全球定位系統（GPS）具有高精度，全天候觀測的特性，已廣泛應用於地殼變動之範疇。本文以1999年8月到2002年12月台灣地區衛星追蹤站的資料，與IGS國際衛星追蹤站進行聯測，藉此獲得國內衛星追蹤站在ITRF框架下的坐標值與速度量，以明瞭台灣地區板塊運動的情形。

為了節省在運算上所花費的時間，使資料在處理過程更有效率，本文以GPS解算軟體Bernese v4.2自動化計算程序，進行資料處理與平差計算。在資料計算方面，每日計算皆以關島（GUAM）作為約制進行最小二乘平差計算，再將每日之成果進行合併求解。在實驗過程中，發現僅是精密星曆與衛星追蹤站框架不一致所造成國內追蹤站坐標的偏差最大可達約±4cm，此數量級對於地殼運動而言不可忽視，因此本文亦考慮精密星曆與衛星追蹤站參考框架之一致性。

最後分析結果顯示，在921地震之後陽明山（YMSM）、北港（PKGM）、馬祖（MZUM）三個衛星追蹤站之水平速度場量值介於30～40 mm/yr，鳳林（FLNM）與墾丁（KDNM）兩個衛星追蹤站之水平速度場約為20 mm/yr，TMAM衛星追蹤站之水平速度場約為4 mm/yr。此外，各站在垂直方向均呈現下降的情況除北港、鳳林與馬祖兩個衛星追蹤站在垂直方向的速度場偏大（約20mm/yr以上）以外，其餘各衛星追蹤站皆小於10 mm/yr。

The Global Position System (GPS) has been widely applied to monitor the crustal movements for a long time due to its high precision and independence of meteorological conditions. In the thesis, the data of the GPS tracking stations of Taiwan and the international tracking stations of IGPS over the interval August 1999-December 2002 are combined for the computation and analysis of crustal movement and for the establishment of velocity field of GPS tracking stations of Taiwan.

To improve the efficiency of computation, an automatic process based on BPE of the GPS software Bernese v4.2 has been developed. In the daily data processing, the GPS tracking station GUAM of IGS is constrained with its displacement velocity in the network adjustment using Bernese. Our experiment has shown that the inconsistence of reference frame between ephemeris and coordinates of constrained tracking stations can cause a difference up to about 4cm in coordinate. Therefore, that inconsistence cannot be neglected, and is taken into consideration carefully in the thesis.

The results of our experiment have shown that the horizontal velocity of YMSM, PKGM and MZUM tracking stations after 921 earthquake is about 30 ~ 40 mm/yr, FLNM and KDNM about 20 mm/yr, and TMAM about 4 mm/yr . Moreover, all tracking stations of Taiwan descend gradually. YMSM, KDNM, MZUM and TMAM have the vertical displacement velocity less than 10 mm/yr , while PKGM and FLNM have the larger vertical displacement velocity over 20 mm/yr.

中研院地球科學研究所：http://www.earth.sinica.edu.tw/

中央氣象局地震測報中心：http://www.cwb.gov.tw/index-f.htm

內政部，（2000）：九二一集集大地震災害後重建一、二等衛星控制點測量工作報告。

內政部（2001）：衛星追蹤站坐標維護及計算分析。工作成果報告書，國立成功大學衛星資訊研究中心，台南。

尤瑞哲（2003)：測量坐標系統，國立成功大學測量工程學系，台南。

李彥弘、楊翼聰（1996）：內政部GPS衛星測量業務規劃。第二屆GPS衛星定位技術研討會，國立成功大學，台南，第89-102頁。

馬國鳳，活躍的地球，地球科學研究所，地球科學文集網址：http://www.earth.sinica.edu.tw/paper/ma.pdf

郭隆晨（2001）：高精度GPS衛星測量在地殼變形觀測之研究。國立交通大學土木工程學系博士論文，新竹。

郭隆晨、余水倍、陳春盛（1994）：GPS連續觀測資料處理方法研究。第十三屆測量學術研討及應用會論文級，國立交通大學，新竹，第 121-131頁。

張嘉強、陳鴻智（2000）：參考框架轉換對衛星追蹤站坐標之影響，測量工程，42（1）,第 53-68頁。

國際地球參考框架：ftp://lareg.ensg.ign.fr/pub/itrf/ITRF.TP
http://large.ensg.ign.fr/ITRF/ITRF97/itrf97-in/itrf97-in.html#solusion
http://lareg.ensg.ign.fr/ITRF/solutions.html
http://schubert.ign.fr/ITRF/ITRF97/results/ITRF97_GPS.SSC
http://schubert.ign.fr/ITRF/ITRF2000/results/ITRF2000_GPS.SSC

國立中央大學應用地質研究所，工程地質與新科技研究室，地理資訊網址：http://140.115.123.30/gis/eq/twactft/loc.htm

鄭煒騰（2000）：台灣地區GPS衛星追蹤網形形變分析，以1999年觀測資料為例。國立成功大學測量研究所碩士論文，台南。

蔡豐吉（2002）：利用短基線場檢測GPS天線相位中心之研究。國立交通大學土木工程學系碩士論文，新竹。

Beutler, G., H., Bock, E. Brockmann, R. Dach, P. Fridez, W. Gurtner, U. Hugentobler, D. Ineichen, J. Johnson, M. Meindl, L. Mervart, M. Rothacher, S. Schaer, T. Springer, and R. Weber（2001）, Bernese GPS Software Version 4.2, Astronomical Institute, University of Berne, Berne, Switzerland.

Blewitt, G., M. B. Heflin, F. H. Webb, U. J. Lindqwister and R. P. Malla, （1992）：Global Coordinates with Centimeter Accuracy in the International Terrestrial Reference Frame Using GPS. Geophysical Reserch Letters, Vol. 19, No.9, pp. 853-856.

Boucher, C., and Z. Altamimi (1996): International Terrestrial Reference Frame, GPS World, September 1996, pp 71-74.

Boucher, C., Z. Altamimi, M. Feissel, and P. Sillard, (1998):Results and Analysis of the ITRF96, IERS Technical Note 24, Observatoire de Paris, Paris, France.

Boucher, C., Z. Altamimi, M. Feissel, and P. Sillard, (1999): The 1997 international terrestrial reference frame (ITRF97), IERS Technical Note 27, Observatoire de Paris, Paris, France.

Brockmann, E., (1997): Combination of Solutions for Geodynamic Applications of the Global Positioning System(GPS). Funfundfunfzigster Band Vol. 55.

Deffontaines, B., Lee, J.C, Angelier, J Carvalho, j. and Rudant, J.P.,(1994): New Geomorphic Data on the Active Taiwan Orogen:A Multi source Approach. J. Geophys. R. Vol. 99, No B10, pp.20243-20266.

Hugentobler, U., Schaer, S., Fridez, P. ,(2001): Bernese GPS Sofware Version 4.2, Astronmical institute university of Berne.

IGS(1999)，Information and Resources, December 1999.

Kotzev, V., Nakov, R., Burchfiel, B.C., Reilinger, R., (2001): GPS study of active tectonics in Bulgaria:result from 1996 to 1998 . Journal of Geodynamics , Vol. 31, pp. 189-200.

Krumm, F., Grafarend E. and Friedrich, D., (2000): Datum free Deformation Analysis ITRF networks. Geoid and Geodynamics 2000,July 31- August 4 , 2000.

Liu, C.C., and Yu, S.B, (1990): Vertical Crustal Movements in Eastern Taiwan and thei Tectonic Implication. Tectonphysics, 183 , pp. 111-119.

Louis, D. B., (1986): World Geodetic System 1984. Proceedings of the Fouth International Geodetic Symosium on Satellite Positioning, Vol. I.1, pp. 69-88.

SOPAC網站：http://sopac.ucsd.edu/cgi-bin/sector.cgi

Willemin, J.H, and Knuepfer, P.L.K., (1994): Kinematics of Arc-Continent Collision in the eastern central Range of Taiwan Inferred from Geomorphic Analysis. J. Geophy. R. , Vol.99, No.B10, pp. 20267-20280.

Yu, S.B., and Lee, C., (1986): Geodetic Measurement of Horizontal Crustal Deformation in Eastern Taiwan. Tectonophysics, Vol. 125, pp. 73-85.

Yu, S.B., Chen, H.Y., and Kuo, L.C., (1997): Velocity field of GPS station in the Taiwan area. Tectonophysics, Vol. 274 , pp. 41-59.