SOLS: A lake database to monitor in the Near Real Time water level and storage variations from remote sensing data

An accurate and continuous monitoring of lakes and inland seas is available since 1993 thanks to the satellite altimetry missions (Topex–Poseidon, GFO, ERS-2, Jason-1, Jason-2 and Envisat). Global data processing of these satellites provides temporal and spatial time series of lakes surface height with a decimetre precision on the whole Earth. The response of water level to regional hydrology is particularly marked for lakes and inland seas in semi-arid regions. A lake data centre is under development at by LEGOS (Laboratoire d’Etude en Géophysique et Océanographie Spatiale) in Toulouse, in coordination with the HYDROLARE project (Headed by SHI: State Hydrological Institute of the Russian Academy of Science). It already provides level variations for about 150 lakes and reservoirs, freely available on the web site (HYDROWEB: http://www.LEGOS.obs-mip.fr/soa/hydrologie/HYDROWEB), and surface-volume variations of about 50 big lakes are also calculated through a combination of various satellite images (Modis, Asar, Landsat, Cbers) and radar altimetry. The final objective is to achieve in 2011 a fully operating data centre based on remote sensing technique and controlled by the in situ infrastructure for the Global Terrestrial Network for Lakes (GTN-L) under the supervision of WMO (World Meteorological Organization) and GCOS (Global Climate Observing System).     

GNSS-R open-loop difference phase altimetry: Results from a bridge experiment

The paper explores a method to obtain accurate lake surface heights using measurements of the Global Navigation Satellite System (GNSS) carrier phase reflected from the lake surface. The method is referred to as Global Navigation Satellite System-Reflection (GNSS-R) open-loop difference phase altimetry method. It consists of two key technologies: one is the open-loop tracking method to track the GNSS-R signals, where the direct GNSS signal’s frequency is used as a reference frequency to obtain the carrier phases of the GNSS-R signals; the other key technology is time difference phase altimetry method to invert the lake surface heights using two or more carrier phases of GNSS-R signals received simultaneously. A validation experiment is carried out on the SANYING bridge over GUANTING lake using a GNSS-R receiver developed by the Center for Space Science and Applied Research (CSSAR), processing the data with GNSS-R open-loop difference phase altimetry method. The lake surface height results are consistent with the height results of GPS dual-frequency differential positioning altimetry. The results show that we can achieve centimeter level height in one minute average, by using 11 minutes carrier phase data of three GNSS-R signals received simultaneously.     

Calibration of Envisat radar altimeter over Lake Issykkul

This study presents the results of calibration/validation (C/V) of Envisat satellite radar altimeter over Lake Issykkul located in Kyrgyzstan, which was chosen as a dedicated radar altimetry C/V site in 2004. The objectives are to estimate the absolute altimeter bias of Envisat and its orbit based on cross-over analysis with TOPEX/Poseidon (T/P), Jason-1 and Jason-2 over the ocean. We have used a new method of GPS data processing in a kinematic mode, developed at the Groupe de Recherche de Geodesie Spatiale (GRGS), which allows us to calculate the position of the GPS antenna without needing a GPS reference station. The C/V is conducted using various equipments: a local GPS network, a moving GPS antenna along the satellites tracks over Lake Issykkul, In Situ level gauges and weather stations. The absolute bias obtained for Envisat from field campaigns conducted in 2009 and 2010 is between 62.1 and 63.4 ± 3.7 cm, using the Ice-1 retracking algorithm, and between 46.9 and 51.2 cm with the ocean retracking algorithm. These results differ by about 10 cm from previous studies, principally due to improvement of the C/V procedure. Apart from the new algorithm for GPS data processing and the orbit error reduction, more attention has been paid to the GPS antenna height calculation, and we have reduced the errors induced by seiche over Lake Issykkul. This has been assured using cruise data along the Envisat satellite track at the exact date of the pass of the satellite for the two campaigns. The calculation of the Envisat radar altimeter bias with respect to the GPS levelling is essential to allow the continuity of multi-mission data on the same orbit, with the expected launch of SARAL/Altika mission in 2012. Implications for hydrology in particular, will be to produce long term homogeneous and reliable time series of lake levels worldwide.     

Prospects of Global Navigation Satellite System (GNSS) reflectometry for geodynamic studies

Innovative processing of satellite radar altimetry over solid Earth has been successfully applied for observing geodynamic process of glacial isostatic adjustment over the former Laurentide Ice Sheet in the present-day Hudson Bay land region. In this contribution, a simulation is conducted to study the prospects of the applications of space-/airborne and land-based Global Navigation Satellite System (GNSS) reflectometry to synoptically observe global-scale geodynamic processes with a vertical accuracy of ∼2 mm/yr.     

Atmospheric correction of satellite altimetry observations and sea-level variability in the NE Atlantic

Satellite altimetry provides continuous and spatially regular measurements of the height of the sea surface. Sea level responds to density changes of the water, to mass changes, due to addition or reduction of water mass, and to changes in the atmosphere above it. The present study examines the influence of atmospheric effects on sea-level variability in the North-East Atlantic. The association between the height of the sea surface and the North Atlantic Oscillation (NAO) is investigated by considering different sets of altimetry measurements for which the atmospheric effects have been handled differently. Altimetry data not corrected for atmospheric effects are strongly anti-correlated with the state of the NAO, reflecting the hydrostatic response of sea-level to the NAO pressure dipole. The application of an atmospheric correction to satellite altimetry observations in the NE Atlantic decreases variability of the height time series by more than 70% and reduces the amplitude of the seasonal cycle by ∼5 cm. Altimetry data for which atmospheric effects are removed via an inverse barometer correction show a non-negligible correlation with the NAO index at some locations suggesting further indirect non-hydrostatic influences of the state of the NAO on sea level variability.     

Using altimetry and seafloor pressure data to estimate vertical deformation offshore: Vanuatu case study

Measuring ground deformation underwater is essential for understanding Earth processes at many scales. One important example is subduction zones, which can generate devastating earthquakes and tsunamis, and where the most important deformation signal related to plate locking is usually offshore. We present an improved method for making offshore vertical deformation measurements, that involve combining tide gauge and altimetry data. We present data from two offshore sites located on either side of the plate interface at the New Hebrides subduction zone, where the Australian plate subducts beneath the North Fiji basin. These two sites have been equipped with pressure gauges since 1999, to extend an on-land GPS network across the plate interface. The pressure series measured at both sites show that Wusi Bank, located on the over-riding plate, subsides by 11 ± 4 mm/yr with respect to Sabine Bank, which is located on the down-going plate. By combining water depths derived from the on-bottom pressure data with sea surface heights derived from altimetry data, we determine variations of seafloor heights in a global reference frame. Using altimetry data from TOPEX/Poseidon, Jason-1, Jason-2 and Envisat missions, we find that the vertical motion at Sabine Bank is close to zero and that Wusi Bank subsides by at least 3 mm/yr and probably at most 11 mm/yr.     

Detection of Envisat RA2/ICE-1 retracked radar altimetry bias over the Amazon basin rivers using GPS

Altimetry is now routinely used to monitor stage variations over rivers, including in the Amazon basin. It is desirable for hydrologic studies to be able to combine altimetry from different satellite missions with other hydrogeodesy datasets such as leveled gauges and watershed topography. One requirement is to accurately determine altimetry bias, which could be different for river studies from the altimetry calibrated for deep ocean or lake applications. In this study, we estimate the bias in the Envisat ranges derived from the ICE-1 waveform retracking, which are nowadays widely used in hydrologic applications. As a reference, we use an extensive dataset of altitudes of gauge zeros measured by GPS collocated at the gauges. The thirty-nine gauges are spread along the major tributaries of the Amazon basin. The methodology consists in jointly modeling the vertical bias and spatial and temporal slope variations between altimetry series located upstream and downstream of each gauge. The resulting bias of the Envisat ICE-1 retracked altimetry over rivers is 1.044 ± 0.212 m, revealing a significant departure from other Envisat calibrations or from the Jason-2 ICE-1 calibration.     

Preliminary gravity recovery from CryoSat-2 data in the Baffin Bay

A number of geophysical phenomenons in the open ocean are still unresolved by conventional altimetry, but could be resolved through the potential improvements offered by Synthetic Aperture Radar (SAR), also called Delay-Doppler, altimetry. The SAR altimeter offers the following benefits with respect to conventional satellite altimetry: factor of 20 improvements in the along-track resolution, the along-track footprint length which does not vary with wave height (sea state), and improved precision in sea surface height measurements or sea surface slope measurements.     

Influence of time variable geopotential models on precise orbits of altimetry satellites,global and regional mean sea level trends

During the last decade a significant progress has been reached in the investigation of the gravity field of the Earth. Besides static, also time variable geopotential models have been recently created. In this paper we investigate the impact of the recent time variable geopotential models on altimetry satellite orbits and such altimetry products based on these orbits, as global and regional mean sea level trends. We show that the modeling of time variable gravity improves the orbit solutions, at least for the GRACE period where time variable gravity is sufficiently accurately observed by this mission. Our analysis includes six geopotential models jointly developed by GFZ German Research Centre for Geosciences and Space Geodesy Research Group (CNES/GRGS) Toulouse: the stationary model EIGEN-GL04S, a stationary version of EIGEN-6S (EIGEN-6S_stat), a corrected version of EIGEN-6S and three enhanced versions of EIGEN-6S called EIGEN-6S2, EIGEN-6S2A and EIGEN-6S2B. By “stationary” we mean “containing periodic parameters such as annual and semi-annual variations, but no secular (drift) terms”. We computed precise orbits for the radar altimetry satellites ERS-1, ERS-2, TOPEX/Poseidon, and Envisat over 20 years between 1991 and 2011. The orbit, single-mission and multi-mission altimetry crossover analyses show that the time variable models EIGEN-6S_corrected, EIGEN-6S2 and its two precursors EIGEN-6S2A/B perform notably better than the stationary models for the GRACE period from 2003 onwards. Thus, using EIGEN-6S2 and EIGEN-6S2A/B we have got 3.6% smaller root mean square fits of satellite laser ranging observations for Envisat, as when using EIGEN-GL04S. However, for the pre-GRACE period 1991–2003, the stationary geopotential models EIGEN-GL04S and EIGEN-6S_stat as well as EIGEN-6S2 having no drift terms for degree 3–50 at this time interval perform superior compared to EIGEN-6S_correct and EIGEN-6S2A/B which contain drifts for this period. We found, that the time variable geopotential models have a low (0.1–0.2 mm/yr) impact on our results for the global mean sea level trend. However, we found strong East/West differences up to 3 mm/yr in the regional mean sea level trends when using orbits of all four satellites based on time variable and stationary geopotential models. We show that these differences are related to the relative drifts of the centers-of-origin between the orbit solutions based on the time variable and stationary geopotential models. From the results of our detailed study, we conclude that the final version of the time variable gravity field model EIGEN-6S2 performs best for the four satellites tested. This model provides the most reliable and mission-consistent sea level estimates for the whole time period from 1992 to 2010. This model is of maximum spherical harmonic degree and order 260 and contains time series for drifts as well as annual and semiannual variations of the spherical harmonic coefficients for degree 2–50.     

Characterizing slope correction methods applied to satellite radar altimetry data: A case study around Dome Argus in East Antarctica

Slope correction is important to improve the accuracy of satellite radar elevation measurements by mitigating the slope-induced error (SE), especially over uneven ground surfaces. Although several slope correction methods have been proposed, guidance in the form of stepwise algorithm on how to implement these methods in processing radar altimetric data at the coding level, and the differences among these methods need to be presented and discussed systematically. In this paper, three existing types of slope correction methods—the direct method (DM), intermediate method (IM), and relocation method (RM, further divided into RM1 and RM2)—are described in detail. In addition, their main differences and features for various scientific applications are analyzed. We conduct a systematic experiment with CryoSat-2 Low Resolution Mode (LRM) data in a physically stable area around Dome Argus in East Antarctica, where in-situ measurements were available for comparison. The slope correction is implemented separately using the three methods, with the latest high-accuracy Reference Elevation Model of Antarctica (REMA) as the a-priori topography model. The bias and precision of the slope-corrected CryoSat-2 data results from the RM2 is ?0.18 ± 0.86 m based on the comparison with the field Global Navigation Satellite System (GNSS) data. The results from the RM2 indicate higher precision compared to those from the RM1. According to the correlation analysis of the slope-corrected CryoSat-2 data results (RM1 and RM2), the bias enlarges and the precision becomes worse when the surface slope increases from 0 to 0.85°. After a comprehensively comparative analysis, we find that the results from the RM1 and RM2 are superior in precision (0.93 m and 0.86 m) with respect to the GNSS data. The relatively low precision (1.22 m) from the IM is due to the potential error from the a-priori digital elevation model (DEM). The DM has the lowest precision (2.66 m). Another experiment over rough topography in West Antarctica is carried out for comparison, especially between the RM1 (precision of 15.27 m) and RM2 (precision of 16.25 m). In general, the RM is recommended for the SE elimination among the three methods. Moreover, the RM2 is firstly considered over smooth topography due to the superior performance in bias and precision, while the RM1 is more suggested over the rough topography because of the slightly smaller bias and better precision. The IM relies much on the accuracy of the a-prior DEM and is not usually recommended, because of the strict requirement in the sampling time between the radar altimetry data and the a-priori DEM to avoid any surface change over time.     

A distant look at the cryosphere

Ninety nine per cent of all the fresh water on the surface of the Earth is in the form of ice. Observations from space have revealed more about the ice than about most other parts of the environment because at the dawn of the satellite era, less was known about it. The cryosphere includes all forms of naturally occurring ice but here we review what space science has done for knowledge of glaciers and ice sheets. Whereas in global terms the cryosphere exists as a response to climate, over large areas it controls climate. While imaging spacecraft systems have proved easiest to interpret, microwave sensors with poor spatial resolution are able to distinguish transient and stable surface features that are invisible to the eye. Imaging radars quite effectively describe sea ice, but precision altimetry is the only practicable method for monitoring changes in the total mass of ice on land.     

The contribution of local gravimetric geoid models to the calibration of satellite altimetry data and an outlook of the latest GOCE GGM performance in Gavdos

The use of geoid heights has been one of the available methodologies utilized for the independent calibration/validation of altimeters on-board satellites. This methodology has been employed for long in the Gavdos dedicated cal/val facility (Crete, Greece), where calibration results for the Jason satellites have been estimated, both for ascending and descending passes. The present work gives a detailed overview of the methodology followed in order to estimate a high-resolution and accuracy gravimetric geoid model for the wider Gavdos area, in support of the on-going calibration work. To estimate the geoid model, the well-known remove-compute-restore method is used while residual geoid heights are estimated through least-squares collocation so that associated errors are determined as well. It is found that the estimated formal geoid errors from LSC along passes 018 and 109 of Jason satellites, used for the bias estimation, range between ±0.8–1.6 cm. The so-derived geoid heights are employed in the determination of the Jason-2 altimeter bias for all available cycles (cycles 1-114, spanning the period from July 2008 to August 2011) together with the RioMed DOT model. From the results acquired the Jason-2 bias has been estimated to be +196.1 ± 3.2 mm for pass 109 and +161.9 ± 5.1 mm for pass 018. Within the same frame, the GOCE/GRACE-based geopotential model GOCO02s has been used to estimate the mean dynamic ocean topography and the steady-state circulation in the area around Gavdos. The so-derived DOT model was used to estimate the Jason-2 bias in an effort to evaluate the performance of satellite-only geoid models and investigate whether their spatial resolution and accuracy provides some improvement w.r.t. traditional local gravimetric geoids. From the results acquired with geoid heights from GOCO02s, the estimated Jason-2 bias deviates significantly from that of the local gravimetric model, which can be attributed to a possible mean offset and the low resolution of GOCE-based GGMs. On the other hand, when the newly estimated GOCE-based DOT was employed with geoid heights from the local gravimetric geoid model, the Jason-2 bias has been estimated to be +185.1 ± 3.2 mm for pass 109 and +130.2 ± 5.1 mm for pass 018.     

Ion composition behaviour in low and mid-latitudes during high solar activity

In the present analysis, the mass spectrometer data from the ISS-b satellite, available in the form of contour plots at an average height of 1100 km for every alternate hour, is used. This analysis showed some interesting results in terms of the diurnal variation of the transition heights; at times dominance of He⁺ over the other ions, and the seasonal variations of different constituents. In the development of ion composition models, it is suggested that this type of result should be taken into consideration.     

Semi-automatic determination of the Azores Current axis using satellite altimetry: Application to the study of the current variability during 1995–2006

Satellite altimetry has been widely used to study the variability of the ocean currents such as the Azores Current (AzC) in the North Atlantic. Most analyses are performed over the region that encloses the current, thus being somehow affected by other oceanographic signals, e.g., eddies. In this study, a new approach for extracting the axis of a zonal current solely based on satellite altimetry is presented. This is a semi-automatic procedure that searches for the maximum values of the gradient of absolute dynamic topography (ADT), using the geostrophic velocity as auxiliary information. The advantage of this approach is to allow the analyses to be performed over a buffer centered on the current axis instead of using a wider region. It is here applied to the AzC for the period June 1995–October 2006.     

Improved kHz-SLR tracking techniques and orbit quality analysis for LEO missions

Satellite gravity field missions such as CHAMP, GRACE and GOCE are designed as low Earth orbiting spacecraft (LEO) with orbit heights of about 250–500 km. The challenging mission objectives require a very precise knowledge of the satellite orbit position in space. For these missions precise orbit information is typically provided by GPS satellite-to-satellite tracking (SST) observations supported by satellite laser ranging (SLR).     

一种基于北斗三号系统的GNSS-R海面干涉测高技术

GNSS-R干涉测高技术可用于中尺度海面高度观测，具有空间分辨率高、测量精度高等优势。与传统的GNSS-R本地码测高技术相比，GNSS-R干涉测高技术可以有效提升高度测量精度。虽然GNSS-R干涉测高技术已有一些研究，但是基于北斗三号的干涉测高应用还很少。本文根据GNSS-R干涉测高技术优势，针对北斗三号系统在干涉测高技术上的应用，研发了支持北斗三号的GNSS-R干涉测高接收机并描述了整体架构及实现。利用所研发的接收机进行水面干涉测高试验，首次获取了北斗三号B1和B2干涉测高波形，与传统GPS L1和北斗B1本地码测高波形进行对比。对两种方法计算出的水面高度进行对比，结果显示北斗三号干涉测高精度明显优于GPS L1和北斗B1传统本地码测高精度。      

Stereoscopic observations from meteorological satellites

The capability of making stereoscopic observations of clouds from meteorological satellites is a new basic analysis tool with a broad spectrum of applications. Stereoscopic observations from satellites were first made using the early vidicon tube weather satellites (e.g., Ondrejka and Conover [1]). However, the only high quality meteorological stereoscopy from low orbit has been done from Apollo and Skylab, (e.g., Shenk $\text{et al.}$ [2] and Black [3], [4]). Stereoscopy from geosynchronous satellites was proposed by Shenk [5] and Bristor and Pichel [6] in 1974 which allowed Minzner $\text{et al.}$ [7] to demonstrate the first quantitative cloud height analysis. In 1978 Bryson [8] and desJardins [9] independently developed digital processing techniques to remap stereo images which made possible precision height measurement and spectacular display of stereograms (Hasler $\text{et al.}$ [10], and Hasler [11]). In 1980 the Japanese Geosynchronous Satellite (GMS) and the U.S. GOES-West satellite were synchronized to obtain stereo over the central Pacific as described by Fujita and Dodge [12] and in this paper. Recently the authors have remapped images from a Low Earth Orbiter (LEO) to the coordinate system of a Geosynchronous Earth Orbiter (GEO) and obtained stereoscopic cloud height measurements which promise to have quality comparable to previous all GEO stereo. It has also been determined that the north-south imaging scan rate of some GEOs can be slowed or reversed. Therefore the feasibility of obtaining stereoscopic observations world wide from combinations of operational GEO and LEO satellites has been demonstrated.Stereoscopy from satellites has many advantages over infrared techniques for the observation of cloud structure because it depends only on basic geometric relationships. Digital remapping of GEO and LEO satellite images is imperative for precision stereo height measurement and high quality displays because of the curvature of the earth and the large angular separation of the two satellites. A general solution for accurate height computation depends on precise navigation of the two satellites. Validation of the geosynchronous satellite stereo using high altitude mountain lakes and vertically pointing aircraft lidar leads to a height accuracy estimate of ± 500 m for typical clouds which have been studied. Applications of the satellite stereo include: 1) cloud top and base height measurements, 2) cloud-wind height assignment, 3) vertical motion estimates for convective clouds (Mack $\text{et al.}$ [13], [14]), 4) temperature vs. height measurements when stereo is used together with infrared observations and 5) cloud emissivity measurements when stereo, infrared and temperature sounding are used together (see Szejwach $\text{et al.}$ [15]).When true satellite stereo image pairs are not available, synthetic stereo may be generated. The combination of multispectral satellite data using computer produced stereo image pairs is a dramatic example of synthetic stereoscopic display. The classic case uses the combination of infrared and visible data as first demonstrated by Pichel $\text{et al.}$ [16]. Hasler $\text{et at.}$ [17], Mosher and Young [18] and Lorenz [19], have expanded this concept to display many channels of data from various radiometers as well as real and simulated data fields.A future system of stereoscopic satellites would be comprised of both low orbiters (as suggested by Lorenz and Schmidt [20], [19]) and a global system of geosynchronous satellites. The low earth orbiters would provide stereo coverage day and night and include the poles. An optimum global system of stereoscopic geosynchronous satellites would require international standarization of scan rate and direction, and scan times (synchronization) and resolution of at least 1 km in all imaging channels. A stereoscopic satellite system as suggested here would make an extremely important contribution to the understanding and prediction of the atmosphere.     

The multi-source data fusion global ionospheric modeling software—IonoGim

We introduce a new global ionospheric modeling software—IonoGim, using ground-based GNSS data, the altimetry satellite and LEO (Low Earth Orbit) occultation data to establish the global ionospheric model. The software is programmed by C++ with fast computing speed and highly automatic degree, it is especially suitable for automatic ionosphere modeling. The global ionospheric model and DCBs obtained from IonoGim were compared with the CODE (Center for Orbit Determination in Europe) to verify its accuracy and reliability. The results show that IonoGim and CODE have good agreement with small difference, indicating that IonoGim owns high accuracy and reliability, and can be fully applicable for high-precision ionospheric research. In addition, through comparison between only using ground-based GNSS observations and multi-source data model, it can be demonstrated that the space-based ionospheric data effectively improve the model precision in marine areas where the ground-based GNSS tracking station lacks.     

利用GNSS反射信号载波测量湖面高度变化

给出一种利用导航卫星反射信号的载波相位测量光滑湖面高度变化的方法.首先建立载波相位观测方程,通过在直射与反射信号的观测量之间求单差,再在历元间求双差的方法,消除了大量误差并得到了直射与反射信号的几何路径差;然后根据各颗卫星的观测值和高度角等信息计算出湖面高度变化,并把各卫星得到的湖面高度度变化值进行加权平均以提高测量系统的精度和可靠性.仿真验证了方法在-15dB信噪比下达到了毫米级的测量精度.     

基于地面移动通信基站的差分气压测高方法

为了解决现有的卫星定位系统高程定位精度较低的问题,提出利用已发展完备的地面移动通信基站作为气压差分测量基准点的差分气压测高的方法.该方法以用户接收机附近的分布密集的地面移动通信基站为气压校正基准点,并利用移动通信基站的传输链路把相关的测量数据传送给用户,用户接收机便可以用移动通信基站观测点的高程值、气压值以及温度值,结合接收机测点处的气压值和温度值,利用高程与气压值之间的对应关系,得到用户接收机处的绝对高程值.实际试验结果表明:使用本方法后,室内高程定位精度约为1 m,室外的高程定位精度优于1 m.测高分辨率约为0.2 m,可以使GNSS(Global Navigation Satellite System)的三维位置定位精度提高约60%.本方法不仅可以弥补GNSS高程精度差的缺陷,作为高程约束参与测量方程解算还会提高水平测量精度,而且还克服了传统的气压测高精度依赖于气象站,而大多数用户接收机离气象站的距离较远,造成气压测高精度偏低的局限性.所以本方法能辅助增强卫星导航系统,提高三维定位精度,用于室内定位可以解决室内楼层的分辨问题.