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.     

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.     

Rates of Deglaciation during the Last Glaciation and Holocene in the Cordillera Vilcanota-Quelccaya Ice Cap Region, Southeastern Perú

Moraine chronology is combined with digital topography to model deglacial rates of paleoglacier volumes in both the Huancané Valley on the west side of the Quelccaya Ice Cap and the Upismayo Valley on the northwest side of the Cordillera Vilcanota. The fastest rates of deglaciation (39×10⁻⁵ to 114×10⁻⁵ km³ yr⁻¹ and 112×10⁻⁵ to 247×10⁻⁵ km³ yr⁻¹ for each valley, respectively) were calculated for the most recent paleoglaciers, corresponding to the last few centuries. These results are consistent with observations in the Venezuelan Andes showing high rates of deglaciation since the Little Ice Age. These rates also fall within the range of 20th century rates of deglaciation measured on the Quelccaya Ice Cap (29×10⁻⁵ to 220×10⁻⁵ km³ yr⁻¹, Brecher and Thompson, 1993; Thompson, 2000). These results imply that rates of deglaciation may fluctuate significantly over time and that high rates of deglaciation may not be exclusive to the late 20th century. Equilibrium line altitude (ELA) depressions for the ice volumes of the last glaciation modeled here were computed as 230 m for the Quelccaya Ice Cap and 170 m for the Cordillera Vilcanota. Maximum ELA depressions are lower than previously published: <500 m for the Cordillera Vilcanota and <400 m for the Quelccaya Ice Cap. These lower values could imply a topographic control over paleoglacier extent.     

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.