A Brondata, bewerking en download

A.1 Zeespiegelstijging

In de DSR Waddenzee is de jaargemiddelde waterstand voor de zes hoofdstations langs de Nederlandse kust overgenomen uit de Zeespiegelmonitor. In deze studie is de stijging van de zeespiegel nauwkeurig onderzocht en berekend. De berekening en aanvullende en interactieve figuren kunnen worden ingezien via de Jupyter Notebook Zeespiegelmonitorberekening.

A.2 Waterstanden

A.2.1 Brondata

Voor de beschikbare stations in de Waddenzee en aangrenzende delen van de Noordzee worden de tijdseries van zowel de gemeten als de berekende waterstand (astronomische waterstand) ingelezen via de webservice van de Data Distributie Laag van Rijkswaterstaat. Merk op dat er een aantal stations missen in de Waddenzee, zoals vermeld in paragraaf 9.2.

De recente waterstandsmetingen, vanaf 1984, hebben een tijdstap van 10 minuten op alle stations. Voor 1984 is met een lagere frequentie gemeten, met een tijdstap van 1 uur, en voor 1966 is dit zelfs een tijdstap van 3 uur, zie figuur A.1). Voor deze lagere frequenties is de hieronder beschreven methode voor het afleiden van de waterstanden te onnauwkeurig. Voor die periode beschikt RWS wel over een dataset van hoog- en laagwaterstanden, die hiervoor kunnen worden gebruikt. Daarvan kan dan het jaarlijkse gemiddelde worden bepaald. Deze data zijn momenteel nog niet online beschikbaar.

Gebruikte tijdstap in de dataset voor alle stations. Vanaf 1984 is voor alle stations data op elke 10 minuten beschikbaar. Tussen 1966 en 1984 (variërend per station), is dit een tijdstap van een uur. In de jaren voor 1966 is dit zelfs 3 uur.

Figuur A.1: Gebruikte tijdstap in de dataset voor alle stations. Vanaf 1984 is voor alle stations data op elke 10 minuten beschikbaar. Tussen 1966 en 1984 (variërend per station), is dit een tijdstap van een uur. In de jaren voor 1966 is dit zelfs 3 uur.

De getijcomponenten zijn berekend op basis van de waterstationsmeetseries, zoals ontsloten via Open Earth Data.

A.2.2 GHW, GLW en getijslag

De jaargemiddelde hoog- en laagwaters (GHW en GLW) zijn bepaald door in de tijdseries van gemeten waterstanden te zoeken naar het maximum en minimum per getijcyclus, met een onderlinge afstand van minimaal 4 uur. In figuur A.2 is daarvan een voorbeeld te zien. Vervolgens wordt hiervan het jaargemiddelde bepaald. De jaargemiddelde getijslag wordt berekend door het verschil tussen het jaargemiddelde GHW en GLW.

Om trends te kunnen bepalen is vervolgens een lineaire trendlijn toegevoegd die varieert met de 18,6-jarige cyclus. Doordat de frequentie van de 18,6-jarige cyclus bekend is en de tijdserie voldoende lang, kan de software zelf de amplitude en de fase van de 18,6-jarige cyclus fitten.

Gemeten waterstand en hoog- (rood) en laagwater (groen) voor eind maart 2018 (inclusief overgang naar zomertijd van 24 op 25 maart), zoals berekend voor de hele tijdserie.

Figuur A.2: Gemeten waterstand en hoog- (rood) en laagwater (groen) voor eind maart 2018 (inclusief overgang naar zomertijd van 24 op 25 maart), zoals berekend voor de hele tijdserie.

A.2.3 Getijcomponenten en -asymmetrie

Op basis van de gemeten waterstanden wordt met het python package ‘hatyan’ de amplitude en fase van de verschillende getijcomponenten bepaald. Op basis van deze amplitudes en fases kan de amplitudeverhouding M4/M2 en het faseverschil (\(2 \phi M_{2} - \phi M_{4}\)) worden berekend https://github.com/Deltares/hatyan. Voor uitgebreide documentatie wordt verwezen naar de gebruikershandleiding.

Getijcomponenten zijn berekend per station en jaar, en kunnen gedownload worden vanaf watersysteemdata.deltares.nl.

A.2.4 Gemiddeld laag laag water (GLLW) en gemiddeld hoog hoog water (GHHW)

GLLW en GHHW worden berekend als de maandminima en -maxima van het berekend getij per station. Het berekend getij wordt bepaald uit de getijcomponenten, zie sectie A.2.3.

A.2.5 Stormopzet en stormvloedhoogte

De stormopzet (en -afzet) wordt berekend als de gemeten waterstand, verminderd met de astronomische waterstand (het getij). In de systeemrapportage worden de 99,5-percentiel per jaar en het maximum per jaar gebruikt als maat voor storm.

Voor de stormvloedhoogte gebruiken we het 99,5-percentiel per jaar van de gemeten waterstand (dus inclusief het getij).

De astronomische waterstand is berekend per station en jaar, en kan gedownload worden vanaf watersysteemdata.deltares.nl.

A.3 Golven

A.3.1 Brondata

De significante golfhoogte, golfperiode (Tm02) en golfrichting worden rechtstreeks ingelezen uit de Data DistributieLaag. De ingelezen set die gebruikt is voor deze rapportage is beschikbaar via onderstaande downloadknop:

A.3.2 Anomalie van de maandgemiddelde golfhoogte

Als maat voor de seizoensfluctuatie van de golfhoogte wordt de afwijking (anomalie) van de maandgemiddelde golfhoogte ten opzichte van het langjarige gemiddelde gebruikt.

De anomalie is berekend aan de hand van de volgende stappen:

  1. Incomplete jaren verwijderen.

  2. Langjarig gemiddelde voor elk station berekenen als gemiddelde van alle maandgemiddelden.

  3. De afwijking per station en maand berekenen als het verschil tussen maandgemiddelde en langjarig gemiddelde.

  4. De afwijking uitzetten tegen de maand in het jaar, en voor alle jaren in een figuur. De kleurenschaal toont het onderscheid tussen de jaren.

A.4 Zoetwaterafvoeren

A.4.1 Brondata

De daggemiddelde waarden van de drie grote afvoeren vanaf Rijkswateren (Den Oever, Kornwerderzand en Lauwersoog) zijn opgevraagd via Servicedesk Data van RWS. De 10-minuutwaarden bevatten veel gaten in de tijdreeksen en zijn daarom niet gebruikt.

Afvoeren van de Eems en Leda zijn opgevraagd bij NLWKN en toegevoegd aan de Datahuis Wadden Geoserver. De data kunnen ook via onderstaande links direct worden gedownloadt als csv. Om niet ongewild grote datasets te downloaden is deze nu nog beperkt tot de eerste 50 records. Wanneer alle data gewenst zijn, moet “&maxFeatures=50” verwijderd worden uit de URL:

Voor de polderafvoeren van de waterschappen zijn gegevens opgevraagd bij het desbetreffende waterschap en voor deze rapportage beschikbaar gemaakt via de Datahuis Wadden Geoserver. Directe links naar de afvoerdata staan hieronder.

Rijkswaterstaat

De etmaalgemiddelde afvoeren die in deze rapportage zijn gebruikt zijn opgevraagd bij de servicedesk data van Rijkswaterstaat. Deze gegevens zijn lokaal opgeslagen als bevroren versie voor deze rapportage en te downloaden via onderstaande link.

NLWKN (Eems en Leda)

Overzichtskaart van hydrologische meetstations in het stroomgebied van de Eems en West-Duitse Waddenzee.

Figuur A.3: Overzichtskaart van hydrologische meetstations in het stroomgebied van de Eems en West-Duitse Waddenzee.

Waterschappen

A.4.2 Bewerking van afvoeren

In de systeemrapportage worden de daggemiddelde, maandgemiddelde en jaargemiddelde afvoeren weergegeven.

Daarnaast wordt voor de grotere afvoeren ook de gemiddelde zomerafvoer en de gemiddelde winterafvoer berekend. Voor de zomer wordt daarbij gemiddeld over de periode 1 april tot en met 30 september. En voor de winter wordt gemiddeld over de periode van 1 oktober tot en met 31 maart.

A.5 Wind

A.5.1 Brondata

De data betreffen de KNMI uurgegevens voor de vier stations De Kooy, Lauwersoog, Hoorn Terschelling, Vlieland.

A.5.2 Windsterkte en windrichting

Als indicatoren voor de windsterkte worden de windrozen per station van de volledige tijdreeks en het jaarlijkse 25, 50, 75 en 95 percentiel van de windsnelheid weergegeven.

Voor het bepalen van de periodes met hoge windsterkte wordt gekeken naar het aantal achtereenvolgende uren dat de windsnelheid boven een bepaalde grenswaarde komt. In de systeemrapportage worden de grenswaardes >10,8 m/s (Beaufort 6) en >13,9 m/s (Beaufort 7) gebruikt. Dit aantal uren wordt uitgezet in de tijd voor alle windrichtingen en uitgesplitst per kwadrant. Daarnaast wordt ook het aantal uren gesommeerd per jaar weergegeven met onderverdeling per kwadrant, zodat de jaren onderling beter vergeleken kunnen worden.

A.6 Waterkwaliteit

A.6.1 Brondata

MWTL metingen van Rijkswaterstaat worden beschikbaar gemaakt via Waterinfo en zijn ook direct via een webservice te downloaden uit de Datadistributielaag (DDL) (documentatie webservice).

Zwevend stof was tijdens de samenstelling van dit rapport niet beschikbaar uit de DDL, en is daarom apart opgevraagd bij de servicedesk data van RWS. Voor het kunnen combineren van deze data met andere waterkwaliteitsdata is een aparte nabewerking nodig (mapping).

A.6.2 Gemiddelde watertemperatuur (maand, jaar en zomer/winter)

Omdat de MWTL metingen een andere opnamefrequentie in de zomer en de winter hebben, is het van belang eerst de ‘maandgemiddelden’ te bepalen, en deze vervolgens te middelen per jaar. Anders worden de jaargemiddelden beïnvloed door de hogere meetfrequentie in de zomer. (Het woord ‘maandgemiddelden’ staat tussen aanhalingstekens, omdat er in de winterperiode vaak maar 1 meting per maand is.)

De seizoensgemiddelden zijn alleen berekend en weergegeven als er minimaal 5 metingen in het seizoen beschikbaar waren. Het zomerseizoen wordt gerekend vanaf 1 april t/m 30 september, en het winterseizoen vanaf 1 oktober t/m 31 maart.

Zomer- en wintergemiddelden zijn berekend om te zien of er trends in de seizoenen zitten die bij jaargemiddelden weg worden gemiddeld.

A.6.3 Gemiddelde saliniteit (maand, jaar) en anomalie

Omdat de MWTL metingen een andere opnamefrequentie in de zomer en de winter hebben, is het van belang eerst de ‘maandgemiddelden’ te bepalen, en deze vervolgens te middelen per jaar. Anders worden de jaargemiddelden beinvloed door de hogere meetfrequentie in de zomer. (Het woord ‘maandgemiddelden’ staat tussen aanhalingstekens, omdat er in de winterperiode vaak maar 1 meting per maand is.)

Het jaargemiddelde is alleen berekend bij meer dan 10 metingen per jaar.

De anomalie van de saliniteit is de afwijking van de maandgemiddelde saliniteit ten opzichte van het langjarig gemiddelde. De anomalie is berekend aan de hand van de volgende stappen:

  1. Incomplete jaren verwijderen.

  2. Langjarig gemiddelde voor elk station berekenen als gemiddelde van alle maandgemiddelden.

  3. De afwijking per station en maand berekenen als het verschil tussen maandgemiddelde en langjarig gemiddelde.

  4. De afwijking uitzetten tegen de tijd.

A.6.4 Gemiddelde zwevende stof gehalte (maand, jaar) en anomalie

Omdat de MWTL metingen een andere opnamefrequentie in de zomer en de winter hebben, is het van belang eerst de ‘maandgemiddelden’ te bepalen, en deze vervolgens te middelen per jaar. Anders worden de jaargemiddelden beinvloed door de hogere meetfrequentie in de zomer. (Het woord ‘maandgemiddelden’ staat tussen aanhalingstekens, omdat er in de winterperiode vaak maar 1 meting per maand is.)

Voor zwevende stof is de anomalie, de afwijking ten opzichte van het langjarig gemiddelde, berekend voor alle jaren en maanden. Bovendien is deze data log-getransformeerd. De stappen voor de afleiding zijn:

  1. Incomplete jaren verwijderen.

  2. Alle meetwaarden log (e-basis) transformeren.

  3. Langjarig gemiddelde berekenen voor elk station als gemiddelde van alle maandgemiddelden.

  4. De afwijking berekenen per station per maand, als het verschil tussen maandgemiddelde en langjarig gemiddelde.

  5. De afwijking uitzetten tegen de tijd, zowel voor de gehele meetperiode als per maand. De kleurenschaal geeft de absolute waarde van de zwevend stofconcentratie aan. Bij het figuur met de gehele meetperiode achter elkaar is een lijn toegevoegd op basis van het 13-maanden (dus minstens 1 jaar) lopend gemiddelde.

A.7 Morfologie

A.7.1 Brondata bodemligging

De verschillende bodemkaarten van de Waddenzee zijn een samenslag van lodingen en hoogtemetingen over de verschillende kombergingsgebieden waarvan er jaarlijks één bemeten wordt. Zo is elke 6 jaar door Rijkswaterstaat een kaart van de volledige Waddenzee geconstrueerd, op basis van de kaarten van deelgebieden over die periode. Deze zogenaamde vaklodingen zijn opgehaald van de Open Earth Dataserver. De database bevat de vaklodingen data in 20x20 m resolutie. De vaklodingen zijn een combinatie van lodingen van de diepe delen aangevuld met LiDAR opnames van platen en waterpassingen/RTK-DGPS van de hoge delen zoals kwelders. Voor meer informatie over deze metingen verwijzen we naar de desbetreffende webpagina van Rijkswaterstaat. De eerste volledige kaart dateert van 1991, met data over de periode 1985-1991. Let wel, deze kennen missende punten en vooral de kwelders en ondiepe gebieden langs de eilandkusten zijn niet volledig. De serie bodems vanaf 1985 zijn begin 2024 door Deltares opnieuw gecombineerd, aangevuld en gecorrigeerd en ontsloten via Datahuis Wadden. Voor de bodems vóór 1985 geldt dat de lodingen minder systematisch plaats vonden en ontsloten werden, waardoor de samengestelde bodems minder duidelijk terug te voeren zijn op een bepaalde periode. De tabel hieronder geeft een overzicht van de gebiedsdekkende kaarten vanuit de volledige vaklodingencycli en de onderliggende jaren van de afzonderlijke zeegaten.

…1 Marsdiep Eierland Vlie Ameland Friesche ZG Groninger Wad ED
Cyclus 1- WadGrid 1985-1991 1985-1988 1987 1988 1989 1987 1989 1990
Cyclus 2- WadGrid 1991-1997 1991 1993 1992 1993 1994 1995 1996-1997
Cyclus 3- WadGrid 1997-2002 1997 1999 1998 1999 2000 2001 2002
Cyclus 4- WadGrid 2003-2008 2003 2005 2004 2005 2006 2007 2008
Cyclus 5- WadGrid 2009-2014 2009 2011 2010 2011 2012 2013 2014
Cyclus 6- WadGrid 2015-2020 2015 2017 2016 2017 2018-2019 2019 2020
Cyclus 7- WadGrid 2021-present 2021 2023 2022 2023
NA

De visualisatie van de bodemligging van de Waddenzee is interactief beschikbaar gemaakt voor deze Digitale Systeemrapportage. De interactiveit wordt beschikbaar gemaakt met behulp van een r shiny app, die in de rapportage getoond wordt. De bathymetrie wordt hier gepresenteerd als een laag bovenop een Openstreetmap achtergrondkaart. De bathymetrielaag bevat gestapelde zgn. “tiles” van vaklodingen tussen de jaren die gekozen zijn met de slider. De bodemligging kan met en zonder “Hillshading” worden bekeken.

De oorspronkelijke data is gemeten en opgeleverd door RWS voor elk jaar3, en opgeslagen op een thredds server van Deltares als netcdf4. Van elke netcdf is één tiff file per jaar gemaakt door een script5 en in een Google Cloud storage (gs) map gevat. Hierna is van gs een Google Earth Engine (gee) Image Collection van vaklodingen data gemaakt. De r shiny app gebruikt een API service6 die de Image Collection leest, op basis van de gegeven input data.

Voor het berekenen van morfologische indicatoren zijn de vaklodingen als rasters (tiff) samengesteld voor elk jaar voor de hele Waddenzee. De rasters bestaan uit data van het betreffende jaar, gecombineerd met data van alle voorgaande jaren, waarbij alleen de meest recente waarde per pixel gebruikt wordt. Nadeel is dat sommige delen van een raster oudere waarden kunnen bevatten. Het voordeel is dat er voor elk jaar een gebiedsdekkend raster voorhanden is.

A.8 Brondata bodemslibgehalte

Er zijn verschillende historische kaarten met het bodemslibgehalte beschikbaar (Lely ~1890, De Glopper ~1960, Sedimentatlas ~1989, SIBES/SUBES sinds 2008/2019).

De sedimentatlas is digitaal beschikbaar en geeft een volledig dekkende kaart van de bodemsamenstelling in de Waddenzee inclusief het Eems estuarium. Deze data is getoond in deze digitale systeemrapportage en beschikbaar in OpenEarthData.

Sinds 2008 zijn door het NIOZ monsters gestoken van het intergetijdengebied en is ook de korrelgrootteverdeling bepaald. Deze data wordt aangeduid met SIBES. De data is niet vrij beschikbaar en daarom vooralsnog niet opgenomen in deze rapportage.

In ongeveer 2019 is de SIBES dataset gecomplementeerd met data van monsters die zijn genomen in de diepe delen van de Waddenzee, waarmee een ruimtelijk volledig dekkende kaart van de Waddenzee kan worden gemaakt. De data uit de diepe delen worden aangeduid met SUBES. Ook deze data is (nog) niet vrij beschikbaar.

A.8.1 Brondata afbakening kombergingsgebieden

Voor de opdeling in kombergingen zijn polygonen van de Wageningen Marine Research Geoserver gebruikt.

A.8.2 Hypsometrische curves

Hypsometrische curves worden bepaald voor subgebieden (kombergingsgebieden) door de oppervlakte per diepte uit te zetten tegen de diepte.

A.8.3 Ruimtelijk dekkende waterstanden

Voor de bepaling van bepaalde morfologische karakteristieken is het ook belangrijk om een ruimtelijk dekkend beeld te hebben van de waterstanden in het bekken. De gemeten waterstanden geven alleen informatie van de waterstanden op de stations (zie paragraaf A.2).

Om van de gemeten waterstand een ruimtelijk dekkende waterspiegel te reconstrueren, moet rekening worden gehouden met de getijvoortplanting. Numerieke modellen die de waterstanden berekenen, hebben moeite om laagwaterstanden goed te reproduceren bij een acceptabele horizontale celresolutie. Dit komt doordat in een numeriek model de bodem wordt geschematiseerd naar grotere gebieden (rekencellen), waardoor de kleinere geultjes en prielen niet in het rekenrooster zijn opgenomen. Deze geultjes hebben een belangrijke functie in de drainage van de intergetijdengebieden, en daarmee ook op de laagwaterstand in die gebieden.

Voor de (morfologische) zonering van de ondiepe gebieden is er voor gekozen om in eerste instantie uit te gaan van vaste begrenzingen van de hoogteklasses. Daarmee worden niet de veranderingen in getijslag door de tijd of over het gebied meegenomen. Dit betekent dat verandering in arealen alleen worden veroorzaakt door veranderingen in de bodemligging, en niet door een verandering in getijvoortplanting. Een indeling aan de hand van dergelijke begrenzingen geeft weliswaar niet het exacte intergetijdengebied weer, maar levert wel een bruikbare proxy hiervoor. Bovendien geven de vaste waterstandsgrenzen inzicht in veranderingen louter als gevolg van bodemveranderingen, en niet als gevolg van veranderingen in de waterstanden. Voor het systeembegrip is dit een voordeel.

Om nauwkeuriger inzicht te krijgen in de areaalveranderingen van verschillende morfologische eenheden van de kombergingsgebieden, kan de hypsometrische curve worden onderverdeeld in vier verschillende diepteklassen. Omdat hierbij wordt gerekend met vaste begrenzingen van de diepteklassen, en veranderingen in getijslag door de tijd of over het gebied niet worden meegenomen, worden de (ondiepere) klassen aangeduid met ‘proxy’.

We onderscheiden de volgende 4 klassen:

  1. Geul (\(z < -5 m\))

  2. Proxy subgetijdengebied ( \(-5 m < z < -1 m\))

  3. Proxy intergetijdengebied ( \(-1 < z < 1 m\))

  4. Proxy supragetijdengebied (\(z > 1 m\))

Waar z de verticale coordinaat t.o.v. NAP betekent.

A.8.4 Plaathoogte en -volume

Voor de berekening van de plaathoogte wordt het volume boven de NAP -1 m grens berekend, en gedeeld door de oppervlakte. In formulevorm:

\[\begin{equation} \text{Plaathoogte} = \frac{\sum_i(h_i \cdot A_i)}{\sum_i A_i} \text{ if } h\geq -1 \tag{A.1} \end{equation}\]

waarin:

  • A = oppervlakte van 1 rasterpixel (400 \(m^2\))
  • h = hoogtewaarde van een pixel
  • i = de index van de pixels in het raster

Voor het plaatvolume wordt het volume boven de NAP -1 m grens berekend:

\[\begin{equation} \text{Plaatvolume} = \sum_i(h_i \cdot A_i) \text{ if } h\geq -1 \tag{A.2} \end{equation}\]

A.8.5 Geuldiepte en -volume

Voor de berekening van de geuldiepte wordt het volume onder de NAP -5 m grens berekend, en gedeeld door de oppervlakte. In formulevorm:

\[\begin{equation} \text{Geuldiepte} = \frac{\sum_i(h_i \cdot A_i)}{\sum_i A_i} \text{ if } h\leq -5 \tag{A.3} \end{equation}\]

waarin:

  • A = oppervlakte van 1 rasterpixel (400 \(m^2\))
  • h = hoogtewaarde van een pixel
  • i = de index van de pixels in het raster

Voor het geulvolume wordt het volume onder de NAP -5 m grens berekend:

\[\begin{equation} \text{Geulvolume} = \sum_i(h_i \cdot A_i) \text{ if } h\leq -5 \tag{A.4} \end{equation}\]

A.8.6 Lengte van de laagwaterlijn

Voor de berekening van de lengte van de laagwaterlijn is in eerste instantie de totale lengte berekend van de -1 m NAP contourlijnen binnen elk bekken.

A.9 Overzicht alle brondata

Tabel A.1: Overzicht van alle gebruikte data
hoofdstuk brondataset bron
Hydrodynamiek waterhoogte data distributielaag en Deltares OPeNDAP server
Hydrodynamiek waterhoogteberekend data distributielaag
Hydrodynamiek golfhoogte data distributielaag
Hydrodynamiek golfrichting data distributielaag
Hydrodynamiek golfperiode data distributielaag
Hydrodynamiek zoetwaterafvoeren via servicedesk data RWS
Hydrodynamiek hoog- en laagwaterstanden aangevraagd bij servicedesk data RWS
Hydrodynamiek waterhoogte additionele stations via servicedesk data RWS
Meteorologie windsterkte https://www.knmi.nl/nederland-nu/klimatologie/uurgegevens
Meteorologie windrichting https://www.knmi.nl/nederland-nu/klimatologie/uurgegevens
Waterkwaliteit watertemperatuur data distributielaag
Waterkwaliteit saliniteit data distributielaag
Waterkwaliteit zwevend stof via servicedesk data RWS
Morfologie vaklodingen https://opendap.deltares.nl/thredds/catalog/opendap/rijkswaterstaat/vaklodingen_new/catalog.html
Morfologie vakken shape files uit project
Morfologie ecotopenkaart https://www.nationaalgeoregister.nl/geonetwork/srv/dut/catalog.search;jsessionid=0B4745908C502557E461385BDB69FE9F#/metadata/291d329a-a421-41f4-b219-a904756e7aa5
Morfologie slibgehalte sedimentatlas https://opendap.deltares.nl/thredds/catalog/opendap/rijkswaterstaat/sedimentatlas_waddenzee/catalog.html
Algemeen achtergrondkaarten openstreetmap.org via leaflet widget

  1. https://waterinfo-extra.rws.nl/monitoring/morfologie/#h53b817e6-e5cd-49d1-873c-9a0f85d8e416↩︎

  2. https://opendap.deltares.nl/thredds/catalog/opendap/rijkswaterstaat/vaklodingen_new/catalog.html↩︎

  3. https://github.com/openearth/eo-bathymetry/blob/master/notebooks/Vaklodingen2EE.ipynb↩︎

  4. https://github.com/openearth/hydro-engine-service/blob/master/hydroengine_service/main.py↩︎