Interpreting tracer breakthrough tailing in a conduit-dominated karstic aquifer

Breakthrough tailing has been observed during dye-tracing recovery tests in the Norville aquifer system (chalk), France. Karst-conduit flow and transport parameters were assessed using two different interpretative methods: the linear graphical method and the Chatwin method (implemented in the Qtracer2 program). The linear graphical method was used to model the observed tailing effects, which was explained by a second smaller delayed breakthrough curve. By comparing the results of tracer-test interpretation for the two methods, it was possible to relate the area of this second curve to the importance of turbulent flow in spring discharge. The more turbulent the flow, the less important the contribution of the second breakthrough curve and the tailing effect. The observed tailing could possibly be controlled by hydrodynamics to a greater extent than usually expected, the tailing effects being mostly attributed to diffusion phenomena. Tailing effects were expected to increase with discharge and the piezometric level, which would have resulted in overpressure in conduits, fissure flooding, etc. Instead, breakthrough tailing tended to disappear with increasing aquifer discharge, which would support the hypothesis of there being mostly hydrodynamic-controlled tailing effects instead of matrix- or fissure-diffusion.An erratum to this article can be found at     

Synergy of multiple geophysical approaches to unravel explosive eruption conduit and source dynamics - A case study from Stromboli

An effective approach to understanding the dynamics of explosive volcanic eruptions and the conduit systems that drive them is through synergy of multiple data sets. Three data sets that lend themselves to ease of integration are seismic, infrasonic and thermal. Although approaches involving these data have been used to record volcanological phenomena since 1862, 1955 and 1965, respectively, their integrated use has only developed since 1999. When combined, these three data sets allow constraint of shallow system geometry and the dynamics of the explosive events that occur within that system. Using Stromboli volcano (Italy) as a case study, we review the complete range of geochemical and geophysical studies that can be applied. In doing so, we aim to show how integration of these diverse studies allows insights into a plumbing system and the dynamics of the eruptive activity that the system feeds. When combined at Stromboli, these data provide constraint of multiple system parameters including chamber depths, gas and magma fluxes, shallow system magma residence times, explosion source depths, and the rise/ejection velocities of ascending gas slugs and ejecta. In turn, these results allow various conduit and eruption dynamic models to be applied and tested.The persistent and repeated mildly explosive events that characterize Stromboli have been modeled in terms of the coalescence of gas within the magma to form large gas slugs that ascend the remaining portion of the conduit to burst at the free surface. Our integrated seismic, infrasonic and thermal data sets indicate that gas coalescence occurs at a depth of ∼260 m, with a typical event frequency of ∼9/h. Infrasonic and thermal data show the explosion source to be located 20-220 m below the vent. Thermal data give emission velocities for the ejected fragments of 8-20 m/s, which converts to gas jet velocities of 23-39 m/s. Tracking these parameters in space and time shows that, although eruptions at Stromboli can be grouped into two characteristic types (simple and complex-each of which characterizes a particular crater, NE and SW, respectively), events within each type show significant short-term variability. The system does, however, appear robust, maintaining its characteristic strombolian eruption style after significant effusive phases and more energetic explosive events.     

Conduit-vent structures and related proximal deposits in the Las Cañadas caldera, Tenerife, Canary Islands

The Las Cañadas caldera wall and the outer slopes of the caldera provide three-dimensional exposures of numerous proximal-welded fallout deposits and have been mapped in detail. As a result, some parts of the Ucanca and Guajara Formations of the stratigraphy of Martí et al. (1994) have been divided into members that correspond to individual eruptions. Mapping has also revealed the occurrence of conduit-vent structures associated with proximal-welded fallout deposits. Conduit-vent structures consist of an upper flaring area and a lower narrow conduit. Conduit-vent geometry and dimensions include cylindrical plugs and eruptive fissures steeply dipping towards the caldera depression and elongated vents. The flaring area can be rather asymmetric and is usually filled by down-vent rheomorphic flow of the proximal fallout deposit. The lower conduits are filled by lava plug, agglutination of juveniles onto conduit walls and dyke intrusion with eventual dome extrusion. The eruption dynamics of welded fallout deposits and magma fragmentation within the conduit are consistent with an evolution from explosive to effusive. In this context conduit flow regimes evolve from turbulent to annular flow in which the conduit is progressively choked, and laminar flow leading to the final conduit closure.     

Crystallization of microlites during magma ascent: the fluid mechanics of 1980–1986 eruptions at Mount St Helens

Eruptions of Mount St Helens (Washington, USA) decreased in intensity and explosivity after the main May 18, 1980 eruption. As the post-May 18 eruptions progressed, albitic plagioclase microlites began to appear in the matrix glass, although the bulk composition of erupted products, the phenocryst compositions and magmatic temperatures remained fairly constant. Equilibrium experiments on a Mount St Helens white pumice show that at 160 MPa water pressure and 900°C, conditions deduced for the 8 km deep magma storage zone, the stable plagioclase is An₄₇. The microlites in the natural samples, which are more albitic, had to grow at lower water pressures during ascent. Isothermal decompression experiments reported here demonstrate that a decrease in water pressure from 160 to 2 MPa over four to eight days is capable of producing the albitic groundmass plagioclase and evolved melt compositions observed in post-May 18 1980 dacites. Because groundmass crystallization occurs over a period of days during and after decreases in pressure, microlite crystallization in the Mount St Helens dacites must have occurred during the ascent of each magma batch from a deep reservoir rather than continuously in a shallow holding chamber. This is consistent with data on the kinetics of amphibole breakdown, which require that a significant portion of magma vented in each eruption ascended from a depth of at least 6.5 km (160 MPa water pressure) in a matter of days. The size and shape of the microlite population have not been studied because of the small size of the experimental samples; it is possible that the texture continues to mature long after chemical equilibrium is approached. As the temperature, composition, crystal content and water content of magma in the deep reservoir remained approximately constant from May 1980 to at least March 1982, the spectacular decrease in eruption intensity during this period cannot be attributed to changes in viscosity or density of the magma. Simple fluld mechanical considerations indicate, however, that the observed changes in mass flux of magma can be modelled by a five-fold decrease in conduit radius from 35 to 7 m, produced perhaps by plating of magma along the conduit walls. The decreased ascent rates which accompanied the decrease in conduit radius can explain the change from closed-system to open-system degassing and the shift from explosive to effusive eruptions during 1980.     

Development of fluid conduits in the auriferous shear zones of the Hutti Gold Mine, India: evidence for spatially and temporally heterogeneous fluid flow

The gold mineralization of the Hutti Mine is hosted by nine parallel, N–S trending, steeply dipping, 2–10 m wide shear zones, that transect Archaean amphibolites. The shear zones were formed after peak metamorphism during retrograde ductile D₂ shearing in the lower amphibolite facies. They were reactivated in the lower to mid greenschist facies by brittle–ductile D₃ shearing and intense quartz veining. The development of a S₂–S₃ crenulation cleavage facilitates the discrimination between the two deformation events and contemporaneous alteration and gold mineralization. Ductile D₂ shearing is associated with a pervasively developed distal chlorite–sericite alteration assemblage in the outer parts of the shear zones and the proximal biotite–plagioclase alteration in the center of the shear zones. D₃ is characterized by development of the inner chlorite-K-feldspar alteration, which forms a centimeter-scale alteration halo surrounding the laminated quartz veins and replaces earlier biotite along S₃. The average size of the laminated vein systems is 30–50 m along strike as well as down-dip and 2–6 m in width.Mass balance calculations suggest strong metasomatic changes for the proximal biotite–plagioclase alteration yielding mass and volume increase of ca. 16% and 12%, respectively. The calculated mass and volume changes of the distal chlorite–sericite alteration (ca. 11%, ca. 8%) are lower. The decrease in δ¹⁸O values of the whole rock from around 7.5‰ for the host rocks to 6–7‰ for the distal chlorite–sericite and the proximal biotite–plagioclase alteration and around 5‰ for the inner chlorite-K-feldspar alteration suggests hydrothermal alteration during two-stage deformation and fluid flow.The ductile D₂ deformation in the lower amphibolite facies has provided grain scale porosities by microfracturing. The pervasive, steady-state fluid flow resulted in a disseminated style of gold–sulfide mineralization and a penetrative alteration of the host rocks. Alternating ductile and brittle D₃ deformation during lower to mid greenschist facies conditions followed the fault-valve process. Ductile creep in the shear zones resulted in a low permeability environment leading to fluid pressure build-up. Strongly episodic fluid advection and mass transfer was controlled by repeated seismic fracturing during the formation of laminated quartz(-gold) veins. The limitation of quartz veins to the extent of earlier shear zones indicate the importance of pre-existing anisotropies for fault-valve action and economic gold mineralization.     

Simulation of the development of gypsum maze caves

The development of gypsum maze caves under artesian conditions has been simulated. The numerical model simulations show that the evolution of maze caves in this type of setting requires structural preferences such as laterally extended fissure networks in a horizon of the gypsum layer. Without any structural preferences vertical shafts rather than maze caves are predicted to develop. The most important stage for the development of horizontal caves under artesian conditions is found to be the initial karstification period. During this period the structure of the mature conduit system is established. The solutional enlargement of conduits is spatially extended, total dissolution rates are higher than the later ones.     

Modelling eruption cycles and decay of mud volcanoes

Recent debates about the eruptive behavior of mud volcanoes and their activation mechanisms have been driven particularly by the LUSI eruption in Indonesia that resulted in huge commercial and cultural damages. Numerical modeling of mud volcanoes, of which few exist, can provide insight into eruptive behavior and contribute to hazard evaluation. In this paper, we present a simple model to describe fluid escape from an underground reservoir through a conduit, extruded as a mud volcano at the surface. The governing equations result in oscillatory behavior, and we study the influence of changes in rheological properties of surrounding rock and fluid characteristics of the mud on extrusion dynamics. We focus on understanding long-term eruption behavior, flow cycles, and decay factors. Model results can be used to estimate the discharge rates and extruded volume from assumptions on the mud reservoir and conduit, or conversely, the reservoir or conduit properties from discharge rates.     

Observations of volcanic tremor during the January–February 2005 eruption of Mt. Veniaminof,Alaska

Mt. Veniaminof, Alaska Peninsula, is a stratovolcano with a summit ice-filled caldera containing a small intracaldera cone and active vent. From January 2 to February 21, 2005, Mt. Veniaminof erupted. The eruption was characterized by numerous small ash emissions (VEI 0 to 1) and accompanied by low-frequency earthquake activity and volcanic tremor. We have performed spectral analyses of the seismic signals in order to characterize them and to constrain their source. Continuous tremor has durations of minutes to hours with dominant energy in the band 0.5–4.0 Hz, and spectra characterized by narrow peaks either irregularly (non-harmonic tremor) or regularly spaced (harmonic tremor). The spectra of non-harmonic tremor resemble those of low-frequency events recorded simultaneously with surface ash explosions, suggesting that the source mechanisms might be similar or related. We propose that non-harmonic tremor at Mt. Veniaminof results from the coalescence of gas bubbles while low-frequency events are related to the disruption of large gas pockets within the conduit. Harmonic tremor, characterized by regular and quasi-sinusoidal waveforms, has duration of hours. Spectra containing up to five harmonics suggest the presence of a resonating source volume that vibrates in a longitudinal acoustic mode. An interesting feature of harmonic tremor is that frequency is observed to change over time; spectral lines move towards higher or lower values while the harmonic nature of the spectra is maintained. Factors controlling the variable characteristics of harmonic tremor include changes in acoustic velocity at the source and variations of the effective size of the resonator.     

On the relation between excess hydraulic head and flow in a conduit imbedded within a porous matrix

This article provides a partial answer to the question “What is the relation between excess hydraulic head and volume flux of water in a conduit within a porous matrix?”, focusing on the case that the forcing is steady. The conduit is modelled as a horizontal circular cylinder, imbedded within a porous matrix of rectangular cross section, having constant head prescribed on the sidewalls and being confined top and bottom. Laminar flow in the matrix is assumed to obey Darcy's law, while turbulent flow in the conduit is quantified using the Darcy–Weisbach equation. Analysis of the latter equation shows that the length scale of variations in the direction of the conduit is large compared with the scale of lateral and vertical variations. This permits separation of the full three-dimensional non-linear problem into a two-dimensional linear problem for head within the matrix and a one-dimensional non-linear problem for head within the conduit. Analytic solutions are obtained for the distribution of head in the matrix and in a conduit of either infinite or finite length. In both cases, the volume flux of water is proportional to the excess head to the 2/3 power, the conduit radius to the 5/3 power, the matrix permeability to the 1/3 power and gravity to the 1/3 power. The scale of variation of head along the conduit is proportional to the excess head to the ?1/3 power, the conduit radius to the 5/3 power, the matrix permeability to the ?2/3 power and gravity to the 1/3 power.     

Phreatic drainage conduits within quartz sandstone: Evidence from the Jurassic Precipice Sandstone, Carnarvon Range, Queensland, Australia

Discrete underground drainage conduits in quartz sandstones are far less common than in limestones. This paper provides field evidence from the quartzose Precipice Sandstone in the Carnarvon Range of south-central Queensland, Australia, for tubular underground drainage networks similar in many ways to limestone conduits. Diameters range from less than 1 or 2 cm to over 1.5 m, most display a near-circular to oval cross-section that seems to suggest phreatic or epiphreatic development, and the internal surfaces of many are case-hardened by secondary silica deposits. A number of the region's perennial springs appear to be fed by such tubes.The dominant vertical jointing of the quartz sandstone and relatively high permeability of the sandstone are important controls on tube formation. Solutional weathering of the sandstone is widespread, and is followed by the removal of loosened sand grains by flowing underground water, the process of ‘arenisation’. Tube development would appear to have been happening for a very long time, and may still be occurring. A model for tube network formation is proposed.These findings highlight our potentially poor understanding of groundwater flow within some quartz sandstones, and may have important groundwater management implications.