最近幾年，台灣的氣候型態呈現明顯的變化，主要表現在極端氣候事件的頻繁發生和降雨型態的劇烈變動。首先，極端氣候事件，如颱風、豪雨和極端高溫，呈現了增加的趨勢，其次，降雨型態的改變也是近年台灣氣候變化的一大特徵，乾旱危機也相繼爆發，造成水資源短缺，尤其對農業灌溉帶來重大影響，因此期望精確評估極端氣候和降雨型態改變對地下 水資源的衝擊，以發展更為有效的水資源管理策略。
本研究選用 MODFLOW 地下水數值模式，對高雄美濃、旗山地區進行地下水模擬，以探討氣候變遷下研究區域的地下水變化情形。本研究的模擬時間為 2015 年至 2023 年，模式再經過率定後採用相關係數及均方根誤差來評估率定結果。擬和結果顯示，相關性高，皆在 0.7 以上;方根誤差 除了吉洋為 1.51m 外，其餘皆在 1m 左右，顯示出本研究的地下水數值模 型能夠準確模擬實際的地下水位變化。
為了分析乾旱對地下水位變化的影響，利用 2020 年的乾旱情況來作為 未來乾旱年模擬的基準，設定初始水位、抽水量及水力傳導係數等參數， 再利用臺灣氣候變遷推估與資訊平台計畫(TCCIP)提供的耦合大氣-海洋- 陸面系統模式(GCMs)模擬未來氣候情境。結果顯示，隨著氣溫上升， 地下水位在不同情境下的變化有所不同，其中 TaiESM1 模型在 4°C 情境 下的地下水位升高，而 EC-Earth3 模型在 2°C 情境下的地下水位升高，顯 示不同模式對地下水補注的預測存在差異。因此，考慮更多的氣候模式及 情境，觀察降雨量的季節性分布及強度變化對地下水補注的影響，有助於 在針對乾旱及極端降雨事件時，能更有效地實施應對策略跟措施，以確保 地下水資源的可持續利用性。

In recent years, Taiwan's climate pattern has shown significant changes, mainly in the frequency of extreme climate events and drastic changes in rainfall patterns. Firstly, extreme climate events, such as typhoons, heavy rainfall, and extreme high temperatures, have shown an increasing trend. Secondly, changes in rainfall patterns have also been a major feature of climate change in Taiwan in recent years. Droughts have erupted one after another, resulting in water scarcity, which has had a significant impact on agricultural irrigation in particular. Therefore, it is hoped that the impacts of the changes in extreme climatic conditions and rainfall patterns on groundwater resources will be assessed accurately, in order to develop more effective water resource management strategies.
In this study, the MODFLOW groundwater numerical model was employed to simulate groundwater changes in the Meinong and Qishan areas of Kaohsiung under climate change conditions. The simulation period spans from 2015 to 2023, and the model calibration was evaluated using correlation coefficients and root-mean-square errors (RMSE). The results showed that the correlation coefficients were all above 0.7, and the RMSE values were around 1 m, except for Jiyang, which had an RMSE of 1.51 m. This demonstrates that the numerical groundwater model can accurately simulate the actual changes in groundwater levels.
In order to analyze the effect of drought on groundwater level changes, the drought in 2020 was used as a benchmark for future drought years, and the initial water level, pumping volume, and hydraulic transmission coefficient were set to simulate future climate scenarios using the coupled atmosphere-ocean-surface system models (GCMs) provided by the Taiwan Climate Change Estimation and Information Platform (TCCIP) program. The results show that the changes in groundwater levels with increasing temperature vary under different scenarios, with the TaiESM1 model showing an increase in groundwater levels under the 4°C scenario and the EC-Earth3 model showing an increase in groundwater levels under the 2°C scenario, which indicates that there is a discrepancy in the prediction of groundwater recharge between different models. Therefore,considering more climate models and scenarios, and observing the impacts of seasonal distribution and intensity changes of rainfall on groundwater recharge can help to implement more effective response strategies and measures in response to droughts and extreme rainfall events, so as to ensure the sustainable utilization of groundwater resources.

1. Arlai, P., Lukjan, A., & Koch, M. (2012). Numerical Investigation of the
Groundwater Balance in the Mae Sai Aquifer, Northern Thailand. Procedia Engineering, 32, 1214-1220. doi:10.1016/j.proeng.2012.02.079.2. Chen, Z., Grasby, S. E., & Osadetz, K. G. J. J. o. H. (2002),Predicting
average annual groundwater levels from climatic variables: an empirical model,260(1-4), 102-117.
3. Eckhardt , K., & Ulbrich, U. J. J. o. H. (2003), Potential impacts of climate change on groundwater recharge and streamflow in a central European low mountain range, 284(1-4), 244-252.
4. Freez,R.A. ( 1979 ) , Groundwater.Prentice Hall,Englewood Cliffs NJ,pp.604.
5. Kirshen, P. H. ( 2002 ) , Potential Impacts of Global Warming on Groundwater in Eastern Massachusetts. Journal of Water Resources Planning
and Management, 128(3), pp.216–226. DOI 10.1061/(asce)0733- 9496 (2002)128:3(216).
6. Patil, N. S., Chetan, N. L., Nataraja, M., & Suthar, S. (2020). Climate
change scenarios and its effect on groundwater level in the Hiranyakeshi watershed. Groundwater for Sustainable Development, 10. doi:10.1016/j.gsd.2019.100323.
7. Ritesh Karki, Puneet Srivastava, Latif Kalin, Subhasis Mitra, Sarmistha Singh.(2021). Assessment of impact in groundwater levels and stream-
aquifer interaction due to increased groundwater withdrawal in the lower Apalachicola- Chattahoochee-Flint (ACF) River Basin using MODFLOW, Journal of Hydrology 100802.
8. Shaban, A., Khawlie M. and Abdallah C.(2006), Use of remote sensing
and GIS to determine recharge potential zone: the case of Occidental Lebanon, Hydrogeology Journal 14: 433-443.
9.Sobeih, M. M., El-Arabi, N. E., Helal, E. E. D. Y., & Awad, B. S. (2017), Management of water resources to control groundwater levels in the southern area of the western Nile delta, Egypt. Water Science, 31(2), pp.137-150. DOI 10.1016/j.wsj.2017.09.001.
10. Taheri Tizro, A., Fryar, A. E., For, M. K., Voudouris, K. S., & Mashhadian,
M. J. (2019). Groundwater conditions related to climate change in the semi-
arid area of western Iran. Groundwater for Sustainable Development, 9. doi:10.1016/j.gsd.2019.100273.
11. Yeh, H.F., Lee, C.H., Hsu, K.C. and Chang, P.H.(2009), GIS for the assessment of groundwater recharge potential zone. Environmental Geology 58, 185- 195.
12. 丁澈士(1997)應用氯離子平衡法推估地下水補注量-屏東平原個案研 究，第二屆地下水資源及水質保護研討會論文集，第 703-714 頁。
13. 田巧玲、張炎銘、簡振源、蔡文豪和陶芳策(1999)，濁水溪沖積扇地 下水水流數值模擬，第三屆地下水資源及水質保護研討會，11~12 頁， 國 立中央大學。
14. 吳雪蘋(2000)，濁水溪沖積扇地區地下水位變化之研究，碩士論文， 國立台灣大學地質學研究所。
15. 吳聲暉(2005)，雪山隧道開鑿滲流量與翡翠水庫進流量相關性之研究， 碩士論文，國立成功大學資源工程學系。
16. 吳偉競(2007)，蘭陽平原整合模式下之 MODFLOW 地下水模擬研究， 碩士論文，國立中正大學地震研究所暨應用地球物理研究所。
17. 吳若穎、胡茵婷、許駿騰(2009)，利用 HEC-HMS 模擬降雨、逕流並 推估 CN 值與 NDVI 之關係，國立屏東科技大學水土保持系。
18. 吳政緯(2012)降雨變異對水文分佈及運移時間影響之研究，碩士論 文，國立成功大學資源工程學系碩博士班。
19. 李心惟(2014)，結合 HEC-RAS 與 MODFLOW 於濁水溪沖積扇地下 水與地層下陷模擬，碩士論文，國立成功大學資源工程學系。
20. 邱永嘉(2002)雲林沿海地區垂直方向地下水補注之研究，碩士論文， 國立臺灣大學地質科學研究所。
21. 林美華(2021)高屏溪流域伏流水資源開發對地下水位及水質影響評 估，碩士論文，國立成功大學水利及海洋工程學系。22. 林聖婷(2012)，濁水溪沖積扇補注量及抽水量空間分佈之模式，碩士 論文，國立臺灣大學土木工程研究所。
23. 林楷文(2023)，美濃、旗山地區農業用水對地下水資源的影響，碩士 論文，國立中央大學土木工程學系。
24. 姜儷安、歐陽湘( 1996 )，雲林地區地下水及水平衡初步分析，濁水 溪沖積扇地下水及水文地質研討會論文集，181~206 頁，國立臺灣大學。
25. 范純志(1998)氣候變遷對臺灣地下水補注量的影響，碩士論文，國 立臺灣大學農業工程學系。
26. 秦啟文、謝文忠、簡傳彬和林癸妙(2000)，嘉義地區地下水安全出水 量，第十一屆水利工程研討會，國立臺灣大學。
27. 張興亞(2001)，灌溉水滲漏之數值模擬，碩士論文，中原大學土木工 程研究所。
28. 國科會(2018)，智慧水管理創新產業計畫-區域地下水智慧管理模式 及技術研發計畫。
29. 陳柏穎(2020)，應用集水區耦合型模式於地下水抽水影響與敏感度之 研究-以屏東平原為例，碩士論文，國立中興大學土木工程學系。
30. 黃智昭(2018)，高屏溪流域上游地下水補注與開發潛能區劃設及可開 發量評估，博士論文，國立成功大學資源工程學系。
31. 黃奕語(2020)不同時空降雨型態對於地下水補注量之探討-以屏東平 原為例，碩士論文，國立中央大學土木工程學系。
32. 曾伃利(2023)以乾旱指標結合概念模式推估地下水補注熱區之研究， 碩士論文，國立中興大學水土保持學系所。
33. 經濟部中央地質調查所(2002)，屏東平原水文地質調查研究總報告。
34. 經濟部中央地質調查所( 2012 )，「臺灣地區地下水區水文地質調查及 地下水資源評估地下水補注潛勢評估與地下水模式建置」。
35. 經濟部中央地質調查所(2014)，地下水補注地質敏感區畫定計畫書 (G0002 屏東平原)。
36. 廖垂麟(2016)地層下陷重點區地下水補注對策之研究，碩士論文， 國立中興大學水土保持學系所。
37. 蔡清研(2007)，濁水溪沖積扇整合模式下之 MODFLOW 地下水模擬 研究，國立中正大學地震研究所暨應用地球物理研究所。
38. 蔡忠遠(2015)探討土地利用和氣候變遷下鳳山溪流域地下水位及流 量的影響，碩士論文，國立中央大學土木工程學系。
39. 賴再壽(2001)，濁水溪沖積扇地下水汲取之研究，碩士論文，國立中 正大學應用地球物理研究所。
40. 賴楷元(2009)，應用河川水位掃瞄評估地下水補注量之研究，碩士論 文，國立成功大學資源工程學系碩博士班。
41. 譚仲哲、童慶斌(2008)氣候變遷對台北地下水補注之衝擊，農業工 程學報，第 54 卷第 1 期，第 1-15 頁。