Differences

This shows you the differences between two versions of the page.

--- wiki:site-protecting_forest [2018/07/12 09:30] – eurac
+++ wiki:site-protecting_forest [2018/07/13 15:14] (current) – eurac
@@ Line 1: / Line 1: @@
 ====== Protection of areas against avalanches, mudslides and rockfalls - Supply, Flow, Demand ======
-====    General description:    ====
+==== General description: ====
-<font inherit/inherit;;windowtext;;inherit>The following model is an approach to delineate areas where the biotic ecosystems (forests) contribute to the mitigation of natural hazards and the protection of human assets from hazardous natural processes. At the Alpine-wide scale, this has been done by combining separate regional models for avalanches, rock-falls and water channel relevant processes. The models use topographical information derived from the 25m EU DEM and statistically derived threshold values identified by earlier projects and assessments (performed in all Alpine countries) to model potential avalanche and rock-fall release and transition zones.</font>**<font inherit/inherit;;windowtext;;inherit>Figure</font>**  **<font inherit/inherit;;windowtext;;inherit>1</font>**<font inherit/inherit;;windowtext;;inherit></font><font inherit/inherit;;windowtext;;inherit>describes</font><font inherit/inherit;;windowtext;;inherit>in</font><font inherit/inherit;;windowtext;;inherit>detail</font><font inherit/inherit;;windowtext;;inherit>the</font><font inherit/inherit;;windowtext;;inherit>calculation</font><font inherit/inherit;;windowtext;;inherit>procedure</font><font inherit/inherit;;windowtext;;inherit>to derive the developed indicators per local administrative units (LAU2) of the Alpine Space.</font>
+The following model is an approach to delineate areas where the biotic ecosystems (forests) contribute to the mitigation of natural hazards and the protection of human assets from hazardous natural processes. At the Alpine-wide scale, this has been done by combining separate regional models for avalanches, rock-falls and water channel relevant processes. The models use topographical information derived from the 25m EU DEM and statistically derived threshold values identified by earlier projects and assessments (performed in all Alpine countries) to model potential avalanche and rock-fall release and transition zones. **The flowchart below** describes in detail the calculation procedure to derive the developed indicators per local administrative units (LAU2) of the Alpine Space.
+==== Input data ====
-====      Input data      ====
+  * DEM (slope, slope-length, flow direction, watershed, plan curvature, contour lines) * Land Cover * River Network * Land Cover * River Network
-  * <font inherit/inherit;;windowtext;;inherit>DEM (slope, slope-length, flow direction, watershed, plan curvature, contour lines)</font>* <font inherit/inherit;;windowtext;;inherit>Land Cover</font>* <font inherit/inherit;;windowtext;;inherit>River Network</font>  * <font 14px/calibri,sans-serif;;windowtext;;inherit>Land Cover</font>
+==== Calculation processes ====
-  * <font 14px/calibri,sans-serif;;windowtext;;inherit>River Network</font>  \\
-====      Calculation processes     ====
+=== Avalanche release and Transition zones ===
-=== **__<font 12.0pt/inherit;;inherit;;inherit>Avalanche release and Transition zones</font>__** ===
 **(1) Calculate potential avalanche release areas:**<font inherit/inherit;;windowtext;;inherit>Select all areas with a slope between 28° and 55°, a plan curvature between -2 and 0.2, a forest land-cover and if there is no better information on potential snow-accumulation available above an altitude threshold (Bauerhansel et al 2009).</font>
@@ Line 22: / Line 20: @@
 **(4) Calculate energy line angle for every release area:**<font inherit/inherit;;windowtext;;inherit>The energy line,</font><font inherit/inherit;;windowtext;;inherit>or</font><font inherit/inherit;;windowtext;;inherit>generalized</font><font inherit/inherit;;windowtext;;inherit>gradient</font><font inherit/inherit;;windowtext;;inherit>, represents the angle of the connecting line from the release point and the outer edge of the run out zone of a mass movement with the horizontal plan. (Bauerhansel et al 2009)</font>
-==== Figure 1 : Bauerhansl, C., Berger, F., Dorren, others, 2010. Development of harmonized indicators and estimation procedures for forests with protective functions against natural ====
+====   ====
 <font inherit/inherit;;windowtext;;inherit>In Arc GIS the calculation of the Energy Line Angle can be done with the following steps:</font>
-  * <font inherit/inherit;;windowtext;;inherit>Creating an integer raster with the elevation of the release areas.</font>* <font inherit/inherit;;windowtext;;inherit>Running the tool Euclidean Allocation and Distance with the integer raster as Input (to speed up the calculation, a maximum distance of 2000m should be set).</font>* <font inherit/inherit;;windowtext;;inherit>Using the Raster Calculator, subtract the result from the DEM. The value you get is the elevation difference from every pixel to the release zones Δh.</font>* <font inherit/inherit;;windowtext;;inherit>The Euclidean Distance calculated in step (b) represents Δl.</font>* <font inherit/inherit;;windowtext;;inherit>The energy line angle can than simply be calculated with ELA = atan( Δh/Δl) (</font><font inherit/inherit;;windowtext;;inherit>à</font><font inherit/inherit;;windowtext;;inherit>Arc</font><font inherit/inherit;;windowtext;;inherit>GIS</font><font inherit/inherit;;windowtext;;inherit>calculates</font><font inherit/inherit;;windowtext;;inherit>angles in radians to get a result in degrees you have to multiply it with 57.2958).</font>
+  * <font 14px/inherit;;windowtext;;inherit>Creating an integer raster with the elevation of the release areas.</font>
+  * <font 14px/inherit;;windowtext;;inherit>Running the tool Euclidean Allocation and Distance with the integer raster as Input (to speed up the calculation, a maximum distance of 2000m should be set).</font>
+  * <font 14px/inherit;;windowtext;;inherit>Using the Raster Calculator, subtract the result from the DEM</font><font inherit/inherit;;windowtext;;inherit>. The value you get is the elevation difference from every pixel to the release zones Δh.</font>* <font 14px/inherit;;windowtext;;inherit>The Euclidean Distance calculated in step (b) represents Δl.</font>
+  * <font 14px/inherit;;inherit;;inherit>The energy line angle can than simply be calculated with ELA = atan( Δh/Δl) (→Arc GIS calculates angles in radians to get a result in degrees you have to multiply it with 57.2958)</font>
 **(****5) Limit avalanche Paths with energy line angle:**<font inherit/inherit;;windowtext;;inherit>Knowing the avalanche path, the energy line angle allows a rough estimation of maximal runout distance. The previously measured events occurred at an ELA of 17° - 47° with the mean 28°. Events with ELA 17° are very improbable. We used this threshold for our approach to consider all possible events. (Bauerhansel et al. 2009, PARAmount Project 2012).</font>
-===     Rock-fall Release and Transition Zones     ===
+=== Rock-fall Release and Transition Zones ===
+**(6) Calculate potential rock-fall start zones:**<font inherit/inherit;;windowtext;;inherit>Similar to the avalanche release areas, this can roughly be done by selecting all areas with bare-rock land-cover and a slope steeper than 43°. In the Alpine-wide approach we had to use CORINE Land Cover as the input dataset where most rock surfaces are not delineated (because of resolution and minimum mapping width). For that reason, we also included all the areas without vegetation. Depending on the DEM resolution, different slope thresholds should be applied (Berger et al 2010).</font>
+|DTM resolution/cell size [m]|Threshold slope gradient [º]|
+|1|55|
+|5|49|
+|10|46|
+|25|43|
+**(7) Create Cost Raster:**A cost raster represents the theoretical costs a falling rock would have to overcome from one area to the next (pixel to pixel). These costs depend on the slope length and the land-cover/surface roughness. The highest costs are found in flat, forested areas and low costs exist in steep areas without forest-cover. A single raster can be created by combining two reclassified slope rasters (one for forest areas and one for all other surfaces). The reclassification should occur with the following values:
+^Slope^Costs no forest^Costs in forest|
+|<5º|10|30|
+|>5-10°|9|27|
+|>10-15°|8|24|
+|>15-20°|7|21|
+|>20-25°|6|18|
+|>25-30°|5|15|
+|>30-35°|4|12|
+|>35-40°|3|9|
+|>40-45°|2|6|
+|>45°|1|3|
+**(8) Calculate energy line angle for every start zone:**Energy line angle is calculated using the same methodology as in step 4.
+**(9) Perform Path Distance ****analysis****to**** calculate maximum runout lengths of falling-rocks**: The threshold of 2250 was used which is approximatly a runout of 250 meters in forested terrain (assuming the same dbh and tree density in all forests). With better information on the tree stands in forests different thresholds could be applied using the Rockfor<sup>NET</sup>   tool, found here: [[http://www.ecorisq.org/rockfor-net-en|http://www.ecorisq.org/rockfor-net-en]] .
+=== Areas relevant for water channels ===
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(6) Calculate potential rock-fall start zones:</font>**<font inherit/calibri,sans-serif;;windowtext;;inherit>Similar to the avalanche release areas, this can roughly be done by selecting all areas with bare-rock land-cover and a slope steeper than 43°. In the Alpine-wide approach we had to use CORINE Land Cover as the input dataset where most rock surfaces are not delineated (because of resolution and minimum mapping width). For that reason, we also included all the areas without vegetation. Depending on the DEM resolution, different slope thresholds should be applied (Berger et al 2010).</font>
+**(10) Delineate Forest areas with a protective effect for water-channels: **  All areas with an on-surface distance to water channels of between 10m (trees too near to the channel can have negative influence) and 200m and a slope steeper than >20° are potential erosion and mudslide risk zones that could clog up the river. Forests within this zone are thus particularly important in preventing erosion and events.
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(7) Create Cost Raster:</font>**<font inherit/calibri,sans-serif;;inherit;;inherit>A cost raster represents the theoretical costs a falling rock would have to overcome from one area to the next (pixel to pixel). These costs depend on the slope length and the land-cover/surface roughness. The highest costs are found in flat, forested areas and low costs exist in steep areas without forest-cover. A single raster can be created by combining two reclassified slope rasters (one for forest areas and one for all other surfaces). The reclassification should occur with the following values:</font>
+=== Supply Indicator creation ===
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(8) Calculate energy line angle for every start zone:</font>**<font inherit/calibri,sans-serif;;inherit;;inherit>Energy line angle is calculated using the same methodology as in step 4.</font>
+**(11) Combine areas where forests have a protective function: **By merging all natural hazard zones that already have a forest land-cover, the sites of protection forests are identified. By definition, the protective function of forests exists in all avalanche release and transition zones, rock-fall transition zones and water channel relevant terrain.
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(9) Perform Path Distance analysis</font>** [[http://www.ecorisq.org/rockfor-net-en|**to calculate maximum runout lengths of falling-rocks**: The threshold of 2250 was used which is approximatly a runout of 250 meters in forested terrain (assuming the same dbh and tree density in all forests).  With better information on the tree stands in forests different thresholds could be applied using the Rockfor<sup>NET</sup> tool, found here: http://www.ecorisq.org/rockfor-net-en]] .
+=== Demand indicator creation ===
-===      Areas relevant for water channels     ===
+**(12) Calculate natural hazard potential:**Merge the spatial extent of the avalanche and rock-fall transition zones calculated above.
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(10) Delineate Forest areas with a protective effect for water-channels:</font>**<font inherit/calibri,sans-serif;;inherit;;inherit>All areas with an on-surface distance to water channels of between 10m (trees too near to the channel can have negative influence) and 200m and a slope steeper than >20° are potential erosion and mudslide risk zones that could clog up the river. Forests within this zone are thus particularly important in preventing erosion and events.</font>
+**(13) Calculate Damage Potential:**Create a raster with all human infrastructure, such as settlements areas, built up areas, buildings, roads and railways.
-===      Supply Indicator creation     ===
+**(14) Calculate protection forest demand index:**Overlay the natural hazard potential and the damage potential areas. This allows one to identify all of the Infrastructure in potential hazard zones that is in need of protection.
-<font inherit/calibri,sans-serif;;inherit;;inherit>**(11) Combine areas where forests have a protective function: **By merging all natural hazard zones that already have a forest land-cover, the sites of protection forests are identified. By definition, the protective function of forests exists in all avalanche release and transition zones, rock-fall transition zones and water channel relevant terrain.</font>
+=== Flow indicator creation ===
-===      Demand indicator creation     ===
+**(15) Backlink intersections of hazard potential and damage potential to release areas:**In the case of avalanches all the transition zones that intersect with infrastructure have to be traced back to the corresponding release area. Using Raster in ArcGIS, this can be done by inverting the DEM and performing a new cost path analysis starting from the intersected damage potential areas. Finally, by combining the avalanche release and transition zones as well as rockfall transition zones that possibly harm infrastructure we are able to identify Object-protecting forest.
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(12) Calculate natural hazard potential:</font>**<font inherit/calibri,sans-serif;;inherit;;inherit>Merge the spatial extent of the avalanche and rock-fall transition zones calculated above.</font>
+<font 26px/inherit;;inherit;;inherit>**Protection Forest**</font>
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(13) Calculate Damage Potential:</font>**<font inherit/calibri,sans-serif;;inherit;;inherit>Create a raster with all human infrastructure, such as settlements areas, built up areas, buildings, roads and railways.</font>
+{{:en:protection_forest_flowchart_f1.jpg?nolink&1643x2001}}
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(14) Calculate protection forest demand index:</font>**<font inherit/calibri,sans-serif;;inherit;;inherit>Overlay the natural hazard potential and the damage potential areas. This allows one to identify all of the Infrastructure in potential hazard zones that is in need of protection.</font>
+{{:en:biomassproductionfromgrassland_page_12_1_.jpg?nolink&500x300}}
-===      Flow indicator creation     ===
+<font 18px/inherit;;inherit;;inherit>**Flowchart depicting the procedures used to derive the supply, flow and demand indicators**</font>
-**<font inherit/calibri,sans-serif;;inherit;;inherit>(15) Backlink intersections of hazard potential and damage potential to release areas:</font>**<font inherit/calibri,sans-serif;;inherit;;inherit>In the case of avalanches all the transition zones that intersect with infrastructure have to be traced back to the corresponding release area. Using Raster in ArcGIS, this can be done by inverting the DEM and performing a new cost path analysis starting from the intersected damage potential areas. Finally, by combining the avalanche release and transition zones as well as rockfall transition zones that possibly harm infrastructure we are able to identify Object-protecting forest.</font>
+\\