Engineering topics
Molod, Andrea
Publications and source records attributed to Molod, Andrea.
Tropical Waves and the Quasi-Biennial Oscillation in a 7-km Global Climate Simulation
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.
Technical Report Series on Global Modeling and Data Assimilation, Volume 43: Initial Evaluation of the Climate - MERRA-2
The years since the introduction of MERRA have seen numerous advances in the GEOS-5 Data Assimilation System as well as a substantial decrease in the number of observations that can be assimilated into the MERRA system. To allow continued data processing into the future, and to take advantage of several important innovations that could improve system performance, a decision was made to produce MERRA-2, an updated retrospective analysis of the full modern satellite era. One of the many advances in MERRA-2 is a constraint on the global dry mass balance; this allows the global changes in water by the analysis increment to be near zero, thereby minimizing abrupt global interannual variations due to changes in the observing system. In addition, MERRA-2 includes the assimilation of interactive aerosols into the system, a feature of the Earth system absent from previous reanalyses. Also, in an effort to improve land surface hydrology, observations-corrected precipitation forcing is used instead of model-generated precipitation. Overall, MERRA-2 takes advantage of numerous updates to the global modeling and data assimilation system. In this document, we summarize an initial evaluation of the climate in MERRA-2, from the surface to the stratosphere and from the tropics to the poles. Strengths and weaknesses of the MERRA-2 climate are accordingly emphasized.
Evaluation of the 7-km GEOS-5 Nature Run
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.
Using Ensemble Short-Term Initialized Coupled NASA GEOS5 Climate Model Integrations to Study Convective Bias Growth
The representation of convective processes, particularly deep convection in the tropics, remains a persistent problem in climate models. In fact structural biases in the distribution of tropical rainfall in the CMIP5 models is hardly different than that of the CMIP3 versions. Given that regional climate change at higher latitudes is sensitive to the configuration of tropical forcing, this persistent bias is a major issue for the credibility of climate change projections. In this study we use model output from integrations of the NASA Global Earth Observing System Five (GEOS5) climate modeling system to study the evolution of biases in the location and intensity of convective processes. We take advantage of a series of hindcast experiments done in support of the US North American Multi-Model Ensemble (NMME) initiative. For these experiments a nine-month forecast using a coupled model configuration is made approximately every five days over the past 30 years. Each forecast is started with an updated analysis of the ocean, atmosphere and land states. For a given calendar month we have approximately 180 forecasts with daily means of various quantities. These forecasts can be averaged to essentially remove "weather scales" and highlight systematic errors as they evolve. Our primary question is to ask how the spatial structure of daily mean precipitation over the tropics evolves from the initial state and what physical processes are involved. Errors in parameterized convection, various water and energy fluxes and the divergent circulation are found to set up on fast time scales (order five days) compared to errors in the ocean, although SST changes can be non-negligible over that time. For the month of June the difference between forecast day five versus day zero precipitation looks quite similar to the difference between the June precipitation climatology and that from the Global Precipitation Climatology Project (GPCP). We focus much of our analysis on the influence of SST gradients, associated PBL baroclinicity enabled by turbulent mixing, the ensuing PBL moisture convergence, and how changes in these processes relate to convective precipitation bias growth over this short period.
Development of the GEOS-5 Atmospheric General Circulation Model: Evolution from MERRA to MERRA2.
The Modern-Era Retrospective Analysis for Research and Applications-2 (MERRA2) version of the GEOS-5 (Goddard Earth Observing System Model - 5) Atmospheric General Circulation Model (AGCM) is currently in use in the NASA Global Modeling and Assimilation Office (GMAO) at a wide range of resolutions for a variety of applications. Details of the changes in parameterizations subsequent to the version in the original MERRA reanalysis are presented here. Results of a series of atmosphere-only sensitivity studies are shown to demonstrate changes in simulated climate associated with specific changes in physical parameterizations, and the impact of the newly implemented resolution-aware behavior on simulations at different resolutions is demonstrated. The GEOS-5 AGCM presented here is the model used as part of the GMAO's MERRA2 reanalysis, the global mesoscale "nature run", the real-time numerical weather prediction system, and for atmosphere-only, coupled ocean-atmosphere and coupled atmosphere-chemistry simulations. The seasonal mean climate of the MERRA2 version of the GEOS-5 AGCM represents a substantial improvement over the simulated climate of the MERRA version at all resolutions and for all applications. Fundamental improvements in simulated climate are associated with the increased re-evaporation of frozen precipitation and cloud condensate, resulting in a wetter atmosphere. Improvements in simulated climate are also shown to be attributable to changes in the background gravity wave drag, and to upgrades in the relationship between the ocean surface stress and the ocean roughness. The series of "resolution aware" parameters related to the moist physics were shown to result in improvements at higher resolutions, and result in AGCM simulations that exhibit seamless behavior across different resolutions and applications.
Constraints on the Profiles of Total Water PDF in AGCMs from AIRS and a High-Resolution Model
Atmospheric general circulation model (AGCM) cloud parameterizations generally include an assumption about the subgrid-scale probability distribution function (PDF) of total water and its vertical profile. In the present study, the Atmospheric Infrared Sounder (AIRS) monthly-mean cloud amount and relative humidity fields are used to compute a proxy for the second moment of an AGCM total water PDF called the RH01 diagnostic, which is the AIRS mean relative humidity for cloud fractions of 0.1 or less. The dependence of the second moment on horizontal grid resolution is analyzed using results from a high-resolution global model simulation.The AIRS-derived RH01 diagnostic is generally larger near the surface than aloft, indicating a narrower PDF near the surface, and varies with the type of underlying surface. High-resolution model results show that the vertical structure of profiles of the AGCM PDF second moment is unchanged as the grid resolution changes from 200 to 100 to 50 km, and that the second-moment profiles shift toward higher values with decreasing grid spacing.Several Goddard Earth Observing System, version 5 (GEOS-5), AGCM simulations were performed with several choices for the profile of the PDF second moment. The resulting cloud and relative humidity fields were shown to be quite sensitive to the prescribed profile, and the use of a profile based on the AIRS-derived proxy results in improvements relative to observational estimates. The AIRS-guided total water PDF profiles, including their dependence on underlying surface type and on horizontal resolution, have been implemented in the version of the GEOS-5 AGCM used for publicly released simulations.
Report of the 4th World Climate Research Programme International Conference on Reanalyses
The 4th WCRP International Conference on Reanalyses provided an opportunity for the international community to review and discuss the observational and modelling research, as well as process studies and uncertainties associated with reanalysis of the Earth System and its components. Characterizing the uncertainty and quality of reanalyses is a task that reaches far beyond the international community of producers, and into the interdisciplinary research community, especially those using reanalysis products in their research and applications. Reanalyses have progressed greatly even in the last 5 years, and newer ideas, projects and data are coming forward. While reanalysis has typically been carried out for the individual domains of atmosphere, ocean and land, it is now moving towards coupling using Earth system models. Observations are being reprocessed and they are providing improved quality for use in reanalysis. New applications are being investigated, and the need for climate reanalyses is as strong as ever. At the heart of it all, new investigators are exploring the possibilities for reanalysis, and developing new ideas in research and applications. Given the many centres creating reanalyses products (e.g. ocean, land and cryosphere research centres as well as NWP and atmospheric centers), and the development of new ideas (e.g. families of reanalyses), the total number of reanalyses is increasing greatly, with new and innovative diagnostics and output data. The need for reanalysis data is growing steadily, and likewise, the need for open discussion and comment on the data. The 4th Conference was convened to provide a forum for constructive discussion on the objectives, strengths and weaknesses of reanalyses, indicating potential development paths for the future.
The GEOS-5 Atmospheric General Circulation Model: Mean Climate and Development from MERRA to Fortuna
This report is a documentation of the Fortuna version of the GEOS-5 Atmospheric General Circulation Model (AGCM). The GEOS-5 AGCM is currently in use in the NASA Goddard Modeling and Assimilation Office (GMAO) for simulations at a wide range of resolutions, in atmosphere only, coupled ocean-atmosphere, and data assimilation modes. The focus here is on the development subsequent to the version that was used as part of NASA s Modern-Era Retrospective Analysis for Research and Applications (MERRA). We present here the results of a series of 30-year atmosphere-only simulations at different resolutions, with focus on the behavior of the 1-degree resolution simulation. The details of the changes in parameterizations subsequent to the MERRA model version are outlined, and results of a series of 30-year, atmosphere-only climate simulations at 2-degree resolution are shown to demonstrate changes in simulated climate associated with specific changes in parameterizations. The GEOS-5 AGCM presented here is the model used for the GMAO s atmosphere-only and coupled CMIP-5 simulations.
A Comparison of FIFE Observation with GEOS Assimilated Data Including a Heterogeneous LSM
Several recent studies have shown that much can be learned by comparing grid-point data from a data assimilation system with in-situ observations from field experiments. While the surface heterogeneity is acknowledged in these studies, they lack quantitative representations of the influence of heterogeneity on the near-surface meteorology and surface hydrologic and energy balance. Here, we use the Betts and Ball FIFE site-averaged data. Standard deviations of the site-average will provide an estimate of the FIFE site heterogeneity. Recently, the Mosaic Land-Surface Model (LSM) has been incorporated into the Goddard Earth Observing System (GEOS) Data Assimilation System (DAS). The Mosaic LSM computes the surface energy and hydrologic balance for nine distinct surface types at each grid-point. Each surface type is proportionally weighted to determine the mean grid point properties. Hence, we can compare modeled and observed grid-point variability in addition to the mean properties. Also, assimilated data sets created with and without the LSM are compared. The results indicate the importance of including quantitative estimates of heterogeneity in the analysis of the land surface hydrology and energy balances in assimilation systems.
Validation of Land-Surface Mosaic Heterogeneity in the GEOS DAS
The Mosaic Land-surface Model (LSM) has been included into the current GEOS Data Assimilation System (DAS). The LSM uses a more advanced representation of physical processes than previous versions of the GEOS DAS, including the representation of sub-grid heterogeneity of the land-surface through the Mosaic approach. As a first approximation, Mosaic assumes that all similar surface types within a grid-cell can be lumped together as a single'tile'. Within one GCM grid-cell, there might be 1 - 5 different tiles or surface types. All tiles are subjected to the grid-scale forcing (radiation, air temperature and specific humidity, and precipitation), and the sub-grid variability is a function of the tile characteristics. In this paper, we validate the LSM sub-grid scale variability (tiles) using a variety of surface observing stations from the Southern Great Plains (SGP) site of the Atmospheric Radiation Measurement (ARM) Program. One of the primary goals of SGP ARM is to study the variability of atmospheric radiation within a G,CM grid-cell. Enough surface data has been collected by ARM to extend this goal to sub-grid variability of the land-surface energy and water budgets. The time period of this study is the Summer of 1998 (June I - September 1). The ARM site data consists of surface meteorology, energy flux (eddy correlation and bowen ratio), soil water observations spread over an area similar to the size of a G-CM grid-cell. Various ARM stations are described as wheat and alfalfa crops, pasture and range land. The LSM tiles considered at the grid-space (2 x 2.5) nearest the ARM site include, grassland, deciduous forests, bare soil and dwarf trees. Surface energy and water balances for each tile type are compared with observations. Furthermore, we will discuss the land-surface sub-grid variability of both the ARM observations and the DAS.
Validation of the Land-Surface Energy Budget and Planetary Boundary Layer for Several Intensive field Experiments
Land-surface processes in a data assimilation system influence the lower troposphere and must be properly represented. With the recent incorporation of the Mosaic Land-surface Model (LSM) into the GEOS Data Assimilation System (DAS), the detailed land-surface processes require strict validation. While global data sources can identify large-scale systematic biases at the monthly timescale, the diurnal cycle is difficult to validate. Moreover, global data sets rarely include variables such as evaporation, sensible heat and soil water. Intensive field experiments, on the other hand, can provide high temporal resolution energy budget and vertical profile data for sufficiently long periods, without global coverage. Here, we evaluate the GEOS DAS against several intensive field experiments. The field experiments are First ISLSCP Field Experiment (FIFE, Kansas, summer 1987), Cabauw (as used in PILPS, Netherlands, summer 1987), Atmospheric Radiation Measurement (ARM, Southern Great Plains, winter and summer 1998) and the Surface Heat Budget of the Arctic Ocean (SHEBA, Arctic ice sheet, winter and summer 1998). The sites provide complete surface energy budget data for periods of at least one year, and some periods of vertical profiles. This comparison provides a detailed validation of the Mosaic LSM within the GEOS DAS for a variety of climatologic and geographic conditions.
Precipitation Recycling in the NASA GEOS Data Assimilation System
Analysis of precipitation recycling can improve the understanding of regional hydrologic anomalies, especially their evolution and maintenance. Diagnostic models of the recycling of precipitation and are applied to 15 years of the NASA Goddard Earth Observing System (GEOS) Data Assimilation System (DAS). Recycled precipitation is defined as the fraction of precipitation within a given region that originated as surface evaporation from the same region. The focus of the present work is on the interannual variability of the central United States hydrologic cycle and precipitation recycling. The extreme years of 1988 (drought) and 1993 (flood) are compared with the 15 year base period mean annual cycle. The results indicate that recycling ratio (the amount of precipitation with a local source relative to the total precipitation) is greater in 1988 than both the base period mean and the 1993 season (with 1993 recycling ratio less than the mean). On the other hand, both the summers of 1988 and 1993 show less total recycled precipitation than the mean. The results also show that precipitation recycling may have been more important in the spring of 1993, when the region was primed for flooding, than the summer, when the sever flooding occurred. The diagnostic approaches to precipitation recycling suffer from some weaknesses. Numerical simulations and assimilation using passive tracers have the potential to provide more accurate calculations of precipitation recycling and the remote sources of water. This ability is being incorporated into the latest GEOS data assimilation system, and some preliminary results will be presented.
The Vertical Extent of the Influence of the Land Surface Boundary on the Atmosphere Above
Heterogeneities in the land surface type on scales smaller than the typical General Circulation Model (GMC) grid size have long presented a challenge to properly modelling the impact of the subgrid scale variability on the grid scale. In particular, at some height above the heterogeneous vegetated surface the atmospheric properties (temperature, humidity, for example) become homogenized. This is an issue with implications for climate modelling as well as for the optimal use of field data for climate studies. There is evidence that the influence of the soil and/or vegetation properties extends well beyond the atmospheric surface layer. The present study is aimed at understanding the implications of this influence for a GCM's description of the coupling at the earth-atmosphere interface.
The Impact of Very Fine Vertical Resolution on the Climate of the GEOS-2 GCM
Typical General Circulation Models (GCM's) that are used for tropospheric climate studies extend vertically with significant resolution only into the lower stratosphere, while GCM's specifically designed for the study of the stratospheric climate do not finely resolve the tropospheric jet or boundary layer regions. The NASA/Goddard Data Assimilation Office (DAO) Goddard Earth Observing System GCM, (GEOS-2 GCM) which is a component of the GEOS-2 Data Assimilation System (DAS), is configured with a sufficient number of vertical levels to provide very fine resolution in the tropospheric and stratospheric regimes. The present study is aimed at understanding the benefits for accurate simulation of the climate seen in the GEOS-2 GCM due to its choice of vertical resolution. Raising the model lid to the lower mesosphere so as to explicitly resolve the stratospheric circulation, in particular the descent over the winter pole, is included in GEOS-2 to address the common GCM deficiency of a stratospheric cold pole bias. Increasing the number of levels is directed at improving the representation of the planetary boundary layer, the troposphere at the altitude of the sub-tropical jet stream, and the stratosphere. Comparisons of the climate of simulations with three different vertical resolution formulations of the GEOS-2 GCM are presented. In addition to the reduction in the climate bias, the impact of the resolution on model variability will be briefly discussed.
Technical report series on global modeling and data assimilation. Volume 6: A multiyear assimilation with the GEOS-1 system: Overview and results
The Data Assimilation Office (DAO) at Goddard Space Flight Center has produced a multiyear global assimilated data set with version 1 of the Goddard Earth Observing System Data Assimilation System (GEOS-1 DAS). One of the main goals of this project, in addition to benchmarking the GEOS-1 system, was to produce a research quality data set suitable for the study of short-term climate variability. The output, which is global and gridded, includes all prognostic fields and a large number of diagnostic quantities such as precipitation, latent heating, and surface fluxes. Output is provided four times daily with selected quantities available eight times per day. Information about the observations input to the GEOS-1 DAS is provided in terms of maps of spatial coverage, bar graphs of data counts, and tables of all time periods with significant data gaps. The purpose of this document is to serve as a users' guide to NASA's first multiyear assimilated data set and to provide an early look at the quality of the output. Documentation is provided on all the data archives, including sample read programs and methods of data access. Extensive comparisons are made with the corresponding operational European Center for Medium-Range Weather Forecasts analyses, as well as various in situ and satellite observations. This document is also intended to alert users of the data about potential limitations of assimilated data, in general, and the GEOS-1 data, in particular. Results are presented for the period March 1985-February 1990.
Technical report series on global modeling and data assimilation. Volume 1: Documentation of the Goddard Earth Observing System (GEOS) General Circulation Model, version 1
This technical report documents Version 1 of the Goddard Earth Observing System (GEOS) General Circulation Model (GCM). The GEOS-1 GCM is being used by NASA's Data Assimilation Office (DAO) to produce multiyear data sets for climate research. This report provides a documentation of the model components used in the GEOS-1 GCM, a complete description of model diagnostics available, and a User's Guide to facilitate GEOS-1 GCM experiments.