本研究使用複數解載(Complex demodulation)、波線追跡法(Ray tracing)與三維衰減構造(Attenuation structure)等方法，以中央氣象局強地動觀測計畫(Taiwan Strong Motion Instrumental Program ,TSMIP )臺中地區、臺南地區、宜蘭地區、花蓮地區與臺東地區測站地震事件記錄，解算震源P波加速度時間函數，為研究震源P波加速度時間函數之新解法。
本研究目標為震源P波加速度時間函數，選取芮氏規模(ML)大於6.0的地震事件，如此震波傳遞至測站方有足夠之初達P波加速度波形紀錄可供分析；並且各地區測站地震紀錄之地震強度大於4級，如此初達P波波形方足夠明顯供分析使用。
由臺中地區、宜蘭地區與臺東地區測站資料解算地震事件I之震源P波加速度時間函數結果可見，各地區解算同一個地震震源加速度變化一致，故推論本研究應為研究震源特性之可信方法。臺南地區測站解算地震事件I之震源P波加速度時間函數結果與其他地區測站解算結果比較可見其加速度變化相對較低，由此地區地下構造推論可能為沖積層造成之場址效應造成，然而更精確的原因尚待進一步的研究。
由宜蘭地區測站資料解算12個地震事件震源P波加速度時間函數，地震事件A-I的P波的震源時間函數具有高線性相關性，震央分布沿著宜蘭平原向東延伸至琉球溝與沖繩槽之間，在大地構造上屬於菲律賓海板塊與歐亞板塊碰撞隱沒區，於中央氣象局地震活動分區中屬於花蓮分區。事件J的P波加速度變化以較低的加速度變化率上升到最大值，並且以較低的加速度變化率下降，震央位於海岸山脈至花東海盆之間，在大地構造上屬於菲律賓海板塊與歐亞板塊弧陸碰撞區，於中央氣象局之地震活動分區中屬於臺東分區。事件K和L具有相似的P波加速度變化，具有高線性相關，震央位於臺灣本島中央，於大地構造上屬於板塊運動形成造山運動之陸地抬升區，於中央氣象局之地震活動分區中屬於臺中分區。
研究結果顯示不同分區之震源加速度時間函數具有不同加速度變化，可知各地區之孕震構造有不同的地震加速度特性，推論位於不同孕震區的震源可能具有不同的錯動過程，在應力釋放時以不同的加速度變化反應其差異。

This study uses complex demodulation, attenuation structure, and ray tracing to compute the source time function (STF) of P-wave acceleration based on the earthquake data recorded at the Taichung, Tainan, Yilan, and Taitung stations that are operational under the Taiwan Strong Motion Instrumental Program (TSMIP) of the Central Weather Bureau. Complex demodulation estimates the variation of acceleration with time at each frequency for the seismic records. The acceleration of source P-waves at each frequency are first calculated by determining the propagation path and travel time of seismic waves using ray tracing. Subsequently, the attenuation factor for each block in the three-dimensional attenuation structure traversed by the waves is computed using the propagation path. The acceleration of source P-waves is calculated based on the travel time, attenuation factors, as well as the data at each frequency. The STF of P-wave acceleration is calibrated by substituting the P-wave amplitude for each frequency into the source time function of acceleration, to obtain the STF of P-wave acceleration. The regional seismic results achieved using this method agree with the geotectonic features analyzed by the Central Weather Bureau. The reliability of the protocol developed in this study is therefore corroborated. Apart from velocity based methods, this study offers a different technique of accurately solving source time functions.

一、中文文獻
工程地震學研究室 (2017) 國立中正大學地球與環境科學系，http://www.eq.ccu.edu.tw/lab/lab520/introduction.html（2016年6月查詢）。
中央氣象局 (2016) 地震活動彙整，中央氣象局，http://www.cwb.gov.tw/V7/earthquake/rtd_eq.htm（2016年3月查詢）。
中央氣象局地球物理資料管理系統 (2017) 自由場強地動觀測網，http://gdms.cwb.gov.tw/gdms-freefield.php（2016年3月查詢）。
中央氣象局地震測報中心 (2009) 中央氣象局地震測報中心20週年專刊，中央氣象局地震測報中心，http://scweb.cwb.gov.tw/Twenty.aspx?ItemId=8&loc=tw（2016年3月查詢）。
王仁聖 (2011) 以時頻分析探討地動加速度譜之異常。國立臺灣師範大學碩士論文，共90頁。
王郁如 (2010) 台灣地殼及頂部地函三維P波、S波衰減模型對於造山帶構造特性與機制之探討。國立中央大學博士論文，共159頁。
何春蓀 (1986) 臺灣地質概論。經濟部中央地質調查所。
何英銳 (2002) 地磁複變解調與地震活動之相關性：以台灣崙坪全磁力資料為例。 臺灣師範大學碩士論文，共54頁。
陳光榮 (1993) 臺灣地區Qp值之空間分佈及其特性。國立臺灣大學博士論文，共153頁。
劉貞伶 (2009) 以複變解調法分析竹子湖強震站最大地動加速度的異常放大現象。國立臺灣師範大學碩士論文，共55頁。
鄭世楠，葉永田，徐明同，與辛在勤 (1999) 台灣十大災害地震圖集。中央氣象局，報告編號 CWB-9-1999-002-9，p. 153-174.

二、英文文獻
Ammon, C.J., Velasco, A.A., and Lay, T. (2006). Rapid estimation of first‐order rupture characteristics for large earthquakes using surface waves: 2004 Sumatra‐Andaman earthquake. Geophysical Research Letters, vol. 33, no. 14, L14314. doi:10.1029/2006GL026303.
Banks, R.J. (1975). Complex Demodulation of Geomagnetic Data and the Estimation of Transfer Functions. Geophysical Journal of the Royal Astronomical Society, vol. 43, no. 1, p. 87-101. doi:10.1111/j.1365-246X.1975.tb00629.x.
Chang, C.P., Angelier, J., and Huang, C.Y. (2000). Origin and evolution of a mélange: the active plate boundary and suture zone of the Longitudinal Valley, Taiwan. Tectonophysics, vol. 325, no. 1–2, p. 43-62. doi:http://dx.doi.org/10.1016/S0040-1951(00)00130-X.
Chen, K.J. (1998). S-wave attenuation structure in the Taiwan area and its correlation to seismicity. Terrestrial Atmospheric and Oceanic Sciences, vol. 9, no. 1, p. 97-118.
Chen, K.J., Chiu, B., Wang, J.S., Lee, C.Y., Lin, C.H., and Chao, K. (2011). Geomagnetic anomaly at Lunping before the 1999 Chi-Chi earthquake (Mw= 7.6) in Taiwan. Natural hazards, vol. 58, no. 3, p. 1233-1252.
Chen, K.J., Ho, Y.R., Chiu, B., Wang, J.S., Sun, R., and Lin, C.H. (2007). Correlation between time change in modulus of short-period geomagnetic variation and seismicity in Taiwan. Terrestrial Atmospheric and Oceanic Sciences, vol. 18, no. 3, p. 577-591.
Chen, K.J., Yeh, Y.H., and Shyu, C.T. (1996). Qp structure in the Taiwan area and its correlation to seismicity. Terrestrial Atmospheric and Oceanic Sciences, vol. 7, no. 4, p. 409-429.
Chen, K.J., Yeh, Y.H., and Tsai, Y.B. (1985). A study on complex demodulation of the geomagnetic total intensity data of Taiwan. Bulletin of the Institute of Earth Sciences, Academia Sinica, vol. 5, p. 125-136.
Chen, Y., and Xu, L. (2000). A time‐domain inversion technique for the tempo‐spatial distribution of slip on a finite fault plane with applications to recent large earthquakes in the Tibetan Plateau. Geophysical Journal International, vol. 143, no. 2, p. 407-416.
Courboulex, F., Virieux, J., Deschamps, A., Gibert, D., and Zollo, A. (1996). Source investigation of a small event using empirical Green's functions and simulated annealing. Geophysical Journal International, vol. 125, no. 3, p. 768-780.
Ho, C.S. (1986). A synthesis of the geologic evolution of Taiwan. Tectonophysics, vol. 125, no. 1, p. 1-16. doi:http://dx.doi.org/10.1016/0040-1951(86)90004-1.
Julian, B., and Gubbins, D. (1977). Three-dimensional seismic ray tracing. Journal of Geophysics, vol. 43, no. 1, p. 95-114.
Kikuchi, M., and Kanamori, H. (1991). Inversion of complex body waves—III. Bulletin of the Seismological Society of America, vol. 81, no. 6, p. 2335-2350.
Lee, W.H.K., and Lahr, J.C. (1972). HYPO71: A computer program for determining hypocenter, magnitude, and first motion pattern of local earthquakes (2331-1258). Retrieved from. 69Open-File Report (No. 6972-224), U.S. Geological Survey, 100p.
Ma, K.F., Mori, J., Lee, S.J., and Yu, S. (2001). Spatial and temporal distribution of slip for the 1999 Chi-Chi, Taiwan, earthquake. Bulletin of the Seismological Society of America, vol. 91, no. 5, p. 1069-1087.
Ma, K.F., Wang, J.H., and Zhao, D. (1996). Three-Dimensional Seismic Velocity Structure of the Crust and Uppermost Mantle beneath Taiwan. Journal of Physics of the Earth, vol. 44, no. 2, p. 85-105.
Rau, R.J., and Wu, F.T. (1995). Tomographic imaging of lithospheric structures under Taiwan. Earth and Planetary Science Letters, vol. 133, no. 3, p. 517-532.
Roecker, S., Yeh, Y., and Tsai, Y. (1987). Three‐dimensional P and S wave velocity structures beneath Taiwan: deep structure beneath an arc‐continent collision. Journal of Geophysical Research: Solid Earth, vol. 92, no. B10, p. 10547-10570.
Roumelioti, Z., Kiratzi, A., and Dreger, D. (2004). The source process of the 2001 July 26 Skyros Island (Greece) earthquake. Geophysical Journal International, vol. 156, no. 3, p. 541-548.
Scherbaum, F. (1990). Combined inversion for the three-dimensional Q structure and source parameters using microearthquake spectra. Journal of Geophysical Research: Solid Earth, vol. 95, no. B8, p. 12423-12438. doi:10.1029/JB095iB08p12423.
Zhang, W., Iwata, T., Irikura, K., Pitarka, A., and Sekiguchi, H. (2004). Dynamic rupture process of the 1999 Chi‐Chi, Taiwan, earthquake. Geophysical research letters, vol. 31, no. 10. doi:10.1029/2004GL019827.
Zollo, A., Capuano, P., and Singh, S.K. (1995). Use of small earthquake records to determine the source time functions of larger earthquakes: an alternative method and an application. Bulletin of the Seismological Society of America, vol. 85, no. 4, p. 1249-1256.