SEARCH · Engineering Papers
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.
In this talk, five specific major GCE improvements: (1) ice microphysics, (2) longwave and shortwave radiative transfer processes, (3) land surface processes, (4) ocean surface fluxes and (5) ocean mixed layer processes are presented. The performance of these new GCE improvements will be examined. Observations are used for model validation.
Cloud-resolving case-study simulations of convective transport and lightning NO production have yielded results which are directly applicable to the design of lightning parameterizations for global chemical transport models. In this work we have used cloud-resolving models (the Goddard Cumulus Ensemble Model (GCE) and MMS) to drive an off-line cloud-scale chemical transport model (CSCTM). The CSCTM, in conjunction with aircraft measurements of NO(x) in thunderstorms and ground-based lightning observations, has been used to constrain the amount of NO produced per flash. Cloud and chemistry simulations for several case studies of storms in different environments will be presented. Observed lightning flash rates have been incorporated into the CSCTM, and several scenarios of NO production per intracloud (IC) and per cloud-to-ground (CG) flash have been tested for each storm. The resulting NO(x) mixing ratios are compared with aircraft measurements taken within the storm (typically the anvil region) to determine the most likely NO production scenario. The range of values of NO production per flash (or per meter of lightning channel length) that have been deduced from the model will be shown and compared with values of production in the literature that have been deduced from observed NO spikes and from anvil flux calculations. Results show that on a per flash basis, IC flashes are nearly as productive of NO as CG flashes. This result simplifies the lightning parameterization for global models (ie., an algorithm for estimating the ICKG ratio is not necessary). Vertical profiles of lightning NOx mass at the end of the 3-D storm simulations have been summarized to yield suggested profiles for use in global models. Estimates of mean NO production per flash vary by a factor of three from one simulated storm to another. When combined with the global flash rate of 44 flashes per second from NASA's Optical Transient Detector (OTD) measurements, these estimates and the results from other techniques yield global NO production rates of 2-9 TgN/year. Simulations of the photochemistry over the 24 hours following a storm has been performed to determine the additional ozone production which can be attributed to lightning NO. Convective transport of HO(x), precursors leads to the generation of a HO(x), plume which substantially aids the downstream ozone production.
Parameterization of lightning NO(x) in global chemical transport models requires a method to specify the temporal and geographic distribution of flash rates, an estimate of NO production per flash, and a method to specify the effective vertical distribution of NO production. Flash rate parameterizations using the upward convective cloud mass flux have been developed and tested for several GCM and assimilated meteorological data sets. This parameterization and others involving the cloud-top height and convective precipitation provide a generally acceptable climatological distribution of flashes after normalization with satellite data. However, considerable improvement in terms of geographic and temporal variability may be achieved in the future using cloud microphysics data available from next-generation GCMs. Cloud-resolving case-study simulations of convective transport and lightning NO production have yielded results which are directly applicable to the NO production per flash and the vertical NO distribution portions of the global model lightning parameterizations. In this work we have used cloud-resolving models (the Goddard Cumulus Ensemble Model (GCE) and MM5) to drive an off-line cloud-scale chemical transport model (CSCTM). The CSCTM, in conjunction with aircraft measurements of NO(x) in thunderstorms and ground-based lightning observations, has been used to constrain the amount of NO produced per flash. Observed lightning flash rates have been incorporated into the CSCTM, and several scenarios of NO production per intracloud (IC) and per cloud-to-ground (CG) flash have been tested in each of several case study storms. The resulting NO(x) mixing ratios are compared with aircraft measurements taken within the storm (typically the anvil region) to determine the most likely NO production scenario. The range of values of NO production per flash (or per meter of lightning channel length) that have been deduced from the model will be shown and compared with values of production in the literature that have been deduced from observed NO spikes and from anvil flux calculations. Results show that on a per flash basis, IC flashes are nearly as productive of NO as CG flashes. Estimates of mean NO production per flash vary by a factor of three from one simulated storm to another. Vertical profiles of lightning NOx mass at the end of the 3-D storm simulations have been summarized to yield suggested profiles for use in global models. These profiles show a larger percentage of the total lightning NO(x) mass in the free troposphere than those suggested by Pickering et al. (1998).
Cloud-resolving case-study simulations of convective transport and lightning NO production have yielded results which are directly applicable to the design of lightning parameterizations for global chemical transport models. In this work we have used cloud-resolving models (the Goddard Cumulus Ensemble Model (GCE) and MMS) to drive an off-line cloud-scale chemical transport model (CSCTM). The CSCTM, in conjunction with aircraft measurements of NO x in thunderstorms and ground-l;>ased lightning observations, has been used to constrain the amount of NO produced per flash. Cloud and chemistry simulations for several case studies of storms in different environments will be presented. Observed lightning flash rates have been incorporated into the CSCTM, and several scenarios of NO production per intracloud (IC) and per cloud-to-ground (CG) flash have been tested for each storm. The resulting NOx mixing ratios are compared with aircraft measurements taken within the storm (typically the anvil region) to determine the most likely NO production scenario. The range of values of NO production per flash (or per meter of lightning channel length) that have been deduced from the model will be shown and compared with values of production in the literature that have been deduced from observed NO spikes and from anvil flux calculations. Results show that on a per flash basis, IC flashes are nearly as productive of NO as CG flashes. This result simplifies the lightning parameterization for global models (ie., an algorithm for estimating the IC/CG ratio is not necessary). Vertical profiles of lightning NOx mass at the end of the 3-D storm simulations have been summarized to yield suggested profiles for use in global models. Estimates of mean NO production per flash vary by a factor of three from one simulated storm to another. When combined with the global flash rate of 44 flashes per second from NASA's Optical Transient Detector (OTD) measurements, these estimates and the results from other techniques yield global NO production rates of2-9 TgN/year. Simulations of the photochemistry over the 24 hours following a storm has been performed to determine the additional ozone production which can be attributed to lightning NO. Convective transport of HOx precursors leads to the generation of a HOx plume which substantially aids the downstream ozone production.
This paper examines the ability of coupled global climate models to predict decadal variability of Arctic and Antarctic sea ice. We analyze decadal hindcasts/predictions of 11 Coupled Model Intercomparison Project Phase 5 (CMIP5) models. Decadal hindcasts exhibit a large multimodel spread in the simulated sea ice extent, with some models deviating significantly from the observations as the predicted ice extent quickly drifts away from the initial constraint. The anomaly correlation analysis between the decadal hindcast and observed sea ice suggests that in the Arctic, for most models, the areas showing significant predictive skill become broader associated with increasing lead times. This area expansion is largely because nearly all the models are capable of predicting the observed decreasing Arctic sea ice cover. Sea ice extent in the North Pacific has better predictive skill than that in the North Atlantic (particularly at a lead time of 3-7 years), but there is a reemerging predictive skill in the North Atlantic at a lead time of 6-8 years. In contrast to the Arctic, Antarctic sea ice decadal hindcasts do not show broad predictive skill at any timescales, and there is no obvious improvement linking the areal extent of significant predictive skill to lead time increase. This might be because nearly all the models predict a retreating Antarctic sea ice cover, opposite to the observations. For the Arctic, the predictive skill of the multi-model ensemble mean outperforms most models and the persistence prediction at longer timescales, which is not the case for the Antarctic. Overall, for the Arctic, initialized decadal hindcasts show improved predictive skill compared to uninitialized simulations, although this improvement is not present in the Antarctic.
This paper provides initial results from a multi-model ensemble analysis based on the volc-pinatubo-full experiment performed within the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The volc-pinatubo-full experiment is based on an ensemble of volcanic forcing-only climate simulations with the same volcanic aerosol dataset across the participating models (the 1991–1993 Pinatubo period from the CMIP6-GloSSAC dataset). The simulations are conducted within an idealized experimental design where initial states are sampled consistently across models from the CMIP6-piControl simulation providing unperturbed preindustrial background conditions. The multi-model ensemble includes output from an initial set of six participating Earth system models (CanESM5, GISS-E2.1-G, IPSL-CM6A-LR, MIROC-E2SL, MPI-ESM1.2-LR and UKESM1). The results show overall good agreement between the different models on the global and hemispheric scales concerning the surface climate responses, thus demonstrating the overall effectiveness of VolMIP's experimental design. However, small yet significant inter-model discrepancies are found in radiative fluxes, especially in the tropics, that preliminary analyses link with minor differences in forcing implementation; model physics, notably aerosol–radiation interactions; the simulation and sampling of El Niño–Southern Oscillation (ENSO); and, possibly, the simulation of climate feedbacks operating in the tropics. We discuss the volc-pinatubo-full protocol and highlight the advantages of volcanic forcing experiments defined within a carefully designed protocol with respect to emerging modelling approaches based on large ensemble transient simulations. We identify how the VolMIP strategy could be improved in future phases of the initiative to ensure a cleaner sampling protocol with greater focus on the evolving state of ENSO in the pre-eruption period.
The possibility of using the information from simultaneous equivalent Large Eddy Simulations (LAS) for improving the subgrid scale modeling is investigated. An ensemble average dynamic model is proposed as an alternative to the usual spatial average versions. It is shown to be suitable independently of the existence of any homogeneity directions, and its formulation is thus universal. The ensemble average dynamic model is shown to give very encouraging results for as few as 16 simultaneous LES's.
An ensemble of land surface models and forcing data was developed to assess variability in SWE estimation over North America. In this study, the ensemble output was used to assess how SWE uncertainty impacts stream flow estimation. The analysis was conducted by major basins of North America over the 2009-2017 time period.
Recent changes in the climate system have led to growing concern, especially in communities which are highly vulnerable to resource shortages and weather extremes. There is an urgent need for better climate information to develop solutions and strategies for adapting to a changing climate. Climate models provide excellent tools for studying the current state of climate and making future projections. However, these models are subject to biases created by structural uncertainties. Performance metrics-or the systematic determination of model biases-succinctly quantify aspects of climate model behavior. Efforts to standardize climate model experiments and collect simulation data-such as the Coupled Model Intercomparison Project (CMIP)-provide the means to directly compare and assess model performance. Performance metrics have been used to show that some models reproduce present-day climate better than others. Simulation data from multiple models are often used to add value to projections by creating a consensus projection from the model ensemble, in which each model is given an equal weight. It has been shown that the ensemble mean generally outperforms any single model. It is possible to use unequal weights to produce ensemble means, in which models are weighted based on performance (called "intelligent" ensembles). Can performance metrics be used to improve climate projections? Previous work introduced a framework for comparing the utility of model performance metrics, showing that the best metrics are related to the variance of top-of-atmosphere outgoing longwave radiation. These metrics improve present-day climate simulations of Earth's energy budget using the "intelligent" ensemble method. The current project identifies several approaches for testing whether performance metrics can be applied to future simulations to create "intelligent" ensemble-mean climate projections. It is shown that certain performance metrics test key climate processes in the models, and that these metrics can be used to evaluate model quality in both current and future climate states. This information will be used to produce new consensus projections and provide communities with improved climate projections for urgent decision-making.
Recent changes in the climate system have led to growing concern, especially in communities which are highly vulnerable to resource shortages and weather extremes. There is an urgent need for better climate information to develop solutions and strategies for adapting to a changing climate. Climate models provide excellent tools for studying the current state of climate and making future projections. However, these models are subject to biases created by structural uncertainties. Performance metrics-or the systematic determination of model biases-succinctly quantify aspects of climate model behavior. Efforts to standardize climate model experiments and collect simulation data-such as the Coupled Model Intercomparison Project (CMIP)-provide the means to directly compare and assess model performance. Performance metrics have been used to show that some models reproduce present-day climate better than others. Simulation data from multiple models are often used to add value to projections by creating a consensus projection from the model ensemble, in which each model is given an equal weight. It has been shown that the ensemble mean generally outperforms any single model. It is possible to use unequal weights to produce ensemble means, in which models are weighted based on performance (called "intelligent" ensembles). Can performance metrics be used to improve climate projections? Previous work introduced a framework for comparing the utility of model performance metrics, showing that the best metrics are related to the variance of top-of-atmosphere outgoing longwave radiation. These metrics improve present-day climate simulations of Earth's energy budget using the "intelligent" ensemble method. The current project identifies several approaches for testing whether performance metrics can be applied to future simulations to create "intelligent" ensemble-mean climate projections. It is shown that certain performance metrics test key climate processes in the models, and that these metrics can be used to evaluate model quality in both current and future climate states. This information will be used to produce new consensus projections and provide communities with improved climate projections for urgent decision-making.
The GEOS-5 atmospheric model is being developed as a weather-and-climate capable model. It must perform well in assimilation mode as well as in weather and climate simulations and forecasts and in coupled chemistry-climate simulations. In developing GEOS-5, attention has focused on the representation of moist processes. The moist physics package uses a single phase prognostic condensate and a prognostic cloud fraction. Two separate cloud types are distinguished by their source: "anvil" cloud originates in detraining convection, and large-scale cloud originates in a PDF-based condensation calculation. Ice and liquid phases for each cloud type are considered. Once created, condensate and fraction from the anvil and statistical cloud types experience the same loss processes: evaporation of condensate and fraction, auto-conversion of liquid or mixed phase condensate, sedimentation of frozen condensate, and accretion of condensate by falling precipitation. The convective parameterization scheme is the Relaxed Arakawa-Schubert, or RAS, scheme. Satellite data are used to evaluate the performance of the moist physics packages and help in their tuning. In addition, analysis of and comparisons to cloud-resolving models such as the Goddard Cumulus Ensemble model are used to help improve the PDFs used in the moist physics. The presentation will show some of our evaluations including precipitation diagnostics.
Increasing temperature trends are expected to impact yields of major field crops by affecting various plant processes, such as phenology and growth. However, future projections, especially at the global scale, do not consider any agronomic adaptation in farming practices. We use an ensemble of seven Global Gridded Crop Models (GGCMs) contributing to the AgMIP-Gridded project to quantify the impacts and adaptation potential of field crops under increasing temperature. We study how uniform warming scenarios up to 6 K affect the productivity and growing period duration of five major crops. The design of the experiment aims at understanding adaptation measures targeted to temperature-driven impacts and therefore leaves out changes in any other climate variable. We find that, without adaptation the dominant effect of temperature increase is to shorten the growing period and to reduce yields and production, consistently across crops and regions. We then test the potential of two agronomic measures to contrast warming-induced yield reduction, assuming that (i) cultivars with adjusted phenology would be used to regain the reference growing period duration; (ii) rainfed systems would be converted to irrigated. Despite substantial uncertainties in model parametrization and little agreement in spatial patterns, at the global aggregation, model results are robust. We find that up to 3 K of temperature increase, production losses could be fully compensated by maintaining the original crop growing period. Irrigation would also compensate production losses, but would not reduce the temperature impacts. Across regions, the model ensemble estimates larger adaptation potentials from unchanged growing periods in continental and temperate regions, than in tropical and arid regions, where also irrigation has respectively little effects and availability.
Decadal trends and interannual variations in the hydroxyl radical (OH), while poorly constrained at present, are critical for understanding the observed evolution of atmospheric methane (CH4). Through analyzing the OH fields simulated by the model ensemble of the Chemistry-Climate Model Initiative (CCMI), we find (1) the negative OH anomalies during the El Niño years mainly corresponding to the enhanced carbon monoxide (CO) emissions from biomass burning and (2) a positive OH trend during 1980-2010 dominated by the elevated primary production and the reduced loss of OH due to decreasing CO after 2000. Both two-box model inversions and variational 4D inversions suggest that ignoring the negative anomaly of OH during the El Niño years leads to a large overestimation of the increase in global CH4 emissions by up to10±3Tg yr-1to match the observed CH4 increase over these years. Not accounting for the increasing OH trends given by the CCMI models leads to an underestimation of the CH4 emission increase by 23±9Tg yr-1from 1986 to 2010. The variational inversion estimated CH4 emissions show that the tropical regions contribute most to the uncertainties related to OH. This study highlights the significant impact of climate and chemical feedbacks related to OH on the top-down estimates of the global CH4 budget.
Early summer (May–June–July; MJJ) droughts over the Central Plains are often caused by atmospheric ridging, but it is uncertain if these events will increase in frequency or if their influence on drought severity will change in a warming world. Here, we use tree-ring based reconstructions (1500–2020 CE) of MJJ ridging and 0–200 cm soil moisture with six CMIP6 model ensembles to investigate the response of Central Plains drought dynamics to a moderate warming scenario (SSP2-4.5). By the end of the 21st century (2071–2100), precipitation increases in most models during the preceding months (February–March–April), especially over the northern part of the Central Plains, while changes during MJJ are non-robust. By contrast, vapor pressure deficit increases strongly in all models, resulting in five of the six models projecting robust median soil moisture drying and all six models projecting more rapid seasonal soil moisture declines during the transition into the summer. Major ridging events increase in frequency in some models, and there is strong agreement across all models that when ridging events do occur, they will cause more severe soil moisture drought and seasonal drying at the end of the 21st century. The median multi-model response also indicates, by the end of the 21st century, that the Central Plains will experience a three-fold increase in the risk of drought events equivalent to the most extreme droughts of the last 500 years. Our results demonstrate that even moderate warming is likely to increase early summer soil moisture drought severity and risk over the Central Plains, even in the absence of robust precipitation declines, and that drought responses to major atmospheric ridging events will be significantly stronger.
Concerns over climate change are motivated in large part because of their impact on human society. Assessing the effect of that uncertainty on specific potential impacts is demanding, since it requires a systematic survey over both climate and impacts models. We provide a comprehensive evaluation of uncertainty in projected crop yields for maize, spring and winter wheat, rice, and soybean, using a suite of nine crop models and up to 45 CMIP5 and 34 CMIP6 climate projections for three different forcing scenarios. To make this task computationally tractable, we use a new set of statistical crop model emulators. We find that climate and crop models contribute about equally to overall uncertainty. While the ranges of yield uncertainties under CMIP5 and CMIP6 projections are similar, median impact in aggregate total caloric production is typically more negative for the CMIP6 projections (+1% to −19%) than for CMIP5 (+5% to −13%). In the first half of the 21st century and for individual crops is the spread across crop models typically wider than that across climate models, but we find distinct differences between crops: globally, wheat and maize uncertainties are dominated by the crop models, but soybean and rice are more sensitive to the climate projections. Climate models with very similar global mean warming can lead to very different aggregate impacts so that climate model uncertainties remain a significant contributor to agricultural impacts uncertainty. These results show the utility of large-ensemble methods that allow comprehensively evaluating factors affecting crop yields or other impacts under climate change. The crop model ensemble used here is unbalanced and pulls the assumption that all projections are equally plausible into question. Better methods for consistent model testing, also at the level of individual processes, will have to be developed and applied by the crop modeling community.
In situ hydrologic monitoring over regions most susceptible to food insecurity can be a challenge in current times due to various socio-economic and political issues in combination with environmental factors such as ongoing famine or drought. Hydrologic monitoring and initializing forecasts based on remotely sensed and analyzed data can contribute significantly to early warning in such regions. Routine hydrologic forecasts, as provided by NASA’s Hydrologic Forecasting and Analysis System (NHyFAS), are a recent addition to early warning systems. A custom instance of NHyFAS, termed FLDAS-Forecast, is used by FEWS NET’s Land Data Assimilation System (FLDAS). The FLDAS-Forecast’s dynamic forecasting component was originally set up with Goddard Earth Observing System (GEOS) forecast inputs and has been recently expanded with precipitation forecast forcing from the North American Multi-Model Ensemble (NMME). This paper describes the improvements in seasonal hydrologic forecasts produced with this updated system. Evaluations in this study focus on soil moisture across southern Africa’s growing season. Soil moisture forecasts are benchmarked and evaluated relative to climatology-based forecasts and historic runs, which are driven by observation-based meteorological forcing fields, and they are verified with remotely sensed observations of soil moisture and vegetation. Through multiple deterministic and probabilistic skill assessments, we show that using the larger ensemble of NMME precipitation inputs in the forecast system results in higher quality hydrologic forecasts than are allowed by climatology- or GEOS-only-based forecasts. Further, the near-real-time NMME-based rootzone soil moisture forecasts were able to correctly predict developing drought conditions over southern Africa through late 2019 and into early 2020.
The Goddard Convective-Stratiform Heating (CSH) algorithm, used to estimate cloud heating in support of the Tropical Rainfall Measuring Mission (TRMM), is upgraded in support of the Global Precipitation Measurement mission (GPM). The algorithm is required to use look-uptables (LUTs) from cloud-resolving model (CRM) simulations from the Goddard Cumulus Ensemble model (GCE). This paper will present the heating retrievals from the Goddard CSH algorithm in the TRMM and GPM using precipitation products (rainfall, radar reflectivity).
The objective of the present study is to evaluate the dynamic procedure in LES of stratocumulus topped atmospheric boundary layer and assess the relative importance of subgrid-scale modeling, cloud microphysics and radiation modeling on the predictions. The simulations will also be used to gain insight into the processes leading to cloud top entrainment instability and cloud breakup. In this report we document the governing equations, numerical schemes and physical models that are employed in the Goddard Cumulus Ensemble model (GCEM3D). We also present the subgrid-scale dynamic procedures that have been implemented in the GCEM3D code for the purpose of the present study.