Search NASA⌕ Search

SEARCH · Search NASA

Results for “atmosphere: 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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
DOE OSTI2020

Surface Chemistry of Liquid Bismuth under Oxygen and Water Vapor Studied by Ambient Pressure X-ray Photoelectron Spectroscopy

Bismuth (Bi) has a fairly low melting point of 544 K making it practical as a liquid metal medium in a number of applications. Under ambient atmospheric conditions the surface of solid Bi oxidizes. While the solid Bi interface has been fairly well characterized with surface science studies upon exposure to oxygen, to date little is known about the molecular level reactivity of the liquid interface. Here, using ambient pressure X-ray photoelectron spectroscopy the liquid-gas interfacial chemistry of Bi was examined upon exposure to oxygen gas and water vapor at 550 K up to a maximum pressure of 0.5 Torr. Water vapor remained unreactive towards the liquid interface over the entire pressure range. Oxygen remined unreactive up to 10 –4 Torr, whereas above this pressure oxidation was observed forming Bi 2 O 3 . The oxidation exposure onset was > 5 x 10 5 Langmuirs, significantly higher than what is required for solid Bi interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

An ultrahigh-resolution E3SM land model simulation framework and its first application to the Seward Peninsula in Alaska

The availability of supercomputers and state-of-science datasets has made it possible to conduct large-scale land simulations at an ultrahigh-resolution. This study reported a computational framework for land surface simulation using the E3SM land model (ELM) at an unprecedented resolution (1 km x 1 km gridcell). The ultrahigh-resolution ELM (uELM) simulation framework includes three parts: (1) high-resolution atmospheric forcing and surface properties dataset generation, (2) massive gridcell-based simulation, and (3) large-scale simulation results analysis. Additionally, we implemented the uELM simulation framework and completed the first 1 km x 1 km terrestrial ecosystem simulation (from 1850 to 2014) over the Seward Peninsula in Alaska (78,000 km 2 ). The experiment contained two phases: a spin-up simulation and a transient simulation, and required five weeks of calculations using 320 cores in a 44-node Linux HPC computer. It created approximately 1.3 TB of data from the transient simulation alone (1850 - present). We selected sample results (monthly and daily simulation outputs) to illustrate the temporal and spatial variations of several variables in high-latitude Arctic ecosystems’ water, energy, and carbon cycles. At last, we summarized the lessons learned and proposed new developments for full-scale uELM simulations over the entire North American continent (approximately 22,000,000 km 2 ).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluating Large‐Storm Dominance in High‐Resolution GCMs and Observations Across the Western Contiguous United States

Abstract Extreme precipitation events are projected to increase in frequency across much of the land‐surface as the global climate warms, but such projections have typically relied on coarse‐resolution (100–250 km) general circulation models (GCMs). The ensemble of HighResMIP GCMs presents an opportunity to evaluate how a more finely resolved atmosphere and land‐surface might enhance the fidelity of the simulated contribution of large‐magnitude storms to total precipitation, particularly across topographically complex terrain. Here, the simulation of large‐storm dominance, that is, the number of wettest days to reach half of the total annual precipitation, is quantified across the western United States (WUS) using four GCMs within the HighResMIP ensemble and their coarse resolution counterparts. Historical GCM simulations (1950–2014) are evaluated against a baseline generated from station‐observed daily precipitation (4,803 GHCN‐D stations) and from three gridded, observationally based precipitation data sets that are coarsened to match the resolution of the GCMs. All coarse‐resolution simulations produce less large‐storm dominance than in observations across the WUS. For two of the four GCMs, bias in the median large‐storm dominance is reduced in the HighResMIP simulation, decreasing by as much as 62% in the intermountain west region. However, the other GCMs show little change or even an increase (+28%) in bias of median large‐storm dominance across multiple sub‐regions. The spread in differences with resolution amongst GCMs suggests that, in addition to resolution, model structure and parameterization of precipitation generating processes also contribute to bias in simulated large‐storm dominance.

Environmental Sciences & Ecology↗
DOE OSTI2022

Stochastic Model for High Temperature Oxidation of Cr–Ni Austenitic Steels Assisted by Spallation

Abstract Cr–Ni austenitic steels offer significant high temperature corrosion protection by forming a surface oxide layer. However, above critical service conditions (temperature, atmosphere, thermal cycling), oxidized surface can experience intensive degradation because of scale spallation, which could be detrimental to the in-service life. To predict the effect of scale spallation on oxidation kinetics, a simulation was implemented using a stochastic model. The model considers topological parameters and intensity of spallation which can occur, while delivering a true oxidation constant. The experimental procedure identified the amount of formed spalled scale and topology of spallation based on the use of element mapping of the surface. This information was used to determine a true kinetic constant for a corresponding spallation intensity in oxidized Cr–Ni austenitic steel. To illustrate the capability of the stochastic model, a parametric analysis was performed. The model verified how the spallation parameters could change the oxidation processes from parabolic growth of an adhered oxide layer without spallation to a mixed linear-parabolic, or with a constant thickness of residual scale at high spallation intensity. The spallation model will be used in a separate article to characterize high temperature surface degradation of several Cr–Ni austenitic steels during harsh oxidation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

NSA Site Science: Use of ARM Observations from Northern Alaska to Evaluate and Improve Prediction Capabilities

The Arctic is warming at a rate nearly double that of the rest of the planet, leading to profound changes in atmospheric, oceanic, and ice processes. The U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility has played a significant role in Arctic research, operating observatories in Alaska's North Slope for over 25 years. These observatories provide a rich, and wide-reaching dataset that offers insight into atmospheric processes in northern Alaska. This report details the results of a nine-year research project (2015–2024) supported by the DOE Atmospheric Systems Research (ASR) program, that leverages data from ARM’s deployment of observing facilities at Utqiaġvik (known as the North Slope of Alaska, or NSA, site) and Oliktok Point, Alaska. The project was conducted in two phases: - Phase 1 (2015–2019): Focused on understanding key atmospheric processes at Oliktok Point, including cloud formation, high-latitude precipitation, aerosol-cloud interactions, and cloud properties. - Phase 2 (2019–2024): Extended the research to the broader North Slope region, using data from both Oliktok Point and Utqiaġvik. Topics explored included surface energy budgets, atmospheric stability, ice nucleation processes, and microphysics in Arctic clouds. The project resulted in numerous research products, including 50 peer-reviewed publications and dissertations, 169 presentations, and 10 data products. These products cover a variety of topics, including: - Cloud Macro- and Microphysical Properties: Arctic clouds play a crucial role in energy transfer, and accurate representation in models is critical. The study explored cloud transitions, ice crystal shapes, and dual-wavelength radar data to understand ice crystal habits and size distributions. - Aerosol Properties and Processes: The team examined aerosol sources in the Arctic, including industrial emissions and natural sources. Observations showed significant spatial gradients in aerosol concentrations due to human activities and wildfire smoke. The influence of aerosols on cloud formation and the surface energy budget was also assessed. - Aerosol-Cloud Interactions: Research revealed that aerosols might suppress cloud ice production, affecting cloud radiative forcing and precipitation. The impact of local industrial emissions on cloud properties was also investigated. - Contextualizing the North Slope of Alaska in the context of the broader Arctic: To understand broader trends, the project evaluated large-scale circulation patterns and the influence of weather systems on the Arctic. Studies indicated that large-scale processes play a significant role in temperature patterns and the timing of snowmelt. - Advancing ARM Observational and Modeling Capabilities: The project developed new radar data products and advanced measurement techniques, including clutter mitigation and drizzle detection. Uncrewed aerial systems (UAS) and tethered balloon systems (TBS) were deployed to gather detailed atmospheric data. Additionally, the project supported 10 early career scientists, providing training and mentorship to undergraduate interns, graduate students, postdoctoral researchers, and early career researchers. These efforts contributed to the advancement of ARM research capabilities and fostered a new generation of scientists skilled in Arctic atmospheric research. Ultimately, this ASR-supported project has provided valuable insights into Arctic atmospheric processes and their broader climate implications. Recommendations for future work include continuing support for long-term observing at Arctic locations to foster additional research, further exploration of aerosol-cloud interactions and the potential impacts of enhanced industrialization of the Arctic, and expanded use of uncrewed systems to gather data in this remote and harsh environment. Additionally, the data products developed by this work, and the data products developed through the ARM infrastructure, leave a treasure-trove of additional information that should be explored for many years to come to gain additional insight into physical processes in the Arctic atmosphere that drive the rapid changes occurring in at high latitudes and their global impact.

58 GEOSCIENCES↗
DOE OSTI2021

University of Colorado and Black Swift Technologies RPAS-based measurements of the lower atmosphere during LAPSE-RATE

Abstract. Between 14 and 20 July 2018, small remotely piloted aircraft systems (RPASs) were deployed to the San Luis Valley of Colorado (USA) together with a variety of surface-based remote and in situ sensors as well as radiosonde systems as part of the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE). The observations from LAPSE-RATE were aimed at improving our understanding of boundary layer structure, cloud and aerosol properties, and surface–atmosphere exchange and provide detailed information to support model evaluation and improvement work. The current paper describes the observations obtained using four different types of RPASs deployed by the University of Colorado Boulder and Black Swift Technologies. These included the DataHawk2, the Talon and the TTwistor (University of Colorado), and the S1 (Black Swift Technologies). Together, these aircraft collected over 30 h of data throughout the northern half of the San Luis Valley, sampling altitudes between the surface and 914 m a.g.l. Data from these platforms are publicly available through the Zenodo archive and are co-located with other LAPSE-RATE data as part of the Zenodo LAPSE-RATE community (https://zenodo.org/communities/lapse-rate/, last access: 27 May 2021). The primary DOIs for these datasets are https://doi.org/10.5281/zenodo.3891620 (DataHawk2, de Boer et al., 2020a, e), https://doi.org/10.5281/zenodo.4096451 (Talon, de Boer et al., 2020d), https://doi.org/10.5281/zenodo.4110626 (TTwistor, de Boer et al., 2020b), and https://doi.org/10.5281/zenodo.3861831 (S1, Elston and Stachura, 2020).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Pressure Transducer Measurement Variability in Deep Wells Screened Across the Water Table

Abstract Automated water level measurements collected using vented pressure transducers in deep wells screened across the water table may exhibit a greater response to barometric pressure changes than the true water level. The cause was hypothesized to be disequilibrium in barometric pressure between the wellbores and land surface due to air exchange with the deep vadose zone. In this study, vented and nonvented pressure transducers were installed and operated simultaneously in two deep wells screened across the water table. A vent tube open to the atmosphere at land surface allowed for barometric compensation of the vented transducers. Two nonvented transducers were installed in each well, one submerged in the water and one above the water surface. The difference in readings allowed for barometric compensation. Manual measurements were also collected. It was confirmed that measurements from the vented transducers exhibited greater variability in response to barometric pressure changes than the nonvented transducers and manual measurements. Comparison of the downhole barometric pressure measurements to values from a nearby meteorology station showed the response in the wells to changes in barometric pressure was time‐lagged and attenuated. Thus, the reference pressure from land surface supplied to the vented transducers was not representative of the air pressure within the wells. This caused fluctuations of the transducer readings in response to barometric pressure changes to be greater than the true water level change. This issue can be resolved by the use of nonvented pressure transducers.

McDonald, John P.↗
DOE OSTI2021

Effect of Humidity on Charge Decay in Varying Atmospheric Gases

Characterization of surface charge decay for varying materials – PTFE, Acrylic, and metals – of 100 mm diameter, radially symmetric bodies was performed. In a gaseous environment, charged materials will lose their charge depending on gas type, humidity, and other environmental factors. Eventually, the charge will be (1) neutralized by ions or electrons in the surrounding gas (2) distributed by conduction on the material surface or in its volume (3) or neutralized by ions or electrons in the material surface and volume. Metal and polymer samples are charged to a known potential in the tens of kilovolts range, and the decaying voltage is monitored with a non-contact probe until the potential reaches near zero. Humidity has an impact on both groups of materials. As a baseline, in air at 42% relative humidity in an enclosed chamber, a positive charge equivalent to 20 kV potential on a metal sphere of 50 mm radius decays to near zero within 25 minutes. In contrast, the same sphere negatively charged decays slightly faster within 20 minutes. Overall, the humidity has a significant effect on the charge decay time, such that an order of magnitude longer times (hours) are observed at low humidity. Effects of chamber size were also studied, drawing a connection to charge carriers in the air available for decay, with chambers of volume(s): 3.3 105 cm3 and 1.4 106 cm3, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessing the Robustness of Arctic Sea Ice Bi‐Stability in the Presence of Atmospheric Feedbacks

Abstract Arctic sea‐ice loss is influenced by multiple positive feedbacks, sparking concerns of accelerated loss in the coming years or even a tipping point, where a sea‐ice equilibrium disappears at a given CO 2 value and sea ice rapidly evolves to a new steady state. Such a tipping point would imply a bi‐stability of the Arctic climate—where multiple steady‐state Arctic climates are possible at the same CO 2 value. Previous works have sought to establish the existence of bi‐stability using a range of models, from zero‐dimensional sea ice thermodynamic models to fully coupled global climate models, with conflicting results. Here, we present a new model of the Arctic that includes both sea‐ice thermodynamics and key atmospheric feedbacks in a simple framework. We exploit the model's simplicity to identify physical mechanisms that control the timing and extent of sea‐ice bi‐stability, and the abruptness of ice loss. We show that longwave radiation feedbacks can have a strong influence on Arctic surface climate from atmospheric temperature increases alone, even without major contributions from clear‐sky moisture or convective clouds suggested previously. While winter sea‐ice bi‐stability is robust to changes in uncertain model parameters in this study, summer sea ice is more sensitive. Finally, our model indicates that positive feedbacks may modulate the CO 2 threshold of sea‐ice loss and the width of bi‐stability much more strongly than the abruptness of loss. These results lead to a comprehensive understanding of the conditions that favor Arctic sea‐ice bi‐stability, particularly the role of atmospheric feedbacks, in both future and past climates.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

Planetary Boundary-Layer Height (PBLHT) Value-Added Product: Remote-Sensing Retrievals

The planetary boundary layer (PBL) is fundamental to numerous atmospheric processes, including aerosol mixing and transport, cloud evolution, and precipitation formation. A critical parameter in these studies is the PBL height (PBLHT). This vertical depth is essential for characterizing PBL structures in numerical simulations and serves as a primary metric for estimating flux exchanges between the Earth’s surface and the atmosphere. Radiosonde (SONDE) observations provide high-vertical-resolution measurements of temperature and moisture profiles and are widely used to estimate PBLHT (Liu and Liang 2010, Seidel et al. 2010). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s PBLHT value-added product (VAP) for radiosonde measurements, known as PBLHTSONDE, applies three commonly used methods—the Heffter (1980) method, the Liu and Liang (2010) method, and the bulk Richardson number approach (Seibert et al. 2000)—to derive PBLHT. The PBLHTSONDE VAP operates routinely at ARM observatories and mobile facilities, with data available from the ARM Data Center shortly after sounding observations are collected (Sivaraman et al. 2013). However, radiosonde observations are limited by their low temporal resolution. Most stations launch soundings only twice daily, which constrains the ability to investigate and characterize the temporal evolution of the PBL using radiosonde data alone. The use of continuous remote-sensing observations provides high temporal resolution of PBLHT estimates. These observations include aerosol lidars (Dang et al. 2019, Su et al. 2020), Doppler lidar (DL; Tucker et al. 2009, Krishnamurthy et al. 2021), and water vapor and/or temperature lidars and radiometers (Turner et al. 2014). These observations provide valuable data on the PBL’s thermodynamic properties (e.g., water vapor and/or temperature lidars and radiometers), dynamic properties (e.g., DL), and distribution of tracer substances (e.g., aerosol lidars), all of which can be used to estimate PBLHT. ARM developed PBLHT estimates from the micropulse lidar (MPL; PBLHTMPL), Doppler lidar (PBLHTDL), and combined Raman lidar (RL)/atmospheric emitted radiance interferometer (AERI) thermodynamic profiles (PBLHTTHERMO). Each estimate captures different physical characteristics of the boundary layer—aerosol tracers, vertical velocity turbulence, and thermodynamic structure—and exhibits distinct strengths and limitations depending on the PBL regime and time of day. In addition, the ARM ceilometer (CEIL) provides three potential PBLHT candidates derived from the vendor's built-in algorithm. Building on these individual retrievals, ARM developed the PBLHTBEML VAP, which combines the four remote-sensing-based estimates with ancillary meteorological variables using the machine learning approach of Zhang et al. (2025) to produce a best-estimate PBLHT at 10-minute resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Systematically tracking the hourly progression of large wildfires using GOES satellite observations

Abstract. In the western United States, prolonged drought, a warming climate, and historical fuel buildup have contributed to larger and more intense wildfires as well as to longer fire seasons. As these costly wildfires become more common, new tools and methods are essential for improving our understanding of the evolution of fires and how extreme weather conditions, including heat waves, windstorms, droughts, and varying levels of active-fire suppression, influence fire spread. Here, we develop the Geostationary Operational Environmental Satellites (GOES)-Observed Fire Event Representation (GOFER) algorithm to derive the hourly fire progression of large wildfires and create a product of hourly fire perimeters, active-fire lines, and fire spread rates. Using GOES-East and GOES-West geostationary satellite detections of active fires, we test the GOFER algorithm on 28 large wildfires in California from 2019 to 2021. The GOFER algorithm includes parameter optimizations for defining the burned-to-unburned boundary and correcting for the parallax effect from elevated terrain. We evaluate GOFER perimeters using 12 h data from the Visible Infrared Imaging Radiometer Suite (VIIRS)-derived Fire Event Data Suite (FEDS) and final fire perimeters from the California's Fire and Resource Assessment Program (FRAP). Although the GOES imagery used to derive GOFER has a coarser resolution (2 km at the Equator), the final fire perimeters from GOFER correspond reasonably well to those obtained from FRAP, with a mean Intersection-over-Union (IoU) of 0.77, in comparison to 0.83 between FEDS and FRAP; the IoU indicates the area of overlap over the area of the union relative to the reference perimeters, in which 0 is no agreement and 1 is perfect agreement. GOFER fills a key temporal gap present in other fire tracking products that rely on low-Earth-orbit imagery, where perimeters are available at intervals of 12 h or longer or at ad hoc intervals from aircraft overflights. This is particularly relevant when a fire spreads rapidly, such as at maximum hourly spread rates of over 5 km h−1. Our GOFER algorithm for deriving the hourly fire progression using GOES can be applied to large wildfires across North and South America and reveals considerable variability in the rates of fire spread on diurnal timescales. The resulting GOFER product has a broad set of potential applications, including the development of predictive models for fire spread and the improvement of atmospheric transport models for surface smoke estimates. The resulting GOFER product has a broad set of potential applications, including the development of predictive models for fire spread and the improvement of atmospheric transport models for surface smoke estimates (https://doi.org/10.5281/zenodo.8327264, Liu et al., 2023).

Geology↗
DOE OSTI2023

Atmospheric Turbulent Intermittency Over the Arctic Sea‐Ice Surface During the MOSAiC Expedition

Turbulent motions in the Arctic stable boundary layer are characterized by intermittency, but they are rarely investigated due to limited observations, in particular over the sea‐ice surface. In the present study, we explore the characteristics of turbulent intermittency over the Arctic sea‐ice surface using data collected during the Multidisciplinary drifting Observation for the Study of Arctic Climate expedition from October 2019 to September 2020. We first develop a new algorithm, which performs well in identifying the spectral gap over the Arctic sea‐ice surface. Then the characteristics of intermittency are investigated. It is found that the strength of intermittency increases under the conditions of light surface wind speed, small surface wind speed gradient, and strong surface air temperature gradient. The momentum flux, sensible heat flux, and latent heat flux calculated by raw eddy‐covariance fluctuations are overestimated by 3%, 10%, and 24%, respectively, because submesoscale motions are included. Furthermore, the characteristics of the atmospheric boundary layer structure under various intermittency conditions reveal that strong low‐level jets are favorable to surface turbulent motions that result in weak intermittency, while strong temperature inversions above the surface layer suppress surface turbulent motions and lead to strong intermittency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Accumulation of microplastics in soil after long-term application of biosolids and atmospheric deposition

Land-applied biosolids can be a considerable source of microplastics in soils. Previous studies reported microplastics accumulation in soils from biosolid application, however, little is known about the contribution of atmospherically deposited microplastics to agricultural soils. In this study, we quantified and characterized microplastics in soils that have been amended with biosolids over the past 23 years. We also collected atmospheric deposition samples to determine the amount and type of plastics added to soils through atmospheric input over a period of about 2 years. Soil samples were taken from a replicated field trial where biosolids have been applied at rates of 0, 4.8, 6.9, and 9.0 t/ha every second crop. The biosolids were anaerobically digested and dewatered, and were applied by spreading onto the soil surface. Soil and atmospheric samples were extracted for microplastics by Fenton's reaction to remove organic matter followed by flotation in a zinc chloride solution to separate plastic from soil particles. Samples were analyzed for microplastics by optical microscopy and Laser Direct Infrared Imaging Analysis (LDIR). The mean number of microplastics identified from biosolids samples was 12,000 particles/kg dry biosolids. The long-term applications of biosolids to the soil led to mean plastics concentrations of 383, 500, and 361 particles/kg dry soil in the 0–10 cm depth for low, medium, and high biosolids application rates, respectively. These plastic concentrations were not significantly different from each other, but significantly higher than those found in non biosolids-amended soil (117 particles/kg dry soil). The dominant plastic types by number found in biosolids were polyurethane, followed by polyethylene, and polyamide. The most abundant plastics in soil samples were polyurethane, polyethylene terephthalate, polyamide, and polyethylene. Atmospheric deposition contributed to 15 particles/kg dry soil per year and was mainly composed of polyamide fibers. In conclusion, this study shows that long-term application of biosolids led to an accumulation of microplastics in soil, but that atmospheric deposition also contributes a considerable input of microplastics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

High-Resolution ESM Projections for Energy Applications Over the CONUS

Assessing energy resources under future scenarios requires high-resolution meteorological information that is physically consistent and suitable for regional-scale analysis. While Earth system model (ESM) projections provide valuable large-scale information, their coarse resolution and systematic biases limit direct applicability for energy system modeling and planning. In this study, we develop a high-resolution dynamical downscaling framework based on the Weather Research and Forecasting (WRF) model to translate global-scale ESM data into energy-relevant regional projections over the contiguous United States (CONUS). The framework identifies an optimized WRF configuration through numerical experiments and evaluates raw and bias-corrected ESM initial and boundary conditions, with soil moisture (SM) and soil temperature (ST) bias correction implemented as an integral part of the bias-corrected ESM forcing to improve land-atmosphere coupling prior to WRF dynamical downscaling. Using an optimized WRF configuration at 4-km resolution, we show that raw ESM forcing introduces systematic dry and cold soil biases that propagate into pronounced warm biases in near-surface air temperature and positive biases in solar irradiance, particularly during summer. Applying bias-corrected atmospheric forcing together with bias-corrected SM and ST substantially reduces these downstream biases and improves the surface energy balance and near-surface atmospheric fields. These results demonstrate that bias-aware treatment of initial conditions is critical for producing high-resolution downscaled projections suitable for energy system modeling and planning applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗