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.     

Mobilization analysis of the dolmen of Dombate (Northwest Spain)

This paper presents an analysis of the mobilization of the dolmen of Dombate, a megalithic monument situated in the Northwest of Spain. The study is founded on simple static principles, with the application of basic concepts of soil mechanics (soil reactions were computed by the reaction modulus method) and rock mechanics (a simple secant approach was adopted to model the contact between orthostats). In addition, a simplified geometric model for dolmens was adopted. On the basis of these concepts, a numerical solver able to provide quick estimations of the system's mobilization was developed. Using this tool, a number of sensitivity analyses were conducted to estimate the mobilization that could be caused by a conservation project currently under consideration. It was found that mobilization might be excessive, and therefore, a safer procedure is proposed. In short, this paper puts forth a new rational approach for archaeological practices which may help to estimate the possible consequences of excavations and rehabilitations on megalithic structures.     

Equilibrium constants of triethanolamine in major seawater salts

Acid–base equilibrium constants of triethanolamine (TEA) have been determined by potentiometric titrations with a glass electrode, at 25 °C. Ionic strength was kept constant with only one electrolyte (using one of these salts: NaCl, KCl, MgCl₂ or CaCl₂), with binary mixtures of MgCl₂ and CaCl₂, and finally, in a solution with a composition approximately similar to that of natural seawater without sulfate. Equilibrium constants have been expressed in function of ionic strength by means of Pitzer equations and interaction parameters proposed in this theory have been obtained. It has been found that acid–base behaviour of TEA depends greatly on the salt used: basicity of TEA is decreased by CaCl₂, while it is increased by the other electrolytes used in this work.     

Spatial Heterogeneity Created by Burrowing Crabs Affects Human Impact on Migratory Shorebirds

We show that, as previously described, there are conflicts between areas used by people and shorebirds. However, we propose that burrowing organisms add complexity to the shorebird–human interaction that should be taken into account for management and conservation. This is because SW Atlantic mudflats are dominated by high densities (up to 60 crabs m⁻²) of a strong bioturbator, the burrowing crab Neohelice granulata (=Chasmagnathus granulatus). These crabs affect the habitat use and foraging performance of shorebird species. The two-banded plover Charadrius falklandicus feed more in burrow areas whereas the yellowlegs Tringa spp. and the white-rumped sandpiper Calidris fuscicollis tend to avoid burrow areas. We evaluated the effect of human recreational activities on shorebirds, taking into account shorebird–crab interaction, during shorebird migratory periods (November to April) of 2001, 2002, and 2003. This study was performed at the Mar Chiquita coastal lagoon (37° 40′ S, Argentina). Results showed that people only rarely used the extensive soft bottom intertidals dominated by burrowing assemblages of N. granulata because the soft sediments make walking difficult. Therefore, human impact is reduced for the two-banded plover. However, for several migratory shorebirds that also avoid burrow areas, these areas act as suboptimal alternative areas when human impact is important in their preferred habitats. When human activity is low, shorebirds remain foraging in the area but they decrease their feeding intakes. The extensions of burrow areas differ among estuaries and are likely to generate between-estuaries differences as stopover sites for shorebirds. Also, other species that form large burrowing assemblages can have a similar impact on shorebird–human interactions.     

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.     

Sediment Organic Carbon in Todos Santos Bay, Baja California, Mexico

Sediment grain size and organic carbon (OC) data collected over the past 50 years, together with δ¹³C values of OC in recently collected samples, were analyzed to improve understanding of sediment OC distribution and abundance in Todos Santos Bay. Sediments in the submarine canyon at the mouth of the bay and in quiet-water locations along the shore are fine grained, high in OC, and have generally low δ¹³C values; sediments in high-energy environments are low in OC and have high δ¹³C values. A bivariate isotopic mixing model indicates that none of the sediments contain >50% terrigenous OC (average ~30%), and that the terrigenous OC content of the sediments is a small proportion of the OC content of local soils. Sediment OC composition is apparently controlled by energy-related sorting and deposition, oxidation of much of the original terrigenous OC, and replacement of some terrigenous OC by marine OC.     

The Geysers - Cobb Mountain Magma System, California (Part 1): U-Pb zircon ages of volcanic rocks, conditions of zircon crystallization and magma residence times

Combined U-Pb zircon and ⁴⁰Ar/³⁹Ar sanidine data from volcanic rocks within or adjacent to the Geysers geothermal reservoir constrain the timing of episodic eruption events and the pre-eruptive magma history. Zircon U-Pb concordia intercept model ages (corrected for initial ²³⁰Th disequilibrium) decrease as predicted from stratigraphic and regional geological relationships (1σ analytical error): 2.47 ± 0.04 Ma (rhyolite of Pine Mountain), 1.38 ± 0.01 Ma (rhyolite of Alder Creek), 1.33 ± 0.04 Ma (rhyodacite of Cobb Mountain), 1.27 ± 0.03 Ma (dacite of Cobb Valley), and 0.94 ± 0.01 Ma (dacite of Tyler Valley). A significant (∼0.2-0.3 Ma) difference between these ages and sanidine ⁴⁰Ar/³⁹Ar ages measured for the same samples demonstrates that zircon crystallized well before eruption. Zircons U-Pb ages from the underlying main-phase Geysers Plutonic Complex (GPC) are indistinguishable from those of the Cobb Mountain volcanics. While this is in line with compositional evidence that the GPC fed the Cobb Mountain eruptions, the volcanic units conspicuously lack older (∼1.8 Ma) zircons from the shallowest part of the GPC. Discontinuous zircon age populations and compositional relationships in the volcanic and plutonic samples are incompatible with zircon residing in a single long-lived upper crustal magma chamber. Instead we favor a model in which zircons were recycled by remelting of just-solidified rocks during episodic injection of more mafic magmas. This is consistent with thermochronologic evidence that the GPC cooled below 350° C at the time the Cobb Mountain volcanics were erupted.