期末報告二：花蓮奇萊平原附近地表起伏與地震剖面等資料

由GMT與Matlab程式應用，來看花蓮奇萊平原地區歷年(1900年~2009年)地震分布在本課程學到不少。另外如使用儀器（3D Scanner）或相關應用軟體（Mapper、Surfer等）描繪課程，與在本課程中以GMT使用twdtm_small.grd檔案繪製立體地形顯示，此grd檔案所做出來的地形效果應僅粗淺顯示地表起伏，而不為高解析度地形分析所可使用，需有更好的grd檔為佳。

以2008年地震資料作出奇萊平原地表剖面起伏與地震剖面來看(如下圖)，顯示地表剖面穿越中央山脈並橫越木瓜溪河谷至奇萊平原(應有穿過米崙台地，惟無起伏變化可見)。地震剖面則集中分布於淺成地震深度範圍，顯示板塊交界處的擠壓作用產生近地表岩層破裂的機制十分頻繁，以下列出下圖Script：

project -C121.25/24 -E121.85/24 -G0.1 -N -Q > temp1.d
gawk "{print $1, $2}" temp1.d > temp2.d
grdtrack temp2.d -Gtwdtm_small.grd -V > temp3.d
gawk "{print $1, $3}" temp3.d > temp4.d
makecpt -Crelief -T0/3000/100 -V -Z > tw_test.cpt
makecpt -Crainbow -T0/400/40 -V -Z > depth.cpt
grdgradient twdtm_small.grd -Gtwdtm_small_temp_shade.grd -A60 -Ne0.9 -V
psbasemap -JM4i -R121/121.85/23.1/24.5 -Ba1f1NwsE -V -K > tw_topo.ps
grdimage twdtm_small.grd -Itwdtm_small_temp_shade.grd -JM4i -R121/121.85/23.1/24.5 -Ctw_test.cpt -V -K -O >> tw_topo.ps
psxy -JM4i -R121/121.85/23.1/24.5 eq_tw_formatted02.txt -Sc -Cdepth.cpt -: -V -K -O >> tw_topo.ps
psxy temp2.d -JM4i -R121/121.85/23.1/24.5 -W10/250 -V -K -O >> tw_topo.ps
pscoast -JM4i -R121/121.85/23.1/24.5 -Dh -W3 -Na -Ia -T121.15/23.85/0.5 -L121.15/23.75/23.75/10k -V -K -O >> tw_topo.ps
psscale -D0.4c/5c/4c/0.25c -Cdepth.cpt -B50 -V -K -O >> tw_topo.ps
psxy temp4.d -JX3/2i -R121.25/121.85/-5000/3000 -Sc0.01 -B1/1000nwSE -X5i -Y3i -V -K -O >> tw_topo.ps
psxy eq1.d -JX3/2i -R121.25/121.85/-100/0 -Sc -B1/25nwSE -X0 -Y-3i -V -O >> tw_topo.ps
pause

rem grd2cpt twdtm_small.grd -Crelief -R121.5/122.5/23.5/24.5 -V -Z > tw_test.cpt



如果將地表地震剖面範圍擴大(如下圖)，則隨著地表變化震源深度也會增加，這是靠近張裂海槽的板塊構造活動特性。(圖示地震規模由ML2~ML6)



以下則為3D震源分佈圖：



以震源機制言，擷取奇萊平原地區在1976年~2008年Havard CMT Project Web資料，可見張應力與壓應力分布與斷層滑動面的符合狀況，並可見奇萊平原有逆斷層帶平移破裂現象，縱谷東側則多為逆斷層破裂，中央山脈東翼則有正斷層破裂機制。
set range=119.5/122.5/21.5/25.5
psbasemap -JM6i -R%range% -P -B0.5nwSE -V -K > focal_mech.ps
psmeca EQ_CMT_Harvard_Over6_Taiwan.dat -JM6i -R%range% -Sd1c/9 -H -V -K -O >> focal_mech.ps
pscoast -JM6i -R%range% -Df -W1 -O >> focal_mech.ps
ps2pdf focal_mech.ps
pause

台灣東部地區2008年地震震源機制圖

在 Windows XP 下裝設Generic Mapping Tools