利用卫星测高、GRACE和GOCE资料估计全球海洋表面地转流

重力恢复和气候试验GRACE(gravity recovery and climate experiment)卫星极大地提高了地球重力场的精度和分辨率,特别是中长波分量,联合卫星测高数据可获得全球海洋表面大尺度洋流循环。另外,新一代地球重力和海洋环流探测卫星GOCE(gravity field and steady-state ocean circulation explorer)于2009年3月成功发射,采用卫星重力梯度测量原理,对重力场的高频部分非常敏感,使其高分辨率监测全球海洋循环成为可能。本文利用1~7年GRACE观测数据确定的重力场模型和18个月GOCE观测数据确定的地球重力场模型GO_CONS_GCF_2_TIM_R3,联合卫星测高确定的平均海面高模型MSS_CNES_CLS_11,分别估计全球海洋表面地转流,并且与实测浮标数据结果进行比较。分析表明GOCE重力卫星确定的重力场模型具有更高的空间分辨率,能够确定高精度和高空间分辨率的全球海洋地转流,如墨西哥湾暖流的细节和特征,并且与实测浮标结果基本一致。而基于1~4年GRACE观测资料的模型不能很好估计全球地转流特征,基于7年GRACE观测资料的重力场模型ITG-Grace2010s确定的全球地转流的精度仍低于18个月GOCE观测数据确定的地球重力场模型GO_CONS_GCF_2_TIM_R3的结果,估计的全球地转流仍含有较大的噪声,不能很好地反应中小尺度地转流细节特征。并计算ITG_Grace2010s和GOCE_TIM3的稳态海面地形和全球平均地转流的内符合精度,结果显示,在全球范围内,GOCE_TIM3的稳态海面地形和全球平均地转流的精度都比ITG_Grace2010s结果的精度有着很大的改善,其中ITG_Grace2010s的稳态海面地形的精度为21.6cm,而GOCE_TIM3的结果则为7.45cm,ITG_Grace2010s的全球平均地转流的精度为40.7cm/s,而GOCE_TIM3的结果则为19.6cm/s。     

DGS AT1M-3海洋重力仪的应用及精度评估

为了评估DGS AT1M-3海洋重力仪的精度,首先通过分割的有效重力测线进行重力异常质量的内符合精度评价,得出该海洋重力仪测量精度符合海洋调查规范要求。然后与Sandwell v23测高卫星重力异常进行比对分析,对重力异常质量进行外符合精度评价,可看出DGS测量异常与Sandwell卫星测高异常在整体变化趋势上基本一致。最后以西南印度洋断桥热液区为例,简要阐述DGS所测重力数据的有效性。通过评估,可知该海洋重力仪在动态环境下的工作性能较好,测量精度较高,为今后的海洋矿产资源勘探提供新的测量工具。     

利用卫星测高垂线偏差计算其他重力场元的特性分析

利用卫星测高技术确定海洋重力场,垂线偏差数据作为导出观测量在实际工作中被普遍采用。利用物理大地测量边值问题的定义以及扰动位在球面边界条件下的解,给出了由垂线偏差计算大地水准面高、重力异常和扰动重力的公式。分析了不同积分计算公式在重力场阶谱表达形式下对垂线偏差误差的抑制作用,也分析了不同积分核函数的变化特性,得出基本结论:在利用卫星测高数据求解海洋重力场时,当以格网化海面垂线偏差数据计算重力场参数时,求解的大地水准面高的有效性和稳定性优于重力异常和扰动重力。     

卫星测高中的垂线偏差法

卫星测高中的垂线偏差法是当前利用卫星测高技术研究海洋重力场的最优方法,包括利用卫星测高数据计算垂线偏差和利用该垂线偏差确定海洋重力场两部分。研究了Sandwell、Olgiafi、Hwang测高垂线偏差的计算方法和Molodenskii、Hwang利用测高垂线偏差确定海洋重力场的基本原理,分析比较了上述方法的异同,为科学地利用卫星测高资料反演海洋重力场提供理论依据。     

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.     

Retracking of SARAL/AltiKa Radar Altimetry Waveforms for Optimal Gravity Field Recovery

The accuracy of the marine gravity field derived from satellite altimetry depends on dense track spacing as well as high range precision. Here, we investigate the range precision that can be achieved using a new shorter wavelength Ka-band altimeter AltiKa aboard the SARAL spacecraft. We agree with a previous study that found that the range precision given in the SARAL/AltiKa Geophysical Data Records is more precise than that of Ku-band altimeter by a factor of two. Moreover, we show that two-pass retracking can further improve the range precision by a factor of 1.7 with respect to the 40 Hz-retracked data (item of range_40 hz) provided in the Geophysical Data Records. The important conclusion is that a dedicated Ka-band altimeter-mapping mission could substantially improve the global accuracy of the marine gravity field with complete coverage and a track spacing of <6 km achievable in ～1.3 years. This would reveal thousands of uncharted seamounts on the ocean floor as well as important tectonic features such as microplates and abyssal hill fabric.     

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.     

利用EGM2008模型快速构建扰动重力梯度基准图

重力梯度张量是重力位二阶导,相比重力异常能够更好反映局部区域的细节特征。因此重力梯度导航理论上能为惯性导航提供更好的辅助。重力梯度导航的关键技术之一是背景基准图的构建,推导了扰动重力梯度张量与扰动位在局部指北坐标系中的关系式,并基于EGM2008地球重力场模型构建了一块范围的海域扰动重力梯度张量基准图。为了快速构建基准图,选取了合适的勒让德函数,并将每一个梯度张量的计算式改变求和顺序来提高同一纬度圈上的计算点的计算速度。最后利用梯度张量对角线上三元素满足拉普拉斯约束条件的原理验证了所得基准图的正确性。     

基于卫星测高与卫星测温的表层流产品生成方法应用评估

提高海洋表层流的精度与分辨率对于相关应用领域至关重要。研究引入海表面温度(sea surface temperature,SST)信息对高度计导出流场的改进效果,计算了2018年每日无间隔的全球海洋表层流速度。通过增加热量守恒方程约束,引入卫星测温产品,在地转流基础上生成表层流产品,并与现场漂流浮标速度比较,评估了使用多源卫星遥感获取海流的质量。研究表明,本文方法在不改变沿等温线的切向流速分量的前提下,沿等温线的法向流速分量得到了改进,方法能充分利用海表温度信息提取流场特征信息;方法不适用于海表面温度梯度较小区域,在梯度较大区域可获得明显改进;优化流场能更好地刻画海面海流变化,证明结合卫星测高与卫星测温可改善海洋表层流动;存在强烈中尺度活动和热梯度区域,增加SST观测能够弥补测高资料的不足。     

Validation of ERS-1 and High-Resolution Satellite Gravity with in-situ Shipborne Gravity over the Indian Offshore Regions: Accuracies and Implications to Subsurface Modeling

Geoid and gravity anomalies derived from satellite altimetry are gradually gaining importance in marine geoscientific investigations. Keeping this in mind, we have validated ERS-1 (168 day repeat) altimeter data and very high-resolution free-air gravity data sets generated from Seasat, Geosat GM, ERS-1 and TOPEX/POSEIDON altimeters data with in-situ shipborne gravity data of both the Bay of Bengal and the Arabian Sea regions for the purpose of determining the consistencies and deviations. The RMS errors between high resolution satellite and ship gravity data vary from 2.7 to 6.0 mGal, while with ERS-1 data base the errors are as high as 16.5 mGal. We also have generated high resolution satellite gravity maps of different regions over the Indian offshore, which eventually have become much more accurate in extracting finer geological structures like 85° E Ridge, Swatch of no ground, Bombay High in comparison with ERS-1satellite-derived gravity maps. Results from the signal processing related studies over two specific profiles in the eastern and western offshore also clearly show the advantage of high resolution satellite gravity compared to the ERS-1 derived gravity with reference to ship gravity data.     

局部重力场球冠谐分析中的导数计算及应用

在建立局部地区的地球重力场模型中,球冠谐分析是一种有效的方法。在球冠坐标系下建立重力场模型时,仅涉及到重力位的垂向分量,不必进行关于极角的微分运算,所以仅用单个正交基来建模,但若要在球冠坐标系下对重力场作进一步分析,如类似于平面处理的导数运算,来突出地质体的边界,要计算极角的微分,利用单个正交基表示则不易收敛,且误差较大。本文中作者采用两个正交基函数在球冠坐标系下对重力场建模,推导计算了径向导数和曲面导数模,并利用棱柱体正演模型进行检验,结果表明双基函数建模更适合曲面上的位场表示。最后计算了冲绳海槽地区卫星测高重力异常数据的曲面导数,径向导数清晰划分了东海陆架上大构造单元,曲面导数梯度带显示了东海陆架外缘隆起带,也指示了陆架盆地内的一些局部构造。这表明双基函数的球冠谐分析是一种有效的曲面位场建模方法。     

Surface available gravitational potential energy in the world oceans

Satellite altimetry observations, including the upcoming Surface Water and Ocean Topography mission, provide snapshots of the global sea surface high anomaly field. The common practice in analyzing these surface elevation data is to convert them into surface velocity based on the geostrophic approximation. With increasing horizontal resolution in satellite observations, sea surface elevation data will contain many dynamical signals other than the geostrophic velocity. A new physical quantity, th...     

Oceanic gravity by analytical inversion of Hotine's formula

An analytical inversion of the Hotine formula is developed using fast Fourier transform techniques. Detailed mathematical derivations are used to explain the concepts behind the inverse transformation. Three modifications of the analytical inversion of the Hotine formula are compared and tested using both synthetic data from the OSU91A geopotential model and real GEOSAT altimetry data from the Exact Repeat Mission. The stability of this inverse Hotine approach is investigated using simulated data, and numerical tests are done to quantify the stability of this approach. The approach seems to be numerically stable without employing any stabilization technique. Estimated gravity information from GEOSTAT altimetry data is compared to marine gravity data from shipboard measurements in the Orphan Knoll area. The standard deviations and mean values of the differences between satellite and marine gravity disturbances are 8.2 and 2.9 mGal for the planar approximation, 9.2 and 3.7 mGal for the spherical approximation, and 9.5 and 1.9 mGal for the Molodenskii‐like approximation, respectively, indicating that latitude‐dependent errors affect the latter two approximations. Such errors could be eliminated by performing the calculations by the rigorous one‐dimensional fast Fourier transform (FFT) technique, and any data noise could be filtered out by utilizing covariance knowledge about the input geoid undulations and their errors. Simulation studies also showed that the accuracy of the techniques (for all approximations) can reach a root‐mean‐square (RMS) level of only a few mGal when proper treatment of FFT edge effects is employed and a rather wide area of results is disregarded around the edges.     

Wind wave characteristics based on visual Observations and satellite altimetry

Joint analysis of wind wave characteristics derived from the Voluntary Observing Ship data (VOS) and satellite altimetry is presented as the first step of the synthesis of different data sources. Global distributions of significant wave heights and periods along with wind speed are constructed using various techniques and empirical parameterizations. Good qualitative and quantitative agreement of VOS and satellite altimetry is found especially for regions with high spatio-temporal density of observations. The problems and prospects of the further development of the study are discussed in the context of global wave climatology and marine safety.     

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

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

太平洋卫星测高重力场与地球动力学特征

通过多卫星测高数据的综合处理，获得西太平洋卫星测高重力场，进行不同尺度、深度构造动力信息的分离，探讨诸边缘海盆的地球动力学问题。测高大地水准面反映了研究区板块相互作用的特点，其高频成分可以刻画各海盆的构造特征。测高空间重力异常也可刻画陆架构造及盆地分布，由其推算出的海底地形含有大量的海底构造信息。各边缘海盆的莫霍面埋深具有往南变浅的趋势，与菲律宾海各海盆的莫霍面埋深大致相当，说明岛弧两侧的构造动力强度基本相似。大尺度地幔流应力场总体上反映了欧亚板块向东南蠕散和太平洋板块向北西扩张的特点；日本海北侧和南海巽他陆架的中尺度上地幔对流与地幔柱之间有着密切关系，西菲律宾海的上地幔对流强化了日本-琉球-台湾-菲律宾岛弧的活动强度；小尺度地幔流主要限于软流圈层内部，在各海盆分散，而在冲绳海槽和马里亚纳海槽则会聚，可与均衡重力异常类比。还讨论了大、中、小地幔流体系的特点及相互之间的关系，籍以阐明海盆及海槽演化的地球动力学过程。     

Geodetic measurements in the ocean

This paper covers the topic of marine geodesy, its goals, and applications. Specifically discussed are position determination and establishment of geodetic control on the ocean bottom, ocean surface, and subsurface, and the determination of the geoid, a vertical reference surface. The various techniques used in position determination (including satellite, airborne, radio, inertial and acoustic techniques) are assessed in terms of accuracy, coverage, and contribution to the solution of specific problems associated with position and control. The results of several marine geodetic control experiments are presented. Classical techniques for the determination of the geoid are discussed and assessed, as are new techniques such as satellite altimetry. The outlook for marine geodetic measurements in the ocean is outlined in terms of what is being planned or considered for the next decade, and several recommendations are made.     

基于验潮站资料的HY-2A测高数据质量评定

验潮站资料为评定卫星测高海面高度观测值的精度提供了有效途径。基于数据编辑准则筛选出HY-2A数据,通过引入NCEP实时大气压模型,解决了HY-2A卫星任务后期测高数据产品中部分干对流层延迟项和大气逆压校正项不可用的问题。在此基础上,将HY-2A海面高观测值与验潮站数据进行时空匹配,在选取的9个站点上进行了相关系数和标准差计算。结果表明,两者相关系数平均约为0.676 9,最优可达0.898 7,平均标准差为0.089 5 m。整体验证结果表明,HY-2A卫星测高数据质量符合设计指标,达到预期水平,为海洋重力场反演等应用研究提供了新的可靠数据源。     

海洋测绘信息处理新技术

着重对当前海洋测绘信息处理技术的几个研究热点进行了介绍,内容包括海陆基准的统一、海洋遥感图像的精校正、卫星测高和遥感测深技术,并对海洋测量信息处理与多元化表示作了简要的讨论。     

Modification of Globwave satellite altimetry database for sea wave field diagnostics

A new database of ocean wave parameters has been created based on satellite altimetry observations. The basis was data from the European Space Agency project GlobWave (www.globwave.org), which was transformed to suit upcoming requirements for global wave analysis. The new database contains additional wave characteristics (altimetry wind speed estimated using different parametric models, steepness, period, and some quality control parameters). It provides up-to-date tools for mass data preprocessing. The new database makes it possible to optimize wave field diagnostics on regional and global scales. Using the Envisat and Jason-1 satellite missions as an the example, we demonstrate the specific features of using the initial GlobalWave data set and the modified database.