Variability of Intrinsic Groundwater Vulnerability to Pollution in River Valley due to Groundwater Depth and Recharge Changes
Abstract
:1. Introduction
2. Material and Methods
2.1. Study Area
2.2. Analysis of Changes in the Depth to GroundwaterTable
2.3. Assessment of Recharge Value
2.4. Assessment of Vulnerability
2.5. Characteristic of DRASTIC Parameters
2.5.1. Depth to Groundwater
2.5.2. Net Recharge
2.5.3. Aquifer Media
2.5.4. Soil Media
2.5.5. Topography
2.5.6. Impact of Vadose Zone
2.5.7. Hydraulic Conductivity
2.6. Parameter Sensitivity Analysis
2.7. Vulnerability Validation
3. Results and Discussion
4. Conclusions
- The groundwater vulnerability assessment in the Vistula Valley was performed for three variants: base variant A—average depth to groundwater from the period of 1999–2013, hydrogeological drought variant B (lowest long-term depth—2003), and flood risk variant C (greatest long-term depth—2011). In variant A, the moderately high and medium vulnerability classes occupy the greatest area. During hydrogeological drought, the largest area is occupied by the medium vulnerability class. During the potential flood risk in the river valleys, the greatest area is occupied by the moderately high vulnerability class. All these variants reveal a spatial variability in the distribution of the individual classes.
- Based on sensitivity analysis, depth to water table was the most effective parameter responsible for highest variation in the vulnerability index. The effective weight of parameter D is much higher than the theoretical value. In all variants, the calculated IPZ index indicates (after optimization of the parameter values) the greatest increase in the area of the moderately high and high classes, while the area of the low class decreased. The greatest changes in relations to the theoretical values occurred in variant B for the lowest groundwater level variant.
- The studies confirm that the assessment of groundwater vulnerability to contamination in areas of shallow depth to groundwater table, which are GDE, should be made for a number of variants, with a particular focus on the highest and lowest observed groundwater levels. Groundwater vulnerability changes are also associated with sensitivity analysis and calculated effective weights.
- This study suggests that DRASTIC is an effective tool for the groundwater vulnerability assessment studies and can be used for the prioritization of the susceptible areas to avoid contamination as well as a frequent and detailed monitoring of pollution potential in the high and medium vulnerability zones. There is also a need for proper planning in areas of shallow groundwater and considerable changes. Vulnerability maps for various variants of input data can be the basis for plans of land management and flood protection, as well as for prediction of climate change effects.
Author Contributions
Funding
Conflicts of Interest
References and Notes
- Aller, L.; Bennett, T.; Lehr, J.H.; Petty, R.J.; Hackett, G. DRASTIC: A Standardized System of Evaluating Ground Water Pollution Potential Using Hydrogeological Settings; EPA no. 600287035; USEPA: Washington, DC, USA, 1987.
- Secunda, S.; Collin, M.L.; Melloul, A.J. Groundwater vulnerability assessment using a composite model combining DRASTIC with extensive agricultural land use in Israel’s Sharon region. J. Environ. Manag. 1998, 54, 39–57. [Google Scholar] [CrossRef]
- Al-Adamat, R.; Foster, I.; Baban, S.M.J. Groundwater vulnerability mapping for the Basaltic aquifer of the Azraq basin of Jordan using GIS, Remote Sensing and DRASTIC. Appl. Geogr. 2003, 23, 303–324. [Google Scholar] [CrossRef]
- Krogulec, E. Vulnerability Assessment of Groundwater Pollution in the River Valley on the Basis of Hydrodynamic Evidence; Warsaw University Print: Warsaw, Poland, 2004. (In Polish) [Google Scholar]
- Krogulec, E. Intrinsic and specific vulnerability of groundwater in a river valley: Assessment, verification and analysis of uncertainty. J. Earth Sci. Clim. Chang. 2013, 4. [Google Scholar] [CrossRef]
- Babiker, I.S.; Mohamed, M.A.A.; Hiyama, T.; Kato, K. A GISbased DRASTIC model for assessing aquifer vulnerability in Kakamigahara Heights, Gifu Prefecture, central Japan. Sci. Total Environ. 2005, 345, 127–140. [Google Scholar] [CrossRef]
- Rahman, A. A GIS based DRASTIC model for assessing groundwater vulnerability in shallow aquifer in Aligarh, India. Appl. Geogr. 2008, 28, 32–53. [Google Scholar] [CrossRef]
- Al Hallaq, A.; Elaish, A. Assessment of aquifer vulnerability to contamination in Khan younis Governorate, Gaza Strip, Palestine, using the DRASTIC model within GIS environment. Arab. J. Geosci. 2011, 1, 1–15. [Google Scholar] [CrossRef]
- Panagopoulos, G.; Antonakos, A.; Lambrakis, N. Optimization of DRASTIC model for groundwater vulnerability assessment, by the use of simple statistical methods and GIS. Hydrogeol. J. 2006. published online. [Google Scholar] [CrossRef]
- Pathak, D.R.; Hiratsuka, A.; Awata, I.; Chen, L. Groundwater vulnerability assessment in shallow aquifer of Kathmandu Valley using GIS-based DRASTIC model. Environ. Geol. 2009, 57, 1569–1578. [Google Scholar] [CrossRef]
- Saidi, S.; Bouri, S.; Dhia, H.B. Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer, Tunisia: A validation study. Hydrol. Sci. J. 2011, 56, 288–304. [Google Scholar] [CrossRef]
- Usman, M.; Liedl, R.; Kavousi, A. Estimation of distributed seasonal net recharge by modern satellite data in irrigated agricultural regions of Pakistan. Environ. Earth Sci. 2015. [Google Scholar] [CrossRef]
- Lobo-Ferreira, J.P.; Delgado-Rodrigues, J. BALSEQ—A model for the estimation of water balances, including aquifer recharge, requiring scarce hydrogeological data. In Estimation of Natural Groundwater Recharge; Simmers, E., Ed.; NATO ASI Series Vol. 222; D. Reidel: Dordrecht, The Netherlands, 1988. [Google Scholar]
- Lobo-Ferreira, J.P.; Oliveira, M.M. DRASTIC ground water vulnerability mapping of Portugal. In Proceedings of the 27th Congress of the International Association for Hydraulic Research, San Francisco, CA, USA, 10–15 August 1997; pp. 132–137. Available online: http://www-dh.lnec.pt/gias/novidades/drastic_e.html (accessed on 28 December 2018).
- Ducci, D. The spatial variability of climatic data in the hydrogeologic budget. In Proceedings of the International Conference on “ModelCARE’99”, Zurich, Switzerland, 20–23 September 1999; pp. 673–678. [Google Scholar]
- Farjad, B.; Mohamed, T.A.; Wijesekara, N.; Pirasteh, S.; Shafri, H.Z.B.M. Groundwater intrinsic vulnerability and risk mapping. Proc. ICE Water Manag. 2012, 165, 441–450. [Google Scholar] [CrossRef]
- Krogulec, E. Evaluating the risk of groundwater drought in groundwater-dependent ecosystems in the central part of the Vistula River Valley, Poland. Ecohydrol. Hydrobiol. 2018, 18, 82–91. [Google Scholar] [CrossRef]
- Directive 2007/60/EC of the European Parliament and of the Council of 23 October 2007 on the assessment and management of flood risks.
- Regulation of the Council of Ministers of January 16, 1959 on the creation of the Kampinos National Park (Journal of Laws 1959, No. 17, item 91).
- UNESCO MaB; since the year 2000. Available online: http://www.unesco.org/new/en/natural-sciences/environment/ecological-sciences/man-and-biosphere-programme/ (accessed on 15 December 2018).
- Available online: http://natura2000.gdos.gov.pl/ (accessed on 15 December 2018).
- Krogulec, E.; Zabłocki, S.; Sawicka, K. Changes in groundwater regime during vegetation period in groundwater dependent ecosystems. Acta Geol. Pol. 2016, 66, 525–540. [Google Scholar] [CrossRef]
- Baraniecka, M.D.; Konecka-Betley, K. Fluvial Sediments of the Vistulian and Holocene in the Warsaw Basin. Evolution of the Vistula River Valley During the last 15 000 Years. Geogr. Stud. 1987, 4, 111–171. (In Polish) [Google Scholar]
- Sarnacka, Z. Stratigraphy of Quaternary sediments: Warsaw and surrounding areas. Pr. Panstw. Inst. Geol. 1992, 138, 1–36. (In Polish) [Google Scholar]
- Krogulec, E. Evaluation of hydrogeological conditions based on modelling: A case study in the Kampinos National Park region, central Poland. Geol. Q. 2003, 47, 83–90. [Google Scholar]
- Krogulec, E.; Marciniak, M. Laboratory and Field Studies of Conductivity in Kampinos National Park Area; Współczesne Problemy Hydrogeologii: GUT Publishing House, Gdansk, 2003; pp. 361–369. (In Polish) [Google Scholar]
- Krogulec, E.; Andrzejewska, A.; Furmankowska, A.; Zabłocki, S. Analysis of groundwater level in wetland areas on the basis of observations in Kampinos National Park. Biul. Panstw. Inst. Geol. 2009, 436, 281–288. (In Polish) [Google Scholar]
- Kleczkowski, A.S. (Ed.) Map of Major Groundwater Reservoirs (GZWP) in Poland Requiring Special Protection, 1:500 000; AGH: Krakov, Poland, 1990. (In Polish) [Google Scholar]
- Krogulec, E. Numerical Analysis of the Structure of the Filtration Flux in the Edge Zone of Blonski Level; Warsaw University Print: Warsaw, Poland, 1997. (In Polish) [Google Scholar]
- Krogulec, E.; Zabłocki, S. Relationship between the environmental and hydrogeological elements characterizing groundwater-dependent ecosystems in central Poland. Hydrogeol. J. 2015, T23, 1587–1602. [Google Scholar] [CrossRef]
- McDonald, M.D.; Harbaugh, A.W. A Modular Three-Dimensional Finite-Difference Flow Model; Book 6; Techniques of Water Resources Investigations of the U.S. Geological Survey: Denver, CO, USA, 1988.
- Gruszczyński, T.; Krogulec, E. Numerical Model of the Filtration Field in the Area of the Kampinos National Park. Result no. 9 of Project “Development of Methods for the Reconstruction of the Primary Water Conditions in Kampinos National Park for Recovering Natural Degradation and Improvement State of Biodiversity”; EEA Grants. PL0268; Warsaw University of Life Sciences Press: Warsaw, Poland, 2010. [Google Scholar]
- Gruszczyński, T.; Krogulec, E. Selection of the renaturalization variant of wetlands in Kampinos National Park based on hydrodynamic model tests. Biul. Inst. Geol. 2012, 451, 45–52. (In Polish) [Google Scholar]
- Rosen, L. A study of the DRASTIC methodology with emphasis on Swedish conditions. Ground Water, 32:278-285 Rupert MG (2001) Calibration of the DRASTIC ground water vulnerability mapping method. Ground Water 1993, 39, 625–630. [Google Scholar]
- Knox, R.C.; Sabatini, D.A.; Canter, L.W. Subsurface Transport and Fate Processes; Lewis Publishers: Boca Raton, FL, USA, 1993. [Google Scholar]
- Barrocu, G.; Biallo, G. Application of GIS for aquifer vulnerability evaluation. In W: Application of Geographical Information Systems in Hydrology and Water Resources Management; HydroGIS93; No. 211; IAHS Publications: Vienna, Austria, 1993; pp. 571–580. [Google Scholar]
- Vrba, J.; Zaporozec, A. (Eds.) Guidebook on Mapping Groundwater Vulnerability; IAH International Contribution to Hydrogeology: Balkema, The Netherlands, 1994; Volume 16. [Google Scholar]
- Kim, Y.J.; Hamm, S.Y. Assessment of the potential for groundwater contamination using the DRASTIC/EGIS technique, Cheongju area, South Korea. Hydrogeol. J. 1999, 7, 227–235. [Google Scholar] [CrossRef]
- Hamza, M.H.; Added, A.; Rodríguez, R.; Abdeljaoued, S.; Mammou, A.B. A GIS-based DRASTIC vulnerability and net recharge reassessment in an aquifer of a semi-arid region (Metline-Ras Jebel-Raf Raf aquifer, Northern Tunisia). J. Environ. Manag. 2007, 84, 12–19. [Google Scholar] [CrossRef]
- Leone, A.; Ripa, M.N.; Uricchio, V.; Deák, L.; Vargay, Z. Vulnerability and risk evaluation of agricultural nitrogen pollution for Hungary’s main aquifer using DRASTIC and GLEAMS models. J. Environ. Manag. 2009, 90, 2969–2978. [Google Scholar] [CrossRef]
- Protection Plan of Kampinos National Park; Arch. KPN: Izabelin, Poland, 1995. (In Polish)
- Przytuła, E.; Cabalska, J.; Modliński, P. Hydrogeological Map of Poland 1:50 000, Sheet Kampinos (0520); Polish Geological Institute: Warsaw, Poland, 2000. (In Polish)
- Słowański, W.; Piechulska-Słowańska, B.; Gogołek, W. Geological Map of Poland 1:200 000, Sheet Warsaw; Polish Geological Institute: Warsaw, Poland, 1994. (In Polish)
- Herbich, P.; Ćwiertniewska, Z.; Fert, M. Methodological Recommendations for the Development of Information Layers of the GIS Database of the Hydrogeological Map of Poland 1:50 000 First Aquifer—Vulnerability and Groundwater Quality; Polish Geological Institute: Gdansk, Poland, 2008. (In Polish)
- Napolitano, P.; Fabbri, A.G. Single-Parameter Sensitivity Analysis for Aquifer Vulnerability Assessment Using DRASTIC and SINTACS; Application of Geographic Information Systems in Hydrology and Water Resources Management (Proceedings of the Vienna Conference, April 1996), IAHS Publ. no.235; IAHS Publications: Vienna, Austria, 1996; pp. 559–566. [Google Scholar]
- Elçi, A. Calibration of groundwater vulnerability mapping using the generalized reduced gradient method. J. Contam. Hydrol. 2017, 207, 39–49. [Google Scholar] [CrossRef]
- Sadat-Noori, M.; Ebrahimi, K. Groundwater vulnerability assessment in agricultural areas using a modified DRASTIC model. Environ. Monit. Assess. 2016, 188, 19. [Google Scholar] [CrossRef]
- Krogulec, E.; Trzeciak, J. Drastic assessment of groundwater vulnerability to pollution in the Vistula floodplain in central Poland. Hydrol. Res. 2017, 48, 1088–1099. [Google Scholar] [CrossRef]
Groundwater Depth Statistics (m b.g.l.) | All Monitoring Points (n = 56) in: | Selected Piezometers in 1999–2013 | ||||
---|---|---|---|---|---|---|
1999–2013 | 2003 | 2011 | P11 | P12 | P14 | |
Average | 1.87 | 2.11 | 1.35 | 1.89 | 0.47 | 0.74 |
Median | 1.92 | 2.09 | 1.35 | 1.94 | 0.43 | 0.75 |
Maximum | 2.46 | 2.46 | 1.67 | 2.43 | 1.15 | 1.42 |
Minimum | 1.11 | 1.59 | 1.11 | 1.15 | −0.45 | −0.10 |
Standard deviation | 0.28 | 0.21 | 0.16 | 0.28 | 0.31 | 0.32 |
Maximum amplitude | 1.35 | 0.87 | 0.56 | 1.28 | 1.60 | 1.52 |
Annual amplitude | 0.68 | 0.87 | 0.56 | 0.73 | 0.88 | 0.78 |
Parameter | Symbol | Weight | Data Source |
---|---|---|---|
Depth to groundwater | D | 5 | Monitoring; numerical groundwater flow modeling |
Net recharge | R | 4 | Numerical groundwater flow modeling |
Aquifer media | A | 3 | Analysis of profiles from borehole database; hydrogeological maps [42]; field works; geophysical data |
Soil media | S | 2 | Soil map of KNP [41]; archival data; direct studies |
Topography | T | 1 | Numerical Model Terrain, 10 × 10 m resolution |
Impact of vadose zone | I | 5 | Analysis of profiles from borehole database; Subsurface sediment map 1:200,000 [43]; field works; geophysical data |
Hydraulic conductivity | C | 3 | Archival field tests in wells; PARAMEX and BAT tests in piezometers and shallow boreholes [26,29] |
Groundwater Depth (m b.g.l.) | Range of Parameter D | Area of Occurrence in Selected Years (km2) | ||
---|---|---|---|---|
1999–2013 | 2003 | 2011 | ||
<1 | 10 | 171.35 | 62.57 | 322.24 |
1–3 | 9 | 266.34 | 277.26 | 176.21 |
3–5 | 8 | 89.39 | 157.38 | 48.6 |
>5 | 7 | 47.26 | 77.14 | 27.29 |
Recharge Rate (mm/year) | Range of Parameter R | Area of Occurrence in Selected Years (km2) | ||
---|---|---|---|---|
1999–2013 | 2003 | 2011 | ||
<−50 | 1 | 24.84 | 10.99 | 120.09 |
−50–0 | 2 | 139.74 | 120.77 | 146.44 |
0–20 | 3 | 69.12 | 107.08 | 57.06 |
20–40 | 4 | 77.47 | 290.69 | 55.69 |
40–60 | 5 | 263.18 | 32.42 | 85.21 |
>60 | 6 | 12.39 | 109.86 |
Parameter | Characteristic | Range | Area of Occurrence (km2) |
---|---|---|---|
Aquifer media A | Surface water | 0 | 1.22 |
Fine sands, silty sands with till inlays | 2 | 23.74 | |
Fine sands with silty alluvial deposits | 3 | 148.70 | |
Sands with inlays of weathered tills | 4 | 36.10 | |
River and fluvioglacial sands | 6 | 213.39 | |
Gravels and eolian sands | 8 | 183.54 | |
Soil media S | Dusty, mineral-muck, silty alluvial soil | 5 | 166.07 |
Brown soil, black soil | 6 | 133.08 | |
Eolian-erosive soil | 7 | 32.24 | |
Silt‒gley, peat soil, muck soil | 8 | 125.21 | |
Ground‒gley, gley‒podsolic, podsolic soil | 9 | 117.52 | |
Topography Slope (%) T | >3.0 | 7.5 | 12.92 |
2.5–3.0 | 8 | 8.81 | |
2.0–2.5 | 8.5 | 17.91 | |
1.5–2.0 | 9 | 23.57 | |
1.0–1.5 | 9.5 | 63.78 | |
0.0–1.0 | 10 | 447.47 | |
Impact of vadose zone I | Surface water | 0 | 1.22 |
Laminated clays, tills | 2 | 23.74 | |
Aggregate mud, silty alluvial deposits | 3 | 148.70 | |
Weathered tills | 4 | 36.10 | |
River and fluvioglacial sands | 6 | 213.39 | |
Gravels, eolian sands | 8 | 183.54 | |
Hydraulic conductivity (m/day) C | <4 | 1 | 2.96 |
4–12 | 2 | 87.15 | |
13–28 | 4 | 265.77 | |
29–40 | 6 | 148.65 | |
41–80 | 8 | 69.67 |
IPZ Index | Vulnerability Class | Vulnerability for Groundwater State Defined in Years: | ||
---|---|---|---|---|
1999–2013 | 2003 | 2011 | ||
Area (km2) | ||||
75–100 | Very low | 1.20 | 1.85 | 1.11 |
101–125 | Low | 138.80 | 160.27 | 141.11 |
126–150 | Medium | 193.21 | 216.26 | 197.99 |
151–175 | Medium high | 231.87 | 190.75 | 221.75 |
176–200 | High | 5.68 | 2.12 | 9.29 |
Mean IPZ | 142.7 | 139.5 | 142.99 | |
Standard deviation | 19.33 | 18.77 | 20.54 | |
IPZ scope | 80.5–181 | 76.5–180 | 80.5–185 |
Parameter | “Theoretical” Weight | “Theoretical” Weight (%) | “Effective” Weight | “Effective” Weight (%) | |||
---|---|---|---|---|---|---|---|
Mean | Minimum | Maximum | Standard Deviation | ||||
D | 5 | 21.74 | 7.4 | 32.28 | 20.47 | 59.52 | 6.63 |
R | 4 | 17.39 | 2.4 | 10.37 | 2.44 | 23.81 | 3.54 |
A | 3 | 13.04 | 2.7 | 11.58 | 0 | 17.33 | 3.13 |
S | 2 | 8.70 | 2.2 | 9.59 | 5.78 | 19.88 | 1.89 |
T | 1 | 4.35 | 1.6 | 6.99 | 4.3 | 11.9 | 1.15 |
I | 5 | 21.74 | 4.4 | 19.3 | 0 | 28.88 | 5.21 |
C | 3 | 13.04 | 2.3 | 9.89 | 1.94 | 23.76 | 3.71 |
Total | 23 | 100 | 23 | 100 |
IPZ Index | Vulnerability Class | Class Area after Optimization (km2) | Class Area Changes according to Original IPZ (%) | ||||
---|---|---|---|---|---|---|---|
Average Groundwater State 1999–2013 | Groundwater State in 2003 | Groundwater State in 2011 | Average State 1999–2013 | Groundwater State in 2003 | Groundwater State in 2011 | ||
<100 | Very low | 0.21 | 0.27 | 0.15 | −0.2 | −0.3 | −0.2 |
100–125 | Low | 5.38 | 18.38 | 2.73 | −23.4 | −24.8 | −24.2 |
125–150 | Medium | 175.67 | 186.68 | 173.38 | −3.1 | −5.2 | −4.3 |
151–175 | Medium‒high | 314.32 | 328.80 | 287.30 | 14.4 | 24.2 | 11.5 |
>175 | High | 75.67 | 37.12 | 107.69 | 12.3 | 6.1 | 17.2 |
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Krogulec, E.; Zabłocki, S.; Zadrożna, D. Variability of Intrinsic Groundwater Vulnerability to Pollution in River Valley due to Groundwater Depth and Recharge Changes. Appl. Sci. 2019, 9, 1133. https://doi.org/10.3390/app9061133
Krogulec E, Zabłocki S, Zadrożna D. Variability of Intrinsic Groundwater Vulnerability to Pollution in River Valley due to Groundwater Depth and Recharge Changes. Applied Sciences. 2019; 9(6):1133. https://doi.org/10.3390/app9061133
Chicago/Turabian StyleKrogulec, Ewa, Sebastian Zabłocki, and Danuta Zadrożna. 2019. "Variability of Intrinsic Groundwater Vulnerability to Pollution in River Valley due to Groundwater Depth and Recharge Changes" Applied Sciences 9, no. 6: 1133. https://doi.org/10.3390/app9061133
APA StyleKrogulec, E., Zabłocki, S., & Zadrożna, D. (2019). Variability of Intrinsic Groundwater Vulnerability to Pollution in River Valley due to Groundwater Depth and Recharge Changes. Applied Sciences, 9(6), 1133. https://doi.org/10.3390/app9061133