Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

Impacts of Indian Ocean SST biases on the Indian Monsoon: as simulated in a global coupled model

In this study, the impact of the ocean–atmosphere coupling on the atmospheric mean state over the Indian Ocean and the Indian Summer Monsoon (ISM) is examined in the framework of the SINTEX-F2 coupled model through forced and coupled control simulations and several sensitivity coupled experiments. During boreal winter and spring, most of the Indian Ocean biases are common in forced and coupled simulations, suggesting that the errors originate from the atmospheric model, especially a dry islands bias in the Maritime Continent. During boreal summer, the air-sea coupling decreases the ISM rainfall over South India and the monsoon strength to realistic amplitude, but at the expense of important degradations of the rainfall and Sea Surface Temperature (SST) mean states in the Indian Ocean. Strong SST biases of opposite sign are observed over the western (WIO) and eastern (EIO) tropical Indian Ocean. Rainfall amounts over the ocean (land) are systematically higher (lower) in the northern hemisphere and the south equatorial Indian Ocean rainfall band is missing in the control coupled simulation. During boreal fall, positive dipole-like errors emerge in the mean state of the coupled model, with warm and wet (cold and dry) biases in the WIO (EIO), suggesting again a significant impact of the SST errors. The exact contributions and the distinct roles of these SST errors in the seasonal mean atmospheric state of the coupled model have been further assessed with two sensitivity coupled experiments, in which the SST biases are replaced by observed climatology either in the WIO (warm bias) or EIO (cold bias). The correction of the WIO warm bias leads to a global decrease of rainfall in the monsoon region, which confirms that the WIO is an important source of moisture for the ISM. On the other hand, the correction of the EIO cold bias leads to a global improvement of precipitation and circulation mean state during summer and fall. Nevertheless, all these improvements due to SST corrections seem drastically limited by the atmosphere intrinsic biases, including prominently the unimodal oceanic position of the ITCZ (Inter Tropical Convergence Zone) during summer and the enhanced westward wind stress along the equator during fall.     

Long-term groundwater resource observatory for Southwestern Madagascar

Madagascar has one of the highest poverty rates in the world and consequently the long-term monitoring of groundwater resources is not a priority for the authorities. However, groundwater is often the only sustainable resource that has a satisfactory quality to supply the population. This is especially true in the south-west of the country, which is a semi-arid region and a global change hot spot (intense land use and climate changes). In response to the lack of data, the Groundwater Resource Observatory for Southwestern Madagascar (GROSoM) was established to monitor piezometry and meteorology over the longer term as part of a humanitarian response. The first site was setup in 2014 in a catchment located over a carbonate plateau; in 2018, a second site was installed in an alluvial setting within a crystalline basement catchment and a third site will be installed in 2020 to monitor groundwater dynamics in a coastal setting. The three sites, located between Toliara and Taolagnaro cities, are complementary and representative of various hydrogeological systems in Southwestern Madagascar. Each site includes a weather station and between 3 and 6 piezometric probes. The monitoring data indicate a strong inter-annual variability in precipitation, which induces a strong variability in aquifers recharge. One of the driest years in 2016 seems to be consistent with strong El Niño – Southern Oscillation (ENSO) effects observed at the global scale, while years with higher recharge appear to be related to cyclones such as Fundi in 2015 and Eketsang in 2019. Preliminary results of cross-disciplinary studies demonstrated a link between groundwater and health issues (i.e., admissions to basic health centres). This observatory aims to produce long-term data and has two objectives: (i) strengthening the early warning system for humanitarian crises in Madagascar; (ii) contributing to a better understanding of the effects of climate change on groundwater resources in this semi-arid region.     

Tectonometamorphic evolution of the Atbashi high‐P units (Kyrgyz CAOB,Tien Shan): Implications for the closure of the Turkestan Ocean and continental subduction–exhumation of the South Kazakh continental margin

The South Tien Shan (STS) belt results from the last collision event in the western Central Asian Orogenic Belt (CAOB). Understanding its formation is of prime importance in the general framework of the CAOB. The Atbashi Range preserves high‐P (HP) rocks along the STS suture, but still, its global metamorphic evolution remains poorly constrained. Several HP units have been identified: (a) a HP tectonic mélange including boudins of mafic eclogites in a sedimentary matrix, (b) a large (>100 km long) high‐P metasedimentary unit (HPMU) and (c) a lower blueschist facies accretionary prism. Raman Spectroscopy on carbonaceous material combined with phengite and chlorite multiequilibria and isochemical phase diagram modelling indicates that the HPMU recorded homogeneous P–T conditions of 23–25 kbar and 560–570°C along the whole unit. ⁴⁰Ar/³⁹Ar dating on phengite from the HPMU ranges between 328 and 319 Ma at regional scale. These ages are interpreted as (re‐) crystallization ages of phengite during T_max conditions at a pressure range of 20–25 kbar. Thermobarometry on samples from the HP tectonic mélange provides similar metamorphic peak conditions. Thermobarometry on the blueschist to lower greenschist facies accretionary prism indicates that it underwent P–T conditions of 5–6 kbar and 290–340°C, highlighting a 17–20 kbar pressure gap between the HPMU‐tectonic mélange units and the accretionary prism. Comparison with available geochronological data suggests a very short time span between the prograde path (340 Ma), HP metamorphic peak (330 Ma), the T_max (328–319 Ma) and the final exhumation of the HPMU (303–295 Ma). Extrusion of the HPMU, accommodated by a basal thrust and an upper detachment, was driven by buoyant forces from 70–75 km up to 60 km depth, which directly followed continental subduction and detachment of the HPMU. At crustal depths, extrusion was controlled by collisional tectonics up to shallow levels. Lithological homogeneity of the HPMU and its continental‐derived character from the North Tien Shan suggest this unit corresponds to the hyper‐extended continental margin of the Kazakh continent, subducted southward below the north continental active margin of the Tarim craton. Integration of the available geological data allows us to propose a general geodynamic scenario for Tien Shan during the Carboniferous with a combination of (a) N‐dipping subduction below the Kazakh margin of Middle Tien Shan until 390–340 Ma and (b) S‐dipping subduction of remaining Turkestan marginal basins between 340 and 320 Ma.     

A geochemical approach to constraining the formation of glassy fallout debris from nuclear tests

Glassy nuclear fallout debris from near-surface nuclear tests is fundamentally reprocessed earth material. A geochemical approach to analysis of glassy fallout is uniquely suited to determine the means of reprocessing and shed light on the mechanisms of fallout formation. An improved understanding of fallout formation is of interest both for its potential to guide post-detonation nuclear forensic investigations and in the context of possible affinities between glassy debris and other glasses generated by high-energy natural events, such as meteorite impacts and lightning strikes. This study presents a large major-element compositional dataset for glasses within aerodynamic fallout from the Trinity nuclear test (“trinitite”) and a geochemically based analysis of the glass compositional trends. Silica-rich and alkali-rich trinitite glasses show compositions and textures consistent with formation through melting of individual mineral grains—quartz and alkali feldspar, respectively—from the test-site sediment. The volumetrically dominant glass phase—called the CaMgFe glass—shows extreme major-element compositional variability. Compositional trends in the CaMgFe glass are most consistent with formation through volatility-controlled condensation from compositionally heterogeneous plasma. Radioactivity occurs only in CaMgFe glass, indicating that co-condensation of evaporated bulk ground material and trace device material was the main mechanism of radioisotope incorporation into trinitite. CaMgFe trinitite glasses overlap compositionally with basalts, rhyolites, fulgurites, tektites, and microtektites but display greater compositional diversity than all of these naturally formed glasses. Indeed, the most refractory CaMgFe glasses compositionally resemble early solar system condensates—specifically, CAIs.     

Influence of damage on pore size distribution and permeability of rocks

The model proposed in this article relates permeability to porosity measurements that can easily be performed in the laboratory. The pore size distribution (PSD) curve is updated with strains and damage. The updated volumetric fractions of natural pores and cracks are introduced in the expression of permeability. Contrary to classical permeability models based on PSD integrations, the model proposed in this article accounts for possible changes in the porosity modes: one mode for undamaged samples and two modes for cracked samples. The proposed approach also accounts for varying states of damage, as opposed to classical fracture network models, in which the cracks pattern is fixed. The only material parameters that are required to describe the microstructure are the lower and upper bounds of the pores size for both natural pores and cracks. All the other PSD parameters involved in the model are related to macroscopic parameters that can easily be determined in the laboratory, such as the initial void ratio. The framework proposed in this article can be used in any damage constitutive model to determine the permeability of a brittle porous medium. Drained triaxial compression tests have been simulated. Before cracks initiation, permeability decreases while the larger natural pores are getting squeezed. After the occurrence of damage, permeability grows due to the increase of cracks density. The model performs well to represent the influence of the confining pressure on damage evolution and permeability variations. Copyright © 2012 John Wiley & Sons, Ltd.