A novel methodological approach for land subsidence prediction through data assimilation techniques

Computational Geosciences - Anthropogenic land subsidence can be evaluated and predicted by numerical models, which are often built over deterministic analyses. However, uncertainties and...     

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.     

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.     

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.     

Robust numerical implementation of a 3D rate-dependent model for reservoir geomechanical simulations

A 3D elasto-plastic rate-dependent model for rock mechanics is formulated and implemented into a Finite Element (FE) numerical code. The model is based on the approach proposed by Vermeer and Neher (A soft soil model that accounts for creep. In: Proceedings of the International Symposium “Beyond 2000 in Computational Geotechnics,” pages 249-261, 1999). An original strain-driven algorithm with an Inexact Newton iterative scheme is used to compute the state variables for a given strain increment.The model is validated against laboratory measurements, checked on a simplified test case, and used to simulate land subsidence due to groundwater and hydrocarbon production. The numerical results prove computationally effective and robust, thus allowing for the use of the model on real complex geological settings.     

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.     

Parameterisation of the Planetary Boundary Layer for Diagnostic Wind Models

The planetary boundary-layer (PBL) parameterization is a key issue for the definition of initial wind flow fields in diagnostic models. However, PBL theories usually treat separately stable, neutral, and convective stability conditions, so that their implementation in diagnostic wind models is not straightforward. In the present paper, an attempt is made to adopt a comprehensive PBL parameterisation, covering stable/neutral and unstable atmospheric conditions, which appears suitable to diagnostic models. This parameterisation is implemented into our diagnostic mass-consistent code. A validation of the consistency between the implemented PBL parameterisations has been checked through an analysis of the sensitivity of the vertical wind profiles to atmospheric stability.     

High-temperature emplacement of the Cerro Galán and Toconquis Group ignimbrites (Puna plateau,NW Argentina) determined by TRM analyses

Estimates of pyroclastic flow emplacement temperatures in the Cerro Galán ignimbrite and Toconquis Group ignimbrites were determined using thermal remanent magnetization of lithic clasts embedded within the deposits. These ignimbrites belong to the Cerro Galán volcanic system, one of the largest calderas in the world, in the Puna plateau, NW Argentina. Temperature estimates for the 2.08-Ma Cerro Galán ignimbrite are retrieved from 40 sites in 14 localities (176 measured clasts), distributed at different distances from the caldera and different stratigraphic heights. Additionally, temperature estimates were obtained from 27 sample sites (125 measured clasts) from seven ignimbrite units forming the older Toconquis Group (5.60–4.51 Ma), mainly outcropping along a type section at Rio Las Pitas, Vega Real Grande. The paleomagnetic data obtained by progressive thermal demagnetization show that the clasts of the Cerro Galán ignimbrite have one single magnetic component, oriented close to the expected geomagnetic field at the time of emplacement. Results show therefore that most of the clasts acquired a new magnetization oriented parallel to the magnetic field at the moment of the ignimbrite deposition, suggesting that the clasts were heated up to or above the highest blocking temperature (T _b) of the magnetic minerals (T _b = 580°C for magnetite; T _b = 600–630°C for hematite). We obtained similar emplacement temperature estimations for six out of the seven volcanic units belonging to the Toconquis Group, with the exception of one unit (Lower Merihuaca), where we found two distinct magnetic components. The estimation of emplacement temperatures in this latter case is constrained at 580–610°C, which are lower than the other ignimbrites. These estimations are also in agreement with the lowest pre-eruptive magma temperatures calculated for the same unit (i.e., 790°C; hornblende–plagioclase thermometer; Folkes et al. 2011b). We conclude that the Cerro Galán ignimbrite and Toconquis Group ignimbrites were emplaced at temperatures equal to or higher than 620°C, except for Lower Merihuaca unit emplaced at lower temperatures. The homogeneity of high temperatures from proximal to distal facies in the Cerro Galán ignimbrite provides constraints for the emplacement model, marked by a relatively low eruption column, low levels of turbulence, air entrainment, surface–water interaction, and a high level of topographic confinement, all ensuring minimal heat loss.     

Expected loss‐based alarm threshold set for earthquake early warning systems

Earthquake early warning systems (EEWS) seem to have potential as tools for real‐time seismic risk management and mitigation. In fact, although the evacuation of buildings requires warning time not available in many urbanized areas threatened by seismic hazard, they may still be used for the real‐time protection of critical facilities using automatic systems in order to reduce the losses subsequent to a catastrophic event. This is possible due to the real‐time seismology, which consists of methods and procedures for the rapid estimation of earthquake features, as magnitude and location, based on measurements made on the first seconds of the P‐waves. An earthquake engineering application of earthquake early warning (EEW) may be intended as a system able to issue the alarm, if some recorded parameter exceeds a given threshold, to activate risk mitigation actions before the quake strikes at a site of interest. Feasibility analysis and design of such EEWS require the assessment of the expected loss reduction due to the security action and set of the alarm threshold. In this paper a procedure to carry out these tasks in the performance‐based earthquake engineering probabilistic framework is proposed. A merely illustrative example refers to a simple structure assumed to be a classroom. Structural damage and non‐structural collapses are considered; the security action is to shelter occupants below the desks. The cost due to a false alarm is assumed to be related to the interruption of didactic activities. Results show how the comparison of the expected losses, for the alarm‐issuance and non‐issuance cases, allows setting the alarm threshold on a quantitative and consistent basis, and how it may be a tool for the design of engineering applications of EEW. Copyright © 2007 John Wiley & Sons, Ltd.