USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Oceans”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 1,693 records · Page 94Linked to original sources
USGS2003

Evidence for a Battle Mountain-Eureka crustal fault zone, north-central Nevada, and its relation to Neoproterozoic-Early Paleozoic continental breakup

Combined evidence from gravity, radiogenic isotope, and magnetotelluric (MT) data indicates a crustal fault zone that coincides with the northwest-trending Battle Mountain-Eureka (BME) mineral trend in north-central Nevada, USA. The BME crustal fault zone likely originated during Neoproterozoic-Early Paleozoic rifting of the continent and had a large influence on subsequent tectonic events, such as emplacement of allochthons and episodic deformation, magmatism, and mineralization throughout the Phanerozoic. MT models show the fault zone is about 10 km wide, 130-km long, and extends from 1 to 5 km below the surface to deep crustal levels. Isotope data and gravity models imply the fault zone separates crust of fundamentally different character. Geophysical evidence for such a long-lived structure, likely inherited from continental breakup, defies conventional wisdom that structures this old have been destroyed by Cenozoic extensional processes. Moreover, the coincidence with the alignment of mineral deposits supports the assertion by many economic geologists that these alignments are indicators of buried regional structures.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Slip rate and earthquake recurrence along the central Septentrional fault, North American-Caribbean plate boundary, Dominican Republic

The Septentrional fault zone (SFZ) is the major North American-Caribbean, strike-slip, plate boundary fault at the longitude of eastern Hispaniola. The SFZ traverses the densely populated Cibao Valley of the Dominican Republic, forming a prominent scarp in alluvium. Our studies at four sites along the central SFZ are aimed at quantifying the late Quaternary behavior of this structure to better understand the seismic hazard it represents for the northeastern Caribbean. Our investigations of excavations at sites near Rio Cenovi show that the most recent ground-rupturing earthquake along this fault in the north central Dominican Republic occurred between A.D. 1040 and A.D. 1230, and involved a minimum of ???4 m of left-lateral slip and 2.3 m of normal dip slip at that site. Our studies of offset stream terraces at two locations, Rio Juan Lopez and Rio Licey, provide late Holocene slip rate estimates of 6-9 mm/yr and a maximum of 11-12 mm/yr, respectively, across the Septentrional fault. Combining these results gives a best estimate of 6-12 mm/yr for the slip rate across the SFZ. Three excavations, two near Tenares and one at the Rio Licey site, yielded evidence for the occurrence of earlier prehistoric earthquakes. Dates of strata associated with the penultimate event suggest that it occurred post-A.D. 30, giving a recurrence interval of 800-1200 years. These studies indicate that the SFZ has likely accumulated elastic strain sufficient to generate a major earthquake during the more than 800 years since it last slipped and should be considered likely to produce a destructive future earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Reservoir depletion at The Geysers geothermal area, California, shown by four-dimensional seismic tomography

Intensive geothermal exploitation at The Geysers geothermal area, California, induces myriads of small-magnitude earthquakes that are monitored by a dense, permanent, local seismometer network. Using this network, tomographic inversions were performed for the three-dimensional Vp and Vp/Vs structure of the reservoir for April 1991, February 1993, December 1994, October 1996, and August 1998. The extensive low-Vp/Vs anomaly that occupies the reservoir grew in strength from a maximum of 9% to a maximum of 13.4% during the 7-year study period. This is attributed to depletion of pore liquid water in the reservoir and replacement with steam. This decreases Vp by increasing compressibility, and increases Vs because of reduction in pore pressure and the drying of argillaceous minerals, e.g., illite, which increase the shear modulus. These effects serendipitously combine to lower Vp/Vs, resulting in a strong overall effect that provides a convenient tool for monitoring reservoir depletion. Variations in the Vp and Vs fields indicate that water depletion is the dominant process in the central part of the exploited reservoir, and pressure reduction and mineral drying in the northwest and southeast parts of the reservoir. The rate at which the Vp/Vs anomaly grew in strength in the period 1991-1998 suggests most of the original anomaly was caused by exploitation. Continuous monitoring of Vp, Vs, and Vp/Vs is an effective geothermal reservoir depletion monitoring tool and can potentially provide information about depletion in parts of the reservoir that have not been drilled.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Three-dimensional crustal structure of Long Valley caldera, California, and evidence for the migration of CO2 under Mammoth Mountain

A temporary network of 69 three-component seismic stations captured a major seismic sequence in Long Valley caldera in 1997. We performed a tomographic inversion for crustal structure beneath a 28 km ?? 16 km area encompassing part of the resurgent dome, the south moat, and Mammoth Mountain. Resolution of crustal structure beneath the center of the study volume was good down to ???3 km below sea level (???5 km below the surface). Relatively high wave speeds are associated with the Bishop Tuff and lower wave speeds characterize debris in the surrounding moat. A low-Vp/Vs anomaly extending from near the surface to ???1 km below sea level beneath Mammoth Mountain may represent a CO2 reservoir that is supplying CO2-rich springs, venting at the surface, and killing trees. We investigated temporal variations in structure beneath Mammoth Mountain by differencing our results with tomographic images obtained using data from 1989/1990. Significant changes in both Vp and Vs were consistent with the migration of CO2 into the upper 2 km or so beneath Mammoth Mountain and its depletion in peripheral volumes that correlate with surface venting areas. Repeat tomography is capable of detecting the migration of gas beneath active silicic volcanoes and may thus provide a useful volcano monitoring tool.

Journal of Geophysical Research B: Solid Earth↗
USGS2014

An empirical approach to modeling methylmercury concentrations in an Adirondack stream watershed

Inverse empirical models can inform and improve more complex process-based models by quantifying the principal factors that control water quality variation. Here we developed a multiple regression model that explains 81% of the variation in filtered methylmercury (FMeHg) concentrations in Fishing Brook, a fourth-order stream in the Adirondack Mountains, New York, a known “hot spot” of Hg bioaccumulation. This model builds on previous observations that wetland-dominated riparian areas are the principal source of MeHg to this stream and were based on 43 samples collected during a 33 month period in 2007–2009. Explanatory variables include those that represent the effects of water temperature, streamflow, and modeled riparian water table depth on seasonal and annual patterns of FMeHg concentrations. An additional variable represents the effects of an upstream pond on decreasing FMeHg concentrations. Model results suggest that temperature-driven effects on net Hg methylation rates are the principal control on annual FMeHg concentration patterns. Additionally, streamflow dilutes FMeHg concentrations during the cold dormant season. The model further indicates that depth and persistence of the riparian water table as simulated by TOPMODEL are dominant controls on FMeHg concentration patterns during the warm growing season, especially evident when concentrations during the dry summer of 2007 were less than half of those in the wetter summers of 2008 and 2009. This modeling approach may help identify the principal factors that control variation in surface water FMeHg concentrations in other settings, which can guide the appropriate application of process-based models.

New York↗
USGS1991

Fractal hierarchies of magma transport in Hawaii and critical self- organization of tremor

A hierarchical model of magma transport in Hawaii is developed from the seismic records of deep (30-60 km) and intermediate-depth (5-15 km) harmonic tremor between January 1, 1962, and December 31, 1983. A tremor model of magma transport is developed from mass balances of percolation that are proportional to tremor durations. It gives reasonable magma fractions and residence times for a vertical drift velocity of 4 km yr -1 and yields patterns of intermittency that are in accord with singularity analyses of the 22-year time series record. It is suggested that spatiotemporal universality extends from small to large scales in Hawaiian and other magmatic systems. The apparent universal scaling of frequencies may be more than 15 decades in time (1 s to about 60 m.y.) and 10 decades in length (0.1 mm to 10 3 km). -from Authors

Hawaii↗
USGS1999

Carbon cycling in boreal wetlands: A comparison of three approaches

Three independent methods were used to measure net ecosystem production (NEP) in four wetlands near Thompson, Manitoba, Canada. The first method calculated NEP by subtracting heterotrophic respiration from net primary productivity, using both measurements and estimates derived from the literature. The second method used radiocarbon data from cores to derive long-term NEP averaged over the past several decades. The third method used direct measurement of NEP combined with a model to fill in for days with no data. The three methods, with their independently derived uncertainties, all show the same magnitude and pattern of NEP variation across four different wetland types. However, direct measurement yielded distinctly lower estimates of NEP in the most productive sites. Highest NEP (31–180 gC m −2 yr −1 ) was observed in the two wetlands with the highest proportion of sedge vegetation. A bog collapse scar and a nutrient-rich fen had NEP values not significantly different from zero. The maximum NEP at sites with intermediate nutrient status is due to slower overall decomposition and is likely associated with greater allocation of production below ground by sedges. The three methods for estimating NEP differ in the effort required, the sources of error, and in the timescale over which they apply. Used in combination, they allow estimation of parameters such as below- ground production and the contribution of heterotrophic decomposition to total soil respiration. Using the radiocarbon method, we also derived estimates of the rate of N accumulation in the four wetland types.

Manitoba↗
USGS2022

Combining eddy covariance and chamber methods to better constrain CO2 and CH4 fluxes across a heterogeneous restored tidal wetland

Tidal wetlands play an important role in global carbon cycling by storing carbon in sediment at millennial time scales, transporting dissolved carbon into coastal waters, and contributing significantly to global CH 4 budgets. However, these ecosystems' greenhouse gas monitoring and predictions are challenging due to spatial heterogeneity and tidal flooding. We utilized eddy covariance and chamber measurements to quantify fluxes of CO 2 and CH 4 at a restored tidal saltmarsh across spatial and temporal scales. Eddy covariance data revealed that the site was a strong net sink for CO 2 (−387 g C-CO 2 m −2 yr −1 , SD = 46) and a small net source of CH 4 (0.7 g C-CH 4 m −2 yr −1 , SD = 0.4). After partitioning net ecosystem exchange of CO 2 into gross primary production and ecosystem respiration, we found that high net uptake of CO 2 was due to low respiration emissions rather than high photosynthetic rates. We also found that respiration rates varied between land covers with increased respiration in mudflats compared to vegetated areas. Daytime soil chamber measurements revealed that the greatest CO 2 emission was from higher elevation mudflat soils (0.5 μmol m −2 s −1 , SE = 1.3) and CH 4 emission was greatest from lower elevation Spartina foliosa soils (1.6 nmol m −2 s −1 , SD = 8.2). Overall, these results highlight the importance of the relationships between wetland plant community and elevation, and inundation for CO 2 and CH 4 fluxes. Future research should include the use of high-resolution imagery, automated chambers, and a focus on quantifying carbon exported in tidal waters.

California↗
USGS2008

Evaluating the effects of historical land cover change on summertime weather and climate in New Jersey: Land cover and surface energy budget changes

The 19th-century agrarian landscape of New Jersey (NJ) and the surrounding region has been extensively transformed to the present-day land cover by urbanization, reforestation, and localized areas of deforestation. This study used a mesoscale atmospheric numerical model to investigate the sensitivity of the warm season climate of NJ to these land cover changes. Reconstructed 1880s-era and present-day land cover data sets were used as surface boundary conditions for a set of simulations performed with the Regional Atmospheric Modeling System (RAMS). Three-member ensembles with historical and present-day land cover were compared to examine the sensitivity of surface air and dew point temperatures, rainfall, and the individual components of the surface energy budget to these land cover changes. Mean temperatures for the present-day landscape were 0.3-0.6??C warmer than for the historical landscape over a considerable portion of NJ and the surrounding region, with daily maximum temperatures at least 1.0??C warmer over some of the highly urbanized locations. Reforested regions, however, were slightly cooler. Dew point temperatures decreased by 0.3-0.6??C, suggesting drier, less humid near-surface air for the present-day landscape. Surface warming was generally associated with repartitioning of net radiation from latent to sensible heat flux, and conversely for cooling. While urbanization was accompanied by strong surface albedo decreases and increases in net shortwave radiation, reforestation and potential changes in forest composition have generally increased albedos and also enhanced landscape heterogeneity. The increased deciduousness of forests may have further reduced net downward longwave radiation. Copyright 2008 by the American Geophysical Union.

New Jersey↗
USGS2019

Basal stress equations for granular debris masses on smooth or discretized slopes

Knowledge of basal stresses is essential for analyzing slope stability and modeling the dynamics and erosive potential of debris flows and avalanches. Here we derive and test new algebraic formulas for calculating the shear stress τ and normal stress σ at the base of variable‐thickness granular debris masses in states of static or dynamic equilibrium on slopes. The formulas include a lateral pressure coefficient κ , but use of a fixed value κ = 0.7 yields predictions of σ that on average err by less than 3% and of τ that on average err by less than 13% in matching basal stresses measured in six large‐scale experiments involving wet debris masses with varying geometries and compositions. Much larger prediction errors result from use of infinite‐slope or shallow‐debris approximations. Specialized versions of the new formulas apply if basal topography is discretized and represented by a “staircase” function in a digital elevation model. Use of these formulas to assess static limiting equilibrium conditions shows that the apparent basal Coulomb friction angle ϕ tread of debris that engages friction acting on the horizontal surfaces (or “treads”) of a staircase sloping at an angle θ is generally described by tan ϕ tread = tan ( ϕ − θ )+ κ tan θ , where ϕ is the true basal friction angle of the same debris in contact with a uniformly sloping bed. Differences between the values of ϕ and ϕ tread can greatly influence the results of numerical simulations that use unsmoothed digital elevation model topography to calculate the stability or dynamics of debris masses on slopes.

Journal of Geophysical Research-Earth Surface↗
USGS2021

Drivers of methane flux differ between lakes and reservoirs, complicating global upscaling efforts

Methane is an important greenhouse gas with growing atmospheric concentrations. Freshwater lakes and reservoirs contribute substantially to atmospheric methane concentrations, but the magnitude of this contribution is poorly constrained. Uncertainty stems partially from whether the sites currently sampled represent the global population as well as incomplete knowledge of which environmental variables predict methane flux. Thus, determining the main drivers of methane flux across diverse waterbody types will inform more accurate upscaling approaches. Here we use a new database of total, diffusive, and ebullitive areal methane emissions from 313 lakes and reservoirs (ranging in surface area from 6 m 2 to 5,400 km 2 ) to identify the best predictors of methane emission. We found that the best predictors of methane emission differed by waterbody type (lakes vs. reservoirs), and that ecosystem morphometric variables (e.g., surface area and maximum depth) were more important predictors in lakes whereas metrics of autochthonous production (e.g., chlorophyll a ) were more important in reservoirs. We also found that productivity strongly predicted methane ebullition, whereas ecosystem morphometry and waterbody type were more important predictors of diffusive methane flux. Finally, we identify several knowledge gaps that limit upscaling efforts. First, we need more methane emission measurements in small reservoirs, large lakes, and both natural and artificial ponds. Additionally, more accurate upscaling efforts require improved global information about waterbody surface area, waterbody type (lake vs. reservoir), ice phenology, and the distribution of productivity‐related predictor variables such as total phosphorus, DOC, and chlorophyll a .

Journal of Geophysical Research-Biogeosciences↗
USGS2021

Temporal seismic velocity variations: Recovery following from the 2019 Mw 7.1 Ridgecrest, California earthquake

We investigated seismic velocity changes ( dv / v ) associated with the 2019 Ridgecrest earthquake sequence with high‐frequency autocorrelations of ambient seismic noise data. Daily autocorrelation functions were computed for the entirety of 2019 and the first quarter of 2020 for broadband stations within the region, including the temporary broadband stations installed during the aftershock deployment. Travel time shifts in the daily autocorrelation functions, relative to the mean autocorrelation waveform, were computed to produce dv / v time series, which are sensitive to the evolving material properties of the shallow crust surrounding the Ridgecrest fault zone (RFZ). A short‐term velocity drop follows the M w 7.1 earthquake at stations in the vicinity of the rupture surface, while those greater than 50 km away showed no such drop. The maximum, absolute changes in seismic velocity are proportional to the logarithm of distance from the fault rupture and to the peak dynamic strain experienced during the earthquake. Near the areas of the highest coseismic slip within the RFZ, seismic velocities recovered over 3 months. However, in the vicinity of the nearby Garlock fault, where triggered slip manifested, and north of the RFZ, seismic velocities recovered within a month. We interpret the seismic velocity changes and their recovery to be largely due to changes in the physical properties of the shallow crust, such as fault zone damage recovery caused by the earthquake rupture process and in response to the large dynamic stresses of passing seismic waves from the mainshock.

California↗
USGS2021

Extreme precipitation across adjacent burned and unburned watersheds reveals impacts of low severity wildfire on debris-flow processes

In steep landscapes, wildfire-induced changes to soil and vegetation can lead to extreme and hazardous geomorphic responses, including debris flows. The wildfire-induced mechanisms that lead to heightened geomorphic responses, however, depend on many site-specific factors including regional climate, vegetation, soil texture, and soil burn severity. As climate and land use change drive changes in fire regime, there is an increasing need to understand how fire alters geomorphic responses, particularly in areas where fire has been historically infrequent. Here, we examine differences in the initiation, magnitude, and particle-size distribution of debris flows that initiated within the area burned by the 2019 Woodbury Fire in central Arizona, USA, and those that initiated in a nearby unburned area. Despite similar rainfall intensities, unburned watersheds were less likely to produce debris flows. Debris flows in unburned areas initiated from both runoff and shallow landslides, while debris flows only initiated from runoff-related processes in the burned area. The grain-size distribution making up the matrix of debris-flow deposits within the burned area generally had a lower ratio of sand to silt relative to debris flows that initiated in the unburned area, though there were no systematic differences in the coarse fraction of debris-flow sediment between burned and unburned areas. Results help expand our ability to predict postwildfire debris-flow activity in a wider range of settings, specifically the Sonoran Desert ecoregion, and provide general insight into the impact of wildfire on geomorphic processes in steep terrain.

Arizona↗
USGS2022

Synthetic evaluation of infrasonic multipole waveform inversion

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3-D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal-to-noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill-posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies can help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

Journal of Geophysical Research-Solid Earth↗
USGS2021

Influence of permafrost type and site history on losses of permafrost carbon after thaw

We quantified permafrost peat plateau and post-thaw carbon (C) stocks across a chronosequence in Interior Alaska to evaluate the amount of C lost with thaw. Macrofossil reconstructions revealed three stratigraphic layers of peat: (1) a base layer of fen/marsh peat, (2) peat from a forested peat plateau (with permafrost) and, (3) collapse-scar bog peat (at sites where permafrost thaw has occurred). Radiocarbon dating revealed that peat initiated within the last 2,500 years and that permafrost aggraded during the Little Ice Age (ca. 250 – 575 years ago) and degraded within the last several decades. The timing of permafrost thaw within each feature was not related to thaw bog size. Their rate of expansion may be more influenced by local factors, such as ground ice content and subsurface water inputs. We found C losses due to thaw over the past century were up to 46% of the C available, but the absolute amount of C lost (kg m -2 ) was over 50% lower than losses previously described in other Alaskan peatland chronosequences. We hypothesize that this difference stems from the process by which permafrost aggraded, with sites that formed permafrost epigenetically (significantly later than most peat accumulation) experiencing less absolute C loss with thaw than sites that formed syngenetically (simultaneously with peat accumulation). Epigenetic peat from our site had lower C:N ratios as compared to Alaskan sites that have syngenetic peat. This difference could help predict the magnitude of C loss with thaw across a range or permafrost types and histories.

Alaska↗
USGS2022

Higher temperature sensitivity of ecosystem respiration in low marsh compared to high elevation marsh ecosystems

Salt marsh habitats contain some of the highest quantities of soil organic carbon (C) per unit area, but increasing anthropogenic stressors threaten their ability to maintain themselves as large C reservoirs in some regions. We quantify rates of C gas exchange (methane [CH 4 ] and carbon dioxide [CO 2 ]) monthly across a 16-month period from a low nitrogen “reference” salt marsh on Cape Cod in New England using static chambers. While the summer period is the most dynamic period of marsh C gas exchange, we observed substantial fluxes in the early summer through late fall, highlighting the importance of including shoulder seasons in studies of marsh C exchange. We estimate annual ecosystem respiration between 108 and 252 g C m −2 yr −1 , which varied based on temperature and elevation. This flux is lower than in other nearby marshes, which we attribute to the frequently inundated, microtidal nature of the site, resulting in the majority of respired CO 2 being exported via lateral, not vertical, fluxes from this marsh. We observed significantly higher temperature sensitivity from the low elevation of the marsh compared to the high marsh. Recent acceleration in the rate of sea level rise is leading to a well-documented expansion of low marsh into high marsh vegetation zones in this marsh system and others in the region. While rates of C burial are higher in the low marsh compared to the high marsh, the higher temperature sensitivity of respiration in the low marsh may diminish the longevity of marsh C stocks with climate warming.

Massachusetts↗
USGS2012

Interactive coupling of a lake thermal model with a regional climate model

A one‐dimensional model of lake temperature, evaporation, and ice has been coupled in an interactive mode with the climate version of the National Center for Atmospheric Research/Pennsylvania State University regional (mesoscale) atmospheric model (MM4). The coupled MM4‐lake model makes possible high‐resolution simulations of climate in the proximity of large water bodies. Atmospheric variables required as input for the lake model are supplied by MM4 and simulated values of lake temperature, evaporation, and ice cover axe supplied to MM4 by the lake model. We have tested the coupled model system with a 60‐day, summertime simulation at Pyramid Lake, Nevada, and with a 10‐day, wintertime simulation of the North American Great Lakes and vicinity. Both simulations were conducted at a 60‐km resolution. Results from these tests indicate the coupled model system produces realistic simulations of lake temperature, evaporation, and ice cover and that the coupled system is applicable to simulations of regional climate change.

Journal of Geophysical Research-Atmospheres↗
USGS2024

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

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.

Alaska↗