Mixing and differentiation in the Oruanui rhyolitic magma, Taupo, New Zealand: evidence from volatiles and trace elements in melt inclusions

Large pyroclastic rhyolites are snapshots of evolving magma bodies, and preserved in their eruptive pyroclasts is a record of evolution up to the time of eruption. Here we focus on the conditions and processes in the Oruanui magma that erupted at 26.5 ka from Taupo Volcano, New Zealand. The 530 km³ (void-free) of material erupted in the Oruanui event is comparable in size to the Bishop Tuff in California, but differs in that rhyolitic pumice and glass compositions, although variable, did not change systematically with eruption order. We measured the concentrations of H₂O, CO₂ and major and trace elements in zoned phenocrysts and melt inclusions from individual pumice clasts covering the range from early to late erupted units. We also used cathodoluminescence imaging to infer growth histories of quartz phenocrysts. For quartz-hosted inclusions, we studied both fully enclosed melt inclusions and reentrants (connecting to host melt through a small opening). The textures and compositions of inclusions and phenocrysts reflect complex pre-eruptive processes of incomplete assimilation/partial melting, crystallization differentiation, magma mixing and gas saturation. ‘Restitic’ quartz occurs in seven of eight pumice clasts studied. Variations in dissolved H₂O and CO₂ in quartz-hosted melt inclusions reflect gas saturation in the Oruanui magma and crystallization depths of ∼3.5–7 km. Based on variations of dissolved H₂O and CO₂ in reentrants, the amount of exsolved gas at the beginning of eruption increased with depth, corresponding to decreasing density with depth. Pre-eruptive mixing of magma with varying gas content implies variations in magma bulk density that would have driven convective mixing. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

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.     

Geographic and temporal trends in proboscidean and human radiocarbon histories during the late Pleistocene

The causes of large animal extinctions at the end of the Pleistocene remain a hotly debated topic focused primarily on the effects of human over hunting and climate change. Here we examine multiple, large radiocarbon data sets for humans and extinct proboscideans and explore how variation in their temporal and geographic distributions were related prior to proboscidean extinction. These data include 4532 archaeological determinations from Europe and Siberia and 1177 mammoth and mastodont determinations from Europe, Siberia, and North America. All span the period from 45,000 to 12,000 calendar years BP. We show that while the geographic ranges of dated human occupations and proboscidean remains overlap across the terminal Pleistocene of the Old World, the two groups remain largely segregated and increases in the frequency of human occupations do not coincide with declines in proboscidean remains. Prior to the Last Glacial Maximum (LGM; ca 21,000 years BP), archaeological ¹⁴C determinations increase slightly in frequency worldwide while the frequency of dated proboscidean remains varies depending on taxon and location. After the LGM, both sympatric and allopatric groups of humans and proboscideans increase sharply as climatic conditions ameliorate. Post-LGM radiocarbon frequencies among proboscideans peak at different times, also depending upon taxon and location. Woolly mammoths in Beringia reach a maximum and then decline beginning between 16,000 and 15,500 years BP, woolly mammoths in Europe and Siberia ca 14,500 and 13,500 BP, and Columbian mammoth and American mastodont only after 13,000 BP. Declines among woolly mammoths appear to coincide with the restructuring of biotic communities following the Pleistocene–Holocene transition.     

Gold ore-forming fluids of the Tanami region, Northern Australia

Fluid inclusion studies have been carried out on major gold deposits and prospects in the Tanami region to determine the compositions of the associated fluids and the processes responsible for gold mineralization. Pre-ore, milky quartz veins contain only two-phase aqueous inclusions with salinities ≤19 wt% NaCl eq. and homogenization temperatures that range from 110 to 410°C. In contrast, the ore-bearing veins typically contain low to moderate salinity (<14 wt% NaCl eq.), H₂O + CO₂ ± CH₄ ± N₂-bearing fluids. The CO₂-bearing inclusions coexist with two-phase aqueous inclusions that exhibit a wider range of salinities (≤21 wt% NaCl eq.). Post-ore quartz and carbonate veins contain mainly two-phase aqueous inclusions, with a last generation of aqueous inclusions being very CaCl₂-rich. Salinities range from 7 to 33 wt% NaCl eq. and homogenization temperatures vary from 62 to 312°C. Gold deposits in the Tanami region are hosted by carbonaceous or iron-rich sedimentary rocks and/or mafic rocks. They formed over a range of depths at temperatures from 200 to 430°C. The Groundrush deposit formed at the greatest temperatures and depths (260–430°C and ≤11 km), whereas deposits in the Tanami goldfield formed at the lowest temperatures (≥200°C) and at the shallowest depths (1.5–5.6 km). There is also evidence in the Tanami goldfield for late-stage isothermal mixing with higher salinity (≤21 wt% NaCl eq.) fluids at temperatures between 100 and 200°C. Other deposits (e.g., The Granites, Callie, and Coyote) formed at intermediate depths and at temperatures ranging from 240 to 360°C. All ore fluids contained CO₂ ± N₂ ± CH₄, with the more deeply formed deposits being enriched in CH₄ and higher level deposits being enriched in CO₂. Fluids from deposits hosted mainly by sedimentary rocks generally contained appreciable quantities of N₂. The one exception is the Tanami goldfield, where the quartz veins were dominated by aqueous inclusions with rare CO₂-bearing inclusions. Calculated δ ¹⁸O values for the ore fluids range from 3.8 to 8.5‰ and the corresponding δD values range from −89 to −37‰. Measured δ ¹³C values from CO₂ extracted from fluid inclusions ranged from −5.1 to −8.4‰. These data indicate a magmatic or mixed magmatic/metamorphic source for the ore fluids in the Tanami region. Interpretation of the fluid inclusion, alteration, and structural data suggests that mineralization may have occurred via a number of processes. Gold occurs in veins associated with brittle fracturing and other dilational structures, but in the larger deposits, there is also an association with iron-rich rocks or carbonaceous sediments, suggesting that both structural and chemical controls are important. The major mineralization process appears to be boiling/effervescence of a gas-rich fluid, which leads to partitioning of H₂S into the vapor phase resulting in gold precipitation. However, some deposits also show evidence of desulfidation by fluid–rock interaction and/or reduction of the ore-fluid by fluid mixing. These latter processes are generally more prevalent in the higher crustal-level deposits.     

Characterizing hydrological processes within kettle holes using stable water isotopes in the Uckermark of northern Brandenburg,Germany

Understanding the hydrologic connectivity between kettle holes and shallow groundwater, particularly in reaction to the highly variable local meteorological conditions, is of paramount importance for tracing water in a hydro(geo)logically complex landscape and thus for integrated water resource management. This article is aimed at identifying the dominant hydrological processes affecting the kettle holes' water balance and their interactions with the shallow groundwater domain in the Uckermark region, located in the north-east of Germany. For this reason, based on the stable isotopes of oxygen (δ¹⁸O ) and hydrogen (δ²H ), an isotopic mass balance model was employed to compute the evaporative loss of water from the kettle holes from February to August 2017. Results demonstrated that shallow groundwater inflow may play the pivotal role in the processes taking part in the hydrology of the kettle holes in the Uckermark region. Based on the calculated evaporation/inflow (E/I) ratios, most of the kettle holes (86.7%) were ascertained to have a partially open, flow-through-dominated system. Moreover, we identified an inverse correlation between E/I ratios and the altitudes of the kettle holes. The same holds for electrical conductivity (EC) and the altitudes of the kettle holes. In accordance with the findings obtained from this study, a conceptual model explaining the interaction between the shallow groundwater and the kettle holes of Uckermark was developed. The model exhibited that across the highest altitudes, the recharge kettle holes are dominant, where a lower ratio of E/I and a lower EC was detected. By contrast, the lowest topographical depressions represent the discharge kettle holes, where a higher ratio of E/I and EC could be identified. The kettle holes existing in between were categorized as flow-through kettle holes through which the recharge takes place from one side and discharge from the other side.     

Chloride source delineation in an urban-agricultural watershed: Deicing agents versus agricultural contributions

Analyses (n = 525) of chloride (Cl⁻), bromide (Br⁻), nitrate as nitrogen (NO₃-N), sodium (Na⁺), calcium (Ca²⁺) and potassium (K⁺) in stream water, tile-drain water and groundwater were conducted in an urban-agricultural watershed (10% urban/impervious, 87% agriculture) to explore potential differences in the signature of Cl⁻ originating from an urban source as compared with an agricultural source. Only during winter recharge events did measured Cl⁻ concentrations exceed the 230 mg/L chronic threshold. At base flow, nearly all surface water and tile water samples had Cl⁻ concentrations above the calculated background threshold of 18 mg/L. Mann–Whitney U tests revealed ratios of Cl⁻ to Br⁻ (p = .045), to NO₃-N (p < .0001), to Ca²⁺ (p < .0001), and to Na⁺ (p < .0001) to be significantly different between urban and agricultural waters. While Cl⁻ ratios indicate that road salt was the dominant source of Cl⁻ in the watershed, potassium chloride fertilizer contributed as an important secondary source. Deicing in watersheds where urban land use is minimal had a profound impact on Cl⁻ dynamics; however, agricultural practices contributed Cl⁻ year-round, elevating stream base flow Cl⁻ concentrations above the background level.     

Granitic plutonism as an indicator of microplates in the Palaeozoic of central and eastern Maine

Oxygen and strontium isotope ratios have been used to characterize source regions for granitic magmas for a transect across the northern Appalachian orogen in central and eastern Maine. The northwestern plutons (Katahdin and Seboeis) have δ¹⁸O values of 10.3–13.3 and initial ⁸⁷Sr/⁸⁶Sr ratios of 0.7083 and 0.7066, respectively. The central plutons (Bottle Lake and Center Pond) have lower δ¹⁸O values (8.2–9.9) and initial ⁸⁷Sr/⁸⁶Sr ratios (0.7043–0.7055). The southeastern plutons (Lucerne and Deblois) have δ¹⁸O values (9.0–11.0) but initial ⁸⁷Sr/⁸⁶Sr ratios (0.7077 and 0.7041, respectively) which are intermediate between the northwestern and central plutons.Source models derived from these results and other petrological and geochemical data reflect the juxtaposition of discrete source regions by transcurrent faulting, which may be related to oblique plate motions. This model illustrates the importance of microplate accretion in the Palaeozoic history of the northern Appalachian orogen.     

Prevalence and distribution of fecal indicator organisms in South Florida beach sand and preliminary assessment of health effects associated with beach sand exposure

Fecal indicator levels in nearshore waters of South Florida are routinely monitored to assess microbial contamination at recreational beaches. However, samples of sand from the surf zone and upper beach are not monitored which is surprising since sand may accumulate and harbor fecal-derived organisms. This study examined the prevalence of fecal indicator organisms in tidally-affected beach sand and in upper beach sand and compared these counts to levels in the water. Since indicator organisms were statistically elevated in sand relative to water, the study also considered the potential health risks associated with beach use and exposure to sand. Fecal coliforms, Escherichia coli, enterococci, somatic coliphages, and F(+)-specific coliphages were enumerated from sand and water at three South Florida beaches (Ft. Lauderdale Beach, Hollywood Beach, and Hobie Beach) over a 2-year period. Bacteria were consistently more concentrated in 100g samples of beach sand (2-23 fold in wet sand and 30-460 fold in dry sand) compared to 100ml samples of water. Somatic coliphages were commonly recovered from both sand and water while F(+)-specific coliphages were less commonly detected. Seeding experiments revealed that a single specimen of gull feces significantly influenced enterococci levels in some 3.1m(2) of beach sand. Examination of beach sand on a micro-spatial scale demonstrated that the variation in enterococci density over short distances was considerable. Results of multiple linear regression analysis showed that the physical and chemical parameters monitored in this study could only minimally account for the variation observed in indicator densities. A pilot epidemiological study was conducted to examine whether the length of exposure to beach water and sand could be correlated with health risk. Logistic regression analysis results provided preliminary evidence that time spent in the wet sand and time spent in the water were associated with a dose-dependent increase in gastrointestinal illness.     

Estuarine shoreline processes in a dynamic low-energy system

Estuarine shorelines are often classified as low-energy coasts and are, therefore, expected to undergo little variation. Port Stephens (SE Australia) is a ria-like drowned river valley microtidal estuary located on a wave-dominated coast. The outer part of the estuary is tide-dominated and has a large shallow flood-tide delta, which is also affected by waves. The northern (predominantly low-energy) shoreline of outer Port Stephens is a continuous stretch of sand comprising areas of high mobility and areas of relative stability terminating in a western extending sand spit. This paper investigates the effects of periodic high-energy conditions during which waves penetrate into the estuary by analysing two types of storms, low to moderate (more frequent type) storms and severe to extreme (low frequency) storms. It is established that low to moderate storms cause generalised erosion over the northern shoreline. On the contrary, severe to extreme storms, while causing erosion on parts of the beach, can transport new sediment across the flood-tide delta and deposit it to build a mobile shore attached sandwave. Long-term (decadal) trends identified in the study area are in agreement with short- and medium-term results. Moving into the estuary are four complementary zones of sediment transport which include: (1) sandwave formation and westward migration; (2) a relatively stable area between the sandwave and an erosion zone; (3) an erosion zone undergoing shoreline retreat and finally (4) a depositional terminus causing westward extension of the sand spit.