On the basis of S4max data retrieved from COSMIC GPS radio occultation measurements, the long-term climatology of the intensity of E s layers is investigated for the period from December 2006 to January 2014. The global maps of E s intensity shows a high spatial resolution geographical distributions and strong seasonal dependence of E s layers.The maximum intensity of E s occurs in the midlatitudes, and its value in summer is 2-3 times larger than that in winter. A relatively strong E s layer is observed at the North and South Poles with a distinct boundary dividing the middle latitudes and high latitudes along 60 • -5 80 • geomagnetic latitude bands. Besides, simulation results shows that the convergence of vertical ion velocity could partially explain the seasonal dependence of E s intensity. Furthermore, some disagreements between the distributions of calculated divergence of vertical ion velocity and observed E s intensity indicate that other processes such as magnetic field effects, meteoric mass influx into the earth's atmosphere and chemical processes of metallic ions should also be considered, which play an important role in the spatial and seasonal variations of E s layers. 10 1 IntroductionThe ionospheric sporadic E (E s ) layers are known as thin-layered structures of intense high electron density at 90-130 km altitudes. Rocket-borne mass spectrometric measurements proved that the E s layer is mostly the ionization of metal atoms such as Fe + , Mg + , and Na + (Kopp, 1997;Grebowsky and Aikin, 2002). The E s layer is mainly at midlatitudes and relatively absent at geomagnetic equator and high latitudes (Whitehead, 1989). It is widely accepted that the mechanism responsible for the E s 15 layer formation at midlatitudes is the windshear theory, in which the zonal and meridional winds provide the vertical windshear convergence nodes. As a result, the long-lived metallic ions are forced to converge towards the wind shear null to form a thin layer of intense metallic ionization (Whitehead, 1961;Macleod, 1966;Whitehead, 1970;Nygren et al., 1984;Whitehead, 1989;Haldoupis, 2012). In the equatoral region, the physical process of E s irregularities is attributed to the gradient-drift instabilities associated with the equatorial electrojet (Tsunoda, 2008). The E s layer generally has a vertical scale of 1 km or less, but its 20 1 Atmos. Chem. Phys. Discuss., https://doi.