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.

Sedimentary Metal Bioavailability Determined By the Digestive Constraints of Marine Deposit Feeders: Gut Retention Time and Dissolved Amino Acids

Contaminant metals bound to sediments are subject to considerable solubilization during passage of the sediments through the digestive systems of deposit feeders. We examined the kinetics of this process, using digestive fluids extracted from deposit feeders Arenicola marina and Parastichopus californicus and then incubated with contaminated sediments. Kinetics are complex, with solubilization followed occasionally by readsorption onto the sediment. In general, solubilization kinetics are biphasic, with an initial rapid step followed by a slower reaction. For many sediment-organism combinations, the reaction will not reach a steady state or equilibrium within the gut retention time (GRT) of the organisms, suggesting that metal bioavailability in sediments is a time-dependent parameter. Experiments with commercial protein solutions mimic the kinetic patterns observed with digestive fluids, which corroborates our previous study that complexation by dissolved amino acids (AA) in digestive fluids leads to metal solubilization (Chen & Mayer 1998b; Environ Sci Technol 32:770-778). The relative importance of the fast and slow reactions appears to depend on the ratio of ligands in gut fluids to the amount of bound metal in sediments. High ligand to solid metal ratios result in more metals released in fast reactions and thus higher lability of sedimentary metals. Multiple extractions of a sediment with digestive fluid of A. marina confirm the potential importance of incomplete reactions within a single deposit-feeding event, and make clear that bioavailability to a single animal is Likely different from that to a community of organisms. The complex kinetic patterns lead to the counterintuitive prediction that toxification of digestive enzymes by solubilized metals will occur more readily in species that dissolve less metals.

The Yanque prospect (Peru): from polymetallic Zn-Pb mineralization to a nonsulfide deposit

The Yanque nonsulfide Pb-Zn deposit (inferred resources 12.5 Mt @ 3.7% Pb and @ 3.5% Zn) is located in the Andahuaylas-Yauri ore province (Cuzco, southern Peru). The deposit occurs within a base metal mineralized district, centered on the medium-sized Dolores porphyry copper. A thorough geological, mineralogical and geochemical study has carried out in order to define: the relationships between the Dolores Cu-porphyry ore and the Yanque Zn-Pb polymetallic mineralization, and the characteristics of the economic nonsulfide concentrations. Both sedimentary and igneous rocks constitute the backbone of the Yanque-Dolores area. The sedimentary lithologies belong to the Soraya, Mara and Ferrobamba Fms. (upper Jurassic-middle Cretaceous). The Yanque orebody is hosted by the Mara Fm., which prevailingly consists of a siliciclastic sedimentary breccia. The original sulfide mineralization consisted of galena, pyrite and sphalerite. The host rock has been affected by a strong hydrothermal alteration, characterized by prevailing sericite/illite, as in the typical porphyry-related phyllic-argillic alteration stage, and by minor kaolinite, dolomite and quartz. Minor element geochemistry, characterized by Sb, As, Mn, Ag and locally also by Cu, points to magmatic-hydrothermal related mineralizing fluids. The Pb isotopic compositions from Dolores and Yanque sulfides are similar, and are typical of the Tertiary magmatically-derived ores in this part of Peru. The hydrothermally altered rocks at Yanque have the same Pb isotopic compositions as the sulfides, thus confirming the hypothesis that the Yanque primary Zn-Pb mineralization may have been produced by hydrothermal circulation related to the emplacement of the Dolores Cu-porphyry, as it is the case of other porphyry Cu systems associated with polymetallic mineralization elsewhere. However, no simple genetic model for the mineralization involving just one fluid circulation episode is able to explain the data. The Yanque economic nonsulfide ore association consists of sauconite, hemimorphite, smithsonite and cerussite, which result from the weathering and alteration of the original sulfide mineralization. Zinc is allocated mainly in sauconite (Zn-smectite), rather than in carbonates: a factor strictly related to the prevailing siliciclastic character of the host rock. Distinctive features of the Yanque orebody are the comparable ore grades for both Pb and Zn (3.5% Zn and 3.7% Pb), and the inverse supergene chemical zoning. In fact, contrary to other supergene ores of this type, zinc prevails in the top zone of the Yanque deposit, whereas lead content increases with depth. Considering the different mobility of the two metals in solution, it may be assumed that most of the primary zinc that was the source for the Yanque nonsulfides was originally located far from the position occupied by the galena mineralization, whose remnants have been observed on site. Zinc sulfides may have been originally contained in the eroded rock volumes that surrounded the actual deposit: the zinc-rich solutions have possibly migrated through the siliciclastic Mara Fm. and precipitated the nonsulfide minerals by porosity filling and replacement processes. In this sense, the Yanque secondary Zn-Pb deposit could be considered as a special type of “Exotic” mineralization.

Spring temperature variability and eutrophication history inferred from sedimentary pigments in the varved sediments of Lake Żabińskie, north-eastern Poland, AD 1907–2008

Varved lake sediments are excellent natural archives providing quantitative insights into climatic and environmental changes at very high resolution and chronological accuracy. However, due to the multitude of responses within lake ecosystems it is often difficult to understand how climate variability interacts with other environmental pressures such as eutrophication, and to attribute observed changes to specific causes. This is particularly challenging during the past 100 years when multiple strong trends are superposed. Here we present a high-resolution multi-proxy record of sedimentary pigments and other biogeochemical data from the varved sediments of Lake Żabińskie (Masurian Lake District, north-eastern Poland, 54°N–22°E, 120 m a.s.l.) spanning AD 1907 to 2008. Lake Żabińskie exhibits biogeochemical varves with highly organic late summer and winter layers separated by white layers of endogenous calcite precipitated in early summer. The aim of our study is to investigate whether climate-driven changes and anthropogenic changes can be separated in a multi-proxy sediment data set, and to explore which sediment proxies are potentially suitable for long quantitative climate reconstructions. We also test if convoluted analytical techniques (e.g. HPLC) can be substituted by rapid scanning techniques (visible reflectance spectroscopy VIS-RS; 380–730 nm). We used principal component analysis and cluster analysis to show that the recent eutrophication of Lake Żabińskie can be discriminated from climate-driven changes for the period AD 1907–2008. The eutrophication signal (PC1 = 46.4%; TOC, TN, TS, Phe-b, high TC/CD ratios total carotenoids/chlorophyll-a derivatives) is mainly expressed as increasing aquatic primary production, increasing hypolimnetic anoxia and a change in the algal community from green algae to blue-green algae. The proxies diagnostic for eutrophication show a smooth positive trend between 1907 and ca 1980 followed by a very rapid increase from ca. 1980 ± 2 onwards. We demonstrate that PC2 (24.4%, Chl-a-related pigments) is not affected by the eutrophication signal, but instead is sensitive to spring (MAM) temperature (r = 0.63, pcorr < 0.05, RMSEP = 0.56 °C; 5-yr filtered). Limnological monitoring data (2011–2013) support this finding. We also demonstrate that scanning visible reflectance spectroscopy (VIS-RS) data can be calibrated to HPLC-measured chloropigment data and be used to infer concentrations of sedimentary Chl-a derivatives {pheophytin a + pyropheophytin a}. This offers the possibility for very high-resolution (multi)millennial-long paleoenvironmental reconstructions.

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.

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

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.