Constraints on the application of carbon isotope thermometry in high- to ultrahigh-temperature metamorphic terranes

Nine marble horizons from the granulite facies terrane of southern India were examined in detail for stable carbon and oxygen isotopes in calcite and carbon isotopes in graphite. The marbles in Trivandrum Block show coupled lowering of δ¹³C and δ¹⁸O values in calcite and heterogeneous single crystal δ¹³C values (? 1 to ? 10‰) for graphite indicating varying carbon isotope fractionation between calcite and graphite, despite the granulite facies regional metamorphic conditions. The stable isotope patterns suggest alteration of δ¹³C and δ¹⁸O values in marbles by infiltration of low δ¹³C–δ¹⁸O‐bearing fluids, the extent of alteration being a direct function of the fluid‐rock ratio. The carbon isotope zonation preserved in graphite suggests that the graphite crystals precipitated/recrystallized in the presence of an externally derived CO₂‐rich fluid, and that the infiltration had occurred under high temperature and low f_O2 conditions during metamorphism. The onset of graphite precipitation resulted in a depletion of the carbon isotope values of the remaining fluid+calcite carbon reservoir, following a Rayleigh‐type distillation process within fluid‐rich pockets/pathways in marbles resulting in the observed zonation. The results suggest that calcite–graphite thermometry cannot be applied in marbles that are affected by external carbonic fluid infiltration. However, marble horizons in the Madurai Block, where the effect of fluid infiltration is not detected, record clear imprints of ultrahigh temperature metamorphism (800–1000 °C), with fractionations reaching <2‰. Zonation studies on graphite show a nominal rimward lowering δ¹³C on the order of 1 to 2‰. The zonation carries the imprint of fluid deficient/absent UHT metamorphism. Commonly, calculated core temperatures are > 1000 °C and would be consistent with UHT metamorphism.     

Far field effects of Alpine plate tectonism in the Iberian microplate recorded by fault-related denudation in the Spanish Central System

Apatite fission track analysis was performed on 56 samples from central Spain to unravel the far field effects of the Alpine plate tectonic history of Iberia. The modelled thermal histories reveal complex cooling in the Cenozoic, indicative of intermittent denudation. Accelerated cooling events occurred across the Spanish Central System (SCS) from the Middle Eocene to Recent. These accelerated cooling events resulted in up to 2.8±0.9 km of denudation in the western Sierra de Gredos and 3.6±1.0 km in the central and eastern Gredos (assuming a paleogeothermal gradient of 28±5 °C and a surface temperature of 10 °C). The greatest amount of denudation (5.0±1.6 km) occurred in the Sierra de Guadarrama. Accompanying rock uplift was 4.7±1.0 and 5.9±1.6 km in the eastern Gredos and Guadarrama, respectively. Most denudation in the Gredos occurred from the Middle Eocene to the Early Miocene and can be related to the N–S stress field, induced by the Pyrenean compression. In the Guadarrama, the greatest denudation was Pliocene to Recent of age and seems related to the ongoing NW–SE Betic compression. The fact that the formation of the E–W trending Gredos coincides with the N–S Pyrenean compression and the creation of the present day morphology of the NE–SW trending Guadarrama with the younger NW–SE Betic compression, indicates that they record the far field effects of Alpine plate tectonics on Iberia. The trend of pre-existing lineaments was of major importance in influencing the style and magnitude of these of far field effects.     

Stability of earthquake moment tensor inversions: effect of the double-couple constraint

We show that spurious large non-double-couple components can be obtained in inversions for the full deviatoric moment tensor for shallow crustal earthquakes due to inaccurate Earth models. The traditional “best double-couple” solution does not in general provide an optimal estimate of a double-couple mechanism, and is only reliable when the non-double-couple component of the full deviatoric solution is small. The inverse problem for the moment tensors of the 1998 Antarctic Plate and 2000 Wharton Basin strike-slip earthquakes is shown in each case to have two well-fitting minima in the misfit function of pure double-couple solutions. Such pairs of solutions are most likely to exist for earthquakes which are close either to vertical strike-slip or to dip-slip on a fault plane dipping at 45°. It is shown theoretically that these pairs of solutions arise from the combination of the pure double-couple constraint and the instability of two elements of the moment tensor. No significant non-double-couple component is found for the shallow thrusting 1996 Biak, Indonesia earthquake.     

New Mexico structural zone—an analogue of the Colorado mineral belt

Updated aeromagnetic maps of New Mexico together with current knowledge of the basement geology in the northern part of the state (Sangre de Cristo and Sandia–Manzano Mountains)—where basement rocks were exposed in Precambrian-cored uplifts—indicate that the northeast-trending Proterozoic shear zones that controlled localization of ore deposits in the Colorado mineral belt extend laterally into New Mexico. The shear zones in New Mexico coincide spatially with known epigenetic precious- and base-metal ore deposits; thus, the mineralized belts in the two states share a common inherited basement tectonic setting. Reactivation of the basement structures in Late Cretaceous–Eocene and Mid-Tertiary times provided zones of weakness for emplacement of magmas and conduits for ore-forming solutions. Ore deposits in the Colorado mineral belt are of both Late Cretaceous–Eocene and Mid-Tertiary age; those in New Mexico are predominantly Mid-Tertiary in age, but include Late Cretaceous porphyry-copper deposits in southwestern New Mexico.The mineralized belt in New Mexico, named the New Mexico structural zone, is 250-km wide. The northwest boundary is the Jemez subzone (or the approximately equivalent Globe belt), and the southeastern boundary was approximately marked by the Santa Rita belt. Three groups (subzones) of mineral deposits characterize the structural zone: (1) Mid-Tertiary porphyry molybdenite and alkaline-precious-metal deposits, in the northeast segment of the Jemez zone; (2) Mid-Tertiary epithermal precious-metal deposits in the Tijeras (intermediate) zone; and (3) Late Cretaceous porphyry-copper deposits in the Santa Rita zone. The structural zone was inferred to extend from New Mexico into adjacent Arizona. The structural zone provides favorable sites for exploration, particularly those parts of the Jemez subzone covered by Neogene volcanic and sedimentary rocks.     

Cosmogenic-nuclide ages for New England coastal moraines, Martha's Vineyard and Cape Cod, Massachusetts, USA

Cosmogenic-nuclide exposure ages for 13 glacially transported boulders atop the Martha's Vineyard moraine, MA, USA, indicate that the southeastern margin of the Laurentide ice sheet reached its maximum extent during the last glaciation 23,200±500 yr ago. Another 10 age determinations from the younger Buzzards Bay moraine near Woods Hole, MA, indicate that this moraine complex was formed 18,800±400 yr ago. These ages correlate approximately with the terminations of cooling cycles recorded in Greenland ice cores and coeval ice-rafting events, suggesting that the marginal position of this sector of the ice sheet was tightly coupled to North Atlantic climate during the last glacial maximum.     

Probability Density Functions as Botanical-Climatological Transfer Functions for Climate Reconstruction

We present a new procedure, the pdf method (pdf=probability density function), for reconstructing Quaternary climate utilizing botanical data. The procedure includes the advantages of the indicator species method by considering the fossil and modern presence and absence of taxa rather than their frequencies, thus avoiding the need for modern analog plant communities. Overcoming the problematic use of absolute limits to describe climate response ranges is the main progress of the pdf method in comparison to the indicator species method. This advantage results from estimating probability density functions (pdfs) for monthly mean January and July temperature conditional on the present day occurrence of single taxa. Gaussian distributions sufficiently approximate pdfs of many, although not all, studied taxa. On the assumption of statistical independence, the procedure calculates a joint pdf as the product of the pdfs of the individual taxa. This algorithm weights each taxon according to the extent of its climate response range expressed by its covariance structure. We interpret the maximum of the resulting pdf as the most likely climate and its confidence interval as the uncertainty range. To avoid an artificial reduction of uncertainty arising from the use of numerous similar pdfs, a preselection method is proposed based on the Mahalanobis distance between pdfs. The pdf method was applied to the Carpinus phase of a profile from Gröbern, Germany, that spans the last interglaciation (Eemian). The reconstructed most probable January and July temperatures of about 0.0°C and 18.4°C barely differ from the modern values of −0.5°C and 18.3°C.     

Crack damage stress as a composite function of porosity and elastic matrix stiffness in dolomites and limestones

Twenty-five uniaxial compression tests were performed to determine stress at onset of dilation, referred to herein as “the crack damage stress,” in heterogeneous dolomites and limestones. A simplified model for crack damage stress (σ_cd) is developed here using porosity, elastic modulus, Poisson's ratio and three empirical coefficients. The model shows that when porosity decreases and elastic modulus increases, σ_cd rapidly increases and approaches its maximum value. On the other hand, when porosity increases and elastic modulus decreases, σ_cd rapidly decreases and approaches its minimum value. The proposed model is validated for six heterogeneous limestone and dolomite formations which are widely distributed in Israel.     

Fluid permeability of sedimentary rocks in a complete stress–strain process

The permeability of sedimentary rocks during triaxial compression tests was investigated to relate it to the complete strain–stress process. It was found that the permeability was not constant, but varied with the stress and strain states in the rocks. Prior to the peak strength, the permeability decreases with increasing load. A dramatic increase in permeability occurs during the strain softening period. In the present study, in situ measurements of fluid flow and pressure in floor strata was carried out in a double longwall mining face in the Yangzhuang colliery. These measurements show that both the strata pressure and the position with respect to the mining face influence the hydrogeologic properties. The permeability increased in the floor strata behind the mining face because those mining induced fractures opened as the strata pressure decreased. To better understand this change in hydraulic behavior around the mining faces, 3-D numerical modeling was carried out. The model provides the general picture of the stress distribution and failure zone both in the floor and roof strata. The field and model results demonstrate the importance of changes in the stress and strain states on the hydrogeology of a site.     

Pressure effect on Ti- or P-rich accessory mineral saturation in evolved granitic melts with differing K2O/Na2O ratios

The solubility of Ti- and P-rich accessory minerals has been examined as a function of pressure and K₂O/Na₂O ratio in two series of highly evolved silicate systems. These systems correspond to (a) alkaline, varying from alkaline to peralkaline with increasing K₂O/Na₂O ratio; and (b) strongly metaluminous (essentially trondhjemitic at the lowest K₂O/Na₂O ratio) and remaining metaluminous with increasing K₂O/Na₂O ratio (to 3). The experiments were conducted at a fixed temperature of 1000 °C, with water contents varying from 5 wt.% at low pressure (0.5 GPa), increasing through 5–10 wt.% at 1.5–2.5 GPa to 10 wt.% at 3.5 GPa. Pressure was extended outside the normal crustal range, so that the results may also be applied to derivation of hydrous silicic melts from subducted oceanic crust.

For the alkaline composition series, the TiO₂ content of the melt at Ti-rich mineral saturation decreases with increasing pressure but is unchanged with increasing K content (at fixed pressure). The P₂O₅ content of the alkaline melts at apatite saturation increases with increased pressure at 3.5 GPa only, but decreases with increasing K content (and peralkalinity). For the metaluminous composition series (termed as “trondhjemite-based series” (T series)), the TiO₂ content of the melt at Ti-rich mineral saturation decreases with increasing pressure and with increasing K content (at fixed pressure). The P₂O₅ content of the T series melts at apatite saturation is unchanged with increasing pressure, but decreases with increasing K content. The contrasting results for P and Ti saturation levels, as a function of pressure in both compositions, point to contrasting behaviour of Ti and P in the structure of evolved silicate melts. Ti content at Ti-rich mineral saturation is lower in the alkaline compared with the T series at 0.5 GPa, but is similar at higher pressures, whereas P content at apatite saturation is lower in the T series at all pressures studied. The results have application to A-type granite suites that are alkaline to peralkaline, and to I-type metaluminous suites that frequently exhibit differing K₂O/Na₂O ratios from one suite to another.     

The Port Mouton Shear Zone: intersection of a regional fault with a crystallizing granitoid pluton

The peraluminous tonalite–monzogranite Port Mouton Pluton is a petrological, geochemical, structural, and geochronological anomaly among the many Late Devonian granitoid intrusions of the Meguma Lithotectonic Zone of southern Nova Scotia. The most remarkable structural feature of this pluton is a 4-km-wide zone of strongly foliated (040/subvertical) monzogranites culminating in a narrow (10–30 m), straight, zone of compositionally banded rocks that extends for at least 3 km along strike. The banded monzogranites consist of alternating melanocratic and leucocratic compositions that are complementary to the overall composition of that part of the pluton, suggesting an origin by mineral–melt and mineral–mineral sorting. Biotite and feldspar are strongly foliated in the plane of the compositional bands. These compositional variations and foliations originated by a process of segregation flow during shearing of the main magma with a crystallinity of 55–75%. Subsequent minor brittle fracturing of feldspars, twinning of microcline, development of blocky sub-grains in quartz, and kinking of micas demonstrate overprinting by a high-temperature deformation straddling the monzogranite solidus. Small folds and late sigmoidal dykes indicate dextral movement on the shear zone. This Port Mouton Shear Zone (PMSZ) is approximately co-linear with the only outcrops of Late Devonian mafic intrusions in the area, two of which are syn-plutonic with well-developed mingling textures in the marginal tonalite of the Port Mouton Pluton. Also closely co-linear with the mafic intrusions are a granitoid dyke that extends well beyond the outer contact of the Port Mouton Pluton, a swarm of large aligned angular xenolithic slabs, a zone of thin wispy schlieren banding, a large Be-bearing pegmatite, and a breccia pipe with abundant garnetiferous metapelitic xenoliths. In various ways, the shear zone may control all of these features. The Port Mouton Shear Zone is parallel to many other NE-trending faults and shear zones in the northern Appalachians, probably related to the docking of the Meguma Zone along the Cobequid–Chedabucto Fault system.