ORIGINAL RESEARCH article

Front. Mar. Sci., 31 October 2025

Sec. Ocean Observation

Volume 12 - 2025 | https://doi.org/10.3389/fmars.2025.1683307

Estimating the Mediterranean sea-level trends from 1993 to 2022 showing the consistency of X-TRACK-L2P coastal altimetry with tide gauge data and trend impacts on Egypt's northern coast

  • 1. State Key Laboratory of Precision Geodesy, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan, China

  • 2. University of Chinese Academy of Sciences, Beijing, China

  • 3. Civil Engineering Department, Benha Faculty of Engineering, Benha University, Benha, Egypt

  • 4. Civil Engineering Department, Faculty of Engineering, Delta University for Science and Technology, Gamasa, Egypt

Abstract

In recent decades, sea level has risen markedly along Mediterranean coastlines, highlighting the growing impacts on vulnerable coastal regions. This study analyzed sea-level variability using tide gauge (TG) records from 45 stations, coastal altimetry data (X-TRACK-L2P and X-TRACK/ALES), and virtual altimetry station (VAS) data. Based on TG records from 1993 to 2022, the results reveal positive linear trends, with 11 stations recording rates well above the global mean. The highest rate, 5.281 mm/year, was observed at TARANTO (Station ID: 2095). Strong consistency was also found between TG and X-TRACK-L2P coastal data, with X-TRACK-L2P trends ranging from −3 to 6 mm/year. In addition, sea-level rise in the Adriatic, Aegean, and Levantine sub-basins exceeded the current global mean. The agreement between TG and X-TRACK/ALES data was more limited, observed at only two sites located about 7 km and 20 km apart from the first valid point on the altimetry track, respectively. Furthermore, the findings highlight the diverse impacts of sea-level rise on Egypt's northern coast, where densely populated areas, productive agricultural lands, and cultural heritage sites are particularly at risk. In conclusion, all types of data used indicate positive linear trends in the Mediterranean Sea. In addition, several sites have rates of sea-level trend that exceed the global mean, underscoring the urgent need for proactive adaptation measures.

1 Introduction

In today's Earth science research, studying and forecasting sea-level changes at both the global and regional scales remain a challenging scientific issue (). The Public is becoming increasingly concerned about how climate change is affecting the rate of sea-level rise (SLR), especially as it affects densely populated coastal areas (; Puthucherril, 2025). Our ability to forecast future sea-level variability is essential for developing plans for coastal protection and adaptation projects under climate change scenarios. The dependability of models, assimilation systems, and current sea-level data all play a role in this ().

Numerous spatiotemporal scales show how sea level fluctuates, reflecting the superposition of large-scale processes such as changes in ocean circulation and global mean thermosteric changes imposed by variations in ocean temperature. It also reflects changes in ocean mass brought on by land-ice melt (). The global oceans exhibit diverse spatiotemporal sea-level structures, and a significant portion of the interannual to decadal sea-level variability is linked to climatic modes (Wang et al., 2021; Rusu, 2025). However, sea-level fluctuations around the coast can be significantly different from those in the open ocean (). This difference can be linked to distinct spatiotemporal features ().

The Mediterranean Sea, which is a semi-enclosed body of water connected to the North Atlantic Ocean through the Gibraltar Strait, has been identified as one of the world's most vulnerable climate zones (Pinardi et al., 2015). This vulnerable position is due to its numerous low-lying, densely populated coastal regions (). Coastal populations and activities may be greatly impacted by SLR and associated future projections (). Because of the local vertical land motion mechanisms, the low elevation may considerably worsen SLR in a number of coastal places throughout the Mediterranean Sea (Uzun and Otay, 2025). Sea-level trends in the Mediterranean deviate from world mean sea-level trends (Pinardi et al., 2014; ). They may even deviate from the values recorded in the adjacent Atlantic Ocean (). That is because non-linear local oceanographic processes can cause sea-level trends in the Mediterranean to deviate from the basin mean ().

The poor sampling of the two observational datasets, tide gauge (TG) and satellite altimetry (ALT) data, is the primary cause of the challenges in accurately measuring sea-level changes (Zhou et al., 2022). TG and ALT observations must first be expressed in the same reference frame in order to be compared (). TG data and continuous global navigation satellite system (GNSS) observations of the TG benchmark displacement can be combined to achieve this comparison (). It is also important to note that TG and ALT measurements are made at different sites (Ray et al., 2021). These may be farther apart than the spatial scales of local oceanic variability. Only up to 10–15 kilometers from the coast can traditional ALT products provide reliable sea-level readings (Woodworth et al., 2019). Such distance is insufficient in many areas to address the intricate ocean dynamics along the shoreline (Woodworth et al., 2019).

Literature addressing sea-level variability and patterns in the Mediterranean, through the comparison of TG and ALT data, typically examines a small number of sub-regional basin tide-gauge stations (e.g., ). Furthermore, related to the type of ALT data, some studies utilized along-track ALT data such as , while Taibi and Haddad (2019); , and employed multi-mission gridded ALT data. In addition, coastal ALT data were utilized to estimate trends (; ).

This study investigates sea-level changes in the Mediterranean Sea, presenting an updated time series of tidal data up to 2022. We assess how effective ALT data are in measuring coastal sea levels and examine the impacts of rising sea levels on Egypt's northern coast. To achieve this, we utilized a regional multi-mission coastal ALT dataset, including X-TRACK-L2P (the first application of this dataset in the Mediterranean) and X-TRACK/ALES data. These were used alongside all available TG records in the Mediterranean Sea throughout thirty years (1993–2022), making it the longest TG time series in the Mediterranean used during the era of ALT. Additionally, we incorporate virtual altimetry station (VAS) data to make up for the lack of TG stations along the African coast. Furthermore, we focused on Egypt's northern coast to evaluate the regional effects of these changes. That is because rising sea levels lead to significant threats, including coastal flooding, erosion of coastlines, flooding of dry land and low areas, and more saltwater entering rivers and groundwater (; ).

Section 2 describes the data used and the processing that went into them. Section 3 discusses the results of estimating the rates of change in sea level and compares different techniques based on coastal data. Section 4 illustrates the effects of SLR on Egypt's northern coast. Finally, Section 5 presents the conclusion.

2 Data and methods

2.1 Dataset selection

The types and sources of the data utilized in this study are described in the section that follows. These data comprise tide gauge records, satellite altimetry datasets, and virtual altimetry station data; each is explained in turn below.

2.1.1 Tide gauge data

We analyzed the TG time series of monthly mean sea level (MSL) utilizing the revised local reference (RLR), which has a documented datum history. The data were sourced from the Permanent Service for Mean Sea Level (PSMSL, ; PSMSL, 2025). The RLR datum for each site is around 7000 mm below the mean sea level. The PSMSL made this discretionary selection to prevent negative values in the resultant RLR monthly sea level-time series, which is recommended for scientific applications (Woodworth et al., 2019).

Only TGs with time series matching the altimeter period (1993–2022) and lasting at least 18 years during this period were chosen. They also had to have more than 75% of the monthly value available compared to the total months in the time series. These thresholds signify a balance between an adequate depiction of the altimeter period and maintaining reasonable geographical coverage of the Mediterranean coastline. Forty-nine tide gauges were picked from over one hundred available gauges. Figure 1 illustrates the IDs and positions of TGs along the Mediterranean coastline and within the principal sub-basins of the Mediterranean.

Figure 1

2.1.2 Vertical land motion data

The GNSS-related VLM was employed to adjust the TG trends. GNSS stations must be co-located with the TGs to ensure that the recorded VLM accurately reflects the measurements obtained from the TGs (). Particularly in areas with notable tectonic and human activity, the risk of different VLM between the GNSS and TG locations increases with their separation (Santamaría-Gómez et al., 2017). The VLM values employed to adjust the TG-derived relative sea-level time series were sourced from the Nevada Geodetic Laboratory (NGL). These data depend on the Median Interannual Difference Adjusted for Skewness (MIDAS, ; ). VLM data were also obtained from the SONEL website, which relies on GNSS data from processing centers at the University of La Rochelle (ULR7, ).

We imposed an arbitrary maximum distance of 30 km for TG-GNSS (about 85% of stations were placed with less than 5 km) and a minimum data series duration of 10 years. This was done to improve the representativeness of the data incorporated in TG measurements. However, for some stations, these conditions were not met, so glacial isostatic adjustment (GIA)-related VLM derived from the gridded ICE-6G_C (VM5a) model by Peltier et al. (2015) was used. These data are accessible on the PSMSL website. Figure 2 illustrates the temporal availability of TG and GNSS-related VLM measurements at the designated sites.

Figure 2

2.1.3 Altimetry data

Traditional satellite altimetry has developed into a technique for measuring ocean sea-level variations. This procedure was done by employing radar altimeters to ascertain sea surface height in relation to a reference ellipsoid (). The coastal altimetry method is a specialized form of traditional satellite altimetry that focuses on assessing sea-level changes along the shoreline (Stammer and Cazenave, 2017; Wang and Huang, 2021; ). We considered two different coastal ALT datasets, hereby identified as X-TRACK-L2P and X-TRACK/ALES data.

The X-TRACK-L2P (1 Hz) data comprises along-track sea level anomalies (SLA) mapped onto reference tracks, with a spatial interval of approximately 6–7 km between points (1 second). The SLA time series has been generated using the X-TRACK processing system, utilizing Level 2 Plus (L2P) products provided by Aviso+. The altimetric corrections implemented are included in this product. The X-TRACK-L2P data relied on the T/P, Jason-1, 2, and 3, T/P interleaved, Jason-1 and 2 interleaved, ERS-1 and 2, ENVISAT, and SARAL/AltiKa missions spanning from 1993 to 2022 (). In coastal zones, this technique works on modifying and post-processing the data. This practice improves the accuracy and completeness of sea surface height (SSH) measured by satellite altimetry ().

The X-TRACK/ALES (2 Hz) data is a product of coastal sea levels along-track, provided for six designated zones by the sea-level project of the ESA's Climate Change Initiative (The Climate Change Coastal Sea Level Team 2020). This dataset is generated by predicting the satellite altimetry range with the "Adaptive Leading Edge Sub-waveform Retracker" (ALES, Passaro et al., 2014) and further processed utilizing the X-TRACK methodology (; ). The X-TRACK/ALES data encompasses the three Jason-1, 2, and 3 missions from January 2002 to May 2018. The spatial resolution is 350 m at a frequency of 20 Hz, while the temporal resolution is 10 days (). Subsequently, this data was employed to investigate the ALT and TG measurement comparison zone of influence (), to observe sea-level changes on the West African coast (), to assess sea-level changes along the coast of New Zealand (), and in the Mediterranean Sea region (; ).

2.1.4 Virtual altimetry stations data

VAS does not correspond to physical locations on Earth, such as TG or meteorological stations. Rather, they represent virtual locations along satellite tracks from which researchers can obtain accurate sea level measurements by ALT (). A VAS is created by the trajectory's closest distance to the seashore of the first valid point. This dataset is derived from reprocessed high-resolution (300 m) along-track ALT data collected during the Jason-1, Jason-2, and Jason-3 missions, which spanned from early 2002 to 2021 (). reported SLA and corresponding coastal sea-level trends at more than 700 VAS situated along the global coasts, including the Mediterranean Sea, all within 6 km of the coast. Therefore, we utilized the VAS data to make up for the lack of TG stations along the African coast.

2.2 Data processing

2.2.1 Tide gauge data

To estimate sea-level trends depending on different techniques and to conduct a robust inter-technique comparison, it is essential to guarantee that the time series of the various observations achieve optimal consistency (Taibi and Haddad, 2019; ). 49 TG stations were chosen after meeting the criteria mentioned in the dataset selection part. Then, the modified seasonal Mann-Kendall Test (MKT) was utilized. This test is a non-parametric statistical test employed to identify monotonic trends (either increasing or decreasing) in time series data. It modifies the calculation of the standard test's variance to account for the presence of serial correlation in the data (; ). The null hypothesis (H0) in this test indicates that no trend exists in the data. This implies that the data are independently dispersed and random (; Taibi and Haddad, 2019).

For a data time series consisting of n points, denoted as x1, x2, …, xn, the MKT statistic (S) is calculated using Equations 1, 2, where S represents the difference between the number of increasing pairs and the number of declining pairs. A significant positive S implies an increasing trend, while a significant negative S denotes a decreasing trend ().

The seasonal MKT is a straightforward adaptation of the MKT that addresses observable seasonal variations in the dataset. This test conducts an independent MKT trend for each of the m distinct seasons, where m is the total number of seasons (). The assessment involves calculating the MKT statistic (Si) and its variation Var (Si) for each season using measurements gathered throughout the years. The cumulative (S') and the cumulative Var (S') are derived by aggregating the S statistic for each season as shown in Equations 3, 4 (; ).

The next step involves the most robust modification for the seasonal MKT, which utilizes a correction factor to account for serial correlation. The modified variance can be represented by Equation 5 (). The effective sample size represents the number of independent observations that are equivalent to the autocorrelated series. It is calculated as Equation 6 (; ).

Where:

  • - n is the actual number of observations.

  • - is the autocorrelation function at lag k of the ranks of the data, not the original data itself. This is a crucial point.

Finally, the computation of the Z' statistics is calculated as Equation 7 (; ).

Evaluate ∣∣ against the crucial value α/2 from the standard normal distribution at a specified level of significance α (e.g., 0.05 for 95% confidence). If ∣∣> α/2, reject H0 (the trend is significant). The p-value can be obtained from Z; if the p-value is less than α, the trend is significant (; ; ).

The modified seasonal MKT was conducted at a 5% level of significance for each dataset of the 49 TG stations. Subsequent to the test application, four time-series datasets were omitted: SOUDHAS station (ID 1232), ZADAR station (ID 1859), PALINURO station (ID 2082), and CATANIA II station (ID 2094), as indicated in Table 1. Consequently, 45 of the 49 stations are deemed acceptable following the test application, with all relevant application information shown in Table A in Supplementary Material.

Table 1

Station ID.Station nameS'Var(S') standardVar(S') modifiedZP-valueNull hypothesis (H0)
1232SOUDHAS37440734628191.4880.13669Not rejected
1859ZADAR471655662970.5800.56213Not rejected
2082PALINURO37440734628191.4880.13669Not rejected
2094CATANIA II641501474110.7320.46429Not rejected

The excluded TG stations after applying the modified seasonal Mann-Kendall test.

S' is the Seasonal Mann-Kendall statistic, Var(S') standard is its variation, and Var(S') modified is the modified variance which uses a correction factor to account for serial correlation. The (Equation 7) is utilized as an indicator of the trend's significance. The Null Hypothesis (H₀) in this test indicates that no trend exists in the data.

Following the application of the modified seasonal MKT to assess trends, the seasonal cycle was eliminated to detrend the time series. It was accomplished using the Multiple Seasonal Trend decomposition based on LOESS (MSTL) (). This method is commonly employed to separate time series of data into three parts: trend, seasonal, and residual, as delineated in Equation 8 (). Subsequently, the original trend was reintroduced into the deseasonalized time series.

where:

  • -Yt is the observed time series,

  • -Tt​ is the trend component,

  • -S1,t,S2,t,…, Sn,t are the seasonal components (one for each seasonal period),

  • -Rt​ is the remainder (residual) component.

In the absence of an error accompanying the data, trends are calculated utilizing the ordinary least squares (OLS) approach. This approach estimates the unknown parameters in a linear regression model (Zdaniuk, 2024). The objective is to identify the line (or hyperplane in higher dimensions) that optimally fits the data. That is done by reducing the total squared differences between the values that were observed and those that the model predicted as per Equation 9 ().

where:

  • -Y is the dependent variable (response),

  • -x1, x2,…, xp are the independent variables (predictors),

  • -β0, β1,…, βp are the regression coefficients,

  • -ϵ is the error term (residuals).

The information about the accepted TG stations is shown in Table B in Supplementary Material. Figure 3 shows the monthly SLA and trends for the different stations over the period 1993–2022.

Figure 3

The TG time series were subsequently adjusted for the long-term impact of VLM, utilizing the trends obtained from GNSS data or GIA data () as per Equation 10.

where:

  • -ASL: Absolute Sea level

  • -RSL: Relative Sea level (Tide gauge data)

  • -VLM: Vertical Land Motion

The results demonstrate that intense VLM significantly influences sea-level trends along some sites on the northeast coast of the Mediterranean Sea. The large number of earthquakes that occurred in the Mediterranean region throughout the study period may account for the robust values of VLM. Seismic activity persists in the Mediterranean region, with 9,712 earthquakes of magnitude 4 or more documented from January 1993 to December 2022, as per the US Geological Survey (USGS) database. The US Geological Survey (USGS) database documented approximately 709 earthquakes with magnitudes of 5 or more during this extensive timeframe, and Figure 4 shows their locations. Consequently, we observe that regions frequently impacted by earthquakes influence sea level rise. This phenomenon is evident in locations such as DUBROVNIK (ID: 760) with a VLM of -1.57 mm/yr, KHALKIS NORTH (ID: 1237), and KHALKIS SOUTH (ID: 1441) with a VLM of 1.14 mm/yr (Figure 4).

Figure 4

2.2.2 The X-TRACK-L2P altimetry data

The X-TRACK data, as detailed in , is utilized in conjunction with AVISO L2P data. Equation 11 illustrates the computation of the SLA associated with the X-TRACK data.

The linear trends were calculated via the least-squares technique. The spatial distribution of sea level trends, shown in Figure 5, offers an explanation for the dynamic variations obtained from about thirty years of monthly mean ALT observations. The trends vary from −3 to 6 mm/year, indicating beneficial variations along the whole coastal zone of the Mediterranean Sea. These results match those conducted by Taibi and Haddad (2019). The highest value, 6 mm/year, recorded in certain areas on the coastlines, surpasses 1.8 times the present global average of 3.3 ± 0.3 mm/year reported by . This value indicates acceleration in SLR. The Adriatic, Aegean, and Levantine subbasins have the largest positive rates (4 mm/yr and higher), which matches the results by and . These high trends occur particularly in areas with recurrent gyres and eddies in circulation, such as the Shikmona Gyre System (SGS), the Pelops Gyre (PG), and the Mersa Matruh Gyre System (MMGS, ; ). This observation is especially alarming as it exceeds worldwide forecasts and may result in accelerated coastal flooding, erosion, and the degradation of essential habitats. Conversely, some portions of the Levantine and Ionian sub-basins note negative trends (reaching -3 mm/yr) at considerable distances from the coastline. These adverse trends indicate limited oceanic phenomena, such as upwelling, alterations in dynamic topography, or sediment redistribution, rather than terrestrial influences ().

Figure 5

2.2.3 The X-TRACK/ALES dataset

produced the X-TRACK/ALES dataset using algorithms and geophysical corrections refined for ALT observations of coastal areas. However, combining these data with in-situ data is not straightforward, as Jason's tracks do not always meet the TG sites. Consequently, TG and X-TRACK/ALES records may denote sites that are tens of kilometers apart. Therefore, a condition was established to select points from the dataset that are at maximum 40 kilometers apart from both the TG and the nearest location on any track. Thus, nine TG stations were chosen, and trends were calculated using the least squares technique. Figure 6 shows the tracks spanning the Mediterranean region, accompanied by the locations of the TGs incorporated in our analysis.

Figure 6

3 Results and discussion

The sea-level trends of all selected TG stations were calculated based on Equation 10. All TG trends corrected to VLM have a statistically significant trend and show an increasing annual trend. Table B in Supplementary Material presents all the information about the TG data and their trend. The TARANTO II station (Station ID: 2095) has the highest yearly trend after the VLM adjustments, with a change rate of 5.281 mm/year. Conversely, the MALAGA II station (Station ID: 1810) exhibits the lowest yearly trend, of 0.696 mm/year. The mean trend across the stations is approximately 2.503 mm/year. This is consistent with the findings of , who reported that the sea level in the Mediterranean has fluctuated at a rate of 2.5 mm/year from 1993 to 2017.

The rate of sea-level rise in some regions exceeds the global average of 3.3 mm/year reported by . These regions are THESSALONIKI (ID: 373) at 4.197 mm/year, PIRAIEVS (ID: 374) at 4.183 mm/year, TARIFA (ID: 488) at 3.319 mm/year, KHALKIS NORTH (ID: 1237) at 3.995 mm/year, HADERA (ID: 1797) at 5.034 mm/year, BARCELONA (ID: 1811) at 3.823 mm/year, CORFU (ID: 1933) at 3.886 mm/year, PORTO TORRES (ID: 2084) at 3.925 mm/year, TARANTO II (ID: 2095) at 5.281 mm/year, REGGIO CALABRIA II (ID: 2142) at 3.651 mm/year, and ALMERÍA (ID: 2435) at 3.819 mm/year. These results are in good agreement with the findings from the common TG stations studied by ; Ramos-Alcántara et al. (2022), and . It should be noted that variations in the time series period lead to differences of approximately 1 mm/year at a few sites. Figure 7 illustrates the linear trend of the TG time series, the values of VLM, and the adjusted trends corresponding to VLM along with their uncertainty.

Figure 7

The trends obtained from the X-TRACK-L2P data were calculated for all 45 TGs. A comparison with the VLM-corrected TG trend was conducted, as illustrated in Table 2 and Figure 8.

Table 2

Station ID.Lat. (in Deg.)Long. (in Deg.)RSL (TG)GNSS station nameDistance of TG-GNSS (km)VLMASL(RSL+VLM) versus ALT
(RSL+VLM)ALT
6143.2795.3542.512MARS a0.007-0.7431.7692.3670.598
35243.50716.4421. 928(GIA)—–-0.1401.7882.8221.034
37340.63322.9354.283THS1 a2.142-0.0864.1973.429-0.768
37437.93723.6274.373(GIA)—–-0.1904.1832.484-1.699
40838.37226.1412.346(GIA)—–-0.0902.2563.2390.983
41038.95920.7573.423LFKD a14.355-0.1803.2432.681-0.562
48836.009-5.6033.989CEU1 b0.009-0.6703.3191.809-1.511
76042.65818.0632.714DUB2 a3.96-1.5671.1472.7221.575
76145.08313.6281.837PORE a16.075-1.1230.7142.9122.198
78843.7297.4211.957EZEV a7.932-0.8111.1462.1170.971
95843.3983.6993.212SETE a0.006-1.3121.92.8470.947
96038.339-0.4811.429ALAC a0-0.3271.1022.0130.911
98043.1135.9153.012(GIA)—–-0.3902.6222.265-0.356
123337.1326.8481.95(GIA)—–-0.1901.762.4400.680
123437.4424.9462.729(GIA)—–-0.2402.4892.320-0.169
123738.47223.5932.856HALK a1.1511.1393.9952.865-1.130
123840.84425.8781.633(GIA)—–-0.0101.6232.9271.304
124037.64521.322.045(GIA)—–-0.2801.7651.9140.149
144138.46123.5891.138HALK a1.2491.1392.2772.8570.580
146843.6967.2862.709NICE a3.5210.0272.7362.125-0.611
173535.8214.5332.248MALT a2.08-1.0761.1721.7570.584
176442.0543.2062.03GIRO a29.117-0.9611.0692.6471.578
179732.4734.8634.795CSAR a3.2180.2395.0343.438-1.597
181036.712-4.4171.451MLGA a1.689-0.7550.6962.0821.386
181141.3422.1664.193(GIA)—–-0.373.8232.635-1.188
181339.442-0.3112.872VCIA a2.894-0.5932.2792.3350.056
181745.54813.7251.985KOPE a0.001-0.3691.6162.9041.288
189239.5522.6392.121MALL b0.875-0.0202.1011.751-0.350
192941.9238.7632.016AJAC a0.5180.2552.2711.946-0.326
193238.9111.454.205IBIZ a0.076-1.4622.7431.812-0.931
193339.62819.9054.556KASI b13.372-0.6703.8862.825-1.061
207843.8788.0192.106EZEV b0.767-0.3601.7462.0930.347
207935.512.6043.767LAMP a0.111-0.4743.2932.238-1.055
208043.54610.2991.845(GIA)—–-0.2301.6152.2630.649
208440.8428.4044.515(GIA)—–-0.5903.9251.756-2.170
208640.67714.7513.507(GIA)—–-0.303.2072.345-0.861
208741.88816.1772.938(GIA)—–-0.2202.7182.591-0.127
208939.219.1142.689CAGL b1.891-0.2502.4391.687-0.752
209044.418.9262.765GENO a1.086-0.6242.1411.988-0.153
209540.47617.2245.541(GIA)—–-0.2605.2813.630-1.651
209640.14718.4972.03GIUR a6.246-0.4571.5732.7501.176
209742.35614.4153.503FRRA a12.214-1.9881.5152.4680.954
209843.62513.5071.47ANCG a2.4650.7782.2482.6190.371
214238.12215.6493.736TGRC a1.498-0.0853.6511.796-1.855
243536.832-2.4853.849ALME b0.002-0.0303.8192.245-1.575

Comparison between trends of corrected TG data and the X-TRACK 1 Hz trend (ALT) at the same locations (mm/year).

a is the VLM rate from NGL(MIDAS), while b is the VLM rate from SONEL (ULR7A).

Figure 8

The computed trends indicate positive linear trends in 1993–2022. Notably, higher rates of sea level increase were recorded in several locations, particularly THESSALONIKI (ID: 373), HADERA (ID: 1797), and TARANTO II (ID: 2095), where both TG records and X-TRACK-L2P data surpassed the world MSL trend. In some areas, such as the Aegean and Ionian seas, the rates of sea level rise significantly exceed the global mean. This pattern is consistent with the findings of , who addressed a high average rate of sea level rise across the entire Mediterranean Sea between 2000 and 2018.

Based on the values shown in Table 2, the statistics ​​were calculated. The TG mean trend was 2.503 mm/year, and the TG standard deviation was 1.132. The ALT mean trend was 2.408 mm/year, and the ALT standard deviation was 0.529. In addition, the absolute discrepancies between rates from ALT and TG measurements ranged from 0.056 to 2.229 mm/year, with a mean of 0.095 mm/year. Finally, the root mean square error (RMSE X-TRACK-L2P) was 1.134. The discrepancies between X-TRACK-L2P and TG data are consistent with those derived from the Copernicus Marine Environment Monitoring Service (CMEMS) and VLM-corrected TG data analyzed by and . These results sufficiently support the rising sea level over the Mediterranean Sea and the consistency between coastal sea level observations and in-situ data.

We also looked at whether the ALT product made for coastal use (X-TRACK/ALES) (20 Hz) matched better with TG observations than the X-TRACK-L2P (1 Hz) data. Therefore, as previously stated, we used the nearest point on the data track to determine the trend values at the TG locations. The sea-level trends were calculated for the selected nine points with approximately a complete time series of data. The minimum ALT-TG distance was 7.2 km at DUBROVNIK (ID: 760), and the maximum was 36.6 km at MARSAXLOKK (ID: 1735). To compare these trends with TG trends, the TG data had to be reanalyzed to match the XTRACK/ALES data period, which ran from 2002 to 2018. Table 3 displays the absolute values of the differences between the estimated trends from XTRACK/ALES data and in situ measurements at the nine sites.

Table 3

Station ID.Lat. (in Deg.)Long. (in Deg.)RSL (TG)GIA-Related VLMASLRSL+VLM versus ALT
(RSL+VLM)ALT (20Hz)
76042.65818.0632.051-0.141.9111.704-0.207
78843.7297.4211.835-0.291.5450.293-1.253
95843.3983.6992.479-0.232.2490.858-1.390
146843.6967.2863.655-0.293.3651.232-2.133
173535.82014.5332.923-0.512.4131.999-0.413
176442.0543.2061.779-0.391.3891.9990.610
207935.50012.6043.047-0.362.6870.288-2.399
208440.8428.4044.01-0.593.4201.036-2.384
208741.88816.1774.045-0.223.8253.030-0.795

Comparison between trends of corrected TG data to GIA-related VLM and the X-TRACK\ALES (ALT 20Hz) at the same locations (mm/year).

The statistics were calculated using data from 2002 to 2018. The TG mean trend was 2.534 mm/year, and the TG standard deviation was 0.862. While the ALT mean trend was 1.382 mm/year, and the ALT standard deviation was 0.893. In addition, the mean value of the differences between the trends from the two datasets was 1.152 mm/year, and the RMSE XTRACK/ALES was 1.517.

When comparing the trend differences between TG trends and both ALT datasets for the same nine sites (Figure 9), the X-TRACK-L2P data better matched the TG trends. This result may be due to the long measurement period of X-TRACK-L2P (from 1993 to 2022) and the corresponding TG records, which could better align with the SLR change (). This conclusion is also supported by its lower RMSE value for the nine stations (RMSEX-TRACK-L2P – 9 stations< RMSEX-TRACK/ALES). Therefore, the RMSE X-TRACK-L2P was recalculated for these nine TG datasets, giving a value of 1.144. However, X-TRACK/ALES exhibited even stronger trend agreement at only two sites: the in-site observations at DUBROVNIK (ID: 760), where the distance between TG and X-TRACK/ALES is smallest (about 7 km), which exactly matches results. The second site is L'ESTARTIT (ID: 1764), which is about 20 km (Figure 9).

Figure 9

When comparing the trend differences derived from XTRACK/ALES and TG data with those concluded by , it became evident that the results are in agreement. The small discrepancies (approximately 0.5 mm/year) arise from applying more stringent selection criteria for TG stations, which require at least 18 years of data and 75% completeness, as well as a TG-GNSS distance of no more than 30 km. In contrast, use less stringent TG selection criteria, requiring at least 14 years of data and 65% completeness.

Twenty-four VASs were used to estimate the change in sea level in the southern and eastern parts of the Mediterranean Sea to make up for the lack of TG stations. The trends of sea level were calculated using the least squares approach, depending on the values of monthly sea level anomalies provided for these stations. Furthermore, all the results show positive linear trends over the Mediterranean Sea, as shown in Figure 10. However, the issue with using the Virtual Altimetry Stations is their distance from the coast. That is because the first reliable points in the track that could be considered a VAS are located within 6 kilometers of the coast. Table C in Supplementary Material represents the information and trends on all the Virtual Altimetry Stations used in the study.

Figure 10

Although several research studies advance our understanding of the dynamics of the Mediterranean Sea, they utilize different approaches and data sources. All of this investigation leads to complementary views about the positive linear trends of the Mediterranean Sea and its impacts on the coastal areas, but with slight differences in some results. However, some research, like Taibi and Haddad (2019), relied on a shorter time frame of the TG data while examining 18 high-quality tide-gauge stations from 1993 to 2015. Several investigations, such as , also utilized other Alt datasets, focusing only on validating the X-TRACK/ALES (20 Hz) product. In addition, Ramos-Alcántara et al. (2022) and used CMEMS datasets. Furthermore, as noted in the dataset selection section, the station selection criteria in this study differ from those in earlier research. The stricter criteria applied here contribute to more accurate outcomes. Moreover, using a more updated dataset—which incorporates the X-TRACK-L2P product alongside X-TRACK/ALES, the longest available TG records (1993–2022) from 45 stations, and the VAS on the southern and eastern Mediterranean coastlines—enables a more robust and comprehensive assessment of three-decade trends across the entire Mediterranean basin.

Finally, it is noted that all the techniques used to calculate SLR rates in the Mediterranean Sea indicate positive linear trends. Some sites exhibit trends that surpass the global mean SLR, underscoring the necessity for proactive measures to address potential impacts of this phenomenon. The results also confirm the consistency of X-TRACK-L2P and TG data in representing the rates of SLR. For the southern and eastern coastline of the Mediterranean Sea, where no suitable tide gauge (TG) time series can be used, the X-TRACK-L2P or the VAS data can be relied upon. The RSL at any location without TG data can be estimated by subtracting the GNSS-related VLM from the ASL obtained from ALT data, provided that the closest GNSS station is used. If no nearby GNSS station is available, the GIA-related VLM can be used to correct the ALT data.

4 Impacts of SLR on Egypt

Rising sea levels, projected to cause flooding, are one of the most important challenges for the Mediterranean countries during this century (; Schuerch et al., 2025). According to scientists, up to 20 million people in the region could be permanently relocated by the end of the century if current trends in global warming and SLR continue (Hilmi et al., 2022). Egypt, located in northeastern Africa, shares its northern border with the Mediterranean Sea, as shown in Figure 11. Coastal regions at risk include Alexandria— a major industrial, commercial, urban center, and UNESCO World Heritage site on the Mediterranean—and the northern Nile Delta, where there is high population density. These areas face significant threats from coastal erosion and flooding driven by SLR (; ).

Figure 11

Egypt's Nile Delta is regarded as the most significant of all, because it supports more than half of Egypt's agricultural land and is home to more than 50% of the country's people (). Delta's low-lying zones make up half of its total area. 2.4% of Egypt's land area is below 2 meters above mean sea level, making it especially susceptible to flooding, storms, and increasing sea levels (Shaltout et al., 2015; ). A assessment from the International Panel on Climate Change (IPCC) states that rising sea levels, a consequence of global warming, are making Egypt's Nile Delta more susceptible to coastal flooding. According to the geospatial study that was done by Vecchio et al. (2024), the Mediterranean region's coastal plains have a total area of approximately 38,529 km² that could be vulnerable to flooding. France (3,681 km²), Italy (10,060 km²), and Egypt (12,879 km²) are the most affected nations. Allan et al. (2023) also note that by 2100, it is expected that 2,660 km² of the northern Delta will be underwater due to flooding (Vecchio et al., 2024).

According to Schuerch et al. (2025), the Nile Delta might lose thousands of acres suitable for investment, agriculture, or human habitation, with projections indicating a total loss of all coastal wetlands, even under optimistic climate change scenarios. showed that sea levels in Egypt have risen by 3.2 mm annually since 2012, which will continue to flood and erode the Delta's northern shore. Our results indicate higher rates of sea level rise, with the Delta's sea level trend measuring 3.88 and 3.70 mm/year, respectively, based on the VAS and X-TRACK-L2P data (Figure 11).

Additionally, sea level rise will worsen saltwater intrusion, which will be hastened by rising temperatures and force saltwater into the groundwater and soil that farmers use for irrigation (Wassef and Schüttrumpf, 2016). The lakes and reclaimed lowlands in the Nile Delta depend on the sand belt, which can be destroyed by rising sea levels (Wassef and Schüttrumpf, 2016). Additionally, the Nile Delta may subside at a rate of 0.5 cm annually, hastening the process of Delta submergence (Rateb and Abotalib, 2020).

A similar issue is present in coastal cities like Alexandria (Radwan et al., 2021). Therefore, the rate of SLR in Alexandria (station ID: 503) was calculated using the PSMSL dataset for the period 1944-2006. After correction for GIA-related VLM, the SLR was estimated at 1.478 mm/year. For the shorter period 1993-2006, the SLR increased to 4.39 mm/year, indicating an acceleration in sea-level rise. In addition, the SLR was calculated from X-TRACK-L2P data, and it is estimated to be 3.00 mm/year (Figure 11), which exceeds the mean SLR in the Mediterranean Sea. According to experts, this is because the soil lacks silt, which solidifies the soil and is being held back more and more by the Nile Dam at Aswan, Upper Egypt. Therefore, everything related to water is at risk since rising sea levels brought on by global warming cause coastal cities to flood ().

emphasize that sea-level rise, in conjunction with vertical land movement, exacerbates tsunami risks in the Mediterranean by elevating relative sea levels. Therefore, in Alexandria, residential and commercial properties are at risk of tidal, wind-driven, and storm surge flooding almost annually, and occasionally multiple times throughout the year, during times of intense rainfall, strong winds, and sea level rise ().

At best, according to other Egyptian and UN studies, "30% of Alexandria will be flooded, 1.5 million or more people will be displaced, 195,000 jobs will be destroyed, and losses in land and construction will reach $30 trillion" if the Mediterranean rises just 50 centimeters (). In addition to its history and its vestiges of the past, Alexandria is also home to the largest port in the country, which is the hub of the Egyptian economy; thus, the catastrophe will have an impact on approximately one hundred million Egyptians ().

The Mediterranean is home to various UNESCO World Heritage Sites (WHS) that face threats from floods and erosion due to rising sea levels (). Reimann et al. (2018) assembled datasets that identify UNESCO cultural World Heritage Sites in the Mediterranean's Low-Elevation Coastal Zone (LECZ) as of 2018. Reimann et al. (2018) observed that several heritage sites, in regions facing sea level rise far beyond the global average, are particularly susceptible. Figure 12a shows the distribution of the UNESCO WHS on the Mediterranean coast. It also shows that several TG and VAS stations documented sea-level rise rates beyond the global mean where they are located in high population density areas. In addition, these stations are in proximity to multiple WHS, heightening the hazards of flooding and the potential loss of invaluable cultural landmarks. Furthermore, the historic coastal city of Alexandria in Egypt is facing increasing dangers that may result in the loss of cultural assets and substantial effects on tourism (). Reimann et al. (2018) provide and illustrate all details about the Mediterranean UNESCO WHS, including their location in urban areas, site elevation, category, and distance from the coast.

Figure 12

Numerous prior investigations, such as ; , and Vousdoukas et al. (2022), demonstrated that Alexandria's coastal heritage sites will be impacted by sea inundation and disturbances in groundwater flow. These sites are low-lying archaeological locations such as the Qaitbay Citadel, Ras El Tin Palace, Chatby Tombs, and Anfoushy Necropolis. The four locations, as shown in Figure 12b, are situated at distances of about 25, 117, 180, and 300 meters from the sea, making them very susceptible to fluctuations in sea level (). The infiltration of seawater at these locations will also impact on the archaeological structures, thus necessitating the prompt implementation of essential protective actions.

5 Conclusion

Our findings emphasize the critical importance of understanding and monitoring sea level rise (SLR) along the Mediterranean coastline. In this region, SLR has been shaped by various processes, including the complex interplay between ocean circulation and climate variability. This study provides comprehensive information on sea-level variations along the Mediterranean coast between 1993 and 2022, utilizing a multifaceted methodology that incorporates both contemporary coastal ALT and TG data. It also demonstrates the agreement between these two techniques and emphasizes the implications of sea level rise (SLR) on Egypt's coastlines, using them as a case study.

We examined TG observations from the ALT era (1993–2022), which represents the longest period used in the analysis of TG data along the Mediterranean coast since satellite altimetry began. The network included 45 TG sites that provide time series spanning at least 18 years, with a minimum of 75% data completeness, and located within 30 kilometers of the nearest GNSS station. A modified seasonal Mann–Kendall test (MKT) was applied to detect significant patterns. Four TG datasets with otherwise suitable time series had to be excluded due to limitations in applying for the test. Sea-level trends, estimated using ordinary least squares (OLS) after removing the seasonal cycle with MSTL, showed rising sea levels at all TG sites with significant trends.

In our dataset, we established a maximum TG-GNSS distance of 30 km. In addition, a minimum data series duration of 10 years was set to improve the representativeness of the data included in TG measurements. Otherwise, we depended on GIA-related VLM. We observed that the sea level trends at certain locations on the northeast coast of the Mediterranean Sea are markedly affected by substantial vertical land movement (VLM). That brought us back to the significant frequency of earthquakes that occurred in the Mediterranean region, particularly in the Adriatic and Aegean sub-basins, during the research period. Upon rectifying the trends of TG to VLM, it was observed that all TG stations exhibit a statistically significant trend, characterized by an upward annual trajectory. The average trend among the stations is roughly 2.503 mm per year. Moreover, the rate of sea level increase in certain areas surpassed the global average of 3.3 mm/year, highlighting the pressing need for preventative actions to lessen the effects of SLR.

This study also confirmed the accuracy of coastal X-TRACK-L2P data by comparing it with TG records at selected sites. We mapped the regional distribution of sea-level trends from the X-TRACK-L2P dataset and their uncertainties. The trends ranged from −3 to 6 mm/year, showing notable variations across the Mediterranean coastline. Interannual variability revealed a strong correlation between coastal X-TRACK-L2P and TG data, with absolute differences of 0.056 to 2.229 mm/year (average of 0.095 mm/year), and the RMSEX-TRACK-L2P of 1.134. In several locations, SLR rates exceeded the global mean, particularly at THESSALONIKI (ID: 373), HADERA (ID: 1797), and TARANTO II (ID: 2095), where both TG and ALT values were above the global mean trend. Similarly, the Adriatic, Aegean, and Levantine subbasins also displayed SLR rates well above the global mean.

We also utilized the X-TRACK/ALES coastal data with only 9 TG stations located within 40 kilometers of the TG locations. The X-TRACK/ALES (20 Hz) trends were calculated for the 9 TGs, and a comparison with the VLM-corrected TG trend was also applied. It is noted that during the period (2002–2018), X-TRACK-L2P data exhibited superior concordance with in-situ data regarding the discrepancies between their trends than X-TRACK/ALES. However, TG and X-TRACK/ALES showed a higher level of agreement at only two locations, where the distances between TG locations and the closest point on the altimetry track were about 7 km and 20 km, respectively.

Concerning the southern and eastern Mediterranean coastline, we depended on the X-TRACK-L2P data, which has exhibited strong agreement with the TG data and the Virtual Altimetry Stations (VAS). Ultimately, all techniques employed to assess SLR rates in the Mediterranean Sea demonstrate positive linear trends, which points to the need for preventative actions to mitigate the SLR's possible effects.

Finally, we investigated the trend of sea level along Egypt's northern coast and noticed that it had risen above the Mediterranean Sea's mean, which threatens the low-lying zones of Egypt's Delta that are ideal for investment, agriculture, and human habitation. In addition, Alexandria's coastal heritage sites, particularly low-lying archaeological locations such as the Qaitbay Citadel, RAS ELTIN Palace, Chatby Tombs, and Anfoushy Necropolis, will be impacted by sea inundations and disturbances in groundwater flow.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. The monthly tide-gauge time series are available from the PSMSL (https://psmsl.org/data/obtaining/). The vertical land motion values are sourced from the Nevada Geodetic Laboratory (NGL, https://geodesy.unr.edu/PlugNPlayPortal.php) using MIDAS velocities https://geodesy.unr.edu/velocities/midas.IGS14.txt), and from the SONEL website (https://www.sonel.org/-Vertical-land-movements-.html?lang=en). The utilization of GIA-related VLM derived from the gridded ICE-6G_C (VM5a) model, accessible on the PSMSL website (https://psmsl.org/train_and_info/geo_signals/gia/peltier/index.php). The X-TRACK-L2P SLA data was obtained from AVISO+ (https://www.aviso.altimetry.fr/en/data/products/sea-surface-height-products/regional/x-track-sla.html), while XTRACK/ALES data v1.1 are available at https://catalogue.ceda.ac.uk/uuid/222cf11f49a94d2da8a6da239df2efc4, https://doi.org/10.5270/esasl_cci-xtrack_ales_sla-200206_201805-v1.1-202006. The network of virtual altimetry stations v2.4 are available at https://www.seanoe.org/data/00631/74354/. The locations of UNESCO cultural World Heritage in the Mediterranean coastal zone are available at https://figshare.com/articles/dataset/UNESCO_cultural_World_Heritage_in_the_Mediterranean_coastal_zone/5759538. Finally, seismic activity locations in the Mediterranean region are obtained from the US Geological Survey (USGS) database (https://earthquake.usgs.gov/earthquakes/search/).

Author contributions

MY: Investigation, Conceptualization, Data curation, Writing – original draft, Visualization, Formal analysis, Methodology, Validation. QL: Writing – review & editing, Visualization, Investigation, Methodology, Data curation. AZ: Validation, Methodology, Data curation, Writing – review & editing. LW: Funding acquisition, Writing – review & editing, Methodology, Validation. LB: Investigation, Conceptualization, Writing – review & editing, Methodology, Supervision, Funding acquisition, Formal analysis.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This work was supported in part by the National Natural Science Foundation of China (Grant Nos. 42192535, 42174102, 42274116, 42404004) and the Natural Science Foundation of Wuhan (2024040701010065).

Acknowledgments

All data suppliers are appreciated by the authors. The authors also would like to thank the editor and the reviewers, whose comments and suggestions helped improve and clarify this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2025.1683307/full#supplementary-material

References

  • 1

    AllanR. P.AriasP. A.BergerS.CanadellJ. G.CassouC.ChenD.et al. (2023). “ Intergovernmental panel on climate change (IPCC). Summary for policymakers,” in Climate change 2021: The physical science basis. Contribution of working group I to the sixth assessment report of the intergovernmental panel on climate change (Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press), 332. doi: 10.1017/9781009157896.001

  • 2

    AndroulidakisY.PytharoulisI. (2025). Variability of marine heatwaves and atmospheric cyclones in the Mediterranean Sea during the last four decades. Environ. Res. Lett.20, 034031. doi: 10.1088/1748-9326/adb505

  • 3

    BandaraK.HyndmanR. J.BergmeirC. (2025). MSTL: A seasonal-trend decomposition algorithm for time series with multiple seasonal patterns. Int. J. Oper. Res. 52 (1), 7998. doi: 10.1504/IJOR.2025.143957

  • 4

    Bellas-ManleyA.NeremR. S.HamlingtonB. D. (2025). Extrapolation of the satellite altimeter record to understand regional variations in future sea level change. J. Geophys. Res.: Oceans130, e2024JC022094. doi: 10.1029/2024JC022094

  • 5

    BirolF.DelebecqueC. (2014). Using high sampling rate (10/20 Hz) altimeter data for the observation of coastal surface currents: A case study over the northwestern Mediterranean Sea. J. Mar. Syst.129, 318333. doi: 10.1016/j.jmarsys.2013.07.009

  • 6

    BirolF.FullerN.LyardF.CancetM.NinoF.DelebecqueC.et al. (2017). Coastal applications from nadir altimetry: Example of the X-TRACK regional products. Adv. Space Res.59, 936953. doi: 10.1016/j.asr.2016.11.005

  • 7

    BirolF.LégerF.PassaroM.CazenaveA.NiñoF.CalafatF. M.et al. (2021). The X-TRACK/ALES multi-mission processing system: New advances in altimetry towards the coast. Adv. Space Res.67, 23982415. doi: 10.1016/j.asr.2021.01.049

  • 8

    BlewittG.KreemerC.HammondW. C.GazeauxJ. (2016). MIDAS robust trend estimator for accurate GPS station velocities without step detection. J. Geophys. Res.: Solid Earth121, 20542068. doi: 10.1002/2015JB012552

  • 9

    BonaduceA.PinardiN.OddoP.SpadaG.LarnicolG. (2016). Sea-level variability in the Mediterranean Sea from altimetry and tide gauges. Climate Dyn.47, 28512866. doi: 10.1007/s00382-016-3001-2

  • 10

    BruniS.FenoglioL.RaicichF.ZerbiniS. (2022). On the consistency of coastal sea-level measurements in the Mediterranean Sea from tide gauges and satellite radar altimetry. J. Geodesy96, 41. doi: 10.1007/s00190-022-01626-9

  • 11

    CalafatF. M.FrederikseT.HorsburghK. (2022). The sources of sea-level changes in the Mediterranean Sea since 1960. J. Geophys. Res.: Oceans127, e2022JC019061. doi: 10.1029/2022JC019061

  • 12

    CazenaveA.CozannetG. L. (2014). Sea level rise and its coastal impacts. Earth's Future2, 1534. doi: 10.1002/2013EF000188

  • 13

    CazenaveA.GouzenesY.BirolF.LegerF.PassaroM.CalafatF. M.et al. (2022). Sea level along the world’s coastlines can be measured by a network of virtual altimetry stations. Commun. Earth Environ.3, 117. doi: 10.1038/s43247-022-00448-z

  • 14

    DiengH. B.CazenaveA.GouzenesY.SowB. A. (2021). Trends and inter-annual variability of altimetry-based coastal sea level in the Mediterranean Sea: Comparison with tide gauges and models. Adv. Space Res.68, 32793290. doi: 10.1016/j.asr.2021.06.022

  • 15

    DurandF.PiecuchC. G.BeckerM.PapaF.RajuS. V.KhanJ. U.et al. (2019). Impact of continental freshwater runoff on coastal sea level. Surveys Geophys.40, 14371466. doi: 10.1007/s10712-019-09536-w

  • 16

    ErkoçM. H.DoğanU. (2023). Monitoring of recently sea level changes on the coast of New Zealand using XTRACK coastal altimetry and tide gauge data. Continent. Shelf Res.269, 105127. doi: 10.1016/j.csr.2023.105127

  • 17

    Fox-KemperB. (2021). “ Ocean, cryosphere and sea level change,” in AGU fall meeting abstracts, vol. 2021. New Orleans, LA, U13BU109. doi: 10.1017/9781009157896.011

  • 18

    FrederikseT.AdhikariS.DaleyT. J.DangendorfS.GehrelsR.LandererF.et al. (2021). Constraining 20th-century sea-level rise in the South Atlantic Ocean. J. Geophys. Res.: Oceans126, e2020JC016970. doi: 10.1029/2020jc01697

  • 19

    FrihyO. E.El-SayedM. K. (2013). Vulnerability risk assessment and adaptation to climate change induced sea level rise along the Mediterranean coast of Egypt. Mitig. Adapt. Strateg. Global Change18, 12151237. doi: 10.1007/s11027-012-9418-y

  • 20

    GravelleM.WöppelmannG.GobronK.AltamimiZ.GuichardM.HerringT.et al. (2023). The ULR-repro3 GPS data reanalysis and its estimates of vertical land motion at tide gauges for sea level science. Earth Syst. Sci. Data15, 497509. doi: 10.5194/essd-15-497-2023

  • 21

    GrezioA.AnzideiM.ArmigliatoA.BaglioneE.MaramaiA.SelvaJ.et al. (2025). “ Tsunami hazard and risk in the Mediterranean Sea,” in Probabilistic Tsunami Hazard and Risk Analysis ( Elsevier), 397415. doi: 10.1016/B978-0-443-18987-6.00015-4

  • 22

    GuérouA.MeyssignacB.PrandiP.AblainM.RibesA.Bignalet-CazaletF. (2023). Current observed global mean sea level rise and acceleration estimated from satellite altimetry and the associated measurement uncertainty. Ocean Sci.19, 431451. doi: 10.5194/egusphere-2022-330

  • 23

    HamedK. H.RaoA. R. (1998). A modified Mann-Kendall trend test for autocorrelated data. J. Hydrol.204, 182196. doi: 10.1016/S0022-1694(97)00125-X

  • 24

    HammondW. C.BlewittG.KreemerC.NeremR. S. (2021). GPS imaging of global vertical land motion for studies of sea level rise. J. Geophys. Res.: Solid Earth126, e2021JB022355. doi: 10.1029/2021JB022355

  • 25

    HamzawyA.Al-MailamM.ArkehJ. (2023). Climate change in Egypt: Opportunities and obstacles ( Carnegie Endowment for International Peace). Available online at: https://coilink.org/20.500.12592/6d61qd (Accessed June 25, 2025).

  • 26

    HemedaS. (2021). Geotechnical modelling of the climate change impact on world heritage properties in Alexandria, Egypt. Heritage Sci.9, 73. doi: 10.1186/s40494-021-00547-8

  • 27

    HilmiN.AliE.Carnicer ColsJ.CramerW.GeorgopoulouE.Le CozannetG.et al. (2022). “ IPCC AR6 WGII cross-chapter paper 4: Mediterranean region,” in EGU General Assembly, Vienna, Austria, 23–27 May 2022, EGU2210590. doi: 10.5194/egusphere-egu22-10590

  • 28

    HirschR. M.SlackJ. R.SmithR. A. (1982). Techniques of trend analysis for monthly water quality data. Water Resour. Res.18, 107121. doi: 10.1029/WR018i001p00107

  • 29

    HolgateS. J.MatthewsA.WoodworthP. L.RickardsL. J.TamisieaM. E.BradshawE.et al. (2013). New data systems and products at the permanent service for mean sea level. J. Coast. Res.29, 493504. doi: 10.2112/JCOASTRES-D-12-00175.1

  • 30

    JinY.ChenM.YanH.WangT.YangJ. (2023). Sea level variation in the Arctic Ocean since 1979 based on ORAS5 data. Front. Mar. Sci.10. doi: 10.3389/fmars.2023.1197456

  • 31

    KamalI.FekriM.Abou El-MagdI.SolimanN. (2021). Mapping the impacts of projected sea-level rise on Cultural heritage sites in Egypt: Case study (Alexandria). J. Fac. Tour5, 120. doi: 10.21608/MFTH.2021.190350

  • 32

    Kemgang GhomsiF. E.ErkoçM. H.RajR. P.PirtiA.BonaduceA.AbiodunB. J.et al. (2025). Sea level trends along the South African coast from 1993 to 2022 using XTRACK altimetry, tide gauges, and GNSS measurements. Sci. Rep.15, 4946. doi: 10.1038/s41598-025-89258-9

  • 33

    LeclercqL.CazenaveA.LegerF.BirolF.NinoF.ToluL.et al. (2025). Coastal sea level rise at altimetry-based virtual stations in the Gulf of Mexico. Adv. Space Res.75, 16361652. doi: 10.1016/j.asr.2024.11.069

  • 34

    LegeaisJ. F.MeyssignacB.FaugèreY.GuerouA.AblainM.PujolM. I.et al. (2021). Copernicus sea level space observations: a basis for assessing mitigation and developing adaptation strategies to sea level rise. Front. Mar. Sci.8. doi: 10.3389/fmars.2021.704721

  • 35

    LiY.FengJ.YangX.ZhangS.ChaoG.ZhaoL.et al. (2024). Analysis of sea level variability and its contributions in the Bohai, Yellow Sea, and East China Sea. Front. Mar. Sci.11. doi: 10.3389/fmars.2024.1381187

  • 36

    LornezhadE.EbrahimiH.RabieifarH. R. (2023). Analysis of precipitation and drought trends by a modified Mann–Kendall method: a case study of Lorestan province, Iran. Water Supply23, 15571570. doi: 10.2166/ws.2023.068

  • 37

    MarcosM.WöppelmannG.MatthewsA.PonteR. M.BirolF.ArdhuinF.et al. (2019). Coastal sea level and related fields from existing observing systems. Surveys Geophys.40, 12931317. doi: 10.1007/s10712-019-09513-3

  • 38

    MartiF.CazenaveA.BirolF.PassaroM.LégerF.NiñoF.et al. (2021). Altimetry-based sea level trends along the coasts of western Africa. Adv. Space Res.68, 504522. doi: 10.1016/j.asr.2019.05.033

  • 39

    Masson-DelmotteV.ZhaiP.PiraniA.ConnorsS. L.PéanC.BergerS.et al. (2021). “ Climate change 2021: the physical science basis,” in Contribution of working group I to the sixth assessment report of the intergovernmental panel on climate change 2, no. 1Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press, 2391. doi: 10.1017/9781009157896

  • 40

    MeliM.CamargoC. M.OlivieriM.SlangenA. B.RomagnoliC. (2023). Sea-level trend variability in the Mediterranean during the 1993–2019 period. Front. Mar. Sci.10. doi: 10.3389/fmars.2023.1150488

  • 41

    MohamedB.AbdallahA. M.Alam El-DinK.NagyH.ShaltoutM. (2019). Inter-annual variability and trends of sea level and sea surface temperature in the Mediterranean Sea over the last 25 years. Pure Appl. Geophys.176, 37873810. doi: 10.1007/s00024-019-02156-w

  • 42

    OelsmannJ.CalafatF. M.PassaroM.HughesC.RichterK.PiecuchC.et al. (2024). Coherent modes of global coastal sea level variability. J. Geophys. Res.: Oceans129, e2024JC021120. doi: 10.1029/2024JC021120

  • 43

    OelsmannJ.PassaroM.DettmeringD.SchwatkeC.SánchezL.SeitzF. (2021). The zone of influence: matching sea level variability from coastal altimetry and tide gauges for vertical land motion estimation. Ocean Sci.17, 3557. doi: 10.5194/os-17-35-202

  • 44

    PandžićK.LiksoT.BiondićR.BiondićB. (2024). A review of the contribution of satellite altimetry and tide gauge data to evaluate sea level trends in the Adriatic Sea within a Mediterranean and Global Context. GeoHazards5, 112141. doi: 10.3390/geohazards5010006

  • 45

    PassaroM.CipolliniP.VignudelliS.QuartlyG. D.SnaithH. M. (2014). ALES: A multi-mission adaptive subwaveform retracker for coastal and open ocean altimetry. Remote Sens. Environ.145, 173189. doi: 10.1016/j.rse.2014.02.008

  • 46

    PeltierW. R.ArgusD. F.DrummondR. (2015). Space geodesy constrains ice age terminal deglaciation: The global ICE-6G_C (VM5a) model. J. Geophys. Res.: Solid Earth120, 450487. doi: 10.1002/2014JB011176

  • 47

    Permanent Service for Mean Sea Level (PSMSL) (2025). Tide Gauge Data. Available online at: http://www.psmsl.org/data/obtaining/ (Accessed June 10, 2025).

  • 48

    PinardiN.BonaduceA.NavarraA.DobricicS.OddoP. (2014). The mean sea level equation and its application to the Mediterranean Sea. J. Climate27, 442447. doi: 10.1175/JCLI-D-13-00139.1

  • 49

    PinardiN.ZavatarelliM.AdaniM.CoppiniG.FratianniC.OddoP.et al. (2015). Mediterranean Sea large-scale low-frequency ocean variability and water mass formation rates from 1987 to 2007: A retrospective analysis. Prog. Oceanogr.132, 318332. doi: 10.1016/j.pocean.2013.11.003

  • 50

    PuthucherrilT. G. (2025). Adapting to sea level rise: Is India on- or off-track? Front. Mar. Sci.12. doi: 10.3389/fmars.2025.1516241

  • 51

    RadwanA. M.MagdyM.RabahM.SaberA.ZakiA. (2021). Sea level analysis using tide gauge observations at the northern Delta coast, Egypt. NRIAG J. Astron. Geophys.10, 361371. doi: 10.1080/20909977.2021.1940742

  • 52

    Ramos-AlcántaraJ.GomisD.JordàG. (2022). Reconstruction of Mediterranean coastal sea level at different timescales based on tide gauge records. Ocean Sci.18, 17811803. doi: 10.5194/os-18-1781-2022

  • 53

    RatebA.AbotalibA. Z. (2020). Inferencing the land subsidence in the Nile Delta using Sentinel-1 satellites and GPS between 2015 and 2019. Sci. Total Environ.729, 138868. doi: 10.1016/j.scitotenv.2020.138868

  • 54

    RayR. D.LoomisB. D.ZlotnickiV. (2021). The mean seasonal cycle in relative sea level from satellite altimetry and gravimetry. J. Geodesy95, 80. doi: 10.1007/s00190-021-01529-1

  • 55

    ReimannL.VafeidisA. T.BrownS.HinkelJ.TolR. S. J. (2018). Mediterranean UNESCO World Heritage at risk from coastal flooding and erosion due to sea-level rise. Nature communications9 (1), 4161. doi: 10.1038/s41467-018-06645-9

  • 56

    RusuL. (2025). Climate change impact on the sea state conditions in the Mediterranean Sea under RCP and SSP emission scenarios. Renewable Energy243, 122616. doi: 10.1016/j.renene.2025.122616

  • 57

    Santamaría-GómezA.GravelleM.DangendorfS.MarcosM.SpadaG.WöppelmannG. (2017). Uncertainty of the 20th century sea-level rise due to vertical land motion errors. Earth Planet. Sci. Lett.473, 2432. doi: 10.1016/j.epsl.2017.05.038

  • 58

    SchuerchM.KieselJ.BoutronO.GuelmamiA.WolffC.CramerW.et al. (2025). Large-scale loss of Mediterranean coastal marshes under rising sea levels by 2100. Commun. Earth Environ.6, 128. doi: 10.1038/s43247-025-02099-2

  • 59

    ShaltoutM.TonbolK.OmstedtA. (2015). Sea-level change and projected future flooding along the Egyptian Mediterranean coast. Oceanologia57, 293307. doi: 10.1016/j.oceano.2015.06.004

  • 60

    StammerD.CazenaveA. (Eds.) (2017). Satellite altimetry over oceans and land surfaces (Boca Raton: CRC Press), 343380. doi: 10.1201/9781315151779

  • 61

    TaibiH.HaddadM. (2019). Estimating trends of the Mediterranean Sea level changes from tide gauge and satellite altimetry data, (1993–2015). J. Oceanol. Limnol.37, 11761185. doi: 10.1007/s00343-019-8164-3

  • 62

    The Climate Change Coastal sea level team. (2020). ESA Sea Level Climate Change Initiative (Sea_Level_cci): Altimeter along-track high resolution sea level anomalies in some coastal regions (2002-2018) from the JASON satellites, v1.1. Centre for Environmental Data Analysis, date of citation. doi: 10.5270/esa-sl_cci-xtrack_ales_sla-200206_201805-v1.1-202005

  • 63

    UzunP.OtayE. N. (2025). Sea-level change in Eastern Mediterranean, Aegean and Black seas analyzed at subregional scale of six CORDEX domains. Regional Stud. Mar. Sci.83, 104099. doi: 10.1016/j.rsma.2025.104099

  • 64

    VecchioA.AnzideiM.SerpelloniE. (2024). Sea level rise projections up to 2150 in the northern Mediterranean coasts. Environ. Res. Lett.19, 014050. doi: 10.1088/1748-9326/ad127e

  • 65

    VousdoukasM. I.ClarkeJ.RanasingheR.ReimannL.KhalafN.DuongT. M.et al. (2022). African heritage sites threatened as sea-level rise accelerates. Nat. Climate Change12, 256262. doi: 10.1038/s41558-022-01280-1

  • 66

    WangJ.ChurchJ. A.ZhangX.ChenX. (2021). Reconciling global mean and regional sea level change in projections and observations. Nat. Commun.12, 990. doi: 10.1038/s41467-021-21265-6

  • 67

    WangH.HuangZ. (2021). Waveform decontamination for improving satellite radar altimeter data over nearshore area: upgraded algorithm and validation. Front. Earth Sci.9. doi: 10.3389/feart.2021.748401

  • 68

    WassefR.SchüttrumpfH. (2016). Impact of sea-level rise on groundwater salinity at the development area western Delta, Egypt. Groundwater Sustain. Dev.2, 85103. doi: 10.1016/j.gsd.2016.06.001

  • 69

    WoodworthP. L.MeletA.MarcosM.RayR. D.WöppelmannG.SasakiY. N.et al. (2019). Forcing factors affecting sea level changes at the coast. Surveys Geophys.40, 13511397. doi: 10.1007/s10712-019-09531-1

  • 70

    ZdaniukB. (2024). “ Ordinary least-squares (OLS) model,” in Encyclopedia of quality of life and well-being research ( Springer International Publishing, Cham), 48674869. doi: 10.1007/978-3-031-17299-1_2008

  • 71

    ZhouD.LiuY.FengY.ZhangH.FuY.LiuY.et al. (2022). Absolute sea level changes along the coast of China from tide gauges, GNSS, and satellite altimetry. J. Geophys. Res.: Oceans127, e2022JC018994. doi: 10.1029/2022JC018994

Summary

Keywords

tide gauge data, X-TRACK dataset, virtual altimetry stations, sea-level trends, the Mediterranean region, Egypt's low-lying coastal zones

Citation

Youssef MM, Li Q, Zaki A, Wu L and Bao L (2025) Estimating the Mediterranean sea-level trends from 1993 to 2022 showing the consistency of X-TRACK-L2P coastal altimetry with tide gauge data and trend impacts on Egypt's northern coast. Front. Mar. Sci. 12:1683307. doi: 10.3389/fmars.2025.1683307

Received

10 August 2025

Accepted

06 October 2025

Published

31 October 2025

Volume

12 - 2025

Edited by

David Alberto Salas de León, National Autonomous University of Mexico, Mexico

Reviewed by

Kreso Pandzic, Retired, Zagreb, Croatia; Muhammad Luqman Ahmad Affandi, Universiti Teknologi Malaysia, Malaysia

Updates

Copyright

*Correspondence: Lifeng Bao,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics