空間大地測量技術在火山形變監測中的應用1

胡亞軒　許建東　劉國明　宋尚武　李煜航

?

空間大地測量技術在火山形變監測中的應用1

胡亞軒1）許建東2）劉國明3）宋尚武1）李煜航1）

1）中國地震局第二監測中心，西安 710054 2）中國地震局地質研究所，北京 100029 3）吉林省地震局，長春 130022

巖漿活動的不同階段引起地表變化不同。地表形變受壓力源大小、形狀、深度及巖漿運移速率等影響；另外火山類型不同，地形不同，形變特征也不同。地表形變幅度范圍很大，為1×10-7—1米量級。火山區形變監測可以了解火山活動狀態，有助于進行噴發危險性的預測預報。形變監測從20世紀60年代的傳統技術逐漸過渡到20世紀90年代發展起來的GNSS和InSAR等大地測量新技術，火山區形變時空監測能力得到提高，同時縮短了預測時間。我國火山形變監測開始較晚，現已在長白山天池、騰沖以及海南等主要火山區開展監測。傳統的連續測量以地傾斜觀測為主；新技術主要以流動GNSS監測為主，連續觀測站少，InSAR技術研究時間密度不夠；目前形變監測還不能實現很好的時空覆蓋。

GNSS InSAR 形變模型 形變監測網

引言

全球陸地上已知活火山超過1500座，火山噴發給人們的生命財產帶來很大威脅。世界范圍內有500萬以上的人口居住在可能受火山噴發嚴重影響的地區，而且這些地方的人口還在不斷增長。20世紀90年代至今，仍有許多火山在噴發。2015年全球因火山噴發產生的難民人數近30000人（張傳杰等，2016）。火山災害預測在日本、美國等一些多火山的國家普遍受到高度重視，建立了完善的火山監測系統，研究工作也比較深入（許建東，2006）。我國于1988年建立五大連池火山觀測站，這些觀測站觀察火山變化（包括蒸汽噴發、植物死亡等），記錄火山地震活動，測量火山形態、地溫、重力的變化，研究活動火山的噴發過程（洪漢凈，2013）。地震監測是觀測站最常用的方法，其次是形變測量（Scarpa等，2001）。火山形變可以反映巖漿的累積或后撤、巖墻的侵入、火山側翼的活動、火山區斷層活動及地震等。形變測量對火山噴發活動的預測預報有很大幫助。如1980年5月18日噴發的美國圣海倫斯火山，噴發前觀測到大范圍的隆升，一些區域隆升大于150m；4月10日觀測到50urad/h的傾斜等明顯形變異常。形變技術應用于監測火山活動已經非常久遠，如日本櫻島火山于1891年就開始開展水準測量，還進行三邊、三角及地傾斜等傳統方法的測量（Dvorak等，1997）。隨著形變觀測技術的不斷發展，火山形變監測方法也在發生變化，全球衛星導航系統（GNSS）和合成孔徑雷達干涉測量（InSAR）等新技術逐漸取代傳統技術。現代新大地測量技術被廣泛應用于火山監測，Fournier等（2010）統計的全球233個活動火山中，93%以上采用了新技術。

中國現代火山有長白山天池、騰沖、海南和五大連池等，具有噴發時距長、噴發頻率低和噴發類型單一等特點，但現在仍不同程度地存在著噴氣、地熱和火山型地震活動等特點（李玉鎖等，1998）。中國從20世紀后期才開始對火山開展系統監測，監測歷史短、數據積累少。與西方多火山的發達國家相比，我國開展系統監測以來還沒經歷過火山噴發事件，反映火山活動強度異常信號的背景基線尚在探索之中，還需要借鑒國外典型火山的監測背景基線進行類比。本文在總結國外火山形變研究結果的基礎上，結合模型分析形變監測技術在中國主要火山區的應用現狀及形變新技術的應用前景。

1 形變監測新技術

巖漿從深部向淺層巖漿房侵入往往引起短暫或永久膨脹；火山噴發期間，伴隨巖漿從地殼向地表的噴溢同樣引起地表形變的變化。隨著巖漿的噴出或向裂隙等的運移，往往引起地表收縮。火山類型不同，地形不同，形變特征也各異。盾狀火山山體通常具有寬闊頂面和緩坡度（小于10°）側翼，如冒納羅亞火山底徑達150km；冰島火山錐體高度很少超過0.1km。層狀/復合火山經過多次噴發形成，由熔巖流和火山碎屑物交替組成，熔巖粘滯性較高，如圣海倫斯火山長15km，高程3km左右。復合火山由多個錐體組成一個大的火山，其中一個主火山構成復合火山的主體。錐形火山（熔巖穹丘）為高粘滯性富硅巖漿緩慢擠出，相對高差0.05—0.75km，坡度15°—30°，如美國亞利桑那州的落日火山長1.5km，高0.3km。Dvorak（1992）研究表明，火山活動在距火山山頂約10km距離范圍內，至少會出現幾厘米的地表形變；如果在火山頂一定距離內布設測點，對于盾狀火山，這種較大范圍的形變相對容易被傳統的形變技術監測到，如夏威夷的基拉韋厄和冒納羅亞，冰島的克拉夫拉火山，以及相對高差不很大的層狀火山，如長谷火山和櫻島火山等；而一些火山由于山坡陡峭，大的地面運動僅發生在距火山頂約lkm的范圍內，如圣海倫斯山，埃特納火山等。常規的形變測量方法通常要求所選測點之間有清楚的視線，故很難觀測到發生在遠距離或測線崎嶇的局部火山運動。

形變監測方法多，技術發展快，不同方法觀測重復周期不同，可以用于不同形變幅度的測量（圖1）。20世紀60年代到80年代，火山區主要采用傳統的形變監測技術，連續形變測量采用傾斜儀和應變儀，應用水準測量和三角測量分別得到垂直形變及水平形變。20世紀80年代，隨著測距儀精度的提高，三邊測量方法逐漸取代三角測量方法。20世紀90年代以來，隨著美國的GPS、俄羅斯的Glonass、歐洲的Galileo、中國的北斗等GNSS，日本QZSS、印度IRNSS等區域系統以及相關增強系統等定位技術的出現和發展，可同時實現多衛星導航系統（Multi-GNSS）的監測。美國夏威夷冒納羅亞火山山頂的2個水準觀測點間距離測量采用了電子激光測距儀（EDM），監測時段內發生1975和1984年2次噴發；2000年后采用GPS連續測量，觀測數據有很好的時間覆蓋，2002年5月基線由縮短轉為伸長，后發生深源地震；2014年中期以后持續伸長約6cm，推測由于巖漿充填了火山山頂下方及西南裂隙引起的。夏威夷納帕烏火山口附近2個距離最近的GPS觀測站基線長度在1997年1月出現36cm的拉長，巖漿侵入到東側裂隙并于1月30日噴發，該變化持續到31號；一半以上的觀測量出現在噴發前8小時，也說明巖漿8小時后運到地表。同時，InSAR技術也迅速崛起，可以提供毫米—厘米級的三維形變。

應用GPS方法，國內外學者已取得了大量的火山水平變形監測資料。冰島已于1986年在火山區開展GPS觀測；Sigmundsson等（1992）監測到海克拉火山裂隙噴發引起的地表形變；Marshall等（1997）從1990—1994年，每年重復GPS觀測資料分析得到長谷火山巖墻張裂引起的地表膨脹。GPS監測資料在地下巖漿壓力變化過程研究判定方面提供了極大幫助。Owen等（2000a，2000b）分析了夏威夷火山1990—1996年70個GPS觀測站的重復觀測資料，得到了水平速率，結果表明1999年1月—2002年5月莫納羅亞火山口周圍幾個觀測點位移速率很小，速率差也很小，推測該階段無明顯巖漿活動；而2002年5月—9月火山區出現明顯膨脹，推測巖漿源發生了較強的壓力增大過程。

GNSS的不同觀測方法各有優缺點，適用于火山不同的活動階段。流動觀測是在一定時間段內對同一監測點的重復測量，能比較快地監測較大區域，不需要長時間的場地設備維護，但是時間密度不夠，可用來監測火山水平變形階段性變化，評價火山危險性，適用于平靜期或初始擾動期的火山變形監測。固定臺站的連續觀測可實時監測火山水平變形的動態變化過程，但通常點位密度不夠，適用于地下巖漿活動強烈、噴發危險性大的火山。1988年3月在伊豆半島建設的GPS連續觀測網為了監測該區地震，是最早建設的連續觀測網之一。觀測資料顯示相距10km的2個觀測點在地震活動期間伸長14.5cm，揭示為巖漿充填裂隙，并在10天后海底火山噴發（Dzurisin，2003）。隨后長谷、蘇弗里耶爾、基拉韋厄等火山也開展連續GPS觀測（CGPS）（Battaglia等，2006）。實時動態測量（The Real-Time Kinematic，RTK）技術也被應用于火山活動的實時監測，其相對費用高，達厘米級精度，監測基線長度范圍小于10—20km，在拉包爾和三宅島等火山都有應用和研究。拉包爾火山采用RTK技術由流動站倒推基本站位置，德國由此發展了基于GPS的在線控制及預警系統。采用多種觀測方法可更好地監測火山活動的不同階段，如夏威夷島布設有16個CGPS觀測站和100多個流動觀測站，用來監測基拉韋厄、莫納羅亞和莫納克亞等火山的活動。

圖1 不同觀測技術重復周期及測量幅度

InSAR技術可以監測大區域變形，彌補火山區沒有形變監測或觀測點少的情況，更適合對盾狀火山和大的破火山口進行監測。Hooper等人（2004）利用ERS影像生成的21幅干涉圖提取了美國加州長谷火山的形變時間序列，與GPS、水準、EDM的監測結果很好地吻合。Briole等（1997）利用1992—1993年ERS-1 SAR影像資料得到埃特納火山1986—1987及1989年中等大小和強度的巖漿流動。Amelung等（2000）從InSAR獲得的高空間分辨率位移圖提供了厄瓜多爾加拉帕戈斯群島火山群地表下巖漿積聚的最初異常，結果表明在伊莎貝拉島和費爾南迪納島上，七分之六的火山在1992—1999年發生了形變。InSAR技術也會受到一定限制。對于層狀/復合火山，接近山頂處存在大的幾何變形和冰雪覆蓋等不利觀測環境，會導致失相干，進而影響數據質量。早期的InSAR研究成果中，一些火山區由于可用的衛星干涉圖少等原因，結果并不理想，如像云仙岳火山于1990—1995年噴發，在1992—1993期間的70天，C波段干涉圖僅在北翼得到小部分的干涉（Zebker等，2000）。另外，短波段信號更容易丟失；長波段更易穿透植被，相干性好，但對電離層延遲更敏感。21世紀以來隨著發射和運行衛星的增多，存在C、L和X等不同波段。長白山火山一年中多數時間被積雪覆蓋，故常采用長波段衛星資料及一些特殊的處理方法來分析地表形變變化。InSAR新技術方法包括永久散射體（PS）技術、人工角反射器（CR）輔助技術、小基線集（SBAS）技術、干涉點目標分析（IPTA）技術等。法國學者Peltier等人利用PSInSAR技術提取了富爾奈斯火山2003—2009年的形變時間序列，與GPS監測結果具有很好的一致性（Peltier等，2010）。Bato等人（2011）利用PSInsAR技術提取了菲律賓馬榮火山的形變速率，與GPS結果一致。此外，未來還會發射Radarsat Constellation、SOACOM-1a/b等升級衛星，可克服InSAR技術應用的缺陷，實現互補，為技術的推廣應用奠定基礎。

2 形變模型及監測網布設

火山區地表形變受壓力源形狀、大小及深度等影響，其變化特征不同，可以用球形、橢球形、巖床和巖墻等模型進行模擬。火山噴發時，壓力源處于淺位，會在相對小的區域引起大的形變（Dvorak等，1997）；噴發間歇期通常由于壓力變化引起更大區域每年幾毫米到幾厘米的位移變化，如基拉韋厄和冰島的一些火山在間歇期形變量主要表現為小于1m。另外，火山類型不同，地形不同，可采用不同的形變技術進行監測。

Mogi模型是模擬地表形變最簡單適用的模型（Kiyoo，1958），是McTigue（1987）公式的近似，對有限大小的球型源適應性較好，一般源深度需大于半徑的3倍以上。橢球形模型是針對壓力源的長短軸近似橢球形。當巖漿充填裂隙、斷層時，會形成水平的巖床或垂直的巖墻，需要用不同模型進行模擬。Mogi模型假設地形為平坦的，當火山坡度超過20°時需要考慮地形的影響。Dvorak等人（1997）統計21座火山巖漿源深度范圍在1.3—13km。火山噴發往往是巖漿向淺層侵入結果，統計結果表明一般壓力源深度在3—5km，5km也是大多數火山的脆—韌轉換帶深度（Hill，1993）。火山噴發規模不同，造成的災害影響不同，對于強烈噴發的火山，噴發體積在0.1—1km3。假定噴發物體積為壓力源體積的變化，則半徑約為287—620m。取壓力源深度為5km，等效體積變換半徑為500m，得到不同高程（）的高斯型火山地表形變，如圖2所示（圖中為徑向距離（km）；為源深（km）；h為水平形變（km）；v為垂直形變（km）；vmax為垂直形變最大量（km）；V(H)=Hexp(-(-)2/22)，代表不同高程的火山，其中：為觀測點至壓力源中心的距離（km），為位置參數，為形狀參數，H分別取0=0、1=500m、2=1500m）。

圖2 Mogi模型引起不同高程高斯型火山的地表形變

由圖2可以看出火山區的形變變化范圍很大。距離壓力源中心越遠，壓力源體積變化越小，形變量越小，需要更多高精度的測量手段。距離壓力源中心的形變受地形影響也很明顯。火山形變監測站布設需要很好的空間分布，Dzurisin（2003）建議10km范圍內至少有1—3個觀測站，10—20km范圍內至少1—2個，這樣一般火山區共需設8—20個觀測站。從歷史記錄來看（Dvorak等，1997），一些火山翼的裂隙噴發，形變變化影響的距離最大可達50km。

對于Mogi模型，如果取壓力源深度分別為5km和2km（圖3），當半徑變化150m，垂向位移在距壓力源中心0—5km處最大，為3—6cm；20km處僅0.1cm左右。對于水平位移，10km處形變量在2—3cm；40km遠處仍存在幾毫米的變形。因此，合理的形變監測網應當滿足足夠的時長，以確保記錄到事件全過程；足夠的儀器靈敏度記錄微小和緩慢的形變變化，足夠的空間覆蓋記錄整個事件（Dvorak等，1997）。考慮到三維變形，對于活動強的火山，可以在火山區不同距離和方向上布設20個左右的CGPS站進行實時形變監測。噴發間歇期的合理監測方法是在近火山區域開展一年1次或2次觀測，而遠離火山區的地方觀測頻率可以少一些（Sigurdsson等，2015）。

圖3 不同深度的壓力源引起的地表形變

火山形變監測的優劣取決于觀測數據的時空分布、觀測精度等。地表位移變化范圍非常大，在穹狀或中酸性巖漿的火山中可達數十米，火山造錐階段形變量可達百米以上。圣海倫斯火山1980年噴發后進入造錐階段，截至1987年生長的穹丘達910m寬、240m高（Harris等，1988），并且在1989—1991年有幾個月的噴發活動。當巖漿運動緩慢或在很深位置運移時，位移的年變化量級僅為10-6m，甚至小于10-7m。新技術的應用可以克服傳統測量不能監測到微小地表形變的缺點。

綜合利用多種形變監測技術可以克服各自缺點，更好地監測火山活動。合理的GNSS網能夠監測到由于深部或分散巖漿源引起的地表形變，GNSS觀測精度高，但需要大量觀測點以確保測量范圍可以覆蓋整個火山。InSAR不需要在地面上配置儀器，可以實現很好的空間覆蓋，但是受植被、冰雪的影響，不易獲得小的地表形變。為了獲得1cm/a的形變（Sentinel-1衛星），在卡斯卡特火山區需要60—130天的觀測時間（Parker等，2015）。另外，所需時間與火山類型有關，不同類型的火山噴發前擾動持續時間不同，所需的連續衛星影像個數會不同，如盾狀火山異常擾動會持續40天，需要10個連續的干涉影像（Sentinel-1衛星）。傾斜儀和應變儀可進行連續觀測，對于小的形變比GPS和InSAR更敏感，但易受溫度、降雨等環境影響，能更快得到形變變化，但需要去除地球潮汐和小的非火山活動信號。綜合多種監測技術可以獲取更好的形變時空分布，如冰島全國火山監測網就是由40臺地震儀、17個GPS固定站和10臺100—400m深的鉆孔應變儀組成。

3 中國主要火山的空間大地測量技術監測

中國主要火山形變監測網建站時間短，觀測站少（許建東，2011）。目前中國僅在長白山天池、騰沖、海南3個火山區建有火山形變監測網，開展了水準、GNSS、InSAR及地傾斜等的觀測與研究工作，監測情況見表1。另外，主要服務于地震預測預報的中國地殼運動觀測網絡工程（CMONOC I）、中國大陸構造環境監測網絡（CMONOC II）是以GNSS為主的國家級地球科學綜合觀測網絡，由260多個連續觀測站（CGPS）及2000多個流動觀測的區域站組成，其中部分站點覆蓋了火山區，還有部分用于地震監測的臺站距離火山區較近，布設有傾斜儀及應變儀等進行定點形變連續觀測，這些監測數據均可服務于火山預測預報。

表1 大地測量新技術應用于中國活動火山監測情況統計

3.1 長白山天池火山

長白山天池火山是一座巨型復合層狀火山錐體，天池西側中國境內的主峰白云峰海拔2691m，錐體半徑20km。天池火山有歷史記錄以來曾發生過多次噴發活動，最近一次噴發活動發生在1903年，目前已經休眠達100余年。洪漢凈等（2007）依據火山活動的危險性將長白山火山歸為Ⅰ類，認為其是位于深震活動區上方的弧后板內火山，仍然具有產生巖漿活動的構造條件。2002—2005年多種監測資料表明火山經歷了一次擾動過程（Xu等，2012）。火山區常規形變監測有水準、GPS和定點形變觀測，見圖4，布設2條水準路線獲取垂直形變（李克等，2009；劉國明等，2011），北坡水準路線長24km，相對高差901m，從2002年開始觀測，監測到2002—2003年最大垂向隆起46.3mm；西坡水準路線全長約30km，相對高差1084m，2006年開始觀測。GPS監測網由15個點組成，其中1999年建立8個點，并于2000年開始觀測，每年夏季進行1期流動觀測；2006年后補充建設7個點，2007年開始每年一次的觀測，全部點位基本圍繞天池火山口周圍。定點連續形變觀測主要是洞體內的傾斜儀和伸縮儀，從1999年12月開始觀測。

從圖4可以看出，在新技術應用方面，15個GNSS火山流動觀測站中，4個分布在10km范圍內、7個分布在10—20km范圍內，點位分布在各方向不均勻。距火山口約40km處，分布有CMONOC流動觀測站和CGPS點位各1個。長白山天池火山地區地形復雜、地勢高峻、積雪期長、植被非常茂密，很容易造成SAR干涉圖像失相干，需要特殊的影像及處理方法。韓宇飛等（2010）選取了L波段的JERS-1 SAR數據進行差分干涉處理（D-InSAR），提取了研究區1994—1998年的地表形變場；Ji等（2013）基于ENVISAT ASAR影像資料，采用永久散射技術（PSInSAR）獲取了天池火山2004—2010年的形變時間序列。可見，要全面監測火山區形變的時空演化，目前已有的形變監測時空密度還不夠，需要更多的CGPS觀測站進行監測，以及長期連續InSAR資料的應用及分析等。

圖4 長白山火山形變監測網（TC-天池）

3.2 云南騰沖火山

騰沖火山最新的一次活動發生于全新世，形成打鷹山、老龜坡、馬鞍山和黑空山4個截頭圓錐狀火山錐。其中打鷹山最高，海拔2595m，火山錐體長2km、寬1.5km、高0.6km（皇甫崗等，2000），1609年仍有噴發活動，該火山區溫泉發育、地震頻繁，表明火山仍有活動性；馬鞍山海拔1793m，火山錐體長0.8km、寬0.6km、高0.11km；黑空山海拔2072m，火山錐體長0.8km、寬0.8km、高0.15km。

火山主要活動區域已開展水準、重力及GPS等觀測。1997年建成了由64個測點組成的水準網（圖5），路線總長約150km，基本覆蓋主要火山，2003年增設了熱海—勐連—騰沖水準測線85km，共30個點，1998—2004年共進行了5期水準測量（胡亞軒等，2003，2007）。在GPS監測方面，分別于2002、2003年對原激光測距網（1997年建立）進行了觀測；2003年建成了由20點組成的新GPS觀測網，并于2003、2004年開展了野外觀測。

從圖5可以看出，GNSS火山觀測站基本分布在20km范圍內；距離火山口約40km范圍內有CMONOC流動觀測站9個和CGPS點位1個（王伶俐等，2015），其中部分點位可用于火山監測。騰沖火山區地貌復雜、所處緯度偏低、植被覆蓋茂密，很容易造成SAR影像失相干。季靈運等（2011）基于JERS-1Ｌ波段SAR影像，利用小基線集—合成孔徑雷達干涉測量（SBAS-DInSAR）技術提取了騰沖火山地區1995—1997年間地表形變時間序列。可見，騰沖火山區目前形變監測時空密度還不夠，也需要更多的CGPS觀測站和InSAR等資料的長期連續監測分析。

圖5 騰沖火山形變監測網（HKS-黑空山；DYS-打鷹山；MAS-馬鞍山）

3.3 海南火山

海南島北部熔巖覆蓋面積達4000余平方千米（胡久常，2002；白志達等，2003）。第四紀火山區的石山、永興一帶大小30多個火山口明顯呈北西方向排列，形成典型的中心式火山群，是中國第四紀以來火山活動最為強烈和頻繁的地區之一，其中雷虎嶺和馬鞍嶺為典型的現代火山口。白志達等（2003）研究認為，馬鞍嶺地區第四紀火山活動具有多期性，不同期次具有不同的火山活動方式、噴發強度及火山結構類型。海口地區火山噴發活動較弱，海拔最高222.8m，火山錐體高度在10—40m，最大底徑不足3km，絕大部分火山底徑在1km以下。

從圖6監測網看出，火山區共有6個GPS觀測站，從2008年開始每年進行一次觀測；另有CMONOC和中國地震局第二監測中心建設的8個GNSS流動觀測站、1個CGPS站可供觀測，但個別觀測站受周圍樹木及（或）建（構）筑物影響，觀測數據質量差（Hu等2016a，2016b）；此外，火山區西部觀測點位較少。受茂密的植被影響，Ji等（2015）采用C波段的ENVISAT ASAR影像不能得到干涉圖，后采用L波段的ALOS PALSAR衛星影像利用干涉圖堆疊技術（stacking-InSAR）獲取了2007—2011年的形變演化圖。目前火山區缺乏對微小形變敏感的連續形變監測手段；雷虎嶺和馬鞍嶺等主要活動火山區，監測點位少，形變監測時空密度還不夠。

圖6 海南火山形變監測網（MAL-馬鞍嶺；LHL-雷虎嶺）

4 結論

火山區地表形變很復雜，地表位移的變化和諸多因素有關，通常由于地下巖漿運動引起，可能與火山地熱系統中的流體壓力與流速變化有關。火山區形變監測可幫助人們了解火山的活動動態，火山的噴發周期可能經歷部分熔化、初始上升、地殼同化、巖漿混合、存儲、部分結晶和最終上升到地表等全部或部分過程，從而引起的形變特征也不同；形變變化在火山噴發前存在很大的差異，從零到十多米直至幾百米，影響范圍也不同。為了獲得可靠的形變時空演化特征，可以應用水準、GPS和InSAR獲取高精度大范圍的觀測資料；應用應變儀、傾斜儀以及CGPS獲取理想的形變隨時間演化特征。多種監測技術在火山區的聯合應用，可以更好地觀測到由于巖漿侵入引起的大范圍非震形變時空特征，有利于精確綜合判斷火山活動狀況，以及對噴發前數天到幾個月的短期預報和災害評估。中國主要火山采用的形變監測技術較少，目前火山區域的GNSS觀測站密度還不能很好滿足監測火山活動的需要，用于小形變的連續觀測站還不夠，監測網還需要進一步優化，此外，InSAR技術在中國火山僅進行了部分時段的研究工作。因此，建議對計劃開展形變監測的火山，在規劃設計時能夠依據火山類型、活動性等考慮所采用的監測技術、建站密度及觀測頻次等。

白志達，徐德斌，魏海泉等，2003. 瓊北馬鞍嶺地區第四紀火山活動期次劃分. 地震地質，25（S1）：12—20.

陳國滸，單新建，Moon W. M. 等，2008. 基于InSAR、GPS形變場的長白山地區火山巖漿囊參數模擬研究．地球物理學報，51（4）：1085—1092.

韓宇飛，宋小剛，單新建等，2010. D-InSAR技術在長白山天池火山形變監測中的誤差分析與應用. 地球物理學報，53（7）：1571—1579.

洪漢凈，吳建平，王慶良等，2007. 中國火山危險性等級與活動性分類. 地震地質，29（3）：447—458.

洪漢凈，2013. 火山預測與預警. 北京：地震出版社，1—278.

胡久常，2002. 海南有座休眠活火山. 防災博覽，（6）：43.

胡亞軒，施行覺，王慶良等，2003. 騰沖火山區地表垂直形變分析. 大地測量與地球動力學，23（2）：37—41.

胡亞軒，王慶良，崔篤信等，2007. 用形變資料分析騰沖火山區巖漿的活動特征. 地震研究，30（2）：164—168.

皇甫崗，姜朝松，2000. 騰沖火山研究. 昆明：云南科技出版社.

季靈運，王慶良，崔篤信等，2011. 利用SBAS-DInSAR技術提取騰沖火山區形變時間序列. 大地測量與地球動力學，31（4）：149—153，159.

季靈運，許建東，趙波等，2013. 利用InSAR技術研究新疆阿什庫勒火山群現今活動性. 地震地質，35（3）：532—541.

李克，劉俊清，盤曉東等，2009. 2000—2007年期間長白山天池火山區地殼變形監測與分析. 地震地質，31（4）：639—646.

李玉鎖，修濟剛，李繼泰等，1998. 火山噴發機制與預報. 北京：地震出版社.

劉國明，孫鴻雁，郭峰，2011. 長白山火山最新監測信息. 巖石學報，27（10）：2905—2911.

Scarpa R.，Tilling R. I.，1996.火山監測與減災. 劉若新等譯，2001. 北京: 地震出版社，285—303.

唐攀攀，單新建，王長林，2014. 基于PSInSAR技術的長白山天池火山形變監測. 地震地質，36（1）：177—185.

王伶俐，邵德盛，洪敏，2015．云南境內陸態網絡GNSS觀測資料數據處理與初步結果．震災防御技術，10（1）：141—150．

許建東，2006. 我國火山災害的主要類型及火山災害區劃圖編制現狀探討. 震災防御技術，1（3）：266—272.

許建東，2011. 中國活動火山監測進展回顧. 礦物巖石地球化學通報，30（4）：390—392.

許建東，趙波，Sindney H. 等，2014. 西昆侖阿什庫勒火山群地質特征和活動分期. 巖石學報，30（12）：3521—3530.

張傳杰，李霓，龔麗文，2016. 2015年全球火山活動綜述. 國際地震動態，（11）：4—11.

Amelung F., Jónsson S., Zebker H., et al., 2000. Widespread uplift and ‘trapdoor’ faulting on galápagos volcanoes observed with radar interferometry. Nature, 407 (6807): 993—996.

Bato M. P., Lagmay A. A., Paguican E. R., 2011. Interferometric SAR persistent scatterer analysis of Mayon volcano, Albay, Philippines. In: American Geophysical Union, Fall Meeting. Washington: American Geophysical Union, G23A-0848.

Battaglia M., Vasco D. W., 2006. The search for magma reservoirs in Long Valley Caldera: single versus distributed sources. Geological Society, London, Special Publications, 269 (1): 173—180.

Briole P., Massonnet D., Delacourt C., 1997. Post-eruptive deformation associated with the 1986-87 and 1989 lava flows of etna detected by radar interferometry. Geophysical Research Letters, 24 (1): 37—40.

Dvorak J. J., 1992. Tracking the movement of Hawaiian volcanoes; Global Positioning System (GPS)measurement. Earthquakes & Volcanoes, 23 (6): 255—267.

Dvorak J. J., Dzurisin D., 1997. Volcano geodesy: the search for magma reservoirs and the formation of eruptive vents. Review of Geophysics, 35 (3): 343—384.

Dzurisin D., 2003. A comprehensive approach to monitoring volcano deformation as a window on the eruption cycle. Reviews of Geophysics, 41 (1): 1001.

Fournier T. J., Pritchard M. E., Riddick S. N., 2010. Duration, magnitude, and frequency of subaerial volcano deformation events: new results from Latin America using InSAR and a global synthesis. Geochemistry Geophysics Geosystems, 11 (1): Q01003.

Harris S. L., 1988. Fire mountains of the west: the cascade and mono lake volcanoes. Missoula: Mountain Press Publishing Company.

Hill D. P., 1993. Temperatures at the base of the seismogenic crust beneath Long Valley caldera, California, and the Phlegrean Fields caldera, Italy. In: Gasparini P., Scarpa R., Aki K., Volcanic Seismology. Berlin, Heidelberg: Springer-Verlag, 432-461.

Hooper A., Zebker H., Segall P., 2004. A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers. Geophysical Research Letters, 31 (23): L23611.

Hu Y. X., Hao M., Ji L. Y., et al., 2016a. Three-dimensional crustal movement and the activities of earthquakes, volcanoes and faults in Hainan Island, China. Geodesy and Geodynamics, 7 (4): 284—294.

Hu Y. X., Cheng L., Wang X., 2016b. Quality analysis of the campaign GPS stations observation in Northeast and North China. Geodesy and Geodynamics, 7 (2): 87—94.

Ji L. Y., Xu J. D., Wang Q. L., et al., 2013. Episodic deformation at Changbaishan Tianchi volcano, northeast China during 2004 to 2010, observed by persistent scatterer interferometric synthetic aperture radar. Journal of Applied Remote Sensing, 7 (1): 073499.

Ji L. Y., Hu Y. X., Wang Q. L., et al., 2015. Large-scale deformation caused by dyke intrusion beneath eastern Hainan Island, China observed using InSAR. Journal of Geodynamics, 88: 52—58.

Kiyoo M., 1958. Relations between the eruptions of various volcanoes and the deformations of the ground surfaces around them. Bulletin of the Earthquake Research Institute, University of Tokyo, 36 (2): 99—134.

Marshall G. A., Langbein J., Stein R. S., et al., 1997. Inflation of Long valley caldera, California, basin and range strain, and possible mono craters dike opening from 1990-94 GPS surveys. Geophysical Research Letters, 24 (9): 1003-1006.

McTigue D. F., 1987. Elastic stress and deformation near a finite spherical magma body: resolution of the point source paradox. Journal of Geophysical Research: Solid Earth, 92 (B12): 12931—12940.

Owen S., Segall P., Lisowski M., et al., 2000a. January 30, 1997 eruptive event on Kilauea volcano, Hawaii, as monitored by continuous GPS. Geophysical Research Letters, 27 (17): 2757-2760.

Owen S., Segall P., Lisowski M., et al., 2000b. Rapid deformation of Kilauea Volcano: global positioning system measurements between 1990 and 1996. Journal of Geophysical Research: Solid Earth, 105 (B8): 18983—18998.

Parker A. L., Biggs J., Walters R J., et al., 2015. Systematic assessment of atmospheric uncertainties for InSAR data at volcanic arcs using large-scale atmospheric models: application to the cascade volcanoes, United States. Remote Sensing of Environment, 170: 102—114.

Peltier A., Bianchi M., Kaminski E., et al., 2010. PSInSAR as a new tool to monitor pre-eruptive volcano ground deformation: Validation using GPS measurements on Piton de la Fournaise. Geophysical Research Letters, 37 (12): L12301.

Sigmundsson F., Einarsson P., Bilham R., 1992. Magma chamber deflation recorded by the Global Positioning System: The Hekla 1991 eruption. Geophysical Research Letters, 19 (14): 1483—1486.

Sigurdsson H., Houghton B., McNutt S., et al., 2015. The encyclopedia of volcanoes. 2nd ed. New York: Acadamic Press, 1101—1123.

Xu J. D., Liu G. M., Wu J. P., et al., 2012. Recent unrest of Changbaishan volcano, northeast China: A precursor of a future eruption? Geophysical Research Letters, 39 (16): L16305.

Zebker H. A., Amelung F., Jonsson S., 2000. Remote sensing of volcano surface and internal processes using radar interferometry. In: Mouginis-Mark P. J., Crisp J. A., Fink J. H. Remote Sensing of Active Volcanism. Washington DC: Wiley, 179—205.

Application of Geodesy Technology in Deformation Monitoring of Volcanoes

Hu Yaxuan1), Xu Jiandong2), Liu Guoming3), Song Shangwu1)and Li Yuhang1)

1）The Second Monitoring and Application Center, China Earthquake Administration, Xi’an 710054, China 2）Insitute of Geology, China Earthquake Administration, Beijing 100029, China 3）Jilin Earthquak Agency, Changchun 130022, China

The surface deformation varies with magmatic activities in different periods. Deformation monitoring can help to understand the activities and to evaluate the potential risk of eruption. Generally, the deformation in volcano area is affected by the magamatic source pressure, size, shape, depth and migration rate and etc. The volcanic types also can cause different deformation on the surface. The range of amplitude is from 1×10-7to 1 meter-scale. The new geodetic technologies developed in the 1990s, such as GNSS and InSAR, gradually replaced the traditional technique of the 1960s. The monitoring capability has been improved in temporal-spatial domain. The deformation monitoring of active volcanoes in China started late in 1990s, which has been used in the volcano monitoring of Tianchi Changbaishan, Tengchong and Hainan volcanoes. However, the continuous deformation measurement is done only with tilting. The application of new geodesy technologies, such as continuous GNSS, is not enough in these volcano regions.

GNSS; InSAR; Deformation model; Deformation monitoring network

胡亞軒，許建東，劉國明，宋尚武，李煜航，2018．空間大地測量技術在火山形變監測中的應用．震災防御技術，13（2）：410—423．

10.11899/zzfy20180216

國家自然科學基金（41372345）

2017-02-16

胡亞軒，女，生于1970年。高級工程師。主要研究方向：地形變機理。E-mail：happy_hu6921@sina.com