Screening the Future Security of Groundwater-Based Water Supply on the Scale of a Regional Water Supplier

This repository contains all data required to reproduce the calculations from the corresponding publication. A brief description of the content is given below.

 

#=== Population ===#

"ASG_ew_monthly.csv"

Historical Population per municipality, monthly resolution, actual number of people supplied by ASG (=/= people living in municipality)

 

"scaled_population_projections.csv"

Population projections, monthly resolution, shows relative change to initial state (2021) until 2100 (after 2070 extrapolation) p(t) = p(0)*(1+x(t))

Based on population projections by the federal statistical office, Germany 

(https://www.destatis.de/DE/Themen/Gesellschaft-Umwelt/Bevoelkerung/Bevoelkerungsvorausberechnung/Publikationen/Downloads-Vorausberechnung/statistischer-bericht-bvb-bundeslaender-2070-5124205219005.html?nn=208696)

 

"EW_HBG.csv"

Population of the city of Herrenberg including the supplied districts, annual resolution

 

#=== Water Demand ===#

"Monthly_processed.csv"

Historical water delivery by ASG per municipality, monthly resolution [m^3/month]

 

"Wasserverbrauch ab 2000.csv"

Historical water delivery by the city of Herrenberg (only parts of the city are supplied by the city), annual resolution [m^3/year]

 

"bwv.csv"

Historical water delivery by the distance water supplier bwv, annual resolution [m^3/year]

 

"total_extraction_asg.csv"

Historical total water extraction between 2015 and 2021 by the ASG [m^3/month], used for the GWUI correlation analysis

 

#=== Spatial Data ===#

"LURCH_GrundFl_v1.zip"

Contains "LURCH_GrundFl_v1.shp"

Spatial information (shp) of the polygons for groundwater recharge calculation (EPSG:31467 - DHDN / 3-degree Gauss-Kruger zone 3)

Shapefile contains 6112 individual polygons corresponding to 6050 patches, column "GRUND_ID" contains identifier for each polygon

Alternatively Polygons can be identified by their x- and y-coordinate.

 

"LURCH_GrundFl_v1.dat"

Parameters/attributes of the polygons 

 

"ASG Recharge Area.zip"

Contains "ASG Recharge Area.shp", Polygons delineating the catchments of the two well fields (EPSG:31467 - DHDN / 3-degree Gauss-Kruger zone 3)

 

#=== Weather and Climate Data ===#

"Lurch01_(...).dat"

Historical weather and seepage data, "Summenwert_MonatEinzel": Monthly sum, "Tageswert_Monateinzel": Monthly average, "Tageswert_HalbjahrEinzel": half-year average

"Niederschlag": Precipitation [mm/month]

"Sickerwasser": Seepage rate (=/= gw recharge), must be multiplied with base flow index ("QbasQges" in "LURCH_GrundFl_v1.dat")

"Temperature": Temperature

 

"empProj2(...).dat"

Empirical climate and seepage projections

 

"Lurch-ProjXX(...).dat"

Processed climate projections, with seepage rate calculations *

 

"ProjModel.csv"

*Each climate projection (XX) is based on a global circulation model, a down scaling model, and a realisation. These are listed here.

 

All files in "Weather and Climate Data" are generated by the program "GWN-BW". Each file has 6050 rows and the column number

equals the number of time steps. Each row represents one polygon from "LURCH_GrundFl_v1.shp". Note that the time steps are listed in the header, and the first two columns contain x and y coordinate of the polygon center. When loaded in "Python" with the "pandas" package, the first column (x) will be the index column and the second column (y) will be in the "TimeStep" column. When loaded in "Excel" or "LibreOffice" columns can be shifted.

To relate the rows to the individual polygons from "LURCH_GrundFl_v1.zip", the "GRUND_ID" column from "LURCH_GrundFl_v1.dat"

can be concatenated with the weather data file as an additional column. The order of patches in "LURCH_GrundFl_v1.dat" is the same 

as in each of the weather data files. 

(patch_properties = pandas.read_csv("LURCH_GrundFl_v1.dat"), file = pandas.read_csv("empProj2(...).dat"), file["Grund_ID"] = patch_properties["Grund_ID"])
Alternatively, polygons can be identified by their x and y coordinate (EPSG:31467 - DHDN / 3-degree Gauss-Kruger zone 3).