Next Article in Journal
A Combined Raindrop Aggregate Destruction Test-Settling Tube (RADT-ST) Approach to Identify the Settling Velocity of Sediment
Previous Article in Journal
Simulating the Hydrologic Impact of Arundo donax Invasion on the Headwaters of the Nueces River in Texas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimizing Groundwater Monitoring Networks Using Integrated Statistical and Geostatistical Approaches

by
Jay Krishna Thakur
Environment and Information Technology Centre—UIZ, Neue Grünstraße 38, Berlin 10179, Germany
Hydrology 2015, 2(3), 148-175; https://doi.org/10.3390/hydrology2030148
Submission received: 11 April 2015 / Accepted: 22 July 2015 / Published: 24 August 2015

Abstract

:
The aim of this work is to investigate new approaches using methods based on statistics and geo-statistics for spatio-temporal optimization of groundwater monitoring networks. The formulated and integrated methods were tested with the groundwater quality data set of Bitterfeld/Wolfen, Germany. Spatially, the monitoring network was optimized using geo-statistical methods. Temporal optimization of the monitoring network was carried out using Sen’s method (1968). For geostatistical network optimization, a geostatistical spatio-temporal algorithm was used to identify redundant wells in 2- and 2.5-D Quaternary and Tertiary aquifers. Influences of interpolation block width, dimension, contaminant association, groundwater flow direction and aquifer homogeneity on statistical and geostatistical methods for monitoring network optimization were analysed. The integrated approach shows 37% and 28% redundancies in the monitoring network in Quaternary aquifer and Tertiary aquifer respectively. The geostatistical method also recommends 41 and 22 new monitoring wells in the Quaternary and Tertiary aquifers respectively. In temporal optimization, an overall optimized sampling interval was recommended in terms of lower quartile (238 days), median quartile (317 days) and upper quartile (401 days) in the research area of Bitterfeld/Wolfen. Demonstrated methods for improving groundwater monitoring network can be used in real monitoring network optimization with due consideration given to influencing factors.

1. Introduction

In regions with frequent water stress and large aquifer systems, groundwater is often used as an additional water resource. If extensive groundwater abstraction exceeds natural groundwater recharge for longer times, overexploitation or persistent groundwater depletion occurs [1].
Population growth and groundwater depletion are two of the most significant dangers to global water stability. It is evident from the trending water scarcity and continual groundwater pollution throughout the world that the existing policies fail to protect this vital resource [2,3].
In almost all major water resource management programs, groundwater monitoring is required [4]. Ground water monitoring involves an array of sequential processes ranging from long term standardized measurement, observation, status and trend evaluation to the final reporting of processed data in order to meet the objectives of monitoring programmes [5,6].
Monitoring networks with sets of strategically located groundwater monitoring wells collecting prerequisite data [7] that, in turn, provide insights for strategic planning and decision making [8]. Since complex infrastructure is crucial for entire data analysis and processing activities, ground water monitoring acts as a significant burden at many ground water contamination sites [9,10].
The last three decades has achieved great improvement in groundwater monitoring strategies, the majority of which have strongly focused on statistical approach for monitoring locations and sample size [11,12,13]. However, relatively little work is in place with optimization considering multiple contaminants in multiple aquifers.
For proper design of monitoring strategy as per management objectives, it is necessary that required information be extracted from measured variables [4,14]. Thereafter, the monitoring strategy can be derived based on optimization of the groundwater monitoring network and concurrent considerations of the uncertainty associated with spatial and temporal variability [15,16].
In the past, various studies have specified various methods for monitoring network optimization [17,18,19,20,21,22] but have failed to incorporate the primary factors that influence the methods for so. This study attempts to thoroughly analyze the influence of various factors on the optimization method, which can have a remarkable role in the decision about redundancy and necessity of new wells in the monitoring network.
The aim of this research is to identify and investigate new approaches using statistical and geostatistical methods for groundwater monitoring network optimization in order to improve groundwater monitoring strategies. The factors influencing groundwater monitoring optimization methods will be analyzed and new and improved methods will be integrated by testing those with groundwater monitoring network optimization at contaminated sites.

2. Study Area and Data Used

2.1. Study Area

In order to address the research questions and test the developed approaches and methods, Bitterfeld/Wolfen, located in the Federal State Saxony-Anhalt, Germany, was selected as a study area (Figure 1). Yet again, to make precise application of improved and developed methods in this study area, an area of about 100 km2 in urbanized zone of Bitterfeld/Wolfen (latitude 51°35′30″–51°41′30″ and longitude 12°14′10″–12°20′0.5″) consisting of a long term monitoring network (LTM) of 357 wells in the Tertiary and 462 wells in the Quaternary aquifer was selected.
Figure 1. Location map of the study area in Bitterfeld, Federal State of Saxony-Anhalt, Germany (a) Saxony-Anhalt in eastern Germany, (b) Federal State of Saxony-Anhalt, (c) study area of 100 km2 used for monitoring network optimization in Bitterfeld/Wolfen showing location of groundwater monitoring wells).
Figure 1. Location map of the study area in Bitterfeld, Federal State of Saxony-Anhalt, Germany (a) Saxony-Anhalt in eastern Germany, (b) Federal State of Saxony-Anhalt, (c) study area of 100 km2 used for monitoring network optimization in Bitterfeld/Wolfen showing location of groundwater monitoring wells).
Hydrology 02 00148 g001
In terms of geography, the western part of the study area is covered by glacial outwash sediments, whereas the flood plain of the Mulde river constitutes the eastern part of research area [23]. The geological setting of the research area reveals the presence of pre-tertiary rocks (overlaid by Cenozoic sediments), separated hydro-geologically from one another by clay layer at the depth of 50–70 m [24]. An upper Quarternary aquifer system, lower Tertiary aquifer system and pre-tertiary basement constitutes a geological cross section in the research area as portrayed by the Bitterfeld/Wolfen model Stollberg, Gossel et al. [25].

2.2. Nature and Availability of Data

An available data set of physicochemical parameters (α-Hexachlorocyclohexane (α-HCH), Monochlorobenzene (MCB), temperature, pH-value, oxidation potential (Eh), sulphate (SO42−), sulphite (SO32−), nitrate (NO3), ammonium (NH4+) and iron (Fe3+)) from the period 2003 to 2009 were used for LTM network optimization. Assessory information like name of the groundwater observation well, coordinates, elevation data of the monitoring well, screen depth, stratigraphic geological layer, stratigraphic horizon (Quaternary (Q), Tertiary (T) and Quaternary-Tertiary (Q-T)), the date and time of sampling and the name of the analysing laboratory (Table 1) were also available. In the table, the 9th column gives the total number of wells and samples during the overall period 2003–2009.
Table 1. Number of wells and samples from 30 September 2003 to 15 December 2009, showing number of wells monitored each year for physicochemical parameters.
Table 1. Number of wells and samples from 30 September 2003 to 15 December 2009, showing number of wells monitored each year for physicochemical parameters.
Year20032004200520062007200820092003–2009
No of well47757949666368252138827
No of sample796749729847711519384389
pH786719720865501519384148
Eh787701703841501519384090
NO3787709719866414519314045
SO32−796703700864424519314037
SO42−795729720866424519314084
NH4+717729720866420519314002
Fe2+0810218000128
Fe3+772709694841475519384048
α-HCH772682698823500519384032
MCB735678739215619542307209710,666

3. Methods

Design and improvement of monitoring strategies requires the description of its components. The components of groundwater monitoring strategy, analyzed for monitoring scenario in the study area, is represented in the logic diagram (Figure 2).
Among the groundwater physicochemical properties and contaminants concentration data, representative contaminants were first identified and then prioritized for inclusion in the monitoring strategy. Their associated uncertainties along with spatial and temporal variance were incorporated.

3.1. Selecting the Representative Variable

Identification of monitoring components in the research area was done through Principal Component Analysis (PCA). PCA of the normalized variables was performed to extract significant PCs. The variables were standardized through z-scale transformation to avoid misclassification due to difference in data dimensionality [26,27]. PCA provides information about the principle variables responsible for spatial variation in groundwater quality [28].
Figure 2. Logic diagram showing components of groundwater monitoring strategies as a circular continuous process.
Figure 2. Logic diagram showing components of groundwater monitoring strategies as a circular continuous process.
Hydrology 02 00148 g002
Measures of descriptive statistics in conjunction with statistical modelling of all available water quality and contaminant concentration data were used to identify the representative contaminants. World Health Organization (WHO) standards for drinking water and maximum contamination limit (MCL) of pollutants and contaminant concentration as well as associated human health risk were considered in the course of identification of the representative contaminants [29,30].

3.2. LTM Network Optimization

Necessary information about contaminants, including their properties, hydro-geochemical properties of aquifers, and other related information, are the pre-requisites for optimization of LTM network. Representative contaminants in the study area were MCB, SO42 and α-HCH, which was identified based on PCA that was used for monitoring network optimization. MCB and SO42− represent organic and inorganic contamination respectively in the research area [31,32]. α-HCH, which is an organochloride and one of the isomers of hexachlorocyclohexane (HCH), represents pesticides [33,34].

3.2.1. Spatial Optimization Methods

Redundant wells were predicted using Geostatistical Temporal-Spatial algorithm (GTS) on a concrete basis that the remaining wells represent the same amount of statistical information about the underlying plume [35,36]. GTS concept considers a well or its group redundant if it plays insignificant contribution to the interpolated map of contaminant plumes. The investigation steps involved in locating essential and redundant wells are shown in Figure 3.
An interpolated contaminant concentration plume was generated with input data of MCB, α-HCH and SO42− at specified locations of the well. The maximum concentration limit (MCL) for each contaminant i.e., 100 μg/L, 0.2 μg/L and 0.25 g/L respectively for MCB, α-HCH and SO42− [30] was given as input indicator limit to the statistical algorithm. The contaminant concentration was converted into an indicator value as
iC(t) = 1 if z(t) > C and 0 if z(t) ≤ C
where, iC(t) is the indictor value, z(t) is the concentration at time t, and C is the maximum contamination limit (MCL).
Figure 3. Research steps to locate essential and redundant wells and to propose new wells in the existing monitoring network (Modified from Cameron and Philip M. Hunter (2010) [36]).
Figure 3. Research steps to locate essential and redundant wells and to propose new wells in the existing monitoring network (Modified from Cameron and Philip M. Hunter (2010) [36]).
Hydrology 02 00148 g003
Once again, to take into account the correlation based on distance and direction between pairs of sampling locations for each contaminant, the empirical variogram was modelled. The empirical spatial variogram, γ, as a function of lag distance can be expressed as:
γ ( Δ h ) = 1 2 N ( Δ h ) i . j | | x j x i | Δ h | | ε [ i C ( x j ) i C ( x i ) ] 2
where Δh is the targeted lag distance and may depend on direction; xi and xj are the ith and jth monitoring well locations; and N(Δh) is the number of indicator pairs contributing to the summation for lag Δh.
Numerical weight was given to the wells based on global kriging and a plume map was generated for the monitoring network. For this process, two intermediate computations were carried out:
(i) the local kriging weights of sampled locations which were processed to generate a “global” interpolation weight for each well [37]; and (ii) the local kriging estimation variance which indicate the relative uncertainty of the local block estimates in comparison to the estimates at other blocks [36].
Spatial-temporal algorithm incorporates series of procedural steps to estimate global interpolation weight and the local kriging estimation variance. At first, the study area is divided into sets of non-overlapping blocks. The set of closest sampled locations is then located at each block with the help of simple search algorithm. Then, local kriging weights (λB) is computed for each well based on the spatial configuration of known indicator value within and surrounding the block as well as spatial correlation between average block location and each known indicator by using modelled spatial covariance function.
A block indicator estimate is then generated with the combination of local weights and indicator values. A block indicator estimate consist of weighted average of the n(xB) indicators located within the search radius of the block (xB), which is estimated as [36]:
i C ( x B ) = i = 1 n ( x B ) λ i B i C ( x i )
for xi ϵ search radius of block xB.
Global interpolation weights (λG) is generated by averaging the local kriging weights to the given well that can be used to estimate the well’s overall contribution to the interpolated map [36], and are given by:
λ G ( x i ) = 1 N B B = 1 N B λ i B ( x i )
where NB is the total estimated number of blocks and xi is the location of the ith sampled well.
The global kriging weights provide relative rankings of well locations in terms of independent spatial information which gives a logical interpretation that the wells with the lowest global kriging weights are predominantly spatially redundant wells. If the temporary removal of the subsets of these wells does not significantly affect the nature of plume map, these wells were permanently removed. This process is repeated until the removal of a subset of wells change the plume map. Performing this process subsequently creates a situation where the removal of subsets significantly changes the plume. At this point in time, the subset is not removed from the network. In order to fix the limitation for such changes in plume map, a two-sided (upper and lower) confidence interval of 95% is assigned while limiting median plume concentration.
Again, the relative spatial uncertainty for the installation of new monitoring wells in the existing LTM network based on the local kriging variance is given by the global kriging variance (kvG), defined as:
kv G   =   1 N B B = 1 N B k v B ( x B )
where xB denotes the location of the Bth block and kvB (xB) is the local kriging variance of the Bth block [36].
The groundwater LTM network has been spatially optimized for MCB, α-HCH and SO42− in 2-D and 2.5-D of the groundwater aquifer. The inherent assumption in 2-D analysis is the existence of all well locations in flat 2-D which is most applicable with a single, fairly uniform and well-connected aquifer. In contrast, the 2.5-D analysis assumes multiple hydrogeological layers in the aquifers without hydraulic interconnection in which LTM is optimized for separate hydrogeological layers. The Quadratic Logistic Regression (QLR) mapping algorithm then uses the data (segregated into subsets representing one chemical of concern for each vertical zone and time slice triplet) from a given subset to map the layer and time frame represented by that given triplet.

3.2.2. Temporal Optimization Method

To optimize groundwater monitoring network temporally, Sen’s method [38] was applied which is based on iterative thinning. The algorithm consists of (i) estimating a trend using the entire time series; (ii) thinning the time series by a randomly selected fraction of the measurements; and then (iii) re-estimating the trend to determine if the slope estimate is still close to the original slope. The random removal is only retained if the measured slope is within the limits of the confidence interval on the slope using the full dataset [36]. To eliminate the effect of statistical assumptions inherent in standard linear regression methods, trend estimation is carried out using Sen’s method [38,39]. Sen’s procedure is non-parametric, readily adapted to non-detect measurements and irregular sampling frequencies [40]. Furthermore, this method precludes the calculation of slope for missing points and can predict the median slope even if the number of undetected measurements is less than (n − 1)/2.
Sen’s estimate (Q) is simply the median value of the resulting list of slopes and is given by:
Q = x j x i j i
Q is the slope between data points; xi and xj are concentrations measured at times i and j; Time j is after time i (j > i).
Q = median slope = Q[(N + 1)/2] if N is odd = (Q[N/2] + Q[(N + 2)/2])/2 if N is even
where N is the number of calculated slopes
A two side (M1, upper and M2, lower) confidence interval for the median slope is estimated using Zstatistic and Mann-Kendall statistic (VAR(S)). If there is a two-sided confidence interval of 95%, the Z(1−0.05/2) = Z0.975 = 1.96. Mann-Kendall statistic (VAR(S)) [41,42] is given by:
V A R ( S ) = 1 18 [ n ( n 1 ) ( 2 n + 5 ) p 1 q t p ( t p 1 ) ( 2 t p + 5 ) ]
where n is number of sampling data points, tp is the number of ties for the pth value and q is the number of tied values. Equation (8) may be used for values between 10 and 40. The range of ranks for the specified confidence interval (Ci) [38] is given as shown below:
C i = Z 1 i / 2 × V A R ( S )
Taking the value of Equation (9), the ranks of the lower (M1) and upper (M2 + 1) confidence limits can be calculated using the following relation:
M 1 = N C i 2   And  M 2 = N + C i 2
The values of lower (M1) and upper (M2 + 1) confidence limits were used to define lower and upper boundaries along the median slope. The temporal optimization of the existing groundwater monitoring network was carried out using α-HCH, MCB and SO42− concentration data set from 2003–2009 for each contaminant separately and together.

3.3. Factors Influencing Monitoring Optimization Results

3.3.1. Block Width and Dimension Dependency

In the monitoring network optimization using geostatistical methods, a plume map of well location is created with the help of global kriging weight [8] computed by interpolation method which depends upon width and number of blocks. Hence, in the optimization process, a block width from 1000 m to 1 m was demarcated in order to discover ambiguities in the method.
The 2-D treatment ignores the sampling depth of monitoring well while in the 2.5-D case, the aquifer is assumed to have the number of hydrogeological layers and thus the sampling depth of the monitoring wells needs to be taken into account. For systematic study of the influence of dimension and block width on monitoring network optimization methods, the optimization process has considered block widths from 1000 m to 1 m in both 2 and 2.5 dimensions.

3.3.2. Contaminants Association

Based on PCA of hydro-geochemical data, contaminant concentration of three representative variables, α-HCH, MCB and SO42− are used individually to monitor network optimization. The likely influence of multiple variables on optimization process is examined using data sets in groups as well as all three variables together. In this case, co-kriging was used to compute co-kriging numerical weights from contaminant concentrations for each monitoring well which in turn were used to create a plume map of the monitoring network which is used for geostatistical optimization of monitoring network as explained in Section 3.2.1.

3.4. Groundwater Flow Direction and Aquifer Homogeneity

In addition to spatio-temporal optimization of ground water monitoring network, direction of ground-water flow and homogeneity of the aquifer was determined. Groundwater flow direction is known to play a prominent role in monitoring network optimization. The directional dependency was studied using directional variograms. Estimation of experimental variogram was carried out using Golden Surfer software [43].
Since variograms based on contaminant concentration data set were anisotropic in nature, valid variogram model incorporating directional dependence were constructed using standard models to model the data set. The best fitting model, considering range, sill, nugget effect and shape of the model, was selected. A variogram was built from a concentration data set of α-HCH, MCB and SO42− for the years 2003–2009, according to seasons (summer: May–October; winter: November–April) and to hydrologic seasons (March–May: High groundwater level, September–November: Low groundwater level). The range and sill were estimated for various directional variogram models.
A homogeneity index (RV index) which numerically estimates homogeneity of the aquifer was defined based on range, sill and variance of the variogram [44] based on the variogram model as given by:
Homogeneity index = range s i l l + v a r i a n c e 2

3.5. Integrating Approaches and Improving Groundwater Monitoring Strategy

New sets of methods were analyzed with new approaches based on integrated statistical and geostatistical methods, integrated in the light of various optimization objectives. However, basic framework and components of groundwater monitoring strategies were duly taken into consideration while subjecting these methods to integration. The new integrated methods, although based on different approaches, could be used as an impeccable tool for groundwater monitoring network optimization and feedback for updating groundwater monitoring strategies. While integrating these methods, the advantages and disadvantages of the methods being used have also been analysed.

4. Results

4.1. Representative Variable in the Research Area

For the purpose of dimension reduction and analysis of major components loading the system, PCA was applied to a matrix of hydro-geochemical data. The first four principal components in this dataset with eigenvalues greater than one explain 66% of the cumulative variance. The significance of the parameters in four components was further depicted by the component matrix extracted with PCA and varimax with the Kaiser Normalization rotation method. At 75% significance level, SO42−, Fe, pH and Eh play a major role in the system whereas α-HCH and MCB are minor components. In the first component, SO42− and Fe have positive loading. In the second component, Eh has positive loading. However, pH is has negative loading because the acidic groundwater environment in the study area arises from a high SO42− concentration. SO42− and Fe have positive loading in the first component whereby Eh has positive loading in the second component. The negative loading of pH could be linked to the acidic groundwater environment in the study area from high SO42− concentration.

4.2. Optimization of LTM Network

4.2.1. Spatial Optimization of LTM Network

Both types of aquifers were considered separately during the spatial optimization of LTM network in accordance with the methods described in Section 3.2. The optimized monitoring network suggests 292 wells out of 462 wells in the Quaternary aquifer and 256 wells out of 357 wells in the Tertiary aquifer for monitoring needs at the suggested temporal interval.
The spatial uncertainties analysis, resulting from the global kriging variance (using Equation (5)) made recommendations for 22 and 41 new monitoring wells, respectively, in the Tertiary and Quaternary aquifers. The spatial distribution of redundant, essential, and proposed new wells in the Quaternary and Tertiary aquifer is depicted in Figure 4. Similar kriging weights in the original and optimized LTM network datasets points out that reduction in observation points does not compromise the quality or resolution of the collected samples as long as the distribution is logically and conveniently designed.
Figure 4. Optimized LTM network map showing location of essential, redundant, and proposed new wells in the monitoring network.
Figure 4. Optimized LTM network map showing location of essential, redundant, and proposed new wells in the monitoring network.
Hydrology 02 00148 g004

4.2.2. Temporal Optimization of LTM Network

Monitoring network was temporally optimized through statistical approach using Sen’s method, with 95% confidence intervals around the slope estimates, as described in Section 3.2.2. The outputs from temporal optimization are the sampling intervals at 3 quartiles (lower, median and upper) for individual monitoring well and each monitoring parameter (Table 2). The optimization results reveal that the optimization differs remarkably while considering pair and multiple contaminants.
Table 2. Temporal optimization of monitoring network in Quaternary and Tertiary aquifer for monochlorobenzene (MCB), α-HCH and SO42−.
Table 2. Temporal optimization of monitoring network in Quaternary and Tertiary aquifer for monochlorobenzene (MCB), α-HCH and SO42−.
V. ZoneCOCPresent Sampling Interval (Days)Recommended Sampling Interval (Days)
Lower QuartileMedian QuartileUpper QuartileLower QuartileMedian QuartileUpper Quartile
QMCB92210254205328503
Qα-HCH138181224298327418
QSO42−179217268325506704
TMCB143221280282443810
Tα-HCH224224224342398553
TSO42−188217278435589836
Qall138210254298328503
Tall188221278342443810
Bothall161217261311420628
Note: Q: Quaternary aquifer; T: Tertiary aquifer; COC: Chemical of Concern.
In order to clearly visualize and precisely describe the temporal optimization results of three contaminants in combination, the recommended median quartile sampling frequency for the monitoring wells was divided into five classes: three months, six months, one year, two years and three years. Recommendations for specific number of monitoring wells with respect to each of the categorized temporal sampling interval are specified in Table 3.
Table 3. Number of monitoring wells for each sampling interval.
Table 3. Number of monitoring wells for each sampling interval.
Sampling Interval3 Months6 Months1 Year2 Years3 YearsTotal
Quaternary34861737693462
Tertiary16691148474357
Total No. of wells50155287160167819
Figure 5 shows the distribution of monitoring wells and their recommended sampling intervals. The highest number of sampling wells is recommended at the yearly sampling interval.
In the study area, the overall temporally optimized sampling interval was recommended as 238 days, 317 days and 401 days in terms of lower, median and upper quartiles respectively.

4.3. Factor Influencing Monitoring Optimization Results

4.3.1. Dimension and Block Width Dependency

With an objective of incorporating dependency of dimension and block width on optimization process, the optimization was carried out for different interpolation block widths (from 1000 m to 1 m) for both aquifers in both 2 and 2.5 D dimensions as described in Section 3.3.1.
Figure 5. Statistically temporally optimized LTM network map showing recommended temporal frequency of the monitoring wells in the monitoring network.
Figure 5. Statistically temporally optimized LTM network map showing recommended temporal frequency of the monitoring wells in the monitoring network.
Hydrology 02 00148 g005

Optimization in 2-D Aquifer

Considering the aquifer as a 2-D plain, groundwater monitoring network was spatially optimized for MCB, α-HCH and SO42. The result unveiled a significant change in number of suggested/essential wells as well as existing redundant wells with changing block width for interpolation (from 1000 m to 1 m). The increase in the relative spatial uncertainties based on local kriging variance with decreasing block width eventually recommends installation of new monitoring wells.
The overall analysis of optimization in 2-D plain unveils that the relative spatial uncertainty increases with decreasing block width from 200 m to 1 m in both aquifers (Figure 6 and Figure 7), thus suggesting the need of installation of 63 and 36 new wells, respectively, in Quaternary and Tertiary aquifers at 1 m block width.
Figure 6. Block width dependency in the LTM network optimization in the Quaternary aquifer for MCB, α-HCH and SO42− in a 2-D plain. Essential, redundant and total number of wells is on left Y-axis and number of new wells is on right Y-axis.
Figure 6. Block width dependency in the LTM network optimization in the Quaternary aquifer for MCB, α-HCH and SO42− in a 2-D plain. Essential, redundant and total number of wells is on left Y-axis and number of new wells is on right Y-axis.
Hydrology 02 00148 g006
Figure 7. Block width dependency in the LTM network optimization in the Tertiary aquifer for MCB, α-HCH and SO42− in a 2-D plain. Essential, redundant and total number of wells are on left Y-axis and the number of new wells is on right Y-axis.
Figure 7. Block width dependency in the LTM network optimization in the Tertiary aquifer for MCB, α-HCH and SO42− in a 2-D plain. Essential, redundant and total number of wells are on left Y-axis and the number of new wells is on right Y-axis.
Hydrology 02 00148 g007

Optimization in a 2.5-D Aquifer

Spatial optimization for MCB, α-HCH and SO42− in groundwater LTM network has been conducted in combination in 2.5-D aquifer. The LTM network is optimized separately for each hydro-stratigraphic layer in the aquifer using data from that layer only.
In the case when all three contaminants are considered together, the change in block width for interpolation (1000 m to 1 m) does not significantly change the number of suggested wells (Figure 8 and Figure 9). The increment of spatial relative uncertainty significantly with decreasing block width from 20 m to 1 m suggests the installation of 58 and 38 new monitoring wells in Quaternary and Tertiary aquifers, respectively.
Figure 8. Block width dependency in the LTM network optimization in the Quaternary aquifer for MCB, α-HCH and SO42− considering 2-D and 2.5-D aquifers. Essential, redundant and total number of wells are on left Y-axis and number of new well is on right Y-axis.
Figure 8. Block width dependency in the LTM network optimization in the Quaternary aquifer for MCB, α-HCH and SO42− considering 2-D and 2.5-D aquifers. Essential, redundant and total number of wells are on left Y-axis and number of new well is on right Y-axis.
Hydrology 02 00148 g008
Figure 9. Block width dependency in the LTM network optimization in the Tertiary aquifer for MCB, α-HCH and SO42− in 2-D and 2.5-D aquifers. Essential, redundant and total number of wells are on the left Y-axis and the number of new wells is on the right Y-axis.
Figure 9. Block width dependency in the LTM network optimization in the Tertiary aquifer for MCB, α-HCH and SO42− in 2-D and 2.5-D aquifers. Essential, redundant and total number of wells are on the left Y-axis and the number of new wells is on the right Y-axis.
Hydrology 02 00148 g009

4.3.2. Contaminants Association

To account for the association between the contaminants, the LTM network was spatially optimized taking 7 m interpolation block width for the contaminants MCB, α-HCH and SO42− in pairs as well as all grouped together, considering the groundwater aquifer separately as 2-D plain and as 2.5-D aquifer. Dependency of contaminants association on network optimization is supported by result that the spatial optimizations of the LTM network for MCB and α-HCH individually show less redundancy and low recommendations for new monitoring wells while that only for SO42− recommends higher number of redundant wells (70%). Such effect of incorporating greater number of contaminants is presented in Table 4 and Table 5. For example, in the Quaternary aquifer considered as a 2-D plain, for α-HCH and SO42−, 73 wells are recommended, for MCB and SO42, 293 wells are recommended and for MCB and α-HCH, 335 wells are recommended. However, 244 wells are recommended when considering all three contaminants MCB, SO42− and α-HCH together (Table 4). Similarly, in the tertiary aquifer considered as a 2.5-D plain, 72, 267, 247 and 213 wells are recommended, respectively, for combination of α-HCH and SO42−; MCB and SO42−; MCB and α-HCH; and MCB, SO42− and α-HCH together (Table 5).
Table 4. Number of essential, redundant and new well location recommendation from 2-D groundwater monitoring network optimization for α-HCH and SO42; MCB and SO42−; MCB and α-HCH; and MCB, SO42− and α-HCH in Quaternary (462 wells) and Tertiary (357 wells) aquifers.
Table 4. Number of essential, redundant and new well location recommendation from 2-D groundwater monitoring network optimization for α-HCH and SO42; MCB and SO42−; MCB and α-HCH; and MCB, SO42− and α-HCH in Quaternary (462 wells) and Tertiary (357 wells) aquifers.
Quaternary Aquifer/Contaminantsα-HCH and SO42−MCB and SO42−MCB and α-HCHMCB, SO42− and α-HCH
Essential well73293335244
Redundant well351169127218
New well location1471414
Tertiary Aquifer
Essential wells52233247220
Redundant well302124110137
New well location0000
Table 5. Number of essential, redundant and new well location recommendation from 2.5-D groundwater monitoring network optimization for α-HCH and SO42; MCB and SO42− MCB and α-HCH; and MCB, SO42− and α-HCH in Quaternary (462 wells) and Tertiary (357 wells) aquifers.
Table 5. Number of essential, redundant and new well location recommendation from 2.5-D groundwater monitoring network optimization for α-HCH and SO42; MCB and SO42− MCB and α-HCH; and MCB, SO42− and α-HCH in Quaternary (462 wells) and Tertiary (357 wells) aquifers.
Quaternary Aquifer/Contaminantsα-HCH and SO42−MCB and SO42−MCB and α-HCHMCB, SO42− and α-HCH
Essential well107334341298
Redundant well317127140164
New well location02114
Tertiary Aquifer
Essential well72267247213
Redundant well28290148144
New well location0385
In the optimization of LTM network for combination of three contaminants, the local kriging weight of each contaminant is averaged, which clearly depicts the dependence of relative spatial uncertainty with the spatial distribution of each contaminants. It has been observed that optimization for each type of aquifer considering individual contaminants and in groups gives recommendations for different numbers and locations of new wells (Table 4 and Table 5).

4.4. Groundwater Contaminants Flow Direction and Aquifer Homogeneity

Geostatistical and flow direction modeling methods have been used to predict the groundwater flow direction and integrate its effect into LTM network optimization. The spatial characterization of the groundwater contamination scenario (degree of spatial co-relation between contaminant concentration values) was observed using an experimental variogram exhibiting the contaminant concentration data of MCB, α-HCH and SO42−.
The geometric anisotrophy as depicted by the experimental variogram was observed for each 30° of lag direction with reference to North direction. Again, the heterogeneity in the nature of geological formation was quantified by using a homogeneity index called RV Index (Equation (11)), which takes its base from spatial variability of contaminant concentration distribution.

4.4.1. Contaminant Wise Variogram Modelling

Based on the results from experimental variogram, geometrical anisotropy was observed in the data set of α-HCH, as range and sill differed in different directions.
The experimental variogram (incorporating each 30° of lag direction) revealed the highest range in Northern direction (0°) which is also coherent with the results from RV index that was highest in the Northern direction for Quaternary aquifer (Table 6). Meanwhile, the variogram showed highest range and RV index at 30° from the Northern direction for the Tertiary aquifer (Table 6). These results specify that the overall conspicuous concentration of α-HCH flowed towards the North in both aquifers. In case of Tertiary aquifer, there was slight deviation in the flow towards east direction.
Table 6. Directional variogram modelling α-HCH concentration in the long term monitoring network (LTM) network in Quaternary and Tertiary aquifers (January 2003 to February 2009).
Table 6. Directional variogram modelling α-HCH concentration in the long term monitoring network (LTM) network in Quaternary and Tertiary aquifers (January 2003 to February 2009).
α-HCH for Quaternary Aquifer (2003–2009)α-HCH for Tertiary Aquifer (2003–2009)
Model: Gaussian, Variance: 0.62Model: Gaussian, Variance: 0.57
DirectionRangeSillRV IndexDirectionRangeSillRV Index
Omni4040.64641.27Omni11820.701862.88
06660.621074.1909600.791416.97
304640.64736.513013400.672161.29
603850.66601.56609990.691591.91
903210.66501.56907100.701118.11
1203600.65566.931207220.701137.01
1504330.63692.801507500.681200.00
1805900.62951.611809600.691529.76
The relatively lower RV index in the Quaternary aquifer than in the Tertiary aquifer shows high heterogeneity in the Quaternary aquifer.

4.4.2. Season Wise Variogram Modelling for α-HCH Data

Seasonal variability was embodied into optimization by using the data set of α-HCH for hydrological summer and winter seasons to model the experimental variogram. Periodic data from May to October (2003–2009) and November to April (2003–2009) were categorized as hydrological summer and winter seasons, respectively. Analysis of the experimental variogram, using α-HCH data for the hydrological summer season, again showed the highest range in the North direction as well as highest RV index in Northern direction for both Quaternary and Tertiary aquifers (Table 7 and Table 8). These results show that preferential groundwater flow and α-HCH concentration movement occurred in the Northern direction.
Table 7. Directional variogram modelling for α-HCH concentration in the LTM network in Quaternary aquifers for summer and winter seasons (May to October in 2003–2009).
Table 7. Directional variogram modelling for α-HCH concentration in the LTM network in Quaternary aquifers for summer and winter seasons (May to October in 2003–2009).
α-HCH for Quaternary Aquifer (May to October in 2003–2009)α-HCH for Quaternary Aquifer in Winter Seasons (Nov. to April from 2003–2009)
Model: Gaussian, Variance: 0.70Model: Gaussian, Variance: 0.68
DirectionRangeSillRV IndexDirectionRangeSillRV Index
Omni4910.78661.99Omni6690.73954.01
04910.70701.4303980.74562.54
303870.75532.54305340.71771.12
603400.77461.52605430.68801.48
903540.78477.28903970.72570.20
1203640.79488.591204110.73585.47
1503330.76456.161504220.74596.89
1803900.65577.781805130.73729.21
Table 8. Directional variogram modelling of α-HCH concentration in the LTM network in Quaternary and Tertiary aquifers for winter seasons (November to April from 2003–2009).
Table 8. Directional variogram modelling of α-HCH concentration in the LTM network in Quaternary and Tertiary aquifers for winter seasons (November to April from 2003–2009).
α-HCH for Tertiary Aquifer (May to October in 2003–2009)α-HCH for Tertiary Aquifer in Winter Seasons (Nov. to April from 2003–2009)
Model: Gaussian, Variance: 0.70Model: Gaussian, Variance: 1.25
DirectionRangeSillRV IndexDirectionRangeSillRV Index
Omni11660.412854.69Omni5851.41439.85
012580.472875.4305791.31452.34
3010620.462455.49305501.25440.00
609390.382392.36606001.33465.12
908880.432122.12908501.26677.29
1207370.431761.2612010501.22850.20
1507930.421934.151508231.25658.40
1808820.442082.891805261.32409.34
In the next round, the experimental variogram was modelled using α-HCH data for the hydrological winter season (November to April, 2003–2009). The experimental variogram showed the highest range in the direction of 60° from North in the Quaternary aquifer (Table 8) whereas the range was highest in the direction of 120° from North in Tertiary aquifer. Further, the RV index is lower in the Quaternary aquifer than the Tertiary aquifer, thus indicating high heterogeneity in Quaternary aquifer (Table 8).

4.4.3. Year Wise Directional Variogram Modelling for α-HCH Data

To acquaint with the geometrical anisotropy of the data set, variogram modelling was carried out for the monitoring data sets of the years 2005 and 2006. To understand the year wise alterations in the homogeneity of aquifer and flow direction, variogram modelling was carried out for the monitoring data sets of the years 2005 and 2006. The experimental variogram showed that the range was highest in the Northern direction, as shown in Table 9 and Table 10.
Again, the range was highest at 30° from North in the Tertiary aquifer. The RV index varied with different lag directions, showing no distinct pattern or peak.
Table 9. Directional variogram modelling α-HCH concentration in the LTM Network in Quaternary and Tertiary aquifers for 2005.
Table 9. Directional variogram modelling α-HCH concentration in the LTM Network in Quaternary and Tertiary aquifers for 2005.
α-HCH for Quaternary Aquifer in 2005α-HCH for Tertiary Aquifer in 2005
Model: Gaussian, Variance: 0.89Model: Gaussian, Variance: 1.05
DirectionRangeSillRV IndexDirectionRangeSillIndex
Omni8500.89955.06Omni7501.33631.31
05000.93550.9609371.30797.45
305330.92588.303011251.35937.50
607580.95826.16609331.30794.04
908560.93939.63907501.30638.30
1208000.94874.321207501.30638.30
1505330.95580.931505621.32474.26
1805000.91555.561806871.33577.31
Table 10. Directional variogram modelling α-HCH concentration in the LTM Network in Quaternary and Tertiary aquifers for 2006.
Table 10. Directional variogram modelling α-HCH concentration in the LTM Network in Quaternary and Tertiary aquifers for 2006.
α-HCH for Quaternary Aquifer in 2006α-HCH for Tertiary Aquifer in 2006
Model: Gaussian, Variance: 0.59Model: Gaussian, Variance: 1.05
DirectionRangeSillRV IndexDirectionRangeSillRV Index
Omni3750.65605.23Omni13330.811904.29
04680.61793.22016840.722405.71
304370.68740.683017270.842467.14
603750.67635.596013330.861904.29
903750.67635.599010580.841511.43
1203750.66635.5912010580.841511.43
1503750.66635.5915010550.681507.14
1803750.53635.5918012770.691824.29

4.5. Integrating Approaches for Improving Groundwater Monitoring

As an important objective of this research, statistical and geostatistical methods were taken as a basis for integrating new approaches to illustrate several sets of methods for different optimization objectives. Figure 10 depicts integration of statistical and geostatistical methods. Such integration was executed to fully understand, evaluate and optimize groundwater monitoring with reference to groundwater monitoring of a contaminated site.
Multivariate and descriptive statistics were the primary element of groundwater monitoring framework used for analysis, classification, modelling and interpretation of large volume of available data set. The second element consisted of optimization of existing LTM networks based on target variables (MCB, α-HCH and SO42−) using a geostatistical spatial optimization algorithm (GTA). GTS itself accounted for dimension dependency, influence of block width for interpolation and the influence of multiple contaminants on the LTM network optimization, whereas variogram modelling estimated the groundwater flow direction and heterogeneity of aquifers. As third and fourth elements of the groundwater monitoring framework, the obtained results were analysed and interpreted incorporating the prominent effects of other influencing factors such as legal requirements and land use changes so that essential, redundant and new monitoring wells in the existing LTM network could be recommended.
Figure 10. Research steps showing the integration of statistical and geostatistical methods for groundwater monitoring.
Figure 10. Research steps showing the integration of statistical and geostatistical methods for groundwater monitoring.
Hydrology 02 00148 g010
Strategic planning of monitoring network optimization is crucial to groundwater monitoring since it is an important component of monitoring strategies [45,46]. With an overall objective to improve the monitoring strategies, the developed and existing methods were applied to the observed and model-based data sets for the mega-contaminated site, Bitterfeld/Wolfen.
The existing monitoring strategies can be improved by incorporating the optimization outcomes from integrated methods along with functionality consideration of expert knowledge pertaining to the study area.

5. Discussion

PCA technique is an effective technique for pattern recognition which attempts to explain the variance of a data set of intercorrelated variables with small set of independent variables [47]. Representative variables for groups of groundwater contaminants (SO42−, α-HCH and MCB) were selected for monitoring network optimization in the study area based on these PCA results.
The optimization results from geostatistical methods suggests the need for monitoring for only 292 out of 462 wells in the Quaternary aquifer and 256 out of 357 wells in the Tertiary aquifer (Figure 4). Furthermore, the spatial uncertainty analysis in terms of kriging variance recommends 22 and 41 new monitoring wells to be added to the existing network in the Tertiary and the Quaternary aquifers, respectively, when three representative contaminants (α-HCH, MCB and SO42−) are considered together. Given the presence of variety of highly mobile and persistent contaminants in the research area [48,49], monitoring network optimization with more number of available contaminants would result in better optimization (as presented in Section 3.3.2. and Section 4.3.2).
Temporal optimization was performed with due consideration of pairs of contaminants as well as three contaminants together; recording the intervals in terms of lower quartile, median quartile and upper quartile for each monitoring well for each parameter. The result shows that the optimization differs markedly for single and multiple contaminants (Table 2).
Optimization in the 2-D plain considering three contaminants together shows the increase in relative spatial uncertainty with decreasing block width form 200 m to 1 m in both aquifers (Figure 6 and Figure 7) suggesting installation of 63 and 36 new wells, respectively, in Quaternary and Tertiary aquifers at 1 m block width. Again, with the same consideration of contaminants in 2.5-D analysis, the change in block width for interpolation (1000 m to 1 m) has an insignificant effect on the number of suggested wells (Figure 8 and Figure 9). In this case, the increment of spatial relative uncertainty significantly with decreasing block width from 20 m to 1 m suggests the installation of 58 and 38 new monitoring wells in Quaternary and Tertiary aquifers, respectively. A 2-D analysis treats all wells in a manner that they exist in a flat 2-D plane which should be more applicable in case of single, fairly uniform and well-connected aquifer [50]. On the other hand, a 2.5-D analysis treats the aquifer as comprising multiple Hydrogeological layers, each of which is considered separately during optimization within the GTS [36]. Again, the analysis reveals notable dependence of optimization of existing monitoring network on block width. With the gradual decrease in the block width, the recommendation of the number of new monitoring wells increase. The incorporation of the relatively influential role of block width for interpolation has not been documented in previous studies based on geostatistical methods [13,16]. Hence, this study has pioneered the incorporation of this factor in the monitoring network optimization process. Previous studies have, however, addressed problems related to scale in hydrology such as the works of Bergström and Graham [51], Dooge [52], Gupta et al. [53], Sivapalan and Kalma [54] and Sivapalan et al. [55]. The results based on the present study gives recommendations of scale size of 100 m to 7 m for monitoring network optimization at local scale.
The role of contaminants association was pronounced when three representative contaminants were considered in different pairs or all together (Table 4 and Table 5). The optimization result in terms of the number of essential, redundant and new wells was found to be different. Spatial variograms for each contaminant was first estimated for the spatial covariance model, which was determined through a non-linear modelling program. This modelling utilized the Levenberg-Marquardt algorithm [56,57]. The fitting algorithm was set up to fit either a combination of up-to-three spherical, exponential and/or Gaussian components or a combination of up-to-three power model structures [36]. As the spatial variogram for each contaminant varies, the combination of two or more contaminants determines the covariance of the contaminants for the kriging interpolation.
The analysis of flow direction (from January 2003 to February 2009) shows that α-HCH flows northwards in the Quaternary aquifer, whereas in the direction 30° from north in the Tertiary aquifer. This distinction in the flow direction with respect to aquifer type is the result of historic shift in flow direction from northwards to eastwards and then back to northwards. This represents source specific spreading direction of the α-HCH contaminant. The seasonal alterations in the groundwater flow system has also been documented in previous studies [58,59]. Seasonal influence was observed in both aquifers since the experimental variogram modelling using the α-HCH data shows north as the preferential flow direction in hydrological summer season (May to October in 2003–2009) while during hydrological winter season, the flow direction was 60° from north in the Quaternary aquifer and 120° from north in the Tertiary aquifer (Table 8). Seasonal fluctuation in the water level in Mulde river and surrounding water bodies [60] and the location of contaminant source strongly control such seasonal variation in the flow direction.
The geometric anisotropy in contaminant concentration data set using experimental variogram was the means of determining groundwater and contaminants flow direction (described in Section 3.4). Table 6, Table 7, Table 8 and Table 9 show the prominent groundwater and contaminant flow direction that are revealed from the data sets of α-HCH, MCB and SO42− concentration. It has been shown that although the prominent groundwater and contaminant flow direction is northwards, various contaminants spread non-uniformly in varied time scale. The RV index points out the hydrogeological homogeneity in the aquifer in the direction where there is the highest range. The likely flow direction and homogeneity analysis in a year wise basis displayed the complete change in flow direction in 2006 which was towards east in Quaternary aquifer, but in the direction of 30° from north in the Tertiary aquifer during 2005. The preferential flow direction was northwards in the Quaternary aquifer and 30° from north in the Tertiary aquifer.
Groundwater monitoring network optimization methods are greatly challenged by other factors besides the problems associated with nature and availability of data. Given the existence of a wide range of methods for LTM network optimization, most of them fail to implicate the factors that influence the optimization methods [7]. New and improved approaches based on integrated statistical and geostatistical methods for optimization of monitoring network along with their possible implications are dealt with in this study. In addition, a range of deterministic factors and their influence on application of these methods in the research area are analyzed. Existing monitoring strategies can be re-shaped by removing the existing redundancy as well as adding requisite newer monitoring wells. A new approach has been devised in this study by integrating statistical and geostatistical methods for spatio-temporal optimization of the LTM network which has been tested in the study area. The outcomes of this research can be taken as a basis for improving the strategies to make the long term monitoring program more efficient and equally effective. However, expert knowledge and site-specific peculiarities are important to take into consideration in order to implement this approach.

6. Conclusions and Recommendation

The results have verified that the geostatistical method, based on the kriging interpolation weight, shows realistic optimization results in terms of recommended essential, redundant and new monitoring wells. For the time being, several factors influencing optimization results such as block width, multiple contaminants association, and contaminant spreading direction were analyzed. The results obtained from such analysis revealed their significant role in the monitoring optimization process. Meanwhile, the network was temporally optimized using a statistical approach with the help of Sen’s method. Spatio-temporal optimization reduces the number of samples and optimizes the monitoring frequency, which could make the monitoring program more effective in terms of resource requirement. The flexibility inherent in both statistical and geostatistical approaches allow them to be operated at various confidence level limit set by user. Furthermore, these methods could be applied to field based real groundwater quality data set.
Since various methods have strengths and weaknesses, the strength of these methods can be integrated on the basis of optimization objectives of groundwater monitoring strategies. For example, in a monitoring network where there is a low density of monitoring wells, statistical and geostatistical methods can be integrated. An understanding of the distribution of contaminants, component analysis can be gained, as well as temporal optimization can be done with the help of statistical methods, whereas spatial optimization of the monitoring network can be done by using geostatistical methods. Temporal optimization of the monitoring network can be carried out using statistical methods. In this manner, these methods can be systematically integrated in monitoring network optimization according to the prime objective of the groundwater monitoring strategies.
Thus, the study has demonstrated that the existing monitoring network could be agreeably optimized using the presented statistical and geostatistical methods without the obvious fear of losing essential information from the current monitoring network. In order to accomplish sound groundwater resource management practices, necessary improvements in groundwater monitoring strategies is the basic key. As such, the efforts to optimize and evaluate the monitoring network will enhance the performance of the water management system. The methods presented in this research are convenient both in the case of inadequate networks, which need new monitoring wells as well as for dense monitoring networks, with redundant wells. The method can be conveniently used to find redundancy in the existing monitoring network as well as to locate necessary new monitoring wells in developing countries, whereas it is convenient to thin out the dense monitoring network (without losing valuable information) in developed countries. The integrated statistical and geostatistical method can be integrated with expert knowledge and field attributes to improve existing monitoring strategies by adding or removing existing monitoring wells in the monitoring network.

Acknowledgment

The presented work is funded by a Helmholtz PhD scholarship of the Helmholtz Centre for Environmental Research—UFZ in collaboration with the Martin Luther University Halle-Wittenberg, Germany. Groundwater monitoring data of Bitterfeld/Wolfen was provided by the LAF Sachsen-Anhalt. The data exchange was embedded in the SAFIRA I and II projects (funding number 02WT0023 & 02WT00412). This work has been possible because of kind co-operation and support from the Department of Hydrogeology and Environmental Geology, Institute of Geosciences and Geography, Martin Luther University, Halle (Saale), Germany.

Conflict of Interest

The authors declare no conflicts of interest.

References

  1. Wada, Y. Water Stress over the Year: Quantitative Analysis of Seasonality and Severity on a Global Scale. Master’s Thesis, Utrecht University, Utrecht, The Netherlands, 2008. [Google Scholar]
  2. Johnson, N.; Revenga, C.; Echeverria, J. Managing water for people and nature. Science 2001, 292, 1071–1072. [Google Scholar] [CrossRef] [PubMed]
  3. Thakur, J.K.; Thakur, R.K.; Ramanathan, A.; Kumar, M.; Singh, S.K. Arsenic contamination of groundwater in nepal—An overview. Water 2011, 3, 1–20. [Google Scholar] [CrossRef]
  4. Van Geer, F.C.; Bierkens, M.F.P.; Broers, H.P. Groundwater monitoring strategies. In Encyclopedia of Hydrological Sciences; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2006. [Google Scholar] [CrossRef]
  5. Erechtchoukova, M.G.; Chen, S.Y.; Khaiter, P.A. Application of optimization algorithms for the improvement of water quality monitoring systems. In Information Technologies in Environmental Engineering; Springer Berlin: Heidelberg, Germany, 2009; pp. 176–188. [Google Scholar]
  6. Thakur, J.K.; Gossel, W.; Weiß, H.; Wycisk, P. Long-Term Groundwater Monitoring Optimization: Improving Performance and Reducing Costs Associated for a 3d Regionalization of Groundwater Contamination—A Case Study of Bitterfeld/Wolfen, East Germany. In Proceedings of the Sustainable Water Resource Management and Treatment Technology, Nagpur, India, 19–22 January 2011; Labhasetwar, P., Ed.; NEERI: Nagpur, India, 2011; p. 20. [Google Scholar]
  7. Loaiciga, H.; Charbeneau, R.; Everett, L.; Fogg, G.; Hobbs, B.; Rouhani, S. Review of ground-water quality monitoring network design. J. Hydraul. Eng. 1992, 118, 11–37. [Google Scholar] [CrossRef]
  8. Thakur, J.K.; Gossel, W.; Weiß, H.; Wycisk, P. Optimizing long-term groundwater monitoring network using geostatistical methods. In Iahs Red Book: Models—Repositories of Knowledge; IAHS: London, UK, 2012; pp. 133–138. [Google Scholar]
  9. Massmann, J.; Freeze, R.A. Groundwater contamination from waste management sites: The interaction between risk-based engineering design and regulatory policy: 1. Methodology. Water Resour. Res. 1987, 23, 351–367. [Google Scholar] [CrossRef]
  10. Parker, J.; Kim, U.; Kitanidis, P.K.; Cardiff, M.; Liu, X. Stochastic cost optimization of multistrategy dnapl site remediation. Ground Water Monit. Remediat. 2010, 30, 65–78. [Google Scholar] [CrossRef]
  11. MacKenzie, M.C.; Palmer, R.N.; Millard, S.P. Analysis of statistical monitoring network design. J. Water Resour. Plan. Manag. 1987, 113, 599–615. [Google Scholar] [CrossRef]
  12. Ben-Jemaa, F.; Marino, M.A.; Loaiciga, H.A. Multivariate geostatistical design of ground-water monitoring networks. J. Water Resour. Plan. Manag. 1994, 120, 505–522. [Google Scholar] [CrossRef]
  13. Nabi, A.; Gallardo, A.; Ahmed, S. Optimization of a groundwater monitoring network for a sustainable development of the Maheshwaram Catchment, India. Sustainability 2011, 3, 396–409. [Google Scholar] [CrossRef]
  14. Thakur, J.; Diwakar, J.; Singh, S. Hydrogeochemical evaluation of groundwater of Bhaktapur Municipality, Nepal. Environ. Earth Sci. 2015. [Google Scholar] [CrossRef]
  15. Datta, B.; Chakrabarty, D.; Dhar, A. Optimal dynamic monitoring network design and identification of unknown groundwater pollution sources. Water Resour. Manag. 2009, 23, 2031–2049. [Google Scholar] [CrossRef]
  16. Chadalavada, S.; Datta, B.; Naidu, R. Uncertainty based optimal monitoring network design for a chlorinated hydrocarbon contaminated site. Environ. Monit. Assess. 2011, 173, 929–940. [Google Scholar] [CrossRef] [PubMed]
  17. Nunes, L.M.; Cunha, M.C.; Ribeiro, L. Groundwater monitoring network optimization with redundancy reduction. J. Water Resour. Plan. Manag. 2004, 130, 33–43. [Google Scholar] [CrossRef]
  18. Dhar, A.; Datta, B. Global optimal design of ground water monitoring network using embedded kriging. Ground Water 2009, 47, 806–815. [Google Scholar] [CrossRef] [PubMed]
  19. Meyer, P.D.; Brill, E.D. A method for locating wells in a groundwater monitoring network under conditions of uncertainty. Water Resour. Res. 1988, 24, 1277–1282. [Google Scholar] [CrossRef]
  20. Cieniawski, S.E.; Eheart, J.W.; Ranjithan, S. Using genetic algorithms to solve a multiobjective groundwater monitoring problem. Water Resour. Res. 1995, 31, 399–409. [Google Scholar] [CrossRef]
  21. Reed, P.M.; Minsker, B.S. Striking the balance: Long-term groundwater monitoring design for conflicting objectives. J. Water Resour. Plan. Manag. 2004, 130, 140–149. [Google Scholar] [CrossRef]
  22. Kollat, J.B.; Reed, P.M. Comparing state-of-the-art evolutionary multi-objective algorithms for long-term groundwater monitoring design. Adv. Water Resour. 2006, 29, 792–807. [Google Scholar] [CrossRef]
  23. Stollberg, R. Groundwater Contaminant Source Zone Identification at an Industrial and Abandoned Mining Site: A Forensic Backward-in-Time Modelling Approach; Martin Luther University: Halle (Saale), Germany, 2013. [Google Scholar]
  24. Heidrich, S.; Schirmer, M.; Weiss, H.; Wycisk, P.; Grossmann, J.; Kaschl, A. Regionally contaminated aquifers—Toxicological relevance and remediation options (Bitterfeld case study). Toxicology 2004, 205, 143–155. [Google Scholar] [CrossRef] [PubMed]
  25. Stollberg, R.; Gossel, W.; Wycisk, P.; Weiß, H. Source and Pathway Identification of Groundwater Contaminants Using a Backward Modeling Technique. In Proceedings of the 2nd International FEFLOW User Conference, Potsdam, Germany, 14–18 September; 2009. [Google Scholar]
  26. Simeonov, V.; Stratis, J.; Samara, C.; Zachariadis, G.; Voutsa, D.; Anthemidis, A.; Sofoniou, M.; Kouimtzis, T. Assessment of the surface water quality in northern greece. Water Res. 2003, 37, 4119–4124. [Google Scholar] [CrossRef]
  27. Kumar, M.; Ramanathan, A.; Keshari, A.K. Understanding the extent of interactions between groundwater and surface water through major ion chemistry and multivariate statistical techniques. Hydrol. Processes 2009, 23, 297–310. [Google Scholar] [CrossRef]
  28. Wold, S.; Esbensen, K.; Geladi, P. Principal component analysis. Chemom. Intell. Lab. Syst. 1987, 2, 37–52. [Google Scholar] [CrossRef]
  29. World Health Organization (WHO). Guidelines for Drinking-Water Quality, 4th ed.; WHO: Geneva, Swizerland, 2011. [Google Scholar]
  30. Environmental Protection Agency (EPA). National Primary Drinking Water Regulations. Available online: http://water.epa.gov/drink/contaminants/index.cfm (accessed on 25 December 2012).
  31. Weiss, H.; Schirmer, M.; Merkel, P. Development of novel reactive barrier technologies at the safira test site, Bitterfeld. In Advanced Groundwater Remediation: Active and Passive Technologies; Thomas Telford Publishing: London, UK, 2002; pp. 159–172. [Google Scholar]
  32. Haase, G.; Händel, L.; Nagel, C.; Opp, C.; Zierath, R. Environmental impacts by input of substances in landscapes of the district (Bezirk) of Leipzig, Germany. GeoJournal 1990, 22, 153–165. [Google Scholar] [CrossRef]
  33. Moisey, J.; Fisk, A.T.; Hobson, K.A.; Norstrom, R.J. Hexachlorocyclohexane (hch) isomers and chiral signatures of α-hch in the arctic marine food web of the northwater polynya. Environ. Sci. Technol. 2001, 35, 1920–1927. [Google Scholar] [CrossRef] [PubMed]
  34. Wycisk, P.; Stollberg, R.; Neumann, C.; Gossel, W.; Weiss, H.; Weber, R. Integrated methodology for assessing the hch groundwater pollution at the multi-source contaminated mega-site bitterfeld/wolfen. Environ. Sci. Pollut. Res. 2013, 20, 1907–1917. [Google Scholar] [CrossRef] [PubMed]
  35. Cameron, K.; Hunter, P.M. Geostatistical Temporal-Spatial (Gts) 1.0.0.; Steps, G., Ed.; MacStat Consulting, Ltd.: Colorado, CO, USA, 2010. [Google Scholar]
  36. Cameron, K.; Hunter, P. Using spatial models and kriging techniques to optimize long-term ground-water monitoring networks: A case study. Environmetrics 2002, 13, 629–656. [Google Scholar] [CrossRef]
  37. Isaaks, E.; Srivastava, R. Applied Geostatistics; Oxford University Press: New York, NY, USA, 1989; p. 561. [Google Scholar]
  38. Gilbert, R.O. Statistical Methods for Environmental Pollution Monitoring; Wiley: Hoboken, NJ, USA, 1987. [Google Scholar]
  39. Sen, P.K. Estimates of the regression coefficient based on Kendall’s tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  40. Brauner, J.S. Non-Parametric Estimation of Slope: Sen’s Method in Environmental Pollution. Environmental sampLing and Monitoring Primer. Available online: http://www.webapps.cee.vt.edu/ewr/environmental/teach/smprimer/sen/sen.html (accessed on 25 June 2012).
  41. Mann, H.B. Nonparametric tests against trend. Econom. J. Econom. Soc. 1945, 13, 245–259. [Google Scholar] [CrossRef]
  42. Kendall, M.G.; Stuart, A. The Advanced Theory of Statistics: Design and Analysis, and Time-Series; Macmillan Pub. Co.: London, UK, 1976; Volume 3. [Google Scholar]
  43. Software, G. Surfer 11. Available online: http://www.goldensoftware.com/products/surfer (accessed on 5 December 2011).
  44. Hubert, T. Methodische Ansätze zur Geologischen 3d-Modellierung Oberflächennaher Schichtenfolgen und Ihre Auswirkung auf die Numerische Grundwassermodellierung; Martin Luther Univesity: Halle, Germany, 2011. (In German) [Google Scholar]
  45. Thakur, J.K. Methods in Groundwater Monitoring: Strategies Based on Statistical, Geostatistical, and Hydrogeological Modelling and Visualization; Martin-Luther University: Halle-Wittenberg, Halle (Saale), Germany, 2013. [Google Scholar]
  46. Thakur, J.K.; Gossel, W.; Weiß, H.; Wycisk, P. Improving groundwater monitoring strategies based on groundwater monitoring network optimization. In TASK—The Centre of Competence for Soil, Groundwater and Site Revitalisation; Helmholtz Centre for Environmental Research - UFZ: Liepzig, Germany, 2011. [Google Scholar]
  47. Singh, S.K.; Singh, C.K.; Kumar, K.S.; Gupta, R.; Mukherjee, S. Spatial-temporal monitoring of groundwater using multivariate statistical techniques in Bareilly district of Uttar Pradesh, India. J. Hydrol. Hydromech. 2009, 57, 45–54. [Google Scholar] [CrossRef]
  48. Weiß, H.; Schirmer, M.; Teutsch, G.; Merkel, P. Sanierungsforschung in regional kontaminierten aquiferen (safira)—2. Projektüberblick und pilotanlage. Grundwasser 2002, 7, 135–139. (In German) [Google Scholar] [CrossRef]
  49. Weiß, H.; Teutsch, G.; Fritz, P.; Daus, B.; Dahmke, A.; Grathwohl, P.; Trabitzsch, R.; Feist, B.; Ruske, R.; Böhme, O. Sanierungsforschung in regional kontaminierten aquiferen (safira)—1. Information zum forschungsschwerpunkt am standort bitterfeld. Grundwasser 2001, 6, 113–122. (In German) [Google Scholar]
  50. Knotters, M.; Brus, D.; Oude Voshaar, J. A comparison of kriging, co-kriging and kriging combined with regression for spatial interpolation of horizon depth with censored observations. Geoderma 1995, 67, 227–246. [Google Scholar] [CrossRef]
  51. Bergström, S.; Graham, L.P. On the scale problem in hydrological modelling. J. Hydrol. 1998, 211, 253–265. [Google Scholar] [CrossRef]
  52. Dooge, J.C. Scale problems in hydrology. Spec. Publ. 1997, 48, 84–145. [Google Scholar]
  53. Gupta, V.K.; Rodríguez-Iturbe, I.; Wood, E.F. Runoff generation and basin response. In Scale Problem in Hydrology; Springer Netherlands: Rotterdam, Netherlands, 1986. [Google Scholar]
  54. Sivapalan, M.; Kalma, J.D. Scale problems in hydrology: Contributions of the robertson workshop. Hydrol. Process. 1995, 9, 243–250. [Google Scholar] [CrossRef]
  55. Sivapalan, M.; Grayson, R.; Woods, R. Scale and scaling in hydrology. Hydrol. Processes 2004, 18, 1369–1371. [Google Scholar] [CrossRef]
  56. Press, W.H.; Flannery, B.P.; Teukolsky, S.A.; Vetterling, W.T. Numerical Recipes in Fortran 77: Volume 1, Fortran Numerical Recipes: The Art of Scientific Computing; Cambridge University Press: Cambridge, UK, 1992; Volume 1. [Google Scholar]
  57. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes 3rd Edition: The Art of Scientific Computing; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  58. Storey, R.G.; Howard, K.W.; Williams, D.D. Factors controlling riffle-scale hyporheic exchange flows and their seasonal changes in a gaining stream: A three-dimensional groundwater flow model. Water Resour. Res. 2003. [Google Scholar] [CrossRef]
  59. Wroblicky, G.J.; Campana, M.E.; Valett, H.M.; Dahm, C.N. Seasonal variation in surface-subsurface water exchange and lateral hyporheic area of two stream-aquifer systems. Water Resour. Res. 1998, 34, 317–328. [Google Scholar] [CrossRef]
  60. Thakur, J.K.; Gossel, W.; Weiß, H.; Wycisk, P. Optimizing Long Term Groundwater Monitoring Network in a Complex Contaminated Site. In Proceedings of the HIGRADE Conference, Leipzig, Germany, 26 March 2012; Helmholtz Centre for Environmental Research (UFZ): Leipzig, Germany, 2012. [Google Scholar]

Share and Cite

MDPI and ACS Style

Thakur, J.K. Optimizing Groundwater Monitoring Networks Using Integrated Statistical and Geostatistical Approaches. Hydrology 2015, 2, 148-175. https://doi.org/10.3390/hydrology2030148

AMA Style

Thakur JK. Optimizing Groundwater Monitoring Networks Using Integrated Statistical and Geostatistical Approaches. Hydrology. 2015; 2(3):148-175. https://doi.org/10.3390/hydrology2030148

Chicago/Turabian Style

Thakur, Jay Krishna. 2015. "Optimizing Groundwater Monitoring Networks Using Integrated Statistical and Geostatistical Approaches" Hydrology 2, no. 3: 148-175. https://doi.org/10.3390/hydrology2030148

APA Style

Thakur, J. K. (2015). Optimizing Groundwater Monitoring Networks Using Integrated Statistical and Geostatistical Approaches. Hydrology, 2(3), 148-175. https://doi.org/10.3390/hydrology2030148

Article Metrics

Back to TopTop