Measurements of a water potential and water content in unsaturated crystalline rock

A water desaturation zone develops around a tunnel in water-saturated rock when the evaporative water loss at the rock surface is larger than the water flow from the surrounding saturated region of restricted permeability. We describe the methods with which such water desaturation processes in rock materials can be quantified. The water retention characteristic theta(psi) of crystalline rock samples was determined with a pressure membrane apparatus. The negative water potential, identical to the capillary pressure, psi, below the tensiometric range (psi < -0.1 MPa) can be measured with thermocouple psychrometers (TP), and the volumetric water contents, theta, by means of time domain reflectometry (TDR). These standard methods were adapted for measuring the water status in a macroscopically unfissured granodiorite with a total porosity of approximately 0.01. The measured water retention curve of granodiorite samples from the Grimsel test site (central Switzerland) exhibits a shape which is typical for bimodal pore size distributions. The measured bimodality is probably an artifact of a large surface ratio of solid/voids. The thermocouples were installed without a metallic screen using the cavity drilled into the granodiorite as a measuring chamber. The water potentials observed in a cylindrical granodiorite monolith ranged between -0.1 and -3.0 MPa; those near the wall in a ventilated tunnel between -0.1 and -2.2 MPa. Two types of three-rod TDR Probes were used, one as a depth probe inserted into the rock, the other as a surface probe using three copper stripes attached to the surface for detecting water content changes in the rock-to-air boundary. The TDR signal was smoothed with a low-pass filter, and the signal length determined based on the first derivative of the trace. Despite the low porosity of crystalline rock these standard methods are applicable to describe the unsaturated zone in solid rock and may also be used in other consolidated materials such as concrete.

Matrix pore water in low-permeable crystalline bedrock: An archive for the palaeohydrogeological evolution of the Olkiluoto Investigation Site

Matrix pore water in the connected inter- and intragranular pore space of low-permeable crystalline bedrock interacts with flowing fracture groundwater predominately by diffusion. Based on the slow exchange between the two water reservoirs, matrix pore water acts as an archive of past changes in fracture groundwater compositions and thus of the palaeohydrological history of a site. Matrix pore water of crystalline bedrock from the Olkiluoto investigation site (SW Finland) was characterised using the stable water isotopes (δ18O, δ2H), combined with the concentrations of dissolved chloride and bromide as natural tracers. The comparison of tracer concentrations in pore water and present-day fracture groundwater suggest for the pore water the presence of old, dilute meteoric water components that infiltrated into the fractures during various warm climate stages. These different meteoric components can be discerned based on the diffusion distance between the two reservoirs and be brought into context with the palaeohydrological evolution of the site.

Assessing Initial Conditions for Chloride Transport Across Low- permeability Argillaceous Rocks, Wellenberg, Switzerland

Information about fluid evolution and solute transport in a low-permeability metamorphic rock sequence has been obtained by comparing chloride concentrations and chlorine isotope ratios of pore water, groundwater, and fluid inclusions. The similarity of d37Cl values in fluid inclusions and groundwater suggests a closed-system evolution during the metamorphic overprint, and signatures established at this time appear to form the initial conditions for chloride transport after exhumation of the rock sequence.

Subcritical Propagation and Coalescence of Oil-Filled Cracks: Getting the Oil Out of Low-Permeability Source Rocks

We use a fracture mechanics model to study subcritical propagation and coalescence of single and collinear oil-filled cracks during conversion of kerogen to oil. The subcritical propagation distance, propagation duration, crack coalescence and excess oil pressure in the crack are determined using the fracture mechanics model together with the kinetics of kerogen-oil transformation. The propagation duration for the single crack is governed by the transformation kinetics whereas the propagation duration for the multiple collinear cracks may vary by two orders of magnitude depending on initial crack spacing. A large amount of kerogen (>90%) remains unconverted when the collinear cracks coalesce and the new, larger cracks resulting from coalescence will continue to propagate with continued kerogen-oil conversion. The excess oil pressure on the crack surfaces drops precipitously when the collinear cracks are about to coalesce, and crack propagation duration and oil pressure on the crack surfaces are strongly dependent on temperature. Citation: Jin, Z.-H., S. E. Johnson, and Z. Q. Fan (2010), Subcritical propagation and coalescence of oil-filled cracks: Getting the oil out of low-permeability source rocks, Geophys. Res. Lett., 37, L01305, doi:10.1029/2009GL041576.

Metamorphism, Transient Mid-Crustal Rheology, Strain Localization and the Exhumation of High-Grade Metamorphic Rocks

We present a series of three-dimensional numerical models investigating the effects of metamorphic strengthening and weakening on the geodynamic evolution of convergent orogens that are constrained by observations from an exposed mid-crustal section in the New England Appalachians. The natural mid-crustal section records evidence for spatially and temporally variable mid-crustal strength as a function of metamorphic grade during prograde polymetamorphism. Our models address changes in strain rate partitioning and topographic uplift as a function of strengthening/weakening in the middle crust, as well as the resultant changes in deformation kinematics and potential exhumation patterns of high-grade metamorphic rock. Results suggest that strengthening leads to strain rate partitioning around the zone and suppressed topographic uplift rates whereas weakening leads to strain rate partitioning into the zone and enhanced topographic uplift rates. Deformation kinematics recorded in the orogen are also affected by strengthening/weakening, with complete reversals in shear sense occurring as a function of strengthening/weakening without changes in plate boundary kinematics.

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.

Chevkinite-Group Minerals from Granulite-Facies Metamorphic Rocks and Associated Pegmatites of East Antarctica and South India

Electron microprobe data are presented for chevkinite-group minerals from granulite-facies rocks and associated pegmatities of the Napier Complex and Mawson Station charnockite in East Antarctica and from the Eastern Ghats, South India. Their compositions conform to the general formula for this group, viz. A(4)BC(2)D(2)Si(4)O(22) where, in the analysed specimens A = (rare-earth elements (REE), Ca, Y, Th), B = Fe(2+) Mg, C = (Al, Mg, Ti, Fe(2+), Fe(3+), Zr) and D = Ti and plot within the perrierite field oftlic total Fe (as FeO) (wt.%) vs. CaO (wt.%) discriminator diagram of Macdonald and Belkin (2002). In contrast to most chevkinite-group minerals, the A site shows unusual enrichment in the MREE and HREE relative to the LREE and Ca. In one sample from the Napier Complex, Y is the dominant cation among the total REE + Y in the A site, the first reported case of Y-dominance in the chevkinite group. The minerals include the most Al-rich yet reported in the chevkinite group (<= 9.15 wt.% Al(2)O(3)), sufficient to fill the C site in two samples. Conversely, the amount of Ti in these samples does not fill the D site. and, thus, some of the Al could be making up the deficiency at D, a situation not previously reported in the chevkinite group. Fe abudances are low, requiring Mg to occupy up to 45% of the B site. The chevkinite-group minerals analysed originated from three distinct parageneses: (1) pegmatites containing hornblende and orthopyroxene or garnet; (2) orthopyroxene-bearing gneiss and granulite; (3) highly aluminous paragneisses in which the associated minerals are relatively magnesian or aluminous. Chevkinite-group minerals from the first two parageneses have relatively high FeO content and low MgO and Al(2)O(3) contents; their compositions plot in the field for mafic and intermediate igneous rocks. In contrast, chevkinite-group minerals from the third paragenesis are notably more aluminous and have greater Mg/Fe ratios.

Crystalline Confinement

We show that exotic phases arise in generalized lattice gauge theories known as quantum link models in which classical gauge fields are replaced by quantum operators. While these quantum models with discrete variables have a finite-dimensional Hilbert space per link, the continuous gauge symmetry is still exact. An efficient cluster algorithm is used to study these exotic phases. The (2+1)-d system is confining at zero temperature with a spontaneously broken translation symmetry. A crystalline phase exhibits confinement via multi stranded strings between chargeanti-charge pairs. A phase transition between two distinct confined phases is weakly first order and has an emergent spontaneously broken approximate SO(2) global symmetry. The low-energy physics is described by a (2 + 1)-d RP(1) effective field theory, perturbed by a dangerously irrelevant SO(2) breaking operator, which prevents the interpretation of the emergent pseudo-Goldstone boson as a dual photon. This model is an ideal candidate to be implemented in quantum simulators to study phenomena that are not accessible using Monte Carlo simulations such as the real-time evolution of the confining string and the real-time dynamics of the pseudo-Goldstone boson.

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.

Crystalline, Mixed-Valence Manganese Analogue of Prussian Blue: Magnetic, Spectroscopic, X-ray and Neutron Diffraction Studies

The compound of stoichiometry Mn(II)3[Mn(III)(CN)6]2·zH2O (z = 12−16) (1) forms air-stable, transparent red crystals. Low-temperature single crystal optical spectroscopy and single crystal X-ray diffraction provide compelling evidence for N-bonded high-spin manganese(II), and C-bonded low-spin manganese(III) ions arranged in a disordered, face-centered cubic lattice analogous to that of Prussian Blue. X-ray and neutron diffraction show structured diffuse scattering indicative of partially correlated (rather than random) substitutions of [Mn(III)(CN)6] ions by (H2O)6 clusters. Magnetic susceptibility measurements and elastic neutron scattering experiments indicate a ferrimagnetic structure below the critical temperature Tc = 35.5 K.

On the Methyl-Transfer Reaction in Crystalline Methyl 2-(Methylthio)benzenesulfonate: a Thermally Induced Non-Topochemical Solid-State Reaction

The rearrangement of methyl 2-(methylthio)benzenesulfonate (1) to the zwitterionic 2-(dimethyl-sulfonium)benzenesulfonate (2) is known to proceed in solution by intermolecular Me transfers. The same rearrangement has been observed to occur in crystalline 1, but the crystal structure shows that the molecular packing is not conducive to intermolecular Me transfer. The reaction has been carried out with mixed crystals composed of 1 and deuteriomethylated (D6)-l. By fast-atom-bombardment mass spectroscopy, it has been shown that the product consists of a 1:2:1 mixture of the non-, tri-, and hexadeuterated species, the mixture expected, if the solid-state reaction proceeds by intermolecular Me transfers. From this result, together with the slower rates of conversion in the single crystal compared with the melt, it can be concluded that the reaction must occur not topochemically but rather at defects such as microcavities, surfaces, and other irregularities in the ordered crystal arrangement.

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.