Crustal melting and recycling: geochronology and sources of Variscan syn-kinematic anatectic granitoids of the Tormes Dome (Central Iberian Zone). A U–Pb LA-ICP-MS study

International Journal of Earth Sciences - In this study, we report U–Pb Laser Ablation ICP-MS zircon and ID-TIMS monazite ages for peraluminous granitoid plutons...     

A Cambrian mixed carbonate–siliciclastic platform in SW Gondwana: evidence from the Western Sierras Pampeanas (Argentina) and implications for the early Paleozoic paleogeography of the proto-Andean margin

International Journal of Earth Sciences - The Western Sierras Pampeanas (WSP) of Argentina record a protracted geological history from the Mesoproterozoic assembly of the Rodinia supercontinent to...     

Site scale modeling of slow-moving landslides,a 3D viscoplastic finite element modeling approach

This paper presents an advanced 3D numerical methodology to reproduce the kinematics of slow active landslides, more precisely, to reproduce the nearly constant strain rate (secondary creep) and the acceleration/deceleration of the moving mass due to hydrological changes. For this purpose, finite element analyses are performed in a large area covering a long time-span (12 years), in order to exhibit different interacting slope movements. First, we perform a stability analysis using the shear strength reduction (SSR) technique with a Mohr-Coulomb failure criteria. It is done in order to compute factors of safety (FS) and to identify two different scenarios, the first one being stable (FS > 1) and the second one being unstable (FS < 1). In the studied test case, the Portalet landslide (Central Spanish Pyrenees), the first scenario corresponds to an initial stable configuration of the slope and the second one to an unstable excavated configuration. Second, taking the first scenario as an initial condition, a time-dependent analysis is performed using a coupled formulation to model solid skeleton and pore fluids interaction, and a simplified ground water model that takes into account daily rainfall intensity. In this case, a viscoplastic constitutive model based on Perzyna’s theory is applied to reproduce soil viscous behavior and the delayed creep deformation due to the excavation. The fluidity parameter is calibrated to reproduce displacements measured by the monitoring systems. Our results demonstrate that 3D analyses are preferable to 2D ones for reproducing in a more realistic way the slide behavior. After calibration, the proposed model is able to simulate successfully short- and medium-term predictions during stages of primary and secondary creep.     

Rainfall thresholds for the triggering of landslides considering previous soil moisture conditions (Asturias,NW Spain)

Given its geological and climatic conditions and its rugged orography, Asturias is one of the most landslide prone areas in the North of Spain. Most of the landslides occur during intense rainfall episodes. Thus, precipitation is considered the main triggering factor in the study area, reaching average annual values of 960 mm. Two main precipitation patterns are frequent: (i) long-lasting periods of moderate rainfall during autumn and winter and (ii) heavy short rainfall episodes during spring and early summer. In the present work, soil moisture conditions in the locations of 84 landslides are analysed during two rainfall episodes, which represent the most common precipitation patterns: October–November 2008 and June 2010. Empirical data allowed the definition of available water capacity percentages of 99–100% as critical soil moisture conditions for the landslide triggering. Intensity-duration rainfall thresholds were calculated for each episode, considering the periods with sustained high soil moisture levels before the occurrence of each analysed landslide event. For this purpose, data from daily water balance models and weather stations were used. An inverse relationship between the duration of the precipitation and its intensity, consistent with published intensity-duration thresholds, was observed, showing relevant seasonal differences.     

Creep behaviour of intact and remoulded fibrous peat

This paper presents the creep behaviour of intact and remoulded specimens of fibrous peat obtained from a field site near Anzac, Alberta, Canada. The creep behaviour was investigated by means of long-term drained and undrained triaxial tests. The development of volumetric, axial, and undrained axial strain and strain rate during drained and undrained creep tests under variable stress conditions is presented. The stress _– strain _– strain rate (p′–ε _v–\(\dot{\varepsilon }_{\text{v}}\)) relationship is found to be unique for different stress and loading durations. The p′–ε _v–\(\dot{\varepsilon }_{\text{v}}\) relationship is analysed and represented by creep isotaches. The applicability of different creep models developed for normally consolidated clay is discussed and applied to define the development of creep strain in fibrous peat under varying isotropic and deviator stresses. The secondary consolidation coefficient for evaluating the volumetric strain rate of peat is found to be applicable with some limits. The drained creep behaviour of remoulded peat specimens differs from the behaviour shown by Shelby tube specimens, whereas the undrained creep behaviour in remoulded and Shelby tube specimens is similar.     

An orthotropic interface damage model for simulating drying processes in soils

The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory.     

Accessing effects and signals of leakage from a CO<Subscript>2</Subscript> reservoir to a shallow freshwater aquifer by reactive transport modelling

CO₂ geological storage is a transitional technology for the mitigation of climate change. In the vicinity of potential CO₂ reservoirs in Hungary, protected freshwater aquifers used for drinking water supplies exist. Effects of disaster events of CO₂ escape and brine displacement to one of these aquifers have been studied by kinetic 1D reactive transport modelling in PHREEQC. Besides verifying that ion concentrations in the freshwater may increase up to drinking water limit values in both scenarios (CO₂ or brine leakage), total porosity of the rock is estimated. Pore volume is expected to increase at the entry point of CO₂ and to decrease at further distances, whereas it shows minor increase along the flow path for the effect of brine inflow. Additionally, electrical conductivity of water is estimated and suggested to be the best parameter to measure for cost-effective monitoring of both worst-case leakage scenarios.     

Bioerosion and encrustation rates in recent mollusk death assemblages on a supratidal siliciclastic setting,Playa Norte,Veracruz State,Mexico

This study presents both qualitative and quantitative data regarding marine mollusk (gastropods and bivalves) shell bioerosion and encrustation based on death assemblages obtained from a recent supratidal environment in Playa Norte, Veracruz State. The objectives of this study were to assess the nature of bioerosion and encrustation processes and to investigate the role of these taphonomic features contributing to the deterioration of natural shell accumulations within a tropical siliciclastic tidal environment. The assemblage comprises 31 species: 13 gastropods and 18 bivalves. The bioerosion and encrustation degrees were low to moderate for both types. The most abundant traces were predatory gastropod structures (Oichnus paraboloides and O. simplex), whereas sponge borings (Entobia isp.), polychaete dwellings (Caulostrepsis taeniola), and echinoid raspings (Gnatichnus isp.) were less frequent. The encrusting organisms include polychaete serpulids, bryozoans, and rare foraminifers (Homotrema rubrum). Because of the low bioerosion and encrustation degrees occurring in this area, accumulation is expected to predominate over biotic destruction. As deposition conditions (richly fossiliferous carbonate sandstone beds) were similar to those prevailing in the Tuxpan Formation during the Miocene (Langhian), it is suggested that this study provides an equivalent reference to interpret mollusk fossil assemblages located in this site.     

Exploratory Factor Analysis of Wireline Logs Using a Float-Encoded Genetic Algorithm

In the paper, a novel inversion approach is used for the solution of the problem of factor analysis. The float-encoded genetic algorithm as a global optimization method is implemented to extract factor variables using open-hole logging data. The suggested statistical workflow is used to give a reliable estimate for not only the factors but also the related petrophysical properties in hydrocarbon formations. In the first step, the factor loadings and scores are estimated by Jöreskog’s fast approximate method, which are gradually improved by the genetic algorithm. The forward problem is solved to calculate wireline logs directly from the factor scores. In each generation, the observed and calculated well logs are compared to update the factor population. During the genetic algorithm run, the average fitness of factor populations is maximized to give the best fit between the observed and theoretical data. By using the empirical relation between the first factor and formation shaliness, the shale volume is estimated along the borehole. Permeability as a derived quantity also correlates with the first factor, which allows its determination from an independent source. The estimation results agree well with those of independent deterministic modeling and core measurements. Case studies from Hungary and the USA demonstrate the feasibility of the global optimization based factor analysis, which provides a useful tool for improved reservoir characterization.