地球物理学报杂志社
分享到:

地球物理学报杂志

《地球物理学报》杂志在全国影响力巨大,创刊于1948年,公开发行的月刊杂志。创刊以来,办刊质量和水平不断提高,主要栏目设置有:学术论文、综述或评述、研究短论、科学论坛、学术动态等。
  • 主管单位:中国科学院
  • 主办单位:中国地球物理学会;中国科学院地质与地球物理研究所
  • 国际刊号:0001-5733
  • 国内刊号:11-2074/P
  • 出版地方:北京
  • 邮发代号:2-571
  • 创刊时间:1948
  • 发行周期:月刊
  • 期刊开本:B5
  • 复合影响因子:1.73
  • 综合影响因子:1.867
相关期刊
服务介绍

地球物理学报 2015年第10期杂志 文档列表

地球物理学报杂志空间物理学·重力与大地测量学

伴随着高速流的磁尾电流片拍动观测研究

摘要:磁尾电流片在磁尾动力学过程中起着重要作用.卫星观测表明磁尾电流片经常处于拍动状态.但磁尾电流片拍动的特性和产生机制至今仍然没有被完全弄清楚.本文主要利用欧洲空间局Cluster卫星数据,研究一个伴随高速离子流的电流片拍动事件.该电流片拍动事件具有很强的周期性.拍动的周期约是2min,磁场振荡幅度约为20nT.能量电子和离子的通量具有周期性增强和减弱的特征.电流密度X和Y分量也具有周期性的振荡,并且振荡周期与磁场振荡周期一致.通过对粒子流速矢量与电流矢量的分析,发现粒子运动具有涡旋的特征.因此可以推断,该磁尾电流片的拍动不是由磁尾等离子体片高速流产生的,而是与局地等离子体不稳定性有关.
3449-3456

外源Sq等效电流体系电流涡中心电流强度与太阳活动关系

摘要:本文选取了INTERMAGNET地磁台网2001年到2012年的地磁数据,对其进行世界时(UT)到地方时(LT)的转换后利用自然正交分量法(NOC)从所选资料中提取出太阳静日变化Sq成分,再通过球谐分析方法建立模型分离内、外源Sq成分,逐日反演出内、外源Sq等效电流体系,并得到外源Sq等效电流体系南北电流涡中心电流强度.本文将外源Sq等效电流体系南北电流涡中心电流强度与同一时期的Dst指数进行了对比分析,研究表明它们之间具有同步变化的规律,且北半球电流涡中心电流强度在磁暴发生时的异常现象远高于南半球.对F10.7cm太阳射电流量与外源Sq等效电流体系南、北半球电流涡中心电流强度的长短周期分析发现,Sq等效电流表现出明显的11年周期特点,与太阳活动周期一致.外源南、北半球电流涡中心电流强度和F10.7cm年均值的相关系数分别达到了0.93和0.90,说明太阳活动是导致外源Sq电流体系变化的最直接也最主要的因素,这可能与电离层电导率受控于太阳的电磁辐射相关.
3457-3468

一种基于LEO卫星信标的电离层层析成像新算法

摘要:LEO卫星信标是电离层监测的重要手段之一.利用电离层层析成像算法,LEO卫星信标能够实现区域电离层电子密度的快速重构.针对LEO卫星信标的特点,本文提出了一种函数基模型与像素基模型组合的电离层层析成像新算法.选择差分相对电离层总电子含量作为输入数据源,先通过函数基模型法获取电离层电子密度初始分布,再利用像素基模型法对初始分布进行二次迭代重构,该方法可有效降低电离层层析成像对背景电离层模型的依赖,同时能够实现电离层小尺度扰动结构的有效反演.利用数值仿真方法及低纬度电离层层析成像网的实测数据的反演结果验证了本文提出的新算法的可行性和可靠性.
3469-3480

利用瑞利激光雷达观测北京地区上平流层地形重力波活动

摘要:本文利用中国科学院空间科学与应用研究中心的瑞利激光雷达首次观测到了平流层地形重力波活动的现象,并结合美国国家环境预报中心(NCEP)的全球预报系统(GFS)的风场数据分析了该地形重力波的基本参数.与惯性重力波相比较,地形重力波的密度扰动没有下传的相位,在同一高度上,其扰动相位保持不变.北京空间科学与应用研究中心瑞利激光雷达自2012年开始观测实验以来,已经观测到多起地形重力波活动事件.本文以2013年11月11日的观测数据为例,研究北京上空的地形重力波活动,并结合GFS风场数据分析了北京上平流层地形重力波的波长、传播方向、传播速度等参量.通过分析得到在2013年11月11日北京上空存在一列传播方向为北偏西52.4°,水平波长为5.5km,平均垂直波长约为6.0km的地形重力波.
3481-3486

模拟分析低低跟踪模式重力卫星反演地球重力场的精度

摘要:本文利用卫星重力反演与模拟软件ANGELS系统(ANalyst of Gravity Estimation with Low-orbit Satellites)对低低跟踪模式的重力卫星的关键载荷精度指标进行了深入分析.模拟结果表明:(1)对短弧长积分法而言,在低低跟踪模式的关键载荷精度指标中,重力场反演精度对星间距离变率精度最为敏感;(2)通过对目前在轨运行GRACE的载荷指标进行分析,发现轨道数据的误差主要影响重力场的低阶部分(约小于25阶),较高阶次部分(约大于26阶)主要受星间距离变率的误差限制;(3)如果下一代低低跟踪模式的重力卫星的目标之一是把重力异常反演精度较GRACE提高约10倍,则在保持轨道高度和GRACE相同的前提下,轨道、星间距离变率和星载加速度计等关键载荷指标需要达到的最低精度分别约为2cm、10nm·s^-1和3.0×10^-10 m·s^-2;(4)轨道精度和混频误差将是影响下一代低低跟踪模式重力卫星重力场恢复能力进一步提高的主要制约因素,距离变率精度和加速度计精度存在盈余.
3487-3495

青藏高原GRACE卫星重力长期变化

摘要:本文采用最新的GRACE(Gravity Recovery and Climate Experiment)(RL05)数据,通过水文模型(Global Land Data Assimilation System,GLDAS与Climate Prediction Center,CPC)扣除土壤水及雪水的影响,利用Paulson提供的冰川模型结果扣除GIA(Glacial Isostatic Adjustment)的影响,采用尺度因子的方法减少数据处理过程中误差的影响,最终基于最小二乘计算方法得到2003—2013中国及周边地区长期性重力异常变化情况.结果发现青藏高原有较为明显的重力上升信号,我们认为该信号可能由印度板块俯冲欧亚板块导致青藏高原地壳增厚所引起.接着依据GPS观测结果和艾黎均衡假说构建了地壳形变模型并通过直立长方体模型予以正演模拟分析.以班公湖—怒江断裂带为界将青藏高原划分为南北两大区块,结果显示青藏高原重力异常大致以0.2μGal·a^-1的速率在递增,小于GRACE得到的0.3±0.08μGal·a^-1的增长速率(对应于地壳增厚速率约3mm·a^-1),剩余未解释部分可能与湖水、冰川因素、冻土因素等有关.该结果对于认识青藏高原隆升动力学有一定参考意义.
3496-3506

基于Forward-Modeling方法的黑河流域水储量变化特征研究

摘要:黑河流域陆地水储量变化对流域下游等周边区域水资源的合理利用以及经济和社会发展等有着重要的意义.本文利用2003年1月至2013年12月的GRACE RL05数据反演了黑河流域陆地水储量长时间序列的变化,并针对重力场模型和数据处理中产生的信号泄漏问题,采用Forward-Modeling方法进行了改正并恢复泄漏信号;将GRACE获得的泄漏信号恢复前后的黑河流域水储量变化结果与全球水文模型GLDAS和CPC进行比较分析,结果表明泄漏信号改正后的结果与水文模型结果的时间序列相关性均有明显提高,从其空间分布结果可以看出Forward-Modeling方法有效地恢复初始信号、增强被湮没的信号,泄漏信号误差减小;通过分析黑河流域水储量变化的长时间序列结果,发现其具有明显的阶段性变化特征,即2003—2006年呈明显下降趋势,约为-0.86cm·a^-1,在2007—2010年趋于平衡状态,而2011—2013年则呈现缓慢上升趋势约为0.14cm·a^-1;联合GRACE数据和GLDAS数据反演了黑河流域地下水储量变化,并与全球降雨数据GPCC进行了比较分析,两者相关性可达到0.88以上.
3507-3516

地球形成和演化过程中的分异能计算方法研究

摘要:地球形成初期,构成地球的物质在组成上是大致均一的.目前地球的地核-地幔-地壳圈层结构,是由分异作用形成的.分异过程释放的能量称为分异能.Sorokhtin和Chilingarian等人从行星吸积的定义出发,导出了基于地球内部密度分布的势能计算公式,计算出的分异能大小为1.698×10^31J.本文采用计算球体势能的思路,导出分异能计算的解析公式和数值计算公式,通过求取原始地球模型与均匀分层模型、PREM模型的势能差计算分异能.两种方法的计算结果分别为1.535×10^31J和1.698×10^31J.前者与Sorokhtin等的结果相近,后者与之相同.本文初步分析了方法间的异同以及造成结果偏差的主要原因.
3530-3539
地球物理学报杂志地球动力学·地震学·地热学·地磁学

上地壳纵横波速度结构相关反演成像方法

摘要:基于纵横波初至走时数据的层析成像方法越来越广泛地被应用于揭示不同构造域壳幔速度结构特征.我们从同一地质体的纵横波速度属性相关这一基本思想出发,提出一种相关反演成像的方法:纵横波速度反演交替进行,在迭代反演过程中每通过一次反演获得相应的纵波速度(或横波速度)结构后,更新相应的纵横波速度比模型以及相应的横波(或纵波)速度反演的初始模型,然后继续开展后续横波(或纵波)速度反演工作.在反演过程中依据纵横波速度的相关性信息和射线路径长度将走时残差以不同权重分配到射线路径经过的单元,依据网格节点周围平均的慢度扰动更新速度模型.正反演过程分别基于有限差分走时计算方法和反投影成像方法.两种典型模型试验表明,该技术应用于上地壳速度结构反演成像过程,可有效提高反演结果的可靠性,在很大程度上避免了常规单独反演纵波和横波速度过程容易带来的畸变和失真.该方法应用于重建青藏高原西部札达—泉水沟深地震测深(DSS)剖面下方的上地壳速度结构,揭示出与青藏高原西缘板块碰撞相关的上地壳速度结构特征.
3553-3570

SsPmp震相地壳探测方法

摘要:SsPmp波是远震S波经地表反射转换的P波在莫霍面发生反射后被地表台站接收得到的震相.震中距在30°~50°之间的远震S波震相经地表反射转换的P波射线参数较大,在莫霍面发生全反射,使得台站接收的SsPmp波具有较强的能量,能够从地震记录中清楚地识别出来,为探测台站附近的莫霍面形态提供新的途径.本文通过合成理论地震图分析了SsPmp震相与地壳厚度、射线参数和Pn波速度之间的关系.结果表明:对于水平界面,地壳厚度只影响SsPmp与Ss波之间的相对到时差;Pn波速度只影响SsPmp的相位;射线参数既对SsPmp波的相对到时有影响,也会引起SsPmp波的相位变化.对于复杂的界面,SsPmp反映的深度与速度梯度最大的深度接近,而反映的Pn波速度与实际的Pn波速度一致.
3571-3582

用接收函数建立区域模型的震源机制反演及其在芦山地震序列研究中的应用

摘要:本文提出并试验了一种基于接收函数建立区域模型进行震源机制反演的方法.选取四川地震台网记录的M≥3且信噪比高的近震波形资料,反演得到了芦山地震序列中74个地震的震源机制.通过对震源深度和震源机制的综合分析,探讨了芦山地震的发震构造和区域应力场状态.采用接收函数方法反演获取了26个台站下方的S波速度结构,对不同区域的台站反演结果进行叠加平均,以此区域平均S波速度作为本文震源机制反演使用的区域模型的S波速度;区域模型的P波速度由经验公式给出.反演稳定性测试表明,使用不同模型或对原始波形记录加入随机噪声的反演结果与原始反演相比,震源深度最大误差为1km,断层面各参数误差水平也很低,且显示的发震类型是一致的,其中随机噪声带来的误差小于模型带来的误差.主震反演得到的震源机制解为:震源深度17km,矩震级6.47;节面Ⅰ走向213°,倾角51°,滑动角98°;节面Ⅱ走向20°,倾角40°,滑动角80°;显示芦山主震可视为纯逆冲型地震,发震构造可能是某个具有较大倾角的逆冲断层,而不是低缓的推覆构造的基底滑脱面.同时本文反演获取的73个M≥3余震的震源机制绝大多数也显示了类似的发震类型,逆冲型地震为67个,占92%,具有绝对优势;走滑型地震为5个,正断型地震为1个.其中5个走滑型地震中的4个均分布在震源区的东北端.整个芦山地震序列深度集中在12~20km,且沿震源区短轴的余震深度剖面有自西向东呈逐步变浅的趋势,呈现清晰的铲形断面结构,结合本地地质构造,可以推断芦山地震序列主要发生在龙门山前山断裂以东的逆冲推覆体内的一个隐伏断裂上.P轴方位角优势方位与区域应力场及汶川震源区南段的相一致,表明芦山序列地震活动主要受区域应力场控制,且汶川震后该区应该不存在应力场变化.P轴仰角随深度分布则显示了孕�
3583-3600

一种地方与区域地震震源机制反演技术:广义极性振幅技术(二)——对实际震例的应用

摘要:为了检验广义极性振幅技术(GPAT)的实用性,我们利用GPAT反演确定了49次实际地震的震源机制、矩震级和震源深度.为了检验GPAT对地方地震、区域地震以及远震的实用性,我们选用了震级范围约为ML0.2~MS7.0之间的地震与震中距范围约在5~8000km之间的观测资料.对反演结果的分析表明,利用GPAT获取的震源机制结果在合理误差范围内是正确的,利用GPAT获取的矩震级结果是可靠的,而利用GPAT获取的震源深度总体上似乎比常规定位深度深约0.6km.总体而言,GPAT在震源机制、矩震级和震源深度的反演方面表现出良好的实用性.
3601-3614

拉分盆地形成机制三维数值模拟——以海原断裂带老龙湾盆地为例

摘要:拉分盆地是走滑断层系中受拉伸作用形成的断陷盆地.一般在两条平行断层控制下发育.盆地形似菱形,几何形态主要受两条主控走滑断层错距和叠接长度影响.本文以青藏高原东北缘海原断裂带老龙湾拉分盆地第四纪所处的构造环境为基础,参考盆地周围断层几何分布,建立了三维有限元数值模型,模拟该拉分盆地的演化过程;进一步分析了断层力学性质、地壳分层结构等各因素对盆地形成和演化的影响.模拟结果显示,盆地地表沉降伴随有下地壳物质的上涌,此上涌对盆地地表沉降存在阻碍作用.各因素的影响具体表现为:(1)断层力学性质(弹性模量和黏滞系数)越弱,其对构造应力较低的传递效率导致盆地两端差异性运动越明显,从而形成较大的盆地地表沉降和明显的上地壳减薄.(2)平行主控断层的叠接长度反映盆地形成的拉伸作用范围,叠接长度越大,相同的差异性运动在单位面积形成的拉伸应力越小,盆地地表沉降较小.(3)下地壳流变性影响其物质的上涌量,下地壳黏滞系数越小,其对上部拉伸作用的响应越明显,上涌量越大,此上涌对上地壳沉降形成的阻碍作用也越明显.根据老龙湾拉分盆地所处的构造格局,将平行断层的叠接长度取20km,当断层黏滞系数取值为周围基岩的1/10,参考该盆地第四纪构造演化历史,模拟得到的盆地第四纪下沉量与盆地内第四系沉积层厚度在规模上近似,下地壳黏滞系数取值在(2.5~5.0)×10^21 Pa·s范围内时,盆地下沉量模拟结果与老龙湾拉分盆地第四系地层厚度吻合较好.
3615-3626

华北地区水位下降是否会减缓气温上升——浅部地温影响的数值模拟分析

摘要:华北地区由于长期持续的地下水过量开采,导致了大面积地下水位大幅下降,引发地面塌陷、地下水质污染等一系列地质环境问题,这些现象早已为人们所熟知和关注.然而地下水位下降还会造成百米量级浅部地温及其梯度的变化,因此即使来自地球深部的大地热流密度没有变化,年度平均的从表浅部位通过地表实际传导进入大气的热流密度会减小,这是中外文献中尚未见讨论过的问题.我们通过数值模拟发现假定大地热流密度不变的条件下,华北数万平方公里地下水位下降会造成百米尺度内的地温降低,从而传入大气的热流密度降低40%以上,且会持续数百年以上的时间.这种长时间大范围的传导入大气的热流密度变化对环境会造成什么影响是一个十分值得关注的问题.这一预测在一定程度上得到了气象站地温观测数据的支持,但由于目前气象观测站只有3.2m深度范围内的地温资料,累计不超过5、60年,中间还有10余年的间断,而且表浅深度地温受地表多种因素的影响也较大,这些资料难以对我们关心的地下水位下降引起流入大气的热流密度变化这一问题提供直接确凿的数据来进行分析,因此今后有必要开展对地下数十乃至数百米地温进行持续精确的监测工作.
3649-3659

川西坳陷鸭子河地区基于多种古温标的钻井热史恢复

摘要:利用镜质体反射率(Ro)、磷灰石裂变径迹(AFT)和伊利石结晶度(IC指数)等古温标恢复了四川盆地川西坳陷的钻井热史,对比了不同温标最高古地温的恢复结果.研究表明,研究区晚白垩世至今总体表现为冷却及抬升剥蚀的过程,地温梯度由约26℃·km^-1降低至约22℃·km^-1,剥蚀量约1.3~1.9km.约80 Ma以来开始抬升剥蚀,40—2.5 Ma经历了一个热平静期,第四纪存在一定的增温,地温梯度增高约5℃·km^-1.三种古地温恢复结果具有较高的一致性,相对于镜质体反射率(Ro)和磷灰石裂变径迹(AFT)等成熟古温标,伊利石结晶度作为有机质成熟度指标和沉积岩古温标的应用处于定性分析阶段,该指标的热演化模型仍需进一步探索.
3660-3670

鄂尔多斯盆地深部热结构特征及其对华北克拉通破坏的启示

摘要:基于二维稳态热传导方程,利用有限元数值模拟方法,选取东西向横穿鄂尔多斯盆地地质与地球物理解释大剖面进行了深部温度场数值模拟研究,得到了华北克拉通西部的鄂尔多斯盆地下伏岩石圈热结构特征.地幔热流变化范围:21.2~24.5mW·m^-2,体现为东高西低特征.壳幔热流比(Qc/Qm)介于1.51~1.84之间,为"热壳冷幔".与华北东部地幔热流对比表明,西部的鄂尔多斯盆地相对处于稳定的深部动力学环境.在岩石圈热结构研究基础上,对克拉通地震岩石圈与热岩石圈厚度差异进行了对比,研究表明:鄂尔多斯盆地西部地震岩石圈与热岩石圈厚度差异约达140km,而东部的汾渭地堑,渤海湾盆地二者差异逐渐减小.华北克拉通自西向东,地震岩石圈厚度与热岩石圈厚度差异不断减小,意味着华北克拉通岩石圈下部的软流圈地幔黏性系数自西向东逐渐降低,本文从地热学角度可能印证了太平洋俯冲脱水作用对华北克拉通的影响.
3671-3686

柴达木盆地现今大地热流与晚古生代以来构造-热演化

摘要:依据钻孔系统稳态测温、静井温度资料与实测热导率数据分析了柴达木盆地地温场分布特征,建立了柴达木盆地热导率柱,新增了17个大地热流数据.柴达木盆地现今地温梯度介于17.1~38.6℃·km^-1,平均为28.6±4.6℃·km^-1,大地热流介于32.9~70.4mW·m^-2,平均55.1±7.9mW·m^-2.盆地不同构造单元地温场存在差异,昆北逆冲带、一里坪坳陷属于"高温区",祁南逆冲带属于"中温区",三湖坳陷、德令哈坳陷及欧龙布鲁克隆起属于"低温区",盆地现今地温场分布特征受控于地壳深部结构、盆地构造等因素.以现今地温场为基础,采用磷灰石、锆石裂变径迹年龄分布特征定性分析与径迹长度分布数据定量模拟相结合,研究了柴达木盆地晚古生代以来的沉积埋藏、抬升剥蚀和热演化史,并结合区域构造背景,对柴达木盆地构造演化过程进行了探讨,研究表明柴达木盆地晚古生代以来经历了六期(254.0—199 Ma,177—148.6 Ma,87—62 Ma,41.1—33.6 Ma,9.6—7.1 Ma,2.9—1.8 Ma)构造运动,六期构造事件与研究区构造演化的动力学背景相吻合.其中白垩纪末期(87—62 Ma)的构造事件导致了柴达木盆地东部隆升并遭受剥蚀,欧龙布鲁克隆起形成雏形,柴达木盆地北缘在弱挤压环境下形成坳陷盆地;中新世末的两期构造事件(9.6—7.1 Ma和2.9—1.8 Ma)使柴达木盆地遭受强烈挤压,盆地快速隆升,构造变形强烈,基本形成现今的构造面貌.
3687-3705

华南二叠系卡匹敦阶高分辨率浮点年代标尺

摘要:基于详细的生物地层学研究,以磁化率为古气候替代指标,对广西来宾铁桥剖面卡匹敦阶地层开展时间序列分析,建立高分辨率浮点年代标尺(FPTS).结果表明,磁化率记录了铁桥剖面中二叠世晚期沉积序列中的米兰科维奇旋回,卡匹敦阶上部磁化率突然增加与峨眉山玄武岩喷发和卡匹敦晚期全球性海退有关,这些事件导致同期沉积物中碎屑物质增加.铁桥剖面瓜德鲁普—乐平统界线附近磁化率和蓬莱滩剖面(乐平统底界GSSP)表现出一致的变化趋势,具可对比性.利用多窗谱法(MTM)和傅里叶变换(FT)从磁化率序列中识别出五个米兰科维奇周期:长偏心率周期(E2,405ka)、短偏心率周期(E1,100ka)、长地轴斜率周期(O2,44.1ka)、长岁差周期(P2,20.95ka)和短岁差周期(P1,17.7ka).对比基于E2周期建立的磁性地层磁化率(MSS)带和标准参考带(SRZ),建立整个沉积序列的高分辨率(200ka)FPTS,提出卡匹敦阶的时限为3.85 Ma(存在+0~0.28 Ma误差),整段沉积序列的平均沉积速率为2.91cm·ka^-1.同时计算出卡匹敦阶内部七个牙形石带的时限,从最短26.6ka到最长2.3 Ma.另外,估算出峨眉山大火成岩省喷发启动时间为262.67 Ma,位于瓜德鲁普—乐平统界线之下1.42 Ma.
3719-3734