New versions of REGIS and GeoTOP

[dbrakenhoff]

Received by email:

Beste gebruiker van DGM, REGIS II en GeoTOP,

Graag informeren we je over een nieuwe versie van de ondergrondmodellen DGM, REGIS II en GeoTOP. In de nieuwe versies worden de stratigrafische en hydrogeologische coderingen in overeenstemming gebracht met die van de Stratigrafische Nomenclator. De verandering van de codes betekent dat programmatuur die gebruik maakt van deze modellen aangepast moet worden. We zijn voornemens om de nieuwe versies (DGM v2.2.1, REGIS II v2.2.3 en GeoTOP v1.6.1) op 18 maart 2025 op BROloket en DINOloket te publiceren.

Nieuwe coderingen

Stratigrafische en hydrogeologische eenheden in de BRO-modellen DGM, REGIS II en GeoTOP worden aangeduid met een alfanumerieke code, bijvoorbeeld ‘BX’ voor de Formatie van Boxtel (stratigrafische eenheid in DGM en GeoTOP) en ‘BXk1’ voor de Formatie van Boxtel, eerste kleiige eenheid (hydrogeologische eenheid in REGIS II). Formaties kunnen verder onderverdeeld zijn in lagen en laagpakketten, zo staat ‘NAWA’ voor Formatie van Naaldwijk, Laagpakket van Walcheren, en ‘NAWAZU’ voor de Zuiderzee Laag binnen het Laagpakket van Walcheren.

Wat tot nu toe niet in de coderingen tot uitdrukking kwam, is dat de formaties tot een groep behoren. Deze groepscode is nu toegevoegd aan de bestaande codes. Zo behoren veel formaties in de ondiepe ondergrond tot de Boven-Noordzee Groep, aangeduid met de code ‘NU’. Deze code ‘NU’ wordt dus toegevoegd aan de bestaande codes. De hierboven genoemde codes worden daarmee respectievelijk ‘NUBX’, ‘NUBXk1’, ‘NUNAWA’ en ‘NUNAWAZU’. Je vindt deze codes ook terug in de Stratigrafische Nomenclator. Door deze wijziging is het makkelijker om toelichtende informatie over een stratigrafische en hydrogeologische eenheid te bekijken in de Nomenclator.

Naast het toevoegen van de groepscode zijn er nog een aantal andere wijzigingen. Zo wijzigt van een aantal formaties de formatiecode, bijvoorbeeld ‘MT’ wordt ‘MA’ (Formatie van Maastricht). Verder zijn enkele modeleenheidcodes in GeoTOP herzien, bijvoorbeeld ‘AEC’ wordt ‘ECga’ (Formatie van Echteld, geulgeneratie A). In bijgaand memo vind je meer informatie en een overzicht van alle wijzigingen in de coderingen.

Aandachtspunten bij gebruik van REGIS II en GeoTOP

Voor gebruik in bijvoorbeeld een grondwaterstromingsmodel is het van belang dat de namen van de rasterbestanden van het lagenmodel gewijzigd zijn. De code maakt namelijk deel uit van de bestandsnaam van de rasters. Programmatuur die de rasters gebruikt zal daarom aangepast moeten worden.

Op het 3D voxelmodel van GeoTOP heeft de wijziging van de coderingen geen effect. De geologische eenheden in het voxelmodel worden namelijk met een numerieke, viercijferige code aangeduid (bijvoorbeeld 3100 voor de Formatie van Boxtel). Deze numerieke code blijft ongewijzigd.

Andere wijzigingen in GeoTOP

Voor DGM en REGIS II is de verandering van de codering de enige aanpassing in de nieuwe release. Voor GeoTOP is ook een inhoudelijke verbetering doorgevoerd en bevat de nieuwe release tevens een verbeterde inschatting van de modelonzekerheid van de gemodelleerde geologische eenheden. Ook is een tweetal terugmeldingen op GeoTOP opgelost, namelijk:

• Meldingsdatum 19-02-2024, omschrijving: Het bestand van BXWISIKO_BCC (lagenmodel onderkant BXWISIKO) is incorrect. Dit bestand geeft de onderkant weer van BXWIKO en niet die van BXWISIKO.
• Meldingsdatum 19-02-2024, omschrijving: In het samengevoegde lagenmodel van GeoTOP v01r6 zijn de grids van eenheid NAWO nog van een oudere versie. Deze grids kruisen daardoor met jongere of oudere afzettingen zoals NAWA, NIHO en NIBA.

In de nieuwe versie van GeoTOP zijn de juiste grids opgenomen.

[rubencalje]

On March 18 2018 the layer codes of REGIS and GeoTOP will change. This is described in this document (in dutch):

Toelichting bij naamgevingsreleases DGM REGIS II en GeoTOP.pdf

These changes will probably break all code in which layer-names are used. As soon as the new model versions become available, I think we should release a new NLmod-version, in which we edit the names used within the code of NLmod itself (for example nlmod.read.regis.get_combined_layer_models).

We can also make a method to rename the layers of regis and geotop back to their current layer-names, as a quick fix for existing models. I will request TNO to deliver the tables with layer-name-changes mentioned in the document above,

[rubencalje]

I received an Excel-file with the code-changes. See the following attachment:

Tabellen bij naamgevingsreleases DGM REGIS II en GeoTOP.xlsx

[rubencalje]

So what will we do? There are the following options:

1. We do nothing. For GeoTOP we have programmed the layer names ourselves. nlmod.read.geotop.get_geotop return stratographic unit ids as intergers, that do not change, so nothing changes. For REGIS, we download the layer names when we download the data. These names will be changed since a previous run of your script, and you script and some methods within nlmod will probably fail.
2. We change the geotop names to the new names in nlmod and change the use of regis-names within nlmod (for example replace "HLc" by "NUhlc" as the default value for the parameter geotop_layers in nlmod.read.regis.get_combined_layer_models.
3. We rename the downloaded regis-names to the old ones, as soon as we have downloaded regis-data (in nlmod.read.regis.get_regis). Existing scripts will behave as before (after you have updated nlmod), and you can continue to use the old names.

I am in favor of option 3 for now. We can then change to the new names of regis and geotop (option 2) in a major version of nlmod, at a moment we can choose ourselves.

[ArtesiaWater]

+1 for option 3, with the option to not rename (through some keyword argument or something)

[martinvonk]

nlmod.read.regis.get_regis(extent=[101000, 120000, 412000, 432000])

results in a KeyError: 'AKc'

[rubencalje]

This issue is solved by PR #424, which was merged just now, by implementing option 3 from the comment above. The fix will be included from nlmod version 0.9.3.

I will leave this issue active for now, so people will find it more easy.