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.     

Adiabatic chaos in the spin–orbit problem

We provide evidences that the angular momentum of a symmetric rigid body in a spin–orbit resonance can perform large scale chaotic motions on time scales which increase polynomially with the inverse of the oblateness of the body. This kind of irregular precession appears as soon as the orbit of the center of mass is non-circular and the angular momentum of the body is far from the principal directions with minimum (maximum) moment of inertia. We also provide a quantitative explanation of these facts by using the theory of adiabatic invariants, and we provide numerical applications to the cases of the 1:1 and 1:2 spin–orbit resonances.     

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.     

鄂尔多斯盆地西南部三叠纪末抬升剥蚀事件及热年代学记录

三叠纪末期大型鄂尔多斯盆地遭受了中生代成盆以来首次较大规模的抬升剥蚀,显著改造了中晚三叠世延长期盆地面貌,并控制了侏罗纪早期沉积格局和油藏分布,对盆地演化及矿产资源分布产生了重要影响。本文利用地质及大量钻井资料揭示了该期构造事件对盆地的剥蚀改造特征,盆地及周邻地区磷灰石裂变径迹年代学记录并约束了此次构造抬升的时限与过程;综合周邻区域构造研究成果,探讨了其发育的动力学背景。结果表明,鄂尔多斯盆地三叠纪末期的剥蚀具西南强、东北弱的特点,西南部大范围内延长组地层残缺不全,剥蚀量最大可达1000余米;前侏罗纪沉积古地貌总体呈西南高、北东低的特点;其抬升时间始于205～190Ma,西南部稍早于盆地腹部,抬升速率大于1℃/Myr,可持续至中侏罗世(约160Ma)。该期抬升剥蚀事件范围可涉及至盆地西南缘更广阔的区域,与同期秦岭造山带内出现的快速抬升冷却事件具有较好的时空耦合关系,是对秦岭造山带区域构造环境转变的响应和纪录。该研究丰富和发展了三叠纪末期构造事件在华北克拉通的影响,对该区油气、煤炭资源的进一步勘探和评价提供了新的思路,具有一定的现实意义。     

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.     

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.     

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.