Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial force correlation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Transport coefficient approach for characterizing nonequilibrium dynamics in soft matter

Nonequilibrium states in soft condensed matter require a systematic approach to characterize and model materials, enhancing predictability and applications. Among the tools, X-ray photon correlation spectroscopy (XPCS) provides exceptional temporal and spatial resolution to extract dynamic insight into the properties of the material. However, existing models might overlook intricate details. We introduce an approach for extracting the transport coefficient, denoted as $J(t)$, from the XPCS studies. This coefficient is a fundamental parameter in nonequilibrium statistical mechanics and is crucial for characterizing transport processes within a system. Our method unifies the Green–Kubo formulas associated with various transport coefficients, including gradient flows, particle–particle interactions, friction matrices, and continuous noise. We achieve this by integrating the collective influence of random and systematic forces acting on the particles within the framework of a Markov chain. We initially validated this method using molecular dynamics simulations of a system subjected to changes in temperatures over time. Subsequently, we conducted further verification using experimental systems reported in the literature and known for their complex nonequilibrium characteristics. The results, including the derived $J(t)$ and other relevant physical parameters, align with the previous observations and reveal detailed dynamical information in nonequilibrium states. This approach represents an advancement in XPCS analysis, addressing the growing demand to extract intricate nonequilibrium dynamics. Further, the methods presented are agnostic to the nature of the material system and can be potentially expanded to hard condensed matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Plating and Stripping of Lithium Metal Stabilized by a Block Copolymer Electrolyte: Local Current Density Measurement and Modeling

Here, plating and stripping of lithium protrusions in lithium metal symmetric cells containing a solid block copolymer electrolyte was studied as a function of time in 3D using time-resolved X-ray tomography. These measurements enabled determination of the spatial variation in current densities at the plating and stripping electrodes. The initial interelectrode distance was 27 μm. Correlation functions were calculated to reveal the relationships between current densities at the two electrodes and local electrolyte thickness. Current densities at opposing electrode locations during protrusion growth is uncorrelated until the local interelectrode distance decreases to less than 6 μm, just before the cell shorts. Mass balance was used to determine the area from which lithium ions that form a protrusion were stripped. Computational modeling of the plating and stripping process reveals the interplay between electrochemical and mechanical driving forces and their effect on nonuniform current distribution. Model predictions were compared with experiments without resorting to any adjustable parameters. The computed correlation functions were in qualitative agreement with experiments. Finally, the model was used to calculate contour plots of electrochemical potential within the electrolyte, shedding light on how geometry, salt concentration, interelectrode distance, and mechanical stress influence local rates of electrochemical reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Description of the Three-Dimensional Large-Scale Forcing Data from the 3D Constrained Variational Analysis (VARANAL3D)

This technical report introduces a Three-Dimensional Constrained Variational Analysis (3DCVA) (Tang and Zhang 2015) and its product of three-dimensional large-scale forcing data to drive single-column models (SCM), cloud-resolving models (CRM), and large-eddy simulation (LES) models, and to evaluate model results. The 3DCVA algorithm is an extension of the original 1D constrained variational analysis (1DCVA) (Zhang and Lin 1997, Zhang et al. 2001). The three-dimensional structure of the forcing data allows studies of spatial variation of the large-scale forcing fields and tests of physical parameterizations across scales. In the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility, the 3D forcing data are assigned the datastream name varanal3d. In this technical report, 3DCVA will be used to refer to the algorithm, while VARANAL3D will be used to refer to the data product.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

System Identification of Damped Truss-Like Space Structures

A spacecraft payload flown on a launch vehicle experiences dynamic loads. The dynamic loads are caused by various phenomena ranging from the start-up of the launch vehicle engine to wind gusts. A spacecraft payload should be designed to meet launch vehicle dynamic loads. One of the major steps taken towards determining the dynamic loads is to correlate the finite element model of the spacecraft with the test results of a modal survey test. A test-verified finite element model of the spacecraft should possess the same spatial properties (stiffness, mass, and damping) and modal properties (frequencies and mode shapes) as the test hardware representing the spacecraft. The test-verified and correlated finite element model of the spacecraft is then coupled with the finite element model of the launch vehicle for analysis of loads and stress. Modal survey testing, verification of a finite element model, and modification of the finite element model to match the modal survey test results can easily be accomplished if the spacecraft structure is simple. However, this is rarely the case. A simple structure here is defined as a structure where the influence of nonlinearity between force and displacement (uncertainty in a test, for example, with errors in input and output), and the influence of damping (structural, coulomb, and viscous) are not pronounced. The objective of this study is to develop system identification and correlation methods with the focus on the structural systems that possess nonproportional damping. Two approaches to correct the nonproportional damping matrix of a truss structure were studied, and have been implemented on truss-like structures such as the National Aeronautics and Space Administration's space station truss. The results of this study showed nearly 100 percent improvement of the correlated eigensystem over the analytical eigensystem. The first method showed excellent results with up to three modes used in the system identification process. The second method could handle more modes, but required more computer usage time, and the results were less accurate than those of the first method.

Armand, Sasan↗
DOE OSTI2022

Gaussian approximation of dispersion potentials for efficient featurization and machine-learning predictions of metal–organic frameworks

Energy-related descriptors in machine learning are a promising strategy to predict adsorption properties of metal–organic frameworks (MOFs) in the low-pressure regime. Interactions between hosts and guests in these systems are typically expressed as a sum of dispersion and electrostatic potentials. The energy landscape of dispersion potentials plays a crucial role in defining Henry’s constants for simple probe molecules in MOFs. To incorporate more information about this energy landscape, we introduce the Gaussian-approximated Lennard-Jones (GALJ) potential, which fits pairwise Lennard-Jones potentials with multiple Gaussians by varying their heights and widths. The GALJ approach is capable of replicating information that can be obtained from the original LJ potentials and enables efficient development of Gaussian integral (GI) descriptors that account for spatial correlations in the dispersion energy environment. GI descriptors would be computationally inconvenient to compute using the usual direct evaluation of the dispersion potential energy surface. Here, we show that these new GI descriptors lead to improvement in ML predictions of Henry’s constants for a diverse set of adsorbates in MOFs compared to previous approaches to this task.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

The Mexican Drought Atlas: Tree-Ring Reconstructions of the Soil Moisture Balance During the Late Pre-Hispanic, Colonial, and Modern Eras

Mexico has suffered a long history and prehistory of severe sustained drought. Drought over Mexico is modulated by ocean-atmospheric variability in the Atlantic and Pacific, raising the possibility for long-range seasonal climate forecasting, which could help mediate the economic and social impacts of future dry spells. The instrumental record of Mexican climate is very limited before 1920, but tree-ring chronologies developed from old-growth forests in Mexico can provide an excellent proxy representation of the spatial pattern and intensity of past moisture regimes useful for the analysis of climate dynamics and climate impacts. The Mexican Drought Atlas (MXDA) has been developed from an extensive network of 252 climate sensitive tree-ring chronologies in and near Mexico. The MXDA reconstructions extend from 1400 CE-2012 and were calibrated with the instrumental summer (JJA) self-calibrating Palmer Drought Severity Index (scPDSI) on a 0.5deg latitude/longitude grid extending over land areas from 14 to 34degN and 75-120degW using Ensemble Point-by-Point Regression (EPPR) for the 1944-1984 period. The grid point reconstructions were validated for the period 1920-1943 against instrumental gridded scPDSI values based on the fewer weather station observations available during that interval. The MXDA provides a new spatial perspective on the historical impacts of moisture extremes over Mexico during the past 600-years, including the Aztec Drought of One Rabbit in 1454, the drought of El Ano de Hambre in 1785-1786, and the drought that preceded the Mexican Revolution of 1909-1910. The El Nino/Southern Oscillation (ENSO) is the most important ocean-atmospheric forcing of moisture variability detected with the MXDA. In fact, the reconstructions suggest that the strongest central equatorial Pacific sea surface temperature (SST) teleconnection to the soil moisture balance over North America may reside in northern Mexico. This ENSO signal has stronger and more time-stable correlations than computed for either the Atlantic Multidecadal Oscillation or Pacific Decadal Oscillation. The extended Multivariate ENSO Index is most highly correlated with reconstructed scPDSI over northern Mexico, where warm events favor moist conditions during the winter, spring, and early summer. This ENSO teleconnection to northern Mexico has been strong over the past 150 years, but it has been comparatively weak and non-stationary in the MXDA over central and southern Mexico where eastern tropical Pacific and Caribbean/tropical Atlantic SSTs seem to be more important. The ENSO teleconnection to northern Mexico is weaker in the available instrumental PDSI, but analyses based on the millennium climate simulations with the Community Earth System Model suggest that the moisture balance during the winter, spring, and early summer over northern Mexico may indeed be particularly sensitive to ENSO forcing. Nationwide drought is predicted to become more common with anthropogenic climate change, but the MXDA reconstructions indicate that intense "All Mexico" droughts have been rare over the past 600 years and their frequency does not appear to have increased substantially in recent decades.

El Nino↗
DOE OSTI2024

Optical tuning of the diamond Fermi level measured by correlated scanning probe microscopy and quantum defect spectroscopy

Quantum technologies based on quantum point defects in crystals require control over the defect charge state. Here we tune the charge state of shallow nitrogen-vacancy and silicon-vacancy centers by locally oxidizing a hydrogenated surface with moderate optical excitation and simultaneous spectral monitoring. The loss of conductivity and change in work function due to oxidation are measured in atmosphere using conductive atomic force microscopy and Kelvin probe force microscopy (KPFM). We correlate these scanning probe measurements with optical spectroscopy of the nitrogen-vacancy and silicon-vacancy centers created via implantation 15–25 nm beneath the diamond surface and annealing. The observed charge state of the defects as a function of optical exposure demonstrates that laser oxidation provides a way to precisely tune the Fermi level over a range of at least 2.00 eV. We also observe a significantly larger oxidation rate for implanted surfaces compared to unimplanted surfaces under ambient conditions. Here, combined with knowledge of the electron affinity of a surface, these results suggest KPFM is a powerful, high-spatial-resolution technique to advance surface Fermi level engineering for charge stabilization of quantum defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗
DOE OSTI2020

Controlling protein adsorption modes electrostatically

Protein adsorption on surfaces is ubiquitous in biology and in biotechnology. There are various forces required for controlling protein adsorption. Here, we introduce an explicit ion coarse-grained molecular dynamics simulation approach for studying the effects of electrostatics on protein adsorption, and 2D protein assembly on charged surfaces. Our model accounts for the spatial distribution of protein charges. We use catalase as our model protein. We find that the preferential adsorption mode of proteins at low protein concentration on a charged surface is “standing up”. When the protein concentration in a solution increases to reach a critical density on the surface, the adsorption mode switches from “standing up” to a mixed state “flat on” and “standing up”, which increases the lateral correlations among the adsorbed proteins. As such, the changes in the adsorption mode arise from the protein adsorption that cancel the surface charge and the protein–protein repulsion. This correlated surface structure melts as the salt concentration increases because the charged surface is cancelled by the salt ions and the proteins de-adsorb. For the case of strongly charged surfaces the “standing up” conformation remains more favorable even at high protein adsorption at low salt concentrations since in that conformation the surface charge is cancelled more effectively, generating an even more laterally correlated structure. We elucidate the effects of parameters such as surface charge density, salt concentration, and protein charges on the different adsorption modes and the structure and organization of proteins on the charged surfaces. Finally, this study provides a guide for controlling protein assembly on surfaces.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Core/Shell-Like Localized Emission at Atomically Thin Semiconductor–Au Interface

Localized emission in atomically thin semiconductors has sparked significant interest as single-photon sources. Despite comprehensive studies into the correlation between localized strain and exciton emission, the impacts of charge transfer on nanobubble emission remains elusive. Here, in this paper, we report the observation of core/shell-like localized emission from monolayer WSe 2 nanobubbles at room temperature through near-field studies. By altering the electronic junction between monolayer WSe 2 and the Au substrate, one can effectively adjust the semiconductor to metal junction from a Schottky to an Ohmic junction. Through concurrent analysis of topography, potential, tip-enhanced photoluminescence, and a piezo response force microscope, we attribute the core/shell-like emissions to strong piezoelectric potential aided by induced polarity at the WSe 2 –Au Schottky interface which results in spatial confinement of the excitons. Our findings present a new approach for manipulating charge confinement and engineering localized emission within atomically thin semiconductor nanobubbles. These insights hold implications for advancing the nano and quantum photonics with low-dimensional semiconductors.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1981

Planetary-scale wave structure in the Martian atmosphere

Wave-like perturbations are found in the Mariner 9 IRIS atmospheric temperature data during late Northern Hemisphere winter in a latitude band between 45 deg N and 65 deg N. The nature of the data base prevents a unique separation of spatial and temporal behavior, but Fourier analysis of the data constrains the waves to discrete combinations of planetary wavenumber and period. One major spectral component possesses a meridional amplitude cross section with a maximum near the 1 mbar level at 60 deg N and is strongly correlated with the circumpolar jet observed in thermal winds calculated from the mean meridional temperature cross section. This feature is consistent with the low-wavenumber baroclinic waves observed in Viking Lander data, and the vertical structure reflects the behavior anticipated for a vertically penetrating quasi-geostrophic disturbance. Other possible origins for the wave cannot be ruled out, however. Among these is a stationary wave forced by wavenumber-2 topographic relief.

Conrath, B. J.↗
DOE OSTI2020

Correlative nm-Scale Nonuniformity of Active Charge Carriers and Electrical Potential Along Both the Plane-View and Depth Directions in Group-V-Doped CdTe Thin Films: Preprint

We report nanometer-scale imaging on inhomogeneous distributions of active carrier and electrical potential in an As-doped CdTe film along both plane-view and film-depth directions. Despite Se grading, the SCM imaging does not show a clear variation of carrier concentration along the depth of the film. Instead, we observe carrier concentration variations of about 1 order of magnitude (high 1015 to low 1017/cm3) with inhomogeneous spatial regions ranging from a few hundred nm to a few ?m. This nonuniformity is distributed randomly in both the film lateral and vertical directions, independent of grain structure and GBs. We further mapped the surface potential using Kelvin probe force microscopy (KPFM). Higher potential was found on GBs, illustrating positive GB charging but not GB-specific carrier concentration. The results indicate that this suite of techniques can help identify nonuniform carrier concentration and potential fluctuations that can contribute to Voc deficits in GrV-doped CdTe devices.

CdTe thin film photovoltaics↗
DOE OSTI2022

Correlative piezoresponse and micro-Raman imaging of CuInP 2 S 6 –In 4/3 P 2 S 6 flakes unravels phase-specific phononic fingerprint via unsupervised learning

Characterizing the novel properties of layered van der Waals materials is key for their application in functional devices. A better understanding of this type of material requires correlative imaging of diverse nanoscale material properties. Within this class of materials, CuInP 2 S 6 (CIPS) has received a significant degree of interest due to its ionically mediated room temperature ferroelectricity. Moreover, it is possible to form stable self-assembled heterostructures of ferroelectric CuInP 2 S 6 (CIPS) and non-ferroelectric (i.e., lacking Cu) In 4/3 P 2 S 6 (IPS) phases, by controlling the targeted composition and kinetics of synthesis. In this work, we present a correlative nanometric imaging study of the phononic modes and piezoelectricity of the phase-separated thin heteroepitaxial CIPS/IPS flakes. Here, we show that it is possible to isolate the different phononic modes of the two phases by spatially correlating them with their distinct ferroelectric behavior. The coupling of our experimental data with unsupervised learning statistical methods enables unraveling specific Raman peaks that are characteristic of each chemical phase (CIPS and IPS) present in the composite sample, discarding the less significant ones.

correlative microscopy↗
DOE OSTI2020

Statistical properties of homogeneous and isotropic turbulence in He II measured via particle tracking velocimetry

Despite being a quantum two-fluid system, superfluid helium-4 (He II) is observed to behave similarly to classical fluids when a flow is generated by mechanical forcing. This similarity has brought up the feasibility of utilizing He II for high Reynolds number classical turbulence research, considering the small kinematic viscosity of He II. However, it has been suggested that the nonclassical dissipation mechanism in He II at small scales may alter its turbulent statistics and intermittency. In this work, we report our study of a nearly homogeneous and isotropic turbulence (HIT) generated by a towed grid in He II. We measure the velocity field using particle tracking velocimetry with solidified deuterium particles as the tracers. By correlating the velocities measured simultaneously on different particle trajectories or at different times along the same particle trajectory, we are able to conduct both Eulerian and Lagrangian flow analyses. Spatial velocity structure functions obtained through the Eulerian analysis show scaling behaviors in the inertial subrange similar to that for classical HIT but with enhanced intermittency. The Lagrangian analysis allows us to examine the flow statistics down to below the dissipation length scale. Interestingly, abnormal deviations from the classical scaling behaviors are observed in this regime. Lastly, we discuss how these deviations may relate to the motion of quantized vortices in the superfluid component in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2010

High-Speed PLIF Imaging of Hypersonic Transition over Discrete Cylindrical Roughness

In two separate test entries, advanced laser-based instrumentation has been developed and applied to visualize the hypersonic flow over cylindrical protrusions on a flat plate. Upstream of these trips, trace quantities of nitric oxide (NO) were seeded into the boundary layer. The protuberances were sized to force laminar-to-turbulent boundary layer transition. In the first test, a 10-Hz nitric oxide planar laser-induced fluorescence (NO PLIF) flow visualization system was used to provide wide-field-of-view, high-resolution images of the flowfield. The images had sub-microsecond time resolution. However these images, obtained with a time separation of 0.1 sec, were uncorrelated with each other. Fluorescent oil-flow visualizations were also obtained during this test. In the second experiment, a laser and camera system capable of acquiring NO PLIF measurements at 1 million frames per second (1 MHz) was used. This system had lower spatial resolution, and a smaller field of view, but the images were time correlated so that the development of the flow structures could be observed in time.

Danehy, P. M.↗
DOE OSTI2024

Distilling the Evolving Contributions of Anthropogenic Aerosols and Greenhouse Gases to Large‐Scale Low‐Frequency Surface Ocean Changes Over the Past Century

Abstract Anthropogenic aerosols (AER) and greenhouse gases (GHG)—the leading drivers of the forced historical change—produce different large‐scale climate response patterns, with correlations trending from negative to positive over the past century. To understand what caused the time‐evolving comparison between GHG and AER response patterns, we apply a low‐frequency component analysis to historical surface ocean changes from CESM1 single‐forcing large‐ensemble simulations. While GHG response is characterized by its first leading mode, AER response consists of two distinct modes. The first one, featuring long‐term global AER increase and global cooling, opposes GHG response patterns up to the mid‐twentieth century. The second one, featuring multidecadal variations in AER distributions and interhemispheric asymmetric surface ocean changes, appears to reinforce the GHG warming effect over recent decades. AER thus can have both competing and synergistic effects with GHG as their emissions change temporally and spatially.

Dong, Yue↗
DOE OSTI2024

Local-scale heterogeneity of soil thermal dynamics and controlling factors in a discontinuous permafrost region

In permafrost regions, the strong spatial and temporal variability in soil temperature cannot be explained by the weather forcing only. Understanding the local heterogeneity of soil thermal dynamics and their controls is essential to understand how permafrost systems respond to climate change and to develop process-based models or remote sensing products for predicting soil temperature. In this study, we analyzed soil temperature dynamics and their controls in a discontinuous permafrost region on the Seward Peninsula, Alaska. We acquired one-year temperature time series at multiple depths (at 5 or 10 cm intervals up to 85 cm depth) at 45 discrete locations across a 2.3 km 2 watershed. We observed a larger spatial variability in winter temperatures than that in summer temperatures at all depths, with the former controlling most of the spatial variability in mean annual temperatures. We also observed a strong correlation between mean annual ground temperature at a depth of 85 cm and mean annual or winter season ground surface temperature across the 45 locations. We demonstrate that soils classified as cold, intermediate, or warm using hierarchical clustering of full-year temperature data closely match their co-located vegetation (graminoid tundra, dwarf shrub tundra, and tall shrub tundra, respectively). We show that the spatial heterogeneity in soil temperature is primarily driven by spatial heterogeneity in snow cover, which induces variable winter insulation and soil thermal diffusivity. These effects further extend to the subsequent summer by causing variable latent heat exchanges. Finally, we discuss the challenges of predicting soil temperatures from snow depth and vegetation height alone by considering the complexity observed in the field data and reproduced in a model sensitivity analysis.

54 ENVIRONMENTAL SCIENCES↗