Relative Effects of the Greenhouse Gases and Stratospheric Ozone Increases on Temperature and Circulation in the Stratosphere over the Arctic

Using a stratosphere-resolving general circulation model, the relative effects of stratospheric ozone and greenhouse gases (GHGs) increase on the temperature and circulation in the Arctic stratosphere are examined. Results show that stratospheric ozone or GHGs increase alone could result in a cooling and strengthening extratropical stratosphere during February, March and April. However, the contribution of stratospheric ozone increases alone on the cooling and strengthening Arctic stratosphere is approximately 2 fold that of the GHGs increase alone. Model simulations suggested that the larger responses of the Arctic stratosphere to the ozone increase alone are closely related to the wave fluxes in the stratosphere, rather than the wave activity in the stratosphere. In response to the ozone increase, the vertical propagation of planetary waves from the troposphere into the mid-latitude stratosphere weakens, mainly contributed by its wavenumber-1 component. The impeded planetary waves tend to result from the larger zonal wind shear and vertical gradient of the buoyancy frequency. The magnitudes of anomalies in the zonal wind shear and buoyancy frequency in response to GHGs increase alone are smaller than in response to the ozone increase, which is in accordance with the larger contribution of stratospheric ozone to the temperature and circulation in the Arctic stratosphere.

Driving Influences of the Doppler Flash Observed by SuperDARN HF Radars in Response to Solar Flares

Sudden enhancement in high-frequency absorption is a well-known impact of solar flare-driven Short-Wave Fadeout (SWF). Less understood, is a perturbation of the radio wave frequency as it traverses the ionosphere in the early stages of SWF, also known as the Doppler flash. Investigations have suggested two possible sources that might contribute to it's manifestation: first, enhancements of plasma density in the D-and lower E-regions; second, the lowering of the F-region reflection point. Our recent work investigated a solar flare event using first principles modeling and Super Dual Auroral Radar Network (SuperDARN) HF radar observations and found that change in the F-region refractive index is the primary driver of the Doppler flash. This study analyzes multiple solar flare events observed across different SuperDARN HF radars to determine how flare characteristics, properties of the traveling radio wave, and geophysical conditions impact the Doppler flash. In addition, we use incoherent scatter radar data and first-principles modeling to investigate physical mechanisms that drive the lowering of the F-region reflection points. We found, (a) on average, the change in E- and F-region refractive index is the primary driver of the Doppler flash, (b) solar zenith angle, ray's elevation angle, operating frequency, and location of the solar flare on the solar disk can alter the ionospheric regions of maximum contribution to the Doppler flash, (c) increased ionospheric Hall and Pedersen conductance causes a reduction of the daytime eastward electric field, and consequently reduces the vertical ion-drift in the lower and middle latitude ionosphere, which results in lowering of the F-region ray reflection point.

Impact of Increased Vertical Resolution in WACCM on the Climatology of Major Sudden Stratospheric Warmings

Sudden stratospheric warmings (SSWs) are a major mode of variability of the winter stratosphere. In recent years, climate models have improved their ability to simulate SSWs. However, the representation of the frequency and temporal distribution of SSWs in models depends on many factors and remains challenging. The vertical resolution of a model might be one such factor. Therefore, here we analyse the impact of increased vertical resolution on the simulation of major sudden stratospheric warmings (SSWs) in the Whole Atmosphere Community Climate Model (WACCM). We compare two versions of the model, WACCM3.5 and WACCM4. We find that the frequency of occurrence of SSWs is improved in the newer version and closer to that obtained using reanalysis. Furthermore, simulations with a coupled ocean best reproduce the behaviour of temperature during these events. Increasing vertical resolution increases the number of occurrences; however, it does not produce significantly different results than standard resolution. WACCM4 also does not reproduce vortex split events well, generating far fewer of these than observed. Finally, the ratio between polar vortex splits and displacement events in the model is slightly better for non-ocean-coupled simulations. We conclude that, at least for WACCM4, the use of the high vertical resolution configuration is not cost-effective for the study of SSWs.

Intra-Annual Variation of Eddy Diffusion (k zz ) in the MLT, From SABER and SCIAMACHY Atomic Oxygen Climatologies

Atomic oxygen (O) in the mesosphere and lower thermosphere (MLT) results from a balance between production via photo-dissociation in the lower thermosphere and chemical loss by recombination in the upper mesosphere. The transport of O downward from the lower thermosphere into the mesosphere is preferentially driven by the eddy diffusion process that results from dissipating gravity waves and instabilities. The motivation here is to probe the intra-annual variability of the eddy diffusion coefficient (k _zz ) and eddy velocity in the MLT based on the climatology of the region, initially accomplished by Garcia and Solomon (1985, https://doi.org/10.1029/JD090iD02p03850). In the current study, the intra-annual cycle was divided into 26 two-week periods for each of three zones: the northern hemisphere (NH), southern hemisphere (SH), and equatorial (EQ). Both 16 years of SABER (2002-2018) and 10 years of SCIAMACHY (2002-2012) O density measurements, along with NRLMSIS^® 2.0 were used for calculation of atomic oxygen eddy diffusion velocities and fluxes. Our prominent findings include a dominant annual oscillation below 87 km in the NH and SH zones, with a factor of 3-4 variation between winter and summer at 83 km, and a dominant semiannual oscillation at all altitudes in the EQ zone. The measured global average k _zz at 96 km lacks the intra-annual variability of upper atmosphere density data deduced by Qian et al. (2009, https://doi.org/10.1029/2008JA013643). The very large seasonal (and hemispherical) variations in k _zz and O densities are important to separate and isolate in satellite analysis and to incorporate in MLT models.

Transport and Variability of Tropospheric Ozone over Oceania and Southern Pacific during the 2019–20 Australian Bushfires

The present study contributes to the scientific effort for a better understanding of the potential of the Australian biomass burning events to influence tropospheric trace gas abundances at the regional scale. In order to exclude the influence of the long-range transport of ozone precursors from biomass burning plumes originating from Southern America and Africa, the analysis of the Australian smoke plume has been driven over the period December 2019 to January 2020. This study uses satellite (IASI, MLS, MODIS, CALIOP) and ground-based (sun-photometer, FTIR, ozone radiosondes) observations. The highest values of aerosol optical depth (AOD) and carbon monoxide total columns are observed over Southern and Central Australia. Transport is responsible for the spatial and temporal distributions of aerosols and carbon monoxide over Australia, and also the transport of the smoke plume outside the continent. The dispersion of the tropospheric smoke plume over Oceania and Southern Pacific extends from tropical to extratropical latitudes. Ozone radiosonde measurements performed at Samoa (14.4°S, 170.6°W) and Lauder (45.0°S, 169.4°E) indicate an increase in mid-tropospheric ozone (6–9 km) (from 10% to 43%) linked to the Australian biomass burning plume. This increase in mid-tropospheric ozone induced by the transport of the smoke plume was found to be consistent with MLS observations over the tropical and extratropical latitudes. The smoke plume over the Southern Pacific was organized as a stretchable anticyclonic rolling which impacted the ozone variability in the tropical and subtropical upper-troposphere over Oceania. This is corroborated by the ozone profile measurements at Samoa which exhibit an enhanced ozone layer (29%) in the upper-troposphere. Our results suggest that the transport of Australian biomass burning plumes have significantly impacted the vertical distribution of ozone in the mid-troposphere southern tropical to extratropical latitudes during the 2019–20 extreme Australian bushfires.

Global climate disruption and regional climate shelters after the Toba supereruption

The Toba eruption ∼74,000 y ago was the largest volcanic eruption since the start of the Pleistocene and represents an important test case for understanding the effects of large explosive eruptions on climate and ecosystems. However, the magnitude and repercussions of climatic changes driven by the eruption are strongly debated. High-resolution paleoclimate and archaeological records from Africa find little evidence for the disruption of climate or human activity in the wake of the eruption in contrast with a controversial link with a bottleneck in human evolution and climate model simulations predicting strong volcanic cooling for up to a decade after a Toba-scale eruption. Here, we use a large ensemble of high-resolution Community Earth System Model (CESM1.3) simulations to reconcile climate model predictions with paleoclimate records, accounting for uncertainties in the magnitude of Toba sulfur emissions with high and low emission scenarios. We find a near-zero probability of annual mean surface temperature anomalies exceeding 4 °C in most of Africa in contrast with near 100% probabilities of cooling this severe in Asia and North America for the high sulfur emission case. The likelihood of strong decreases in precipitation is low in most of Africa. Therefore, even Toba sulfur release at the upper range of plausible estimates remains consistent with the muted response in Africa indicated by paleoclimate proxies. Our results provide a probabilistic view of the uneven patterns of volcanic climate disruption during a crucial interval in human evolution, with implications for understanding the range of environmental impacts from past and future supereruptions.

On the response of the middle atmosphere to anthropogenic forcing

Anthropogenic forcing of the atmosphere by greenhouse gases (GHG) and ozone-depleting substances has provided an unintended test of the robustness of current understanding of the physics and chemistry of the middle atmosphere, that is, the stratosphere and mesosphere. We explore this topic by examining how well anthropogenic changes can be simulated by modern, comprehensive numerical models. Specifically, we discuss the simulations of trends in global mean temperature; the development of the ozone hole and its impact on the dynamics of the Southern Hemisphere, both in the stratosphere and troposphere; trends in the stratospheric Brewer-Dobson circulation; and the response of the quasi-biennial oscillation (QBO) to increasing burdens of CO₂ . We find that, in most of these cases, numerical simulation is able to reproduce observed changes and provide physical insights into the relevant mechanisms. Simulation of the QBO is on a less firm footing. Although many numerical models can now generate realistic QBOs, future projections of its behavior under the increasing burdens of GHG are inconsistent and even contradictory.

Transport of Nitric Oxide Via Lagrangian Coherent Structures Into the Top of the Polar Vortex

The energetic particle precipitation (EPP) indirect effect (IE) refers to the downward transport of reactive odd nitrogen (NOx = NO + NO2) produced by EPP (EPP-NOx) from the polar winter mesosphere and lower thermosphere to the stratosphere where it can destroy ozone. Previous studies of the EPP IE examined NOx descent averaged over the polar region, but the work presented here considers longitudinal variations. We report that the January 2009 split Arctic vortex in the stratosphere left an imprint on the distribution of NO near the mesopause, and that the magnitude of EPP-NOx descent in the upper mesosphere depends strongly on the planetary wave (PW) phase. We focus on an 11-day case study in late January immediately following the 2009 sudden stratospheric warming during which regional-scale Lagrangian coherent structures (LCSs) formed atop the strengthening mesospheric vortex. The LCSs emerged over the north Atlantic in the vicinity of the trough of a 10-day westward traveling planetary wave. Over the next week, the LCSs acted to confine NO-rich air to polar latitudes, effectively prolonging its lifetime as it descended into the top of the polar vortex. Both a whole atmosphere data assimilation model and satellite observations show that the PW trough remained coincident in space and time with the NO-rich air as both migrated westward over the Canadian Arctic. Estimates of descent rates indicate five times stronger descent inside the PW trough compared to other longitudes. This case serves to set the stage for future climatological analysis of NO transport via LCSs.

The Diurnal Variation in Stratospheric Ozone from MACC Reanalysis, ERA-Interim, WACCM, and Earth Observation Data: Characteristics and Intercomparison

In this study, we compare the diurnal variation in stratospheric ozone of the MACC (Monitoring Atmospheric Composition and Climate) reanalysis, ECMWF Reanalysis Interim (ERA-Interim), and the free-running WACCM (Whole Atmosphere Community Climate Model). The diurnal variation of stratospheric ozone results from photochemical and dynamical processes depending on altitude, latitude, and season. MACC reanalysis and WACCM use similar chemistry modules and calculate a similar diurnal cycle in ozone when it is caused by a photochemical variation. The results of the two model systems are confirmed by observations of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) experiment and three selected sites of the Network for Detection of Atmospheric Composition Change (NDACC) at Mauna Loa, Hawaii (tropics), Bern, Switzerland (midlatitudes), and Ny-Ålesund, Svalbard (high latitudes). On the other hand, the ozone product of ERA-Interim shows considerably less diurnal variation due to photochemical variations. The global maxima of diurnal variation occur at high latitudes in summer, e.g., near the Arctic NDACC site at Ny-Ålesund, Svalbard. The local OZORAM radiometer observes this effect in good agreement with MACC reanalysis and WACCM. The sensed diurnal variation at Ny-Ålesund is up to 8% (0.4 ppmv) due to photochemical variations in summer and negligible during the dynamically dominated winter. However, when dynamics play a major role for the diurnal ozone variation as in the lower stratosphere (100–20 hPa), the reanalysis models ERA-Interim and MACC which assimilate data from radiosondes and satellites outperform the free-running WACCM. Such a domain is the Antarctic polar winter where a surprising novel feature of diurnal variation is indicated by MACC reanalysis and ERA-Interim at the edge of the polar vortex. This effect accounts for up to 8% (0.4 ppmv) in both model systems. In summary, MACC reanalysis provides a global description of the diurnal variation of stratospheric ozone caused by dynamics and photochemical variations. This is of high interest for ozone trend analysis and other research which is based on merged satellite data or measurements at different local time.

Persisting volcanic ash particles impact stratospheric SO2 lifetime and aerosol optical properties

Volcanic ash is often neglected in climate simulations because ash particles are assumed to have a short atmospheric lifetime, and to not participate in sulfur chemistry. After the Mt. Kelut eruption in 2014, stratospheric ash-rich aerosols were observed for months. Here we show that the persistence of super-micron ash is consistent with a density near 0.5 g cm-3, close to pumice. Ash-rich particles dominate the volcanic cloud optical properties for the first 60 days. We also find that the initial SO2 lifetime is determined by SO2 uptake on ash, rather than by reaction with OH as commonly assumed. About 43% more volcanic sulfur is removed from the stratosphere in 2 months with the SO2 heterogeneous chemistry on ash particles than without. This research suggests the need for re-evaluation of factors controlling SO2 lifetime in climate model simulations, and of the impact of volcanic ash on stratospheric chemistry and radiation.

Homogeneity of the Temperature Data Series from ERA5 and MERRA2 and Temperature Trends

The stratosphere and its dynamics are a very important part of atmospheric circulation. We need to analyze its climatology, as well as long-term trends. A long-term trend study needs homogenous datasets without significant artificial discontinuities. The analysis is based on the two newest released reanalyses, Modern Era-Retrospective Analysis (MERRA2) and European Center for Medium-Range Weather Forecast Reanalysis (ERA5). The aim of this study is to detect discontinuities in the temperature time series from the above reanalyses with the help of the Pettitt homogeneity test for pressure layers above 500 hPa up to 1 hPa in January and February, and show a comparison of temperature trends from the studied reanalyses and GPS radio occultation (GPS RO). We search for individual grid points where these discontinuities occur, and also for the years when they occur (geographical and temporal distribution). As expected, the study confirms better results for the Northern Hemisphere due to the denser data coverage. A high number of grid points with jumps on the Southern Hemisphere is found, especially at higher pressure levels (from 50 hPa). The spatial and vertical distribution of discontinuities is also presented. The vertical distribution reveals the reduction of the number of jumps around 10 hPa, especially for ERA5 reanalysis. The results show that ERA5 has significantly less jumps than MERRA2. We also study temperature trends from reanalyses and GPS RO and our analysis shows that the agreement between the reanalyses and observations are very good for the period 2006–2018.

Observation and Attribution of Temperature Trends Near the Stratopause From HALOE

This study considers time series of temperature versus pressure, T(p), from the Halogen Occultation Experiment (HALOE) across the stratopause region, where the effects of radiative forcings from the greenhouse gases (CO2 and H2O) and from ozone are most pronounced. Trend analyses of HALOE T(p) values for 1993-2005 are for six levels from 2.0 to 0.3 hPa with a vertical resolution of about 4 km and for eight latitude zones from 65°S to 65°N. The analyses account for the forcing effects from the 11-yr solar cycle. HALOE trends at 2.0 hPa are of the order of -1.0 K/decade across the tropics and subtropics but then become smaller (-0.5 K/decade) at the middle latitudes. Near-global T(p) trends are of order -0.5 K/decade but have a minimum of -0.2 K/decade at 1.0 hPa; they are clearly negative in the southern but slightly positive in the northern hemisphere. The combined radiative forcings from CO2, H2O, and ozone vary between -0.4 and -0.7 K/decade for 1993-2005 and are hemispherically symmetric. The HALOE temperature trend and total radiative cooling profiles differ from those reported from observations and calculations for 1980-2000, mainly because the ozone trends changed from clearly negative in the 1980s through mid-1990s to slightly positive during the time of HALOE. Trends at low latitudes for the tracer, methane (CH4), increase from 2% to 4%/decade from 50 to 10 hPa and then to ~6%/decade by 5 hPa. Analyses of time series of CH4 across the stratopause reveal subseasonal scale variability within the northern hemisphere that reduces the significance of the T(p) trends.

Observed responses of mesospheric water vapor to solar cycle and dynamical forcings

This study focuses on responses of mesospheric water vapor (H2O) to the solar cycle flux at Lyman-α wavelength and to wave forcings according to the multivariate ENSO index (MEI). The zonal-averaged responses are for latitudes from 60°S to 60°N and pressure-altitudes from 0.01 to 1.0 hPa, as obtained by multiple linear regression (MLR) analyses of time series of H2O from the Halogen Occultation Experiment (HALOE) for July 1992 to November 2005. The solar responses change from strong negative H2O values in the upper mesosphere to very weak, positive values in the tropical lower mesosphere. Those response profiles at the low latitudes agree reasonably with published results for H2O from the Microwave Limb Sounder (MLS). The distribution of seasonal H2O amplitudes corresponds well with that for temperature and is in accord with the seasonal net circulation. In general, the responses of H2O to MEI are anti-correlated with those of temperature. H2O responses to MEI are negative in the upper mesosphere and largest in the northern hemisphere; responses in the lower mesosphere are more symmetric with latitude. The H2O trends from MLR for the lower mesosphere agree with those reported from time series of microwave observations at two ground-based network stations.

Long-term climate change in the D-region

Controversy exists over the potential effects of long-term increases in greenhouse gas concentrations on the ionospheric D-region at 60–90 km altitudes. Techniques involving in-situ rocket measurements, remote optical observations, and radio wave reflection experiments have produced conflicting results. This study reports a novel technique that analyses long-distance subionospheric very low frequency radiowave observations of the NAA 24.0 kHz transmitter, Cutler, Maine, made from Halley Station, Antarctica, over the period 1971–2016. The analysis is insensitive to any changes in the output power of the transmitter, compensates for the use of different data logging equipment, and can confirm the accuracy of the timing systems operated over the 45 year long record. A ~10% reduction in the scale size of the transmitter nighttime interference fringe pattern has been determined, taking into account the quasi-11 year solar cycle. Subionospheric radiowave propagation modeling suggests that the contraction of the interference fringe pattern about the mid-latitude NAA transmitter is due to a 3 km reduction in the effective height of the nighttime ionospheric D-region over the last 45 years. This is consistent with the effect of enhanced infra-red cooling by increasing greenhouse gases.

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.

On transient climate change at the Cretaceous-Paleogene boundary due to atmospheric soot injections

Climate simulations that consider injection into the atmosphere of 15,000 Tg of soot, the amount estimated to be present at the Cretaceous-Paleogene boundary, produce what might have been one of the largest episodes of transient climate change in Earth history. The observed soot is believed to originate from global wildfires ignited after the impact of a 10-km-diameter asteroid on the Yucatán Peninsula 66 million y ago. Following injection into the atmosphere, the soot is heated by sunlight and lofted to great heights, resulting in a worldwide soot aerosol layer that lasts several years. As a result, little or no sunlight reaches the surface for over a year, such that photosynthesis is impossible and continents and oceans cool by as much as 28 °C and 11 °C, respectively. The absorption of light by the soot heats the upper atmosphere by hundreds of degrees. These high temperatures, together with a massive injection of water, which is a source of odd-hydrogen radicals, destroy the stratospheric ozone layer, such that Earth's surface receives high doses of UV radiation for about a year once the soot clears, five years after the impact. Temperatures remain above freezing in the oceans, coastal areas, and parts of the Tropics, but photosynthesis is severely inhibited for the first 1 y to 2 y, and freezing temperatures persist at middle latitudes for 3 y to 4 y. Refugia from these effects would have been very limited. The transient climate perturbation ends abruptly as the stratosphere cools and becomes supersaturated, causing rapid dehydration that removes all remaining soot via wet deposition.

Solar cycle response and long-term trends in the mesospheric metal layers

The meteoric metal layers (Na, Fe, and K)-which form as a result of the ablation of incoming meteors-act as unique tracers for chemical and dynamical processes that occur within the upper mesosphere/lower thermosphere region. In this work, we examine whether these metal layers are sensitive indicators of decadal long-term changes within the upper atmosphere. Output from a whole-atmosphere climate model is used to assess the response of the Na, K, and Fe layers across a 50 year period (1955-2005). At short timescales, the K layer has previously been shown to exhibit a very different seasonal behavior compared to the other metals. Here we show that this unusual behavior is also exhibited at longer timescales (both the ~11 year solar cycle and 50 year periods), where K displays a much more pronounced response to atmospheric temperature changes than either Na or Fe. The contrasting solar cycle behavior of the K and Na layers predicted by the model is confirmed using satellite and lidar observations for the period 2004-2013.

Equivalent latitude computation using regions of interest (ROI)

This paper introduces a novel algorithm to compute equivalent latitude by applying regions of interest (ROI). The technique is illustrated using code written in Interactive Data Language (IDL). The ROI method is compared with the "piecewise-constant" method, the approach commonly used in atmospheric sciences, using global fields of atmospheric potential vorticity. The ROI method is considerably more accurate and computationally faster than the piecewise-constant method, and it also works well with irregular grids. Both the ROI and piecewise-constant IDL codes for equivalent latitude are included as a useful reference for the research community.