瞭解地下水對降雨事件的反應有助於瞭解地下水之傳輸，然而，較少利用長時間且連續現地觀測之雨水及地下水資料評估兩者之間的關係。本研究分析(1) 1990至2020年坪林測站歷史雨量資料、2012年5月至2020年5月坪林國小的雨水同位素資料和(2)雪山隧道導坑內2006年7月至2018年6月測量之地下水流量、2012年10月至2019年5月蒐集的地下水同位素資料以及(3)2018年10月至2020年5月於坪林一號豎井附近地下水井之地下水位及其同位素資料，以探討本區雨水－地下水之關連。
研究結果顯示雨水的同位素值分布具明顯季節差異，在導坑及井等不同深度採集的地下水樣本其同位素有明顯空間差異，而所有採樣的地下水同位素值均落在區域天水線上，反映出地下水是雨水的混合。根據前述成果，本研究利用端源混合分析、水文Φ指數及同位素Φ指數評估雨水補注地下水的情形，兩種Φ指數的分析結果顯示分別有31.5% ~ 63.5%和85.2% ~ 91%的年雨量入滲。比較上述三種入滲推估方法中乾季及濕季月份雨水對地下水的相對貢獻，以端源混合分析（乾季19.1% ~ 34.1%、濕季65.9% ~ 80.9%）及同位素Φ指數（乾季20.2% ~ 37.2%、濕季62.8% ~ 79.8%）的分析結果和坪林全年雨量分布（乾季31%、濕季69%）較相似，而根據泰平及坪林氣象站雨量及坪林流量站計算的水文Φ指數則顯示相對平均且不同的貢獻度分布（乾季45.9% ~ 58.7%，濕季41.3% ~ 54.1%）。為進一步探討地下水對降雨事件的反應，首先分析了雨量和井地下水位、雨量和導坑流量的相關性，結果顯示井地下水位和導坑流量對降雨事件反應的延遲時間分別為1.5 ~ 118.8小時和6 ~ 33天，且井地下水位的變化量也和降雨量呈高度正相關。其次，運用週期性回歸分析分別對降雨和地下水中δ18O的季節性趨勢進行建模，擬合結果顯示井地下水δ18O的變化週期約1.9年，導坑地下水為1.9至3.4年，雨水δ18O變化週期則為1年，再根據地下水中δ18O振幅減少的程度計算井地下水的平均停留時間約405至993天；導坑地下水則為342至1691天。
根據上述研究結果，可以觀察到雪山隧道地區地下水對降雨事件的反應特徵相當複雜，透過流量反應（井地下水位、導坑流量）及示蹤劑反應（同位素）所估算之反應時間的尺度相當不同，說明地下水傳輸包含波速(Celerity)與速度(Velocity)兩種概念。研究結果拓展我們對雪山隧道地區地表與地下水關連性的理解，並提供對研究區地下水補注機制的初步想像，期望可作為未來水資源管理的參考。

Understanding the response of groundwater to rainfall events helps to realize groundwater transport process. However, long-term and continuous on-site observation of hydrological data is rarely used to assess the relationship between rainfall and groundwater. This study collected historical rainfall amount data derived from Pinglin Meteorological Station from 1990 to 2020 and isotope data collected in Pinglin Elementary School from May 2012 to May 2020, as well as the deeper groundwater isotope (from October 2012 to May 2019) and discharge (from July 2006 to June 2018) collected from Hsueh-Shan Tunnel. Besides, the shallow groundwater isotope and groundwater level data from the monitoring well near by Pinglin No. 1 shaft from October 2018 to May 2020 were also analyzed to explore the relationship between rainfall and groundwater in Pinglin.
The results demonstrated that the isotope of rainfall had obvious seasonal differences; and the isotope of groundwater collected at distinct depths (i.e. pilot tunnel and well) also had significant difference. All isotopic data of groundwater at each sampling plot approximately fit the local meteoric water line, reflecting that groundwater is the mixture of rainfall. This study applied the end member mixing analysis, hydrological Φ index, and isotopic Φ index to evaluate the amount of rainfall recharging to groundwater. The results demonstrated that there were 31.5% ~ 63.5% and 85.2% ~ 91% of annual rainfall recharging to groundwater, computed by hydrological Φ index and isotopic Φ index, respectively. Comparing with the contribution of rainfall recharge into groundwater in the dry season and wet season, the end member mixing analysis (19.1% ~ 34.1% for the dry season, 65.9% ~ 80.9% for the wet season) and isotopic Φ index (20.2% ~ 37.2% for the dry season and 62.8% ~ 79.8% for the wet season) were similar with the annual rainfall distribution in Pinglin, and hydrological Φ index calculated based on rainfall amount data of Taiping and Pinglin stations showed a relatively average and different distribution of contribution (45.9% ~ 58.7% for the dry season and 41.3% ~ 54.1% for the wet season). To explore the groundwater response to rainfall events, we calculated the correlation between time series of rainfall, groundwater level, and pilot tunnel discharge. The correlation results indicated that the time-lag response of groundwater level and pilot tunnel discharge to rainfall was about 1.5 ~ 118.8 hours and 6 ~ 33 days, respectively; and the water table changes were also highly correlated with rainfall amounts. In addition, the periodic regression analysis with sinusoidal wave is used to fit the periodic trends of δ18O in rainfall and groundwater, indicating that the change period of δ18O in the well, pilot tunnel, and rainfall was about 1.9 years, 1.9 ~ 3.4 years and 1 year, respectively. Based on the damping of the amplitude of δ18O in groundwater, the mean residence time of the well and pilot tunnel was about 405 ~ 993 days and 342 ~ 1691 days, respectively.
It can summarize that the groundwater response to rainfall events in the Hsueh-Shan tunnel was quite complex. The estimated response time of hydrologic responses (groundwater level, pilot tunnel discharge) and tracer response (isotope) were quite different, indicating that groundwater transmission includes two concepts: Celerity and Velocity. The research results expand our understanding of the relationship between surface flow and groundwater in the Hsueh-Shan tunnel area, and provide the preliminary imagination of the groundwater recharge mechanism in the study area, which can be used as a reference for future water resources management.

中文文獻
江昭輝、簡光佑、莊秉潔、黨美齡、李育棋、洪景山、郭珮萱、蔡徵霖（民104）。台灣區域土壤含水率觀測網之建置與資料分析。大氣科學，43(2)，133-150。
江崇榮、汪中和（民91）。以氫氧同位素組成探討屏東平原之地下水補注源。經濟部中央地質調查所彙刊，15，49-67。
李光敦（民94）。水文學。臺北市：五南圖書出版股份有限公司。
佳美檢驗科技股份有限公司（民106）。交通部臺灣區國道高速公路局北區工程處一號豎井附近地下水監測井設置完井報告書。
邱祈榮、梁玉琦、賴彥任、黃名媛（民93）。臺灣地區氣候分區與應用之研究。臺灣地理資訊學刊，1，41-62。
邱琳濱、曾大仁、侯秉承、李民政、張博翔（民97）。雪山隧道營運中地下水監測與分析成果。地下空間與工程學報，4(4)，776-780。
國家災害防救科技中心（民106）。臺灣氣候變遷科學報告2017－物理現象與機制，臺灣氣候變遷推估資訊與調適知識平台。https://tccip.ncdr.nat.gov.tw/publish_01_one.aspx?bid=20171220135820。（2021/01/05瀏覽）
簡銘成、杜永昌、汪中和、丁澈士（民100）。應用氫氧穩定同位素分析地下水補注之研究。農業工程學報，57(3)，61-74。

 
英文文獻
Arnold, J. G., Allen, P. M., Muttiah, R., & Bernhardt, G. (1995). Automated base flow separation and recession analysis techniques. Groundwater, 33(6), 1010-1018.
Beven, K. (1981). Kinematic subsurface stormflow. Water Resources Research, 17(5), 1419-1424.
Bliss, C. I. (1970). Statistics in biology. Vol. 2.
Böhlke, J. K., Jurgens, B. C., Uselmann, D. J., & Eberts, S. M. (2014). Educational webtool illustrating groundwater age effects on contaminant trends in wells. Ground water, 52(Suppl 1), 8-9.
Burgman, J. O., Calles, B., & Westman, F. (1987). Conclusions from a ten year study of oxygen-18 in precipitation and runoff in Sweden. In Isotope techniques in water resources development (pp. 579-590).
Chiu, Y. C., & Chia, Y. (2012). The impact of groundwater discharge to the Hsueh-Shan tunnel on the water resources in northern Taiwan. Hydrogeology journal, 20(8), 1599-1611.
Craig, H. (1961). Isotopic variations in meteoric waters. Science, 133(3465), 1702-1703.
Dansgaard, W. (1964). Stable isotopes in precipitation. Tellus, 16(4), 436-468.
DeWalle, D. R., Edwards, P. J., Swistock, B. R., Aravena, R., & Drimmie, R. J. (1997). Seasonal isotope hydrology of three Appalachian forest catchments. Hydrological Processes, 11(15), 1895-1906.
Evaristo, J., Jasechko, S., & McDonnell, J. J. (2015). Global separation of plant transpiration from groundwater and streamflow. Nature, 525(7567), 91-94.
Gabrielli, C. P., & McDonnell, J. J. (2018). No direct linkage between event‐based runoff generation and groundwater recharge on the MaiMai Hillslope. Water Resources Research, 54(11), 8718-8733.
Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167.
Hou, P. C. & Hsiao, T. H. (2005). Geological exploration of the Hsuehshan tunnel －Review and discussion. In International symposium on design, construction and operation of long tunnels, 105-112.
Hughes, C. E., & Crawford, J. (2012). A new precipitation weighted method for determining the meteoric water line for hydrological applications demonstrated using Australian and global GNIP data. Journal of Hydrology, 464, 344-351.
Kendall, C., & Coplen, T.B. (2001). Distribution of oxygen-18 and deuterium in river waters across the United states. Hydrological Processes, 15, 1363-1393.
Kirchner, J. W. (2003). A double paradox in catchment hydrology and geochemistry. Hydrological Processes, 17(4), 871-874.
Kirchner, J. W., Feng, X., & Neal, C. (2000). Fractal stream chemistry and its implications for contaminant transport in catchments. Nature, 403(6769), 524-527.
Lapworth, D. J., MacDonald, A. M., Krishan, G., Rao, M. S., Gooddy, D. C., & Darling, W. G. (2015). Groundwater recharge and age‐depth profiles of intensively exploited groundwater resources in northwest India. Geophysical Research Letters, 42(18), 7554-7562.
Lis, G., Wassenaar, L. I., & Hendry, M. J. (2008). High-precision laser spectroscopy D/H and 18O/16O measurements of microliter natural water samples. Analytical chemistry, 80(1), 287-293.
Liu, T. K., Cheng-Hong, C., Yang, T. F., & Lee, M. (2005). Tritium concentrations and radiocarbon ages of gushing groundwater from Hsuehshan Tunnel, Northern Taiwan. TAO: Terrestrial, Atmospheric and Oceanic Sciences, 16(4), 909-917.
Lu, H. Y., Peng, T. R., Liu, T. K., Wang, C. H., & Huang, C. C. (2006). Study of stable isotopes for highly deformed aquifers in the Hsinchu-Miaoli area, Taiwan. Environmental Geology, 50(6), 885-898.
McDonnell, J. J. (2014). The two water worlds hypothesis: ecohydrological separation of water between streams and trees?. Wiley Interdisciplinary Reviews: Water, 1(4), 323-329.
McDonnell, J. J., & Beven, K. (2014). Debates—The future of hydrological sciences: A (common) path forward? A call to action aimed at understanding velocities, celerities and residence time distributions of the headwater hydrograph. Water Resources Research, 50(6), 5342-5350.
Meehan, T.D., Giermakowski, J.T., & Cryan, P.M. (2004) GIS-based model of stable hydrogen isotope ratios in North American growing-season precipitation for use in animal movement studies. Isotopes in Environmental and Health Studies, 40(4), 291- 300.
Muñoz‐Villers, L. E., & McDonnell, J. J. (2012). Runoff generation in a steep, tropical montane cloud forest catchment on permeable volcanic substrate. Water Resources Research, 48(9).
Nathan, R. J., & McMahon, T. A. (1990). Evaluation of automated techniques for base flow and recession analyses. Water resources research, 26(7), 1465-1473.
Peng, T. R., Chen, K. Y., Zhan, W. J., Lu, W. C., & Tong, L. T. J. (2015). Use of stable water isotopes to identify hydrological processes of meteoric water in montane catchments. Hydrological Processes, 29(23), 4957-4967.
Peng, T. R., Liu, K. K., Wang, C. H., & Chuang, K. H. (2011). A water isotope approach to assessing moisture recycling in the island‐based precipitation of Taiwan: A case study in the western Pacific. Water Resources Research, 47(8).
Peng, T. R., Wang, C. H., Huang, C. C., Fei, L. Y., Chen, C. T. A., & Hwong, J. L. (2010). Stable isotopic characteristic of Taiwan's precipitation: A case study of western Pacific monsoon region. Earth and Planetary Science Letters, 289(3-4), 357-366.
Quenet, M., Celle-Jeanton, H., Voldoire, O., Albaric, J., Huneau, F., Peiry, J. L., Allain, E., Garreau, A., Beauger, A. (2019). Coupling hydrodynamic, geochemical and isotopic approaches to evaluate oxbow connection degree to the main stream and to adjunct alluvial aquifer. Journal of Hydrology, 577, 123936.
Rodgers, P., Soulsby, C., Waldron, S., & Tetzlaff, D. (2005). Using stable isotope tracers to assess hydrological flow paths, residence times and landscape influences in a nested mesoscale catchment. Hydrology and Earth System Sciences Discussions, 9(3), 139-155.
Rozanski, K., Araguás‐Araguás, L., & Gonfiantini, R. (1993). Isotopic patterns in modern global precipitation. Geophysical Monograph-American Geophysical Union, 78, 1-1.
Scaini, A., Hissler, C., Fenicia, F., Juilleret, J., Iffly, J. F., Pfister, L., & Beven, K. (2018). Hillslope response to sprinkling and natural rainfall using velocity and celerity estimates in a slate-bedrock catchment. Journal of hydrology, 558, 366-379.
Sharp, Z. (2017). Principles of stable isotope geochemistry (2nd ed.). doi: https://doi.org/10.25844/h9q1-0p82
Taylor, R. G., Scanlon, B., Döll, P., Rodell, M., Van Beek, R., Wada, Y., ... & Treidel, H. (2013). Ground water and climate change. Nature climate change, 3(4), 322-329.