Isotopic constraints on growth conditions of multiphase calcite–pyrite–barite concretions in Carboniferous mudstones

Carbonate concretions in the Lower Carboniferous Caton Shale Formation contain diagenetic pyrite, calcite and barite in the concretion matrix or in different generations of septarian fissures. Pyrite was formed by sulphate reduction throughout the sediment before concretionary growth, then continued to form mainly in the concretion centres. The septarian calcites show a continuous isotopic trend from δ¹³C=?28·7‰ PDB and δ¹⁸O=?1·6‰ PDB through to δ¹³C=?6·9‰ PDB and δ¹⁸O=?14·6‰ PDB. This trend arises from (1) a carbonate source initially from sulphate reduction, to which was added increasing contributions of methanogenic carbonate; and (2) burial/temperature effects or the addition of isotopically light oxygen from meteoric water. The concretionary matrix carbonates must have at least partially predated the earliest septarian cements, and thus used the same carbonate sources. Consequently, their isotopic composition (δ¹³C=?12·0 to ?10·1‰ PDB and δ¹⁸O=?5·7 to ?5·6‰ PDB) can only result from mixing a carbonate cement derived from sulphate reduction with cements containing increasing proportions of carbonate from methanogenesis and, directly or indirectly, also from skeletal carbonate. Concretionary growth was therefore pervasive, with cements being added progressively throughout the concretion body during growth. The concretions contain barite in the concretion matrix and in septarian fissures. Barite in the earlier matrix phase has an isotopic composition (δ³⁴S=+24·8‰ CDT and δ¹⁸O=+16·4‰ SMOW), indicating formation from near‐surface, sulphate‐depleted porewaters. Barites in the later septarian phase have unusual isotopic compositions (δ³⁴S=+6 to +11‰ CDT and δ¹⁸O=+8 to +11‰ SMOW), which require the late addition of isotopically light sulphate to the porewaters, either from anoxic sulphide oxidation (using ferric iron) or from sulphate dissolved in meteoric water. Carbon isotope and biomarker data indicate that oil trapped within septarian fissures was derived from the maturation of kerogen in the enclosing sediments.     

Analogue study of particle segregation in pyroclastic density currents, with implications for the emplacement mechanisms of large ignimbrites

Abstract Analogue flume experiments were conducted to investigate the transport and sedimentation behaviour of turbulent pyroclastic density currents. The experimental currents were scaled approximately to the natural environment in three ways: (1) they were fully turbulent; (2) they had a very wide range of particle sizes and associated Rouse numbers (the ratio of particle settling velocity to effective turbulent eddy velocity in the current); and (3) they contained particles of two different densities. Two sets of surge‐type experiments were conducted in a 5 m long, water‐filled lock‐exchange flume at five different volumetric particle concentrations from 0·6% to 23%. In one set (one‐component experiments), the currents contained just dense particles; in the other set (two‐component experiments), they contained both light and dense particles in equal volume proportions. In both sets of experiments, the population of each component had a log‐normal size distribution. In the two‐component experiments, the size range of the light particle population was selected in order to be in hydrodynamic equivalence with that of the dense particles. Dense particles were normally graded, both vertically and downstream, in the deposits from both sets of experiments. The mass loading (normalized to the initial mass of the suspension) and grain size of the dense component in the deposits decreased with distance from the reservoir and were insensitive to initial total particle concentration in the currents. On the other hand, in the two‐component experiments, the light particles were extremely sensitive to concentration. They were deposited in hydrodynamic equivalence with the dense particles from dilute currents, but were segregated efficiently at concentrations higher than a few per cent. With increasing particle concentration, the large, light particles were carried progressively further down the flume because of buoyancy effects. Deposits from the high‐concentration currents exhibited reverse vertical grading of the large, light particles. Efficient segregation of the light component was observed even if the bulk density of the current was less than that of the light particles. In both sets of experiments, marked inflexions in the rate of downstream decline in mass loading and maximum grain size of the dense component can be attributed to the presence of two different particle settling regimes in the flow: (1) particles with Rouse numbers >2·5, which did not respond to the turbulence and settled rapidly; and (2) particles with Rouse numbers <2·5, which followed the turbulent eddies and settled slowly. The results are applied to the transport and sedimentation dynamics of pyroclastic density currents that generate large, widespread ignimbrites. Field data fail to reveal significant departures from aerodynamic equivalence between pumice and lithic clasts in three such ignimbrites: the particulate loads of some large ignimbrites are transported principally in turbulent suspensions of low concentration. In some ignimbrites, the well‐developed inflexions in curves of maximum lithic (ML) size vs. distance can be attributed to the existence of distinct high and low Rouse number particle settling regimes that mark the transition from an overcharged state to one in which the residual particulate load is transported more effectively by turbulence.     

High-T, low-P metamorphism in the Palaeoproterozoic Halls Creek Orogen, northern Australia: the middle crustal response to a mantle-related transient thermal pulse

High‐T, low‐P metamorphic rocks of the Palaeoproterozoic central Halls Creek Orogen in northern Australia are characterised by low radiogenic heat production, high upper crustal thermal gradients (locally exceeding 40 °C km^?1) sustained for over 30 Myr, and a large number of layered mafic‐ultramafic intrusions with mantle‐related geochemical signatures. In order to account for this combination of geological and thermal characteristics, we model the middle crustal response to a transient mantle‐related heat pulse resulting from a temporary reduction in the thickness of the mantle lithosphere. This mechanism has the potential to raise mid‐crustal temperatures by 150–400 °C within 10–20 Myr following initiation of the mantle temperature anomaly, via conductive dissipation through the crust. The magnitude and timing of maximum temperatures attained depend strongly on the proximity, duration and lateral extent of the thermal anomaly in the mantle lithosphere, and decrease sharply in response to anomalies that are seated deeper than 50–60 km, maintained for <5 Myr in duration and/or have half‐widths <100 km. Maximum temperatures are also intimately linked to the thermal properties of the model crust, primarily due to their influence on the steady‐state (background) thermal gradient. The amplitudes of temperature increases in the crust are principally a function of depth, and are broadly independent of crustal thermal parameters. Mid‐crustal felsic and mafic plutonism is a predictable consequence of perturbed thermal regimes in the mantle and the lowermost crust, and the advection of voluminous magmas has the potential to raise temperatures in the middle crust very quickly. Although pluton‐related thermal signatures significantly dissipate within <10 Myr (even for very large, high‐temperature intrusive bodies), the interaction of pluton‐ and mantle‐related thermal effects has the potential to maintain host rock temperatures in excess of 400–450 °C for up to 30 Myr in some parts of the mid‐crust. The numerical models presented here support the notion that transient mantle‐related heat sources have the capacity to contribute significantly to the thermal budget of metamorphism in high‐T, low‐P metamorphic belts, especially in those characterised by low surface heat flow, very high peak metamorphic geothermal gradients and abundant mafic intrusions.     

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.     

Basement-controlled faulting of Paleozoic strata in southern Ontario, Canada: new evidence from geophysical lineament mapping

Basement fault reactivation is now recognized as an important control on sedimentation and fault propagation in intracratonic basins. In southern Ontario, the basement consists of complexly structured mid-Proterozoic (ca. 1.2 Ga) crystalline rocks and metasedimentary rocks that are overlain by up to 1500 m of Paleozoic sedimentary strata. Reactivation of basement structures is suspected to control the location of Paleozoic fault and fracture systems, but evaluation has been hindered by a limited understanding of the regional structural characteristics of the buried basement. New aeromagnetic- and gravimetric-lineament mapping presented in this paper better resolves the location of basement discontinuities and provides further evidence for basement controls on the distribution of Paleozoic fault and fracture systems. Lineament mapping was facilitated by reprocessing and digital image enhancement (micro-levelling, regional residual separation, derivative filtering) of existing regional gravity and aeromagnetic datasets. Reprocessed images identify new details of the structural fabric of the basement below southern Ontario and delineate several previously unrecognized aeromagnetic and gravity lineaments and linear zones. Linear zones parallel the projected trends of mid-Proterozoic terrane boundaries identified by field mapping on the exposed shield to the north of the study area, and are interpreted as zones of shearing and basement faulting. Mapped aeromagnetic and gravity lineaments show similar trends to Paleozoic faults and fracture networks and broad zones of seismicity in southern Ontario. These new data support an ‘inheritance model’ for Paleozoic faulting, involving repeated reactivation and upward propagation of basement faults and fractures into overlying cover strata.     

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.     

Crustal structure of the Ruby Mountains and southern Carlin Trend region, Nevada, from magnetotelluric data

We have collected about 150 magnetotelluric (MT) soundings in northeastern Nevada in the region of the Ruby Mountains metamorphic core complex uplift and southern Carlin mineral trend, in an effort to illuminate controls on core complex evolution and deposition of world-class gold deposits. The region has experienced a broad range of tectonic events including several periods of compressional and extensional deformation, which have contributed to the total expression of electrical resistivity. Most of the soundings reside in three east–west profiles across increasing degrees of core uplift to the north (Bald Mountain, Harrison Pass, and Secret Pass latitudes). One short cross-line was also taken to assess an east–west structure to the north of the northern profile. Model resistivity cross-sections were derived from the MT data using a 2-D inversion algorithm, which damps departures of model parameters from an a priori structure. Geological interpretation of the resistivity combines previous seismic, potential field and isotope models, structural and petrological models for regional compression and extension, and detailed structural/stratigraphic interpretations incorporating drilling for petroleum and mineral exploration. To first order, the resistivity structure is one of a moderately conductive, Phanerozoic sedimentary section fundamentally disrupted by intrusion and uplift of resistive crystalline rocks. Late Devonian and early Mississippian shales of the Pilot and Chainman Formations together form an important conductive marker sequence in the stratigraphy and show pronounced increases in conductance (conductivity–thickness product) from east to west. These increases are attributed to graphitization caused by Elko–Sevier era compressional shear deformation and possibly by intrusive heating. The resistive crystalline central massifs adjoin the host stratigraphy across crustal-scale, steeply dipping fault zones. The zones provide pathways to the lower crust for heterogeneous, upper crustal induced, electric current flow. Resistive core complex crust appears steeply bounded under the middle of the neighboring grabens and not to deepen at a shallow angle to arbitrary distances to the west. The numerous crustal breaks imaged with MT may contribute to the low effective elastic thickness (T_e) estimated regionally for the Great Basin and exemplify the mid-crustal, steeply dipping slip zones in which major earthquakes nucleate. An east–west oriented conductor in the crystalline upper crust spans the East Humboldt Range and northern Ruby Mountains. The conductor may be related to nearby graphitic metasediments, with possible alteration by middle Tertiary magmatism. Lower crustal resistivity everywhere under the profiles is low and appears quasi one-dimensional. It is consistent with a low rock porosity (<1 vol.%) containing hypersaline brines and possible water-undersaturated crustal melts, residual to the mostly Miocene regional extension. The resistivity expression of the southern Carlin Trend (CT) in the Pinon Range is not a simple lineament but rather a family of structures attributed to Eocene intrusion, stratal deformation, and alteration/graphitization. Substantial reactivation or overprinting by core complex uplift or Basin–Range extensional events seems likely. We concur with others that the Carlin Trend may result in part from overlap of the large Eocene Northeast Nevada Volcanic Field with Precambrian–Paleozoic deep-water clastic source rocks thickening abruptly to the west of the Pinon Range, and projecting to the north–northwest.     

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.