孙和平1,2 周江存1 徐建桥1 陈晓东1 雷湘鄂1崔小明1 牛晓伟3
1 中国科学院精密测量科学与技术创新研究院 大地测量与地球动力学国家重点实验室
2 中国科学院大学
3 中国科学院精密测量科学与技术创新研究院 武汉大地测量国家野外科学观测研究站
1985 年以来,武汉大地测量国家野外科学观测研究站积累了近 35 年的高精度重力观测数据,为国家精密测绘、全球大尺度地球动力学和空间技术提供新的理论指导和背景材料。
地球重力场的时间变化综合反映了来自其他天体的引潮力变化和地球系统内不同圈层质量迁移情况。而这些质量在空间分布的调整又与地球自转、地表负荷、地震甚至地球内核运动等地球动力学现象密切相关。地球是由地表流体层、地壳、地幔、液态外核和固态内核等组成的分层旋转椭球体,在不同圈层内部和各圈层之间耦合的动力学过程中会产生多种自由运动简正模。它们是与地球深内部结构物性参数和内部动力学等地球科学基础研究领域的前沿热点问题密切相关的,并且会在地表及外部空间的重力场变化中得以体现。因此,对重力场的时间变化特征的分析与研究有助于加深我们对地球各圈层内部活动及各圈层之间的耦合等现象的认识。
目前,观测精度达到 10-8 m/s2(微伽)的重力仪就属于高精度重力仪的范畴。而超导重力仪的精度要比微伽标准高 1—2 个量级,具有超高的观测精度、极低的漂移率、非常高的稳定性和灵敏度。仪器主要部件封装在绝缘的充满液氦的杜瓦瓶内,一个超导小球悬浮在磁场中,磁力正好精确地平衡重力。对上部和下部线圈的电流的比例进行调整,可以使磁力梯度(等同于弹簧常数)非常微弱。因此,非常小的重力变化,会引起超导球大的移动,从而具有非常高的灵敏度。零电阻的超导性能,允许产生磁场的电流可以稳恒流动,只要保持在临界温度的超导状态下,就永远没有电阻带来的损耗,因此具有极强的稳定性。采用磁悬浮测试感应体,从根本上解决了机械弹簧重力仪的零漂问题。因此,超导重力仪是目前性能最好的相对重力仪,是研究重力固体潮及重力长期变化最理想的观测仪器,能有效检测地球表面和内部物质迁移,以及深内部动力学效应。
武汉大地测量国家野外科学观测研究站(以下简称“武汉大地测量站”)超导重力观测始于 1985 年,其超导重力仪在 1997 年经过了仪器厂家的升级更新。后来与比利时、英国和德国等国的重力仪进行了长期的同址比对观测,建立了武汉国际重力潮汐基准,对早期建立的基准值进行了校正,确定了更为精密的潮汐模型,为空间和地面大地测量以及地球物理观测提供精密的潮汐改正,它也是亚洲唯一的国际重力潮汐基准。为研究青藏高原地下结构及隆升等区域构造动力学现象,武汉大地测量站在与武汉基本处于同一纬度区的云南丽江和西藏拉萨分别安装了新型超导重力仪,三地构成一条超导重力长期连续观测的东西链条,为相关研究提供了重要的数据保障。
在对超导重力观测资料的长期分析与研究过程中,也促进了相关理论的发展。得益于全球地球动力学计划(GGP)的实施,基于全球超导重力仪观测资料共享,台站科研人员开展了重力固体潮、大气负荷和海潮与陆地水负荷、地球自由核章动、地球自由振荡、内核平动振荡等全球动力学信号的研究与检测;同时,还开展了地下水重力效应的理论模拟研究。
1
建立了武汉大地测量站的背景噪声模型
台站背景噪声水平是评价观测数据质量和台站观测环境优良水平的一个关键指标。武汉大地测量站分频段(如地震、亚地震、潮汐频段)定量研究了超导重力仪台站的背景噪声水平,发现在低频部分超导重力仪比传统地震仪的背景噪声水平要低,同时也证实了超导台站具有较低的噪声水平、较高的观测质量。定量评估了“中国大陆构造环境监测网络”30个 gPhone 连续重力台站背景噪声水平,探讨了背景噪声的地域分布特征,初步建立了中国大陆重力场背景噪声模型。
2
构建了精密潮汐改正模型
武汉大地测量站利用全球长期连续超导重力观测资料,获得了精密的重力潮汐因子和相位滞后,构建了全球固体潮实验模型。通过对比实验模型与实测结果,发现该实验模型中 O1、K1、M2 和 S2 4 个主要潮波的差异平均优于 0.2%。该实验模型接近于非弹性、非流体静压力平衡初始状态的地球模型的结果。其中,周日潮波部分与 Mathews的结果比较一致,而半日波部分与 Dehant 等的结果更加一致。
利用不同的理论固体潮模型和上述的实验固体潮模型,结合由 NAO99b 全球海潮模型和中国近海海潮模型,基于格林函数方法计算的海潮负荷效应,构建了中国大陆精密重力潮汐改正模型。结果表明,采用不同的固体潮模型或实验模型会对重力潮汐结果产生相对变化幅度小于 0.06% 的差异;在沿海地区,海潮负荷的影响约为整个潮汐的 4%,而中部地区约为1%;其中,中国近海潮汐模型的影响约占整个海潮负荷的 10%,内插或外推潮波的负荷约占海潮负荷的3%。通过与超导重力观测资料进行比较表明,重力潮汐改正模型的精度远远优于 5×10-9 m/s2,这说明构建的改正模型具有很好的实用性,可为中国大陆高精度重力测量提供精密的改正模型。
3
改善了地表质量负荷重力效应的理论模拟
武汉大地测量站利用标准大气模型构造了大气重力格林函数,克服了传统的计算大气负荷效应时模型复杂、计算量大的缺点,获得了精度较高的大气负荷效应。根据大气重力信号随区域分布而变化的特征,计算了大气重力导纳值。基于弹性地球表面负荷理论,利用准三维大气模型,建立了一套大气负荷效应的计算系统。采用地表气压和温度资料计算了大气负荷效应和大气重力导纳值,为快速计算大气负荷效应奠定了基础。
采用由高空间分辨率的数字高程模型确定的海陆边界来模拟沿海和岛屿的海潮负荷效应,构建了我国沿海和岛屿的重力潮汐模型。显示了考虑高精度海陆边界信息在计算沿海和岛屿海潮负荷的重要性,同时也指出准确的高程值在计算沿海重力海潮负荷中的重要性。
利用全球 20 个台站高精度重力潮汐观测序列,深入研究了目前海潮和固体潮模型的适定性及海潮负荷改正的有效性。顾及不同潮波振幅特征,提出了计算各台站平均观测残差和剩余残差矢量的“非等权均值法”,获得了经海潮负荷改正后,全球台站主波平均振幅因子与理论模型间的差异小于 0.3% 和仪器标定误差不超过 0.5% 的重要结论;用地表重力实测数据证实了 Mathews 等关于周日潮 K1 和 O1 波相位滞后是正值的理论计算结果。
利用全球陆地水储量变化模型,基于格林函数方法,模拟了全球陆地水储量变化对重力场的影响;基于重力位理论,联合地表绝对重力测量和空间重力测量结果反演了区域的陆地水储量变化,并用井水位资料验证了反演结果的正确性。从地下水渗透过程的物理机制出发,采用一维水动力学模拟,利用井水位和降雨数据模拟了武汉大地测量站附近的土壤含水率变化,在此基础上估计了地下水变化导致的重力效应 。
4
探测了与地球动力学相关的信号
自由核章动
地球自由核章动(FCN)是由于地幔与液态外核的瞬时旋转轴不一致而在椭球形核幔边界产生的一种自转简正模,它与液核的动力学扁率、核幔边界附近的黏滞系数、电导率、磁感应强度等密切相关。
武汉大地测量站利用全球超导重力观测数据,确定了 FCN 的本征周期为 430 天左右,精密测定了 FCN 的共振参数。估算了核幔边界的黏滞系数和电磁耦合参数,数值结果表明液核顶部的黏滞系数应该在 6.6×102—2.6×103 Pa·s之间,这与根据地球章动、液核自由章动及日长变化等的实际观测得到的黏滞系数结果非常吻合。结果还表明,地幔底部的电导率需要达到 2.6×106—1.0×107 S/m 才能符合实际观测的 FCN品质因子量级,耗散耦合对 FCN 本征周期的影响仅为 1—2 恒星日。
利用模拟数据研究了周日潮波对 FCN 周期变化的敏感程度,根据 K1 和 ψ1 波的变化幅度发现 FCN 周期的变化幅度为 2.5—4 恒星日,与实际甚长基线干涉(VLBI)的观测结果非常符合,进一步验证了 FCN的时变特征。在此基础之上,研究了 FCN 周期的变化,成功提取到高信噪比的地球液态地核动力学信号并发现与核幔边界地磁急变的相关性(图 1),从而说明地磁急变可通过核幔边界电磁耦合影响液核与地幔相对运动。
图 1 地磁急变时间和 FCN 周期变化(实线)关系
自由震荡
1. 提出了解算液核共振参数的三频谱线法。成功检测到 2001 年秘鲁里氏(Ms)8.2 级地震激发的从 0S0 到 0S48 的基频震型,并且观测到了 0S2 和 0S3 的谱峰分裂现象。观测到秘鲁地震与阿拉斯加地震激发的0S2 振型周期的差异,该结果有可能反映了地球内核各向异性,发现自转方向和逆自转方向上 1S2 振型谱线分裂不对称因子之间存在较大差异(图 2)。
2. 检测到 2004 年 12 月 26 日苏门答腊—安达曼特大地震激发的球型基频振型和 2 个径向振型和环型自由振荡的耦合现象,以及谱峰分裂现象。通过自由震荡本征模和谱峰分裂参数,有效地约束了该特大地震的破裂机制(图 2)。
图 2 由武汉超导重力仪检测到的地球简正模谱峰分裂
3. 检测到 2—5 mHz 的背景自由振荡信号,分析了其季节性和年际性变化特征,并结合海洋和大气数据对其激发机制进行了定性探讨。提取了 2—7 mHz 频段范围的长周期面波和背景自由振荡信号,并提取到面波群速度频散曲线,证明了超导重力仪背景噪声观测资料具有反演上地幔结构的潜力。
内核平动振荡
探测了内核平动振荡 Slichter 模三重谱线存在的可能性。结果表明,对于超导重力仪来说,目前的精度条件下有能力探测到 Slichter 模。发现了一组可能来自于内核平动震荡的三重谱线,其周期分别为 5.310、4.995 和 4.344 h。尽管如此,Slichter 模是否存在仍然是一个争议的话题,需要进行更进一步的深入研究。因此,在深内部结构如何影响内核平动振荡的本征周期方面,也开展了相关的理论模拟工作,相关结果可为 Slichter 模的研究、探测及其对地球深内部结构的约束提供理论依据。
5
促进了地球形变理论及数值模拟的发展
固体潮、地表质量负荷等效应在研究地球动力学现象过程中需要作为干扰信号加以扣除,这涉及到地球变形的理论模拟。武汉大地测量站研究人员发展了地球受不同力源作用下的理论变形模拟,开展了模拟方法的研究。改进了谱元法计算负荷和位错勒夫数 ,提出了计算位错勒夫数的解析解方法。发现改正后的残差中有明显与太阳加热相关的周日信号。
特别是,最近发现青藏高原拉萨站的潮汐因子与云南丽江和湖北武汉的潮汐因子之间存在非常明显的差异。该差异是否与青藏高原非常大的昼夜温差引起的热负荷有关?这需要进一步加以澄清。因此,我们顾及球形地球的分层结构,考虑地球的自引力效应、地球内部介质及热力学参数的各向异性,建立了球型分层地球在地表热负荷作用下的形变理论,计算了表征地球变形的热负荷勒夫数。这些新的数值计算方法和理论将会促进超导重力观测研究,特别是观测信号的地球动力学解释。
6
结语
35 多年来,武汉大地测量站研究人员在利用超导重力仪观测技术研究相关科学问题方面,取得了一大批原创性的成果。例如:完善和发展了固体潮理论、海潮负荷计算理论,建立了地表热负荷理论,为相关的研究提供了理论指导;构建了全球固体潮实验模型和中国大陆精密重力潮汐改正模型,为基于地表的大地测量和地球物理的观测提供了必要的潮汐改正模型;对目前国际上的一些热点和前沿问题,如自由核章动、核幔耦合、地磁急变、内核平动振荡、地球简正模及谱峰分裂现象等,进行了深入的研究,得到了一些具有显示度的成果。
上述成果先后获得国家自然科学基金重点项目、中国科学院知识创新工程项目、国家自然科学基金创新群体项目、国家重点研发计划、中国科学院战略性先导科技专项及其他项目的持续资助。获得了测绘科技进步奖一等奖、湖北省自然科学奖一、二等奖。在国际和国内重要期刊上发表了超过百篇高水平研究论文(含 Science 论文 1 篇)。
孙和平 中国科学院院士。中国科学院精密测量科学与技术创新研究院研究员,中国科学院大学教授。1995年毕业于比利时鲁汶大学地球物理专业,获博士学位。长期从事地球动力学与微小形变的高精度重力信号检测、理论模拟和力学机制解释等工作。国际大地测量协会(IAG)会士,国际IAG固体潮与地球动力学委员会委员。曾任中国科学院测量与地球物理研究所所长,中国测绘学会与湖北省测绘学会副理事长。国家杰出青年科学基金获得者,国家自然科学基金委员会创新群体项目与重点项目负责人。Journal of Geodesy、Earth and Planetary Physics,以及《地球物理学报》《测绘学报》和《地震学报》等杂志编委。获 2007 年全国测绘科技进步奖一等奖和 2012 年湖北省自然科学奖一等奖。
文章源自:
孙和平,周江存,徐建桥,陈晓东,雷湘鄂,崔小明,牛晓伟.高精度超导重力观测与研究为国家精密测绘和全球地球动力学提供理论基础.中国科学院院刊,2021,36(2):216-223.
总监制:杨柳春
责任编辑:张帆
助理编辑、校对:PAN
排版:曹艺璇
联系客服