本研究為了瞭解對於玉山國家公園在九二一地震發生後環境變遷狀況,特別是地震直接造成的土地覆蓋變遷,運用遙感探測(Remote Sensing:RS)技術來進行玉山國家公園全區地震前、後期的土地覆蓋及變遷偵測,並判定地震受災情況以及潛在危險區域。
本研究分別使用衛星影像來偵測玉山國家公園全區在九二一地震發生後土地覆蓋變遷的大致狀況,以及使用航空攝影照片來針對主要受災地區---玉山國家公園西北部園區,即東埔、玉山地區,進行較為細緻的崩塌地判釋與分析工作。
由於玉山國家公園幅員遼闊,面積達105,490公頃,且地理位置處於中央山脈地帶,崇山峻嶺、溪谷深邃,許多地區人跡都無法到達,在九二一地震發生後,玉管處僅能針對主要步道及設施所在地,派員加以巡視調查,然而所能涵蓋面積卻極為有限,為了能夠有效率地針對整個園區,進行全面性的受災調查,本研究首先選定運用衛星影像為工具,透過遙感探測技術中的影像判釋相關技術,來進行環境變遷偵測,由於衛星影像拍攝位置高,因此地形起伏效應的影像較小,同時影像的涵蓋範圍大,都有助於減少資料處理的繁雜度,此外衛星通過同一地點的週期短(從數天到一個月),能夠有效地提供緊鄰地震發生前、後時間的地表覆蓋情形。
遙測衛星影像資料來源有許多,有最早1972年的Landsat Multi-spectral Scanner System (MSS)、1982年與1984年的Landsat 30公尺精度七個頻譜之Thematic Mapper(TM )、1986年至1993年的法國SPOT衛星之20公尺精度三個頻譜或全色態(黑白,panchromatic)10公尺精度、氣象衛星(Meteorological Satellite)、美國NASA太空船之HyperSpectral Imager (HIS)、……等。在運用遙測衛星觀測地球時偵測方式係利用不同頻率(波長)之波段,常用的波段整理如表3.1:
|
名稱 |
波長範圍 |
輻射來源 |
偵測方式 |
|
可見光(Visible,V) |
0.4-0.7 m m |
太陽能 |
反射 |
|
近紅外光(Near InfraRed, NIR) |
0.7-1.1 m m |
太陽能 |
反射 |
|
短波紅外光 (Short Wave InfraRed, SWIR) |
1.1-1.35 m m 1.4-1.8 m m 2-2.5 m m |
太陽能 |
反射 |
|
中波紅外光 (Mid Wave InfraRed, MWIR) |
3-4 m m 4.5-5 m m |
太陽能 熱能 |
反射 溫度 |
|
溫度的紅外光 (Thermal InfraRed, TIR) |
8-9.5 m m 10-14 m m |
熱能 |
溫度 |
|
微波(Microwave) 雷達波(radar) |
1 mm-1 m |
熱能(被動式) 人工製造(主動式) |
溫度(被動式) 粗糙度(主動式) |
本研究所使用的衛星影像是由國立中央大學太空及遙測研究中心所代為接收的法國衛星SPOT 影像,該影像的具有多光譜態(XS)及全色態(PAN)兩種能力,為了有效地針對地表覆蓋進行分類,本研究選用多光譜態影像,亦即由綠光段(0.50μm~0.59μm)、紅光段(0.61μm~0.68μm)及近紅外光段(0.78μm~0.89μm)所組成,而影像解析度在經過處理後可以像元大小達到12.5公尺*12.5公尺,而由於SPOT衛星具有傾斜攝影的能力,因此平均約每隔五天就可以對同一地點進行觀測。
原始的衛星影像資料,係由衛星系統拍攝,由於系統、大氣環境及地表起伏等因素都會對影響產生雜訊及扭曲現象,所以需要透過影像的雜訊消除、輻射糾正、幾何校正等前處理動作後,才得以獲得無幾何扭曲、幅射偏差且具有大地座標的影像,據此的影像判釋及分析成果,才能夠正確地反映出真實狀況,目前中大太遙中心有提供針對衛星影像代為前處理的服務,即所謂Level 10的產品,本研究為了有效率地進行後續研究工作,直接選用經過處理後的SPOT衛星影像來作為玉山國家公園境內土地覆蓋判釋的主要素材。
根據玉管處委託逢甲大學地理資訊系統研究中心所做的研究(周天穎,2000),針對賀伯風災後(1996年)及九二一地震後(1999年)的土地覆蓋,利用衛星影像加以分類所得的結果,如表3.2。
表3.2 玉山國家公園土地覆蓋面積變遷表(周天穎,2000)
|
1996年(賀伯風災後) |
1999年(九二一地震後) |
|||
|
面積(公頃) |
百分比(%) |
面積(公頃) |
百分比(%) |
|
|
草生地 |
12877.2 |
12.37 |
13238.1 |
12.72 |
|
裸露地 |
9356.33 |
8.99 |
8543.47 |
8.41 |
|
林地 |
79303.5 |
76.2 |
78818.6 |
75.73 |
|
水體(河川地) |
1816.08 |
1.74 |
1089.73 |
1.05 |
|
雲 |
- |
- |
2326.42 |
2.24 |
|
總面積 |
104077.4 |
100 |
104077.4 |
100 |
從研究中的結果發現,九二一地震後玉山國家公園園區整體而言土地覆蓋類型比例並未有重大的變化,但從較細微的觀點來看,從衛星影像上所判釋出的林地面積略有減小,但裸露地卻不增反減,根據檢視過1996年SPOT影像過後,發現風災過後所拍攝的影像,因為山區大量積水,土壤及樹冠潮濕的狀況,使得地表覆蓋物的光譜反射值受到水份的影響而產生變動,進而造成判釋上的誤差,此外由於兩個時期的衛星影像相隔三年,其間影響土地覆蓋變異的因子相當地多,無法單純地區隔出九二一地震所導致的土地覆蓋變異程度。
因此本研究另行選用更接近九二一地震發生前的SPOT衛星影像(1999年3月至8月間),來重新進行土地覆蓋分類,以其更精準地掌握九二一地震所造成土地覆蓋變遷狀況,同時藉以探討土地覆蓋變遷的空間分佈趨勢。茲將研究中所使用的影像資料表列如下:
|
編號 |
任務 |
日期 |
涵蓋範圍 |
備註 |
|
1 |
賀伯風災後 |
1996/08/18 |
全區 |
* |
|
2 |
九二一地震 |
1998/11/13 |
西北園區區 |
|
|
3 |
九二一地震前 |
1999/03/06 |
西北園區區 |
|
|
4 |
九二一地震前 |
1999/04/09 |
西南園區 |
|
|
5 |
九二一地震前 |
1999/08/16 |
東部園區 |
|
|
6 |
九二一地震後 |
1999/10/12 |
全區 |
* |
<*為玉山國家公園委託逢甲大學地理資訊研究中心研究案中所購置>
本研究所使用的分類方式上,採取監督式分類方法來加以操作,所謂監督式分類方法是利用光譜類型辨別以地面真實資訊,選定分析者所欲分析資料的類別,如不同類型的森林、不同的土地使用型態,然後選取訓練樣本區(training area),即,以代表每一類別的光譜特質,再從訓練區中導出相關統計值作為分類依據。此外,選用Purdue大學所發展的遙測影像分類軟體MultiSpec7.0作為分類工具,並以監督式分類方法下的最大概似法 (Maximum likelihood method)來達成土地覆蓋類型分類目的。
分類目標主要為了找出因為九二一地震使得裸露地的增長面積及分佈位置,因此設定分類類別時僅以裸露地、草生地、林地(向陽面)、林地(陰影面)及雲等五項類別作為分類之標的,之所以不那納入一般分類時會採用的人工設施(建築、道路)及水體的原因,是玉山國家公園境內人為開發少,因此人工設施所佔的比例極小,同時建築、道路、河床的光譜特徵與裸露地十分接近,區分不易,而道路及河川都有GIS圖形資料可供參考,可以進行完監督式分類之後再加以區隔出來,由於地形的效應導致太陽光照射強度之不同,林地在向陽面及陰影面亦會呈現不同的明暗對比,為求分類的精準度,在建立訓練樣本時將其視作獨立類別,待分類結束後,再將者兩各類別加以合併成單一林地類型。
將分類過後的影像以抽樣的方式檢核分類結果的正確性,在整體正確度Kappa統計值皆在85%以上,接受其分類結果。則兩個時期的土地覆蓋類別分類結果如表3.4
|
1999年3∼8月(九二一地震前) |
1999年(九二一地震後) |
|||
|
面積(公頃) |
百分比(%) |
面積(公頃) |
百分比(%) |
|
|
草生地 |
9,061.5 |
8.60 |
8790.7 |
8.35 |
|
裸露地 |
4,393.7 |
4.17 |
4886.0 |
4.64 |
|
林地 |
89,207.56 |
84.70 |
90275.2 |
85.7 |
|
雲 |
2659 |
2.71 |
1363.8 |
1.29 |
|
總面積 |
105321.7 |
100 |
105315.7 |
100 |
根據土地覆蓋類別分類結果,可以得知九二一地震並未對玉山國家公園境內的土地覆蓋類別造成巨大的變化,僅在裸露地部分略微增加約500公頃,但整體面積的0.47%,主要可能原因,來自於新增的崩塌地以及就崩塌地的擴張所累積出來的,而在草生地及林地上亦各有增減,前後時期的差異,主要原因,可能來季節差異造成植物生長狀況差異,導致產生判釋上之模糊地帶,同時由於兩個時期都有1~2%的雲量,阻礙了對地表的觀測,亦會對於面積統計上造成變動。
將地震前、後時期的土地覆蓋類別分類結果以圖形方式,分別呈現於圖3.1及圖3.2。
圖3.1 九二一地震前土地覆蓋分類圖
圖3.2九二一地震後土地覆蓋分類圖
由圖3.1及3.2相互比較,主要的土地覆蓋變異都發生在玉山國家公園的西北部,即東埔、玉山地區,在這個區域內地震過後的裸露地及草生地有都顯著增加的狀況,至於東部及西南部園區,扣除因為雲層及雲層陰影所造成的阻隔及判釋錯誤,則地表覆蓋類別及面積分佈皆無太大的變化,即便直接比對前後時期兩張衛星影像,仍不容易發現大面積的新崩塌地,這樣的結果與九二一地震發生時西北部地區所處的震度規模較大,相互契合。
根據GIS圖層套疊的概念將九二一地震前的裸露地分佈圖套疊在九二一地震後的衛星影像上,如圖3.3所示,底圖為九二一地震後的Spot衛星影像,青色的九二一地震前的裸露地判釋結果,同時可以比對九二一地震前的Spot衛星影像,於是很容易就可以以人工的方式找出主要的新增崩塌地位置,如圖所框選位置內即為新增的崩塌地。此外套疊上行政界線、地名、山名、水系及道路等輔助的空間資料,並能得知這些新的崩塌所處的位置,作為後續的航空照片判釋及田野調查的主要調查對象。將新的裸露地分佈位置標示如圖3.4,針對個別的裸露地區,呈現前、後期的衛星影像,並說明相關背景資訊如表3.5、表3.6、表3.7及表3.8。
圖3.3 九二一地震後衛星影像套疊地震前裸露地分佈圖
|
編號 |
位置 |
崩塌面積 |
崩塌類型 |
|
A |
陳有蘭溪右岸,雲龍瀑布旁,崩塌地阻斷日據八通關古道。 |
約37.3公頃 |
舊崩塌地擴張 |
|
B |
陳有蘭溪左岸,樂樂溫泉上方250公尺處。 |
約13.9公頃 |
新崩塌地 |
|
C |
東埔山西北方3公里處坡地上,出水溪右岸。 |
約13.1公頃 |
新崩塌地 |
|
D |
東埔山西方2.5公里處坡地上,出水溪左岸。 |
約11.8公頃 |
新崩塌地 |
|
E |
玉山西南方3公里處坡地,楠梓仙溪右岸,崩塌地阻斷塔塔加往玉山主峰步道。 |
約18.3公頃 |
舊崩塌地擴張 |
|
F |
南玉山西南方800公尺處,位於南玉山稜線左右兩側。 |
約56.9公頃 |
舊崩塌地擴張 |
|
位置:陳有蘭溪兩岸 |
坡度:大於55% |
|
|
面積:A約37.3公頃、B約13.9公頃 |
地震前土地覆蓋狀況:舊崩塌地與植生 |
|
|
描述:A、B兩崩塌地分屬陳有蘭溪兩岸,A崩塌地原本就已崩塌,地震後崩塌面積則更加擴大,阻斷東埔往玉山主峰的日據八通關越嶺道。B崩塌地原本多為植生,清八通關古道正通過此一位置,需進一步調查受損狀況。 |
||
|
|
||
|
地震發生前SPOT衛星影像 |
地震發生後SPOT衛星影像 |
|
|
位置:東埔山以西,新中橫公路下方 |
坡度:大於55% |
|
|
面積:C約13.1公頃、D約11.8公頃 |
地震前土地覆蓋狀況:植生 |
|
|
描述:此兩各崩塌地位於新中橫公路下方,分屬出水溪左、右兩岸,在地震之前地表都是屬於植生覆蓋狀態,所以是地震直接造成的新崩塌地。
|
||
|
|
||
|
地震發生前SPOT衛星影像 |
地震發生後SPOT衛星影像 |
|
|
位置:玉山主峰西南方 |
坡度:大於55% |
|
|
面積:E約18.3公頃、F約56.9公頃 |
地震前土地覆蓋狀況:崩塌地及草生地 |
|
|
描述:E、F兩崩塌地主要在南北向的玉山山脊線西方接近稜線的地方發生崩塌,原先都已存在崩塌現象,地震後崩塌面積則更加擴大。其中E崩塌地,阻斷塔塔加往玉山主峰的玉山主線步道,需額外留意該崩塌地是否持續發展。 |
||
|
|
||
|
地震發生前SPOT衛星影像 |
地震發生後SPOT衛星影像 |
|
在九二一地震後玉山國家公園,由衛星影像中所自動判釋出來的新增崩塌地約增加400~500公頃,再透過人工比對可以確切找出六個主要的新增崩塌地點,總計約為為150公頃,而從空間分佈上來看,這些新增的崩塌地都分佈在東埔、玉山群峰地區,其結果與九二一發生時該地區地震規模達到六級有密切關係。針對這些新增的崩塌地,將震災前的衛星影像及數值地形模型進行分析,可以發現九二一地震造成新的崩塌地,多位於坡度達55%以上的坡地,且震災前多已存在有舊的崩塌現象,而九二一地震的發生使得舊有的崩塌持續擴大。
這樣的結果與國家地震研究中心,以調查統計方式針對於全台九二一地震造成坡地破壞進行之分析(羅俊雄,2000),結果發現地震所造成得坡地破壞處皆發生於山脈經水流切割地形上坡度較陡且地表PGA值較高(即震度較高)之地方,兩者相符。由此可知,地震災害所造成的山崩、地滑現象除了與地震強度、地震延時特性相關外,山坡地的穩定度亦直接相關。因此針對全園區進行坡地災害敏感區的分析與分區劃定,並採取有效的經營管理措施,將會有助於減少地震所帶來的災害。
由於衛星影像的解析度僅12.5公尺,因此對於小的崩塌地並無法正確地判釋出來,因此有必要針對於東埔、玉山地區再進一步使用影像解析度更高的航空攝影照片,作更細緻的崩塌地判釋工作。
雖然衛星影像能夠有效率地提供全面性的環境變遷大致面貌,但是由於本研究使用的SPOT衛星影像的空間解析度僅12.5公尺,且在崩塌地判釋精確度上會受地形起伏、太陽入射角度及當時的土壤水分等因素所影響,同時對於於崩塌的類型及崩塌可能造成的原因,並無法直接提供有效資訊,仍須倚賴相關的輔助資訊,如數值高程資料、地質資料、水文資料等。針對於變遷程度大或環境敏感地區,為求更精細地掌握環境變遷的面貌,則有必要使用航空攝影照片來進行更進一步的判釋工作,但由於航空攝影照片往往無法配合災害發生後立即拍攝(研究區內的東埔玉山地區直到89年10月份才完成攝影工作),因此在時效性上仍有其缺點,同時航空攝影照片的購買成本高,因此本研究僅針對研究區內受到九二一地震導致土地覆蓋變遷程度較大的東埔玉山地區,進一步利用航空照片,再作更細緻的崩塌地判釋與分析,試圖找出該區域內直接因為由九二一地震所導致新的崩塌地之位置,一方面作為衛星影像自動化判釋的比對資料,一方面也提供田野調查的參考資訊。
在資料收集方面,本研究自農委會林務局農林航空測量所購得相關資料,資料內容如表3.9及圖3.5、圖3.6所示,同時針對東埔玉山地區的崩塌地判釋工作,程延年等在1987年時受玉管處所委託的《玉山國家公園東埔玉山區地質調查暨解說規劃研究報告》中,已經有具體的成果,因此本研究根據該研究的成果作為基礎,進行九二一前、後後期航空照片之分析與研判,將新增的崩塌地描繪於程延年所繪製的九幅的崩坍分析圖上(原圖亦轉繪自比例尺為一萬分之一的像片基本圖)。
表 3.9 東埔玉山地區九二一地震前、後期航照及像片基本圖資料清單
|
資料名稱 |
比例尺 |
年代 |
數量 |
資料編號 |
|
|
像片基本圖 |
1:10000 |
1980 |
15幅 |
9520-2-18 和社 9520-2-20 郡大山 9520-2-22 阿里山 9520-2-23 同富山 9520-2-24 沙里仙 9520-2-25 觀高 9519-1-02 自忠山 |
9519-1-03 鹿林山 9519-1-04 玉山 9519-1-05 八通關 9519-1-07 霞山 9519-1-08 南玉山西部 9519-1-09 南玉山 9519-1-10 南山 |
|
玉山林區像片基本圖 |
1:5000 |
1986 |
17幅 |
9519-1-004 ~ 008 9519-1-014 ~ 018 |
9519-1-025 ~ 028 9519-1-036 ~ 038 |
|
台大實驗林區像片基本圖 |
1:5000 |
1996 |
13幅 |
9520-2-075 ~ 078 9520-2-085 ~ 088 |
9520-2-094 ~ 098 |
|
航空照片(地震前) |
1996~1997 |
62幅 |
|||
|
航空照片(地震後) |
1999~2000 |
65幅 |
|
|
|
圖3.5 東埔玉山地區九二一地震前、後期航空照片位置圖
圖3.6 東埔玉山地區崩塌分析圖圖幅接合表
根據程延年.,等(1987年)的製圖基礎,加上本研究所判釋九二一地震後新增的崩坍地、國家公園界線及重要標註說明等圖例,將圖例說明如下:
|
|
山嶺線 |
|
裸露地,上面仍未有固結土質(1987年所判釋) |
|
|
新路 |
|
重要崩塌地,提供進一步說明 |
|
|
舊路 |
|
玉山國家公園界線 |
|
|
河流、實點表示河床為現代沈積物充塞 |
|
裸露地,上面仍有未固結土質(2000年所判釋) |
|
|
裸露地,上面仍有未固結土質(1987年所判釋) |
|
裸露地,上面仍未有固結土質(2000年所判釋) |
a.由於研究區內多屬於地形起伏大的區域,需選取影像重疊度達75%以上的相片對,最為航照立體判釋的對象。
b.利用立體鏡並適當調整相片位置,以獲得清晰的立體影像。並開始研判區域性的地形、地質與山坡地崩塌狀況。
c.判釋項目主要以崩塌地為主,並將判釋出的崩塌地範圍轉繪製崩坍分析圖上,如果崩塌裸露地上已有間雜著植生,則可判定為固結土質之崩塌裸露地,若崩塌裸露地上無生長任何植物,且有大量土石堆積狀況,則可判定為上方未固結土質之崩塌裸露地。
d.此為針對於聚落、步道或國家公園內相關設施上、下方所出現的崩塌地,在崩坍分析圖加以註記說明,以利後續更深入探討。
三、崩塌判釋結果
將崩坍分析圖依圖幅分別說明判釋結果如下:
#1.和社幅(9520-II-18)
#2.東埔幅(9520-II-19)
#3.郡大山幅(9520-II-20)
#4.同富山幅(9520-II-23)
#5.沙里仙幅(9520-II-24)
#6.觀高幅(9520-II-25)
#7.鹿林山幅(9519-I-03)
#8.玉山幅(9519-I-04)
#9.八通關幅(9519-I-05)


四、小結
除了且解析度的因素外,透過航照立體對的觀察方式,能夠呈現出立體的地形起伏狀態,有助於崩塌地判釋工作之準確度,因此相對於衛星影像所能判釋出來崩塌地,藉由航空照片判釋的結果,可以找,更多,且面積更小的崩塌地來。
由於時間與經費的限制,僅針對東埔、玉山地區實行航空照片判釋工作,在各區域內,新增的崩塌地有很多約在50~60各左、右,其中以東埔、同富山、玉山及八通關圖幅所增加的面積最大。而新崩塌地多出現於陡峭的坡地以及集水區上游坡地。總結歸納這些新增的崩塌地與幾個原因有關:1. 可能是路基開挖時上方地層坡角被砍除與下方棄土堆積有關;2.地質結構破碎;3. 斷層通過之區域;4.未穩定之舊崩塌地持續崩落;5.坡度陡峻地區;6.屬於溪流集水區上方坡地。
對於這些崩塌地,除了對於園區內步道與相關設施的直接破壞外,崩落的土石以會淤積於河道上或坡面上,造成大雨過後發生土石流或是的土石崩落之危害。因此在進行針園區內災後復原時,需對於新、舊崩塌地加以檢視評估,是否上方邊坡存有未固結的岩塊、以及和河道上是否由土石淤積現象,並針對土石崩解及淤積情況,採取適當的警告或預防措施。