Surface velocities of Upsala glacier, Southern Patagonian Andes, estimated using cross-correlation satellite imagery: 2013-2014 period

: In this study we present surface velocities estimation for the Upsala glacier catchment, South Patagonian Ice Field (SPI) during the summer season of years 2013 (January-March) and 2014 (March-April), including the Bertacchi, Cono, and Murallón tributaries using satellite images from Advanced Spaceborne ermal Emission and Reﬂection Radiometer (ASTER). e Cross-Correlation method was applied by COSI-Corr technique with sub-pixel accuracy. In general, it should be noted that the SPI glaciers, and Upsala glacier in particular, are fast-ﬂowing ice bodies, which makes the technique works properly. Results of surface velocities estimation ranged from 0.22 to 2.93 md-1 for January-March 2013 and 0.12 to 5.8 md-1 for March-April 2014. In summary, COSI-Corr can achieved accurate and reliable results for glacier displacements and surface velocities estimation, also contributing in the better knowledge of the velocities change processes in time, taking into account Upsala is one of the most dynamic temperate glaciers of the SPI.


Introduction
Glaciers worldwide stand out as one of the best climate change indicators, such as the fluctuations like a function of height (Oerlemans, 2005), and the calving glaciers are particularly sensitive due to the loss of mass in the terminus (Badino and Romeo, 2005).In recent decades, the globally widespread trend of glacier area loss has accelerated (Rignot and Kanagaratnam, 2006;Ding et al., 2006;Bolch et al., 2012;Zemp et al., 2015;Solomina et al., 2016;Yang et al., 2016;Ragettli et al., 2016) and the dynamics of glaciers has consequently readjusted (Bolch et al., 2012).In the Southern Andes, the Southern Patagonian Ice Field (SPI) constitute the largest masses of glacier ice in the Southern Hemisphere (Aniya, 2013).Glacier thinning in the SPI has sped up since 1968 (Rignot et al., 2003;Glasser et al., 2011;Willis et al., 2012), thus making an increasing Table 1 e main goal of the present study is to estimate surface velocities for Upsala glacier and its tributaries Bertacchi, Cono and Murallón as well as the update of the last period studied of Upsala glacier's surface speeds to 2014, since detailed information was available until 2011.We use a temporal series of optical satellite images acquired by the Advanced Spaceborne ermal Emission and Reflection Radiometer (ASTER) sensor (15 m resolution, www.lpdaac.usgs.gov (last visit 22-03-2016)) during the summer periods of the years 2013 (January-March) and 2014 (March-April).e cross-correlation technique with COSI-Corr soware was applied for this purpose.e investigation is a continuation of our previous effort (Moragues et al., 2016).Additionally, a review of the previously published literature on surface velocities of Upsala glacier for the last seventy years, has been used for reconstructing the recent evolution of glacier velocities.

Study area
e study area is located at the eastern margin and latitudinally centered in the SPI (Fig. 1a), and encompasses the Upsala glacier tongue (49°50'-73°15') and its major tributaries (Bertacchi, Cono and Murallón glaciers).e main Upsala outlet glacier and the aforementioned tributaries constitute a compound basin with an area of 838 km2 (calculated from Landsat satellite images, November 2016).e main (western) glacier tongue calves into Argentino Lake and a secondary (eastern) tongue discharges in the Guillermo proglacial lake (Fig. 1b) (Skvarca, 2002).e Holocene recession of Upsala glacier exposed Brazo Cristina more than 10,115±100 calyrs BP. e Upsala glacier readvanced at least seven times, the first being a relatively minor expansion documented only in stratigraphic sections between 7,730±50 and 7,210±45 calyrs BP (Strelin et al., 2014).e recent glacier retreat has exposed a number of proglacial features, such as striated pavements, vast moraine deposits, and rochemoutonnées, among others.
Climate in this region is characterized by the westerly flow that carries humidity from the Pacific Ocean.is system hits the area more pronouncedly during the austral winter months (July-September) (Garreaud, 2009), producing precipitation maxima of up to 10,000 mm over the Southern Patagonian Icefield (Villalba et al., 2003;Garreaud et al., 2013).Due to the strong W-E orographic precipitation gradient, precipitation over the eastern slopes of the Argentinean Andes decrease to less than 300 mm within 100 km of the main Andean axis (Villalba et al., 2003).

Data Sets
e present study is based on four 15 m spatial resolution ASTER images (3N band) (Table 2).e ASTER images were obtained from the United States Geological Survey (USGS) Earth Explorer portal, and there were acquired during the ablation period years 2013 and 2014 (January-April), contain no significant seasonal snow, and are devoid of clouds.Note that time separation between processed scenes is crucial for the calculation of surface velocities and their uncertainties (Berthier et al., 2005), i.e., separation in days must be long enough to show the velocity signal in the slower sectors of the glacier, but short enough to preserve the morphological features that are tracked down in the correlation process.
e Shuttle Radar Topography Mission (SRTM), carrying interferometric synthetic aperture radar (SAR) sensor in C-band, acquired altitudinal data between 11 and22 February 2000 (www.jpl.nasa.gov/srtm, (last visit 16-04-2016)).We used the SRTM DEM with processing level of 1 Arc-Second, 30 m resolution, obtained from the USGS, data pool to resolve the orthorectification processes.To provide an optimal coregistration between the ASTER imagery and the SRTM (Ayoub et al., 2015), all data was projected to the UTM zone 18S projection, and re-sampled to 30 m spatial resolution by bilinear interpolator.
In addition, Corona, Hexagon and Landsat images were used to map the historical reconstruction of the frontal positions since 1968 until 2016 (Table 2).e images were selected according to their availability in the archive.Corona and Hexagon scenes (10 m resolution) (Fowler, 2013) were geometrically treated to eliminate the produced deformations and were georeferenced through the incorporation of ground control points (GCPs) from Landsat images, and re-sampled to the same spatial resolution of Landsat (Masiokas et al., 2015).e horizontal average errors were 1.2 pixels.To more detail of the overall processing, see Hollingsworth et al. (2012).All of the process described below was done using SOPI GIS soware.We use a manual delineation to detect the changes on the front positions (Nuth et al., 2013).An axial profile (A) was drawn trough the main flux of Upsala glacier to monitor the fluctuations.e origin (zero) of the axis was taken in the high zone of the central flow.

ASTER images pair cross correlation
In the first stage of COSI-Corr soware, geometric artifacts such as topographic were removed by the images orthorectification to mitigate those influences in the final processes (Ayoub et al., 2015).en, the image orthorectification processes were carried out selecting 30 tie-points (X,Y coordinates) over stable areas (mainly rock outcrops), areas outside the glaciers or with water, shade and snow cover.e tie points are defined either by taking points between the image and a georeferenced referenceimage (i.e., orthorectified image, shaded DEM, high resolution digitized map, etc.) (Ayoub et al., 2015).us, in this work tie points were taken between the reference image (January, 2013) and the shaded SRTM DEM.ese tie-points are optimized by sub-pixel correlation and converted into Ground Control Points (GCPs) because manual tie points selection leads to unavoidable mis-registrations between the orthorectified referenceimage and the orthorectified slaveimage.However, these mis-registrations can be reduced aerwards by optimizing the image with slave GCPs (Ayoub et al., 2009).Because GCPs have known ground coordinates, as they improve the georeferencing of the referenceand slave images.It can be assumed that GCPs can be located on the ASTER images with a±1 pixel accuracy (Berthier et al., 2005).Once the images were orthorectified, the scenes were spatially trimmed to the Upsala glacier area, tributary glaciers and surrounding terrain to avoid the presence of clouds and streamline the process.
In the subsequent, all of the ASTER images were co-registered to the reference2013 image.In a similar process to the previous one, 30 homologous points were selected and optimized between the reference2013 image and the three remaining ASTER scenes.is step assumes that the referenceimage georeferencing as faultless.Aer this, the referenceand slaveimages should be co-registered.
During this stage, the ASTER scenes were correlated to determine the offset in Δx and Δy coordinates.Once the multi-temporal images were orthorectified, the cross-correlation was performed in a region of interest (ROI) containing the glacier area and the immediate surroundings.e correlation was carried out in a two-step process.e first step estimates the image offset and in the second one, the offsets were obtained with sub-pixels accuracy (Ayoub et al., 2009), in this case less than 15 m.e horizontal displacements were calculated implementing the frequency correlator based on the Fourier method that is more sensitive to noise, and is therefore recommended for good quality optical images as those in the present study (Leprince et al., 2007;Ayoub et al., 2015).In the Fourier domain, the modulus of a white noise remains constant, and assuming that the images are degraded with some additive white noise, the phase information is then most likely to be biased in the high frequencies (Leprince et al., 2007a).
Low coherence values, caused by cloud cover or seasonal snow, may cause image de-correlation.Before extraction of ground displacements, areas with de-correlation were manually removed.An example of this can be found at Upsala glacier front in 2013 (Fig. 2c).e gaps were filled using the nearest neighbor interpolator.Nearest Neighbor interpolation finds the closest subset of input samples to a query point and applies weights to them based on proportionate areas (Sibson, 1981).

Extraction of glacier displacements
To obtain the values of large glacier displacements in the velocity field, multiple scale correlation analysis was performed with an initial search window size of 64 to 32 pixels, with steps of 8 pixels.e COSI-Corr cross correlation process output consists in the horizontal Δx (E-W) and Δy (N-S) displacements, and an additional SNR (Signal to Noise Ratio) one that numerically evaluates the quality of the measurement.Note that in a range of 0 to 1, a higher SNR indicates a better correlation.In some areas, the correlation is lost, then SNR reach very small atypical values close to 0 (Ayoub et al., 2009).e resultant glacier displacements are obtained by using e daily velocity field calculation (in md-1) was obtained by dividing each of the displacement values by the number of days elapsed between the referenceand slave images (Fig. 3).
Axial and transversal profiles were identified to extract the velocities field in the Upsala, Bertacchi, Cono, and Murallón glaciers (Fig. 4).In the axial profiles drawn on the tributaries (Bertacchi glacier, AB, Cono, AC and Murallón, AM), the origin (zero) was taken from the accumulation zone to the frontal zone.In the case of Upsala (AU), the origin was taken according to the border of the selected images.e transverse ones in Upsala (TU1, TU2) were traced from the West margin to the East margin.e cross-wises in the tributaries (Bertacchi glacier, TB1, TB2, Cone, TC1, TC2 and Murallón, TM1, TM2) were drawn in two areas, one in the lower part, near the front; and the other in the upper part on changes of topography, traced from South to North.In all cases the distances were expressed in meters.

Uncertainty estimation
e cross correlation is limited to a certain degree by the surface changes due to the different illumination/ incidence angles among the different satellite images (Berthier et al., 2005).is cause changes in the shadow extent and apparent displacements in shadow areas between the referenceand slave scenes, unrelated to glacier flow.
e random error (σr) was calculated on stable terrain and we choose three test areas (TA1, TA2, TA3) (Fig. 1b), where no terrain displacements is to be expected (Berthier et al., 2005;Scherler et al., 2008).Noise in displacements may resemble the Additive White Gaussian Noise.Since x and y are independent variables and follow a normal distribution.e value of the modules obeys the Rayleigh distribution and have a variance and a mean (Herman et al., 2011).For the 2013 calculation, 644 pixels were considered, whereas a total of 695 were utilized for the 2014 computation.3a-b).Also, SNR assessing the quality of the measure shows the areas with low correlation in the terminal portion of Upsala glacier in 2013 with values between 0 and 0.2 due clouds covering (Fig. 2c).
Considering the entire Upsala basin as a whole, glacier ice presents surface velocities with a range of 0.22 to 2.93±0.06md-1 in 2013 (Fig. 3a), and reached values of 0.12 to 5.8±0.1 md-1 in 2014 (Fig. 3b).e lower velocities correspond to the lateral margins of glaciers, whereas the maximum ones were always observed inthe middle part of the Upsala glacier.e pattern of glacier direction velocities corresponds to the orientation of the ice flow (Fig. 3a-b).us, it can be easily noted that Upsala glacier shows an increase in velocity at a faster rate in March-April 2014 (Fig. 3b) compared to January-March 2013 (Fig. 3a).Also, the main Upsala glacier flows at higher speed than any of its tributaries.e highest velocities were recorded towards the front of Upsala glacier, coinciding with the main N-S flow direction, reaching flow rates of 2.93 md-1 in 2013, where the values can be observed in a small set of pixels of the central zone of the glacier.Due to de-correlation, the flow acceleration at the glacier front cannot be valued in figure 3a.e maximum surface velocities yielded 5.8 md-1 in 2014; the area registering these values is wider than in 2013 and represents an increase of 25% from 2013 to 2014.
e mean velocities recorded across the axial profile AU were 1.87 md-1 and 2.61 md-1 for 2013 and 2014, respectively (Table 3).In 2013, velocities across the transversal profiles of the glacier differ to a lesser degree, with a maximum of 2.3 md-1 (TU1) and 1.7 md-1 (TU2), respectively.On the other hand, in 2014, transversal TU1 and TU2 profiles also show velocity differences between the outer margin of the glacier and the middle zone, where surface velocity reached a maximum of 2.5 md-1 (TU1) and 1.8 md-1 (TU2).Figure 4b shows the axial and transversal profiles on Upsala glacier for the 2013-2014 period, blue color line corresponds to year 2013 and the green color to 2014, respectively.Note that data from the 2013 image the pixels with de-correlation were eliminated in the frontal area.e 2013 velocity variability is due to the presence of clouds on the glacier (Fig. 4b).
On average, and similarly to Upsala glacier, both the Bertacchi (1.73 md-1 and 2.33 md-1) and Cono (1.07 md-1 and 1.12 md-1) glaciers show higher speeds in 2014 compared to 2013, whereas Murallón glacier slowed down from 0.4 md-1 to 0.16 md-1 (Table 3).e second fastest glacier in the Upsala catchment was Bertacchi with maxima of 2.8 md-1 (2013) and 4.9 md-1 (2014), followed by Cono (1.6 md-1 and 1.8 md-1) and Murallón (1 md-1 and 2.2 md-1).Interestingly, despite the fact that axial profiles show increased velocities in 2014, the transversal profiles AB, AC and AM at the tributary glaciers terminus show only minor variability in flow rates.On the contrary, increases in glacier velocity in these tributaries have concentrated in the upper portions of the ablation areas.

Discussion
Image quality is one of the factors limiting the accuracy of glacier surface speed estimation in a multitemporal analysis of satellite images.Satellite based optical sensors may oen have observation limiting factors, such as cloud cover, especially in mountain regions (Gleitsmann and Kappas, 2006) such as in the SPI region.In this sense, the 2013 ASTER pair, used in this study, had severe issues with areas of image de-correlation compared to the 2014 pair.We associated the low correlation areas with the presence of variable seasonal snow, clouds and changing illumination conditions between the images.Although this situation, it was possible to estimate accurate results of surface speed motion for the Upsala glacier in the 2013-2014 period.
As shown in table 1, the lowest surface velocities known for the Upsala glacier were estimated by Aniya and Skvarca (1992), between 1968 and 1970, with a velocity of 700 my-1 (1.9 md-1).Mouginot and Rignot (2015) presented a comprehensive study of Upsala glacier surface velocities in the NPI and SPI between 1993 and 2014, revealing a velocities increase between 1993 (0.30 md-1) and 2001 (2.75 md-1), but later decreased by 20% between 2001 and 2005.In the same study authors report that between 2005 and 2010, the glacier sped up yet again by 50% before decreasing again in 2014 (3.56 md-1).In addition, Sakakibara and Sugiyama (2014) found velocities higher than 3,000 my1 (8.22 md1) near the calving front for 2010 year, when most freshwater glaciers in the SPI had a general deceleration trend.
e acceleration of surface velocities found in the period 2013-2014 has been related to the calving processes and the interaction between the glacier itself and the Argentino Lake (Sugiyama et al., 2016).On calving glaciers, surface velocities reach a maximum at the glacier front, due to the pulling effect of high calving rates (Sakakibara and Sugiyama, 2014).Upsala glacier has been in fact described as one of the most representative glaciers in this regard (Sakakibara et al., 2013).e high calving rates may be further enhanced when near-buoyancy conditions are reached at the front of a glacier calving into deep waters (Rivera et al., 2012).In the case of Upsala, this is confirmed by the recently bathimetric survey of the lake, where in 2016 it was investigated until approximately 600 m near the front of the glacier (Sugiyama et al., 2016).e glacier retreat and surface velocity of large calving glaciers of the SPI are driven by the bedrock topography (Naruse and Skvarca, 2000;Skvarca et al., 2002), which in turn is related to lithology (Kraemer and Riccardi, 1997), and may be unrelated to the climate variability to a great extent (Yde and Paasche, 2010;Post et al., 2011;Wilson et al., 2016).e great depth of the Argentino Lake and the Upsala channel in particular might be related to the abundance of easily erodible rocks such as the slates of Pizarras Formation in the Upsala basin (Lo Vecchio et al., 2016).is noticeable acceleration of glaciers at their terminus has been previously described for other outlet glaciers in Patagonia such as Jorge Montt (Rivera et al., 2012;Sakakibara and Sugiyama, 2014;Bown, 2015), where calving is driven at least in part by water depths near the front.

Frontal position and surface velocities variation
ak Sakakibara and Sugiyama (2014) studied the glacier frontal variations of the outlet glaciers in the SPI.ese authors reported a general recession trend for most glaciers in the area, though spatially and temporally varied.e ice speed acceleration implied a substantial increase in frontal ablation (calving plus subaqueous frontal melting) (Sugiyama et al., 2016).is statistically consistent and significant recession suggests the existence of warming trend in the region.However, there is not yet enough statistical evidence of climate data within the ice field to state that a regional atmospheric warming is the cause of glacier retreat (Skvarca and De Angelis, 2002).
We examined the relationship between the frontal positions and variations of the surface speeds of the Upsala glacier (Fig. 5).For this figure, we took into consideration the different techniques and acquisition dates for each of the studies listed in table 1.To make the comparison easy, we refer to the glacier surface velocities as near as possible to the glacier terminus (e.g., Sakakibara et al., 2013;Muto and Furuya, 2013), and the original my-1 values were converted into md-1 rates.We detected at least two episodes where a surface velocity increase coincided with enhanced glacier retreat.e first one is observed for the 1993-2002 period, when the glacier reached velocities up to 4.4 md-1 in the year 1993 (Skvarca et al., 1995), whereas between 2000 and 2002 maximum velocities exceeded 4.7 md-1 (Sakakibara and Sugiyama, 2014).e surface velocity of Upsala glacier increased significantly starting in the year 2008, coinciding with the rapid retreat of the frontal position.is velocity increase (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008) coincides with an enhanced -90 my-1 retreat rate.On the other hand, during the 2009-2010 period, Sakakibara and Sugiyama (2014) reported a maximum glacier velocity of 9.3 md-1 (historical maximum registered so far for Upsala glacier), coinciding with an abrupt retreat rate of 1,000 my-1 approx.
Figure 5b shows a rough estimate of Upsala glacier terminus positions from early 20th century based on Warren et al.(1995Warren et al.( ) (1914Warren et al.( -1945 period) period), and additional mapping on CORONA (1963CORONA ( -1979) ) and Landsat imagery .At least three stages of varied retreat rates can be recognized.During the first stage , Upsala glacier retreated at -193 my-1.In the second stage , the glacier was fairly stable, retreating merely 30 my-1.roughout this second phase, Upsala glacier occupied the shallower part of the Upsala channel south of the Vacas valley, (Skvarca et al., 2002), where extensive rock exists.e combination of these two factors may have contributed to the stabilization of the glacier front during this period.Aer this relatively stable period, the glacier terminus begun a steady and fast retreat phase in 1978 (Naruse and Skvarca, 2000) that is presently ongoing, and which represented approximately 15% of the total mass loss of the ice field between February 2000 and August 2011 (Willis et al., 2012).is stage is characterized by retreat rates of -263 my-1, which are in agreement with those found for Jorge Montt glacier (-240 my-1) for the 2000-2012 period (Sakakibara and Sugiyama, 2014).In terms of glacier elevation changes, Willis et al. (2012) determined thinning rates of 24 my-1 at the glacier snout for the 2000-2012 period.Mean rates of frontal recoil of 400 my-1 (Skvarca et al., 2003), with estimated maxima higher than 800 my-1 (Naruse et al., 1997)

Bertacchi glacier retreat
Due to the fast retreat of Upsala glacier, the tributary Bertacchi tongue is no longer restrained at its terminus (Fig. 6).us, the glacier may be subject to longitudinal stress due to a new, faster flow rate, and this longitudinal stretching may lead to increased velocities in certain glacier areas.Since glacier movement responds to the viscous-elastic flow laws and flows in response to gravity, the glacier areas with steeper slopes should be the most affected by this process (Fig. 6).Higher glacier velocities (4 md-1 in the upper ablation area and 5.5 md-1 at the terminus) coinciding with abrupt slope areas can be observed in figure 6. Presently, the terminus of Bertacchi glacier is now calving into Argentino Lake and may be undergoing similar glacierlake interactions as in the large outlet glaciers (Fig. 3d).

Conclusions
Upsala glacier is one of the temperate glaciers of high scientific interest in the SPI as it presents a combination of processes and mechanisms such as rapid frontal retreating, surface velocity variations, ablation and mass balance changes, high calving rates, thinning dynamic features, that make it a highly active glacier.From the late 20th century to the present, it had a high rate of frontal retraction and an acceleration of exceptional surface velocities, compared to other glaciers studies in the SPI.
e present study estimates surface velocities of Upsala glacier and tributary glaciers, by means of crosscorrelation of ASTER images with the COSI-Corr soware.e technique proved to be appropriate for the detection and estimation of glacier displacements and velocities, as it offered an operative sub-pixel methodology for the measurements of the horizontal flow, taking in consideration a fast-flowing glacier such as Upsala.In addition, the investigation makes a significant contribution to the understanding of the dynamic of Upsala glacier, since it adds new data to the previously published studies focusing on glacier velocity, extending the surveyed period until 2014.Results of surface velocities estimations ranged from 0.22 to 2.93 md-1 for january-march 2013 and 0.12 to 5.8 md-1 for march-april 2014.e available glacier velocities were put into context of glacier fluctuations since 1914 and a tight correspondence between glacier velocities and retreat rates was implied.
FIG. 1 FIG. 1. a. Study area located in the Southern Patagonian Ice Field (SPI), Santa Cruz, Argentina (yellow lines); b.Upsala glacier and tributary glaciers Bertacchi, Cono and Murallón and stable test areas (red squares).

FIG. 2 .
FIG. 2. FIG. 2. Net horizontal glacier displacements for January 25 to March 14, 2013 and March 17 to April 2, 2014.Negative values indicate the direction of the N-S flow (b), while the positive values correspond to the flow W-E (a).SNR values for low correlation tend to be 0 (c).Negative N-S values indicate northerly flow (b), whereas positive W-E values correspond to Easterly flow (a).SNR values of low correlation tend to 0 (c).

FIG. 3 FIG
FIG. 3 FIG.3. Flux and ice surface velocities, expressed in meter/day, on Upsala glacier and its tributaries, derived from the crosscorrelation method with ASTER images (a 2013; b 2014).(c-d) Front evolution of Bertacchi glacier in 2013 and 2014 where in 2014 (right) it shows the glacier front and its interaction with the lake.Vectors in black: velocities derived by cross-correlation.
Figure2aand b shows the horizontal displacements in meters of the glaciers in the Upsala catchment: the Upsala glacier moves from north to south towards Argentino Lake, whereas the tributaries have displacements from west to east.e predominance of N-S displacement values in Upsala glacier ranges from 25 m to 200 m (2013-2014), with the negative sign being the N-S direction of displacement.For 2013, the maximum values of displacements, between 50 m to 200 m, are located in the middle of the main flow with a separation of time between images of 48 days (Fig.2a), and for 2014 are situated near to the front (50-200 m) with a separation between scenes of 16 days (Fig.2b).For tributary glaciers, the direction of horizontal displacement to indicate positive values in both periods, ranging from 5 m to 200 m displaced in a W-E direction.In general, in the mentioned periods, the displacements on valley glaciers have maximum values at the center of the valley, with 50 m to 200 m (W-E) for both periods.However, for the direction N-S the displacements were 50-200 m in 2013, 25-50 m in 2014, decreasing to a minimum at the margins, 0-50 m (W-E) and 0 m to 15 m (N-S) (Fig.3a-b).Also, SNR assessing the quality of the measure shows the areas with low correlation in the terminal portion of Upsala glacier in 2013 with values between 0 and 0.2 due clouds covering (Fig.2c).Considering the entire Upsala basin as a whole, glacier ice presents surface velocities with a range of 0.22 to 2.93±0.06md-1 in 2013 (Fig.3a), and reached values of 0.12 to 5.8±0.1 md-1 in 2014 (Fig.3b).e lower velocities correspond to the lateral margins of glaciers, whereas the maximum ones were always observed inthe middle part of the Upsala glacier.e pattern of glacier direction velocities corresponds to the orientation of the ice flow (Fig.3a-b).us, it can be easily noted that Upsala glacier shows an increase in velocity at a faster rate in March-April 2014 (Fig.3b) compared to January-March 2013 (Fig.3a).Also, the main Upsala glacier flows at higher speed than any of its tributaries.e highest velocities were recorded towards the front of Upsala glacier, coinciding with the main N-S flow direction, reaching flow rates of 2.93 md-1 in 2013, where the values can be observed in a small set of pixels of the central zone of the glacier.Due to de-correlation, the flow acceleration at the glacier front cannot be valued in figure3a.e maximum surface velocities yielded 5.8 md-1 in 2014; the area registering these values is wider than in 2013 and represents an increase of 25% from 2013 to 2014.e mean velocities recorded across the axial profile AU were 1.87 md-1 and 2.61 md-1 for 2013 and 2014, respectively (Table3).In 2013, velocities across the transversal profiles of the glacier differ to a lesser degree, with a maximum of 2.3 md-1 (TU1) and 1.7 md-1 (TU2), respectively.On the other hand, in 2014, transversal TU1 and TU2 profiles also show velocity differences between the outer margin of the glacier and the middle zone, where surface velocity reached a maximum of 2.5 md-1 (TU1) and 1.8 md-1 (TU2).Figure4bshows the axial and transversal profiles on Upsala glacier for the 2013-2014 period, blue color line corresponds to year 2013 and the green color to 2014, respectively.Note that data from the 2013 image the pixels with de-correlation were eliminated in the frontal area.e 2013 velocity variability is due to the presence of clouds on the glacier (Fig.4b).On average, and similarly to Upsala glacier, both the Bertacchi (1.73 md-1 and 2.33 md-1) and Cono (1.07 md-1 and 1.12 md-1) glaciers show higher speeds in 2014 compared to 2013, whereas Murallón glacier slowed down from 0.4 md-1 to 0.16 md-1 (Table3).e second fastest glacier in the Upsala catchment was Bertacchi with maxima of 2.8 md-1 (2013) and 4.9 md-1 (2014), followed by Cono (1.6 md-1 and 1.8 md-1) and Murallón (1 md-1 and 2.2 md-1).Interestingly, despite the fact that axial profiles show increased velocities in 2014, the transversal profiles AB, AC and AM at the tributary glaciers terminus show only minor variability in flow rates.On the contrary, increases in glacier velocity in these tributaries have concentrated in the upper portions of the ablation areas.
FIG. 5 FIG. 5. a. Temporal evolution of Upsala glacier frontal position, from 1914 to 2016 and the registration of studies carried out by different authors from the end of the 20th century to the present day on the surface velocities of the frontal area or the nearest sector recorded, expressed in md-1 (upper); b.Front fluctuations of Upsala glacier since 1914 at present (lower), the profile (A) by which the glacier retreated was calculated in kilometers, from the accumulation zone towards the frontal zone.
FIG. 6. FIG. 6. Isolines surface velocities of the Bertacchi Glacier in the year 2014, expressed in md-1 (upper).Topographic profile of the bed rock glacier (black line), dotted lines delimit the areas with steeper slopes (average).Panoramic photograph of the front of Bertacchi glacier, in contact with the vertical front of the Upsala glacier to the right (lower).