中国地震遥感应用研究与地球物理场探测卫星计划[PDF全文]
申旭辉, 张学民, 崔静, 周新, 姜文亮, 龚丽霞, 李永生, 刘芹芹
摘要: 中国地震遥感应用研究始于20世纪70年代中期,从简单的多光谱遥感图像解译线性活动构造开始,迄今已走过了40余年的历程。其中20世纪70、80年代以可见光遥感图像应用为特征,90年代为红外遥感、GNSS遥感和InSAR技术试验应用阶段,20世纪末到21世纪前10年主要表现为综合遥感应用和卫星电磁、卫星重力等地球物理场探测卫星计划预研,2011年以来侧重于集成前期研究成果,实施地球物理场卫星计划和卫星地震应用系统建设运行。目前可见光遥感技术在活动构造探查和灾害评估中已实现业务化应用,卫星红外、卫星电磁、InSAR技术、卫星重力和高光谱气体地球化学探测技术在地震监测领域越来越凸显出其应用效能。作为国家地球物理场探测卫星计划首发星,张衡一号电磁监测试验卫星于2018年2月2日进入预定轨道,目前在轨运行正常并产出了中国首幅全球地磁图,填补了中国在全球地球物理场获取能力方面的空白。
关键词: 地震遥感     地球物理场探测卫星     卫星电磁     卫星重力     InSAR     红外遥感     光学遥感     GNSS遥感    
DOI: 10.11834/jrs.20188337    
收稿日期: 2018-06-08
中图分类号: TP73    文献标识码: A    
作者简介: 申旭辉,1965年生,男,研究员,研究方向为遥感地学应用和活动构造。E-mail:shenxh@seis.ac.cn
通信作者: 刘芹芹,1987年生,女,博士后,研究方向为地球物理和高光谱遥感。E-mail:liuqq2008@163.com.
基金项目: 政府间国际科技创新合作重点专项(编号:2016YFE0122200);卫星团队建设项目
Remote sensing application in earthquake science research and geophysical fields exploration satellite mission in China
SHEN Xuhui, ZHANG Xuemin, CUI Jing, ZHOU Xin, JIANG Wenliang, GONG Lixia, LI Yongsheng, LIU Qinqin
1. The Institute of Crustal Dynamics, China Earthquake Administration, Beijing 100085, China
2. School of Environment, Tsinghua University, Beijing 100084, China
3. Institute of Earthquake Forecasting, China Earthquake Administration, Beijing 100036, China
Abstract: The application of seismic remote sensing in China began in mid 1970s. More than 40 years have passed since the beginning of the simple linear active structure of multi spectral image interpretation. In the 1970s and 1980s, it was characterized by the application of visible light remote sensing images, and in the 1990s, the infrared remote sensing, GNSS remote sensing, and InSAR technology were at the experimental application stage. From the late 20th century to the first decade of the 21st century, it was mainly manifested as integrated remote sensing application and pre-research of the satellite program of geophysical field exploration, such as electromagnetic field and gravity satellites. Since 2011, more attentions were played to integrating previous research results, implementing the geophysical field satellite plan, and the constructing, as well as operating the satellite seismic application system. After more than 40 years of development, the application of visible-light remote sensing technology in active tectonic exploration and disaster assessment has been realized. Satellite infrared, satellite electromagnetism, InSAR technology, satellite gravity, and hyperspectral gas geochemical detection technology have become increasingly prominent in the field of seismic monitoring. As the first satellite of the national geophysical field exploration satellite program, China Seismic-Electromagnetic Satellite, named ZH-1, entered its present orbit on February 2, 2018. At present, China’s first global geomagnetic map was produced after the normal operation in orbit, which fills the gap in the acquisition capacity of the global geophysical field in China.
Key Words: seismic remote sensing    geophysical field detection satellites    satellite electromagnetic    satellite gravity    InSAR    infrared remote sensing    optical remote sensing    GNSS remote sensingremote sensing    
1、引 言

中国是世界上地震灾害十分严重的国家之一。据《后汉书·张衡传》记载:早在公元132年,张衡制造了世界第1台可感知地震方向的候风地动仪,并准确记录了公元138年陇西大地震的方向和时间。自19世纪末期现代地震学的诞生以来,中国逐渐建成了世界上规模最大、覆盖学科门类最广的地面地震观测系统。经历了1975年海城地震成功预报的喜悦,1976年唐山地震和2008年汶川地震不期而至的打击,地震预报这一世界性科学难题至今仍然未能被攻克。而伴随着社会经济高速发展的需要,防震减灾工作也由逐渐整合形成为监测预测、震害防御和应急救援的3大工作体系。在防震减灾3大工作体系业务能力建设和科技支撑水平发展过程中,遥感技术以其高动态、大观测、全天时、全天候等特点,展现出显著的技术优势。卫星遥感等空间技术手段,空间对地观测各种业务手段被大量用于防震减灾各业务体系、各业务环节,且在防震减灾各业务领域中展现了广阔的应用前景,成为推进地震监测预测和应急救援能力提升的重要技术手段(申旭辉 等,2007)。

2、地震遥感应用研究历程的简要回顾

1974年,中国著名地震学家顾功叙院士和丁国瑜院士率团访问美国地质调查局(USGS),获赠了一套中国多光谱卫星遥感影像图,由此拉开了中国遥感应用研究的序幕(丁国瑜和李永善,1979)。经历40多年的发展,目前地质遥感和地貌遥感方法已经成为地震科学研究的基本手段,灾害遥感技术发展为地震应急救援辅助决策和灾害评估发挥了重要作用,以卫星红外、卫星电磁、InSAR、GNSS为代表的卫星遥感技术已经在地震监测预测研究中发挥重要作用。

    (2.1) 地震遥感应用研究起步阶段(20世纪70年代中期到80年代末期)

从1974年开始,中国地震学家开展了全国性活动构造遥感调查研究工作,并随即引进了国外大型遥感图像分色和融合设备,光学遥感信息解译成为全国活动构造普查、线性构造分析和全国地震区划图编制的重要手段,推进了地震构造研究的快速发展。中国著名地震学家丁国瑜院士(1980)认为:“由于大量新资料的积累,特别是宇航遥感资料的急速增长,得到了许多以前在近距离所看不到的现象和新认识,对线性构造在地壳表面大范围的分布特征及其活动发展动态的研究进入了一个新阶段。”,黄培华等人(1982)编辑出版《地震地质学基础》,明确陆地卫星图像应用是活动构造研究的5种主要方法之一。

在此期间,以遥感技术应用为基础,编制了《中国活动构造典型卫星影像集》(国家地震局地震研究所 等,1982)和《中国卫星影象地震构造判读图》(国家地震局地质研究所,1985),初步查明了中国活动构造和地震构造分布,为地震区划和区域活动构造研究奠定了重要基础。

    (2.2) 多种遥感手段全面探索推进阶段(20世纪80年代末期到20世纪末期)

1988年,苏联科学家发现了地震前红外辐射增温现象。中国学者迅速跟踪国际前沿,开始了这方面的研究工作。围绕震前红外辐射增强现象及其统计特征、地震红外辐射机理以及红外辐射信号提取技术等方面开展了大量的研究工作,在地震系统建设了多个NOAA/MODIS卫星数据接收站,以保障红外遥感数据的时效性和连续性,并在历年地震监测预测实践中发挥相应的作用。

1993年,欧洲空间局专家反演获取了1992年美国LANDERS地震同震形变,中国学者王超等人(2000)单新建(2002)张景发(2002)等随后利用国外SAR遥感数据反演获得了中国历次7级以上地震的同震变形并分析了相应的地震变形特征。为进一步提高InSAR形变测量精度,中国学者分别在青藏高原东北缘地区、京西北地区、云南丽江和鲜水河地区以及西藏当雄—羊八井地区,布设了人工角反射器台阵,发展了毫米级形变测量和反演技术。

这一阶段,高分辨率遥感信息越来越丰富,王晓青等人(2013)利用高分辨率光学遥感信息发展了地震灾害遥感评估技术并在历次破坏性地震中得到检验。谢广林(1991)申旭辉等人(2000)利用高分辨率遥感数据定量分析断层位移规模,并识别和寻找历史地震形变带,申旭辉等人还利用SAR数据的穿透性能在西北干旱地区开展了隐伏构造提取试验研究,多源遥感应用在地震行业得到长足发展并全面覆盖防震减灾3大业务领域。

    (2.3) 地球物理场遥感探测与多源遥感应用技术协同推进阶段(21世纪初)

20世纪末至21世纪初,国际学术界围绕地震能否预报展开了一场大讨论(Geller 等,1997Wyss 等,1997陈运泰,2009),同时一系列7级乃至8级以上强烈破坏性地震袭击了世界各地,困惑中的地震预报与巨大的地震灾难触动了关于地震预报方法论的思考。在此背景下,中国地震局等正式启动地震卫星规划论证工作,并相继发布了《国家防震减灾规划(2006—2020)》、《民用航天“十二五”规划》、《防震减灾“十二五”规划》。地震遥感应用研究也由分散性课题研究开始整合,形成了以国家科技支撑计划、863计划、国际合作交流计划和民用航天科研计划为依托的系统研究分工合作新格局。

此阶段地震遥感应用的发展,一是地球物理场遥感探测技术得到高度关注,卫星电磁和卫星重力探测、反演与应用技术研究和研发取得突破性进展;二是灾害遥感应用进入业务化,2004年伽师震群后,遥感技术作为一种重要的技术手段进入应急救援与灾害评估业务领域,且在汶川地震、玉树地震灾后应急救援业务发挥了重要作用;三是多源遥感综合应用不断深入和拓展,高光谱气体地球化学、LiDAR技术以及红外多角度遥感等新技术开始进入防震减灾领域。

    (2.4) 遥感应用业务化示范与地球物理场遥感探测系统发展阶段(2011年以来)

2010年,“高分地震遥感监测与应急应用示范工程”正式启动,旨在建设具备业务化运行能力的地震遥感示范应用系统,2013年国务院正式批准电磁监测试验卫星工程研制,以高分重大科技专项和电磁监测试验卫星工程为依托,建成了卫星地震应用示范系统,并具备业务化运行能力。2018年2月2日,中国地球物理场探测卫星计划首发星张衡一号顺利发射入轨,目前在轨运行正常,并首次获取了中国自主知识产权的全球地磁图,填补了中国在地球物理探测卫星和全球地球物理场战略资源获取能力的空白。

2008年汶川地震后,国家启动了623专项论证工作,并于2015年正式发布了《国家民用空间基础设施中长期发展规划(2015—2025)》,按照规划和《卫星对地观测数据产品分类分级规则标准》,正式提出和明确了地球物理场探测卫星的概念。目前,除电磁监测试验卫星之外,其他地球物理探测卫星计划同步论证推进中。

3、地震遥感应用研究现状与趋势     (3.1) 遥感技术在地震监测预测领域的应用研究

地震遥感应用现状及其总体趋势归纳总结见图1,主要从InSAR形变反演技术、卫星红外遥感技术、卫星电磁探测技术以及其他遥感技术应用等方面进行概述总结。

图 1 地震遥感应用现状 Figure 1 Application status of seismic remote sensing
         3.1.1. InSAR形变反演技术

1993年,欧洲科学家利用ERS-1卫星数据,首次获得了1992年美国兰德斯7.3级地震同震形变场,引起国际科学界轰动(Massonnet 等,1993Zebker 等,1994),InSAR技术作为一种低成本、高分辨率的空间对地测量手段越来越得到学术界的认可。

(1) 同震形变提取。D-InSAR技术在中国地震同震形变中的应用研究始于1997年玛尼7.5级地震(单新建 等,2002),此后,1998年张北Ms6.2级地震,2001年昆仑山口西Mw8.1 级地震(单新建 等,2004),2008年改则Mw6.4级地震(孙建宝 等,2008),2010年玉树Mw7级地震(屈春燕 等,2013),D-InSAR均给出了快速可靠的同震形变信息,为强震考察和相关科学研究提供了重要依据。

(2) 3维形变反演。洪顺英等人(2010)利用InSAR、GPS、像素偏移Offset-tracking方法以及形变模拟获取了2008年1月9日和16日改则Mw6.4和Mw5.9双震,3月21日于田Mw7.1地震和2003年12月26日伊朗巴姆Mw6.9地震3维形变特征。此外,随着SAR卫星技术的不断发展,GNSS观测能力的不断提高,基于GNSS和InSAR技术结合的地壳3维形变场观测,将是未来地震地壳形变监测的核心手段。

(3) 缓慢微小变形监测。由于雷达波在通过电离层时的大气效应以及受数字高程模型精度的影响,D-InSAR技术形变结算精度总体上在厘米量级,距离长期缓慢地壳形变测量需求存在明显不足。为了最大限度提高D-InSAR地壳形变解算精度,Ferretti等人(2001) 提出了基于永久干涉体PS(Permanent Scatters)测量方法以获取毫米级变形精度。中国先后在海原、六盘山、西秦岭北缘、羊八井地区、丽江断裂等地区布设了永久角反射器,基于InSAR技术的震前形变监测技,在当雄—羊八井以及玛尼地区的试验均取得了良好效果。

         3.1.2. 卫星红外遥感技术

20世纪80年代末期,苏联科学家在分析中亚地区地震活动时,发现了震前红外辐射增强现象(Gorny 等,1988)。此后许多科学家集中开展了地震红外异常分析研究,徐秀登等人(2000)邓明德(1997)初步研究了地震红外辐射的机理,刘德富和康春丽等人(2005)康春丽等人(2003)开展了大量震例分析研究。

(1)红外遥感机理研究。国内外学者通过空气电场实验、岩石加载实验等手段开展了积极探索研究,提出了多种机理解释。主要有:气热说(Gorny 等,1988Qiang 等,1995)、应力致热(Wu 等,2000吴立新,2001)、地下流体作用(Tronin 等,2002)、“P-hole”理论(Freund,2007)、氡气衰变释放潜热(Pulinets,2007)、地下热传导(蔡永恩 等,1987)、断层蠕动摩擦生热(郭子祺,1998)、微震声波振动促使土壤热通量变化(王纯莹 等,2005)等假说。

(2)红外遥感前兆统计特征分析。近年来,利用卫星热红外遥测技术研究断层活动以及地震前兆受到国内外学者的广泛关注。美国NASA利用EOS/MODIS卫星热红外数据研究2001年1月26日印度古吉拉特邦Ms7.6地震时,发现地震前5 d在印度西部震中及外围出现了大面积热红外异常。中国学者的系统研究结果表明,震前孕震所处地域存在明显的热辐射增强现象,通过卫星遥感技术所产出的数据产品—长波辐射场强度(OLR值)和亮度温度(TBB)的变化识别出这种增强,对这种震前红外遥感异常特征进行了分析。近年发生的几次显著性地震事件,如2005年江西九江地震、2006年河北文安地震、2007年云南普洱地震、2008年汶川8级地震等,卫星红外遥感信息都有明显异常显示。

目前,利用卫星红外遥感资料进行地震监测预报方面已取得了一定的进展,对1975年以来中国全部6级以上地震震例进行了全面整理,取得了可喜的成果。其中部分研究成果已应用于日常的地震预报工作之中,成为日常地震监测预报的重要手段之一。

(3)数据分析与处理技术研发。地震热红外的识别算法和应用也逐渐成熟,数据源也从单一影像发展到多源遥感相结合(Qin 等,2012),分析也从水平方向向垂直空间领域发展,异常识别方法也逐步发展。异常提取方法上主要包括4大类:1)目视解译(徐秀登 等,2000);2)基于差值分析的异常提取算法,典型算法包括震前震后亮温差值法、断裂带内外亮温差值法、涡度法等(Saraf和Choudhury,2005);3)基于信号分析的异常提取算法,典型算法有小波、功率谱、夜间热梯度(NTG)等(Saradjian和Akhoondzadeh,2011);4)基于背景场分析的异常提取算法,典型算法包括历年同期偏移指数法、RST算法等(Pergola 等,2010)。在多参数异常分析应用上,Wu等人(2012a2012b)提出了LCA(Lithosphere Coversphere Atmosphere)耦合模型以揭示多参数多圈层准同步变化之间相互关系。卫星红外遥感应用于地震相关信息研究,需考虑到红外遥感信息影响因素多样,以及采用不同波段的信息如中红外、远红外与微波辐射,进行综合分析研究。

         3.1.3. 卫星电磁探测技术

1965年,俄罗斯科学家在分析Alouette卫星电磁信号时,发现卫星记录到地震低频电磁辐射前兆现象。2004年6月,法国成功发射DEMETER卫星,专门用于研究与地震、火山喷发等有关的电磁异常和电离层扰动以及与人类活动有关的全球电磁环境。中国科学家重点基于中法、中俄和中意合作,开展了系统的卫星电磁技术研发和科学研究工作。

(1) 典型地震震例解剖。2008年汶川地震以来,中国科学家围绕地震电离层前兆问题开展了大量研究工作,并在单参量、单手段震例分析的基础上逐渐向多参量、天地联合等方向快速发展。针对2008年汶川特大地震进行了重点剖析,发现强震前卫星观测多频段电磁场、电离层等离子体密度、温度等多参量异常扰动,并同步在地基同类观测中出现类似扰动增强现象(张学民 等,2009Zhang 等,2009, 2010a)。对于国外发生的破坏性地震也开展了响应研究和震例积累,如2010年海地7.0级地震(Zhang 等,2014)、2010年智利7.9级地震(Zhang 等,2010b)、2005年印尼8.5级地震等(Zhang 等,2011),发现了震前1—7 d内电场与电子密度等多参量同步扰动现象。在卫星多参量综合分析方面,张学民等人(2014)利用地基电离层垂测、GPS TEC和DEMETER卫星探测资料分析了2010年玉树7.1级地震前的扰动现象,发现震前1 d多参量同步增强;同时针对这次地震,Shen等人(2017)对比卫星与地基记录的VLF电波信号,发现在震前人工源电波信号强度降低。Liu等人(2016)结合GPS TEC、DEMETER和DMSP 3种卫星资料,分析了印尼一次7.2级地震前后3种探测技术在相同时空内的响应关系。研究结果显示,多参量同步异常对于解析异常在不同高度的时空分布特征及其传播耦合关系具有重要的参考价值。

(2) 地震电离层前兆信息的统计分析。徐彤等人(2012)分析了中国14个大地震临界频率(f0F2),认为85.7 %的地震在震前1周存在异常,并且几种出现在地方时11:00—17:00。Le等人(2011)利用GPS TEC数据分析了2002年—2010年全球 736次 6.0级以上地震,表明GPS TEC的变化随震级、震源深度、发震时刻而变化。Yao等人(2012)研究了2010年全球7级以上地震GPS TEC,发现多出现于震前0—2 d和地方时12:00—20:00。刘静和万卫星(2014)利用GPS TEC对1998年—2010年中国6级以上强震进行了分析,结果显示震中周围多个方向均在震前3—5 d观测到异常扰动。Zhang等人(2012, 2014)利用DEMETER卫星资料分析了2005年—2010年全球7级以上地震前电场异常,结果表明:板块边缘地震更易于引起电离层扰动,伴随着ULF/ELF频段电场功率谱增加1—2个数量级,电场强度增加3—15 mV/m;He等人(2011)利用DEMETER卫星资料研究了2006年—2009年全球约7000 次5级以上地震前后电子浓度的变化,发现震前异常集中在震中附近;Li和Parrot(2013)统计分析了DEMETER卫星运行6年期间记录的离子密度参量与地震的关联,发现异常数目随地震震级增加有明显增多,且有地震当天异常数量最多;Zhang等人(2013)总结了2005年—2010年全球7级强震前的高能粒子沉降特征,发现70—330 keV的低能电子沉降与地震存在更好的关联性。

(3) 地震电离层耦合机理。前期研究结果显示(Pulinets和Boyarchuk,2004),地震电离层耦合主要以3种模式从震源传播到地表,穿过大气层向电离层传播(图2),其中包括化学—电场效应、声重波及电磁波辐射传播效应。赵庶凡等人(2011)Shen等人(2011)发展了甚低频电磁波传播的全波计算模型,并开展了岩石层—大气层—电离层耦合机制以及地震电磁辐射向电离层传播的模拟;Zhou等人(2017)优化发展了附加直流电场在电离层中传播耦合过程,揭示了不同方向电场在不同纬度的时空演化图像;Hao等人(2012)研究了2011年日本9级地震时的地震波、地磁场、TEC、次声波及多普勒多种观测技术的同震信号,发现地震面波激发的声重波向电离层传播扩散的特征。观测信息和理论模型研究方面的阶段性成果为发展地震电磁立体监测体系提供了科学支撑(张学民 等,2016)。

图 2 主要地震电磁电离层效应 Figure 2 Main electromagnetic ionospheric effects
         3.1.4. 其他遥感技术应用

(1) GNSS(Global Navigation Satellite System)遥感。1988年开始的滇西实验场GPS观测工作及随后开展的攀登计划,充分展示了GPS技术在监测现今地壳运动的前景。在此基础上,“九五”期间,由中国地震局牵头,中国建立了以地震预报为主要目的、以GPS观测技术为主的“中国地壳运动观测网络”(一期工程),已经建成覆盖全国的基本站28个,基准站56个,区域站1000个。目前中国地壳运动观测网络(CMONOC)二期工程已经得到国家发改委批准,即将启动实施。

目前,GNSS已经成为地震监测预报的重要手段。利用覆盖全国GNSS观测网络,已经能够获得一年尺度的全国地壳运动分布图,并通过国际合作,分析板块运动和地壳运动变化信息,为地震年度趋势预测提供重要的决策依据。利用GNSS连续观测站提供的信息开展地震前兆跟踪研究也取得了实质性进展,取得了震前地壳形变特征和电离层扰动信息的大量震例研究结果。如2001年11月14日发生在新疆维吾尔自治区与青海交界处的昆仑山口西M8.1大地震,利用GNSS取得了震前极为珍贵的地壳运动和电离层变化信息。

(2) 卫星重力探测技术。低轨道的重力卫星GRACE(Gravity Recovery and Climate Experiment)和GOCE (Gravity Field and Stead-state Ocean Circulation Explore)的发射、运行和观测从空间上提供了前所未有的高精度、高分辨率的全球重力场和大地水准面模型。GRACE卫星采用低—低卫星跟踪测距模式,能够获取大型地震产生的重力场变化(Sun和Okubo,2004)。目前GRACE能够检测出2004年苏门答腊Mw 9.3、2010年智利Mw 8.8和2011年日本东北Mw 9.0地震同震和震后黏滞性调整产生的重力场变化(Han 等,2006Panet 等,2007Cambiotti 等,2011Zhou 等,2012)。利用GRACE观测到的空间重力变化能够约束地震事件的震源参数(Han 等,2011Cambiotti 等,2011),也能够约束断层的同震滑移分布(Zhou 等,2018)。

2009年欧空局发射的GOCE卫星轨道约为255 km,比GRACE卫星轨道更低,搭载的重力梯度计能够通过测量地球的重力梯度恢复精度更高的中高阶(约50—250)重力场模型。GOCE卫星在轨运行期间(2009年3月—2013年11月)经历了全球两次特大地震——2010年智利Mw 8.8和2011年日本东北Mw 9.0。GOCE卫星可以当作空间地震仪,能够观测到大地震引起的大气密度波动信号 (Garcia 等,2013)。重力卫星GRACE和GOCE资料为震源提供了地球质量迁移的约束,加深了人们对震源物理的理解。

(3) 高光谱遥感技术应用。气体地球化学作为地震前兆研究的重要内容,自20世纪60年代中期以来引起了广大学者的重视,大量的地面观测数据和震例研究表明断层气体(CH4、CO2、He、Rn和Hg等)与地震活动、构造有关,对地震响应敏感, 并在地震监测预测应用中取得一定成果(车用太 等,1998)。近一二十年来,国内外学者将高光谱痕量气体探测技术引进到了地震监测领域,常见的异常分析方法包含背景差值法、异常指数法、标准差法、时间序列分析法等,研究表明利用现有的卫星高光谱数据可以监测到地震前后大气中的CO2、CH4、CO、O3、水汽、气溶胶等异常。中国学者利用AIRS卫星数据研究了2004年12月26日苏门答腊8.9级(孙玉涛 等,2014)、2010年下加利福尼亚7.2级(崔月菊 等,2011)、2008年汶川8.0级(Cui 等,2017)、2013年芦山Ms7.0级(王杰 等,2013)地震前后的CO和CH4总量变化,发现地震前后大气中的CO、CH4浓度变化与地震和构造活动有一定的关系。基于卫星遥感的气溶胶监测也发现地震前后气溶胶光学厚度(AOD)有一定的异常,如2008年5月汶川Ms8.0地震7天前AOD增加(Qin 等,2014)。利用高光谱遥感气体地球化学异常,有助于分析其与地震发生的关系,有助于实现孕震过程的反演与理解。

    (3.2) 遥感技术在地震构造调查方面的应用研究

活动构造是大地震孕育发生的重要条件,也是地震危险区划分的重要依据,是震害防御领域的重要工作内容。目前,在活动构造调查研究方面多源遥感影像数据应用越来越多,《活动断层探察遥感调查》技术标准也已经发布(张景发 等,2017),一方面有效提高了工作效率,同时降低了工作成本。“兵马未动,遥感先行”,已经成为活动构造探查工作的基本模式。

         3.2.1. 区域活动构造格局分析

作为活动构造研究的第一步,需要了解大范围活动构造格局。中高分辨率的MSS、TM、ETM、HJ-1-A/B等卫星遥感信息成为获取大范围活动构造基本状态的重要手段。实际上早在20世纪70年代,许多学者利用卫星影像对亚洲地区的活动构造格局进行了宏观的研究,取得了青藏高原和中亚地区活动构造分布状态,并在此基础上,提出了青藏高原形成演化的构造模式,引领了一个时代世界第三极的研究(Molnar和Tapponnier, 1975, 1977)。中国学者利用陆地资源卫星图像分析揭示了地表广泛分布的不同方向的线性破裂网络(丁国瑜和李永善,1979),解译完成的《中国活动构造典型卫星影像集》和《中国卫星影象地震构造判读图》,给出了中国大陆活动构造分布特征,在此基础上进一步开展的《IGPB-206中国活断层图集》和《中国大陆岩石圈动力学图集》,迄今仍然是地震区划、重点监视防御区圈定和震害防御工作的重要基础性材料。

         3.2.2. 活动构造位移量和精细结构获取

在区域活动构造研究基础上,利用高分辨率遥感影像,包括SPOT、GF-1/2、ZY-3等卫星,以及相应的航空遥感信息,可以对活动断层及其控制的构造微地貌进行全面深入的定量研究,分析断层几何形态、结构特征与构造变形,进而明确断层运动学特征及力学机制。例如,丁国瑜院士(1982)利用MSS卫星影像定量解译了宁夏海原断裂带的水系变形现象;申旭辉等人(2000)在1995年利用区域TM影像资料,发现了四川安宁河东侧凉山活动构造带大规模水平位移证据,修正了鲜水河—安宁河—则木河—小江断裂带终端位移亏损问题;洪顺英等人(2006)研究了阿尔泰山脉蒙古一侧的地震构造,定量给出了活动位移,并初步圈定了历史地震破裂带,为全面认识俄中蒙边界地区阿尔泰地震带的活动提供了有力证据。近年来,随着中国航空航天技术的发展,国产卫星数据也逐渐应用于活动断层定量调查(鲁恒新 等,2017)。此外,无人机或LiDAR技术在活动断层精细结构与定量调查方面也得到广泛应用(Jiao 等,2016Jiang 等,2017),可以提升对活动断裂带精细几何结构与大震复发特征的认识。

         3.2.3. 地震破裂带特征提取

由于大地震造成明显的地表破裂带,应用遥感资料研究地表破裂带的展布、结构甚至位错量的工作也得到普遍应用。如1997年11月8日发生的西藏自治区玛尼(简称西藏玛尼)地震,由于玛尼地区荒无人烟,交通不便,对地震断层的研究十分不便,中巴地球资源卫星CBERS-1图像在该地区的影纹清晰,清楚显示线性形迹长达220 km,清晰分辨出线性刻蚀痕迹,地震专家判释该断层至今活动十分强烈。图像表明,在玛单盖茶卡湖和朝阳湖一带,流经断层带的一系列河流发生同步左向扭曲,玛尼发震断层所造成地表破裂带长110 km,最大左旋位移量达6—7 m。CBERS-1资料对发震条件信息提取和对西藏玛尼地震震源确定、灾情评估提供了科学依据。

         3.2.4. 其他方面

利用雷达数据具有一定穿透能力的特点探测隐伏活动构造的实验研究,在郑州、南京、银川、海口等城市活断层探测工程中发挥积极的作用(洪顺英 等,2007)。基于DEM和多光谱遥感影像丰富的影像特征,可以给出相关区域3维地震构造特征(付碧宏 等,2008),成为系统分析活动构造特征的有效工具。

    (3.3) 卫星遥感技术应用于灾害快速评估

1975年海城地震和1976年唐山地震期间,中国采用遥感飞机对灾区进行了专门飞行任务,为中央救灾决策提供了重要支持,在此后的地震应急实践中,应急遥感技术也迅速发展起来(陈鑫连,1995)。空间对地观测方法的出现和发展为地震调查提供了一种经济、快速、安全的评估方式。在震前应急准备(预案)阶段,遥感图像为风险评估提供参考;在灾后应急阶段,遥感图像用于灾害监测和震害快速评估;在损失评估和重建阶段,遥感图像辅助震害调查、烈度评估等工作。

         3.3.1. 应急遥感信息提取技术研发

应急遥感技术的发展,总体上基于灾害引致的地物状态改变以及由此带来的地物影像特征变化,遥感震害信息提取主要是应用遥感影像的变化检测和图像分类技术。变化检测技术通过对震前和震后遥感影像进行对比分析,检测影像地震前和遭受地震破坏后出现的地物变化及其分布,以此获取地震灾害信息。根据变化检测的信息层次,变化检测方法可分为基于像元的变化检测(张景发 等,2002)、基于特征的变化检测(李成龙和张景发,2013)和基于目标的变化检测(龚丽霞 等,2013)。图像分类技术则直接利用震后遥感影像通过地物分类提取震害信息。传统的分类大多基于像元,利用像元的光谱特征和纹理特征识别震害信息。随着高分辨率遥感影像资料急剧增长,目标呈现更多的对象特征,现在大多采用面向对象的分类方法,或者结合像元与对象来识别震害信息(Gong 等,2016)。

         3.3.2. 应急遥感应急系统研制

“九五”以来,通过“利用卫星遥感技术建立地震灾害预估、评估及救灾决策的运行系统”、“基于遥感图像的震害快速评估技术研究”等项目实施,以卫星遥感信息为主,辅以航空遥感和实地调查资料,首次将卫星遥感技术应用于震害研究中,并初步建立了地震灾情快速获取与损失评估遥感以及通讯技术系统,开发了“地震遥感信息分析与处理系统”并在2003年新疆维吾尔自治区伽师(简称新疆伽师)地震和2005年江西九江地震震灾应急工作中得到检验(王晓青 等,2003)。2011年启动实施的高分地震遥感监测与应急应用示范工程,研制了应急遥感示范系统,该系统集成了建筑物、道路交通、地震地质灾害3个主要地震灾害识别功能,能全面地反映灾区宏观灾情,在地震应急遥感业务工作具有重要作用(窦爱霞 等,2013)。

         3.3.3. 应急遥感应用实践

应急遥感技术业务化应用,始于2003年—2004年新疆巴楚—伽师地震震群活动,利用大量的航空与卫星遥感资料,快速处理获取了地震造成的建筑物震害与地质灾害的遥感图像特征,并根据以往震害遥感影像统计经验与本次地震震害遥感特征,提出了遥感震害分级分类标准和地震烈度划分标准,进而得到基于震害遥感影像的伽师地震等震线图,为应急救援决策提供了重要依据(王晓青 等,2003)。2008年汶川特大地震后,在缺乏应有的地面震害调查资料情况下,通过遥感手段获取灾区震后遥感影像,图像快速处理和震害信息提取,对地震灾害及其损失进行快速评估,将弥补地面调查的不足,为地震应急指挥和救援决策提供重要依据(王晓青 等,2009)。

2008年汶川地震之后,遥感灾害评估工作加快步伐,提高了其实用性,在2013年芦山地震和2017年九寨沟等地震中都有良好表现,震后24 h内,就利用高分辨率遥感数据对地震损毁程度进行了初步评估,用于应急救援决策支持,并向民众发布大致灾情。

4、地球物理场探测卫星计划与电磁监测试验卫星     (4.1) 地球物理场探测卫星计划概况

早在1977年,中国科学家就提出卫星地球物理探测将在地震研究中发挥重要作用(中国科学院地质研究所,1977)。20世纪末至21世纪初,国际学术界围绕地震能否预报展开了一场大讨论(Geller 等,1997Wyss,1997)。中国学者在系统回顾和总结中国地震监测预测发展历程和取得的主要成果基础上,倾向认为地震作为一个自然科学问题,应该可以预测(陈运泰,2009)。而空间对地观测技术以其大范围、高动态、全天时、全天候、不受地面条件限制等技术优势,成为现阶段提高地震监测能力,推进地震预测科学探索的重要途径。在充分论证基础上,提出了中国地震立体观测体系发展的总体框架和发展思路(申旭辉 等,2007),地震立体观测系统概念框图如图3所示。

图 3 地震立体观测系统概念框图 Figure 3 Conceptual diagram of seismic stereo observation system

地震预测的主要困难有:地震震源的不可入性、破坏性地震发生的小概率性和地震监测预测技术手段的局限性。空间对地观测能够显著弥补现有地震观测能力不足,从而为推进防震减灾能力建设提供十分重要的信息保障。

根据统计,中国平均每3年发生2次7级以上地震,每年发生4次6级以上地震,而全球平均每年发生7级以上地震17—18次,6级以上地震150次,通过卫星观测,能够获取全球震例,从而显著弥补破坏性地震样本不足的问题,为科学家研究和发展地震预测方法提供支持。与此同时,地面监测系统只能给出孤立点的地球物理场信息,地震监测预测研究需要完整的3维地球物理场动态信息可以通过地球物理场卫星技术获取,从而有效提高地震监测能力,推进地震预测科学探索。

2015年,国务院批准《国家民用空间基础设施中长期发展规划(2015—2025)》,正式明确发展国家地球物理场探测卫星计划。总体上,地球物理场探测卫星主要包含电磁卫星和重力卫星两大类。其中,电磁监测卫星以获取地球基本磁场、电磁场和电磁波以及电离层等离子体结构和原位参数为主要任务;重力场测量卫星以获取地球重力场和重力梯度场为主要任务目标。

通过获取地球磁场、电场、重力场、电离层模型及其变化,可直接用于支撑导航、惯性制导和通信环境管理。与此同时,通过3维地球物理场环境及其变化状态获取,将为有利推进空间地球物理和空间物理学发展,发展地球多圈层联合探测技术,支撑地球系统科学发展和地震监测预测科学突破奠定重要的基础。

    (4.2) 电磁监测试验卫星

作为中国地球物理场探测卫星计划首发星,电磁监测试验卫星(张衡一号)于2003年初开始规划预研、2013年启动工程研制工作,于2018年2月2日顺利进入预定轨道。

张衡一号卫星的主要科学目标是:获取全球电磁场、等离子体、高能粒子观测数据,进行卫星在轨期间全球7级以上、中国6级以上地震卫星电磁信息分析研究,总结地震电离层扰动特征;对中国全境实施准实时观测,开展试验性地震短临监测,研究大震短临预测新方法;开展电磁卫星技术试验;为相关领域应用与国际合作提供电磁数据信息服务。

卫星采用三轴稳定太阳同步极轨圆轨道,轨道高度507 km,倾角97.4°,降交点地方时14:00,卫星重访周期设定为5 d,卫星设计寿命5 a。按照设计,张衡一号卫星采用CAST2000通用小卫星平台,对地三轴稳定,星上搭载3类8种科学载荷,包括:电磁场测量载荷、高精度磁强计、感应式磁力仪和电场仪;电离层原位测量载荷:朗缪尔探针、等离子体分析仪和高能粒子探测器;电离层结构探测载荷:GNSS掩星接收机和三频信标发射机(卫星构型如图4所示)。

图 4 中国首颗电磁监测试验卫星构型 Figure 4 China’s first configuration of electromagnetic monitoring test satellite

张衡一号卫星在轨期间,已产出全球地磁场和大气电场分频段波形和频谱数据,电离层原位电子和离子的密度、温度、离子漂移速度及高能粒子沉降数据,电离层等离子体2维和3维结构探测数据,提供地震科学、固体地球科学和空间科学研究使用,同时服务于地震监测、空间天气监测预警以及电波环境管理等任务。

5、思考与展望

到目前为止,空间技术应用到防震减灾事业发展和相关科学研究已经走过了40余年的历程。空间技术应用已经渗透到防震减灾工作各个工作环节,为地震监测预测、震害防御及应急救援科技水平提升和业务能力建设提供了重要支撑。从防震减灾业务能力建设需求出发,考虑到地震遥感应用现状和空间对地观测技术发展,未来地震遥感应用将重点在如下方向推进:

(1) 以地球物理场探测卫星计划发展为基础,发展综合性、天空地一体化、有效协同的地震立体观测体系。当前的重点是电磁卫星和重力卫星等地球物理场探测卫星系列的发展,在此基础上,构建融合光学和微波遥感平台于一体的虚拟化业务卫星星座,并推进临近空间艇载地球物理综合观测平台和地面比测校验系统发展。

(2) 加快推进遥感应用业务化水平。目前,应急遥感应用已经具备业务化能力,卫星红外遥感已经成为日常地震监测预测的重要手段,基于GNSS与InSAR相结合的地壳形变立体监测网络正在逐步形成,电磁监测试验卫星在轨运行。多源空间信息应用与业务系统的密切结合是地震遥感应用的一个重要发展方向。

(3) 积极拓展空间信息综合分析与大数据技术应用,发展有效的融合技术与分析方法,提高数据处理能力,发展强背景和复杂环境下地震微弱信号提取技术,加强新技术研发和应用。

(4) 在传统的板块和块体相互作用模型基础上,充分利用立体观测系统特点和优势,加强地球多圈层、多物理场耦合与相互作用特征与机制研究,推进地震科学理论取得突破。

参考文献
[1] 蔡永恩, 殷有泉, 王仁. 地震前断层蠕动与地温异常的探讨[J]. 地震学报, 1987, 9 (2) : 167 –175. Cai Y N, Yin Y Q and Wang R. A study on preseismic fault creep and ground temprature anomaly[J]. Acta Seismologica Sinica, 1987, 9 (2) : 167 –175.
[2] Cambiotti G, Bordoni A, Sabadini R and Colli L. GRACE gravity data help constraining seismic models of the 2004 Sumatran earthquake[J]. Journal of Geophysical Research, 2011, 116 (B10) : B10403 . DOI: 10.1029/2010JB007848
[3] 车用太, 王基华, 林元武, 朱清钟, 张培仁. 张北-尚义地震前的地下流体异常及其跟踪预报[J]. 地震地质, 1998, 20 (2) : 99 –104. Che Y T, Wang J H, Lin Y W, Lin Y W, Zhu Q Z and Zhang P R. The anomalies of subsurface fluids before the Zhangbei Shangyi earthquake and trace prediction of the earthquake[J]. Seismology and Geology, 1998, 20 (2) : 99 –104.
[4] 陈鑫连. 1995. 地震灾害的航空遥感信息快速评估与救灾决策. 北京: 科学出版社 Chen X L. 1995. Quick Evaluation and Decision-making for Earthquake Disaster by Airborne Remote Sensing Information. Beijing: Science Press
[5] 陈运泰. 地震预测: 回顾与展望[J]. 中国科学D辑: 地球科学, 2009, 39 (12) : 1633 –1658. Chen Y T. Earthquake prediction: review and prospect[J]. Science China Series D: Earth Science, 2009, 39 (12) : 1633 –1658.
[6] Cui Y, Ouzounov D, Hatzopoulos N, Sun K, Zou Z and Du J. Satellite observation of CH4 and CO anomalies associated with the Wenchuan MS 8.0 and Lushan MS 7.0 earthquakes in China [J]. Chemical Geology, 2017, 469 : 185 –191. DOI: 10.1016/j.chemgeo.2017.06.028
[7] 崔月菊, 杜建国, 周晓成, 陈志, 李营, 刘雷, 谢超, 张炜斌. 墨西哥下加利福尼亚Mw7.2地震前后CO遥感地球化学异常[J]. 矿物岩石地球化学通报, 2011, 30 (4) : 458 –464. Cui Y J, Du J G, Zhou X C, Chen Z, Li Y, Liu L, Xie C and Zhang W B. Geochemical anomaly of CO remote sensing associated with Baja California Mw 7.2 Earthquake in Mexico[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2011, 30 (4) : 458 –464. DOI: 10.3969/j.issn.1007–2802.2011.04.014
[8] 邓明德, 耿乃光, 崔承禹, 支毅乔, 樊正芳, 籍全权. 岩石应力状态改变引起岩石热状态改变的研究[J]. 中国地震, 1997, 13 (2) : 179 –185. Deng M D, Geng N G, Cui C Y, Zhi Y Q, Fan Z F and Ji Q Q. The study on the variation of thermal state of rocks caused by the variation of stress state of rocks[J]. Earthquake Resea Rch in China, 1997, 13 (2) : 179 –185.
[9] 丁国瑜, 李永善. 中国地震活动与地壳现代破裂网络[J]. 地质学报, 1979, 53 (1) : 22 –34. Ding G Y and Li Y S. Seismicity and the recent fracturing pattern of the earth crust in China[J]. Acta Geological Sinica, 1979, 53 (1) : 22 –34.
[10] 窦爱霞, 丁香, 王龙, 王晓青, 单新建, 张景发, 马海建. 高分光学遥感地震应急应用原型系统[J]. 地震, 2013, 33 (2) : 11 –18. Dou A X, Ding X, Wang L, Wang X Q, Shan X J, Zhang J F and Ma H J. A prototype system for emergency earthquake damage analysis and processing based on optical high resolution images[J]. Earthquake, 2013, 33 (2) : 11 –18. DOI: 10.3969/j.issn.1000–3274.2013.02.002
[11] Ferretti A, Prati C and Rocca F. Permanent scatterers in SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39 : 8 –20. DOI: 10.1109/36.898661
[12] Freund F T. Pre-earthquake signals-Part II: flow of battery currents in the crust[J]. Natural Hazards and Earth System sciences, 2007, 7 (5) : 543 –548. DOI: 10.5194/nhess–7–543–2007
[13] 付碧宏, 二宫芳树, 董彦芳, 时丕龙, 张之武. 三维卫星遥感图像生成技术及其在第四纪构造地貌研究中的应用[J]. 第四纪研究, 2008, 28 (2) : 189 –196. Fu B H, Yoshiki N, Dong Y F, Shi P L and Zhang Z W. Generation of 3-Dimensional perspective satellite imagery and its application on quaternary geomorphological analysis[J]. Quaternary Sciences, 2008, 28 (2) : 189 –196. DOI: 10.3321/j.issn:1001–7410.2008.02.001
[14] Garcia R F, Bruinsma S, Lognonné P, Doornbos E and Cachoux F. GOCE: the first seismometer in orbit around the Earth[J]. Geophysical Research Letters, 2013, 40 (5) : 1015 –1020. DOI: 10.1002/grl.50205
[15] Geller R J, Jackson D D, Kagan Y Y and Mulargia F. Earthquakes cannot be predicted[J]. Science, 1997, 275 (5306) : 1616 . DOI: 10.1126/science.275.5306.1616
[16] 龚丽霞, 李强, 张景发, 曾琪明, 刘明众, 李成龙. 面向对象的房屋震害变化检测方法[J]. 地震, 2013, 33 (2) : 109 –114. Gong L X, Li Q, Zhang J F, Zeng Q M, Liu M Z and Li C L. Object-oriented detection of earthquake building damages[J]. Earthquake, 2013, 33 (2) : 109 –114. DOI: 10.3969/j.issn.1000–3274.2013.02.014
[17] Gong L X, Wang C, Wu F, Zhang J F, Zhang H and Li Q. Earthquake-induced building damage detection with post-event Sub-Meter VHR terraSAR-X staring spotlight imagery[J]. Remote Sensing, 2016, 8 (11) : 887 . DOI: 10.3390/rs8110887
[18] Gorny V I, Salman A G, Tronin A A and Shilin B B. Terrestrial outgoing infrared radiation as an indicator of seismic activity[J]. Proceeding of the Academy of Sciences of the USSR, 1988, 30 (1) : 67 –69.
[19] 郭子祺. 航天遥感技术进行短临地震预报评述[J]. 遥感信息, 1998, 2 : 13 –15. Guo Z Q. Space remote sensing technology for short - impending earthquake prediction[J]. Remote Sensing Information, 1998, 2 : 13 –15.
[20] Han S C, Shum C K, Bevis M, Ji C and Kuo C Y. Crustal dilatation observed by GRACE after the 2004 Sumatra-Andaman Earthquake[J]. Science, 2006, 313 (5787) : 658 –662. DOI: 10.1126/science.1128661
[21] Han S C, Sauber J, and Riva R. Contribution of satellite gravimetry to understanding seismic source processes of the 2011 Tohoku-Oki earthquake[J]. Geophysical Research Letters, 2011, 38 : 24312 . DOI: 10.1029/2011GL049975
[22] Hao Y Q, Xiao Z and Zhang D H. Multi-instrument observation on co-seismic ionospheric effects after great Tohoku earthquake[J]. Journal of Geophysical Research, 2012, 117 (A2) : A02305 . DOI: 10.1029/2011JA017036
[23] He Y, Yang D, Qian J and Parrot M. Response of the ionospheric electron density to different types of seismic events[J]. Natural Hazards and Earth System Sciences, 2011, 11 (8) : 2173 –2180. DOI: 10.5194/nhess–11–2173–2011
[24] 洪顺英, 申旭辉, 赖木收, 陈正位, 荆凤, 欧阳新艳, RogoZhin E A. 阿尔泰山东缘主要活动断裂影像特征分析[J]. 地震地质, 2006, 28 (1) : 119 –128. Hong S Y, Shen X H, Lai M S, Chen Z W, Jing F, Ouyang X Y and RogoZhin E A. Analyses of image feature of major active faults along the eastern margin of Altai mountains[J]. Seismology and Geology, 2006, 28 (1) : 119 –128. DOI: 10.3969/j.issn.0253–4967.2006.01.012
[25] 洪顺英, 申旭辉, 单新建, 刘智荣, 戴娅琼, 荆凤. 基于升降轨ASAR的于田Ms 7.3级地震同震形变场信息提取与分析 [J]. 国土资源遥感, 2010, 22 (4) : 98 –102. Hong S Y, Shen X H, Shan X J, Liu Z R, Dai Y Q and Jing F. The calculation and analysis of the co-seismic deformation field of Yutian Ms 7.3 earthquake basing on the ascending and descending orbit ASAR data [J]. Remote Sensing For Land and Resources, 2010, 22 (4) : 98 –102. DOI: 10.6046/gtzyyg.2010.04.20
[26] 洪顺英, 张红英, 申旭辉, 胡友健, 张丽娜. 南京市活动断层卫星遥感图像的研究[J]. 遥感信息, 2007 (4) : 45 –50. Hong S Y, Zhang H Y, Shen X H, Hu Y J and Zhang L N. The research on remote sensing images of Nanjing active falts[J]. Remote Sensing Information, 2007 (4) : 45 –50. DOI: 10.3969/j.issn.1000–3177.2007.04.010
[27] 黄培华, 金凤英, 丁宝田, 韩健. 1982. 地震地质学基础. 北京: 地震出版社 Huang P H, Jin F Y, Ding B T and Han J. 1982. Foundation of Seismic Geology. Beijing: Seismological Press
[28] 国家地震局地质研究所. 1985. 遥感地震地质文集. 北京: 地震出版社, 55-125 Institute of Geology, China Earthquake Administration. 1985. Remote Sensing Seismic Geology Anthology. Beijing: Seismological Press, 55-125
[29] 中国科学院地质研究所. 1977. 中国地震地质概论. 北京: 科学出版社, 18-22 Institute of Geology, Chinese Academy of Geological Sciences. 1977. Introduction to China’s Seismic Geology. Beijing: Science Press, 18-22
[30] Jiang W L, Han Z J, Guo P, Zhang J F, Jiao Q S and Kang S. Slip rate and recurrence intervals of the east Lenglongling fault constrained by morphotectonics: tectonic implications for the northeastern Tibetan Plateau[J]. Lithosphere, 2017, 9 (3) : 417 –430. DOI: 10.1130/L597.1
[31] Jiao Q, Jiang W L, Zhang J F, Jiang H B, Luo Y and Wang X. Identification of paleoearthquakes based on geomorphological evidence and their tectonic implications for the southern part of the active Anqiu–Juxian fault, eastern China[J]. Journal of Asian Earth Sciences, 2016, 132 : 1 –8. DOI: 10.1016/j.jseaes.2016.10.012
[32] 康春丽, 陈正位, 陈立泽, 田勤俭, 刘德富. 昆仑山口西8.1级地震的卫星热红外前兆特征分析[J]. 西北地震学报, 2003, 25 (1) : 12 –15. Kang C L, Chen Z W, Chen L Z, Tian Q J and Liu D F. Analysis on the satellite infrared anomaly feature before west to Kunlun mountain pass M8.1 earthquake [J]. Northwestern Seismological Journal, 2003, 25 (1) : 12 –15. DOI: 10.3969/j.issn.1000–0844.2003.01.003
[33] Le H J, Liu J Y and Liu L. A statistical analysis of ionospheric anomalies before 736 M6.0+ earthquakes during 2002-2010 [J]. Journal of Geophysical Research, 2011, 116 (A2) : A02303 . DOI: 10.1029/2010JA015781
[34] 李成龙, 张景发. 基于主成分分析的遥感震害变化检测方法与应用[J]. 地震, 2013, 33 (2) : 103 –108. Li C L and Zhang J F. Method and application of earthquake damage change detection from remote sensing images based on principal component analysis[J]. Earthquake, 2013, 33 (2) : 103 –108. DOI: 10.3969/j.issn.1000–3274.2013.02.013
[35] Li M and Parrot M. Statistical analysis of an ionospheric parameter as a base for earthquake prediction[J]. Journal of Geophysical Research, 2013, 118 (6) : 3731 –3739. DOI: 10.1002/jgra.50313
[36] 刘德富, 康春丽. 苏门答腊岛8.7级大地震前的辐射异常现象[J]. 国际地震动态, 2005, 1 : 36 –39. Liu D F and Kang C L. Radiation anomalism before Sumatra Island M8.7 earthquake[J]. Recent Dvelopment in World Seismology, 2005, 1 : 36 –39. DOI: 10.3969/j.issn.0253–4975.2005.05.007
[37] 刘静, 万卫星. 中国6.0级以上地震临震电离层扰动时空分布特征研究[J]. 地球物理学报, 2014, 57 (7) : 2181 –2189. Liu J and Wan W X. Spatial-temporal distribution of the ionospheric perturbations prior to Ms≥6.0 earthquakes in China main land [J]. Chinese Journal of Geophysics, 2014, 57 (7) : 2181 –2189. DOI: 10.6038/CJG20140713
[38] Liu J, Zhang X M, Novikov V and Shen X H. Variations of ionospheric plasma at different altitudes before the 2005 Sumatra Indonesia Ms 7.2 earthquake [J]. Journal of Geophysical Research, 2016, 121 (9) : 9179 –9187. DOI: 10.1002/2016JA022758
[39] 鲁恒新, 徐岳仁, 陈立泽, 张学民, 申旭辉, 崔静. 国产GF-2和ZY-3卫星数据在活动断裂定量研究中的应用——以哈思山南麓断裂带为例[J]. 地震, 2017, 37 (1) : 121 –133. Lu H X, Xu Y R, Chen L Z, Zhang X M, Shen X H and Cui J. Application of GF-2 and ZY-3 satellite data in quantitative research of active fault: a case study of the south Hasi mountain fault zone[J]. Earthquake, 2017, 37 (1) : 121 –133. DOI: 10.3969/j.issn.1000-3274.2017.01.013
[40] Massonnet D, Rossi M, Carmona C, Adragna F, Peltzer G, Feigl K and Rabaute T. The displacement field of the Landers earthquake mapped by radar interferometry[J]. Nature, 1993, 364 (6433) : 138 –142. DOI: 10.1038/364138a0
[41] Molnar P and Tapponnier P. Cenozoic Tectonics of Asia: effects of a Continental Collision: features of recent continental tectonics in Asia can be interpreted as results of the India-Eurasia collision[J]. Science, 1975, 189 (4201) : 419 –426. DOI: 10.1126/science.189.4201.419
[42] Molnar P and Tapponnier P. The collision between India and Eurasia[J]. Scientific American, 1977, 236 (4) : 30 –41. DOI: 10.1038/scientificamerican0477–30
[43] Panet I, Mikhailov V, Diament M, Pollitz F, King G, De Viron O, Holschneider M, Biancale R and Lemoine J M. Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity[J]. Geophysical Journal International, 2007, 171 (1) : 177 –190. DOI: 10.1111/j.1365–246X.2007.03525.x
[44] Pergola N, Aliano C, Coviello I, Filizzola C, Genzano N, Lacava T, Lisi M, Mazzeo G and Tramutoli V. Using RST approach and EOS-MODIS radiances for monitoring seismically active regions: a study on the 6 April 2009 Abruzzo earthquake[J]. Natural Hazards and Earth System Sciences, 2010, 10 (2) : 239 –249. DOI: 10.5194/nhess–10–239–2010
[45] Pulinets S A and Boyarchuk K A. 2004. Ionospheric Precursors of Earthquakes. Berlin: Springer: 1–289
[46] Pulinets S A. Natural radioactivity, earthquakes, and the ionosphere[J]. EOS, Transactions American Geophysical Union, 2007, 88 (20) : 217 –218. DOI: 10.1029/2007EO200001
[47] Qiang Z J, Kong L C, Guo M H, San G P, Lin C G and Zhao Y. Laboratory research on mechanism of satellite infrared anomaly[J]. Chinese Science Bulletin, 1995, 40 (16) : 1403 –1404. DOI: 10.1360/sb1995–40–16–1403
[48] Qin K, Wu L X, Zheng S, Bai Y and Lv X. Is there an abnormal enhancement of atmospheric aerosol before the 2008 Wenchuan earthquake?[J]. Advances in Space Research, 2014, 54 (6) : 1029 –1034. DOI: 10.1016/j.asr.2014.04.025
[49] 屈春燕, 张桂芳, 单新建, 张国宏, 刘云华, 宋小刚. 2010年青海玉树地震同震—震后形变场特征及演化过程[J]. 地球物理学报, 2013, 56 (7) : 2280 –2291. Qu C Y, Zhang G F, Shan X J, Zhang G H, Liu Y H, and Song X G. Coseismic and postseismic deformation fields of the 2010 Yushu, Qinghai MS7.1 earthquake and their evolution processes [J]. Chinese Journal of Geophysics, 2013, 56 (7) : 2280 –2291. DOI: 10.6038/cjg20130715
[50] Saradjian M R and Akhoondzadeh M. Thermal anomalies detection before strong earthquakes (M>6.0) using interquartile, wavelet and Kalman filter methods [J]. Natural Hazards and Earth System Sciences, 2011, 11 (4) : 1099 –1108. DOI: 10.5194/nhess–11–1099–2011
[51] Saraf A K and Choudhury S. Thermal remote sensing technique in the study of pre-earthquake thermal anomalies[J]. The Journal of Indian Geophysical Union, 2005, 9 (3) : 197 –207.
[52] 单新建, 柳稼航, 马超. 2001年昆仑山口西8.1级地震同震形变场特征的初步分析[J]. 地震学报, 2004, 26 (5) : 474 –480. Shan X J, Liu J H and Ma C. Preliminary analysis on characteristics of coseismic deformation associated with Ms=8.1 western Kunlunshan pass earthquake in 2001 [J]. Acta Seismologica Sinica, 2004, 26 (5) : 474 –480. DOI: 10.3321/j.issn:0253–3782.2004.05.003
[53] 单新建, 马瑾, 王长林, 柳稼航. 利用差分干涉雷达测量技术(D-InSAR)提取同震形变场[J]. 地震学报, 2002, 24 (4) : 413 –420. Shan X J, Ma J, Wang C L and Liu J H. Extracting coseismic deformation of the 1997 Mani earthquake with differen-tial interferometric SAR[J]. Acta Seismologica Sinica, 2002, 24 (4) : 413 –420. DOI: 10.3321/j.issn:0253–3782.2002.04.009
[54] 申旭辉, 陈正位, 许任德, 张宁. 凉山活动构造带晚新生代变形特征与位移规模[J]. 地震地质, 2000, 22 (3) : 232 –238. Shen X H, Chen Z W, Xu R D and Zhang N. Deformation characteristics and displacement amount of the Liangshan active fault zone in late Cenozoic era[J]. Seismology and Geology, 2000, 22 (3) : 232 –238. DOI: 10.3969/j.issn.0253–4967.2000.03.004
[55] 申旭辉, 吴云, 单新建. 地震遥感应用趋势与中国地震卫星发展框架[J]. 国际地震动态, 2007 (8) : 38 –45. Shen X H, Wu Y and Shan X J. Remote sensing application in earthquake science and general proposal for earthquake satellite project in China[J]. Recent Developments in World Seismology, 2007 (8) : 38 –45. DOI: 10.3969/j.issn.0253-4975.2007.08.011
[56] Shen X H, Zhang X M, Wang L W, Chen H R, Wu Y, Yuan S G, Shen J F, Zhao S F, Quan J D and Ding J H. The earthquake-related Disturbances in ionosphere and project of the First China seismo-electromagnetic satellite[J]. Earthquake Science, 2011, 24 (6) : 639 –650. DOI: 10.1007/s11589–011–0824–0
[57] Shen X H, Zhima Z, Zhao S F, Qian G, Ye Q and Ruzhin Y. VLF radio wave anomalies associated with the 2010 Ms 7.1 Yushu earthquake [J]. Advances in Space Research, 2017, 59 (10) : 2636 –2644. DOI: 10.1016/j.asr.2017.02.040
[58] 孙建宝, 梁芳, 沈正康, 徐锡伟. 汶川Ms 8.0地震InSAR形变观测及初步分析 [J]. 地震地质, 2008, 30 (3) : 769 –759. Sun J B, Liang F, Shen Z K and Xu X W. InSAR deformation observation and preliminary analysis of the Ms 8 Wenchuan earthquake [J]. Seismology and Geology, 2008, 30 (3) : 769 –759.
[59] Sun W K and Okubo S. Co-seismic deformations detectable by satellite gravity missions: a case study of Alaska (1964, 2002) and Hokkaido (2003) earthquakes in the spectral domain[J]. Journal of Geophysical Research, 2004, 109 (B4) : B04405 . DOI: 10.1029/2003JB002554
[60] 孙玉涛, 崔月菊, 刘永梅, 杜建国, 张炜斌, 张冠亚. 苏门答腊2004、2005年两次大地震前后CO和O3遥感信息 [J]. 遥感信息, 2014, 29 (2) : 49 –55. Sun Y T, Cui Y J, Liu Y M, Du J G, Zhang W B and Zhang G Y. Remote sensing anomalies of CO and O3 related to two giant Sumatra earthquakes occurred in 2004 and 2005 [J]. Remote Sensing Information, 2014, 29 (2) : 49 –55. DOI: 10.3969/j.issn.1000–3177.2014.02.009
[61] Tronin A A, Hayakawa M and Molchanov O A. Thermal IR satellite data application for earthquake research in Japan and China[J]. Journal of Geodynamics, 2002, 33 (4/5) : 519 –534. DOI: 10.1016/S0264–3707(02)00013–3
[62] 王超, 刘智, 张红, 单新建. 张北-尚义地震同震形变场雷达差分干涉测量[J]. 科学通报, 2000, 46 (6) : 514 –518. Wang C, Liu Z, Zhang H and Shan X J. The coseismic displacement field of the Zhangbei-Shangyi earthquake mapped by differ- ential radar interferometry[J]. Chinese Science Bulletin, 2000, 46 (6) : 514 –518. DOI: 10.3321/j.issn:0023–074X.2000.23.017
[63] 王纯莹, 朱子豪, 刘英毓. 利用时空内插及数据正规化辅助方法监测地震前热异常[J]. 国土资源遥感, 2005 (4) : 11 –15. Wang C Y, Zhu Z H and Liu Y Y. The application of spatial –temporal interpolation and data normalization methods to detecting the pre- earthquake abnormal temperature increasing phenomenon[J]. Remote Sensing for Land and Resources, 2005 (4) : 11 –15. DOI: 10.3969/j.issn.1001–070X.2005.04.003
[64] 王杰, 张雄, 潘黎黎, 曾佐勋. 芦山地震(Ms7.0)前甲烷释放与大气增温异常 [J]. 地学前缘, 2013, 20 (6) : 29 –35. Wang J, Zhang X, Pan L L and Zeng Z X. Anomalies of temperature increase and methane release before Lushan earthquake (Ms7.0) [J]. Earth Science Frontiers, 2013, 20 (6) : 29 –35.
[65] 王晓青, 窦爱霞, 孙国清, 丁香, 王龙, 袁小祥. 基于综合震害指数的玉树地震烈度遥感评估研究[J]. 地震, 2013, 33 (2) : 1 –10. Wang X Q, Dou A X, Sun G Q, Ding X, Wang L and Yuan X X. Intensity assessment of the 2010 Yushu Ms7.1 earthquake based on synthetic seismic damage index [J]. Earthquake, 2013, 33 (2) : 1 –10. DOI: 10.3969/j.issn.1000–3274.2013.02.001
[66] 王晓青, 王龙, 章熙海, 丁香, 邱虎, 董彦芳, 窦爱霞, 张飞宇, 刘在涛. 汶川8.0级地震震害遥感定量化初步研究——以都江堰城区破坏为例[J]. 地震, 2009, 29 (1) : 174 –181. Wang X Q, Wang L, Zhang X H, Ding X, Qiu H, Dong Y F, Dou A X, Zhang F Y and Liu Z T. Primary quantitative study on earthquake damage extracted from remote sensing imagery—a case study of Doujiangyan due to the Wenchuan M8.0 earthquake [J]. Earthquake, 2009, 29 (1) : 174 –181. DOI: 10.3969/j.issn.1000–3274.2009.01.021
[67] 王晓青, 魏成阶, 苗崇刚, 张景发, 单新建, 马庆尊. 震害遥感快速提取研究——以2003年2月24日巴楚—伽师6.8级地震为例[J]. 地学前缘, 2003, 10 (S1) : 285 –291. Wang X Q, Wei C J, Miao C G, Zhang J F, Shan X J and Ma Q Z. The extraction of seismic damage from remote sensing images—a case study of Bachu-Jiashi earthquake with Ms=6.8 occurred on Feb 24, 2003 [J]. Earth Science Frontiers, 2003, 10 (S1) : 285 –291. DOI: 10.3321/j.issn:1005–2321.2003.z1.039
[68] 吴立新. 遥感岩石力学及其新近进展与未来发展[J]. 岩石力学与工程学报, 2001, 20 (2) : 139 –146. Wu L X. Remote sensing rock mechanics and its recent achievements and future development[J]. Chinese Journal of Rock Mechanics and Engineering, 2001, 20 (2) : 139 –146. DOI: 10.3321/j.issn:1000–6915.2001.02.001
[69] Wu L X, Cai C Y, Geng N G and Wang J Z. Remote sensing rock mechanics (RSRM) and associated experimental studies[J]. International Journal of Rock Mechanics and Mining Sciences, 2000, 37 (6) : 879 –888. DOI: 10.1016/S1365–1609(99)00066–0
[70] Wu L X, Qin K and Liu S J. GEOSS-based thermal parameters analysis for earthquake anomaly recognition[J]. Proc IEEE, 2012a, 100 (10) : 2891 –2907. DOI: 10.1109/JPROC.2012.2184789
[71] Wu L X, Qin K, Liu S J, De S A and Cianchini G. 2012b. Importance of lithosphere–coversphere–atmosphere coupling to earthquake anomaly recognition, geoscience and remote sensing symposium (IGARSS), 2012 IEEE International IEEE, 3532-3535 [DOI: 10.1109/IGARSS.2012.6350657]
[72] Wyss M, Aceves R L, Park S K, Geller R J, Jackson D D, Kagan Y Y and Mulargia F. 1997. Cannot earthquakes be predicted? Science, 278(5337): 487–490 [DOI: 10.1126/science.278.5337.487]
[73] 谢广林. 1991. 应用卫星图像分析中国东部断裂活动的几点特征. 地质学报, (2): 116–126 Xie G L. 1991. Analysis of some characteristics of the fault activities in eastern china by satellite images. Acta Geological Sinica, (2): 116–126
[74] 徐彤, 胡艳莉, 吴健, 李春斌, 吴振森, 索玉成, 冯健. 中国大陆14次强震前电离层异常统计分析[J]. 电波科学学报, 2012, 27 (3) : 507 –512. Xu T, Hu Y L, Wu J, Li C B, Wu Z S, Suo Y C and Feng J. Statistical analysis of seismo-ionospheric perturbation before 14 Ms≥7.0 strong earthquakes in Chinese subcontinent [J]. Chinese Journal of Radio Science, 2012, 27 (3) : 507 –512. DOI: 10.13443/j.cjors.2012.03.005
[75] 徐秀登, 徐向民, 王煜. 台湾南投7.6级地震前卫星红外异常[J]. 地震学报, 2000, 22 (6) : 666 –669. Xu X D, Xu X M and Wang Y. Satellite infrared anomaly before Nantou Ms=7.6 earthquake in Taiwan, China [J]. Acta Seismologica Sinica, 2000, 22 (6) : 666 –669. DOI: 10.3321/j.issn:0253–3782.2000.06.014
[76] Yao Y B, Chen P, Zhang S, Chen J J, Yan F and Peng W F. Analysis of pre-earthquake ionospheric anomalies before the global M=7.0+ earthquakes in 2010 [J]. Natural Hazards and Earth System Sciences, 2012, 12 (3) : 575 –585. DOI: 10.5194/nhess–12–575–2012
[77] 张景发, 刘钊. InSAR技术在西藏玛尼强震区的应用[J]. 清华大学学报(自然科学版), 2002, 42 (6) : 847 –850. Zhang J F and Liu Z. Application of InSAR technology the strong earthquake in Mani, Tibet[J]. Journal of Tsinghua University (Science and Technology), 2002, 42 (6) : 847 –850. DOI: 10.3321/j.issn:1000–0054.2002.06.037
[78] 张景发, 谢礼立, 陶夏新. 建筑物震害遥感图像的变化检测与震害评估[J]. 自然灾害学报, 2002, 11 (2) : 59 –64. Zhang J F, Xie L L and Tao X X. Change detection of remote sensing image for earthquake-damaged buildings and its application in seismic disaster assessment[J]. Journal of Natural Disasters, 2002, 11 (2) : 59 –64. DOI: 10.3969/j.issn.1004–4574.2002.02.010
[79] 张景发, 徐锡伟, 姜文亮, 申旭辉, 田甜, 田勤俭, 何宏林, 杨晓平, 秦其明, 王鑫, 罗毅, 李永生, 焦其松, 李强. 2017. DB/T 69-2017. 活动断层探察 遥感调查. 北京: 中国地震局 Zhang J F, Xu X W, Jiang W L, Shen X H, Tian T, Tian Q J, He H L, Yang X P, Qin Q M, Wang X, Luo Y, Li Y S, Jiao Q S and Li Q. 2017. DB/T 69-2017. Active fault exploration and remote sensing investigation. China Earthquake Administration.
[80] Zhang X, Shen X, Liu J, Ouyang X, Qian J and Zhao S. Analysis of ionospheric plasma perturbations before Wenchuan earthquake[J]. Natural Hazards and Earth System Sciences, 2009, 9 (4) : 1259 –1266. DOI: 10.5194/nhess–9–1259–2009
[81] Zhang X, Shen X, Zhao S, Yao L, Ouyang X and Qian J. The characteristics of quasistatic electric field perturbations observed by DEMETER satellite before large earthquakes[J]. Journal of Asian Earth Sciences, 2014, 79 : 42 –52. DOI: 10.1016/j.jseaes.2013.08.026
[82] Zhang X, Zeren Z, Parrot M, Battiston R, Qian J and Shen X. ULF/ELF ionospheric electric field and plasma perturbations related to Chile Earthquakes[J]. Advances in Space Research, 2011, 47 (6) : 991 –1000. DOI: 10.1016/j.asr.2010.11.001
[83] 张学民, 丁鉴海, 申旭辉, 王敏, 刘静, 余素荣, 王亚丽, 欧阳新艳. 汶川8级地震前电磁扰动与电磁立体监测体系[J]. 电波科学学报, 2009, 24 (1) : 1 –8. Zhang X M, Ding J H, Shen X H, Wang M, Liu J, Yu S R, Wang Y L and Ouyang X Y. Electromagnetic perturbations before Wenchuan M8 earthquake and stereo electromagnetic observation system[J]. Chinese Journal of Radio Science, 2009, 24 (1) : 1 –8. DOI: 10.3969/j.issn.1005–0388.2009.01.001
[84] Zhang X M, Fidani C, Huang J P, Shen X H, Zeren Z, and Qian J. Burst increases of precipitating electrons recorded by the DEMETER satellite before strong earthquakes[J]. Natural Hazards and Earth System Sciences, 2013, 13 : 197 –209. DOI: 10.5194/nhess–13–197–2013
[85] Zhang X M, Liu J, Shen X H, Parrot M, Qian J D, Ouyang X Y, Zhao S F and Huang J P. Ionospheric perturbations associated with the M8.6 Sumatra earthquake on 28 March 2005 [J]. Chinese Journal of Geophysics, 2010b, 53 (3) : 567 –575. DOI: 10.3969/j.issn.0001–5733.2010.03.010
[86] 张学民, 刘静, 赵必强, 徐彤, 申旭辉, 姚璐. 玉树地震前的电离层异常现象分析[J]. 空间科学学报, 2014, 34 (6) : 822 –829. Zhang X M, Liu J, Zhao B Q, Xu T, Shen X H and Yao L. Analysis on ionospheric perturbations before Yushu earthquake[J]. Chinese Journal of Space Science, 2014, 34 (6) : 822 –829. DOI: 10.11728/cjss2014.06.822
[87] Zhang X M, Shen X H, Parrot M, Zeren Z, Ouyang X, Liu J, Qian J, Zhao S and Miao Y. Penomena of electrostatic perturbations before strong earthquakes (2005-2010) observed on DEMETER[J]. Natural Hazards and Earth System Sciences, 2012, 12 : 75 –83. DOI: 10.5194/nhess–12–75–2012
[88] Zhang X M, Shen X H, Liu J, Ouyang X Y, Qian J D and Zhao S F. Ionospheric perturbations of electron density before the Wenchuan earthquake[J]. International Journal of Remote Sensing, 2010a, 31 (13) : 3559 –3569. DOI: 10.1080/01431161003727762
[89] 张学民, 申旭辉, 赵庶凡, 刘静, 欧阳新艳, 娄文宇, 泽仁志玛, 何建辉, 钱庚. 地震电离层探测技术及其应用研究进展[J]. 地震学报, 2016, 38 (3) : 356 –375. Zhang X M, Shen X H, Zhao S F, Liu J, Ouyang X Y, Lou W Y, Zeren Z, He J H and Qian G. The seismo-ionospheric monitoring technologies and their application research development[J]. Acta Seismologica Sinica, 2016, 38 (3) : 356 –375.
[90] 赵庶凡, 申旭辉, 潘威炎, 张学民. VLF波及大气层到低电离层的传输特性分析[J]. 空间科学学报, 2011, 31 (2) : 194 –200. Zhao S F, Shen X H, Pan W Y and Zhang X M. Penetration characteristics of VLF wave from atmosphere into the lower ionosphere[J]. China Journal of Space Sciences, 2011, 31 (2) : 194 –200.
[91] Zebker H A, Rosen P A, Goldstein R M, Gabriel A and Werner C L. On the derivation of coseismic displacement fields using differential radar interferometry: The Landers earthquake[J]. Journal of Geophysical Research: Solid Earth banner, 1994, 99 (19) : 19617 –19634. DOI: 10.1029/94JB01179
[92] Zhou C, Liu Y, Zhao S F, Liu J, Zhang X M, Huang J P, Shen X H, Ni B B and Zhao Z Y. An electric field penetration model for seismo-ionospheric research[J]. Advances in Space Research, 2017, 60 (10) : 2217 –2232. DOI: 10.1016/j.asr.2017.08.007
[93] Zhou X, Cambiotti G, Sun W K and Sabadini R. Co-seismic slip distribution of the 2011 Tohoku (Mw 9.0) earthquake inverted from GPS and space-borne gravimetric data [J]. Earth and Planetary Physics, 2018, 2 (2) : 120 –138. DOI: 10.26464/epp2018013
[94] Zhou X, Sun W K, Zhao B, Fu G Y, Dong J and Nie Z S. Geodetic observations detecting coseismic displacements and gravity changes caused by the Mw=9.0 Tohoku-Oki earthquake[J]. Journal of Geophysical Research, 2012, 117 (B5) : B05408 . DOI: 10.1029/2011JB008849