Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

False waves in wave records

A methodology to define discrete waves from free sea surface elevation time series is presented. The method allows an objective discrimination of false waves among small waves, avoiding the use of arbitrary criteria associated with zero-up-crossing, zero-down-crossing, or other definitions. The method uses the Hilbert Transform and the representation of the time series in the complex plane. A discrete wave corresponds to a 2π phase-advance in the complex plane. The waves between zero crossings which do not show 2π phase advance are considered false waves. Wave rider buoy records, measured off the west coast of Portugal, were employed in assessing some of the statistical implications of this methodology.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

Deformability characteristics of brazilian laterites

This paper summarizes the resilient and permanent deformation characteristics of laterite gravels, determined by cyclic loading in triaxial tests of compacted specimens. A study of laterite pavements was carried out on eight paved road sections in the state of Mato Grosso. Also reported is the influence of fine sand admixture, specimen size and soaking in two laterite samples from Roraima and Rondônia in the North of Brazil. The experimental relationships of resilient moduli against applied stresses are interpreted empirically for models what may be incorporated in a structural analysis of full depth of laterite pavement.     

Surface-water quality for the Salí River watershed in NW Argentina

Spatial and temporal variations of Na, Mg, Si, K, Ca, SO₄^2-, Cl^-, HCO₃^-, Ag, Al, As, Au, B, Ba, Be, Bi, Br, Ce, Cd, Co, Cr, Cs, Cu, Dy, Er, Eu, F, Fe, Ga, Ge, Gd, Hf, Hg, Ho, I, La, Li, Mn, Mo, Nb, Nd, Ni, P, Pb, Pr, Rb, Sb, Sm, Sn, Sr, Ta, Tb, Te, Th, Ti, Tl, U, V, W, Y, Yb, Zn, and Zr were monitored through 37 sampling stations to determine the main aspects influencing the surface-water quality in the Salí River watershed (Tucumán Province, NW Argentina). The influence of the regional geological setting on water chemistry allows to distinguish three sub-basins. The interaction with sedimentary rocks was found to be dominant in the northern and central-eastern sub-basin as well as in the southern sub-basin, whereas the metamorphic-granitic basement of the Pampean Ranges was noted in the central sub-basin. In addition, anthropogenic activities affect the spatial variation of K, P, Mn, Rb, and Pb as well as dissolved oxygen concentrations and Eh. Temporal water-quality variation is related to the spatial distribution of precipitation and to the seasonal character of the main local industries (sugar cane, alcohol, citrus), increasing P and K concentrations and decreasing dissolved oxygen concentration and Eh in winter. Cl^-, Na, SO₄^2-, Al, As, B, Fe, Mn, Se and U concentrations exceed the regulated drinking-water thresholds at several sampling stations.     

Hydrogeology and sustainable future groundwater abstraction from the Agua Verde aquifer in the Atacama Desert,northern Chile

The hyper-arid conditions prevailing in Agua Verde aquifer in northern Chile make this system the most important water source for nearby towns and mining industries. Due to the growing demand for water in this region, recharge is investigated along with the impact of intense pumping activity in this aquifer. A conceptual model of the hydrogeological system is developed and implemented into a two-dimensional groundwater-flow numerical model. To assess the impact of climate change and groundwater extraction, several scenarios are simulated considering variations in both aquifer recharge and withdrawals. The estimated average groundwater lateral recharge from Precordillera (pre-mountain range) is about 4,482 m³/day. The scenarios that consider an increase of water withdrawal show a non-sustainable groundwater consumption leading to an over-exploitation of the resource, because the outflows surpasses inflows, causing storage depletion. The greater the depletion, the larger the impact of recharge reduction caused by the considered future climate change. This result indicates that the combined effects of such factors may have a severe impact on groundwater availability as found in other groundwater-dependent regions located in arid environments. Furthermore, the scenarios that consider a reduction of the extraction flow rate show that it may be possible to partially alleviate the damage already caused to the aquifer by the continuous extractions since 1974, and it can partially counteract climate change impacts on future groundwater availability caused by a decrease in precipitation (and so in recharge), if the desalination plant in Taltal increases its capacity.     

Comparison of long-, medium- and short-term variations of beach profiles with and without submerged geological control

Victoria Beach (Cadiz, Spain) comprises a rocky flat outcrop in its northern zone and a sand-rich southern zone. These natural features allowed for a 5-year monitoring period and subsequent analysis of two different profiles (one in each zone) based on differences in bottom contours. Topo-bathymetric data were analysed using empirical orthogonal functions (EOFs) to determine changes over the short-, medium- and long-term. Several morphologic phenomena were identified (generalised erosion, seasonal or summer–winter tilting of the profile around different hinge points, berm development and its posterior destruction, etc.) in terms of their importance in explaining the variability of the collected data for both profiles. It is worth mentioning that both profiles undergo parallel regression in the medium-term. Thus, the 1st eigenfunction enabled us to identify the true regression of the beach shoreline, independent of seasonal or summer–winter slope changes. Reconstruction of profiles using EOF components demonstrated that though accretion periods in the medium-term were similar for both types of profiles, the accretion speed was much faster in the sand-rich profile than in the reef-protected profile (1.01 m³/day versus 0.33 m³/day). Moreover, the seasonal erosion rate and the subsequent shoreline retreat for the sand-rich profile were much larger than for the reef-protected profile (121 m³/year versus 29 m³/year). Analysis in the short-term (changes induced by a single day's storm) showed an instantaneous tilting of the profile, with the mobilised sand volume being much greater for the sand-rich than for the reef-protected profile (68 m³/m versus 12 m³/m).     

Effects of wave–current interaction on the current profile

Wave–current laboratory experiments have shown that the logarithmic current profile observed in pure current flows is modified due to the presence of waves. When waves propagate opposite the current, an increase in the current intensity is achieved near the mean water level, while a reduction is obtained for following waves and currents. With the aim of analyzing these nonlinear effects along the whole water column, an Eulerian wave–current model is presented. In contrast to previously presented wave–current models, the present is able to include the variation of the free surface elevation due to the wave motion and the effect of a non hydrostatic pressure field. Therefore it does not restrict its application to waves in shallow waters. Moreover, the model is able to simulate all the possible angles between waves and currents.     

Acceleration response modification factors for nonstructural components attached to inelastic moment‐resisting frame structures

A statistical analysis of the peak acceleration demands for nonstructural components (NSCs) supported on a variety of stiff and flexible inelastic regular moment‐resisting frame structures with periods from 0.3 to 3.0 s exposed to 40 far‐field ground motions is presented. Peak component acceleration (PCA) demands were quantified based on the floor response spectrum (FRS) method without considering dynamic interaction effects. This study evaluated the main factors that influence the amplification or decrease of FRS values caused by inelasticity in the primary structure in three distinct spectral regions namely long‐period, fundamental‐period, and short‐period region. The amplification or decrease of peak elastic acceleration demands depends on the location of the NSC in the supporting structure, periods of the component and building, damping ratio of the component, and level of inelasticity of the supporting structure. While FRS values at the initial modal periods of the supporting structure are reduced due to inelastic action in the primary structure, the region between the modal periods experiences an increase in PCA demands. A parameter denoted as acceleration response modification factor (R_acc) was proposed to quantify this reduction/increase in PCA demands. Copyright © 2007 John Wiley & Sons, Ltd.