Strength criterion for cross-anisotropic sand under general stress conditions

By incorporating the fabric effect and Lode’s angle dependence into the Mohr–Coulomb failure criterion, a strength criterion for cross-anisotropic sand under general stress conditions was proposed. The obtained criterion has only three material parameters which can be specified by conventional triaxial tests. The formula to calculate the friction angle under any loading direction and intermediate principal stress ratio condition was deduced, and the influence of the degree of the cross-anisotropy was quantified. The friction angles of sand in triaxial, true triaxial, and hollow cylinder torsional shear tests were obtained, and a parametric analysis was used to detect the varying characteristics. The friction angle becomes smaller when the major principal stress changes from perpendicular to parallel to the bedding plane. The loading direction and intermediate principal stress ratio are unrelated in true triaxial tests, and their influences on the friction angle can be well captured by the proposed criterion. In hollow cylinder torsional shear tests with the same internal and external pressures, the loading direction and intermediate principal stress ratio are related. This property results in a lower friction angle in the hollow cylinder torsional shear test than that in the true triaxial test under the same intermediate principal stress ratio condition. By comparing the calculated friction angle with the experimental results under various loading conditions (e.g., triaxial, true triaxial, and hollow cylinder torsional shear test), the proposed criterion was verified to be able to characterize the shear strength of cross-anisotropic sand under general stress conditions.     

Identifying influence areas with connectivity analysis – application to the local perturbation of heterogeneity distribution for history matching

Numerical representations of a target reservoir can help to assess the potential of different development plans. To be as predictive as possible, these representations or models must reproduce the data (static, dynamic) collected on the field. However, constraining reservoir models to dynamic data – the history-matching process – can be very time consuming. Many uncertain parameters need to be taken into account, such as the spatial distribution of petrophysical properties. This distribution is mostly unknown and usually represented by millions of values populating the reservoir grid. Dedicated parameterization techniques make it possible to investigate many spatial distributions from a small number of parameters. The efficiency of the matching process can be improved from the perturbation of specific regions of the reservoir. Distinct approaches can be considered to define such regions. For instance, one can refer to streamlines. The leading idea is to identify areas that influence the production behavior where the data are poorly reproduced. Here, we propose alternative methods based on connectivity analysis to easily provide approximate influence areas for any fluid-flow simulation. The reservoir is viewed as a set of nodes connected by weighted links that characterize the distance between two nodes. The path between nodes (or grid blocks) with the lowest cumulative weight yields an approximate flow path used to define influence areas. The potential of the approach is demonstrated on the basis of 2D synthetic cases for the joint integration of production and 4D saturation data, considering several formulations for the weights attributed to the links.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

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...     

Pico- and Nanophytoplankton Dynamics in Two Coupled but Contrasting Coastal Bays in the NW Mediterranean Sea (France)

Due to its ecological context, the Toulon bay represents a site of scientific interest to study temporal plankton distribution, particularly pico- and nanophytoplankton dynamics. A monthly monitoring was performed during a two-year cycle (October 2013–December 2015) at two coupled sampling sites, referred to as Little and Large bays, which had different morphometric characteristics and human pressures. Flow cytometry analyses highlighted the fact that pico- and nanophytoplankton were more abundant in the eutrophic Little bay. Furthermore, it evidenced two community structures across the Toulon bays: at times, a co-dominance of picoeukaryotes, nanoeukaryotes, Synechococcus 1-like cells and Prochlorococcus-like cells was found, and at other times, a Synechococcus 1-like dominated community existed. The alternation of one structure or the other can be explained by a combined action of temperature regime, nutrient conditions and degree of contamination. This study showed that pico- and nanophytoplankton dynamics were mainly driven by temperature in both sites, as in other temperate Mediterranean regions. Thus, the community was mainly composed of picoeukaryotes and Prochlorococcus-like cells in the winter (<?15 °C), while it was dominated by Synechococcus 1-like cells in the summer (>?20 °C). Additionally, the multiple human stressors in the Little bay seemed to affect the increase in abundance of Synechococcus 1-like cells as they were preferentially observed in the Large bay.     

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.     

Assessing rockfall susceptibility in steep and overhanging slopes using three-dimensional analysis of failure mechanisms

Rockfalls strongly influence the evolution of steep rocky landscapes and represent a significant hazard in mountainous areas. Defining the most probable future rockfall source areas is of primary importance for both geomorphological investigations and hazard assessment. Thus, a need exists to understand which areas of a steep cliff are more likely to be affected by a rockfall. An important analytical gap exists between regional rockfall susceptibility studies and block-specific geomechanical calculations. Here we present methods for quantifying rockfall susceptibility at the cliff scale, which is suitable for sub-regional hazard assessment (hundreds to thousands of square meters). Our methods use three-dimensional point clouds acquired by terrestrial laser scanning to quantify the fracture patterns and compute failure mechanisms for planar, wedge, and toppling failures on vertical and overhanging rock walls. As a part of this work, we developed a rockfall susceptibility index for each type of failure mechanism according to the interaction between the discontinuities and the local cliff orientation. The susceptibility for slope parallel exfoliation-type failures, which are generally hard to identify, is partly captured by planar and toppling susceptibility indexes. We tested the methods for detecting the most susceptible rockfall source areas on two famously steep landscapes, Yosemite Valley (California, USA) and the Drus in the Mont-Blanc massif (France). Our rockfall susceptibility models show good correspondence with active rockfall sources. The methods offer new tools for investigating rockfall hazard and improving our understanding of rockfall processes.     

Subduction-related mafic to felsic magmatism in the Malayer–Boroujerd plutonic complex,western Iran

The Malayer–Boroujerd plutonic complex (MBPC) in western Iran, consists of a portion of a magmatic arc built by the northeast verging subduction of the Neo-Tethys plate beneath the Central Iranian Microcontinent (CIMC). Middle Jurassic-aged felsic magmatic activity in MBPC is manifested by I-type and S-type granites. The mafic rocks include gabbroic intrusions and dykes and intermediate rocks are dioritic dykes and minor intrusions, as well as mafic microgranular enclaves (MMEs). MBPC Jurassic-aged rocks exhibit arc-like geochemical signatures, as they are LILE- and LREE-enriched and HFSE- and HREE-depleted and display negative Nb–Ta anomalies. The gabbro dykes and intrusions originated from metasomatically enriched garnet-spinel lherzolite [Degree of melting (f_mel) ~ 15%] and exhibit negative Nd and positive to slightly negative ε_Hf(T) (+ 3.0 to ? 1.6). The data reveal that evolution of Middle Jurassic magmatism occurred in two stages: (1) deep mantle-crust interplay zone and (2) the shallow level upper crustal magma chamber. The geochemical and isotopic data, as well as trace element modeling, indicate the parent magma for the MBPC S-type granites are products of upper crustal greywacke (f_mel: 0.2), while I-type granites formed by partial melting of amphibolitic lower crust (f_mel: 0.25) and mixing with upper crustal greywacke melt in a shallow level magma chamber [Degree of mixing (f_mix): 0.3]. Mixing between andesitic melt leaving behind a refractory dense cumulates during partial crystallization of mantle-derived magma and lower crustal partial melt most likely produced MMEs (f_mix: 0.2). However, enriched and moderately variable ε_Nd(T) (? 3.21 to ? 4.33) and high (⁸⁷Sr/⁸⁶Sr)_i (0.7085–0.7092) in dioritic intrusions indicate that these magmas are likely experienced assimilation of upper crustal materials. The interpretations of magmatic activity in the MBPC is consistent with the role considered for mantle-derived magma as heat and mass supplier for initiation and evolution of magmatism in continental arc setting, elsewhere.