Skip to main content

Main menu

  • Home
    • Journal home
    • Lyell Collection home
    • Geological Society home
  • Content
    • Online First
    • Current issue
    • Past issues
    • Collections
    • Supplementary publications
    • Open Access
  • Subscribe
    • GSL fellows
    • Institutions
    • Corporate
    • Other member types
  • Info for
    • Authors
    • Librarians
    • Readers
    • Access for GSL Fellows
    • Access for other member types
    • Press office
    • Help
  • Alerts
    • eTOC alerts
    • Online First alerts
    • RSS feeds
    • Newsletters
    • GSL blog
  • Submit
  • Geological Society of London Publications
    • Engineering Geology Special Publications
    • Geochemistry: Exploration, Environment, Analysis
    • Journal of Micropalaeontology
    • Journal of the Geological Society
    • Lyell Collection home
    • Memoirs
    • Petroleum Geology Conference Series
    • Petroleum Geoscience
    • Proceedings of the Yorkshire Geological Society
    • Quarterly Journal of Engineering Geology and Hydrogeology
    • Quarterly Journal of the Geological Society
    • Scottish Journal of Geology
    • Special Publications
    • Transactions of the Edinburgh Geological Society
    • Transactions of the Geological Society of Glasgow
    • Transactions of the Geological Society of London

User menu

  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Scottish Journal of Geology
  • Geological Society of London Publications
    • Engineering Geology Special Publications
    • Geochemistry: Exploration, Environment, Analysis
    • Journal of Micropalaeontology
    • Journal of the Geological Society
    • Lyell Collection home
    • Memoirs
    • Petroleum Geology Conference Series
    • Petroleum Geoscience
    • Proceedings of the Yorkshire Geological Society
    • Quarterly Journal of Engineering Geology and Hydrogeology
    • Quarterly Journal of the Geological Society
    • Scottish Journal of Geology
    • Special Publications
    • Transactions of the Edinburgh Geological Society
    • Transactions of the Geological Society of Glasgow
    • Transactions of the Geological Society of London
  • My alerts
  • Log in
  • My Cart
  • Follow gsl on Twitter
  • Visit gsl on Facebook
  • Visit gsl on Youtube
  • Visit gsl on Linkedin
Scottish Journal of Geology

Advanced search

  • Home
    • Journal home
    • Lyell Collection home
    • Geological Society home
  • Content
    • Online First
    • Current issue
    • Past issues
    • Collections
    • Supplementary publications
    • Open Access
  • Subscribe
    • GSL fellows
    • Institutions
    • Corporate
    • Other member types
  • Info for
    • Authors
    • Librarians
    • Readers
    • Access for GSL Fellows
    • Access for other member types
    • Press office
    • Help
  • Alerts
    • eTOC alerts
    • Online First alerts
    • RSS feeds
    • Newsletters
    • GSL blog
  • Submit

Use of non-parametric statistics as a tool for the hydraulic and hydrogeochemical characterization of hard rock aquifers

David Banks, Geir Morland and Bjørn Frengstad
Scottish Journal of Geology, 41, 69-79, 1 April 2005, https://doi.org/10.1144/sjg41010069
David Banks
Holymoor Consultancy, 8 Heaton Street, Brampton, Chesterfield, Derbyshire, S40 3AQ, UK (e-mail:david@holymoor.co.uk)
  • Find this author on Google Scholar
  • Search for this author on this site
Geir Morland
NorChip AS, Industriveien 8, N-3490 Klokkarstua, Norway
  • Find this author on Google Scholar
  • Search for this author on this site
Bjørn Frengstad
Norges Geologiske Undersøkelse, N-7491 Trondheim, Norway
  • Find this author on Google Scholar
  • Search for this author on this site
PreviousNext
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Synopsis

The new EU Water Framework Directive (Directive 2000/60/EC) implies that hard rock units will need to be considered as aquifers, and to be hydrogeologically and hydrogeochemically characterized. This is not a straightforward task, given that such aquifers are heterogeneous, discontinuous and strongly three-dimensional in the distributions of both hydraulic and hydrochemical properties. The Geological Survey of Norway has recently carried out characterizations of Norway’s hard rock aquifers using simple non-parametric statistical techniques, focusing on median and percentile values for aquifer properties and concentrations of chemical parameters. Such characterizations reveal far larger parameter ranges within given aquifer lithological units than between different units. Indeed, for most lithologies, median well yields and concentrations of major chemical parameters (e.g. pH, major ions) are remarkably similar. For certain elements, however, systematic differences can be detected. For example, elements such as uranium, radon and fluoride are significantly enriched in Precambrian granites and gneisses, and depleted in anorthosites.

Introduction

This paper considers some aspects of the philosophy of siting successful wells in hard rock aquifers. It also suggests means of characterizing such heterogeneous aquifers, both in terms of their hydraulic parameters and their hydrogeochemistry. In this paper, the term ‘hard rock’ aquifer is taken to mean the crystalline metamorphic or igneous rock (basalt, slate, granite, schist, gneiss etc.) underlying large areas of Scandinavia, Scotland and Wales. It is not taken to include deeply weathered basement terrain, with significant saprolite development, such as is typical of much of tropical India or Africa (Lloyd 1999).

Approaches to the science of locating successful water supply boreholes in hard rock aquifers typically focus on the identification of fracture zones or lineaments, by means of field surveys and study of maps, aerial photographs or satellite images. Geophysical techniques (ground penetrating radar, resistivity, electromagnetic and magnetic techniques) are often used to similar ends. Some authors have argued that certain fractures are likely to be more productive than others, depending on their orientation in past or current in-situ stress fields, resulting in Larsson's (1972) ‘Hydrotectonic’ theory (further discussed and criticized by Banks et al. 1994; 1996).

The authors do not wish to suggest that fracture zones never represent groundwater flow pathways, nor do they wish to suggest that geophysics are of little use or that stress orientation is not of significance. The relative lack of success of hydrogeologists in locating productive wells in hard rock terrain (and the continued apparent willingness of customers to use dowsers in preference to hydrogeologists) does, however, indicate that current techniques and theories lack predictive power. In particular, the following points should be noted:

  1. Not all fracture zones are water transmissive. For example, Banks et al. (1992) clearly demonstrated that water leakages to Hvaler subsea tunnel (southern Norway, Fig. 1) were derived from small fractures in relatively massive granite, rather than from major fracture zones and crush zones (which were often ‘tightened’ by smectite-rich fault gouge).

  2. Remote sensing and most geophysical techniques are unable to distinguish between water-transmissive fractures and clay/gouge-filled fractures (both result in lineaments in the terrain, and high conductivity electromagnetic responses).

  3. Most ‘hydrotectonic’ theories are simply too plausible: too many tectonic/stress scenarios can be envisaged to explain the yield (or lack thereof) of a set of boreholes. Such theories have explanatory power, but lack predictive power (Banks et al. 1996).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Water leakages in the Hvaler tunnel (after Banks et al. (1992), reproduced with permission of the Geological Society Publishing House). Note that zones of water leakage and injection grouting do not necessarily coincide with fracture zones as detected by geophysics and in the tunnel itself.

Many hydrogeologists approach well siting in hard rock terrain as a game of chess. By using the right tools and the right strategies, a positive outcome ought to be assured. In reality, this is often a delusion. The resources at the disposal of most groundwater supply investigations are insufficient to adequately characterize the complexity of most hard rock aquifers. The inability of hydrogeologists to predict the hydraulic properties of hard rock aquifers has its roots in the fractured, discontinuous and extremely heterogeneous nature of such aquifers. In a similar manner, the chemical quality of water from a single prospective borehole cannot be adequately predicted, and also has its roots in heterogeneous and discontinuous fracturing and fracture mineralization. As will be seen later in this article, hydrochemical parameters exhibit wide variations within a single lithology, largely reflecting heterogeneity in fracture residence times and hydrodynamics.

The authors would argue that well siting in hard rock terrain is more akin to a game of bingo. The outcome of drilling is largely determined by chance. Yet, hydrogeologists are still required to make prognoses for clients. Likewise, the EU Water Framework Directive (Directive 2000/60/EC) requires the hydrogeological and hydrochemical characterization of hard rock aquifers. It is clearly problematic to approach such characterization in a deterministic manner, but can we approach it in a probabilistic manner? Can we use statistics to assist in decision-making, when drilling for water? Can we statistically characterize hard rock aquifers?

The statistical characterization of well yield from hard rock aquifers

Some form of statistical technique has been used to characterize aquifers in most of the hard-rock-dominated nations of Europe, including Sweden, Finland, the Czech Republic and Norway. Indeed, the IAH hydrogeological map of Europe uses cumulative frequency diagrams for such purposes (e.g. Persson et al. 1985).

Let us take a data set from the borehole archives of the Geological Survey of Norway, for example, 24 boreholes from the island group of Hvaler, in the Precambrian Iddefjord Granite (Table 1). In this data set, yields range from 40001/h to dry (nominal 1 1/h entered into archive). We can calculate an arithmetic mean yield of 7641/h, with a standard deviation of 998 1/h, but is this a meaningful characterization of the Iddefjord Granite of Hvaler? Figure 2 suggests not: parameters such as the arithmetic mean and standard deviation only adequately characterize the central values of normally (or at least symmetrically) distributed data sets. Figure 2 demonstrates the distribution of yields to be highly skewed towards low yields (typical of all hard rock yield distributions). To characterize non-normally distributed data sets, we need to use ‘non-parametric’ statistics. The median value 400 1/h is taken as a central value, while the interquartile range (787 1/h) between the 25th percentile value (138 1/h) and the 75th percentile value (925 1/h) is a non-parametric measure of variance.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Histogram showing the frequency distribution of the borehole yield data set for 24 boreholes in the Iddefjord Granite of Hvaler (Table 1).

View this table:
  • View inline
  • View popup
TABLE 1

Yield (from short-term testing) of 24 water supply boreholes in the Iddefjord Granite of the Norwegian Hvaler Islands. The yield of 1 1/h is probably a nominal figure assigned to a ‘dry’ borehole.

In a sense, these statistics, either as raw numbers or as a cumulative frequency distribution (c.f.d.) curve (Fig. 3), represent a characterization of the Iddefjord Granite aquifer (provided that a great enough data set is available). Although neither boxplots nor c.f.d. diagrams make any specific assumptions about the type or nature of the distribution of yield, it is intriguing that, for many lithologies, the distribution of borehole yields often approximates to log-normality (Banks 1998). In such a case, the distribution can be characterized by only two parameters, the median (which will be identical to the geometric mean) and the standard deviation of the log values. The explanation for this phenomenon lies in the hierarchical evolution of fractures within the lithosphere and the observed lognormal distribution of fracture apertures (Banks et al. 1996).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Cumulative frequency distribution diagram for the borehole yield data set for 24 boreholes in the Iddefjord Granite of Hvaler (Table 1). The median is located at cumulative probability = 0.5, the 25th percentile at cumulative probability = 0.25 etc.

Furthermore, it can be argued that such statistics are meaningful in providing a first-order characterization of the hydraulic properties of an aquifer, given that well yield will be closely related to specific capacity (yield per metre drawdown) and that specific capacity is approximately proportional to aquifer transmissivity. Czech and Swedish hydrogeologists have used borehole yield as a surrogate measure for aquifer transmissivity when producing hydrogeological maps (Jetel & Krasny 1968; Krasny 1975; Carlsson & Carlstedt 1977), while Wright (1997) has used a plot of well yield against specific capacity (the ‘QSC diagram’) to characterize aquifers in Ireland. Banks (1992) has specifically examined the published relationships between well yield and transmissivity and concluded that:Embedded Imagewhere Ta - apparent transmissivity (m2/s) and F = specific capacity (m3/s/m drawdown)

Application of statistical characterizations to hydrogeological problem solving

Cumulative frequency curves can be constructed for each specific lithology in an area. Although yield distributions are often very similar, some systematic differences can be distinguished (Fig. 4). From such curves, one may offer a client the possibility of making an informed decision as to the chances of obtaining a given yield. For example, the chances of obtaining a yield of 12001/h in the Iddefjord Granite are only around 28%. However, the chances of obtaining a yield of 100 1/h in almost all lithologies is very good, often around 90%. Such low yields are in fact adequate for many domestic properties, summer cottages or small farms, and this statistic helps to explain the apparent relative success of dowsers, who are typically employed by clients with such properties. This statistic also explains why many Scandinavian drillers are able to offer a ‘yield guarantee’ to small clients, absorbing the cost of the small risk of failure into their pricing structure.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Cumulative frequency diagram showing yield distribution curves for Norwegian water boreholes in the Precambrian Iddefjord Granite, Cambro-Silurian metasediments of the Norwegian Caledonian terrain, and Precambrian gneisses. The '1001/h' guideline shows the approximate 10% yield for most lithologies (i.e. 90% of wells yield better than this figure).The 50% and 72% guidelines show the median yield and the 72% yield for the Iddefjord Granite. After Banks & Robins (2002) and based on data from Morland (1997). Reproduced with permission from the Geological Survey of Norway.

Furthermore, yield distributions, once characterized, can be combined to permit an estimation of the chances of obtaining a given total yield from a given number of multiple wells. Banks (1998) used a spreadsheet method of combining frequency curves in order to estimate yield distributions from multiple wells and points out that, for single wells, the 50% likely yield (F50) is given by Y5Q=Mm, where Mm is the median yield of single borehole.

For a non-interfering wellfield of n wells, the total 50% probability yield approaches: Y50(tota\)=nMa, where Ma is the arithmetic mean yield of a single well (n = large). For smaller wellfields of n wells: nMm < F50(total)<nMa.

Gustafson (2002) takes a similar approach and uses equal density sampling of cumulative frequency distribution curves to estimate the likely total yield of the best i boreholes out of a total drilled number b boreholes (assuming that (b — i) boreholes will not be commissioned; Fig. 5). Such approaches are not deterministic, but they do allow a client to make an informed cost- risk-benefit analysis of drilling. The economic aspects of such decision-making are clearly discussed further by Gustafson (2002).

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Cumulative plots of total predicted borehole yields from the i best wells out of b drilled. These plots are based on a Swedish data set of (rather good) boreholes, sited on the basis of Very Low Frequency (VLF) geophysical measurements (see Müllern (1980) for a description of the VLF technique), and with a median yield of 36001/h for single wells (note that the median yield of all Swedish hard rock wells, whether sited by VLF methods or not, is only some 6001/h (n = 59 000)). After Gustafson (2002) and reproduced with the permission of the Geological Survey of Norway.

The use of non-parametric statistical presentations also allows us to test hydrogeological hypotheses. For example, Morland (1997) examined the entire groundwater well data set held by the Geological Survey of Norway, to try to establish whether yield distributions differed significantly between differing hard rock lithologies. To this end, he used a presentation termed the ‘boxplot’ (Fig. 6). Here, the central ‘box’ contains the interquartile range of data, with a horizontal line at the median value. The parentheses on either side of the median represent a robust 95% confidence interval on the median. The ‘whiskers’ represent the extraquartile data, with outlying data (as defined by a statistical algorithm) represented by squares. Such boxplots are highly useful tools, as they do not make any assumptions about the nature of the data distribution, nor do they place especial emphasis on extreme low or high data values which may represent experimental error or ‘abnormal’ values (the emphasis is rather on the central mass of data).

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Boxplots illustrating the distribution of borehole yields within different lithological subgroups in Norway (after Morland 1997). The numbers on the left vertical axis refer to the lithological group (Table 2), those on the right give the total number of boreholes (#) in each subgroup. The 'box' of the boxplot contains the central interquartile range of data, with a line at the median value (and square parentheses giving the 95% confidence interval on the median). The 'whiskers' show the extraquartile data range, with squares showing extreme outlying data.

View this table:
  • View inline
  • View popup
TABLE 2

Lithological groupings used in connection with Norwegian hydrogeological and hydrochemical data, based on the legend to Sigmond's (1992) bedrock map of Norway at scale 1:3 000 000.

Morland (1997) found that, for a total of 12 757 hard rock boreholes, the median yield was 600 1/h (±17 1/h at a 95% confidence level), the median borehole depth was 56 m (±0.58 m) and the median normalized yield was 12.04 1/h/drilled m (±0.38 1/h/m). Intriguingly, a similar investigation of some 59 000 Swedish hard rock boreholes (Gustafson 2002; Fagerlind 1988) also arrived at an identical median of 6001/h (25th percentile = c. 2301/h, 75th percentile = c. 15001/h, arithmetic mean = 1643 1/h).

Morland (1997) divided the Norwegian data set into subsets according to lithology (Table 2; as defined in Sigmond's (1992) geological map), and produced the boxplots shown in Figure 6. From this diagram, lithologies whose median yield differs significantly from each other can be identified by the fact that the 95% confidence intervals do not overlap. For example, the median yield of lithology 57 (Permian volcanics in the Oslo Rift, 25001/h) is significantly higher than any other lithology, while that of 79 does not differ significantly (at a 95% confidence level) from that of lithology 87. The lowest median yield (lithology 72, 2501/h) is from Caledonian gabbros, diorites and ultramafic rocks. It will be noted that, while the spread of yields within any given lithology is large, the difference in median yield between lithologies is modest. This is presumed to reflect the similar geomechanical properties of many hard rock lithologies.

Morland (1997) also demonstrated an intriguing inverse relationship between median yield and median well depth (Fig. 7), presumably reflecting the fact that successful boreholes are terminated at shallower depths than poorly yielding boreholes, which are drilled deeper in the hope of encountering a water strike (see Ronka (1993) for a similar phenomenon in Finland).

Fig. 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 7.

Morland's (1997) data set of Norwegian bedrock boreholes, with median normalized yield (1/h per drilled metre) plotted against lithological subgroup (see Table 2) in descending order from left to right. Also shown are median total yield (1/h per borehole) and median borehole depth (m).

Morland (1997) then divided a section of Precambrian gneissic terrain in western-central Norway into regions of differing post-glacial isostatic uplift. The median normalized yield value for each subregion was calculated, and Figure 8 was produced, appearing to confirm a relationship first tentatively identified by Rohr-Torp (1994), between borehole yield and postglacial neotectonic activity. Henriksen (2003) subsequently attempted to continue this study further east, into the Swedish Precambrian terrain. He failed to find such a strong degree of correlation as Morland (1997), however, although he concludes that isostatic rebound might be able to explain around 9% of the observed variation in borehole yield along some profiles.

Fig. 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 8.

Correlation between well yield and total postglacial isostatic uplift; the diagram shows the median yield of boreholes in Norwegian Precambrian rocks per metre drilled depth as a function of annual land uplift. After Morland (1997); Banks & Robins (2002), and reproduced with the permission of the Geological Survey of Norway.

A word of caution

Naturally, the utility of the statistical techniques described above depends on the quality of the data set under consideration. In this context, the following caveats should be noted:

  1. There will often be a tendency towards overestimation of yield distributions, because drillers in most countries tend not to report dry or failed wells to the relevant authorities.

  2. The yields reported to well databases are often short-term yields measured over a period of hours. These yields may well not be representative of long-term sustainable yields (Van Tonder et al. 1998).

Hydrogeochemistry

In the same way that borehole yield varies considerably from borehole to borehole over distances of only a few metres within a given hard rock aquifer, depending on whether a transmissive fracture is encountered, so does the hydrogeochemistry. Different fractures contain groundwater of differing chemistry, depending, amongst other things, on the degree of opportunity for water rock interaction on fracture walls (surface area relative to water volume, residence time) and fracture mineralogies (which may differ from whole-rock mineralogy). However, it is possible to use the same statistical techniques as applied to well yield (above) to quantify the chances of obtaining water of a given quality within a given aquifer, and also to reveal statistically significant differences in the water quality distributions of differing lithological units.

For example, Banks et al. (1998b) systematically analysed water samples from some 1604 Norwegian bedrock groundwater boreholes. They observed:

  • that bedrock groundwater chemistry, like borehole yield, could not be predicted at a given site in a deterministic manner;

  • that bedrock groundwaters were hydrochemically more mature than both surface waters and water from superficial drift aquifers, e.g. they contained higher dissolved ionic contents, higher pH and alkalinity etc.;

  • as regards water type, that bedrock groundwaters were typically of calcium bicarbonate composition, although a significant number could be classed as high pH, sodium bicarbonate waters;

  • that a weak trend to a more ‘mature’ hydrochemistry could be discerned with increasing borehole depth (Frengstad & Banks 2003);

  • that within any given lithology there was typically a very wide range of concentrations of most parameters;

  • that the systematic differences between different lithologies were, for most parameters (including major ions and pH) rather low. In fact, variation in concentration between boreholes within a given lithology was often orders of magnitude greater than the variation in median values between different lithologies.

Figure 9 shows a series of boxplots documenting the distribution of groundwater pH amongst differing lithologies in Norway (see Table 2). The wide range in pH values should be noted, as should the fact that the median pH in almost all lithologies is around 8.2 to 8.3. Banks et al. (1998b, 2000) and Frengstad & Banks (2000) suggested pH increases with progressive carbonate and silicate weathering and that this median pH represents an interval of buffering during calcite saturation and precipitation from the water. Ca+++HC03 – =CaC03+H+

Fig. 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 9.

Boxplots displaying distribution of laboratory-determined pH values in Norwegian hard rock groundwater, sorted according to lithological groups (see Table 2). For comparison, boxplots are displayed on the extreme left showing: the entire data set for all hard rock groundwaters (‘Bedrock’, n - 1604), for groundwaters from superficial drift aquifers (‘Drift’, n = 72) and for n = 1 surface water works (‘Surface water’, n = 1). See Banks et al. (1998b, 2000).

The variation in pH within the various lithologies is ascribed largely to varying degrees of hydrochemical maturity (i.e. water-rock interaction) which would probably depend on factors such as residence time, fracture shape etc., which ultimately depend more on the geomechanical properties of the aquifer than on their detailed lithochemistry. As the chemistry of so many other major ions and trace elements depends ultimately on pH and the carbonate system (Frengstad et al. 2001), the concentrations of the majority of hydrochemical parameters in hard rock aquifers can be ascribed to:

  • hydrodynamic factors (residence time, fracture shape)

  • universal hydrochemical equilibria; specifically the calcium carbonate system

and not, to any great extent, to the detailed geochemical or mineralogical composition of the host rock.

A few elements do, however, exhibit a discernible lithological dependence. Prominent among these are uranium and radon. Uranium is highly soluble under oxidizing conditions, while radon is ultimately dependent on the uranium and radium content of the solid phases in contact with water. These parameters are thus dependent to a lesser extent on, for example, water pH and the degree of hydrochemical evolution of the water relative to the carbonate system, and more on the uranium content of the host rock aquifer. Indeed, uranium concentrations in bedrock groundwater in Norway vary across six orders of magnitude (and, along with several other elements, exhibit a distribution approaching log-normality; Frengstad et al. 2000). Uranium (and also radon and fluoride) concentrations can also be seen to be clearly elevated in granites (group 92, in Fig. 10) and depleted in groundwaters from anorthosites (group 93). The apparently elevated concentrations in group 72 are, on detailed inspection, also related to minor granitic inliers within that group (Frengstad et al. 2002). Uranium concentrations of up to 2 mg/1 have been recorded elsewhere in Norway (Morland et al. 1997), while over 14 mg/1 have been found in Finland (Asikainen & Kahlos 1979). Norwegian radon concentrations were found to range up to 31 900 Bq/1 (Banks et al. 2000).

Fig. 10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 10.

Boxplots displaying distribution of uranium concentrations (determined by ICP-MS, in jig/1) in Norwegian hard rock groundwater, sorted according to lithological groups (see Table 2). For comparison, the extreme right-hand boxplot shows the entire data set for all hard rock groundwaters (‘Bedrock’, n - 476). See Frengstad et al. (2000); Banks et al. (2000). The apparently elevated concentrations in group 72 are, on detailed inspection, also related to minor granitic inliers within that group (Frengstad et al. 2002). Data below the analytical detection limit are plotted at a value of half the detection limit for the purposes of graphical presentation.

Boxplots also enable the comparison of groundwater composition in similar aquifers from different climatic conditions. For example, Banks et al. (1998a) compared groundwater from the granitic Isles of Scilly, with groundwaters from the granitic Hvaler islands of Norway. Huge hydrochemical differences were noted (Fig. 11): the Scilly waters were far less hydrochemically mature (lower pH, alkalinity, silicon, Na/Cl ratio, calcium content etc.), lower in fluoride, uranium and radon, but much richer in nitrate and potassium. The Hvaler waters appeared to be more chemically reducing in nature, with higher iron and manganese concentrations. The differences were ascribed to the following factors:

  • the Scilly Isles are currently submergent in nature, while the Hvaler Isles are emergent (post-glacial isostatic rebound);

  • Scilly thus has a long subaerial weathering history, while Hvaler has recently been glacially scoured;

  • the Hvaler Isles are thus composed of relatively fresh, unweathered granite; Scilly's granite is likely to be more weathered and possibly depleted in uranium and basic minerals along the main groundwater flow pathways;

  • the Scilly Isles have higher topographic gradients than Hvaler (faster groundwater flow);

  • the Scilly granite may be more permeable than that of Hvaler;

  • Scilly is subject to relatively intense floriculture and agriculture, compared to Hvaler; soils are sandy and unfrozen in winter and may be more prone to nitrate leaching.

Fig. 11.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 11.

Boxplots comparing distribution of selected parameters in granitic groundwater from the Hvaler Islands (n = 11) with granitic groundwater from the Isles of Scilly (n = 10). Based on data published by Banks et al. (1998a). Data below the analytical detection limit are plotted at a value of half the detection limit for the purposes of graphical presentation.

Shand & Frengstad (2002) and Frengstad (2002) have continued this study further, comparing a range of Norwegian hard rock lithologies with Welsh Lower Palaeozoic argillite aquifers, the Cornish Permian Carnmenellis Granite and the Scottish Devonian Sandstone aquifers. These studies confirmed the generally higher pH, higher alkalinity nature of the Norwegian groundwaters and their generally lower nitrate concentrations.

Conclusions

Hydrogeological prognosis

How might we restore the self respect of the hydrogeologist in hard rock terrain vis-à-vis the dowser? Should we continue to play our games of hydrogeological chess against granitic opponents, marshalling our resources of geophysics and structural analysis upon a playing board of aerial photos and LandSat images?

Alternatively, should we simply accept that hard rock hydrogeology is a game of bingo, and focus our energies on providing clients with realistic estimates of their chances of obtaining a well of adequate yield and chemistry, via statistical techniques such as those described above?

The authors would argue that the correct analogy is possibly a game of poker. The outcome is largely dependent on chance. However, the astute poker player will employ sensible strategies (i.e. using geophysics or lineament analysis, where appropriate) and will have an excellent understanding of the probabilities of various outcomes (using statistical analysis). The best poker players, not being content to accept the probabilities, will probably also have an ace up their sleeves, with which to cheat Mother Nature. One example of such a hydrogeological ‘ace’ would be the use of explosives or hydraulic fracturing (hydrofraccing). The latter is a technique of creating or opening subsurface fractures by injection of water at pressures of several MPa. During the past years, hydrofraccing techniques have advanced considerably (Less & Anderson 1994; Lloyd 1999), although even three decades ago, studies performed in Sweden suggested that hydrofraccing typically results in a median increase in bedrock borehole yield of around 2001/h, with increases of up to 18001/h (Müllern & Eriksson 1977; Fagerlind 1979; 1982).

Aquifer vulnerability and characterization

The EU Water Framework Directive requires the characterization of hard rock aquifers. The Nordic nations have progressed a long way on the road to this objective by the publication of groundwater atlases and maps (Lahermo et al. 1990; Aastrup et al. 1995) or by the analysis of large databases of hydraulic and hydrochemical information (Morland 1997; Banks et al. 1998b; Frengstad et al. 2000; 2001). The hard rock areas of the UK (including Scotland, Wales, SW England and parts of Northern Ireland) need to complete this task within a relatively short time frame. Currently, the British Geological Survey is contracted to characterize aquifer units to produce ‘aquifer property’ maps of Scotland. This current characterization involves placing each aquifer unit in one of several ‘classes’ based on its productivity (typical long-term well yield), groundwater chemistry (pH, mineralization, redox condition), depth of flow, length of flow pathway, transmissivity, porosity and degree of confinement (Fitzsimons et al. 2004).

Additionally, the vulnerability of Scottish groundwater is being assessed in a GIS environment, using a methodology described by O'Dochartaigh et al. (2005). This methodology considers, however, only the vertical contaminant pathway from the surface to the water table. The vulnerability of this pathway, according to O'Dochartaigh et al.’s methodology, does not depend on the properties or lithology of fractured bedrock but is largely determined by the nature of the soils and superficial deposits. However, once any contamination reaches the water table in a fractured bedrock aquifer, it may migrate laterally down-gradient towards specific receptors such as wells and springs. In this case, one would expect the vulnerability of the receptor to depend to some extent on the hydraulic properties of the fractured aquifer.

It is argued that the analysis of available well yield and hydrochemical data sets using simple non-parametric statistical methods would be a straightforward and inexpensive way of characterizing the hydraulic properties and hydrochemistry of fractured, hard rock aquifers. Authors such as Shand & Frengstad (2002) have already attempted the statistical characterization of the hydrochemistry of some British hard rock aquifers. Available data from Scotland are sparse and possibly somewhat non-representative, however, and those available are presented in publications by Robins (1999) and Shand & Frengstad (2002). It is likely that locally operating drillers may hold additional information on the yields and water chemistry of smaller private and farm wells in the hard rock areas of Scotland. There is a case for the systematic collation of such data, via the development of mutually trusting and beneficial relationships between drillers and regulators, in order to establish larger data sets to support statistical analysis. If this could be achieved, it can be argued that the characterization of a Scottish granite aquifer is as simple as presenting a data table such as Table 3.

View this table:
  • View inline
  • View popup
TABLE 3

Example of the possible characterization of, for example, a granite aquifer using non-parametric statistical analysis - Norwegian Precambrian Granite (Lithological Group 92, see Table 2) (data derived from studies of Morland 1997 and Banks et al. 1998b).

  • © 2005 Scottish Journal of Geology

References

  1. ↵
    1. AASTRUP, M.,
    2. THUNHOLM, B.,
    3. JOHNSON, J.,
    4. BERTILLS, U. &
    5. BERNTELL, A.
    1995. Groundwater chemistry in Sweden. Report 4416, Naturvårdsverket Förlag.
  2. ↵
    1. ASIKAINEN, M. &
    2. KAHLOS, H.
    1979. Anomalously high concentrations of uranium, radium and radon in water from drilled wells in the Helsinki region. Geochimica et Cosmochimica Acta, 43, 1681–1686.
    OpenUrlCrossRefWeb of Science
  3. ↵
    1. BANKS, D.
    1992. Estimation of apparent transmissivity from capacity-testing of boreholes in bedrock aquifers. Applied Hydrogeology. 4/92, 5–19.
    OpenUrl
  4. ↵
    1. BANKS, D.,
    2. SOLBJØRG, M.-L. &
    3. ROHR-TORP, E.
    1992. Permeability of fracture zones in a Precambrian granite. Quarterly Journal of Engineering Geology, 25, 377–388.
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. BANKS, D.,
    2. ROHR-TORP, E. &
    3. SKARPHAGEN, H.
    1994. Groundwater resources in hard rock; experiences from the Hvaler study, Southeastern Norway. Applied Hydrogeology, 2/94, 33–42.
    OpenUrl
  6. ↵
    1. BANKS, D.,
    2. ODLING, N.E.,
    3. SKARPHAGEN, H. &
    4. ROHR-TORP, E.
    1996. Permeability and stress in crystalline rocks. Terra Nova, 8, 223–235.
    OpenUrlCrossRefWeb of Science
  7. ↵
    1. BANKS, D.
    1998. Predicting the probability distribution of yield from multiple boreholes in crystalline bedrock. Ground Water, 36, 269–274.
    OpenUrlCrossRefWeb of Science
  8. ↵
    1. BANKS, D. &
    2. ROBINS, N.
    2002. An Introduction to Groundwater in Crystalline Bedrock, Norges geologiske undersøkelse, Trondheim, Norway.
  9. ↵
    1. BANKS, D.,
    2. REIMANN, C. &
    3. SKARPHAGEN, H.
    1998a. The comparative hydrochemistry of two granitic island aquifers: the Isles of Scilly, U.K. and the Hvaler Islands, Norway. The Science of the Total Environment, 209, 169–183.
    OpenUrlCrossRefPubMed
  10. ↵
    1. BANKS, D.,
    2. FRENGSTAD, B.,
    3. MIDTGÅRD, A.K.,
    4. KROG, J.R. &
    5. STRAND, T.
    1998b. The chemistry of Norwegian groundwaters: I. the distribution of radon, major and minor elements in 1604 crystalline bedrock groundwaters. The Science of the Total Environment, 222, 71–91.
    OpenUrlCrossRefPubMed
  11. ↵
    1. BANKS, D.,
    2. FRENGSTAD, B.,
    3. SKREDE, A.K.,
    4. KROG, J.R.,
    5. STRAND, T.,
    6. SIEWERS, U. &
    7. LIND, B.
    2000. Groundwater – not just water (in Norwegian) Grunnvann – ikke bare vann, Gråsteinen, 6, Trondheim.
  12. ↵
    1. CARLSSON, L. &
    2. CARLSTEDT, A.
    1977. Estimation of transmissivity and permeability in Swedish bedrock. Nordic Hydrology, 8, 103–116.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. FAGERLIND, T.
    1979. Högtryckspumpning, sprängning och gradborrning. Information från Brunnsarkivet (Sveriges geologiska undersökning, Uppsala), 2/1979, 4–8. Sveriges geologiska undersökning, Uppsala.
    1. FAGERLIND, T.
    1982. Tryckning av brunnar. Information från Brunnsarkivet (Sveriges geologiska undersökning, Uppsala), 2/1982, 16–17.
    OpenUrl
  14. ↵
    1. FAGERLIND, T.
    1988. Urberget – en undervärderad grundvattenakvifer? Grundvattenförsörjning – Geohydrologi i praktiken; Symposium, VIAK AB and Swedish Hydrological Council.
  15. ↵
    1. FITZSIMONS, V.,
    2. CLEWS, J. &
    3. SIMPSON, E.
    2004. Environmental Characterisation of Scotland’s Water Environment: Supporting Information for Groundwaters. Report, 29/6/04, Scottish Environment Protection Agency.
  16. ↵
    1. FRENGSTAD, B. &
    2. BANKS, D.
    2000. Evolution of high-pH Na-HCO3 groundwaters in anorthosites: silicate weathering or cation exchange? In Sililo, et al. (ed.) Groundwater: Past Achievements and Future Challenges. Proceedings of XXXIInd Congress of the International Association of Hydrogeologists, Cape Town, South Africa, Balkema, Rotterdam, 493–498.
  17. ↵
    1. FRENGSTAD, B.,
    2. MIDTGÅRD SKREDE, A.K.,
    3. BANKS, D.,
    4. KROG, J.R. &
    5. SIEWERS, U.
    2000. The chemistry of Norwegian groundwaters: III. The distribution of trace elements in 476 crystalline bedrock groundwaters, as analysed by ICP-MS techniques. The Science of the Total Environment, 246, 21–40.
    OpenUrlCrossRefPubMed
  18. ↵
    1. FRENGSTAD, B.,
    2. BANKS, D. &
    3. SIEWERS, U.
    2001. The chemistry of Norwegian groundwaters: IV. The pH-dependence of element concentrations in crystalline bed-rock groundwaters. The Science of the Total Environment, 277, 101–117.
    OpenUrlCrossRefPubMed
  19. ↵
    1. FRENGSTAD, B.
    2002. Groundwater quality of crystalline bedrock aquifers in Norway. Dr. Ing. Thesis, 2002.53 Norwegian University of Science and Technology.
  20. ↵
    1. FRENGSTAD, B. &
    2. BANKS, D.
    2003. Groundwater chemistry related to depth of shallow crystalline bedrock boreholes in Norway. In Krásný, J., et al. (ed.) Proceedings of the International Association of Hydrogeologists’ Conference on Groundwater in Fractured Rocks, 15th–19th September 2003, Prague, Czech Republic. Extended abstracts / IHP-VI, Series on groundwater, 7, 203–204.
  21. ↵
    1. FRENGSTAD, B.,
    2. BANKS, D.,
    3. MIDTGARD SKREDE, A.K.,
    4. KROG, J.R.,
    5. SIEWERS, U. &
    6. STRAND, T.
    2002. The hydrochemistry of crystalline bedrock groundwater in Norway. Norges Geologiske Undersokelse Bulletin, 439, 87–98.
    OpenUrl
  22. ↵
    1. GUSTAFSON, G.
    2002. Strategies for groundwater prospecting in hard rocks: probabilistic approach. Norges Geologiske Undersøkelse Bulletin, 439, 21–25.
    OpenUrl
  23. ↵
    1. HENRIKSEN, H.
    2003. The role of some regional factors in the assessment of well yields from hard-rock aquifers of Fennoscandia. Hydrogeology Journal, 11, 628–645.
    OpenUrlCrossRef
    1. JETEL, J. &
    2. KRÁSNY, J.
    1968. Approximate aquifer characteristics in regional hydrogeological study. Ustredniho ustavu geologickeho (Prague), 43, 459–461.
    OpenUrl
    1. KRÁSNY, J.
    1975. Variation in transmissivity of crystalline rocks in southern Bohemia. Ustredniho ustavu geologickeho (Prague), 50, 207–216.
    OpenUrl
  24. ↵
    1. LAHERMO, P.,
    2. ILMASTI, M.,
    3. JUNTUNEN, R. &
    4. TAKA, M.
    1990. The geochemical atlas of Finland, Part 1: The hydrogeochemical mapping of Finnish groundwater, Geological Survey of Finland, Espoo.
    1. LARSSON, I.
    1972. Groundwater in granite rocks and tectonic models. Nordic Hydrology, 3, 111–129.
    OpenUrlAbstract/FREE Full Text
    1. LESS, C. &
    2. ANDERSEN, N.
    1994. Hydrofracture – state of the art in South Africa. Applied Hydrogeology, 2 (2), 59–63.
    OpenUrlCrossRef
  25. ↵
    1. LLOYD, J.W.
    1999. Water resources in hard rock aquifers in arid and semi-arid areas. Studies and Reports in Hydrogeology, 58, UNESCO.
  26. ↵
    1. MORLAND, G.
    1997. Petrology, lithology, bedrock structures, glaciation and sea level. Important factors for groundwater yield and composition of Norwegian bedrock boreholes (2 vols). Report, 97.122, Norges geologiske undersøkelse.
  27. ↵
    1. MORLAND, G.,
    2. REIMANN, C.,
    3. STRAND, T.,
    4. SKARPHAGEN, H.,
    5. BANKS, D.,
    6. BJORVATN, K.,
    7. HALL, G.M. &
    8. SIEWERS, U.
    1997. The hydrogeochemistry of Norwegian bedrock groundwater – selected parameters (pH, F-, Rn, U, Th, B, Na, Ca) in samples from Vestfold and Hordaland, Norway. Norges geologiske undersøkelse Bulletin, 432, 103–118.
    OpenUrl
    1. MÜLLERN, C.F. &
    2. ERIKSSON, A.
    1977. Kapacitetsökning hos bergborrade brunnar genom högtrckspumpar. Vannet i Norden, 2.
    1. MÜLLERN, C.F.
    1980. Airborne geophysical measurements used for hydrogeological mapping. Proceedings of 6th Nordic Hydrological Conference, Vemdalen (10th–16th August 1980), 135–142.
    1. O’DOCHARTAIGH, B.E.,
    2. BALL, D.F.,
    3. MACDONALD, A.M.,
    4. LILLY, A.,
    5. FITZSIMONS, V.,
    6. DEL RIO, M. &
    7. AUTON, C.A.
    2005. Mapping groundwater vulnerability in Scotland: a new approach for the Water Framework Directive. Scottish Journal of Geology, 41, 21–30.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. PERSSON, G.,
    2. DEGEER, J.,
    3. EGOROV, S.V.,
    4. KIRKHUSMO, L.A.,
    5. OLSSON, T.,
    6. PANOVA, Z.P.,
    7. TAKA, M. &
    8. WIKNER, T.
    1985. Explanatory notes for the hydrogeological map of Europe, scale 1:1,500,000. Sheet D2 Haparanda. IAH/Commission for the Geological Map of the World, Bundesanstalt für Geowissenschaften und Rohstoffe and UNESCO, Hannover and Paris.
  29. ↵
    1. ROBINS, N.S.
    1999. Groundwater occurrence in the Lower Palaeozoic and Precambrian rocks of the UK: implications for source protection. Journal of the Chartered Institution of Water and Environmental Management, 13, 447–453.
    OpenUrlCrossRefWeb of Science
  30. ↵
    1. ROHR-TORP, E.
    1994. Present uplift rates and groundwater potential in Norwegian hard rocks. Norges geologiske undersøkelse Bulletin, 426, 47–52.
    OpenUrl
    1. RÖNKA, E.
    1993. Increased depth of drilled wells – advantage or drawback?. In Banks, S.B. & Banks, D. (eds.) Hydrogeology of Hard Rocks. Memoir of the 24th Congress of the International Association of Hydrogeologists. Memoir, Norges Geologiske Undersøkelse, Trondheim, 772–779.
    1. SHAND, P. &
    2. FRENGSTAD, B.
    2002. Baseline groundwater quality. A comparison of Selected British and Norwegian Aquifers. Groundwater Systems and Water Quality Programme Internal Report, IR/01/177, British Geological Survey.
    1. SIGMOND, E.M.O.
    1992. Bedrock map of Norway and adjacent ocean areas. Scale 1:3 million, Norges geologiske undersøkelse, Trondheim, Norway.
  31. ↵
    1. VAN TONDER, G.,
    2. KUNSTMANN, H. &
    3. XU, Y.
    1998. Estimation of the sustainable yield of a borehole including boundary information, drawdown derivatives and uncertainty propagation, Institute for Groundwater Studies, Bloemfontein, RSA.
  32. ↵
    1. WRIGHT, G.
    1997. QSC graphs – a tool to assist in aquifer classification. GSI (Geological Survey of Ireland) Groundwater Newsletter, 32, 11–20.
    OpenUrl
PreviousNext
Back to top

In this issue

Scottish Journal of Geology
Volume 41, Issue 1
April 2005
  • Table of Contents
  • Table of Contents (PDF)
  • Index by author
  • Back Matter (PDF)
  • Front Matter (PDF)
Alerts
Sign In to Email Alerts with your Email Address
Citation tools

Use of non-parametric statistics as a tool for the hydraulic and hydrogeochemical characterization of hard rock aquifers

David Banks, Geir Morland and Bjørn Frengstad
Scottish Journal of Geology, 41, 69-79, 1 April 2005, https://doi.org/10.1144/sjg41010069
David Banks
Holymoor Consultancy, 8 Heaton Street, Brampton, Chesterfield, Derbyshire, S40 3AQ, UK (e-mail:david@holymoor.co.uk)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Geir Morland
NorChip AS, Industriveien 8, N-3490 Klokkarstua, Norway
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bjørn Frengstad
Norges Geologiske Undersøkelse, N-7491 Trondheim, Norway
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Permissions
View PDF
Share

Use of non-parametric statistics as a tool for the hydraulic and hydrogeochemical characterization of hard rock aquifers

David Banks, Geir Morland and Bjørn Frengstad
Scottish Journal of Geology, 41, 69-79, 1 April 2005, https://doi.org/10.1144/sjg41010069
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Email to

Thank you for sharing this Scottish Journal of Geology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Use of non-parametric statistics as a tool for the hydraulic and hydrogeochemical characterization of hard rock aquifers
(Your Name) has forwarded a page to you from Scottish Journal of Geology
(Your Name) thought you would be interested in this article in Scottish Journal of Geology.
Print
Download PPT
  • Tweet Widget
  • Facebook Like
  • Google Plus One
  • Article
    • Synopsis
    • Introduction
    • The statistical characterization of well yield from hard rock aquifers
    • Application of statistical characterizations to hydrogeological problem solving
    • A word of caution
    • Hydrogeochemistry
    • Conclusions
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

Related Articles

Similar Articles

Cited By...

More in this TOC Section

  • Numerical testing of conceptual models of groundwater flow: a case study using the Dumfries Basin aquifer
  • Agriculture and diffuse pollution: groundwater nitrate vulnerable zones in Scotland
Show more: THEMATIC SET: Hydrogeology of Scotland
  • Most read
  • Most cited
Loading
  • A sauropod-dominated tracksite from Rubha nam Brathairean (Brothers’ Point), Isle of Skye, Scotland
  • Discussion of ‘The Highland Boundary Fault and the Highland Border Complex’ by B. J. Bluck, Scottish Journal of Geology, 46, 113–124
  • Reply to the discussion by Tanner on ‘The Highland Boundary Fault and the Highland Border Complex’ Bluck (2010) Scottish Journal of Geology, 46, 113–124
  • Book Review
  • Introducing Geomorphology: A Guide to Landforms and Processes
More...

Scottish Journal of Geology

  • About the journal
  • Editorial Board
  • Submit a manuscript
  • Author information
  • Supplementary Publications
  • Subscribe
  • Pay per view
  • Alerts & RSS
  • Copyright & Permissions
  • Activate Online Subscription
  • Feedback
  • Help

Lyell Collection

  • About the Lyell Collection
  • Lyell Collection homepage
  • Collections
  • Open Access Collection
  • Open Access Policy
  • Lyell Collection access help
  • Recommend to your Library
  • Lyell Collection Sponsors
  • MARC records
  • Digital preservation
  • Developing countries
  • Geofacets
  • Manage your account
  • Cookies

The Geological Society

  • About the Society
  • Join the Society
  • Benefits for Members
  • Online Bookshop
  • Publishing policies
  • Awards, Grants & Bursaries
  • Education & Careers
  • Events
  • Geoscientist Online
  • Library & Information Services
  • Policy & Media
  • Society blog
  • Contact the Society

 

EGS logo

 

 

Published by The Geological Society of London, registered charity number 210161

Print ISSN 
0036-9276
Online ISSN 
2041-4951

Copyright © 2019 Scottish Journal of Geology