Региональная оценка ресурсов подземных вод и их защищенности от загрязнения на Равнине Пинг-Тонг, Тайвань тема диссертации и автореферата по ВАК РФ 25.00.36, кандидат геолого-минералогических наук Цао Хао-Вэй

• Цао Хао-Вэй
• кандидат геолого-минералогических наук
• 2005

Оглавление диссертации кандидат геолого-минералогических наук Цао Хао-Вэй

Введение1.

Глава 1. Современные методы качественной оценки уязвимости подземных вод к загрязнению. 1.1 Общие положения; основные понятия и определения оценки уязвимости подземных вод.5.

1.2 Анализ методов оценки уязвимости подземных вод к загрязнению13.

Глава 2. Природные условия Равнины Пинг-Тонг

2.1 Физико-географические условия.35.

2.2 Основные водоносные горизонты.42.

Глава 3. Характеристики подземного стока и естественных ресурсов подземных

3.1 Количественная оценка модулей и коэффициентов естественных ресурсов46.

3.2 Ориентировочная оценка водного баланса на территории Равнины Пинг-Тонг58.

Глава 4. Разработка комбинированной модели оценки уязвимости подземных вод и ее применения на территории Пинг-Тонг.

4.1 Общее положение и специфика загрязнения нитратов74.

4.2 Оценка уязвимости подземных вод по методу DRASTIC

4.3 Оценка времени проникновения ЗВ до уровня ГВ

4.4 Оценка уязвимости подземных вод по комбинированному методу

Введение диссертации (часть автореферата) на тему «Региональная оценка ресурсов подземных вод и их защищенности от загрязнения на Равнине Пинг-Тонг, Тайвань»

Известно, что подземные воды играют большую роль в водообеспечении коммунального хозяйства. Чтобы с достаточной точностью оценить и определить перспективы использования подземных вод, необходимо проводить основательное гидрогеологическое исследование в интересующих нас районах. Это предполагает моделирование подземных потоков, питания и расходов в целом.

На Тайване в настоящее время увеличивается потребность в изучении подземных вод, так как интенсивная их эксплуатация влияет на окружающую среду. Однако только в последнее время начались исследования по оценке подземных вод в региональном масштабе. В предполагаемой диссертации анализируется район «Равнина Пинг-Тонг», где недавно были проведены исследования для определения перспектив управления подземными водами. В этих исследованиях проводились расчеты питания подземных вод с использованием методов баланса почвенной влагоемкости и массо-равновесия хлорида, с учетом различных элементов уравнения водного баланса. Кроме того, учитываются фильтрации из искусственных водоемов.

В данной диссертации проведен всесторонний анализ распространения подземных вод и стоков в районе исследуемой территории с помощью методов, разработанных в России. Эти методы делятся на четыре основные группы: гидродинамические, балансовые, гидрометрические, аналоговые. Причем каждый из них дополняет друг друга. Помимо этого, рассматривается уязвимость подземных вод к загрязнению в региональных масштабах. Таким образом, достигается практическое сравнение настоящих методов оценки с существующими. Важнейшие задачи, которые должны быть решены - это исследование интенсивности испарения и искусственного питания подземных вод от орошения.

В результате исследования предполагается усовершенствование методики региональной оценки естественных ресурсов и защищенности подземных вод применительно к гидрогеологическим и геологическим условиям Равнины Пинг-Тонг, и выявление участков, перспективных для проведения поисково-разведочных работ.

Цель и задачи работы. Целью данной работы является оценка защищенности грунтовых вод Равнины Пинг-Тонг от загрязнения.

Для достижения поставленной цели были сформулированы следующие задачи:

1. Охарактеризовать природные условия формирования ресурсов подземных вод.

2. Изучить природные факторы защищенности грунтовых вод Равнины Пинг-Тонг.

3. Усовершенствовать методику оценки уязвимости подземных вод применительно к природным условиям Равнины Пинг-Тонг.

4. Провести качественную и количественную оценку уязвимости подземных вод на Равнине Пинг-Тонг.

5. Выделить районы, перспективные для использования подземных вод с точки зрения обеспеченности ресурсами и защищенности от загрязнения.

Фактический материал и методика исследования.

В основу данной работы положены результаты исследований совместного Тайваньско-российского проекта «Карты ресурсов подземных вод Тайваня». При выполнении работы использовались геологические и гидрогеологические данные, собранные Гидравлической Лабораторией Тайнана - Университет Ченг-гонг. Данные по геологическим разрезам получены по результатам проекта «Сеть мониторинга подземных вод». Основным объектом исследований были грунтовые воды водоносного горизонта Равнины Пинг-Тонг, сложенного песчаниками и гравием.

В работе также использованы литературные данные.

Методика обработки и интерпретации информации основана на анализе и использовании существующих методов оценки уязвимости подземных вод к загрязнению. Фактический материал обработан автором, в результате чего проведена оценка уязвимости подземных вод к загрязнению и построены карты защищенности подземных вод Равнины

Пинг-Тонг.

Научная новизна работы.

1. Проведена оценка уязвимости подземных вод к загрязнению различными методами; выполнены анализ достоинств и недостатков этих методов для региональных оценок.

2. Предложен новый способ оценки влияния зоны аэрации, и ее учет при оценке уязвимости методом DRASTIC.

3. Обоснован комбинированный подход к оценке уязвимости, основанный на семи факторах DRASTIC.

4. Выявлены особенности природной защищенности подземных вод Равнины Пинг-Тонг в зависимости от характеристик природных условий.

Практическая значимость работы.

Результаты выполненных исследований могут быть использованы для:

- рационального выбора безопасных мест для размещения водозаборных сооружений;

- обоснования планов размещения и развития крупных промышленных проектов по захоронению опасных отходов;

- обоснования базы данных для будущих исследований;

- гидрогеологического обоснования различных водоохранных мероприятий.

Защищаемые положения.

1. Метод оценки уязвимости подземных вод к загрязнению, применительно к исследуемой территории.

2. Выявление преимуществ и ограничений существующих методов региональной оценки защищенности подземных вод от загрязнения.

3. Количественная оценка времени проникновения загрязняющих веществ в подземные воды Равнины Пинг-Тонг

4. Рекомендации по учету защищенности подземных вод от загрязнения при обосновании водозаборных мероприятий и оценке перспектив использования подземных вод Равнины Пинг-Тонг.

Апробация работы

Основные положения и результаты диссертационной работы докладывались на научных семинарах Института водных проблем РАН, на семинарах Гидравлической лаборатории Тайнана (на Тайване) и на следующих конференциях: Groundwater Vulnerability Assessment and Mapping (Ustron, Poland 2004), Экватек 2004 (Москва, 2004).

По теме диссертации опубликованы 3 работы.

Структура работы

Диссертация содержит 118 страниц текста, состоит из введения, 4-х глав, заключения и выключает 14 таблиц, 37 рисунков и список литературы, состоящий из 103 наименований.

Заключение диссертации по теме «Геоэкология», Цао Хао-Вэй

Заключение

Карта уязвимости подземных вод к загрязнению показывает нам чувствительность подземных вод к различным условиям их залегания. Эта карта может быть использована для разработки стратегий использования и защиты подземных вод в районах с различной природной уязвимостью. Сочетание этой карты с картой ресурсов подземных вод позволяет определить зоны, которые имеют значительные ресурсы подземных вод и в тоже время требуют более детальных исследований по их защищенности от загрязнения. В результате такого процесса обосновывается конкретный выбор зоны, где необходимо проводить 2-ой этап оценки уязвимости подземных вод. 2-ой этап оценки уязвимости подземных вод заключается в оценке времени, в течение которого загрязнитель достигает уровня грунтовых вод.

На Равнине Пинг-Тонг при сравнении карты уязвимости и ресурсов подземных вод видно, что зоны, где имеются обильные ресурсы подземных, являются уязвимыми к загрязнению. Зоны, которые не имеют большие ресурсы подземных вод, наоборот не подвержены воздействию загрязнения. Поскольку факторы, которые определяют режим ресурсов и уязвимости подземных вод на Равнине Пинг-Тонг имеют одинаковые закономерности, например в зоне илистых пород степень уязвимости - низкая и ресурсы подземных вод - небольшие, а в зоне крупнообломочных пород степень уязвимости -высокая и ресурсы подземных вод - большие, можно сделать вывод, что при строении водозаборных строений в предгорных районах на Равнине Пинг-Тонг необходимо проводить детальные расследования времени достижения различных видов ЗВ.

В зависимости от источников загрязнения оценки времени добегания ЗВ до уровня ГВ можно использовать для обоснования водозащитных мероприятий и дальнейших мероприятий по улучшению экологических условий Равнины Пинг-Тонг. При этом выделяются участки, которые являются перспективными для эксплуатации и уязвимыми к загрязнению. Эти участки показываются на карте (рис. 1). В основном эти участки имеют мощность зоны аэрации от 20-30 метров и сложены крупнообломочными породами. Наличие этих хорошо фильтрирующих пород и отсутствие глин означает, что там есть хорошее условие для формирования ресурсов подземных вод. Модуль подземного стока в этих участках оценивается в пределах 20 - 40 л/с*км2. Эти участки расположены в аллювиальных конусах выноса равнины недалеко от сброса Цао-цо. На Равнине Пинг-Тонг не существуют участки, где имеются большие ресурсы подземных вод и низкая степень уязвимости к загрязнению. Где подземные воды защищенные, ресурсы подземных вод небольшие.

В целом анализ влияния различных природных факторов на защитные свойства безнапорных водоносных горизонтов и формирование ресурсов подземных вод в безнапорном горизонте Равнины Пинг-Тонг позволяет делать вывод, что чем благоприятнее условие для формирования ресурсов подземных вод, чем выше уязвимости подземных вод к загрязнению.

На участках, где уязвимость и ресурсы подземных вод большие, рекомендируется проводить второй этап исследования. При более детальных исследованиях и сборе данных о сорбционных свойствах зоны аэрации зоны и водоносного горизонта и о миграционных параметрах, можно рассчитать время проникновения ЗВ. Исследование влияния загрязнителей зависит также от расположения конкретных загрязняющих объектов. В настоящее время существуют многие математические модели и компьютерные программы, например Modflow, СНЕМ, PATH, HydrogeoCHEM, которые могут быт использованы для этой задачи. участки с самой высокой уязвимостью и большими ресурсами подземных вод

Е2 участки с высокой уязвимостью и большими ресурсами подземных вод участки с самой высокой уязвимостью и минимальными ресурсами подземных вод

ЕЕ] участки с высокой уязвимостью и минимальными ресурсами подземных вод □ участки со средней уязвимостью и минимальными ресурсами подземных вод

ZD участки с низкой уязвимостью и минимальными ресурсами подземных вод О

О с=иг=

I 21,200 м

Рис. 1. Карта, где отражаются ресурсы и уязвимости подземных вод Равнины Пинг-Тонг.

В заключении нужно отметить, что региональная оценка и картирование защищенности подземных вод основаны на анализе и обработке только имеющихся геологических, гидрологических и гидрогеологических данных без проведения специальных дорогостоящих буровых, опытно-фильтрационных или лабораторных работ. Выбор метода оценки уязвимости подземных вод при этих условиях во многом определяется количеством имеющихся фактических данных. На Равнине Пинг-Тонг имеются детальные данные по ее геологическим строению и разрезам и на основе этого использовались методы, предложенные в данной работе.

Решение применить комбинированный метод для оценки уязвимости подземных вод основано на том, что качественные и количественные подходы к оценке уязвимости имеют свои недостатки и совместный подход может дать более обоснованную оценку уязвимости подземных вод. Результаты количественных методов, используемых в работе, показывают, что в зависимости от использованных уравнений расчеты времени проникновения ЗВ значительно отличаются. Поэтому нельзя использовать чисто результаты, полученные из математических уравнений, которые предполагают многие допущения в региональной оценке уязвимости подземных вод. Качественный метод оценки уязвимости подземных вод позволяет определить какую территорию лучше защищена от загрязнения на основе баллов, которые характеризуют степень влияния каждого фактора на защищенности данной территории. Этот подход - обычно применяется для оценки уязвимости в региональном масштабе, но часто вызывает критику, посязшенную выборам бальных значений каждого фактора.

Карта уязвимости подземных вод позволяет сравнить различные территории в отношении защищенности подземных вод от загрязнения и решить, где существует большая опасность для загрязнения скважин, эксплуатирующих подземную воду для водоснабжения, где в первую очередь необходимы водоохранные мероприятия. При совместном использовании карты уязвимости и ресурсов подземных вод можно выделить участки, на которых необходимо проводить более детальные полевые исследования. Этот процесс помогает специалистам решить какие районы требуют больше финансовым затратам для водоохранных и водозаборных мероприятий.

Результаты оценки уязвимости подземных вод на Равнине Пинг-Тонг показывают, что закономерности формирования аллювиальных конусов выноса равнины играют большую роль в определении условий защищенности этой территории. Эти закономерности сильно отражаются на карте уязвимости подземных вод Равнины Пинг-Тонг. Уязвимость подземных вод на Равнине постепенно уменьшается от предгорных районов до прибережных зон.

При сочетании карты уязвимости подземных вод с картой модули подземного стока можно делать следующее выводы:

1. Питание подземных вод в предгорных зонах значительно влияет на степень уязвимости подземных вод.

2. Закономерности в формировании аллювиальных конусов выноса исследуемой территории определяют и интенсивность питания подземных вод, и степень уязвимости подземных вод.

3. Необходимо организировать систему мониторинга за ГВ в предгорных зонах Равнины Пинг-Тонг. Необходимо проводить специальные геофильтрационные и геомиграционные исследования с целью определения параметров фильтрации и миграции ЗВ.

4. В дальнейших рекомендируется провести оценку уязвимости подземных вод по времени достижения ПДК для нитратов, как практически не сортируемого загрязнителя.

Список литературы диссертационного исследования кандидат геолого-минералогических наук Цао Хао-Вэй, 2005 год

1. Арипов С.А. О миграции пестицидов в зоне аэрации и влияние их на загрязнение грунтовых вод//Гидрогеология ноосферы. Ташкент, 1976. Выпуск. 1. с. 42-48.

2. Белоусова А.П. Качество подземных вод, современные подходы к оценке. М.: Наука 2001 с. 339

3. Белоусова А.П. Основные принципы и рекомендации по оценке и картированию защищенности подземных вод от загрязнения. Водные ресурсы, 2003, том 33, №6, с. 667-677.

4. By М.Ч. , Чеиг Ч.Ш., Кау Р.Ч., Цао Х.В., Ценг Ч.М., JIu Ю.П. Оценка уязвимости подземных вод в регионах на основе составной модели пример региона Пинг-Тонг, Тайвань. // Шестой международный конгресс ВОДА: Экология и Технология ЭКВАТЭК-2004. С. 223.

5. Гольдберг В.М. Оценка условий защищенности подземных вод и построение карт защищенности // Гидрогеологические основы охраны подземных вод. М.: Недра, 1984. с. 171-177.

6. Данило Ж. Маркович. Социальная Экология. Изд-во «Российский Университет Дружбы народов», 1997.

7. Данилов-Данильян В.И., Лосев КС. Экологический вызов и устойчивое развитие. Изд-во, 2000.

8. Зекцер И.С. Подземные воды как компонент окружающей среды. М.: Научный мир 2001 с. 327.

9. Зекцер И.С. Сколько воды под землей? Изд-во «Знание» 1986.

10. Зекцер И.С., О.А. Каримова, Бужуоли Ж., Буччи М., Региональная оценка уязвимости пресных подземных вод: методологические аспекты и пратическое применение. Водные ресурсы, 2004, том 31, №6, с. 645-650.

11. Зелинский И.Л. и др. К методике оценки защищенности подземных вод от загрязнения // Тез. Докл. Всесоюзного совещания по подземным водам Востока СССР. Иркутск, 1988. с. 171-172.

12. Карагодина Н.В. Геологический аспект составлении региональных схем охраны природы // Тез. Докл. Первого Всесоюзного съезда инженеров-геологов, гидрогеологов и геокриологов. Киев, 1988.4.5. с. 62-63.

13. Мироненко В.А., Румынии В.Г. Оценка защитных свойств зоны аэрации // Инженерная геология. 1990. No. 2. с. 3-18.

14. Рогачевская Л.М. Современное состояние региональных исследований защищенности и уязвимости подземных вод к радионуклидному загрязнению // ВИНИТИ проблемы окружающей среды и природных ресурсов No.5,2001.

15. Шаманже Ж. Оценка защищенности грунтовых вод бассейна р. Бенуэ от загрязнения нитратами соединениями. Проблемы окружающей среды и природных ресурсов. Обзорная информация. ВИНИТИ, выпуск No. 9, 2004 г. с. 80-88.

16. Шестаков В.М., Поздняков С.П. Геогидрология // ИКЦ «АКАДЕМКНИГА» 2003.

17. Штритер Е.Е. и Ахметьева Н.П. Protection of groundwater against nitrogen compound // ICHE-98, Германия, 1998, с. 112. (SD. 112-118).

18. Al-Zabet T. Evaluation of aquifer vulnerability to contamination potential using the DRASTIC method. DOI 10.1007/s00254-002-0645-5 Environmental Geology (2002) 43:203-208.

19. Anderson L.J. and E. Gosk. Applicability of vulnerability maps // Vulnerability of Soil and groundwater to pollutants. 1987 Proceedings and Information no. 38. p. 321-332.

20. Artuso, E, Oliveira, M.M., & Lobo Ferreira, J.P. (2002) "Assessment of Groundwater Vulnerability to Pollution using six different Methods: AVI, GOD, DRASTIC, SI, EPPNA and SINTACS. Application to the Evora Aquifer", LNEC, Rel. 184/02-GIAS, June 2002.

21. Arya L.M., Leij F.J., Shouse P. J., and Genuchten M. Relationship between the Hydraulic Conductivity Function and the Particle-Size Distribution. Soil Sci. Soc. Am. J. (1999) 63:1063-1070.

22. Bachmat Y., Collin M. Mapping to assess groundwater vulnerability to pollution // Vulnerability of Soil and groundwater to pollutants. 1987 Proceedings and Information no. 38. p. 297-306.

23. Bishop G.D., Church M.R., Aber J.D., Neilson R.P., Ollinger S.V., Daly C. A comparison of mapped estimates of long-term runoff in the northeast United States Journal of Hydrology 206(1998) 176-190.

24. Chang, C.C. Water Management and the Outlook for Water Adequacy in Taiwan's Agricultural Sector. Industry of Free China, Nov. 2001. p. 47-58.

25. Chen C. and Payne W.A. Measured and Modelled Unsaturated Hydraulic Conductivity of a Walla Walla Silt Loam Soil Sci. Soc. Am. J. (2001) 65:1385-1391.

26. Chen CJ, et al. Malignant neoplasms among residents of a blackfoot disease-endemic area in Taiwan: high-arsenic artesian well water and cancers. Cancer research, 1985, 45:5895-5899.

27. Chen X., Ни Q. Groundwater influences on soil moisture and surface evaporation Journal of Hydrology 297 (2004) 285-300

28. Chiang, C.J., Huang, C.C., Chen, J.E., The groundwater recharge areas of Pingtung Plain. Bulletin of the Central Geological Survey, (2002) No. 15, p. 17-47. (in Chinese)

29. Church M.R., Bishop G.D., Cassell D.L. Maps of regional evapotranspiration and runoff precipitation ratios in the northeast United States Journal of Hydrology 168 (1995) 283-298.

30. Civita M. (2003). Assessing and mapping groundwater vulnerability to contamination the Italian "combined" approach, Aquifer Vulnerability and Risk. Proceedings of the First International Workshop, May 28-30,2003, Salamanca, Mexico. Vol. 1. pp. 3-19.

31. Civita, M. (1994) Le carte della vulnerability degli acquiferi all'inquinamento. Studi sulla vulnerability degli acquiferi. Teoria & Pratica. Pitagora Editrice, Bologna, 325 pp.

32. Corniello A., Ducci D. and Monti G.M. Aquifer pollution vulnerability in the Sorrento peninsula, southern Italy, evaluated by SINTACS method Geofisica Internacional (2004), Vol. 43, Num. 4, pp. 575-581.

33. Doerfliger N. Jeannin P.J. and Zwahlen F. Water vulnerability assessment in karst environments: a new method of defining protection areas using a multi-attribute approach and GIS tools (EPIK method). Environmental Geology (1999) 39(2): 165-176.

34. Droogers P. & Allen R.G. Estimating reference evapotranspiration under inaccurate data conditions Irrigation and Drainage Systems 16: 33-45, 2002.

35. Fetter C. W. Applied Hydrogeology 2001 Prentice-Hall Inc., Upper Saddle River, N.J. pp.598

36. Foster S. Aquifer protection. Are the vulnerability maps a useful tool? Aquifer Vulnerability and Risk. Proceedings of the First International Workshop, May 28-30, 2003, Salamanca, Mexico. Vol. 1. pp. 325-338.

37. Freeze R.A. and Cherry J.A. Groundwater 1979 Prentice-Hall Inc., Englewood Cliffs, NJ. pp.604

38. Genuchten, M.Th. Van. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. (1980) 44:892-898.

39. Gogu R.C. and Dassargues A. Current trends and future challenges in groundwater vulnerability assessment using overlay and index methods. Environmental Geology 39 (6) April 2000

40. Gogu R.C. and Dassargues A. Sensitivity analysis for the EPIK method of vulnerability assessment in a small karstic aquifer, southern Belgium. Hydrogeology Journal (2000) 8A 337-345.

41. Groundwater Resources Map of Taiwan Modulus of Groundwater Flow- Published by the Water Resources Agency, Ministry of Economic Affairs, in Taiwan 2003.

42. Holmes R.M. and Robertson G.W. A modulated soil moisture budget. March monthly weather review. 1959 pp. 101-105

43. Hsu S.K. Plan for a groundwater monitoring network in Taiwan. Hydrogeology Journal (1998) 6:405-415

44. Kim Y.J. & Hamm S. Y. Assessment of the potential for groundwater contamination using the DRASTIC/EGIS technique, Cheongju area, South Korea. Hydrogeology Journal (1999) 7:227-235.

45. Kondoh, A., Harto A.B., Eleonora R. and Kojiri T. Hydrological regions in monsoon Asia Hydrol. Process. (2004) vol. 18, 3147-3158.

46. Kuo, C.H. and Wang, C.H. The implication of the delineation on the hydrogeologic framework of the shallow Pingtung Plain aquifer. Western Pacific Earth Sciences, 2001 V.l, No.4, pp 459-472.

47. Le Grand, H. 1983. A standardized system for evaluating waste disposal sites. NWWA, Worthington, OH: 49 pp.

48. Leal J.A.R. and Castillo R.R. Aquifer vulnerability mapping in the Turbio river valley, Mexico: A validation study. Geoflsica Internacional (2003), Vol. 42, Num. 1, pp. 141-156.

49. Li G.C., Wong S.S. and Tsai M.C. Safety Evaluation and Regulatory Control of Pesticide Residues in Taiwan Journal of Food and Drug Analysis, Vol. 10, No. 4, 2002, Pages 269-277.

50. Lyon S. W., Walter M. Т., G'erard-Marchant P. andSteenhuis T.S. Using a topographic index to distribute variable source area runoff predicted with the SCS curve-number equation Hydrol. Process. 18, 2757-2771 (2004) DOI: 10.1002/hyp.l494.

51. Maidment, David R. and Reed, Seann M. Water Balance of Africa Exercises on GIS in Water Resources for Africa FAO/UNESCO.

52. Margat J., Suais-Parascandola M.F. Mapping the vulnerability of groundwater to pollution. Some lessons of experience in France // Vulnerability of Soil and groundwater to pollutants. 1987 Proceedings and Information no. 38. p .433-436.

53. Махе L. and Johansson P. Assessing groundwater vulnerability using travel time and specific surface area as indicators. Hydrogeology Journal (1998) 6:441—449

54. Methods of Rating Unsaturated Zone and Watershed Characteristics of Public Water Supplies in North Carolina U.S. GEOLOGICAL SURVEY Water-Resources Investigations Report 99—4283

55. Moges, S.A., Katambara, Z. and Bashar, K. Decision Support System for Estimation of Potential Evapotranspiration in Pangani Basin 3rd WaterNet/Warfsa Symposium 'Water Demand Management for Sustainable Development1, Dar es Salaam, 30-31 October 2002.

56. NRC, 1993. Groundwater vulnerability assessment; contamination potential under conditions of uncertainty. National Academy Press (Washington DC-USA) 198pp.

57. Poianil K.A., Bedford B.L. and Merrid M.D. A GIS-based index for relating landscape characteristics to potential nitrogen leaching to wetlands Landscape Ecology vol. 11 no. 4 pp 237-255 (1996)

58. Scanlon B.R. and Goldsmith R.S. Field study of spatial variability in unsaturated flow beneath and adjacent to playas. Water Resources Research 1997, vol. 33, No. 10, pp. 22392252.

59. Scanlon BR, Healy RW, Cook PG. Choosing appropriate techniques for quantifying groundwater recharge. Hydrogeo J.(2002) DOI 10.1007/s 10040-001-0176-2.

60. Scanlon, B. R., and Goldsmith, R. S., 1997, Field study of spatial variability in unsaturated flow beneath and adjacent to playas: Water Resources Research, v. 33, no. 10, p. 2239-2252.

61. Schaap M.G. andLeij F.J. Improved Prediction of Unsaturated Hydraulic Conductivity with the Mualem-van Genuchten Model. Soil Sci. Soc. Am. J. (2000) 64:843-851.

62. Schoups G. and Hopmans J. W. Analytical Model for Vadose Zone Solute Transport with Root Water and Solute Uptake Vadose Zone Journal 1:158-171 (2002).

63. Sophocleous M. From safe yield to sustainable development of water resources—the Kansas experience. Journal of Hydrology 235 (2000) 27-43.

64. Sophocleous, Marios. Interactions between groundwater and surface water: the state of science. Hydrogeology Journal (2002) 10:52-67.

65. Snow J. On the mode of communication of cholera// London Medical Gazette XLIV. — 1849.-p. 730-732.

66. Sun H. Estimating The Fate of Precipitation From Stream Discharge:A Case Study in New Jersey Bull. N.J. Acad. Sci. (2004), 49(2), pp. 9-15.

67. Taiwan Groundwater Observation Network-Final Report for first research term (hydrogeological conditions of Pingtung Plain). Central Geological Survey, 2002.(in Chinese)

68. Thornthwaite, C. W. An approach toward a rational classification of climate. Geographic Review (1948) 38:55-94

69. Ting C.S., Kerh T.F., Liao C.J. Estimation of groundwater recharge using the chloride mass-balance method, Pingtung Plain, Taiwan. Hydrogeology Journal (1998) 6:282-292

70. Ting, Cheh-Shyh. 1997. Groundwater Resources Evaluation and Management for Pingtung Plain, Taiwan. Amsterdam: Free University, pp. 199.

71. Tsai C.C., Chen Z.S., Duh C.T. and Horng F.W. Prediction of Soil Depth Using a Soil-Landscape Regression Model: a Case Study on Forest Soils in Southern Taiwan. Proc. Natl. Sci. Counc. ROC(B) Vol. 25, No. 1, 2001. pp. 34-39

72. Van Stempvoort, D., Ewert, L. & Wassenaar, L. (1993), Aquifer Vulnerability Index: A GIS compatible method for groundwater vulnerability mapping. Canadian Water Resources Journal. Vo. 18, No. 1, pp.25-37.

73. Voigt H.J., Heinkele Т., Jahnke C. and Wolter R. Characterization of groundwater vulnerability to fulfill requirements of the water framework directive of the European Union. Geofisica Internacional (2004), Vol. 43, Num. 4, pp. 567-574

74. Voigt, Hans-Jiirgen; Heinkele, Thomas; Jahnake, Christoph; Wolter, Riidiger. (2003). Characterization of Groundwater Vulnerability. Aquifer Vulnerability and Risk. First International Workshop, May 28-30, 2003, Salamanca, Mexico. Vol. 2. pp. 189-201.

75. Vorosmarty C.J., Federer C.A., Schloss A.L., Potential evaporation functions compared on US watersheds: Possible implications for the global scale water balance and terrestrial ecosystem modeling. Journal of Hydrology 207 (1998) pp/ 147 169.

76. Vrba J., Zaporozec A. Guidebook on Mapping Groundwater Vulnerability. Int. Association of Hydrogeologists.1994. Hannover, Heise. Vol.16

77. Wang C. Deterioration of groundwater quality in the coastal Pingtung Plain, southern Taiwan // Impacts of Urban Growth on Surface Water and Groundwater Quality. 1999. Proceedings of IUGG 99 Symposium HS5. IAHS Publ. No. 259, 1999.

78. Wang C.H. and Liu C. Dissipation of Organochlorine Insecticide Residues in the Environment of Taiwan, 1973-1999 Journal of Food and Drug Analysis, (2000) Vol. 8, No. 3.

79. Water Resources Bureau, (1998) Hydrologic yearbook of Taiwan. Ministry of Economic Affairs, Taipei, Taiwan.

80. Welsch D.L., Kroll C. N., McDonnell J.J. and Burns D.A. Topographic controls on the chemistry of subsurface stormflow Hydrol. Process. 15, 1925-1938 (2001) DOI: 10.1002/hyp.247

81. Wosten, Henk; Nemes, Attila; Pachepsky, Yakov; Rawls, Walter and Zeiliguer, Anatole. Translating available basic soil data into missing soil hydraulic characteristics. 17th WCSS, 14-21 August 2002 Thailand.

82. Wu M.C., King S.H., Zektser I.S., Dzhamalov R.G., Kao R.C., Cheng C.L, Chen C.K. & Chen C.K, 2002, Inventory of the regional natural groundwater resources in Pingtung region of southern Taiwan, IAH Australia Conference, Юр.

83. Wu M.C.; Cheng C.L. & Shih KM. 1996, Paleofluminolobic Settings and Artesian Aquifer Within Pingtung Alluvial Plain A Preliminary Study, 1996 Annual Meeting, Geological Society of China, pp. 93-97. (In Chinese)

84. Wu M.C.; Cheng C.L. Chen C.K. & Chen C.K, 2001, Preliminary Discussion for Regional Lateral Recharge Assessment of Pingtung Plain, The 4th Conference for Protection of Groundwater Resources and Water Quality, pp. 217-222. (In Chinese)

85. Zektser Igor S., Pozdniakov Sergey P., Szpakiewicz Michael andRogachevskaya Liliya M. Regional assessment of groundwater vulnerability in the Snake river plain aquifer basin, USA Geofisica Intemacional (2004), Vol. 43, Num. 4, pp. 697-705.

86. Zhang R„ Hamerlinck J.D., Gloss S.P. and Munn L. Determination of Nonpoint-Source Pollution Using GIS and Numerical Models Journal of Environmental Quality (1996) 25(3):411-418.

87. Zhou Y.X., Ting C.S., Liu C.W. Design of groundwater monitoring networks with case ■f study of Pingtung Plain, Taiwan. Taiwan, Wu-Nan book inc. pp. 117.