Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “2018 data”

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,117 records · Page 62

SECARB Post Injection Assessment Report

The Southern States Energy Board (SSEB) managed the overall Southeast Regional Carbon Sequestration Partnership (SECARB) Phase III project and The Bureau of Economic Geology (BEG), at the University of Texas at Austin, managed all activities associated with the Early Test field site near Cranfield, Mississippi (the Early Test or Cranfield Project). The SECARB project at Cranfield was focused downdip of the oil rim and on demonstrating the stacked nature of the storage potential at this location: enhanced oil recovery (EOR) up-dip and saline downdip. Because the injection for EOR continued after the SECARB project ended, no actual post-injection period occurred. However, information collected during the SECARB-monitored injection period, including pressure and CO2 plume position, was utilized in an EOR carbon lifecycle assessment. The carbon lifecycle assessment provided insight on tracking injection fronts and production efficiencies across the field. To complete the post injection assessment, taking into account the continuing EOR period after closeout of SECARB field activities, this report provides data and analysis done under another DOE-funded project. The follow up project focused on the carbon lifecycle of EOR and included an assessment of EOR at as-operated and under possible future scenarios. This work was published by Hossieni et al. (2018) and Hosseininoosheri (2019).

01 COAL, LIGNITE, AND PEAT↗

Real Time Data Analysis With the ATLAS Trigger at the LHC in Run-2

The trigger selection capabilities of the ATLAS detector have been significantly enhanced for the LHC Run- 2 in order to cope with the higher event rates and with the large number of simultaneous interactions (pile-up) per protonproton bunch crossing. A new hardware system, designed to analyse real time event-topologies at Level-1 came to full use in 2017. A hardware-based track reconstruction system, expected to be used real-time in 2018, is designed to provide track information to the high-level software trigger at its full input rate. The high-level trigger selections are largely relying on offline-like reconstruction techniques, and in some cases multivariate analysis methods. Despite the sudden change in LHC operations during the second half of 2017, which caused an increase in pile-up and therefore also in CPU usage of the trigger algorithms, the set of triggers (so called trigger menu) running online has undergone only minor modifications thanks to the robustness and redundancy of the trigger system, and the use of a levelling luminosity scheme in agreement with LHC and other experiments. Here, this presentation gives a brief yet comprehensive review of the real-time performance of the ATLAS trigger system in 2017. Considerations will be presented on the most relevant parameters of the trigger (efficiency to collect signal and output data rate) as well as details on some aspects of the algorithms which are run real-time on the High Level Trigger CPU farm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-Fidelity Analysis of ULF Wave Mode Structure Following Interplanetary Shock Compression of the Dayside Magnetopause Using MMS Multi-Point Observations

During the 9 March 2018 event with two consecutive interplanetary shocks compressing the dayside magnetosphere, the azimuthal mode structure and frequency spectrum of ultra low frequency magnetic pulsations are resolved using a cross-spectral analysis based on high-fidelity multi-probe Magnetospheric Multiscale Mission (MMS) magnetometer data. The results based on the MMS 4 and MMS 3 pair of measurements show that shock arrival leads to low mode (|m| ≤ 3) magnetic fluctuations in the Pc4-5 regimes, and smaller spatial scale fluctuations implied by the dominant high mode numbers are observed after both shock signatures hit and passed the magnetosphere. Detailed evolution of the mode structure is also shown for the first shock to reveal the development of high mode structure from a bump-on-tail distribution at m≈20 to a dominant peak at m≈50 in about 10 min. In addition, an interesting change of sign in m from negative to positive is observed as MMS crosses ~11 MLT pre-noon, which is consistent with the picture of wave generation by dayside magnetopause compression and then anti-sunward propagation. Finally, for both shocks, the contribution of higher frequency waves (Pc-4 range compared with Pc-5) to the total wave power is found to be negligible before and after the shock impact, but it becomes more significant during the shock impact.

79 ASTRONOMY AND ASTROPHYSICS↗

BeyondPlanck: XII. Cosmological parameter constraints with end-to-end error propagation

We present cosmological parameter constraints estimated using the Bayesian BEYONDPLANCK analysis framework. This method supports seamless end-to-end error propagation from raw time-ordered data onto final cosmological parameters. As a first demonstration of the method, we analyzed time-ordered Planck LFI observations, combined with selected external data (WMAP 33–61 GHz, Planck HFI DR4 353 and 857 GHz, and Haslam 408 MHz) in the form of pixelized maps that are used to break critical astrophysical degeneracies. Overall, all the results are generally in good agreement with previously reported values from Planck 2018 and WMAP, with the largest relative difference for any parameter amounting about 1σ when considering only temperature multipoles between 30 ≤ ℓ ≤ 600. In cases where there are differences, we note that the BEYONDPLANCK results are generally slightly closer to the high-ℓ HFI-dominated Planck 2018 results than previous analyses, suggesting slightly less tension between low and high multipoles. Using low-ℓ polarization information from LFI and WMAP, we find a best-fit value of τ = 0.066 ± 0.013, which is higher than the low value of τ = 0.052 ± 0.008 derived from Planck 2018 and slightly lower than the value of 0.069 ± 0.011 derived from the joint analysis of official LFI and WMAP products. Most importantly, however, we find that the uncertainty derived in the BEYONDPLANCK processing is about 30 % greater than when analyzing the official products, after taking into account the different sky coverage. We argue that this uncertainty is due to a marginalization over a more complete model of instrumental and astrophysical parameters, which results in more reliable and more rigorously defined uncertainties. We find that about 2000 Monte Carlo samples are required to achieve a robust convergence for a low-resolution cosmic microwave background (CMB) covariance matrix with 225 independent modes, and producing these samples takes about eight weeks on a modest computing cluster with 256 cores.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the Frequency of Extreme Precipitation Events at Los Alamos National Laboratory

The Environmental Protection and Compliance (EPC) Storm Water Permitting and Compliance Program is currently updating its stormwater compliance standards for the Los Alamos National Laboratory (LANL) Engineering Standards Manual (2015) and is addressing the required current meteorological data needed for storm water applications. Information on the frequency of extreme precipitation events at LANL is needed. In response to this request, LANL Meteorology & Air Quality (MAQ) was tasked to estimate the frequency of extreme precipitation events. An earlier study by Foley (2012) which estimated extreme precipitation rates at LANL for a set of duration periods and average recurrence intervals (return periods) was evaluated. However, it was determined that this study would not be helpful since its technical support documentation was unavailable and the meteorological data and the methodology employed were not well understood. This led to the need to develop a new peer-reviewed methodology. This report documents the development of the new methodology and the application of a temporally representative precipitation dataset from February 1990 through December 2018, monitored at the official LANL meteorological station, Technical Area (TA)-6, to develop the frequency of extreme precipitation events.

54 ENVIRONMENTAL SCIENCES↗

A Review of Value of Solar Studies In Theory and In Practice

This brief summarizes a collection of state- and utility-commissioned value-of-solar (VoS) studies and related literature, with a focus on who commissioned the study, which value and cost categories were discussed and/or quantified, and the methods used. Our objective is to compile information on prior VoS studies to inform state regulators and other stakeholders that may pursue related studies or integrate findings into rate design. The brief is organized into three parts: 1) an introduction to distributed solar photovoltaic (DPV) compensation; 2) a review of theoretical research on VoS; and 3) a review of VoS studies. The vast majority of VoS studies have served an informational role of quantifying the net benefits of PV. Three studies were commissioned in states or utility service territories that subsequently implemented VoS tariffs in California, New York, and Austin, Texas. When applied as a tariff, VoS aims to compensate PV output as efficiently as possible by doing so at rates that reflect the marginal benefits and costs of PV through value and cost categories that may vary temporally and/or geographically. This could lead to higher compensation in locations and times where more PV output is more valuable and consequently drive adoption in those locations to provide more societal benefits. Value and cost factors can be broadly grouped into five categories: generation, transmission, distribution, other utility, and other social categories. Those conducting VoS studies must weigh various tradeoffs when deciding which categories to include and quantify. Tradeoffs include prioritizing values based on their magnitude of value or cost impact, as well as taking into account the feasibility of data collection and accurate quantification. Values of higher magnitude and estimation feasibility are quantified in the majority of studies, including the earliest of studies conducted in the 2000s and 2010s. Additionally, some values of higher magnitude but low feasibility in the earliest of studies have become quantifiable in recent years. There are some values with low average system-wide levels but very high magnitude in specific locations or hours. The value magnitude in some cases can be tied to DPV penetration with low value in areas with little congestion and/or low penetration and vice versa. In these cases, values that are easier to quantify are often incorporated, while those that are more difficult are often addressed via a placeholder value. The placeholder value is paired with a discussion around data needs and methods to improve future estimates, as well as a conversation about when these value categories may increase in magnitude and necessitate more rigorous quantification. This brief summarizes findings from two meta-analyses of VoS studies that took place between 2005 and 2018, as well as findings from four additional studies published from 2018 to 2023. Table ES-1 summarizes the various value and cost categories included in each respective study and whether they were quantified, discussed, or omitted. Values such as avoided energy, capacity, transmission capacity, line losses, and avoided environmental costs are quantified in every study. Some categories were deemed harder to quantify and less impactful at the time of the study, so they were discussed but not quantified (e.g., ancillary services). Other categories, including many at the distribution level, were very locationally and/or temporally specific and dependent on high DPV penetration. These were sometimes quantified and at other times discussed. Notably, when it came to utility costs, integration costs were discussed in all cases, though they were deemed to have a small impact. Other utility costs were omitted for the most part; however, the utility-commissioned study (by NorthWestern Energy in Montana) included both lost utility revenue and programmatic/administrative cost categories. While there are some similarities across studies, each had fairly unique methods that are detailed in the body of this brief.

14 SOLAR ENERGY↗

Dark energy survey year 3 results: Cosmology with peaks using an emulator approach

ABSTRACT We constrain the matter density Ωm and the amplitude of density fluctuations σ8 within the ΛCDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We take into account the uncertainty from the remaining, largely unconstrained ΛCDM parameters (Ωb, ns, and h). We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv ~\sigma _8\sqrt{\Omega _{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68 per cent confidence limits), with a precision of 1.8 per cent, an improvement of 38 per cent compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in S8 is recorded. We find a mild tension of $1.5 \, \sigma$ between our study and the results from Planck 2018, with our analysis yielding a lower S8. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment AIA and S8, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance of the CMS high-level trigger during LHC Run 2

The CERN LHC provided proton and heavy ion collisions during its Run 2 operation period from 2015 to 2018. Proton-proton collisions reached a peak instantaneous luminosity of 2.1× 10 34 cm -2 s -1 , twice the initial design value, at √(s)=13 TeV. The CMS experiment records a subset of the collisions for further processing as part of its online selection of data for physics analyses, using a two-level trigger system: the Level-1 trigger, implemented in custom-designed electronics, and the high-level trigger, a streamlined version of the offline reconstruction software running on a large computer farm. This paper presents the performance of the CMS high-level trigger system during LHC Run 2 for physics objects, such as leptons, jets, and missing transverse momentum, which meet the broad needs of the CMS physics program and the challenge of the evolving LHC and detector conditions. Sophisticated algorithms that were originally used in offline reconstruction were deployed online. Highlights include a machine-learning b tagging algorithm and a reconstruction algorithm for tau leptons that decay hadronically.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Role of the Hubble scale in the weak lensing versus CMB tension

Here, we explore a reparametrization of the lensing amplitude tension between weak lensing (WL) and cosmic microwave background (CMB) data and its implications for a joint resolution with the Hubble tension. Specifically, we focus on the lensing amplitude over a scale of 12 Mpc in absolute distance units using a derived parameter S 12 and show its constraints from recent surveys in comparison with Planck 2018. In WL alone, we find that the absolute distance convention correlates S 12 with H 0 . Accounting for this correlation in the 3D space S 12 ×w m ×ℎ reproduces the usual levels of 2~3⁢σ tension inferred from S 8 ×Ω m . Additionally, we derive scaling relations in the S 8 ×ℎ and S 12 ×ℎ planes that are allowed by Λ⁢CDM and extrapolate target scalings needed to solve the H 0 and lensing-amplitude tensions jointly in a hypothetical beyond-Λ⁢CDM model. As a test example, we quantify how the early dark energy scenario compares with these target scalings. Useful fitting formulas for S 8 and S 12 as a function of other cosmological parameters in Λ⁢CDM are provided, with 1% precision.

79 ASTRONOMY AND ASTROPHYSICS↗

Flaming Moe

Current clamp measurements collected on various small electronic devices. Details on the data set can be found in J. M. Vann, T. P. Karnowski, R. Kerekes, C. D. Cooke and A. L. Anderson, A Dimensionally Aligned Signal Projection for Classification of Unintended Radiated Emissions, in IEEE Transactions on Electromagnetic Compatibility, vol. 60, no. 1, pp. 122-131, Feb. 2018, doi: 10.1109/TEMC.2017.2692962.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Regolith, rock and fluid distributions at the Upper Colorado River Basin via a multicomponent seismic imaging approach

In October 2018, we conducted an exploratory seismic survey of the East River Watershed Function Science Focus Area, near Crested Butte, CO. Our survey goals were to map seismic bedrock, regolith and sediment properties and relate them to geomorphic and hydrologic processes. We acquired ~12 km of seismic data along roads with a 72-channel, 1.25m spaced streamer, an average shot spacing of 2.5m and an average acquisition rate of 3 km/day. We also obtained ~3km of planted geophone data with geophone and shot spacing between 2.5 and 5 meters and an acquisition rate of 0.5 km/day. Vp results show a sharp transition between regolith and bedrock along hillslopes suggesting a similarly abrupt change in porosity. Vp anisotropy in the Mancos Shale is observable in the lower reaches of the watershed, but is absent in regions near crystalline intrusions where contact metamorphism has likely modified pore structure within the shale. Alluvium thickness varies along the river axis from 0 to > 20 meters reflecting variations in depositional processes as function of stream gradient. We show that when strong discontinuities at the regolith base are present, inhomogeneous waves that travel along subsurface interfaces are generated and can be measured on vertical component geophones. We present a summary manuscript in the report appendix.

58 GEOSCIENCES↗

Data and scripts associated with a manuscript analyzing ELM-FATES parameter sensitivity under pre-fire and postfire scenarios using machine learning

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Fire Severity-Dependent Shifts in Vegetation Parameter Sensitivity: A Pre- and Post-Fire Analysis Using ELM-FATES and Explainable AI” submitted to Journal of Advances in Modeling Earth Systems (Zahura et al. 2026). The study examines vegetation physiological parameters controlling pre-fire and post-fire vegetation dynamics. To support this analysis, 73 vegetation parameters in Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (Fisher et al., 2018) , which is coupled with E3SM (Energy Exascale Earth System Model) land model (ELM, ELM-FATES), were perturbed using a Sobol sequence to generate 1,024 ensemble members for two plant functional types: needleleaf evergreen extratropical trees (NEET) and C3 grass. Simulations were conducted for the pre-fire period (2016) and post-fire period (2018–2023). Burn severity was represented by modifying the Nesterov index in FATES to 75,000, 150,000, and 300,000 for low, moderate, and high severity, respectively. A no-fire scenario was also included. Simulations were performed for 16 grid cells in the American River Watershed across different burn severities and plant functional types. XGBoost (eXtreme Gradient Boosting) models were trained using the parameter ensembles and ELM-FATES-simulated outputs, including leaf area index (LAI), gross primary productivity (GPP), aboveground biomass, vegetation evaporation, transpiration, and soil evaporation. Models were trained separately for each year and burn severity, followed by SHAP (SHapley Additive exPlanations) analysis to identify changes in dominant parameters after fire disturbance. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package contains the ELM-FATES simulation data. The scripts and data related to the analysis will be added later. The inputs and outputs from ELM-FATES are inside the “FATES” folder. “FATES_domain_surface” contains the domain and surface netcdfs that were used to run ELM-FATES in the study area. “FATES_parameters” contains the 1024 ensembles that were generated using Sobol sequence. “FATES_outputs” folder contains ELM-FATES simulated variables. All files are .csv and .nc (NetCDF).

Aboveground biomass↗

Spatial and temporal variation in forest transpiration across a forested boreal peatland complex

Transpiration is a globally important component of evapotranspiration. Careful upscaling of transpiration from point measurements is thus crucial for quantifying water and energy fluxes. In spatially heterogeneous landscapes common across the boreal biome, upscaled transpiration estimates are difficult to determine due to variation in local environmental conditions (e.g., basal area, soil moisture, permafrost). Here, we sought to determine stand-level attributes that influence transpiration scalars for a forested boreal peatland complex consisting of sparsely treed wetlands and densely treed permafrost plateaus as land cover types. The objectives were to quantify spatial and temporal variability in stand-level transpiration, and to identify sources of uncertainty when scaling point measurements to the stand-level. Using heat ratio method sap flow sensors, we determined sap velocity for black spruce and tamarack for 2-week periods during peak growing season in 2013, 2017 and 2018. We found greater basal area, drier soils, and the presence of permafrost increased daily sap velocity in individual trees, suggesting that local environmental conditions are important in dictating sap velocity. When sap velocity was scaled to stand-level transpiration using gridded 20 × 20 m resolution data across the ~10 ha Scotty Creek ForestGEO plot, we observed significant differences in daily plot transpiration among years (0.17–0.30 mm), and across land cover types. Daily transpiration was lowest in grid-cells with sparsely treed wetlands compared to grid-cells with well-drained and densely treed permafrost plateaus, where daily transpiration reached 0.80 mm, or 30% of the daily evapotranspiration. When transpiration scalars (i.e., sap velocity) were not specific to the different land cover types (i.e., permafrost plateaus and wetlands), scaled stand-level transpiration was overestimated by 42%. To quantify the relative contribution of tree transpiration to ecosystem evapotranspiration, we recommend that sampling designs stratify across local environmental conditions to accurately represent variation associated with land cover types, especially with different hydrological functioning as encountered in rapidly thawing boreal peatland complexes.

54 ENVIRONMENTAL SCIENCES↗

Hidden delays of climate mitigation benefits in the race for electric vehicle deployment

Although battery electric vehicles (BEVs) are climate-friendly alternatives to internal combustion engine vehicles (ICEVs), an important but often ignored fact is that the climate mitigation benefits of BEVs are usually delayed. The manufacture of BEVs is more carbon-intensive than that of ICEVs, leaving a greenhouse gas (GHG) debt to be paid back in the future use phase. Here we analyze millions of vehicle data from the Chinese market and show that the GHG break-even time (GBET) of China’s BEVs ranges from zero (i.e., the production year) to over 11 years, with an average of 4.5 years. 8% of China’s BEVs produced and sold between 2016 and 2018 cannot pay back their GHG debt within the eight-year battery warranty. We suggest enhancing the share of BEVs reaching the GBET by promoting the effective substitution of BEVs for ICEVs instead of the single-minded pursuit of speeding up the BEV deployment race.

33 ADVANCED PROPULSION SYSTEMS↗

Leaf mass per area and leaf water content measurements from field survey in association with NEON AOP survey, East River, CO 2018

Leaf area, wet, and dry weights for leaf samples from sites where field sampling was conducted during the 2018 National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) imaging spectroscopy and lidar surveys in Gunnison County, Colorado. The sampling sites were located across East River, Washington Gulch, Slate River, and Coal Creek watersheds and contained a mixture of meadow, shrub, and tree sampling sites. This data package contains the leaf area and wet and dry weights for leaves from tree and shrub sampling sites, as well as samples for meadow species present within a given sampling area. Associated data packages contain species composition information at the plot level (DOI: 10.15485/1618130), as well as additional biogeochemical, microbial, and geophysical data as they become available. For full documentation, please see associated reference.This package contains zipped folders of all scans that were used for area calculations for both sampling sites (tree/shrub, site_scans.zip) as well as representative samples from a daily sampling area of meadow species that were found in the 1 m plots (meadow_area_scans.zip). It also contains the macro code that was used in ImageJ for processing of the area calculations. Finally, it includes excel files with the area (cm2), wet weight (g), dry weight (g), leaf mass per area (LMA, g/m2), and leaf water content (%) for each set of scans.

54 ENVIRONMENTAL SCIENCES↗

Generator Interconnection Costs to the Transmission System in non-ISO Balancing Authorities [Slides]

Electric transmission system operators—including Independent System Operators (ISOs), Regional Transmission Organizations (RTOs), and utilities—require proposed power plants to undergo a series of interconnection studies before connecting to the grid. These studies assess what transmission upgrades or new infrastructure may be necessary and assign the associated costs to the project. Lawrence Berkeley National Laboratory has compiled, aggregated, and cleaned interconnection cost data, originally for ISOs/RTOs, and now for five non-ISO Balancing Authorities: PacifiCorp, Bonneville Power Authority, Duke Energy Progress, Duke Energy Carolinas and Duke Energy Florida. Insufficient transparency in interconnection cost data may contribute to rapidly expanding interconnection queues, with active queue capacities tripling between 2020 and 2024 in the studied BAs. Most projects withdraw after receiving high interconnection cost estimates. Interconnection costs have increased since the early 2000s, with average costs for "complete" projects reaching $194/kW between 2018 and 2024. Active queue projects and withdrawn projects incur substantially higher costs, primarily due to rising network upgrade costs. Recent interconnection costs in non-ISO balancing authorities are higher than in ISO regions, potentially due to a greater willingness to pay among developers. Utility-scale solar, wind, and storage projects have interconnection costs that exceed those for natural gas. However, when focusing on projects that do not withdraw from the queue, the interconnection costs for these technologies are more similar to natural gas projects. Other key findings include: (1) Larger generation projects benefit from lower proportional interconnection costs, (2) capacity transmission service (NRIS) often requires additional network investments, and (3) projects with high network upgrade costs are often clustered geographically. The dataset includes results from 2,104 interconnection studies conducted between 2000 and 2024, covering projects that are operational, withdrawn, or still progressing through the study process. The Excel file contains (a) the complete project-level interconnection cost dataset, and (b) seven additional tabs summarizing cost metrics across dimensions such as time, market structure, cost category (point of interconnection vs. broader network upgrades), fuel type, service type (ERIS vs. NRIS), generator size, and geography.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil" in the Journal of Geophysical Research-Biogeosciences. We simulated biogeochemical cycling in arctic soils using the PFLOTRAN geochemical model combined with measurements from previous NGEE Arctic incubations of polygonal permafrost soils in northern Alaska (Zheng et al., 2018). Simulated iron cycling, carbon dioxide production, and methane production were compared with incubation measurements and the parameterized model was then used to simulate coupled iron and carbon cycling over repeated oxic-anoxic cycles at different levels of carbon substrate availability and pH. The most recent data version (2.0) in the archive incorporates changes to the model and simulations as suggested by reviewers during the manuscript review process. These changes include an updated parameterization of the model; a new set of simulations omitting the iron cycle for direct evaluation of how iron cycle processes affect modeled outcomes; and a set of simulations testing different scenarios of carbon substrate availability in addition to scenarios of initial soil pH. This archive contains simulation code, model output, and analysis code for PFLOTRAN simulations. All scripts are python except the batch script for submitting multiprocessor jobs. Note that the model also requires compiled versions of the Alquimia interface and the NGEE Arctic fork of the PFLOTRAN geochemical simulator (see the README_INSTALL document for basic instructions). The Output directory contains eight data files in netCDF format generated by the model. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗