SEARCH · Engineering Papers

Results for “variability study”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Controlled Shifts of X-ray Emission Lines Measured with Transition Edge Sensors at the Advanced Photon Source

The measurement of shifts in the energy of X-ray emission lines is important for understanding the electronic structure and physical properties of materials. In this study, we demonstrate a method using a synchrotron source to introduce controlled eV-scale shifts of a narrow line in between fixed-energy fluorescence lines. We use this to characterize the ability of a hard X-ray superconducting Transition Edge Sensor (TES) array to measure line shifts. Fixed fluorescence lines excited by higher harmonics of the monochromatic X-ray beam are used for online energy calibration, while elastic scattering from the primary harmonic acts as the variable energy emission line under study. We use this method to demonstrate the ability to track shifts in the energy of the elastic scattering line of magnitude smaller than the TES energy resolution, and find we are ultimately limited by our calibration procedure. The method can be applied over a wide X-ray energy range and provides a robust approach for the characterization of the ability of high-resolution detectors to detect X-ray emission line shifts, and the quantitative comparison of energy calibration procedures.

Guruswamy, Tejas (ORCID:0000000241656765)↗

DOE OSTI2022

XMM–Newton study of the Sculptor dwarf spheroidal galaxy

ABSTRACT We analysed observations with XMM–Newton in the field of the Sculptor dwarf spheroidal galaxy (dSph). The aim of the study was the classification of X-ray binaries and accreting white dwarfs belonging to the Sculptor dSph. Using different methods of X-ray timing and spectral analyses, together with an extensive multiwavelength study of the optical and infrared counterparts of the X-ray sources, we classified the sources detected with XMM–Newton in the field of the Sculptor dSph. The long-term variability of the sources has been studied over two XMM–Newton observations. None of the members of the Sculptor dSph show significant long-term variability over these two observations. We also searched for periodicity and pulsation using the Lomb–Scargle and Rayleigh Z$^{2}_{n}$ techniques. No signals of pulsation or periodicity have been found for the X-ray sources. The results show the presence of a noticeable number of background X-ray sources in the field of this galaxy. We classified 43 sources as active galactic nuclei, galaxies, and galaxy candidates. Three Galactic foreground stars have been identified in the field of the Sculptor dSph, and one of them is an M-dwarf candidate. Moreover, we classified four symbiotic-star candidates and three quiescent low-mass X-ray binary candidates in the Sculptor dSph. The luminosity of these X-ray sources is ∼1033−35 erg s−1.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Parameter-Varying Hydrodynamic Model of a Single Vane in a Variable-Geometry Oscillating Surge Wave Energy Converter

This paper presents a preliminary study on the parameter-varying hydrodynamic modeling of an individual vane of a variable-geometry oscillating surge wave energy converter (VGOSWEC). The WEC design incorporates controlled surfaces that can modify their orientation relative to the wave motion, reducing hydrodynamic pressure and related loads. This research focuses on characterizing the behavior of the oscillating WEC using a simplified model and three methods for achieving a continuous parameter-varying model: coarse hydrodynamic parameters, interpolation of hydrodynamic parameters, and a fitting function. The results of this study contribute to the understanding of parameter-varying hydrodynamic effects in variable geometry oscillating WECs. The findings provide insights into the potential for reducing structural loads and improving the overall performance of such devices. Further research and development in this area could lead to advancements in WEC technologies, enabling their integration into the competitive energy market.

analytical models↗

DOE OSTI2023

Probabilistic feasibility assessment of sequestration reliance for climate targets

Countries worldwide are transforming their energy systems to achieve Carbon-Neutrality. Investing in renewable resources-based technologies and implementing Carbon Capture and Storage (CCS) are common strategies to achieve higher sequestration levels. Negative emissions through Bioenergy with CCS are expected to play an essential role in the transition to full decarbonization. On top of that, biomass is a limited resource that depends on environmental factors, which create uncertainties related to the amount that can be sustainably provided to the energy system. Here this study emphasizes the relevance of variability in carbon sequestration for achieving climate targets by 2050. Here this paper proposes a probabilistic approach that integrates the Global Change Analysis Model for Latin America (GCAM-LA) with a chance constraint approach. GCAM-LA is used to assess the impact on the energy sector of different limits of sustainable biomass and carbon budget scenarios. The risk associated with exceeding the sequestration capacity of a given region is modeled via Chance Constraint. Results show that electrification is an appropriate long-term decarbonization strategy. It smoothes the effects of uncertainty in sequestration capacity and responds to end-user demands. For this case study, higher levels of electrification are obtained at likelihood levels >66% for end-use sectors.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Design Principles for the Synthesis of Self-Pillared ZSM-5 Zeolite Nanosheets

The design of next-generation materials for emerging energy and environmental applications heavily relies on empirical approaches to direct nonclassical nucleation and crystal growth pathways, polymorphism, intercrystalline transformations, and seed-assisted growth processes. A long-standing obstacle to nanoporous materials design is the complexity of their crystallization, which hinders the development of predictive models and/or physical descriptors that can guide their synthesis. In this study, we use a combination of state-of-the-art synthesis, characterization, and computational design to prepare hierarchical MFI-type zeolites, which we couple with benchmark catalytic testing to assess structure-property-performance relationships. These hierarchical materials are intergrowths of two commercially relevant zeolite frameworks, MFI and MEL, prepared as self-pillared pentasil (SPP) zeolites through seed-assisted, organic-free syntheses for which little theoretical guidance existed. Here, by comparing a large library of zeolite seeds with different pore sizes, dimensions, and structural composite building units, we determined the relative impact of seed and silica source selection, among other synthesis variables. Combined experimental and computational studies are used to test several hypotheses in literature to rationalize the choice of seed structure and establish a more robust selection criteria for seed-assisted synthesis of zeolites. Specifically, we show that a data-driven approach to develop structural descriptors correlates to new, facile routes to rationally design SPP zeolites, addressing knowledge gaps in the fundamental understanding of (non)classical crystal growth mechanisms that are characteristic of nanoporous aluminosilicates.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Cohesive and adhesive degradation in PET-based photovoltaic backsheets subjected to ultraviolet accelerated weathering

Delamination between layers in photovoltaic (PV) backsheets is often reported in the literature, causing voids that can collect moisture, diminish module backside heat dissipation, and reduce the backsheet’s effectiveness as a physical barrier. While backsheets with weathering-resistant fluoropolymer outer layers have traditionally been used in modules, more recent backsheets using non-fluoropolymer outer layers, such as polyethylene terephthalate (PET), have been developed. These backsheets have shown signs of premature degradation, and their adhesion degradation, in particular, has not been widely studied. In the present work, the single cantilever beam (SCB) adhesion test was used to quantify the adhesion energy in two commercially available PET-based backsheets. To study the effect of minor changes in formulation, the backsheets were obtained from same manufacturer and product line but during different years. To study the effect of environmental variables on adhesion degradation, the backsheets were subjected to artificial weathering at controlled temperature, humidity, and ultraviolet (UV) radiation in an indoor weathering chamber, and the adhesion energy was quantified at several intervals of exposure time. Layering structure, composition, and adhesion failure mode were compared between the backsheets, using Raman and infrared spectroscopy and thermogravimetric analysis. The results show a large difference in initial (unexposed) adhesion energy between the backsheets, despite very similar structures and compositions. Following exposure, adhesion energy dropped significantly, primarily due to thermo-hydrolytic degradation of the polyurethane (PU)-based adhesive layers. Significant UV-induced adhesion degradation of the PET outer layer surface was also observed. The study represents an important step in understanding adhesion degradation in PET-based backsheets, suggesting ways in which adhesion integrity – and, correspondingly, module service life – can be improved.

14 SOLAR ENERGY↗

DOE OSTI2020

Electric Power Grid and Natural Gas Network Operations and Coordination

In this report, we present findings from three studies related to the coordination of natural gas and electricity system operations. We first propose and demonstrate a modeling platform for examining the interdependence of natural gas and electricity networks based on a direct current unit-commitment and economic dispatch model for the power system and a transient hydraulic gas model for the gas system. We use this platform to analyze the value of day-ahead coordination of power and natural gas network operations and to show the importance of considering gas system constraints when analyzing power systems operation with high penetration of gas generators and variable renewable energy sources. In the second study, we utilize our modeling platform to consider the U.S. Federal Energy Regulatory Commission (FERC) Order 809, issued in 2015 to improve day-ahead and intraday coordination of power and gas systems. Finally, in the third study we expand our modeling platform to focus on market-based coordination of electricity and natural gas system operations for a real system, namely a subset of the power and gas networks in the Front Range region of Colorado. We use real system data to evaluate the benefits of coordination operations under different conditions, including different levels of renewable penetration and the use of time-variant, shaped flow nominations. Our results indicate that coordination at various timescales can contribute to a reduction in curtailed gas in high-stress periods (such as those with large ramps in gas offtakes) and a reduction in energy consumption of gas compressor stations. We find that intraday coordination can reduce total power system production costs and natural gas deliverability constraints, yielding cost and reliability benefits. We observe these benefits for the test system as well as in the Colorado case study, where we find that coordination and shaped flows may provide additional value for systems with high penetration of variable renewable energy. Together, these three studies demonstrate a pathway for integrated gas and electricity grid modeling and for studying the benefits of coordinated operations of these increasingly interdependent energy systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2024

A survey study on arsenic speciation in coal fly ash and insights into the role of coal combustion conditions

Coal fly ashes (CFAs) are the low-density byproducts of the coal combustion process. Improper or uncontrolled CFA disposal poses significant environmental and health concerns due to the potential leaching of toxic heavy metals such as arsenic (As). Previous studies have investigated the content and speciation of As in different CFA samples, yet systematic information on As speciation in CFA with representative coal source and combustion conditions is still missing. Based on a recent survey study on the typical coal sources and combustion conditions across the U.S., this study selected 19 representative CFA samples to systematically investigate As speciation and potential correlations with these parameters. The composition, morphology, mineralogy, and As speciation of these CFA samples were characterized by complementary analytical, microscopic, and spectroscopic techniques. Synchrotron X-ray spectroscopy and microscopy analyses revealed the dominant As oxidation state to be As(V) and with strong associations to Ca, with the exception of 3 samples that had 19–51% As(III), likely due to the use of selective catalytic reduction (SCR) process. Principal component analysis was conducted to identify potential correlations of As concentration and oxidation state with parameters such as major element content, loss on ignition (LOI), average particle size, coal source, and combustion condition. Al 2 O 3 and FeO content were found to capture a majority of the variability. Further, results from this study provide fundamental basis for understanding the correlations between coal source, combustion conditions, CFA characteristics, and As speciation, and providing insights for downstream beneficial utilization or disposal management.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2025

Key Environmental and Ecological Variables of Wetland CH 4 and CO 2 Fluxes Change With Warming

Wetlands are important ecosystems for the global carbon cycle, impacting regional and global methane (CH 4 ) and carbon dioxide (CO 2 ) budgets. This study examines how environmental and ecological variables impact wetland CH 4 flux and net ecosystem exchange of CO 2 (NEE) across 17 sites globally. We also quantified the importance of variables for each wetland type and site at monthly scale under normal and warm temperatures using dominance analysis. We identified soil and air temperature (TS, TA, respectively) as key variables influencing wetland CH4, and latent heat (LE) and shortwave radiation (SW) for NEE under normal and warm conditions. However, the importance of some variables shifted with warming. For predicting the variability of wetland CH4 flux under warming, gross primary productivity (GPP) and LE, replacing wind direction (WD), were dominant variables for tropical swamps, while NEE was important for high-latitude fens and bogs under warm temperatures. For wetland NEE, the role of TA and TS decreased across all wetland types with warming, while vapor pressure deficit (VPD) became more important for mid and high-latitude wetlands. Our results reveal the complex responses of wetland carbon flux to environmental and ecological variables with warming and provide new insights into improving wetland models by incorporating additional variables and accounting for the changing roles of variables in carbon flux under warming.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A thermodynamic analysis of hypervelocity impacts on metals

Hypervelocity impacts on solid surfaces generate and propagate shock waves and create extreme density, pressure, and temperature conditions. Understanding hypervelocity impacts requires a multi-physics approach at the intersection of fluid dynamics, solid mechanics, and quantum theory, which has the capability of spanning multiple different dominating physics regimes. These strong shock conditions found in hypervelocity impacts are studied in the context of high-speed collisions between spacecraft and meteoroids and debris. We study the problem of hypervelocity impacts on metals with the lens of understanding the processes that lead to the creation of an impact generated plasma. In order to characterize the impacting particle as well as its threat to spacecraft from the subsequent radiating plasma an analysis was performed that highlights how the kinetic energy of the impactor is deposited into a target, how the shock wave and shocked material in the target evolve, and how these effects can lead towards the creation of an impact plasma through heating. The analysis suggests that at low velocities the majority of the energy from a hypervelocity impact is stored in reversible elastic energy and does not contribute to strong heating. At higher impact velocities, nucleus and electron thermal components dominate, increasing the energy available for ionization. Overall, this work provides the complete set of state variable initial conditions for any study of hypervelocity impact plasma expansion and its effects.

42 ENGINEERING↗

DOE OSTI2021

Quantifying the influence of natural climate variability on in situ measurements of seasonal total and extreme daily precipitation

Abstract While various studies explore the relationship between individual sources of climate variability and extreme precipitation, there is a need for improved understanding of how these physical phenomena simultaneously influence precipitation in the observational record across the contiguous United States. In this work, we introduce a single framework for characterizing the historical signal (anthropogenic forcing) and noise (natural variability) in seasonal mean and extreme precipitation. An important aspect of our analysis is that we simultaneously isolate the individual effects of seven modes of variability while explicitly controlling for joint inter-mode relationships. Our method utilizes a spatial statistical component that uses in situ measurements to resolve relationships to their native scales; furthermore, we use a data-driven procedure to robustly determine statistical significance. In Part I of this work we focus on natural climate variability: detection is mostly limited to DJF and SON for the modes of variability considered, with the El Niño/Southern Oscillation, the Pacific–North American pattern, and the North Atlantic Oscillation exhibiting the largest influence. Across all climate indices considered, the signals are larger and can be detected more clearly for seasonal total versus extreme precipitation. We are able to detect at least some significant relationships in all seasons in spite of extremely large (> 95%) background variability in both mean and extreme precipitation. Furthermore, we specifically quantify how the spatial aspect of our analysis reduces uncertainty and increases detection of statistical significance while also discovering results that quantify the complex interconnected relationships between climate drivers and seasonal precipitation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Impact of Dust-Cloud-Radiation-Precipitation Dynamical Feedback on Subseasonal-to-Seasonal Variability of the Asian Summer Monsoon in Global Variable-Resolution Simulations With MPAS-CAM5

In this study, we investigate the effects of increased dust emission from the Middle East deserts on subseasonal-to-seasonal (S₂S) variability of the Asian summer monsoon (ASM). Numerical experiments are performed using the Model for Prediction Across Scales (MPAS) coupled with the Community Atmosphere Model (CAM5) physics, with regional refinement at 30 km grid resolution over South Asia and the surrounding regions. Result shows that increased dust emission and transport from the Middle East/West Asia region induces a strong dust-cloud-radiation-precipitation-circulation feedback, resulting in a colder surface over the desert regions and the western Tibetan Plateau, but warmer and moister troposphere with enhanced cloudiness and precipitation over the Pakistan/Northwest India (PNWI) region. The latter changes are amplified by the dust aerosol induced Elevated Heat Pump (EHP) mechanism along the West Himalayas/Iranian Plateau foothill regions, most pronounced during May-June. During July-August, cloud radiation feedback further enhances the warming of the upper troposphere, and cooling of the land surface over the PNWI and adjacent regions over West Asia. The upper tropospheric heating and increased precipitation over PNWI spur a large-scale anomalous Rossby wavetrain and a northward displacement of the subtropical jetstream, manifesting in a contraction of the South Asian High and westward shift of the Western Pacific Subtropical High. As a result, the entire ASM precipitation-cloud system is displaced westward. Precipitation and cloudiness are intensified over northwest and western India and west Asia, but suppressed over southern and central East Asia. Analyses of the S₂S variability of the upper level vorticity balance suggests that heating by Middle East dust plays an important role in exciting, and anchoring a teleconnection pattern through interactions among dust-cloud radiation, precipitation heating over the PNWI, the development of an upper level Rossby wavetrain, and the northward shift of the boreal summer jetstream over Eurasia.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Salt Sample Statistical Study

Historically there has, at times, been large variability in measured uranium concentration amongst salt samples taken at the same time from the electrorefiner with relative standard deviations ranging from 0.59% up to a high of 114%. This has led to uncertainty in the actual uranium content of the salt. A series of surrogate experiments were performed with gadolinium to study the effects of Li2O, Gd2O3, and Gd metal additions to an LiCl-KCl-GdCl3 salt on the gadolinium content of the salt. These experiments showed no increase in salt sample variability following the additions of the mentioned species, indicating that the carryover of these species from the oxide reduction process may not lead to the historic variability issues.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Does the leaf economic spectrum hold within plant functional types? A Bayesian multivariate trait meta-analysis

The leaf economic spectrum is a widely-studied axis of plant trait variability that defines a trade-off between leaf longevity and productivity. While this has been investigated at the global scale, where it is robust, and at local scales, where deviations from it are common, it has received less attention at the intermediate scale of PFTs. We investigated whether global leaf economic relationships are also present within the scale of plant functional types (PFTs) commonly used by Earth System models, and the extent to which this global-PFT hierarchy can be used to constrain trait estimates. We developed a hierarchical multivariate Bayesian model that assumes separate means and covariance structures within and across PFTs and fit this model to seven leaf traits from the TRY database related to leaf longevity, morphology, biochemistry, and photosynthetic metabolism. Although patterns of trait covariation were generally consistent with the leaf economic spectrum, we found three approximate tiers to this consistency. Relationships among morphological and biochemical traits (SLA, N, P) were the most robust within and across PFTs, suggesting that covariation in these traits is driven by universal leaf construction trade-offs and stoichiometry. Relationships among metabolic traits (R_d, V_{c,max}, J_max) were slightly less consistent, reflecting in part their much sparser sampling (especially for high-latitude PFTs), but also pointing to more flexible plasticity in plant metabolistm. Finally, relationships involving leaf lifespan were the least consistent, indicating that leaf economic relationships related to leaf lifespan are dominated by across-PFT differences and that within-PFT variation in leaf lifespan is more complex and idiosyncratic. Across all traits, these covariance were an important source of information, as evidenced by the improved imputation accuracy and reduced predictive uncertainty in multivariate models compared to univariate models. Ultimately, our study reaffirms the value of studying not just individual traits but the multivariate trait space and the utility of hierarchical modeling for studying the scale dependence of trait relationships.

biodiversity, ecological modeling, ecology, leaf b↗

DOE OSTI2025

A global dataset of terrestrial biological nitrogen fixation

Biological nitrogen fixation (BNF) is the main natural source of new nitrogen inputs in terrestrial ecosystems, supporting terrestrial productivity, carbon uptake, and other Earth system processes. We assembled a comprehensive global dataset of field measurements of BNF in all major N-fixing niches across natural terrestrial biomes derived from the analysis of 376 BNF studies. The dataset comprises 32 variables, including site location, biome type, N-fixing niche, sampling year, quantification method, BNF rate (kg N ha −1 y −1 ), the percentage of nitrogen derived from the atmosphere (%N dfa ), N fixer or N-fixing substrate abundance, BNF rate per unit of N fixer abundance, and species identity. Overall, the dataset combines 1,207 BNF rates for trees, shrubs, herbs, soil, leaf litter, woody litter, dead wood, mosses, lichens, and biocrusts, 152 herb %N dfa values, 1,005 measurements of N fixer or N-fixing substrate abundance, and 762 BNF rates per unit of N fixer abundance for a total of 424 species across 66 countries. This dataset facilitates synthesis, meta-analysis, upscaling, and model benchmarking of BNF fluxes at multiple spatial scales.

Reis Ely, Carla R. [Oregon State Univ., Corvallis,↗

DOE OSTI2025

Transition Zones at the Changing Coastal Terrestrial‐Aquatic Interface

Coastal soils are a significant but highly uncertain component of global biogeochemical cycles. These systems experience spatial and temporal variability in biogeochemical processes, driven by marsh‐to‐upland gradients and hydrological fluctuations. These fluctuations make it difficult to understand and predict biogeochemical processes in these highly dynamic systems. We studied coastal soil biogeochemistry and its variability (a) at regional scales and (b) across transects from upland forest to marsh, in two contrasting regions—Lake Erie, a freshwater lacustrine system, and Chesapeake Bay, a saltwater estuarine system. Salinity‐related analytes were a key source of variability in soil biogeochemistry, not just in the saltwater system, but surprisingly, also in the freshwater system. We had hypothesized linear trends in biogeochemical parameters along the TAI—however, contrary to expectations, transition soils were not consistently intermediate between upland and marsh endmembers; the non‐monotonic trends of C, P, Fe along our transects suggest that these do not behave as expected and may be difficult to model and predict—thus these are key analytes to study in our regions. Rapidly changing soil factors across coastal gradients (e.g., Ca, K, CEC, and TS) may act as precursors to ecosystem shifts. Our comprehensive soil characterization represents a snapshot of a single timepoint of surface soils and provides essential data for mechanistic modeling of ecosystem dynamics across coastal transects.

Patel, Kaizad F. [Pacific Northwest National Labor↗

DOE OSTI2021

Empirical estimation of weather-driven yield shocks using biophysical characteristics for U.S. rainfed and irrigated maize, soybeans, and winter wheat

Agricultural yields are highly susceptible to changes in weather system patterns, including both annual and sub-annual changes in temperature and precipitation. Understanding the impacts of future changes in meteorological variables on crop yields have been widely studied in both empirical and process-based models. This work presents a structural econometric approach that combines historical weather data with information about the biophysical growth cycles of maize, winter wheat, and soy to predict the year-to-year yield responses for rainfed crops. Modeled soil moisture and temperature are taken as key predictors, which allows testing the fitted rainfed model’s ability to predict irrigated yields by assuming irrigation produces the numerically optimal level of soil moisture. This approach is grounded in known biophysical processes, which improves confidence in the model predictions. It produces estimates of interannual variability in yields, which enable study of the effects of extreme events on agricultural productivity. Finally, it enables prediction of the potential impacts of changing weather patterns on irrigated crops in areas that are currently primarily rainfed. We present the results of the empirical model, fitted with rainfed data; out-of-sample validation on irrigated crops; and projections of yield shocks under multiple future climate scenarios. Under a bias-corrected GFDL RCP8.5 scenario, this approach predicts average yield change, relative to 2006-2020 average yields, across U.S. counties in the 2040-2060 period of -16.9% and -13.4% for rainfed and irrigated maize, of -19.1% and -18.4% for rainfed and irrigated soy, and of -4.0% and -2.0% for rainfed and irrigated winter wheat. And in the 2070-2090 period of -30.3% and -28.8 % for rainfed and irrigated maize, of -33.1% and -32.8% for rainfed and irrigated soy, and of -6.5% and -5.2%for rainfed and irrigated winter wheat.

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↗