Search NASA⌕ Search

SEARCH · Search NASA

Results for “Representativity”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cluster analysis of occupancy schedules in residential buildings in the United States

The energy performance of residential buildings significantly depends on the building occupants’ behavior, which can be highly variable. When the heating, ventilation and air conditioning (HVAC) system is controlled based on the presence or absence of occupants in a building, occupant behavior is of even further importance to its energy performance. In current practice, building energy simulation tools generally use a single occupancy profile to represent the building’s occupancy schedule, the schedule of which is considered to be the same, regardless of the type of household being modeled. Thus, there is significant potential for improvement to allow for more flexibility and accuracy in calculation of occupancy. The objective of this study is to assess the variations in the typical types of occupancy schedules followed by the U.S. population using cluster analysis. American Time Use Survey data, which statically represents the overall U.S. population’s activities, across 12 years (2006–2017), is used. The ATUS data is segregated into smaller groups based on age and weekday/weekend, then divided into activities that are considered “at home” and “away from home”, which are mapped to the presence or non-presence of occupants in the home. Cluster analysis is then used to identify common types of occupancy schedule patterns for each age group. Three main types of patterns are obtained from cluster analysis for each age group, which together represent approximately 88% of people in the United States. The output of the cluster analysis is further analyzed to evaluate the variation in characteristics, including the number of times leaving home, time of day when leaving the home, and the timespan of absence from the home. The results of this study provide detailed insights on how typical occupants in the United States spend their time in residential spaces which can be used to create occupancy profiles for residential buildings. Finally, these occupancy profiles could be utilized inform an assessment of the energy use impact of occupancy-based controls of energy consuming systems and technologies.

42 ENGINEERING↗
DOE OSTI2020

The timing of lunar solidification and mantle overturn recorded in ferroan anorthosite 62237

Ferroan anorthosite suite (FAS) rocks are widely interpreted to represent primordial lunar crust. Despite their importance in pinpointing the timing of lunar crust formation, robust chronological investigations for this rock type are scarce. Here, we report the Ar-Ar, Rb-Sr, and Sm-Nd isotopic systematics for the FAS troctolitic anorthosite 62237. The Ar-Ar isotopic system has been reset by a thermal event at 3710 ± 48 Ma, and the Rb-Sr isotopic systematics has been disturbed such that a Rb-Sr isochron age cannot be determined. However, an internal isochron for the Sm-Nd isotopic system has yielded an age of 4350 ± 73 Ma (MSWD = 2.0) with an initial ε 143 Nd CHUR of -0.53 ± 0.26. The mineral and whole-rock fractions of 62237 plot on the same internal isochron as FAS sample 60025. The combined datasets define an age of 4372 ± 35 Ma (MSWD = 4.0) with an initial ε 143 Nd CHUR of -0.17 ± 0.22. Literature Sm-Nd data for FAS and Mg-suite whole-rocks also plot on the 60025-62237 isochron. The coherence of data from both FAS and Mg-suite rocks examined thus far suggests that both rock suites formed contemporaneously from identical, or nearly identical, sources. In addition, the concordance of FAS and Mg-suite ages suggests that primordial crust solidification either involved both magmatic suites, or that Mg-suite magmatism was contemporaneous with FAS magmatism within resolution of the Sm-Nd chronometer. The ages for FAS and Mg-suite also coincide with the formation ages of the mare basalt source regions and urKREEP. Ferroan anorthosite suite rocks and urKREEP are thought to represent primordial LMO solidification products, whereas Mg-suite and the mare basalt source regions are argued to represent mixtures of various LMO crystallization products that were formed during density-driven overturn of the LMO. The concordance of ages implies that the 4372 ± 35 Ma Sm-Nd isochron records the age of mantle overturn, and that overturn occurred during, or shortly after, solidification of the LMO.

58 GEOSCIENCES↗
DOE OSTI2020

Forensic species identification of elasmobranchs landed in Costa Rican artisanal fisheries

Elasmobranchs represent one of the most imperiled taxa, yet few nations have developed catch limits and no international or bilateral catch limits exist. One barrier to establishing catch limits is a lack of accurate species-specific extraction rates due to many species looking similar, distinguishing characteristics (fins and head) being removed, or grouping of all elasmobranchs together in fisheries data. To better understand the species-specific impacts fisheries are having on elasmobranchs we used forensic genetic techniques to identify the types and quantities of shark species landed in Costa Rican artisanal small-scale and mid-scale fisheries. We collected elasmobranch tissue samples from fishermen in three locations along Costa Rica’s Pacific Coast: Coyote, Bejuco, and Ojochal, from April 2013 to September 2014. We used DNA barcoding techniques, utilizing a 1044 base pair region of the NADH dehydrogenase subunit 2 gene, to identify to species 275 of the 416 samples collected. Within these fisheries, seven species of shark and one ray were identified, with the scalloped hammerhead shark (Sphyrna lewini) accounting for ~750-80% of all sharks landed. Observed total lengths (TL) for the scalloped hammerhead shark (10.9–99 cm) in these fisheries suggests the sharks sampled were all either juveniles or neonates. Based on these data it is likely that the areas from which these samples were collected represent nursery grounds for the scalloped hammerhead shark. The quantity of young scalloped hammerhead sharks removed from these areas, and ones similar to them, by these artisanal fisheries could be of great importance when determining future conservation measures for this species. We conclude our findings also indicate that sharks caught in these artisanal fisheries do not appear to represent those sold in markets around San Jose.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Constrained Turret Defense with Fixed Final Time

In this paper, we extend existing turret defense differential game formulations involving a turn-constrained turret and mobile agent to include specified final time and a constraint. For the purposes of this analysis, the specified final time may represent some exogenous input, perhaps representing the time at which some other event will take place. As for the constraint, it represents a no-fly zone for the mobile agent. The scenario is formulated as a two-player, zero-sum differential game and solved via the method of characteristics (i.e., back-propagation of equilibrium trajectories). Three different trajectory types make up the solution: trajectories that end with the turret aligned with the mobile agent, trajectories that end with the mobile agent on the constraint boundary, and regular trajectories.

differential game, game theory↗
DOE OSTI2020

Equilibrium vapor pressure properties for absorbent and adsorbent materials

Sorption devices are important tools for the efficient utilization of fuels and waste heat. Amid a tremendous diversity of cycles and applications, all sorption systems have an equilibrium vapor pressure that depends on the sorbent temperature and composition. The vapor pressure properties of working fluids are reported in the literature in a variety of ways, which impedes wide-ranging cross comparisons or screening studies for novel applications. This work compiles equilibrium vapor pressure properties for 123 liquid absorbents with 31 absorbates and 139 solid adsorbents with 16 adsorbates. The adsorption pairs are represented with six functional forms. Additionally, most of the absorption pairs are represented with 10 functional forms, plus several that are represented with custom empirical equations. Because the functional forms used in the literature vary widely, in this work each functional form was generalized. This paper is designed to facilitate comparisons of working fluid properties for screening studies, provide a quick reference to existing research, and present a framework for standardizing the reportage of vapor pressure equilibrium data for existing and future working pairs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Proper orthogonal descriptors for efficient and accurate interatomic potentials

Here, we present the proper orthogonal descriptors for efficient and accuracy representation of the potential energy surface. The potential energy surface is represented as a many-body expansion of parametrized potentials in which the potentials are functions of atom positions and parameters. The proper orthogonal decomposition is employed to decompose the parametrized potentials into a set of proper orthogonal descriptors (PODs). Because of the rapid convergence of the proper orthogonal decomposition, relevant snapshots can be sampled exhaustively to represent the atomic neighborhood environment accurately with a small number of descriptors. The proper orthogonal descriptors are used to develop interatomic potentials by using a linear expansion of the descriptors and determining the expansion coefficients from a weighted least-squares regression against a density functional theory (DFT) training set. We present a comprehensive evaluation of the POD potentials on previously published DFT data sets comprising Li, Mo, Cu, Ni, Si, Ge, and Ta elements. The data sets represent a diverse pool of metals, transition metals, and semiconductors. The accuracy of the POD potentials are comparable to that of state-of-the-art machine learning potentials such as the spectral neighbor analysis potential (SNAP) and the atomic cluster expansion (ACE).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Feature engineering for machine learning enabled early prediction of battery lifetime

Accurate battery lifetime estimates enable accelerated design of novel battery materials and determination of optimal use protocols for longevity in deployments. Unfortunately, traditional battery testing may take years to reach thousands of cycles. Recent studies have shown that machine learning (ML) tools can predict lithium-ion battery lifetimes from 100 or fewer preliminary cycles, representing only a few weeks of cycling. Until now, conclusions about the efficacy and broad applicability of these predictions across a variety of cathode chemistries have been limited by available experimental information. In this work, we leverage a battery cycling dataset representing six cathode chemistries (NMC111, NMC532, NMC622, NMC811, HE5050, and 5Vspinel), multiple electrolyte/anode compositions, and 300 total carefully prepared pouch batteries to explore feature selection and battery chemistry's role in ML battery lifetime predictions. Here, a mean absolute error (MAE) of 78 cycles in prediction was seen for a chemistry-spanning test set from 100 preliminary cycles. Furthermore, an MAE of 103 cycles was seen when using only the first cycle. This study represents an in-depth investigation of strategies for feature selection for battery lifetime prediction, ML models' generalization across multiple battery chemistries, and predictions beyond the training set in the chemical space.

25 ENERGY STORAGE↗
DOE OSTI2021

Material Interactions in Severe Accidents – Benchmarking the MELCOR V2.2 Eutectics Model for a BWR-3 Mark-I Station Blackout: Part II – Uncertainty Analysis

Single case comparisons between severe accident simulations can provide detailed insights into severe accident model behavior, however, they cannot offer insights into model uncertainty, sensitivity to uncertain parameters, or underlying model biases.Here in this analysis, the single case benchmark comparison of the MELCOR material interaction models for a station blackout (SBO) scenario of a boiling water reactor (BWR) using representative Fukushima Daiichi Unit 1 boundary conditions is expanded to include an uncertainty analysis. As part of this uncertainty analysis, 1200 simulations are performed for each material interaction model (2400 total), with random sampling of 14 uncertain MELCOR input parameters. Input parameters are selected for their impact on models representing core degradation processes. These include candling, fuel rod failure, debris quenching and dryout. The analysis performed here is not a traditional “best-estimate” uncertainty analysis that uses best-estimate parameters or identifies best-estimate figure of merit distributions. Instead, it is an exploratory uncertainty analysis that identifies and interrogates underlying model form biases of the two material interaction models (eutectics and interactive materials models). Uniform distributions are applied to all uncertain parameters to ensure coverage of the model parameter uncertainty space. Key findings from this study include underlying model form biases exhibited by material interaction models, and notable differences in accident progression outcomes between the material interaction models. This uncertainty study extends and confirms the conclusions from the first part of this study, which compared the impact of material interaction modeling on simulation of a short-term station blackout scenario with representative Fukushima Daiichi Unit I boundary conditions. In particular, this study confirms that the eutectics model generally exhibits accelerated degradation and failure of fuel components, the core plate, and the lower head. The eutectics model also has a tendency to exhibit a greater degree of core degradation, greater debris mass formation, and larger debris mass ejection. Finally, the eutectics model exhibits higher maximum temperatures for fuel, cladding, particulate debris, oxidic molten pool, and metallic molten pool components than the interactive materials model; interactive materials model simulations exhibit a soft “limitation” on maximum temperatures that is related to the temperature at which material relocation occurs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Scaling methodologies and similarity analysis for thermal hydraulics test facility development for water-cooled small modular reactor

Small modular reactors (SMRs) represent a promising option for providing clean and sustainable energy due to their potential for enhanced safety, reduced capital costs, and increased siting flexibility. However, new reactor systems require the development and operation of representative scaled-down test facilities to support the verification and validation of system computer codes and models. Here this study reviews the research on scaling methodologies and similarity principles pivotal in developing non-nuclear integral effects test and separate effects test facilities for water-cooled SMRs. The study focuses on a review of the scaling methods, similarity approaches, and possible challenges posed by the unique and compact design features of integral-pressurized water reactor-type SMRs, and their representative test facilities. This study also reviews previous research related to scaling and similarity methodologies and provides insights into design considerations for achieving prototypic conditions in test facilities. The findings and recommendations emphasize the broader impact of appropriate scaling and similarity principles to ensure meaningful and transferable results from non-nuclear test facilities to accelerate the safe and efficient deployment of next-generation water-cooled SMRs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Effect of electrode rinse solutions on the electrodialysis of concentrated salts

The data presented indicate that modest attention to the chemistry of the electrode rinse solution can yield as much as 56% improvement in process efficiency. When operating electrodialysis (ED) with salt concentrations greater than 10,000 mg/l (as NaCl) the amperage can be relatively high even at low voltages. To maximize the overall rate of ion transport from diluate to concentrate, there is a need to minimize resistances in the electrode cells since these can represent as much as 30% of the total resistances in the entire process and in this particular ED unit, represented resistance roughly equivalent to that within the membrane stack. This work describes the performance of a 10-cell pair, 200 cm 2 (0.02 m 2 ) per membrane, pilot ED unit operated in batch mode at 5 V potential. The feedstock to the stack was varied from 0.5% to 6% NaCl. Results from Volt-amp profiles (2.0–15 V) were used as the rationale for choosing the preferred electrode rinse solution, 90 g/l (kg/m 3 ) disodium sulfate at pH 12.5. The recommended solution of 30 g/l (kg/m 3 ) disodium sulfate at neutral pH was tested against stronger solutions (60, 90, and 120 g/l (kg/m 3 ) disodium sulfate) at neutral pH and with 1 g/l (kg/m 3 ) sodium hydroxide added to yield solutions at pH 12.5. There were stark differences in the performance of the ED unit between the solutions at pH 7 and those at pH 12.5. The slopes of the Volt-amp profiles improved indicating less resistance to ion flow at pH 12.5. Most of the difference can be attributed to the greater conductivity of sodium hydroxide compared to disodium sulfate, where even a modest addition of 1 g/l (kg/m 3 ) sodium hydroxide increased the conductivity of the solution between 8 and 20% depending on the concentration of disodium sulfate. Most notably, the initiation voltage was 2.39 ± 0.08 V with solutions at pH 7 and 1.99 V ± 0.04 V when the rinse solutions were around pH 12.5, a difference of around 0.42 V. Since the full batch runs were performed at 5.0 V, this 0.42 V differential represented a loss of 8% in process efficiency. A series of tests showed a shift in initiation voltage occurred between pH 11.5 and 12.5. Nernst equations for water electrolysis are coupled with a simple flux model that estimates the concentration of hydroxide and hydronium ions at the electrode surfaces. This model indicates that the potential demanded by the electrodes is directly linked to the current. Furthermore, there is a steep drop in the potential demanded by the anode when there is sufficient hydroxide available to generate oxygen from hydroxide as opposed to generation from water. The voltage demand is minimal at pH 12.5 and reflective of the reference standard voltage (just over 1.23 V predicted at 0.1 amps). However, the flux model predicts that the voltage demanded by the electrodes at pH 7 is 1.61 V at 0.1 amps (the lowest amperage measurable with the power pack). In conclusion, comparing this result to the reference standard (E 0 = -1.23 V) yields a satisfactory explanation for the observed difference below reference standard of 0.42 V.

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

Revising the dynamic energy budget theory with a new reserve mobilization rule and three example applications to bacterial growth

Dynamic energy budget (DEB) theory has been applied to model a wide range of organisms, including microbes. In the standard DEB model, biomass is partitioned into reserve and structural compartments, where reserve biomass is mobilized in a pseudolinear manner (while the reserve biomass density, defined as the ratio between reserve and structural biomass, decays linearly) to drive maintenance and the growth of structural biomass (and extracellular enzyme production if it is considered). However, the linear dynamics of the reserve biomass density makes the standard DEB model incapable of explaining the slowdown of microbial growth at high reserve density that is caused by macromolecular crowding effect which reduces biochemical reaction rates (a typical situation occurs when microbes are experiencing severe moisture stress) and is inconsistent with the observation that intracellular enzymatic reactions generally follow non-linear kinetics. By partitioning biomass into reserve, kinetic, and structural compartments, we show here that the Equilibrium Chemistry Approximation (ECA) kinetics can be used to represent enzymatically catalyzed reserve biomass mobilization that can then drive the kinetic and structural biomass synthesis. This revised DEB model better represents the tradeoff in ribosome allocation for structural growth and internal enzyme production, is structurally compatible with metabolic models of cell individuals, and includes the standard DEB model and the popular compromise model as special cases for representing population growth. We then applied the revised DEB model to interpret components of bacterial respiration, their dependence on substrate availability, and emergent microbial carbon use efficiency dynamics for an exponentially growing population. We found that the revised DEB model enables a better understanding of bacterial substrates use (carbon in our examples) than that can be derived from a few other models in the literature. In particular, the revised DEB model explains why carbon use efficiency may first increase, then plateau, and finally decrease with growth rate (and substrate uptake rate), as a function of proteomics. Additionally, the revised DEB model explains why the kinetic biomass compartment needs to be divided to reasonably incorporate proteomic control of microbial growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2026

Ultrafast population and structural dynamics of a Ni-bipyridine photoredox catalyst reveal a significant deactivation pathway

The ultrafast excited state pathways and dynamics of NiII-bipyridine complexes influence the yield of photochemical processes involved in their catalytic cross-coupling reactions. Here we present ultrafast Ni K x-ray emission spectroscopy (XES) and x-ray solution scattering (XSS) of a NiII-bipyridine aryl halide complex, [Ni(t-Bubpy)(o-tol)Br], to quantify the excited state population dynamics and structural changes of the pre-catalyst. Due to the local spin-sensitivity of XES, the population dynamics of metal-to-ligand charge transfer (MLCT) and metal-centered (MC) excited states is established. A rapid ground state recovery pathway is newly identified, representing a significant deactivation pathway during photocatalysis. Furthermore, the pseudotetrahedral structure of the long-lived MC excited state is unambiguously identified and refined by XSS. The results advance our understanding of the ultrafast relaxation mechanisms that impact the photocatalytic mechanism and yield for NiII-bipyridine aryl halide cross-coupling catalysts.

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

Chemical Differences between Phenolic Secondary Organic Aerosol Formed through Gas-Phase and Aqueous-Phase Reactions

Phenolic compounds, which are significant emissions from biomass burning (BB), undergo rapid photochemical reactions in both gas and aqueous phases to form secondary organic aerosol, namely, gasSOA and aqSOA, respectively. The formation of gasSOA and aqSOA involves different reaction mechanisms, leading to different product distributions. In this study, we investigate the gaseous and aqueous reactions of guaiacol a representative BB phenol to elucidate the compositional differences between phenolic aqSOA and gasSOA. Aqueous-phase reactions of guaiacol produce higher SOA yields than gas-phase reactions (e.g., roughly 60 vs 30% at one half-life of guaiacol). These aqueous reactions involve more complex reaction mechanisms and exhibit a more gradual SOA evolution than their gaseous counterparts. Initially, gasSOA forms with high oxidation levels (O/C > 0.82), while aqSOA starts with lower O/C (0.55–0.75). However, prolonged aqueous-phase reactions substantially increase the oxidation state of aqSOA, making its bulk chemical composition closer to that of gasSOA. Additionally, aqueous reactions form a greater abundance of oligomers and high-molecular-weight compounds, alongside a more sustained production of carboxylic acids. AMS spectral signatures representative of phenolic gasSOA have been identified, which, together with tracer ions of aqSOA, can aid in the interpretation of field observation data on aerosol aging within BB smoke. The notable chemical differences between phenolic gasSOA and aqSOA highlighted in this study also underscore the importance of accurately representing both pathways in atmospheric models to better predict the aerosol properties and their environmental impacts.

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

Mathematical Reconstruction of Land Carbon Models From Their Numerical Output: Computing Soil Radiocarbon From C Dynamics

Radiocarbon ( 14 C) is a powerful tracer of the global carbon cycle that is commonly used to assess carbon cycling rates in various Earth system reservoirs and as a benchmark to assess model performance. Therefore, it has been recommended that Earth System Models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 report predicted radiocarbon values for relevant carbon pools. However, a detailed representation of radiocarbon dynamics may be an impractical burden on model developers. Here, we present an alternative approach to compute radiocarbon values from the numerical output of an ESM that does not explicitly represent these dynamics. The approach requires computed 12 C stocks and fluxes among all carbon pools for a particular simulation of the model. From this output, a time-dependent linear compartmental system is computed with its respective state-transition matrix. Using transient atmospheric 14 C values as inputs, the state-transition matrix is then applied to compute radiocarbon values for each pool, the average value for the entire system, and component fluxes. We demonstrate the approach with ELMv1-ECA, the land component of an ESM model that explicitly represents 12 C, and 14 C in 7 soil pools and 10 vertical layers. Results from our proposed method are highly accurate (relative error <0.01%) compared with the ELMv1-ECA 12 C and 14 C predictions, demonstrating the potential to use this approach in CMIP6 and other model simulations that do not explicitly represent 14 C.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Incorporating Microtopography in a Land Surface Model and Quantifying the Effect on the Carbon Cycle

Abstract Northern peatlands are a terrestrial carbon store, with an annual sink of 0.1 Pg C yr −1 and a total storage estimate of 547 Pg C. Northern peatlands are also major contributors of atmospheric methane. Most land surface models do not accurately represent peatland carbon emissions, partly because they do not represent the hydrologic cycle and/or microtopography adequately. Interactions between water table depth and microtopography in peatlands influence decomposition and modulate CO 2 and CH 4 fluxes. A modified version of the land surface component of the Energy Exascale Earth System Model, was recently created to represent the microtopography and hydrology of a raised dome bog in northern Minnesota, USA. In this study, three microtopographic parameters are analyzed in the modified version: hummock height, hummock‐hollow spacing, and percent hollow. Terrestrial laser scanning observations are used to set uncertainty bounds for these parameters. Our model experiment results suggest that carbon‐related quantities of interest (QOI) were typically the most sensitive to hummock height, and those QOI (especially net ecosystem exchange, NEE) were sensitive to interactions between parameters. Furthermore, NEE was most relatively influenced by microtopographic parameters in the model, varying by 35%. We found that increasing hummock height resulted in more C being stored in plant tissue and less in soil organic matter. This coincided with decreases in Sphagnum and increases in Picea and shrub net primary production. These results suggest that future studies may consider extending prognostic capabilities of carbon cycling by incorporating hummock hollow microtopography into earth system models.

58 GEOSCIENCES↗
DOE OSTI2022

A Decentralized Approach for Modeling Organized Convection Based on Thermal Populations on Microgrids

Abstract In this study, a spectral model for convective transport is coupled to a thermal population model on a two‐dimensional horizontal “microgrid,” covering the typical gridbox size of general circulation models. The goal is to explore new ways of representing impacts of spatial organization in cumulus cloud fields. The thermals are considered the smallest building block of convection, with thermal life cycle and movement represented through binomial functions. Thermals interact through two simple rules, reflecting pulsating growth and environmental deformation. Long‐lived thermal clusters thus form on the microgrid, exhibiting scale growth and spacing that represent simple forms of spatial organization and memory. Size distributions of cluster number are diagnosed from the microgrid through an online clustering algorithm, and provided as input to a spectral multiplume eddy‐diffusivity mass flux scheme. This yields a decentralized transport system, in that the thermal clusters acting as independent but interacting nodes that carry information about spatial structure. The main objectives of this study are (a) to seek proof of concept of this approach, and (b) to gain insight into impacts of spatial organization on convective transport. Single‐column model experiments demonstrate satisfactory skill in reproducing two observed cases of continental shallow convection. Metrics expressing self‐organization and spatial organization match well with large‐eddy simulation results. We find that in this coupled system, spatial organization impacts convective transport primarily through the scale break in the size distribution of cluster number. The rooting of saturated plumes in the subcloud mixed layer plays a key role in this process.

54 ENVIRONMENTAL SCIENCES↗