NASA Satellite Measurements Show Global-Scale Reductions in Free Tropospheric Ozone in 2020 and Again in 2021 During COVID-19.

NASA satellite measurements show that ozone reductions throughout the Northern Hemisphere (NH) free troposphere reported for spring-summer 2020 during the COronaVIrus Disease 2019 pandemic have occurred again in spring-summer 2021. The satellite measurements show that tropospheric column ozone (TCO) (mostly representative of the free troposphere) for 20°N-60°N during spring-summer for both 2020 and 2021 averaged ∼3 Dobson Units (DU) (or ∼7%-8%) below normal. These ozone reductions in 2020 and 2021 were the lowest in the 2005-2021 record. We also include satellite measurements of tropospheric NO2 that exhibit reductions of ∼10%-20% in the NH in early spring-to-summer 2020 and 2021, suggesting that reduced pollution was the main cause for the low anomalies in NH TCO in 2020 and 2021. Reductions of TCO ∼2 DU (7%) are also measured in the Southern Hemisphere in austral summer but are not associated with reduced NO2.

Reanalysis intercomparisons of stratospheric polar processing diagnostics.

We compare herein polar processing diagnostics derived from the four most recent full-input reanalysis datasets: the National Centers for Environmental Prediction Climate Forecast System Reanalysis / Climate Forecast System, version 2 (CFSR/CFSv2), the European Centre for Medium-Range Weather Forecasts Interim Reanalysis (ERA-Interim), the Japanese Meteorological Agency's Japanese 55-year Reanalysis (JRA-55), and the National Aeronautics and Space Administration's Modern Era Retrospective-analysis for Research and Applications version 2 (MERRA-2). We focus on diagnostics based on temperatures and potential vorticity (PV) in the lower to middle stratosphere that are related to formation of polar stratospheric clouds (PSCs), chlorine activation, and the strength, size, and longevity of the stratospheric polar vortex. Polar minimum temperatures (Tmin) and the area of regions having temperatures below PSC formation thresholds (APSC) show large persistent differences between the reanalyses, especially in the southern hemisphere (SH), for years prior to 1999. Average absolute differences of the reanalyses from the reanalysis ensemble mean (REM) in Tmin are as large as 3 K at some levels in the SH (1.5 K in the NH), and absolute differences of reanalysis APSC from the REM up to 1.5% of a hemisphere (0.75% of a hemisphere in the NH). After 1999, the reanalyses converge toward better agreement in both hemispheres, dramatically so in the SH: Average Tmin differences from the REM are generally less than 1 K in both hemispheres, and average APSC differences less than 0.3% of a hemisphere. The comparisons of diagnostics based on isentropic PV for assessing polar vortex characteristics, including maximum PV gradients (MPVG) and the area of the vortex in sunlight (or sunlit vortex area, SVA), show more complex behavior: SH MPVG showed convergence toward better agreement with the REM after 1999, while NH MPVG differences remained largely constant over time; differences in SVA remained relatively constant in both hemispheres. While the average differences from the REM are generally small for these vortex diagnostics, understanding such differences among the reanalyses is complicated by the need to use different methods to obtain vertically-resolved PV for the different reanalyses. We also evaluated other winter season summary diagnostics, including the winter mean volume of air below PSC thresholds, and vortex decay dates. For the volume of air below PSC thresholds, the reanalyses generally agree best in the SH, where relatively small interannual variability has led to many winter seasons with similar polar processing potential and duration, and thus low sensitivity to differences in meteorological conditions among the reanalyses. In contrast, the large interannual variability of NH winters has given rise to many seasons with marginal conditions that are more sensitive to reanalysis differences. For vortex decay dates, larger differences are seen in the SH than in the NH; in general the differences in decay dates among the reanalyses follow from persistent differences in their vortex areas. Our results indicate that the transition from the reanalyses assimilating Tiros Operational Vertical Sounder (TOVS) data to Advanced TOVS and other data around 1998 - 2000 resulted in a profound improvement in the agreement of the temperature diagnostics presented (especially in the SH) and to a lesser extent the agreement of the vortex diagnostics. We present several recommendations for using reanalyses in polar processing studies, particularly related to the sensitivity to changes in data inputs and assimilation. Because of these sensitivities, we urge great caution for studies aiming to assess trends derived from reanalysis temperatures. We also argue that one of the best ways to assess the sensitivity of scientific results on polar processing is to use multiple reanalysis datasets.

Recent decline in extratropical lower stratospheric ozone attributed to circulation changes.

1998-2016 ozone trends in the lower stratosphere (LS) are examined using the Modern-Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2) and related NASA products. After removing biases resulting from step-changes in the MERRA-2 ozone observations, a discernible negative trend of -1.67±0.54 Dobson units per decade (DU/decade) is found in the 10-km layer above the tropopause between 20°N and 60°N. A weaker but statistically significant trend of -1.17±0.33 DU/decade exists between 50°S and 20°S. In the Tropics, a positive trend is seen in a 5-km layer above the tropopause. Analysis of an idealized tracer in a model simulation constrained by MERRA-2 meteorological fields provides strong evidence that these trends are driven by enhanced isentropic transport between the tropical (20°S-20°N) and extratropical LS in the past two decades. This is the first time that a reanalysis dataset has been used to detect and attribute trends in lower stratospheric ozone.

Chemical Mechanisms and Their Applications in the Goddard Earth Observing System (GEOS) Earth System Model.

NASA's Goddard Earth Observing System (GEOS) Earth System Model (ESM) is a modular, general circulation model (GCM), and data assimilation system (DAS) that is used to simulate and study the coupled dynamics, physics, chemistry, and biology of our planet. GEOS is developed by the Global Modeling and Assimilation Office (GMAO) at NASA Goddard Space Flight Center. It generates near-real-time analyzed data products, reanalyses, and weather and seasonal forecasts to support research targeted to understanding interactions among Earth System processes. For chemistry, our efforts are focused on ozone and its influence on the state of the atmosphere and oceans, and on trace gas data assimilation and global forecasting at mesoscale discretization. Several chemistry and aerosol modules are coupled to the GCM, which enables GEOS to address topics pertinent to NASA's Earth Science Mission. This paper describes the atmospheric chemistry components of GEOS and provides an overview of its Earth System Modeling Framework (ESMF)-based software infrastructure, which promotes a rich spectrum of feedbacks that influence circulation and climate, and impact human and ecosystem health. We detail how GEOS allows model users to select chemical mechanisms and emission scenarios at run time, establish the extent to which the aerosol and chemical components communicate, and decide whether either or both influence the radiative transfer calculations. A variety of resolutions facilitates research on spatial and temporal scales relevant to problems ranging from hourly changes in air quality to trace gas trends in a changing climate. Samples of recent GEOS chemistry applications are provided.

The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2).

The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) is the latest atmospheric reanalysis of the modern satellite era produced by NASA's Global Modeling and Assimilation Office (GMAO). MERRA-2 assimilates observation types not available to its predecessor, MERRA, and includes updates to the Goddard Earth Observing System (GEOS) model and analysis scheme so as to provide a viable ongoing climate analysis beyond MERRA's terminus. While addressing known limitations of MERRA, MERRA-2 is also intended to be a development milestone for a future integrated Earth system analysis (IESA) currently under development at GMAO. This paper provides an overview of the MERRA-2 system and various performance metrics. Among the advances in MERRA-2 relevant to IESA are the assimilation of aerosol observations, several improvements to the representation of the stratosphere including ozone, and improved representations of cryospheric processes. Other improvements in the quality of MERRA-2 compared with MERRA include the reduction of some spurious trends and jumps related to changes in the observing system, and reduced biases and imbalances in aspects of the water cycle. Remaining deficiencies are also identified. Production of MERRA-2 began in June 2014 in four processing streams, and converged to a single near-real time stream in mid 2015. MERRA-2 products are accessible online through the NASA Goddard Earth Sciences Data Information Services Center (GES DISC).

Reanalysis comparisons of upper tropospheric/lower stratospheric jets and multiple tropopauses.

The representation of upper tropospheric/lower stratospheric (UTLS) jet and tropopause characteristics is compared in five modern high-resolution reanalyses for 1980 through 2014. Climatologies of upper tropospheric jet, subvortex jet (the lowermost part of the stratospheric vortex), and multiple tropopause frequency distributions in MERRA (Modern Era Retrospective Analysis for Research and Applications), ERA-I (the ECMWF interim reanalysis), JRA-55 (the Japanese 55-year Reanalysis), and CFSR (the Climate Forecast System Reanalysis) are compared with those in MERRA-2. Differences between alternate products from individual reanalysis systems are assessed; in particular, a comparison of CFSR data on model and pressure levels highlights the importance of vertical grid spacing. Most of the differences in distributions of UTLS jets and multiple tropopauses are consistent with the differences in assimilation model grids and resolution: For example, ERA-I (with coarsest native horizontal resolution) typically shows a significant low bias in upper tropospheric jets with respect to MERRA-2, and JRA-55 a more modest one, while CFSR (with finest native horizontal resolution) shows a high bias with respect to MERRA-2 in both upper tropospheric jets and multiple tropopauses. Vertical temperature structure and grid spacing are especially important for multiple tropopause characterization. Substantial differences between MERRA and MERRA-2 are seen in mid- to high-latitude southern hemisphere winter upper tropospheric jets and multiple tropopauses, and in the upper tropospheric jets associated with tropical circulations during the solstice seasons; some of the largest differences from the other reanalyses are seen in the same times and places. Very good qualitative agreement among the reanalyses is seen between the large scale climatological features in UTLS jet and multiple tropopause distributions. Quantitative differences may, however, have important consequences for transport and variability studies. Our results highlight the importance of considering reanalyses differences in UTLS studies, especially in relation to resolution and model grids; this is particularly critical when using high-resolution reanalyses as an observational reference for evaluating global chemistry climate models.

Assessment of upper tropospheric and stratospheric water vapor and ozone in reanalyses as part of S-RIP.

Reanalysis data sets are widely used to understand atmospheric processes and past variability, and are often used to stand in as "observations" for comparisons with climate model output. Because of the central role of water vapor (WV) and ozone (O3) in climate change, it is important to understand how accurately and consistently these species are represented in existing global reanalyses. In this paper, we present the results of WV and O3 intercomparisons that have been performed as part of the SPARC (Stratosphere-troposphere Processes and their Role in Climate) Reanalysis Intercomparison Project (S-RIP). The comparisons cover a range of timescales and evaluate both inter-reanalysis and observation-reanalysis differences. We also provide a systematic documentation of the treatment of WV and O3 in current reanalyses to aid future research and guide the interpretation of differences amongst reanalysis fields. The assimilation of total column ozone (TCO) observations in newer reanalyses results in realistic representations of TCO in reanalyses except when data coverage is lacking, such as during polar night. The vertical distribution of ozone is also relatively well represented in the stratosphere in reanalyses, particularly given the relatively weak constraints on ozone vertical structure provided by most assimilated observations and the simplistic representations of ozone photochemical processes in most of the reanalysis forecast models. However, significant biases in the vertical distribution of ozone are found in the upper troposphere and lower stratosphere in all reanalyses. In contrast to O3, reanalysis estimates of stratospheric WV are not directly constrained by assimilated data. Observations of atmospheric humidity are typically used only in the troposphere, below a specified vertical level at or near the tropopause. The fidelity of reanalysis stratospheric WV products is therefore mainly dependent on the reanalyses' representation of the physical drivers that influence stratospheric WV, such as temperatures in the tropical tropopause layer, methane oxidation, and the stratospheric overturning circulation. The lack of assimilated observations and known deficiencies in the representation of stratospheric transport in reanalyses result in much poorer agreement amongst observational and reanalysis estimates of stratospheric WV. Hence, stratospheric WV products from the current generation of reanalyses should generally not be used in scientific studies.

Evaluation of the Ozone Fields in NASA's MERRA-2 Reanalysis.

We describe and assess the quality of the assimilated ozone product from the Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) produced at NASA's Global Modeling and Assimilation Office (GMAO) spanning the time period from 1980 to present. MERRA-2 assimilates partial column ozone retrievals from a series of Solar Backscatter Ultraviolet (SBUV) radiometers on NASA and NOAA spacecraft between January 1980 and September 2004; starting in October 2004 retrieved ozone profiles from the Microwave Limb Sounder (MLS) and total column ozone from the Ozone Monitoring Instrument on NASA's EOS Aura satellite are assimilated. We compare the MERRA-2 ozone with independent satellite and ozonesonde data focusing on the representation of the spatial and temporal variability of stratospheric and upper tropospheric ozone and on implications of the change in the observing system from SBUV to EOS Aura. The comparisons show agreement within 10 % (standard deviation of the difference) between MERRA-2 profiles and independent satellite data in most of the stratosphere. The agreement improves after 2004 when EOS Aura data are assimilated. The standard deviation of the differences between the lower stratospheric and upper tropospheric MERRA-2 ozone and ozonesondes is 11.2 % and 24.5 %, respectively, with correlations of 0.8 and above, indicative of a realistic representation of the near-tropopause ozone variability in MERRA-2. The agreement improves significantly in the EOS Aura period, however MERRA-2 is biased low in the upper troposphere with respect to the ozonesondes. Caution is recommended when using MERRA-2 ozone for decadal changes and trend studies.

Multi-year composite view of ozone enhancements and stratosphere-to-troposphere transport in dry intrusions of northern hemisphere extratropical cyclones.

We examine the role of extratropical cyclones in stratosphere-to-troposphere (STT) exchange with cyclone-centric composites of O3 retrievals from the Microwave Limb Sounder (MLS) and the Tropospheric Emission Spectrometer (TES), contrasting them to composites obtained with the Modern-Era Retrospective-analysis for Research and Applications (MERRA and MERRA-2) reanalyses and the GEOS-Chem chemical transport model. We identify 15,978 extratropical cyclones in the northern hemisphere (NH) for 2005-2012. The lowermost stratosphere (261 hPa) and middle troposphere (424 hPa) composites feature a 1,000 km-wide O3 enhancement in the dry intrusion (DI) airstream to the southwest of the cyclone center, coinciding with a lowered tropopause, enhanced potential vorticity, and decreased H2O. MLS composites at 261 hPa show that the DI O3 enhancements reach a 210 ppbv maximum in April. At 424 hPa, TES composites display maximum O3 enhancements of 27 ppbv in May. The magnitude and seasonality of these enhancements are captured by MERRA and MERRA-2, but GEOS-Chem is a factor of two too low. The MERRA-2 composites show that the O3-rich DI forms a vertically aligned structure between 300 and 800 hPa, wrapping cyclonically with the warm conveyor belt. In winter and spring DIs, O3 is enhanced by 100 ppbv or 100-130% at 300 hPa, with significant enhancements below 500 hPa (6-20 ppbv or 15-30%). We estimate that extratropical cyclones result in a STT flux of 119±56 Tg O3 yr-1, accounting for 42±20 % of the NH extratropical O3 STT flux. The STT flux in cyclones displays a strong dependence on westerly 300 hPa wind speeds.

Frequency and Impact of Summertime Stratospheric Intrusions over Maryland during DISCOVER-AQ (2011): New Evidence from NASA's GEOS-5 Simulations.

Aircraft observations and ozonesonde profiles collected on July 14 and 27, 2011, during the Maryland month-long DISCOVER-AQ campaign, indicate the presence of stratospheric air just above the planetary boundary layer (PBL). This raises the question of whether summer stratospheric intrusions (SIs) elevate surface ozone levels and to what degree they influence background ozone levels and contribute to ozone production. We used idealized stratospheric air tracers, along with observations, to determine the frequency and extent of SIs in Maryland during July 2011. On 4 of 14 flight days, SIs were detected in layers that the aircraft encountered above the PBL from the coincidence of enhanced ozone, moderate CO, and low moisture. Satellite observations of lower tropospheric humidity confirmed the occurrence of synoptic scale influence of SIs as do simulations with the GEOS-5 Atmospheric General Circulation Model. The evolution of GEOS-5 stratospheric air tracers agree with the timing and location of observed stratospheric influence and indicate that more than 50% of air in SI layers above the PBL had resided in the stratosphere within the previous 14 days. Despite having a strong influence in the lower free troposphere, these events did not significantly affect surface ozone, which remained low on intrusion days. The model indicates similar frequencies of stratospheric influence during all summers from 2009-2013. GEOS-5 results suggest that, over Maryland, the strong inversion capping the summer PBL limits downward mixing of stratospheric air during much of the day, helping to preserve low surface ozone associated with frontal passages that precede SIs.

Structure and Dynamics of the Quasi-Biennial Oscillation in MERRA-2.

The structure, dynamics, and ozone signal of the Quasi-Biennial Oscillation produced by the 35-year NASA MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications) reanalysis are examined based on monthly mean output. Along with the analysis of the QBO in assimilation winds and ozone, the QBO forcings created by assimilated observations, dynamics, parameterized gravity wave drag, and ozone chemistry parameterization are examined and compared with the original MERRA system. Results show that the MERRA-2 reanalysis produces a realistic QBO in the zonal winds, mean meridional circulation, and ozone over the 1980-2015 time period. In particular, the MERRA-2 zonal winds show improved representation of the QBO 50 hPa westerly phase amplitude at Singapore when compared to MERRA. The use of limb ozone observations creates improved vertical structure and realistic downward propagation of the ozone QBO signal during times when the MLS ozone limb observations are available (October 2004 to present). The increased equatorial GWD in MERRA-2 has reduced the zonal wind data analysis contribution compared to MERRA so that the QBO mean meridional circulation can be expected to be more physically forced and therefore more physically consistent. This can be important for applications in which MERRA-2 winds are used to drive transport experiments.