Next Article in Journal
Behavior of Circular CFST Columns Subjected to Different Lateral Impact Energy
Next Article in Special Issue
Treatment of Organic Matter and Tetracycline in Water by Using Constructed Wetlands and Photocatalysis
Previous Article in Journal
Research and Development of Automatic Monitoring System for Livestock Farms
Previous Article in Special Issue
Role of Wetland Plants and Use of Ornamental Flowering Plants in Constructed Wetlands for Wastewater Treatment: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variability of Intrinsic Groundwater Vulnerability to Pollution in River Valley due to Groundwater Depth and Recharge Changes

Faculty of Geology, University of Warsaw, Żwirki i Wigury 93, 02-089 Warszawa, Poland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(6), 1133; https://doi.org/10.3390/app9061133
Submission received: 10 February 2019 / Revised: 7 March 2019 / Accepted: 11 March 2019 / Published: 18 March 2019

Abstract

:
Groundwater management can be effectively implemented by assessment of groundwater vulnerability to pollution. The research was carried out in the Vistula River valley (Poland) in an area of shallow groundwater occurrence, defined as a Groundwater-Dependent Ecosystem area. The goal of this study was to identify the average, maximum, and minimum depths of the groundwater table for variant assessment of groundwater vulnerability to contamination. The variants correspond to the average vulnerability, the vulnerability during hydrogeological drought, and the vulnerability during the flood risk period. Theoretical and effective weights of vulnerability parameters were calculated using the DRASTIC method. Vulnerability maps constructed for the various vulnerability variants and by using different parameter weights show the spatial variability of the individual vulnerability classes. Due to the specifics of this protected area, the expected dependency between vulnerability index and chloride concentrations in the monitoring points was not found. Comparison showed strong dependency of water chemistry and the value of recharge, lithology of aquifer, and unsaturated zone. The research results confirm the need for variant vulnerability assessment to protect against floods or predict the effects of climate change.

1. Introduction

Maps of groundwater vulnerability to contamination, which are the main tool in spatial planning and prediction of the effects of aquifer pollution, are compiled using many methods and procedures. Parameters and data taken for the vulnerability assessment are represented by average values and therefore are characterized by an “average” groundwater vulnerability to contamination. In the assessment of natural groundwater vulnerability to contamination, DRASTIC is the rating method commonly used. It is dedicated to areas of porous aquifers [1]. In the DRASTIC method, like in other methods of groundwater vulnerability assessment, the depth to groundwater table (D) and the groundwater recharge by infiltration (R) depend on climate, weather, hydrological, and environmental changes and are the main parameters of vulnerability assessment, especially in shallow groundwater areas. The value of D is taken for the vulnerability assessment based on statistical data from long-term monitoring [2,3,4,5,6,7,8], and the value of R is based on averaged precipitation data, lithologies of the near-surface zone, or retention data [4,6,9,10,11,12]. The pollution vulnerability can be defined as ‘the sensitivity of groundwater quality to an imposed contaminant load’ depending only on the ‘proper characteristic of the aquifer’. The vulnerability is relatively static, excluding some variation over time due, e.g., to the piezometric level changes [13,14,15,16].
The aim of the studies was the identification of characteristic groundwater levels in the Vistula River valley for the variant vulnerability assessment. During the studies, the variants of groundwater vulnerability to contamination were calculated by taking into account the average, maximum, and minimum values of D. These variants correspond to the average vulnerability, vulnerability in conditions of hydrogeological drought [17], and vulnerability in the conditions of floods in river valleys [18]. The groundwater vulnerability assessment was carried out for three variants of infiltration: for the average groundwater levels from the period 1999–2013, and for the lowest and highest (average annual) values of groundwater levels. Sensitivity analysis and map removal parameter for all the variants, interpretation of vulnerability variations due to groundwater and recharge changes, and vulnerability validation with hydrochemical data usage were carried out.

2. Material and Methods

2.1. Study Area

The analysis of vulnerability assessment changes was performed on a well-examined area located in the Middle Vistula Valley, Poland. It has been covered by legal protection since 1959 as the Kampinos National Park [19], Biosphere Reservoir [20] and Natura 2000—Kampinos Forest [21]. The national park covers an area of 385 km2. The vulnerability assessment was carried out in the part of the Vistula Valley that covers 571 km2. Valuable plant communities whose existence is ensured by the shallow groundwater table cause areas to be classified as groundwater-dependent ecosystems (GDE), which cover approximately 62% of the park acreage [22]. The main factor determining proper functioning of these GDEs is dependency: precipitation-infiltration-ecosystem type-surface water drainage. In the GDE areas, about 33% is covered by typical swamp vegetation, 6% is represented by wet communities, and 19% are dry communities, where the humidity of the aeration zone and the depth of the groundwater table are insufficient to maintain the hydrophilic vegetation [22].
In the Vistula Valley, Quaternary sediments are the collector of groundwater. The surface of tertiary loams, constituting a floor of a quaternary water-bearing layer, in the Vistula valley in the region of KNP occurs from 2 to 54 m a.s.l. The top part of the water-bearing layer has a sandy and sandy‒gravel character; the bottom is created by sandy-silt sediments, in places changing to sandy clay and dusty clay [23,24,25]. The hydrogeological conditions of KNP have been the subject of several research projects, summarized in several publications [4,5,26,27]. Since its resources are significant, this is the main usable aquifer and the major groundwater basin No. 222 [28]. The regional groundwater flow occurs to the north, towards the Vistula River. In some regions, there is also a local groundwater circulation, where recharge areas are represented by dunes of higher elevations, whereas swamp areas with their minor watercourses and canals are the drainage base [29]. The study area comprises characteristic longitudinal zones of similar geological structure and hydrogeological conditions [4]: the Vistula flood plain terrace, northern and southern marsh zones, northern and southern dune zones, and the accumulative-erosional Warsaw‒Błonie terrace (Figure 1).

2.2. Analysis of Changes in the Depth to GroundwaterTable

The Vistula River Valley is a typical shallow groundwater area. The groundwater levels are under regular monitoring within a stationary monitoring network. Groundwater level measurements have been conducted since 1999 on a fortnight basis in 56 piezometers in KPN. An analysis of changes of the groundwater table level has provided data on the quantity of groundwater resources. They include the determination of: changes over a long period [4,25], seasonal changes of resources [30], changes in the amount of resources during the vegetation period and the periods and sites of possible moisture deficit in the GDE areas [22], and duration of and depth to groundwater table during hydrogeological drought [17].
Statistical analysis of groundwater levels in Vistula valley included basic statistics for the 1999–2013 period and particular years of period were analyzed and the scope of groundwater levels was defined. The statistics allow for determining the range of long-term changes in the depth to groundwater table and seasonal changes, as well as defining the highest and lowest levels in a long period. The average depth to groundwater table in the piezometers of the monitoring network during the period of 1999–2013 was 1.87 m b.g.l., and the range of changes varied from 1.11 to 2.46 m. The individual piezometers record typical seasonal level fluctuations. The years 2003 and 2011 stand out for a distinctly different trend of seasonal variations in the groundwater table depth. This phenomenon was observed in all zones of the KNP—marsh, dune, and floodplain. In 2003, the average depth was 2.11 m b.g.l., and the amplitude was 0.87 m. The highest groundwater levels were recorded in 2011, when the average depth to groundwater table was 1.35 m b.g.l., and the amplitude was 0.56 m (Table 1).
In the period of 1999–2013, in marsh zones, the groundwater table occurred at average depths ranging 0–1 m in the central parts (25% of the total area) to 1–3 m in the peripheries (46%) (Figure 2). The zones where the groundwater table was above the ground surface comprised ca. 6.5% of the study area. The greatest depth to groundwater table (>5 m) was found in dune zones and at the boundary with the upland area, attaining a maximum value of 13.5 m.
There are distinct differences between the groundwater level for the long period and the lowest (2003) and highest (2011) levels in the period of 1999–2013. In 2003, the groundwater table in the marsh zones was at depths of 0–1 m (12%) or 1–2 m (49% of the total area), which indicates drying of organic soils whose maximum thickness does not exceed 1.5 m [27]. We have also observed that the areas where the groundwater table occurs at depths of >5 m doubled, and the maximum groundwater table depth was determined to be 14.5 m.
In 2011, the groundwater table was ca. 0.54 m above the ground surface on average, and locally up to 2 m. This phenomenon was observed in almost the entire area of the marsh zones, i.e., ca. 29% of the total study area. The depth intervals of 0–1 and 1–3 m account for 29% of the area each. The calculated maximum depth was 11.89 m.

2.3. Assessment of Recharge Value

Hydrodynamic modeling studies in the Vistula valley were performed to provide a quantitative description of the components of groundwater balance within the Quaternary aquifer and create spatial distribution of groundwater depth. The modelling was carried out using the Visual ModFLOW 4.2. software [31]. The model was built under the steady state, defined as the average annual hydrodynamic state in the period of 1999–2008 and it concerned variants related to renaturalization of wetland areas [32,33]. Lateral boundaries of this hydrogeological system were delimited along the main rivers of the region: in the north the Vistula River, in the west the Bzura River, in the south along the edge of the accumulative-erosional level, and in the east along the border of Warsaw. The lower boundary was set as the top part of poorly permeable Neogene clays.
The relationship between the aquifer system and its environment was projected as boundary conditions: RIVER type—for Bzura, Vistula and small canals; General Head type—along south boundary with groundwater flow into Vistula valley from Warsaw‒Blonie terrace side. The discretization step was set to 100 × 100 m. The calibration process consisted of applying various values of infiltration recharge and underground evaporation with simultaneous verification of flow rates in rivers and canals [32]. Then, the hydrodynamic model was updated by applying three calculation variants of infiltration: modeling for average groundwater levels from the period of 1999–2013, for the lowest level (in 2003), and for the highest level (in 2011) [30]. The correlation statistics obtained at 49 monitoring points (located in the model area) and calculated head were as follows for all variants: standard error of the estimate 0.023–0.044 m; root mean squared 0.158–0.308 m; normalized RMS 0.953–1.849%.
The rate of infiltration recharge was calculated from the hydrodynamic model for the period of 1999–2013 (average value) and for the lowest and highest groundwater levels (Figure 3). The greatest recharge rate values are found predominantly in dune areas. In marsh zones, which are characterized by the highest evaporimetric potential (hydrogenic soils), due to the shallow depth to groundwater, intensification of evapotranspiration-related losses causes a so-called negative recharge in the annual cycle. An analysis of seasonal changes [22] has shown that the effective recharge in the marsh areas usually occurs in springtime due to thawing. This is why pollutants can migrate with the percolating water, although, in the annual cycle, the observed losses are higher than the recharge rate. This phenomenon was the reason for covering the areas of negative annual balance by the vulnerability analysis. The recharge rate varies spatially within the range of −89–65 mm/year, and in the individual model variants, the differences in the average values are 6–12 mm/year. For the average groundwater level, the average recharge rate in this area was approximately 21 mm/year, with standard deviation of 35 mm/year. In 2003, the evaporation rate from the groundwater table in the marsh areas was higher than the recharge rate, and remained at a similar level to the average (−50–0 mm/year). In the dune areas, the recharge rate was clearly lower, by ca. 20 mm/year. The average recharge rate for this period was 15 mm/year, with an excursion up to 25 mm/year. In 2011, the recharge rate was 3 mm/year. In the dune areas, high recharge rates dominated (>60 mm/year), whereas in the marsh areas, low values prevailed (<−50 mm/year), which was related to intense evaporation from the groundwater table occurring at the shallowest depths as compared to the other variants.

2.4. Assessment of Vulnerability

The parametric method DRASTIC (indexation of parameters) was used to assess the natural vulnerability to contamination, dedicated to porous water-bearing media. The most important assumptions of this method include the advective migration of conservative pollutants from the ground surface through the vadose zone. DRASTIC is a method often used to assess groundwater vulnerability to a wide range of potential contaminants. It is a standard tool for assessing groundwater vulnerability to contamination, and is applied in many countries [2,10,11,34,35,36,37,38,39,40].
Calculations of the DRASTIC index (IPZ) are conducted according to the following formula [1]:
IPZ = DR Dw + RRRw + ARAw + SRSw + TRTw + IRIw + CRCw,
where R—parameter rating (Table 2); W—parameter weight (Table 2); D—depth to groundwater; R—net recharge; A—aquifer media; S—soil media; T—topography; I—impact of vadose zone; C—hydraulic conductivity of the aquifer.
An advantage of this method is that the used parameters comprehensively characterize the analyzed multi-aquifer system by providing information about its both typical hydrogeological elements (groundwater table depth, recharge rate, aquifer lithology, hydraulic conductivity) and environmental elements that include terrain topography and the type of soil profile and its parent lithology of the vadose zone.
Because the circulation conditions are very well known in the study area, the analysis was based on various, highly detailed data (Table 2). They comprised the analysis of sections of research and hydrogeological boreholes, geophysical surveys, maps, including those compiled within the framework of the Protection plan of Kampinos National Park [41].
To calculate the index of vulnerability (IPZ index), a high-resolution (10 × 10 m) digital terrain model was also applied. A separate issue was the analysis of spatial variability of hydraulic conductivity, whose values were determined based on grain-size analyses, pumping tests, and direct studies using the PARAMEX [26] and BAT tests [29]. For the calculations of groundwater vulnerability to contamination, the original weight values were taken ([1]; Table 2).
The assessment of groundwater vulnerability to contamination was performed for three variants of the depth to groundwater table and infiltration recharge: variant A—average groundwater levels from the period of 1999–2013, treated as long-term average values; variant B—the lowest long-term levels from 2003; and variant C—the highest long-term levels from 2011. The vulnerability calculations were performed using GIS software. Continuous spatial data (i.e., soil types, aquifer media etc.) were discretized using ArcGIS 10.3 and the Raster calculator tool.

2.5. Characteristic of DRASTIC Parameters

2.5.1. Depth to Groundwater

The depth to groundwater table was determined for the period of 1999–2013 and for the lowest (in 2003) and highest (in 2011) recorded groundwater levels.
In 1999–2013, the largest area was that where the groundwater table was present at 1–3 m depth (more than 260 km2) and where it occurred at shallow depths down to 1 m b.g.l. (more than 170 km2). In 2003, the groundwater table was found at greater depths; the area where the groundwater table occurred at depths of 1–5 m was 430 km2, and the area where it remained below a depth of 5 m was greater than 75 km2. In 2011, the depth to groundwater table occurred most frequently at depths down to 1 m (more than 320 km2). Relevant parameter rating values were attributed to the individual depth intervals (Table 3).

2.5.2. Net Recharge

The division of the R parameter into 6 ranks was made on the basis of the recharge variability in particular periods of analysis (Table 4). In 1999–2013, the largest area was that where the recharge was in the range of 40–60 mm/year (over 263 km2), and mostly represented dune area. Second area with negative recharge from −50 to 0 mm/year (over 139 km2) occurred in the marsh zones. In 2003, the recharge took lower values with 290 km2 area represented by range of 20–40 mm/year. Respectively lower values of the negative recharge were a result of a deeper groundwater occurring. In 2011, the recharge was again bimodal with a clearer selection of the highest values, i.e., exceeding 60 mm/year (109 km2) and lowest values, i.e., below −50 mm/year (120 km2).

2.5.3. Aquifer Media

The aquifer lithology is generally inhomogeneous. Slight variations are observed as regards thin till interlayers, weathered tills or alluvial muds in sands. The aquifer lithology is represented by three main sediment types (Table 5): fluvial and glaciofluvial sands (37% of the area), gravels and eolian sands (32%), and fine-grained sands with alluvial muds (26%). Their distribution is related to the extents of the individual generations of the Vistula River terraces, as well as to the near-surface deposits at the shallow groundwater table (Figure 4).

2.5.4. Soil Media

The marsh zones are characterized by the greatest heterogeneity with respect to soil types. In this area, there are peaty soils, black earth soils, as well as alluvial soils composed of the finest grain fractions. The dune areas are covered mainly with podzols and gleyed podzols with poorly developed soil horizons, especially the humic layer. The soil ranking scheme (Table 5) was based on defining their protective ability against the migration of pollutants into groundwater, resulting from their field water capacity [44].

2.5.5. Topography

The area is characterized by low relief gradients, ranging predominantly up to 2.5%. Gradients higher than 3% were found on dune slopes; they account for only 2.2% of the area. The most common areas are those favoring infiltration, where the land gradient does not exceed 1% (78% of the area), with the T parameter rank of 10 (Table 5).

2.5.6. Impact of Vadose Zone

The I parameter has similar classes as the aquifer lithology (Table 5). The shallow groundwater table occurs within the same lithology, causing the parameters A and I to be very similar. Their high mutual correlation can be the basis for concluding on the necessity of modifying the DRASTIC assessment due to repeated information in the input layers (Table 2).

2.5.7. Hydraulic Conductivity

The hydraulic conductivity values range from 13 to 28 m/day in most of the area (more than 265 km2), with an average of ca. 25 m/day in 46% of the area. Slightly higher values, 29–40 m/day, are reported for the Vistula floodplain (26%). Zones of increased water permeability occur in the western, central and north-eastern parts of the area, covering ca. 70 km2 (Table 5). Low values, below 4 m/day, are found locally in the southeastern part of the area, occupying 2.96 km2.

2.6. Parameter Sensitivity Analysis

To determine the effect of weights of the individual parameters applied in the DRASTIC method, a single parameter sensitivity analysis was performed [3,6,45]. This analysis determines, by comparing the theoretical and effective weights, the possibility of optimizing the results depended on the input data and the relations between them. The “effective” weight of a parameter is obtained with regard to the other parameters of DRASTIC assessment. To assess the “effective” weight of each parameter—W was calculated using the following formula:
W = n = 1 7 ( parameter   weight   ×   parameter   rating / DRASTIC   index ) × 100 ,
where W is the effective parameter weight; the other symbols are as in Equation (1).

2.7. Vulnerability Validation

Validation of vulnerability is carried out through correlation with chemical substances, which can be commonly found in groundwater, and could be treated as an indicator of groundwater pollution [11,46,47]. In order to perform the validation of variant assessments of intrinsic vulnerability, a correlation analysis was carried out. Correlation analysis consisted in searching for a simple (linear) dependence between datasets and for dependence determined as a result of graphical interpretation of a point cloud. As an indicator of validation, the averaged chloride concentrations from 44 piezometers were used. The chlorides were determined every quarter in the period of 2006–2009. The validation of vulnerability variant of the average hydrodynamic state for theoretical and effective weights (the period 1999–2013) was made.

3. Results and Discussion

The DRASTIC index was calculated for three variants of parameter D: variant A—for the average depth to groundwater table from 1999 to 2013; variant B—for the lowest groundwater levels in a long period in 2003; and variant C—for the highest groundwater levels in a long period in 2011.
The groundwater vulnerability to contamination, determined for the average level from the period 1999–2013, is characterized by a high proportion of the area with a moderately high vulnerability class (IPZ = 151–175), accounting for 41% of the study area, i.e., 231.87 km2. The moderately high vulnerability occurs mainly in the dune zones (Figure 5). Despite the shallow depth to groundwater table in the marsh zones, the vulnerability is lower by 1–2 points, ranging from medium (IPZ = 126–150) in the margins of the zone, to low (IPZ = 101–125) in its central parts. The low vulnerability class is determined for the Vistula floodplain. The total area classified as having low vulnerability covers 138.80 km2. The high vulnerability class (IPZ > 176) is identified in the area of 5.68 km2 in some patches of the southern zone and in the suburbs of Warsaw, in the eastern part of the study area.
The assessment of groundwater vulnerability to contamination for the average depth and recharge rates in the years 1999–2013 was the basis for the analysis of vulnerability changes caused by the fluctuations in the depth to groundwater table and recharge rates. The vulnerability determined for average groundwater levels from a long-term period has been defined as the base state. Vulnerability maps, constructed based on average input data values in various vulnerability assessment methods and procedures, are the basis for making spatial development and groundwater protection decisions. In the case of hydrogeological systems with unconfined groundwater table and significant changes in its dynamics, it is necessary to analyze vulnerability changes for various variants.
The vulnerability determined for the lowest average levels has been defined as the vulnerability during hydrogeological drought. In the conditions of hydrogeological drought, the vulnerability is distinctly lower throughout the area. The medium vulnerability class occupies the greatest part of the study area, i.e., 216.26 km2 (338%), mainly in the dune zones (Figure 5). The moderately high vulnerability covers an area of 190.75 km2. Despite the shallow depth to groundwater table in the marsh zones, the vulnerability score is lower by 1–2, ranging from medium in the margins of the zone to low in its central parts. The low vulnerability class covers the Vistula River floodplain. The total area classified as having low vulnerability is 138.80 km2. The high vulnerability class is determined for the total area of 2.12 km2 in a few patches of the southern dune zone and in the suburbs of Warsaw, in the eastern part of the study area.
The assessment results were the basis for analyzing the differences in the areas occupied by the individual vulnerability classes determined for the base state and for the conditions of hydrogeological drought. In general, the class distribution pattern has not changed significantly. The most distinct changes refer to the size of highly vulnerable areas (a decrease from 0.99% to 0.37%) and to the increase in the very low class area (from 0.21% to 0.32%). The analysis of spatial variations in the IPZ index values shows that its greatest decrease in relation to the average value was 15, and the average decrease was 3.19. Considering this change with regard to the size of the classes (IPZ = 25), the changes can be deemed insignificant.
In variant C, for the highest groundwater levels in 2011, the vulnerability is close to the base level. The greatest area is covered by the moderately high vulnerability class, which is 221.75 km2 of the study area, predominantly in the dune zones (Figure 5). Highly vulnerable areas occupy 197.99 km2. The low vulnerability class occurs chiefly in the Vistula floodplain, and the total area classified as having low vulnerability is 141.11 km2. The high vulnerability class occupies the total area of 9.29 km2, in a few patches of the southern dune zone and in the suburbs of Warsaw—in the eastern part of the study area. When analyzing the differences between the base vulnerability and the flood risk vulnerability (Table 6), a slight increase in the area of low vulnerability class and nearly a twice increase in the area of high vulnerability class (from 0.99% to 1.62%; Figure 5) can be noticed. The analysis of spatial differences in the IPZ values has shown an average increase by merely 0.3, with a maximum increase by 14.5, although there are areas where the IPZ decreased (by up to 13).
A single parameter sensitivity analysis was carried out and, after determination of new weights, the IPZ index was recalculated for the above-mentioned three variants dependent on the groundwater levels and recharge rates (Figure 2 and Figure 3). The greatest increase in the effective weight in comparison with the theoretical weight refers to parameter D; the effective value exceeds the theoretical one by more than 50% (Table 7; Figure 6). Similar values of effective weights were also obtained in the sensitivity analysis concerning Vistula flood plain itself [48], which means that this dependence is characteristic for the whole area covered by the analysis. For parameter I, which has the theoretical weight of 5, like parameter D, the calculated effective weight is lower and determined at 4.4. For parameter R, which is characterized by the theoretical weight of 4, the effective weight is much lower and it was calculated on 2.4. Lower effective weights have also been obtained for parameter A, from 3 to 2.7, and for parameter C, from 3 to 2.3 (Table 7). The calculations made by the single parameter sensitivity analysis confirm the role of the depth to groundwater table in the assessment of groundwater vulnerability to contamination using the DRASTIC method.
In relation to the vulnerability classes determined earlier, there was an increase in the IPZ index in all cases (Figure 7; Table 8). The percentage of the area of the lowest vulnerability classes, in particular the low class, decreased by ca. 23–24%. The area occupied by the moderately high and very high classes increased by several percent. An increase in the very high class area is lower (by 6.1%) only for the low groundwater level variant (in 2003). Moderately high and very high classes dominate in the dune areas, whereas the axial parts of the marsh areas are covered by the medium class (locally passing into the low class), with the moderately high class at the margins. The largest area of the low and very low classes occurs in the south, along the margin of the Błonie Level. The change in the DRASTIC index value for the average level during the period was from −1 to 26, at the average value indicating an increase, after optimization, by 14.4. In 2003, the index change varied from 2 to 27, with an average of 13.9, while in 2011, the range of changes was 0–27, with an average of 16. The greatest increase in the index value was found for the marsh areas, which means that the contribution of parameter D is underestimated when the theoretical weights remain at the same level, and the obtained vulnerability classes are lower than those inferred from the actual effect of the depth on the ultimate vulnerability assessment. No changes are observed in the dune areas, where the vulnerability change is low (Figure 7).
The study area is characterized by low concentrations of chlorides, which reach values from 4.2 to 111 mg/L, with an average of 28.2 mg/L (median: 17.6 mg/L). The lowest concentrations occur on dune areas, while the highest are in drainage zones near the canals. Validation of vulnerability as an expected high dependence of chloride concentration and IPZ index was found as the reverse. For low values of Cl ions, the highest IPZ indexes were calculated (Figure 8). This results from the lack of source of groundwater contamination and intensity of infiltration. Intense recharge in the area of dune belts causes a slight transformation of precipitation water. In drainage zones, due to the long residence time of water in the system, processes such as evaporation result in much higher chloride concentrations.
The results of the correlation analysis indicate a better adjustment of the trend line to the results of the DRASTIC method with the original weights, where lesser importance is attributed to depth, and higher to the recharge and lithology of the aquifer and the vadose zone.

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

Conceptualization, E.K.; methodology, E.K. and S.Z.; software, S.Z. and D.Z.; validation, S.Z.; formal analysis, S.Z.; investigation, E.K., S.Z. and D.Z.; resources, E.K.; data curation, E.K.; writing—original draft preparation, E.K.; writing—review and editing, E.K., S.Z. and D.Z.; visualization, D.Z.; supervision, E.K.; project administration, E.K.; funding acquisition, E.K.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References and Notes

  1. 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.
  2. 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]
  3. 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]
  4. 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]
  5. 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]
  6. 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]
  7. 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]
  8. 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]
  9. 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]
  10. 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]
  11. 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]
  12. 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]
  13. 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]
  14. 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).
  15. 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]
  16. 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]
  17. 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]
  18. Directive 2007/60/EC of the European Parliament and of the Council of 23 October 2007 on the assessment and management of flood risks.
  19. 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).
  20. 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).
  21. Available online: http://natura2000.gdos.gov.pl/ (accessed on 15 December 2018).
  22. 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]
  23. 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]
  24. Sarnacka, Z. Stratigraphy of Quaternary sediments: Warsaw and surrounding areas. Pr. Panstw. Inst. Geol. 1992, 138, 1–36. (In Polish) [Google Scholar]
  25. 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]
  26. 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]
  27. 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]
  28. 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]
  29. 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]
  30. 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]
  31. 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.
  32. 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]
  33. 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]
  34. 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]
  35. Knox, R.C.; Sabatini, D.A.; Canter, L.W. Subsurface Transport and Fate Processes; Lewis Publishers: Boca Raton, FL, USA, 1993. [Google Scholar]
  36. 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]
  37. Vrba, J.; Zaporozec, A. (Eds.) Guidebook on Mapping Groundwater Vulnerability; IAH International Contribution to Hydrogeology: Balkema, The Netherlands, 1994; Volume 16. [Google Scholar]
  38. 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]
  39. 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]
  40. 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]
  41. Protection Plan of Kampinos National Park; Arch. KPN: Izabelin, Poland, 1995. (In Polish)
  42. 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)
  43. 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)
  44. 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)
  45. 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]
  46. Elçi, A. Calibration of groundwater vulnerability mapping using the generalized reduced gradient method. J. Contam. Hydrol. 2017, 207, 39–49. [Google Scholar] [CrossRef]
  47. 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]
  48. 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]
Figure 1. 3D view of study area location (vertical exaggeration 5.0; view from west; gridlines every 2 km).
Figure 1. 3D view of study area location (vertical exaggeration 5.0; view from west; gridlines every 2 km).
Applsci 09 01133 g001
Figure 2. Groundwater depth (kriging): (a) lowest (2003); (b) highest (2011); (c) average (1999–2013).
Figure 2. Groundwater depth (kriging): (a) lowest (2003); (b) highest (2011); (c) average (1999–2013).
Applsci 09 01133 g002
Figure 3. Annual recharge rate taken from numerical groundwater flow modeling: (a) 2003; (b) 2011; (c) average in 1999–2013.
Figure 3. Annual recharge rate taken from numerical groundwater flow modeling: (a) 2003; (b) 2011; (c) average in 1999–2013.
Applsci 09 01133 g003
Figure 4. Parameter ranges of: (a) aquifer media (A); (b) soil media (S); (c) topography (T); (d) impact of vadose zone (I); (e) hydraulic conductivity (C).
Figure 4. Parameter ranges of: (a) aquifer media (A); (b) soil media (S); (c) topography (T); (d) impact of vadose zone (I); (e) hydraulic conductivity (C).
Applsci 09 01133 g004
Figure 5. IPZ index for groundwater state defined in years: (a) 2003; (b) 2011; (c) 1999–2013.
Figure 5. IPZ index for groundwater state defined in years: (a) 2003; (b) 2011; (c) 1999–2013.
Applsci 09 01133 g005
Figure 6. Comparison of original and effective weights of parameters.
Figure 6. Comparison of original and effective weights of parameters.
Applsci 09 01133 g006
Figure 7. IPZ index after optimization by single parameter sensitivity analysis for groundwater state defined in years: (a) 2003; (b) 2011; (c) 1999–2013.
Figure 7. IPZ index after optimization by single parameter sensitivity analysis for groundwater state defined in years: (a) 2003; (b) 2011; (c) 1999–2013.
Applsci 09 01133 g007
Figure 8. Correlation of chloride concentration and IPZ index in selected piezometers.
Figure 8. Correlation of chloride concentration and IPZ index in selected piezometers.
Applsci 09 01133 g008
Table 1. Main statistics of groundwater depth based on monitoring observations in all the piezometers and in selected piezometers of dune (P11), marsh (P12), and Vistula flood plain (P14).
Table 1. Main statistics of groundwater depth based on monitoring observations in all the piezometers and in selected piezometers of dune (P11), marsh (P12), and Vistula flood plain (P14).
Groundwater Depth Statistics (m b.g.l.)All Monitoring Points (n = 56) in:Selected Piezometers in 1999–2013
1999–201320032011P11P12P14
Average1.872.111.351.890.470.74
Median1.922.091.351.940.430.75
Maximum2.462.461.672.431.151.42
Minimum1.111.591.111.15−0.45−0.10
Standard deviation0.280.210.160.280.310.32
Maximum amplitude1.350.870.561.281.601.52
Annual amplitude0.680.870.560.730.880.78
Table 2. DRASTIC parameters and sources of data used.
Table 2. DRASTIC parameters and sources of data used.
ParameterSymbolWeightData Source
Depth to groundwaterD5Monitoring; numerical groundwater flow modeling
Net rechargeR4Numerical groundwater flow modeling
Aquifer mediaA3Analysis of profiles from borehole database; hydrogeological maps [42]; field works; geophysical data
Soil mediaS2Soil map of KNP [41]; archival data; direct studies
TopographyT1Numerical Model Terrain, 10 × 10 m resolution
Impact of vadose zoneI5Analysis of profiles from borehole database; Subsurface sediment map 1:200,000 [43]; field works; geophysical data
Hydraulic conductivityC3Archival field tests in wells; PARAMEX and BAT tests in piezometers and shallow boreholes [26,29]
Table 3. Groundwater depth and its range in selected years.
Table 3. Groundwater depth and its range in selected years.
Groundwater Depth (m b.g.l.)Range of Parameter DArea of Occurrence in Selected Years (km2)
1999–201320032011
<110171.3562.57322.24
1–39266.34277.26176.21
3–5889.39157.3848.6
>5747.2677.1427.29
Table 4. Net recharge and its range in selected years.
Table 4. Net recharge and its range in selected years.
Recharge Rate (mm/year)Range of Parameter RArea of Occurrence in Selected Years (km2)
1999–201320032011
<−50124.8410.99120.09
−50–02139.74120.77146.44
0–20369.12107.0857.06
20–40477.47290.6955.69
40–605263.1832.4285.21
>606 12.39109.86
Table 5. Characteristics and range of A, S, T, I, C parameters.
Table 5. Characteristics and range of A, S, T, I, C parameters.
ParameterCharacteristicRangeArea of Occurrence (km2)
Aquifer media
A
Surface water01.22
Fine sands, silty sands with till inlays223.74
Fine sands with silty alluvial deposits3148.70
Sands with inlays of weathered tills436.10
River and fluvioglacial sands6213.39
Gravels and eolian sands8183.54
Soil media
S
Dusty, mineral-muck, silty alluvial soil5166.07
Brown soil, black soil6133.08
Eolian-erosive soil732.24
Silt‒gley, peat soil, muck soil8125.21
Ground‒gley, gley‒podsolic, podsolic soil9117.52
Topography
Slope (%)
T
>3.07.512.92
2.5–3.088.81
2.0–2.58.517.91
1.5–2.0923.57
1.0–1.59.563.78
0.0–1.010447.47
Impact of vadose zone
I
Surface water01.22
Laminated clays, tills223.74
Aggregate mud, silty alluvial deposits3148.70
Weathered tills436.10
River and fluvioglacial sands6213.39
Gravels, eolian sands8183.54
Hydraulic conductivity
(m/day)
C
<412.96
4–12287.15
13–284265.77
29–406148.65
41–80869.67
Table 6. Statistics of IPZ index and occurrence of vulnerability classes on research area.
Table 6. Statistics of IPZ index and occurrence of vulnerability classes on research area.
IPZ IndexVulnerability ClassVulnerability for Groundwater State Defined in Years:
1999–201320032011
Area (km2)
75–100Very low1.201.851.11
101–125Low138.80160.27141.11
126–150Medium193.21216.26197.99
151–175Medium high231.87190.75221.75
176–200High5.682.129.29
Mean IPZ142.7139.5142.99
Standard deviation19.3318.7720.54
IPZ scope80.5–18176.5–18080.5–185
Table 7. Original and effective weight of parameters.
Table 7. Original and effective weight of parameters.
Parameter“Theoretical” Weight “Theoretical” Weight (%)“Effective” Weight“Effective” Weight (%)
MeanMinimum Maximum Standard Deviation
D521.74 7.432.2820.4759.526.63
R417.392.410.372.4423.813.54
A313.042.711.58017.333.13
S28.702.29.595.7819.881.89
T14.351.66.994.311.91.15
I521.744.419.3028.885.21
C313.042.39.891.9423.763.71
Total2310023100
Table 8. IPZ changes after optimization by single parameter sensitivity analysis.
Table 8. IPZ changes after optimization by single parameter sensitivity analysis.
IPZ IndexVulnerability ClassClass Area after Optimization (km2)Class Area Changes according to Original IPZ (%)
Average Groundwater State 1999–2013Groundwater State in 2003Groundwater State in 2011Average State 1999–2013Groundwater State in 2003Groundwater State in 2011
<100Very low0.210.270.15−0.2−0.3−0.2
100–125Low5.3818.382.73−23.4−24.8−24.2
125–150Medium175.67186.68173.38−3.1−5.2−4.3
151–175Medium‒high314.32328.80287.3014.424.211.5
>175High75.6737.12107.6912.36.117.2

Share and Cite

MDPI and ACS Style

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

AMA Style

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 Style

Krogulec, 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 Style

Krogulec, 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

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop