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.

The 1452 or 1453 A.D. Kuwae Eruption Signal Derived from Multiple Ice Core Records: Greatest Volcanic Sulfate Event of the Past 700 Years

We combined 33 ice core records, 13 from the Northern Hemisphere and 20 from the Southern Hemisphere, to determine the timing and magnitude of the great Kuwae eruption in the mid-15th century. We extracted volcanic deposition signals by applying a high-pass loess filter to the time series and examining peaks that exceed twice the 31 year running median absolute deviation. By accounting for the dating uncertainties associated with each record, these ice core records together reveal a large volcanogenic acid deposition event during 1453 - 1457 A. D. The results suggest only one major stratospheric injection from the Kuwae eruption and confirm previous findings that the Kuwae eruption took place in late 1452 or early 1453, which may serve as a reference to evaluate and improve the dating of ice core records. The average total sulfate deposition from the Kuwae eruption was 93 kg SO4/km(2) in Antarctica and 25 kg SO4/km(2) in Greenland. The deposition in Greenland was probably underestimated since it was the average value of only two northern Greenland sites with very low accumulation rates. After taking the spatial variation into consideration, the average Kuwae deposition in Greenland was estimated to be 45 kg SO4/km(2). By applying the same technique to the other major eruptions of the past 700 years our result suggests that the Kuwae eruption was the largest stratospheric sulfate event of that period, probably surpassing the total sulfate deposition of the Tambora eruption of 1815, which produced 59 kg SO4/km(2) in Antarctica and 50 kg SO4/km(2) in Greenland.

Impact of solar versus volcanic activity variations on tropospheric temperatures and precipitation during the Dalton Minimum

The aim of this work is to elucidate the impact of changes in solar irradiance and energetic particles versus volcanic eruptions on tropospheric global climate during the Dalton Minimum (DM, AD 1780–1840). Separate variations in the (i) solar irradiance in the UV-C with wavelengths λ < 250 nm, (ii) irradiance at wavelengths λ > 250 nm, (iii) in energetic particle spectrum, and (iv) volcanic aerosol forcing were analyzed separately, and (v) in combination, by means of small ensemble calculations using a coupled atmosphere–ocean chemistry–climate model. Global and hemispheric mean surface temperatures show a significant dependence on solar irradiance at λ > 250 nm. Also, powerful volcanic eruptions in 1809, 1815, 1831 and 1835 significantly decreased global mean temperature by up to 0.5 K for 2–3 years after the eruption. However, while the volcanic effect is clearly discernible in the Southern Hemispheric mean temperature, it is less significant in the Northern Hemisphere, partly because the two largest volcanic eruptions occurred in the SH tropics and during seasons when the aerosols were mainly transported southward, partly because of the higher northern internal variability. In the simulation including all forcings, temperatures are in reasonable agreement with the tree ring-based temperature anomalies of the Northern Hemisphere. Interestingly, the model suggests that solar irradiance changes at λ < 250 nm and in energetic particle spectra have only an insignificant impact on the climate during the Dalton Minimum. This downscales the importance of top–down processes (stemming from changes at λ < 250 nm) relative to bottom–up processes (from λ > 250 nm). Reduction of irradiance at λ > 250 nm leads to a significant (up to 2%) decrease in the ocean heat content (OHC) between 0 and 300 m in depth, whereas the changes in irradiance at λ < 250 nm or in energetic particles have virtually no effect. Also, volcanic aerosol yields a very strong response, reducing the OHC of the upper ocean by up to 1.5%. In the simulation with all forcings, the OHC of the uppermost levels recovers after 8–15 years after volcanic eruption, while the solar signal and the different volcanic eruptions dominate the OHC changes in the deeper ocean and prevent its recovery during the DM. Finally, the simulations suggest that the volcanic eruptions during the DM had a significant impact on the precipitation patterns caused by a widening of the Hadley cell and a shift in the intertropical convergence zone.

Volcanic Influence on European Summer Precipitation through Monsoons: Possible Cause for “Years without Summer”*

Strong tropical volcanic eruptions have significant effects on global and regional temperatures. Their effects on precipitation, however, are less well understood. Analyzing hydroclimatic anomalies after 14 strong eruptions during the last 400 years in climate reconstructions and model simulations, a reduction of the Asian and African summer monsoons and an increase of south-central European summer precipitation in the year following the eruption was found. The simulations provide evidence for a dynamical link between these phenomena. The weaker monsoon circulations weaken the northern branch of the Hadley circulation, alter the atmospheric circulation over the Atlantic–European sector, and increase precipitation over Europe. This mechanism is able to explain, for instance, the wet summer in parts of Europe during the “year without a summer” of 1816, which up to now has not been explained. This study underlines the importance of atmospheric teleconnections between the tropics and midlatitudes to better understand the regional climate response to stratospheric volcanic aerosols.

Ice Core Evidence for a Second Volcanic Eruption Around 1809 in the Northern Hemisphere

A volcanic signal observed in ice cores from both polar regions six years prior to Tambora is attributed to an unknown tropical eruption in 1809. Recovery of dacitic tephra from the 1809 horizon in a Yukon ice core ( Eclipse) that is chemically distinct from andesitic 1809 tephra found in Antarctic ice cores indicates a second eruption in the Northern Hemisphere at this time. Together with the similar magnitude and timing of the 1809 volcanic signal in the Arctic and Antarctic, this could suggest a large tropical eruption produced the sulfate and Antarctic tephra and a minor Northern Hemisphere eruption produced the Eclipse tephra. Nonetheless, the possibility that there were coincidental eruptions of similar magnitude in both hemispheres, rather than a single tropical eruption, should not be discounted. Correctly attributing the source of the 1809 volcanic signal has important implications for modeling the magnitude and latitudinal distribution of volcanic radiative forcing.

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.

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

A tephra lattice for Greenland and a reconstruction of volcanic events spanning 25–45 ka b2k

Tephra layers preserved within the Greenland ice-cores are crucial for the independent synchronisation of these high-resolution records to other palaeoclimatic archives. Here we present a new and detailed tephrochronological framework for the time period 25,000 e 45,000 a b2k that brings together results from 4 deep Greenland ice-cores. In total, 99 tephra deposits, the majority of which are preserved as cryptotephra, are described from the NGRIP, NEEM, GRIP and DYE-3 records. The major element signatures of single glass shards within these deposits indicate that 93 are basaltic in composition all originating from Iceland. Specifically, 43 originate from Grimsv € otn, 20 are thought to be sourced from the Katla volcanic system and 17 show affinity to the Kverkfj € oll system. Robust geochemical characterisations, independent ages derived from the GICC05 ice-core chronology, and the stratigraphic positions of these deposits relative to the Dansgaard-Oeschger climate events represent a key framework that provides new information on the frequency and nature of volcanic events in the North Atlantic region between GS-3 and GI-12. Of particular importance are 19 tephra deposits that lie on the rapid climatic transitions that punctuate the last glacial period. This framework of well-constrained, time-synchronous tie-lines represents an important step towards the independent synchronisation of marine, terrestrial and ice-core records from the North Atlantic region, in order to assess the phasing of rapid climatic changes during the last glacial period.

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.