Global Aquifers Dominated by Fossil Groundwaters but Wells Vulnerable to Modern Contamination
Total Page:16
File Type:pdf, Size:1020Kb
ARTICLES PUBLISHED ONLINE: 25 APRIL 2017 | DOI: 10.1038/NGEO2943 Global aquifers dominated by fossil groundwaters but wells vulnerable to modern contamination Scott Jasechko1*, Debra Perrone2,3, Kevin M. Befus4, M. Bayani Cardenas5, Grant Ferguson6, Tom Gleeson7, Elco Luijendijk8, Jerey J. McDonnell9,10,11, Richard G. Taylor12, Yoshihide Wada13,14 and James W. Kirchner15,16,17 The vulnerability of groundwater to contamination is closely related to its age. Groundwaters that infiltrated prior to the Holocene have been documented in many aquifers and are widely assumed to be unaected by modern contamination. However, the global prevalence of these ‘fossil’ groundwaters and their vulnerability to modern-era pollutants remain unclear. Here we analyse groundwater carbon isotope data (12C, 13C, 14C) from 6,455 wells around the globe. We show that fossil groundwaters comprise a large share (42–85%) of total aquifer storage in the upper 1 km of the crust, and the majority of waters pumped from wells deeper than 250 m. However, half of the wells in our study that are dominated by fossil groundwater also contain detectable levels of tritium, indicating the presence of much younger, decadal-age waters and suggesting that contemporary contaminants may be able to reach deep wells that tap fossil aquifers. We conclude that water quality risk should be considered along with sustainable use when managing fossil groundwater resources. lobal groundwater is an immense resource, storing ∼100 (12C, 13C) data were used to correct for the dissolution of carbonate times more water than all the world's lakes1,2, supplying rocks, which are devoid of radiocarbon19 and thus would otherwise G∼40% of the water for global irrigated agriculture3, and distort 14C-based fossil groundwater calculations. We estimated providing drinking water to billions of people around the world. fossil groundwater fractions in wells around the world using a Recent research has evaluated the global depths of both the recently developed radiocarbon endmember mixing model, which groundwater table4 and modern groundwater recharged within accounts for both radioactive decay and carbonate dissolution for the past ∼50 years1, but the global prevalence and distribution pre- and post-Holocene recharge20. Our approach, which estimates of `fossil groundwater' remain unclear. Here we define fossil the fraction of fossil groundwater in a water sample rather than groundwater as groundwater recharged by precipitation more than the sample's average age, is designed to be less vulnerable to the ∼12,000 years ago, prior to the beginning of the Holocene epoch; aggregation errors that are known to bias mean groundwater age we define prevalence as the frequency with which regional well calculations21–23. We plotted depth profiles of fossil groundwater waters contain fossil groundwater. Understanding the global extent for aquifers around the world and calculated the depth below and depth of fossil groundwater resources is important because which fossil groundwater becomes common (>50% of wells pump of their distinctive susceptibility to overdraft5, presumed isolation some fossil groundwater) or dominant (>50% of wells pump from surface-borne pollutants6,7, potential vulnerability to geogenic >50% fossil groundwater; Methods). Where tritium (3H) data contaminants8, and isolation from modern climate variability9. were available (n D 5,661 well water samples), we determined the To calculate the prevalence of fossil groundwater in well waters, fraction of the groundwater sample that recharged more recently we compiled a groundwater carbon isotope (12C, 13C, 14C) database than around 1953 by relating groundwater 3H concentrations to of 6,455 well water samples from around the globe. The continental historical precipitation 3H time series24 (Methods). The threshold USA and Europe are over-represented in our compilation, which year 1953 was selected because widespread thermonuclear testing is also inevitably biased towards sedimentary basins where in subsequent years increased precipitation tritium levels by ∼5 groundwater use is common (Supplementary Figs 1,2 and Methods). to ∼500 times above local natural background concentrations24, Radiocarbon (14C) has a half-life of 5,730 years and has been widely providing a tracer of recently recharged groundwater (for example, used to identify fossil groundwaters10–18. Stable carbon isotope ref. 15). For samples with both radiocarbon and tritium data, we 1Department of Geography, University of Calgary, Calgary, Alberta T2N 4H7, Canada. 2Water in the West, Stanford University, Stanford, California 94305, USA. 3Civil and Environmental Engineering, Stanford University, Stanford, California 94305, USA. 4Civil and Architectural Engineering, University of Wyoming, Laramie, Wyoming 82071, USA. 5Department of Geological Sciences, The University of Texas at Austin, Austin, Texas 78712, USA. 6Department of Civil and Geological Engineering, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5A9, Canada. 7Department of Civil Engineering and School of Earth and Ocean Sciences, University of Victoria, Victoria, British Columbia V8P 5C2, Canada. 8Geoscience Centre, University of Göttingen, Göttingen 37077, Germany. 9Global Institute for Water Security, and School of Environment and Sustainability, University of Saskatchewan, Saskatoon, Saskatchewan S7N 3H5, Canada. 10School of Geosciences, University of Aberdeen, Aberdeen AB24 3FX, UK. 11Department for Forest Engineering, Resources and Management, Oregon State University, Corvallis, Oregon 97330, USA. 12Department of Geography, University College London, London WC1E 6BT, UK. 13International Institute for Applied Systems Analysis, Schlossplatz 1, Laxenburg A-2361, Austria. 14Department of Physical Geography, Utrecht University, Utrecht 80115, The Netherlands. 15Department of Environmental System Sciences, ETH Zürich, Universitätstrasse 16, CH-8092, Switzerland. 16Swiss Federal Research Institute WSL, Birmensdorf CH-8903, Switzerland. 17Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA. *e-mail: [email protected] NATURE GEOSCIENCE | VOL 10 | JUNE 2017 | www.nature.com/naturegeoscience 425 © 2017 Macmillan Publishers Limited, part of Springer Nature. All rights reserved. ARTICLES NATURE GEOSCIENCE DOI: 10.1038/NGEO2943 a TR−Kazan Trona TN−Djerid Basin TH−Bangkok Basin EE−Cambrian−Vendian FR−Bathonian−Bajocian AU−Ngalia−Amadeus AU−Western Port DE−Benkerstein ID−Jakarta Basin CN−Laizhou Bay FR−Aquitaine Basin DZ−Grand Erg Oriental BW−Kalahari Desert CN−Songnen Plain IL−Arava Valley HU−Great Hungarian Plain BD−Bengal Basin IN−Tiruvadanai Aquifer DK−Ribe Formation HU−Pannonian Basin CN−North China Plain FR−Lorraine Sandstone LY−Sirte Basin PL−Malm Limestones PL−Mazovian Basin PL−Southern Carbonates N G−Chad Basin SN−Diass Aquifer IT−Emilia−Romagna Plain JP−Kimitsu−Yoro Aquifer LY−Kufra Basin US−Tampa Bay (FL) US−Tertiary Aquifer (MT) US−Jacksonville Aquifer (FL) US−Aquia−Patapsco (MD) US−Atlantic Coast (VA) US−Castle-Hayne (NC) US−Palouse Basin (WA) US−Surficial Aquifers (FL) US−West Embayment (LA) US−Eureka Coast (CA) US−Mississippi Alluvial (AR) US−Houston Coast (TX) US−Wharton Tract (NJ) US−Floridan Aquifer (GA) US−Central High Plains (KS) US−Bay Area, East (CA) US−Denver Basin (CO) US−Livermore Valley (CA) US−North High Plains (NE) C N−Hexi Corridor US−San Diego Basin (CA) US−Bay Area, Monterey (CA) US−South High Plains (TX) US−South Central Valley (CA) US−Temecula−Murrieta (CA) US−Columbia Basin (WA) US−Cambrian−Ordovician (WI) US−Silurian−Devonian (IL) US−Mahomet Aquifer (IL) US−East Embayment (MS) US−Los Angeles Basin (CA) 0 123 4 5 678 91011 12 13 100 14 15 16 17 18 19 20 21 22 23 24 25 26 27 200 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 300 48 49 50 51 52 53 54 55 400 56 57 Depth below land surface (m) 58 59 60 500 61 62 600 bc2 55 48 17 13 53 26 Aquifer 3 14 35 4 56 Fossil groundwater: 28 39 40 50 30 51 58 20 31 7 43 52 27 57 Present (>0%) in most wells below here 18 37 5 38 36 59 33 60 24 32 47 9 29 22 Dominant (>50%) in most wells below here 49 4419 42 11 46 16 62 45 61 12 21 41 54 23 25 8 34 6 1 15 10 Figure 1 | Prevalence of fossil groundwater in global aquifers. a, Depth to the fossil groundwater transition in 62 aquifers. The shallow depth (top of peach bar) represents a depth below which most wells (>50%) contain detectable fossil groundwater (minimum fossil groundwater fraction >0%). The deeper depth (top of red bar) represents a depth below which most wells (>50%) are dominated by fossil groundwater (minimum fossil groundwater fraction is >50%). Fossil groundwater becomes dominant at a median depth of 200 m (lower–upper quartiles are 115–290 m). The lower limit of our graph (600 m) does not necessarily represent the lower boundary for our study aquifers, nor do the depths covered by red and orange bars imply that groundwater is of high quality. b,c, Locations of the 62 aquifers. calculated the three fractions of groundwater that recharged more North China Plain (at depths >200 m), the southern Central Valley than ∼12 thousand years ago (fossil groundwater); recharged of California (at depths >260 m), the north, central and south High more recently than the year 1953 (post-1953 groundwater); and Plains aquifers of the central USA (at depths >120–280 m), Italy's is of an intermediate age, having recharged more recently than Emilia–Romagna Plain (at depths >100–300 m) and Hungary's ∼12 thousand years ago, but before the year 1953. Pannonian Basin (at depths >160–300 m). Among our 62 study aquifers (Fig.1), we find the range of depths below which fossil Fossil groundwater in global aquifers groundwaters dominate well waters (that is, fossil groundwaters Our global compilation of radiocarbon data shows that fossil comprise >50% of the water pumped from more than half of all groundwater is not an anomaly in the upper 1 km of the crust, deeper samples) has a median of 200 m, a lower–upper quartile but instead is common in wells drilled to depths of more than range of 115–290 m, and a 10th–90th percentile range of 70–430 m.