Institute of Oceanology, Chinese Academy of Sciences
Article Information
 TAIBI Hebib, HADDAD Mahdi
 Estimating trends of the Mediterranean Sea level changes from tide gauge and satellite altimetry data (19932015)
 Journal of Oceanology and Limnology, 37(4): 11761185
 http://dx.doi.org/10.1007/s0034301981643
Article History
 Received Jun. 12, 2018
 accepted in principle Aug. 12, 2018
 accepted for publication Oct. 12, 2018
The knowledge of sea level is a longstanding concern to the human being. The first motivation, which led man to focus closely on the observation of sea level changes, was probably navigation. This necessity has led man to imagine more and more sophisticated technical means to accurately measure sea level fluctuations; first, by simply installing tide board near the ports; then, by designing instruments of a high technical level, such as pressure sensors or acoustic tide gauges; finally, in recent years, by using artificial altimetry satellites. The tide gauge instruments are employed in coastal areas for the measurement of sea level changes relative to fixed benchmarks on the land, while the satellite altimetry allows measuring in the open ocean surfaces the absolute change of sea level with regard to the surface of the ellipsoid (geocentric reference frame). Therefore, the tide gauge records and space altimetry measurements are complementary for observing the sea level changes. It must be noted, that satellite altimetry is inadequate near the coastlines due to the bad spatial and temporal sampling of altimetry measurements (White et al., 2005) and also to the instability in the altimeter instrument during its passage from the land to the marine environment (FenoglioMarc, 2002).
Tide gauge being fixed to the ground provides a measure of the sea level relative to the coast: it records indeed such Vertical Land Motion (VLMs) of the coast, in addition, to the variations of the sea level. This relative character of the tide gauge compared to the coast is critical for the study of longterm sea level changes because the rates of these changes and the VLMs are of the same order of magnitude (a few millimeters per year). In contrast to tide gauges, altimeters directly measure the absolute sea level, relative to the Earth's center of mass.
Among the main phenomena causing vertical movements of the Earth's crust are: effects of Earth tides, tectonic effects (co and post seismic deformations), postglacial rebound, load effects (vertical crustal deformations under effect of air mass, oceanic and continental waters), subsidence effects related to sediment weight, pumping of groundwater, gas and oil, etc. (Wöppelmann, 1997). The best studied physical process causing VLM of the Earth's crust is certainly the postglacial rebound. This is the only geological process that can be modeled globally, even though large uncertainties remain. For example, Peltier's successive models, most recently: ICE6G (Argus et al., 2014; Peltier et al., 2015) are the most commonly used to correct tide gauge measurements.
At present, sea level change draws much attention as it represents a significant consequence of global climate warming. Global Mean Sea Level (GMSL) is the object of clear increase since the late nineteenth century. More recently, GMSL increased at a rate close to 2 mm/a in the twentieth century (Douglas and Peltier, 2002; Bindoff et al., 2007; Nicholls and Cazenave, 2010; Church and White, 2011) and more recently, since 1993 (satellite altimetry era) by a rate near to 3.2 mm/a (Haddad et al., 2013a; Taibi et al., 2013). The regional variability in the GMSL trend is due to largescale changes in ocean density in response to forcing factors (eg, wind forcing, heat, and freshwater exchange to the airsea interface) and their consequences in the ocean circulation. Thus, the greatest regional variabilities of the sea level are mainly due to changes in ocean temperature (i.e. nonuniform thermal expansion), but in some regions, ocean salinity variations are also important. A recent study based on altimetry data has revealed the major contributions, to the GMSL rise during altimetry era, are from the Indian Ocean with (3.78±0.08 mm/a) (Taibi et al., 2013).
The Mediterranean Sea being a landlocked sea is one of the most vulnerable regions in the world to the impact of climate change, estimating an increase in average surface temperatures in the range of 2.2℃ for the period of 1980–2000 and 5.1℃ for the period of 2080–2100. That further indicates an increasing heat content of this basin that can directly translate into an increase in the mean sea level by the thermal expansion (IPCC, 2013). A possible impact of sealevel rise is enhanced further in the Mediterranean Sea by the presence, particularly along the coasts, of densely populated areas and a highvalue economic activity. There are several studies focusing on sealevel variability and trends in the Mediterranean Sea based on tidegauge and satellite altimetry data analysis (Vignudelli et al., 2005; Pascual et al., 2007; FenoglioMarc et al., 2012; Birol and Delebecque, 2014; Bonaduce et al., 2016). Analyzing satellite altimetry data reported that the Mediterranean Sea level has risen with an average trend of 2.44 mm/a for the period 1993–2013 (Haddad et al., 2013b; Bonaduce et al., 2016).
The aim of this study is to reveal the information stored in the Permanent Service of Mean Sea Level (PSMSL) published datasets of monthly sea level height of 18 tide gauge stations situated along the Mediterranean Sea coast: Tarifa, Malaga Ⅱ, Barcelona, L'Estartit, Sete, Marseille, Toulon, Nice, Marsaxlokk (Formally Valletta), Preveza, Levkas, Katakolon, Kalamai, Piraievs, Khios, Siros, Leros and Ceuta stations. By means of statistical methods, we evaluate the presence or not of significant trends and their magnitude of change by year in these sea level time series. The used statistical methods in this study are the seasonal MannKendall test (Mann, 1945; Kendall, 1975; Hirsch et al., 1982) and the seasonal Sen's method for slope estimates (Sen, 1968; Hirsch et al., 1982). We are motivated to use the MannKendall test because it is a nonparametric test where data must not obey a particular distribution. The second advantage of this test is its low sensitivity to changepoints or abrupt breaks in the time series (Tabari et al., 2011). To estimate the magnitude of trend change per unit time, Sen's method uses a linear model and the variance of the residuals should be constant in time (Salmi et al., 2002). These tests are incorporated in XLSTAT (2017) statistical analysis software that offers a wide variety of functions for data analysis. Prompted by this analysis, the rates of sea level changes from tide gauge measurements, corrected for VLM using the ICE6G model, are compared with those obtained using satellite altimetry data. The main aim of this comparison is to check if the altimetry data are representative of a local sea level trend, as revealed by tide gauge measurements.
The rest of the paper is organized as follows. Section 2 describes the tide gauge and satellite altimetry data used in this study. The seasonal MannKendall test used for trend detection and the Sen's slope approach used to estimate the magnitude of trends are explained briefly in Section 3. The obtained results are given in Section 4. These results are followed by a comparison between trends from tide gauge records and satellite altimetry measurements. Section 5 concludes.
2 DATA 2.1 Tide gauge dataThe study of the trend of sea level variability does not necessarily require direct manipulation of hourly tide gauge data. Indeed, the mechanism for calculating of average hourly data for one month or one year filters the fluctuations of short period observed in tide gauge records, that they are of irregular nature (waves of the storm, tidal wave…) or periodical (diurnal, tidal waves, tides…). The time series of monthly or annual averages are consequently suitable under investigation of the sea level variability. In this context, we have chosen within the framework of this study to analyze the available monthly averages sea level series. In the Mediterranean basin, there are nearly 100 tide gauge stations with Revised Local Reference (RLR) datasets covering the period from 1993 to 2015 (satellite altimetry era) (Holgate et al., 2013; PSMSL, 2018). Less than half of these have complete time series. Accordingly, we selected 18 tide gauge stations that have at least 80% monthly data present during the analysis period: Tarifa, Malaga Ⅱ, Barcelona, L'Estartit, Sete, Marseille, Toulon, Nice, Marsaxlokk (Formally Valletta), Preveza, Levkas, Katakolon, Kalamai, Piraievs, Khios, Siros, Leros and Ceuta stations. Figure 1 shows the spatial distribution of these tide gauge stations over the Mediterranean Sea.
The corresponding monthly tidal heights series are obtained from the PSMSL. These data are of the "RLR" category, i.e., those checked in terms of quality, continuity, and local stability of the tide gauge reference. A description of data checks routinely made by the PSMSL is given in (Woodworth et al., 1990) and (IOC, 1992). Examination of sea level time series from these 18 stations is very important to perform and interpret the trend analysis. Table 1 presents the descriptive statistics of the considered series for the entire period of study (1993–2015) as well as their completeness expressed in a number of observations without missing data. The magnitude of change and its standard deviation for each tide gauge are estimated using the classical least squares method that finds the line of best fit for a dataset.
The RLR datum at each station is approximately 7 000 mm below mean sea level. This arbitrary choice was made by the PSMSL in order to avoid negative values in the resulting RLR monthly sea level time series (see Table 1, Columns from 7 to 9). The standard deviation in Column 10 of Table 1, that quantifies the dispersion of data, shows that the sea level at each station is object of strong variations due probably to the annual signal present in the Mediterranean Sea level variations (standard deviation from 56 to 95 mm) (Haddad et al., 2013a). As can be seen in Column 11 of Table 1, all the sea level series show a positive trend. However, concerning the Marsaxlokk (Formally Valletta) station, the reliability of its slope is rather suspicious in view of its error margin that is of the same order of magnitude (1.56±1.03 mm/a). This may be due to the multiple dysfunctions and repositioning of the float gauge thus inducing bias on the zero datum of the Marsaxlokk station. The history of this station is available from the Permanent Service for Mean Sea Level (PSMSL) web site (http://www.psmsl.org/data/obtaining/stations/1735.php). Here, we must take into account the problem of the presence or not of significant monotonic trends in time series before any estimation of the trend.
2.2 Satellite altimetry dataIn this study, we also used the monthly gridded Maps of Sea Level Anomaly (MSLA) data that cover the Mediterranean basin. These data are produced by AVISO (Archiving, Validation, and Interpretation of Satellite Oceanographic data) and are available through the website: http://www.aviso.oceanobs.com/. These maps cover the period from January 1993 to December 2015, with a spatial resolution of 1/8×1/8 degrees in latitude and longitude.
Sea Level Anomaly (SLA), that represents the variation of the sea surface height with respect to a severalyear mean or a mean sea surface, is generally used as a precious and main indicator for the analysis of ocean variability. The main input data for processing MSLA grids are the Geophysical Data Records produced from different space altimetry missions that are therefore of the highest quality, notably in terms of satellite's orbit determination (position and altitude). All of the standard corrections to the altimeter data were applied including tropospheric and ionospheric corrections, removal of ocean tides, sea state bias correction, and inverted barometer correction (CLS, 2018). These maps of SLA are available in NetCDF format. This format is a generic format and multiplatform data storage. An exhaustive presentation of the NetCDF conventions is available on the Unidata website: http://www.unidata.ucar.edu/. The necessary routines for using the NetCDF are available on this website.
To highlight the spatial variations of the Mediterranean Sea level, a map (30°N–46°N, 6°W–37°E), showing the sea level variation rates (drift or velocity in mm/a), was performed from the monthly maps of SLA (1993–2015) by estimating the linear slopes for each cell of this map (Fig. 2). As shown in Fig. 2, the Mediterranean Sea level is characterized by an EastWest differentiation. While the Eastern Mediterranean basin is the object of a clear increasing trend, an inverse situation occurs at the Ionian Sea where the sea level is falling. The observed trends are in agreement with the results of Haddad et al. (2013b) obtained by performing a singular spectrum analysis (SSA) with satellite altimetry data to estimate the seasonal cycle and trend, in almost the same analysis period.
3 METHOD 3.1 MannKendall testThe MannKendall test (Mann, 1945; Kendall, 1975) is a nonparametric test that is used commonly to detect the presence or not of significant trends in environmental, climate, or hydrological time series. The null hypothesis, H_{0}, is that the data come from a population with independent realizations that does not follow any trend. The alternative hypothesis, H_{1}, is that the data follow a clear monotonic trend.
Considering a time series of n data points and t_{i} and t_{j} as two subsets of data where i=1, 2, 3, …, n–1 and j=i+1, i+2, i+3, …, n, the MannKendall statistic S is calculated by analyzing all possible pairs of measurements in the dataset. If the earlier measurement is less in magnitude than a later one, S is incremented by one. On the other hand, if an earlier value is greater in magnitude than a later sample, S is decremented by one. For two equivalent measurement values, S remains unchanged:
with:
The variance for the Sstatistic is defined by:
where n: size of the series, g: represents the number of groups of ties in the time series, and t_{j}: the number of ties in the j^{th} group of ties.
The statistic S is approximately normal distributed. The standardized Zstatistic for an increasing (or decreasing) trend is computed using the equation:
The statistic Z is used as an indicator of the significance of the trend. If Z is greater than Z_{α/2}, where α represents the chosen significance level (eg: 5% with Z_{0.025}=1.96), then the null hypothesis H_{0} is invalid implying that the trend is significant (Motiee and McBean, 2009). A positive value of Z indicates an increasing trend and a negative value indicates a decreasing trend. Using Pvalue calculated for Z, H_{0} is rejected if P < α.
The seasonal MannKendall (Hirsch et al., 1982) is a simple modification to the MannKendall test that accounts for apparent seasonal fluctuations in the series. It is suitable for seasonal data with a moderate level of autocorrelation between the different components of the series and may be used even in the case of presence of missing data in the series. This seasonal approach runs a separate Mann Kendall trend test on each of m seasons separately, where m is the number of separate seasons (Spring, Summer etc.). Seasonal can also refer to other time periods, like hours, days, or months.
The test consists of computing the MannKendall statistic S_{i} and its variance Var(S_{i}) for each season with data collected over years. For monthly time series, January would only be compared with January and February would only be compared with February, etc. The overall seasonal MannKendall statistic S' is calculated from summing each season's Kendall S statistic, and the Z' statistic is computed:
The hypothesis H_{0} is rejected in favor of the alternative hypothesis of a downward trend if Z' is negative and If Z' is greater than Z'_{α/2}. Using Pvalue calculated for Z', H_{0} is rejected if P < α. For further information, the reader is referred to Hipel and McLeod (1994) and Hirsch et al. (1982).
3.2 Sen's slope estimatorSen's slope estimation (Sen, 1968) is a nonparametric method that is used for analysis of the magnitude of the trend. The advantage of using Sen's slope is that it is not affected by the possible presence of outliers and data errors in the dataset (Salmi et al., 2002). This method assumes that the trend line is a linear function in the time series:
where Q is the slope and B is a constant.
The slopes of all data pairs t_{j} and t_{k} for j and k times of a period, where j>k, are calculated as follow:
Median is computed from n observations of the slope to estimate the Sen's Slope estimator.
The simple Sen's Slope estimator can be modified to account for seasonality in the observations (Hirsch et al., 1982). It is calculated by the following formula:
For each (t_{ij}, t_{ik} pair i=1, 2, …, m, where 1≤k≤j≤n_{i} and n_{i} is the number of known values in the i^{th} season. The seasonal slope estimator is the median of the Q_{ijk} values.
Software used in this study for performing the statistical seasonal MannKendall test and seasonal Sen's slope estimator is Addinsoft's XLSTAT (2017).
4 RESULT AND DISCUSSIONThe trend analysis of the sea level series has been done using the data of Tarifa, Malaga Ⅱ, Barcelona, L'Estartit, Sete, Marseille, Toulon, Nice, Marsaxlokk (Formally Valletta), Preveza, Levkas, Katakolon, Kalamai, Piraievs, Khios, Siros, Leros and Ceuta tide gauge stations. Here, we used a Pvalue approach to test the null hypothesis, H_{0}, of no trend versus the alternative hypothesis, H_{1}, of either an upward or downward trend. Table 2 gives the results of the running of the seasonal MannKendall test on tidal heights data for the considered 18 stations. The hypothesis H_{0} is rejected when Pvalue is less than the considered significance level α of 0.05. The reject of the null hypothesis meaning the presence of a trend in the series is statistically significant. Table 2 indicates that acceptance of the null hypothesis (H_{0}) occurred only for the Marsaxlokk (Formally Valletta) station (the Pvalue of the test of 0.374 is greater than the significance level α of 0.05). For the remainder, the results show a significant positive trend (Pvalue < α). We can conclude that relative sea level changes along the Mediterranean Sea coast are the object of increasing behavior during 1993–2015.
Using the seasonal Sen's approach, the rate of change over time (in mm per year) of each sea level time series, that exhibits a statistically significant monotonic trend, is computed as the median of all slopes between data pairs within the same season. The obtained trend values are given in Table 3, Column 5. Corrections for VLM due to postglacial rebound (Table 3, Column 6) are derived from predictions of the ICE6G model of Peltier et al. (2015), which is available from the University of Toronto web site:www.atmosp.physics.utoronto.ca/~peltier/data.php.
An increasing annual trend is detected for all stations that exhibit a statistically significant trend. After VLM corrections, the largest annual trend is noticed for the Levkas station with a rate of change of 8.72 mm/a, whereas the low annual trend is observed for the Leros station with a value of 1.48 mm/a (Table 3, Column 7). These results provide further evidence, if needed, the rise in the costal Mediterranean Sea level.
With the aim to assess whether offshore altimetry signals were representative in term of temporal variability of the local sea level along the coasts of the Mediterranean basin, as seen from the tide gauges records, the magnitude of change at the location of each tide gauge station was estimated from the drift map (Fig. 2) using the nearestneighbor interpolation method.
The obtained rates of change, estimated from the analysis of insitu tidal gauge and satellite altimetry data (Table 3, Columns 7 and 8), show positive linear trends during 1993–2015. These results are sufficient to admit the rise in the coastal sea level of the Mediterranean Sea. The differences between the estimated rates from insitu and satellite measurements vary, in absolute value, from 0.49 to 5.57 mm/a in an average of 1.94 mm/a (Table 3, Column 9). These differences corroborate the agreement and the complementarity of the two techniques of measurement (satellite altimetry and tide gauge) to cover the Earth's ocean surfaces, just as well in the open ocean as in coastal areas.
Nevertheless, note that the main reasons for the difference between the trends estimated from these two techniques are due to the lack of altimetry data in the coastal areas of the fact that the observations are confined typically every 10 days along the ground tracks, which decreased the accuracy of the corrective terms within this dynamic environment operating between land and sea. Tide gauges are fundamentally different which measuring in situ sea level each second to hour accurately at selected coastal locations.
For illustrative purposes, the rates of change estimated from tide gauges and satellite altimetry measurements for the period from 1993 to 2015 are illustrated in Fig. 3. The highest tide gauge sealevel trend is estimated at the Levkas station to 8.72 mm/a, on the Western coast of Greece. This value is 5.57 mm/a higher than the absolute trend estimated at the nearest satellite grid point, 3.15 mm/a (see Fig. 2), suggesting a strong subsidence rate at this site.
5 CONCLUSIONThe main objective of this study is to take advantage of tide gauge and satellite altimetry techniques to monitor sea level variability along the Mediterranean Sea coast. To that end, we have analyzed 18 tidal height series that are available over the satellite altimetry era period (1993–2015). The seasonal MannKendall test is performed to identify statistically significant trends of these series. The magnitudes of sea level change, estimated by the use of seasonal Sen's approach, show that there is an increase in sea level for all tide gauge stations that exhibit a significant trend. After correction for vertical land motions (VLM) due to postglacial rebound, the magnitudes of changes are from 1.48 to 8.72 mm/a during 1993–2015.
A second analysis was carried out in order to assess whether the altimetry data can be representative of the local sea level along the Mediterranean coast. For this purposes, monthly maps of sea level anomaly derived from different space altimetry missions at the Mediterranean scale were analyzed. The results indicate that there is a clear differentiation between the behavior of the Eastern and the Western parts of the Mediterranean Sea. While a greater rise was observed in the Eastern Mediterranean Sea, the Ionian Sea level fell during 1993–2015.
These rates from the altimetry measurements were exploited to provide the magnitude of change at the location of each tide gauge station, using the nearestneighbor interpolation method. The differences between the estimated rates of sea level change from insitu tide gauge and satellite measurements vary, in absolute value, from 0.49 to 5.57 mm/a with an average of 1.94 mm/a. One of the reasons for these differences can be attributed to the fact that the altimetry measurements are not accurate near the coasts. Nevertheless, all local trends, estimated from the set of insitu tide gauge and satellite altimetry data, show positive rates during 1993–2015. This finding is sufficient to admit firstly the coastal sea level rise of the Mediterranean basin, and secondly, the complementarity of the two measurement techniques (satellite altimetry and tide gauges).
In the face of climate change phenomenon, this experimentation shows the efficacy of statistical tests in detecting the sea level trends derived from the tidal gauge measurements. The obtained results may appear primary; however, these trends are essential for climatologists to determine what correlations exist between climatic events that the Earth is undergoing. To conclude, evidence from this study confirms that the coastal Mediterranean Sea level has risen during the period of the altimetry era (1993–2015).
6 DATA AVAILABILITY STATEMENTThe tide gauge data that support this study are available from the Permanent Service of Mean Sea Level (PSMSL) Website: http://www.psmsl.org. The altimetry data that support this study are available from AVISO (Archiving, Validation, and Interpretation of Satellite Oceanographic data) through the website: http://www.aviso.oceanobs.com/.
7 ACKNOWLEDGEMENTThe authors are enormously grateful to PSMSL for providing sea level time series and to AVISO for providing altimetry data. The authors greatly thank the anonymous reviewers for their valuable and constructive comments.
Argus D F, Peltier W R, Drummond R, Moore A W. 2014. The Antarctica component of postglacial rebound model ICE6G_C (VM5a) based on GPS positioning, exposure age dating of ice thicknesses, and relative sea level histories. Geophysical Journal International, 198(1): 537563.
DOI:10.1093/gji/ggu140 
Bindoff N L, Willebrand J, Artale V, Cazenave A, Gregory J, Gulev S, Hanawa K, Le Quéré C, Levitus S, Nojiri Y, Shum C K, Talley L D, Unnikrishnan A. 2007.Observations: oceanic climate change and sea level. In: Climate Change 2007: The Physical Science Basis.Cambridge University Press, Cambridge.

Birol F, Delebecque C. 2014. Using high sampling rate (10/20Hz) altimeter data for the observation of coastal surface currents:a case study over the northwestern Mediterranean Sea. Journal of Marine Systems, 129: 318333.
DOI:10.1016/j.jmarsys.2013.07.009 
Bonaduce A, Pinardi N, Oddo P, Spada G, Larnicol G. 2016. Sealevel variability in the Mediterranean Sea from altimetry and tide gauges. Climate Dynamics, 47(910): 28512866.
DOI:10.1007/s0038201630012 
Church J A, White N J. 2011. Sealevel rise from the late 19th to the early 21st century. Surveys in Geophysics, 32(45): 585602.
DOI:10.1007/s1071201191191 
CLS. 2018. SSALTO/DUACS Experimental Product Handbook, SALPMUPEA23172CLS, 49 pp. https://www.aviso.altimetry.fr/fileadmin/documents/data/tools/hdbk_duacs_experimental.pdf. Accessed on 24 Sep. 2018.

Douglas B C, Peltier W R. 2002. The puzzle of global sealevel rise. Physics Today, 55(3): 3540.
DOI:10.1063/1.1472392 
FenoglioMarc L, Braitenberg C, Tunini L. 2012. Sea level variability and trends in the Adriatic Sea in 19932008 from tide gauges and satellite altimetry. Physics and Chemistry of the Earth, Part A/B/C, 4041: 4758.
DOI:10.1016/j.pce.2011.05.014 
FenoglioMarc L. 2002. Longterm sea level change in the Mediterranean Sea from multisatellite altimetry and tide gauges. Physics and Chemistry of the Earth, Part A/B/C, 27(3234): 14191431.
DOI:10.1016/S14747065(02)000840 
Haddad M, Hassani H, Taibi H. 2013b. Sea level in the Mediterranean Sea:seasonal adjustment and trend extraction within the framework of SSA. Earth Science Informatics, 6(2): 99111.
DOI:10.1007/s1214501301146 
Haddad M, Taibi H, Mohammed Arezki S M. 2013a. On the recent global mean sea level changes:trend extraction and El Niño's impact. Comptes Rendus Geoscience, 345(4): 167175.
DOI:10.1016/j.crte.2013.03.002 
Hipel K W, McLeod A I. 1994. Time Series Modelling of Water Resources and Environmental Systems. Elsevier, Amsterdam.

Hirsch R M, Slack J R, Smith R A. 1982. Techniques of trend analysis for monthly water quality data. Water Resources Research, 18(1): 107121.
DOI:10.1029/WR018i001p00107 
Holgate S J, Matthews A, Woodworth P L, Rickards L J, Tamisiea M E, Bradshaw E, Foden P R, Gordon K M, Jevrejeva S, Pugh J. 2013. New data systems and products at the permanent service for mean sea level. Journal of Coastal Research, 29(3): 493504.

IOC. 1992. Joint IAPSOIOC Workshop on Sea Level Measurements and Quality Control. Intergovernmental Oceanographic Commission, Workshop Report No. 81, IOC, Paris.

IPCC. 2013. The Physical Science Basis. http://www.climatechange2013.org. Accessed on 17 July 2018.

Kendall M G. 1975. Rank Correlation Methods. 4^{th} edn.Charles Griffin, London.

Mann H B. 1945. Nonparametric tests against trend. Econometrica, 13(3): 245259.
DOI:10.2307/1907187 
Motiee H, McBean E. 2009. An assessment of long term trends in hydrologic components and implications for water levels in lake superior. Hydrology Research, 40(6): 564579.
DOI:10.2166/nh.2009.061 
Nicholls R J, Cazenave A. 2010. Sealevel changes and their impacts on coastal zones. Science, 328(5985): 15171520.
DOI:10.1126/science.1185782 
Pascual A, Pujol M I, Larnicol G, Traon P Y L, Rio M H. 2007. Mesoscale mapping capabilities of multisatellite altimeter missions:first results with real data in the Mediterranean Sea. Journal of Marine Systems, 65(14): 190211.
DOI:10.1016/j.jmarsys.2004.12.004 
Peltier W R, Argus D F, Drummond R. 2015. Space geodesy constrains ice age terminal deglaciation:the global ICE6G_C (VM5a) model. Journal of Geophysical Research:Solid Earth, 120(1): 450487.
DOI:10.1002/2014JB011176 
Permanent Service for Mean Sea Level  PSMSL. 2018. http://www.psmsl.org, Accessed on 20180218.

Salmi T, Määttä A, Antilla P, RuohoAirola T, Amnell T. 2002. Detecting Trends of Annual Values of Atmospheric Pollutants by the MannKendall Test and Sen's Slope Estimates  the Excel Template Application Makesens.Finnish Meteorological Institute, Helsinki, Finland, 35p.

Sen P K. 1968. Estimates of the regression coefficient based on Kendall'S tau. Journal of the American Statistical Association, 63(324): 13791389.
DOI:10.1080/01621459.1968.10480934 
Tabari H, Marofi S, Ahmadi M. 2011. Longterm variations of water quality parameters in the Maroon River, Iran. Environmental Monitoring and Assessment, 177(14): 273287.
DOI:10.1007/s106610101633y 
Taibi H, Kahlouche S, Haddad M, Rami A. 2013. Trends in global and regional sea level from satellite altimetry within the framework of AutoSSA. Arabian Journal of Geosciences, 6(12): 45754584.
DOI:10.1007/s1251701207762 
Vignudelli S, Cipollini P, Roblou L, Lyard F, Gasparini G P, Manzella G, Astraldi M. 2005. Improved satellite altimetry in coastal systems:case study of the Corsica channel (Mediterranean Sea). Geophysical Research Letters, 32(7): L07608.

White N J, Church J A, Gregory J M. 2005. Coastal and global averaged sea level rise for 1950 to 2000. Geophysical Research Letters, 32(1): L01601.

Woodworth P L, Spencer N E, Alcock G A. 1990. On the availability of European mean sea level data. International Hydrographic Review, 67(1): 131146.

Wöppelmann G. 1997. Rattachement géodésique des marégraphes dans un système de référence mondial par techniques de géodésie spatiale. Paris Observatory, France.

XLSTAT. 2017. Data Analysis and Statistical Solution for Microsoft Excel. Addinsoft, Paris, France.
