Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Location”

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 703 records · Page 39

Global Life Cycle and Techno-Economic Assessment of Algal-Based Biofuels

Techno-economic analyses (TEAs) and life cycle assessments (LCAs) of algal biofuels often focus on locations in suboptimal latitudes for algal cultivation, which can under-represent the sustainability potential of the technology. This study identifies the optimal global productivity potential, environmental impacts, and economic viability of algal biofuels by using validated biophysical and sustainability modeling. The biophysical model simulates growth rates of Scenedesmus obliquusbased on temperature, photoinhibition, and respiration effects at 6685 global locations. Region-specific labor costs, construction factors, and tax rates allow for spatially resolved TEA, while the LCA includes regional impacts of electricity, hydrogen, and nutrient markets across ten environmental categories. The analysis identifies optimal locations for algal biofuel production in terms of environmental impacts and economic viability which are shown to follow biomass yields. Modeling results highlight the global variability of productivity with maximum yields ranging between 24.8 and 27.5 g m -2 d -1 in equatorial regions. Environmental impact results show favorable locations tracked with low-carbon electricity grids, with the well-to-wheels global warming potential (GWP) ranging from 31 to 45 g CO 2eq MJ -1 in South America and Central Africa. When including direct land use change impacts, the GWP ranged between 44 and 55 g CO 2eq MJ -1 in these high-productivity regions. Low-carbon electricity also favors air quality and eutrophication impacts. The TEA shows that minimum algal fuel prices of $\$1.89$-$\$2.15$ per liter of gasoline-equivalent are possible in southeast Asia and Venezuela. Furthermore, this discussion focuses on the challenges and opportunities to reduce fuel prices and the environmental impacts of algal biofuels in various global regions.

09 BIOMASS FUELS↗

Unlocking Solutions: Innovative Approaches to Identifying and Mitigating the Environmental Impacts of Undocumented Orphan Wells in the United States

In the United States, hundreds of thousands of undocumented orphan wells have been abandoned, leaving the burden of managing environmental hazards to governmental agencies or the public. These wells, a result of over a century of fossil fuel extraction without adequate regulation, lack basic information like location and depth, emit greenhouse gases, and leak toxic substances into groundwater. For most of these wells, basic information such as well location and depth is unknown or unverified. Addressing this issue necessitates innovative and interdisciplinary approaches for locating, characterizing, and mitigating their environmental impacts. Our survey of the United States revealed the need for tools to identify well locations and assess conditions, prompting the development of technologies including machine learning to automatically extract information from old records (95%+ accuracy), remote sensing technologies like aero-magnetometers to find buried wells, and cost-effective methods for estimating methane emissions. Notably, fixed-wing drones equipped with magnetometers have emerged as cost-effective and efficient for discovering unknown wells, offering advantages over helicopters and quadcopters. Efforts also involved leveraging local knowledge through outreach to state and tribal governments as well as citizen science initiatives. These initiatives aim to significantly contribute to environmental sustainability by reducing greenhouse gases and improving air and water quality.

54 ENVIRONMENTAL SCIENCES↗

Spectroscopic Online Monitoring: Using a Multi-Track Visible Spectrometer to Facilitate a Mass Balance Study in a Simulated TALSPEAK Process

Nuclear energy is a promising low-carbon energy candidate to meet the increased demand for green energy, where the integration of fuel recycling can have significant benefits for material usage and waste reduction. Utilizing in situ monitoring tools can provide ample opportunities to better control and safeguard nuclear material recycle processes while also offering knowledge and insight into real-time solution properties. The simultaneous measurement of analytical targets in multiple process locations can enable real-time mass balance and material accountancy calculations. This is demonstrated here with a mass balance study of Nd 3+ on countercurrent aqueous/organic metal extraction within a single centrifugal contactor. The Nd 3+ concentration was simultaneously monitored at the inlets and outlets of both aqueous and organic phases using a visible absorbance detector that allowed for the simultaneous measurement of up to six locations. The Nd 3+ concentration was calculated by using chemical data science algorithms, where model training sets were collected on a single track of the detector. The discussion includes addressing the challenges of using a model collected on a single track and applying it as a model across the other tracks on the detector. Each track of the detector corresponds to one measurement location on the contactor. The difference in the integrated moles of Nd 3+ between the inlet and outlet at the end of the experiment was near zero, indicating that the mass balance of this experiment was maintained. Overall, the online spectroscopic monitoring was able to follow changing solution conditions and accurately measure the concentration of Nd 3+ in different locations within the contactor system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

macroMS: Image-Guided Analysis of Random Objects by Matrix-Assisted Laser Desorption/Ionization Time-of-Flight Mass Spectrometry

Mass spectrometry imaging is well-suited to characterizing sample surfaces for their chemical content in a spatially resolved manner. However, when the surface contains small objects with significant empty spaces between them, more efficient approaches to sample acquisition are possible. Image-guided mass spectrometry (MS) enables high-throughput analysis of a diverse range of sample types, such as microbial colonies, liquid microdroplets, and others, by recognizing and analyzing selected location targets in an image. Here, we describe an imaging protocol and macroMS, an online software suite that can be used to enhance MS measurements of macroscopic samples that are imaged by a camera or a flatbed scanner. Furthermore, the web-based tool enables users to find and filter targets from the optical images, correct optical distortion issues for improved spatial location of selected targets, input the custom geometry files into an MS device to acquire spectra at the selected locations, and finally, perform limited data analysis and use visualization tools to aid locating samples containing compounds of interest. Using the macroMS suite, an enzyme mutant library of Saccharomyces cerevisiae and nL droplet arrays of Escherichia coli and Pseudomonas fluorescens have been assayed at a rate of ~2 s/sample.

47 OTHER INSTRUMENTATION↗

River Plume Liftoff Dynamics and Surface Expressions

The water surface expression of liftoff and its dependence on discharge are examined using numerical simulations with the Regional Ocean Modeling System (ROMS). Liftoff is the process by which buoyant river water separates from the bed and flows over denser saltwater. During low-discharge conditions liftoff occurs in the river and is accompanied by a change in the surface slope. During high-discharge conditions liftoff occurs outside the mouth and generates a ridge on the water surface. The location and height of the ridge can be described by analytical equations in terms of discharge, shelf slope, and river mouth aspect ratio. The offshore distance and height of the ridge are proportional to the river discharge and vary inversely with river mouth aspect ratio. For steep shelf slopes liftoff occurs close to the river mouth and generates a large ridge. The ridge is modified, but not eliminated, by the presence of tides. The water surface slope change at the ridge peak is large enough to be detected by the upcoming Surface Water and Ocean Topography (SWOT) altimeter and can be used to identify the liftoff location during high discharge. However, during low discharge the water surface slope change at the liftoff location is too small to be detected by SWOT. These results indicate that remote measurements of the presence or absence of the ridge may be useful to distinguish between low and high flows, and remote measurements of the ridge location or height could be used to estimate freshwater discharge.

42 ENGINEERING↗

Exploring the Constraints on Simulated Aerosol Sources and Transport across the North Atlantic with Island-Based Sun Photometers

Atmospheric aerosol over the North Atlantic Ocean impact regional clouds and climate. In this work, we use a set of sun photometer observations of Aerosol Optical Depth (AOD) located on the Graciosa and Cape Verde islands, along with the GEOS-Chem chemical transport model to investigate the sources of these aerosol and their transport over the North Atlantic Ocean. At both locations, the largest simulated contributor to aerosol extinction is the local source of sea-salt aerosol. In addition to this large source, we find that signatures consistent with long-range transport of anthropogenic, biomass burning, and dust emissions are apparent throughout the year at both locations. Model simulations suggest that this signal of long-range transport in AOD is more apparent at higher elevation locations; the influence of anthropogenic and biomass burning aerosol extinction is particularly pronounced at the height of Pico Mountain, near the Graciosa Island site. Using a machine learning approach, we further show that simulated observations at these three sites (near Graciosa, Pico Mountain, and Cape Verde) can be used to predict the simulated background aerosol imported into cities on the European mainland, particularly during the local winter months, highlighting the utility of background AOD monitoring for understanding downwind air quality.

54 ENVIRONMENTAL SCIENCES↗

Continental-Scale Geographic Trends in Barometric-Pumping Efficiency Potential: A North American Case Study

Barometric pumping is a gas transport mechanism that has important implications for many applications involving subsurface gas seepage processes. Here, this study provides the first continental-scale analysis of barometric-pumping efficiency potential based on meteorology. We quantified the barometric-pumping efficiency potential at 1,257 locations across the continental US and Canada. The results provide continental-scale geographic dependencies of barometric-pumping efficiency potential, indicating a significant correlation with latitude and a nonlinear dependence on longitude. The analysis also indicates that variability in barometric-pumping efficiency potential decreases with distance from the coast and as elevation increases. Locations far from the coastline are more likely to have upper mid-range potentials, while higher elevation locations are more likely to have low potentials. The highest barometric-pumping efficiency potentials are mostly found around the Gulf of St. Lawrence around 50°N. Locations along the Atlantic coast exhibit large-scale variations in potentials with a clear increasing trend with latitude.

58 GEOSCIENCES↗

Probing the Magnetosheath Boundaries Using Interstellar Boundary Explorer (IBEX) Orbital Encounters

Inside the magnetosheath, the IBEX-Hi energetic neutral atom (ENA) imager measures a distinct background count rate that is more than 10 times the typical heliospheric ENA emissions observed when IBEX is outside the magnetosheath. The source of this enhancement is magnetosheath ions of solar wind (SW) origin that deflect around the Earth's magnetopause (MP), scatter and neutralize from the anti-sunward part of the IBEX-Hi sunshade, and continue into the instrument as neutral atoms, behaving indistinguishably from ENAs emitted from distant plasma sources. While this background pollutes observations of outer heliospheric ENAs, it provides a clear signature of IBEX crossings over the magnetospheric boundaries. In this study, we investigate IBEX encounters with the magnetosheath boundaries using ~8 yr of orbital data, and we determine the MP and bow shock (BS) locations derived from this background signal. We find 280 BS crossings from X GSE ~ 11 R e to X GSE ~ –36 R e and 241 MP crossings from X GSE ~ 6 R e to X GSE ~ –48 R e . We compare IBEX BS and MP crossing locations to those from IMP-8, Geotail, Cluster, Magion-4, ISEE, and Magnetospheric Multiscale Mission, and we find that IBEX crossing locations overlap with the BS and MP locations inferred from these other data sets. In this paper, we demonstrate how IBEX can be used to identify magnetosheath crossings, and extend boundary observations well past the terminator, thus further constraining future models of magnetosheath boundaries. Furthermore, we use the IBEX data set to show observational evidence of near-Earth magnetotail squeezing during periods of strong interplanetary magnetic field B y .

79 ASTRONOMY AND ASTROPHYSICS↗

Spatiotemporal Graph Convolutional Networks for Earthquake Source Characterization

Abstract Accurate earthquake location and magnitude estimation play critical roles in seismology. Recent deep learning frameworks have produced encouraging results on various seismological tasks (e.g., earthquake detection, phase picking, seismic classification, and earthquake early warning). Many existing machine learning earthquake location methods utilize waveform information from a single station. However, multiple stations contain more complete information for earthquake source characterization. Inspired by recent successes in applying graph neural networks (GNNs) in graph‐structured data, we develop a Spatiotemporal Graph Neural Network (STGNN) for estimating earthquake locations and magnitudes. Our graph neural network leverages geographical and waveform information from multiple stations to construct graphs automatically and dynamically by adaptive message passing based on graphs' edges. Using a recent graph neural network and a fully convolutional neural network as baselines, we apply STGNN to earthquakes recorded by the Southern California Seismic Network from 2000 to 2019 and earthquakes collected in Oklahoma from 2014 to 2015. STGNN yields more accurate earthquake locations than those obtained by the baseline models and performs comparably in terms of depth and magnitude prediction, though the ability to predict depth and magnitude remains weak for all tested models. Our work demonstrates the potential of using GNNs and multiple stations for better automatic estimation of earthquake epicenters.

58 GEOSCIENCES↗

Characterizing Mesoscale Cellular Convection in Marine Cold Air Outbreaks With a Machine Learning Approach

Abstract During marine cold‐air outbreaks (MCAOs), when cold polar air moves over warmer ocean, a well‐recognized cloud pattern develops, with open or closed mesoscale cellular convection (MCC) at larger fetch over open water. The Cold‐Air Outbreaks in the Marine Boundary Layer Experiment provided a comprehensive set of ground‐based in situ and remote sensing observations of MCAOs at a coastal location in northern Norway. MCAO periods that unambiguously exhibit open or closed MCC are determined. Individual cells observed with a profiling Ka‐band radar are identified using a watershed segmentation method. Using self‐organizing maps (SOMs), these cells are then objectively classified based on the variability in their vertical structure. The SOM nodes contain some information about the location of the cell transect relative to the center of the MCC. This adds classification noise, requiring numerous cell transects to isolate cell dynamical information. The SOM‐based classification shows that comparatively intense convection occurs only in open MCC. This convection undergoes an apparent lifecycle. Developing cells are associated with stronger updrafts, large spectrum width, larger amounts of liquid water, lower surface precipitation rates, and lower cloud tops than mature and weakening cells. The weakening of these cells is associated with the development of precipitation‐induced cold pools. The SOM classification also reveals less intense convection, with a similar lifecycle. More stratiform vertical cloud structures with weak vertical motions are common during closed MCC periods and are separated into precipitating and non‐precipitating stratiform cores. Convection is observed only occasionally in the closed MCC environment.

Meteorology & Atmospheric Sciences↗

Aerial imagery dataset of lost oil wells

Orphaned wells are wells for which the operator is unknown or insolvent. The location of hundreds of thousands of these wells remain unknown in the United States alone. Cost-effective techniques are essential to locate orphaned wells to address environmental problems. In this paper, we present a dataset consisting of 120,948 aerial images of recently documented orphan wells. Each of these 512 × 512 images is paired with segmentation masks that indicate the presence or absence of such well. These images, sourced from the National Agriculture Imagery Program, cover the continental United States with spatial resolutions ranging from 30 centimeters to 1 meter. Additionally, we included negative examples by selecting locations uniformly across the United States. Accompanying metadata includes the IDs and spatial resolution of the original images, which are available for free through the United States Geological Survey, and the pixel coordinates of documented orphaned wells identified in these images. This dataset is intended to support the development of deep-learning models that can help locating undocumented orphan wells from such imagery, thereby blunting the environmental damage they do.

Climate-change mitigation↗

Analysis of solid uranium particulates on cotton swipes with an automated microextraction-ICP-MS system

Here, an automated microextraction method coupled to an inductively coupled plasma – mass spectrometer (ICP-MS) was developed for the direct analysis of solid uranium particulates on the surface of cotton swipes. The microextraction probe extracts particulates from the sample surface, in a flowing solvent, and directs the removed analyte to an ICP-MS for isotopic determination. The automated system utilizes a mechanical XY stage that is software controlled with the capability of saving and returning to specific locations and a camera focused to the swipe surface for optimal viewing of the extracted locations (i.e., material present). Here, particulates (n = 135) were extracted and measured by ICP-MS, including 35 depleted uranyl nitrate hexahydrate (UN) (used for mass bias corrections), 50 uranyl fluoride (UO 2 F 2 ), and 50 uranyl acetate (UAc) particulates. Blank extractions were performed on the cotton swipes between triplicate sample analyses. Between each swipe extraction, the probe was sent between two wells containing 10% and 5% HNO 3 to clean the probe head and to eliminate any analyte carryover between particulates. The measured 235 U/ 238 U and 234 U/ 238 U isotope ratios for the UO 2 F 2 particulates were 0.00725(8) and 0.000054(4), a percent relative difference (% RD) of –0.041% and –1.7% from the reference isotope ratios determined in-lab through multi-collector ICP-MS analysis of dissolved aliquots of the U material. The UAc samples had a measured 235 U/ 238 U isotope ratio of 0.00206(7), a –0.96% relative difference from the reference value of 0.00208(1). The 234 U/ 238 U and 236 U/ 238 U isotope ratios were 0.000008(1) and 0.000031(4), –5.1% RD and –4.3% RD, respectively. The automated sample stage enabled seamless and rapid particle analysis, leading to a significant increase in throughput versus what was previously possible. Additionally, the saved location capability reduced user sampling error as sampling locations were easily stored and recalled. Analysis of U particles on the swipe surface – including blanks, mass bias, and triplicate extractions – was completed in less than an hour without any sample preparation necessary.

47 OTHER INSTRUMENTATION↗

Explaining an unusual electromigration behavior—A comprehensive experimental and theoretical analysis using finite element method

In metallic interconnects, it is generally assumed that electromigration (EM) failure location is independent of the applied electrical current and always occurs at the highest-current-density area. Our experiments show otherwise. We designed an Al interconnect that alters its failure location by only varying the applied current density. The failure occurs near the high for a current above 2 × 10 7 A/cm 2 , but at a location with 59% of the maximum for lower current densities. Thermoreflectance thermal imaging is employed to gather time-dependent high-resolution spatial temperature distributions of the Al interconnect during EM. More importantly, we propose a computationally inexpensive 2D finite element method that tracks EM evolution in time and matches well with the observations from different experimental conditions. A detailed analysis covering the major driving forces of EM is carried out to understand the complex physics behind EM. The atomic depletion rate contributed by each force is quantitatively studied. By examining the results from every tested experimental condition, the model reveals that the temperature gradient is the key reason causing atomic depletion near the failure location. Graphical illustrations and qualitative analysis are provided to intuitively show the key findings of our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Field validation of an actuated convergent-beam acoustic Doppler profiler for high resolution flow mapping

Abstract Energetic river and tidal flow environments feature complex hydrodynamic conditions. Conventional acoustic Doppler profiling instrumentation typically requires assumptions of flow homogeneity over the spatial scales of the divergent beam separation. This removes the ability to measure spatio-temporal variability within the flow. However, velocity variability within these spatial scales is often important in the design of structures and devices exposed to such flow, informing dynamic and peak hydrodynamic load predictions. The research presented outlines the development and testing of a flow measurement instrument consisting of multiple spatially-separated single beam acoustic Doppler profilers converging on a remote focal point, the location of which can be programmatically adjusted through actuation. This increases the spatial resolution at which remote field measurements can be made in energetic flow environments. Field testing of the instrument was conducted in a tidal channel at Sequim Bay Inlet, WA, USA. Results are compared with independent reference measurements made by a co-located, motion-corrected acoustic Doppler velocimeter. This comparison, across tidal velocities of 0.4–0.7 m s −1 , showed a mean velocity error of 0.5%–13.2% across nine sample locations within a 3.00 m × 2.25 m plane. To the authors’ knowledge, this is the first publication of flow velocity results from non-intrusive acoustic measurements at an off-axis, focal point location that features 3D positional control.

13 HYDRO ENERGY↗

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neural networks for estimation of divertor conditions in DIII-D using C III imaging

Deep learning approaches have been applied to images of C III emission in the lower divertor of DIII-D to develop models for estimating the level of detachment and magnetic configuration (X-point location and strike point radial location). The poloidal distance from the target to the C III emission front is used to represent the level of detachment. The models perform well on a test dataset not used in training, achieving $F_1$ scores as high as 0.99 for detachment state classification and root mean squared error (RMSE) as low as 2cm for front location regression. Predictions for shots with intermittent reattachment are studied, with class activation mapping used to aid in interpretation of the model predictions. Based on the success of these models, a third model was trained to predict the X-point location and strike point radial position from C III images. Though the dataset covers only a small range of possible magnetic configurations, the model shows promising results, achieving RMSE around 1cm for the test data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Canopy cover and microtopography control precipitation-enhanced thaw of ecosystem-protected permafrost

Abstract Northern high-latitudes are projected to get warmer and wetter, which will affect rates of permafrost thaw and mechanisms by which thaw occurs. To better understand the impact of rain, as well as other factors such as snow depth, canopy cover, and microtopography, we instrumented a degrading permafrost plateau in south-central Alaska with high-resolution soil temperature sensors. The site contains ecosystem-protected permafrost, which persists in unfavorable climates due to favorable ecologic conditions. Our study (2020–2022) captured three of the snowiest years and three of the four wettest years since the site was first studied in 2015. Average thaw rates along an across-site transect increased nine-fold from 6 ± 5 cm yr −1 (2015–2020) to 56 ± 12 cm yr −1 (2020–2022). This thaw was not uniform. Hummock locations, residing on topographic high points with relatively dense canopy, experienced only 8 ± 9 cm yr −1 of thaw, on average. Hollows, topographic low points with low canopy cover, and transition locations, which had canopy cover and elevation between hummocks and hollows, thawed 44 ± 6 cm yr −1 and 39 ± 13 cm yr −1 , respectively. Mechanisms of thaw differed between these locations. Hollows had high warm-season soil moisture, which increased thermal conductivity, and deep cold-season snow coverage, which insulated soil. Transition locations thawed primarily due to thermal energy transported through subsurface taliks during individual rain events. Most increases in depth to permafrost occurred below the ∼45 cm thickness seasonally frozen layer, and therefore, expanded existing site taliks. Results highlight the importance of canopy cover and microtopography in controlling soil thermal inputs, the ability of subsurface runoff from individual rain events to trigger warming and thaw, and the acceleration of thaw caused by consecutive wet and snowy years. As northern high-latitudes become warmer and wetter, and weather events become more extreme, the importance of these controls on soil warming and thaw is likely to increase.

54 ENVIRONMENTAL SCIENCES↗