Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data reduction”

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 1,351 records · Page 75

Throughfall-Reduction Drying Effects on Gravimetric Soil Moisture at Two Depths in Four Lowland Panamanian Forests from 2015-2022

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We also surveyed all tree species and identified to species in 2018 and 2019.Results: We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests to 20cm depths, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Datasets included here include .csv and .xls files for gravimetric soil moisture (weight/weight). Soil moisture was collected on a quarterly basis from 0-10 cm and 10-20 cm depths using hand-held constant-volume soil corers. Data are on ~3 month timescales from 2015-2022 with some gaps. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗

Discrimination of electronic recoils from nuclear recoils in two-phase xenon time projection chambers

We present a comprehensive analysis of electronic recoil vs. nuclear recoil discrimination in liquid/gas xenon time projection chambers, using calibration data from the 2013 and 2014-16 runs of the Large Underground Xenon (LUX) experiment. We observe strong charge-to-light discrimination enhancement with increased event energy. For events with S1 = 120 detected photons, i.e. equivalent to a nuclear recoil energy of ~100 keV, we observe an electronic recoil background acceptance of <10 –5 at a nuclear recoil signal acceptance of 50%. We also observe modest electric field dependence of the discrimination power, which peaks at a field of around 300 V/cm over the range of fields explored in this study (50-500 V/cm). In the WIMP search region of S1 = 1-80 phd, the minimum electronic recoil leakage we observe is (7.3 ± 0.6) 10 –4 , which is obtained for a drift field of 240-290 V/cm. Pulse shape discrimination is utilized to improve our results, and we find that, at low energies and low fields, there is an additional reduction in background leakage by a factor of up to 3. We develop an empirical model for recombination fluctuations which, when used alongside the Noble Element Scintillation Technique (NEST) simulation package, correctly reproduces the skewness of the electronic recoil data. We use this updated simulation to study the width of the electronic recoil band, finding that its dominant contribution comes from electron-ion recombination fluctuations, followed in magnitude of contribution by fluctuations in the S1 signal, fluctuations in the S2 signal, and fluctuations in the total number of quanta produced for a given energy deposition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Active learning for the design of polycrystalline textures using conditional normalizing flows

Generative modeling has opened new avenues for solving previously intractable materials design problems. However, these new opportunities are accompanied by a drastic increase in the required amount of training data. This is in stark juxtaposition to the high expense and difficulty in curating such large materials datasets. In this work, we propose a novel framework for integrating generative models within an active learning loop. Further, this enables the training of generative models with datasets significantly smaller than what has previously been demonstrated, providing a direct route for their application in data constrained environments. The functionality of this framework is then demonstrated by addressing the challenge of designing polycrystalline textures associated with target anisotropic mechanical properties. The developed protocol exhibited a cost reduction between 14 to 18 times over a randomly sampled experimental design.

36 MATERIALS SCIENCE↗

A fast and accurate domain decomposition nonlinear manifold reduced order model

Here, this paper integrates nonlinear-manifold reduced order models (NM-ROMs) with domain decomposition (DD). NM ROMs approximate the full order model (FOM) state in a nonlinear-manifold by training a shallow, sparse autoencoder using FOM snapshot data. These NM-ROMs can be advantageous over linear-subspace ROMs (LS-ROMs) for problems with slowly decaying Kolmogorov n-width. However, the number of NM-ROM parameters that need to be trained scales with the size of the FOM. Moreover, for “extreme-scale” problems, the storage of high-dimensional FOM snapshots alone can make ROM training expensive. To alleviate the training cost, this paper applies DD to the FOM, computes NM-ROMs on each subdomain, and couples them to obtain a global NM-ROM. This approach has several advantages: Subdomain NM-ROMs can be trained in parallel, involve fewer parameters to be trained than global NM-ROMs, require smaller subdomain FOM dimensional training data, and can be tailored to subdomain specific features of the FOM. The shallow, sparse architecture of the autoencoder used in each subdomain NM-ROM allows application of hyper-reduction (HR), reducing the complexity caused by nonlinearity and yielding computational speedup of the NM-ROM. This paper provides the first application of NM-ROM (with HR) to a DD problem. In particular, this paper details an algebraic DD reformulation of the FOM, training a NM-ROM with HR for each sub domain, and a sequential quadratic programming (SQP) solver to evaluate the coupled global NM-ROM. Theoretical convergence results for the SQP method and a priori and a posteriori error estimates for the DD NM-ROM with HR are provided. The proposed DD NM-ROM with HR approach is numerically compared to a DD LS-ROM with HR on the 2D steady-state Burgers’ equation, showing an order of magnitude improvement in accuracy of the proposed DD NM-ROM over the DD LS-ROM.

97 MATHEMATICS AND COMPUTING↗

Monitoring Fracture Saturation With Internal Seismic Sources and Twin Neural Networks

Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply scattered or reflected waves. However, a challenge arises when multiply scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural or induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (TNN also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The TNN with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. Furthermore, these results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.

58 GEOSCIENCES↗

Strong Sensitivity of Simulated Biomass Burning Aerosol Transport and Radiative Effects Over the South Atlantic to Carbonaceous Aerosol Aging and Particle Density

Biomass burning aerosol (BBA) impacts climate through aerosol‐cloud‐radiation interactions, but models disagree on the sign and magnitude of BBA radiative effects. We quantify the sensitivity of BBA radiative effects and transport to three BBA‐relevant processes and properties: parameterized oxidative aging of organic aerosol (OA), a combined change to black carbon (BC) density and the method for calculating aerosol refractive index, and reduction in OA density. We evaluate Unified Model simulations against two aircraft campaigns from summer 2017 over the Southeast Atlantic. The model generally performs well, such that discrepancies between the observational data sets may sometimes limit the precision of the evaluation. Our newly developed aging parameterization reproduces observed OA:BC mass ratios well and allows modeled OA:BC to decrease with smoke age, but increases bias in aerosol extinction and changes the BBA radiative effect little (+0.12 W m -2 ). We calculate aerosol refractive index using either a volume‐weighted component average or the Maxwell‐Garnett (MG) mixing assumption, which represents BC as small inclusions in a host material. Compared to MG mixing, the volume‐weighted average refractive index and reduced BC density increase aerosol absorption, substantially increasing the total BBA radiative effect (+2.66 W m -2 ) and amount of BBA transported across the ocean through BC self‐lofting. Reducing OA density to better match literature values changes the total BBA radiative effect by −1.96 W m -2 . Changes to direct radiative effects exceed changes to cloud radiative effects. Our findings emphasize the sensitivity of aerosol radiative effects and transport to these processes and properties, which we suggest could be improved in climate models.

Southeast Atlantic↗

Can ferric-oxyl excited states explain elongated iron-oxygen bonds in heme peroxidase catalytic intermediates?

The use of X-ray structures to determine and interpret the ferryl iron-oxygen bond order in molecular oxygen-activating heme enzymes has, in the past, been controversial. This has mainly stemmed from the susceptibility of ferryl species to X-ray-induced electronic state changes. In this work we establishe using time-resolved serial femtosecond X-ray crystallography (tr-SFX) on a dye-decolourising peroxidase that the ferryl intermediate species (Compounds I and II) captured following in situ mixing of microcrystals with H 2 O 2 have single, rather than the double bond character expected. X-ray emission validated tr-SFX data with quantum refinement, time-dependent-DFT calculations and QM/MM geometry optimizations together support the concept that the single iron-oxygen bond character is not an indication of ferryl reduction or a protonated form (Fe IV -OH) but is instead attributed to the existence of accessible excited states possessing ferric-oxyl (Fe III –O •– ) character. Such states offer insight into the nature of ferryl heme.

Williams, Lewis J. [University of Essex, Colcheste↗

Charging-management And Infrastructure-planning (cmip) Model

CMIP model explores various charging infrastructure network designs to serve a free-floating car-sharing fleet and determine the charging downtime experienced by the fleet for each design. Development of the CMIP model had two major steps: (1) describing modeling assumptions and (2) developing an integer program (IP) that jointly optimizes decisions about locations to install DC fast chargers and EV-to-charger assignments. The CMIP model integrates an EV charging model, EV energy consumption model, and heterogeneous, real-world vehicle use data with an integer programming optimization model to identify optimal location of new charging stations and calculate vehicle downtime for charging. The CMIP model can be applied to understand: (a) the reduction of EV fleet downtime if an additional fast-charging station is added to the current infrastructure and (b) to what extent total vehicle downtime would be sensitive to additional charging infrastructure.

Roni, MohammadS↗

Effect of fluctuating redox conditions on trace metal release from wetland soils from South Carolina and stream sediments from Tennessee

Natural aquatic systems undergo fluctuating redox conditions due to microbial activity, varying water saturation levels, and nutrients dynamics. With fluctuating oxic and anoxic conditions, trace metals can be mobilized or sequestered in response to changes in iron and sulfur speciation and the concentrations and lability of organic carbon. We conducted a systematic laboratory-based microcosm study to examine the effect of redox fluctuations on trace metal mobility in samples collected from two different natural aquatic systems: riparian wetlands and a stream. We incubated water-saturated soils under three cycles of anoxic-oxic conditions (τanoxic:τoxic = 3) spanning 24 days and monitored the change in dissolved and bioavailable metal (copper (Cu), nickel (Ni), zinc (Zn), cobalt (Co), iron (Fe), and manganese (Mn)) concentrations. The dataset includes the results of the variation in metal concentrations due to redox fluctuations in wetland soils and stream sediments. Additionally, the variation in the bioavailable concentration of metals obtained using passive samplers based on diffuse gradient in thin-films technique is included in the dataset. We have also reported the change in the pH values and dissolved oxygen concentrations in the microcosms during the entire study of 24 days. The dataset also contains the dissolved organic carbon and sulfate concentrations, and the data on fractions of metals bound to reducible and oxidizable phases. For materials from both the wetlands and the stream, anoxic conditions favored Co and Zn release, which corresponded with the reductive dissolution of iron oxides. In contrast, dissolved Cu concentrations increased under oxic conditions for both sites and correlated positively with the release of sulfate. In wetland soils, dissolution of Fe (hydr)oxides increased Ni solubility; however, in stream sediments, Ni release occurred when sulfides or organic matter were oxidized. For stream sediments, each subsequent redox cycle increased the bioavailability of trace metals. Upon redox fluctuations in wetland soils, the bioavailability of Zn and Cu increased, whereas bioavailable Ni and Co decreased. This study illustrates that different trace metals display distinct bioavailability patterns during oxic-anoxic fluctuations in natural environments. The biogeochemical cycling of carbon and nitrogen in systems with redox fluctuations may be influenced by these patterns in trace metal availability in addition to the availability of electron donors and acceptors.

54 ENVIRONMENTAL SCIENCES↗

Bridging the length scales on mechanical property evaluation (Final Report)

The development of small-scale mechanical testing in combination with microstructural investigation is of great interest to the nuclear materials community for both materials development and monitoring applications. Dramatically reducing the sample sizes to reduce radioactivity and obtaining mechanical properties of irradiated samples is truly intriguing. Moreover, such studies promise a range of benefits including cost reduction, fundamental insight in structure-property relationships, increased statistics on less sample material, and reinvestigation of prior irradiated and tested reactor samples while simultaneously enabling the generation of mechanical test data on ion beam irradiated materials with limited penetration depths. Small scale materials testing on sub-sized samples has been studied for several decades, though it has only been after the development of micro-testing based on Focused Ion Beam (FIB) sample manufacturing in that orders of magnitude smaller samples could really be investigated in a quantitative manner. In recent years, small scale mechanical testing techniques at a number of length scales has been developed for both unirradiated and irradiated (ion and neutron) materials. Technological advances made in this field have enabled ex-situ and in-situ transmission electron microscopy (TEM) and scanning electron microscopy (SEM) examination, thus leading to more accurate measurements as well as additional mechanistic information. A recent review of the benefits of these techniques show that these techniques are at a stage to tackle multi-scale ranges of materials investigations and can be utilized to obtain fundamental science-based understanding of nuclear materials. Considering the tremendous advances made, one can see how small-scale mechanical testing techniques combined with modeling can enable true small scale to bulk scale mechanical property correlations. However, for the engineering community to adapt this approach fully, one needs to demonstrate that a) that these techniques can produce results with high fidelity and reproducibility, b) generate engineering stress-strain data that one can utilize to understand bulk behavior, and c) generate new insight into relevant phenomena fostering the true understanding of radiation damage and microstructure in materials for nuclear applications. It is the objective of this proposal to bridge the length scale between macro- and micro- scale mechanical testing of unirradiated and irradiated materials. This involves the development and demonstration of procedures for multi-scale mechanical testing that enable high fidelity reproducibility of data and the generation bulk property data from small-scale mechanical tests. Through this, the proposal aims to enhance the confidence in the obtained data at the smaller length scales and enhance the insight provided from these techniques for bulk scale applications on both unirradiated and irradiated materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Transformational Nano-confined Ionic Liquid Membrane for Greater than or Equal to 97 Percent Carbon Dioxide Capture from Natural Gas Combined Cycle Flue Gas

A transformational process based on nano-confined ionic liquid (NCIL) membranes was developed for capturing ≥97% CO 2 from natural gas combined cycle (NCCC) flue gas. The NCIL membranes were prepared by loading amino acid ionic liquid into a framework composed of single-walled carbon nanotube mesh filled with graphene oxide quantum dots. The membranes exhibited CO 2 permeance as high as 2,000 GPU with a CO 2 /N2 selectivity of 2,300 for a typical NGCC flue gas composition. When H 2 O vapor sweep was applied in the permeate side, 96.6% CO2 dry-basis purity and 97.6% CO 2 capture rate were achieved for a simulated NGCC flue gas with single stage. In the process design, a highly H 2 O-selective membrane would be needed to recover majority of the H 2 O vapor, and the recovered H 2 O vapor could be recycled to the permeate side of the NCIL membrane. Sulfonated poly(ether ether ketone) membranes were successfully developed for this purpose. These membranes exhibited H 2 O permeance great than 11,000 GPU and H 2 O/CO 2 selectivity greater than 1,000 at 70ºC for a feed mixture consisting of 14.5 vol% H2O and balanced CO 2 . A standalone membrane model using MATLAB platform was developed for process simulation. The model was validated with experimental data. Techno-economic analysis based on the testing data collected during the current program suggests this transformational membrane process can achieve 97% CO 2 capture efficiency with a cost of $47.8/tonne of CO 2 , which is a 21% reduction versus DOE’s reference case B31B.97.

03 NATURAL GAS↗

New Insights Into Lignification via Network and Multi-Omics Analyses of Arogenate Dehydratase Knock-Out Mutants in Arabidopsis thaliana

Multiple Arabidopsis arogenate dehydratase (ADT) knock-out (KO) mutants, with phenotypes having variable lignin levels (up to circa 70% reduction), were studied to investigate how differential reductions in ADTs perturb its overall plant systems biology. Integrated "omics" analyses (metabolome, transcriptome, and proteome) of wild type (WT), single and multiple ADT KO lines were conducted. Transcriptome and proteome data were collapsed into gene ortholog (GO) data, with this allowing for enzymatic reaction and metabolome cross-comparisons to uncover dominant or likely metabolic biosynthesis reactions affected. Network analysis of enzymes-highly correlated to stem lignin levels-deduced the involvement of novel putative lignin related proteins or processes. These included those associated with ribosomes, the spliceosome, mRNA transport, aminoacyl tRNA biosynthesis, and phosphorylation. While prior work helped explain lignin biosynthesis regulation at the transcriptional level, our data here provide support for a new hypothesis that there are additional post-transcriptional and translational level processes that need to be considered. These findings are anticipated to lead to development of more accurate depictions of lignin/phenylpropanoid biosynthesis models in situ, with new protein targets identified for further biochemical analysis and/or plant bioengineering. Additionally, using KEGG defined functional categorization of proteomics and transcriptomics analyses, we detected significant changes to glucosinolate, a-linolenic acid, nitrogen, carotenoid, aromatic amino acid, phenylpropanoid, and photosynthesis-related metabolic pathways in ADT KO mutants. Metabolomics results also revealed that putative carotenoid and galactolipid levels were generally increased in amount, whereas many glucosinolates and phenylpropanoids (including flavonoids and lignans) were decreased in the KO mutants.

Arabidopsis, arogenate dehydratases, lignin, multi↗

“Intelligent” Pt Catalysts Based on Thin LaCoO3 Films Prepared by Atomic Layer Deposition

LaCoO3 films were deposited onto MgAl2O4 powders by atomic layer deposition (ALD) and then used as catalyst supports for Pt. X-ray diffraction (XRD) showed that the 0.5 nm films exhibited a perovskite structure after redox cycling at 1073 K, and scanning transmission electron microscopy and elemental mapping via energy-dispersive X-ray spectroscopy (STEM/EDS) data demonstrated that the films covered the substrate uniformly. Catalysts prepared with 3 wt % Pt showed that the Pt remained well dispersed on the perovskite film, even after repeated oxidations and reductions at 1073 K. Despite the high Pt dispersion, CO adsorption at room temperature was negligible. Compared with conventional Pt on MgAl2O4, the reduced forms of the LaCoO3-containing catalyst were highly active for the CO oxidation and water gas shift (WGS) reactions, while the oxidized catalysts showed much lower activities. Surprisingly, the reduced catalysts were much less active than the oxidized catalysts for toluene hydrogen. Catalysts prepared from thin films of Co3O4 or La2O3 exhibited properties more similar to Pt/MgAl2O4. Possible reasons for how LaCoO3 affects properties are discussed.

42 ENGINEERING↗

Bias correcting regional scale Earth system model projections: novel approach using empirical mode decomposition

Bias correction is a crucial step in using Earth system model outputs for assessments, as it adjusts systematic errors by comparing the model to observations. However, standard methods – ranging from mean-based linear scaling to distribution-based quantile mapping typically treat bias correction as a single-scale process, overlooking the fact that biases can manifest differently across daily, seasonal, and annual timescales. In this study, we propose a novel, timescale-aware bias-correction approach built on Empirical Mode Decomposition. By decomposing the meteorological signal into multiple oscillatory components and aggregating them to represent distinct timescales, we apply targeted corrections to each component, thereby preserving both short- and long-term structure in the data. Experimental illustrations show that the timescale-aware EMDBC framework matches the performance of conventional quantile-delta mapping (QDM) at the native daily scale and achieves progressively larger bias reductions at bi-weekly, seasonal, and annual scales. As a result, the proposed approach offers a more robust path to accurate and reliable Earth system projections, strengthening their utility for resilience and adaptation planning.

Ganguli, Arkaprabha [Argonne National Laboratory (↗

Observations and simulations of a wind farm modifying a thunderstorm outflow boundary

Abstract. On 18 June 2019, National Weather Service (NWS) radar reflectivity data indicated the presence of thunderstorm-generated outflow propagating east-southeastward near Lubbock, Texas. A section of the outflow boundary encountered a wind farm and then experienced a notable reduction in ground-relative velocity, suggesting that interactions with the wind farm impacted the outflow boundary progression. We use the Weather Research and Forecasting model and its wind farm parameterization to address the extent to which wind farms can modify the near-surface environment of thunderstorm outflow boundaries. We conduct two simulations of the June 2019 outflow event: one containing the wind farm and one without. We specifically investigate the outflow speed of the section of the boundary that encounters the wind farm and the associated impacts on near-surface wind speed, moisture, temperature, and changes to precipitation features as the storm and associated outflow pass over the wind farm domain. The NWS radar and nearby West Texas Mesonet surface stations provide observations for validation of the simulations. The presence of the wind farm in the simulation clearly slows the progress of the outflow boundary by over 20 km h -1 , similar to what was observed. Simulated perturbations of surface wind speed, temperature, and moisture associated with outflow passage were delayed by up to 6 min when the wind farm was present in the simulation compared to the simulation without the wind farm. However, impacts on precipitation were localized and transient, with no change to total accumulation across the domain.

17 WIND ENERGY↗

Performance Enhancement of the Transcritical CO2 Cycle with a Near-Isothermal Liquid Piston Compression by an Ejector

This study investigates the use of an ejector in a transcritical CO2 cycle employing a near-isothermal liquid piston compressor, with the objective of mitigating refrigerant degassing in mineral oil. While the liquid piston compressor demonstrated high isothermal efficiency (˜90%), overall cycle performance was limited by CO2 degassing during chamber pressure fluctuations. A thermodynamic model was developed to compare baseline and ejector-assisted cycles, incorporating solubility correlations, entrainment ratio characteristics, and a degassing correction factor. Results indicate that the ejector increased suction pressure and reduced CO2–oil solubility variation, which improved refrigerant intake and reduced degassing losses. Despite these benefits, the ejector-assisted cycle showed a lower coefficient of performance (COP) than the baseline. Possible contributing factors include the limited marginal benefit of pressure reduction in a highly efficient isothermal compressor, reduced effectiveness of the separator due to the use of a suction-line heat exchanger in both cycles, and uncertainties in the solubility data under transient conditions. These findings suggest that degassing mitigation is essential for realizing the potential of liquid piston compression, but ejector integration alone is unlikely to enhance COP. Future work should focus on ejector optimization, improved oil–refrigerant management, introducing less-soluble working liquids, and experimental validation under realistic operating conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

ReEDS Performance Improvement

The Regional Energy Deployment System (ReEDS) is an open-source, spatially explicit, long-term capacity expansion model for the bulk electric power system of the contiguous United States, encompassing multiple scenarios with technological and political assumptions (see https://github.com/NREL/ReEDS-2.0). With the increased needs for capabilities, higher temporal and spatial resolutions to model the evolution of the power system with modern technologies and low-carbon pathways, ReEDS' model solution times have increased significantly from 4-6 hours in 2018 to 18-48+ hours in 2023 . Also, the model size for commonly-run ReEDS scenarios reached 22 and 28 million equations and variables, respectively. These runtimes can be especially challenging under certain scenario settings (e.g., very high temporal or spatial resolution) or with limited computational power. In this presentation, we will discuss several methods we used to improve model runtime, including data preparation, model modification, and solver tuning. The implementation of these methods shrank the model size to 7.2 and 7.3 million equations and variables, respectively. Furthermore, this led to a 77% reduction in the model's run time for commonly-run ReEDS scenarios. We will discuss the process of identifying areas for solve time improvements and how the specific enhancements for the ReEDS model might be applied to other similar large-scale models.

ENERGY PLANNING, POLICY, AND ECONOMY,MATHEMATICS A↗

Throughfall-Reduction Drying Effects on Soil Respiration, Moisture, and Temperature in Four Distinct Lowland Panamanian Forests from 2015-2022

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that throughfall reduction (i.e. chronic drying) would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We conducted a field study in four distinct lowland Panamanian forests that vary in background rainfall and soil fertility. In each forest we constructed throughfall reduction structures that reduce moisture inputs to plots by ~50% (see methods). In drying and paired control plots we measured soil CO2 fluxes, soil moisture (volume/volume), soil temperature (°C), air temperature (°C), and forest floor biomass during wet-dry cycles (2015 – 2022). We also collected litterfall biomass and surveyed all tree species for diameter at breast height (dbh), and identified tree species in 2018 and 2019.Results:We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks. Datasets included: Datasets included here include .csv and .xls files for soil respiration (carbon dioxide flux) collected using the LiCor (see methods), air temperature °C , soil moisture (vl/vl), and soil temperature °C. Soil moisture and soil temperature were collected at 0-10cm depths using hand-held probes (see methods). Data are on ~3 week timescales from 2015-2022 with some gaps. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗