Debris Flows Occurrence in the Semiarid Central Andes under Climate Change Scenario

: This review paper compiles research related to debris flows and hyperconcentrated flows in the central Andes (30°–33° S), updating the knowledge of these phenomena in this semiarid region. Continuous records of these phenomena are lacking through the Andean region; intense precipitations, sudden snowmelt, increased temperatures on high relief mountain areas, and permafrost degradation are related to violent flow discharges. Documented catastrophic consequences related to these geoclimatic events highlight the need to improve their understanding in order to prepare the Andean communities for this latent danger. An amplified impact is expected not only due to environmental changes potentially linked to climate change but also due to rising exposure linked to urban expansion toward more susceptible or unstable areas. This review highlights as well the need for the implementation of preventive measures to reduce the negative impacts and vulnerability of the Andean communities in the global warming context.


Introduction
As in other high mountain ranges around the world, landslide processes are widespread in the Central Andes of Chile and Argentina [1][2][3][4][5]. Geomorphological studies have described paleo-landslides in the Central Andes [6][7][8][9]; however, studies of the most dangerous and impacting events such as debris flows are rare. These sudden events have generated many fatalities, destruction along international transportation corridors, road damages, house destructions, traffic interruptions, infrastructures damages, and complications on drinking water supply with underestimated regional economical loses [10][11][12]. Despite the high economic costs, the hydraulic characterization of these destructive events is not well known, thus preventing the adoption of adequate structural preventive measures for mitigating risks. At the same time, non-systematic monitoring of these phenomena impedes evaluation on the recurrence or frequency of events. The focus of attention on emergencies instead of preparedness or the adoption of preventive measures is the most common approach. Landslide hazard and risk assessment studies are starting to gradually be included in land-use planning of government offices.
Current climate trends and future projections show in addition to temperature rise a decrease in winter rainfall [13], increased summer rainfall (this is better documented in the Argentine side), and a rise in the torrential nature of storms [14]. Snow precipitation has been reduced since 2010 on the Central Andes [15,16]. This period of 10 consecutive dry years has been called around the region the Mega Drought [16]. All these natural environmental conditions are correlated with landslide activity and the occurrence of debris flows of the Central Andes, showing a positive temporal trend [7,[17][18][19]. Global warming has led to an increase in the elevation of the 0 °C isotherm [20] inducing glacier retreat, which enhances the debutressing of hillslopes and mass wasting processes that might result in a double phenomenon. First, larger amounts of erodible material are generated and available for remobilization as flows; second, new glacial lakes are formed, constituting unstable water reservoirs in high mountain areas [21]. At the same time, growing pressure from human activities is expanding toward unstable areas and there is a continuing high rate of deforestation due to logging, burning, and development, which are factors that increase severity of debris flows activity. Assessment of debris flow hazard is fundamental to develop preventive measures that reduce the severity of this potential natural hazard.
This review gathers research conducted on debris flows and hyperconcentrated flows in the Central Andes (30°-33° S), which is a region underrepresented in the international literature. Local reports and articles written in Spanish were included in this work to bring these cases to a broader audience. The main scope of our paper is to update the knowledge of these phenomena in this semiarid region trying to understand their behavior in the context of documented climate change [20]. Furthermore, this review highlights research gaps and a certain lack of local measurement data that are limiting a better understanding of debris flows in the Central Andes.

Semiarid Central Andes
The Central or Subtropical Andes (CA) extend from 31° to 35° south latitude comprising Western Chilean and Eastern Argentinean hillslopes [22]. The extent of CA has been defined differently [23], but this range system is considered as a transition zone between the hyperarid and semiarid conditions of the northern arid Andes and the southern humid Patagonian Andes. This climate transition segment is forced by the Arid Diagonal that separates summer precipitation patterns driven from the tropics, from winter precipitation driven by the westerlies to the south [24] (Figure 1). The summer precipitation cycle east of the Andes is originated by convective storms influenced by continental air masses [25,26], whereas winter precipitation on the highest elevation of the Andes is a consequence of the interaction between the moisture fluxes coming from the synoptic systems of the middle latitudes of the South Pacific, the topography of the Andes, the occurrence of cut-off low during the cold season [27], and the passage of cold fronts [28]. Topography and latitude strongly force temperature and precipitation distribution. Mountain ranges elevation decreases toward the south, favoring westerlies incoming the eastern side of the Andes (Argentina). The precipitation increases to the south due to the rising moisture flux from the Pacific Ocean and the weakening of the Southeast Pacific Subtropical Anticyclone. Between 30°-32° S and 71°-68° W, the average precipitation reaches 100-400 mm per year [30], including solid precipitation in the highest mountain areas. The snow coverage shows seasonal and inter-annual variations. Maximum solid precipitation occurs during the winter season, reaching 58% of the annual average snow coverage. Snow coverage reaches 12% during dry inter-annual periods against the highest incidence of snowfall periods with an average cover of 32% [31]. Anomalous winter precipitation has been associated with the El Niño-Southern Oscillation (ENSO) phenomenon [32]. Above mean precipitations are recorded on the Andes range during the warm phase. Opposite conditions prevail during the cold phase [24]. Likewise, the positive phase of the Pacific Decadal Oscillation (PDO) and the negative phase of the Antarctic Oscillation (AAO) match with increased winter precipitations [27]. Even though the ENSO impact on summer precipitation is rather weak [33], the precipitation variability is associated with equatorial symmetric circulation anomalies linked to ENSO-like warmer conditions after 1976/77 [34].
At present, this region is suffering a large drought receiving 30% less rainfall than normal over the past decade due to a negative phase of the Pacific Decadal Oscillation (PDO) [15]. In 2019, central Chile had an 80% rainfall deficit [16] (Figure 2). Climate models project a decrease in snowfall [35] and an increase of temperature [36] in the mountain areas of Central Andes. Meanwhile, the Mega-drought [16] has strongly forced the fire regime in the region in the last decade [37]. . Brown indicates vegetation that is less abundant and healthy than normal for this time. Gray highlights snowy and icy areas, or those with minimal vegetation. Green corresponds to abundant vegetation areas in valleys near glaciers and snow cover. The map was done by Lauren Dauphin, using MODIS data from NASA EOSDIS/LANCE and GIBS/Worldview and topographic data from the Shuttle Radar Topography Mission (SRTM).

Geographical Context
At these latitudes, the CA comprises different morpho-tectonic units from west to east, the Coastal Cordillera, the Main Cordillera, the Frontal Cordillera, and the Precordillera with very diverse lithologies. The Pampean flat slab segment, between 28° and 32° S, is related to an intense interplate and cortical seismicity, which has been related to landslides as well [39][40][41]. In this active tectonic segment of the CA, the highest peaks lay near 7 km (Aconcagua Mt 6960.8 m and Mercedario Mt 6780 m). This abrupt topography environment with low temperatures favored the presence of glacial and periglacial landforms. Permafrost predominates above 3200 m asl, which is indicated by the presence of debrisrock glaciers [42].
The main valleys of the CA run at lower altitudes fed by other streams and creeks with intermittent flows that generate alluvial fans in their mouth. In this type of environment, mountain villages have been developed, which occasionally are severely affected by channeled debris flows or floods. Major routes connecting both sides of the Andes are located along the main valleys, which often cut alluvial fans and creeks subject to debris and hyperconcentrated flows. These routes constitute the main trade corridors of South America, linking the Pacific and Atlantic Oceans through terrestrial transport. Furthermore, the Andean piedmont, where all these fluvial systems discharge, coexist with the settlement of main capital cities: La Serena-Coquimbo (30° S), Santiago (33° S), San Juan (31° S), and Mendoza (32°30′ S), which are historically affected by floods that are not considered in this review [43][44][45][46].

Debris and Hyperconcentrated Flows
The active tectonics and near-freezing temperatures of the high Andean mountains exacerbate sediment production feeding extensive valleys with steep slopes and long path gullies. This environmental context is associated with an arid-semiarid climate, where soil development and plant cover are low, which promotes erosion. The most efficient and common mechanism to remove these sediments is through aqueous fluidization. Andean catchments accumulate sediments commonly removed by sudden intense precipitations originating these debris/hyperconcentrated flows [39] mainly known as "aluviones" or "huaycos" in Latin America (Spanish or native terms).
Debris flows (DFs) are very rapid to extremely rapid (3 to 5 m/sec) gravity-driven flows composed of a mixture of debris and water [47,48]. They are recurrent phenomena within their path, since they usually follow pre-existing channels [49], and each event may include one single surge or several surges characterized by a coarse-grained head followed by a more liquid and fine-grained tail [48,50]. In CA, DFs are normally initiated on the basin headwaters of intermittent streamflows or dry creek by water saturation of hillslope material that moves downslope until the end of the alluvial system and deposit on fan landform. Whereas, hyperconcentrated flows (HFs) are defined as a two-phase turbulent flow composed of sediment-water mixtures that flow and behave as an intermediate flux between a normal streamflow or flood (mainly water with low fine sediments load) and a DF (higher proportion of coarse sediments than water) [50]. Even though HF was first defined according to a sediment concentration threshold of 20 to 60% [51] or 40-80% [52], more recent literature argue that these boundary values are arbitrary and not enough to define HF [50]. According to this author [50], HFs are unequivocal because of their high fluid viscosity due to suspended fine sediments that allow an intermittent dynamic suspension of coarse sediments (sand and fine gravel). Moreover, mean flow velocity and sediment loads transported by HF are greater than for water flow at analogous conditions of depth and slope.
When DFs lose coarse sediments by dilution and selective deposition, they may derive into HF. HF can also be generated when water flow entrains suspended sediments through valley erosion and entrainment. HFs are more related to streamflows saturated by sediments from upstream landslides or DFs or from natural or man-made dam collapse.
The failure of natural or artificial dams by overtopping, piping, or by mixing reservoir water with dam material may also generate outburst floods from naturally dammed lakes. In the CA, several glacial lake outburst floods (GLOFs) or landslide-dammed lakes outburst floods (LLOFs) have been reported [53][54][55][56]. A severe LLOF along the Santa Cruz River (32° S) in 2005 damaged a dam under construction and caused problems on water supply in San Juan city [57][58][59]. These low-frequency phenomena provide massive volumes of water flowing at high-velocities that can evolve as DF or HF due to erosion and sediment entrainment processes [50,60]. Since lakes producing GLOFs and LLOFs are most of the time located in remote areas, the resulting outburst flood can be misjudged downstream as a DF or HF, given that they were not originated by intense rainfall or the thaw of snow or frozen soil.
Finally, volcanic DFs, known as lahars, have not been reported in the study region, where Quaternary volcanic activity practically does not exist, starting south of 33° S. Lahars are generated by the mobilization of tephra deposits laying on hillslopes due to intense rainfall, snow, or ice melting or geothermal heating. Nonetheless, the role in DFs occurrence from the mobilization of laying pyroclastic material ejected from large historical eruptions further south (33°-36° S-e.g., Quizapu volcano) has not been analyzed in the study area yet.

Debris Flows, Hyperconcentrated Flows, and Debris Flood Records
During the 19th century, landslides were rarely reported by foreign travelers crossing the Andes along the Mendoza Valley (32° S) [61][62][63][64][65]. Santiago Molina, a muleteer, was killed by a falling block at Cortaderas locality in 1790 according to the communication of Miers [61] during his voyage in 1819. Afterward, landslide studies began to gain importance after the construction of mountain roads construction during the late 19th century [66] and with the Transandino railway operation along the Aconcagua (Chile) and Mendoza (Argentina) valleys [67]. Fatalities caused by falls, extreme cold, and snow avalanches along these corridors are reported by Osculati in 1834 (In: [68]), and serious damages on the Transandino railway were reported in 1890 during its inauguration (Los Andes, 1890 in [69]). Facilities and villages growing along the Andean valleys escalated landslides and snow avalanches records in the Argentinian side. Verdaguer [70] mentioned damage caused by rain-triggered DFs in the Precordillera during the summer of 1824 along Mendoza Valley (32° S).
On the other hand, historical floods and flows in the Elqui River (30° S) have been recorded since 1827 [71,72] (Figure 3). At least seven flood events were recorded from 1833 to 1915 by historian José Varela [73] and ≈71% of total floods (373 floods) recorded in between 1981 and 1990 correspond to the Elqui valley [45]. Particularly, several damaging debris flows were reported in the Elqui Valley after a hard colder winter in 1880. Violent DFs in Leiva and Uchumí creeks destroyed the railway along the Elqui valley in 1988, isolating the Vicuña village. The Algarrobal DF was an extraordinary event that occurred on May 21, 1934, which was triggered after 120 mm of rainfall in 72 h [73]. A bridge from Algarrobal creek (Rivadavia village) was dragged all the way to Diaguitas village more than 10 km downstream, where it destroyed the Andallito bridge. The Vicuña bridge and the railway were also completely damaged, leaving these mountain communities isolated. During the same rainstorm, a debris flow channelized in the Uchumí creek and dammed the Elqui River, generating a 10 m high lagoon that was fortunately quickly emptied 2 h later [73]. The most violent events were reported in Rivadavia village, where DFs coming from the San Juan and Cementerio creeks generated serious loses downstream. The material also blocked the Elqui River in the Cerda bridge. The dammed lagoon reached 500 m long and 10 m high, collapsing drastically with catastrophic consequences downstream [73]. The Diaguitas DF (30° S) struck the Vicuña village on 22 April 2004, during the night, after intense rain (90 mm in 12 h). Several tons of muddy material with blocks of different size, trunks, and branches came down from the Puyalles creek. This catastrophic event affected 150 families, injured 60 people, destroyed some houses, eroded farmland, razed grapevines, and devastated the irrigation system. The village´s school was covered by a muddy layer of 0.70 m [74]. Losses were estimated between six and eight million dollars [75].
The most catastrophic DF of central Chile occurred in November 1987 on the Colorado River catchment (33° S). The Parraguirre rock slide-avalanche [76] generated a DF that destroyed the fields, roads, bridges, and damaged the under-construction Alfalfal hydroelectric power plant along the Colorado River. This event caused 41 fatalities [1,2,77]. The DF reached an estimated volume ranging from 2.5 to 15 × 10 6 cubic meters [2]. A summary of the main historic DF and HF events in central Chile is presented in Table 1.  [11] Estimated volume greater than 10,000 m 3 , splash marks up to 3-4 m high On the Argentinean side of the Andes, an extensive record of more than 300 historical DFs exists for the Mendoza Valley (32°30′ S) [7,17,19]. The most dangerous reported DFs were channelized in the Seca and Camino creeks in 1976, 1980, 1999, and 2000. One of these events dragged a car in 1968. Another DF from Avispas Creek dragged out a lorry that fell into the Mendoza River with one fatality. Violent DFs have been reported as well in the Polcura (known as 60 creek) and Guido locality in 1999, 2001, and 2005, swamping the international road [17,84]. Based on these records, a recurrence interval of 4 years was estimated for the Frontal Cordillera, while this interval is higher (15 years) for the Precordillera [7].
Along the Blanco Valley (33° S), many damaging DFs were recorded in the main creeks as Salto, Angostura, Mulas, and Colorada, impacting adversely the local population [10]. Damaging torrents were reported in January 1946, destroying a bridge in Chacritas, and in January 1947, affecting a school at Potrerillos village. On December 1954 and 1955, an intense rainstorm generated several debris falls along the Colorada creek, disallowing the access to the Vallecitos ski tourist complex. A DF cut again the access on 7th February 1967. Finally, the 7th February 2013 another DF channeled in Las Mulas creek destroyed a small bridge and isolated people in the Valle del Sol locality. This event severely affected the drinking water supply of these communities for at least four days [85].
After an intense rainfall affecting the central-west part of Argentina on February 8th, 2013, a DF affected the Agua Negra international road (30°20′ S). This road connects the Agua Negra valley, San Juan (Argentina) with Elqui valley (Road CH40), and Coquimbo-La Serena (Chile). Unfortunately, this event led to the death of a 10-year-old boy when the vehicle in which he was traveling lost control trying to pass the muddy moving mass material. The material volume was estimated at 50,000 m 3 with an average speed of 35 km/h. As a result, the international road was closed for several days and isolated 11 families at the Chilean customs [86].
Impacts of DF and HF generated by intensive summer rainstorms have been also documented for the lowlands. Polanski [87] studied these processes proposing a classification for DF and floods based on their deposit characteristics in this arid region.

Debris Flows Characterization
Both DF and HF constitute hazardous processes due to their great velocity, long runout distance (several km), and high impact loads capable of moving large boulders and trees, implying significant destructive power [88]. The estimated volume of DF in the CA fluctuates from 250 m 3 [11] to 15 × 10 6 m 3 [76]; these values are possibly underestimated, as the real water volume is rarely known. Figure 4 shows that greater recorded volumes correspond to GLOF and LLOF. The transported material can involve boulders up to 4 m in diameter. The DF recorded in Juncal valley (32° S, Valparaíso region) on 18th August 1987 removed a 450 t stone. Boulder and block material tends to be thicker in Chile hillslopes than in Argentina side with more gentle slope basins [11]. Nonetheless, many mud flows have been described in Elqui valley (30° S) (Table 1). Less common are woodladen flow (WLF) in this semiarid region, even though a WLF event was described in February 2017 in the Pocuro stream (Central Chile, 33° S) [89]. Different flood or flow pulses are often generated because the sediments form dams that contain the flow for a variable time, and when these give in, a new flood pulse is generated [76,80]. Flood heights of 0.7 to 35 m depositing 1 to 5 m of sediments [7,10,11] have been reported.
The larger sizes of Andean catchments influence a large debris supply, while steep slopes exacerbate the flows velocity. The DF/HF affecting this region are relatively highspeed events with velocities exceeding 30 km/h [76,90]. In the Mendoza Valley (32°50' S), some DFs reach 126-180 km/h in very steep ravines [41]. Historical DFs caused by rainstorm during summer seasons in Anchoris (33° S) reach 1-2 m·s −1 in velocity, remobilizing a great amount of fine material [87]. The velocity of 4.7 m·s −1 for the first phase of the event was estimated from videos recording in Pocuro [89]. Curiously, higher speeds were indirectly established for the Amarilla creek and Paraguirre DF, whereas lower velocities were estimated for historical GLOF (Figure 4).
Channeled alluvial DFs mainly affect smaller catchments traveling short to moderate distances (up to ≈15 km), but when they end into a main river streamflow, they may travel larger distances in order of 100 km or more. The travel distance of Andean DF/HF is very variable, which was probably linked to the catchment areas, slopes, and nature of mobilized material (Figure 4). GLOF and LLOF events have longer runouts (up to ≈250 km). Still, outburst wave dissipation is strongly forced by the dammed lake volume. The Plomo GLOF event on 1934 reached 170 km, developing a frontal wave of variable height (7 to 12 m) according to the narrowing/widening of the valley. This event destroyed several kilometers of road and railways, seven bridges, a hydroelectrical power plant, and a hotel, affecting 750 km 2 of cultivated land and killed 18 people [53]. Particularly, the HF generated by the Los Erizos LLOF event recorded the longer runout with a travel distance of 250 km (Figure 4).
Flow lubrication by ice input can propel excessive travel distances as well. This was the case of the 1987 Parraguirre-Alfalfal DF in the Colorado River (33° S) that traveled 50 km downstream arriving close to Santiago, the capital city of Chile [11]. Out of the study area, the Santa Lucia DF (43°22′ S, 72°22′ W) on 16th December 2017 was generated by a 122 mm rainfall that triggered a rock slide that in turn caused partial melting of a receding glacier. The ice-debris mass traveled 10 km with an average speed of 72 km/h, mobilizing 7.2 × 10 6 m³. Two million cubic meters of sediments were deposited on the alluvial cone where the Santa Lucía village is settled, affecting 536 ha and resulting in 21 fatalities with a missing person [74].

Spatial Distribution
The reconstruction of debris flows occurrence, similar to another natural process in mountain areas, is strongly constrained by historical sources. The availability of chronological documents is rare in remote regions, and event reports are concentrated exclusively on the inhabited sectors along the Andean valleys ( Figure 3). For this reason, inventories from remote sensors are becoming more significant in the CA and have partly overcome this limitation. From these inventories, the occurrence of DFs in high Andean altitudes dominated by a glacier environment has been ruled out [39,91,92] (Figure 5). However, those slopes located in lower positions related to periglacial domains where permafrost conditions could be affected by degradation take a key role [93][94][95]. ). Light blue polygons represent an isolated DF and dark blue ones correspond to debris flows area where small DFs has been gathered. A limited spatial distribution of DFs is clearly observed in the highest mountain zones due to glacial domain. Along the main valley, several DF geoforms match with historical records (see Figure 3).
These findings match as well with reported DFs where 75% of these events occurred in the Frontal Cordillera with altitudes between 1880 and 2900 m asl along the Mendoza Valley (32° S). Only 16% of the complex events happened in the higher Main Cordillera (heights between 3000 and almost 7000). The spatial distribution of historical DFs is constrained as well by the precipitations conditions. Snowfall over the higher catchments promotes snow avalanches during winter season, while DFs occur in the springtime associated with snow or ice thawing [18,19,39] and rain induced-DFs occur in lower elevations areas of catchments during the summer season [17][18][19].

Precipitation
In this arid-semiarid region, DF/HF occurrence has been mostly associated with intense summer rainfall [11,17,79,86,97,98]. Most of them are convective storms associated with multiple synoptic conditions, such as the negative anomaly of the zonal component of the wind in the middle layers or, in certain cases, the presence of a cut-off low promotes the formation of convective storms on the eastern slopes of the Andes mountain range [99]. Meanwhile, in the western side of the Andes, DFs tend to occur not only linked to summer rainstorms but also to intense and persistent precipitation during winter.
Total day precipitation threshold values are relatively low compared to humid-climate environments. In the wetter Chilean side, daily precipitation of more than 60 mm/24 h may trigger DFs in the CA (33° S) [100]. However, DFs have been generated with 12 mm/h in the Macul Creek (Santiago, 33° S) [80,101]. In the driest eastern Argentinean side, a range between 6.5 and 12.9 mm/24 h was determined for the historical DFs in the Mendoza Valley (32°30′ S) [17]. A very low threshold precipitation of 5.5 mm/24 h was estimated for the 8F event in 2013 that occurred in the Amarilla creek (30° S, San Juan) [86] ( Figure 6a). This threshold value was estimated using satellite rainfall data (CMORPH), as weather stations in the surrounding area are lacking. Nonetheless, the influence of antecedent precipitation is not dismissed, as accumulated precipitation during previous days was of 13.2 mm in this case. The role of precipitation during the previous days was denoted as well for the events studied in the Mendoza River valley (32° S), as 60% of the cases recorded precipitation until 3 or 5 previous days [7,17]. In Chile, the accumulated rainfall influence during 20-25 days before the event was found significant in the Andes mountain front at the latitude of Santiago (33° S) [93]. Hence, not only the daily rainfall, but the accumulated rainfall and the altitude of the zero isotherm are crucial when analyzing the causes of DF in this region [46].  [70]; (c,d) DFs in Picheuta area affecting the international road on 7th Feb 2013 [11]; (e-g) events recorded in Seca creek on 23th January, one roofed a bulldozer during the night [11]; and (h) event occurred in Picheuta region on 7th February 2015 with an estimated volume ~ 115.000 m 3 [11].
In the middle portion of the Elqui Valley (30° S), where the probability of permafrost is limited to the highest summits, 1 h of precipitation is enough to trigger DFs and HFs [79]. The accumulated precipitation necessary for the triggering of DFs varies between 4.5 and 10 mm depending on zero-isotherm elevation, since for this area, the precipitation occurs fundamentally in winter (see Section 2). The triggering of HFs requires about two millimeters less for equal variations of the zero-isotherm elevation, but also for these events, their occurrence was not recorded for values less than 4.5 mm/24 h. Considering that HFs are generally lower energy phenomena (lower peak discharge, speed and size of suspended sediment than DF), it is logical that their liquid precipitation thresholds tend to be lower as well [50,101]. Unlike the southern CA, here, it was calculated that soil moisture preceding the rainfall burst does not condition the potential subsequent DFs occurrence, which could be due to different failure mechanisms and/or the more arid conditions of this sector that generate greater evapotranspiration [79].
The perspective of future temperature in the CA reveals a general rise [14]. A still series of summer precipitation of the Mendoza river valley shows a negative trend in the 1905-2018 period; this tendency varies. In fact, temperature trends will depend on the analyzed time period. Since 1965, precipitation has increased, matching with a sustained rise of DFs occurrence (Figure 7a). Although a greater number of DF have been recorded in the last decades, there is still more information to be collected to attribute a relationship between such events with climate change.

Temperatures and Zero-Isotherm
Temperature is a key parameter for DFs generation in the mountain environment, as it influences the cryosphere melting [18,76,78]. Analysis of temperature series  of the Mendoza river valley shows a rising trend. This increase in temperatures could be correlated with more activity of DFs in this valley (Figure 7b).
Even though DFs from the melt of seasonal snow and ice within the active layer has lower frequency with respect to other triggering causal factors, DFs related to this melting processes have been recorded in some places (e.g., [12,18,103]). The Aconcagua Park (31° S, Argentina) visited annually by more than 120,000 international tourists has been affected by violent DFs at the beginning of the warm season [12,92,104]. A DF coming from Durazno creek was reported in 1985, and another event was recorded in the Vacas valley in 2008 [105]. Moreover, two climbers reported an ice-core moraine expelling saturated fine material in the Vieja Alta Creek near the Plaza Guanacos camp in February 2004. This event partially blocked the Vacas River (Figure 8c) [105]. An amateur video recorded the DF on January 2016 in the Blanco alluvial fan after several days of successive increases in air temperature (video: https://www.youtube.com/watch?v=jNVwspunNn0). This DF was generated from a creek with heavy snow accumulation coming from the Almacenes Sur Mount (5410 m). At least 3 pulses were recorded in this event (Figure 8a,b).  [94]; (b) the DF occurred in January 2016 when tourists were crossing the Blanco creek [93]; (c) DF generated from an ice-core moraine in the Vieja Alta Creek near the Plaza Guanacos camp in February 2004 [104]; (d) DFs generated from a debris-rock glacier melting in the Negro creek on 28-29 December 2015 with an estimated volume of 30,018 m 3 [12], the circle pointed out with a car; and (e) DFs caused by the melting of a debris-rock glacier in La Salada creek (Mendoza Valley 32°50′ S) with an estimated volume of 238,891 m 3 [12], circle corresponds to a car.
In December 2015, the successive climbing of the zero-isotherm produced the melting of a debris-rock glacier located in La Salada creek (Mendoza Valley 32° S) generating four successive events that blocked the international road (Figure 8e). Due to this associated substrate warming, also, a debris-rock glacier located at the headwater of the Negro creek melted, triggering another three DFs during the same period [12] (Figure 8d). Those DFs generated by the melting of glacier tongues have generally performed successive events during several days, unlike in DFs triggered by rainfall, which have a sudden discharge in a single event, even though several pulses may occur.
Although direct measurement on soil temperatures is critical, the temperature is established by correlation with air temperatures or by insolation exposition [106]. The study of the thermo-radiative conditions previous, synchronous, and subsequent to debris flows from the melt of seasonal snow and ice within the active layer is complex in this region due to the scarce meteorological monitoring stations in the headwaters of the hazardous watersheds. In the upper basin of the Elqui Valley (30° S), these events could be investigated in greater detail due to the good meteorological network of the area, with stations at hourly resolution and located at high elevations [78,79]. The debris flows have several thermo-radiative patterns such as (a) greater correlation with the daily accumulation of heat available for melting than with the daily maximum in this heating, (b) occurrence in clusters of 3-5 days during the early summer, (c) initiation of clusters at the warmest day until that moment of the summer trimester, and (d) an accelerated increase in temperature over the days previous to the beginning of the clusters. However, the amount of water generated by summer snowmelt or melting of the active permafrost layer is still uncertain in the CA [78].

Climate Phenomena-ENSO Role
The warm phase of the El Niño Southern Oscillation (ENSO) [32], linked to above average precipitation, promotes DF activity in the CA [7,[106][107][108]. The largest event occurred in the Elqui valley took place in 1997 during the most intensive El Niño event [45]. The peak of DFs activity has been documented during the El Niño warm phase in the Mendoza Valley (32° S) [7]. At least 55 events were exceptionally reported during 29 consecutive days in the 2015/2016 South Hemisphere summer when the GODZILLA event was hitting [12]. In the highest Andean sectors, westerlies intensification during El Niño set off a denser snow precipitation that promotes DFs escalation on the following warmer season [18].
On the western side, the 1993 event on the 2nd and 3rd of May that affected seriously the city of Santiago de Chile was generated by heavy rainfall associated with warm conditions in the troposphere linked to a moderate ENSO event [76]. During this single event, a maximum intensity rainfall of 12 mm/h was reported, which is a high value expected only every 25 years [109].

Environmental Impacts
Large economic losses resulting from unexpected violent DFs have been underestimated in the CA, with the environmental consequences even less considered. In turn, these effects might be important for catchment dynamics.

Alteration on Stream Dynamics
The DFs are generated in valley hillslopes and move downstream, removing material along its path. This clastic and organic material is discharged along alluvial channels, creeks, or streamflows controlling sediment budgets in the catchment [110,111]. The sediments mobilization from hillslopes switches further slope erosion, while sediments transferred downstream supply removable material for future flows. This debris material transfer certainly drives basin dynamics.
The DFs provide important sediment loads, and bank undercutting increases its erosion power. This sediment input derivates on changes in geomorphological aspect of valleys. As an example, overflows, as a consequence of several DFs channeled in the Blanco Valley (33° S), were recorded in 1942, 1945, and 1957. These powerful erosive events caused changes on valley configuration, and bridges and roads were swept away.

Effects of Landslide Dams on Streams and Valleys
The long-term contribution of sediment material forces the progression of alluvial fans toward the main drained streamflow, bringing local channel constriction and riverbed adjustments. Whereas, fast DFs have generated partial or complete river blockages that breed impounded lakes. The Elqui Valley (30° S) was dammed in two sectors during the Algarrobal DF in 1934 [73]. The Turbio River, upstream in the Elqui Valley, was partly blocked during 2015 [79] (and 2017 events (Table 1)).
The Mendoza River (32° S) was partially obstructed by the Soltera DF on 7th February 2016. Historical DFs have also obstructed partially the Vacas Valley such as during the events reported in 2004 and 2008 [104].

Impact on Ecosystems
Stream flows with great sediment load may also generate negative effects on the coastal ecosystem. The sediment discharge into the sea should generate a great mortality over banks of filtering organism, such as clams, generating environmental and economic problems. In the Coquimbo region, this issue was really important before the operation of the Puclaro reservoir that was built in 1997.
In addition, streamflow with great sediment load generates troubles in the water treatment process for consumption. Commonly, these sediment-laden flows produce damage to the filters occasioning reduction on drinking water in the Andean communities.

Discussion Concerning Climate Change and Final Remarks
Until now, precipitation threshold values for triggering DF have been mainly established from total day precipitation, as records of precipitation intensity are lacking in the most catchments of the study area. The occurrence of DFs in the CA is mainly related to rainfall during 1-2 days, thus understanding that the synoptic situation associated with these events may improve prediction and system alerts. A small positive trend in rainfalldriven landslides during summer in the eastern Argentinean piedmont persists since the mid-20th century (32° S), matching with an increase on Atlantic monsoon precipitation.
Less common events are those related to the release of melt water from the periglacial environment located in the upper sectors of Andean catchments. The warming and ensuing rise of the zero-isotherm in a context of global change will affect these environments, increasing the frequency of DFs. In this sense, the behavior of melt-driven events could vary for the coming decades depending on the type and the water source. The future situation of the DFs triggered from the melt of seasonal snow and ice within the active layer is not clear; on the one hand, the solid precipitation that feeds the seasonal snowpack and active layer will decrease, but on other hand, an increase is expected in the maximum potential melting rate due to the rise in maximum temperatures during the ablation period [78].
Higher temperatures are raising the zero-isotherm elevation; that stimulates permafrost degradation with increased suitable water-saturated areas during rainfall. The intensification of the warm phase of El Niño-ENSO linked to climate change could force the distribution and frequency of DFs as well as their intensity and volume whether ice glacial input is involved.
Overall, the connection between climate change and debris flow occurrence is not fully proved. It is feasible that concentrated episodes of heavy rain may trigger more frequent DFs, even though a decreasing winter precipitation is predict by global models. In the same way, a weakened permafrost and glacier cover might produce more DFs. Longterm data series are needed to establish a causal link with climate, which need to be properly analyzed with statistical methods. As far as we are aware, this kind of analysis has not been carried out yet. Records of DFs in the Mendoza river basin show a rising trend in recent decades, but the cause of this pattern is still uncertain. Therefore, this issue is worthy of further research and would be a natural continuation of this line of research in semiarid zones.
Unanswered research inquiries remain about the behavior of the DFs. The role of seismic shaking on debris supply in the DF-generating slopes has been little explored in this active tectonic region. Studies focused on the role of cineritic levels (tephra) as facilitators of DF events are also needed. The consequences of ensuring fires due to the megadrought on slope behavior are still not known. For example, landslides were induced after the wildfires of California in 2020 [110], and this still has to be studied in the area.
Finally, the most severe and impacting DFs in the CA are those channelized in alluvial creeks that may feed into river streamflow; communities located on alluvial fans will be more vulnerable. Specific geomorphological and hydraulic studies taking care of this alluvial environment dynamics are particularly required to reduce potential damages. Hence, the urgent implementation of mainly structural preventive measures is required based on detailed technical studies. Even though vulnerability assessment is a useful tool on a certain scale, it does not reduce the exposure of Andean communities to risk per se. Institutional Review Board Statement: Ethical review and approval were waived for this study because review article do not require follow this type of procedure.
Informed Consent Statement: Not applicable. This study did involve humans.
Data Availability Statement: Data reported in this article are available in published papers and PhD thesis of some authors of this review.