This study describes the modifications made to the Goddard Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) to conserve atmospheric dry‐air mass and to guarantee that the net ...source of water from precipitation and surface evaporation equals the change in total atmospheric water. The modifications involve changes to both the atmospheric model and the analysis procedure. In the model, sources and sinks of water are included in the continuity equation; in the analysis, constraints are imposed to penalize (and thus minimize) analysis increments of dry‐air mass. Finally, changes are also required to the Incremental Analysis Update (IAU) procedure. The effects of these modifications are evaluated separately in free‐running and assimilation experiments. Results are also presented from a multiyear reanalysis (Version 2 of the Modern Era Retrospective‐Analysis for Research and Applications: MERRA‐2) that uses the modified system.
Many studies have shown that recent reanalyses (e.g., NASA's Global Modeling and Assimilation Office (GMAO) Modern Era Retrospective‐Analysis for Research and Applications: MERRA, and European Centre for Medium‐Range Weather Forecasts: ERA‐Interim) contain spurious trends in globally integrated dry‐air mass, as well as large variations in globally integrated total precipitation associated with changes in the satellite observing systems. In this study we examine modifications to GMAO's Goddard Earth Observing System‐Version 5 (GEOS‐5) Data Assimilation System (DAS) which preserves total dry‐air mass and maintains the natural global balance between total precipitation and surface evaporation. Results are shown from targeted historical periods, as well as from GMAO's 35‐year MERRA‐2 Reanalysis.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
A recent attempt to downscale the 50 km MERRA-2 analyses to 7 km revealed an instability associated with the Incremental Analysis Update (IAU) procedure that has thus far gone unnoticed. A ...theoretical study based on a simple damped harmonic oscillator with complex frequency provides the framework to diagnose the problem and suggests means to avoid it. Three possible approaches to avoid the instability are to: (i) choose an "ideal" ratio of the lengths of the Predictor and Corrector steps of IAU based on a theoretical stability diagram; (ii) time average the background fields used to construct the IAU tendencies with given frequency; or (iii) apply a digital filter modulation to the IAU tendencies. All these are shown to control the instability for a wide range of resolutions when doing up- or down-scaling, experiments with the NASA/GMAO atmospheric general circulation model. Furthermore, it is found that combining IAU with the ensemble re-centering step typical of hybrid ensemble-variational approaches, also results in an instability based on the same mechanisms in the members of the ensemble. An example of such occurrence arises in an experiment performed with the GMAO 12.8 km hybrid 4D-EnVar system. Modulation of the ensemble IAU tendencies with a digital filter is shown to avoid the instability. In addition, the stability of certain 4DIAU implementations is analyzed and a suggestion is made to improve its results, though a complete study of this subject is postponed to a follow up work.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
Offline chemical transport models (CTMs) have traditionally been used to perform studies of atmospheric chemistry in a fixed dynamical environment. An alternative to using CTMs is to constrain the ...flow in a general circulation model using winds from meteorological analyses. The Goddard Earth Observing System (GEOS) “replay” approach involves reading in analyzed fields every 6 h and recomputing the analysis increments, which are applied as a forcing to the meteorology at every model time step. Unlike in CTM, all of the subgrid‐scale processes are recalculated online so that they are consistent with the large‐scale analysis fields, similar in spirit to “nudged” simulations, in which the online meteorology is relaxed to the analysis. Here we compare the transport of idealized tracers in different replay simulations constrained with meteorological fields taken from The Modern‐Era Retrospective Analysis for Research and Applications, Version 2 (MERRA‐2). We show that there are substantial differences in their large‐scale stratospheric transport, depending on whether analysis fields or assimilated fields are used. Replay simulations constrained with the instantaneous analysis fields produce stratospheric mean age values that are up to 30% too young relative to observations; by comparison, simulations constrained with the time‐averaged assimilated fields produce more credible stratospheric transport. Our study indicates that care should be taken to correctly configure the model when the replay technique is used to simulate stratospheric composition.
Key Points
GEOS replay simulations produce credible large‐scale stratospheric and tropospheric transport for use in studies of atmospheric composition
Simulations constrained with analysis fields produce stratospheric mean ages that are too young, versus when assimilated fields are used
By comparison, large‐scale tropospheric transport properties are relatively insensitive to whether analysis or assimilated fields are used
Full text
Available for:
DOBA, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, SIK, UILJ, UKNU, UL, UM, UPUK
The age of air has recently emerged as a diagnostic of atmospheric transport unaffected by chemical parameterizations, and the features in the age distributions computed in models have been ...interpreted in terms of the models' large-scale circulation field.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The development of and results obtained with a variable-resolution stretched-grid GCM for the regional climate simulation mode are presented. A global variable-resolution stretched grid used in the ...study has enhanced horizontal resolution over the United States as the area of interest. The stretched-grid approach is an ideal tool for representing regional- to global-scale interactions. It is an alternative to the widely used nested-grid approach introduced over a decade ago as a pioneering step in regional climate modeling. The major results of the study are presented for the successful stretched-grid GCM simulation of the anomalous climate event of the 1988 U.S. summer drought. The straightforward (with no updates) 2-month simulation is performed with 60-km regional resolution. The major drought fields, patterns, and characteristics, such as the time-averaged 500-hPa heights, precipitation, and the low-level jet over the drought area, appear to be close to the verifying analyses for the stretched-grid simulation. In other words, the stretched-grid GCM provides an efficient downscaling over the area of interest with enhanced horizontal resolution, in spite of degradation of skill over the coarser resolution far away from the area of interest. It is also shown that the stretched-grid GCM skill is sustained over the area of interest throughout the simulation extended to 1 yr. The stretched-grid GCM, developed and tested in a simulation mode, is a viable tool for regional and subregional climate studies and applications.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
Abstract
Multiyear (1987–97) limited ensemble integrations using a stretched-grid GCM, previously developed and experimented with by the authors, are employed for U.S. regional climate simulations. ...The ensemble members (six in total) are produced at two different regional resolutions: three members with 60-km and the other three members with 10-km regional resolution. The use of these two finer and coarser regional resolution ensemble members allows one to examine the impact of resolution on the overall quality of the simulated regional fields. For the multiyear ensemble simulations, an efficient regional downscaling to realistic mesoscales has been obtained. The ensemble means of the midtroposphere prognostic variables (height and meridional wind) show an overall good resemblance to the global reanalysis, especially for summer. Low-level features like the warm season Great Plains low-level jet are well represented in the simulations. During winter the 100-km simulations develop a southward wind east of the Rockies that is present neither in the reanalyses nor in the 60-km simulations. The analysis of the annual mean precipitation and its variance reveals that the ensemble simulations reproduce many of the observed features of a high-resolution rain gauge dataset analyzed on a 0.5° × 0.5° grid. Signal-to-noise ratios are larger than 1.5 s over a major part of the United States, especially over the Midwest and also over the mountainous regions like the Rockies and the Appalachians, suggesting that the orographic forcing is contributing to a larger signal. The ratios are smaller toward the eastern and western U.S. coastlines. This result could be attributed, at least in part, to limits in the representation of the land–sea contrasts.
For comparison purposes, an additional simulation has been performed using a global uniform 2° × 2.5° grid with the same number of global grid points as those of the above stretched grids. The stretched-grid GCM ensemble means show, overall, a better regional depiction of features than those of the uniform-grid GCM.
The results of the study show that even using limited ensemble integrations with a state-of-the-art stretched-grid GCM is beneficial for reducing the uncertainty of the multiyear regional climate simulation, especially when using finer 60-km regional resolution.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
In other words, the longer-term integrations suffer from an inaccurate interaction at the boundaries of the area of interest. ...a nested-grid approach does not allow one to perform continuous ...long-term integrations without updating initial and boundary conditions. ...the fine-resolution fields may provide only a marginal to moderate impact on the coarse-resolution fields within a buffer zone around the area of interest.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
The new stretched‐grid (SG) design with multiple (four) areas of interest, one at each global quadrant, is implemented into both a SG general circulation model (GCM) and a SG data assimilation system ...(DAS). The four areas of interest include the United States/northern Mexico, the El Niño area/central South America, India/China, and the eastern Indian Ocean/Australia. Both SG‐GCM and SG‐DAS annual (November 1997 to December 1998) integrations are performed with 50‐km regional resolution. The efficient regional downscaling to mesoscales is obtained for each of the four areas of interest, while the consistent interactions between regional and global scales and the high quality of global circulation are preserved. This is the advantage of the SG approach. The global variable‐resolution DAS incorporating the SG‐GCM has been developed and tested as an efficient tool for producing regional analyses and diagnostics with enhanced mesoscale resolution. The anomalous regional climate events of 1998 that occurred over the United States, Mexico, South America, China, India, African Sahel, and Australia are investigated in both simulation and data assimilation modes. The assimilated products are also used, along with gauge precipitation data, for validating the simulation results. The obtained results show that the SG‐GCM and SG‐DAS are capable of producing realistic high‐quality simulated and assimilated products at mesoscale resolution for regional climate studies and applications.
Derived Meteorological Products (DMPs, including potential temperature, potential vorticity (PV), equivalent latitude (EqL), horizontal winds and tropopause locations) from several meteorological ...analyses have been produced for the locations and times of measurements taken by several solar occultation instruments and the Aura Microwave Limb Sounder (MLS). MLS and solar occultation data are analyzed using DMPs to illustrate sampling issues that may affect interpretation and comparison of data sets with diverse sampling patterns and to provide guidance regarding the kinds of studies that benefit most from analyzing satellite data in relation to meteorological conditions using the DMPs. Using EqL or PV as a vortex‐centered coordinate does not alleviate all sampling problems, including those in studies using “vortex averages” of solar occultation data and in analyses of localized features (such as polar stratospheric clouds) and other fields that do not correlate well with PV. Using DMPs to view measurements with respect to their air mass characteristics is particularly valuable in studies of transport of long‐lived trace gases, polar processing in the winter lower stratosphere, and distributions and transport of O3 and other trace gases from the upper troposphere through the lower stratosphere. The comparisons shown here demonstrate good agreement between MLS and solar occultation data for O3, N2O, H2O, HNO3, and HCl; small biases are attributable to sampling effects or are consistent with detailed validation results presented elsewhere in this special section. The DMPs are valuable for many scientific studies and to facilitate validation of noncoincident measurements.