Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high-latitudes”

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

Flux Cube Reconstruction from Slitless Spectroscopy

Slitless spectroscopy enables efficient, large-area surveys without target preselection, yet it faces challenges from source blending, higher noise, and lost spatial–spectral information. We present an advanced, nonparametric, data-driven algorithm that leverages multiple dispersion angles to reconstruct three-dimensional flux distributions, providing low-resolution integral field unit capabilities from slitless data. By treating each pixel as an independent element, our method naturally handles source confusion without requiring prior assumptions regarding redshifts, templates, or model libraries. We validate the algorithm using simulated Roman Space Telescope wide-field slitless spectroscopy images that are equivalent to what is expected from the High-Latitude Time-Domain Survey. First, we demonstrate that a host-galaxy model reconstructed from multiple dispersion angles can be used to accurately subtract host light from a transient, recovering a Type Ia supernova spectrum with minimal bias. Second, we showcase a high-fidelity flux-cube reconstruction of a complex galaxy, successfully measuring the redshift and recovering continuum, emission, and absorption features. This approach highlights the potential of multi-dispersion-angle slitless data to provide spatially resolved spectral information in a nonparametric way, which is traditionally accessible only with integral field spectroscopy, opening a new window into large, unbiased, and spatially resolved studies of galaxy evolution.

Griggio, M. [Space Telescope Science Institute, Ba↗
DOE OSTI2023

Stellar-reddening-based Extinction Maps for Cosmological Applications

Cosmological surveys must correct their observations for the reddening of extragalactic objects by Galactic dust. Existing dust maps, however, have been found to have spatial correlations with the large-scale structure of the Universe. Errors in extinction maps can propagate systematic biases into samples of dereddened extragalactic objects and into cosmological measurements such as correlation functions between foreground lenses and background objects and the primordial non-Gaussianity parameter f NL . Emission-based maps are contaminated by the cosmic infrared background, while maps inferred from stellar reddenings suffer from imperfect removal of quasars and galaxies from stellar catalogs. Thus, stellar-reddening-based maps using catalogs without extragalactic objects offer a promising path to making dust maps with minimal correlations with large-scale structure. We present two high-latitude integrated extinction maps based on stellar reddenings, with a point-spread functions of FWHMs 6'. 1 and 15'. We employ a strict selection of catalog objects to filter out galaxies and quasars and measure the spatial correlation of our extinction maps with extragalactic structure. Our galactic extinction maps have reduced spatial correlation with large-scale structure relative to most existing stellar-reddening-based and emission-based extinction maps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Spectral Properties of Globally Distributed ENA Fluxes across Diverse Regions of the Heliosphere

This study analyzes energetic neutral atom (ENA) spectral properties across distinct regions of globally distributed flux (GDF) sky maps, using Interstellar Boundary Explorer data from a full solar cycle, corrected for time dispersion. By time-shifting the data to the heliosheath using GDF source distances from D. B. Reisenfeld et al., we achieve a more accurate representation of heliosheath GDF energy spectra. We quantify ENA spectral characteristics, heliosheath line-of-sight-integrated proton pressure, and heliosheath proton temperature, comparing these to solar wind properties at 1 au and interplanetary scintillation-derived solar wind data. Our findings show that the spectral index is generally anticorrelated with heliosheath proton temperature and pressure, except in the central tail, where a partial positive correlation is observed. The lowest spectral index values occur when high-latitude heliosheath regions are dominated by fast solar wind from polar coronal holes. The south pole exhibits the flattest energy spectra due to plasma heating from both fast solar wind and a late-2014 pressure pulse. The central tail shows shorter variability (5–6 yr) for spectral index and heliosheath proton temperature, while proton pressure follows the 11 yr solar cycle. Most spectral shapes exhibit a “knee” distribution, peaking during solar maximum, with an “ankle” shape observed only at the south pole during solar cycle transitions. Asymmetry in proton pressure in the lobes is driven by the draping effect of the local interstellar magnetic field. This study provides insights into the energetic properties of GDF across the heliosphere, enhancing our understanding of the heliospheric environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Impact of Cosmic Filaments on Galaxy Morphological Evolution and Predictions of Early Cosmic Web Structure for Roman

We leverage the IllustrisTNG cosmological simulations to test how the large-scale cosmic web shapes galaxy morphology and to forecast the early cosmic web structure that the Nancy Grace Roman Space Telescope will reveal. In the hydrodynamic TNG50 and N-body TNG50-Dark runs, we reconstruct the cosmic web at redshifts z = 0, 0.5, 1, 2, 3, and 4 with the Monte Carlo Physarum Machine density estimator and the DisPerSE structure identification framework. We confirm that dark matter halos start out predominantly prolate (elongated), and their shapes are aligned with their nearest filaments; prolate galaxies retain strong shape alignment with their outer halos to later times. At z ≥ 1, the fraction of prolate (spheroidal) halos increases (decreases) toward lower stellar mass, higher redshift, and lower filament density. At z < 1, more spheroidal (oblate) stellar structures preferentially reside in higher-density (lower-density) filaments. We also find that higher-density filaments favor extended rotationally supported disks, whereas lower-density filaments more often host smaller dispersion-supported systems. Then, generating mock galaxy samples from TNG100 and TNG50, we predict the early cosmic web accessible to Roman. We find that the spectroscopic emission-line depth planned for the High-Latitude Wide-Area Survey (HLWAS) yields a highly incomplete galaxy sample that does not accurately trace the z = 1 cosmic web. A survey ≥2.5× deeper over a few square degrees would enable a proper reconstruction and reveal qualitatively correct filament–galaxy morphology relationships. Nevertheless, the planned HLWAS Deep field should still identify most galaxy overdensities; targeted deeper spectroscopy of these regions would efficiently and adequately map the early filamentary structure.

Hasan, Farhanul [Space Telescope Science Institute↗
DOE OSTI2026

Global Characteristics of the IBEX Ribbon and Its Temporal Variation: Implications for Relationships with Solar Cycle Variability and the Interstellar Magnetic Field

The Interstellar Boundary Explorer (IBEX) has revealed the existence of the energetic neutral atom (ENA) ribbon, a narrow and nearly circular band of enhanced flux that enables probing plasma and magnetic field conditions outside the heliopause (HP). In this study, we investigate the global structure and temporal variation of the IBEX ribbon using 14 yr of IBEX-Hi ENA observations. We derive geometric parameters of the ribbon: center location, angular radius, eccentricity, and major axis orientation. We find that the ribbon exhibits systematic variations tied to the solar cycle. The ribbon center shows a solar cycle–modulated latitudinal motion at low energies. The angular radius increases with ENA energy and displays a pronounced enlargement between 2013 and 2018, coinciding with changes in high-latitude fast solar wind from polar coronal holes. The eccentricity remains consistently around 0.3 with some temporal evolution, while the major axis rotates gradually over time, becoming increasingly perpendicular to the ecliptic meridional plane. A notable hemispheric asymmetry has emerged over time: the ribbon in the northern hemisphere remains remarkably stable in position, whereas the southern ribbon exhibits temporal migration. By combining these results with solar wind observations, global heliosphere, and ribbon modeling, we suggest that this asymmetry is explained by a coupled effect of magnetic field draping and the latitude-dependent timing of fast solar wind from polar coronal holes. These results establish the most detailed long-term geometric characterization of the IBEX ribbon to date and provide observational support for models of secondary ENA production and magnetic field draping outside the HP.

Noh, Sung Jun [New Mexico Consortium, Los Alamos, ↗
DOE OSTI2023

Physicochemical characterization and source apportionment of Arctic ice-nucleating particles observed in Ny-Ålesund in autumn 2019

Abstract. Ice-nucleating particles (INPs) initiate primary ice formation in Arctic mixed-phase clouds (MPCs), altering cloud radiative properties and modulating precipitation. For atmospheric INPs, the complexity of their spatiotemporal variations, heterogeneous sources, and evolution via intricate atmospheric interactions challenge the understanding of their impact on microphysical processes in Arctic MPCs and induce an uncertain representation in climate models. In this work, we performed a comprehensive analysis of atmospheric aerosols at the Arctic coastal site in Ny-Ålesund (Svalbard, Norway) from October to November 2019, including their ice nucleation ability, physicochemical properties, and potential sources. Overall, INP concentrations (NINP) during the observation season were approximately up to 3 orders of magnitude lower compared to the global average, with several samples showing degradation of NINP after heat treatment, implying the presence of proteinaceous INPs. Particle fluorescence was substantially associated with INP concentrations at warmer ice nucleation temperatures, indicating that in the far-reaching Arctic, aerosols of biogenic origin throughout the snow- and ice-free season may serve as important INP sources. In addition, case studies revealed the links between elevated NINP and heat lability, fluorescence, high wind speeds originating from the ocean, augmented concentration of coarse-mode particles, and abundant organics. Backward trajectory analysis demonstrated a potential connection between high-latitude dust sources and high INP concentrations, while prolonged air mass history over the ice pack was identified for most scant INP cases. The combination of the above analyses demonstrates that the abundance, physicochemical properties, and potential sources of INPs in the Arctic are highly variable despite its remote location.

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

Natural marine cloud brightening in the Southern Ocean

Abstract. The number of cloud droplets per unit volume (Nd) is a fundamentally important property of marine boundary layer (MBL) liquid clouds that, at constant liquid water path, exerts considerable controls on albedo. Past work has shown that regional Nd has a direct correlation to marine primary productivity (PP) because of the role of seasonally varying, biogenically derived precursor gases in modulating secondary aerosol properties. These linkages are thought to be observable over the high-latitude oceans, where strong seasonal variability in aerosol and meteorology covary in mostly pristine environments. Here, we examine Nd variability derived from 5 years of MODIS Level 2-derived cloud properties in a broad region of the summer eastern Southern Ocean and adjacent marginal seas. We demonstrate latitudinal, longitudinal and temporal gradients in Nd that are strongly correlated with the passage of air masses over high-PP waters that are mostly concentrated along the Antarctic Shelf poleward of 60∘ S. We find that the albedo of MBL clouds in the latitudes south of 60∘ S is significantly higher than similar liquid water path (LWP) clouds north of this latitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Surface-based observations of cold-air outbreak clouds during the COMBLE field campaign

Abstract. Cold-air outbreaks (CAOs) are characterized by extreme air–sea energy exchanges and low-level convective clouds over large areas in the high-latitude oceans. As such, CAOs are an important component of the Earth's climate system. The CAOs in the Marine Boundary Layer Experiment (COMBLE) deployment of the US Department of Energy Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) provided the first comprehensive view of CAOs using a suite of ground-based observations at the northern coast of Norway. Here, cloud and precipitation observations from 13 CAO cases during COMBLE are analyzed. A vertical air motion retrieval technique is applied to the Ka-band ARM Zenith-pointing Radar (KAZR) observations. The CAO cumulus clouds are characterized by strong updrafts with magnitudes between 2–8 m s−1, vertical extents of 1–3 km, and horizontal scales of 0.25–3 km. A strong relationship between our vertical air velocity retrievals and liquid water path (LWP) measurements is found. The LWP measurements exceed 1 kg m−2 in strong updraft areas, and the vertical extent of the updraft correlates well with the LWP values. The CAO cumulus clouds exhibit eddy dissipation rate values between 10−3 and 10−2 m2 s−3 in the lowest 10 km of the atmosphere, and using a radar Doppler spectra technique, evidence of secondary ice production is found during one of the cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Validation of the Cloud_CCI (Cloud Climate Change Initiative) cloud products in the Arctic

The role of clouds in the Arctic radiation budget is not well understood. Ground-based and airborne measurements provide valuable data to test and improve our understanding. However, the ground-based measurements are intrinsically sparse, and the airborne observations are snapshots in time and space. Passive remote sensing measurements from satellite sensors offer high spatial coverage and an evolving time series, having lengths potentially of decades. However, detecting clouds by passive satellite remote sensing sensors is challenging over the Arctic because of the brightness of snow and ice in the ultraviolet and visible spectral regions and because of the small brightness temperature contrast to the surface. Consequently, the quality of the resulting cloud data products needs to be assessed quantitatively. In this study, we validate the cloud data products retrieved from the Advanced Very High Resolution Radiometer (AVHRR) post meridiem (PM) data from the polar-orbiting NOAA-19 satellite and compare them with those derived from the ground-based instruments during the sunlit months. The AVHRR cloud data products by the European Space Agency (ESA) Cloud Climate Change Initiative (Cloud_CCI) project uses the observations in the visible and IR bands to determine cloud properties. The ground-based measurements from four high-latitude sites have been selected for this investigation: Hyytiälä (61.84°N, 24.29°E), North Slope of Alaska (NSA; 71.32°N, 156.61°W), Ny-Ålesund (Ny-Å; 78.92°N, 11.93°E), and Summit (72.59°N, 38.42°W). The liquid water path (LWP) ground-based data are retrieved from microwave radiometers, while the cloud top height (CTH) has been determined from the integrated lidar–radar measurements. The quality of the satellite products, cloud mask and cloud optical depth (COD), has been assessed using data from NSA, whereas LWP and CTH have been investigated over Hyytiälä, NSA, Ny-Å, and Summit.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Using atmospheric observations to quantify annual biogenic carbon dioxide fluxes on the Alaska North Slope

The continued warming of the Arctic could release vast stores of carbon into the atmosphere from high-latitude ecosystems, especially from thawing permafrost. Increasing uptake of carbon dioxide (CO 2 ) by vegetation during longer growing seasons may partially offset such release of carbon. However, evidence of significant net annual release of carbon from site-level observations and model simulations across tundra ecosystems has been inconclusive. To address this knowledge gap, we combined top-down observations of atmospheric CO 2 concentration enhancements from aircraft and a tall tower, which integrate ecosystem exchange over large regions, with bottom-up observed CO 2 fluxes from tundra environments and found that the Alaska North Slope is not a consistent net source nor net sink of CO 2 to the atmosphere (ranging from ₋6 to +6 Tg C yr -1 for 2012–2017). Our analysis suggests that significant biogenic CO 2 fluxes from unfrozen terrestrial soils, and likely inland waters, during the early cold season (September–December) are major factors in determining the net annual carbon balance of the North Slope, implying strong sensitivity to the rapidly warming freeze-up period. At the regional level, we find no evidence of the previously reported large late-cold-season (January–April) CO 2 emissions to the atmosphere during the study period. Despite the importance of the cold-season CO 2 emissions to the annual total, the interannual variability in the net CO 2 flux is driven by the variability in growing season fluxes. During the growing season, the regional net CO 2 flux is also highly sensitive to the distribution of tundra vegetation types throughout the North Slope. This study shows that quantification and characterization of year-round CO 2 fluxes from the heterogeneous terrestrial and aquatic ecosystems in the Arctic using both site-level and atmospheric observations are important to accurately project the Earth system response to future warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hemispherically symmetric strategies for stratospheric aerosol injection

Stratospheric aerosol injection (SAI) comes with a wide range of possible design choices, such as the location and timing of the injection. Different stratospheric aerosol injection strategies can yield different climate responses; therefore, understanding the range of possible climate outcomes is crucial to making informed future decisions on SAI, along with the consideration of other factors. Yet, to date, there has been no systematic exploration of a broad range of SAI strategies. This limits the ability to determine which effects are robust across different strategies and which depend on specific injection choices. This study systematically explores how the choice of SAI strategy affects climate responses in one climate model. Here, we introduce four hemispherically symmetric injection strategies, all of which are designed to maintain the same global mean surface temperature: an annual injection at the Equator (EQ), an annual injection of equal amounts of SO 2 at 15° N and 15° S (15N+15S), an annual injection of equal amounts of SO 2 at 30° N and 30° S (30N+30S), and a polar injection strategy that injects equal amounts of SO 2 at 60° N and 60° S only during spring in each hemisphere (60N+60S). We compare these four hemispherically symmetric SAI strategies with a more complex injection strategy that injects different quantities of SO 2 at 30° N, 15° N, 15° S, and 30° S in order to maintain not only the global mean surface temperature but also its large-scale horizontal gradients. All five strategies are simulated using version 2 of the Community Earth System Model with the middle atmosphere version of the Whole Atmosphere Community Climate model, version 6, as the atmospheric component, CESM2(WACCM6-MA), with the global warming scenario, Shared Socioeconomic Pathway (SSP)2-4.5. We find that the choice of SAI strategy affects the spatial distribution of aerosol optical depths, injection efficiency, and various surface climate responses. In addition, injecting in the subtropics produces more global cooling per unit injection, with the EQ and the 60N+60S cases requiring, respectively, 59 % and 50 % more injection than the 30N+30S case to meet the same global mean temperature target. Injecting at higher latitudes results in larger Equator-to-pole temperature gradients. While all five strategies restore Arctic September sea ice, the high-latitude injection strategy is more effective due to the SAI-induced cooling occurring preferentially at higher latitudes. These results suggest trade-offs wherein different strategies appear better or worse, depending on which metrics are deemed important.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Machine-learning-based estimates of global natural vegetated wetland methane emissions (2000–2025)

Wetlands are the largest natural source of atmospheric methane (CH 4 ), yet comprehensive global budgets are typically delayed by years, preventing a timely understanding of CH 4 sources, sinks, and trends. To reduce this delay, we present a model emulator-driven framework and accompanying workflow that enable timely, continuous emission updates using a machine-learning emulator to reconstruct spatially explicit monthly emission fields at 1° × 1° resolution. We apply this framework to a global dataset of natural vegetated wetland CH 4 emissions to extend the most recent Global Methane Budget (GMB; Saunois et al., 2025) record that covers the 2000–2020 emissions through 2025. In the test data (∼ 30 % of the total dataset), the emulator achieved a global R 2 of 0.65 ± 0.003 (mean ± 95 % CI, hereafter) and an RMSE of 5.49 ± 0.12×10 -3 Tg CH 4 yr −1 . The emulator is trained on 35 GMB model estimates, including 22 process-based models and 13 atmospheric inversions, paired with 10 ensemble realizations of 11 gridded climate predictor variables from atmospheric reanalyses. Our results show that the global mean predicted wetland CH 4 emissions for 2021–2025 (157.8 ± 2.4 Tg CH 4 yr −1 ) are not significantly higher (∼ 0.05 Tg CH 4 yr −1 ) than the 2000–2020 baseline. However, this stability masks a significant hemispheric redistribution of emissions. We detect an increase in Northern Hemisphere (NH) emissions in 2021–2025, with mid- and high-latitudes increasing by 0.76 ± 0.07 and 0.35 ± 0.03 Tg CH 4 yr −1 , respectively, while the tropics and Southern Hemisphere (SH) extratropics show offsetting negative trends (−0.95 ± 0.19 and -0.11 ± 0.02 Tg CH 4 yr −1 , respectively). The predicted emissions are able to capture the low emissions in 2023 in South America linked to El Niño-related drought, as reported by recent studies (Ciais et al., 2026; Quinn et al., 2025). Furthermore, we identify a distinct seasonal amplification of global emission trends that peaks in late boreal summer. This new modeled dataset and operational framework bridge the gap between the latest updated budgets and low-latency monitoring, providing a scalable capacity to frequently update global emission estimates and critical early warnings of regional wetland feedback loops. The data are publicly available at https://doi.org/10.5281/zenodo.18870108 (Li et al., 2026).

Li, Mengze [National University of Singapore (Sing↗
DOE OSTI2026

Benchmarking soil moisture and its relationship to ecohydrologic variables in Earth System Models

Soil moisture (SM) is a key regulator of ecosystem biogeophysics, influencing plant water relations and land-atmosphere energy exchanges. We evaluate the representation of SM in 16 Earth System Models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) using the International Land Model Benchmarking (ILAMB) framework, focusing on surface (0–5, 0–10 cm) and rootzone (0–100 cm) depths, as well as key ecohydrological variables like gross primary productivity (GPP), leaf area index (LAI), and evapotranspiration (ET), and their coupling. Models are benchmarked against multiple observational and assimilated datasets to assess both state variables and cross-variable relationships. Surface SM is generally well represented (r > 0.87), while rootzone SM variability is systematically overestimated (normalized standard deviation > 1). ET shows strong agreement with observations (r > 0.9), whereas GPP and LAI exhibit larger inter-model spread. Skill in individual variables does not guarantee realistic SM–ecohydrology coupling, which varies strongly across models and depends on the reference dataset. Köppen-based regional analyses reveal strong regime dependence, with several models performing well in Tropical and Temperate regions but degrading in Continental (high-latitude) zones. Across both global and regional benchmarks, models cluster by land surface framework, indicating that structural choices in soil hydrology and soil–plant coupling exert a first-order control on performance. These results provide process-relevant benchmarks and suggest that improving the representation of vertical soil structure, rooting depth distributions, and soil–plant hydraulic coupling will be central to advancing soil moisture realism in next-generation Earth system models.

CMIP6↗
DOE OSTI2024

Wetting and drying trends in the land–atmosphere reservoir of large basins around the world

Abstract. Global change is altering hydrologic regimes worldwide, including large basins that play a central role in the sustainability of human societies and ecosystems. The basin water budget is a fundamental framework for understanding these basins' sensitivity and future dynamics under changing forcings. In this budget, studies often treat atmospheric processes as external to the basin and assume that atmosphere-related water storage changes are negligible in the long term. These assumptions are potentially misleading in large basins with strong land–atmosphere feedbacks, including terrestrial moisture recycling, which is critical for global water distribution. Here, we introduce the land–atmosphere reservoir (LAR) concept, which includes atmospheric processes as a critical component of the basin water budget and use it to study long-term changes in the water storage of some of the world's largest basins. Our results show significant LAR water storage trends over the last 4 decades, with a marked latitudinal contrast: while low-latitude basins have accumulated water, high-latitude basins have been drying. If they continue, these trends will disrupt the discharge regime and compromise the sustainability of these basins, resulting in widespread impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Climatology and surface impacts of atmospheric rivers on West Antarctica

Abstract. Atmospheric rivers (ARs) transport large amounts of moisture from the mid- to high-latitudes and they are a primary driver of the most extreme snowfall events, along with surface melting, in Antarctica. In this study, we characterize the climatology and surface impacts of ARs on West Antarctica, focusing on the Amundsen Sea Embayment and Marie Byrd Land. First, we develop a climatology of ARs in this region, using an Antarctic-specific AR detection tool combined with the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) and the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5) atmospheric reanalyses. We find that while ARs are infrequent (occurring 3 % of the time), they cause intense precipitation in short periods of time and account for 11 % of the annual surface accumulation. They are driven by the coupling of a blocking high over the Antarctic Peninsula with a low-pressure system known as the Amundsen Sea Low. Next, we use observations from automatic weather stations on Thwaites Eastern Ice Shelf with the firn model SNOWPACK and interferometric reflectometry (IR) to examine a case study of three ARs that made landfall in rapid succession from 2 to 8 February 2020, known as an AR family event. While accumulation dominates the surface impacts of the event on Thwaites Eastern Ice Shelf (> 100 kg m−2 or millimeters water equivalent), we find small amounts of surface melt as well (< 5 kg m−2). The results presented here enable us to quantify the past impacts of ARs on West Antarctica's surface mass balance (SMB) and characterize their interannual variability and trends, enabling a better assessment of future AR-driven changes in the SMB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Daily Arctic Lightning Strokes

In recent decades, lightning activity at high latitudes has increased. Tall thunderstorm clouds affect radiation balance directly, as well as indirectly through lightning-initiated fires and the resulting smoke. One can remotely sense lightning strokes over the globe through their VLF radio emission, and, with multiple receivers, it is possible to precisely locate lightning strokes. This technique makes it possible to continuously monitor arctic lightning--a capability not possible by other means. In this research effort, a World-Wide Lightning Location Network station, consisting of a VLF receiver, signal processing hardware, and analysis software, were installed at the North Slope of Alaska (NSA) facility and planned to operate for several years. This far north location is expected to improve the network's high-latitude detection efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗