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.     

Stability analysis of tunnels driven in jointed rocks by means of a homogenized limit analysis approach

Conceived as a potential alternative to the classical design methods employed for analyzing the stability of underground structures driven in jointed rocks, the homogenization approach stems from the heuristic idea that, from a macroscopic point of view, a rock mass cut by a network of joints may be perceived as a homogenized continuum. The strength properties of the latter can be theoretically obtained from the failure conditions of its individual constituents: rock matrix and joint interfaces. At the material level, the limit analysis reasoning is used in the context of homogenization to formulate the homogenized strength criterion of a jointed rock mass in the particular situation of a single set of parallel joints. As it could be expected, the obtained closed‐form expressions show the strength anisotropy induced by joint preferential orientation. The support functions (π functions) associated with the homogenized strength criterion are also determined in both plane strain and three‐dimensional cases. This criterion is then applied to the investigation of stability analysis of a tunnel excavated in a jointed rock mass. Upper bounds estimated of the stability factor are derived from the implementation of the kinematic approach directly on the homogenized underground structure. Finally, the approach is applied to analyze and discuss the collapse of the Pinheiros subway station (São Paulo, Brazil). Copyright © 2014 John Wiley & Sons, Ltd.     

Elastoplastic model with loading memory surfaces (LMS) for monotonic and cyclic behaviour of geomaterials

A new elastoplastic model called loading memory surface based on the critical state concept and the multi‐surface framework is proposed for geomaterials. The model uses a hypoelastic formulation and two plastic mechanisms. The formulations of the model are made in three‐dimensional stress–strain space and work under both monotonic and cyclic loadings. A newly introduced formalism makes it possible to obtain the cyclic response directly from the monotonic loading one. This formalism gives a three‐dimensional generalization of the well‐known Masing rule. The model has been validated against test results of Hostun sand under several conditions: monotonic and cyclic, drained and undrained, tests in compression and in extension, and at different confining pressures and different densities. Copyright © 2014 John Wiley & Sons, Ltd.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

Segmentation and kinematics of the North America‐Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high‐resolution multibeam echo‐sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate‐boundary structures are a series of strike‐slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre‐existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike‐slip regime. Along the most recent trace of the SOFZ, we measured a strike‐slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS‐derived motion of 9.8 ± 2 mm a^?1 has remained stable during the entire Quaternary.     

Unraveling complex hydrogeological processes in Andean basins in south‐central Chile: An integrated assessment to understand hydrological dissimilarity

Groundwater storage, drainage, and interbasin water exchange are common hydrological processes but often difficult to quantify due to a lack of local observations. We present a study of three volcanic mountainous watersheds located in south‐central Chile (~36.9 ° S) in the Chillán volcanic complex (Chillán, Renegado, and Diguillín river basins). These are neighboring basins that are similar with respect to the metrics normally available for characterization everywhere (e.g., precipitation, temperature, and land cover). In a hydrological sense, similar (proportional) behavior would be expected if these catchments would be characterized with this general information. However, these watersheds show dissimilar behavior when analyzed in detail. The surface water balance does not fit for any of these watersheds individually; however, the water balance of the whole system can be explained by likely interbasin water exchanges. The Renegado river basin has an average annual runoff per unit of area on the order of 60–65% less than those of the Diguillín and Chillán rivers, which is contradictory to the hydrological similarity among the basins. To understand the main processes that control streamflow generation, two analyses were performed: (a) basin metrics (land cover, geologic, topographic, and climatological maps) and hydro‐meteorological data analyses and (b) a water balance model approach. The analyses contribute to a plausible explanation for the hydrogeological processes in the system. The soils, topography, and geology of the Chillán–Renegado–Diguillín system favor the infiltration and groundwater movements from the Renegado river basin, mainly to the neighboring Diguillín basin. The interbasin water exchanges affect hydrological similarity and explain the differences observed in the hydrological processes of these three apparently similar volcanic basins. The results highlight the complexity of hydrological processes in volcanic mountainous systems and suggest that a simple watershed classification approach based on widely available data is insufficient. Simple local analyses such as specific flow analysis with a review of the geology and morphology can contribute to a better understanding of the hydrology of volcanic mountainous areas.     

Soil cracking effects on hydrological and erosive processes: a study case in Mediterranean cultivated vertisols

Shrink–swell soils, such as those in a Mediterranean climate regime, can cause changes in terms of hydrological and erosive responses because of the changing soil water storage conditions. Only a limited number of long‐term studies have focused on the impacts on both hydrological and erosive responses and their interactions in an agricultural environment. In this context, this study aims to document the dynamics of cracks, runoff and soil erosion within a small Mediterranean cultivated catchment and to quantify the influence of crack processes on the water and sediment supplied to a reservoir located at the catchment outlet using water and sediment measurements at a cultivated field outlet as baseline. Detailed monitoring of the presence of topsoil cracks was conducted within the Kamech catchment (ORE OMERE, Tunisia), and runoff and suspended sediment loads were continuously measured over a long period of time (2005–2012) at the outlets of a field (1.32 ha) and a catchment (263 ha). Analysis of the data showed that topsoil cracks were open approximately half of the year and that the rainfall regime and water table level conditions locally control the seasonal cracking dynamics. Topsoil cracks appeared to seriously affect the generation of runoff and sediment concentrations and, consequently, sediment yields, with similar dynamics observed at the field and catchment outlets. A similar time lag in the seasonality between water and sediment delivery was observed at these two scales: although the runoff rates were globally low during the presence of topsoil cracks, most sediment transport occurred during this period associated with very high sediment concentrations. This study underlines the importance of a good prediction of runoff during the presence of cracks for reservoir siltation considerations. In this context, the prediction of cracking effects on runoff and soil erosion is a key factor for the development of effective soil and water management strategies and downstream reservoir preservation. Copyright © 2016 John Wiley & Sons, Ltd.     

New major element analyses of proximal tephras from the Azores and suggested correlations with cryptotephras in North-West Europe

The Azores Archipelago is one of the most active volcanic areas in the North Atlantic region. Approximately 30 eruptions have been reported over the last 600 years with some major VEI 5 (Volcanic Explosivity Index) eruptions further back in time. The geochemical composition of associated tephra-derived glass, however, is not well characterized. An Azorean origin of cryptotephras found in distal areas such as North Africa, the British Isles and Greenland has been suggested, but proximal data from the Azores are scarce and the correlations have only been tentative. These tephras have a traychtic composition, which excludes an Icelandic origin. In a previous study, we presented major element analyses of proximal tephra-derived glass from five Holocene eruptions on the Azores Islands. There is a striking geochemical similarity between tephras from volcanoes on São Miguel and Irish cryptotephras, and especially with eruptives from the Furnas volcano. Here we present new analyses of proximal tephras that confirm and strengthen a link between Furnas and cryptotephras found in south-west Ireland. We also suggest a correlation between a previously unsourced tephra found in a Swedish bog with an eruption of the Sete Cidades volcano c. 3880 a cal BP.