Atmospheric Corrections for Altimetry Studies over Inland Water

Originally designed for applications over the ocean, satellite altimetry has been proven to be a useful tool for hydrologic studies. Altimeter products, mainly conceived for oceanographic studies, often fail to provide atmospheric corrections suitable for inland water studies. The focus of this paper is the analysis of the main issues related with the atmospheric corrections that need to be applied to the altimeter range to get precise water level heights. Using the corrections provided on the Radar Altimeter Database System, the main errors present in the dry and wet tropospheric corrections and in the ionospheric correction of the various satellites are reported. It has been shown that the model-based tropospheric corrections are not modeled properly and in a consistent way in the various altimetric products. While over the ocean, the dry tropospheric correction (DTC) is one of the most precise range corrections, in some of the present altimeter products, it is the correction with the largest errors over continental water regions, causing large biases of several decimeters, and along-track interpolation errors up to several centimeters, both with small temporal variations. The wet tropospheric correction (WTC) from the on-board microwave radiometers is hampered by the contamination on the radiometer measurements of the surrounding lands, making it usable only in the central parts of large lakes. In addition, the WTC from atmospheric models may also have large errors when it is OPEN ACCESS Remote Sens. 2014, 6 4953 provided at sea level instead of surface height. These errors cannot be corrected by the user, since no accurate expression exists for the height variation of the WTC. Alternative and accurate corrections can be computed from in situ data, e.g., DTC from surface pressure at barometric stations and WTC from Global Navigation Satellite System permanent stations. The latter approach is particularly favorable for small lakes and reservoirs, where GNSS-derived WTC at a single location can be representative of the whole lake. For non-timely critical studies, for consistency and stability, model-derived tropospheric corrections from European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis ERA Interim, properly computed at surface height, are recommended. The instrument-based dual-frequency ionospheric correction may have errors related with the land contamination in the Ku and C/S bands, making it more suitable to use a model-based correction. The most suitable model-based ionospheric correction is the Jet Propulsion Laboratory (JPL) global ionosphere map (GIM) model, available after 1998, properly scaled to the altimeter height. Most altimeter products provide the GIM correction unreduced for the total electron content extending above the altitude of these satellites, thus overestimating the ionospheric correction by about 8%. Prior to 1998, the NIC09 (NOAA Ionosphere Climatology 2009) climatology provides the best accuracy.


Introduction
Although originally designed for ocean and ice studies, satellite radar altimetry has been successfully used in the monitoring of continental water surfaces.In recent years, numerous studies have been devoted to different fields of continental hydrology based on satellite altimetry (e.g., [1][2][3][4]).These studies have proven that radar altimetry is currently an essential technique for various applications over inland water, such as: altimeter calibration in the continental domain (e.g., [5][6][7]), study of the hydrological water balance [2,8], assessment of lake-level variation [9,10], studies of anthropogenic impact on lake water storage [11] and climate impacts of lake level fluctuations at the regional scale [12].
The sea or lake surface height, , above a reference ellipsoid, is obtained as: (1) where is the computed satellite height above a reference ellipsoid provided by a precise orbit solution, referred to an International Terrestrial Reference Frame (ITRF), is the observed altimeter range corrected for all instrument effects and is the altimeter range corrected for all instrument, range and geophysical effects.The term includes all range and geophysical corrections and is given by Equation ( 1): (2) The first four terms in Equation ( 2) are the range corrections required to account for the interaction of the radar signal with the atmosphere and with the sea surface.The fifth and sixth terms in Equation ( 2) refer to geophysical phenomena, which must be accounted for in order to separate them from the signals of interest.Using the terminology adopted in the Radar Altimeter Database System (RADS), the last term in Equation ( 2) is called the reference frame offset, , and is only required when multi-mission data are used.This term accounts for all instrument and any systematic effects due to, for example, the adoption of different gravity field models or ITRF in the orbit computation.When applied, the resulting multi-mission time series should be inter-calibrated and continuous.A full description of the range and geophysical corrections can be found, e.g., in [13,14].
In this paper, all corrections are defined as they are usually provided in the altimeter products, that is, the values that must be added to the measured range and therefore subtracted from the sea or lake surface height to get the corresponding corrected values.
Equation (1) shows that in parallel with the improvement of the range estimates through dedicated retracking, a correct modeling and computation of the range and geophysical corrections are of major importance for hydrological studies.For these studies using satellite altimetry, the most relevant range and geophysical corrections are those concerning the effects of the dry troposphere ( ), wet troposphere ( ), the ionosphere ( ) and tides ( ), solid earth, ocean tide loading, pole tide and lake tide, if applicable).The inverse barometer or dynamic atmospheric correction ( ) is not applied, because the lakes/reservoirs are closed systems, and the sea state bias ( ) is usually also not applied, because wind effects tend to be averaged out along-track [15].
Primarily aimed for ocean studies, satellite altimetry products often fail to provide valid corrections over inland water regions.Amongst the required corrections, the effect of the atmosphere in the measured altimeter range is of particular relevance, because of the specificity of these corrections over continental water regions.The atmosphere reduces the speed of the radar pulse, bending its trajectory and, therefore, causing a "path delay" of the altimeter signal, which would, if unaccounted for, result in a too long observed range and, therefore, too short an elevation measurement with respect to a reference ellipsoid.This effect of the atmospheric refraction is usually modeled into three different parts: the delay due to the dry gasses in the troposphere (mainly nitrogen and oxygen), the delay due to the water vapor in the troposphere and the delay due to the free electrons in the upper atmosphere.
The focus of this study is the analysis of the major issues associated with the atmospheric corrections that need to be applied to satellite altimeter range measurements over inland water bodies, often requiring range improvement through retracking, to get centimeter-level measurements of the water surface height.This includes a survey of the main errors present in the current standard altimeter products and the identification of the correct procedures to improve data accuracy.Although many authors have reported various problems in each of these corrections, e.g., [5][6][7]15], particularly on the dry and wet tropospheric corrections, this study shows that the agencies did not take these reports into account, and the various altimeter products still fail to provide atmospheric corrections appropriate for use in these regions.In spite of this, many successful studies have been published either by using dedicated altimeter products or by, e.g., using in situ data sources to derive the corrections, such as surface pressure (for the dry correction) and wet tropospheric corrections derived from Global Navigation Satellite Systems (GNSS) stations, e.g., [5,7,16].
Figure 1 illustrates the spatial coverage of inland water measurements for Envisat selected from RADS using the criteria explained in Section 2.1.2.To increase the spatiotemporal resolution of the observations, most studies over continental waters make use of multiple mission data, from at least two complementary missions, e.g., one 10-day and one 35-day repeat mission.Unfortunately, with the loss of Envisat in April 2012, well before the launch of the follow up ESA mission, Sentinel-3 (which will be on a 27-day repeat cycle), and until the launch of SARAL in February 2013, for the first time after nearly 20 years, this two-mission data set was discontinued.This increased the use of CryoSat-2 data for applications for which the satellite had not been designed, including studies over ocean and inland water [17,18].
At present, considering the whole error budget (range measurement, orbit, instrumental, range and geophysical corrections; see Equations ( 1) and ( 2)), altimetry measurements over ocean have an accuracy of a few (1-3) centimeters.Over inland water, this accuracy ranges from a few centimeters (e.g., Great Lakes, USA, or Lake Victoria, Africa) to tens of centimeters, depending on the size of the water body and wind conditions [19,20].Comparisons with tide gauge data indicate a root mean square (RMS) accuracy of e.g., 2.6 cm for Lake Victoria [15,19] and a ~3-5-cm RMS for the Great Lakes [19,21].This accuracy decreases to several decimeters for narrow inland water bodies in rough terrain [1,15].Studies performed over small water bodies usually involve the retracking of the altimetric waveforms, which considerably increases the amount of data involved and significantly improves the accuracy of water level time series [22].
In each section of this paper, whenever applicable, the specificities of the corrections for each mission will be discussed, with emphasis on the missions mostly used in altimetry over inland waters: T/P, Jason-1, Jason-2, Envisat, CryoSat-2, SARAL and, to a somewhat lesser extent, GFO and ERS.To facilitate the reference to all corrections and missions discussed in this study, the most relevant information relative to the products present in RADS and in the Geophysical Data Records (GDR) of three of the most relevant altimetric missions (Envisat and Jason-1/2) is presented in Table 1.For the remaining missions, the information will be given in the text as appropriate.The information relative to the T/P GDR is not provided in Table 1, since this dataset is expected to be replaced soon by an updated version.In the Level 2 altimeter products, Interim Geophysical Data Records (IGDR) and GDR, the agencies usually provide measurement data and all corrections at a 1-Hz rate, while some fields (range and satellite position-latitude, longitude and the height above the reference ellipsoid-with respect to an adopted ITRF) are provided at higher rates: 10 Hz (T/P), 40 Hz (SARAL) or 20 Hz (all other missions).The corrections are always provided at 1 Hz and need to be interpolated to the higher data rates by the user.
Compared with T/P and Jason-2, Jason-1 reveals a significant data loss over continental water regions, both due to on-board processing (data filtering) and the way the 1-Hz measurement is built from the 20-Hz measurements [15].For this reason, Jason-1 data are usually restricted to larger lakes.The same happens to GFO, due to its poor performance over narrow inland water regions [15,19].
Unlike the NASA/CNES/EUMETSAT missions, which are primarily aimed at open-ocean, the ESA missions have multidisciplinary objectives, and their methods of tracking are thus more sophisticated.The radar altimeters on board ERS-1 and ERS-2 had two tracking modes, ocean-tracking mode (for sea surfaces) and ice-tracking mode (for ice sheets and sea ice).While ERS-1 alternated between these two modes over land, ERS-2 spent the whole mission lifetime in ice-mode over land.The changing orbit configurations in ERS-1, having been during its life time in three different types of orbits, with different ground tracks, make derivation of time series over inland water impossible over all but a small number of very large targets, whereas the dataset from ERS-2 provides a long time-span over inland water, with the "ice mode" utilized over land greatly facilitating the acquisition of inland water targets at the cost of some vertical precision [23].For Envisat, because of the automatic change of resolution over rougher terrain and the availability of various retrackers (see Table 1), data recovery over continental waters is generally good [15].
River and Lake provides near-real-time water-level products derived from Envisat and Jason-2, using appropriate waveform retracking for each echo shape [23].The RLA product (River/Lake Altimetry), as opposed to RLH (River/Lake Hydrology, intended for non-specialist users), provides all corrections separately, so that the expert user can replace them with local or regional data, if available and necessary.Apparently, the available or applied corrections are those present in the original products, and no dedicated corrections are provided, nor is the source of the original altimeter products given [24].
Two main products are available from LEGOS over land, those from the CTOH and the Hydroweb databases.
The Hydroweb database contains time series over water levels of ~250 large rivers, ~100 lakes and wetlands around the world.For lakes and reservoirs, the water levels are based on merged 1-Hz measurements from T/P, Jason, ERS, ENVISAT and GFO data provided by ESA, NASA and CNES data centers.The classical corrections (ionospheric and tropospheric corrections, polar and solid Earth tides and sea state bias (SSB)) are applied, but it is not specified which corrections (model or measurements) are chosen.For rivers and floodplains, the time series are constructed using the most upgraded T/P GDRs made available by AVISO (Archiving, Validation and Interpretation of Satellite Oceanographic data).The basic T/P 10-Hz measurements are used, corrected for the ionosphere (based on total electron content (TEC) from Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS)), dry tropospheric correction (DTC) from the European Centre for Medium-Range Weather Forecasts (ECMWF) surface pressure fields, and wet tropospheric correction (WTC) from the National Centers for Environmental Prediction (NCEP) air temperature and specific humidity fields [25], solid Earth tide and for on-board instrumental drifts and bias.Further details on the Hydroweb products can be found in [26].CTOH provides along-track Merged Geophysical Data Records (MGDR) for the T/P, Jason-1, Jason-2, GFO and Envisat missions.The MGDR have both 1-Hz and raw data at higher rate (10 Hz or 20 Hz, according to the mission) measurements and up-to-date corrections.The user can extract along-track data both using standard, most up-to-date CTOH corrections, or choose "à la carte" combinations of data and corrections.Dedicated CTOH up-to-date corrections include a WTC at the continents' altitude, interpolated from the gridded NCEP levels [25] and two different ionospheric corrections (IC): a dual-frequency-based one, smoothed with a 20-point median filter, adapted with a variable length approaching the coasts to maintain maximum data coverage and linearly interpolating for small gaps (<5 points); and another from GIM [27].CTOH altimetry products also use radar waveform retracking, allowing the choice between one of four retracking algorithms (ocean, ice1, ice2 or sea ice) for most of the missions available, in particular for Jason-2, T/P and Envisat.From all analyzed dedicated inland water products, CTOH seems to provide the most elaborate dataset, although it is not entirely global, providing data for a set of about 100 large lakes and about 250 sites (called virtual stations) on large rivers.
GRLM provides time-series of water level variations for some of the world's largest lakes and reservoirs; large (>100 km 2 ) lakes in important agricultural regions are the main targets.The altimetric datasets currently exploited are: T/P (AVISO and PO.DAAC (Physical Oceanography Distributed Active Archiving Center) GDR-C), Jason-1 (IGDR datasets available via ftp at podaac.jpl.nasa.gov),Jason-2 (IGDR from both AVISO and NOAA (National Oceanic and Atmospheric Administration), GFO (GDR from NOAA) and Envisat (unspecified data version).These products primarily use 1-Hz data or averaged data of several individual 10-Hz or 20-Hz data without radar waveform retracking.The applied corrections include the DTC, the radiometer-based WTC when valid (and the model-derived correction, when not valid), the DORIS IC, Earth tide and the elastic-ocean and ocean-loading (the last two only applied to the Caspian Sea) [15].Most of the above-mentioned altimeter-derived products are advanced, ready-to-use products in whose generation the user has no intervention.Only for some of them detailed information on how they are generated and the type of range and geophysical corrections that have been applied to the data is easily available, as, for example, the CTOH database, which provides along-track data and detailed information about the available range and geophysical corrections.
Apart from the various products provided by the space agencies, near-real time (NRT), IGDRs and GDRs for general purposes and research, including hydrological studies, there are several data centers providing multi-mission altimetry, such as AVISO, PO.DAAC and RADS.The only limitation of these datasets, when compared with the IGDR/GDR, is the fact that, at present, they all provide data only at a 1-Hz rate.
The Radar Altimeter Database System was developed first in the late 1990s at the Delft University of Technology to create a harmonized, up-to-date database of all then flying altimeters, ironing out known errors in the GDRs and inconsistencies between the various GDR products available.Further developments took place at NOAA, adopting the Network Common Data Format (netCDF), efficient additional software, operational updating of the database and upgrading it to a prototype Climate Data Record for sea level.RADS keeps being developed jointly with EUMETSAT.
For many studies using multi-mission altimetry, RADS provides the most complete, harmonized data set, including a wide collection of, e.g., orbit solutions, mean sea surface (MSS), geoid models, topographic/bathymetric models and range and geophysical corrections (e.g., various ionospheric models [27][28][29][30] and pole tide [31] properly takeing into account different Love numbers for ocean, land and lakes).In addition, when selecting the RADS default parameters, the sea-level time series of the various satellites will be aligned due to the application of the reference frame offset; see Equations ( 1) and (2).For studies requiring higher data rates, only available at (I)GDR, advanced users can still make use of updated corrections, orbits, MSS, geoid and topographic models available in RADS and merge them with the 10-Hz/20-Hz/40-Hz (I)GDR ranges.The same applies to developers of improved ranges through retracking, who can also merge their ranges with the remaining parameters available in RADS.Table 1 provides a summary of the main corrections provided in RADS and which are relevant to this study.An update description of RADS can be found in [32].In the next sections, the most relevant issues related to the atmospheric corrections required for altimetric studies over continental waters will be addressed.Some of the identified issues will be illustrated with examples using RADS data.The RADS data set was selected due to the large collection of range and geophysical corrections that are provided.
The tropospheric corrections (dry and wet) are addressed in Section 2 (2.1 and 2.2, respectively), while the ionospheric correction is addressed in Section 3. In each section, the procedures used in the estimation of each correction are reviewed, the main error sources identified in the current altimeter products are analyzed and examples illustrating the impact of these errors on water level measurements are given.Section 4 presents a summary of the main conclusions and recommendations.

DTC Estimation
The path delay, , due to the dry neutral gases in the atmosphere, the dry tropospheric correction (DTC), in absolute terms, is the largest range correction in satellite altimetry.With an absolute value of about 2.3 m at sea level and a range of about 0.2 m, with almost linear height dependence, it can be estimated with an accuracy of a few millimeters from surface atmospheric pressure using the modified Saastamoinen model [33]: ( where is the surface pressure in hPa, is the geodetic latitude, is the surface height above the geoid (in meters) and results in meters.
Since, in Equation ( 3), is the total atmospheric pressure, i.e., the partial pressure of dry air plus the pressure due to water vapor, this expression actually gives the zenith path delay caused by the hydrostatic component of air and not just by the dry component of atmospheric pressure, the effective "dry delay".The authors of [33] suggested combining the "dry" term part of the zenith path delay and the same effect of water vapor, thus yielding the hydrostatic term that can be calculated from the total atmospheric pressure, data that are more routinely available.Given the small difference between the hydrostatic and dry components of the tropospheric path delay, the term "dry tropospheric delay" is usually used within the altimetry community to refer to the hydrostatic tropospheric path delay.According to [13], the expression Equation ( 3), when combined with the corresponding expression for the wet path delay given by [33], is accurate to better than 0.2%, i.e., for a total tropospheric delay (dry + wet) of 2.5 m, the error is < 5 mm.
The most common sources of atmospheric pressure are the atmospheric models from the European Centre for Medium-Range Weather Forecasts (ECMWF) [34] and the National Centers for Environmental Prediction (NCEP), [35].All models provide global grids of sea level pressure (SLP) and surface pressure (SurfP), every six hours at different spatial samplings.Two ECMWF products are of particular interest: the operational model at 0.125 × 0.125° (regular grid) or about 16 km (Gaussian grid) and the latest reanalysis product, ERA Interim [36], which is provided at 0.75° × 0.75° (regular grid) or about 80-km (Gaussian grid) resolution (e.g., http://www.ecmwf.int/products/forecasts/guide/user_guide.pdf).For simplicity, throughout the text, in particular, in the figures, whenever the term, ECMWF, is used, it refers to the ECMWF operational model.NCEP provides similar products at various spatial samplings, and for example, corrections present in RADS have been derived from 2.5° × 2.5° grids (see Table 1).This choice for NCEP is historical, since it is the model used as far back as Geosat and remained consistent over time [35], thus avoiding discontinuities and drifts present in the time series of the ECMWF operational model.Each of these models provides grids at six hourly intervals.
Since most satellite altimetry applications are over ocean, the most common approach adopted in the computation of the dry tropospheric correction at each satellite measurement location is to perform the computation in two steps.First, the correction is computed at sea level from SLP grids, using Equation (3) and setting equal to zero.The correction is then reduced to surface height by performing a suitable height reduction.
The height reduction of the DTC can be performed by means of the height dependence of atmospheric pressure, which, according to the International Standard Atmosphere [37], is modeled as: (4) where is the normal lapse rate of temperature (-0.0065 ), is mean sea-level temperature (288.15K), is the mean gravity (9.80665 ) and is the specific constant for dry air (287.053 ).By using all constants as above, Equation ( 4) can be written as: (5) A similar expression, with slightly different coefficients, given by [38], has been used by many authors (e.g., [39]).These two expressions lead to surface pressure values that differ by less than 0.3%.
In Equation (5), and are the atmospheric surface pressures (in hPa) at heights and (in m), respectively.Since is the sea level pressure, then and is the surface orthometric height.
Equations ( 4) and ( 5) are simplified expressions for the height dependence of atmospheric pressure in which any variation of sea-level temperature is neglected.More precise expressions can be found in the literature, most of them equivalent to the formula by [40], adopted by authors such as [41][42][43]: - In Equation ( 6), is the specific constant for dry air, is the mean temperature (in K) of the layer between heights and and is the mean gravity, as given by: (7) can be estimated as the mean value of temperatures and at heights and , respectively, obtained, for example, from the values of at mean sea level given by the Global Pressure and Temperature (GPT) model [41] and considering a value of -0.0065 for the normal lapse rate of temperature with height (for the estimation of ).
In [43], the effect of the height reduction on the estimation of the dry path delay, performed either by Equations ( 5) or ( 6), was evaluated, for a set of GNSS coastal sites with heights up to 1000 m. Results show that the global grids of sea level pressure provided by the ECMWF operational model, either at 0.25° or 0.125° spacing, or the ERA Interim reanalysis product at 1.5°, allow the estimation of the hydrostatic component of the tropospheric delay with an accuracy of 1 to 3 mm at the global scale, provided an adequate model for the height dependence of atmospheric pressure such as given by Equations ( 6) and ( 7), is adopted.It was also shown that expressions such as Equation ( 5), will induce seasonal signals with amplitudes of several millimeters, due to the seasonal variation of pressure with temperature.For higher altitudes, as in several continental water regions, larger effects can be expected.
The six-hour SLP or SurfP global grids commonly available from ECMWF or NCEP do not permit proper resolution of the tidal variations in atmospheric pressure, the air tides, particularly the semi-diurnal solar tide, S2 [44].As a consequence of the six-hour temporal sampling being the S2 Nyquist frequency, S2 appears as a standing wave rather than a known westward propagating wave, with much underestimated amplitudes near longitudes where sampling times happen to coincide with times of the S2 zero values or nodes [45].A strategy for handling tidal signals in (re)analysis products, such as the six-hour SLP grids from ECMWF or NCEP, for routine geodetic and oceanographic data processing, was developed by [46].A similar procedure has been implemented in RADS, so that the air tides are properly accounted for in the computed DTC (see Table 1).However, this effect is very small, only up to ±3 mm, and could be neglected for most inland water applications.
The DTC has strong height dependence, expressed by Equations ( 5) or (6).Since satellite altimetry is primarily designed for ocean applications, for which no such dependence exists, some altimeter products fail to provide the DTC appropriate for inland water studies.This problem has been pointed out by several authors [5][6][7]15] who reported large differences in the dry tropospheric correction of some satellite products.

Analysis of DTC Errors Present on Altimetric Products
In this section, we aim to inspect how the dry tropospheric correction is presently handled on inland water measurements, in the products of the various satellite missions, using the RADS products as examples.For this purpose, the values present in RADS derived from the three following models were analyzed: the ECMWF operational model, ERA interim and NCEP.The three DTC data sets present in RADS for each of the above-mentioned models have been summarized in Table 1 and are further detailed below: (1) For ERA Interim and NCEP, for all altimetric missions, the SLP is first interpolated in space and time from the available model grids and then reduced to surface height using Equation ( 6), while setting = 0.The surface height has been extracted from the DTU10 topographic data set, except for the Caspian Sea, where a constant value of −27 m was adopted.DTU10 are upgrades of the corresponding DNSC08 models [47].The DTU10 topographic data set contains topographic information as follows: over land and most inland water regions, it contains the surface height above the geoid as modeled by a digital elevation model (DEM), while over oceans and over the Caspian Sea, it possesses bathymetric information.
(2) For the ECMWF operational model, the DTC values are those extracted from the GDR products of each altimeter mission, without any modification, except for CryoSat-2, for which the correction has been computed from the model grids using the same procedure as in (1).For Geosat and GFO, only the DTC from ERA Interim and NCEP models are available.For T/P, J1 and J2, the AVISO GDRs have been used.
For use in the subsequent analysis, 1-Hz altimeter measurements over inland water have been extracted from RADS using the following criteria: points over inland waters only using the appropriate data flags; points for which the RMS of the 10-/20-/40-Hz measurements (depending on the mission) is less than 0.4 m; points for which the absolute difference between the measured altimeter height above the geoid and the corresponding height given by the DTU10 topography model is less than 500 meters.For this purpose the orthometric height above the EGM2008 geoid [48] ( ) was computed by subtracting the EGM2008 geoid height above a reference ellipsoid (N) from the water level height above the same ellipsoid (h) given by Equation (1).The last two criteria aim to remove points with the largest errors in the altimeter measurements.With this selection, we aim to have a representative set of measurements for the various conditions that affect the atmospheric corrections to be inspected.Figure 1 illustrates the location of these points for Envisat.
Figures 2 to 7 illustrate how the three DTC datasets present in RADS altimetric products handle the height dependence of the correction and the dramatic impact of this procedure on the products' accuracy over inland waters.
We find that the DTC has a mean variation of 2.5 cm per each 100 m (see Figures 2 and 3).While at sea level, the correction is about −2.3 m, at a height of 4000 m it only amounts up to −1.4 m.
Figure 2 represents the DTC from three models versus surface height above the geoid for various altimetric missions, showing that the various products handle the height dependence of the correction in different ways.Figure 3 illustrates the three DTC data sets for a portion of Jason-1 Cycle 62.The x-axis is the along-track point number (1-Hz measurements), but showing only inland water measurements as extracted according to the criteria described above.In this figure, all points correspond to inland water measurements, shown continuously, as if there were no land or sea in between these points.Various satellite passages over some of the largest lakes, such as the Great Lakes and Lake Victoria, Africa, can be identified.
It can be observed that the DTC derived from ERA and NCEP correctly follow the surface topography, but the correction derived from the ECMWF operational model, extracted directly from the GDR products of the various missions, is not provided in a consistent way.
For Jason-1, Jason-2, Envisat (GDR version 2.1), CryoSat-2 and SARAL, the ECMWF-derived correction is provided at the surface height (Figure 2(top)), while for T/P and for all ERS-2 and ERS-1 analyzed cycles, the correction seems to be provided at sea level without any height reduction.This is evident in the flat red points in the bottom graph of Figure 2.For these satellites, the correction was provided by the program CORIOTROP (CORrections IOnosphé riques et TROPosphé riques Franç aises), running at Mé té o-France.Unfortunately, this process was not part of the T/P algorithm specifications, so we cannot verify our assumption.For T/P and J1, the plots refer to all cycles of Phase A of each mission, i.e., the period when the satellite was in the so-called reference orbit.Together with Figure 3, Figure 4 helps to illustrate another type of DTC error present in some products.Figure 4 shows how the DTC is handled in Jason-1 GDRs, for Pass 222 over Lake Tanganyika, located in central Africa, with a mean height of 769 m.It shows that interpolation errors arise from the fact that the DTC has been computed using surface pressure grids.Surface pressure is provided by ECMWF at the level of the atmospheric model orography, which is similar to a smoothed DEM, and can significantly depart from the actual surface height by values up to ~500 m.Consequently, for regions, such as Lake Tanganyika, where the difference between lake height and the model orography amounts up to 300 m, errors up to 7-8 cm can occur.Figure 4 is representative of the interpolation errors, which occur on Jason-1 cycles up to 149 (see the further explanation below and the bottom panel of Figure 5).For all of these cycles, very similar plots are obtained for the same pass over Lake Tanganyika.This type of error also explains why in Figures 5 and 6, the ECMWF DTC is always above (is less negative than) the corresponding ERA and NCEP values.
Figure 5 shows the mean value of the dry correction provided by each of the three models for each cycle of T/P Phase A (top) and Jason-1 Phase A (middle), respectively.Figure 6 shows the equivalent results for two versions of the whole set of Envisat GDRs. Figure 5 shows that while for all Jason-1 cycles, the DTC from the three models is provided at the surface height (see Table 1), for all T/P cycles, ECMWF-derived DTC are given at sea level.Those from ERA and NCEP are given at the surface height for both satellites (see Table 1).
Figure 6 (top) illustrates that the DTC was not consistently modeled over the various cycles of Envisat GDR version 1: while for Cycles 6 to 10, 48 and all cycles after 63, the correction was provided at surface height, for Cycles 11 to 62 (except Cycle 48), the correction appeared with most of the values at sea level (−2.3 m), but surprisingly, values significantly different were also present (up to −1.9 m).While it is known that different types of processing have been applied in the generation of the GDRs for the various Envisat version 1 cycles, the reason for having mixed DTC values (at sea level and at the surface) on the same cycle is not understood.The bottom plot of Figure 6 shows that in the most recent Envisat data (GDR version 2.1), the DTC values are now consistently provided at the surface height (see Table 1).
The bottom plot of Figure 5 represents the mean cycle values of the DTC for Jason-1 Pass 222 over Lake Tanganyika, showing that different processing has been applied up to Cycle 149 and from Cycle 150 onwards.Until 31 January 2006 (the end of J1 Cycle 149), the correction is provided at the surface height, derived from surface pressure grids, as shown in Figure 4. On 1 February 2006, ECMWF introduced a new processing version, which was likely paired with a software update at Mé té o-France to use the proper elevation from a DEM, thus eliminating the errors associated with the previous approach.Unfortunately, no documentation of the processing at Mé té o-France is available, nor was that process redone during subsequent updates of the Jason-1 GDRs in 2009.
In view of these results, we concluded that the ECMWF-derived DTC is not provided homogeneously in the GDRs of the various missions.On the contrary, the RADS DTC from the ERA and NCEP models are homogeneous corrections for all missions, correctly estimated at the water surface level (Table 1), though limited by the resolution of the underlying DEM.
The DTC errors, over inland water regions, present in the current altimeter products, can be grouped into three types: (1) The DTC has been computed from SLP grids, with an expression similar to Equation (3) with no further reduction to the water surface level.For each location, this error is mainly a bias, which, for inland water at high altitudes, such as the Lake Titicaca (~3800 m), can reach 80 cm.Apart from a constant part, this error also has a seasonal component, since, according to Equation ( 6), an accurate height dependence of surface pressure can only be modeled by accounting for the seasonal variations of atmospheric pressure with temperature [43].
(2) The DTC is estimated from surface pressure grids at the level of the model orography.This will induce interpolation errors, due to the fact that the model orography might depart significantly from the actual surface height.In regions of large variations in terrain height, such as a river in a deep valley or a lake surrounded by high mountains, even if the model orography correctly fits the terrain heights, there may be a large difference between the DTC at the points over the mountains and the DTC at the center of lake or river.This may cause along-track interpolation errors at the decimeter level, with a linear decrease/increase in the DTC as the points approach the center of the lake, reported as "V shape" errors by, e.g., [5], and illustrated in Figures 3 and 4 for Lake Tanganyika.
(3) The DTC has been computed from SLP grids, with an expression, such as Equation ( 3), with further reduction to the surface water level, for each along-track altimeter location, using an appropriate height reduction.Studies by [43] show that, for heights up to 1000 m, Equation ( 6) provides the DTC with an accuracy of a few millimeters if an accurate DEM is adopted and the GPT climatological model is used to model the seasonal dependence of atmospheric pressure on temperature.Further research is required to find out how these results can be extended to larger heights, up to 4000 or 5000 m, since a few continental water regions of interest exist at these altitudes.Preliminary results indicate that for these altitudes, errors due to the inaccurate modeling of both the height and the temperature dependence of surface pressure are expected.
Figure 7 further illustrates the importance of using correct height information in the required height reduction of the DTC.This figure illustrates the DTC (top, in m) and surface topography from DTU10 (bottom, in m) versus the point number by ascending time order, over the Caspian Sea, for SARAL Cycles 1 and 2. The ECMWF Operational DTC values are those present in RADS, extracted from the (I)GDRs.The blue line represents ERA DTC computed using a mean height of −27 m, while the black line shows the ERA DTC values obtained using the DTU10 topography model in the computation of the height reduction.The difference between the red and blue lines (a mean value of 7 mm) reveals that the ECMWF operational DTC, present in the (I)GDRs, was computed using a mean height of 0 m.
In summary, to avoid most of the reported errors over inland water regions, the DTC for satellite altimetry should be computed as follows: (1) For each atmospheric model SLP grid node, compute the correction at sea level using Equation (3) and setting = 0.
(3) Apply height reduction using Equation ( 6) and a model topography, such as DTU10 (apart from the Caspian Sea, where a mean sea height must be adopted).Users may also use a local and more precise DEM where available.For studies over lakes or reservoirs, for example, it will be preferable to use the mean lake height instead.In the implementation of Equation ( 6), the surface temperature can be obtained from the values of the temperature at mean sea level given by the GPT climatologic model and considering a mean normal lapse rate of temperature with height of −0.0065 K/m.For non-time-critical inland water studies, the authors recommend the use of the correction from the ERA Interim model, which is more stable in time, although both the ECMWF operational and NCEP models provide enough accuracy, if the appropriate computation procedure is adopted.Alternatively, whenever measurements of surface pressure are available, from reliable and well-monitored barometric stations, these can be used to compute the DTC from Equation ( 3) with an accuracy of about 1 cm [7,43].

WTC Estimation
The path delay due to the presence of water vapor in the atmosphere, the wet tropospheric correction (WTC), is one of the major error sources in satellite altimetry.With an absolute value less than 50 cm, it is highly variable, both in space and time.Due to this high variability, over the ocean, the most accurate way to model this effect is through the measurements of microwave radiometers (MWR) on board the altimetric missions.
The passive MWR on board altimeter missions retrieve the water vapor content from the instantaneous measured brightness temperatures, at the nadir, in channels inside and outside the water vapor absorption line centered at 22.235 GHz.An analysis of the errors associated with the WTC estimation from on-board MWR of various missions can be found, e.g., in [49].
Three-band radiometers have been used in the NASA/CNES/NOAA/EUMETSAT missions (T/P, Jason-1 and Jason-2): one at the water vapor line at 18-18.7 GHz (low sensitivity to clouds), one in the absorption line (21-23.8GHz) and one in the 34-37 GHz band (as sensitive to the surface as the low frequency one, but more sensitive to cloud liquid water).
In the ESA missions (ERS-1, ERS-2 and Envisat) two-band radiometers have been used: 23.8 GHz and 36.5 GHz.GFO also carried a dual frequency MWR at 22 GHz and 37 GHz [50].SARAL also carries a two-band radiometer (23.8 and 36.8GHz) similar to that of Envisat, sharing the AltiKa antenna and allowing a smaller footprint [51].
CryoSat-2 does not carry a radiometer, thus relying on models, such as those from ECMWF.New approaches are being developed to determine the WTC for CryoSat-2, based on the combination of all available data types: from imaging radiometers on board remote sensing satellites, from models and GNSS-derived wet delays from inland stations [43,52].
The algorithms used to retrieve the WTC from the measured brightness temperatures of the various MWR channels assume a constant surface ocean or calm waters emissivity.Therefore, in the presence of surfaces with different emissivity, the measurements become invalid.This is the case of land, vegetation or ice surfaces.
The problems associated with the computation of the wet tropospheric correction in the coastal regions, where the MWR measurements also become invalid, in the context of missions possessing an on-board MWR, has been addressed by several authors [53][54][55][56].These authors have derived algorithms for improving the WTC in the coastal regions, but these algorithms have not been customized for inland water regions.
Over continental waters, the wet path delay retrieval from passive microwave radiometers is hampered by the contamination on the radiometer measurements of the surrounding lands (e.g., in lakes or rivers), mixed vegetation in wetlands or ice at high latitudes.Because of its large footprint, the MWR measurements over an oceanic or lake surface become invalid, due to emerged lands, within distances of about 20-40 km of the water area limit, depending on the frequency.Consequently, the MWR correction provided in the GDRs is usable only in the central parts of very large lakes.In the case of smaller lakes (<2000 km 2 ), the footprint of the radiometer always encompasses surrounding lands or islands, and then, the MWR-derived WTC is either not available or invalid.Thus, MWR-derived WTC cannot be used in small lakes or over river channels.In these regions, alternative sources must be used.
For regions possessing permanent GNSS stations, GNSS-derived wet path delays can be derived with an accuracy close to 1 cm, the same accuracy as MWR-derived WTC (e.g., [55]).This approach is particularly favorable for small lakes and reservoirs, where the measurements at a single location can be representative of the whole lake, the approach followed by, e.g., [7].However, unlike the coastal regions, where GNSS data have been used in combination with other available wet path delay sources (valid MWR measurements in the vicinity of the point and model values) [55], so far, this procedure has not been exploited over inland water regions.
In the absence of any other data source, the WTC from meteorological models must be used.Although the overall accuracy of WTC from meteorological models is worse than the corresponding values from MWR, the quality of the recent models has been increasing significantly [34], particularly for the latest reanalysis product from ECMWF, the ERA Interim model [36,57].
For use in satellite altimetry, the WTC can be calculated from global grids of two single-level parameters provided by global atmospheric models, the total column water vapor (TCWV, expressed in mm or, the equivalent, kg/m 2 ) and near-surface air temperature (two-meter temperature, ), from the following expression [58,59]: (8) where is the mean temperature of the troposphere, which can be modeled from , according to, e.g., [60,61]: (9) Equations ( 8) and ( 9) provide the WTC at the level of the atmospheric model orography.This orography may depart from the actual surface height by hundreds of meters.If accurate surface heights are available, the corresponding value at the surface can be obtained by an appropriate height reduction.
Unlike DTC, which has a relatively well-known height dependence (see Section 2.1), the height dependence of water vapor is not easy to model, due to its large variability.In spite of this, Kouba [39] proposed the following empirical expression: (10) where and are the ellipsoidal heights of the model orography and surface, respectively.Thus, height errors of, e.g., 100 and 500 m, will induce errors of 5% and 28%, respectively (1 cm and 5.6 cm for a correction of 20 cm).Although these errors are considerable, in the altimetric products, the model-derived WTCs are commonly referred to the model orography or given at sea level, since in some regions it is difficult to have reliable DEM.According to [39], Equation (10) should only be used to performed WTC height reductions up to 1000 m.
In summary, for each model grid node, the correction is: (1) computed at the orography height using Equation (8), followed, or not, by a height reduction to the surface height by Equation ( 10); (2) interpolated at the altimetric measurement location; and (3) interpolated at the measurement time using the two closest model grids, six hours apart.
Alternative approaches for the computation of the wet path delay from other parameters provided by the global atmospheric models can be adopted, e.g., from 3D parameters (specific humidity, near-surface air temperature and surface pressure; see, e.g., Equation (1) in [62]).In this case, the integration over height can be performed from the top of the troposphere until the surface.In [62], a computation algorithm is proposed, where the effective thickness of the atmosphere column is deduced from the altimetric measurement itself.

Analysis of WTC Errors Present on Altimetric Products
Similarly to the analysis performed in Section 2.1.2for the DTC, in this section, we present an evaluation of the WTC available on altimeter products, using the values provide in RADS as examples.
As summarized in Table 1, for each mission, RADS provides the wet tropospheric correction from the on-board MWR (if available) and from three models: ECMWF operational, ERA Interim and NCEP (apart from Geosat and GFO, where the values from the ECMWF operational model are not available, as they are not present in the corresponding GDR).In addition, as mentioned before, the MWR correction is not available for CryoSat-2.
As for the DTC, the WTC from ERA Interim and NCEP present in RADS have been calculated from the model grids, at 80 km and 2.5° spatial sampling and six-hour intervals.The values have been computed at the model orography height using Equations ( 8) and (9).For the ECMWF operational model, the values from each mission GDR are given without any change, except for CryoSat-2, for which the correction has also been computed directly from the model grids.
In the following analysis, a selection of points over inland water was made based on the same criteria as described in Section 2.1.2.In addition, in the analyses involving the measurements from the on-board MWR of each mission, only points for which the radiometer flag indicates no corruption by land were considered.
Figure 8 shows the mean WTC from each of the three analyzed models, for each 10-day cycle of T/P (Phase A), Jason-1 (Phase A) and Jason-2.The dominant feature is the strong seasonal signal of this correction, with a 12-15-cm peak-to-peak amplitude.Figure 9 shows the mean (top) and standard deviation (bottom) of the differences between the WTC from ERA and each of the other two models, for each 10-day cycle of T/P (Phase A), Jason-1 (Phase A) and Jason-2.Both figures show that the WTC from ECMWF is not provided in a consistent way in the GDRs of T/P, J1 and J2, the so-called reference missions.While for T/P, it is provided at sea level, for Jason-1 and Jason-2, they are correctly determined by integrating water vapor from the surface upwards (Table 1).These results also show that, when computed at the same level as is done for Jason-1 and Jason-2, the differences between the ERA and ECMWF operational model decrease with time, with a mean of 0.4 cm and a standard deviation of 1.5 cm for the whole Jason-1 and Jason-2 period.These values are in agreement with the results obtained for the global comparison of the two models (not shown here), which shows that the ECMWF operational model has been updated and improved several times over the years.The global difference between the two models has a standard deviation, for each six-hour grid, that goes from 2.5 cm in 1995 to about 1.1 cm after 2004.The corresponding values shown in Figure 9 are somewhat larger than the global results and might be due to the different sampling and to the fact that the WTC from the two models might have been computed using different formulae.Results also show that the differences between ERA and NCEP are larger than the corresponding differences between the ERA and ECMWF operational model (considering only these last ones for the Jason-1 and Jason-2 period).The differences between ERA and NCEP are stable in time, with a mean value of −0.6 cm and a standard deviation of 2.3 cm for the whole set of Jason-1 and Jason-2 cycles.
Figure 10.Wet tropospheric correction (WTC) for all T/P measurements over Lake Victoria (Africa) for five years (1993)(1994)(1995)(1996)(1997)(1998).Only measurements over inland water with valid microwave radiometer (MWR) correction were considered.  2 and 3 show the comparison between the WTC from three models and from the on-board MWR of each mission for Lake Victoria (Africa) and the Caspian Sea.These are examples of two regions with significantly different conditions for various range corrections, including the WTC.
Lake Victoria is one of the world's largest lakes, located in Africa, with a mean height of about 1135 m.Over this lake, the WTC has a weak seasonal signal, with values between 10 cm and 30 cm all over the year (Figure 10).Over this lake, all models capture very well the main wet path delay variations revealed by the MWR correction, except for some wetter conditions.
Results show that the T/P correction from the ECMWF operational model is too negative by about 10 cm (Figures 10 and 11 and Table 2), since it is provided at sea level instead of lake height.The authors of [5] reported a similar behavior of the ECMWF WTC on the T/P GDR, which also provided the correction at sea level, derived from global 3D grids.
From the three models, ERA is the one providing the best results, with a mean difference with respect to MWR, considering the whole analyzed period from 1992 to 2013, of 0.0 cm and a standard deviation of 2.3 cm.The corresponding statistics for the difference between MWR and NCEP are −3.4 cm for the mean and 2.9 cm for the standard deviation.Thus, although in terms of variability, ERA and NCEP seem to be similar over Lake Victoria, NCEP has a bias of about 3 cm with respect to ERA and MWR.The Caspian Sea has a mean height of −27 m.The most striking feature of the WTC over this closed sea is its strong seasonal signal with a peak-to-peak amplitude close to 15 cm (Figures 12 and 13).

Models
Regarding the ECMWF operational model, due to the small sea level height, the fact that the correction is provided at sea level does not introduce any significant error.
As for Lake Victoria, from the three models, ERA is the one that provides the best results (Figure 13 and Table 3), with a mean difference with respect to MWR of 0.3 cm and a standard deviation of 1.7 cm, considering the period from 1992 to 2013.The corresponding statistical parameters for the difference between MWR and NCEP are −0.5 cm for the mean and 3.2 cm for the standard deviation.Figure 13 shows that the amplitude of the differences of NCEP with respect to MWR is considerably larger than those for the other two models, showing that NCEP is not able to model the large seasonal variations present in the WTC of the Caspian Sea.The results obtained for the ECMWF operational model are of the same order of magnitude of those obtained for ERA, showing that when provided at the same level and for the period after 2004, these two models provide similar results.
Concerning other altimetric missions, no systematic differences were found between the WTC from the three analyzed models for Envisat, CryoSat-2 and SARAL, thus indicating that these corrections are provided at surface height and should have an accuracy close to that of Jason-1 and Jason-2.As shown in Tables 2 and 3, for Envisat and, generally, for all 35-day missions, the standard deviations of the differences between the MWR correction and each model are larger than the corresponding values for the reference missions by a factor of 1.5-2.0.This can be explained by the fact that these missions have many more passes over these regions, thus spanning a larger range of conditions and mostly because the MWR correction of these missions is much more affected by land effects than the corresponding corrections for the reference missions.Thus, these larger standard deviations cannot be just attributed to model errors, but also to MWR land contaminated measurements that have not been properly edited by just using the radiometer land flag.In spite of the consistency between the three models for all these missions, for Envisat, before 2004, the accuracy of the correction from the ECMWF operational model gets worse, while ERA Interim and NCEP have the same accuracy for the whole missions lifetimes.Table 3. Statistical parameters (in cm) of the comparison of the various WTC for the Caspian Sea.Values refer to the differences for the following periods of each mission: T/P (Phase A), J1 (Phase A) and J2 (Cycles 1-183) and Envisat (EN, Cycles 10-93).For all inland water measurements of the ERS-2 and ERS-1 analyzed cycles over Lake Victoria, a bias around 10 cm was found between the ECMWF operational model and ERA Interim, as reported for T/P, suggesting that for these satellites, the WTC from ECMWF operational model is provided at sea level.

Models
In summary, considering the best available atmospheric model, ERA Interim, results based on the differences between the WTC from the on-board MWR and the model-based wet path delay indicate that the WTC can be estimated over continental waters with an overall accuracy (one standard deviation) of 1 to 3 cm, depending on the region.Moreover, like all models, ERA may have local biases of 1-3 cm.For studies using data only after 2004, the best option is the ECMWF operational model, provided the correction is properly estimated at surface height, due to its higher spatial resolution compared to ERA Interim.The limitation of ERA Interim is its spatial resolution (~80 km), which may not resolve short wavelength variations in the WTC, particularly critical in regions of strong water vapor variability, such as the tropics.
The ECMWF operational model provided in the GDR products of the various missions is not provided consistently and must be used with caution, particularly on T/P and ERS products.Here, the ECMWF-derived WTC is provided at sea level, and the correction by the user to surface height using Equation ( 7) and a suitable DEM does not provide enough accuracy for heights above 1000 m [39].
For inland water studies, the authors recommend the use of the MWR-derived wet tropospheric corrections, using appropriate editing, where available (over large lakes and closed seas), or GNSS-derived wet path delays whenever GNSS stations exist in these regions.
In all other cases, in particular over rivers, small lakes and wetlands, the ERA Interim model, computed at the model orography, as described in this section, can be used.Alternatively, as suggested by [62], the correction can be estimated at surface height using 3D parameters.In this case, the integration can be performed over height, from the surface to the sensor, and more realistic surface heights can be used, for example those derived from the altimeter measurements.This procedure has the advantage that the height dependence of the WTC is modeled with higher accuracy, but requires a larger computational effort, and altimeter-derived surface heights must be used with caution, with an appropriate editing criterion.

Examples of Case Studies
Aiming to illustrate the effective impact of the tropospheric correction errors on altimeter water level measurements, in this section, two case studies are presented: for Lake Victoria and for Lake Tanganyika.Since the major tropospheric errors were found in the T/P products, data from T/P, J1 and J2 have been used in these examples.
For the so-called reference orbit, each mission has mainly only one pass, Pass 120 over Lake Victoria and Pass 222 over Lake Tanganyika.Since only 1-Hz measurements were used from the RADS database, the relatively small number of available measurements limits the results.For this reason, in the water level time series, instead of the mean, the median values are represented.
For each of the mentioned lakes, altimeter data have been extracted from RADS using the criteria described in Section 2.1.1.For each measurement point, three sets of water level heights above the EGM2008 geoid were computed, using the default RADS orbit and a common set of corrections: ionosphere (JPL_GIM), solid Earth tide, load tide (GOT4.8)and pole tide.Then, for each data set, the dry and wet tropospheric corrections were computed from each of the three analyzed atmospheric models: ECMWF operational, ERA Interim and NCEP. Figure 14.Median water level anomalies (m) for Lake Victoria, Africa, derived using tropospheric (dry + wet) corrections from three models, for each cycle of T/P Phase A, Jason-1 Phase A and Jason-2.Water level heights are with respect to the EGM2008 geoid.A mean lake height of 1135 m has been removed.

Figure 15.
Median water level anomalies (m) for Lake Tanganyika, derived using tropospheric (dry + wet) corrections from three models, for each cycle of T/P Phase A, Jason-1 Phase A and Jason-2.Water level heights are with respect to the EGM2008 geoid.A mean lake height of 769 m has been removed.Figures 14 and 15 illustrate the impact of the corrections on the determination of water level time series.As discussed in the previous sections, the fact that in the T/P GDRs, the tropospheric corrections from ECMWF are provided at sea level causes a bias (39 cm for Lake Victoria and 26 cm for Lake Tanganyika).The actual errors are not exactly constant, since, e.g., the height dependence of the wet tropospheric correction is a function of the correction itself, Equation (10); therefore, for wetter conditions, the errors are larger than for dryer conditions.However, when computing the mean for, e.g., a 10-day cycle, the errors nearly average out.When deriving time series from various missions, for which the corrections are not consistent, the impact on the water level time series is evident, with a significant impact on the long-term water level trend.This effect is larger for Lake Victoria, since it is at a higher elevation than Lake Tanganyika.  Figure 16 shows the impact of providing the corrections at sea level instead of at lake height, as happens for the ECMWF corrections from T/P GDR.Apart from a large bias of ~26 cm, the difference between the water level height computed, e.g., from ECMWF and ERA Interim, has a pattern related to the ERA Interim orography used in the WTC computations from this model.Results for NCEP are not shown, since they are very similar to those from ERA Interim.
Figure 17 illustrates the errors associated with the computations of the DTC from surface pressure grids.In this case, the pattern of the differences is due to the difference between the ECMWF model orography and the actual mean lake height.The amplitude of these errors can reach the decimeter level (see also Figure 4).Figure 18 shows the impact of the changes introduced in the J1 processing after the end of Cycle 149, largely correcting the previously mentioned errors (see also Figure 3, bottom plot).In summary, Figures 16 to 18 show that the errors introduced by the inadequate computation of the tropospheric corrections are not just a bias and may have large along-track variations.These types of errors may have a strong impact on comparisons of altimeter data with in situ data, such as, e.g., in altimeter calibration sites over lakes.

Estimation of the Ionospheric Correction
As part of the range estimation for satellite altimetry, which enables the inland water level (or stage) to be estimated, the refraction caused by the presence of electrically-charged particles in the atmosphere has to be accounted for.Most of those charged particles occur in the so-called ionosphere, the part of the atmosphere between an altitude of ~50 and 2000 km, where ions are produced by the photoionization of atomic and molecular gasses [63].
The total delay affecting the altimetric radar pulse along its path through the ionosphere, , can be estimated from the integral over height, from the target to the sensor, of the ionospheric refractivity , which is, in turn, proportional to the density of electrons in the ionosphere, .Unlike its hydrostatic and wet counterparts, the ionospheric refraction is frequency-dependent, i.e., the ionosphere is a dispersive medium for the propagation of the radar pulse, the ionospheric refractivity being inversely proportional to the square of the frequency, , of the propagating signal (e.g., [13]): (11) where , is expressed in electrons/m 3 , in Hz and results in meters.The last integral in Equation (11) represents the atmospheric columnar electron density, i.e., the number of electrons per unit area in a column extending from the surface of the Earth to the altimeter.This total electron content (TEC) is usually expressed in TECU (TEC units), with 1 TECU = 10 16 electrons/m 2 , allowing the altimetric ionospheric path delay to be more practically written as in Equation ( 12): (12) where , is expressed in GHz, and is the total electron content below the altimeter (in TECU), yielding in meters.For altimeters operating in the Ku-band (approximately 13.6 GHz), this comes out to 2.18 mm of path delay per TECU (e.g., [64]).Operating on two different frequencies, recent high-precision radar altimeters flying on Envisat, Jason-1, Jason-2 and TOPEX allow, in principle, direct estimation of from the difference between the measured ranges on the two frequencies [28]: (13) The ranges, and , in Equation ( 13) are those measured on the primary Ku-band and on the secondary C-band (5.3 GHz for TOPEX and Jason) or S-band (3.2 GHz for Envisat), respectively, neglecting the effects of atmospheric refraction.The above formulation has been the primary source of estimation for the ionospheric path delay for ocean altimetry performed by dual-frequency altimeters, as the necessary correction can be retrieved simultaneously, i.e., coincident in space and time, with the range measurements.Unfortunately, the range measurements at the secondary frequencies are inherently less precise than those at the Ku-band frequency, such that the dual-frequency path delay correction would create a much noisier sea level determination.Therefore, in order to reduce the cumulative error arising from the instrumental noise at the two radar frequencies, the dual-frequency retrieved ionospheric correction is usually smoothed with the use of an along-track low-pass filter over about 140-200 km (approximately 20-35 s), as first suggested by [65].
Single-frequency altimeters (such as CryoSat-2, ERS-1, ERS-2, Geosat, GFO and Poseidon) must rely on independent estimates of for the ionospheric range correction, as those available from observational or climatological global models.AltiKa, on board SARAL, is a special case, as unlike all other single-frequency satellite altimeters, it is operating in the Ka-band (35.75 GHz) and is thus about seven times less affected by the ionosphere than Ku-band altimeters, approximately 0.31 mm of path delay per TECU, as determined by Equation (12).
The first attempts to model the ionosphere were based on the development of climatologies, e.g., the Bent model [29,30] and the International Reference Ionosphere (IRI) (e.g., [66][67][68]).Neither of these climatologies ever included GNSS data.However, just like the dual-frequency altimeters, the dual-frequency GNSS measurements provide an abundance of TEC observations along the slant-ranges, and since 1998, institutions like the Jet Propulsion Laboratory (JPL) and the University of Berne started the operational production of global ionosphere maps (GIM) of vertical TEC based on those GNSS slant-range observations.At present, JPL GIM are produced from a global network of about 170 stations and 31 GNSS satellites that contribute their data through the International GNSS Service (IGS) and other institutes [69,70].Since 28 August 1998, two-hourly GIM maps are available from JPL with a single-layer gridded solution of 2.5° spacing in (geocentric) latitude and 5° in longitude [70,71].More recently, NIC09 (NOAA Ionosphere Climatology 2009), a global climatological model based on JPL GIM for the period 1998-2008, has been developed by [28], which can provide global TEC maps with the same spatial and temporal resolutions as JPL GIM, allowing, in addition, the estimation of ionospheric path delays prior to 1998.While GIM uses actual measurements of the total electron content (TEC) along slant-paths through the atmosphere and then models the vertical TEC from that, the NIC09 uses a proxy of the solar activity (global mean electron content or a combination of solar flux and geomagnetic index), a single value for a given time, and then models the global distribution of the vertical TEC based on season and time of day.
In order to use GNSS-derived TEC to estimate the ionospheric correction for altimeter range measurements, besides the necessary interpolation in space and time to the altimeter ground track, additional altitude scaling has to be performed [27], since GNSS measurements integrate TEC to an altitude of about 20,200 km, whereas the altimetric satellites fly only at 1350 km (T/P, Jason-1 and Jason-2) or 800 km (most others).The approach by [27], based on the climatological model, IRI95 [72], suggests that there are essentially no free electrons above 1400 km, and thus, no altitude scaling is needed for the high flying altimeters (T/P and Jason).As for the low flying altimeters, the scaling factor to be applied in the retrieval of TEC alt from the GNSS-derived TEC is the ratio of TEC IRI95 < 800 km to TEC IRI95 < 1400 km .However, [28] demonstrated, by comparing dual-frequency altimeter estimates (from TOPEX, Jason-1 and Envisat) and GIM TEC data, that about 8% of the TEC extends above 1350 km and about 14% above 800 km, thus suggesting the adoption of scaling factors of 0.925 for altimeters around 1350 km, and 0.856 for those around an altitude of 800 km.For CryoSat-2, flying around an altitude of 730 km, this factor is 0.844.The vertical TEC that determines the ionospheric correction to be applied to the altimeter range measurement is known to exhibit large spatial and temporal variations.As illustrated on the global maps of Figure 19, the ionospheric delay presents a latitudinal dependence, varying by a factor of four and with maxima occurring along two bands paralleling the geomagnetic equator.Towards the poles, both the mean and the variation become smaller.
Despite the relative homogeneity of its temporal mean, local ionospheric delays are prone to exhibit large variations, depending on season of the year, time of day and solar activity.The annual variation is mainly due to the changing orientation of the sun with respect to the geomagnetic field, being in general relatively small.However, seasonal variations are very significant in polar regions, where summer levels of TEC can be five times higher than the winter levels.Diurnal TEC variation is quite pronounced, generally peaking around 14:00 local time, with the highest level up to ten times the lowest, occurring around 02:00 during the night.The dependency on solar activity is mostly related to the 11-year cycle in solar radio flux, with levels for high activity periods going up to about five times those of low activity (see Figure 20).For the whole ionosphere, TEC values can range from one TECU during the night time to about 180 TECU in daytime during peak periods of solar activity [28].In a comparison recently performed by [28], the RMS error of TEC models was shown to be more or less proportional to the TEC itself, with the RMS error of 14% to GIM, 18% to NIC09 and 35% to IRI2007.Therefore, from 1998 onwards, whenever dual-frequency measurements are not available or cannot be used to derive reliable ionospheric range corrections, properly interpolated and scaled GIM TEC data provide the more accurate alternative.Prior to that, one has to rely on less accurate climatological models, with increased error for periods of high solar activity, as ionosphere climatologies, tuned to model the lower TEC conditions, were shown to have a tendency to underestimate higher TEC values, IRI2007 significantly more so than NIC09 [28].

Analysis of Issues Regarding the Computation of the Ionospheric Correction over Inland Waters
As the ionosphere is insensitive to coastlines or surface topography, its effect on the propagation of the altimetric radar pulse, unlike those of the troposphere (wet and hydrostatic) does not depend on the altitude of the target, nor is it sensitive to any near-surface atmospheric currents.However, originally designed for ocean applications, the Ku and C (or S) radar altimeter bands have their range measurements differently affected by land.As a consequence, for inland waters applications, the accuracy of the estimation of the ionospheric path delay from the dual-frequency range difference, as detailed in Equation ( 13), is constrained both by the contamination of the returned radar pulse by land within the altimeter footprint and by the ability to perform the necessary along-track smoothing of the difference of the two range measurements (e.g., [73]).
Pioneering inland water altimetry applications (e.g., [1,74]) have relied mostly on climatologies or data derived from DORIS to compensate for the effect of the ionosphere, while others (e.g., [21]) relied on the smoothed dual-frequency altimeter data.With the extension of the time series of altimetric data and the consequent availability of GNSS-derived TEC, the more accurate GIM data became the preferred data source for the ionospheric correction of the range measurement (e.g., [3,[75][76][77]), and the DORIS-derived ionospheric correction is no longer reported on the altimeter data products due to its relatively reduced accuracy.Although most of the authors using altimetry data for inland waters studies do not rely on altimeter dual-frequency-derived ionospheric correction, as part of a calibration study performed over Lake Issykkul, Cretaux et al. [5] compared the use of the altimeter dual-frequency ionospheric correction with that derived from GIM data, for Envisat and Jason-1.Their analysis is based on the comparison of the absolute bias of the altimeters (as deduced from GPS campaigns) when each of the ionospheric corrections available in the respective GDRs is applied.In that case study, GIM was shown to produce lower absolute Ku-band range measurement biases than dual-frequency (46.9 vs. 48.1 cm for Envisat and 5.4 vs. 7.0 cm for Jason-1) for the short, and low solar activity, period tested (September 2004, and September 2005).Studies by [16], as part of their validation exercise of Jason-2 altimeter range data over inland water bodies, show that for large unfrozen lakes, the dual-frequency instrument-based correction presents a much higher variability than that derived from GIM.As part of the same study, the particular case of Lake Ontario is also referred to as a case in which the use of the GIM-derived ionospheric improves the RMS of the estimated lake height (when compared to ground-based gauge data) by ~1.5 cm over the application of dual-frequency ionospheric correction.In their absolute calibration campaign of the Jason-1 and Jason-2 altimeters over Lake Issykkul, Cretaux et al. [6] further compared the use of the dual-frequency ionospheric corrections and those derived from GIM. Considering only altimeter data 10 km off the coastline and further filtering for outliers (differences larger than 5 cm), the authors found an agreement relative to GIM of about 1.1 cm for Jason-1 [7] and 0.4 ± 1.1 cm for Jason-2, over the period of August 2008 to August 2010.Moreover, due to the occurrence of errors of several centimeters in the retrieval of the dual-frequency ionospheric correction, affecting about 10% of Jason-2 data and even more for Jason-1, the GIM correction was selected as the most appropriate for their study.
In order to make a more definitive statement about the preferred ionospheric correction for inland water studies, we compared lake level measurements from Envisat and Jason-1, while applying various ionospheric corrections: JPL-GIM data scaled for the altitude of the satellite orbits as recommended by [28], dual-frequency altimeter data and the same dual-frequency data smoothed by averaging valid lake returns over the entire track crossing the lake.The best way to gauge the relative errors in each of the approaches consists in differencing lake levels at cross-over locations on the satellite ground tracks of the same altimeter, shortly enough spaced in time so that the lake level will not have substantially changed in the meantime.That time lapse was limited to half the length of the repeat cycle: five days in the case of Jason-1; 17.5 days in the case of Envisat.The measured lake level differences on these crossing points, known as crossover height differences, were used as a measure of the quality of the measurements: the lower the RMS crossover height difference, the more precise the measurements.Unfortunately, only the Great Lakes are a large enough body to have a substantive amount of crossovers to do this.An oceanic region with a similar size and latitudinal domain was also selected for comparison, and the results are presented in Table 4 for a two-year period of median solar activity (2003)(2004)).The water level in the Great Lakes varies roughly annually, with an amplitude of about 50 cm.Given the general distribution of the time difference of the crossovers, this contributes about 3 cm to the RMS figures presented for Envisat; less for Jason-1.According to the presented results, the conclusions inferred for inland waters are not, in essence, different from what was already known (and again shown) for ocean applications.Whenever dual-frequency derived ionospheric correction can be retrieved and smoothed, its use yields a slightly better range estimation (lower RMS error) than what is obtained with GIM.However, over smaller lakes, and certainly over rivers, it will not be possible to smooth the dual-frequency correction sufficiently to bring the noise down.In addition, land-contamination of the C-(or S-) band signal is of concern, as well, similar to coastal regions [14].Thus, the impact of smoothing the dual-frequency ionospheric correction over inland waters is less pronounced than over ocean.Only over very sizeable water bodies, like the Great Lakes in the example, smoothing can be expected to be effective enough to bring the crossover RMS error down below that obtained with the GIM-derived correction.Even over the Great Lakes, the differences in the standard deviation of the crossover height differences are small: 0.1 and 0.2 mm for Jason-1 and Envisat, respectively.These results largely confirm what was suggested earlier [5,16], that the ionospheric correction based on JPL GIM was to be preferred for inland waters.However, we like to stress that for large bodies of inland water (200 km and larger, like the Great Lakes), the smoothed dual-frequency ionospheric correction can still provide better results.
Another possibility that could be further investigated is the use of regional GNSS-derived models for the ionospheric correction of inland water altimetry.With the enhanced GNSS coverage in Europe and North America, regional ionosphere maps (RIMs) can be produced at a much finer spatial and temporal resolution in these areas.The RIMs produced since 2004 for Europe (30°-70°N; 15°W-40°E) by the SPECTRE project (Service Products for ionospheric Electron Content and Troposphere Refractivity over Europe) were shown to better resolve small-scale variations of the ionosphere than the JPL GIMs (e.g., [78]).
In conclusion, in agreement with previous studies, the authors recommend using the GIM-derived ionospheric correction over inland water bodies, for dual-and single-frequency altimeters alike.This avoids the possibility of land contamination in the C-or S-band footprint, as well as the limited number of points over which to smooth the dual-frequency correction.For the correction of altimeter measurements prior to the availability of GIM data (1998), the use of NIC09 is recommended.One should be aware that the GIM-derived ionospheric correction provided on most TOPEX, Jason-1 and Jason-2 GDR products is not reduced for the TEC extending above the altitude of these satellites, thus overestimating the ionospheric correction by about 8%.

Discussion and Conclusions
For inland water (IW) studies using satellite radar altimetry, in addition to the atmospheric corrections, the most relevant corrections are those from solid Earth, ocean tide loading and lake tide, if applicable.The dynamic atmospheric correction is not applied, because the lakes/reservoirs are closed systems.The SSB is usually also not applied, because wind effects tend to be averaged out along-track [15], although recent studies made by the authors, outside the scope of this paper, indicate that for regions, such as the Great Lakes, the application of the SSB correction slightly reduces the lake level height anomaly variance at the crossovers.
It has been shown that the tropospheric corrections are not modeled and computed properly and in a consistent way in the various altimetric products, as reported by many authors.If not appropriately corrected, they may be the main source of error for IW applications.In contrast to the ocean domain, IW studies require the modeling of the height dependence of tropospheric corrections.Since most radar altimeter products have been designed for oceanographic purposes, they often fail to provide valid corrections over continental water regions.
While over the ocean, the DTC is one of the most precise range corrections (better than 1 cm), in some of the present altimeter products, it is the correction with the largest errors for IW studies, up to several decimeters.It has been shown that the DTC determined from model grids may possess mainly two types of errors: (1) large biases if it is provided at the sea instead of at the surface level; and (2) interpolation errors with a strong along-track signature, when it is given at the level of the model orography and has been computed from surface pressure grids.
Users must be aware that the DTC derived from the ECMWF operational model is not provided homogeneously on the GDR of the various missions.When provided at sea level, it can be reduced to surface height by applying an appropriate height reduction.However, when provided at the model orography, it is difficult to correct for the interpolation errors that might occur.Therefore, for IW studies, the DTC should be derived by first computing the correction at each altimeter measurement location by interpolation from SLP model grids and then performing a suitable height reduction, preferably using Equation (6), which accounts for the seasonal variations of pressure with temperature.Due to model stability, the use of the DTC from ERA Interim, as provided in RADS, is recommended.
Due to its large spatiotemporal variability, the WTC may be the second largest source of error over IW regions.Being the best correction over the ocean, MWR-derived WTC should also be adopted whenever available, over the largest lakes, such as the Great Lakes and Lake Victoria.In regions possessing GNSS permanent stations, GNSS-derived WTC can be a very valuable and accurate data source, particularly, e.g., for small lakes for which the data from a single station is representative of the local wet path delay conditions.In the absence of any of the previous data types, over small lakes and rivers, model-derived corrections from the ERA Interim, computed at surface height, provide the highest accuracy (1-3 cm).
The WTC from ECMWF operational model provided in the GDR products of the various missions is not given consistently and must be used with caution, particularly on T/P and ERS products.Here, the ECMWF-derived WTC is computed at sea level, and the correction to surface height using Equation (11) does not give enough accuracy for heights above 1000 m [39].When properly given at the model orography or at surface height, for the latest years (since 2004), the ECMWF operational model provides similar or better results compared with ERA Interim, while the NCEP model available in RADS is less accurate than any of the other two ECMWF products.
While over the ocean, smoothed dual-frequency instrument-derived ionospheric correction is the most accurate, the fact that the terrain has different effects on Ku-and C-/S-bands and the difficulty in performing an efficient smoothing over inland water regions make this correction unsuitable for use in most of these regions.Thus, from 1998 onwards, whenever dual-frequency measurements are not available or cannot be used to derive reliable IC, properly interpolated and scaled (as recommended by [28]) GIM TEC data provide the best source for the ionospheric path delay.Prior to that, the climatological model NIC09 can be used, with increased errors for periods of high solar activity.
This paper has addressed the main issues in the atmospheric corrections present in the products of the current altimetric missions, for studies over continental water surfaces, with emphasis on the missions most suitable for these studies.
In spite of the fact that CryoSat-2 does not carry an on-board MWR, the unique characteristics of its data makes them of great interest for inland water studies due to the fact that it operates in three modes [79].For IW studies, the synthetic aperture radar (SAR) mode is of particular relevance, providing measurements with lower noise and improved spatial resolution of 300 m in the along-track direction, over pre-defined regions of interest, including IW regions such as the Amazon basin.
Concerning future missions, an altimeter similar to that of CryoSat-2 is planned for Sentinel-3, currently planned to be operating continuously in SAR mode, while Jason-CS (Jason Continuity of Service) is expected to operate for the first time in the so-called interleaved mode, where the satellite is capable of performing simultaneous and continuous measurements both in the traditional and SAR modes.These missions will increase the capability of acquiring measurements over IW surfaces with increasing accuracy and spatial resolution.A correct handling of all range and geophysical corrections, in particular those from the atmosphere, is crucial in the exploitation of these datasets.
Hydrological studies will reach a new era with the launch of the Surface Water and Ocean Topography (SWOT), a Ka-band wide-swath radar interferometer that will provide sea surface heights (SSH) and terrestrial water heights over a 120 km-wide swath with a gap of approximately 10 km around the nadir track.Over the deep oceans, it will provide SSH with a resolution of 2 m in along-track directions, while in the cross-track direction, the resolution will range from 10 m in the far swath to 60 m in the near swath.Over land, it will produce a water mask able to resolve 100-m rivers and 1-km 2 lakes, wetlands or reservoirs.Associated with this mask, water level elevations with an accuracy of 10 cm and a slope accuracy of 1 cm/km will be available.This high spatial resolution sets new requirements for the range and geophysical corrections.For this purpose, SWOT will carry a new type of radiometer with additional high frequency channels, to allow for a smaller footprint and improved water vapor retrieval in the coastal regions and over inland waters.

Figure 1 .
Figure 1.Location of inland water measurements for a typical Envisat cycle extracted from Radar Altimeter Database System (RADS) using the criteria explained in Section 2.1.2.The color scale represents the water level height above the Earth Gravitational Model 2008 (EGM2008) geoid in meters.
altimeters only: based on difference in the Ku-band and S-or C-band range.Recomputed taking into account range biases in secondary channel.An along-track smoothed version is also provided.Based on difference in the Ku-band and S-or C-band range.Does not take range biases into account.Not available on Envisat since the failure of the S-band on 27 January 2008.

Figure 2 .
Figure 2. (Top) Dry tropospheric correction from three models as present in RADS: ECMWF operational (red), ERA interim (blue) and NCEP (green), in meters, versus surface height above the EGM2008 geoid (m) for Jason-1 (all points for all Phase A cycles, 1 to 259).Similar plots are obtained for Jason-2, Envisat, CryoSat-2 and SARAL.(Bottom) The same as in the top plot for TOPEX/Poseidon (T/P) (all cycles for Phase A, 1 to 364).Similar plots are obtained for all ERS-2 and ERS-1 analyzed cycles.See the text for further details.

Figure 3 .
Figure 3. Dry tropospheric correction (top, in meters) and surface topography from DTU10 (bottom, in meters) versus the point number by ascending time order, for Jason-1 Cycle 62. Since NCEP and ERA corrections are very similar, when the green line cannot be seen, it is behind the blue line.

Figure 4 .
Figure 4. Illustration of dry tropospheric correction (DTC) errors for Jason-1 Pass 222 over Lake Tanganyika.(a) Pass location over ECMWF orography (contours, contour interval = 100 m) and surface pressure (color map, in hPa) of the closest in time ECMWF operational model grid; (b) Pass location over a DEM; (c) DTC (top, m) for Pass 222, Cycle 62; the surface height from ECMWF orography (bottom, in brown, m) and the lake level height above the EGM2008 geoid (as measured by Jason-1, in blue, m).

Figure 5 .
Figure 5. Mean cycle values of the dry tropospheric correction (mean value of the DTC for all selected inland water points of each cycle, ~10 days) from the ECMWF operational model (red), ERA interim (blue) and NCEP (green), in meters, versus the cycle number, for T/P Cycles 1 to 364 (top), for Jason-1 Cycles 1 to 259 (middle) and for Jason-1 Cycles 1 to 259 Pass 222 over lake Tanganyika (bottom).See the main text for further details.

Figure 6 .
Figure 6.The same as in Figure 5, for Envisat Cycles 10 to 93 for GDR data version 1 (top) and version 2.1 (bottom).Cycles 1 to 9 are not shown, since Cycles 1 to 5 are not available in RADS and Cycles 6 to 9 contain a small number of points.

Figure 7 .
Figure 7. DTC (top, in m) and surface topography from DTU10 (bottom, in m) versus the point number by ascending time order, over the Caspian Sea, for SARAL cycles 1 and 2. ECMWF Operational DTC values are those extracted from the IGDRs present in RADS.The blue line represents ERA DTC computed using a mean height of −27 m, while the black line shows the ERA DTC values obtained using the DTU10 topography model in the computation of the height reduction.The difference between the red and blue lines (mean value of 7 mm) reveals that the ECMWF operational DTC was computed using a mean height of 0 m.The DTC from NCEP computed using a mean height of −27 m (not shown) is very similar to the corresponding ERA values (blue curve).

Figure 8 .
Figure 8. Global mean WTC for each 10-day cycle from three models as present in RADS, for all inland water points extracted as described in Section 2.1.2:ECMWF operational (red), ERA Interim (blue) and NCEP (green), for T/P Phase A, Jason-1 Phase A and Jason-2.

Figure 9 .
Figure 9. Mean cycle differences (top) and standard deviation of the differences (bottom) between the WTC from ERA Interim and from the ECMWF operational model (red) and from NCEP (green) for T/P Phase A, Jason-1 Phase A and Jason-2.

Figures 10 -
Figures 10-13 and Tables2 and 3show the comparison between the WTC from three models and from the on-board MWR of each mission for Lake Victoria (Africa) and the Caspian Sea.These are examples of two regions with significantly different conditions for various range corrections, including the WTC.Lake Victoria is one of the world's largest lakes, located in Africa, with a mean height of about 1135 m.Over this lake, the WTC has a weak seasonal signal, with values between 10 cm and 30 cm all over the year (Figure10).Over this lake, all models capture very well the main wet path delay variations revealed by the MWR correction, except for some wetter conditions.Results show that the T/P correction from the ECMWF operational model is too negative by about 10 cm (Figures10 and 11and Table2), since it is provided at sea level instead of lake height.The authors of[5] reported a similar behavior of the ECMWF WTC on the T/P GDR, which also provided the correction at sea level, derived from global 3D grids.From the three models, ERA is the one providing the best results, with a mean difference with respect to MWR, considering the whole analyzed period from 1992 to 2013, of 0.0 cm and a standard deviation of 2.3 cm.The corresponding statistics for the difference between MWR and NCEP are −3.4 cm for the mean and 2.9 cm for the standard deviation.Thus, although in terms of variability, ERA and NCEP seem to be similar over Lake Victoria, NCEP has a bias of about 3 cm with respect to ERA and MWR.

Figure 11 .
Figure11.Mean cycle differences between the MWR-derived WTC and the corresponding values from three models, over Lake Victoria, for T/P Phase A, Jason-1 Phase A and Jason-2.Only measurements over inland water with valid MWR correction were considered.

Figure 13 .Table 2 .
Figure 13.Mean cycle differences between the MWR-derived WTC and the corresponding values from three models, over the Caspian Sea, for T/P Phase A, Jason-1 Phase A and Jason-2.Only measurements over inland water with valid MWR correction were considered.

Figures 16 to 18
Figures 14 and 15 illustrate the impact of the corrections on the determination of water level time series.As discussed in the previous sections, the fact that in the T/P GDRs, the tropospheric corrections from ECMWF are provided at sea level causes a bias (39 cm for Lake Victoria and 26 cm for Lake Tanganyika).The actual errors are not exactly constant, since, e.g., the height dependence of the wet tropospheric correction is a function of the correction itself, Equation(10); therefore, for wetter conditions, the errors are larger than for dryer conditions.However, when computing the mean for, e.g., a 10-day cycle, the errors nearly average out.When deriving time series from various missions, for which the corrections are not consistent, the impact on the water level time series is evident, with a significant impact on the long-term water level trend.This effect is larger for Lake Victoria, since it is at a higher elevation than Lake Tanganyika.Figures 16 to 18 illustrate the impact of different types of tropospheric errors in the determination of mean lake profiles.Three examples are given for Pass 222 over Lake Tanganyika.

Figure 16 .
Figure 16.Mean water level anomalies (m) along T/P Pass 222 over Lake Tanganyika, derived using tropospheric (dry and wet) corrections from two models, for the whole period of T/P Phase A. Water level heights are with respect to the EGM2008 geoid.A mean lake height of 769 m has been removed.Results for NCEP are not shown, since they are very similar to those from ERA Interim.

Figure 17 .
Figure 17.Mean water level anomalies (m) along T/P Pass 222 over Lake Tanganyika, derived using tropospheric (dry and wet) corrections from two models, for the period of J1 Cycles 1-149.Water level heights are with respect to the EGM2008 geoid.A mean lake height of 769 m has been removed.Results for NCEP are not shown, since they are very similar to those from ERA Interim.

Figure 18 .
Figure 18.Mean water level anomalies (m) along T/P Pass 222 over Lake Tanganyika, derived using tropospheric (dry and wet) corrections from two models, for the period of J1 Cycles 150-259.Water level heights are with respect to the EGM2008 geoid.A mean lake height of 769 m has been removed.Results for NCEP are not shown, since they are very similar to those from ERA Interim.

Figure 19 .
Figure 19.Mean (top) and standard deviation (bottom) of the Ku-band ionospheric path delay based on 13 years of JPL global ionosphere map (GIM) TEC maps (1999-2011).

Figure 20 .
Figure 20.Variation of solar 10.7 cm radio flux during the last three solar cycles.The duration of the altimeter missions is shown by rectangles.Outlined rectangles of TOPEX, Jason-1, Jason-2 and Envisat indicate the availability of dual-frequency measurements.The checkered patterns indicate limited data coverage or extended missions.CS2, CryoSat-2; SA, SARAL/AltiKa.Global TEC models are shown at the bottom.

Table 1 .
Detailed description and comparison of the data fields contained in the Geophysical Data Records (GDR) products and RADS for the altimeters, Jason-1, Jason-2 and Envisat.ECMWF, Centre for Medium-Range Weather Forecasts; NCEP, National

Table 4 .
Standard deviations (in mm) of the crossover height differences for Jason-1 and Envisat during the years 2003-2004, using JPL GIM, non-smoothed dual-frequency range measurements and smoothed dual-frequency range measurements for the correction of the ionospheric refraction.The N-W Atlantic Ocean region (52°-40°W, 41°-49°N) is shown for comparison.