Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble prediction”

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.

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Combining Emergent Constraints for Climate Sensitivity

A method is proposed for combining information from several emergent constraints into a probabilistic estimate for a climate sensitivity proxy Y such as equilibrium climate sensitivity (ECS). The method is based on fitting a multivariate Gaussian PDF for Y and the emergent constraints using an ensemble of global climate models (GCMs); it can be viewed as a form of multiple linear regression of Y on the constraints. The method accounts for uncertainties in sampling this multidimensional PDF with a small number of models, for observational uncertainties in the constraints, and for overconfidence about the correlation of the constraints with the climate sensitivity. Its general form (Method C) accounts for correlations between the constraints. Method C becomes less robust when some constraints are too strongly related to each other; this can be mitigated using regularization approaches such as ridge regression. An illuminating special case, Method U, neglects any correlations between constraints except through their mutual relationship to the climate proxy; it is more robust to small GCM sample size and is appealingly interpretable. These methods are applied to ECS and the climate feedback parameter using a previously published set of 11 possible emergent constraints derived from climate models in the Coupled Model Intercomparison Project (CMIP). The ±2σ posterior range of ECS for Method C with no overconfidence adjustment is 4.3 ± 0.7 K. For Method U with a large overconfidence adjustment, it is 4.0 ± 1.3 K. This study adds confidence to past findings that most constraints predict higher climate sensitivity than the CMIP mean.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO S2S Forecast System

Sea surface salinity (SSS) observations from space allow us to investigate if improved estimates of near-surface density stratification and associated mixing can positively impact seasonal to interannual variability of tropical Pacific Ocean dynamics as well as dynamical ENSO forecasts. For the first part of the presentation, we utilize our intermediate-complexity coupled model. Baseline experiments assimilate satellite sea level (multi-satellite gridded AVISO, 2013), SST (Reynolds et al., 2004), and in situ subsurface temperature and salinity observations (GTSPP NODC, 2006). These baseline experiments are then compared with experiments that additionally assimilate Aquarius (V5.0 Lilly and Lagerloef, 2008) and SMAP (V4.0 Fore et al., 2016) SSS. Twelve-month forecasts are initialized for each month from September 2011 to September 2017. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. This pattern enhances air/sea interaction and amplifies the equatorial Kelvin wave signal. We find that including satellite SSS significantly improves NINO3.4 sea surface temperature anomaly validation over most forecast lead times. We next assess how different satellite SSS products impact the validation of ENSO forecasts. SMAP V4 reduces the salty bias in the western Pacific and so is an improvement upon SMAP V2 and SMOS V2 (Boutin et al., 2017) has similar validation characteristics as a combination of Aquarius and SMAP V4. Next we shift to present results from the NASA GMAO Sub-seasonal to seasonal (S2S_v2.1) production coupled model (i.e. the same model that contributes ENSO forecasts to the North American Multi-Model Ensemble Experiment). From March to June 2015, the availability of two overlapping satellite SSS instruments, Aquarius and SMAP, allows a unique opportunity to compare and contrast forecasts initialized with the benefit of these two satellite SSS observation types. We assess the impact of satellite sea surface salinity (SSS) observations on dynamical ENSO forecasts for the big 2015 El Nino event. We will present distinct experiments for the overlap period that include 1) freely evolving SSS (i.e. no satellite SSS as the production system), 2) Aquarius, and 3) SMAP initialization. Our results show that using Aquarius slightly improves validation of the reanalysis (including sea level and temperature statistics). Our production system without SSS assimilation generated too warm forecasts for the 2015 El Nino from March initial conditions. Incorporating Aquarius into initialization of the coupled system leads to a deeper, more realistic MLD that acts to damp the downwelling Kelvin signal and slightly cool NINO3.4 SST. With Aquarius the forecasts better match the observed amplitude of the 2015 event. On the other hand, SMAP V2 relaxation generally degrades validation statistics. At forecast initialization, SMAP is much too salty within 10o of the equator, leading to deeper MLD east of 165W. This deeper MLD leads to over-damping of the downwelling signal (i.e. relative upwelling), in turn leading to relatively too cool ENSO forecasts.

Hackert, E.↗
DOE OSTI2021

Data-Driven Specification of Earth's Magnetic Field

As our society is increasingly dependent on technologies in space, predicting space weather becomes ever more important. We have come to rely on technologies such as cellphones, GPS, the Internet, and other commercial and military assets. In order to protect these systems a high priority is to develop better understanding of the harmful conditions of space weather encountered in the near-Earth radiation environment which are mostly dictated by the dynamics of energetic charged particles. The main objective of our proposal was to improve the specification and forecast of the particle distribution in the Earth’s inner magnetosphere using data assimilation together with machine learning (ML) algorithms in our ring current-atmosphere interactions model with self-consistent magnetic and electric fields (RAM-SCBE). Specifically, we explored the application of machine learning algorithms to emulate and improve, through infusion of observational data, the magnetic field (SCBE) of our ring current model. Our work advanced towards the characterization of the RAM-SCBE particle flux field and its potentially different simulation scenarios. We explored how the model behaves in a number of simulation scenarios in order to better understand how the model reacts to perturbations in the solar wind. We were able to capture this behavior through the simulation of an ensemble of RAM-SCB model solutions. The ensemble showed the important correlations between the particle fluxes and the magnetic field. We were able to start the training of a neural network (Autoencoder) to identify the dominant features of the magnetic field, and the main correlations (linear and non-linear) between the particle fluxes and the magnetic field. Our next steps consists of training a neural network to capture the time evolution of the dominant features and obtain a consistent surrogate model for the RAM-SCB.

42 ENGINEERING↗
DOE OSTI2020

Radiative heating rate profiles over the southeast Atlantic Ocean during the 2016 and 2017 biomass burning seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the first Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1) in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96° S, 14.35° W). Using observed profiles of temperature, humidity, and clouds from the field campaign alongside aerosol optical properties from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM), profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Idealized experiments using the RRTM with and without aerosols and a range of values for the single-scattering albedo (SSA) demonstrate that shortwave (SW) heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K d -1 depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA, and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave (LW) cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single-scattering albedo assumptions in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Surface Passivation by Excess Sulfur for Controlled Synthesis of Large, Thin SnS Flakes

Tin sulfide (SnS) is part of the group IV chalcogenides (SnX, GeX; X: S, Se), a family of anisotropic layered materials considered for thin film photovoltaics, optoelectronics, and valleytronics, and predicted narrow bandgap multiferroic materials. Large ultrathin SnS flakes, suitable for a variety of applications, are challenging to synthesize because of the enhanced surface reactivity due to the open layer structure of SnS which gives rise to a competition between lateral and vertical growth. Here we investigate the effects of added sulfur on passivating the surface and modifying the balance between lateral expansion and thickening of SnS flakes in chemical vapor transport from a SnS precursor. In addition, we investigate the growth of SnS flakes and compare the results of synthesis from a pure SnS precursor with growth, in which a slight excess of sulfur is supplied in different ways. Our results demonstrate that small amounts of excess sulfur can profoundly affect the size and thickness distributions of SnS flakes. The largest and thinnest flakes are obtained if (i) traces of sulfur are added and (ii) the sulfur source consists either of small Sx fragments released at high temperatures from the reactor walls, or atomic S supplied by sublimation from SnS 2 . The likely mechanism for the observed growth modifications is a transient surface passivation of SnS flakes during growth, which reduces the reactivity of the top facet of the flakes, limits vertical growth, and thus gives rise to ensembles with increased lateral size and reduced thickness.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Investigating the linkage between simulated precipitation climatology and ENSO-related precipitation anomaly based on multi-model and perturbed parameter ensembles

It remains a challenge for climate models to correctly capture the relationship between precipitation and ENSO. This studyexamines the linkage between the simulated precipitation climatology and ENSO-related precipitation anomaly during borealwinter based on the multi-model ensemble from the Atmospheric Model Intercomparison Project Phase 5 (AMIP5) and perturbedparameter ensemble (PPE) with the Beijing Climate Center (BCC) atmospheric model. The AMIP5 models have large biases insimulating the tropical precipitation anomaly during El Niño, such as the shifts of the Inter-Tropical Convergence Zone (ITCZ)and South Pacific Convergence Zone (SPCZ). The inter-model differences show that the precipitation change in response to seasurface temperature (SST) change increases with enhanced precipitation climatology. The ENSO-related precipitation anomalycan also be related to the spatial distribution of the mean-state precipitation. The simulated ITCZ/SPCZ displacements aresignificantly correlated with the spatial precipitation-SST relationship in the mean state. Models with stronger mean-stateprecipitation-SST relationship also produce stronger SPCZ/ITCZ displacements. The mean-state tropical precipitation also hasstrong impacts on the ENSO-related precipitation anomalies over East Asia. In the BCC PPE, the connection between the mean-state precipitation and ENSO-related precipitation anomaly is overall consistent with that in AMIP5. Parameters associated withthe low-cloud and deep convection processes are the most influential ones for the precipitation simulations in BCC. Comparedwith the version in AMIP5, the new BCC model can better simulate the precipitation climatology and the relationship betweenENSO and precipitation over southern China. These results have important implications for model development and model-basedclimate predictions.

Yang, Ben↗
DOE OSTI2022

Inference of bipolar neutrino flavor oscillations near a core-collapse supernova based on multiple measurements at Earth

Neutrinos in compact-object environments, such as core-collapse supernovae, can experience various kinds of collective effects in flavor space, engendered by neutrino-neutrino interactions. These include “bipolar” collective oscillations, which are exhibited by neutrino ensembles where different flavors dominate at different energies. Considering the importance of neutrinos in the dynamics and nucleosynthesis in these environments, it is desirable to ascertain whether an Earth-based detection could contain signatures of bipolar oscillations that occurred within a supernova envelope. To that end, we, in this study, continue examining a cost-function formulation of statistical data assimilation (SDA) to infer solutions to a small-scale model of neutrino flavor transformation. SDA is an inference paradigm designed to optimize a model with sparse data. Our model consists of two monoenergetic neutrino beams with different energies emanating from a source and coherently interacting with each other and with a matter background, with radially varying interaction strengths. We attempt to infer flavor transformation histories of these beams using simulated measurements of the flavor content at locations “in vacuum” (that is, far from the source), which could in principle correspond to Earth-based detectors. Within the scope of this small-scale model, we found that: (i) based on such measurements, the SDA procedure is able to infer whether bipolar oscillations had occurred within the protoneutron star envelope, and (ii) if the measurements sample the full amplitude of the neutrino oscillations in vacuum, then the amplitude of the prior bipolar oscillations is well predicted. This result intimates that the inference paradigm can well complement numerical integration codes, via its ability to infer flavor evolution at physically inaccessible locations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1997

Theories of Simplification and Scaling of Spatially Distributed Processes

The problem of scaling is at the heart of ecological theory, the essence of understanding and of the development of a predictive capability. The description of any system depends on the spatial, temporal, and organizational perspective chosen; hence it is essential to understand not only how patterns and dynamics vary with scale, but also how patterns at one scale are manifestations of processes operating at other scales. Evolution has shaped the characteristics of species in ways that result in scale displacement: Each species experiences the environment at its own unique set of spatial and temporal scales and interfaces the biota through unique assemblages of phenotypes. In this way, coexistence becomes possible, and biodiversity is enhanced. By averaging over space, time, and biological interactions, a genotype filters variation at fine scales and selects the arena in which it will face the vicissitudes of nature. Variation at finer scales is then noise, of minor importance to the survival and dynamics of the species, and consequently of minor importance in any attempt at description. In attempting to model ecological interactions in space, contributors throughout this book have struggled with a trade-off between simplification and "realistic" complexity and detail. Although the challenge of simplification is widely recognized in ecology, less appreciated is the intertwining of scaling questions and scaling laws with the process of simplification. In the context of this chapter simplification will in general mean the use of spatial or ensemble means and low-order moments to capture more detailed interactions by integrating over given areas. In this way, one can derive descriptions of the system at different spatial scales, which provides the essentials for the extraction of scaling laws by examination of how system properties vary with scale.

Levin, Simon A.↗
NASA NTRS2012

Impact of Satellite Atmospheric Motion Vectors In the GMAO GEOS-5 Global Data Assimilation System

The WMO and THORPEX co-sponsored fifth Workshop on the Impact of Various Observing Systems on Numerical Weather Prediction will be organized by the Expert Team on the Evolution of the Global Observing System in Sedona, Arizona, United States, from 22 to 25 May 2012. Participants are expected to come from all the major NWP centres which are active in the area of impact studies. The workshop will be conducted in English. As for the first four workshops it is planned to produce a workshop report to be published as a WMO Technical Report that will include the papers submitted by the participants. The previous four workshops in this series took place in Geneva {April 1997), Toulouse (March 2000), Alpbach (March 2004) and Geneva (May 2008). Results from Observing System Experiments (OSEs), both with global and regional aspects were presented and conclusions were drawn concerning the contributions of the various components of the observing system to the large scale forecast skill at short and medium range (Workshop Proceedings were published as WMO World Weather Watch Technical Reports TD No. 868, 1034, 1228 and 1450). Since then, some significant changes and developments have affected the global observing system and more efforts have been devoted to meso-scale observing and assimilation systems. There has also been a trend toward using techniques other than OSEs to document data impact, such as adjoint-based sensitivity to observations or ensemble-based sensitivity. Field experiments have been carried out, in particular through the THORPEX project, and the use of targeted data has been assessed.

Gelaro, Ronald↗
NASA NTRS2024

Dust Model Sensitivity to Dust Source Mask, Sandblasting Efficiency, Air Density, and Land Use: Implications for Model Improvement

This study compares dust storm simulations using two commonly adopted methods for representing four important dust emission parameters. Compared to a dynamic dust source mask based on land use and vegetation cover, a static mask based solely on land use overestimates dust concentration and optical depth by a factor of 2, besides generating spurious emissions. The results reinforce that seasonal variations in vegetation cover can significantly affect dust emissions. For sandblasting efficiency, a clay-dependent semiempirical expression produces 10 times more dust than a physics-based expression. Simulations using model-predicted versus a fixed constant for air density differ by only 8%. However, this difference could range between 12 and 22% for annual simulations over global dust source regions. Simulations with updated versus old land use data, using the same dust source mask, differ twofold, indicating the significant impact of land use change on regional dust emission in central Arizona. The differences in the pairs of these simulations are generally larger than the uncertainty due to meteorology. The simulations align better with observation when using the dynamic dust source mask, the physics-based sandblasting efficiency, and the up-to-date land use data. Given the high sensitivity of dust to surface conditions, the results discussed have implications for improving the dust cycle in weather and climate models and for interpreting model intercomparisons.

Dust-emission source↗
NASA NTRSDate not supplied

Reconstructing Subgrid Cloud Variability Guided by CloudSat/CALIPSO Observations

Predicting realistic cloud subgrid variability remains a challenge for Global Climate Models (GCMs) even though it can play an important role for proper representation of processes pertaining to cloud microphysics, precipitation, and radiation, but also for comparisons with satellite observations which are of much higher resolution than model grids. The diagnostic approach of subgrid cloud variability is commonly handled by subcolumn cloud generators. For a specific GCM, one ideally wants the subgrid variability used for comparisons with satellite observations to be created by the same generator and with the same rules as the one used for model integration. With this in mind, we have embarked in an effort to test and improve cloud subcolumn generators appropriate for GCMs. For this purpose, we use cloud (hydrometeor) products from active observations by the CloudSat radar (CPR) and the CALIPSO lidar (CALIOP). Cloud products from active sensors while suffering significant sampling and coverage drawbacks have the advantage of resolving both horizontal and vertical variability. The main question is: given a profile of cloud condensate mean and variance, can we create a subgrid cloud field that is statistically similar to the observed subgrid cloud field? By “statistically”, we suggest that we do not aspire to reproduce “well” each individual subgrid cloud field of a GCM-scale region, but that our generator performs well for a large ensemble of cases. We simulate radar, passive imager, and radiation flux fields from the observed 2D cloud fields and create one-point statistics; we use the profiles of cloud condensate mean and variance as input to the generator to create subgrid cloud fields; we compare the statistics; we adjust the rules of the generator to create the best possible agreement between the radar, imager and radiation field statistics. In this process, the active observations have actually a dual role: they provide the actual subgrid cloud field which can be reduced to a profile of mean and variance used by the subcolumn generator, but they also provide the rules needed by the generator, such as measures (e.g, decorrelation length) of the vertical overlap of cloud fraction and of the condensate horizontal variability. This presentation will show our progress using Cloudsat and CALIPSO products in this dual fashion.

cloud↗
NASA NTRSDate not supplied

Reconstructing Subgrid Cloud Variability Guided by CloudSat/CALIPSO Observations

Predicting realistic cloud subgrid variability remains a challenge for Global Climate Models (GCMs) even though it can play an important role for proper representation of processes pertaining to cloud microphysics, precipitation, and radiation, but also for comparisons with satellite observations which are of much higher resolution than model grids. The diagnostic approach of subgrid cloud variability is commonly handled by subcolumn cloud generators. For a specific GCM, one ideally wants the subgrid variability used for comparisons with satellite observations to be created by the same generator and with the same rules as the one used for model integration. With this in mind, we have embarked in an effort to test and improve cloud subcolumn generators appropriate for GCMs. For this purpose, we use cloud (hydrometeor) products from active observations by the CloudSat radar (CPR) and the CALIPSO lidar (CALIOP). Cloud products from active sensors while suffering signi􀂦cant sampling and coverage drawbacks have the advantage of resolving both horizontal and vertical variability. The main question is: given a profile of cloud condensate mean and variance, can we create a subgrid cloud field that is statistically similar to the observed subgrid cloud field? By “statistically”, we suggest that we do not aspire to reproduce “well” each individual subgrid cloud field of a GCM-scale region, but that our generator performs well for a large ensemble of cases. We simulate radar, passive imager, and radiation flux fields from the observed 2D cloud fields and create one-point statistics; we use the profiles of cloud condensate mean and variance as input to the generator to create subgrid cloud fields; we compare the statistics; we adjust the rules of the generator to create the best possible agreement between the radar, imager and radiation field statistics. In this process, the active observations have actually a dual role: they provide the actual subgrid cloud field which can be reduced to a profile of mean and variance used by the subcolumn generator, but they also provide the rules needed by the generator, such as measures (e.g, decorrelation length) of the vertical overlap of cloud fraction and of the condensate horizontal variability. This presentation will show our progress using Cloudsat and CALIPSO products in this dual fashion.

clouds↗
NASA NTRS2012

Projections of Declining Surface-Water Availability for the Southwestern United States

16 of the CMIP5 models had all the data needed for this work for at least one simulation that was continuous from 1950 to 2040. Details of the models analyzed here are provided in Table S1. The model data analyzed here are available at http://strega.ldeo.columbia.edu:81/expert/home/.naomi/.AR5/.v2/.historical:rcp85/.mmm16/ a. Assessing the climatology of the models Despite increases in horizontal resolution of many models compared to their CMIP3 counterparts none of these models can adequately resolve the topography of the south west United States, such as the Sierra Nevada and Rocky Mountains and the associated orographic precipitation. This requires that caution be used when interpreting the results presented here. To assess the ability of the models to simulate the current hydroclimate, in Figure S1 we show the observed (from the Global Precipitation Climatology Centre gridded rain gauge data, (1)) monthly climatology of precipitation and the same for all the models and the multimodel mean for the California-Nevada, Colorado headwaters and Texas regions. The GPCC data uses rain gauges only and interpolates to regular grids of which we used the 1◦ by 1◦ one. Details of the data set can be found in (2). While the models apparently overestimate precipitation in California and Nevada the seasonal cycle with wet winters and dry summers is very well represented. It is also possible that the rain gauge observations are biased low by inadequately sampling the higher mountain regions. How ever the models might also be expected to underestimate orographic precipitation due to inadequate horizontal resolution. The 25 models are also too wet in the Colorado headwaters region but correctly represent the quite even distribution though the year. The bimodal distribution of precipitation in Texas, with peaks in May and September, and the absolute amounts, are well modeled but with the September peak too weak. The positive precipitation bias translates into a positive runoff bias for the Colorado headwaters as also shown in Figure S1. Here the observed runoff values are taken from simulations of the Variable Infiltration Capacity (VIC) land surface-hydrology model (3) forced by observed meteorology (5) that were conducted as part of the North American Land Data Assimilation System project phase 2 ( (NLDAS-2), http://www.emc.ncep.noaa.gov/mmb/nldas/. Runoff for California-Nevada is better simulated but there is a positive bias over Texas despite no strong precipitation bias. To check whether regional climate models better simulate P and runoff in these regions we analyzed the historical simulation with the Regional Climate Model version 3 driven by the National Centers for Environmental Prediction-Department of Energy Reanalysis 2 available from the North American Regional Climate Change Assessment Program (http://www.narccap.ucar.edu). This model configuration retained these biases in P and runoff although they were reduced in amplitude. Given these varying biases we plot P and P − E changes in actual values but apply the simplest bias correction possible to the runoff and soil moisture values and show the modeled changes in terms of percentages of the 20th Century model climatologies. A thorough assessment of the simulation of North American climate in CMIP5 models is conducted in Sheffield at al. (North American Climate in CMIP5 Experiments. Part I: Evaluation of 20th Century Continental and Regional Climatology, manuscript submit ted to J. Climate, available at http://www.climate.noaa.gov/index.jsp?pg=./cpo pa/ mapp/cmip5 publications.html). Sheffield et al. analyze the climatology of precipitation, surface air temperature, low level winds, moisture fluxes, runoff etc. and conclude that the main features of the hydrological cycle, including characteristics of the atmospheric moisture balance and its seasonality, are captured in the CMP5 models subject to biases in total precipitation amounts. We chose to use all available models instead of selecting some and rejecting others based on an assessment of model realism. This is in accord with the suggestions of Mote et al. for CMIP3 (4) but future work needs to revisit this matter for the case of the CMIP5 ensemble.

Seager, Richard↗
DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Attribution of Chemistry-Climate Model Initiative (CCMI) Ozone Radiative Flux Bias from Satellites

The top-of-atmosphere (TOA) outgoing longwave flux over the 9.6-μm ozone band is a fundamental quantity for understanding chemistry-climate coupling. However, observed TOA fluxes are hard to estimate as they exhibit considerable variability in space and time that depend on the distributions of clouds, ozone (O3), water vapor (H2O), air temperature (Ta), and surface temperature (Ts). Benchmarking present day fluxes and quantifying the relative influence of their drivers is the first step for estimating climate feedbacks from ozone radiative forcing and predicting radiative forcing evolution. To that end, we constructed observational instantaneous radiative kernels (IRKs) under clear-sky conditions, representing the sensitivities of the TOA flux in the 9.6-μm ozone band to the vertical distribution of geophysical variables, including O3, H2O, Ta, and Ts based upon the Aura Tropospheric Emission Spectrometer (TES) measurements. Applying these kernels to present-day simulations from the Chemistry-Climate Model Initiative (CCMI) project as compared to a 2006 reanalysis assimilating satellite observations, we show that the models have large differences in TOA flux, attributable to different geophysical variables. In particular, model simulations continue to diverge from observations in the tropics, as reported in previous studies of the Atmospheric Chemistry Climate Model Inter-comparison Project (ACCMIP) simulations. The principal culprits are tropical mid and upper tropospheric ozone followed by tropical lower tropospheric H2O. Five models out of the eight studied here have TOA flux biases exceeding 100 mWm-2 attributable to tropospheric ozone bias. Another set of five models have flux biases over 50 mWm-2 due to H2O. On the other hand, Ta radiative bias is negligible in all models (no more than 30 mWm-2). We found that AM3 and CMAM have the lowest TOA flux biases globally but are a result of cancellation of opposite biases due to difference processes. Overall, the multi-model ensemble mean bias is –133±98 mWm-2, indicating that they are too atmospherically opaque due to trapping too much radiation in the atmosphere by overestimated tropical tropospheric O3 and H2O. Having too much O3 and H2O in the troposphere would have different impacts on the sensitivity of TOA flux to O3 and these competing effects add more uncertainties on the ozone radiative forcing. We find that the inter-model TOA outgoing longwave radiation (OLR) difference is well anti-correlated with their ozone band flux bias. This suggests that there is significant radiative compensation in the calculation of model outgoing longwave radiation.

Aura Tropospheric Emission Spectrometer (TES) meas↗
DOE OSTI2026

How Well Can CMIP6 Models Represent the Observed Influence of the Pacific and Indian Oceans on the Indian Summer Monsoon Rainfall?

This study evaluates the ability of CMIP6 climate models to simulate the observed effects of tropical Pacific and Indian Ocean sea surface temperature anomalies (SSTAs) on Indian summer monsoon rainfall (ISMR) variability. Using observational data and the large ensemble historical simulations of seven CMIP6 models from 1950 to 2014, we applied a cyclostationary linear inverse model (CS-LIM) to isolate the impacts of tropical Pacific SSTAs, Indian Ocean SSTAs and their interaction on the interannual variability of ISMR. Overall, CMIP6 models well reproduced the observed enhanced (reduced) ISMR variability from Pacific SSTAs (Indian Ocean SSTAs and the Indo-Pacific interaction), but with varying spatial patterns and magnitudes. While CESM2 and E3SM-2-0 showed the best agreement with observations for the effects of Pacific SSTAs and the Indo-Pacific interaction, respectively, CMIP6 models showed mixed results for the impacts from Indian Ocean SSTAs. Composite analysis of ISMR anomalies during the developing phases of pure and co-occurring El Niño-Southern Oscillation (ENSO) and Indian Ocean dipole (IOD) events revealed that the impacts from Pacific SSTAs were captured reasonably well by E3SM-2-0, CESM2, MIROC6, and MPI-ESM1-2-LR, while E3SM-2-0 also showed the best agreement with observations for the effects from the Indo-Pacific interaction. However, all models showed substantial biases in simulating the Indian Ocean SSTA impacts on ISMR, especially for pure El Niño events. Overall, this study provides new insights into how individual CMIP6 models simulate the isolated impacts from the tropical Pacific and Indian Oceans, which has important applications for improving ISMR predictions and interpreting ISMR future projections.

monsoon↗