Glacier melt is an important source of water for high Andean rivers in central Chile, especially in dry years, when it can be an important contributor to flows during late summer and autumn. However, few studies have quantified glacier melt contribution to streamflow in this region. To address this shortcoming, we present an analysis of meteorological conditions and ablation for Universidad Glacier, one of the largest valley glaciers in the central Andes of Chile at the head of the Tinguiririca River, for the 2009–2010 ablation season. We used meteorological measurements from two automatic weather stations installed on the glacier to drive a distributed temperature-index and runoff routing model. The temperature-index model was calibrated at the lower weather station site and showed good agreement with melt estimates from an ablation stake and sonic ranger, and with a physically based energy balance model. Total modelled glacier melt is compared with river flow measurements at three sites located between 0.5 and 50 km downstream. Universidad Glacier shows extremely high melt rates over the ablation season which may exceed 10 m water equivalent in the lower ablation area, representing between 10 and 13 % of the mean monthly streamflow at the outlet of the Tinguiririca River Basin between December 2009 and March 2010. This contribution rises to a monthly maximum of almost 20 % in March 2010, demonstrating the importance of glacier runoff to streamflow, particularly in dry years such as 2009–2010. The temperature-index approach benefits from the availability of on-glacier meteorological data, enabling the calculation of the local hourly variable lapse rate, and is suited to high melt regimes, but would not be easily applicable to glaciers further north in Chile where sublimation is more significant.

Recent evidence shows that most Patagonian glaciers are receding rapidly. Due to the lack of in situ long-term meteorological observations, the understanding of how glaciers are responding to changes in climate over this region is extremely limited, and uncertainties exist in the glacier surface mass balance model parameterizations. This precludes a robust assessment of glacier response to current and projected climate change. An issue of central concern is the accurate estimation of precipitation phase. In this work, we have assessed spatial and temporal patterns in snow accumulation in both the North Patagonia Icefield (NPI) and South Patagonia Icefield (SPI). We used a regional climate model, RegCM4.6 and four Phase Partitioning Methods (PPM) in addition to short-term snow accumulation observations using ultrasonic depth gauges (UDG). Snow accumulation shows that rates are higher on the west side relative to the east side for both icefields. The values depend on the PPM used and reach a mean difference of 1,500 mm w.e., with some areas reaching differences higher than 3,500 mm w.e. These differences could lead to divergent mass balance estimations depending on the scheme used to define the snow accumulation. Good agreement is found in comparing UDG observations with modeled data on the plateau area of the SPI during a short time period; however, there are important differences between rates of snow accumulation determined in this work and previous estimations using ice core data at annual scale. Significant positive trends are mainly present in the autumn season on the west side of the SPI, while on the east side, significant negative trends in autumn were observed. Overall, for the rest of the area and during other seasons, no significant changes can be determined. In addition, glaciers with positive and stable elevation and frontal changes determined by previous works are related to areas where snow accumulation has increased during the period 2000–2015. This suggests that increases in snow accumulation are attenuating the response of some Patagonian glaciers to warming in a regional context of overall glacier retreat.

Glaciar Jorge Montt (48°20' S/73°30' W), one of the main tidewater glaciers of the Southern Patagonian Icefield (SPI), has experienced the greatest terminal retreat observed in Patagonia during the past century, with a recession of 19.5 km between 1898 and 2011. This retreat has revealed trees laying subglacially until 2003. These trees were dated using radiocarbon, yielding burial ages between 460 and 250 cal yrs BP. The presence of old growth forest during those dates indicates that Glaciar Jorge Montt was upvalley of its present position before the commonly recognized Little Ice Age (LIA) period in Patagonia. The post-LIA retreat was most likely triggered by climatically induced changes during the 20th century; however, Glaciar Jorge Montt has responded more dramatically than its neighbours. The retreat of Jorge Montt opened a 19.5 km long fjord since 1898, which reaches depths in excess of 390 m. The bathymetry is well correlated with glacier retreat rates, suggesting that dynamic responses of the glacier are at least partially connected to near buoyancy conditions at the ice front, resulting in high calving fluxes, accelerating thinning rates and rapid ice velocities.

Tidewater calving glaciers can undergo large fluctuations not necessarily in direct response to climate, but rather owing to complex ice–water interactions at the glacier termini. One example of this process in Chilean Patagonia is Glaciar Jorge Montt, where two cameras were installed in February 2010, collecting up to four glacier photographs per day, until they were recovered on 22 January 2011. Ice velocities were derived from feature tracking of the geo-referenced photos, yielding a mean value of 13 ±4 md^–1 for the whole lower part of the glacier. These velocities were compared to satellite-imagery-derived feature tracking obtained in February 2010, resulting in similar values. During the operational period of the cameras, the glacier continued to retreat (1 km), experiencing one of the highest calving fluxes ever recorded in Patagonia (2.4 km³ a^–1). Comparison with previous data also revealed ice acceleration in recent years. These very high velocities are clearly a response to enhanced glacier calving activity into a deep water fjord.

The Northern Patagonian Icefield (NPI) and the Southern Patagonian Icefield (SPI) have increased their ice mass loss in recent decades. In view of the impacts of glacier shrinkage in Patagonia, an assessment of the potential future surface mass balance (SMB) of the icefields is critical. We seek to provide this assessment by modelling the SMB between 1976 and 2050 for both icefields, using regional climate model data (RegCM4.6) and a range of emission scenarios. For the NPI, reductions between 1.5 m w.e. (RCP2.6) and 1.9 m w.e. (RCP8.5) were estimated in the mean SMB during the period 2005–2050 compared to the historical period (1976–2005). For the SPI, the estimated reductions were between 1.1 m w.e. (RCP2.6) and 1.5 m w.e. (RCP8.5). Recently frontal ablation estimates suggest that mean SMB in the SPI is positively biased by 1.5 m w.e., probably due to accumulation overestimation. If it is assumed that frontal ablation rates of the recent past will continue, ice loss and sea-level rise contribution will increase. The trend towards lower SMB is mostly explained by an increase in surface melt. Positive ice loss feedbacks linked to increasing in meltwater availability are expected for calving glaciers.

Glaciers in the central Andes of Chile are fundamental freshwater sources for ecosystems and communities. Overall, glaciers in this region have shown continuous recession and down-wasting, but long-term glacier mass balance studies providing precise estimates of these changes are scarce. Here, we present the first long-term (1955–2013/2015), region-specific glacier elevation and mass change estimates for the Maipo River Basin, from which the densely populated metropolitan region of Chile obtains most of its freshwater supply. We calculated glacier elevation and mass changes using historical topographic maps, Shuttle Radar Topography Mission (SRTM), TerraSAR-X add-on for Digital Elevation Measurements (TanDEM-X), and airborne Light Detection and Ranging (LiDAR) digital elevation models. The results indicated a mean regional glacier mass balance of −0.12 ± 0.06 m w.e.a⁻¹, with a total mass loss of 2.43 ± 0.26 Gt for the Maipo River Basin between 1955–2013. The most negative glacier mass balance was the Olivares sub-basin, with a mean value of −0.29 ± 0.07 m w.e.a⁻¹. We observed spatially heterogeneous glacier elevation and mass changes between 1955 and 2000, and more negative values between 2000 and 2013, with an acceleration in ice thinning rates starting in 2010, which coincides with the severe drought. Our results provide key information to improve glaciological and hydrological projections in a region where water resources are under pressure.