SEARCH · Geology Reports

Results for “The Cryosphere”

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

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5Linked to original sources

USGS2016

Simulated high-latitude soil thermal dynamics during the past four decades

Soil temperature (Ts ) change is a key indicator of the dynamics of permafrost. On seasonal and inter-annual time scales, the variability of Ts determines the active layer depth, which regulates hydrological soil properties and biogeochemical processes. On the multi-decadal scale, increasing T 5 s not only drives permafrost thaw/retreat, but can also trigger and accelerate the decomposition of soil organic carbon. The magnitude of permafrost carbon feedbacks is thus closely linked to the rate of change of soil thermal regimes. In this study, we used nine process-based ecosystem models with permafrost processes, all forced by different observation-based climate forcing during the period 1960–2000, to characterize the warming rate of Ts 10 in permafrost regions. There is a large spread of Ts trends at 20 cm depth across the models, with trend values ranging from 0.010 ± 0.003 to 0.031 ± 0.005 ◦C yr−1 . Most models show smaller increase in Ts with increasing depth. Air temperature (Ta ) and longwave downward radiation (LWDR) are the main drivers of Ts trends, but their relative contributions differ 15 amongst the models. Different trends of LWDR used in the forcing of models can explain 61 % of their differences in Ts trends, while trends of Ta only explain 5 % of the differences in Ts trends. Uncertain climate forcing contributes a larger uncertainty in Ts trends (0.021 ± 0.008 ◦C yr−1 , mean ± SD) than the uncertainty of model structure (0.012 ± 0.001 ◦C yr−1 ), diagnosed from the range of response between different mod- 20 els, normalized to the same forcing. In addition, the loss rate of near-surface permafrost area, defined as total area where the maximum seasonal active layer thickness (ALT) is less than 3 m loss rate is found to be significantly correlated with the magnitude of the trends of Ts at 1 m depth across the models (R = −0.85, P = 0.003), but not with the initial total near-surface permafrost area (R = −0.30, P = 0.438). The sensitivity of the total boreal near-surface permafrost area to T 25 s at 1 m, is estimated to be of −2.80 ± 0.67 million km2 ◦C −1 . Finally, by using two long-term LWDR datasets and relationships between trends of LWDR and Ts across models, we infer an observationconstrained total boreal near-surface permafrost area decrease comprised between 39 ± 14 × 103 and 75 ± 14 × 103 km2 yr−1 from 1960 to 2000. This corresponds to 9– 18 % degradation of the current permafrost area.

Cryosphere Discussions↗

USGS2016

Diagnostic and model dependent uncertainty of simulated Tibetan permafrost area

We perform a land-surface model intercomparison to investigate how the simulation of permafrost area on the Tibetan Plateau (TP) varies among six modern stand-alone land-surface models (CLM4.5, CoLM, ISBA, JULES, LPJ-GUESS, UVic). We also examine the variability in simulated permafrost area and distribution introduced by five different methods of diagnosing permafrost (from modeled monthly ground temperature, mean annual ground and air temperatures, air and surface frost indexes). There is good agreement (99 to 135  ×  10 4 km 2 ) between the two diagnostic methods based on air temperature which are also consistent with the observation-based estimate of actual permafrost area (101  × 10 4 km 2 ). However the uncertainty (1 to 128  ×  10 4 km 2 ) using the three methods that require simulation of ground temperature is much greater. Moreover simulated permafrost distribution on the TP is generally only fair to poor for these three methods (diagnosis of permafrost from monthly, and mean annual ground temperature, and surface frost index), while permafrost distribution using air-temperature-based methods is generally good. Model evaluation at field sites highlights specific problems in process simulations likely related to soil texture specification, vegetation types and snow cover. Models are particularly poor at simulating permafrost distribution using the definition that soil temperature remains at or below 0 °C for 24 consecutive months, which requires reliable simulation of both mean annual ground temperatures and seasonal cycle, and hence is relatively demanding. Although models can produce better permafrost maps using mean annual ground temperature and surface frost index, analysis of simulated soil temperature profiles reveals substantial biases. The current generation of land-surface models need to reduce biases in simulated soil temperature profiles before reliable contemporary permafrost maps and predictions of changes in future permafrost distribution can be made for the Tibetan Plateau.

Tibeta Plateau↗

USGS2017

Observationally constrained surface mass balance of Larsen C ice shelf, Antarctica

The surface mass balance (SMB) of the Larsen C ice shelf (LCIS), Antarctica, is poorly constrained due to a dearth of in situ observations. Combining several geophysical techniques, we reconstruct spatial and temporal patterns of SMB over the LCIS. Continuous time series of snow height (2.5–6 years) at five locations allow for multi-year estimates of seasonal and annual SMB over the LCIS. There is high interannual variability in SMB as well as spatial variability: in the north, SMB is 0.40 ± 0.06 to 0.41 ± 0.04 m w.e. year −1 , while farther south, SMB is up to 0.50 ± 0.05 m w.e. year −1 . This difference between north and south is corroborated by winter snow accumulation derived from an airborne radar survey from 2009, which showed an average snow thickness of 0.34 m w.e. north of 66° S, and 0.40 m w.e. south of 68° S. Analysis of ground-penetrating radar from several field campaigns allows for a longer-term perspective of spatial variations in SMB: a particularly strong and coherent reflection horizon below 25–44 m of water-equivalent ice and firn is observed in radargrams collected across the shelf. We propose that this horizon was formed synchronously across the ice shelf. Combining snow height observations, ground and airborne radar, and SMB output from a regional climate model yields a gridded estimate of SMB over the LCIS. It confirms that SMB increases from north to south, overprinted by a gradient of increasing SMB to the west, modulated in the west by föhn-induced sublimation. Previous observations show a strong decrease in firn air content toward the west, which we attribute to spatial patterns of melt, refreezing, and densification rather than SMB.

Larsen C Ice Shelf↗

USGS2018

Local topography increasingly influences the mass balance of a retreating cirque glacier

Local topographically driven processes – such as wind drifting, avalanching, and shading – are known to alter the relationship between the mass balance of small cirque glaciers and regional climate. Yet partitioning such local effects from regional climate influence has proven difficult, creating uncertainty in the climate representativeness of some glaciers. We address this problem for Sperry Glacier in Glacier National Park, USA, using field-measured surface mass balance, geodetic constraints on mass balance, and regional climate data recorded at a network of meteorological and snow stations. Geodetically derived mass changes during 1950–1960, 1960–2005, and 2005–2014 document average mass change rates during each period at −0.22 ± 0.12, −0.18 ± 0.05, and −0.10 ± 0.03 m w.e. yr −1 , respectively. A correlation of field-measured mass balance and regional climate variables closely (i.e., within 0.08 m w.e. yr −1 ) predicts the geodetically measured mass loss from 2005 to 2014. However, this correlation overestimates glacier mass balance for 1950–1960 by +1.20 ± 0.95 m w.e. yr −1 . Our analysis suggests that local effects, not represented in regional climate variables, have become a more dominant driver of the net mass balance as the glacier lost 0.50 km 2 and retreated further into its cirque.

Northern Rocky Mountains↗

USGS2018

Interannual snow accumulation variability on glaciers derived from repeat spatially extensive ground-penetrating radar surveys

There is significant uncertainty regarding the spatiotemporal distribution of seasonal snow on glaciers, despite being a fundamental component of glacier mass balance. To address this knowledge gap, we collected repeat, spatially extensive high-frequency ground-penetrating radar (GPR) observations on two glaciers in Alaska for five consecutive years. GPR measurements showed steep snow water equivalent (SWE) elevation gradients at both sites; continental Gulkana Glacier’s SWE gradient averaged 115 mm 100 m–1 and maritime Wolverine Glacier’s gradient averaged 440 mm 100 m–1 (over >1000 m). We extrapolated GPR point observations across the glacier surface using terrain parameters derived from digital elevation models as predictor variables in two statistical models (stepwise multivariable linear regression and regression trees). Elevation and proxies for wind redistribution had the greatest explanatory power, and exhibited relatively time-constant coefficients over the study period. Both statistical models yielded comparable estimates of glacier-wide average SWE (1 % average difference at Gulkana, 4 % average difference at Wolverine), although the spatial distributions produced by the models diverged in unsampled regions of the glacier, particularly at Wolverine. In total, six different methods for estimating the glacier-wide average agreed within ± 11 %. We assessed interannual variability in the spatial pattern of snow accumulation predicted by the statistical models using two quantitative metrics. Both glaciers exhibited a high degree of temporal stability, with ~85 % of the glacier area experiencing less than 25 % normalized absolute variability over this five-year interval. We found SWE at a sparse network (3 stakes per glacier) of long-term glaciological stake sites to be highly correlated with the GPR-derived glacier-wide average. We estimate that interannual variability in the spatial pattern of SWE is only a small component (4–10 % of glacier-wide average) of the total mass balance uncertainty and thus, our findings support the concept that sparse stake networks effectively measure interannual variability in winter balance on glaciers, rather than some spatially varying pattern of snow accumulation.

Alaska↗

USGS2020

The firn meltwater Retention Model Intercomparison Project (RetMIP): Evaluation of nine firn models at four weather station sites on the Greenland ice sheet

Perennial snow, or firn, covers 80 % of the Greenland ice sheet and has the capacity to retain surface meltwater, influencing the ice sheet mass balance and contribution to sea-level rise. Multilayer firn models are traditionally used to simulate firn processes and estimate meltwater retention. We present, intercompare and evaluate outputs from nine firn models at four sites that represent the ice sheet's dry snow, percolation, ice slab and firn aquifer areas. The models are forced by mass and energy fluxes derived from automatic weather stations and compared to firn density, temperature and meltwater percolation depth observations. Models agree relatively well at the dry-snow site while elsewhere their meltwater infiltration schemes lead to marked differences in simulated firn characteristics. Models accounting for deep meltwater percolation overestimate percolation depth and firn temperature at the percolation and ice slab sites but accurately simulate recharge of the firn aquifer. Models using Darcy's law and bucket schemes compare favorably to observed firn temperature and meltwater percolation depth at the percolation site, but only the Darcy models accurately simulate firn temperature and percolation at the ice slab site. Despite good performance at certain locations, no single model currently simulates meltwater infiltration adequately at all sites. The model spread in estimated meltwater retention and runoff increases with increasing meltwater input. The highest runoff was calculated at the KAN_U site in 2012, when average total runoff across models ( ±2 σ ) was 353±610 mm w.e. (water equivalent), about 27±48 % of the surface meltwater input. We identify potential causes for the model spread and the mismatch with observations and provide recommendations for future model development and firn investigation.

The Cryosphere↗

USGS2021

The Holocene dynamics of Ryder Glacier and ice tongue in north Greenland

The northern sector of the Greenland Ice Sheet is considered to be particularly susceptible to ice mass loss arising from increased glacier discharge in the coming decades. However, the past extent and dynamics of outlet glaciers in this region, and hence their vulnerability to climate change, are poorly documented. In the summer of 2019, the Swedish icebreaker Oden entered the previously unchartered waters of Sherard Osborn Fjord, where Ryder Glacier drains approximately 2 % of Greenland's ice sheet into the Lincoln Sea. Here we reconstruct the Holocene dynamics of Ryder Glacier and its ice tongue by combining radiocarbon dating with sedimentary facies analyses along a 45 km transect of marine sediment cores collected between the modern ice tongue margin and the mouth of the fjord. The results illustrate that Ryder Glacier retreated from a grounded position at the fjord mouth during the Early Holocene ( > 10.7±0.4 ka cal BP) and receded more than 120 km to the end of Sherard Osborn Fjord by the Middle Holocene ( 6.3±0.3 ka cal BP), likely becoming completely land-based. A re-advance of Ryder Glacier occurred in the Late Holocene, becoming marine-based around 3.9±0.4 ka cal BP. An ice tongue, similar in extent to its current position was established in the Late Holocene (between 3.6±0.4 and 2.9±0.4 ka cal BP) and extended to its maximum historical position near the fjord mouth around 0.9±0.3 ka cal BP. Laminated, clast-poor sediments were deposited during the entire retreat and regrowth phases, suggesting the persistence of an ice tongue that only collapsed when the glacier retreated behind a prominent topographic high at the landward end of the fjord. Sherard Osborn Fjord narrows inland, is constrained by steep-sided cliffs, contains a number of bathymetric pinning points that also shield the modern ice tongue and grounding zone from warm Atlantic waters, and has a shallowing inland sub-ice topography. These features are conducive to glacier stability and can explain the persistence of Ryder's ice tongue while the glacier remained marine-based. However, the physiography of the fjord did not halt the dramatic retreat of Ryder Glacier under the relatively mild changes in climate forcing during the Holocene. Presently, Ryder Glacier is grounded more than 40 km seaward of its inferred position during the Middle Holocene, highlighting the potential for substantial retreat in response to ongoing climate change.

Ryder Glacier↗

USGS2021

Invited perspective: What lies beneath a changing Arctic?

As permafrost thaws in the Arctic, new subsurface pathways open for the transport of groundwater, energy, and solutes. We identify different ways that these subsurface changes are driving observed surface consequences, including the potential for increased contaminant transport, modification to water resources, and enhanced rates of infrastructure (e.g. buildings and roads) damage. Further, as permafrost thaws it allows groundwater to transport carbon, nutrients, and other dissolved constituents from terrestrial to aquatic environments via progressively deeper subsurface flow paths. Cryohydrogeology, the study of groundwater in cold regions, should be included in northern research initiatives to account for this hidden catalyst of environmental and societal change.

The Cryosphere↗

USGS2022

The effect of changing sea ice on wave climate trends along Alaska's central Beaufort Sea coast

Diminishing sea ice is impacting the wave field across the Arctic region. Recent observation- and model-based studies highlight the spatiotemporal influence of sea ice on offshore wave climatologies, but effects within the nearshore region are still poorly described. This study characterizes the wave climate in the central Beaufort Sea coast from 1979 to 2019 by utilizing a wave hindcast model that uses ERA5 winds, waves, and ice concentrations as input. The spectral wave model SWAN (Simulating Waves Nearshore) is calibrated and validated based on more than 10 000 in situ time point measurements collected over a 13-year time period across the region, with friction variations and empirical coefficients for newly implemented empirical ice formulations for the open-water and shoulder seasons. Model results and trends are analyzed over the 41-year time period using the non-parametric Mann–Kendall test, including an estimate of Sen's slope. The model results show that the reduction in sea ice concentration correlates strongly with increases in average and extreme wave conditions. In particular, the open-water season extended by ∼96 d over the 41-year time period (∼2.4 d yr−1), resulting in a 5-fold increase in the yearly cumulative wave power. Moreover, the open-water season extends later into the year, resulting in relatively more open-water conditions during fall storms with high wind speeds. The later freeze-up results in an increase in the annual offshore median wave heights of 1 % yr−1 and an increase in the average number of rough wave days (defined as days when maximum wave heights exceed 2.5 m) from 1.5 in 1979 to 13.1 d in 2019. Trends in the nearshore areas deviate from the patterns offshore. Model results indicate a saturation limit for high wave heights in the shallow areas of Foggy Island Bay. Similar patterns are found for yearly cumulative wave power.

Alaska↗

USGS2022

Assessing the seasonal evolution of snow depth spatial variability and scaling in complex mountain terrain

Dynamic natural processes govern snow distribution in mountainous environments throughout the world. Interactions between these different processes create spatially variable patterns of snow depth across a landscape. Variations in accumulation and redistribution occur at a variety of spatial scales, which are well established for moderate mountain terrain. However, spatial patterns of snow depth variability in steep, complex mountain terrain have not been fully explored due to insufficient spatial resolutions of snow depth measurement. Recent advances in uncrewed aerial systems (UASs) and structure from motion (SfM) photogrammetry provide an opportunity to map spatially continuous snow depths at high resolutions in these environments. Using UASs and SfM photogrammetry, we produced 11 snow depth maps at a steep couloir site in the Bridger Range of Montana, USA, during the 2019–2020 winter. We quantified the spatial scales of snow depth variability in this complex mountain terrain at a variety of resolutions over 2 orders of magnitude (0.02 to 20 m) and time steps (4 to 58 d) using variogram analysis in a high-performance computing environment. We found that spatial resolutions greater than 0.5 m do not capture the complete patterns of snow depth spatial variability within complex mountain terrain and that snow depths are autocorrelated within horizontal distances of 15 m at our study site. The results of this research have the potential to reduce uncertainty currently associated with snowpack and snow water resource analysis by documenting and quantifying snow depth variability and snowpack evolution on relatively inaccessible slopes in complex terrain at high spatial and temporal resolutions.

Montana↗

USGS2024

Observing glacier elevation changes from spaceborne optical and radar sensors – an inter-comparison experiment using ASTER and TanDEM-X data

Observations of glacier mass changes are key to understanding the response of glaciers to climate change and related impacts, such as regional runoff, ecosystem changes, and global sea level rise. Spaceborne optical and radar sensors make it possible to quantify glacier elevation changes, and thus multi-annual mass changes, on a regional and global scale. However, estimates from a growing number of studies show a wide range of results with differences often beyond uncertainty bounds. Here, we present the outcome of a community-based inter-comparison experiment using spaceborne optical stereo (ASTER) and synthetic aperture radar interferometry (TanDEM-X) data to estimate elevation changes for defined glaciers and target periods that pose different assessment challenges. Using provided or self-processed digital elevation models (DEMs) for five test sites, 12 research groups provided a total of 97 spaceborne elevation-change datasets using various processing approaches. Validation with airborne data showed that using an ensemble estimate is promising to reduce random errors from different instruments and processing methods but still requires a more comprehensive investigation and correction of systematic errors. We found that scene selection, DEM processing, and co-registration have the biggest impact on the results. Other processing steps, such as treating spatial data voids, differences in survey periods, or radar penetration, can still be important for individual cases. Future research should focus on testing different implementations of individual processing steps (e.g. co-registration) and addressing issues related to temporal corrections, radar penetration, glacier area changes, and density conversion. Finally, there is a clear need for our community to develop best practices, use open, reproducible software, and assess overall uncertainty to enhance inter-comparison and empower physical process insights across glacier elevation-change studies.

The Cryosphere↗

USGS2024

Brief communication: Recent estimates of glacier mass loss for western North America from laser altimetry

Glaciers in western North American outside of Alaska are often overlooked in global studies because their potential to contribute to changes in sea level is small. Nonetheless, these glaciers represent important sources of freshwater, especially during times of drought. Differencing recent ICESat-2 data from a digital elevation model derived from a combination of synthetic aperture radar data (TerraSAR-X/TanDEM-X), we find that over the period 2013–2020, glaciers in western North America lost mass at a rate of -12.3 ± 3.5 Gt yr −1 . This rate is comparable to the rate of mass loss (-11.7 ± 1.0 Gt yr −1 ) for the period 2018–2022 calculated through trend analysis using ICESat-2 and Global Ecosystems Dynamics Investigation (GEDI) data.

western North America↗

USGS2025

Automated snow cover detection on mountain glaciers usingspaceborne imagery and machine learning

Tracking the extent of seasonal snow on glaciers over time is critical for assessing glacier vulnerability and the response of glacierized watersheds to climate change. Existing snow cover products do not reliably distinguish seasonal snow from glacier ice and firn, preventing their use for glacier snow cover detection. Despite previous efforts to classify glacier surface facies using machine learning on local scales, currently there is no published comparison of machine learning models for classifying glacier snow cover across different satellite image products. We present an automated snow detection workflow for mountain glaciers using supervised machine-learning-based image classifiers and Landsat 8 and 9, Sentinel-2, and PlanetScope satellite imagery. We develop the image classifiers by testing numerous machine learning algorithms with training and validation data from the U.S. Geological Survey Benchmark Glacier Project glaciers. The workflow produces daily to twice monthly time series of several glacier mass balance and snowmelt indicators (snow-covered area, accumulation area ratio, and seasonal snow line) from 2013 to present. Workflow performance is assessed by comparing automatically classified images and snow lines to manual interpretations at each glacier site. The image classifiers exhibit overall accuracies of 92%–98%, K scores of 84%–96%, and F scores of 93%–98% for all image products. The median difference between automatically and manually delineated median snow line altitudes is 31m (IQR of 73to0m)across all image products. The Sentinel-2 classifier (support vector machine) produces the most accurate glacier mass balance and snowmelt indicators and distinguishes snow from ice and f irn the most reliably. Although they are less accurate, the Landsat- and PlanetScope-derived estimates greatly enhance the temporal coverage of observations. The transient accumulation area ratio produces the least noisy time series, making it the most reliable indicator for characterizing seasonal snow trends. The temporally detailed accumulation area ratio time series reveal that the timing of minimum snow cover conditions varies by up to a month between Arctic (63°N) and midlatitude (48°N) sites, underscoring the potential for bias when estimating glacier minimum snow cover conditions from a single late-summer image. Widespread application of our automated snow detection workflow has the potential to improve regional assessments of glacier mass balance, land ice representations within Earth system models, water resources, and the impacts of climate change on snow cover across broad spatial scales.

The Cryosphere↗

USGS2025

Brief communication: Not as dirty as they look, flawed airborne and satellite snow spectra

Key to the success of spaceborne missions is understanding snowmelt in our warming climate, as this has implications for nearly 2 billion people. An obstacle is that surface reflectance products over snow show an erroneous hook with decreases in the visible wavelengths, causing per-band and broadband reflectance errors of up to 33 % and 11 %, respectively. This hook is sometimes mistaken for soot or dust but can result from three artifacts: (1) background reflectance that is too dark, (2) an assumption of level terrain, or (3) differences in optical constants of ice. Sensor calibration and directional effects may also contribute. Solutions are being implemented.

The Cryosphere↗

USGS2026

Evaluating snow depth measurements from ground-penetrating radar and airborne lidar in boreal forest and tundra environments during the NASA SnowEx 2023 campaign

Snow is a vital component of high-latitude terrestrial systems, but environmental factors (e.g., permafrost) and complex vegetation challenge the accurate measurement of key snowpack properties. We evaluated local-scale ground-penetrating radar (GPR) and large-scale airborne lidar retrievals of snow depth collected during the NASA SnowEx 2023 campaign in tundra and boreal forest environments in Alaska along 44 short (3–12 m) transects. Compared to in situ observations, we identified modest biases for GPR snow depths (bias <0.03 m in tundra, +0.06 m in boreal forests) and larger biases for lidar snow depths in the boreal forests (–0.16 m). At the Upper Kuparuk-Toolik tundra site, lidar snow depths exhibited a small bias (–0.02 m), whereas the bias was much larger at the Arctic Coastal Plain tundra site (+0.19 m). For most sites, biases were primarily related to sub-snow vegetation, tussocks, and seasonally dynamic ground. However, we identified vertical alignment issues with the Arctic Coastal Plain lidar snow depth dataset that likely contributed to the higher bias. The complex ground surface and sub-snow vegetation in these environments present a challenge to established snow depth measurement methods, which needs to be considered when evaluating novel remote sensing approaches.

Alaska↗

USGS2017

In situ nuclear magnetic resonance response of permafrost and active layer soil in boreal and tundra ecosystems

Characterization of permafrost, particularly warm and near-surface permafrost which can contain significant liquid water, is critical to understanding complex interrelationships with climate change, ecosystems, and disturbances such as wildfires. Understanding the vulnerability and resilience of permafrost requires an interdisciplinary approach, relying on (for example) geophysical investigations, ecological characterization, direct observations, remote sensing, and more. As part of a multi-year investigation into the impacts of wildfires to permafrost, we have collected in situ measurements of the nuclear magnetic resonance (NMR) response of active layer and permafrost in a variety of soil conditions, types, and saturations. In this paper, we summarize the NMR data and present quantitative relationships between active layer and permafrost liquid water content and pore sizes. Through statistical analyses and synthetic freezing simulations, we also demonstrate that borehole NMR can image the nucleation of ice within soil pore spaces.

Alaska↗

USGS2016

Evaluation of air-soil temperature relationships simulated by land surface models during winter across the permafrost region

A realistic simulation of snow cover and its thermal properties are important for accurate modelling of permafrost. We analyze simulated relationships between air and near-surface (20 cm) soil temperatures in the Northern Hemisphere permafrost region during winter, with a particular focus on snow insulation effects in nine land surface models and compare them with observations from 268 Russian stations. There are large across-model differences as expressed by simulated differences between near-surface soil and air temperatures, (Δ T ), of 3 to 14 K, in the gradients between soil and air temperatures (0.13 to 0.96°C/°C), and in the relationship between Δ T and snow depth. The observed relationship between Δ T and snow depth can be used as a metric to evaluate the effects of each model's representation of snow insulation, and hence guide improvements to the model’s conceptual structure and process parameterizations. Models with better performance apply multi-layer snow schemes and consider complex snow processes. Some models show poor performance in representing snow insulation due to underestimation of snow depth and/or overestimation of snow conductivity. Generally, models identified as most acceptable with respect to snow insulation simulate reasonable areas of near-surface permafrost (12–16 million km 2 ). However, there is not a simple relationship between the quality of the snow insulation in the acceptable models and the simulated area of Northern Hemisphere near-surface permafrost, likely because several other factors such as differences in the treatment of soil organic matter, soil hydrology, surface energy calculations, and vegetation also provide important controls on simulated permafrost distribution.

The Cryosphere↗

USGS2016

Presence of rapidly degrading permafrost plateaus in south-central Alaska

Permafrost presence is determined by a complex interaction of climatic, topographic, and ecological conditions operating over long time scales. In particular, vegetation and organic layer characteristics may act to protect permafrost in regions with a mean annual air temperature (MAAT) above 0 °C. In this study, we document the presence of residual permafrost plateaus in the western Kenai Peninsula lowlands of south-central Alaska, a region with a MAAT of 1.5 ± 1 °C (1981–2010). Continuous ground temperature measurements between 16 September 2012 and 15 September 2015, using calibrated thermistor strings, documented the presence of warm permafrost (−0.04 to −0.08 °C). Field measurements (probing) on several plateau features during the fall of 2015 showed that the depth to the permafrost table averaged 1.48 m but at some locations was as shallow as 0.53 m. Late winter surveys (augering, coring, and GPR) in 2016 showed that the average seasonally frozen ground thickness was 0.45 m, overlying a talik above the permafrost table. Measured permafrost thickness ranged from 0.33 to > 6.90 m. Manual interpretation of historic aerial photography acquired in 1950 indicates that residual permafrost plateaus covered 920 ha as mapped across portions of four wetland complexes encompassing 4810 ha. However, between 1950 and ca. 2010, permafrost plateau extent decreased by 60.0 %, with lateral feature degradation accounting for 85.0 % of the reduction in area. Permafrost loss on the Kenai Peninsula is likely associated with a warming climate, wildfires that remove the protective forest and organic layer cover, groundwater flow at depth, and lateral heat transfer from wetland surface waters in the summer. Better understanding the resilience and vulnerability of ecosystem-protected permafrost is critical for mapping and predicting future permafrost extent and degradation across all permafrost regions that are currently warming. Further work should focus on reconstructing permafrost history in south-central Alaska as well as additional contemporary observations of these ecosystem-protected permafrost sites south of the regions with relatively stable permafrost.

Alaska↗