利用GPS三频数据计算电离层垂直总电子含量变化率

研究了通过采样间隔为1 s的GPS三频载波相位观测数据计算1 Hz电离层垂直总电子含量变化率(RVTEC)的方法,推导了直接采用双频载波相位观测量计算RVTEC的公式,对一般空间环境与电离层活动剧烈两种条件下L1L2、L1L5、L2L5三种载波相位组合计算的RVTEC进行了一致性分析,给出了通过三频载波相位观测数据计算电离层RVTEC的修正方法,比较了X射线太阳耀斑期间RVTEC与由传统方法计算的总电子含量变化率(ROT)响应,在双星串飞编队测高模式下对计算方法进行了应用.结果表明,L1L2、L1L5两种组合之间的一致性较好,由L1L2、L1L5组合计算的RVTEC的中误差约为0.004 TECu/s,RVTEC较ROT更能体现TEC变化的细节,在双星串飞编队测高模式下海面高梯度的计算中,电离层延迟之差可以忽略.     

Improving gravity anomalies over China marginal sea from retracked Geosat and ERS-1 data

The quality of altimeter data and ocean tide model is critical to the recovery of coastal gravity anomalies. In this contribution, three retracking methods (threshold, improved threshold and Beta-5) are investigated with the aim of improving the altimeter data over a shallow water area. Comparison indicates that the improved threshold is the best retracking method over China Sea. Two ocean tide models, NAO99b and CSR4.0, are analyzed. Results show that different tide models used in the processing of altimeter data may result in differences more than 10 mGal in recovered coastal gravity anomalies. Also, NAO99b is more suitable than CSR4.0 over the shallow water area of China Sea. Finally, gravity anomalies over China Sea are calculated from retracked Geosat/GM and ERS-1/GM data by least squares collocation. Comparison with shipborne gravimetry data demonstrates that gravity anomalies from retracked data are significantly superior to those from non-retracked data. Our results have the same order as the other two altimeter-derived gravity models: Sandwell&Smith(V16) and DNSC08.     

基于星载激光测高估计南极冰盖高程变化

利用NSIDC公布的ICESat运行期间的19个任务期的观测数据, 采用重复轨迹分析方法估计近年来南极冰盖高程的时空变化。结果表明,南极大陆冰盖整体上呈现消融趋势,基于不同的GIA模型和ICESat激光测高数据的南极大陆冰盖高程变化的趋势约为 -1.17~-1 cm/a。     

四种南极数字高程模型的精度比较与分析

南极数字高程模型（DEM）是南极冰盖研究的基础数据,目前国际通用的全南极DEM数据主要有JLB97 DEM、RAMPv2 DEM、ICESat DEM以及Bamber 1km DEM. 利用DEM对DEM验证的方式对四种DEM的精度进行比较分析. 结果表明：Bamber 1km DEM和ICESat DEM之间的高程差异最小,平均高程差小于1.8 m,二者均有较高的可靠性. RAMPv2 DEM与Bamber 1km DEM的高程差大于1.9 m,在81.5° S以南和坡度较大的区域,高程差异更为明显,高程可靠性较低. JLB97 DEM与上述三种DEM的偏差超过10 m,高程可靠性最低.     

卫星测高混合边值问题的球谐级数解法

研究了球界面下卫星测高问题的解法,利用有限逼近方法得到了下列结论：若陆地部分是球冠,则卫星测高问题的解可以转换成关于球谐级数位系数的线性方程组,并且位系数的阶和次是以分离形式出现的,从而确保该解法具有实用意义. 利用重力场360阶模型进行模拟计算的结果表明：该解法得到的位系数的相对精度达到了10－１１. 同时证明了常用的Stokes问题、Dirichlet问题、Neumann问题可以看成卫星测高问题的特殊情况.     

Jason-1 Precision Orbit Determination by Combining SLR and DORIS with GPS Tracking Data

With the implementation of the Jason-1 satellite altimeter mission, the goal of reaching the 1-cm level in orbit accuracy was set. To support the Precision Orbit Determination (POD) requirements, the Jason-1 spacecraft carries receivers for DORIS (Doppler Orbitography and Radiopositioning Integrated by Satellite) and GPS (Global Positioning System), as well as a retroreflector for SLR (Satellite Laser Ranging). The overall orbit accuracy for Jason will depend on the quality and the relative weighting of the available tracking data. In this study, the relative importance of the SLR, DORIS, and GPS tracking data is assessed along with the most effective parameterization for accounting for the unmodeled accelerations through the application of empirical accelerations. The optimal relative weighting for each type of tracking data was examined. It is demonstrated that GPS tracking alone is capable of supporting a radial orbit accuracy for Jason-1 at the 1-cm level, and that including SLR tracking provides additional benefits. It is also shown that the GRACE (Gravity Recovery and Climate Experiment) gravity model GGM01S provides a significant improvement in the orbit accuracy and reduction in the level of geographically correlated orbit errors.     

Impact of Ocean Mean Dynamic Topography on Satellite Data Assimilation

The response of an eddy-permitting ocean model to changes imposed by the use of different mean dynamic topographies (MDT) is analyzed in a multivariate assimilation context, allowing the evaluation of this impact, not only on the surface circulation, but also on the interior ocean representation. The assimilation scheme is a reduced-order sequential Kalman filter (SEEK). In a first set of experiments, high resolution sea surface temperature, along-track sea surface height and sea surface salinity from climatology are assimilated into a 1/3° resolution North and Tropical Atlantic version of the HYCOM model. In a second experiment, in situ profile data are assimilated in addition to the surface measurements.

The first set of experiments illustrates that important differences in the representation of the horizontal model circulation pattern are related to differences in the MDT used. The objective of assimilation is to improve the representation of the 3D ocean state. However, the imperfect representation of the mean dynamic topography appears to be an important limiting factor with regard to the degree of realism obtained in the simulated flow.

Vertical temperature and salinity profiles are key observations to drive a general circulation ocean model toward a more realistic state. The second set of experiments shows that assimilating them in addition to sea surface measurements is a far from trivial exercise. A specific difficulty is due to inconsistencies between the dynamic topography diagnosed from in situ observations and that diagnosed from sea surface height. These two fields obtained from different data sources do not contain exactly the same information. In order to overcome this difficulty, a strategy is proposed and validated.     

Accuracy Assessment of the TOPEX/Poseidon Ionosphere Measurements

We conducted an assessment of the TOPEX dual-frequency nadir ionosphere observations in the TOPEX/Poseidon (T/P) GDR by comparing TOPEX with the Center for Orbit Determination in Europe (CODE) Global Ionosphere Map (GIM), the climatological model IRI2001, and the DORIS (onboard T/P) relative ionosphere delays. We investigated the TOPEX (TOPEX Side A and TOPEX Side B altimeters, TSA and TSB, respectively) ionosphere observations for the time period 1995–2001, covering periods of low, intermediate, and high solar activity. Here, we use absolute path delays (at Ku-band frequency of the TOPEX altimeter and with positive signs) rather than Total Electron Content (TEC). We found significant biases between GIM and TOPEX (GIM–TOPEX) nadir ionosphere path delays: ?8.1 ± 0.4 {mm} formal uncertainties and equivalent to 3.7 TECu) and ?9.0 ± 0.7 {mm} (4.1 TECu) for TSA and TSB, respectively, indicating that the TOPEX path delay is longer (or with higher TECu) than GIM. The estimated relative biases vary with latitude and with daytime or nighttime passes. The estimated biases in the path delays (DORIS–TOPEX) are: ?10.9 ± 0.4 {mm} (5.0 TECu) and ?14.8 ± 0.6 {mm} (6.7 TECu), for TSA and TSB, respectively. There is a distinct jump of the DORIS path delays (?3.9 ± 0.7 {mm}, TSA delays longer than TSB delays) at the TSB altimeter switch in February 1999, presumably due to inconsistent DORIS processing. The origin of the bias between GIM (GPS, L-band) and TOPEX (radar altimeter, Ku-band) is currently unknown and warrants further investigation. Finally, the estimated drift rates between GIM and TSA, DORIS and TSA ionosphere path delays for the 6-year study span are ?0.4 mm/yr and ?0.8 mm/yr, respectively, providing a possible error bound for the TOPEX/Poseidon sea level observations during periods of low and intermediate solar activity.     

Towards a Seamless Transition from TOPEX/Poseidon to Jason-1

The Jason-1 verification phase has proven to be a unique and successful calibration experiment to quantify the agreement with its predecessor TOPEX/Poseidon. Although both missions have met prescribed error budgets, comparison of the mean and time-varying sea surface height profiles from near simultaneous observations derived from the missions' Geophysical Data Records exhibit significant basin scale differences. Several suspected sources causing this disagreement are identified and improved upon, including (a) replacement of TOPEX and Jason project POE with enhanced orbits computed at GSFC within a consistent ITRF2000 terrestrial reference frame, (b) application of waveform retracking corrections to TOPEX significant wave height and sea surface heights, (c) resultant improved efficacy of the TOPEX sea state bias estimation from the value added sea surface height, and (d) estimation of Jason-1 sea state bias employing dual TOPEX/Jason crossover and collinear sea surface height residuals unique to the validation mission. The resultant mean sea surface height comparison shows improved agreement at better than 60 percent level of variance reduction with a standard deviation less then 0.5 cm.