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

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

GOCE卫星数据确定的全球稳态海面地形和海表地转流

利用GOCE卫星235天的观测数据恢复了200阶次的重力场模型SWJTU-GO01S,结合欧空局提供的最新GOCE重力场模型和CNES-CLS 2011平均海面高模型,计算了全球稳态海面地形和海表地转流,并采用GRACE模型、多源数据同化模型和海洋浮标观测数据对GOCE模型的计算结果进行对比分析.结果表明:由于重力场模型精度和分辨率较高,GOCE计算结果所需的滤波半径小于GRACE结果;GOCE和GRACE模型的计算结果与CNES-CLS09稳态海面地形差异的RMS分别为6 cm和7 cm左右;与海洋浮标实测数据对比发现,GOCE和GRACE的计算结果与实测数据差异明显,但GOCE的计算结果优于GRACE结果,而SWJTU-GO01S与DIR-R4和TIM-R4模型在全球范围内具有较好的一致性.整体而言,GOCE比GRACE数据的计算结果可以反映更小尺度地转流,且计算的精度更高;海洋环流结果和水准数据的对比表明SWJTU-GO01S与DIR-R4和TIM-R4模型的精度符合性较好,三者计算的地转流精度基本相当.     

GRACE卫星资料确定的稳态海面地形及其谱特征

介绍由新一代卫星重力得到的重力场模型GGM01C, 并据此计算构造新的全球高精度高分辨率的稳态海面地形. 利用信号的频谱分析方法对以不同阶次GGM01C模型求得的剩余海面高进行谱分析, 并与EGM96重力场模型结果作相应比较, 结果表明: 海面地形的谱结构非常稳定, 主要由长波构成, 中波部分占据一定比例, 但经纬向表现了不同的频谱特征；两个不同来源的重力场模型存在较大差异, 解释了Tapley据此两模型结果计算的洋流状态不同的主要原因.     

利用卫星测量技术和小波滤波方法探测表层地转流

根据海面地形的空域频谱特征, 提出用小波滤波器对海面地形进行滤波降噪的方法. 通过对全球及黑潮流系区域的海面地形进行小波滤波和高斯低通滤波降噪后所确定的地转流的比较, 显示前者较之后者更能表现地转流的局部特征. 依据CG01C卫星重力场模型和EGM96重力场模型, 分别联合由卫星测高确定的KMSS04平均海面高模型构造海面地形. 将据小波滤波降噪后所推算出的全球平均地转流、黑潮和热带太平洋区域地转流与海洋学结果的比较, 以及据此两不同重力场模型推算的大地水准面误差对计算地转流流速精度的对比都表明: 重力卫星确定的地球重力场模型较之以前存在的重力场模型在长波部分精度有较大提高. 联合现有卫星重力和卫星测高数据探测的全球平均地转流, 在大、中尺度上与海洋学结果相一致, 这表明从大地测量(空间)角度来研究洋流已达到较高的精度.     

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

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

利用卫星大地测量技术研究海面动力地形及地转流的进展

随着卫星大地测量技术的发展,平均海平面高度与大地水准面资料的分辨率和精度得到大幅提高,这为确定高精度高分辨率的海面动力地形提供了条件,进而也促进了海面动力地形及地转流研究的发展.本文从数据资料、计算方法、空间尺度分析、误差分析及多源数据融合五个方面总结了目前利用卫星大地测量技术研究海面动力地形及地转流的进展.在此基础上,对目前该研究中存在的问题进行了分析,并对该研究今后的重点发展方向进行了展望.     

南海海面高度和输运流函数: 全球 变网格模式结果

为了研究南海环流结构和变异及其与外部水域的关系, 我们建立了一个嵌套于全球大洋环流模式中的高分辨率中国近海环流数值模式. 给出模拟所得南海月平均以及年平均的海面高度和流函数分布, 与TOPEX/POSEIDON资料比较表明, 所得海面高度距平与观测十分一致. 基于这些结果, 讨论了南海的环流结构, 尤其是上层环流结构. 结果表明: 对于表层海水来说, 黑潮在冬、春和秋季均通过吕宋海峡入侵南海, 夏季则表层没有入侵. 但对于整个海水而言, 全年均有海水从太平洋通过吕宋海峡进入南海. 这一差异表明, 在夏季, 太平洋的海水是在次表层和中层入侵南海的.南海北部陆坡附近全年受气旋式环流控制. 夏季的南海南部反气旋流圈、越南东南离岸流和冬季的南海南部气旋流圈都得到了很好的再现. 南海海面高度和海面高度距平之间的差异明显. 表明, 在利用卫星高度计资料研究南海的上层环流时, 长期平均海面高度的空间分布有重要意义.     

GOCE卫星监测的中国区域重力异常

介绍了利用最新一代重力卫星GOCE(Gravity field and steady-state Ocean Circulation Explorer)重力场模型数据计算重力异常的原理和方法,采用最新发布的GOCE重力场模型数据(2009年10月-2010年7月)计算了中国区域(70°～130°E,15°～55°N)的重力异常,为了进一步分析重力异常与区域地质构造及地震活动性的对应关系,将计算结果与中国区域地形、地震活动区域等资料进行了对比分析,结果表明利用GOCE重力场模型数据计算的重力异常能够较好的反映区域地质构造分布特征,而且强震震中通常位于重力异常变化剧烈的高梯度带上.     

南海海面高度、动力地形和环流的周年变化——TOPEX/Poseidon卫星测高应用研究

用1993~1999年的TOPEX/Poseidon卫星测高资料, 分析了南海海面高度距平场(SSHA)的平均周年变化; 结合历史水文资料, 反演了多年平均的逐月海面动力地形; 探讨了南海动力地形及其所反映的上层环流季节特征和演变规律. 分析表明, 南海大尺度环流的周年演替可分为4个阶段. 冬季(11~2月)南海环流表现为以北部气旋环流为主的气旋型双圈结构, 相关的特征还包括吕宋海峡的黑潮入侵和加里曼丹岛西北外海的东北向离岸流. 春季(3~4 月)气旋型双圈结构解体, 北部的气旋型环流依然维持, 南部环流则向反气旋型演变, 大尺度环流结构呈现偶极子特征. 夏季(5~7 月)和秋季(8~10月)海盆内部不存在明显封闭的大尺度环流, 环流以西南-东北流向的季风急流为主要特征, 但夏、秋流态有较大差别. 5~7 月季风急流贴中南半岛北上, 在海南岛东南18°N附近沿地形折向东形成反气旋型弯曲, 穿越南海后再折向东北. 8~10 月季风急流在13°N附近即离开中南半岛海岸进入海盆中部, 其流态转变为气旋型. 在春、夏、秋三季, 黑潮的入侵都不明显. 上述变化规律显示南海环流的动力调整在季风盛期过后就已经开始.     

利用高精度复测GPS网进行中国大陆区域地壳运动特征分析

本文利用高精度复测GPS网,用多面函数拟合的方法计算出中国大陆地壳运动速度场模型,得到了有实用价值的中国地壳平面运动整体速度场图像和垂直速率等值线图像,同时计算出ITRF97参考框架下中国大陆地壳运动的区域背景场,该局域背景场与NNR-NUVELlA、Sillard、Larson 等全球板块运动模型所提供的欧亚板块的欧拉极方向和旋转速率在量级和大趋势上是一致的.据此对全球框架下的中国大陆区域地壳运动特征进行分析.     

Performance of GOCE and GRACE-derived mean dynamic topographies in resolving Antarctic Circumpolar Current fronts

Presently, two satellite missions, Gravity Recovery and Climate Experiment (GRACE) and Gravity field and steady-state Ocean Circulation Explorer (GOCE), are making detailed measurements of the Earth’s gravity field, from which the geoid can be obtained. The mean dynamic topography (MDT) is the difference between the time-averaged sea surface height and the geoid. The GOCE mission is aimed at determining the geoid with superior accuracy and spatial resolution, so that a more accurate MDT can be estimated. In this study, we determine the mean positions of the Antarctic Circumpolar Current fronts using the purely geodetic estimates of the MDT constructed from an altimetric mean sea surface and GOCE and GRACE geoids. Overall, the frontal positions obtained from the GOCE and GRACE MDTs are close to each other. This means that these independent estimates are robust and can potentially be used to validate frontal positions obtained from sparse and irregular in situ measurements. The geodetic frontal positions are compared to earlier estimates as well as to those derived from MDTs based on satellite and in situ measurements and those obtained from an ocean data synthesis product. The position of the Sub-Antarctic Front identified in the GOCE MDT is found to be in better agreement with the previous estimates than that identified in the GRACE MDT. The geostrophic velocities derived from the GOCE MDT are also closer to observations than those derived from the GRACE MDT. Our results thus show that the GOCE mission represents an improvement upon GRACE in terms of the time-averaged geoid.     

Impact of GOCE Level 1b data reprocessing on GOCE-only and combined gravity field models

The reprocessing of Gravity field and steady-state Ocean Circulation Explorer (GOCE) Level 1b gradiometer and star tracker data applying upgraded processing methods leads to improved gravity gradient and attitude products. The impact of these enhanced products on GOCE-only and combined GOCE+GRACE (Gravity Recovery and Climate Experiment) gravity field models is analyzed in detail, based on a two-months data period of Nov. and Dec. 2009, and applying a rigorous gravity field solution of full normal equations. Gravity field models that are based only on GOCE gradiometer data benefit most, especially in the low to medium degree range of the harmonic spectrum, but also for specific groups of harmonic coefficients around order 16 and its integer multiples, related to the satellite’s revolution frequency. However, due to the fact that also (near-)sectorial coefficients are significantly improved up to high degrees (which is caused mainly by an enhanced second derivative in Y direction of the gravitational potential — V_YY), also combined gravity field models, including either GOCE orbit information or GRACE data, show improvements of more than 10% compared to the use of original gravity gradient data. Finally, the resulting gradiometry-only, GOCE-only and GOCE+GRACE global gravity field models have been externally validated by independent GPS/levelling observations in selected regions. In conclusion, it can be expected that several applications will benefit from the better quality of data and resulting GOCE and combined gravity field models.     

Detecting surface geostrophic currents using wavelet filter from satellite geodesy

According to the features of spatial spectrum of the dynamic ocean topography (DOT),wavelet filter is proposed to reduce short-wavelength and noise signals in DOT. The surface geostrophic currents calculated from the DOT models filtered by wavelet filter in global and Kuroshio regions show more detailed information than those from the DOT models filtered by Gaussian filter. Based on a satellite gravity field model (CG01C) and a gravity field model (EGM96),combining an altimetry-derived mean sea surface height model (KMSS04),two mean DOT models are estimated. The short-wavelength and noise signals of these two DOT models are removed by using wavelet filter,and the DOT models asso-ciated global mean surface geostrophic current fields are calculated separately. Comparison of the surface geostrophic currents from CG01C and EGM96 model in global,Kuroshio and equatorial Pacific regions with that from oceanography,and comparison of influences of the two gravity models errors on the precision of the surface geostrophic currents velocity show that the accuracy of CG01C model has been greatly improved over pre-existing models at long wavelengths. At large and middle scale,the surface geostrophic current from satellite gravity and satellite altimetry agrees well with that from oceanography,which indicates that ocean currents detected by satellite measurement have reached relatively high precision.     

GRACE/GOCE扩展重力场模型确定我国1985高程基准重力位的精度分析

高精度高程基准重力位的确定往往依赖于高精度全球重力场模型,其对全球和区域高程基准的高精度统一非常关键,GRACE、GOCE卫星重力计划极大地提高了全球重力场模型中长波的精度.本文首先对GRACE/GOCE卫星重力场模型的内符合和外符合精度进行讨论分析,结果说明卫星重力模型的截断误差影响可达到分米级水平,在确定高程基准重力位时该影响不可忽略.利用EGM2008模型扩展GRACE/GOCE卫星重力场模型至2190阶,可有效减弱卫星重力模型的截断误差影响,但不同模型扩展时的最优拼接阶次不同,其中DIR-1、DIR-5模型对应的最优拼接阶次分别为180阶和220阶,以GPS水准数据检验,扩展模型在中国区域的精度均优于18cm.最后,基于最优拼接阶次获得的扩展重力场模型对我国1985高程基准重力位进行了估计,DIR-5和TIM-5模型对应数值分别为62636853.47m~2·s~(-2)和62636853.49m~2·s~(-2),精度均为1.51m~2·s~(-2);发现在中国区域模型大地水准面与GPS/水准数据的差值存在微弱的系统性倾斜,东西向倾斜约为9cm,南北向倾斜约为1.4cm,考虑倾斜改正后基于DIR-5和TIM-5模型估计我国1985高程基准重力位的精度提高了0.16m~2·s~(-2).     

Gravity field contribution analysis of GOCE gravitational gradient components

A gravity field model is computed from the four accurate gravitational gradient components of GOCE (Gravity field and steady-state Ocean Circulation Explorer), combined with the analysis of the kinematic orbits, and some moderate constraint (or stabilization) in the polar areas where no observation from GOCE is available due to the orbit geometry. The normal matrix of each component is computed individually in order to study its contribution to the combined solution. The results show that the contribution of V_zz is the largest, with an average value of 32.74% of the total solution; the second and the third largest are V_zz and V_yy, with average values of 28.04% and 26.08%, respectively; the component V_xz contributes 11.81%. Validation with external data shows that each component has its characteristic value and that the information content of the component V_xz is not negligible and should be included for gravity field recovery. The orbit part as derived from high-low satellite-to-satellite tracking (SST-hl) to the GPS contributes mostly to the coefficients below degree and order (d/o) 20, and to non-zonal coefficients from d/o 20 to 80. The mean value of the contribution of the polar stabilization is the smallest with a value of 0.22%, nevertheless it is important. In addition to the contribution analysis in terms of the normal matrices, each individual component of the gradiometer has been combined with SST and polar stabilization, to give a set of single component gravity field models. These partially combined solutions are compared to the fully combined solution in terms of geoid differences. They show that the partially combined solution with V_zz is closest to the complete solution. Even closer is a combination with V_xx and V_yy. In addition to the GOCE-only solution, a GOCE-GRACE (Gravity Recovery And Climate Experiment) combined gravity field model is derived and the information content of GOCE and an available set of normal equations of GRACE are investigated. Results show that, as expected, GRACE dominates the solution below degree 90 and GOCE above degree 140.     

GOCE, Satellite Gravimetry and Antarctic Mass Transports

In 2009 the European Space Agency satellite mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) was launched. Its objectives are the precise and detailed determination of the Earth’s gravity field and geoid. Its core instrument, a three axis gravitational gradiometer, measures the gravity gradient components V _xx , V _yy , V _zz and V _xz (second-order derivatives of the gravity potential V) with high precision and V _xy , V _yz with low precision, all in the instrument reference frame. The long wavelength gravity field is recovered from the orbit, measured by GPS (Global Positioning System). Characteristic elements of the mission are precise star tracking, a Sun-synchronous and very low (260 km) orbit, angular control by magnetic torquing and an extremely stiff and thermally stable instrument environment. GOCE is complementary to GRACE (Gravity Recovery and Climate Experiment), another satellite gravity mission, launched in 2002. While GRACE is designed to measure temporal gravity variations, albeit with limited spatial resolution, GOCE is aiming at maximum spatial resolution, at the expense of accuracy at large spatial scales. Thus, GOCE will not provide temporal variations but is tailored to the recovery of the fine scales of the stationary field. GRACE is very successful in delivering time series of large-scale mass changes of the Antarctic ice sheet, among other things. Currently, emphasis of respective GRACE analyses is on regional refinement and on changes of temporal trends. One of the challenges is the separation of ice mass changes from glacial isostatic adjustment. Already from a few months of GOCE data, detailed gravity gradients can be recovered. They are presented here for the area of Antarctica. As one application, GOCE gravity gradients are an important addition to the sparse gravity data of Antarctica. They will help studies of the crustal and lithospheric field. A second area of application is ocean circulation. The geoid surface from the gravity field model GOCO01S allows us now to generate rather detailed maps of the mean dynamic ocean topography and of geostrophic flow velocities in the region of the Antarctic Circumpolar Current.     

融合GOCE和GRACE卫星数据的无约束重力场模型Tongji-GOGR2019S

本文在法方程层面融合GOCE卫星的V_xx、V_yy、V_zz和V_xz重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用ⅡR带通滤波器在5~41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明：引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型.     

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.     

Band-Limited Analysis of Current Satellite-to-Satellite Tracking, Gradiometry and Combined Earth Gravity Models

Several satellite-only gravity models based on the analysis of satellite-to-satellite tracking (SST) data have become available in the course of the last decade. The realization of the satellite missions CHAllenging Minisatellite Payload (CHAMP) and Gravity Recovery And Climate Experiment (GRACE) enabled the practical implementation of two modes of the SST principle, namely the high–low and the low–low SST. Though similar in their fundamental idea, which is the indirect observation of the gravity field based on the position of two satellites orbiting the Earth, the different architecture and geometrical layout of these techniques capture different fingerprints of the observed field. In the last few years, satellite-only gravity models based on the analysis of satellite gravity gradiometry (SGG) data became available and led to a new insight into the gravity field. The implementation of the SGG principle became possible after the launch of Gravity field and steady-state Ocean Circulation Explorer (GOCE), the first gravitational gradiometry mission. Based on the principle of differential accelerometry, GOCE provides the gravitational gradients which can be used in gravity field retrieval as primary observations of the field at satellite altitude. In the present study, we consider some of the current satellite-only and combined gravity models based on the analysis of CHAMP, GRACE, GOCE, gravimetry and altimetry data. In order to perform a thorough analysis of the models, we present an overview of tools for their quality assessment both in an absolute and relative sense in terms of computing spectral quantities, such as correlation or smoothing coefficients per degree and per order, attempting to demonstrate possible non-isotropic features in the models. Furthermore, typical geodetic measures in computing second-order derivatives, such as degree and order variances and difference variances, have been also evaluated for the same models, using the combined model EGM2008 as reference. Apart from these standard spectral assessment quantities, a systematic spatial representation of the second derivatives at satellite altitude has been performed. The combination of the two analysis steps (spectral and spatial) permits a first detailed assessment of the models, focusing especially on the identification of characteristic interpretable bandwidths.