Coastal groundwater systems: mapping chloride distribution from borehole and geophysical data

Information on chloride (Cl) distribution in aquifers is essential for planning and management of coastal zone groundwater resources as well as for simulation and validation of density-driven groundwater models. We developed a method to derive chloride concentrations from borehole information and helicopter-borne electromagnetic (HEM) data for the coastal aquifer in the Elbe-Weser region where observed chloride and electrical conductivity data reveal that the horizontal distribution of salinity is not uniform and does not correlate with the coastline. The integrated approach uses HEM resistivity data, borehole petrography information, grain size analysis of borehole samples as well as observed chloride and electrical conductivity to estimate Cl distribution. The approch is not straightforward due to the complex nature of the geology where clay and silt are present. Possible errors and uncertainties involved at different steps of the method are discussed.


Introduction
Worldwide coastal aquifers are increasingly endangered by saline water intrusion from oceans and rivers (van Weert et al. 2009). Saline water intrusion is triggered by several driving forces, natural and anthropogenic, that act on different spatial and temporal scales. Changes in rainfall patterns and sea level rise driven by climate change interact with increased groundwater exploitation, land drainage, and urbanization (White and Kaplan 2017). The effects of these driving forces make sustainable coastal zone groundwater management very complicated. Accurate groundwater resource assessment and information on the spatial distribution of salinity is required (Barlow and Reichard 2010). Groundwater resource assessment allows the planning of drinking water supplies which involves considering present and future water demand (Cosgrove and Loucks 2015). Understanding the three-dimensional distribution of chloride supports the siting of groundwater extraction wells and the definition of the proper filter depth. Groundwater models can support this planning action (Dogrul et al. 2016) but for density-dependent groundwater modelling, information on existing freshwater head and chloride distribution is required at sufficient resolution.
Direct and indirect methods are available to estimate chloride (Cl) content in the groundwater (Goes et al. 2009). Direct methods involve groundwater sampling and measurement of Cl concentration in the laboratory. This analysis is more expensive than the indirect measurement of groundwater electrical conductivity (ECw) that can be measured in the field by a conductivity meter (Peinado-Guevara et al. 2012). It is common practice to install monitoring stations for groundwater levels and salinity measurements, and in some cases, data availability is adequate and fulfils the requirement for groundwater modelling. However, as the chloride concentration varies with depth, multilevel moni-toring stations would be necessary (Margane 2004), but are not often realized due to cost issues.
Some indirect survey techniques can support the required data density for mapping the chloride concentration in a coastal aquifer (Table 1). These indirect techniques are mostly based on ECw or the bulk electrical resistivity (Rb). Both can be translated into chloride concentration by empirical equations that require aquifer and water properties, such as the formation factor or porosity.
A number of studies apply geophysical surveys for salinity estimation or fresh-saline groundwater interface delineation. In particular, airborne electromagnetic (AEM) methods enable spatial mapping of Rb with high resolution (in the order of tens of meters to hundreds of meters) along the line and acceptable resolution (hundreds of meters) across the lines ). Vertical resolution decreases with depth (from meters to tens of meters) and depends on system parameters and the subsurface conductivity distribution, which limits the depth of investigation (Christiansen and Auken 2012). Large-scale groundwater surveys increasingly apply AEM methods to investigate large areas in reasonable time and at relatively low costs (Paine and Minty 2005;Siemon et al. 2009). AEM results have been used not only for mapping purposes, but also as base-line data estimation for geological and groundwater modelling (e.g. De Louw et al. 2011;Gunnink et al. 2012;Sulzbacher et al. 2012).
Ground-based resistivity methods like Vertical Electric Sounding (VES) for fresh-saline water interface mapping were used by Kumar et al. (2015 (India)), supported by water quality parameter estimation, and Goes et al. (2009 (the Netherlands)), supported by ECPT, borehole logs, ECw and chloride measurements. Paine (2003 (Texas, USA)) used AEM and borehole induction logs to map the resistivity distribution in the groundwater.
The spatial distribution of chloride in groundwater was mapped by Delsman et al. (2018 (Zeeland, the Netherlands)) using frequency-domain helicopter-borne electromagnetic (HEM) and Peinado-Guevara et al. (2012 (Mexico)) using VES; both studies were supported by chloride measurements.
In this paper an integrated approach is presented to estimate the formation factor distribution from HEM resistivity data, borehole petrography information, and grain size analysis of borehole samples. Cl concentrations are then calculated from observed chloride and electrical conductivity (EC) values. A final map is produced showing the Cl distribution in the aquifer at several depth intervals. The method is applied to a coastal area of Lower Saxony (Ger-many) between the Weser and Elbe rivers. The obtained Cl distribution is used to model its future development under climate change conditions (González et al. 2020). The method is general and applicable to other regions with comparable geological conditions. We critically discuss possible errors and uncertainties involved at different steps of the method.

Study area
The study area is situated in the coastal area of Lower Saxony, Germany, between the Elbe and Weser rivers (Fig. 1). This study is performed as part of the project TOPSOIL in the framework of the EU INTERREG VB North Sea Region Programme. The analysis covers an area of about 1700 km 2 , and the depth range covered by the HEM survey is between the surface and 30-160 m depending on the bulk electrical conductivity distribution in the subsurface. The study area is bounded by the North Sea in the west (ca. 37 km coastline) and the Elbe River in the north and east (80 km long reach). Another major river within the area is the Oste (reach ca. 50 km). The topography of the area is relatively flat (marshlands) with some highlands (Geest)

Geology and hydrogeology of the area
The geology of Elbe-Weser region is relatively complex and characterized by sedimentation in Quaternary glacial and interglacial periods as well as by the Holocene sea level rise. During the Elsterian and the Saalian glaciations, the glacial maximum reached the Central German Basin; thus the Elbe-Weser region was covered entirely by ice several times (Ehlers 2011). During this period, massive glaciofluvial sediment bodies were deposited, mainly consisting of sand and gravel with an average thickness of 20 m. In addition, during the Elsterian glaciation, tunnel valleys up to 500 m deep were formed by subglacial erosion and filled predominantly by glaciofluvial sand. The valleys are covered by so called 'Lauenburger Clay' (Kuster and Meyer 1979). Due to sea level rise in the following Eem interglacial, the sea extended into the region bounded by the moraine ridges of Land Hadeln (Höfle et al. 1985). The maximum of the Weichselian glaciation did not reach the Elbe-Weser region (Streif and Köster 1978) so that the area has only been influenced by proglacial processes. Brackishmarine sediments of the marshlands were deposited in the Holocene (Streif and Köster 1978).
Pleistocene and partly Pliocene sands, and the occasional presence of gravel, form the upper aquifer of the Elbe-Weser region with local variations. The uppermost groundwater storage in the Elbmarsch consists of Weichselian fluvial sand and glaciofluvial sand of the Drenthe glacial (Saalian). In the Unterweser Marsch and the Geest area, Pleistocene sand forms the upper aquifer. The tunnel valley aquifer consists of Elsterian glaciofluvial sediments. The hydraulic contact between the tunnel valley aquifer and the surrounding Pliocene sand is discontinuous. The Lauenburger Clay in the region is the most important dividing layer between the different aquifers (Elbracht et al. 2016). Hydraulic conductivity values of the Saalian aquifer at the region range between 10 -4 and 10 -3 m/s (Reutter 2011). Clayey Holocene sediments in the marshlands work as a protective layer for the groundwater, whereas the protection potential of groundwater in the moraine ridges is very low. Groundwater in the marshlands has a high salinity in contrast to the moraine ridges, where groundwater is fresh (Elbracht et al. 2016).
The study area lies between the moraine ridges of Altenwalde in the west and the Elbe River in the east, including the Hadelner Marsch. The area is characterized by low groundwater recharge (between 51 mm/year and 150 mm/ year) in the marshlands and high groundwater recharge (between 101 mm/year and 400 mm/year) with a relatively high groundwater level (> 10 m ASL) (NIBIS ® Kartenserver 2018) in the moraine ridges. The groundwater level in the study area follows the topography. In general, the groundwater level varies between -4 m and 28 m ASL.

Source of salinity in the aquifer
The primary source of salinity in the aquifer is saline water intrusion from the North Sea (average Cl concentration of North Sea water close to the coast is about 16,000 mg/l) and from the Elbe River (average Cl concentration between 9000 mg/l and 40 mg/l (WSV 2019)). The chloride concentration of the Elbe River decreases upstream. Salt also accumulated in the shallow aquifer below Unterweser Marsch, Hadelner Marsch, and Elbmarsch due to Holocene transgressions (leading to saline water infiltrating vertically from the flooded area before the dykes were built). At Stade, a salt structure exists that contains a different geochemical signature than the saline water originating from the sea (Rahman et al. 2018). The contribution of this salt structure to the groundwater salinity at Stade is yet to be quantified.

Methodology
The following steps were performed to estimate the Cl distribution from HEM data: (1) analysis of borehole information, (2) formation factor estimation, (3) assigning petrography classes to HEM resistivity at different depths, (4) estimation of groundwater electrical conductivity (ECw), (5) establishment of a relation between ECw and Cl, (6) 2-D  (7) manual inspection, validation and adjustment of Cl distribution. All data were obtained from several secondary sources ( Table 2). The data were checked for possible errors, consistency, then were processed after corrections. Python codes were written for step 1 to step 7 (except step 6) to perform the tasks quickly and efficiently. The codes are general and, therefore, applicable to any other comparable HEM survey area (unconsolidated sedimentary rocks) after minor modification. Golden Software Surfer 8 was used for 2D interpolation. A brief description of each step is given below.

Analysis of borehole information
Borehole description is considered as one of the primary sources of information to understand the horizontal and vertical distribution of the petrography at the study area. A total of 8622 boreholes with depths in the range of 1-2000 m were obtained from Geological Survey of Lower Saxony (LBEG) for the study area. Borehole ID, coordinates, depth, and corresponding petrographic symbols were extracted from the database. The petrographic symbols were converted to petrography descriptions using the geology key-symbol database of LBEG (LBEG 2018), e.g. T-U (clay to silt) was converted to 'silty clay'. Altogether, 13 petrography classes were defined this way and listed in Table 3, together with values for porosity and, if available, formation factor (FF). Afterwards, the petrography information was assigned to every meter depth of the borehole.

Formation factor estimation
The straightforward way to estimate the formation factor FF would result from its definition. FF is the ratio of ECw to the bulk conductivity of the fully saturated rock (Archie 1942) or, expressed in terms of resistivity, the quotient of the resistivity of the fully saturated rock (Rb) and the resistivity of the pore water (1/ECw). As information on ECw is not spatially available, we considered estimating FF from grain size data and porosity. FF was estimated using the following empirical equation (Taylor Smith 1971): where Φ is the porosity of the sediment sample. This equation is for clay-free sediment, and is valid for Φ ≤ 0.6, which is within the limit for the porosities of our samples (except peat). Sediment porosity was estimated using the following equation (Vukovic and Soro 1992): where U is the coefficient of uniformity and defined as the ratio of d60 to d10 from grain size analysis. A total of 647 grain size analysis data at different depths from 83 boreholes were obtained from the LBEG, BGR, and LIAG databases. The spatial distribution of grain size data is mainly limited to the Geest area in the western part of the study area, only a few samples are available in the marshlands. FF values of silt, clay, coarse sand, gravel and peat were obtained from the literature (e.g., De Louw et al. 2011; Goes et al. 2009). In addition, if petrography classes whose formation factor or porosity were not included in our database or the literature, we considered FF values from similar petrography classes (see details in Table 3).
Assigning formation factor to HEM data at different depths BGR conducted airborne surveys covering the entire study area, simultaneously recording electromagnetic, magnetic, and radiometric data (Siemon et al. 2020). The electromagnetic data used here were acquired in 2000 for the Cuxhaven area , in 2001 for the Bremerhaven area (Siemon et al. ), in 2004 for the Hadelner Marsch area (Siemon et al. 2017) and in 2008-2010 for the Glückstadt area (Steuer et al. 2013), all available at Produktcenter BGR (2017). After inversion of the frequency domain electromagnetic data by BGR using Levenberg-Marquardt inversion techniques (Sengpiel andSiemon 2000, Siemon et al. 2009), the layered vertical distribution of the bulk resistivity Rb was obtained at sampling distances of 4 m along the flight line. The number of layers varies between 5 and 6, however for the Glückstadt survey, smooth multilayer mod-els were available. Considering the relative stationarity of the groundwater Cl distribution, time differences between the surveys are acceptable (De Louw et al. 2011).
For each borehole in the study area with petrographic information, the nearest HEM data point was identified and the corresponding petrography class and FF were assigned at each meter depth of the HEM resistivity model (Rb data).

Estimation of groundwater electrical conductivity (EC w )
Using the information of Rb [Ωm] and FF at each meter depth of the selected HEM data points, ECw [S/m] was derived by: This relation is valid for clay-free sediments. Drilling results show that clayey sediments mainly occur in the Geest area with the exception of the clayey near-surface layer in the marshlands.

Relation between EC w and Cl
ECw is usually used as a salinity indicator. Since chloride ions are the main constituents in groundwater and in saline sediments that directly affect ECw, chloride concentrations in groundwater can be an indicator of the degree of saline intrusion (Peinado-Guevara et al. 2012).
Lower Saxony Water Management, Coastal Protection and Nature Conservation Agency (NLWKN) provides a database of measured Cl content and ECw values mea-sured at 25°C. The dataset used in this paper contains 1936 ECw-Cl pairs from 67 groundwater-monitoring stations. We applied the ECw temperature compensation relation (Sorensen and Glass 1987) defined by: Here, ECwt is the ECw at 10°C, the mean groundwater temperature at the study area. In Eq. 4, α is a temperature compensation factor which was recommended as 0.0191 by Clesceri et al. (1998).
In many studies (e.g., De Louw et al. 2011) a strong correlation between ECw and Cl has been found. We applied a linear regression analysis to determine the correlation between ECw and chloride concentration Cl (mg/l) from groundwater samples: where m and c are regression coefficients. The high electrical conductivity of groundwater in the study area (> 2000 µS/cm) is mostly dominated by chloride and sodium ions, although the piper diagram presented by Rahman et al. (2018) shows that Ca(HCO3)2 and Mg(HCO3)2 are also present in the groundwater. To avoid uncertainties involved in bicarbonate measurement and to consider that the electrical conductivity is not linearly related to the Cl concentration (De Louw et al. 2011), we separated the entire data set into three classes. Therefore, three different regression analyses were performed. The classifications based on the range of groundwater conductivity were freshwater (EC25 < 500 μS/cm), brackish water (500 µS/cm < EC25 < 2000 µS/cm) and saline water (EC25 > 2000 µS/cm). Using the regression analysis results, we converted all estimated ECw (µS/cm), derived from HEM resistivity ( m) and borehole data, to Cl concentration (mg/l).
After the regression analysis, the relations were validated against the measured EC25 and Cl from 350 samples from direct push tests in 2017 and 2018 provided by LBEG.

2-D interpolation of chloride data
Instead of interpolating FF values, we chose chloride values to interpolate all data in the entire area because the chloride distribution in the groundwater is smoother than the FF distribution which is rather heterogeneous due to spatial and vertical variability of the petrography. The interpolation was based on a 100 m by 100 m grid using ordinary Kriging method with Auto Fit option. The Cl concentration was interpolated at different depths by Surfer 16.0 in the 64-bit version 16.0.330 (www.goldensoftware.com).

Manual inspection, validation and adjustment of Cl distribution
The groundwater in the Geest area is considered freshwater, i.e., the Cl concentration in groundwater is less than 250 mg/l, the threshold set for drinking water in Germany (BMJV 2001).
The initial interpolation results showed relatively higher Cl concentration in some areas due to clayey material, e.g. the Lauenburger Clay on top of a tunnel valley. To avoid such visible and unrealistic errors, we performed some manual validation (e.g., checking resistivity and borehole petrography information) and adjustment for the Geest area. The fresh saline water interface defined by the resistivity limit of 30 Ωm (Repsold 1990;Deus et al. 2015) is found at the Geest-Marsch boundary. First, we deleted resistivity values below Rb = 30 Ωm in the Geest area. Second, we decided to set the maximum Cl concentration to 250 mg/l in the Geest area. Finally, the new dataset was re-interpolated. The interpolation was done for each meter depth below the surface of the vertical extent covered by the survey area.

Results and discussion
Formation factor for the petrography classes Table 3 shows the porosity and formation factor of each petrography class. For completeness, formation factors of coarse sand, silt, sandy clay, and clayey sand were taken from De Louw et al. (2011) or Goes et al. (2009) who investigated comparable sites in the Netherlands. In addition, because of the invalidity of Eq. 1 for cohesive materials (in which the porosity must be below 0.6, and surface conductivity which is common to clayey material is neglected), the apparent formation factor of clay and peat were also taken from the literature. Due to a lack of available information on the formation factor of silty clay and clayey silt, we assigned formation factors of clay and silt, respectively.

Resistivity and petrography classes
Measured resistivities (Rb) from HEM data were grouped in a box plot, which shows that the median resistivity value of most petrography classes increases with the increase of grain size, except in the case of peat (no grains) and gravel (Fig. 2a). The highest median resistivity occurs in coarse sand with only 70 Ωm, and the lowest resistivity is in clay, around 20 Ωm. Medium sand shows the largest spectrum of resistivity values ranging between 10 Ωm and 140 Ωm, while silty clay has the smallest range, between 15 Ωm and 35 Ωm. The resistivity occurrence of some selected classes shows the dominance of sandy material (Fig. 2b). The most Fig. 2 a Box plot of resistivity from HEM data for each petrography class (increase of grain size on x-axis); b resistivity occurrence from HEM data for selected petrography classes (sand comprises fine, medium and coarse sand) Abb. 2 a Box-Plot des spezifischen elektrischen Widerstands aus HEM-Daten für jede Petrografieklasse (auf X-Achse mit zunehmender Korngröße aufgetragen); b Widerstandshäufigkeit für ausgewählte Petrografieklassen abgeleitet aus HEM-Daten (Sand beinhaltet Fein-, Mittelund Grobsand) frequently occurring resistivities for the sand fraction are in the range of 80 m, which corresponds to the freshwatersaturated sands in the Geest areas. The accumulated values in the lower resistivity range are interpreted as saline watersaturated sands in the marshlands; the minimum at about 8 Ωm is suggested to define the threshold to the saline affected sands. This relatively straightforward interpretation of the histogram is consistent with Wiederhold and Ullmann (2017); Siemon et al. (2019) also discuss such low resistivity threshold values for fresh-saline groundwater interfaces. We used Eqs. 6-8 to convert ECw derived from HEM resistivity to Cl concentration. Fig. 4 shows the different steps to obtain the chloride distribution from the Rb distribution, and the resistivity chloride conversion, omitting resistivities below 30 Ωm in the Geest area and additionally setting the maximum chloride value to 250 mg/l in the Geest area. Fig. 5 shows the final maps for different elevations down to -60 m ASL. The high chloride concentration in the groundwater occurs mainly in the marshlands close to the seaside; the high salinity is also distributed along the Elbe River and Oste River. The saline area expands with the decrease in elevation, as higher chloride concentrations are also present in the lower part of the seaside. However, below -50 m ASL, the resistivity information from the HEM survey as well as lithology information is scarce and cannot provide good interpolation results.

Discussion on errors and uncertainties
We identified the following uncertainties that might contribute to the erroneous Cl distribution.

Analysis of borehole information
The deduction of petrography classes and formation factors from well information depends on the quality of the sediment profiles. Different drilling methods and thus associated uncertainties in sediment description need to be considered. In hydraulic circulation drillings, for example, the fine grain content is often underestimated. Another problem can be the simplification in the description of the sediment profiles, which is based on the subjective view of the responsible geologist.

Formation factor estimation
The sampling procedure for grain size analysis can lead to errors and uncertainties in follow-up estimations. Grain size analyses provide a quick overview of the grainsize distribution and the porosity. Because of the very high spatial variation in the composition of Quaternary sediments, the results are only valid for the sampling location. In addition, the method does not take the alignment and the compaction of the sediment into account, and results could differ from real conditions. Porosity is calculated from grain size analysis using the coefficient of uniformity (Eq. 2). If the content of fine grain material is underestimated, as some of the fine material might be lost in the mud, the uniformity (d60/d10) will be too low leading to an overestimate of porosity. This results in an underestimate of the formation factor. Additionally, the porosity exponent in Eq. 1 is set to 1.5 following Taylor Smith (1971), but for unconsolidated sediments it can be in the range of 1.3-2. Values from Goes et al. (2009) imply porosity exponents (cementation factors) of 1.9 for coarse sand and 1.7 for gravel (Table 3), while other authors, e.g., Repsold (1989), used the modified Humble formula (FF = 0.81 ϕ -2 ).
Comparing formation factors for a porosity of 0.4 calculated after Taylor Smith, Goes and Repsold leads to an uncertainty of about 25%. Therefore, the uncertainty in poros-ity estimation from uniformity can be neglected against the uncertainty resulting from porosity exponent errors.
For drillings where no porosity data are available, the formation factors were assigned to the layers according to the particular petrography class. Here the same uncertainty must be taken into account.
The applied method for formation factor estimation is valid for clay-free material only. In this study, areas with the possibility of clay occurrence (e.g. in the Geest area) were excluded. Even in the marshlands the clay content is low compared to the sand fraction as shown in the histograms of Fig. 2b.

HEM data and inversion
As measured data always have errors, even after thorough data processing, the derived physical parameters are also biased, particularly if they are the result of an optimization process being prone to over-interpretation, e.g., the inversion of HEM data to resistivity-depth models. The resulting 1D models used in this study represent only one solution explaining the data. Such a solution is always a simplified image of the subsurface. Comparison studies (e.g., Delsman et al. 2018 andKing et al. 2018), however, demonstrated that the majority of models were able to reveal the principle layered resistivity features. A further challenge occurs if the subsurface is strongly heterogeneous in the lateral direction. In this case, 1D models may fail to explain the data, and 2D/3D modelling or inversion is required. Fortunately, this point is not essential here.

Assigning petrography class to HEM data at different depths
At the HEM resistivity model location, we assigned petrography from the nearest available (proximity max. 100 m) borehole. Considering the geologic and petrographic heterogeneity at the study area, assigned petrography might not be accurate. Therefore, the formation factor, which is related to petrography, might also be inaccurate. Subsequently, erroneous results can be generated. In a regional aquifer system, it is challenging to resolve the spatial and vertical heterogeneity in an automated process. Manual inspection definitely reduces the uncertainty but time and resources would be vital.

Estimation of groundwater electrical conductivity (EC w )
ECw is estimated from the bulk resistivity and formation factor. Since both quantities are linearly related, the 25% uncertainty of the formation factor leads to the same amount of uncertainty in ECw. Additionally, a uniform groundwater temperature of 10 ºC is assumed that also corresponds to the average temperature of groundwater samples at the study area (average value of 204 samples is 10.1 ºC and ε is 0.92). The spatial variation in temperature is between 13.1 ºC and 8.4 ºC, which is non-uniform and might contribute uncertainty in the ECw estimation by +8% to -4%.

Relation between EC w and Cl
We considered that Na+ and Cl-dominate the electrical conductivity through the entire range of ECw data available, although the linear relation may only be valid at higher ECw (> 2000 µS/cm) (De Louw et al. 2011). Below this limit, other ions, such as bicarbonate ion (HCO3 -) due to the presence of Ca(HCO3)2 and Mg(HCO3)2, also might influence the ECw (Goes et al. 2009). The Cl concentration might therefore be overestimated.

Interpolation error
In addition to the above problems, interpolation error can also be considered as one of the main issues that generates error in the chloride distribution. In our case study, the data points were neither evenly nor densely distributed, and therefore, the interpolation methods might produce some deviation from the true chloride distribution.

Conclusions
Despite many uncertainties, the 3D chloride distribution in the aquifer of the study area was determined by HEM data in combination with borehole data, grain size analysis and laboratory results on the relation of CL content and electrical conductivity from groundwater samples. A threshold value of 8 Ωm for saline-affected sandy glacial sediments was found by a statistical analysis of HEM resistivity data.
The database was substantially enlarged by the HEM data, nonetheless, the full data density of HEM has not yet been exploited. HEM resistivity data were taken only to determine formation factors in the surrounding boreholes. The chloride concentrations obtained at the borehole locations were interpolated. Resistivity variations from between boreholes as seen in the HEM data, indicating lateral changes in the geology, were not taken into account. For future studies, one has to consider interpolation methods that take the geology into account to be able to use the full HEM data set for interpolation of Cl concentrations.
The estimated ECw can be useful to generate a densitydriven groundwater model. The ECw data could be used as an initial condition when the generation of initial condition considering a historical change in natural recharge and other boundaries (such as seaside boundary, historical flooding) is not possible.
The fresh-saline groundwater interface derived in this study can be validated and investigated in detail with densely-spaced surface measurements, e.g. with 2D or 3D ERT.
The grouping of borehole data in terms of petrography classes as applied in this study was used to determine a depth profile of porosities and formation factors. This method can be extended to other petrophysical characteristics such as, for example, hydraulic conductivity and verified by borehole NMR. This can lead to an extended database for groundwater modelling.
Acknowledgements The TOPSOIL project is co-funded by the Interreg North Sea Region Programme (J No 38 2 27 15). We thank NL-WKN for providing the database of measured Cl content and ECw values. The manuscript was improved by the comments of two reviewers.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4. 0/.