Modeling of middle atmosphere dynamics with LIMA

• Published in: Journal of atmospheric and solar-terrestrial physics 2008-06, Vol.70 (8), p.1170-1200
• Main Author: Berger, U.
• Format: Article
• Language: English
• Subjects: Earth, ocean, space; Exact sciences and technology; External geophysics; Middle atmosphere dynamics; Numerical modeling; Physics of the ionosphere; Physics of the magnetosphere; Summer mesopause temperatures
• ISSN: 1364-6826; 1879-1824
• DOI: 10.1016/j.jastp.2008.02.004

This paper describes a new circulation model of the middle atmosphere called LIMA (Leibniz-Institute middle atmosphere model) which especially aims to model the thermal structure around mesopause altitudes. LIMA is a fully non-linear, global, and 3-d Eulerian grid-point model which extends from the ground to the lower thermosphere (0–150 km) taking into account major processes of radiation, chemistry, and transport. The major improvements of LIMA compared to its predecessor COMMA/IAP are the implementation of a triangular horizontal grid structure with 41 804 grid points in every horizontal layer ( Δ x ∼ Δ y ∼ 110 km ), and the assimilation of tropospheric and lower stratospheric data from ECMWF/ERA-40. Applying LIMA we are able to simulate in situ wind reversals of the mesospheric jets and cold summer mesopause states without applying any gravity wave parametrization. In the summer season the mesospheric wind reversal of the westward zonal wind is in the order of + 5 –10 m/s which implies that the resolved eddy fields and hence EP-flux divergences, determining the momentum budget of the mesopause region, are underestimated by approximately a factor of two. On the other hand the thermal structure of the summer mesopause, and hence at least the balance of the energy budget in the summer mesopause region, is reproduced. E.g. during early July the observed mesopause temperature from rocket (falling spheres) experiments over ALOMAR ( 69 ∘ N , 16 ∘ E ) is 131 K at 88 km compared to corresponding model temperature of 133 K at 88 km. The results of our numerical simulations suggest (1) that the intrinsic variability of the mesosphere is indeed determined by the lower atmosphere variability, and (2) that the simulated eddy fields may be interpreted as some fraction of inertia gravity waves resolved by the model ( λ xy ⩾ 500 km , 2 h ⩽ T ⩽ 12 h ).