Search NASA⌕ Search

SEARCH · Search NASA

Results for “Latent variable”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Performance assessment of a real water source heat pump within a hardware-in-the-loop (HIL) testing environment

Over the last decade, the global fight against climate change through electrification has led to an increase in research on building heating, ventilation, and air conditioning (HVAC) systems that utilize intelligent control algorithms to provide demand-side grid service while also maintaining the thermal comfort of building occupants. As the pivotal point between building electricity consumption and indoor thermal comfort, high-efficiency electrical vapor-compression heat pumps are at the center of these emerging studies, and various grid-interactive and occupant-comfort control algorithms have been developed for them. The impact of these algorithms on heat pump operation and performance when subjected to different weather conditions, building loads, and grid requests calls for investigation and verification via experimental testing with actual heat pumps integrated with real-time building and grid responses. This study introduces a Water-Source Heat Pump (WSHP) Hardware-in-The-Loop (HIL) Test Facility that is the first of its kind. This testbed utilizes a 2-ton variable speed water-to-air heat pump that is capable of interacting with a virtual environment currently comprised of an EnergyPlus (E+) building simulation, an agent-based occupant behavioral model, and a single U-tube ground-loop heat exchanger (GLHE) model. Detailed descriptions of the testbed’s physical design and operation, virtual environment, as well as their mutual communication is provided. An uncertainty analysis is also performed under manufacturer specified heating and cooling design conditions. This analysis shows that the total load across the WSHP’s demand side heat exchanger, i.e., the sum of its latent and sensible components, can be measured with a relative uncertainty of ± 10.4% and ± 3.6% in cooling and heating mode respectively. The WSHP’s coefficient of performance (COP) can be measured with relative uncertainties of ± 10.4% in cooling mode, and ± 3.7% in heating mode. A preliminary 24-h experimental demonstration is then performed utilizing the DOE prototype small commercial office building model in E+. The simulation takes place in Atlanta, GA on the date of 08/26/15 from 12:00 AM to 11:59 PM using TMY3 weather data. Here, the results from this demonstration show that over the course of this experiment the simulated outputs of zone dry-bulb temperature, zone humidity ratio, and WSHP inlet water temperature can be tracked by testbed emulators up to a root mean squared error (RMSE) of ± 0.27 °C, ± 0.376 g/kg, and ± 0.85 °C respectively. The WSHP’s dynamic behavioral characteristics and performance are also captured, and correspond well with the authors’ previous understanding of heat pump efficiency as a function of evaporator and condenser fluid inlet conditions respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Influence of Tundra Polygon Type and Climate Variability on CO 2 and CH 4 Fluxes Near Utqiagvik, Alaska

Arctic tundra has the potential to generate significant climate feedbacks, but spatial complexity makes it difficult to quantify the impacts of climate on ecosystem-atmosphere fluxes, particularly in polygonal tundra comprising wetter and drier polygon types on the scale of tens of meters. We measured CO 2 , CH 4 , and energy fluxes using eddy covariance for 7 yr (April to November, 2013–2019) in polygonal tundra near Utqiagvik, Alaska. This period saw the earliest snowmelt, latest snow accumulation, and hottest summer on record. To estimate fluxes by polygon type, we combined a polygon classification with a flux-footprint model. Methane fluxes were highest in the summer months but were also large during freeze-up and increased with the warming trend in August–November temperatures. While CO 2 respiration had a consistent, exponential relationship with temperature, net ecosystem exchange was more variable among years. CO 2 and CH 4 exchange (June–September) ranged between -0.83 (Standard error [SE] = 0.03) and -1.32 (SE = 0.04) μmol m -2 s-1 and 13.92 (SE = 0.26)—23.42 (SE = 0.45) nmol m -2 s -1 , respectively, and varied interannually (p ≤ 0.05). The maximum-influence method effectively attributed fluxes to polygon types. Areas dominated by low-centered polygons had higher CO 2 fluxes except in 2016–2017. Methane fluxes were highest in low-centered polygons 2013–2015 and in flat-centered polygons in subsequent years, possibly due to increasing temperature and precipitation. Sensible and latent heat fluxes also varied significantly among polygon types. Accurate characterization of Arctic fluxes and their climate dependencies requires spatial disaggregation and long term observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluating the Conservation of Energy Variables in Simulations of Deep Moist Convection

Abstract It is often assumed in parcel theory calculations, numerical models, and cumulus parameterizations that moist static energy (MSE) is adiabatically conserved. However, the adiabatic conservation of MSE is only approximate because of the assumption of hydrostatic balance. Two alternative variables are evaluated here: MSE − IB and MSE + KE, wherein IB is the path integral of buoyancy ( B ) and KE is kinetic energy. Both of these variables relax the hydrostatic assumption and are more precisely conserved than MSE. This article quantifies the errors that result from assuming that the aforementioned variables are conserved in large-eddy simulations (LES) of both disorganized and organized deep convection. Results show that both MSE − IB and MSE + KE better predict quantities along trajectories than MSE alone. MSE − IB is better conserved in isolated deep convection, whereas MSE − IB and MSE + KE perform comparably in squall-line simulations. These results are explained by differences between the pressure perturbation behavior of squall lines and isolated convection. Errors in updraft B diagnoses are universally minimized when MSE − IB is assumed to be adiabatically conserved, but only when moisture dependencies of heat capacity and temperature dependency of latent heating are accounted for. When less accurate latent heat and heat capacity formulae were used, MSE − IB yielded poorer B predictions than MSE due to compensating errors. Our results suggest that various applications would benefit from using either MSE − IB or MSE + KE instead of MSE with properly formulated heat capacities and latent heats.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Indications of a Decrease in the Depth of Deep Convective Cores with Increasing Aerosol Concentration during the CACTI Campaign

Abstract An aerosol indirect effect on deep convective cores (DCCs), by which increasing aerosol concentration increases cloud-top height via enhanced latent heating and updraft velocity, has been proposed in many studies. However, the magnitude of this effect remains uncertain due to aerosol measurement limitations, modulation of the effect by meteorological conditions, and difficulties untangling meteorological and aerosol effects on DCCs. The Cloud, Aerosol, and Complex Terrain Interactions (CACTI) campaign in 2018–19 produced concentrated aerosol and cloud observations in a location with frequent DCCs, providing an opportunity to examine the proposed aerosol indirect effect on DCC depth in a rigorous and robust manner. For periods throughout the campaign with well-mixed boundary layers, we analyze relationships that exist between aerosol variables (condensation nuclei concentration > 10 nm, 0.4% cloud condensation nuclei concentration, 55–1000-nm aerosol concentration, and aerosol optical depth) and meteorological variables [level of neutral buoyancy (LNB), convective available potential energy, midlevel relative humidity, and deep-layer vertical wind shear] with the maximum radar-echo-top height and cloud-top temperature (CTT) of DCCs. Meteorological variables such as LNB and deep-layer shear are strongly correlated with DCC depth. LNB is also highly correlated with three of the aerosol variables. After accounting for meteorological correlations, increasing values of the aerosol variables [with the exception of one formulation of aerosol optical depth (AOD)] are generally correlated at a statistically significant level with a warmer CTT of DCCs. Therefore, for the study region and period considered, increasing aerosol concentration is mostly associated with a decrease in DCC depth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sea Ice and Cloud Processes Mediating Compensation between Atmospheric and Oceanic Meridional Heat Transports across the CMIP6 Preindustrial Control Experiment

Abstract Bjerknes compensation (BJC) refers to the anticorrelation observed between atmospheric and oceanic heat transport (AHT/OHT) variability, particularly on decadal to longer time scales that may be important to the predictability of the climate system. This study investigates the spread in BJC across fully coupled simulations of phase 6 of the Coupled Model Intercomparison Project (CMIP6) and critical processes (particularly related to sea ice and clouds) that may contribute to that spread. BJC on decadal to longer time scales is confirmed across all the simulations evaluated, and it is strongest in the Northern Hemisphere (NH) between 60° and 70°N. At these latitudes, BJC appears to be primarily driven by the exchange of turbulent fluxes (sensible and latent) in the Greenland, Iceland, and Barents Seas. Metrics to break down how sea ice and clouds uniquely modify the radiative balance of the polar atmosphere during anomalous OHT events are presented. These metrics quantify the impacts of sea ice and clouds on surface and top of atmosphere (latent, sensible, longwave, and shortwave radiative) energy fluxes. Cloud responses tend to counter the clear sky impacts over the Marginal Ice Zone (MIZ). It is further shown that the degree of BJC present in a simulation at high latitudes is heavily influenced by the sensitivity of the sea ice to OHT, which is most influential over the MIZ. These results are qualitatively robust across models and explain the intermodel spread in NH BJC in the preindustrial control experiment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Advanced HVAC Humidity Control for Hot-Humid Climates

During this project we develop and validate a cost-effective, integrated control solution to improve humidity control and comfort for energy-efficient homes in hot-humid climates. This study focuses on developing a strategy that is effective, field tested, and practical for builders to install with minimal disruption to standard practices. A successful solution would simplify the transition to high-performance humidity control and be the basis for design and installation guidance. By relying on the central system as a starting point, the strategy employed minimizes system complexity and cost for builders, while improving comfort and operating cost for homeowners. The solution strategy was to coordinate the cooling, dehumidification, and ventilation functions of central, ducted HVAC systems to better control indoor humidity, improve occupant thermal comfort, and capture energy savings. The primary strategic goals were to: (1) optimize dehumidification by the central air-conditioning system, particularly during part-load conditions, using conventional equipment with modified control settings and lower system airflows; (2) maximize ventilation during heating/cooling on-cycles, to “bank” and condition outdoor ventilation air, and minimize ventilation during off-cycles; (3) quantify the effectiveness and energy impact of the dehumidification and ventilation strategies, while identifying a metric that would be useful to evaluate latent effectiveness. For the test houses in our study, located in Richmond Hill, Georgia; Houston, Texas; and Monroe, Louisiana we observed: (1) the indoor humidity did not exceed 60% RH during the monitored cooling season for 99% of the time in Richmond Hill, 96% of the time in Houston, and 90% of the time in Monroe; (2) the dehumidification strategy improved the steady-state latent capacity of the HVAC system at design conditions by 16% to 49% at the Houston test house and by 28% to 71% at the Monroe test house, depending on which mode the system was operating in; and (3) the good results at the test houses were primarily due to the amount of time the air-conditioning system operated in ramping or dehumidification modes, or both, particularly during the early cooling season. This study demonstrates that air conditioners or heat pumps with a single-stage compressor can provide good humidity control without the need for a two-stage or variable-stage compressor system. The airflow and control settings for ramping and dehumidification modes are critical to control indoor humidity in hot-humid climates, particularly during part-load and shoulder season conditions. The dehumidification strategy used in this study did not jeopardize the mechanical reliability of the cooling equipment. The strategies used in this study are applicable across various equipment brands, models, and efficiency levels, and also applicable to a broad range of homes in hot-humid climates. Results will vary by specific equipment, location, and house configuration and construction.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Recharge Demand Mitigation for Latent Heat Thermal Energy Storage at Off-peak Hours

The utilization of integrated heat pump thermal energy storage (HP-TES) systems for space conditioning applications can reduce electrical power demand and achieve utility savings resulting from peak shaving and load shifting. However, issues arise when the energy required to recharge the TES exceeds the shifted energy demand, especially in extreme climates where high-temperature lifts between the outdoor air conditions and phase-change materials (PCMs) during recharging may occur. In this work, a single room-temperature (22°C) PCM-TES is integrated into a single-speed heat pump system using a secondary hydronic loop where the TES is recharged using the outdoor unit during off-peak hours. Three recharge strategies were investigated as potential pathways to reduce recharging power demand: (i) considering different heat source/sink temperatures, (ii) using a variable-speed compressor, and (iii) using a variable-speed pump. Simulations were conducted using Modelica for HP-TES recharging at different heat source/sink temperatures from -30 – 20℃ in heating mode and 0 – 30℃ in cooling mode. Results show that the most effective recharge demand reduction strategy is an HP-TES system using a single-speed compressor with a variable-speed pump for the hydronic loop run at preferable heat source/sink temperatures when the temperature lifts were reduced. Utilizing a variable-speed compressor for the HP-TES system reduced recharge time in extreme cases at the cost of additional compressor power input. The findings motivate potential control strategies that minimize recharge energy and maximize peak energy savings, hence reducing overall operating costs and annual energy consumption for HP-TES systems for space conditioning applications.

25 ENERGY STORAGE↗
DOE OSTI2023

End-Use Savings Shapes Measure Documentation: Variable Refrigerant Flow with Heat Recovery and Dedicated Outdoor Air System

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock™ and ComStock™ models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility measure applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on a single heating, ventilation, and air-conditioning (HVAC) end-use savings shape measure - a variable refrigerant flow with heat recovery (VRF HR) heating and cooling system coupled with a dedicated outdoor air system (DOAS) for ventilation. This measure replaces existing multi-zone variable air volume (VAV) systems or single-zone rooftop units (RTU) with a VRF HR system coupled with a DOAS that includes an energy/heat recovery ventilator (E/HRV). The measure covers 53% of exisiting building stock's floor area and is not applicable to HVAC system types using district heating or cooling or buildings/spaces that include high-ventilation spaces such as kitchens where the amount of exhaust air is large. A DOAS with E/HRV is used to provide required outdoor ventilation air to spaces since ventilation air is generally not supplied by a VRF HR system. An exhaust air energy recovery ventilator (ERV ) with sensible and latent heat exchange is added to humid climate zones while a heat recovery ventilator (HRV ) with sensible only exchange is added to drier climate zones. The ERV is modeled as a fixed membrane plate counterflow heat exchanger, while the HRV is modeled as a sensible-only fixed aluminum plate counterflow heat exchanger. Both systems include a bypass (for temperature control and economizer lockout) and minimum exhaust temperature control for frost prevention.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗
DOE OSTI2025

Diagnostic-free onboard battery health assessment

Diverse usage patterns induce complex and variable aging behaviors in lithiumion batteries, complicating accurate health diagnosis and prognosis. Separate diagnostic cycles are often used to untangle the battery’s current state of health from prior complex aging patterns. However, these same diagnostic cycles alter the battery’s degradation trajectory, are time-intensive, and cannot be practically performed in onboard applications. Here, in this work, we leverage portions of operational measurements in combination with an interpretable machine learning model to enable rapid, onboard battery health diagnostics and prognostics without offline diagnostic testing and the requirement of historical data. We integrate mechanistic constraints within an encoder-decoder architecture to extract electrode states in a physically interpretable latent space and enable improved reconstruction of the degradation path. The health diagnosis model framework can be flexibly applied across diverse application interests with slight fine-tuning.

battery aging reconstruction↗
DOE OSTI2024

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗
DOE OSTI2025

Urban weather modeling using WRF: linking physical assumptions, code implementation, and observational needs

The Weather Research and Forecasting (WRF) model includes urban schemes that simulate the influence of urban surfaces on the atmosphere using parameterizations for flux, and radiative exchanges. Three core schemes – the Bulk urban parameterization, Single-Layer Urban Canopy Model (SLUCM), and Multi-Layer Urban Canopy Model (MLUCM) – represent increasing levels of complexity. Although the parameterizations within these urban schemes are described in the literature, their specific implementation remains poorly documented, thus slowing down model development efforts. This manuscript provides a roadmap to the three urban schemes in WRF version 4.5.2, presenting equations using the same symbols as in the model code, along with references to code lines, and including graphics and explanations that connect the code to its physical foundations. Our thorough review of the urban parameterizations implemented in WRF version 4.5.2 highlighted a handful of parameters that may introduce discontinuities in simulations: (i) in the SLUCM, a 1 mm h −1 rain rate threshold is employed to switch between two minimum moisture availability parameterizations, thus impacting latent heat flux calculations; (ii) in the SLUCM a threshold is used to partition shortwave radiation into direct and diffuse components; (iii) in all three urban schemes, the bulk Richardson number is employed to select the similarity function, which influences the vertical distribution of heat and momentum. We also identified a highly simplified treatment of the radiative balance on roof surfaces. The implications of these simplifications can be assessed through targeted observations across relevant conditions, including varying precipitation rates, cloud cover, and transitions between stability regimes. Furthermore, the widespread application of the Monin-Obukhov similarity theory in these urban schemes warrants model evaluation under highly stable and unstable conditions and in heterogeneous urban settings with variable land cover and building heights on scales finer than model resolution. To address these challenges, we offer guidance on observational strategies, emphasizing the need for multi-parametric measurements to capture potential compensating biases and multi-height measurements that align with the levels where quantities are diagnostic and prognosed in the model (i.e., the lowest atmospheric level of the WRF model). Finally, our inspection of the code revealed implementation bugs that have now been corrected in WRF versions 4.6.0 and 4.6.1. Sensitivity tests over the Atlanta urban area show that these corrections affect surface temperatures, underscoring the importance of performing rigorous documentation and verification of the implementation of parameterizations in model code.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Opinion: A critical evaluation of the evidence for aerosol invigoration of deep convection

Deep convective updraft invigoration via indirect effects of increased aerosol number concentration on cloud microphysics is frequently cited as a driver of correlations between aerosol and deep convection properties. Here, we critically evaluate the theoretical, modeling, and observational evidence for warm- and cold-phase invigoration pathways. Though warm-phase invigoration is plausible and theoretically supported via lowering of the supersaturation with increased cloud droplet concentration in polluted conditions, the significance of this effect depends on substantial supersaturation changes in real-world convective clouds that have not been observed. Much of the theoretical support for cold-phase invigoration depends on unrealistic assumptions of instantaneous freezing and unloading of condensate in growing, isolated updrafts. When applying more realistic assumptions, impacts on buoyancy from enhanced latent heating via fusion in polluted conditions are largely canceled by greater condensate loading. Many foundational observational studies supporting invigoration have several fundamental methodological flaws that render their findings incorrect or highly questionable. Thus, much of the evidence for invigoration has come from numerical modeling, but different models and setups have produced a vast range of results. Furthermore, modeled aerosol impacts on deep convection are rarely tested for robustness, and microphysical biases relative to observations persist, rendering many results unreliable for application to the real world. Without clear theoretical, modeling, or observational support, and given that enervation rather than invigoration may occur for some deep convective regimes and environments, it is entirely possible that the overall impact of cold-phase invigoration is negligible. Substantial mesoscale variability of dominant thermodynamic controls on convective updraft strength coupled with substantial updraft and aerosol variability in any given event are poorly quantified by observations and present further challenges to isolating aerosol effects. Observational isolation and quantification of convective invigoration by aerosols is also complicated by limitations of available cloud condensation nuclei and updraft speed proxies, aerosol correlations with meteorological conditions, and cloud impacts on aerosols. Furthermore, many cloud processes, such as entrainment and condensate fallout, modulate updraft strength and aerosol–cloud interactions, varying with cloud life cycle and organization, but these processes remain poorly characterized. Considering these challenges, recommendations for future observational and modeling research related to aerosol invigoration of deep convection are provided.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Groundwater regulates interannual variations in evapotranspiration in a riparian semiarid ecosystem

Dryland ecosystems are an important component of Earth’s carbon, water, and energy budgets. Changes in hydroclimatic conditions can significantly alter land surface fluxes in these ecosystems. In regions with groundwater-river water exchange, land surface fluxes are strongly modulated by groundwater table fluctuations in response to river stage variations; however, interactions between fluxes and groundwater remain poorly understood in drylands. Here, we use eddy covariance data collected in paired upland and riparian semiarid ecosystems during two years with contrasting river flow conditions to examine interactions among river stage conditions, water table dynamics, meteorological drivers, and ecosystem fluxes. Elevated groundwater availability resulting from groundwater-surface water exchange at the riparian site enhanced the latent heat flux (LE, the energy form of evapotranspiration), which was sustained longer in the growing season during both years. At the upland site, mean LE during the dry season differed by less than 2% between the two years. In contrast, at the riparian site contrasting river flow conditions resulted in a 41% difference in peak monthly LE between the two years. The two-year mean LE in the summer at the riparian site was twice that at the upland site. Carbon uptake (NEP) and LE were well correlated at the upland site during both years, but at the riparian site enhanced summertime LE during the high river flow year did not correspond t enhanced NEP. Our results demonstrate that enhanced groundwater availability modulated by river stage variability is a critical driver of land surface fluxes in semiarid riparian ecosystems.

Missik, Justine E.↗
DOE OSTI2025

Variational AutoEncoders Reveal Intensifying GPP Extremes in Continental United States based on CESM2 Simulations

Climate extremes significantly impact terrestrial carbon cycle dynamics, necessitating robust methods for detecting and analyzing anomalous behavior in plant productivity. This study presents a novel application of variational autoencoders (VAE) for identifying extreme events in gross primary productivity (GPP) from Community Earth System Model version 2 simulations across four AR6 regions in the Continental United States. We compare VAE-based anomaly detection with traditional singular spectral analysis (SSA) methods across three time periods: 1850-80, 1950-80, and 2050-80 under SSP5-8.5 scenario. The VAE architecture employs three dense layers and a latent space with input sequence length of 12 months, training on normalized GPP time series to reconstruct the GPP and identify anomalies based on reconstruction errors. Extreme events are defined using 5th percentile thresholds applied to both VAE and SSA anomalies. Results demonstrate strong regional agreement between VAE and SSA methods in spatial patterns of extreme event frequencies, despite VAE consistently producing higher threshold values (179-756 GgC for VAE vs. 100-784 GgC for SSA across regions and periods). Both methods reveal increasing magnitudes and frequencies of negative carbon cycle extremes toward 2050-80, particularly in Western and Central North America. The VAE approach shows comparable performance to established SSA techniques while offering computational advantages and enhanced capability for capturing non-linear temporal dependencies in carbon cycle variability. This research demonstrates the potential of deep learning approaches for extremes detection and provides a foundation for improved understanding of future carbon cycle risks under future conditions.

Sharma, Bharat [ORNL] (ORCID:0000000266982487)↗
DOE OSTI2025

An integrated integral projection model ( IPM 2 ) to disentangle size‐structured harvest and natural mortality

Abstract Body size is one of the most important traits governing individual‐level demographic rates and modulating population‐level processes. Multiple size‐dependent demographic rates can simultaneously change population structure, so distinguishing their individual contributions to overall population dynamics remains a challenge. Disentangling size‐dependent harvest rates from other demographic rates is critical for assessing the impact of removal on populations of invasive species. Inference about invasive populations can be difficult, however, as observations are often collected opportunistically as part of removal programs, rather than experimentally designed. Yet accurate inference is essential for understanding the feasibility of population suppression and optimising management decisions. We develop an integrated integral projection model (IPM 2 ) that leverages the strengths of the integrated population model and integral projection model to enable inference about complex, size‐structured demographic rates from imperfect observations. We apply the IPM 2 in the context of invasive European green crab ( Carcinus maenas ), a species for which individual body size strongly regulates both the observation‐generating process and latent, population dynamics. The IPM 2 facilitates the distinct estimation of green crab size‐structured harvest and natural mortality rates, parameters for which no explicit data is collected and that are unidentifiable in component datasets of the integrated population model. The model represents how the green crab population changes over time, providing the first estimates of size‐structured abundance of this high‐priority species. By forecasting the stable size distribution and equilibrium population size under varying removal efforts, we demonstrate that extremely high levels of removal effort can reduce the equilibrium green crab population size. Yet these high mortality rates also shift the stable size distribution and increase the equilibrium abundance of smaller crabs, since size‐selective removal alters intraspecific interactions. The ecological outcome of this shift in size structure will be variable, as green crab size modulates only some of its interactions with other species. These results highlight the value of the IPM 2 framework for inferring complex population dynamics with information needs that outpace information in individual observational datasets, providing a path forward for accurate assessment of conservation programs.

Keller, Abigail G. [Department of Environment Scie↗
DOE OSTI2025

Is There a Scalar Atmospheric Surface Layer Within a Convective Boundary Layer? Implications for Flux Measurements

Top‐down entrainment shapes the vertical gradients of sensible heat, latent heat, and CO 2 fluxes, influencing the interpretation of eddy covariance (EC) measurements in the unstable atmospheric surface layer (ASL). Using large eddy simulations for convective boundary layer flows, we demonstrate that decreased temperature gradients across the entrainment zone increase entrainment fluxes by enhancing the entrainment velocity, amplifying the asymmetry between top‐down and bottom‐up flux contributions. These changes alter scalar flux profiles, causing flux divergence or convergence and leading to the breakdown of the constant flux layer assumption (CFLA) in the ASL. As a result, EC‐measured fluxes either underestimate or overestimate “true” surface fluxes during divergence or convergence phases, contributing to energy balance non‐closure. The varying degrees of the CFLA breakdown are a fundamental cause for the non‐closure issue. These findings highlight the underappreciated role of entrainment in interpreting EC fluxes, addressing non‐closure, and understanding site‐to‐site variability in flux measurements.

eddy covariance fluxes of scalars↗
DOE OSTI2019

Eddy Flux and Meteorology over Deciduous Forest, Prairie, and Soybean Ecosystems in Missouri, USA, during the Total Solar Eclipse of 2017

Eddy fluxes and meteorology data are reported at high temporal resolution for three flux tower sites in Mid-Missouri before, during, and after the Total Solar Eclipse of August 21, 2017. Mid-Missouri experienced up to 2 minutes and 40 seconds of totality at around solar noon during the eclipse (Fig. 1). The sites are located in deciduous oak-hickory forest, native prairie, and soybean ecosystems (Table 1 and Fig. 3). Eddy fluxes were computed using a wavelet-based approach that permitted the calculation of two-minute mean fluxes without losing low frequency flux contributions (Fig. 2). During the eclipse, standard meteorological variables were sampled and recorded at 5 second intervals. The two-minute means of fluxes and meteorology are reported. This data product includes three data files in NetCDF (*.nc) format, three data dictionary files and a file-level metadata file in comma separate (*.csv) format, and a user guide in PDF (*.pdf) format.

54 ENVIRONMENTAL SCIENCES↗