Gravitational natural hazards such as snow avalanches, rockfalls,
shallow landslides and volcanic activity represent a risk to mountain
communities around the world. In particular, where documentary records about
these processes are rare, decisions on risk management and land-use planning
have to be based on a variety of other sources including vegetation,
tree-ring data and natural hazard process models. We used a combination of
these methods in order to evaluate dynamics of natural hazards with a focus
on snow avalanches at Valle Las Trancas, in the Biobío region in
Chile. Along this valley, natural hazards threaten not only the local human
population, but also the numerous tourists attracted by outdoor recreational
activities. Given the regional scarcity of documentary records, tree-ring
methods were applied in order to reconstruct the local history of snow
avalanches and debris flow events, which are the most important
weather-related processes at respective tracks. A recent version of the model
Rapid Mass MovementS (RAMMS), which includes influences of forest structure,
was used to calculate different avalanche parameters such as runout distances
and maximum pressures, taking into consideration the presence or absence of
forest along the tracks as well as different modeled return periods. Our
results show that local
In mountainous areas, settlements and infrastructure are commonly endangered by natural hazards such as snow avalanches, rockfalls, shallow landslides or volcanic activities. The documentation of past events of these natural hazards, including a detailed description of their properties (e.g., magnitude, affected area, type of damage) as well as their effects on and interactions with ecosystems, are important for assessing risk (Stoffel and Huggel, 2012; Gądek et al., 2017). In mountain ranges with a long history of settlement such as the European Alps, documentary records of gravitational natural hazards have been compiled for more than five centuries due to the frequent contact between mountain inhabitants and these phenomena (Laely, 1984). Yet, only more contemporary records have provided suitable information for accurate risk assessments (Wisner, 2004). In many other mountainous regions, human settlement has been rather limited until the 20th century and has substantially increased only during the last decades (Riebsame et al., 1996; Romero and Ordenes, 2004; Su et al., 2011). It is thus particularly important in more remote mountain areas, in which human populations and infrastructure are expanding, to learn more about (i) different natural hazard processes, (ii) how mountain ecosystems can be used to reduce risk and (iii) how people can be affected by these processes.
In regions where documentary records of past events of gravitational hazards are limited, dendrochronological methods provide a valuable and accurate tool for their spatiotemporal reconstruction (Braam et al., 1987; Shroder, 1980; Pop et al., 2016). The maximum length of this type of reconstruction depends primarily on the maximum age of the trees available for sampling at the study area. These ages, in turn, depend on two main factors: (i) the site history related to human interventions such as induced forest fires, logging and grazing (Mundo et al., 2013) and (ii) the frequency of events including avalanches and mass movements (i.e., more frequent events typically result in lower tree ages; Corominas and Moya, 2010). Documented damage produced by gravitational natural hazards (in particular snow avalanches) to old buildings, such as churches, has made it possible to couple long-term dendrochronological reconstructions with very precise information on partial and total destruction of infrastructure, thus validating the use of tree-ring methods for such studies (Fliri, 1998).
Forests can protect infrastructure from natural hazards in mountainous areas
(Brang et al., 2006; Getzner et al., 2017). In the case of snow avalanches,
which are main the focus of our investigation, mountain forests are known
for stabilizing the snowpack in potential avalanche release areas and thus
help to prevent avalanche initiation (Weir, 2002; Bebi et al., 2009,).
Regarding small to medium-scale snow avalanches (i.e., < ca. 10 000 m
Study area. We focused our studies on natural hazards and protection forest along the road (red line) connecting Las Trancas with Nevados de Chillán. The north-facing slope adjacent to the road is much less vegetated compared to the south-facing slope, which is mostly explained by the lower radiation received on the latter and thus longer remaining snow cover that provides water in drier periods of the year (see also Fig. 2).
Eco-DRR has been defined as “the sustainable management, conservation and restoration of ecosystems to reduce disaster risk, with the aim of achieving sustainable and resilience development” (Estrella and Saalismaa, 2013). Some of the most commonly applied measures to disaster risk reduction in Latin American countries are nature-based solutions such as protection forest, bio-engineering and eco-construction measures (McBreen, 2016). Nature-based solutions are particularly attractive due to their low costs, the use of local materials and their minimal impact on the landscape. A further advantage is that they have the potential to simultaneously reduce natural hazard risks and to provide benefits such as clean drinking water, increasing biodiversity, recreation, wood production and a wide range of other ecosystem services (Baró et al., 2015). In the Southern Andes, human settlement in villages and cities as we know them nowadays started only about 150 years ago, coinciding with one of the largest European immigration waves and the expansion of populations towards the slopes of the Andes (Germani, 1966). This process resulted in a sudden exposure of human populations and infrastructure to relatively well-known natural hazards but in a new environment. Since then and up to the present day, knowledge transfer has been made from (particularly) the European Alps to the Southern Andes to better understand the triggers and dynamics of natural hazards (Casteller et al., 2008; Stoffel et al., 2012). Yet, further major efforts are still needed to fulfill these goals.
Native forests in Chile provide several ecosystem services that directly or
indirectly benefit human wellbeing, including the provision of food (e.g.,
The overreaching goal of this investigation is to evaluate the role of mountain forest ecosystems as an Eco-DRR measure against natural hazards (primarily snow avalanches) at Valle Las Trancas, Nevados de Chillán, in the Biobío region of Chile. Within this region, avalanche-related casualties have been reported from the ski resort Nevados de Chillán and at the road connecting the village of Las Trancas with infrastructure located along the valley ending at the ski resort (La Discusión, 2011). Yet, the risk of snow avalanches in this area is largely unknown and, besides road temporal closures, no other passive or active measures (except within the limits of the ski resort Nevados de Chillán) are in place to reduce avalanche-related risks.
Besides snow avalanches, other natural hazard processes such as volcanic activity, rockfall, debris flows and shallow landslides also affect our study area. We thus complemented our detailed study on snow avalanche and vegetation interactions with investigations on other natural hazard processes that occur in the area. This allowed us to evaluate the importance of native forests in relation to different natural hazards in the study area and to provide a suitable basis for proposing adequate Eco-DRR measures.
This study is part of the EPIC (Ecosystems Protecting Infrastructure and Communities) project, which has been implemented by the International Union for Conservation of Nature and several local partners and aims to demonstrate how improved ecosystem management can reduce the risk of disasters and help to adapt to climate change with different study cases around the world. The study is focused on answering the following questions. (i) What is the effect of avalanches and debris flows on the structure of affected forests? (ii) What is the influence of forest on avalanche runout and avalanche hazard? (iii) What assessment approach is most suitable to assess avalanche–forest interactions where information on past events is limited? (iv) What are the potential Eco-DRR measures for reducing overall risks in the study area?
The study area is located in the Biobío region of Chile in northern
Patagonia. We focused our investigations at Valle Las Trancas, an
east–west facing valley that is located within the limits of the transition
zone of the biosphere reserve Corredor Biológico Nevados de
Chillán – Laguna del Laja (Fig. 1). This biosphere reserve was opened by UNESCO on 29 June 2011. The presence of the northernmost
population of the Patagonian huemul (
In recent decades Valle Las Trancas has experienced high economic development due to its popular tourist attractions, mainly related to mountain sports and thermal springs (Pfanzelt et al., 2008). The municipality of Pinto, in which the Valle Las Trancas is located, lacks official land-use plans and regulations, leading to uncontrolled establishment of new buildings in the area. This development has also increased the demand of water for local consumption and of wood for heating.
A growing demand along Valle Las Trancas of new tourist infrastructure, in an area lacking appropriate hazard maps, can result in increased exposure to gravitational processes with serious risk for people. In addition, the ski resort Nevados de Chillán, located in the basin headwater, is one important local actor of the private sector authorized to cut patches of native forests to clear the land for new ski tracks. This potentially results in tradeoffs with the forests' provision of protection against snow avalanches.
The investigated area is located between the temperate and Mediterranean
macrobioclimates (Luebert and Pliscoff, 2006). It is characterized by dry
summers and cold winters; more than 75 % of the annual precipitation fall
in winter between April and September (Donoso, 1993). Local climate data
obtained from a weather station located at Las Trancas
(36
A high diversity of geological substrates caused by high volcanic activity
characterizes the study region. Adjacent to the study area is the
Nevados de Chillán volcano, which consists of the Cerro Blanco and Las
Termas subcomplex (formed about 40 000 years ago). The most recent eruptions
from the Nevados de Chillán volcano occurred in the years 1973, 2003 and
2016; due to the presence of tourist facilities, it is categorized as one of
the volcanos with the highest damage potential in Chile (Dixon et al.,
1999). These diverse geological substrates are overlaid with characteristic
soils classified as Inceptisols, which develop from these volcanic
products (Freiberg, 1984; Soil Survey Staff, 2014). The high diversity of
geological substrates leads in turn to 11 different vegetation units
(Pflanzelt et al., 2008). Close to the riverbed in humid localities,
Identification of gravitational natural hazards within the study area was
conducted in an early stage of this investigation and complemented with subsequently obtained
geographical information. We first focused on
spatially demarcating the avalanche tracks that reached (according to
available records) infrastructure along the road of Valle Las Trancas.
Spatial analysis and interpretation of aerial images were made using ArcMap
10.2 (ESRI, 2013). Google Earth images were used to display
the information and acted as a basis for the interpretation of aerial images.
Additionally, available aerial images of the study area (from the years 1978,
1979, 1982, 1998 and 2000) were stereoscopically analyzed in order to
identify chronological modifications in the forest structure due to natural
disturbances following the approaches of Bebi et al. (2001) and Panayotov et al. (2011).
For the assessment of the current forest structure in selected
avalanche tracks, we complemented these analyses with high-resolution
remote-sensing data (including a digital terrain model (DTM) and cover
image) obtained in February 2015 with an eBee drone (resolution: 8 cm pixel
We complemented these data with repeated field surveys. Trees that were either sampled or simply identified for features of interest were individually positioned using a portable navigation system (model: GPSMAP® 64s). Forests in which the spatial extent of past avalanches was visible in the field were mapped using the same method, in particular with borders and avalanche runouts. Photographs were taken of avalanche release zones (which are mostly above the treeline and on steep slopes) from different angles and locations (Fig. 2). These photographs subsequently helped for spatial demarcation of these release zones using a Geographic Information System (GIS).
Snow avalanche tracks, debris flow tracks and forest plots (fp) at the study area. Tracks depicted in blue correspond to snow avalanches, whereas those in red correspond to both snow avalanches and debris flows. A total of 18 tracks were identified in the study area and 19 forest plots were defined.
Between the village of Las Trancas and the base of the ski resort Nevados de
Chillán (36
Two different approaches were applied to reconstruct a spatiotemporal
occurrence of gravitational natural hazards. The first approach aimed to
evaluate the forest age and structure in specific plots. The locations of
these forest plots were selected according to the following criteria: (i) plots
are placed at each natural hazard patch (i.e., avalanche/debris flow
track) with the presence of forest; and (ii) at each forest stand was in between
two natural hazard patches. Each plot had a radius of 10 m, thus covering a
surface of 314 m
The second approach focused on identifying particular trees at selected avalanche tracks with the largest amount of visible damage (which is typically associated with repeated events), such as broken stems and/or branches, decapitation, uprooting and exposed scars. From the total of 18 avalanche tracks identified, we applied this approach (hereafter called “selective sampling”) to two of them (defined as av6 and av9; Table 1). We selected these two tracks for the following reasons: av6 had one recorded avalanche event in 2011 (which reportedly caused one casualty) that could be used for calibration purposes, whereas av9 presented a very clear avalanche runout zone and borders with numerous trees evidencing recent events. Selected trees were sampled using increment borers (Haglöfs; width 0.5 cm) as well as handsaws, and in all cases unnecessary damage to living trees was avoided in this protection forest (Bebi et al., 2001; Dorren et al., 2004). Morphological features of sampled trees, such as the height of the damage along the main stem, size of the scars and inclination of the stems were registered and complemented by drawings and pictures.
Samples collected at different tracks and forest plots.
All tree-ring samples were prepared following standard dendrochronological techniques (Stokes and Smiley, 1968) and subsequently ring-width measured using a LINTAB system from RINNTECH, Germany. The software COFECHA was utilized to detect (and if necessary to correct) dating errors, using chronological tree samples from adjacent undisturbed sectors of the forest as a reference (Holmes, 1983). Further, every tree-ring sample was carefully examined for the presence of scars, growth releases and growth decreases, as well as peculiar features in the wood such as frozen rings or wood compartmentalization. Growth reactions were grouped into three different categories following Dubé et al. (2004). Both the expert (Stoffel and Corona, 2014) and indices (Shroder, 1978) approaches were then applied in a comparative way to achieve a final reconstruction of events at each avalanche/debris flow track.
In regions where the precipitation patterns determine the occurrence of debris flows during the vegetation period of trees, it is possible to differentiate between the occurrence of snow avalanches and debris flows by identifying the position of the injuries within the tree rings (i.e., dormancy for avalanches and active growth for debris flows; Stoffel et al., 2012). In our study area, where most precipitation occurs in winter, we could determine that most of the studied debris flows occurred in spring due to rain-on-snow events.
Numerical simulation models have been used to reconstruct snow avalanches in
various mountain regions including the Andes (Valero et al., 2016), central
Europe (Blahut et al., 2017), Japan (Takeuchi et al., 2011) and Scandinavia
(Gauer and Kristensen, 2016). The software RAMMS was used to conduct snow
avalanche simulations at Valle Las Trancas. RAMMS is a two-dimensional
software to calculate mass movement dynamics in a three-dimensional terrain
(Christen et al., 2010). We applied RAMMS::AVALANCHE version 1.6.25© (from WSL-SLF) with its new forest applications. Compared to standard
simulations with RAMMS::AVALANCHE, this version differentiates among the
effects of different forest structural characteristics on avalanche runout.
Therefore, spatial information on crown coverage, roughness and DBH had to
be defined from the forest plot survey and remote-sensing data for the
simulations at each track. Based on the forest type, crown coverage and
roughness in the area, we subsequently deduced explicit spatial information
on the detrainment coefficient
At avalanche tracks 6 and 9 (av6 and av9), specific information regarding runout distances from past events was available (both from tree-ring data and documentary records). Out of the results of our tree-ring analyses, trees with synchronous reactions to avalanches at these run-out zones were used to delineate the affected area of past (tree-ring dated) events. We then compared this information with the spatial data provided by the numerical simulations with RAMMS. To further evaluate the effect of local forests on avalanche dynamics, different simulations considering both the presence and absence of forest were conducted. Further, avalanche simulations were conducted at a selected track (av6) in order to quantify the improved protective capacity of the forest by conducting additional afforestation, with the assumption that the additionally considered forests have similar functional properties to the current forest.
Our analyses show that avalanches and debris flows play a role in
influencing the structure of the forest at the study area. Canopy densities
of the investigated plots range between 20 and 90 % (average value of
60 %; Fig. 4). By grouping plots in different categories (i.e., no
disturbance, debris flows and snow avalanches), we observe that avalanches
reduce the canopy density by about 20 % in relation to undisturbed
forest, whereas for debris flow this reduction is 10 % (Fig. 5); these
results are, however, not significant (
Canopy density distribution (%) per plot.
Canopy density distribution (%) per disturbance.
Reconstructed borders and runouts at avalanche track 6 using tree-ring data from the first registered event in 1942 up to the sampling time in 2015. The green circles indicate sampled trees.
Reconstructed avalanche and debris flow events at the investigation area. Fp3 and fp5 show event years in the debris flow tracks (df3 and df4), whereas fp8, av6, fp9, fp12, av9, fp18 and fp19 show event years at different avalanche tracks (av6, av6, av7, av8, av9, av17 and av18).
The results of our tree-ring analyses allowed us to reconstruct past snow avalanche and debris flows years at the investigated tracks. Scars and abrupt growth changes provided the most reliable indicators of past events. Table 2 contains a list of all years during which, according to our reconstruction, either snow avalanches or debris flows occurred at the investigated tracks. Only forest plots located in avalanche or debris flows patches are considered here (see Fig. 3 for reference), in addition to tracks av6 and av9 in which the selective sampling was conducted.
The first year during which an avalanche event was detected in our tree-ring data is 1942. The main factors that prevent the extension of our reconstruction beyond that year are, on the one hand, the young age of the local forests and, on the other hand, wood decay. At avalanche track 6 (av6; using selective sampling), a forest plot (fp8) was also selected in order to compare the event reconstruction using both approaches. We observed that only two events have been captured at fp8 since 1942, whereas the selective sampling allowed a reconstruction of 23 events (Table 2). During 1995 and 2000 a particularly high occurrence of events was registered at many of the investigated tracks.
By making a linkage between trees showing synchronic avalanche reactions we can reconstruct five different runouts at av6, and four of them correspond to more than one year (Fig. 6). Events with a shorter runout compared to the runout of 2011 could not be reconstructed due to the absence of trees.
Avalanche simulations at avalanche track 9 with a 10-year return
period:
In order to evaluate the influence of the
Avalanche simulations at avalanche track 9 with a 100-year return
period:
Influence of additional afforestation in avalanche track 6 with a
10-year return period:
Simulations of avalanches considering artificial afforestation scenarios were conducted for av6, which shows the largest potential impact on infrastructure along the road (see also Fig. 6), considering a 10-year return period scenario. Our results show that the runout of avalanches would not be significantly shortened (5 %) due to afforestation at the track, but that there would be a reduction in maximum avalanche pressure at the intersection with the road of 24 % (Fig. 9).
Snow avalanches and debris flows are natural hazards that have not been
thoroughly studied in Chile (Espinosa et al., 1983). In this investigation
we identified the main natural hazards representing a risk to local
communities at Valle Las Trancas within the biosphere reserve Nevados de
Chillán – Laguna del Laja. Besides the continuous growth of its human
population, the valley receives a large and increasingly growing number of
visitors attracted by leisure activities and the magnificent mountain
landscape throughout the year. Yet, the lack of land-use planning and
regulation on natural hazards such as snow
avalanches and debris flows contribute to risk to people and
infrastructure. This study assessed the role of native mountain forests,
represented here by different species of
The lack of documentary records of spatiotemporal patterns of gravitational natural hazards in the Southern Andes is a recognized limitation to estimating related risks. Efforts have been made in the past to apply tree-ring methods to this region in order to substitute for or complement documentary records, which, even if available, tend to be fragmentary (e.g., Mundo et al., 2007; Casteller et al., 2009, 2011, 2015). Methodological limitations exist, however: for example, a reconstructed avalanche event year implies the occurrence of at least one event during that particular winter, but it is not possible to determine whether more than one event occurred at the same track during the same winter. In contrast, distinguishing among events occurring in different times of the year is possible by identifying the position of reaction wood within tree rings. In regions where debris flows occur in summer, it is normally possible to distinguish their associated tree-ring reactions (from early wood to late wood) from those originating from snow avalanche impacts (the reaction can be found in this case at the boundary of two rings; Stoffel et al., 2006). In our study area, we were able to distinguish between avalanche and debris flow events using tree-ring methods as snow avalanches occur in winter and debris flows, in turn, typically occur during the early spring when rain events melt the snow from the winter remaining on the catchments. In addition, available records, types of damage in the forest, topographical features (including channel geometry) and process modeling helped us to distinguish one process from another.
As we have shown, different sampling techniques can lead to different results when analyzing temporal patterns of snow avalanches and debris flows. The selection of single trees with widespread visible damage from past events at the borders (and for snow avalanches also at the runouts) of the tracks can provide longer event reconstructions compared to sampling in plots adjacent to disturbed areas. Yet, sampling plots, as we have shown, can also provide valuable data related to forest structure, from which it is possible to infer past event dates. In addition, forest plots are a valuable source of information and can be analyzed with remote-sensing tools and for larger areas. A combination of both sampling approaches was found to be the most useful.
The role of mountain forests as protective measures against natural hazards has been acknowledged for centuries in the European Alps, where mountain settlements are in close contact with phenomena such as snow avalanches and landslides (Price and Thompson, 1997; Getzner et al., 2017. The use of the new model RAMMS::AVALANCHE (which considers forest effects), developed and first tested in the Alps, was used in this investigation to reconstruct past avalanche events and to determine potential avalanche scenarios. These simulations, with a few exceptions (Casteller et al., 2008), are practically nonexistent in the Andes. This type of multi-methodological approach has proven to be most effective (Gądek et al, 2017), in particular in regions where documentary data on past events is scarce. In addition to the methods and data used in this study, the availability of detailed data, e.g., on the size and precise location of avalanche release zones (including snow depth), types of snow and snow entrainment along the track would allow even more robust event modeling.
Our results indicate that local native broadleaf forests may influence avalanche dynamics, particularly in small to medium-sized events if the distance between the release area of the avalanche and the upper timberline is less than ca. 200 m (Weir, 2002; Teich et al., 2012). The effect of broadleaf forests on avalanche releases is certainly less than that of evergreen forests due to lower interception rates during snowfall (Lundberg et al., 2004; Bebi et al., 2009). However, their effect on reducing runout distances may still be considerable, in particular for small avalanches that do not cause widespread tree breakage (Feistl et al., 2014). In such cases, a shortening of runout distances and reduction of impact pressure can largely be explained by the extraction of snow mass from the avalanche volume and thus the reduction of the momentum of the avalanche flow by the trees. These results suggest that appropriate and spatially well-targeted afforestation may potentially reduce risks related to snow avalanches in our study area. In contrast to grey measures such as snow sheds (Campbell et al., 2007), afforestation would clearly have a less detrimental effects on the provision of other ecosystem services, in particular if they are conducted with regional tree species and locally limited to appropriate site conditions (Schönenberger, 2001). In line with the expectations of the local community (Cortés Donoso, 2017), strict controls of further forest cutting within the biosphere reserve are needed in order to allow sustainable protection against natural hazards without major impacts on the forest structure and ecosystem dynamics (Bebi et al., 2017).
As indicated by gathered data at the study site (Cortés Donoso, 2017), the two most relevant risks for the local community at Valle Las Trancas are snow avalanches and volcanic activity. The large volcanic complex has recently reached over 18 months of visible activity. In addition, the occurrence of forest fires in Chile has increased exponentially during recent years, representing a major national hazard (González et al., 2011). While the potential effects of Eco-DRR measures against volcano activity are very limited, they represent a window of opportunity in terms of protection against snow avalanches. The promotion of services provided by the native forest and the participation of the local community in conserving and sustainably managing the forest are issues to be closely considered.
Our results show that both avalanches and debris flows can affect the
structure of
Our research data will be uploaded in the coming weeks to the new data portal
Envrionmental Data Portal EnviDat (ENVIDAT
The authors declare that they have no conflict of interest.
The authors thank the Ministry of Environment of Chile (MMA) for its political support throughout the implementation of the project and the International Climate Initiative (ICI) of the German Federal Ministry for Environment, Nature Conservation, Building and Nuclear Safety (BMUB) for funding the project. Gabriela Jara Aburto and Gabriel Orozco from SERNAGEOMIN are acknowledged for providing support and valuable information from the beginning to the end of the project. Field assistants Ernesto Corvalán, Jorge Silva, Eugenia Marcotti, Carolina Córdova, Viola Debus, Álvaro Gutiérrez and Paul Szejner permitted field activities to be successfully conducted. Likewise, our gratitude goes to Jaime Soto and Leandro Olivares, who conducted drone flights at the study area and shared this information with us for the avalanche modeling. Last but not least, we thank Marc Christen, Yves Buehler and Lukas Stoffel at the SLF Institute for their valuable insight related to the avalanche modeling. Edited by: Thom Bogaard Reviewed by: two anonymous referees