Geology topics

Christopher J. McNeil

Publications and source records attributed to Christopher J. McNeil.

14 recordsLinked to original sources

USGS2024

GNSS reflectometry from low-cost sensors for continuous in situ contemporaneous glacier mass balance and flux divergence

Recent advances in remote sensing have produced global glacier surface elevation change data. Parsing these elevation change signals into contributions from the climate (i.e. climatic mass balance) and glacier dynamics (i.e. flux divergence) is critical to enhance our process-based understanding of glacier change. In this study, we evaluate three approaches for direct, continuous measurements of the climatic mass balance, flux divergence, and elevation change at a site on Gulkana Glacier in Alaska using low-cost GNSS sensors, GNSS interferometric reflectometry (GNSS-IR), banded ablation stakes with time-lapse cameras, and combinations thereof. Cumulative climatic mass balance over the season was 4.85 m and the three approaches were within 0.08 m through early July before the snowpack melted, and within 0.28 m through mid-August. The flux divergence increased from 0.52 ± 0.03 cm d -1 before June 3 to roughly 0.73 cm d -1 after June 27. We demonstrate a single GNSS system fixed atop an ablation stake can measure contemporaneous climatic mass balance, flux divergence, and elevation change based on the antenna’s position and GNSS-IR techniques. The ability of these systems to measure glacier mass balance and flux divergence offers unique opportunities for year-round observations on mountain glaciers in the future.

Alaska↗

USGS2024

Direct measurements of firn-density evolution from 2016 to 2022 at Wolverine Glacier, Alaska

Knowledge of snow and firn-density change is needed to use elevation-change measurements to estimate glacier mass change. Additionally, firn-density evolution on glaciers is closely connected to meltwater percolation, refreezing and runoff, which are key processes for glacier mass balance and hydrology. Since 2016, the U.S. Geological Survey Benchmark Glacier Project has recovered firn cores from a site on Wolverine Glacier in Alaska's Kenai Mountains. We use annual horizons in repeat cores to track firn densification and meltwater retention over seasonal and interannual timescales, and we use density measurements to quantify how the firn air content (FAC) changes through time. The results suggest the firn is densifying due primarily to compaction rather than refreezing. Liquid-water retention in the firn is transient, likely due to gravity-fed drainage and irreducible-water-content decreases that accompany decreasing porosity. We show that the uncertainty (±60 kg m −3 ) in the commonly used volume-to-mass conversion factor of 850 kg m −3 is an underestimation when glacier-wide FAC variability exceeds 12% of the glacier-averaged height change. Our results demonstrate how direct measurements of firn properties on mountain glaciers can be used to better quantify the uncertainty in geodetic volume-to-mass conversions.

Alaska↗

USGS2024

Accelerating glacier volume loss on Juneau Icefield driven by hypsometry and melt-accelerating feedbacks

Globally, glaciers and icefields contribute significantly to sea level rise. Here we show that ice loss from Juneau Icefield, a plateau icefield in Alaska, accelerated after 2005 AD. Rates of area shrinkage were 5 times faster from 2015–2019 than from 1979–1990. Glacier volume loss remained fairly consistent (0.65–1.01 km 3 a −1 ) from 1770–1979 AD, rising to 3.08–3.72 km 3 a −1 from 1979–2010, and then doubling after 2010 AD, reaching 5.91 ± 0.80 km 3 a −1 (2010–2020). Thinning has become pervasive across the icefield plateau since 2005, accompanied by glacier recession and fragmentation. Rising equilibrium line altitudes and increasing ablation across the plateau has driven a series of hypsometrically controlled melt-accelerating feedbacks and resulted in the observed acceleration in mass loss. As glacier thinning on the plateau continues, a mass balance-elevation feedback is likely to inhibit future glacier regrowth, potentially pushing glaciers beyond a dynamic tipping point.

Alaska↗

USGS2023

How to handle glacier area change in geodetic mass balance

Innovations in geodesy enable widespread analysis of glacier surface elevation change and geodetic mass balance. However, coincident glacier area data are less widely available, causing inconsistent handling of glacier area change. Here we quantify the bias introduced into meters water equivalent (m w.e.) specific geodetic mass balance results when using a fixed, maximum glacier area, and illustrate the bias for five North American glaciers. Sites span latitudes from the northern U.S. Rocky Mountains (48°N) to the Central Alaska Range (63°N) between 1948 and 2021. Results show that fixed (maximum) area treatment subdues the m w.e. mass change signal, underestimating mass balance by up to 19% in our test cases. This bias scales with relative glacier area change and the mass balance magnitude. Thus, the bias for specific geodetic mass balances will be most pronounced across rapidly deglaciating regions. Our analysis underscores the need for temporally resolved glacier area in geodetic mass balance studies.

Journal of Glaciology↗

USGS2023

Inventory of glaciers and perennial snowfields of the conterminous USA

This report summarizes an updated inventory of glaciers and perennial snowfields of the conterminous United States. The inventory is based on interpretation of mostly aerial imagery provided by the National Agricultural I magery Program, US Department of Agriculture, with some satellite imagery in places where aerial imagery was not suitable. The inventory includes all perennial snow and ice features ≥ 0.01 km 2 . Due to aerial survey schedules and seasonal snow cover, imageries acquired over a number of years were required. The earliest date is 2013 and the latest is 2020, but more than 73 % of the outlines were acquired from 2015 imagery. The inventory is compiled as shapefiles within a geographic information system that includes feature classification, area, and location. The inventory identified 1331 (366.52 ± 14.34 km 2 ) glaciers, 1176 (31.01 ± 9.30 km 2 ) perennial snowfields, and 35 (3.57 km 2 ± no uncertainty) buried-ice features. The data including both the shapefiles and tabulated results are publicly available at https://doi.org/10.15760/geology-data.03 (Fountain and Glenn, 2022).

Earth System Science Data↗

USGS2023

Historical Structure from Motion (HSfM): Automated processing of historical aerial photographs for long-term topographic change analysis

Precisely measuring the Earth’s changing surface on decadal to centennial time scales is critical for many science and engineering applications, yet long-term records of quantitative landscape change are often temporally and geographically sparse. Archives of scanned historical aerial photographs provide an opportunity to augment these records with accurate elevation measurements that capture the historical state of the Earth surface. Structure from Motion (SfM) photogrammetry workflows produce high-quality digital elevation models (DEMs) and orthoimage mosaics from these historical images, but time-intensive tasks like manual image preprocessing (e.g., fiducial marker identification) and ground control point (GCP) selection impede processing at scale. We developed an automated method to process historical images and generate self-consistent time series of high-resolution (0.5–2 m) DEMs and orthomosaics, without manual GCP selection. The method relies on SfM to correct camera interior and exterior orientation and a robust multi-stage co-registration approach using modern reference terrain datasets for geolocation refinement. We demonstrate the method using scanned images from the North American Glacier Aerial Photography (NAGAP) archive collected between 1967 and 1997. We present results for two sites with variable photo acquisition geometry and overlap — Mount Baker and South Cascade Glacier in Washington State, USA. The automated method corrects initial camera position errors of several kilometers and produces accurately georeferenced, high-resolution DEMs and orthoimages, regardless of camera configuration, acquisition geometry, terrain characteristics, and reference DEM properties. The average RMS reprojection error following bundle adjustment optimization was 0.67 px (0.15 m) for the 261 images contributing to 10 final DEM mosaics between 1970 and 1992 at Mount Baker, and 0.65 px (0.13 m) for the 243 images contributing to 18 individual DEMs between 1967 and 1997 at South Cascade Glacier. The relative accuracy of elevation values in the historical time series stacks was 0.68 m at Mount Baker and 0.37 m at South Cascade Glacier. Our products have reduced systematic error and improved accuracy compared to DEM products generated using SfM with manual GCP selection. Final elevation change measurement precision was ∼ "> ∼ 0.7–1.0 m over a 30-year period, enabling the study of processes with rates as low as ∼ "> ∼ 1-3 cm/yr. Our results demonstrate the potential of this scalable method to rapidly process archives of historical imagery and deliver new quantitative insights on long-term geodetic change and Earth surface processes.

Washington↗

USGS2022

Beyond glacier-wide mass balances: Parsing seasonal elevation change into spatially resolved patterns of accumulation and ablation at Wolverine Glacier, Alaska

We present spatially distributed seasonal and annual surface mass balances of Wolverine Glacier, Alaska, from 2016 to 2020. Our approach accounts for the effects of ice emergence and firn compaction on surface elevation changes to resolve the spatial patterns in mass balance at 10 m scale. We present and compare three methods for estimating emergence velocities. Firn compaction was constrained by optimizing a firn model to fit three firn cores. Distributed mass balances showed good agreement with mass-balance stakes (RMSE = 0.67 m w.e., r = 0.99, n = 41) and ground-penetrating radar surveys (RMSE = 0.36 m w.e., r = 0.85, n = 9024). Fundamental differences in the distributions of seasonal balances highlight the importance of disparate physical processes, with anomalously high ablation rates observed in icefalls. Winter balances were found to be positively skewed when controlling for elevation, while summer and annual balances were negatively skewed. We show that only a small percent of the glacier surface represents ideal locations for mass-balance stake placement. Importantly, no suitable areas are found near the terminus or in elevation bands dominated by icefalls. These findings offer explanations for the often-needed geodetic calibrations of glaciological time series.

Alaska↗

USGS2022

Topographic controls on ice flow and recession for Juneau Icefield (Alaska/British Columbia)

Globally, mountain glaciers and ice caps are losing dramatic volumes of ice. The resultant sea-level rise is dominated by contributions from Alaska. Plateau icefields may be especially sensitive to climate change due to the non-linear controls their topography imparts on their response to climate change. However, Alaskan plateau icefields have been subject to little structural glaciological or regional geomorphological assessment, which makes the controls on their present and former mass balance difficult to ascertain. We inventoried 1050 glaciers and 368 lakes in the Juneau Icefield region for the year 2019. We found that 63 glaciers had disappeared since the 2005 inventory, with a reduction in glacier area of 422 km 2 (10.0%). We also present the first structural glaciological and geomorphological map for an entire icefield in Alaska. Glaciological mapping of >20 800 features included crevasses, debris cover, foliation, ogives, medial moraines and, importantly, areas of glacier fragmentation, where glaciers either separated from tributaries via lateral recession ( n = 59), or disconnected within areas of former icefalls ( n = 281). Geomorphological mapping of >10 200 landforms included glacial moraines, glacial lakes, trimlines, flutes and cirques. These landforms were generated by a temperate icefield during the Little Ice Age (LIA) neoglaciation. These data demonstrate that the present-day outlet glaciers, which have a similar thermal and ice-flow regime, have undergone largely continuous recession since the LIA. Importantly, disconnections occurring within glaciers can separate accumulation and ablation zones, increasing rates of glacier mass loss. We show that glacier disconnections are widespread across the icefield and should be critically taken into consideration when icefield vulnerability to climate change is considered.

Alaska↗

USGS2022

Uncertainty of ICESat-2 ATL06- and ATL08-derived snow depths for glacierized and vegetated mountain regions

Seasonal snow melt dominates the hydrologic budget across a large portion of the globe. Snow accumulation and melt vary over a broad range of spatial scales, preventing accurate extrapolation of sparse in situ observations to watershed scales. The lidar onboard the Ice, Cloud, and land Elevation, Satellite (ICESat-2) was designed for precise mapping of ice sheets and sea ice, and here we assess the feasibility of snow depth-mapping using ICESat-2 data in more complex and rugged mountain landscapes. We explore the utility of ATL08 Land and Vegetation Height and ATL06 Land Ice Height differencing from reference elevation datasets in two end member study sites. We analyze ∼3 years of data for Reynolds Creek Experimental Watershed in Idaho's Owyhee Mountains and Wolverine Glacier in southcentral Alaska's Kenai Mountains. Our analysis reveals decimeter-scale uncertainties in derived snow depth and glacier mass balance at the watershed scale. Both accuracy and precision decrease as slope increases: the magnitudes of the median and median of the absolute deviation of elevation errors (MAD) vary from ∼0.2 m for slopes <5° to >1 m for slopes >20°. For glacierized regions, failure to account for intra- and inter-annual evolution of glacier surface elevations can strongly bias ATL06 elevations, resulting in under-estimation of the mass balance gradient with elevation. Based on these results, we conclude that ATL08 and ATL06 observations are best suited for characterization of watershed-scale snow depth and mass balance gradients over relatively shallow slopes with thick snowpacks . In these regions, ICESat-2 elevation residual-derived snow depth and mass balance transects can provide valuable watershed scale constraints on terrain parameter- and model-derived estimates of snow accumulation and melt.

Alaska, Idaho↗

USGS2021

The imminent calving retreat of Taku Glacier

Along the rugged Southeast Alaska coast, 30 kilometers northeast of the state capital Juneau, a tidewater glacier has largely defied global trends by steadily advancing for most of the past century while most glaciers on Earth retreated. This 55-kilometer-long and nearly 1,500-meter-thick tidewater glacier, named Taku Glacier, or T'aaḵú Ḵwáan Sít'i in the language of the Indigenous Tlingit people, has been the focus of continuous scientific study for more than 70 years. Some records even extend back to the mid-18th century. With this long observation record and the glacier’s year-round accessibility and proximity to Juneau and adjacent research facilities, Taku provides an unparalleled locale to study tidewater glaciers and their response to Earth’s rapidly changing climate.

Eos, American Geophysical Union↗

USGS2020

Explaining mass balance and retreat dichotomies at Taku and Lemon Creek Glaciers, Alaska

We reanalyzed mass balance records at Taku and Lemon Creek Glaciers to better understand the relative roles of hypsometry, local climate and dynamics as mass balance drivers. Over the 1946–2018 period, the cumulative mass balances diverged. Tidewater Taku Glacier advanced and gained mass at an average rate of +0.25±0.28 m w.e. a–1, contrasting with retreat and mass loss of –0.60±0.15 m w.e. a-1 at land-terminating Lemon Creek Glacier. The uniform influence of regional climate is demonstrated by strong correlations among annual mass balance and climate data. Regional warming trends forced similar statistically significant decreases in surface mass balance after 1989: –0.83 m w.e. a–1 at Taku Glacier and –0.81 m w.e. a–1 at Lemon Creek Glacier. Divergence in cumulative mass balance arises from differences in glacier hypsometry and local climate. Since 2013 negative mass balance and glacier-wide thinning prevailed at Taku Glacier. These changes initiated terminus retreat, which could increase dramatically if calving begins. The future mass balance trajectory of Taku Glacier hinges on dynamics, likely ending the historic dichotomy between these glaciers.

Alaska↗

USGS2019

Reanalysis of the U.S. Geological Survey Benchmark Glaciers: Long-term insight into climate forcing of glacier mass balance

Mountain glaciers integrate climate processes to provide an unmatched signal of regional climate forcing. However, extracting the climate signal via intercomparison of regional glacier mass balance records can be problematic when methods for extrapolating and calibrating direct glaciological measurements are mixed or inconsistent. To address this problem, we reanalyzed and compared long-term mass balance records from the U.S. Geological Survey Benchmark Glaciers. These five glaciers span maritime and continental climate regimes of the western United States and Alaska. Each glacier exhibits cumulative mass loss since the mid-20th century, with average rates ranging from –0.58 to –0.30 m water equivalent (w.e.) a-1. We produced a set of solutions using different extrapolation and calibration methods to inform uncertainty estimates, which range from 0.22–0.44 m w.e. a-1. Mass losses are primarily driven by increasing summer warming. Continentality exerts a stronger control on mass loss than latitude. Similar to elevation, topographic shading, snow redistribution, and glacier surface features often exert first-order control on mass balance. The reanalysis underscores the value of geodetic calibration to resolve mass balance magnitude, as well as the irreplaceable value of direct measurements in contributing to process-based understanding of glacier mass balance.

Journal of Glaciology↗

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↗

USGS2015

End-of-winter snow depth variability on glaciers in Alaska

A quantitative understanding of snow thickness and snow water equivalent (SWE) on glaciers is essential to a wide range of scientific and resource management topics. However, robust SWE estimates are observationally challenging, in part because SWE can vary abruptly over short distances in complex terrain due to interactions between topography and meteorological processes. In spring 2013, we measured snow accumulation on several glaciers around the Gulf of Alaska using both ground- and helicopter-based ground-penetrating radar surveys, complemented by extensive ground truth observations. We found that SWE can be highly variable (40% difference) over short spatial scales (tens to hundreds of meters), especially in the ablation zone where the underlying ice surfaces are typically rough. Elevation provides the dominant basin-scale influence on SWE, with gradients ranging from 115 to 400 mm/100 m. Regionally, total accumulation and the accumulation gradient are strongly controlled by a glacier's distance from the coastal moisture source. Multiple linear regressions, used to calculate distributed SWE fields, show that robust results require adequate sampling of the true distribution of multiple terrain parameters. Final SWE estimates (comparable to winter balances) show reasonable agreement with both the Parameter-elevation Relationships on Independent Slopes Model climate data set (9–36% difference) and the U.S. Geological Survey Alaska Benchmark Glaciers (6–36% difference). All the glaciers in our study exhibit substantial sensitivity to changing snow-rain fractions, regardless of their location in a coastal or continental climate. While process-based SWE projections remain elusive, the collection of ground-penetrating radar (GPR)-derived data sets provides a greatly enhanced perspective on the spatial distribution of SWE and will pave the way for future work that may eventually allow such projections.

Alaska↗