Abstract. The methods used to identify coherent ocean eddies are either Eulerian or Lagrangian in nature, and nearly all existing eddy datasets are based on
the Eulerian method. In this study, millions of Lagrangian particles are advected by satellite-derived surface geostrophic velocities over the
period of 1993–2019. Using the method of Lagrangian-averaged vorticity deviation (LAVD), we present a global Lagrangian eddy dataset (GLED
v1.0, Liu and Abernathey, 2022, https://doi.org/10.5281/zenodo.7349753). This open-source dataset contains not only the general features (eddy
center position, equivalent radius, rotation property, etc.) of eddies with lifetimes of 30, 90, and 180 d, but also the trajectories of
particles trapped by coherent eddies over the lifetime. We present the statistical features of Lagrangian eddies and compare them with those of the
most widely used sea surface height (SSH) eddies, focusing on generation sites, size, and propagation speed. A remarkable feature is that Lagrangian
eddies are generally smaller than SSH eddies, with a radius ratio of about 0.5. Also, the validation using Argo floats indicates that coherent
eddies from GLED v1.0 exist in the real ocean and have the ability to transport water parcels. Our eddy dataset provides an additional option for
oceanographers to understand the interaction between coherent eddies and other physical or biochemical processes in the Earth system.