SEARCH · Search NASA

Results for “Association study”

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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Computational schemes for the Magnus expansion of the in-medium similarity renormalization group

The in-medium similarity renormalization group (IMSRG) is a popular many-body method used for computations of nuclei. It solves the many-body Schrödinger equation through a continuous unitary transformation of the many-body Hamiltonian. The IMSRG transformation is typically truncated at the normal-ordered two-body level, the IMSRG(2), but recently several approaches have been developed to capture the effects of normal-ordered three-body operators, the IMSRG(3). In particular, a factorized approximation to the IMSRG(3) proposes to capture the leading effects of three-body operators at the same computational cost as the IMSRG(2) approximation. This approach often employs an approximate scheme for solving the IMSRG equations, the so-called hunter-gatherer scheme. In this work, I study the uncertainty associated with this scheme. I find that the hunter-gatherer scheme differs by up to 7MeV for ground-state energies and 0.5MeV for excitation energies from standard IMSRG(2) approaches. These differences are in some cases comparable to the expected size of IMSRG(3) corrections.

39 ≤ A ≤ 58↗

DOE OSTI2021

Atacama Cosmology Telescope: Combined kinematic and thermal Sunyaev-Zel’dovich measurements from BOSS CMASS and LOWZ halos

The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel’dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope DR5 and Planck in combination with the CMASS (mean redshift $\langle z \rangle$ = 0.55 and host halo mass $\langle M_{vir} \rangle$ = 3 × 10 13 M ⊙ ) and LOWZ ( $\langle z \rangle$ = 0.31 , $\langle M_{vir} \rangle$ = 5 × 10 13 M ⊙ ) galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothesis at 6.5σ, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally > 90 σ for the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e., the electron thermal pressure profile of the same CMASS objects, and reject the no-tSZ hypothesis at 10σ. We combine tSZ and kSZ measurements to estimate the electron temperature to 20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ) signal at 2.9 (13.9) σ , and leave their interpretation to future work. This paper and the companion paper demonstrate that current CMB experiments can detect and resolve gas profiles in low mass halos and at high redshifts, which are the most sensitive to feedback in galaxy formation and the most difficult to measure any other way. They will be a crucial input to cosmological hydrodynamical simulations, thus improving our understanding of galaxy formation. These precise gas profiles are already sufficient to reduce the main limiting theoretical systematic in galaxy-galaxy lensing: baryonic uncertainties. Future such measurements will thus unleash the statistical power of weak lensing from the Rubin, Euclid and Roman observatories. Finally, our stacking software ThumbStackis publicly available and directly applicable to future Simons Observatory and CMB-S4 data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Analytical and numerical studies of dark current in radiofrequency structures for short-pulse high-gradient acceleration

High-gradient acceleration is a key research area that could enable compact linear accelerators for future colliders, light sources, and other applications. In the pursuit of high-gradient operation, rf breakdown limits the attainable accelerating gradient in normal-conducting rf structures. Recent experiments at the Argonne Wakefield Accelerator suggest a promising approach: using short rf pulses with durations of a few nanoseconds. Experimental studies show that these 𝒪⁡(1⁢ ⁢ns) rf pulses can mitigate breakdown limitations, resulting in higher gradients. For example, an electric field of nearly 400 MV/m was achieved in an 𝑋-band photoemission gun driven by 6-ns-long rf pulses, with rapid rf conditioning and low dark current observed. Despite these promising results, the short-pulse regime remains an underexplored parameter space, and rf breakdown physics under nanosecond-long pulses requires further investigation. In this paper, we present analytical and numerical simulations of dark current dynamics in accelerating cavities operating in the short-pulse regime. We study breakdown-associated processes spanning different time scales, including field emission, multipacting, and plasma formation, using simulations of the 𝑋-band photogun cavities. The results reveal the advantages of using short rf pulses to reduce dark current and mitigate rf breakdown, offering a path toward a new class of compact accelerators with enhanced performance and reduced susceptibility to breakdown.

Linear accelerators↗

DOE OSTI2022

Rapid changes in coastal ocean microbiomes uncoupled with shifts in environmental variables

Disturbances, here defined as events that directly alter microbial community composition, are commonly studied in host–associated and engineered systems. Notably, in spite of global change both altering environmental averages and increasing extreme events, there has been relatively little research into the causes, persistence and population–level impacts of disturbance in the dynamic coastal ocean. Here, we utilize 3 years of observations from a coastal time series to identify disturbances based on the largest week–over–week changes in the microbiome (i.e. identifying disturbance as events that alter the community composition). In general, these microbiome disturbances were not clearly linked to specific environmental factors and responsive taxa largely differed, aside from SAR11, which generally declined. However, several disturbance metagenomes identified increased phage–associated genes, suggesting that unexplained community shifts might be caused by increased mortality. Furthermore, a category 1 hurricane, the only event that would likely be classified a priori as an environmental disturbance, was not an outlier in microbiome composition, but did enhance a bloom in seasonally abundant phytoplankton. Thus, as extreme environmental changes intensify, assumptions of what constitutes a disturbance should be re–examined in the context of ecological history and microbiome responses.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

A new mechanism for void-cascade interaction from nondestructive depth-resolved atomic-scale measurements of ion irradiation–induced defects in Fe

The nondestructive investigation of single vacancies and vacancy clusters in ion-irradiated samples requires a depth-resolved probe with atomic sensitivity to defects. The recent development of short-pulsed positron beams provides such a probe. Here, we combine depth-resolved Doppler broadening and positron annihilation lifetime spectroscopies to identify vacancy clusters in ion-irradiated Fe and measure their density as a function of depth. Despite large concentrations of dislocations and voids in the pristine samples, positron annihilation measurements uncovered the structure of vacancy clusters and the change in their size and density with irradiation dose. When combined with transmission electron microscopy measurements, the study demonstrates an association between the increase in the density of small vacancy clusters with irradiation and a remarkable reduction in the size of large voids. This, previously unknown, mechanism for the interaction of cascade damage with voids in ion-irradiated materials is a consequence of the high porosity of the initial microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Vasomotion in Human Fingers

Introduction: We describe methods by which vasomotion can be recorded in awake and anesthetized human subjects without significant interference from other spontaneous vascular oscillations. Methods: In three separate studies, we used photoplethysmography (PPG) to record vasomotion in fingertips. In Study 1, we induced chemical sympathectomy in the studied hand of 11 awake subjects who received intravenous dexmedetomidine infusions. In Study 2, we administered four progressively increasing intravenous dexmedetomidine infusions to 16 awake volunteers. In Study 3, we recorded vasomotion simultaneously from 6 fingers of 7 patients who were under dexmedetomidine-based anesthesia. Five-minute epochs of PPG recordings that displayed slow vascular oscillations were analyzed for frequency and amplitude. Results: In Study 1, vasomotion frequencies were 0.025 ± 0.008 Hz. In Study 2, vasomotion frequencies were 0.033 ± 0.006 Hz, and 0.032 ± 0.008 Hz during the two highest dexmedetomidine infusion steps. In Study 3, vasomotion frequencies ranged from 0.020 to 0.037 Hz and were observed in all 6 fingers, with no synchrony between the six fingers. Conclusion: The vascular oscillations we observed without significant interference from other spontaneous oscillations are independent of neural activity (Study 1), local in nature (Study 3), and associated with alpha-2-adrenoceptor activation, consistent with known properties of vasomotion.

Cardiovascular System & Cardiology↗

DOE OSTI2020

Psychological Distress and Alzheimer’s Disease Mortality in the United States: Results from the 1997–2014 National Health Interview Survey-National Death Index Record Linkage Study

Objective: This study examines the association between psychological distress and Alzheimer’s disease mortality among US adults aged ≥45. Methods: We analyzed the Kessler 6-item psychological distress scale as a risk factor for Alzheimer’s mortality using the pooled 1997–2014 National Health Interview Survey (NHIS)- National Death Index (NDI) database ( N = 265,089). Cox regression was used to model mortality as a function of psychological distress and sociodemographic and behavioral covariates. Results: The Alzheimer’s mortality risk was 97% higher (HR = 1.97; 95% confidence interval [CI] = 1.37, 2.84) in adults with serious psychological distress compared with those without psychological distress, controlling for sociodemographic covariates. The relative mortality risk remained statistically significant (HR = 1.49; 95% CI = 1.04, 2.13) after additional adjustment for smoking, alcohol consumption, health status, activity limitation, and body mass index. Discussion: US adults had significantly higher risks of Alzheimer’s disease mortality at higher psychological distress levels. These findings underscore the significance of addressing psychological well-being as a strategy for reducing Alzheimer’s disease mortality.

Singh, Gopal K.↗

DOE OSTI2024

Differential patterns of reproductive and lifestyle risk factors for breast cancer according to birth cohorts among women in China, Japan and Korea

Background: The birth cohort effect has been suggested to influence the rate of breast cancer incidence and the trends of associated reproductive and lifestyle factors. We conducted a cohort study to determine whether a differential pattern of associations exists between certain factors and breast cancer risk based on birth cohorts. Methods: This was a cohort study using pooled data from 12 cohort studies. We analysed associations between reproductive (menarche age, menopause age, parity and age at first delivery) and lifestyle (smoking and alcohol consumption) factors and breast cancer risk. We obtained hazard ratios (HRs) with 95% confidence intervals (CIs) using the Cox proportional hazard regression analysis on the 1920s, 1930s, 1940s and 1950s birth cohorts. Results: Parity was found to lower the risk of breast cancer in the older but not in the younger birth cohort, whereas lifestyle factors showed associations with breast cancer risk only among the participants born in the 1950s. In the younger birth cohort group, the effect size was lower for parous women compared to the other cohort groups (HR [95% CI] 0.86 [0.66–1.13] compared to 0.60 [0.49–0.73], 0.46 [0.38–0.56] and 0.62 [0.51–0.77]). Meanwhile, a higher effect size was found for smoking (1.45 [1.14–1.84] compared to 1.25 [0.99–1.58], 1.06 [0.85–1.32] and 0.86 [0.69–1.08]) and alcohol consumption (1.22 [1.01–1.48] compared to 1.10 [0.90–1.33], 1.15 [0.96–1.38], and 1.07 [0.91–1.26]). Conclusion: We observed different associations of parity, smoking and alcohol consumption with breast cancer risk across various birth cohorts.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Urbanization and malaria have a contextual relationship in endemic areas: A temporal and spatial study in Ghana

In West Africa, malaria is one of the leading causes of disease-induced deaths. Existing studies indicate that as urbanization increases, there is corresponding decrease in malaria prevalence. However, in malaria-endemic areas, the prevalence in some rural areas is sometimes lower than in some peri-urban and urban areas. Therefore, the relationship between the degree of urbanization, the impact of living in urban areas, and the prevalence of malaria remains unclear. This study explores this association in Ghana, using epidemiological data at the district level (2015–2018) and data on health, hygiene, and education. We applied a multilevel model and time series decomposition to understand the epidemiological pattern of malaria in Ghana. Then we classified the districts of Ghana into rural, peri-urban, and urban areas using administratively defined urbanization, total built areas, and built intensity. We converted the prevalence time series into cross-sectional data for each district by extracting features from the data. To predict the determinant most impacting according to the degree of urbanization, we used a cluster-specific random forest. We find that prevalence is impacted by seasonality, but the trend of the seasonal signature is not noticeable in urban and peri-urban areas. While urban districts have a slightly lower prevalence, there are still pockets with higher rates within these regions. These areas of high prevalence are linked to proximity to water bodies and waterways, but the rise in these same variables is not associated with the increase of prevalence in peri-urban areas. The increase in nightlight reflectance in rural areas is associated with an increased prevalence. We conclude that urbanization is not the main factor driving the decline in malaria. However, the data indicate that understanding and managing malaria prevalence in urbanization will necessitate a focus on these contextual factors. Finally, we design an interactive tool, ’malDecision’ that allows data-supported decision-making.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost: Model Inputs and Outputs

This dataset contains model inputs and simulation outputs associated with the study “Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost.” The ecosys model was extended to represent thermokarst processes in laterally coupled two-dimensional transects and applied to the Scotty Creek basin (Northwest Territories, Canada). The ecosys model and documentation is available for download at https://github.com/jinyun1tang/ECOSYS. The dataset contains the inputs used to run the model for this study (meteorological forcing data (1950–2018) and site-specific soil and vegetation parameters) and model outputs used to generate Figure 3 and Figure 4 in the paper (simulated ground thermal and hydrologic states, vegetation dynamics, and carbon fluxes).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Historic climate, cosmogenic 10Be, denudation-rate, and geospatial datasets from the Pikes Peak region, Colorado, USA

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of elevation-dependent denudation rates on Pikes Peak in the Front Range of the Rocky Mountains, Colorado, USA. The package includes GIS layers used to produce the study-area map, including sample locations, sample watershed boundaries, the Pikes Peak batholith, Pleistocene glacier extent, weather station locations, and elevation and hillshade rasters, together with comma-separated value (CSV) tables and matching CSV data dictionaries. These mapped layers provide the geographic framework for interpreting denudation patterns across the Pikes Peak region and for relating sample locations to watershed geometry, bedrock setting, glacial history, and nearby climate stations. The first group of tables reports climate and geospatial context for the study area. These files include station-based temperature and precipitation data used to characterize elevational gradients in mean annual climate and monthly climate seasonality, sample locations, denudation-rate and topographic metrics, fixed frost-cracking model parameters, frost-cracking intensity and precipitation-frequency metrics, and stream-power inversion results. Together, these data provide the basis for evaluating how denudation varies with elevation, climate, and landscape form across sampled catchments on Pikes Peak. The second group of tables reports cosmogenic nuclide and erosion-model results used in the denudation analysis. Included files contain accelerator mass spectrometry (AMS) measurements for in situ-produced cosmogenic beryllium-10 (10Be), including sample identifiers, measured 10Be:9Be ratios, analytical uncertainties, carrier mass, quartz mass, blank corrections, blank-group statistics, and calculated 10Be concentrations and uncertainties. Additional tables summarize stream-power-law inversion results for sampled catchments, including optimized model parameters, predicted erosion rates, residual metrics, channel-pixel counts, and convergence status, as well as regression equations and summary statistics used to evaluate relationships among elevation, climate, frost cracking, precipitation forcing, and denudation rate. The package contains GIS files, comma-separated value files (.csv), Microsoft Excel files (.xlsx), CSV data dictionaries, a file-level metadata table, and a readme text file.

10Be cosmogenic nuclides↗

DOE OSTI2023

The Effect of Prostate-Specific Antigen (PSA) Test on Radiation Risk Estimate for Prostate Cancer Incidence among Atomic-Bomb Survivors

Following our previous report on the radiation doseresponse for prostate cancer incidence rates in the Life Span Study (LSS) cohort of atomic bomb survivors, we reevaluated the radiation-related risk adjusting for differences in baseline cancer incidence rates among three subsets of the LSS cohort defined by the timing of their first participation in biennial health examinations offered to the Adult Health Study (AHS) sub-cohort members and prostate-specificantigen (PSA) testing status for AHS participants: 1. nonAHS participants, 2. AHS participants before receiving PSA test, and 3. AHS participants after receiving PSA test. We found a 2.9-fold increase in the baseline incidence rates among AHS participants after receiving PSA test. After adjusting for the PSA-testing-status effects on the baseline rates the estimated excess relative risk (ERR) per Gy was 0.54 (95% CI: 0.15, 1.05), which was almost identical to the previously reported unadjusted ERR estimate (0.57, 95% CI: 0.21, 1.00). The current results confirmed that, while the PSA testing among AHS participants increased the baseline incidence rates, it did not impact the radiation risk estimate, strengthening the previously reported dose-response relationship for prostate cancer incidence in the LSS. As the use of PSA tests continue in screening and medical settings, analyses of possible effects of PSA testing should be an important aspect of future epidemiological studies of the association between radiation exposure and prostate cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2020

Hydropower Geotechnical Foundations: Executive Summary on Current Practice and Innovation Opportunities for Low-Head Applications

The following content provides a companion executive summary of Hydropower Geotechnical Foundations: Current Practice and Innovation Opportunities for Low-Head Applications (DeNeale et al. 2020), a technical report developed by Oak Ridge National Laboratory (ORNL) and Knight Piésold Consulting, with funding from the US Department of Energy (DOE) Water Power Technologies Office (WPTO). This executive summary provides a brief synopsis of the broader information contained in the technical report. This study (and the associated technical report): (1) covers background information on hydropower; (2) describes various characteristics relevant to hydropower foundations for undeveloped US streams; (3) presents the current state of practice in foundation development; (4) provides a representative assessment of conventional hydropower foundation costs and timelines; (5) addresses key challenges facing conventional hydropower foundations; and (6) presents opportunity areas and examples for innovative hydropower foundation technologies while highlighting some advances in non-hydropower industries.

13 HYDRO ENERGY↗

DOE OSTI2024

PyPop: a mature open-source software pipeline for population genomics

Python for Population Genomics (PyPop) is a software package that processes genotype and allele data and performs large-scale population genetic analyses on highly polymorphic multi-locus genotype data. In particular, PyPop tests data conformity to Hardy-Weinberg equilibrium expectations, performs Ewens-Watterson tests for selection, estimates haplotype frequencies, measures linkage disequilibrium, and tests significance. Standardized means of performing these tests is key for contemporary studies of evolutionary biology and population genetics, and these tests are central to genetic studies of disease association as well. Here, we present PyPop 1.0.0, a new major release of the package, which implements new features using the more robust infrastructure of GitHub, and is distributed via the industry-standard Python Package Index. New features include implementation of the asymmetric linkage disequilibrium measures and, of particular interest to the immunogenetics research communities, support for modern nomenclature, including colon-delimited allele names, and improvements to meta-analysis features for aggregating outputs for multiple populations.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Nitrogen fixation and mucilage production on maize aerial roots is controlled by aerial root development and border cell functions

Exploring natural diversity for biological nitrogen fixation in maize and its progenitors is a promising approach to reducing our dependence on synthetic fertilizer and enhancing the sustainability of our cropping systems. We have shown previously that maize accessions from the Sierra Mixe can support a nitrogen-fixing community in the mucilage produced by their abundant aerial roots and obtain a significant fraction of their nitrogen from the air through these associations. In this study, we demonstrate that mucilage production depends on root cap and border cells sensing water, as observed in underground roots. The diameter of aerial roots correlates with the volume of mucilage produced and the nitrogenase activity supported by each root. Young aerial roots produce more mucilage than older ones, probably due to their root cap’s integrity and their ability to produce border cells. Transcriptome analysis on aerial roots at two different growth stages before and after mucilage production confirmed the expression of genes involved in polysaccharide synthesis and degradation. Genes related to nitrogen uptake and assimilation were up-regulated upon water exposure. Altogether, our findings suggest that in addition to the number of nodes with aerial roots reported previously, the diameter of aerial roots and abundance of border cells, polysaccharide synthesis and degradation, and nitrogen uptake are critical factors to ensure efficient nitrogen fixation in maize aerial roots.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Prompt Searches for Very-high-energy γ -Ray Counterparts to IceCube Astrophysical Neutrino Alerts

The search for sources of high-energy astrophysical neutrinos can be significantly advanced through a multimessenger approach, which seeks to detect the γ-rays that accompany neutrinos as they are produced at their sources. Multimessenger observations have so far provided the first evidence for a neutrino source, illustrated by the joint detection of the flaring blazar TXS 0506+056 in high-energy (E > 1 GeV) and very-high-energy (VHE; E > 100 GeV) γ-rays in coincidence with the high-energy neutrino IceCube-170922A, identified by IceCube. Imaging atmospheric Cherenkov telescopes (IACTs), namely FACT, H.E.S.S., MAGIC, and VERITAS, continue to conduct extensive neutrino target-of-opportunity follow-up programs. These programs have two components: follow-up observations of single astrophysical neutrino candidate events (such as IceCube-170922A), and observation of known γ-ray sources after the identification of a cluster of neutrino events by IceCube. Here we present a comprehensive analysis of follow-up observations of high-energy neutrino events observed by the four IACTs between 2017 September (after the IceCube-170922A event) and 2021 January. Our study found no associations between γ-ray sources and the observed neutrino events. We provide a detailed overview of each neutrino event and its potential counterparts. Furthermore, a joint analysis of all IACT data is included, yielding combined upper limits on the VHE γ-ray flux.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Analysis of genomic signatures associated with Variovorax endosphere colonization

This repository contains the analysis code and supporting datasets associated with the study “Genomic signatures in Variovorax enabling colonization of the Populus endosphere.” Beals DG, Carper DL, Hochanadel LH, Jawdy SS, Klingeman DM, Piatkowski BT, Weston DJ, Doktycz MJ, Pelletier DA. 2026. Genomic signatures in Variovorax enabling colonization of the Populus endosphere. mSystems 11:e01605-25. https://doi.org/10.1128/msystems.01605-25 The scripts are organized sequentially (01–07) and document the workflows used for: Sequence-read alignment and feature counting Orthogroup and KEGG Ortholog annotation Count normalization Statistical analysis and aggregation Generation of manuscript figures and tables Repository contents The uncompressed files are the finalized, formatted datasets used to generate the figures and tables reported in the study, including the supplemental CSV files referenced in the manuscript. The accompanying ZIP archive contains the complete codebase and example data_input/ and data_output/ directories illustrating the organization and execution of the analytical workflow. Individual scripts identify the corresponding manuscript analyses and figure panels. Raw sequencing data Raw sequencing reads are available through the NCBI Sequence Read Archive under BioProject accession PRJNA1322484.

Beals, Delaney [ORNL] (ORCID:0000000306274574)↗

DOE OSTI2026

Data and Code for: Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits

This repository contains the simulation outputs and processing scripts associated with the study of winter wheat traits across the United States, utilizing the Ecosys agroecosystem model. The dataset includes model results for both rainfed and irrigated winter wheat systems, supporting the findings presented in the manuscript titled "Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits." Data includes the original Ecosys simulation outputs (archived in .db format within the compressed .zip files) and extracted analysis data (stored in .pkl files for efficient processing). Python code for data processing and figure generation is provided in a Jupyter notebook. External Observational Datasets should refer to the following official repositories for the input and validation data used in this study. The eddy covariance data from the AmeriFlux network (https://ameriflux.lbl.gov/). Climate-forcing data of NLDAS-2 from NASA LDAS (https://ldas.gsfc.nasa.gov/nldas/nldas-2-forcing-data). Soil data from the Gridded Soil Survey Geographic Database (gSSURGO), available at (https://www.nrcs.usda.gov/resources/data-and-reports/gridded-soil-survey-geographic-gssurgo-database). Crop yields, planting and harvest dates from the USDA public databases (https://quickstats.nass.usda.gov/; https://webapp.rma.usda.gov/apps/actuarialinformationbrowser/CropCriteria.aspx). Satellite-derived SLOPE GPP data from ORNL DAAC (https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1786). Land use and crop progress information from the USDA Crop Data Layer and Crop Progress and Condition Gridded Layers (https://www.nass.usda.gov/Research_and_Science/). The Ecosys model code is available online at https://github.com/jinyun1tang/ECOSYS.

Wheat↗