Construction of lunar DEMs based on reflectance modelling

Existing lunar DEMs obtained based on laser altimetry or photogrammetric image analysis are characterised by high large-scale accuracies while their lateral resolution is strongly limited by noise or interpolation artifacts. In contrast, image-based photometric surface reconstruction approaches reveal small-scale surface detail but become inaccurate on large spatial scales. The framework proposed in this study therefore combines photometric image information of high lateral resolution and DEM data of comparably low lateral resolution in order to obtain DEMs of high lateral resolution which are also accurate on large spatial scales. Our first approach combines an extended photoclinometry scheme and a shape from shading based method. A novel variational surface reconstruction method further increases the lateral resolution of the DEM such that it reaches that of the underlying images. We employ the Hapke IMSA and AMSA reflectance models with two different formulations of the single-particle scattering function, such that the single-scattering albedo of the surface particles and optionally the asymmetry parameter of the single-particle scattering function can be estimated pixel-wise. As our DEM construction methods require co-registered images, an illumination-independent image registration scheme is developed. An evaluation of our framework based on synthetic image data yields an average elevation accuracy of the constructed DEMs of better than 20 m as long as the correct reflectance model is assumed. When comparing our DEMs to LOLA single track data, absolute elevation accuracies around 30 m are obtained for test regions that cover an elevation range of several thousands of metres. The proposed illumination-independent image registration method yields subpixel accuracy even in the presence of 3D perspective distortions. The pixel-wise reflectance parameters estimated simultaneously with the DEM reflect compositional contrasts between different surface units. Specifically, the detected variations of the parameter of the single-particle scattering function indicate small-scale variations of the regolith particle size, possibly as a result of differences in soil maturity.     

Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

Inter-agency comparison of TanDEM-X baseline solutions

TanDEM-X (TerraSAR-X add-on for Digital Elevation Measurement) is the first Synthetic Aperture Radar (SAR) mission using close formation flying for bistatic SAR interferometry. The primary goal of the mission is to generate a global digital elevation model (DEM) with 2 m height precision and 10 m ground resolution from the configurable SAR interferometer with space baselines of a few hundred meters. As a key mission requirement for the interferometric SAR processing, the relative position, or baseline vector, of the two satellites must be determined with an accuracy of 1 mm (1D RMS) from GPS measurements collected by the onboard receivers. The operational baseline products for the TanDEM-X mission are routinely generated by the German Research Center for Geosciences (GFZ) and the German Space Operations Center (DLR/GSOC) using different software packages (EPOS/BSW, GHOST) and analysis strategies. For a further independent performance assessment, TanDEM-X baseline solutions are generated at the Astronomical Institute of the University of Bern (AIUB) on a best effort basis using the Bernese Software (BSW).     

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.     

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.     

Coastal sea level changes in Europe from GPS,tide gauge,satellite altimetry and GRACE, 1993–2011

Sea level changes are threatening the human living environments, particularly along the European Coasts with highly dense population. In this paper, coastal sea level changes in western and southern Europe are investigated for the period 1993–2011 using Global Positioning System (GPS), Tide Gauge (TG), Satellite Altimetry (SA), Gravity Recovery and Climate Experiment (GRACE) and geophysical models. The mean secular trend is 2.26 ± 0.52 mm/y from satellite altimetry, 2.43 ± 0.61 mm/y from TG+GPS and 1.99 ± 0.67 mm/y from GRACE mass plus steric components, which have a remarkably good agreement. For the seasonal variations, annual amplitudes of satellite altimetry and TG+GPS results are almost similar, while GRACE Mass+Steric results are a little smaller. The annual phases agree remarkably well for three independent techniques. The annual cycle is mainly driven by the steric contributions, while the annual phases of non-steric (mass component) sea level changes are almost a half year later than the steric sea level changes.     

DORIS and GPS monitoring of the Gavdos calibration site in Crete

Due to its specific geographical location as well as its geodetic equipment (DORIS, GNSS, microwave transponder and tide gauges), the Gavdos station in Crete, Greece is one of the very few sites around the world used for satellite altimetry calibration. To investigate the quality of the Gavdos geodetic coordinates and velocities, we analyzed and compared here DORIS and GPS-derived results obtained during several years of observations. The DORIS solution is the latest ignwd11 solution at IGN, expressed in ITRF2008, while the GPS solution was obtained using the GAMIT software package. Current results show that 1–2 mm/yr agreement can be obtained for 3-D velocity, showing a good agreement with current geophysical models. In particular, the agreement obtained for the vertical velocity is around 0.3–0.4 mm/yr, depending on the terrestrial reference frame. As a by-product of these geodetic GPS and DORIS results, Zenith Tropospheric Delays (ZTDs) estimations were also compared in 2010 between these two techniques, and compared to ECMWF values, showing a 6.6 mm agreement in dispersion without any significant difference between GPS and DORIS (with a 97.6% correlation), but with a 13–14 mm agreement in dispersion when comparing to ECMWF model (with only about 90% correlation for both techniques). These tropospheric delay estimations could also provide an external calibration of the tropospheric correction used for the geophysical data of satellite altimetry missions.     

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.     

Study of long path VLF signal propagation characteristics as observed from Indian Antarctic station,Maitri

To examine the quality and propagation characteristics of the Very Low Frequency (VLF) radio waves in a very long propagation path, Indian Centre for Space Physics, Kolkata, participated in the 27th Indian scientific expedition to Antarctica during 2007–2008. One Stanford University made AWESOME VLF receiving system was installed at the Indian Antarctic station Maitri and about five weeks of data were recorded successfully from the Indian transmitter VTX and several other transmitting stations worldwide. The quality of the signal from the VTX transmitter was found to be very good, consistent and highly stable in day and night. The signal shows the evidences of the presence of the 24 h solar radiation in the Antarctic region during local summer. Here we report the both narrow band and broadband VLF observations from this site. The diurnal variations of VTX signal (18.2 kHz) are presented systematically for Antarctica path and also compared the same with the variations for a short propagation path (VTX-Kolkata). We compute the spatial distribution of the VTX signal along the VTX-Antarctica path using the most well-known LWPC model for an all-day and all-night propagation conditions. The calculated signal amplitudes corresponding to those conditions relatively corroborate the observations. We also present the attenuation rate of the dominant waveguide modes corresponding to those propagation conditions where the effects of the Antarctic polar ice on the attenuation of different propagating waveguide modes are visible.     

Towards development of a consistent orbit series for TOPEX,Jason-1, and Jason-2

The TOPEX/Poseidon, Jason-1 and Jason-2 set of altimeter data now provide a time series of synoptic observations of the ocean that span nearly 17 years from the launch of TOPEX in 1992. The analysis of the altimeter data including the use of altimetry to monitor the global change in mean sea level requires a stable, accurate, and consistent orbit reference over the entire time span. In this paper, we describe the recomputation of a time series of orbits that rely on a consistent set of reference frames and geophysical models. The recomputed orbits adhere to the IERS 2003 standards for ocean and earth tides, use updates to the ITRF2005 reference frame for both the SLR and DORIS stations, apply GRACE-derived models for modeling of the static and time-variable gravity, implement the University College London (UCL) radiation pressure model for Jason-1, use improved troposphere modeling for the DORIS data, and apply the GOT4.7 ocean tide model for both dynamical ocean tide modeling and for ocean loading. The new TOPEX orbits have a mean SLR fit of 1.79 cm compared to 2.21 cm for the MGDR-B orbits. These new TOPEX orbits agree radially with independent SLR/crossover orbits at 0.70 cm RMS, and the orbit accuracy is estimated at 1.5–2.0 cm RMS over the entire TOPEX time series. The recomputed Jason-1 orbits agree radially with the Jason-1 GDR-C orbits at 1.08 cm RMS. The GSFC SLR/DORIS dynamic and reduced-dynamic orbits for Jason-2 agree radially with independent orbits from the CNES and JPL at 0.70–1.06 cm RMS. Applying these new orbits, and using the latest altimeter corrections for TOPEX, Jason-1, and Jason-2 from September 1992 to May 2009, we find a global rate in mean sea level of 3.0 ± 0.4 mm/yr.     

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.     

Searching for life in extreme environments relevant to Jovian’s Europa: Lessons from subglacial ice studies at Lake Vostok (East Antarctica)

The objective was to estimate the genuine microbial content of ice samples from refrozen water (accretion ice) from the subglacial Lake Vostok (Antarctica) buried beneath the 4-km thick East Antarctic ice sheet. The samples were extracted by heavy deep ice drilling from 3659 m below the surface. High pressure, a low carbon and chemical content, isolation, complete darkness and the probable excess of oxygen in water for millions of years characterize this extreme environment. A decontamination protocol was first applied to samples selected for the absence of cracks to remove the outer part contaminated by handling and drilling fluid. Preliminary indications showed the accretion ice samples to be almost gas free with a low impurity content. Flow cytometry showed the very low unevenly distributed biomass while repeated microscopic observations were unsuccessful.     

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

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

Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.     

Estimation and analysis of GPS receiver differential code biases using KGN in Korean Peninsula

The total electron content (TEC) estimation by the Global Positioning System (GPS) can be seriously affected by the differential code biases (DCB), referred to as inter-frequency biases (IFB), of the satellite and receiver so that an accuracy of GPS–TEC value is dependent on the error of DCBs estimation. In this paper, we proposed the singular value decomposition (SVD) method to estimate the DCB of GPS satellites and receivers using the Korean GPS network (KGN) in South Korea. The receiver DCBs of about 49 GPS reference stations in KGN were determined for the accurate estimation of the regional ionospheric TEC. They obtained from the daily solution have large biases ranging from +5 to +27 ns for geomagnetic quiet days. The receiver DCB of SUWN reference station was compared with the estimates of IGS and JPL global ionosphere map (GIM). The results have shown comparatively good agreement at the level within 0.2 ns. After correction of receiver DCBs and knowing the satellite DCBs, the comparison between the behavior of the estimated TEC and that of GIMs was performed for consecutive three days. We showed that there is a good agreement between KASI model and GIMs.     

The use of ionospheric tomography and elevation masks to reduce the overall error in single-frequency GPS timing applications

Signals from Global Positioning System (GPS) satellites at the horizon or at low elevations are often excluded from a GPS solution because they experience considerable ionospheric delays and multipath effects. Their exclusion can degrade the overall satellite geometry for the calculations, resulting in greater errors; an effect known as the Dilution of Precision (DOP). In contrast, signals from high elevation satellites experience less ionospheric delays and multipath effects. The aim is to find a balance in the choice of elevation mask, to reduce the propagation delays and multipath whilst maintaining good satellite geometry, and to use tomography to correct for the ionosphere and thus improve single-frequency GPS timing accuracy. GPS data, collected from a global network of dual-frequency GPS receivers, have been used to produce four GPS timing solutions, each with a different ionospheric compensation technique. One solution uses a 4D tomographic algorithm, Multi-Instrument Data Analysis System (MIDAS), to compensate for the ionospheric delay. Maps of ionospheric electron density are produced and used to correct the single-frequency pseudorange observations. This method is compared to a dual-frequency solution and two other single-frequency solutions: one does not include any ionospheric compensation and the other uses the broadcast Klobuchar model. Data from the solar maximum year 2002 and October 2003 have been investigated to display results when the ionospheric delays are large and variable. The study focuses on Europe and results are produced for the chosen test site, VILL (Villafranca, Spain). The effects of excluding all of the GPS satellites below various elevation masks, ranging from 5° to 40°, on timing solutions for fixed (static) and mobile (moving) situations are presented. The greatest timing accuracies when using the fixed GPS receiver technique are obtained by using a 40° mask, rather than a 5° mask. The mobile GPS timing solutions are most accurate when satellites at lower elevations continue to be included: using a mask between 10° and 20°. MIDAS offers the most accurate and least variable single-frequency timing solution and accuracies to within 10 ns are achieved for fixed GPS receiver situations. Future improvements are anticipated by combining both GPS and Galileo data towards computing a timing solution.     

Mapping mean lake surface from satellite altimetry and GPS kinematic surveys

Lake water height is a key variable in water cycle and climate change studies, which is achievable using satellite altimetry constellation. A method based on data processing of altimetry from several satellites has been developed to interpolate mean lake surface (MLS) over a set of 22 big lakes distributed on the Earth. It has been applied on nadir radar altimeters in Low Resolution Mode (LRM: Jason-3, Saral/AltiKa, CryoSat-2) in Synthetic Aperture Radar (SAR) mode (Sentinel-3A), and in SAR interferometric (SARin) mode (CryoSat-2), and on laser altimetry (ICESat). Validation of the method has been performed using a set of kinematic GPS height profiles from 18 field campaigns over the lake Issykkul, by comparison of altimetry’s height at crossover points for the other lakes and using the laser altimetry on ICESat-2 mission. The precision reached ranges from 3 to 7 cm RMS (Root Mean Square) depending on the lakes. Currently, lake water level inferred from satellite altimetry is provided with respect to an ellipsoid. Ellipsoidal heights are converted into orthométric heights using geoid models interpolated along the satellite tracks. These global geoid models were inferred from geodetic satellite missions coupled with absolute and regional anomaly gravity data sets spread over the Earth. However, the spatial resolution of the current geoid models does not allow capturing short wavelength undulations that may reach decimeters in mountaineering regions or for rift lakes (Baikal, Issykkul, Malawi, Tanganika). We interpolate in this work the geoid height anomalies with three recent geoid models, the EGM2008, XGM2016 and EIGEN-6C4d, and compare them with the Mean Surface of 22 lakes calculated using satellite altimetry. Assuming that MLS mimics the local undulations of the geoid, our study shows that over a large set of lakes (in East Africa, Andean mountain and Central Asia), short wavelength undulations of the geoid in poorly sampled areas can be derived using satellite altimetry. The models used in this study present very similar geographical patterns when compared to MLS. The precision of the models largely depends on the location of the lakes and is about 18 cm, in average over the Earth. MLS can serve as a validation dataset for any future geoid model. It will also be useful for validation of the future mission SWOT (Surface Water and Ocean Topography) which will measure and map water heights over the lakes with a high horizontal resolution of 250 by 250 m.     

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

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

The impact of microwave absorber and radome geometries on GNSS measurements of station coordinates and atmospheric water vapour

We have used microwave absorbing material in different geometries around ground-based Global Navigation Satellite System (GNSS) antennas in order to mitigate multipath effects on the estimates of station coordinates and atmospheric water vapour. The influence of a hemispheric radome – of the same type as in the Swedish GPS network SWEPOS – was also investigated. Two GNSS stations at the Onsala Space Observatory were used forming a 12 m baseline. GPS data from October 2008 to November 2009 were analyzed by the GIPSY/OASIS II software using the Precise Point Positioning (PPP) processing strategy for five different elevation cutoff angles from 5° to 25°. We found that the use of the absorbing material decreases the offset in the estimated vertical component of the baseline from ∼27 mm to ∼4 mm when the elevation cutoff angle varies from 5° to 20°. The horizontal components are much less affected. The corresponding offset in the estimates of the atmospheric Integrated Water Vapour (IWV) decreases from ∼1.6 kg/m² to ∼0.3 kg/m². Changes less than 5 mm in the offsets in the vertical component of the baseline are seen for all five elevation cutoff angle solutions when the antenna was covered by a hemispheric radome. Using the radome affects the IWV estimates less than 0.4 kg/m² for all different solutions. IWV comparisons between a Water Vapour Radiometer (WVR) and the GPS data give consistent results.