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.     

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.     

Hydraulic Characteristics of Rough Fractures in Linear Flow under Normal and Shear Load

Summary A hydro-mechanical testing system, which is capable of measuring both the flow rates and the normal and shear displacement of a rock fracture, was built to investigate the hydraulic behaviour of rough tension fractures. Laboratory hydraulic tests in linear flow were conducted on rough rock fractures, artificially created using a splitter under various normal and shear loading. Prior to the tests, aperture distributions were determined by measuring the topography of upper and lower fracture surfaces using a laser profilometer. Experimental variograms of the initial aperture distributions were classified into four groups of geostatistical model, though the overall experimental variograms could be well fitted to the exponential model. The permeability of the rough rock fractures decayed exponentially with respect to the normal stress increase up to 5 MPa. Hydraulic behaviours during monotonic shear loading were significantly affected by the dilation occurring until the shear stress reached the peak strength. With the further dilation, the permeability of the rough fracture specimens increased more. However, beyond shear displacement of about 7 to 8 mm, permeability gradually reached a maximum threshold value. The combined effects of both asperity degradation and gouge production, which prohibited the subsequent enlargement of mean fracture aperture, mainly caused this phenomenon. Permeability changes during cyclic shear loading showed somewhat irregular variations, especially after the first shear loading cycle, due to the complex interaction from asperity degradations and production of gouge materials. The relation between hydraulic and mechanical apertures was analyzed to investigate the valid range of mechanical apertures to be applied to the cubic law. Received June 12, 2001; accepted February 26, 2002 Published online September 2, 2002     

Trace element mobility during sub-seafloor alteration of basaltic glass from Ocean Drilling Program site 953 (off Gran Canaria)

Trace element concentrations of altered basaltic glass shards (layer silicates) and zeolites in volcaniclastic sediments drilled in the volcanic apron northeast of Gran Canaria during Ocean Drilling Program (ODP) leg 157 document variable element mobilities during low-temperature alteration processes in a marine environment. Clay minerals (saponite, montmorillonite, smectite) replacing volcanic glass particles are enriched in transition metals and rare earth elements (REE). The degree of retention of REE within the alteration products of the basaltic glass is correlated with the field strength of the cations. The high field-strength elements are preferentially retained or enriched in the alteration products by sorption through clay minerals. Most trace elements are enriched in a boundary layer close to the interface mineral-altered glass. This boundary layer has a key function for the physico-chemical conditions of the subsequent alteration process by providing a large reactive surface and by lowering the fluid permeability. The release of most elements is buffered by incorporation into secondary precipitates (sodium-rich zeolites, phillipsite, Fe- and Mn-oxides) as shown by calculated distribution coefficients between altered glasses and authigenic minerals. Chemical fluxes change from an open to a closed system behavior during prograde low-temperature alteration of volcaniclastic sediments with no significant trace metal flux from the sediment to the water column.     

Occurrence of naphthalenesulfonates and their condensates with formaldehyde in a landfill leachate and their transport behavior in groundwater of the Upper Rhine Valley, Germany

. Naphthalenesulfonates and their condensation products with formaldehyde are manufactured in the chemical industry for broad application as tanning agents, dispersing agents, and as superplasticizers for concrete. With the disposal of waste they have found their way into the aquatic environment. Five monomeric naphthalenesulfonates and a group of dimeric sulfonated naphthalene-formaldehyde condensates (SNFC) were identified in the leachate and in groundwater observation wells downstream in the plume of the Karlsruhe-West landfill, situated in the Upper Rhine Valley in southwest Germany. For the main monomer of technical SNFC products, 2-naphthalenesulfonate (2-NS), the concentration in the leachate was 170 µg l^–1, which decreased to 9.5 µg l^–1 at a distance of about 300 m from the boundary of the landfill. In the Merdingen tracer test field, the transport of these compounds was studied under more controlled conditions in two tracer experiments with (1) two model compounds, and (2) a technical SNFC product and uranine as a reference tracer. In these experiments, the monomeric naphthalenesulfonates and the SNFC with n=2 behave as conservative tracers. Thus, the findings of the landfill study were supported by these results. Higher molecular SNFC were strongly retarded, which is attributed to adsorption to soil particles. The results of the second tracer experiment suggest a degradation of 2-NS and 2,6-naphthalenedisulfonate (2,6-NDS) after adaptation of the microorganisms in the groundwater aquifer as a consequence of the first tracer experiment.     

Water quality variation and clam growth: Is pH really a non-issue in estuaries?

A tandem deployment system was used to critically evaluate relationships between important water chemistry parameters (pH, salinity, dissolved oxygen) and biotic performance based on clam growth. The effects of environmental conditions on growth of juvenile clams,Mercenaria mercenaria, were determined after 7-day field deployments in cages at reference sites from 1998 to 2000. Continuous measurements of the overlying water chemistry parameters were monitored by deploying an in situ water quality instrument (Hydrolab Datasonde) at the same time. While salinity was identified as an important determinant of clam growth over wide salinity ranges (10–35‰), pH was also found to be a very important parameter, especially in low-salinity regimes (<25‰). Average pH measurements ranged from 7.2 to 7.8; minimal pHs ranged from 6.9 to 7.6. The results indicated that when average pH levels fell below 7.5 or minimum pH levels fell below 7.2, growth rates were <50% that of clams deployed under higher pH conditions. Estuarine systems are generally perceived as being well-buffered so pH is frequently assumed to be unimportant, but our results suggest that pH levels can decline in estuarine systems to levels that can adversely affect biological responses. The potential impacts on biological resources of even moderate decreases in pH, particularly in systems that naturally tend to have lower pH conditions, may be more important than previously realized.     

Impact of suboxia on sinking particulate organic carbon: Enhanced carbon flux and preferential degradation of amino acids via denitrification

Fluxes of particulate organic carbon (POC) through the oxygen deficient waters in the eastern tropical North Pacific were found to be relatively less attenuated with depth than elsewhere in the eastern North Pacific. The attenuation coefficient (b) for the flux was found to be 0.40 versus the composite value of 0.86 determined by Martin et al. (1987). To examine this further, sinking POC was collected using sediment traps and allowed to degrade in oxic and suboxic experiments. Using a kinetic model, it was found that degradation proceeded at similar rates (roughly 0.8 day⁻¹) under oxic and suboxic conditions, but a greater fraction of bulk POC was resistant to degradation in the suboxic experiments (61% vs. 23%). Amino acids accounted for 37% of POC collected at 75m, but following degradation the value dropped to 17% and 16% in the oxic and suboxic experiments respectively. POC collected from 500m was 10% amino acids. The non-AA component of POC collected at 75m was not degraded under suboxic conditions, while under oxic conditions it was. These results suggest that microbes degrading OC under suboxic conditions via denitrification preferentially utilize nitrogen-rich amino acids. This preferential degradation of amino acids suggests that 9% more nitrogen may be lost via water column denitrification than is accounted for when a more “Redfieldian” stoichiometry for POC is assumed.