Magma-Driven Multiple Dike Propagation and Fracture Toughness of Crustal Rocks

Dike swarms consisting of tens to thousands of subparallel dikes are commonly observed at Earth's surface, raising the possibility of simultaneous propagation of two or more dikes at various stages of a swarm's development. The behavior of multiple propagating dikes differs from that of a single dike owing to the interacting stress fields associated with each dike. We analyze an array of parallel, periodically spaced dikes that grow simultaneously from an overpressured source into a semi-infinite, linear elastic host rock. To simplify the analysis, we assume steady state (constant velocity) magma flow and dike propagation. We use a perturbation method to analyze the coupled, nonlinear problem of multiple dike propagation and magma transport. The stress intensity factor at the dike tips and the opening displacements of the dike surfaces are calculated. The numerical results show that dike spacing has a profound effect on the behavior of dike propagation. The stress intensity factors at the tips of parallel dikes decrease with a decrease in dike spacing and are significantly smaller than that for a single dike with the same length. The reduced stress intensity factor indicates that, compared to a single dike, propagation of parallel dikes is more likely to be arrested under otherwise the same conditions. It also implies that fracture toughness of the host rock in a high confining pressure environment may not be as high as inferred from the propagation of a single dike. Our numerical results suggest fracture toughness values on the order of 100 MPa root m. The opening displacements for parallel dikes are smaller than that for a single dike, which results in higher magma pressure gradients in parallel dikes and lower flux of magma transport.

Timing of Tertiary Extension in the Railroad Valley Pioche Transect, Nevada: Constraints from Ar-40/Ar-39 Ages of Volcanic Rocks

Time-space relations of extension and volcanism place critical constraints on models of Basin and Range extensional processes. This paper addresses such relations in a 130-km-wide transect in the eastern Great Basin, bounded on the east by the Ely Springs Range and on the west by the Grant and Quinn Canyon ranges. Stratigraphic and structural data, combined with 40Ar/39Ar isotopic ages of volcanic rocks, document a protracted but distinctly episodic extensional history. Field relations indicate four periods of faulting. Only one of these periods was synchronous with nearby volcanic activity, which implies that volcanism and faulting need not be associated closely in space and time. Based on published dates and the analyses reported here, the periods of extension were (1) prevolcanic (pre-32 Ma), (2) early synvolcanic (30 to 27 Ma), (3) immediately postvolcanic (about 16 to 14 Ma), and (4) Pliocene to Quaternary. The break between the second and third periods is distinct. The minimum gap between the first two periods is 2 Ma, but the separation may be much larger. Temporal separation of the last two periods is only suggested by the stratigraphic record and cannot be rigorously demonstrated with present data. The three younger periods of faulting apparently occurred across the entire transect. The oldest period is recognized only at the eastern end of the transect, but appears to correlate about 150 km northward along strike with extension in the Northern Snake Range-Kern Mountains area. Therefore the oldest period also is regional in extent, but affected a different area than that affected by younger periods. This relation suggests that distinct extensional structures and master detachment faults were active at different times. The correlation of deformation periods of a few million years duration across the Railroad Valley-Pioche transect suggests that the scale of active extensional domains in the Great Basin may be greater than 100 km across strike.

Revisiting the Nd-142 deficits in the 1.48 Ga Khariar alkaline rocks, India

The 146Sm–142Nd system plays a central role in tracing the silicate differentiation of the Earth prior to 4.1 Ga. After this time, given its initial abundance, the 146Sm can be considered to be extinct. Upadhyay et al. (2009) reported unexpected negative 142Nd anomalies in 1.48 Ga rocks of the Khariar nepheline syenite complex (India) and inferred that an early enriched, low-Sm/Nd reservoir must have contributed to the mantle source rocks of the Khariar complex. As 146Sm had been effectively extinct for about 2.6 billion years before the crystallisation of the Khariar samples, this Nd signature should have remained isolated from the convective mantle for at least that long. It was thus suggested that the source rock of Khariar samples had been sequestered in the lithospheric root of the Indian craton. Using a different chemical separation method, and a different Thermal Ionization Mass Spectrometry (TIMS) analysis protocol, the present study attempted to replicate these negative 142Nd anomalies, but none were found. To determine which data set is correct, we investigated three possible sources of bias between them: imperfect cancellation of Faraday collector efficiencies during multidynamic TIMS analysis, rapid sample fractionation between the sequential measurement of 146Nd/144Nd and 142Nd/144Nd, and non-exponential law behaviour resulting from so-called “domain mixing.” Incomplete cancellation of collector efficiencies was found unlikely to cause resolvable biases at the estimated level of variation among collector efficiencies. Even in the case of highly variable efficiency and resolvable biases, there is no reason to suspect that they would reproducibly affect only four rocks out of 10 analysed by Upadhyay et al. (2009). Although domain mixing may explain apparent “reverse” fractionation trends observed in some TIMS analyses, it cannot be the cause of the apparent negative anomalies in the study of Upadhyay et al. (2009). It was determined that rapid mass fractionation during the course of a multidynamic TIMS analysis can bias all measured Nd ratios. After applying an approximate correction for this effect, only one rock from Upadhyay et al. (2009) retained an apparent negative 142Nd anomaly. This, in conjunction with our new, anomaly-free data set measured at fractionation rates too low to cause bias, leads to the conclusion that the anomalies reported by Upadhyay et al. (2009) are a subtle and reproducible analytical artefact. The absence of negative 142Nd anomalies in these rocks relaxes the need for a mechanism (other than crust formation) that can isolate a Nd reservoir from the convective mantle for billions of years.

Precambrian to Paleozoic zircon record in the Siviez-Mischabel basement (western Swiss Alps)

U–Pb zircon analyses from three meta-igneous and two metasedimentary rocks from the Siviez-Mischabel nappe in the western Swiss Alps are presented, and are used to derive an evolutionary history spanning from Paleoarchean crustal growth to Permian magmatism. The oldest components are preserved in zircons from metasedimentary albitic schists. The oldest zircon core in these schists is 3.4 Ga old. Detrital zircons reveal episodes of crustal growth in the Neoarchean (2.7–2.5 Ga), Paleoproterozoic (2.2–1.9 Ma) and Neoproterozoic (800–550 Ma, Pan-African event). The maximum age of deposition for the metasedimentary rocks is given by the youngest detrital zircons within both metasedimentary samples dated at ~490 Ma (Cambrian-Ordovician boundary). This is in the age range of two granitoid samples dated at 505 ± 4 and 482 ± 7 Ma, and indicates sedimentation and magmatism in an extensional setting preceding an Ordovician orogeny. The third felsic meta-igneous rock gives a Permian age of intrusion, and is part of a long-lasting Variscan to post-Variscan magmatic activity. The zircons record only minor disturbance of the U–Pb system during the Alpine orogeny.

Matrix Porewater In Crystalline Rocks: Extraction and Analysis

The chemical and isotopic characterization of porewater residing in the inter- and intragranular pore space of the low-permeability rock matrix is an important component with respect to the site characterization and safety assessment of potential host rocks for a radioactive waste disposal. The chemical and isotopic composition of porewater in such low permeability rocks has to be derived by indirect extraction techniques applied to naturally saturated rock material. In most of such indirect extraction techniques – especially in case of rocks of a porosity below about 2 vol.% – the original porewater concentrations are diluted and need to be back-calculated to in-situ concentrations. This requires a well-defined value for the connected porosity – accessible to different solutes under in-situ conditions. The derivation of such porosity values, as well as solute concentrations, is subject to various perturbations during drilling, core sampling, storage and experiments in the laboratory. The present study aims to demonstrate the feasibility of a variety of these techniques to charac-terize porewater and solute transport in crystalline rocks. The methods, which have been de-veloped during multiple porewater studies in crystalline environments, were applied on four core samples from the deep borehole DH-GAP04, drilled in the Kangerlussuaq area, Southwest Greenland, as part of the joint NWMO–Posiva–SKB Greenland Analogue Project (GAP). Potential artefacts that can influence the estimation of in situ porewater chemistry and isotopes, as well as their controls, are described in detail in this report, using specific examples from borehole DH-GAP04

Magnetic anisotropy in clinopyroxene and orthopyroxene single crystals

Pyroxenes constitute an important component in mafic igneous and metamorphic rocks. They often possess a prismatic habit, and their long axis, the crystallographic c axis, helps define a lineation in a textured rock. Anisotropy of magnetic susceptibility (AMS) serves as a fabric indicator in igneous and metamorphic rocks. If a rock’s AMS is carried by pyroxenes, it can be related to their crystallographic preferred orientation and degree of alignment. This requires knowing the intrinsic AMS of pyroxene single crystals. This study provides a comprehensive low-field and high-field AMS investigation of chemically diverse orthopyroxene and clinopyroxene crystals in relation to crystal structure, chemical composition, oxidation state of Fe, and the possible presence of ferromagnetic inclusions. The paramagnetic anisotropy, extracted from high-field data, shows clear relationships to crystallographic directions and Fe concentration both in clinopyroxene and orthopyroxene. In the diopside-augite series, the intermediate susceptibility is parallel to b, and the maximum is at 45° to the c axis. In aegirine, the intermediate axis remains parallel to b, while the maximum susceptibility is parallel to c. The AMS of spodumene depends on Fe concentration. In enstatite, the maximum susceptibility aligns with c and the minimum with b, and in the case of hypersthene, the maximum susceptibility is normal to the exsolution lamellae. Magnetite inclusions within augite possess a ferromagnetic anisotropy with consistent orientation of the principal susceptibilities, which dominates the low-field anisotropy. These results provide better understanding of magnetic anisotropy in pyroxenes and form a solid basis for interpretation of magnetic fabrics in pyroxene-bearing rocks.

Predictive modelling of mineral scaling, corrosion and the performance of solute geothermometers in a granitoid-hosted, enhanced geothermal system

Experience is lacking with mineral scaling and corrosion in enhanced geothermal systems (EGS) in which surface water is circulated through hydraulically stimulated crystalline rocks. As an aid in designing EGS projects we have conducted multicomponent reactive-transport simulations to predict the likely characteristics of scales and corrosion that may form when exploiting heat from granitoid reservoir rocks at ∼200 °C and 5 km depth. The specifications of an EGS project at Basel, Switzerland, are used to constrain the model. The main water–rock reactions in the reservoir during hydraulic stimulation and the subsequent doublet operation were identified in a separate paper (Alt-Epping et al., 2013b). Here we use the computed composition of the reservoir fluid to (1) predict mineral scaling in the injection and production wells, (2) evaluate methods of chemical geothermometry and (3) identify geochemical indicators of incipient corrosion. The envisaged heat extraction scheme ensures that even if the reservoir fluid is in equilibrium with quartz, cooling of the fluid will not induce saturation with respect to amorphous silica, thus eliminating the risk of silica scaling. However, the ascending fluid attains saturation with respect to crystalline aluminosilicates such as albite, microcline and chlorite, and possibly with respect to amorphous aluminosilicates. If no silica-bearing minerals precipitate upon ascent, reservoir temperatures can be predicted by classical formulations of silica geothermometry. In contrast, Na/K concentration ratios in the production fluid reflect steady-state conditions in the reservoir rather than albite–microcline equilibrium. Thus, even though igneous orthoclase is abundant in the reservoir and albite precipitates as a secondary phase, Na/K geothermometers fail to yield accurate temperatures. Anhydrite, which is present in fractures in the Basel reservoir, is predicted to dissolve during operation. This may lead to precipitation of pyrite and, at high exposure of anhydrite to the circulating fluid, of hematite scaling in the geothermal installation. In general, incipient corrosion of the casing can be detected at the production wellhead through an increase in H2(aq) and the enhanced precipitation of Fe-bearing aluminosilicates. The appearance of magnetite in scales indicates high corrosion rates.

P-T-t estimation of syn-kinematic strain in low-grade quartz-feldspar bearing rocks using thermodynamic modeling and 40Ar/39Ar dating techniques: example of the Plan-de-Phasy shear zone unit (Briançonnais Zone, Western Alps)

Pressure–Temperature–time (P–T–t) estimates of the syn-kinematic strain at the peak-pressure conditions reached during shallow underthrusting of the Briançonnais Zone in the Alpine subduction zone was made by thermodynamic modelling and 40Ar/39Ar dating in the Plan-de-Phasy unit (SE of the Pelvoux Massif, Western Alps). The dated phengite minerals crystallized syn-kinematically in a shear zone indicating top-to-the-N motion. By combining X-ray mapping with multi-equilibrium calculations, we estimate the phengite crystallization conditions at 270 ± 50 °C and 8.1 ± 2 kbar at an age of 45.9 ± 1.1 Ma. Combining this P–T–t estimate with data from the literature allows us to constrain the timing and geometry of Alpine continental subduction. We propose that the Briançonnais units were scalped on top of the slab during ongoing continental subduction and exhumed continuously until collision.

Microscale mapping of alteration conditions and potential biosignatures in basaltic-ultramafic rocks on early Earth and beyond

Subseafloor environments preserved in Archean greenstone belts provide an analogue for investigating potential subsurface habitats on Mars. The c. 3.5-3.4 Ga pillow lava metabasalts of the mid-Archean Barberton greenstone belt, South Africa, have been argued to contain the earliest evidence for microbial subseafloor life. This includes candidate trace fossils in the form of titanite microtextures, and sulfur isotopic signatures of pyrite preserved in metabasaltic glass of the c. 3.472 Ga Hooggenoeg Formation. It has been contended that similar microtextures in altered martian basalts may represent potential extraterrestrial biosignatures of microbe-fluid-rock interaction. But despite numerous studies describing these putative early traces of life, a detailed metamorphic characterization of the microtextures and their host alteration conditions in the ancient pillow lava metabasites is lacking. Here, we present a new nondestructive technique with which to study the in situ metamorphic alteration conditions associated with potential biosignatures in mafic-ultramafic rocks of the Hooggenoeg Formation. Our approach combines quantitative microscale compositional mapping by electron microprobe with inverse thermodynamic modeling to derive low-temperature chlorite crystallization conditions. We found that the titanite microtextures formed under subgreenschist to greenschist facies conditions. Two chlorite temperature groups were identified in the maps surrounding the titanite microtextures and record peak metamorphic conditions at 315 ± 40°C (XFe3+(chlorite) = 25-34%) and lower-temperature chlorite veins/microdomains at T = 210 ± 40°C (lower XFe3+(chlorite) = 40-45%). These results provide the first metamorphic constraints in textural context on the Barberton titanite microtextures and thereby improve our understanding of the local preservation conditions of these potential biosignatures. We suggest that this approach may prove to be an important tool in future studies to assess the biogenicity of these earliest candidate traces of life on Earth. Furthermore, we propose that this mapping approach could also be used to investigate altered mafic-ultramafic extraterrestrial samples containing candidate biosignatures.

Vapnikite Ca3UO6 – a new double-perovskite mineral from pyrometamorphic larnite rocks of the Jabel Harmun, Palestinian Autonomy, Israel

The new mineral species vapnikite, Ca3UO6, was found in larnite pyrometamorphic rocks of the Hatrurim Formation at Jabel Harmun in the Judean desert, Palestinian Autonomy, Israel. Vapnikite is an analogue of the synthetic ordered double-perovskite β-Ca3UO6 and is isostructural with the natural fluorperovskite – cryolite Na3AlF6. Vapnikite Ca3UO6 (P21/n,Z = 2, a = 5.739(1), b = 5.951(1), c = 8.312(1) Å, β = 90.4(1)°, V = 283.9(1) Å3) forms yellow-brown xenomorphic grains with a strong vitreous lustre. Small grains up to 20 – 30 m m in size are wedged between larnite, brownmillerite and ye'elimite. Vapnikite has irregular fracture, cleavage and parting were not observed. The calculated density is 5.322 g cm–3, the microhardness is VHN25 = 534 kg mm–2 (mean of seven measurements) corresponding to the hardness of ~5 on the Mohs scale. The crystal structure of vapnikite Ca3UO6 differs from that of its synthetic analogue β-Ca3UO6 by having a larger degree of Ca, U disorder. Vapnikite formed at the high-temperature retrograde stage of pyrometamorphism when larnite rocks were altered by fluids/melts of high alkalinity.

Granular flow in polymineralic rocks bearing sheet silicates: New evidence from natural examples

Natural deformation in carbonate mylonites bearing sheet silicates occurs via a complex interaction of granular flow and solution transfer processes and involves continuous cycles of dissolution, grain boundary diffusion, nucleation and growth. In this way, new sheet silicates (a) nucleate within voids formed by grain boundary sliding of calcite grains. (b) grow, and (c) rotate towards the shear plane. As a consequence, small mica grains show a wide range of orientations with respect to the shear plane, but moderate to large grains are subparallel both to each other and to the shear plane. Increases of average grain sizes with increasing temperature of sheet silicates in mica-rich layers is more pronounced than in mica-poor layers. In the calcitic matrix however, sheet silicates can only grow via solution-precipitation and mass transfer processes. Therefore, the observed grain size variability indicates drastic differences in mass transfer behavior between the individual layers, which might be related to differences in the fluid flux. Based on these observations, a conceptual model for the microfabric evolution in sheet silicate bearing mylonites is presented. © 2001 Elsevier Science B.V. All rights reserved.

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

39Ar-40Ar geochronology of mono- and polymetamorphic basement rocks

The 39Ar-40Ar technique is often used to date the metamorphic evolution of basement rocks. The present review article examines systematic aspects of the K-Ar decay system in different mineral chronometers frequently found in mono- and polymetamorphic basements (amphibole, biotite, muscovite/phengite, K-feldspar). A key observation is that the measured dissolution rate of silicates in aqueous fluids is many orders of magnitude faster, and has a much lower activation energy, than the rate of Fickian diffusion of Ar. The effects of this inequality are patchy age zonations, very much like those observed in many U-Pb chronometers, unaccompanied by intra-crystalline bell¬shaped Ar loss profiles. Recognizing the importance of the respective rate constants in field situations leads to re-evaluating the ages and the interpretive paradigms in classic examples such as the Central Alpine "Lepontine" amphibolite event and the Western Alpine eclogitic event.