Caspian Sea Level Change Observed by Satellite Altimetry

: We analyze satellite altimeter observed Caspian Sea level (CSL) changes over the period January 1993 to December 2021 using the lake level series from the Hydroweb project and global sea level anomalies (SLA) grids provided by the Copernicus Marine Environment Monitoring Service (CMEMS). The two altimeter-based CSL series agree well at interannual and longer time scales, but show signiﬁcantly large discrepancies at seasonal and shorter time scales. The large discrepancies are found to be introduced by the approximately inverted barometer (IB) correction applied to the CMEMS SLA over the Caspian Sea. The IB correction over the Caspian Sea or any enclosed lakes needs to be treated separately from the ocean by using the correct reference mean pressure. The actual IB effects over the Caspian Sea are signiﬁcantly smaller than those applied in the CMEMS SLA grids. After applying an improved IB correction using the global mean sea level pressure ﬁelds from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis atmospheric model, the two CSL series agree remarkably well. Altimeter observed CSL series show a signiﬁcant decreasing trend on top of strong seasonal variations. The estimated linear trends for the Hydroweb and CMEMS CSL series are − 5.37 ± 0.11 and − 5.40 ± 0.11 cm/yr, respectively. Annual amplitudes are 17.03 ± 1.33 vs. 15.79 ± 1.30 cm, with nearly the same phases. The CSL change shows notable acceleration in the decreasing trend since around 2005, and the estimated trends have increased to − 8.86 ± 0.10 and − 8.81 ± 0.10 cm/yr, respectively for the two-altimeter CSL series.


Introduction
The Caspian Sea is the largest lake on Earth with an area of about 371,000 km 2 .Over the past few hundred years, the Caspian Sea level (CSL) has experienced substantial fluctuations up to several meters [1][2][3].The Caspian Sea is located within an endorheic basin between Europe and Asia, and the CSL change is mainly controlled by water mass exchange between the Caspian Sea and Caspian drainage basin (see Figure 1) via river discharge, precipitation, and evaporation.Over the past three decades, the CSL shows a large and steady decreasing trend on top of strong seasonal fluctuations, which is believed to be driven by imbalanced water fluxes or increased evaporation over the Caspian Sea due to global warming [3][4][5].The CSL is projected to fall by 9-18 m by the end of this century under the medium to high emissions scenarios [6], which is expected to cause catastrophic impacts on the Caspian Sea coastal regions, especially in the northern part where most of the water depths are less than only 5-10 m [5].If this were to happen, the Kara-Bogaz-Gol (KBG) Bay would dry up completely.
. Accurate measurements of CSL change play a critical role in understanding iability at different temporal scales and connections with driving forces in the gl regional climate system and predicting future CSL change and potential impac ecosystem in the Caspian Sea region.Tide gauge observations are also important itoring CSL change, especially at long-term time scales (from multi-decadal to a years).However, due to special data policies in the surrounding countries, mo all) tide gauge observations along the Caspian Sea coasts are not available to the public.Other challenges in using and interpreting tide gauge observations inc inconsistent temporal coverages, data qualities, and vertical land movements of gauge stations.
Satellite remote sensing has become an advanced observational tool to mo global ocean since the late 70s.Satellite-borne Synthetic Aperture Radar techn been used since 1978 for sea surface monitoring, as well as for mapping applica 9].Satellite radar altimeter experiments were also started at about the same time, accuracies of the early altimeter sea level observations at only meter level, inadeq meaningful sea level change studies.
Since the launch of the TOPEX/Poseidon (T/P) high-accuracy satellite radar a mission in August 1992, modern satellite altimeters have provided continuous m Accurate measurements of CSL change play a critical role in understanding CSL variability at different temporal scales and connections with driving forces in the global and regional climate system and predicting future CSL change and potential impacts on the ecosystem in the Caspian Sea region.Tide gauge observations are also important for monitoring CSL change, especially at long-term time scales (from multi-decadal to a hundred years).However, due to special data policies in the surrounding countries, most (if not all) tide gauge observations along the Caspian Sea coasts are not available to the general public.Other challenges in using and interpreting tide gauge observations include the inconsistent temporal coverages, data qualities, and vertical land movements of the tide gauge stations.
Satellite remote sensing has become an advanced observational tool to monitor the global ocean since the late 70s.Satellite-borne Synthetic Aperture Radar technique has been used since 1978 for sea surface monitoring, as well as for mapping applications [7][8][9].Satellite radar altimeter experiments were also started at about the same time, with the accuracies of the early altimeter sea level observations at only meter level, inadequate for meaningful sea level change studies.
Since the launch of the TOPEX/Poseidon (T/P) high-accuracy satellite radar altimeter mission in August 1992, modern satellite altimeters have provided continuous measurements of global sea level (including large lakes water level) changes with unprecedented accuracy for over 30 years [10,11].Altimeter-observed global sea level changes are commonly presented as global gridded sea level anomalies (SLA), defined as residual sea level changes after effects from ocean tides and atmospheric loading have been removed [12,13].
Estimating large lake water level change from a satellite altimeter involves some additional challenges due to the small spatial scales and different dynamical processes compared to the ocean [14].Specially calibrated water level changes for over hundreds of large lakes in the world are available from the Hydroweb project (http://hydroweb.theia-land.fr/;accessed on 24 October 2022) developed by LEGOS (Laboratoire d'Etudes en Géophysique et Océanographie Spatiales) and CNES (Centre National d'Etudes Spatiales) [14,15].
In the present study, we carry out an analysis of CSL variations over the period of January 1993 to December 2021 using multi-satellite altimeter measurements.Two monthly CSL change series are included in the analysis.One is the specially calibrated CSL change series provided by the Hydroweb project and the other is estimated from the global SLA grids distributed by the Copernicus Marine Environment Monitoring Service (CMEMS, https://marine.copernicus.eu/;accessed on 27 July 2022).The main goals are to (1) study CSL variations at different time scales during the past 30 years; (2) compare CSL changes from the two data sources; and (3) investigate potential discrepancies between the two altimeter-based CSL estimates and identify the causes for the discrepancies.

Altimeter Observed CSL Changes
The Hydroweb project has been providing satellite altimeter-observed water level time series of selected major lakes in the world in near real-time via an automated algorithm by combing observations from multi-satellites altimeter missions, including T/P, Jason-1, -2, and -3 series, Geosat Follow-On, Envisat, and others [14,15].The original Hydroweb CSL series is irregularly sampled (with intervals from ~10 days in early years to less than a day in recent years).The Hydroweb CSL series has been resampled to monthly intervals by taking the mean of all estimates for a given month.Water level variations for the Caspian Sea and KBG are separately provided from Hydroweb.KBG is connected to the Caspian Sea via a narrow channel, and KBG shows similar but notably different variabilities compared to the Caspian Sea.In the present analysis, we focus on the Caspian Sea, with KBG only included in the inverted barometer (IB) related analysis.The monthly Hydroweb CSL series over the period of January 1993 to December 2021 is shown in Figure 2 (blue curves).
Global monthly mean SLA estimates given on 0.25 • × 0.25 • grids and covering the period January 1993 to December 2021 are provided by CMEMS.The CSL anomaly estimates are also included in the CMEMS SLA grids.The CMEMS global SLA grids are also derived from multi-satellite altimeter missions (similar to the Hydroweb lake level data).Monthly CMEMS CSL changes are computed from the CMEMS SLA grids over the Caspian Sea with the cosine of latitude as weighting.To be consistent with the Hydroweb CSL series, KBG is also excluded in the CMEMS CSL computation.As mentioned above, effects from ocean tides and atmospheric loading (i.e., IB -Inverted Barometer-) have been removed from the CMEMS SLA grids.Since the Caspian Sea is an enclosed water body, in theory, the ocean tide and IB corrections should not affect the estimated mean CSL change.The CMEMS monthly CSL change series (red curves) is compared with the Hydroweb estimates (blue curves) in Figure 2.Both CSL change series exhibit similar decreasing trends and other long-term var bilities over the 30-year period, with the trends appearing to accelerate since around 20 However, the CMEMS CSL series show significantly smaller seasonal variability, as co pared to the Hydroweb estimates.To better understand the discrepancy, we estimate nual and semiannual variations using least-squares fit and remove them from both C series and compare the two nonseasonal CSL series in Figure 3.It becomes more obvio -the two altimeter CSL series agree remarkably well at interannual time scales.In anoth word, it confirms that the discrepancy is mostly at seasonal time scales.The estimat linear trends for the Hydroweb and CMEMS CSL series are -5.37±0.11and -5.40±0 cm/yr, respectively (virtually the same).However, the estimated mean annual amplitu of the Hydroweb CSL series is 17.03±1.33cm, which is twice as large as that (8.72±1.30c from the CMEMS CSL series, although the annual phases are nearly the same (266±4 265±9 degree).Table 1 summarizes the comparisons of annual and semiannual amplitud and phases and linear trends from the two-altimeter CSL estimates (see rows 1-2).Annual and semiannual amplitudes and phases and linear trends from the twoaltimeter CSL estimates are estimated from least-squares fit using the following equation, In which, H and t represent the altimeter monthly CSL series and time epochs; a 1 and a 1 represent the annual and semiannual amplitudes; and ϕ 1 and ϕ 2 represent the annual and semiannual phases.r represents the linear trend and b is the mean.
Both CSL change series exhibit similar decreasing trends and other long-term variabilities over the 30-year period, with the trends appearing to accelerate since around 2005.However, the CMEMS CSL series show significantly smaller seasonal variability, as compared to the Hydroweb estimates.To better understand the discrepancy, we estimate annual and semiannual variations using least-squares fit and remove them from both CSL series and compare the two nonseasonal CSL series in Figure 3.It becomes more obvious-the two altimeter CSL series agree remarkably well at interannual time scales.In another word, it confirms that the discrepancy is mostly at seasonal time scales.The estimated linear trends for the Hydroweb and CMEMS CSL series are −5.37 ± 0.11 and −5.40 ± 0.11 cm/yr, respectively (virtually the same).However, the estimated mean annual amplitude of the Hydroweb CSL series is 17.03 ± 1.33 cm, which is twice as large as that (8.72 ± 1.30 cm) from the CMEMS CSL series, although the annual phases are nearly the same (266 ± 4 vs. 265 ± 9 degree).Table 1 summarizes the comparisons of annual and semiannual amplitudes and phases and linear trends from the two-altimeter CSL estimates (see rows 1-2).

Atmospheric Loading Effects on CSL Change
Given the particularly large magnitude (8-9 cm) and frequency band (seasonal) of the discrepancy between the two altimeter CSL change estimates, one likely cause may be the inconsistent treatments of atmospheric loading effects on CSL change through the IB correction.In the CMEMS altimeter global SLA grids, sea level change due to atmospheric loading is removed using the IB assumption, by adjusting the sea level to changes in bar-

Atmospheric Loading Effects on CSL Change
Given the particularly large magnitude (8-9 cm) and frequency band (seasonal) of the discrepancy between the two altimeter CSL change estimates, one likely cause may be the inconsistent treatments of atmospheric loading effects on CSL change through the IB correction.In the CMEMS altimeter global SLA grids, sea level change due to atmospheric loading is removed using the IB assumption, by adjusting the sea level to changes in barometric pressure as [16,17], ∆H = −0.9948× P − P ocn (2) in which, ∆H is the IB correction to sea level in units of cm, P is sea surface pressure in units of mb (mBar) at a given grid point, and P ocn (in mb) is the global mean sea surface pressure over the oceans, also called reference mean pressure.In the case of full adjustment, an increase in barometric pressure of 1 mb corresponds to a fall in sea level of ~1 cm.
In the IB correction, this pressure-driven sea level change needs to be removed from altimeter observations [12].Equation (2) shows that the integration of IB effects over the global ocean is zero, which means the IB effect only changes the spatial viability of the sea (or water) level, but not the global mean sea level change (or mean CSL change).The same IB Equation ( 2) can be applied to large lakes as well (for small lakes, the IB effect is negligible).However, for large lakes, the IB correction needs to be implemented separately, i.e., treating each enclosed lake as an individual ocean surface.The reference mean pressure P needs to be defined as the mean surface pressure over the specific lake (e.g., Caspian Sea) as In the Hydroweb altimeter lake-level data processing, the IB effect is not considered [14].Considering IB or not will give the same result, as the derived mean lake level change is not affected by the IB effect.This means that if the IB correction was implemented correctly in the CMEMS SLA grids, the atmospheric loading effect should not be responsible for the observed discrepancy between the Hydroweb and CMEMS CSL series.
Upon the clarifications from the CMEMS Sea Level Thematic Center (Françoise Mertz, Mercator-Ocean, personal communications), during the CMEMS SLA data processing, an approximate IB correction was applied for the Caspian Sea, in which the reference mean pressure used was the same as for the ocean, i.e., global mean sea level pressure over the ocean (P ocn ).The difference between the approximate reference mean pressure (P ocn ) and the actual reference mean pressure (P casp ) is expected to lead to a significant difference in the IB corrections and thus SLA estimates.

Improved Quantification of IB Effects on CSL Change
To better understand the potential impacts of the approximate IB correction on CMEMS CSL estimates, we use the global mean sea level pressure fields from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis [18] to quantify the CMEMS approximate and actual IB effects on Caspian Sea SLA and CSL change.ERA5 is the latest generation (5th) of the ECMWF reanalysis system developed through the Copernicus Climate Change Service [19].The ERA5 monthly mean sea level pressure fields on 0.25 • × 0.25 • grids over the period January 1993 to December 2021 are provided by the Copernicus Climate Data Store (https://cds.climate.copernicus.eu/;accessed on 24 November 2022).Monthly mean sea level pressure changes over the global ocean and the Caspian Sea are computed from the ERA5 sea level pressure fields and compared in Figure 4.Not surprisingly, the two mean sea level pressure series or the two reference mean pressure series (P ocn and P casp ) for the IB correction differ substantially.The seasonal variability of the global mean sea level pressure over the ocean is much smaller than that over the Caspian Sea.The large differences (up to over 10 mb) in the two mean sea level pressure series will be translated into up to over 10 cm differences in the derived SLA.For each grid point over the Caspian Sea, the actual and CMEMS approximate IB corrections are computed from the ERA5 sea level data using  ̅  and  ̅  as the reference mean pressure, respectively.Annual amplitudes at each grid point from the two IB corrections are estimated using an unweighted least-squares fit and shown in Figures 5a  and 5b (note the different color scales used for the two maps).Annual amplitudes of the actual IB corrections over the Caspian Sea are mostly less than 1 cm, while the annual amplitudes of the CMEMS approximate IB corrections are about 7-8 cm.Some spatial variability exists in both cases, reflecting regional pressure variations over the Caspian Sea.Two grid points, A (51°E, 46°N) and B (52°E, 38°N) in the north and south Caspian Sea are selected and marked by the black "+" in Figure 5a for IB affect time series comparisons shown in Figure 6.The selections of the two points are arbitrary, to simply illustrate the slightly spatial variability of atmospheric surface pressure over the Caspian Sea (between the north and south in this case).For each grid point over the Caspian Sea, the actual and CMEMS approximate IB corrections are computed from the ERA5 sea level data using P casp and P ocn as the reference mean pressure, respectively.Annual amplitudes at each grid point from the two IB corrections are estimated using an unweighted least-squares fit and shown in Figure 5a,b (note the different color scales used for the two maps).Annual amplitudes of the actual IB corrections over the Caspian Sea are mostly less than 1 cm, while the annual amplitudes of the CMEMS approximate IB corrections are about 7-8 cm.Some spatial variability exists in both cases, reflecting regional pressure variations over the Caspian Sea.Two grid points, A (51 • E, 46 • N) and B (52 • E, 38 • N) in the north and south Caspian Sea are selected and marked by the black "+" in Figure 5a for IB affect time series comparisons shown in Figure 6.The selections of the two points are arbitrary, to simply illustrate the slightly spatial variability of atmospheric surface pressure over the Caspian Sea (between the north and south in this case).After applying the improved or actual IB corrections to the CMEMS SLA grids over the Caspian Sea, we recompute the CMEMS CSL change series and compare it with the Hydroweb CSL series in Figure 7.The two altimeter CSL change results agree remarkably well and the discrepancies noted earlier are almost completely gone.The annual and semiannual amplitudes and phases and linear trend of the recomputed CMEMS CSL change series are estimated from unweighted least-squares fit and listed in Table 1 (row 3) for comparisons.With the improved IB correction, the mean annual amplitude of the CMEMS CSL is increased from 8.72±1.30cm to 15.79±1.30cm, and agrees well with the Hydroweb estimates (17.03±1.33 cm).The RMS of the residuals between the Hydroweb and CMEMS CSL estimates is reduced from 6.88 cm to 3.14 cm.As mentioned above, the IB correction, if implemented correctly, does not affect the CSL change estimates.To compute CMEMS CSL change, we can equivalently undo or remove the approximate IB correction applied to the CMEMS SLA grids.
At seasonal scales, the CMEMS approximate IB corrections (see the blue curves in Figures 6a & 6b) are in phase with the observed CSL changes (see Figure 2), all reaching the minimums around the start of the year.This is due to the fact the mean atmospheric pressure over the Caspian Sea is coincidentally out of phase with the observed CSL change, which also explains that the approximate or mis-modeled IB corrections only affect the amplitudes, but not the phase of the CSL change.After applying the improved or actual IB corrections to the CMEMS SLA grids over the Caspian Sea, we recompute the CMEMS CSL change series and compare it with the Hydroweb CSL series in Figure 7.The two altimeter CSL change results agree remarkably well and the discrepancies noted earlier are almost completely gone.The annual and semiannual amplitudes and phases and linear trend of the recomputed CMEMS CSL change series are estimated from unweighted least-squares fit and listed in Table 1 (row 3) for comparisons.With the improved IB correction, the mean annual amplitude of the CMEMS CSL is increased from 8.72 ± 1.30 cm to 15.79 ± 1.30 cm, and agrees well with the Hydroweb estimates (17.03 ± 1.33 cm).The RMS of the residuals between the Hydroweb and CMEMS CSL estimates is reduced from 6.88 cm to 3.14 cm.As mentioned above, the IB correction, if implemented correctly, does not affect the CSL change estimates.To compute CMEMS CSL change, we can equivalently undo or remove the approximate IB correction applied to the CMEMS SLA grids.

Conclusions and Discussion
We analyze satellite altimeter observed CSL changes over the period of January 1993 to December 2021 using the specially calibrated lake level data from the Hydroweb project and the global SLA grids provided by CMEMS.The two CSL series agree very well at interannual and longer time scales, but show significantly large discrepancies at seasonal and shorter time scales.The discrepancies are found to be introduced by the approximate IB corrections applied to the CMEMS SLA grids over the Caspian Sea.After applying the correct IB corrections using the ECMWF ERA5 sea surface pressure fields, the two altimeter-based CSL series agree remarkably well.
Over the 30-year period, satellite altimeter observed CSL series show a significant decreasing trend on top of strong seasonal variations.Decadal variations are also evident in CSL change.The estimated linear trends are −5.37±0.11and −5.40±0.11cm/yr for the Hydroweb and CMEMS CSL series, respectively.The annual amplitudes and phases of the two estimates are 17.03±1.33cm/266±4 degrees vs. 15.79±1.30cm/267±5 degrees.The CSL series show notable acceleration in the decreasing trend.For the period since 2005, the trends have increased to −8.86±0.10 and −8.81±0.10cm/yr, respectively for the two altimeter CSL series.
The Caspian Sea is an enclosed water body separated from the ocean.The IB correction over the Caspian Sea or any enclosed lakes needs to be treated separately from the ocean by using the correct reference mean pressure.The actual IB effects over the Caspian Sea are significantly smaller than the approximate estimate in the CMEMS SLA grid data product.The IB effects only affect spatial variability of sea level (or lake water level), but At seasonal scales, the CMEMS approximate IB corrections (see the blue curves in Figure 6a,b) are in phase with the observed CSL changes (see Figure 2), all reaching the minimums around the start of the year.This is due to the fact the mean atmospheric pressure over the Caspian Sea is coincidentally out of phase with the observed CSL change, which also explains that the approximate or mis-modeled IB corrections only affect the amplitudes, but not the phase of the CSL change.

Conclusions and Discussion
We analyze satellite altimeter observed CSL changes over the period of January 1993 to December 2021 using the specially calibrated lake level data from the Hydroweb project and the global SLA grids provided by CMEMS.The two CSL series agree very well at interannual and longer time scales, but show significantly large discrepancies at seasonal and shorter time scales.The discrepancies are found to be introduced by the approximate IB corrections applied to the CMEMS SLA grids over the Caspian Sea.After applying the correct IB corrections using the ECMWF ERA5 sea surface pressure fields, the two altimeter-based CSL series agree remarkably well.
Over the 30-year period, satellite altimeter observed CSL series show a significant decreasing trend on top of strong seasonal variations.Decadal variations are also evident in CSL change.The estimated linear trends are −5.37 ± 0.11 and −5.40 ± 0.11 cm/yr for the Hydroweb and CMEMS CSL series, respectively.The annual amplitudes and phases of the two estimates are 17.03 ± 1.33 cm/266 ± 4 degrees vs. 15.79 ± 1.30 cm/267 ± 5 degrees.The CSL series show notable acceleration in the decreasing trend.For the period since 2005, the trends have increased to −8.86 ± 0.10 and −8.81 ± 0.10 cm/yr, respectively for the two altimeter CSL series.
The Caspian Sea is an enclosed water body separated from the ocean.The IB correction over the Caspian Sea or any enclosed lakes needs to be treated separately from the ocean by using the correct reference mean pressure.The actual IB effects over the Caspian Sea are significantly smaller than the approximate estimate in the CMEMS SLA grid data product.The IB effects only affect spatial variability of sea level (or lake water level), but not the average sea or water level change of the global ocean or entire water body.
Satellite altimetry has been an established space geodetic technique for measuring and monitoring the global sea level change for over three decades since the launch of T/P, the first modern satellite altimeter mission.However, accurate quantification of SLA or water level changes in lakes using a satellite altimeter is still challenging.For large lakes, such as the Caspian Sea, altimeter-derived water level changes (or anomalies) are subject to larger uncertainties due to relatively smaller spatial scales, difficulties in modeling the tides, and special implementations of atmospheric loading corrections.Climate change is taking a big toll on the Caspian Sea and the surrounding region.A long record of continuous and accurate measurements of CSL change from satellite altimetry is important for people to better understand CSL change and its connections with climate change.

Figure 1 .
Figure 1.The Caspian Sea and the Caspian Sea drainage basin (circled by the red boundar with five surrounding countries (Russia, Azerbaijan, Iran, Turkmenistan, and Kazakhstan

Figure 1 .
Figure 1.The Caspian Sea and the Caspian Sea drainage basin (circled by the red boundary curves), with five surrounding countries (Russia, Azerbaijan, Iran, Turkmenistan, and Kazakhstan) labeled.

Figure 2 .
Figure 2. Monthly Caspian Sea level changes (in units of cm) observed by satellite altimeter over period of January 1993 to December 2021, provided by the Hydroweb project (https:// droweb.theia-land.fr/;Accessed on October 24 2022) and estimated from the CMEMS global S product (https://data.marine.copernicus.eu/;Accessed on July 27 2022).The two lines (red and bl represent the linear trends estimated from least-squares fit (the two lines are mostly overlapp with each other, as the estimated trends are too close).The estimated trends, −5.36 cm/yr and −5 cm/yr are labeled in the plot with corresponding colors.

Figure 2 .
Figure 2. Monthly Caspian Sea level changes (in units of cm) observed by satellite altimeter over the period of January 1993 to December 2021, provided by the Hydroweb project (https://hydroweb.theia-land.fr/;accessed on 24 October 2022) and estimated from the CMEMS global SLA product (https://data.marine.copernicus.eu/;accessed on 27 July 2022).The two lines (red and blue) represent the linear trends estimated from least-squares fit (the two lines are mostly overlapped with each other, as the estimated trends are too close).The estimated trends, −5.36 cm/yr and −5.41 cm/yr are labeled in the plot with corresponding colors.

12 Figure 3 .
Figure 3. Nonseasonal Caspian Sea level changes (in units of cm) observed by satellite altimeter over the period of January 1993 to December 2021, provided by the Hydroweb project (https://hydroweb.theia-land.fr/;Accessed on October 24 2022) and CMEMS global sea level anomaly product (https://data.marine.copernicus.eu/;Accessed on July 27 2022).The two lines (red and blue) represent the linear trends estimated from least-squares fit (the two lines are mostly overlapped with each other, as the estimated trends are too close).The estimated trends, −5.36 cm/yr, and −5.41 cm/yr are labeled in the plot with corresponding colors.

Figure 3 .
Figure 3. Nonseasonal Caspian Sea level changes (in units of cm) observed by satellite altimeter over the period of January 1993 to December 2021, provided by the Hydroweb project (https:// hydroweb.theia-land.fr/;accessed on 24 October 2022) and CMEMS global sea level anomaly product (https://data.marine.copernicus.eu/;accessed on 27 July 2022).The two lines (red and blue) represent the linear trends estimated from least-squares fit (the two lines are mostly overlapped with each other, as the estimated trends are too close).The estimated trends, −5.36 cm/yr, and −5.41 cm/yr are labeled in the plot with corresponding colors.

Figure 4 .
Figure 4. Monthly mean sea level pressure (in units of mb) over the global ocean (blue curve) and Caspian Sea (red curve) over the period January 1993 to December 2021, estimated from the ECMWF ERA5 mean sea level pressure data.

Figure 4 .
Figure 4. Monthly mean sea level pressure (in units of mb) over the global ocean (blue curve) and Caspian Sea (red curve) over the period January 1993 to December 2021, estimated from the ECMWF ERA5 mean sea level pressure data.

Figure 5 .
Figure 5. (a) Mean annual amplitudes (in cm) of the actual IB effects in Caspian Sea level ch over the period January 1993 to December 2021 based on the ECMWF ERA5 atmospheric data; (b) Mean annual amplitudes of the CMEMS approximate IB effects in Caspian Sea changes over the period January 1993 to December 2021 (also based on the ECMWF ERA5 pheric model data).Note the different color scales used for the two maps.Two grid points A 46°N) and B (52°E, 38°N) are marked by the black "+" in (a) for later comparisons.

Figure 5 .
Figure 5. (a) Mean annual amplitudes (in cm) of the actual IB effects in Caspian Sea level changes over the period January 1993 to December 2021 based on the ECMWF ERA5 atmospheric model data; (b) Mean annual amplitudes of the CMEMS approximate IB effects in Caspian Sea level changes over the period January 1993 to December 2021 (also based on the ECMWF ERA5 atmospheric model data).Note the different color scales used for the two maps.Two grid points A (51 • E, 46 • N) and B (52 • E, 38 • N) are marked by the black "+" in (a) for later comparisons.

Figure 6 .
Figure 6.Comparisons of the actual and CMEMS approximate IB effects in two selected points (a) A and (b) B in the Caspian Sea over the period January 1993 to December 2021 based on the ECMWF ERA5 atmospheric model data (Points A and B are marked by the "+" in Figure 5a).

Figure 6 .
Figure 6.Comparisons of the actual and CMEMS approximate IB effects in two selected points (a) A and (b) B in the Caspian Sea over the period January 1993 to December 2021 based on the ECMWF ERA5 atmospheric model data (Points A and B are marked by the "+" in Figure 5a).

Figure 7 .
Figure 7. Monthly Caspian Sea level changes (in units of cm) observed by satellite altimeter over the period of January 1993 to December 2021, provided by the Hydroweb project (https://hydroweb.theia-land.fr/;Accessed on October 24 2022) and estimated from the CMEMS global sea level anomaly product (https://data.marine.copernicus.eu/;Accessed on July 27 2022) with the CMEMS approximate IB correction removed and replaced by the appropriate one.The two lines (red and blue) represent the linear trends estimated from least-squares fit (the two lines are mostly overlapped with each other, as the estimated trends are too close).The estimated trends, −5.36 cm/yr, and −5.41 cm/yr are labeled in the plot with corresponding colors.

Figure 7 .
Figure 7. Monthly Caspian Sea level changes (in units of cm) observed by satellite altimeter over the period of January 1993 to December 2021, provided by the Hydroweb project (https: //hydroweb.theia-land.fr/;accessed on 24 October 2022) and estimated from the CMEMS global sea level anomaly product (https://data.marine.copernicus.eu/;accessed on 27 July 2022) with the CMEMS approximate IB correction removed and replaced by the appropriate one.The two lines (red and blue) represent the linear trends estimated from least-squares fit (the two lines are mostly overlapped with each other, as the estimated trends are too close).The estimated trends, −5.36 cm/yr, and −5.41 cm/yr are labeled in the plot with corresponding colors.

Table 1 .
Annual and semiannual amplitudes and phases and trends of the Hydroweb and CMEMS CSL altimetry series for January 1993 to December 2021 estimated from unweighted least-squares fit.

Table 1 .
Annual and semiannual amplitudes and phases and trends of the Hydroweb and CMEMS CSL altimetry series for January 1993 to December 2021 estimated from unweighted least-squares fit.