{'fraction_non_alphanumeric': 0.027016414042252416, 'fraction_numerical': 0.009777742872850074, 'mean_word_length': 4.965448582806278, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Modelling woody material transport and deposition in alpine rivers

Recent ﬂood events in Switzerland and Western Austria in 2005 were characterised by an increase in impacts and associated losses due to the transport of woody material. As a consequence, protection measures and bridges suffered considerable dam-ages. Furthermore, cross-sectional obstructions due to woody material entrapment caused unexpected ﬂood plain inundations resulting in severe damage to elements at risk. Until now, the transport of woody material is neither sufﬁciently taken into account nor systematically considered, leading to prediction inaccuracies during the procedure of hazard mapping. To close this gap, we propose a modelling approach that (1) allows the estimation of woody material recruitment from wood-covered banks and ﬂood plains; (2) allows the evaluation of the disposition for woody material entrainment and transport to selected critical conﬁgurations along the stream and that (3) enables the delineation of hazard process patterns at these critical conﬁgurations. Results from a case study suggest the general applicability of the concept. This contribution to woody material transport analysis reﬁnes ﬂood hazard assessments due to the consideration of woody material transport scenarios.


Introduction
Socio-economic developments in European mountain environments and related forelands are reflected in increasing settlement and economic activities in areas affected by natural hazards (Fuchs and Holub 2007). Consequently, considerable economic losses have resulted in recent years from events (Mitchell 2003;Autonome Brovinz Bozen-Südtirol 2006;Oberndorfer et al. 2007) despite the efforts made towards the mitigation of flood hazards and the reduction of specific risks (Fuchs and McAlpin 2005). At critical stream geometry configurations in mountain streams (e.g., bridge cross sections), a remarkable increase in process intensities could be attributed to woody material transport phenomena (e.g., Diehl 1997;Lyn et al. 2007). Due to the general necessity of assessing natural hazards and risks in a reproducible manner, guidelines for hazard mapping were defined in European countries (e.g., BUWAL 1998;Autonome Provinz Bozen-Südtirol 2006), thereby providing milestones with respect to the quality of integral risk management. The major starting point for managing risk from an integral point of view is the deduction and systematic construction of consistent and reliable scenarios. However, subjective assumptions on relevant impact variables such as woody material transport intensities on the system-loading side and response mechanisms at critical configurations often cause biases and inaccuracies in the results. Considering the hazardous effects of woody material transport, clear indications emerged from the analyses of the debris flow and flood events that occurred recently in several Alpine regions (e.g., Bänziger 1990; Rickli and Bucher 2006). At critical stream geometry configurations (e.g., bridge locations), the transported woody material is repeatedly entrapped. In addition to increasing the loading conditions on the structural components of the bridges (e.g., piers, abutments and superstructure), overflowing is likely to occur more frequently. In order to assess these phenomena, this paper aims at contributing to a systematic investigation of woody material recruitment processes, an evaluation of the propensity for woody material entrainment and transport to critical configurations, and a detection of hazard process patterns at these critical configurations during extreme floods. The analysis of such elements is indispensable for comprehensive flood hazard assessments and for optimizing forest management strategies. Moreover, knowledge about the quantity of woody material, the main woody material pathways in the stream channel and the main places of deposition is fundamental for the design of resilient protection systems and for optimised emergency planning. In order to approach problems emerging from in-stream structures, a detailed study of the damage potential has been carried out by Diehl (1997). The relation between the flow resistance due to the presence of large woody debris (LWD) and increased inundation frequency has been analysed in detail by Shields et al. (2001). New insights into the dynamics of wood transport in streams have been achieved by flume experiments (Rickenmann 1997;Braudrick et al. 1997;Grant 2000, 2001;Degetto 2000;Haga et al. 2002;Curran and Wohl 2003;Bocchiola et al. 2006). The interaction of woody material transport with protection measures has been investigated with a special focus on check dams (Bezzola et al. 2004;Lange and Bezzola 2006) and on rope nets (Rimböck and Strobl 2002).
Acknowledging the fact that hazard impacts at critical configurations along the stream could be interpreted as effects of a complex process interaction field, the main objective of this paper is to propose a modelling concept for the analysis of the following key aspects in considering woody material transport in flood hazard mapping: 1. Disposition: Wood stand productivity and dead wood production in the recruitment areas are important factors that determine the disposition for woody transport phenomena in mountain streams. 2. Intensity of wood-flood interaction: Recruitment processes directly connected to the dynamics of wood-flood interaction become relevant. Recruitment processes due to wood-flood interaction can be attributed directly to hydrodynamic pressure loading and subsequent breakage of the stems. Moreover, morphodynamics such as stream bed erosion and aggradation and side erosion plays a relevant role (e.g., Abbe et al. 1997;Hildebrand et al. 1997;Gurnell et al. 2000;Keim et al. 2002;Montgomery and Piegay 2003;Comiti et al. 2006). In fact, the erosive action of the current is responsible for the scouring of root wads which in turn induces tree toppling. 3. Entrainment and transport of the woody material: The intensity of the flood process, in terms of flow depths and flow velocities, has to be considered as a critical parameter. 4. Interaction phenomena at critical channel geometry configurations: Significant in this context are woody material entrapment and the related consequences, e.g., bridge failures due to hydrostatic and hydrodynamic overloading. Because of the complexity of the involved process chains, magnitude-frequency-related considerations deserve close attention.
Given earlier, the developed conceptual structure comprises: (1) criteria for the localization and classification of woody material recruitment areas as well as the assessment of the woody material recruitment volumes; (2) a computational procedure for woody material entrainment processes; (3) a computational scheme for woody material transport, deposition and remobilization dynamics and (4) an analysis procedure of the interaction phenomena involving transported woody material occurring at critical stream configurations.
Data were implemented into a GIS environment. Subsequently, we tested the GIS application for reconstructing woody material recruitment, transport and deposition patterns during a design event with a reoccurrence period of 300 years and compared it with a flash-flood event occurred in 1987 in the Passirio/Passer River in South Tyrol (Northern Italy).

Theoretical background
Throughout the paper, we will consider a wood-covered flood plain region as a system, X, as shown in Fig. 1. The system is confined at the downstream side by the outflow boundaries, C out , and at the upstream side by the inflow boundaries, C in . To simplify matters, these boundaries are assumed to be invariant. Furthermore, it is supposed that the material flux exchanges (e.g., discharge, sediment rates and wood material amounts) within the environment are taking place at these boundaries. At the margins, the system is confined to lateral flood plain boundaries, namely the slopes of the mountains. The system consists of stocks or storage compartments and flows or fluxes. Three storage compartments are defined as (1) sediment storage, (X X 1 ) t , (2) water storage, (X X 2 ) t and (3) wood material storage, (X X 3 ) t . The corresponding fluxes within the system and at the inflow boundaries are: (1) sediment fluxes, (D X 1 ) t and (D C 1 ) t ; (2) water fluxes, (D X 2 ) t and (D C 2 ) t ; and (3) woody material fluxes, (D X 3 ) t and (D C 3 ) t . Each flood can be intended as a disturbance of the system, and the effects of a certain flood event depend also on the settings given by the previous floods. As a consequence at the beginning of each flood event, t i 0 with i = 1, … , N, a specific set of initial conditions, X j X À Á t0 i and D j X À Á t0 i , with j = 1, 2, 3 and inflow boundary conditions, D j C À Á t i , have to be considered.
The simulations of the material fluxes taking place within the considered system during successive extreme events would require the solution of distinct boundary conditions and initial value problems. Between successive extreme events, the morphodynamical development as well as the growth of the wood stocks has to be monitored in order to correctly assess both the initial and boundary conditions. In the following subsections, a concise summary of the main findings regarding woody material recruitment, entrainment and transport are given. Rickli and Bucher (2006) identify the following relevant recruitment processes from hill slopes:

Woody material recruitment
1. Bank erosion: through the shear stress exerted on the wetted perimeter of the cross section of the channel, erosion processes occur along the banks and alter the static equilibrium of the trees. Synchronisms between debris flows, flood events and the earlier mentioned woody material recruitment processes are quite probable. 2. Wind-throw: strong wind conditions can either destabilise the trees that consequently fall as a whole into the channel or lead to the recruitment of the epigeous parts if their stems break under the wind loading. May and Gresswell (2003) point out that falling trees with a horizontal distance to the channel that exceeds their height can exert a destabilizing action on other trees (e.g., knock-on effect).  Fig. 1 Wood, sediment and fluxes and system dynamics 3. Snow loading: the pressure exerted by snow, in particular wet snow in spring and autumn, can cause stability problems to broad-leaved trees, increasing bending moments and shear forces. Through the cleaving action of ice, crowns of trees are susceptible to break off. 4. Landslides and other slope processes convey standing and lying woody material towards the channel. 5. Avalanches are likely to incorporate and convey large volumes of woody material within the process perimeter.
Only a few exhaustive experimental results currently exist regarding woody material production and recruitment from different wood structures which are exposed to hydrodynamic impacts. Extensive experimental investigations have been carried out in order to study the hydraulic impact on vegetated riverbanks, e.g., the experiments carried out in the soil bioengineering test flume of the Vienna river (Rauch 2005). Furthermore, a series of 3D numerical simulations of vegetated Compound Channel Flows have been performed (Wilson et al. 2004). A second impact factor influencing woody material production is the wood structure itself. An analysis of recent flood events documented in the Province of Bolzano showed that wood vegetation in the river bed has a higher relative tendency to produce woody material (depending on wood stand volume) compared to wood vegetation on stream banks. In absolute terms, the amount of woody material recruited from the stream banks exceeded the volume recruited from the riverbank (Mazzorana et al. 2009). This fact can be attributed to the wood stand volume being on average larger on the bank slopes than in the river bed. In relation to the wood stand volume, wood vegetation located on the stream banks delivers more woody material due to lateral erosion phenomena than wood vegetation of the flood plain.

Woody material entrainment and transport
Theoretical models for woody material entrainment based on the balance of hydrodynamic (F) and resistance forces (R) on individual large woody material pieces have been developed by Braudrick et al. (1997) and Braudrick and Grant (2000). Assuming the shape of the woody debris pieces being cylindrical, and neglecting the influence of buoyancy, the hydrodynamic force can be expressed as follows (Haga et al. 2002): where C d , drag coefficient for the woody material element in water; q, density of the water; l, length of the woody material element; h, flow depth given by hydrodynamic simulations; A sub , submerged area of the log perpendicular to length; h, angle of the element axis relative to the main flow direction; U, flow velocities given by hydrodynamic simulations as: . The resistance forces can be expressed as follows: where d, diameter of the woody material element; l, friction coefficient between the element and the channel bed; r, density of the woody material element; a, channel bed slope; g, gravity acceleration.
Expressing the submerged area of the log perpendicular to its length can be defined as follows: Expressing the non-dimensional force W ¼ F R in terms of Eqs. 1 and 2, the following expression can be obtained (modified from Haga et al. 2002): Analysing Eq. 4, the dynamics of a single woody material element with known dimensions can be described within a simplified scheme as follows (Haga et al. 2002): Floating condition: Resting condition: Sliding or rolling condition: Field observations are in good agreement with the earlier outlined theory. According to those evidences, the entrainment condition for smooth wood logs with an approximate cylindrical form is 1 h d 1:2. The ratio, h d , increases for wood logs with branches up to 1.5 and for wood logs with root wads up to 1.7.
The earlier outlined theoretical model does not consider the effects of morphodynamics. After Lange and Bezzola (2006), the entrainment is facilitated in the case of sediment transport due to the fact that the movable stream bed layer acts like a roller-bearing. In these conditions, the required flow depths, ceteris paribus, are 20-30% lower compared to the case without sediment transport.
The procedure for evaluating woody material dynamics from a hazard-related perspective relies on these theoretical principles and is described in the next section.

Modelling approach
This section provides a description of the developed modelling approach, which has been implemented into an ArcGIS-Esri Ó environment. The criteria for the identification, localization and classification of the woody material recruitment areas are explained, followed by the procedure for woody material volume assessment. Subsequently, the relevant aspects of woody material transport dynamics, namely the entrainment, transport and deposition and remobilization processes, are thoroughly discussed and a computational scheme is proposed. Finally, an analysis procedure of interactions at critical stream configurations is introduced.

Identification, localization and classification of woody material recruitment areas
Woody material recruitment areas are identified by the interpretation of aerial stereo images (photographs). According to the findings of Rauch (2005), an innovative ad hoc classification of the typologies of alluvial forests, lowland riparian forests and riverside woodland is proposed. The forested areas in the influence zone of the stream are classified into seven structural typologies according to their different behaviour when exposed to hydrodynamic loadings (Table 1). The classification criteria take into account the response of different vegetation and forest typologies to the hydraulic forces and impacts of flood processes. Figure 2 shows the interdependencies between woody material production (interpreted as ''distance'' from the point of origin in the coordinate system), the level of hydrodynamic impact forces (on the vertical axis), the position within the riparian zone (on the horizontal axis) and the vegetation structure (oblique axis). Table 1 also shows the potentially available volumes of woody material for each vegetation structure typology. In addition, Table 1 provides the description of the corresponding typical response mechanisms. The geomorphologic characteristics of the flood process areas have to be determined. These areas are classified as (1) stream bed, (2) stream bank and (3) process area of an extreme flood event.
The classification is made on the basis of digital elevation models, aerial photographs and a Lidar-based digital terrain model with a high spatial resolution. An output data set is obtained with the morphology of the stream influence zone.

Assessment of the recruited woody material volume
Three distinct computational procedures for the assessment of woody material recruitment volumes are presented for: the recruitment from (1) hill slopes, the recruitment from (2) tributaries, and (3) within the maximum extent of the flooded area.

Assessment of the recruitment volumes from hill slopes
In this subsection, a procedure to assess the recruited woody material volumes from hill slopes is outlined. This assessment involves six steps: Step 1 Determination of the perimeter of the flooded areas of the considered extreme event from flow depth or velocity raster data sets. The flow depth or velocity data sets corresponding to the different time steps are given as output files of the hydrodynamic simulations performed with the 2D numerical model Sobek Rural (WL/Delft Hydraulics 2004). These raster data sets are overlaid in order to derive the perimeter identifying the maximum extent of the flood (see Fig. 3 for details).
Step 2 A screening for the availability of wood-covered areas outside of the maximum flood boundaries. These wood-covered areas within a buffer width are categorised as active wood buffer (AWB) areas. This buffer width corresponds either to the potential tree height derived from the forest typology map (Klosterhuber et al. 2007) or to the real tree height retrieved from the Lidar-based digital surface model and the digital terrain model (Mazzorana et al. 2009). The difference between the surface model and the terrain model can give hints about the reference tree height and therefore about the width of the active wood buffer.
Step 3 Identification of the recruitment wood buffer strips. Inside the delimitation polygon of the maximum extent flood area, buffer strips for the wood volumes, potentially recruited form the AWB areas, are identified in direct proximity to the AWB areas (see Fig. 3 for details). The width of these strips is half of the tree height of the adjacent AWB area and corresponds to the area of possible location of the centre of gravity of the recruited wood logs. Step 4 Computation of the length of the contact boundary between the AWB areas and the buffer strips (identified in step 3). The user has to define the volume of recruited woody material per unit length of contact boundary for each AWB area. An estimation of recruitment volume scenarios is made on the basis of experimental data provided by inventories (Rickli and Bucher 2006). It should be noted that uncertainty is not negligible without detailed forest inventories and a precise estimation of dead wood volumes and forest stand stability conditions of the forest population. It is recommended to define a plurality of scenarios of recruitment volumes per unit length.
Step 5 Specification of the dimension and number of the model wood logs. In this step, the dimensions of the ''model wood log'' (diameter and log length) are specified and the determination of the recruited number of model wood logs for each strip is straightforward, once known the volumes from the recruitment scenarios defined in the previous step. The position of the centre of gravity of the ''model wood logs'' is either chosen at random (see Fig. 3) or assessed by field investigations.
Step 6 Determination of the starting position of the ''model wood logs''. The volume of each recruited ''model wood log'' is assigned to the corresponding flow cell depending on the position of its centre of gravity (compare Fig. 3).

Assessment of the recruitment volumes from tributaries
This assessment involves the following two steps: Step 1 The recruitment scenarios of wood material volumes from the tributaries can be reliably defined only if a detailed documentation exists for events. If information is missing, possible scenarios of recruited woody material volumes can be defined on the basis of empirical equations (Rickenmann 1997).
Step 2 The number of ''model wood logs'' is calculated. Again the problem arises of the positioning of the ''model wood logs'' within the maximum extent flood area. In the case of accurate event documentations, the ''model wood logs'' can be placed within the geo-referenced depositional areas (see Fig. 3). Otherwise, engineering judgement is required to position the ''model wood logs'' by choosing, e.g., probable deposition sites.

Assessment of the recruitment volumes within the flooded area
Forests of different wood typologies produce determined amounts of dead wood ready for transport depending on the current stage in their ''life cycle''. During the flood event itself, additional dead wood is produced as a consequence of either direct hydrodynamic impact on trees or destabilization of the trees' anchorage through erosion phenomena (root wad scouring). The vegetation structures are mapped and the morphology of the riverbank is overlaid with the outputs of the respective flood simulation. A woody material recruitment indicator is computed based on an impact-response assessment approach, which considers either the morphological characteristics or the flood intensity. This impact-response assessment approach has been developed on the basis of the findings of Rauch (2005) and Hübl et al. (2008). Experimental studies should further investigate the interplay of the following key factors: (1) Hydrodynamic impact: The hydrodynamic impact acts, on the one hand, through static and dynamic pressure forces on wood vegetation and, on the other hand, by yielding stress on the soil, weakening the root-soil anchorages.
(2) Wood stand volume: A positive correlation between the recruited woody material volume and the wood stand volume is postulated. Average wood stand volume estimations for different wood stand structures are reported in Table 1. (3) Wood vegetation resistance-resilience: The influence of the aforementioned key factors cannot be understood and quantitatively estimated without assessing wood vegetation resistance-resilience mechanisms. These depend on both wood structural and species-specific characteristics. For the necessary accuracy level, the analysis is limited to the first type of characteristics underlining that resistance-resilience against hydrodynamic impact depends significantly on the flexibility of the wood and on recovery capacity. In addition, a very flexible vegetation structure protects the soil from erosion, while inflexible-old growth-population structures are weakened by erosion mechanisms and are also unstable due to an unfavourable slenderness (h/d) ratio. These interaction phenomena are described in Table 1.
The assessment scheme is shown in Fig. 2, where the energy indicator C i,j describes hydraulic impact on vegetation structures (see Eq. 8). The energy indicator is composed of a hydrostatic and a hydrodynamic pressure term and it is calculated as follows (Egli 2008, Holub and: where i,j, cell indices; h i,j , flow depth [m]; U i,j , flow velocity [m/s]; Q i,j , pressure load; c d , drag coefficient; q, density of the water. Dividing both terms on the right hand side of Eq. 8a by qg and assuming c d & 1, we obtain the specific load: The flow velocity is calculated from the velocities in direction x and y: where u xi;j , flow velocity in x direction (m/s); u yi;j , flow velocity in y direction (m/s). By a qualitative analysis of documented flood events (Mazzorana et al. 2009), the following assessment procedure for woody material recruitment could be established (see Fig. 4): (1) Identification and spatial delimitation of possible woody material recruitment areas, (2) assignment of these recruitment areas to the appropriate river morphology categories (e.g., stream bed, river banks, flood plain), (3) determination of the wood structure characteristics and definition of the respective structure typologies (see Table 1), (4) calculation of an indicator describing the hydraulic impact (Eq. 8a, 8b), (5) estimation of the recruitment volumes (recruited wood volume per hectare/stand volume per hectare), here referred to as woody material recruitment indicator (SVI) on the basis of the scheme shown in Fig. 4 and (6) quantification (see Eq. 10) of the absolute volume of recruited woody material (V SH ). The parameter of the absolute volume of recruited woody material indicates the maximum volume of woody material that could be ripped at a given location by the flood process.
Finally, the recruited woody material volume for each cell is calculated as Since the wood stand volume is given as a parameter showing values per hectares, Eq. 11 is used to calculate the woody material volume.
where V cell,i,j , wood stand volume for the cell i,j; V ha , wood stand volume referred to an area of 1 ha; A i,j , area of the cell i,j in m 2 . The resulting output from this step is the calculated maximal volume of woody material that could be ripped out from each cell by the flood.
The estimations of the recruited wood volumes per hectare and the quantification of the absolute volume of recruited woody material (step 5 and 6 of the earlier outlined procedure, respectively) can be computed also following a slightly modified procedure, which explicitly accounts for the subdivision of the flood duration in time steps Dt. Two essential requirements are: a. Knowledge about the dead wood material available in each cell at time t 0 . b. Knowledge about the amount of greenwood, which, through either stem breakage or ''uprooting'' (and toppling) due to erosion, becomes dead wood ready for transport at the beginning of the successive time step.
Consequently, maximum estimated dead wood amounts for different forest typologies and predictions of the dead wood amounts produced in a given time step are required. Plausible values of dead wood volumes for different wood typologies range from 3 to 40.2 m 3 /ha. In Table 1, estimations of possible dead wood volumes for different wood typologies are reported.
The prediction of the dead wood volumes, which are produced in a given time step in a determined forest typology exposed to determined hydrodynamic impacts and to determined river bed dynamics, is conduced with the following linear relationship: where V DW,i,j

Dt
, dead wood produced within the time step Dt in the cell i,j; SVI i,j t , response class [1/s] (determined in full analogy to the first theory) for cell i,j; V GW,i,j t , living wood volume in the cell at time step t [m 3 ] in cell i,j; Dt, time step (s); k, empirical constant, to be assessed on the basis observed extreme floods.
The amount of living wood at time step t ? 1 is computed as follows: It should be noted that the influence of morphodynamics is not considered in its full complexity in this approach (e.g., no mobile bed or shear stress computations are performed). Morphodynamics are considered indirectly and in a simplified way in the estimation procedure of the woody material recruitment indicator (SVI i,j t ) by assigning a location attribute (i.e., stream bed, stream bank or flood plain) to each cell (compare Figs. 2, 4).

Woody material transport dynamics
Having defined and quantified the recruitment areas, the transport of material from these areas to the defined weak points or critical river locations and cross sections can be modelled. The modelling procedure is able to consider two different approaches for modelling the transport of woody material. If the main objective is the estimation of the pathways and the possible deposition zones, the transport of woody material is calculated on a cell-by-cell basis. This approach allows analysing the transportation and deposition dynamics either on a regional level or on a more detailed level. With this modelling approach, only one time step of the hydraulic modelling results can be considered. If the main objective is to study the interactions of transported woody material with obstacles such as bridges, the transport of woody material is calculated following an object-oriented approach. This approach allows for a consideration of more time steps of the hydraulic modelling and the deposition of woody material on sand banks during the falling limb of a flood hydrograph. In such a way, it is possible to keep track of the positions of the woody material elements from time step to time step.
The basics for the calculation of the transport dynamics are the same in the two approaches and are based on the following simplified model. This model delineates the possible pathways for woody material transport and computes for each of them the entrainment and transport conditions based on the theory outlined in Sect. 2.2 and on the following method. The input data are the raster results from a hydrodynamic 2D simulation of the design event (1 in 300 years return period) for water depths h i,j and flow velocities in x and y direction, u xi;j and u yi;j . These rasters represent the state of the hydraulic simulation of different time steps. For each cell in the affected flood area, the flow direction is calculated according to the respective flow velocity in x and y directions, and given this flow direction, the source cell is moved until it is fitted into the next two neighbouring cells. The amount of woody material arriving from the source cell is divided into two portions according to the overlapping area of the moved source cell with the two neighbouring cells (Fig. 5). The particular portion of each subsequent cell is added, giving the total amount in this cell. The volumes passed through each cell are cumulated, and a grid data set of the woody material volume passed through each cell of the river influence zone is produced. This approach is based on a further simplification of Eqs. 1 and 2 in order to precisely determine the woody material transport analysis conducted on a hazard index level (Petraschek and Kienholz 2003). Assuming that the woody material elements are positioned perpendicular to the flow direction (h = 90°) and that the length of each element is expressed as a multiple of its diameter, l = nd, with n C 1; Eq. 14 is obtained from Eq. 1.
Taking the hypothesis that the density of the wood element can be expressed as r = n Á q with n % 1 and that the local slope angle at the wood element locations is very small, a % 0°, Eq. 2 can be rewritten as: Based on Eqs. 3 and 4 can be rewritten as: Subsequently, the velocity corresponding to W ¼ F R ¼ 1, here named as threshold velocity, U lim , for the movement of the wood element, is determined accordingly by Given these results and the conditions stated for expressions 5, 6 and 7, the cell-based transport inhibition parameter is defined as follows: Case 1: If h i,j C d, the wood material element is floating and the associated specific transport inhibition parameter is c i,j * = 0. Case 2: If h i,j \ d and ' s or 0\U i;j U limi;j a condition of resting is imposed to the wood material element. The associated specific transport inhibition parameters c i,j * = 1. Case 3: If h i,j \ d and U i;j [ U lim i;j a condition of either sliding or rolling is imposed to the wood material element. The value of the associated transport inhibition parameter is expressed by: where c i,j * , transport inhibition parameter of the cell i,j (non-dimensional).
Observations noted by Diehl (1997) and Ng and Richardson (2001) indicate that woody material in ideal conditions is transported on the surface as individual pieces aligned with the flow and travelling at about the same velocity as the average water velocity at the surface.
Using the average velocity instead of the surface velocity as reference velocity for the moving woody material for a wide range of flow conditions, velocity along the transport trajectory for each moving woody material model log is estimated as follows: where U wood i;j , velocity of a wood log in the cell i,j. An analysis of the expression for the wood log velocity (Eq. 19) reveals that if the transport inhibition parameter tends to 1, coherently with resting or deposition conditions the velocity of the wood log tends to 0, whereas if the transport inhibition parameter tends to 0, the velocity of the wood log tends to U i,j . Through this method, it is possible to describe woody material transport pathways under unsteady flow conditions. A deposited woody material log can be remobilised in a successive time step under changed flow depth and flow velocity conditions.

Potential hazard impacts at critical stream configurations
The object-oriented approach for modelling the woody material transport dynamics as outlined in chapter 3.3 allows for an assessment of potential hazard impacts at critical stream configurations. This approach considers the model wood logs as points, representing the centre points of the model wood log. Each single object has information about log diameter, log length, diameter of root wads and impact-resistance characteristics as outlined in Fig. 2. The flowing of the model wood logs is computed following the procedure as proposed in Chap. 3.3.
A simplified assessment procedure for entrapment and deposition phenomena at special obstacles (e.g., bridges) is outlined. In Fig. 6a, a stream section with a crossing bridge is shown. Along the flow path, woody material can potentially interact with: (1) the superstructure of the bridge, (2) a single bridge pier and/or (3) two or more bridge piers.
Two types of obstacles are defined in order to model the interaction between these obstacles and the floating woody material: In-stream-obstacles Obstacles standing in the water such as piers or abutments. On these obstacles, transported logs can get entrapped upon collision at any stage of the flood. Crossing obstacles Obstacles crossing the stream at a given height on which the logs can get entrapped with their root plate when the flow depth approaches the object's height i.e., lower chord of the superstructure of a bridge, see Fig. 6a.
Within the computational procedure, these obstacles are represented by polygonal objects. The following attributes are required to comprehensively describe these obstacles: 1. Retention probability (p ret ): Probability for each colliding log to get entrapped at the considered obstacle. This gives the expert the possibility to consider the geometry (e.g., profile) of an obstacle (e.g., inappropriately shaped piers) and estimate a proper retention probability 2. Obstacle height (h obst ): Height of the lower chord of a crossing obstacle above the initial water level. The obstacle type is defined using the obstacle height (Fig. 6a). In case the obstacle height equals zero, an object is treated as in-stream-obstacle, otherwise the object is handled as a crossing obstacle A collision with an obstacle occurs when the flow path of a woody material log intersects an obstacle. If it is a crossing obstacle, the flow depth needs to reach a critical value. Above this value, the root plate can contact the lower chord of the obstacle. Depending on the retention probability of the obstacle, specified by expert judgement, the woody material logs may either be entrapped or flushed through (see Fig. 6b).
If an object is floating between two in-stream-obstacles, e.g., two piers, it is checked whether the length of the log exceeds the shortest distance (d obst ) between the two piers. If so, a spanning blockage can occur, provided that the log is unfavourably oriented. The corresponding probability (p ent ) is estimated through expert judgement.
The entrapped wood logs occlude part of the available flow section and become a part of the obstacle for wood logs subsequently approaching the critical configuration (see Fig. 6c).

(b)
(c) (a) Fig. 6 a Example of a bridge as an obstacle consisting of two piers as in-stream-obstacles and the superstructure as a crossing obstacle d obst indicates the minimum distance between the two piers and h obst the height of the superstructure. b Example of a collision with an in-stream-object. The flow path intersects the shape of a pier. In this case, the floating log is entrapped at the first intersection point between the flow path and the shape of the pier. c Example of a woody material object floating between two piers of a bridge. The flow path does not cross any of the piers; however, the length (l log ) exceeds the minimum distance between the two pillars (d obst ) and the possibility of entrapment occurs

Study site and test application
The Passer/Passirio River in the Autonomous Province of Bolzano-South Tyrol, Italy was chosen as test site for the application of the procedure described in the previous sections (Fig. 7). The Passer/Passirio River drains a catchment area of approximately 415 km 2 and opens to the receiving watercourse Etsch/Adige River near the city of Meran/Merano, Italy.
The study area comprises the river bed and the extent of the simulated flooded areas of the Passer/Passirio River during a flood event with a return period of 1 in 300 years. A relevant volume of woody material arriving from the upper parts of the catchment is entrapped at the open check dam located in the upstream river reach. The upstream boundary of the simulation area is located at this open check dam near the community of St. Leonhard in Passeier/S. Leonardo in Passiria. A high-magnitude flash flood occurred in the Passer/ Passirio River in 1987 (Fig. 8), resulting in severe damage to regionally important bridges and roads. The downstream boundaries were defined near the locality of Saltaus/Saltusio at a bridge location (Fig. 7b). The cumulative volume of the transported woody material was calculated here. Within the studied river reach, the channel bed is characterised by twelve tributaries with relevant input of woody material. The material supplied by the tributaries was considered in this study. Figure 7c summarises the extent of the system. The woody material volumes delivered by the tributaries were assessed and quantified by analysing recent debris flow events. Table 2 shows the assumed values for the available woody material inputs. The interceptors (e.g., bridges and check dams) within the river influence zones were geo-referenced. For each of these critical configurations, the geometrical characteristics were assessed and a retention probability was assumed. This retention factor accounted for the woody material volumes retained by the obstacle. The simulation procedure was repeated with and without consideration of these weak points. The inflow hydrograph (return period of 1 in 300 years) at the upstream boundaries of the study area was calculated with the GIS-based hydrologic modelling system BaSIn 30 (AIDI 2005). The flood propagation computations were carried out with the hydrodynamic simulation model SOBEK Rural (WL/Delft Hydraulics 2004), which is capable of computing the full numerical solution of the shallow water equations. For a hazard indication analysis level, a pure 2D overland flow simulation was performed. Outputs included the flow depths and the flow velocities in x and y direction for the different time steps (e.g., 30 min for computations at the hazard index level). The flood simulation was performed on the basis of a digital elevation model delineated by airborne laser scanning technique with an original resolution of 2.5 m, upscaled to a resolution of 10 m.
For the calculation of the woody material transport dynamics, a reference diameter of the wood elements of d = 0.3 m was assumed. A drag coefficient of C d = 0.8 and a friction coefficient of l = 1.0 were used during the sets of calculation. For an analysis conducted at a hazard index level, a value of 1.4 9 10 -4 was chosen for the parameter k 1 in Eq. 12. A detailed analysis of transport dynamics at critical bridge locations has been carried out.

Simulation results
The application of the method outlined in Sect. 3.1 resulted in a map of the vegetation structures, a map of the morphological characteristics of the river influence zone and computed maximal volumes of recruited woody material within the considered system. Figure 9a shows an extract of the mapped vegetation structures within a channel section, indicating that a substantial part of the river influence zone in the study area is covered by vegetation. Thus, the recruitment of considerable amount of woody material is plausible. In Fig. 9b, the mapped geomorphologic classification of the river influence zone is shown, and in Fig. 10a, detail of the cell-by-cell simulation results is provided. The results show an increased concentration of woody material transport in the centre of the streamline. Additionally, the potential deposition areas of woody material are given; deposition primarily took place in flooded areas with low flow depths or low flow velocities outside of the river channel. Within the river channel, deposition of woody material was modelled at the waterside slopes, at the waterside of river bends and at calm loops. The modelled woody material deposition areas corresponded accurately with the mapped potential deposition areas. The maximum of transported woody material was calculated as 632 m 3 in   Table 1. This underlines that woody material transport and related phenomena are not negligible in the elaboration of the flood hazard zone map for the Passer/Passirio river. Additionally, the modelling results show the increase of transported driftwood along the river. At the upper reach of the river, small amounts of driftwood were computed; whereas downstream of each flooded wood stand with large dead wood volumes, the volume of transported driftwood increased. As expected, after a considerable deposition of driftwood, the volume of transported driftwood within the river channel decreased. Therefore, the modelling results lead to the identification of river reaches with significant driftwood transport. Since time dependency was not considered, the modelled design event had virtually an infinite duration, meaning that the entire amount of potentially removable wood stand volume was removed and mobilised. Thus, the computed woody material volumes represent the upper threshold in terms of potentially maximum values. The cell-by-cell-based modelling procedure does not consider the rising and falling limb of the flood hydrograph as does the object-oriented modelling approach. Therefore, the results of the two approaches differ slightly. In general, the modelled deposition areas of both approaches are the same, but the computed deposition volumes differ slightly. Since the object-oriented approach considers different time steps of the flood process, it considers different process areas during the flood event. Therefore, the deposited volumes differ in situations where dead wood is re-mobilised after deposition due to the increase in the wetted perimeter or the increase in flow depth and/or velocities. The object-oriented approach outlines the track of each single model wood log and as such, the origin of the entrapped or deposited log could be assessed. Additionally, the consideration of more time steps of the flood simulation enables the tracking of the development of dead wood and the transport of the model wood logs in time. It was shown that most of the model wood logs flow repeatedly out from some of the flood forest areas. Only the flooding of a large area with the same wood stand characteristics and with high flood intensities led to a simultaneous burst of woody material into the river reach. The modelling procedure allows for the observation of these situations and the assessment of the consequences of the evolvement of log jams due to a simultaneous over-flooding. Figure 10b shows the woody material deposition areas computed by the object-oriented approach. Each red dot represents one model wood log. The results of the investigation of woody material dynamics at the selected weak point (i.e., bridge location) are shown in Fig. 11. The computed interaction process of the model wood logs showed that the log jams increased constantly at obstacles during the modelled flood event. After the modelled flood process, the net volume of the log jam at the selected bridge between Oberpsairer and Mörre was around 3 m 3 . Figure 11 shows the computed log jam at the bridge, and Fig. 12 show the bridge destroyed by the flood event in 1987. The picture does not allow reconstructing the total amount of driftwood entrapped at the bridge, as the bridge was over-flooded after the break of the left pillar and the entrapped woody material was transported away. However, Fig. 11 confirms the entrapment of driftwood with this obstacle. The procedure for modelling the interaction of the model wood logs with critical stream configurations calculates the volume of retained driftwood at every obstacle. The retention volume depends on the characteristics of the obstacle itself, and trapped driftwood necessarily reduces the volume which is transported downstream of the obstacle. This method allows for the identification of the systematically most relevant weak points in the system.

Discussion and conclusions
By the application of the developed method, valuable insights were provided regarding woody material recruitment processes and the propensity for entrainment and transport to critical configurations during extreme flood events. In particular, an overview of the Fig. 11 Entrapment at a critical configuration. The graphic at the right side shows the entrapped driftwood on a more detailed scale Fig. 12 Bridge destroyed across the Passer/Passirio river soon after the flood event in July 18th, 1987. The bridge was occluded with woody material and was destroyed after the outbreak of the river on the left side. The figure shows only the driftwood at the bridge that remained after the event maximal amount of transported and deposited woody material in a defined river reach was given. Moreover, the main pathways of the woody material in the river channel were traced and subsequently the potential depositional areas were identified. Furthermore, the procedure provided the order of dimension of the woody material volume expected to pass through bridges during a flood event. Such quantified knowledge on woody material volume approaching a weak point location is crucial for a reliable scenario definition during hazard assessment, in particular with respect to a possible log jam formation, clogging or similar obstruction phenomena. Consequently, the transparency of hazard mapping procedure and the quality of the results are increased, and the supervision of the entire hazard mapping procedure by the respective public authorities is facilitated. The outcomes can also be used in order to support the definition of policies in riparian forest management and for defining particular measures such as thinning and other forest management actions. Nevertheless, the procedure presented earlier might still have certain limitations. The spatial resolution of the flood simulation has to be adapted and refined if the goal is a more detailed scale. With a grid size of 10 m, the river channel morphology is approximated with limited accuracy, resulting in a less precise calculation of the woody material flow paths. The procedure was found to be sensitive to the assumed wood stand volumes of different vegetation structures studied; however, a very detailed investigation on the wood structure might improve the procedure significantly. Additionally, the exposure time of the different vegetation structures to the flood event might be another relevant factor. During flood events with a relatively short duration, significantly smaller volumes are expected compared to long-duration flood events.
Despite these limitations, the developed method is of particular relevance to mitigate flood risk. The results of mapped recruitment areas with the respective forest structure typologies and identified transport paths with the transport dynamics for a determined critical configuration are essential indicators for hazard assessment.
The time-dependent modelling of woody material dynamics indicates whether or not the formation of log jams due to a simultaneous flooding of a large forested area with high flood intensities is plausible in the studied river reach. The assessment of the interaction processes induced at the considered critical configuration refines the hazard analysis and provides significantly more detailed input for the subsequent risk assessment. As shown during previous events, the critical configuration was repeatedly located at bridges; hence, the hypothetical debris accumulation for the entire bridge and the resulting consequences can be calculated according to the indications provided by Diehl (1997). This is fundamental information for effective risk mitigation strategies. Key elements of these strategies include: (1) the removal of those parts of the critical configuration that induce woody material accumulation, (2) the reconfiguration of the weak point to improve flow conditions and (3) the adherence to silvicultural measures in the respective recruitment areas in order to reduce the hazard source.
An application of the proposed model would allow enhanced emergency planning and preparations (e.g., excavating identified places in order to diminish accumulations at bridges; installation of temporary protection measures; implementation of local structural protection). These activities would introduce redundancies and buffer capacities into the system, thereby achieving an increased resilience of elements at risk. anonymous reviewers and to Sven Fuchs and Francesco Comiti for valuable comments on an earlier version of the manuscript, to Emily Procter for improving the English manuscript and to Claudio Volcan for performing the hydrodynamic simulations.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.