OR/18/026 Water quality
R S Ward1, P L Smedley1, G Allen2, B J Baptie1, M R Cave1, Z Daraktchieva3, R Fisher5, D Hawthorn1, D G Jones1, A Lewis4, D Lowry5, R Luckett1, B P Marchant1, R M Purvis4 and S Wilde4. 2018. Environmental baseline monitoring: phase III final report (2017–2018). British Geological Survey Internal Report, OR/18/026.
With contributions from: |
Introduction
Monitoring of both surface water and groundwater in the Vale of Pickering has continued over the grant period from April 2017 to March 2018. The water monitoring network (Smedley et al., 2015[1]) was initially set up for sampling in September 2015 and augmented with sampling sites provided by new BGS purpose-drilled boreholes from March 2016. We have also been sampling Third Energy’s boreholes at KM8 since April 2016 (Figure 1). New multi-level samplers for groundwater sampling, installed some 2 km east of KM8, have also been monitored frequently since May 2017. Two multi-level samplers have been installed: one with a depth range of 4–20 m (Quaternary extends to depth 7 m) and another with a range of 8–75 m (Kimmeridge Clay extends from 7 m to ca.200 m). The sampler ranges overlap slightly but are designed to sample mainly from the Quaternary and mainly from the Kimmeridge Clay respectively.
The monitoring network aims to establish baseline chemical characteristics of groundwater and surface water, the former from a Superficial aquifer (Quaternary ± topmost Kimmeridge Clay; depicted as Upper Jurassic clays in Figure 1) and the Corallian Limestone aquifer; the latter from low-order streams dispersed across the Vale of Pickering. The monitoring infrastructure and sampling protocols have been described in detail by Ward et al. (2017)[2]. At the time of reporting (March 2018), monitoring of the initial network has been carried out 21 times (with additional samples at other times), the BGS boreholes 21 times, Third Energy Boreholes 20 times and the BGS multi-level samplers 7 times. Frequency of sampling of the network sites and BGS boreholes was initially monthly for 12 months, according to the requirements (for methane) of the 2015 Infrastructure Act. Thereafter, sampling reduced to quarterly, until the point where hydraulic fracturing looked likely to start (October 2017). Monthly sampling was conducted in each from October 2018–January 2018, with a final sampling for the reporting period in March 2018. Laboratory data are here reported mostly up to January 2018.
In addition, summary data from the five sondes, installed to collect real-time quality and physical-properties data in the BGS monitoring boreholes, are reported for the period up to March 2018.

Monitoring network data
Data for the sites in the monitoring network from the Superficial aquifer are shown in Figure 2–Figure 4. Groundwater pH varies spatially from 7.0–8.8 with broad consistency in values for individual sites over time. Redox potential is dominantly between -200–+200 mV, indicating the mildly to strongly reducing conditions across the Superficial aquifer. Specific electrical conductance (SEC) shows a large range spatially from 500–3000 µS/cm but with a notable consistency at individual sites over time. Groundwater temperature commonly reflects a dampened response to ambient air temperature but with greatest response tending to be in the shallowest boreholes. In the boreholes with the most extreme temperature fluctuations, some artefact from the sampling procedure is possible, including changes in water level and flow induced by pumping.

Data for major and minor ions (Figure 3, Figure 4) also indicate a range of compositions in the groundwater from the Superficial aquifer but most have a notable consistency over time. The data confirm the Na-HCO3 dominance of the water types in the aquifer, with reducing conditions supported by high concentrations of dissolved NH4, Fe and Mn in many. By far the largest variability is seen in groundwater sampled from Site 15. This borehole is on a site with numerous nearby boreholes and was disused before being incorporated into the monitoring network. Changes over time could reflect changes induced by pumping of the borehole itself, effects of pumping neighbouring boreholes, or a combination of the two. The borehole was pumped more regularly during February–March 2016, November–December 2016 and January 2017 (pump rates unknown). The changes in water chemistry at Site 15 over time imply some stratification of water with depth, with draw-in from different horizons depending on local pumping rates and regime. These observations highlight the care required in designing/selecting and operating groundwater monitoring points. It also illustrates the need for multiple monitoring points to ensure that results are representative of groundwater quality in the area and that future changes that might arise from shale gas site operations can be detected.
Concentrations of dissolved CH4 are repeatedly high in groundwater from many of the boreholes (Figure 4), several exceeding 1 mg/L and two showing groundwater with consistently greater than 5 mg/L and maxima to 50 mg/L, albeit with substantial variation between sampling visits. Data for dissolved CO2 show a strong inverse relationship with CH4 in Site 15 (Figure 4), suggesting varying redox status over time, in further response to changing pumping rates and inflow horizons. In combination with observations for the major/minor this site is displaying very unusual behaviour, when compared to other nearby sites. This will be investigated further and the outcomes used to inform future site selection, operation and sampling guidance.

Data for groundwater from the new BGS boreholes in the Superficial aquifer (Figure 5–Figure 7) show similar features to the pre-existing boreholes, with Na-HCO3 dominance, pH varying from 6.5–8.5, Eh -100 to +200 mV, SEC 1000–3500 µS/cm, and temperature varying seasonally. The data show a greater temporal variability but still with broad consistency in chemical composition. This reflects the greater number of monitoring points and data now available.
Groundwater in the BGS boreholes is reducing with respect to Fe, Mn, NH4 and CH4, and CH4 concentrations reach up to 40 mg/L in line with values observed in the Superficial aquifer groundwater monitoring network.




Groundwater in the Corallian Limestone aquifer from the margins of the Vale of Pickering (Figure 1) contrasts with that from the Superficial clay-dominated aquifer, being of Ca-HCO3 composition, with a pH range of 7.2–8.0, Eh from 0–500 mV reflecting a large range of redox conditions from oxic to anoxic, and SEC values of 300–750 µS/cm (lower than in the Superficial aquifer). Temperature is mainly in the 10–11°C range, although temperatures for Sites 7 and 16 are more variable (Figure 8). Redox-sensitive minor-element concentrations reflect the range of redox conditions, with some showing high Fe and Mn concentrations and NH4 usually <0.2 mg/L but reaching up to 0.5 mg/L (Figure 9, Figure 10). These very occasional peaks may reflect incomplete purging of boreholes before sampling. Trends in CH4 concentration over time are spiky, but with much lower concentrations (and range) than observed in the Superficial aquifer (up to 0.7 mg/L in the Corallian; Figure 10). Compositions of many of the solutes in the Corallian groundwater are relatively stable over time compared to observations from the new BGS boreholes in the Superficial aquifer.



In the monitored streams, pH values are typically 7.0–8.5 with SEC mainly between 200 and 800 µS/cm and occasionally up to 1400 µS/cm (Figure 11, Figure 8). The stream site with the highest SEC values (Site 29) also periodically has the highest concentrations of Na and Cl as well as the greatest variability in these ions over time (Figure 12). For other major ions and many of the minor ions, compositions are less variable. Occasional spikes in concentrations of Fe may be from material present as a fine colloidal suspension reflecting higher flow rates.



Monitoring statistics
Statistical summary data for selected important analytes in groundwater from the Superficial aquifer and for streams, each collected between September 2015 and January 2018, are provided as violin plots in Figure 14, Figure 15 and Figure 16. These display median, interquartile range and range in the box plots as well as probability density functions (the outer envelope showing relative likelihood of a concentration occurring). The data are divided into: all groundwater data for the period; data in a defined API (‘area of potential impact’, categorised as within an approximate 2 km radius around KM8), control area outside this zone; data for 2016; data for 2017. Data for streams also include all data for the period as well as 2016 and 2017 data. The divisions investigate baseline chemical distributions both spatially and temporally; the time differences have been included to investigate similarities and differences between the distributions for two separate 12-month intervals of monitoring. The summary distributions for groundwater are given in Table 1 and for streams in Table 2. Once a complete baseline dataset is available, the statistics will be re-run and a range of thresholds identified, above which any detection in groundwater will suggest that a deviation from baseline has occurred with an associated level of confidence.



Analyte | Units | Group | n | mean | sd | median | mad | min | max | range | skew |
Ca | mg/L | All data | 493 | 80.4 | 86.8 | 45.7 | 46.1 | 1.40 | 370 | 369 | 1.57 |
Mg | mg/L | All data | 493 | 9.91 | 7.95 | 7.85 | 5.44 | 0.50 | 40.0 | 39.5 | 1.71 |
Na | mg/L | All data | 493 | 249 | 195 | 190 | 227 | 14.9 | 834 | 819 | 0.69 |
K | mg/L | All data | 493 | 2.88 | 1.30 | 2.56 | 1.04 | 0.96 | 6.72 | 5.76 | 0.88 |
Cl | mg/L | All data | 493 | 104 | 94.2 | 71.8 | 53.9 | 18.7 | 443 | 425 | 1.76 |
SO4 | mg/L | All data | 493 | 63.0 | 91.5 | <0.25 | 857 | ||||
NH4 | mg/L | All data | 473 | 0.95 | 0.71 | 0.83 | 0.73 | 0.010 | 3.12 | 3.11 | 0.87 |
Ba | µg/L | All data | 493 | 198 | 293 | 90.0 | 92.8 | 7.60 | 1380 | 1370 | 2.44 |
Sr | µg/L | All data | 493 | 737 | 759 | 499 | 340 | 86.0 | 4230 | 4145 | 2.43 |
Li | µg/L | All data | 493 | 35.7 | 27.3 | 30.0 | 26.7 | 3.00 | 120 | 117 | 0.73 |
CH4 | mg/L | All data | 474 | 0.20 | 0.30 | <0.001 | 48.8 | ||||
Ca | mg/L | API | 375 | 84.8 | 97.4 | 34.0 | 36.2 | 1.40 | 370 | 369 | 1.34 |
Mg | mg/L | API | 375 | 10.39 | 8.90 | 7.60 | 6.12 | 0.50 | 40.0 | 39.5 | 1.46 |
Na | mg/L | API | 375 | 296 | 199 | 248 | 201 | 19.6 | 834 | 814 | 0.34 |
K | mg/L | API | 375 | 3.01 | 1.41 | 2.87 | 1.47 | 0.96 | 6.72 | 5.76 | 0.65 |
Cl | mg/L | API | 375 | 124 | 100 | 101 | 74.2 | 19.2 | 444 | 424 | 1.42 |
SO4 | mg/L | API | 375 | 67.0 | 99.0 | <0.25 | 857 | ||||
NH4 | mg/L | API | 361 | 1.04 | 0.75 | 0.88 | 0.78 | 0.010 | 3.12 | 3.11 | 0.75 |
Ba | µg/L | API | 375 | 127 | 146 | 70.9 | 69.1 | 7.60 | 753 | 745 | 2.29 |
Sr | µg/L | API | 375 | 808 | 845 | 514 | 418 | 86.0 | 4230 | 4140 | 2.09 |
Li | µg/L | API | 375 | 41.1 | 27.8 | 36.0 | 31.1 | 5.00 | 120 | 115 | 0.55 |
CH4 | mg/L | API | 361 | 0.29 | 0.43 | <0.001 | 48.8 | ||||
Ca | mg/L | Control | 118 | 66.5 | 33.3 | 57.3 | 47.1 | 15.9 | 129.5 | 114 | 0.20 |
Mg | mg/L | Control | 118 | 8.36 | 3.03 | 9.00 | 3.12 | 2.95 | 13.6 | 10.7 | -0.45 |
Na | mg/L | Control | 118 | 99.3 | 67.7 | 94.3 | 92.7 | 14.9 | 283 | 268 | 0.49 |
K | mg/L | Control | 118 | 2.47 | 0.72 | 2.14 | 0.39 | 1.69 | 5.04 | 3.35 | 1.29 |
Cl | mg/L | Control | 118 | 42.3 | 12.8 | 41.4 | 9.41 | 18.7 | 97.7 | 79.0 | 1.39 |
SO4 | mg/L | Control | 118 | 50.3 | 35.1 | <0.25 | 108 | ||||
NH4 | mg/L | Control | 112 | 0.69 | 0.47 | 0.46 | 0.45 | 0.15 | 1.76 | 1.61 | 0.43 |
Ba | µg/L | Control | 118 | 427 | 472 | 124 | 97.0 | 36.2 | 1380 | 1340 | 0.81 |
Sr | µg/L | Control | 118 | 512 | 274 | 384 | 119 | 252 | 1160 | 905 | 1.21 |
Li | µg/L | Control | 118 | 18.4 | 16.4 | 10.0 | 8.90 | 3.00 | 50.4 | 47.4 | 0.74 |
CH4 | mg/L | Control | 113 | 0.10 | 0.15 | <0.01 | 24.0 | ||||
Ca | mg/L | 2016 | 238 | 80.6 | 79.6 | 51.7 | 52.5 | 1.50 | 299 | 297 | 1.31 |
Mg | mg/L | 2016 | 238 | 10.1 | 7.71 | 8.39 | 5.52 | 0.53 | 40.0 | 39.4 | 1.67 |
Na | mg/L | 2016 | 238 | 252 | 197 | 199 | 234 | 15.2 | 834 | 818 | 0.68 |
K | mg/L | 2016 | 238 | 2.97 | 1.35 | 2.69 | 1.22 | 0.96 | 6.44 | 5.48 | 0.71 |
Cl | mg/L | 2016 | 238 | 110 | 101 | 72.3 | 53.7 | 19.2 | 443 | 424 | 1.69 |
SO4 | mg/L | 2016 | 238 | 139 | 184 | 61.3 | 89.9 | <0.25 | 855 | ||
NH4 | mg/L | 2016 | 218 | 1.00 | 0.74 | 0.84 | 0.78 | 0.01 | 3.12 | 3.11 | 0.78 |
Ba | µg/L | 2016 | 238 | 212 | 298 | 94.2 | 103 | 8.00 | 1270 | 1270 | 2.14 |
Sr | µg/L | 2016 | 238 | 780 | 777 | 505 | 381 | 87.6 | 4230 | 4140 | 2.25 |
Li | µg/L | 2016 | 238 | 39.3 | 29.5 | 33.5 | 30.4 | 4.00 | 120 | 116 | 0.63 |
CH4 | mg/L | 2016 | 236 | 0.20 | 0.30 | <0.01 | 44.2 | ||||
Ca | mg/L | 2017 | 181 | 88.1 | 101 | 40.4 | 40.3 | 1.40 | 370 | 369 | 1.44 |
Mg | mg/L | 2017 | 181 | 9.84 | 8.18 | 7.70 | 5.43 | 0.50 | 40.0 | 39.5 | 1.69 |
Na | mg/L | 2017 | 181 | 251 | 196 | 213 | 246 | 14.9 | 676 | 661 | 0.60 |
K | mg/L | 2017 | 181 | 2.79 | 1.21 | 2.55 | 0.99 | 1.09 | 6.72 | 5.63 | 1.07 |
Cl | mg/L | 2017 | 181 | 100 | 80.8 | 73.3 | 54.9 | 18.8 | 365 | 347 | 1.62 |
SO4 | mg/L | 2017 | 181 | 66.8 | 97.2 | <0.25 | 783 | ||||
NH4 | mg/L | 2017 | 181 | 0.91 | 0.68 | 0.83 | 0.72 | 0.02 | 2.90 | 2.88 | 0.83 |
Ba | µg/L | 2017 | 181 | 189 | 295 | 85.1 | 82.3 | 7.60 | 1380 | 1370 | 2.69 |
Sr | µg/L | 2017 | 181 | 690 | 705 | 494 | 309 | 86.0 | 3860 | 3770 | 2.64 |
Li | µg/L | 2017 | 181 | 33.5 | 23.7 | 32.0 | 28.2 | 3.00 | 87.0 | 84.0 | 0.55 |
CH4 | mg/L | 2017 | 181 | 0.19 | 0.28 | <0.01 | 44.9 |
Notes: All data: samples collected September 2015 to January 2018; API: area of potential impact (ca.2 km radius of KM8); control: Superficial aquifer beyond 2 km radius; data for 2016 and 2017 only are also given. Min, max: minimum and maximum observed values; sd: standard deviation; mad: median absolute deviation. Values for mean, standard deviation, range and skew are not computed for data with non-detects to avoid assumptions about their distributions
Analyte | Units | n | mean | sd | median | mad | min | max | range | Skew |
Ca | mg/L | 205 | 95.4 | 31.0 | 98.3 | 21.7 | 13.8 | 219 | 205 | 0.07 |
Mg | mg/L | 205 | 5.89 | 2.22 | 5.63 | 2.56 | 1.99 | 16.8 | 14.8 | 0.66 |
Na | mg/L | 205 | 25.8 | 37.7 | 15.9 | 5.34 | 9.10 | 304 | 295 | 5.20 |
K | mg/L | 205 | 3.54 | 11.8 | 1.86 | 0.86 | 0.89 | 169 | 168 | 13.3 |
Cl | mg/L | 205 | 42.1 | 42.6 | 30.8 | 8.19 | 11.8 | 340 | 329 | 4.64 |
SO4 | mg/L | 205 | 45.3 | 23.7 | 43.1 | 15.8 | 3.23 | 121 | 118 | 1.45 |
NO3 | mg/L | 205 | 26.9 | 17.8 | <0.15 | 63.3 | ||||
NH4 | mg/L | 195 | 0.04 | 0.04 | <0.004 | 3.73 | ||||
Ba | µg/L | 205 | 75.8 | 32.6 | 70.4 | 17.9 | 41.5 | 311 | 269 | 3.42 |
Sr | µg/L | 205 | 189 | 79.8 | 161 | 67.7 | 48.6 | 467 | 418 | 0.71 |
Li | µg/L | 205 | 3.00 | 1.48 | <1 | 36.0 |
For groundwater from the Superficial aquifer, data from the API and control areas (all dates) have been compared using the non-parametric Wilcoxon signed rank test of significance (Table 3). Statistically significant differences are observed for Na, Cl, SO4, NH4 and Li between the API and control data sets. Significant differences are not observed for Ca, Mg, K, Sr and CH4. For several indicators of salinity therefore, baseline concentrations in the groundwater closer to KM8 (i.e. in the central part of the Vale of Pickering Superficial aquifer) are significantly higher than in the control area (median and maximum values respectively for Na in the API being 248 mg/L/834 mg/L, for control area being 94.3 mg/L/283 mg/L; median and maximum values respectively for Cl in the API being 101 mg/L/444 mg/L, for control area being 41.4 mg/L/97.7 mg/L). It is important to have established this distinction under baseline conditions before considering regional variations in chemistry that may occur under post- baseline (operational) conditions.
Test | Ca | Mg | Na | K | Cl | SO4 | NH4 | Ba | Sr | Li | CH4 |
Two-sided | NS | NS | S | NS | S | S | S | S | NS | S | NS |
One-sided | NS | NS | S | NS | S | S | S | NS | NS | S | NS |
Table 4 and Table 5 show the comparison of results for 2016 and 2017 data from groundwater in the Superficial aquifer and from streams respectively. For the analytes considered, no statistically significant differences (at 99% confidence level) were observed. From the sample sites and the analytes investigated, there is no indication that significant differences occurred in the water chemistry from the two 12-month intervals of sampling. This suggests that a 12-month period of monthly monitoring ahead of any site operations may be sufficient to characterise the baseline.
Test | Ca | Mg | Na | K | Cl | SO4 | NH4 | Ba | Sr | Li | CH4 |
Two-sided | NS | NS | NS | NS | NS | NS | NS | NS | NS | NS | NS |
One-sided | NS | NS | NS | NS | NS | NS | NS | NS | NS | NS | NS |
Test | Ca | Mg | Na | K | Cl | SO4 | NH4 | Ba | Sr | Li | CH4 |
Two-sided | NS | NS | NS | NS | NS | NS | NS | NS | NS | NS | NS |
One-sided | NS | NS | NS | NS | NS | NS | NS | NS | NS | NS | NS |
Real-time groundwater data
Results for hourly monitoring of groundwater quantity and quality (pH, SEC, temperature and groundwater level) from five sondes over the last year of monitoring (March 2017 to March 2018) are displayed in Figure 17. Short-term spikes in the time series indicate intervals when sondes were removed for maintenance or calibration and have been retained in the record to indicate frequency of intervention. We experienced significant problems with EBM8 over the monitored interval despite attempts to rectify. For this sonde, some values for pH, SEC and temperature from September 2017 to March 2018 are spurious. The sonde has since been repaired. The pH electrode in EBM9 also gave spurious data in early September 2017, since repaired. Other fluctuations over the year are believed to be real variations; pH shows variations typically of 0.5 pH units (up to 1 pH unit), SEC of up to 1000 µS/cm and water level of <1 m (Figure 17). Temperature has been the most stable parameter, varying between 10–11°C for the boreholes observed. Hourly monitoring is continuing.


Figure 18 shows a comparison of pH and SEC values determined during spot sampling in the monitoring rounds from the BGS boreholes (Superficial aquifer) with values determined downhole from the installed sondes at the approximate time of sampling. Values are broadly comparable for each parameter and provide confidence in the data produced. Data for Site BGS 45 (EBM9) are an exception, with both pH and SEC showing systematically higher values in the sonde data than the monitoring samples, albeit only by 0.5 units for pH. Spot sampling involves recording of field readings while pumping until stable values are obtained. In the case of Site BGS 45, on most sampling occasions, SEC values started at around 3500 µS/cm but dropped during the course of monitoring (periods of up to 1 hour). This suggests that low-flow pumping of the borehole over time induces flow with a differing composition from different horizons in the aquifer sediments. This may point to salinity stratification with depth in this borehole, and will be investigated further using depth profiling. Measurements of SEC made during spot sampling are reported after stable readings have been obtained and therefore reflect the lower SEC values.
Profiles from multi-level samplers
Profiles of groundwater from the two multi-level samplers indicate some notable consistency over the successive sampling rounds between June 2017 and February 2018 (Figure 19). Where depths of the two sampler systems overlap, differences at a given depth are likely in part due to spatial heterogeneity. However, systematic differences in redox status are also apparent, with more oxidising conditions indicated (higher SO4, CO2, lower CH4) in the shallower CMT system and more strongly reducing conditions in the Waterloo MLS system.
Low SO4 concentrations in the MLS groundwater, accompanied by a notable sulphide smell, are consistent with dissimilatory SO4 reduction. The strongly reducing conditions in the MLS system are also indicated by high concentrations of CH4 and C2H6 in the groundwater (up to 80 mg/L and 800 µg/L respectively). The simplest explanation is that these hydrocarbons are generated in situ within the Kimmeridge Clay, rather than deriving from extraneous sources.

Summary
A high-quality groundwater baseline dataset is being compiled. It is providing important insights into the importance of establishing robust information on the conditions before shale gas operations start and also some of the challenges in achieving it. Monthly/quarterly monitoring has shown the greatest consistency in chemical compositions occurring in groundwater from sites in the original monitoring network (both Superficial and Corallian aquifers). Streams also showed a broad consistency over time. Greatest variability was seen in groundwater from the BGS new borehole sites. This suggests that these boreholes, which have only been pumped under low-flow conditions due to the small water yields, may not yet have been flushed and reached a steady state. This is an important finding for future monitoring borehole design and operation.
Wilcoxon signed rank tests indicate that for groundwater from the Superficial aquifer, salinity (Na, Cl, SO4) as well as NH4 and Li differ significantly between the aquifer around KM8 (API) and that in the control area. This could be due to differences in e.g. Quaternary lithology (clay representation) and thickness, groundwater flow regimes and redox conditions between the central part of the vale and the more peripheral areas. Comparing statistical summary data for 2016 versus 2017 for groundwater from the Superficial aquifer and for streams showed that sample populations were not significantly different in either case. These statistical distributions will continue to be investigated as more baseline monitoring data become available.
Real-time sensor data for five boreholes monitored hourly over the year of reporting show greatest fluctuations for pH and SEC. Although some of the variability has been due to analytical artefacts, underlying fluctuations are likely due to real variability in flow and chemistry within the boreholes. Fluctuations have been greater in EBM8 and EBM9 than the other boreholes. EBM9 in particular is found to vary significantly in its discharge chemistry over a period of pumping. The real-time data are providing a further measure of baseline groundwater chemistry in the Vale of Pickering at times not sampled by the monitoring campaigns and will continue to provide evidence of any material change in water chemistry should the state of monitoring change from baseline to operational conditions.
Groundwater data from multi-level samplers installed into both Quaternary and Kimmeridge Clay formations provide much greater spatial resolution on groundwater chemical compositions. Groundwater from the Kimmeridge Clay is strongly reducing (sulphate-reducing and methanogenic). In-situ formation of methane (and ethane) within this shallow formation can explain the high concentrations of these gases observed in the shallow groundwater elsewhere in the Vale of Pickering.
References
- ↑ Smedley, P L, Ward, R S, Allen, G, Baptie, B, Daraktchieva, Z, Jones, D G, Jordan, C J, Purvis, R M, and Cigna, F. 2015. Site selection strategy for environmental monitoring in connection with shale-gas exploration, Vale of Pickering, Yorkshire and Fylde, Lancashire. BGS Report, OR/15/067.
- ↑ Ward, R S, Smedley, P L, Allen, G, Baptie, B J, Daraktchieva, Z, Horleston, A, Jones, D G, Jordan, C J, Lewis, A, Lowry, D, Purvis, R M, Rivett, M O. 2017. Environmental Baseline Monitoring Project. Phase II, final report. BGS Report OR/17/049.