Evolution of the Performances of Radar Altimetry Missions from ERS-2 to Sentinel-3A over the Inner Niger Delta

Radar altimetry provides unique information on water stages of inland hydro-systems. In this study, the performance of seven altimetry missions, among the most commonly used in land hydrology (i.e., European Remote-Sensing Satellite-2 (ERS-2), ENVIronment SATellite (ENVISAT), Satellite with Argos and ALtika (SARAL), Jason-1, Jason-2, Jason-3 and Sentinel-3A), are assessed using records from a dense in situ network composed of 19 gauge stations in the Inner Niger Delta (IND) from 1995 to 2017. Results show an overall very good agreement between altimetry-based and in situ water levels with correlation coefficient (R) greater than 0.8 in 80% of the cases and Root Mean Square Error (RMSE) lower than 0.4 m in 48% of cases. Better agreement is found for the recently launched missions such as SARAL, Jason-3 and Sentinel-3A than for former missions, indicating the advance of the use of the Ka-band for SARAL and of the Synthetic-aperture Radar (SAR) mode for Sentinel-3A. Cross-correlation analysis performed between water levels from the same altimetry mission leads to time-lags between the upstream and the downstream part of the Inner Niger Delta of around two months that can be related to the time residence of water in the drainage area.


Introduction
Surface waters, which are part of the continental branch of the terrestrial water cycle, play an essential role in supplying fresh water for basic human and economic needs.They are strongly impacted by climate changes and anthropogenic pressures caused by population growth and changes in agricultural practices [1][2][3].Despite the importance of their monitoring for addressing integrated water resource management, use in operational flood forecasting or disaster mitigation, reliable in situ measurements of water stage and discharge has become increasingly scarce information due to either the disappearance of the gauge networks or the difficulty to get access to data [4,5].
Satellite radar altimetry, initially developed for the measurement of the ocean surface topography through the measurement of the distance between the Earth's surface and the spaceborne radar altimeter [6], has demonstrated its efficiency for deriving water levels of inland water bodies (see Crétaux et al. [7] for a recent review).Radar altimetry, was initially used over land to retrieve water levels over homogeneous surfaces such as large lakes and enclosed seas [8,9], but also at cross-sections between rivers and altimetry ground-tracks of several kilometers of width in large river basins [10,11].These early results were obtained using Geosat and Topex/Poseidon (T/P) ranges (i.e., the distance between the satellite and the surface) derived from the Ocean retracking algorithm.Root Mean Square Errors (RMSE) lower than 0.05 m and 1.1 m through comparisons with in situ water stages were obtained over lakes and rivers respectively.The comparisons are performed at the so-called Virtual Stations.Virtual stations (VS) are defined as the cross-sections of an altimetry ground-track and a water body (i.e., lake, reservoir river channel, floodplain, or wetland) where the temporal variations of the height from one cycle to the next can be associated with changes in water level [7].With the launch of ENVISAT in 2002, ranges processed using other retracking algorithms were included in the Geophysical Data Records (GDR) made available by space agencies.Among them, the Offset Center Of Gravity (OCOG, also known as Ice-1) was found to provide, most of the time, the most accurate estimate of river water levels (with RMSE generally lower than 0.3 m and correlation coefficient R greater than 0.9) [12].Combined with the availability of land-dedicated corrections of the ionosphere, wet troposphere delays and improvements in the data processing, this allowed the generalization of the use of radar altimetry for the monitoring of inland waters [7,13].
Thanks to these different improvements and the use of high-frequency data (10,18,20 or 40 Hz depending on the altimetry mission) instead of 1 Hz data (~7 km of sampling along the track) as over the open ocean water bodies of a few or below one hundred meters of width can now be monitored with very good accuracy (e.g., [14][15][16]).As radar altimetry data have global coverage and are freely available, they are now commonly used in a wide range of hydrological applications ( see Crétaux et al. [7] for a recent review), and, even in support for the management of in situ networks [17].Until now, no study provided a systematic assessment of the performance over rivers of different altimetry missions that were operating since the beginning of the high-precision altimetry era, which started with the launch of Topex/Poseidon (T/P) in 1992, contrary to that done over lakes [18,19].
The goal of this study is to evaluate the quality of altimetry-based water levels for all missions in repetitive orbits whose data contained in the GDR were processed using the OCOG retracking algorithm.For this purpose, comparisons between altimetry-based water stages and in situ measurements from a dense gauges network were performed in the Inner Niger Delta (IND).
Several studies already used altimetry-based water levels to better understand spatio-temporal dynamics of the flood in this region [20][21][22] and to estimate river discharges [23,24].The IND was chosen as study area as it is densely covered with (i) in situ gauge stations whose records are available over the whole high-precision altimetry era and (ii) cross-sections between altimetry ground-tracks from different missions and rivers of various widths.
A dense network of VS was built in the IND, composed of 52, 63, 623 VS for European Remote-Sensing Satellite-2 (ERS-2), ENVIronment SATellite (ENVISAT), Satellite with Argos and ALtika (SARAL) respectively, 31 for Sentinel-3A and 8, 8, 9 for Jason-1, Jason-2 and Jason-3 respectively.At each virtual station, time variations of river levels from radar altimetry are constructed.
In this study, comparisons between altimetry-based water stages derived from acquisitions of Jason-1, Jason-2 and Jason-3, ERS-2, ENVISAT, SARAL and Sentinel-3A and in situ water levels from 19 gauge stations located in the IND are presented in terms of RMSE and R. Intra-mission results consistency were also assessed through cross-correlations between virtual stations along the river.

Principle of Altimetry Measurement
The principle of radar altimetry is the following: a radar altimeter emits an electromagnetic wave in the nadir direction and measures its round-trip time.The distance between the satellite and the Earth surface-the altimeter range (R 0 )-is derived with a precision of a few centimeters.The satellite altitude (H) referred to an ellipsoid is determined from precise orbitography technique with accuracy better than 2 cm.Taking into account propagation corrections caused by delays resulting from interactions of electromagnetic wave with the atmosphere, and geophysical corrections, the height of the reflecting surface (h) with reference to an ellipsoid can be estimated as [25,26]: where H is the height of the center of mass of the satellite above the ellipsoid estimated using precise orbit determination (POD) technique, R 0 is the nadir altimeter range from the center of mass of the satellite to the sea surface taking into account instrumental corrections.
where ∆R ion is the atmospheric refraction range correction due to the free electron content associated with the dielectric properties of the ionosphere, ∆R dry is the atmospheric refraction range correction due to the dry gas component of the troposphere, ∆R wet is the atmospheric refraction range correction due to the water vapor and the cloud liquid water content of the troposphere.
∑ ∆R geophysical = ∆R solid Earth + ∆R pole (3) where ∆R solid Earth and ∆R pole are the corrections respectively accounting for crustal vertical motions due to the solid Earth and pole tides.

Time Variations of River Levels from Radar Altimetry Measurements
In this study, the Multi-mission Altimetry Processing Software (MAPS), developed by Frappart et al. [27] was used to visualize and process the altimetry data over land [22,[28][29][30] and ocean [31,32] to build the VS in the IND.Data processing is composed of three main steps: (i) a coarse delineation of the VS using Google Earth; (ii) a refined selection of the valid altimetry data based on visual inspection; and (iii) the computation of the time series of water level.The altimetry-based water level is computed for each cycle using the median of the selected altimetry heights, along with their respective deviation (i.e., mean absolute deviation).This process is repeated each cycle to construct the water level time series at the virtual stations and illustrated in Figure 1.
Altimetry datasets are referenced either to WGS84 ellipsoid or to Topex/Poseidon ellipsoid.A datum conversion from T/P ellipsoid to WGS84 is automatically performed using Equation (4) adapted from Jekeli et al. [33] and implemented in the version of MAPS used in Salameh et al. [34]: where ∆h is the variation of height at latitude ϕ due to the change of ellipsoid from T/P to WGS84 datum, a = 6, 378, 137 m and e = 0.081819190842621 are the semi-major axis and the eccentricity of the WGS84 datum, a' = 6,378,136.3m and e' = 0.081819221456 are the semi-major axis and the eccentricity of the T/P datum.Some along-track altimetry profiles exhibit a parabolic shape caused by non-nadir reflections known as hooking effect (see Figure 2).Hooking effect is corrected as follows: where s is the along-track coordinate, h(s0) is the altimeter height at nadir, Rcorr(s0) the altimeter range at nadir corrected from the geophysical and environmental effects, s0 the location of the nadir along the altimeter track, si the coordinates of the slant measurements, / the rate of altitude variation of the satellite along the orbital segment, and ds the along track difference between s0 and si.Then, the altimeter height at nadir is computed using the summit of the parabola representing the actual water level: (6) Figure 1.The different steps of the altimetry data using Multi-mission Altimetry Processing Software (MAPS).First, a rough selection of the altimetry data (represented with crosses of color) is performed: (a) all altimetry data located a few kilometers away from the center of the river are selected (left panel), the major topographic features, such hills, river banks, tributaries, etc. can be identified (central panel), temporal variations of the altimetry signal cannot be related to any hydrological signal (right panel).Then, a more accurate selection is made; (b) the number of data is decreasing and outliers are removed (left panel), the shape of the river and its temporal variations in width (central panel) and height (left panel) becomes clearer.This process is repeated until final selection is achieved; (c) all subfigures are derived from the MAPS Graphical User Interface (GUI).Some along-track altimetry profiles exhibit a parabolic shape caused by non-nadir reflections known as hooking effect (see Figure 2).Hooking effect is corrected as follows: where s is the along-track coordinate, h(s 0 ) is the altimeter height at nadir, R corr (s 0 ) the altimeter range at nadir corrected from the geophysical and environmental effects, s 0 the location of the nadir along the altimeter track, s i the coordinates of the slant measurements, ∂H/∂s the rate of altitude variation of the satellite along the orbital segment, and ds the along track difference between s 0 and s i .Then, the altimeter height at nadir is computed using the summit of the parabola representing the actual water level: where a, b and c are parabola coefficients calculating using a least-square fitting of the altimeter data affected by hooking.Finally, s 0 and h(s 0 ) are defined as follow: Remote Sens. 2018, 10, x FOR PEER REVIEW 5 of 28 where a, b and c are parabola coefficients calculating using a least-square fitting of the altimeter data affected by hooking.Finally, s0 and h(s0) are defined as follow:

Validation of the Altimetry-Based Water Levels
Validation of the altimetry-based water levels is performed against records from the closest in situ gauge stations.The along-stream distance between a VS and the closest in situ gauge stations is generally lower than 100 km (see Tables S1-S5).Root Mean Square Error (RMSE), R and R 2 values were estimated between altimetry-based water levels and in situ ones measured the same day using the classical formulas:

Validation of the Altimetry-Based Water Levels
Validation of the altimetry-based water levels is performed against records from the closest in situ gauge stations.The along-stream distance between a VS and the closest in situ gauge stations is generally lower than 100 km (see Tables S1-S5).Root Mean Square Error (RMSE), R and R 2 values were estimated between altimetry-based water levels and in situ ones measured the same day using the classical formulas: where h alti and h in situ are the altimetry-based and the in situ water stages respectively, t i is the measurement time and n the number of common observations.The average of a variable x is written <x>.
As the in situ gauge stations are leveled against a reference unavailable to us, no bias estimates were computed between the in situ and the altimetry-based water levels, but they were between the different missions in the same orbit as follows: while h alti1 is the more recent mission in the orbit and h alti2 is the older one.
The consistency of the intra-mission altimetry-based water levels as well as likely time-lag between water stages in the IND were estimated using the maximum of the cross-correlation function R hh and the argument of the maximum: where τ is the time displacement.

Study Area
The IND is an extensive Sahelian floodplain located between longitudes 3-5 • W and latitudes 13-17 • N in Central Mali (Figure 3a).It is encompassed between the in situ gauge stations of Macina (−5.37 • W, 13.95 • N), on the Niger River, and Douna (−5.9 • W, 13.22 • N), on the Bani River, upstream, and Diré (−3.38 • W, 16.27 • N), downstream.Its drainage area represents a surface of 73,000 km 2 [35].The flooded area extent depends on the intensity of the West African Monsoon and can reach 35,000 km 2 during the wettest rainy seasons [36][37][38][39][40].The flooding period ranges from August to December and during the dry season, from March to May, the area dries out with the exception of the rivers mainstem and the permanent lakes3. 2

. Radar Altimetry Data
The data used in this study come from the acquisitions of the following radar altimetry missions in their nominal orbit: Jason-

Missions with a 10-Day Repeat Period (Jason-1, Jason-2 and Jason-3)
Jason-1, Jason-2 and Jason-3 orbit at an altitude of 1336 km, with an inclination of 66°, on a 10day repeat cycle, providing observations of the Earth surface (ocean and land) from 66° latitude North to 66° latitude South, with an equatorial ground-track spacing of about 315 km.This orbit was formerly used by Topex/Poseidon mission.Jason-1, Jason-2 and Jason-3 orbit at an altitude of 1336 km, with an inclination of 66 • , on a 10-day repeat cycle, providing observations of the Earth surface (ocean and land) from 66 • latitude North to 66 • latitude South, with an equatorial ground-track spacing of about 315 km.This orbit was formerly used by Topex/Poseidon mission.
Jason-1 mission was launched on 7 December 2001 by a cooperation between CNES and National Aeronautics and Space Administration (NASA).Jason-1 sensors are based on the former Topex/Poseidon missions, composed of the Poseidon-2 altimeter which is a two-frequency altimeter with C (5.3 GHz) and Ku (13.575GHz) -bands.Its payload is also composed of the Jason Microwave Radiometer from NASA and a triple system for precise orbit determination: DORIS instrument from the CNES, Black Jack Global Positioning System receiver from NASA and a LRA from NASA/Jet Propulsion Laboratory (JPL) [43].Jason-1 remained in its nominal orbit until 26 January 2009 and was decommissioned on 21 June 2013.
Jason-2 mission was launched on 20 June 2008 as a cooperation between CNES, the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT), NASA and the National Oceanic and Atmospheric Administration (NOAA).Its payload is mostly composed of the Poseidon-3 radar altimeter from CNES, the Advanced Microwave Radiometer (AMR) from JPL/NASA, and a triple system for precise orbit determination: the real-time tracking system DIODE of DORIS instrument from CNES, a Global Navigation Satellite System (GNSS) receiver and a LRA from NASA/JPL.Poseidon-3 radar altimeter is a two-frequency solid-state altimeter that measures accurately the distance between the satellite and the surface (range) and provides ionospheric corrections over the ocean [44].It operates at Ku and C bands.Raw data are processed by SSALTO (Segment Sol multimissions d'ALTimétrie, d'Orbitographie).Jason-2 remained in its nominal orbit until 3 July 2016.
Jason-3 mission was launched on 17 January 2016 as cooperation between CNES, EUMETSAT, NASA and NOAA.This satellite is composed of Poseidon-3B radar altimeter with a Precise Orbit Determination (POD) package with a Global Positioning System (GPS) receiver, DORIS and a LRA from NASA/JPL.

Mission with a 27-Day Repeat Period (Sentinel-3A)
Sentinel-3A mission was launched on 16 February 2016 by ESA to an orbit of altitude 814 km.The satellite caries one altimeter radar called SRAL (SAR Radar ALtimeter), a dual-frequency SAR altimeter (Ku-band at 13.575 GHz and C-band at 5.41 GHz).Its payload comprises also a Microwave Radiometer (MWR) instrument for wet path delay measurements and a triple system for precise orbit determination: a POD including a GPS receiver, a LRA and a DORIS instrument [45].
All this information is summarized in Table 1.The data used in this study are summarized in Table 2. Ranges used to derive altimeter heights and backscattering coefficients are those processed with OCOG/Ice-1/Ice retracking algorithm [46].Previous studies showed that Ice-1-derived altimetry heights are the more suitable for hydrological studies in terms of accuracy of water levels and availability of the data (e.g., [12,47,48]) among the commonly available retracked data present in the GDRs.

In Situ Water Levels
Daily stage records from 19 in situ gauge stations located in the IND were used in this study to validate altimetry-based water levels (see Table 3 for their names, locations and periods of data availability and Figure 3c for their locations).Measurements were acquired at 12:00 a.m.local time.They are made available by the Malian water agency (Direction Nationale de l'Hydraulique-DNH).

Direct Validation of the Altimetry-Based Water Stages
The nominal altimetry ground-tracks from ERS-2, ENVISAT, SARAL (35-day repeat orbit), Sentinel-3A (27-day repeat orbit), Jason-1, Jason-2 and Jason-3 (10-day repeat orbit) missions present a large number of cross-sections with river streams and floodplains in the IND (see Figure 3b).A dense network of virtual stations from different missions was defined in the IND (see Table 4).Virtual station (VS) locations in the IND are presented in Figure 3c.Altimetry-based water levels were compared to water stage records from close in situ gauge.These comparisons were performed for VS located on the rivers and not on the floodplains for distances between the in situ gauge and the VS lower than 100 km.In situ gauge records from 19 stations were used to perform the 89 following comparisons: The complete results of these comparisons (distance between the in situ gauge and the VS, number of data used for comparisons (N), RMSE, R and R 2 ) are presented in Tables S1 to S5 for ERS-2, ENVISAT, SARAL, Jason-1/2/3 and Sentinel-3A missions (in supplementary information).The results of these comparisons are also presented as maps in Figure 4 for the altimetry missions that were launched before 2010 (ERS-2, ENVISAT, Jason-1 and 2) and in Figure 5 for the most recent missions.The number of data used for the comparison is, most of the time, statistically significant, except for 15 comparisons against ENVISAT (less than 20 common observations), 4 against SARAL (less than 15 common observations) as well as the whole comparisons against Sentinel-3A as only 16  Very good agreements were generally found between altimetry-based and in situ water stages for all the missions over a total of 89 comparisons performed.Values of R greater than 0.95 were obtained 41 times (45%), between 0.95 and 0.9, 18 times (20%), and between 0.8 and 0.9, 18 times (20%).Correlation coefficients R lower than 0.7 were obtained only 4 times (4%) (Figure 6a).The minimum R value is 0.57.RMSE lower than 0.3 m were obtained 12 times (13%), between 0.3 and 0.5 m, 17 times (19%), between 0.5 and 0.75 m, 29 times (32%) and above 1 m, 19 times (21%) (Figure 6b).
Better agreement was found for the recent missions such as SARAL, Jason-3 and Sentinel-3A than for the older ones (ERS-2, ENVISAT, Jason-1 and Jason-2).Focusing on the correlations, for instance R was greater than 0.9 in 10 out of 19 comparisons for ERS-2 (52%), in 19 out of 32 (59%) for ENVISAT, in 10 out of 14 (71%) for SARAL, in 14 out of 16 (88%) for Sentinel-3A, 1 (R = 0.89) out of 3 (33%) for Jason-1, 3 out of 3 for Jason-2 (100%) and 3 out of 3 (100%) for Jason-3 (Tables S1-S5).The first results of the use of Sentinel-3A are very encouraging as only very few cycles were available, keeping in mind that among the three correlation coefficients lower than 0.9, three were already higher than 0.85. the older ones (ERS-2, ENVISAT, Jason-1 and Jason-2).Focusing on the correlations, for instance R was greater than 0.9 in 10 out of 19 comparisons for ERS-2 (52%), in 19 out of 32 (59%) for ENVISAT, in 10 out of 14 (71%) for SARAL, in 14 out of 16 (88%) for Sentinel-3A, 1 (R = 0.89) out of 3 (33%) for Jason-1, 3 out of 3 for Jason-2 (100%) and 3 out of 3 (100%) for Jason-3 (Tables S1-S5).The first results of the use of Sentinel-3A are very encouraging as only very few cycles were available, keeping in mind that among the three correlation coefficients lower than 0.9, three were already higher than 0.85.Combining altimetry-based time series of water levels from missions in the same nominal orbits (i.e., ERS-2, ENVISAT and SARAL, and Jason-1, Jason-2 and Jason-3) or at inter-mission cross-overs (ERS-2/ENVISAT/SARAL with Sentinel-3A ground-tracks), multi-mission time series of water levels were obtained.In this latter case, the maximum difference in time between SARAL and Sentinel-3A acquisitions is half the length of the shortest repeat period of the two satellites (i.e., 13 days considering the 27 days of Sentinel-3A repeat period).Inter-mission biases were removed: Low water periods were chosen, rather than high water periods or the complete hydrological cycle, to minimize the effect of the difference in temporal sampling, assuming that water levels are more stable during low water stages.
Examples of multi-mission time series of water levels are presented in Figures 7 and 8 for ERS-2/ENVISAT/SARAL and Sentinel-3A (when there is a cross-over) and for Jason-1, Jason-2 and Jason-3 respectively along with in situ gauge station gauge records of Diré upstream part of the IND, Mopti, central part and Macina, downstream part, over 1995-2017 (Figure 7a-c respectively), and of Macina, over 2002-2017 (Figure 8).
In examples presented in Figure 7, a very good agreement is found between altimetry-based and in situ water stages.There is a gap in the time series between November 2010 and January 2013 as no altimetry mission was in the 35-day repeat orbit during this period.In spite of the distance between the VS and the two first in situ gauge stations (77 and 40 km, with Diré and Mopti stations respectively, but only 1 km from the Macina station under the ERS-2/ENVISAT and SARAL groundtracks, and 8 km from Macina station under Sentinel-3A ground-tracks), better results were found in the first examples than in the latter one, with higher R 2 and lower RMSE for ERS-2 and ENVISAT and similar ones for SARAL and Sentinel-3A.An underestimation of the annual amplitude of the water levels is observed during the ERS-2 observation period.In the example presented in Figure 8, the VS and the Macina in situ gauge station are separated by only 1 km.The quality of the water stage retrieval is increasingly better from Jason-1 to Jason-3.It is important to mention that Jason-1 data contained in the GDR E, released in May 2016, allow the accurate and continuous estimation of water stages over the IND contrary to the previous GDR versions that contained few useful data over land (except over large lakes, see [18,51,52] for instance).Combining altimetry-based time series of water levels from missions in the same nominal orbits (i.e., ERS-2, ENVISAT and SARAL, and Jason-1, Jason-2 and Jason-3) or at inter-mission cross-overs (ERS-2/ENVISAT/SARAL with Sentinel-3A ground-tracks), multi-mission time series of water levels were obtained.In this latter case, the maximum difference in time between SARAL and Sentinel-3A acquisitions is half the length of the shortest repeat period of the two satellites (i.e., 13 days considering the 27 days of Sentinel-3A repeat period).Inter-mission biases were removed: • using the acquisitions made during tandem phases when two missions were in the same orbit a few seconds or minutes apart from each other (e.g., Jason-1 and Jason-2, Jason-2 and Jason-3, ERS-2 and ENVISAT); Low water periods were chosen, rather than high water periods or the complete hydrological cycle, to minimize the effect of the difference in temporal sampling, assuming that water levels are more stable during low water stages.
Examples of multi-mission time series of water levels are presented in Figures 7 and 8 for ERS-2/ENVISAT/SARAL and Sentinel-3A (when there is a cross-over) and for Jason-1, Jason-2 and Jason-3 respectively along with in situ gauge station gauge records of Diré upstream part of the IND, Mopti, central part and Macina, downstream part, over 1995-2017 (Figure 7a-c respectively), and of Macina, over 2002-2017 (Figure 8).
In examples presented in Figure 7, a very good agreement is found between altimetry-based and in situ water stages.There is a gap in the time series between November 2010 and January 2013 as no altimetry mission was in the 35-day repeat orbit during this period.In spite of the distance between the VS and the two first in situ gauge stations (77 and 40 km, with Diré and Mopti stations respectively, but only 1 km from the Macina station under the ERS-2/ENVISAT and SARAL ground-tracks, and 8 km from Macina station under Sentinel-3A ground-tracks), better results were found in the first examples than in the latter one, with higher R 2 and lower RMSE for ERS-2 and ENVISAT and similar ones for SARAL and Sentinel-3A.An underestimation of the annual amplitude of the water levels is observed during the ERS-2 observation period.In the example presented in Figure 8, the VS and the Macina in situ gauge station are separated by only 1 km.The quality of the water stage retrieval is increasingly better from Jason-1 to Jason-3.It is important to mention that Jason-1 data contained in the GDR E, released in May 2016, allow the accurate and continuous estimation of water stages over the IND contrary to the previous GDR versions that contained few useful data over land (except over large lakes, see [18,51,52] for instance).

Intermission Water Stage Comparison
Comparisons between water levels estimated by altimetry missions during their tandem phase (i.e., few cycles during which two missions were orbiting a few minutes apart from one another) were performed between ERS-2 and ENVISAT from June 2002 to July 2003 (11 cycles in common), Jason-1 and Jason-2 from July 2008 to January 2009 (21 cycles in common), Jason-2 and Jason-3 from January to September 2016 (23 cycles in common).They allow the increase of the number of comparisons in the IND, not only on the rivers but also on the wetlands that are monitored using in situ gauges.

Intermission Water Stage Comparison
Comparisons between water levels estimated by altimetry missions during their tandem phase (i.e., few cycles during which two missions were orbiting a few minutes apart from one another) were performed between ERS-2 and ENVISAT from June 2002 to July 2003 (11 cycles in common), Jason-1 and Jason-2 from July 2008 to January 2009 (21 cycles in common), Jason-2 and Jason-3 from January to September 2016 (23 cycles in common).They allow the increase of the number of comparisons in the IND, not only on the rivers but also on the wetlands that are monitored using in situ gauges.
Better agreement is found between altimetry-based water stages from Jason-2 and 3. Correlation coefficients R greater than 0.95 were obtained in all the cases except one (86%) for which R equals to 0.84 (Figure 12b).Values of RMSE lower than 0.30 m were obtained in all the cases except one (86%) for which RMSE equals to 0.40 m (Figure 12c).Biases between Jason-1 and 2 are very similar, ranging from −0.36 to −0.14 m (−0.27 ± 0.08 m) (Figure 12d).A total of 5 and 7 comparisons between Jason-1 and Jason-2, and Jason-2 and Jason-3-based water levels were performed in the IND with several samples (N) varying from 5 to 15 for Jason-1/Jason-2, and from 13 to 22 for Jason-2/Jason-3.Three of these comparisons were performed on more than 10 samples (60%) for Jason-1 (Figure 11a).Very good agreement was generally found between altimetry-based water stages from Jason-1 and 2. Values of R greater than 0.95 were obtained 3 times (60%).In the two other cases, R equal to 0.80 and 0.69 were found (Figure 11b).Values of RMSE lower than 0.3 m were obtained 3 times (60%).In the two other cases, RMSE lower than 0.35 m were found (Figure 11c).Biases between Jason-1 and 2 ranged between 0.34 and 1.06 m (0.75 ± 0.28 m) (Figure 11d).
Better agreement is found between altimetry-based water stages from Jason-2 and 3. Correlation coefficients R greater than 0.95 were obtained in all the cases except one (86%) for which R equals to 0.84 (Figure 12b).Values of RMSE lower than 0.30 m were obtained in all the cases except one (86%) for which RMSE equals to 0.40 m (Figure 12c).Biases between Jason-1 and 2 are very similar, ranging from −0.36 to −0.14 m (−0.27 ± 0.08 m) (Figure 12d).A total of 5 and 7 comparisons between Jason-1 and Jason-2, and Jason-2 and Jason-3-based water levels were performed in the IND with several samples (N) varying from 5 to 15 for Jason-1/Jason-2, and from 13 to 22 for Jason-2/Jason-3.Three of these comparisons were performed on more than 10 samples (60%) for Jason-1 (Figure 11a).Very good agreement was generally found between altimetry-based water stages from Jason-1 and 2. Values of R greater than 0.95 were obtained 3 times (60%).In the two other cases, R equal to 0.80 and 0.69 were found (Figure 11b).Values of RMSE lower than 0.3 m were obtained 3 times (60%).In the two other cases, RMSE lower than 0.35 m were found (Figure 11c).Biases between Jason-1 and 2 ranged between 0.34 and 1.06 m (0.75 ± 0.28 m) (Figure 11d).
Better agreement is found between altimetry-based water stages from Jason-2 and 3. Correlation coefficients R greater than 0.95 were obtained in all the cases except one (86%) for which R equals to 0.84 (Figure 12b).Values of RMSE lower than 0.30 m were obtained in all the cases except one (86%) for which RMSE equals to 0.40 m (Figure 12c).Biases between Jason-1 and 2 are very similar, ranging from −0.36 to −0.14 m (−0.27 ± 0.08 m) (Figure 12d).Examples of time series combining ERS-2 and ENVISAT, Jason-1 and 2, and Jason-2 and 3 are presented on Figure 13a-c.Virtual stations from close locations were chosen: ERS-2/ENVISAT VS 0545-a and Jason-1/Jason-2/Jason-3 VS 046-d.They are located in the upstream part of the IND at an approximate along-stream distance of 45 km.Seasonal amplitudes of between 4 and 5 m are observed during the common period of availability of the different altimetry-based water levels.The agreement is better for the recent missions than for the older ones: R increases from 0.92 to 0.99 whereas RMSE decreases from 0.62 to 0.22 m as well as the bias from −1.35 to 0.36 m.Lower deviations are generally observed on the time series from the more recent missions than other the older ones: 0.13 ± 0.11, 0.13 ± 0.15, 0.19 ± 0.22, 0.21 ± 0.16, 0.15 ± 0.11 m were obtained averaging the mean absolute deviation from individual cycles for ERS-2, ENVISAT, Jason-1, 2 and 3 respectively.Please note that, for readability purpose, biases between the time series were removed in Figure 13.

Multi-Mission Time Series on Floodplains
Floodplains and wetlands are generally not monitored using in situ gauges.Radar altimetry is a unique tool for the long-term observations of the changes in water levels over inundated areas [53][54][55][56][57].Among the VS defined in the IND, 16 were built in floodplains under ERS-2/ENVISAT/SARAL ground-tracks and 9 under Sentinel-3A ground-tracks, but none under Jason-1, Jason-2 and Jason-3 ground-tracks.Their consistency was checked during the tandem phase between ERS-2 and ENVISAT (Figure 14).Values of R greater than 0.95 were obtained 10 times (63%), between 0.95 and 0.9, twice (12%), and between 0.8 and 0.9, 3 times (19%).However, R lower than 0.5 was obtained once (6%) (Figure 14a).Values of RMSE lower than 0.3 m were obtained 7 times (44%), between 0.3 and 0.5 m, 5 times (31%), between 0.5 and 0.75 m, 3 times (19%) and above 1 m, once (6%) (Figure 14b).The agreement is better for the recent missions than for the older ones: R increases from 0.92 to 0.99 whereas RMSE decreases from 0.62 to 0.22 m as well as the bias from −1.35 to 0.36 m.Lower deviations are generally observed on the time series from the more recent missions than other the older ones: 0.13 ± 0.11, 0.13 ± 0.15, 0.19 ± 0.22, 0.21 ± 0.16, 0.15 ± 0.11 m were obtained averaging the mean absolute deviation from individual cycles for ERS-2, ENVISAT, Jason-1, 2 and 3 respectively.Please note that, for readability purpose, biases between the time series were removed in Figure 13.

Consistency of the Altimetry-Based Water Levels in the Inner Niger Delta (IND)
A consistency check was performed between a reference VS located either in Niger River mainstem or in the Bani major tributary and the other VS from the same altimetry mission located upstream and downstream on the same river course estimating the maximum of correlation using the cross-correlation function and the associated time-lag, similarly as in Bogning et al. [29].Crosscorrelation function maxima and associated time-lags are presented for ERS-2 and ENVISAT in Figures 16 and S1

Consistency of the Altimetry-Based Water Levels in the Inner Niger Delta (IND)
A consistency check was performed between a reference VS located either in Niger River mainstem or in the Bani major tributary and the other VS from the same altimetry mission located upstream and downstream on the same river course estimating the maximum of correlation using the cross-correlation function and the associated time-lag, similarly as in Bogning et al. [29].Cross-correlation function maxima and associated time-lags are presented for ERS-2 and ENVISAT in Figure 16 and Figure S1 respectively and for SARAL and Sentinel-3A in Figure 17  Due to the changes of the river features (slope, depth, width, etc.), biases and RMSEs were not computed between the time series of altimetry-based water levels.Correlation coefficients (R) higher than 0.7 are generally observed for ERS-2 and ENVISAT, except in the upstream part of the delta for ERS-2, with correlation ranging from 0.5 to 0.7 (Figure 16c).Correlations higher than 0.8 are observed in the IND downstream and central parts for ENVISAT (Figure 16).Very high correlations, generally higher than 0.85, especially in the downstream part, except on a few locations, were observed for SARAL and, especially, for Sentinel-3A (Figure 17).No time-lag was observed over the Bani River for none of the altimetry missions (Figures S1a, d and S2a,d).Time lags of plus or minus one cycle were respectively found for the upper and lower stations along the Niger mainstem (Figures S1 and S2).Correlation coefficients (R) higher than 0.7 are generally observed for ERS-2 and ENVISAT, except in the upstream part of the delta for ERS-2, with correlation ranging from 0.5 to 0.7 (Figure 16c).Correlations higher than 0.8 are observed in the IND downstream and central parts for ENVISAT (Figure 16).Very high correlations, generally higher than 0.85, especially in the downstream part, except on a few locations, were observed for SARAL and, especially, for Sentinel-3A (Figure 17).No time-lag was observed over the Bani River for none of the altimetry missions (Figures S1a, d and S2a,d).Time lags of plus or minus one cycle were respectively found for the upper and lower stations along the Niger mainstem (Figures S1 and S2).

Discussion
The coverage of altimetry tracks from missions of the 35-day (ERS-2/ENVISAT/SARAL) and the 27-day (Sentinel-3A) repeat orbits allows the construction of a dense network of VS on both rivers and floodplains in the IND, completed by a few VS from missions in the 10-day repeat orbit (Figure 3c).Yet, a larger number (more than 15%) of ENVISAT (6) and SARAL-based (62) VS than ERS-2-based (52) ones were built on the ground-tracks of the 35-day repeat orbit.These VS with no valid ERS-2-based water levels are located on the upstream part of the IND for 5 of them, central part for 3 others and downstream part for the last 3 ones.
Due to the relative flatness of the IND, the no-construction of VS cannot be attributed to data losses caused by tracking issues, but to incorrect range estimates.They were likely to be caused by the narrow width of the river streams, especially in the upstream part, where they are generally lower than 300 m.They can also be accounted for the possible complexity of the waveforms (e.g., multi-peaked over areas with several bright targets encompassed in the altimeter footprint such as several river streams that can reach 1 km of width as in the downstream part of the IND, or several floodplain lakes as in the central part).
Comparisons with a high number of in situ gauge records show an overall very good agreement between in situ and altimetry-based water levels, increasingly better for the most recent missions than the older ones.Contrary to what was found in other study areas with variations of topography (e.g., [28,29]), the data used to build the VS were mostly acquired at 330 and 320 MHz Ku chirp bandwidth operation mode by both ERS-2 and ENVISAT (i.e., with the better range resolution in more than 95% of the cases for any VS-Figure 18).This is the reason quite similar results when comparing ERS-2 and ENVISAT-based water stages with in situ water levels were found in terms of R, with generally higher RMSE for ERS-2 (see Figure 4, Tables S1 and S2), and also a good agreement when comparing ERS-2 and ENVISAT-based water levels during their tandem phase, except in terms of bias (Figure 10).In the same orbit, better results were found using SARAL data (Figure 5 and Table S3), as already observed in other river basins, resulting from the use of the Ka-band with its smaller footprint and higher chirp frequency of 500 MHz [29,58].
Very similar numbers of VS were built using Jason-1, Jason-2 and Jason-3 data.The release of the Jason-1 GDR E, which contained valid data over land and not mostly only over large lakes as in the previous versions, allowed the extension of the duration of the time series of water levels over rivers for the missions in the 10-day repeat orbit down to 2001, at the expense of lower accuracy than using Jason-2 and Jason-3 data (Figures 4 and 5 and Table S4).This lower accuracy can be accounted for the small river widths, most of the time lower than 500 m, except in some locations in the downstream part of the IND.
On the contrary, Jason-2 and Jason-3 confirm their great capability for detecting and accurately estimating water levels over narrow rivers, especially over flat areas [28,30].Comparisons performed during the tandem phase logically confirm the lower agreement between Jason-1 and Jason-2 than between Jason-2 and Jason-3 (Figures 11 and 12).In terms of bias, if high and variable biases were found between Jason-1 and Jason-2-based time series of water levels, much lower ones with low variability were found between Jason-2 and Jason-3 ones during their tandem period.Poseidon-3, on-board Jason-3, is nominally operating in open loop or DIODE/DEM tracking mode over land surfaces, meaning that the reception window is controlled by an a-priori elevation from an along-track DEM loaded in the altimeter [30].
Very good results were also found using one year and a half of measurements from Sentinel-3A, the first altimeter to operate in SAR over all types of surfaces.During this short time period, the results obtained are almost as good as the ones obtained using SARAL, confirming the strong potential of this technique for monitoring inland water stages.This is likely due to the sharper waveform obtained in SAR mode compared with Low Resolution Mode (LRM) [59] as OCOG is particularly well adapted for very specular echoes from a single reflector (e.g., a river stream).On the contrary, less accurate estimates can be expected from multi-peak waveforms.The consistency of intra-mission water levels derived from radar altimetry measurements was analyzed using cross-correlations for ERS-2, ENVISAT, SARAL and Sentinel-3A.In this case, too, higher correlations were found for the most recent missions operating in (LRM) at Ka band (SARAL) and in SAR at Ku band (Sentinel-3A).Time-lags corresponding to the maximum of cross-correlation (Figures S1 and S2) also showed that there is a time-lag of around two months between the upstream and the downstream part of the IND, in accordance with the water residence time estimated using satellite images [39].The consistency of intra-mission water levels derived from radar altimetry measurements was analyzed using cross-correlations for ERS-2, ENVISAT, SARAL and Sentinel-3A.In this case, too, higher correlations were found for the most recent missions operating in (LRM) at Ka band (SARAL) and in SAR at Ku band (Sentinel-3A).Time-lags corresponding to the maximum of cross-correlation (Figures S1 and S2) also showed that there is a time-lag of around two months between the upstream and the downstream part of the IND, in accordance with the water residence time estimated using satellite images [39].

Conclusions
An extensive assessment of the performance of almost all the missions put in orbit from ERS-2 to Sentinel-3A was performed over the IND.Around 90 comparisons between in situ and altimetry-based were performed as well as more than 75 intermission comparisons between ERS-2 and ENVISAT, Jason-1 and Jason-2, and Jason-2 and Jason-3.Results of these comparisons show (i) a better agreement between altimetry-derived and in situ water stages over rivers of width varying from a few hundreds of meters to ~1.5 km and floodplains, (ii) an increase in accuracy in water level estimates for the most recent missions such as SARAL, Jason-3 and Sentinel-3A benefiting from the use of the Ka-band for SARAL and of the SAR acquisition mode for Sentinel-3A.
Due to the small number of VS defined on Jason-3 ground-tracks, the DIODE/DEM was not evaluated.Intra-mission consistency check performed on ERS-2, ENVISAT, SARAL and Sentinel-3A missions exhibits a time-lag for the maximum of around 2 months that can be related to the residence time of the water in the IND.
Very good performance of the recently launched altimeter SARAL onboard Sentinel-3A ensures the continuity of the monitoring of the IND whose data density should increase with the launch, in 2018, of Sentinel-3B, in the same orbit as Sentinel-3A but Sentinel-3B flies at ±140 • out of phase with Sentinel-3A.
This network of VS will present a strong interest for (i) assessing the impacts of climate variability and human effects (e.g., dam operation and rice production), and (ii) validating the measurement of the future Surface Water and Ocean Topography (SWOT) mission, the first mission to operate in close nadir interferometry SAR for providing elevations in two swaths, to be launched in 2021, on both rivers and floodplains.S1: Results of the comparisons between in situ and altimetry water levels derived from ERS-2 mission: number of sample (N), correlation coefficient (R), determination coefficient (R 2 ) and Root Mean Square Error (RMSE), Table S2: Results of the comparisons between in situ and altimetry water levels derived from ENVISAT mission: number of sample (N), correlation coefficient (R), determination coefficient (R 2 ) and Root Mean Square Error (RMSE), Table S3: Results of the comparisons between in situ and altimetry water levels derived from SARAL mission: number of sample (N), correlation coefficient (R), determination coefficient (R 2 ) and Root Mean Square Error (RMSE), Table S4: Results of the comparisons between in situ and altimetry water levels derived from Jason-1, Jason-2 and Jason-3 missions: number of sample (N), correlation coefficient (R), determination coefficient (R 2 ) and Root Mean Square Error (RMSE), Table S5: Results of the comparisons between in situ and altimetry water levels derived from Sentinel-3A mission: number of sample (N correlation coefficient (R), determination coefficient (R 2 ) and Root Mean Square Error (RMSE).
Author Contributions: C.N., F.F., A.T.D. and E.M. designed the study.C.N. processed the altimetry data to produce the time series of water levels.A.T.D. and N.B. provided the in situ data and gave important information on their quality.All of the authors of the present work contributed to the discussion of the results, as well as the writing of the manuscript.
Funding: This study was supported by Centre National d'Etudes Spatiales (CNES) through "Centre de Topographie des Océans et de l'Hydrosphère-CTOH" grant.Cassandra Normandin is supported by Ministère de l'Enseignement Supérieur et de la Recherche through a PhD grant.

Figure 1 .
Figure 1.The different steps of the altimetry data using Multi-mission Altimetry Processing Software (MAPS).First, a rough selection of the altimetry data (represented with crosses of color) is performed:(a) all altimetry data located a few kilometers away from the center of the river are selected (left panel), the major topographic features, such hills, river banks, tributaries, etc. can be identified (central panel), temporal variations of the altimetry signal cannot be related to any hydrological signal (right panel).Then, a more accurate selection is made; (b) the number of data is decreasing and outliers are removed (left panel), the shape of the river and its temporal variations in width (central panel) and height (left panel) becomes clearer.This process is repeated until final selection is achieved; (c) all subfigures are derived from the MAPS Graphical User Interface (GUI).

Figure 2 .
Figure 2. Example of the hooking correction.(a) On a rough selection, parabolic profiles in different cycles were identified.(b) Zooming on them, it appears that they are responsible for deviation of several tenths of centimeters of the river levels.(c)Once the correction of the hooking effect is applied, the deviation is reduced to a maximum of a couple of tenths cm).

Figure 2 .
Figure 2. Example of the hooking correction.(a) On a rough selection, parabolic profiles in different cycles were identified.(b) Zooming on them, it appears that they are responsible for deviation of several tenths of centimeters of the river levels.(c)Once the correction of the hooking effect is applied, the deviation is reduced to a maximum of a couple of tenths cm).
developed for measuring height over ocean, land and ice caps.ERS-2 data are available from 17 May 1995 to 9 August 2010.After 22 June 2003, the dataset coverage is limited to ground station visibility.ENVISAT mission was launched on 1 March 2002 by ESA.It carried 10 instruments including the advanced radar altimeter (RA-2).It was based on the heritage of the sensor on-board the ERS-1 and 2 satellites.RA-2 was a nadir-looking pulse-limited radar altimeter operating at two frequencies at Ku-(13.575GHz), as ERS-1 and 2, and S-(3.2GHz) bands [41].ENVISAT remained in its nominal orbit until October 2010 and its mission ended on 8 April 2012.RA-2 stopped operating correctly at S-band in January 2008.Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 28 orbit until October 2010 and its mission ended on 8 April 2012.RA-2 stopped operating correctly at S-band in January 2008.

Figure 3 .
Figure 3. (a) Location of the IND in Africa; (b) Altimetry tracks over the IND from European Remote-Sensing Satellite-2 (ERS-2), ENVIronment SATellite (ENVISAT), Satellite with Argos and ALtika (SARAL) (blue dots), Jasons-1/Jason-2/Jason-3 (yellow dots) and Sentinel-3A (red dots); and (c) Location of virtual stations to calculate water levels (colored circles) using altimetry data and in situ gauge stations (black diamonds) in IND.SARAL mission was launched on 25 February 2013.SARAL is a new collaboration between Centre National d'Etudes Spatiales (CNES) and Indian Space Research Organization (ISRO).Its payload comprises the AltiKa radar altimeter and bi-frequency radiometer, and a triple system for precise orbit determination: the real-time tracking system Détermination Immédiate d'Orbite par Doris embarqué (DIODE) of the Doppler Orbitography and Radio-positioning Integrated by Satellite (DORIS) instrument, a Laser Retroflector Array (LRA), and the Advance Research and Global Observation Satellite (ARGOS-3).AltiKa radar altimeter is a solid-state mono-frequency altimeter that provides accurate range measurements.It is the first altimeter to operate in the Ka-band (35.75 GHz) [42].It has been put in a drifting orbit since July 2016.

Figure 4 .
Figure 4. Comparisons between in situ and altimetry water levels for several missions (a) ERS-2 and Jason-1 data and (b) ENVISAT and Jason-2.For each comparison, the number of samples (N), correlation (R), and RMSE is presented.Diamonds points correspond to Jason data and circles with black contours to ERS-2/ENVISAT.

Figure 4 .
Figure 4. Comparisons between in situ and altimetry water levels for several missions (a) ERS-2 and Jason-1 data and (b) ENVISAT and Jason-2.For each comparison, the number of samples (N), correlation (R), and RMSE is presented.Diamonds points correspond to Jason data and circles with black contours to ERS-2/ENVISAT.

Figure 5 .
Figure 5. Comparisons between in situ and altimetry water levels for several missions (a) SARAL and Jason-3 data and (b) Sentinel-3A.For each comparison, the number of samples (N), correlation (R), and RMSE is presented.Diamonds points correspond to Jason data and circles with black contours to SARAL/Sentinel-3A.

Figure 5 .
Figure 5. Comparisons between in situ and altimetry water levels for several missions (a) SARAL and Jason-3 data and (b) Sentinel-3A.For each comparison, the number of samples (N), correlation (R), and RMSE is presented.Diamonds points correspond to Jason data and circles with black contours to SARAL/Sentinel-3A.

Figure 6 .
Figure 6.Histograms of (a) correlations and (b) RMSE between altimetry-based and in situ water stages for all missions.


using the acquisitions made during tandem phases when two missions were in the same orbit a few seconds or minutes apart from each other (e.g., Jason-1 and Jason-2, Jason-2 and Jason-3, ERS-2 and ENVISAT);  averaging the acquisitions made during the common period of observations at low water stages (April-May-June) for Sentinel-3A and SARAL;  averaging the acquisitions made during low water periods (April-May-June) on different years for ENVISAT (2003-2010) and SARAL (2013-2016).

Figure 6 .
Figure 6.Histograms of (a) correlations and (b) RMSE between altimetry-based and in situ water stages for all missions.
• averaging the acquisitions made during the common period of observations at low water stages (April-May-June) for Sentinel-3A and SARAL; • averaging the acquisitions made during low water periods (April-May-June) on different years for ENVISAT (2003-2010) and SARAL (2013-2016).

Figure 8 .
Figure 8. Altimetry-based water levels from 2002 to 2017 from Jason-1/Jason-2/Jason-3 (red/black/blue respectively) and in situ water levels (light green) data (left).Scatter plots of water levels anomalies from radar altimetry and in situ gauge stations at Macina (right figure, downstream IND).

Figure 8 .
Figure 8. Altimetry-based water levels from 2002 to 2017 from Jason-1/Jason-2/Jason-3 (red/black/blue respectively) and in situ water levels (light green) data (left).Scatter plots of water levels anomalies from radar altimetry and in situ gauge stations at Macina (right figure, downstream IND).

Figure 8 .
Figure 8. Altimetry-based water levels from 2002 to 2017 from Jason-1/Jason-2/Jason-3 (red/black/blue respectively) and in situ water levels (light green) data (left).Scatter plots of water levels anomalies from radar altimetry and in situ gauge stations at Macina (right figure, downstream IND).

Figure 13 .
Figure 13.Intermission quality time series: (a) for the VS 0545-a located in the upstream part in the IND between ERS-2/ENVISAT; (b) for the VS 046-d (upstream part) between Jason-1/Jason-2; and (c) for the VS 046-d between Jason-2/Jason-3.The two VS are separated by 45 km.

Figure 14 .
Figure 14.Histograms of (a) R and (b) RMSE for ERS-2/ENVISAT intermission water stages comparisons on floodplains virtual stations.

Figure 13 .
Figure 13.Intermission quality time series: (a) for the VS 0545-a located in the upstream part in the IND between ERS-2/ENVISAT; (b) for the VS 046-d (upstream part) between Jason-1/Jason-2; and (c) for the VS 046-d between Jason-2/Jason-3.The two VS are separated by 45 km.

Figure 13 .
Figure 13.Intermission quality time series: (a) for the VS 0545-a located in the upstream part in the IND between ERS-2/ENVISAT; (b) for the VS 046-d (upstream part) between Jason-1/Jason-2; and (c) for the VS 046-d between Jason-2/Jason-3.The two VS are separated by 45 km.

Figure 14 .
Figure 14.Histograms of (a) R and (b) RMSE for ERS-2/ENVISAT intermission water stages comparisons on floodplains virtual stations.Figure 14.Histograms of (a) R and (b) RMSE for ERS-2/ENVISAT intermission water stages comparisons on floodplains virtual stations.

Figure 14 .
Figure 14.Histograms of (a) R and (b) RMSE for ERS-2/ENVISAT intermission water stages comparisons on floodplains virtual stations.Figure 14.Histograms of (a) R and (b) RMSE for ERS-2/ENVISAT intermission water stages comparisons on floodplains virtual stations.

Figure 15 .
Figure 15.Water levels time series on floodplains derived from radar altimetry multi-mission (ERS-2/ENVISAT and SARAL) from 1995 to 2016 at 3 different locations in the Inner Niger Delta: the center of the delta (a); the north eastern (b) and the southern (c) parts.
respectively and for SARAL and Sentinel-3A in Figures 17 and S2 respectively for three different VS locations in the IND: near Douna, over the Bani tributary, near Mopti, in the downstream part of the IND, and near Diré in the upstream part of the IND.As the repeat period of these missions is either 35 or 27 days in their nominal orbit, only time-lags of plus or minus one repeat period and a half (i.e., 53 and 41 days) were considered due to the relatively small scale of the IND.

Figure 15 .
Figure 15.Water levels time series on floodplains derived from radar altimetry multi-mission (ERS-2/ENVISAT and SARAL) from 1995 to 2016 at 3 different locations in the Inner Niger Delta: the center of the delta (a); the north eastern (b) and the southern (c) parts.
and Figure S2 respectively for three different VS locations in the IND: near Douna, over the Bani tributary, near Mopti, in the downstream part of the IND, and near Diré in the upstream part of the IND.As the repeat period of these missions is either 35 or 27 days in their nominal orbit, only time-lags of plus or minus one repeat period and a half (i.e., 53 and 41 days) were considered due to the relatively small scale of the IND.Due to the changes of the river features (slope, depth, width, etc.), biases and RMSEs were not computed between the time series of altimetry-based water levels.Remote Sens. 2018, 10, x FOR PEER REVIEW 20 of 28

Figure 16 .
Figure 16.Maxima of cross-correlation between time series of water levels derived from ERS-2 (left) and ENVISAT (right) data in the IND for three different VS near Douna (a,d); Mopti (b,e) and Diré (c,f).Red crosses symbolize the VS chosen as reference (auto-correlation of 1).

Figure 16 .
Figure 16.Maxima of cross-correlation between time series of water levels derived from ERS-2 (left) and ENVISAT (right) data in the IND for three different VS near Douna (a,d); Mopti (b,e) and Diré (c,f).Red crosses symbolize the VS chosen as reference (auto-correlation of 1).

Figure 17 .
Figure 17.Maxima of cross-correlation between time series of water levels derived from SARAL (left) and Sentinel-3A (right) data in the IND for three different VS near Douna (a,d); Mopti (b,e) and Diré (c,f).Red crosses symbolize the VS chosen as reference (auto-correlation of 1).

Figure 17 .
Figure 17.Maxima of cross-correlation between time series of water levels derived from SARAL (left) and Sentinel-3A (right) data in the IND for three different VS near Douna (a,d); Mopti (b,e) and Diré (c,f).Red crosses symbolize the VS chosen as reference (auto-correlation of 1).

Figure 18 .
Figure 18.Percentage of acquisition in the different Ku chirp bandwidth operation mode for each VS in the IND for ENVISAT at 320 MHz (a); at 80 MHz; (b) and at 20 MHz; (c) ERS-2 at 330 MHz (d) and at 82.5 MHz (e).

Figure 18 .
Figure 18.Percentage of acquisition in the different Ku chirp bandwidth operation mode for each VS in the IND for ENVISAT at 320 MHz (a); at 80 MHz; (b) and at 20 MHz; (c) ERS-2 at 330 MHz (d) and at 82.5 MHz (e).

Supplementary Materials:
The following are available online at http://www.mdpi.com/2072-4292/10/6/833/s1,FigureS1: Time-lag corresponding to maxima of cross-correlation between time series of water levels derived from SARAL (left) and SENTINEL-3A (right) data in the IND for three different VS near Douna (a and d), Mopti (b and e) and Diré (c and f).Light blue crosses symbolize the VS chosen as reference (auto-correlation of 1 and time-lag null), Figure S2: Time-lag corresponding to maxima of cross-correlation between time series of water levels derived from SARAL (left) and SENTINEL-3A (right) data in the IND for three different VS near Douna (a and d), Mopti (b and e) and Diré (c and f).Light blue crosses symbolize the VS chosen as reference (auto-correlation of 1 and time-lag null), Table

Table 1 .
Major characteristics of the high-precision radar altimetry missions used in this study.

Table 2 .
Major characteristics of the high-precision radar altimetry missions used in this study.

Table 3 .
List of in situ gauge stations in the IND used this study.

Table 4 .
Number of virtual stations defined in the IND for each mission.