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.     

Study of the Interaction Mechanism in the Biosorption of Copper(II) Ions onto Posidonia oceanica and Peat

A systematic approach was used to characterize the biosorption of copper(II) onto two biosorbents, Posidonia oceanica and peat, focusing on the interaction mechanisms, the copper(II) sorption–desorption process and the thermal behavior of the biosorbents. Sorption isotherms at pH 4–6 were obtained and the experimental data were fitted to the Langmuir model with a maximum uptake (q_max) at pH 6 of 85.78 and 49.69 mg g^?1, for P. oceanica and peat, respectively. A sequential desorption (SD) with water, Ca(NO₃)₂, and EDTA was applied to copper‐saturated biosorbents. Around 65–70% copper(II) were desorbed with EDTA, indicating that this heavy metal was strongly bound. The reversibility of copper(II) sorption was obtained by desorption with HCl and SD. Fourier transform IR spectroscopy (FTIR) analysis detected the presence of peaks associated with OH groups in aromatic and aliphatic structures, CH, CH₂, and CH₃ in aliphatic structures, COO^? and COOH groups and unsaturated aromatic structures on the surface of both biosorbents, as well as peaks corresponding to Si? O groups on the surface of peat. The results of SEM‐EDX and FTIR analysis of copper‐saturated samples demonstrated that ion exchange was one of the mechanisms involved in copper(II) retention. Thermal analysis of biosorbent samples showed that copper(II) sorption–desorption processes affected the thermal stability of the biosorbents.     

Elasticity of antigorite,seismic detection of serpentinites,and anisotropy in subduction zones

Serpentinization of the mantle wedge is an important process that influences the seismic and mechanical properties in subduction zones. Seismic detection of serpentines relies on the knowledge of elastic properties of serpentinites, which thus far has not been possible in the absence of single-crystal elastic properties of antigorite. The elastic constants of antigorite, the dominant serpentine at high-pressure in subduction zones, were measured using Brillouin spectroscopy under ambient conditions. In addition, antigorite lattice preferred orientations (LPO) were determined using an electron back-scattering diffraction (EBSD) technique. Isotropic aggregate velocities are significantly lower than those of peridotites to allow seismic detection of serpentinites from tomography. The isotropic V_P/V_S ratio is 1.76 in the Voigt–Reuss–Hill average, not very different from that of 1.73 in peridotite, but may vary between 1.70 and 1.86 between the Voigt and Reuss bonds. Antigorite and deformed serpentinites have a very high seismic anisotropy and remarkably low velocities along particular directions. V_P varies between 8.9 km s^? 1 and 5.6 km s^? 1 (46% anisotropy), and 8.3 km s^? 1 and 5.8 km s^? 1 (37%), and V_S between 5.1 km s^? 1 and 2.5 km s^? 1 (66%), and 4.7 km s^? 1 and 2.9 km s^? 1 (50%) for the single-crystal and aggregate, respectively. The V_P/V_S ratio and shear wave splitting also vary with orientation between 1.2 and 3.4, and 1.3 and 2.8 for the single-crystal and aggregate, respectively. Thus deformed serpentinites can present seismic velocities similar to peridotites for wave propagation parallel to the foliation or lower than crustal rocks for wave propagation perpendicular to the foliation. These properties can be used to detect serpentinite, quantify the amount of serpentinization, and to discuss relationships between seismic anisotropy and deformation in the mantle wedge. Regions of high V_P/V_S ratios and extremely low velocities in the mantle wedge of subduction zones (down to about 6 and 3 km.s^?1 for V_P and V_S, respectively) are difficult to explain without strong preferred orientation of serpentine. Local variations of anisotropy may result from kilometer-scale folding of serpentinites. Shear wave splittings up to 1–1.5 s can be explained with moderately thick (10–20 km) serpentinite bodies.     

Montecarlo DLA-type simulations of wetting effects in fluid displacement in porous media

In this work, we report diffusion-limited aggregation (DLA)-type Montecarlo computations of a stochastic model of displacement of a viscous fluid by another that preferentially wets a porous medium, for the case when both fluids are immiscible in the absence of buoyancy forces. The model has the aim to simulate cooperative invasion processes found in experiments of immiscible wetting displacement. The model considers the nonlocal effects of the Laplacian pressure field and the capillary forces via hydrodynamic equations in the Darcy regime with a boundary condition for the pressure at the interface. The boundary condition contains two different types of disorder: the capillary term, which constitutes an additive random disorder, and a term containing an effective random surface tension, which couples to a curvature (it constitutes a multiplicative random term that carries nonlocal information of the whole pressure). We generate different displacement patterns for different setting of the parameters of the model. We analyze these patterns by studying the scaling properties of the interface that separate the two fluids and calculating the fractal dimension of the interface. The results show the existence of three distinct regimes of scaling. One regime at the smallest-length scales is due to the multiplicative random disorder together with the nonlocal coupling; it reveals itself in a roughness exponent α ≈ 0.80. Additionally, we find a DLA-type scaling regime with a roughness exponent α ≈ 0.60 at the largest scales and intermediate scaling regime with α ≈ 0.70 corresponding to invasion percolation with trapping. Each regime has definite scaling ranges that depend on the capillary number and the relative wetting tendency of the fluids. The behavior of the fractal dimensions of the interfaces of the aggregates constitutes a further confirmation of the existence of three scaling regimes and the multi-self-affinity of the perimeter of the interface boundaries.     

Stepwise selection of functional covariates in forecasting peak levels of olive pollen

High levels of airborne olive pollen represent a problem for a large proportion of the population because of the many allergies it causes. Many attempts have been made to forecast the concentration of airborne olive pollen, using methods such as time series, linear regression, neural networks, a combination of fuzzy systems and neural networks, and functional models. This paper presents a functional logistic regression model used to study the relationship between olive pollen concentration and different climatic factors, and on this basis to predict the probability of high (and possibly extreme) levels of airborne pollen, selecting the best subset of functional climatic variables by means of a stepwise method based on the conditional likelihood ratio test.     

Daily precipitation concentration in Central Coast Vietnam

Theoretical and Applied Climatology - Empirical frequency distribution of daily precipitation amounts can be fitted by a negative exponential distribution, because anywhere there are many small...     

Structure and nuclear density distribution in the cheralite—CaTh(PO4)2: studies of its behaviour under high pressure (36?GPa)

The crystal structure of the cheralite—CaTh(PO₄)₂—has been revisited by neutron diffraction and its behaviour under high pressure investigated by X-ray diffraction up to 36?GPa. The neutron diffraction data at ambient pressure gave a more accurate determination of the Ca/Th cation position than previous XRD data, taking advantage that the neutron scattering lengths of calcium and thorium are of same order of magnitude. The nuclear density distribution was also determined using the maximum entropy method (MEM) confirming that the two cations are not located at the same position in the unit cell but are slightly displaced from one another along a specific direction in order to minimize the electrostatic repulsion with the surrounding phosphorus atoms. At high pressure, the compound did not show any phase transition or amorphization. From the evolution of the unit-cell volume as a function of the pressure, the zero-pressure bulk modulus B₀ and its pressure derivative B₀ ^′ have been determined by fitting the experimental compressibility curve to the Birch–Murnaghan equation of state. The results are B₀?=?140(2) GPa and B ₀ ^′ ?=?4.4(4) GPa.