Estimation of Sea Level Rise on the West African Coasts: Case of Senegal, Mauritania and Cape Verde

The advent of technology in the world has led to enormous, varying econom-ic and social changes. Those changes have influenced the climate balance, which leads to global temperature rise, hence oceanic expansion and melting glaciers. These various phenomena caused the variation in the mean sea level. Thus, in this study we have evaluated that variation from the instrumental recordings of the tide gauges available in the region, which led us to choose the tide gauges of Dakar (Senegal), Nouakchott (Mauritania) and Palmeira (Cape Verde). The analysis of the tidal series began with the test of presence or absence of significant trends in the series using the Man-Kendall method followed by the decomposition of these series using the SSA, ARIMA and SLT techniques, which would result in the trend, seasonality and residues. The results of this analysis show a variation of 1.52 mm/year, 4.51 mm/year and 20.08 mm/year respectively for the Dakar, Palmeira and Nouakchott tide gauges.


Introduction
The sea level has changed on various time scales. During the last ice age, around 20,000 years ago, the sea level was on average 120 m lower than today. With the melting of the ice caps that then covered northern America and Europe, the world's fastest sea level rise occurred at an average rate of around 10 mm/year between 15,000 and 6000 years before our era [1].
According to geological data, eustatic sea level rise has subsequently occurred at an average rate of 0.5 mm per year for the past 6000 years, and 0.1 to 0.2 mm per year for the past 3000 years. This rate is about 10 times lower than that observed in the 20th century. During the last 1000 years and before the 20 th century, the rate of change in sea level would have been less than 0.2 mm/year [2] [3].
While for 1000 years the sea level has changed little, observations made from core samples and tide gauge measurements show us that at the level of modern man, the rise in sea level has significantly started over the last century. Indeed, for the improvement of living conditions, developed countries used energy resources for their industrialization at the end of the 19 th century [4] [5]. The use of such resources has led to the emission of greenhouse gases which indirectly affect the sea level by warming our planet.
While the mean sea level has risen very little over the past 1000 years, recent analyses using measurements from historical tide gauges indicate a much higher global sea level rise, of around 1.5 to 2 mm/year for the 20 th century [6]. That increase, which is 10 times greater than (those) in previous centuries, is subject to considerable interest in the current context of global warming observed in recent decades. Since the end of 1992, spatial altimetry has also provided absolute (geocentric) measurements of the global mean level of the oceans, on a continuous basis, calculated from the altimetric data. Recent estimates based on the latter measurement technique show a sea level rise of 3.1 ± 0.4 mm/year between November 1992 and January 2018 [7]. The reasons for this rise result from two major processes [8]: the thermal expansion of the ocean caused by variations in the thermal content of the ocean and the inflow of fresh water into the ocean caused by the melting of the polar caps and mountain glaciers, as well as potential exchanges with continental water reservoirs (climatic and anthropogenic effects).
On average over the period 1993-2017, the warming of the oceans in the sea level balance is 42%. The rate of average sea level rise due to thermal expansion was estimated at 1.3 ± 0.4 mm/year between 1993 and 2017. The contribution of continental glaciers and ice caps (Greenland and Antarctica) between 1993 and 2011 has roughly the same contribution in the sea level balance, at 21% and 23%, respectively. The average contribution of continental glaciers and the two polar caps was 0.65 ± 0.15 mm/year between 1993 and 2017. The contribution of ice caps for the same period is 0.73 ± 0.10 mm/year. Inland waters have played little role in long-term sea level change over the past decades [8]. This rise in global sea level will have a lot of consequences for many populations living on the coasts. The physical phenomena linked to the rise in sea level in coastal areas fall into five categories: worsening of the submersion of low-lying coasts (deltas, coral reefs, some islands), acceleration of the erosion of beaches and cliffs, increased salinization of aquifers and surface water, greater amplitude of high water, and finally the effects of natural disasters such as storms and floods [9].
If the overall trend of mean level is increasing, regional causes are mainly due to large-scale changes in ocean density in response to forcing factors (e.g. wind forcing, heat exchange and freshwater at the air-sea interface) and their consequences in ocean circulation. Thus, the greater regional variabilities of the mean level are mainly due to changes in ocean temperature (i.e. from non-uniform thermal expansion), but in some regions variations in ocean salinity are also significant.
The main objective of this project relates to the study of the evolution of the local sea level on the West African coasts (Senegal, Mauritania and Cape Verde) where few scientific studies on this topic have been carried out. As such, specific objectives have been defined: To objectively test whether the three tidal series recorded at the level of the three ports of Dakar (Senegal), Nouakchott (Mauritania) and Palmeira (Cape Verde) show or not significant monotonic trends upward or downward. The statistical test used is the Mann-Kendall nonparametric test which takes into account the presence of seasonality in the series; the case of tidal series.
To carry out a decomposition of the three series of measurements into a set of several independent components, such that each component in this sum can be identified as being the "trend" component, the "seasonality" component and the The rest of this article is organized as follows: The "Tide data used" section presents the tide series considered in this study. The "Methodology" section introduces a pedagogical and exhaustive description of the approaches to analysis of tidal series. The results of detection or not of significant trends within tide gauge series as well as the results of decomposition of these series according to the different techniques (SSA, ARIMA and STL) are presented in the "Experimental results" section. Finally, the conclusions of this study are presented in the last section.

Tide Data Used
Tide gauges are the first instruments to measure variations in sea level. This measurement was taken in the 18 th century by simply reading the water level 124 International Journal of Geosciences on a graduated scale located at certain ports. Mechanical tide gauges were then developed in the 19 th century, allowing automatic recording of the water heights measured by the tide gauge. The first mechanical instruments were float tide gauges. The measurement principle is based on the vertical displacement of a float in a tube, the stilling well, which mechanically filters out low wavelength oscillations from the water surface. The vertical movements of the float are plotted on a roll of paper which then provides an automatic recording of water heights as a function of time [10]. Since then, techniques have been considerably modernized: ultrasonic tide gauges measure the round trip propagation time of an acoustic wave emitted above the body of water, while pressure tide gauges measure the pressure exerted by the column of water above the submerged device. They are digital and remotely searchable instruments for data collection. West Africa appears as the plunging hull in the Atlantic Ocean, despite this peculiarity it has the fewest tide gauges among the world's coasts. The tidal data used in this study (Dakar, Nouakchott and Palmeira) with a monthly resolution are those taken from the database of the global sea level observing system GLOSS (Global Sea Level Observing System).
GLOSS was established by the Intergovernmental Oceanographic Commission (IOC) of the United Nations Educational, Scientific and Cultural Organization (UNESCO) in 1985. GLOSS's main mission is to oversee and coordinate regional and global networks observation of sea level to support scientific research. In the course of its mission, we witness the creation of the CGN in 1990, made up of about 300 tide gauges spread all over the world and most of which are now operational ( Figure 1). In 1997, the second GLOSS implementation plan incorporated tide gauges at a number of selected sites, mainly islands, for the calibration of altimeters (the GLOSS-ALT set) [11].
Dakar station: Dakar (Senegal) has one of the longest sea level records in Africa, with a multi-decadal time series. This station is maintained locally by the Phares et Balises sub-division of the Dakar Port Authority. The year 1902 was associated with the first measurements with the old Dakar tide gauge. The latter operated for a period of less than eighteen years. It was not until 22 years, in 1942, that was experimented a second series of measurements with this tide gauge. This series, seen in 1966, were followed by a long period of absence [12]. Thus a new series has been updated since 1992 and will be subject to analysis in this study. Palmeira resort: Cape Verde is an island state, made up of an archipelago of ten volcanic islands. The volcanic activity of these islands is explained by the presence of a hot spot associated with a bathymetric swell phenomenon, which would have formed these islands. Despite this volcanic activity, it was not until 2000 that a tide gauge was installed in Palmeira. The latter is located in the Atlantic Ocean, off the coast of Senegal.
Nouakchott station: The Autonomous Port of Nouakchott known as "Port de l'Amitié" has had a tide gauge installed at the western end of the quay since May 21, 1999. However, the data available for download is from 2008; this is because the Port de l'Amitié tide gauge is an OTT R 16 type vertical float limigraph. The recording was carried out directly on graph paper. Figure 2 represents the study area and the geographical location of the three tide stations (Dakar, Palmeira and Nouakchott). Figure 3 represents the three tidal series. Table 1 indicates for each tide gauge the geographical position, the range of data used, the percentage of existing data within the series and the rate of change in mm/year estimated by linear regression.

Mann-Kendall Test
The main goal of trend analysis is to find out, objectively, whether there is a significant systematic change in the time series or not. The Mann-Kendall trend test is a non-parametric test widely used to identify a trend in a series, even if there is a seasonal component in the series. This nonparametric trend test is the result of an improvement of the test studied first by [13] then taken up by [14] and finally optimized by [15] so as to take into account a seasonal component.
where ( ) f t is an increasing or decreasing monotonic continuous function.
The residuals i ε , with the same distribution, have a zero mean. The dispersion (or variance) of the distribution is stationary. The data can be annual, monthly, weekly… for a given site. For a series that does not have a seasonal component, the Mann-Kendall statistic is calculated using the following formulation x j and x k are the monthly values of months j and k, with j > k, with the expectation E(S) = 0, and ( ) The variance of S is given by: where t k represents the number of equalities (number of ties) involving k values. The value of S and Var(S) are used to calculate the statistical test Z and is written in the following form If the values of the series tend to increase (respectively fall) over time, Z should be positive (respectively negative).
The probability value p of the statistical test S can be estimated using a cumulative normal distribution function: Thus, if p ≤ 0.05, then the local trend is statistically significant at the 5% risk threshold.
To take into account the seasonality of the data, we define Var(S') as With: ( )  (2).
Using the normal approximation, Z' is defined as: If Z' > 0, the trend is monotonic and increasing. Otherwise, it is decreasing monotonously

Technical Decomposition
Singular Spectral Analysis-SSA Singular Spectral Analysis is a relatively recent technique that consists of breaking down a signal into periodic components, trend and noise. The main idea of SSA is to decompose into the sum of a reduced number of components thanks to the singular value decomposition (SVD) of a specific matrix built from the data ( [16] [17] [18]). The SSA algorithm has two steps [18].
Decomposition: this first step corresponds to performing the SVD of a Hankel matrix called a "trajectory matrix". This is constructed by associating with a signal S formed by N samples (K = NL + 1) column vectors of dimension L, the k th vector being given by ( ) By SVD [8], this matrix of rank r ≤ L breaks down into r elementary matrices X i of rank 1. Let 1 , , r σ σ  be the non-zero singular values of X classified in decreasing order and (U 1 , V 1 ); (U r , V r ) the associated left and right singular vectors Reconstruction: signals 1 , , r X X  of length N are constructed by averaging the elements of the anti-diagonals of the matrices 1 , , r X X  . A grouping into components is then carried out, generally referring to the analysis of the singular values and/or the behavior of the eigenvectors and/or the correlation rate between the X i . This results in m components Y j , each being associated with a group I j This algorithm is thus particularly effective for identifying sinusoidal signals (possibly modulated in amplitude by an exponential) or polynomial functions of time, but these components must be clearly distinguished for this, which leads to the notion of separability within the meaning of the SSA.
ARIMA ARIMA (Autoregressive Integrated Moving Average) [20] is a technique commonly used to fit time series data and to make forecasts. ARIMA models allow three types of temporal processes to be combined: autoregressive processes In the most general case, an ARIMA model combines the three types of random processes, the contribution of each of them being specified by the notation ARIMA (p, d, q), where p is the order of the autoregressive process AR (p), d the degree of integration of a process I (d), and q the order of the moving average MA (q). From a practical point of view, the first parameter to choose is the degree of differentiation d. To do so, we first check the stationarity of the data. If the data is stationary, then d = 0. If the data shows a trend, we take the first difference and check the stationarity again. We keep taking the differences until we get a stationary output. The number of difference required is the number d. We can also calculate the empirical autocorrelations (autocorrelation function ACF) and analyze the rate of decrease towards 0 of this function. If this decrease is slow, we can suspect a non-stationarity. Note that in practice the case d > 2 rarely happens. After determining d, we can use the partial auto correlation function (PACF) to obtain the order p. If the PACF stops after a few shifts, this number is the order of AR. Finally, the q order of MA is determined by examining the autocorrelation (ACF) of the differentiated data. If the ACF stops after a few shifts, this number is the order of MA Seasonal-Trend Decomposition Based on Loess-STL STL initially proposed by [21] and developed by [22] is a non-parametric method which, from an additive model, decomposes a time series (Y t ) into three components: trend (C t ), season (S t ) and residues (ε t ).
STL is an iterative procedure that uses weighted local regression functions, called LOcally wEighted regreSsion Smoother-LOESS. Each point (x 0 , y 0 ) in the dataset is replaced by a linear regression on the points (x i , y i ) of its neighborhood, assigned a weighting depending on the distance |x i − x 0 | with respect to the point (x 0 , y 0 ). This is a nonparametric smoothing determined by the number of points participating in local regressions. The parameters to be defined are the size of the neighborhood and the degree of the polynomial for the adjustment function (constant, linear or quadratic). Due to its ability to manage outliers or missing data in a time series, the Seasonal-Trend Decomposition Based on Loess-STL technique provides a reliable and robust estimate of trend and seasonality. Moreover, according to [23] and [24] its flexibility, simplicity and computing performance offer the possibility of easily processing a large volume of data

Trend Tests
The results of this section relating to the Mann-Kendall trend test are obtained through the "trend" package [25] developed in the R environment. This package is a set of functions dedicated to the analysis of time series. The "trend" functions can be used as a basis for detecting trends, estimating rates of change, de- In the case of the Mann-Kendall test with seasonality, the seasonality of the series is taken into account. In other words, for monthly data with a seasonality of 12 months, we will not seek to know if there is an overall growth over the series, but simply whwther from January to January, from February to the next, so on and so forth, there is a trend. The variance of the statistics can be calculated by assuming that the series are independent (for example the values for the months of January and the months of February are independent).
The results of the application of the Mann-Kendall test with seasonality on the three tide gauge series (Dakar, Palmeira and Nouakchott) are given in Table 2. The p-values calculated from one month to another are lower than the significance level α = 0.05. We therefore reject the null hypothesis H0 (There is no significant trend in the series for the twelve months, and we retain the alternative hypothesis H1 (There is a trend in the series). The three tide series present, from month to month, a significant monotonically increasing trend: S' > 0.

Breakdown of Series
The objective is to decompose the sea level series into a sum of series, such that each component in this sum can be identified as being the trend component, the seasonality component and the noise component.
The results of this section relating to the SSA analysis are obtained through the RSSA package [8] For a good use of RSSA, examples with code fragments under the R environment, dedicated to the explanation of the methodology, are presented in [26].
The width of the "covariance-lag L" window is the only parameter in the decomposition step of the SSA technique. According to [16], the covariance-lag acts as a window on the signal, it must be long enough to encompass the period of the sought signal (exactly as in Fourier spectral analysis…). On this recommendation, we have chosen for our application and for each tide gauge series L close to N/2 where N is the dimension of the series. In our case, L is taken equal to 144 for the Dakar series, 108 for the Palmeira series and 48 for the Nouakchott series.
In order to validate the quality of the SSA reconstruction of each sea level series, two further decompositions of the tidal series were performed using the ARIMA technique and the Seasonal-Trend Decomposition Based on Loess-STL

procedure.
Regarding the decomposition by the ARIMA method, the different steps are implemented in the R programming environment by the functions: "diff" (for the differentiation), "acf" (for the calculation of the function autocorrelation) and "PACF" (partial autocorrelation function). The adjustment of the data according to the possible ARIMA (p, d, q) models is done through the function "arima".
The best ARIMA adjustment model (p, d, q) is selected among the different The tidal measurement is a measurement relating to the earth's crust. The observation of sea level by tide gauges is therefore affected by the vertical movement (in English, Vertical Land Motion-VLM) of the earth. The possible geological origins of this movement are diverse. The physical process affecting tide gauges and the best studied is certainly the postglacial rebound (glacial isostatic adjustment, GIA): it is the only geological process that can be modeled in a global way, even if large uncertainties remain. Successive models of [27] (latest: ICE-6G_C (VM5a), available from the University of Toronto website: site: www.atmosp.physics.utoronto.ca/~peltier/data.php) are most commonly used to correct tide gauge measurements of postglacial rebound. Such models are relatively reliable, but studies show that the estimated velocity values are still too sensitive to some parameters of the earth model, such as the thickness of the li- thosphere or the mantle viscosity, and therefore too uncertain [19]. Table 3 gives the sea level variation rates estimated by applying the linear regression method to the extracted trend components. The columns of this table indicate: the station, the contribution of the trend component in the initial series, the rate of change of the trend component estimated in mm/year and the standard error (SSA, Arima and STL), the correction due vertical movements of the earth's crust (VLM) according to the ICE-6G_C model (VM5a) and the corrected rate of change of VLM (Arima). As can be seen in column 8 of Table 3, all series of tide measurements show increasing rates of change. The highest trend estimated at the Nouakchott station at 20.08 mm/year could be due to a high local subsidence rate.

Conclusions
The study of sea level variability is a complex one that requires powerful statistical tools. The use of these techniques enabled us to analyze the trend of variations in sea level at the level of the Dakar, Nouakchott and Palmeira tide stations. The seasonal Mann-Kendall test is carried out to objectively identify the presence or absence of significant trends in the sea level series. All tide stations show a significant positive trend.
A second analysis was carried out in order to decompose these three series into several components, which can be identified as a non-linear trend, seaso- Due to the phenomenon of climate change, this study shows the effectiveness of statistical tests in detecting trends in tidal series. The results obtained may seem primary, but it is these trends which, for specialists in climatology, are of paramount importance for the search for correlations between climatic events that the Earth undergoes.
In conclusion and in the light of the results obtained, the following recommendations emerge: The need to search for forgotten historical measurements in order to extend the data series as far as possible in the past; The need to conserve the scientific heritage consists of observations of the sea level of the past decades.