屏東平原地表水資源不足，居民大量抽取地下水使用，造成沿海地區發生嚴重之地下水水位大幅下降，進而導致海水入侵與地層下陷問題。本研究整合地下水水質與水量條件劃定屏東平原地下水適合飲用之區域，在地下水水質評估方面將使用多變量指標克利金法根據飲用水水質標準分析適合飲用水質之空間變異，並機率劃定適合飲用之地下水區域；在地下水水量評估方面則以地下水水流模式MODFLOW模擬屏東平原地下水流動與抽水量，模式經過率定後，分析各鄉鎮市抽水與洩降關係，以決定抽水對含水層衝擊性較低之區域；最後結合地下水水質與水量條件評估屏東平原地下水可飲用之區域。本研究結果顯示就水質與水量觀點屏東平原地下水最適合飲用之鄉鎮為鹽埔鄉，而不適合作為飲用水之區域有大寮鄉、東港鄉、南州鄉、崁頂鄉、新園鄉、枋寮鄉、萬巒鄉、佳冬鄉、林邊鄉、林園鄉。

Surface water is considerably lack in the Pingtung Plain.Thus, many residents massively extract groundwater to use, resulting that rapid decreases in groundwater levels, seawater intrusion and land subsidence. This study is to combine water quality variability and flow simulation models to assess potential zones of groundwater used for drinking in the Pingtung Plain. In water quality, multivariate indicator Kriging and water quality standard are used to spatially analyze groundwater quality in the Pingtung plain and to probabilistically delineate zones of proper drinking. In water quantity, the groundwater flow model, MODFLOW, was adopted to simulate groundwater drawdown for specified pumping rates. To reduce impacts on land subsidence by pumping, zones of pristine groundwater quality with minimizing pumping drawdown were the optimal locations to extract groundwater for drinking. The results of this study could offer agricultural administrators a sound guideline for obtaining clean and low-impact groundwater resources for drinking.

Adamsson, M., Dave, G., Forsberg, L., Guterstam, B., 1998. Toxicity identification evaluation of ammonia, nitrite and heavy metals at the Stensund Wastewater Aquaculture plant, Sweden. Water Science and Technology 38(3): 151-157.
AERC (Agriculture Engineering Research Center), 2008. Analysis and Evaluation of the Groundwater Quality Survey in Taiwan Areas. Water Resources Agency, Ministry of Economic Affairs, Executive Yuan, Taiwan, 487pp.
Al-Salamah, L.S., Ghazaw, Y.M., Ghumman, A.R., 2011. Groundwater modeling of Saq Aquifer Buraydah Al Qassim for better water management strategies. Environmental Monitoring and Assessment 173 (1–4), 851–860.
Allan, G.L., Maguire, G.B., 1995. Effect of sediment on growth and acute ammonia toxicity for the school prawn, Metapenaeus macleayi (Haswell). Aquaculture 131(1-2): 59-71.
Appelo, C.A.J., Postma, D., 1993. Geochemistry, Groundwater and Pollution. A.A. Balkema Publishers, VT, USA. pp. 239-295
Bauder, T.A., Waskom, R.M., Davis, J.G., 2007. Irrigation Water Quality Criteria. Available at http://www.ext.colostate.edu/pubs/crops/00506.html.
Bradley, C., 2002. Simulation of the annual water table dynamics of a floodplain wetland, Narborough Bog, UK. Journal of Hydrology 261 (1–4), 150–172.
Chiang, W. H. and W. Kinzelbach, 1998. Processing Modflow for Windows(Version 5.06), Hamburg, Germany.
D''Alessandro, W., Bellomo, S., Parello, F., Brusca, L., Longo, M., 2008. Survey on fluoride, bromide and chloride contents in public drinking water supplies in Sicily (Italy). Environmental Monitoring and Assessment 145(1-3): 303-313.
Deutsch, C.V., Journel, A.G., 1998. GSLIB: Geostatistical Software Library and User’s Guide, second ed. Oxford University Press, New York.
Diodato, N., Ceccarelli, M., 2004. Multivariate indicator kriging approach using a GIS to classify soil degradation for Mediterranean agricultural lands. Ecological Indicators 4(3), 177-187.
Eastoe, C.J., Hutchison, W.R., Hibbs, B.J., Hawley, J., Hogan, J.F., 2010. Interaction of a river with an alluvial basin aquifer: stable isotopes, salinity and water budgets. Journal of Hydrology 395 (1–2), 67–78.
Goovaerts, P., 1997. Geostatistics for Natural Resources Evaluation. Oxford University Press, New York. pp. 259–368.
Grieve, C.M., Poss, J.A., Amrhein, C., 2006. Response of Matthiola incana to irrigation with salinity wastewaters. Hortscience 41(1): 119-123.
Halvorson, J.J., Smith, J.L., Papendick, R.I., 1996. Integration of multiple soil parameters to evaluate soil quality: A field example. Biology and Fertility of Soils 21(3): 207–214.
Harbaugh, A.W., McDonald, M.G., 1996. Programmer documentation for MODFLOW-96, an update to the US Geological Survey modular finitedifference groundwater flow model. USGS Open-File Report, pp. 96–486.
Huang, Y.K., Lin, K.H., Chen, H.W., Chang, C.C., Liu, C.W., Yang, M.H., Hsueh, Y.M., 2003. Arsenic species contents at aquaculture farm and in farmed mouthbreeder (Oreochromis mossambicus) in blackfoot disease hyperendemic areas. Food and Chemical Toxicology 41: 1491-1500.
Huq, S.M.I., Naidu, R., 2005. Arsenic in groundwater and contamination of the food chain: Bangladesh scenario. Natural Arsenic in Groundwater Occurrence, Remediation and Management. The 32nd International Geological Congress, Florence, Italy, 95-101.
Isaaks, E.H.; Srivastava, R.M., 1989. An Introduction to Applied Geostatics. Oxford University Press, New York,. pp278-322.
Jang, C.S., Chen, J.S., 2009. Probabilistic assessment of groundwater mixing with surface water for agricultural utilization. Journal of Hydrology 376: 188–199
Jang, C.S., Chen, S.K., Kuo, Y.M., 2012. Establishing an irrigation management plan of sustainable groundwater based on spatial variability of water quality and quantity. Journal of Hydrology 414-415: 201-210
Jang, C.S., Chen, S.K., Lin, C.C., 2008. Using multiple-variable indicator kriging to assess groundwater quality for irrigation in the aquifers of the Choushui River alluvial fan. Hydrological Processes 22(22): 4477-4489.
Jang, C.S., Liu, C.W., Cuo, Y.L., 2012. Assessment of groundwater emergency utilization in Taipei Basin during drought. Journal of Hydrology 414-415: 405-412
Jang, C.S., Liu, C.W., Lin, K.H., Huang, F.M., Wang, S.W., 2006. Spatial analysis of potential carcinogenic risks associated with ingesting arsenic in aquacultural tilapia (Oreochromis mossambicus) in blackfoot disease hyperendemic areas. Environmental Science & Technology 40: 1707-1713.
Jang, C.S., Liou, Y.T., Liang, C.P., 2010. Probabilistically determining roles of groundwater use in aquacultural fishponds. Journal of Hydrology 388: 491–500
Karaivazoglou, N.A., Papakosta, D.K., Divanidis, S., 2005. Effect of chloride in irrigation water and form of nitrogen fertilizer on Virginia (flue-cured) tobacco. Field Crops Research 92(1): 61-74.
Karaman, M.R., Saltali, K., Ersahin, S., Gulec, H., Rifat, D.M., 2005. Modeling nitrogen uptake and potential nitrate leaching under different irrigation programs in nitrogen-fertilized tomato using the computer program NLEAP. Environmental Monitoring and Assessment 101 (1-3), 249-259.
Knobeloch, L., Salna, B., Hogan, A., Postle, J., Anderson, H., 2000. Blue babies and nitrate-contaminated well water. Environmental Health Perspectives 108(7): 675-678.
Lee, J.J., Liu, C.W., Jang, C.S., Liang, C.P., 2008. Zonal management of multi-purpose use of water from arsenic-affected aquifers by using a multi-variable indicator kriging approach. Journal ofHydrology 359(3-4): 260-273.
Liao, C.M., Jau, S.F., Chen, W.Y., Lin, C.M., Jou, L.J., Liu, C.W., Liao, V.H.C., Chang, F.J., 2008.Acute toxicity and bioaccumulation of arsenic in freshwater clam Corbicula fluminea. Environmental Toxicology 23(6): 702-711.
Liu, C.W., Lin, C.N., Jang, C.S., Chen, C.P., Chang, J.F., Fan, C.C., Lou, K.H., 2006. Sustainable groundwater management in Kinmen Island. Hydrological Processes 20, 4363–4372.
Liu, C.W., Wang, S.W., Jang, C.S., Lin, K.H., 2006. Occurrence of arsenic in groundwater of the Choshui river alluvial fan, Taiwan. Journal of Environmental Quality 35: 68-75.
Liu, C.W., Chou, Y.L., Lin, S.T., Lin, G.J., Jang, C.S., 2010. Management of high groundwater level aquifer in the Taipei Basin. Water Resources Management 24 (13), 3513–3525.
Mittelstet, A.R., Smolen, M.D., Fox, G.A., Adams, D.C., 2011. Comparison of aquifer sustainability under groundwater administrations in Oklahoma and Texas. Journal of the American Water Resources Association 47 (2), 424–431.
Oyedele, D.J., Amusan, A.A., Obi, A.O., 1996. The use of multiple-variable indicator kriging technique for assessment of the suitability of an acid soil for maize. Tropical Agriculture 73(4): 259 –263.
Packham, R.F., 1992. Public health and regulatory aspects of inorganic nitrogen compounds in drinking water. Water Supply 10(3): 1-6.
Rejani, R., Jha, M.K., Panda, S.N., Mull, R., 2008. Simulation modeling for efficient groundwater management in balasore coastal basin, India. Water Resources Management 22 (1), 23–50.
Sedki, A., Ouazar, D., 2011. Simulation-optimization modeling for sustainable groundwater development: a Moroccan coastal aquifer case study. Water Resources Management 25 (11), 2855–2875.
Smith, J.L., Halvorson, J.J., Papendick, R.I., 1993. Using multiple-variable indicator kriging for evaluating soil quality. Soil Science Society of America Journal 57: 743 –749.
Taiwan Sugar Company, 2005. Groundwater Quality Monitoring and Analysis for Taiwan Groundwater Monitoring Network (2/2). Water Resources Agency, Ministry of Economic Affairs, Executive Yuan, Taiwan, pp.10-1–10-74.
Taiwan WRA (Water Resources Agency), 2004. The Inquiry and Supply System of Water Resources Data. Water Resources Agency, Ministry of Economic Affairs, Executive Yuan, Taiwan. <http://gweb.wra.gov.tw/wrwebeng/index.htm>.
Tsai, J.W., Liao, C.M., Liao, V.H.C., 2006. A biologically based damage assessment model to enhance aquacultural water quality management. Aquaculture 251(2-4): 280-294.
Tseng, W.P., 1977. Effects and dose-response relationships of skin cancer and blackfoot disease with arsenic. Environ Health Perspect 19:109-119
Wang, S.W., Liu, C.W., Jang, C.S., 2007. Factors responsible for high arsenic concentrations in two groundwater catchments at Taiwan. Applied Geochemistry 22: 460–476.
Yevtushenko, N.Y., 1998. Accumulation of microelements in organs and tissues of fishes differing in feeding specialization under conditions of cultivation in fishponds with heated water. Hydrobiological Journal 34(4-5): 112-124.
Zabin, S., Foaad, M.A., Al-Ghamdi, A., 2008. Non-carcinogenic risk assessment of heavy metals and fluoride in some water wells in the Al-Baha region, Saudi Arabia. Human and Ecological Risk Assessment 14(6): 1306-1317.
Zaidel, J., Markham, B., Bleiker, D., 2010. Simulating seepage into mine shafts and tunnels with MODFLOW. Ground Water 48 (3), 390–400.
中央地質調查所（2001）台灣地區地下水觀測網第一期計畫-屏東平原水文地質調查研究總報告。
吳銘志、鄭鉅霖（1997）屏東沖平原內自流井之分佈及其賦存層之地層組構。第二屆地下水資源及水質保護研討會論文集，國立成功大學，第947-954頁。
徐鐵良（1961）臺灣南部屏東各地之自升地下水系。中國地質學會會刊，4卷，第73-81頁。
張良正、蔡威平、陳字文，1999，屏東地區地下水補注量推估及分級，第三屆地下水資源及水質保護研討會論文集，pp. 65-76。
陸挽中，陳瑞娥，黃智昭，江崇榮，費立沅（2007），屏東地下水文地質與地下水補注，地質環境與資源研討會論文集，經濟部中央地質調查所，第141-148頁。
經濟部水資源統一規劃委員會（1995）臺灣地區之水資源，共334頁。
經濟部水資源局（2000）中華民國八十八年臺灣水文年報。行政院經濟部水資源局，共702頁。
環境保護署（1998）環境白皮書。八十七年版，行政院環境保護署，共517頁。
經濟部水利署（2005）台灣地下水資源情勢之研究。