Abstract. This study investigated the statistics of eddy splitting and merging in the global oceans based on 23 years’ altimetry data. Multicore structures were identified using an improved threshold-free closed-contour algorithm of sea surface height. Splitting and merging events were discerned from continuous time series maps of sea level anomalies. Multicore structures represent an intermediate stage in the process of eddy evolution, similar to the generation of multiple nuclei in a cell as a preparatory phase for cell division. Generally, splitting or merging events can change substantially (by a factor of two or more) the eddy scale, amplitude, and eddy kinetic energy. Specifically, merging (splitting) generally causes an increase (decrease) of eddy properties. Multicore eddies were found to tend to split into two eddies with different intensities. Similarly, eddy merging is not an interaction of two equal-intensity eddies, and that it tends to manifest as a strong eddy merging with a weaker one. A hybrid tracking strategy based on the eddy overlap ratio, considering both multicore and single-core eddies, was used to confirm splitting and merging events globally. The census revealed that eddy splitting and merging do not always occur most frequently in eddy-rich regions, e.g., their frequencies of occurrence in the Antarctic Circumpolar Current and western boundary currents were found obviously higher than mid-latitude regions (20°–35°) north and south. Eddy splitting and merging are caused primarily by an unstable configuration of multicore structures due to obvious current– or eddy–topography interaction, strong current variation, and eddy–mean flow interaction.