comparing empirical models, 2d- and 3d process-based ...€¦ · different modeling results in...

1
Geological Survey of Austria - Depart. of Engineering Geology Neulinggasse 38, A-1030 Vienna phone: +43-1-7125674-0 Fax: +43-1-7125674-56 393 [email protected] www.geologie.ac.at European Geosciences Union General Assembly 2010 Austria, Vienna, 02- 07 May 2010 Comparing empirical models, 2D- and 3D process-based models for delineating maximum rockfall runout distances Sandra Melzner*, Michael Mölk**, Luuk Dorren° and Richard Bäk°° *Geological Survey of Austria,**Austrian Service for Torrent and Avalanche Control,°ecorisQ,°°Federal State Government of Carinthia N Austria Carinthia Vienna Local Extent Regional Extent Slope Extent (Green line= energy angle lines: red dots= maximum runout distance; blue polygone= possible conflict) (Red Polygon= detailed parameter maps Rockyfor3D; blue lines= detailed input parameters Rockfall 7.1 ) (Red lines= Rockyfor3D trajectories; blue lines=output Rockfall 7.1; red dots= maximum runout distance; black dots= possible endangered area) Energy line approach (“Geometrisches Gefälle”) “Fahrböschung” File name (raster) Describtion E_mean.asc Mean kinetic energy (translational & rotational; inkJ) per cell E_95.asc The 95% confidence level of all kinetic energy values (in kJ) in each cell Ph_mean.asc Mean passage height (in m; measured in normal direction to the slope surface) per cell Ph_95.asc The 95% confidence level of all passage height values in each cell Nr_passage.asc The number of blocks passed through each cell Nr_deposited.asc The number of blocks passed through each cell Rvol_deposit.asc The maximum block volume (in m³) stopped in each cell EL_angles.asc Mean energy line angles per cell (in-) OUTPUT “Rockyfor3D” OUTPUT “Rockfall 7.1” * analysing model sensitivity to the accuracy of input data * comparing model results at regional, local and slope scale * comparing Rockyfor3D with the 3D model STONE (Guzzetti et al. 2002) Future activities Strategy We gratefully thank Hannes Burger (Austrian Service for Torrent and Avalanche control), Fausto Guzzetti & Paola Reichenbach (CNR-IRPI, Perugia, Italy) and the township of Bad Bleiberg for their support in carrying out this work. Acknowledgments References Dorren (2009): Rockyfor3D revealed. Description of the complete 3D rockfall model, ecorisQ-publication. Guzzetti et al. (2002): STONE: a computer programm for the three- dimensional simulation of rock falls, Com.& Geosc. 28: 1079-1093. Mölk (2008):Steinschlagschutz Bad Bleiberg. Steinschlagsimulation als Basis für die Dimensionierung von Schutzmaßnahmen. Unveröfftl. WLV Bericht, Innsbruck, Österreich Melzner et al. (2009a):Erstellung einer detaillierten Gefahrenhinweiskarte und ingenieurgeologische Ausarbeitung relativer Gefahrenklassen für ausgewählte Gebiete Kärntens im Maßstab 1:10.000. Testgebiet Bad Bleiberg- Prozessgruppe Fallen/Stürzen. Unveröfftl. GBA Endbericht, Wien, Österreich. Melzner et al. (2009b): Entwurf zur Entwicklung einer Bearbeitungsstrategie für die Einschätzung des Gefahrenpotentials durch Sturzprozesse im Rahmen von Regionalstudien.Unveröfftl. GBA Bericht, Wien, Österreich. Step 1: Identification of potential conflicts within large regions (regional extent) Step 2: Specification of a potential conflict in some parts of the region (local extent) Step 3: Proposal of protective measures for a possible conflict for some slopes (slope extent) To define maximum rockfall runout distances at regional scale, usually simple methods are used (e.g. energy line principle, shadow angle or “Fahrböschung”), due to the lack of high resolution and accurate slope cover and type information for larger areas. A regional simulation with Rockyfor3D was tested, using greatly generalised parameter maps, which were prepared by calculating mean values from samples of a representative number of punctual information of mapped parameters in predefined homogeneous classes. Overlaying these results with areas of habitations and infrastructure could identify potential conflict areas. Comparing different regionally applied methods for delineating potentially endangered areas on local scale. These potentially endangered areas are then investigated in more detail to pin down the potential conflict (if, where and to what extent a conflict exists). For the 2D simulation the necessary input parameters were acquired on representative (adverse) profiles in previously determined homogeneous areas. The input parameters for Rockyfor were acquired by spatially continuous field mapping. Comparing different locally applied methods for delineating potentially endangered areas on slope extent. In case these results reveal an imminent danger for certain areas, more detailed investigations should follow up. Only on basis of such detailed studies at slope scale a realistic quantitative hazard assessment can be performed. This is important when it comes to planning protective measures (e.g. what type and dimension are required) or adaption/designation of hazard zones. : Comparing different models regarding their capacity for planning purposes and dimensioning of preventive measures. Goal: Goal: Goal Introduction An important step in rockfall hazard analysis is the delineation of potential maximum runout distances. Different approaches exist to determine the runout zones of rock falls varying from emperical models to 2D and 3D process- based models. Depending on the research goal, spatial scale, and quality of available input data, the appropriate model approach has to be applied according to the requested data validity domain of the output data. The present work compares the capacity of empirical models and a 2D resp. 3D process-based model (”Rockfall 7.1” resp. “Rockyfor3D”) to predict rockfall runouts at different spatial scales. Special focus is on the analysis of the effects of the quality of the model input data on the modelling results. This is conducted by a step-wise downscaling from the region scale to the slope scale, which brings about an increasing quality of the accompanying model input maps. A highly rock fall prone area in Carinthia (Austria) was chosen for a model comparison, also because of the fact that a large amount of field data about past rockfall events and rock fall determining factors was already available (see Mölk 2008, Melzner et al. 2009a&b). This allows an evaluation of the different model types and a more accurate validation of the different modeling results in respect to the real conditions of the study area. Models/ Input parameters Input parameter „Fahrböschung“ & „Geometrisches Gefälle“ „Rockfall 7.1“ (see Dr. Spang GmbH 2008) „Rockyfor3D“ (see Dorren 2009) Source areas v v v DEM v -- v Topography -- v -- Rock density -- v v Shape & size of the falling block -- v v Surface roughness -- v v Surface elasticity -- v v Dynamic friction coefficient for rolling -- v -- Vegetation -- (v) v Starting velocity -- v v Initial falling heights -- -- v First results Photo by S. Melzner, GBA dimension: 2000 kJ; heights: 4m black dots= mapped boulders white squares=houses Output “E_mean” Output “Ph_mean” Output “Nr_deposited” red dots= mapped boulders white squares=houses green lines= dimensioned protective measure Photo by S. Melzner, GBA Methods level of detail/quality of input data/ data validity domain Spatial scale Regional extent low medium high Slope extent Local extent e.g. empirical models e.g. „Rockyfor3D“ e.g. „Rockyfor3D“ & „Rockfall 7.1“ downscaling downscaling

Upload: others

Post on 18-Oct-2020

1 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Comparing empirical models, 2D- and 3D process-based ...€¦ · different modeling results in respect to the real conditions of thestudyarea. Models/ Input parameters Input parameter

Geological Survey of Austria - Depart. of Engineering GeologyNeulinggasse 38, A-1030 Vienna

phone: +43-1-7125674-0 Fax: +43-1-7125674-56 [email protected] www.geologie.ac.at

European Geosciences UnionGeneral Assembly 2010

Austria, Vienna, 02- 07 May 2010

Comparing empirical models, 2D- and 3D process-based modelsfor delineating maximum rockfall runout distancesSandra Melzner*, Michael Mölk**, Luuk Dorren° and Richard Bäk°°

*Geological Survey of Austria,**Austrian Service for Torrent and Avalanche Control,°ecorisQ,°°Federal State Government of Carinthia

N

Austria

Carinthia

Vienna

Lo

calE

xte

nt

Reg

ion

alE

xte

nt

Slo

pe

Exte

nt

(Green line= energy angle lines: reddots= maximum runout distance; bluepolygone= possible conflict)

(Red Polygon= detailed parameter mapsRockyfor3D; blue lines= detailed inputparameters Rockfall 7.1 )

(Red lines= Rockyfor3D trajectories;blue lines=output Rockfall 7.1; reddots= maximum runout distance;black dots= possible endangeredarea)

Energy line approach(“Geometrisches Gefälle”)

“Fahrböschung”

File name

(raster)

Describtion

E_mean.asc Mean kinetic energy (translational & rotational; inkJ) percell

E_95.asc The 95%confidence level of all kinetic energy values (in kJ) in

each cell

Ph_mean.asc Mean passage height (in m; measured in normaldirection to the

slope surface) percell

Ph_95.asc The 95%confidence level of all passage height values in each

cell

Nr_passage.asc The number of blocks passed through each cell

Nr_deposited.asc The number of blocks passed through each cell

Rvol_deposit.asc The maximum block volume (in m³) stoppedin each cell

EL_angles.asc Mean energy line anglesper cell (in-)

OUTPUT “Rockyfor3D”

OUTPUT “Rockfall 7.1”

* analysing model sensitivity to the accuracy of input data

* comparing model results at regional, local and slope scale

* comparing Rockyfor3D with the 3D model STONE(Guzzetti et al. 2002)

Future activities

Strategy

We gratefully thank Hannes Burger (Austrian Service forTorrent and Avalanche control), Fausto Guzzetti & PaolaReichenbach (CNR-IRPI, Perugia, Italy) and the township ofBad Bleiberg for their support in carrying out this work.

AcknowledgmentsReferencesDorren (2009): Rockyfor3D revealed. Description of the complete 3D rockfall model, ecorisQ-publication.Guzzetti et al. (2002): STONE: a computer programm for the three- dimensional simulation of rock falls, Com.& Geosc. 28:1079-1093.Mölk (2008):Steinschlagschutz Bad Bleiberg. Steinschlagsimulation als Basis für die Dimensionierung vonSchutzmaßnahmen. Unveröfftl. WLV Bericht, Innsbruck, ÖsterreichMelzner et al. (2009a):Erstellung einer detaillierten Gefahrenhinweiskarte und ingenieurgeologische Ausarbeitungrelativer Gefahrenklassen für ausgewählte Gebiete Kärntens im Maßstab 1:10.000. Testgebiet Bad Bleiberg-Prozessgruppe Fallen/Stürzen. Unveröfftl. GBAEndbericht, Wien, Österreich.Melzner et al. (2009b): Entwurf zur Entwicklung einer Bearbeitungsstrategie für die Einschätzung des Gefahrenpotentialsdurch Sturzprozesse im Rahmen von Regionalstudien.Unveröfftl. GBABericht, Wien, Österreich.

Step 1: Identification of potential conflicts within largeregions (regional extent)

Step 2: Specification of a potential conflict in some partsof the region (local extent)

Step 3: Proposal of protective measures for a possibleconflict for some slopes (slope extent)

To define maximum rockfall runout distances at regionalscale, usually simple methods are used (e.g. energy lineprinciple, shadow angle or “Fahrböschung”), due to the lackof high resolution and accurate slope cover and typeinformation for larger areas. A regional simulation withRockyfor3D was tested, using greatly generalised parametermaps, which were prepared by calculating mean values fromsamples of a representative number of punctual informationof mapped parameters in predefined homogeneous classes.Overlaying these results with areas of habitations andinfrastructure could identify potential conflict areas.

Comparing different regionally applied methods fordelineating potentially endangered areas on local scale.

These potentially endangered areas are then investigated inmore detail to pin down the potential conflict (if, where and towhat extent a conflict exists).For the 2D simulation the necessary input parameters wereacquired on representative (adverse) profiles in previouslydetermined homogeneous areas. The input parameters forRockyfor were acquired by spatially continuous fieldmapping.

Comparing different locally applied methods fordelineating potentially endangered areas on slope extent.

In case these results reveal an imminent danger for certainareas, more detailed investigations should follow up. Only onbasis of such detailed studies at slope scale a realisticquantitative hazard assessment can be performed. This isimportant when it comes to planning protective measures(e.g. what type and dimension are required) oradaption/designation of hazard zones.

: Comparing different models regarding their capacity forplanning purposes and dimensioning of preventivemeasures.

Goal:

Goal:

Goal

IntroductionAn important step in rockfall hazard analysis is thedelineation of potential maximum runout distances. Differentapproaches exist to determine the runout zones of rock fallsvarying from emperical models to 2D and 3D process- basedmodels. Depending on the research goal, spatial scale, andquality of available input data, the appropriate modelapproach has to be applied according to the requested datavalidity domain of the output data.The present work compares the capacity of empirical modelsand a 2D resp. 3D process-based model (”Rockfall 7.1” resp.“Rockyfor3D”) to predict rockfall runouts at different spatialscales. Special focus is on the analysis of the effects of thequality of the model input data on the modelling results. Thisis conducted by a step-wise downscaling from the regionscale to the slope scale, which brings about an increasingquality of the accompanying model input maps.Ahighly rock fall prone area in Carinthia (Austria) was chosenfor a model comparison, also because of the fact that a largeamount of field data about past rockfall events and rock falldetermining factors was already available (see Mölk 2008,Melzner et al. 2009a&b). This allows an evaluation of thedifferent model types and a more accurate validation of thedifferent modeling results in respect to the real conditions ofthe study area.

Models/ Input parameters

Input parameter „Fahrböschung“

&

„Geometrisches

Gefälle“

„Rockfall 7.1“

(see Dr. Spang

GmbH 2008)

„Rockyfor3D“

(see Dorren

2009)

Source areas v v v

DEM v -- v

Topography -- v --

Rock density -- v v

Shape & size of the

falling block-- v v

Surface roughness -- v v

Surface elasticity -- v v

Dynamic friction

coefficient for

rolling

-- v --

Vegetation -- (v) v

Starting velocity -- v v

Initial falling

heights-- -- v

First results

Photo by S. Melzner, GBA

dimension: 2000 kJ; heights: 4m

black dots= mapped boulderswhite squares=houses

Output “E_mean”

Output “Ph_mean”

Output “Nr_deposited”

red dots= mapped boulderswhite squares=housesgreen lines= dimensionedprotective measurePhoto by S. Melzner, GBA

Methods

level of detail/quality ofinput data/

data validity domain

Sp

ati

al

sc

ale

Regionalextent

low medium high

Slopeextent

Localextent

e.g.empiricalmodels

e.g.„Rockyfor3D“

e.g.„Rockyfor3D“& „Rockfall

7.1“

downscaling

downscaling