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.
We present an open-source Python-based framework that can be used for objective evaluation of Atmospheric Rivers (ARs) in climate models. With a comprehensive pool of pre-defined metrics, this package is useful for systematic AR evaluation across a range of model simulations, such as the MIPs in CMIP5/6 models and large ensemble experiments.
This study applies an idealized modeling framework, alongside a Bayesian Markov chain Monte Carlo (MCMC) algorithm, to explore which combinations of upstream environmental conditions and cloud microphysical parameter values can produce a particular precipitation distribution over an idealized two-dimensional, bell-shaped mountain. Simulations focus on orographic precipitation produced when an atmospheric river interacts with topography. MCMC-based analysis reveals that different combinations of parameter values produce a similar precipitation distribution, with the most influential parameters being relative humidity (RH), horizontal wind speed (U), surface potential temperature (θ sfc ), and the snow fall speed coefficient (A s ). RH, U, and A s exhibit interdependence: changes in one or more of these factors can be mitigated by compensating changes in the other(s) to produce similar orographic precipitation rates. The results also indicate that the parameter sensitivities and relationships can vary for spatial subregions and given different environmental conditions. In particular, high θ sfc values are more likely to produce the target precipitation rate and spatial distribution, and thus the ensemble of simulations shows a preference for liquid precipitation at the surface. The results presented here highlight the complexity of orographic precipitation controls, and have implications for flood and water management, observational efforts, and climate change.
Future water scarcity is a global concern with significant impacts on the energy, water, and land (EWL) sectors. Countries in Latin America and the Caribbean (LAC) are significant producers of agricultural goods consumed around the world, so disruptions resulting from land and water scarcity in LAC are a source of global supply chain risk. Understanding where water scarcity in LAC could occur and what could drive that scarcity to emerge is critical for strategic resource management and planning, both regionally and globally. Assessing future water scarcity impacts in LAC is challenging given the multisector dynamic interactions among the EWL sectors and the multiple uncertainties acting across different spatial scales. We use scenario discovery to explore the drivers of future water scarcity, considering assumptions related to climate, socioeconomics, and the EWL sectors. Understanding what factors drive uncertainty across outcomes can help stakeholders in targeting long-term solutions and future data collection efforts. To illuminate these dynamics, we use a database consisting of a large ensemble of scenarios representing diverse worlds simulated using the Global Change Analysis Model. We quantify future water scarcity and its economic impacts across the scenarios using three metrics: (1) physical water scarcity, (2) water price, and (3) crop profit. We find that physical water scarcity and water price are driven primarily by reservoir storage capacity assumptions, highlighting the importance of strategic water infrastructure development in maintaining future water availability and accessibility. Crop profit is driven by a combination of water supply and demand assumptions, which emphasizes the complex nature of EWL multisector dynamics. While most of LAC is poised to have abundant land and water resources available for development, water basins in Mexico and along the Pacific coast of South America experience high levels of severity and uncertainty across the future scenarios for at least one of the metrics. We find that the drivers of extreme scarcity vary spatially and across the metrics, further highlighting the heterogeneity of the region and the importance of considering multiple metrics to assess water scarcity vulnerability.
Mountains play an outsized role in water resource availability, and the amount and timing of water they provide depend strongly on temperature. To that end, we ask the question: How well are atmospheric models capturing mountain temperatures? We synthesize results showing that high-resolution, regionally relevant climate models produce 2-m air temperature (T2m) measurements colder than what is observed (a “cold bias”), particularly in snow-covered midlatitude mountain ranges during winter. We find common cold biases in 44 studies across global mountain ranges, including single-model and multimodel ensembles. We explore the factors driving these biases and examine the physical mechanisms, data limitations, and observational uncertainties behind T2m. Our analysis suggests that the biases are genuine and not due to observation sparsity or resolution mismatches. Cold biases occur primarily on mountain peaks and ridges, whereas valleys are often warm biased. Our literature review suggests that increasing model resolution does not clearly mitigate the bias. By analyzing data from the Surface Atmosphere Integrated Field Laboratory (SAIL) field campaign in the Colorado Rocky Mountains, we test various hypotheses related to cold biases and find that local wind circulations, longwave (LW) radiation, and surface-layer parameterizations contribute to the T2m biases in this particular location. We conclude by emphasizing the value of coordinated model evaluation and development efforts in heavily instrumented mountain locations for addressing the root cause(s) of T2m biases and improving predictive understanding of mountain climates.
Mountains play an outsized role in water resource availability, and the amount and timing of water they provide depend strongly on temperature. To that end, we ask the question: How well are atmospheric models capturing mountain temperatures? We synthesize results showing that high-resolution, regionally relevant climate models produce 2-m air temperature (T2m) measurements colder than what is observed (a “cold bias”), particularly in snow-covered midlatitude mountain ranges during winter. We find common cold biases in 44 studies across global mountain ranges, including single-model and multimodel ensembles. We explore the factors driving these biases and examine the physical mechanisms, data limitations, and observational uncertainties behind T2m. Our analysis suggests that the biases are genuine and not due to observation sparsity or resolution mismatches. Cold biases occur primarily on mountain peaks and ridges, whereas valleys are often warm biased. Our literature review suggests that increasing model resolution does not clearly mitigate the bias. By analyzing data from the Surface Atmosphere Integrated Field Laboratory (SAIL) field campaign in the Colorado Rocky Mountains, we test various hypotheses related to cold biases and find that local wind circulations, longwave (LW) radiation, and surface-layer parameterizations contribute to the T2m biases in this particular location. We conclude by emphasizing the value of coordinated model evaluation and development efforts in heavily instrumented mountain locations for addressing the root cause(s) of T2m biases and improving predictive understanding of mountain climates.
Earth's hydrological cycle is expected to intensify in response to global warming, with a “wet-gets-wetter, dry-gets-drier” response anticipated over the ocean. Subtropical regions (∼15°–30°N/S) are predicted to become drier, yet proxy evidence from past warm climates suggests these regions may be characterized by wetter conditions. Here we use an integrated data-modeling approach to reconstruct global and zonal-mean rainfall patterns during the early Eocene (∼56–48 million years ago). The Deep-Time Model Intercomparison Project (DeepMIP) model ensemble indicates that the mid- (30°–60°N/S) and high-latitudes (>60°N/S) are characterized by a thermodynamically dominated hydrological response to warming and overall wetter conditions. The tropical band (0°–15°N/S) is also characterized by wetter conditions, with several DeepMIP models simulating narrowing of the Inter-Tropical Convergence Zone. However, the latter is not evident from the proxy data. The subtropics are characterized by negative precipitation-evaporation anomalies (i.e., drier conditions) in the DeepMIP models, but there is surprisingly large inter-model variability in mean annual precipitation (MAP). Intriguingly, we find that models with weaker meridional temperature gradients (e.g., CESM, GFDL) are characterized by a reduction in subtropical moisture divergence, leading to an increase in MAP. These model simulations agree more closely with our new proxy-derived precipitation reconstructions and other key climate metrics and imply that the early Eocene was characterized by reduced subtropical moisture divergence. If the meridional temperature gradient was even weaker than suggested by those DeepMIP models, circulation-induced changes may have outcompeted thermodynamic changes, leading to wetter subtropics. This highlights the importance of accurately reconstructing zonal temperature gradients when reconstructing past rainfall patterns.
Climate models often exhibit biases in their precipitation predictions, particularly underestimating high-intensity events and overestimating low precipitation. Deep learning approaches offer promising solutions, but their epistemic uncertainty associated with a deep learning–based bias correction method has not previously been quantified for reliable downstream climate impact studies. While methods for capturing the epistemic uncertainty in deep learning frameworks exist, there is currently no consensus on the best method. In this work, we compare three uncertainty quantification (UQ) methods—Deep Ensembles (DEns), Monte Carlo Dropout (MCD), and Flipout—by assessing the reliability of their uncertainty estimates using standard measures such as sharpness and calibration. These UQ methods are applied to an existing deep learning precipitation bias correction model known as UFNet: a coupled U-Net and fully connected neural network. The methods utilized to assess the models’ uncertainties are 1) calibration, which ensures that the expected probabilities of the model align with reality and 2) sharpness, which is a measure of the precision of the model’s probabilistic predictions. Of the three UQ methods evaluated, the DEns and MCD methods demonstrated the best-calibrated performance (expected calibration error of 0.36 and 0.35, respectively), compared to Flipout (0.58). In contrast, Flipout had the sharpest predictions and the highest metric performance in bias correcting precipitation—especially for higher-order moments such as kurtosis with a spatial correlation of 72% compared to 32% and 55% spatial correlation for DEns and MCD, respectively. Of the three UQ methods, MCD was found to be the most suitable method for UQ purposes based on its calibration, sharpness, and computational requirements.
The objective of the North American Land Data Assimilation System (NLDAS) is to provide best available estimates of near-surface meteorological conditions and soil hydrological status for the continental United States. To support the ongoing efforts to develop data assimilation (DA) capabilities for NLDAS, the results of Gravity Recovery and Climate Experiment (GRACE) DA implemented in a manner consistent with NLDAS development are presented. Following previous work, GRACE terrestrial water storage (TWS) anomaly estimates are assimilated into the NASA Catchment land surface model using an ensemble smoother. In contrast to many earlier GRACE DA studies, a gridded GRACE TWS product is assimilated, spatially distributed GRACE error estimates are accounted for, and the impact that GRACE scaling factors have on assimilation is evaluated. Comparisons with quality-controlled in situ observations indicate that GRACE DA has a positive impact on the simulation of unconfined groundwater variability across the majority of the eastern United States and on the simulation of surface and root zone soil moisture across the country. Smaller improvements are seen in the simulation of snow depth, and the impact of GRACE DA on simulated river discharge and evapotranspiration is regionally variable. The use of GRACE scaling factors during assimilation improved DA results in the western United States but led to small degradations in the eastern United States. The study also found comparable performance between the use of gridded and basin averaged GRACE observations in assimilation. Finally, the evaluations presented in the paper indicate that GRACE DA can be helpful in improving the representation of droughts.
Extreme weather events have significant consequences, dominating the impact of climate on society. While high-resolution weather models can forecast many types of extreme events on synoptic timescales, long-term climatological risk assessment is an altogether different problem. A once-in-a-century event takes, on average, 100 years of simulation time to appear just once, far beyond the typical integration length of a weather forecast model. Therefore, this task is left to cheaper, but less accurate, low-resolution or statistical models. But there is untapped potential in weather model output: despite being short in duration, weather forecast ensembles are produced multiple times a week. Integrations are launched with independent perturbations, causing them to spread apart over time and broadly sample phase space. Collectively, these integrations add up to thousands of years of data. We establish methods to extract climatological information from these short weather simulations. Using ensemble hindcasts by the European Center for Medium-range Weather Forecasting archived in the subseasonal-to-seasonal (S2S) database, we characterize sudden stratospheric warming (SSW) events with multi-centennial return times. Consistent results are found between alternative methods, including basic counting strategies and Markov state modeling. By carefully combining trajectories together, we obtain estimates of SSW frequencies and their seasonal distributions that are consistent with reanalysis-derived estimates for moderately rare events, but with much tighter uncertainty bounds, and which can be extended to events of unprecedented severity that have not yet been observed historically. These methods hold potential for assessing extreme events throughout the climate system, beyond this example of stratospheric extremes.
Abstract. A new set of standalone parameterizations is presented for simulating the injection, evolution, and radiative forcing by stratospheric volcanic aerosols against an idealized Held–Suarez–Williamson (HSW) atmospheric background in the Energy Exascale Earth System Model version 2 (E3SMv2). In this model configuration (HSW with enabled volcanism, HSW-V), sulfur dioxide (SO2) and ash are injected into the atmosphere with a specified profile in the vertical, and they proceed to follow a simple exponential decay. The SO2 decay is modeled as a perfect conversion to a long-living sulfate aerosol which persists in the stratosphere. All three species are implemented as tracers in the model framework and are transported by the dynamical core's advection algorithm. The aerosols contribute simultaneously to local heating of the stratosphere and cooling of the surface by a simple plane-parallel Beer–Lambert law applied on two zonally symmetric radiation broadbands in the longwave and shortwave ranges. It is shown that the implementation parameters can be tuned to produce realistic temperature anomaly signatures of large volcanic events. In particular, results are shown for an ensemble of runs that mimic the volcanic eruption of Mt. Pinatubo in 1991. The design requires no coupling to microphysical subgrid-scale parameterizations and thus approaches the computational affordability of prescribed aerosol forcing strategies. The idealized simulations contain a single isolated volcanic event against a statistically uniform climate, where no background aerosols or other sources of externally forced variability are present. HSW-V represents a simpler-to-understand tool for the development of climate source-to-impact attribution methods.
In the domain of predicting land surface fluxes, models are used to bring data from large observation networks and satellite remote sensing together to make predictions about present and future states of the Earth. Characterizing the uncertainty about such predictions is a complex process and one that is not yet fully understood. Uncertainty exists about initialization, measurement and interpolation of input variables; model parameters; model structure; and mixed spatial and temporal supports. Multiple models or structures often exist to describe the same processes. Uncertainty about structure is currently addressed by running an ensemble of different models and examining the distribution of model outputs. To illustrate structural uncertainty, a multi-model ensemble experiment we have been conducting using the Terrestrial Observation and Prediction System (TOPS) will be discussed. TOPS uses public versions of process-based ecosystem models that use satellite-derived inputs along with surface climate data and land surface characterization to produce predictions of ecosystem fluxes including gross and net primary production and net ecosystem exchange. Using the TOPS framework, we have explored the uncertainty arising from the application of models with different assumptions, structures, parameters, and variable definitions. With a small number of models, this only begins to capture the range of possible spatial fields of ecosystem fluxes. Few attempts have been made to systematically address the components of uncertainty in such a framework. We discuss the characterization of uncertainty for this approach including both quantifiable and poorly known aspects.
Abstract To improve the representation of microphysical processes in convective clouds and their interaction with aerosol and stratiform clouds, a two‐moment convective microphysics parameterization (CMP) scheme developed by Song and Zhang (2011, https://doi.org/10.1029/2010jd014833 ) is upgraded and implemented in E3SM. The new developments include: (a) implementing a parameterization for graupel to enhance the representation of ice‐phase microphysical processes; (b) representing the impact of spatial inhomogeneity of cloud droplets in cumulus ensembles on autoconversion and accretion processes to improve the representation of warm‐rain microphysical processes; (c) implementing a comprehensive Bergeron process parameterization to better represent mixed‐phase microphysical processes; and (d) representing the interactions between ice‐phase microphysics and cloud thermodynamics. Simulations show that the cloud microphysical properties simulated by the CMP are generally in good agreement with observations. It reasonably simulates the changes in droplets effective radius related to precipitation formation in convective clouds, as identified from satellite observations. It also successfully simulates the contrast in these processes between maritime and continental clouds, demonstrating its capability to simulate the impact of aerosols on convection. Analyses of the impact of CMP on climate mean state simulation demonstrate that the CMP slightly improves the simulations of precipitation, cloud macrophysical properties, longwave cloud radiative forcing, zonal wind, and temperature. However, a degradation in shortwave cloud radiative forcing occurs.
Predicting how regional precipitation will respond to future warming is among the most challenging undertaking in climate change projection. Despite sustained efforts to improve modeling and understanding of precipitation, the overall uncertainty in projecting regional precipitation has not been reduced substantially. In this work, the potential for more robust precipitation projections is demonstrated through the use of discriminating metrics to subsample the multimodel ensemble. Using a two-dimensional metric of precipitation and its relationship with large-scale circulation indices in East Asia, 31 models in the Coupled Model Intercomparison Project Phase 5 (CMIP5) are classified into three groups. Models in the top performing group projected statistically significant increasing trends in precipitation and the regional precipitation patterns are more similar to each other than to the patterns in the bottom performing group. In contrast, models in the bottom performing group projected diverse responses, with overall small drying or no significant trends in precipitation.
Aerosol-cloud interactions (ACI) in warm clouds are the primary source of uncertainty in effective radiative forcing (ERF) during the historical period and, by extension, inferred climate sensitivity. The ERF due to ACI (ERFaci) is composed of the radiative forcing due to changes in cloud microphysics and cloud adjustments to microphysics. Here, we examine the processes that drive ERFaci using a perturbed parameter ensemble (PPE) hosted in CAM6. Observational constraints on the PPE result in substantial constraints in the response of cloud microphysics and macrophysics to anthropogenic aerosol, but only minimal constraint on ERFaci. Examination of cloud and radiation processes in the PPE reveal buffering of ERFaci by the interaction of precipitation efficiency and radiative susceptibility.
Abstract Biomass burning is a significant source of aerosol emissions in some regions and has a considerable impact on regional climate. Earth system model simulations indicate that increased biomass burning aerosol emissions contributed to statistically significant decreases in tropical precipitation over the 20th century. In this study, we use the Community Earth System Model version 1 Large Ensemble (CESM1-LENS) experiment to evaluate the mechanisms by which biomass burning aerosol contributed to decreased tropical precipitation, with a focus on South America and Southeast Asia. We analyze the all-but-one forcing simulations in which biomass burning aerosol emissions are held constant while other forcings (e.g., greenhouse gas concentrations) vary throughout the 20th century. This allows us to isolate the influence of biomass burning aerosol on processes that contribute to decreasing precipitation, including cloud microphysics, the radiative effects of absorbing aerosol particles, and alterations in regional circulation. We also show that the 20th century reduction in precipitation identified in the CESM1-LENS historical and biomass burning experiments is consistent across Coupled Model Intercomparison Project Phase 5 models with interactive aerosol schemes and the CESM2 single-forcing experiment. Our results demonstrate that higher concentrations of biomass burning aerosol increases the quantity of cloud condensation nuclei and cloud droplets, limiting cloud droplet size and precipitation formation. Additionally, absorbing aerosols (e.g., black carbon) contribute to a warmer cloud layer, which promotes cloud evaporation, increases atmospheric stability, and alters regional circulation patterns. Corresponding convectively coupled circulation responses, particularly over the tropical Andes, contribute to further reducing the flow of moisture and moisture convergence over tropical land. These results elucidate the processes that affect the water cycle in regions prone to biomass burning and inform our understanding of how future changes in aerosol emissions may impact tropical precipitation over the 21st century.
Aerosol perturbations over selected land regions are imposed in Version-4 of the Goddard Earth Observing System (GEOS-4) general circulation model (GCM) to assess the influence of increasing aerosol concentrations on regional circulation patterns and precipitation in four selected regions: India, Africa, and North and South America. Part 1 of this paper addresses the responses to aerosol perturbations in India and Africa. This paper presents the same for aerosol perturbations over the Americas. GEOS-4 is forced with prescribed aerosols based on climatological data, which interact with clouds using a prognostic scheme for cloud microphysics including aerosol nucleation of water and ice cloud hydrometeors. In clear-sky conditions the aerosols interact with radiation. Thus the model includes comprehensive physics describing the aerosol direct and indirect effects on climate (hereafter ADE and AIE respectively). Each simulation is started from analyzed initial conditions for 1 May and was integrated through June-July-August of each of the six years: 1982 1987 to provide a 6-ensemble set. Results are presented for the difference between simulations with double the climatological aerosol concentration and one-half the climatological aerosol concentration for three experiments: two where the ADE and AIE are applied separately and one in which both the ADE and AIE are applied. The ADE and AIE both yield reductions in net radiation at the top of the atmosphere and surface while the direct absorption of shortwave radiation contributes a net radiative heating in the atmosphere. A large net heating of the atmosphere is also apparent over the subtropical North Atlantic Ocean that is attributable to the large aerosol perturbation imposed over Africa. This atmospheric warming and the depression of the surface pressure over North America contribute to a northward shift of the inter-Tropical Convergence Zone over northern America, an increase in precipitation over Central America and the Caribbean, and an enhancement of convergence in the North American monsoon region.
South Africa currently possesses no operational ocean forecasting system for the purpose of predicting ocean state variables including temperature,salinity and velocity. Substantial initial efforts towards this goal have been made and resulted in a system using a regional Hybrid Coordinate Ocean Model (HYCOM) along with the Ensemble Optimal Interpolation (EnOI)assimilation scheme. Assimilating only sea surface temperature (SST) observations from the Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) product into the system resulted in a degraded forecast. Aiming to address this, Climate Change Initiative (CCI) SSTs are assimilated into the system in an effort to improve the forecast skill. Observation errors in the assimilated product are used in the EnOI to determine whether more confidence should be placed in the model or observations in producing the analysis, but overconfidence in observations can shock the model and result in failure. To tweak the impact of the assimilation, a scaling factor is applied in the assimilation code. A scaling factor of 25 was found to produce a favourable result with lowest mean root mean square error (RMSE;1.098°C) between the model and observations over time. Postulating the error to be overconfident, a floor value is introduced in order to set a minimum value for the observation error thereby reducing confidence in the observations. These experiments fared less favourably with a floor value of 0.5 and a scaling factor of 15 producing the best mean RMSE (1.118°C).
In 2022, the United Nations launched the Early Warnings for All (EW4ALL) Program to establish global early warning systems by 2027. To assess the impact of the substantial $3.1 billion annual investment over five years, EW4ALL will consider factors that will require national coordination for the data needed for these assessments. In 2023, Ecuador was identified as one of the world's most climate-vulnerable countries, emphasizing the need to enhance its early warning systems. In 2020, the SERVIR Amazonia hub implemented the GEOGLOWS streamflow forecast service in collaboration with Ecuador's national meteorological agency (INAMHI). GEOGloWS provides 15-day ensemble forecasts and 80 years of historical streamflow data for every river worldwide through a free web service. The World Meteorological Organization has recognized this initiative as essential in contributing to the UN's call to ensure an 'Early Warning for All' by 2027. In 2023, as part of NASA's continuous efforts to fund research for Policy-Relevant Implementations, an economic impact assessment (EIA) was performed to understand the potential socioeconomic benefits of Early streamflow predictions in Ecuador using the GEOGLOWS service. Preliminary findings highlighted that gaps remain in effectively integrating socioeconomic and Earth observation (EO) data to capture the total value of these predictions. Implementing GEOGLOWS has led to valuable hydrological forecasts; however, the total economic benefits have yet to be documented. This study addresses the gaps and makes recommendations for future work that should focus on capturing the socio-economic benefits and costs associated with these forecasts, including their impact on decision-making at national and local levels. Despite the daily use of GEOGLOWS by key figures, including the President of Ecuador, the need for comprehensive recommendations and assessments is urgent.