Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble integration”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Future changes in extreme precipitation over the San Francisco Bay Area: Dependence on atmospheric river and extratropical cyclone events

Extreme precipitation poses a major challenge for local governments, including the City and County of San Francisco, California, as flooding can damage and destroy infrastructure and property. As the climate continues to warm, reliable future precipitation projections are needed to provide the best possible information to decision makers. However, future changes in the magnitude of extreme precipitation are uncertain, as current state-of-the-art global climate models are typically run at relatively coarse horizontal resolutions that require the use of convective parameterization and have difficulty simulating observed extreme rainfall rates. Here, we performed ensembles of convection-permitting regional climate model simulations to investigate how five historically impactful extreme precipitation events over the San Francisco Bay Area could change if similar events occurred in future climates. We found that changes in storm-total precipitation depend strongly on storm type. Precipitation associated with an atmospheric river (AR) accompanied by an extratropical cyclone (ETC) is projected to increase at a rate exceeding (by up to 1.5 times) the theoretical Clausius Clapeyron scaling of 6–7% per °C warming. On the other hand, future precipitation changes are weak or negative for events characterized by an AR only, despite increases in precipitable water and integrated vapor transport that are similar to those of the co-occurring AR and ETC events. The differences in the sign of future precipitation change between AR-only events and co-occurring AR and ETC events is instead linked with changes in mid-tropospheric vertical velocity. Given that the majority of observed ARs are associated with an ETC, this research has important implications for future precipitation impacts over the Bay Area, as it indicates that storm-total precipitation associated with the most common type of storm event may increase by up to 26–37% in 2100 relative to historical.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Rock Physics-Based Data Assimilation of Integrated Continuous Active-Source Seismic and Pressure Monitoring Data during Geological Carbon Storage

Summary There has been substantial controversy concerning the role of geological carbon storage (GCS) in sequestering anthropogenic carbon emissions to mitigate climate change and global warming. Arguments center on the inability to monitor a geological storage site precisely and continuously, especially highlighting the associated costs and spatiotemporal trade-offs when using conventional subsurface monitoring techniques (well logs, core samples, chemical tracers, and 4D seismics). Active surveillance of GCS sites is essential for managing and mitigating potential leaks but is also required by regulation. With the goal of enhancing the monitoring capability at GCS sites, we present a rock physics-based joint data assimilation model to study a popular GCS site at Cranfield, Mississippi, USA. Synthetic continuous active-source seismic monitoring (CASSM) data (in the form of Vp and Qp measurements) and wellbore pressure monitoring data are assimilated with an ensemble of reservoir realizations to monitor gas saturation and reservoir pressure changes over a period of 100 years. Synthetic seismic attributes are generated using rock physics models (RPMs) and wellbore pressure monitoring data are extracted from the ground truth. Two assimilation methods, ensemble Kalman filter (EnKF) and ensemble Kalman smoother (EnKS), are tested in an observation system simulation experiment (OSSE) environment to assess the prediction accuracy of the individual and composite observation systems. The joint monitoring system achieves more accurate estimates of gas saturation and pressure, across the time span from start of injection to end of forecast, as compared to a single type of monitoring tool and irrespective of data assimilation algorithm choice. These results indicate that jointly assimilated data from two types of sensors (in this case, crosswell seismic and downhole pressure) may lead to a more risk-reducing monitoring design. One would expect that more data, vis-à-vis inclusion of a new sensor type, will improve the accuracy of any GCS monitoring system. However, from a practical standpoint, one important question is whether such a gain in accuracy is worth the additional cost associated with the new sensor. This paper focuses on quantifying the gain in accuracy, such that a practitioner can answer this question.

Engineering↗
NASA NTRS1994

A comparison of CO(J = 1 yields 0) and CO(J = 2 yields 1) emission in the Milky Way molecular ring

We have carried out a CO (J = 2 right arrow 1) survey of the Scutum Arm region of the Milky Way molecular ring. Our goals are to compare CO (J = 2 right arrow 1) maps of individual Galactic clouds with the large-scale CO (J = 2 right arrow 1) emission from the Galactic plane, and to predict the CO (J = 2 right arrow 1) appearance of a Galactic cloud ensemble in an external galaxy. The angular resolution and spatial coverage of our survey are compatible with the existing CO (J = 1 right arrow 0) survey of this region by Sanders et al., which we use for comparison. We identify 34 molecular clouds in our map region; their relationships between size and line width and between virial mass and luminosity are consistent with the relationships seen in CO (J = 1 right arrow 0) emission. However, we note that previous studies have shown considerable variation in these relationships; we attribute much of this variation to differences in their cloud definition algorithm and statistical method. We find that the median ratio of integrated emission in the two lowest CO transitions for the clouds in our sample is I(2 right arrow 1)/ I(1 right arrow 0) = 0.69, implying that the typical emitting region in the line of sight contains cold gas that is not actively star-forming. Our conclusion that the molecular ring emission is not dominated by star-forming regions is consistent with other large-scale studies of the Milky Way. Our distributions of cloud size and temperature also imply that relatively massive molecular clouds that lac k star formation, such as Maddalena's cloud, are not rare in the inner Galaxy.

Chiar, Jean E.↗
DOE OSTI2022

Deep Generative Modeling of Periodic Variable Stars Using Physical Parameters

The ability to generate physically plausible ensembles of variable sources is critical to the optimization of time domain survey cadences and the training of classification models on data sets with few to no labels. Traditional data augmentation techniques expand training sets by reenvisioning observed exemplars, seeking to simulate observations of specific training sources under different (exogenous) conditions. Unlike fully theory-driven models, these approaches do not typically allow principled interpolation nor extrapolation. Moreover, the principal drawback of theory-driven models lies in the prohibitive computational cost of simulating source observables from ab initio parameters. In this work, we propose a computationally tractable machine learning approach to generate realistic light curves of periodic variables capable of integrating physical parameters and variability classes as inputs. Our deep generative model, inspired by the transparent latent space generative adversarial networks, uses a variational autoencoder (VAE) architecture with temporal convolutional network layers, trained using the OGLE-III optical light curves and physical characteristics (e.g., effective temperature and absolute magnitude) from Gaia DR2. A test using the temperature–shape relationship of RR Lyrae demonstrates the efficacy of our generative “physics-enhanced latent space VAE” (PELS-VAE) model. Such deep generative models, serving as nonlinear nonparametric emulators, present a novel tool for astronomers to create synthetic time series over arbitrary cadences.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2012

NASA Center for Climate Simulation (NCCS) Presentation

The NASA Center for Climate Simulation (NCCS) offers integrated supercomputing, visualization, and data interaction technologies to enhance NASA's weather and climate prediction capabilities. It serves hundreds of users at NASA Goddard Space Flight Center, as well as other NASA centers, laboratories, and universities across the US. Over the past year, NCCS has continued expanding its data-centric computing environment to meet the increasingly data-intensive challenges of climate science. We doubled our Discover supercomputer's peak performance to more than 800 teraflops by adding 7,680 Intel Xeon Sandy Bridge processor-cores and most recently 240 Intel Xeon Phi Many Integrated Core (MIG) co-processors. A supercomputing-class analysis system named Dali gives users rapid access to their data on Discover and high-performance software including the Ultra-scale Visualization Climate Data Analysis Tools (UV-CDAT), with interfaces from user desktops and a 17- by 6-foot visualization wall. NCCS also is exploring highly efficient climate data services and management with a new MapReduce/Hadoop cluster while augmenting its data distribution to the science community. Using NCCS resources, NASA completed its modeling contributions to the Intergovernmental Panel on Climate Change (IPCG) Fifth Assessment Report this summer as part of the ongoing Coupled Modellntercomparison Project Phase 5 (CMIP5). Ensembles of simulations run on Discover reached back to the year 1000 to test model accuracy and projected climate change through the year 2300 based on four different scenarios of greenhouse gases, aerosols, and land use. The data resulting from several thousand IPCC/CMIP5 simulations, as well as a variety of other simulation, reanalysis, and observationdatasets, are available to scientists and decision makers through an enhanced NCCS Earth System Grid Federation Gateway. Worldwide downloads have totaled over 110 terabytes of data.

Webster, William P.↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
DOE OSTI2026

Generating entangled steady states in multistable open quantum systems via initial state control

Entanglement underpins the power of quantum technologies, yet it is fragile and typically destroyed by dissipation. Paradoxically, the same dissipation, when carefully engineered, can drive a system toward robust entangled steady states. However, this engineering task is nontrivial, as dissipative many-body systems are complex, particularly when they support multiple steady states. Here, we derive analytic expressions that predict how the steady state of a system evolving under a Lindblad equation depends on the initial state, without requiring integration of the dynamics. These results extend Refs. [V. V. Albert and L. Jiang, Phys. Rev. A 89, 022118 (2014); V. V. Albert et al., Phys. Rev. X 6, 041031 (2016)], showing that while the steady-state manifold is determined by the Liouvillian kernel, the weights within it depend on both the Liouvillian and the initial state. We identify a special class of Liouvillians for which the steady state depends only on the initial overlap with the kernel. Our framework provides analytical insight and a computationally efficient tool for predicting steady states in open quantum systems. As an application, we propose schemes to generate metrologically useful entangled steady states in spin ensembles via balanced collective decay.

Dissipative dynamics↗
DOE OSTI2022

Distilling physical origins of hardness in multi-principal element alloys directly from ensemble neural network models

Abstract Despite a plethora of data being generated on the mechanical behavior of multi-principal element alloys, a systematic assessment remains inaccessible via Edisonian approaches. We approach this challenge by considering the specific case of alloy hardness, and present a machine-learning framework that captures the essential physical features contributing to hardness and allows high-throughput exploration of multi-dimensional compositional space. The model, tested on diverse datasets, was used to explore and successfully predict hardness in Al x Ti y (CrFeNi) 1- x - y , Hf x Co y (CrFeNi) 1- x - y and Al x (TiZrHf) 1- x systems supported by data from density-functional theory predicted phase stability and ordering behavior. The experimental validation of hardness was done on TiZrHfAl x . The selected systems pose diverse challenges due to the presence of ordering and clustering pairs, as well as vacancy-stabilized novel structures. We also present a detailed model analysis that integrates local partial-dependencies with a compositional-stimulus and model-response study to derive material-specific insights from the decision-making process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Discreteness and integrality in Conformal Field Theory

Various observables in compact CFTs are required to obey positivity, discreteness, and integrality. Positivity forms the crux of the conformal bootstrap, but understanding of the abstract implications of discreteness and integrality for the space of CFTs is lacking. We systematically study these constraints in two-dimensional, non-holomorphic CFTs, making use of two main mathematical results. First, we prove a theorem constraining the behavior near the cusp of integral, vector-valued modular functions. Second, we explicitly construct non-factorizable, non-holomorphic cuspidal functions satisfying discreteness and integrality, and prove the non-existence of such functions once positivity is added. Application of these results yields several bootstrap-type bounds on OPE data of both rational and irrational CFTs, including some powerful bounds for theories with conformal manifolds, as well as insights into questions of spectral determinacy. We prove that in rational CFT, the spectrum of operator twists t ≥ c/12 is uniquely determined by its complement. Likewise, we argue that in generic CFTs, the spectrum of operator dimensions Δ > c–1/12 is uniquely determined by its complement, absent fine-tuning in a sense we articulate. Finally, we discuss implications for black hole physics and the (non-)uniqueness of a possible ensemble interpretation of AdS 3 gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Improved measurement of the strong-phase difference $\delta _D^{K\pi }$ in quantum-correlated $D{\overline{D}}$ decays

The decay D → K - π + is studied in a sample of quantum-correlated D$\overline{D}$ pairs, based on a data set corresponding to an integrated luminosity of 2.93 fb -1 collected at the ψ(3770) resonance by the BESIII experiment. The asymmetry between CP-odd and CP-even eigenstate decays into K - π + is determined to be A Kπ = 0.132 ± 0.011 ± 0.007, where the first uncertainty is statistical and the second is systematic. This measurement is an update of an earlier study exploiting additional tagging modes, including several decay modes involving a $K^{0}_{L}$ meson. The branching fractions of the $K^{0}_{L}$ modes are determined as input to the analysis in a manner that is independent of any strong phase uncertainty. Using the predominantly CP-even tag D → π + π - π 0 and the ensemble of CP-odd eigenstate tags, the observable $A^{πππ^0}_{Kπ}$ is measured to be 0.130 ± 0.012 ± 0.008. The two asymmetries are sensitive to $r^{Kπ}_{D}$ cos $δ^{Kπ}_{D}$ where $r^{Kπ}_{D}$ and $δ^{Kπ}_{D}$ are the ratio of amplitudes and phase difference, respectively, between the doubly Cabibbo-suppressed and Cabibbo-favoured decays. In addition, events containing D → K - π + tagged by D → $K^{0}_{S,L}$π + π - are studied in bins of phase space of the three-body decays. This analysis has sensitivity to both $r^{Kπ}_{D}$ cos $δ^{Kπ}_{D}$ and $r^{Kπ}_{D}$, sin $δ^{Kπ}_{D}$. A fit to A Kπ , $A^{πππ^0}_{Kπ}$ and the phase-space distribution of the D → $K^{0}_{S,L}$π + π - tags yields $δ^{Kπ}_{D}$ = ($187.6^{+8.9+5.4}_{-9.7-6.4}$)°, where external constraints are applied for $r^{Kπ}_{D}$ and other relevant parameters. This is the most precise measurement of $δ^{Kπ}_{D}$ in quantum-correlated D$\overline{D}$ decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗
NASA NTRS2004

Formations of Tethered Spacecraft as Stable Platforms for Far IR and Sub-mm Astronomy

In this paper we describe current research in tethered formations for interferometry, and a roadmap to demonstrating the required key technologies via on-ground and in-orbit testing. We propose an integrated kilometer-size tethered spacecraft formation flying concept which enables Far IR and Sub-mm astronomy observations from space. A rather general model is used to predict the dynamics, control, and estimation performance of formations of spacecraft connected by tethers in LEO and deep space. These models include the orbital and tethered formation dynamics, environmental models, and models of the formation estimator/controller/commander. Both centralized and decentralized control/sensing/estimation schemes are possible, and dynamic ranges of interest for sensing/control are described. Key component/subsystem technologies are described which need both ground-based and in-orbit demonstration prior to their utilization in precision space interferometry missions using tethered formations. Defining an orbiting formation as an ensemble of orbiting spacecraft performing a cooperative task, recent work has demonstrated the validity of the tethering the spacecraft to provide both the required formation rigidity and satisfy the formation reconfiguration needs such as interferometer baseline control. In our concept, several vehicles are connected and move along the tether, so that to reposition them the connecting tether links must vary in length. This feature enables variable and precise baseline control while the system spins around the boresight. The control architecture features an interferometer configuration composed of one central combiner spacecraft and two aligned collector spacecraft. The combiner spacecraft acts as the formation leader and is also where the centralized sensing and estimation functions reside. Some of the issues analyzed with the model are: dynamic modes of deformation of the distributed structure, architecture of the formation sensor, and sources of dynamical perturbation that need to be mitigated for precision operation in space. Examples from numerical simulation of an envisioned scenario in heliocentric orbit demonstrate the potential of the concept for space interferometry.

Quadrelli, Marco B.↗
NASA NTRS2013

Advancing Drought Understanding, Monitoring and Prediction

Having the capacity to monitor droughts in near-real time and providing accurate drought prediction from weeks to seasons in advance can greatly reduce the severity of social and economic damage caused by drought, a leading natural hazard for North America. The congressional mandate to establish the National Integrated Drought Information System (NIDIS; Public Law 109-430) in 2006 was a major impulse to develop, integrate, and provide drought information to meet the challenges posed by this hazard. Significant progress has been made on many fronts. On the research front, efforts by the broad scientific community have resulted in improved understanding of North American droughts and improved monitoring and forecasting tools. We now have a better understanding of the droughts of the twentieth century including the 1930s "Dust Bowl"; we have developed a broader array of tools and datasets that enhance the official North American Drought Monitor based on different methodologies such as state-of-the-art land surface modeling (e.g., the North American Land Data Assimilation System) and remote sensing (e.g., the evaporative stress index) to better characterize the occurrence and severity of drought in its multiple manifestations. In addition, we have new tools for drought prediction [including the new National Centers for Environmental Prediction (NCEP) Climate Forecast System, version 2, for operational prediction and an experimental National Multimodel Ensemble] and have explored diverse methodologies including ensemble hydrologic prediction approaches. Broad NIDIS-inspired progress is influencing the development of a Global Drought Information System (GDIS) under the auspices of the World Climate Research Program. Despite these advances, current drought monitoring and forecasting capabilities still fall short of users' needs, especially the need for skillful and reliable drought forecasts at regional and local scales. To tackle this outstanding challenging problem, focused and coordinated research efforts are needed, drawing from excellence across the broad drought research community. To meet this challenge, National Oceanic and Atmospheric Administration (NOAA)'s Drought Task Force was established in October 2011 with the ambitious goal of achieving significant new advances in the ability to understand, monitor, and predict drought over North America. The Task Force (duration of October 2011-September 2014) is an initiative of NOAA's Climate Program Office Modeling, Analysis, Predictions, and Projections (MAPP) program in partnership with NIDIS. It brings together over 30 leading MAPP-funded drought scientists from multiple academic and federal institutions [involves scientists from NOAA's research laboratories and centers, the National Aeronautics and Space Administration (NASA), U.S. Department of Agriculture, National Center for Atmospheric Research (NCAR), and many universities] in a concerted research effort that builds on individual MAPP research projects. These projects span the wide spectrum of drought research needed to make fundamental advances, from those aimed at the basic understanding of drought mechanisms to those aimed at testing new drought monitoring and prediction tools for operational and service purposes (as part of NCEP's Climate Test Bed). The Drought Task Force provides focus and coordination to MAPP drought research activities and also facilitates synergies with other national and international drought research efforts, including those by the GDIS.

GDIS↗
DOE OSTI2021

A violation of global symmetries from replica wormholes and the fate of black hole remnants

We show that the presence of replica wormholes in the Euclidean path integral of gravity leads to a non-perturbative violation of charge conservation for any global symmetry present in the low-energy description of quantum gravity. Explicitly, we compute the scattering probability between different charged states in several two-dimensional models of quantum gravity and find a non-vanishing answer. Here, this suggests that the set of all charged states is typically over-complete, which has drastic consequences for the fate of black hole remnants that could carry a global symmetry charge. In the holographic context, we argue that the presence of such a symmetry in the effective description of the bulk should appear on the boundary as an emergent global symmetry after ensemble averaging.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

Self-Consistent Magnetosphere-Ionosphere Coupling and Associated Plasma Energization Processes

Magnetosphere-Ionosphere (MI) coupling and associated with this process electron and ion energization processes have interested scientists for decades and, in spite of experimental and theoretical research efforts, are still ones of the least well known dynamic processes in space plasma physics. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is large scale ring current (RC) electrodynamic coupling that includes calculation of the magnetospheric electric field that is consistent with the ring current (RC) distribution. A general scheme for numerical simulation of such large-scale magnetosphere-ionosphere coupling processes has been presented earlier in many works. The mathematical formulation of these models are based on "modified frozen-in flux theorem" for an ensemble of adiabatically drifting particles in the magnetosphere. By tracking the flow of particles through the inner magnetosphere, the bounce-averaged phase space density of the hot ions and electrons can be reconstructed and the magnetospheric electric field can be calculated such that it is consistent with the particle distribution in the magnetosphere. The new a self-consistent ring current model has been developed that couples electron and ion magnetospheric dynamics with calculation of electric field. Two new features were taken into account in addition to the RC ions, we solve an electron kinetic equation in our model, self-consistently including these results in the solution. Second, using different analytical relationships, we calculate the height integrated ionospheric conductances as the function of precipitated high energy magnetospheric electrons and ions as produced by our model. This results in fundamental changes to the electric potential pattern in the inner magnetosphere, with a smaller Alfven boundary than previous potential formulations would predict but one consistent with recent satellite observations. This leads to deeper penetration of the plasma sheet ions and electrons into the inner magnetosphere and more effective ring current ions and electron energization.

Khazanov, G. V.↗
NASA NTRSDate not supplied

Openet: Applications of Satellite-Based Evapotranspiration Data for Water Resources Management in the Western United States

Advancing water security in overallocated river basins globally requires consistent and reproducible information on consumptive use of water that can anchor the development of data-driven solutions to the challenge of balancing water supply and demand. OpenET is a fully automated system for field-scale (30 m), satellite-based mapping of evapotranspiration (ET) at daily, monthly and annual timesteps. OpenET currently provides spatially contiguous data throughout the 23 westernmost states in the continental US, and includes both current information as well as multi-year timeseries of ET. The OpenET consortium has implemented an ensemble of satellite-based ET models (ALEXI/DisALEXI, eeMETRIC, PT-JPL, geeSEBAL, SIMS and SSEBop) on Google Earth Engine, which provides a shared computing platform for collaboration on processing of data from Landsat and other satellites, land cover and meteorological inputs, leading to increased consistency and accuracy across the ensemble of models. Earth Engine also facilitates hosting and distribution of data via open data collections and an application programming interface. We provide updates on the OpenET framework, open data services and data access tools, approach to geographic expansion, recent accuracy assessments, and describe how a user-driven design approach has facilitated successful applications of OpenET data for a wide range of water resource management activities. Applications to date include: use of ET data to improve quantification of ET and consumptive use in Oregon, Utah and the Upper Colorado River Basin; streamlining of water use reporting requirements in the California Delta; support for calculation of water budgets for the implementation of the Sustainable Groundwater Management Act in California; and integration into decision support tools for irrigation management. The use cases demonstrate how satellite-derived ET data that are easily accessed and seen as broadly accepted can accelerate adoption of innovative water management practices at scale, and support advances in the sustainability of water supplies. Uptake and use of data by the OpenET science community has also led to advances in our understanding of the impacts of landcover change, irrigation intensification and wildfire events on hydrology and the water security.

Applications↗
DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2022

Skillful Seasonal Forecasts of Land Carbon Uptake in Northern Mid- and High Latitudes

Here we present a first look at the Gross Primary Production (GPP) forecast skill levels achievable with a state-of-the-art subseasonal-to-seasonal (S2S) forecast system. Using NASA’s retrospective S2S ensemble forecast in conjunction with a terrestrial biosphere model, and using an independent, remote sensing-based dataset for validation, we demonstrate an ability to accurately forecast spring-summer carbon uptake at multi-month leads. Averaged across mid-and high latitudes of the Northern Hemisphere land, the GPP forecast initialized on January 1 produces statistically significant skill through summer. The skill achieved, however, is spatially variable, with some regions appearing to extract skill from accurate forecasts of snowpack removal and others extracting skill from the initialization of carbon and nitrogen states. Our results reveal some heretofore unexplored facets of climate predictability and provide a look at what might be possible with future S2S forecast systems that are fully integrated with biogeochemical cycles.

Eunjee Lee↗