Changes in zooplankton communities along a mercury contamination gradient in a coastal lagoon (Ria de Aveiro,Portugal)

The main objective of this paper was to evaluate the impact of mercury on the zooplankton communities’ structure and functioning and their bioaccumulation patterns along a contamination gradient in a temperate coastal lagoon. Our results demonstrated that total abundance was not negatively affected by Hg contamination, since the most contaminated areas presented the highest values, being the copepod Acartia tonsa the dominant species, which means that it is a very well adapted and tolerant species to mercury. Nevertheless, negative effects were observed in terms of species diversity, since the most contaminated areas presented the lowest values of species richness, evenness and heterogeneity. Moreover, the spatial mercury gradient was reflected on the bioaccumulation patterns of the zooplankton communities. This reinforces the idea that zooplankton can be considered as an important vehicle of mercury transfer through the food pelagic web since it constitutes a primordial food resource for several commercial fish species.     

Long‐term change of daily and multi‐daily precipitation in southern Sweden

Daily rain series from southern Sweden with records dating back to the 1870s have been analysed to investigate the trends of daily and multi‐day precipitation of different return periods with emphasis on the extremes. Probabilities of extreme storms were determined as continuously changing values based on 25 years of data. An extra set of data was used to investigate changes in Skåne, the southernmost peninsula of Sweden. Another 30‐year data set of more than 200 stations of a dense gauge network in Skåne was used to investigate the relation between very large daily rainfall and annual precipitation. The annual precipitation has increased significantly all over southern Sweden due to increased winter precipitation. There is a trend of increasing maximum annual daily precipitation at only one station, where the annual maximum often occurs in winter. The number of events with a short return period is increasing, but the number of more extreme events has not increased. Daily and multi‐daily design storms of long return periods determined from extreme value analysis with updating year by year are not higher today than during the last 100 years. The largest daily storms are not related to stations with annual rainfall but seem to occur randomly. Copyright © 2013 John Wiley & Sons, Ltd.     

Rainfall,run‐off,and sediment transport dynamics in a humid mountain badland area: Long‐term results from a small catchment

In this study, we investigated rainfall, run‐off, and sediment transport dynamics (414 run‐off events and 231 events with sediment information) of a humid mountain badland area—the Araguás catchment (Central Pyrenees, Spain)—from October 2005 to September 2016. Use of this long‐term database allows characterization of the hydrological response, which consist of low‐magnitude/high‐frequency events and high‐magnitude/low‐frequency events, and identification of seasonal dynamics and rainfall‐run‐off thresholds. Our results indicate that the Araguás catchment, similarly to other humid badlands, had high hydrological responsiveness (mean annual run‐off coefficient: 0.52), a non‐linear relationship of rainfall with run‐off (common in Mediterranean environments), and seasonal hydrological and sedimentological dynamics. We created and validated a multivariate regression model to characterize the hydrological variables (stormflow and peak discharge) and sedimentological variables (mean and maximum suspended sediment concentrations and total suspended sediment load). In summer and at the beginning of autumn, the response was mainly related to rainfall intensity, suggesting a predomination of Hortonian flows. In contrast, in spring and winter, the responses were mainly related to the antecedent conditions (previous rainfall and baseflow), suggesting the occurrence of saturated excess flow processes, and the contribution of neighbouring vegetated areas. The multivariate analysis also showed that total sediment load is better predicted by a multivariate regression model that integrates pre‐event, rainfall, and run‐off variables. In general, our models provided more accurate predictions of small‐magnitude/high‐frequency events than high‐magnitude/low‐frequency events. This study highlights the high inter‐ and intra‐annual variability response in humid badland areas and that long‐term records are needed to reduce the uncertainty of hydrological and sedimentological responses in Mediterranean badland areas.     

Model simulations of tidal changes in a coastal lagoon,the Ria de Aveiro (Portugal)

Tidal changes in the Ria de Aveiro lagoon, Portugal, show an average increase of 0.245 m in M₂ amplitude and an average 17.4° decrease in M₂ phase, over 16 years (1987–2004). This paper investigates the causes of these changes using: (a) an analytical model and (b) hydrodynamic (vertically integrated) model of the Ria de Aveiro lagoon.     

Long‐term discharge prediction for the Turnu Severin station (the Danube) using a linear autoregressive model

The long‐term trends of yearly discharge time series and runoff variability at seven stations along the River Danube are identified. The results of statistical analysis of discharge time series indicate the period around the year 1860 was the driest decade in central and eastern Europe since 1840. In these years, the mean annual air temperature in central Europe was lower by about 1 °C compared with the 1990s. It is important to notice that the two driest decades (around 1860s and 1990s) of the instrumental era occurred in very different temperature conditions. The 28–31 years; 20–21 years; 14 years, as well as 4·2, 3·6, and 2·4 years fluctuations of annual discharge in the River Danube were found. Also, the long‐term streamflow prediction based on stochastic modelling methods is treated. Harmonic models and the Box–Jenkins methods were used. The predictions of yearly River Danube discharge time series were made for two decades ahead. From the stochastic models it follows that the annual discharge in the Danube at Turnu Severin station should reach its local maximum within the years 2004–06. The period 2015–19 should be dry. Copyright © 2005 John Wiley & Sons, Ltd.     

Long‐term investigations of the snow cover in a subalpine semi‐forested catchment

To improve spring runoff forecasts from subalpine catchments, detailed spatial simulations of the snow cover in this landscape is obligatory. For more than 30 years, the Swiss Federal Research Institute WSL has been conducting extensive snow cover observations in the subalpine watershed Alptal (central Switzerland). This paper summarizes the conclusions from past snow studies in the Alptal valley and presents an analysis of 14 snow courses located at different exposures and altitudes, partly in open areas and partly in forest. The long‐term performance of a physically based numerical snow–vegetation–atmosphere model (COUP) was tested with these snow‐course measurements. One single parameter set with meteorological input variables corrected to the prevailing local conditions resulted in a convincing snow water equivalent (SWE) simulation at most sites and for various winters with a wide range of snow conditions. The snow interception approach used in this study was able to explain the forest effect on the SWE as observed on paired snow courses. Finally, we demonstrated for a meadow and a forest site that a successful simulation of the snowpack yields appropriate melt rates. Copyright © 2006 John Wiley & Sons, Ltd.     

Seasonal Liza aurata tissue-specific DNA integrity in a multi-contaminated coastal lagoon (Ria de Aveiro, Portugal)

In this study, the DNA integrity of golden grey mullet (Liza aurata) collected in differently contaminated sites of a coastal lagoon, Ria de Aveiro (Portugal), was assessed, over the period of 1 year, using the DNA alkaline unwinding assay, in four different tissues (gill, kidney, liver and blood) and compared to a reference site. The four tissues displayed different DNA integrity basal levels, clearly affected by seasonal factors. Gill and kidney were, respectively, the most and least sensitive tissues. All sites demonstrated the capacity to interfere with DNA integrity. The sites displaying the highest and lowest DNA damage capability were, respectively, Barra (subject to naval traffic) and Vagos (contaminated with polycyclic aromatic hydrocarbons). In terms of seasonal variability, autumn seems to be the more critical season (more DNA damage) unlike summer when no DNA damage was found in any tissue. Data recommend the continued monitoring of this aquatic system.     

Long‐wave radiation and heat flux estimates within a small tributary in Catamaran Brook (New Brunswick,Canada)

River water temperature is an important water quality parameter that also influences most aquatic life. Physical processes influencing water temperature in rivers are highly complex. This is especially true for the estimation of river heat exchange processes that are highly dependent on good estimates of radiation fluxes. Furthermore, very few studies were found within the stream temperature dynamic literature where the different radiation components have been measured and compared at the stream level (at microclimate conditions). Therefore, this study presents results on hydrometeorological conditions for a small tributary within Catamaran Brook (part of the Miramichi River system, New Brunswick, Canada) with the following specific objectives: (1) to compare between stream microclimate and remote meteorological conditions, (2) to compare measured long‐wave radiation data with those calculated from an analytical model, and (3), to calculate the corresponding river heat fluxes. The most salient findings of this study are (1) solar radiation and wind speed are parameters that are highly site specific within the river environment and play an important role in the estimation of river heat fluxes; (2) the incoming, outgoing, and net long‐wave radiation within the stream environment (under the forest canopy) can be effectively calculated using empirical formula; (3) at the study site more than 80% of the incoming long‐wave radiation was coming from the forest; (4) total energy gains were dominated by solar radiation flux (for all the study periods) followed by the net long‐wave radiation (during some periods) whereas energy losses were coming from both the net long‐wave radiation and evaporation. Conductive heat fluxes have a minor contribution from the overall heat budget (<3·5%); (5) the reflected short‐wave radiation at the water surface was calculated on average as 3·2%, which is consistent with literature values. Results of this study contribute towards a better understanding of river heat fluxes and water temperature models as well as for more effective aquatic resources and fisheries management. Copyright © 2011 John Wiley & Sons, Ltd.     

Seasonal fluctuations of tissue mercury contents in the European shore crab Carcinus maenas from low and high contamination areas (Ria de Aveiro, Portugal)

The main objective was to study the seasonal variation of mercury concentrations in different tissues (muscle, hepatopancreas and gills) of Carcinus maenas from low and high Hg contaminated areas, a valuable resource in temperate estuaries and a possible pathway for human uptake. Individuals of two size classes (around 35 and 55 mm cephalothorax wide) were captured monthly between March 1999 and May 2000 in two areas of Ria de Aveiro: in the main navigation channel that connects the lagoon to the sea, and in the inner lagoon area heavily contaminated by mercury (maximum Hg in sediments of 5.4 microg g(-1)). Pronounced decreases in salinity and temperature and reduced food availability in winter seemed to be the responsible for the decline of the crab condition index (0.75-0.45) in larger individuals. Muscle and hepatopancreas exhibited higher mercury concentrations than gills, with concentrations in the contaminated site ranging from 0.03 to 0.63 microg g(-1) and 0.02 to 0.34 microg g(-1), respectively. Linear regressions between muscle and hepatopancreas (r=0.94, p<0.001) and muscle and gills (r=0.97, p<0.001) suggested a rapid redistribution of mercury inside the organism. During winter, a rapid elimination of mercury was found in the three analysed tissues followed by uptake. Larger crabs presented elimination rates from 18 to 34 ng g(-1) per week, while the smaller crabs showed lower elimination rates (10-24 ng g(-1) per week). The uptake was similar in both size classes (11-15 ng g(-1) and 8.1-15 ng g(-1) per week, respectively for large and small crabs). Our results suggest that C. maenas harvested in the contaminated areas must be considered with caution, since Hg concentrations were found to exceed the threshold concentration allowed for human consumption (0.5 microg g(-1)).     

Modelling long‐term (300 ka) upland catchment response to multiple lava damming events

Landscapes respond in complex ways to external drivers such as base level change due to damming events. In this study, landscape evolution modelling was used to understand and analyse long‐term catchment response to lava damming events. PalaeoDEM reconstruction of a small Turkish catchment (45 km²) that endured multiple lava damming events in the past 300 ka, was used to derive long‐term net erosion rates. These erosion rates were used for parameter calibration and led to a best fit parameter set. This optimal parameter set was used to compare net erosion landscape time series of four scenarios: (i) no uplift and no damming events; (ii) no uplift and three damming events; (iii) uplift and no damming events; and (iv) uplift and three damming events. Spatial evolution of net erosion and sediment storage of scenario (iii) and (iv) were compared. Simulation results demonstrate net erosion differences after 250 000 years between scenarios with and without dams. Initially, trunk gullies show less net erosion in the scenario with damming events compared with the scenario without damming events. This effect of dampened erosion migrates upstream to smaller gullies and local slopes. Finally, an intrinsic incision pulse in the dam scenario results in a higher net erosion of trunk gullies while decoupled local slopes are still responding to the pre‐incision landscape conditions. Sediment storage differences also occur on a 100 ka scale. These differences behaved in a complex manner owing to different timings of the migration of erosion and sediment waves along the gullies for each scenario. Although the specific spatial and temporal sequence of erosion and deposition events is sensitive to local parameters, this model study shows the manner in which past short‐lived events like lava dams have long‐lasting effects on catchment evolution. Copyright © 2014 John Wiley & Sons, Ltd.     

Long‐term morphological impacts of the opening of a new inlet on a multiple inlet system

In this study, the artificial opening of a new tidal inlet in an existing multiple inlet system is shown to significantly modify the adjacent nearshore and backbarrier morphology, as well as both updrift and downdrift shorelines. The study focuses on the dominant Faro‐Olhão and Armona inlets in the Ria Formosa barrier island system of southern Portugal. The equilibrium state and future evolution of the system are inferred using a range of morphological and hydrodynamic indicators, including the evolution of the inlet cross‐section, changes in tidal prism, and changes in the dimensions (length and area) of barrier islands. The results reveal how the morphology of an interconnected two‐inlet bay system and the adjacent coastlines has evolved following the artificial opening and stabilization of Faro‐Olhão inlet since 1929. A clear relationship between barrier island size, inlet cross‐section/width, and tidal prism is demonstrated. Decadal time‐scale changes in the tidal prism of the two interconnected inlets are shown to be the main mechanism responsible for morphological change, and have resulted in the remobilization of ebb‐tidal delta sediments deposited during previous hydraulic configurations. These changes, in turn, have contributed to a narrowing of Armona inlet and an increase in the size of Culatra Island. The work highlights the importance of ebb‐tidal deltas both as sand reservoirs and as conduits through which sand exchange between estuaries or lagoons and the open coast is regulated. It also shows the pivotal role of ebb‐tidal deltas in trapping longshore‐transported sediment and releasing it again during periods of increased wave activity. The findings have implications regarding the accurate assessment of the stability of multiple inlet systems. Copyright © 2011 John Wiley & Sons, Ltd.     

Detecting the flow relationships between deep and shallow aquifers in an exploited groundwater system,using long‐term monitoring data and quantitative hydrogeology: the Acque Albule basin case (Rome,Italy)

Five years of hydrogeological monitoring and field activities performed in the complex hydrogeological system of the Acque Albule basin (AAB) were conducted to define the hydrogeological setting, the relationship between deep and shallow aquifers and a conceptual groundwater flow model of this exploited area using conventional quantitative techniques. The basin, which is located close to Rome (Italy) on the west side of the Apennine chain and just north of the Colli Albani volcano, subsided after development of a north–south fault system (about 115 000 y bp). The AAB experiences intense hydrothermal activity, which has produced a large travertine deposit (80‐m thick). The travertine deposit constitutes a fractured aquifer that is the final destination of more than 5 m³ s^‐1 of water and is strongly dewatered by quarry activities. The complex hydrogeology of this basin was investigated, revealing two main hydraulically connected aquifers, one thermalised and partly confined into the limestone bedrock and one unconfined in the travertine. The two aquifers are separated by a non‐continuous clayey aquiclude. The hydrogeological survey and geological characterisation contributed to the development of the groundwater flow conceptual model. Analysis and comparison of the monitored levels highlighted the pattern of flow between the deep and shallow parts of the flow system. Copyright © 2012 John Wiley & Sons, Ltd.     

Integrative assessment of organotin contamination in a southern European estuarine system (Ria de Aveiro, NW Portugal): tracking temporal trends in order to evaluate the effectiveness of the EU ban

Organotin (OT) compounds have been used as biocide agents in antifouling paints since the mid 1960s and are now ubiquitous in the marine environment. Due to their high toxicity to non-target species they were banned from antifouling paints in the European Union in 2003 (2002/62/EC directive). The aim of the present work is to assess any obvious decline of the OT environmental levels at Ria de Aveiro (NW Portugal) after the ban. The organotins - monobutyltin (MBT), dibutyltin (DBT), tributyltin (TBT), monophenyltin (MPT), diphenyltin (DPT), triphenyltin (TPT), monoctyltin (MOT), dioctyltin (DOT) and trioctyltin (TOT) - were quantified in the gastropod Nassarius reticulatus, in the mussel Mytilus galloprovincialis and in sediments. Imposex (imposition of male characters on females of gonochorist gastropods) in N. reticulatus was additionally used as a biomarker of TBT pollution. Time comparisons show a slight decrease of imposex between 2003 and 2005 probably as a consequence of the EU ban, though in some cases this trend seems to have started earlier since 2000. The fraction of TBT relatively to its metabolites has been decreasing over the last years but still remains high suggesting that there are still recent inputs of this compound into the study area.     

The impact of exceptional events on erosion,bedload transport and channel stability in a step‐pool channel

Sediment transport in the Erlenbach, a small stream with step‐pool morphology in the canton of Schwyz, Switzerland, has been monitored for more than 20 years. During this time three exceptional events (events with high sediment yield and long return times that have a large effect on channel morphology) have impacted the stream and partly or completely rearranged the existing step‐pool morphology. In the aftermath of the events, sediment transport rates at a given discharge and total sediment yield remained elevated for about a year or longer. For the last event, dated on the 20 June 2007, observations of boulder mobility and step destruction were used to interpret channel stability. Boulders with median diameters of up to 135 cm and estimated weights of more than 2·5 tons have moved during the 2007 event. Using hydraulic observations and shear stress calculations boulders up to 65 cm in diameter were predicted to have been fully mobile in peak conditions, even if form resistance and increased critical stresses needed for the initiation of motion in steep streams were taken into account. For two of the events, estimated peak shear stresses at the bed exceeded 1000 Pa, calculated both from observations of the flow hydraulics and from boulder mobility. This suggests that highly energetic flows occur relatively frequently in small, steep streams and that large boulders can be transported by fluvial processes in such streams. The observations have potential significance for hazard risk mitigation, stream engineering and restoration. Copyright © 2009 John Wiley & Sons, Ltd.     

Estimating stream solute loads from fixed frequency sampling regimes: the importance of considering multiple solutes and seasonal fluxes in the design of long‐term stream monitoring networks

Reduced sampling frequency is known to increase the error associated with estimates of stream solute load. However, the extent to which the magnitude of error differs among commonly measured solutes and across seasons is unclear. In this study, a high sampling frequency data set from two forested streams (one upland‐draining and one wetland‐draining stream) in south‐central Ontario was systematically sub‐sampled to simulate weekly, fortnightly and monthly fixed frequency sampling regimes for 12 stream solutes. We found that solutes which had a higher degree of temporal variation in concentration (i.e. higher %RSD) had poorer precision (C_v) in estimates of annual load relative to solutes with a lower %RSD. In addition, the magnitude and direction of bias varied considerably among solutes and were related to differences in spring concentration‐discharge relationships (m_{[spring Q vs C]}) among the 12 solutes. Solutes which decreased in concentration with increases in spring flow (i.e. m_{[spring Q vs C]} <0) exhibited positive bias in annual load while solutes which increased in concentration with increases in spring flow (i.e. m_{[spring Q vs C]} >0) were negatively biased. In terms of differences between seasonal and annual load errors, precision was generally lower for estimates of seasonal load relative to annual load while bias varied in both magnitude and direction among seasons. When the root mean square error (RMSE) of load estimates was compared to a threshold of acceptable error (<15%), the proportion of solutes attaining acceptable levels of uncertainty ranged from 11/12 for annual load estimates at a weekly sampling frequency to only 4/12 at a monthly frequency when both annual and seasonal loads were considered. Our results demonstrate that commonly measured solutes do not behave uniformly in response to changes in sampling frequency and that estimates of seasonal loads are often less accurate than estimates of annual load. These findings provide important insights into the design of stream monitoring programs and the evaluation of existing long‐term data sets. Copyright © 2015 John Wiley & Sons, Ltd.     

Long‐term effects of management practices on water‐driven soil erosion in an intense agricultural sub‐watershed: monitoring and modelling

The overarching objective of this research was to provide an improved understanding of the role of land use and associated management practices on long‐term water‐driven soil erosion in small agricultural watersheds by coupling the established, physically based, distributed parameter Water Erosion Prediction Project (WEPP) model with long‐term hydrologic, land use and soil data. A key step towards achieving this objective was the development of a detailed methodology for model calibration using physical ranges of key governing parameters such as effective hydraulic conductivity, critical hydraulic shear stress and rill/inter‐rill erodibilities. The physical ranges for these governing parameters were obtained based on in situ observations within the South Amana Sub‐Watershed (SASW) (～26 km²) of the Clear Creek, IA watershed where detailed documentation of the different land uses was available for a period of nearly 100 years. A quasi validation of the calibrated model was conducted through long‐term field estimates of water and sediment discharge at the outlet of SASW and also by comparing the results with data reported in the literature for other Iowa watersheds exhibiting similar biogeochemical properties. Once WEPP was verified, ‘thought experiments’ were conducted to test our hypothesis that land use and associated management practices may be the major control of long‐term erosion in small agricultural watersheds such as SASW. Those experiments were performed using the dominant 2‐year crop rotations in the SASW, namely, fall till corn–no till bean (FTC‐NTB), no till bean–spring till corn (NTB‐STC) and no till corn–fall till bean (NTC‐FTB), which comprised approximately 90% of the total acreage in SASW. Results of this study showed that for all crop rotations, a strong correspondence existed between soil erosion rates and high‐magnitude precipitation events during the period of mid‐April and late July, as expected. The magnitude of this correspondence, however, was strongly affected by the crop rotation characteristics, such as canopy/residue cover provided by the crop, and the type and associated timing of tillage. Tillage type (i.e. primary and secondary tillages) affected the roughness of the soil surface and resulted in increases of the rill/inter‐rill erodibilities up to 35% and 300%, respectively. Particularly, the NTC‐FTB crop rotation, being the most intense land use in terms of tillage operations, caused the highest average annual erosion rate within the SASW, yielding quadrupled erosion rates comparatively to NTB‐STC. The impacts of tillage operation were further exacerbated by the timing of the operations in relation to precipitation events. Timing of operations affected the ‘life‐time’ of residue cover and as a result, the degree of protection that residue cover offers against the water action on the soil surface. In the case of NTC‐FTB crop rotation, dense corn residue stayed on the ground for only 40 days, whereas for the other two rotations, corn residue provided a protective layer for nearly 7 months, lessening thus the degree of soil erosion. The cumulative effects of tillage type and timing in conjunction with canopy/residue cover led to the conclusion that land management practices can significantly amplify or deamplify the impact of precipitation on long‐term soil erosion in small agricultural watersheds. Copyright © 2009 John Wiley & Sons, Ltd.     

Trends in timing of low stream flows in Canada: impact of autocorrelation and long‐term persistence

The annual timing of river flows might indicate changes that are climate related. In this study, trends in timing of low flows for the Reference Hydrometric Basin Network were investigated under three different hypotheses namely: independence, short‐term persistence (STP) and long‐term persistence (LTP). Both summer and winter time series were characterized with scaling behaviour providing strong evidence of LTP. The Mann–Kendall trend test was modified to account for STP and LTP, and used to detect trends in timing of low flows. It was found that considering STP and LTP resulted in a significant decrease in the number of detected trends. Numerical analysis showed that the timing of summer 7‐day low flows exhibited significant trends in 16, 9 and 7% of stations under independence, STP and LTP assumptions, respectively. Timing of summer low flow shifted toward later dates in western Canada, whereas the majority of stations in the east half of the country (except Atlantic Provinces) experienced a shift toward earlier dates. Timing of winter low flow experienced significant trends in 20, 12, and 6% of stations under independence, STP and LTP assumptions, respectively. Shift in timing of winter low flow toward earlier dates was dominant all over the country where it shifted toward earlier dates in up to 3/4 of time series with significant trends. There are local patterns of upward significant/insignificant trends in southeast, southwest and northern Canada. This study shows that timing of low flows in Canada is time dependent; however, addressing the full complexity of memory properties (i.e. short term vs long term) of a natural process is beyond the scope of this study. Copyright © 2009 John Wiley & Sons, Ltd.     

Long‐term Holocene groundwater fluctuations in a chalk catchment: evidence from Rock‐Eval pyrolysis of riparian peats

The depositional history of peat‐dominated wetlands can be used to understand palaeoclimate and palaeohydrology and also constrain the impacts of future climate change. However, in chalkland valleys, seasonal water table fluctuations and a high alkalinity have diminished key environment indicators such as pollen, and there is a need for alternative investigative techniques. The method of Rock‐Eval pyrolysis can track changes in organic matter source and degradation, potentially relating to historic changes in vegetation cover. This is the first Rock‐Eval on cores from a groundwater‐dependent riparian chalk valley wetland combined with radiocarbon dating. The dating showed that the cores represented approximately 4000 years of depositional history. Changes in hydrocarbon chemistry including normal alkane composition of the peat indicated shifts of around 500 to 1000 years between terrestrial and more aquatic species, relating to periods of climate wetness. These climatic shifts are broadly consistent with other evidence from ombrotrophic peatland and lacustrine sediments across northwest Europe. However, the connection between climate wetness and groundwater dependent chalkland wetlands is complicated by external anthropogenically driven factors relating to land use and vegetation cover changes in the catchment. Nonetheless, this study suggests that Rock‐Eval pyrolysis is a useful and cost‐effective tool that can provide evidence for long‐term Holocene groundwater fluctuations. Copyright © 2016 British Geological Survey.