Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “balance function”

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 289 records · Page 16

MICCO: An Enhanced Multi-GPU Scheduling Framework for Many-Body Correlation Functions

Calculation of many-body correlation functions is one of the critical kernels utilized in many scientific computing areas, especially in Lattice Quantum Chromodynamics (Lattice QCD). It is formalized as a sum of a large number of contraction terms each of which can be represented by a graph consisting of vertices describing quarks inside a hadron node and edges designating quark propagations at specific time intervals. Due to its computation- and memory-intensive nature, real-world physics systems (e.g., multi-meson or multi-baryon systems) explored by Lattice QCD prefer to leverage multi-GPUs. Different from general graph processing, many-body correlation function calculations show two specific features: a large number of computation-/data-intensive kernels and frequently repeated appearances of original and intermediate data. The former results in expensive memory operations such as tensor movements and evictions. The latter offers data reuse opportunities to mitigate the data-intensive nature of many-body correlation function calculations. However, existing graph-based multi-GPU schedulers cannot capture these data-centric features, thus resulting in a sub-optimal performance for many-body correlation function calculations. To address this issue, this paper presents a multi-GPU scheduling framework, MICCO, to accelerate contractions for correlation functions particularly by taking the data dimension (e.g., data reuse and data eviction) into account. This work first performs a comprehensive study on the interplay of data reuse and load balance, and designs two new concepts: local reuse pattern and reuse bound to study the opportunity of achieving the optimal trade-off between them. Based on this study, MICCO proposes a heuristic scheduling algorithm and a machine-learning-based regression model to generate the optimal setting of reuse bounds. Specifically, MICCO is integrated into a real-world Lattice QCD system, Redstar, for the first time running on multiple GPUs. The evaluation demonstrates MICCO outperforms other state-of-art works, achieving up to 2.25× speedup in synthesized datasets, and 1.49× speedup in real-world correlation functions.

Wang, Qihan↗

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Immunophenotypic characterization and therapeutics effects of human bone marrow- and umbilical cord-derived mesenchymal stromal cells in an experimental model of sepsis

Highlights: • Mesenchymal stromal cells (MSCs) may be a therapeutic option for sepsis. • Bone marrow (BM) vs. umbilical cord (UC) -MSCs induce different M1/M2 gene profiles. • In a CLP model of sepsis, only BM-MSCs improved survival and bacterial clearance. Sepsis is a complicated multi-system disorder characterized by a dysregulated host response to infection. Despite substantial progress in the understanding of mechanisms of sepsis, translation of these advances into clinically effective therapies remains challenging. Mesenchymal Stromal Cells (MSCs) possess immunomodulatory properties that have shown therapeutic promise in preclinical models of sepsis. The therapeutic effects of MSCs may vary depending on the source and type of these cells. In this comparative study, the gene expression pattern and surface markers of bone marrow-derived MSCs (BM-MSCs) and umbilical cord-derived MSCs (UC-MSCs) as well as their therapeutic effects in a clinically relevant mouse model of polymicrobial sepsis, cecal ligation and puncture (CLP), were investigated. The results showed remarkable differences in gene expression profile, surface markers and therapeutic potency in terms of enhancing survival and pro/anti-inflammatory responses between the two MSC types. BM-MSCs improved survival concomitant with an enhanced systemic bacterial clearance and improved inflammatory profile post CLP surgery. Despite some improvement in the inflammatory profile of the septic animals, treatment with UC-MSCs did not enhance survival or bacterial clearance. Overall, the beneficial therapeutic effects of BM-MSCs over UC-MSCs may likely be attributed to their pro-inflammatory function, and to some extent anti-inflammatory features, reflected in their gene expression pattern enhancing macrophage polarization to M1/M2 phenotypes resulting in a balanced pro- and anti-inflammatory response against polymicrobial sepsis.

60 APPLIED LIFE SCIENCES↗

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗

Metabolic Redox Coupling Controls Methane Production in Permafrost‐Affected Peatlands Through Organic Matter Quality‐Dependent Energy Allocation

ABSTRACT Permafrost thaw represents one of Earth's largest climate feedback risks, potentially releasing vast carbon (C) stores as greenhouse gases (GHG). However, our ability to predict emissions remains limited by poor understanding of how changing organic matter (OM) composition affects microbial carbon processing. We test a metabolism‐centered redox framework, which views microbial processes as coupled oxidative‐reductive reactions, to mechanistically explain how organic matter metabolite quality controls greenhouse gas production in permafrost‐affected peatland ecosystems. Rather than relying solely on geochemical redox measurements, our approach examines how microbes balance electron flow through metabolic pathways. Using active layer peat (9–19 cm) from contrasting environments (bog and fen), we employed multi‐omics approaches, including metabolomics, metagenomics, and metatranscriptomics, to link OM chemistry to microbial function. Our results reveal distinct dissolved organic matter metabolite composition, with fen systems enriched in compounds with higher substrate quality (low molecular weight (MW) sugars with high H:C ratios and low aromaticity) and bog systems dominated by compounds with lower substrate quality (high MW phenols with lower H:C ratios and higher aromaticity). In fen samples, these sugar‐like compounds correlated with higher oxidative metabolism and methanogenesis, supported by increased glycolysis gene expression. Initially, electrons from increased oxidative metabolism were balanced through nitrate and sulfate reduction, but as these electron acceptors were depleted, methanogenesis increased to maintain redox balance. Fen samples showed rapid degradation of both high‐ and low‐substrate‐quality compounds, suggesting sufficient energy for efficient C cycling. Conversely, bog samples exhibited more polyphenolic compounds, lower glycolysis activity, and higher stress‐related gene expression, suggesting energy was diverted towards cell maintenance under acidic conditions rather than C processing. This approach suggests that predicting greenhouse gas emissions requires an understanding of how organic matter quality shapes microbial energy allocation strategies, providing a mechanistic framework for improving emission predictions from permafrost‐affected peatlands and similar ecosystems.

Biodiversity & Conservation↗

Device Performance of Emerging Photovoltaic Materials (Version 2)

Following the 1st release of the "Emerging photovoltaic (PV) reports", the best achievements in the performance of emerging photovoltaic devices in diverse emerging photovoltaic research subjects are summarized, as reported in peer-reviewed articles in academic journals since August 2020. Updated graphs, tables, and analyses are provided with several performance parameters, e.g., power conversion efficiency, open-circuit voltage, short-circuit current density, fill factor, light utilization efficiency, and stability test energy yield. These parameters are presented as a function of the photovoltaic bandgap energy and the average visible transmittance for each technology and application and are put into perspective using, e.g., the detailed balance efficiency limit. The 2nd installment of the "Emerging PV reports" extends the scope toward tandem solar cells and presents the current state-of-the-art in tandem solar cell performance for various material combinations.

36 MATERIALS SCIENCE↗

Device Performance of Emerging Photovoltaic Materials (Version 3)

Abstract Following the 2nd release of the “Emerging PV reports,” the best achievements in the performance of emerging photovoltaic devices in diverse emerging photovoltaic research subjects are summarized, as reported in peer‐reviewed articles in academic journals since August 2021. Updated graphs, tables, and analyses are provided with several performance parameters, e.g., power conversion efficiency, open‐circuit voltage, short‐circuit current density, fill factor, light utilization efficiency, and stability test energy yield. These parameters are presented as a function of the photovoltaic bandgap energy and the average visible transmittance for each technology and application, and are put into perspective using, e.g., the detailed balance efficiency limit. The 3rd installment of the “Emerging PV reports” extends the scope toward triple junction solar cells.

14 SOLAR ENERGY↗

Adaptive Sampling-Based Bi-Fidelity Stochastic Trust Region Method for Stochastic Derivative-Free Optimization

Bi-fidelity stochastic optimization has gained increasing attention as an efficient approach to reduce computational costs by leveraging a low-fidelity (LF) model to optimize an expensive high-fidelity (HF) objective. In this paper, we propose ASTRO-BFDF, an adaptive sampling trust-region method specifically designed for unconstrained bi-fidelity stochastic derivative-free optimization problems. In ASTRO-BFDF, the LF function serves two purposes: (i) to identify better iterates for the HF function when the optimization process indicates a high correlation between them and (ii) to reduce the variance of the HF function estimates using bi-fidelity Monte Carlo (BFMC). The algorithm dynamically determines sample sizes while adaptively choosing between crude Monte Carlo and BFMC to balance the trade-off between optimization and sampling errors. We prove that the iterates generated by ASTRO-BFDF converge to a first-order stationary point almost surely. Additionally, we demonstrate the effectiveness of the proposed algorithm through numerical experiments on synthetic benchmarks and simulation optimization problems involving discrete event systems.

97 MATHEMATICS AND COMPUTING↗

Impacts of a Prescribed Fire on Air Quality in Central New Mexico

A short-duration but high-impact air quality event occurred on 28 November 2018 along the Rio Grande Valley of New Mexico. This fire occurred outside the typical wildfire season, and greatly impacted the air quality in Socorro, NM, and the surroundings. Measurements were taken during the event using an aerosol light scattering technique (integrating nephelometer) and a particulate mass concentration monitor (DustTrak PM optical monitor). The instruments sampled the ambient air during the event on the campus of the New Mexico Institute of Mining and Technology in Socorro, New Mexico. The peak values on a 5-min basis of light scattering and the PM mass concentration reached 470 Mm -1 and 270 µg/m 3 , respectively. We examined the meteorological context of the event using local meteorological data and back trajectories using the NOAA HYSPLIT model to determine atmospheric transport and possible sources. Several fires, both prescribed and wildfires, occurred in the region including a prescribed burn at Bosque del Apache National Wildlife Refuge (17 km south-southeast of the receptor site). The data suggest that the prescribed burn at Bosque del Apache was the dominant contributor due to transport evidence and the event’s narrow spatiotemporal extent. The increasing importance of restoring ecosystem function using prescribed fire in wildland fire management will likely lead to more frequent air quality impacts and sets up policy tradeoffs that require a balance between these public goals. This study examines the evidence of the effects of a prescribed fire in a protected wildland area impacting the air quality in a nearby populated area.

54 ENVIRONMENTAL SCIENCES↗

Computational prediction of ferromagnetic 𝐴⁢𝑇 6 ⁢𝑋 6 kagome compounds

We present a systematic high-throughput density-functional-theory investigation of the structural and magnetic stability of 312 substitutional compounds in the magnetic kagome 𝐴⁢𝑇 6⁢ 𝑋 6 family. Our screening confirms the stability of many previously reported structures and predicts several additional stable candidates. Within collinear spin configurations, we find that Fe-based systems predominantly adopt antiferromagnetic ground states, whereas Mn-based analogs exhibit a more balanced distribution between ferromagnetic and antiferromagnetic order. For compounds exhibiting several nearly degenerate collinear configurations, we analyze the nature of their magnetic ground states, assess the possible emergence of noncollinear order, and discuss the limitations and uncertainties inherent to standard density-functional approaches. Our electronic structure analysis further reveals that predicted ferromagnetic kagome systems display characteristic features of topological metals, with rich magnetic configurations that can be tuned by chemical substitution. Altogether, these ferromagnetic kagome compounds constitute a broad and still largely unexplored materials platform for the emergence of exciting magnetotransport phenomena.

Density functional calculations↗

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plating and Stripping of Lithium Metal Stabilized by a Block Copolymer Electrolyte: Local Current Density Measurement and Modeling

Here, plating and stripping of lithium protrusions in lithium metal symmetric cells containing a solid block copolymer electrolyte was studied as a function of time in 3D using time-resolved X-ray tomography. These measurements enabled determination of the spatial variation in current densities at the plating and stripping electrodes. The initial interelectrode distance was 27 μm. Correlation functions were calculated to reveal the relationships between current densities at the two electrodes and local electrolyte thickness. Current densities at opposing electrode locations during protrusion growth is uncorrelated until the local interelectrode distance decreases to less than 6 μm, just before the cell shorts. Mass balance was used to determine the area from which lithium ions that form a protrusion were stripped. Computational modeling of the plating and stripping process reveals the interplay between electrochemical and mechanical driving forces and their effect on nonuniform current distribution. Model predictions were compared with experiments without resorting to any adjustable parameters. The computed correlation functions were in qualitative agreement with experiments. Finally, the model was used to calculate contour plots of electrochemical potential within the electrolyte, shedding light on how geometry, salt concentration, interelectrode distance, and mechanical stress influence local rates of electrochemical reaction.

36 MATERIALS SCIENCE↗

Electrically controllable chirality in a nanophotonic interface with a two-dimensional semiconductor

Chiral nanophotonic interfaces enable propagation direction-dependent interactions between guided optical modes and circularly dichroic materials. Electrical tuning of interface chirality would aid active, switchable non-reciprocity in on-chip optoelectronic and photonic circuitry, but remains an outstanding challenge. Here, in this paper, we report electrically controllable chirality in a nanophotonic interface with atomically thin monolayer tungsten diselenide (WSe 2 ). Titanium dioxide waveguides are directly fabricated on the surface of low-disorder, boron nitride-encapsulated WSe 2 . Following integration, photoluminescence from excitonic states into the waveguide can be electrically switched between balanced and directionally biased emission. The operational principle leverages the doping-dependent valley polarization of excitonic states in WSe 2 . Furthermore, the nanophotonic waveguide can function as a near-field source for diffusive exciton fluxes, which display valley and spin polarizations that are inherited from the interface chirality. Our versatile fabrication approach enables the deterministic integration of photonics with van der Waals heterostructures and could provide optical control over their excitonic and charge-carrier behaviour.

42 ENGINEERING↗

Nonintuitive Surface Self-Assembly of Functionalized Molecules on Ag(111)

The fabrication of nanomaterials involves selfordering processes of functional molecules on inorganic surfaces. To obtain specific molecular arrangements, a common strategy is to equip molecules with functional groups. However, focusing on the functional groups alone does not provide a comprehensive picture. Especially at interfaces, processes that govern self-ordering are complex and involve various physical and chemical effects, often leading to unexpected structures, as we showcase here on the example of a homologous series of quinones on Ag(111). Naively, one could expect that such quinones, which all bear the same functionalization, form similar motifs. In salient contrast, our joint theoretical and experimental study shows that profoundly different structures are formed. Using a machine-learning-based structure search algorithm, we find that this is due to a shift of the balance of three antagonizing driving forces: adsorbate-substrate interactions governing adsorption sites, adsorbate-adsorbate interactions favoring close packing, and steric hindrance inhibiting certain otherwise energetically beneficial molecular arrangements. The theoretical structures show excellent agreement with our experimental characterizations of the organic/ inorganic interfaces, both for the unit cell sizes and the orientations of the molecules within. The nonintuitive interplay of similarly important interaction mechanisms will continue to be a challenging aspect for the design of functional interfaces. With a detailed examination of all driving forces, we are, however, still able to devise a design principle for self-assembly of functionalized molecules.

structure prediction↗

Model-based Cross-correlation Search for Gravitational Waves from the Low-mass X-Ray Binary Scorpius X-1 in LIGO O3 Data

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO and Advanced Virgo. This is a semicoherent search that uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25 to 1600 Hz, as well as ranges in orbital speed, frequency, and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100 and 200 Hz, correspond to an amplitude h $_{0}$ of about 10$^{−25}$ when marginalized isotropically over the unknown inclination angle of the neutron star’s rotation axis, or less than 4 × 10$^{−26}$ assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically marginalized upper limits are close to the predicted amplitude from about 70 to 100 Hz; the limits assuming that the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40 to 200 Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500 Hz or more.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic response of photorespiration in fluctuating light environments

Abstract Photorespiration is a dynamic process that is intimately linked to photosynthetic carbon assimilation. There is a growing interest in understanding carbon assimilation during dynamic conditions, but the role of photorespiration under such conditions is unclear. In this review, we discuss recent work relevant to the function of photorespiration under dynamic conditions, with a special focus on light transients. This work reveals that photorespiration is a fundamental component of the light induction of assimilation where variable diffusive processes limit CO2 exchange with the atmosphere. Additionally, metabolic interactions between photorespiration and the C3 cycle may help balance fluxes under dynamic light conditions. We further discuss how the energy demands of photorespiration present special challenges to energy balancing during dynamic conditions. We finish the review with an overview of why regulation of photorespiration may be important under dynamic conditions to maintain appropriate fluxes through metabolic pathways related to photorespiration such as nitrogen and one-carbon metabolism.

Plant Sciences↗

Filling the Gaps: A Bayesian Mixture Model for Imputing Missing Soil Water Content Data

ABSTRACT Soil water content (SWC) data are central to evaluating how soil moisture varies over time and space and influences critical plant and ecosystem functions, especially in water‐limited drylands. However, sensors that record SWC at high frequencies often malfunction, leading to incomplete timeseries and limiting our understanding of dryland ecosystem dynamics. We developed an analytical approach to impute missing SWC data, which we tested at six eddy flux tower sites along an elevation gradient in the southwestern United States. We impute missing data as a mixture of linearly interpolated SWC between the observed endpoints of a missing data gap and SWC simulated by an ecosystem water balance model (SOILWAT2). Within a Bayesian framework, we allowed the relative utility (mixture weight) of each component (linearly interpolated vs. SOILWAT2) to vary by depth, site and gap characteristics. We explored “fixed” weights versus “dynamic” weights that vary as a function of cumulative precipitation, average temperature, and time since the start of the gap. Both models estimated missing SWC data well ( R 2 = 0.70–0.88 vs. 0.75–0.91 for fixed vs. dynamic weights, respectively), but the utility of linearly interpolated versus SOILWAT2 values depended on site and depth. SOILWAT2 was more useful for more arid sites, shallower depths, longer and warmer gaps and gaps that received greater precipitation. Overall, the mixture model reliably gap‐fills SWC, while lending insight into processes governing SWC dynamics. This approach to impute missing data could be adapted to accommodate more than two mixture components and other types of environmental timeseries.

Ogle, Kiona [School of Informatics, Computing, and↗