Velocities of Major Outlet Glaciers of the Patagonia Icefield Observed by TerraSAR-X

The capabilities of TerraSAR-X data for feature tracking by amplitude correlation over glacier surfaces are investigated. Methodical aspects of the amplitude correlation approach are described. The TerraSAR-X based velocity fields are compared with former InSAR derived velocities and field measurements on three outlet glaciers on the South Patagonia ice field

Distributed ice thickness and glacier volume in southern South America

South American glaciers, including those in Patagonia, presently contribute the largest amount of meltwater to sea level rise per unit glacier area in the world. Yet understanding of the mechanisms behind the associated glacier mass balance changes remains unquantified partly because models are hindered by a lack of knowledge of subglacial topography. This study applied a perfect-plasticity model along glacier centre-lines to derive a first-order estimate of ice thickness and then interpolated these thickness estimates across glacier areas. This produced the first complete coverage of distributed ice thickness, bed topography and volume for 617 glaciers between 41°S and 55°S and in 24 major glacier regions. Maximum modelled ice thicknesses reach 1631 m ± 179 m in the South Patagonian Icefield (SPI), 1315 m ± 145 m in the North Patagonian Icefield (NPI) and 936 m ± 103 m in Cordillera Darwin. The total modelled volume of ice is 1234.6 km3 ± 246.8 km3 for the NPI, 4326.6 km3 ± 865.2 km3 for the SPI and 151.9 km3 ± 30.38 km3 for Cordillera Darwin. The total volume was modelled to be 5955 km3 ± 1191 km3, which equates to 5458.3 Gt ± 1091.6 Gt ice and to 15.08 mm ± 3.01 mm sea level equivalent (SLE). However, a total area of 655 km2 contains ice below sea level and there are 282 individual overdeepenings with a mean depth of 38 m and a total volume if filled with water to the brim of 102 km3. Adjusting the potential SLE for the ice volume below sea level and for the maximum potential storage of meltwater in these overdeepenings produces a maximum potential sea level rise (SLR) of 14.71 mm ± 2.94 mm. We provide a calculation of the present ice volume per major river catchment and we discuss likely changes to southern South America glaciers in the future. The ice thickness and subglacial topography modelled by this study will facilitate future studies of ice dynamics and glacier isostatic adjustment, and will be important for projecting water resources and glacier hazards

Quantifying mass balance processes on the Southern Patagonia Icefield

Artículo de publicación ISIWe present surface mass balance simulations of the Southern Patagonia Icefield (SPI) driven by downscaled reanalysis data. The simulations were evaluated and interpreted using geodetic mass balances, measured point balances and a complete velocity field of the icefield for spring 2004. The high measured accumulation of snow of up to 15.4 m w.e. yr−1 (meters water equivalent per year) as well as the high measured ablation of up to 11 m w.e. yr−1 is reproduced by the model. The overall modeled surface mass balance was positive and increasing during 1975–2011. Subtracting the surface mass balance from geodetic balances, calving fluxes were inferred. Mass losses of the SPI due to calving were strongly increasing from 1975–2000 to 2000– 2011 and higher than losses due to surface melt. Calving fluxes were inferred for the individual glacier catchments and compared to fluxes estimated from velocity data. Measurements of ice thickness and flow velocities at the glaciers’ front and spatially distributed accumulation measurements can help to reduce the uncertainties of the different terms in the mass balance of the Southern Patagonia Icefield.FONDECYT 3140135 European Union 22637

A decreasing glacier mass balance gradient from the edge of the Upper Tarim Basin to the Karakoram during 2000-2014

In contrast to the glacier mass losses observed at other locations around the world, some glaciers in the High Mountains of Asia appear to have gained mass in recent decades. However, changes in digital elevation models indicate that glaciers in Karakoram and Pamir have gained mass, while recent laser altimetry data indicate mass gain centred on West Kunlun. Here, we obtain results that are essentially consistent with those from altimetry, but with two-dimensional observations and higher resolution. We produced elevation models using radar interferometry applied to bistatic data gathered between 2011 and 2014 and compared them to a model produced from bistatic data collected in 2000. The glaciers in West Kunlun, Eastern Pamir and the northern part of Karakoram experienced a clear mass gain of 0.043 ± 0.078~0.363 ± 0.065 m w.e. yr−1. The Karakoram showed a near-stable mass balance in its western part (−0.020 ± 0.064 m w.e. yr−1), while the Eastern Karakoram showed mass loss (−0.101 ± 0.058 m w.e. yr−1). Significant positive glacier mass balances are noted along the edge of the Upper Tarim Basin and indicate a decreasing gradient from northeast to southwest

On the frontal ablation of Alaska tidewater glaciers

Thesis (Ph.D.) University of Alaska Fairbanks, 2013Sea level rise is a major problem that society will face in the coming century. One of the largest unknown components of sea level rise is frontal ablation (the sum of mass loss through calving and subaqueous melting) from glaciers and ice sheets. Using estimates of ice thickness, rates of glacier length change, and glacier velocities, we present a record of frontal ablation over the period 1985-2012 for 20 Alaska tidewater glaciers. We also present a new method for estimating ice thickness by solving the continuity equation between adjacent flowlines. Because of the wealth of data available, we apply this method to Columbia Glacier, Alaska. The mean ice thickness and volume of Columbia Glacier were approximately halved over the period 1957-2007, from 281 m to 143 m, and from 294 km�_ to 134 km�_, respectively. Using bedrock slope and considering how waves of thickness change propagate through the glacier, we conclude that the rapid portion of this tidewater glacier's retreat is likely nearing an end. We present a 64 year record of length change for 50 Alaska tidewater glaciers, over the period 1948-2012. Most (31) glaciers retreated over the period. Examination of the onset of glacier retreats indicates a correlation between high summer sea surface temperature and the triggering of retreat. Finally, we present a 27 year record of surface velocity for 20 Alaska tidewater glaciers derived from Landsat imagery. Surface velocities vary by as much as 80% throughout the year, indicating that using measurements from one time of year may bias estimates of frontal ablation. The total mean rate of frontal ablation for these 20 glaciers over the period 1985-2012 is 16.2 � 6.5 Gt a����_. Extending this to the remaining 30 Alaska tidewater glaciers yields an estimate of frontal ablation of 18.3 � 7.3 Gt a����_, approximately 50% of the climatic mass balance of the region. This indicates the important, non-negligible role frontal ablation can play in regional mass balance, even where tidewater glaciers cover a small fraction of the total area.Chapter 1. Introduction -- Chapter 2. Using surface velocities to calculate ice thickness and bed topography: a case study at Columbia Glacier, Alaska -- Chapter 3. Alaska tidewater glacier terminus positions, 1948-2012 -- Chapter 4. Alaska tidewater glacier velocities and frontal ablation, 1985-2012 -- Chapter 5. Conclusions

Crustal deformation across the Southern Patagonian Icefield observed by GNSS

Geodetic GNSS observations at 43 sites well distributed over the Southern Patagonian Icefield region yield site velocities with a mean accuracy of 1 mm/a and 6 mm/a for the horizontal and vertical components, respectively. These velocities are analyzed to reveal the magnitudes and patterns of vertical and horizontal present-day crustal deformation as well as their primary driving processes. The observed vertical velocities confirm a rapid uplift, with rates peaking at 41 mm/a, causally related to glacial-isostatic adjustment (GIA). They yield now an unambiguous preference between two competing GIA models. Remaining discrepancies between the preferred model and our observations point toward an effective upper mantle viscosity even lower than 1.6.10¹⁸ and effects of lateral rheological heterogeneities. An analysis of the horizontal strain and strain-rate fields reveals some complex superposition, with compression dominating in the west and extension in the east. This deformation field suggests significant contributions from three processes: GIA, a western interseismic tectonic deformation field related to plate subduction, and an extensional strain-rate field related to active Patagonian slab window tectonics.Facultad de Ciencias Astronómicas y Geofísica

Using a new generation of remote sensing to monitor Peru’s mountain glaciers

Remote sensing technologies are integral to monitoring mountain glaciers in a warming world. Tropical glaciers, of which around 70% are located in Peru, are particularly at risk as a result of climate warming. Satellite missions and field-based platforms have transformed understanding of the processes driving mountain glacier dynamics and the associated emergence of hazards (e.g. avalanches, floods, landslides), yet are seldom specialised to overcome the unique challenges of acquiring data in mountainous environments. A ‘new generation’ of remote sensing, marked by open access to powerful cloud computing and large datasets, high resolution satellite missions, and low-cost science-grade field sensors, looks to revolutionise the way we monitor the mountain cryosphere. In this thesis, three novel remote sensing techniques and their applicability towards monitoring the glaciers of the Peruvian Cordillera Vilcanota are examined. Using novel processing chains and image archives generated by the ASTER satellite, the first mass balance estimate of the Cordillera Vilcanota is calculated (-0.48 ± 0.07 m w.e. yr-1) and ELA change of up to 32.8 m per decade in the neighbouring Cordillera Vilcabamba is quantified. The performance of new satellite altimetry missions, Sentinel-3 and ICESat-2, are assessed, with the tracking mode of Sentinel-3 being a key limitation of the potential for its use over mountain environments. Although currently limited in its ability to extract widespread mass balance measurements over mountain glaciers, other applications for ICESat-2 in long-term monitoring of mountain glaciers include quantifying surface elevation change, identifying large accumulation events, and monitoring lake bathymetry. Finally, a novel low-cost method of performing timelapse photogrammetry using Raspberry Pi camera sensors is created and compared to 3D models generated by a UAV. Mean difference between the Raspberry Pi and UAV sensors is 0.31 ± 0.74 m, giving promise to the use of these sensors for long-term monitoring of recession and short-term warning of hazards at glacier calving fronts. Together, this ‘new generation’ of remote sensing looks to provide new glaciological insights and opportunities for regular monitoring of data-scarce mountainous regions. The techniques discussed in this thesis could benefit communities and societal programmes in rapidly deglaciating environments, including across the Cordillera Vilcanota

The glaciers climate change initiative: Methods for creating glacier area, elevation change and velocity products

Glaciers and their changes through time are increasingly obtained from a wide range of satellite sensors. Due to the often remote location of glaciers in inaccessible and high-mountain terrain, satellite observations frequently provide the only available measurements. Furthermore, satellite data provide observations of glacier character- istics that are difficult to monitor using ground-based measurements, thus complementing the latter. In the Glaciers_cci project of the European Space Agency (ESA), three of these characteristics are investigated in detail: glacier area, elevation change and surface velocity. We use (a) data from optical sensors to derive glacier outlines, (b) digital elevation models from at least two points in time, (c) repeat altimetry for determining elevation changes, and (d) data from repeat optical and microwave sensors for calculating surface velocity. For the latter, the two sensor types provide complementary information in terms of spatio-temporal coverage. While (c) and (d) can be generated mostly automatically, (a) and (b) require the intervention of an analyst. Largely based on the results of various round robin experiments (multi-analyst benchmark studies) for each of the products, we suggest and describe the most suitable algorithms for product creation and provide recommendations concerning their practical implementation and the required post-processing. For some of the products (area, velocity) post-processing can influence product quality more than the main-processing algorithm

Spatio-temporal variability in Southern Hemisphere glacier snowline altitudes from 2000-2020

The glacierised Southern Hemisphere is vulnerable to continued shrinkage under climate change, but representation of these mountainous regions in climate research is limited by hemispheric and altitudinal scarcity of meteorological observations. End-of-summer snowline altitude (SLAEOS) indicates glacier response to climatic forcing, though has been estimated with low spatio-temporal coverage for the Southern Hemisphere. This study presents the first Southern Hemisphere-wide quantification of SLAEOS, with analysis of regional and intra-regional trends. An automated approach was implemented in Google Earth Engine, in which glacier snow cover was classified in Landsat scenes using Otsu image segmentation and SLAEOS was estimated as the lowest altitude from which snow cover ratio was continuously > 0.5. Results encompassed 6485 glaciers of the Southern Alps, Andes, and Antarctic Peninsula, with trends calculated from 2000-2020. Snowlines underwent widespread retreat in this period; mean rates of SLAEOS rise were between 2.19 and 6.28 m yr-1 for regions, between 1.63 and 7.55 m yr-1 for east/west sub-regions, and were mostly accelerated for the recent decade (2010-2020). Mean SLAEOS lowering (-30 to -1 m yr-1) indicated stability in the southernmost Andes, contrasting to rapid SLAEOS rise (10 to 30 m yr-1) in the southern Central Chilean Andes, and eastern slopes generally experienced increased rates of SLAEOS rise compared to western slopes. SLAEOS variability was reflected in periods of summer warming and reductions in summer snowfall, though correlation with these variables was not consistently identified. East-west and north-south disparities in absolute SLAEOS and rates of SLAEOS change were linked to spatial variability in terrain elevation and prevailing moisture transport, with the latter evidencing the variability and impact of large-scale climatic modes. Given implications of observed trends for glacier mass loss, continued research may involve developing an annually-updated global dataset, investigating additional drivers of SLAEOS variability, and estimating glacier response times