Single explosions at Stromboli in 2002: Use of clast microtextures to map physical diversity across a fragmentation zone

A combination of fortuitous eruptive and wind conditions at Stromboli volcano in September 2002 enabled the collection of samples of multiple lapilli from individual observed explosions. These ejecta present the first opportunity to analyze the vesicularity of material ejected in a single Strombolian explosion. Samples of between 40 and 92 lapilli were collected from each of six sequential explosions on 30 September 2002, and 28 lapilli were obtained from a single explosion  24 h later (1 October). Density measurements and microtextural observations show that considerable heterogeneity existed within each of the seven samples. Centimeter to millimeter size bubble-rich and bubble-poor zones are present and are, in places, mingled together. These data confirm that the shallow conduit at Stromboli is a texturally diverse environment at the instant of a single explosion, and that a similar range of heterogeneity can persist through closely-spaced sequences of explosions on a timescale of hours and, probably, days, despite the dynamic processes of ascent of melt and decoupled gas phases.     

Pb,Hf and Nd isotope compositions of the two Réunion volcanoes (Indian Ocean): A tale of two small-scale mantle “blobs”?

Pb, Hf, Nd and Sr isotopes of basaltic lavas from the two Réunion Island volcanoes are reported in order to examine the origin of the sources feeding these volcanoes and to detect possible changes through time. Samples, chosen to cover the whole lifetime of the two volcanoes (from 2 Ma to present), yield a chemically restricted (compared to OIB lavas) but complex distribution. Réunion plume isotopic characteristics have been defined on the basis of the composition of uncontaminated shield-building lavas from the Piton de la Fournaise volcano. The average ?_Nd, ?_Hf, ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb isotope ratios calculated for this component are + 4.4, + 9.1, 0.70411, 18.97, 15.59 and 39.03, respectively. In Pb–Pb isotope space, each volcano defines a distinct linear trend but slight variations are also detected within the various volcanic sequences. The Piton des Neiges volcano yields a distinct and significantly more scattered isotopic distribution than Piton de la Fournaise for both Pb, Hf and Nd isotope tracers. A principal component analysis of the Pb isotope data from Piton de la Fournaise reveals a major contribution of the C and EM-1 components (with a clear Dupal flavor) as main components for the modern Réunion plume. The same components have been identified for Piton des Neiges but with a stronger participation of a depleted mantle component and a weaker EM-1 contribution. The compositional change of the lavas erupted by the Piton des Neiges and Piton de la Fournaise volcanoes is attributed to the impingement of two small-scale blobs of plume material at the base of the Réunion lithosphere. Compared to other hot-spots worldwide, in particular Hawaii and Kerguelen, magmas beneath Réunion are generated from a considerably more homogeneous, compositionally more primitive plume higher in ²⁰⁶Pb. Although shallow-level contamination processes have been locally detected they did not alter significantly the composition of the plume magmas. This is tentatively attributed to mantle dynamics producing small, high-velocity blobs that ascend rapidly through the lithosphere, and to the lack of a well-developed magma chamber at depth in the lithosphere.     

The impact of storm-triggered landslides on sediment dynamics and catchment-wide denudation rates in the southern Central Range of Taiwan following the extreme rainfall event of Typhoon Morakot

Extreme erosion events can produce large short-term sediment fluxes. Such events complicate erosion rates estimated from cosmogenic nuclide concentrations in river sediment by providing sediment with a concentration different from the long-term basin average. We present a detrital ¹⁰Be study in southern Taiwan, with multiple samples obtained in a time sequence bracketing the 2009 Typhoon Morakot, to assess the impact of landslide sediment on ¹⁰Be concentrations (N¹⁰Be) in river sediment. Sediment samples were collected from 13 major basins, two or three times over the last decade, to observe the temporal variation of N¹⁰Be. Landslide inventories with time intervals of 5–6 years were used to quantify sediment flux changes. A negative correlation between N¹⁰Be and landslide areal density indicates dilution of N¹⁰Be by landslide sediment. Denudation rates estimated from the diluted N¹⁰Be can be up to three times higher than the lowest rate derived from the same basins. Observed increases imply that, 3 years after the passage of Typhoon Morakot, fluvial channels still contain a considerable amount of sediment provided by hillslope landslides during the event. However, higher N¹⁰Be in 2016 samples indicate that the contribution from landslide sediment at the sampled grain size has decreased with time. The correlation between changes in N¹⁰Be and landslide area and volume is not strong, likely resulting from the stochastic and complex nature of sediment transport. To simultaneously evaluate the volume of landslide-derived sediment and estimate the background denudation rate, associated with less impulsive sediment supply, we constructed a sediment-mixing model with the time series of N¹⁰Be and landslide inventories. The spatial pattern of background erosion rate in southern Taiwan is consistent with the regional tectonic framework, indicating that the landscape is evolving mainly in response to the tectonic forcing, and this signal is modified, but not obscured by impulsive sediment supply. © 2019 John Wiley & Sons, Ltd.     

An Analysis Platform for Multiscale Hydrogeologic Modeling with Emphasis on Hybrid Multiscale Methods

One of the most significant challenges faced by hydrogeologic modelers is the disparity between the spatial and temporal scales at which fundamental flow, transport, and reaction processes can best be understood and quantified (e.g., microscopic to pore scales and seconds to days) and at which practical model predictions are needed (e.g., plume to aquifer scales and years to centuries). While the multiscale nature of hydrogeologic problems is widely recognized, technological limitations in computation and characterization restrict most practical modeling efforts to fairly coarse representations of heterogeneous properties and processes. For some modern problems, the necessary level of simplification is such that model parameters may lose physical meaning and model predictive ability is questionable for any conditions other than those to which the model was calibrated. Recently, there has been broad interest across a wide range of scientific and engineering disciplines in simulation approaches that more rigorously account for the multiscale nature of systems of interest. In this article, we review a number of such approaches and propose a classification scheme for defining different types of multiscale simulation methods and those classes of problems to which they are most applicable. Our classification scheme is presented in terms of a flowchart (Multiscale Analysis Platform), and defines several different motifs of multiscale simulation. Within each motif, the member methods are reviewed and example applications are discussed. We focus attention on hybrid multiscale methods, in which two or more models with different physics described at fundamentally different scales are directly coupled within a single simulation. Very recently these methods have begun to be applied to groundwater flow and transport simulations, and we discuss these applications in the context of our classification scheme. As computational and characterization capabilities continue to improve, we envision that hybrid multiscale modeling will become more common and also a viable alternative to conventional single‐scale models in the near future.     

The effect of riparian tree roots on the mass‐stability of riverbanks

Plants interact with and modify the processes of riverbank erosion by altering bank hydrology, flow hydraulics and bank geotechnical properties. The physically based slope stability model GWEDGEM was used to assess how changes in bank geotechnical properties due to the roots of native Australian riparian trees affected the stability of bank sections surveyed along the Latrobe River. Modelling bank stability against mass failure with and without the reinforcing effects of River Red Gum (Eucalyptus camaldulensis) or Swamp Paperbark (Melaleuca ericifolia) indicates that root reinforcement of the bank substrate provides high levels of bank protection. The model indicates that the addition of root reinforcement to an otherwise unstable bank section can raise the factor of safety (F _s) from F _s = 1·0 up to about F _s = 1·6. The addition of roots to riverbanks improves stability even under worst‐case hydrological conditions and is apparent over a range of bank geometries, varying with tree position. Trees growing close to potential failure plane locations, either low on the bank or on the floodplain, realize the greatest bank reinforcement. Copyright © 2000 John Wiley & Sons, Ltd.     

Site-dependent models of earthquake ground motion

Random vibration analyses of structural systems subjected to seismic loading are dependent upon the characterization of earthquake ground motion as a stochastic process. The response of structural systems to earthquakes is dependent strongly on the local geological conditions, which should be incorporated into seismological models of ground motion. In the study presented herein, three previously developed ground-motion models are adapted to incorporate site-dependent characteristics. Records obtained from two recording stations in California are used as a basis for the ground-motion models. Single-degree-of-freedom (SDOF) oscillators are subjected to ensembles of accelerograms generated from these models, and both elastic and inelastic response are considered. Response statistics are compared to those generated by the analysis of structural response to ensembles of recorded motion from the two sites. The important features of the ground motion for effective reproduction of response statistics are identified, and observations are made on the sensitivity of specific response parameters to site-dependent characteristics of the ground motion.     

Comparison of amplitude decay rates in reflection,refraction, and local earthquake records

Three types of seismic data recorded near Coalinga, California were analyzed to study the behavior of scattered waves: 1) aftershocks of the May 2, 1983 earthquake, recorded on verticalcomponent seismometers deployed by the USGS; 2) regional refraction profiles using large explosive sources recorded on essentially the same arrays above; 3) three common-midpoint (CMP) reflection surveys recorded with vibrator sources over the same area. Records from each data set were bandpassed filtered into 5 Hz wide passbands (over the range of 1–25 Hz), corrected for geometric spreading, and fit with an exponential model of amplitude decay. Decay rates were expressed in terms of inverse codaQ (Q _c ^–1 ).Q _c ^–1 values for earthquake and refraction data are generally comparable and show a slight decrease with increasing frequency. Decay rates for different source types recorded on proximate receivers show similar results, with one notable exception. One set of aftershocks shows an increase ofQ _c ^–1 with frequency.Where the amplitude decay rates of surface and buried sources are similar, the coda decay results are consistent with other studies suggesting the importance of upper crustal scattering in the formation of coda. Differences in the variation ofQ _c ^–1 with frequency can be correlated with differences in geologic structure near the source region, as revealed by CMP-stacked reflection data. A more detailed assessment of effects such as the depth dependence of scattered contributions to the coda and the role of intrinsic attenuation requires precise control of source-receiver field geometry and the study of synthetic seismic data calculated for velocity models developed from CMP reflection data.     

Fluids in the lithosphere, 2. Experimental constraints on CO2 transport in dunite and quartzite at elevated P-T conditions with implications for mantle and crustal decarbonation processes

In a series of experiments at 0.5–1.3 GPa and 1050–1200°C we have monitored the transport, via crack propagation, of CO₂ into well-annealed olivine and quartz aggregates. The objectives were to determine (1) the extent and rate of fluid penetration; (2) the effect of varying both P-T conditions and microstructure; and (3) the fluid penetration pathways. Experiments on CO₂ penetration into dunite annealed in the absence of MgO indicate rapid and pervasive fluid transport on a grain-dimension scale, but a limited penetration distance ( 1 mm). Additional experiments on dunite annealed in the presence of MgO (either dispersed or present at both ends), however, resulted in CO₂ penetration that was both pervasive on the scale of individual grains and almost always completely through the 5 mm long samples. The abundance of fine (10 μm) grains in the MgO-free dunite, in contrast to the much larger grain sizes of the samples annealed with MgO present, suggests the difference in fluid penetration behavior may arise because the strength variation in dunite scales with the grain size. Effects arising from changes in olivine point defect chemistry, however, are an additional possibility. The response of synthetic quartzite to CO₂ overpressure is distinct from that of dunite: Quartzite experiences rapid and complete penetration of CO₂, via a macroscopically visible system of transgranular fractures, over the range of P-T conditions investigated.The small amount of porosity ( 2–3%) present in most rock samples fabricated for this study, lacks three-dimensional connectivity, thus precluding any enhanced fluid penetration via porous flow. Pores could possibly enhance fluid penetration as the result of a small reduction in resistance to fracture, but the probable abundance of strength-controlling flaws in natural rocks is likely to produce similar behavior.The results of our experiments on olivine and olivine + MgO suggest that the transport of pressurized CO₂ in very olivine-rich mantle environments will be pervasive on the scale of individual grains and its extent may be dependent on rock microstructure and/or crystal chemical effects. Such pervasive fluid transport, perhaps associated with magma decarbonation, may have interesting implications for both magma transport and local LREE enrichment of adjacent mantle wall-rock. The ease with which quartzite is penetrated by CO₂ at the conditions of our experiments underscores the possible role of decarbonation reactions in crustal permeability-enhancement processes.