SEARCH · Search NASA

Results for “complex ensemble”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Modeling aerosol transmission spectra from n(λ) and k(λ) infrared optical constants measurements of organic liquids and solids

The effects of light scattering and refraction play significantly different roles for aerosols than for bulk materials, making it challenging to identify aerosolized chemicals using traditional spectral methods or spectral reference libraries. Due to a potentially infinite number of particle morphologies, sizes, and compositions, constructing a database of laboratory-measured aerosol spectra is not a practical solution. Here, as an alternative approach, the measured n / k optical vectors of two example organic materials (diethyl phthalate and D-mannitol) are used in combination with particle absorption / scattering theory (Mie theory and FDTD) and the Beer-Lambert law to generate a series of synthetic infrared transmission / scattered light spectra. The synthetic spectra show significant differences versus simple slab transmission spectra, even for small changes in particle size (e.g., 5 vs. 10 µm) for both single particles and ensembles, potentially serving as useful reference data for aerosol sensing. For spherical single particles with diameters of 1 to 10 µm, FDTD simulations predict changes in the magnitudes of spectral shifts and the shapes of the peaks vs. particle size with only small deviations from Mie theory predictions, yet reliably capture the direction of the shifts. Typical spectral peak shifts in the longwave infrared correspond to Δλ ∼0.20 µm (∼34 cm -1 ) when compared to corresponding slab transmission spectra. Additionally, synthetic spectra generated from the n / k values derived using two different methods (KBr pellet transmission and single-angle reflectance) are compared using the Mie theory model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Attribution of the Australian bushfire risk to anthropogenic climate change

Abstract. Disastrous bushfires during the last months of 2019 and January 2020 affected Australia, raising the question to what extent the risk of these fires was exacerbated by anthropogenic climate change. To answer the question for southeastern Australia, where fires were particularly severe, affecting people and ecosystems, we use a physically based index of fire weather, the Fire Weather Index; long-term observations of heat and drought; and 11 large ensembles of state-of-the-art climate models. We find large trends in the Fire Weather Index in the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) Atmospheric Reanalysis (ERA5) since 1979 and a smaller but significant increase by at least 30 % in the models. Therefore, we find that climate change has induced a higher weather-induced risk of such an extreme fire season. This trend is mainly driven by the increase of temperature extremes. In agreement with previous analyses we find that heat extremes have become more likely by at least a factor of 2 due to the long-term warming trend. However, current climate models overestimate variability and tend to underestimate the long-term trend in these extremes, so the true change in the likelihood of extreme heat could be larger, suggesting that the attribution of the increased fire weather risk is a conservative estimate. We do not find an attributable trend in either extreme annual drought or the driest month of the fire season, September–February. The observations, however, show a weak drying trend in the annual mean. For the 2019/20 season more than half of the July–December drought was driven by record excursions of the Indian Ocean Dipole and Southern Annular Mode, factors which are included in the analysis here. The study reveals the complexity of the 2019/20 bushfire event, with some but not all drivers showing an imprint of anthropogenic climate change. Finally, the study concludes with a qualitative review of various vulnerability and exposure factors that each play a role, along with the hazard in increasing or decreasing the overall impact of the bushfires.

van Oldenborgh, Geert Jan (ORCID:0000000268989535)↗

NASA NTRSDate not supplied

Satellite-Assisted Particulate Matter (SAPM) for the Models, In situ, and Remote Sensing of Aerosols (MIRA) Working Group

Models, In situ, and Remote sensing of Aerosols (MIRA) is an international working group that provides a forum for collaborations amongst these three atmospheric aerosol communities. MIRA consists of a collection of interdisciplinary and independently funded Topic/Project groups with clear goals and generally characterized by requests for additional scientific data. The Satellite-Assisted Particulate Matter (SAPM) Topic group focuses on the study of fine particulate matter (PM2.5), as it is a major contributor to air pollution and negatively impacts human health. Our group aims to provide intercomparisons of various methods and techniques for estimating surface PM2.5 assisted by satellite remote sensors (passive and active), global aerosol models, and in situ aerosol measurements. The overall motivation of our work is to enhance PM2.5 coverage beyond in situ ground stations, which can be limited. In this presentation, we provide an overview of various PM2.5 estimation approaches by current SAPM team members, as each approach has its own strengths and limitations. One approach uses near-surface aerosol extinction retrievals from the spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and an assumed value of mass extinction efficiency to derive PM2.5 concentrations (Toth et al 2019). A study using this method found a promising agreement (R=0.60; slope=0.89) with U.S. Environmental Protection Agency (EPA) in situ PM2.5 measurements over the contiguous United States (CONUS) for a 12-year (2007-2018) period (Toth et al 2022). In a different approach combining spaceborne lidar and a model, the Cloud Aerosol Transport System (CATS) lidar and Goddard Earth Observing System (GEOS) model are blended in an ensemble variational assimilation scheme to retrieve aerosol extinction and convert speciated mass concentrations to total PM2.5. A recent study (Matus et al in prep.) found agreement within 2 μg/m3 over the CONUS on the 2016 annual mean when compared to EPA ground-based measurements. Another study explored trends in city aerosol optical depth (AOD) from a spaceborne passive remote sensor, Moderate Resolution Imaging Spectroradiometer (MODIS), and found that these trends agree well with trends in city surface PM2.5, which provides confidence in the use of city AOD for city PM2.5 trend studies (Vohra et al 2021). In a follow-on study, steep and significant trends (~2.5 to 7.8% a-1) in AOD are found for several South Asian cities (e.g., Bangalore and Hyderabad), suggesting rapid growth in PM2.5 (Vohra et al 2022). Also in this presentation, we report results on our current collective SAPM study of India, a country that exhibits high levels of PM2.5 pollution. We show the temporal (i.e., seasonal) and spatial variability of PM2.5 over India using the CALIOP and CATS+model approaches, comparisons between the estimated PM2.5 from the retrieval schemes, and validation with in situ PM2.5 ground-based observations from the Central Pollution Control Board in India. This region provides us with an excellent case study to test our PM2.5 retrievals in heavily polluted scenes and over a complex and varying topography. The SAPM Topic group actively seeks international participants/collaborators in our group, including those working with in situ aerosol measurements (e.g., ACTRIS). We are interested in aerosol datasets in order to either improve our PM2.5 estimates from various approaches (e.g., using in situ mass scattering/absorption coefficient and aerosol hygroscopic properties for various aerosol species) and/or validate the PM2.5 estimates (e.g., using in situ ground-based PM2.5 concentrations).

Travis D Toth↗

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

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↗

DOE OSTI2024

Chemical properties and single-particle mixing state of soot aerosol in Houston during the TRACER campaign

Abstract. A high-resolution soot particle aerosol mass spectrometer (SP-AMS) was used to selectively measure refractory black carbon (rBC) and its associated coating material using both the ensemble size-resolved mass spectral mode and the event trigger single particle (ETSP) mode in Houston, Texas, in summer 2022. This study was conducted as part of the Department of Energy Atmospheric Radiation Measurement (ARM) program's TRacking Aerosol Convection interactions ExpeRiment (TRACER) field campaign. The study revealed an average (±1σ) rBC concentration of 103 ± 176 ng m−3. Additionally, the coatings on the BC particles were primarily composed of organics (59 %; 219 ± 260 ng m−3) and sulfate (26 %; 94 ± 55 ng m−3). Positive matrix factorization (PMF) analysis of the ensemble mass spectra of BC-containing particles resolved four distinct types of soot aerosol, including an oxidized organic aerosol (OOABC,PMF) factor associated with processed primary organic aerosol, an inorganic sulfate factor (SO4,BC,PMF), an oxidized rBC factor (O-BCPMF), and a mixed mineral dust–biomass burning aerosol factor with significant contribution from potassium (K-BBBC,PMF). Additionally, K-means clustering analysis of the single-particle mass spectra identified eight different clusters, including soot particles enriched in hydrocarbon-like organic aerosol (HOABC,ETSP), sulfate (SO4,BC,ETSP), two types of rBC, OOA (OOABC,ETSP), chloride (ClBC,ETSP), and nitrate (NO3,BC,ETSP). The single-particle measurements demonstrate substantial variation in BC coating thickness with coating-to-rBC mass ratios ranging from 0.1 to 100. The mixing state index (χ), which denotes the degree of homogeneity of the soot aerosol, varied from 4 % to 94 % with a median of 40 %, indicating that the aerosol population lies in between internal and external mixing but has large temporal and source type variability. In addition, a significant fraction of BC-containing particles, a majority enriched with oxidized organics and sulfate, exhibit sufficiently high κ values and diameters conducive to activation as cloud nuclei under atmospherically relevant supersaturation conditions. This finding bears significance in comprehending the aging processes of rBC-containing particles and their activation into cloud droplets. Our analysis highlights the complex nature of soot aerosol and underscores the need to comprehend its variability across different environments for accurate assessment of climate change.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2001

Quasi-Equilibrium States in the Tropics Simulated by a Cloud-Resolving Model: Specific Features and Budget Analysis - Part 1

A series of long-term integrations using the two-dimensional Goddard Cumulus Ensemble (GCE) model were performed by altering imposed environmental components to produce various quasi-equilibrium thermodynamic states. Model results show that the genesis of a warm/wet quasi-equilibrium state is mainly due to either strong vertical wind shear (from nudging) or large surface fluxes (from strong surface winds), while a cold/dry quasi-equilibrium state is attributed to a remarkably weakened mixed-wind shear (from vertical mixing due to deep convection) along with weak surface winds. In general, latent heat flux and net large-scale temperature forcing, the two dominant physical processes, dominate in the beginning stage of the simulated convective systems, then considerably weaken in the final stage, which leads to quasi-equilibrium states. A higher thermodynamic regime is found to produce a larger rainfall amount, as convective clouds are the leading source of rainfall over stratiform clouds even though the former occupy much less area. Moreover, convective clouds are more likely to occur in the presence of strong surface winds (latent heat flux), while stratiform clouds (especially the well-organized type) are favored in conditions with strong wind shear (large-scale forcing). The convective systems, which consist of distinct cloud types due to the variation in horizontal winds, are also found to propagate differently. Accordingly, convective systems with mixed-wind shear generally propagate in the direction of shear, while the system with strong (multidirectional) wind shear propagates in a more complex way. Based on the results from the temperature (Q1) and moisture (Q2) budgets, cloud-scale eddies are found to act as a hydrodynamic 'vehicle' that cascades the heat and moisture vertically. Several other specific features such as atmospheric stability, CAPE, and mass fluxes are also investigated and found to be significantly different between diverse quasi-equilibrium states. Detailed comparisons between the various states are presented.

Shie, C.-L.↗

DOE OSTI2020

The shared socio-economic pathway (SSP) greenhouse gas concentrations and their extensions to 2500

Anthropogenic increases in atmospheric greenhouse gas concentrations are the main driver of current and future climate change. The integrated assessment community has quantified anthropogenic emissions for the shared socio-economic pathway (SSP) scenarios, each of which represents a different future socio-economic projection and political environment. Here, we provide the greenhouse gas concentrations for these SSP scenarios using the reduced-complexity climate–carbon-cycle model MAGICC7.0. We extend historical, observationally based concentration data with SSP concentration projections from 2015 to 2500 for 43 greenhouse gases with monthly and latitudinal resolution. CO 2 concentrations by 2100 range from 393 to 1135 ppm for the lowest (SSP1-1.9) and highest (SSP5-8.5) emission scenarios, respectively. We also provide the concentration extensions beyond 2100 based on assumptions regarding the trajectories of fossil fuels and land use change emissions, net negative emissions, and the fraction of non-CO 2 emissions. By 2150, CO 2 concentrations in the lowest emission scenario are approximately 350 ppm and approximately plateau at that level until 2500, whereas the highest fossil-fuel-driven scenario projects CO 2 concentrations of 1737 ppm and reaches concentrations beyond 2000 ppm by 2250. We estimate that the share of CO 2 in the total radiative forcing contribution of all considered 43 long-lived greenhouse gases increases from 66% for the present day to roughly 68% to 85% by the time of maximum forcing in the 21st century. For this estimation, we updated simple radiative forcing parameterizations that reflect the Oslo Line-By-Line model results. In comparison to the representative concentration pathways (RCPs), the five main SSPs (SSP1-1.9, SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5) are more evenly spaced and extend to lower 2100 radiative forcing and temperatures. Performing two pairs of six-member historical ensembles with CESM1.2.2, we estimate the effect on surface air temperatures of applying latitudinally and seasonally resolved GHG concentrations. We find that the ensemble differences in the March–April–May (MAM) season provide a regional warming in higher northern latitudes of up to 0.4 K over the historical period, latitudinally averaged of about 0.1 K, which we estimate to be comparable to the upper bound (~5 % level) of natural variability. In comparison to the comparatively straight line of the last 2000 years, the greenhouse gas concentrations since the onset of the industrial period and this studies' projections over the next 100 to 500 years unequivocally depict a “hockey-stick” upwards shape. The SSP concentration time series derived in this study provide a harmonized set of input assumptions for long-term climate science analysis; they also provide an indication of the wide set of futures that societal developments and policy implementations can lead to – ranging from multiple degrees of future warming on the one side to approximately 1.5°C warming on the other.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

The Influence of Anomalous Biomass Emissions on ENSO in CESM2

Abstract The influence of biomass burning (BB) aerosols arising from wildfires and agricultural fires on the transient coupled evolution of El Niño–Southern Oscillation (ENSO) is explored in Community Earth System Model, version 2 (CESM2). For both El Niño and La Niña, two 20-member ensembles are generated from initial states that are predisposed to evolve into ENSO events. For each ENSO phase, one ensemble is forced with the observed BB emissions during satellite-era ENSO events while the other is forced with a climatological annual cycle, with the responses to anomalous BB emissions estimated from interensemble differences. It is found that the regional responses to anomalous BB emissions occur mainly during boreal fall, which is also the time of the climatological seasonal maximum in emissions. Transient responses are identified in precipitation, clouds, and radiation in both the tropics and extratropics. At the onset of El Niño, these include increased precipitation in the northern branch of the intertropical convergence zone (ITCZ) and an enhancement of cloud albedo and amount across the Maritime Continent and eastern subtropical Pacific Ocean. Additional responses are identified through the course of El Niño and successive La Niña events, the net effect of which is to strengthen sea surface temperature (SST) anomalies in the eastern Pacific Ocean during El Niño and warm the tropical Pacific Ocean during La Niña. These responses improve the simulation of ENSO power, diversity, and asymmetry in CESM2. Significance Statement Biomass burning emissions from both wildfires and agricultural fires during ENSO act as a climate feedback, as they are both driven by fire weather anomalies tied to ENSO and influence clouds, radiation, and precipitation on a global scale. While the largest biomass burning (BB) emissions emanate mainly from Indonesia, other regions also exhibit systematic responses. When simulated in CESM2, the emissions collectively impact the flows of energy across the tropical Pacific Ocean in a spatially and temporally complex fashion, with the net effect of improving CESM2’s simulation of ENSO events.

Fasullo, John T. [National Center for Atmospheric ↗

NASA NTRS2008

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗

NASA NTRS2016

Polarized Radiative Transfer of a Cirrus Cloud Consisting of Randomly Oriented Hexagonal Ice Crystals: The 3 x 3 Approximation for Non-Spherical Particles

The reflection and transmission of polarized light for a cirrus cloud consisting of randomly oriented hexagonal columns were calculated by two very different vector radiative transfer models. The forward peak of the phase function for the ensemble-averaged ice crystals has a value of order 6 x 10(exp 3) so a truncation procedure was used to help produce numerically efficient yet accurate results. One of these models, the Vectorized Line-by-Line Equivalent model (VLBLE), is based on the doubling- adding principle, while the other is based on a vector discrete ordinates method (VDISORT). A comparison shows that the two models provide very close although not entirely identical results, which can be explained by differences in treatment of single scattering and the representation of the scattering phase matrix. The relative differences in the reflected I and Q Stokes parameters are within 0.5 for I and within 1.5 for Q for all viewing angles. In 1971 Hansen showed that for scattering by spherical particles the 3 x 3 approximation is sufficient to produce accurate results for the reflected radiance I and the degree of polarization (DOP), and he conjectured that these results would hold also for non-spherical particles. Simulations were conducted to test Hansen's conjecture for the cirrus cloud particles considered in this study. It was found that the 3 x 3 approximation also gives accurate results for the transmitted light, and for Q and U in addition to I and DOP. For these non-spherical ice particles the 3 x 3 approximation leads to an absolute error 2 x 10(exp -6) for the reflected and transmitted I, Q and U Stokes parameters. Hence, it appears to be an excellent approximation, which significantly reduces the computational complexity and burden required for multiple scattering calculations.

Polarization↗

NASA NTRSDate not supplied

Medium-Range River Flood Forecasts Using a Long Short-Term Memory Network

River flooding and the impacts are a concern for decision makers throughout the United States. Accurate medium-range forecasts (~3-7 days) are critical for providing advanced outlooks to emergency management officials. Unfortunately, accurately forecasting rainfall-runoff and the subsequent rise and fall within rivers remain a challenge in hydrological modeling. While complex physical modeling systems are the standard for representing the hydrological processes, they are computationally demanding and can require extensive calibration. Further, uncertainties remain in the model parameters and input data. The use of machine learning can reduce some of the computational demand while maintaining high accuracy. Therefore, this project makes use of a Long Short-Term Memory (LSTM) network which explicitly accounts for the time-dependent nature of rainfall-runoff modeling. The developed LSTM was trained to predict river gauge height, or stage height, based on time-lagged input features which include: gauge height to initialize the model, the NASA Short-term Prediction Research and Transition Center’s instance of the Land Information System (SPoRT-LIS) relative soil moisture to describe the rainfall infiltration rate, and 6-hr Multi-Radar Multi-Sensor quantitative precipitation estimate (MRMS QPE). The developed LSTM based system is then used to produce 7-day forecasts with a 6-hr temporal resolution using three different quantitative precipitation forecasts (QPF) from the NWS’s Weather Prediction Center (WPC), the NCEP Global Forecast System (GFS) model and the National Blend of Models (NBM). This trained modeling system has been implemented as an experimental product at over 100 different rivers in collaboration with at multiple National Weather Service (NWS) Forecast Offices and River Forecast Centers (RFC) across the eastern half of the United States. The developed LSTM model achieved average Nash-Sutcliffe efficiency (NSE) 0.89 higher than the equivalent medium-range National Water Model ensemble member forecast over a 7-day forecast. In addition to the initial development and evaluation, this project has continued to expand. While the initial model was developed for precipitation dominated basins, expansion of the project has taken it to basins effected by snow melt. This presentation will provide an overview of the project with focus on recent developments on incorporating snow melt processes into the model.

Andrew T. White↗

DOE OSTI2025

Conformational heterogeneity in the dGsw purine riboswitch: role of Mg²⁺ and 2’-dG in aptamer folding

Recent advancements in RNA structural biology have focused on unraveling the complexities of non-coding mRNA elements like riboswitches. These cis-acting regulatory regions undergo structural changes in response to specific cellular metabolites, leading to up or downregulation of downstream genes. The purine riboswitch family regulates many prokaryotic genes involved in purine degradation and biosynthesis. They feature an aptamer domain organized around a 3-way helical junction, where ligand encapsulation occurs at the junctional core. In our study, we chemically probed the aptamer domain of the 2’-dG-sensing purine riboswitch from Mesoplasma florum (dGsw) under various solution conditions to understand how Mg²⁺ and 2’-dG influence riboswitch folding. Here, we find that efficient 2’-dG binding strongly depends on Mg²⁺, indicating that Mg²⁺ is essential for priming dGsw for ligand interactions. We identified a previously undescribed sequence in the 5’ tail of dGsw that is complementary to a conserved helix. The inclusion of this region in a construct led to intramolecular competition between the alternate helix, Palt, and P1. Mutational analysis confirmed that 5’ flanking end of the aptamer domain forms an alternate helix in the absence of ligand. Molecular dynamics simulations revealed that this alternative conformation is stable. This helix may, therefore, facilitate the formation of an anti-terminator helix by opening the 3-way junction surrounding the 2’-dG binding site. Our study further establishes the importance of a closed terminal P1 helix conformation for metabolite binding and suggests that the delicate interplay between P1 and Palt may fine-tune downstream gene regulation. These insights offer a new perspective on riboswitch structure and enhance our understanding of the role that a conformational ensemble plays in riboswitch activity and regulation.

Biochemistry & Molecular Biology↗

DOE OSTI2020

Mechanistic understanding of catalytic conversion of ethanol to 1-butene over 2D-pillared MFI zeolite

Ethanol is an important C2 platform molecule for producing value-added chemicals and distillate hydrocarbon fuels (e.g., jet and diesel). Among these, catalytic upgrading of ethanol to butenes can generate valuable commodity chemicals (e.g., 1-butene) and provide C4 olefin intermediates that can be further upgraded to jet/diesel fuels. 2D zeolites offer hierarchical mesoporous structures, leading to improved mass transport and reduced diffusion length, which can help to address the coking challenges faced by ethanol conversion to hydrocarbons over 3D zeolites. In this study, we investigate the acid-catalyzed conversion of ethanol to 1-butene over the Brønsted acid sites (BAS) in 2D-pillared MFI zeolite (2D-PMFI) using ab initio molecular dynamics (AIMD) simulations, in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and calorimetric measurements. A detailed thermodynamic analysis, using quasi-harmonic approximation (QHA), on the Gibbs free energy pathway of ethanol conversion shows that the consideration of entropy is critical to accurately capture the detailed thermodynamic profiles. The formation of bulky ethoxy-ethene complex is found to be a potential rate-determining step (RDS), proceeding via stepwise mechanism. The reactivity of 2D-PMFI can be further tuned by manipulating RDS through a careful control of number of BAS and operating temperatures. The accuracy of computational model is confirmed through the comparison of calculated vibrational density of states (VDOS) against the experimental DRIFTS measurements. Overall, our study provides mechanistic insights into ethanol upgrading over the 2D-PMFI and shows the importance of evaluating entropic effects in such confined system.

Ethanol conversion, 2D-pillared zeolite, MFI, Anha↗

NASA NTRS2005

Electronic Transport through Self Assembled Thiol Molecules: Effect of Monolayer Order, Dynamics and Temperature

We present the charge transport and tunneling conductance of self assembled organic thiol molecules and discuss the influence of order and dynamics in the monolayer on the transport behavior and the effect of temperature. Conjugated thiol molecular wires and organometals such as terpyridine metal complexes provide a new platform for molecular electronic devices and we study their self assembly on Au(111) substrates by the scanning tunneling microscope. Determining the organization of the molecule and the ability to control the nature of its interface with the substrate is important for reliable performance of the molecular electronic devices. By concurrent scanning tunneling microscopy and spectroscopy studies on SAMs formed from oligo (phenelyne ethynelyne) monolayers with and without molecular order, we show that packing and order determine the response of a self assembled monolayer (SAM) to competing interactions. Molecular resolution STM imaging in vacuum shows that the OPES adopt an imcommensurate SAM structure on Au(111) with a rectangular unit cell. Tunneling spectroscopic measurements were performed on the SAM as a function of junction resistance. STS results show that the I-Vs are non linear and asymmetric due to the inherent asymmetry in the molecular structure, with larger currents at negative sample biases. The asymmetry increases with increasing junction resistance due to the asymmetry in the coupling to the leads. This is brought out clearly in the differential conductance, which also shows a gap at the Fermi level. We also studied the effect of order and dynamics in the monolayer on the charge transport and found that competing forces between the electric field, intermolecular interactions, tip-molecule physisorption and substrate-molecule chemisorption impact the transport measurements and its reliability and that the presence of molecular order is very important for reproducible transport measurements. Thus while developing new electronic platforms based on molecules, it is important to have a good control of the molecule-substrate interface, for the devices to perform reliably. While such a control would minimize fluctuations and dynamics in the ensemble, the real challenge is to develop device architectures that are tolerant to fluctuations, since they cannot be totally eliminated in these low dimensional soft systems. Results of temperature dependent STS measurements will also be discussed.

Dholakia, Geetha↗

DOE OSTI2021

AI-driven multiscale simulations illuminate mechanisms of SARS-CoV-2 spike dynamics

We develop a generalizable AI-driven workflow that leverages heterogeneous HPC resources to explore the time-dependent dynamics of molecular systems. We use this workflow to investigate the mechanisms of infectivity of the SARS-CoV-2 spike protein, the main viral infection machinery. Our workflow enables more efficient investigation of spike dynamics in a variety of complex environments, including within a complete SARS-CoV-2 viral envelope simulation, which contains 305 million atoms and shows strong scaling on ORNL Summit using NAMD. We present several novel scientific discoveries, including the elucidation of the spike’s full glycan shield, the role of spike glycans in modulating the infectivity of the virus, and the characterization of the flexible interactions between the spike and the human ACE2 receptor. We also demonstrate how AI can accelerate conformational sampling across different systems and pave the way for the future application of such methods to additional studies in SARS-CoV-2 and other molecular systems.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Simulation of ENSO Teleconnections to Precipitation Extremes over the United States in the High-Resolution Version of E3SM

Abstract We evaluate the simulated teleconnection of El Niño–Southern Oscillation (ENSO) to winter season precipitation extremes over the United States in a long (98 years) 1950 control high-resolution version (HR; 25-km nominal atmosphere model horizontal resolution) of the U.S. Department of Energy’s (DOE) Energy Exascale Earth System Model version 1 (E3SMv1). The model bias and spatial pattern of ENSO teleconnections to mean and extreme precipitation in HR overall are similar to the low-resolution model’s (LR; 110 km) historical simulation (four-member ensemble, 1925–59). However, over the southeastern United States (SE-U.S.), HR produces stronger El Niño–associated extremes, reducing LR’s model bias. Both LR and HR produce weaker than observed increase in storm track activity during El Niño events there, but HR improves the ENSO-associated variability of moisture transport over SE-U.S. During El Niño, stronger vertical velocities in HR produce stronger large-scale precipitation, causing larger latent heating of the troposphere that pulls in more moisture from the Gulf of Mexico into the SE-U.S. This positive feedback also contributes to the stronger mean and extreme precipitation response in HR. Over the Pacific Northwest, LR’s bias of stronger than observed La Niña associated extremes is amplified in HR. Both models simulate stronger than observed moisture transport from the Pacific Ocean into the region during La Niña years. The amplified HR bias there is due to stronger orographically driven vertical updrafts that create stronger large-scale precipitation, despite weaker La Niña–induced storm track activity. Significance Statement New high-resolution Earth system models (ESMs) solve mathematical equations of fluid flow at much smaller spatial scales than prevalent ESMs, and thus are prohibitively expensive to compute. However, they can be useful for simulating accurate details of regional climate extremes that are driven by naturally occurring climate oscillations like El Niño–Southern Oscillation (ENSO). Here, we evaluate the simulation of ENSO-driven precipitation extremes over the United States in the high-resolution version of the U.S. Department of Energy’s new Energy Exascale Earth System Model version 1. We find that the high-resolution model improves upon its low-resolution counterpart over the southeastern United States by producing a better transport of moisture into the region from the Gulf of Mexico during El Niño. Over the U.S. Pacific Northwest, the high-resolution model simulates the atmospheric flow in more detail over the complex mountainous terrain. However, it also brings in more moisture from the Pacific Ocean just like the low-resolution model. This causes it to produce precipitation extremes during La Niña years there that are stronger than that observed in the real world.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗