基于Grace月重力场模型的稳态地球重力场模型分析

基于106个月的Grace（gravity recovery and climate experiment）月重力场模型（120阶次）,消除了月重力场的月、季度及年度变化,得到了稳态的地球重力场模型（Grace_sta）。在2～120阶次之间,Grace_sta与已有高阶重力场EGM2008及EGM96三个模型的阶方差是一致的。在2～100阶次之间,Grace_sta模型误差阶方差要小于EGM2008与EGM96误差阶方差。在全球范围内,Grace_sta重力场的大地水准面与EGM2008相应阶次的大地水准面标准差约为3cm,与EGM96模型大地水准面差异则高达52cm。结果表明,Grace_sta足可以取代EGM2008重力场模型2～100阶次的低阶部分,新得到的稳态重力场模型可为海面地形分析提供了可靠的参考场。     

基于Stokes边值理论的大地水准面计算模型改化及分析检验

针对源自经典Stokes边值理论的大地水准面计算模型适用性问题,提出了3种实用化的Stokes积分模型分阶段改化方法,分析讨论了数据截断误差和观测噪声对大地水准面计算结果的影响,设计了改化计算模型的数值检验方案,使用超高阶位模型EGM2008作为数值模拟标准场,对3种改化计算模型进行了数值计算实验和精度检核,同时开展了数据观测噪声影响评估检验,得出了一些有参考和实际应用价值的研究结论。在一定条件下,使用改化Stokes方法计算大地水准面可获得1 cm的内符合度。     

卫星测高技术在监测海啸现象中的研究及应用

利用卫星测高技术,通过共线平差和沿轨、垂直轨迹测高大地水准面梯度的计算,得到了海啸发生2个小时后的海面波高和海面梯度异常现象,为利用卫星测高技术监测海啸、模拟海啸现象提供了依据。     

多源数据推求的西太平洋区域海面动力地形比较分析

简述利用空间大地测量观测数据和海洋水文数据推求海面动力地形的方法。基于EGM96重力场模型和卫星重力恢复的重力场模型GL04C,联合卫星测高平均海面高模型分别推算西太平洋海域的平均海面动力地形,并与根据海洋水文数据推算之结果进行比较分析。结果表明:卫星重力场模型GL04C更好地表现了海面地形的细节特征。卫星重力和卫星测高的联合应用将成为确定海面动力地形的有效途径之一。     

Comparison of Satellite Altimetric Gravity and Global Geopotential Models with Shipborne Gravity in the Red Sea

The determination of high-resolution geoid for marine regions requires the integration of gravity data provided by different sources, e.g. global geopotential models, satellite altimetry, and shipborne gravimetric observations. Shipborne gravity data, acquired over a long time, comprises the short-wavelengths gravitation signal. This paper aims to produce a consistent gravity field over the Red Sea region to be used for geoid modelling. Both, the leave-one-out cross-validation and Kriging prediction techniques were chosen to ensure that the observed shipborne gravity data are consistent as well as free of gross-errors. A confidence level equivalent to 95.4% was decided to filter the observed shipborne data, while the cross-validation algorithm was repeatedly applied until the standard deviation of the residuals between the observed and estimated values are less than 1.5 mGal, which led to the elimination of about 17.7% of the shipborne gravity dataset. A comparison between the shipborne gravity data with DTU13 and SSv23.1 satellite altimetry-derived gravity models is done and reported. The corresponding results revealed that altimetry models almost have identical data content when compared one another, where the DTU13 gave better results with a mean and standard deviation of ?2.40 and 8.71 mGal, respectively. A statistical comparison has been made between different global geopotential models (GGMs) and shipborne gravity data. The Spectral Enhancement Method was applied to overcome the existing spectral gap between the GGMs and shipborne gravity data. EGM2008 manifested the best results with differences characterised with a mean of 1.35 mGal and a standard deviation of 11.11 mGal. Finally, the least-squares collocation (LSC) was implemented to combine the shipborne gravity data with DTU13 in order to create a unique and consistent gravity field over the Red Sea with no data voids. The combined data were independently tested using a total number of 95 randomly chosen shipborne gravity stations. The comparison between the extracted shipborne gravity data and DTU13 altimetry anomalies before and after applying the LSC revealed that a significant improvement is procurable from the combined dataset, in which the mean and standard deviation of the differences dropped from ?3.60 and 9.31 mGal to ?0.39 and 2.04 mGal, respectively.     

High Precision Geoid Models HKGEOID-2000 for Hong Kong and SZGEOID-2000 for Shenzhen, China

High precision geoid models HKGEOID-2000 for Hong Kong and SZGEOID-2000 for Shenzhen, China, have been developed with a hybrid approach of so-called sequential processing, using high precision GPS/leveling data, land and sea gravity anomalies, and digital terrain models. These two local geoid models have the same 1-km resolution. The estimated accuracy (external accuracy) is better than 1.7 cm for HKGEOID-2000 and 1.4 cm for SZGEOID-2000. Some common areas are covered by HKGEOID-2000 and SZGEOID-2000. So these two geoid models, along with high quality GPS/leveling data collected on the overlapping areas, can be used to detect the systematic bias between HKGEOID-2000 and SZGEOID-2000, as well as the difference between Hong Kong Principal Datum and 1956 yellow sea height datum of China, yielding RMS errors of 1.011 m and 1,003 m, respectively. Moreover, HKGEOID-2000, along with GPS ellipsoidal heights, is employed to determine the errors of the “orthometric heights” from purely trigonometric heighting, yielding an RMS error of 0.102 m. The combination of SZGEOID-2000 and GPS ellipsoidal heights has been used to replace the traditional spirit leveling and mapping, called GPS mapping.     

High Precision Geoid Models HKGEOID-2000 for Hong Kong and SZGEOID-2000 for Shenzhen,China

High precision geoid models HKGEOID-2000 for Hong Kong and SZGEOID-2000 for Shenzhen, China, have been developed with a hybrid approach of so-called sequential processing, using high precision GPS/leveling data, land and sea gravity anomalies, and digital terrain models. These two local geoid models have the same 1-km resolution. The estimated accuracy (external accuracy) is better than 1.7 cm for HKGEOID-2000 and 1.4 cm for SZGEOID-2000. Some common areas are covered by HKGEOID-2000 and SZGEOID-2000. So these two geoid models, along with high quality GPS/leveling data collected on the overlapping areas, can be used to detect the systematic bias between HKGEOID-2000 and SZGEOID-2000, as well as the difference between Hong Kong Principal Datum and 1956 yellow sea height datum of China, yielding RMS errors of 1.011 m and 1,003 m, respectively. Moreover, HKGEOID-2000, along with GPS ellipsoidal heights, is employed to determine the errors of the “orthometric heights” from purely trigonometric heighting, yielding an RMS error of 0.102 m. The combination of SZGEOID-2000 and GPS ellipsoidal heights has been used to replace the traditional spirit leveling and mapping, called GPS mapping.     

基于Stokes-Helmert边值理论的大地水准面计算模型改化及分析检验

针对源自Stokes-Helmert边值理论的大地水准面计算模型适用性问题,首先详细解析了Stokes-Helmert方法的基本原理和计算模型,研究探讨了该方法与其他现有方法的关联性和差异性,分析比较了Stokes-Helmert计算模型和相同类型模型的技术特点,提出了Stokes-Helmert计算模型的实用改化方法,同时设计了4种分步改化模型的实验验证方案,使用超高阶位模型EGM2008作为数值模拟标准场,对分步改化模型的计算效果进行了数值精度检核,并开展了数据观测噪声影响评估检验,得出了一些有参考和实际应用价值的研究结论。在一定条件下,使用改化Stokes-Helmert方法计算(似)大地水准面可获得优于2cm的内符合度。     

利用T/P 卫星高度计资料调和分析南海潮汐信息

利用j,v模型调和分析1992～2002年共10 a的TOPEX/Poseidon(T/P)海面高度距平资料,提取了南海K1,O1,P1,Q1,M2,S2,N2和K2等8个主要分潮的潮汐调和常数。分析比较了卫星上下行轨道的19个交叉点的振幅和迟角,其中M2,S2,K1和O1的平均向量均方根偏差分别是1.5,1.1,2.5和1.4 cm;将交叉点的调和常数与TPXO7.2模式的结果进行了比较,结果表明M2,S2,K1和O1分潮振幅的绝对平均误差均小于3 cm,迟角的最大绝对平均误差为7.8°。选取了与卫星轨道较近的8个验潮站,对验潮站的实测数据调和常数和本文所得调和常数进行了比较,结果显示K1分潮的向量均方根偏差为4.7 cm,M2分潮的向量均方根偏差为3.7 cm。论文结果表明利用j,v模型调和分析方法对南海海域卫星高度计资料进行潮汐信息提取是可靠的,并可为局部重力场的研究提供海洋潮汐改正数据,有一定的参考价值。     

Comparison of geosat and gravimetric geoid profiles in the North Sea

Sea surface height profiles derived from 2‐year, repeat track, Geosat altimeter data have been compared with a regional gravimetric geoid in the western North Sea, computed using a geopotential model and terrestrial gravity data. The comparison encompasses 18 Geosat profiles covering a 750 × 850 km area of the North Sea. After a second‐order polynomial was used to model the long‐wavelength differences which cannot be clearly separated over an area of this size, results show agreement to better than ±3 cm for wavelengths between approximately 20 and 750 km. In regions where terrestrial gravity data were not available to improve the geoid, similar comparisons with the OSU91A geopotential model alone show differences of up to ±6 cm. This illustrates the importance of incorporating local gravity data in regional geoid computations, and partly validates the regional gravimetric geoid solution and Geosat sea surface profiles in the western North Sea. It is concluded that, in marine areas where the sea surface topography is known to be small in magnitude, Geosat sea surface profiles can act as an independent control on gravimetric geoids in the medium‐wavelength range.     

Estimating the sea surface dynamic topography from Geosat altimetry data

Optimal interpolation method is applied to Geosat altimetry data both to remove orbit error and to separate temporal mean sea surface dynamic topography (SSDT) from temporal fluctuations around the mean. Loss of long-wavelength oceanic signals at orbit error reduction procedure is smaller in this method than that in conventional collinear methods, but the areal average height over the study domain is still removed as the orbit error. The fluctuation SSDT is quantitatively evaluated by sea level data from tide gauge stations at Japanese islands. The correlation coefficient of the two sea-level variations is 0.83 when the loss of the areal average is compensated by the seasonal variation of the areal average height determined from the climatological monthly-mean SSDT. In addition, the improvement of the geoid model by combined use of Seasat altimetry data and hydrographic data is validated through the estimated temporal mean SSDT. In a local area where hydrographic data contemporary with the Seasat mission exist, the geoid model has been significantly improved so that the absolute SSDT can be determined from combination of the altimetry data and geoid model; the absolute SSDT describes the onset event of a quasi-stationary large meander of the Kuroshio south of Japan very well. Outside this local area, however, errors of several tens of centimeters still remain in the improved geoid model.     

南桑威奇俯冲带大地水准面异常及其地幔对流特征分析

以南桑威奇俯冲带为例,根据EGM2008超高阶地球重力场模型、卫星重力数据为基础,利用移去-恢复原理计算了研究区大地水准面,实现了研究区不同场源深度大地水准面异常信息的分离,根据Runcorn模型计算了研究区小尺度地幔对流应力场,并结合天然地震空间展布和前人研究成果,对俯冲带结构特征与地幔对流模式进行了探讨。结果表明:南桑威奇俯冲带具有俯冲倾角较大、地震震级较低、弧前侵蚀明显等典型的马里亚纳型俯冲带特征,俯冲带南北部俯冲深度存在明显差异,中段偏北俯冲深度可达500 km;受到软流圈与上地幔上部物质密度差异的控制,东斯科舍海脊下存在沿海脊轴向南流动强地幔流;俯冲带结构与小尺度地幔对流应力场具有很强的相关性。本研究对于搞清南桑威奇俯冲带深部构造特征,理解俯冲运动、地幔对流方向及其动力控制机制提供了新的研究思路和方法。     

Combined satellite altimetry and shipborne gravimetry data processing

The recovery of quantities related to the gravity field (i.e., geoid heights and gravity anomalies) is carried out in a test area of the central Mediterranean Sea using 5' × 5' marine gravity data and satellite altimeter data from the Geodetic Mission (GM) of ERS‐J. The optimal combination of the two heterogeneous data sources is performed using (1) the space‐domain least‐squares collocation (LSC) method, and (2) the frequency‐domain input‐output system theory (IOST). The results derived by these methods agree at the level of 2 cm in terms of standard deviation in the case of the geoid height prediction. The gravity anomaly prediction results by the same methods vary between 2.18 and 2.54 mGal in terms of standard deviation. In all cases, the spectral techniques have a much higher computational efficiency than the collocation procedure. In order to investigate the importance of satellite altimetry for gravity field modeling, a pure gravimetric geoid solution, carried out in a previous study for our lest area by the fast collocation approach (FCOL), is used in comparison with the combined geoid models. The combined solutions give more accurate results, at the level of about 15 cm in terms of standard deviation, than the gravimetric geoid solution, when the geoid heights derived by each method are compared with TOPEX altimeter sea surface heights (SSHs). Moreover, nonisotropic power spectral density functions (PSDs) can be easily used by IOST, while LSC requires isotropic covariance functions. The results show that higher prediction accuracies are always obtained when using a priori nonisotropic information instead of isotropic information.     

Assessment of Radar Waveform Retracked Jason-2 Altimetry Sea Surface Heights Near Taiwan Coastal Ocean

This study focuses on assessing the accuracy of 20-Hz waveform retracked Jason-2 (J-2) altimetry sea surface heights (SSHs) in the vicinity of Taiwan by comparisons with the TOPEX/Poseidon (T/P) 10-Hz SSHs and sea level data from the Anping tide gauge. The study areas exhibit high, medium, and low amplitudes of ocean tides and contain diverse bathymetries with depths of 0–4000 m. The performance of Offset Center of Gravity (OCOG), threshold, modified threshold, and ice retrackers was examined by comparing the retracked SSHs with Earth Gravitational Model 2008 (EGM08) geoid via the use of the improvement percentages (IMPs). The results indicate that both altimetry measurements are significantly improved by waveform retracking techniques, with a maximum IMP of 46.6% for T/P and 82.0% for J-2, and the optimal achievement of retrackers is influenced by the characteristics of the study areas. In addition, valid retracked J-2 SSHs are much closer to shorelines than T/P. A comparison of retracked J-2 data with Anping tide gauge records reveals that applying the optimal retracking algorithms reduces the root mean squares of differences and increases the number of valid measurements.     

RTM Gravity Forward-Modeling Using Topography/Bathymetry Data to Improve High-Degree Global Geopotential Models in the Coastal Zone

We apply the residual terrain modeling (RTM) technique for gravity forward-modeling to successfully improve high-resolution global gravity fields at short spatial scales in coastal zones. The RTM scheme is combined with the concept of rock-equivalent topography, allowing to use a single uniform constant mass-density in the RTM forward-modeling, both at land and sea. SRTM30_PLUS bathymetry is merged with higher-resolution SRTM V4.1 land topography, and expanded into spherical harmonics to degree 2160, yielding a new and consistent high-degree RTM reference surface. The forward-modeling performance is demonstrated in coastal zones of Greece and Canada using ground-truth vertical deflections, gravity from land and shipborne gravimetry, and geoid heights from GPS/leveling, with improvements originating from bathymetry clearly identified. We demonstrate that the SRTM30_PLUS bathymetry carries information on gravity field structures at spatial scales less than 5 arc minutes, which can be used to augment EGM2008 in (rugged) coastal zones, both over land and marine areas. This may be of value (i) to partially reduce the signal omission error in EGM2008/GOCE-based height transfer in areas devoid of dense gravity data, (ii) to fill the gap between land gravity and shipborne gravity along rugged coastlines, and (iii) for the development of next-generation altimetric gravity fields.     

Adjustment and filtering of seasat altimetry with the least squares response technique

The least squares response technique has been applied to SEASAT altimetry data for the adjustment and filtering of individual profiles without the use of standard crossover analyses. Based on a multichannel input-single output model, the technique allows a straightforward multi-correlation of the altimetry information with different types of forcing inputs. In addition, it can account for other known constituents such as datum biases, linear trends, periodic constituents of known periods and arbitrarily defined user-specified functions. In this paper, results of this approach are compared with geoid profiles computed from standard geopotential models such as GEM-10B which show a consistently good agreement with rms differences found at the level of the errors at long wavelengths of the geoid model used.     

东印度洋海域最优深度基准面模型构建

利用东印度洋海域周边长期验潮站实测数据、TOPEX/Poseidon等系列卫星测高反演结果,评估了DTU10,EOT11a,FES2014,GOT4.8,OSU12和TPXO8六种全球潮汐模型精度,根据卫星测高结果给出了浅水分潮改正量和长周期分潮改正量的经验模型,又在此基础上分析并构建了研究区域精度最优的深度基准面模型。考虑到全球潮汐模型在近岸的影响因素及验潮站位置,将13个验潮站分成开阔海域与近海海域两类,与潮汐模型的对比,结果表明,DTU10和FES2014模型分别在开阔海域和近海海域精度最优。根据潮汐模型在不同分潮处的精度,如EOT11a模型在O1和K1分潮处精度较高,DTU10在N2,M2,S2和K2分潮处精度较高等,分别构建了开阔海域与近海海域的组合深度基准面模型,计算得知误差分别为11.33和20.95 cm,其精度显著提高。     

Improvement in the Determination of the Marine Geoid by Estimating the Bathymetry from Altimetry and Depth Soundings

Abstract

The contribution of bathymetry to the estimation of gravity field related quantities is investigated in an extended test area in the Mediterranean Sea. The region is located southwest of the island of Crete, Greece, bounded between 33^? ≤ ? ≤ 35^? and 15^? ≤ λ ≤ 25^?. Gravity anomalies from the KMS99 gravity field and shipborne depth soundings are used with a priori statistical characteristics of depths in a least-squares collocation procedure to estimate a new bathymetry model. Two different global bathymetry models, namely JGP95E and Sandwell and Smith V8, are used to derive the depth a priori statistical information, while the estimated model is compared against both the global ones and the shipborne depth soundings to assess whether there is an improvement. Various marine geoid models are estimated using ERS1 and GEOSAT Geodetic Mission altimetry and shipborne gravity data. In that process, the effect of the bathymetry is computed using both the estimated and the original depths through a residual terrain modeling reduction. The TOPEX/Poseidon Sea Surface Heights, known for their high accuracy and precision, and the GEOMED solution for the geoid in the Mediterranean are used as control for the validation of the new geoid models and to assess the improvement that the estimated depths offer to geoid modeling. The results show that the newly estimated bathymetry agrees better (by about 30 to 300 m) with the shipborne depth soundings and provides smoother residual geoid heights and gravity anomalies (by about 8–20%) than those from global models. Finally, the achieved accuracy in geoid modeling ranges between 6 and 10 cm (1σ).     

Using advanced soft computing techniques for regional shoreline geoid model estimation and evaluation

This study aims at evaluating the global geoid model for a regional shoreline fitting using advanced soft computing techniques and global navigation satellite system/leveling measurements. Artificial neural networks, fuzzy logic, and least square support vector machine models are developed and used to fit the global geoid model for the north coastal Egyptian line. In addition, a novel estimation geoid model is designed and evaluated based on the latest global geoid models. The results of the three estimation models show that they can be used to correct the shoreline geoid model, in terms of root mean square error that ranges from 1.7 to 8.5?cm. Moreover, it is found that the least square vector machine model is a competitive approach with certain advantage in solving complex problems represented by missing data.     

Leveling the Sea Surface Using a GPS-Catamaran Special Issue: Jason-1 Calibration/Validation

In the framework of the TOPEX/Poseidon and Jason-1 CNES-NASA missions, two probative experiments have been conducted at the Corsica absolute calibration site in order to determine the local marine geoid slope under the ascending TOPEX/Poseidon and Jason-1 ground track (No. 85). An improved determination of the geoid slope was needed to better extrapolate the offshore (open-ocean) altimetric data to on-shore tide-gauge locations. This in turn improves the overall precision of the calibration process. The first experiment, in 1998, used GPS buoys. Because the time required to cover the extended area with GPS buoys was thought to be prohibitive, we decided to build a catamaran with two GPS systems onboard. Tracked by a boat at a constant speed, this innovative system permitted us to cover an area of about 20 km long and 5.4 km wide centered on the satellites' ground track. Results from an experiment in 1999 show very good consistency between GPS receivers: filtered sea-surface height differences have a mean bias of ?0.2 cm and a standard deviation of 1.2 cm. No systematic error or distortions have been observed and crossover differences have a mean value of 0.2 cm with a standard deviation of 2.7 cm. Comparisons with tide gauges data show a bias of 1.9 cm with a standard deviation of less than 0.5 cm. However, this bias, attributable in large part to the effect of the catamaran speed on the waterline, does not affect the geoid slope determination which is used in the altimeter calibration process. The GPS-deduced geoid slope was then incorporated in the altimeter calibration process, yielding a significant improvement (from 4.9 to 3.3 cm RMS) in the agreement of altimeter bias determinations from repeated overflight measurements.