Finite element analysis of land subsidence above depleted reservoirs with pore pressure gradient and total stress formulations

The solution of the poroelastic equations for predicting land subsidence above productive gas/oil fields may be addressed by the principle of virtual works using either the effective intergranular stress, with the pore pressure gradient regarded as a distributed body force, or the total stress incorporating the pore pressure. In the finite element (FE) method both approaches prove equivalent at the global assembled level. However, at the element level apparently the equivalence does not hold, and the strength source related to the pore pressure seems to generate different local forces on the element nodes. The two formulations are briefly reviewed and discussed for triangular and tetrahedral finite elements. They are shown to yield different results at the global level as well in a three‐dimensional axisymmetric porous medium if the FE integration is performed using the average element‐wise radius. A modification to both formulations is suggested which allows to correctly solve the problem of a finite reservoir with an infinite pressure gradient, i.e. with a pore pressure discontinuity on its boundary. Copyright © 2001 John Wiley & Sons, Ltd.     

Two-dimensional dam break flooding simulation: a GIS-embedded approach

In the twenty-first century, around 200 notable dam and reservoir failures happened worldwide causing massive fatalities and economic costs. In order to reduce the losses, managers usually define mitigation strategies identifying flooding area due to dam break by using standalone hydrodynamic models and then importing the results within a GIS to perform risk analysis. This two-step procedure is time expensive, error prone due to export/import requirements and not user friendly. For this reason with this work, a new numerical model for the solution of the two-dimensional dam break problem has been implemented in the GRASS GIS with a GIS-embedded approach. The model solves the conservative form of the 2D shallow water equations using a finite volume method; the intercell flux is computed by one-side upwind conservative scheme extended to a two-dimensional problem. The newly developed GIS module, among others outputs, allows to derive maximum intensity maps that can be directly used for risk assessment. Finally, the model has been (1) tested against two standard synthetic problems referenced in literature showing differences in estimated water depth of 2, 3 and 15% and (2) verified against official flooding map of an existing dam (Verzasca) detecting 75% of similarity. The problem formulation, the new GRASS module and its validation is presented.     

Macroseismic effects highlight site response in Rome and its geological signature

A detailed analysis of the earthquake effects on the urban area of Rome has been conducted for the L??Aquila sequence, which occurred in April 2009, by using an online macroseismic questionnaire. Intensity residuals calculated using the mainshock and four aftershocks are analyzed in light of a very accurate and original geological reconstruction of the subsoil of Rome based on a large amount of wells. The aim of this work is to highlight ground motion amplification areas and to find a correlation with the geological settings at a subregional scale, putting in evidence the extreme complexity of the phenomenon and the difficulty of making a simplified model. Correlations between amplification areas and both near-surface and deep geology were found. Moreover, the detailed scale of investigation has permitted us to find a correlation between seismic amplification in recent alluvial settings and subsiding zones, and between heard seismic sound and Tiber alluvial sediments.     

IFKIS-Hydro: an early warning and information system for floods and debris flows

IFKIS-Hydro is an information and warning system for hydrological hazards in small- and medium-scale catchments. The system collects data such as weather forecasts, precipitation measurements, water level gauges, discharge simulations and local observations of event-specific phenomena. In addition, IFKIS-Hydro incorporates a web-based information platform, which serves as a central hub for the submission and overview of data. Special emphasis is given to local information. This is accomplished particularly by human observers. In medium-scale catchments, discharge forecast models have an increasing importance in providing valuable information. IFKIS-Hydro was developed in several test regions in Switzerland and the first results of its application are available now. The system is constantly extended to additional regions and may become the standard for warning systems in smaller catchments in Switzerland.     

Toward a rapid probabilistic seismic vulnerability assessment using satellite and ground-based remote sensing

Natural hazards such as earthquakes threaten millions of people all around the world. In a few decades, most of these people will live in fast-growing, inter-connected urban environments. Assessing risk will, therefore, be an increasingly difficult task that will require new, multidisciplinary approaches to be tackled properly. We propose a novel approach based on different imaging technologies and a Bayesian information integration scheme to characterize exposure and vulnerability models, which are among the key components of seismic risk assessment.     

Lava flow hazard at Nyiragongo volcano,D.R.C.

The 2002 eruption of Nyiragongo volcano constitutes the most outstanding case ever of lava flow in a big town. It also represents one of the very rare cases of direct casualties from lava flows, which had high velocities of up to tens of kilometer per hour. As in the 1977 eruption, which is the only other eccentric eruption of the volcano in more than 100 years, lava flows were emitted from several vents along a N–S system of fractures extending for more than 10 km, from which they propagated mostly towards Lake Kivu and Goma, a town of about 500,000 inhabitants. We assessed the lava flow hazard on the entire volcano and in the towns of Goma (D.R.C.) and Gisenyi (Rwanda) through numerical simulations of probable lava flow paths. Lava flow paths are computed based on the steepest descent principle, modified by stochastically perturbing the topography to take into account the capability of lava flows to override topographic obstacles, fill topographic depressions, and spread over the topography. Code calibration and the definition of the expected lava flow length and vent opening probability distributions were done based on the 1977 and 2002 eruptions. The final lava flow hazard map shows that the eastern sector of Goma devastated in 2002 represents the area of highest hazard on the flanks of the volcano. The second highest hazard sector in Goma is the area of propagation of the western lava flow in 2002. The town of Gisenyi is subject to moderate to high hazard due to its proximity to the alignment of fractures active in 1977 and 2002. In a companion paper (Chirico et al., Bull Volcanol, in this issue, 2008) we use numerical simulations to investigate the possibility of reducing lava flow hazard through the construction of protective barriers, and formulate a proposal for the future development of the town of Goma.     

The thermal evolution of Corsica as recorded by zircon fission-tracks

New zircon fission-track (ZFT) ages from Corsica record multiple thermal events that can be tied to the structural evolution of the western Mediterranean region. The Corsican zircons have a wide scatter of ZFT grain ages (243–14 Ma), which together define several age domains. Western Corsica consists largely of stable Hercynian basement characterized by ZFT ages in the range 161–114 Ma. We interpret these ages (Late Jurassic–Early Cretaceous) as the product of a long-lived Tethyan thermal event related to continental rifting and subsequent drifting during the separation of the European and African plates and the formation of the Liguro–Piemontese ocean basin. In contrast to Hercynian Corsica, Alpine Corsica (northeast Corsica) experienced widespread deformation and metamorphism in Late Cretaceous(?)–Tertiary time. Dated samples from Alpine Corsica range in age from 112 to 19 Ma and all are reset or partially reset by one or more Alpine thermal events. The youngest ZFT grain ages are from the northernmost Alpine Corsica and define an age population at  24 Ma that indicates cooling after Tertiary thermal events associated with the Alpine metamorphism and the opening of the Liguro–Provençal basin. A less well-defined ZFT age population at  72 Ma is present in both Alpine Corsica and Hercynian basement rocks. The thermal history of these rocks is not clear. One interpretation is that the ZFT population at  72 Ma reflects resetting during a Late Cretaceous event broadly synchronous with the early Alpine metamorphism. Another interpretation is that this peak is related to variable fission-track annealing and partial resetting during the Tertiary Alpine metamorphic event across central to north-eastern Corsica. This partial age resetting supports the presence of a fossil ZFT partial annealing zone and limits the peak temperature in this area below 300 °C, for both the affected pre-Alpine and Alpine units.     

Markovian modeling of seismic damage accumulation

Analysis of civil structures at the scale of life‐cycle requires stochastic modeling of degradation. Phenomena causing structures to degrade are typically categorized as aging and point‐in‐time overloads. Earthquake effects are the members of the latter category this study deals with in the framework of performance‐based earthquake engineering (PBEE). The focus is structural seismic reliability, which requires modeling of the stochastic process describing damage progression, because of subsequent events, over time. The presented study explicitly addresses this issue via a Markov‐chain‐based approach, which is able to account for the change in seismic response of damaged structures (i.e. state‐dependent seismic fragility) as well as uncertainty in occurrence and intensity of earthquakes (i.e. seismic hazard). The state‐dependent vulnerability issue arises when the seismic hysteretic response is evolutionary and/or when the damage measure employed is such that the degradation increment probabilistically depends on the conditions of the structure at the time of the shock. The framework set up takes advantage also of the hypotheses of classical probabilistic seismic hazard analysis, allowing to separate the modeling of the process of occurrence of seismic shocks and the effect they produce on the structure. It is also discussed how the reliability assessment, which is in closed‐form, may be virtually extended to describe a generic age‐ and state‐dependent degradation process (e.g. including aging and/or when aftershock risk is of interest). Illustrative applications show the options to calibrate the model and its potential in the context of PBEE. Copyright © 2015 John Wiley & Sons, Ltd.     

Exhumation patterns along shallow low‐angle normal faults: an example from the Altotiberina active fault system (Northern Apennines,Italy)

A multi‐method approach (palaeothermal and thermochronological analyses; thermal modelling) is applied to reconstruct the exhumation history of the Altotiberina Fault (ATF), a representative example of crustal‐scale active low‐angle normal faulting in the Northern Apennines (Italy). Thermal maturity and thermochronological data yield similar burial histories but different exhumation patterns for the sedimentary successions in the hangingwall and the footwall of the ATF. Since 3.8 Ma, the ATF footwall has exhumed at rates of 0.90 mm a^?1. Exhumation led to bending and deactivation of the ATF uppermost portion as a result of tectonic unloading and isostatic adjustment, followed by migration of extension and the development of a set of domino‐like, east‐dipping normal faults, rooting on the buried portion of the ATF. ATF activity and isostatic rebound exhumed Triassic rock units from depths of about 4 km. We suggest that isostatic instability is accommodated at shallow crustal levels, in a similar way to what is observed on larger structures at mid‐low crustal levels.