Scholarly article on topic 'Assessing Connectivity Between an Overlying Aquifer and a Coal Seam Gas Resource Using Methane Isotopes, Dissolved Organic Carbon and Tritium'

Assessing Connectivity Between an Overlying Aquifer and a Coal Seam Gas Resource Using Methane Isotopes, Dissolved Organic Carbon and Tritium Academic research paper on "Earth and related environmental sciences"

Share paper
Academic journal
Sci. Rep.

Academic research paper on topic "Assessing Connectivity Between an Overlying Aquifer and a Coal Seam Gas Resource Using Methane Isotopes, Dissolved Organic Carbon and Tritium"


Received: 22 July 2015 Accepted: 07 October 2015 Published: 04 November 2015

Assessing Connectivity Between an Overlying Aquifer and a Coal Seam Gas Resource Using Methane Isotopes, Dissolved Organic Carbon and Tritium

Charlotte P. Iverach1, Dioni I. Cendon1'2, Stuart I. Hankin2, David Lowry3, Rebecca E. Fisher3, James L. France4, Euan G. Nisbet3, Andy Baker1 & Bryce F. J. Kelly1

Coal seam gas (CSG) production can have an impact on groundwater quality and quantity in adjacent or overlying aquifers. To assess this impact we need to determine the background groundwater chemistry and to map geological pathways of hydraulic connectivity between aquifers. In southeast Queensland (Qld), Australia, a globally important CSG exploration and production province, we mapped hydraulic connectivity between the Walloon Coal Measures (WCM, the target formation for gas production) and the overlying Condamine River Alluvial Aquifer (CRAA), using groundwater methane (CH4) concentration and isotopic composition (513C-CH4), groundwater tritium (3H) and dissolved organic carbon (DOC) concentration. A continuous mobile CH4 survey adjacent to CSG developments was used to determine the source signature of CH4 derived from the WCM. Trends in groundwater 513C-CH4 versus CH4 concentration, in association with DOC concentration and 3H analysis, identify locations where CH4 in the groundwater of the CRAA most likely originates from the WCM. The methodology is widely applicable in unconventional gas development regions worldwide for providing an early indicator of geological pathways of hydraulic connectivity.

Unconventional gas production, which is increasingly important to the global energy industry, is the focus of major environmental questions. Debates about the impact of gas production can hinge on claims about CH4 leaks and emission, but methodologies to fingerprint gas sources are as yet weak. The production of unconventional gas typically requires the joint extraction of very large quantities of groundwater per day, which can affect groundwater levels in neighbouring aquifers1'2. Another environmental concern is the impact of gas migration on shallow groundwater resources1-5. Recent studies around shale gas production sites in the USA report elevated CH4 in aquifers up to two kilometres away from the production wells6,7. More recent studies have argued that this may not be a result of shale gas production8,9. Although research is largely associated with shale gas developments in the USA1-9 similar issues associated with aquifer hydraulic connectivity also accompany CSG developments worldwide. It is clear that a methodology to quantify potential impacts of unconventional gas production needs to be developed, as well as an understanding of existing pathways of hydraulic connectivity, prior to significant gas production.

■"■Connected Water Initiative Research Centre, UNSW Australia, UNSW Sydney, NSW, 2052, Australia. 2Australian Nuclear Science and Technology Organisation, New Illawarra Rd, Lucas Heights, NSW, 2234, Australia. 3Royal Holloway, University of London, Egham Hill, Egham, Surrey TW20 0EX, United Kingdom. 4School of Environmental Sciences, University of East Anglia, Norwich, Norfolk, NR4 7TJ, United Kingdom. Correspondence and requests for materials should be addressed to C.P.I. (email:







Figure 1. Site map showing the extent of the region investigated within the Condamine Catchment, south-east Queensland. The shaded polygon delineates the extent of the irrigation district53. The location of the CSG production and exploration wells are shown as blue54 circles and the sampled irrigation bores are shown as white circles. The red line highlights the continuous mobile CH4 survey route around the irrigation district. The black square shows the map area in Fig. 3a (Map created in QGIS; data and imagery: MapQuest, Open Street Map and contributors, CC-BY-SA55. Modified with Corel Painter 201556).

This research tests the hypothesis that CH4 can be used to map zones and identify pathways of hydraulic connectivity between a gas production site and an adjacent freshwater aquifer. Natural pathways of connectivity may be via faults, fracture networks and permeable sedimentary formations10. There are also potential pathways of hydraulic connectivity via abandoned exploration wells and faulty well casings11,12.

Previous studies in the USA near shale gas production have made use of stable isotopes to identify sources of CH4 in shallow groundwater13,14. This study, for the first time, tests the suitability of jointly using 613C-CH4, DOC concentration [DOC] and 3H activity in the groundwater to assess hydraulic connectivity. These three measurements provide considerable insights into pathways of groundwater and gas movement because:

- 3 H activities provide information on groundwater residence times and recharge pathways,

- [DOC] provides a measure of the carbon inputs, either from the river recharge or the upward migration of CH4 from a coal bed, and

- 613C-CH4 can be used to characterise the potential sources of the CH4 within an aquifer.

The suitability of using 3H, [DOC] and 613C-CH4 as key parameters for identifying locations where there is hydraulic connectivity between a freshwater aquifer and an underlying unconventional gas reservoir is described in detail below. The methodology is applied to a hydraulic connectivity investigation between the target WCM (on an international scale, a large unconventional gas development) and the overlying CRAA, which supplies water for irrigated agricultural regions, producing a billion dollars worth of commodities. The methodologies presented here are applicable to the many large and geologically similar CSG resources worldwide, including those from the USA, India, China and South Africa15, because we measure geochemical parameters that are common to all sedimentary systems.

Study Area

In Australia, CSG exploration and production have expanded rapidly throughout south-east Qld over the last decade (Fig. 1)15. This has generated considerable public concern about the impact on adjacent aquifers used to support irrigated agriculture, stock and domestic water supplies. The proximity of gas extraction to aquifers used for irrigation or domestic water supply is common to many CSG production sites globally. In this study we investigate groundwater from the unconfined CRAA, which has

Stratigraphy Index

Condamine River Alluvial Aquifer

I Main Range Volcanics

□ Kumbarilla Beds

■200 m

® Walloon Coal Measures

Marburg/Hutton Sandstone

Base (Great Artesian Basin)

Figure 2. Geological cross section along A-A' in Fig. 1 (adapted from KCB Final Report26).

historically supplied 81.4 GL/year of water for irrigated agriculture16. This aquifer is located immediately east of the expanding Daandine and Kogan North CSG developments, which began producing in 200617. Currently, within 50 km of the western boundary of the CRAA, there are around 500 producing wells (Fig. 1). The target gas resource for production wells (WCM) starts at approximately 300 m below the ground surface18.

There are insufficient baseline CH4 data of near-surface ambient air and dissolved gas in groundwater within the Condamine Catchment19-21, and our understanding of the extent of hydraulic connectivity between the WCM and the CRAA is limited. Owen et at.20. used principal component analysis (PCA) to investigate connectivity using historical geochemical major ion data from 879 wells within the CRAA and adjacent Surat Basin. They concluded "in general, no relationships were observed between CSG groundwater in the WCM and the alluvial groundwater.20" Major ions migrate between aquifers via advection and there would need to be significant movement of water between formations for geochemical mixing to be apparent. They did not examine [DOC], CH4 concentration [CH4] or isotopic composition (613C-CH4) within the groundwater. Another study in the region, using 3H, dissolved inorganic carbon isotopes (613C-DIC) and Sr isotopes (87Sr/86Sr) found that there were limitations in using conventional tracers to assess aquifer connectivity22. The authors concluded that "other innovative tracers need to be tested in order to confirm inter-aquifer interactions in CSG environments"22.

Hydrogeological Setting. The entire study area sits within the Surat Basin, which is part of the Great Artesian Basin (GAB) in south-east Qld (Fig. 1). The units of the GAB, including the WCM, vary locally between semi-confined and confined23 and the WCM in places immediately underlie the CRAA

The environment of deposition for the Surat Basin was mainly fluvio-lacustrine during the late Triassic-Jurassic and shallow marine and coastal in the Cretaceous, similar to associated GAB units25. The middle Jurassic WCM are a low-rank CSG resource in the Surat Basin25. They consist of very fine- to medium-grained sandstone, siltstone, mudstone and coal, with minor calcareous sandstone, impure limestone and ironstone26. The coal consists of numerous discontinuous thin lenses separated by sediments of low per-meability27. The WCM are up to 700 m thick, however the thickness of the coal makes up less than 10% of the total thickness of the unit. The WCM dip gently (~4°) to the west, consistent with the general trend of the Surat Basin in this region. In the region of study, the WCM is thicker (150 m to 350 m) along the western margin of the CRAA, and thins to approximately 50 m in the east, where it outcrops26. Also, along the eastern margin of the valley, the CRAA is bounded by the Main Range Volcanics24. The WCM overly the Eurombah Formation (a conglomeratic sandstone with minor siltstones and mudstone beds) and underlies the Kumbarilla beds (predominantly sandstone, with lesser mudstone, siltstones and conglomerates)26. The unconfined CRAA fills a palaeovalley that was carved through the GAB formations. These reworked eolian and alluvial sediments were deposited throughout the mid-Miocene to the present23,28. The valley-filling sediments, which have a maximum thickness of 134 m near Dalby23,24, consist of gravels and fine- to coarse-grained channel sands interbedded with floodplain clays and, on the margins, colluvial deposits.

Sources of Methane

Within the Condamine Catchment there are many potential sources of CH4 that could contribute to the measured [CH4] and 613C-CH4 in the groundwater and air. These include CH4 from:

- the upward migration of gas from the WCM,

- biological activity in the saturated zone beneath rivers and wetlands,

- biological activity in the saturated zone of the CRAA,

- biological activity in the vadose zone (CH4 sink),

- biological activity within the borehole, and

- anthropogenic inputs into the atmosphere.

(Fig. 2)24.

Each of the above CH4 sources potentially has a unique 613C-CH4 signature29. Coal commonly has a portion of the CH4 derived from thermogenic origins and therefore contains more 13C, resulting in a less negative 613C-CH4 signature. In contrast, all the listed biological sources of CH4 will be depleted in 13C relative to thermogenic sources and atmospheric background, resulting in a more negative 613C-CH4 signature11,30. However, there can be complications in interpreting the isotopes, as various processes can alter the concentration or isotopic composition of the gas. Oxidation of CH4 in the vadose zone and atmosphere31, groundwater flow fractionation within boreholes32, and the mixing of various sources of CH433 can alter the isotopic signature. Therefore, improved knowledge of the sources of CH4 is required if we are to attribute measured changes in concentration in the groundwater and air to the appropriate sources.

Biological activity in saturated zones produces CH4 as a result of methanogenesis. This can occur via either acetate fermentation or reduction of carbon dioxide (CO2) (equation (1) and (2) respectively):31

CH3 COOH ^ CH4 + CO2 (1)

CO2 + 8H + + 8e- ^ CH4 + 2H2 O (2)

Anaerobic oxidation of CH4 (AOM), via reduction of sulphate (SO42—) or denitrification of nitrate (NO3—) may also occur in the saturated zone34-36. One of the conditions required for AOM is elevated concentrations of SO42— or NO3— in the groundwater. Therefore, testing the groundwater for SO42— or NO3— can give an indication as to the potential for AOM to occur. Microbes in the deep aquifer are also responsible for AOM37. However, no member of the methanotrophic Euryarchaeota group responsible for anaerobic oxidation has been cultured yet and the various pathways by which AOM occurs are still largely unknown37.

Biological activity within the vadose zone yields a net consumption (oxidation) of CH4. Oxidation of CH4 results in an enrichment of 13C, as the methanotrophs preferentially consume 12C31. This process can go through different reaction intermediates before producing CO2 (equation (3)):31

CH4 ^ CH3OH, CHOH, HCOO- ^ CO2 (3)

Methanotrophic processes can mask the full extent of the upward migration of CH4 because they deplete the concentration of CH4 in the subsurface. Oxidation processes in the near surface can also alter the source isotopic signature via fractionation. This allows us to detect methanotrophic activity in a 12C rich environment, such as the vadose zone, because the isotopic signature of CH4 will be depleted in 12C. This is evident in the results presented below.

Within the WCM most of the CH4 is adsorbed to the coal, however a small portion can exist in a free state15. Methane in its free state is buoyant, and over time may migrate upwards. Since the 1960s groundwater levels within the CRAA have declined24 by 25 m and this may have enhanced the upwards migration of CH4. The reduction in weight immediately overlying the WCM has the potential to unload the coal interval, releasing a portion of the free state CH4. This CH4 will have an isotopic signature that is different from near-surface microbiologically produced CH4.

Using Isotopes of Carbon and Hydrogen to Attribute the Source of Methane

To assess the potential of CH4 as a marker for hydraulic connectivity, we measured [CH4] and 613C-CH4 in both the groundwater and air of the Condamine Catchment. The [DOC] was measured in each irrigation borehole to investigate pathways of groundwater flow. Tritium activity was measured in the groundwater to provide insights about groundwater residence times, recharge processes and potential groundwater mixing. Anions (SO42— and NO3—) were also measured to provide insights into anaerobic processes deep in the aquifer and 613C of dissolved inorganic carbon (613C-DIC) was measured to give insights into the source of the inorganic carbon in the groundwater.

Atmospheric background concentrations of CH4 change with latitude due to changing sources and sinks38. The background concentration of CH4 at the Cape Grim Baseline Air Pollution Station in Tasmania, Australia (40.683°S), was 1.754 ± 0.002 ppm in 201339. Continuous measurements from the Commonwealth Scientific and Industrial Research Organisation (CSIRO) and Geoscience Australia station at Arcturus, central Queensland (24.027°S) currently range between 1.750 ppm and 1.770 ppm40. The isotopic signature value (613C-CH4) commonly reported for ambient background air is — 47%o41.

Accepted threshold values and isotopic ranges allow us to attribute a source from the measured 613C-CH4 obtained for each source of CH4 encountered throughout the region41. Australian CSG has a 613C-CH4 value of less than -60% for biologically sourced gas (biogenic CH4) and greater than —50% for gas sourced from coal at depth (thermogenic CH4)42. The intermediate zone is classified as being of mixed source. Previous studies have measured the CH4 from the WCM using core samples from gas wells. The 613C-CH4 readings for these cored samples ranged from —58.5% to — 45.3%25,43,44, which indicates that the gas within the WCM is secondary biogenic in origin, with a minor thermogenic component25.

150°55'E 150°57'E 150°59'E


Figure 3. (a) A map of CSG developments and co-produced water storage reservoirs. The blue rectangles represent the location of the water storage reservoirs and the red line indicates the route of the mobile survey. The A-B represents the mobile CH4 measurement traverse lines, where Tedlar bag samples were collected downwind of the co-produced water reservoir to establish the isotopic composition of CH4 extracted from the WCM (Map created in QGIS; data and imagery: MapQuest, Open Street Map and contributors, CC-BY-SA55. Modified with Corel Painter 201556). (b) CH4 concentration measured along the line A-B in the plume downwind of the co-produced water storage reservoir. The horizontal red line indicates the background CH4 concentration of 1.774 ppm.

Results and Discussion

Background ambient air. We established the background [CH4] and isotopic signature of air in rural New South Wales (NSW) and Qld by collecting 8 ambient air samples in areas of native vegetation and dryland farming between Narrabri (NSW) and Dalby (Qld). The average background [CH4] was 1.774 ± 0.002 ppm and the isotopic signature was -47.0 ± 0.05%o (n = 8; Supplementary Table S3 online). These values are slightly higher than those reported for the Cape Grim Baseline Air Pollution Station in Tasmania, and at Arcturus in central Qld39,40.

Co-produced water storage reservoir. We were unable to measure [CH4] directly from CSG well-heads; however, at each production well large quantities of water are co-produced. This water is stripped of most of the CH4 and then stored in large holding reservoirs. We were able to indirectly measure the gas from the WCM by analysing the emissions from one such reservoir with dimensions 800 m by 500 m. To determine the signature of the gas being produced from the WCM we surveyed immediately adjacent to and downwind (east) of the reservoir (Fig. 3a). At the closest point the measurements were within 50 m of the reservoir. We detected an area of elevated [CH4] that had a peak of 2.107 ppm and a width of 2.3 km.

The [CH4] measured during the near-reservoir traverses are plotted in Fig. 3b. Each coloured curtain represents a single survey run in the vehicle. For each run, there are slight differences in the recorded [CH4] at the same location. We attribute this to the varied mixing processes associated with the fluctuating winds (both direction and speed).

The isotopic composition of a mixed air sample, which consists of a point source (613C-CH4 (s); [CH4 (s)]) added to background air (613C-CH4 (b); [CH4 (b)]), is determined from:45-47

0.2 0.3 0.4

1/[CH4] (ppm)

Figure 4. Co-produced water reservoir downwind plume mixing plot. Background air is shown as a red dot, and the downwind plume samples are the black dots. The blue line is the line of best fit for the mixing model (equation (4)) and the 90% confidence interval bands are shown in orange. Source intercept determined for the mixing model is -50.8%o (90% CI, -55.7%o to -45.8%o).

where 613C-CH4 (a) and [CH4 (a)] are the isotopic value and concentration measured in the ambient air sample, respectively.

For 613C-CH4 (b) and [CH4 (b)] we used the average of the air samples collected in areas of natural vegetation for the background values (613C-CH4 (b) = -47.0%; [CH4 (b)] = 1.774 ppm). The isotopic signature of the off-gassing co-produced water reservoir was then determined by fitting equation (4) to the data using the nonlinear model fitting function in Mathematica48 (Fig. 4; Supplementary Table S6 online). To determine the source 613C-CH4 value for the WCM three air samples collected along line A-B (Fig. 3), when the wind was blowing from the west directly across the reservoir, were combined with two other samples downwind of the CSG region. Figure 4 shows the mixing model (equation (4)) line of best fit for the combined samples. This yielded a value of -50.8% (90% CI, -55.7% to -45.8%) for 613C-CH4 (s). The isotopic signature of the CH4 off-gassing from the reservoir indicates that the gas has mixed thermogenic/biogenic origin. The broad confidence interval (-55.7% to -45.8%) is due to a combination of the sample size and measurement error. These 90% confidence interval bounds sit within the 613C-CH4 range reported for the WCM43. We therefore use the isotopic signature of - 50.8% as our reference value for the WCM. In the discussion below we use this isotopic signature for the WCM as part of attributing the source of CH4 within the aquifers.

The irrigation district. The irrigation district, to the south of Dalby and the east of Cecil Plains, has a zone of CH4 production along the Condamine River, with [CH4] peaks as high as 1.930 ppm. However, most of the region is undergoing methanotrophic processes, resulting in a slight decrease in [CH4] in the irrigation district. In and around the irrigation district, the continuous CH4 survey (Fig. 1) recorded low [CH4] values, with the lowest concentration measured being 1.764 ppm. The average [CH4] recorded in the continuous survey throughout the district was 1.771 ± 0.007 ppm (n = 8954), slightly lower than the regional atmospheric background (613C-CH4 (b)=-47.0%; [CH4 (b)] = 1.774 ppm).

We were able to further examine the processes occurring in the unsaturated and saturated zones of the aquifer by analysing the [CH4] and isotopic signatures within irrigation boreholes in the region.

Irrigation boreholes. The second data set focused on the groundwater and gas within private irrigation boreholes, details of which are provided in Supplementary Table S4 online.

Our 613C-DIC and 613C-CH4 measurements from all irrigation boreholes display no correlation (R2 = 0.04), with the values falling outside of the methanogenesis or sulfate reduction zones7 (Fig. 5; Supplementary Table S5 online). The DIC in this groundwater is in part from the dissolution of rego-lith carbonates formed in the vertosols. These calcareous soils are often distributed close to areas with Ca-rich bedrock, such as limestone and basalts, such as the Main Range Volcanics in the study area49. Therefore the lack of a correlation between 613C-DIC and 613C-CH4 implies that the CH4 is from a different source. In other studies microbial CH4 from shallow depths has a strong positive correlation between 613C-DIC and 613C-CH47. Results obtained for 613C-DIC are 18% to 22% lighter (more negative) than values expected for DIC sourced from microbial CH4 in shallow groundwater.

Figure 5. 513C-DIC against 613C-CH4 showing no correlation (R2 = 0.04).

In this study we use the absence of 3H and [DOC] above the detection limit to provide insights into hydraulic connectivity because this combination implies that, in the old groundwater (>70 years), there is another source of [DOC]. Recorded 3H activities in the groundwater from irrigation boreholes were generally higher in the Condamine River corridor, and lower compared to modern rainfall values in the central portion and east of the catchment.

To interpret the CH4 measurements from the boreholes we plotted 613C-CH4 versus 1/[ CH4] (Fig. 6; Supplementary Table S5 online (Fig. 6b highlights the cluster of irrigation borehole samples around 0.55 ppm-1 in Fig. 6a.)). Two distinct distributions of points are shown in Fig. 6. These relate to the measured activity of 3H, [DOC] > 0.1 mg/L, [CH4] and 613C-CH4. We assigned samples to subsets for fitting equation (4) to the data (allowing the background values to float) based on the presence or absence of 3H activity and [DOC]. Only samples 9, 16, 17 and 19 have 3H activity below the quantification limit and [DOC] above the detection limit, and these samples were assigned to subset A. The remaining samples have a mixture of 3H activities and detectable [DOC] values, and were assigned to subset B.

Under the assumption that there are only two end members, an estimate of their isotopic signatures can be obtained from a least squares regression fit of equation (4). Subset A has a /-axis intercept for the regression line of — 55.9%o (90%o CI, — 58.3%o to — 53.4%o). This is the source signature of the CH4 in these boreholes. The isotopic value falls within the range reported for gas from the WCM and is also similar to our recorded value for the WCM from the co-produced water reservoir. The other end member indicated from the line of best fit has a value of — 40.7% (90% CI, — 40.75% to —40.74%), which is consistent with mixing from an oxidised source (for example the vadose zone). The absence of 3H activity in the subset A samples indicates that recharge from the river corridor takes > 70 years to reach these locations. The elevated [DOC] in subset A cannot be attributed to recharge. This is because in the time that it takes for the recharge to reach these boreholes (> 70 years), it is reasonable to assume that biological processes would have already consumed the available DOC. Therefore another source is supplying the detected DOC. Upwards migration of CH4 from the WCM would be the most likely source.

Subset A boreholes all had elevated concentrations of DO; including the highest three values recorded in our study (Supplementary Table S5 online). Microbial activity would consume this oxygen50. We therefore infer that there is minimal biological activity occurring at depth within the aquifer at these locations. In addition, the SO42— and NO3— concentrations in the groundwater across the irrigation boreholes were too low to support AOM (Supplementary Table S5 online). Anaerobic oxidation requires one of these pathways to occur and therefore it is reasonable to assume that oxidation deep within the aquifer is minimal.

In summary for subset A, the absence of 3H activity in irrigation boreholes with detectable [DOC], in association with the isotopic value of the regression line within the range of the WCM, suggests local connectivity between the CRAA and the WCM in these areas. The detected CH4 is attributed to the upward migration of gas from the WCM.

Subset B has a y-axis intercept for the regression line of — 69.1% (90% CI, — 73.2% to —65.0%), which indicates a biological source. The regression line indicates a second end member with an isotopic signature of —29.3% (90% CI, —29.33% to — 29.32%), which represents mixing with a highly oxidised background (the vadose zone). All samples located in a cluster at 0.55 ppm—1 have [DOC] below the detection limit, with the exception of irrigation borehole 8, despite many of these irrigation boreholes being near the river. However, the 3H activity indicates that samples 3, 5, 8, and 18 are all sourcing water less than 70 years old, which indicates that these irrigation boreholes are extracting water influenced by near surface processes. Samples 1, 2, 4, 6, 7, 10, 11, 12, and 14 have no recorded 3H activity and the [DOC] is below the detection threshold. This indicates that these samples are accessing water > 70 years old and the DOC has already been utilised by biological processes.

0.4 0.6 0.8 1.0 1.2

1/[CH4] (ppm)

1/[CH4] (ppm)

: Tritium, [DOC] < Detection Limit s Tritium, [DOC] > Detection Limit No Tritium, [DOC] < Detection Limit # No Tritium , [DOC] > Detection Limit

Figure 6. (a) A combined mixing plot for the irrigation bores. The regression lines represent either the bores with detectable [DOC] (intercept: -55.9% (90% CI, - 58.3% to -53.4%)) or the bores with no detectable [DOC] (intercept: - 69.1% (90% CI, -73.2% to - 65.0%)). Tritium activity is indicated by the marker color. All boreholes with no 3H activity and detectable [DOC] sit on the regression line with intercept of -55.9%, indicating gas from the WCM. (b) The cluster of bores around 0.55 ppm-1 highlights a mixing of parameters at the intercept of the regression lines. This indicates that the absence of 3H activity can be used to attribute source but it must be used in tandem with detectable [DOC] values. Mixing model lines of best fit are shown in blue, and the 90% confidence interval bands in orange.

We attribute the source signature of these samples to biological activity within and near the irrigation boreholes. The low [DOC] and DO measurements in the presence of biologically produced CH4 is possibly due to microbiological activity after recharge consuming both the DOC and DO. This biological activity would have produced CH4 as a by-product, and this source of CH4 dominates the isotopic signature measured in the groundwater at these locations50.

Measuring 613C-CH4 and [CH4] is not enough to assign source, there is also a need to measure both the 3H activity and DOC, which is highlighted by samples 5 and 17. Samples 5 and 17 both had [DOC] above the detection limit, however sample 17 had no recorded 3H activity, whereas sample 5 did. Sample 13 had no detectable [DOC] and no 3H activity measured in the groundwater. We propose that biological processes have consumed the DOC. Sample 5 groundwater has been affected by near surface biological activity, whilst sample 17 shows evidence of upward migration of CH4 from the WCM. The difference in these three samples despite their isotopic similarities highlights the need for this combination of measurements in attributing source. The combined method presented here is more robust than the individual measurements because 613C-CH4 is difficult to interpret individually in environments where both thermogenic and biogenic methanogenesis is occurring along with methanotrophy.

Evidence for aquifer connectivity. At our case study site in the Condamine Catchment, the isotopic value of -55.9% from the irrigation boreholes with detectable [DOC] is not as 13C enriched in CH4 as expected from a classical thermogenic source, due to biological processes occurring in situ. However, it falls within the range of the isotopic signature from the WCM and is significantly more enriched in 13C than the samples from irrigation boreholes with [DOC] below the detection limit. Irrigation borehole samples with no 3H activity and detectable [DOC] all sat on the isotopic regression line that fell within

the range of values reported for the WCM (—58.5% to — 45.3%)25,43,44. The value of this regression had a 5.1% difference from our own isotopic signature from the co-produced water reservoir. Thus the isotopic signature from off-gassing air samples collected from irrigation boreholes 9, 16, 17 and 19 (refer Fig. 6), the plume from the co-produced water reservoir, and the literature reported values for gas extracted from coal core samples all sit within the mixed thermogenic/biologic range. This is a strong indicator that the CH4 sampled from irrigation boreholes 9, 16, 17 and 19 is from the underlying WCM.

This research has demonstrated that three conditions need to be present to infer that pathways of geological hydraulic connectivity exist in areas of CSG development. The 613C-CH4 versus inverse [CH4] data must plot on a mixing plot regression line with a y-axis intercept value indicative of thermogenic to mixed thermogenic/biogenic gas. However, at high concentrations (low 1/[CH4]) it is difficult to assign samples to a specific source-determining regression line. Detectable [DOC] can indicate either river corridor recharge or gas migration upwards from depth. Further information about the probable origin of the CH4 and DOC is provided by the 3H activity. Detectable [DOC] along with 3H activity above the quantification limit indicates relatively young water and near-surface biologically sourced CH4. Conversely, [DOC] above the detection limit and 3H activity below the quantification limit strongly suggests CH4 dominated by the upward migration of gas from the WCM. Our results show that measurements of the isotopic composition of CH4, [DOC] and 3H activity in the groundwater and CH4 in the air can be used as an initial assessment of pathways of geological hydraulic connectivity where an alluvial aquifer overlies coal measures targeted for CSG production.


From 20 January 2014 to 1 February 2014 we collected CH4 samples from 19 irrigation boreholes. The irrigation boreholes had a gas outlet, which allowed us to pump the gas directly into 3 L Tedlar bags using an SKC 222-2301 battery-operated air pump. This gas was being stripped from the groundwater within the borehole and is representative of the gas within the groundwater. The groundwater was extracted from 35.4 m to 199.9 m within the aquifer. The depth to the water table ranged from 5 to 20 m51. No purging of these boreholes was required, as the pumps had been running for 2-3 months.

We also sampled groundwater from all 19 irrigation boreholes in the same period. A sampling tube was installed 2 m inside the pump outlet to minimise air/water interactions. Physio-chemical parameters (pH, Eh, temperature, electrical conductivity and dissolved oxygen (DO)) were checked for stability prior to sampling, with all samples fully stable within minutes.

Sample collection involved an in-line, 0.45 |m, high-volume filter, which was connected to the pump outlet. Samples for 613C-DIC and [DOC] were further filtered through a 0.22 |m membrane and collected in 12 mL glass vials (Exetainers) and 60 mL high-density polyethylene (HDPE) bottles, respectively. Samples for 613C-DIC were refrigerated at 4 °C and [DOC] samples were frozen within 12 hours of collection until they were analysed. Samples for 3H analysis were collected in 2 L HDPE bottles. Samples for the SO42—and NO3— analyses were collected in 125 mL HDPE bottles with no further treatment.

A mobile CH4 survey was conducted from 12 March 2014 to 18 March 2014 (Fig. 1). The route covered the irrigation district and CSG developments around Dalby and Cecil Plains to the south, and the location of the irrigation boreholes sampled in January. We sampled the air through Teflon tubing connected to an inlet mounted on top of the vehicle (2.2 m above ground level). This air was passed directly into a Picarro G2301 Cavity Ring-Down Spectrometer (CRDS) (Picarro, Inc., Santa Clara, CA), which measured the [CH4] every 5 seconds as we drove along the road at 40 km/h. The precision of the Picarro G2301 is <0.001 ppm over 5 seconds. Simultaneously, GPS location was measured using a Hemisphere GPS (model A21), whilst wind speed and direction data were measured with an S2 Climatronics 2D sonic anemometer, which was also attached to the top of the vehicle. The setup is a modification of methods previously used to determine gas leakage around Boston, USA52.

Where we discovered an interval of elevated [CH4] individual samples of ambient air were pumped into 3 L Tedlar bags using a KNF NMP 830 KNDC B Micro Diaphragm Gas Pump. The vehicle was stationary and the engine switched off when we pumped air into the Tedlar bags for later analysis in the laboratory (see Supplementary Information online).


1. Vidic, R. D., Brantley, S. L., Vandenbossche, J. M., Yoxtheimer, D. & Abad, J. D. Impact of Shale Gas Development on Regional Water Quality. Science. 340, 1235009-1-9, doi: 10.1126/science.1235009 (2013).

2. Roy, J. W. & Ryan, M. C. Effects of unconventional gas development on groundwater: a call for total dissolved gas pressure field measurements. Groundwater. 51(4), 480-482, doi: 10.1111/gwat.12065 (2013).

3. Vengosh, A., Jackson, R. B., Warner, N., Darrah, T. H. & Kondash, A. A Critical Review of the Risks to Water Resources from Unconventional Shale Gas Development and Hydraulic Fracturing in the United States. Environ. Sci. Technol 48, 8334-8348, doi: 10.1021/es405118y (2014).

4. Moritz, A. et al. Methane Baseline Concentrations and Sources in Shallow Aquifers from the Shale Gas-Prone Region of the St. Lawrence Lowlands (Quebec, Canada). Environ. Sci. Technol. Article ASAP, doi: 10.1021.acs.est.5b00443 (2015).

5. Heilweil, V. M. et al. Stream Measurements Locate Thermogenic Methane Fluxes in Groundwater Discharge in an Area of Shale-Gas Development. Environ. Sci. Technol. Article ASAP, doi: 10.1021/es503882b (2015).

6. Osborn, S. G., Vengosh, A., Warner, N. R. & Jackson, R. B. Methane contamination of drinking water accompanying gas-well drilling and hydraulic fracturing. Proc. Natl. Acad. Sci. USA. 108, 8172-8176, doi: 10.1073/pnas.1100682108 (2011).

7. Jackson, R. B. et al. Increased stray gas abundance in a subset of drinking water wells near Marcellus shale gas extraction. Proc. Natl. Acad. Sci. USA 110(28), 11250-11255, doi: 10.1073/pnas.1221635110 (2013).

8. Lamb, B. K. et al. Direct Measurements Show Decreasing Methane Emissions from Natural Gas Local Distribution Systems in the United States. Environ. Sci. Technol. Article ASAP, doi: 10.1021/es505116p (2015).

9. Siegel, D. I., Azzolina, N. A., Smith, B. J., Perry, A. E. & Bothun, R. L. Methane Concentrations in Water Wells Unrelated to Proximity to Existing Oil and Gas Wells in Northeastern Pennsylvania. Environ. Sci. Technol. Article ASAP, doi: 10.1021/ es505775c (2015).

10. Ward, C. R. & Kelly, B. F. J. Background Paper on New South Wales Geology. University of New South Wales, Sydney, Australia


11. Barker, J. F. & Fritz, P. The occurrence and origin of methane in some groundwater flow systems. Can. J. Earth. Sci. 18, 1802-1816, doi: 10.1139/e81-168 (1981).

12. Fontenot, B. E. et al. An evaluation of water quality in private drinking water wells near natural gas extraction sites in the Barnett Shale Formation. Environ. Sci. Technol. 47, 10032-10040, doi: 10.1021/es4011724 (2013).

13. Hakala, J. A. Use of stable isotopes to identify sources of methane in Appalachian Basin shallow groundwaters: a review. Environ. Sci. Processes Impacts. 16, 2080-2086, doi: 10.1039/c4em00140k (2014).

14. Sharma, S., Mulder, M. L., Sack, A., Schroeder, K. & Hammack, R. Isotope Approach to Assess Hydrologic Connections During Marcellus Shale Drilling. Groundwater. 1-10, doi: 10.1111/gwat.12083 (2013).

15. Moore, T. A. Coalbed methane: a review. Int. J. Coal Geol. 101, 36-81, doi: 10.1016/j.coal.2012.05.011 (2012).

16. Upper Condamine Alluviums groundwater system. Queensland Government, Dept. of Natural Resources and Mines, Toowoomba, Queensland, Australia, 2012.

17. Arrow Energy, Our Operations (2015). Available at: (Accessed: 27th March 2015).

18. Surat Gas Project - Environmental Impact Statement, Arrow Energy Pty Ltd, Queensland, Australia, 2012.

19. Maher, D. T., Santos, I. R. & Tait, D. R. Mapping Methane and Carbon Dioxide Concentrations and 613C Values in the Atmosphere of Two Australian Coal Seam Gas Fields. Water Air Soil Pollut. 225, 2216 (1-9). doi: 10.1007/s11270-014-2216-2 (2014).

20. Owen, D. D. R. & Cox, M. E. Hydrochemical evolution within a large alluvial groundwater resource overlying a shallow coal seam gas reservoir. Sci. Total. Environ. 523, 233-252, doi: 10.1016/j.scitotenv.2015.03.115 (2015).

21. Day, S., Connell, L., Etheridge, D., Norgate, T. & Sherwood, N. Fugitive greenhouse gas emissions from coal seam gas production in Australia. CSIRO, Australia (2012).

22. Duvert, C. et al. Hydrochemical processes in a shallow coal seam gas aquifer and its overlying stream-alluvial system: Implications for recharge and inter-aquifer connectivity. Appl. Geochem. 61, 146-159, doi: 10.1016/j.apgeochem.2015.05.021 (2015).

23. Kelly, B. F. J. & Merrick, N. P. Groundwater Knowledge and Gaps in the Condamine Alliance Area for the Cotton Catchment Communities CRC, University of Technology, Sydney - National Center for Groundwater Management Report, NCGM 2006/9d, 49 pages (2007).

24. Dafny, E. & Silburn, D. M. The hydrogeology of the Condamine River Alluvial Aquifer, Australia: a critical assessment. Hydrogeology Journal. 22, 705-727, doi: 10.1007/s10040-013-1075-z (2014).

25. Hamilton, S. K., Esterle, J. S. & Golding, S. D. Geological interpretation of gas content trends, Walloon Subgroup, eastern Surat Basin, Queensland, Australia. Int. J. Coal Geol. 101, 21-35, doi: 10.1016/j.coal.2012.07.001 (2012).

26. Conceptualisation of the Walloon Coal Measures beneath the Condamine Alluvium - Final Report. KCB (Klohn Crippen Berger), Dept. of Environment and Resource Management, Toowoomba, Queensland, Australia, 2011.

27. Hillier, J. R. Groundwater connections between the Walloon Coal Measures and the Alluvium of the Condamine River. Central Downs Irrigators Limited, Bribie Island, Queensland, Australia, 2010.

28. Kelly, B. F. J. et al. A reassessment of the Lower Namoi Catchment aquifer architecture and hydraulic connectivity with reference to climate drivers. Aust. J. Earth Sci. 61(3), 501-511, doi: 10.1080/08120099.2014.900647 (2014).

29. Fisher, R. E. et al. Arctic methane sources: Isotopic evidence for atmospheric inputs. Geophys. Res. Letters. 38 (21), L21803 (1-6), doi: 10.1029/2011GL049319 (2011).

30. Quay, P. et al. The isotopic composition of atmospheric methane. Global Biogeochem. Cycles. 13, 445-461, doi: 10.1029/1998GB900006 (1999).

31. Whiticar, M. J. Carbon and hydrogen isotope systematics of bacterial formation and oxidation of methane. Chem. Geol. 161, 291-314, doi: 10.1016/S0009-2541(99)00092-3 (1999).

32. Qin, S., Tang, X., Song, Y. & Wang, H. Distribution and fractionation mechanism of stable carbon isotope of coalbed methane. Sci. in China Series D: Earth Sci. 49(12), 1252-1258, doi: 10.1007/s11430-006-2036-3 (2006).

33. Cuna, S. et al. Separating contributions from natural and anthropogenic sources in atmospheric methane from the Black Sea region, Romania. Appl. Geochem. 23, 2871-2879, doi: 10.1016/j.apgeochem.2008.04.019 (2008).

34. Valentine, D. L. & Reeburgh, W. S. New perspectives on anaerobic methane oxidation. Env. Microbiol. 2, 477-484, doi: 10.1046/j.1462-2920.2000.00135.x (2000).

35. Antler, G. Sulfur and oxygen isotope tracing of sulfate driven anaerobic methane oxidation in estuarine sediments. Estuarine, Coastal and Shelf Sci. 142, 4-11, doi: 10.1046/j.1462-2920.2000.00135.x (2014).

36. Green-Saxena, A., Dekas, A. E., Dalleska, N. F. & Orphan, V. J. Nitrate-based niche differentiation by distinct sulfate-reducing bacteria involved in the anaerobic oxidation of methane. The ISME J. 8, 150-163, doi: 10.1038/ismej.2013.147 (2014).

37. Knittel, K. & Boetius, A. Anaerobic Oxidation of Methane: Progress with an Unknown Process. Annu. Rev. Microbiol. 63, 311-334, doi: 10.1146/annurev.micro.61.080706.093130 (2009).

38. Nisbet, E. G., Dlugokencky, E. J. & Bousquet, P. Methane on the rise - again. Science. 343, 493-494, doi: 10.1126/science.1247828


39. Stalker, L. Methane Origins and Behaviour. Commonwealth Scientific and Industrial Research Organisation, Australia, 2013;

40. Integrating a remotely deployed stable isotope Picarro CRDS with a global monitoring network. Centre for Australian Weather and Climate Research, CSIRO Marine and Atmospheric Research, Victoria, Australia, 2013

41. Dlugokencky, E. J., Nisbet, E. G., Fisher, R. & Lowry, D. Global atmospheric methane: budget, changes and dangers. Phil. Trans. R. Soc. A. 369, 2058-2072, doi: 10.1098/rsta.2010.0341 (2011).

42. Kinnon, E. C. P., Golding, S. D., Boreham, C. J., Baublys, K. A. & Esterle, J. S. Stable isotope and water quality analysis of coal bed methane production waters and gases from the Bowen Basin, Australia. Int. J. Coal Geol. 82, 219-231, doi: 10.1016/j. coal.2009.10.014 (2010).

43. Hamilton, S. K., Golding, S. D., Baublys, K. A. & Esterle, J. S. Stable isotopic and molecular composition of desorbed coal seam gases from the Walloon Subgroup, eastern Surat Basin, Australia. Int. J. Coal Geol. 122, 21-36, doi: 10.1016/j.coal.2013.12.003 (2014).

44. Papendick, S. L. et al. Biogenic methane potential for Surat Basin, Queensland coal seams. Int. J. Coal Geol. 88, 123-134, doi: 10.1016/j.coal.2011.09.005 (2011).

45. Keeling, C. D. The concentration and isotopic abundances of atmospheric carbon dioxide in rural areas. Geochem. Cosmochim. Acta. 13, 322-334, doi: 10.1016/0016-7037(58)90033-4 (1958).

46. Keeling, C. D. The concentration and isotopic abundances of carbon dioxide in rural and marine air. Geochem. Cosmochim. Acta. 24, 277-298, doi: 10.1016/0016-7037(61)90023-0 (1961).

47. Pataki, D. E. et al. The application and interpretation of Keeling plots in terrestrial carbon cycle research. Global Biogeochem. Cycles. 17, 1022-1037, doi: 10.1029/2001GB001850 (2003).

48. Wolfram Research, Inc., Mathematica, Version 9.0, Champagne, IL, 2012.

49. Chen, X. Y., Lintern, M. J. & Roach, I. C. Calcrete: characteristics, distribution and use in mineral exploration. Cooperative Research Centre for Landscape Environments and Mineral Exploration, 160 pages (2002).

50. Chapelle, F. H. in Groundwater Microbiology and Geochemistry (John Wiley & Sons, 2001).

51. Central Condamine Alluvium, Stage II - Conceptual hydrogeology study. KCB (Klohn Crippen Berger), Dept. of Environment and Resource Management, Toowoomba, Queensland, Australia, 2010.

52. Phillips, N. G. et al. Mapping urban pipeline leaks: Methane leaks across Boston. Environ. Poll. 173, 1-4, doi: 10.1016/j. envpol.2012.11.003 (2013).

53. Queensland Government Queensland Spatial Catalogue, Land use mapping - South East Queensland NRM region (2015). Available at:{C3CD23AD-7101-4765-8BED-AFC2C5DF1F5F} (First accessed: 24 August 2015).

54. Queensland Government Data, Coal seam gas well locations - Queensland (2013). Available at: queensland-borehole-series/resource/e57de19e-e471-4fb9-b87e-c970d71a4ef0 (First accessed: 2013).

55. QGIS 2.8.2 Wien, 2015, MapQuest and Open Street Map licensed under Creative Commons Attribution - ShareAlike 3.0 license (CC-BY-SA).

56. Corel Corporation, Corel Painter Education, Version (2015).


This research was funded by the Cotton Research and Development Corporation and the National Centre for Groundwater Research and Training (funded by the Australian Research Council and the National Water Commission). We also acknowledge the National Collaborative Research Infrastructure Strategy groundwater funding used to purchase the Picarro G2132-i.

Author Contributions

Experimental conceptualisation and design was carried out by D.I.C., E.G.N. and B.F.J.K. Fieldwork was conducted by C.P.I., D.I.C., S.I.H., R.E.F., D.L., J.L.F. and B.F.J.K. Geochemical analyses were conducted by D.I.C. and S.I.H. Funding for the work was won by B.F.J.K., D.I.C., A.B. and E.G.N. The manuscript was written by C.P.I. and B.F.J.K with input from all authors.

Additional Information

Supplementary information accompanies this paper at

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Iverach, C. P. et al. Assessing Connectivity Between an Overlying Aquifer and a Coal Seam Gas Resource Using Methane Isotopes, Dissolved Organic Carbon and Tritium. Sci. Rep. 5, 15996; doi: 10.1038/srep15996 (2015).

l/gv 0 I This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit