ALBERT

Description: Variations in the solar spectral irradiance (SSI) with the 11-year sunspot cycle have been shown to have a significant impact on temperatures and the mixing ratios of atmospheric constituents in the stratosphere and mesosphere. Uncertainties in modelling the effects of SSI variations arise from uncertainties in the empirical models reconstructing the prescribed SSI data set as well as from uncertainties in the chemistry–climate model (CCM) formulation. In this study CCM simulations with the ECHAM/MESSy Atmospheric Chemistry (EMAC) model and the Community Earth System Model 1 (CESM1)–Whole Atmosphere Chemistry Climate Model (WACCM) have been performed to quantify the uncertainties of the solar responses in chemistry and dynamics that are due to the usage of five different SSI data sets or the two CCMs. We apply a two-way analysis of variance (ANOVA) to separate the influence of the SSI data sets and the CCMs on the variability of the solar response in shortwave heating rates, temperature, and ozone. The solar response is derived from climatological differences of time slice simulations prescribing SSI for the solar maximum in 1989 and near the solar minimum in 1994. The SSI values for the solar maximum of each SSI data set are created by adding the SSI differences between November 1994 and November 1989 to a common SSI reference spectrum for near-solar-minimum conditions based on ATLAS-3 (Atmospheric Laboratory of Applications and Science-3). The ANOVA identifies the SSI data set with the strongest influence on the variability of the solar response in shortwave heating rates in the upper mesosphere and in the upper stratosphere–lower mesosphere. The strongest influence on the variability of the solar response in ozone and temperature is identified in the upper stratosphere–lower mesosphere. However, in the region of the largest ozone mixing ratio, in the stratosphere from 50 to 10 hPa, the SSI data sets do not contribute much to the variability of the solar response when the Spectral And Total Irradiance REconstructions-T (SATIRE-T) SSI data set is omitted. The largest influence of the CCMs on variability of the solar responses can be identified in the upper mesosphere. The solar response in the lower stratosphere also depends on the CCM used, especially in the tropics and northern hemispheric subtropics and mid-latitudes, where the model dynamics modulate the solar responses. Apart from the upper mesosphere, there are also regions where the largest fraction of the variability of the solar response is explained by randomness, especially for the solar response in temperature.

Description: In response to global warming, the Brewer–Dobson circulation in the stratosphere is expected to accelerate and the mean transport time of air along this circulation to decrease. This would imply a negative stratospheric age of air trend, i.e. an air parcel would need less time to travel from the tropopause to any point in the stratosphere. Age of air as inferred from tracer observations, however, shows zero to positive trends in the northern mid-latitude stratosphere and zonally asymmetric patterns. Using satellite observations and model calculations we show that the observed latitudinal and vertical patterns of the decadal changes of age of air in the lower to middle stratosphere during the period 2002–2012 are predominantly caused by a southward shift of the circulation pattern by about 5°. After correction for this shift, the observations reveal a hemispherically almost symmetric decrease of age of air in the lower to middle stratosphere up to 800 K of up to −0.25 years over the 2002–2012 period with strongest decrease in the northern tropics. This net change is consistent with long-term trends from model predictions.

Description: We compare simulations from three high-top (with upper lid above 120 km) and five medium-top (with upper lid around 80 km) atmospheric models with observations of odd nitrogen (NOx  =  NO + NO2), temperature, and carbon monoxide from seven satellite instruments (ACE-FTS on SciSat, GOMOS, MIPAS, and SCIAMACHY on Envisat, MLS on Aura, SABER on TIMED, and SMR on Odin) during the Northern Hemisphere (NH) polar winter 2008/2009. The models included in the comparison are the 3-D chemistry transport model 3dCTM, the ECHAM5/MESSy Atmospheric Chemistry (EMAC) model, FinROSE, the Hamburg Model of the Neutral and Ionized Atmosphere (HAMMONIA), the Karlsruhe Simulation Model of the Middle Atmosphere (KASIMA), the modelling tools for SOlar Climate Ozone Links studies (SOCOL and CAO-SOCOL), and the Whole Atmosphere Community Climate Model (WACCM4). The comparison focuses on the energetic particle precipitation (EPP) indirect effect, that is, the polar winter descent of NOx largely produced by EPP in the mesosphere and lower thermosphere. A particular emphasis is given to the impact of the sudden stratospheric warming (SSW) in January 2009 and the subsequent elevated stratopause (ES) event associated with enhanced descent of mesospheric air. The chemistry climate model simulations have been nudged toward reanalysis data in the troposphere and stratosphere while being unconstrained above. An odd nitrogen upper boundary condition obtained from MIPAS observations has further been applied to medium-top models. Most models provide a good representation of the mesospheric tracer descent in general, and the EPP indirect effect in particular, during the unperturbed (pre-SSW) period of the NH winter 2008/2009. The observed NOx descent into the lower mesosphere and stratosphere is generally reproduced within 20 %. Larger discrepancies of a few model simulations could be traced back either to the impact of the models' gravity wave drag scheme on the polar wintertime meridional circulation or to a combination of prescribed NOx mixing ratio at the uppermost model layer and low vertical resolution. In March–April, after the ES event, however, modelled mesospheric and stratospheric NOx distributions deviate significantly from the observations. The too-fast and early downward propagation of the NOx tongue, encountered in most simulations, coincides with a temperature high bias in the lower mesosphere (0.2–0.05 hPa), likely caused by an overestimation of descent velocities. In contrast, upper-mesospheric temperatures (at 0.05–0.001 hPa) are generally underestimated by the high-top models after the onset of the ES event, being indicative for too-slow descent and hence too-low NOx fluxes. As a consequence, the magnitude of the simulated NOx tongue is generally underestimated by these models. Descending NOx amounts simulated with medium-top models are on average closer to the observations but show a large spread of up to several hundred percent. This is primarily attributed to the different vertical model domains in which the NOx upper boundary condition is applied. In general, the intercomparison demonstrates the ability of state-of-the-art atmospheric models to reproduce the EPP indirect effect in dynamically and geomagnetically quiescent NH winter conditions. The encountered differences between observed and simulated NOx, CO, and temperature distributions during the perturbed phase of the 2009 NH winter, however, emphasize the need for model improvements in the dynamical representation of elevated stratopause events in order to allow for a better description of the EPP indirect effect under these particular conditions.

Description: We analyze the impact of energetic particle precipitation on the stratospheric nitrogen budget, ozone abundances and net radiative heating using results from three global chemistry-climate models considering solar protons and geomagnetic forcing due to auroral or radiation belt electrons. Two of the models cover the atmosphere up to the lower thermosphere, the source region of auroral NO production. Geomagnetic forcing in these models is included by prescribed ionization rates. One model reaches up to about 80 km, and geomagnetic forcing is included by applying an upper boundary condition of auroral NO mixing ratios parameterized as a function of geomagnetic activity. Despite the differences in the implementation of the particle effect, the resulting modeled NOy in the upper mesosphere agrees well between all three models, demonstrating that geomagnetic forcing is represented in a consistent way either by prescribing ionization rates or by prescribing NOy at the model top.Compared with observations of stratospheric and mesospheric NOy from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) instrument for the years 2002–2010, the model simulations reproduce the spatial pattern and temporal evolution well. However, after strong sudden stratospheric warmings, particle-induced NOy is underestimated by both high-top models, and after the solar proton event in October 2003, NOy is overestimated by all three models. Model results indicate that the large solar proton event in October 2003 contributed about 1–2 Gmol (109 mol) NOy per hemisphere to the stratospheric NOy budget, while downwelling of auroral NOx from the upper mesosphere and lower thermosphere contributes up to 4 Gmol NOy. Accumulation over time leads to a constant particle-induced background of about 0.5–1 Gmol per hemisphere during solar minimum, and up to 2 Gmol per hemisphere during solar maximum. Related negative anomalies of ozone are predicted by the models in nearly every polar winter, ranging from 10–50 % during solar maximum to 2–10 % during solar minimum. Ozone loss continues throughout polar summer after strong solar proton events in the Southern Hemisphere and after large sudden stratospheric warmings in the Northern Hemisphere. During mid-winter, the ozone loss causes a reduction of the infrared radiative cooling, i.e., a positive change of the net radiative heating (effective warming), in agreement with analyses of geomagnetic forcing in stratospheric temperatures which show a warming in the late winter upper stratosphere. In late winter and spring, the sign of the net radiative heating change turns to negative (effective cooling). This spring-time cooling lasts well into summer and continues until the following autumn after large solar proton events in the Southern Hemisphere, and after sudden stratospheric warmings in the Northern Hemisphere.

Description: We present a recently developed emission module for the ICON (ICOsahedral Non-hydrostatic)-ART (Aerosols and Reactive Trace gases) modelling framework. The emission module processes external flux data sets and increments the tracer volume mixing ratios in the boundary layer accordingly. The performance of the emission module is illustrated with simulations of acetone, using a simplified chemical depletion mechanism based on a reaction with OH and photolysis only. In our model setup, we calculate a tropospheric acetone lifetime of 33 days, which is in good agreement with the literature. We compare our results with ground-based as well as with airborne IAGOS-CARIBIC measurements in the upper troposphere and lowermost stratosphere (UTLS) in terms of phase and amplitude of the annual cycle. In all our ICON-ART simulations the general seasonal variability is well represented but uncertainties remain concerning the magnitude of the acetone mixing ratio in the UTLS region. In addition, the module for online calculations of biogenic emissions (MEGAN2.1) is implemented in ICON-ART and can replace the offline biogenic emission data sets. In a sensitivity study we show how different parametrisations of the leaf area index (LAI) change the emission fluxes calculated by MEGAN2.1 and demonstrate the importance of an adequate treatment of the LAI within MEGAN2.1. We conclude that the emission module performs well with offline and online emission fluxes and allows the simulation of the annual cycles of emissions-dominated substances.

Description: We compare simulations from three high-top (with upper lid above 120 km) and five medium-top (with upper lid around 80km) atmospheric models with observations of odd nitrogen (NOx = NO + NO2), temperature, and carbon monoxide from seven satellite instruments (ACE-FTS on SciSat, GOMOS, MIPAS, and SCIAMACHY on Envisat, MLS on Aura, SABER on TIMED, and SMR on Odin) during the Northern Hemisphere (NH) polar winter 2008/2009. The models included in the comparison are the 3d Chemistry Transport model (3dCTM), the ECHAM5/MESSy Atmospheric Chemistry (EMAC) model, FinROSE, the Hamburg Model of the Neutral and Ionized Atmosphere (HAMMONIA), the Karlsruhe Simulation Model of the Middle Atmosphere (KASIMA), the modeling tools for SOlar Climate Ozone Links studies (SOCOL and CAO-SOCOL), and the Whole Atmosphere Community Climate Model (WACCM4). The comparison focuses on the energetic particle precipitation (EPP) indirect effect, that is, the polar winter descent of NOx largely produced by EPP in the mesosphere and lower thermosphere. A particular emphasis is given to the impact of the sudden stratospheric warming (SSW) in January 2009 and the subsequent elevated stratopause (ES) event associated with enhanced descent of mesospheric air. The chemistry climate model simulations have been nudged toward reanalysis data in the troposphere and stratosphere while being unconstrained above. An odd nitrogen upper boundary condition obtained from MIPAS observations has further been applied to medium-top models. Most models provide a good representation of the mesospheric tracer descent in general, and the EPP indirect effect in particular, during the unperturbed (pre-SSW) period of the NH winter 2008/2009. The observed NOx descent into the lower mesosphere and stratosphere is generally reproduced within 20 %. Larger discrepancies of a few model simulations could be traced back either to the impact of the models' gravity wave drag scheme on the polar wintertime meridional circulation or to a combination of prescribed NOx mixing ratio at the uppermost model layer and low vertical resolution. In March–April, after the ES event, however, modelled mesospheric and stratospheric NOx distributions deviate significantly from the observations. The too fast and early downward propagation of the NOx tongue, encountered in most simulations, coincides with a temperature high bias in the lower mesosphere (0.2–0.05 hPa) being likely caused by an overestimation of descent velocities. On the other hand, upper mesospheric temperatures (at 0.05–0.001 hPa) are generally underestimated by the high-top models after the onset of the ES event, being indicative for too slow descent and hence too low NOx fluxes. As a consequence, the magnitude of the simulated NOx tongue is generally underestimated by these models. Descending NOx amounts simulated with medium-top models are on average closer to the observations but show a large spread of up to several hundred percent. This is primarily attributed to the different vertical model domains in which the NOx upper boundary condition is applied. In general, the intercomparison demonstrates the ability of state-of-the-art atmospheric models to reproduce the EPP indirect effect in dynamically and geomagnetically quiescent NH winter conditions. The encountered differences between observed and simulated NOx, CO, and temperature distributions during the perturbed phase of the 2009 NH winter, however, emphasize the need for model improvements in the dynamical representation of elevated stratopause events in order to allow for a better description of the EPP indirect effect under these particular conditions.

Description: We analyze the impact of energetic particle precipitation on the stratospheric nitrogen budget, ozone abundances and net radiative heating using results from three global chemistry-climate models considering solar protons and geomagnetic forcing due to auroral or radiation belt electrons. Two of the models cover the atmosphere up to the lower thermosphere, the source region of auroral NO production. Geomagnetic forcing in these models is included by prescribed ionization rates. One model reaches up to about 80 km, and geomagnetic forcing is included by applying an upper boundary condition of auroral NO mixing ratios parameterized as a function of geomagnetic activity. Despite the differences in the implementation of the particle effect, the resulting modeled NOy in the upper mesosphere agrees well between all three models, demonstrating that geomagnetic forcing is represented in a consistent way either by prescribing ionization rates or by prescribing NOy at the model top. Compared with observations of stratospheric and mesospheric NOy from the MIPAS instrument for the years 2002–2010, the model simulations reproduce the spatial pattern and temporal evolution well. However, after strong sudden stratospheric warmings, particle induced NOy is underestimated by both high-top models, and after the solar proton event in October 2003, NOy is overestimated by all three models. Model results indicate that the large solar proton event in October 2003 contributed about 1–2 Gmol (109 mol) NOy per hemisphere to the stratospheric NOy budget, while downwelling of auroral NOx from the upper mesosphere and lower thermosphere contributes up to 4 Gmol NOy. Accumulation over time leads to a constant particle-induced background of about 0.5–1 Gmol per hemisphere during solar minimum, and up to 2 Gmol per hemisphere during solar maximum. Related negative anomalies of ozone are predicted by the models nearly in every polar winter, ranging from 10–50 % during solar maximum to 2–10 % during solar minimum. Ozone loss continues throughout polar summer after strong solar proton events in the Southern hemisphere and after large sudden stratospheric warmings in the Northern hemisphere. During mid-winter, the ozone loss causes a reduction of the infrared radiative cooling, i.e., a positive change of the net radiative heating (effective warming), in agreement with analyses of geomagnetic forcing in stratospheric temperatures which show a warming in the late winter upper stratosphere. In late winter and spring, the sign of the net radiative heating change turns to negative (effective cooling). This spring-time cooling lasts well into summer and continues until the following autumn after large solar proton events in the Southern hemisphere, after sudden stratospheric warmings in the Northern hemisphere.