Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initialization”

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 109 records · Page 6

The Global Spectra-Trait Initiative: A database of paired leaf spectroscopy and functional traits associated with leaf photosynthetic capacity

Accurate assessment of leaf functional traits is crucial for a diverse range of applications from crop phenotyping to parameterizing global climate models. Leaf reflectance spectroscopy offers a promising avenue to advance ecological and agricultural research by complementing traditional, time-consuming gas exchange measurements. However, the development of robust hyperspectral models for predicting leaf photosynthetic capacity and associated traits from reflectance data has been hindered by limited data availability across species and environments. Here we introduce the Global Spectra-Trait Initiative (GSTI), a collaborative repository of paired leaf hyperspectral and gas exchange measurements from diverse ecosystems. The GSTI repository currently encompasses over 7500 observations from 397 species and 41 sites gathered from 36 published and unpublished studies, thereby offering a key resource for developing and validating hyperspectral models of leaf photosynthetic capacity. The GSTI database is developed on GitHub (https://github.com/plantphys/gsti, last access: 4 January 2026) and published to ESS-DIVE https://doi.org/10.15485/2530733, Lamour et al., 2025). It includes gas exchange data, derived photosynthetic parameters, and key leaf traits often associated with traditional gas exchange measurements such as leaf mass per area and leaf elemental composition. By providing a standardized repository for data sharing and analysis, we present a critical step towards creating hyperspectral models for predicting photosynthetic traits and associated leaf traits for terrestrial plants.

Lamour, Julien [Université of Toulouse (France); U↗

Initial Engineering Framework for Species Tracking in MSRs Involving Fuel Salt and Structures

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track leaching and plating of materials between the fuel salt and structures in Molten Salt Reactors (MSRs). Specifically, this framework is demonstrated in modeling various processes within MSRs including hot and cold leg corrosion, redox potential driven corrosion, colloid precipitation and deposition, and noble metal plating and decay heat modeling. The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the evolving chemical behavior of the fuel salt in interactions with structural materials. Initial findings suggest that hot-cold leg corrosion in MSR systems may become significant if the redox potential of the fuel salt is not sufficiently controlled. Additionally, the impact of noble metal decay heat deposition in large power reactors is explored. The resulting engineering framework for species tracking will be applied in the future toward specific validation work with ongoing experimental efforts focusing on both chloride and fluoride salt loop experiments with in-situ corrosion and redox control instrumentation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Positron emission tomography harmonization in the Alzheimer's Disease Neuroimaging Initiative: A scalable and rigorous approach to multisite amyloid and tau quantification

Abstract INTRODUCTION A key goal of the Alzheimer's Disease NeuroImaging Initiative (ADNI) positron emission tomography (PET) Core is to harmonize quantification of β‐amyloid (Aβ) and tau PET image data across multiple scanners and tracers. METHODS We developed an analysis pipeline (Berkeley PET Imaging Pipeline, B‐PIP) for ADNI Aβ and tau PET images and applied it to PET data from other multisite studies. Steps include image pre‐processing, refacing, magnetic resonance imaging (MRI)/PET co‐registration, visual quality control (QC), quantification of tracer uptake, and standardization of Aβ and tau standardized uptake value ratios (SUVrs) across tracers. RESULTS Measurements from 10,105 cross‐sectional and longitudinal Aβ and tau PET scans acquired in several studies between 2010 and 2024 can be processed, harmonized, and directly merged across tracers and cohorts. DISCUSSION The B‐PIP developed in ADNI is a scalable image harmonization approach used in several observational studies and clinical trials that facilitates rigorous Aβ and tau PET quantification and data sharing. Highlights Quantitative results from ADNI Aβ and tau PET data are generated using a rigorous, scalable image processing pipeline This pipeline has been applied to PET data from several other large, multisite studies and trials Quantitative outcomes are harmonizable across studies and are shared with the scientific community

Neurosciences & Neurology↗

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗

Initiation of dusting corrosion in high-temperature alloys under CO exposure

Carbon monoxide is commonly encountered in energy systems, yet its reactivity with structural alloys—critical heat-resistant components in these systems—has been largely overlooked compared to the well-documented effects of oxidizing gases. In contrast, we demonstrate the high-temperature reaction of CO with NiAl using in-situ low-energy electron microscopy and X-ray photoemission electron microscopy. Our results show that CO dissociates into atomic oxygen and carbon, resulting in two concurrent reactions: selective oxidation of aluminum to form Al 2 O 3 and the initiation of dusting corrosion through carbon dissolution into the alloy and subsequent carbon deposition on the surface. These reactions produce spatially distinct surface products, preventing the formation of a continuous protective Al oxide layer. These results reveal a preference for the dissociative pathway of CO over the classic Boudouard disproportionation reaction that forms CO 2 . These insights not only advance our understanding of CO-induced alloy degradation but also highlight the practical implications for managing alloy stability and optimizing catalysis in carbon-rich environments, such as those in petrochemical processing and hydrocarbon combustion.

36 MATERIALS SCIENCE↗

Computational analysis of flame initiation, quenching, and re-ignition in a prechamber natural gas engine under varying EGR-dilution levels

The on-road natural-gas (NG) fueled transportation relies on stoichiometric spark-ignition engines for the advantages of simple after-treatment system despite the efficiency penalty relative to lean-burn combustion strategies. Exhaust gas recirculation (EGR) has the potential to reduce this efficiency gap at low to moderate loads without the need for complex lean-exhaust aftertreatment systems. However, EGR dilution leads to reduced combustion stability and increased cycle-to-cycle variability. A promising technology that has the potential to achieve reliable operation under diluted conditions is the prechamber ignition (or turbulent jet ignition) which uses chemically active turbulent jets generated from combustion inside a prechamber to initiate, stabilize and accelerate combustion of the mixture inside the main chamber. The present work focusses on developing a RANS-based CFD approach to accurately reproduce in-cylinder phenomena in a stoichiometric NG prechamber-assisted heavy-duty engine without relying on complex combustion models that account for turbulence-chemistry interactions. This is necessary because reactive prechamber jets at high EGR dilution tend to extinguish while emerging into the main chamber, which is followed by a phase of re-ignition — a phenomenon that conventional G-equation or well-stirred reactor combustion models cannot reproduce. With addition of a damping multiplier to the well-stirred reactor model, the predictions are seen to show good agreement with experimental pressure evolution and combustion images acquired from a single cylinder Cummins N-14 optical diesel engine retrofitted with a prechamber ignition system. Model predictions of local heat release in the flame and temperature evolution inside the flame are used to investigate combustion dynamics in the prechamber and the main chamber. It is seen that the well-stirred reactor model with the inclusion of damping is able to reproduce the temporary reduction in heat release within the flame, which can be considered equivalent to quenching of jets, and the subsequent re-ignition of the flame inside the main chamber. The delay between quenching and re-ignition depends on the amount of dilution, as explained by an illustration of flame evolution in a Borghi diagram.

Prechamber ignition↗

Influence of initial conditions on data-driven model identification and information entropy for ideal mhd problems

Data-driven methods of model identification are able to discern governing dynamics of a system from data. Such methods are well suited to help us learn about systems with unpredictable evolution or systems with ambiguous governing dynamics given our current understanding. Many plasma problems of interest fall into these categories as there are a wide range of models that exist, however each model is only useful in a certain regime and often limited by computational complexity. To ensure data-driven methods align with theory, they must be consistent and predictable when acting on data whose governing dynamics are known. Weak Sparse Identification of Nonlinear Dynamics (WSINDy) is a recently developed data-driven method that has shown promise in learning governing dynamics from data with high noise levels [1]. This work examines how WSINDy acts on ideal MHD test problems as the initial conditions are varied and specifies limiting requirements for successful equation identification. Furthermore, it is hard to recover the governing dynamics from data that emphasize a single dominant behavior. In these low information cases, Shannon information entropy is able to pick up on the redundancies in the data that affect recoverability.

97 MATHEMATICS AND COMPUTING↗

Dual roles of stepped surfaces: Catalytic initiators and stabilizers of oxygen-induced reconstructions

Understanding surface restructuring under reactive conditions is crucial for designing next-generation catalysts with enhanced activity and selectivity. Here, we employ in situ transmission electron microscopy to directly observe the dynamic behavior of Cu(100) and Cu(410) surfaces under both oxidizing and vacuum annealing conditions, revealing a complex interplay among surface crystallography, local oxygen coverage, Cu atom mobility, and step-edge reactivity. The stepped Cu(410) surface acts as an active site for O 2 dissociation, triggering the oscillatory transformation of the c(2 × 2)–O phase into the more stable (2$\sqrt2$ ×$\sqrt2$)R45°–O missing-row (MR) structure on the adjacent flat Cu(100) terrace. Under subsequent vacuum annealing, this same Cu(410) facet exhibits remarkable structural resilience, preserving the MR reconstruction and chemisorbed oxygen. In contrast, the Cu(100) surface undergoes reversible transitions from the MR structure back to the c(2 × 2)–O phase. These results highlight the critical role of surface morphology in directing both the formation and stability of oxygen-induced reconstructions, demonstrating that stepped surfaces serve dual roles as both catalytic initiators and structural stabilizers. Furthermore, this work offers atomic-level insights into the environment-responsive behavior of copper surfaces, establishing a mechanistic basis for designing Cu-based catalysts through facet-specific control of surface reactivity.

36 MATERIALS SCIENCE↗

Initial design concepts for solid boron injection in ITER

As part of ITER’s consideration to change its first wall material from beryllium to tungsten, the ITER Organization has proposed studying the feasibility of real-time solid boron injection (SBI) into the plasma to coat the walls and divertor to supplement glow discharge boronization (GDB). Boron deposits getter oxygen and reduce sputtering of tungsten from plasma facing components (PFCs). Particularly in areas with significant plasma wall interactions, boron coatings are expected to be short-lived under high performance plasma conditions. The proposed SBI system aims to maintain boron layers in these areas to avoid excessive radiation from tungsten in the plasma as a risk mitigation to ensure ITER will be able to reach and sustain Q = 10 conditions. The system will be used sparingly, as redeposition of boron can lead to significant tritium retention, which must be minimized in ITER to comply with nuclear safety concerns. SBI is proposed to limit and precisely control the amount of boron injected in real-time during plasma operation. Here, some of the design requirements and initial concepts for an SBI system in ITER are presented based on previous results carried out with SBI systems on a number of tokamaks and stellarators around the world. Previous results using SBI systems installed by PPPL have injected boron particles from 5 µm–2 mm diameter at calibrated rates of 2–200 mg/s in real-time during plasma operation on AUG, DIII-D, EAST, KSTAR, LHD, TFTR, WEST, and W7-X, leading to improved wall conditions with reduced plasma impurity concentrations and radiated power and improved plasma performance. The boron is ionized in the plasma edge and then deposited on plasma-wetted surfaces. On AUG, EAST and WEST reduced tungsten sputtering sources were observed following several discharges with SBI. Extrapolation of these SBI results are presented to estimate the amount of boron needed for wall conditioning in ITER. Real-time SBI control requirements and plasma operation scenarios for ITER are also described.

36 MATERIALS SCIENCE↗

Impact of recent ENDF nuclear data update, high initial enrichment and high burnup fuel on critical experiments applicability determination via the integral index c k for burnup credit validation

In 2012, NUREG/CR-7109 reported on the validation of burnup credit calculations involving major and minor actinides and major fission products which was investigated for pressurized and boiling water reactor (PWR and BWR) fuel enrichments up to 5 wt% 235 U and assembly-average burnups up to 60 GWd/MTU. Recently, there has been interest in increasing the maximum enrichment used in PWR fuel as high as 8 wt% 235 U and correspondingly increasing the maximum assembly-average burnups to approximately 75 GWd/MTU. These proposed increases in enrichment and burnup necessitate reinvestigation of the validation basis for k eff calculations for this expanded application space. Additionally, the 2012 study was performed by using the Evaluated Nuclear Data File (ENDF)/B-VII.0 nuclear data with the SCALE 6 covariance library, and the effects of using the newly released ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data and covariance libraries should be evaluated. In this work, published in NUREG/CR-7309 in 2025, the validation assessment was performed consistently with NUREG/CR-7109: modeling irradiated fuel assemblies in the Generic Burnup Credit (GBC)-32 cask defined in NUREG/CR-6747. The TSUNAMI-3D sequence was used to generate sensitivity data for the application model, and the data were compared with sensitivity data from select benchmark models. The integral parameter c k is the metric of similarity used in this study and is consistent with NUREG/CR-7109, where a c k value in excess of 0.8 indicates sufficient similarity for use in validation. A new set of benchmark experiments with sensitivity data has been assembled for this effort. The number of experiments with available sensitivity data is now 2,104, compared to 474 in NUREG/CR-7109. This increase was facilitated by the efforts of the Nuclear Energy Agency to generate sensitivity data for a majority of the experiments in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook to supplement the data available in the Oak Ridge National Laboratory (ORNL) Verified, Archived Library of Inputs and Data (VALID). The complete set of benchmarks considered here includes experiments for low-enriched uranium (LEU), intermediate enriched uranium (IEU), and a mixture of uranium and plutonium (MIX) from the ICSBEP Handbook and VALID, as well as ORNL models of the Haut Taux de Combustion (HTC) experiments and other potentially relevant models not included in VALID. The updated similarity study shows that none of the extended burnup and higher enrichment combinations considered show a significant decrease in the number of potentially applicable experiments, meaning sufficient critical experiments exist for the validation of BUC criticality safety calculations, with initial enrichments up to 8 wt% 235 U and burnups up to 80 GWd/MTU. Additionally, both the ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data libraries can be used for validation since the number of critical experiments applicable for validation increases for most cases with the most recent nuclear data compared to the previous one. As in previous BUC validation studies, the French HTC experiments are the most similar in a majority of the application cases studied, especially from representative discharge burnups ranging from 40 to 80 GWd/MTU. In conclusion, these results match the conclusions presented in NUREG/CR-7109 regarding validation of the primary actinides in BUC analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

How can an ecosystem approach support integrated management of marine renewable energy? An initial assessment from an environmental point of view

With the increasing installation of marine renewable energy (MRE) devices in areas already subject to multiple anthropogenic activities and environmental changes, it is necessary to develop tools and methods for the integrated management of marine ecosystems. The ecosystem approach is a holistic environmental management method that considers all components of an ecosystem. The ecosystem approach has demonstrated utility in the application to various anthropogenic activities and is relevant for consideration within the context of MRE. Indeed, many of the effects observed on marine ecosystems from those other activities are also applicable to MRE development. This review is an initial assessment where we summarize the potential effects of MRE development on marine ecosystems and propose schematic frameworks for applying the ecosystem approach to MRE. We also provide a non-exhaustive list of commonly used models pertinent to the ecosystem approach and associated with several reference studies. An outline of core questions that can currently be answered using available modeling tools central to the ecosystem approach is provided, along with recommendations for the application of this approach to the MRE context. Further, we identify key knowledge gaps and areas that require additional investigation for meaningful application of the ecosystem approach to MRE development. Our recommendations mainly concern the current limitations of applying the ecosystem approach to concrete cases, such as consolidating knowledge of the effects of MRE on the local environment, the need to obtain fine-scale data, considering effects at different spatiotemporal scales, and, finally, the need for an interdisciplinary vision.

16 TIDAL AND WAVE POWER↗

Insights from Initial Engineering Designs of Point Source Capture at Industrial Facilities

Initial engineering design studies examining the application of state-of-the-art carbon capture technology at industrial plants contain generally overlooked real-world design considerations for near-term deployment of point source capture (PSC). The implementation of PSC across a wide range of industrial applications presents unique challenges associated with fluctuating CO<sub>2</sub> concentrations, flue gas composition, and utility and land availability. In this article, seven recent industrial retrofit PSC projects are reviewed to investigate the impact of site-specific factors on project design and cost. Across the seven projects, three capture technology classes and four industrial applications are considered, allowing insight into industry-specific opportunities for PSC technology synergy. Common challenges across projects and proposed design solutions are highlighted to propagate ideas and solutions to close the technology gaps and accelerate learning rates.

FEED Studies↗

Investigation of O( 3 P) Initiated Oxidation Products of 2,3-Dimethylfuran Using Synchrotron Photoionization

Biofuels are an appealing alternative to augment conventional fossil fuels, given that traditional fossil fuels are finite in nature despite being the highest-consumed energy source. Furans have been studied as potential biofuels due to their renewability, high energy density, and ability to both blend with existing fossil fuels and serve as a stand-alone fuel. Although promising, more research efforts on the combustion of furan-based biofuels are needed. Here, in this study, the oxidation of 2,3-dimethylfuran (2,3-DMF) initiated by O( 3 P) was investigated using vacuum-ultraviolet synchrotron radiation from the Advanced Light Source at the Lawrence Berkeley National Laboratory. The reaction was studied at 550 K and 7 Torr. Major reaction products include 3-methyl-3-buten-2-one, 4,5-dimethyl-3(2H)-furanone, and 3-methyl-4-oxo-2(Z)-pentenal, with 3-methyl-3-buten-2-one being the most abundant.

CBS-QB3↗

Initial Laboratory Measurements Probing Hydrogen Interactions with Eagle Ford Shale and Pyrite: Potential Implications for Subsurface Hydrogen Storage

Hydrogen (H 2 ) has the potential to be a transformative technology as an enabler to a low-carbon future and promoter of renewable energy. When H 2 is injected and stored in the subsurface, it has the potential to interact with the caprock (usually shale) which overlies and seals the storage reservoir. Here, this study examines geochemical reactions or changes in surface morphology to Eagle Ford Shale, a proxy for caprock, upon exposure to H 2 at 50°C and 10.3 MPa. Reactions were also performed with N 2 to provide an experimental control. Scanning Electron Microscopy with Energy Dispersive Spectroscopy (SEM-EDS), Atomic Force Microscopy (AFM), and Optical Photothermal-Infrared (O-PTIR) Spectroscopy were applied to quantify changes on the microscale and nanoscale level. Fluid chemistry changes were monitored with Ion Chromatography (IC) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Exposure of Eagle Ford Shale to H 2 gas alone did not result in any alterations to the shale chemically or any changes in the surface morphology. Exposure of Eagle Ford Shale to both H 2 and water as well as N 2 and water resulted in changes to the surface morphology because of gypsum dissolution and reprecipitation, thus indicating that H 2 is not necessary to promote changes. Pure pyrite was the most reactive with H 2 possibly resulting in a reduction to pyrrhotite. These initial studies suggest that the extent of reactions activated by hydrogen with caprock are minor under the temperature and pressure conditions that would represent underground hydrogen storage.

08 HYDROGEN↗

Thermally Initiated Formation of Criegee Intermediate CH 2 OO in the Oxidation of Ethane

Criegee intermediates (CIs) play an important role in atmospheric chemistry as a transient source of the OH radical through their formation by the ozonolysis of unsaturated organic compounds. Here, we report thermally initiated formation of the smallest CI (CH 2 OO) in the oxidation of ethane (CH 3 CH 3 ) that may be relevant to combustion and flames. The SiO 2 /SiC oxidation microreactor is heated to 1800 K and has a short residence time of ∼100 μs. The CH 2 OO we observe is likely formed in a lower-temperature region near the microreactor’s exit. Plausible mechanisms for CH 2 OO formation and retention under these conditions mediated by methylperoxy (CH 3 OO) radicals are discussed. Pure rotational spectra of CH 2 OO and other intermediates (HO 2 , CH 3 CHO, CH 2 CHOH, c-CH 2 OCH 2 , CH 3 CH 2 CHO, CH 3 OOH, and HCOOH) are detected with a chirped-pulse Fourier transform millimeter-wave spectrometer operating in the frequency range of 60–90 GHz. In conclusion, detection occurs in a molecular beam, where the species are supersonically cooled to 5 K.

Hydrocarbons↗

Operando Infrared Nanospectroscopy of the Silicon/Electrolyte Interface during Initial Stages of Solid-Electrolyte-Interphase Layer Formation

The solid electrolyte interphase (SEI) is a critical component in Li-ion batteries; however, its nanoscale structure and composition and unstable nature make it difficult to characterize and ascertain primary functional mechanisms. We use operando nanoscale Fourier transform infrared spectroscopy (nano-FTIR) with a broadband synchrotron IR source to study the SEI formation on a thin-film Si electrode at nanometer-scale spatial resolution as a function of time and voltage. By probing the Si/carbonate electrolyte interface through a 25 nm-thick amorphous Si window/electrode, we detect molecular vibrational modes within a 10s of nanometers region adjacent to the Si surface and observe that PF6 – anions react to form LiF at 0.5 V. Spatially resolved nano-FTIR spectra showcase subtle nanoscale heterogeneities in the initial solid/liquid interface and the resulting deposited LiF. With its nanoscale resolution and high chemical specificity, operando nano-FTIR provides unique insights into the dynamics and heterogeneous formation of SEIs and opens opportunities for connecting nanoscale interfacial properties to bulk performance metrics.

Dopilka, Andrew↗

Electrochemically Initiated Depolymerization of Poly(Methyl Methacrylate)

Efficient depolymerization of polymers with all-carbon backbones under mild conditions would be valuable in chemically recycling commodity plastics. Poly(methyl methacrylate) (PMMA) is a commodity thermoplastic that is currently depolymerized under temperatures in excess of 400 °C. Herein, we lower the temperatures needed to achieve depolymerization of PMMA by performing radical generation and depropagation with orthogonal stimuli. This first demonstration of electrochemically initiated PMMA depolymerization relies on reduction of phthalimide esters that, upon subsequent decarboxylation, generate polymer-centered radicals. These radicals then spontaneously unzip the polymer back to its monomeric constituents at temperatures as low as 105 °C. We studied the mechanism and efficiency of this transformation as a function of phthalimide ester placement, incorporation density, and polymer molecular weight. We found that chain-end activation is effective for modest molecular weights but suffers diminished efficiency at higher degrees of polymerization. In contrast, pendent-group activation is more effective for depolymerizing higher molecular weight species. Integrating higher molar amounts of phthalimide ester pendants leads to more effective depolymerization, with >95% depolymerization in copolymers with 5 mol% phthalimide ester incorporation. We leveraged this understanding to create a custom electro-distillation apparatus that allowed us to simultaneously electrochemically depolymerize PMMA and directly distill methyl methacrylate in >22% yield, which could ultimately be repolymerized. These findings establish electrochemistry as a versatile and orthogonal stimulus for vinyl polymer depolymerization and provide a foundation for closed-loop electrochemical recycling of widely used plastics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

fluxfinder: An R Package for Reproducible Calculation and Initial Processing of Greenhouse Gas Fluxes From Static Chamber Measurements

Fluxes of greenhouse gases are a critical component of the earth's natural climate, but anthropogenic emissions have created an imbalance and resulted in global climate change. Quantifying the emission of these gases is vital to our understanding of their sources and sinks, both natural and anthropogenic. The static chamber method, in which a system of interest is enclosed, and gas concentrations are measured over time, is widely used to estimate fluxes of greenhouse gases. With the development of instruments such as infrared gas analyzers (IRGAs) supporting high-frequency concentration data, there is a growing need for open-source workflows to calculate fluxes. Here we present fluxfinder, an R package designed to support reproducible calculations and processing of greenhouse gas fluxes measured with the static chamber method. The package includes raw data file parsing from widely used IRGAs, metadata matching, unit conversion, flux estimations, and initial quality assurance/quality control (QA/QC). Diagnostic graphical plots provide a transparent way to differentiate between measurement issues and nonlinear behavior. The package is also designed to be easily integrated with the gasfluxes package for further fitting of nonlinear concentration-time models, allowing alternative or additional flux QA/QC. The fluxfinder package offers a flexible workflow that is easily adaptable to promote open and reproducible greenhouse gas flux estimations.

Wilson, Stephanie J.↗