Geochronology of granulites from the south Itabuna-Salvador-Curaçá Block,São Francisco Craton (Brazil): Nd isotopes and U–Pb zircon ages

This work provides five new U–Pb zircon dating and the corresponding Nd isotope data for felsic granulites from the south Itabuna-Salvador-Curaçá Block (ISCB), in the São Francisco Craton, Brazil. Three major sets of felsic granulites can be recognised. The oldest set is tonalitic in composition and of TTG affinity. It is Archaean in age with magmatic zircon cores dated at 2675 ± 11 Ma by LA-ICPMS and up to ca 2.7–2.9 Ga by SHRIMP on an other sample. It exhibits epsilon Nd values between ?8 and ?11 at 2.1 Ga. This Nd signature is similar to that of granulites found in the western Archaean Jequié Block. Cartographically, this set of Archaean terrains represents at least 50% of the granulites in the studied area. The second set corresponds to a Palaeoproterozoic calc-alkaline tonalitic suite with zircon ages from 2019 ± 19 Ma to 2191 ± 10 Ma and epsilon Nd values between ?3 and ?4 at 2.1 Ga, corresponding partially to a newly formed crust. The third set of granulites is also Palaeoproterozoic. It is shoshonitic to monzonitic in composition and synchronous with the high grade metamorphism dated by metamorphic zircons at 2086 ± 7 Ma (average of five samples). The Nd isotope signature for this alkaline set is similar to that of the Palaeoproterozoic calc-alkaline one. Nd isotopes appear to be a very efficient tool to distinguish Archaean from Palaeoproterozoic felsic protoliths in granulitic suites of the Itabuna-Salvador-Curaçá Block (ISCB). Finally, the southern part of the ISCB is composed of a mixture of Archaean and Palaeoproterozoic protoliths, in similar amounts, suggesting that it was probably an active margin between 2.1 and 2.2 Ga located on the eastern border of the Archaean Jequié Block. A major crustal thickening process occurred at ca 2.09 Ga in the ISCB and seems significantly younger towards the west, in the Jequié granulites, where an average of 2056 ± 9 Ma is determined for the high grade event.     

Differential response at the seafloor during Palaeocene and Eocene ocean warming events at Walvis Ridge,Atlantic Ocean (ODP Site 1262)

The Latest Danian Event (LDE, c. 62.1 Ma) is an early Palaeogene hyperthermal or transient (<200 ka) ocean warming event. We present the first deep‐sea benthic foraminiferal faunal record to study deep‐sea biotic changes together with new benthic (Nuttallides truempyi) stable isotope data from Walvis Ridge Site 1262 (Atlantic Ocean) to evaluate whether the LDE was controlled by similar processes as the minor early Eocene hyperthermals. The spacing of the double negative δ¹³C and δ¹⁸O excursion and the slope of the δ¹⁸O–δ¹³C regression are comparable, strongly suggesting a similar orbital control and pacing of eccentricity maxima as well as a rather homogeneous carbon pool. However, in contrast to early Eocene hyperthermals, the LDE exhibits a remarkable stability of the benthic foraminiferal fauna. This lack of benthic response could be related to the absence of threshold‐related circulation changes or better pre‐adaptation to elevated deep‐sea temperatures, as the LDE was superimposed on a cooling trend, in contrast to early Eocene warming.     

Plio-Pleistocene mammals from Mille-Logya,Ethiopia, and the post-Hadar faunal change

We describe the non-primate mammalian fauna from the late Pliocene to earliest Pleistocene deposits of Mille-Logya in the Lower Awash Valley, Ethiopia, dated to c. 2.9–2.4 Ma, and divided into three successive units: Gafura, Seraitu, and Uraitele. We identify 41 mammalian taxa (including rodents), the most diverse group being the Bovidae, with 17 taxa. While the Gafura assemblage still resembles those from the earlier Hadar Formation, the younger Seraitu assemblage documents a major turnover. While there is little change in the species present across this interval, the relative abundances of various taxa change dramatically, with suids being largely replaced by open-country bovids (Alcelaphini and Antilopini). We interpret this faunal change as reflective of an environmental shift, contemporaneous with the replacement of Australopithecus afarensis by Homo in the area.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part 1: Saturated state

The behavior of granular materials is known to depend on its loose or dense nature, which in turns depends both on density and confining pressure. Many models developed in the past require the use of different sets of constitutive parameters for the same material under different confining pressures. The purpose of this paper is to extend a basic generalized plasticity model for sands proposed by Pastor, Zienkiewicz and Chan by modifying the main ingredients of the model flow—rule, loading–unloading discriminating direction and plastic modulus—to include a dependency on the state parameter. The proposed model is tested against the available experimental data on three different sands, using for each of them a single set of material parameters, finding a reasonably good agreement between experiments and predictions. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Likely effects of climate change on groundwater availability in a Mediterranean region of Southeastern Spain

Groundwater resources are typically the main fresh water source in arid and semi‐arid regions. Natural recharge of aquifers is mainly based on precipitation; however, only heavy precipitation events (HPEs) are expected to produce appreciable aquifer recharge in these environments. In this work, we used daily precipitation and monthly water level time series from different locations over a Mediterranean region of Southeastern Spain to identify the critical threshold value to define HPEs that lead to appreciable aquifer recharge in this region. Wavelet and trend analyses were used to study the changes in the temporal distribution of the chosen HPEs (≥20 mm day^?1) over the observed period 1953–2012 and its projected evolution by using 18 downscaled climate projections over the projected period 2040–2099. The used precipitation time series were grouped in 10 clusters according to similarities between them assessed by using Pearson correlations. Results showed that the critical HPE threshold for the study area is 20 mm day^?1. Wavelet analysis showed that observed significant seasonal and annual peaks in global wavelet spectrum in the first sub‐period (1953–1982) are no longer significant in the second sub‐period (1983–2012) in the major part of the ten clusters. This change is because of the reduction of the mean HPEs number, which showed a negative trend over the observed period in nine clusters and was significant in five of them. However, the mean size of HPEs showed a positive trend in six clusters. A similar tendency of change is expected over the projected period. The expected reduction of the mean HPEs number is two times higher under the high climate scenario (RCP8.5) than under the moderate scenario (RCP4.5). The mean size of these events is expected to increase under the two scenarios. The groundwater availability will be affected by the reduction of HPE number which will increase the length of no aquifer recharge periods (NARP) accentuating the groundwater drought in the region. Copyright © 2016 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.