Next Article in Journal
Rain, Snow and Frozen Soil: Open Questions from a Porescale Perspective with Implications for Geohazards
Next Article in Special Issue
Magnetic Expression of Hydrothermal Systems Hosted by Submarine Calderas in Subduction Settings: Examples from the Palinuro and Brothers Volcanoes
Previous Article in Journal
Analysis of the Northern Hemisphere Atmospheric Circulation Response to Arctic Ice Reduction Based on Simulation Results
Previous Article in Special Issue
Age of Initial Submarine Volcanism in the Paleo-Tsushima Basin and Implications for Submarine Volcanism in the Opening Stage of the Japan Sea in Northern Kyushu
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Volcanic Thermohaline Fluctuations of Tagoro Submarine Volcano (El Hierro Island, Canary Islands, Spain)

by
Anna Olivé Abelló
1,2,
Beatriz Vinha
2,
Francisco Machín
3,
Francesco Zerbetto
4,
Evangelos Bakalis
4,* and
Eugenio Fraile-Nuez
2,*
1
Departament d’Oceanografia Física i Tecnològica, Institut de Ciències del Mar, CSIC, 08003 Barcelona, Spain
2
Centro Oceanográfico de Canarias, Instituto Español de Oceanografía (IEO), Consejo Superior de Investigaciones Científicas (CSIC), 38180 Santa Cruz de Tenerife, Spain
3
Departamento de Física, Universidad de Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
4
Dipartimento di Chimica “G. Ciamician”, Università di Bologna, 40126 Bologna, Italy
*
Authors to whom correspondence should be addressed.
Geosciences 2021, 11(9), 374; https://doi.org/10.3390/geosciences11090374
Submission received: 29 June 2021 / Revised: 26 August 2021 / Accepted: 28 August 2021 / Published: 4 September 2021
(This article belongs to the Special Issue Present and Past Submarine Volcanic Activity)

Abstract

:
Temperature and conductivity fluctuations caused by the hydrothermal emissions released during the degasification stage of the Tagoro submarine volcano (Canary Islands, Spain) have been analysed as a robust proxy for characterising and forecasting the activity of the system. A total of 21 conductivity-temperature-depth time series were gathered on a regular high-resolution grid over the main crater of Tagoro volcano. Temperature and conductivity time series, as manifestations of stochastic events, were investigated in terms of variance and analysed by the Generalised Moments Method (GMM). GMM provides the statistical moments, the structure functions of a process whose shape is an indicator of the underlying stochastic mechanisms and the state of activity of the submarine volcano. Our findings confirm an active hydrothermal process in the submarine volcano with a sub-normal behaviour resulting from anti-persistent fluctuations in time. Its hydrothermal emissions are classified as multifractal processes whose structure functions present a crossover between two time scales. In the shorter time scale, findings point to the multiplicative action of two random processes, hydrothermal vents, which carries those fluctuations driving the circulation over the crater, and the overlying aquatic environment. Given that both temperature and conductivity fluctuations are nonstationary, Tagoro submarine volcano can be characterised as an open system exchanging energy to its surroundings.

1. Introduction

Hydrothermal vents, found in submarine volcanic islands, release important emissions of hot fluids that rise in a buoyant turbulent plume over the source [1]. Those hydrothermal emissions range from intense high-temperature plumes to diffuse low temperature diluted discharges emanating directly from the seafloor [2]. The fluids mix with the surrounding seawater and rise to a level where their density matches the water outside the plume being able to alter the temperature and salinity fields specific to the vertical seawater structure [3]. These vertical motions generate a horizontal flow, resulting in an anticyclonic vortex of sufficient magnitude to trap water, minerals, tracers, and organisms in a local recirculation [4]. As a result, temperature and conductivity fields may fluctuate, influencing the local circulation. At the same time, fluctuations are the outcome of random perturbations acting on a field, noise sources caused by the short-term circulation over the volcano [5]. Decoding fluctuations with stochastic analysis tools can shed light and benefit understanding the complexity of hydrothermal emissions [6,7,8].
In the Mediterranean, a newly developed mathematical technique was applied to the time series of temperature and conductivity fluctuations, recorded by remotely operated vehicles (ROV), to understand the nature of the stochastic processes that drive the activity of volcanoes [6,7]. These studies suggest that temperature fluctuations can be used as an indicator to determine whether the volcano operates as an open or closed system, i.e., if it exchanges energy with the surroundings or not. Open stands for nonstationary temperature time series and closed for stationary ones. On the other hand, conductivity, the ability of a water sample to conduct an electrical current, fluctuations can provide information about its current state of activity. During unrest (active) degassing periods, conductivity and temperature follow a universal multifractal behaviour, i.e., a multiplicative action of random processes, with continuous energy dissipation. Otherwise, conductivity behaves as a universal multifractal but with a stationary pattern [6,7].
El Hierro is the youngest and southwesternmost island of the Canary Archipelago characterised by a three-armed rift system and originated by hotspot volcanism [9,10] (Figure 1). On 10 October 2011, a submarine volcano eruption took place 1.8 km south of the coast of El Hierro Island, forming the shallow submarine volcano Tagoro, located at 27°37.1160′ N, 017°59.4660′ W [11]. When the eruption finished in March 2012, hydrothermal manifestations were characterised by the release of hot fluids, gases, inorganic nutrients and metals as a result of the new degasification phase of the post-eruptive stage of the Tagoro submarine volcano [12]. This degasification process is exclusively associated with hydrothermalism linked with shallow submarine volcanoes, resulting in water properties changes [13,14].
Since 2011, the Tagoro submarine volcano has been deeply monitored by more than 28 different multidisciplinary cruises in order to understand submarine volcanic processes better and how the volcano evolves in terms of different physical-chemical seawater parameters and the relationship of these with the surrounding marine ecosystem [8,11,12,15,16,17,18,19,20,21,22,23,24,25,26,27]. However, there is a big gap in knowledge in terms of understanding how fluctuations in temperature and conductivity can provide information regarding the hydrothermal system and its state of activity.
The main goal of this manuscript is, using stochastic calculus tools, to assess and classify the state of activity (rest/unrest) and characterise the nature of stochastic mechanisms pertaining to around and over the main crater of the Tagoro submarine volcano in its degassing stage. This objective will be evaluated by applying the mathematical methodology proposed in [6] and using the temperature and conductivity measurements collected by a CTD sensor (Conductivity-Temperature-Depth) mounted on a rosette over the Tagoro submarine volcano. In addition, this will be the first time that a cost-effective oceanographical instrument (rosette) is implemented to register statically hydrothermal emissions in the Atlantic Ocean instead of the ROV, widely used in the Mediterranean Sea [6,26,28]. Although ROVs are a valuable tool for collecting in situ samples and recording HD images, compared with the rosette sampler, with a significantly higher measurement frequency, they are much more limited.

2. Materials and Methods

In October 2016, the multidisciplinary cruise VULCANO-II-1016 on board of R/V Ángeles Alvariño was carried out. One of the survey objectives consisted of deploying up to 21 oceanographic stations with a CTD sensor mounted on the rosette sampler in a regular 10.5-m high-resolution grid above the main crater, main crest and surroundings (Figure 1d).
Each station records a unique time series of conductivity and temperature fluctuations, referred to here with the Ec and θ parameters, between a time length from 9 to 13 min, with the exceptions of CTDs 41, 43 and 47, as a result of discarding some periods of fluctuations because the rosette was not static one meter above the seafloor (Table 1). Thus, each CTD time series was sampled with a frequency of 24 Hz or 0.0417 s distance between two consecutive measurements maintaining the rosette sampler, through the altimeter sensor, at approximately 1 m above the seabed. At the same time, the ship remained stationary on the surface through its Dynamic Positioning System. Data were acquired using an SBE 911-plus CTD equipped with dual temperature and conductivity sensors, calibrated at the SeaBird laboratory before and after the cruise, with accuracies of 0.001 °C and 0.0003 S/m, respectively.
The time series of temperature (θ) and conductivity (Ec) recorded at each location inside the grid are subject to fluctuations and have been analysed as the result of stochastic processes. In this sense, an initial evaluation of the dynamics of the hydrothermal emissions fluctuations driven by the activity of the volcano, x ( t ) , can be made by examining the time dependence of its variance, W ( t ) = x 2 ( t ) x ( t ) 2 . If the process, x ( t ) , corresponds to uncorrelated random events (normal distribution), then the variance of x ( t ) grows linearly in time (Gaussian behaviour). Any departure from linearity, i.e., anomalous behaviour, might indicate the existence of anti-correlated/correlated events and/or systems where the environment fluctuates at similar time scales as the random variable x ( t ) [29]. In the anomalous regime, the behaviour is sub-normal/super-normal for growth slower/faster than linear, and the variance reads [30]:
W ( t ) = K γ t γ Γ ( 1 + γ )
where K γ is a generalised coefficient expressed in proper units, e.g., for temperature the units are °C2 sγ, and Γ () is the Gamma function. According to the exponent γ, the stochastic processes can be classified as sub-normal (anti-persistent process) for values from 0 < γ < 1, Brownian or normal for γ = 1, super-normal (persistent process) for 1 < γ < 2, ballistic for γ = 2, and stationary for γ = 0. It is worth mentioning that Equation (1) does not apply to instantaneous changes of enormous amplitude (Lévy flights) where the first moment diverges. For discrete time series, Equation (1) is expressed as time average and reads:
W ( Δ , T ) = 1 M i = 1 M 1 T Δ n = 1 T Δ ( x i ( n + Δ ) x i ( n ) ) 2 ( 1 M i = 1 M 1 T Δ n = 1 T Δ ( x i ( n + Δ ) x i ( n ) ) ) 2
where T = N τ is the total time, N is the total number of measurements, and τ = 0.0417 s is the reciprocal of the maximum sampling rate. M represents the number of time series monitored under the same conditions, the first summation being omitted when M = 1. The parameter Δ is the lag time and plays the role of the actual time in the analysis; it takes values in the range τ Δ T / 10 . Temperature or conductivity time series are represented by x i ( n ) with n = 1,2, 3, …, N and i = θ or Ec.
Equation (2) provides an adequate description of when the process is normal, γ = 1. Any departure from normality requires the estimation of more moments, including fractional ones. We apply the Generalised Moments Methods (GMM) following these required steps:
  • Construction of N/10 new time series each of length T Δ . The new time series, y n ( Δ ) , contains the absolute change between two values of the original series that are set apart by Δ, so that y n ( Δ ) = | x ( n τ + Δ ) x ( n τ ) | for n = 1 ,   2 , , ( ( T Δ ) / τ ) and for Δ = m τ with m = 1, 2, ..., N/10 where τ and N τ / 10 define the minimum and maximum time lags.
  • Estimation of the statistical moments, order of q, for each one of the new time series, y n ( Δ ) , is carried out according to:
    ρ ( q , Δ ) = 1 T Δ n = 1 T Δ ( y n ( Δ ) ) q
    where only positive moments, q , are taken into account. Moments between 0 < q ≤ 2 are responsible for the core of the probability density function (pdf), while moments q > 2 contribute to the tails of the pdf (Bakalis et al., 2017 and the references therein).
  • The moments will scale according to:
    ρ m ( q , Δ ) Δ z ( q )
    where z ( q ) is the structure function and its shape gives information about the stochastic mechanisms that govern the process. When the structure function fits a linear form, that is, z ( q ) = H q , then there is a direct relation between the scaling exponent γ of Equation (1) and the Hurst exponent, H, γ = 2 H , being the process characterised as monofractal. The parameter H is the mean fluctuating scaling exponent. For 0 < H < 0.5 the process is anti-persistent (sub-normal), normal for H = 0.5 , and super-normal or persistent for 0.5 < H < 1 . If z ( q ) presents any convex shape, then the process is multifractal. Among multifractals, universal multifractals are ubiquitous, and their structure function reads [31]:
    z ( q ) = H q C a 1 ( q a q )
    where the Lévy index, a , 0 ≤ a ≤ 2, indicates the class to which the probability distribution belongs and shows how fast the inhomogeneity rises. Inhomogeneity refers to the degree of mixing of various stochastic mechanisms shaping the final process. For a = 0, z ( q ) = H q , the process is monofractal, and the field homogeneous. For a = 1, z ( q ) = H q C q l o g ( q ) , the process draws steps from a Cauchy–Lorentz distribution. For a = 2, z ( q ) = H q C ( q 2 q ) , the process is the multiplication of two random processes and is described by a log-normal or Kolmogorov distribution. For 1 < a < 2, the multifractal character is the multiplicative result of more than two random processes. Finally, C is the intermittency parameter that measures the mean inhomogeneity of the field and always takes positive values.

3. Results

Implementing a high-resolution grid of temporal series around and over the main crater reveals temperature and conductivity fluctuations in the whole domain, suggesting that the fields are disturbed by emissions from the volcano and may be far from a thermodynamic equilibrium. Using non-affected vertical profiles of θ (°C) and Ec (mS/cm) measured simultaneously on the submarine volcano surroundings (see Figure 1c), we obtained the average temperature and conductivity reference profiles and their standard deviation down to 127 m (shown in Figure 2). The high variations of both properties during the whole time series over the seabed, concerning the reference profile’s standard deviation, confirm the significance of the fluctuations due to the continuous release of hydrothermal emissions. These fluctuations as a function of time and their correlations for the randomly selected CTD30, one of the 21 time series collected with a time interval of 11 min and 21 s, are illustrated in detail in Figure 3.
The variances of temperature and conductivity have been obtained as a function of lag time for all CTD time series. Figure 4 shows the variance of both temperature (Figure 4a) and conductivity (Figure 4b) of the CTD30 sampling point. Both properties present a quite similar scaling with a crossover point at Δ ≈ 3.5 s, defined as the time where the scaling exponent changes. The crossover point divides the time series into two time scales or regimes: the shorter time scale corresponds to the first regime and small lag times, and the longer time scale considers all the lag times and therefore, we will refer to it as the whole regime. We fitted each curve with Equation (1), W(Δ) = bΔγ, obtaining the early (small lag times) and the long-time (all lag times) behaviour for θ and Ec. For small lag times, 0.042 < Δ 3.5   s , temperature and conductivity pose similar weak sub-normal behaviour, with scaling exponents 0.84 and 0.85, respectively. For the long-time behaviour, the whole lag times range has been considered. Again, both θ and Ec scale similarly, with values of exponents γ = 0.75 and γ = 0.77, respectively. For both small and large times, temperature and conductivity follow sub-normal behaviour, anti-persistent variations, which become stronger for longer times and lesser when scaling exponents. The absence of stationary behaviour (γ = 0) in temperature confirms the presence of a mechanism that causes the alterations (see also below), suggesting the existence of a remarkable hydrothermal vent field, consistent with what has been described so far in the area [8,21].
The rest of the CTD time-series exhibit crossover points whose positional values range from 1.8 to 5.04 s. Their variances for the whole regime are consistent with the CTD30 sub-normal behaviour, always with scaling exponent in the range 0 < γ < 1 (Table 2). In the short time scale (first regime), although most of the sampling points show a similar pattern in terms of variance and can be defined as sub-normal, the temperature variances for CTDs 37, 38, and 50 and conductivity variances for CTDs 32, 38, 40, 48, and 50 suggest in those sampled locations processes that can be characterised as Brownian (γ ≈ 1), since 0.95 ≤ γ. The latter highlights the necessity of a closer investigation of the analysis of the recorded data which extends beyond the standard variance.
Generalised moments for both temperature and conductivity measurements up to the fourth order were determined by using Equation (3) and the moments for the CTD30 time series are illustrated in Figure 5. A vertical solid line corresponds to the crossover point dividing the two time scales, and its choice is based on the turnover point of the variance. For both θ and Ec, each moment shows dependence on the time lags.
Each moment has been fitted by Equation (4), ρ ( q , Δ ) Δ z ( q ) , and the obtained exponent is equal to the value of the structure function for this specific moment. In this way, the values of z ( q ) for the moments defined in the range 0.25 q 4 will be fitted later by Equation (5) for the estimate of the structure function. For the CTD30, the structure function for each property is represented in Figure 6. The value z ( q ) has a convex shape as a function of the order of the moment, q, for both short and long times. This convexity indicates the multifractality of temperature and conductivity time series. For the first regime (short times), the Lévy stability index, a , of the structure function is equal to 2 and points to the existence of a log-normal or Kolmogorov distribution, which describes the multiplicative effect of two random processes. In the simplest scenario, the first process is likely the volcano hydrothermal vents that would drive the second one, a short-term circulation over the crater that contributes to the rise of the temperature and conductivity variables. For longer times, the process becomes even more complicated. The stability index takes values a = 1.61 and a = 1.58 for temperature and conductivity, respectively, highlighting an increase in complexity as the observation time increases. And yet, this increase in complexity may reflect the existence of more than the previously mentioned two random processes that shape the overall multifractal character. The same is true for all CTD time series analysed in the vicinity of the crater because the structure function for short times always presents a log-normal distribution with Lévy stability index equal to 2 and for long times corresponds to a universal multifractal, Equation (5), whose Lévy stability index takes value in the range 1.04 ≤ a < 1.90, see Table 2.
For short times (Figure 6), the H and C parameters of the structure function of CTD30 take values of 0.662 and 0.055 for θ and 0.669 and 0.054 for Ec, respectively. Considering the structure functions obtained for all the CTD time series for the first regime ( a = 2), listed in Table 2, the H exponent is varied between 0.540 and 0.695, implying super-normal processes if the intermittency parameter, parameter C, were zero. However, the latter is not true, and the non-zero intermittent effects establish conductivity and temperature mean fields as non-conservative. In the long-time limit, the structure-function of CTD30 provides values of H and C of 0.424 and 0.059 for θ and 0.434 and 0.060, for Ec, respectively. For long times, the Hurst exponents are lower than 0.5, implying stronger sub-normal processes and an increasing complexity reflected in the value of the alpha-stable index, which is different now for temperature and conductivity. For short times, the structure functions of both temperature and conductivity have a high similarity indicating a strong correlation of one another. Indeed, the normalised cross-correlation coefficient of the two, rθ,Ec, is higher than 0.95 for short times and decreases as the time increases (Figure 3c). Strong correlation does not necessarily imply the existence of causality between temperature and conductivity since this is the result of a broader mechanism affecting both properties equally and simultaneously, as the hydrothermal field does. Actually, by differentiating, with respect to time, the time series of temperature, Δθ, and conductivity, ΔEc, (Figure 3d,e), respectively, the normalised cross-correlation coefficient of the two, rΔθ, ΔEc, goes rapidly to zero (it is not exactly delta correlated) and the crossing of the abscissa is at 0.6 s. After that point and for a short window, before fluctuations start around zero (see Figure 3f), Δθ and ΔEc are anticorrelated to one another. These findings support the argument that the hydrothermal field is the reason for the strong correlation between temperature and conductivity since it operates as a vehicle carrying material to a sampling point.
The findings show that for both short and long time scales as defined above, the volcano area operates as an open system releasing a continuous flux of energy into the seawater, establishing thus temperature and conductivity to be non-conservative mean fields as a result of local swirls and/or vortexes. Alternatively, a non-conservative conductivity mean field might also be due to the presence of a magnetic field generated by the ions discharged from the crust of the volcanic hydrothermal chimneys to the water [6].
To understand how the environment and seawater above Tagoro main crater are affected by those mechanisms governing the stochastic process and characterising the system, the spatial distribution of CTD time series according to its stochastic behaviour in the first regime (short times) is shown in Figure 7. There is not a significant trend in stochastic processes behaviour concerning the main crater distance, which is located over and between the CTDs 30, 34, 35 and 36. Analysing together the scaling of the mean field of temperature and its level of complexity, we found that the entire area, and not only the area above the crater, presents a higher activity and non-conservative field acting as an open system far from equilibrium. Near the crater, a less non-conservative region is found, characterised by an H smaller than 0.57. Additionally, in the middle of the high-resolution grid, the conductivity field can be considered an almost conservative one, C ~ 0 (Figure 7). These findings likely reflect local structures which do not allow strong mixing with the overlying surrounding environment, or it might mirror a stronger discharge of the fluid near the crater that weakens the effect of the marine short-term circulation at the sampling point in the window of observation, or it can be a combination of the two. The overall result reduces the intermittency parameter. Located mainly in the grid edges, the CTDs 40, 48, and 50, whose short time behaviour in terms of variance points to Brownian motion, present a multifractal character as well and match with a lower intermittency parameter, C, shown in Table 2.
In the extremes of the grid, the intermittency parameter of the conductivity field increases, highlighting higher complexity and displays the opposite behaviour for the temperature field. Possible contributions to the conductivity field at each sampling point can come directly from the volcano, and/or from the crust, and/or from diffusive ions from nearby points, and/or from the drifting of physical-chemical plumes, and/or from entrainment processes. If volcanic activity were the only source that affects the conductivity field, higher values of intermittency, C, would describe a non-continuous outgassing process reminiscent of short volcano breaths followed by longer ones, whereas low or even zero values would describe a continuous flow of material whose H value provides the risk of eruption (high risk for H ~1, ballistic behaviour, and low risk for H ~0.5, random behaviour). The hypothesis of diffusivity as a source that affects the conductivity field is also discarded because it is not able to wash out local inhomogeneities (e.g., a typical diffusion coefficient of 10−9 m2/s for ions moving in water during our observations results only in a displacement of ~1.3 mm) and to create a homogeneous one in such a reduced time of observation. Dynamic drifting is always present in open oceans, and its contributions might differ at different sampling points, thus mirroring the structure of the submarine volcano. The existence of drift is compatible with the variations of the temperature field. At sampling points where intermittency is almost zero for the conductivity field, its corresponding value for the temperature field is high indicating a strong mixing process. Certainly, there are points where intermittency takes high values for both fields mirroring strong mixing/competition of various random processes at these locations.

4. Discussion

The non-conservative mean fields are rather a consequence of the action of multiplicative noise, which can affect a physical system in many ways. To mention some of them, multiplicative noise can (i) stabilise a transient state [32], (ii) affect growth rates [33], and (iii) turn stationary distributions to non-Gaussian ones [34]. Furthermore, the strength of a multiplicative noise can modify the response of a system and unravel hidden periodicities, in the case that they exist, due to stochastic resonance [35]. Moreover, it is demonstrated that the effect of the buoyant seawater plume on the surrounding environment, together with the entrainment effect from the background waters and the general convective circulation over the hydrothermal sources, could potentially induce rapid changes in the temperature and conductivity mean fields [4]. This demonstrates that those physical-chemical anomalies or organisms remain near the source around the Tagoro system. Other physical oceanic processes, such as waves, tides, internal tides, etc., influence the variability of the temperature and conductivity mean fields at the Tagoro system on a larger time scale; however, a detailed study to understand such potential processes will be considered for futures studies.
The use of a high-resolution grid (Figure 7) allows us (i) to unveil the sub-normal behaviour of both temperature and conductivity fields in the entire area, and not only inside the main crater whose alterations are detected even several meters above the seafloor, (ii) to designate potential mechanisms for the functionality of the volcano, such as intermittency describing the outgassing ions and drift terms redistributing this material in the entire area, and (iii) to characterise the volcano as an active one, since multifractal properties of the conductivity field show that the volcano continuously provides material to the marine environment, and that it yet operates as an open system because of the non-conservative temperature field. The results confirm that the Tagoro submarine volcano is an active vent system significantly altering the properties of the seawater surrounding the main crater. Such physical-chemical changes have often been documented as changes in the temperature, salinity, pH, and ORP [8] in the pH and carbonate system [12,36], in the emission of reduced species and O2 [11], in nutrient enrichment [27], and also in biological changes in terms of species composition [11,16,26], mainly related with bacterial and planktonic groups [20,22]. In agreement with the high fluctuations observed in the CTD profiles and the emerging gasses and particles documented by ROV videos, our results show that the system is far from both equilibrium and homogeneity.
The observed multifractality needs a closer look at the causes that produce it. A system driven by noise sources, multiplicative and/or additive, is usually described by a general class of stochastic differential equations (see Equation (1) of [37]). Alternatively, complex systems can be described by (i) fractional calculus and generalised fractional Langevin equations [38,39], (ii) linear Langevin equation subject to multiplicative and/or to additive noise [34], and (iii) stochastic population models [40,41,42]. A stochastic approach and description of the present system should take into account the findings of fluctuation analysis. We leave this task for forthcoming work. It is worth mentioning that the method (GMM) analytically described in the present work is not the only one that can be used for fluctuations analysis. Methods such as detrended fluctuation analysis (DFA) [43], multifractal detrended fluctuation analysis (MFDFA) [44], diffusion entropy analysis (DEA) [45], or even the classical rescaled analysis (RA) [46] can be used if the corresponding time series is stationary.
For the first time, our study applied the methodology described by [6] to a shallow submarine volcano in the Atlantic Ocean. Previous studies had used time series data collected by a remotely operated vehicle (ROV) in Avyssos (Nisyros Island) [7], and in Kolumbo (Santorini Island) [6], a shallow submarine volcano, whose last eruption happened in 1650 [47]. The results of the study [6] showed that Kolumbo is a far-from-equilibrium system following sub-normal behaviour, similar to the results obtained in the present study for the Tagoro’s hydrothermal system. For Kolumbo the crossover point, where the scaling exponent changes, is 5 s, similar to 3.5 s for the present study; however, a difference in the dynamics between the two (Kolumbo and Tagoro) might not be excluded due to closed/open sea conditions.
Overall, our results show that the GMM methodology can be applied to different systems and that the use of other oceanographical instruments for data collection is consistent with the effectiveness of the methodology. This can be useful since it allows choice of the instrument to be employed, although the ROV represent higher cost-effectiveness related to a rosette sampler since the utilisation of ROVs in oceanographical campaigns presents a notably higher cost.
The cost-effective methodology applied in this study may be exploited for forecasting future eruptions. By using this sampling methodology in time and using GMM for method analysis, we can characterise the level of activity (values of H and C for both temperature and conductivity) and monitor the evolution of such a physical system, predicting the natural hazard according to its potentiality. A combination of a high value of H and a value of C close to zero (almost homogeneous process, which means that the hydrothermal material dominates over the other processes) underlines an intense activity of the volcano, possibly according to a high hazard. On the other hand, for values of H close to 0.5 and non-vanishing values of intermittency (C parameter), the hazard is low since random environmental processes can affect the direct hydrothermal vent turning. The validity of the behaviour of H and C as prediction indexes should be tested in future experiments, which will monitor the dynamics of the El-Hierro submarine complex. These values of C and H are likely of importance for the existence of aquatic life in volcano surroundings, but extensive simulations in terms of fractional Langevin equations should be done to see how these parameters, especially C, affect the system. Such an analysis will deliver the memory kernel of the temperature and conductivity, examine possible repeated patterns, and help to study the future behaviour of Tagoro submarine volcano, see [6].

5. Conclusions

The Tagoro shallow submarine volcano, located south of El Hierro Island at the Canary Archipelago, has been deeply monitored in the last ten years since its origin. Its emissions, recorded with high sensitivity temperature and conductivity sensors mounted on an oceanographic rosette, exhibit important fluctuations in the surrounding of the main crater. The 21 θ-Ec time series have been analysed by the GMM method in a short window of observation, delivering a wealth of findings. GMM returns structure functions with convex shapes indicating that both temperature and conductivity fluctuations correspond to multiplicative processes. This behaviour indicates the direct competition between the overlying marine short-term circulation and the degassing material of the volcano leading to long-lived non-equilibrium states. The nonstationary nature of the temperature field points to an open system that releases energy continuously to the environment. The estimated H and C parameters showed that the risk of eruption is low at the corresponding time window. However, continuous monitoring of the volcanic activity is necessary to further compare and classify the future volcano activity. The latter can be done by using the metrics (moments) analysed in the present work, which can be part of a machine learning approach aiming at forecasting the future activity of the Tagoro volcano.

Author Contributions

A.O.A., F.M. and E.F.-N.: data processing, initial mathematical analysis and graphical representations. E.B. and F.Z.: mathematical methodology verification. B.V., A.O.A., E.B., E.F.-N. and F.M.: writing–original draft preparation. B.V., A.O.A., F.M., E.B., F.Z. and E.F.-N.: writing–draft contribution. E.F.-N., F.M. and E.B.: supervision. E.F.-N.: project administration and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FEDER and the Spanish Ministry of Economy and Competitiveness through the VULCANO-II project (CTM2014-51837-R) and for the Spanish Institute of Oceanography through the VULCANA-I (IEO-2015-2017) project, both projects led by the Spanish Institute of Oceanography (IEO), and by Fundación CajaCanarias-LaCaixa through project Cyclovent (2019SP18). Moreover, the Swiss company ORIS Höltsein 1904 funded the scholarship of one of the authors in the Spanish Institute of Oceanography facilities making possible the data preparation, data assessment and training.

Data Availability Statement

Data available under request.

Acknowledgments

We would like to thank the crew of the R/V Ángeles Alvariño and R/V Ramón Margalef for their assistance during the multiple oceanographic cruises.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Edmond, J.M.; Von Damm, K.L.; McDuff, R.E.; Measures, C.I. Chemistry of hot springs on the East Pacific Rise and their effluent dispersal. Nat. Cell Biol. 1982, 297, 187–191. [Google Scholar] [CrossRef]
  2. Baker, E.T.; Massoth, G.J.; Walker, S.L.; Embley, R.W. A method for quantitatively estimating diffuse and discrete hydrothermal discharge. Earth Planet. Sci. Lett. 1993, 118, 235–249. [Google Scholar] [CrossRef]
  3. Speer, K.G.; Rona, P.A. A model of an Atlantic and Pacific hydrothermal plume. J. Geophys. Res. Space Phys. 1989, 94, 6213–6220. [Google Scholar] [CrossRef]
  4. Speer, K.G. A forced baroclinic vortex around a hydrothermal plume. Geophys. Res. Lett. 1989, 16, 461–464. [Google Scholar] [CrossRef]
  5. Fiasconaro, A.; Valenti, D.; Spagnolo, B. Noise in ecosystems: A short review. Math. Biosci. Eng. 2004, 1, 185–211. [Google Scholar] [CrossRef] [PubMed]
  6. Bakalis, E.; Mertzimekis, T.J.; Nomikou, P.; Zerbetto, F. Breathing modes of Kolumbo submarine volcano (Santorini, Greece). Sci. Rep. 2017, 7, 46515. [Google Scholar] [CrossRef] [Green Version]
  7. Bakalis, E.; Mertzimekis, T.J.; Nomikou, P.; Zerbetto, F. Temperature and Conductivity as Indicators of the Morphology and Activity of a Submarine Volcano: Avyssos (Nisyros) in the South Aegean Sea, Greece. Geosciences 2018, 8, 193. [Google Scholar] [CrossRef] [Green Version]
  8. Fraile-Nuez, E.; Santana-Casiano, J.M.; González-Dávila, M.; Vázquez, J.T.; Fernández-Salas, L.M.; Sánchez-Guillamón, O.; Palomino, D.; Presas-Navarro, C. Cyclic Behavior Associated with the Degassing Process at the Shallow Submarine Volcano Tagoro, Canary Islands, Spain. Geosciences 2018, 8, 457. [Google Scholar] [CrossRef] [Green Version]
  9. Carracedo, J.-C. The Canary Islands: An example of structural control on the growth of large oceanic-island volcanoes. J. Volcanol. Geotherm. Res. 1994, 60, 225–241. [Google Scholar] [CrossRef]
  10. Carracedo, J.C.; Day, S.; Guillou, H.; Badiola, E.R.; Canas, J.A.; Perez-Torrado, F.J. Hotspot volcanism close to a passive continental margin: The Canary Islands. Geol. Mag. 1998, 135, 591–604. [Google Scholar] [CrossRef] [Green Version]
  11. Nuez, E.F.; Gonzalez-Davila, M.; Casiano, J.M.S.; Arístegui, J.; Alonso-González, I.J.; Hernández-León, S.; Blanco, M.J.; Rodriguez-Santana, A.; Hernández-Guerra, A.; Gelado-Caballero, M.; et al. The submarine volcano eruption at the island of El Hierro: Physical-chemical perturbation and biological response. Sci. Rep. 2012, 2, 486. [Google Scholar] [CrossRef] [Green Version]
  12. Santana-Casiano, J.M.; Gonzalez-Davila, M.; Nuez, E.F.; De Armas, D.; González, A.G.; Domínguez-Yanes, J.F.; Escánez, J. The natural ocean acidification and fertilization event caused by the submarine eruption of El Hierro. Sci. Rep. 2013, 3, srep01140. [Google Scholar] [CrossRef]
  13. Resing, J.; Baker, E.T.; Lupton, J.E.; Walker, S.L.; Butterfield, D.; Massoth, G.J.; Nakamura, K.-I. Chemistry of hydrothermal plumes above submarine volcanoes of the Mariana Arc. Geochem. Geophys. Geosyst. 2009, 10, 1–23. [Google Scholar] [CrossRef]
  14. Hawkes, J.A.; Connelly, D.P.; Rijkenberg, M.J.A.; Achterberg, E.P. The importance of shallow hydrothermal island arc systems in ocean biogeochemistry. Geophys. Res. Lett. 2014, 41, 942–947. [Google Scholar] [CrossRef] [Green Version]
  15. Rivera, J.; Lastras, G.; Canals, M.; Acosta, J.; Arrese-González, B.; Hermida-Jiménez, N.; Micallef, A.; Tello-Antón, M.; Amblas, D. Construction of an oceanic island: Insights from the El Hierro (Canary Islands) 2011–2012 submarine volcanic eruption. Geology 2013, 41, 355–358. [Google Scholar] [CrossRef]
  16. Ariza, A.; Kaartvedt, S.; Røstad, A.; Garijo, J.C.; Arístegui, J.; Fraile-Nuez, E.; Hernández-León, S. The Submarine Volcano Eruption off El Hierro Island: Effects on the Scattering Migrant Biota and the Evolution of the Pelagic Communities. PLoS ONE 2014, 9, e102354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Coca, J.; Ohde, T.; Redondo, A.; García-Weil, L.; Santana-Casiano, M.; González-Dávila, M.; Arístegui, J.; Nuez, E.F.; Ramos, A.G. Remote sensing of the El Hierro submarine volcanic eruption plume. Int. J. Remote Sens. 2014, 35, 6573–6598. [Google Scholar] [CrossRef] [Green Version]
  18. Eugenio, F.; Martin, J.; Marcello, J.; Fraile-Nuez, E. Environmental monitoring of El Hierro Island submarine volcano, by combining low and high resolution satellite imagery. Int. J. Appl. Earth Obs. Geoinf. 2014, 29, 53–66. [Google Scholar] [CrossRef]
  19. Blanco, M.J.; Fraile-Nuez, E.; Felpeto, A.; Santana-Casiano, J.M.; Abella, R.; Fernández-Salas, L.M.; Almendros, J.; Díaz-del-Río, V.; Domínguez Cerdeña, I.; García-Cañada, L.; et al. Comment on “Evidence from acoustic imaging for submarine volcanic activity in 2012 off the west coast of El Hierro (Canary Islands, Spain)” by Pérez, N.M.; Somoza, L.; Hernández, P.A.; González de Vallejo, L.; León, R.; Sagiya, T.; Biain, A.; González, F.J.; Medialdea, T.; Barrancos, J.; Ibáñez, J.; Sumino, H.; Nogami, K.; Romero, C. (2014). Bull. Volcanol. 2015, 76, 882–896. [Google Scholar] [CrossRef]
  20. Ferrera, I.; Arístegui, J.; Gonzalez, J.; Montero, M.F.; Fraile-Nuez, E.; Gasol, J.M. Transient Changes in Bacterioplankton Communities Induced by the Submarine Volcanic Eruption of El Hierro (Canary Islands). PLoS ONE 2015, 10, e0118136. [Google Scholar] [CrossRef]
  21. Santana-Casiano, J.M.; Nuez, E.F.; Gonzalez-Davila, M.; Baker, E.T.; Resing, J.; Walker, S.L. Significant discharge of CO2 from hydrothermalism associated with the submarine volcano of El Hierro Island. Sci. Rep. 2016, 6, 25686. [Google Scholar] [CrossRef] [PubMed]
  22. Danovaro, R.; Canals, M.; Tangherlini, M.; Dell’Anno, A.; Gambi, C.; Lastras, G.; Amblas, D.; Sanchez-Vidal, A.; Frigola, J.; Calafat, A.; et al. A submarine volcanic eruption leads to a novel microbial habitat. Nat. Ecol. Evol. 2017, 1, 144. [Google Scholar] [CrossRef] [PubMed]
  23. Santana-González, C.; Santana-Casiano, J.M.; González-Dávila, M.; Fraile-Nuez, E. Emissions of Fe(II) and its kinetic of oxidation at Tagoro submarine volcano, El Hierro. Mar. Chem. 2017, 195, 129–137. [Google Scholar] [CrossRef]
  24. Álvarez-Valero, A.; Burgess, R.; Recio, C.; de Matos, V.; Guillamón, O.S.; Gómez-Ballesteros, M.; Recio, G.; Nuez, E.F.; Sumino, H.; Flores, J.; et al. Noble gas signals in corals predict submarine volcanic eruptions. Chem. Geol. 2018, 480, 28–34. [Google Scholar] [CrossRef] [Green Version]
  25. Santana-Casiano, J.M.; González-Dávila, M.; Fraile-Nuez, E. The Emissions of the Tagoro Submarine Volcano (Canary Islands, Atlantic Ocean): Effects on the Physical and Chemical Properties of the Seawater. In Volcanoes: Geological and Geophysical Setting, Theoretical Aspects and Numerical Modeling, Applications to Industry and Their Impact on the Human Health; IntechOpen: London, UK, 2018; Volume 53. [Google Scholar] [CrossRef] [Green Version]
  26. Sotomayor-García, A.; Rueda, J.L.; Sánchez-Guillamón, O.; Urra, J.; Vázquez, J.T.; Palomino, D.; Fernández-Salas, L.M.; López-González, N.; González-Porto, M.; Santana-Casiano, J.M.; et al. First Macro-Colonizers and Survivors Around Tagoro Submarine Volcano, Canary Islands, Spain. Geosciences 2019, 9, 52. [Google Scholar] [CrossRef] [Green Version]
  27. González-Vega, A.; Fraile-Nuez, E.; Santana-Casiano, J.M.; González-Dávila, M.; Escánez-Pérez, J.; Gómez-Ballesteros, M.; Tello, O.; Arrieta, J.M. Significant Release of Dissolved Inorganic Nutrients from the Shallow Submarine Volcano Tagoro (Canary Islands) Based on Seven-Year Monitoring. Front. Mar. Sci. 2020, 6, 829. [Google Scholar] [CrossRef] [Green Version]
  28. Sotomayor-Garcia, A.; Rueda, J.L.; Sánchez-Guillamón, O.; Vázquez, J.T.; Palomino, D.; Fernández-Salas, L.M.; Lopez-Gonzalez, N.; González-Porto, M.; Urra, J.; Santana-Casiano, J.M.; et al. Geomorphic features, main habitats and associated biota on and around the newly formed Tagoro submarine volcano, Canary Islands. In GeoHab Atlas of Seafloor Geomorphic Features and Benthic Habitats; Elsevier: Amsterdam, The Netherlands, 2020; Volume 51, pp. 835–846. [Google Scholar] [CrossRef]
  29. Wang, B.; Anthony, S.; Bae, S.C.; Granick, S. Anomalous yet Brownian. Proc. Natl. Acad. Sci. USA 2009, 106, 15160–15164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Metzler, R.; Klafter, J. The restaurant at the end of the random walk: Recent developments in the description of anomalous transport by fractional dynamics. J. Phys. A Math. Gen. 2004, 37, R161–R208. [Google Scholar] [CrossRef]
  31. Moss, M.E.; Bárdossy, A.; Kotwicki, V.; Strupczewski, W.G.; Mitosek, H.T.; Lian, G.S.; Kemblowski, M.W.; Wen, J.-C.; Lovejoy, S.; Schertzer, D.; et al. New Uncertainty Concepts in Hydrology and Water Resources; Cambridge University Press: Cambridge, UK, 1995; Volume 270, pp. 61–103. [Google Scholar]
  32. Spagnolo, B.; Dubkov, A.A.; Agudov, N.V. Enhancement of stability in randomly switching potential with metastable state. Eur. Phys. J. B 2004, 40, 273–281. [Google Scholar] [CrossRef] [Green Version]
  33. Giuffrida, A.; Valenti, D.; Ziino, G.; Spagnolo, B.; Panebianco, A. A stochastic interspecific competition model to predict the behaviour of Listeria monocytogenes in the fermentation process of a traditional Sicilian salami. Eur. Food Res. Technol. 2009, 228, 767–775. [Google Scholar] [CrossRef] [Green Version]
  34. Biro, T.; Jakovác, A. Power-Law Tails from Multiplicative Noise. Phys. Rev. Lett. 2005, 94, 132302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Caruso, A.; Gargano, M.E.; Valenti, D.; Fiasconaro, A.; Spagnolo, B. Cyclic fluctuations, climatic changes and role of noise in planktonic foraminifera in the mediterranean sea. Fluct. Noise Lett. 2005, 5, L349–L355. [Google Scholar] [CrossRef] [Green Version]
  36. Longpré, M.-A.; Stix, J.; Klügel, A.; Shimizu, N. Mantle to surface degassing of carbon- and sulphur-rich alkaline magma at El Hierro, Canary Islands. Earth Planet. Sci. Lett. 2017, 460, 268–280. [Google Scholar] [CrossRef] [Green Version]
  37. Anteneodo, C.; Riera, R. Additive-multiplicative stochastic models of financial mean-reverting processes. Phys. Rev. E 2005, 72, 026106. [Google Scholar] [CrossRef] [Green Version]
  38. West, B.J. Colloquium: Fractional calculus view of complexity: A tutorial. Rev. Mod. Phys. 2014, 86, 1169–1186. [Google Scholar] [CrossRef]
  39. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  40. Ciuchi, S.; de Pasquale, F.; Spagnolo, B. Nonlinear relaxation in the presence of an absorbing barrier. Phys. Rev. E 1993, 47, 3915–3926. [Google Scholar] [CrossRef]
  41. Ciuchi, S.; de Pasquale, F.; Spagnolo, B. Self-regulation mechanism of an ecosystem in a non-Gaussian fluctuation regime. Phys. Rev. E 1996, 54, 706–716. [Google Scholar] [CrossRef]
  42. Spagnolo, B.; La Barbera, A. Role of the noise on the transient dynamics of an ecosystem of interacting species. Phys. A Stat. Mech. Appl. 2002, 315, 114–124. [Google Scholar] [CrossRef]
  43. Peng, C.-K.; Buldyrev, S.; Havlin, S.; Simons, M.; Stanley, H.E.; Goldberger, A.L. Mosaic organization of DNA nucleotides. Phys. Rev. E 1994, 49, 1685–1689. [Google Scholar] [CrossRef] [Green Version]
  44. Kantelhardt, J.W.; Zschiegner, S.A.; Koscielny-Bunde, E.; Havlin, S.; Bunde, A.; Stanley, H. Multifractal detrended fluctuation analysis of nonstationary time series. Phys. A Stat. Mech. Appl. 2002, 316, 87–114. [Google Scholar] [CrossRef] [Green Version]
  45. Scafetta, N.; Grigolini, P. Scaling detection in time series: Diffusion entropy analysis. Phys. Rev. E 2002, 66, 036130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Hurst, H.E. Long-Term Storage Capacity of Reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 770–799. [Google Scholar] [CrossRef]
  47. Hübscher, C.; Ruhnau, M.; Nomikou, P. Volcano-tectonic evolution of the polygenetic Kolumbo submarine volcano/Santorini (Aegean Sea). J. Volcanol. Geotherm. Res. 2015, 291, 101–111. [Google Scholar] [CrossRef]
Figure 1. (a) Map of El Hierro Island with the triple rift system marked (dashed line) and the location of Tagoro submarine volcano (red triangle). (b) Location of El Hierro Island (red square box) in the Canary Islands Archipelago. (c) High-resolution bathymetric map of the Tagoro submarine volcano with the location of the main crater of 15 m in diameter (red dot), the reference profiles (grey dots) and (d) the distribution of the 21 CTD time series recorded along the crest of the seamount.
Figure 1. (a) Map of El Hierro Island with the triple rift system marked (dashed line) and the location of Tagoro submarine volcano (red triangle). (b) Location of El Hierro Island (red square box) in the Canary Islands Archipelago. (c) High-resolution bathymetric map of the Tagoro submarine volcano with the location of the main crater of 15 m in diameter (red dot), the reference profiles (grey dots) and (d) the distribution of the 21 CTD time series recorded along the crest of the seamount.
Geosciences 11 00374 g001
Figure 2. Mean vertical profile of (a) temperature, θ, and (b) conductivity, Ec, with their standard deviation used as a reference from outlying stations not affected by the submarine volcano emissions. The time series fluctuations measured around and over the main crater, in orange for temperature and dark green for conductivity, are located close to 127 m depth.
Figure 2. Mean vertical profile of (a) temperature, θ, and (b) conductivity, Ec, with their standard deviation used as a reference from outlying stations not affected by the submarine volcano emissions. The time series fluctuations measured around and over the main crater, in orange for temperature and dark green for conductivity, are located close to 127 m depth.
Geosciences 11 00374 g002
Figure 3. (a) Temperature, θ, and (b) conductivity, Ec, time series of CTD30 collected as part of the high-resolution grid. The similarity between both variable fluctuation patterns as well as the occasional distinct peaks is remarkable, and (c) their normalised cross-correlation coefficient rθ,Εc showing a strong correlation for short times (t < 3.5 s), which remains large even for long times. (d) The differentiation with respect to time of the temperature time series, Δθ, (e) the differentiation with respect to time of the conductivity time series, ΔEc, and (f) the normalised cross-correlation coefficient rΔθ,ΔEc of the differentiation with respect to time of temperature and conductivity time series going rapidly to zero, at the crossing point 0.6 s. In grey is indicated the 95% confidence bounds, where the time series are uncorrelated.
Figure 3. (a) Temperature, θ, and (b) conductivity, Ec, time series of CTD30 collected as part of the high-resolution grid. The similarity between both variable fluctuation patterns as well as the occasional distinct peaks is remarkable, and (c) their normalised cross-correlation coefficient rθ,Εc showing a strong correlation for short times (t < 3.5 s), which remains large even for long times. (d) The differentiation with respect to time of the temperature time series, Δθ, (e) the differentiation with respect to time of the conductivity time series, ΔEc, and (f) the normalised cross-correlation coefficient rΔθ,ΔEc of the differentiation with respect to time of temperature and conductivity time series going rapidly to zero, at the crossing point 0.6 s. In grey is indicated the 95% confidence bounds, where the time series are uncorrelated.
Geosciences 11 00374 g003
Figure 4. Logarithmic plots of the variance for (a) temperature, θ, (orange dots) and (b) conductivity, Ec, (dark green dots) with respect to the lag time for the CTD30 time series. The variance for both temperature and conductivity has been fitted by Equation (1). In each figure, the short black solid line represents the best fit for short times (first regime), and the long one the best fit for long times (whole regime). Notice that Δ (s) refers to the unit of time in seconds.
Figure 4. Logarithmic plots of the variance for (a) temperature, θ, (orange dots) and (b) conductivity, Ec, (dark green dots) with respect to the lag time for the CTD30 time series. The variance for both temperature and conductivity has been fitted by Equation (1). In each figure, the short black solid line represents the best fit for short times (first regime), and the long one the best fit for long times (whole regime). Notice that Δ (s) refers to the unit of time in seconds.
Geosciences 11 00374 g004
Figure 5. Generalised moments of CTD30 for (a) temperature, θ, and (b) conductivity, Ec, as function of the lag time, Δ, for different moment orders, q , ranging from 0.25 to 4 with steps of 0.25. The line represents the turning point where the convex shape of the short lag times ends. Notice that Δ (s) refers to the unit of time in seconds.
Figure 5. Generalised moments of CTD30 for (a) temperature, θ, and (b) conductivity, Ec, as function of the lag time, Δ, for different moment orders, q , ranging from 0.25 to 4 with steps of 0.25. The line represents the turning point where the convex shape of the short lag times ends. Notice that Δ (s) refers to the unit of time in seconds.
Geosciences 11 00374 g005
Figure 6. The structure function z ( q ) versus the order of the moment, q , for (a) the first regime (short times) and (b) the whole regime (long times) for temperature, θ, and conductivity, Ec, time series and for the CTD30 data.
Figure 6. The structure function z ( q ) versus the order of the moment, q , for (a) the first regime (short times) and (b) the whole regime (long times) for temperature, θ, and conductivity, Ec, time series and for the CTD30 data.
Geosciences 11 00374 g006
Figure 7. High-resolution grid distribution of (a) temperature, θ, and (b) conductivity, Ec, measurements above of the main crater. The black dots indicate the position where the CTD time series were registered over the seafloor. The filled contours represent the C parameter, while the white contours show the values of the Hurst parameter. The main crater (dashed orange line) is located over the CTDs 30, 34, 35, and 36.
Figure 7. High-resolution grid distribution of (a) temperature, θ, and (b) conductivity, Ec, measurements above of the main crater. The black dots indicate the position where the CTD time series were registered over the seafloor. The filled contours represent the C parameter, while the white contours show the values of the Hurst parameter. The main crater (dashed orange line) is located over the CTDs 30, 34, 35, and 36.
Geosciences 11 00374 g007
Table 1. The location, duration and date of the 21 CTD time series recorded during the multidisciplinary cruise VULCANO-II-1016. The initial and final times correspond to the entire time in which the oceanographic rosette sampler remained in the water. The duration at the seabed (in minutes and seconds) corresponds to the length of the CTD time series fluctuations without the descending and ascending track. All the minutes until the oceanographic rosette reach the volcano edifice and remain steady one meter above the seafloor were discarded, resulting in significatively shorter time series for CTDs 41, 43 and 47.
Table 1. The location, duration and date of the 21 CTD time series recorded during the multidisciplinary cruise VULCANO-II-1016. The initial and final times correspond to the entire time in which the oceanographic rosette sampler remained in the water. The duration at the seabed (in minutes and seconds) corresponds to the length of the CTD time series fluctuations without the descending and ascending track. All the minutes until the oceanographic rosette reach the volcano edifice and remain steady one meter above the seafloor were discarded, resulting in significatively shorter time series for CTDs 41, 43 and 47.
Nº CTDLatitudeLongitudeInitial TimeFinal TimeDuration
CTD3027°37.1788′ N17°59.5831′ W28 October 2016 12:5928 October 2016 13:1711 min. 21 s
CTD3127°37.1793′ N17°59.5893′ W28 October 2016 13:3828 October 2016 13:5612 min. 20 s
CTD3227°37.1793′ N17°59.5944′ W28 October 2016 14:3628 October 2016 14:549 min. 37 s
CTD3327°37.1788′ N17°59.6008′ W28 October 2016 15:1428 October 2016 15:3210 min. 44 s
CTD3427°37.1795′ N17°59.5775′ W28 October 2016 15:5728 October 2016 16:1511 min. 40 s
CTD3527°37.1725′ N17°59.5773′ W28 October 2016 16:3628 October 2016 16:5511 min. 49 s
CTD3627°37.1732′ N17°59.5836′ W28 October 2016 17:1828 October 2016 17:3712 min. 32 s
CTD3727°37.1729′ N17°59.5895′ W28 October 2016 17:5728 October 2016 18:1612 min. 17 s
CTD3827°37.1729′ N17°59.5257′ W28 October 2016 18:2828 October 2016 18:4713 min. 19 s
CTD3927°37.1731′ N17°59.6008′ W28 October 2016 19:4228 October 2016 19:5911 min. 25 s
CTD4027°37.1849′ N17°59.5878′ W28 October 2016 20:1628 October 2016 20:3412 min. 1 s
CTD4127°37.1670′ N17°59.5904′ W28 October 2016 20:5528 October 2016 21:135 min. 0 s
CTD4227°37.1848′ N17°59.5949′ W29 October 2016 07:1029 October 2016 07:2811 min. 35 s
CTD4327°37.1849′ N17°59.5936′ W29 October 2016 07:5629 October 2016 08:136 min. 6 s
CTD4427°37.1674′ N17°59.5948′ W29 October 2016 08:4229 October 2016 08:5911 min. 28 s
CTD4527°37.1680′ N17°59.5792′ W29 October 2016 09:3329 October 2016 09:5010 min. 45 s
CTD4627°37.1704′ N17°59.5822′ W29 October 2016 10:5129 October 2016 11:1111 min. 25 s
CTD4727°37.1701′ N17°59.5895′ W29 October 2016 11:3329 October 2016 11:526 min. 38 s
CTD4827°37.1616′ N17°59.5887′ W29 October 2016 12:2329 October 2016 12:4011 min. 28 s
CTD4927°37.1674′ N17°59.6005′ W29 October 2016 13:0129 October 2016 13:189 min. 53 s
CTD5027°37.1909′ N17°59.5891′ W29 October 2016 13:4829 October 2016 14:0610 min. 46 s
Table 2. γ, α, H, and C parameters of the temperature and conductivity for both first and whole regimes (short and long times) and for each one of the recorded CTD time series. Notice that for short times the alpha index is 2.
Table 2. γ, α, H, and C parameters of the temperature and conductivity for both first and whole regimes (short and long times) and for each one of the recorded CTD time series. Notice that for short times the alpha index is 2.
Temperature θConductivity Ec
First RegimeWhole RegimeFirst RegimeWhole Regime
γαHCγαHCγαHCγαHC
CTD300.8420.6620.0550.751.6100.4240.0590.8520.6690.0540.771.5780.4340.060
CTD310.8020.6200.0410.501.0450.3520.0990.8120.6280.0450.511.0780.3600.100
CTD320.9120.7380.0460.551.1590.3160.0850.9920.7300.0520.611.4280.3300.068
CTD330.7820.5980.0320.851.790.4030.0530.9020.6200.0330.871.7930.3910.052
CTD340.7820.5570.0460.441.0760.3160.0650.9020.5710.0480.461.1040.3240.063
CTD350.7420.5640.0250.421.1360.3820.1240.7520.5780.0280.421.1880.3840.117
CTD360.8020.6010.0280.571.1150.4460.1800.7520.5930.0190.621.2080.4560.168
CTD370.9620.6190.0560.741.6220.4050.0480.8220.6310.0610.791.6330.4220.047
CTD380.9620.6220.0470.551.1900.3900.0670.9620.6440.0550.661.1710.4060.065
CTD390.7320.5680.0210.801.1990.4690.1090.8620.5880.0240.831.2440.4830.110
CTD400.8320.5640.0040.391.3950.2950.1051.0120.5890.0180.421.4800.3060.105
CTD410.8520.6040.0580.261.2750.3070.1020.8420.6270.0580.291.3140.3240.100
CTD420.9320.6140.0210.591.7960.4300.0470.9420.6310.0170.711.8210.4460.055
CTD430.8520.6230.050.621.3670.4260.1260.8420.6420.0560.631.3910.4380.123
CTD440.7720.5980.0640.611.2100.3560.0510.8520.6230.0750.531.2230.3680.055
CTD450.8820.6200.0130.881.9020.4160.0510.9020.6460.0220.901.8470.4300.056
CTD460.7720.5400.0390.522.0000.3450.0170.8220.5730.0480.772.0000.3580.018
CTD470.7320.5590.0320.571.5470.3560.0550.7720.5990.0390.591.5400.3770.061
CTD480.8820.6190.0340.431.4310.3930.1620.9920.6520.0450.211.4450.4050.163
CTD490.7720.5840.0360.591.0000.5360.1300.7920.6000.0430.621.0000.5450.134
CTD500.9720.6620.0370.150.9410.2370.1221.0720.6950.0440.140.9280.2830.131
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Olivé Abelló, A.; Vinha, B.; Machín, F.; Zerbetto, F.; Bakalis, E.; Fraile-Nuez, E. Analysis of Volcanic Thermohaline Fluctuations of Tagoro Submarine Volcano (El Hierro Island, Canary Islands, Spain). Geosciences 2021, 11, 374. https://doi.org/10.3390/geosciences11090374

AMA Style

Olivé Abelló A, Vinha B, Machín F, Zerbetto F, Bakalis E, Fraile-Nuez E. Analysis of Volcanic Thermohaline Fluctuations of Tagoro Submarine Volcano (El Hierro Island, Canary Islands, Spain). Geosciences. 2021; 11(9):374. https://doi.org/10.3390/geosciences11090374

Chicago/Turabian Style

Olivé Abelló, Anna, Beatriz Vinha, Francisco Machín, Francesco Zerbetto, Evangelos Bakalis, and Eugenio Fraile-Nuez. 2021. "Analysis of Volcanic Thermohaline Fluctuations of Tagoro Submarine Volcano (El Hierro Island, Canary Islands, Spain)" Geosciences 11, no. 9: 374. https://doi.org/10.3390/geosciences11090374

APA Style

Olivé Abelló, A., Vinha, B., Machín, F., Zerbetto, F., Bakalis, E., & Fraile-Nuez, E. (2021). Analysis of Volcanic Thermohaline Fluctuations of Tagoro Submarine Volcano (El Hierro Island, Canary Islands, Spain). Geosciences, 11(9), 374. https://doi.org/10.3390/geosciences11090374

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop