許多地震發生之後，都可以偵測到重力的改變，這是因為岩層受到應力作用而變形或破裂，使得地表地形及地下構造改變，這些現象都會造成地表重力值的永久性變化。透過精密的重力測量資料，我們可以了解地震對台灣的影響。本研究大致可分為兩個部分。第一部分，利用新竹市十八尖山超導重力儀 (Superconducting Gravimeter，簡稱SG)，所記錄到的實際重力觀測資料 (2006/4-2011/9)，去除各種環境引起的訊號，如：日、月引力產生的固體潮、海潮以及地下水位、氣壓、極移等因素，分析地震前後，因板塊或斷層位移所造成的重力變化。第二部分，使用台灣寬頻地震網 (Broadband Array in Taiwan for Seismology，簡稱BATS)，2006年4月至2011年9月，台灣本島附近，深度小於300公里，地震矩規模 (Mw)大於5的地震資料，配合Fu and Sun (2008)的彈性地球模型，計算各地震所產生的永久重力變化量。此模型可計算地下構造和地表受到地震影響而產生的重力變化。本研究最後比較儀器記錄數據與理論模型計算結果，發現兩者之間結果正負號相同但大小不同，主要原因可能是台灣的地質構造過於複雜，導致模型無法準確估計地震對台灣的影響。

After earthquakes, we can detect the Gravity changes. This is because the rock will deform or brake by tectonic stresses, both them will cause permanent gravity change. Therefore, through the data of precision gravimeter, we can understand the influence of the earthquake on Taiwan. This research is divided into two parts. In the first part, we took gravity data (2006/4-2011/9) from Hsinchu superconducting gravimeter (SG) station, and removed those data we didn’t want such as solid earth tide, ocean tide, polar motion, atmospheric, continental water storage then we had got residual gravity data. Then we used the residual gravity data to compute permanent gravity changes caused by quakes. In the second part, we used earthquakes data which is start from 2006/4-2011/9, near Taiwan, depth around 0 to 300 km, moment magnitude (Mw) greater than 5 from the Broadband Array in Taiwan for Seismology (BATS). Than used Fu and Sun (2008) numerical earth model with earthquakes data from BATS, to compute permanent gravity changes caused by earthquakes.
Finally, we analyzed and compared SG data result with the model computing result. The result show that the sign of two results is the same but the value’s order is quite different. The main reason why maybe is the complicated tectonic in Taiwan. It cause the model could not compute the influence of the earthquake precisely. Complicated tectonic also influence the quality of SG data, although we could get clear signals by many data processes, we could not sure that it is meaning the gravity changed by the earthquake.

〔1〕Baker, T. F., and M. S. Bos ,“Validating Earth and ocean tide models using tidal gravity measurements”, Geophysical Journal International, 152, 468–485, 2003.
〔2〕Boy, J.P., and J. Hinderer,“Study of the seasonal gravity signal in superconducting gravimeter data”, Journal of Geodynamics, 41(1-3), 227-233,2006.
〔3〕Chao, B. F., and R. s. Gross ,“Changes in the Earth's rotation and low-degree gravitational field induced by earthquakes”, Journal of Geophysical Research, 91, 569-596,1987.
〔4〕Chao, B. F., and R. Eanes,“Global gravitational changes due to atmospheric mass redistribution”, Geophysical Journal International, 122, 755-764,1995.
〔5〕Chao, B. F., R. S. GrossD., and N. Dong, “Changes in the Earth's rotational energy induced by earthquakes”, Geophysical Journal International, 122, 784-789,1995.
〔6〕Crossley, D., et al., “Network of Superconducting Gravimeters Benefits a Number of Disciplines”, EOS Trans, 80(11), 121-126, 1999.
〔7〕D. Dong, T. A. Herring, and R. W. King,“Estimating regional deformation from a combination of space and terrestrial geodetic data”, Journal of Geodsy, 72, 200-214,1998.
〔8〕Dragert, H., T. S. James, and A. Lambert,“Ocean loading corrections for continuous GPS: A case study at the Canadian coastal site Holberg”, Geophysical Research Letters, 27(14), 2045-2048,2000.
〔9〕Fu, G., and W. Sun,“Surface coseismic gravity changes caused by dislocations in a 3-D heterogeneous earth”, Geophysical Journal International, 172(2), 479-503,2008.
〔10〕Goodkind, J. M.,“Continuous measurement of nontidal variations of gravity”, Journal of Geophysical Research, 91(B9), 9125,1986.
〔11〕Goodkind, J. M., and R. J. Warburton,“Superconductivity applied to gravimeter”, IEEE Transactions on Magnetics, 708-711,1975.
〔12〕Heki, K., and K. Matsuo,“Coseismic gravity changes of the 2010 earthquake in central Chile from satellite gravimetry”, Geophysical Research Letters, 37(24), 2010.
〔13〕Hwang, C., R. Kao, C.-C. Cheng, J.-F. Huang, C.-W. Lee, and T. Sato ,“Results from parallel observations of superconducting and absolute gravimeters and GPS at the Hsinchu station of Global Geodynamics Project, Taiwan”, Journal of Geophysical Research, 114(B7) ,2009.
〔14〕Imanishi, Y., T. Sato, T. Higashi, W. Sun, and S. Okubo ,“A network of superconducting gravimeters detects submicrogal coseismic gravity changes”, Science, 306(5695), 476-478,2004.
〔15〕Karato, S. i.,“Importance of anelasticity in the interpretation of seismic tomography”, Geophysical Research Letters, 1623-1626,1993.
〔16〕Lo, C. L., and S. K. Hsu, “Earthquake-induced gravitational potential energy change in the active Taiwan orogenic belt”, Geophysical Journal International, 162(1), 169-176, 2005.
Lundgren, P.,“Source model for the 2001 flank eruption of Mt. Etna volcano”, Geophysical Research Letters, 30(7), 2003.
〔17〕Matsumoto, K., T. Takanezawa, and M. Ooe,“Ocean tide models developed by assimilating TOPEX/POSIDON altimeter data into hydrodynamical model A global model and a regional model around Japan”, Journal of Oceanography, 56, 567–581,2000.
〔18〕Meichior, P., and B. Ducarme,“Detection of inertial gravity oscillations in the Earth's core with a superconducting gravimeter at Brussels’’, 42, 149-134, 1986.
〔19〕Okada, Y.,“Surface deformation to shear and tensile faults in a half-space, Bulletin of the Seismological Society of America’’, 75, 1135-1154,1985.
〔20〕Okubo, S.,“Potential and gravity changes raised by point dislocations’’, Geophys. J. Int., 105, 573-586,1991.
〔21〕Okubo, S.,“ Potential and gravity changes due to shear and tensilefaults in a half-space’’, J. geophys. Res., 97, 7137-7144,1992.

〔22〕Prothero, W. A., and J. M. Goodkind ,“Earth-tide measurements with the superconducting gravimeter’’, Journal of Geophysical Research, 77(5), 926-937,1972.
〔23〕Sun, W., and S. Okubo, “Surface potential and gravity changes due to internal dislocations in a spherical earth-I. Theory for a point dislocation’’, Geophysical Journal International, 114, 569-592, 1993.
〔24〕Sun, W., and S. Okubo, “Surface potential and gravity changes due to internal dislocations in a spherical earth-II. Application to a finite fault’’, Geophysical Journal International, 132, 79-88,1998.
〔25〕Sun, W., and S. Okubo, “Effects of earth's spherical curvature and radial heterogeneity in dislocation studies—for a point dislocation’’, Geophysical Research Letters, 29(12), 2002.
〔26〕Wang, K., “A revised dislocation model of interseismic deformation of the Cascadia subduction zone’’, Journal of Geophysical Research, 108(B1), 2003.
〔27〕Yu, S.B., H.Y. Chen, and L.C. Kuo , “Velocity field of GPS stations in the Taiwan area,’’ Tectonophysics, 274, 41-59, 1997.
〔28〕Yen, Y.T., and K.F. Ma , “Source-Scaling Relationship for M 4.6–8.9 Earthquakes, Specifically for Earthquakes in the Collision Zone of Taiwan’’, Bulletin of the Seismological Society of America, Vol. 101, No. 2, pp. 464–481,2011.
〔29〕孫文科，地震位錯理論，中國，2011。
〔30〕蘇洵頡，利用超導重力儀訊號偵測地震造成的永久性重力變化，國立交通大學，碩士論文，民國99年。