Search NASA⌕ Search

Engineering topics

Fang, Yilin

Publications and source records attributed to Fang, Yilin.

15 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A unified ensemble soil moisture dataset across the continental United States

Abstract A unified ensemble soil moisture (SM) package has been developed over the Continental United States (CONUS). The data package includes 19 products from land surface models, remote sensing, reanalysis, and machine learning models. All datasets are unified to a 0.25-degree and monthly spatiotemporal resolution, providing a comprehensive view of surface SM dynamics. The statistical analysis of the datasets leverages the Koppen-Geiger Climate Classification to explore surface SM’s spatiotemporal variabilities. The extracted SM characteristics highlight distinct patterns, with the western CONUS showing larger coefficient of variation values and the eastern CONUS exhibiting higher SM values. Remote sensing datasets tend to be drier, while reanalysis products present wetter conditions. In-situ SM observations serve as the basis for wavelet power spectrum analyses to explain discrepancies in temporal scales across datasets facilitating daily SM records. This study provides a comprehensive soil moisture data package and an analysis framework that can be used for Earth system model evaluations and uncertainty quantification, quantifying drought impacts and land–atmosphere interactions and making recommendations for drought response planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An ML-based terrestrial data fusion and augmentation framework to enable advanced understanding of the terrestrial carbon and water interactions

Soil moisture is essential to the terrestrial carbon and water cycles and land–atmosphere interactions. There are various types of soil moisture data, and each type has the distinct spatiotemporal strengths and limitations, depending on the diverse applications and retrieval methodologies of different data types (Li et al., in review; The PNNL-82151 FY23 Report). However, the limitations of different soil moisture data in terms of accuracy and spatiotemporal coverage hinder our ability to further understand the soil moisture dynamics across scales. To have a gap free soil moisture data product with a fine spatiotemporal coverage and vertical profiles, we train extreme gradient boosting (XGBoost) models by using (1) in-situ soil moisture measurements from the International Soil Moisture Network (ISMN), (2) soil moisture from the ECMWF reanalysis (ERA) at the 9 km and sub-daily spatiotemporal resolution, (3) the Daymet meteorological fields, and (4) data products that characterize surface conditions, including soil texture, organic content, topography, vegetation type, and rooting depth. We use the trained XGBoost models that have consistent performance across seven soil layers, i.e., 0–5 cm, 5–10 cm, 10–20 cm, 20–40 cm, 40–60 cm, 60–100 cm, and 100–200 cm, and the gridded model predictors to generate a soil moisture data at the 1 km and daily spatiotemporal resolution for the Continental United States (CONUS) from 2001–2020. This dataset can be broadly used for Earth system model benchmark, monitoring extreme weathers, making informed decisions regarding agriculture, water resource management, climate change mitigation, and ecosystem preservation.

58 GEOSCIENCES↗
DOE OSTI2024

On the transferability of residence time distributions in two 10-km long river sections with similar hydromorphic units

Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface and their residence time distributions (RTDs) in the subsurface are important for managing the water quality and ecosystem health in dynamic river corridors. However, direct simulating high-spatial resolution HEFs and RTDs can be time-consuming, especially for watershed-scale modeling. Efficient surrogate models linking RTDs to hydromorphic units (HUs) can be alternatives for simulating RTDs in large-scale models. A common concern of these surrogate models, though, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this issue, this work evaluates the HEFs and resulting RTD-HU relationships for two 10-km long river corridors along the Columbia River leveraging a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework we previously developed. Applying such a framework at the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. Finally, our comparison shows that the similarity and transferability of the RTD-HU relationship is very low for the two investigated river sections, which suggests that devising a general algorithm to estimate RTDs based solely on surface water hydrodynamics and short-distance river channel topography data, as well as HU classification, might be nearly impossible.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying the Impacts of Land-Cover Change on the Hydrologic Response to Hurricane Ida in the Lower Mississippi River Basin

Abstract The Lower Mississippi River basin (LMRB) has experienced significant changes in land cover and is one of the most vulnerable regions to hurricanes in the United States. Here, we study the impacts of land-cover change on the hydrologic response to Hurricane Ida in LMRB. By using an integrated surface–subsurface hydrologic model, Energy Exascale Earth System Model (E3SM) Land Model coupled with the three-dimensional ParFlow subsurface flow model (ELM-ParFlow), we simulate the effects of land-cover change on the flood volume and peak timing induced by rainfall from Hurricane Ida. The results show that land-cover changes from 1850 to 2015, which resulted in a smoother surface and less vegetation, exacerbated both flood peak time and volume induced by Hurricane Ida. The effects of land-cover changes can be decomposed into two mechanisms: a smoother surface routes more water faster to a watershed outlet and less vegetation allows more water to contribute to surface runoff. By comparing scenarios in which the two mechanisms were isolated, we found that changes in soil moisture due to vegetation cover change have more dominant effects on floods in the southern part and changes in Manning’s coefficient have the largest effect on floods in the northern part of the LMRB. The study provides important insights into the complex relationship between land-use, land-cover, and hydrologic processes in coastal regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data and scripts associated with a manuscript on residence time distribution simulation in two 10-kilometer long river sections

This data package is associated with the publication “On the Transferability of Residence Time Distributions in Two 10-km Long River Sections with Similar Hydromorphic Units” submitted to the Journal of Hydrology (Bao et al. 2024).Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface, along with their residence time distributions (RTDs) in the subsurface, is crucial for managing water quality and ecosystem health in dynamic river corridors. However, directly simulating high-spatial resolution HEFs and RTDs can be a time-consuming process, particularly for watershed-scale modeling. Efficient surrogate models that link RTDs to hydromorphic units (HUs) may serve as alternatives for simulating RTDs in large-scale models. One common concern with these surrogate models, however, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this, we evaluated the HEFs and the resulting RTD-HU relationships for two 10-kilometer-long river corridors along the Columbia River, using a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework that we previously developed. Applying this framework to the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. This data package includes the model inputs files and the simulation results data. This data package contains 10 folders. The modeling simulation results data are in the folders 100H_pt_data and 300area_pt_data, for the study domain Hanford 100H and 300 area respectively. The remaining eight folders contain the scripts and data to generate the manuscript figures. The file-level metadata file (Bao_2024_Residence_Time_Distribution _flmd.csv) includes a list of all files contained in this data package and descriptions for each. The data dictionary file (Bao_2024_Residence_Time_Distribution _dd.csv) includes column header definitions and units of all tabular files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Northern Hemisphere Snow Drought in Earth System Model Simulations and ERA5‐Land Data in 1980–2014

Abstract Low snow levels over the past few decades and predictions of a low‐to‐no snow future have spurred research into snow droughts, which pose a threat to water security and management. Systematic data‐model comparisons of snow drought have been lacking, hindering our understanding of the drivers of snow drought in the past. To address this gap, we analyzed snow drought events using standardized snow water equivalent index derived from monthly results of four numerical experiments using the E3SM Land Model (ELM) and ERA5‐Land data during the period of 1980–2014. Additionally, we compared snow drought duration calculated from models with those from the ERA5‐Land data during selected El Niño‐Southern Oscillation (ENSO) years. The numerical experiments were conducted with ELM driven by two prescribed atmospheric forcings, and with the coupled land‐atmosphere configuration of E3SM with and without plant hydraulics scheme feedback. Analysis reveals that 20%–30% of snow droughts occur due to factors other than above‐normal temperature and low snowfall, such as low soil moisture, warm soil temperature, and low relative humidity, etc., especially in high latitudes (50° North). Furthermore, our study highlights the exacerbating effect of ENSO events on snow drought conditions in various regions, despite some discrepancies between model and ERA5‐Land results. We also identified limitations of the coupled land‐atmosphere models in our current configuration in capturing the spatial patterns of snow droughts. This study underscores the challenge of predicting and mitigating snow drought and the need for a comprehensive understanding of the factors contributing to snow drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluation of Multi-Fidelity Soil Moisture Products Across the Continental United States

We have aggregated the most recent soil moisture datasets from a diverse range of sources, encompassing the Continental United States (CONUS). These sources encompass gridded data from remote sensing products, reanalysis products, machine learning-based projects, and land surface modeling products. Additionally, we have obtained and processed in-situ soil moisture observations from the International Soil Moisture Network. The collected datasets exhibit variations in both temporal and spatial resolutions. Among the 20 datasets, six are available at a spatial resolution of 0.25 degrees, while three are at a coarser spatial resolution of 25 km. To minimize spatial interpolation, we conducted data uncertainty evaluations at the 0.25-degree spatial resolution. For our data evaluations, we maintained a monthly temporal resolution, which effectively captures soil moisture seasonality and interannual variability. Our data processing strategy preserves the raw data and interpolated data at their original temporal resolutions. Datasets with higher temporal resolutions, including daily, three-hourly, and hourly datasets, are set aside for subsequent analyses. These analyses will delve into topics such as soil moisture changes and recovery during extreme weather events. Furthermore, we have processed auxiliary data to enhance our evaluation, leveraging tools such as Google Earth Engine. This includes incorporating topography data, land use land cover data, Köppen-Geiger climate classification, and more to provide a comprehensive assessment from multiple sources.

Li, Lingcheng↗
DOE OSTI2023

Resolving dynamic mineral-organic interactions in the rhizosphere by combining in-situ microsensors with plant-soil reactive transport modeling

Associations between minerals and organic matter represent one of the most important carbon storage mechanisms in soils. Plant roots are major sources of soil carbon, and resolving the dynamics and dominance of microbial consumption versus mineral sorption of root-derived carbon is critical to understanding soil carbon storage. Here we integrate in-situ rhizosphere microsensor and plant physiological measurements with a 3-D plant-soil reactive transport model to explore the fate of dissolved organic carbon (DOC) in the rhizosphere, particularly its microbial consumption and interaction with Fe oxide minerals. Over several days, a microdialysis probe sampling pore water at the root-soil interface of growing Vicia faba roots in live soil, revealed clear diel patterns of DOC concentration. Daytime DOC spikes coincided with peaks in leaf-level photosynthesis rates and were accompanied by declining redox potential and dissolved oxygen as well as increasing pH in the rhizosphere. Incorporating microsensor data into our modeling framework showed that the measured rapid loss of DOC after each mid-day spike could not be explained by consumption via aerobic respiration, nor via anaerobic respiration dominated by Fe oxide reduction. Rather, in the model, a large fraction of rhizosphere DOC was rapidly immobilized each day by adsorption to Fe oxides. Further, modeled microbial Fe reduction (fueled by DOC) did not mobilize significant organic carbon from Fe oxides during the day. Instead, the model predicted equilibrium desorption of organic carbon from Fe oxides at night. This new mechanistic modeling framework, which couples aboveground plant physiological measurements with non-destructive high-resolution monitoring of rhizosphere processes, has great potential for exploring the dynamics and balance of the various microbial reactions and mineral interactions controlling carbon transformations and storage in soils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
DOE OSTI2023

Modeling the mechanisms of conifer mortality under seawater exposure

Summary Relative sea level rise (SLR) increasingly impacts coastal ecosystems through the formation of ghost forests. To predict the future of coastal ecosystems under SLR and changing climate, it is important to understand the physiological mechanisms underlying coastal tree mortality and to integrate this knowledge into dynamic vegetation models. We incorporate the physiological effect of salinity and hypoxia in a dynamic vegetation model in the Earth system land model, and used the model to investigate the mechanisms of mortality of conifer forests on the west and east coast sites of USA, where trees experience different form of sea water exposure. Simulations suggest similar physiological mechanisms can result in different mortality patterns. At the east coast site that experienced severe increases in seawater exposure, trees loose photosynthetic capacity and roots rapidly, and both storage carbon and hydraulic conductance decrease significantly within a year. Over time, further consumption of storage carbon that leads to carbon starvation dominates mortality. At the west coast site that gradually exposed to seawater through SLR, hydraulic failure dominates mortality because root loss impacts on conductance are greater than the degree of storage carbon depletion. Measurements and modeling focused on understanding the physiological mechanisms of mortality is critical to reducing predictive uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Part I: Predicting performance of Purolite A532E resins for remediation of comingled contaminants in groundwater

Ion exchange (IX) resins are used in pump-and-treat (P&T) facilities to remove soluble groundwater contaminants. However, natural anions present at concentrations orders of magnitude higher than contaminants can compete for IX sites and impact resin lifecycles. Here, the Hanford Site’s 200 West Area P&T facility (Washington State, USA) was selected as a case study because it currently uses two IX resins: Purolite® A532E (A532E) to remove pertechnetate (TcO 4 - ) and DOWEX 21K (DOWEX) to remove uranium from groundwater. Nitrate (NO 3 - ), sulfate (SO 4 2- ), chloride (Cl - ), and carbonate (CO 3 2- ) anions have been identified to potentially compete for A532E and DOWEX IX sites. Hanford-relevant anion groundwater concentrations were used to design a series of laboratory-scale batch experiments to evaluate the impact of competing anions on resin performance and potential kinetic effects. These data are then modeled to obtain Cl--normalized equilibrium exchange coefficients (K) needed to predict IX resin performance. The work is presented in two parts, with IX performance evaluated for A532E in Part I and DOWEX in Part II. Part I results demonstrate that TcO 4 - uptake is not impacted by NO 3 - , SO 4 2- , Cl - , CO 3 2- (as HCO 3 - ) and U(VI) carbonate anions, with K TcO4-/Cl- > 4,000, likely due to the high selectivity of A532E trihexylammonium sites for the large, weakly hydrated TcO 4 - anion. Other anion K values were K NO3-/Cl- = 20, K SO4--/Cl- = 0.2, K HCO3-/Cl- = 0.09, K U/Cl- = 370–1000. These K values provide conservative parameters for predicting A532E performance, and demonstrate that, under these test conditions, A532E will remove TcO 4 - from current and future influent streams to meet groundwater treatment objectives.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Ecohydrological decoupling under changing disturbances and climate

Terrestrial disturbances are increasing in frequency and severity, perturbing the hydrologic cycle by altering vegetation-mediated water use and microclimate. Changes in the hydrological cycle alter vegetation succession, which is additionally influenced by changing climate and atmospheric CO 2 . Here, the interacting drivers and feedbacks between disturbances, climate, and vegetation cause uncertainty regarding the sustainable provision of freshwater through streamflow. Here we synthesize the literature on post-disturbance ecohydrological coupling, in particular, including the mechanistic relationship between vegetation and streamflow, under changing disturbance regimes, atmospheric CO 2 , and climate. Disturbance can cause decoupling between transpiration and streamflow by altering the connectivity, size, availability, and spatial distribution of their source pools. Decoupling increases when soil moisture is low due to reduced connectivity, and when soil moisture is high and Q and T are independent. Successional trajectories regulate the physiological and physical features influencing influence the dynamics of source water partitioning. Changing climate and disturbance regimes can alter succession and prolong decoupling, with potential moderation by CO 2 . Increasing rates, severity, and spread of disturbances along with warming could promote greater decoupling globally. From this review emerges a framework of testable hypotheses that identify the critical processes regulating ecohydrological coupling and provides a roadmap for future research. Accurate prediction of post-disturbance coupling and subsequent water availability for human consumption requires understanding the degree of hydraulic connectivity between source water pools for transpiration and streamflow and their response to succession under changing disturbance and climate regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Part II: Predicting performance of $\mathrm{DOWEX 21K}$ resin for remediation of comingled contaminants in groundwater

The selectivity of ion exchange (IX) resins for aqueous contaminant removal can be impacted by changing concentrations of competing natural groundwater ions. In a two-part investigation, the Hanford Site 200 West Area pump-and-treat (P&T) facility in Washington State, USA is used as a case study to evaluate the performance of two IX resins for groundwater treatment: Purolite® A532E for pertechnetate (TcO 4 - ) removal, explored in Part I, and DOWEX 21K (DOWEX) for uranium (U) removal. In Part II, DOWEX selectivity for U, as uranyl carbonate species, and uptake kinetics is quantified in a series of laboratory-scale aqueous batch experiments containing Hanford-relevant concentrations of competing anions nitrate (NO 3 - ), sulfate (SO 4 2- ), chloride (Cl - ), and carbonate (CO 3 2- ), as well as co-mingled contaminant TcO 4 - . Here the results demonstrate that DOWEX trimethylammonium functional groups are highly selective for U carbonate species (85–100 % uptake) under all conditions investigated. Only NO 3 - concentrations of 100 mM were shown to decrease U removal, with the extent (85–99 %) depending on competing anion concentrations present in solution. However, at the highest NO 3 - concentrations reported for groundwaters treated at the P&T facility (25 mM), the effect on U uptake is minimal. The batch sorption results are modeled to obtain chloride normalized equilibrium exchange coefficients (K) for predicting DOWEX performance: K SO4--/Cl- = 2.0, K NO3-/Cl- = 5.0, K HCO3-/Cl- = 1.5, K TcO4-/Cl- = 2,000, and K U/Cl- = 50,000. These K values predict little effect of current and future influent chemistries on U removal by DOWEX, where both uranyl carbonate species and TcO 4 - are removed such that effluent concentrations meet groundwater treatment requirements.

54 ENVIRONMENTAL SCIENCES↗