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.     

An experimental study of grain scale melt segregation mechanisms in two common crustal rock types

Creation of pathways for melt to migrate from its source is the necessary first step for transport of magma to the upper crust. To test the role of different dehydration‐melting reactions in the development of permeability during partial melting and deformation in the crust, we experimentally deformed two common crustal rock types. A muscovite‐biotite metapelite and a biotite gneiss were deformed at conditions below, at and above their fluid‐absent solidus. For the metapelite, temperatures ranged between 650 and 800 °C at P_c=700 MPa to investigate the muscovite‐dehydration melting reaction. For the biotite gneiss, temperatures ranged between 850 and 950 °C at P_c=1000 MPa to explore biotite dehydration‐melting under lower crustal conditions. Deformation for both sets of experiments was performed at the same strain rate (ε.) 1.37×10^?5 s^?1. In the presence of deformation, the positive ΔV and associated high dilational strain of the muscovite dehydration‐melting reaction produces an increase in melt pore pressure with partial melting of the metapelite. In contrast, the biotite dehydration‐melting reaction is not associated with a large dilational strain and during deformation and partial melting of the biotite gneiss melt pore pressure builds more gradually. Due to the different rates in pore pressure increase, melt‐enhanced deformation microstructures reflect the different dehydration melting reactions themselves. Permeability development in the two rocks differs because grain boundaries control melt distribution to a greater extent in the gneiss. Muscovite‐dehydration melting may develop melt pathways at low melt fractions due to a larger volume of melt, in comparison with biotite‐dehydration melting, generated at the solidus. This may be a viable physical mechanism in which rapid melt segregation from a metapelitic source rock can occur. Alternatively, the results from the gneiss experiments suggest continual draining of biotite‐derived magma from the lower crust with melt migration paths controlled by structural anisotropies in the protolith.     

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.     

Cation distribution in amphiboles synthesized along the Ga analogue of the tremolite–aluminotschermakite join

 Amphiboles were synthesized from bulk compositions prepared along the join Ca_1.8Mg_5.2Si₈O₂₂(OH)₂–Ca_1.8Mg₃Ga₄Si₆O₂₂(OH)₂ hydrothermally at 750–850 °C and 1.0–1.8 GPa, and along the join Ca₂Mg₅Si₈O₂₂F₂–Ca₂Mg₃Ga₄Si₆O₂₂F₂, anhydrously at 1000 °C and 0.7 GPa to document how closely the tschermak-type substitution is obeyed in these analogues of aluminous amphiboles. Electron-microprobe analyses and Rietveld X-ray diffraction structure refinements were performed to determine cation site occupancies. The extent of Ga substitution was found to be limited in both joins, but with the fluorine series having about twice the Ga content (0.6 atoms per formula unit, apfu) of the hydroxyl-series amphiboles (0.3 apfu). The tschermak-type substitution was followed very closely in the hydroxyl series with essentially equal partitioning of Ga between tetrahedral and octahedral sites. The fluorine-series amphiboles deviated significantly from the tschermak-type substitution and, instead, appeared to follow a substitution that is close to a Ca-pargasite substitution of the type: ^[6]Ga³⁺+2^[4]Ga³⁺+1/2^[A] Ca²⁺ = ^[6]Mg²⁺+2^[4]Si⁴⁺+1/2^[A]□. Infrared spectroscopy revealed an inverse correlation between the intensity of the OH-stretching bands and the Ga content for the hydroxyl- and fluorine-series amphiboles. The direct correlation between the Ga and F content and inverse relationship between the Ga and OH content may be a general phenomenon present in other minerals and suggests, for example, that high F contents in titanite are controlled by the Al content of the host rock and that there may be similar direct Al–F correlations in tschermakitic amphiboles. Evidence for the possibility that Al (Ga) might substitute onto only a subset of the tetrahedral sites in tschermakitic amphiboles was sought but not observed in this study. Received: 5 March 2001 / Accepted: 31 July 2001     

Incorporation of Cr<Superscript>3+</Superscript> in dickite: a spectroscopic study

 We have investigated a well-ordered sample of natural Cr-bearing dickite from Nowa Ruda (Lower Silesia, Poland) using electron paramagnetic resonance (EPR) at X- and Q-band frequencies (9.42 and 33.97 GHz, respectively) and optical diffuse reflectance spectroscopy. The observation of the spin-forbidden transitions at 15500 and 14690 cm⁻¹ allows us to unambiguously identify the major contribution of octahedrally coordinated Cr³⁺ ions in the optical spectrum. The X- and Q-band EPR spectra show two superposed Cr³⁺ signals. The corresponding fine-structure parameters were determined at room temperature and 145 K. These results suggest the substitution of Cr³⁺ for Al³⁺ in equal proportions in the two unequivalent octahedral sites of the dickite structure. In kaolin group minerals, the distortion around Cr³⁺ ions (λ≈ 0.2–0.4) in Al sites is significantly less rhombic than that observed around Fe³⁺ ions (λ≈ 0.6–0.8). Received: 29 June 2001 / Accepted: 22 October 2001     

Low-temperature thermodynamic properties of disordered zeolites of the natrolite group

 The heat capacity of paranatrolite and tetranatrolite with a disordered distribution of Al and Si atoms has been measured in the temperature range of 6–309 K using the adiabatic calorimetry technique. The composition of the samples is represented with the formula (Na_1.90K_0.22Ca_0.06)[Al_2.24Si_2.76O₁₀]·nH₂O, where n=3.10 for paranatrolite and n=2.31 for tetranatrolite. For both zeolites, thermodynamic functions (vibrational entropy, enthalpy, and free energy function) have been calculated. At T=298.15 K, the values of the heat capacity and entropy are 425.1 ± 0.8 and 419.1 ±0.8 J K⁻¹ mol⁻¹ for paranatrolite and 381.0 ± 0.7 and 383.2 ± 0.7 J K⁻¹ mol⁻¹ for tetranatrolite. Thermodynamic functions for tetranatrolite and paranatrolite with compositions corrected for the amount of extraframework cations and water molecules have also been calculated. The calculation for tetranatrolite with two water molecules and two extraframework cations per formula yields: C _p (298.15)=359.1 J K⁻¹ mol⁻¹, S(298.15) −S(0)=362.8 J K⁻¹ mol⁻¹. Comparing these values with the literature data for the (Al,Si)-ordered natrolite, we can conclude that the order in tetrahedral atoms does not affect the heat capacity. The analysis of derivatives dC/dT for natrolite, paranatrolite, and tetranatrolite has indicated that the water- cations subsystem within the highly hydrated zeolite may become unstable at temperatures above 200 K. Received: 30 July 2001 / Accepted: 15 November 2001     

A three-dimensional parametric study of the use of soil nails for stabilising tunnel faces

Applying an effective nailing system at a tunnel heading, not only improves the stability of the tunnel heading and limits deformation at the tunnel face, but it also reduces volume loss during excavation and hence reduces ground surface settlement. The effectiveness of a soil nail system is affected by many factors such as the diameter and stiffness of the nails. In this paper, a systematic parametric study was conducted to study the axial rigidity of a nail, E_nA_n, for improving the stability of tunnel headings and reducing ground movements in stiff clay. The parametric study involved a series of three-dimensional elasto-plastic coupled-consolidation finite element analyses. The stability of the tunnel face is improved with increasing E_nA_n. For a given nail density applied at the tunnel face, an optimum axial rigidity of the nail (E_nA_n)_opt can be identified. The efficiency of the nailing system diminishes when (E_nA_n)_opt is reached. The use of a soil nailing system reduces the magnitude of stress relief at the tunnel heading during excavation. Thus, this reduction of stress relief minimises the amount of soil yielding and excess pore water pressure generated in the soil around the tunnel heading.