基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

张超, 姚华建, 童平, 刘沁雅, 雷霆. 2020. 基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究. 地球物理学报, 63(11): 4065-4079, doi: 10.6038/cjg2020O0181
引用本文: 张超, 姚华建, 童平, 刘沁雅, 雷霆. 2020. 基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究. 地球物理学报, 63(11): 4065-4079, doi: 10.6038/cjg2020O0181
ZHANG Chao, YAO HuaJian, TONG Ping, LIU QinYa, LEI Ting. 2020. Joint inversion of linear array ambient noise surface-wave and teleseismic body-wave data based on an adjoint-state method. Chinese Journal of Geophysics (in Chinese), 63(11): 4065-4079, doi: 10.6038/cjg2020O0181
Citation: ZHANG Chao, YAO HuaJian, TONG Ping, LIU QinYa, LEI Ting. 2020. Joint inversion of linear array ambient noise surface-wave and teleseismic body-wave data based on an adjoint-state method. Chinese Journal of Geophysics (in Chinese), 63(11): 4065-4079, doi: 10.6038/cjg2020O0181

基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

  • 基金项目:

    江苏省自然科学基金项目(BK20190499),河海大学中央高校基本科研业务费项目(2019B00714),国家自然科学基金项目(42004037),中国自然科学基金委重大项目课题(41790464)联合资助

详细信息
    作者简介:

    张超, 男, 1988年生, 河海大学讲师, 主要从事地震层析成像、结构反演等研究.E-mail:czhang18@hhu.edu.cn

    通讯作者: 姚华建, 男, 1979年生, 教授, 博士生导师, 主要从事地震波和背景噪声成像、岩石圈结构与变形、大地震破裂过程等领域的研究工作.E-mail:hjyao@ustc.edu.cn
  • 中图分类号: P315

Joint inversion of linear array ambient noise surface-wave and teleseismic body-wave data based on an adjoint-state method

More Information
  • 伴随层析成像(Adjoint Tomography)通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.

  • 加载中
  • 图 1 

    基于伴随方法的背景噪声面波经验格林函数和远震体波联合成像流程图

    Figure 1. 

    The workflow for joint inversion of ambient noise surface-wave empirical Green′s functions and teleseismic body-wave data based on the adjoint-state method

    图 2 

    三种类型伴随层析成像第1次迭代的S波速度敏感核函数及反演模型对比. (a)真实模型:地壳低速体模型,黑色虚线表示Moho面,上下层分别为地壳和地幔,模型参数参见表 1;(b)-(d)分别表示三种类型反演第1次迭代对应的S波速度敏感核函数;(e)-(g)分别表示三种类型反演第1次迭代产生的S波速度模型

    Figure 2. 

    The S-wave velocity sensitivity kernel and inverted model of the 1st iteration for 3 types of adjoint inversion. (a) True model: crust-over-mantle model with low-velocity anomaly in the crust; the black dashed line gives the Moho interface; the model parameters are given in Table 1; (b)-(d) The S wave velocity sensitivity kernel of the 1st iteration for 3 types of inversion; (e)-(g) The inverted S wave velocity model of the 1st iteration for 3 types of inversion.

    图 3 

    三种类型伴随层析成像最终反演模型对比. (a)-(c)不含噪声合成数据联合反演第1,4,8次迭代更新获得的S波速度模型;(d)不含噪声的体波反演的最终速度模型;(e)不含噪声的面波反演的最终S波速度模型; (f)含噪声合成数据联合反演成像结果:其中面波、体波噪声水平分别为2%和5%; (g) (a)-(f)速度模型中X=200 km处的S波速度垂直剖面

    Figure 3. 

    The final S wave velocity models for 3 types of inversion. (a)-(c) The inverted model after 1、4、5 and 8 iterations for joint inversion using synthetic data without noise; (d) The final inverted model for body wave tomography without noise; (e) The final inverted model for surface wave tomography without noise; (f) The inverted model for joint inversion with 2% random noise in surface waves and 5% random noise in body waves; (g) The vertical profile of S wave velocity located at X=200 km in the inverted model in (a)-(f)

    图 4 

    联合反演(不含噪声数据)第1,4,8次迭代波形拟合及数据误差下降曲线. (a)(c)(e)面波波形拟合; (b)(d)(f)体波波形拟合, 蓝、红线分别表示观测和理论波形;(g)数据误差下降曲线, 黑色、蓝色和红色五角星分别代表体波、面波和联合反演

    Figure 4. 

    Waveform fitting and reduction of data misfit for joint inversion without noise in the data. (a)(c)(e) Waveform fitting for surface waves; (b)(d)(f) Waveform fitting for teleseismic body waves; The blue and red waveforms are the observed and synthetic waveforms, respectively; (g) Reduction of data misfit values as a function of iteration number; The black, blue, and red stars represent results for body wave, surface wave, and joint inversion, respectively

    图 5 

    含起伏Moho界面的地壳-地幔模型联合反演. (a)真实模型,黑色实线表示Moho面;(b)-(d)分别为第1,4,7次迭代获得的S波速度模型; (e)-(g)分别为第9,12,15次迭代获得的S波速度模型;白色虚线表示更新的Moho界面

    Figure 5. 

    Joint inversion for the crust-mantle model with an undulating Moho interface. (a) True model, with the black solid line showing the location of Moho interface; (b)-(d) The inverted S wave velocity model after 1、4 and 7 iterations; (e)-(g) The inverted S wave velocity model after 9、12 and 15 iterations; The white dashed line denotes the updating Moho interface

    图 6 

    联合反演第一阶段和第二阶段反演波形拟合.蓝色实线表示观测数据,红色实线表示理论正演数据

    Figure 6. 

    Waveform fitting of the 1st and 2nd step for joint inversion. The blue solid line denotes the observed data, while the red solid line denotes the synthetic data

    图 7 

    联合反演数据误差下降曲线.蓝色和红色五角星分别代表第一阶段和第二阶段反演

    Figure 7. 

    Reduction of data misfit values as a function of iteration number for the 1st step (blue star) and the 2nd step (red star)

    图 8 

    华北克拉通地壳-地幔模型联合反演. (a)真实模型,中下地壳中存在一个低速体,相对速度扰动6%,低速体下方在X=200 km和X=350 km之间出现Moho面突降;为了分析联合反演方法对初始模型的依赖性程度,测试了两种不同的初始模型; (b)含起伏Moho面的地壳-地幔双层模型(简称model-1);(e)含水平Moho面的地壳-地幔双层模型(简称model-2); (c) (f)分别为基于model-1和model-2的联合反演结果;(d) (g)分别为基于model-1和model-2的面波单独反演结果

    Figure 8. 

    Joint inversion for the typical crust-mantle model of North China Craton. (a) True model, where a low velocity anomaly locates in the mid-lower crust with 6% velocity perturbation and there is a Moho depression between X=200 km and X=350km; In order to test the impact of initial model selection on the final model, we set up two different initial models: (b) the crust-mantle model with undulating Moho (model-1); (e) the crust-mantle model with flatting Moho (model-2); (c) (f) The final model for joint inversion based on model-1 and model-2; (d) (g) The final model for surface wave tomography based on model-1 and model-2

    图 9 

    联合反演第一阶段和第二阶段反演波形拟合.蓝色实线表示观测数据,红色实线表示理论正演数据

    Figure 9. 

    Waveform fitting of the 1st and 2nd step for joint inversion. The blue solid line denotes the observed data, while the red solid line denotes the synthetic data

    图 10 

    数据误差下降曲线.蓝色和红色五角星分别代表第一阶段和第二阶段反演

    Figure 10. 

    Reduction of data misfit values as a function of iteration number for the 1st step (blue star) and the 2nd step (red star)

    表 1 

    地壳-地幔双层模型参数

    Table 1. 

    Parameters in the two layer crust-mantle model

    密度
    /(kg·m-3)
    纵波
    /(m·s-1)
    横波
    /(m·s-1)
    厚度/km
    地壳(crust) 2600 5800 3198 60
    地幔(mantle) 3380 8080 4485 40
    下载: 导出CSV

    表 2 

    体波和面波反演中震源参数设置

    Table 2. 

    Source parameter settings for body wave and surface wave tomography

    体波反演 面波反演
    远震事件 入射角 子波类型 虚拟震源 震源位置(x, z) 子波类型
    Event 1
    Event 2
    Event 3
    Event 4
    Event 5
    Event 6
    Event 7
    Event 8
    16°
    20°
    24°
    28°
    332°
    336°
    340°
    344°

    Ricker
    子波
    (主频:
    0.4 Hz)
    Source 1
    Source 2
    Source 3
    Source 4
    Source 5
    Source 6
    Source 7
    Source 8
    (100 km, 0 km)
    (125 km, 0 km)
    (150 km, 0 km)
    (175 km, 0 km)
    (225 km, 0 km)
    (250 km, 0 km)
    (275 km, 0 km)
    (300 km, 0 km)

    高斯
    子波
    (主频:
    0.2 Hz)
    下载: 导出CSV

    表 3 

    数值模拟参数及反演计算时间统计

    Table 3. 

    Parameters of numerical modeling and computational time for inversion

    反演类型 联合反演
    (model-1)
    面波反演
    (model-1)
    联合反演
    (model-2)
    面波反演
    (model-2)
    模型大小
    (NX×NZ)
    X方向:500 km; Z方向: 100 km;
    NX=200; NZ=40
    正演参数 步数:12000, 步长: 0.02 s
    震源个数 32 16 32 16
    CPU耗费 128 64 128 64
    计算时间/
    迭代次数
    4.92 h
    /18次
    3.65 h
    /14次
    5.45 h
    /18次
    3.80 h
    /14次
    反演模型
    RMS误差
    5.41% 12.63% 6.68% 15.64%
    下载: 导出CSV
  •  

    Ammon C J, Randall G E, Zandt G. 1990. On the nonuniqueness of receiver function inversions. Journal of Geophysical Research:Solid Earth, 95(B10):15303-15318. doi: 10.1029/JB095iB10p15303

     

    Bao X W, Sun X X, Xu M J, et al. 2015. Two crustal low-velocity channels beneath SE Tibet revealed by joint inversion of Rayleigh wave dispersion and receiver functions. Earth and Planetary Science Letters, 415:16-24. doi: 10.1016/j.epsl.2015.01.020

     

    Beller S, Monteiller V, Operto S, et al. 2017. Lithospheric architecture of the South-Western Alps revealed by multiparameter teleseismic full-waveform inversion. Geophysical Journal International, 212(2):1369-1388. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1050446023331f84902384b4892bfb37

     

    Bozdaǧ E, Peter D, Lefebvre M, et al. 2016. Global adjoint tomography:first-generation model. Geophysical Journal International, 207(3):1739-1766. doi: 10.1093/gji/ggw356

     

    Brocher T M. 2005. Empirical relations between elastic wavespeeds and density in the Earth's crust. Bulletin of the Seismological Society of America, 95(6):2081-2092. doi: 10.1785/0120050077

     

    Chen L, Wei Z G, Cheng C. 2010. Significant structural variations in the Central and Western North China craton and its implications for the craton destruction. Earth Science Frontiers (in Chinese), 17(1):212-228. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dxqy201001018

     

    Chen M, Huang H, Yao H, et al. 2014. Low wave speed zones in the crust beneath SE Tibet revealed by ambient noise adjoint tomography. Geophysical Research Letters, 41:334-340. doi: 10.1002/2013GL058476

     

    Chen M, Niu F L, Tromp J, et al. 2017. Lithospheric foundering and underthrusting imaged beneath Tibet. Nature Communications, 8:15659. doi: 10.1038/ncomms15659

     

    Chen P, Zhao L, Jordan T H. 2007. Full 3D tomography for the crustal structure of the Los Angeles region. Bulletin of the Seismological Society of America, 97(4):1094-1120. doi: 10.1785/0120060222

     

    Chen Y, Zhou H W, Ge H K. 2005. Seismic array in North China. Journal of Geodesy and Geodynamics (in Chinese), 25(4):1-5. http://en.cnki.com.cn/Article_en/CJFDTOTAL-DKXB200504000.htm

     

    Chen Y L, Niu F L. 2016. Joint inversion of receiver functions and surface waves with enhanced preconditioning on densely distributed CNDSN stations:Crustal and upper mantle structure beneath China. Journal of Geophysical Research:Solid Earth, 121(2):743-766. doi: 10.1002/2015JB012450

     

    Fang H J, Zhang H J, Yao H J, et al. 2016. A new algorithm for three-dimensional joint inversion of body wave and surface wave data and its application to the Southern California plate boundary region. Journal of Geophysical Research:Solid Earth, 121(5):3557-3569. doi: 10.1002/2015JB012702

     

    Fichtner A, Van Herwaarden D P, Afanasiev M, et al. 2018. The collaborative seismic earth model:generation 1. Geophysical Research Letters, 45(9):4007-4016. doi: 10.1029/2018GL077338

     

    Galetti E, Curtis A. 2012. Generalised receiver functions and seismic interferometry. Tectonophysics, 532-535:1-26. doi: 10.1016/j.tecto.2011.12.004

     

    Julià J, Ammon C J, Herrmann R B, et al. 2000. Joint inversion of receiver function and surface wave dispersion observations. Geophysical Journal International, 143(1):99-112. doi: 10.1046/j.1365-246x.2000.00217.x

     

    Komatitsch D, Liu Q Y, Tromp J, et al. 2004. Simulations of ground motion in the Los Angeles Basin based upon the Spectral-Element Method. Bulletin of the Seismological Society of America, 94(1):187-206. doi: 10.1785/0120030077

     

    Li J T, Song X D, Zhu L P, et al. 2017a. Joint inversion of surface wave dispersions and receiver functions with P velocity constraints:Application to Southeastern Tibet. Journal Geophysical Research:Solid Earth, 122(9):7291-7310. doi: 10.1002/2017JB014135

     

    Li Y H, Wang X C, Zhang R Q, et al. 2017b. Crustal structure across the NE Tibetan Plateau and Ordos Block from the joint inversion of receiver functions and Rayleigh-wave dispersions. Tectonophysics, 705:33-41. doi: 10.1016/j.tecto.2017.03.020

     

    Lin F C, Moschetti M P, Ritzwoller M H. 2008. Surface wave tomography of the western United States from ambient seismic noise:Rayleigh and Love wave phase velocity maps. Geophysical Journal International, 173(1):281-298. doi: 10.1111/j.1365-246X.2008.03720.x

     

    Liu C L, Chen H P, Xie J. 2018. Progress in the studies of the joint inversion of surface wave dispersion and receiver functions. Progress in Geophysics (in Chinese), 33(2):479-488, doi:10.6038/pg2018BB0189.

     

    Liu Q, Gu Y J. 2012. Seismic imaging:From classical to adjoint tomography. Tectonophysics, 566-567:31-66. doi: 10.1016/j.tecto.2012.07.006

     

    Liu Q Y, Li Y, Chen J H, et al. 2010. Joint inversion of receiver function and ambient noise based on Bayesian theory. Chinese Journal of Geophysics (in Chinese), 53(11):2603-2612. http://onlinelibrary.wiley.com/doi/10.1002/cjg2.1566/full

     

    Liu S L, Yang D H, Dong X P, et al. 2017. Element-by-element parallel spectral-element methods for 3-D teleseismic wave modeling. Solid Earth, 8(5):969-986. doi: 10.5194/se-8-969-2017

     

    Monteiller V, Chevrot S, Komatitsch D, et al. 2013. A hybrid method to compute short-period synthetic seismograms of teleseismic body waves in a 3-D regional model. Geophysical Journal International, 192(1):230-247. doi: 10.1093/gji/ggs006

     

    Nocedal J. 1980. Updating quasi-Newton matrices with limited storage. Mathematics of Computation, 35(151):773-782. doi: 10.1090/S0025-5718-1980-0572855-7

     

    Owens T J, Zandt G, Taylor S R. 1984. Seismic evidence for an ancient rift beneath the Cumberland Plateau, Tennessee:A detailed analysis of broadband teleseismic P waveforms. Journal of Geophysical Research:Solid Earth, 89(B9):7783-7795. doi: 10.1029/JB089iB09p07783

     

    Shapiro N M, Campillo M. 2004. Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise. Geophysical Research Letters, 31(7):L07614, doi:10.1029/2004GL019491.

     

    Shapiro N M, Campillo M, Stehly L, et al. 2005. High-resolution surface-wave tomography from ambient seismic noise. Science, 307(5715):1615-1618. doi: 10.1126/science.1108339

     

    Shen W S, Ritzwoller M H, Schulte-Pelkum V. 2013. A 3-D model of the crust and uppermost mantle beneath the Central and Western US by joint inversion of receiver functions and surface wave dispersion. Journal of Geophysical Research:Solid Earth, 118(1):262-276. doi: 10.1029/2012JB009602

     

    Tao K, Grand S P, Niu F L. 2018. Seismic structure of the upper mantle beneath eastern asia from full waveform seismic tomography. Geochemistry, Geophysics, Geosystems, 19(8):2732-2763. doi: 10.1029/2018GC007460

     

    Tape C, Liu Q Y, Maggi A, et al. 2009. Adjoint tomography of the southern California crust. Science, 325(5943):988-992. doi: 10.1126/science.1175298

     

    Tarantola A, Valette B. 1982. Generalized nonlinear inverse problems solved using the least squares criterion. Reviews of Geophysics, 20(2):219-232. doi: 10.1029/RG020i002p00219

     

    Tong P, Chen C W, Komatitsch D, et al. 2014. High-resolution seismic array imaging based on an SEM-FK hybrid method. Geophysical Journal International, 197, 369-395. doi: 10.1093/gji/ggt508

     

    Tromp J, Tape C, Liu Q Y. 2005. Seismic tomography, adjoint methods, time reversal and banana-doughnut kernels. Geophysical Journal International, 160(1):195-216. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f03b16f523e68d6e1164cec62546b458

     

    Virieux J, Operto S. 2009. An overview of full-waveform inversion in exploration geophysics. Geophysics, 74(6):WCC1-WCC26. doi: 10.1190/1.3238367

     

    Wang K, Yang Y J, Basini P, et al. 2018. Refined crustal and uppermost mantle structure of southern California by ambient noise adjoint tomography. Geophysical Journal International, 215(2):844-863. doi: 10.1093/gji/ggy312

     

    Wang Y, Chevrot S, Monteiller V, et al. 2016. The deep roots of the western pyrenees revealed by full waveform inversion of teleseismic P waves. Geology, 44(6):475-478. doi: 10.1130/G37812.1

     

    Wu R S, Luo J R, Wu B Y. 2014. Seismic envelope inversion and modulation signal model. Geophysics, 79(3):WA13-WA24. doi: 10.1190/geo2013-0294.1

     

    Yang Y J, Ritzwoller M H, Lin F C, et al. 2008. Structure of the crust and uppermost mantle beneath the western United States revealed by ambient noise and earthquake tomography. Journal of Geophysical Research:Solid Earth, 113(B12):B12310, doi:10.1029/2008JB005833.

     

    Yao H J, Van Der Hilst R D, De Hoop M V. 2006. Surface-wave array tomography in SE Tibet from ambient seismic noise and two-station analysis-I. Phase velocity maps. Geophysical Journal International, 166(2):732-744. doi: 10.1111/j.1365-246X.2006.03028.x

     

    Yao H J, Beghein C, Van Der Hilst R D. 2008. Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis-Ⅱ. Crustal and upper-mantle structure. Geophysical Journal International, 173(1):205-219. doi: 10.1111/j.1365-246X.2007.03696.x

     

    Zhang C, Yao H J, Liu Q Y, et al. 2018. Linear array ambient noise Adjoint tomography reveals intense crust-mantle interactions in North China Craton. Journal of Geophysical Research:Solid Earth, 123(1):368-383. doi: 10.1002/2017JB015019

     

    Zhang C. 2018. Linear array ambient noise and body wave waveform adjoint tomography: methodology and applications[Ph. D. thesis] (in Chinese). Hefei: University of Science and Technology of China.

     

    Zhao L, Jordan T H, Olsen K B, et al. 2005. Fréchet kernels for imaging regional earth structure based on three-dimensional reference models. Bulletin of the Seismological Society of America, 95(6):2066-2080. doi: 10.1785/0120050081

     

    Zheng C, Ding Z F, Song X D. 2018. Joint inversion of surface wave dispersion and receiver functions for crustal and uppermost mantle structure beneath the northern north-south seismic zone. Chinese Journal of Geophysics (in Chinese), 61(4):1211-1224, doi:10.6038/cjg2018L0443.

     

    Zheng T Y, Duan Y H, Xu W W, et al. 2017. A seismic model for crustal structure in North China Craton. Earth and Planetary Physics, 1(1):26-34. doi: 10.26464/epp2017004

     

    Zhu H J, Komatitsch D, Tromp J. 2017. Radial anisotropy of the North American upper mantle based on adjoint tomography with USArray. Geophysical Journal International, 211(1):349-377. doi: 10.1093/gji/ggx305

     

    Zhu H J. 2018. Crustal wave speed structure of North Texas and Oklahoma based on ambient noise cross-correlation functions and adjoint tomography. Geophysical Journal International, 214(1):716-730. doi: 10.1093/gji/ggy169

     

    Zhu L P, Kanamori H. 2000. Moho depth variation in southern California from teleseismic receiver functions. Journal of Geophysical Research:Solid Earth, 105(B2):2969-2980. doi: 10.1029/1999JB900322

     

    Zhu R X, Chen L, Wu F Y, et al. 2011. Timing, scale and mechanism of the destruction of the North China Craton. Science China Earth Sciences, 54(6):789-797. doi: 10.1007/s11430-011-4203-4

     

    Zhu R X, Xu Y G, Zhu G, et al. 2012. Destruction of the North China craton. Science China Earth Sciences, 55(10):1565-1587. doi: 10.1007/s11430-012-4516-y

     

    陈凌, 危自根, 程骋. 2010.从华北克拉通中、西部结构的区域差异性探讨克拉通破坏.地学前缘, 17(1):212-228. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dxqy201001018

     

    陈顒, 周华伟, 葛洪魁. 2005.华北地震台阵探测计划.大地测量与地球动力学, 25(4):1-5. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz200504001

     

    刘成林, 陈浩朋, 谢军. 2018.面波频散与体波接收函数联合反演研究回顾及展望.地球物理学进展, 33(2):479-488, doi:10.6038/pg2018BB0189.

     

    刘启元, 李昱, 陈九辉等. 2010.基于贝叶斯理论的接收函数与环境噪声联合反演.地球物理学报, 53(11):2603-2612, doi:10.3969/j.issn.0001-5733.2010.11.008. http://www.geophy.cn//CN/abstract/abstract3413.shtml

     

    张超. 2018.基于线性台阵背景噪声与体波波形伴随成像方法研究及应用[博士论文].合肥: 中国科学技术大学.

     

    郑晨, 丁志峰, 宋晓东. 2018.面波频散与接收函数联合反演南北地震带北段壳幔速度结构.地球物理学报, 61(4):1211-1224, doi:10.6038/cjg2018L0443. http://www.geophy.cn//CN/abstract/abstract14445.shtml

     

    朱日祥, 陈凌, 吴福元等. 2011.华北克拉通破坏的时间、范围与机制.中国科学:地球科学, 41(5):583-592. http://www.cnki.com.cn/Article/CJFDTotal-JDXK201105001.htm

     

    朱日祥, 徐义刚, 朱光等. 2012.华北克拉通破坏.中国科学:地球科学, 42(8):1135-1159. http://www.cnki.com.cn/Article/CJFDTotal-JDXK201208002.htm

  • 加载中

(10)

(3)

计量
  • 文章访问数: 
  • PDF下载数: 
  • 施引文献:  0
出版历程
收稿日期:  2020-05-14
修回日期:  2020-08-05
上线日期:  2020-11-05

目录