Resources

1.
Bain, R. et al. Fecal contamination of drinking-water in low- and middle-income countries: a systematic review and meta-analysis. PLoS Med. 11, e1001644 (2014).
Article  Google Scholar 
2.
Charles, K. & Pond, K. Drinking water quality regulations. Routledge Handb. Water Health https://doi.org/10.4324/9781315693606 (2015).
Article  Google Scholar 

3.
Rahman, Z., Crocker, J., Chang, K., Khush, R. & Bartram, J. A comparative assessment of institutional frameworks for managing drinking water quality. J. Water Sanit. Hyg. Dev. 1, 242–258 (2011).
Article  Google Scholar 

4.
Ashbolt, N. J., Grabow, W. O. K. & Snozzi, M. Indicators of microbial water quality. in Water Quality: Guidelines, Standards and Health: Assessment of Risk and Risk Management for Water-Related Infectious Disease (eds Fewtrell, L. & Bartram, J. K.) (World Health Organization, London, UK, 2001).

5.
HMSO. The Bacteriological Examination of Water Supplies. Reports on Public Health and Medical Subjects, 1st edn, No. 71 (1934).

6.
Leclerc, H., Mossel, D. A. A., Edberg, S. C. & Struijk, C. B. Advances in the Bacteriology of the Coliform Group: Their Suitability as Markers of Microbial Water Safety. (2001).

7.
Bonde, G. J. Bacteriological methods for estimation of water pollution. Health Lab. Sci. 3, 124–128 (1966).
CAS  Google Scholar 

8.
Gordon, D. M. The ecology of Escherichia coli. in Escherichia coli: Pathotypes and Principles of Pathogenesis (ed Donnenberg, M.S.) 3–20 (Academic Press, 2013).

9.
Medema, G. J. et al. Safe drinking water: an ongoing challenge. in Assessing Microbial Safety of Drinking Water: Improving Approaches and Methods (eds O. E. C. D. & WHO) (2003).

10.
WHO. Guidelines for Drinking-water Quality: Third edition, incorporating the first and second Addenda, Volume 1 Recommendations. (World Health Organisation, 2008).

11.
WHO. Guidelines for Drinking Water Quality 2nd Edition. Volume 1. Recommendations. (World Health Organisation, 1993).

12.
Anon. Council Directive 98/83/EC of 3 November 1998 on the quality of water intended for human consumption. Off. J. Eur. Commun. 5.12.98, L330/32–L330/53 (1998).

13.
Gray, N. & Gleeson, C. The Coliform Index and Waterborne Disease. https://doi.org/10.4324/9780203476888 (1996).

14.
Hamilton, W. P., Kim, M. & Thackston, E. L. Comparison of commercially available Escherichia coli enumeration tests: implications for attaining water quality standards. Water Res. 39, 4869–4878 (2005).
CAS  Article  Google Scholar 

15.
WHO & UNICEF. Rapid Assessment of Drinking-Water Quality: A Handbook for Implementation. (2012).

16.
Mahmud, Z. H. et al. Occurrence of Escherichia coli and faecal coliforms in drinking water at source and household point-of-use in Rohingya camps, Bangladesh. Gut Pathog. 11, 1–11 (2019).
Article  Google Scholar 

17.
Allen, M. A. The Public Health Significance of Bacterial Indicators in Drinking Water. in Coliforms and E. coli, Problem or Solution? (eds Kay, D. & Fricker, C.) 176–181 (Royal Society of Chemistry, 1997).

18.
Geldreich, E. E. Reinventing microbial regulations for safer water supplies. in Coliforms and E. coli, Problem or Solution? (eds Kay, D. & Fricker, C.) 218–234 (Royal Society of Chemistry, 1997).

19.
WHO. Guidelines for Drinking Water Quality. First Edition. Vol 1 Recommendations. (World Health Organisation, 1984).

20.
Pepper, D. Bangladeshis poisoned by arsenic sue British organisation. Lancet 367, 199–200 (2006).
Article  Google Scholar 

21.
George, C. M. The arsenic crisis in Bangladesh. in Routledge Handbook of Water and Health (ed Bartram, J.) 685–692 https://doi.org/10.4324/9781315693606-82 (Routledge, 2015).

22.
WHO. Arsenic in Drinking-water Background Document for Development of WHO Guidelines for Drinking-Water Quality. (World Health Organisation, 2011).

23.
WHO. Guidelines for Drinking-water Quality 4th ed., WHO, Geneva, p. 340. World Health Organization https://doi.org/10.1016/S1462-0758(00)00006-6 (2011).

24.
Anon. Water, England and Wales The Water Supply (Water Quality) Regulations SI 2016/614. Water, England and Wales 52 (2016).

25.
Waite, W. M. Drinking water quality regulation—a European perspective. in Coliforms and E. coli, Problem or Solution? (eds Kay, D. & Fricker, C.) 208–217 (Royal Chemistry Society, 1997).

26.
Geldreich, E. E. Coliforms: a new beginning to an old problem. in Coliforms and E. coli, Problem or Solution? (eds Kay, D. & Fricker, C.) 3–11 (Royal Society of Chemistry, 1997).

27.
Mossel, D. A. A. Index and indicator organisms: a current assessment of their usefulness and significance. Food Technol. Aust. 30, 212–219 (1978).
Google Scholar 

28.
Wright, J. et al. Water quality laboratories in Colombia: a GIS-based study of urban and rural accessibility. Sci. Total Environ. 485–486, 643–652 (2014).
Article  Google Scholar 

29.
Brennan, F. P., O’Flaherty, V., Kramers, G., Grant, J. & Richards, K. G. Long-term persistence and leaching of Escherichia coli in temperate maritime soils. Appl. Environ. Microbiol. 76, 1449–1455 (2010).
CAS  Article  Google Scholar 

30.
Brennan, F. P., Abram, F., Chinalia, F. A., Richards, K. G. & O’Flaherty, V. Characterization of environmentally persistent Escherichia coli isolates leached from an Irish soil. Appl. Environ. Microbiol. 76, 2175–2180 (2010).
CAS  Article  Google Scholar 

31.
Haller, L., Poté, J., Loizeau, J. L. & Wildi, W. Distribution and survival of faecal indicator bacteria in the sediments of the Bay of Vidy, Lake Geneva, Switzerland. Ecol. Indic. https://doi.org/10.1016/j.ecolind.2008.08.001 (2009).

32.
Pote, J. et al. Persistence and growth of faecal culturable bacterial indicators in water column and sediments of Vidy Bay, Lake Geneva, Switzerland. J. Environ. Sci. 21, 62–69 (2009).
CAS  Article  Google Scholar 

33.
Ashbolt, N. J., Dorsch, M. R., Cox, P. T. & Banens, B. Blooming E. Coli, what do they mean? in Coliforms and E. coli, Problem or Solution? (eds Kay, D. & Fricker, C.) (The Royal Society of Chemistry, Cambridge, England, 1997).

34.
Block, J. C., Muteaux, L., Gatel, D. & Reasoner, D. J. Survival and growth of E. coli in drinking water distribution systems. in Coliforms and E. coli, Problem or Solution? (eds Kay, D. & Fricker, C.) 157–167 (Royal Society of Chemistry, 1997).

35.
Ferguson, A. S. et al. Hand-pumps as reservoirs for microbial contamination of well water. J. Water Health 9, 708–717 (2011).
CAS  Article  Google Scholar 

36.
LeChavallier, M. W. Conditions favouring coliform and HPC bacterial growth in drinking-water and on water contact surfaces. in Heterotrophic Plate Counts and Drinking-water Safety: The Significance of HPCs for Water Quality and Human Health (eds Bartram, J., Cotruvo, J., Exner, M., Fricker, C. & Glasmacher, A.) 177–198 (WHO, 2003).

37.
Benz, S. A., Bayer, P. & Blum, P. Global patterns of shallow groundwater temperatures. Environ. Res. Lett. 12, 34005 (2017).
Article  Google Scholar 

38.
McDonough, L. K. et al. Changes in global groundwater organic carbon driven by climate change and urbanization. Nat. Commun. 11, 1–10 (2020).
Article  Google Scholar 

39.
Leclerc, H., Mossel, D. A. A., Edberg, S. C. & Struijk, C. B. Advances in the bacteriology of the coliform group: their suitability as markers of microbial water safety. Annu. Rev. Microbiol. 55, 201–234 (2001).
CAS  Article  Google Scholar 

40.
Osborn, M. J. et al. Indicators for Waterborne Pathogens. (The National Academies Press, 2004).

41.
Murphy, H. Persistence of Pathogens in Sewage and Other Water Types. in Global Water Pathogen Project (eds Rose, J. & Jiménez-Cisneros, B.) https://doi.org/10.14321/waterpathogens.51 (2017).

42.
Taylor, R., Cronin, A., Pedley, S., Barker, J. & Atkinson, T. The implications of groundwater velocity variations on microbial transport and wellhead protection—review of field evidence. FEMS Microbiol. Ecol. 49, 17–26 (2004).
CAS  Article  Google Scholar 

43.
Payment, P. & Locas, A. Pathogens in water: value and limits of correlation with microbial indicators. Ground Water 49, 4–11 (2011).
CAS  Article  Google Scholar 

44.
Bivins, A. W. et al. Estimating infection risks and the global burden of diarrheal disease attributable to intermittent water supply using QMRA. Environ. Sci. Technol. 51, 7542–7551 (2017).
CAS  Article  Google Scholar 

45.
WHO. Surveillance of Drinking Water Quality. (WHO, 1976).

46.
Taylor, E. W., Beale, J. F., Suckling, E. V. & Thresh, J. C. The Examination of Waters and Water Supplies. (Churchill, 1958).

47.
Setty, K. E. et al. Time series study of weather, water quality, and acute gastroenteritis at Water Safety Plan implementation sites in France and Spain. Int. J. Hyg. Environ. Health 221, 714–726 (2018).
Article  Google Scholar 

48.
Hunter, P. R., Zmirou-Navier, D. & Hartemann, P. Estimating the impact on health of poor reliability of drinking water interventions in developing countries. Sci. Total Env. 407, 2621–2624 (2009).
CAS  Article  Google Scholar 

49.
Kirchman, D., Sigda, J., Kapuscinski, R. & Mitchell’, R. Statistical analysis of the direct count method for enumerating bacteria. Appl. Environ. Microbiol. 44 (1982).

50.
Setty, K. et al. Assessing operational performance benefits of a Water Safety Plan implemented in Southwestern France. Perspect. Public Health 138, 270–278 (2018).
CAS  Article  Google Scholar 

51.
Prescott, S. C. & Winslow, C.-E. A. Elements of Water Bacteriology, with Special Reference to Sanitary Analysis. (Wiley, 1931).

52.
Kelly, E. R., Cronk, R., Kumpel, E., Howard, G. & Bartram, J. How we assess water safety: a critical review of sanitary inspection and water quality analysis. Sci. Total Environ. 718, 137237 (2020).
CAS  Article  Google Scholar 

53.
Gunnarsdottir, M. J. Benefits of Water Safety Plans: Microbiology, Compliance, and Public Health. https://doi.org/10.1021/es300372h (2012).

54.
WHO. Water Safety Plan Manual: Step-by-step risk management for drinking-water suppliers. WHO Library. https://doi.org/10.1111/j.1752-1688.1970.tb00528.x (2009).

55.
Nowicki, S., Koehler, J. & Charles, K. J. Including water quality monitoring in rural water service provision: why safe water requires challenging the quantity versus quality dichotomy. npj Clean Water 3, 1–9 (2020).
Article  Google Scholar 

56.
Bradley, D. J. & Bartram, J. K. Domestic water and sanitation as water security: monitoring, concepts and strategy. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 371, 20120420 (2013).
Article  Google Scholar 

57.
Anon. Metadata on 6.1.1. (2017). Available at: https://unstats.un.org/sdgs/metadata/files/Metadata-06-01-01.pdf. Accessed 1 May 2020.

58.
WHO. Guidelines for Drinking-Water Quality: Fourth Edition Incorporating the First Addendum. (WHO Press, 2017).

59.
Damania, R. et al. Uncharted Waters. The New Economics of Water Scarcity and Variability. (2017).

60.
Bartram, J., Elliott, M. & Chuang, P. Getting wet, clean, and healthy: why households matter. Lancet 380, 85–86 (2012).
Article  Google Scholar 

61.
UN-WATER. Kenya National Water Development Report: Case Study. A WWAP Case Study Prepared for the 2nd UN World Water Development Report: Water, a Shared Responsibility. (2006).

62.
Mendes, D. M. & Paglietti, L. Kenya Irrigation Market Brief. (FAO and IFC, 2015).

63.
Muller, M. Cape Town’s drought: don’t blame climate change. Nature 559, 174–176 (2018).
CAS  Article  Google Scholar 

64.
Howard, G. et al. Securing 2020 vision for 2030: climate change and ensuring resilience in water and sanitation services. J. Water Clim. Chang. https://doi.org/10.2166/wcc.2010.205 (2010).
Article  Google Scholar 

65.
Taye, M. T., Dyer, E., Hirpa, F. A. & Charles, K. Climate change impact on water resources in the Awash basin, Ethiopia. Water 10, 1–16 (2018).
Article  Google Scholar 

66.
Auld, H., MacIver, D. & Klaassen, J. Heavy rainfall and waterborne disease outbreaks: the walkerton example. J. Toxicol. Environ. Heal. Part A 67, 1879–1887 (2004).
CAS  Article  Google Scholar 

67.
Brown, G. H., Tranter, M. & Sharp, M. J. Experimental investigations of the weathering of suspended sediment by alpine glacial meltwater. Hydrol. Process. 10, 579–597 (1996).
Article  Google Scholar 

68.
Chowdhury, F. J., Ahmad, Z. U. & Aalderink, H. Protecting the Meghna river. A Sustainable Water Resource for Dhaka. (Asian Development Bank, 2019).

69.
UNICEF & WHO. Progress on household drinking water, sanitation and hygiene 2000–2017. Special focus on inequalities. Progress on Drinking Water, Sanitation and Hygiene 2000–2017 (2019). More

1.
Nicholls, R. J. & Cazenave, A. Sea-level rise and its impact on coastal zones. Science 328, 1517–1520 (2010).
CAS  Google Scholar 
2.
Bamber, J. L., Oppenheimer, M., Kopp, R. E., Aspinall, W. P. & Cooke, R. M. Ice sheet contributions to future sea-level rise from structured expert judgment. Proc. Natl Acad. Sci. USA 166, 11195–11200 (2019).
Google Scholar 

3.
Spencer, T. et al. Global coastal wetland change under sea-level rise and related stresses: the DIVA wetland change model. Glob. Planet. Change 139, 15–30 (2016).
Google Scholar 

4.
Moftakhari, H. R. et al. Increased nuisance flooding due to sea-level rise: past and future. Geophys. Res. Lett. 42, 9846–9852 (2015).
Google Scholar 

5.
Vitousek, S. Doubling of coastal flooding frequency within decades due to sea-level rise. Sci. Rep. 7, 1399 (2017).
Google Scholar 

6.
Church, J. A. et al. in Climate Change 2013: The Physical Science Basis (eds Stocker, T. F. et al.) 1137–1216 (IPCC, Cambridge Univ. Press, 2013).

7.
Neumann, B., Vafeidis, A. T., Zimmermann, J. & Nicholls, R. J. Future coastal population growth and exposure to sea-level rise and coastal flooding—a global assessment. PLoS ONE 10, e0118571 (2015).
Google Scholar 

8.
Diaz, D. B. Estimating global damages from sea level rise with the Coastal Impact and Adaptation Model (CIAM). Clim. Change 137, 143–156 (2016).
Google Scholar 

9.
Barnard, P. L. et al. Dynamic flood modeling essential to assess the coastal impacts of climate change. Sci. Rep. 9, 4309 (2019).
Google Scholar 

10.
Rotzoll, K. & Fletcher, C. H. Assessment of groundwater inundation as a consequence of sea-level rise. Nat. Clim. Change 3, 477–481 (2013).
Google Scholar 

11.
Webb, M. D. & Howard, K. W. F. Modeling the transient response of saline intrusion to rising sea-levels. Ground Water 49, 560–569 (2011).
CAS  Google Scholar 

12.
Werner, A. D. & Simmons, C. T. Impact of sea-level rise on sea water intrusion in coastal aquifers. Ground Water 47, 197–204 (2009).
CAS  Google Scholar 

13.
Michael, H. A., Russoniello, C. J. & Byron, L. A. Global assessment of vulnerability to sea-level rise in topography-limited and recharge-limited coastal groundwater systems. Water Resour. Res. 49, 2228–2240 (2013).
Google Scholar 

14.
Masterson, J. P. et al. Effects of sea-level rise on barrier island groundwater system dynamics—ecohydrological implications. Ecohydrology 7, 1064–1071 (2014).
Google Scholar 

15.
Kirwan, M. L. & Gedan, K. B. Sea-level driven land conversion and the formation of ghost forests. Nat. Clim. Change 9, 450–457 (2019).
Google Scholar 

16.
Hummel, M. A., Berry, M. S. & Stacey, M. T. Sea level rise impacts on wastewater treatment systems along the U.S. coasts. Earth’s Future 6, 622–633 (2018).
Google Scholar 

17.
Liu, T., Su, X. & Prigiobbe, V. Groundwater–sewer interaction in urban coastal areas. Water 10, 1774 (2018).
CAS  Google Scholar 

18.
Knott, J. F., Daniel, J. S., Jacobs, J. M. & Kirshen, P. Adaptation planning to mitigate coastal-road pavement damage from groundwater rise caused by sea-level rise. Transp. Res. Rec. 2672, 11–22 (2018).
Google Scholar 

19.
Myers, N. Environmental refugees: a growing phenomenon of the 21st century. Phil. Trans. R. Soc. Lond. B 357, 609–613 (2002).
Google Scholar 

20.
Nicholls, R. J. et al. Sea-level rise and its possible impacts given a ‘beyond 4 °C world’ in the twenty-first century. Phil. Trans. R. Soc. A 369, 161–181 (2011).
Google Scholar 

21.
Abarca, E., Karam, H., Hemond, H. F. & Harvey, C. F. Transient groundwater dynamics in a coastal aquifer: the effects of tides, the lunar cycle and the beach profile. Water Resour. Res. 49, 2473–2488 (2013).
Google Scholar 

22.
Nielsen, P. Tidal dynamics of the water table in beaches. Water Resour. Res. 26, 2127–2134 (1990).
Google Scholar 

23.
Ketabchi, H., Mahmoodzadeh, D., Ataie-Ashtiani, B. & Simmons, C. T. Sea-level rise impacts on seawater intrusion in coastal aquifers: review and integration. J. Hydrol. 535, 235–255 (2016).
Google Scholar 

24.
Masterson, J. P. & Garabedian, S. P. Effects of sea-level rise on ground water flow in a coastal aquifer system. Ground Water 45, 209–217 (2007).
CAS  Google Scholar 

25.
Werner, A. D. et al. Vulnerability indicators of sea water intrusion. Ground Water 50, 48–58 (2012).
CAS  Google Scholar 

26.
Burnett, W. C., Bokuniewicz, H., Huettel, M., Moore, W. S. & Taniguchi, M. Groundwater and pore water inputs to the coastal zone. Biogeochemistry 66, 3–33 (2003).
CAS  Google Scholar 

27.
Hoover, D. J., Odigie, K. O., Swarzenski, P. W. & Barnard, P. Sea-level rise and coastal groundwater inundation and shoaling at select sites in California, USA. J. Hydrol. Reg. Stud. 11, 234–249 (2017).
Google Scholar 

28.
Plane, E., Hill, K. & May, C. A rapid assessment method to identify potential groundwater flooding hotspots as sea levels rise in coastal cities. Water 11, 2228 (2019).
CAS  Google Scholar 

29.
Lu, C., Werner, A. D. & Simmons, C. T. Threats to coastal aquifers. Nat. Clim. Change 3, 605 (2013).
Google Scholar 

30.
Harbaugh, A. W. MODFLOW-2005: The U.S. Geological Survey Modular Ground-Water Model—the Ground-Water Flow Process Techniques and Methods No. 6-A16 (US Geological Survey, 2005).

31.
Topologically Integrated Geographic Encoding and Referencing (TIGER) Database (US Census Bureau, 2016).

32.
Gleeson, T., Moosdorf, N., Hartmann, J. & van Beek, L. P. H. A glimpse beneath Earth’s surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophys. Res. Lett. 41, 3891–3898 (2014).
Google Scholar 

33.
Glover, R. E. The pattern of fresh-water flow in a coastal aquifer. J. Geophys. Res. 64, 457–459 (1959).
Google Scholar 

34.
Kopp, R. E. et al. Probabilistic 21st and 22nd century sea-level projections at a global network of tide-gauge sites. Earth’s Future 2, 383–406 (2014).
Google Scholar 

35.
Sweet, W. V. et al. Global and Regional Sea Level Rise Scenarios for the United States Technical Report NOS CO-OPS (NOAA, 2017).

36.
State of California Sea-Level Rise Guidance—2018 Update (California Ocean Protection Council, 2018).

37.
Duvall, A., Kirby, E. & Burbank, D. Tectonic and lithologic controls on bedrock channel profiles and processes in coastal California. J. Geophys. Res. 109, F03002 (2004).
Google Scholar 

38.
Vitousek, S., Barnard, P. L., Limber, P., Erikson, L. & Cole, B. A model integrating longshore and cross-shore processes for predicting long-term shoreline response to climate change. J. Geophys. Res. Earth Surf. 122, 782–806 (2017).
Google Scholar 

39.
Erikson, L. H., O’Neill, A., Barnard, P. L., Vitousek, S. & Limber, P. Climate change-driven cliff and beach evolution at decadal to centennial time scales. In Proc. Coastal Dynamics 2017 (eds Aagaard, T. et al.) 125–136 (2017).

40.
Limber, P. W., Barnard, P. L., Vitousek, S. & Erikson, L. H. A model ensemble for projecting multidecadal coastal cliff retreat during the 21st century. J. Geophys. Res. Earth Surf. 123, 1566–1589 (2018).
Google Scholar 

41.
Knott, J. F., Elshaer, M., Daniel, J. S., Jacobs, J. M. & Kirshen, P. Assessing the effects of rising groundwater from sea level rise on the service life of pavements in coastal road infrastructure. Transp. Res. Rec. 2639, 1–10 (2017).
Google Scholar 

42.
Habel, S., Fletcher, C. H., Rotzoll, K. & El-Kadi, A. I. Development of a model to simulate groundwater inundation induced by sea-level rise and high tides in Honolulu, Hawaii. Water Res. 114, 122–134 (2017).
CAS  Google Scholar 

43.
Hughes, J. D. & White, J. T. Hydrologic Conditions in Urban Miami-Dade County, Florida, and the Effect of Groundwater Pumpage and Increased Sea Level on Canal Leakage and Regional Groundwater Flow Scientific Investigations Report No. 2014–5162 (US Geological Survey, 2014).

44.
Guha, H. & Panday, S. Impact of sea level rise on groundwater salinity in a coastal community of South Florida. J. Am. Water Resour. Assoc. 48, 510–529 (2012).
Google Scholar 

45.
Sukop, M. C., Rogers, M., Guannel, G., Infanti, J. M. & Hagemann, K. High temporal resolution modeling of the impact of rain, tides, and sea level rise on water table flooding in the Arch Creek basin, Miami-Dade County Florida USA. Sci. Total Environ. 616–617, 1668–1688 (2018).
Google Scholar 

46.
Bakker, M. et al. Scripting MODFLOW model development using Python and FloPy. Groundwater 54, 733–739 (2016).
CAS  Google Scholar 

47.
Reitz, M., Sanford, W. E., Senay, G. B. & Cazenas, J. Annual Estimates of Recharge, Quick-Flow Runoff, and ET for the Contiguous US Using Empirical Regression Equations, 2000–2013 (US Geological Survey, 2017).

48.
Reitz, M., Sanford, W. E., Senay, G. B. & Cazenas, J. Annual estimates of recharge, quick-flow runoff, and evapotranspiration for the contiguous U.S. using empirical regression equations. J. Am. Water Resour. Assoc. 53, 961–983 (2017).
Google Scholar 

49.
Hanson, R. T., Martin, P. & Koczot, K. M. Simulation of Ground-Water/Surface-Water Flow in the Santa Clara-Calleguas Ground-Water Basin, Ventura County, California Water-Resources Investigations Report No. 2002-4136 (US Geological Survey, 2003).

50.
Hanson, R. T., Schmid, W., Faunt, C. C., Lear, J. & Lockwood, B. Integrated Hydrologic Model of Pajaro Valley, Santa Cruz and Monterey Counties, California Scientific Investigations Report No. 2014-5111 (US Geological Survey, 2014).

51.
Reichard, E. G. et al. Geohydrology, Geochemistry, and Ground-Water Simulation-Optimization of the Central and West Coast Basins, Los Angeles County, California Water-Resources Investigations Report No. 03-4065 (US Geological Survey, 2003).

52.
Nishikawa, T. A Simulation-Optimization Model for Water-Resources Management, Santa Barbara, California Water-Resources Investigations Report No. 97-4246 (US Geological Survey, 1998).

53.
Farrar, C. D., Metzger, L. F., Nishikawa, T., Koczot, K. M. & Reichard, E. G. Geohydrological Characterization, Water-Chemistry, and Ground-Water Flow Simulation Model of the Sonoma Valley Area, Sonoma County, California Scientific Investigations Report No. 2006-5092 (US Geological Survey, 2006).

54.
Bright, D. J., Nash, D. B. & Martin, P. Evaluation of Ground-Water Flow and Solute Transport in the Lompoc Area, Santa Barbara County, California Water-Resources Investigations Report No. 97-4056 (US Geological Survey, 1997).

55.
Knott, J. F., Jacobs, J. M., Daniel, J. S. & Kirshen, P. Modeling groundwater rise caused by sea-level rise in coastal New Hampshire. J. Coast. Res. 35, 143–157 (2019).
Google Scholar 

56.
Huscroft, J., Gleeson, T., Hartmann, J. & Börker, J. Compiling and mapping global permeability of the unconsolidated and consolidated Earth: GLobal HYdrogeology MaPS 2.0 (GLHYMPS 2.0). Geophys. Res. Lett. 45, 1897–1904 (2018).
Google Scholar 

57.
Gleeson, T. et al. Mapping permeability over the surface of the Earth. Geophys. Res. Lett. 38, L02401 (2011).
Google Scholar 

58.
Zamrsky, D., Oude Essink, G. H. P. & Bierkens, M. F. P. Estimating the thickness of unconsolidated coastal aquifers along the global coastline. Earth Syst. Sci. Data 10, 1591–1603 (2018).
Google Scholar 

59.
Tyler, D. J. & Danielson, J. J. Topobathymetric Model for the Southern Coast of California and the Channel Islands, 1930 to 2014 (US Geological Survey, 2018).

60.
Danielson, J. J. et al. Topobathymetric elevation model development using a new methodology: coastal national elevation. Database J. Coast. Res. 76, 75–89 (2016).
Google Scholar 

61.
Tyler, D. J., Danielson, J. J., Poppenga, S. K. & Gesch, D. B. Topobathymetric Model for the Central Coast of California, 1929 to 2017 (US Geological Survey, 2018).

62.
Tarboton, D. G. Terrain Analysis Using Digital Elevation Models (TauDEM) (Utah State Univ., 2005).

63.
Estimation of Vertical Uncertainties in VDatum (National Oceanic and Atmospheric Administration, 2018).

64.
National Oceanic Data Center (Levitus) World Ocean Atlas (National Oceanic and Atmospheric Administration, 1994).

65.
Schraga, T. S. & Cloern, J. E. Water quality measurements in San Francisco Bay by the U.S. Geological Survey, 1969-2015. Sci. Data 4, 170098 (2017).
CAS  Google Scholar 

66.
Post, V., Kooi, H. & Simmons, C. Using hydraulic head measurements in variable-density ground water flow analyses. Ground Water 45, 664–671 (2007).
CAS  Google Scholar 

67.
Befus, K. M. kbefus/ca_gw_slr Zenodo https://doi.org/10.5281/zenodo. 3897502 (2020).

68.
Befus, K. M., Hoover, D., Barnard, P. L. & Erikson, L. H. California Coastal Groundwater Projected Response with Sea-Level Rise (US Geological Survey, 2020); https://doi.org/10.5066/P9H5PBXP

69.
Befus, K. M., Barnard, P. L., Hoover, D. J., Finzi Hart, J. A. & Voss C. California saline groundwater wedge footprint model results. HydroShare https://doi.org/10.4211/hs.1c95059edcf041a0959e0b4a1f05478c (2020).

70.
Badon Ghyben, W. Nota in Verband Met de Voorgenomen Putboring Nabil Amsterdam. Tijdschr. K. Inst. Ing. 9, 8–22 (1888).
Google Scholar 

71.
Herzberg, A. Die wasserversorgung einiger Nordseebader. J. Gasbeleucht. Wasserversorg. 44, 815–819 (1901).
Google Scholar 

72.
Feistel, R. A Gibbs function for seawater thermodynamics for −6 to 80 °C and salinity up to 120 g kg−1. Deep Sea Res. I 55, 1639–1671 (2008).
Google Scholar 

73.
Kuan, W. K. et al. Tidal influence on seawater intrusion in unconfined coastal aquifers. Water Resour. Res. 48, W02502 (2012).
Google Scholar 

74.
Ataie-Ashtiani, B., Volker, R. E. & Lockington, D. A. Tidal effects on sea water intrusion in unconfined aquifers. J. Hydrol. 216, 17–31 (1999).
Google Scholar 

75.
Pool, M., Post, V. E. A. & Simmons, C. T. Effects of tidal fluctuations and spatial heterogeneity on mixing and spreading in spatially heterogeneous coastal aquifers. Water Resour. Res. 51, 1570–1585 (2015).
Google Scholar 

76.
Werner, A. D. et al. Seawater intrusion processes, investigation and management: recent advances and future challenges. Adv. Water Res. 51, 3–26 (2013).
Google Scholar 

77.
Yu, X. & Michael, H. A. Mechanisms, configuration typology, and vulnerability of pumping-induced seawater intrusion in heterogeneous aquifers. Adv. Water Resour. 128, 117–128 (2019).
Google Scholar 

78.
Strack, O. D. L. & Ausk, B. K. A formulation for vertically integrated groundwater flow in a stratified coastal aquifer. Water Resour. Res. 51, 6756–6775 (2015).
Google Scholar  More

Design, synthesis, characterization, and adsorption tests
We sought to design a HOIF through a new, simple, and scalable synthetic approach. In particular, the reaction of Zn(NO3)2, 2-methylimidazole (Hmim), and the copper chelator salicylaldoxime (H2salox) in water around 50 °C results in the formation of a zinc imidazole salicylaldoxime supramolecule, hereafter referred to as ZIOS. The material forms rapidly as small crystals in relatively high yield (76%) following the combination of separate mixtures of H2salox/Zn(NO3)3 and H2salox/mim. Synchrotron single-crystal X-ray diffraction (SC-XRD) data (Fig. 1a) revealed that the ZIOS structure features Zn3(Hsalox)4(mim)2 (i.e., Zn3(C6H4CHNOHO)4(CH3C3H2N2)2) trinuclear units, wherein H2salox and 2-methylimidazole are deprotonated and bonded directly to Zn2+ nodes through tetrahedral and pentahedral coordination by oxygen and nitrogen, respectively. Intermolecular hydrogen bonds between the hydroxyimino oxygen of H2salox and the pyrrolic nitrogen of 2-methylimidazole lead to the formation of a two-dimensional supramolecular network. The ZIOS structure is essentially nonporous, in contrast to MOFs and ZIFs28, with a Brunauer–Emmett–Teller (BET) surface area of ~12–14 m2 g−1 confirmed by both BET measurements at different activation conditions, e.g. 100 and 200 °C and grand canonical Monte Carlo simulations (Fig. 2c and Supplementary Fig. 2; c.f. BET surface area ~1200–1492 m2 g−1 for ZIF-8 determined in this work, see the Supporting Information). Even still, we find that the Cu2+ uptake capacity of ZIOS is 40% higher than that of a sample of ZIF-8 prepared in our hands, and the supramolecular structure exhibits unprecedented rapid adsorption kinetics (via infra). This surprising performance suggests the role of a chemisorption-based uptake mechanism, based on the coordination of Cu2+ by the aldoxime and histidine-inspired imidazole pendant groups present in the host structure. Powder X-ray diffraction characterization of ZIOS revealed that the bulk, as-synthesized material is highly crystalline (Fig. 2d). Thermogravimetric analysis data (Supplementary Fig. 6) collected under a N2 atmosphere revealed minimal ( More

1.
Barnett, T. P., Adam, J. C. & Lettenmaier, D. P. Potential impacts of a warming climate on water availability in snow-dominated regions. Nature 438, 303–309 (2005).
CAS  Article  Google Scholar 
2.
Huss, M. & Hock, R. Global-scale hydrological response to future glacier mass loss. Nat. Clim. Change 8, 135–140 (2018).
Article  Google Scholar 

3.
Immerzeel, W. W., van Beek, L. P. H. & Bierkens, M. F. P. Climate change will affect the Asian water towers. Science 328, 1382–1385 (2010).
CAS  Article  Google Scholar 

4.
Kaser, G., Grosshauser, M. & Marzeion, B. Contribution potential of glaciers to water availability in different climate regimes. Proc. Natl Acad. Sci. USA 107, 20223–20227 (2010).
CAS  Article  Google Scholar 

5.
Pritchard, H. D. Asia’s shrinking glaciers protect large populations from drought stress. Nature 569, 649–654 (2019).
CAS  Article  Google Scholar 

6.
Milner, A. M. et al. Glacier shrinkage driving global changes in downstream systems. Proc. Natl Acad. Sci. USA 114, 9770–9778 (2017).

7.
Marshall, S. J. et al. Glacier water resources on the eastern slopes of the Canadian Rocky Mountains. Can. Water Resour. J. 36, 109–134 (2011).
Article  Google Scholar 

8.
Clarke, G. K. C., Jarosch, A. H., Anslow, F. S., Radić, V. & Menounos, B. Projected deglaciation of western Canada in the twenty-first century. Nat. Geosci. 8, 372–377 (2015).
CAS  Article  Google Scholar 

9.
Ebrahimi, S. & Marshall, S. J. Parameterization of incoming longwave radiation at glacier sites in the Canadian Rocky Mountains. J. Geophys. Res. Atmos. 120, 12536–12556 (2015).
Article  Google Scholar 

10.
Fitzpatrick, N., Radić, V. & Menounos, B. Surface energy balance closure and turbulent flux parameterization on a mid-latitude mountain glacier, Purcell Mountains, Canada. Front. Earth Sci. 5, 67 (2017).
Article  Google Scholar 

11.
Gascoin, S. et al. Glacier contribution to streamflow in two headwaters of the Huasco River, Dry Andes of Chile. Cryosphere 5, 1099–1113 (2011).
Article  Google Scholar 

12.
Bliss, A., Hock, R. & Radić, V. Global response of glacier runoff to twenty-first century climate change. J. Geophys. Res. Earth Surf. 119, 717–730 (2014).
Article  Google Scholar 

13.
Comeau, L. E. L., Pietroniro, A. & Demuth, M. N. Glacier contribution to the North and South Saskatchewan Rivers. Hydrol. Process. 23, 2640–2653 (2009).
Article  Google Scholar 

14.
Jost, G., Moore, R. D., Menounos, B. & Wheate, R. Quantifying the contribution of glacier runoff to streamflow in the upper Columbia River Basin, Canada. Hydrol. Earth Syst. Sci. 16, 849–860 (2012).
Article  Google Scholar 

15.
Naz, B. S., Frans, C. D., Clarke, G. K. C., Burns, P. & Lettenmaier, D. P. Modeling the effect of glacier recession on streamflow response using a coupled glacio-hydrological model. Hydrol. Earth Syst. Sci. 18, 787–802 (2014).
Article  Google Scholar 

16.
Soruco, A. et al. Contribution of glacier runoff to water resources of La Paz city, Bolivia (16° S). Ann. Glaciol. 56, 147–154 (2015).
Article  Google Scholar 

17.
Greve, P. et al. Global assessment of water challenges under uncertainty in water scarcity projections. Nat. Sustain. 1, 486–494 (2018).
Article  Google Scholar 

18.
Flörke, M., Schneider, C. & McDonald, R. I. Water competition between cities and agriculture driven by climate change and urban growth. Nat. Sustain. 1, 51–58 (2018).
Article  Google Scholar 

19.
Hoekstra, A. Y. Water scarcity challenges to business. Nat. Clim. Change 4, 318–320 (2014).
Article  Google Scholar 

20.
Farinotti, D., Usselmann, S., Huss, M., Bauder, A. & Funk, M. Runoff evolution in the Swiss Alps: projections for selected high-alpine catchments based on ENSEMBLES scenarios. Hydrol. Process. 26, 1909–1924 (2012).
Article  Google Scholar 

21.
Hagg, W., Hoelzle, M., Wagner, S., Mayr, E. & Klose, Z. Glacier and runoff changes in the Rukhk catchment, upper Amu-Darya basin until 2050. Glob. Planet. Change 110, 62–73 (2013).
Article  Google Scholar 

22.
Schindler, D. W. & Donahue, W. F. An impending water crisis in Canada’s western prairie provinces. Proc. Natl Acad. Sci. USA 103, 7210–7216 (2006).
CAS  Article  Google Scholar 

23.
Downing, D. & Pettapiece, W. Natural Regions and Subregions of Alberta (Natural Regions Committee, 2006).

24.
Dee, D. P. et al. The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 137, 553–597 (2011).
Article  Google Scholar 

25.
Marshall, S. J. Meltwater run-off from Haig Glacier, Canadian Rocky Mountains, 2002–2013. Hydrol. Earth Syst. Sci. 18, 5181–5200 (2014).
Article  Google Scholar 

26.
Demuth, M. & Keller, R. in Peyto Glacier: One Century of Science (eds Demuth, M. et al.) 83–132 (Environment Canada, 2006).

27.
RGI Consortium Randolph Glacier Inventory (RGI)—A Dataset of Global Glacier Outlines (GLIMS, 2017); https://doi.org/10.7265/N5-RGI-60

28.
Bash, E. A. & Marshall, S. J. Estimation of glacial melt contributions to the Bow River, Alberta, Canada, using a radiation–temperature melt model. Ann. Glaciol. 55, 138–152 (2014).
Article  Google Scholar 

29.
Moore, R. D. et al. Glacier change in western North America: influences on hydrology, geomorphic hazards and water quality. Hydrol. Process. 23, 42–61 (2009).
CAS  Article  Google Scholar 

30.
Human Activity and the Environment: Freshwater in Canada. Section 2: Freshwater Supply and Demand (Statistics Canada, 2017).

31.
Kohonen, T. Self-organized formation of topologically correct feature maps. Biol. Cybern. 43, 59–69 (1982).
Article  Google Scholar 

32.
Shen, C. A transdisciplinary review of deep learning research and its relevance for water resources scientists. Water Resour. Res. 54, 8558–8593 (2018).
Article  Google Scholar 

33.
Fairfield, J. & Leymarie, P. Drainage networks from grid digital elevation models. Water Resour. Res. 27, 709–717 (1991).
Article  Google Scholar 

34.
Farr, T. G. et al. The Shuttle Radar Topography Mission. Rev. Geophys. 45, RG2004 (2007).
Article  Google Scholar 

35.
Hsieh, W. W. Machine Learning Methods in the Environmental Sciences: Neural Networks and Kernels (Cambridge Univ. Press, 2009).

36.
Unglert, K., Radić, V. & Jellinek, A. M. Principal component analysis vs. self-organizing maps combined with hierarchical clustering for pattern recognition in volcano seismic spectra. J. Volcanol. Geotherm. Res. 320, 58–74 (2016).
CAS  Article  Google Scholar 

37.
Steiger, M. et al. Explorative analysis of 2D color maps. In Proc. WSCG 2015 Conference on Computer Graphic, Visualization, and Computer Vision (eds Gavrilova, M. & Skala, V.) 151–160 (Union Angency, 2015).

38.
Vesanto, J., Himberg, J., Alhoniemi, E. & Parhankangas, J. SOM Toolbox for Matlab 5 Report A57 (Helsinki Univ. Technol., 2000).

39.
Strahler, A. N. Quantitative analysis of watershed geomorphology. EOS 38, 913–920 (1957).
Google Scholar 

40.
Jiskoot, H., Curran, C. J., Tessler, D. L. & Shenton, L. R. Changes in Clemenceau Icefield and Chaba Group glaciers, Canada, related to hypsometry, tributary detachment, length–slope and area–aspect relations. Ann. Glaciol. 50, 133–143 (2009).
Article  Google Scholar 

41.
Silverman, B. Density Estimation for Statistics and Data Analysis (Chapman & Hall, 1986).

42.
Water Survey of Canada HYDAT Data (Environment Canada, 2018); https://wateroffice.ec.gc.ca/mainmenu/historical_data_index_e.html

43.
Anderson, S. Alberta municipal water supply overview. Zenodo https://doi.org/10.5281/zenodo.3266447 (2019).

44.
Anderson, S. andersonsam/pca_som_streamflow: first release. Zenodo https://doi.org/10.5281/zenodo.3742162 (2020). More

1.
Misdan N, Lau WJ, Ismail AF. Seawater reverse osmosis (SWRO) desalination by a thin-film composite membrane—current development, challenges and future prospects. Desalination. 2012;287:228–37.
CAS  Article  Google Scholar 
2.
Li D, Wang H. Recent developments in reverse osmosis desalination membranes. J Mater Chem. 2010;20:4551–66.
CAS  Article  Google Scholar 

3.
Shintani T, Matsuyama H, Kurata N. Development of a chlorineresistant polyamide reverse osmosis membrane. Desalination. 2007;207:340–8.
CAS  Article  Google Scholar 

4.
Daer S, Kharraz J, Giwa A, Hasan SW. Recent applications of nanomaterials in water desalination: a critical review and future opportunities. Desalination. 2015;367:37–48.
CAS  Article  Google Scholar 

5.
Mayyahi AA. Important approaches to enhance reverse osmosis (RO) thin film composite (TFC) membranes performance. Membranes. 2018;8:68. https://doi.org/10.3390/membranes8030068.
Article  PubMed Central  Google Scholar 

6.
Xue SM, Ji CH, Xu ZL, Tang YJ, Li RH. Chlorine resistant TFN nanofiltration membrane incorporated with octadecylamine-grafted GO and fluorine-containing monomer. J Membr Sci. 2018;545:185–95.
CAS  Article  Google Scholar 

7.
Aljundi IH. Desalination characteristics of TFN-RO membrane incorporated with ZIF-8 nanoparticles. Desalination. 2017;420:12–20.
CAS  Article  Google Scholar 

8.
Xu R, Wang J, Kanezashi M, Yoshioka T, Tsuru T. Development of robust organosilica membranes for reverse osmosis. Langmuir. 2011;27:13996–9.
CAS  Article  Google Scholar 

9.
Yamamoto K, Ohshita J. Bridged polysilsesquioxane membranes for water desalination. Polym J. 2019;51:1103–16.
CAS  Article  Google Scholar 

10.
Xu R, Kanezashi M, Yoshioka T, Okuda T, Ohshita J, Tsuru T. Tailoring the Affinity of organosilica membranes by introducing polarizable ethenylene bridges and aqueous ozone modification. ACS Appl Mater Interfaces. 2013;5:6147–54.
CAS  Article  Google Scholar 

11.
Xu R, Ibrahim SM, Kanezashi M, Yoshioka T, Ito K, Ohshita J, et al. New Insights into the microstructure-separation properties of organosilica membranes with ethane, ethylene, and acetylene bridges. ACS Appl Mater Interfaces. 2014;6:9357–64.
CAS  Article  Google Scholar 

12.
Zhao L, Ho WSW. Novel reverse osmosis membranes incorporated with a hydrophilic additive for seawater desalination. J Membr Sci. 2014;455:44–54.
CAS  Article  Google Scholar 

13.
Waki M, Mizoshita N, Ohsuna T, Tani T, Inagaki S. Crystal-like periodic mesoporous organosilica bearing pyridine units within the framework. Chem Commun. 2010;46:8163–5.
CAS  Article  Google Scholar 

14.
Li Y, Li S, Zhang K. Influence of hydrophilic carbon dots on polyamide thin film nanocomposite reverse osmosis membranes. J Membr Sci. 2017;537:42–53.
CAS  Article  Google Scholar 

15.
Ohshita J, Muragishi H, Yamamoto K, Mizumo T, Kanezashi M, Tsuru T. Preparation and separation properties of porous norbornane-bridged silica membrane. J Sol Gel Sci Technol. 2015;73:365–70.
CAS  Article  Google Scholar 

16.
Mizumo T, Muragishi H, Yamamoto K, Ohshita J, Kanezashi M, Tsuru T. Preparation and separation properties of oxalylurea-bridged silica membranes. Appl Organomet Chem. 2015;29:433–8.
CAS  Article  Google Scholar 

17.
Yamamoto K, Koge S, Sasahara K, Mizumo T, Kaneko Y, Kanezashi M, et al. Preparation of bridged polysilsesquioxane membranes from bis [3-(triethoxysilyl) propyl] amine for water desalination. Bull Chem Soc Jpn. 2017;90:1035–40.

18.
Yamamoto K, Muragishi H, Mizumo T, Gunji T, Kanezashi M, Tsuru T, et al. Diethylenedioxane-bridged microporous organosilica membrane for gas and water separation. Sep Purif Technol. 2018;207:370–6.
CAS  Article  Google Scholar 

19.
Yamamoto K, Kanezashi M, Tsuru T, Ohshita J. Preparation of bridged polysilsesquioxane-based membranes containing 1, 2, 3-triazole moieties for water desalination. Polym J. 2017;49:401–6.
CAS  Article  Google Scholar 

20.
Liu Y, Tan J, Choi W, Hsu JH, Han DS, Han A, et al. Influence of nanoparticle inclusions on the performance of reverse osmosis membranes. Environ Sci Water Res. 2018;4:411–20.
CAS  Article  Google Scholar 

21.
Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, et al. Gaussian 09, revision A.02. Wallingford, CT: Gaussian, Inc; 2009.
Google Scholar 

22.
Guo M, Kanezashi M, Nagasawa H, Yu L, Yamamoto K, Gunji T, et al. Tailoring the microstructure and permeation properties of bridged organosilica membranes via control of the bond angles. J Membr Sci. 2019;84:56–65.
Article  Google Scholar 

23.
Asay DB, Kim SHJ. Evolution of the adsorbed water layer structure on silicon oxide at room temperature. Phys Chem B. 2005;109:16760–3.
CAS  Article  Google Scholar 

24.
Marsmann HC. Silicon-29 NMR. eMagRes, Wiley 2011. https://doi.org/10.1002/9780470034590.emrstm0505.pub2.

25.
Wang J, Kanezashi M, Yoshioka T, Tsuru T. Effect of calcination temperature on the PV dehydration performance of alcohol aqueous solutions through BTESE-derived silica membranes. J Membr Sci. 2012;415–416:810–5.
Article  Google Scholar  More

1.
Margat, J. F. & Gun, Jvd. Groundwater Around the World: A Geographic Synopsis (CRC Press, Boca Raton, 2013).
Google Scholar 
2.
Siebert, S. et al. Groundwater use for irrigation—a global inventory. Hydrol. Earth Syst. Sci. 14, 1863–1880 (2010).
ADS  Google Scholar 

3.
Scanlon, B. R. et al. Groundwater depletion and sustainability of irrigation in the US High Plains and Central Valley. Proc. Natl. Acad. Sci. USA 109, 9320–9325 (2012).
ADS  CAS  PubMed  Google Scholar 

4.
Dalin, C., Wada, Y., Kastner, T. & Puma, M. J. Groundwater depletion embedded in international food trade. Nature 543, 700–704 (2017).
ADS  CAS  PubMed  Google Scholar 

5.
Gleeson, T., Wada, Y., Bierkens, M. F. P. & van Beek, L. P. H. Water balance of global aquifers revealed by groundwater footprint. Nature 488, 197–200 (2012).
ADS  CAS  PubMed  Google Scholar 

6.
Aeschbach-Hertig, W. & Gleeson, T. Regional strategies for the accelerating global problem of groundwater depletion. Nat. Geosci. 5, 853–861 (2012).
ADS  CAS  Google Scholar 

7.
Richey, A. S. et al. Uncertainty in global groundwater storage estimates in a Total Groundwater Stress framework. Water Resour. Res. 51, 5198–5216 (2015).
ADS  PubMed  PubMed Central  Google Scholar 

8.
Richey, A. S. et al. Quantifying renewable groundwater stress with GRACE. Water Resour. Res. 51, 5217–5238 (2015).
ADS  PubMed  PubMed Central  Google Scholar 

9.
Rodell, M. et al. Emerging trends in global freshwater availability. Nature 557, 651 (2018).
ADS  CAS  PubMed  PubMed Central  Google Scholar 

10.
Long, D. et al. Global analysis of spatiotemporal variability in merged total water storage changes using multiple GRACE products and global hydrological models. Remote Sens. Environ. 192, 198–216 (2017).
ADS  Google Scholar 

11.
Ashraf, B. et al. Quantifying anthropogenic stress on groundwater resources. Sci. Rep. 7, 12910 (2017).
ADS  PubMed  PubMed Central  Google Scholar 

12.
Rodell, M., Velicogna, I. & Famiglietti, J. S. Satellite-based estimates of groundwater depletion in India. Nature 460, 999–1002 (2009).
ADS  CAS  PubMed  Google Scholar 

13.
Famiglietti, J. S. et al. Satellites measure recent rates of groundwater depletion in California’s Central Valley. Geophys. Res. Lett. 38, L03403 (2011).
ADS  Google Scholar 

14.
Longuevergne, L., Scanlon, B. R. & Wilson, C. R. GRACE hydrological estimates for small basins: evaluating processing approaches on the High Plains Aquifer, USA. Water Resour. Res. 46, 11517 (2010).
ADS  Google Scholar 

15.
Voss, K. A. et al. Groundwater depletion in the Middle East from GRACE with implications for transboundary water management in the Tigris-Euphrates-Western Iran region. Water Resour. Res. 49, 904–914 (2013).
ADS  PubMed  PubMed Central  Google Scholar 

16.
Feng, W. et al. Evaluation of groundwater depletion in North China using the Gravity Recovery and Climate Experiment (GRACE) data and ground-based measurements. Water Resour. Res. 49, 2110–2118 (2013).
ADS  Google Scholar 

17.
Huang, Z. Y. et al. Subregional-scale groundwater depletion detected by GRACE for both shallow and deep aquifers in North China Plain. Geophys. Res. Lett. 42, 1791–1799 (2015).
ADS  Google Scholar 

18.
Liu, C. M., Yu, J. J. & Kendy, E. Groundwater exploitation and its impact on the environment in the North China Plain. Water Int. 26, 265–272 (2001).
Google Scholar 

19.
Chaussard, E., Wdowinski, S., Cabral-Cano, E. & Amelung, F. Land subsidence in central Mexico detected by ALOS InSAR time-series. Remote Sens. Environ. 140, 94–106 (2014).
ADS  Google Scholar 

20.
Döll, P., Müller Schmied, H., Schuh, C., Portmann, F. T. & Eicker, A. Global-scale assessment of groundwater depletion and related groundwater abstractions: combining hydrological modeling with information from well observations and GRACE satellites. Water Resour. Res. 50, 5698–5720 (2014).
ADS  Google Scholar 

21.
Famiglietti, J. S. The global groundwater crisis. Nat. Clim. Change 4, 945–948 (2014).
ADS  Google Scholar 

22.
Wada, Y., Wisser, D. & Bierkens, M. F. P. Global modeling of withdrawal, allocation and consumptive use of surface water and groundwater resources. Earth Syst. Dynam. 5, 15–40 (2014).
ADS  Google Scholar 

23.
Chen, M. et al. Imaging land subsidence induced by groundwater extraction in Beijing (China) using satellite radar interferometry. Remote Sens. 8, 468 (2016).
ADS  Google Scholar 

24.
Gleeson, T. et al. Groundwater sustainability strategies. Nat. Geosci. 3, 378–379 (2010).
ADS  CAS  Google Scholar 

25.
Condon, L. E. & Maxwell, R. M. Simulating the sensitivity of evapotranspiration and streamflow to large-scale groundwater depletion. Sci. Adv. 5, eaav4574 (2019).
ADS  PubMed  PubMed Central  Google Scholar 

26.
Beijing Water Authority. Beijing Water Resources Bulletin 2003–2016 (Beijing Water Authority, 2016).

27.
Russo, T. A. & Lall, U. Depletion and response of deep groundwater to climate-induced pumping variability. Nat. Geosci. 10, 105–108 (2017).
ADS  CAS  Google Scholar 

28.
Cuthbert, M. O. et al. Observed controls on resilience of groundwater to climate variability in sub-Saharan Africa. Nature 572, 230–234 (2019).
ADS  CAS  PubMed  Google Scholar 

29.
Hartmann, A., Gleeson, T., Wada, Y. & Wagener, T. Enhanced groundwater recharge rates and altered recharge sensitivity to climate variability through subsurface heterogeneity. Proc. Natl. Acad. Sci. USA 114, 2842–2847 (2017).
ADS  CAS  PubMed  Google Scholar 

30.
Taylor, R. G. et al. Ground water and climate change. Nat. Clim. Change 3, 322–329 (2013).
ADS  Google Scholar 

31.
Trenberth, K. E. Changes in precipitation with climate change. Clim. Res. 47, 123–138 (2011).
Google Scholar 

32.
Su, B. D. et al. Drought losses in China might double between the 1.5 degrees C and 2.0 degrees C warming. Proc. Natl. Acad. Sci. USA 115, 10600–10605 (2018).
ADS  CAS  PubMed  Google Scholar 

33.
Kang, S. & Eltahir, E. A. B. North China Plain threatened by deadly heatwaves due to climate change and irrigation. Nat. Commun. 9, 2894 (2018).
ADS  PubMed  PubMed Central  Google Scholar 

34.
Scanlon, B. R., Reedy, R. C., Faunt, C. C., Pool, D. & Uhlman, K. Enhancing drought resilience with conjunctive use and managed aquifer recharge in California and Arizona. Environ. Res. Lett. 11, 035013 (2016).
ADS  Google Scholar 

35.
Schiermeier, Q. Purification with a pinch of salt. Nature 452, 260–261 (2008).
ADS  CAS  PubMed  Google Scholar 

36.
Muller, M. Lessons from Cape Town’s drought. Nature 559, 174–176 (2018).
ADS  CAS  PubMed  Google Scholar 

37.
Tapley, B. D., Bettadpur, S., Ries, J. C., Thompson, P. F. & Watkins, M. M. GRACE measurements of mass variability in the Earth system. Science 305, 503–505 (2004).
ADS  CAS  PubMed  Google Scholar 

38.
Tapley, B. D. et al. Contributions of GRACE to understanding climate change. Nat. Clim. Chang 9, 358–369 (2019).
ADS  Google Scholar 

39.
Xia, J., Wang, Q., Zhang, X., Wang, R. & She, D. X. Assessing the influence of climate change and inter-basin water diversion on Haihe River basin, eastern China: a coupled model approach. Hydrogeol. J. 26, 1455–1473. (2018).
ADS  Google Scholar 

40.
Li, X., Ye, S. Y., Wei, A. H., Zhou, P. P. & Wang, L. H. Modelling the response of shallow groundwater levels to combined climate and water-diversion scenarios in Beijing-Tianjin-Hebei Plain, China. Hydrogeol. J. 25, 1733–1744 (2017).
ADS  Google Scholar 

41.
Zhang, M. L., Hu, L. T., Yao, L. L. & Yin, W. J. Numerical studies on the influences of the South-to-North Water Transfer Project on groundwater level changes in the Beijing Plain, China. Hydrol. Process. 32, 1858–1873 (2018).
ADS  Google Scholar 

42.
Vorosmarty, C. J., Douglas, E. M., Green, P. A. & Revenga, C. Geospatial indicators of emerging water stress: an application to Africa. Ambio 34, 230–236 (2005).
PubMed  Google Scholar 

43.
Beijing Government. in 13th Five-Year Development Plan of Water Resources in Beijing (ed. People’s Government of Beijing Municipality) (Beijing Government, Beijing, 2016).

44.
Beijing Municipal Commission for City Planning and Land Resources Management. Beijing City Overall Planning (2016–2035) (Beijing Municipal Commission for City Planning and Land Resources Management, 2017).

45.
Döll, P. Vulnerability to the impact of climate change on renewable groundwater resources: a global-scale assessment. Environ Res Lett 4, 035006 (2009).
ADS  Google Scholar 

46.
Barnett, J., Rogers, S., Webber, M., Finlayson, B. & Wang, M. Transfer project cannot meet China’s water needs. Nature 527, 295–297 (2015).
ADS  CAS  PubMed  Google Scholar 

47.
AghaKouchak, A., Feldman, D., Hoerling, M., Huxman, T. & Lund, J. Water and climate: recognize anthropogenic drought. Nature 524, 409–411 (2015).
ADS  CAS  PubMed  Google Scholar 

48.
Chen, D. et al. The impact of water transfers from the lower Yangtze River on water security in Shanghai. Appl. Geogr. 45, 303–310 (2013).
Google Scholar 

49.
Webber, M. et al. Impact of the Three Gorges Dam, the South-North Water Transfer Project and water abstractions on the duration and intensity of salt intrusions in the Yangtze River estuary. Hydrol. Earth Syst. Sci. 19, 4411–4425 (2015).
ADS  Google Scholar 

50.
Wang, Y. G., Zhang, W. S., Zhao, Y. X., Peng, H. & Shi, Y. Y. Modelling water quality and quantity with the influence of inter-basin water diversion projects and cascade reservoirs in the Middle-lower Hanjiang River. J. Hydrol. 541, 1348–1362 (2016).
ADS  Google Scholar 

51.
Niswonger, R. G., Panday S. & Ibaraki, M. MODFLOW-NWT: A Newton Formulation for MODFLOW-2005 (US Geological Survey, 2005).

52.
Cao, G., Zheng, C., Scanlon, B. R., Liu, J. & Li, W. Use of flow modeling to assess sustainability of groundwater resources in the North China Plain. Water Resour. Res. 49, 159–175 (2013).
ADS  Google Scholar 

53.
Burek, P. et al. Development of the Community Water Model (CWatM v1.04): a high-resolution hydrological model for global and regional assessment of integrated water resources management. Geosci Model Dev Discuss 2019, 1–49. (2019).
Google Scholar 

54.
Kendy, E. et al. A soil−water−balance approach to quantify groundwater recharge from irrigated cropland in the north China plain. Hydrol. Process. 17, 2011–2031 (2003).
ADS  Google Scholar 

55.
Wei, M. et al. An efficient soil water balance model based on hybrid numerical and statistical methods. J. Hydrol. 559, 721–735 (2018).
Google Scholar 

56.
Fortin, F.-A., Rainville, F.-M. D., Gardner, M.-A., Parizeau, M. & Gagné, C. DEAP: evolutionary algorithms made easy. J. Mach. Learn. Res. 13, 2171–2175 (2012).
MathSciNet  Google Scholar 

57.
Beijing Water Authority. Beijing Water Resource Statistics Year Book 2005–2018 (Beijing Water Authority, 2018).

58.
Beijing Hydrology Bureau. Annual Beijing Water Monitoring Bulletins 2005–2019 (Beijing Hydrology Bureau, 2019).

59.
Dai, Y. et al. Development of a China dataset of soil hydraulic parameters using pedotransfer functions for land surface modeling. J. Hydrometeorol. 14, 869–887 (2013).
ADS  Google Scholar 

60.
Guo, Y. Q., Zhang, X., Yu, X. & Zou, Z. The increasing effects in energy and GHG emission caused by groundwater level declines in North China’s main food production plain. Agric. Water Manag. 203, 138–150 (2018).
Google Scholar 

61.
Gong, P. et al. Finer resolution observation and monitoring of global land cover: First mapping results with Landsat TM and ETM+ data. Int. J. Remote Sens. 34, 2607–2654 (2013).
ADS  Google Scholar 

62.
Wu, L. Digital Elevation Model (1 km Spatial Resolution) Dataset of China (National Basic Geographic Information Center. China 1km resolution digital elevation model dataset. National Tibetan Plateau Data Center, 2014).

63.
Messager, M. L., Lehner, B., Grill, G., Nedeva, I. & Schmitt, O. Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nat. Commun. 7, 13603 (2016).
ADS  CAS  PubMed  PubMed Central  Google Scholar 

64.
Gleeson, T., Moosdorf, N., Hartmann, J. & Van Beek, L. A glimpse beneath earth’s surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophys. Res. Lett. 41, 3891–3898 (2014).
ADS  Google Scholar 

65.
Andreadis, K. M., Schumann, G. J. P. & Pavelsky, T. A simple global river bankfull width and depth database. Water Resour. Res. 49, 7164–7168 (2013).
ADS  Google Scholar 

66.
Shen, Y. & Xiong, A. Validation and comparison of a new gauge‐based precipitation analysis over mainland China. Int. J. Climatol. 36, 252–265 (2016).
CAS  Google Scholar 

67.
Meng, X., Wang, H., Shi, C., Wu, Y. & Ji, X. Establishment and evaluation of the China Meteorological Assimilation Driving Datasets for the SWAT model (CMADS). Water 10, 1555 (2018).
Google Scholar 

68.
Okamoto, K.i., Ushio, T., Iguchi, T., Takahashi, N. & Iwanami, K. The global satellite mapping of precipitation (GSMaP) project. Proceedings. 2005 IEEE International Geoscience and Remote Sensing Symposium, 2005. IGARSS ’05. 2005, pp. 3414-3416 https://doi.org/10.1109/IGARSS.2005.1526575 (2004).

69.
Weedon, G. P. et al. The WFDEI meteorological forcing data set: WATCH Forcing Data methodology applied to ERA‐Interim reanalysis data. Water Resour. Res. 50, 7505–7514 (2014).
ADS  Google Scholar 

70.
Berrisford, P. et al. ERA-Interim Arch., version 2, 0 (2011).
Google Scholar 

71.
Lange, S. EartH2Observe, WFDEI and ERA-Interim Data Merged and Bias-Corrected for ISIMIP (EWEMBI), V.1.1 (GFZ Data Services, 2019).

72.
Department of Civil and Environmental Engineering at Princeton University. in Global Meteorological Forcing Dataset for Land Surface Modeling (ed. Princeton University) (Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory, Boulder, 2006).

73.
Kalnay, E. et al. The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc. 77, 437–472 (1996).
ADS  Google Scholar  More

Characterization
Evidence for the successful incorporation of alumina into the PES membranes after SIS treatment can be seen in X-ray photoelectron spectroscopy (XPS) measurements (Fig. 1b), which show an absence of alumina before treatment and a strong intensity Al 2p peak after treatment consistent with Al2O3. Alumina incorporation was further confirmed by differences in membrane weight before (71.6 ± 0.1 g) and after (85.9 ± 0.5 g) SIS treatment, which suggests the treated films are 17% alumina by weight. The alumina content as determined by thermogravimetric analysis (TGA) was slightly higher at 28% (Supplementary Fig. 1a). This discrepancy could be explained by mass loss during the reaction of TMA with the PES membrane during SIS, making the actual loading higher than the mass change would suggest. Overall, however, the chemical composition of the original PES polymer is unchanged in the resulting organic–inorganic composite membrane, as confirmed by FTIR measurements (Supplementary Fig. 1b) and in agreement with the previous work37.
The LIG pyrolysis mechanism has been studied previously for other polymers, though not specifically for PES, and is thought to involve thermal decomposition of C–O and C–N bonds due to the rapid increase in temperature, followed by growth of ring clusters forming graphitic structures38. In the case of PES (and other sulfonated polymers), the transformation results in clusters with some insertion of sulfur in the graphene skeleton in the form of C-S-S and -C=S bonds22. Here, the conversion of PES to LIG is observed as a visual color change in the membranes from white to black (Supplementary Fig. 2) and is verified using Raman spectroscopy, which shows the presence of D, G and 2D bands characteristic of graphene containing materials (Fig. 1c). Optical images of the membranes lased at increasing laser powers (Supplementary Fig. 2) show that a critical laser power is required for conversion to LIG, as has been seen previously for other polymers13. The images also show that without the SIS treatment, the membranes transition through two regimes: at low but non-zero power, the membranes soften at their glass transition and become transparent (Supplementary Fig. 2); at a critical laser power, they exhibit the expected visual color change, though to not nearly as dark a color as the SIS-treated membranes and at higher laser powers than those required to convert SIS-treated membranes. SIS-treated membranes do not soften but instead graphitize directly, giving the first indication that the incorporation of alumina allows the membranes to resist changes to their nanoscale features during graphitization. Although the critical laser power that is required for LIG formation depends on lasing parameters such as the spacing between the laser scanning lines and laser speed, this same trend was observed at all setting tested, including the minimum spacing between scanning lines.
Cross-sectional scanning electron microscopy (SEM) images (Fig. 1d–f) demonstrate the powerful impact of SIS treatment on the membrane’s physical stability during lasing. Figure 1d shows the cross-sectional structure of a treated PES microfiltration membrane before lasing, along with a higher magnification image showing the micropores (Fig. 1g). When PES membranes without SIS treatment are lased, the membrane’s structure collapses into a dense bottom layer (Fig. 1f, k) and an exfoliated top layer (Fig. 1j). The membrane’s total cross-sectional thickness decreases to 38–58 μm compared to the initial 140 μm (Fig. 1d). In contrast, after lasing, SEM images of SIS-treated membranes have a thickness ranging from 90 to 132 μm (Fig. 1e, Supplementary Figs. 3–6) depending on the plane of the cross-section and the laser power used, showing that most of the membrane thickness is retained. Both the top 40.2 ± 1.0 μm of the membrane (high contrast in Fig. 1e, h), which is taken to be the lased region, and the remaining 85 ± 1.0 μm bottom layer (Fig. 1e, i) maintain an open and porous structure very similar to the starting membrane (see SI and Supplementary Fig. 7 for more detail about the membrane’s surface and cross-section structure). This indicates that graphitization happens in place without any macroscale deformation.
Performance
SIS-treated LIG membranes maintain the same permeability within uncertainty before and after lasing (872 Lm−2 h−1 bar−1) (Fig. 2a). Without SIS treatment, membranes show a dramatic decrease (from 1124 Lm−2 h−1 bar−1 to 35 Lm−2 h−1 bar−1) in permeability, owing to pore-closure in the subsurface.
Fig. 2: Performance of conductive membranes.

a Permeability of PES membranes (with and without SIS treatment) before and after lasing. Error bars represent the range of data from repeated measurements. b Sheet resistance of lased membranes with and without SIS treatment as a function of the laser power used. c Applied potential required to maintain a reducing current of 10 mA cm−2 using SIS-PES-LIG electrodes over 14 days.

Full size image

In addition to maintaining their permeability, SIS-treated membranes also exhibit relatively low sheet resistance. Conductivity measurements of the LIG-coated membranes with and without SIS treatment (measured by the Van der Pauw method) as a function of laser power (Fig. 2b) show that the SIS-treated membranes achieve a sheet resistance of 37.7 ± 0.7 Ω □−1 or a conductivity of 660 S/m, based on the thickness of the conductive region seen in Fig. 1e. This sheet resistance is slightly higher than LIG formed from polyimide polymer, which show sheet resistance values as low as 15 Ω □−1, but is comparable to CNT membrane coatings which exhibit similar sheet resistances3,39. In contrast, the sheet resistance of lased membranes without SIS treatment (PES-LIG), which only graphitize and become conductive at higher laser powers (Supplementary Fig. 2), is an order of magnitude higher, in excess of 1000 Ω □−1. Their conductivity is also highly anisotropic: the two-point probe conductivity of these membranes is much higher along the direction of the laser path compared to perpendicular to the laser path.
Beyond high electrical conductivity, during operation, conductive membranes must be electrochemically stable enough to sustain either a capacitive voltage to electrostatically repel foulants or a faradaic current to electrochemically degrade foulants, strip scalants, and generate gas bubbles to remove accumulated contaminants. To verify their electrochemical stability, SIS-treated, lased membranes were subjected a reducing current of 10 mA/cm2 sufficient to perform water electrolysis40. Testing was performed on strips of the membrane surface that were dipped into 0.1 M NaCl electrolyte solution. To prevent contact between the electrolyte and the electrical contact wire due to wicking, only a small fraction of the lased area was dipped into the solution, far from the electrical contact point. In this configuration, membranes showed no loss in performance up to at least 14 days of continuous operation (Fig. 2c). An initial drop in the voltage required to drive 10 mA cm−2 of current density is attributed to the wicking of water further up the test strip during the first day of testing. Thus, laser-scribed SIS-treated membranes make excellent candidates for separations that require conductive membranes, or any technology where templated conductive structures are required.
Mechanism of improved stability
Given these favorable results, we explored the mechanisms behind the improved stability of SIS-treated membranes during lasing (i.e. suppression of deformation) and the improved conductivity of the LIG formed. One possible mechanism behind the structural stability during lasing is a change in the glass transition temperature (Tg) of the membrane due to the addition of alumina. However, DSC analysis (Fig. 3a) shows only a minor increase in the glass transition temperature (Tg) from 212 °C to 230 °C after SIS treatment. These values are consistent with previously reported Tg values for PES membranes41,42 and further indicate that the inclusion of alumina does not prevent the polymer in the membranes from undergoing a glass transition. SEM images of the membranes after the DSC measurements (i.e. after heating them above their glass transitions) (Fig. 3b, c) show that the PES completely loses its original porosity, while the alumina infiltrated sample remains mostly porous, with slight deformation. This suggests that, despite the similarity in Tg values, the presence of alumina alters the rheological properties of the SIS-PES membranes and stabilizes the membrane structure under elevated temperatures above the Tg of PES (the laser irradiation increases the PES temperature well above 230 °C)13. To test this hypothesis, we performed dynamic mechanical analysis (DMA) to measure the tensile storage (G′) and loss (G′′) modulus of PES and SIS-PES as a function of temperature (Fig. 3d). While the PES sample undergoes a full transition from a glassy state at 200 °C to terminal flow (i.e. liquid-like) behavior at 275 °C with only a mild entanglement plateau, SIS-PES shows a slight relaxation above 230 °C, but then exhibits a prolonged plateau in G′ (around 0.1 GPa) up to the instrumental limit of 400 °C. The solid-like properties of the SIS-PES sample are also reflected by the minimal sample elongation during testing, especially compared to PES, which yields over 100% and prematurely ends testing (Supplementary Fig. 8). The dramatic difference in mechanical properties between the two samples and the extended plateau in G′ of SIS-PES suggests that the alumina has formed a continuous network that stabilizes the original membrane structure well beyond the Tg of PES43.
Fig. 3: Mechanism of structural resilience during laser pyrolysis.

a DSC scans of PES membranes with and without SIS treatment showing the similar glass transition temperatures of the polymer in the two membranes. SEM image of the b PES (scale bar 2 µm) and c SIS-treated PES membranes (scale bar 2 µm) after DSC measurement showing the different pore structure between the two membranes after heating. d Storage and loss modulus measurements of PES and SIS-treated PES membranes. e Cross-section SEM image of lased, SIS-treated membrane (scale bar 4 µm) and f its EDX line-scan along the yellow arrow shown in the SEM image. Lower magnification top surface SEM images of g lased PES (scale bar 100 µm) and h SIS-treated PES membranes (scale bar 100 µm). i D to G band intensity ratios and full width at half max of the D band of Raman spectra of SIS-treated PES membranes lased at increasing laser powers. Error bars represent the standard deviation of repeated measurements.

Full size image

While the infiltrated alumina is responsible for the mechanical structural resilience of the membrane during lasing, it is unlikely to be responsible for the improved conductivity and electrochemical stability of the LIG formed. Cross-sectional SEM and energy-dispersive X-ray spectroscopy (EDX) cross-sections of SIS-PES (Fig. 3e, f) show that after lasing, the top of the film is absent of aluminum introduced by SIS, likely due to sublimation under the high temperatures induced by the laser. At the LIG/PES interface, there is a small region where the alumina appears to have ripened into nanoparticles coating the polymer/LIG film but the conductive region of the resulting film is completely absent of the infiltrated alumina in its original crosslinked structure. Note that Fig. 3e is an SEM image of a membrane lased at lower laser power (14%) than the SEM image shown in Fig. 1e, which is why the thickness of the conductive layer is different. Further analysis of the PES-LIG and SIS-PES-LIG using XPS does not reveal significant chemical compositional differences. Both materials show primary C 1s peaks at a binding energy of 284.4 eV in XPS fine scans, consistent with sp2-bonded carbon13, with some additional higher binding energy peaks (Supplementary Fig. 9). While the PES-LIG shows a greater intensity of higher binding energy carbon peaks (286.2 eV), suggesting more ether carbon remains in the films after laser treatment, it is unlikely this difference in composition would lead to such a drastic difference in sheet resistance. Al 2p fine scans of the SIS-PES-LIG (Supplementary Fig. 10) also show that the aluminum within the film has remained in an oxide form, eliminating the possibility for the formation of Al metal. Raman spectra of PES-LIG and SIS-PES-LIG also do not exhibit significant enough differences to account for differences in conductivity. The D to G peak ratios in Raman spectra of carbonaceous materials are often used to make qualitative statements about the nature of LIG such as the crystallite size of graphitic clusters44. The peak ratios of the Raman spectra for PES-LIG and SIS-PES-LIG (Fig. 1c) suggest similar crystallite sizes ranging from 7–10 nm. However, LIG from untreated PES have broader peaks indicating increased disorder in these films relative to the SIS-treated samples45. In addition, some spots on the untreated PES films showed significant fluorescence (Supplementary Fig. 11a), indicating the presence of regions with minimal LIG coverage. Surface SEM images of the lased PES membrane without SIS show a heavily exfoliated structure with order-of-magnitude larger features (Supplementary Fig. 11a), and continuous regions of non-porous polymer underneath and in between, which is consistent with these Raman spectra. In contrast, spectra from SIS-treated PES after lasing showed consistent bands (Supplementary Fig. 11b) and SEM images of these membranes show a more homogeneous structure, similar to the pore-structure of the starting membrane.
The drastic differences in conductivity are therefore likely due to the less homogenous coverage of the lased surface by LIG. While membranes with and without SIS treatment show similar LIG thickness after lasing (Fig. 1e, f), lower magnification top-down SEM images of PES-LIG (Fig. 3g) reveal that after laser scribing, the laser creates linear trench structures of graphitized regions, separated by large gaps of non-porous polymer. These gaps between LIG structures also explain the high degree of anisotropy in conductivity mentioned earlier. A structural anisotropy is also visible in SIS-treated membranes (Fig. 3h). Previous work has shown that differences in lasing conditions can lead to vastly different LIG morphologies, in part due to differences in the resulting polymer temperature and anneal rate14. Here, regions outside of the direct laser path are still graphitized without any loss of porosity, but they likely do not reach the same annealing temperature and therefore exhibit different morphology.
Ultimately, the temperature that is achieved during laser irradiation plays a critical role in the resulting pyrolysis process and is affected by a combination of factors, such as the total absorption of IR light by the polymer and alumina (when present), the heat capacity of the membranes with and without alumina (which can be qualitatively deduced from the DSC data (Fig. 3a)), the structural changes due to polymer softening, the loss of porosity (without alumina), the rate of heat dissipation, and the total mass of polymer irradiated. While the relative impact of each of these factors was outside the scope of this work, it can be reasonably concluded that the presence of alumina increases the resulting membrane temperature during laser irradiation, since there is a reduction in critical laser power required for graphitization (Supplementary Fig. 2): the SIS-PES membranes can reach a higher temperature with lower laser energy and thus require less laser power to achieve the temperature required for graphitization.
The saturation in sheet resistance of the SIS-treated membranes above a laser power of 20–22% is also partially explained by the anisotropy in lased area. Initially, increasing the laser power is correlated with improved conductivity (Fig. 2b), as a greater fraction of the top surface of the membrane is converted to LIG and the laser continues to penetrate deeper into the membrane bulk. At 14% power, the laser is able to convert only parts of the membrane to LIG, which is seen as dark regions with low alumina content in surface SEM images, EDS maps, and EDS line scans (Supplementary Fig. 12). Cross-sectional images at 14% power show that the dark regions are part of a hemispherical path scribed by the laser (Supplementary Fig. 12). As the laser power is increased, the radius of this path grows and neighboring paths overlap until full surface coverage is achieved at which point the conductivity also saturates. We analyzed Raman spectra of the lased surfaces at different powers to see if higher laser powers changed the graphitic nature of the formed LIG. Between 14 and 18% laser power, there is an increase in the D to G peak intensity and area ratio, a narrowing of the D band full width at half max (Fig. 3i), and an emergence of a 2D peak (Supplementary Fig. 13). These changes plateau at higher powers, mirroring the conductivity measurements. These trends seen in the Raman spectra indicate that although initially, higher laser power leads to less disorder in the type of defects found on the LIG and increased stacking of graphitic clusters, defects found on the graphitic regions persist even at high powers46. Thus, continued improvement in the sheet resistance of the membranes will require alternative approaches, such as other chemical treatments besides alumina.
Interfaces and interfacial properties play a central role in many technologies other than membranes at the water and energy nexus2. PES membranes are taken as a prototypical example of porous polymers where nanosized features need to be preserved during laser scribing. The method of combining SIS of organometallic precursors into polymers with LIG formation described in this paper can be generalized to other applications where maintaining micro and nano-sized features of polymers at temperatures well above their Tg is desired. An even broader advantage of this approach is the improved mechanical and chemical properties that are observed at these high temperatures, without the need to change the chemistry of the underlying polymer28. Other modification approaches that could potentially be used, such as crosslinking of the polymer, would alter the polymer chemistry, potentially altering the material’s LIG forming properties. Since stability of LIG formed from polymers is a known concern47, the method described herein could be beneficial for all other polymers used for LIG formation as well.
In summary, we present a simple, solvent-free process for making conductive membrane coatings without altering underlying polymer structure. We demonstrate how infiltration with alumina stabilizes the PES membrane against deformation above the glass transition temperature37, allowing it to maintain its structure during laser treatment. These membranes are shown to be more conductive than LIG formed directly from the bare polymer, are electrochemically stable and maintain their permeability after lasing. These results demonstrate the immense versatility of hybrid polymer-ceramic materials as a promising class of materials to be used in conjunction with the LIG process. More

The HYSETS database contains a multitude of datasets that, when combined, provide a rich and comprehensive environment for testing in various applications such as in hydrological modelling. Four categories of data are combined to provide this sandbox environment: Hydrometric, watershed delineation, meteorological and physiographic data. The methods to extract, validate and combine these sources of data are presented according to each category.
Hydrometric data
The daily hydrometric data were collected independently for the three countries covered in this database, namely Canada, the United States and Mexico. The Canadian hydrometric data were provided through the Environment and Climate Change Canada (ECCC) Water Survey Canada (WSC) National Water Data Archive (HYDAT), available at: https://collaboration.cmc.ec.gc.ca/cmc/hydrometrics/www/. The data are downloadable in a single Microsoft Access table and include station metadata, such as station location, drainage area and flow regime, as well as the actual daily flow data. A filter was applied to select only the stations whose flow regime is natural, i.e. is not regulated by man-made structures over the study period. This was performed to ensure that the hydrometric data and meteorological data were as naturally correlated as possible.
Hydrometric data for the United States were collected from the United States Geological Survey (USGS) National Water Information Service (NWIS) web portal23, available at: https://waterdata.usgs.gov/nwis/uv/?referred_module=sw. Data were batch-downloaded and processed to obtain full time-series for each station. While some metadata is made available such as drainage area, station coordinates and statistics on historical flow, there is no information on the flow regime or on the presence of regulation structures. An alternative method to filter regulated structures was thus devised. To determine which sites were regulated or were affected by regulation, the streamflow dataset was cross-checked against the peak flow statistics database from NWIS. Stations whose peaks were influenced by regulation works at least once in the archive were removed from the HYSETS dataset. Therefore, all regulated, or partly regulated stations were excluded to keep the data as close to natural as possible. While this method seems to have produced the desired results, it is possible that the filtering method is not perfectly accurate and as such, users should verify individual hydrographs if any doubts arise on the river’s regulation status.
Finally, hydrometric data in Mexico were collected from the “Banco Nacional de Datos de Aguas Superficiales” (BANDAS) produced and maintained by the National Water Commission (CONAGUA) from Mexican Ministry of Environment. The data were downloaded from the BANDAS web portal and was filtered according to data quality (visual inspection, hydrological model calibration performance) and time series length. The BANDAS data are available at: www.conagua.gob.mx/CONAGUA07/Contenido/Documentos/Portada%20BANDAS.htm.
Metadata is available for every station in the database and includes station coordinates, drainage area and information about the presence/absence of regulation structures. The stations included in the HYSETS database are all located on basins exempt from any regulation structures. Finally, a filter was applied to ensure that all stations had at least one year of recorded streamflow data. All hydrometric data were converted into units of cubic meters per second (m3/s).
Watershed delineation
The watershed delineation boundaries are a critical component of the dataset as all meteorological data need to be extracted according to those limits for each watershed. For most of the watersheds, the water management agencies provide the official boundaries directly in the form of a shapefile or geodatabase. The Canadian data are available at: http://donnees.ec.gc.ca/data/water/products/national-hydrometric-network-basin-polygons/?lang=en and the United States boundaries are available at the following website: water.usgs.gov/GIS/metadata/usgswrd/XML/streamgagebasins.xml.
The drainage area was made available for most hydrometric gauges by the water management agencies that collate and curate those sources of data. However, a filter was applied to remove all stations that did not have an official drainage area value at the hydrometric gauge, as the value is key in determining if the watershed bounds are acceptable or not. The drainage areas were validated using the watershed delineation boundaries as described above in the geospatial analysis software QGIS 3.4. However, in some instances, the water management agencies did not provide watershed boundary files as they had not been produced or made available publicly. In those cases, estimated watershed contours were taken from the Global Streamflow Indices and Metadata (GSIM) project21,22 where available. For catchments where GSIM boundaries were kept for the data extraction a flag (“flag_GSIM_boundaries”) was set to 1 to inform users that the boundaries are from GSIM and not from the official agencies. The GSIM-derived area is also identified in those cases in the dataset, under the “Drainage_Area_GSIM_km2” heading. For catchments smaller or equal to 50 km2 in size according to the official gauge, a bounding box equal to the surface area around the catchment outlet was provided as the contour of the catchment as at those scales catchment delineations are difficult due to the resolution and hydrometric gauge accuracy. These catchments are represented by the “flag_artificial_boundaries” indicator in the dataset files. Furthermore, weather data and other catchment attributes are coarser than the area of the catchments in most cases.
In the HYSETS database, all boundaries are provided in a WGS84-projected ESRI shapefile and include the following properties: Watershed ID (to link to the data in the netCDF files), source of the data, name of the hydrometric station, official ID of the hydrometric station and the flag to identify if the boundaries were derived from GSIM. All drainage areas are in km2.
Meteorological data
The HYSETS database contains meteorological data from five sources. The following sections detail the methods applied to integrate the data into the database at the catchment scale. All precipitation data are provided in millimetres per day (mm/d), and temperature values are in degrees Celsius (°C).
Station data
Three weather station products were used to cover the North-America domain: The Environment and Climate Change Canada (ECCC) weather stations for Canada, available at: https://climate.weather.gc.ca/, the Global Historical Climate Network Daily (GHCND) station database24,25 for the United States and Mexico, available at: https://www.ncdc.noaa.gov/ghcnd-data-access and the station-based serially complete dataset for North America (SCDNA)26, available at: https://zenodo.org/record/3735534. The ECCC and GHCND datasets were combined to provide a North American weather station dataset for raw observations with a potential incomplete coverage for the desired 1950–2018 period. The SCDNA was also added to provide another dataset of stations with complete records between the 1979–2018 period.
The daily historical ECCC weather data were downloaded from the ECCC web portal for the years 1950–2018. Data include daily precipitation, maximum and minimum temperature for over 8578 stations across Canada, with varying levels of data completeness and record length. While the GHCND database contained Canada’s weather stations as well, the ECCC database was more complete and was preferred. Over Mexico, a total of 5249 precipitation and 5071 temperature stations are available in the GHCND dataset. Similarly, there are 55693 precipitation and 16011 temperature stations across the United States. While these stations undergo some quality control (QC), both data with and without QC were extracted and provided in the HYSETS database.
These stations were combined into a 69520-station dataset for precipitation and 29960-station dataset for temperature. The extraction process was performed twice for station data: Once using the quality controlled GHCND database and another using all available data, even the non-quality-controlled data.
The SCDNA provides daily precipitation, maximum and minimum temperature for 27280 stations over North America. A strict quality-control was performed on the original station data by the SCDNA authors. Missing data were infilled/reconstructed using information from neighbouring stations as well as three reanalyses products (ERA5, JRA-55 and MERRA-2). Strategies based on quantile mapping, statistical interpolation and machine learning were used to implement the corrections. Overall, this dataset was shown to provide a better agreement to station observations compared to four gridded products.
All three station datasets (ECCC/GHCND with and without QC and SCDNA) were weighted separately using Thiessen polygons for each watershed. The stations contributing to the Thiessen polygon calculation were those located within an artificial boundary defined as the real watershed boundaries extended by a buffer of 1° of latitude and longitude. This step was performed in order to exclude stations that would be too distant to represent the watershed conditions. Due to the highly variable nature of gauge-data quality and station longevity in the case of the ECCC/GHCND combined datasets, Thiessen polygons were computed for each day, using the available data for each day. Therefore, there are some discontinuities when a station is added or removed, or when a station temporarily has no record for a given period of time. There are also cases where no data at all were available for a given period. In those cases, the meteorological data fields are set to NaN, and the user is encouraged to replace those data using other means (either manual replacement or replacement with one of the gridded products as described below). The 1950–1978 period for the SCDNA dataset was also set to NaN.
All three sets of catchment-averaged data are available in the HYSETS database, opening the possibility of performing quantitative assessments of the impacts of using more (but perhaps less reliable) meteorological data in impact studies.
Natural resources canada gridded climate data for canada
The Natural Resources Canada (NRCan) gridded climate data product was made available from Natural Resources Canada’s Canadian Forest Service and covers the entirety of Canada up to approximately 84°N latitude27,28,29. It includes daily precipitation, maximum and minimum temperature data on a daily scale on a ~10 km spatial grid. Data cover the period 1950–2010 inclusively. Data points falling within the catchment boundaries were averaged to obtain a single time-series of continuous data as the NRCan dataset contains no missing data. When catchments were too small to contain a data point, the closest data point to the catchment centroid was used to populate the time series for that watershed. Information on the NRCan dataset can be found at the following website: https://cfs.nrcan.gc.ca/projects/3/4.
Livneh gridded climate data for continental USA, Mexico and southern Canada
The Livneh database includes interpolated precipitation and temperature data on a regular 0.0625 × 0.0625° grid over the continental United States, Mexico and southern Canada30,31. The data cover the period 1915–2015, although only the portion 1950–2015 was used in this dataset. It includes daily precipitation, maximum and minimum temperatures for the entire period without any missing data. Some discontinuities are present at the United States/Mexico border as station density and quality differ and influence the interpolation process. The same is also present but at a smaller scale on the United States/Canada border. The catchment-averaging process was the same as for the NRCan dataset. The Livneh data were provided by the NOAA/OAR/ESRL PSD, Boulder, Colorado, USA, from their Web site at https://www.esrl.noaa.gov/psd/.
ERA5 and ERA5-Land reanalyses
The ERA532 and ERA5-Land33 reanalyses are hourly products developed by the European Center for Medium-Range Weather Forecasting (ECMWF). The reanalyses provide estimates of a multitude of hydrometeorological and atmospheric variables including precipitation and temperature on regular grids covering the entire surface of the Earth. ERA5 was first implemented with a 0.25° x 0.25° spatial resolution and data are available from 1979–2019, although the HYSETS database stops in 2018 to remain consistent with the other data sources. The ERA5-Land reanalysis is a refined version of ERA5 with a spatial resolution of approximately 9 km. It was driven by the ERA5 reanalysis and a mask was applied such that only land masses are modelled in the refined domain. ERA5-Land covers the period 1981-onwards, and as such the years 1981–2018 are available at present in the HYSETS database.
As both the ERA5 and ERA5-Land products are hourly and data are archived without UTC offsets, it was required to shift the data according to the grid point locations. The longitude of the watershed centroid was used to assess to which time zone it belongs. Based on this time zone, the hourly data were shifted from the same number of hours to realign the daily cycle between 00:00 and 24:00. For instance, a station located in the time zone -7 will have the whole time series shifted by 7 hours to match its proper daily cycle. The ERA5 and ERA5-Land data were downloaded from the Copernicus Climate Data Store, available at: https://climate.copernicus.eu/climate-reanalysis. Note that the HYSETS dataset contains modified ERA5 and ERA5-Land reanalysis data from the Copernicus Climate Change Service Information and that neither the European Commission nor ECMWF is responsible for any use that may be made of the Copernicus Information or Data it contains.
Once the data processing was complete to bring it to the daily scale, the extraction followed the same process as for the NRCan and Livneh datasets, however, for the reanalysis products, the data are available for all watersheds given that reanalyses cover the entire globe.
Snow Water Equivalent (SWE) data
Two snow-water equivalent databases are provided in HYSETS. The first is a 9-year (2010–2018) time series of watershed-averaged daily high-resolution (roughly 1 km) data provided by the Snow Data Assimilation System (SNODAS) analysis34 available at: https://nsidc.org/data/g02158. The SNODAS data were averaged at the catchment scale using points within the watershed boundaries, or the closest point for the smallest of watersheds that did not contain a point within their limits. SWE data units are millimetres (mm) and represent the value expected on the ground for that day. Any missing values are replaced by NaNs. SNODAS data incorporate multiple sources of data and provide the best possible estimate of SWE at each point location. HYSETS simply averages those values at the catchment scale for users to easily interpret and analyse the data with respect to the rest of the data. SNODAS’ main limitation is that its spatial coverage includes the continental USA as well as the lower portions of Canada below 54°N latitude. This means that many snowy catchments in northern Canada and Alaska are not covered by SNODAS.
The second dataset is the ERA5-Land reanalysis product, which covers the period 1981–2018. It was extracted using the same method as for the ERA5-Land precipitation and temperature data and was also averaged at the daily scale using an hourly UTC offset depending on the time zone. One main advantage of the ERA5-Land reanalysis SWE is that it covers the entire globe and thus is available for all catchments, even those above 54°N that are not covered by SNODAS.
Physiographic data
One of the strengths of the dataset is the inclusion of a multitude of properties to describe and characterize each of the watersheds. The process is similar for all of the properties, but some variables required slightly more complex operations than others.
The first set of data was based on geographic and topographic properties and was derived from the EarthEnv-DEM90 digital elevation model35 available at: https://www.earthenv.org/DEM.html. The process was performed in a meta-software called PAVICS-Hydro (Power Analytics and Visualization for Climate Science – hydrological modelling toolbox) being developed for this purpose, available at: https://pavics-sdi.readthedocs.io/. This set includes mean watershed elevation (meters), slope (degrees), aspect (degrees), Gravelius index (unitless) and perimeter (kilometers). The elevation and perimeter are self-explanatory. The slope is the average slope when considering the individual elevation differences between tiles and can be seen as an indicator of the catchment relief, with higher slopes indicating more mountainous regions. The aspect is the main orientation of the catchment, i.e. where the average slope points towards. The Gravelius index is the ratio of the perimeter of the watershed compared to the perimeter of a circle of the same area. Higher values indicate more elongated or less compact catchments. All the DEM points falling within a watershed boundary were used to compute these characteristics.
The database also provides elevation band data for each catchment in 100-meter intervals. This information was extracted from the EarthEnv 90 m DEM by applying a zonal histogram to a reclassified DEM based on 100-meter intervals in the QGIS software. The data are provided in a separate.csv file “HYSETS_elevation_bands.csv” and represent the percentage of catchment area lying below that elevation. Therefore, the curves are cumulative sums of these areas. This will allow users to provide information to more complex routines and models such as the well-known CemaNeige snow accounting and melt model. The elevation bands can also be used to adjust precipitation rates and temperatures per precipitation bands based on the average catchment elevation and the elevations in the bands, depending on the desired lapse rates and correction methods.
The land use percentages reflect which fraction of the watershed is covered in the different classification categories. The North American Land Change Monitoring System (NALCSM) imagery data from 2010 was used for this purpose36,37. NALCMS was developed by Canada, the United States and Mexico to track the evolution of land use over time. A static dataset for 2010 is available and contains 19 land use classes values that were combined to form 7 meta-categories: forests, shrubs, croplands, wetlands, water, urban and permanent snow/ice. For example, coniferous, deciduous and mixed forests were combined into the “forest” category. The original data are available at the Commission for Environmental Cooperation website:

For each watershed, the NALCMS raster dataset was queried through the PAVICS zonal statistics toolbox for each of the 19 original categories. The values were then aggregated to the 7 categories used in the HYSETS dataset and the relative fraction of each was computed.
The dataset used to characterize catchment geology is the GLobal HYdrogeology MaPS (GLHYMPS) of subsurface permeability and porosity38. This dataset provides quantitative estimates of permeability and porosity below the soil horizon. Catchment averages of these two variables were calculated, by considering the contribution of each spatial polygon being weighted by the fraction of a catchment it covers. The arithmetic mean was used for porosity, but for permeability, the geometric mean was taken. The same process as for NALCMS was performed. However, a pre-processing of the GLHYMPS vector data was performed to transform it to a raster format. The same zonal statistics tools from PAVICS were used to extract the average values of both variables for each catchment. The permeability units are in m2 whereas the porosity is archived as a fraction. The GLHYMPS dataset is available at: https://dataverse.scholarsportal.info/dataset.xhtml?persistentId=doi:10.5683/SP2/DLGXYO.
For all the above-mentioned properties, if for any reason the extraction process could not be performed (watershed out of product boundaries, unavailability at a given location) the values are replaced by NaNs and a flag was set in the metadata and properties file.
Monthly hydrometeorological data
One aspect that must be noted is that the data are all averaged temporally at the daily scale and spatially at the catchment scale. This means that for large catchments, it is possible that the daily data are not very representative of the localized precipitation and runoff events. For this reason, the HYSETS database also includes monthly-aggregated data for all catchments, which will allow evaluating products in terms of bias and mass balance over longer periods. The daily data are still made available for all catchments and the users are invited to consider which timescale is more appropriate for their use-case.
The monthly hydrometeorological data include data from the seven temperature and precipitation data sources as well as for the streamflow. SWE values were not provided at the monthly scale as they are typically not as variable as other variables. All monthly data are combined into a single netCDF file named “HYSETS_2020_monthly_meteorolgical_data.nc”. More