2008年3月4日17時31分47秒(UT)在高雄市桃源區發生芮氏規模5.2之地震。該地震伴隨著前震與大量之餘震,且餘震主要集中發生於主震後8小時之內,為一個完整之地震序列。地震序列發生的時間適逢為TAIGER 計畫執行期間,該計畫在全台各地設置大量的地震儀以記錄人工爆震源與天然地震。因此該地震序列之波形資料,被鄰近密佈之地震觀測儀完整地記錄下來,提供了多方向研究之可能。
首先,本研究利用近震源之高密度地震網,對此地震序列進行完整之初步定位,並使用兩階段重新定位,第一階段先引入測站修正值重新定位,提高其相對定位精準度,第二階段使用hypoDD進行第二次重新定位。由最後之定位結果,得以了解該地震斷層之形貌,確定其斷層走向為東北-西南方向,傾角約為 45°往東南方向傾斜,且地震深度隨著地震走向往西南方向逐漸變深。同時此地震序列之空間分布位置與形貌,以及主震之震源機制,皆與鄰近之內英斷層近似,因此推論此地震序列應與內英斷層有關。
本研究進一步再利用tomoDD逆推此地震序列發生之區域,以及鄰近地區之三維速度構造。在Vp與Vs之速度變化以及Vp/Vs分布上,可以看出地震多分佈在高值與低值區域之交界處,而主震發生區域緊臨一個相對低速區域。推斷此區域可能為破裂面,或為斷層滑動之界面所在。同時在地震分布並延伸至地表的方向上,存在一個傾角為45°的低速帶,且此低速帶與地表之內英斷層的位置相當。因此由速度構造中低速帶之存在、震源之空間分布及內英斷層的地表位置,三者間在空間中之關係可佐證內英斷層與地震序列之間的相關性。同時在鄰近區域之速度變化上,西邊之速度也較東邊之速度低,代表地層由更新世的新沉積層,轉變為中央山脈近地表
之第三紀變質岩帶。
最後,透過使用時間項之方法,利用TAIGER計畫之人工震源資料並聯合天然地震資料,讀取台灣北部測站之 Pn 波到時,推求北台灣下方之深部構造變化之時間項變化。時間項之分佈特徵為,台灣北部至東北部存在一個時間項之相對低區,在台北盆地、宜蘭清水以及宜蘭盆地,皆表現出時間項之相對高區。由透過不同之折射層P波速度、地殼P波速度以及折射層深度的組合,所得到之理論到時與觀測到時的殘差值為最小時,可得到合理的折射層P波速度、地殼P波波速以及地殼厚度組合為 7.6 km/sec、5.6 km/sec 以及35公里。但沉積層或地殼中之異常的速度構造,將對時間項產生影響,使得地殼厚度之估算造成誤差。因此,在東西向之剖線上進行沉積層修正,所求得之莫荷面較淺且變化較為平緩,但變化的趨勢與前
人研究一致。其變化為西部平原與麓山帶約30公里,進入中央山脈後深度
加深至約34公里,再往東部深度再度變淺為約28公里。

The March 4, 2008, a moderate earthquake (ML 5.2) occurred in the southwestern part of Taoyuan District of Kaohsiung City in southern Taiwan. It was followed by numerous aftershocks in the next 48 hours, including three events with magnitude lager than 4. This sequence occurred during the TAIGER (TAiwan Integrated GEodynamics Research) controlled-source seismic experiment, several seismic networks were deployed in the Taiwan area at this time
and many stations recorded this earthquake sequence in the near-source region.
We processed the near-source observations and determined 2,340 events in preliminary result. Then we relocated events with station corrections first and relocated again by applying the double-difference relocation algorithm with the
differential arrival times obtained by the waveform cross-correlation. The resulting aftershocks are extended along the NE-SW direction and located on a 45° SE-dipping plane. This plane agrees to one of the nodal planes of the focal
mechanism determined by BATS (strike=N 37°E, dip=48°SE and rake=96°).
Furthermore, we used double-difference tomography to relocate earthquake locations and determine P and S wave velocity structure around the sequence of Taoyuan earthquake. A total of 3,034 events were selected from preliminary results and CWBSN catalog. The relocation results are similar to the prior results and the 3D velocity structures show that obvious velocity variations exist in the crust beneath the Ne-In fault. The sequence occurred on the interface between the low and high velocity area. Based on the hypocentral distribution pattern, focal mechanism and velocity structure variations, we suggest that the Taoyuan earthquake was probably caused by Ne-In fault.
The P arrivals of 431 ray paths with the epicentral distances longer than 160 km from 10 explosive sources and 10 earthquakes in Taiwan were analyzed according to the time term method, a simple travel time inversion for seismic re-
fraction crustal studies. We used the data to investigate the variation of time term values in northern Taiwan. The results show that relatively high time term values are obtained in the areas of Taipei basin and Chingshui area. The relative high values locate at the northernmost and northeastern area in Taiwan. We also determine the optimal combinations of Pn velocity, crust velocity and Moho depth. The values of 7.6 km/s, 5.6 km/s and 35 km are obtained, respectively.

Ai Y., Q. F. Chen, F. Zeng, X. Hong, and W. Ye, The crust and upper mantle structure beneath southeastern China, Earth Planet. Sci. Lett. , 260, 549-563, 2007.
Aki, K. and W. H. K. Lee, Determination of three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes: 1. A homogeneous initial model, J. Geophys. Res. , 81, 4381-4399, 1976.
Angerlier, J., Geodynamics of the Eurasia-Philippine Sea Plate boundary: Preface, Tectonophys. , 125, IX-X, 1986.
Berry, M. J. and G. F. West, An interpretation of the first-arrival data of the Lake Superior experiment by the time-term method, Bull. Seism. Soc. Am. , 56, 141-171, 1966.
Biq, C., Dual-trench structure in the Taiwan-Luzon region, J. Geol. Soc. China , 15, 65-75, 1972.
Castagna, J. P., M. L. Batzle, and R. L. Eastwood, Reltionships between compressional-wave and shear-wave velocities in clastic silicate rocks, Geophysics , 50, 571-581, 1985.
Chemenda, A. I., R. K. Yang, C. H. Hsieh, and A. L. Groholsky, Evolutionary model for the Taiwan collision based on physical modeling, Tectonophys. , 274, 253-274, 1997.
Chen, C. H., Y. H. Chen, H. Y. Yen, and G. K. Yu, Lateral variations of Pn velocity and anisotropy in Taiwan from travel-time tomography, Earth, Planets and Space , 55, 223-230, 2003.
Davis, D., J. Suppe, and F. A. Dahlen, Mechanics of fold-and-thrust belts and accretionary wedges, J. Geophys. Res. , 88, 1153-1172, 1983.
Dewey, J. W., Seismicity studies with the method of joint hypocenter determination, Ph. D. Thesis, University of California at Berkeley, 166, 1971.
Freedman, H. W., Seismological measurements and measurement error, Bull. Seism. Soc. Am. , 58, 1261-1271, 1968.
Frohlich, C., An efficient method for joint hypocenter determination for large groups of earthquakes, Comp. Geosci. , 5, 387-389, 1979.
Got, J. L., J. Frèchet, and F. W. Klein, Deep fault plane geometry inferred from multiplet relative relocation beneath the south flank of Kilauea, J. Geophys. Res. , 99, 15375-15386, 1994.
Ho, C. S., An introduction to the geology of Taiwan explanatory text of the geologic map of Taiwan, Bulletin of the Central Geological Survey , Taiwan, 2006.
Hsu, S. K. and J. C. Sibuet, Is Taiwan the result of arc-continent or arc-arc collision?, Earth Planet. Sci. Lett. , 136, 315-324, 1995.
Huang, B. S., K. C. Chen, and H. Y. Yen, Velocity of Pn-waves in the Taiwan Strait and its surrounding area from regional earthquakes, Terr. Atmos. Ocean. Sci. , 9, 473-486, 1998.
Kao, H. and R. J. Rau, Detail structures of the subducted Philippine Sea plate beneath northeast Taiwan: A new type of double seismic zone, J.Geophys. Res. , 104, 1015-1033, 1999.
Kim K. H., J. M. Chiu, H. Kao, Q. Liu, and Y. H. Yeh, A preliminary study of crustal structure in Taiwan region using receiver function analysis, Geophys. Int. J. , 159, 146-164, 2004.
Kim K. H., J. M. Chiu, J. Pujol, K. C. Chen., B. S. Huang., Y. H. Yeh and P. Shen, Three-dimensional Vp and Vs structure models associated with the active subduction and collision tectonics in the Taiwan region, Geophys. Int. J. , 162, 204-220, 2005.
Kuo-Chen H., Francis T. Wu, and S. W. Roecker, Three-dimensional P velocity structures of the lithosphere beneath Taiwan from the analysis of TAIGER and related seismic data sets, J. Geophys. Res. , 117, B06306, 2012.
Lee, C. T., C. T. Cheng, C. W. Liao, and Y.B. Tsai, Site classification of Taiwan free-field strong-motion stations, Bull. Seism. Soc. Am. , 91, 1283-1297, 2001.
Lee, C. T., and B. R. Tsai, Mapping V S30 in Taiwan, Terr. Atmos. Ocean. Sci. , 19, 671-682, 2008.
Lin, A. T., A. B. Watts and S. P. Hesselbo, Cenozoic stratigraphy and subsidence history of the South China Sea margin in the Taiwan region, Basin Research , 15, 453-478, 2003.
Lin, C. H., Y. H. Yeh, H. Y. Yen, K. C. Chen, B. S. Huang, S. W. Roecker, and J. M. Chiu, Three-dimensional elastic wave velocity structures of the Hualien region of Taiwan : Evidence of active crustal exhumation, Tectonics , 17, 89-103, 1998.
Lin C. H. and Y. H. Yeh, Three-dimensional P- and S-wave velocity structures of the Chingshui-Tuchang geothermal area in northeastern Taiwan, Western Pacific Earth Sciences, 1, 73-84, 2001.
Lin, C. H., Active continental subduction and exhumation: Taiwan orogeny, Terra Nova , 14, 281-287, 2002.
Lin, C.H., Temporal b-value variations throughout seismic faulting process: the 2008 Taoyuan earthquake in Taiwan, Terr. Atmos. Ocean. Sci. , 21, 229-234, 2010.
Lin, Y. P., L. Zhao, and S. H. Hung., First-arrival times from an active-source experiment in the TAIGER project and their application to the assessment of tomography models, Bull. Seism. Soc. Am. , 101, 866-880, 2011.
Ludwig, W. J., N. Kumar, and R. E. Houtz, Profiler-sonobuoy
measurements in the South China Sea Basin, J. Geophys. Res. , 84, 3505-3518, 1979.
Ma, K. F., J. H. Wang, and D. Zhao, Three-dimensional seismic velocity structure of the crust and uppermost mantle beneath Taiwan, J. of Phy. Earth , 44, 85-105, 1996.
Malavieille, J., S. Lallemand, S. Dominguez, A. Deschamps, C. Y. Lu, C. S. Liu, P. Schnqrle, and the ACT scientific crew, Arc-continent collision in Taiwan: new marine observations and tectonic evolution, Geol. Soc. Am. , Spec. Pap. 358, 187-211, 2002.
McCollon, R. L. and R. S. Crosson, An array study of upper mantle velocity in Washington State, Bull. Seism. Soc. Am. , 65, 467-482, 1975.
Murdock, J. N. and L. H. Jaksha, Time-term solutions and corresponding data for the crustal structure of north-central New Mexico, U.S. Geological Survey, Open-file Report 80-2014, 39pp, 1980.
Pavlis, G. L., F. Vernon, D. Harvey, and D. Quinlan, The generalized earthquake location (GENLOC) package: an earthquake location library, Comp. Geosci. , 30, 1079-1091, 2004.
Pujol, J., Comments on the joint determination of hypocenters and station corrections, Bull. Seism. Soc. Am. , 78, 1179-1189, 1988.
Pujol, J., Application of the JHD technique to the Loma Prieta, California, mainshock-aftershock sequence and implications for earthquake location, Bull. Seism. Soc. Am ., 85, 129-150, 1995.
Rau, R. J., and F. T. Wu, Tomographic imaging of lithospheric structures under Taiwan, Earth Planet. Sci. Lett. , 133, 517-532, 1995.
Schaff, D. P., G. H. R. Bokelmann, W. L. Ellsworth, E. Zanzerkia, F. Waldhauser, and G. C. Beroza, Optimizing correlation techniques for improved earthquake location, Bull. Seism. Soc. Am. , 94, 705-721, 2004.
Scheidegger, A. E. and P. L. Willmore, The use of a least squares method for the interpretation of data from seismic surveys, Geophysics , 22, 9-22, 1957.
Shearer, P. M., Introduction to seismology , Cambridge University Press. 260pp, 1999.
Suppe, J., Imbricated structure of western Foothills belt,
southcentral Taiwan, Petrol. Geol. Taiwan , 17, 1-16, 1980.
Suppe, J., Mechanics of mountain building and metamorphism in Taiwan, Geological Society of China , 4, 67-89, 1981.
Teng, L. S., Geotectonic evolution of late Cenozoic arc-collision in Taiwan. Tectonophys. , 183, 57-76, 1990.
Teng, L. S. and C. T. Lee, and Y. B. Tsai, An integrated lithospheric modal of Taiwan: Proc. 7 th Taiwan Symposium in Geophy., 519-527, 1998.
Teng, L. S., C. T. Lee, Y. B. Tsai, and L. Y. Hsiao, Slab breakoff as a mechanism for flipping of subduction polarity in Taiwan, Geology , 28, 155-158, 2000.
Thurber C., S. W. Roecker, W. Ellsworth, Y. Chen, W. Lutter, and R. Sessions, Two-dimensional seismic image of the San Andreas Fault in the Northern Gabilan Range, central California: evidence for fluids in the fault zone, Geophys. Res. Lett. , 24, 1591-1594, 1997.
Toomey, D. R. and G. R. Foulger, Application of tomographic inversion to local earthquake data from the Hengill-Grensdalur central volcano complex, Iceland, J. Geophys. Res. , 94, 17497-17510, 1989.
Tsan, S. F. and W. P. Keng, The Neogene rocks and major structural features of southwestern Taiwan: Proc. Geol. Soc. China ., 11, 45-4, 1968.
Walck, M. C., Three-dimensional Vp/Vs variations for the Coso region, California, J. geophys. Res. , 93, 2047-2052, 1988.
Waldhauser, F., and W. L. Ellsworth, A double-difference earthquake location algorithm: method and application to the northern Hayward fault, California, Bull. Seism. Soc. Am. , 90, 1353-1368, 2000.
Wang C. Y., Y. H. Lee, M. L. Ger, and Y. L. Chen, Investing Subsurface Structure and P- and S-wave Velocities in the Taipei Basin, Terres. Atmos. Ocean. Sci. , 15, 609-627, 2004.
Willmore, P. L. and A. M. Bancroft, The time term approach to refraction seismology, Geophys. J. Res. , 3419-3432, 1960.
Wu, F. T., R. J. Rau, and D. Salzberg, Taiwan orogeny: thin-skinned or lithospheric collision?, Tectonophys. , 274, 191-220, 1997.
Wu, Y. M., C. H. Chang, L. Zhao, J. Bruce H. Shyu, Y. G. Chen, Kerry Sieh, and J. P. Avouac, Seismic tomography of Taiwan: Improved constraints from a dense network of string motion stations, J. Geophy. Res. , 112, B08312, 2007.
Yeh, Y. H. and Y. B. Tsai, Crustal structure of central Taiwan from inversion of P-wave arrival times, Bull. Inst. Earth Sci., Academia Sinica , 1, 83-102, 1981.
Yeh, Y. H., E. Barrier, C. H. Lin, and J. Angelier, Stress tensor analysis in the Taiwan area from focal mechanisms of earthquakes, Tectonophys. , 200, 267-280, 1991.
Yeh, Y. H., R. C. Shih, C. H. Lin, C. C. Liu, H. Y. Yen, B. S. Huang, C. S. Liu, P. Z. Chen, C. S. Huang, C. J. Wu, and F. T. Wu, Onshore/offshore wide-angle deep seismic profiling in Taiwan, Terres. Atmos. Ocean. Sci. , 9, 301–316, 1998.
Yen, H. Y. and Y. H. Yeh, Two-dimensional crust structure of Taiwan from gravity data, Tectonics , 17, 104-111, 1998.
Zhang, H. and C. H. Thurber, Double-difference tomography: The method and its application to the Hayward fault, California, Bull. Seism. Soc. Am. , 93, 1875-1889, 2003.
大井上義近、千谷好之助、飯塚保五郎、佐藤戈止和村山賢一, 台灣油田地質概查圖及報告 ,比例尺 1:200,000,日本海軍省,第 115-146 頁, 1928 年。
林啟文,張徽正,盧詩丁,石同生,黃文正, 台灣活動斷層概論 ,第二版, 經濟部中央地質調查所,共 122 頁,2000 年。
林啟文、陳文山、劉彥求與陳柏村, 中央地質調查所特刊第 23 號:台灣東部與南部的活動斷層 ,共 178 頁,2009 年。
何美儀, 台灣西部地區三維速度構造 ,國立中央大學地球物理研究所碩士論文,共 108 頁,1994 年。
何春蓀, 台灣地質概論:台灣地質圖說明書 ,第二版,經濟部中央地質調查所,共 164 頁,1986 年。
陳文山、陳于高、楊小青, 地震地質調查及活動斷層資料庫建置-槽溝開挖與古地震研究計畫(4/5) ,經濟部中央地質調查所研究報告 94-7 號,共 165 頁,2005 年。
陳燕玲, 台灣地區三維速度構造與隱沒構造之相關探討 ,國立中央大學地球物理研究所,碩士論文,1995 年。