Computations for the improvement of global atmospheric models
Global atmospheric models computations based on treating molecular and eddy diffusion as dynamic processes
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Global atmospheric models computations based on treating molecular and eddy diffusion as dynamic processes
The new generation of heterogeneous CPU/GPU computer systems offer much greater computational performance but are not yet widely used for climate modeling. One reason for this is that traditional climate models were written before GPUs were available and would require an extensive overhaul to run on these new machines. In addition, even conventional “high–resolution” simulations don't currently provide enough parallel work to keep GPUs busy, so the benefits of such overhaul would be limited for the types of simulations climate scientists are accustomed to. The vision of the Simple Cloud-Resolving Energy Exascale Earth System (E3SM) Atmosphere Model (SCREAM) project is to create a global atmospheric model with the architecture to efficiently use GPUs and horizontal resolution sufficient to fully take advantage of GPU parallelism. After 5 years of model development, SCREAM is finally ready for use. In this paper, we describe the design of this new code, its performance on both CPU and heterogeneous machines, and its ability to simulate real-world climate via a set of four 40 day simulations covering all 4 seasons of the year.
This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40-day prescribed-sea-surface-temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be ~ 6x the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long-standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.
When a new, better-formulated physical parameterization is introduced into a global atmospheric model, aspects of the global model solutions are sometimes degraded. Then, in order to use the new global model to address science questions, there is an incentive to restore its accuracy. Oftentimes this restoration is achieved by tuning of model parameter values. Unfortunately, the retuning process is expensive because characterizing the parameter dependence requires numerous time-consuming global simulations. To reduce the cost of tuning, this manuscript introduces a “poor man's” model tuner, “QuadTune”. QuadTune carves the globe into regions and approximates the model parameter dependence through the use of an uncorrelated quadratic emulator (i.e., response surface). The simplicity of the emulator reduces the required number of global model simulations and aids explainability of tuner behavior. Tuning removes parametric error but leaves behind model structural error. Structural error manifests itself as regional residual biases, such as stubborn biases and tuning trade-offs. To visualize these residual biases, QuadTune's software includes a set of diagnostic plots. This paper illustrates the use of the plots for characterizing residual biases with an example tuning problem.
The first global model of meteoric iron in the atmosphere (WACCM-Fe) has been developed by combining three components: the Whole Atmosphere Community Climate Model (WACCM), a description of the neutral and ion-molecule chemistry of iron in the mesosphere and lower thermosphere (MLT), and a treatment of the injection of meteoric constituents into the atmosphere. The iron chemistry treats seven neutral and four ionized iron containing species with 30 neutral and ion-molecule reactions. The meteoric input function (MIF), which describes the injection of Fe as a function of height, latitude, and day, is precalculated from an astronomical model coupled to a chemical meteoric ablation model (CABMOD). This newly developed WACCM-Fe model has been evaluated against a number of available ground-based lidar observations and performs well in simulating the mesospheric atomic Fe layer. The model reproduces the strong positive correlation of temperature and Fe density around the Fe layer peak and the large anticorrelation around 100 km. The diurnal tide has a significant effect in the middle of the layer, and the model also captures well the observed seasonal variations. However, the model overestimates the peak Fe+ concentration compared with the limited rocket-borne mass spectrometer data available, although good agreement on the ion layer underside can be obtained by adjusting the rate coefficients for dissociative recombination of Fe-molecular ions with electrons. Sensitivity experiments with the same chemistry in a 1-D model are used to highlight significant remaining uncertainties in reaction rate coefficients, and to explore the dependence of the total Fe abundance on the MIF and rate of vertical transport.
Abstract Can the current successes of global machine learning‐based weather simulators be generalized beyond 2‐week forecasts to stable and accurate multiyear runs? The recently developed AI2 Climate Emulator (ACE) suggests this is feasible, based upon 10‐year simulations with a network trained on output from a physics‐based global atmosphere model using a grid spacing of approximately 110 km and forced by a repeating annual cycle of sea‐surface temperature. Here we show that ACE, without modification, can be trained to emulate another major atmospheric model, EAMv2, run at a comparable grid spacing for at least 10 years with similarly small climate biases—a prerequisite to wider applicability. With an analysis that combines multiple temporal, spatial, and frequency domain perspectives, we show that ACE faithfully represents the spatiotemporal structure of EAMv2 precipitation and related variables. Finally, we show that a pretrained ACE network is able to adapt to a new global climate model simulation data set with 10 fewer training steps than when starting from random initialization, all while still maintaining low levels of climate bias. Further analysis of these fine‐tuning experiments reveal ACE's intriguing ability to interpolate between distinct global climate models.
The Global Reference Atmospheric Model (GRAM) program includes the capability for simulating pseudo-random perturbations in density, temperature, pressure, or wind components along a simulated reentry trajectory or other path through the atmosphere. Some concerns were expressed by GRAM users, however, that the mean-square perturbation gradients may be too large for small values of the vertical separation Delta z. The present GRAM perturbation simulations, based on a one-step autoregressive model, yield a power spectrum versus wavenumber k which is proportional to k sup -2 at high wavenumbers. This feature also produces mean-square perturbation differences which are directly proportional to Delta z, and mean-square perturbation gradients which are inversely proportional to Delta z. Thus, root-mean-square gradients, (Delta f/Delta z) sub rms, increase with decreasing Delta a as Delta z sup -1/2. A simple modification to GRAM is suggested which overcomes this problem, i.e., which produce root-mean-square gradient that remain bound as Delta z approaches zero. Possible applications of more sophisticated simulation approaches, e.g., second order autoregressive models, or fractal model techniques, were also explored briefly but found to yield improvements which appear too small to justify their considerable added complexity for use in the GRAM programs.
Venus Global Reference Atmospheric Model (Venus-GRAM) is an engineering-level atmospheric model developed by MSFC that is widely used for diverse mission applications including: Systems design; Performance analysis; Operations planning for aerobraking, Entry, Descent and Landing, and aerocapture; Is not a forecast model; Outputs include density, temperature, pressure, wind components, and chemical composition; Provides dispersions of thermodynamic parameters, winds, and density; Optional trajectory and auxiliary profile input files Has been used in multiple studies and proposals including NASA Engineering and Safety Center (NESC) Autonomous Aerobraking and various Discovery proposals; Released in 2005; Available at: https://software.nasa.gov/software/MFS-32314-1.
Surface spherical harmonics were used to analyze the horizontal fields of various quantities generated by a global climate model. Also, the computed monthly mean forecast fields were compared with the corresponding observed fields.
Progress in better determining CO2 sources and sinks will almost certainly rely on utilization of more extensive and intensive CO2 and related observations including those from satellite remote sensing. Use of advanced data requires improved modeling and analysis capability. Under NASA Carbon Cycle Science support we seek to develop and integrate improved formulations for 1) atmospheric transport, 2) terrestrial uptake and release, 3) biomass and 4) fossil fuel burning, and 5) observational data analysis including inverse calculations. The transport modeling is based on meteorological data assimilation analysis from the Goddard Modeling and Assimilation Office. Use of assimilated met data enables model comparison to CO2 and other observations across a wide range of scales of variability. In this presentation we focus on the short end of the temporal variability spectrum: hourly to synoptic to seasonal. Using CO2 fluxes at varying temporal resolution from the SIB 2 and CASA biosphere models, we examine the model's ability to simulate CO2 variability in comparison to observations at different times, locations, and altitudes. We find that the model can resolve much of the variability in the observations, although there are limits imposed by vertical resolution of boundary layer processes. The influence of key process representations is inferred. The high degree of fidelity in these simulations leads us to anticipate incorporation of realtime, highly resolved observations into a multiscale carbon cycle analysis system that will begin to bridge the gap between top-down and bottom-up flux estimation, which is a primary focus of NACP.
Machine learning (ML)-based models have demonstrated high skill and computational efficiency, often outperforming conventional physics-based models in weather and subseasonal predictions. While prior studies have assessed their fidelity in capturing synoptic-scale atmospheric dynamics, their performance across timescales and under out-of-distribution forcing, such as +3K or +4K uniform-warming forcings, and the sources of biases remain elusive, to establish the model's reliability for Earth science. Here, we design three sets of experiments targeting synoptic-scale phenomena, interannual variability, and out-of-distribution uniform-warming forcings. We evaluate the Neural General Circulation Model (NeuralGCM), a hybrid model integrating a dynamical core with ML-based component, against observations and physics-based Earth system models (ESMs). At the synoptic scale, NeuralGCM captures the evolution and propagation of extratropical cyclones with performance comparable to ESMs. At the interannual scale, when forced by El Niño-Southern Oscillation sea surface temperature (SST) anomalies, NeuralGCM successfully reproduces associated teleconnection patterns but exhibits deficiencies in capturing nonlinear response. Under out-of-distribution uniform-warming forcings, NeuralGCM simulates similar responses in global-average temperature and precipitation and reproduces large-scale tropospheric circulation features similar to those in ESMs. Notable weaknesses include overestimating the tracks and spatial extent of extratropical cyclones, biases in the teleconnected wave train triggered by tropical SST anomalies, and differences in upper-level warming and stratospheric circulation responses to SST warming compared to physics-based ESMs. The causes of these weaknesses were explored. Despite the noted weaknesses, NeuralGCM reproduces responses across experiments reasonably and performs comparably to ESMs. By integrating a dynamical core with ML, NeuralGCM shows potential for developing ML-based ESMs.
Abstract Clouds are parameterized in climate models using quantities on the model grid‐scale to approximate the cloud cover and impact on radiation. Because of the complexity of processes involved with clouds, these parameterizations are one of the key challenges in climate modeling. Differences in parameterizations of clouds are among the main contributors to the spread in climate sensitivity across models. In this work, the clouds in three generations of an atmosphere model lineage are evaluated against satellite observations. Satellite simulators are used within the model to provide an appropriate comparison with individual satellite products. In some respects, especially the top‐of‐atmosphere cloud radiative effect, the models show generational improvements. The most recent generation, represented by two distinct branches of development, exhibits some regional regressions in the cloud representation; in particular the southern ocean shows a positive bias in cloud cover. The two branches of model development show how choices during model development, both structural and parametric, lead to different cloud climatologies. Several evaluation strategies are used to quantify the spatial errors in terms of the large‐scale circulation and the cloud structure. The Earth mover's distance is proposed as a useful error metric for the passive satellite data products that provide cloud‐top pressure‐optical depth histograms. The cloud errors identified here may contribute to the high climate sensitivity in the Community Earth System Model, version 2 and in the Energy Exascale Earth System Model, version 1.
An empirical 4-D atmospheric model is described, whose ultimate goal is to give pressure, temperature, and density variables and their structure as a function of latitude, longitude, altitude, and the time domain (seasonal and perhaps time of day) over the altitude range from sea level to 185 km. The computer program combines the 4-D model of Siegler and Fowler (1972) for the range 0-25 km altitude, the modified Groves (1971) model for the range 30-90 km, and the Jacchia (1970) model for the range 115-185 km. Between 90 and 115 km, a smooth transition between the Groves values and Jacchia values was obtained by a fairing technique, while between 25 and 30 km an interpolation scheme was used between the 4-D results and the Groves values. The program is designed to give atmospheric parameter values along a simulated trajectory and a profile, such as a vertical profile at a single location. A set of trajectory positions for any desired trajectory, ascent or re-entry, can be input to the program for evaluation of atmospheric parameters for design purposes.
Substitution of observed monthly mean sea-surface temperatures (SSTs) as lower boundary conditions, in place of climatological SSTs, failed to improve the model simulations. While the impact of SST anomalies on the model output is greater at sea level than at upper levels the impact on the monthly mean simulations is not beneficial at any level. Shifts of one and two days in initialization time produced small, but non-trivial, changes in the model-generated monthly mean synoptic fields. No improvements in the mean simulations resulted from the use of either time-averaged initial data or re-initialization with time-averaged early model output. The noise level of the model, as determined from a multiple initial state perturbation experiment, was found to be generally low, but with a noisier response to initial state errors in high latitudes than the tropics.
Global aerosol models are used to study the distribution and properties of atmospheric aerosol particles as well as their effects on clouds, atmospheric chemistry, radiation, and climate. The present article provides an overview of the basic concepts of global atmospheric aerosol modeling and shows some examples from a global aerosol simulation. Particular emphasis is placed on the simulation of aerosol particles and their effects within global climate models.
The Venus Global Reference Atmospheric Model (Venus-GRAM) is an engineering-level atmospheric model applicable for engineering design analyses, mission planning, and operational decision making. Missions to Venus have generated a wealth of atmospheric data, however, Venus-GRAM has not been updated since its development and release in 2005. GRAM upgrades and maintenance have depended on inconsistent and waning project-specific support. The NASA Science Mission Directorate (SMD) has agreed to provide funding support in Fiscal Year 2018 and 2019 to upgrade the GRAMs. This presentation will provide an overview of Venus-GRAM and the objectives, tasks, and milestones related to the GRAM upgrades.
The status of the Global Reference Atmospheric Model (GRAM) and the Mars Global Reference Atmospheric Model (MARS-GRAM) is reviewed. The wavelike perturbations observed in the Viking 1 and 2 surface pressure data, in the Mariner 9 IR spectroscopy data, and in the Viking 1 and 2 lander entry profiles were studied and the results interpreted.
Global Reference Atmospheric Model (GRAM-99) is an engineering-level model of the Earth's atmosphere. It provides both mean values and perturbations for density, temperature, pressure, and winds, as well as monthly- and geographically-varying trace constituent concentrations. From 0-27 km, thermodynamics and winds are based on National Oceanic and Atmospheric Administration Global Upper Air Climatic Atlas (GUACA) climatology. Above 120 km, GRAM is based on the NASA Marshall Engineering Thermosphere (MET) model. In the intervening altitude region, GRAM is based on Middle Atmosphere Program (MAP) climatology that also forms the basis of the 1986 COSPAR Intemationa1 Reference Atmosphere (CIRA). MAP data in GRAM are augmented by a specially-derived longitude variation climatology. Atmospheric composition is represented in GRAM by concentrations of both major and minor species. Above 120 km, MET provides concentration values for N2, O2, Ar, O, He, and H. Below 120 km, species represented also include H2O, O3, N2O, CO, CH, and CO2. Water vapor in GRAM is based on a combination of GUACA, Air Force Geophysics Laboratory (AFGL), and NASA Langley Research Center climatologies. Other constituents below 120 km are based on a combination of AFGL and h4AP/CIRA climatologies. This report presents results of comparisons between GRAM Constituent concentrations and those provided by the Naval Research Laboratory (NRL) climatology of Summers (NRL,/MR/7641-93-7416, 1993). GRAM and NRL concentrations were compared for seven species (CH4, CO, CO2, H2O, N2O, O2, and O3) for months January, April, July, and October, over height range 0-115 km, and latitudes -90deg to + 90deg at 10deg increments. Average GRAM-NRL correlations range from 0.878 (for CO) to 0.975 (for O3), with an average over all seven species of 0.936 (standard deviation 0.049).