Search

The ability to seasonally forecast the Quasi-Biennial Oscillation (QBO) was examined using NASA S2S (Sub-seasonal to Seasonal), 9-month, retrospective forecasts. Validation of these forecasts showed that the S2S retrospective QBO forecasts improved skill in predicting the QBO amplitude and phase over a simple QBO phase propagation model at forecast lead times of 1 to 3 months. Results from an initial assessment of whether more accurate QBO forecasts can improve Northern Hemisphere winter sea level pressure forecasts showed no significant forecast improvement at a 1-month lead time, indicating the need for improved stratosphere-troposphere QBO coupling metrics and pathway identification. Overall, these results suggest that future improvements in representing the QBO in global models can increase the ensemble fraction of valid 1 to 3 month QBO forecasts and potentially extend useful QBO forecasts beyond 3 months.

The eruption of Mount Pinatubo in June 1991 had dramatic effects on the global climate, introducing a sudden and extreme forcing on the radiative budget. In this presentation we analyze the effects of the Mt. Pinatubo eruption on the seasonal forecast performed with Goddard Earth Observing System (GEOS), an Earth System Model that includes a bulk aerosol model coupled to radiation and an interactive ocean.We performed three ten-member ensembles of 12-month simulations (June 1991-May 1992). These simulations were performed with a 0.5 longitude by 0.5 latitude horizontal resolution. Out of the three ensembles, one excludes the eruption, representing the control experiment. The other two ensembles include the eruption of Mt. Pinatubo, varying in the effective radius of the volcanic sulfate: one assumed the effective radius to be equal to 0.35 micron (as for tropospheric aerosol) and the other has the effective radius set to 0.6 micron (closer to natural observation). Through the analysis of the aerosol forcing results derived from the two assumptions, we can show how this forcing acts on the seasonal forecast system. In particular, we will focus on the impacts to the surface and ocean temperatures and precipitation.

The major stratospheric sudden warming in February 2018 occurred in unison with changes throughout the troposphere in the northern hemisphere. Changes in large-scale dynamics and regional weather patterns are described in this study, along with features in major features of the Earth system, such as sea ice. The focus is on analyses using the Goddard Earth Observing System (GEOS) analyses, weather prediction, and seasonal forecasts. Specific topics to discuss include: - The morphology of the sudden warming, including planetary waves in the troposphere and stratosphere.- Relationships with weather over North America. - Warming over the Arctic and anomalous sea-ice distributions. - Predictability on the one-to-four week timescales. The predictions will be based on routine products generated using the GEOS systems, as well as targeted experiments that focus on the roles of different mechanisms (such as ice-atmosphere interactions).

Project Loon has an overall goal of providing worldwide internet coverage using a network of long-durationsuper-pressure balloons. Since 2013, Loon has launched over 1600 balloons from multiple tropical and middlelatitude locations. These GPS tracked balloon trajectories provide lower stratospheric wind information overthe oceans and remote land areas where traditional radiosonde soundings are sparse, thus providing uniquecoverage of lower stratospheric winds. To fully investigate these Loon winds we: 1) compare the Loon windsto winds produced by a global data assimilation system (DAS: NASA GEOS) and 2) assimilate the Loon windsinto the same comprehensive DAS. Results show that in middle latitudes the Loon winds and DAS winds agreewell, and the Loon wind assimilation has only a minor impact on the forecasts. However, in the Tropics, thereis often a substantial difference between the assimilated winds and the observed Loon winds, of 8 m/s or morein magnitude. In these cases, assimilating the Loon winds significantly improves the meteorological analysesand subsequently the forecasts of the Loon winds. By highlighting cases where the Loon and DAS winds differ,these results can lead to improved understanding of stratospheric winds, especially in the tropics, as well asaiding analyses of the representation of dynamical forcing mechanisms in the GEOS model.

There occurred a unique and significant disruption of the Quasi-Biennial Oscillation (QBO) during the Northern Hemisphere winter of 2015-16. Here we document the return of the QBO to its normal downward phase speed and period based on Singapore soundings, MERRA-2 re-analysis (Modern Era Retrospective Reanalysis for Research and Applications), and a simple QBO model. Daily averaged zonal winds from 100-10 hPa are used to characterized the behavior of the QBO's amplitude and phase as seen in the first two EOFs (Empirical Orthogonal Functions). These EOFs capture the QBO structure and evolution during the pre-disruption, disruption, and post-disruption times. Results show that the amplitude and phase returned to normal by June 2016, however the post-disruption QBO phase was delayed relative to the pre-disruption phase by four tenths of a QBO cycle (~11 months). A rapid, seasonal, phase shift of this magnitude is shown to be unique in the QBO observational record.

The 1998-2016 ozone trends in the lower stratosphere are examined using the Modern-Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2) and related National Aeronautics and Space Administration 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 deg N and 60 deg N. A weaker but statistically significant trend of 1.17 +/- 0.33 DU/decade exists between 50 deg S and 20 deg 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 deg S-20 deg N) and extratropical lower stratosphere in the past two decades. This is the first time that a reanalysis data set has been used to detect and attribute trends in lower stratospheric ozone.

Long-term changes of ozone in the UTLS (Upper Troposphere / Lower Stratosphere) reflect the response to decreases in the stratospheric concentrations of ozone-depleting substances as well as changes in the stratospheric circulation induced by climate change. To date, studies of UTLS ozone changes and variability have relied mainly on satellite and in-situ observations as well as chemistry-climate model simulations. By comparison, the potential of reanalysis ozone data remains relatively untapped. This is despite evidence from recent studies, including detailed analyses conducted under SPARC (Scalable Processor Architecture) Reanalysis Intercomparison Project (S-RIP), that demonstrate that stratospheric ozone fields from modern atmospheric reanalyses exhibit good agreement with independent data while delineating issues related to inhomogeneities in the assimilated observations. In this presentation, we will explore the possibility of inferring long-term geographically and vertically resolved behavior of the lower stratospheric (LS) ozone from NASA's MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications -2) reanalysis after accounting for the few known discontinuities and gaps in its assimilated input data. This work builds upon previous studies that have documented excellent agreement between MERRA-2 ozone and ozonesonde observations in the LS. Of particular importance is a relatively good vertical resolution of MERRA-2 allowing precise separation of tropospheric and stratospheric ozone contents. We also compare the MERRA-2 LS ozone results with the recently completed 37-year simulation produced using Goddard Earth Observing System in "replay" mode coupled with the GMI (Global Modeling Initiative) chemistry mechanism. Replay mode dynamically constrains the model with the MERRA-2 reanalysis winds, temperature, and pressure. We will emphasize the areas of agreement of the reanalysis and replay and interpret differences between them in the context of our increasing understanding of model transport driven by assimilated winds.

Project Loon has an overall goal of providing worldwide internet coverage using a network of long-duration super-pressure balloons. Beginning in 2013, Loon has launched over 1600 balloons from multiple tropical and middle latitude locations. These GPS tracked balloon trajectories provide lower stratospheric wind information over the oceans and remote land areas where traditional radiosonde soundings are sparse, thus providing unique coverage of lower stratospheric winds. To fully investigate these Loon winds we: 1) compare the Loon winds to winds produced by a global data assimilation system (DAS: NASA GEOS) and 2) assimilate the Loon winds into the same comprehensive DAS. Results show that in middle latitudes the Loon winds and DAS winds agree well and assimilating the Loon winds have only a small impact on short-term forecasting of the Loon winds, however, in the tropics the loon winds and DAS winds often disagree substantially (8 m/s or more in magnitude) and in these cases assimilating the loon winds significantly improves the forecast of the loon winds. By highlighting cases where the Loon and DAS winds differ, these results can lead to improved understanding of stratospheric winds, especially in the tropics.

The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) is an atmospheric reanalysis, spanning 1980 through near-realtime, that uses state-of-the-art processing of observations from the continually evolving global observing system. The effectiveness of any reanalysis is a function not only of the input observations themselves, but also of how the observations are handled in the assimilation procedure. Relevant issues to consider include, but are not limited to, data selection, data preprocessing, quality control, bias correction procedures, and blacklisting. As the assimilation algorithm and earth system models are fundamentally fixed in a reanalysis, it is often a change in the character of the observations, and their feedbacks on the system, that cause changes in the character of the reanalysis. It is therefore important to provide documentation of the observing system so that its discontinuities and transitions can be readily linked to discontinuities seen in the gridded atmospheric fields of the reanalysis. With this in mind, this document provides an exhaustive list of the input observations, the context under which they are assimilated, and an initial assessment of selected core observations fundamental to the reanalysis.

This study investigates tropical waves and their role in driving a quasi-biennial oscillation (QBO)-like signal in stratospheric winds in a global 7-km-horizontal-resolution atmospheric general circulation model. The Nature Run (NR) is a 2-year global mesoscale simulation of the Goddard Earth Observing System Model, version 5 (GEOS-5). In the tropics, there is evidence that the NR supports a broad range of convectively generated waves. The NR precipitation spectrum resembles the observed spectrum in many aspects, including the preference for westward-propagating waves. However, even with very high horizontal resolution and a healthy population of resolved waves, the zonal force provided by the resolved waves is still too low in the QBO region and parameterized gravity wave drag is the main driver of the NR QBO-like oscillation (NRQBO). The authors suggest that causes include coarse vertical resolution and excessive dissipation. Nevertheless, the very-high-resolution NR provides an opportunity to analyze the resolved wave forcing of the NR-QBO. In agreement with previous studies, large-scale Kelvin and small-scale waves contribute to the NRQBO driving in eastward shear zones and small-scale waves dominate the NR-QBO driving in westward shear zones. Waves with zonal wavelength,1000 km account for up to half of the small-scale (,3300 km) resolved wave forcing in eastward shear zones and up to 70% of the small-scale resolved wave forcing in westward shear zones of the NR-QBO.

This report documents an evaluation by the Global Modeling and Assimilation Office (GMAO) of a two-year 7-km-resolution non-hydrostatic global mesoscale simulation produced with the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model. The simulation was produced as a Nature Run for conducting observing system simulation experiments (OSSEs). Generation of the GEOS-5 Nature Run (G5NR) was motivated in part by the desire of the OSSE community for an improved high-resolution sequel to an existing Nature Run produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), which has served the community for several years. The intended use of the G5NR in this context is for generating simulated observations to test proposed observing system designs regarding new instruments and their deployments. Because NASA's interest in OSSEs extends beyond traditional weather forecasting applications, the G5NR includes, in addition to standard meteorological components, a suite of aerosol types and several trace gas concentrations, with emissions downscaled to 10 km using ancillary information such as power plant location, population density and night-light information. The evaluation exercise described here involved more than twenty-five GMAO scientists investigating various aspects of the G5NR performance, including time mean temperature and wind fields, energy spectra, precipitation and the hydrological cycle, the representation of waves, tropical cyclones and midlatitude storms, land and ocean surface characteristics, the representation and forcing effects of clouds and radiation, dynamics of the stratosphere and mesosphere, and the representation of aerosols and trace gases. Comparisons are made with observational data sets when possible, as well as with reanalyses and other long model simulations. The evaluation is broad in scope, as it is meant to assess the overall realism of basic aspects of the G5NR deemed relevant to the conduct of OSSEs. However, because of the relatively short record and other practical considerations, these comparisons cannot provide a definitive, statistically sound assessment of all model deficiencies, or guarantee the G5NR's suitability for all OSSE applications. Differences between the observed and simulated behavior also must be judged in the context of basic internal atmospheric variability which can introduce variations that are not necessarily controlled by the prescribed sea surface temperatures used in generating the G5NR. The results show that the G5NR performs well as measured by the majority of metrics applied in this evaluation. Particular benefits derived from the 7-km resolution of G5NR include realistic representations of extreme weather events in both the tropics and extratropics including tropical cyclones, Nor'easters and mesoscale convective complexes; improved representation of the diurnal cycle of precipitation over land; well-resolved surface-atmosphere interactions such as katabatic wind flows over Antarctica and Greenland; and resolution of orographically generated gravity waves that propagate into the upper atmosphere and influence the large scale circulation. Obvious deficiencies in the G5NR include a "splitting" of the inter-tropical convergence zone, which leads to a weaker-than-observed Hadley circulation and related deficiencies in the depiction of stationary wave patterns. Also, while the G5NR captures global cloud features and radiative effects well in general, close comparison with observations reveals higher-than-observed cloud brightness, likely due to an overabundance of cloud condensate; less distinct cloud minima in subtropical subsidence zones, consistent with a weak Hadley circualtion; and too few near-coastal marine stratocumulus clouds.

Meteorological reanalyses provide multi-year gridded datasets that describe the evolution of the atmosphere. Such products use a data assimilation system, comprising of an atmospheric model, a broad suite of observations, and an analysis system that optimally combines the model forecast with the observations, using an algorithm that includes information about model and data accuracy. The mixture of observations is of central importance to the quality of the assimilated datasets. The Modern-era Retrospective Analysis for Research and Applications (MERRA) included constraints on the thermal structure of the middle atmosphere from nadir sounders on the NOAA polar-orbiting platforms (Stratospheric Sounding Units and Advanced Microwave Sounding Units). These instruments have peak sensitivities that occur well below the stratopause. As such, the radiance measurements do not provide strong constraints on stratopause temperature. The new MERRA-2 reanalysis is using EOS-MLS temperature retrievals after they are available: it will be demonstrated that these data lead to a more realistic stratopause structure in MERRA-2 than in MERRA. Similarly, the work demonstrates the improvements in lower stratospheric ozone in MERRA-2 than in MERRA, for the period when EOS-MLS ozone data are assimilated. This improvement occurs because of the ozone profile information offered by MLS in the low stratosphere, in contrast to the SBUV/2 data used for the rest of MERRA-2. The impacts of choosing to use the EOS-MLS datasets are discussed in context of the continuity of the data record in MERRA- 2.

We present a new method to diagnose the middle atmosphere climate sensitivity by extending the Climate Feedback-Response Analysis Method (CFRAM) for the coupled atmosphere-surface system to the middle atmosphere. The Middle atmosphere CFRAM (MCFRAM) is built on the atmospheric energy equation per unit mass with radiative heating and cooling rates as its major thermal energy sources. MCFRAM preserves the CFRAM unique feature of an additive property for which the sum of all partial temperature changes due to variations in external forcing and feedback processes equals the observed temperature change. In addition, MCFRAM establishes a physical relationship of radiative damping between the energy perturbations associated with various feedback processes and temperature perturbations associated with thermal responses. MCFRAM is applied to both measurements and model output fields to diagnose the middle atmosphere climate sensitivity. It is found that the largest component of the middle atmosphere temperature response to the 11-year solar cycle (solar maximum vs. solar minimum) is directly from the partial temperature change due to the variation of the input solar flux. Increasing CO2 always cools the middle atmosphere with time whereas partial temperature change due to O3 variation could be either positive or negative. The partial temperature changes due to different feedbacks show distinctly different spatial patterns. The thermally driven globally averaged partial temperature change due to all radiative processes is approximately equal to the observed temperature change, ranging from 0.5 K near 70 km from the near solar maximum to the solar minimum.

We examine the major stratosphere sudden warming (SSW) that occurred on 6 January 2013, using output from the NASA Global Modeling and Assimilation Office (GMAO) GEOS-5 (Goddard Earth Observing System) near-real-time data assimilation system (DAS). Results show that the major SSW of January 2013 falls into the vortex splitting type of SSW, with the initial planetary wave breaking occurring near 10 hPa. The vertical flux of wave activity at the tropopause responsible for the SSW occurred mainly in the Pacific Hemisphere, including the a pulse associated with the preconditioning of the polar vortex by wave 1 identified on 23 December 2012. While most of the vertical wave activity flux was in the Pacific Hemisphere, a rapidly developing tropospheric weather system over the North Atlantic on 28 December is shown to have produced a strong transient upward wave activity flux into the lower stratosphere coinciding with the peak of the SSW event. In addition, the GEOS-5 5-day forecasts accurately predicted the major SSW of January 2013 as well as the upper tropospheric disturbances responsible for the warming. The overall success of the 5-day forecasts provides motivation to produce regular 10-day forecasts with GEOS-5, to better support studies of stratosphere-troposphere interaction.

The Global Modeling Initiative (GMI) science team is developing a three dimensional chemistry and transport model (CTM) to be used in assessment of the atmospheric effects of aviation. Requirements are that this model be documented be validated against observations, use a realistic atmospheric circulation, and contain numerical transport and photochemical modules representing atmospheric processes. The model must also retain computational efficiency to be tractable to use for multiple scenarios and sensitivity studies. To meet these requirements, a facility model concept was developed in which the different components of the CTM are evaluated separately. The first use of the GMI model will be to evaluate the impact of the exhaust of supersonic aircraft on the stratosphere. The assessment calculations will depend strongly on the wind and temperature fields used by the CTM. Three meteorological data sets for the stratosphere are available to GMI: the National Center for Atmospheric Research Community Climate Model (CCM2), the Goddard Earth Observing System Data Assimilation System (GEOS DAS), and the Goddard Institute for Space Studies general circulation model (GISS). Objective criteria were established by the GMI team to identify the data set which provides the best representation of the stratosphere. Simulations of gases with simple chemical control were chosen to test various aspects of model transport. The three meteorological data sets were evaluated and graded based on their ability to simulate these aspects of stratospheric measurements. This paper describes the criteria used in grading the meteorological fields. The meteorological data set which has the highest score and therefore was selected for GMI is CCM2. This type of objective model evaluation establishes a physical basis for interpretation of differences between models and observations. Further, the method provides a quantitative basis for defining model errors, for discriminating between different models, and for ready re-evaluation of improved models. These in turn will lead to a higher level of confidence in assessment calculations.

Winds derived from a stratospheric and tropospheric data assimilation system (STRATAN) are compared with balance winds derived from National Meteorological Center/Climate Analysis Center (NMC/CAC) heights. At middle latitudes in the lower stratosphere, the results show that STRATAN winds are comparable to the balance winds. In addition STRATAN winds provide useful horizontal divergence analyses, and hence, vertical velocity fields. More generally, the STRATAN winds are useful in a more extended domain than the balanced winds. In particular, they are useful in the Tropics and the upper stratosphere where the balanced winds fail. The assimilation also captures the quasi-biennial oscillation, but does not do a good job of representing tropical waves.

Data from the mesosphere-stratosphere-troposphere radar at Jicamarca, Peru, together with other available data, are used to diagnose the mean structure of winds and gravity-wave momentum fluxes from the surface to 90 km during two-ten-day campaigns in June and August of 1987. In the stratosphere a layer of maximum eastward flow associated with the quasi-biennial oscillation (QBO) was seen to strengthen and descend rapidly from June to August, overlying persistent westward flow. A layer of enhanced signal return, suggestive of a turbulent layer, was observed just above the descending QBO eastward maximum. Notable zonal asymmetries were present during this transition and the local meridional circulation departed from zonal-mean QBO theory. A substantial northeastward momentum flux was found below 25 km, which may be related to topographic gravity waves excited by southeastward flow across the Andes. In the lower mesosphere a relatively weak 'second' mesopause semiannual oscillation is confirmed. In both the lower stratosphere and mesosphere, body forces were comparable in magnitude to inferred Coriolis torques.

Density fluctuations obtained along seven Space Shuttle reentry tracks are used to examine the horizontal structure and the vertical distribution of density variance in the mesosphere and lower thermosphere. The tracks lie primarily over open ocean at middle and low latitudes and represent the only measurements of horizontal atmospheric structure at these heights available to date. The density fluctuations are interpreted in terms of gravity-wave motions and reveal significant density (and velocity) variance at horizontal scales ranging from about 10 to 1000 km. Fluctuation amplitudes are used to infer corresponding velocity perturbations and characteristic vertical scales and frequencies of the wave spectrum. Results suggest that the mean velocity variance is smaller over the Pacific ocean than over major land masses, and that the variance increases with height in a manner consistent with that expected in the present of wave saturation processes.