Comparison of 3D surfaces from historical aerial images and UAV acquisitions to understand glacier dynamics: The Aneto glacier changes in 40 years

Pyrenean glaciers have shown a marked area and thickness decrease in the last century, especially in the last decades, and currently are highly threatened by climate change. Out of the 39 glaciers existing in the Pyrenees in 1984, 23 very small glaciers remain in this mountain range, from which only four have more than 10 ha. Probably, the most emblematic glacier of these four is Aneto glacier as it is located in the North-East face of the highest summit in the Pyrenees, the Aneto peak (3404 m a.s.l.). This work presents the Aneto glacier surface reconstruction from aerial images obtained in 1981, and its comparison with the glacier surface obtained in 2021 with Unmanned Aerial Vehicles (UAV) images.

The 1981 and 2021 images have been processed with Structure from Motion (SfM) algorithms to reconstruct the Digital Surface Model (DSM) of the glacier and nearby terrain. Taking advantage of the accurate geolocation of the UAV images in 2021 (GPS with RTK/PPK surveying), the DSM obtained has a precise representation of the glacier surface. Oppositely the aerial images of 1981 lack precise geolocation and thus require a post-processing analysis. The aerial images of the '80s have been firstly geolocated with Ground Control Points (GCPs) of known coordinates within the study area (summits, crests, and rock blocks with unaltered position). After this initial geolocation, the DSM of 1981 was generated with SfM algorithms. Nevertheless, this DSM still lacks a geolocation accuracy. To allow a comparison between the 1981 and the 2021 DSMs, the glacier surface in 1981 was registered to the 2021 surface with an Iterative Close Point (ICP) routine in the surrounding area of the glacier. The technique described in this work may be applicable to other historical aerial images, which may allow studying glacier evolutions all over the world for dates without field observations.

The surface comparison generated with images that have a temporal difference of 40 years has shown the dramatic area and thickness loss of this glacier, with areas decreasing more than 68 m, and an average thickness reduction of 31.5 m. In this period, the glacier has reduced its extent by about a 60%. There is a recent acceleration in the rate of shrinkage if we compare these data with the obtained for the period 2011-2021, in which area loss reaches 15% and thickness reduction almost reaches 10 m. During the 1981-2021 period the shrinkage rate is 0.78 m thickness/year and 1.5% area/year, meanwhile, during the 2011-2021 period the shrinkage rate is 0.99 m thickness/year and 2.7% area/year.