Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface parameters”

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

Savannah River Site Climatology for Research Activities Conducted Onsite and in P-Area (1994-2023)

Executive Summary This report presents a climatology of meteorological variables at the Savannah River Site (SRS) using measurements from the P-Area and Central Climatology meteorological towers to support research activities at the Savannah River Site. Key values included in this report are wind, surface pressure, temperature, relative humidity, soil moisture, cloud fraction, cloud base, lightning, atmospheric stability, and boundary layer parameters. The Central Climatology tower is useful in that it is instrumented at four levels, allowing for an examination of height dependence on some variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Reaction Diffusion Modelling of 3D Pillar Electrodes in Single-Catalyst CO 2 Reduction Cascades

Effective electrochemical CO 2 reduction to liquid fuels requires that the local catalytic environment facilitates the desired reactivity, yet a microscopic understanding of this environment is difficult to achieve from experiment alone. In this work, a 3D reaction-diffusion model was developed to explore the effects of electrode surface area and local geometry on the performance of a heterogeneous catalyst that performs a two-step CO 2 reduction cascade reaction to CO and then CH 3 OH under aqueous conditions. Kinetic parameters for the model were inspired by experimental results using a cobalt phthalocyanine (CoPc) catalyst. Three-dimensional architectures composed of arrays of square pillars with varying dimensions and either smooth or periodically modulated surfaces were tested, revealing the extent to which geometry modulates the performance of the cascade reactions. Although structural variations modulate local concentration gradients, we find that electrochemically active surface area predominantly governs the overall cascade reaction. Moreover, the results suggest that supersaturation of CO, with concentrations up to ten-fold higher than the equilibrium solubility limit, might be critical for more efficient conversion to CH 3 OH. For any given geometry, the spatially averaged ratio of [CO] to [CO 2 ] is dictated by the electrochemically active surface area and determines the yield of CH 3 OH. For a fixed surface area, geometries that spatially confine the electrolyte yield moderate local [CO] to [CO 2 ] ratios within small volumes. In contrast, less confining geometries result in a broader distribution of local ratios spread over larger volumes, with both configurations yielding the same spatially averaged [CO] to [CO 2 ] ratio. These insights provide valuable design principles—highlighting the critical importance of surface area and possibly CO supersaturation—for engineering advanced electrode architectures that leverage intermediate trapping and CO supersaturation to enhance overall performance in tandem CO 2 reduction systems.

COMSOL↗
DOE OSTI2026

BNF C-band Scanning ARM Precipitation Radar 2nd Generation (CSAPR-2) Extracted Radar Columns and In-Situ Sensors (RadCLss)

Corrected Moments in Antenna Coordinates (CMAC) calculates quantitative precipitation estimates (QPE) from empirical relationships based on equivalent radar reflectivity factor, specific differ- ential phase, and specific attenuation. To evaluate these empirical relationships, the Extracted Radar Columns and In-Situ Sensors (RADclss) product was developed. Utilizing Py-ART, RAD- clss extracts CMAC radar columns above various ARM and partner locations. These columns are then spatiotemporally synced with in-situ observations at the surface utilizing the Atmospheric data Community Toolkit (ACT; Theisen et al. 2025), allowing direct comparison of radar parameters with rain gauges and laser disdrometers for further investigation.

bankhead↗
DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗
DOE OSTI2024

A Nonintrusive Optical Approach to Characterize Heliostats in Utility-Scale Power Tower Plants: Camera Position Sensitivity Analysis

Optics plays a major role in the effectiveness of concentrating solar power (CSP) technologies. The nonintrusive optical (NIO) approach developed by the National Renewable Energy Laboratory uses uncrewed aircraft system (UAS)-based imaging to survey heliostats in a commercial-scale power tower CSP plant and characterize their optical errors. The image processing algorithm uses photogrammetry to calculate the camera position for each image frame, and the accuracy of the estimated optical errors is highly sensitive to the calculated camera position accuracy. In this study, we simulate a series of case studies in python to examine the impact of different parameters of the sensitivity of the camera calculation, including the number of facet corners used as control points for the photogrammetric calculation, precision error in the detected pixel locations of the facet corners in the image, and precision error of the canting and mounting positions of the facets of the heliostat. The case studies consider heliostat geometry based on three commercial designs to serve as representative examples of different possible sizes of heliostats that the NIO method could be applied to. The results show that increasing the number of control points can improve accuracy for heliostats with many facets, pixel precision has a significantly larger impact on camera calculation accuracy than facet canting and mounting errors, and the camera distance and focal length must be chosen to ensure adequate pixel accuracy on the heliostat surface depending on the size of heliostat. In conclusion, based on the results, recommendations for the allowable values of each parameter are provided to achieve the required NIO optical error estimation accuracy depending on the size of heliostat.

14 SOLAR ENERGY↗
DOE OSTI2025

An automated integrated web-based smart tool for open stope design

The Stability Graph is a widely used tool for the design of open stopes in underground mining. Many users of the Stability Graph still apply this design method manually. Although the manual approach has benefits, using multiple graphs and stability number computation charts for each stope surface is time-consuming, even for the experienced mining engineer. Current practice in the use of the method also limits data sharing. This paper presents a StopeSoft web-based tool for open stope stability prediction that is developed on the basis of the Stability Graph method and is available at openstope.com. StopeSoft incorporates flexibility in terms of Stability Graph options and incorporates additional critical factors often overlooked. As a web-based tool, StopeSoft encourages and makes data sharing possible globally, focused on expanding the database and improving the current limitations of the Stability Graph to provide practical, reliable solutions for mining engineers, consultants, and academics. The StopeSoft automated process facilitates the process of open stope stability prediction, saving time and minimizing potential human errors. Statistical treatment of the data accounts for the variability of input parameters to emphasize the probabilistic nature of the Stability Graph method. The probabilistic interpretation of the stability states of stope surfaces eliminates the false feeling of absolute stope performance based on its location on the Stability Graph , as implied by the deterministic approach.

58 GEOSCIENCES↗
DOE OSTI2024

Dynamic changes in evapotranspiration, canopy photosynthesis and expansive growth in open field under rapid fluctuating radiation

Abstract In open fields, plants experience dynamic changes in environment, particularly radiation, temperature, wind, and humidity but their short-term responses have not been adequately characterized under natural conditions. In this study, we assessed causal effects of rapid radiation fluctuation on seven plant parameters in open fields of cotton and sweet corn canopies. The parameters are evapotranspiration ( ET ), canopy photosynthesis indicated by CO2 flux from air to canopy ( FCO2 ), sensible heat energy flux ( H ), canopy conductance for water vapor ( gcw ), canopy surface temperature ( Ts ), shoot and leaf elongation rate, and stem diameter change. The energy and CO2 fluxes were measured with Bowen ratio/energy balance/CO2 gradient (BREB+) technique, using averaging time of 5 min. Shoot + leaf elongation and stem diameter change were monitored with position transducers using averaging time of 1 min. All parameters were all found to respond to change in radiation and transpiration within minutes or sooner. While radiation effects on canopy gas exchanges are expected, illuminating are the indirect but immediate effects on shoot + leaf growth and stem diameter change through radiation effects on transpiration and plant water status. A novel finding is that gcw also responded within minutes or sooner to radiation fluctuations and that FCO2 was related almost linearly to gcw . Results are discussed in terms of soil-plant-atmosphere continuum, and interpreted in terms of dynamic interactions between transpiration and plant water status. The clear inverse relationship between ET and elongation rate or stem diameter changes provides additional evidence supporting the validity of 5-min averaging for the BREB + technique.

Hsiao, Theodore C.↗
DOE OSTI2025

Extending Badger's rule. I. The relationship between energy and structure in hydrogen bonds

Here, we derive a new expression for the strength of a hydrogen bond (V HB ) in terms of the elongation of the covalent bond of the donor fragment participating in the hydrogen bond (Δr HB ) and the intermolecular coordinates R (separation between the heavy atoms) and θ (deviation of the hydrogen bond from linearity). The expression includes components describing the covalent D–H bond of the hydrogen bond donor via a Morse potential, the Pauli repulsion, and electrostatic interactions between the constituent fragments using a linear expansion of their dipole moment and a quadratic expansion of their polarizability tensor. We fitted the parameters of the model using ab initio electronic structure results for six hydrogen bonded dimers, namely, NH 3 –NH 3 , H 2 O–H 2 O, HF–HF, H 2 O–NH 3 , HF–H 2 O, and HF–NH 3 , and validated its performance for extended parts of their potential energy surfaces, resulting in a mean absolute error ranging from 0.07 to 0.31 kcal/mol. The derived expression describes the energy–structure relationship in terms of a single structural parameter, namely, the elongation of the donor’s covalent bond (Δr HB ), and suggests the novel relationship of 8.0 kcal/mol pm −1 (or 0.8 kcal/mol per 0.001 Å elongation). This structural parameter is easily obtained from theory and can serve as the single descriptor of the strength of individual hydrogen bonds.

Santis, Garrett D. [Univ. of Washington, Seattle, ↗
DOE OSTI2025

A comparison of surrogate constitutive models for viscoplastic creep simulation of HT-9 steel

Mechanistic microstructure-informed constitutive models for the mechanical response of polycrystals are a cornerstone of computational materials science. However, as these models become increasingly more complex – often involving coupled differential equations describing the effect of specific deformation modes – their associated computational costs can become prohibitive, particularly in optimization or uncertainty quantification tasks that require numerous model evaluations. To address this challenge, surrogate constitutive models that balance accuracy and computational efficiency are highly desirable. Data-driven surrogate models, that learn the constitutive relation directly from data, have emerged as a promising solution. In this work, we develop two local surrogate models for the viscoplastic response of a steel: a piecewise response surface method and a mixture of experts model. These surrogates are designed to adapt to complex material behavior, which may vary with material parameters or operating conditions. The surrogate constitutive models are applied to creep simulations of HT-9 steel, an alloy of considerable interest to the nuclear energy sector due to its high tolerance to radiation damage, using training data generated from viscoplastic self-consistent (VPSC) simulations. In conclusion, we define a set of test metrics to numerically assess the accuracy of our surrogate models for predicting viscoplastic material behavior, and show that the mixture of experts model outperforms the piecewise response surface method in terms of accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2024

Updimensioning strategy derived from synthetic equiaxed grain structures for approximating 3D grain size distributions from 2D visualizations with 1D parameters

We generated synthetic equiaxed grain structures using computer graphics software to explore the relationship between various grain size determination methods and true three-dimensional (3D) grain diameters. Mirroring grain measurement techniques, the synthetic 3D grain structures are imaged as 2D micrographs which are measured to yield 1D grain size parameters. Synthetic grain structures provide data at a mass scale and permit exploration of both polished and fractured surface micrographs, revealing one-to-one correspondence between exposed 2D grain cross-sections and individual 3D grains. Analysis of this correspondence yielded a procedure to approximate 3D equiaxed grain size and volume distributions based on the mode of the 2D fractograph grain size distribution. The 3D approximation procedure is shown to be less susceptible to different imaging conditions that affect small, undiscernible grains compared to the standard planimetric and linear intercept methods, which by design also tend to underestimate the 3D grain diameter. The procedure requires larger sample sizes to lower variance and a deeper analysis which could become more practical with machine learning (ML) models for grain boundary segmentation, which synthetic grain structures can help train. This work lays the foundation for analyzing other grain distributions such as columnar and composite grains in similar depth.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

High-efficiency, high-fidelity charge initialization of shallow nitrogen-vacancy centers in diamond

Nitrogen-vacancy (N-𝑉) centers in diamond exhibit long spin-coherence times, optical initialization, and optical-spin readout under ambient conditions, making them excellent quantum sensors. However, the conventional scheme for charge-state initialization based on off-resonant green excitation results in significant state-preparation errors, typically around 30%. One method for improving charge-state initialization fidelity is to use multicolor excitation, which has been demonstrated to achieve a near-unity preparation fidelity for bulk N-𝑉 centers by using a few milliseconds of near-infrared (NIR) (5-mW) and green (10-μ⁢W) excitation. The translation of such schemes to N-𝑉 centers near the diamond surface with higher-efficiency optical pumping would enable new applications in nanoscale sensing. Here, we demonstrate a protocol for efficient charge initialization of shallow N-𝑉 centers between 5 nm and 15 nm from the diamond surface. By carefully studying the charge dynamics of shallow N-𝑉 centers, we identify a region of parameter space that allows for near-unity (95%) charge initialization within 300 μ⁢s of NIR (905-nm, 1-mW) and green (520-nm, 10-μ⁢W) excitation. The time to 90% charge initialization can be as fast as 10 μ⁢s for 4 mW of NIR and 39 μ⁢W of green illumination. This fast, efficient charge initialization protocol will especially benefit nanoscale sensing applications in which state-preparation errors currently prohibit scaling, such as measuring higher-order multipoint correlators.

infrared techniques↗
DOE OSTI2024

Contactless Production Testing of Silicon Solar Cells

Critical cost reductions in silicon solar cells are made by minimizing silver usage in their metal contacts. The reason why is both clear and unavoidable—silver is not only expensive, but also a potential limiting resource for scale-up (the PV industry alone used 10.3 % of the global silver supply in 2020 [1]). The result is solar cells with vanishing electrical contact area that are extremely difficult to test prior to module manufacture, when they must be sorted for quality to maximize power output and reliability of the subsequent modules. This project developed a new measurement instrument that enables the transition to solar cells with very low silver content, including cell designs with both minimal busbars and no busbars at all to electrically contact. The resulting instrument at project completion demonstrated the ability to sort cells with comparable—or better—results than existing technology. Compared to existing cell-test instruments, the new tool has minimal contacting requirements, which lowers the maintenance costs and use of consumable parts. We demonstrated innovative, yet pragmatic, solutions for reporting a comprehensive set of measurement results useful for both cell sorting (going forward in the line for module manufacture) and process control (looking backward in the line to identify cell manufacturing issues from wafer to cell-test). In addition to the traditional current-voltage characterization at cell test—short circuit current, open-circuit voltage, power, efficiency and fill factor—the new tool maintains all the advanced parameter characterization of our existing product line such as a substrate doping measurement, carrier recombination (lifetime) analysis, and surface recombination analysis.

14 SOLAR ENERGY↗
DOE OSTI2024

Thermographic phosphor digital image correlation

The invention advances the ratio method of phosphor thermography using two machine vision cameras for full-field temperature measurements of a solid surface. As an example, algorithms from digital image correlation (DIC) can be used to determine the stereoscopic imaging system intrinsic and extrinsic parameters, and accurately register material points on the sample to subpixel locations in each image with 0.07 px or better accuracy. A phosphor calibration sample fabricated using aerosol deposition can be used for in situ determination of the temperature-versus-intensity ratio relationship. The full calibration methodology and several improvements on two-color phosphor thermography open the door for full-field temperature measurements in dynamic tests with deforming test specimens. In particular, thermographic phosphors can be combined with stereo digital image correlation in a novel diagnostic, TP+DIC, to measure time-resolved full-field surface strains and temperatures simultaneously, thereby enabling testing in combined thermo-mechanical environments.

Winters, Caroline↗
DOE OSTI2025

Observation of Vortex Stripes in UTe 2

In conventional type-II superconductors, quantum vortices determine the magnetic response of the system and tend to form regular lattices. UTe 2 is a recently discovered heavy Fermion superconductor exhibiting many anomalous macroscopic behaviors, but whether it has a multicomponent order parameter remains open. Here, in this work, we study the vortices of UTe 2 by employing scanning superconducting quantum interference device microscopy. While a small out-of-plane magnetic field produces typical isolated vortices, a higher field generates vortex stripe patterns that evolve with the vortex density. The stripes form at different locations and along different directions in the surface plane when the field is along the b or c axis. The behavior is reproduced by our simulation based on an anisotropic two-component order parameter. This study shows that UTe 2 has a nontrivial disparity of multiple length scales, placing constraints on the multicomponent superconductivity.

UTe2↗
DOE OSTI2024

Tuning Surface Adhesion Using Grayscale Electron-beam Lithography

Surface texturing of manufactured products tailors their properties, such as friction, adhesion, biocompatibility, or fluid interactions. However, advancements in this area are largely the result of trial-and-effort testing and generally lack a science-guided framework for determining the surface topography that will optimize performance. The present investigation explores grayscale electron-beam lithography as a means to create multiscale surface patterns to control surface performance. Here, we created and characterized a set of surface textures on a silicon wafer; the textures were superpositions of sine waves of varying wavelengths and amplitudes. First, the multiscale topography of the patterned surface was characterized, using profilometry and atomic force microscopy, to understand its fidelity to the designed-in pattern. The results of this analysis demonstrated how grayscale lithography accurately controlled the lateral size of features but was less precise on the vertical height of the surface, and also introduced inherent roughness below the scale of patterning. Second, a micromechanical tester was used to characterize the adhesion of the surfaces with large-scale polished silicon spheres. The results showed that adhesion could be tailored, with significant contribution from all of the designed-in length scales of topography. The strength of adhesion did not correlate with conventional roughness parameters but could be accurately modeled using simple numerical integration. Taken together, this investigation demonstrates the promise and challenges of grayscale e-beam lithography with multiscale patterns as a method for the tailoring of surface performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Adsorption↗
DOE OSTI2024

Soil Carbon and Nitrogen Stocks Across Hillslopes Underlain by Continuous Permafrost in the Northern Arctic Foothills, Alaska, United States

Constraining the variability of soil organic carbon (SOC) and total nitrogen (TN) stocks across hillslopes in Low Arctic permafrost-affected landscapes remains a significant challenge for improving global estimates of permafrost SOC stocks. We investigated SOC and TN stocks across hillslopes at two sites in the Arctic Foothills of Alaska, United States (Happy Valley and Sagwon Hills). Average SOC and TN stocks for the 0–1-m depth interval were high (52.0 ± 15.1 kg C m −2 and 2.74 ± 0.82 kg N m −2 ) and linearly related ( R 2 = 0.74, p < 0.0001). Unlike soils of other permafrost and nonpermafrost landscapes, variability was greatest within rather than between hillslope positions. Furthermore, SOC and TN stocks in the surface 1 m did not exhibit strong patterns by hillslope position and were only weakly associated with major geomorphic parameters that typically predict SOC and TN stocks well in other landscapes. Although sampling at upper hillslope positions was largely limited to depths of less than 1.5 m due to the presence of coarse fragments in reworked glacial till, deeper observations at lower hillslope positions (footslopes, toeslopes, and basins) revealed significantly larger SOC stocks (92.0 ± 18.0 kg C m −2 at 2 m; 117.1 ± 10.4 kg C m −2 at 3 m). The unique small-scale variability in ice content, cryoturbation, patterned ground, and organic layer thickness on these broad, Low Arctic sites contributes to the relatively homogeneous distribution of SOC and TN stocks across hillslope positions in the top 1 m, but a future focus on deeper sampling may reveal greater differences in SOC and TN stocks.

54 ENVIRONMENTAL SCIENCES↗