Reply to ‘Comment on “The beginnings of hydrous mantle wedge melting” by Till et al.’ by Green, Rosenthal and Kovacs

The comment of Green et al. debates the interpretation of the temperature of the H₂O-saturated peridotite solidus and presence of silicate melt in the experiments of Till et al. (Contrib Mineral Petrol 163:669–688, 2012) at <1,000?°C. The criticisms presented in their comment do not invalidate any of the most compelling observations of Till et al. (Contrib Mineral Petrol 163:669–688, 2012) as discussed in the following response, including the changing minor element and Mg# composition of the solid phases with increasing temperature in our experiments with 14.5?wt% H₂O at 3.2?GPa, as well as the results of our chlorite peridotite melting experiments with 0.7?wt% H₂O. The point remains that Till et al. (Contrib Mineral Petrol 163:669–688, 2012) present data that call into question the H₂O-saturated peridotite solidus temperature preferred by Green (Tectonophysics 13(1–4):47–71, 1972; Earth Planet Sci Lett 19(1):37–53, 1973; Can Miner 14:255–268, 1976); Millhollen et al. (J Geol 82(5):575–587, 1974); Mengel and Green (Stability of amphibole and phlogopite in metasomatized peridotite under water-saturated and water-undersaturated conditions, Geological Society of Australia Special Publication, Blackwell, pp 571-581, 1989); Wallace and Green (Mineral Petrol 44:1–19, 1991) and Green et al. (Nature 467(7314):448–451, 2010).     

Applicability of different ground-motion prediction models for northern Iran

A total of 163 free-field acceleration time histories recorded at epicentral distances of up to 200 km from 32 earthquakes with moment magnitudes ranging from M _w 4.9 to 7.4 have been used to investigate the predictive capabilities of the local, regional, and next generation attenuation (NGA) ground-motion prediction equations and determine their applicability for northern Iran. Two different statistical approaches, namely the likelihood method (LH) of Scherbaum et al. (Bull Seismol Soc Am 94:341–348, 2004) and the average log-likelihood method (LLH) of Scherbaum et al. (Bull Seismol Soc Am 99:3234–3247, 2009), have been applied for evaluation of these models. The best-fitting models (considering both the LH and LLH results) over the entire frequency range of interest are those of Ghasemi et al. (Seismol 13:499–515, 2009a) and Soghrat et al. (Geophys J Int 188:645–679, 2012) among the local models, Abrahamson and Silva (Earthq Spectra 24:67–97, 2008) and Chiou and Youngs (Earthq Spectra 24:173–215, 2008) among the NGA models, and finally Akkar and Bommer (Seism Res Lett 81:195–206, 2010) among the regional models.     

EPR study of chromium-doped forsterite crystals: Cr3+(M1) with associated trivalent ions Al3+ and Sc3+

Electron paramagnetic resonance (EPR) study of single crystals of forsterite co-doped with chromium and scandium has revealed, apart from the known paramagnetic centers Cr³⁺(M1) and Cr³⁺(M1)– $ V_{{{\text{Mg}}^{2 + } }} $ (M2) (Ryabov in Phys Chem Miner 38:177–184, 2011), a new center Cr³⁺(M1)– $ V_{{{\text{Mg}}^{2 + } }} $ (M2)–Sc³⁺ formed by a Cr³⁺ ion substituting for Mg²⁺ at the M1 structural position with a nearest-neighbor Mg²⁺ vacancy at the M2 position and a Sc³⁺ ion presumably at the nearest-neighbor M1 position. For this center, the conventional zero-field splitting parameters D and E and the principal g values have been determined as follows: D?=?33,172(29) MHz, E?=?8,482(13) MHz, g?=?[1.9808(2), 1.9778(2), 1.9739(2)]. The center has been compared with the known ion pair Cr³⁺(M1)–Al³⁺ (Bershov et al. in Phys Chem Miner 9:95–101, 1983), for which the refined EPR data have been obtained. Based on these data, the known sharp M1″ line at 13,967?cm^?1 (with the splitting of 1.8?cm^?1), observed in low-temperature luminescence spectra of chromium-doped forsterite crystals (Glynn et al. in J Lumin 48, 49:541–544, 1991), has been ascribed to the Cr³⁺(M1)–Al³⁺ center. It has been found that the concentration of the new center increases from 0 up to 4.4?×?10¹⁵?mg^?1, whereas that of the Cr³⁺(M1) and Cr³⁺(M1)– $ V_{{{\text{Mg}}^{2 + } }} $ (M2) centers quickly decreases from 7.4?×?10¹⁵?mg^?1 down to 3?×?10¹⁵?mg^?1 and from 2.7?×?10¹⁵?mg^?1 down to 0.5?×?10¹⁵?mg^?1, i.e., by a factor of 2.5 and 5.4, respectively, with an increase of the Sc content from 0 up to 0.22 wt?% (at the same Cr content 0.25 wt?%) in the melt. When the Sc content exceeds that of Cr, the concentration of the new center decreases most likely due to the formation of the Sc³⁺(M1)– $ V_{{{\text{Mg}}^{2 + } }} $ (M2)–Sc³⁺ complex instead of the Cr³⁺(M1)– $ V_{{{\text{Mg}}^{2 + } }} $ (M2)–Sc³⁺ center. The formation of such ordered neutral complex is in agreement with the experimental results, concerning the incorporation of Sc into olivine, recently obtained by Grant and Wood (Geochim Cosmochim Acta 74:2412–2428, 2010).     

Paleogeographical restoration and ramp tectonic evidence in Tunisian Tellian domain: Ain El Bey-Bou Awen area

The imbrication’s area in northern Tunisia is the most external segment of Alpine range, where several associated folds types with thrust ramps are recognized within imbricate units beneath Numidian front slope. Their presence help to understand thrusting mechanisms installation through studied area. In fact, this zone was considered as a result of Paleogene gravitary slop (Kujawski (Ann Miner Géol Tunis (24):281, 1969); Caire (Ann Min Géol Tunis 26:87–110, 1973); Rouvier 1977), which is proved to be affected by major deep decollement, given rise to various structures, some are propagation folds, specific of foreland front, limited to this area, and those in more external position: Tunisian Atlas (Creusot et al. (C R Acad Sci Paris 314(Sér II):961–965, 1992); Ouali and Mercier (PII: S0191-8141(97):00048-5, 1997); Ouali 1984; Ahmadi et al. (J Struct Geol 28:721–728, 2006)). Various categories of fold ramps could be identified: frontal folds ramp NE–SW and others as lateral or oblique ramp with NW–SE trend (Aridhi et al. (C R Geosci 343:360–369, 2011)). The relation between various structures has been used as recognition tools of thrusting sequences and to propose a new deformation chronology. Delimited outcropping of these structures between two both parallel faults strikes with regional displacement, leads to interpret these faults as cogenetic tear faults of propagation thrusts; this fault separates two domains with different deformation styles from each other side.     

High-temperature and high-pressure equation of state for the hexagonal phase in the system NaAlSiO4 – MgAl2O4

Thermal equation of state of an Al-rich phase with Na_1.13Mg_1.51Al_4.47Si_1.62O₁₂ composition has been derived from in situ X-ray diffraction experiments using synchrotron radiation and a multianvil apparatus at pressures up to 24 GPa and temperatures up to 1,900 K. The Al-rich phase exhibited a hexagonal symmetry throughout the present pressure–temperature conditions and the refined unit-cell parameters at ambient condition were: a=8.729(1) Å, c=2.7695(5) Å, V ₀=182.77(6) Å³ (Z=1; formula weight=420.78 g/mol), yielding the zero-pressure density ρ₀=3.823(1) g/cm³ . A least-square fitting of the pressure-volume-temperature data based on Anderson’s pressure scale of gold (Anderson et al. in J Appl Phys 65:1534–543, 1989) to high-temperature Birch-Murnaghan equation of state yielded the isothermal bulk modulus K ₀=176(2) GPa, its pressure derivative K ₀ ^′ =4.9(3), temperature derivative (?K _T /?T) _P =?0.030(3) GPa K^?1 and thermal expansivity α(T)=3.36(6)×10^?5+7.2(1.9)×10^?9 T, while those values of K ₀=181.7(4) GPa, (?K _T /?T) _P =?0.020(2) GPa K^?1 and α(T)=3.28(7)×10^?5+3.0(9)×10^?9 T were obtained when K ₀ ^′ was assumed to be 4.0. The estimated bulk density of subducting MORB becomes denser with increasing depth as compared with earlier estimates (Ono et al. in Phys Chem Miner 29:527–531 2002; Vanpeteghem et al. in Phys Earth Planet Inter 138:223–230 2003; Guignot and Andrault in Phys Earth Planet Inter 143–44:107–128 2004), although the difference is insignificant (<0.6%) when the proportions of the hexagonal phase in the MORB compositions (～20%) are taken into account.     

Oxygen fugacity control in piston-cylinder experiments: a re-evaluation

Jakobsson (Contrib Miner Petrol 164(3):397–407, 2012) investigated a double capsule assembly for use in piston-cylinder experiments that would allow hydrous, high-temperature, and high-pressure experiments to be conducted under controlled oxygen fugacity conditions. Using a platinum outer capsule containing a metal oxide oxygen buffer (Ni–NiO or Co–CoO) and H₂O, with an inner gold–palladium capsule containing hydrous melt, this study was able to compare the oxygen fugacity imposed by the outer capsule oxygen buffer with an oxygen fugacity estimated by the AuPdFe ternary system calibrated by Barr and Grove (Contrib Miner Petrol 160(5):631–643, 2010). H₂O loss or gain, as well as iron loss to the capsule walls and carbon contamination, is often observed in piston-cylinder experiments and often go unexplained. Only a few have attempted to actually quantify various aspects of these changes (Brooker et al. in Am Miner 83(9–10):985–994, 1998; Truckenbrodt and Johannes in Am Miner 84:1333–1335, 1999). It was one of the goals of Jakobsson (Contrib Miner Petrol 164(3):397–407, 2012) to address these issues by using and testing the AuPdFe solution model of Barr and Grove (Contrib Miner Petrol 160(5):631–643, 2010), as well as to constrain the oxygen fugacity of the inner capsule. The oxygen fugacities of the analyzed melts were assumed to be equal to those of the solid Ni–NiO and Co–CoO buffers, which is incorrect since the melts are all undersaturated in H₂O and the oxygen fugacities should therefore be lower than that of the buffer by 2 log $a_{{{\text{H}}_{ 2} {\text{O}}}}$ .     

Discussion of “The Bushveld Complex,South Africa: formation of platinum–palladium,chrome- and vanadium-rich layers via hydrodynamic sorting of a mobilized cumulate slurry in a large,relatively slowly cooling,subsiding magma chamber” by Maier et al. (2013) Miner Deposita 48:1–56

We believe the hypothesis presented by Maier et al. (Miner Deposita 48:1–56, 2012) for the formation of the various ore bodies in the Bushveld Complex to be overly simplistic, and we raise concerns that some of our work, used in support of this hypothesis, has been misrepresented. The formation of both diverse metalliferous layers (platinum-group element (PGE) reefs and Ti-magnetite layers) and some discordant (pipe) ore deposits has been ascribed by Maier et al. to the single unifying process of hydrodynamic sorting. The problem faced by authors of universal hypotheses for the Bushveld Complex is the sheer size and complexity of the intrusion. We disagree with many aspects of the overall Maier et al. model and have also identified several minor errors on maps and photographs, although some of these do not have a material effect on the model. The nature and origin of the layering is, however, too complex a topic to deal with in the context of this commentary, and we restrict ourselves to noting that our preferred hypothesis, namely the incremental buildup of layering from numerous episodes of replenishment, by different magma lineages, is consistent with field relationships. Our hypothesis for the origin of the ultramafic-hosted PGE-rich reefs, i.e., lateral mixing, is applicable to economically mineralized reefs (Mitchell and Scoon, Econ Geol 102:971–1009, 2007) and poorly mineralized layers such as the Pseudoreef harzburgite (Scoon and De Klerk, Canad Mineral 25:51–77, 1987) and the chromitite layers below the UG2 (Scoon and Teigler, Econ Geol 89:1094–1121, 1994).     

Low-temperature evolution of OH bands in synthetic forsterite, implication for the nature of H defects at high pressure

We performed in situ infrared spectroscopic measurements of OH bands in a forsterite single crystal between ?194 and 200 °C. The crystal was synthesized at 2 GPa from a cooling experiment performed between 1,400 and 1,275 °C at a rate of 1 °C per hour under high silica-activity conditions. Twenty-four individual bands were identified at low temperature. Three different groups can be distinguished: (1) Most of the OH bands between 3,300 and 3,650 cm^?1 display a small frequency lowering (<4 cm^?1) and a moderate broadening (<10 cm^?1) as temperature is increased from ?194 to 200 °C. The behaviour of these bands is compatible with weakly H-bonded OH groups associated with hydrogen substitution into silicon tetrahedra; (2) In the same frequency range, two bands at 3,617 and 3,566 cm^?1 display a significantly anharmonic behaviour with stronger frequency lowering (42 and 27 cm^?1 respectively) and broadening (~30 cm^?1) with increasing temperature. It is tentatively proposed that the defects responsible for these OH bands correspond to H atoms in interstitial position; (3) In the frequency region between 3,300 and 3,000 cm^?1, three broad bands are identified at 3,151, 3,178 and 3,217 cm^?1, at ?194 °C. They exhibit significant frequency increase (~20 cm^?1) and broadening (~70 cm^?1) with increasing temperature, indicating moderate H bonding. These bands are compatible with (2H)_Mg defects. A survey of published spectra of forsterite samples synthesized above 5 GPa shows that about 75 % of the incorporated hydrogen belongs to type (1) OH bands associated with Si substitution and 25 % to the broad band at 3,566 cm^?1 (type (2); 3,550 cm^?1 at room temperature). The contribution of OH bands of type (3), associated to (2H)_Mg defects, is negligible. Therefore, solubility of hydrogen in forsterite (and natural olivine compositions) cannot be described by a single solubility law, but by the combination of at least two laws, with different activation volumes and water fugacity exponents.     

Octahedral cation ordering in Mg, Fe2+-olivine: an optical absorption spectroscopic study

Polarized optical absorption spectra of natural olivine, Fa₁₀Fo₉₀, were measured before and after annealing/quenching experiments performed at 650, 800, 1,000 and 1,200°C under controlled oxygen fugacity. It was found that the annealing induces weak but definite changes in the olivine spectra. The intensity of the spin-allowed Z > X-polarized band at 9,560 cm^?1 and shoulder at ～8,300 cm^?1 attributed to Fe²⁺(M2), continuously decreases with annealing temperature, whereas a weaker band at ～11,600 cm^?1 assigned to electronic spin-allowed transitions of Fe²⁺(M1), increases. This evidently shows that annealing treatments cause a redistribution of Fe²⁺ from M2 to M1. The fractionation increases with increasing temperature. This observation is in good correspondence with many diffraction structural studies of natural and synthetic olivines, as well as with recent Raman and Mössbauer investigations by Kolesov and Geiger (Mitt Österr Mineral Ges 149:48, 2004) and Morozov et al. (Eur J Mineral 17:495–500, 2005) evidencing a weak tendency of Fe to order into the M1 site with increasing temperature. However, this deduction is incompatible with the results of the in situ neutron power diffraction study of synthetic FeMgSiO₄ by Redfern et al. (Phys Chem Minerals 27:630–637, 2000). Polarization properties of the UV absorption edge, attributed to ligand-to-metal charge-transfer transitions in Fe³⁺, changes from Y > X ? Z in natural samples to a weak Y ≥ X ≥ Z-pleochroism in annealed ones. This may be due to redistribution of a small content of Fe³⁺ among M1 and M2 structural sites.     

Convergence of iterative coupling for coupled flow and geomechanics

In this paper, we study solving iteratively the coupling of flow and mechanics. We demonstrate the stability and convergence of two widely used schemes: the undrained split method and the fixed stress split method. To our knowledge, this is the first time that such results have been rigorously obtained and published in the scientific literature. In addition, we propose a new stress split method, with faster convergence rate than known schemes. These results are specially important today due to the interest in hydraulic fracturing (Dean and Schmidt SPE J. 14:707–714, 2009; Ji et al. SPE J. 14:423–430, 2009; Samier and De Gennaro 2007; Settari and Maurits SPE J. 3:219–226, 1998), in oil and gas shale reservoirs.     

Comment on “A Comparison of Modified Fuzzy Weights of Evidence,Fuzzy Weights of Evidence,and Logistic Regression for Mapping Mineral Prospectivity” by Daojun Zhang,Frits Agterberg,Qiuming Cheng,and Renguang Zuo Math Geosci DOI 10.1007/s11004-013-9496-8

Despite a missing definition of equivalence of mathematical models or methods by Zhang et al. (Math Geosci, 2013), an “equivalence” (Zhang et al., Math Geosci, 2013, p. 6,7,8,14) of modified weights-of-evidence (Agterberg, Nat Resour Res 20:95–101, 2011) and logistic regression does not generally exist. Its alleged proof is based on a previously conjectured linear relationship between weights of evidence and logistic regression parameters (Deng, Nat Resour Res 18:249–258, 2009), which does not generally exist either (Schaeben and van den Boogaart, Nat Resour Res 20:401–406, 2011). In fact, an extremely simple linear relationship exists only if the predictor variables are conditionally independent given the target variable, in which case the contrasts, i.e., the differences of the weights, are equal to the logistic regression parameters. Thus, weights-of-evidence is the special case of logistic regression if the predictor variables are binary and conditionally independent given the target variable.     

A phase-by-phase upstream scheme that converges to the vanishing capillarity solution for countercurrent two-phase flow in two-rock media

We discuss the convergence of the upstream phase-by-phase scheme (or upstream mobility scheme) towards the vanishing capillarity solution for immiscible incompressible two-phase flows in porous media made of several rock types. Troubles in the convergence were recently pointed out by Mishra and Jaffré (Comput. Geosci. 14, 105–124, 2010) and Tveit and Aavatsmark (Comput. Geosci. 16, 809–825, 2012). In this paper, we clarify the notion of vanishing capillarity solution, stressing the fact that the physically relevant notion of solution differs from the one inferred from the results of Kaasschieter (Comput. Geosci. 3, 23–48, 1999). In particular, we point out that the vanishing capillarity solution depends on the formally neglected capillary pressure curves, as it was recently proven in by Andreianov and Cancès (Comput. Geosci. 17, 551–572, 2013). Then, we propose a numerical procedure based on the hybridization of the interfaces that converges towards the vanishing capillarity solution. Numerical illustrations are provided.     

Assessment of social vulnerability to natural disasters: a comparative study

The purpose of this study is to examine and compare the methodologies being developed in assessing social vulnerability to natural disasters. Existing vulnerability literature shows that two methods have been used in developing social vulnerability indexes: (1) a deductive approach based on a theoretical understanding of relationships and (2) an inductive approach based on statistical relationships (Adger et al. in New indicators of vulnerability and adaptive capacity. Tyndall Centre for Climate Change Research, Norwich, 2004). Two techniques were also utilized in aggregating social vulnerability indicators: (1) a deductive approach using standardization techniques such as z scores or linear scaling (Wu et al. in Clim Res 22:255?C270, 2002; Chakraborty et al. in Nat Hazards Rev 6(1):23?C33, 2005) and (2) an inductive approach using data-reduction techniques such as factor analysis (Clark et al. in Mitig Adapt Strateg Glob Change 3(1):59?C82, 1998; Cutter et al. Soc Sci Quart 84(2):242?C261, 2003). This study empirically compares deductive and inductive index development and indicator aggregation methods in assessing social vulnerability to natural disasters in the Gulf of Mexico and Atlantic coastal areas. The aggregated social vulnerability index is used to examine a relationship with disaster losses in the Gulf of Mexico and Atlantic coastal areas. The results show that coastal counties with more vulnerability in terms of social achieved status are positively associated with disaster damages, while variations in the development of the index using deductive and inductive measurement approaches produce different outcomes.     

Impedance Functions of Slab Foundations with Rigid Piles

The earthquake resistant design of structures depends upon the soil-structure interaction during seismic excitation. The dynamic behavior of surface foundations and deep foundations has been investigated by several authors (Gazetas in Soil Dyn Earthq Eng 2(1):2–42, 1983; Pecker in Dynamique des sols. Presses de l’Ecole Nationale des Ponts et Chaussées, Paris, 1984; Sieffert and Cevaer in Manuel des fonctions d’impédance. Ouest-Editions, Nantes, 1992, etc.) but the dynamic behavior of pile-reinforced soils has not been sufficiently studied yet. The design of pile-reinforced soils comprises the design of the rigid piles and the design of the earth-platform. The foundation system studied in this article consists of an earth-platform over a soft ground reinforced by deep piles. In order to understand the dynamic behavior of a rigid pile, a series of dynamic tests is conducted on an experimental site. The vertical and horizontal impedances of the slab foundation are obtained with and without rigid piles. The numerical models are developed to interpret these dynamic tests. The numerical and experimental impedance functions are compared in both the vertical and the horizontal directions. A sensitivity analysis on the influence of the physical and geometrical properties of rigid piles on the impedance functions is presented.     

The Use of Rock Mass Classification Systems to Estimate the Modulus and Strength of Jointed Rock

Three-dimensional, elastic and elasto-plastic finite element (FE) programs have permitted calculation of the displacements and the factor of safety (FOS) for the excavation for a tower, 132.70 m high (above foundation) on the island of Tenerife. The tower is supported by a 2 m thick reinforced concrete slab on jointed, vesicular and weathered basalt and scoria. The installation of rod extensometers at different depths below the slab has permitted comparison between measured and calculated displacements and the estimation of in situ deformation modulus. The moduli deduced from the simple empirical equations proposed by Hoek et al. (In: NARMS-TAC, 2002) and Gokceoglu et al. (Int J Rock Mech Min Sci 40:701–710, 2003) as a function of GSI, and Nicholson and Bieniawski (Int J Min Geol Eng 8:181–202, 1990) as a function of RMR, provide an acceptable fit with the measured settlements in this type of rock. Good correlation is also obtained with the empirical equation presented by Verman et al. (Rock Mech Rock Eng 30(3):121–127, 1997) that incorporates the influence of confining stress in the deformation modulus. The FOS obtained from different correlations with geomechanical classifications is within a relatively narrow range. These results increase our confidence in the use of classification schemes to estimate the deformation and stability in jointed rock.     

Infrared spectroscopic properties of goethite: anharmonic broadening,long-range electrostatic effects and Al substitution

The infrared spectrum and its temperature dependence (20–320 K) were collected on a synthetic goethite sample (α-FeOOH). In addition, the infrared powder absorption spectrum of goethite and aluminum-substituted goethite was computed using ab initio quantum mechanical calculations based on density functional theory. This combined experimental and theoretical approach allows (1) the unequivocal assignment of absorption bands to the corresponding vibrational modes, (2) separate identification of the effects of the particle shape and of the aluminum substitution on the infrared spectrum, and (3) a discussion of the anharmonic properties and the origin of the line broadening in goethite. In particular, the two well-resolved OH bending absorption bands show different temperature evolution. Their detailed analysis suggests that the broadening of the band at ~800 cm^?1 cannot be described solely by a usual three-phonon process. The strong anharmonic behavior of this mode implies the addition of a four-phonon process, such as a pure dephasing process. In our calculations, the effect of the Hubbard U correction is also investigated and found to be most visible on the OH stretching and bending modes, in relation to the associated structural relaxation. The OH stretching frequencies decrease, leading to a better agreement with experimental frequencies, while the OH bending frequencies increase.     

Investigation of the variability of strong ground motions from Vrancea earthquakes

A systematic investigation of the applicability of several ground motion prediction models for Vrancea intermediate-depth seismic source is conducted in this research. Two ground motion prediction models recommended by previous evaluations (Vacareanu et al. in Bull Earthq Eng 11(6):1867–1884, 2013a; Pavel et al. in Earthq Struct 6(1):1–18, 2014), as well as two new state-of-the-art ground motion prediction equations (Vacareanu et al. in J Earthq Eng, 2013b; Earthq Struct 6(2):141–161, 2014) are tested using an increased strong ground motion database consisting of 150 recordings from Vrancea subcrustal earthquakes. The evaluation is performed by using several goodness-of-fit parameters from the literature. Moreover, the applicability of the single-station sigma method is also investigated by using the same strong ground motion database recorded in 30 seismic stations from southern and eastern Romania. The influence of the soil conditions on the numerical results obtained in this study is investigated and discussed using the results provided by the analysis of variance method. The impact of the single-station standard deviation on the levels of seismic hazard is also assessed in this study, and the results show, in the analyzed cases, significant reductions of the hazard levels.     

Evaluating the Poroelastic Effect on Anisotropic, Organic-Rich, Mudstone Systems

Understanding the poroelastic effect on anisotropic organic-rich mudstones is of high interest and value for evaluating coupled effects of rock deformation and pore pressure, during drilling, completion and production operations in the oilfield. These applications include modeling and prevention of time-dependent wellbore failure, improved predictions of fracture initiation during hydraulic fracturing operations (Suarez-Rivera et al. Presented at the Canadian Unconventional Resources Conference held in Calgary, Alberta, Canada, 15–17 November 2011. CSUG/SPE 146998 2011), improved understanding of the evolution of pore pressure during basin development, including subsidence and uplift, and the equilibrated effective in situ stress (Charlez, Rock mechanics, vol 2 1997; Katahara and Corrigan, Pressure regimes in sedimentary basins and their prediction: AAPG Memoir, vol 76, pp 73–78 2002; Fjær et al. Petroleum related rock mechanics. 2nd edn 2008). In isotropic rocks, the coupled poro-elastic deformations of the solid framework and the pore fluids are controlled by the Biot and Skempton coefficients. These are the two fundamental properties that relate the rock framework and fluid compressibility and define the magnitude of the poroelastic effect. In transversely isotropic rocks, one desires to understand the variability of these coefficients along the directions parallel and longitudinal to the principal directions of material symmetry (usually the direction of bedding). These types of measurements are complex and uncommon in low-porosity rocks, and particularly problematic and scarce in tight shales. In this paper, we discuss a methodology for evaluating the Biot’s coefficient, its variability along the directions parallel and perpendicular to bedding as a function of stress, and the homogenized Skempton coefficient, also as a function of stress. We also predict the pore pressure change that results during undrained compression. Most importantly, we provide values of transverse and longitudinal Biot’s coefficients and the homogenized Skempton coefficient for two important North American, gas-producing, organic-rich mudstones. These results could be used for petroleum-related applications.     

Crystal reaming during the assembly, maturation, and waning of an eleven-million-year crustal magma cycle: thermobarometry of the Aucanquilcha Volcanic Cluster

Phenocryst assemblages of lavas from the long-lived Aucanquilcha Volcanic Cluster (AVC) have been probed to assess pressure and temperature conditions of pre-eruptive arc magmas. Andesite to dacite lavas of the AVC erupted throughout an 11-million-year, arc magmatic cycle in the central Andes in northern Chile. Phases targeted for thermobarometry include amphibole, plagioclase, pyroxenes, and Fe–Ti oxides. Overall, crystallization is documented over 1–7.5 kbar (~25 km) of pressure and ~680–1,110 °C of temperature. Pressure estimates range from ~1 to 5 kbar for amphiboles and from ~3 to 7.5 kbar for pyroxenes. Pyroxene temperatures are tightly clustered from ~1,000–1,100 °C, Fe–Ti oxide temperatures range from ~750–1,000 °C, and amphibole temperatures range from ~780–1,050 °C. Although slightly higher, these temperatures correspond well with previously published zircon temperatures ranging from ~670–900 °C. Two different Fe–Ti oxide thermometers (Andersen and Lindsley 1985; Ghiorso and Evans 2008) are compared and agree well. We also compare amphibole and amphibole–plagioclase thermobarometers (Ridolfi et al. 2010; Holland and Blundy 1994; Anderson and Smith 1995), the solutions from which do not agree well. In samples where we employ multiple thermometers, pyroxene temperature estimates are always highest, zircon temperature estimates are lowest, and Fe–Ti oxide and amphibole temperature estimates fall in between. Maximum Fe–Ti oxide and zircon temperatures are observed during the middle stage of AVC activity (~5–3 Ma), a time associated with increased eruption rates. Amphibole temperatures during this time are relatively restricted (~850–1,000 °C). The crystal record presented here offers a time-transgressive view of an evolving, multi-tiered subvolcanic reservoir. Some crystals in AVC lavas are likely to be true phenocrysts, but the diversity of crystallization temperatures and pressures recorded by phases in individual AVC lavas suggests erupting magma extensively reams and accumulates crystals from disparate levels of the middle to upper crust.