Spectral analysis of quasi-stationary sea surface topography from GRACE mission

The mean sea surface height (MSSH) refers to the average of the long-term sea height. The quasi-sea surface topography (QSST) is usually defined as the height difference between the MSSH and the geoid. As to 100 years of time yardstick of geodesy, the time that satellite altimetry data sets spanned is relatively shorter, in this paper, the QSST refers to the residual sea surface height (RSSH) that shows the height dif-ference between MSSH derived from altimetry and the geoid[1]. As w…     

Assimilation of geodetic dynamic ocean topography using ensemble based Kalman filter

Estimation of ocean circulation is investigated via assimilation of satellite measurements of the dynamic ocean topography (DOT) into the global finite-element ocean model (FEOM). The DOT was obtained by means of a geodetic approach from carefully cross-calibrated multi-mission altimeter data and GRACE gravity fields. The spectral consistency was achieved by consistently filtering both, the sea surface and the geoid. The filter length is determined by the spatial resolution of the gravity field and corresponds to approximately 241 km half width for the GRACE-based gravity field model ITG-Grace03s.The assimilation of the geodetic DOT was performed by employing a local singular evolutive interpolated Kalman (SEIK) filter in combination with the method of weighting of observations. It is shown that this approach leads to a successful assimilation technique that reduced the RMS difference between the model and the data from 16 cm to 5 cm during one year of assimilation. The ocean model returns an optimized mean dynamic ocean topography. The effects of assimilation on transport estimates across several hydrographic World Ocean Circulation Experiment (WOCE) sections show improvements compared to the FEOM run without data assimilation. As a result of the assimilation, DOT estimates are available in the polar or coastal regions where the geodetic estimates from satellite data alone are not adequate. Furthermore, more realistic features of the ocean can be seen in these areas compared to those obtained using the filtered data fields.     

用GRACE卫星跟踪数据反演地球重力场

利用141天GRACE卫星观测资料,包括K波段、星载加速度和卫星轨道数据,反演了80阶地球重力场模型IGGGRACE01S,该模型在半波长为500km的空间分辨率上,确定大地水准面的精度约为0012m,中长波(<80阶)精度优于重力卫星发射以前研制的重力场模型. 与EIGEN_GRACE02S、EIGEN_CHAMP03S和EGM96模型的位系数相比,该模型系数最接近于EIGEN_GRACE02S,与另两个模型差异较大. 比较几种模型确定的全球重力异常和大地水准面起伏,结果发现IGGGRACE01S与EIGEN_GRACE02S模型的计算结果比较接近,与EGM96模型结果差异较大,差别较大地区主要在南极等地区. 对于中国大陆,比较IGGGRACE01S模型（前72阶）计算的重力异常和NIMA重力异常数据（25°×25°网格）,两者之间的标准偏差为48mGal.     

由GOCE引力场模型和CNES-CLS2010平均海面高计算的稳态海面地形

利用欧空局发布的三组GOCE引力场模型及CNES-CLS 2010平均海面高数据,计算得到了全球的稳态海面地形,进而得到了全球地转流速度图.在此基础上重点对黑潮进行了对比分析.结果表明:GOCE不同组解的稳定性较好,所计算的稳态海面地形的差异基本在厘米量级内,这间接表明了GOCE引力场模型提供的大地水准面的精度达到了厘米量级.此外,通过将GOCE与GRACE相应结果进行对比发现,GOCE可提供更多的局部信息,特别是对于流速快、水流窄的边界流,如黑潮、墨西哥湾流等,GOCE所得结果更加清晰,速度也更精确.     

基于区域滤波的GOCE稳态海面动力地形和地转流

基于频域法,利用最新的GOCE卫星重力场模型和卫星测高数据计算了稳态海面动力地形.结合海洋表层漂流浮标的观测结果,对稳态海面动力地形进行了最优空间滤波尺度分析,给出了区域、纬度带和全球稳态海面动力地形的最优空间滤波尺度因子.在此基础上,给出了全球和区域地转流.结果表明:在中高纬度和全球区域,可以分别获得空间尺度优于102km和127km的稳态海面动力地形信息.与海洋表层漂流浮标对比可知,在强流区域,采用稳态海面动力地形得到的地转流速可以解释观测浮标流速的70%;在中高纬度区域,由GOCE重力场得到的地转流略优于对应的GRACE结果;在近赤道区域,由GOCE重力场得到的地转流精度略低于对应的GRACE结果;在北大西洋和阿古拉斯强流区域,由GOCE得到的地转流场明显优于对应的GRACE结果,其精度分别提高了16%和24%.     

基于卫星测高交叉点的海洋表面地转流速度

在流体静力平衡状态下,海洋Coriolis力和压力梯度平衡就形成地转流,世界上大多数海流都近似为地转流.本文利用卫星测高交叉点方法计算海洋表面地转流速度,分析了利用测高交叉点计算地转流速度的不确定性,上升和下降弧段的海面倾斜在分辨率50 km上可以达到10^-7量级,才可能获得优于10 cm/s的地转流速度.在低纬度或者纬度接近卫星轨道倾角的地区,由交叉点方法计算的地转流速度精度低于中纬度地区.以中国台湾东部黑潮为试验区,利用最新的中国台湾周边海域大地水准面模型参考场计算高精度的大地水准面高,利用TOPEX/Poseidon和Jason-1的GDR数据(2002~2005年)计算海面高,然后计算交叉点的动力高,确定交叉点的地转流速度,结果与中国台湾NCOR(National Center for Ocean Research)的流速基本一致.     

Comparison of global geopotential models from the champ and grace missions for regional geoid modelling in Turkey

The continuous efforts on establishment and modernization of the geodetic control in Turkey include a number of regional geoid models that have been determined since 1976. The recently released gravimetric Geoid of Turkey, TG03, is used in geodetic applications where GPS-heights need to be converted to the local vertical datum. To reach a regional geoid model with improved accuracy, the selection of the appropriate global geopotential model is of primary importance. This study assesses the performance of a number of recent satellite-only and combined global geopotential models (GGMs) derived from CHAMP and GRACE missions’ data in comparison to the older EGM96 model, which is the underlying reference model for TG03. In this respect, gravity anomalies and geoid heights from the global geopotential models were compared with terrestrial gravity data and low-pass filtered GPS/levelling data, respectively. Also, five new gravimetric geoid models, computed by the Fast Fourier Transform technique using terrestrial gravity data and the geopotential models, were validated at the GPS/levelling benchmarks. The findings were also compared with the validation results of the TG03 model. The tests showed that as it was expected any of the high-degree combined models (EIGEN-CG03C, EIGEN-GL04C, EGM96) can be employed for determining the gravity anomalies over Turkey. In the west of Turkey, EGM96 and EIGEN-CHAMP03S fit the GPS/levelling surface better. However, all the tested GGMs revealed equal performance when they were employed in gravimetric geoid modelling after de-trending the gravimetric geoid model with corrector surface fitting. The new geoid models have improved accuracy (after fit) compared to TG03.     

地球重力场恢复中的位旋转效应

分析了地球自转引起的位旋转效应公式中采用近似速度的影响. 对一组GFZ的快速科学轨道、一组TUM的约化动力法轨道以及一组GFZ的事后科学轨道,计算了星历提供的速度与只有地球引力场对卫星产生作用时的卫星速度的差值,其中参考重力场模型分别采用EGM96、EIGEN2和EIGEN_CG01C. 通过比较得出：轨道数据与EIGEN2地球重力场模型的自恰性优于EGM96和EIGEN_CG01C地球重力场模型. 速度差各分量的变化具有很明显的周期性且与卫星轨道的运行周期相吻合. 当要求在卫星轨迹处获得1m2/s2精度的扰动位时,也即要求位旋转效应公式中卫星速度的近似精度小于2mm/s时,GFZ的快速科学轨道、TUM的约化动力法轨道只需要剔除那些速度精度不满足要求的卫星轨迹点;当要求在卫星轨迹处获得05m2/s2精度的扰动位时,应当重新估算上述轨道的速度信息,或采用精度更高的GFZ事后科学轨道.     

Direct BEM for high-resolution global gravity field modelling

The paper presents a high-resolution global gravity field modelling by the boundary element method (BEM). A direct BEM formulation for the Laplace equation is applied to get a numerical solution of the linearized fixed gravimetric boundary-value problem. The numerical scheme uses the collocation method with linear basis functions. It involves a discretization of the complicated Earth’s surface, which is considered as a fixed boundary. Here 3D positions of collocation points are simulated from the DNSC08 mean sea surface at oceans and from the SRTM30PLUS_V5.0 global topography model added to EGM96 on lands. High-performance computations together with an elimination of the far zones’ interactions allow a very refined integration over the all Earth’s surface with a resolution up to 0.1 deg. Inaccuracy of the approximate coarse solutions used for the elimination of the far zones’ interactions leads to a long-wavelength error surface included in the obtained numerical solution. This paper introduces an iterative procedure how to reduce such long-wavelength error surface. Surface gravity disturbances as oblique derivative boundary conditions are generated from the EGM2008 geopotential model. Numerical experiments demonstrate how the iterative procedure tends to the final numerical solutions that are converging to EGM2008. Finally the input surface gravity disturbances at oceans are replaced by real data obtained from the DNSC08 altimetryderived gravity data. The ITG-GRACE03S satellite geopotential model up to degree 180 is used to eliminate far zones’ interactions. The final high-resolution global gravity field model with the resolution 0.1 deg is compared with EGM2008.     

Tide model comparison over the Southwestern Atlantic Shelf

Sea surface height (SSH) as measured by satellites has become a powerful tool for oceanographic and climate related studies. Whereas in the open ocean good accuracy has been achieved, more energetic dynamics and a number of calibration problems have limited applications over continental shelves and near the coast. Tidal ranges in the Southwestern Atlantic (SWA) continental shelf are among the highest in the world ocean, reaching up to 12 m at specific locations. This fact highlights the relevance of the accuracy of the tidal correction that must be applied to the satellite data to be useful in the region. In this work, amplitudes and phases of tidal constituents are extracted from five global tide models and three regional models and compared to the corresponding harmonics estimated from coastal tide gauges (TGs) and satellite altimetry data. The Root Sum Square (RSS) of the misfit of the common set of the five tidal constituents solved by the models (M₂, N₂, S₂, K₁ and O₁) is higher than 18 cm close to the coast for two of the regional models and higher than 24.5 cm for the rest of the models considered. Both values are too high to provide an accurate estimation of geostrophic non-tidal currents from satellite altimetry in the coastal region. On the other hand, the global model with the highest spatial resolution has a RSS lower than 4.5 cm over the continental shelf even when the non-linear M₄ overtide is considered. Comparison with in-situ current measurements suggests that this model can be used to de-tide altimetry data to compute large-scale patterns of SSH and associated geostrophic velocities. It is suggested that a local tide model with very high resolution that assimilates in-situ and satellite data should meet the precision needed to estimate geostrophic velocities at a higher resolution both close to the coast and over the Patagonian shelf.     

A comparison of GOCE gravitational models with EGM2008

Four new gravity field models from GOCE, two of them combined with GRACE, are compared here with EGM2008. The objectives are to look into the differences in consecutive ranges of the spherical harmonic expansion globally as well as in selected geographical regions and in the regions of the various data sources used for EGM2008. In general, GOCE is able to contribute to improved global gravity models in the spherical harmonic range between 120 and 200 (and above). The agreement between EGM2008 and the GOCE models is very good in well-surveyed regions such as North America, Europe and Australia, with geoid RMS-differences on the order of 4–6 cm. In other regions, where the surface gravity data available for the development of EGM2008 were poor, such as South America, Africa, South-East Asia or China the RMS-differences are on a level of 30 cm. Here GOCE leads to a significant improvement. These findings are confirmed by the analysis of the areas of the various EGM2008 data sources. In the regions of the so-called “fill-in” data of EGM2008 RMS-geoid height differences are high. In Antarctica GOCE also gives important improvements in terms of spatial resolution and accuracy. In general, the agreement between EGM2008 and the GOCE-models up to degree and order (d/o) 200 is good, with a global (excluding the polar gaps of GOCE orbits, throughout) geoid difference RMS of 11 cm, in the ocean areas 8 cm and 20 cm in the continental areas. GOCE models are better suited for ocean circulation studies because no prior ocean information enters into the data reduction process, as it is the case when deducing gravity anomalies from an altimetric mean sea surface. On the other hand, the good consistency between GOCE-models and EGM2008 in ocean areas very likely indicates that the influence of ocean circulation information on EGM2008 is rather small. The four tested GOCE models behave similarly except at the highest latitudes where GOCE lacks data due to its orbit inclination of 96.5° and some form of regularization which has to be applied.     

One technique for refining the global Earth gravity models

The results of the theoretical and experimental research on the technique for refining the global Earth geopotential models such as EGM2008 in the continental regions are presented. The discussed technique is based on the high-resolution satellite data for the Earth’s surface topography which enables the allowance for the fine structure of the Earth’s gravitational field without the additional gravimetry data. The experimental studies are conducted by the example of the new GGMplus global gravity model of the Earth with a resolution about 0.5 km, which is obtained by expanding the EGM2008 model to degree 2190 with the corrections for the topograohy calculated from the SRTM data. The GGMplus and EGM2008 models are compared with the regional geoid models in 21 regions of North America, Australia, Africa, and Europe. The obtained estimates largely support the possibility of refining the global geopotential models such as EGM2008 by the procedure implemented in GGMplus, particularly in the regions with relatively high elevation difference.     

Gravity Anomaly Map over The Black Sea Using Corrected Sea Surface Heights from ERS1, ERS2 and Topex/Poseidon Satellite Altimetry Missions

Approximately nine-year data from ERS1, ERS2 and TOPEX/POSEIDON (T/P) satellite altimetry missions have been used for the recovery of gravity anomalies over the Black Sea. The Corrected Sea Surface Height product of Aviso/Altimetry has been proven to be homogeneous after a cross-over adjustment. The Least Squares Collocation method was applied in a so-called remove-restore procedure. The residual geoid heights, obtained by subtracting EGM96 geoid heights from cross-over adjusted sea surface heights, were inverted to recover residual gravity anomalies in a grid structure over the Black Sea. Finally, EGM96 free air gravity anomalies were added to the predicted residual gravity anomalies to obtain the free air gravity anomalies. In order to check the consistency with respect to an external source, these computed free air gravity anomalies were compared to ship gravity observations, and to alternative satellite altimetry derived gravity anomalies. Comparisons with the observed gravity data yielded that external consistency of the gravity anomalies computed in this study is about 3 mGal for individual ship tracks. Overall external consistency in the test area is 4.8 mGal. Comparison with other satellite altimetry derived gravity anomalies presented a good agreement.     

Regional gravimetric quasi-geoid model and transformation surface to national height system for Turkey (THG-09)

Turkish regional geoid models have been developed by employing a reference earth gravitational model, surface gravity observations and digital terrain models. The gravimetric geoid models provide a ready transformation from ellipsoidal heights to the orthometric heights through the use of GPS/leveling geoid heights determined through the national geodetic networks. The recent gravimetric models for Turkish territory were computed depending on OSU91 (TG-91) and EGM96 (TG-03) earth gravitational models. The release of the Earth Gravitational Model 2008 (EGM08), the collection of new surface gravity observations, the advanced satellite altimetry-derived gravity over the sea, and the availability of the high resolution digital terrain model have encouraged us to compute a new geoid model for Turkey. We used the Remove-Restore procedure based on EGM08 and applied Residual Terrain Model (RTM) reduction of the surface gravity data. Fast Fourier Transformation (FFT) was then used to obtain the residual quasigeoid from the reduced gravity. We restored the individual contributions of EGM08 and RTM to the whole quasi-geoid height (TQG-09). Since the Helmert orthometric height system is adopted in Turkey, the quasi-geoid model (TQG-09) was then converted to the geoid model (TG-09) by making use of Bouguer gravity anomalies and digital terrain model. After all we combined a gravimetric geoid model with GPS/leveling geoid heights in order to obtain a hybrid geoid model (THG-09) (or a transformation surface) to be used in GPS applications. The RMS of the post-fit residuals after the combination was found to be ± 0.95 cm, which represents the internal precision of the final combination. And finally, we tested the hybrid geoid model with GPS/leveling data, which were not used in the combination, to assess the external accuracy. Results show that the external accuracy of the THG-09 model is ± 8.4 cm, a precision previously not achieved in Turkey until this study.     

基于重力卫星几何轨道线性化的地球重力场反演方法

传统动力学法的观测方程以6个初始轨道参数和先验力模型为初值进行线性化,其线性化误差随积分弧长拉长而增大.本文直接以重力卫星的几何观测轨道为初值进行线性化,其线性化误差与轨道弧长无关,且不需要初始重力场模型和初始轨道参数.导出了基于卫星轨道观测值反演重力场模型的相关公式,利用JPL公布的RL02版本2008年全年的GRACE双星轨道数据和加速度计数据解算了90阶次的地球重力场模型TJGRACE01S,并以EGM2008模型为基准与其他模型进行了比较分析,结果表明:TJGRACE01S模型直到90阶次的大地水准面累积误差为17.6 cm,优于同阶次的EIGEN-CHAMP03S和EIGEN-CHAMP05S模型,前27阶位系数整体精度优于EIGEN-GRACE01S,前15阶位系数整体精度与EIGEN-GRACE02S模型精度大致相当.利用美国8221个GPS水准点数据的分析结果也表明,本文模型也优于同阶次的EIGEN-CHAMP03S和EIGEN-CHAMP05S模型.     

海潮对卫星重力场恢复的影响

本文讨论了海潮对卫星重力测量的影响问题. 首先介绍了海潮对卫星重力测量影响的基本理论;采用FES02和TPXO6海潮模型计算了海潮负荷对卫星重力结果前60阶的影响;并用两个模型之间的差异作为海潮模型精度的估计量,据此计算了海潮模型误差对卫星重力结果的影响. 与GRACE恢复的重力场精度的比较说明：海潮对重力场40阶以下的影响都超过了目前重力场恢复精度;尽管由于卫星测高技术的发展,海潮模型的精度有了很大的提高,但目前的全球海潮模型用于GRACE重力场恢复的前12阶的改正还是不够精确. 另外,我们也利用中国东海和南海潮汐资料以及FES02海潮模型讨论了中国近海潮汐效应对GRACE观测的影响. 结果说明该影响与海潮模型的误差相当. 这反映了当前海潮模型的不确定度,因此通过结合全球验潮站资料有望提高海潮对卫星重力测量的改正精度.     

Combined Regional Gravity Model of the Andean Convergent Subduction Zone and Its Application to Crustal Density Modelling in Active Plate Margins

The Central Andean subduction system is one of the most active geological structures on Earth. Although there have been a few previous studies, the structure and dynamics of the system are still not well understood. In the present study, we determine a combined regional gravity model of the Andean convergent subduction region for constraining lithospheric models. After a thorough validation and cleaning of the terrestrial gravity and height databases, the method of Least Squares Collocation was applied to consistently combine terrestrial and satellite gravity data, putting much emphasis on the stochastic modelling of the individual data components. As a result, we computed the first high-resolution regional gravity model of the study region that includes GOCE satellite gravity information. The inclusion of GOCE is an essential distinction from the independent global gravity model EGM2008. Validation against EGM2008 reveals that our regional solution is very consistent in regions where terrestrial gravity data are available, but shows systematic differences in areas with terrestrial data gaps. Artefacts in the EGM2008 of up to 150 mGal could be identified. The new combined regional model benefits from the very homogeneous error characteristics and accuracy of GOCE gravity data in the long-to-medium wavelengths down to 80–100 km. Reliable density modelling became possible also in the region of Central Andes, which lacks terrestrial gravity data. Finally, density models were adapted to fit the new regional gravity field solution. The results clearly demonstrate the capabilities of GOCE to better constrain lithospheric models.     

近四年全球海水质量变化及其时空特征分析

本文利用卫星重力、卫星测高和海洋温盐数据反演计算全球海水质量变化,并分析其时空变化特征.卫星重力数据利用2003年1月~2006年12月的GRACE月时变重力场球谐系数,同时考虑替换一阶项和C₂₀项,并进行了相关误差滤波、高斯滤波和陆地水文信号泄漏改正,计算得到海洋等效水高变化;利用相同时间跨度的卫星测高数据和海洋温度、盐度水文观测数据,计算全球海平面变化和比容海平面变化,反演得到海水质量变化.反演的两种海水质量变化的年际变化特征一致性较好.三种数据得到的长期趋势变化,与1993~2003年的结果相比,可以看出,海水质量变化加速,并已成为全球海平面上升的主要因素.     

应用多种来源重力异常编制中国海陆及邻区空间重力异常图及重力场解读

本文通过分析陆地实测空间重力异常数据、海洋船载测量空间重力异常数据、卫星测高重力异常,布格重力异常数据、EGM2008地球重力模型数据等多种来源数据的性质和精度,并对相关数据进行对比,研究了编制1:500万中国海陆空间重力异常图的数据使用方案和技术方法.在地形较为平坦、实测数据分布均匀的陆区,使用实测数据,在地形复杂,实测数据稀少以及没有实测数据的陆区或岛屿,利用布格重力异常反推空间异常的方法合成平均空间重力数据,西藏地区的数据对比实验证明合成平均空间重力异常数据是一种有效的数据补充.利用三观测列方差分解法在南海地区对船载测量空间重力数据和美国SS系列及丹麦DNSC08GRA卫星重力数据进行了方差分解计算,结果表明不同来源的卫星测高重力数据具有很大的一致性,数据精度较以往有了很大的提高.海区空间重力数据使用原则是在船载重力测量数据校准下,全面使用卫星测高重力数据进行编图.海陆过渡区的异常处理应以EGM2008地球重力模型重力场为基准参考场,实现海陆异常平缓过渡,无缝连接.对中国海陆空间重力异常场进行了小波变换处理,对空间重力异常场进行了解读,勾画出三横四竖的一级重力梯级带及其所围限的8个一级重力异常区,并划分了二级重力异常区和梯级带,为块体构造学体系中大地构造格架的建立提供了地球物理证据.     

Seasonal water storage change of the Yangtze River basin detected by GRACE

1 Introduction Large-scale mass redistribution, or temporal varia- tion of mass within the Earth system, the driving force of interactions between solid Earth and geophysical fluids envelope (i.e., atmosphere, ocean, and hydro- sphere), is an important geophysical process critical to human life. Most of the interactions between solid Earth and the atmosphere/oceans happen at seasonal and inter-annual time scales. One important contribu- tor of mass redistribution at seasonal and inter-annual …