臺灣博碩士論文加值系統

六甲─木屐寮逆衝斷層是台灣西南部沿海平原區至麓山帶的交界斷層之一，在六甲─木屐寮斷層以東的區域有比周圍更高的壓縮應變率，此應變率反映了該區域內較高的應變累積量以及所對應的地震潛勢，此應變累積量可能是淺部斷層所反應或由深部的盲斷層系統所累積，而此區域恰好與1862年台南地震的震央位置有重合，而前人的研究工作也指出1862年台南地震的震源應位於六甲斷層深處，可以看出目前該斷層正處於地震循環中的間震累積能量階段，而1862年的台南地震為上一次的同震能量釋放事件。本研究使用了2002至2016年間之135個GPS移動站、2003至2016年間之57個GPS連續站資料，以Bernese v.5.0解算得到各測站之座標每日解，並以最小二乘法求得各站速度；本研究也使用了地調所2003至2014年間之精密水準成果，另外還使用了由歷史文獻所推估之1862年台南地震的震度分佈資料，從同震及間震兩個時期的資料進行研究，來探討該地震循環的地震潛能與地震週期。本研究中有三個主要的假設，第一為該地震循環中能量釋放之破裂面會與能量累積之斷層面相符，亦即1862年台南地震之破裂面會與現今能量累積之位置相符，也會是下次地震之破裂面，第二為該地震循環中每次地震時釋放之能量都相同，第三為現今斷層面上之能量累積都是由滑移虧損率所貢獻。在這三個假設之下，本研究先設定合理之斷層幾何參數，透過使用經驗式來模擬同震之地表最大加速度值，並參考氣象局之震度分級表繪製為震度圖，再與震度分布資料做比對，以找出最佳發震斷層之斷層深度；本研究同時也使用了由現今之大地測量資料所計算得之地表水平及垂直速度場，採用基線反演模型進行逆推，並設定合理之斷層幾何參數，以求得在研究區域內各個斷層面之滑移虧損率的分佈。由震度圖與震度資料分析的結果顯示，1862年台南地震之發震斷層可能為深部之六甲斷層，破裂面之上緣深度為6公里、下緣深度為16公里，該面積對應之發震規模為MW 6.8；在部分台南地區可以觀察到地表最大加速度之估計值會低於觀測值，而此低估之現象與2016年美濃地震造成的現象很類似，可能是由於地動衰減經驗式低估了台南地區的場址效應所造成。由基線反演模型之成果可看出，於六甲斷層南段深部有一能量累積面，其深度約為9至19公里，與由震度圖與震度資料推估的破裂面相吻合，可呼應本研究之假設，而經計算後得知其平均滑移虧損率約為1.7 mm/yr。在1862年台南地震的破裂面與現今斷層面上能量累積的位置吻合的情況下，若下次地震時的破裂面仍與此面積相同，且能量的累積完全是由平均滑移虧損率所提供，則可計算出地震週期約為677年。

關鍵詞: 六甲─木屐寮斷層、GPS、精密水準、基線反演模型、等震度圖

SUMMARY
The Liuchia–Muchiliao reverse fault is a major seismogenic fault in SW Taiwan, which is also proposed as the source fault of the 1862 Tainan earthquake. In order to characterize the seismic potential in SW Taiwan, I collect the geodetic data during 2002 to 2016 and the intensity distribution of the 1862 Tainan earthquake, then use empirical formulas to simulate the isoseismal maps for optimal fault geometry and corresponding moment magbitude of the source area of 1862 event, and also use baseline inversion model to invert the present–day velocity fields for the distribution of slip deficit rates on each fault plane, finally can estimate the recurrence interval. The results show that the optimal fault depth of the source area is about 6 – 16 km in depth, the corresponding moment magbitude is about MW 6.8, and also show that the current energy accumulating area is similar to the source area, the average slip deficit rate is about 1.7 mm/yr, so the estimated recurrence interval is 677 years. And I also found that there are some PGA underestimates in part of Tainan area, this underestimate may be caused by the ground motion attenuation formula used in this study, and this condition is in good agreement with the phenomenon in 2016 Meinong event.

Key words: Liuchia-Muchiliao fault, GPS, Precise leveling, Baseline inversion model, Isoseismal maps
INTRODUCTION
Taiwan is located at the plate boundary between the Eurasian plate and the Philippine Sea plate, where the convergence rate between two plates is about 82 mm/yr. The high convergence rate results in the generation of fold-and-thrust belt in SW Taiwan. Previous studies based on GPS observations have found that the area to the east of the Liuchia–Muchiliao fault has high contraction strain rate, which means high strain accumulation in the area and the corresponding seismic potential. Interestingly, the location of the epicenter of 1862 Tainan earthquake is covered in the area which has the high strain accumulation, and previous study also indicate that the source area of 1862 Tainan earthquake is in the depth of the Liuchia fault. So I consider that the 1862 event is the previous one coseismic event in the earthquake cycle, and the fault is undergoing interseismic stage now. In order to realize the seismic potential and the recurrence interval, there are three main assumptions in this study:
(1) In the earthquake cycle, the source area which release energy in coseismic event is similar to the area which accumulate energy in interseismic period, which means the source area of 1862 Tainan event is in good agreement with present–day locked area on the fault plane.
(2) In the earthquake cycle, the amount of the release energy is the same, which means the magnitude of 1862 Tainan event will be equal to the magnitude of next event.
(3) The energy accumulation of present–day locked area on the fault plane will fully be attributed to the distribution of slip deficit rates.
Under the assumptions, I collect present–day geodetic data as interseismic data, and intensity distribution of 1862 Tainan earthquake as coseismic data. And then use empirical formulas to simulate the isoseismal maps for optimal fault depth and corresponding moment magbitude of the source area of 1862 event, and also use baseline inversion model to invert the present–day velocity fields for the distribution of slip deficit rates on each fault plane.
MATERIALS ＆ METHODS
I process the present–day geodetic data for the interseismic surface velocity fields. Then set proper fault geometry parameters from previous studies, and use empirical formulas to simulate the peak ground acceleration, next plot the isoseismal maps by comparing with the intensity scale from Central Weather Bureau. After that I can compare the isoseismal maps with the intensity distribution data to find out the optimal fault parameters of the source area. And also invert the velocity fields by using the baseline inversion model, which do Delaunay triangulation to the coordinates of each station for making baselines and calculate the elongation rates of each baseline with vector difference in the baseline direction, then we can get the distribution of slip deficit rates. The fault parameters of the source area can be used to calculate the corresponding moment magnitude, and the slip deficit rates can be further used to estimate the recurrence interval.
RESULTS ＆ DISCUSSION
After comparing the isoseismal maps with the intensity distribution data, I found that the optimal depth of the source area seems to be 6 km in depth of top edge and 16 km in depth of bottom edge (Figure 1). In this figure, the black triangles represent the location where the magnitude is 6, the gray triangles represent the location where the magnitude is 5 and the white triangles represent the location where the magnitude is 4. But even the optimal result still not fit intensity distribution data very well, there still have some observed intensity values are larger than the calculated values in Tainan urban area. This underestimate may be caused by the ground motion attenuation formula used in this research, and this condition is in good agreement with the phenomenon in 2016 Meinong event from previous study. Modeling results show that there is a energy accumulating area which is in the depth of the Liuchia fault as 9 – 19 km, this result is similar to the source area of 1862 Tainan earthquake and also can reflect on the first assumption in this research. So in this study, the baseline inversion model is used to estimate the distribution of slip deficit rates for optimal fault geometry of the source area, the average slip deficit rate of depth as 9 – 19 km is about 1.7 mm/yr.
CONCLUSION
In this study, under the assumptions, I found that:
(1) The result of the isoseismal map analysis shows that the source area is about 6 – 16 km in depth of the Liuchia fault, and the corresponding moment magnitude is about MW 6.8.
(2) The isoseismal map of the optimal source area indicate that there are some PGA observed values larger than calculated values in part of Tainan area, and this underestimate may be caused by the ground motion attenuation formula which is lack of considering the site effect in Tainan area, this condition is in good agreement with the phenomenon in 2016 Meinong event.
(3) The result of the baseline inversion model shows that there is a energy accumulating area in the depth as 9 – 19 km of Liuchia fault, this area is similar to the source area from the result of the isoseismal map analysis. And the average slip deficit rate is about 1.7 mm/yr.
(4) On condition that the source area is similar to the energy accumulating area, which can reflect on the first assumption of this research. If the rupture area of next earthquake event is the same with 1862 event, and the energy accumulation is totally attributed by average slip deficit rate, then the estimated recurrence interval is about 677 years.

摘要 I
Abstract III
誌謝 VII
目錄 VIII
表目錄 X
圖目錄 XI
第一章 緒論 1
1.1 研究動機與目的 1
1.2 研究方法及步驟 2
第二章 研究區域地質概述 7
2.1 斷層 8
2.1.1 觸口斷層 9
2.1.2 木屐寮斷層 10
2.1.3 六甲斷層 10
2.1.4 新化斷層 13
2.1.5 左鎮斷層 15
2.1.6 後甲里斷層 16
2.2 1862年台南地震 17
第三章 測量資料來源及處理 19
3.1 GPS觀測資料 19
3.1.1 GPS觀測資料來源 19
3.1.2 GPS解算策略 20
3.1.3 速度擬合 21
3.2 水準測量資料 22
3.3地表間震速度場 22
3.3.1水平速度場 22
3.3.2垂直速度場 23
第四章 基線應變模型 27
4.1 斷層模型建立方法與參數設定 27
4.1.1 模型理論 27
4.1.2 斷層幾何參數設定 29
4.1.3 滑移虧損率分布粗糙度係數測試 33
4.1.4 斷層淺部之滑移虧損率分布約制測試 33
4.2 模型模擬結果 35
第五章 1862年台南地震震度分析 38
5.1 震度分布資料 38
5.2 震度圖之模擬 38
5.2.1斷層幾何參數設定 38
5.2.2 發震規模 39
5.2.3 最大地表加速度 40
5.3 震度圖模擬結果 41
第六章 討論 44
第七章 結論 48
參考資料 49
附錄一 各測站之時間序列 58
附錄二 模型模擬成果 155
附錄三 震度圖模擬成果 175

王宏霖（2006），以反射震測法研究木屐寮斷層及其附近構造，國立中央大學地球物理所碩士論文，共134頁。
中央氣象局（2000），地震震度分級表，http://www.cwb.gov.tw/V7/earthquake/quake_preparedness.htm，最近檢索時間2017年6月。
中央地質調查所（2013），活動斷層分布圖，http://fault.moeacgs.gov.tw/TaiwanFaults_2009/PageContent.aspx?type=C＆id=5，最近檢索時間2017年7月。
中央地質調查所（2014），斷層活動性觀測研究第三階段─斷層整合性觀測與潛勢分析（2/4）（第二年度），經濟部中央地質調查所執行期末報告書。
中央地質調查所（2015），斷層活動性觀測研究第三階段─斷層整合性觀測與潛勢分析（3/4）（第三年度），經濟部中央地質調查所執行期末報告書。
中央氣象局（2016），地震專欄，http://www.cwb.gov.tw/V7/knowledge/planning/seismological.htm，最近檢索時間2017年6月。
中國石油公司台灣油礦探勘總處（1989），十萬分之一地質圖幅No.6─台南。
水利署（2016），地層下陷狀況，http://www.wra.gov.tw/ct.asp?xItem=45261＆ctNode=9915＆comefrom=lp#9915，最近檢索時間2017年7月。
石再添、鄧國雄、張瑞津、石慶得、楊貴三（1986），臺灣活斷層的地形學研究，師大地理研究報告，第12期，第1–44頁。
石瑞銓、王乾盈、謝昭輝（1998），雲嘉南地區活動斷層震測調查，經濟部中央地質調查所活動斷層地球物理探勘計畫八十七年度報告，共130頁。
吉田要（1931），所謂觸口斷層，臺灣地學記事，第2號，第27–29頁。
辛在勤（1998），台灣地區地震預警之初探訪，氣象學報42卷2期，第118–134頁。
李錫堤、鄭錦桐、紀立民、廖啟雯、廖卿妃（2000），槽溝開挖探查成果簡介，第三屆台灣活動斷層調查研討會資料。
李錫堤 （2002），後甲里斷層槽溝探查工作成果報告，國立中央大學應用地質研究所，共40頁。
李秋賢（2016），利用GPS觀測資料及塊體模型來探討台灣的地殼變形，國立中央大學碩士論文，共170頁。
何春蓀（1986），台灣地質概論，經濟部中央地質調查所，共118頁。
邱朝聰（1997），台灣西南地區之震測研究─西南海域之大陸棚及新化斷層為例，國立臺灣海洋大學應用地球物理研究所碩士論文，共47頁。
林朝棨（1957），臺灣地形，臺灣省文獻委員會，共424頁。
林啟文、張徽正、盧詩丁、石同生、黃文正（2000a），台灣活動斷層概論：五十萬分之一台灣活動斷層分析圖說明書，經濟部中央地質調查所，共121頁。
林啟文、石同生、盧詩丁（2000b），左鎮斷層，經濟部中央地質調查所活動斷層調查報告。
林燕慧、劉彥求、李明書、林偉雄、林啟文（2004），後甲里斷層，經濟部中央地質調查所活動斷層調查報告。
林耕霈（2012），利用永久性散射體差分干涉法探討台南地區之地殼形變，國立中央大學地球物理研究所碩士論文，共81頁。
耿文溥（1986），台灣中部竹山及嘉義間之地質，經濟部中央地質調查所彙刊，第1–26頁。
徐明同（1980），台灣地震目錄（西元1644年至1979年），國立台灣大學地震工程研究中心，CEER S69-01，共77頁。
徐明同（1983），明清時代破壞性大地震規模及震度之評估，明清時代破壞性大地震規模及震度之評估，氣象學報第 29 卷第四期，第1–18頁。
孫習之（1970），台灣省新營至嘉義平原區域航照地質之研究，台灣石油地質，第7號，第133–144頁。
孫習之（1971），台灣省新營至嘉義平原區域航照地質之研究，台灣石油地質，第8號，第65–76頁。
陳文山、李錫堤、陳于高（2002），槽溝開挖與古地震研究計畫（1/5）（上），經濟部中央地質調查所地震地質調查與活動斷層資料庫建置計畫研究報告，共107頁。
陳文山、李錫堤、石瑞銓、楊小青、楊志成、顏一勤、劉力豪、張徽正、侯進雄（2004），新化斷層的構造特性與古地震研究，經濟部中央地質調查所特刊，第15號，第111–119頁。
陳文山、陳于高、楊小青（2006），地震地質調查及活動斷層資料庫建置─槽溝開挖與古地震研究計畫（5/5）與總報告，經濟部中央地質調查所報告95–08號，共133頁。
陳心怡（2006），台灣南部地殼變形分析：1995–2005年GPS觀測，國立成功大學地球科學系碩士論文，共144頁。
陳延宗（2007），以反射震測法研究左鎮斷層及其附近構造，國立中央大學地球物理研究所碩士論文，共90頁。
陳建宏（2007），以反射震測法研究六甲斷層及其附近構造，國立中央大學地球物理研究所碩士論文，共80頁。
陳盈蒨（2016），以2016年美濃地震為例探討如何適當地應用TELES進行台南地區地震災害評估，成功大學地球科學系碩士在職專班學位論文，共144頁。
郭炫佑 （1999），後甲里斷層及其附近構造，國立中央大學地球物理研究所碩士論文，共83頁。
郭素妙（2008），利用GPS與地震資料探討台灣西南部地震潛能，國立成功大學地球科學系碩士論文，共66頁。
張徽正、林啟文、盧詩丁、陳勉銘（1998），臺灣活動斷層概論：五十萬分之一臺灣活動斷層分布圖說明書，經濟部中央地質調查所，共103頁。
國土測繪中心（2012），大地基準及一九九七座標系統2010年成果工作總報告，共137頁。
張麗旭、周敏、陳培源（1947），民國35年12月5日台南之地震，台灣省地質調查所彙刊，第1號，第11–18頁。
張徽正、林啟文、陳勉銘、盧詩丁（1998），台灣活動斷層分布圖說明書，經濟部中央地質調查所特刊，第十號，共103頁。
張憲卿（2008），台灣地質圖說明書，嘉義圖幅（五萬分之一），經濟部中央地質調查所。
許雅儒（2004），集集地震之震前、同震及震後變形模式研究，國立中央大學地球物理研究所博士論文，共133頁。
許京穎（2008），台灣潛勢地震之發生機率估計，國立中央大學地球物理研究所碩士論文，共133頁。
黃鑑水、張憲卿、劉桓吉（1994），台灣南部觸口斷層之調查與探勘，經濟部中央地質調查所彙刊，第9號，第51–76頁。
黃旭燦（2003），台灣中南部褶皺逆衝斷層帶地質構造特徵分析，國立中央大學地球物理研究所博士論文，共130頁。
黃旭燦、楊耿明、吳榮章、丁信修、李長之、梅文威、徐祥宏（2004），斷層活動性觀測與地震潛勢評估調查研究，台灣陸上斷層帶地質構造與地殼變形調查研究（5/5）─台灣西部麓山帶地區地下構造綜合分析，經濟部中央地質調查所報告93–13號，共59頁。
黃勝群、陳彥宏、石瑞銓、張徽正、林啟文、陳文山（2004），新化斷層帶之淺部構造特徵，經濟部中央地質調查所特刊，第15號，第99–110頁。
楊耿明、洪日豪、吳榮章、黃旭燦（2001），斷層活動性觀測與地震潛勢評估調查研究：台灣陸上斷層帶地質構造與地殼變形調查研究（1/5）─西南地區（觸口斷層），經濟部中央地質調查所報告，第89–4號，共93頁。
楊志成、顏一勤、宋時驊、黃能偉、陳勇全、陳文山、陳于高、吳樂群、張徽正、侯進雄、林啓文（2005），六甲斷層近萬年來滑移速率之探討，經濟部中央地質調查所特刊，第16號，第1–16頁。
劉俊忠（2005），六甲斷層之淺層反射震測與野外地質分析，國立中正大學應用地球物理研究所碩士論文，共74頁。
鄭世楠、葉永田（1989），西元1604年至1988年台灣地區地震目錄，中央研究院地球科學研究所，IES-R-661，共255頁。
鄭世楠、葉永田、徐明同、辛在勤（1999），台灣十大災害地震圖集，中央氣象局及中央研究院地球科學研究所，共289頁。
鄭世楠、王子賓、林祖慰、江嘉豪（2010），台灣地區地震目錄的建置，中央氣象局地震技術報告彙編，第54卷，第575–605頁。
鄭世楠、江嘉豪、陳燕玲（2012），台灣地區歷史地震資料的建置，交通部中央氣象局地震技術報告彙編，第60卷，第427–448頁。
鄭世楠、邵承芬（2013），1862年臺南地震的回顧與探討，2013土木工程與防災研討會，第1–4頁。
鄭世楠 （2014），歷史地震第二講：1862年臺南地震，台灣地震科學中心─歷史地震典藏及推廣，共51頁。
盧詩丁、石同生、林啟文（2000a），木屐寮斷層，經濟部中央地質調查所活動斷層調查報告。
盧詩丁、石同生、林啟文（2000b），六甲斷層，經濟部中央地質調查所活動斷層調查報告。
饒瑞鈞、景國恩、謝宗訓、余致義、侯進雄、李元希、胡植慶、詹瑜璋、李建成、洪日豪（2003），台南台地的地表變形與地震潛能，經濟部中央地質調查所特刊，第14號，第147–156頁。
饒瑞鈞（2014），GPS觀測之雜訊分析研究II，交通部中央氣象局地震測報中心。
饒瑞鈞、景國恩（2016），利用大地測量資料探討台灣西南部之地震潛能，交通部中央氣象局委託研究計畫期中成果報告，共29頁。
Bonilla, M. G.（1975）, A review of recently active faults in Taiwan, US Geological Survey Open-File Report, No. 75–41, 58.
Ching, K. E., M. L. Hsieh, K. M. Johnson, K. H. Chen, R. J.Rau, and M. Yang (2011), Modern vertical deformation rates and mountain building in Taiwan from precise leveling and continuous GPS observations, 2000-2008, J. Geophys. Res., 116(B8).
Dach, R., U. Hugentobler, P. Fridez, and M. Meindl (Eds.) (2007), Bernese GPS Software Version 5.0, Astronomical Institute, University of Berne, 612.
Harris, R. A., and P. Segall, (1987), Detection of a locked zone at depth on the Parkfield, California, segment of the San Andreas fault, Journal of Geophysical Research: Solid Earth, 92(B8), 7945-7962.
Ho, C. S. (1976), Foothills tectonics of Taiwan, Bull. Geol. Surv. Taiwan, 25, 9–28.
Huang, S. T., K. M. Yang, J. H. Hung, J. C. Wu, H. H. Ting, W. W. Mei, S. H. Hsu, and M. Lee (2004), Deformation front development at the northeast margin of the Tainan basin, Tainan–Kaohsiung area, Taiwan, Marine Geophysical Research, 25(1), 139–156.
Huang, M. H., J. C. Hu, K. E. Ching, R. J. Rau, C. S. Hsieh, E. Pathier, B. Fruneau, and B. Deffontaines (2009), Active deformation of Tainan tableland of southwestern Taiwan based on geodetic measurements and SAR interferometry, Tectonophysics, 466(3–4), 322–33410.
Hung, J. H. (1999), Structure and motion of the southwestern Taiwan fold and thrust belt, TAO, Vol.10, No. 3, 543–568.
Jamison, W. R. (1993), Mechanical stability of the triangle zone: the backthrust wedge, Journal of Geophysical Research: Solid Earth, 98(B11), 20015–20030.
Jean, W. Y., Y. W. Chang, K. L. Wen, and C. H. Loh (2006), Early estimation of seismic hazard for strong earthquakes in Taiwan, Natural hazards, 37(1–2), 39–53.
Jean, W. Y., and C. H. Loh (2001), A study on the classification of site effects and its application to the seismic hazard and microzonation, In 10th International Conference on Soil Dynamics and Earthquake Engineering, SDEE, 82.
Johnson, K. H. (2013), High Moment Accumulation Rates on Southern California Faults Inferred from GPS Data, Indiana University, personal communication, 1–21.
Lacombe, O., F. Mouthereau, B. Deffontaines, J. Angelier, H. T. Chu, C. T. Lee (1999), Geometry and quatemary kinematics of fold-and-thrust units of southwestern Taiwan: Tectonics, 18, 1198–1223.
Menke, W. (1989), Geophysical data analysis: Discrete Inverse Theory, 2nd Edition, Academic press, New York, 289.
Okada, Y. (1985), Surface deformation due to shear and tensile faults in a half-space, Bulletin of the seismological society of America, 75(4), 1135–1154.
Okada, Y. (1992), Internal deformation due to shear and tensile faults in a half-space, Bulletin of the Seismological Society of America, 82(2), 1018–1040.
Price, N. J., and J. W. Cosgrove (1990), Analysis of Geological Structures, Cambridge Univ. Press, Cambridge, 502.
Rothacher, M., and L. Mervart (1996), Bernese GPS Sofrware version 4.0, Astronomical Institude, University of Berne, 418.
Savage, J. C. (1983). A dislocation model of strain accumulation and release at a subduction zone, Journal of Geophysical Research: Solid Earth, 88(B6), 4984–4996.
Shyu, J. B. H., Y. R. Chuang, Y. L. Chen, Y. R. Lee, and C. T. Cheng (2016), A New On-Land Seismogenic Structure Source Database from the Taiwan Earthquake Model (TEM) Project for Seismic Hazard Analysis of Taiwan, Terrestrial, Atmospheric ＆ Oceanic Sciences, 27(3).
Sun, S. C. (1964), Photogeologic study of the Tainan-Kaohsiung Coastal plain area, Taiwan, Petrol. Geol. Taiwan, 3, 39–51.
Sun, S. C. (1971), Photogeologic study of Hsinying-Chiayi coastal plain, Taiwan, Petrol. Geol. Taiwan, 65–75.
Suppe, J. (1980), Imbricated structure of western foothills belt, south-central Taiwan, Petrol. Geol. Taiwan, 17, 1–16.
Tsai, Y. B. (1985), A study of disastrous earthquakes in Taiwan, 1683-1895, Bull. Inst. Earth Sc., Acad. Sin., 5, 1–44.
Wang, C. Y., G. P. Chen, and D. T. Jong (1994), The detection of active faults on Taiwan using shallow reflection seismics, Terr. Atmos. Ocean. Sci., 5, 277–293.
Wang, S. (1976), ERTS–1 satellite imagery and its application in regional geologic study of southwestern Taiwan, Petrol. Geol. Taiwan, 13, 37–57.
Wdowinski, S., Y. Bock, J. Zhang, P. Fang, and J. Genrich (1997), Southern California permanent GPS geodetic array: Spatial filtering of daily positions for estimating coseismic and postseismic displacements induced by the 1992 Landers earthquake, Journal of Geophysical Research: Solid Earth, 102(B8), 18057–18070.
Wells, D. L., and K. J. Coppersmith (1994), New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement, Bulletin of the seismological Society of America, 84(4), 974–1002.
Yang, K. M., S. T. Huang, J. C. Wu, H. H. Ting, W. W. Mei, and J. B. Wang (2010), Structural Transfer Zones in the Foothills Belt, Western Taiwan, European Geosciences Union, General Assembly 2010, 10776.
Youngs, R. R., and K. J. Coppersmith (1985), Implications of fault slip rates and earthquake recurrence models to probabilistic seismic hazard estimates, Bulletin of the Seismological society of America, 75(4), 939–964.
Yu, S. B., H. Y. Chen, and L. C. Kuo (1997), Velocity field of GPS stations in the Taiwan area, Tectonophysics, 274, 41–59.