Journal cover Journal topic
Hydrology and Earth System Sciences An interactive open-access journal of the European Geosciences Union
Journal topic
Hydrol. Earth Syst. Sci., 23, 3603–3629, 2019
https://doi.org/10.5194/hess-23-3603-2019
Hydrol. Earth Syst. Sci., 23, 3603–3629, 2019
https://doi.org/10.5194/hess-23-3603-2019

Review article 04 Sep 2019

Review article | 04 Sep 2019

# Error in hydraulic head and gradient time-series measurements: a quantitative appraisal

Error in hydraulic head and gradient time-series measurements: a quantitative appraisal
Gabriel C. Rau1,2, Vincent E. A. Post3, Margaret Shanafield4, Torsten Krekeler3, Eddie W. Banks4, and Philipp Blum1 Gabriel C. Rau et al.
• 1Karlsruhe Institute of Technology (KIT), Institute of Applied Geosciences (AGW), Karlsruhe, Germany
• 2Connected Waters Initiative Research Centre (CWI), the University of New South Wales (UNSW), Sydney, Australia
• 3Federal Institute for Geosciences and Natural Resources (BGR), Hanover, Germany
• 4National Centre for Groundwater Research and Training (NCGRT), College of Science and Engineering, Flinders University, Adelaide, Australia

Correspondence: Gabriel C. Rau (gabriel.rau@kit.edu)

Abstract

Hydraulic head and gradient measurements underpin practically all investigations in hydrogeology. There is sufficient information in the literature to suggest that head measurement errors can impede the reliable detection of flow directions and significantly increase the uncertainty of groundwater flow rate calculations. Yet educational textbooks contain limited content regarding measurement techniques, and studies rarely report on measurement errors. The objective of our study is to review currently accepted standard operating procedures in hydrological research and to determine the smallest head gradients that can be resolved. To this aim, we first systematically investigate the systematic and random measurement errors involved in collecting time-series information on hydraulic head at a given location: (1) geospatial position, (2) point of head, (3) depth to water, and (4) water level time series. Then, by propagating the random errors, we find that with current standard practice, horizontal head gradients $<{\mathrm{10}}^{-\mathrm{4}}$ are resolvable at distances 170 m. Further, it takes extraordinary effort to measure hydraulic head gradients $<{\mathrm{10}}^{-\mathrm{3}}$ over distances <10 m. In reality, accuracy will be worse than our theoretical estimates because of the many possible systematic errors. Regional flow on a scale of kilometres or more can be inferred with current best-practice methods, but processes such as vertical flow within an aquifer cannot be determined until more accurate and precise measurement methods are developed. Finally, we offer a concise set of recommendations for water level, hydraulic head and gradient time-series measurements. We anticipate that our work contributes to progressing the quality of head time-series data in the hydrogeological sciences and provides a starting point for the development of universal measurement protocols for water level data collection.

1 Introduction

Water level and hydraulic head time series are critical for understanding water flow-related processes and properties in both surface and subsurface aquatic environments. At the surface, water levels are important for understanding relationships between water level and flow and for estimating surface water–groundwater interactions (e.g., Kalbus et al.2006; McCallum et al.2014). In the subsurface, measurements of hydraulic head are used to determine groundwater flow, estimate aquifer properties, and investigate aquifer processes such as the response to pumping or groundwater recharge . While it has been confirmed in several studies that the accuracy of water level measurements is a limiting factor for drawing conclusions about hydrogeological processes , measurement errors are not always properly recognised .

Pressure transducers (PTs) have been used since the 1960s to measure water level , and the development and availability of a wide variety of commercial instruments has made collection of high-temporal-resolution water level time series common practice. This has been a major advancement in our capability to study hydrological processes, but the proper use of automated sensors means that researchers need to have a good understanding of instrument technology and operating procedures. This is by no means trivial and certainly much more complex than collecting manual measurements. Knowledge is already required during the procurement phase, as there are many brands and logger types available, and the specific research objectives of a project determine which sensors are suitable and which are not . The same is true for modern positioning and levelling instruments, needed to establish the horizontal and vertical position of the monitoring well (e.g.,  Brinker1995; Hegarty2017). The storage and quality assurance of the large volume of time-series data are not straightforward either and can require programming skills to process data in an efficient manner. All things considered, modern hydrologists and hydrogeologists require a broad skill set, which is typically too extensive to be comprehensively covered in standard textbooks and water-related educational programmes.

Yet water level measurement lies at the heart of all hydrogeological investigation, and knowledge of the measurement error associated with modern instruments is fundamental to the collection of reliable time-series data. Studies on the topic published in the literature focus mainly on the instruments themselves. One of the first estimates of PT drift was published by , who showed how these errors would have led to incorrect interpretation of water levels at several sites. More recently, examined the accuracy and drift of different brands of PTs and found that the manufacturers' specifications were not met during field deployment. The effect of temperature on sensor performance has also received some attention . These studies concluded that strong temperature fluctuations such as those that occur under field conditions affect PTs of all types.

More comprehensive treatments of the subject tend to be published as reports by national research organisations. Prime examples include and , who not only discussed sensor technology but also provided technical procedures for collecting water levels and some of the errors involved. Moreover, some relevant works were published in a non-English language or as conference proceedings so that, despite their usefulness, their findings did not permeate the indexed international literature.

When collecting hydraulic head time series in the field, many different factors apart from instrument drift influence the stability of the measurement set-up . These include cable stretch, well clogging, sensor fouling, variable-density effects, and even changes in the vertical position of the observation well. This requires regular field site maintenance, recalibration and record-keeping. However, without knowledge of the magnitude of the water level error caused by such effects, there is no general guidance to develop adequate systematic field procedures. Unrecognised and unaccounted for systematic errors can accumulate (or cancel), leading to unquantifiable inaccuracies, while the random errors increase the uncertainty. contended that the propagation of measurement errors can result in ±100 % uncertainty in calculated flow velocities and that the uncertainty of the head gradient may be of a similar magnitude as that of the hydraulic conductivity.

We anticipate that our analysis is helpful to field practitioners at all levels and can be used as an educational resource. By providing a concise list of best practice recommendations at the end of the paper, we intend to provide a starting point for the development of comprehensive and universal international standard procedures, which are currently lacking.

2 Review of measurements and error terminology

Figure 1Overview of the four individual measurements (enumerated as steps and marked in red) required to calculate time series of hydraulic head (one location) and gradient (two locations) using two different types of groundwater monitoring infrastructure (GMI). Location 1 shows a cased borehole that is open to the atmosphere (open GMI), whereas Location 2 illustrates a fully grouted-in piezometer (closed GMI). The boreholes are drawn at an angle to highlight the importance of errors caused by inclination during construction of the borehole.

## 2.1 From measurements to heads

In this work we use the term groundwater-monitoring infrastructure (GMI) as an umbrella term for open and cased boreholes, wells, and standpipe or grouted-in piezometers (Sect. 4). The most typical GMI in hydrogeology consists of boreholes equipped with a standpipe piezometer, where the standing water level is in contact with the atmosphere (open GMI; see Location 1 in Fig. 1) and therefore readily accessible for measurement. Fully grouted-in piezometers contain a single or multi-array string of PTs and are closed to the atmosphere (closed GMI; see Location 2 in Fig. 1) and are often used in mining and geotechnical engineering (e.g., McKenna1995; Mikkelsen and Green2003).

GMI allows access to measuring depth to water or groundwater pressure from which the hydraulic head can be calculated (terminology is illustrated in Fig. 1). The hydraulic head is defined as (e.g., Hubbert1940; Freeze and Cherry1979)

$\begin{array}{}\text{(1)}& \begin{array}{rl}h\left(x,\phantom{\rule{0.125em}{0ex}}y,\phantom{\rule{0.125em}{0ex}}z,\phantom{\rule{0.125em}{0ex}}t\right)& ={z}_{\mathrm{h}}\left(x,\phantom{\rule{0.125em}{0ex}}y\right)+\frac{p\left(x,\phantom{\rule{0.125em}{0ex}}y,\phantom{\rule{0.125em}{0ex}}z,\phantom{\rule{0.125em}{0ex}}t\right)-{p}_{\mathrm{b}}\left(x,\phantom{\rule{0.125em}{0ex}}y,\phantom{\rule{0.125em}{0ex}}z,\phantom{\rule{0.125em}{0ex}}t\right)}{\mathit{\rho }\left(x,\phantom{\rule{0.125em}{0ex}}y,\phantom{\rule{0.125em}{0ex}}z,\phantom{\rule{0.125em}{0ex}}t\right)g}\\ & ={z}_{\mathrm{h}}\left(x,\phantom{\rule{0.125em}{0ex}}y\right)+{h}_{\mathrm{p}}\left(x,\phantom{\rule{0.125em}{0ex}}y,\phantom{\rule{0.125em}{0ex}}z,\phantom{\rule{0.125em}{0ex}}t\right),\end{array}\end{array}$

where $\left(x,\phantom{\rule{0.125em}{0ex}}y,\phantom{\rule{0.125em}{0ex}}z\right)$ are the Cartesian coordinates (m) of the measurement point, t is time (s), zh is referred to as elevation head (m), p is the total groundwater pressure (Pa) and pb is the barometric pressure (Pa), ρ is the groundwater density (kg m−3) across the water column above zh, and g is the gravitational constant (≈9.81 m s−2). The term hp (m) is the pressure head.

The dependence of the variables in Eq. (1) on $\left(x,\phantom{\rule{0.125em}{0ex}}y,\phantom{\rule{0.125em}{0ex}}z,\phantom{\rule{0.125em}{0ex}}t\right)$ has been deliberately emphasised to stress the point that their magnitude varies in space and time. Determining hydraulic head time series requires four measurements (hereafter also referred to as steps), which are conceptualised in Fig. 1 and can be summarised as follows:

1. geo-positioning or relative positioning of the GMI, i.e. determining its location at the Earth's surface, ${s}_{g}=\left({x}_{g},\phantom{\rule{0.125em}{0ex}}{y}_{g},\phantom{\rule{0.125em}{0ex}}{z}_{g}\right)$ (Sect. 3),

2. establishing the point of (or location representative of) head measurement ${s}_{\mathrm{h}}=\left({x}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}{y}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}{z}_{\mathrm{h}}\right)={s}_{g}+\mathbf{\Delta }{\mathbit{s}}_{p}$, with $\mathbf{\Delta }{\mathbit{s}}_{p}=\left(\mathrm{\Delta }{x}_{p},\mathrm{\Delta }{y}_{p},\mathrm{\Delta }{z}_{p}\right)$ being the vector that represents the location offset sh with respect to sg (Sect. 4),

3. measurement of the water depth below the top of casing dw(tj) at a discrete times tj (open GMI only; Sect. 5),

4. automated pressure measurements at PT location ${s}_{\mathrm{pt}}=\left({x}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{y}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{z}_{\mathrm{pt}}\right)$ of ppt(spt, ti) at discrete times ti (Sect. 6).

There are two methods to obtain $h\left({x}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}{y}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}{z}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}t\right)=h\left({s}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}t\right)$ based on field measurements.

Method 1 (only for open GMI). When only dw has been measured in the field (for example, by taking regular manual water level measurements) the hydraulic head simply follows from

$\begin{array}{}\text{(2)}& h\left({s}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}{t}_{j}\right)={z}_{g}-{d}_{\mathrm{w}}\left({t}_{j}\right),\end{array}$

where tj is the distinct time at which the water level measurement was made.

Hydraulic head time series are nowadays commonly determined from the pressure readings of a transducer located at elevation zpt (m). The head is then calculated using

$\begin{array}{}\text{(3)}& \begin{array}{rl}h\left({s}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right)& ={z}_{g}-{d}_{\mathrm{w}}\left({t}_{j}\right)\\ & +\left[{h}_{\mathrm{pt}}\left({s}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right)-{h}_{\mathrm{pt}}\left({s}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{t}_{j}\right)\right],\end{array}\end{array}$

where hpt is the transducer pressure head, i.e., the pressure recorded by the PT expressed as a water column height (e.g., Hölting and Coldewey2013):

$\begin{array}{}\text{(4)}& \begin{array}{rl}{h}_{\mathrm{pt}}\left({s}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right)& =\frac{{p}_{\mathrm{pt},\phantom{\rule{0.125em}{0ex}}\mathrm{abs}}\left({s}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right)-{p}_{\mathrm{b}}\left({s}_{\mathrm{b}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right)}{{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}\left({t}_{i}\right)g}\\ & =\frac{{p}_{\mathrm{pt}}\left({s}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right)}{{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}\left({t}_{i}\right)g},\end{array}\end{array}$

where ppt, abs and ppt are, respectively, the absolute and relative transducer recorded pressures, and ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ is the average density (kg m−3) across the water column above the transducer's elevation zpt. Application of Eq. (4) is referred to as barometric compensation. The location of the barometric pressure measurement ${s}_{\mathrm{b}}=\left({x}_{\mathrm{b}},{y}_{\mathrm{b}},{z}_{\mathrm{b}}\right)$ must be chosen so that it is representative of the barometric pressure experienced by the PT .

Method 2 (for open and closed GMI). For a PT installed at location spt=sh, the hydraulic head follows from

$\begin{array}{}\text{(5)}& h\left({s}_{\mathrm{h}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right)={z}_{\mathrm{pt}}+\frac{{p}_{\mathrm{pt}}\left({s}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right)}{{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}\left({t}_{i}\right)g}={z}_{\mathrm{pt}}+{h}_{\mathrm{pt}}\left({s}_{\mathrm{pt}},\phantom{\rule{0.125em}{0ex}}{t}_{i}\right).\end{array}$

This is the only way by which heads can be measured in closed GMI for which dw cannot be determined.

For open GMI ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ can be measured. For closed GMI, however, ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ is the average density of the groundwater above zpt, which has to be estimated in the absence of direct measurements. Because the PT is at elevation zpt=zh, ${p}_{\mathrm{pt}}=p-{p}_{\mathrm{b}}$, and Eq. (5) is identical to Eq. (1) when ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}=\mathit{\rho }$. These considerations have important implications when density effects influence the pressure–head relationship of GMI (Sect. 6.4.2).

Hydraulic head is a scalar quantity, and the gradient of the head field in combination with hydraulic conductivity enables quantification of groundwater flow rates using Darcy's law. In three dimensions the hydraulic head gradient (or simply head gradient) is a vector defined as (e.g., Domenico and Schwartz1997)

$\begin{array}{}\text{(6)}& \mathrm{\nabla }h=\mathbit{i}\frac{\partial h}{\partial x}+\mathbit{j}\frac{\partial h}{\partial y}+\mathbit{k}\frac{\partial h}{\partial z},\end{array}$

where the bold italic i, j and k symbols denote the unit vectors in the x, y, and z direction, respectively. Since h and h are continuous field variables, and, in practice, h can only be measured at discrete points sh, head measurements can only be used to approximate h. Moreover, it is rare for field studies to determine h in three dimensions. Therefore, for the purpose of error propagation (Sect. 7), we consider the horizontal (in the xy plane) and vertical components (indicated by either a superscript h or v, respectively) separately by

$\begin{array}{}\text{(7)}& {\left(\frac{\mathrm{d}h}{\mathrm{d}s}\right)}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}\approx \frac{\mathrm{\Delta }h}{\mathrm{\Delta }{s}_{\mathrm{h}}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}},\end{array}$

where the term on the left-hand side represents the rate of head change per unit of distance s, which is approximated by the ratio of Δh, the head difference between two points of measurement, over

$\begin{array}{}\text{(8)}& \mathrm{\Delta }{s}_{\mathrm{h}}^{\mathrm{h}}=\sqrt{\mathrm{\Delta }{x}_{\mathrm{h}}^{\mathrm{2}}+\mathrm{\Delta }{y}_{\mathrm{h}}^{\mathrm{2}}},\end{array}$

or

$\begin{array}{}\text{(9)}& \mathrm{\Delta }{s}_{\mathrm{h}}^{\mathrm{v}}=\mathrm{\Delta }{z}_{\mathrm{h}},\end{array}$

where Δxh, Δyh and Δzh are the distances between two points of head measurement in the x, y and z direction, respectively.

It must be emphasised that considering the horizontal head difference between two points is only meaningful when they are located along the direction of the maximum rate of head change, i.e. perpendicular to the contour planes of equal head (assuming isotropic and constant-density conditions). Hydraulic head measurements from at least three different locations, which are best arranged in the form of an equilateral triangle, are required to determine the head gradient in two dimensions (e.g., Freeze and Cherry1979) or four locations in three dimensions . Even more locations are required for head contour maps (e.g., Ohmer et al.2017). For accurate vertical gradients it is important to use short screens that are within a single hydrogeological unit.

## 2.3 Barometric effects

The following discussion is only applicable to open GMI (i.e., open to the atmosphere; Location 1 in Fig. 1). Air pressure changes are transmitted instantaneously to the water column in open GMI. In contrast, the formation response is more complex because air pressure changes must propagate through the subsurface to the point of measurement, which can result in a delay. Barometric pressure can change as part of the local weather (e.g., the passing of high- and low-pressure systems) by as much as the 1.5 m water level equivalent for the most extreme weather events. If a barometric pressure change propagates through the unsaturated zone of an unconfined system without delay, the water level in an open GMI is a direct representation of the groundwater pressure. However, since the unsaturated zone can resist air movement, for example under low (air) permeability or variably saturated conditions (e.g., Weeks1979), there can be a time lag between barometric pressure changes and the associated GMI water level response (e.g., Rasmussen and Crawford1997). This can be quantified using the barometric response function, which can change over time .

In addition to this, the response to air pressure changes of an open GMI's water level is fundamentally different than the response of the hydraulic head due to the elastic storage behaviour of the subsurface. This can be understood by considering that an increase in barometric pressure raises the total stress acting on both the GMI's water column and the subsurface. The additional stress is borne exclusively by the water column inside the GMI, whereas it is shared between the water and the formation in the surrounding subsurface . As a result, the pressure increase inside the GMI is larger than the groundwater pressure increase, which induces water flow from the GMI into the formation, thus leading to a lowering of the measured water level. The result is an inverse relationship between changes in water level inside open GMI and the changing barometric pressure (e.g., Meinzer1939; Gonthier2003). This relationship can be exploited to detect aquifer confinement but also necessitates the correction of water levels measured in open GMI to faithfully infer the hydraulic head in the formation.

The barometric efficiency (BE) expresses the ratio between the water level change in a GMI Δhpt and the barometric pressure change Δpb causing it :

$\begin{array}{}\text{(10)}& \mathrm{BE}=-\frac{\mathrm{\Delta }{h}_{\mathrm{pt}}}{\mathrm{\Delta }{p}_{\mathrm{b}}}{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}g=\frac{\mathrm{\Delta }h}{\mathrm{\Delta }{p}_{\mathrm{b}}}{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}g=\frac{n\mathit{\beta }}{n\mathit{\beta }+\mathit{\alpha }},\end{array}$

where n is the total porosity of the formation (–), β is the compressibility of water ($\approx \mathrm{4.59}×{\mathrm{10}}^{-\mathrm{10}}$ Pa−1) and α is the (undrained) compressibility of the formation (Pa−1). The minus sign is due to the discussed inverse relationship between hpt and pa.

The BE quantifies the partitioning of the total stress change between the formation and the groundwater . If the subsurface is assumed to be incompressible (α=0 so BE=1, an often-made assumption), the inverse relationship between water level measured in the GMI and hydraulic head in the subsurface is most pronounced. However, the majority of geological materials are more compressible than water (β>α), so realistically $\mathrm{0}<\mathrm{BE}<\mathrm{1}$ . Methods for reducing barometric effects on hydraulic head measurements were suggested in the literature and are referred to as barometric correction (not to be confused with barometric compensation; Eq. 4; e.g., Hubbell et al.2004; Toll and Rasmussen2007; Noorduijn et al.2015). This discussion highlights that the BE of a formation is an important property, and ignoring it can have significant implications when hydraulic heads or gradients are derived from water level measurements with the aim to interpret groundwater processes (Spane2002).

Avoiding barometric effects requires GMI with a specific design. suggested a sealed well and showed that their design reduced barometric pressure effects by an order of magnitude, especially for sites with deep vadose zones. Furthermore, a laboratory study by demonstrated that measured total pressure recorded in sealed and unsealed wells is equal assuming that barometric pressure is also measured; water levels can be accurately measured in either sealed or unsealed standpipes. This is convenient for fluvial environments, where long standpipes are subject to the forces of river flows, which can be quite violent in ephemeral streams especially (e.g., Shanafield and Cook2014).

## 2.4 Clarification of error terminology

Figure 2Possible combinations of accuracy, precision and resolution illustrated in a matrix, when 1000 measurements of the same head (h=1 m) are made. Measurements are (a) inaccurate and imprecise, (b) inaccurate and precise, (c) accurate and imprecise, and (d) accurate and precise. Examples are illustrated with two different values of accuracy, precision and resolutions (equal to bin width in histogram).

Figure 3Illustration of the influence that the instrument resolution has on the measurement error: a continuous, time-variable head is measured at discrete time intervals by instruments with analogue-to-digital conversion resolution of 5, 2 and 0.1 mm.

Despite their importance, the terms related to measurement error are often mixed up or used ambiguously. Thus, before proceeding, it is crucial to clarify their meanings within the context of head measurement.

• Accuracy is a measure of how closely the mean of the measured head corresponds to the real head. The deviation between the true value and the mean of its measurements is the systematic (or absolute) error (Fig. 2).

• Precision is the spread of the measured heads around their mean value. When the measurements are normally distributed, it can be expressed by the standard deviation of a Gaussian distribution. It is also referred to as random error (Fig. 2).

• Resolution is the smallest numerical separation at which the change of real value can be distinguished.

• Range is the difference between the minimum and maximum value an instrument can measure.

Electronic measurements use analogue-to-digital converters (ADCs), which convert continuous analogue signals into discrete (digital) values. ADCs generally have limited steps (resolution bins; Fig. 3), leading to an inverse relationship between the measurement range and resolution. Consequently, the larger the range of measurement, the coarser the resolution. For example, a 12 bit ADC has 212=4096 resolution bins, which equates to a theoretical resolution of 2.4 mm, when the range is 10 m, or a resolution of 12.2 mm, when the range is 50 m. As Fig. 3 demonstrates, the difference between the continuous and instrument-reported (quantised) head, and thus the measurement error, decreases with increasing resolution.

3 Geo-spatial positioning of groundwater monitoring infrastructures

There are two ways to determine a GMI's position (sg; Fig. 1). The first is surveying, which is the determining of the three-dimensional distance between points of interest. The second is to use navigation satellites. This section briefly summarises both. More details on surveying can be found in and on satellite system technology and applications in , , and .

## 3.1 Relative positioning using traditional surveying

Determining the horizontal and vertical distances to a reference point (known as trigonometric levelling) can be done using a total station theodolite. These are equipped with a precision telescope that can rotate in the horizontal and vertical direction, allowing visual adjustment of the telescope to points of interest. Precise optical sensors can pinpoint a barcode on the staff and digitise the angle and azimuth readings from which the horizontal and vertical distances are calculated using a built-in computer. They further include an electronic distance measurement (EDM) device, based on the travel time of laser pulses reflecting off a target, and have satellite receivers to determine geo-coordinates (Sect. 3).

Levelling, the technique of measuring vertical distances (heights) relative to a known survey benchmark, can be conducted using optical or light-based instruments operating from a tripod. The latest generations of optical levelling instruments use a rotating precision telescope to magnify the scale printed onto a levelling rod (staff) that is held vertically on top of a point of interest. The telescope is used to read the vertical distance above the point of interest of a laser beam rotating in a horizontal plane. The levelling rod is equipped with a receiver that can be moved vertically until it detects the beam.

The maximum measurement distance of digital levels or total stations is limited to hundreds of metres, depending on the telescope, the range of the laser beam and the visibility of the target . Longer distances are surveyed by leap-frogging survey devices along multiple points (traversing; Brinker1995). Measurement error is a function of distance, and accuracy and precision of leap-frog surveys tend to be poorer than surveys where the instrument does not require moving. An indication of the measurement error can be obtained by returning to the starting location of the survey and determining the difference between the recorded positions at the start and end. When GMI locations are to be referenced with respect to a national datum, the accuracy is further dependent on the quality of the known benchmarks that provide the link between the local survey to the national datum (Fig. 1).

It is difficult to determine the accuracy of high-precision surveying because this must be compared to a more accurate benchmark method. The measurement error for state-of-the-art survey devices depends on many factors, including instrument set-up, calibration, sun position, temperature elevation gradient, battery level and, most importantly, operator's expertise . The literature contains very few peer-reviewed investigations that test manufacturers' specifications. However, one assessment has illustrated that digital levelling can reach an accuracy of 2 mm km−1 with precision of 1 mm + 1 mm km−1 . Leap-frogging using 150 m distance steps found an elevation precision of 1.9 mm km−1 .

Estimating the positioning errors of total stations is even more complicated due to the combination of EDM and angle sensors . thoroughly investigated the accuracy and precision of industry-standard EDM devices over a well-calibrated distance of 40 m. They found that the accuracy varied from 0–4 mm, with some devices showing dependence on the measurement distance. We use the precision of 0.5 mm stated by for our error analysis (Table 1).

## 3.2 Navigation satellite positioning

Global navigation satellite systems (GNSSs) currently available include the widely used Global Positioning System (GPS; USA) and Globalnaya Navigazionnaya Sputnikovaya Sistema or Global Navigation Satellite System (GLONASS; Russia) as well as Galileo (European Union) and BeiDou (China), which are currently being deployed. Additionally there are local systems such as the Indian Regional Navigation Satellite System (IRNSS; India) and the Quasi-Zenith Satellite System (QZSS; Japan). Each system type consists of a network of satellites that orbit the Earth at 18 000–25 000 km altitude.

The network satellites transmit their location and absolute, synchronised time, encoded in radio signals with at least two different frequencies. A GNSS receiver can decode these signals and calculate the distance to multiple satellites using the signal arrival times. In the case of global systems, the intersect of distances from at least four individual satellites enables a GNSS device to calculate a location in geo-coordinates via trilateration. Single-point positioning (SPP) requires only one GNSS receiver (Hegarty2017). The horizontal positioning accuracy is at best within 5–8.5 m , and the vertical accuracy is poorer still. This is because the visible satellites are more closely aligned in a horizontal plane and the Earth shields the signals from remaining satellites, which would provide more vertical information. Recent developments have focussed on eliminating the need for multiple GNSS receivers and speeding up the time required to achieve accurate positioning .

Measuring locations relies on a reference system (georeferencing) that is Earth-centred and Earth-fixed. A catalogue of 3-D positions is given by the International Terrestrial Reference Frame (ITRF). This falls to within ±1 m of the 1984 (WGS84) and is therefore used as the common reference frame for geo-positioning . The International Hydrographic Organization mandates the use of WGS84 as the horizontal reference for hydrographic mapping (Rizos2017).

Geo-positioning is based on the geographical coordinate system, which delivers the spherical coordinates of latitude, longitude and height (geoidal geometry as a global reference point). Measuring lengths and areas in spherical coordinates is not straightforward. For the purpose of hydrogeological investigations, these coordinate points are transformed into a projected coordinate system, a 2-D representation of the Earth's surface. Although, there is some uncertainty as to the origin of this projection , the most commonly used projection is the Universal Transverse Mercator (UTM) system, which divides the Earth into 60 zones and 20 latitude bands. Each zone is then assumed to be planar, and coordinates are expressed in metres as northing, easting and elevation (projected from the geoid to a flat surface with the local zone as the reference point). Note that height (vertical distance above the ground surface) and elevation (vertical distance above sea level) should not be confused.

Differential global navigation satellite system (DGNSS) positioning can provide much better accuracy and precision than GNSS. This approach requires at least two GNSS receivers, one of which is stationary and located at a known point (base station). The base station uses single-point positioning in conjunction with its known location to calculate an error correction. The second, mobile GNSS receiver (rover) uses the GNSS signals in conjunction with the error correction to calculate its distance from the base station. The error correction is determined from signal phase observations at both stations (Remondi1985). This can be achieved offline by post-processing the stored satellite signals in both receivers or in real time through a radio link between the rover and the base.

Recent developments in many countries have resulted in continuous operating reference stations (CORSs) at strategic locations whose error corrections can be accessed via mobile data networks as long as there is network coverage. The most sophisticated GNSS devices can nowadays provide positioning with millimetre horizontal and sub-centimetre vertical accuracy . However, these innovations have yet to make it into commercial receivers. More typically, best-achievable horizontal accuracy and precision are 15 and 10 mm, respectively, whereas vertical accuracy and precision are 30 and 40 mm, respectively . These numbers have been adopted for the purpose of our error analysis (Table 1). Interestingly, found that accuracy and precision did not show any dependence on the distance to the base station within their test area of about 11 km. It should be noted that these accuracies are achievable only when there are a sufficient number of visible satellites (for both receivers). When points of interest are near or under vegetation, the geo-positioning accuracy is significantly degraded .

When traditional surveying is undertaken there appears to be a horizontal or vertical distance dependent error, whereas for DGNSS this is not the case (Table 1). Using the random error estimates, a horizontal cut-off distance at which the precision from state-of-the-art DGNSS is better than that of a total station theodolite is ≈700 m. For vertical distances, DGNSSs become more precise than digital levelling when two locations are further apart than ≈15 km in the horizontal direction. In this case it is meaningless to derive vertical head gradients. Consequently, the surveying approach should be chosen according to the distance between the locations.

Table 1Summary of precisions for the four different steps and methods required to calculate hydraulic heads and gradients. For a graphical explanation, see Fig. 1. Values are best possible estimates or collated from the literature.

4 Point of head measurement

## 4.1 Representative point of measurement

However, larger errors can be expected when vertical gradients are higher than the example value used, which may be the case near groundwater discharge zones, under pumped conditions, and in formations of low permeability. found that for aquitard hydraulic conductivity tests, where the propagation of an induced head drop in a piezometer is recorded as a function of time, the representative point of measurement was biased towards the bottom of the screen and that this significantly influenced the outcomes of the parameters to be determined. Moreover, as the gradients changed in time, so did the representative point of measurement. In layered aquifer systems, the water level in wells with long screens was found to depend on the transmissivities of the layers intersected by the wells (Sokol1963). These findings highlight the need for using short screens. However, the finite screen length of standpipe piezometers means that some uncertainty remains about the representative vertical position of the head measurement.

## 4.2 Borehole verticality and screen location

Despite best efforts in many countries on the mandatory requirement to report accurate information on the drilling and completion of GMIs, construction details are often reported at a precision of decimetres (not centimetres) and prone to significant systematic error. Due to the variety of different field and environmental conditions as well as the different qualifications and experience of drillers, we assume that the vertical screen locations can be estimated from driller's logs, with a precision no better than about 0.5 m (Table 1).

The deviation from the vertical of a GMI further results in uncertainty about sh (NUDLC2011). The importance of borehole deviation surveys is critical in other industries such as oil and gas, where errors in the observed inclination angle and other parameters of the monitored fracture-system geometry impact the monitoring and interpretation of hydraulic fracturing . Yet in hydrogeology borehole verticality is typically ignored, in particular when calculating heads gradients. Poorly aligned boreholes impact significantly the integrity of casing and hence increase the risk of flow short-circuiting and water column density stratification (Sect. 6.4.2).

A thorough investigation of borehole deviation was conducted by , who applied a correction factor to water level data when a borehole deviation survey indicated a change of more than 0.06 m between the measured borehole length and the true vertical depth. From the 177 boreholes surveyed, correction factors to the historical water levels of these wells ranged from 0.06 to 1.8 m, and inclination angles ranged from 1.6 to 16. A comprehensive examination of borehole deviation was conducted in more than 100 boreholes drilled (up to 1000 m deep) at the Swedish nuclear repository site by the Swedish Nuclear Fuel and Waste Management Company (SKB; Nilsson and Nissen2007). Their investigation provided an uncertainty of deviation measurement of the inclination of the borehole (up to 3) as well as an elevation uncertainty at the bottom of the borehole (up to 15 m for the boreholes measured).

Guidelines for drilling and water bore construction for plumbness and straightness is generally a “do the best you can” approach within practical limits using appropriate equipment and drilling operation (e.g. drilling centralisers, correct collar and feed pressure) for the geological conditions . Drillers consider angles of less than 5 to be acceptable . Hence, the horizontal positioning error of a 10 m long borehole would become $\mathrm{sin}\left(\mathrm{5}{}^{\circ }\right)\cdot \mathrm{10}$ m ≈872 mm, whereas the vertical error is $\left[\mathrm{1}-\mathrm{cos}\left(\mathrm{5}{}^{\circ }\right)\right]\cdot \mathrm{10}$ m ≈38 mm. In the absence of more literature reporting on borehole deviations, we use this reported figure as the random error when determining the point of head (Table 1).

Since this uncertainty is much greater than the achievable accuracy of the GMI's geo-position, sg, the verticality (plumbness) of a borehole should be measured using downhole geophysical tools such as verticality probes or inclinometers. This includes a gyroscope or an accelerometer to measure the vertical angle combined with a magnetometer to provide the probe's rotational position around the vertical axis. Both measurements can be logged continuously while lowering the probe. For example, assuming an industry-standard precision of 0.5 (e.g., Verticality Sonde by GeoVista, UK) and an otherwise straight 10 m deep borehole, the resulting precision in identifying the horizontal screen offset would be ≈87 mm, which is an improvement over the crude guess based on a 5 angle according to best drilling practice. The borehole deviation survey should be combined with a downhole camera to determine the position of the screen relative to the top of the GMI. We estimate the depth measurement precision of a typical system to be approximately 20 mm (Table 1).

5 Depth-to-water measurements

There are a number of different ways to measure depth to water (dw, Eq. 3). It is commonly done by hand and involves the use of a measurement tape . Most groundwater projects today use electric water level meters, colloquially called dip meters, which provide an audible or visible signal when a sensor touches the water surface. When the acoustic signal is not electronic but an audible noise is generated mechanically, for example by lifting and dropping a hollow brass cylinder just touching the water surface, the instrument is called a plopper. Another inexpensive method uses a steel tape that is covered with chalk .

Depth-to-water measurements should be performed frequently (at a minimum every 3 months) for checking the performance, and adjustment, of automatic sensors. Good-quality measuring tapes are marked every 1 mm (metric) or every 0.01 ft (imperial). The chalked-tape method can potentially deliver a precision that corresponds to the resolution of the graduated steel tape (Table 1 and Fig. 2), whereas dip meters and ploppers are generally read to the nearest centimetre. This may involve human measurement errors such as the switching of digits (e.g., noting 57 instead of 75) or reading to the wrong decimetre or metre marking on the tape .

There is only minor information in the literature about the errors associated with manual head measurements. The lack of assessment is surprising given that manual measurement is the most important link which ties automated pressure time series to a benchmark (Eq. 3). Some controlled experiments have been conducted though, most recently in the Netherlands by . Sixteen operators, with varying degrees of experience, each took a reading in a total of 16 standpipes. Half of the readings were done with an electronic dip meter, and the other half were done with a plopper. After discarding the obvious mistakes from the data set, the errors were fitted to a normal distribution with a mean and standard deviation of 5 mm and ±8.4 mm, respectively (Table 1), for the electronic dip meter, versus 0.3 mm and ±9.5 mm for the plopper. The measurements by were representative of very shallow water tables. A poorer precision (0.05 ft =15 mm) was reported by (cited in ) for water levels more than 120 m below the surface measured by two observers using the same instrument within a short time period. conducted a rather comprehensive experiment themselves but reported the errors as percentages, which makes the figures difficult to compare to the other studies.

Figure 4Manual measurement of nine different groundwater depths using eight different dip meters: difference from mean and standard deviation (precision) over mean depth.

noted that the graduations on some of the tapes that were used showed noticeable differences and that this caused systematic measurement error. compared three water level probes to a reference probe and found differences that increased with depth, reaching a maximum value of 0.1 ft =0.03 m at a depth of 61 ft =19 m. Comparable findings based on our own experiments are shown in Fig. 4, which summarises variations in manual measurements using several commercial electric dip meters of various lengths to measure water levels at nine depths between 5 and 90 m. Electric water level loggers by the same manufacturer differed by several centimetres, with differences of up to 0.12 m observed overall (one person taking the measurements sequentially for each borehole, using the same location at the top of casing for each bore). The differences increased with depth to the water table for several instruments, confirming the observations made by 25 years ago. Discrepancies of this magnitude preclude the use of data for accurately identifying small head gradients and call for replacement of the measuring instrument.

Wear (e.g. kinks and tears) on electric water level tapes causes additional discrepancies in measurements over time. detail procedures for calibrating electric water level devices before each use; this involves measuring the electric tape against a steel measuring tape kept in the office only for this purpose. For consistent time series, it is extremely important that a datum on the casing is used to always measure water level from the same point . This seems obvious, but there can be confusion when different operators are involved, and regular maintenance is required to make sure the mark stays visible. Repeating the measurement a number of times can avoid tape-reading errors and ensures proper functioning of the electronic dip meter, which may give inconsistent readings sometimes .

In areas prone to vertical land surface movement a regular check of the elevation of the casing (zg) is necessary. The causes for such movements can be manifold and include tectonic processes, slope instability, freezing and thawing cycles , or clay swelling and shrinking with changing moisture conditions. In peat areas, subsidence by compaction, oxidation or drying is a well-known cause for movement of the well casing . Moreover, damage can occur to standpipes either by vandalism or natural processes, such as ice expansion when the water inside the GMI freezes .

6 Automated water level time-series measurement

## 6.1 Automated measurements

Automated measurement of water levels or pressures in GMI requires electronic devices capable of time keeping and sensing. Many commercial instruments combine a stand-alone clock, a sensor, an ADC unit, memory and a power supply in a single housing. There are also instruments that house only the sensor and are connected to a data logger that converts the sensor signal and stores the quantised readings. The focus of this section is on systematic errors that occur during automated collection of water level data in the field.

Automated instruments have the capacity to record unattended for a long period of time. Data loss as a result of logger chip or battery failure can be prevented by the automated transmission of instrument readings to a receiving data management system via radio, infrared signals, a GSM network or satellites . This is referred to as telemetry. The expansion of cellular-network-provider coverage and the reduced cost of data-only plans in the recent past have made telemetry systems a more accessible and viable option for remote hydrogeological monitoring. Nowadays, transmitted data are stored on network servers which can be accessed in real time via computer or smart phone. Telemetry systems allow remote modification of sensor settings and identification of sensor problems, early identification of logger failure as a safeguard against data loss, and re-synchronisation of clocks to avoid time-based errors (Sect. 6.6). However, the deployment of a telemetered system does not avoid certain errors such as sensor drift (Sect. 6.4.3). Consequently, to ensure the accuracy of automated water level measurements, frequent site inspections and manual measurements are still required.

## 6.2 Types of devices

At present, water level time series are typically determined from pressure measured using submersible or grouted-in PTs (Fig. 1). This relies on Eqs. (3) and (4) to determine the water level. A comprehensive overview of PTs can be found in or .

The most popular type of PT consists of a piezo-resistive crystal made from silicone or ceramics, which acts as a strain gauge as it deforms under pressure. The deformation causes the electrical resistance of a Wheatstone bridge to change, which is gauged by recording the changing voltage due to a constant current. Vented PTs are connected to a venting tube that connects the air chamber of the submerged PT to the atmosphere. The measured pressure is the relative pressure (ppt in Eq. 4), and there is no need for barometric compensation. We estimate a best-possible precision for vented PTs to be 1.5 mm (Table 1). This value is based on water level time series recorded with three vented PTs (with a range of either 10 or 20 m) inside a standpipe piezometer in Hanover (Germany) over a period of more than 15 months. Using one logger as a reference and calculating the difference with the remaining two loggers resulted in standard deviations of 1.4 and 1.7 mm (n=11 279), thereby demonstrating excellent performance. However, readings can be influenced when the venting tube does not remain dry. A desiccant capsule is therefore attached to the tube, and this can causes some practical difficulties when measuring in riverbeds or other areas subject to flooding as well as when freezing occurs .

Non-vented PTs measure absolute pressure ppt, abs, which are converted to relative pressure ppt by subtracting the barometric pressure pb (Eq. 4; barometric compensation), typically measured with a barometric PT near the GMI. The subtraction of barometric pressure from absolute pressure results in a loss of water level measurement precision because the two instrument measurement errors accumulate (Sect. 7). Because of this, we estimate the best possible water level measurement precision as 3 mm (double that of vented PTs; Table 1).

A second type of PT is the so-called vibrating wire piezometer (VWP) that uses electromagnetic coils to excite a wire exposed to differing strain resulting from pressure changes. The square of the resonant frequency is linearly proportional to the pressure . VWPs are designed for long-term stability and are therefore used for closed GMI when the instrument is fully grouted-in (Fig. 1, Location 2). However, recalibration becomes impossible once installed . One strategy to verify PT performance in that case (if the budget permits it) is to install three instruments at the same depth. VWPs of the non-vented and vented type exist, and some models have a pressure range of 10 MPa (equivalent to about 1 km of water) or sometimes higher. We estimate the best possible water level measurement precision of VWPs to 7 mm (Table 1).

Wet–wet pressure transducers measure the pressure difference between two points that are both exposed to water . Such devices are ideal for obtaining small head gradients, such as is required for measuring surface water–groundwater interactions, because they eliminate the uncertainties arising from barometric correction or the spatial positioning of two individual measurement points.

There are also electronic water level measurement devices that emit a laser pulse and determine the depth of water from the time it takes the pulse to reach the water table and return to the sensor (known as lidar: light detection and ranging). When connected to a time-keeping data logger, this technology is suitable for time-series collection . A recent development and test of a lidar-based system demonstrated an outstanding precision of 0.5 mm (Table 1Benjamin and Kaplan2017), but condensation and, in groundwater studies, borehole non-verticality can interfere with the light reaching the water surface.

Another type of water level sensor is based on electronic capacitance measurement. It consists of two electrically isolated plates or wires that are aligned in parallel at close proximity. Submergence of the wires in water creates a contrast in electrical capacitance compared to air, with values that are proportional to the submerged length. Their range (typically 1 to 2 m of water level) is smaller than piezo-resistive PTs (which can be used in water depths of 100 m or even more). An important advantage of capacitance probes is that they are rugged and can withstand overload, drying and freezing. In contrast to the measurement-tape and lidar techniques, which measure dw from the top downward, the capacitance probe sits in the water column and records water levels on a data logger (similar to PTs).

## 6.3 Instrument range and resolution

In their instrument specifications, manufacturers typically provide the accuracy of a PT as a percentage of the range, or full scale (FS). Unfortunately, this number is not defined unambiguously. Typically, it may comprise a combination of a sensor's non-linearity (the relationship between p and voltage V not being a straight line), hysteresis (differing pV relations during p increases or decreases) and repeatability (the closeness of measured p values for the same V), and thermal artefacts (influence of temperature on the pV relation). These are non-adjustable errors and are therefore not related to accuracy in the sense that they can be corrected by applying a simple offset to calibrate the instrument to a known value. Moreover, since there are different ways to quantify the instrument's deviation from the ideal pV relation, the number specified as the instrument's accuracy can have a different meaning depending on a manufacturer's definition. This can even mean that an instrument with 0.5 % FS accuracy is as accurate as an instrument with 0.1 % FS accuracy .

Figure 5(a) Depth to water level (WL) measured by three different instrument types, illustrating the influence of precision and resolution on head time series: capacitance PT (blue line; high precision and low resolution), vented PT (orange line; high precision and high resolution) and non-vented PT (corrected for barometric pressure; green line; low precision and high resolution). The examples are from (a) Ti Tree basin in the Northern Territory (Australia) and (b) a farm dam in South Australia (note the effect of a 0.2 mm rainfall event on 5 January, which is visible in the orange line but not in the green line). The low precision achieved by the non-vented PT is specific to the particular instrument used and not representative of all PTs of this type.

As a practical example of precision and resolution (see Sect. 2.4), Fig. 5a and b show several days of automated depth to water level measurements made using different logger types. The difference between the graphs is the vertical scale, with the water levels in Fig. 5a fluctuating at the millimetre scale and those in Fig. 5b showing an overall decline of a few centimetres. The curves recorded by different water level measurement devices illustrate the discrete time and magnitude nature of automated measurement. Here, the blue line illustrates high-precision and low-resolution values, and the orange line shows high-precision and high-resolution measurements, whereas the green line represents low-resolution and low-precision data.

Current standard practice allows the quantification of subsurface processes and properties from time changes in heads due to either natural or induced causes. Common examples such as aquifer tests rely on large head changes that are sufficiently resolved using off-the-shelf instruments. However, more recent research advances have demonstrated that subtle signals in hydraulic heads can also be used to passively quantify hydrogeological processes and properties. For example, Fig. 5 demonstrates sub-centimetre diel (i.e., daily) fluctuations that originate from phreatophyte evapotranspiration (e.g. Gribovszki et al.2013) or Earth and atmospheric tides . Such subtle signals can only be detected with appropriately high sensor resolution. Currently, it is advisable to deploy vented transducers to minimise errors resulting from clock differences and imprecisions due to barometric compensation (refer to Sect. 2.1). For such intentions the measurement range must be minimised in favour of maximum resolution, which reduces the measurement error (Fig. 3). It should be noted that when PTs are used to calculate gradients, the readings from non-vented PTs may be used directly without compensating for atmospheric pressure changes as long as the PTs all experience the same atmospheric pressure change.

## 6.4 Issues related to pressure transducers

### 6.4.1 Temperature effects

Figure 6(a) Barometric pressure and temperature and (b) water level versus time. The data in (a) were collected using a logger that was exposed to significant temperature fluctuations prior to 16 November 2013. The water head time series in (b) was derived by subtracting the measured pressures shown in (a) from the total pressures measured using a non-vented PT. The artefacts caused by the temperature variations prior to 16 November 2013 are clearly reflected by diurnal oscillations of the water levels (grey shaded area). Manual dips are indicated by blue dots.

The response of piezo-resistive sensors to pressure changes is a function of temperature; hence most instruments record temperature alongside pressure and use this to compensate the readings. Nevertheless, the operation of PTs in transient temperature environments has been found to affect water levels that are calculated from pressure readings. For example, showed that when PTs are exposed to direct sunlight, thermal effects add noise to water level measurements. noted that temperature compensation often significantly compromises the accuracy of pressure readings.

For non-vented PTs, it is especially important to consider placement of the barometric PT to prevent adding noise from thermal effects into water levels during barometric compensation , which is demonstrated in Fig. 6. The graph of water level versus time shows a clear diurnal variation in the time series recorded by a non-vented and barometric PT pair. The data were converted to water levels by subtracting the recorded barometric pressures from the total pressure recorded by the non-vented PT, which was suspended in a surface water pond. Inspection of the diurnal temperature variations of the barometric PT shows daily temperature variations of 10 C or more prior to 16 November, on which date the PT was placed in a more constant temperature environment. As a result, the periodicity that characterised the water level data before that date disappears from the water level time series. argued that such thermal effects should not affect vented PTs, which therefore should be suitable for fine-scale (e.g. sub-daily) measurements as required for evapotranspiration calculations. However, found that rapid changes in temperature on a sub-daily timescale can cause the air in the line to expand or contract and that the relationship between temperature fluctuation and logger error varies between loggers.

### 6.4.2 Water column density changes

Figure 7(a) Time series of water level pressure recorded by two transducers, a shallow (superscript s) and deep (superscript d) transducer, located in the same well. (b) Difference between the deep and shallow pressure measurements. (c) Water level in an observation well showing offset in mean water level after the well was purged for hydrochemical sampling (AHD is Australian Height Datum). Note that due to the high temporal resolution, the discrete measurement points are not shown and the time series are drawn as lines instead.

Figure 8(a) Temperature and average density as a function of depth for an observation well in Japan that experienced significant warming due to urbanisation between 1993 and 2003. Temperature data from . (b) Theoretical head increase due to the temperature-related density decrease and head increase that would be perceived if the PT's vertical position is lowered due to thermal expansion of the steel suspension wire.

For internally consistent hydraulic head time series, it is imperative that the average density across the water column ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ in Eq. (4) stays constant in time. Strictly speaking, this is never the case, and the impact of the changes of ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ represents a systematic measurement error that has to be assessed and, when not negligible, corrected for. The effects are largest in groundwater systems with changes in salinity, such as coastal aquifers . When the density varies within the water column and with time, ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ is given by (e.g.,  Post and von Asmuth2013)

$\begin{array}{}\text{(11)}& {\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}\left(t\right)=\frac{\mathrm{1}}{{h}_{\mathrm{pt}}}\underset{\mathrm{0}}{\overset{{h}_{\mathrm{pt}}}{\int }}\mathit{\rho }\left(z,\phantom{\rule{0.125em}{0ex}}t\right)\mathrm{d}z,\end{array}$

where ρ(z, t) is the density as a function of the vertical dimension and time, and hpt is the vertical distance between the top of the water column and the PT (Eq. 4).

Application of Eq. (11) requires knowledge of the density distribution across the length of the water column at multiple times, which is seldom collected. Pressure-to-head conversion errors due to unknown knowledge of ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}\left(t\right)$ are therefore probably one of the most overlooked issues in head time-series measurement. Some instruments provide a correction function based on the change of the electrical conductivity and temperature measured by sensors housed in the same instrument as the PT, but this is only meaningful if the density of the water column above the PT is constant. It is important to distinguish the effects discussed here from the head corrections that must be applied when studying flow in variable-density groundwater systems .

A subtle effect of the change of ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ with time is shown in Fig. 7. Figure 7a shows the pressures recorded during an experiment whereby two PTs were hanging inside the same standpipe piezometer. One was located just beneath the air–water interface, and one was just above the bottom of the piezometer. The latter case corresponds to the way the pressures are recorded when heads are calculated using Method 2 (Sect. 2.1), whereas the former is representative of Method 1. Because of the well's vicinity to the sea, the recorded pressures vary with the tide. The difference between the PT readings is shown in Fig. 7b. In a well with a constant ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$, the pressure difference would be constant in time. Clearly, this is not the case here, and two effects are notable: (i) a linear trend (grey dashed line in Fig. 7b), causing the pressure difference to become smaller, and (ii) oscillations that are superimposed on this linear trend.

The linear trend was due to leaking casing joints, which led to the ingress of fresh groundwater in the upper parts of the piezometer, as a result of which a salinity stratification developed. As more freshwater seeped in with time, ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ decreased by an amount $\mathrm{\Delta }{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ per unit of time Δt. In fact, the slope of the linear trend line is equal to $\frac{\mathrm{\Delta }{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}g{h}_{\mathrm{pt}}}{\mathrm{\Delta }t}=\mathrm{200}$ Pa d−1, which for hpt≈50 m gives $\frac{\mathrm{\Delta }{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}}{\mathrm{\Delta }t}=$ 0.04 kg m−3 d−1, and this is roughly consistent with the estimate of $\frac{\mathrm{\Delta }{\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}}{\mathrm{\Delta }t}=\mathrm{0.03}$ kg m−3 d−1 derived from consecutive downhole probe measurements on 1 July and 5 August 2015. The superimposed tidal oscillations were caused by the change of the density stratification inside the piezometer standpipe with the tide: as the tide rose, groundwater with an ambient, high salinity entered across the well screen, and this caused more saltwater to stand above the deepest PT. The shallow PT, however, remained in the freshwater part of the stratified water column. Both PTs experienced the same increase in water column height above the sensor, but because this added height consisted of freshwater for the shallow PT, it sensed a smaller pressure change than the deeper PT. Correcting for these effects shows that the pressure difference becomes virtually constant although not all fluctuations disappear. The fluctuations around the mean difference decrease to become around 0.1 kPa (1 mm of water column height). The cause of the remaining fluctuations is not clear; they may be due to clock synchronisation issues.

While the previous example showed a subtle trend of relatively low magnitude, the time series in Fig. 7c show the potentially large magnitude of an abrupt change of ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$. In this case it was caused by the purging of the well for hydrochemical sampling. Prior to sampling, the water inside the well had a non-constant salinity because it had not been properly developed at the time of construction. After sampling, the well was filled with water with the same salinity as the groundwater at the well screen. As a result, ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ increased from 1006.7 to 1015.1 kg m−3. Based on these density values, the length of the water column inside the well would have changed from 72 m to $\frac{\left(\mathrm{72}\cdot \mathrm{1006.7}\right)}{\mathrm{1015.1}}=\mathrm{71.4}$ m, i.e. a decrease of 0.6 m, which corresponds to the measured change of 0.67 m; the additional difference may be due to the removal of silt and other fouling material from the well screen by the pumping.

An example of the effect of temperature-related density changes on the head error is shown in Fig. 8. In this example, the change in temperature was caused by urbanisation, resulting in a noticeable warming of the upper 75 m of the subsurface. This caused the density of the water column to decrease, and hence a longer column of water is required to balance the pressure at the screen. Figure 8b shows the magnitude of this effect as a function of depth, which in this example remains limited to less than 1 cm. Another effect that plays a role is the lengthening of the logger's metal suspension wire as it warms. Assuming a linear expansion coefficient for steel of $\mathrm{11}×{\mathrm{10}}^{-\mathrm{6}}$ K−1, the increase in wire length as a function of depth is shown in Fig. 8c. The magnitude of this effect is less than 1 mm. Because this example was chosen to represent a case of relatively strong temperature increase for groundwater, it is expected that these values represent the upper bounds for thermal expansion effects, which thereby represent relatively small errors under typical groundwater conditions. Larger effects could occur though near-aquifer thermal storage facilities, geothermal areas, or in very deep wells.

### 6.4.3 Measurement drift

Sensor drift is one of the most common errors in automated hydraulic head measurements. Here it is expressed as Δdw, which is defined as the depth below the top of casing (TOC) measured manually with a dip meter minus the depth measured by the PT. tested 14 different transducer brands commonly used in hydrogeological studies. For PTs with a range <15 m H2O, the drift was observed to be $-\mathrm{8}\le \mathrm{\Delta }{d}_{\mathrm{w}}\le \mathrm{27}$ mm after 99 d in the field, but the models with a greater range showed up to 5 times more drift. Data available to the present study from Syria, where 11 observation wells were equipped with vented PTs in January 2009 and were not inspected until June 2010, showed $-\mathrm{199}\le \mathrm{\Delta }{d}_{\mathrm{w}}\le \mathrm{153}$ mm, with only one of the PTs showing a negative Δdw value.

In an extensive study of 473 piezometers, all equipped with the same brand logger and inspected every 3 months for a total of 2 years, statistically evaluated Δdw values were based on a data set of 5583 measurements. For 144 piezometers, a statistically significant linear trend could be identified. The slope of the trend line was negative for 95 and positive for 49 of the piezometers. The drift was reported (as the median of the trend line slopes) to be −3.6 and 4.4 cm yr−1 for the negative and positive trends, respectively.

Apart from technical reasons that cause PT drift, fouling of instruments is a well-known problem that affects the quality of head time series. This can be due to the formation of mineral precipitates by hydrochemical processes . Biological processes often build biofilms of microorganisms, or larger organisms such as snails attach themselves to a sensor. Biofouling filters consist of copper coiled wire that can slow down these effects, but regular inspection and cleaning are a requirement to prevent measurements from being compromised. Moreover, improper suspension cables may stretch, hence causing the logger to sit deeper below the water surface, or frequent removal of loggers for downloading may cause the wire length to change due to kinks. The cables of vented PTs may be large relative to the well diameter, and sometimes there is little room for the desiccation unit at the top, which may mean that the logger is not always returned to exactly the same position after the GMI was accessed for maintenance or other kinds of measurements (e.g. water sampling).

Drift introduces errors of unknown magnitude that remain unnoticed unless identified by frequent checks using an independent measurement . The examples of field-observed drift in and show that drift is not generally linear. The rate of change can vary in time, sometimes suddenly, and even reverse direction. Frequent, independent dw measurements by manual dipping provides the only means to correct for drift. Drift correction involves removal of the linear trend between manual measurements, which introduces uncertainty because of the linearity assumption. Drift corrections must be carefully documented, and at all times the original data must be stored alongside the corrected data. Data downloaded at different times must be stored separately to ensure that a drift correction applicable to a particular block of data is not inadvertently applied to other blocks of data.

Figure 9Water level versus time for a piezometer in the Netherlands that had a clogged well screen until it was rehabilitated in June 1996 . The temporal dynamics of the head in the aquifer were not registered by the piezometer until that time. Data were obtained from http://www.dinoloket.nl (last access: 16 January 2019).

Another form of drift, which is related to the GMI and not the PT itself, occurs if the conditions in the GMI change such that the relationship between the recorded water level hpt in the well and the groundwater pressure p in the aquifer (Eq. 4) is not constant over time. When the change of ${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$ with time (Sect. 6.4.2) is responsible for this, it will cause Δdw≠0. However, Δdw≠0 cannot be used to detect measurement errors due to clogging of the well screen by suspended sediment particles, geochemical processes (e.g., iron oxidation) or biofilm growth. An example of the detrimental effect on time series because of the latter phenomenon is illustrated in Fig. 9, which shows the water level in a piezometer as a function of time. The temporal dynamics remained very much subdued until the well screen was mechanically rehabilitated in June 1996 . As soon as the hydraulic connection between the piezometer and the aquifer was restored, the temporal variability of the head in the aquifer became apparent.

## 6.5 Clock drift

Automated water level and pressure recorders use an autonomous or external clock which relies on crystal oscillators (commonly made of quartz) for a counter-measuring process that forms the basis for digital time keeping. The crystal oscillators are highly accurate, yet small deviations in their oscillation frequency, which changes with time (a phenomenon known as ageing), can add up over long measurement periods. This results in a gradual drift of the internal clock in relation to the real time. Elimination of this form of systematic error can be achieved by synchronising the transducer clock with a more accurate clock, such as in a field laptop that is always set to the same time zone and does not update to daylight savings time. Clock stability is an important consideration when using multiple instruments. Examples include the barometric correction of absolute pressure measurements from a non-vented transducer or the calculation of hydraulic gradients using two different time series.

Table 2Example of an assessment of clock stability for eight different standard PTs (AEDT is Australian eastern daylight time). None of the PTs complied with the clock stability of ±1 min yr−1, as specified by the manufacturer.

Figure 10The influence of clock stability on calculating a vertical head gradient: (a) pressure heads measured in a streambed using two PTs (note that the two lines are too close together to distinguish) with clocks that are in sync and the calculated vertical head gradient. (b) Erroneous vertical head gradients arising from time differences Δt due to out-of-sync instrument clocks caused by clock drift (the data were synthetically shifted).

The clock stability of eight PTs was assessed during a long-term surface water–groundwater exchange monitoring programme in the arid zone of Australia (Fowlers Creek at Fowlers Gap, New South Wales, Australia), where streams are dry for most parts of the year but flow if there is enough rainfall . Monitoring stream flow under such conditions relies on long-term and accurate resolution hydraulic gradients. To monitor the spatial and temporal dynamics of stream flow we used streambed arrays similar to those reported in . Before deploying the PTs, the internal clock of the field laptop was synchronised to an online time server. This ensured that all loggers had the same time stamp. The transducers were set up to start logging on 21 October 2014 at 18:00 AEDT (Australian eastern daylight time), with a sampling interval of 30 min.

Due to the remoteness of the field site, monitoring continued for over 2 years. After removal and disassembly of the streambed arrays, the internal clock of each PT was compared to that of a synchronised computer. The findings demonstrate that the majority of the PTs did not comply with the manufacturers' specifications of ±1 min yr−1, with most of the clocks running slower and the worst clock drift being +7.5 min yr−1 (Table 2). Such deviations are unfortunately not unusual for commercial PTs .

The influence of clock stability on measuring hydraulic head gradients is illustrated in Fig. 10. In this example, a vertical PT array similar to that used in was deployed in a streambed at Maules Creek (New South Wales). Resolving vertical head gradients over small distances is a significant challenge. Both PTs were calibrated against each other by placing the array inside a water bath overnight.

Figure 10a shows the pressure heads as well as the vertical head gradient during the experiment. Figure 10b illustrates the outcome if either one of the PTs was synchronised to a different time or as a result of clock drift. It is clear that the largest error arises during fastest head changes with time, where the gradient disregarding time errors could be interpreted as either gaining or losing conditions with different magnitudes. Similar to this example, showed how clock drift led to erroneous flow estimates in a coastal aquifer subject to ocean tides. Hydrological processes could be fundamentally misinterpreted if time-related monitoring errors are ignored, which is not always properly recognised.

## 6.6 Miscellaneous errors

Figure 11Water level in an observation well responding to irregular and frequent pumping from a nearby production bore automatically measured at two different sampling intervals (12 h or 30 min). Note that a manual dip which is taken at a time that falls between samples of automated measurement (e.g., the time period indicated by the grey box) can result in a significant measurement error in dynamic hydrogeological systems.

The importance of setting the logger to the appropriate time resolution (sampling rate) is illustrated by Fig. 11. Both lines show the same water level time series, but the red line shows how the curve would look if the measurement frequency had been set to twice daily, whereas the blue line shows the data as measured using an interval of 30 min. Obviously, the short-term variations caused by the operations of a nearby production bore are not captured when an inappropriate measurement interval is chosen. Similar issues can arise in aquifers affected by ocean tides or river stage fluctuations.

When unresolved temporal head fluctuations occur between two consecutive automated measurement intervals ti, a large discrepancy can arise between a manual measurement taken at time tj and the nearest measurement at time ti. When field personnel take manual measurements during their regular site visits, the timing of which is usually not determined by the logger recording settings but by logistical factors instead, considerable differences can arise from the fact that tjti. Using the data in Fig. 11 as an example, a manual head measurement taken at time (tj) of 09:11 on 5 September 2015 would be 0.83 m higher than the closest automated reading of the logger (set to a 12 h sampling interval) at the time (ti) of 12:00. The difference is unrelated to any instrument error and is solely due to unresolved temporal variability. A manual dip taken at any time between the two sampling times would result in an error that falls within the grey box in Fig. 11. While this hypothetical example represents an extreme case of this effect, misalignment of ti and tj is very common, and it supports the contention by that unrecognised hydrological processes are a form of measurement noise. To avoid such errors, a suitable measurement interval must be chosen upon initial logger deployment, which depends on the hydrogeological conditions at the measurement location. Only when it becomes clear that there is no temporal variability at this timescale can the sampling interval be increased to avoid unnecessary data handling and storage requirements.

Open GMI may suffer from a time delay in the water level response to changes in subsurface pore pressure, either because the well screen is partially clogged or improperly sized or because the well volume is so large that water cannot flow fast enough through the surrounding porous media and well screen to allow equilibration of the water level inside the well with the groundwater pressure (e.g., Hvorslev1951). In low-permeability materials like compacted peat or fine-grained sediments, time lags can be on the order of hours or even longer, which precludes the registration of the response to rapid processes such as for example river flooding events but also to water level changes induced by pumping, ocean tides or atmospheric pressure changes (e.g., Bredehoeft1967). Observation wells may also take appreciable time to readjust after the water level inside was raised by water displaced by inserting measurement instruments.

There are a variety of reasons why PTs do not always accurately record the water level in open GMI, many of which can be prevented by proper installation. When suspension cables are attached to well caps, the logger may not always be in exactly the same position after having been removed from the GMI. Some lightweight pressure transducers may experience buoyancy, especially in saltwater, and hence their vertical position is not constant in time. As a consequence of suspension cables being too short, PTs may end up being suspended above the water surface inside the GMI when the water level falls and hence record the atmospheric pressure . Air bubbles that become entrapped in the PT after the water level rises again can cause inaccurate readings and must be removed.

When open GMI becomes artesian, which can occur when water levels rise higher than the standpipe's top, the PT no longer indicates the true water level. When the PT is too deep to withstand the pressure of the water column (the so-called burst pressure, usually about twice the measurement range), the sensor may become damaged and the logger will malfunction. Freezing of the water column and lightning strikes can also cause damage to the PT . Sometimes PTs show erratic readings for no apparent reason, which can be due to the overheating of electronics. Temperatures in the standpipe sticking up above the land surface can easily exceed 40 C, the upper temperature threshold for correct functioning of electronic parts, due to sun exposure. Shading or ventilation measures are therefore also an important part of GMI.

One issue not discussed in the literature is the considerable confusion that arises due to clock adjustments related to daylight saving time (DST). Perhaps this is because it is considered a trivial point but an important one nonetheless that must be specifically addressed in the measurement protocol. Some devices automatically adjust to daylight savings time, whereas others do not. When they do, the instrument's clock setting depends on the time of year it was set up. Some manufacturers apply DST corrections only when the recorded data are exported to a file, and this depends on the computer's DST settings. The same data readings therefore even also end up with different time stamps if multiple computers are used.

7 Random error propagation

Figure 12(a) Visual comparison of horizontal and vertical random errors based on precision values in Table 1 (note that some errors are distance dependent) for the different steps (Fig. 1) and method options. Minimum achievable relative random head gradient error in the horizontal (b) and vertical (c) direction calculated using standard-practice measurements (highlighted by the black frame: step 1, options A, B and C; step 2, options B and C; step 3, option A; step 4, option A). Note that HHG and VHG values with errors exceeding 100 % are blanked out in (b) and (c). Please note the example error calculation in the main text.

Figure 1 illustrates steps required to calculate a head gradient. In what follows it is assumed that all systematic errors have been eliminated from this process in a way that the only error remaining is the random error. For the sake of simplicity and in the absence of further information, we also assume that all random errors are normally distributed and not correlated (Table 1). Furthermore, horizontal errors stated in Table 1 are isotropic, i.e. do not vary in the x and y directions. Note that all error quantification in this subsection assumes a standard practice that includes the following (see Table 1 and Fig. 12a):

1. Horizontal distance is measured using a total station (step 1, option B), and vertical distance is measured using digital levelling (step 1, option C). Distance errors are limited to the respective errors from DGNSS (step 1, option A). While GNSS surveying with a single receiver is useful for mapping, the precision of coordinates is not high enough to determine the distances between GMI for head gradient calculations.

2. Point of head is measured using a downhole camera for the vertical (step 2, option C) and verticality for horizontal error (step 2, option B), assuming a 10 m deep well.

3. The manual water level is measured using a dip meter (step 3, option A), and we assume that there is no depth dependency of the random error.

4. The automated pressure is measured using a vented transducer (step 4, option A).

Figure 12a graphically compares the random measurement error magnitudes for the different steps and methods summarised in Table 1. We stress that the adopted values reflect the absolute best-case scenario from current standard field practice.

The random error associated with head differences arises from steps 1, 3 and 4 and can be expressed as

$\begin{array}{}\text{(12)}& \mathit{\delta }\mathrm{\Delta }h=\mathrm{2}\sqrt{\left(\mathit{\delta }{z}_{g}{\right)}^{\mathrm{2}}+\left(\mathit{\delta }{d}_{\mathrm{w}}{\right)}^{\mathrm{2}}+\left(\mathit{\delta }{h}_{\mathrm{pt}}{\right)}^{\mathrm{2}}}.\end{array}$

We use the random errors associated with standard practice in this equation to estimate the minimum achievable precision (combined random error) when calculating head differences as δΔh=0.017 m. This error is somewhat higher than the findings by but lower than the field-based values (δΔh=0.022 m) reported by . Measured head differences that are smaller than this value will not allow much confidence in detecting the direction of groundwater flow. To improve this precision, approaches to reduce the achievable random errors when measuring steps 1, 3 and 4 must be found (Fig. 1), likely resulting in greater effort and a higher cost than what is currently standard practice.

Using the measurements illustrated in Fig. 1, the horizontal hydraulic head gradient (HHG) is calculated as

$\begin{array}{}\text{(13)}& \mathrm{\nabla }{h}^{\mathrm{h}}=\frac{{\left[{z}_{g}-{d}_{\mathrm{w}}+{h}_{\mathrm{pt}}\right]}_{\mathrm{2}}-{\left[{z}_{g}-{d}_{\mathrm{w}}+{h}_{\mathrm{pt}}\right]}_{\mathrm{1}}}{\mathrm{\Delta }{s}_{\mathrm{h}}^{\mathrm{h}}},\end{array}$

where the numeric subscripts depict the two locations. Analogously, the vertical hydraulic head gradient (VHG) can be determined as

$\begin{array}{}\text{(14)}& \mathrm{\nabla }{h}^{\mathrm{v}}=\frac{{\left[{z}_{g}-{d}_{\mathrm{w}}+{h}_{\mathrm{pt}}\right]}_{\mathrm{2}}-{\left[{z}_{g}-{d}_{\mathrm{w}}+{h}_{\mathrm{pt}}\right]}_{\mathrm{1}}}{{\left[{z}_{g}+\mathrm{\Delta }{z}_{p}\right]}_{\mathrm{2}}-{\left[{z}_{g}+\mathrm{\Delta }{z}_{p}\right]}_{\mathrm{1}}}.\end{array}$

A propagation of random errors accounts for the errors involved in measuring the different variables explained in Fig. 1. The relative error for the head gradients is

$\begin{array}{}\text{(15)}& \frac{\mathit{\delta }\mathrm{\nabla }{h}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}}{|\mathrm{\nabla }{h}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}|}=\sqrt{{\left(\frac{\mathit{\delta }\mathrm{\Delta }h}{\mathrm{\Delta }h}\right)}^{\mathrm{2}}+{\left(\frac{\mathit{\delta }\mathrm{\Delta }{s}_{\mathrm{h}}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}}{\mathrm{\Delta }{s}_{\mathrm{h}}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}}\right)}^{\mathrm{2}}},\end{array}$

where

$\begin{array}{}\text{(16)}& \mathit{\delta }\mathrm{\Delta }{s}_{\mathrm{h}}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}=\mathrm{2}\sqrt{{\left(\mathit{\delta }{s}_{g}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}\right)}^{\mathrm{2}}+{\left(\mathit{\delta }\mathrm{\Delta }{s}_{p}^{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}\right)}^{\mathrm{2}}}.\end{array}$

For the VHG case, $\mathit{\delta }\mathrm{\Delta }{s}_{\mathrm{h}}^{\mathrm{v}}=\mathit{\delta }\mathrm{\Delta }{z}_{\mathrm{h}}$, $\mathit{\delta }{s}_{g}^{\mathrm{v}}=\mathit{\delta }{z}_{g}$ and $\mathit{\delta }\mathrm{\Delta }{s}_{p}^{\mathrm{v}}=\mathit{\delta }\mathrm{\Delta }{z}_{p}$. These equations quantify the vertical positioning of the GMI and point of measurement resulting from an non-vertical borehole (steps 1 and 2). We use Eq. (15) to calculate the minimum achievable random relative error for HHGs and VHGs as a function of horizontal or vertical distance between two points of head measurement (Fig. 12b and c).

Figure 12 clearly demonstrates the relationship between HHGs or VHGs and distance between the points of head. In general, the greater the distance between screens, the smaller the relative head gradient error. For example, the random error of determining a HHG or VHG of 10−2 at a 10 m horizontal or vertical point of head distance is ≈17 % (see examples in Fig. 12). Figure 12b further illustrates that measuring a HHG $<{\mathrm{10}}^{-\mathrm{4}}$ with an error less than 100 % requires a distance 170 m between points of head. VHGs of 10−4 are unresolvable within the considered maximum vertical distance of 100 m (Fig. 12c). We stress that these errors are the best-case scenario, as in reality there is a likeliness of additional systematic errors contained in the measurements. The errors calculated here are thus unlikely to be achieved in practice. Extraordinary effort must be put towards improving the precision of measurements when head gradients less than 10−2 are to be calculated for distances smaller than 10 m (Fig. 12b). Note that in order to additionally determine the direction of the gradient, a minimum area between GMI is required . In practice, heterogeneity of the hydraulic conductivity will further add to the uncertainty of groundwater flow estimates.

8 Concluding remarks

Reliable water level measurements are at the core of every hydrogeological investigation, and the measurement error determines which processes or properties can be resolved. We have analysed unpublished and published data to quantify the best possible accuracy and precision of hydraulic head measurements using commonly available, state-of-the-art commercial instruments. By propagating the random errors, we find that with current standard practice, horizontal head gradients $<{\mathrm{10}}^{-\mathrm{4}}$ are only resolvable at distances 170 m and that it takes extraordinary effort to measure hydraulic head gradients $<{\mathrm{10}}^{-\mathrm{3}}$ over distances <10 m. However, we consider these estimates very optimistic, as they assume that systematic errors are absent or that systematic error corrections do not introduce additional error.

The magnitude of systematic errors tends to be much larger than that of random errors, and hence failure to recognise systematic errors can seriously compromise the outcomes of an investigation. It is difficult to establish if systematic errors are accumulating or cancelling, and hence they must be either avoided or identified and corrected. In part, systematic errors are due to the measurement conditions in the field, which are not easy to control and negatively affect instrument performance. But other factors play a role too, including improper instrument use, faulty or unsuitable GMI (e.g., long well screens), and the lack of measurement protocols that pay due consideration to all sources of error. Some measurement techniques have not seen performance improvement in decades, and there does not seem to be the same quest for measurement error reduction in hydrogeology as there is in other fields of science, where the smallest of dimensions are measured with ever-better accuracy and precision and advances in measurement technology are pushing the frontiers of science.

We acknowledge that the measurement error with available technology could already be sufficiently small (e.g., within a few centimetres) for a lot of practical applications. However, the quantification and reporting of measurement error does not seem to be commonplace yet. Existing standards like those of or contain useful guidelines for a maximum error as follows: (1) ±3 mm (0.01 ft) for general applications, (2) 0.1 % of expected water level changes and (3) 0.01 % for cases where the depth to water exceeds ≈30 m (100 ft). Such standards must see wider uptake, and the development of more sophisticated or site-specific standards, suited for a particular study area or research objective, would be important steps towards better hydrogeological data quality and consistency. Moreover, technological advances are necessary for enabling the measurement of vertical flow within an aquifer, the subtle temporal head fluctuations related to tidal cycles. Increased sensor performance and sensitivity would underpin new developments, such as the use of the groundwater response to Earth and atmospheric tides to characterise the degree of groundwater confinement (e.g. Butler et al.2011; Acworth et al.2017) and quantify compressible subsurface properties . Such advances highlight the need to innovate beyond standard practice to support research in the hydrogeological sciences. We believe that researchers and industry should work together and find ways to increase instrument performance.

The following list of recommendations synthesises the findings from our study and focuses on aspects that could considerably improve the current practice of hydraulic head measurement. These are as follows.

• Elimination of systematic errors. Our estimation of the minimum achievable random error across all measurements presumes the absence of systematic errors (Table  1; Fig. 12). Not all systematic errors (e.g., sensor drift) can be eliminated, but to minimise human error, measurements should be conducted exclusively by personnel that has received formal training. Moreover, a detailed measurement protocol must be designed and periodically evaluated, which outlines the procedures for measurement, maintenance, note keeping (using standardised field data sheets), and data storage and handling.

• Point of measurement. Our review of the literature demonstrates that GMI can significantly deviate from the vertical (Sect. 4). The resulting error in the point of head measurement is larger in the horizontal compared to the vertical direction (step 2 in Fig. 1). While this potentially introduces one of the largest errors, it is generally ignored when calculating head gradients (Table 1 and Fig. 12). If investigations necessitate the detection of small HHGs, we recommend measuring the borehole verticality using downhole profiling tools in open GMI. The best possible precision in measuring the point of head is achieved by combining a verticality sonde with a downhole optical camera. Geophysical logging (Keys2017) or flow meter measurements can identify GMI construction errors and ageing issues, such as casing leaks. For open GMI, joint interpretation of barometric pressure and water level time series is required to determine hydraulic heads (barometric correction; Sect. 2.3). For closed GMI, the point of head measurement accuracy depends on the details contained in the original drilling report if inclinometer casing is not used in the installation .

• Geo-spatial positioning. Our error propagation analysis (Fig. 12) clearly demonstrates that precise measurement of the horizontal and vertical distances between GMI is paramount to resolve the small hydraulic gradients inherent to groundwater investigations. This is particularly important when the GMI locations are in close proximity. Geo-spatial data from the single-receiver GNSS should only be used for mapping but not to calculate distances. Traditional surveying techniques deliver more precise results for horizontal distances <700 m compared to DGNSS. Vertical distances should only be calculated using data from digital levelling and not DGNSS. If possible, leap-frogging the survey device should be avoided (Sect. 3).

• Automated head measurements. The widespread use of automated PTs for hydraulic head and gradient measurement has perhaps led to the impression that manual measurements have become less important. Our analysis demonstrates that regular, frequent manual water level measurement remains essential, as it is the only way to verify that PTs are accurately recording the correct water level (Sect. 6.4.3). It also improves the precision of the automated measurements by averaging out the error introduced from manual dipping (Table 1). Given the significance of manual measurement, it is surprising to note that commercially available dip meters show as much error today (Fig. 4) as a quarter-century ago (Plazak1994). Telemetry does not obviate field site visits but only offers the convenience of not having to download devices and the advantage of being able to detect potential problems remotely, albeit at a higher cost of installation and maintenance (e.g., data service and connection problems).

• Time-related errors. Automated transducers rely on the stability of their internal or external time base once synchronised with the clock of the device that is used to set up the logging protocol. We demonstrate that clocks can drift significantly (Sect. 6.5), which leads to silent measurement errors and false interpretations (Fig. 10), especially for highly dynamic systems where uninterrupted long-term monitoring is required. When non-vented PTs are used to assess barometric effects, the clock stability error is a function of two clocks. We recommend that the clock is re-synchronised as frequently as possible or, where this is impossible, careful documentation of the device's internal clock status when monitoring is finished. Such practice is not always supported by off-the-shelf devices, and the limitations of the software and the instrument have to be trialled before deployment. Good time-keeping practice also includes the use of one and the same field laptop, which is regularly synchronised with a time server. Moreover, we recommend the use of an absolute time base, for example universal time coordinated (UTC), to avoid systematic errors arising from daylight savings time confusion.

• Density and temperature effects. We demonstrate that automated pressure measurement and accurate conversion into water levels necessitates knowledge of the average density of water inside the borehole (${\stackrel{\mathrm{‾}}{\mathit{\rho }}}_{\mathrm{w}}$; Eq. 11). Since water density depends on the amount of dissolved substances as well as temperature, there is a need for measuring water temperature and electric conductivity across the length of the water column to establish their potential influence. If the water density inside the open GMI is not constant, the best solution is to position the PT such that it measures ppt at elevation zh (Method 2 in Sect. 2.1), although this may come at the cost of greater measurement error due to a larger PT range. Further, PT readings are often affected by temperature despite internal compensation (Sect. 6.4.1). While subsurface temperatures beyond 2–3 m depth are generally roughly stable, avoiding temperature effects can be a significant problem when measuring water levels or barometric pressure at or near the surface.

• Type of pressure transducer. The type of PT is an important consideration that should be made according to the purpose of the investigation. For general groundwater monitoring away from topographic depressions and waterways (no risk of borehole over-topping), we recommend vented PTs as long as there is no problem with keeping the venting tube dry. Because vented PTs measure a relative pressure instead of an absolute pressure, they have a smaller range and do not require a separate instrument to simultaneously record the atmospheric pressure. As such, they have better accuracy, precision and resolution. Also, there is less risk of human error, and their use avoids the problems that are introduced with two PTs (i.e., sensor and clock drift). Nevertheless, barometric pressure must still be acquired in order to perform a barometric correction. For reliably resolving head gradients and flow direction at small vertical distances, for example when assessing surface water–groundwater interactions, we recommend the use of wet–wet differential pressure sensors (e.g., Cuthbert et al.2011). Quartz oscillator PTs are much more accurate than the commonly used strain gauge-type PT. However, they have hardly been used in groundwater studies to date, probably because of their higher cost.

• Technical specifications. In the technical specification of PTs, much of the focus is on accuracy as a percentage of the full-scale range. We noted that this value is not consistently defined between manufacturers and may contain adjustable (i.e., errors that can be corrected using manual depth-to-water measurements) as well as non-adjustable errors (hysteresis, repeatability and non-linearity). Before purchasing, it is wise to approach manufacturers and enquire about the various technical details. Consideration must be paid to minimising the measurement range in favour of maximum possible resolution (Sect. 2.4). Practice has shown that PTs have high failure rates; hence reliability is also an important selection criterion.

As a final remark, documentation of measurement procedure is critical for data validation. Without any assessment of the measurement uncertainty it is impossible to assign a quality label to the data, which severely limits their worthiness for consideration in public databases. In addition to data-collection protocols, quality-control procedures must be in place to ensure the reliability of the distributed water-level data. The development of such procedures should be considered in future work.

Data availability
Data availability.

Some of the ground- water data are available through the NCRIS Groundwater Database: http://groundwater.anu.edu.au (last access: 2 September 2019). Additional data are available at https://doi.org/10.13140/RG.2.2.28337.94569 (Shanafield et al., 2019).

Author contributions
Author contributions.

The idea for this paper was conceived by GCR and VEAP, and they were the lead authors. Everyone contributed ideas, data and writing. GCR coordinated the effort and produced the figures and tables.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

The authors thank Nick White for collecting the data shown in Fig. 4 and Michael Teubner for the measurements shown in Fig. 6.

Financial support
Financial support.

This research was supported by the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant (no. 835852). Gabriel C. Rau acknowledges support from the Australian NSW State Government’s Research Attraction and Acceleration Program in 2019. Margaret Shanafield was supported by the Australian Research Council (DE150100302). Eddie W. Banks acknowledges support from the Australian Research Council (ARC) Linkage Project (LP150100588). We acknowledge support by the publication fund of the Karlsruhe Institute of Technology (KIT). Some of the data used in this analysis were collected with equipment provided by the National Collaborative Research Infrastructure Strategy (NCRIS) financed by the Australian federal government.

The article processing charges for this open-access
publication were covered by a Research
Centre of the Helmholtz Association.

Review statement
Review statement.

This paper was edited by Brian Berkowitz and reviewed by Donald Rosenberry and one anonymous referee.

References

Acworth, R. I., Rau, G. C., McCallum, A. M., Andersen, M. S., and Cuthbert, M. O.: Understanding connected surface-water/groundwater systems using Fourier analysis of daily and sub-daily head fluctuations, Hydrogeol. J., 23, 143–159, https://doi.org/10.1007/s10040-014-1182-5, 2015. a

Acworth, R. I., Halloran, L. J. S., Rau, G. C., Cuthbert, M. O., and Bernardi, T. L.: An objective frequency domain method for quantifying confined aquifer compressible storage using Earth and atmospheric tides, Geophys. Res. Lett., 43, 611–671, https://doi.org/10.1002/2016GL071328, 2016a. a, b

Acworth, R. I., Rau, G. C., Cuthbert, M. O., Jensen, E., and Leggett, K.: Long-term spatio-temporal precipitation variability in arid-zone Australia and implications for groundwater recharge, Hydrogeol. J., 24, 905–921, https://doi.org/10.1007/s10040-015-1358-7, 2016b. a

Acworth, R. I., Rau, G. C., Halloran, L. J. S., and Timms, W. A.: Vertical groundwater storage properties and changes in confinement determined using hydraulic head response to atmospheric tides, Water Resour. Res., 53, 2983–2997, https://doi.org/10.1002/2016WR020311, 2017. a, b

Atwood, D. and Lamb, B.: Resolution problems with obtaining accurate ground water elevation measurement in a hydrogeologic site investigation, in: Proceedings, First National Outdoor Action Conference on Aquifer Restoration, Ground Water Monitoring, and Geophysical Methods, National Water Well Association, Westerville, OH, 185–193, 1987. a, b

Bailey, D.: Preface, in: Practical Radio Engineering and Telemetry for Industry, edited by: Bailey, D., p. xi, Elsevier, Oxford, https://doi.org/10.1016/B978-075065803-4/50013-5, 2003. a

Bakuła, M., Oszczak, S., and Pelc-Mieczkowska, R.: Performance of RTK Positioning in Forest Conditions: Case Study, J. Surv. Eng., 135, 125–130, https://doi.org/10.1061/(ASCE)0733-9453(2009)135:3(125), 2009. a

BDA: Guidance For The Operation Of Cable Percussion Rigs And Equipment (rev 1.4), Tech. rep., British Drilling Association (BDA), Pinxton, 2017. a

Benjamin, J. and Kaplan, D.: Development of a laser-based water level sensor for fine-scale ecohydrological measurements, in: 2017 IEEE Conference on Technologies for Sustainability (SusTech), 1–8, IEEE, 12–14 November 2017, Phoenix, AZ, USA, https://doi.org/10.1109/SusTech.2017.8333522, 2017. a, b, c

Beshr, A. A. and Abo Elnaga, I. M.: Investigating the accuracy of digital levels and reflectorless total stations for purposes of geodetic engineering, Alexandria Engineering Journal, 50, 399–405, https://doi.org/10.1016/j.aej.2011.12.004, 2011. a

Bitelli, G., Roncari, G., Tini, M. A., and Vittuari, L.: High-precision topographical methodology for determining height differences when crossing impassable areas, Measurement: Journal of the International Measurement Confederation, 118, 147–155, https://doi.org/10.1016/j.measurement.2018.01.013, 2018. a, b, c

Bock, Y. and Melgar, D.: Physical applications of GPS geodesy: A review, Rep. Prog. Phys., 79, 106801, https://doi.org/10.1088/0034-4885/79/10/106801, 2016. a, b

Bouma, J., Maasbommel, M., and Schuurman, I.: Handboek meten van grondwaterstanden in peilbuizen (Translated title: Handbook on measurement of groundwater levels in piezometers), Rapport/STOWA; 2012-50, STOWA, Amersfoort, 2012 (in Dutch). a

Braun, J., Štroner, M., Urban, R., and Dvořáček, F.: Suppression of Systematic Errors of Electronic Distance Meters for Measurement of Short Distances, Sensors, 15, 19264–19301, https://doi.org/10.3390/s150819264, 2015. a, b, c

Bredehoeft, J. D.: Response of well-aquifer systems to Earth tides, J. Geophys. Res., 72, 3075–3087, https://doi.org/10.1029/JZ072i012p03075, 1967. a

Brinker, R. C.: The Surveying Handbook, Springer US, Boston, MA, https://doi.org/10.1097/BCR.0b013e3182a2ad17, 1995. a, b, c

Buchroithner, M. F. and Pfahlbusch, R.: Geodetic grids in authoritative maps–new findings about the origin of the UTM Grid, Cartogr. Geogr. Info. Sci., 44, 186–200, https://doi.org/10.1080/15230406.2015.1128851, 2017. a

Bulant, P., Eisner, L., Pšenčík, I., and Calvez, J. L.: Importance of borehole deviation surveys for monitoring of hydraulic fracturing treatments, Geophys. Prospect., 55, 891–899, https://doi.org/10.1111/j.1365-2478.2007.00654.x, 2007. a, b

Butler, J. J., Jin, W., Mohammed, G. A., and Reboulet, E. C.: New insights from well responses to fluctuations in barometric pressure, Ground Water, 49, 525–533, https://doi.org/10.1111/j.1745-6584.2010.00768.x, 2011. a, b

Cain, S. F., Davis, G. A., Loheide, S. P., and Butler, J. J.: Noise in Pressure Transducer Readings Produced by Variations in Solar Radiation, Ground Water, 42, 939–944, https://doi.org/10.1111/j.1745-6584.2004.t01-12-.x, 2004. a

Ceylan, A. and Baykal, O.: Precise Height Determination Using Leap-Frog Trigonometric Leveling, J. Surv. Eng., 132, 118–123, https://doi.org/10.1061/(ASCE)0733-9453(2006)132:3(118), 2006. a, b

Clark, W. E.: Computing the barometric efficiency of a well, J. Hydr. Eng. Div.-ASCE, 93, 93–98, 1967. a

Contreras, I. A., Grosser, A. T., and Ver Strate, R. H.: The Use of the Fully-Grouted Method for Piezometer Installation, in: 7th FMGM 2007, 1–20, American Society of Civil Engineers, Reston, VA, https://doi.org/10.1061/40940(307)67, 2007. a

Cuevas, J., Calvo, M., Little, C., Pino, M., and Dassori, P.: Are diurnal fluctuations in streamflow real?, J. Hydrol. Hydromech., 58, 149–162, https://doi.org/10.2478/v10098-010-0014-0, 2010. a, b

Cunningham, W. L. and Schalk, C. W.: Groundwater Technical Procedures of the U. S. Geological Survey Techniques and Methods 1 – A1, Tech. rep., U.S. Geological Survey, Reston, Virginia, USA, 2016. a, b, c, d, e

Cuthbert, M., Greswell, R., and Mackay, R.: A Wet/Wet Differential Pressure Sensor for Measuring Vertical Hydraulic Gradient, Ground Water, 49, 781–782, https://doi.org/10.1111/j.1745-6584.2011.00789.x, 2011. a, b

Devlin, J. F. and McElwee, C. D.: Effects of Measurement Error on Horizontal Hydraulic Gradient Estimates, Ground Water, 45, 62–73, https://doi.org/10.1111/j.1745-6584.2006.00249.x, 2007. a, b, c, d, e, f

Domenico, P. A. and Schwartz, F. W.: Physical and Chemical Hydrogeology, John Wiley & Sons, Inc., New York, 2nd edn., 1997. a, b, c, d

Drexler, J. Z., Bedford, B. L., Scognamiglio, R., and Siegel, D. I.: Fine-scale characteristics of groundwater flow in a peatland, Hydrol. Process., 13, 1341–1359, https://doi.org/10.1002/(SICI)1099-1085(19990630)13:9<1341::AID-HYP810>3.0.CO;2-5, 1999. a

Dunnicliff, J. and Green, G. E.: Geotechnical Instrumentation for Monitoring Field Performance, A Wiley-Interscience publication, Wiley, New York, 1993. a

El-Ashmawy, K. L. A.: Accuracy, time cost and terrain independence comparisons of levelling techniques, Geodesy and Cartography, 40, 133–141, https://doi.org/10.3846/20296991.2014.962727, 2014. a

Elçi, A., Flach, G. P., and Molz, F. J.: Detrimental effects of natural vertical head gradients on chemical and water level measurements in observation wells: identification and control, J. Hydrol., 281, 70–81, https://doi.org/10.1016/S0022-1694(03)00201-4, 2003. a

Freeman, L. A., Carpenter, M. C., Rosenberry, D. O., Rousseau, J. P., Unger, R., and McLean, J. S.: Use of Submersible Pressure Transducers in Water-Resources Investigations, Tech. rep., U.S. Geological Survey, Reston, Virginia, USA, 2004. a, b, c, d

Freeze, R. A. and Cherry, J. A.: Groundwater, Prentice Hall, Inc., Upper Saddle River, NJ 07458, 1979. a, b, c, d, e

Garrido, M. S., Giménez, E., de Lacy, M. C., and Gil, A. J.: Surveying at the limits of local RTK networks: Test results from the perspective of high accuracy users, Int. J. Appl. Earth Obs., 13, 256–264, https://doi.org/10.1016/j.jag.2010.12.001, 2011. a, b

Gonthier, G.: A Graphical Method for Estimation of Barometric Efficiency from Continuous Data – Concepts and Application to a Site in the Piedmont, Air Force Plant 6, Marietta, Georgia, Tech. rep., US Geological Survey, Reston, Virginia, USA, https://doi.org/10.3133/sir20075111, 2003. a

Gribovszki, Z., Kalicz, P., and Szilágyi, J.: Does the accuracy of fine-scale water level measurements by vented pressure transducers permit for diurnal evapotranspiration estimation?, J. Hydrol., 488, 166–169, https://doi.org/10.1016/J.JHYDROL.2013.03.001, 2013. a, b

Hanschke, T. and Baird, A. J.: Time-lag errors associated with the use of simple standpipe piezometers in wetland soils, Wetlands, 21, 412–421, https://doi.org/10.1672/0277-5212(2001)021[0412:TLEAWT]2.0.CO;2, 2006. a

Hegarty, C. J.: The Global Positioning System (GPS), in: Springer Handbook of Global Navigation Satellite Systems, 197–218, Springer International Publishing, Cham, https://doi.org/10.1007/978-3-319-42928-1_7, 2017. a, b, c

Hölting, B. and Coldewey, W. G.: Hydrogeologie, Spektrum Akademischer Verlag, Heidelberg, https://doi.org/10.1007/978-3-8274-2354-2, 2013. a, b

Hubbell, J. M., Sisson, J. B., Nicholl, M. J., and Taylor, R. G.: Well Design to Reduce Barometric Pressure Effects on Water Level Data in Unconfined Aquifers, Vadose Zone J., 3, 183–189, https://doi.org/10.2113/3.1.183, 2004. a, b

Hubbert, M. K.: The theory of ground-water motion, Eos, Transactions American Geophysical Union, 21, 648–648, https://doi.org/10.1029/TR021i002p00648-1, 1940. a

Hvorslev, M.: Time Lag and Soil Permeability in Ground-water Observations, Tech. rep., Army Engineer Waterways Experiment Station, Vicksburg, Mississippi, USA, 1951. a

Jacob, C. E.: On the flow of water in an elastic artesian aquifer, Eos, Transactions American Geophysical Union, 21, 574–586, 1940. a

Kalbus, E., Reinstorf, F., and Schirmer, M.: Measuring methods for groundwater – surface water interactions: a review, Hydrol. Earth Syst. Sci., 10, 873–887, https://doi.org/10.5194/hess-10-873-2006, 2006. a

Keys, W. S.: A Practical Guide to Borehole Geophysics in Environmental Investigations, Routledge, Abingdon, https://doi.org/10.1201/9781315136417, 2017. a

Kim Sun, G. O. and Gibbings, P.: How well does the virtual reference station (VRS) system of gps base stations perform in comparison to conventional RTK?, J. Spat. Sci., 50, 59–73, https://doi.org/10.1080/14498596.2005.9635038, 2005. a

Knotters, M., Meij, T. v. d., and Pleijter, M.: Nauwkeurigheid van handmatig gemeten grondwaterstanden en stijghoogtes: verslag van een veldexperiment, Wageningen, Tech. rep., 2013. a, b, c, d, e

Kouba, J., Lahaye, F., and Tétreault, P.: Precise Point Positioning, in: Springer Handbook of Global Navigation Satellite Systems, chap. E(25), 724–751, Springer International Publishing, Basel, https://doi.org/10.1007/978-3-319-42928-1_7, 2017. a

Li, X., Ge, M., Dai, X., Ren, X., Fritsche, M., Wickert, J., and Schuh, H.: Accuracy and reliability of multi-GNSS real-time precise positioning: GPS, GLONASS, BeiDou, and Galileo, J. Geodesy., 89, 607–635, https://doi.org/10.1007/s00190-015-0802-8, 2015. a

Liu, Z. and Higgins, C. W.: Does temperature affect the accuracy of vented pressure transducer in fine-scale water level measurement?, Geosci. Instrum. Method. Data Syst., 4, 65–73, https://doi.org/10.5194/gi-4-65-2015, 2015. a, b, c, d

Lusczynski, N. J.: Head and flow of ground water of variable density, J. Geophys. Res., 66, 4247–4256, https://doi.org/10.1029/JZ066i012p04247, 1961. a

Mäkinen, R. and Orvomaa, M.: Experiences and recommendations on automated groundwater monitoring, in: 20th International Northern Research Basins Symposium and Workshop – Kuusamo, Finland, 16–21 August 2015, edited by: Korhonen, J. and Kuusisto, E., 71–75, Finnish Environment Institute (SYKE), Kuusamo, Finland, 2015. a, b

McCallum, A. M., Andersen, M. S., Rau, G. C., Larsen, J. R., and Acworth, R. I.: River-aquifer interactions in a semiarid environment investigated using point and reach measurements, Water Resour. Res., 50, 2815–2829, https://doi.org/10.1002/2012WR012922, 2014. a, b, c

McKenna, G. T.: Grouted-in installation of piezometers in boreholes, Can. Geotech. J., 32, 355–363, https://doi.org/10.1139/t95-035, 1995. a, b

McLaughlin, D. L. and Cohen, M. J.: Thermal artifacts in measurements of fine-scale water level variation, Water Resour. Res., 47, https://doi.org/10.1029/2010WR010288, 2011. a, b

McMillan, T. C., Rau, G. C., Timms, W. A., and Andersen, M. S.: Utilizing the Impact of Earth and Atmospheric Tides on Groundwater Systems: A Review Reveals the Future Potential, Rev. Geophys., 57, 2018RG000630, https://doi.org/10.1029/2018RG000630, 2019. a, b

Meinzer, O. E.: Ground water in the United States, a summary of ground-water conditions and resources, utilization of water from wells and springs, methods of scientific investigation, and literature relating to the subject, Tech. rep., U.S. G.P.O., Reston, Virginia, USA, https://doi.org/10.3133/wsp836D, 1939. a

Mikkelsen, P. E. and Green, G. E.: “Piezometers in fully grouted boreholes”, in:, Proceedings of the Sixth International Symposium on Field Measurements in Geomechanics, edited by: Myrvoll, F., CRC Press, London, pp. 545-553, 2003. a, b

Misra, P. and Enge, P.: Global Positioning System: Signals, Measurements, and Performance, Revised Second Edition, Ganga-Jamuna Press, Lincoln, 2010. a

Morgenschweis, G.: Messung des Wasserstands (Translated title: Measurement of water level), in: Hydrometrie (Translated title: Hydrometry), edited by: Morgenschweis, G., 25–114, Springer Berlin Heidelberg, Berlin, Heidelberg, https://doi.org/10.1007/978-3-662-55314-5_3, 2018 (in German). a, b

Nielsen, D. M. and Nielsen, G.: The Essential Handbook of Ground-Water Sampling, Cambridge University Press, Cambridge, 2006. a, b, c

Nilsson, G. and Nissen, J.: Revision of borehole deviation measurements in Forsmark. Forsmark site investigation, Tech. rep., Swedish Nuclear Fuel and Waste Management Company, Stockholm, 2007. a

Noorduijn, S. L., Cook, P. G., Wood, C., and White, N.: Using Sealed Wells to Measure Water Levels Beneath Streams and Floodplains, Groundwater, 53, 872–876, https://doi.org/10.1111/gwat.12287, 2015. a, b

NUDLC: Minimum Construction Requirements for Water Bores in Australia, Tech. rep., National Uniform Drillers Licensing Committee, Australian Government National Water Commission, Canberra, 2011. a, b

Ohmer, M., Liesch, T., Goeppert, N., and Goldscheider, N.: On the optimal selection of interpolation methods for groundwater contouring: An example of propagation of uncertainty regarding inter-aquifer exchange, Adv. Water Resour., 109, 121–132, https://doi.org/10.1016/j.advwatres.2017.08.016, 2017. a

Plazak, D.: Differences Between Water- Level Probes, Ground Water Monit. R., 14, 84, https://doi.org/10.1111/j.1745-6592.1994.tb00091.x, 1994. a, b, c

Pleijter, M., Hamersveld, L., and Knotters, M.: Systematische fouten in metingen van grondwaterstanden met drukopnemers: verslag van een data-analyse (Alterra-rapport: 2666), 1775, Alterra, Wageningen-UR, available at: http://edepot.wur.nl/360149 (last access: 2 September 2019), 2015. a, b

Post, V. E. A. and von Asmuth, J. R.: Review: Hydraulic head measurements–new technologies, classic pitfalls, Hydrogeol. J., 21, 737–750, https://doi.org/10.1007/s10040-013-0969-0, 2013. a, b, c, d, e, f

Post, V. E. A., Banks, E., and Brunke, M.: Groundwater flow in the transition zone between freshwater and saltwater: a field-based study and analysis of measurement errors, Hydrogeol. J., 26, 1821–1838, https://doi.org/10.1007/s10040-018-1725-2, 2018. a, b, c, d

Rasmussen, T. C. and Crawford, L. A.: Identifying and Removing Barometric Pressure Effects in Confined and Unconfined Aquifers, Ground Water, 35, 502–511, https://doi.org/10.1111/j.1745-6584.1997.tb00111.x, 1997. a, b

Rau, G. C., Acworth, R. I., Halloran, L. J. S., Timms, W. A., and Cuthbert, M. O.: Quantifying Compressible Groundwater Storage by Combining Cross-Hole Seismic Surveys and Head Response to Atmospheric Tides, J. Geophys. Res.-Earth, 123, 1910–1930, https://doi.org/10.1029/2018JF004660, 2018. a, b

Remondi, B. W.: Performing Centimeter-Level Surveys in Seconds with GPS Carrier Phase: Initial Results, Navigation, 32, 386–400, https://doi.org/10.1002/j.2161-4296.1985.tb00918.x, 1985. a

Ritzema, H. P., Heuvelink, G. B. M., Heinen, M., Bogaart, P. W., Bolt, F. J. E. v. d., Hack-ten Broeke, M. J. D., Hoogland, T., Knotters, M., Massop, H. T. L., and Vroon, H. R. J.: Meten en interpreteren van grondwaterstanden: analyse van methodieken en nauwkeurigheid, Tech. rep., Alterra, Wageningen-UR, 309, 2012. a

Rizos, C.: Surveying, in: Springer Handbook of Global Navigation Satellite Systems, chap. F(35), 1011–1038, Springer International Publishing, Basel, https://doi.org/10.1007/978-3-319-42928-1_7, 2017. a

Rosenberry, D. O.: Effect of Sensor Error on Interpretation of Long-Term Water-Level Data, Ground Water, 28, 927–936, https://doi.org/10.1111/j.1745-6584.1990.tb01729.x, 1990. a, b

Rosenberry, D. O., LaBaugh, J. W., and Hunt, R. J.: Use of Monitoring Wells, Portable Piezometers, and Seepage Meters to Quantify Flow Between Surface Water and Ground Water, in: Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water, chap. 2 of Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water, U.S. Geological Survey, Reston, Virginia, USA, 39–70, 2008. a, b

Rowe, R. K. and Nadarajah, P.: Evaluation of the hydraulic conductivity of aquitards, Int. J. Rock. Mech. Min., 31, 781–800, https://doi.org/10.1016/0148-9062(94)90614-9, 1994. a

Saines, M.: Errors in interpretation of ground-water level data, Ground Water Monit. R., 1, 56–61, https://doi.org/10.1111/j.1745-6592.1981.tb00798.x, 1981. a

Shanafield, M. and Cook, P. G.: Transmission losses, infiltration and groundwater recharge through ephemeral and intermittent streambeds: A review of applied methods, J. Hydrol., 511, 518–529, https://doi.org/10.1016/j.jhydrol.2014.01.068, 2014. a

Shanafield, M., Post, E. A. V., Rau, G., Banks, E., Blum, P., and Krekeler, T.: Fig5a rawdata, https://doi.org/10.13140/RG.2.2.28337.94569, 2019.

Siejka, Z.: Validation of the Accuracy and Convergence Time of Real Time Kinematic Results Using a Single Galileo Navigation System, Sensors, 18, 2412, https://doi.org/10.3390/s18082412, 2018. a

Silliman, S. E. and Mantz, G.: The Effect of Measurement Error on Estimating the Hydraulic Gradient in Three Dimensions, Ground Water, 38, 114–120, https://doi.org/10.1111/j.1745-6584.2000.tb00208.x, 2000. a, b, c, d

Simeoni, L.: Laboratory tests for measuring the time-lag of fully grouted piezometers, J. Hydrol., 438–439, 215–222, https://doi.org/10.1016/j.jhydrol.2012.03.025, 2012. a, b

Sokol, D.: Position and fluctuations of water level in wells perforated in more than one aquifer, J. Geophys. Res., 68, 1079–1080, https://doi.org/10.1029/JZ068i004p01079, 1963. a

Sorensen, J. P. and Butcher, A. S.: Water Level Monitoring Pressure Transducers-A Need for Industry-Wide Standards, Ground Water Monit. R., 31, 56–62, https://doi.org/10.1111/j.1745-6592.2011.01346.x, 2011. a, b, c, d, e

Spane, F. A.: Considering barometric pressure in groundwater flow investigations, Water Resour. Res., 38, 14–1, https://doi.org/10.1029/2001WR000701, 2002. a, b

Spane, F. J. and Mercer, R.: HEADCO: a program for converting observed water levels and pressure measurements to formation pressure and standard hydraulic head, Tech. rep., Rockwell International Corporation, Richland, WA, USA, 1985. a

STS Sensors: Genauigkeitsangaben bei Drucksensoren richtig deuten (Translated title: Correct interpretation of accuracy specifications of pressure transducers), available at: https://campaign.stssensors.com/blog/de/druckmesstechnik-genauigkeitsangaben-drucksensor (24 January 2019), 2017 (in German). a

Sweet, H., Rosenthal, G., and Atwood, D.: Water level monitoring – Achievable accuracy and precision, in: Ground Water and Vadose Zone Monitoring, edited by Nielsen, D. and Johnson, A., ASTM STP 1053, 178–192, ASTM, 1990. a, b, c, d

Toll, N. J. and Rasmussen, T. C.: Removal of Barometric Pressure Effects and Earth Tides from Observed Water Levels, Ground Water, 45, 101–105, https://doi.org/10.1111/j.1745-6584.2006.00254.x, 2007. a

Treskatis, C.: Bewirtschaftung von Grundwasserressourcen: Planung, Bau und Betrieb von Grundwasserfassungen, Inst. WAR, Darmstadt, Germany, 2006. a

Twining, B. V. B.: Borehole deviation and correction factor data for selected wells in the eastern Snake River Plain aquifer at and near the Idaho National Laboratory, Idaho, Tech. rep., U.S. Geological Survey, Reston, Virginia, USA, https://doi.org/10.3133/sir20165163, 2016. a

van der Kamp, G. and Gale, J. E.: Theory of earth tide and barometric effects in porous formations with compressible grains, Water Resour. Res., 19, 538–544, https://doi.org/10.1029/WR019i002p00538, 1983.  a

Walker, J. and Awange, J. L.: Total Station Differential Levelling, in: Surveying for Civil and Mine Engineers, 93–101, Springer International Publishing, Cham, https://doi.org/10.1007/978-3-319-53129-8_6, 2018. a

Weeks, E. P.: Barometric fluctuations in wells tapping deep unconfined aquifers, Water Resour. Res., 15, 1167–1176, https://doi.org/10.1029/WR015i005p01167, 1979. a

Willemsen, J.: Basismeetnet grondwater (Translated title: Primary monitoring network groundwater), Tech. rep., Waterschap Hollandse Delta, Ridderkerk, the Netherlands 2006 (in Dutch). a, b

Yamano, M., Goto, S., Miyakoshi, A., Hamamoto, H., Lubis, R. F., Monyrath, V., and Taniguchi, M.: Reconstruction of the thermal environment evolution in urban areas from underground temperature distribution, Sci. Total Environ., 407, 3120–3128, https://doi.org/10.1016/j.scitotenv.2008.11.019, 2009. a

Zandbergen, P. A. and Barbeau, S. J.: Positional accuracy of assisted GPS data from high-sensitivity GPS-enabled mobile phones, J. Navigation, 64, 381–399, https://doi.org/10.1017/S0373463311000051, 2011. a

Zarriello, P. J.: Accuracy, Precision, and Stability of a Vibrating‐Wire Transducer Measurement System to Measure Hydraulic Head, Ground Water Monit. R., 15, 157–168, https://doi.org/10.1111/j.1745-6592.1995.tb00527.x, 1995. a, b