留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

卫星大地测量成像地震周期形变研究综述

郭汝梦 杨浩哲 汤雄伟 张文婷 徐晓雪 刘德川 孙和平

郭汝梦, 杨浩哲, 汤雄伟, 张文婷, 徐晓雪, 刘德川, 孙和平. 卫星大地测量成像地震周期形变研究综述[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 799-806. doi: 10.13203/j.whugis20220166
引用本文: 郭汝梦, 杨浩哲, 汤雄伟, 张文婷, 徐晓雪, 刘德川, 孙和平. 卫星大地测量成像地震周期形变研究综述[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 799-806. doi: 10.13203/j.whugis20220166
GUO Rumeng, YANG Haozhe, TANG Xiongwei, ZHANG Wenting, XU Xiaoxue, LIU Dechuan, SUN Heping. A Review on Satellite Geodesy Applied to Image the Earthquake Cycle Deformation[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 799-806. doi: 10.13203/j.whugis20220166
Citation: GUO Rumeng, YANG Haozhe, TANG Xiongwei, ZHANG Wenting, XU Xiaoxue, LIU Dechuan, SUN Heping. A Review on Satellite Geodesy Applied to Image the Earthquake Cycle Deformation[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 799-806. doi: 10.13203/j.whugis20220166

卫星大地测量成像地震周期形变研究综述

doi: 10.13203/j.whugis20220166
基金项目: 

中国科学院战略性先导科技专项 XDB41000000

武汉引力与固体潮国家野外科学观测研究站开放研究基金 WHYWZ202111

详细信息
    作者简介:

    郭汝梦,博士,副研究员,研究方向为地震周期断层形变。guorm@apm.ac.cn

    通讯作者: 孙和平,博士,研究员,中国科学院院士。heping@whigg.ac.cn
  • 中图分类号: P315; P22

A Review on Satellite Geodesy Applied to Image the Earthquake Cycle Deformation

Funds: 

The Strategic Priority Research Program of Chinese Academy of Sciences XDB41000000

Open Fund of Wuhan, Gravitation and Solid Earth Tides, National Observation and Research Station WHYWZ202111

More Information
    Author Bio:

    GUO Rumeng, PhD, associate professor, specializes in fault deformation associated with earthquake cycle. E-mail: guorm@apm.ac.cn

    Corresponding author: SUN Heping, PhD, professor, Academician of Chinese Academy of Sciences.E-mail: heping@whigg.ac.cn
  • 摘要: 近年来,卫星大地测量技术的快速发展为精确测量地壳形变和断层行为提供了前所未有的多维观测数据。它与地震学的结合使得地震周期形变监测的时空分辨率大大提升,为更加深入地研究地震周期过程和机理提供了一个窗口,地震大地测量学应运而生。它能够对地壳运动进行定量描述、对断层活动进行精准建模,从而为洞悉整个地震周期过程的应力应变演化提供科学依据,同时为评估地震危险性、实现地震预测预警提供科学指导。以卫星大地测量观测探究断层形变为主线,分析了断层处于地震周期不同阶段的运动学特征(震间、同震和震后),回顾了地震大地测量学在震源物理方面的一些重要发现。研究表明,利用卫星大地测量数据判定断层所处地震周期的阶段是实现地震预测的可行思路。
  • 图  1  地震周期断层形变图(改自文献[9])

    Figure  1.  Map of Earthquake-Cycle Fault Deformation (Modified from Reference [9])

    图  2  走滑断层不同孕震阶段形变模式(改自文献[16])

    Figure  2.  Deformation Pattern of Transform Fault inDifferent Stages(Modified from Reference [16])

    图  3  基于不同数据有限断层反演结果的时间、空间以及地震矩分辨率(改自文献[37])

    Figure  3.  Temporal, Spatial, and Moment Resolutions of Fault Inversion Results by Using Different Datasets (Modified from Reference [37])

    图  4  震后形变概念示意图

    Figure  4.  Schematic Diagram of Postseismic Deformation

  • [1] Segall P, Davis J L. GPS Applications for Geodynamics and Earthquake Studies[J]. Annual Review of Earth and Planetary Sciences, 1997, 25: 301-336 doi:  10.1146/annurev.earth.25.1.301
    [2] Bürgmann R, Rosen P A, Fielding E J. Synthetic Aperture Radar Interferometry to Measure Earth's Surface Topography and Its Deformation[J]. Annual Review of Earth and Planetary Sciences, 2000, 28: 169-209 doi:  10.1146/annurev.earth.28.1.169
    [3] 沈正康. 卫星大地测量用于东亚大陆地球动力学与地震学研究回顾[J]. 地球物理学报, 2021, 64(10): 3514-3520 doi:  10.6038/cjg2021P0541

    Shen Zhengkang. Satellite Geodesy Applied to Geodynamic and Seismological Studies in East Asia: A Review[J]. Chinese Journal of Geophysics, 2021, 64(10): 3514-3520 doi:  10.6038/cjg2021P0541
    [4] Elliott J R, Walters R J, Wright T J. The Role of Space-Based Observation in Understanding and Responding to Active Tectonics and Earthquakes[J]. Nature Communications, 2016, 7: 13844 doi:  10.1038/ncomms13844
    [5] 瞿伟, 高源, 陈海禄, 等. 利用GPS高精度监测数据开展青藏高原现今地壳运动与形变特征研究进展[J]. 地球科学与环境学报, 2021, 43(1): 182-204 https://www.cnki.com.cn/Article/CJFDTOTAL-XAGX202101013.htm

    Qu Wei, Gao Yuan, Chen Hailu, et al. Review on Characteristics of Present Crustal Tectonic Movement and Deformation in Qinghai-Tibet Plateau, China Using GPS High Precision Monitoring Data[J]. Journal of Earth Sciences and Environment, 2021, 43(1): 182-204 https://www.cnki.com.cn/Article/CJFDTOTAL-XAGX202101013.htm
    [6] Zhang W T, Ji L Y, Zhu L Y, et al. Current Slip and Strain Rate Distribution Along the Ganzi-Yushu-Xianshuihe Fault System Based on InSAR and GPS Observations[J]. Frontiers in Earth Science, 2022, 10: 821761 doi:  10.3389/feart.2022.821761
    [7] 季灵运, 朱良玉, 李宁, 等. 基于大地测量观测的断层运动研究综述[J]. 大地测量与地球动力学, 2017, 37(8): 771-776 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201708001.htm

    Ji Lingyun, Zhu Liangyu, Li Ning, et al. Review of Fault Movement Based on Geodetic Observations[J]. Journal of Geodesy and Geodynamics, 2017, 37(8): 771-776 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201708001.htm
    [8] Guo R M, Zheng Y, Xu J Q. Stress Modulation of the Seismic Gap Between the 2008 Ms 8.0 Wenchuan Earthquake and the 2013 Ms 7.0 Lushan Earthquake and Implications for Seismic Hazard[J]. Geophysical Journal International, 2020, 221(3): 2113-2125 doi:  10.1093/gji/ggaa143
    [9] Feigl K L, Thatcher W. Geodetic Observations of Post-Seismic Transients in the Context of the Earthquake Deformation Cycle[J]. Comptes Rendus Geoscience, 2006, 338(14/15): 1012-1028
    [10] Zheng Y, Guo R M. Earthquake Potential of the Seismic Gap Between the Wenchuan and Lushan Earthquakes: Current Status, Thoughts, and Challenges[J]. Science China Earth Sciences, 2021, 64(3): 503-506 doi:  10.1007/s11430-020-9719-y
    [11] Lindsey E O, Mallick R, Hubbard J A, et al. Slip Rate Deficit and Earthquake Potential on Shallow Megathrusts[J]. Nature Geoscience, 2021, 14(5): 321-326 doi:  10.1038/s41561-021-00736-x
    [12] Fialko Y. Interseismic Strain Accumulation and the Earthquake Potential on the Southern San Andreas Fault System[J]. Nature, 2006, 441(7096): 968-971 doi:  10.1038/nature04797
    [13] 邹镇宇, 江在森, 武艳强, 等. 利用带倾角断层形变公式研究川滇块体东边界断裂带形变特征[J]. 武汉大学学报·信息科学版, 2018, 43(11): 1688-1695 doi:  10.13203/j.whugis20160415

    Zou Zhenyu, Jiang Zaisen, Wu Yanqiang, et al. Deformation Characteristics of the Eastern Boundary Fault Zone of Sichuan-Yunnan Block Using the Deformation Formula of Seismic Fault with Dip Angle[J]. Geomatics and Information Science of Wuhan University, 2018, 43(11): 1688-1695 doi:  10.13203/j.whugis20160415
    [14] Meade B J. Block Models of Crustal Motion in Southern California Constrained by GPS Measurements[J]. Journal of Geophysical Research, 2005, 110: B03403
    [15] Wang K L, Zhu Y J, Nissen E, et al. On the Relevance of Geodetic Deformation Rates to Earthquake Potential[J]. Geophysical Research Letters, 2021, 48(11): e2021GL093231
    [16] Diao F Q, Xiong X, Wang R J, et al. Slip Rate Variation Along the Kunlun Fault (Tibet): Results from New GPS Observations and a Viscoelastic Earthquake-Cycle Deformation Model[J]. Geophysical Research Letters, 2019, 46(5): 2524-2533 doi:  10.1029/2019GL081940
    [17] Guo R M, Zheng Y, Tian W, et al. Locking Status and Earthquake Potential Hazard Along the Middle-South Xianshuihe Fault[J]. Remote Sensing, 2018, 10(12): 2048 doi:  10.3390/rs10122048
    [18] Daout S, Jolivet R, Lasserre C, et al. Along-Strike Variations of the Partitioning of Convergence Across the Haiyuan Fault System Detected by InSAR[J]. Geophysical Journal International, 2016, 205(1): 536-547 doi:  10.1093/gji/ggw028
    [19] Jolivet R, Lasserre C, Doin M P, et al. Shallow Creep on the Haiyuan Fault (Gansu, China) Revealed by SAR Interferometry[J]. Journal of Geophysical Research: Solid Earth, 2012, 117(B6): e2000JB008732
    [20] Guo R M, Yang H F, Li Y, et al. Complex Slip Distribution of the 2021 Mw 7.4 Maduo, China, Earthquake: An Event Occurring on the Slowly Slipping Fault[J]. Seismological Research Letters, 2022, 93(2A): 653-665 doi:  10.1785/0220210226
    [21] Meade B J, Klinger Y, Hetland E A. Inference of Multiple Earthquake-Cycle Relaxation Timescales from Irregular Geodetic Sampling of Interseismic Deformation[J]. Bulletin of the Seismological Society of America, 2013, 103(5): 2824-2835 doi:  10.1785/0120130006
    [22] Wang K L, Hu Y, He J H. Deformation Cycles of Subduction Earthquakes in a Viscoelastic Earth[J]. Nature, 2012, 484(7394): 327-332 doi:  10.1038/nature11032
    [23] Hilley G E, Johnson K M, Wang M, et al. Earthquake-Cycle Deformation and Fault Slip Rates in Northern Tibet[J]. Geology, 2009, 37(1): 31-34 doi:  10.1130/G25157A.1
    [24] Chuang R Y, Johnson K M. Reconciling Geologic and Geodetic Model Fault Slip-Rate Discrepancies in Southern California: Consideration of Nonsteady Mantle Flow and Lower Crustal Fault Creep[J]. Geo-logy, 2011, 39(7): 627-630
    [25] Bruhat L, Segall P. Deformation Rates in Northern Cascadia Consistent with Slow Updip Propagation of Deep Interseismic Creep[J]. Geophysical Journal International, 2017, 211(1): 427-449 doi:  10.1093/gji/ggx317
    [26] Lay T, Kanamori H. An Asperity Model of Large Earthquake Sequences[M]//Earthquake Prediction: An International Review. Washington DC: American Geophysical Union, 1981: 579–592
    [27] Guo R M, Zheng Y, An C, et al. The 2018 Mw 7.9 Offshore Kodiak, Alaska, Earthquake: An Unusual Outer Rise Strike-Slip Earthquake[J]. Journal of Geophysical Research: Solid Earth, 2020, 125(5): e2019JB019267
    [28] Yue H, Lay T, Koper K D. En Échelon and Orthogonal Fault Ruptures of the 11 April 2012 Great Intraplate Earthquakes[J]. Nature, 2012, 490(7419): 245-249 doi:  10.1038/nature11492
    [29] An C, Yue H, Sun J B, et al. The 2015 Mw 8.3 Illapel, Chile, Earthquake: Direction-Reversed Along-Dip Rupture with Localized Water Reverberation[J]. Bulletin of the Seismological Society of America, 2017, 107(5): 2416-2426 doi:  10.1785/0120160393
    [30] Yao H J, Gerstoft P, Shearer P M, et al. Compressive Sensing of the Tohoku-Oki Mw 9.0 Earthquake: Frequency-Dependent Rupture Modes[J]. Geophysical Research Letters, 2011, 38.
    [31] Tang X W, Guo R M, Xu J Q, et al. Probing the Fault Complexity of the 2017 Ms 7.0 Jiuzhaigou Earthquake Based on the InSAR Data[J]. Remote Sensing, 2021, 13(8): 1573 doi:  10.3390/rs13081573
    [32] Guo R M, Yang H F, Zhu Y F, et al. Narrow Rupture of the 2020 Mw 7.4 La Crucecita, Mexico, Earthquake[J]. Seismological Research Letters, 2021, 92(3): 1891-1899 doi:  10.1785/0220200328
    [33] 许才军, 邓长勇, 周力璇. 利用方差分量估计的地震同震滑动分布反演[J]. 武汉大学学报·信息科学版, 2016, 41(1): 37-44 doi:  10.13203/j.whugis20150500

    Xu Caijun, Deng Changyong, Zhou Lixuan. Coseismic Slip Distribution Inversion Method Based on the Variance Component Estimation[J]. Geomatics and Information Science of Wuhan University, 2016, 41(1): 37-44 doi:  10.13203/j.whugis20150500
    [34] Ji C. Source Description of the 1999 Hector Mine, California, Earthquake, Part Ⅰ: Wavelet Domain Inversion Theory and Resolution Analysis[J]. Bulletin of the Seismological Society of America, 2002, 92(4): 1192-1207 doi:  10.1785/0120000916
    [35] Hayes G P. The Finite, Kinematic Rupture Properties of Great-Sized Earthquakes Since 1990[J]. Earth and Planetary Science Letters, 2017, 468: 94-100 doi:  10.1016/j.epsl.2017.04.003
    [36] Hartzell S, Mendoza C, Ramirez-Guzman L, et al. Rupture History of the 2008 Mw 7.9 Wenchuan, China, Earthquake: Evaluation of Separate and Joint Inversions of Geodetic, Teleseismic, and Strong-Motion Data[J]. Bulletin of the Seismological Society of America, 2013, 103(1): 353-370 doi:  10.1785/0120120108
    [37] Yue H, Zhang Y, Ge Z X, et al. Resolving Rupture Processes of Great Earthquakes: Reviews and Perspective from Fast Response to Joint Inversion[J]. Science China Earth Sciences, 2020, 63(4): 492-511 doi:  10.1007/s11430-019-9549-1
    [38] Xu C, Liu Y, Wen Y, et al. Coseismic Slip Distribution of the 2008 Mw 7.9 Wenchuan Earthquake from Joint Inversion of GPS and InSAR Data[J]. Bulletin of the Seismological Society of America, 2010, 100(5B): 2736-2749 doi:  10.1785/0120090253
    [39] Shen Z K, Sun J, Zhang P, et al. Slip Maxima at Fault Junctions and Rupturing of Barriers During the 2008 Wenchuan Earthquake[J]. Nature Geoscience, 2009, 2(10): 718-724 doi:  10.1038/ngeo636
    [40] Guo R M, Zheng Y, Diao F Q, et al. Rupture Model of the 2013 Mw 6.6 Lushan (China) Earthquake Constrained by a New GPS Data Set and Its Effects on Potential Seismic Hazard[J]. Earthquake Science, 2018, 31(3): 117-125 doi:  10.29382/eqs-2018-0117-1
    [41] Yang H F, Wang D, Guo R M, et al. Rapid Report of the 8 January 2022 ​Ms 6.9 Menyuan Earthquake, Qinghai, China[J]. Earthquake Research Advances, 2022, 2(1): 100113 doi:  10.1016/j.eqrea.2022.100113
    [42] Bürgmann R, Dresen G. Rheology of the Lower Crust and Upper Mantle: Evidence from Rock Mechanics, Geodesy, and Field Observations[J]. Annual Review of Earth and Planetary Sciences, 2008, 36: 531-567 doi:  10.1146/annurev.earth.36.031207.124326
    [43] Barbot S, Fialko Y. A Unified Continuum Representation of Post-Seismic Relaxation Mechanisms: Semi-Analytic Models of Afterslip, Poroelastic Rebound and Viscoelastic Flow[J]. Geophysical Journal International, 2010, 182(3): 1124-1140 doi:  10.1111/j.1365-246X.2010.04678.x
    [44] 李志才, 许才军, 张鹏, 等. 顾及地壳粘弹性结构的地震断层震后形变反演分析[J]. 武汉大学学报·信息科学版, 2014, 39(12): 1477-1481 http://ch.whu.edu.cn/article/id/3144

    Li Zhicai, Xu Caijun, Zhang Peng, et al. Post-Seismic Deformation Inversion of Seismic Fault Considering the Crustal Viscoelastic Structure[J]. Geomatics and Information Science of Wuhan University, 2014, 39(12): 1477-1481 http://ch.whu.edu.cn/article/id/3144
    [45] Bedford J, Moreno M, Li S Y, et al. Separating Rapid Relocking, Afterslip, and Viscoelastic Relaxation: An Application of the Postseismic Straightening Method to the Maule 2010 CGPS[J]. Journal of Geophysical Research: Solid Earth, 2016, 121(10): 7618-7638 doi:  10.1002/2016JB013093
    [46] Marone C J, Scholtz C H, Bilham R. On the Mechanics of Earthquake Afterslip[J]. Journal of Geophysical Research, 1991, 96(B5): 8441 doi:  10.1029/91JB00275
    [47] Guo R M, Zheng Y, Xu J Q, et al. Seismic and Aseismic Fault Slip Associated with the 2017 Mw 8.2 Chiapas, Mexico, Earthquake Sequence[J]. Seismological Research Letters, 2019, 90(3): 1111-1120 doi:  10.1785/0220180262
    [48] He L J, Feng G C, Wu X X, et al. Coseismic and Early Postseismic Slip Models of the 2021 Mw 7.4 Maduo Earthquake (Western China) Estimated by Space-Based Geodetic Data[J]. Geophysical Research Letters, 2021, 48(24): e2021GL095860
    [49] Barbot S, Fialko Y, Bock Y. Postseismic Deformation Due to the Mw 6.0 2004 Parkfield Earthquake: Stress-Driven Creep on a Fault with Spatially Variable Rate-and-State Friction Parameters[J]. Journal of Geophysical Research, 2009, 114: B07405
    [50] Tian Z, Freymueller J T, Yang Z Q. Postseismic Deformation Due to the 2012 Mw 7.8 Haida Gwaii and 2013 Mw 7.5 Craig Earthquakes and Its Implications for Regional Rheological Structure[J]. Journal of Geophysical Research: Solid Earth, 2021, 126(2): e2020JB020197
    [51] Wang K, Fialko Y. Observations and Modeling of Coseismic and Postseismic Deformation Due to the 2015 Mw 7.8 Gorkha (Nepal) Earthquake[J]. Journal of Geophysical Research: Solid Earth, 2018, 123(1): 761-779 doi:  10.1002/2017JB014620
    [52] Zhao B, Bürgmann R, Wang D Z, et al. Dominant Controls of Downdip Afterslip and Viscous Relaxation on the Postseismic Displacements Following the Mw 7.9 Gorkha, Nepal, Earthquake[J]. Journal of Geophysical Research: Solid Earth, 2017, 122(10): 8376-8401 doi:  10.1002/2017JB014366
    [53] Johnson K M, Fukuda J, Segall P. Challenging the Rate-State Asperity Model: Afterslip Following the 2011 M 9 Tohoku-Oki, Japan, Earthquake[J]. Geophysical Research Letters, 2012, 39(20): 2012GL052901
    [54] Shrivastava M N, González G, Moreno M, et al. Coseismic Slip and Afterslip of the 2015 Mw 8.3 Illapel (Chile) Earthquake Determined from Continuous GPS Data[J]. Geophysical Research Letters, 2016, 43: 710-719
    [55] Sun T, Wang K L, Iinuma T, et al. Prevalence of Viscoelastic Relaxation After the 2011 Tohoku-Oki Earthquake[J]. Nature, 2014, 514(7520): 84-87 doi:  10.1038/nature13778
    [56] Huang M H, Bürgmann R, Freed A M. Probing the Lithospheric Rheology Across the Eastern Margin of the Tibetan Plateau[J]. Earth and Planetary Science Letters, 2014, 396: 88-96 doi:  10.1016/j.epsl.2014.04.003
    [57] Guo R M, Zheng Y, Xu J Q, et al. Transient Viscosity and Afterslip of the 2015 Mw 8.3 Illapel, Chile, Earthquake[J]. Bulletin of the Seismological Society of America, 2019, 109(6): 2567-2581 doi:  10.1785/0120190114
    [58] Freed A M. Afterslip (and only Afterslip) Following the 2004 Parkfield, California, Earthquake[J]. Geophysical Research Letters, 2007, 34(6): L06312
    [59] Jónsson S, Segall P, Pedersen R, et al. Post-Earthquake Ground Movements Correlated to Pore-Pressure Transients[J]. Nature, 2003, 424(6945): 179-183 doi:  10.1038/nature01776
    [60] Peltzer G, Rosen P, Rogez F, et al. Poroelastic Rebound Along the Landers 1992 Earthquake Surface Rupture[J]. Journal of Geophysical Research: Solid Earth, 1998, 103(B12): 30131-30145 doi:  10.1029/98JB02302
    [61] Wang R J, Kümpel H J. Poroelasticity: Efficient Modeling of Strongly Coupled, Slow Deformation Processes in a Multilayered Half-Space[J]. Geophysics, 2003, 68(2): 705-717 doi:  10.1190/1.1567241
    [62] Yang H Z, Guo R M, Zhou J C, et al. Transient Poroelastic Response to Megathrust Earthquakes: A Look at the 2015 Mw 8.3 Illapel, Chile, Event[J]. Geophysical Journal International, 2022, 230(2): 908-915 doi:  10.1093/gji/ggac099
    [63] Hu Y, Bürgmann R, Freymueller J T, et al. Contributions of Poroelastic Rebound and a Weak Volcanic Arc to the Postseismic Deformation of the 2011 Tohoku Earthquake[J]. Earth, Planets and Space, 2014, 66: 106 doi:  10.1186/1880-5981-66-106
  • [1] 王啸, 邹蓉, 李瑜, 王琪.  北斗二代导航卫星系统地壳运动监测能力 . 武汉大学学报 ● 信息科学版, 2022, 47(6): 934-945. doi: 10.13203/j.whugis20210316
    [2] 江颖, 刘子维, 张晓彤, 张丽娜, 韦进.  2021年青海玛多Mw 7.4地震前后b值的变化特征研究 . 武汉大学学报 ● 信息科学版, 2022, 47(6): 907-915. doi: 10.13203/j.whugis20220071
    [3] 高源, 瞿伟, 张勤, 王庆良, 郝明.  GNSS揭示的汾渭盆地及周缘现今地壳运动与应变差异 . 武汉大学学报 ● 信息科学版, 2021, 46(7): 1063-1070. doi: 10.13203/j.whugis20190365
    [4] 顾国华, 王武星.  GPS测得的2018年夏威夷6.9级地震与火山喷发地壳运动 . 武汉大学学报 ● 信息科学版, 2019, 44(8): 1191-1197, 1204. doi: 10.13203/j.whugis20180463
    [5] 袁鹏, 孙宏飞, 秦昌威, 张理想.  安徽CORS参考站三维速度场分析 . 武汉大学学报 ● 信息科学版, 2016, 41(4): 535-540. doi: 10.13203/j.whugis20140552
    [6] 李冲, 李建成, 黄瑞金, 谭理.  考虑欧拉矢量的BP神经网络模型建立区域地壳运动速率场 . 武汉大学学报 ● 信息科学版, 2014, 39(3): 362-366. doi: 10.13203/j.whugis20120046
    [7] 丁开华, 许才军, 邹蓉, 胡友健.  利用GPS分析川滇地区活动地块运动与应变模型 . 武汉大学学报 ● 信息科学版, 2013, 38(7): 822-827.
    [8] YinMyoMinHtwe, 申文斌, 孙榕.  缅甸Yangon及其周边地区地震危险性评估 . 武汉大学学报 ● 信息科学版, 2010, 35(4): 463-466.
    [9] 顾国华, 王武星, 孟国杰, 徐岳仁.  GPS测得的汶川大地震前后的地壳运动 . 武汉大学学报 ● 信息科学版, 2009, 34(11): 1336-1339.
    [10] 郭东美, 李军, 熊熊, 刁法启.  地壳运动速度场的数据融合研究 . 武汉大学学报 ● 信息科学版, 2008, 33(2): 212-215.
    [11] 陈俊勇, 郭春喜, 程传录, 张世娟.  利用GPS技术测定珠峰高程变化及其北坡地壳形变的初步分析 . 武汉大学学报 ● 信息科学版, 2007, 32(1): 1-5.
    [12] 吴云, 孙建中, 乔学军, 王辉.  GPS在现今地壳运动与地震监测中的初步应用 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 79-82,136.
    [13] 蒋志浩, 张鹏, 秘金钟.  中国新疆及其天山地区地壳运动特征分析 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 142-145.
    [14] 蒋志浩, 张鹏, 李毓麟, 秘金钟.  基于新疆GPS联测资料分析青藏高原北缘地壳运动规律 . 武汉大学学报 ● 信息科学版, 2003, 28(4): 417-421.
    [15] 王泽民, 刘经南.  球面坐标系非连续变形分析的数学模型 . 武汉大学学报 ● 信息科学版, 2001, 26(3): 209-212.
    [16] 刘经南, 施闯, 许才军, 姜卫平.  利用局域复测GPS网研究中国大陆块体现今地壳运动速度场 . 武汉大学学报 ● 信息科学版, 2001, 26(3): 189-195.
    [17] 许才军.  大地测量联合反演理论和方法研究进展 . 武汉大学学报 ● 信息科学版, 2001, 26(6): 555-561.
    [18] 王泽民.  非连续变形分析与现代地壳运动研究 . 武汉大学学报 ● 信息科学版, 2001, 26(2): 122-126.
    [19] 刘经南, 施闯, 陈俊勇.  92'、96国家高精度GPS A级网整体平差结果分析与我国块体运动模型研究 . 武汉大学学报 ● 信息科学版, 1998, 23(4): 314-319.
    [20] 宋成骅, 许才军, 刘经南, 蔡宏翔.  青藏高原块体相对运动模型的GPS方法确定与分析 . 武汉大学学报 ● 信息科学版, 1998, 23(1): 21-25.
  • 加载中
图(4)
计量
  • 文章访问数:  322
  • HTML全文浏览量:  95
  • PDF下载量:  93
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-11
  • 刊出日期:  2022-06-05

卫星大地测量成像地震周期形变研究综述

doi: 10.13203/j.whugis20220166
    基金项目:

    中国科学院战略性先导科技专项 XDB41000000

    武汉引力与固体潮国家野外科学观测研究站开放研究基金 WHYWZ202111

    作者简介:

    郭汝梦,博士,副研究员,研究方向为地震周期断层形变。guorm@apm.ac.cn

    通讯作者: 孙和平,博士,研究员,中国科学院院士。heping@whigg.ac.cn
  • 中图分类号: P315; P22

摘要: 近年来,卫星大地测量技术的快速发展为精确测量地壳形变和断层行为提供了前所未有的多维观测数据。它与地震学的结合使得地震周期形变监测的时空分辨率大大提升,为更加深入地研究地震周期过程和机理提供了一个窗口,地震大地测量学应运而生。它能够对地壳运动进行定量描述、对断层活动进行精准建模,从而为洞悉整个地震周期过程的应力应变演化提供科学依据,同时为评估地震危险性、实现地震预测预警提供科学指导。以卫星大地测量观测探究断层形变为主线,分析了断层处于地震周期不同阶段的运动学特征(震间、同震和震后),回顾了地震大地测量学在震源物理方面的一些重要发现。研究表明,利用卫星大地测量数据判定断层所处地震周期的阶段是实现地震预测的可行思路。

English Abstract

郭汝梦, 杨浩哲, 汤雄伟, 张文婷, 徐晓雪, 刘德川, 孙和平. 卫星大地测量成像地震周期形变研究综述[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 799-806. doi: 10.13203/j.whugis20220166
引用本文: 郭汝梦, 杨浩哲, 汤雄伟, 张文婷, 徐晓雪, 刘德川, 孙和平. 卫星大地测量成像地震周期形变研究综述[J]. 武汉大学学报 ● 信息科学版, 2022, 47(6): 799-806. doi: 10.13203/j.whugis20220166
GUO Rumeng, YANG Haozhe, TANG Xiongwei, ZHANG Wenting, XU Xiaoxue, LIU Dechuan, SUN Heping. A Review on Satellite Geodesy Applied to Image the Earthquake Cycle Deformation[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 799-806. doi: 10.13203/j.whugis20220166
Citation: GUO Rumeng, YANG Haozhe, TANG Xiongwei, ZHANG Wenting, XU Xiaoxue, LIU Dechuan, SUN Heping. A Review on Satellite Geodesy Applied to Image the Earthquake Cycle Deformation[J]. Geomatics and Information Science of Wuhan University, 2022, 47(6): 799-806. doi: 10.13203/j.whugis20220166
  • 卫星大地测量技术主要包括两类,一类是全球卫星导航系统(global navigation satellite system,GNSS),其进一步演化出精密卫星定位,如全球定位系统(global positioning system,GPS),在20世纪80年代开始用于地壳形变监测[1];另一类是合成孔径雷达(synthetic aperture radar,SAR),其进一步发展为合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)技术,在20世纪90年代开始用于地壳形变场追踪[2-3]。在过去十多年中,随着SAR卫星数量的激增,特别是Sentinel-1 SAR卫星的发射,数据量迅速增加;与此同时,InSAR数据处理技术的快速发展大大提高了精确捕捉长波长和小幅度地壳形变信号的能力,能够从大气延迟、轨道误差、非构造形变以及解缠误差等众多误差源中有效分离出毫米级构造变形[4]。此外,不同机构布设的GPS站点也越来越多,台站间隔越来越小[5],且不同大地测量观测融合技术越来越成熟,使得获取全天候、高精度、大范围、高时空分辨率的3D地壳形变场成为可能[6]

    中国是地震频发国家,强震危险性的可靠预测一直是防震减灾事业的急迫需求,了解地震周期各个阶段(震间、同震和震后)构造活动的基本过程有助于减轻日益城市化的人口暴露于此类灾害的风险[7-8]。文献[9]指出强震的孕育、发生和发展始终伴随着地表形变现象,如图 1所示。卫星大地测量技术能够定量、实时或准实时地捕获地震周期过程的地壳形变信号,精确追踪断层活动的时空演化过程,有助于理解构造断裂带的分段运动特征,洞悉整个地震周期过程的应力应变累积-释放机制,同时为评估地震危险性、实现地震预警预报提供科学指导[10]

    图  1  地震周期断层形变图(改自文献[9])

    Figure 1.  Map of Earthquake-Cycle Fault Deformation (Modified from Reference [9])

    • 震间形变过程相对比较缓慢,可以持续几十年、几百年甚至上千年。在此过程中,构造应力持续加载,地壳相关物理量如位移、重力等发生缓慢变化。当应变积累导致应力超过断层摩擦所能承受的临界点时,地震发生。通过震间过程的研究,不仅可以对地震的孕育机制、发震模式有更深刻的认识,为中长期地震潜力评估提供科学依据,还可以对地壳结构特征、地球交互模式有更加深入的了解,从而为相关地球物理现象提供科学的解释和有力的证据[11]

      利用长期的卫星大地测量数据解析断层的运动状态及其所处孕震周期的阶段是实现地震预测的可行思路[12-13]。已有观测和研究表明,在震间的不同阶段,断层两侧的形变模式往往不同。走滑断层不同孕震阶段形变模式如图 2所示,其中t表示上次强震距今时间,T表示地震复发间隔,t/T表示走滑断层所处的震间阶段。由图 2可以看出,早期(t/T=0.05)断层两侧的形变梯度较大,而中(t/T=0.10)、后期(t/T=0.95)则比较平缓。因此断层周边的地壳速率场被广泛用于推测断层的地震潜力[14-16]。目前大多数研究都是基于经典的弹性螺旋位错理论估计断层的滑移率以及闭锁深度[17-19]。文献[20]构建了一个弹性的双断层模型,采用长期的GPS震间形变解析2021年玛多Mw 7.4地震孕震断层的运动学特征,利用贝叶斯蒙特卡洛反演方法同时估计了东昆仑断裂和昆仑山口-江错断裂的滑移率分别为1.7±0.8 mm/a和7.1±0.3 mm/a。这些简单的模型至今仍被广泛使用,主要原因在于其简单性和有效性,也归因于以往震间数据较大的不确定性[21]。目前已经广泛达成的共识是地球在地震周期过程中为粘弹性的[21-22],因此,忽略软流圈粘弹性松弛的效应以及地震周期性的载荷将导致在地震周期不同阶段断层滑移率被系统性地低估或高估,并产生其他运动学参数估计的偏差[15-16]。当前,大面积、毫米级精度的地壳形变数据集具备前所未有的空间尺度和分辨率,为构建更加复杂的、符合真实地球介质的粘弹性地震周期模型提供了一个窗口,各种粘弹性地震周期模型被提出并用于解析不同断层系统的运动学与动力学特征[23-25]。文献[15]构建了一个弹性的上地壳覆盖一个粘弹性基底的双层粘弹性地震周期模型,用于探索中国大陆阿尔金断裂和沂沭断裂的运动特点以及地震潜力。文献[16]则采用了一个弹性板片覆盖粘弹性下地壳和上地幔的三层粘弹性地震周期模型,阐明了东昆仑断层系统托索湖段的滑移率为5.5±0.7 mm/a。

      图  2  走滑断层不同孕震阶段形变模式(改自文献[16])

      Figure 2.  Deformation Pattern of Transform Fault inDifferent Stages(Modified from Reference [16])

    • 地震的本质是断层面上长期的弹性应变累积被快速释放,导致断层两侧的块体发生错动。该过程一般持续几秒、几十秒甚至几百秒,覆盖范围一般从十几千米至数百千米。大地震的时空破裂过程通常指断层面上有限个凹凸体或障碍体的先后破裂[26],破裂行为有时十分复杂,如多条共轭断层先后破裂[27]、破裂在断层之间跳跃[28]、俯冲带地震的破裂折返[29]以及震源频率特征随深度变化[30]等。解析精确的同震破裂过程可以为地面震动预测、震后机制识别以及应力传输估计等提供驱动力,也可以有效地破译孕震断层结构[31-33],其依赖于有限断层成像理论的快速发展。通常的做法是在空间上将主震事件离散成许多尺度较小的子事件,即将发震断层面离散成多个子断层(矩形或三角形),反演子断层的震源时间函数,综合所有子断层的震源时间函数结果,考虑破裂时间上的连续性以及破裂速度的限制,反演得到整个地震在时间和空间(断层面)上的破裂特性,从而描述震源的时空破裂行为[34]

      美国地质调查局长期致力于追踪全球强震破裂过程,采用远震体波和面波进行快速有限断层反演,通常在震后1天之内发布同震滑移模型,对于灾后快速救援起着重要的指导作用[35]。然而,远震体波对于破裂时间更敏感,对于破裂位置却不敏感[36];面波频率较低,其滑移分辨率十分有限[37]。为了提升破裂过程的空间解析度,GPS和InSAR被广泛应用于有限断层反演中。与远震观测不同,大地测量捕捉的地壳形变数据记录的是震源区地震事件导致的静态形变,可以很好地约束断层几何、滑移大小以及方向。GPS台站密度有时比较稀疏,在有限断层反演过程中通常需要加入比较强的平滑权重来提升稳定性,而InSAR空间分辨率比较高,但观测精度略低,因此同时采用两种数据可以达到最优的空间和地震矩解析度[38]。文献[39]基于GPS和InSAR数据联合反演了2008年汶川Mw 7.9地震的几何形态和破裂空间分布,揭示了其分段梯次破裂的特征。尽管这些基于卫星大地测量资料反演得到的震源滑移模型与余震分布、地表破裂的结果保持了较好的一致性,有力地证明了其良好的位移解析度,但是大地测量数据只包含同震形变信息,没有任何时间解析度[40-41]。另外,由于受地理位置、植被覆盖等的影响,大地测量观测方位角的分布有时并不理想(如俯冲型地震的大地测量数据通常位于陆地上,只对俯冲板块边界下部的位错比较敏感),甚至在某些区域没有大地测量观测。图 3为基于不同数据有限断层反演结果的时间、空间以及地震矩分辨率,表明地震数据可以很好地弥补这些缺陷。因此联合地震波形数据和大地测量资料可以精确解析大地震的时空破裂过程,提高模型的稳健性,同时恢复破裂的空间位置、时间演化以及位移量3个维度的属性。

      图  3  基于不同数据有限断层反演结果的时间、空间以及地震矩分辨率(改自文献[37])

      Figure 3.  Temporal, Spatial, and Moment Resolutions of Fault Inversion Results by Using Different Datasets (Modified from Reference [37])

    • 大地震同震破裂会造成震源区岩石圈尺度的应力场扰动,导致震后持续的地壳形变,主要包括断层面上持续的无震滑动(余滑)、下地壳和上地幔的粘弹性松弛以及浅地壳流体运动引起的孔隙弹性响应[42-44]图 4为震后形变概念示意图。由图 4可以看出,大地震发生以后,多种震后效应耦合在一起,但不同成因的震后形变具有不同的时空分布特征,取决于地震的破裂类型、地壳结构和断层面的岩性等。大地测量方法可以有效观测到这些形变造成的地表位移,并用来分离不同机制的震后贡献。通过震后形变的研究可以了解断层面的摩擦性质、软流圈的流变结构以及孔隙介质属性,也可以加深对断层重新闭锁以及应力应变重新累积的理解[45]

      图  4  震后形变概念示意图

      Figure 4.  Schematic Diagram of Postseismic Deformation

      震后余滑是同震应力变化造成的无震滑移,通常发生在同震破裂区上倾或下倾具有速度强化特性且未破裂的断层面[46]。余滑主要控制震后数年近场的快速形变,对于理解断层摩擦行为具有重要指导作用。文献[47]利用震后半年GPS数据反演了2017年墨西哥Mw 8.2地震的震后余滑,发现其累积释放的地震矩相当于Mw 8.0事件,且与同震高滑移区域呈互补的模式,揭示了余震的时间演化过程可能主要由余滑控制。文献[48]通过2021年玛多Mw 7.4地震震后两个月的InSAR观测,得出了类似的结论。但这种直接反演的运动学余滑模型通常忽略了物理的合理性,可能会出现难以解释的凹凸体和过度拟合现象。针对该问题,基于速率-状态摩擦本构关系的应力驱动余滑模型应运而生[49],并成功解释了许多地震震后余滑的时空演化[50]。文献[51-52]基于大地测量观测资料研究了2015年Gorkha Mw 7.8地震早期的震后余滑,指出余滑主要发生在同震破裂面下方的速度强化区。理论和数值模拟表明,速度弱化区域通常不会表现出无震余滑,而速度强化区域在震间阶段不会累积应力。但是文献[53-54]发现限制余滑发生在非凹凸体的范围不能完美地拟合震后形变,提出了对“速率-状态凹凸体模型”的质疑,并且越来越多的学者发现余滑与大地震周期形变研究同震的高滑移区(速度弱化区)存在重叠。

      粘弹性松弛主要来自较弱的下地壳或上地幔缓慢释放同震应力的加载而产生的地壳形变,该阶段持续时间比较长,一般为几年、几十年甚至几百年,在空间上主要作用于地震的中远场[42]。在不同的构造条件下,基于M≥7地震事件震后数月乃至数年的大地测量观测可以推断下地壳和上地幔的粘滞系数,为了解研究区域的断层耦合及流变特性提供关键基础,同时对认识岩石圈的介质结构属性以及地壳演化历程提供重要依据。文献[55]证明了上地幔中的粘弹性松弛在海沟附近产生向陆地运动的重要作用,推断地幔楔瞬态的Kelvin粘度为2.5×1017 Pa·s,稳态的Maxwell粘度为1.8×1018 Pa·s。文献[56]基于2008年汶川地震接近2年的震后观测发现近场的地表形变主要由余滑控制,而中远场的形变主要来源于粘弹性松弛的贡献,并约束了青藏高原瞬态的Kelvin粘度为7.9×1017 Pa·s,稳态的Maxwell粘度为1×1018 Pa·s。文献[57]给出了2015年Illapel Mw 8.3地震震后一个半月下地壳和上地幔的粘度分别为6.0×1017 Pa·s和9.0×1017 Pa·s,指出了以前忽略粘弹性效应的纯弹性余滑模型对破裂区的上倾余滑估计不足,对下倾余滑估计过度。

      孔隙弹性回弹是同震破裂导致断层附近孔隙压力梯度的改变后,流体压力达到再平衡过程引起的地壳形变,其持续时间较短,变形相对较小,主要分布在断裂带近场上地壳[58-59]。在大多数震后机制的研究中,孔弹性效应经常被忽略,可以通过在弹性位错模型中使用泊松比的未排水值和排水值预测的位移场之间的差异来估算。文献[60]利用泊松比在排水(0.27)和未排水(0.31)状态下的差异估算了1992年Landers Mw 7.3地震的震后孔隙回弹机制的贡献,发现余滑和孔隙回弹的综合作用可以很好地解释其震后形变。虽然这种方法给出了孔隙弹性位移的大小和方向,但其复杂性被忽略,孔隙弹性形变的时空演化也未知。文献[61]发展了另外一种方法来模拟由地震或抽注水引起的时间相关的孔隙回弹形变。在正演模型中,格林函数基于一阶多层的弹性结构和深度相关的水力扩散系数生成,然后使用同震滑移模型来计算地震引起的孔隙弹性松弛。文献[62]利用该方法模拟了2015年Illapel地震震后的孔弹性位移,发现忽略孔弹性效应的纯余滑模型高估了同震破裂带南北两侧的余滑,低估了同震破裂带中部的余滑,局部误差高达50%。此外,文献[63]通过构建有限元模型研究了孔隙回弹机制对2011年Tohoku Mw 9.0地震震后形变的影响,研究表明孔隙回弹机制对断层破裂周围区域的变形有着不可忽略的贡献。

    • 地壳形变是断层运动的直接表现形式之一,利用地壳速率场的动态演化和地震震源的过程跟踪是实现地震成功预测的可行思路。人类正在进入一个激动人心的地球观测数据新时代,卫星大地测量技术可以提供毫米级的全球范围尺度的大规模数据。如何利用大尺度形变场分布、活动块体运动特征以及断层动力学信息,识别出紧迫程度相对较高的有限个潜在强震源,判定断层所处的孕震阶段,是强震预测亟待解决的问题。通过对全球范围内活跃断层长期地壳速率场的建模,人们在基本断层过程、地震周期形变的理解方面也取得了较大进步。震源物理机制的研究正从弹性向粘弹性发展,从地震周期内某一阶段的研究向地震孕育、发生、调整的全过程解析及模拟发展。虽然本文专注于卫星大地测量观测平台,但其与地震观测、实地勘测以及航空测绘等有很强的互补联系,如地震观测在全球范围内覆盖了所有的中强地震,可以提供必要的高分辨率时间采样,以观察地震破裂的快速演变,这是静态测量不可能实现的。因此,为了加强对地震灾害和构造过程的理解,通常需要将卫星大地测量和其他技术相结合,进而提升防震减灾的能力。

      以上的回顾仅为卫星大地测量技术在成像地震周期断层形变的一些成果,并不全面,有一定主观成分,但是从中可以反映出地震大地测量学在过去几十年的发展与进步。

参考文献 (63)

目录

    /

    返回文章
    返回