Characterization of juvenile pyroclasts from the Kos Plateau Tuff (Aegean Arc): insights into the eruptive dynamics of a large rhyolitic eruption

Silicic pumices formed during explosive volcanic eruptions are faithful recorders of the state of the magma in the conduit, close to or at the fragmentation level. We have characterized four types of pumices from the non-welded rhyolitic Kos Plateau Tuff, which erupted 161,000 years ago in the East Aegean Arc, Greece. The dominant type of pumice (>90 vol.%) shows highly elongated tubular vesicles. These tube pumices occur throughout the eruption. Less common pumice types include: (1) “frothy” pumice (highly porous with large, sub-rounded vesicles), which form 5–10 vol.% of the coarsest pyroclastic flow deposits, (2) dominantly “microvesicular” and systematically crystal-poor pumices, which are found in early erupted, fine-grained pyroclastic flow units, and are characterized by many small (<50 μm in diameter) vesicles and few mm-sized, irregular voids, (3) grey or banded pumices, indicating the interaction between the rhyolite and a more mafic magma, which are found throughout the eruption sequence and display highly irregular bubble shapes. Except for the grey-banded pumices, all three other types are compositionally identical and were generated synchronously as they are found in the same pyroclastic units. They, therefore, record different conditions in the volcanic conduit leading to variable bubble nucleation, growth and coalescence. A total of 74 pumice samples have been characterized using thin section observation, SEM imagery, porosimetry, and permeametry. We show that the four pumice types have distinct total and connected porosity, tortuosity and permeability. Grey-banded pumices show large variations in petrophysical characteristics as a response to mingling of two different magmas. The microvesicular, crystal-poor, pumices have a bimodal bubble size distribution, interpreted as reflecting an early heterogeneous bubble nucleation event followed by homogeneous bubble nucleation close to fragmentation. Finally, the significant differences in porosity, tortuosity and permeability in compositionally identical tube and frothy pumices are the result of variable shear rates in different parts of the conduit. Differential shear rates may be the result of either: (1) pure shear, inducing a vertical progression from frothy to tube and implying a relatively thick fragmentation zone to produce both types of pumices at the same time or (2) localized simple shear, inducing strongly tubular vesicles along the wall and near-spherical bubbles in the centre of the conduit and not necessarily requiring a thick fragmentation zone.     

Isotopic and geochemical characterization of salinization in the shallow aquifers of a reclaimed subsiding zone: The southern Venice Lagoon coastland

The coastal plain bordering the southern Venice Lagoon is a reclaimed lowland characterized by high subsidence rate, and ground level and water-table depth below sea level. In this agricultural region, where the surface hydrologic network is entirely artificially controlled by irrigation/drainage canals, salinization problems have long been encountered in soils and groundwaters. Here we use isotopic and geochemical tracers to improve our understanding of the origin of salinization and mineralization of the semi-confined aquifer (0–40 m), and the freshwater inputs to this hydrological system. Water samples have been collected at different seasons in the coastal Adriatic Sea, lagoon, rivers and irrigation canals, as well as in the semi-confined aquifer at depths between 12 and 35 m (14 boreholes), and in the first confined aquifer (three boreholes drilled between 40 and 80 m depth). Stable isotopes (δ¹⁸O and δD) and conductivity profiles show that direct saline intrusion from the sea or the lagoon is observed only in a restricted coastal strip, while brackish groundwaters are found over the entire topographic and piezometric depression in the centre of the study area. Fresh groundwaters are found only in the most western zone. The sharp isotopic contrast between the western and central regions suggests disconnected hydrological circulations between these two parts of the shallow aquifer. The border between these two regions also corresponds to the limits of the most strongly subsiding zone.Our results can be interpreted in terms of a four end-member mixing scheme, involving (1) marine water from the lagoon or the open sea, (2) alpine and pre-alpine regional recharge waters carried either by the main rivers Adige, Bacchiglione and Brenta (irrigation waters) or by the regional groundwater circulation, (3) local precipitation, and (4) evaporated waters infiltrated from the surface. Infiltration from the surface is also revealed by the stratification of the electrical conductivity profiles, showing that the brackish groundwaters are overlain by a shallow layer of less saline water all over the central depression. In the first confined aquifer, the groundwaters have isotopic compositions similar to the deep groundwaters of the Venetian confined aquifers (40–400 m depth). The isotopic data and the Br/Cl ratio show that the origin of the salinization of the phreatic aquifer can be ascribed to seawater intrusion alone, with no indication of the involvement of deep brines (identified at 450 m depth) in the process.The chemical composition of the saline and brackish groundwaters is characterized by an excess of sodium and a deficit of calcium compared to conservative mixing between fresh groundwaters and seawater. This suggests that the phreatic aquifer is progressively freshening, as a consequence of the beneficial influence of the extensive irrigation/drainage network, including raised canals acting as a hydraulic barrier along the coast. This freshening tendency may have been lasting since the reclamation in the mid-twentieth century, and has probably been accelerated by the ban on groundwater abstraction since the 1970s.     

Numerical modelling of strain in lava tubes

Received: 5 November 1998 / Accepted: 17 November 1999     

Environmental controls on clay mineralogy of an Early Jurassic mudrock (Blue Lias Formation,southern England)

We present a dataset including clay mineralogy, gamma-ray spectrometry, organic matter content and magnetic susceptibility of the Hettangian to lowest Sinemurian successions of Dorset and Somerset, southern UK (Blue Lias Formation, Bristol Channel and Wessex basins). In both areas, the clay assemblages comprise predominantly detrital illite, kaolinite and illite/smectite mixed layers. Clays probably originated from the erosion of the Hercynian massifs, the relative proportions of kaolinite and illite being modulated by arid-humid climatic fluctuations. The organic matter (OM) content (types II to IV) ranges up to 12% in both areas. A clear stratigraphical trend in clay mineral assemblages is apparent in Somerset, whereas in Dorset sharp contrasts between adjacent horizons and a greater dilution by carbonate mask the long-term evolution. Correlations between both areas based on similar vertical trends in clay mineral abundance support the suggestion of a hiatus within the angulata Zone of the Dorset succession. As expected, the kaolinite/illite ratio correlates with the Th/K ratio deduced from gamma-ray spectrometry. However, significant departures from the correlation occur in OM-rich intervals, suggesting that Th may be partly adsorbed on to OM particles. Surprisingly, high magnetic susceptibility correlates with abundant kaolinite, not with Fe-rich clays, indicating either that kaolinite is accompanied by a soil-inherited magnetisable phase (possibly iron oxide) or that illite-rich rocks are more strongly diluted by carbonate than are kaolinite-rich strata.     

The petrologic evolution and pre-eruptive conditions of the rhyolitic Kos Plateau Tuff (Aegean arc)

The Kos Plateau Tuff is a large (>60 km³) and young (160 k.y.) calc-alkaline, high-SiO₂ rhyolitic ignimbrite from the active Kos-Nisyros volcanic center in the Aegean arc (Greece). Combined textural, petrological and geochemical information suggest that (1) the system evolved dominantly by crystal fractionation from (mostly unerupted) more mafic parents, (2) the magma chamber grew over ≥ 250 000 years at shallow depth (～1.5-2.5 kb) and was stored as a H₂O-rich crystalline mush close to its solidus (～670-750°C), (3) the eruption occurred after a reheating event triggered by the intrusion of hydrous mafic magma at the base of the rhyolitic mush. Rare banded pumices indicate that the mafic magma only mingled with a trivial portion of resident crystal-rich rhyolite; most of the mush was remobilized following partial melting of quartz and feldspars induced by advection of heat and volatiles from the underplated, hotter mafic influx.     

Micro- and nano-textural evidence of Ti(–Ca–Fe) mobility during fluid–rock interactions in carbonaceous lawsonite-bearing rocks from New Zealand

Understanding the mobility of chemical elements during fluid–rock interactions is critical to assess the geochemical evolution of a rock undergoing burial and metamorphism and, more generally, to constrain the geochemical budget of the subduction factory. In particular, determining the behavior and mobility of Ti in aqueous fluids constitutes a great challenge that is still under scrutiny. Here, we study plant fossils preserved in blueschist metasedimentary rocks from the Marybank Formation (New Zealand). Using scanning and transmission electron microscopies (SEM and TEM), we show that the carbonaceous material (CM) composing the fossils contains abundant nano-inclusions of Ti- and Fe-oxides. These nanocrystals are mainly anatase, rutile, and Fe–Ti oxides. The mineral composition observed within the fossils is significantly different from that detected in the surrounding rock matrix. We propose that Ti and Fe might have been mobilized by the alteration of a detrital Ti–Fe-rich protolith during an early diagenetic event under acidic and reducing conditions. Aqueous fluids rich in organic ligands released by the degradation of organic matter may have been involved. Moreover, using mass balance and petrological observations, we show that the contrasted mineralogy between the rock matrix and the fossil CM might be the consequence of the chemical isolation of fossil CM during the prograde path of the rock. Such an isolation results from the early formation of quartz and Fe-rich phyllosilicate layers enclosing the fossil as characterized by SEM and TEM investigations. Overall, this study shows that investigating minerals associated with CM down to the nanometer scale in metamorphic rocks can provide a precious record of early prograde geochemical conditions.     

Experimental Investigation of Turbulent Momentum Transfer in a Neutral Boundary Layer over a Rough Surface

The turbulent characteristics of the neutral boundary layer developing over rough surfaces are not well predicted with operational weather-forecasting models. The problem is attributed to inadequate mixing-length models, to the anisotropy of the flow and to a lack of controlled experimental data against which to validate numerical studies. Therefore, in order to address directly the modelling difficulties for the development of a neutral boundary layer over rough surfaces, and to investigate the turbulent momentum transfer of such a layer, a set of hydraulic flume experiments were carried out. In the experiments, the mean and turbulent quantities were measured by a particle image velocimetry (PIV) technique. The measured velocity variances and fluxes \({(\overline{{u_{i}^{\prime}}{u_{j}^{\prime}}})}\) in longitudinal vertical planes allowed the vertical and longitudinal gradients (?/?z and ?/?x) of the mean and turbulent quantities (fluxes, variances and third-order moments) to be evaluated and the terms of the evolution equations for ?e/?t, \({\partial \overline{u^{\prime 2}}/\partial t}\), \({\partial \overline{w^{\prime 2}}/\partial t}\) and \({\partial \overline{{u^{\prime}}{w^{\prime}}}/\partial t}\) to be quantified, where e is the turbulent kinetic energy. The results show that the pressure-correlation terms allow the turbulent energy to be transferred equitably from \({\overline{{u^{\prime}}^{2}}}\) to \({\overline{{w^{\prime}}^{2}}}\). It appears that the repartition between the constitutive terms of the budget of e, \({\overline{{u^{\prime}}^{2}}}\), \({\overline{{w^{\prime}}^{2}}}\) and \({\overline{{u^{\prime}}{w^{\prime}}}}\) is not significantly affected by the development of the rough neutral boundary layer. For the whole evolution, the transfers of energy are governed by the same terms that are also very similar to the smooth-wall case. The PIV measurements also allowed the spatial integral scales to be computed directly and to be compared with the dissipative and mixing length scales, which were also computed from the data.     

Climate change impacts on the energy system: a review of trends and gaps

Major transformation of the global energy system is required for climate change mitigation. However, energy demand patterns and supply systems are themselves subject to climate change impacts. These impacts will variously help and hinder mitigation and adaptation efforts, so it is vital they are well understood and incorporated into models used to study energy system decarbonisation pathways. To assess the current state of understanding of this topic and identify research priorities, this paper critically reviews the literature on the impacts of climate change on the energy supply system, summarising the regional coverage of studies, trends in their results and sources of disagreement. We then examine the ways in which these impacts have been represented in integrated assessment models of the electricity or energy system.Studies tend to agree broadly on impacts for wind, solar and thermal power stations. Projections for impacts on hydropower and bioenergy resources are more varied. Key uncertainties and gaps remain due to the variation between climate projections, modelling limitations and the regional bias of research interests. Priorities for future research include the following: further regional impact studies for developing countries; studies examining impacts of the changing variability of renewable resources, extreme weather events and combined hazards; inclusion of multiple climate feedback mechanisms in IAMs, accounting for adaptation options and climate model uncertainty.     

Isotope hydrology of a tropical coffee agroforestry watershed: Seasonal and event‐based analyses

Stable isotope variations are extremely useful for flow partitioning within the hydrologic cycle but remain poorly understood throughout the tropics, particularly in watersheds with rapidly infiltrating soils, such as Andisols in Central America. This study examines the fluctuations of stable isotope ratios (δ¹⁸O and δ²H) in the hydrologic components of a tropical coffee agroforestry watershed (~1 km²) with Andisol soils in Costa Rica. Samples were collected in precipitation, groundwater, springs, and stream water over 2 years. The local meteoric water line for the study site was δ²H = 8.5 δ¹⁸O + 18.02 (r² = 0.97, n = 198). The isotope ratios in precipitation exhibited an enriched trend during the dry season and a notable depletion at the beginning of the wet season. The δ¹⁸O compositions in groundwater (average = ?6.4‰, σ = 0.7) and stream water (average = ?6.7‰, σ = 0.6) were relatively stable over time, and both components exhibited more enriched values in 2013, which was the drier year. No strong correlation was observed between the isotope ratios and the precipitation amount at the event or daily time‐step, but a correlation was observed on a monthly scale. Stream water and base flow hydrograph separations based on isotope end‐member estimations showed that pre‐event water originating from base flow was prevalent. However, isotope data indicate that event water originating from springs appears to have been the primary driver of initial rises in stream flow and peak flows. These results indicate that isotope sampling improves the understanding of water balance components, even in a tropical humid location, where significant variations in rainfall challenge current modelling efforts. Further research using fine‐scale hydrometric and isotopic data would enhance understanding the processes driving spring flow generation in watersheds.