Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global sensitivity analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A non-intrusive optical (NIO) approach to characterize heliostats in utility-scale power tower plants: Methodology and in-situ validation

No reliable in-situ optical characterization tool is available that can efficiently and accurately measure various optical errors of individual heliostats in a utility-scale power-tower plant. Here, an innovative non-intrusive optical (NIO) approach is proposed to measure mirror surface slope error, mirror facet canting error, and heliostat tracking error based on reflection images using the natural target—the tower—in a heliostat field. The NIO approach adopts various techniques in photogrammetry, reflectometry, and geometrical optics to first determine the camera position and the actual tower position with respect to the heliostat. Then, Snell’s law is applied to find the surface normal at the reflection point or line of a reflection image and derive the surface slope error, mirror facet canting error, and heliostat tracking error accordingly. In addition, an in-situ validation is performed to demonstrate the feasibility of the NIO approach. Analysis with multiple sets of reflection images on one single heliostat of interest demonstrates that the NIO method can achieve a measurement uncertainty of 0.22 mrad for mirror facet slope error (root-mean-square) at the given test conditions. The sensitivity of measurement uncertainty to various parameters will be presented in a separate work. In the future, state-of-the-art unmanned aerial systems (UAS) are planned to perform the efficient image collection with a goal to complete the survey across a utility-scale field of over 10,000 heliostats within one or a few days depending on the number of UASs. The development of the NIO approach is expected to fill the critical gap for successful solar-field operation and maintenance, which is necessary for the global growth of the solar power-tower technology. The NIO approach can also be applied to any relevant industry requiring high-precision mirror surface shape characterization.

14 SOLAR ENERGY↗
DOE OSTI2023

Measurement of the top quark-antiquark pair charge asymmetry in events with highly boosted top quarks in proton-proton collisions at 13 TeV with the CMS detector

The top quark is the heaviest elementary particle, making it a unique tool to search for new physics. In this talk, I will present a new measurement of the top quark pair charge asymmetry for highly boosted top quarks decaying to a single lepton, missing transverse momentum and jets. The analysis is performed on 13TeV proton-proton collision data recorded by the CMS experiment during Run 2. We have defined a dedicated phase space that selects top quark-antiquark pairs with invariant mass greater than 750GeV in a semileptonic final state where the lepton is not necessarily isolated. This highly boosted sample is enhanced in valence quark production and thus expected to be more sensitive to deviations in the charge asymmetry caused by BSM processes. Dedicated tagging techniques are used to identify the decay products of hadronic top quarks and W bosons. An unfolding procedure is used to correct for detector resolution and acceptance, and inefficiencies in the event reconstruction. The result is presented in the full phase space at parton level and can be used as input to global EFT interpretations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The Cost of Floating Offshore Wind Energy in California Between 2019 and 2032

California’s energy planning is centered around meeting the emissions reduction and renewable energy requirements of Senate Bill 350 by 2030. However, state power system planning is expected to eventually address California’s requirement to achieve 100% of total retail electricity sales from renewable energy and zero-carbon resources by 2045, as mandated by Senate Bill 100. To comply with these directives, California needs to investigate the further development of energy efficiency, storage, and a diverse range of renewable energy, zero-carbon emission, and transmission resources, including offshore wind. Wind resources off the coast of California have the potential to generate a significant portion of the state’s electric energy as it moves toward a zero-carbon economy and can help diversify its energy mix. Floating offshore wind technology, which is suitable for the deep waters along the California coast, is currently in a precommercial phase, with approximately 84 megawatts (MW) installed worldwide at the end of 2019. Globally there are over 7,000 MW in planning and permitting phases of development, with the first commercial-scale projects expected to be operational in 2024. This study provides site-specific cost and performance data for floating offshore wind to inform California’s long-term energy planning. The identification of new resources to meet California’s policy goals at least cost is part of the Integrated Resource Planning (IRP) process, which is coordinated by the California Public Utilities Commission (CPUC). In 2019–2020 IRP modeling, offshore wind was included for the first time as a candidate resource in some sensitivity cases (CPUC 2019). The data and information presented in this report can be used to update offshore wind inputs in future IRP cycles. The authors conducted a geospatial cost analysis over portions of the offshore wind resource area of California. The analyzed spatial domain includes sites with a mean wind speed of at least 7 meters per second and water depths between 40 meters (m) and 1,300 m. Costs and energy production vary across this analysis domain. We calculated these parameters on a grid layout with over 750 sites, with each site representing a 1,000-MW commercial offshore wind power plant. Levelized cost of energy (LCOE) was calculated at each site over the analysis domain. The resulting variation in LCOE across the analysis domain is illustrated through heat maps in this report. Five study areas were selected within the analysis domain where more detailed cost analysis was conducted and cost parameters, such as annual energy production, capital cost expenditures (CapEx), operational cost expenditures (OpEx), and net capacity factors are reported. These five study areas include Morro Bay, Diablo Canyon, Humboldt, Cape Mendocino and Del Norte.

17 WIND ENERGY↗
DOE OSTI2022

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Using onset times from frequent seismic surveys to understand fluid flow at the Peace River Field, Canada

With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.

58 GEOSCIENCES↗
DOE OSTI2026

IAM-FIRE: a Climate Emulator–Based Framework to Project Wildfire Impacts and Risks for Integrated Assessment Models

Most Integrated Assessment Models (IAMs) underrepresent dynamic feedbacks from climate-driven disturbances such as wildfires, potentially overestimating the permanence of land-based carbon sinks. In particular, representing the impacts of forest fires is becoming increasingly important, as these are expected to intensify in the coming years. We introduce IAM-FIRE (Integrated Assessment Model – Fire Impacts & Risks Emulator), a novel framework that enables the projection of wildfire burned area (BA) and carbon emissions (CE) directly from IAM outputs. IAM-FIRE combines a spatial climate emulator, land-use downscaling, vegetation productivity modelling, and an empirical fire model to generate global annual wildfire impacts for arbitrary socioeconomic and emissions scenarios at 0.5° resolution for the period 2020–2100. Calibrated against GFEDv5 observations and using inputs from the Global Change Analysis Model (GCAM), we report projections BA and CE derived from IAM-FIRE for four scenarios: SSP1-2.6, SSP2-4.5, SSP3-6.6 and SSP5-7.6. The model reproduces historical global trends for total BA, including the observed global decline since the early 2000s, and for forest BA. Projected fire trajectories differ strongly among scenarios: total BA range from declines under SSP1-2.6 (-3.36 Mha yr-1) to increases under SSP3-6.6 (+1.6 Mha yr-1). Corresponding total CE show a similar divergence ranging from -15 to +10.6 TgC yr-1. Socioeconomic development exerts a dominant suppressing effect on wildfire impacts while climate change and CO2-driven increases in vegetation productivity amplify fire risk, particularly under high-emissions pathways. Compared with CMIP6 fire-enabled Earth System Models, IAM-FIRE exhibits greater sensitivity to radiative forcing and a stronger role for human-driven fire suppression, highlighting substantial structural uncertainties in future fire projections. By providing a computationally efficient and internally consistent approach to represent wildfire impacts within IAMs, IAM-FIRE enables systematic exploration of fire–climate–land feedbacks and supports improved assessments of mitigation permanence and climate risks in future integrated scenarios.

Rouhette, Theo↗
DOE OSTI2023

ETG turbulent transport in the Mega Ampere Spherical Tokamak (MAST) pedestal

The results of a gyrokinetic analysis of turbulent transport driven by the electron temperature gradient (ETG) in the MAST pedestal are presented. Local nonlinear simulations from the gyrokinetic code, GENE, show that heat flux produced by ETG turbulence is 10%–30% of the total applied heating power in the upper pedestal and pedestal top during both the pre-edge-localized mode (ELM) (80%–99% inter-ELM period) and post-ELM (0%–20%) periods. Increasing strongly with the ETG, the ETG transport appears to be stiff. Considering radiation losses, ion thermal transport, and the strong sensitivity of the transport to the ETG, we propose that ETG transport is a plausible mechanism mediating the inter-ELM temperature profile on MAST. Cognizant of the possibility that sharp variations in the pedestal profiles may violate the local approximation, we conducted global nonlinear simulations; the results are in good agreement with local simulations except near the pedestal top, where extended radial structures and high transport levels (far beyond experimental) develop in the local simulations. We quantify and parameterize the discrepancy between local and global simulations by calculating the ratio of the radial correlation length to a length scale representative of the profile curvature. When this ratio is sufficiently small, local and global simulations agree as expected.

ETG↗
DOE OSTI2020

Argonne’s global critical materials agent-based model (GCMat)

Several studies have identified rare earths as critical materials. Although reasonably abundant in the Earth’s crust, rare earths typically occur in low concentrations of mined ores. Processes for recovering rare earth concentrates from these ores are complex and capital intensive. Further, lead times for deposit development, licensing, and construction are long, with reports of 10- 15 years. China is a major player in the rare earths supply chain, both in production capacity and technology innovation. In 2017, China supplied more than 80% of global rare earth oxide demand. Rare earth elements (REEs) have unique magnetic, catalytic, and phosphorescent properties that significantly improve performance of a wide range of technologies. These technologies span aerospace, energy, telecommunications, electronics, transportation, defense, and other diverse applications. Consequently, disruptions in rare earth supply can have a significant societal impact. Estimating that impact requires understanding of the dynamics across the supply chain, from rare earth oxide extraction to end use application. This report documents Argonne’s Global Critical Materials model (GCMat), first described by Riddle et al. 2015. GCMat provides capabilities to explore supply chain dynamics and uncertainty under scenarios of demand growth or shrinkage, technology adoption, supply disruptions, and trade policies and mitigation strategies of new supply sources, product substitution, consumer thrifting, and stockpiling. Supply chain participants from rare earth mining through final demand are modeled as interacting agents who make market decisions independently as time progresses. Since the version documented in Riddle et al. 2015, GCMat has been expanded to cover additional REEs, derived products and supply chains that use these REEs, and includes new agent behaviors and modeling capabilities. Section 2 provides a summary of the GCMat model design, including the structure of the model and key assumptions, section 3 summarizes methods used for model calibration and sensitivity analysis, and section 4 provides examples of model results.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Understanding neutrinos with accelerator beams and liquid argon time-projection chambers: ICARUS and DUNE

A global program of experiments has worked towards characterizing neutrino oscillation over the past few decades. However, important parameters remain to be measured, and mysteries remain to be elucidated. Current and upcoming experiments are targeting the open questions and probing the consistency of the neutrino oscillation paradigm. Likewise, the liquid argon (LAr) time-projection chamber (TPC) has emerged as a sensitive particle detection technology for neutrino experiments. A current generation of LAr TPC detectors are being used to study neutrinos while also gaining important experience in operating and analyzing with this technology. SBND, MicroBooNE, and ICARUS have collected or are collecting data from beams at Fermilab (near Chicago) to explore the possibility of a sterile neutrino and/or other beyond Standard Model (BSM) physics. SBND and ICARUS will be used to conduct a two-detector analysis as part of the Short Baseline Neutrino (SBN) Program. Additionally, these detectors are enabling important neutrino interaction studies necessary to prepare for the next generation of oscillation experiments. One such oscillation experiment that will come online over the next years is the Deep Underground Neutrino Experiment (DUNE), which will install multiple 10 kiloton LAr TPCs underground in South Dakota (south of Saskatchewan) to conduct oscillation measurements with neutrinos originating in a beamline at Fermilab. A detector complex will be installed at Fermilab as well, to study the beam before the expected flavour oscillations. This “near detector” will also employ a LAr TPC, with a segmented and pixel-based design, as well as other technologies to constrain uncertainties in the oscillation measurement by characterizing the beam and neutrino interactions. This talk will discuss the ICARUS and DUNE experiments, the LAr TPC detector technology, and the efforts to realize and leverage these experiments to better understand the properties of neutrinos.

Howard, Bruce L. [Fermilab]↗
DOE OSTI2026

Life cycle analysis of dedicated energy crops for fuel production in the United States

Dedicated energy crops are promising feedstocks to make biofuels including jet fuels. This study applies life cycle analysis (LCA) to estimate direct well-to-wake (WTW) greenhouse gas (GHG) emissions (g CO 2 e/MJ) for jet fuel derived from five energy crops-biomass sorghum, miscanthus, switchgrass, poplar, and willow-via Fischer-Tropsch-to-Jet (FTJ) and Ethanol-to-Jet (ETJ) pathways. The WTW boundary includes direct emissions from biomass production, fuel production, and fuel combustion. The R&D GREET model is expanded to conduct the LCA, using national average biomass yields and farming inputs from the 2023 Billion-Ton Study. In addition, this study estimates emissions from market-mediated effects, including induced land-use change, induced other crop (non-feedstock) production changes, and induced livestock production changes using global economic and emissions factor models. On a per-dry U.S. ton basis, cultivation and harvest emissions are lowest for willow (51,565 g CO 2 e) and highest for biomass sorghum (104,488 g CO 2 e). Per-acre results show similarly high emissions for sorghum and lowest values for poplar and willow. Direct WTW emissions are substantially lower for FTJ (biomass sorghum: 5.5; miscanthus: 10.3; switchgrass: 11.7; poplar: 11.9; and willow: 8.7 g CO 2 e/MJ) than ETJ (33.2; 33.8; 34.8; 36.2; and 31.7 g CO 2 e/MJ, respectively). When market-mediated emissions are included, miscanthus exhibits the lowest total emissions across energy crop pathways. Although results are sensitive to modeling assumptions, they indicate that high-yielding perennial and woody crops, particularly when planted on marginal land, could significantly reduce WTW emissions for bio-jet fuels by combining low direct emissions with soil carbon gains and favorable market-mediated effects.

Billion-Ton Study↗
DOE OSTI2024

Quantification of locked mode instability triggered by a change in confinement

Abstract This work presents the first analysis of the disruptive locked mode (LM) triggered by the dynamics of a confinement change. It shows that, under certain conditions, the LM threshold during the transient is significantly lower than expected from steady states. We investigate the sensitivity to a controlled n = 1 error field (EF) activated prior to the L-H transition in the COMPASS tokamak, at q 95 ∼ 3, β N ∼ 1, and using EF coils on the high-field side of the vessel. A threshold for EF penetration subsequent to the L-H transition is identified, which shows no significant trend with density or applied torque, and is an apparent consequence of the reduced intrinsic rotation of the 2/1 mode during this transient phase. This finding challenges the assumption made in theoretical and empirical works that natural mode rotation can be predicted by global plasma parameters and urges against using any parametric EF penetration scaling derived from steady-state experiments to define the EF correction strategy in the entire discharge. Furthermore, even at EFs below the identified penetration threshold, disruptive locking of sawtooth-seeded 2/1 tearing modes is observed after about 30% of L-H transitions without external torque.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Molten Salt Aerosol Spectroscopy Approach to Online Process Monitoring in Molten Salt Reactors

Recent research and development in molten salt reactors (MSRs) has led to many advanced reactor designs, many of which have demonstrations planned in the next decade. However, a significant challenge that must be overcome before liquid fueled MSRs become viable globally is nuclear material accountancy (NMA) and safeguards. Liquid-fueled MSRs are problematic to safeguards because the nuclear material is in the bulk phase, inventories are constantly changing, and there are potentially many outlet pathways that must be monitored, especially in designs with online chemical processing. At the Idaho National Laboratory, we have developed a process monitoring tool to circulate molten salt through a nebulizer and optical cell, where laser-induced breakdown spectroscopy (LIBS) and ultraviolet-visible (UV-VIS) are applied to track the elemental and oxidation states of the constituents in the salt aerosol. This approach has the potential to provide information on the actinides, fission products, and corrosion in the reactor or process. Preliminary testing in aqueous praseodymium chloride (PrCl3) and neodymium chloride (NdCl3) showed good sensitivity detection limits n developed calibration curves. Additional tests and analysis or ongoing in aqueous solutions while equipment is being scaled up for high temperature molten salt testing. Results for continued testing as well as a status of the molten salt design and anticipated commissioning will be presented.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

An Adaptive Tracking Algorithm for Convection in Simulated and Remote Sensing Data

A robust and computationally efficient object tracking algorithm is developed by incorporating various tracking techniques. Physical properties of the objects, such as brightness temperature or reflectivity, are not considered. Therefore, the algorithm is adaptable for tracking convection-like features in simulated data and remotely sensed two-dimensional images. In this algorithm, a first guess of the motion, estimated using the Fourier phase shift, is used to predict the candidates for matching. A disparity score is computed for each target–candidate pair. The disparity also incorporates overlapping criteria in the case of large objects. Then the Hungarian method is applied to identify the best pairs by minimizing the global disparity. The high-disparity pairs are unmatched, and their target and candidate are declared expired and newly initiated objects, respectively. They are tested for merger and split on the basis of their size and overlap with the other objects. The sensitivity of track duration is shown for different disparity and size thresholds. The report highlights the algorithm’s ability to study convective life cycles using radar and simulated data over Darwin, Australia. The algorithm skillfully tracks individual convective cells (a few pixels in size) and large convective systems. The duration of tracks and cell size are found to be lognormally distributed over Darwin. The evolution of size and precipitation types of isolated convective cells is presented in the Lagrangian perspective. This algorithm is part of a vision for a modular platform [viz., TINT is not TITAN (TINT) and Tracking and Object-Based Analysis of Clouds (tobac)] that will evolve into a sustainable choice to analyze atmospheric features.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ~$\frac{1}{3}$ of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2% for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluating the impact of peat soils and snow schemes on simulated active layer thickness at pan-Arctic permafrost sites

Abstract Permafrost stability is significantly influenced by the thermal buffering effects of snow and active-layer peat soils. In the warm season, peat soils act as a barrier to downward heat transfer mainly due to their low thermal conductivity. In the cold season, the snowpack serves as a thermal insulator, retarding the release of heat from the soil to the atmosphere. Currently, many global land models overestimate permafrost soil temperature and active layer thickness (ALT), partially due to inaccurate representations of soil organic matter (SOM) density profiles and snow thermal insulation. In this study, we evaluated the impacts of SOM and snow schemes on ALT simulations at pan-Arctic permafrost sites using the Energy Exascale Earth System Model (E3SM) land model (ELM). We conducted simulations at the Circumpolar Active Layer Monitoring (CALM) sites across the pan-Arctic domain. We improved ELM-simulated site-level ALT using a knowledge-based hierarchical optimization procedure and examined the effects of precipitation-phase partitioning methods (PPMs), snow compaction schemes, and snow thermal conductivity schemes on simulated snow depth, soil temperature, ALT, and CO 2 fluxes. Results showed that the optimized ELM significantly improved agreement with observed ALT (e.g. RMSE decreased from 0.83 m to 0.15 m). Our sensitivity analysis revealed that snow-related schemes significantly impact simulated snow thermal insulation levels, soil temperature, and ALT. For example, one of the commonly used snow thermal conductivity schemes (quadratic Sturm or SturmQua) generally produced warmer soil temperatures and larger ALT compared to the other two tested schemes. The SturmQua scheme also amplified the model’s sensitivity to PPMs and predicted deeper ALTs than the other two snow schemes under both current and future climates. The study highlights the importance of accurately representing snow-related processes and peat soils in land models to enhance permafrost dynamics simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Changes in Annual Extremes of Daily Temperature and Precipitation in CMIP6 Models

This study presents an analysis of daily temperature and precipitation extremes with return periods ranging from 2 to 50 years in phase 6 of the Coupled Model Intercomparison Project (CMIP6) multimodel ensemble of simulations. Judged by similarity with reanalyses, the new-generation models simulate the present-day temperature and precipitation extremes reasonably well. In line with previous CMIP simulations, the new simulations continue to project a large-scale picture of more frequent and more intense hot temperature extremes and precipitation extremes and vanishing cold extremes under continued global warming. Changes in temperature extremes outpace changes in global annual mean surface air temperature (GSAT) over most landmasses, while changes in precipitation extremes follow changes in GSAT globally at roughly the Clausius–Clapeyron rate of ~7% °C –1 . Changes in temperature and precipitation extremes normalized with respect to GSAT do not depend strongly on the choice of forcing scenario or model climate sensitivity, and do not vary strongly over time, but with notable regional variations. Over the majority of land regions, the projected intensity increases and relative frequency increases tend to be larger for more extreme hot temperature and precipitation events than for weaker events. To obtain robust estimates of these changes at local scales, large initial-condition ensemble simulations are needed. Appropriate spatial pooling of data from neighboring grid cells within individual simulations can, to some extent, reduce the needed ensemble size.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Discovering the 3D Structure and Dynamics of the Sun-Interstellar Medium System on a Global Scale Los Alamos LDRD Report

We discover the 3D plasma structure of the Sun-interstellar medium interaction region, called the heliosheath. This is enabled by (1) the completion of a full solar cycle of observations of hydrogen energetic neutral atoms (ENAs) emanating from this region and detected by the Los Alamos-led IBEX-Hi ENA imager on NASA’s Interstellar Boundary Explorer (IBEX) mission, and (2) the development of tools to fully exploit all data from the IBEX-Hi imager, leading to a significant increase in the effective instrument sensitivity. This effort prepares us for optimizing the design and data exploitation of the IMAP-Hi imager, the Los Alamos-led next generation imager that will fly on the Interstellar Mapping and Acceleration Probe (IMAP) mission, currently under production. This project realized its goals as evidenced by two main outcomes: (1) The publication of the first-ever empirical 3D map of the interstellar boundary, and (2) the development and release of a new validated ENA data set that has tripled the amount of data available for scientific analysis of the outer heliosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗