Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatiotemporal analysis”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1

Abstract. Groundwater plays a key role in meeting water demands, supplying over 40 % of irrigation water globally, with this role likely to grow as water demands and surface water variability increase. A better understanding of the future role of groundwater in meeting sectoral demands requires an integrated hydro-economic evaluation of its cost and availability. Yet substantial gaps remain in our knowledge and modeling capabilities related to groundwater availability, recharge, feasible locations for extraction, extractable volumes, and associated extraction costs, which are essential for large-scale analyses of integrated human–water system scenarios, particularly at the global scale. To address these needs, we developed Superwell, a physics-based groundwater extraction and cost accounting model that operates at sub-annual temporal and at the coarsest 0.5° (≈50 km × 50 km) gridded spatial resolution with global coverage. The model produces location-specific groundwater supply–cost curves that provide the levelized cost to access different quantities of available groundwater. The inputs to Superwell include recent high-resolution hydrogeologic datasets of permeability, porosity, aquifer thickness, depth to water table, recharge, and hydrogeological complexity zones. It also accounts for well capital and maintenance costs, as well as the energy costs required to lift water to the surface. The model employs a Theis-based scheme coupled with an amortization-based cost accounting formulation to simulate groundwater extraction and quantify the cost of groundwater pumping. The result is a spatiotemporally flexible, physically realistic, economics-based model that produces groundwater supply–cost curves. We show examples of these supply–cost curves and the insights that can be derived from them across a set of scenarios designed to explore model outcomes. The supply–cost curves produced by the model show that most (90 %) nonrenewable groundwater in storage globally is extractable at costs lower than USD 0.57 m−3, while half of the volume remains extractable at under USD 0.108 m−3. The global unit cost is estimated to range from a minimum of USD 0.004 m−3 to a maximum of USD 3.971 m−3. We also demonstrate and discuss examples of how these cost curves could be used by linking Superwell's outputs with other models to explore coupled human–environmental system challenges, such as water resources planning and management, or broader analyses of multisectoral feedbacks.

Global Change Analysis Model (GCAM)↗
DOE OSTI2023

Insights from application of a hierarchical spatio-temporal model to an intensive urban black carbon monitoring dataset

Existing regulatory pollutant monitoring networks rely on a small number of centrally located measurement sites that are purposefully sited away from major emission sources. While informative of general air quality trends regionally, these networks often do not fully capture the local variability of air pollution exposure within a community. Recent technological advancements have reduced the cost of sensors, allowing air quality monitoring campaigns with high spatial resolution. The 100×100 black carbon (BC) monitoring network deployed 100 low-cost BC sensors across the 15 km 2 West Oakland, CA community for 100 days in the summer of 2017, producing a nearly continuous site-specific time series of BC concentrations which we aggregated to one-hour averages. Leveraging this dataset, we employed a hierarchical spatio-temporal model to accurately predict local spatio-temporal concentration patterns throughout West Oakland, at locations without monitors (average cross-validated hourly temporal R 2 =0.60). Using our model, we identified spatially varying temporal pollution patterns associated with small-scale geographic features and proximity to local sources. In a sub-sampling analysis, here we demonstrated that fine scale predictions of nearly comparable accuracy can be obtained with our modeling approach by using ~30% of the 100×100 BC network supplemented by a shorter-term high-density campaign.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Using Explainable Artificial Intelligence to Predict Perovskite Solar Cell Electrical Metastability from Operando Photoluminescence Images in Accelerated Stress Testing

Metal halide perovskite (MHP) solar cells exhibit a metastable response to bias governed by coupled ionic–electronic processes, complicating the conventional reciprocity relation between luminescence intensity and device open-circuit voltage (V oc ). This limits the use of luminescence as a diagnostic for device screening or accelerated stress testing, motivating new approaches that can interpret photoluminescence (PL) signals under nonequilibrium conditions. From the artificial intelligence perspective, we develop an explainable deep learning framework that integrates convolutional neural networks (CNN), long short-term memory (LSTM) layers, and an attention mechanism to learn spatiotemporal features from operando photoluminescence PL image sequences. The model achieves a mean absolute error of ±0.027 V in predicting open-circuit voltage transients and reduces extreme-tail errors by up to 78% compared to physics-based reciprocity calculations. Gradient-weighted Class Activation Mapping (Grad-CAM) provides interpretability by highlighting physically meaningful regions such as electrode edges and emergent defect features. From the engineering application perspective, this framework enables accurate, contactless prediction of device V oc and identification of degradation-relevant features during accelerated aging of perovskite solar cells. This approach demonstrates how explainable AI can enhance operando diagnostics and reliability analysis in photovoltaic devices under nonequilibrium conditions.

14 SOLAR ENERGY↗
DOE OSTI2021

Evaluating Diurnal and Semi-Diurnal Cycle of Precipitation in CMIP6 Models Using Satellite- and Ground-Based Observations

We report the diurnal and semi-diurnal cycle of precipitation simulated from CMIP6 models during 1996-2005 are evaluated globally between 60°S and 60°N, as well as at ten selected locations representing three categories of diurnal cycle of precipitation: (1) afternoon precipitation over land, (2) early morning precipitation over ocean, and (3) nocturnal precipitation over land. Three satellite-based and two ground-based rainfall products are used to evaluate the climate models. Globally, the ensemble mean of CMIP6 models shows a diurnal phase of 3 to 4 hours earlier over land and 1 to 2 hours earlier over ocean, when compared with the latest satellite products. These biases are in line with what were found in previous versions of climate models but reduced compared to the CMIP5 ensemble mean. Analysis at the selected locations complimented with in-situ measurements further reinforces these results. Several CMIP6 models have shown a significant improvement in the diurnal cycle of precipitation compared to their CMIP5 counterparts, notably on delaying afternoon precipitation over land. This can be attributed to the use of more sophisticated convective parameterizations. Most models are still unable to capture the nocturnal peak associated with elevated convection and propagating mesoscale convective systems, with a few exceptions that allow convection to be initiated above the boundary layer to capture nocturnal elevated convection. We also quantify an encouraging consistency between the satellite- and ground-based precipitation measurements despite differing spatiotemporal resolutions and sampling periods, which provides confidence in using them to evaluate the diurnal and semi-diurnal cycle of precipitation in climate models.

58 GEOSCIENCES↗
DOE OSTI2024

Characterizing model uncertainties in simulated coast-to-offshore wind over the northeast U.S. using multi-platform measurements from the TCAP field campaign

Numerical weather prediction (NWP) models, such as the Weather Research and Forecasting (WRF) model, are widely used to provide estimates of the offshore wind energy resource owing to their large spatial coverage compared to available observations. Nevertheless, spatiotemporal distribution of model biases is highly dependent on factors including model configuration, location, and the interplay of multi-scale physical processes. Here, in this study, we focus on the characterization of model uncertainties in simulated coast-to-offshore winds over the northeast U.S., by varying sea surface temperature (SST) forcings, surface layer (SL) and planetary boundary layer (PBL) parameterizations, as well as identifying biases that may be directly passed from initial and boundary conditions. Multiple measurements, including aircraft data collected during the U.S. Department of Energy's Two-Column Aerosol Project (TCAP) experiment, are used to constrain the model results and facilitate quantitative comparisons. Our analysis indicates while SST forcing has notable impacts on simulated air temperature and moisture within PBL, the modeled winds are in general more sensitive to the choices of SL and PBL physics than to SST. The model’s forcing data not only controls the vertical dependence of wind speed errors, but also alters regional variability in wind speed’s spatial correlation. Bias comparisons between ERA5 reanalysis and ensemble simulations revealed significant similarity, particularly in wind speed biases during winter, underscoring their dependency on initial and boundary conditions. Coastal and offshore near-surface wind speed biases tend to exhibit much higher similarity in winter than in summer due to the presence of much stronger and more persistent synoptic wind conditions. This study highlights the importance of accurate atmospheric forcing and parameterization choices in improving wind forecasts and suggests the potential for extrapolating coastal wind biases to offshore locations, aiding wind energy forecasting and informing the Wind Forecast Improvement Project-3 (WFIP3).

17 WIND ENERGY↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Acoustic cavities in 2D heterostructures

Two-dimensional (2D) materials offer unique opportunities in engineering the ultrafast spatiotemporal response of composite nanomechanical structures. In this work, we report on high frequency, high quality factor ( Q ) 2D acoustic cavities operating in the 50–600 GHz frequency ( f ) range with f × Q up to 1 × 10 14 . Monolayer steps and material interfaces expand cavity functionality, as demonstrated by building adjacent cavities that are isolated or strongly-coupled, as well as a frequency comb generator in MoS 2 /h-BN systems. Energy dissipation measurements in 2D cavities are compared with attenuation derived from phonon-phonon scattering rates calculated using a fully microscopic ab initio approach. Phonon lifetime calculations extended to low frequencies (<1 THz) and combined with sound propagation analysis in ultrathin plates provide a framework for designing acoustic cavities that approach their fundamental performance limit. These results provide a pathway for developing platforms employing phonon-based signal processing and for exploring the quantum nature of phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine Learning for Daily Forecasts of Arctic Sea Ice Motion: An Attribution Assessment of Model Predictive Skill

Physics-based simulations of Arctic sea ice are highly complex, involving transport between different phases, length scales, and time scales. Resultantly, numerical simulations of sea ice dynamics have a high computational cost and model uncertainty. We employ data-driven machine learning (ML) to make predictions of sea ice motion. The ML models are built to predict present-day sea ice velocity given present-day wind velocity and previous-day sea ice concentration and velocity. Models are trained using reanalysis winds and satellite-derived sea ice properties. We compare the predictions of three different models: persistence (PS), linear regression (LR), and a convolutional neural network (CNN). We quantify the spatiotemporal variability of the correlation between observations and the statistical model predictions. Additionally, we analyze model performance in comparison to variability in properties related to ice motion (wind velocity, ice velocity, ice concentration, distance from coast, bathymetric depth) to understand the processes related to decreases in model performance. Results indicate that a CNN makes skillful predictions of daily sea ice velocity with a correlation up to 0.81 between predicted and observed sea ice velocity, while the LR and PS implementations exhibit correlations of 0.78 and 0.69, respectively. The correlation varies spatially and seasonally: lower values occur in shallow coastal regions and during times of minimum sea ice extent. LR parameter analysis indicates that wind velocity plays the largest role in predicting sea ice velocity on 1-day time scales, particularly in the central Arctic. Regions where wind velocity has the largest LR parameter are regions where the CNN has higher predictive skill than the LR.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Constraining Black Carbon Aging in Global Models to Reflect Timescales for Internal Mixing

The radiative effects of black carbon depend critically on its atmospheric lifetime, which is controlled by the rate at which freshly emitted combustion particles become internally mixed with other aerosol components. Global aerosol models strive to represent this process, but the timescale for aerosol mixing is not easily constrained using observations. In this study, we apply a timescale parameterization derived from particle‐resolved simulations to quantify, in a global aerosol model, the timescale for internal mixing. We show that, while highly variable, the average timescale for internal mixing is approximately 3 hr, which is much shorter than the 24‐hr aging timescale traditionally applied in bulk aerosol models. We then use the mixing timescale to constrain the aging criterion in the Modal Aerosol Module. Our analysis reveals that, to best reflect timescales for internal mixing, modal models should assume that particles transition from the hydrophobic (fresh) to the hydrophilic (aged) class once they accumulate a coating thickness equal to four monolayers of sulfuric acid, as opposed to the model's current aging criterion of eight monolayers. We show that, in remote regions like the Arctic and Antarctic, predictions of black carbon loading and its seasonal variation are particularly sensitive to the model representation of aging. By constraining aging in global models to reflect mixing timescales simulated by the particle‐resolved model, we eliminate one of the free parameters governing black carbon's long‐range transport and spatiotemporal distribution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Classical Non-Markovian Noise in Symmetry-Preserving Quantum Dynamics

In quantum dynamics, symmetries are vital for identifying and assessing conserved quantities that govern the evolution of a quantum system. When promoted to the open quantum system setting, dynamical symmetries can be negatively altered by system-environment interactions, thus, complicating their analysis. Previous work on noisy symmetric quantum dynamics has focused on the Markovian setting, despite the ubiquity of non-Markovian noise in a number of widely used quantum technologies. Here, in this Letter, we develop a framework for quantifying the impact of non-Markovian noise on symmetric quantum evolution via root space decompositions and the filter function formalism. We demonstrate analytically that symmetry-preserving noise maintains the symmetric subspace, while nonsymmetric noise leads to highly specific leakage errors that are block diagonal in the symmetry representation. We support our findings with numerical studies of a transverse-field Ising model and a quantum error detecting code subject to spatiotemporally correlated multiaxis noise. Our results are broadly applicable, providing new analytic insights into the control and characterization of open quantum system dynamics.

decoherence↗
DOE OSTI2020

Landfalling Droughts: Global Tracking of Moisture Deficits From the Oceans Onto Land

Abstract Droughts threaten food, energy, and water security, causing death and displacement of millions of people and billions of dollars in damages. However, there are still important gaps in the understanding of drought mechanisms and behaviors, inhibiting the accuracy of early‐warning systems designed to protect communities worldwide. We use an object‐tracking algorithm to track clusters of precipitation‐minus‐evaporation moisture deficits across land and ocean areas of the globe from 1981–2018. This analysis reveals a new type of “landfalling drought” that originates over the ocean and “migrates” onto land. We find that 16% of droughts that affected the continents worldwide from 1981–2018 were landfalling droughts. These droughts were significantly larger (220–425%) and more intense (4–30%)—and grew (253–285%) and intensified (9–28%) faster—than droughts that developed solely over the land or ocean. To identify potential underlying mechanisms, we analyze moisture transport associated with landfalling droughts over western North America. We find that landfalling droughts in this region are associated with anomalously anticyclonic atmospheric pressure patterns that reduce moisture fluxes over the Pacific Ocean toward the continent. By advancing understanding of the spatiotemporal evolution of droughts, our findings offer the potential to improve seasonal‐scale prediction and long‐term projection of global drought risks.

Herrera‐Estrada, Julio E.↗
DOE OSTI2023

Picosecond magneto-optic thermometry measurements of nanoscale thermal transport in AlN thin films

The thermal conductivity Λ of wide bandgap semiconductor thin films, such as AlN, affects the performance of high-frequency devices, power devices, and optoelectronics. However, accurate measurements of Λ in thin films with sub-micrometer thicknesses and Λ > 100 W m -1 K -1 is challenging. Widely used pump/probe metrologies, such as time–domain thermoreflectance (TDTR) and frequency–domain thermoreflectance, lack the spatiotemporal resolution necessary to accurately quantify thermal properties of sub-micrometer thin films with high Λ. In this work, we use a combination of magneto-optic thermometry and TiN interfacial layers to significantly enhance the spatiotemporal resolution of pump/probe thermal transport measurements. We use our approach to measure Λ of 100, 400, and 1000 nm AlN thin films. We coat AlN thin films with a ferromagnetic thin-film transducer with the geometry of (1 nm-Pt/0.4 nm-Co) x3 /(2 nm-TiN). This PtCo/TiN transducer has a fast thermal response time of <50 ps, which allows us to differentiate between the thermal response of the transducer, AlN thin film, and substrate. For the 100, 400, and 1000 nm thick AlN films, we determine Λ to be 200 ± 80, 165 ± 35, and 300 ± 70 W m -1 K -1 , respectively. We conclude with an uncertainty analysis that quantifies the errors associated with pump/probe measurements of thermal conductivity, as a function of transducer type, thin-film thermal conductivity, and thin-film thickness. Time resolved magneto-optic Kerr effect experiments can measure films that are three to five times thinner than is possible with standard pump/probe metrologies, such as TDTR. This advance in metrology will enable better characterization of nanoscale heat transfer in high thermal conductivity material systems like wide bandgap semiconductor heterostructures and devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spatiotemporal Variability of Electric System Reliability Metrics

As urbanization continues to grow, urban density also rises, placing significant stress on critical infrastructure. Cities often face challenges in investing in new infrastructure or upgrading existing ones. One area of infrastructure that has come under strain is the electric grid, leading to prolonged power outages. According to data from the U.S. Energy Information Administration, the average duration of power outages has steadily doubled from 2013 to 2021, primarily due to extreme weather events. The IEEE has developed a series of reliability metrics to assess the reliability of the electric system. However, these metrics are only applied at the utility level, making it challenging to comprehend the variations in these indices at finer spatial resolutions and, consequently, evaluate reliability at the non-utility level. In this study, we leveraged electric outage data collected at 15-minute intervals at the county level to estimate electric system reliability over a six-year period using a novel data analysis approach. Our findings reveal a gradual decline in reliability in the U.S. Southwest compared to national averages. Furthermore, the duration of outages per customer per event has significantly decreased in large portions of California and Florida. These discoveries can be valuable for utility companies as they seek to pinpoint regional anomalies and plan for prioritized remediation efforts.

Singh, Nagendra↗
DOE OSTI2022

Exploratory Precipitation Metrics: Spatiotemporal Characteristics, Process-Oriented, and Phenomena-Based

Abstract Precipitation sustains life and supports human activities, making its prediction one of the most societally relevant challenges in weather and climate modeling. Limitations in modeling precipitation underscore the need for diagnostics and metrics to evaluate precipitation in simulations and predictions. While routine use of basic metrics is important for documenting model skill, more sophisticated diagnostics and metrics aimed at connecting model biases to their sources and revealing precipitation characteristics relevant to how model precipitation is used are critical for improving models and their uses. This paper illustrates examples of exploratory diagnostics and metrics including 1) spatiotemporal characteristics metrics such as diurnal variability, probability of extremes, duration of dry spells, spectral characteristics, and spatiotemporal coherence of precipitation; 2) process-oriented metrics based on the rainfall–moisture coupling and temperature–water vapor environments of precipitation; and 3) phenomena-based metrics focusing on precipitation associated with weather phenomena including low pressure systems, mesoscale convective systems, frontal systems, and atmospheric rivers. Together, these diagnostics and metrics delineate the multifaceted and multiscale nature of precipitation, its relations with the environments, and its generation mechanisms. The metrics are applied to historical simulations from phases 5 and 6 of the Coupled Model Intercomparison Project. Models exhibit diverse skill as measured by the suite of metrics, with very few models consistently ranked as top or bottom performers compared to other models in multiple metrics. Analysis of model skill across metrics and models suggests possible relationships among subsets of metrics, motivating the need for more systematic analysis to understand model biases for informing model development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Ensemble learning-iterative training machine learning for uncertainty quantification and automated experiment in atom-resolved microscopy

Deep learning has emerged as a technique of choice for rapid feature extraction across imaging disciplines, allowing rapid conversion of the data streams to spatial or spatiotemporal arrays of features of interest. However, applications of deep learning in experimental domains are often limited by the out-of-distribution drift between the experiments, where the network trained for one set of imaging conditions becomes sub-optimal for different ones. This limitation is particularly stringent in the quest to have an automated experiment setting, where retraining or transfer learning becomes impractical due to the need for human intervention and associated latencies. Here we explore the reproducibility of deep learning for feature extraction in atom-resolved electron microscopy and introduce workflows based on ensemble learning and iterative training to greatly improve feature detection. This approach allows incorporating uncertainty quantification into the deep learning analysis and also enables rapid automated experimental workflows where retraining of the network to compensate for out-of-distribution drift due to subtle change in imaging conditions is substituted for human operator or programmatic selection of networks from the ensemble. This methodology can be further applied to machine learning workflows in other imaging areas including optical and chemical imaging.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spatial and Temporal Interpolation Analysis Process of Shock Loading

The following document serves to describe the current process that is utilized by W-13 analysts to approximate the pressure-time boundary condition seen by a test object from shock tube loading. The intent of this report is to capture the current capability and to encourage continued growth and development. Information on the required input variables, mapping process, and generated output information is covered in the following sections. Throughout the document, gaps and weaknesses of the mapping process are noted to encourage future development efforts.

42 ENGINEERING↗
DOE OSTI2024

Bias Corrected NOAA HRRR Wind Resource Data for Grid Integration Applications

To address the need for regularly updated wind resource data, NREL has processed the High-Resolution Rapid Refresh (HRRR) outputs for use in grid integration modeling. The HRRR is an hourly-updated operational forecast product produced by the National Oceanic and Atmospheric Administration (NOAA) (Dowell et al., 2022). Several barriers have prevented the HRRR's widespread proliferation in the wind energy industry: missing timesteps (prior to 2019), challenging file format for wind energy analysis, limited vertical height resolution, and negative bias versus legacy WIND Toolkit data (2007-2013). NREL has applied re-gridding, interpolation, and bias-correction to the native HRRR data to overcome these limitations. This results in the now-publicly-available bias corrected and interpolated HRRR (BC-HRRR) dataset for weather years 2015 to 2023. Bias correction is necessary for wind resource consistency across weather years to be used simultaneously in planning-focused grid integration studies alongside the original WIND Toolkit data. We show that quantile mapping with the WIND Toolkit as a historical baseline is an effective method for bias correcting the interpolated HRRR data: the BC-HRRR has reduced mean bias versus comparable gridded wind resource datasets (+0.12 m/s versus Vortex) and has very low mean bias versus ground measurement stations (+0.01 m/s) (Buster et al., 2024). BC-HRRR's consistency with the legacy WIND Toolkit allows NREL to extend grid integration analysis to 15+ weather years of wind data with low-overhead extensibility to future years as they are made available by NOAA. As with historical datasets like the WIND Toolkit, BC-HRRR is intended for use in grid integration modeling (e.g., capacity expansion, production cost, and resource adequacy modeling) both independently and alongside the legacy WIND Toolkit.

Array↗
DOE OSTI2025

Machine learning and process-based modeling of spatiotemporal changes in active layer thickness across Alaska

Permafrost degradation poses a growing threat to infrastructure stability and ecosystem resilience in the rapidly warming Arctic. We investigated the spatiotemporal dynamics of active layer thickness (ALT) across Alaska by integrating field observations, environmental datasets, a physically based Stefan model, and machine learning (ML) techniques. Using weather projections from the Coupled Model Intercomparison Project Phase 6 under two Shared Socioeconomic Pathways (SSP 2-4.5 and SSP 5-8.5), we assessed ALT sensitivity to projected future weather conditions. The random forest (RF) model outperformed the Stefan approach in predicting ALT on the training dataset (R² = 0.84 vs. 0.53) but demonstrated lower generalizability on the test dataset (R² = 0.24 vs. 0.54). The root mean square error (RMSE) for the RF model for training and testing ranged from 14 to 22 cm, compared to 17 and 18 cm for the Stefan model. Variable importance analysis revealed that mean annual temperature and slope angle were the strongest predictors of ALT, accounting for 19% and 18% of the variance, respectively, followed by sediment transport index (14%) and stream power index (11%). Comparative analysis of baseline ALT predictions showed the Stefan model tended to project a thicker active layer (mean ± SD: 65 ± 16 cm), compared to the RF model (mean ± SD: 59 ± 8.8) cm). Both models indicated a latitudinal gradient in ALT, with shallower depths at higher latitudes. Projected ALT increases by 2100 were estimated at 3.3 ± 2.2 cm under SSP 2-4.5 and 5.9 ± 4.0 cm under SSP 5-8.5 for the ML model, whereas the Stefan model projected substantially larger increases of 13 ± 2.6 cm (SSP 2-4.5) and 28 ± 4.4 cm (SSP5-8.5). Spatial analysis showed the greatest ALT increases in northern Alaska, with relatively smaller changes in southern regions. These findings highlight the complex, multifactorial nature of ALT dynamics and the value of hybrid modeling approaches. As rising temperatures accelerate permafrost thaw, changes in ALT can disrupt ecosystems, damage infrastructures, and enhance the release of stored soil carbon, highlighting the urgent need for improved predictive capabilities to inform adaptation strategies in the Arctic.

Climate sciences↗