Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Land surface”

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

Citizen science coupled with machine learning to quantify green-blue infrastructure cooling potential in Maricopa County, Arizona

Here, this study investigates the spatiotemporal cooling performance of green and blue infrastructure (GBI) in the Dobson Ranch urban neighborhood in Phoenix, Arizona. We leveraged citizen science near-surface (2 m) air temperature (Tair) measurements to train a highly accurate Tair predicting LightGBM machine learning model (R 2 : 0.986, MAE: 0.251 °C, RMSE: 0.585 °C). On June 16, 2024, the park area exhibited approximately 1 °C cooling effect (relative to the neighborhood mean) during both day and night. In contrast, the nearby artificial lake exhibited a stronger cooling effect of 2.4 °C during the day but a slight warming of 0.3 °C at night. At 00:00, locations 50 m downwind of the park were 0.3 °C warmer than the park, while locations 50 m upwind were 0.8 °C warmer. At 11:00, we observed that the downwind area is 0.8 °C cooler and the upwind area is 0.6 °C warmer—at the same 50 m distances relative to the park. We also observed 1 °C cooler and warmer effects respectively at the same 50 m downwind and upwind locations at 19:00 on June 17, 2024. Our data-driven analysis highlights potential limitations of car-traverse measurements, showing that failure to account for temporal variations during the traverse can lead to overestimation of Tair at night and underestimation during the day. Our analysis also showed only a weak correlation (coefficient: 0.48) between Landsat-derived land surface temperature (LST) and model predicted Tair at the time of the local Landsat overpass (∼11.00). This highlights the potential error of relying solely on LST for human thermal exposure analysis—particularly within the heterogenous built-environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A coupled equilibrium boundary layer model with stable water isotopes and its application to local water recycling

The contribution of evapotranspiration (ET) to regional precipitation, known as “local water recycling”, is a key process in the water cycle that can affect water management. However, the amount of planetary boundary layer (PBL) moisture that arises from ET is highly uncertain due to complex atmosphere and land surface conditions. In this study, an idealized two-layer equilibrium boundary layer model was coupled with a stable water isotope module including HDO and H 2 18 O to constrain PBL growing season water transport processes. The model was validated using turbulent heat fluxes and isotope ratios of water vapor (δ m ) and precipitation (δ P ) measured at a cropland site and a nearby tall tower in the Upper Midwest, United States. The results show that the PBL equilibrium features of δ m and δ P are well-constrained by thermal and moisture equilibrium in the PBL. For this study region, the summer values of rain evaporation fraction (f) and local water recycling ratio (LRR) are estimated to be 0.09 and 0.29 ± 0.12, respectively. The summer LRR values for the years 2006–2010 were 0.35, 0.36, 0.17, 0.29, and 0.29, respectively. The small value of LRR in 2008 corresponded to a drought condition with the lowest precipitation and second lowest ET among the five years. The summer magnitude of the amount effect is –2.8‰ (mm day –1 ) –1 and –0.8‰ (mm day –1 ) –1 for HDO and H 2 18 O, respectively. The local water recycling is identified as a significant factor influencing the continental isotope effect. Cropland has likely changed the regional LRR by –7.6 to 19.5% under different pre-agriculture land use scenarios. In conclusion, the feedback processes revealed here indicate that local water recycling is expected to be weakened under drought conditions, but it will be enhanced if irrigation is applied more intensely with more frequent drought events as the climate continues to warm.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Assessing impacts of selective logging on water, energy, and carbon budgets and ecosystem dynamics in Amazon forests using the Functionally Assembled Terrestrial Ecosystem Simulator

Abstract. Tropical forest degradation from logging, fire, and fragmentation not only alters carbon stocks and carbon fluxes, but also impacts physical land surface properties such as albedo and roughness length. Such impacts are poorly quantified to date due to difficulties in accessing and maintaining observational infrastructures, as well as the lack of proper modeling tools for capturing the interactions among biophysical properties, ecosystem demography, canopy structure, and biogeochemical cycling in tropical forests. As a first step to address these limitations, we implemented a selective logging module into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) by mimicking the ecological, biophysical, and biogeochemical processes following a logging event. The model can specify the timing and aerial extent of logging events, splitting the logged forest patch into disturbed and intact patches; determine the survivorship of cohorts in the disturbed patch; and modifying the biomass and necromass (total mass of coarse woody debris and litter) pools following logging. We parameterized the logging module to reproduce a selective logging experiment at the Tapajós National Forest in Brazil and benchmarked model outputs against available field measurements. Our results suggest that the model permits the coexistence of early and late successional functional types and realistically characterizes the seasonality of water and carbon fluxes and stocks, the forest structure and composition, and the ecosystem succession following disturbance. However, the current version of FATES overestimates water stress in the dry season and therefore fails to capture seasonal variation in latent and sensible heat fluxes. Moreover, we observed a bias towards low stem density and leaf area when compared to observations, suggesting that improvements are needed in both carbon allocation and establishment of trees. The effects of logging were assessed by different logging scenarios to represent reduced impact and conventional logging practices, both with high and low logging intensities. The model simulations suggest that in comparison to old-growth forests the logged forests rapidly recover water and energy fluxes in 1 to 3 years. In contrast, the recovery times for carbon stocks, forest structure, and composition are more than 30 years depending on logging practices and intensity. This study lays the foundation to simulate land use change and forest degradation in FATES, which will be an effective tool to directly represent forest management practices and regeneration in the context of Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Electrical Resistivity Tomography Monitoring of In Situ Soil Flushing at the Hanford 100-K East Area: 100KE Soil Flushing Monitoring

Former operations in the 100-KE Area at the Hanford Site resulted in vadose zone hexavalent chromium that provides a source of groundwater contamination. Previous efforts to remediate vadose zone chromium involved excavation and offsite treatment of contaminated soils. Although much of the chromium was removed, contamination still exists in native soil between the water table and the bottom of the pit, which has since been backfilled. In-situ soil flushing was tested at the 100-KE Area to accelerate the removal of the remaining hexavalent chromium in the vadose zone. Soil flushing works by applying clean flush water at the surface, which mobilizes and transports chromium as it migrates downward to the water table. Once the chromium reaches the water table, it is hydraulically contained, extracted, and treated through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. Therefore, performance monitoring requires and understanding of flush water migration through the vadose zone. It is challenging to comprehensively assess soil flushing performance via borehole access due to the limited volume of investigation afforded by a given wellbore. As an alternative, 3D time-lapse electrical resistivity tomography (ERT) was tested as a method of monitoring the distribution of flush water over time within the vadose zone. ERT is a method of imaging the bulk electrical conductivity of the subsurface, which is highly dependent on water saturation levels in unconsolidated and unsaturated sediments. Thus, the timing and location of changes in bulk conductivity caused by soil flushing can be used to infer flush water migration pathways. In this report we present results of 3D time-lapse ERT imaging during two separate soil flushing campaigns conducted during the spring and summer months of 2022 and 2023. Results show generally that: 1) Pit backfill materials appear to nominally have larger permeability than native soils. Consequently, the boundary between pit backfill and native soil had a significant impact on flush water migration, causing some flush water to migrate along the pit boundary to the bottom of the pit. 2) Non-uniform flows, likely caused by variations in hydrogeologic properties, developed in the pit backfill materials, resulting in uneven flush water distribution on the southern margin of the soil flushing zone. 3) Redistribution of water at the interface between backfill materials and Hanford formation materials likely facilitated elevated and uniform distribution of flush water within native soils beneath the deeper parts of the pit boundary beneath the center of the flush zone, which presumably overlie soils with elevated chromium contamination. These areas appear to have been infiltrated by higher volumes of flush water than the northern and southern margins of the flush zone. 4) In comparison to 2022, high flush water application rates significantly improved flush water distribution throughout the target flushing zone. 5) Imaging resolution was limited to a depth of approximate 20 meters below the land surface, due primarily to limitations on the lateral extent of the surface ERT array. The ~10 m region of the vadose zone between approximately 130 meters elevation and the water table at approximately 120 meters elevation was unresolved.

100 Area↗
DOE OSTI2022

Variability in observed stable water isotopes in snowpack across a mountainous watershed in Colorado

In this study, isotopic information from 81 snowpits was collected over a 5-year period in a large, Colorado watershed. Data spans gradients in elevation, aspect, vegetation, and seasonal climate. They are combined with overlapping campaigns for water isotopes in precipitation and snowmelt, and a land-surface model for detailed estimates of snowfall and climate at sample locations. Snowfall isotopic inputs, describe the majority of δ18O snowpack variability. Aspect is a secondary control, with slightly more enriched conditions on east and north facing slopes. This is attributed to preservation of seasonally enriched snowfall and vapour loss in the early winter. Sublimation, expressed by decreases in snowpack d-excess in comparison to snowfall contributions, increases at low elevation and when seasonal temperature and solar radiation are high. At peak snow accumulation, post-depositional fractionation appears to occur in the top 25 ± 14% of the snowpack due to melt-freeze redistribution of lighter isotopes deeper into the snowpack and vapour loss to the atmosphere during intermittent periods of low relative humidity and high windspeed. Relative depth of fractionation increases when winter daytime temperatures are high and winter precipitation is low. Once isothermal, snowpack isotopic homogenization and enrichment was observed with initial snowmelt isotopically depleted in comparison to snowpack and enriching over time. The rate of δ18O increase (d-excess decrease) in snowmelt was 0.02‰ per day per 100-m elevation loss. Isotopic data suggests elevation dictates snowpack and snowmelt evolution by controlling early snow persistence (or absence), isotopic lapse rates in precipitation and the ratio of energy to snow availability. Hydrologic tracer studies using stable water isotopes in basins of large topographic relief will require adjustment for these elevational controls to properly constrain stream water sourcing from snowmelt.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sensitivity of Fine‐Resolution Urban Heat Island Simulations to Soil Moisture Parameterization

ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.

Talebpour, Mahdad↗
DOE OSTI2020

Plant‐mediated methane transport in emergent and floating‐leaved species of a temperate freshwater mineral‐soil wetland

Abstract Methane flux from freshwater mineral‐soil (FWMS) wetlands and its variability among sites is largely modulated by plant‐mediated transport. However, plant‐mediated transport processes are rarely resolved in land surface models and are poorly parametrized for plants commonly found in FWMS wetlands. Here, relationships between methane flux and CO 2 uptake, as well as plant conductance of methane were evaluated for three plant species and two characteristic functional types: emergent (narrow‐leaved cattail) and floating‐leaved (American lotus and water lily). We found significant but contrasting correlations between methane flux and CO 2 uptake in cattails ( r 2 = 0.51, slope = −0.16, during morning) and water lily ( r 2 = 0.32, slope = 0.064, after midday). This relationship was not significant in American lotus, showing that stomata regulation of methane fluxes is species‐specific and not generalizable across the floating‐leaved plant functional type. Conductance of methane per leaf area showed distinct seasonal dynamics across species. Conductance was similar among the floating‐leaved species (6.2 × 10 −3 m d −1 in lotus and 7.2 × 10 −3 m d −1 in water lily) and higher than conductance in the emergent species (2.7 × 10 −3 m d −1 ). Our results provide direct observations of plant conductance rates and identify the vegetation parameters (leaf area, stomatal conductance) that modify them. Our results further suggest that models of methane emissions from FWMS should parameterize plant‐mediated transport in different plant functional types, scaled by leaf area and with variable seasonal phenological dynamics, and consider possible species‐specific mechanisms that control methane transport through plants.

Villa, Jorge A.↗
DOE OSTI2026

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗
DOE OSTI2021

Heterogeneous multiphase flow properties of volcanic rocks and implications for noble gas transport from underground nuclear explosions

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2021

Developing a strategy for the national coordinated soil moisture monitoring network

Soil moisture is a critical land surface variable, affecting a wide variety of climatological, agricultural, and hydrological processes. Determining the current soil moisture status is possible via a variety of methods, including in situ monitoring, remote sensing, and numerical modeling. Although all of these approaches are rapidly evolving, there is no cohesive strategy or framework to integrate these diverse information sources to develop and disseminate coordinated national soil moisture products that will improve our ability to understand climate variability. The National Coordinated Soil Moisture Monitoring Network initiative has developed a national strategy for network coordination with NOAA's National Integrated Drought Information System. The strategy is currently in review within NOAA, and work is underway to implement the initial milestones of the strategy. This update reviews the goals and steps being taken to establish this national-scale coordination for soil moisture monitoring in the United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Flexible User-Defined Domain Decomposition in Kilometer-Scale E3SM Land Model Simulation

The Energy Exascale Earth System Model (E3SM) Land Model (ELM) has been extended to kilometer-scale (km-ELM) resolutions, enabling high-fidelity simulations of terrestrial processes at 1 km x 1 km grid spacing. In ELM, domain decomposition partitions the computational domain across processors, ensuring efficient parallel execution. Currently, round-robin decomposition is applied, providing a straightforward way to distribute computational workload. As ELM continues evolving at the kilometer-scale (km-scale), particularly with integrating lateral flow modeling, decomposition strategies must also account for the increased workload and data movement. This paper introduces a flexible user-defined domain decomposition framework, allowing users to customize domain partitioning based on application requirements. The impact of different decomposition strategies is evaluated across various applications concerning computation, communication, and I/O. Results demonstrate that while 1D partitioning yields superior I/O performance, k-nearest neighbors (KNN) clustering effectively reduces inter-process communication overhead. This study lays the groundwork for scalable partitioning in large-scale land surface simulations, enhancing next-generation Earth system modeling.

Wang, Dali [ORNL] (ORCID:0000000168065108)↗
DOE OSTI2022

Improving Simulations of Vegetation Dynamics over the Tibetan Plateau: Role of Atmospheric Forcing Data and Spatial Resolution

The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC m –2 yr –1 to 241.21 gC m –2 yr –1 with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Whole-Ecosystem Warming Increases Plant-Available Nitrogen and Phosphorus in an Ombrotrophic Bog

Warming is expected to increase the net release of carbon from peatland soils, contributing to future warming. This positive feedback may be moderated by the response of peatland vegetation to rising atmospheric [CO 2 ] or to increased soil nutrient availability. We asked whether a gradient of whole-ecosystem warming (from + 0 °C to + 9 °C) would increase plant-available nitrogen and phosphorus in an ombrotrophic bog in northern Minnesota, USA, and whether elevated [CO 2 ] would modify the nutrient response. We tracked changes in plant-available nutrients across space and through time and in comparison with other nutrient pools, and assessed whether nutrient warming responses were captured by a point version of the land-surface model, ELM-SPRUCE. We found that warming exponentially increased plant-available ammonium and phosphate, but that nutrient dynamics were unaffected by elevated [CO 2 ]. The warming response increased by an order of magnitude between the first and fourth year of the experimental manipulation, perhaps because of dramatic mortality of Sphagnum mosses in the surface peat of the warmest treatments. However, neither the magnitude nor the temporal dynamics of the responses were captured by ELM-SPRUCE. Here, relative increases in plant-available ammonium and phosphate with warming were similar, but the response varied across raised hummocks and depressed hollows and with peat depth. Plant-available nutrient dynamics were only loosely correlated with inorganic and organic porewater nutrients, likely representing different processes. Future predictions of peatland nutrient availability under climate change scenarios must account for dynamic changes in nutrient acquisition by plants and microbes, as well as microtopography and peat depth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Microbial abundances and carbon use under ambient temperature or experimental warming in a southern boreal peatland

Abstract Organic peat soils occupy relatively little of the global land surface area but store vast amounts of soil carbon in northern latitudes where climate is warming at a rapid pace. Warming may result in strong positive feedbacks of carbon loss and global climate change driven by microbial processes if warming alters the balance between primary productivity and decomposition. To elucidate effects of warming on the microbial communities mediating peat carbon dynamics, we explored the abundance of broad microbial groups and their source of carbon (i.e. old carbon versus more recently fixed photosynthate) using microbial lipid analysis (δ 13 C PLFA) of peat samples under ambient temperatures and before/after initiation of experimental peat warming (+ 2.25, + 4.5, + 6.75, and + 9 °C). This analysis occurred over a profile to 2 m depth in an undrained, ombrotrophic peat bog in northern Minnesota. We found that the total microbial biomass and individual indicator lipid abundances were stratified by depth and strongly correlated to temperature under ambient conditions. However, under experimental warming, statistically significant effects of temperature on the microbial community were sporadic and inconsistent. For example, 3 months after experimental warming the relative abundance of Gram-negative bacterial indicators across depth combined and > 50 cm depth and Gram-positive bacterial indicators at 20–50 cm depth showed significant positive relationships to temperature. At that same timepoint, however, the relative abundance of Actinobacterial indicators across depth showed a significant negative relationship to temperature. After 10 months of experimental warming, the relative abundance of fungal biomarkers was positively related to temperature in all depths combined, and the absolute abundance of anaerobic bacteria declined with increasing temperature in the 20–50 cm depth interval. The lack of observed response in the broader microbial community may suggest that at least initially, microbial community structure with peat depth in these peatlands is driven more by bulk density and soil water content than temperature. Alternatively, the lack of broad microbial community response may simply represent a lag period, with more change to come in the future. The long-term trajectory of microbial response to warming in this ecosystem then could either be direct, after this initial lag time, or indirect through other physical or biogeochemical changes in the peat profile. These initial results provide an important baseline against which to measure long-term microbial community and carbon-cycling responses to warming and elevated CO 2 .

Felice, Mark↗
DOE OSTI2025

Evaluating ecosystem water use efficiency and recovery dynamics during flash droughts: insights from observations and model simulations

Flash droughts (FD), rapidly emerging in a warming future, disrupt ecosystems, agriculture, and water security. Ecosystem water use efficiency (WUE), the ratio of gross primary production (GPP) to actual evapotranspiration (AET), balances carbon assimilation and water loss. FD rapidly disrupts this balance, making WUE critical for assessing plant stress and recovery. Here, this study investigates the dynamics of landscape-scale WUE, and the components of GPP and AET under FD utilizing both observed data from the Missouri Ozark AmeriFlux site (US-MOz) and version 2 of the U.S. Department of Energy’s Earth, Energy, Exascale System Model (E3SM) Land Model (ELMv2). Observations and simulations reveal GPP as dominant for WUE during earlier FD events (2005, 2007, 2012), shifting to AET in recent events (2014, 2018). This agreement indicates that the ELM can capture the shifting dynamics of GPP and AET in regulating WUE under FD conditions. However, the ELM systematically underestimates both GPP and AET and does so in a manner that does not preserve their ratio. As a result, WUE is also underestimated, suggesting that GPP is more strongly underestimated than AET. Furthermore, the ELM also underestimates the speed of GPP recovery, producing an artificially prolonged GPP recovery time following FD events. Observed environmental drivers such as vapor pressure deficit (VPD), soil moisture (SM), and predawn leaf water potential (PLWP) effectively predict WUE, but ELM primarily highlights SM, underestimating VPD’s role. This study demonstrates that relying solely on soil moisture fails to capture the rapid hydraulic recovery observed in PLWP, underscoring the necessity of integrating plant hydraulics into land surface models to improve flash drought predictability.

Evapotranspiration↗
DOE OSTI2021

Daily grass reference evapotranspiration with Meteosat Second Generation shortwave radiation and reference ET products

Here we assess the accuracy of estimating daily grass reference evapotranspiration (PM-ET o ) using daily shortwave radiation (R s ) and reference evapotranspiration (ET REF ) products provided by the Meteosat Second Generation (MSG) geostationary satellite delivered by the Satellite Applications Facility on Land Surface Analysis (LSA-SAF) framework. The accuracy of using reanalysis ERA5 shortwave radiation data (R s ERA5 ) provided by the European Center for Medium-Range Weather Forecasts (ECMWF) is also evaluated. The assessments were performed using observed weather variables at 37 weather stations distributed across continental Portugal, where climate conditions range from semi-arid to humid, and 12 weather stations located in Azores islands, characterized by humid, windy and often cloudy conditions. This study&#x27;s use of data from a variety of climate conditions contributed to a unique and innovative assessment of the usability of LSA-SAF and ERA5 products for ET o estimation. The first assessment focused on comparing LSA-SAF estimates of R s (R s LSA-SAF ) against ground stations (R s ground ). The results showed a good matching between the two R s data sets for continental Portugal but a tendency for R s LSA-SAF to under-estimate R s ground in the cloudy islands of Azores. ET o values computed using R s LSA-SAF data and observed temperature, humidity and wind speed (ET o LSA-SAF ) were then compared with PM-ET o estimates with ground-based data, which were used as benchmark; input data of temperature and humidity needed for PM-ET o were quality checked for surface aridity effects. It was observed that ET o LSA-SAF is strongly correlated with PM-ET o (R 2 > 0.97) for most locations in continental Portugal, with regression coefficient of a linear regression forced to the origin ranging between 0.95 and 1.05, mean root mean square error (RMSE) of 0.13 mm d –1 , and Nash and Sutcliff efficiency of modeling (EF) above 0.95. For most Azores locations, ET o LSA-SAF over-estimated PM-ET o . This is likely a consequence of the high spatio-temporal heterogeneity of weather conditions that occur in these oceanic islands together with the different footprints of satellite (averaged over the pixel) and station observations. Reanalysis ERA5 shortwave radiation data presented similar behavior to the LSA-SAF products, however with slightly lower accuracy. The daily LSA-SAF ET REF product (ET REF LSA-SAF ) was assessed and results have shown a good accuracy of this product, with acceptable RMSE and high EF values, for continental Portugal but a low accuracy for the Azores islands. A simplified bias correction approach was shown to improve both ET o derived from the LSA-SAF products, namely for Azores stations, which seem to be representative of smaller areas. The use of the FAO-PM temperature approach (PMT) was also assessed using the R s LSA-SAF and R s ERA5 data, which showed a superiority of the LSA-SAF product for ET o estimations (ET o PMT LSA-SAF ). No significant differences (p < 0.05) were observed in terms of the median value of the RMSE when adopting ET o PMT and ET REF LSA-SAF . Differently, results showed that using the R s LSA-SAF in the PMT approach (ET o PMT LSA-SAF ) produces significantly better RMSE results than ET o PMT and ET REF LSA-SAF . Overall, the performed assessment allows concluding that the use of R s LSA-SAF, and to a lesser extent the use of the R s ERA5 , highly improves the accuracy of computation of ET o when R s observations are not available, including when only temperature data are accessible. The use of the ET REF LSA-SAF product is a good alternative when observed weather data are not available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mass-conserving implicit–explicit methods for coupled compressible Navier–Stokes equations

Earth system models are composed of coupled components that separately model systems such as the global atmosphere, ocean, and land surface. While these components are well developed, coupling them in a single system can be a significant challenge. Computational efficiency, accuracy, and stability are principal concerns. In this study we focus on these issues. In particular, implicit–explicit (IMEX) tight and loose coupling strategies are explored for handling different time scales. For a simplified model for the air–sea interaction problem, we consider coupled compressible Navier–Stokes equations with an interface condition. Under the rigid-lid assumption, horizontal momentum and heat flux are exchanged through the interface. Several numerical experiments are presented to demonstrate the stability of the coupling schemes. Here, we show both numerically and theoretically that our IMEX coupling methods are mass conservative for a coupled compressible Navier–Stokes system with the rigid-lid condition

42 ENGINEERING↗
DOE OSTI2023

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗