We introduce the first catchment dataset for large sample studies in Chile. This dataset includes 516 catchments; it covers particularly wide latitude (17.8 to 55.0° S) and elevation (0 to 6993 m a.s.l.) ranges, and it relies on multiple data sources (including ground data, remote-sensed products and reanalyses) to characterise the hydroclimatic conditions and landscape of a region where in situ measurements are scarce. For each catchment, the dataset provides boundaries, daily streamflow records and basin-averaged daily time series of precipitation (from one national and three global datasets), maximum, minimum and mean temperatures, potential evapotranspiration (PET; from two datasets), and snow water equivalent. We calculated hydro-climatological indices using these time series, and leveraged diverse data sources to extract topographic, geological and land cover features. Relying on publicly available reservoirs and water rights data for the country, we estimated the degree of anthropic intervention within the catchments. To facilitate the use of this dataset and promote common standards in large sample studies, we computed most catchment attributes introduced by Addor et al. (2017) in their Catchment Attributes and MEteorology for Large-sample Studies (CAMELS) dataset, and added several others. We used the dataset presented here (named CAMELS-CL) to characterise regional variations in hydroclimatic conditions over Chile and to explore how basin behaviour is influenced by catchment attributes and water extractions. Further, CAMELS-CL enabled us to analyse biases and uncertainties in basin-wide precipitation and PET. The characterisation of catchment water balances revealed large discrepancies between precipitation products in arid regions and a systematic precipitation underestimation in headwater mountain catchments (high elevations and steep slopes) over humid regions. We evaluated PET products based on ground data and found a fairly good performance of both products in humid regions (r > 0.91) and lower correlation (r < 0.76) in hyper-arid regions. Further, the satellite-based PET showed a consistent overestimation of observation-based PET. Finally, we explored local anomalies in catchment response by analysing the relationship between hydrological signatures and an attribute characterising the level of anthropic interventions. We showed that larger anthropic interventions are correlated with lower than normal annual flows, runoff ratios, elasticity of runoff with respect to precipitation, and flashiness of runoff, especially in arid catchments. CAMELS-CL provides unprecedented information on catchments in a region largely underrepresented in large sample studies. This effort is part of an international initiative to create multi-national large sample datasets freely available for the community. CAMELS-CL can be visualised from http://camels.cr2.cl and downloaded from https://doi.pangaea.de/10.1594/PANGAEA.894885.

As glaciers adjust their size in response to climate variations, long-term changes in meltwater production can be expected, affecting the local availability of water resources. We investigate glacier runoff in the period 1955–2016 in the Maipo River basin (4843 km², 33.0–34.3^∘ S, 69.8–70.5^∘ W), in the semiarid Andes of Chile. The basin contains more than 800 glaciers, which cover 378 km² in total (inventoried in 2000). We model the mass balance and runoff contribution of 26 glaciers with the physically oriented and fully distributed TOPKAPI (Topographic Kinematic Approximation and Integration)-ETH glacio-hydrological model and extrapolate the results to the entire basin. TOPKAPI-ETH is run at a daily time step using several glaciological and meteorological datasets, and its results are evaluated against streamflow records, remotely sensed snow cover, and geodetic mass balances for the periods 1955–2000 and 2000–2013. Results show that in 1955–2016 glacier mass balance had a general decreasing trend as a basin average but also had differences between the main sub-catchments. Glacier volume decreased by one-fifth (from 18.6±4.5 to 14.9±2.9 km³). Runoff from the initially glacierized areas was 177±25 mm yr⁻¹ (16±7 % of the total contributions to the basin), but it shows a decreasing sequence of maxima, which can be linked to the interplay between a decrease in precipitation since the 1980s and the reduction of ice melt. Glaciers in the Maipo River basin will continue retreating because they are not in equilibrium with the current climate. In a hypothetical constant climate scenario, glacier volume would reduce to 81±38 % of the year 2000 volume, and glacier runoff would be 78±30 % of the 1955–2016 average. This would considerably decrease the drought mitigation capacity of the basin.

We present a field‐data rich modelling analysis to reconstruct the climatic forcing, glacier response, and runoff generation from a high‐elevation catchment in central Chile over the period 2000–2015 to provide insights into the differing contributions of debris‐covered and debris‐free glaciers under current and future changing climatic conditions. Model simulations with the physically based glacio‐hydrological model TOPKAPI‐ETH reveal a period of neutral or slightly positive mass balance between 2000 and 2010, followed by a transition to increasingly large annual mass losses, associated with a recent mega drought. Mass losses commence earlier, and are more severe, for a heavily debris‐covered glacier, most likely due to its strong dependence on snow avalanche accumulation, which has declined in recent years. Catchment runoff shows a marked decreasing trend over the study period, but with high interannual variability directly linked to winter snow accumulation, and high contribution from ice melt in dry periods and drought conditions. The study demonstrates the importance of incorporating local‐scale processes such as snow avalanche accumulation and spatially variable debris thickness, in understanding the responses of different glacier types to climate change. We highlight the increased dependency of runoff from high Andean catchments on the diminishing resource of glacier ice during dry years.

Glaciers and rock glaciers play an important role in the hydrology of the semi-arid northern Chile. Several studies show that glaciers have rapidly lost mass in response to climate change during the last decades. The response of rock glaciers to climate change in this region is, however, less known. In this study we use a combination of historical aerial photography, stereo satellite imagery, airborne lidar, and the Shuttle Radar Topography Mission (SRTM) DEM to report glacier changes for the Tapado Glacier–rock glacier complex from the 1950s to 2020 and to report mass balances for the glacier component of the complex, Tapado Glacier. Furthermore, we examine high-resolution elevation changes and surface velocities between 2012 and 2020 for 35 rock glaciers in the La Laguna catchment. Our results show how Tapado Glacier has shrunk by $-\mathrm{25.2}\pm \mathrm{4.6}$ % between 1956 and 2020, while the mass balance of Tapado Glacier has become steadily more negative, from being approximately in balance between 1956 and 1978 ($-\mathrm{0.04}\pm \mathrm{0.08}$ m w.e. a⁻¹) to showing increased losses between 2015 and 2020 ($-\mathrm{0.32}\pm \mathrm{0.08}$ m w.e. a⁻¹). Climatological (re-)analyses reveal a general increase in air temperature, decrease in humidity, and variable precipitation since the 1980s in the region. In particular, the severe droughts starting in 2010 resulted in a negative mass balance of $-\mathrm{0.54}\pm \mathrm{0.10}$ m w.e. a⁻¹ between 2012 and 2015. The rock glaciers within the La Laguna catchment show heterogenous changes, with some sections of landforms exhibiting pronounced elevation changes and surface velocities exceeding that of Tapado Glacier. This could be indicative of high ice contents within the landforms and also highlights the importance of considering how landforms can transition from more glacial landforms to more periglacial features under permafrost conditions. As such, we believe high-resolution (sub-metre) elevation changes and surface velocities are a useful first step for identifying ice-rich landforms.