D-region ion–neutral coupled chemistry (Sodankylä Ion Chemistry, SIC) within the Whole Atmosphere Community Climate Model (WACCM 4) – WACCM-SIC and WACCM-rSIC

This study presents a new ion–neutral chemical model coupled into the Whole Atmosphere Community Climate Model (WACCM). The ionospheric D-region (altitudes ∼  50–90 km) chemistry is based on the Sodankylä Ion Chemistry (SIC) model, a one-dimensional model containing 307 ion–neutral and ion recombination, 16 photodissociation and 7 photoionization reactions of neutral species, positive and negative ions, and electrons. The SIC mechanism was reduced using the simulation error minimization connectivity method (SEM-CM) to produce a reaction scheme of 181 ion–molecule reactions of 181 ion–molecule reactions of 27 positive and 18 negative ions. This scheme describes the concentration profiles at altitudes between 20 km and 120 km of a set of major neutral species (HNO3, O3, H2O2, NO, NO2, HO2, OH, N2O5) and ions (O2+, O4+, NO+, NO+(H2O), O2+(H2O), H+(H2O), H+(H2O)2, H+(H2O)3, H+(H2O)4, O3−, NO2−, O−, O2, OH−, O2−(H2O), O2−(H2O)2, O4−, CO3−, CO3−(H2O), CO4−, HCO3−, NO2−, NO3−, NO3−(H2O), NO3−(H2O)2, NO3−(HNO3), NO3−(HNO3)2, Cl−, ClO−), which agree with the full SIC mechanism within a 5 % tolerance. Four 3-D model simulations were then performed, using the impact of the January 2005 solar proton event (SPE) on D-region HOx and NOx chemistry as a test case of four different model versions: the standard WACCM (no negative ions and a very limited set of positive ions); WACCM-SIC (standard WACCM with the full SIC chemistry of positive and negative ions); WACCM-D (standard WACCM with a heuristic reduction of the SIC chemistry, recently used to examine HNO3 formation following an SPE); and WACCM-rSIC (standard WACCM with a reduction of SIC chemistry using the SEM-CM method). The standard WACCM misses the HNO3 enhancement during the SPE, while the full and reduced model versions predict significant NOx, HOx and HNO3 enhancements in the mesosphere during solar proton events. The SEM-CM reduction also identifies the important ion–molecule reactions that affect the partitioning of odd nitrogen (NOx), odd hydrogen (HOx) and O3 in the stratosphere and mesosphere

Observed propagation route of VLF transmitter signals in the magnetosphere

Signals of powerful ground transmitters at various places have been detected by satellites in near‐Earth space. The study on propagation mode, ducted or nonducted, has attracted much attentions for several decades. Based on the statistical results from Van Allen Probes (data from October 2012 to March 2017) and DEMETER satellite (from January 2006 to December 2007), we present the ground transmitter signals distributed clearly in ionosphere and magnetosphere. The observed propagation route in the meridian plane in the magnetosphere for each of various transmitters from the combination of DEMETER and Van Allen Probes data in nighttime is revealed for the first time. We use realistic ray tracing simulation and compare simulation results against Van Allen Probes and DEMETER observation. By comparison we demonstrate that the observed propagation route, with partial deviation from the field lines corresponding to ground stations, provides direct and clear statistical evidence that the nonducted propagation mode plays a main role, although with partial contribution from ducted propagation. The propagation characteristics of VLF transmitter signals in the magnetosphere are critical for quantitatively assessing their contribution to energetic electron loss in radiation belts

Observations of Radiation Belt Losses Due to Cyclotron Wave-Particle Interactions

Electron loss to the atmosphere plays a critical role in driving dynamics of the Earths Van Allen radiation belts and slot region. This is a review of atmospheric loss of radiation belt electrons caused by plasma wave scattering via Doppler-shifted cyclotron resonance. In particular, the focus is on observational signatures of electron loss, which include direct measurements of precipitating electrons, measured properties of waves that drive precipitation, and variations in the trapped population resulting from loss. We discuss wave and precipitation measurements from recent missions, including simultaneous multi-payload observations, which have provided new insight into the dynamic nature of the radiation belts

Examination of Radiation Belt Dynamics during Substorm Clusters: Activity Drivers and Dependencies of Trapped Flux Enhancements

Dynamical variations of radiation belt trapped electron fluxes are examined to better understand the variability of enhancements linked to substorm clusters. Analysis is undertaken using the Substorm Onsets and Phases from Indices of the Electrojet substorm cluster algorithm for event detection. Observations from low earth orbit are complemented by additional measurements from medium earth orbit to allow a major expansion in the energy range considered, from medium energy energetic electrons up to ultra-relativistic electrons. The number of substorms identified inside a cluster does not depend strongly on solar wind drivers or geomagnetic indices either before, during, or after the cluster start time. Clusters of substorms linked to moderate (100 nT < AE ≤ 300 nT) or strong AE (AE ≥ 300 nT) disturbances are associated with radiation belt flux enhancements, including up to ultra-relativistic energies by the strongest substorms (as measured by strong southward Bz and high AE). These clusters reliably occur during times of high speed solar winds streams with associated increased magnetospheric convection. However, substorm clusters associated with quiet AE disturbances (AE ≤ 100 nT) lead to no significant chorus whistler mode intensity enhancements, or increases in energetic, relativistic, or ultra-relativistic electron flux in the outer radiation belts. In these cases the solar wind speed is low, and the geomagnetic Kp index indicates a lack of magnetospheric convection. Our study clearly indicates that clusters of substorms occurring outside of high speed wind streams are not by themselves sufficient to drive acceleration, which may be due to the lack of pre-cluster convection

Solar Forcing for CMIP6 (v3.1)

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