ALBERT

Description: While observed mesospheric polar nitric acid enhancements have been attributed to energetic particle precipitation through ion cluster chemistry in the past, this phenomenon is not reproduced in current whole-atmosphere chemistry-climate models. We investigate such nitric acid enhancements resulting from energetic electron precipitation events using a recently developed variant of the Whole Atmosphere Community Climate Model (WACCM) that includes a sophisticated ion chemistry tailored for the D-layer of the ionosphere (50–90 km), namely, WACCM-D. Using the specified dynamics mode, that is, nudging dynamics in the troposphere and stratosphere to meteorological reanalyses, we perform a 1-year-long simulation (July 2009–June 2010) and contrast WACCM-D with the standard WACCM. Both WACCM and WACCM-D simulations are performed with and without forcing from medium-to-high energy electron precipitation, allowing a better representation of the energetic electrons penetrating into the mesosphere. We demonstrate the effects of the strong particle precipitation events which occurred during April and May 2010 on nitric acid and on key ion cluster species, as well as other relevant species of the nitrogen family. The 1-year-long simulation allows the event-related changes in neutral and ionic species to be placed in the context of their annual cycle. We especially highlight the role played by medium-to-high energy electrons in triggering ion cluster chemistry and ion-ion recombinations in the mesosphere and lower thermosphere during the precipitation event, leading to enhanced production of nitric acid and raising its abundance by 2 orders of magnitude from 10−4 to a few 10−2 ppb. ©2018. American Geophysical Union. All Rights Reserved.

Description: The radiation belts are regions in the near-Earth space where solar wind electrons are captured by the Earth's magnetic field. A portion of these electrons is continuously lost into the atmosphere where they cause ionization and chemical changes. Driven by the solar activity, the electron forcing leads to ozone variability in the polar stratosphere and mesosphere. Understanding the possible dynamical connections to regional climate is an ongoing research activity which supports the assessment of greenhouse-gas-driven climate change by a better definition of the solar-driven variability. In the context of the Coupled Model Intercomparison Project Phase 6 (CMIP6), energetic electron and proton precipitation is included in the solar-forcing recommendation for the first time. For the radiation belt electrons, the CMIP6 forcing is from a daily zonal-mean proxy model. This zonal-mean model ignores the well-known dependency of precipitation on magnetic local time (MLT), i.e. its diurnal variability. Here we use the Whole Atmosphere Community Climate Model with its lower-ionospheric-chemistry extension (WACCM-D) to study effects of the MLT dependency of electron forcing on the polar-ozone response. We analyse simulations applying MLT-dependent and MLT-independent forcings and contrast the resulting ozone responses in monthly-mean data as well as in monthly means at individual local times. We consider two cases: (1) the year 2003 and (2) an extreme, continuous forcing. Our results indicate that the ozone responses to the MLT-dependent and the MLT-independent forcings are very similar, and the differences found are small compared to those caused by the overall uncertainties related to the representation of electron forcing in climate simulations. We conclude that the use of daily zonal-mean electron forcing will provide an accurate ozone response in long-term climate simulations.

Description: Atmospheric effects of solar proton events (SPEs) have been studied for decades, because their drastic impact can be used to test our understanding of upper stratospheric and mesospheric chemistry in the polar cap regions. For example, odd hydrogen and odd nitrogen are produced during SPEs, which leads to depletion of ozone in catalytic reactions, such that the effects are easily observed from satellites during the strongest events. Until recently, the complexity of the ion chemistry in the lower ionosphere (i.e., in the D region) has restricted global models to simplified parameterizations of chemical impacts induced by energetic particle precipitation (EPP). Because of this restriction, global models have been unable to correctly reproduce some important effects, such as the increase in mesospheric HNO3 or the changes in chlorine species. Here we use simulations from the WACCM-D model, a variant of the Whole Atmosphere Community Climate Model, to study the statistical response of the atmosphere to the 66 strongest SPEs which occurred in the years 1989–2012. Our model includes a set of D-region ion chemistry, designed for a detailed representation of the atmospheric effects of SPEs and EPP in general. We use superposed epoch analysis to study changes in O3, HOx (OH + HO2), Clx (Cl + ClO), HNO3, NOx (NO + NO2) and H2O. Compared to the standard WACCM which uses an ion chemistry parameterization, WACCM-D produces a larger response in O3 and NOx and a weaker response in HOx and introduces changes in HNO3 and Clx. These differences between WACCM and WACCM-D highlight the importance of including ion chemistry reactions in models used to study EPP.

Description: Most of our understanding of the atmosphere is based on observations and their comparison with model simulations. In middle atmosphere studies it is common practice to use an approach, where the model dynamics are at least partly based on temperature and wind fields from an external meteorological model. In this work we test how closely satellite measurements of a few central trace gases agree with this kind of model simulation. We use collocated vertical profiles where each satellite measurement is compared to the closest model data. We compare profiles and distributions of O3, NO2 and NO3 from the Global Ozone Monitoring by Occultation of Stars instrument (GOMOS) on the Envisat satellite with simulations by the Whole Atmosphere Community Climate Model (WACCM). GOMOS measurements are from nighttime. Our comparisons show that in the stratosphere outside the polar regions differences in ozone between WACCM and GOMOS are small, between 0 and 6%. The correlation of 5-day time series show a very high 0.9–0.95. In the tropical region 10° S–10° N below 10 hPa WACCM values are up to 20 % larger than GOMOS. In the Arctic below 6 hPa WACCM ozone values are up to 20 % larger than GOMOS. In the mesosphere between 0.04 and 1 hPa the WACCM is at most 20 % smaller than GOMOS. Above the ozone minimum at 0.01 hPa (or 80 km) large differences are found between WACCM and GOMOS. The correlation can still be high, but at the second ozone peak the correlation falls strongly and the ozone abundance from WACCM is about 60 % smaller than that from GOMOS. The total ozone columns (above 50 hPa) of GOMOS and WACCM agree within ±2 % except in the Arctic where WACCM is 10 % larger than GOMOS. Outside the polar areas and in the validity region of GOMOS NO2 measurements (0.3–37 hPa) WACCM and GOMOS NO2 agree within −5 to +25 % and the correlation is high (0.7–0.95) except in the upper stratosphere at the southern latitudes. In the polar areas, where solar particle precipitation and downward transport from the thermosphere enhance NO2 abundance, large differences up to −90 % are found between WACCM and GOMOS NO2 and the correlation varies between 0.3 and 0.9. For NO3, we find that the WACCM and GOMOS difference is between −20 and 5 % with a very high correlation of 0.7–0.95. We show that NO3 values strongly depend on temperature and the dependency can be fitted by the exponential function of temperature. The ratio of NO3 to O3 from WACCM and GOMOS closely follow the prediction from the equilibrium chemical theory. Abrupt temperature increases from sudden stratospheric warmings (SSWs) are reflected as sudden enhancements of WACCM and GOMOS NO3 values.

Description: The Ku-band microwave frequencies (10.70–14.25 GHz) overlap emissions from ozone (O3) at 11.072 GHz and hydroxyl radical (OH) at 13.441 GHz. These important chemical species in the polar middle atmosphere respond strongly to high-latitude geomagnetic activity associated with space weather. Atmospheric model calculations predict that energetic electron precipitation (EEP) driven by magnetospheric substorms produces large changes in polar mesospheric O3 and OH. The EEP typically peaks at geomagnetic latitudes of ∼65∘ and evolves rapidly with time longitudinally and over the geomagnetic latitude range 60–80∘. Previous atmospheric modelling studies have shown that during substorms OH abundance can increase by more than an order of magnitude at 64–84 km and mesospheric O3 losses can exceed 50 %. In this work, an atmospheric simulation and retrieval study has been performed to determine the requirements for passive microwave radiometers capable of measuring diurnal variations in O3 and OH profiles from high-latitude Northern Hemisphere and Antarctic locations to verify model predictions. We show that, for a 11.072 GHz radiometer making 6 h spectral measurements with 10 kHz frequency resolution and root-mean-square baseline noise of 1 mK, O3 could be profiled over 8×10-4–0.22 hPa (∼98–58 km) with 10–17 km height resolution and ∼1 ppmv uncertainty. For the equivalent 13.441 GHz measurements with vertical sensor polarisation, OH could be profiled over 3×10-3–0.29 hPa (∼90–56 km) with 10–17 km height resolution and ∼3 ppbv uncertainty. The proposed observations would be highly applicable to studies of EEP, atmospheric dynamics, planetary-scale circulation, chemical transport, and the representation of these processes in polar and global climate models. Such observations would provide a relatively low-cost alternative to increasingly sparse satellite measurements of the polar middle atmosphere, extending long-term data records and also providing “ground truth” calibration data.

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: This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models.For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %).We compare the climatic effects of the CMIP6 solar forcing dataset to its CMIP5 predecessor by using time-slice experiments of two chemistry–climate models and a reference radiative transfer model. The differences in the long-term mean SSI in the CMIP6 dataset, compared to CMIP5, impact on climatological stratospheric conditions (lower shortwave heating rates of −0.35 K day−1 at the stratopause), cooler stratospheric temperatures (−1.5 K in the upper stratosphere), lower ozone abundances in the lower stratosphere (−3 %), and higher ozone abundances (+1.5 % in the upper stratosphere and lower mesosphere). Between the maximum and minimum phases of the 11-year solar cycle, there is an increase in shortwave heating rates (+0.2 K day−1 at the stratopause), temperatures ( ∼  1 K at the stratopause), and ozone (+2.5 % in the upper stratosphere) in the tropical upper stratosphere using the CMIP6 forcing dataset. This solar-cycle response is slightly larger, but not statistically significantly different from that for the CMIP5 forcing dataset.CMIP6 models with a well-resolved shortwave radiation scheme are encouraged to prescribe SSI changes and include solar-induced stratospheric ozone variations, in order to better represent solar climate variability compared to models that only prescribe TSI and/or exclude the solar-ozone response. We show that monthly-mean solar-induced ozone variations are implicitly included in the SPARC/CCMI CMIP6 Ozone Database for historical simulations, which is derived from transient chemistry–climate model simulations and has been developed for climate models that do not calculate ozone interactively. CMIP6 models without chemistry that perform a preindustrial control simulation with time-varying solar forcing will need to use a modified version of the SPARC/CCMI Ozone Database that includes solar variability. CMIP6 models with interactive chemistry are also encouraged to use the particle forcing datasets, which will allow the potential long-term effects of particles to be addressed for the first time. The consideration of particle forcing has been shown to significantly improve the representation of reactive nitrogen and ozone variability in the polar middle atmosphere, eventually resulting in further improvements in the representation of solar climate variability in global models.