臺灣位於菲律賓海板塊及歐亞大陸板塊的聚合邊界，複雜的板塊運動使得臺灣有相當高的變形速率及常有斷層活動造成的地震。此外，西南部地區由於農業灌溉及養殖魚業所需大量用水，超抽地下水造成地層下陷，進而危害到建築物、公路及管線的安全。因此，若能有效地大範圍地偵測地表變形，可有助於減少災害的發生。合成孔徑雷達差分干涉技術(Differential Interferometric Synthetic Aperture Radar, DInSAR) 常被利用來偵測地表變形，在過去許多研究中皆證明它能觀測地表大範圍且精準的三維變動。由歐洲太空總署(European Space Agency, ESA)在2014年發射的Sentinel-1A提供C波段合成孔徑雷達(SAR)影像，並使用Interferometric Wide Swath (IW)模式作為陸地上的觀測模式，IW 模式擁有長170公里及寬250公里的大範圍覆蓋面積，而且Sentinel-1A的短再訪週期每12天可提供一幅臺灣地區的影像。因此利用Sentinel-1A的大覆蓋面積及再訪週期間隔短的特性，可快速且有效地取得整個臺灣的變形圖。
然而，臺灣位於電離層赤道異常地區(Equatorial Ionization Anomaly, EIA)，不同時間的全電子含量(Total Electron Content, TEC)會隨著EIA強弱而改變，兩幅影像中若帶有TEC差異將會在干涉圖(interferogram)上產生額外的相位條帶，進而影響相位解纏及變形圖的結果。此研究比較了全球電離層模式(Global Ionosphere Map, GIM)及臺灣區域電離層模式(Taiwan Regional Ionospheric Map, TRIM)之TEC，模擬電離層在干涉圖中產生的額外相位，並在原始的干涉圖中進行修正，得到電離層修正後的干涉圖。利用2016及2017年Sentinel-1A升軌影像所產生的47對干涉圖中，有2對明顯受到電離層效應並可以利用TRIM的TEC進行消除。經本研究統計，若兩幅影像TEC值在350公里穿透層於南北向的梯度差超過16 TECU (TECU=1016電子數/m2)，則會在干涉圖中產生明顯的電離層效應。最後，利用Global Positioning System (GPS)地面站的時間序列變化比較電離層修正前及修正後的結果，結果顯示經電離層修正後的變形量與地面觀測較為一致，顯示電離層確實會造成變形量的誤差。經修正後，地表垂直向變動的均方根誤差(Root Mean Square Error, RMSE)由修正前的42公厘減少為11公厘。

Taiwan is located at the convergent boundary between the Eurasian Plate and the Philippine Sea Plate. This tectonically active region produces high surface deformation rate and fault activities. In addition, over-pumping of groundwater for agriculture and aquaculture needs in the southwestern part of Taiwan has induced extra land subsidence in the past 20 years. A systematic monitoring method for land deformation using advanced remote sensing technology is thus needed. The Differential Interferometric Synthetic Aperture Radar (DInSAR) technique has been proven useful to measure surface deformation in an effectively way. Sentinel-1A satellite launched by the European Space Agency (ESA) is one of them acquiring images with various pre-defined scenarios to serve different thematic domains. It operates mostly in the Interferometric Wide Swath (IW) mode with dual polarization in east Asia, which provides a broad coverage with 250 km in range and 170 km in azimuth direction. It orbit characteristics allow users to produce interferograms for most of terrestrial areas every 12 days. However, Taiwan is situated on the north edge of the Equatorial Ionization Anomaly (EIA) Region. The Total Electron Content (TEC) irregularity between snapshots causes ionospheric effect on the interferograms. The ionospheric residuals induce extra fringes in the interferograms that further affect unwrapping and displacement result.
In this study, we utilize global and regional ionospheric vertical TEC maps, namely the Global Ionosphere Map (GIM) released by the International GNSS Service and the Taiwan Regional Ionospheric Map (TRIM) provided by the Ionospheric Radio Science Laboratory of National Central University, to compensate ionospheric phase in the interferograms. Our results show that TRIM is capable to reduce most of extra fringes caused by the ionospheric effect. Comparing these two types of TEC map, TRIM presents more accurate and detailed ionospheric patterns than GIM. In a series of 47 interferograms generated from Sentinel-1A ascending mode in 2016–2017, we conclude that when the TEC gradient in south-north direction exceeds 16 TECU (1016 electrons/m2), it will cause extra fringes in the interferogram. Finally, the Global Positioning System (GPS) continuous station measurements are used to compare with DInSAR displacement for the performance of ionospheric correction. The results show that after correction the displacement is more consistent with ground truth. It also confirms that the ionosphere may cause serious error in the displacement estimates. After ionospheric correction, the root-mean-square error (RMSE) of estimated vertical motion is reduced from 42 mm to 11 mm.

[1] Gomba, G., González, F. R., & De Zan, F. (2017). Ionospheric phase screen compensation for the Sentinel-1 TOPS and ALOS-2 ScanSAR modes. IEEE Transactions on Geoscience and Remote Sensing, 55(1), 223-235.
[2] Yu, S. B., Chen, H. Y., & Kuo, L. C. (1997). Velocity field of GPS stations in the Taiwan area. Tectonophysics, 274(1-3), 41-59.
[3] Ma, K. F., Lee, C. T., Tsai, Y. B., Shin, T. C., & Mori, J. (1999). The Chi‐Chi, Taiwan earthquake: Large surface displacements on an inland thrust fault. Eos, Transactions American Geophysical Union, 80(50), 605-611.
[4] Wen, S., Yeh, Y. L., Chang, Y. Z., & Chen, C. H. (2017). The seismogenic process of the 2016 Meinong earthquake, southwest Taiwan. Terrestrial, Atmospheric & Oceanic Sciences, 28(5).
[5] Hung, W. C., Hwang, C., Chen, Y. A., Chang, C. P., Yen, J. Y., Hooper, A., & Yang, C. Y. (2011). Surface deformation from persistent scatterers SAR interferometry and fusion with leveling data: A case study over the Choushui River Alluvial Fan, Taiwan. Remote Sensing of Environment, 115(4), 957-967.
[6] Information for Prevention and Monitoring of Land Subsidence in Taiwan: http://www.lsprc.ncku.edu.tw/index.php
[7] Chang, C. P., Angelier, J., & Huang, C. Y. (2000). Origin and evolution of a melange: the active plate boundary and suture zone of the Longitudinal Valley, Taiwan. Tectonophysics, 325(1-2), 43-62.
[8] Muja, R., Vasile, G., Boudon, R., D'Urso, G., & Boldo, D. (2012). Multidimensional Very High Resolution SAR Signal Processing for Monitoring Energetic Structures (Doctoral dissertation, GIPSA-LAB).
[9] Zebker, H., Werner, C., Rosen, P., & Hensley, S. (1993). Accuracy of topographic maps derived from ERS-1 interferometric radar.
[10] Feng, G., Li, Z., Shan, X., Zhang, L., Zhang, G., & Zhu, J. (2015). Geodetic model of the 2015 April 25 M w 7.8 Gorkha Nepal Earthquake and M w 7.3 aftershock estimated from InSAR and GPS data. Geophysical journal international, 203(2), 896-900.
[11] Wang, X., Liu, G., Yu, B., Dai, K., Zhang, R., Chen, Q., & Li, Z. (2014). 3D coseismic deformations and source parameters of the 2010 Yushu earthquake (China) inferred from DInSAR and multiple-aperture InSAR measurements. Remote sensing of environment, 152, 174-189.
[12] Herrera, G., Tomás, R., Monells, D., Centolanza, G., Mallorquí, J. J., Vicente, F., ... & Mulas, J. (2010). Analysis of subsidence using TerraSAR-X data: Murcia case study. Engineering Geology, 116(3-4), 284-295.
[13] Tomás, R., Romero, R., Mulas, J., Marturià, J. J., Mallorquí, J. J., Lopez-Sanchez, J. M., ... & Duque, S. (2014). Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain. Environmental earth sciences, 71(1), 163-181.
[14] Chen, Q., Cheng, H., Yang, Y., Liu, G., & Liu, L. (2014). Quantification of mass wasting volume associated with the giant landslide Daguangbao induced by the 2008 Wenchuan earthquake from persistent scatterer InSAR. Remote sensing of environment, 152, 125-135.
[15] Singleton, A., Li, Z., Hoey, T., & Muller, J. P. (2014). Evaluating sub-pixel offset techniques as an alternative to D-InSAR for monitoring episodic landslide movements in vegetated terrain. Remote Sensing of Environment, 147, 133-144.
[16] Hooper, A., Zebker, H., Segall, P., & Kampes, B. (2004). A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers. Geophysical research letters, 31(23).
[17] RADAR and SAR Glossary: https://earth.esa.int/handbooks/asar/CNTR5-2.html
[18] InSAR Technology: http://www.insar.sk/en/insar_technology/
[19] Ferretti, A., Monti Guarnieri, A., Prati, C., Rocca, F., & Massonnet, D. (2007). INSAR Principles A. ESA publications.
[20] Kelly, M. C. (1989). The Earth’s ionosphere. Int. Geophys. Ser, 43, 71.
[21] Liu, L., Wan, W., Ning, B., Pirog, O. M., & Kurkin, V. I. (2006). Solar activity variations of the ionospheric peak electron density. Journal of Geophysical Research: Space Physics, 111(A8).
[22]Solar Minimum; Solar Maximum: https://www.nasa.gov/mission_pages/sunearth/news/solarmin-max.html
[23] Solar cycle: https://en.wikipedia.org/wiki/Solar_cycle
[24] SOLAR CYCLE PROGRESSION: https://www.swpc.noaa.gov/products/solar-cycle-progression
[25] Kumar, S., Singh, A. K., & Lee, J. (2014). Equatorial Ionospheric Anomaly (EIA) and comparison with IRI model during descending phase of solar activity (2005–2009). Advances in Space Research, 53(5), 724-733.
[26] Appleton, E. V. (1946). Two anomalies in the ionosphere. Nature, 157(3995), 691.
[27] Rishbeth, H., Müller-Wodarg, I. C. F., Zou, L., Fuller-Rowell, T. J., Millward, G. H., Moffett, R. J., ... & Aylward, A. D. (2000, August). Annual and semiannual variations in the ionospheric F2-layer: II. Physical discussion. In Annales Geophysicae (Vol. 18, No. 8, pp. 945-956). Springer-Verlag.
[28] Susi, M. (2017). Improved receiver tracking models for scintillation monitoring (Doctoral dissertation, University of Nottingham).
[29] Nikolaidis, R. (2002). Observation of geodetic and seismic deformation with the Global Positioning System.
[30] Fattahi, H., Simons, M., & Agram, P. (2017). InSAR Time-Series Estimation of the Ionospheric Phase Delay: An Extension of the Split Range-Spectrum Technique. IEEE Transactions on Geoscience and Remote Sensing, 55(10), 5984-5996.
[31] Belcher, D. P. (2008). Theoretical limits on SAR imposed by the ionosphere. IET Radar, Sonar & Navigation, 2(6), 435-448.
[32] Gomba, G., Parizzi, A., De Zan, F., Eineder, M., & Bamler, R. (2016). Toward operational compensation of ionospheric effects in SAR interferograms: The split-spectrum method. IEEE Transactions on Geoscience and Remote Sensing, 54(3), 1446-1461.
[33] Jung, H. S., Lee, D. T., Lu, Z., & Won, J. S. (2013). Ionospheric correction of SAR interferograms by multiple-aperture interferometry. IEEE Transactions on Geoscience and Remote Sensing, 51(5), 3191-3199.
[34] Meyer, F., Bamler, R., Jakowski, N., & Fritz, T. (2006). The potential of low-frequency SAR systems for mapping ionospheric TEC distributions. IEEE Geoscience and Remote Sensing Letters, 3(4), 560-564.
[35] Brcic, R., Parizzi, A., Eineder, M., Bamler, R., & Meyer, F. (2011, July). Ionospheric effects in SAR interferometry: An analysis and comparison of methods for their estimation. In Geoscience and Remote Sensing Symposium (IGARSS), 2011 IEEE International(pp. 1497-1500). IEEE.
[36] European Space Agency. (2013). Sentinel-1 User Handbook.
[37] Sowter, A., Amat, M. B. C., Cigna, F., Marsh, S., Athab, A., & Alshammari, L. (2016). Mexico City land subsidence in 2014–2015 with Sentinel-1 IW TOPS: Results using the Intermittent SBAS (ISBAS) technique. International Journal of Applied Earth Observation and Geoinformation, 52, 230-242.
[38] INSAR TECHNICAL BACKGROUND: https://skygeo.com/insar-technical-background/ 38
[39] NASA’s Archive of Space Geodesy Data ftp://cddis.gsfc.nasa.gov/gnss/products/ionex/
[40] Jee, G., Lee, H. B., Kim, Y. H., Chung, J. K., & Cho, J. (2010). Assessment of GPS global ionosphere maps (GIM) by comparison between CODE GIM and TOPEX/Jason TEC data: Ionospheric perspective. Journal of Geophysical Research: Space Physics, 115(A10).
[41] IGS network: http://www.igs.org/network?network=multi-GNSS
[42] Space Weather Operational Office: http://swoo.cwb.gov.tw/V1/index.htm
[43] GMTSAR: http://topex.ucsd.edu/gmtsar/
[44] Chen, C. W., & Zebker, H. A. (2000). Network approaches to two-dimensional phase unwrapping: intractability and two new algorithms. JOSA A, 17(3), 401-414.
[45] Chen, C. W., & Zebker, H. A. (2001). Two-dimensional phase unwrapping with use of statistical models for cost functions in nonlinear optimization. JOSA A, 18(2), 338-351.
[46] Chen, C. W., & Zebker, H. A. (2002). Phase unwrapping for large SAR interferograms: Statistical segmentation and generalized network models. IEEE Transactions on Geoscience and Remote Sensing, 40(8), 1709-1719.
[47] SNAPHU: Statistical-Cost, Network-Flow Algorithm for Phase Unwrapping: http://nova.stanford.edu/sar_group/snaphu/
[48] Near-Real-Time Map of the F2-Layer Height Maximum: http://www.spacew.com/www/hmf2.html
[49] 饒瑞鈞(2017)，2016年高雄美濃地震- 震後科學調查，科技部