利用FFT技术计算大地水准面高若干问题研究

本文首先从普通二维 FFT算法的基本定义和要求出发 ,深入分析和比较了大地水准面二维平面和二维球面 FFT算法的特点和差异 ,找出了影响二维球面 FFT计算精度的主要误差源 ,并给出了二维球面 FFT计算公式的改进形式。在此基础上 ,本文详细分析讨论了参考场选取、积分球冠半径确定以及核函数改化对大地水准面计算结果的作用和影响。最后 ,本文利用全国 5′× 5′实测重力异常对我国海陆大地水准面进行了试算 ,并将计算结果分别同陆上GPS水准和海上卫星测高观测值作了比较。     

The reduction of aliasing in gravity anomalies and geoid heights using digital terrain data

Observations of gravity can be aliased by virtue of the logistics involved in collecting these data in the field. For instance, gravity measurements are often made in more accessible lowland areas where there are roads and tracks, thus omitting areas of higher relief in between. The gravimetric determination of the geoid requires mean terrain-corrected free-air anomalies; however, anomalies based only on the observations in lowland regions are not necessarily representative of the true mean value over the topography. A five-stage approach is taken that uses a digital elevation model, which provides a more accurate representation of the topography than the gravity observation elevations, to reduce the unrepresentative sampling in the gravity observations. When using this approach with the Australian digital elevation model, the terrain-corrected free-air anomalies generated from the Australian gravity data base change by between 77.075 and −84.335 mgal (−0.193 mgal mean and 2.687 mgal standard deviation). Subsequent gravimetric geoid computations are used to illustrate the effect of aliasing in the Australian gravity data upon the geoid. The difference between 'aliased' and 'non-aliased' gravimetric geoid solutions varies by between 0.732 and −1.816 m (−0.058 m mean and 0.122 m standard deviation). Based on these conceptual arguments and numerical results, it is recommended that supplementary digital elevation information be included during the estimation of mean gravity anomalies prior to the computation of a gravimetric geoid model.     

辽宁省精细化暴雨洪涝灾害风险评估与预评估

基于辽宁省61个国家气象站1961—2020年和998个区域自动气象观测站建站至2020年逐小时、逐日降水资料,分析了辽宁省暴雨洪涝灾害主要致灾因子,计算了暴雨洪涝孕灾环境指标,完成了辽宁省暴雨洪涝灾害危险性评估。结果表明：暴雨洪涝高危险性地区主要位于丹东;暴雨洪涝灾害人口高风险区主要位于沈阳和大连市区;经济高风险区主要位于大连和盘锦市区;水稻、玉米高风险区主要位于锦州、盘锦和丹东。利用辽宁省无缝隙智能网格预报数据对2022年7月28—29日的暴雨过程灾害风险进行了预评估,发现暴雨灾害危险性高值区域主要分布在朝阳、葫芦岛以及辽宁中部。暴雨灾害可能造成的人口、经济高风险区域主要位于辽宁西部和中部地区;暴雨灾害可能造成的水稻和玉米高风险区主要位于沈阳、铁岭和朝阳北部等地区。预计高风险区主要影响人口约为449万人,经济损失约为1432万元,受影响的水稻面积约为1.028万公顷、玉米面积约为1.798万公顷。通过灾后效果检验,发现预评估模型效果良好,可在实际的暴雨洪涝灾害风险评估业务中使用。     

海洋重力大地水准面的确定

美国海洋卫星测高仪的出现,使应用Hotine积分确定海洋大地水准面成为现实。本文通过对Hotine积分及垂线偏差的计算公式进行改进,较好地改善了求和项的收敛性,减小了截断误差影响,并提出了利用Hotine函数和重力异常确定海洋大地水准面的方法。 实际计算表明:海洋重力大地水准面的精度在1米以内;卫星测高大地水准面间存在0.5米系统差;它和海底地形有一定的相关性,能较好地反映出海底地形的宏观特性。     

Detection of major geological structures in the Eastern offshore of India using Geosat altimeter data

Satellite altimetry can be used to infer subsurface geological structures analogous to gravity anomaly maps generated through ship-borne survey. The Eastern offshore was taken up for analysis using Geosat Exact Repeat Mission (ERM) altimeter data. A methodology is developed to use altimeter data as an aid to offshore hydrocarbon exploration. Processing of altimeter data involves corrections for various atmospheric and oceanographic effects, stacking and averaging of repeat passes, cross-over correction, removal of deeper earth and bathymetric effects, spectral analysis and conversion into free-air gravity anomaly. The final processed results were derived for Eastern offshore in the form of prospecting geoid and gravity anomaly maps and their spectral components. The highs and lows observed in those maps were derived in terms of a number of prominent megastructures e.g., gravity linears, 85°E and 90°E ridges, the Andaman trench complex etc. Satellite-derived gravity profiles along 12°N latitude match well with the existing structures.     

利用垂线偏差计算大地水准面中央区效应的改进方法

为提高利用Molodensky公式反演测高大地水准面中央区效应的精度,视中央区为矩形域,将垂线偏差分量表示成双二次多项式插值形式,引入非奇异变换,推导出了大地水准面的计算公式。垂线偏差理论模型下的分析表明本文导出公式误差为零,而传统公式的误差与纬度以及垂线偏差子午分量与卯酉分量之间的比值有关;以中纬度区域分辨率为2'*2'的垂线偏差数据为背景场进行了实际计算,结果表明在反演计算点本身所在的1个网格对大地水准面的贡献时,传统公式与本文导出公式计算结果差值的最大值达数厘米。本文导出公式可为测高大地水准面的高精度反演提供理论依据。     

The relationship between gravity and bathymetry in the Pacific Ocean

Summary. Surface-ship and satellite derived data have been compiled in new free-air gravity anomaly, bathymetry and geoid anomaly maps of the Pacific Ocean basin and its margin. The maps are based on smoothed values of the gravity anomaly, bathymetry and geoid interpolated on to a 90 × 90 km grid. Each smoothed value was obtained by Gaussian filtering measurements along individual ship and subsatellite tracks. The resulting maps resolve features in the gravity, bathymetry and geoid with wavelengths that range from a few hundred to a few thousand kilometres. The smoothed values of bathymetry and geoid anomaly have been corrected for age. The resulting maps show the Pacific ocean basin is associated with a number of ENE–WSW-trending geoid anomaly highs with amplitudes of about ± 5 m and wavelengths of about 3000 km. The most prominent of these highs correlate with the Magellan seamounts–Marshall Gilbert Islands–Magellan rise and the Hess rise–Hawaiian ridge regions. The correlation between geoid anomaly and bathymetry cannot be explained by models of static compensation, but is consistent with a model in which the geoid anomaly and bathymetry are supported by some form of dynamic compensation. We suggest that the dynamic compensation, which characterizes oceanic lithosphere older than 80 Myr, is the result of mantle convection on scales that are smaller than the lithospheric plates themselves.