Search NASA⌕ Search

SEARCH · Search NASA

Results for “ENERGY BALANCE”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗
DOE OSTI2022

Proteomic Profiling of Intra-Islet Features Reveals Substructure-Specific Protein Signatures

Despite their diminutive size, islets of Langerhans play a large role in maintaining systemic energy balance in the body. New technologies have enabled us to go from studying the whole pancreas, to isolated whole islets, to partial islet sections, and now to islet substructures isolated from within the islet. Using a microfluidic nanodroplet-based proteomics platform coupled with laser capture microdissection (LCM), we present an in-depth investigation of protein profiles specific to regions within the islet. These regions studied include the vascular tissue containing interface boundary, micro-vascular tissue internal to the islet, isolated endocrine cells, islet sections with vasculature intact, and finally acinar tissue from around the islet. Unique protein signatures observed in the inner vasculature potentially indicate increased innervation and intra-islet neuron-like crosstalk compared to external vasculature. We also demonstrate the utility of these data for isolating localized structure-specific drug-target interactions using existing protein/drug binding databases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Metabolic and functional connections between cytoplasmic and chloroplast triacylglycerol storage

Neutral lipids in the form of triacylglycerol (TAG) have emerged as critical regulators of cellular energy balance, lipid homeostasis, growth, development and stress response in organisms ranging from plants to yeast. Although TAGs are mostly recognized as the main storage component in cytoplasmic lipid droplets (LDs), TAG-rich LDs with similar structural and functional characteristics to those found in the cytoplasm also exist in chloroplasts of microalgae and higher plants. Chloroplasts contain up to 70% of total lipids in photosynthetic cells, yet how organisms maintain chloroplast lipid homeostasis remains an under-investigated area of research. In this work, we summarize the current state of knowledge about the metabolism of TAG and its function in chloroplasts, with a focus on the enzymes catalyzing the final steps of TAG assembly and the role of TAG synthesis in protection against lipotoxicity. We also discuss emerging data regarding connections between cytoplasmic and chloroplast TAG metabolism and the role of autophagy in the degradation of chloroplast storage lipids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Methodology for inferring reactor core power distribution from an optical fiber based gamma thermometer array

Here, a 4-step data analytic methodology has been devised for the purpose of inferring the distribution of power in a reactor core, based on the response of an array of optical fiber based gamma thermometers (OFBGTs). This data analytic methodology is crucial for the development of a system of OFBGTs for the purpose of calibrating local power range monitors in boiling water reactors. Such a system would be an improvement to the present calibration system in boiling water reactors, in terms of safety, efficiency, and permanence. The first step of this methodology is to establish an energy balance method. In this method, one uses MCNP to determine response functions, which allow one to convert from gamma thermometer response to power. The gamma thermometers and the reactor core are segmented, such that each gamma thermometer segment provides an estimate of power for each reactor core segment. The estimates of power for the reactor core segments (hereafter referred to as fuel assembly segments) are calculated based on the measured response of the gamma thermometers. The second step of the methodology is to employ a weighting scheme to combine the estimates of the power of the various fuel assembly segments, based on the response of the various OFGBT segments and the incremental dose rates of the various fuel assembly segments to the OFBGT segments. The third step is to iteratively calculate the estimates of the power of the fuel assembly segments, until a convergence criterion is met, which indicates that the calculation has converged. The fourth and final step of the methodology is to estimate the uncertainty of the power for each fuel assembly segment. The mathematical basis for this step is not the focus of this paper. We have used a 3-D homogeneous reactor model to demonstrate the data analytic methodology; and have found that the data analytic methodology operates as intended.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A radiative transfer model for solar induced fluorescence using spectral invariants theory

Solar Induced chlorophyll Fluorescence (SIF) shows promise as an approach for estimating gross primary production (GPP) remotely. However, sun-target-sensor geometry and within-canopy absorption of SIF can alter the relationship between measured SIF and GPP, because sensors can only retrieve some unknown fraction of the total emitted SIF. Radiative transfer models that allow for variation in canopy structure and sensor angles are therefore needed to properly interpret SIF measurements. Spectral invariants allow decoupling of the wavelength-independent canopy structure and the wavelength-dependent leaf and soil spectrum in the radiative transfer process. Here we develop a simple analytical Fluorescence Radiative Transfer model based on Escape and Recollision probability (FluorRTER) to investigate the impact of canopy structure and sun-target-sensor geometry on SIF emissions. SIF simulations using the FluorRTER model agreed well the one-dimensional Soil-Canopy Observation of Photochemistry and Energy balance (SCOPE) model and the three-dimensional Fluorescence model with Weighted Photon Spread (FluorWPS) model. The fractional vegetation cover (FVC) and clumping effect have a large influence the SIF emission of 3D discontinuous canopies. For a moderate solar zenith angle (30°) and a clumped canopy (FVC?=?0.6), the difference between the directional observed SIF of a 3D discontinuous canopy and a 1D homogeneous canopy was as large as 43.2% and 38.4% for Photosystem I?+?II fluorescence at 685?nm and at 740?nm, respectively. By bridging the gap between observed SIF and total emitted SIF over 3D heterogeneous vegetation canopies, the FluorRTER model can assist with the angular normalization of SIF measurements and enable the more robust interpretation of how variations in SIF from directional and hemispherical in-situ, airborne and satellite observations relate to leaf and whole-canopy physiological processes.

Zeng, Yelu↗
DOE OSTI2022

What lies beneath: Vertical temperature heterogeneity in a Mediterranean woodland savanna

As the availability of satellite and airborne thermal infrared remote sensing (TIR-RS) data increases and their spatial, temporal, and spectral resolutions improve, researchers are finding diverse applications for TIR-RS measurements. TIR-RS is now commonly applied in regional- and continental-scale analyses, such as those focused on fire and surface energy balance. However, its application lags in plant physiology and ecology, for which a finer-scale understanding of plant canopy temperatures would be useful to elucidate plant water dynamics, for example. In particular, while methods to disaggregate TIR-RS pixels in horizontal space have advanced, possible vertical stratification of plant canopy temperature and its implications for understanding the correspondence between TIR-RS and finer-scale, field-based thermal measurements (e.g. made with a thermal camera) remain unexplored. Here, we use data from a thermal camera deployed concurrently with the recent ECOSTRESS mission to quantify vertical temperature gradients within tree canopies and temperatures of over- vs. under-story plants in a Mediterranean woodland savanna. We then leverage diverse ancillary data to maximize the geometric comparability of ECOSTRESS and thermal camera measurements, in order to assess the extent to which the two forms of thermal measurements correspond. Specifically, we ask: (1) What are the patterns of intra-canopy and over- vs. under-story vertical temperature in a Mediterranean woodland savanna?, and (2) How can vertically-resolved, but spatially-limited field-based temperature measurements be reconciled with spatially-extensive, but surface-only, temperature measurements of a space-borne remote sensor? Further, we found consistent patterns of vertical thermal heterogeneity both within tree canopies and between ecosystem over- and under-stories. The daytime difference between the top and bottom thirds of blue oak canopies was, on average, 0.48 ° C – and sometimes several times larger. Notably, canopy tops are cooler, likely associated with the under-story grass reaching daytime temperatures often exceeding over-story temperatures by 10° C. Given the consistency of the intra-canopy temperature gradients, we expected the ECOSTRESS sensor would be in better agreement with camera measurements of canopy tops than bulk canopies or canopy bottoms. However, within-canopy gradients were overwhelmed by other sources of disagreement between the measurements, in part associated with upscaling camera measurements across space. Overall, thermal camera and ECOSTRESS measurements were largely in agreement at night (pixel RMSE = 1.1°C), but they were more divergent during times of low (but >0 W/m2) and high incoming solar radiation (daytime pixel RMSE = 3.5°C).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Enhanced Collisional Losses from a Magnetic Mirror Using the Lenard-Bernstein Collision Operator

Collisions are crucial in governing particle and energy transport in plasmas confined in a magnetic mirror trap. Modern gyrokinetic codes model transport in magnetic mirrors, but some utilize approximate model collision operators. This study focuses on a Pastukhov-style method of images calculation of particle and energy confinement times using a Lenard-Bernstein model collision operator. Prior work on parallel particle and energy balances used a different Fokker-Planck plasma collision operator. The method must be extended in non-trivial ways to study the Lenard-Bernstein operator. To assess the effectiveness of our approach, we compare our results with a modern finite element solver. Our findings reveal that the particle confinement time scales like a exp( a 2 ) using the Lenard-Bernstein operator, in contrast to the more accurate scaling that the Coulomb collision operator would yield a 2 exp( a 2 ), where a 2 is approximately proportional to the ambipolar potential. We propose that codes solving for collisional losses in magnetic mirrors utilizing the Lenard-Bernstein or Dougherty collision operator scale their collision frequency of any electrostatically confined species. This study illuminates the collision operator’s intricate role in the Pastukhov-style method of images calculation of collisional confinement.

fusion plasma↗
DOE OSTI2020

Molecular Specificity and Proton Transfer Mechanisms in Aerosol Prenucleation Clusters Relevant to New Particle Formation

Atmospheric aerosol particles influence Earth’s radiative energy balance and cloud properties either directly through interacting with the solar radiation or indirectly through acting as cloud condensation nuclei (CCN), thus impacting Earth’s climate change, air quality, and human health. Regarding both the important scientific and overarching practical implications of aerosols, the past two decades have seen an extensive research effort, with emphasis on the chemical compositions and underlying mechanisms of aerosol formation. New particle formation (NPF) has been recognized as an important pathway, contributing up to 50% of atmospheric aerosols. Nowadays, the general consensus in NPF is that it proceeds via two distinct stages—the nucleation from gaseous precursors to form critical nuclei of sub-1–2 nanometer size, and the subsequent growth into large particles. However, a fundamental understanding of both the NPF process and the molecular-level characterization of the critical size aerosol clusters is still largely missing, hampering the development of reliable and predictive aerosol nucleation and climate models.

Hou, Gao-Lei↗
DOE OSTI2023

Toggling Stereochemical Activity through Interstitial Positioning of Cations between 2D V 2 O 5 Double Layers

The 5/6s 2 lone-pair electrons of p-block cations in their lower oxidation states are a versatile electronic and geometric structure motif that can underpin lattice anharmonicity and often engender electronic and structural instabilities that underpin the function of active elements in nonlinear optics, thermochromics, thermoelectrics, neuromorphic computing, and photocatalysis. In contrast to periodic solids where lone-pair-bearing cations are part of the structural framework, installing lone-pair-bearing cations in the interstitial sites of intercalation hosts provides a means of a systematically modulating electronic structure through the choice of the group and the period of the inserted cation while preserving the overall framework connectivity. The extent of stereochemical activity and the energy positioning of lone-pair-derived mid-gap states depend on the cation identity, stoichiometry, and strength of anion hybridization. V 2 O 5 polymorphs are versatile insertion hosts that can accommodate a broad range of s-, p-, and d-block cations. However, the insertion of lone-pair-bearing cations remains largely underexplored. In this article, we examine the implications of varying the 6s 2 cations situated in interlayer sites between condensed [V 4 O 10 ] n double layers. Systematic modulations of lattice distortions, electronic structure, and magnetic ordering are observed with increasing strength of stereochemical activity from group 12 to group 14 cations. We compare and contrast p-block-layered M x V 2 O 5 (M = Hg, Tl, and Pb) compounds and map the significance of local off-centering arising from the stereochemical activity of lone-pair cations to the emergence of filled antibonding lone-pair 6s 2 –O 2p-hybridized mid-gap states mediated by second-order Jahn–Teller distortions. Crystallographic studies of cation coordination environments and the resulting modulation of V–V interactions have been used in conjunction with variable-energy hard X-ray photoelectron spectroscopy measurements, first-principles electronic structure calculations, and crystal orbital Hamilton population analyses to decipher the origins of stereochemical activity. Magnetic susceptibility measurements reveal antiferromagnetic signatures for all the three compounds. However, the differences in V–V interactions significantly affect the energy balance of the superexchange interactions, resulting in an ordering temperature of 160 and 260 K for Hg 0.5 V 2 O 5 and δ-Tl 0.5 V 2 O 5 , respectively, as compared to 7 K for δ-Pb 0.5 V 2 O 5 . In δ-Pb 0.5 V 2 O 5 , the strong stereochemical activity of electron lone pairs and the resulting electrostatic repulsions enforce superlattice ordering, which strongly modifies the electronic localization patterns along the [V 4 O 10 ] slabs, resulting in disrupted magnetic ordering and an anomalously low ordering temperature. The results demonstrate a versatile strategy for toggling the stereochemical activity of electron lone pairs to modify the electronic structure near the Fermi level and to mediate superexchange interactions.

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

New Particle Formation and Growth Dynamics for α-Pinene Ozonolysis in a Smog Chamber and Implications for Ambient Environments

Secondary organic aerosol (SOA) plays an important role in the Earth’s energy balance and air quality. Understanding the formation of SOA accompanying the oxidation of volatile organic compounds (VOCs) remains incomplete because of the complex reactions and partitioning mechanisms. In this study, we investigated new particle formation (NPF) of SOA derived from α-pinene ozonolysis in the Harvard Environmental Chamber using the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the condensing species as two products, one of lower volatility (LVOC, 0.01 μg m –3 ) and one of higher volatility (SVOC, 1.0 μg m –3 ). A module for the classical nucleation theory was included. The surface tension (σ), bulk diffusivity (D b ), and mass accommodation coefficient (α) were varied until convergence between the simulated and observed particle number size distributions. α of LVOC was constrained as 0.03–0.1 with the reported σ of α-pinene-ozonolysis-derived SOA and Tolman surface correction (~23.0–27.5 mN m –1 ) and Db from 10 –15.5 to >10 –11 cm 2 s –1 . The observed particle number size distribution of α-pinene ozonolysis was well-simulated to derive the physical parameters, including σ, Db, and α, suggesting the possible application to group the reported complex oxidation products and further application for other atmospheric-related oxidations for quantifying NPF.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mesoscale Frank–Kasper Crystal Structures from Dendron Assembly by Controlling Core Apex Interactions

Single-component polymeric materials open up a great potential for self-assembly into mesoscale complex crystal structures that are known as Frank-Kasper (FK) phases. Predicting the packing structures of the soft-matter spheres, however, remains a challenge even when the molecular design is precisely known. Here, we investigate the role of the molecules' enthalpic interaction in determining the low-symmetry crystal structures. To this end, we synthesize architecturally asymmetric dendrons by varying their apex functionalities and examine the packing structures of the second-generation (G2) dendritic wedges. Our work shows that weakening the hydrogen bonding of the dendron apex makes the particles softer and smaller, and leads to the formation of various FK structures at lower temperatures, including the new observation of a FK C14 phase in the cone-shaped dendron systems. As a consequence of the free energy balance between the particle's interfacial tension and the chain's stretching, various packing structures are mainly tuned by designing the hydrogen bonding interaction.

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

Energy Budget Constraints on the Time History of Aerosol Forcing and Climate Sensitivity

Abstract An observationally constrained time series of historical aerosol effective radiative forcing (ERF) from 1750 to 2019 is developed in this study. We find that the time history of aerosol ERFs diagnosed in CMIP6 models exhibits considerable variation and explore how the time history of aerosol forcing influences the probability distributions of present‐day aerosol forcing and emergent metrics such as climate sensitivity. Using a simple energy balance model, trained on CMIP6 climate models and constrained by observed near‐surface warming and ocean heat uptake, we derive estimates for the historical aerosol forcing. We find 2005–2014 mean aerosol ERF to be −1.1 (−1.8 to −0.5) W m −2 relative to 1750. Assuming recently published historical emissions from fossil fuel and industrial sectors and biomass burning emissions from SSP2‐4.5, aerosol ERF in 2019 is −0.9 (−1.5 to −0.4) W m −2 . There is a modest recovery in aerosol forcing (+0.025 W m −2 decade −1 ) between 1980 and 2014. This analysis also gives a 5%–95% range of equilibrium climate sensitivity of 1.8°C –5.1°C (best estimate 3.1°C) with a transient climate response of 1.2°C –2.6°C (best estimate 1.8°C).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Description and Demonstration of the Coupled Community Earth System Model v2 – Community Ice Sheet Model v2 (CESM2‐CISM2)

Abstract Earth system/ice‐sheet coupling is an area of recent, major Earth System Model (ESM) development. This work occurs at the intersection of glaciology and climate science and is motivated by a need for robust projections of sea‐level rise. The Community Ice Sheet Model version 2 (CISM2) is the newest component model of the Community Earth System Model version 2 (CESM2). This study describes the coupling and novel capabilities of the model, including: (1) an advanced energy‐balance‐based surface mass balance calculation in the land component with downscaling via elevation classes; (2) a closed freshwater budget from ice sheet to the ocean from surface runoff, basal melting, and ice discharge; (3) dynamic land surface types; and (4) dynamic atmospheric topography. The Earth system/ice‐sheet coupling is demonstrated in a simulation with an evolving Greenland Ice Sheet (GrIS) under an idealized high CO 2 scenario. The model simulates a large expansion of ablation areas (where surface ablation exceeds snow accumulation) and a large increase in surface runoff. This results in an elevated freshwater flux to the ocean, as well as thinning of the ice sheet and area retreat. These GrIS changes result in reduced Greenland surface albedo, changes in the sign and magnitude of sensible and latent heat fluxes, and modified surface roughness and overall ice sheet topography. Representation of these couplings between climate and ice sheets is key for the simulation of ice and climate interactions.

Muntjewerf, Laura↗
DOE OSTI2022

Multidisciplinary Constraints on the Thermal‐Chemical Boundary Between Earth's Core and Mantle

Abstract Heat flux from the core to the mantle provides driving energy for mantle convection thus powering plate tectonics, and contributes a significant fraction of the geothermal heat budget. Indirect estimates of core‐mantle boundary heat flow are typically based on petrological evidence of mantle temperature, interpretations of temperatures indicated by seismic travel times, experimental measurements of mineral melting points, physical mantle convection models, or physical core convection models. However, previous estimates have not consistently integrated these lines of evidence. In this work, an interdisciplinary analysis is applied to co‐constrain core‐mantle boundary heat flow and test the thermal boundary layer (TBL) theory. The concurrence of TBL models, energy balance to support geomagnetism, seismology, and review of petrologic evidence for historic mantle temperatures supports Q CMB ∼15 TW, with all except geomagnetism supporting as high as ∼20 TW. These values provide a tighter constraint on core heat flux relative to previous work. Our work describes the seismic properties consistent with a TBL, and supports a long‐lived basal mantle molten layer through much of Earth's history.

58 GEOSCIENCES↗
DOE OSTI2022

Physical Controls on the Hydrology of Perennially Ice-Covered Lakes, Taylor Valley, Antarctica (1996–2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. Here, we apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.

58 GEOSCIENCES↗
DOE OSTI2022

Disentangling the Relative Drivers of Seasonal Evapotranspiration Across a Continental-Scale Aridity Gradient

Evapotranspiration (ET) is a significant ecosystem flux, governing the partitioning of energy at the land surface. Understanding the seasonal pattern and magnitude of ET is critical for anticipating a range of ecosystem impacts, including drought, heat-wave events, and plant mortality. Here, in this study, we identified the relative controls of seasonal variability in ET, and how these controls vary among ecosystems. We used overlapping AmeriFlux and PhenoCam time series at a daily timestep from 20 sites to explore these linkages (# site-years >100), and our study area covered a broad climatological aridity gradient in the U.S. and Canada. We focused on disentangling the most important controls of bulk surface conductance (G s ) and evaporative fraction (EF = LE/[H + LE]), where LE and H represent latent and sensible heat fluxes, respectively. Specifically, we investigated how vegetation phenology varied in importance relative to meteorological variables (vapor pressure deficit and antecedent precipitation) as a driver of G s and EF using path analysis, a framework for quantifying and comparing the causal linkages among multiple response and explanatory variables. Our results revealed that the drivers of G s and EF seasonality varied significantly between energy- and water-limited ecosystems. Specifically, precipitation had a much higher effect in water-limited ecosystems, while seasonal patterns in canopy greenness emerged as a stronger control in energy-limited ecosystems. Given that phenology is expected to shift under future climate, our findings provide key information for understanding and predicting how phenology may impact 21st-century hydroclimate regimes and the surface-energy balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Drivers of Decadal Carbon Fluxes Across Temperate Ecosystems

Long-running eddy covariance flux towers provide insights into how the terrestrial carbon cycle operates over multiple timescales. Here, we evaluated variation in net ecosystem exchange (NEE) of carbon dioxide (CO 2 ) across the Chequamegon Ecosystem-Atmosphere Study AmeriFlux core site cluster in the upper Great Lakes region of the USA from 1997 to 2020. The tower network included two mature hardwood forests with differing management regimes (US-WCr and US-Syv), two fen wetlands with varying levels of canopy sheltering and vegetation (US-Los and US-ALQ), and a very tall (400 m) landscape-level tower (US-PFa). Together, they provided over 70 site-years of observations. The 19-tower Chequamegon Heterogenous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 campaign centered around US-PFa provided additional information on the spatial variation of NEE. Decadal variability was present in all long-term sites, but cross-site coherence in interannual NEE in the earlier part of the record became weaker with time as non-climatic factors such as local disturbances likely dominated flux time series. Average decadal NEE at the tall tower transitioned from carbon source to sink to near neutral over 24 years. Respiration had a greater effect than photosynthesis on driving variations in NEE at all sites. Declining snowfall offset potential increases in assimilation from warmer springs, as less-insulated soils delayed start of spring green-up. Higher CO 2 increased maximum net assimilation parameters but not total gross primary productivity. Stand-scale sites were larger net sinks than the landscape tower. Clustered, long-term carbon flux observations provide value for understanding the diverse links between carbon and climate and the challenges of upscaling these responses across space.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ground Heat Flux Reconstruction Using Bayesian Uncertainty Quantification Machinery and Surrogate Modeling

Ground heat flux (G 0 ) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G 0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G 0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G 0 data (measured or modeled) for snow-free period as a reference. When observed G 0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G 0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

54 ENVIRONMENTAL SCIENCES↗