Applying Fractal Dimensions and Energy-Budget Analysis to Characterize Fracturing Processes During Magma Migration and Eruption: 2011–2012 El Hierro (Canary Islands) Submarine Eruption

The impossibility of observing magma migration inside the crust obliges us to rely on geophysical data and mathematical modelling to interpret precursors and to forecast volcanic eruptions. Of the geophysical signals that may be recorded before and during an eruption, deformation and seismicity are two of the most relevant as they are directly related to its dynamic. The final phase of the unrest episode that preceded the 2011–2012 eruption on El Hierro (Canary Islands) was characterized by local and accelerated deformation and seismic energy release indicating an increasing fracturing and a migration of the magma. Application of time varying fractal analysis to the seismic data and the characterization of the seismicity pattern and the strain and the stress rates allow us to identify different stages in the source mechanism and to infer the geometry of the path used by the magma and associated fluids to reach the Earth’s surface. The results obtained illustrate the relevance of such studies to understanding volcanic unrest and the causes that govern the initiation of volcanic eruptions.     

Constraining Holocene sea levels using U‐Th ages of phreatic overgrowths on speleothems from coastal caves in Mallorca (Western Mediterranean)

Phreatic overgrowths on speleothems (POS) are carbonate formations deposited at the water table of caves in unique karstic coastal settings having morphologies that can be directly related to sea level at the time of formation. The U‐Th ages of calcite and aragonite overgrowths collected from the modern water table in coastal caves on Mallorca (Cova de Cala Varques A and Cova des Pas de Vallgornera) were determined using high‐precision MC‐ICPMS techniques. U‐Th ages indicate that phreatic carbonate deposition occurred between ca 2·8 and at least 0·6 ka BP and are in accord with an archeologically estimated age of 3·7–3·0 ka BP for a drowned prehistoric construction at a depth of 1 m below current sea level in a cave from the same area. Speleothem δ¹³C and δ¹⁸O and chemical composition of cave pools provide supportive evidence that POS reflect mixing between seawater and brackish water table. Copyright © 2010 John Wiley & Sons, Ltd.     

Methodology for the computation of volcanic susceptibility: An example for mafic and felsic eruptions on Tenerife (Canary Islands)

A new method to calculate volcanic susceptibility, i.e. the spatial probability of vent opening, is presented. Determination of volcanic susceptibility should constitute the first step in the elaboration of volcanic hazard maps of active volcanic fields. Our method considers different criteria as possible indicators for the location of future vents, based on the assumption that these locations should correspond to the surface expressions of the most likely pathways for magma ascent. Thus, two groups of criteria have been considered depending on the time scale (short or long term) of our approach. The first one accounts for long-term hazard assessment and corresponds to structural criteria that provide direct information on the internal structure of the volcanic field, including its past and present stress field, location of structural lineations (fractures and dikes), and location of past eruptions. The second group of criteria concerns to the computation of susceptibility for short term analyses (from days to a few months) during unrest episodes, and includes those structural and dynamical aspects that can be inferred from volcano monitoring. Thus, a specific layer of information is obtained for each of the criteria used. The specific weight of each criterion on the overall analysis depends on its relative significance to indicate pathways for magma ascent, on the quality of data and on their degree of confidence. The combination of the different data layers allows to create a map of the spatial probability of future eruptions based on objective criteria, thus constituting the first step to obtain the corresponding volcanic hazards map. The method has been used to calculate long-term volcanic susceptibility on Tenerife (Canary Islands), and the results obtained are also presented.     

An investigation of soil erosion and redistribution in a Mediterranean lowland agricultural catchment using caesium-137

This study sought to contribute to an improved understanding of soil erosion and redistribution on Mediterranean agricultural land, where traditional soil conservation practices have been applied over millennia to provide effective protection of cultivated land. The study was undertaken in the Na Borges catchment, a groundwater-dominated lowland limestone basin (319 km2), located in the northeastern part of Mallorca, Spain. The average sediment yield from the basin, based on river sediment load data, is 1 t/km2·yr. The 137Cs technique was used to quantify soil redistribution rates over the past 40 years and to identify the key factors involved in soil erosion and redistribution processes. To estimate erosion and deposition rates and to elucidate the main factors affecting soil redistribution, samples were collected from six slope transects representative of the local land use and slope gradients and the presence or absence of soil conservation practices. A mass balance and a profile distribution conversion model were used for cultivated areas and areas of natural vegetation, respectively, to derive point estimates of the soil redistribution rates from the 137Cs inventories measured for individual soil bulk cores. In areas without soil conservation practices, the estimated mean soil erosion rates ranged from 12.7 to 26.4 t/ha·yr, which correspond to the slight and moderate erosion classes. The erosivity of Mediterranean climatic conditions combined with the influence of agricultural practices and slope gradient on soil erosion, represent the main factors responsible for the variation of soil losses documented for the cultivated land located in downslope areas, in the absence of soil conservation practices. Deposition dominated for those transects affected by soil conservation practices, with rates ranging between 18.8 and 96.6 t/ha·yr. However, this situation does not mean that soil conservation measures retain all the sediment, but rather that agriculture and urbanization (i.e. new rural paths and stone boundaries) modified the micro-topography and diverted sediment from other upslope zones towards the slopes where sampling transects were located.     

Tuning the second-generation SDI: theoretical aspects and real use cases

Spatial data infrastructure (SDI) actors have great expectations for the second-generation SDI currently under development. However, SDIs have many implementation problems at different levels that are delaying the development of the SDI framework. The aims of this article are to identify these difficulties, in the literature and based on our own experience, in order to determine how mature and useful the current SDI phenomena are. We can then determine whether a general reconceptualization is necessary or rather a set of technical improvements and good practices needs to be developed before the second-generation SDI is completed. This study is based on the following aspects: metadata about data and services, data models, data download, data and processing services, data portrayal and symbolization, and mass market aspects. This work aims to find an equilibrium between user-focused geoportals and web service interconnection (the user side vs. the server side). These deep reflections are motivated by a use case in the healthcare area in which we employed the Catalan regional SDI. The use case shows that even one of the best regional SDI implementations can fail to provide the required information and processes even when the required data exist. Several previous studies recognize the value of applying Web 2.0 and user participation approaches but few of these studies provide a real implementation. Another objective of this work is to show that it is easy to complement the classical, international standard-based SDI with a participative Web 2.0 approach. To do so, we present a mash-up portal built on top of the Catalan SDI catalogues.     

Formation of U-depleted rhyolite from a basanite at El Hierro, Canary Islands

Phonolite and trachyte are the felsic magmas of the alkaline magma suites, which characterize the Canary Islands. The October 2011 submarine eruption off El Hierro, the westernmost island, nevertheless, produced a small volume of rhyolitic magma. The rhyolite occurred as highly vesicular, white coloured pumices enveloped in and mingled with darker coloured basanitic pumice. The basanitic pumice is relatively crystal poor with a few euhedral olivines (mostly Fo_77–79), clinopyroxenes and Fe-rich spinels, whereas very rare olivine of same composition is found together with equally rare Fe-sulphide and FeTi-rich oxides in the rhyolite. The Fe–Mg exchange equilibrium in the oxides permits to calculate an equilibrium temperature of 970–890 °C for the rhyolite, in agreement with quartz-melt equilibrium at ca. 930 °C. A striking mineralogical feature of the rhyolite is the presence of rounded to contorted grains of milky quartz, which are xenocrysts incorporated and partly dissolved into the magma. Analyses of residual volatile concentrations in the glasses show that the rhyolite melt was highly degassed, whereas the basanitic glass still has important halogen concentrations. Trace element patterns of the mafic glasses and their elevated incompatible element concentrations are typical of the western Canary Island basanites. In contrast, the trace element composition of the rhyolite shows surprisingly low concentrations for all elements except the most incompatible ones (e.g. Rb, Ba, K and Th). All other measured LILE, HFSE and REE have significantly lower concentration than the basanitic counterpart that can be explained by fractionation of accessory phases (1 % apatite, 1 % sphene and 0.1 % zircon). Surprisingly, low U concentration is presumably related to elevated oxygen fugacity in the rhyolite, causing U to be in a hexavalent state, and fluxing of F-rich gas leading to volatilization of UF₆, known to emanate at low temperature. The results suggest that a gas-rich basanitic melt remobilized a small volume of stagnant rhyolitic melt formed by incorporation of approximately 10 % quartz-rich sediment into a late differentiate of trachytic composition. Sediments at the interface of an old oceanic crust adjacent to a continental shield and younger volcanic island are likely to act as magma traps were sediment assimilation may alter the mantle-derived magma composition. Quartz assimilation thus explains the production of rhyolite magma in a volcanic island characterized by an alkaline magma series from primitive basanites to trachytes.     

Evaluation of sampling devices for the determination of polycyclic aromatic hydrocarbons in surface microlayer coastal waters

The sea surface microlayer (SML) may play an important role on the transport and fate of persistent organic pollutants in the marine environment. In order to evaluate the appropriateness of a number of sampling devices for the analysis of 14 parent polycyclic aromatic hydrocarbons (3-5 aromatic rings), marine SML waters were sampled using a glass plate, a rotating drum/roller, a metal screen and a surface slick sampler. The underlying waters were also sampled for the determination of the corresponding enrichment factors (EF = [C](microlayer)/[C](underlying water)). The EFs were phase dependent, ranging from 1 to 3 for the dissolved phase and between 4 and 7 for the particulate phase. In order to better assess the performance of the different sampling methods, in terms of phase partitioning, the truly dissolved and colloidal phases were also estimated. Generally, no significant differences were found for the enrichment factors provided by the different methods, due to the observed large variability in concentrations that can be attributed to small-scale coastal processes. However, the metal screen is recommended as the most efficient sampling method for the study of PAHs taking into account the amount of water collected versus time.     

Heterogeneous characteristics of streambed saturated hydraulic conductivity of the Touchet River,south eastern Washington,USA

Traditionally a streambed is treated as a layer of uniform thickness and low saturated hydraulic conductivity (K) in surface‐ and ground‐water studies. Recent findings have shown a high level of spatial heterogeneity within a streambed and such heterogeneity directly affects surface‐ and ground‐water exchange and can have ecological implications for biogeochemical transformations, nutrient cycling, organic matter decomposition, and reproduction of gravel spawning fish. In this study a detailed field investigation of K was conducted in two selected sites in Touchet River, a typical salmon spawning stream in arid south eastern Washington, USA. In‐stream slug tests were conducted to determine K following the Bouwer and Rice method. For the upper and lower sites, each 50 m long and 9 m wide and roughly 20 m apart, a sampling grid of 5 m longitudinally and 3 m transversely was used. The slug tests were performed for each horizontal coordinate at 0·3–0·45, 0·6–0·75, 0·9–1·05 and 1·2–1·35 m depth intervals unless a shallower impenetrable obstruction was encountered. Additionally, water levels were measured to obtain vertical hydraulic gradient (VHG) between each two adjacent depth intervals. Results indicated that K ranged over three orders of magnitude at both the upper and lower sites and differed between the two sites. At the upper site, K did not differ significantly among different depth intervals based on nonparametric statistical tests for mean, median, and empirical cumulative distribution, but the spatial pattern of K varied among different depth intervals. At the lower site, K for the 0·3–0·45 m depth interval differed statistically from those at other depth intervals, and no similar spatial pattern was found among different depth intervals. Zones of upward and downward water flow based on VHG also varied among different depth intervals, reflecting the complexities of the water flow regime. Detailed characterization of the streambed as attempted in this study should be helpful in providing information on spatial variations of streambed hydraulic properties as well as surface‐ and ground‐water interaction. Copyright © 2009 John Wiley & Sons, Ltd.