DOI: 10.5281/zenodo.5902880
title: shallowgroundwater: estimated zones with shallow groundwater in the Flemish region of Belgium

description: <strong>General</strong> The data source <code>shallowgroundwater</code> is a geospatial dataset of multipolygons that represent the estimated areas, in the Flemish region of Belgium, where the mean lowest groundwater level (MLW; Knotters & Van Walsum, 1997; Van Heesen, 1970) is less than approximately 2 m below soil surface (hence, “shallow” groundwater). We expect groundwater dependent species and communities to be present within these areas. Outside these areas we assume they are groundwater independent. We combined several data sources in order to estimate these areas. <strong>Compilation of the data source</strong> We compiled the dataset through an iterative process of adding specific data sources (referred in the description of attributes below), followed by validation steps based on both the actual presence of groundwater dependent habitat types or regionally important biotopes (Natura 2000 habitat map) and in situ measurements of groundwater levels (Watina+ database). The coverage of these validation data by <code>shallowgroundwater</code> was 96.9% and 98.6% respectively. Most steps to compile the data source were done manually using QGIS. Final steps were done in R; see R-code in the GitHub repository 'n2khab-preprocessing' at commit 1b004e1. <strong>Detailed properties</strong> The data source is a GeoPackage with a single spatial multipolygon layer <code>shallowgroundwater</code> in the ‘Belge 1972 / Belgian Lambert 72’ coordinate reference system (EPSG-code 31370). All attributes are boolean (true/false), each indicating if a polygon was selected from the corresponding data source by applying a set of criteria. Multiple attributes can be true for a given polygon. In order to reduce file size, polygons were dissolved by each unique combination of the values of all attributes. Hence the dataset consists of multipolygons (also multipart polygons) rather than single (part) polygons. The different attributes of this dataset reveal for each polygon (1) the data source(s) we relied on and (2) the selection criteria we applied to judge if the mean lowest groundwater level is less than approximately 2 m below soil surface. As far as possible, we reference each used data source in the description of attributes below. If one is interested in the original polygons of each datasource, selections can be made by consulting the referenced data sources and applying the specified criteria. These are the attributes: <code>geomorph_wcoast</code>: source: Cosyns et al. (2019) description: polygon belonging to geomorphological entities that are expected to harbour groundwater dependent types, and thus to exhibit shallow groundwater levels. It mainly concerns dune slacks, mud flats and salt marshes. selection: <code>"Code" IN ('ms', 'msl', 'sm', 'ss', 'ys', 'ysl', 'yfs') OR ("Code" = 't' and "Subtype" like '%vlakte%')</code> with ms = medium old dune slack, msl = leveled medium old dune slack, sm = mud flat, ss = salt marsh, ys = young dune slack, ysl = leveled young dune slack, yfs = young frontal dune (intruded), (<code>"Code" = 't' AND "Subtype" like '%vlakte%'</code>) = fossil beach <code>anthrop_gwdep:</code> source: <code>soilmap_simple</code> (10.5281/zenodo.3732903) as derived from the digital soil map of the Flemish Region (version <code>soilmap_2017-06-20</code>; 10.5281/zenodo.3387008); <code>habitatmap_terr</code> (10.5281/zenodo.3468948) as derived from the Natura 2000 habitat map of Flanders (10.5281/zenodo.3354381) description: zones located within a 100 m buffer around (almost) everywhere groundwater dependent habitat types (or regionally important biotopes) ánd situated within zones classified as “anthropogenic” areas within the soil map. Within the zones of the soil map that are designated as “anthropogenic”, we lack information on soil characteristics. However, (almost) everywhere groundwater dependent types are present in these zones according to the Natura 2000 habitat map of Flanders, implying shallow groundwater levels. By including 100 m buffer zones around these types, restricted to the anthropogenic zones of the soil map, we consider the combined areas to have (potentially) shallow groundwater levels. So in practice, we first select the habitatmap polygons with (almost) everywhere groundwater dependent types that intersect the anthropogenic soil polygons, we buffer them, and then clip the result by the anthropogenic soil polygons. selection: see https://github.com/inbo/n2khab-preprocessing/pull/61 for the adopted workflow from <code>soilmap_simple</code>: <code>bsm_mo_soilunitype starts with 'O'</code> from <code>habitatmap_terr</code>: (almost) everywhere groundwater dependent types only; list of this category of types is available through n2khab R-package. <code>narrowanthrop_gwdep</code>: source: <code>soilmap_simple</code> (10.5281/zenodo.3732903) as derived from the digital soil map of the Flemish Region (version <code>soilmap_2017-06-20</code>; 10.5281/zenodo.3387008); <code>habitatmap_terr</code> (10.5281/zenodo.3468948) as derived from the Natura 2000 habitat map of Flanders (10.5281/zenodo.3354381) description: narrow zones classified as “anthropogenic” areas within the soil map that include (almost) everywhere groundwater dependent habitat types (or regionally important biotopes). Regarding the anthropogenic soil type polygons, it appears that the narrow ones containing (almost) everywhere groundwater dependent types are interesting to include as a whole as zones with supposed shallow groundwater levels. Hence we select them as a whole instead of selecting buffers around the polygons with (almost) everywhere groundwater dependent types (cfr. <code>anthrop_gwdep</code> and <code>dunes_gwdep</code>). An appropriate algorithm selects meaningful polygons, based on a “thinness” criterion and the fraction of (almost) everywhere groundwater dependent types that are present within the polygons. selection: see https://github.com/inbo/n2khab-preprocessing/pull/61 <code>drainage</code>: source: <code>soilmap_simple</code> (10.5281/zenodo.3732903) as derived from the digital soil map of the Flemish Region (version <code>soilmap_2017-06-20</code>; 10.5281/zenodo.3387008) description: Drainage classification is based on a combination of groundwater depth, soil permeability, presence of impermeable layers, soil depth and topography (see Van Ranst & Sys, 2000). selection: <code>bsm_mo_drain in ('c-d', 'd', 'e', 'f', 'g', 'h', 'i', 'e-f', 'e-i', 'h-i')</code>; these are soils that are at least moderately gleyic or wet. <code>dunes_gwdep</code>: source: <code>soilmap_simple</code> (10.5281/zenodo.3732903) as derived from the digital soil map of the Flemish Region (version <code>soilmap_2017-06-20</code>; 10.5281/zenodo.3387008); <code>habitatmap_terr</code> (10.5281/zenodo.3468948) as derived from the Natura 2000 habitat map of Flanders (10.5281/zenodo.3354381) description: zones located within a 100 m buffer around (almost) everywhere groundwater dependent habitat types (or regionally important biotopes) ánd situated within zones classified as “dunes” areas within the soil map. Within the zones of the Belgian soil map that are designated as “dunes”, we lack information on soil characteristics. However, (almost) everywhere groundwater dependent types are present in these zones according to the Natura 2000 habitat map of Flanders, implying shallow groundwater levels. By including 100 m buffer zones around these types, restricted to the “dunes” of the soil map, we consider the combined areas to have shallow groundwater levels. So in practice, we first select the habitatmap polygons with (almost) everywhere groundwater dependent types that intersect the “dunes” polygons, we buffer them, and then clip the result by the “dunes” polygons. selection: see https://github.com/inbo/n2khab-preprocessing/pull/61 from <code>soilmap_simple</code>: <code>bsm_mo_soilunitype = 'X'</code> from <code>habitatmap_terr</code>: (almost) everywhere groundwater dependent types only; list of this category of types available through n2khab R-package. <code>peat_profile</code>: source: <code>soilmap_simple</code> (10.5281/zenodo.3732903) as derived from the digital soil map of the Flemish Region (version <code>soilmap_2017-06-20</code>; 10.5281/zenodo.3387008) description: variant of the soil profile indicates a superficial peaty cover, mostly on gleyic or permanently water saturated soil with or without profile development (‘(v)’), eventually combined with strong anthropogenic influence (‘(o)’) selection: <code>bsm_mo_profvar in ('(o)(v)', '(v)')</code> <code>peat_substr</code>: source: <code>soilmap_simple</code> (10.5281/zenodo.3732903) as derived from the digital soil map of the Flemish Region (version <code>soilmap_2017-06-20</code>; 10.5281/zenodo.3387008) description: soil substrate (layer underlying superficial layer, and lithologically diverging from it) consists of peat material starting at small (less than 75 cm; ‘v’) or moderate depths (75-125 cm; ‘(v)’), or a combination of the previous (‘v-’) selection: <code>bsm_mo_substr in ('(v)', 'v', 'v-')</code> <code>peat_parentmat</code>: source: <code>soilmap_simple</code> (10.5281/zenodo.3732903) as derived from the digital soil map of the Flemish Region (version <code>soilmap_2017-06-20</code>; 10.5281/zenodo.3387008) description: parent material contains a mixture of at least 30% of peaty material selection: <code>bsm_mo_parentmat = 'v'</code> <code>peat_texture</code>: source: <code>soilmap_simple</code> (10.5281/zenodo.3732903) as derived from the digital soil map of the Flemish Region (version <code>soilmap_2017-06-20</code>; 10.5281/zenodo.3387008) description: soil consists of plain peat material selection: <code>bsm_mo_tex in ('V-E', 'V')</code> <code>phys_system</code>: source: Fysische systeemkaart - Gegeneraliseerde bodemkaart voor Vlaanderen (here available at geopunt.be). Lhermitte & Honnay (1994). description: polygons designated as seepage areas where groundwater is supposed to gather after having infiltrated elsewhere (infiltration areas) and being transported through the landscape (passage areas). The classification of infiltration/passage/seepage areas is entirely based on the soil map polygons of Flanders (10.5281/zenodo.3387008). Each polygon is attributed to one of these classes by considering soil, topographical and hydrological characteristics. In the seepage areas, we consider the groundwater to be shallow below topography. selection: <code>INKWEL = 3</code> <code>zwin</code>: source: Open Street Map, consulted 2021-10-06 by QGIS plugin QuickOSM description: the contour of the Zwin Nature Reserve in the most eastern part of the Flemish coastal area, as mapped in the Open Streep Map on 2021-10-06 selection: <code>leisure = nature_reserve</code> <code>habitat_1130</code>: source: Natura 2000 habitat map of Flanders (10.5281/zenodo.3354381) description: polygons located within estuaries (habitat type <code>1130</code>); these areas are located outside dike infrastructure and harbour several types that rely on tidal dynamics or brackish water, and at the same time are considered (almost) everywhere groundwater dependent, a.o. brackish pioneer habitats, salt marshes and meadows (habitat types <code>1310</code>, <code>1320</code>, <code>1330</code> resp.), reed beds (<code>6430_mr</code>, <code>rbbmr</code>) and riverine Salix forests and shrubs (<code>91E0_sf</code>, <code>rbbsf</code>). selection: <code>Hab1 = '1130' or Hab2 = '1130' or Hab3 = '1130' or Hab4 = '1130' or Hab5 = '1130'</code> <code>gwdepth_coast</code>: source: interpolated groundwater level of point measurements along the Flemish coast, taken from the WATINA database of the Research Institute for Nature and Forest (INBO) description: locations with estimated average lowest groundwater table less than 2.5 m below soil surface (shallow groundwater), based on interpolation of measured groundwater levels in areas along the Flemish coast with sufficient gauge densities; in these coarse sandy areas, groundwater levels are supposed to behave independently of topography, allowing simple interpolation among point measurements selection: calculation of mean yearly spring groundwater level at each gauge (GVG; i.e. average of three groundwater head measurements in meter Ostend height, EPSG:5710, closest to 1st April, but with a minimum time lapse of 14 days, averaged over years) calculation of mean yearly amplitude of groundwater level at each gauge (i.e. difference between min and max water level per hydrological year, averaged over years) inverse distance weighted interpolation of both these datasets (IDW; factor 3, resolution 5 m) for areas with sufficient gauge densities calculation of relative GVG (depth below soil surface) by subtracting soil surface level (Flemish DTM) from GVG, both in meter Ostend height (EPSG:5710) subtract average amplitude from relative GVG as an estimate of lowest groundwater level (GLG) select locations with estimated <code>GLG &gt;= -2.5</code> m located within the coastal dunes ecodistrict (see Flemish ecoregions and ecodistricts) <code>gwdepth_local</code>: source: Batelaan et al. (2012) description: mean lowest groundwater level less than 2 m below soil surface (MLW) for the present day situation. The modelled area is a large military training site for which no information is available in the soil map of Flanders. selection: <code>GLG &lt;= 2</code> m <code>seepage</code>: source: Batelaan & De Smedt (1994) description: area with modelled seepage fluxes of at least 0.8 mm/day. In these areas we inherently consider the groundwater level to be shallow. The extent of the regional groundwater model encompasses the valleys of the rivers Demer, Dijle en Nete in the central part of eastern Flanders selection: <code>gridcode &gt;= 3</code> <code>peat_survey</code>: source: peat pricks executed or compiled by the Research Institute for Nature and Forest (INBO) layer with interpreted presence of peaty zones evaluated by simple measurements of the depth of plain and superficial peaty soil layers at regular intervals (about 20 m), sampled by manually stabbing the peaty layers with a simple stick until the underlying mineral layer is reached (i.e. resistance is felt). Th
description: {"references": ["Batelaan O., De Becker P., El-Rawy M., Herr C., Schneidewind, U. (2012). Doorrekenen van maatregelen voor herstel van vochtige heidevegetaties op het Schietveld van Houthalen-Helchteren via grondwatermodellering. Vrije Universiteit Brussel/Instituut voor Natuur en Bosonderzoek.", "Batelaan O., De Smedt F. (1994). Regionale grondwaterstroming rond een aantal kwelafhankelijke natuurgebieden, Vrije Universiteit Brussel.", "Cosyns E., Bollengier B., Provoost S. (2019). Masterplan en juridische basis voor grensoverschrijdende samenwerking en bescherming als een transnationaal natuurpark van de duinen tussen Dunkerque (Frankrijk) en Westende (Belgi\u00eb). Partim Masterplan. Rapport in opdracht van Agentschap Natuur en Bos, Conservatoire de l'espace littoral et des rivages lacustres, Conseil G\u00e9n\u00e9ral D\u00e9partement du Nord. https://www.natuurenbos.be/sites/default/files/inserted-files/masterplan_flandre_ned20200210def.pdf", "De Saeger S., Guelinckx R., Oosterlynck P., De Bruyn A., Debusschere K., Dhaluin P., Erens R., Hendrickx P., Hennebel D., Jacobs I., Kumpen M., Op De Beeck J., Spanhove T., Tamsyn W., Van Oost F., Van Dam G., Van Hove M., Wils C. &, Paelinckx D. (2020). Biologische Waarderingskaart en Natura 2000 Habitatkaart, uitgave 2020. (Rapporten van het Instituut voor Natuur- en Bosonderzoek; Nr. 35). Instituut voor Natuur- en Bosonderzoek (INBO). https://doi.org/10.21436/inbor.18840851.", "Knotters M. & van Walsum P.E.V. (1997). Estimating fluctuation quantities from time series of water-table depths using models with a stochastic component. Journal of Hydrology 197 (1): 25\u201346. https://doi.org/10.1016/S0022-1694(96)03278-7.", "Lhermitte K., Honnay O. (1994). \"Kartering van het Fysisch Systeem en de Ruimtelijke structuren in Vlaanderen op schaal 1 : 50 000\", Stichting Plattelandsbeleid 1994, Boerenbond, Leuven.", "Provoost S., Van Gompel W. & Vercruysse E. (2020). Beheerevaluatie kust. Eindrapport 2015-2019. Rapporten van het Instituut voor Natuur- en Bosonderzoek 2020 (18). Instituut voor Natuur- en Bosonderzoek, Brussel. https://doi.org/10.21436/inbor.18039583.", "Van Heesen H.C. (1970). Presentation of the seasonal fluctuation of the water table on soil maps. Geoderma 4 (3): 257\u2013278. https://doi.org/10.1016/0016-7061(70)90006-6.", "Van Ranst E., Sys C. (2000). Eenduidige legende van de digitale bodemkaart van Vlaanderen (schaal 1: 20000). Universiteit Gent, Laboratorium voor Bodemkunde, Ghent, 361 p. http://hdl.handle.net/1854/LU-125899."]}