Revisiting the mystery of recent stratospheric temperature trends.

Simulated stratospheric temperatures over the period 1979-2016 in models from the Chemistry-Climate Model Initiative (CCMI) are compared with recently updated and extended satellite observations. The multi-model mean global temperature trends over 1979- 2005 are -0.88 ± 0.23, -0.70 ± 0.16, and -0.50 ± 0.12 K decade-1 for the Stratospheric Sounding Unit (SSU) channels 3 (~40-50 km), 2 (~35-45 km), and 1 (~25-35 km), respectively. These are within the uncertainty bounds of the observed temperature trends from two reprocessed satellite datasets. In the lower stratosphere, the multi-model mean trend in global temperature for the Microwave Sounding Unit channel 4 (~13-22 km) is -0.25 ± 0.12 K decade-1 over 1979-2005, consistent with estimates from three versions of this satellite record. The simulated stratospheric temperature trends in CCMI models over 1979-2005 agree with the previous generation of chemistry-climate models. The models and an extended satellite dataset of SSU with the Advanced Microwave Sounding Unit-A show weaker global stratospheric cooling over 1998-2016 compared to the period of intensive ozone depletion (1979-1997). This is due to the reduction in ozone-induced cooling from the slow-down of ozone trends and the onset of ozone recovery since the late 1990s. In summary, the results show much better consistency between simulated and satellite observed stratospheric temperature trends than was reported by Thompson et al. (2012) for the previous versions of the SSU record and chemistry-climate models. The improved agreement mainly comes from updates to the satellite records; the range of simulated trends is comparable to the previous generation of models.

Stratospheric water vapor affecting atmospheric circulation.

Water vapor plays an important role in many aspects of the climate system, by affecting radiation, cloud formation, atmospheric chemistry and dynamics. Even the low stratospheric water vapor content provides an important climate feedback, but current climate models show a substantial moist bias in the lowermost stratosphere. Here we report crucial sensitivity of the atmospheric circulation in the stratosphere and troposphere to the abundance of water vapor in the lowermost stratosphere. We show from a mechanistic climate model experiment and inter-model variability that lowermost stratospheric water vapor decreases local temperatures, and thereby causes an upward and poleward shift of subtropical jets, a strengthening of the stratospheric circulation, a poleward shift of the tropospheric eddy-driven jet and regional climate impacts. The mechanistic model experiment in combination with atmospheric observations further shows that the prevailing moist bias in current models is likely caused by the transport scheme, and can be alleviated by employing a less diffusive Lagrangian scheme. The related effects on atmospheric circulation are of similar magnitude as climate change effects. Hence, lowermost stratospheric water vapor exerts a first order effect on atmospheric circulation and improving its representation in models offers promising prospects for future research.

Response to Comment on "A global environmental crisis 42,000 years ago".

Our study on the exact timing and the potential climatic, environmental, and evolutionary consequences of the Laschamps Geomagnetic Excursion has generated the hypothesis that geomagnetism represents an unrecognized driver in environmental and evolutionary change. It is important for this hypothesis to be tested with new data, and encouragingly, none of the studies presented by Picin et al. undermine our model.

Response to Comment on "A global environmental crisis 42,000 years ago".

Our paper about the impacts of the Laschamps Geomagnetic Excursion 42,000 years ago has provoked considerable scientific and public interest, particularly in the so-called Adams Event associated with the initial transition of the magnetic poles. Although we welcome the opportunity to discuss our new ideas, Hawks' assertions of misrepresentation are especially disappointing given his limited examination of the material.

A global environmental crisis 42,000 years ago.

Geological archives record multiple reversals of Earth's magnetic poles, but the global impacts of these events, if any, remain unclear. Uncertain radiocarbon calibration has limited investigation of the potential effects of the last major magnetic inversion, known as the Laschamps Excursion [41 to 42 thousand years ago (ka)]. We use ancient New Zealand kauri trees (Agathis australis) to develop a detailed record of atmospheric radiocarbon levels across the Laschamps Excursion. We precisely characterize the geomagnetic reversal and perform global chemistry-climate modeling and detailed radiocarbon dating of paleoenvironmental records to investigate impacts. We find that geomagnetic field minima ~42 ka, in combination with Grand Solar Minima, caused substantial changes in atmospheric ozone concentration and circulation, driving synchronous global climate shifts that caused major environmental changes, extinction events, and transformations in the archaeological record.

The representation of ionospheric potential in the global chemistry-climate model SOCOL.

In this paper, we present the first results of the ionospheric potential (IP) calculations with the chemistry-climate model (CCM) SOCOL (Solar Climate Ozone Links). For the study, we exploit a parameterization of the difference in electric potential between Earth's surface and lower boundary of the ionosphere as a function of thunderstorm and electrified cloud properties. The model shows a good enough agreement with the IP obtained by balloon soundings. The simulated UT variation of IP exhibits a maximum at 20 Universal time (UT) and minimum at about 2 UT which agree with the UT cycle of the lightning activity. The obtained results allow understanding of IP variability pattern at diurnal, seasonal and annual timescales. We also compare our results with the IP simulated with the climate model INMCM4 using similar IP parameterization. The comparison shows a good agreement of UT cycles especially before 12 UT. Simulated IP annual cycle reaches its maximum in late spring in both models. However, the comparison also reveals some differences in amplitudes of IP variability on different time scales. The large deviations occur after 12 UT for all seasons except summer where the maximum of both results happens before 12 UT. The UT cycle of IP simulated with CCM SOCOL is in a better agreement with observations after 12 UT in terms of phase with similar timing of maximum values. The calculation of IP using climate models can help to fill up the gaps when the observed IP is not available. The interactive calculation of IP is also a step forward in coupling atmospheric and ionospheric processes.

Spectra of high energy electron precipitation and atmospheric ionization rates retrieval from balloon measurements.

The bremsstrahlung from high and relativistic energy electron precipitation (HEEP) measured with balloon based instruments provides information on energy spectra and fluence of the precipitating energetic electrons allowing calculations of the atmospheric ionization. HEEP from the outer radiation belt at the subauroral region causes an increase in the ionization rates down to about 20 km altitudes. We study the variability in the ionization rate using the balloon observations of secondary bremsstrahlung initiated by HEEP. For the first time the changes of atmospheric ionization rates on an hourly and minute time scale at different altitudes was retrieved from balloon observations. These new highlights are important for atmospheric electricity that is sensitive to the local condition in the atmosphere including the local ionization rate.

The influence of mixing on stratospheric age of air changes in the 21st century.

Climate models consistently predict an acceleration of the Brewer-Dobson circulation (BDC) due to climate change in the 21st century. However, the strength of this acceleration varies considerably among individual models, which constitutes a notable source of uncertainty for future climate projections. To shed more light upon the magnitude of this uncertainty and on its causes, we analyze the stratospheric mean age of air (AoA) of 10 climate projection simulations from the Chemistry Climate Model Initiative phase 1 (CCMI-I), covering the period between 1960 and 2100. In agreement with previous multi-model studies, we find a large model spread in the magnitude of the AoA trend over the simulation period. Differences between future and past AoA are found to be predominantly due to differences in mixing (reduced aging by mixing and recirculation) rather than differences in residual mean transport. We furthermore analyze the mixing efficiency, a measure of the relative strength of mixing for given residual mean transport, which was previously hypothesized to be a model constant. Here, the mixing efficiency is found to vary not only across models, but also over time in all models. Changes in mixing efficiency are shown to be closely related to changes in AoA and quantified to roughly contribute 10% to the long-term AoA decrease over the 21st century. Additionally, mixing efficiency variations are shown to considerably enhance model spread in AoA changes. To understand these mixing efficiency variations, we also present a consistent dynamical framework based on diffusive closure, which highlights the role of basic state potential vorticity gradients in controlling mixing efficiency and therefore aging by mixing.

Ultraviolet radiation modelling using output from the Chemistry Climate Model Initiative.

We have derived values of the Ultraviolet Index (UVI) at solar noon using the Tropospheric Ultraviolet Model (TUV) driven by ozone, temperature and aerosol fields from climate simulations of the first phase of the Chemistry-Climate Model Initiative (CCMI-1). Since clouds remain one of the largest uncertainties in climate projections, we simulated only the clear-sky UVI. We compared the modelled UVI climatologies against present-day climatological values of UVI derived from both satellite data (the OMI-Aura OMUVBd product) and ground-based measurements (from the NDACC network). Depending on the region, relative differences between the UVI obtained from CCMI/TUV calculations and the ground-based measurements ranged between -5.9% and 10.6%. We then calculated the UVI evolution throughout the 21st century for the four Representative Concentration Pathways (RCPs 2.6, 4.5, 6.0 and 8.5). Compared to 1960s values, we found an average increase in the UVI in 2100 (of 2-4%) in the tropical belt (30°N-30°S). For the mid-latitudes, we observed a 1.8 to 3.4 % increase in the Southern Hemisphere for RCP 2.6, 4.5 and 6.0, and found a 2.3% decrease in RCP 8.5. Higher increases in UVI are projected in the Northern Hemisphere except for RCP 8.5. At high latitudes, ozone recovery is well identified and induces a complete return of mean UVI levels to 1960 values for RCP 8.5 in the Southern Hemisphere. In the Northern Hemisphere, UVI levels in 2100 are higher by 0.5 to 5.5% for RCP 2.6, 4.5 and 6.0 and they are lower by 7.9% for RCP 8.5. We analysed the impacts of greenhouse gases (GHGs) and ozone-depleting substances (ODSs) on UVI from 1960 by comparing CCMI sensitivity simulations (1960-2100) with fixed GHGs or ODSs at their respective 1960 levels. As expected with ODS fixed at their 1960 levels, there is no large decrease in ozone levels and consequently no sudden increase in UVI levels. With fixed GHG, we observed a delayed return of ozone to 1960 values, with a corresponding pattern of change observed on UVI, and looking at the UVI difference between 2090s values and 1960s values, we found an 8 % increase in the tropical belt during the summer of each hemisphere. Finally we show that, while in the Southern Hemisphere the UVI is mainly driven by total ozone column, in the Northern Hemisphere both total ozone column and aerosol optical depth drive UVI levels, with aerosol optical depth having twice as much influence on the UVI as total ozone column does.

Tropospheric ozone in CCMI models and Gaussian process emulation to understand biases in the SOCOLv3 chemistry-climate model.

Previous multi-model intercomparisons have shown that chemistry-climate models exhibit significant biases in tropospheric ozone compared with observations. We investigate annual-mean tropospheric column ozone in 15 models participating in the SPARC/IGAC (Stratosphere-troposphere Processes and their Role in Climate/International Global Atmospheric Chemistry) Chemistry-Climate Model Initiative (CCMI). These models exhibit a positive bias, on average, of up to 40-50% in the Northern Hemisphere compared with observations derived from the Ozone Monitoring Instrument and Microwave Limb Sounder (OMI/MLS), and a negative bias of up to ~30% in the Southern Hemisphere. SOCOLv3.0 (version 3 of the Solar-Climate Ozone Links CCM), which participated in CCMI, simulates global-mean tropospheric ozone columns of 40.2 DU - approximately 33% larger than the CCMI multi-model mean. Here we introduce an updated version of SOCOLv3.0, "SOCOLv3.1", which includes an improved treatment of ozone sink processes, and results in a reduction in the tropospheric column ozone bias of up to 8 DU, mostly due to the inclusion of N2O5 hydrolysis on tropospheric aerosols. As a result of these developments, tropospheric column ozone amounts simulated by SOCOLv3.1 are comparable with several other CCMI models. We apply Gaussian process emulation and sensitivity analysis to understand the remaining ozone bias in SOCOLv3.1. This shows that ozone precursors (nitrogen oxides (NOx), carbon monoxide, methane and other volatile organic compounds) are responsible for more than 90% of the variance in tropospheric ozone. However, it may not be the emissions inventories themselves that result in the bias, but how the emissions are handled in SOCOLv3.1, and we discuss this in the wider context of the other CCMI models. Given that the emissions data set to be used for phase 6 of the Coupled Model Intercomparison Project includes approximately 20% more NOx than the data set used for CCMI, further work is urgently needed to address the challenges of simulating sub-grid processes of importance to tropospheric ozone in the current generation of chemistry-climate models.

No Robust Evidence of Future Changes in Major Stratospheric Sudden Warmings: A Multi-model Assessment from CCMI.

Major stratospheric sudden warmings (SSWs) are the largest instance of wintertime variability in the Arctic stratosphere. Due to their relevance for the troposphere-stratosphere system, several previous studies have focused on their potential response to anthropogenic forcings. However, a wide range of results have been reported, from a future increase in the frequency of SSWs to a decrease. Several factors might explain these contradictory results, notably the use of different metrics for the identification of SSWs, and the impact of large climatological biases in single-model studies. Here we revisit the question of future SSWs changes, using an identical set of metrics applied consistently across 12 different models participating in the Chemistry Climate Model Initiative. From analyzing future integrations we find no statistically significant change in the frequency of SSWs over the 21st century, irrespective of the metric used for the identification of SSWs. Changes in other SSWs characteristics, such as their duration and the tropospheric forcing, are also assessed: again, we find no evidence of future changes over the 21st century.

Stratospheric ozone loss over the Eurasian continent induced by the polar vortex shift.

The Montreal Protocol has succeeded in limiting major ozone-depleting substance emissions, and consequently stratospheric ozone concentrations are expected to recover this century. However, there is a large uncertainty in the rate of regional ozone recovery in the Northern Hemisphere. Here we identify a Eurasia-North America dipole mode in the total column ozone over the Northern Hemisphere, showing negative and positive total column ozone anomaly centres over Eurasia and North America, respectively. The positive trend of this mode explains an enhanced total column ozone decline over the Eurasian continent in the past three decades, which is closely related to the polar vortex shift towards Eurasia. Multiple chemistry-climate-model simulations indicate that the positive Eurasia-North America dipole trend in late winter is likely to continue in the near future. Our findings suggest that the anticipated ozone recovery in late winter will be sensitive not only to the ozone-depleting substance decline but also to the polar vortex changes, and could be substantially delayed in some regions of the Northern Hemisphere extratropics.

Atmospheric impacts of the strongest known solar particle storm of 775 AD.

Sporadic solar energetic particle (SEP) events affect the Earth's atmosphere and environment, in particular leading to depletion of the protective ozone layer in the Earth's atmosphere, and pose potential technological and even life hazards. The greatest SEP storm known for the last 11 millennia (the Holocene) occurred in 774-775 AD, serving as a likely worst-case scenario being 40-50 times stronger than any directly observed one. Here we present a systematic analysis of the impact such an extreme event can have on the Earth's atmosphere. Using state-of-the-art cosmic ray cascade and chemistry-climate models, we successfully reproduce the observed variability of cosmogenic isotope 10Be, around 775 AD, in four ice cores from Greenland and Antarctica, thereby validating the models in the assessment of this event. We add to prior conclusions that any nitrate deposition signal from SEP events remains too weak to be detected in ice cores by showing that, even for such an extreme solar storm and sub-annual data resolution, the nitrate deposition signal is indistinguishable from the seasonal cycle. We show that such a severe event is able to perturb the polar stratosphere for at least one year, leading to regional changes in the surface temperature during northern hemisphere winters.

Deriving Global OH Abundance and Atmospheric Lifetimes for Long-Lived Gases: A Search for CH3CCl3 Alternatives.

An accurate estimate of global hydroxyl radical (OH) abundance is important for projections of air quality, climate, and stratospheric ozone recovery. As the atmospheric mixing ratios of methyl chloroform (CH3CCl3) (MCF), the commonly used OH reference gas, approaches zero, it is important to find alternative approaches to infer atmospheric OH abundance and variability. The lack of global bottom-up emission inventories is the primary obstacle in choosing a MCF alternative. We illustrate that global emissions of long-lived trace gases can be inferred from their observed mixing ratio differences between the Northern Hemisphere (NH) and Southern Hemisphere (SH), given realistic estimates of their NH-SH exchange time, the emission partitioning between the two hemispheres, and the NH versus SH OH abundance ratio. Using the observed long-term trend and emissions derived from the measured hemispheric gradient, the combination of HFC-32 (CH2F2), HFC-134a (CH2FCF3, HFC-152a (CH3CHF2), and HCFC-22 (CHClF2), instead of a single gas, will be useful as a MCF alternative to infer global and hemispheric OH abundance and trace gas lifetimes. The primary assumption on which this multispecies approach relies is that the OH lifetimes can be estimated by scaling the thermal reaction rates of a reference gas at 272 K on global and hemispheric scales. Thus, the derived hemispheric and global OH estimates are forced to reconcile the observed trends and gradient for all four compounds simultaneously. However, currently, observations of these gases from the surface networks do not provide more accurate OH abundance estimate than that from MCF.

Formaldehyde in the Tropical Western Pacific: Chemical sources and sinks, convective transport, and representation in CAM-Chem and the CCMI models.

Formaldehyde (HCHO) directly affects the atmospheric oxidative capacity through its effects on HOx. In remote marine environments, such as the Tropical Western Pacific (TWP), it is particularly important to understand the processes controlling the abundance of HCHO because model output from these regions is used to correct satellite retrievals of HCHO. Here, we have used observations from the CONTRAST field campaign, conducted during January and February 2014, to evaluate our understanding of the processes controlling the distribution of HCHO in the TWP as well as its representation in chemical transport/climate models. Observed HCHO mixing ratios varied from ~500 pptv near the surface to ~75 pptv in the upper troposphere. Recent convective transport of near surface HCHO and its precursors, acetaldehyde and possibly methyl hydroperoxide, increased upper tropospheric HCHO mixing ratios by ~33% (22 pptv); this air contained roughly 60% less NO than more aged air. Output from the CAM-Chem chemistry transport model (2014 meteorology) as well as nine chemistry climate models from the Chemistry-Climate Model Initiative (free-running meteorology) are found to uniformly underestimate HCHO columns derived from in situ observations by between 4 and 50%. This underestimate of HCHO likely results from a near factor of two underestimate of NO in most models, which strongly suggests errors in NOx emissions inventories and/or in the model chemical mechanisms. Likewise, the lack of oceanic acetaldehyde emissions and potential errors in the model acetaldehyde chemistry lead to additional underestimates in modeled HCHO of up to 75 pptv (~15%) in the lower troposphere.