臺灣半動態基準之建立與展望

Establishment and Further Development of Taiwan Semi-dynamic Datum

景國恩(Kuo-En Ching)；楊名(Ming Yang)；陳鶴欽(He-Chin Chen)；林文勇(Wen-Yung Lin)；梁旭文(Hsu-Wen Liang)；劉正倫(Jeng-Lun Liu)

半動態基準 ； 地表變形模式 ； 克立金空間內插 ； 斷層錯位模型 ； 電子化全球衛星即時動態定位系統 ； Semi-dynamic datum ； Surface deformation model ； Kriging spatial interpolation ； Fault dislocation model ； e-GNSS

國土測繪與空間資訊

5卷2期（2017 / 07 / 01）

83 - 109

繁體中文

臺灣位於板塊邊界的區域，會因板塊間的相對運動使得靜態國家坐標系統隨時間的推進而漸失去其原有的精度。然而，臺灣現行的TWD97[2010]仍屬於靜態的坐標系統，故如何有效延續此國家坐標系統的精確性乃為國家測繪工作的一項重要任務，而建立臺灣地區的半動態基準則扮演了相當重要的角色。半動態基準包含了靜態大地基準(即TWD97[2010])與地表變形。地表變形模式則包含了「速度網格模型」和「位移網格模型」。在建立速度網格模型方面，本研究完成全臺2000 年01 月至2016 年09 月437 個GPS連續站資料之計算，獲得全臺的水平速度場，並結合785個經濟部中央地質調查所2002 年至2015 年之GPS水平速度場，利用克立金空間內插法建立臺灣地區的水平速度網格模型。在建立位移網格模型方面，本研究藉由比較利用不同精度星曆計算得到的同震位移場，顯示所有計算結果近乎一致，亦即利用超快速星曆求解GPS資料即可在最短時間內建立高精度的近即時地表同震位移場。同時在建置震源模型方面，亦指出單一斷層滑移量模型所獲得的位移網格模型較空間內插模型所得到的結果更加合理。本研究進一步利用半動態基準的地表變形模型檢核內政部國土測繪中心2013 年與2016 年基本控制點檢測成果。採用半動態基準可將2013 年之坐標轉換至2016 年的計算坐標，再將其與2016 年實際觀測的檢測坐標成果進行差異比較發現，兩者有相當高的一致性，有91%的測站坐標差異量在6 公分以內。而剩餘未能被適當估算的基本控制點多分布在山區，可能是由大規模崩塌等區域活動所造成。本研究也建議半動態基準可以配合內政部國土測繪中心現行之e-GNSS系統，利用e-GNSS基準站的坐標時間序列協助檢驗地表變形模式之速度模型精度變化，幫助決定地表速度場更新時機之判斷。最後，就高程半動態基準的發展而言，除了解析長期穩定的震間垂直速度場之外，未來還需要掌握例如地層下陷等非自然因素造成的地表垂直位移，以達到高程基準更新維護之目的。

The accuracy of static national coordinate system in any country, located at the present-day plate boundary such as Taiwan, will be decreased due to the relative motions among different plates. However, the current national coordinate system in Taiwan, TWD97[2010], is still a static datum. Therefore, maintaining the accuracy of national coordinate system in high level efficiently is the main issue in Taiwan. One solution of this issue is to establish a semi-dynamic datum. The semi-dynamic datum is consisted of a static geodetic datum, i.e., TWD97[2010] in Taiwan and a surface deformation model, containing the velocity grid model and the displacement grid model. In this study, observations of the 437 continuous GPS sites from Jan. 2001 to Sep. 2016 were adapted to estimate the horizontal velocity field. We also integrated the horizontal velocity field between 2002 and 2015 from 785 campaign-mode GPS sites surveyed by the Central Geological Survey to derive the horizontal grid velocity model in Taiwan by using the Kriging spatial interpolation method. Furthermore, the coseismic displacements of 2016 Meinong earthquake estimated using ultra-rapid ephemerides are consistent with the results using final ephemerides. In other words, the near real-time coseismic displacement field would be solved using ultra-rapid ephemerides in the future. In addition, the displacement grid model inferred from the uniform slip model is more suitable than the model derived from the spatial interpolation. The resurveying results of the basic control points from the National Land Surveying and Mapping Center (NLSC) in 2013 and 2016, respectively, were used for the exterior checking on the accuracy of semi-dynamic datum in this study. The calculated coordinates in 2016, derived from observations in 2013 using the semi-dynamic datum, are nearly corresponding with the observed coordinates in 2016. The coordinate differences in 91% stations are less than 6 cm. The outliers of comparisons are mostly located in the mountain range, which may be caused by the deep-seated landslides. In addition, the e-GNSS established by the NLSC is proposed to be used for checking the accuracy of velocity model and for the timing decision of velocity model renewal. Finally, in order to develop the vertical semi-dynamic datum, we have to estimate the vertical velocity field and collect land subsidence information simultaneously for improving and maintaining the vertical datum in a high accuracy level.

1. Bos, A. G.,Spakman, W.,Nyst, M. C. J.(2013).Surface deformation and tectonic setting of Taiwan inferred from a GPS velocity field.J. Geophys. Res.,108,2458.
2. Chang, C.-P.,Chang, T.-Y.,Angelier, J.,Kao, H.,Lee, J.-C.,Yu, S.-B.(2003).Strain and stress field in Taiwan oblique convergent system: Constraints from GPS observations and tectonic data.Earth Planet. Sci. Lett.,214,115-127.
3. Ching, K.-E.,Chen, K.-H.(2015).Tectonic effect for establishing a semi-dynamic datum in Southwest Taiwan.Earth Planet. Sci. Lett.,67(1),1-14.
4. Ching, K.-E.,Johnson, K. M.,Rau, R.-J.,Chuang, R. Y.,Kuo, L.-C.,Leu, P.-L.(2011).Inferred fault geometry and slip distribution of the 2010 Jiashian, Taiwan, earthquake is consistent with a thick-skinned deformation model, Earth Planet.Earth Planet. Sci. Lett.,301,78-86.
5. Chiu, Y.-H.,Shih, P. T.-Y.(2014).National Datum Uncertainty due to Reference Frame Transformation: Case Study for the Geodetic Datum of Taiwan.Journal of Surveying Engineering,140(3),05014002.
6. Grant, D. B.,Pearse, M. B.(1995).Proposal for a dynamic national geodetic datum for New Zealand.IUGG XXI General Assembly,Boulder Colorado:
7. Grant, D.B.,Blick, G. H.,Pearse, M. B.,Beavan, R. J.,Morgan, P. J.(1999).The development and implementation of New Zealand Geodetic Datum 2000.IUGG99 General Assembly,Birmingham UK:
8. Han, J. Y.(2010).Noniterative approach for solving the indirect problems of linear reference frame transformations.Journal of Surveying Engineering,136(4),150-156.
9. Han, J. Y.,Yu, S. W.,van Gelder, B. H. W.(2011).Time-variant reference frame transformations in a deforming area.Survey Review,43(321),284-295.
10. Hsu, Y.-J.,Yu, S.-B.,Simons, M.,Kuo, L.-C.,Chen, H.-Y.(2009).Interseismic crustal deformation in the Taiwan plate boundary zone revealed by GPS observations, seismicity, and earthquake focal mechanisms.Tectonophysics,479,4-18.
11. Lin, A.-T.,Watts, A. B.,Hesselbo, S. P.(2003).Cenozoic stratigraphy and subsidence history of the South China Seamargin in theTaiwan region.Basin Res.,15,453-478.
12. Nikolaidis, R.(2002).San Diego,University of California.
13. Okada, Y.(1985).Surface deformation due to shear and tensile faults in a half-space.Bull. Seism. Soc. Am.,75,1135-1154.
14. Okada, Y.(1992).Internal deformation due to shear and tensile faults in a half-space.Bull. Seism. Soc. Am.,82,1018-1040.
15. Sagiya, T.,Miyazaki, S.,Tada, T.(2000).Continuous GPS array and present-day crustal deformation of Japan.PAGEOPH,157,2303-2322.
16. Tregoning, P.,Jackson, R.(1999).The need for dynamic datums.Geomatics Research Australasia,71,87-102.
17. Ward, S. N.(1994).A multidisciplinary approach to seismic hazard in southern California.Bull. Seism. Soc. Am.,84,1293-1309.
18. Yang, M.,Chen, K. H.,Shiao, S. W.(2003).A new height reference network in Taiwan.Surv. Rev.,37,260-268.
19. Yang, M.,Tseng, C. L.,Yu, J. Y.(2001).Establishment and maintenance of Taiwan geodetic datum 1997.J. Surv. Eng.,127,119-132.
20. Yoshiyuki, T.,Hiroaki, S.,Jun, S.,Kazumi, I.,Tomoo, T.,Hideaki, H.,Tamotsu, K.,Shigeru, M.,Yuki, H.,Mikio, T.,Yuki, K.,Tetsuro, I.(2007).Efficient maintenance of the Japanese geodetic datum 2000 using crustal deformation models - PatchJGD & semi-dynamic datum.Bulletin of the Geographical Survey Institute,54.
21. Yu, S.-B.,Chen, H.-Y.(1994).Global positioning system measurements of crustal deformation in the Taiwan arc-continent collision zone.Terr. Atmos. Ocean. Sci.,5,477-498.
22. Yu, S.-B.,Chen, H.-Y.,Kuo, L.-C.(1997).Velocity field of GPS Stations in the Taiwan area.Tectonophysics,274,41-59.
23. Yu, S.-B.,Kuo, L.-C.(2001).Present-day crustal motion along the Longitudinal Valley Fault, eastern Taiwan.Tectonophysics,333,199-217.
24. Yu, S.-B.,Kuo, L.-C.,Hsu, Y.-J.,Su, H.-H.,Liu, C.-C.,Hou, C.-S.,Lee, J.-F.,Lai, T.-C.,Liu, C.-C.,Liu, C.-L.,Tseng, T.-F.,Tsai, C.-S.,Shin, T.-C.(2001).Preseismic deformation and coseismic displacements associated with the 1999 Chi-Chi, Taiwan earthquake.Bull. Seism. Soc. Am.,91,995-1012.
25. Yu, S.-B.,Kuo, L.-C.,Punongbayan, R.-S.,Ramos, E.-G.(1999).GPS observation of crustal motion in the Taiwan-Luzon region.Geophys. Res. Lett.,26,923-926.
26. 內政部國土測繪中心(2013)。,國立成功大學。
27. 內政部國土測繪中心(2015)。,國立臺北大學。
28. 內政部國土測繪中心(2014)。,國立臺北大學。
29. 內政部國土測繪中心(2016)。,國立臺北大學。
30. 內政部國土測繪中心(2011)。,國立交通大學。
31. 內政部國土測繪中心(2012)。,內政部國土測繪中心。
32. 邱元宏(2016)。博士論文(博士論文)。國立交通大學土木工程學系。

1. 陳國華,王昱凡(2020)。臺灣水平地表變形模式精度分析與應用。國土測繪與空間資訊，8(2)，61-77。
2. 賴宗裕,徐名洪(2023)。台灣土地變動特性與地籍相關問題之探討。測量工程，60，58-81。
3. 蕭詩涵,景國恩,陳鶴欽,陳國華,林文勇,周逸屏(2019)。臺灣半動態基準地表變形模式精進研究－同震格點位移模型新增時機之探討。國土測繪與空間資訊，7(1)，33-66。