SEARCH · Search NASA

Results for “Forest Modeling”

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 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗

DOE OSTI2023

Towards fast, accurate predictions of RF simulations via data-driven modeling: Forward and lateral models

Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modeling and real-time control applications. More accurate simulations with fast electron diffusion are even slower, requiring multiple hours of run time with parallel processing. The machine learning models use a database of 16,000+ GEN-RAY/CQL3D simulations for training, validation, and testing. Latin hypercube sampling methods implemented in πScope ensure that the database covers the range of 9 input parameters (n e0 , T e0 , I p , B t , R 0 , n ∥︀ , Z e f f , V loop , P LHCD ) with sufficient density in all regions of parameter space. The surrogate models reduce the computation time from minutes-hours to ms with high accuracy across the input parameter space. Data-driven surrogate models also allow for solving inverse and “lateral” problems. A surrogate model for the inverse problem maps from a desired current drive or power deposition profile to a set of input parameters that would result in such a profile, while a surrogate model for the lateral problem maps from a measured experimental quantity such as hard x-ray emission to a current drive or power deposition profile. In conclusion, the πScope database creation workflow is flexible and applicable to other RF simulation codes such as TORIC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Towards Fast and Accurate Predictions of Radio Frequency Power Deposition and Current Profile via Data-driven Modeling

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Assimilation of repeated woody biomass observations constrains decadal ecosystem carbon cycle uncertainty in aggrading forests

Forest carbon sink strengths are governed by plant growth, mineralization of dead organic matter, and disturbance. Across landscapes, remote sensing can provide information about aboveground states of forests and this information can be linked to models to estimate carbon cycling in forests close to steady state. For aggrading forests this approach is more challenging and has not been demonstrated. Here in this paper we apply a Bayesian approach, linking a simple model to a range of data, to evaluate their information content, for two aggrading forests. We compare high information content analyses using local observations with retrievals using progressively sparser remotely sensed information (repeated, single, and no woody biomass observations). The net biome productivity of both forests is constrained to be a net sink with <2 Mg C ha -1 yr -1 variation across the range of inputs. However, the sequestration of particular carbon pool(s) varies with assimilated biomass information. Assimilation of repeated biomass observations reduces uncertainty and/or bias in all ecosystem C pools not just wood, compared to analyses using single or no stock information. As verification, our repeated biomass analysis explains 78–86% of variation in litter dynamics at one forest, while at the second forest total dead organic matter estimates are within observational uncertainty. The uncertainty of retrieved ecosystem traits in the repeated biomass analysis is reduced by up to 50% compared to analyses with less biomass information. This study quantifies the importance of repeated woody observations in constraining the dynamics of both wood and dead organic matter, highlighting the benefit of proposed remote sensing missions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Near‐future forest vulnerability to drought and fire varies across the western United States

Abstract Recent prolonged droughts and catastrophic wildfires in the western United States have raised concerns about the potential for forest mortality to impact forest structure, forest ecosystem services, and the economic vitality of communities in the coming decades. We used the Community Land Model (CLM) to determine forest vulnerability to mortality from drought and fire by the year 2049. We modified CLM to represent 13 major forest types in the western United States and ran simulations at a 4‐km grid resolution, driven with climate projections from two general circulation models under one emissions scenario (RCP 8.5). We developed metrics of vulnerability to short‐term extreme and prolonged drought based on annual allocation to stem growth and net primary productivity. We calculated fire vulnerability based on changes in simulated future area burned relative to historical area burned. Simulated historical drought vulnerability was medium to high in areas with observations of recent drought‐related mortality. Comparisons of observed and simulated historical area burned indicate simulated future fire vulnerability could be underestimated by 3% in the Sierra Nevada and overestimated by 3% in the Rocky Mountains. Projections show that water‐limited forests in the Rocky Mountains, Southwest, and Great Basin regions will be the most vulnerable to future drought‐related mortality, and vulnerability to future fire will be highest in the Sierra Nevada and portions of the Rocky Mountains. High carbon‐density forests in the Pacific coast and western Cascades regions are projected to be the least vulnerable to either drought or fire. Importantly, differences in climate projections lead to only 1% of the domain with conflicting low and high vulnerability to fire and no area with conflicting drought vulnerability. Our drought vulnerability metrics could be incorporated as probabilistic mortality rates in earth system models, enabling more robust estimates of the feedbacks between the land and atmosphere over the 21st century.

Buotte, Polly C.↗

DOE OSTI2022

The Latent Dirichlet Allocation model applied to airborne LiDAR data: A case study on mapping forest degradation associated with fragmentation and fire in the Amazon region

1. LiDAR data are being increasingly used to provide a detailed characterization of the vertical profile of forests. This characterization enables the generation of new insights on the influence of environmental drivers and anthropogenic disturbances on forest structure as well as on how forest structure influences important ecosystem functions and services. Unfortunately, extracting information from LiDAR data in a way that enables the spatial visualization of forest structure, as well as its temporal changes, is challenging due to the high dimensionality of these data. 2. In this study, we show how the Latent Dirichlet Allocation model applied to LiDAR data (LidarLDA) can be used to identify forest structural types and how the relative abundance of these forest types changes throughout the landscape. The code to fit this model is made available through the open-source r package LidarLDA in github. We illustrate the use of LidarLDA both with simulated data and data from a large-scale fire experiment in the Brazilian Amazon region. 3. Using simulated data, we demonstrate that LidarLDA accurately identifies the number of forest types as well as their spatial distribution and absorptance probabilities. For the empirical data, we found that LidarLDA detects both landscape-level patterns in forest structure as well as the strong interacting effect of fire and forest fragmentation on forest structure based on the experimental fire plots. More specifically, LidarLDA reveals that proximity to forest edge exacerbates the impact of fires, and that burned forests remain structurally different from unburned areas for at least 7 years, even when burned only once. Importantly, LidarLDA generates insights on the 3D structure of forest that cannot be obtained using more standard approaches that just focus on top-of-the-canopy information (e.g. canopy height models based on LiDAR data). 4. By enabling the mapping of forest structure and its temporal changes, we believe that LidarLDA will be of broad utility to the ecological research community.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A machine learning approach targeting parameter estimation for plant functional type coexistence modeling using ELM-FATES (v2.0)

Tropical forest dynamics play a crucial role in the global carbon, water, and energy cycles. However, realistically simulating the dynamics of competition and coexistence between different plant functional types (PFTs) in tropical forests remains a significant challenge. This study aims to improve the modeling of PFT coexistence in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a vegetation demography model implemented in the Energy Exascale Earth System Model (E3SM) land model (ELM), ELM-FATES. Specifically, we explore (1) whether plant trait relationships established from field measurements can constrain ELM-FATES simulations and (2) whether machine learning (ML)-based surrogate models can emulate the complex ELM-FATES model and optimize parameter selections to improve PFT coexistence modeling. We conducted three ensembles of ELM-FATES experiments at a tropical forest site near Manaus, Brazil. By comparing the ensemble experiments without (Exp-CTR) and with (Exp-OBS) consideration of observed trait relationships, we found that accounting for these relationships slightly improves the simulations of water, energy, and carbon variables when compared to observations but degrades the simulation of PFT coexistence. Using ML-based surrogate models trained on Exp-CTR, we optimized the trait parameters in ELM-FATES and conducted another ensemble of experiments (Exp-ML) with these optimized parameters. The proportion of PFT coexistence experiments significantly increased from 21 % in Exp-CTR to 73 % in Exp-ML. After filtering the experiments that allow for PFT coexistence to agree with observations (within 15 % tolerance), 33 % of the Exp-ML experiments were retained, which is a significant improvement compared to the 1.4 % in Exp-CTR. Exp-ML also accurately reproduces the annual means and seasonal variations in water, energy, and carbon fluxes and the field inventory of aboveground biomass. This study represents a reproducible method that utilizes machine learning to identify parameter values that improve model fidelity against observations and PFT coexistence in vegetation demography models for diverse ecosystems. Our study also suggests the need for new mechanisms to enhance the robust simulation of coexisting plants in ELM-FATES and has significant implications for modeling the response and feedbacks of ecosystem dynamics to climate change.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1990

Multiband radar characterization of forest biomes

The utility of airborne and orbital SAR in classification, assessment, and monitoring of forest biomes is investigated through analysis of orbital synthetic aperature radar (SAR) and multifrequency and multipolarized airborne SAR imagery relying on image tone and texture. Preliminary airborne SAR experiments and truck-mounted scatterometer observations demonstrated that the three dimensional structural complexity of a forest, and the various scales of temporal dynamics in the microwave dielectric properties of both trees and the underlying substrate would severely limit empirical or semi-empirical approaches. As a consequence, it became necessary to develop a more profound understanding of the electromagnetic properties of a forest scene and their temporal dynamics through controlled experimentation coupled with theoretical development and verification. The concatenation of various models into a physically-based composite model treating the entire forest scene became the major objective of the study as this is the key to development of a series of robust retrieval algorithms for forest biophysical properties. In order to verify the performance of the component elements of the composite model, a series of controlled laboratory and field experiments were undertaken to: (1) develop techniques to measure the microwave dielectric properties of vegetation; (2) relate the microwave dielectric properties of vegetation to more readily measured characteristics such as density and moisture content; (3) calculate the radar cross-section of leaves, and cylinders; (4) improve backscatter models for rough surfaces; and (5) relate attenuation and phase delays during propagation through canopies to canopy properties. These modeling efforts, as validated by the measurements, were incorporated within a larger model known as the Michigan Microwave Canopy Scattering (MIMICS) Model.

Dobson, M. Craig↗

NASA NTRS1995

Estimating release of carbon from 1990 and 1991 forest fires in Alaska

An improved method to estimate the amounts of carbon released during fires in the boreal forest zone of Alaska in 1990 and 1991 is described. This method divides the state into 64 distinct physiographic regions and estimates areal extent of five different land covers: two forest types, peat land, tundra, and nonvegetated. The areal extent of each cover type was estimated from a review of topographic maps of each region and observations on the distribution of foreat types within the state. Using previous observations and theoretical models for the two forest types found in interior Alaska, models of biomass accumulation as a function of stand age were developed. Stand age distributions for each region were determined using a statistical distribution based on fire frequency, which was from available long-term historical records. Estimates of the degree of biomass combusted were based on recent field observations as well as research reported in the literature. The location and areal extent of fires in this region for 1990 and 1991 were based on both field observations and analysis of satellite (advanced very high resolution radiometer (AVHRR)) data sets. Estimates of average carbon release for the two study years ranged between 2.54 and 3.00 kg/sq m, which are 2.2 to 2.6 times greater than estimates used in other studies of carbon release through biomass burning in boreal forests. Total average annual carbon release for the two years ranged between 0.012 and 0.018 Pg C/yr, with the lower value resulting from the AVHRR estimates of fire location and area.

Kaisischke, Eric S.↗

DOE OSTI2016

Explicit modeling of abiotic and landscape factors reveals precipitation and forests associated with aphid abundance

Abstract Increases in natural or noncrop habitat surrounding agricultural fields have been shown to be correlated with declines in insect crop pests. However, these patterns are highly variable across studies suggesting other important factors, such as abiotic drivers, which are rarely included in landscape models, may also contribute to variability in insect population abundance. The objective of this study was to explicitly account for the contribution of temperature and precipitation, in addition to landscape composition, on the abundance of a widespread insect crop pest, the soybean aphid ( Aphis glycines Matsumura), in Wisconsin soybean fields. We hypothesized that higher soybean aphid abundance would be associated with higher heat accumulation (e.g., growing degree days) and increasing noncrop habitat in the surrounding landscape, due to the presence of the overwintering primary hosts of soybean aphid. To evaluate these hypotheses, we used an ecoinformatics approach that relied on a large dataset collected across Wisconsin over a 9‐year period (2003–2011), for an average of 235 sites per year ( n = 2,110 fields total). We determined surrounding landscape composition (1.5‐km radius) using publicly available satellite‐derived land cover imagery and interpolated daily temperature and precipitation information from the National Weather Service COOP weather station network. We constructed linear mixed models for soybean aphid abundance based on abiotic and landscape explanatory variables and applied model averaging for prediction using an information theoretic framework. Over this broad spatial and temporal extent in Wisconsin, we found that variation in growing season precipitation was positively related to soybean aphid abundance, while higher precipitation during the nongrowing season had a negative effect on aphid populations. Additionally, we found that aphid populations were higher in areas with proportionally more forest but were lower in areas where minor crops, such as small grains, were more prevalent. Thus, our findings support our hypothesis that including abiotic drivers increases our understanding of crop pest abundance and distribution. Moreover, by explicitly modeling abiotic factors, we may be able to explore how variable climate in tandem with land cover patterns may affect current and future insect populations, with potentially critical implications for crop yields and agricultural food webs.

Stack Whitney, Kaitlin↗

DOE OSTI2020

Using a coupled dispersion model to estimate depletion of a tritium oxide plume by a forest

Tritium processing facilities may release tritium oxide (HTO) to the atmosphere which poses potential health risks to exposed co-located workers and to offsite individuals. Most radiological consequence analyses determine HTO dose by applying Gaussian plume models to simulate the transport of HTO. Within these models, deposition velocity is used to assess the sum of all deposition processes acting on the plume. While this may account for vegetative and soil uptake or respiration processes, it may currently lack inclusion of the complex interactions within heterogeneous forested environments. In this complex morphology, dispersion patterns are significantly altered by changing flow regimes above and below the forest canopy and by the transfer of plume material across the canopy boundary. To determine the effects of a heterogeneous forest canopy on an airborne HTO plume, a Gaussian plume model coupled with an advection-diffusion plume model was applied to estimate transport in the free atmosphere above the forest and within the forest canopy and understory. During 2012, wind speed and wind direction measurements taken at 5 heights, ranging from 2-m to 28-m, on an instrumented meteorological tower located in a loblolly pine forest at the Department of Energy (DOE) Savannah River Site (SRS), near Aiken, SC. From these measurements, model predictions were made over a full spectrum of meteorological conditions. Deposition and resuspension velocities were calculated based on the model-predicted flux of plume material across the top of the forest canopy. Additionally, net deposition velocity of the plume material was calculated as the difference between the deposition and resuspension velocities. The 1st and 5th percentile net deposition velocities were estimated to be 0.7 cm s-1 and 1.2 cm s-1, respectively.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

NASA NTRS1992

Modeling the effect of land use on carbon storage in the forests of the Pacific Northwest

There is concern as to how the balance of carbon in the terrestrial ecosystem will change in response to a variety of land use practices. A study is described in which a methodology is being developed to help narrow this uncertainty for the temperate forets of the Pacific Northwest region of the US. A carbon storage model is being developed to respond to forest harvesting, the dominant use of land in the region. By linking the carbon model to satellite imagery and a climate simulation model, the current amount of carbon stored in the forests of the Pacific northwest is estimated. The archive of Landsat multispectral scanner (MSS) images permits a 20-year historical perspective of land use changes in the region. With these data, the recent impact of regional land use in forest carbon stores is assessed.

Cohen, Warren B.↗

DOE OSTI2025

Harvest legacies and climate change interact to shape forest structure and biomass through time

Forest structural change driven by climate trends has been observed worldwide and is expected to increase in the future. Management of forest structure has been an important tool for mitigating the impacts of climate change but forest structure may shift independently of management goals as it interacts with climate change. Here, we investigated the long-term impacts of harvest-based management strategies on structure and resistance to climate-induced biomass loss using a process-based ecosystem model for a midwestern USA hardwood forest. We identified aboveground biomass loss events and compared the cumulative number of these events following a five-year period of active management under four management strategies and two climate change scenarios. Management legacy had the clearest impact on climate-driven biomass loss over the mid-term (∼25 years) with the shelterwood scenario experiencing no loss events during this period. However, by the end of the century legacy effects faded and climate change severity became the driver of differences as greater warming scenarios experienced twice the loss events, and end-of-century loss events were 3–10 times more frequent than mid-century events. We found that while structure was distinct among harvest scenarios through the mid-century, differences were negligible by the end-of-century; identical to when management correlated with loss event frequency. Here, we found that loss events were preceded by a drop in precipitation two years prior, while no specific forest structure preceded a loss event. However, the structures preceding a biomass loss event were distinct across different legacies, implying additional influences of past management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Towards mapping biodiversity from above: Can fusing lidar and hyperspectral remote sensing predict taxonomic, functional, and phylogenetic tree diversity in temperate forests?

Abstract Aim Rapid global change is impacting the diversity of tree species and essential ecosystem functions and services of forests. It is therefore critical to understand and predict how the diversity of tree species is spatially distributed within and among forest biomes. Satellite remote sensing platforms have been used for decades to map forest structure and function but are limited in their capacity to monitor change by their relatively coarse spatial resolution and the complexity of scales at which different dimensions of biodiversity are observed in the field. Recently, airborne remote sensing platforms making use of passive high spectral resolution (i.e., hyperspectral) and active lidar data have been operationalized, providing an opportunity to disentangle how biodiversity patterns vary across space and time from field observations to larger scales. Most studies to date have focused on single sites and/or one sensor type; here we ask how multiple sensor types from the National Ecological Observatory Network’s Airborne Observation Platform (NEON AOP) perform across multiple sites in a single biome at the NEON field plot scale (i.e., 40 m × 40 m). Location Eastern USA. Time period 2017–2018. Taxa studied Trees. Methods With a fusion of hyperspectral and lidar data from the NEON AOP, we assess the ability of high resolution remotely sensed metrics to measure biodiversity variation across eastern US temperate forests. We examine how taxonomic, functional, and phylogenetic measures of alpha diversity vary spatially and assess to what degree remotely sensed metrics correlate with in situ biodiversity metrics. Results Models using estimates of forest function, canopy structure, and topographic diversity performed better than models containing each category alone. Our results show that canopy structural diversity, and not just spectral reflectance, is critical to predicting biodiversity. Main conclusions We found that an approach that jointly leverages spectral properties related to leaf and canopy functional traits and forest health, lidar derived estimates of forest structure, fine‐resolution topographic diversity, and careful consideration of biogeographical differences within and among biomes is needed to accurately map biodiversity variation from above.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

An observational constraint on stomatal function in forests: evaluating coupled carbon and water vapor exchange with carbon isotopes in the Community Land Model (CLM4.5)

Abstract. Land surface models are useful tools to quantify contemporary and future climate impact on terrestrial carbon cycle processes, provided they can be appropriately constrained and tested with observations. Stable carbon isotopes of CO2 offer the potential to improve model representation of the coupled carbon and water cycles because they are strongly influenced by stomatal function. Recently, a representation of stable carbon isotope discrimination was incorporated into the Community Land Model component of the Community Earth System Model. Here, we tested the model's capability to simulate whole-forest isotope discrimination in a subalpine conifer forest at Niwot Ridge, Colorado, USA. We distinguished between isotopic behavior in response to a decrease of δ13C within atmospheric CO2 (Suess effect) vs. photosynthetic discrimination (Δcanopy), by creating a site-customized atmospheric CO2 and δ13C of CO2 time series. We implemented a seasonally varying Vcmax model calibration that best matched site observations of net CO2 carbon exchange, latent heat exchange, and biomass. The model accurately simulated observed δ13C of needle and stem tissue, but underestimated the δ13C of bulk soil carbon by 1–2 ‰. The model overestimated the multiyear (2006–2012) average Δcanopy relative to prior data-based estimates by 2–4 ‰. The amplitude of the average seasonal cycle of Δcanopy (i.e., higher in spring/fall as compared to summer) was correctly modeled but only when using a revised, fully coupled An − gs (net assimilation rate, stomatal conductance) version of the model in contrast to the partially coupled An − gs version used in the default model. The model attributed most of the seasonal variation in discrimination to An, whereas interannual variation in simulated Δcanopy during the summer months was driven by stomatal response to vapor pressure deficit (VPD). The model simulated a 10 % increase in both photosynthetic discrimination and water-use efficiency (WUE) since 1850 which is counter to established relationships between discrimination and WUE. The isotope observations used here to constrain CLM suggest (1) the model overestimated stomatal conductance and (2) the default CLM approach to representing nitrogen limitation (partially coupled model) was not capable of reproducing observed trends in discrimination. These findings demonstrate that isotope observations can provide important information related to stomatal function driven by environmental stress from VPD and nitrogen limitation. Future versions of CLM that incorporate carbon isotope discrimination are likely to benefit from explicit inclusion of mesophyll conductance.

58 GEOSCIENCES↗

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 OSTI2025

Quantifying Improvements in Feedstock Performance Resulting from Forest Concepts' Preprocessing Technology (CRADA Final Report)

We will develop structural mechanics models that capture the evolution of feedstock attributes during the Forest Concepts comminution process. Models will be validated by data provided by Forest Concepts and experiments performed by synergistic projects within the Feedstock Conversion Interface Consortium.

09 BIOMASS FUELS↗

DOE OSTI2026

Characterizing Turbulence at a Forest Edge: Comparing Sub-Filter Scale Turbulence Models in Simulations of Flow over a Canopy

In wildfires, atmospheric turbulence plays a major role in the transfer of turbulent kinetic energy. Understanding how turbulence feeds back into a dynamical system is important, down to the varying small scales of fuel structures (i.e. pine needles, grass). Large eddy simulations (LES) are a common way of numerically representing turbulence. The Smagorinsky model (1963) serves as one of the most studied sub-grid scale representations in LES. In this investigation, the Smagorinsky model was implemented in HIGRAD/FIRETEC, LANL’s coupled fire-atmosphere model. This study was motivated by the need to quantitatively investigate the vorticity budget equation in HIGRAD/FIRETEC. The Smagorinsky turbulent kinetic energy (TKE) was compared to FIRETEC’s 1.5-order TKE eddy-viscosity subgrid-scale model, known as the Linn turbulence model. This was done in simulations of flow over flat terrain with a homogeneous, cuboidal canopy in the center of the domain. Examinations of the modeled vertical TKE profile and turbulent statistics at the leading edge, and throughout the canopy, show that the Smagorinsky model provides comparable results to that of the original closure model posed in FIRETEC.

58 GEOSCIENCES↗