Magma convection and mixing dynamics as a source of Ultra-Long-Period oscillations

Many volcanic eruptions are shortly preceded by injection of new magma into a pre-existing, shallow (<10 km) magma chamber, causing convection and mixing between the incoming and resident magmas. These processes may trigger dyke propagation and further magma rise, inducing long-term (days to months) volcano deformation, seismic swarms, gravity anomalies, and changes in the composition of volcanic plumes and fumaroles, eventually culminating in an eruption. Although new magma injection into shallow magma chambers can lead to hazardous event, such injection is still not systematically detected and recognized. Here, we present the results of numerical simulations of magma convection and mixing in geometrically complex magmatic systems, and describe the multiparametric dynamics associated with buoyant magma injection. Our results reveal unexpected pressure trends and pressure oscillations in the Ultra-Long-Period (ULP) range of minutes, related to the generation of discrete plumes of rising magma. Very long pressure oscillation wavelengths translate into comparably ULP ground displacements with amplitudes of order 10⁻⁴–10⁻² m. Thus, new magma injection into magma chambers beneath volcanoes can be revealed by ULP ground displacement measured at the surface.     

Identifying the spatial pattern and importance of hydro‐geomorphic drainage impairments on unpaved roads in the northeastern USA

Roads have been widely studied as sources of runoff and sediment and identified as pollutant production sources to receiving waters. Despite the wealth of research on logging roads in forested, upland settings, little work has been conducted to examine the role of extensive networks of rural, low‐volume, unpaved roads on water quality degradation at the catchment scale. We studied a network of municipal unpaved roads in the northeastern US to identify the type and spatial extent of ‘hydro‐geomorphic impairments’ to water quality. We mapped erosional and depositional features on roads to develop an estimate of pollutant production. We also mapped the type and location of design interventions or best management practices (BMPs) used to improve road drainage and mitigate water quality impairment. We used statistical analyses to identify key controls on the frequency and magnitude of erosional features on the road network, and GIS to scale up from the survey results to the catchment scale to identify the likely importance of unpaved roads as a pollutant source in this setting. An average of 21 hydro‐geomorphic impairments were mapped per kilometer of road, averaging 0.3 m³ in volume. Road gradient and slope position were key controls on the occurrence of these features. The presence of BMPs effectively reduced erosion frequency. Scaled up to the watershed and using a conservative estimate of road–stream connectivity, our results for the Winooski River watershed in the northeastern US suggest that roughly 16% and 6% of the average annual sediment and phosphorus flux, respectively, of the Winooski River may be derived from unpaved roads. Our study identifies an under‐appreciated source of water quality degradation in rural watersheds, provides insights into identifying ‘hot spots’ of pollutant production associated with these networks, and points to effectiveness of design interventions in mitigating these adverse impacts on water quality. Copyright © 2017 John Wiley & Sons, Ltd.     

Reconstruction of the eruptive activity on the NE sector of Stromboli volcano: timing of flank eruptions since 15 ka

A multidisciplinary geological and compositional investigation allowed us to reconstruct the occurrence of flank eruptions on the lower NE flank of Stromboli volcano since 15 ka. The oldest flank eruption recognised is Roisa, which occurred at ~15 ka during the Vancori period, and has transitional compositional characteristics between the Vancori and Neostromboli phases. Roisa was followed by the San Vincenzo eruption that took place at ~12 ka during the early stage of Neostromboli period. The eruptive fissure of San Vincenzo gave rise to a large scoria cone located below the village of Stromboli, and generated a lava flow, most of which lies below sea level. Most of the flank eruptions outside the barren Sciara del Fuoco occurred in a short time, between ~9 and 7 ka during the Neostromboli period, when six eruptive events produced scoria cones, spatter ramparts and lava flows. The Neostromboli products belong to a potassic series (KS), and cluster in two differently evolved groups. After an eruptive pause of ~5,000 years, the most recent flank eruption involving the NE sector of the island occurred during the Recent Stromboli period with the formation of the large, highly K calc-alkaline lava flow field, named San Bartolo. The trend of eruptive fissures since 15 ka ranges from N30°E to N55°E, and corresponds to the magma intrusions radiating from the main feeding system of the volcano.     

New evidence of mantle heterogeneity beneath the Hyblean Plateau (southeast Sicily,Italy) as inferred from noble gases and geochemistry of ultramafic xenoliths

We analyzed major and trace elements, Sr and Nd isotopes in ultramafic xenoliths in Miocenic age Hyblean diatremes, along with noble gases of CO₂-rich fluid inclusions hosted in the same products. The xenoliths consist of peridotites and pyroxenites, which are considered to be derived from the upper mantle. Although the mineral assemblage of peridotites and their whole-rock abundance of major elements (e.g., Al₂O₃ = 0.8–1.5 wt.%, TiO₂ = 0.03–0.08 wt.%) suggest a residual character of the mantle, a moderate enrichment in some incompatible elements (e.g., La_N/Yb_N = 9–14) highlights the presence of cryptic metasomatic events. In this context a deep silicate liquid is considered the metasomatizing agent, which is consistent with the occurrence of pyroxenites as veins in peridotites. Both the Zr/Nb and ¹⁴³Nd/¹⁴⁴Nd ratios of the investigated samples reveal two distinct compositional groups: (1) peridotites with Zr/Nb ≈ 4 and ¹⁴³Nd/¹⁴⁴Nd ≈ 0.5129, and (2) pyroxenites with Zr/Nb ≈ 20 and ¹⁴³Nd/¹⁴⁴Nd ≈ 0.5130. The results of noble-gas analyses also highlight the difference between the peridotite and pyroxenite domains. Indeed, the ³He/⁴He and ⁴He/⁴⁰Ar* ratios measured in the fluid inclusions of peridotites (respectively 7.0–7.4 ± 0.1 Ra and 0.5–8.2, where Ra is the atmospheric ³He/⁴He ratio of 1.38 × 10^? 6) were on average lower than those for the pyroxenites (respectively 7.2–7.6 Ra and 0.62–15). This mantle heterogeneity is interpreted as resulting from a mixing between two end-members: (1) a peridotitic layer with ³He/⁴He ≈ 7 Ra and ⁴He/⁴⁰Ar* ≈ 0.4, which is lower than the typical mantle ratio (~ 1–4) probably due to melt extraction events, and (2) metasomatizing mafic silicate melts that gave rise to pyroxenites characterized by ³He/⁴He ≈ 7.6 Ra, with a variable ⁴He/⁴⁰Ar* due to degassing processes connected with the ascent of magma at different levels in the peridotite wall rock. The complete geochemical data set also suggests two distinct mantle sources for the xenolithic groups highlighted above: (1) a HIMU (high-μ)-type source for the peridotites and (2) a DM (depleted mantle)-type source for the pyroxenites.     

Forecasting vertical ground surface movement from shrinking/swelling soils with artificial neural networks

Artificial neural networks (ANNs) are used to estimate vertical ground surface movement when soils expand and contract due to changes in soil moisture content caused by changing climate conditions. Several counterpropagation ANN test cases were investigated to map climate data (i.e. temperature and rainfall) to vertical ground surface movement at field sites in Texas and Australia. Three of the four ANN test cases use a historical time series of climate data to forecast ground surface elevation relative to a specified datum. The fourth ANN test case predicts the rate of ground surface movement, and requires post‐processing of the predicted rates to calculate ground surface elevation relative to a specified datum. The counterpropagation network has demonstrated a successful mapping of temperature and rainfall data to vertical ground surface movement at a field site when it is trained with a subset of data from the same field site (test cases 1 and 2). The results of training an ANN on one field site and testing it on another field site (test cases 3 and 4) demonstrate the ability of the ANN to capture trends in vertical ground surface movement. When compared with the predictions from a physics‐based method (shrink test‐water content method) that requires measurements/estimates of changes in soil water content, the ANN‐based predictions (based on climatic changes) captured the trends in the field measurements of shrinking–swelling soil surface movements equally well. These findings are promising and merit further investigation with data from additional field sites. Copyright © 2007 John Wiley & Sons, Ltd.     

On the spatio-temporal distribution of M 7.0+ worldwide seismicity with a non-parametric statistics

The aim of this paper is to provide some constrains on the time behavior of earthquake generation mechanism, through the usage of a non-parametric statistics that leads up to the empirical estimation of the hazard function. The results indicate that the most characterizing temporal feature for large (M 7.0+) worldwide shallow earthquake occurrence is a clustering lasting few years, indicating that the probability of earthquake occurrence is higher immediately after the occurrence of an event. After that, the process becomes almost time independent, as in a Poisson process. Remarkably, this time clustering is very similar to what previously found for different spatio-magnitude windows, and it does not seem to depend on the tectonic style of the region. This may support the hypothesis of an universal law for earthquake occurrence.