Search NASA⌕ Search

SEARCH · Search NASA

Results for “Timing Recovery”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Faster recovery of North Atlantic tropical cyclone-induced cold wakes in recent decades

Intense winds associated with tropical cyclones (TCs) generate surface ocean cooling in their wakes, which can persist for several weeks in their aftermath. While multi-decadal observations of the sea surface have shown a substantial warming of the ocean, long-term changes in cold wake recovery time remain largely unknown. Here we find a trend toward faster recovery of TC cold wakes in the Atlantic main development region (MDR) since 2001. This is due primarily to a decrease in the strength of the North Atlantic trade winds, which reduces evaporative cooling of the ocean. The faster damping of TC cold wakes has led to a significant increase in the intensification of subsequent TCs that encounter lingering wakes from prior TCs, with a magnitude that is about 9% of that from long-term warming of the ocean. Finally, earth system model simulations indicate that the observed decrease in the cold wake recovery time will likely continue into the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of Subetadex-α-methyl, a Polyanionic Cyclodextrin Scaffold, as a Medical Countermeasure against Fentanyl and Related Opioids

Subetadex-α-methyl (SBX-Me), a modified, polyanionic cyclodextrin scaffold, has been evaluated for its utilization as a medical countermeasure (MCM) to neutralize the effects of fentanyl and related opioids. Initial in vitro toxicity assays demonstrate that SBX-Me has a nontoxic profile, comparable to the FDA-approved cyclodextrin-based drug Sugammadex. Pharmacokinetic analysis showed rapid clearance of SBX-Me with an elimination half-life of ~7.4 h and little accumulation in major organs. SBX-Me was also evaluated for its ability to counteract the effects of fentanyl, carfentanil, and remifentanil in rats. Recovery times in rats exposed to sublethal fentanyl doses were found to be shorter when treated with SBX-Me after opioid exposure. The recovery times were reduced from ~35 to ~17 min for fentanyl, ~172 to ~59 min for carfentanil, and ~18 to ~12 min for remifentanil. SBX-Me increased the elimination half-life for fentanyl and remifentanil from 5.37 to 6.42 h and 8.24 to 9.74 h, respectively. These data support SBX-Me as a solid platform from which further research can be launched for the development of a MCM against the effects of fentanyl and its analogs. Furthermore, the data suggests that SBX-Me and other analogs are attractive candidates as broad spectrum opioids targeting MCMs.

Chemistry↗
DOE OSTI2024

Regulation of NMDAR activation efficiency by environmental factors and subunit composition

NMDA receptors (NMDAR) convert the major excitatory neurotransmitter glutamate into a synaptic signal. A key question is how efficiently the ion channel opens in response to the rapid exposure to presynaptic glutamate release. Here, we applied glutamate to single channel outside-out patches and measured the successes of channel openings and the latency to first opening to assay the activation efficiency of NMDARs under different physiological conditions and with different human subunit compositions. For GluN1/GluN2A receptors, we find that various factors, including intracellular ATP and GTP, can enhance the efficiency of activation presumably via the intracellular C-terminal domain. Notably, an energy-based internal solution or increasing the time between applications to increase recovery time improved efficiency. However, even under these optimized conditions and with a 1-s glutamate application, there remained around 10–15% inefficiency. Channel activation became more inefficient with brief synaptic-like pulses of glutamate at 2 ms. Of the different NMDAR subunit compositions, GluN2B-containing NMDARs showed the lowest success rate and longest latency to first openings, highlighting that they display the most distinct activation mechanism. In contrast, putative triheteromeric GluN1/GluN2A/GluN2B receptors showed high activation efficiency. Despite the low open probability, NMDARs containing either GluN2C or GluN2D subunits displayed high activation efficiency, nearly comparable with that for GluN2A-containing receptors. These results highlight that activation efficiency in NMDARs can be regulated by environmental surroundings and varies across different subunits.

He, Miaomiao (ORCID:0000000203179136)↗
DOE OSTI2026

Optimizing Heat Recovery with Storage: Control Validation and Sensitivity Analysis of the Time-Independent Energy Recovery Plant Using Modelica

Heat recovery in large building central plants saves energy but traditionally requires simultaneous heating and cooling. The Time-Independent Energy Recovery (TIER) plant shifts this paradigm by integrating thermal energy storage (TES) to enable heat recovery regardless of concurrent demand, offering a highly efficient, space-saving solution to achieve California’s energy goals. However, its integration of heat recovery chillers, cooling-only chillers, cooling towers, and trim air-source heat pumps (ASHPs) creates growing control and sizing complexity. To overcome this, this study employs high-fidelity Modelica dynamic simulation to validate TIER control sequences and optimize equipment sizing. We translated the written Sequences of Operation into executable Control Description Language (CDL) to test logic against sub-hourly loads. This verification workflow successfully identified and resolved critical vulnerabilities, such as thermal storage freezing and equipment short-cycling, in a virtual environment prior to physical deployment. Then, the study analyzes TIER plant performance across three simulated building types in three locations, and a real building load profile, ensuring variety in heating and cooling loads, and simultaneity factors and explores sizing rules for the TES and ASHP capacity. The analysis shows that the TIER plant operates equipment efficiently leading to a plant SCOP of around 7.5 across all scenarios, higher than a traditional ASHP plant, and a viable pathway to de-risk complex system design and control through simulation to identify optimal designs that maximize energy efficiency, minimize operational costs, and ensure robust operation in varied environmental conditions, thereby facilitating the broader adoption of such a solution for large buildings.

Zanetti, Ettore↗
DOE OSTI2024

An Impedance Humidity Sensor Based on CVD Grown WSe 2 2D Films

2D materials-based planar devices have been sought after for gas and humidity sensor applications owing to their high sensitivity due to their ultrathin nature and increased surface area. Herein, 2D WSe 2 films-based planner devices were fabricated to evaluate their performance for potential application in humidity sensors. The structure and morphology of the WSe 2 sensing films were studied by Raman spectroscopy, photoluminescence, optical microscopy, and scanning electron microscopy (SEM) techniques. The relative humidity-dependent sensing performances of the 2D WSe 2 films were evaluated using current–voltage (I–V), relative impedance, and complex total impedance (electrochemical impedance spectroscopy, EIS) techniques. This sensor showed a very stable and repeatable performance investigated over the period of 18 months. The response and recovery times of the WSe 2 -based impedance sensors were found to be 40 and 30 s, respectively. The WSe 2 sensor exhibits a 2 ± 1% minimum hysteresis at a lower humidity level (50% RH) and around 8 ± 2% maximum hysteresis at relative humidity level (70% RH). This value of hysteresis is in the acceptable range during the initial investigation of any material to explore its potential for application in humidity sensors. The density functional theory (DFT) analysis was performed with the interaction from water molecules of the WSe 2 sample for a better understanding of the sensing mechanism of devices.

2D- WSe2↗
DOE OSTI2026

Nanocrystal-to-Ligand Interfacial Thermal Transport: Plasmonic Indium Tin Oxide to Perylene Diimide Adsorbates

The generation and dissipation of heat in nanocrystals upon optical excitation is a process which may either limit or enhance their performance in certain applications. Using transient absorption spectroscopy, we tracked the flow of heat from plasmonic tin-doped indium oxide (ITO) nanocrystals to surface adsorbed perylenediimide (PDI) molecules upon excitation of the ITO near-infrared plasmon resonance. Here, we rationalize the derivative line shapes observed in the transient absorption bleach features are the result of thermal transfer from the ITO core to surface ligands based on temperature-dependent static absorption studies of the PDI molecules adsorbed on plasmonic ITO. Through a series of pump power dependent measurements, we demonstrate that the PDI heating time is largely fluence independent at the powers measured, while the overall recovery time increases. Elucidating thermal transfer rates at the nanocrystal and organic ligand interface is key for applications such as plasmon mediated photocatalysis, where heating may obfuscate hot carrier transfer processes.

heat dissipation↗
DOE OSTI2026

Electrochemical Nutrient Recovery for the Food–Energy–Water Nexus at Municipal Wastewater Facilities: Multivariate Analyses of Seasonal Sampling and Reactor Performance

Digester-equipped municipal wastewater facilities generate recycle streams with high nutrient loads that increase energy consumption and can cause environmental pollution. The reduction of these loads through electrochemical nutrient recovery (ENR) could enhance the food–energy–water nexus by producing fertilizer (struvite). This study investigated the recovery process through a 1 year sampling of recycle streams and the implementation of nutrient recovery. Time series analyses showed that P (as orthophosphate) concentration was time-variant in digester effluent streams, while N (as ammonia) concentration was time-variant in only the aerobic system. Furthermore, these two nutrient concentrations did not correlate in any of the recycle streams. Subsequent multivariate screening analyses identified anode type, NH 4 + concentration, cathodic potential, P concentration, and temperature as most significant for ENR. Finally, the optimum conditions of cathodic potential, anode area-to-volume ratio, and temperature applied to a real recycle stream resulted in 95% P recovery with 0.03 kWh/kg P. This energy consumption is significantly lower than process energy for conventional P fertilizers (1.1 kWh/kg P) and chemical recovery processes at scale (1.7–12.9 kWh/kg P). Overall, this study recommended specific process controls for nutrient recovery, expanded the variables evaluated for ENR, and demonstrated the ability to significantly impact energy demand associated with P-based fertilizers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

HL-2A's ELM cycle simulations by integrating BOUT++'s drift MHD and transport code

A new integrating model has been developed to couple tokamak edge multiscale magnetohydrodynamic (MHD) events and transport simulations, such as edge-localized mode (ELM) cycles. As a proof of principle, we first start from a set of three-field two-fluid model equations, which includes the pressure, current, and vorticity. Here, the equations are separated into the slowly evolving part of the axisymmetric component by taking a time average of the axisymmetric component. The time-averaged fluxes, which are quadratic in fluctuating quantities, act as driven terms for the time-averaged axisymmetric quantities that determine the plasma transport, and therefore the large-scale evolution of the plasma profiles. Then the HL-2A's ELM cycles are simulated using the model. Good agreements of ELM size and pedestal recovery time have been achieved for the solutions obtained from the coupled simulation compared with experiment. For one ELM cycle simulation, the coupled code can achieve a speedup of a factor of up to 30 over standalone code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Excitons in Hematite Fe 2 O 3 : Short-Time Dynamics from TD-DFT and Non-Adiabatic Dynamics Theories

We present a first-principles study of the short-time dynamics of excitons in hematite Fe 2 O 3 . We used time-dependent density functional theory (TD-DFT) with an underlying DFT+U treatment of electron interactions to characterize the electronic structure of excitons and nonadiabatic molecular dynamics theory (NA-MD) to determine their recombination (electronic ground-state recovery) and relaxation dynamics. Decoherence-corrected trajectory surface hopping approaches in NA-MD simulations yielded recovery times of ~1.1 to 1.8 ns and “higher-lying” exciton relaxation times of ~60 to 70 fs, in accord with experimentally derived lifetimes. With hematite phonons in the range of ~100 to 700 cm –1 , higher-lying excitons relax within one or two oscillations of the phonons before getting trapped into an electron–hole pair Exc-3 structure on the first excited state potential energy surface. This structure resembles already a pair of polarons (electron polaron plus hole polaron) with associated lattice distortions three (3) basal planes away. On longer time scales, the electron–hole bipolaronic pair hops to structures Exc-5, then Exc-7, then Exc-9, ... with the electron polaron and hole polaron separated by 5, 7, 9, ... basal planes in a process of charge separation. The largest frequency phonon ~672 cm –1 for the Exc-3 exciton structure is associated with the electron polaron moiety of the exciton. This phonon is a good candidate for giving rise to the recently observed and reported postexcitation transient IR absorption peak.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pedestal formation via different trajectories in the stability space in response to the timing scan of neutral beam heating in DIII-D

The frequency of type-I ELMs decreases as the initiation of the neutral beam injection (NBI) heating is delayed with respect to the time when plasma current (I p ) reaches flat-top in the ITER Baseline Scenario discharges in DIII-D. Henceforth, the time gap between the NBI initiation and I p flat-top will be referred to as “heating delay.” As the heating delay is modified, pedestal formation follows different trajectories in the edge current density–pedestal pressure gradient (j edge -∇p e ped ) space from the L-H transition toward the first ELM event. During the stationary phase after the first ELM, the ELM frequency (f ELM ) decreases by a factor of ~2 as the heating delay is increased. A longer pedestal recovery time in the inter-ELM period is observed for the low f ELM discharges as compared to the high f ELM discharges. Both low and high f ELM discharges show nearly identical profiles of electron density and temperature and have a similar MHD stability just before an ELM crash. However, a marked difference is observed in the magnetic spectrogram of the high and low f ELM discharges in response to the variation in the heating delay. The main difference is in the 200–400 kHz range of the magnetic spectra. A quasi-coherent mode (QCM) at 220 kHz and weaker broadband fluctuations are observed in the high f ELM discharges, while only strong broadband fluctuations are prevalent in the low f ELM discharges. ELM-synchronized analysis shows that the time evolution of these modes is different for the high and low f ELM discharges. The localization of both these modes is confirmed at the maximum gradient region of the pedestal. We hypothesize that these modes cause important pedestal transport and that the difference in the pedestal recovery of the high and low f ELM discharges is a result of the difference in transport driven by these modes, as they change with changes in the heating delay. It is demonstrated experimentally for the first time that discharges with similar pedestal parameters can carry the history of the heating delay into the stationary phase and that changes in turbulent-driven transport are a likely cause of changes in f ELM observed with variations of heating delay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evaluating ecosystem water use efficiency and recovery dynamics during flash droughts: insights from observations and model simulations

Flash droughts (FD), rapidly emerging in a warming future, disrupt ecosystems, agriculture, and water security. Ecosystem water use efficiency (WUE), the ratio of gross primary production (GPP) to actual evapotranspiration (AET), balances carbon assimilation and water loss. FD rapidly disrupts this balance, making WUE critical for assessing plant stress and recovery. Here, this study investigates the dynamics of landscape-scale WUE, and the components of GPP and AET under FD utilizing both observed data from the Missouri Ozark AmeriFlux site (US-MOz) and version 2 of the U.S. Department of Energy’s Earth, Energy, Exascale System Model (E3SM) Land Model (ELMv2). Observations and simulations reveal GPP as dominant for WUE during earlier FD events (2005, 2007, 2012), shifting to AET in recent events (2014, 2018). This agreement indicates that the ELM can capture the shifting dynamics of GPP and AET in regulating WUE under FD conditions. However, the ELM systematically underestimates both GPP and AET and does so in a manner that does not preserve their ratio. As a result, WUE is also underestimated, suggesting that GPP is more strongly underestimated than AET. Furthermore, the ELM also underestimates the speed of GPP recovery, producing an artificially prolonged GPP recovery time following FD events. Observed environmental drivers such as vapor pressure deficit (VPD), soil moisture (SM), and predawn leaf water potential (PLWP) effectively predict WUE, but ELM primarily highlights SM, underestimating VPD’s role. This study demonstrates that relying solely on soil moisture fails to capture the rapid hydraulic recovery observed in PLWP, underscoring the necessity of integrating plant hydraulics into land surface models to improve flash drought predictability.

Evapotranspiration↗
DOE OSTI2024

Decoupling of peeling and ballooning thresholds for pedestal stability and reduction in ELM frequency via enhanced turbulence with edge electron cyclotron heating in DIII-D

Abstract The edge localized mode (ELM) frequency ( f ELM ) decreased by 63% when electron cyclotron heating (ECH) deposition location is shifted from ρ = 0.4 to ρ = 0.8 in DIII-D discharges where the power ratio between neutral beam injection (NBI) and ECH ( P NBI / P ECH ) is kept at ∼1. The performance of the pedestal in the ECH heated case is compared with a pure NBI reference discharge while keeping the total input power constant. All these discharges are performed at balanced input torque conditions. Furthermore, in the pure NBI discharge a strong decoupling of the peeling–ballooning (PB) thresholds is observed. The PB decoupling is preserved when the ECH is deposited at ρ = 0.8 and P NBI / P ECH ∼ 1, while the thresholds manifest a closed stability boundary when the ECH is deposited at ρ = 0.4. The inter-ELM pedestal recovery time is considerably larger for the ECH at ρ = 0.8 case. Increased pedestal turbulence is observed in beam emission spectroscopy (BES), Doppler backscattering and magnetic diagnostics for the ECH at the ρ = 0.8 case. Strong growth of a TEM-like mode is observed in BES and the mode growth is correlated with the decrease in f ELM . In view of these observations, the increased pedestal turbulence seems to be the plausible reason behind the delayed pedestal recovery following an ELM event in the ECH at ρ = 0.8 case, and the preservation of PB decoupling through temperature pedestal profile widening. TRANSP interpretative simulations show that the ECH at the ρ = 0.8 case is more susceptible to ITG/TEM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shallow Soil Response to a Buried Chemical Explosion With Geophones and Distributed Acoustic Sensing

Abstract Shallow sediments can respond non‐linearly to large dynamic strains and undergo a subsequent healing phase as the material gradually recovers following the passing of seismic waves. This study focuses on the physical changes in the subsurface caused by the shaking from a buried chemical explosion detonated in a borehole in Nevada, USA, as a part of the Source Physics Experiment Phase II. The explosion damaged the shallow subsurface and modified the frequency content recorded by 491 geophones and 2240 Distributed Acoustic Sensing (DAS) channels within 2.5 km from surface ground zero. We observe a gradual shift of resonance frequencies in the 10–25 Hz frequency band in the hours following the explosion and develop a method to characterize the related logarithm‐type healing process of the shallow (i.e., upper ∼25 m) subsurface. We find that stronger levels of ground motion increase the relative degree of damage and duration of the subsurface healing; with the spall region exhibiting the largest degree of damage and longest healing recovery time. We observe coherent spatial patterns of damage with the region located to the southeast of the explosion exhibiting more damage than the southwest region. This study demonstrates that both DAS and co‐located geophones capture similar temporal changes associated with the physical processes occurring in the subsurface, with the high‐density sampling of DAS measurements enabling a new capability to monitor the fine‐scale changes of the Earth's shallow subsurface following the detonation of a buried explosion.

Viens, Loïc↗
DOE OSTI2025

Mechanistic and Mitigation-Strategy Insights into NaCl and CaCl 2 Contamination of Proton-Exchange-Membrane Water Electrolysis Using Continuum Modeling

Cationic contaminants are detrimental to proton-exchange-membrane water electrolyzers (PEMWEs). To obtain insight, a 1-D, nonisothermal, multiphase continuum cell model including cationic contamination is developed. Simulations of steady-state cell performance predict decreased performance due to an increase in kinetic overpotential associated with the hydrogen-evolution reaction, which was attributed to decreased protonic-activity within the cathode catalyst layer from proton supplantation with contaminant cations. The accumulation and extent of cation exchange in the cathode catalyst layer depends on the operating current density due to migration. Simulations of cell recovery of potential suggest that a contaminated cell can recover approximately 78% (450 mV) with 24 h of constant current density operation at 2 A cm –2 , with higher current densities accelerating reduced recovery times. Parametric studies show that anode-side acidification at lower current densities inhibit cation contaminant adsorption, and cathode-side acidification at larger current densities facilitate the expulsion of adsorbed cations; for a cathode-side pH of 6 and 5, the cell can recover an additional 10% and 100% performance, respectively. Overall, the model serves as a framework for modeling other aspects of PEMWE systems to address durability and performance aspects, which can assist in improving the viability of the technology.

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

Community Resilience Through Rapid Restoration Leveraging Distributed Energy Resources (DERs) and Low-Cost Sensors

Equitable and automated bottoms-up power restoration following an extreme event will be demonstrated at a site in Puerto Rico. To do so, the team will develop enhanced grid situational awareness techniques integrating behind-the-meter (BTM) distributed energy resources (DER) discovery, impedance sweeping based outage boundary detection, and feasible restoration path identification algorithms. Resilience metric will be developed and incorporated along with situational awareness information in a distributed Model Predictive Control (MPC)-based restoration optimization algorithm to control and mobilize grid assets. These algorithms will be validated through power hardware-in-the-loop experiments and ultimately, a site demonstration to show that outage recovery time and total recovered load could be improved by >20% over the baseline.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Understanding the temperature and humidity dependence of the SiPM characteristics

A miniature version of the ICAL experiment at the India-based Neutrino Observatory, the mini-ICALis in operation at the IICHEP, Madurai. The commissioning work of a Cosmic Muon Veto detector (CMVD)on top of the mini-ICAL is continued using extruded plastic scintillators with embeddedWLS fibers and the SiPM as a photo-transducer. The CMVD is being built to study thefeasibility of a cosmic muon veto for a shallow-depth neutrino experiment. The SiPM is calibratedusing an ultrafast LED driver andthe same calibration techniques will be used during the data taking period with a simplified LED in thesame readout PCB. An experimental setup wasdesigned to characterise the SiPMs in a temperature controlled environment. The readout electronicsinvolves trans-impedance amplifiers and a voltage followers of combined gain 1.24 mV/$\mu$A and adigital storage oscilloscope for the data collection with a minimal distortion of SiPM signal.Various characteristics of the Hamamatsu SiPM (S13360-2050VE), e.g. signal shape, optically correlatedand uncorrelated noise, recovery time etc were studied as a function of $V_{ov}$, number of photoelectrons,the ambient temperature and the humidity. This paper will cover the details of those results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Optimization of electroporation method and promoter evaluation for type-1 methanotroph, Methylotuvimicrobium alcaliphilum

Methanotrophic bacteria are promising hosts for methane bioconversion to biochemicals or bioproducts. However, due to limitations associated with long genetic manipulation timelines and, lack of choice in genetic tools required for strain engineering, methanotrophs are currently not employed for bioconversion technologies. In this study, a rapid and reproducible electroporation protocol is developed for type 1 methanotroph, Methylotuvimicrobium alcaliphilum using common laboratory solutions, analyzing optimal electroshock voltages and post-shock cell recovery time. Successful reproducibility of the developed method was achieved when different replicative plasmids were assessed on lab adapted vs. wild-type M. alcaliphilum strains (DASS vs. DSM19304). Overall, a ∼ 3-fold decrease in time is reported with use of electroporation protocol developed here, compared to conjugation, which is the traditionally employed approach. Additionally, an inducible (3-methyl benzoate) and a constitutive (sucrose phosphate synthase) promoter is characterized for their strength in driving gene expression.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Ecosystem leaf area, gross primary production, and evapotranspiration responses to wildfire in the Columbia River basin

Wildfires impact vegetation mortality and productivity and are increasing in intensity, frequency, and spatial area in the western United States. The rates of vegetation recovery after fires play a major role in the reestablishment of biomass and ecosystem functioning (e.g., structure, resilience, and productivity), but such recovery rates are poorly understood. Here we use remotely sensed data products from the Moderate Resolution Imaging Spectroradiometer (MODIS) to quantify the resistance and resilience of leaf area index (LAI), gross primary production (GPP), and evapotranspiration (ET) to 138 wildfires of various burn severity across the Columbia River basin (CRB) of the Pacific Northwest in 2015. Increasing burn severity caused lower resistance and resilience for all three variables. Resistance and resilience are highest in grasslands, intermediate in savanna, and lowest in needleleaf evergreen forests, consistent with the adaptation of these vegetation types to fire. LAI has consistently lower resistance and resilience than GPP and ET, which is consistent with physical and physiological mechanisms that compensate for reduced LAI. Resilience is influenced by precipitation, vapor pressure deficit (VPD), and burn severity across all three vegetation types; however, burn severity plays a more minor role in grasslands. Increasing wildfire severity will reduce the resistance and resilience and lengthen the recovery time of vegetation structure and fluxes with climate change, with significant consequences for the provision of ecosystem functioning and implications for model predictions.

54 ENVIRONMENTAL SCIENCES↗