Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “R”

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

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

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

SDSS-IV MaNGA: The kinematic-morphology of galaxies on the mass versus star-formation relation in different environments

ABSTRACT We study the link between the kinematic-morphology of galaxies, as inferred from integral-field stellar kinematics, and their relation between mass and star formation rate. Our sample consists of ∼3200 galaxies with integral-field spectroscopic data from the MaNGA survey (Mapping Nearby Galaxies at Apache Point Observatory) with available determinations of their effective stellar angular momentum within the half-light radius $\lambda _{R_e}$. We find that for star-forming galaxies, namely along the star formation main sequence (SFMS), the $\lambda _{R_e}$ values remain large and almost unchanged over about two orders of magnitude in stellar mass, with the exception of the lowest masses $\mathcal {M}_{\star }\lesssim 2\times 10^{9} \, \mathcal {M}_{\odot }$, where $\lambda _{R_e}$ slightly decreases. The SFMS is dominated by spiral galaxies with small bulges. Below the SFMS, but above the characteristic stellar mass $\mathcal {M}_{\rm crit}\approx 2\times 10^{11} \, \mathcal {M}_{\odot }$, there is a sharp decrease in $\lambda _{R_e}$ with decreasing star formation rate (SFR): massive galaxies well below the SFMS are mainly slow-rotator early-type galaxies, namely genuinely spheroidal galaxies without discs. Below the SFMS and below $\mathcal {M}_{\rm crit}$ the decrease of $\lambda _{R_e}$ with decreasing SFR becomes modest or nearly absent: low-mass galaxies well below the SFMS, are fast-rotator early-type galaxies, and contain fast-rotating stellar discs like their star-forming counterparts. We also find a small but clear environmental dependence for the massive galaxies: in the mass range $10^{10.9}\!-\!10^{11.5} \, \mathcal {M}_{\odot }$, galaxies in rich groups or denser regions or classified as central galaxies have lower values of $\lambda _{R_e}$. While no environmental dependence is found for galaxies of lower mass. We discuss how the above results can be understood as due to the different star formation and mass assembly histories of galaxies with varying mass.

Wang, Bitao↗
DOE OSTI2024

Statistical properties of filaments in the cosmic web

ABSTRACT In the context of the cosmological and constrained Exploring the Local Universe with the reConstructed Initial Density field (ELUCID) simulation, this study explores the statistical characteristics of filaments within the cosmic web, focussing on aspects such as the distribution of filament lengths and their radial density profiles. Using the classification of the cosmic web environment through the Hessian matrix of the density field, our primary focus is on how cosmic structures react to the two variables $R_{\rm s}$ and $\lambda _{\rm th}$. The findings show that the volume fractions of knots, filaments, sheets, and voids are highly influenced by the threshold parameter $\lambda _{\rm th}$, with only a slight influence from the smoothing length $R_{\rm s}$. The central axis of the cylindrical filament is pinpointed using the medial-axis thinning algorithm of the COsmic Web Skeleton (COWS) method. It is observed that median filament lengths tend to increase as the smoothing lengths increase. Analysis of filament length functions at different values of $R_{\rm s}$ indicates a reduction in shorter filaments and an increase in longer filaments as $R_{\rm s}$ increases, peaking around $2.5R_{\rm s}$. The study also shows that the radial density profiles of filaments are markedly affected by the parameters $R_{\rm s}$ and $\lambda _{\rm th}$, showing a valley at approximately $2R_{\rm s}$, with increases in the threshold leading to higher amplitudes of the density profile. Moreover, shorter filaments tend to have denser profiles than their longer counterparts.

Zhang, Youcai (ORCID:0000000319674091)↗
DOE OSTI2022

Interannual consistency of migration phenology is season- and breeding region-specific in North American Golden Eagles

Abstract Interannual consistency (an indicator of the strength of adjustments) in migration phenology of Golden Eagles (Aquila chrysaetos) in North America is most strongly associated with the breeding region, the season, and with late-season temperature on breeding and wintering grounds. Consistency was greatest in boreal spring migration and the breeding regions of eastern Canada. Using multi-year GPS tracks of 83 adults breeding in 3 spatially distant regions (Alaska, northeast Canada, and southeast Canada), we quantified the interannual consistency of migration phenology and wintering latitude within and among individuals tracked across multiple years and the repeatability (r) by breeding regions and seasons. By comparing regions and seasons, we found that consistency was highest (r > 0.85) for boreal spring migration in eastern Canada while Alaska had the lowest value (r < 0.15). Because seasonal consistency of migration phenology was only detected in eastern Canada, we conclude that seasonal features are not a primary constraint. While regional differences in consistency were not related to differences in migratory distances, they could be the result of genetic or habitat differences. We also found that temperatures warmer than the decadal average at the region of departure delayed the start of boreal spring migration by ~10 days and advanced boreal autumn migration by ~20 days. These results suggest that warmer temperatures would reduce residence time on breeding grounds, contrary to expectations and trends found in other studies. Wide variations in migratory strategies across a species distribution can add to the list of challenges for conservation but may give migrants the capacity to acclimate to environmental changes.

Zoology↗
DOE OSTI2021

Electronic scattering off a magnetic hopfion

We study scattering of itinerant electrons off a magnetic hopfion in a three-dimensional metallic magnet described by a magnetization vector S(r). A hopfion is a confined topological soliton of S(r) characterized by an emergent magnetic field Bγ(r)≡ε αβγ S∙(∇ α S×∇ β S)/4≠0 with vanishing average value [B(r)]=0. We evaluate the scattering amplitude in the opposite limits of large and small hopfion radius R using the eikonal and Born approximations, respectively. In both limits, we find that the scattering cross section contains a skew-scattering component giving rise to the Hall effect within a hopfion plane. That conclusion contests the popular notion that the Hall effect in noncollinear magnetic structures necessarily implies [B(r)]≠0. In the limit of small hopfion radius pR<<1, we expand the Born series in powers of momentum p and identify different expansion terms corresponding to the hopfion anisotropy, toroidal moment, and skew scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Charge order and superconductivity in a two-band model for infinite-layer nickelates

The recent discovery of superconductivity in infinite-layer nickelates has drawn considerable attention; however, a consensus on the fundamental building blocks and common ingredients necessary to understand and describe their ground states and emergent properties is lacking. A series of experimental and theoretical studies have suggested that an effective two-band Hubbard model with Ni 3d x 2 -y 2 and rare-earth (R) 5d character may describe the low-energy physics. In this report we study the ground state properties of this two-band model on four-leg cylinders using the density-matrix renormalization group (DMRG) technique to better grasp whether such a simple model can embody the essential physics. A key difference compared to single-band Hubbard materials is that the system is self-doped: even at overall half-filling, the R-band acts as an electron reservoir, hole-doping the Ni-layer and fundamentally altering the physics expected from an undoped antiferromagnet. On the four-leg cylinder, the ground state is consistent with a Luttinger liquid, with anti-phase modulations of the charge density in the Ni- and R-layers having corresponding wavevectors that lock together. Light hole doping away from 1/2 filling releases the locking between the Ni and the R charge modulations, as the electron density in the R-band decreases and eventually becomes exhausted at a hole doping concentration that depends sensitively on the effective splitting between the Ni and the R orbitals. The ground state of the doped system is consistent with a Luther-Emery liquid, possessing quasi-long-range superconducting correlations in the Ni layer, similar to the single-band Hubbard model. Our results are consistent with experimental observations and may help to reveal the microscopic mechanism for pairing and other emergent properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamical defects in a two-dimensional Wigner crystal: Self-doping and kinetic magnetism

We study the quantum dynamics of interstitials and vacancies in a two-dimensional Wigner crystal (WC) using a semi-classical instanton method that is asymptotically exact at low density, i.e., in the r s → ∞ limit. Here, the dynamics of these point defects mediates magnetism with much higher energy scales than the exchange energies of the pure WC. Via exact diagonalization of the derived effective Hamiltonians in the single-defect sectors, we find the dynamical corrections to the defect energies. The resulting expression for the interstitial (vacancy) energy extrapolates to 0 at r s = r mit ≈ 70 (r s ≈ 30), suggestive of a self-doping instability to a partially melted WC for some range of r s below r mit . We thus propose a “metallic electron crystal” phase of the two-dimensional electron gas at intermediate densities between a low density insulating WC and a high density Fermi fluid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamical model of J/ψ photoproduction on the nucleon

Here, a dynamical model based on a phenomenological charm quark-nucleon (c – N) potential v cN and the Pomeronexchange mechanism is constructed to investigate the J/ψ photoproduction on the nucleon from threshold to invariant mass W = 300 GeV. The J/ψ – N potential, V J/ψN (r), is constructed by folding v cN into the wave function φ J/ψ ($c\bar{c}$) of J/ψ within a constituent quark model (CQM) of Segovia et al. [Int. J. Mod. Phys. E 22, 1330026 (2013)]. A photoproduction amplitude is also generated by v cN by a $c\bar{c}$–loop integration over the γ → $c\bar{c}$ vertex function and φ J/ψ ($c\bar{c}$). No commonly used vector meson dominance assumption is used to define this photoproduction amplitude which is needed to describe the data near the threshold. The c – N potential v cN (r) is parameterized in a form such that the predicted V J/ψN (r) at large distances has the same Yukawa potential form extracted from a lattice QCD (LQCD) calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)]. The parameters of vcN are determined by fitting the total cross-section data of Jefferson Laboratory (JLab) by performing calculations that include J/ψ – N final-state interactions (FSI). The resulting differential cross sections dσ /dt are found in good agreements with the data. It is shown that the FSI effects dominate the cross section in the very near-threshold region, allowing for sensitive testing of the predicted J/ψ – N scattering amplitudes. By imposing the constraints of J/ψ – N potential extracted from the LQCD calculation of Kawanai and Sasaki, [Phys. Rev. D 82, 091501(R) (2010)], we have obtained three J/ψ – N potentials which fit the JLab data equally well. The resulting J/ψ – N scattering lengths are in the range of a = [-0.05, -0.25] fm. With the determined v cN (r) and the wave functions generated from the same CQM, the constructed model is used to predict the cross sections of photoproduction of η c (1S) and ψ(2S) mesons for future experimental tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Hyper Suprime-Cam Year 3 results: Cosmology from galaxy clustering and weak lensing with HSC and SDSS using the minimal bias model

Here, we present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering about 416 deg 2 , and the SDSS DR11 spectroscopic galaxies spanning the redshift range [0.15, 0.70]. We subdivide the SDSS galaxies into three luminosity-cut, and therefore nearly volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, w p ⁡(R). We also use the measurements of the galaxy-galaxy weak-lensing signal Δ⁢Σ⁡(R) for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, ξ ± ⁡($\vartheta$) measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat Λ⁢CDM cosmology. We use conservative scale cuts, R >12 and 8⁢ℎ -1 Mpc for Δ⁢Σ and w p , respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. We present various validation tests of our model as well as analysis methods. Our baseline analysis yields S 8 = $0.77⁢5^{+0.043}_{-0.038}$ (68% C.I.) for the Λ⁢CDM model, after marginalizing over uncertainties in other parameters. Our value of S 8 is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. We show that various internal consistency tests based on different splits of the data are passed. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model, using Δ⁢Σ⁡(R) and w p ⁡(R) down to R >3 and 2⁢ℎ -1 Mpc, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Anomalous Nematic States in High Half-Filled Landau Levels

It is well established that the ground states of a two-dimensional electron gas with half-filled high ( N ≥ 2 ) Landau levels are compressible charge-ordered states, known as quantum Hall stripe (QHS) phases. The generic features of QHSs are a maximum (minimum) in a longitudinal resistance R x x ( R y y ) and a nonquantized Hall resistance R H . In this letter, we report on emergent minima (maxima) in R x x ( R y y ) and plateaulike features in R H in half-filled N ≥ 3 Landau levels. Remarkably, these unexpected features develop at temperatures considerably lower than the onset temperature of QHSs, suggestive of a new ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Low Thermal Resistance (0.5 K/W) Ga 2 O 3 Schottky Rectifiers With Double-Side Packaging

The low thermal conductivity of Ga 2 O 3 has arguably been the most serious concern for Ga 2 O 3 power and RF devices. Despite many simulation studies, there is no experimental report on the thermal resistance of a large-area, packaged Ga 2 O 3 device. This work fills this gap by demonstrating a 15-A double-side packaged Ga 2 O 3 Schottky barrier diode (SBD) and measuring its junction-to-case thermal resistance (R θJC ) in the bottom-side- and junction-side-cooling configurations. The R θJC characterization is based on the transient dual interface method, i.e., JEDEC 51-14 standard. The R θJC of the junction- and bottom-cooled Ga 2 O 3 SBD was measured to be 0.5 K/W and 1.43 K/W, respectively, with the former R θJC lower than that of similarly-rated commercial SiC SBDs. This low R θJC is attributable to the heat extraction directly from the Schottky junction instead of through the Ga 2 O 3 chip. Here, the R θJC lower than that of commercial SiC devices proves the viability of Ga 2 O 3 devices for high-power applications and manifest the significance of proper packaging for their thermal management.

42 ENGINEERING↗
DOE OSTI2021

On the use of air temperature and precipitation as surrogate predictors in soil respiration modelling

Soil respiration (R S ), the soil-to-atmosphere CO 2 flux that is a major component of the global carbon cycle, is strongly influenced by local soil temperature (T soil ) and water content (SWC). Regional to global-scale R S modelling thus requires this information at local scales, but few high-quality, wall-to-wall (global) T soil and SWC data exist. As a result, such modelling efforts commonly use air temperature (T air ) and monthly precipitation (P m ) as surrogate predictors, but their site-scale accuracy and potential bias are unknown. In this report we used monthly data from 880 sites across a wide variety of different environmental conditions (i.e., climate, ecosystem type, elevation, vegetation leaf habit and drainage conditions) to determine the suitability of T air as a surrogate for T soil , and data from 507 sites to examine the suitability of P m as a surrogate for SWC. Site-specific linear and second-order exponential non-linear models were compared using model evaluation metrics (i.e., slope, p-value of slope, root mean square error [RMSE], index of agreement and model efficiency). We found that T soil and T air are highly correlated and explain similar R S variability. In contrast, P m is not a good surrogate for SWC, even though P m explains a similar amount of R S variability to SWC. The wide variability in the site-specific relationships between R S and SWC means that no single relationship can be used for large-scale modelling. The results from this study support the use of T air in continental-to-global scale R S models, and highlight the urgent need for continental-to-global scale SWC datasets for the modelling and evaluation of future soil carbon dynamics under global climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Photoluminescence spectroscopy of Cr 3+ in β-Ga 2 O 3 and (Al 0.1 Ga 0.9 ) 2 O 3

Alloying β-Ga 2 O 3 with Al 2 O 3 to create (Al x Ga 1–x ) 2 O 3 enables ultra-wide bandgap materials suitable for applications deep into ultraviolet. In this work, photoluminescence (PL) spectra of Cr 3+ were investigated in monoclinic single crystal β-Ga 2 O 3 , and 10 = mol. % Al 2 O 3 alloyed with β-Ga 2 O 3 , denoted β-(Al 0.1 Ga 0.9 ) 2 O 3 or AGO. Temperature-dependent PL properties were studied for Cr 3+ in AGO and β-Ga 2 O 3 from 295 to 16 K. For both materials at room temperature, the red-line emission doublet R 1 and R 2 occurs at 696 nm (1.78 eV) and 690 nm (1.80 eV), respectively, along with a broad emission band at 709 nm (1.75 eV). Here, the linewidths for AGO are larger for all temperatures due to alloy broadening. For both materials, the R-lines blue-shift with decreasing temperature. The (lowest energy) R 1 line is dominant at low temperatures due to the thermal population of the levels. For temperatures above ~50 K, however, the ratio of R 2 to R 1 peak areas is dominated by nonradiative combination.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Maximum Size Intersecting Families of Bounded Minimum Positive Co-degree

Let $\mathcal{H}$ be an $r$-uniform hypergraph. The minimum positive co-degree of $\mathcal{H}$, denoted by $\delta_{r-1}^+(\mathcal{H})$, is the minimum $k$ such that if $S$ is an $(r-1)$-set contained in a hyperedge of $\mathcal{H}$, then $S$ is contained in at least $k$ hyperedges of $\mathcal{H}$. For $r\geq k$ fixed and $n$ sufficiently large, we determine the maximum possible size of an intersecting $r$-uniform $n$-vertex hypergraph with minimum positive co-degree $\delta_{r-1}^+(\mathcal{H}) \geq k$ and characterize the unique hypergraph attaining this maximum. This generalizes the Erd\Hos--Ko--Rado theorem which corresponds to the case $k=1$. Furthermore, our proof is based on the delta-system method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

What Controls the Entrainment Rate of Dry Buoyant Thermals with Varying Initial Aspect Ratio?

Here, this study uses theory and numerical simulations to analyze the nondimensional spreading rate a (change in radius with height) of buoyant thermals as they rise and entrain surrounding environmental fluid. A focus is on how a varies with initial thermal aspect ratio A r , defined as height divided by width of the initial buoyancy perturbation. An analytic equation for thermal ascent rate w t that depends on a is derived from the thermal-volume-averaged momentum bud get equation. The thermal top height when w t is maximum, defining a critical height z c , is inversely proportional to a. The height z c also corresponds to the thermal top height when buoyant fluid along the thermal’s vertical axis is fully replaced by entrained nonbuoyant environmental fluid rising from below the thermal. The time scale for this process is controlled by the vertical velocity of parcels rising upward through the thermal’s core. This parcel vertical velocity is approximated from Hill’s analytic spherical vortex, yielding an analytic inverse relation between a and A r . Physically, this a–A r relation is connected to changes in circulation as A r is modified. Numerical simulations of thermals with A r varied from 0.5 to 2 give a values close to the analytic theoretical relation, with a factor of ~3 decrease in a as A r is increased from 0.5 to 2. The theory also explains why a of initially spherical thermals from past laboratory and modeling studies is about 0.15. Overall, this study provides a theoretical underpinning for understanding the entrainment behavior of thermals, relevant to buoyantly driven atmospheric flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reproductive and leaf litterfall fluxes in forest ecosystem sites globally (1950-2022)

Forest allocation of net primary productivity (NPP) to reproduction is poorly quantified globally, despite its critical role in forest regeneration and a well-supported trade-off with allocation to growth. Although field measurements of total NPP are rare, our work finds that a proxy for reproductive carbon allocation constructed from leaf (L) and reproductive (R) litterfall fluxes, R/(R+L), is strongly correlated with R/NPP, facilitating analysis across a wide range of sites where biometric estimates of NPP are not available (R² = 0.85; Hanbury-Brown et al., 2022, Ward et al., in prep). To investigate relationships between ecosystem-scale reproductive allocation (RA) and climate, soil fertility, and stand age gradients, we conducted a literature search and synthesized 824 observations of annual average leaf and reproductive litterfall fluxes across forest sites globally. The zip file includes 1) a folder Data/ containing the litterfall data ("GlobalForestRA_data.csv") and metadata ("GlobalForestRA_metadata.doc") files. The data file includes geographic coordinates, long-term mean annual temperature and precipitation (1970-2000, extracted from WorldClim2.1), leaf and reproductive litterfall fluxes, sampling interval and protocols, forest characteristics (dominant leaf morphology, information pertaining to forest age and successional stage, and disturbance history) and soil properties (% sand, %silt, %clay, total phosphorus (P), nitrogen (N), cation exchange capacity (CEC) and pH) extracted from SoilGrids250 and from on-site measurements, where available. The metadata file contains information about each variable reported in the data file, including data sources, processing methods, and all references. The Data folder contains two additional files used to create Figure 1; these are described in greater detail in the README.2) R scripts GloalForestRA_analysis.r and GlobalForestRA_SI.r and a folder /Functions used to produce results, figures, and tables in the manuscript Ward et al. (in press)3) a README file describing how the data and R scripts can be used to reproduce statistical results, figures, and tables found in the manuscript. Ward et al. (in press)This repository can also be found at: https://github.com/r-ward/Global_Analysis_ForestRA.Ward, R.E., Zhang-Zheng, H. Aernethy, K., Adu-Bredu, S., Arroyo, L., Bailey, A. et al. (in press). Forest age rivals climate to explain reproductive allocation patterns in forest ecosystems globally. Ecology Letters. Hanbury-Brown, A.R., Ward, R.E. & Kueppers, L.M. (2022). Forest regeneration within Earth system models: current process representations and ways forward. New Phytol., 235, 20–40.Ward et al. (2025), Forest age rivals climate to explain reproductive allocation patterns in forest ecosystems globally, in prep.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Topsoil bulk geochemical compositions - An updated harmonized global dataset

Mineral weathering is a key biogeochemical process because of the capacity of minerals to stabilize organic matter. However, predicting soil weathering status across large spatial areas still isn’t possible due to a lack of global data and theoretical frameworks. To address this knowledge gap, multiple global datasets of bulk topsoil geochemical compositions have been harmonized using R. These datasets document topsoil bulk geochemical compositions across five continents (n = ~16,000 observations). Source data for these observations include the EuroGEOSurveys Geochemical Baseline Database (FOREGS), the US Geological Survey National Geochemical Database (NASGLP), the Geochemical Atlas of Australia (GAA), the US Geological Survey Alaska Geochemical Database (AGD84), the National Cooperative Soil Survey (NCSS), the European Geochemical Mapping of Agricultural Soil (GEMAS), Ecorespira-Amazon (ERA), the New Zealand Geochemical Baseline Survey (NZ_GBS), and the African Soil Information Service (AFSIS). Major elements observed include Aluminum (Al), Calcium (Ca), Iron (Fe), Potassium (K), Magnesium (Mg), Sodium (Na), Titanium (Ti), Manganese (Mn), Phosphorus (P), Carbon (C), and Sulfur (S). This data package includes the harmonized dataset itself, and the R scripts necessary to harmonize these datasets, in addition to metadata that describes all columns, files, and databases used in this project. Methods & Sampling Step 1 – Databases of geochemical data identified This study aimed to leverage existing measurements of topsoil geochemical data. Databases were first identified and deemed appropriate for inclusion if they were measuring soils and performed these measurements on the <2mm soil fraction. Databases such as NCSS and AGD84 needed more post processing to include in the database and this was done using the NCSS_datamerge_031626 R file and Alaska_USGSmerge_031626 R file, respectively. Step 2 – Database harmonization Once appropriate databases were identified, they were harmonized for ease of analysis using the R script Database_Harmonization_031826. This included removing columns from original datasets that would not be used in analysis (removed columns are noted in the code). Then, data cleaning procedures specific to each dataset were undertaken. This includes standardizing columns to include units and adding metadata columns regarding procedures for analyzing specific elements. Functions for standardizing measurements and units are outline in R files: calculate element_mg_kg_031626, calculate_oxide_wt_perc_031626, change_oxide_caps_031626, and conv_2_numeric_031626. This also included adding a unique identifier for each sample to identify it with its respective database (see CD_ID in data dictionary). Geographic information: Data reflect a compilation of datasets collected globally. Geographic areas covered by each of the datasets include: - EuroGEOSurveys Geochemical Baseline Database (FOREGS) - European continent - North American Soil Geochemical Landscapes (NASGLP) - continental United States and limited parts of Canada (see database key for more details) - National Geochemical Survey of Australia (GAA) - Australia - Alaska geochemical database (AGDB4) - Alaska - National Cooperative Soil Survey (NCSS) - Global measurements, but concentrated in the continental United States - Geochemical data for arable land and land under permanent grass cover in continental Europe (GEMAS) - continental Europe - Ecorespira-Amazon (ERA) - Geochemical data from the Amazon basin - Geochemical baseline data for New Zealand (NZGBS) - New Zealand - Geochemical data collected across continental Africa (AfSIS) - Measurements across Africa

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2021

Soil Water Potentials (1990-2018) from a calibrated ELM-FATES, and rooting depth analyses scripts, PA-BCI, Panama

This dataset has ELM-FATES simulation outputs and scripts that accompany and reproduce the manuscript, Chitra-Tarak et al. (2021) Hydraulically-vulnerable trees survive on deep-water access during droughts in a tropical forest. New Phytologist. https: //doi.org/10.1111/nph.17464”. This dataset is organized into four parts: (1) source code for the model, Energy Exascale Earth System Land Model-Functionally Assembled Advanced Terrestrial simulator (ELM-FATES) (2) R scripts associated with ELM-FATES simulations: to prepare files for parameter ensemble simulation runs, calibrate the model with observed hydrological fluxes observed at BCI, re-run the model with tuned parameters and obtain long-term soil matric potentials by depth (3) ELM-FATES hydrological simulation outputs at PA-BCI, Panama (1990-2018, daily), and (4) R scripts for inverse rooting depth modeling and post-processing: prepare datasets and inversely model tree rooting depths at BCI, test hypotheses about the relationship of rooting depths with plant hydraulic traits and drought-induced mortality, and reproduce the manuscript itself. n.b. 1. Data needed to run R scripts in (2) and (4) and thus to reproduce the manuscript are sourced from various public datasets. Refer to the data availability section of the accompanying manuscript. 2. This dataset version (prepared on May 12, 2021) matches the revised and published version of the accompanying manuscript (published May 16, 2021) replacing a previous version. In the publication, this dataset is cited as follows: Chitra-Tarak R, Xu C, Aguilar S, Anderson-Teixeira K, Chambers J, Detto M, Faybishenko B, Fisher R, Knox R, Koven C et al. 2020. Soil water potentials (1990–2018) from a calibrated ELM-FATES, and rooting depth analyses scripts, PA-BCI, Panama. 2.0. NGEE Tropics Data Collection. doi: 10.15486/ngt/1696806

54 ENVIRONMENTAL SCIENCES↗