Next Article in Journal
Fractional View Analysis System of Korteweg–de Vries Equations Using an Analytical Method
Previous Article in Journal
Fractional-Order LCL Filters: Principle, Frequency Characteristics, and Their Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Impacts of the Global Sea Level Dynamics

1
Faculty of Biology, Shenzhen MSU-BIT University, Shenzhen 518172, China
2
Climate Research Group, Division of Environmental Physics and Meteorology, Faculty of Physics, National and Kapodistrian University of Athens, 15784 Athens, Greece
3
Department of General Ecology and Hydrobiology, Lomonosov Moscow State University, Leninskiye Gory, 1, Moscow 119991, Russia
4
Severtsov Institute of Ecology and Evolution, Russian Academy of Sciences, Leninskiy Ave. 33, Moscow 117071, Russia
5
Section of Condensed Matter Physics, Department of Physics, National and Kapodistrian University of Athens, 15784 Zografos, Greece
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(1), 39; https://doi.org/10.3390/fractalfract8010039
Submission received: 19 November 2023 / Revised: 23 December 2023 / Accepted: 24 December 2023 / Published: 5 January 2024

Abstract

:
The temporal evolution of the global mean sea level (GMSL) is investigated in the present analysis using the monthly mean values obtained from two sources: a reconstructed dataset and a satellite altimeter dataset. To this end, we use two well-known techniques, detrended fluctuation analysis (DFA) and multifractal DFA (MF-DFA), to study the scaling properties of the time series considered. The main result is that power-law long-range correlations and multifractality apply to both data sets of the global mean sea level. In addition, the analysis revealed nearly identical scaling features for both the 134-year and the last 28-year GMSL-time series, possibly suggesting that the long-range correlations stem more from natural causes. This demonstrates that the relationship between climate change and sea-level anomalies needs more extensive research in the future due to the importance of their indirect processes for ecology and conservation.

1. Introduction

Many studies have argued that observational data show that the global climate is warming up and that this phenomenon is accompanied by rising sea levels, which obviously has a significant impact on island and coastal communities. It is believed that the cause of sea level rise is mainly thermal expansion, to which the loss of both glaciers and ice sheets in Greenland and Antarctica contributes significantly. In this regard, Thomas et al. [1] suggested that the combination of ocean thermal expansion and increased loss of ice mass from Greenland is causing the acceleration of sea level rise since the 1970s. To date, there is no model or series of models that perform robust simulations in all modes of climate variability and its impacts (e.g., [2,3,4,5,6,7]).
Although a consensus has been achieved between several studies on the likely increase in GMSL at an accelerated rate, the magnitude of this acceleration needs further investigation because it is strongly dependent on the future of the two large ice sheets. In particular, one of the current “hot topics” is the projections of the GMSL rise by the year 2100. Hu and Bates [8], studying the global average and the regional sea level, found a statistically significant reduction in sea level rise between 2061 and 2080.
They have also found that there are areas where the reduction is insignificant (such as the Philippines and west of Australia) owing to ocean dynamics and the intensification of internal variability due to external forcings.
In connection with this, Vousdoukas et al. [9] predicted a likely increase of the global average extreme sea level in the period 2000–2100 by 34–76 cm. Nerem et al. [10] estimated that the GMSL could rise 65 ± 12 cm by 2100 compared to 2005, which is roughly in line with the IPCC 5th Assessment Report (AR5) model projections.
It is noted that earlier, Nerem et al. [11] discussed the correlation between the ENSO (El Niño Southern Oscillation) event in the period 1997–1998 and the extreme fluctuation of GMSL. They also observed similar sudden variations for each major ENSO event since 1981, indicating ENSO’s contribution to the GMSL variability.
Church and White [12] studied the linear trends of the GMSL using satellite altimeter data (from 1993–2009) and coastal and island sea-level measurements (from 1880–2009). Upward trends were detected in the first and second datasets, with sizes of 3.2 ± 0.4 mm per year and 2.8 ± 0.8 mm per year, respectively. This survey has also shown a remarkable variation in the rate of the sea level rise and significant acceleration since 1880.
Gregory et al. [13] explored the important issue of the observed sea-level increase exceeding the sum of identifiable contributions in size (especially in the twentieth century). It has been found that the largest contributions come from the thermal expansion of the ocean and the melting of glaciers and ice caps. They proposed several reasons for this, such as the failure of climate models to consider the thermal expansion of the oceans caused by volcanic forcing, the rise of glaciers melting compared to previous estimates, and the almost equal and opposite contribution of groundwater depletion and impoundment.
Marzeion et al. [14] evaluated the anthropogenic and natural effect on the global glacier mass loss, suggesting that anthropogenic forcing, in 1851–2010, was responsible for only 25 ± 35% of the global glacier mass loss, while in 1991–2010 it increased to 69 ± 24%.
Jevrejeva et al. [15] studied the natural and anthropogenic effects on the sea-level histories (in the last millennium) extracted from a statistical model. As it was derived, the natural signal was the main cause of sea-level variability observed until 1800. In contrast, the observed sea-level rise of the twentieth century was outside the boundaries of its natural variability, and the observed increase in greenhouse gas concentration appears to account for 75% of reported global sea-level trends.
Similar results were derived from other studies. Slangen et al. [16] investigated the variability of sea level based on a climate model, suggesting that natural effects combined with the response to past climatic change are responsible for 67 ± 23% of the observed sea-level increase before 1950 and only 9 ± 18% after 1970. In contrast, the anthropogenic signal was responsible for only 15 ± 55% of the rise in sea level before 1950 but reached 72 ± 39% in 2000.
Along these lines are the results obtained from the analysis carried out by Marcos et al. [17], where the sea level change (observed in the 20th century) was attributed almost exclusively to anthropogenic causes.
In addition, a few recent studies have attempted to explore scaling dynamics at sea level (i.e., the possibility of long-range dependence, self-similarity, and fractal behavior in the sea-level time series). The term “long-range dependence” describes the property of a quantity so that there are values that remain residually correlated to each other even after many years. The terms “self-similarity” and “fractal” have the meaning that a time-series is exactly or almost like a part of itself. From this point of view, Fraedrich and Blender [18] studied temporal correlations in the surface air temperature over the oceans (i.e., a basic thermodynamic parameter considered to be the main cause of dynamic sea-level behavior). A scaling exponent almost equal to unity was detected, revealing long-term dependence.
Persistent long-range correlations (LRC) were also proposed by Monetti et al. [19], who investigated the scaling properties of the sea surface temperature time series in the Atlantic and Pacific Oceans. The suggested persistence displayed two different types: the short-time type (with a time scale of τ < 10 months), which was characterized by non-stationary behavior for both oceans and the long-term correlation decay. About this, Dangendorf et al. [20] performed a mono-fractal analysis suggesting that sea levels exhibit LRC over time scales up to several decades. They also found long-term correlations to the mass loss from glaciers and ice caps.
Becker et al. [21] studied sea level change over the last century and found that it is beyond its natural internal variability, exhibiting power law long-term correlations on a global and regional scale. Additionally, a power-law scaling behavior was also detected in the long-term sea level variability by Tomasicchio et al. [22]. In this context, Gao et al. [23] recently showed that sea level anomalies exhibit multifractal rather than monofractal behavior.
Multifractal systems are considered a generalization of fractal systems in the sense that their dynamics cannot be described by a single exponent but by a continuous spectrum of exponents.
In conclusion, whereas there’s scientific agreement that global and regional mean sea levels have increased since the late 19th century, the relative contribution of natural and man-made forcing remains unclear.
With this in mind, we are trying in the present study to investigate the scaling dynamics in the temporal evolution of the GMSL using two well-known techniques, detrended fluctuation analysis (DFA) and multifractal DFA (MF-DFA), to contribute to the improvement of the future predictions of the sea-level forecasting models, taking into account any systematic rise [24].

2. Materials and Methods

For the present analysis, we use two sets of GMSL data (in mm). The first contains monthly mean GMSL measurements over the period 1993–2020, obtained from the TOPEX/Poseidon, Jason-1, Jason-2, and Jason-3 satellite altimeter missions http://www.cmar.csiro.au/sealevel/sl_hist_last_decades.html (accessed on 14 September 2023). The second, which includes monthly mean GMSL values from 1880 to 2013, was downloaded from the website http://www.cmar.csiro.au/sealevel/sl_data_cmar.html (accessed on 14 September 2023)., which uses a list of stations with coastal and island sea level measurements. Church et al. [25] presented in detail the careful selection and editing criteria used. In particular, the study of Church and White [12] should be considered along with the Calafat et al. [26] paper, which investigates how well methods based on empirical orthogonal functions can reconstruct global mean sea level [27].
For the study of the LRC and multifractal features of both time series (i.e., the satellite altimeter dataset-SAD and the reconstructed dataset-RD), we employ the DFA and MF-DFA techniques, respectively [28,29,30,31,32,33,34,35,36].
To analyze the multifractal properties of the time series data, the method we used is MF-DFA, whose steps are presented below:
(1)
The first step is to integrate the time series y(i) over time by calculating the differences of the N observations y(i) from their average.
(2)
The next step is to divide the integrated time series, x(i), into completely separate boxes of equal length, τ, and repeat the same algorithm starting this time from the end of the profile, thus obtaining 2Nτ boxes (where Nτ is the integer part of the number N/τ).
(3)
The third step is to calculate the polynomial least-square fit (of order l) in each box and the corresponding variance obtained from the below formulas (see a more detailed description in [30]):
  • for each box j = 1, …, Nτ:
    F 2 ( j , τ ) = 1 τ i = 1 τ x ( ( j 1 ) τ + i ) t ( i ) 2
  • for each box j = Nτ + 1, …, 2Nτ:
    F 2 ( j , τ ) = 1 τ i = 1 τ x ( ( N j N τ ) τ + i ) t ( i ) 2
    where t(i) is a locally best polynomial fitted trend (of second degree) to the τ data.
(4)
In the following, the q-th order fluctuation function is estimated by averaging the variances over all boxes:
F q ( τ ) = 1 2 N τ j = 1 2 N τ F 2 ( j , τ ) q / 2 1 / q
where q is the variable moment.
In case of q→0, the Equation (3) becomes as follows:
F 0 ( τ ) = exp 1 2 N τ j = 1 2 N τ ln F 2 ( j , τ )
(5)
The last step is to depict Fq(τ) vs. τ (in log-log plot) for different values of q and in the case of multi-scaling behavior, a power-law behavior for Fq(τ) is observed:
F q ( τ ) ~ τ h ( q )
where h(q) stands for the generalized Hurst exponent (the slope of the regression line).
Another way to characterize a multifractal series is the singularity spectrum f(n) calculated from h(q) using the modified Legendre transform, where n is the singularity strength or Hölder exponent (see [37]).
It should be remembered that the MF-DFA technique originates from the DFA tool, which evaluates the features of the mono-fractal scaling of a time series [28].
In the case of the original DFA tool, the fluctuation function Fd(τ) is calculated using Equation (3) for q = 2 and first-degree fitted trend (without repeating the algorithm beginning at the end of the profile), i.e.,
F d ( τ ) = 1 N τ j = 1 N τ F 2 ( j , τ ) 1 / 2
A power-law behavior of Fd(τ) (i.e., Fd(τ)~τa) is expected for a fractal series, where a is the monofractal exponent. In the case that the a-exponent belongs to the interval (0, 0.5), power-law anticorrelations (antipersistence) are detected, while an a-exponent value between 0.5 and 1.5 denotes long-range power-law correlations (persistence). For a = 0.5, the series is white noise, while for a = 1, the series is the 1/f noise (1/f noise or pink noise is a process with a frequency spectrum such that the power spectral density is inversely proportional to the frequency of the signal).
It is known that many noisy signals in real systems display trends that deform the scaling results obtained from the DFA method. In this regard, Hu et al. [38] have systematically studied the effects of trends on the DFA results. To avoid such trend interference with our analysis, both GMSL time series (derived from the satellite altimeter dataset and the reconstructed dataset) were initially detrended by applying a polynomial best fit (of sixth-order) and then they were deseasonalized by applying the Wiener filter that produces an estimate of a random process by linear time-invariant filtering [39]. It is worth noting that the sixth-order polynomial fitting gave the most significant results (at a 95% confidence level) without, however, removing the relevant long-term oscillations from the studied time series.

3. Results

3.1. GMSL Derived from Satellite Altimeter Data

Figure 1 illustrates the temporal evolution of the monthly mean GMSL values obtained from SAD during the period 1993–2020. According to Figure 1, a clear upward trend is revealed, with a magnitude of 3.54 ± 0.03 mm per year.
To investigate in-depth, the temporal course of GMSL as derived from SAD, it is important to examine its possible intrinsic scaling features by detecting whether the sea level at different times is correlated. For this reason, we applied the DFA method to the GMSL (from SAD) time series, and the results obtained are presented in detail in the next section.

3.1.1. Application of the DFA Method to the GMSL Derived from SAD

Figure 2a shows the detrended and deseasonalized (D&D) time series of the monthly mean GMSL values as deduced from SAD. The scaling exponent derived from DFA was a = 0.77 ± 0.02, thus indicating the persistent LRC (see Figure 2b).
It should be noted that 0.02 is the standard error of a which is calculated by the equation: s a = i = 1 n [ log F d ( τ i ) y i ] 2 ( n 2 ) i = 1 n [ log τ i log τ ¯ ] 2 and log τ ¯ is the mean value of logτi, i = 1, …, n. It turns out that the t-statistic t = a s a obeys the Student distribution with (n − 2) degrees of freedom. Thus, the 95% confidence interval of a is a s a t n 2 , 0.025   ,   a + s a t n 2 , 0.025 (see Section 7.12.1 of [40]).
Figure 2. (a) GMSL time series (derived from SAD) that have been detrended and deseasonalized over the period 1993–2020. (b) Root-mean-square fluctuation function Fd(τ) (derived from the DFA technique) vs. time scale τ (in months), in a log-log plot, for the time series described in (a) (with the best-fit equation: y = 0.77x − 0.42 and R2 = 0.99). (c) Local slopes a(τ) vs. logτ estimated within a window of 10 and 8 points (circles and crosses). The 2σa(τ) intervals around the mean value a ¯ = 0.77 of the local slopes are shown with the grey line. (d) The fit of the power spectral density vs. frequency f (dotted line) with the power-law line (black) and exponential line (grey).
Figure 2. (a) GMSL time series (derived from SAD) that have been detrended and deseasonalized over the period 1993–2020. (b) Root-mean-square fluctuation function Fd(τ) (derived from the DFA technique) vs. time scale τ (in months), in a log-log plot, for the time series described in (a) (with the best-fit equation: y = 0.77x − 0.42 and R2 = 0.99). (c) Local slopes a(τ) vs. logτ estimated within a window of 10 and 8 points (circles and crosses). The 2σa(τ) intervals around the mean value a ¯ = 0.77 of the local slopes are shown with the grey line. (d) The fit of the power spectral density vs. frequency f (dotted line) with the power-law line (black) and exponential line (grey).
Fractalfract 08 00039 g002
To reject the hypothesis that the above-mentioned scaling dynamics could come from random noise, we used DFA in a 900-time series of random (white) noise (i.e., Monte Carlo simulations) of the same size as the SAD one. By this, we will set the 95% confidence interval of a ¯ (i.e., the average value of a). The confidence interval for the white noise had to be constructed to test (with statistical accuracy) the hypothesis H0: a ¯ = 0.77 vs. a ¯ ≠ 0.77. According to the non-parametric Kolmogorov-Smirnov test, the extracted a-exponents were found to obey a Gaussian distribution at a 95% confidence level with an average a ¯ = 0.53 and a standard deviation σa = 0.09 [41,42]. Therefore, the 95% confidence interval of a ¯ is:
a ¯ 1.96 900 σ a   ,   a ¯ + 1.96 900 σ a = 0.52   ,   0 . 54
Since the DFA exponent of the GMSL time series (from SAD) (i.e., 0.77) was found to be higher than the upper limit of the interval shown in (7), persistent characteristics were disclosed (see Figure 2b). It is worth mentioning that the interval given in Equation (7) does not contain 0.5, indicating that the DFA exponent for the random noise time series is approximately (but not exactly) equal to the value 0.5. Persistent LRC means that sea-level fluctuations, from short to longer intervals (up to 6.3 years or log τ = 1.88), are positively correlated throughout the entire time series. It is noteworthy that the upper limit of 6.3 years (which resulted from the DFA formula and is attributed to the total length of the time series) is very close to the upper limit of irregular intervals (3–7 years) in which the natural tropical phenomenon ENSO occurs. It is an oceanic response to purely stochastic atmospheric forcing having climatological impacts in regions far away from the tropical Pacific (i.e., teleconnections) and may be linked to extreme weather conditions (e.g., floods and droughts), changes in the incidence of epidemic diseases (e.g., malaria), severe coral bleaching, civil conflicts, etc. [31,43,44].
However, to establish the aforementioned LRC in the GMSL time series, we should examine the hypothesis that the power-law decay (i.e., y = 13.53x−0.80 with a coefficient of determination R2 = 0.52) could describe the power spectral density profile more accurately than the exponential decay (i.e., y = 23.1e−0.43x with R2 = 0.40) (Maraun et al. [45]). The power-law decay means that a relative change in frequency results in a relative change in spectral density proportional to the negative power of the change, regardless of the initial magnitude of these quantities. Indeed, the F-test showed that the hypothesis H0: R2power law > R2exponential vs. H1: R2power lawR2exponential could not be rejected at the 95% confidence level. Additionally, another criterion should be used to establish the scaling dynamics in the GMSL time series. This is the stability of “local slopes” in a specified range. For this reason, we applied the DFA method to 1000 series of fractional Gaussian noise with a = 0.8 (i.e., Monte Carlo simulations) to compute local slopes-a(τ) vs. logτ for each of them in two separate window sizes (of 8 and 10 points) that were sequentially shifted to all computed scales τ. Thus, we have determined an interval R = a ¯ 1.96 σ a ( τ )   ,   a ¯ + 1.96 σ a ( τ ) where a ¯ is the mean value and σa(τ) is the standard deviation of the 1000 estimated local slopes-a(τ). Figure 2c shows that in the case of the GMSL time series, the local slopes of logFd(τ) vs. logτ revealed stability, and indeed, the local slopes after logτ = 1 appeared to belong to the R range derived from Monte Carlo simulations.
To have more information about the plausible existence of LRC in the GMSL time series, the next step in our analysis was to investigate the features of the spectrum of singularities for the D&D GMSL time series using the MF-DFA technique.

3.1.2. Application of the MF-DFA Method on the GMSL Derived from SAD

The MF-DFA technique provides an estimate of the fluctuation function Fq(τ) of the q-th order for various moments q. The application of this technique to the D&D GMSL time series (derived from SAD) revealed the expected power-law scaling behavior (i.e., Fd(τ)~τa), which corresponds to straight lines in the log-log plot, on large scales τ > 8 months (τ > 12 months) for all the selected positive (negative) moments q (see Figure 3a) [30,44].
Then, the generalized Hurst exponent h(q) was depicted as a function of q-values to confirm the multifractality of the time series examined (see Figure 3b). The fact that the exponent h(q) varies with q and the h(q) values were higher than 0.5 reveals multifractal behavior and persistent LRC for the GMSL time series. It was also noticed that positive q-values (i.e., large fluctuations) corresponded to lower h(q)-values, which is in line with the usual features of multifractal time series.
Figure 3. (a) The MF-DFA fluctuation factor Fq(τ) as a function of the time scale τ for various moments q for the D&D GMSL time series (SAD) during 1993–2020. The straight lines in the log-log plot, at large scales τ > 8 months (τ > 12 months) for all selected positive (negative) moments q, revealed the expected power-law scaling behavior (i.e., Fd(τ)~τa). (b) The dependence of the generalized Hurst exponent h(q) on q-values for the data shown in (a) (for scales 6 months ≤ τ ≤ 6.3 years). The equation of the best fit is h(q) = 7 × 10−5q3 + 0.001q2  0.03q + 0.83, with R2 = 1.00. (c) The dependence of the singularity spectrum f(n) on the singularity strength n. The equation of the best fit is f(n) = −1.54n3 − 4.27n2 + 10.55n − 3.94, with R2 = 0.99.
Figure 3. (a) The MF-DFA fluctuation factor Fq(τ) as a function of the time scale τ for various moments q for the D&D GMSL time series (SAD) during 1993–2020. The straight lines in the log-log plot, at large scales τ > 8 months (τ > 12 months) for all selected positive (negative) moments q, revealed the expected power-law scaling behavior (i.e., Fd(τ)~τa). (b) The dependence of the generalized Hurst exponent h(q) on q-values for the data shown in (a) (for scales 6 months ≤ τ ≤ 6.3 years). The equation of the best fit is h(q) = 7 × 10−5q3 + 0.001q2  0.03q + 0.83, with R2 = 1.00. (c) The dependence of the singularity spectrum f(n) on the singularity strength n. The equation of the best fit is f(n) = −1.54n3 − 4.27n2 + 10.55n − 3.94, with R2 = 0.99.
Fractalfract 08 00039 g003
The feature of multifractality detected in GMSL might substantially contribute to current scientific knowledge. To provide more information on this feature, we also plotted the singularity spectrum f(n) versus the singularity strength n [30]. The maximum f(n) value corresponds to q = 0 and the f(n) values on both sides of the maximum value correspond to positive or negative moments (see Figure 3c).
It should be clarified and highlighted that the combination of the MF-DFA tool along with the two Maraun criteria [45] is the suitable mathematical methodology to unravel the fractal nature of the evolution of geophysical parameters such as GMSL.

3.2. GMSL Derived from Reconstructed Data

To verify the above results, we repeated the same analysis for the mean monthly GMSL values obtained from RD during a longer period (1880–2013). By studying the temporal evolution of this dataset, an apparent upward trend was detected, with a magnitude of 1.6 ± 0.008 mm per year (see Figure 4a). Figure 4b depicts the mean monthly GMSL values derived from the SAD over the period 1993–2013 compared to the corresponding part of Figure 4a for the same period.
We then attempted to investigate the intrinsic properties of the GMSL (from RD) time series, applying the DFA technique.

3.2.1. Application of the DFA Method on the GMSL Data Derived from RD

The time series of monthly GMSL mean values (from RD) was initially D&D by applying the techniques described in Section 2 (see Figure 5a).
The results of applying the DFA technique to the D&D GMSL time series are shown on a log-log graph in Figure 5b with the corresponding best-fit equation: y = 0.76x + 0.002 and R2 = 0.97. This gives a scaling exponent of a = 0.76 ± 0.02. Moreover, the existence of a power-law scaling and LRC in the GMSL (from RD) time series was established by detecting the type of power spectral density and the stability of the local slopes on long-term scales [45].
Figure 5c depicts the local slopes of logFd(τ) vs. logτ, separately for two separate window sizes of 10 and 12 points, which were shifted successively to all computed scales τ. It is worth noting that these two window sizes were chosen to give the most significant results for the estimated scaling exponents (at the 95% confidence level). As can be seen, the entire local slopes (after logτ = 1.00) are within the boundary of the R range (determined by the 1000 Monte Carlo simulations), indicating sufficient stability. Figure 5d shows the profile of the power spectral density for the D&D GMSL (from RD) time series, suggesting the rejection of exponential decay (y = 165.1e−1.18x with R2 = 0.57) compared to the power law (y = 29.35∙x−1.96 with R2 = 0.65), as the hypothesis H0: R2power law > R2exponential vs. H1: R2power lawR2exponential could not be rejected (using F-test at 95% confidence level).
Therefore, the GMSL (from RD) time series meets both criteria proposed by Maraun et al. [45]. However, to study the multifractality and the power-law long-term persistence of the time series considered, we also applied the MF-DFA technique. The results obtained are given in the following section.

3.2.2. Application of the MF-DFA Method on the GMSL Data Derived from RD

The MF-DFA technique was applied to the D&D GMSL time series (derived from RD) and indicated the expected power-law scaling behavior (i.e., Fd(τ)~τa) which corresponds to straight lines in the log-log plot, on large scales τ > 17 months for all the selected positive and negative moments q (see Figure 6a).
Additionally, the plot of the generalized Hurst exponent h(q) as a function of q (Figure 6b) showed h(q) values clearly higher than 0.5 and dependence of h(q) on q, thus verifying multifractality and the persistent LRC of the time series examined. Also, it is noteworthy that the lower values of h(q) for the positive moments (compared to the negative ones) were consistent with the usual features of the multifractal time series (see Figure 6b).
Finally, we plotted the singularity spectrum f(n) versus the singularity strength n for the specific time series. According to Figure 6c, the maximum f(n) value is around q = 0 and the f(n) values on both sides of the maximum value correspond to positive or negative moments (see Figure 6c).
Figure 6. As in Figure 3, but for the case of the detrended and deseasonalized GMSL time series (derived from RD) during 1880–2013 (for scales 17 months < τ ≤ 33 years) (a). The equation of the best-fitted line in (b) is h(q) = −2 × 10−4·q3 + 5 × 10−4q2 − 0.04q + 0.94, with R2 = 1.00, while in (c) is: f(n) = 2.58n3 − 14.43n2 + 20.06n − 7.18, with R2 = 0.96.
Figure 6. As in Figure 3, but for the case of the detrended and deseasonalized GMSL time series (derived from RD) during 1880–2013 (for scales 17 months < τ ≤ 33 years) (a). The equation of the best-fitted line in (b) is h(q) = −2 × 10−4·q3 + 5 × 10−4q2 − 0.04q + 0.94, with R2 = 1.00, while in (c) is: f(n) = 2.58n3 − 14.43n2 + 20.06n − 7.18, with R2 = 0.96.
Fractalfract 08 00039 g006

3.2.3. Application of the MultiFractal Centered Moving Average (MFCMA) Method

As an additional check on the findings presented above, we employ an additional modern multifractal methodology termed MultiFractal Centered Moving Average (MFCMA). MFCMA was introduced by Schumann and Kantelhardt [46], and the details of its implementation are elaborated in their Section 2.1.
MFCMA, which is less computationally demanding than MFDFA, is based on the Centered Moving Average (CMA) method [47,48]. CMA is slightly better than DFA in the limits of small (τ < 10) and large scales (τ > Ν/4) which, according to Schumann and Kantelhardt [46], makes MFCMA suitable for short time series, which is our present case. Of course, MFCMA performs better for time series in the absence of trends, and for this reason, it was applied to the cases of the detrended and deseasonalized time series shown in Figure 2a and Figure 5a for SAD and RD, respectively.
The results obtained (see, e.g., Figure 7 for RD) lead to singularity spectra comparable to those obtained by MFDFA, with maxima of the singularity spectra at approximately n = 0.83 ± 0.02 and 0.86 ± 0.01, for SAD and RD, respectively. The generalized Hurst exponents h(2) result in 0.80 ± 0.02 and 0.81 ± 0.04, for SAD and RD, respectively. We observe that the results found by DFA and MFDFA are validated by CMA and MFCMA.

3.3. GMSL Derived from Reconstructed and Satellite Data, during the Common Period

The final step of this study was to examine whether the historical data set could be used in any way to validate the SAD values. For this purpose, we focused on the common time period of the two datasets (i.e., January 1993–December 2013), and we compared their trends and their scaling properties.
The GMSL time series from SAD and RD showed statistically significant upward trends of 3.18 ± 0.04 mm per year and 3.56 ± 0.06 mm per year, respectively (see Figure 4b).
Additionally, the application of the DFA technique (during the aforementioned common time period) to the D&D GMSL obtained from SAD and RD gave persistent scaling with exponents a = 0.73 ± 0.03 and a = 0.78 ± 0.03, respectively. These scaling features were established using the criteria of rejection of the exponential decay of the autocorrelation function and the stability of local slopes. Furthermore, the two a-values indicate similar scaling properties, suggesting that fluctuations in sea level, from short time intervals to longer ones (up to 5 years), are positively correlated.
As mentioned in the introduction, the long-term persistence of GMSL has already been evaluated using similar approaches [20,21]. However, as noted by Dangendorf et al. [49], none of the current GMSL reconstructions properly represent the temporal variability on time scales up to decades.
To elaborate more on the interpretation of the results presented above, it should be emphasized that MSL variability is closely related to that of sea surface temperature (SST), which displays persistence [19]. MSL also exhibits long-term correlations [20] with changes that are beyond its natural internal variability [20,21]. However, the relationship between SST and MSL is complex because SST is one of many factors that affect MSL, including winds, currents, river discharges, and gravity fluctuations.
Furthermore, the analysis showed that the long-term correlations are mainly due to natural causes, due to the almost identical scaling characteristics for both the 134-year time series and the latest 28-year GMSL time series. This conclusion is only partially consistent with sea-level change being unnatural in two-thirds of the longest tidal records [21].
It is important to emphasize at this point that although time series with strong persistence may exhibit a large upward trend for natural reasons, we cannot rule out anthropogenic forcing (e.g., in the case of a natural downward trend) [21].
Therefore, as also noted by Baker et al. [21], the investigation of the responsible possible combination of external factors for the observed sea level changes requires further research and, indeed, an increase in the available experimental data.

4. Conclusions

The main findings of the analysis presented above are as follows:
  • Applying the DFA technique to the D&D GMSL time series from the satellite altimeter dataset (reconstructed dataset) during the period 1993–2020 (1880–2013) gives a scaling exponent a = 0.77 ± 0.02 (a = 0.76 ± 0.02), thus revealing that the fluctuations in mean sea-level values from short to longer time intervals are positively correlated.
  • The application of the MF-DFA technique to both GMSL time series used suggested the power-law scaling behavior of Fd(τ) on large scales τ for all the selected positive and negative moments. Additionally, the generalized Hurst exponent h(q) appears to depend on q, and the h(q) values were higher than 0.5, revealing multifractality and persistent long-range correlations.
  • A comparison of the trends and scaling properties of both GMSL time series was carried out for the common period (i.e., January 1993–December 2013). Similar scaling properties were revealed for the two-time series, thus suggesting that the historic data set could be used in any way to validate the satellite altimeter dataset.
The above-mentioned multifractality features detected in the GMSL may contribute to the hot topic of projections of global and regional mean sea-level rise and help to integrate a holistic insight for devising ad hoc strategies and mitigating inevitable impacts [50,51]. In relation to this, under accelerating sea-level rise, it is currently impossible to determine the future evolution of climate, marshes and mangroves as, for example, salt marshes are capable of laterally expanding, contracting, and vertically accumulating in response to sea-level rise [52,53,54].

Author Contributions

Conceptualization, C.V. and Y.M.; methodology, N.V.S., M.E. and C.V.; software, N.V.S. and M.E.; validation, C.V., N.V.S. and M.E.; formal analysis, M.E. and N.V.S.; investigation, C.V., Y.M. and M.E.; resources, Y.M. and D.S.; data curation, M.E.; writing—original draft preparation, C.V. and M.E.; writing—review and editing, Y.M., N.V.S. and D.S.; supervision, C.V. and Y.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thomas, Z.A.; Jones, R.T.; Turney, C.S.; Golledge, N.; Fogwill, C.; Bradshaw, C.J.; Menviel, L.; McKay, N.P.; Bird, M.; Palmer, J.; et al. Tipping elements and amplified polar warming during the Last Interglacial. Quat. Sci. Rev. 2020, 233, 106222. [Google Scholar] [CrossRef]
  2. Chattopadhyay, S.; Jhajharia, D.; Chattopadhyay, G. Univariate modelling of monthly maximum temperature time series over northeast India: Neural network versus Yule–Walker equation based approach. Meteorol. Appl. 2011, 18, 70–82. [Google Scholar] [CrossRef]
  3. Ebel, A.; Memmesheimer, M.; Jakobs, H.J. Chemical perturbations in the planetary boundary layer and their relevance for chemistry transport modelling. In Atmospheric Boundary Layers; Baklanov, A., Grisogono, B., Eds.; Springer: New York, NY, USA, 2007; pp. 109–122. [Google Scholar] [CrossRef]
  4. Yang, A. Long-term global sea-level change due to dynamic topography since 410 MA. Deep Sea Res. Part I Oceanogr. Res. Pap. 2023, 191, 103944. [Google Scholar] [CrossRef]
  5. Varotsos, C.A.; Golitsyn, G.S.; Xue, Y.; Efstathiou, M.; Sarlis, N.; Voronova, T. On the relation between rain, clouds, and cosmic rays. Remote Sens. Lett. 2023, 14, 301–312. [Google Scholar] [CrossRef]
  6. Xue, Y.; Ai, J.; Wan, W.; Guo, H.; Li, Y.; Wang, Y.; Guang, J.; Mei, L.; Xu, H. Grid-enabled high-performance quantitative aerosol retrieval from remotely sensed data. Comput. Geosci. 2011, 37, 202–206. [Google Scholar] [CrossRef]
  7. Reese, A. The Atlantic is being ringed with tiny sensors in an ambitious effort to track climate change. Science 2023, 382, 6670. [Google Scholar] [CrossRef]
  8. Hu, A.; Bates, S.C. Internal climate variability and projected future regional steric and dynamic sea level rise. Nat. Commun. 2018, 9, 1068. [Google Scholar] [CrossRef]
  9. Vousdoukas, M.I.; Mentaschi, L.; Voukouvalas, E.; Verlaan, M.; Jevrejeva, S.; Jackson, L.P.; Feyen, L. Global probabilistic projections of extreme sea levels show intensification of coastal flood hazard. Nat. Commun. 2018, 9, 2360. [Google Scholar] [CrossRef]
  10. Nerem, R.S.; Beckley, B.D.; Fasullo, J.T.; Hamlington, B.D.; Masters, D.; Mitchum, G.T. Climate-change–driven accelerated sea-level rise detected in the altimeter era. Proc. Natl. Acad. Sci. USA 2018, 115, 2022–2025. [Google Scholar] [CrossRef]
  11. Nerem, R.S.; Chambers, D.P.; Leuliette, E.W.; Mitchum, G.T.; Giese, B.S. Variations in global mean sea level associated with the 1997-1998 ENSO event: Implications for measuring long term sea level change. Geophys. Res. Lett. 1999, 26, 3005–3008. [Google Scholar] [CrossRef]
  12. Church, J.A.; White, N.J. Sea-level rise from the late 19th to the early 21st century. Surv. Geophys. 2011, 32, 585–602. [Google Scholar] [CrossRef]
  13. Gregory, J.M.; White, N.J.; Church, J.A.; Bierkens, M.F.; Box, J.E.; Van den Broeke, M.R.; Cogley, J.G.; Fettweis, X.; Hanna, E.; Huybrechts, P.; et al. Twentieth-century global-mean sea level rise: Is the whole greater than the sum of the parts? J. Clim. 2013, 26, 4476–4499. [Google Scholar] [CrossRef]
  14. Marzeion, B.; Cogley, J.G.; Richter, K.; Parkes, D. Attribution of global glacier mass loss toanthropogenic and natural causes. Science 2014, 345, 919–921. [Google Scholar] [CrossRef] [PubMed]
  15. Jevrejeva, S.; Grinsted, A.; Moore, J.C. Anthropogenic forcing dominates sea level rise since 1850. Geophys. Res. Let. 2009, 36, L20706. [Google Scholar] [CrossRef]
  16. Slangen, A.B.; Church, J.A.; Agosta, C.; Fettweis, X.; Marzeion, B.; Richter, K. Anthropogenic forcing dominates global mean sea-level rise since 1970. Nat. Clim. Chang. 2016, 6, 701–705. [Google Scholar] [CrossRef]
  17. Marcos, M.; Marzeion, B.; Dangendorf, S.; Slangen, A.B.A.; Palanisamy, H.; Fenoglio-Marc, L. Internal Variability versus Anthropogenic Forcing on Sea Level and Its Components. In Integrative Study of the Mean Sea Level and Its Components; Springer: Berlin/Heidelberg, Germany, 2017; pp. 337–356. [Google Scholar] [CrossRef]
  18. Fraedrich, K.; Blender, R. Scaling of Atmosphere and Ocean Temperature Correlations in Observations and Climate Models. Phys. Rev. Let. 2003, 90, 108501. [Google Scholar] [CrossRef] [PubMed]
  19. Monetti, A.R.; Havlin, S.; Bunde, A. Long-term persistence in the sea surface temperature fluctuations. Phys. A Stat. Mech. Its Appl. 2003, 320, 581–589. [Google Scholar] [CrossRef]
  20. Dangendorf, S.; Rybski, D.; Mudersbach, C.; Müller, A.; Kaufmann, E.; Zorita, E.; Jensen, J. Evidence for long-term memory in sea level. Geophys. Res. Let. 2014, 41, 5530–5537. [Google Scholar] [CrossRef]
  21. Becker, M.; Karpytchev, M.; Lennartz-Sassinek, S. Long-term sea level trends: Natural or anthropogenic? Geophys. Res. Lett. 2014, 41, 5571–5579. [Google Scholar] [CrossRef]
  22. Tomasicchio, G.R.; Lusito, L.; D’Alessandro, F.; Frega, F.; Francone, A.; De Bartolo, S. A direct scaling analysis for the sea level rise. Stoch. Environ. Res. Risk A 2018, 32, 3397–3408. [Google Scholar] [CrossRef]
  23. Gao, M.; Zhang, A.; Zhang, H.; Pang, Y.; Wang, Y. Multifractality of global sea level heights in the satellite altimeter-era. Phys. A Stat. Mech. Its Appl. 2022, 604, 127923. [Google Scholar] [CrossRef]
  24. Varotsos, C.A.; Cracknell, A.P. Remote Sensing Letters contribution to the success of the Sustainable Development Goals-UN 2030 agenda. Remote Sens. Lett. 2020, 11, 715–719. [Google Scholar] [CrossRef]
  25. Church, J.A.; White, N.J.; Coleman, R.; Lambeck, K.; Mitrovica, J.X. Estimates of the regional distribution of sea level rise over the 1950–2000 period. J. Clim. 2004, 17, 2609–2625. [Google Scholar] [CrossRef]
  26. Calafat, F.M.; Chambers, D.P.; Tsimplis, M.N. On the ability of global sea level reconstructions to determine trends and variability. J. Geophys. Res. Oceans 2014, 119, 1572–1592. [Google Scholar] [CrossRef]
  27. Woodworth, P.L.; Player, R. The permanent service for mean sea level: An update to the 21st century. J. Coast. Res. 2003, 19, 287–295. [Google Scholar]
  28. Peng, C.K.; Buldyrev, S.V.; 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] [PubMed]
  29. Weber, R.O.; Talkner, P. Spectra and correlations of climate data from days to decades. J. Geophys. Res. 2001, 106, 20131–20144. [Google Scholar] [CrossRef]
  30. Kantelhardt, J.W.; Zschiegner, A.; Koscielny-Bunde, E.; Havlin, S.; Bunde, A.; Stanley, H.E. Multifractal detrended fluctuation analysis of nonstationary time series. Phys. A Stat. Mech. Its Appl. 2002, 316, 87–114. [Google Scholar] [CrossRef]
  31. Varotsos, C.A.; Tzanis, C.G.; Sarlis, N.V. On the progress of the 2015–2016 El Niño event. Atmos. Chem. Phys. 2016, 16, 2007–2011. [Google Scholar] [CrossRef]
  32. Efstathiou, M.N.; Varotsos, C.A. Intrinsic properties of Sahel precipitation anomalies and rainfall. Theor. Appl. Climatol. 2012, 109, 627–633. [Google Scholar] [CrossRef]
  33. Efstathiou, M.N.; Tzanis, C.; Cracknell, A.P.; Varotsos, C.A. New features of land and sea surface temperature anomalies. Int. J. Remote Sens. 2011, 32, 3231–3238. [Google Scholar] [CrossRef]
  34. Mintzelas, A.; Sarlis, N.V.; Christopoulos, S.-R.G. Estimation of multifractality based on natural time analysis. Phys. A Stat. Mech. Its Appl. 2018, 512, 153–164. [Google Scholar] [CrossRef]
  35. Nikolopoulos, D.; Moustris, K.; Petraki, E.; Cantzos, D. Long-memory traces in PM 10 time series in Athens, Greece: Investigation through DFA and R/S analysis. Meteorol. Atmos. Phys. 2021, 133, 261–279. [Google Scholar] [CrossRef]
  36. Nikolopoulos, D.; Alam, A.; Petraki, E.; Yannakopoulos, P.; Moustris, K. Multifractal Patterns in 17-Year PM10 Time Series in Athens, Greece. Environments 2022, 10, 9. [Google Scholar] [CrossRef]
  37. Kantelhardt, J.W.; Berkovits, R.; Havlin, S.; Bunde, A. Are the phases in the Anderson model long-range correlated? Phys. A Stat. Mech. Its Appl. 1999, 266, 461–464. [Google Scholar] [CrossRef]
  38. Hu, K.; Ivanov, P.C.; Chen, Z.; Carpena, P.; Stanley, H.E. Effect of trends on detrended fluctuation analysis. Phys. Rev. E 2001, 64, 011114. [Google Scholar] [CrossRef]
  39. Wiener, N. Extrapolation, Interpolation and Smoothing of Stationary Time Series; MIT Technology Press; John Wiley and Sons: New York, NY, USA, 1950. [Google Scholar]
  40. Zwillinger, D. CRC Standard Mathematical Tables and Formulae, 31st ed.; Chapman & Hall/CRC: Boca Raton, FL, USA, 2002. [Google Scholar]
  41. Anderson, T.W.; Darling, D.A. A test of goodness of fit. J. Am. Stat. Assoc. 1954, 49, 765–769. [Google Scholar] [CrossRef]
  42. Stephens, M.A. EDF statistics for goodness of fit and some comparisons. J. Am. Stat. Assoc. 1974, 69, 730–737. [Google Scholar] [CrossRef]
  43. Bowen, K.J.; Ebi, K.L.; Woodward, A.; McIver, L.; Tukuitonga, C.; Nayna Schwerdtle, P. Human health and climate change in the Pacific: A review of current knowledge. Clim. Dev. 2023, 1–15. [Google Scholar] [CrossRef]
  44. Ihlen, E.A.F. Introduction to multifractal detrended fluctuation analysis in Matlab. Front. Physiol. 2012, 3, 141. [Google Scholar] [CrossRef]
  45. Maraun, D.; Rust, H.W.; Timmer, J. Tempting long-memory—On the interpretation of DFA results. Nonlinear Process. Geophys. 2004, 11, 495–503. [Google Scholar] [CrossRef]
  46. Schumann, A.Y.; Kantelhardt, J.W. Multifractal moving average analysis and test of multifractal model with tuned correlations. Phys. A Stat. Mech. Its Appl. 2011, 390, 2637–2654. [Google Scholar] [CrossRef]
  47. Alvarez-Ramirez, J.; Rodriguez, E.; Carlos Echeverría, J. Detrending fluctuation analysis based on moving average filtering. Phys. A Stat. Mech. Its Appl. 2005, 354, 199–219. [Google Scholar] [CrossRef]
  48. Bashan, A.; Bartsch, R.; Kantelhardt, J.W.; Havlin, S. Comparison of detrending methods for fluctuation analysis. Phys. A Stat. Mech. Its Appl. 2008, 387, 5080–5090. [Google Scholar] [CrossRef]
  49. Dangendorf, S.; Marcos, M.; Müller, A.; Zorita, E.; Riva, R.; Berk, K.; Jensen, J. Detecting anthropogenic footprints in sea level rise. Nat. Commun. 2015, 6, 7849. [Google Scholar] [CrossRef] [PubMed]
  50. Frederikse, T.; Landerer, F.; Caron, L.; Adhikari, S.; Parkes, D.; Humphrey, V.W.; Dangendorf, S.; Hogarth, P.; Zanna, L.; Cheng, L.; et al. The causes of sea-level rise since 1900. Nature 2020, 584, 393–397. [Google Scholar] [CrossRef] [PubMed]
  51. Vahsen, M.L.; Blum, M.J.; Megonigal, J.P.; Emrich, S.J.; Holmquist, J.R.; Stiller, B.; Todd-Brown, K.E.O.; McLachlan, J.S. Rapid plant trait evolution can alter coastal wetland resilience to sea level rise. Science 2023, 379, 393–398. [Google Scholar] [CrossRef]
  52. Fagherazzi, S.; Mariotti, G.; Leonardi, N.; Canestrelli, A.; Nardin, W.; Kearney, W.S. Salt marsh dynamics in a period of accelerated sea level rise. J. Geophys. Res. Earth Surf. 2020, 125, e2019JF005200. [Google Scholar] [CrossRef]
  53. Saintilan, N.; Khan, N.S.; Ashe, E.; Kelleway, J.J.; Rogers, K.; Woodroffe, C.D.; Horton, B.P. Thresholds of mangrove survival under rapid sea level rise. Science 2020, 368, 1118–1121. [Google Scholar] [CrossRef]
  54. Steele, R.E.; Reinhardt, E.G.; Devos, F.; Meacham, S.; LeMaillot, C.; Gabriel, J.J.; Rissolo, D.; Vera, C.A.; Peros, M.C.; Kim, S.T.; et al. Evidence of recent sea-level rise and the formation of a classic Maya canal system inferred from Boca Paila cave sediments, Sian Ka’an biosphere, Mexico. Quat. Sci. Rev. 2023, 310, 108117. [Google Scholar] [CrossRef]
Figure 1. Monthly mean values of GMSL from 1993 to 2020 as estimated from SAD.
Figure 1. Monthly mean values of GMSL from 1993 to 2020 as estimated from SAD.
Fractalfract 08 00039 g001
Figure 4. (a) GMSL mean monthly values from 1880 to 2013 as derived from RD, (b) the 1993–2013 part of (a) (red line), and the GMSL mean monthly values derived from SAD for the same period (black line).
Figure 4. (a) GMSL mean monthly values from 1880 to 2013 as derived from RD, (b) the 1993–2013 part of (a) (red line), and the GMSL mean monthly values derived from SAD for the same period (black line).
Fractalfract 08 00039 g004
Figure 5. As in Figure 2, but for the case of the detrended and deseasonalized GMSL time series (derived from RD) during 1880–2013 (a). Panels (c,d) are the same as in Figure 2, i.e., the local slopes a(τ) vs. logτ estimated within a window of 10 and 8 points (circles and crosses) with the grey line depicting the 2σa(τ) intervals around the mean value = 0.76 and the fit of the power spectral density vs. frequency f (dotted line) with the power-law line (black) and exponential line (grey), respectively. The best-fit equation in (b) is: y = 0.76x + 0.002 and R2 = 0.97. The local slopes of logFd(τ) vs. logτ shown in (c) are calculated in the windows of 12 and 10 points (circles and crosses).
Figure 5. As in Figure 2, but for the case of the detrended and deseasonalized GMSL time series (derived from RD) during 1880–2013 (a). Panels (c,d) are the same as in Figure 2, i.e., the local slopes a(τ) vs. logτ estimated within a window of 10 and 8 points (circles and crosses) with the grey line depicting the 2σa(τ) intervals around the mean value = 0.76 and the fit of the power spectral density vs. frequency f (dotted line) with the power-law line (black) and exponential line (grey), respectively. The best-fit equation in (b) is: y = 0.76x + 0.002 and R2 = 0.97. The local slopes of logFd(τ) vs. logτ shown in (c) are calculated in the windows of 12 and 10 points (circles and crosses).
Fractalfract 08 00039 g005
Figure 7. As in Figure 3, after applying MFCMA to the detrended and deseasonalized GMSL time series (derived from RD). The equation of the best-fitted line in (a) is h(q) = 8 × 10−5q3 + 0.000q2 − 0.026q + 0.866 with R2 = 1.00, while in (b) is: f(n) = 2.14n3 − 16.63n2 + 24.12n − 8.780 with R2 = 0.99.
Figure 7. As in Figure 3, after applying MFCMA to the detrended and deseasonalized GMSL time series (derived from RD). The equation of the best-fitted line in (a) is h(q) = 8 × 10−5q3 + 0.000q2 − 0.026q + 0.866 with R2 = 1.00, while in (b) is: f(n) = 2.14n3 − 16.63n2 + 24.12n − 8.780 with R2 = 0.99.
Fractalfract 08 00039 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Varotsos, C.; Mazei, Y.; Sarlis, N.V.; Saldaev, D.; Efstathiou, M. On the Impacts of the Global Sea Level Dynamics. Fractal Fract. 2024, 8, 39. https://doi.org/10.3390/fractalfract8010039

AMA Style

Varotsos C, Mazei Y, Sarlis NV, Saldaev D, Efstathiou M. On the Impacts of the Global Sea Level Dynamics. Fractal and Fractional. 2024; 8(1):39. https://doi.org/10.3390/fractalfract8010039

Chicago/Turabian Style

Varotsos, Costas, Yuri Mazei, Nicholas V. Sarlis, Damir Saldaev, and Maria Efstathiou. 2024. "On the Impacts of the Global Sea Level Dynamics" Fractal and Fractional 8, no. 1: 39. https://doi.org/10.3390/fractalfract8010039

APA Style

Varotsos, C., Mazei, Y., Sarlis, N. V., Saldaev, D., & Efstathiou, M. (2024). On the Impacts of the Global Sea Level Dynamics. Fractal and Fractional, 8(1), 39. https://doi.org/10.3390/fractalfract8010039

Article Metrics

Back to TopTop