skip to main content


Title: Quantifying Geodetic Mass Balance of the Northern and Southern Patagonian Icefields Since 1976
Southern Andean glaciers contribute substantially to global sea-level rise. Unfortunately, mass balance estimates prior to 2000 are limited, hindering our understanding of the evolution of glacier mass changes over time. Elevation changes over 1976/1979 to 2000 derived from historical KH-9 Hexagon imagery and NASADEM provide the basis for geodetic mass balance estimates for subsets of the Northern Patagonian Icefield (NPI) and the Southern Patagonian Icefield (SPI), extending current mass balance observations by ∼20 years. Geodetic mass balances were −0.63 ± 0.03 m w.e. yr −1 for 63% of the NPI and −0.33 ± 0.05 m w.e. yr −1 for 52% of the SPI glacierized areas for this historical period. We also extend previous estimates temporally by 25% using NASADEM and ASTER elevation trends for the period 2000 to 2020, and find geodetic mass balances of −0.86 ± 0.03 m w.e. yr −1 for 100% of the NPI and −1.23 ± 0.04 m w.e. yr −1 for 97% of the SPI glacierized areas. 2000–2020 aggregations for the same areas represented in the 1976/1979 to 2000 estimates are −0.78 ± 0.03 m w.e. yr −1 in the NPI and −0.80 ± 0.04 m w.e. yr −1 on the SPI. The significant difference in SPI geodetic mass balance in the modern period for 100% vs. 52% of the glacierized area suggests subsampling leads to significant biases in regional mass balance estimates. When we compare the same areas in each time period, the results highlight an acceleration of ice loss by a factor of 1.2 on the NPI and 2.4 on the SPI in the 21st century as compared to the 1976/1979 to 2000 period. While lake-terminating glaciers show the most significant increase in mass loss rate from 1976/1979–2000 to 2000–2020, mass balance trends are highly variable within glaciers of all terminus environments, which suggests that individual glacier sensitivity to climate change is dependent on a multitude of morphological and climatological factors.  more » « less
Award ID(s):
1853881 2035535 2218664
NSF-PAR ID:
10329085
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Frontiers in Earth Science
Volume:
10
ISSN:
2296-6463
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    With a unique biogeophysical signature relative to other freshwater sources, meltwater from glaciers plays a crucial role in the hydrological and ecological regime of high latitude coastal areas. Today, as glaciers worldwide exhibit persistent negative mass balance, glacier runoff is changing in both magnitude and timing, with potential downstream impacts on infrastructure, ecosystems, and ecosystem resources. However, runoff trends may be difficult to detect in coastal systems with large precipitation variability. Here, we use the coupled energy balance and water routing model SnowModel‐HydroFlow to examine changes in timing and magnitude of runoff from the western Juneau Icefield in Southeast Alaska between 1980 and 2016. We find that under sustained glacier mass loss (−0.57 ± 0.12 m w. e. a−1), several hydrological variables related to runoff show increasing trends. This includes annual and spring glacier ice melt volumes (+10% and +16% decade−1) which, because of higher proportions of precipitation, translate to smaller increases in glacier runoff (+3% and +7% decade−1) and total watershed runoff (+1.4% and +3% decade−1). These results suggest that the western Juneau Icefield watersheds are still in an increasing glacier runoff period prior to reaching “peak water.” In terms of timing, we find that maximum glacier ice melt is occurring earlier (2.5 days decade−1), indicating a change in the source and quality of freshwater being delivered downstream in the early summer. Our findings highlight that even in maritime climates with large precipitation variability, high latitude coastal watersheds are experiencing hydrological regime change driven by ongoing glacier mass loss.

     
    more » « less
  2. ABSTRACT Glaciers spanning large altitudinal ranges often experience different climatic regimes with elevation, creating challenges in acquiring mass-balance and climate observations that represent the entire glacier. We use mixed methods to reconstruct the 1991–2014 mass balance of the Kahiltna Glacier in Alaska, a large (503 km 2 ) glacier with one of the greatest elevation ranges globally (264–6108 m a.s.l.). We calibrate an enhanced temperature index model to glacier-wide mass balances from repeat laser altimetry and point observations, finding a mean net mass-balance rate of −0.74 mw.e. a −1 ( ± σ = 0.04, std dev. of the best-performing model simulations). Results are validated against mass changes from NASA's Gravity Recovery and Climate Experiment (GRACE) satellites, a novel approach at the individual glacier scale. Correlation is strong between the detrended model- and GRACE-derived mass change time series ( R 2 = 0.58 and p ≪ 0.001), and between summer ( R 2 = 0.69 and p = 0.003) and annual ( R 2 = 0.63 and p = 0.006) balances, lending greater confidence to our modeling results. We find poor correlation, however, between modeled glacier-wide balances and recent single-stake monitoring. Finally, we make recommendations for monitoring glaciers with extreme altitudinal ranges, including characterizing precipitation via snow radar profiling. 
    more » « less
  3. Abstract. During the concluding phase of the NASA OperationIceBridge (OIB), we successfully completed two airborne measurementcampaigns (in 2018 and 2021, respectively) using a compact S and C band radarinstalled on a Single Otter aircraft and collected data over Alaskanmountains, ice fields, and glaciers. This paper reports seasonal snow depthsderived from radar data. We found large variations in seasonalradar-inferred depths with multi-modal distributions assuming a constantrelative permittivity for snow equal to 1.89. About 34 % of the snowdepths observed in 2018 were between 3.2 and 4.2 m, and close to 30 % of thesnow depths observed in 2021 were between 2.5 and 3.5 m. We observed snowstrata in ice facies, combined percolation and wet-snow facies, and dry-snow facies fromradar data and identified the transition areas from wet-snow facies to icefacies for multiple glaciers based on the snow strata and radarbackscattering characteristics. Our analysis focuses on the measured strataof multiple years at the caldera of Mount Wrangell (K'elt'aeni) to estimate the localsnow accumulation rate. We developed a method for using our radar readingsof multi-year strata to constrain the uncertain parameters of interpretationmodels with the assumption that most of the snow layers detected by theradar at the caldera are annual accumulation layers. At a 2004 ice core and2005 temperature sensor tower site, the locally estimated average snowaccumulation rate is ∼2.89 m w.e. a−1 between the years2003 and 2021. Our estimate of the snow accumulation rate between 2005 and2006 is 2.82 m w.e. a−1, which matches closely to the 2.75 m w.e. a−1 inferred from independent ground-truth measurements made the sameyear. The snow accumulation rate between the years 2003 and 2021 also showeda linear increasing trend of 0.011 m w.e. a−2. This trend iscorroborated by comparisons with the surface mass balance (SMB) derived forthe same period from the regional atmospheric climate model MAR (ModèleAtmosphérique Régional). According to MAR data, which show anincrease of 0.86 ∘C in this area for the period of 2003–2021, thelinear upward trend is associated with the increase in snowfall and rainfallevents, which may be attributed to elevated global temperatures. Thefindings of this study confirmed the viability of our methodology, as wellas its underlying assumptions and interpretation models.

     
    more » « less
  4. Abstract

    The Patagonian slab window has been proposed to enhance the solid Earth response to ice mass load changes in the overlying Northern and Southern Patagonian Icefields (NPI and SPI, respectively). Here, we present the first regional seismic velocity model covering the entire north‐south extent of the slab window. A slow velocity anomaly in the uppermost mantle indicates warm mantle temperature, low viscosity, and possibly partial melt. Low velocities just below the Moho suggest that the lithospheric mantle has been thermally eroded over the youngest part of the slab window. The slowest part of the anomaly is north of 49°S, implying that the NPI and the northern SPI overlie lower viscosity mantle than the southern SPI. This comprehensive seismic mapping of the slab window provides key evidence supporting the previously hypothesized connection between post‐Little Ice Age anthropogenic ice mass loss and rapid geodetically observed glacial isostatic uplift (≥4 cm/yr).

     
    more » « less
  5. Himalayan glaciers supply meltwater to densely populated catchments in South Asia, and regional observations of glacier change over multiple decades are needed to understand climate drivers and assess resulting impacts on glacier-fed rivers. Here, we quantify changes in ice thickness during the intervals 1975–2000 and 2000–2016 across the Himalayas, using a set of digital elevation models derived from cold war–era spy satellite film and modern stereo satellite imagery. We observe consistent ice loss along the entire 2000-km transect for both intervals and find a doubling of the average loss rate during 2000–2016 [−0.43 ± 0.14 m w.e. year −1 (meters of water equivalent per year)] compared to 1975–2000 (−0.22 ± 0.13 m w.e. year −1 ). The similar magnitude and acceleration of ice loss across the Himalayas suggests a regionally coherent climate forcing, consistent with atmospheric warming and associated energy fluxes as the dominant drivers of glacier change. 
    more » « less