Search NASA⌕ Search

SEARCH · Search NASA

Results for “index models”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Data and Scripts Associated with "Modeling Ecohydrological Responses of Vegetation to Urban Microclimates Using the E3SM Land Model"

This dataset supports the study of vegetation ecohydrological responses to urban microclimates using the land component of the Energy Exascale Earth System Model (ELM) at four urban sites in Knoxville, Tennessee, USA. It includes the model inputs, simulation outputs, and associated scripts for running ELM simulations and analyzing the resulting data. The Model_Inputs folder includes static surface data, satellite-derived phenology (i.e., leaf area index), and atmospheric forcing data used to drive ELM simulations. Detailed descriptions of these datasets are provided in Section 2.3.2 of the associated manuscript. The Model_Outputs folder contains simulation results for the baseline, treatment, and ensemble experiments. Outputs from the baseline and treatment simulations are provided as raw ELM NetCDF files. Because the raw outputs from the 4,000-member ensemble are prohibitively large, the ensemble results are provided as summarized CSV files, which also serve as the source data for Figure 5 of the associated manuscript. The Scripts folder contains three components: E3SM, the core codebase of the Energy Exascale Earth System Model (E3SM); elm-olmt, the Offline Land Model Testbed (OLMT) used to perform the simulations; and knoxville_elm, which contains the analysis scripts used to process model outputs and generate the figures and results presented in the associated manuscript. Additional information is provided in Scripts_readme.txt within the Scripts directory.

Lu, Xiaoman [ORNL] (ORCID:0000000306698780)↗
DOE OSTI2024

Observational benchmarks inform representation of soil organic carbon dynamics in land surface models

Abstract. Representing soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon–climate feedbacks. Machine learning models can help identify dominant environmental controllers and establish their functional relationships with SOC stocks. The resulting knowledge can be integrated into ESMs to reduce uncertainty and improve predictions of SOC dynamics over space and time. In this study, we used a large number of SOC field observations (n=54 000), geospatial datasets of environmental factors (n=46), and two machine learning approaches (namely random forest, RF, and generalized additive modeling, GAM) to (1) identify dominant environmental controllers of global and biome-specific SOC stocks, (2) derive functional relationships between environmental controllers and SOC stocks, and (3) compare the identified environmental controllers and predictive relationships with those in models used in Phase 6 of the Coupled Model Intercomparison Project (CMIP6). Our results showed that the diurnal temperature, drought index, cation exchange capacity, and precipitation were important observed environmental predictors of global SOC stocks. While the RF model identified 14 environmental factors that describe climatic, vegetation, and edaphic conditions as important predictors of global SOC stocks (R2=0.61, RMSE = 0.46 kg m−2), current ESMs oversimplify the relationships between environmental factors and SOC, with precipitation, temperature, and net primary productivity explaining > 96 % of the variability in ESM-modeled SOC stocks. Further, our study revealed notable disparities among the functional relationships between environmental factors and SOC stocks simulated by ESMs compared with observed relationships. To improve SOC representations in ESMs, it is imperative to incorporate additional environmental controls, such as the cation exchange capacity, and refine the functional relationships to align more closely with observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Revisiting the Anisotropic Complex Refractive Indices of Sodium Nitrate for Interpretation of the Reflectance Spectrum of Pressed Pellets

Reflectance spectroscopy is notoriously confounding in that the measured spectral response is highly dependent upon the morphology of the sample. Fortunately, all such perturbations are neatly encoded by the complex refractive index of the sample. Herein, we seek to quantitatively model the measured infrared reflectance spectrum of a specularly flat pressed pellet sample of the birefringent compound, sodium nitrate. Single crystals of sodium nitrate were synthesized via a slow evaporation process and spectroscopically analyzed using polarization-dependent infrared single-angle reflectance spectroscopy. The anisotropic complex refractive index was measured from 7500 to 300 cm-1 (1.33 to 33.33 µm). The deduced anisotropic optical constants were found to be consistent with those previously reported. Once measured and validated, the optical constants were applied to model the pressed pellet reflectance spectrum. It was evident that an average of the anisotropic refractive indices was insufficient to account for the measured pellet reflectance. In order to account for contributions of all possible microcrystalline orientations within the pellet, the Python package PYELLI was used to calculate a basis set of orientation-dependent reflection spectra spanning the distinct ? and ? Euler rotations of the uniaxial crystal. When the population of orientations was allowed to vary freely in a spectral fit analysis, the fit-deduced orientations were tightly clustered along f = 45º, hinting at residual anisotropy in the pressed pellet sample. Conversely, an equally valid spectral fit (with marginally worse fit metric) was obtained when the population was constrained to an isotropic distribution of orientations. Subsequent non-zero cross-polarization reflectance measurements likewise suggested anisotropy in the pellet. However, both grazing incidence wide-angle x-ray scattering and scanning electron microscopy measurements revealed that the microcrystal orientations at the surface of the pressed pellet sample were isotropically distributed (and that the average crystallite size was larger than ?/10). Application of the measured complex refractive indices for modeling the reflectance spectrum of the pressed pellet, and rectification of these seemingly contradictory observations will be discussed.

Wilhelm, Michael J.↗
DOE OSTI2023

Long-term survival and second malignant tumor prediction in pediatric, adolescent, and young adult cancer survivors using Random Survival Forests: a SEER analysis

Abstract Survival and second malignancy prediction models can aid clinical decision making. Most commonly, survival analysis studies are performed using traditional proportional hazards models, which require strong assumptions and can lead to biased estimates if violated. Therefore, this study aims to implement an alternative, machine learning (ML) model for survival analysis: Random Survival Forest (RSF). In this study, RSFs were built using the U.S. Surveillance Epidemiology and End Results to (1) predict 30-year survival in pediatric, adolescent, and young adult cancer survivors; and (2) predict risk and site of a second tumor within 30 years of the first tumor diagnosis in these age groups. The final RSF model for pediatric, adolescent, and young adult survival has an average Concordance index (C-index) of 92.9%, 94.2%, and 94.4% and average time-dependent area under the receiver operating characteristic curve (AUC) at 30-years since first diagnosis of 90.8%, 93.6%, 96.1% respectively. The final RSF model for pediatric, adolescent, and young adult second malignancy has an average C-index of 86.8%, 85.2%, and 88.6% and average time-dependent AUC at 30-years since first diagnosis of 76.5%, 88.1%, and 99.0% respectively. This study suggests the robustness and potential clinical value of ML models to alleviate physician burden by quickly identifying highest risk individuals.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Cloud drop number concentrations over the western North Atlantic Ocean: seasonal cycle, aerosol interrelationships, and other influential factors

Cloud drop number concentrations (N d ) over the western North Atlantic Ocean (WNAO) are generally highest during the winter (DJF) and lowest in summer (JJA), in contrast to aerosol proxy variables (aerosoloptical depth, aerosol index, surface aerosol mass concentrations, surface cloud condensation nuclei (CCN) concentrations) that generally peak inspring (MAM) and JJA with minima in DJF. Using aircraft, satellite remote sensing, ground-based in situ measurement data, and reanalysis data, we characterize factors explaining the divergent seasonal cycles and furthermore probe into factors influencing N d on seasonal timescales. The results can be summarized well by features most pronounced in DJF, including features associated with cold-air outbreak (CAO) conditions such as enhanced values of CAO index, planetary boundary layer height (PBLH),low-level liquid cloud fraction, and cloud-top height, in addition to winds aligned with continental outflow. Data sorted into high- and low-N d days in each season, especially in DJF, revealed that all of these conditions were enhanced on the high-N d days, including reduced sea level pressure and stronger wind speeds. Although aerosols may be more abundant in MAM and JJA, the conditions needed to activate those particles into cloud droplets are weaker than in colder months, which is demonstrated by calculations of the strongest (weakest) aerosol indirect effects in DJF (JJA) based on comparing N d to perturbations in four different aerosol proxy variables (total and sulfate aerosol optical depth, aerosol index, surface mass concentration of sulfate). We used three machine learning models and up to 14 input variables to infer about most influential factors related to N d for DJF and JJA, with the best performance obtained with gradient-boosted regression tree (GBRT) analysis. The model results indicated that cloud fraction was the most important input variable, followed by some combination (depending on season) of CAO index and surface mass concentrations of sulfate and organic carbon. Future work is recommended to further understand aspects uncovered here such as impacts of free tropospheric aerosol entrainment on clouds, degree of boundary layer coupling, wet scavenging, and giant CCN effects on aerosol–N d relationships, updraft velocity, and vertical structure of cloud properties such as adiabaticity that impact the satellite estimation of N d .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improved Constraints on Primordial Gravitational Waves using Planck , WMAP, and BICEP/ Keck Observations through the 2018 Observing Season

Here, we present results from an analysis of all data taken by the BICEP2, Keck Array, and BICEP3 CMB polarization experiments up to and including the 2018 observing season. We add additional Keck Array observations at 220 GHz and BICEP3 observations at 95 GHz to the previous 95/150/220 GHz dataset. The Q/U maps now reach depths of 2.8, 2.8, and 8.8 μK CMB arcmin at 95,150, and 220 GHz, respectively, over an effective area of ≈ 600 square degrees at 95 GHz and ≈ 400 square degrees at 150 and 220 GHz. The 220 GHz maps now achieve a signal-to-noise ratio on polarized dust emission exceeding that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz and evaluate the joint likelihood of the spectra versus a multicomponent model of lensed ΛCDM + r + dust + synchrotron + noise . The foreground model has seven parameters, and no longer requires a prior on the frequency spectral index of the dust emission taken from measurements on other regions of the sky. This model is an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint r 0.05 < 0.036 at 95% confidence. Running maximum likelihood search on simulations we obtain unbiased results and find that σ(r) = 0.009 . These are the strongest constraints to date on primordial gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Data for Yield from Iowa’s first commercial miscanthus fields: implications of spatial variability for productivity and sustainability beyond research plots

This dataset contains biomass yield measurements and associated vegetation index data collected from commercial Miscanthus × giganteus fields in eastern Iowa during the 2022–2023 growing seasons. The data support the analyses presented in the article: “Yield From Iowa's First Commercial Miscanthus Fields: Implications of Spatial Variability for Productivity and Sustainability Beyond Research Plots.” We collected 105 ground-truth biomass samples from four mature commercial fields (>4 years old) covering 92.81 ha. Samples were taken from 3 m² quadrats that were hand-harvested in alignment with commercial harvest timing. Stem biomass (excluding leaves) was weighed, moisture-corrected, and converted to dry-matter yield expressed in Mg DM ha⁻¹. Sampling locations were selected to capture spatial variability visible in aerial imagery and were recorded using RTK GPS. Each biomass observation was paired with vegetation indices derived from high-resolution PlanetScope satellite imagery (3 m resolution). Images were acquired throughout the growing season, and indices were calculated to evaluate their ability to predict end-of-season biomass yield. Statistical and machine learning approaches were used to identify key predictors, and a linear regression model based on end-of-July Green Normalized Difference Vegetation Index (GNDVI) was developed and evaluated. This repository includes the data used in that modeling workflow. Management practices, economic data, full imagery time series, and additional methodological details are described in the associated publication and are not included here. The dataset consists of three comma-separated value (CSV) files: 1. Combine_Groundtruth_Yield_VI_22_23.csv This file contains ground-truth biomass yield measurements and associated key vegetation index values collected during the 2022 and 2023 growing seasons. Rows: 105 observations Columns: Year — Year of observation (2022 or 2023) Field — Field location identifier Sample_number — Unique sample identifier GNDVI_End_Jul — Green Normalized Difference Vegetation Index calculated at end of July GNDVI_End_Aug — Green Normalized Difference Vegetation Index calculated at end of August NDRE_End_Aug — Normalized Difference Red Edge index calculated at end of August Biomass_Stem_Yield_MgDM/ha — Measured stem biomass yield (megagrams dry matter per hectare) 2. trainData_GNDVI.csv This file contains the subset of observations used to train the predictive relationship between July GNDVI and biomass yield. Rows: 76 observations Columns: Unnamed: 0 — Row index retained from the original data processing workflow GNDVI_End_Jul — GNDVI at end of July Stem_Yield_MgDM/ha — Observed stem biomass yield (Mg DM ha⁻¹) 3. testData_GNDVI.csv This file contains the test dataset used to evaluate model performance. Rows: 29 observations Columns: Unnamed: 0 — Row index retained from the original data processing workflow GNDVI_End_Jul — GNDVI at end of July Predicted_Yield_MgDM/ha — Model-predicted stem biomass yield (Mg DM ha⁻¹) Observed_Yield_MgDM/ha — Measured stem biomass yield (Mg DM ha⁻¹)

Potential yield, yield gap, in-field management, y↗
DOE OSTI2024

Unraveling the Nature of Lasing Emission from Hybrid Silicon Nitride and Colloidal Nanocrystal Photonic Crystals with Low Refractive Index Contrast

Silicon nitride is used for its low optical loss and high thermal stability, making it a suitable platform for visible-light applications in integrated photonic devices. However, its application has been limited due to inefficient light emission, a problem addressed by integrating various types of light emitters onto the platform. In particular, the integration of solution-processable colloidal nanocrystals (NCs) as optical gain materials onto the silicon nitride platform is a promising route but requires a more solid theoretical footing. By leveraging 2D surface-emitting photonic crystal structures combined with NCs, we effectively confine and manipulate light to achieve lasing from green to red. Building on this, we model the light–matter interactions of the low index contrast NC/nitride platform, validated by extensive experimental validations through Fourier imaging techniques, revealing the full photonic band structure and showing clear mode congestion. Finally, these comprehensive studies confirm the potential of hybrid NC-based structures for fully integrated on-chip laser applications and indicate routes for further improvement.

BIC Lasing↗
DOE OSTI2026

Resilience Measurement Framework For Post-deployment Artificial Intelligence (ai) Integrated Systems

Resilience is largely defined as the ability to adapt or recover from adverse conditions, stresses, attacks, or compromises on systems that use or are enabled by digital resources. In Artificial Intelligence Management and Research for Advanced Networked Testbed Hub (AMARANTH), resilience is measured in the amount of time it took from the beginning of a testing period for the model to reach predictions outside of the original 95% confidence interval or using the Kullback-Leibler (KL) divergence theorem, the Population Stability Index (PSI), and traditional methods such as root mean squared error (RMSE) threshold. Artificial Intelligence (AI) model drift is of significant concern when deploying AI-integrated systems into critical and/or secure environments. Drift can impact resilience of the AI-integrated system post-deployment and requires consistent maintenance and upkeep to ensure the model is accurate and precise. To quantify model drift and predict the point when a model's drift becomes unacceptable, we describe using Kullback-Leibler (KL) divergence, Population Stability Index (PSI) and/or confidence interval width estimations to determine the point of failure and time to failure of a model post-deployment. Through simple code functions, the KL-divergence, PSI, confidence interval, and root mean squared (RMSE) point of failures can be used to derive when a model needs to be maintained as well as the impact of adversarial action through statistical means.

Yockey, Patience [Idaho National Laboratory (INL),↗
DOE OSTI2022

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗
DOE OSTI2026

Urban Parameters Washington, DC 100m

132 Urban parameters based on building physical dimensions and location were generated for Washington, District of Columbia at 100m resolution using the NATURF model. To use the binary file with WRF, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way NUDAPT44 would be accessed.

Dumas, Melissa [ORNL] (ORCID:0000000233190846)↗
DOE OSTI2022

Machine learning prediction on the fractional free volume of polymer membranes

Fractional free volume (FFV) characterizes the microstructural level features of polymers and affects their properties including thermal, mechanical, and separation performance. Experimental measurements and theoretical analyses have been used to quantify the FFV of polymers, but challenges remain because of their limitations. Experimental measurements are laborious and based on semi empirical equations, while Bondi’s group contribution theory involves ambiguities like the determination of van der Waals volume and the choice of factor values in the theoretical equation. To efficiently evaluate the FFV of polymers, this study utilizes high-throughput molecular dynamics (MD) simulations to build a large dataset regarding polymer’s FFV. Based on this large dataset, we further build machine learning (ML) models to establish the composition-structure relation. Inspired by group contribution theory which correlates polymer’s functional groups to FFV, our ML models correlate polymer’s substructures or physico-chemical indexes to FFV. Here, our study first benchmarks the MD simulation protocol to obtain reliable FFV of polymers and then carries out high-throughput MD simulations for more than 6,500 homopolymers and 1,400 polyamides. Such a large and diverse dataset makes the well-trained ML models more generalizable, compared with the group contribution theory. The efficiency of a feed forward neural network model is further demonstrated by applying it to a hypothetical polyimide dataset of more than 8 million chemical structures. The predicted FFVs of hypothetical polyimides are further validated by MD simulations. The obtained FFVs of the 8 million polymers, plus their previously reported gas separation performances, demonstrate the promising capability of ML virtual screening for the discovery of polymer membranes with exceptional permeability/selectivity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pantropical Indo-Atlantic temperature gradient modulates multi-decadal AMOC variability in models and observations

Abstract Interconnections between ocean basins are recognized as an important driver of climate variability. Recent modeling evidence suggests that the North Atlantic climate can respond to persistent warming of the tropical Indian Ocean sea surface temperature (SST) relative to the rest of the tropics (rTIO). Here, we use observational data to demonstrate that multi-decadal changes in pantropical ocean temperature gradients lead to variations of an SST-based proxy of the Atlantic Meridional Overturning Circulation (AMOC). The largest contribution to this temperature gradient-AMOCconnection comes from gradients between the Indian and Atlantic Oceans. TherTIOindex yields the strongest connection of this tropical temperature gradient to theAMOC. Focusing on the internally generated signal in three observational products reveals that an SST-basedAMOCproxy index has closely followed low-frequency changes ofrTIOtemperature with about 26-year lag since 1870. Analyzing the pre-industrial control simulations of 44 CMIP6 climate models shows that theAMOCproxy index lags simulated mid-latitudeAMOCvariations by 4 ± 4 years. These model simulations reveal the mechanism connectingAMOCvariations to pantropical ocean temperature gradients at a 27 ± 2 years lag, matching the observed time lag in 28 out of the 44 analyzed models. rTIO temperature changes affect the North Atlantic climate through atmospheric planetary waves, impacting temperature and salinity in the subpolar North Atlantic, which modifies deep convection and ultimately the AMOC. Through this mechanism, observed internalrTIOvariations can serve as a multi-decadal precursor ofAMOCchanges with important implications forAMOCdynamics and predictability.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

Urban Parameters Los Angeles County 100m version 2

132 Urban parameters based on building physical dimensions and location were generated for the city of Los Angeles at 100m resolution using the NATURF model. To use the binary file with WRF, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way NUDAPT44 would be accessed. The kmz file can be visualized on Google Earth.

Sweet-Breu, Levi [Baylor University]↗
DOE OSTI2022

The 6-minute walk test predicts mortality in a pulmonary nontuberculous mycobacteria-predominant bronchiectasis cohort

Background: Bronchiectasis is a chronic lung condition frequently associated with nontuberculous mycobacteria pulmonary (NTM) disease. Persons with these conditions are at increased risk of mortality. Patient reported outcome (PRO) instruments and the 6-minute walk test (6MWT) have been shown to predict mortality for several lung conditions, but these measures have not been fully evaluated for bronchiectasis and NTM. Methods: We conducted a retrospective cohort study among adult patients enrolled in a natural history study of bronchiectasis at the National Heart, Lung, and Blood Institute. Electronic medical records were queried for demographic, clinical, microbiologic, radiographic, and PRO instrument data: St. George’s Respiratory Questionnaire (SGRQ), Medical Research Council Dyspnea Scale, and the Pulmonary Symptom Severity Score (PSSS). The study baseline date was defined as the patient’s first visit after January 1st, 2015 with a SGRQ or 6MWT completed. Follow-up was defined as the interval between the study baseline visit and date of death or December 31st, 2019. Sex-stratified Cox proportional-hazards regression was conducted to identify predictors of mortality. Separate models were run for each PRO and 6MWT measure, controlling for age, body mass index (BMI), fibrocavitary disease status, and M. abscessus infection. Results: In multivariable Cox proportional-hazards regression models, the PSSS-severity (aHR 1.29, 95% CI 1.04–1.59), the 6MWT total distance walked (aHR 0.938, 95% CI 0.896–0.981) and distance saturation product (aHR 0.930, 95% CI 0.887–0.974) independently predicted mortality. In addition, BMI was significantly predictive of mortality in all models. Conclusions: The 6MWT and a PRO instrument capturing symptom severity are independently predictive of mortality in our cohort of bronchiectasis patients.

6-minute walk test↗
DOE OSTI2021

A screening model for predicting injection well pressure buildup and plume extent in CO 2 geologic storage projects

We present a new screening model for predicting injection well pressure buildup and CO 2 plume migration for CO 2 geologic sequestration projects. The model requires only limited information and is quite accurate when compared to detailed simulation results. Such models can assist project developers during the early days of project planning (e.g., for 45Q related projects), and also help regulators perform simple checks against detailed numerical models. The screening model consists of two correlations: one for injectivity index in terms of the slope of the CO 2 fractional flow curve, and the second for total storage efficiency within the footprint of plume as a function of gravity number and slope of the CO 2 fractional flow curve. Using these two correlations, as well as a knowledge of some basic reservoir characteristics and estimates of fluid properties from standard correlations, the injection-well pressure buildup and CO 2 plume extent in the formation can be readily estimated. Finally, the new correlations show a good match with the results of the underlying simulations, and also provide good agreement with independent calculations for two example problems.

42 ENGINEERING↗
DOE OSTI2023

Multi-regime mixing modeling for local extinction and re-ignition in turbulent non-premixed flame by using LES/FDF method

Local extinction and re-ignition occur in turbulent non-premixed combustion when the Damköhler number is not large enough and the combustion is not fully mixing controlled. The occurrence of local extinction introduces locally extinguished flame holes followed by a premixed flame propagation toward the hole center to potentially reignite the flame. The co-existence of the non-premixed and premixed combustion regimes complicates the modeling since traditional combustion models are mostly for a single regime. In this work, we examine the effect of multi-regime mixing modeling in the transported filtered density function (FDF) method on the predictions of local extinction and reignition. Predictions of local extinction and reignition remain a challenge for the FDF method despite the progress made in the past. To account for the multi-regime combustion, two different mixing timescale models for non-premixed and premixed combustion are combined. A flame index based on the gradients of fuel and oxidizer is used to define a weighting factor to blend the two mixing timescale models. A turbulent jet non-premixed flame with substantial local extinction, the Sydney piloted jet flame L, is adopted as a test case to examine the performance of the multi-regime model in large-eddy simulation/FDF modeling. It is found that the traditional non-premixed mixing timescale model when combined with the modified Curl mixing leads to global extinction for the Sydney flame L without the presence of the premixed combustion regime. After accounting for the multi-regime combustion with proper detection of the different combustion regimes, the predictions for the flame statistics and the amount of local extinction are significantly improved. Furthermore, it suggests the need of including multi-regime combustion for the predictions of local extinction and re-ignition in turbulent non-premixed combustion configurations.

42 ENGINEERING↗
DOE OSTI2022

Development of an open-source regional data assimilation system in PEcAn v. 1.7.2: application to carbon cycle reanalysis across the contiguous US using SIPNET

Abstract. The ability to monitor, understand, and predict the dynamics of the terrestrial carbon cycle requires the capacity to robustly and coherently synthesize multiple streams of information that each provide partial information about different pools and fluxes. In this study, we introduce a new terrestrial carbon cycle data assimilation system, built on the PEcAn model–data eco-informatics system, and its application for the development of a proof-of-concept carbon “reanalysis” product that harmonizes carbon pools (leaf, wood, soil) and fluxes (GPP, Ra, Rh, NEE) across the contiguous United States from 1986–2019. We first calibrated this system against plant trait and flux tower net ecosystem exchange (NEE) using a novel emulated hierarchical Bayesian approach. Next, we extended the Tobit–Wishart ensemble filter (TWEnF) state data assimilation (SDA) framework, a generalization of the common ensemble Kalman filter which accounts for censored data and provides a fully Bayesian estimate of model process error, to a regional-scale system with a calibrated localization. Combined with additional workflows for propagating parameter, initial condition, and driver uncertainty, this represents the most complete and robust uncertainty accounting available for terrestrial carbon models. Our initial reanalysis was run on an irregular grid of ∼ 500 points selected using a stratified sampling method to efficiently capture environmental heterogeneity. Remotely sensed observations of aboveground biomass (Landsat LandTrendr) and leaf area index (LAI) (MODIS MOD15) were sequentially assimilated into the SIPNET model. Reanalysis soil carbon, which was indirectly constrained based on modeled covariances, showed general agreement with SoilGrids, an independent soil carbon data product. Reanalysis NEE, which was constrained based on posterior ensemble weights, also showed good agreement with eddy flux tower NEE and reduced root mean square error (RMSE) compared to the calibrated forecast. Ultimately, PEcAn's new open-source regional data assimilation framework provides a scalable workflow for harmonizing multiple data constraints and providing a uniform synthetic platform for carbon monitoring, reporting, and verification (MRV) as well as accelerating terrestrial carbon cycle research.

54 ENVIRONMENTAL SCIENCES↗