地球物理大地测量

大型科学

计算平台

综合多种异质离散残差观测扰动场元（包括扰动重力、高程异常、空间异常、扰动重力梯度或垂线偏差），采用球面径向基函数逼近方法，估计给定计算面上全要素残差场元格网，从而实现重力场及大地水准面的统一建模。    

由离散残差扰动重力（mGal）、高程异常（m）、空间异常（mGal）、扰动重力梯度（径向，E）或垂线偏差（南向/西向，″）中某种单一类型观测数据，选择球面点质量核函数、Poisson核函数、m次径向多级子核或Poisson小波核函数中的一种球面径向基函数SRBF进行重力场逼近，估计大地水准面外部残差扰动重力、残差高程异常、残差空间异常、残差扰动重力梯度或残差垂线偏差。

由地面或地球外部非等位的边界面上空间异常或扰动重力（mGal）格网、边界面模型高程异常格网、边界面大地高格网和参考等位面（Stokes边值归算面）大地高格网，计算边界面上离散点空间异常或扰动重力的Molodensky一阶项边值改正数，从而将Molodensky边值问题，转换为Stokes边值问题。

给定时间跨度和采样间隔，计算地球自转参数（自转极移和日长变化）的长周期和短周期（周日/半日）潮汐效应预报时间序列。此功能调用IERS协议2010中收录的模型和Fortran代码模块实现，这些代码可从IERS网站获得。

由边界面大地高格网（m）及其面上残差扰动场元格网，按Poisson严密数值积分算法，计算大地水准面及其外部空间的残差扰动场元。Poisson积分是第一边值问题解，只要求边界面连续可微，不要求边界面是严格的重力等位面。算法的边缘效应小，可有效抑制短波信号衰减，适合向上向下解析延拓，一次积分性能接近梯度法解析延拓效果。

由等位面大地高格网及其面上残差扰动场元格网，按严密径向梯度积分公式，计算该等位面上计算点的残差场元径向梯度值（/km）。扰动场元梯度积分由Stokes边值问题解导出，要求被积扰动场元位于等位面。 径向方向由地心指向地球外部。

由各种大地测量监测量时间序列，按球面径向基函数SRBF负荷形变场逼近算法，计算地表负荷等效水高(cm)，高程异常(大地水准面mm)、地面重力(μGal)、扰动重力(μGal)、地倾斜(mas)、垂线偏差(mas)、地面水平位移(mm)、大地高(mm)、正(常)高(mm)、扰动重力梯度(mE)与水平重力梯度(E)负荷形变效应格网时间序列。    

由各种大地测量监测量时间序列，以负荷格林函数积分为动力学约束，计算地表负荷等效水高(cm)，高程异常(大地水准面mm)、地面重力(μGal)、扰动重力(μGal)、地倾斜(mas)、垂线偏差(mas)、地面水平位移(mm)、大地高(mm)、正(常)高(mm)、扰动重力梯度(mE)与水平重力梯度(E)负荷形变效应格网时间序列。   

由边界面大地高格网（m）及其面上残差高程异常（m）格网，计算大地水准面或地球外部的残差空间异常（mGal）、残差扰动重力（mGal）与残差垂线偏差向量（″，SW）。（可从别人大地水准面成果中反演高精度局部重力场，偷笑）。程序采用Poisson积分与场元微分组合算法，实现高程异常反运算，不要求边界面是重力等位面。

由等位面大地高格网（m）及其面上残差高程异常（m）或垂线偏差向量（S,W，″）格网，按逆Stokes、Hotine或Vening-Meinesz运算数值积分方法，计算等位面上其他类型残差扰动重力场元。
扰动重力场元逆运算积分属Stokes边值问题，要求被积高程异常或垂线偏差向量（S,W）位于等位面，通常用于海洋卫星测高重力场反演计算。

对全球地面/海面大气压分潮、海洋潮高分潮的调和常数球坐标格网进行球谐分析，生成规格化大气压潮、海潮负荷球谐系数模型（IERS中FES2004格式）和负荷潮一阶项同相异相幅值文件。负荷潮球谐系数模型用于预报任意时刻全球地面大气压潮或海面潮高，预报任意时刻固体地球外部全要素大地测量负荷潮效应。负荷潮一阶项同相异相幅值用于预报地球质心变化负荷潮效应。

对地面/海面大气压变化（hPa）、大陆水变化（cm）和海平面变化（cm）等全球非潮汐负荷球坐标格网时间序列进行球谐分析，生成地表环境负荷地球质心变化（mm）时间序列与规格化负荷球谐系数模型（m）时间序列。负荷球谐系数模型时间序列可进一步用于计算地面及固体地球外部全空间全要素大地测量非潮汐负荷效应时间序列。