格陵兰岛地区居里面深度.pdf
《格陵兰岛地区居里面深度.pdf》由会员分享,可在线阅读,更多相关《格陵兰岛地区居里面深度.pdf(11页珍藏版)》请在咨信网上搜索。
1、褚伟,徐亚,张健,张倩文,路书鹏,胡琪鑫.2023.格陵兰岛地区居里面深度.地震学报,45(3):471481.doi:10.11939/jass.20220167.ChuW,XuY,ZhangJ,ZhangQW,LuSP,HuQX.2023.CuriedepthinGreenlandarea.Acta Seismologica Sinica,45(3):471481.doi:10.11939/jass.20220167.格陵兰岛地区居里面深度*褚伟1)徐亚1,2),张健1)张倩文1,2)路书鹏1,2)胡琪鑫1,2)1)中国北京100049中国科学院大学地球与行星科学学院2)中国北京10002
2、9中国科学院地质与地球物理研究所hbhbhmhcQshmQshc摘要基于全球磁异常数据,利用质心法计算格陵兰岛地区的居里面深度,并结合最新的热流资料、地壳结构和大洋年龄进行分析,加深了对研究区热结构的认识。首先,总结了磁性体满足不同分布时,居里面深度计算的理论公式和校正因子的选择。在格陵兰岛地区应用分形分布,结果表明,研究区居里面深度在 40km 以内,平均值约为 20km,估算误差约为 2.7km。校正得到地形到居里面的距离和沉积层底界面到居里面的距离,并根据不同的地壳类型对地表实测热流数据校正得到。对比结果显示与基本成负相关,计算结果合理。进一步分析表明,研究区居里面总体上浅于莫霍面,但在
3、格陵兰岛北部部分区域居里面深于莫霍面,在大洋地区的随洋壳年龄的增大而增大。关键词 格陵兰岛居里面深度分形分布热流doi:10.11939/jass.20220167中图分类号:P314.2文献标识码:ACurie depth in Greenland areaChuWei1)XuYa1,2),ZhangJian1)ZhangQianwen1,2)LuShupeng1,2)HuQixin1,2)1)College of Earth and Planetary Sciences,University of Chinese Academy of Sciences,Beijing 100049,Chi
4、na2)Institute of Geology and Geophysics,Chinese Academy of Sciences,Beijing 100029,ChinahbhbhmhcAbstract:ThispaperestimatestheCuriedepthinGreenlandareabasedontheglobalmagneticanomalydatabyusingthecentroidmethod.Forfurtheranalysisweintroducethelatestheatflowdata,crustalstructureandoceanage,whichhelps
5、todeepentheunderstandingofthermalstructureinthestudiedarea.Whenthemagneticbodysatisfiesdifferentdistributionassump-tions,thedifferenttheoreticalformulasandcorrectionfactorsforcalculatingtheCuriedeptharesummarized.Amongthesedistributionassumptions,weapplythefractaldistributionassump-tiontotheGreenlan
6、darea.TheresultshowsthattheCuriedepth()inthestudiedareaiswithin40km,withanaveragedepthofabout20kmandanestimationerrorofabout2.7km.Basedon,wecalculatethedistancefromtopographytoCuriedepthandthedistance*基金项目国家基金委项目(42074092)和中国科学院青年创新促进会(2016064)联合资助.收稿日期2022-09-09 收到初稿,2022-10-26 决定采用修改稿.作者简介褚伟,在读博士研
7、究生,主要从事海洋地球物理方面的研究,e-mail:;徐亚,博士,副研究员,主要从事海洋地球物理、综合地球物理研究,e-mail:第45卷第3期地震学报Vol.45,No.32023年5月(471481)ACTASEISMOLOGICASINICAMay,2023QsQshmhbhbhcfromthetopcrystallinebasementtoCuriedepth.Besides,wecorrectsurfaceheatflowtoob-tainaccordingtodifferentcrustaltypes.Wefindthatisapproximatelyinverselypropor
8、-tionalto,whichprovesthattheresultisreasonable.FurtheranalysisshowsthatinthestudiedareaisgenerallyshallowerthanMohodepth,whileinpartofthenorthernGreenland,isdeeperthanMoho.Intheocean,increaseswiththeageoftheoceancrust.Key words:Greenland;Curiedepth;thefractaldistribution;geothermalheatflow 引言格陵兰岛位于北
9、美地区,北侧为北冰洋,西侧为巴芬湾、戴维斯海峡,西南方向与加拿大隔海相望,南侧为拉布拉多海,东侧为格陵兰海、丹麦海峡,东南方向与冰岛和欧洲大陆隔海相望。格陵兰岛 80%的面积被冰雪覆盖,形成的格陵兰冰盖是地球上第二大淡水水库。格陵兰岛附近的海平面变化与沿海经济规划甚至人类生存息息相关,对于全球气候变化具有重要的指示意义(江伟伟等,2011)。由于岛上冰雪覆盖,冰盖之下的地形地貌和地质情况相对难以获得。前人对岛上地形地貌的研究表明,格陵兰岛冰盖之下存在一条从中部延伸至西北部的长 750km 的大峡谷(Bamberet al,2013),岛西北部的冰层之下埋藏着一个直径约为31km 的撞击坑(Ma
10、cGregoret al,2019)。冰雪覆盖进一步增加了地质研究的难度,该地区只有海岸线附近才能看到岩石露头,因此重、磁、热、地震等地球物理方法是揭示格陵兰岛地区地壳结构、岩石圈结构和热结构的有效手段。前人对格陵兰岛地区的地壳和岩石圈结构作了大量研究。其中对莫霍面深度的研究表明:莫霍面在格陵兰岛内陆较深,约为 3550km;在大陆边缘较浅,约为 2535km;南部莫霍面较深,可达到 55km;北部莫霍面较浅,只有 35km;平均莫霍面深度约为 40km(Steffenet al,2017;Toyokuniet al,2020)。对于地壳结构,重力和地震给出的平均地壳厚度分别约为 40km 和
11、 36km,其中格陵兰岛东部的地壳最厚,中部的 NW-SE 一线上也存在厚地壳(Martoset al,2018)。对于岩石圈结构,格陵兰岛西南部岩石圈厚约 200270km,向中部方向薄至 180190km,NW-SE 一带岩石圈异常薄,约 100120km,可能标志着冰岛热点轨迹(Artemieva,2019)。研究区的热结构,特别是格陵兰岛的热点迁移问题引起了众多学者的广泛讨论。有学者认为冰岛地幔柱在约 8020Ma 经过格陵兰板块,导致了该地区地热通量的变化(Rogozhinaet al,2016),许多研究结果都验证了这一假设。例如,研究区存在岩石圈变薄和高密度镁铁物质的显著底侵,地
12、震层析成像研究显示岛中部莫霍面上方存在一个细长的剪切波高速区,可能对应高密度区,地震波低速带则反映了冰岛地幔柱的余热(Martoset al,2018;Toyokuniet al,2020a,b)。居里面深度表示磁性层底部的深度,大致相当于 580 的等温线所在深度(Rosset al,2006),因此对于揭示区域地热构造具有重要意义。Martos 等(2018)估计格陵兰岛居里面深度介于 2348km 之间,平均值为 37km;格陵兰岛西部、北部以及埃尔斯米尔岛的内部居里面最深,格陵兰岛东部向冰岛热点的方向居里面深度较浅。Lu 等(2022)认为格陵兰岛北部的居里面深度大于 30km,南部约
13、为 2030km,巴芬湾约为 1837km,扬马延岛微大陆的居里面比其周围洋盆更深。因此,对格陵兰岛地区的居里面结构及其热状态的认识有待进一步深化。本文从全球磁异常数据出发,拟利用质心法计算格陵兰岛地区的居里面深度,并与最新的热流数据结合,472地震学报45卷以期给出整个格陵兰岛地区的热结构。1 地质和地球物理背景格陵兰岛是一块存在了约 40 亿年的稳定陆地,主要由前寒武纪省组成,在早元古代发生碰撞,之后大陆几乎再未发生构造改造。格陵兰岛东侧为大西洋中脊(图 1),自南向北分布着冰岛热点、扬马延岛热点和斯瓦尔巴群岛高地热地区,因此格陵兰岛地区的热结构引发了国内外学者的广泛讨论。居里面深度表示磁
14、性层底部的深度,这一深度大致相当于 580 的等温线,该温度也称为居里面温度,高于该温度,磁性材料将失去铁磁性,因此利用磁异常估算居里面深度能够独立揭示区域地热构造。本研究使用的磁异常数据 来 源 于 EMAG2(EarthMagneticAnomalyGrid2-arc-minuteresolution)(Mauset al,2009),磁异常分布如图 2 所示,数据分辨率为 2 弧分,观测面位于大地水准面之上 4km。此外,可靠的热流测量数据是约束岩石圈热结构的重要依据,而格陵兰岛大陆被冰层覆盖,因此几乎没有热流观测数据,现有的测量数据主要集中在活动的海洋盆地。面对这一问题,不同学者给出了
15、利用已有数据求解热流值的方法。Martos 等(2018)基于全球磁异常数据,利用谱分析和求解热传导方程的方法计算格陵兰岛热流值及其不确定性,给 出 的 热 流 值 变 化 范 围 为 5075mW/m2,平均值为 60mW/m2,东格陵兰岛显示高热流值,西格陵兰岛显示低热流值。Artemieva(2019)结合地震、地形和冰厚数据计算得到的格陵兰岛西南部热流值约为 40mW/m2,NW-SE 带的热流值约为 6070mW/m2,东格陵兰岛热流值大于 100mW/m2。Colgan 等(2022)整理了格陵兰岛所有可获得的热流测量数据共419 个(图 3a),其中海上热流数据 317 个,陆上
16、热流数据 102 个,按照数据来源,其中 290 个为 2018 年国际热流协会已发表的数据,129 个为新整理的数据。基于上述热流测量数据,并543210123海拔/km格陵兰岛冰岛埃尔斯米尔岛斯瓦尔巴群岛扬马延岛大 西 洋 北极圈巴芬湾拉布拉多海90 W6030075N 60中 脊戴维斯海峡丹麦海峡格陵兰海图1研究区地形图Fig.1Thetopographymapofthestudiedarea磁异常/nT90 W6030075N 6030025020015010050050100150200250图2研究区磁异常图Fig.2Themapofmagneticanomalyinthestud
17、iedarea3期褚伟等:格陵兰岛地区居里面深度473结合其它地球物理数据集,Colgan 等(2022)使用机器学习方法得到了格陵兰岛及其周围区域的热流图(图 3b)。本文使用这一最新热流数据集结合居里面深度计算结果分析格陵兰岛地区的热结构。2 研究方法TAT计算居里面深度时,通常需要假设磁性体场源满足特定的磁化分布,从而使用相应的理论计算公式,计算方法包括质心法、谱峰法、非线性参数估计法等(Kelemework et al,2021)。本文使用最常用的质心法,假设磁性体场源满足分形分布,则磁异常总场的径向平均功率谱的对数可以写为lnATC1|2k|htln|2k|12ln1exp|2k|(
18、hbht),(1)C1kht和hb式中,为常数,为波数,分别为磁性体的顶、底界面深度,为分形校正因子。对于各向同性的 3-D 分形磁化强度分布,在波长小于大约两倍的磁性层厚度(中到高波数)时,式(1)简化为ln(|2k|12AT)C2|2k|ht,(2)在低波数时,式(1)则简化为ln(|2k|32AT)C3|2k|h0,(3)C2和C3h0式中,为常数,为磁性层的质心深度。Thth0hb2h0ht根据上述径向平均功率谱的对数表达式可以看出,磁性体的顶界面深度和质心深度均可以表示为磁异常总场的径向平均功率谱对数与波数的线性回归,因此由磁异常功率谱斜率可分别求出场源的顶界面深度和质心深度,再利用
19、求出场源的底界面深50100150热流值/(mWm-1)50100150热流值/(mWm-1)90 W6030090 W6030075N 6075N 60(a)(b)图3研究区热流实测数据(a)和机器学习预测数据(b)(数据来源于 Colganet al,2022)Fig.3Geothermalheatfluxinthestudiedareabymeasuring(a)andpredictingusingmachinelearningmethod(b)(datafromColganet al,2022)474地震学报45卷hb2h0hth0hth0和ht度。同样地,的误差使用计算得到,其中和分
20、别为的数据拟合误差。AT=1=3.9 3,5除了假设磁性体场源满足分形分布之外,还可以假设场源满足随机不相关分布或者统计分布(Kelemeworket al,2021)。经过整理发现,不同的磁化分布假设实际上具有相同形式的理论计算公式,只是需要对作不同形式的校正,即选择不同的:随机不相关分布假设,统计分布假设,分形分布假设。其中,在分形磁化强度分布假设下,可以选择常数,也可以基于数值迭代求解,但是变化的校正因子通常导致过校正(Liet al,2019),常数校正因子虽然会给居里面深度的估计带来额外的误差,但是在相似的地质单元中,一般不会有太大的变化,因此本文对取常数。具体计算时,使用滑动窗口法
21、,根据每一个窗口内的磁异常径向功率谱计算居里面深度。窗口大小的选择也很重要,由于理论模型在水平方向无限延伸,因此窗口过小会限制窗口内可以捕获的最大波长,窗口过大则会导致不同构造单元的效应混淆在一起,造成结果分辨率低。一般来说,窗口大小可以取 50200km。此外,为了减小窗口尺寸大小的影响、压制随机噪声同时提高分辨率,可以选择对不同大小窗口的结果取平均(Liet al,2013)。TTwinTwinATATln(|2k|12AT)ln(|2k|32AT)hth0kbinhth0hth0hb本文计算居里面深度的算法流程如下:输入磁异常数据,使用滑动窗口分割数据,每一个滑动窗口内的磁异常数据记为;
22、对作傅里叶变换得到径向平均功率谱,结果记为;对作分形校正并取对数,得到结果和,分别用于求取和;将位于同一波数段(波数段宽度取)的数据归为一组,用这组数据的均值代表该波数段,即可绘制散点图;对高波数段和低波数段的频谱数据作线性拟合,斜率分别为和,数据拟合误差分别为和,可相应计算得到与误差,并将结果作为该窗口中心点的结果;所有滑动窗口内的磁异常数据完成上述 25 步计算。3 计算结果及分析 3.1 计算结果=3kbin0.006 km1截取研究区磁异常数据并投影到笛卡尔坐标系,通过质心法计算每一个滑动窗口内的居里面深度,选择窗口大小 200km,步长 100km,分形校正因子,波数段宽度,磁性体的
23、顶界面和质心深度的拟合波数段分别为 0.0390.081km1和 0.0030.033km1。某一滑动窗口内磁异常径向平均功率谱对数与波数的线性拟合情况如图 4 所示,其中图 4a 和 4b 分别为高波数段和低波数段拟合所得顶界面深度和质心深度的示例。123452460020.20.40.60.81.21.0ln(|2k|AT)lnAT2k/(2km-1)00.20.40.60.81.21.02k/(2km-1)ht=6.5 kmh0=12.6 km(a)(b)hth0图4某一窗口在高波数段估计(a)和低波数段估计(b)的示例hth0Fig.4Exampleofestimating athig
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 格陵兰岛 地区 居里 深度
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。