abcdefgh

Morfologische gevolgen van rivierverruimende maatregelen langs de Merweden

Opzet en toepassing van een SOBEK-model

Afstudeeropdracht A.C.S. Mol

abcdefgh

Morfologische gevolgen van rivierverruimende maatregelen langs de Merweden

Opzet en toepassing van een SOBEK-model

8 december 2003

Afstudeeropdracht Arjan Mol

Afstudeercommissie: - Dr. Ir. J.S. Ribberink - Dr. Ir. C. M. Dohmen-Janssen - Dr. Ir. Z.B. Wang - Ir. E. Snippen

Voorwoord

......

Dit verslag is het resultaat van mijn afstudeeronderzoek voor de studie Civiele Techniek & Management aan de Universiteit Twente. Het betreft een onderzoek naar de morfologische effecten van rivierverruimende maatregelen in de Merweden. De studie is uitgevoerd aan het Rijksinstituut voor Integraal Zoetwaterbeheer en afvalwaterbehandeling (RIZA) te Dordrecht.

Ik wil een aantal mensen bedanken voor hun steun tijdens dit onderzoek. In de eerste plaats mijn begeleider bij het RIZA, Edwin Snippen, bij wie ik dagelijks kon aankloppen voor vragen, problemen en adviezen. Verder gaat dank uit naar Marjolein Dohmen-Janssen en Zheng Bing Wang voor de hulp en adviezen via de e-mail en naar Jan Ribberink voor de adviezen tijdens de vergaderingen met de examencommissie. Tenslotte wil ik nog de medewerkers van RIZA Dordrecht bedanken voor de hulp die zij mij boden bij problemen. De afstudeerstage bij dit instituut verliep aangenaam dankzij de informele werksfeer, maar ook dankzij de expertise binnen dit instituut.

Dordrecht, November 2003

Arjan Mol

i Morfologische gevolgen van rivierverruimende maatregelen

ii Morfologische gevolgen van rivierverruimende maatregelen

Inhoudsopgave

......

Voorwoord i

Inhoudsopgave iii

Samenvatting vii

1 Inleiding 1 1.1 Achtergrond 1 1.2 Doelstelling 2 1.3 Onderzoeksvragen 2 1.4 Opbouw rapport 3

2 Gebiedsomschrijving 5 2.1 Inleiding 5 2.2 Waterbeweging NDB 6 2.3 Sedimenttransport NDB 7 2.4 Morfologie NDB 8

3 Modelopzet 11 3.1 Inleiding 11 3.2 Morfologische Modellering in SOBEK 12 3.3 Schematisatie 13 3.3.1 Uitbreiding Baseline schematisatie 14 3.3.2 Aanpassing bodemprofielen 15 3.4 Invoer gebiedsspecifieke data 16 3.4.1 Invoer gegevens korreldiameters 16 3.4.2 Invoer gegevens baggerwerk 17 3.4.3 Randvoorwaarden 17 3.5 Modelinstellingen 18 3.6 Initiële Run 19

4 Calibratie 21 4.1 Inleiding 21 4.2 Calibratiestrategie 21 4.2.1 Calibratiedata 21 4.2.2 Calibratievariabelen 22 4.3 Effect aanpassing transportformule 25 4.3.1 Engelund & Hansen 25 4.3.2 Meyer-Peter-Müller 25 4.3.3 Discussie 26 4.3.4 Keuze 27 4.4 Effect aanpassingen calibratiefactor 27 4.4.1 Modelaanpassingen 27 4.4.2 Resultaten 27 4.4.3 Discussie 28 4.4.4 Keuze 29 4.5 Effect aanpassing korreldiameter 29

iii Morfologische gevolgen van rivierverruimende maatregelen

4.5.1 Modelaanpassingen 30 4.5.2 Resultaten 30 4.5.3 Discussie 31 4.5.4 Keuze 31 4.6 Effect aanpassing Chézy waarden 31 4.6.1 Modelaanpassingen 31 4.6.2 Resultaten 32 4.6.3 Discussie 33 4.6.4 Keuze 33 4.7 Uiteindelijke keuzes 33 4.7.1 Modelinstellingen 33 4.7.2 Resultaten 33 4.7.3 Discussie 34 4.7.4 Mogelijkheden van het model 36

5 Effecten rivierverruimende maatregelen 37 5.1 Inleiding 37 5.2 Maatregelen 37 5.2.1 Steurgat 37 5.2.2 Rivierverruiming knelpunt Gorinchem 38 5.3 Theoretische analyse maatregelen 40 5.3.1 Kwalitatieve bepaling effecten 40 5.3.2 Kwantificering effecten 43 5.3.3 Conclusies 45 5.4 Resultaten SOBEK 45 5.4.1 Inleiding 45 5.4.2 Nulsituatie 46 5.4.3 Steurgat 48 5.4.4 Rivierverruiming knelpunt Gorinchem 49

6 Conclusies en aanbevelingen 55 6.1 Conclusies 55 6.2 Aanbevelingen 57

Referenties 59

Bijlage I Modellering waterbeweging SOBEK 65

Bijlage II SOBEK schematisatie takken NDB 67

Bijlage III Beschrijving Baseline 69

Bijlage IV Hydraulische calibratie model 71

Bijlage V Korreldiameters NDB 75

Bijlage VI Baggergegevens onderzoekstakken 79

Bijlage VII Bodemliggingen onderzoekstakken 81

Bijlage VIII Aanpassingen korreldiameter 83

Bijlage IX Controle hydraulische calibratie 85

iv Morfologische gevolgen van rivierverruimende maatregelen

Bijlage X Korte termijn effecten maatregelen 87

Bijlage XI Bodemontwikkeling nulsituatie 91

v Morfologische gevolgen van rivierverruimende maatregelen

vi Morfologische gevolgen van rivierverruimende maatregelen

Samenvatting

......

Het Noordelijk Deltabekken (NDB) is het estuarium van Rijn en Maas. Voor het NDB liggen er tal van beheersvragen op het gebied van ecologisch herstel, veiligheid en morfologische ontwikkeling van het systeem. Het huidige kennisniveau voor wat de morfologische modellering van het NDB betreft is onvoldoende om aan te kunnen geven wat de consequenties voor morfologische ontwikkelingen zijn bij een toekomstige verandering in het waterbeheer en/of bij nieuwe inrichtingsvraagstukken. Het laatste aspect speelt in dit onderzoek een rol. Het project ‘Ruimte voor de Rivier’ is opgezet om ten eerste de veiligheid in het Nederlandse rivierengebied te waarborgen. Toekomstige klimatologische veranderingen stellen nieuwe eisen aan de riviersystemen. Herinrichtingsplannen van bepaalde riviertrajecten zijn opgesteld om het rivierensysteem aan de geëiste veiligheid te laten voldoen. Naast het veiligheidsaspect speelt ook de ruimtelijke kwaliteit van het rivierengebied een belangrijke rol in de voorgestelde maatregelen. Aan de daadwerkelijke uitvoering van de herinrichtingsplannen gaat een lang besluitvormingsproces vooraf. Hierin worden onder meer de maatregelen beoordeeld op bijkomende effecten, waarbij dankbaar gebruik gemaakt wordt van allerhande modellen. Voor wat de morfologische effecten ten gevolge van ingrepen in het NDB betreft schiet het huidige modelinstrumentarium tekort. In 2001 is een poging ondernomen door Wang [Wang, 2001] om een morfologisch SOBEK model van het NDB te calibreren. Een betrouwbare calibratie was echter niet mogelijk door het gebrek aan voldoende gebiedsspecifieke gegevens en tekortkomingen in de schematisaties.

De doelstelling van deze studie luidt: Het opzetten en calibreren van het morfologische NDB model met behulp van SOBEK teneinde een tweetal maatregelen in het kader van het project ‘Ruimte voor de Rivier’, die worden toegepast in de Nieuwe en Boven Merwede, te onderzoeken op hun grootschalige morfologische effecten.

De calibratie van het model is afgebakend tot de Merwede takken en de Waal. Het gebruikte model is een morfologisch SOBEK model van het gehele NDB, dat gebaseerd is op het hydraulische SOBEK model voor het NDB. Tijdens de opzet van het model zijn gebiedsspecifieke data ingevoerd, zoals korrelgegevens en baggergegevens. Verder zijn nieuwe verbeterde bodemschematisaties ingevoerd voor een aantal takken, waaronder de Nieuwe en Beneden Merwede. Tevens zijn de randvoorwaarden gespecificeerd die gebruikt zijn bij de calibratie. Aan de benedenranden zijn morfologisch bepalende getijden opgelegd, aan de bovenranden van het model daggemiddelde afvoeren. Als morfologische randen zijn gefixeerde bodemliggingen gebruikt.

vii Morfologische gevolgen van rivierverruimende maatregelen

Calibratie van het model heeft plaatsgevonden over de periode 1990 – 2000. Hierbij is ten doel gesteld het model zodanig te calibreren dat het een juiste voorspelling betreffende de grootschalige morfologische ontwikkeling gerealiseerd kon worden. Onder grootschalige morfologische ontwikkeling wordt verstaan de totale volumeverandering van de bodem tijdens de periode 1990 – 2000 per riviertak. Aangezien er in de gebruikte baggerdata een grote onnauwkeurigheid schuil gaat, zijn de marges waarbinnen het model de grootschalige morfologie moet beschrijven afhankelijk gesteld van de onnauwkeurigheid in de baggercijfers. De gebruikte calibratievariabelen zijn: keuze transportformule, calibratiefactor transportformule, korreldata en Chézy waarden. De onderzochte transportformules zijn Engelund & Hansen en Meyer- Peter-Müler. Na een serie modelruns waarin beide formules gebruikt zijn, is een keuze gemaakt tussen de twee transportformules. Hieruit kwam naar voren dat de formule van Engelund & Hansen de beste keuze was voor verdere calibratie. Tijdens de verdere calibratie is de ingevoerde korreldiameter aangepast op de Waal en de Merwede takken. Verder zijn de Chézy waarden op de Beneden Merwede enigszins aangepast. Het aldus verkregen model is in staat om de grootschalige morfologische ontwikkeling conform de vooraf gestelde eisen te beschrijven, uitgezonderd voor de Beneden Merwede. De beste resultaten zijn gehaald voor de Nieuwe Merwede.

Bij het voorspellen van de effecten van de voorgestelde maatregelen in het kader van ‘Meer ruimte voor de Rivier’ geldt dat deze altijd in relatieve zin beschouwd zijn, dat wil zeggen in vergelijking met de autonome situatie, waarin geen maatregelen worden toegepast. In het onderzoek zijn twee maatregelen bestudeerd. De eerste betreft het openzetten van het Steurgat, waarbij water tijdens extreme afvoer vanuit de Nieuwe Merwede via het Steurgat afgevoerd kan worden richting de Biesbosch. De andere maatregel heeft betrekking op de vernauwing in de Boven Merwede ter hoogte van Gorinchem. Hier is een zomerbedverbreding in combinatie met de optimalisatie van een nevengeul aan de noordoever toegepast.

Beide maatregelen zijn geschematiseerd in SOBEK en de effecten zijn afzonderlijk bepaald. De gehanteerde simulatieperiode bedraagt 50 jaar. Aangezien het model gecalibreerd is met gebruikmaking van baggerwerkzaamheden is ook tijdens de bepaling van de effecten de baggeroptie toegepast. Hiervoor is een over 10 jaar gemiddeld baggercijfer gebruikt. Verder zijn dezelfde instellingen en randvoorwaarden toegepast als tijdens de calibratie.

De maatregel ‘Steurgat’ brengt geen grote morfologische effecten met zich mee. De maatregel is namelijk zodanig ontworpen dat deze slechts bij extreme afvoeren (met een overschrijdingsfrequentie van één maal per jaar) actief is. Op de Waal, de Boven Merwede en de Beneden Merwede zijn de effecten verwaarloosbaar. Op de Nieuwe Merwede zijn ten opzichte van de nulsituatie sedimentatiegolven zichtbaar met toppen van zo’n 12 cm. In vergelijking met de bodemontwikkelingen die zich in dezelfde tijdsperiode voordoen in de nulsituatie, zijn deze effecten aanzienlijk kleiner. De maatregel zorgt verder voor een

viii Morfologische gevolgen van rivierverruimende maatregelen

bodemstijging op lange termijn (50 jaar) van zo’n 4 cm op een groot deel van de Nieuwe Merwede.

De maatregelen nabij Gorinchem hebben veel grotere effecten op de morfologie tot gevolg. Dit is hoofdzakelijk te wijten aan de zomerbedverbreding. Deze maatregel is bij alle afvoeren merkbaar, terwijl de maatregel ‘Steurgat’ alleen bij extreme afvoeren ingezet zal worden. De bodem zal ter plaatste van het tracé van de zomerbedverbreding 2 meter boven de bodem in de nulsituatie komen te liggen. Verder ontstaat een erosiegolf na de toepassing van de ingrepen, die zich in benedenstroomse richting voortplant. De golf heeft aanvankelijk een diepte van -2 m ten opzichte van de nulsituatie, maar naar gelang de golf zich voortplant neemt de diepte af en de breedte toe. De golf plant zich ook voort op de Beneden Merwede en op de Nieuwe Merwede, waar de golf verder afzwakt door diffusie. Na 50 jaar heeft de golf het eind van beide takken nog niet bereikt. Door de grote effecten van alleen de zomerbedverbreding zijn de effecten van de nevengeul te verwaarlozen.

ix Morfologische gevolgen van rivierverruimende maatregelen

x Morfologische gevolgen van rivierverruimende maatregelen

1 Inleiding

......

1.1 Achtergrond

Het Noordelijk Deltabekken (NDB) is het estuarium van Rijn en Maas. Voor het NDB liggen er tal van beheersvragen op het gebied van ecologisch herstel, veiligheid en morfologische ontwikkeling van het systeem.

Huidige ontwikkelingen ten gevolge van klimaatveranderingen [KNMI, 2003] zoals zeespiegelstijging en verhoogde rivierafvoeren stellen nieuwe eisen aan de riviersystemen. Verdere dijkverbetering past echter niet goed in de integrale benadering van de hoogwaterproblematiek. Meer ruimte voor de rivier is een belangrijk aspect van de hedendaagse visie op waterbeheer. Tegelijkertijd worden er ook steeds hogere eisen aan de ruimtelijke kwaliteit van het Nederlandse rivierengebied gesteld. Met de nieuwe aanpak van de veiligheidsproblematiek in het rivierengebied ontstaat een koppeling met het beleid op het gebied van de ruimtelijke ordening, waarin gepleit wordt voor een verbetering van de ruimtelijke kwaliteit van het rivierengebied. De nadruk wordt hierbij gelegd op natuurontwikkeling in het rivierengebied.

Bovengenoemde aspecten vragen om een herinrichting van de huidige riviergebieden. In het kader van de PKB (Planologische Kern Beslissing) ‘Ruimte voor de rivier’ worden de mogelijkheden van rivierverruiming onderzocht [Startnotitie MER, 2002]. Voor het NDB ligt de prioriteit hierbij in en rond de Boven Merwede, de Beneden Merwede en de Nieuwe Merwede.

Inmiddels zijn tal van rivierverruimende maatregelen ontworpen. Deze maatregelen worden naast de hydraulische effectiviteit ook onderzocht op effecten op de morfologie van het riviersysteem. Inschatting van deze effecten wordt gedaan op basis van expert-judgement. Model berekeningen zouden deze inschattingen kunnen versterken. Echter, een adequaat morfologisch model is voor het NDB nog niet beschikbaar.

In juni 2000 is een begin gemaakt met het opzetten van twee morfologische modellen voor het NDB: SOBEK-MOR en ESTMORF [Wang, 2000]. Dit laatste model is een semi-empirisch 1D model, het SOBEK model daarentegen is proces gebaseerd 1D model. Uit dit onderzoek kwamen naast een aantal problemen met beide modellen ook de mogelijkheden voor morfologische modellering van het NDB naar voren. In 2001 is het onderzoek vervolgd door een eerste calibratie uit te voeren met het SOBEK-MOR model voor het NDB. Problemen met het model en het ontbreken van gegevens maakte een volledige en betrouwbare calibratie onmogelijk [Wang, 2001].

1 Morfologische gevolgen van rivierverruimende maatregelen

Inmiddels zijn aanvullende gegevens beschikbaar gekomen. Hiermee is een calibratie van het model wel mogelijk. Deze studie onderzoekt de morfologische effecten van rivierverruimende maatregelen in de Boven Merwede en de Nieuwe Merwede met behulp van een morfologisch SOBEK model, dat hiervoor eerst gecalibreerd zal worden.

1.2 Doelstelling

De doelstelling van deze studie luidt:

Het opzetten en calibreren van het morfologische NDB model met behulp van SOBEK teneinde een tweetal maatregelen in het kader van het project ‘Ruimte voor de Rivier’, die worden toegepast in de Nieuwe en Boven Merwede, te onderzoeken op hun grootschalige morfologische effecten.

De maatregelen zijn gepland rond het splitsingspunt van de Merweden. De morfologische effecten van deze maatregelen zullen zich voornamelijk manifesteren op de omringende takken. De calibratie beperkt zich dan ook tot de Waal en de Merwede takken, die in het vervolg aangeduid worden als de onderzoekstakken. Wel wordt met de hele NDB SOBEK schematisatie gerekend. Dit met het oog op toekomstige uitbreiding van de calibratie. Hierbij wordt de meeste nadruk gelegd op de Boven en de Nieuwe Merwede. De rivierverruimende maatregelen die in dit onderzoek bestudeerd worden zijn immers op deze takken gepland. Met grootschalige morfologische effecten wordt bedoeld de totale volumeverandering van de rivierbodem voor één bepaalde onderzoekstak, ofwel de sedimentbalans per tak. Hier zal het model ook op worden gecalibreerd volgens bepaalde eisen.

1.3 Onderzoeksvragen

De onderzoeksvragen bij bovenstaande doelstelling luiden: 1. Wat zijn de morfologische effecten van een tweetal maatregelen in het kader van het project ‘Ruimte voor de Rivier’? 2. Wat zijn de karakteristieke instellingen en waarden van parameters voor het SOBEK-MOR model van het Noordelijk Deltabekken?

Deelvragen hierbij zijn: 1.1 Wat is de autonome ontwikkeling van het gebied? 1.2 Welke termijn wordt onderzocht? 1.3 Welke maatregelen worden onderzocht? 1.4 Wat zijn de effecten van de maatregelen volgens SOBEK? 1.5 Zijn de door SOBEK voorspelde effecten overeenkomstig de verwachting?

2.1 Waarop wordt het model afgeregeld? 2.2 Met welke calibratieparameters wordt het model afgeregeld?

2 Morfologische gevolgen van rivierverruimende maatregelen

2.3 Welke calibratie-eis wordt gehanteerd? 2.4 Voor welke takken is een betrouwbare calibratie mogelijk? 2.5 Welke problemen treden er op tijdens de calibratie? 2.6 Wat is hiervan de vermoedelijke oorzaak?

1.4 Opbouw rapport

Dit rapport is als volgt ingedeeld. In het volgende hoofdstuk is een beschrijving gegeven van het onderzoeksgebied. De waterbeweging, sedimentbeweging en morfologische ontwikkeling komen hierin naar voren. In hoofdstuk 3 wordt de opbouw van het model beschreven. Het huidige model wordt hier uitgebreid met aanvullende gegevens zoals gebiedsspecifieke data. Na de opbouw van het model wordt het model gecalibreerd. Hoofdstuk 4 geeft een beschrijving van deze calibratie en de uitkomsten daarvan. In hoofdstuk 5 wordt het gecalibreerde model toegepast in het kader van het project ‘Ruimte voor de rivier’. In dit hoofdstuk worden een tweetal rivierverruimende maatregelen onderzocht op de morfologische effecten in de onderzoekstakken. Naast een analytische analyse van deze effecten zal het model gebruikt worden om deze effecten te beschrijven. In het laatste hoofdstuk zijn enkele conclusies en aanbevelingen opgesomd naar aanleiding van deze studie.

3 Morfologische gevolgen van rivierverruimende maatregelen

4 Morfologische gevolgen van rivierverruimende maatregelen

2 Gebiedsomschrijving

......

2.1 Inleiding

Het huidige Noordelijk Deltabekken (figuur 2.1) is het estuarium van Rijn en Maas. De oostelijke begrenzing van het gebied wordt gevormd door de plaatsen tot waar de getij-invloed merkbaar is (Hagestein, Tiel, Lith)...... Figuur 2.1 Het Noordelijk Deltabekken

Noordzee Gouda el Hoek van Holland ss IJ e Hagestein N h ie c uw s e d W n a a te ll Lek rw o e Rotterdam H H g a r te lk an aa ieuwe l N Maas Tiel N o o r Haringvlietsluizen d n M neden ove erw # Be Me B ede W ui rwede aa Sp O l u de Maas H A ar Dordrecht e in D d f gv e g Lith lie o w e t r r d d e Biesbosch am

t s d M e c e

h M w aas e u as ie Ma K N

i rgsch olla l Be e H ndsch Maa Diep s Amer

#

Volkeraksluizen

k kera Vol

0 5 10 15 Kilometers N Schaal (A4) 1 : 350.000

De westelijke begrenzing ligt in zee, daar waar het zoete rivierwater volledig is vermengd met het zeewater (Maasmond en buiten de Haringvlietsluizen). Van een bekken is pas sprake sinds de afsluiting van het Hollandsch Diep en Haringvliet door de aanleg van de Volkeraksluizen (1969) en de Haringvlietsluizen (1970).

Kenmerkend voor het huidige Noordelijk Deltabekken zijn de vele vertakkingen en kortsluitingen die de riviertakken met elkaar maken. In het algemeen hebben de riviertakken een weinig variërende, regelmatige en min of meer rechthoekige dwarsdoorsnede met een eveneens regelmatige bodemgeometrie zonder plaat/geulconfiguraties. In het Hollandsch Diep en Haringvliet komen nog wel geulen en platen voor, relicten uit de tijd van voor de afsluiting. In de meer bovenstrooms gelegen delen zijn er de typische verschijningsvormen van zomerbedding en uiterwaard; hierbij treden verschillen op tussen stroomvoerende en kombergende dwarsdoorsneden. Globaal de helft van het totale kombergende oppervlak bevindt zich in het Hollandsch Diep/Haringvliet. Daarnaast kan in het gebied komberging in de vorm

5 Morfologische gevolgen van rivierverruimende maatregelen

van “dode takken” en havens onderscheiden worden (Calandkanaal, Hollandsche IJssel, havens in Europoort, langs de Nieuwe Maas en de Afgedamde Maas). Een overzichtskaart van de onderzoekstakken is gegeven in figuur 2.2. Hierin is ook de kilometrering aangegeven.

...... Figuur 2.2 Overzichtskaart onderzoekstakken

2.2 Waterbeweging NDB

De waterbeweging in het gebied wordt bepaald door het binnendringen van het getij vanuit zee, de aanvoer van water door de rivieren en het lozingsprogramma van de Haringvlietsluizen. Door de dichtheidsverschillen tussen het zoete rivierwater en het zoute zeewater komt in delen van het gebied een situatie met een zoute onderlaag en zoete bovenlaag voor, waarvan de stroomrichtingen soms tegengesteld zijn. Ook de wind heeft invloed op de waterbeweging in het gebied. Opwaaiing op zee kan de afvoer van water naar zee bemoeilijken en ook binnen het gebied kan opwaaiing de waterbeweging beïnvloeden.

Het getij is het sterkst op de riviertakken langs de noordrand van het gebied en op het westelijk deel van de Oude Maas. Van daaruit dringt het door naar het zuidelijk en oostelijk deel. In een groot deel van het gebied treedt in het algemeen viermaal per dag een omkering van de stroomrichting op. In het overgangsgebied is dit proces echter afhankelijk van de rivierafvoer. Bij hoge afvoeren zal geen omkering van de stroomrichting optreden. Slechts op de bovenstroomse riviertakken is sprake van een uniforme stroomrichting. De rivierafvoeren van de Rijn en de Maas variëren met de seizoenen. De afvoer van de Maas is sterker weersafhankelijk dan die van de Bovenrijn en qua grootte geringer (10 à 20 % van de Bovenrijnafvoer). Wel zijn beide afvoeren redelijk aan elkaar gecorreleerd vanwege de correlatie in regenval in de stroomgebieden. Er kan wel een verschil van een aantal dagen optreden in het tijdstip waarop de afvoertoppen Nederland bereiken.

De restafvoerverdeling in het Noordelijk Deltabekken wordt met behulp van het spuiregime van de Haringvlietsluizen beïnvloed. In het algemeen wordt gestreefd naar een constante restafvoer over de Waterweg van ongeveer 1500 m3/s (de restafvoer is het gemiddelde

6 Morfologische gevolgen van rivierverruimende maatregelen

debiet over één of meerdere getijperioden). De wind (met name tijdens stormcondities) heeft invloed op de waterstand en kan daarmee de afvoerverdeling in het Noordelijk Deltabekken verstoren.

2.3 Sedimenttransport NDB

De sedimenten die in het gebied worden aangetroffen zijn afkomstig van de rivieren en vanuit de zee. Voor het sedimenttransport vanuit de rivieren is vooral de Waal van belang. De Lek en de Maas hebben respectievelijk bij Hagestein en Lith een stuw waardoor sedimenttoevoer bij lage afvoeren voor een groot deel wordt tegengehouden. Bij hogere rivierafvoeren worden deze stuwen geheven. Daardoor neemt vooral de bijdrage van de Lek aan het sedimenttransport belangrijk toe. In tabel 2.1 worden de door van Dreumel (1995) berekende hoeveelheden slib en zand die gemiddeld jaarlijks via de rivieren worden aangevoerd, gegeven op basis van de sedimentbalans 1982 – 1992.

...... Waal Lek Maas Tabel 2.1 Hoeveelheden slib en zand die gemiddeld Slib (106 ton/j) 2,40 0,55 0,45 per jaar op de rivierranden worden 6 aangevoerd. Zand (10 ton/j) 1,05 0,30 0,20

Meer dan de helft van het via de rivieren aangevoerde sediment komt terecht in het zuidelijke gedeelte van het Noordelijk Deltabekken (Nieuwe Merwede, Amer, Hollandsch Diep en Haringvliet), een deel sedimenteert in het noordelijk gebiedsdeel en het restant wordt afgevoerd naar zee (tabel 2.2). Vanuit zee wordt, via de Maasmond, jaarlijks gemiddeld 3,50.106 ton slib aangevoerd en 2,30.106 ton zand. Dit mariene sediment wordt voornamelijk afgezet in de Maasmond, op de Nieuwe Waterweg en in de Rotterdamse havens.

...... Maasmond Haringvliet Zeeranden Tabel 2.2 Hoeveelheden zand en slib die gemiddeld + - netto sluizen totaal jaarlijks op de zeeranden worden Slib (106 ton/j) 3,50 0,70 +2,80 -0,21 2,59 aangevoerd (+) cq afgevoerd (-), alsmede 6 de resulterende aanvoer. Zand (10 ton/j) 2,30 0,30 +2,00 -0,05 1,95

De sedimenten zijn ruwweg te onderscheiden in het grofkorrelige zand en het fijnkorrelige cohesieve slib. Vanwege een relatieve grote valsnelheid wordt zand voornamelijk als bodemtransport getransporteerd. Bij hogere snelheden (als gevolg van hogere afvoeren) kan zand echter ook in gesuspendeerde vorm getransporteerd worden. Tussen de waterbeweging en het zandtransport bestaat een vrij directe relatie met de korreldiameter als belangrijke parameter.

Slib bestaat voor een belangrijk deel uit organische stof. De hoeveelheid organische stof kan echter sterk in de tijd variëren, waardoor eveneens de slibeigenschappen veranderen. Onder invloed van fysisch/chemische en/of biologische oorzaken kunnen slibdeeltjes samenklonteren (flocculeren) en als grotere deeltjes gaan deelnemen aan de transportprocessen. Vanwege de vrij lage valsnelheid kan het

7 Morfologische gevolgen van rivierverruimende maatregelen

slib vrij gemakkelijk in suspensie blijven en kan daardoor horizontaal over grote afstanden worden verplaatst.

Het netto sedimenttransport van zand en slib is het samenspel van erosie, verticale en horizontale verplaatsing, flocculatie, sedimentatie en consolidatie waarbij de processen elkaar onderling kunnen beïnvloeden. Zowel ten aanzien van zand als slib geldt dat erosie en sedimentatie kunnen leiden tot veranderingen in de waterdiepte en de hydraulische ruwheid en daarmee tot veranderingen in de waterbeweging in het systeem.

2.4 Morfologie NDB

Het verbeteringsprogramma voor de Waterweg (verdieping van de Waterweg en aanleg havens, 1957 - 1987) en de afsluiting van het Volkerak (1969) en het Haringvliet (1970) betekenden een trendbreuk in de morfologische ontwikkeling van het Noordelijk Deltabekken. Het Noordelijk Deltabekken is na deze werken een door Rijn en Maas gevoed reservoir, waarin het getij slechts kan doordringen via de Nieuwe Waterweg (Nieuwe Maas) en Oude Maas. De morfologische ontwikkeling van het Noordelijk Deltabekken in de eerste 10 jaar naar de afsluiting van het Haringvliet (1971-1982) is samengevat door Allersma (1988). Zowel in de noordrand als in de zuidrand traden in deze periode afzettingen van slib en zand op. De sedimentatie van zand en slib in de noordrand (Lek, Nieuwe Maas, Nieuwe Waterweg, Hollandsche IJssel) werd gecompenseerd door nautisch baggerwerk. De sedimentatie in de zuidrand (Bergsche Maas, Amer, Nieuwe Merwede, Haringvliet, Hollandsch Diep) zorgde echter voor verdere verontdieping. Hierbij werd zowel zand als slib afgezet. In het westelijk middendeel van het Noordelijk Deltabekken (Spui, Oude Maas en Noord) was echter sprake van erosie. Van Dreumel (1995) beschrijft de morfologische ontwikkeling van het Noordelijk Deltabekken in de periode1982-1992. De bodem in het grootste gedeelte van de zuidrand is in die periode door slibafzettingen verder verondiept. In de rest van het bekken is de bodem of verdiept (Nieuwe Waterweg en het middengedeelte-west als gevolg van erosie en extra baggerwerk) of nauwelijks veranderd (door middel van baggerwerkzaamheden op diepte houden van de rest van de noordrand en het middengedeelte-oost).

De scheiding tussen zandige en slibrijke bodems kan in het Noordelijke Deltabekken per locatie sterk verschillen. In het algemeen wordt zand voornamelijk in gebieden met een sterke stroming van het water aangetroffen. Dit zijn onder andere de relatief snelstromende rivieren als de Lek en de verschillende Merweden, maar ook de monden van de beide estuaria en oude riviertakken in het middengebied zoals de Oude Maas, Spui, Dordtsche Kil en Noord, De bodem bestaat voornamelijk uit slib in gebieden die gekenmerkt worden door relatief lage stroomsnelheden, zoals achter in de Hollandsche IJssel, in havens, de Biesbosch, de Amer, het Haringvliet, het Hollandsch Diep en de Afgedamde Maas.

8 Morfologische gevolgen van rivierverruimende maatregelen

...... Tabel 2.3 Riviertak D50 bovenstrooms D50 benedenstrooms D50 bovenstrooms en benedenstrooms per tak NBD [mm] [mm] (2002). Maas (Lith)-Bergsche Maas 0,60 0,48 Amer 0,37 0,13 Hollandsch Diep-Haringvliet 0,11 0,12 Waal (Rossum)-Boven Merwede 1,01 0,63 Beneden Merwede 0,83 0,35 Nieuwe Merwede 0,41 0,28 Dordtsche Kil 0,46 0,24 Oude Maas 0,32 0,29 Spui 0,33 0,17 Noord 0,41 0,40 Lek (vanaf Hagestein) 0,47 0,41 Nieuwe Maas 0,43 0,03 Nieuwe Waterweg 0,18 0,49

In 2002 zijn van vrijwel alle takken van het gehele Noordelijk Deltabekken bodemmonsters onderzocht om zodoende de korrelgrootte verdelingen te bepalen. Hieruit zijn verschillende kenmerkende korreldiameters te bepalen, zoals de D50, D10, D90, D35 en de D65. In tabel 2.3 is een overzicht gegeven van de mediane korreldiameter (D50) voor de belangrijkste takken in het NDB. Op vrijwel alle takken vindt er in benedenstroomse richting verfijning van de korrels plaats, het zogenaamde ‘downstream fining’. Op het traject Hollandsch Diep – Haringvliet blijft de D50 vrijwel constant en op de Nieuwe Waterweg neemt de D50 zelfs toe, als gevolg van zand dat vanuit zee wordt afgezet benedenstrooms in de Nieuwe Waterweg.

Ook in de tijd kan de sedimentsamenstelling van een waterbodem veranderen. Kenmerkend was bijvoorbeeld de grote hoeveelheden slib die gedurende de beginperiode na sluiting van het Haringvliet in de Nieuwe Merwede werden afgezet. In deze periode was vanwege het wegvallen van het getij het dwarsprofiel van de Nieuwe Merwede te groot. Hierdoor trad een gemiddeld lage stroomsnelheid in de Nieuwe Merwede op, waardoor slib in de bedding tot afzetting kon komen. Nadat het dwarsprofiel van de Nieuwe Merwede als gevolg van deze slibsedimentatie verkleinde en de gemiddelde stroomsnelheid toenam, werden deze sliblagen met zand afgedekt.

In het algemeen kan worden gesteld dat de bodem in het zuidelijk deel van het Noordelijk Deltabekken ten opzichte van de rest van het gebied stijgt, waardoor het bekken kantelt. De rivierbeddingen in grote delen van het Noordelijk Deltabekken zijn zich nog steeds aan het aanpassen aan de gewijzigde situatie. Dit zal ook nog gedurende langere tijd zo blijven. Met name de zuidrand zal nog verdere verontdieping laten zien. In gedeelten van de Nieuwe Maas, de Hollandsche IJssel, de Nieuwe Merwede, de Bergsche Maas en de Afgedamde Maas is de bodemligging redelijk stabiel. Deze ontwikkeling zal in de toekomst gevolgen kunnen hebben voor de verdeling van water in het Noordelijk Deltabekken en daarmee voor het sedimenttransport en de morfologische ontwikkeling.

9 Morfologische gevolgen van rivierverruimende maatregelen

10 Morfologische gevolgen van rivierverruimende maatregelen

3 Modelopzet

......

3.1 Inleiding

In 2000 is een eerste poging gedaan om een morfologisch SOBEK model van het NDB te calibreren [Wang, 2000]. Hieruit kwam naar voren dat de calibratie niet mogelijk was door het ontbreken van voldoende gebiedsspecifieke gegevens en door onvolkomenheden in de schematisaties van bepaalde takken. Het model is daarom nog niet geschikt als hulpmiddel bij het bepalen van de morfologische effecten ten gevolge van rivierverruimende maatregelen.

Dit hoofdstuk beschrijft de opzet en verbetering van het model door de invoer van de nieuwe data. Onderstaand schema laat zien welke nieuwe data er zijn ingevoerd en welke gegevens en instellingen er uit het model van Wang [2001] zijn overgenomen.

Nieuwe invoer: Model Wang: - schematisering Nieuwe - basis schematisatie Merwede, Beneden Merwede - morfologische en Biesbosch randvoorwaarden - korreldata - benedenstroomse - baggerdata hydraulische - bovenstroomse hydraulische randvoorwaarden randvoorwaarden - numerieke parameters - instellingen sedimentverdeling splitsingspunten - keuze transportformule

nieuw model

Dit hoofdstuk is als volgt opgebouwd: in paragraaf 3.2 wordt SOBEK- MOR beschreven. Ter verduidelijking van de werking van SOBEK- MOR, worden de vergelijkingen beschouwd waarop het morfologische model gebaseerd is. In paragraaf 3.3 wordt de modelschematisatie van het NDB toegelicht. Sinds het onderzoek van Wang [2000] zijn nieuwe schematisaties voor enkele belangrijke takken ontworpen. De invoer van deze nieuwe schematisaties is een belangrijke verbetering van het model, aangezien de verouderde schematisaties in het model van Wang [2001] voor problemen zorgden. Naast de nieuwe schematisaties zijn ook de resultaten beschikbaar van een grootschalig onderzoek in het NDB naar korreldiameters van het bodemmateriaal. Verder zijn voor een groot aantal takken de baggerwerkzaamheden in de afgelopen 20 jaar bekend. Al deze nieuwe gegevens vormen een waardevolle aanvulling op het morfologische

11 Morfologische gevolgen van rivierverruimende maatregelen

model van Wang, zodat een morfologische calibratie nu wel mogelijk is. Paragraaf 3.4 geeft een overzicht van deze gegevens. Vervolgens zijn in paragraaf 3.5 de tijdens de calibratie gehanteerde modelinstellingen gegeven. Afsluitend is in paragraaf 3.6 een beschrijving gegeven van de eerste run met het model. Hierin worden enkele problemen, die zorgen voor instabiliteit, geconstateerd en opgelost.

3.2 Morfologische Modellering in SOBEK

SOBEK is een computerprogramma waarmee de water- en zoutbeweging en daaraan gerelateerde processen als sedimenttransport en morfologie in een systeem van open waterlopen één-dimensionaal kunnen worden gemodelleerd. In deze paragraaf wordt een korte beschrijving gegeven van de werking van het morfologische deel van SOBEK. De modellering van de waterbeweging is opgenomen in bijlage I. De informatie is afkomstig uit de SOBEK technical reference manual [SOBEK Help Desk, 2000].

Morfologische veranderingen worden bepaald met behulp van de continuïteitsvergelijking voor bodemmateriaal: ∂A ∂S s −=−s (vergelijking 3.1) ∂∂tx lat

2 met: As Sedimenttransporterende dwarsdoorsnede [m ]; S Sedimenttransport door de dwarsdoorsnede inclusief porievolume [m3/s]; 2 slat Laterale sedimenttoevoer inclusief porievolume [m /s].

Tezamen met de continuïteitsvergelijking en de bewegingsvergelijking voor water vormt deze vergelijking een set die door SOBEK numeriek opgelost wordt. Hierbij hanteert SOBEK de volgende werkwijze: eerst wordt de waterbeweging bepaald op basis van een vaste dwarsdoorsnede die tijdens de voorgaande rekenstap is bepaald. Vervolgens worden de morfologische aanpassingen bepaald op basis van de waterbeweging en wordt de dwarsdoorsnede aangepast.

Het sedimenttransport S wordt als volgt bepaald:

SW= scf s (vergelijking 3.2)

met: Ws Sedimenttransporterende breedte [m];

cf Calibratiefactor transportformule (paragraaf 4.2.2) [-]; s Sedimenttransport per eenheid van breedte [m2/s].

Met behulp van een aantal mogelijke transportformules kan s bepaald worden. Deze transportformules beschrijven het sedimenttransport op basis van lokale omstandigheden zoals bodemruwheid, stroomsnelheid en korrelgrootte. Concreet bepaalt SOBEK als volgt de morfologische verandering: SOBEK bepaalt het verschil in sedimenttransport tussen twee rekenpunten (gridpunten). Rekening houdend met de laterale toe- of afvoer van sediment, bepaalt SOBEK op basis hiervan de verandering

12 Morfologische gevolgen van rivierverruimende maatregelen

van sedimenttransporterende dwarsdoorsnede, wat neerkomt op een verandering in bodemhoogte.

Voor de verdeling van sediment bij splitsingen bestaat de mogelijkheid om hiervoor verschillende relaties te gebruiken. Naar aanleiding van het vooronderzoek [Wang, 2000] is hier gekozen om de machtsrelatie voorgesteld door Wang et al. [1995] te gebruiken. k 1−k S1 ⎛ Q1 ⎞ ⎛ B1 ⎞ = ⎜ ⎟ ⎜ ⎟ (vergelijking 3.3) S 2 ⎝ Q2 ⎠ ⎝ B2 ⎠ met: 3 Sn sedimenttransport op tak n [m /s] 2 Qn debiet op tak n [m /s]

Bn breedte op tak n [m]

De waarde van k bepaalt de stabiliteit van splitsingen. Bij kleine waarden van k is een splitsing instabiel en zal er op den duur alleen één benedenstroomse tak overblijven. Bij grote waarden van k is een splitsing stabiel en zullen alle benedenstroomse takken open blijven.

3.3 Schematisatie

In bijlage II is een overzicht opgenomen van alle takken en knooppunten in het NDB model. De geometrie van de takken is geschematiseerd door het definiëren van bodemprofielen. In de huidige modelschematisatie is gebruik gemaakt van twee typen schematisaties. Het merendeel van de takken is geschematiseerd gelijk aan het hydraulische model voor het NDB, dat in 2000 is ontwikkeld [van der Linden en van Zetten, 2002]. Hierbij is er elke 2 km een bodemprofiel gedefinieerd. Voor enkele bovenstroomse takken, zoals Waal, Boven Merwede en Lek, is echter gebruik gemaakt van nieuwe schematisaties (op Baseline1 gebaseerd), met bodemprofielen om de 500 m. Uit het onderzoek van Wang [2001] is gebleken dat deze nieuwe schematisaties zorgen voor stabielere en nauwkeurigere modelberekeningen. Uitbreiding van de Baseline-schematisaties naar de benedenstroomse takken is één van de aanbevelingen die Wang doet naar aanleiding van zijn onderzoek. Gezien de scope van dit onderzoek is het vooral voor de Nieuwe Merwede en de Beneden Merwede wenselijk om hier ook met Baseline-schematisaties te werken. In bijlage III is een korte beschrijving opgenomen van Baseline. Figuur 3.1 geeft een overzicht van de stroomvoerende oppervlakken voor de onderzoekstakken gedurende een bepaalde afvoer. Vooral op de Beneden Merwede en Nieuwe Merwede nemen de stroomvoerende oppervlakken aanzienlijk toe in benedenstroomse richting.

1 Voor de aanmaak van dwarsprofielen voor het SOBEK model wordt Baseline gebruikt. Baseline is een softwarepakket waarmee schematisaties tbv de invoerbestanden voor SOBEK gestructureerd kunnen worden uitgevoerd. Baseline is een ARC/INFO-database waarmee geografische informatie van een rivierengebied wordt verwerkt tot een netwerkstructuur met doorstroomprofielen, bergend oppervlak, zomerkaden, enz.

13 Morfologische gevolgen van rivierverruimende maatregelen

...... Figuur 3.1 4000 )

Stroomvoerende oppervlakken langs de 2 m

vier onderzoekstakken bij een Waal-afvoer (

3 k a van ca 1500 m /s. l 3000 v r e p p

d o 2000 en oer v 1000 oom r t S 0 913 923 933 943 953 963 973 983 kmr

Beneden Merw ede Nieuw e Merw ede Waal en Boven Merw ede

3.3.1 Uitbreiding Baseline schematisatie Voor de Merwede takken, de Biesbosch takken en de Bergsche Maas, de Amer en een deel van het Hollandsch Diep zijn inmiddels ook baseline schematisaties beschikbaar [Hartman et al, 2001], gebaseerd op bodempeilingen uit 1999.

...... Figuur 3.2 Uitbreiding Baseline schematisatie. De SOBEK-vakken zijn weergegeven en ook de hoogteligging. Op basis hiervan maakt Baseline de nieuwe schematisaties voor elk vak.

Figuur 3.2 geeft een overzicht van de nieuwe schematisaties van deze takken. Deze schematisaties zijn ingevoerd in het morfologische model. Vervolgens is het model opnieuw hydraulisch gecalibreerd omdat de aanpassingen in de bodemschematisaties ook de waterbeweging beïnvloeden. De calibratie is beschreven in bijlage IV. Hieruit blijkt dat aan de vooraf gestelde doelstellingen voor de hydraulische calibratie niet kon worden voldaan. De effecten op de morfologische modellering zijn echter verwaarloosbaar. Bovendien zijn in vergelijking met het oorspronkelijke hydraulische model de calibratieresultaten niet verslechterd. Het gebruik van het nieuwe hydraulische model in plaats van het oude model als basis voor het morfologische model is hiermee gerechtvaardigd.

14 Morfologische gevolgen van rivierverruimende maatregelen

3.3.2 Aanpassing bodemprofielen De bodemprofielen in het huidige model zijn gebaseerd op bodempeilingen uit 1990. Door de uitbreiding van de baseline schematisatie op enkele takken zijn nu op die takken de bodemprofielen overeenkomstig de peilingen uit 1999. Voor de morfologische calibratie is het noodzakelijk om deze bodemprofielen aan te passen, om overal een gelijke startbodem te verkrijgen. Wegens gebrek aan bodempeilingen van de takken in de Biesbosch, kan voor deze takken geen correctie worden uitgevoerd.

...... Figuur 3.3 Nieuwe Merwede Afwijking beginbodem SOBEK ten opzichte van meting voor de Nieuwe Merwede/ 0,75 ) m ( g

in 0,00 ijk w f a

nieuw e afw ijking oorspronkelijke afw ijking -0,75 961 963 965 967 969 971 973 975 977 979 981 kmr

...... Figuur 3.4 Beneden Merwede Afwijking beginbodem SOBEK ten opzichte 0,75 van meting voor de Beneden Merwede. ) m ng (

i 0,00 k j i w af

nieuw e afw ijking oorspronkelijke afw ijking -0,75 961 963 965 967 969 971 973 975 kmr

Voor de andere takken is de correctie uitgevoerd door een sedimentonttrekking of toevoeging te bepalen, gelijk aan het volumeverschil van de bodem tussen 1990 en 1999. Voor elk SOBEK hoofdgeulvak (hoofdgeulvakken zijn weergegeven in figuur 3.2) is dit volumeverschil bepaald door de bodemhoogtegrids (op basis van metingen) uit 1990 en 1999 van elkaar af te trekken en de zo berekende verschillen per gridcel te sommeren over alle cellen binnen een SOBEK-vak. Deze onttrekking of toevoeging kan worden toegepast in de eerste tijdstap van een modelrun, zodat de bodem ook in één tijdstap wordt aangepast. Figuur 3.3 en figuur 3.4 laten respectievelijk voor de Nieuwe en de Beneden Merwede de resultaten van deze aanpassing zien. De SOBEK beginbodem komt goed overeen met de

15 Morfologische gevolgen van rivierverruimende maatregelen

bodemmetingen uit 1990. De afwijking is veel kleiner geworden dan de oorspronkelijke afwijking van de beginbodem (op basis van peilingen uit 1999 was geschematiseerd) in vergelijking met de peiling uit 1990.

3.4 Invoer gebiedsspecifieke data

3.4.1 Invoer gegevens korreldiameters In 2002 is een grootschalig bodemmonster onderzoek gedaan voor vrijwel het gehele NDB [Fugro, 2002]. Op basis hiervan kunnen allerhande korrelgroottes bepaald worden. Er is gekozen om in SOBEK 5 km gemiddelde korrelgroottes in te voeren. De ruwe gegevens vertonen vaak abrupte schommelingen per plaats. Deze uitschieters kunnen bijdragen aan de instabiliteit van het model. Ook fysisch gezien kan een opeenvolging van grove en fijne trajecten leiden tot irreële verhoging of verlaging van de transportcapaciteit. Immers, in werkelijkheid zal er veel meer sprake zijn van geleidelijke overgangen.

...... Figuur 3.5 Waal - Boven Merwede - Beneden Merwede Mediane korreldiameter langs Waal – Boven Merwede – Beneden Merwede. 2,5 ruw e data

2 5 km gemiddeld )

m 1,5 m ( 0

D5 1

0,5

0 926 931 936 941 946 951 956 961 966 971 976 kmr

...... Figuur 3.6 Waal - Boven Merwede - Nieuwe Merwede Mediane korreldiameter langs Waal – Boven Merwede – Nieuwe Merwede. 2,5 ruw e data 5 km gemiddeld 2

) 1,5 m m 50 (

D 1

0,5

0 926 931 936 941 946 951 956 961 966 971 976 981 kmr

Figuur 3.5 en figuur 3.6 geven een overzicht van de ruwe en de 5 km gemiddelde waarden voor respectievelijk de takken Waal-Boven Merwede-Beneden Merwede en Waal – Boven Merwede - Nieuwe

16 Morfologische gevolgen van rivierverruimende maatregelen

Merwede. Ook al ligt de focus van deze studie op de Waal en de Merwede takken, de korrelgegevens van de overige takken zullen ook worden ingevoerd in het model om zodoende een goed ‘basismodel’ te verkrijgen voor het gehele NDB (bijlage V).

3.4.2 Invoer gegevens baggerwerk Voor een deel van het NDB zijn baggergegevens over de periode 1990- 2000 beschikbaar gesteld [RWS-DZH, 2002]. Deze gegevens, waaronder de baggercijfers voor de Merwede takken en de Waal, zijn ingevoerd in SOBEK. Voor de vier onderzoekstakken zijn de baggerhoeveelheden weergegeven in bijlage VI. De baggerdata voor de overige takken zijn niet in dit rapport opgenomen. Voor de onderzoekstakken geldt dat de data vrij gedetailleerd zijn ingevoerd, namelijk per kilometer. Voor de overige takken zijn de gegevens doorgaans per riviervak opgegeven.

3.4.3 Randvoorwaarden Op alle randen die het NDB model kent, zijn zowel hydraulische als morfologische voorwaarden opgelegd. Het model kent drie bovenstroomse randen, namelijk Lith voor de Maas, Tiel voor de Waal en Hagestein voor de Lek. Verder zijn er drie benedenstroomse randen, namelijk twee in het Haringvliet (er zijn hier namelijk twee geulen) en de monding van de Nieuwe Waterweg, ook wel de Maasmonding genoemd. Alle andere randen in het model betreffen doodlopende takken en/of havens. Aan de randen van deze doodlopende takken en/of havens zijn als morfologische randvoorwaarde een sedimenttransport van 0 en als hydraulische voorwaarde een debiet van 0 opgegeven.

In tabel 3.1 zijn de randvoorwaarden gegeven, die tijdens de calibratie voor de overige randen zijn toegepast.

...... rand Hydraulische Morfologische Tabel 3.1 Hydraulische en morfologische randvoorwaarde randvoorwaarde randvoorwaarden SOBEK-model. Maasmond getij, periode 12u25m bodemhoogte: -25,5 Haringvliet getij, periode 12u25m bodemhoogte: -4,73 Haringvliet-20 getij, periode 12u25m bodemhoogte: -4,45 Maas (Lith) daggemiddelde afvoer bodemhoogte: -5,12 Waal (Tiel) daggemiddelde afvoer bodemhoogte: -1,66 Lek (Hagestein) daggemiddelde afvoer bodemhoogte: -4,93

Het betreft hier dus aan de bovenstroomse randen daggemiddelde debieten [DONAR2] gebaseerd op de werkelijk gemeten debieten in de periode 1990 – 2000 (figuur 3.7). Meer gedetailleerde afvoeren zijn niet voor alle drie de takken voorhanden. Bovendien zullen bij een tijdschaal van 10 jaar de effecten op de morfologie bij het gebruik van meer gedetailleerde afvoergegevens verwaarloosbaar klein zijn. De benedenstroomse hydraulische randvoorwaarden zijn gebaseerd op een

2 DONAR (Data Opslag Natte Rijkswaterstaat) is de centrale database van Rijkswaterstaat, waarin al haar fysische, chemische, biologische en morfologische gegevens worden opgeslagen.

17 Morfologische gevolgen van rivierverruimende maatregelen

morfologisch bepalend getij (figuur 3.8). Het afvoerregime door de Haringvlietsluizen is volgens het afvoerregime LPH843. Deze benedenstroomse hydraulische randvoorwaarden en alle morfologische randvoorwaarden zijn overgenomen uit het model van Wang [2001].

...... Figuur 3.7 8000 Hydraulische randvoorwaarde voor de ) s Waal. / 7000 3 m

( 6000 el i

T 5000 oer v 4000 de af l 3000 de d 2000 emi

agg 1000 D 0 0 1 2 3 4 5 6 7 8 9 0 199 199 199 199 199 199 199 199 199 199 200 Tijd

...... Figuur 3.8 Hydraulische randvoorwaarde voor de 1,2 Maasmond. ) . P

. 0,8 A . N

m 0,4

mond ( 0 as a 0:00 2:00 4:00 6:00 8:00 10:00 12:00 14:00 16:00 18:00 20:00 22:00 0:00

nd M -0,4 ta s r e t

a -0,8 W

-1,2 Tijd

3.5 Modelinstellingen

Tijdens de runs met het model is de optie ‘river’ ingeschakeld, aangezien de optie ‘estuarine’ nog niet goed blijkt te werken, zoals in het onderzoek van Wang [2001] naar voren kwam. De optie ‘salt intrusion’ is niet toegepast, aangezien dit onderzoek zich richt op het bovenstroomse deel van het NDB. Zoutindringing beperkt zich slechts tot de Nieuwe Waterweg, de Nieuwe Maas en de Oude Maas [Houwing en Blom, 2001]. De instellingen betreffende de numerieke

3 Huidig beheersprogramma Haringvlietsluizen. Dit beheer biedt veiligheid tegen overstroming, handhaaft de zoetwaterhuishouding in de zuidrand en minimaliseert de kans op verzilting via de noordrand. De afvoerverdeling tussen de noord- en zuidrand van het NDB wordt geregeld door de opening van de Haringvlietsluizen.

18 Morfologische gevolgen van rivierverruimende maatregelen

parameters worden overgenomen uit het model van Wang. Veel van deze instellingen zorgen voor de stabiliteit van een model run.

Verder worden voor alle takken de volgende instellingen gebruikt: • transportformule: Engelund & Hansen [Wang, 2001]; • calibratiefactor transportformule: 1; • sedimentverdeling splitsingspunten volgens de machtsrelatie Wang (paragraaf 3.2) met k = 2 [Wang, 2001]. Deze instellingen zijn tijdens de eerste runs van de calibratie gehanteerd, maar zijn later aangepast in verband met de afregeling van het model (hoofdstuk 4).

De rekentijdstap wordt op 1 uur gezet. Een lagere waarde zou leiden tot een verhoging van de rekentijd. Een hogere waarde kan gezien de periode van de benedenstroomse hydraulische randvoorwaarden niet gebruikt worden om een realistische berekening te maken. De rekenperiode tijdens de calibratie is 10 jaar: van 1990 tot en met 1999.

3.6 Initiële Run

Met de hierboven omschreven modelinstellingen is een initiële run uitgevoerd, om zodoende het model te controleren op stabiliteit. Het bleek niet mogelijk te zijn om het model voor de gehele 10-jarige periode te laten draaien. De oorzaak hiervan is de schematisering van enkele havens en plassen in SOBEK, namelijk door middel van een extra stroomvoerend oppervlak in de uiterwaard (achter de zomerdijk). Wanneer er water door de uiterwaard stroomt beschouwt SOBEK dit als een debietonttrekking uit de hoofdgeul op de betreffende locatie. Hierdoor is het verloop van de stroomsnelheid in de hoofdgeul niet continu, wat een enorme sedimentatie tot gevolg heeft. De sedimentatie wordt versterkt doordat er alleen binnen de sedimentvoerende breedte morfologische veranderingen kunnen plaatsvinden en deze breedte is per definitie gelijk aan of kleiner dan de breedte van de hoofdgeul. Beter zou zijn om dit gelijk te stellen aan de stroomvoerende breedte, zodat de sedimentatie verspreid zou worden over een veel groter oppervlak. Dit is echter niet mogelijk in SOBEK. Er is nu voor gekozen om de aanwezige plassen en havens te schematiseren door middel van extra bergend oppervlak in plaats van stroomvoerend oppervlak. Hierdoor vindt er geen debietonttrekking plaats en dus blijft de enorme sedimentatie achterwege. Deze aanpassing lijkt ook realistischer aangezien dit soort bekkens vooral invloed heeft op de bergende breedte. Na de aanpassing van de probleemprofielen is de test run geslaagd.

19 Morfologische gevolgen van rivierverruimende maatregelen

20 Morfologische gevolgen van rivierverruimende maatregelen

4 Calibratie

......

4.1 Inleiding

In dit hoofdstuk wordt de morfologische calibratie van het model beschreven. Nu het model is aangevuld met gebiedsspecifieke data en vernieuwde schematisaties, is een betere calibratie mogelijk dan voor het model van Wang [Wang, 2001]. Na een geslaagde calibratie wordt het model toegepast in het kader van ‘ruimte voor de rivier’ (hoofdstuk 5). Uit de doelstelling volgt dat alleen de vier onderzoekstakken worden gecalibreerd. Verder volgt hieruit dat er alleen gecalibreerd wordt op grootschalige morfologische ontwikkeling, dat wil zeggen de totale bodemvolume verandering per tak gedurende de calibratieperiode (1990 – 2000). De data aan de hand waarvan het model wordt afgeregeld zijn hier calibratiedata genoemd. In paragraaf 4.2 wordt de aanpak van de calibratie beschreven. Hierin komen de calibratiedata naar voren, maar ook de calibratievariabelen; dat zijn de keuzemogelijkheden of waarden van parameters in SOBEK- MOR die aangepast worden om zodoende het model af te regelen. In de vier paragrafen die daarop volgen zijn de effecten onderzocht van wijzigingen in de calibratievariabelen. In paragraaf 4.7 zijn de uiteindelijke keuzes voor de calibratievariabelen gegeven.

4.2 Calibratiestrategie

4.2.1 Calibratiedata In bijlage VII zijn figuren opgenomen waarin de bodemliggingen in 1990 en 2000 voor de onderzoekstakken zijn gegeven. Het model is afgeregeld op basis van de grootschalige morfologische veranderingen gedurende 1990 – 2000. De waarden hiervan staan vermeld in tabel 4.1. Het betreft alleen de morfologische veranderingen binnen de hoofdgeul (main channel), want SOBEK berekent ook alleen hier veranderingen. De totale bodemveranderingen volgens de metingen worden vergeleken met de resultaten van SOBEK. Hiervoor zijn de veranderingen van de dwarsprofielen gebruikt die SOBEK als output genereert. Van elk dwarsprofiel wordt de verandering in dwarsdoorsnede vermenigvuldigd met de bijbehorende SOBEK vaklengte en vervolgens wordt gesommeerd over alle dwarsprofielen per tak om de totale bodemverandering per tak te bepalen.

...... 3 Tabel 4.1 Tak Totale bodemverandering [1000 m ] Totale volume verandering bodem Waal -631 hoofdgeul per onderzoekstak tussen 1990 en 2000. Boven Merwede -652 Nieuwe Merwede -1071 Beneden Merwede -30

21 Morfologische gevolgen van rivierverruimende maatregelen

Wanneer de afwijking in modelresultaten kleiner is dan de totale absolute onnauwkeurigheid in de gebruikte data, is de calibratie geslaagd. Het probleem hierbij is dat de onnauwkeurigheid in de data niet bekend is. Uiteraard is niet alleen de onnauwkeurigheid van de bodempeilingen van belang, maar ook de onnauwkeurigheid in baggerdata. De baggercijfers zijn namelijk mede bepalend voor de morfologische ontwikkeling die door het model wordt voorspeld. In figuur 4.1 is een vergelijking gemaakt tussen de geobserveerde morfologische verandering en de baggerhoeveelheden. Het verschil hiertussen is de natuurlijke bodemverandering. Uit de figuur volgt dat de baggerwerkzaamheden een zeer significante invloed op de grootschalige morfologie hebben. Door de relatief grote onnauwkeurigheid die vaak met baggercijfers samengaat en de betrekkelijk grote hoeveelheden baggerwerk, zal de afwijking van de voorspelde grootschalige morfologische ontwikkeling ten opzichte van de werkelijke ontwikkeling gerelateerd worden aan de onnauwkeurigheid in baggerdata. De calibratie-eis luidt dan dat deze afwijking kleiner moet zijn dan de onnauwkeurigheid in de baggerdata plus een marge om te corrigeren voor overige datafouten. Echter deze onnauwkeurigheid van de gebruikte baggerdata is niet bekend. Volstaan moet worden met gegevens uit de literatuur en op basis van ‘expert judgement’. Tijdens de calibratie van het Rijntakken model [Jesse en Kroekenstoel, 2001] zijn onnauwkeurigheden gebruikt variërend tussen de 15 en zelfs meer dan 100%. Andere bronnen melden onnauwkeurigheden rond de 20%. Voor de calibratie wordt hier gestreefd om de afwijking in grootschalige morfologische ontwikkeling kleiner dan 30% van de betreffende baggervolumes te laten zijn.

...... Figuur 4.1 Werkelijke, natuurlijke en kunstmatige bodemverandering tussen 1990 en 2000 De invloed van baggeren op de bodemverandering in de Merwede taken. 2 ] m3 n

e 1 o j l i m [ k a t

0 per g n i -1 der n a r e -2 mev u l o V

-3 Boven Merw ede Beneden Merw ede Nieuw e Merw ede

Werkelijke bodemverandering Verandering door baggeren Natuurlijke verandering

4.2.2 Calibratievariabelen Calibratievariabelen zijn variabelen die worden aangepast om zodoende het model af te regelen. In dit onderzoek zijn de keuzemogelijkheden of waarden van parameters in SOBEK-MOR hiervoor gebruikt. Voor het afregelen van het SOBEK-MOR model zijn verschillende variabelen beschikbaar. Bij het aanpassen van een variabele verdient het de voorkeur om de waarde ervan binnen de - door de literatuur gestelde

22 Morfologische gevolgen van rivierverruimende maatregelen

of op basis van metingen bepaalde - grenzen te houden. Dit beperkt niet alleen het aantal variatiemogelijkheden maar ook het aantal calibratievariabelen zelf, aangezien sommige parameters vastgestelde waarden hebben. De volgende calibratievariabelen zijn gebruikt:

• Transportformule In SOBEK zijn een aantal transportformuleringen voorhanden. In de voorstudie is reeds gebruik gemaakt van de formule van Engelund & Hansen (1967). Het sedimenttransport per meter breedte wordt als volgt berekend: 0.05u5 s = m (vergelijking 4.1) 32 (1 −∆ε ) gC d D50 met:

um gemiddelde stroomsnelheid in de hoofdgeul [m/s] C Chézy coëfficiënt [m½/s] ε porositeit bodemmateriaal [-]

D50 korreldiameter waarbij 50% van het materiaal fijner is [m]

∆d relatieve dichtheid sediment [-]

Deze laatste term is als volgt gedefinieerd: ρ − ρ ∆= s (vergelijking 4.2) d ρ met: 3 ρs dichtheid van sediment [kg/m ] ρ dichtheid van water [kg/m3]

Deze formule berekent het totale transport (exclusief washload). De formule is geverifieerd voor D50 korreldiameters die liggen tussen de

0.19 en 0.93 mm. De D50 die op de Merwede takken worden aangetroffen liggen doorgaans binnen dit bereik. Alleen op de Waal komen diameters boven de 0.93 mm voor. Een totaal transportformule lijkt het meest geschikt voor de Merwede takken. Van Ledden concludeerde in zijn onderzoek naar zand-slib-segregatie op de Nieuwe Merwede, dat bodemtransport van zand een ondergeschikte rol blijkt te spelen op de Nieuwe Merwede [Van Ledden, 1999].

Naast bovenstaande totaaltransport formule wordt in deze studie ook de formule van Meyer-Peter-Müller onderzocht. In tegenstelling tot de vorige formule berekent deze formule alleen bodemtransporten. Deze formule wordt dan ook vaak toegepast op riviertakken met grof bodemmateriaal. Meyer-Peter-Müller wordt hier onderzocht omdat deze formule (in aangepaste vorm) in de calibratie van het Rijntakken model [Jesse en Kroekenstoel, 2001] de beste resultaten gaf.

De formule van Meyer-Peter-Müller luidt:

3/2 3 ⎛⎞2 8 gD∆dm u s =⎜m −0.047⎟ (vergelijking 4.3) (1 −ε ) ⎜⎟3 ⎝⎠CC90 ∆dmD

23 Morfologische gevolgen van rivierverruimende maatregelen

met:

C90 Chézy coëfficiënt gerelateerd aan de D90 van het het bodemmateriaal [m½/s]

D90 korreldiameter waarbij 90% van het materiaal fijner is [m]

Dm gemiddelde korreldiameter [m]

Bij het onderzoek naar de twee transportformules is er voor gekozen om telkens dezelfde formule voor alle takken op te geven. Dit voorkomt plotselinge en onrealistische variaties in sedimenttransport tussen twee takken.

• Calibratiefactor van de transportformule

Deze calibratiefactor (cf,vergelijking 3.2) is een vermenigvuldigingsfactor, waarmee het eigenlijke sedimenttransport wordt berekend. Het eigenlijke sedimenttransport is het transport volgens de transportformule vermenigvuldigd met deze factor. Standaard is de factor ingesteld op 1. Een verandering van deze factor op een bepaalde tak heeft niet alleen invloed op de tak zelf, maar ook op de benedenstroomse takken. Zo zal een verlaging van de factor op een bepaalde tak leiden tot een lager sedimenttransport richting de takken benedenstrooms ervan. Deze factor is in SOBEK ingebouwd voor calibratiedoeleinden, maar kan ook gebruikt worden om de rekentijd te verkorten door de factor groter dan één te kiezen. In dit onderzoek is de factor alleen gebruikt als calibratiemiddel.

• Korreldiameter bodemmateriaal De in het basismodel ingevoerde korreldiameters betreffen 5 km- gemiddelde waarden uit bodembemonsteringen in 2002. Uiteraard zit er in het gebruik van bepaalde korreldiameters een grote onzekerheid. Een bodembemonstering betreft namelijk altijd een momentopname, korreldiameters zullen in de loop der tijd veranderen. Ook uit de metingen blijkt dat er een behoorlijke marge zit in het verloop van de korrelgrootte in rivierrichting (figuur 3.5 en figuur 3.6). Wanneer bepaalde trajecten met de eerder genoemde calibratievariabelen niet goed afgeregeld kunnen worden, kan er voor gekozen worden om veranderingen aan te brengen in de korreldiameter. Een verandering van korreldiameters moet echter wel redelijk binnen de uit de metingen naar voren gekomen marges liggen, zodat de aanpassingen niet onrealistisch zijn.

• Bodemruwheid De bodemruwheid is reeds bepaald gedurende de hydraulische calibratie van het model. De variabele gaat gepaard met een grote onzekerheid. Uit vergelijking 4.1 en 4.3 volgt dat de bodemruwheid van invloed is op het sedimenttransport en dus op de morfologische ontwikkeling. Tijdens de morfologische calibratie is de ruwheid ook gebruikt als calibratievariabele. De consequentie van het aanpassen van de bodemruwheid is dat dit ook het hydraulische model beïnvloedt en dat daarom het hydraulische model gecontroleerd moet worden.

24 Morfologische gevolgen van rivierverruimende maatregelen

4.3 Effect aanpassing transportformule

Met de twee voorgestelde transportformules zijn afzonderlijke runs uitgevoerd. De overige calibratievariabelen zijn hierin nog niet aangepast en zijn dus allen op de standaardwaarden ingesteld. Voor beide formules worden de resultaten in de volgende subparagrafen behandeld.

4.3.1 Engelund & Hansen De totale bodemveranderingen per tak zijn in figuur 4.2 weergegeven. Tabel 4.2 geeft een overzicht van de afwijkingen in relatie tot de baggervolumes. De resultaten van deze run tonen over het algemeen een behoorlijke sedimentatie, terwijl uit de metingen erosie volgt. Alleen op de Waal wordt erosie berekend, waar dat ook volgt uit de metingen. Voor de Merwede takken zijn de resultaten slecht. De voorspellingen wijken veel af van de metingen.

...... Figuur 4.2 Totale bodemverandering 1990 - 2000 Vergelijking totale bodemveranderingen per tak gedurende de periode 1990-2000 3000 van de model runs met E&H en MPM met E&H de waargenomen verandering. 2000 MPM Meting 1000 m3] 1000 ng [ i er

and 0 r e

demv -1000 o B

-2000 Waal Boven Nieuw e Beneden Merw ede Merw ede Mer w ede

...... Tabel 4.2 Tak Afwijking meting Relatief t.o.v. Afwijkingen van de waargenomen [1000 m3] baggervolume bodemveranderingen voor Engelund & Hansen, zowel absoluut als in relatie tot de Waal -577 35% baggervolumes. Boven Merwede 3633 143% Nieuwe Merwede 3045 432% Beneden Merwede 2515 541%

4.3.2 Meyer-Peter-Müller De berekening met Meyer-Peter-Müller laat voor de Merwede takken op het eerste gezicht veel betere resultaten zien dan de vorige berekening (figuur 4.2 en tabel 4.3). De afwijkingen op de Boven en de Beneden Merwede zijn slechts 20% van de gebaggerde hoeveelheden. Dit is ruim binnen de vooraf aangenomen foutmarge van de baggercijfers. Ook voor de Nieuwe Merwede is de voorspelling verbeterd, zij het nog niet voldoende. Alleen voor de Waal is de voorspelling juist verslechterd; het model voorspelt nu sedimentatie in plaats van de geobserveerde erosie.

25 Morfologische gevolgen van rivierverruimende maatregelen

...... Tak Afwijking meting Relatief tov Tabel 4.3 3 Afwijkingen van de [1000 m ] baggervolume waargenomen Waal 1574 94% bodemveranderingen voor MPM, zowel absoluut als in Boven Merwede 543 21% relatie tot de baggervolumes. Nieuwe Merwede 741 105% Beneden Merwede -100 22%

4.3.3 Discussie Gelet op de voorspelde morfologische ontwikkeling die SOBEK met de twee formules liet zien, lijkt MPM in de eerste plaats de beste optie. Deze formule laat namelijk op de Merwede takken erosie zien, wat ook uit de metingen naar voren is gekomen. E&H daarentegen voorspelde een behoorlijke sedimentatie. Opmerkelijk, aangezien transportformules, ondanks de verschillen ertussen, globaal hetzelfde beeld zouden moeten voorspellen. Nader onderzoek wijst uit dat beide transportformules in principe wel hetzelfde beeld laten zien, namelijk sedimentatie, maar doordat de hoeveelheid sedimentatie die MPM berekent veel kleiner is dan de baggerhoeveelheden, is het netto resultaat toch erosie. Om dit beeld te controleren worden de berekende sedimenttransporten met de twee formules vergeleken met waarden uit de literatuur (figuur 4.3). De gebruikte waarden zijn afkomstig uit de sedimentbalansen voor de Rijntakken [Ten Brinke, et al., 2001] en voor het NDB [van Dreumel, 1995]. Vergelijkingen worden gemaakt met jaargemiddelde zandtransporten op de Boven Merwede (kmr 961), de Beneden Merwede (kmr 961 en kmr 976) en de Nieuwe Merwede (kmr 961 en kmr 980).

...... Figuur 4.3 Vergelijking sedimenttransporten rond splitsingspunt Merwede (kmr 961.3) en benedenstrooms van de Nieuwe Merwede (kmr 980) en de Beneden Merwede (kmr 975) voor E&H (boven) en MPM (beneden).

Splitsingspunt Merweden Benedenstrooms Nieuwe en Beneden Merwede ]

1,0E+06 j 5,0E+05 3/ ] j m 3/

8,0E+05 [ 4,0E+05 t m [ t 6,0E+05 por 3,0E+05 por ans s n 4,0E+05 tr 2,0E+05 a r t

2,0E+05 ment 1,0E+05 ment edi di S

e 0,0E+00 0,0E+00 S Nieuw e Merw ede Beneden Merw ede Boven Merw ede Nieuw e Merw ede Beneden Merw ede

Metingen Initiële run E&H Metingen Initiële run E&H

Splitsingspunt Merweden Benedenstrooms Nieuwe en Beneden Merwede 1,0E+06 ]

j 5,0E+05 ] 3/ j 3/ 8,0E+05 m

[ 4,0E+05 m t [ t

6,0E+05 por 3,0E+05 s por n s a r

4,0E+05 t an 2,0E+05 r t

2,0E+05 ment 1,0E+05 ment edi S edi 0,0E+00 0,0E+00 S Nieuw e Merw ede Beneden Merw ede Boven Merw ede Nieuw e Merw ede Beneden Merw ede

Metingen Initiële run MPM Metingen Initiële run MPM

26 Morfologische gevolgen van rivierverruimende maatregelen

De transporten die met MPM worden berekend zijn veel kleiner dan met E&H, wat tot gevolg heeft dat de gradiënt in transport ook veel kleiner is en dus ook de bodemveranderingen. Uit de figuur volgt dat met E&H veel betere sedimenttransporten worden berekend. Alleen benedenstrooms op de Nieuwe en Beneden Merwede kloppen de berekende transporten niet meer, deze zijn veel kleiner dan in werkelijkheid. Dit verklaart ook de te grote sedimentatie op beide takken.

4.3.4 Keuze Uit bovenstaande kan geconcludeerd worden dat de formule van Engelund & Hansen de voorkeur geniet boven die van Meyer-Peter- Müller. Verder is gebleken dat het bodemtransport een ondergeschikte rol blijkt te spelen op de Nieuwe Merwede [van Ledden, 2001] en het gebruik van een totaaltransport formule dus reëler is. Het gebruik van MPM zou door het veel te laag ingeschatte sedimenttransport zeer onrealistisch zijn.

4.4 Effect aanpassingen calibratiefactor

In deze paragraaf is de invloed van de calibratiefactor van de transportformule op de morfologische ontwikkeling onderzocht. Op basis van de voorgaande paragraaf is gerekend met E&H. Er zijn drie runs opgesteld waarin verschillende waarden aan deze factor toegekend. De aanpassingen zijn er op gericht de in de voorgaande paragraaf genoemde probleempunten met Engelund & Hansen op te lossen.

4.4.1 Modelaanpassingen Drie runs zijn getest met de volgende instellingen: • Run 1: op de Merwede takken en op de Waal is de calibratiefactor van de transportformule verlaagd tot 0.8. Dit zou de grote veranderingen in bodemligging kunnen verminderen. Een lager transport leidt immers tot een kleinere gradiënt in de transporten, waardoor de bodemveranderingen kleiner zullen uitpakken. • Run 2: in een tweede run is geprobeerd om de toevoer van sediment vanuit de Waal en Boven Merwede te verkleinen door de alleen hier de calibratiefactor van de transportformule te verkleinen naar 0.8. Hiermee wordt getracht de grote sedimentatie op de Beneden en Nieuwe Merwede op te lossen. • Run 3: een andere mogelijkheid om de grote sedimentatie op de Nieuwe en Beneden Merwede op te lossen is middels een run waarin de calibratiefactor van de takken benedenstrooms van deze takken is verhoogd naar 1.2. Verdere instellingen waren: calibratiefactor op de Waal en Boven Merwede 0.7 en op de Beneden en Nieuwe Merwede 0.9. Op deze manier wordt de afvoer van sediment bevorderd en de toevoer juist geremd.

4.4.2 Resultaten De resultaten voor grootschalige morfologische ontwikkeling (figuur 4.4) tonen aan dat de aanpassingen bijdragen aan een betere voorspelling hiervan voor de Merwede takken. Voor de Merwede

27 Morfologische gevolgen van rivierverruimende maatregelen

takken geeft run 3 de beste resultaten. Dit gaat echter ten koste van de resultaten op de Waal. Run 1 laat voor alle takken een verbetering zien.

...... Figuur 4.4 Totale bodemverandering 1990 - 2000 Totale bodemveranderingen per tak gedurende de periode 1990-2000 volgens 4000 de meting, de initiële run met E&H en de ] drie runs met aangepaste 3 3000 calibratiefactoren.

00 m 2000 Initiële run E&H

g [10 Run 1

n 1000 i

0 Run 2 ander er

v Run 3 -1000 m

de Meting o -2000 B

-3000 Waal Boven Nieuw e Beneden Merw ede Merw ede Mer w ede

Het blijkt onmogelijk, ondanks de verbeterde voorspellingen van run 3, om een goed calibratieresultaat te verkrijgen. Zelfs de aard van de grootschalige morfologische ontwikkeling voor de Merwede takken, erosie dan wel sedimentatie, wordt nog niet goed voorspeld door welke run dan ook. De afwijkingen van de voorspelde grootschalige morfologische ontwikkeling ten opzichte van de baggervolumes laten ook zien dat nog niet aan de vooraf gestelde eisen voldaan kan worden (figuur 4.5).

...... Figuur 4.5 Afw ijking grootschalige morfologie Afwijkingen in voorspelling grootschalige morfologie gerelateerd aan de 600% baggervolumes voor de drie runs alsmede voor de initiële run met E&H. 500% Initiële run 400% Run 1 g n i Run 2 k j

i 300%

w Run 3

Af 200% Eis

100%

0% Waal Boven Nieuwe Beneden Merwede Merwede Merwede

4.4.3 Discussie De aanpassingen in de calibratiefactoren leveren verbeteringen op in de voorspelling van grootschalige morfologische ontwikkeling. Echter, met de voorgestelde aanpassingen kan nog niet aan de calibratie-eis worden voldaan. Verdere aanpassing van deze factoren lijkt niet zinvol; de problemen op bepaalde trajecten zijn zodanig dat er in andere oplossingsrichtingen gezocht moet worden. Bovendien zou een calibratiefactor van de transportformule die erg veel van 1 afwijkt niet meer reëel zijn.

28 Morfologische gevolgen van rivierverruimende maatregelen

4.4.4 Keuze Gekozen is om verder te gaan met de calibratiefactoren van de transportformules op alle onderzoekstakken ingesteld op 0.8 (volgens run 1). Deze simpele aanpassing geeft naast de verbeterde morfologische voorspelling namelijk ook een betere voorspelling van de sedimenttransporten nabij het splitsingpunt ten opzichte van de initiële run.

4.5 Effect aanpassing korreldiameter

Tot nu toe zijn alleen calibratievariabelen gebruikt die betrekking hebben op gehele takken. Transportformule en calibratiefactor zijn per tak constant. Korreldiameters kunnen binnen een tak variëren. Bij gebruikmaking van deze calibratievariabele is het dus van belang te weten waar de grootste probleemtrajecten zich bevinden.

...... Figuur 4.6 Waal en Boven Merwede Gemeten bodemverandering (blauw)versus 2,00 de door SOBEK berekende (paars) bodemverandering tussen 1990 en 2000 ) m voor run 1 uit paragraaf 4.4 (cf = 0.8). (

l i h c s r 0,00 e v odem b

-2,00 913 923 933 943 953 963 kmr

Nie uw e M e r w e de 2,00 ) m (

l i h c s r 0,00 e v em bod

-2,00 961 965 969 973 977 981 kmr

Beneden Merwede 2,00 ) m (

l i h c s r 0,00 e v em bod

-2,00 961 964 967 970 973 976 kmr

29 Morfologische gevolgen van rivierverruimende maatregelen

Figuur 4.6 toont de berekende bodemveranderingen tussen 1990 en 2000 in vergelijking met de gemeten veranderingen voor de onderzoekstakken. De volgende opmerkingen kunnen hierbij worden geplaatst: • Waal: Op de Waal zijn er bovenstrooms locaties waar teveel erosie wordt berekend (rond kmr 918, 927 en 930). Verder benedenstrooms vindt er juist teveel sedimentatie plaats (kmr 938). • Boven Merwede: op de gehele Boven Merwede berekent het model sedimentatie in plaats van erosie. • Nieuwe Merwede: op deze tak valt vooral de enorme sedimentatie vanaf kmr 972 op. Verder wordt aan het begin van de tak ook veel sedimentatie berekend. • Beneden Merwede: een behoorlijke sedimentatie in het middendeel (rond kmr 968) valt op en ook de sedimentatie vanaf kmr 970.

4.5.1 Modelaanpassingen Een fout in morfologische voorspelling op een bepaald traject wordt veroorzaakt door een fout in de transportgradiënt (vergelijking 3.1) en dus ook in het transport zelf. Met behulp van vergelijking 3.1 en de gebruikte transportformule (E&H, vergelijking 4.1) zijn ruwe schattingen gemaakt van de benodigde aanpassingen in de gradiënt van de korreldiameter voor een goede morfologische voorspelling. Vervolgens zijn er enkele runs gedraaid waarin verschillende aanpassingen zijn onderzocht.

De aanpassingen in gradiënt zijn gedaan op bovengenoemde probleemtrajecten. Hierbij is er voor gekozen om de korreldiameters rondom het splitsingspunt van de Merweden onaangetast te laten, om de sedimentverdeling niet negatief te beïnvloeden. Omdat vanuit deze locaties met onveranderde korreldiameters is gewerkt, is het wel zo dat op trajecten waar geen aanpassingen nodig waren, wel andere korreldiameters zijn toegepast, echter de gradiënt op die locaties is niet veranderd.

4.5.2 Resultaten Niet alleen op basis van de resultaten van de runs is een beste oplossing gekozen. Ook de mate waarin de aanpassingen passen binnen de marges zoals die uit de metingen naar voren zijn gekomen tellen mee in deze beoordeling. Op basis daarvan zijn de aanpassingen in korreldiameter die getoond zijn in bijlage VIII gehanteerd voor het vervolgonderzoek. De aanpassingen in korrelgrootte liggen binnen de reële marges (visueel, met behulp van de figuren uit bijlage VIII) die volgen uit de metingen. Alleen voor de Nieuwe Merwede lijken de veranderingen enigszins buiten deze marges te vallen.

De resultaten met betrekking tot de grootschalige morfologie (figuur 4.7) laten zien dat voor drie van de vier takken aan de eis is voldaan. Hier is de afwijking namelijk kleiner dan 30% van de baggervolumes. Alleen voor de Beneden Merwede is de afwijking nog te groot. Overigens is wel een behoorlijke verbetering ten opzichte van de basisrun behaald.

30 Morfologische gevolgen van rivierverruimende maatregelen

...... Figuur 4.7 Afw ijking grootschalige morfologie Resultaten calibratie na aanpassing korreldiameter. 400% 350% 300%

g 250% in

ijk 200% w f 150% A 100% 50% 0% Waal Boven Merw ede Nieuw e Merw ede Beneden Merw ede

Uitgangssituatie (run 1 par. 4.4) Aanpassing korreldiameters Eis

4.5.3 Discussie Op alle takken heeft de aanpassing in korreldiameter geleid tot een behoorlijke verbetering in de resultaten. De aangepaste korreldiameters vallen grotendeels binnen de uit de metingen naar voren gekomen marges. Alleen op de laatste 10 km van de Nieuwe Merwede zijn de aanpassingen behoorlijk. Deze grove aanpassing is hier nodig om de sedimentatie tegen te gaan. Het model is blijkbaar niet goed in staat hier de werkelijke situatie na te bootsen. Wellicht speelt de aanwezigheid van slib een rol. Uit een recent onderzoek blijkt dat op de laatste 10 km van de Nieuwe Merwede slib voorkomt [Koomans, 2003]. De slibpercentages op dit traject kunnen op lopen tot plaatselijk 60 – 70 %. Dit traject komt vrijwel overeen met het probleem traject, waar nu de forse aanpassing in korreldiameter is toegepast. In paragraaf 4.7.3 worden enkele mogelijke verklaringen gegeven.

4.5.4 Keuze Voor de Waal en Boven Merwede is gekozen voor de aanpassing in korreldiameter, zoals weergegeven in bijlage VIII. De keuzen voor de Nieuwe en Beneden Merwede hangen af van de resultaten die met de aanpassingen in Chézy waarden zijn behaald en worden derhalve in paragraaf 4.6.4 behandeld.

4.6 Effect aanpassing Chézy waarden

Uit vergelijking 4.1 volgt dat de Chézy-coëfficiënt ook invloed heeft op het sedimenttransport (S ~ 1/C2). Een verandering van deze waarde heeft echter niet alleen gevolgen voor de morfologie, maar ook voor de waterbeweging. Daarom is na elke aanpassing van Chézy waarden gecontroleerd of het hydraulische model niet negatief beïnvloed is. Hiervoor is dezelfde hydraulische run gedraaid als tijdens de hydraulische calibratie en zijn de uitkomsten vergeleken met de uitkomsten van de hydraulische calibratie.

4.6.1 Modelaanpassingen Als uitgangspunt zijn weer de modelinstellingen genomen van run 1 uit

paragraaf 4.4 (E&H, cf = 0.8 op onderzoekstakken en de niet aangepaste korreldiameters). Zodoende kan het effect van

31 Morfologische gevolgen van rivierverruimende maatregelen

aanpassingen in Chézy waarden goed worden vergeleken met het effect van aanpassingen in korreldiameter. Alleen op de Beneden en Nieuwe Merwede zijn aanpassingen in Chézy waarden onderzocht; op de Beneden Merwede omdat hier met de aanpassingen in korreldiameter geen goede resultaten gehaald zijn en op de Nieuwe Merwede omdat hier een aanzienlijke korrelaanpassing nodig is voor een goed resultaat. Aanpassing in de Chézy- coëfficiënt kunnen de benodigde aanpassing in korreldiameter verkleinen. In het model zijn op beide takken slechts op twee locaties Chézy waarden opgegeven, op het beginpunt en eindpunt. Bovendien zijn de waarden aan het eind van de takken hoger dan die aan het begin (tabel 4.4). Dit zou ook een oorzaak kunnen zijn van de sterke afname van sedimenttransport in benedenstroomse richting op beide takken, aangezien het sedimenttransport afneemt bij hogere Chézy waarden. Op beide takken zijn de Chézy waarden zodanig aangepast dat de waarden ervan over de gehele takken constant is. De constante waarden zijn gelijk aan de waarden die in het oorspronkelijke model aan het begin van de takken zijn opgegeven.

...... Beneden Nieuwe Tabel 4.4 Chézy waarden in het oorspronkelijke Merwede Merwede model. waterdiepte [m] begin eind begin eind 0 45 50 50 60 1 45 50 55 60 2 50 55 60 60 3 55 55 60 60 4 55 55 60 60 5 60 60 65 65

4.6.2 Resultaten De resultaten zijn samengevat in figuur 4.8. Een duidelijke verbetering ten opzichte van de uitgangsituatie is zichtbaar. De verbeteringen op de Nieuwe Merwede en de Beneden Merwede zijn niet zo goed dat aan de calibratie-eis is voldaan...... Figuur 4.8 Afwijking grootschalige morfologie Resultaten calibratie na aanpassing Chézy waarden. 400% 350% 300%

g 250% n i

ijk 200% w f

A 150% 100% 50% 0% Nieuw e Merw ede Beneden Merw ede

Uitgangssituatie (run 1 par. 4.4) Aanpassing Chézy w aarden Eis

Onderzoek naar de invloed op het hydraulische model (bijlage IX) laat zien dat een aanpassing van alleen de Chézy waarden op de Beneden

32 Morfologische gevolgen van rivierverruimende maatregelen

Merwede geen verslechtering in het hydraulische model teweegbrengt. Wanneer echter ook aanpassingen op de Nieuwe Merwede worden doorgevoerd, verslechtert het hydraulische model wel.

4.6.3 Discussie Zowel een aanpassing in korreldiameter als een aanpassing van de Chézy waarden heeft niet geleid tot een bevredigend resultaat voor de Beneden Merwede. Wel leiden beide aanpassingen tot een behoorlijke verbetering ten opzichte van de uitgangssituatie. Dit pleit voor een aanpassing van zowel de korreldiameter als de Chézy waarden in het uiteindelijke model. Voor de Nieuwe Merwede is een aanpassing in korreldiameter wel voldoende geslaagd, de aanpassing in Chézy waarden levert niet het gewenste resultaat op. Bovendien zorgt deze aanpassing voor een verslechtering van het hydraulische model, wanneer ook de Chézy waarden op de Beneden Merwede worden aangepast. De aanpassing op de Beneden Merwede is gezien de resultaten van de aanpassing in korreldiameter noodzakelijk.

4.6.4 Keuze Op basis van bovengenoemde bevindingen is er voor gekozen om op de Beneden Merwede wel een aanpassing in Chézy waarden toe te passen, naast de aanpassing in korreldiameter. Hierbij zijn over de gehele tak constante Chézy waarden gebruikt. De waarden zijn gelijk aan de waarden op het begin van de tak in het oorspronkelijke model (tabel 4.4). Op de Nieuwe Merwede is niet gekozen voor een aanpassing in Chézy waarde, wel voor de aanpassing in korreldiameter. Ten eerste omdat een aanpassing in korrelgrootte al voldoende resultaten oplevert, maar ook omdat een verslechtering van het hydraulische model ongewenst is. In de volgende paragraaf zijn de uiteindelijke keuzes opgesomd en zijn de resultaten gegeven die met de gezamenlijk keuzes behaald zijn.

4.7 Uiteindelijke keuzes

4.7.1 Modelinstellingen Op basis van voorgaande paragrafen is gekomen tot de volgende set van calibratievariabelen: • Transportformule Engelund & Hansen;

• Calibratiefactor transportformule (cf) 0.8 op de onderzoekstakken; • Korreldiameters aangepast ten opzichte van de 5 km gemiddelde waarden (bijlage VIII); • Aangepaste Chézy waarden op de Beneden Merwede, waarden beginpunt oorspronkelijk model nu op gehele tak toegepast (tabel 4.4, kolom 2).

4.7.2 Resultaten In figuur 4.9 zijn de resultaten van de calibratie getoond. Hieruit volgt dat voor drie van de vier takken aan de vooraf gestelde eis is voldaan. Alleen voor de Beneden Merwede is de calibratie niet volledig geslaagd.

33 Morfologische gevolgen van rivierverruimende maatregelen

...... Figuur 4.9 Grafiektitel Resultaten calibratierun met uiteindelijke keuzes. 140% 120%

100%

g 80% in ijk w

f 60% A 40%

20% 0% Waal Boven Merw ede Nieuw e Merw ede Beneden Merw ede

Afw ijking grootschalige morfologie Eis

In figuur 4.10 zijn voor de onderzoekstakken de met bovenstaande instellingen berekende bodemveranderingen vergeleken met de metingen. Een groot deel van de eerdergenoemde probleemtrajecten (paragraaf 4.5) is nu opgelost.

4.7.3 Discussie Gezien de doelstellingen van dit onderzoek, is besloten om met bovenstaande modelinstellingen de maatregelen in het kader van ‘ruimte voor de rivier’ door te rekenen. De twee takken waar de meeste interesse naar uit gaat zijn namelijk de Boven Merwede en de Nieuwe Merwede en hier zijn ook de beste resultaten behaald met de calibratie.

Vanzelfsprekend is het mogelijk een model zodanig te calibreren (bijvoorbeeld door een sterke aanpassing van de korreldiameter), dat het gewenste resultaat wel bereikt wordt. Echter, de aanpassingen van de calibratievariabelen die daarvoor nodig zijn zouden ver buiten de reële marges vallen en zijn dus niet wenselijk. Tijdens de calibratie van dit model is er bewust voor gekozen om niet het model als een black box model te beschouwen, waarbij alleen de uitkomst van het model van belang is, maar om de zich binnen het model afspelende processen zo reëel mogelijk te beschrijven.

Eerder is al ter sprake gekomen dat de aanzienlijke aanpassing in korreldiameter op de Nieuwe Merwede veroorzaakt zou kunnen worden door de aanwezigheid van slib. Het proces zand-slib-segregatie [van Ledden, 2003], waar SOBEK geen rekening mee houdt, lijkt hiervan niet de oorzaak te zijn. Dit proces zorgt er namelijk voor dat de erodeerbaarheid van een bodem afneemt. SOBEK zou dus extra erosie moeten berekenen, maar dat is niet het geval.

De toename van het slibgehalte op het tweede deel van de Nieuwe Merwede [Koomans, 2003] zou een verfijning van korreldiameter teweeg moeten brengen. Dit blijkt echter niet uit de metingen. Uit de metingen volgt dat alleen op het laatste meetpunt op de Nieuwe Merwede de korreldiameter aanzienlijk daalt. In dat opzicht is de

34 Morfologische gevolgen van rivierverruimende maatregelen

toegepaste aanpassing in korreldiameter op de Nieuwe Merwede een beter weergave van de werkelijkheid. Overigens zou ook de onnauwkeurigheid in de gebruikte data een oorzaak kunnen zijn van de problemen. Uit de metingen blijkt dat op het traject aan het eind van de Nieuwe Merwede de bodem erodeert in de periode 1990 – 2000, terwijl aan aangetoonde aanwezigheid van slibjuist doet vermoeden dat er sedimentatie is opgetreden. Verder zou door het slib ook invloed kunnen hebben op de onnauwkeurigheid in baggervolumes.

...... Figuur 4.10 Waal en Boven Merw ede Gemeten bodemverandering (blauw) 2,00 versus de door SOBEK berekende

bodemverandering (paars) tussen 1990 en ) m

2000. (

l i h c s r 0,00 e demv bo

-2,00 913 923 933 943 953 963 kmr

Nie uw e M e r w e de 2,00 ) m ( l i h c s 0,00 er v m de bo

-2,00 961 965 969 973 977 981 kmr

Beneden Merwede 2,00 ) m (

l hi c s 0,00 er v m de bo

-2,00 961 964 967 970 973 976 kmr

Op de Beneden Merwede verliep de calibratie moeizaam. Er is hier niet gekozen voor extra aanpassingen, omdat de focus van dit onderzoek vooral ligt bij de Nieuwe en Boven Merwede. Een mogelijke oorzaak van de problemen op deze tak zou een te grote rekentijdstap kunnen

35 Morfologische gevolgen van rivierverruimende maatregelen

zijn. Hierdoor kan de getijbeweging niet goed worden gemodelleerd, waardoor de berekende transporten te laag uitvallen.

Over het algemeen kan gesteld worden dat versimpeling van de werkelijkheid door een model bijdraagt aan afwijkende resultaten. Het proces zand-slib-segregatie, dat niet in SOBEK is opgenomen, is hier een voorbeeld van, maar ook zou het niet kunnen toepassen van gegradeerd sediment een oorzaak van de problemen kunnen zijn. Het gebruikte SOBEK model was daar niet geschikt voor. In een onderzoek naar gegradeerd sediment bij de Pannerdensche Kop is gebleken dat zulke processen een belangrijke rol kunnen spelen in de morfologische ontwikkeling [van der Werf, 2001]. Verder onderzoek hiernaar is in deze studie niet uitgevoerd.

4.7.4 Mogelijkheden van het model Het gecalibreerde model is in het verdere onderzoek gebruikt bij de bepaling van morfologische effecten van herinrichtingsmaatregelen. De resultaten van de calibratie geven inzicht in de potenties van het model tijdens deze toepassing. Voor alle takken zullen voorspellingen van de effecten op grootschalige morfologie gedaan worden. Voor de Nieuwe Merwede zullen de voorspellingen het meest betrouwbaar zijn. Met de resultaten voor de Beneden Merwede zal voorzichtig moeten worden omgegaan.

Het model is niet gecalibreerd op kleinschalige morfologie, ofwel de bodemontwikkeling binnen een tak. Figuur 4.10 laat zien dat alleen voor de Nieuwe Merwede het model hierop goed presteert. De gemiddelde afwijking bedraagt zo’n 15 cm. Op de Boven Merwede en Waal bedraagt de gemiddelde afwijking zo’n 45 cm. Voor de Beneden Merwede zijn de resultaten slechter, de gemiddelde afwijking hier bedraagt meer dan 50 cm. Met deze afwijkingen moet rekening gehouden worden bij het interpreteren van de effecten van de maatregelen, die door SOBEK worden berekend.

Overigens geldt voor de voorspellingen van de effecten van de voorgestelde maatregelen dat deze altijd in relatieve zin bekeken moeten worden, dat wil zeggen in relatie tot de autonome situatie. Dit is de nulsituatie waarin geen maatregelen worden geïmplementeerd. Alleen op die manier kunnen de morfologische ontwikkelingen bepaald worden die puur het gevolg zijn van de toegepaste maatregelen. Hierdoor wordt de eventuele fout die SOBEK in de nulsituatie berekent niet meegenomen in de resultaten van de berekeningen met de maatregelen.

36 Morfologische gevolgen van rivierverruimende maatregelen

5 Effecten rivierverruimende maatregelen

......

5.1 Inleiding

Het morfologisch SOBEK model voor het NDB, dat in het eerste deel van deze studie is gecalibreerd, is toegepast ter bepaling van de effecten van een tweetal maatregelen in het kader van het project ‘Ruimte voor de Rivier’. Een gedetailleerde beschrijving van deze maatregelen en de schematisering ervan in SOBEK is gegeven in paragraaf 5.2. In paragraaf 5.3 is vervolgens aan de hand van enkele simpele theoretische formuleringen een analyse van de effecten van de maatregelen gegeven. In paragraaf 5.4 zijn vervolgens de resultaten die SOBEK genereert besproken en beoordeeld.

5.2 Maatregelen

5.2.1 Steurgat Algemene beschrijving maatregel Het Steurgat vormt de verbinding tussen de Nieuwe Merwede en de Brabantse Biesbosch. In de huidige situatie wordt het Steurgat afgescheiden van de Nieuwe Merwede door de Biesbosch-sluis. De maatregel ‘Steurgat’ bestaat uit het afvoeren van water vanuit de Nieuwe Merwede via het Steurgat richting de Biesbosch. Er is een aantal varianten ontwikkeld. In deze studie wordt slechts de meest kansrijke variant onderzocht, namelijk ‘Steurgat – aanpassing bedrijventerrein’ [van der Linden, 2003].

...... Figuur 5.1 Maatregel Steurgat, variant ‘aanpassing bedrijventerrein’.

37 Morfologische gevolgen van rivierverruimende maatregelen

Deze variant bestaat uit het verwijderen van de huidige Biesboschsluis en het ten zuiden van de plaats van de huidige sluis bouwen van een nieuwe sluis in combinatie met een overlaat (figuur 5.1). Het geheel wordt zodanig ontworpen dat het Steurgat 1 keer per jaar stroomvoerend zal zijn. De sluis zal in die situaties ook geheel open gezet worden. Het aangrenzende bedrijventerrein en het havengebied moeten zodanig worden aangepast dat bij hoogwater het water onbelemmerd het Steurgat in kan stromen. Tussen de oude en de nieuwe sluis is aan de oostzijde uitbreidingsruimte voor watergebonden bedrijvigheid. In dat geval dient de jachthaven verplaatst te worden naar de nieuwe sluis.

Schematisering in SOBEK In de huidige SOBEK schematisatie is het Steurgat niet verbonden met de Nieuwe Merwede. In de nieuwe situatie komt een nieuw knooppunt in de Nieuwe Merwede, waarmee het Steurgat door middel van een nieuwe tak (Steurgat_a) wordt aangesloten (figuur 5.2). Hiervoor is de tak ‘Nime’ gesplitst in Nime_a en Nime_b. Het knooppunt is gelegd op 1403 meter vanaf het beginpunt van de Nieuwe Merwede. De sluis is in de schematisatie opgenomen als een ‘general structure’ met breedte van 14 m en een drempel op –3.65 NAP. De sluis wordt helemaal open gezet, zodra de waterstand boven de 2.60 m NAP komt en weer gesloten als de waterstand beneden de 2.40 m NAP zakt. De overlaat is geschematiseerd als een ‘weir’ met een drempelhoogte van 2.60 m NAP en een breedte van 300 m. De constructies liggen op 1500 m van het begin van de tak Steurgat_a.

...... Figuur 5.2 Schematisatie verbinding Steurgat met Nieuwe Merwede (Nime). De overlaat en de sluis zijn aangegeven (met *) in de nieuwe tak Steurgat_a.

5.2.2 Rivierverruiming knelpunt Gorinchem Algemene beschrijving In de Boven Merwede bij Gorinchem bevindt zich een versmalling van de rivier, waardoor er problemen ontstaan bij extreme afvoeren. Om dit knelpunt op te lossen zijn er verscheidene varianten in ontwikkeling. In het maatregelenprogramma ruimte voor de rivier [2003] zijn een tweetal varianten genoemd.

38 Morfologische gevolgen van rivierverruimende maatregelen

...... Figuur 5.3 Rivierverruiming knelpunt Gorinchem, variant ‘industrie Avelingen’.

Het onderzoek naar andere varianten is nog bezig en het is daarom niet geheel duidelijk welke maatregel uiteindelijk toegepast gaat worden. In samenspraak met het projectbureau ‘Ruimte voor de Rivier’ [T. van der Linden] is besloten om een onlangs opgestelde variant te onderzoeken, aangezien deze een van de meest waarschijnlijke varianten is van deze maatregel. Dit is de variant ‘industrie Avelingen’ [Lodder en Gao, 2003]. Deze variant betreft aanpassingen aan de noordzijde van de Boven Merwede ter hoogte van het industrieterrein Avelingen (figuur 5.3). De voorlanden worden hier afgegraven en de havens gedempt. Tevens wordt een deel van de industrie verplaatst. Door deze maatregel krijgt het industrieterrein een lange loskade in plaats van het drietal haventjes die er nu aanwezig zijn. Deze maatregel is in feite een zomerbedverbreding met extra een loskade om tegemoet te komen aan de wensen van de industrie. Om de afvoer van water na dit traject zo goed mogelijk te laten verlopen wordt ook een geul gegraven in de uiterwaard, die bij bepaalde afvoeren stroomvoerend zal zijn.

Schematisering in SOBEK Het afgraven van het voorland is in SOBEK gemodelleerd als een verbreding van de hoofdgeul. De mate van verbreding per profiel is af te lezen uit de topografische kaart. In totaal zijn vier profielen aangepast (tabel 5.1).

...... Dwarsprofiel Verbreding [m] Tabel 5.1 Aanpassing profielen bij maatregelen Waal__2 956.2 110 Gorinchem. Waal__2 956.7 180 Waal__2 957.2 190 Waal__2 957.7 290

39 Morfologische gevolgen van rivierverruimende maatregelen

De nevengeul is geschematiseerd als parallelle tak. Om de nevengeul via knooppunten te verbinden met de Waal is de tak ‘Waal__2’ in drie delen opgesplitst (figuur 5.4): 1. Waal__2 (lengte: 43935 m) 2. Waal__3 (lengte: 1498 m) 3. Waal__4 (lengte: 1328 m)

De nevengeul is nu gesitueerd parallel aan tak Waal__3, dit is van kmr 958.2 tot 959.7. De lengte van de nevengeul bedraagt 1750 m. De bodem van de geul ligt op –4 m NAP en is 100 m breed. De breedte van de nevengeul neemt toe met de hoogte, tot 155 m op +1.25 m NAP. Hierboven houdt de geul dezelfde breedte.

...... Figuur 5.4 Schematisatie nevengeul

Om te voorkomen dat de geul te vaak stroomvoerend is, is een inlaatdrempel ontworpen, met de drempelhoogte gelijk aan de hoogte van het maaiveld, die +1.5 m NAP bedraagt.

5.3 Theoretische analyse maatregelen

In deze paragraaf zijn de maatregelen onderworpen aan een theoretische analyse. Met behulp van enkele simpele formuleringen zijn in kwalitatieve en kwantitatieve zin de lange-termijn effecten van de voorgestelde maatregelen bepaald. Deze lange-termijn effecten zijn de wijziging in evenwichtstoestand na toepassing van de maatregel. De aannamen bij deze analyse zijn: • voor de ingrepen zijn de rivieren in hydraulisch en morfologisch evenwicht; • de waterspiegel ver benedenstrooms van de plaats van ingrijpen is constant.

In bijlage X is tevens een analyse opgenomen van de korte termijn effecten. Deze effecten zullen in mindere mate optreden omdat de meeste maatregelen alleen effect hebben bij hoogwater.

5.3.1 Kwalitatieve bepaling effecten Waterbeweging, sedimenttransport en bodemligging zijn onderling gekoppeld en vormen een dynamisch systeem. Dit systeem streeft

40 Morfologische gevolgen van rivierverruimende maatregelen

onder bepaalde condities naar een evenwichtstoestand, waarin in de tijd geen veranderingen meer plaatsvinden. Deze evenwichtstoestand kan worden bepaald met de volgende vergelijkingen [Ribberink et al., 2002a en 2002b]: • de evenwichtswaterstand 1 − n ⎛⎞seq heq = ⎜⎟q (vergelijking 5.1) ⎝⎠m • het evenwichtsverhang 3 n 1 ⎛⎞seq i = (vergelijking 5.2) eq 2 ⎜⎟ Cq⎝⎠m

In bovenstaande formules zijn de volgende symbolen gebruikt: h waterdiepte [m] q debiet per eenheid van breedte [m2/s] i bodemverhang [-] C Chézy coëfficiënt [m1/2/s] s sedimenttransport per eenheid van breedte [m2/s] m, n constanten, hangen af van transportformule

Verder is met het subscript eq de evenwichtstoestand aangeduid. Naast de kwalitatieve beschrijving van de effecten in de volgende twee paragrafen, is in paragraaf 5.3.2 een afschatting gemaakt van de orde van groottes van deze effecten.

Steurgat Alleen in stroomvoerende toestand heeft deze maatregel morfologische effecten tot gevolg. De frequentie waarmee het Steurgat stroomvoerend zal zijn, is dan ook bepalend voor de mate waarin de effecten zullen optreden. Gezien de kwalitatieve aard is in de navolgende analyse uitgegaan van een stroomvoerend Steurgat.

Uit vergelijking 5.1 en 5.2 volgt dat een wateronttrekking tot gevolg heeft dat de waterdiepte afneemt en het verhang toeneemt. Echter, een sedimentonttrekking heeft juist tot gevolg dat de waterdiepte toeneemt en het verhang afneemt. Het is dus noodzakelijk om na te gaan welke van de twee veranderingen de sterkste invloed heeft op de veranderingen van de evenwichtstoestand. Uit de vergelijkingen volgt: 1 − hq∼ , sn

3 iq∼ −1, sn Voor n=5 (is bijvoorbeeld het geval in de transportformule van Engelund & Hansen) blijkt dat een bepaalde relatieve verandering in het debiet een grotere invloed heeft op de waterdiepte en het verhang dan een zelfde relatieve verandering in sediment transport. Bovendien zal de sedimentonttrekking kleiner zijn dan verwacht mag worden op basis van de verandering in debiet. Er geldt niet dat bijvoorbeeld een debietonttrekking van 50% van het oorspronkelijke debiet, een sedimenttonttrekking van 50% van het oorspronkelijke transport met zich meebrengt. Door de aanwezigheid van een drempel

41 Morfologische gevolgen van rivierverruimende maatregelen

bij de inlaat naar het Steurgat zal vooral water onttrokken worden uit de bovenste waterlagen waar de sedimentconcentratie het laagst is. De aanwezigheid van bodemtransport versterkt dit argument.

De nieuwe evenwichtssituatie wordt dus hoofdzakelijk bepaald door de debietonttrekking, wat neerkomt op een afname van waterdiepte en een toename van verhang benedenstrooms van het Steurgat (vergelijking 5.1 en 5.2)...... Figuur 5.5 Aanpassing evenwichtssituatie na h0 , i0 toepassing maatregel Steurgat. De h1 < h0 , i1 > i0 gestippelde lijnen geven de Q0 , S0 oorspronkelijke waterspiegel en bodemligging weer.

gefixeerde waterstand

Steurgat: onttrekking debiet

In figuur 5.5 is een schets gemaakt van de vernieuwde evenwichtssituatie. Hierbij is uitgegaan van een gefixeerde waterstand (ver) benedenstrooms van de ingreep. Uit de figuur volgt dat de effecten benedenstrooms van de ingreep ook gevolgen hebben voor het bovenstroomse gedeelte. Aangezien hier de waterspiegel ook gestegen is, moet de bodem meestijgen om dezelfde waterdiepte te behouden. De stijging van de bodem hier is gelijk aan L·∆i, met L de lengte vanaf het Steurgat tot aan de locatie met de gefixeerde waterstand (bijvoorbeeld de zee) en ∆i het verschil in verhang op dit traject voor en na de toepassing.

Rivierverruiming knelpunt Gorinchem Het voorgenomen ingrijpen bij Gorinchem is feitelijk een pakket van maatregelen. Enerzijds zal het aanleggen van een extra geul in de uiterwaard zorgen voor een debietonttrekking uit de hoofdgeul, anderzijds wordt een deel van de noordelijke uiterwaard afgegraven waardoor het zomerbed wordt verbreed. Beide maatregelen zorgen ervoor dat het debiet per eenheid van breedte (q) afneemt.

Figuur 5.6 toont de vernieuwde evenwichtssituatie. Benedenstrooms van het traject van ingrepen vinden geen veranderingen plaats in de evenwichtstoestand. Hier zijn immers geen wijzigingen in debiet en sedimenttransport. Op het traject zelf daalt de waterdiepte door de afname in q (vergelijking 5.1), wat neerkomt op een stijging van de bodem. Bovendien is het verhang hier groter (vergelijking 5.2), waardoor de bodemstijging stroomopwaarts toeneemt. Bovenstrooms van het traject is het debiet en het sedimenttransport ook niet gewijzigd. De waterdiepte en het verhang blijven hier dus ongewijzigd. Echter, door het grotere bodemverhang op het maatregelentraject, zal de bodem hier stijgen om dezelfde waterdiepte te behouden. De stijging van de bodem hier is gelijk aan L·∆i, met L de lengte van het

42 Morfologische gevolgen van rivierverruimende maatregelen

maatregelentraject (± 5000 m) en ∆i het verschil in verhang op het maatregelentraject voor en na de toepassing...... Figuur 5.6 Aanpassing evenwichtssituatie na toepassing rivierverruimende maatregelen Gorinchem. De gestippelde h0 , i0 h1 < h0 , i1 > i0 lijnen geven de oorspronkelijke waterspiegel en bodemligging weer.

Q0 , S0 h0 , i0

gefixeerde waterstand

Traject maatregelen

5.3.2 Kwantificering effecten De analyse in voorgaande paragrafen is van kwalitatieve aard. In de huidige paragraaf worden met behulp van enkele simpele berekeningen de orde van groottes afgeschat van de effecten. Ook de rol van de morfologische tijd- en ruimteschaal van de voorspelde veranderingen wordt toegelicht.

Steurgat Uit figuur 5.5 volgt dat benedenstrooms van het Steurgat de bodem zal gaan stijgen. Deze stijging vindt plaats over het gehele traject dat tussen het Steurgat en de locatie met de gefixeerde waterstand ligt. Normaal gesproken ligt deze locatie in zee, wanneer verondersteld wordt dat de zeespiegel niet of nauwelijks verandert. In deze situatie kan het Hollandsch Diep als rand gebruikt worden, aangezien het water dat via het Steurgat wordt afgevoerd uiteindelijk via de Biesbosch weer terechtkomt in het Hollandsch Diep. In het Hollandsch Diep is dan ten opzichte van de autonome situatie niets veranderd. In dat geval is de morfologische ruimteschaal zo’n 20 kilometer (lengte Nieuwe Merwede). Nu kan een schatting gemaakt worden van de veranderingen in evenwichtssituatie.

Benedenstrooms van het Steurgat geldt: QQ− ∆ q = 0 (vergelijking 5.3) 1 B

Met ∆Q de debietonttrekking via het Steurgat. Tezamen met vergelijkingen 5.1 en 5.2 volgt: ⎛⎞∆Q hh10=−⎜1 ⎟ (vergelijking 5.4) ⎝⎠Q0

Q0 ii10= (vergelijking 5.5) QQ0 − ∆ Het gemiddelde debiet door de Nieuwe Merwede over de periode 1990-2000 is zo’n 930 m3/s (bepaald uit SOBEK resultaten). Wanneer het Steurgat zodanig wordt ontworpen dat de overschrijdingsfrequentie van de drempel 1 keer per jaar is, verandert dit gemiddeld debiet

43 Morfologische gevolgen van rivierverruimende maatregelen

slechts met 0.11%. Er volgt dan dat h1 = 0.99·h0 en i1 = 1.00·i0. De lange termijn effecten zijn dus zeer minimaal. Met een oorspronkelijke waterdiepte van de orde van grootte van 5 m, levert dit een daling van de waterdiepte op met 0.6 cm op de benedenstroomse rand van de Nieuwe Merwede, ofwel een bodemstijging van 0.6 cm. Er is hier geen rekening gehouden met het feit dat de grootste morfologische veranderingen plaatsvinden bij hoge debieten. De effecten zullen daarom onderschat zijn. Ook zullen er wel enige korte-termijn effecten zichtbaar zijn, elke keer dat het Steurgat stroomvoerend is.

De tijdschaal waarop bovengenoemde lange termijn effecten zich afspelen is moeilijk in te schatten. De bodem zal zijn evenwichtstoestand waarschijnlijk ook nooit bereiken. Door steeds terugkerende aanpassingen in de riviersystemen krijgen de systemen te weinig tijd om evenwichtssituaties te bereiken. Door baggerwerkzaamheden bijvoorbeeld zal de bodemstijging beperkt worden. Het is echter wel zo dat de morfologische tijdschaal waarop de effecten zich afspelen wel informatie verschaft over de benodigde intensivering van de baggerwerkzaamheden om de bodem op gewenst peil te houden. Hierbij kan SOBEK een goed hulpmiddel zijn om deze tijdschaal in te schatten.

Rivierverruiming knelpunt Gorinchem Een soortgelijke analyse is uitgevoerd voor de maatregelen bij Gorinchem. De zomerbedverbreding en de nevengeul zijn apart bekeken. De gemiddelde verbreding van het zomerbed bedraagt 200 m, wat neerkomt op een toename van ongeveer 50%. Bij gelijkblijvend debiet neemt het debiet per eenheid van breedte dus af met een factor 1/1.5 = 2/3. Uit vergelijking 5.1 volgt dan dat de evenwichtswaterdiepte hiermee evenredig afneemt, wat tot uiting komt in een bodemstijging. Bij een huidige gemiddelde waterdiepte van 6 m levert dit een bodemstijging van 2 m. Het evenwichtsverhang neemt toe (vergelijking 5.2) met factor 1.5. Het bodemverhang op dit traject is echter dusdanig klein dat de effecten hiervan verwaarloosbaar zijn. Voor de maatregel ‘nevengeul’ geldt dat met een drempelhoogte van 1.5 m NAP, de geul zo’n 20 dagen per jaar stroomvoerend zal zijn. Dit is het geval bij een debiet door de hoofdgeul van ongeveer 2900 m3/s. Ervan uitgaande dat het extra debiet boven deze 2950 m3/s verdeeld wordt over de hoofdgeul en nevengeul volgens de breedteverhoudingen (400 : 150) van de twee geulen, levert dit een daling van het gemiddeld debiet door de hoofdgeul op van 10 m3/s. Met een oorspronkelijk debiet van zo’n 1500 m3/s volgt uit vergelijking 5.4 dat de waterdiepte afneemt met 4 cm. Vergeleken met het effect door de zomerbedverbreding is dit dus verwaarloosbaar.

Door het zeer kleine verhang in bodemverloop in de Boven Merwede zullen de effecten hier dus beperkt blijven tot het traject waar de maatregelen gepland zijn. Op het traject van de zomerbedverbreding zal een forse bodemstijging plaatsvinden. Verder benedenstrooms ter hoogte van de nevengeul zal de bodemstijging beduidend lager uitvallen. De snelheid waarmee de forse bodemstijging zich manifesteert is afhankelijk van de tijdschaal. Bij een zeer grote tijdschaal zal dit proces

44 Morfologische gevolgen van rivierverruimende maatregelen

zo traag verlopen, dat deze sedimentatie beperkt blijft. Andere ingrepen in het riviersysteem zullen de evenwichtssituatie wederom veranderen. Uiteraard spelen baggerwerkzaamheden hier ook weer een grote rol in. Wordt de sedimentatie consequent weggebaggerd, dan zal deze sedimentatie zich nooit verder kunnen uitbreiden.

Korte termijn effecten Voor beide maatregelen geldt dat er zich ook korte termijn effecten zullen voordoen ten gevolge van de aanpassingen (bijlage X). De morfologische tijdschaal hiervan is korter dan de tijdschaal van de lange termijn effecten. Normaalgesproken zijn de korte termijn effecten al vrij snel merkbaar na toepassing van de ingrepen (<< 1 jaar). Voor het Steurgat en de nevengeul bij Gorinchem geldt dat beide aanpassingen alleen benut worden bij extreme afvoeren. Deze effecten zullen zich dan ook alleen na extreme afvoeren voordoen. Bovendien zullen de effecten na een periode met extreme afvoeren weer min of meer uitgewist worden, doordat het systeem nu weer terugkeert naar de onveranderde situatie. Voor de zomerbedverbreding zullen de korte termijn effecten zich waarschijnlijk veel sterker uiten. De voorspelde bodemgolven zullen zich ook direct na de ingreep ontwikkelen.

SOBEK zal ook hier weer als hulpmiddel gebruikt worden bij het afschatten in welke mate de effecten zich voordoen.

5.3.3 Conclusies Het handmatig afschatten van de voorspelde effecten geeft een globaal inzicht in de verwachte uitkomsten. Met het in de voorgaande hoofdstukken ontwikkelde SOBEK model is zo’n afschatting realistischer, gezien de complexiteit van het rivierensysteem. Hierbij is het wel van belang om een voldoende ruime rekenperiode te hanteren. De effecten zullen zich namelijk vrij langzaam manifesteren aangezien de meeste maatregelen slechts sporadisch effect zullen hebben.

5.4 Resultaten SOBEK

5.4.1 Inleiding Dit hoofdstuk beschrijft de effecten van de maatregelen die met behulp van SOBEK zijn bepaald. Tijdens de SOBEK runs zijn de volgende instellingen gebruikt: • rekenperiode 50 jaar (afweging benodigde rekentijd en gewenste simulatieperiode) met dezelfde initiële morfologische condities als tijdens de calibratie (startbodem 1990); • zelfde hydraulische randvoorwaarden aan de bovenrand als tijdens calibratie, echter nu 5 x achter elkaar. Deze periode is lang genoeg voor een betrouwbare afvoerverdeling; • hydraulische benedenrand is een morfologische bepalend getij, geen zeespiegelrijzing; • baggeren toegepast, middels een over 1990 – 2000 gemiddeld continue sedimentonttrekking (figuur 5.7).

De resultaten van SOBEK zijn in relatie tot de nulsituatie geanalyseerd. De nulsituatie is bepaald door het uitvoeren van een autonome run (paragraaf 5.4.2). De effecten van de maatregel ‘Steurgat’ en

45 Morfologische gevolgen van rivierverruimende maatregelen

‘Rivierverruimende maatregelen Gorinchem’ worden respectievelijk behandeld in paragraaf 5.4.3 en 5.4.4.

...... Figuur 5.7 Baggeren tijdens simulatieperiode Baggerwerkzaamheden die gebruikt zijn tijdens de berekeningen met de rivierverruimende maatregelen. 3.000 ] aar

j 2.500

/10 3 2.000 00 m 1.500 me [10 u l

o 1.000 v

gger 500 a B

0 Waal Boven Merw ede Beneden Merw ede Nieuw e Merw ede

5.4.2 Nulsituatie De grootschalige morfologische ontwikkeling van de vier onderzoekstakken is weergegeven in figuur 5.8. Verder zijn in bijlage XI figuren opgenomen met daarin de kleinschalige morfologische bodemontwikkeling, ofwel de bodemveranderingen per locatie binnen een tak. Het model is hierop niet afgeregeld, dus de resultaten hiervan zijn niet volledig betrouwbaar. Maar juist door het bepalen van de kleinschalige morfologische ontwikkeling van de nulsituatie, kan bij berekeningen met de maatregelen alleen het kleinschalige morfologische effect ten gevolge van de maatregel zelf bestudeerd worden. De onbetrouwbaarheid in de autonome ontwikkeling wordt hierin dus niet meegenomen omdat de resultaten van de berekeningen met de maatregelen altijd ten opzichte van de nulsituatie bekeken worden. Oftewel: effect maatregel = totaal effect – autonome ontwikkeling. Verder geldt dat voor de Nieuwe Merwede het model tijdens de calibratie goede resultaten behaald zijn met betrekking tot de kleinschalige morfologische ontwikkeling. Voor de Beneden Merwede zijn de resultaten van zowel grootschalige als kleinschalige morfologische ontwikkeling het minst betrouwbaar.

Aan de hand van bovengenoemde figuren is hieronder per tak globaal aangegeven hoe de bodemligging zich zal ontwikkelen gedurende 50 jaar, zonder toepassing van maatregelen.

46 Morfologische gevolgen van rivierverruimende maatregelen

...... Figuur 5.8 Nulsituatie De totale bodemverandering per tak voor de situatie waarin geen maatregelen zijn 2.000 toegepast. 1.000 ] 3 0

-1.000 10 ng [1000 m 25 -2.000 50 veranderi -3.000 odem B -4.000

-5.000 Beneden Merwede Nieuwe Merwede Boven Merwede Waal

Waal De resultaten laten over het algemeen erosie zien (figuur 5.8). Gedurende de eerste 10 jaar is al de helft van de uiteindelijke erosie opgetreden. Op bepaalde locaties is deze erosie vrij fors met bodemdalingen van 1 meter na 50 jaar. Ook vindt er op enkele punten sedimentatie plaats, tot maximaal 60 cm. Op het bovenstroomse deel van de Waal (kmr 914 – 934) is de bodem vrij snel in evenwicht gekomen. In de laatste 40 jaar van de simulatie is de bodem hier nauwelijks meer veranderd. Benedenstrooms geldt dat evenwicht later wordt bereikt, hier is de bodemligging de laatste 25 jaar bijna niet meer veranderd.

Boven Merwede: Op deze tak vindt er netto erosie plaats. Wel blijkt dat op de eerste 2.5 km de bodem stijgt met zo’n 20 cm. Verderop daalt de bodem juist met gemiddeld 40 cm. De bodem lijkt nog niet een evenwichtstoestand bereikt te hebben; de verandering in de laatste 25 jaar zijn namelijk vrij groot.

Nieuwe Merwede De bodem van deze tak daalt behoorlijk. Na 25 jaar is de bodem met gemiddeld 35 cm gedaald en na 50 jaar is de totale gemiddelde daling 50 cm. De snelheid waarmee de bodem daalt neemt dus wel af doordat de bodem de evenwichtssituatie nadert, maar de evenwichtsstand lijkt na 50 jaar nog niet gehaald.

Beneden Merwede De enige tak waarop sedimentatie plaatsvindt. Op de eerste 6 km van de Beneden Merwede is de bodem gedurende de 50 jaar lichtelijk geërodeerd, gemiddeld met zo’n 25 cm. Tussen kmr 967 en 969 vindt juist sedimentatie plaats. De bodem is hier tot 80 cm gestegen. Vanaf kmr 970 tot aan het eind is de bodem ook gestegen met gemiddeld 1 meter. Opvallend is dat de bodem gedurende de laatste 25 jaar van de simulatie nauwelijks meer veranderd is. De evenwichtssituatie is dus vrij snel bereikt. Of er werkelijk sedimentatie zal optreden is onzeker doordat de calibratie van deze tak niet goed verliep. Tijdens de calibratie berekende

47 Morfologische gevolgen van rivierverruimende maatregelen

het model namelijk veel sedimentatie, terwijl uit de metingen bleek dat de bodem licht erodeerde. De afwijking van de grootschalige morfologie tijdens de calibratie op deze tak bedroeg 550.000 m3 (figuur 4.9: 123% x baggervolume = 123% x 450.000 m3) in 10 jaar. Uit figuur 5.8 volgt dat een sedimentatie van dezelfde orde van grootte zich voordoet in de eerste 10 jaar van de nulsituatie. Naar alle waarschijnlijkheid zal de werkelijke morfologische veranderingen dus vrij klein zijn. Hieruit blijkt wederom dat, vooral voor de Beneden Merwede, de resultaten van de maatregelen altijd relatief ten opzichte van de nulsituatie bekeken moeten worden.

5.4.3 Steurgat Deze maatregel heeft zeer beperkte morfologische effecten. De effecten op de Beneden Merwede, de Boven Merwede en de Waal zijn verwaarloosbaar klein. Op de Nieuwe Merwede zijn wel enige verschillen te zien in grootschalige morfologie ten opzichte van de nulsituatie (figuur 5.9). De effecten zijn hier echter ook gering.

...... Figuur 5.9 Nieuwe Merwede Totale bodemverandering voor Nieuwe Merwede na toepassen van de maatregel 0 ‘Steurgat’. -500 ]

3 -1000

-1500

1000 m -2000 Nulsituatie -2500 Steurgat -3000 verandering [ -3500

Bodem -4000

-4500

-5000 na 10 jaar na 25 jaar na 50 jaar

Figuur 5.10 geeft een goed beeld van de zich steeds herhalende korte termijneffecten die zich voordoen rond extreme afvoeren, wanneer het Steurgat stroomvoerend is. De sedimentatiegolven die in figuur VIII.1 (bijlage X) zijn weergegeven komen in de SOBEK resultaten ook goed naar voren. De golven verplaatsen zich in benedenstroomse richting. De hoogte van de golven zijn in de figuur zo’n + 12 cm ten opzichte van de run met de nulsituatie. Uiteraard is de figuur een momentopname en zijn hogere golven mogelijk, afhankelijk van de afvoer door het Steurgat. Het maximale verschil in bodemligging ten opzichte van de nulsituatie dat is voorgekomen in de 50-jarige simulatie is + 13 cm op kmr 963 (net na aftakking naar Steurgat).

Verder valt op dat na 50 jaar de bodem van bijna de gehele Nieuwe Merwede met zo’n 2 – 4 cm is gestegen, ruim meer dan voorspeld in paragraaf 5.3.2. Hieruit blijkt weer dat de grootste morfologische veranderingen zich voordoen bij hoge afvoeren.

48 Morfologische gevolgen van rivierverruimende maatregelen

...... Figuur 5.10 Effecten maatregelen Steurgat op Nime Effecten bodemligging Nieuwe Merwede tgv de maatregel ‘Steurgat’ tov de 0,14

] na 50 jaar nulsituatie. m 0,12 [ na 25 jaar e i 0,10

uat na 10 jaar t i s 0,08 nul

v

o 0,06 t l i h c

s 0,04 r e 0,02

demv 0,00 bo -0,02 961 963 965 967 969 971 973 975 977 979 981 kmr

Uit bovenstaande blijkt dat de effecten behoorlijk klein zijn. Figuur 5.11 laat de autonome bodemontwikkeling (grijze lijn) zien en ook de extra effecten die optreden ten gevolge van de maatregel (zwarte lijn). Hieruit blijkt dat zich in de autonome situatie, in een relatief korte tijd van enkele weken, bodemveranderingen van aanzienlijk grotere orde voordoen dan de extra effecten ten gevolge van de maatregel. Het effect ten gevolge van de maatregel na een hoge afvoer bedraagt dan wel zo’n 12 cm ten opzichte van de nulsituatie, de bodemverandering die sowieso al optreedt bedraagt ruim 30 cm.

...... Figuur 5.11 Effect Steurgat t.o.v. autonome bodemontwikkeling Vergelijking bodemontwikkeling nulsituatie, gezien vanaf t=0 (blauwe lijn) met de extra effecten veroorzaakt door de 0,50 maatregel ‘Steurgat’ (gele lijn) op een ] locatie in de Nieuwe Merwede net na de m

[ 0,25 l aftakking naar het Steurgat, waar de korte- i h c

effecten maximaal zijn. s r

e 0,00

odemv -0,25 B

-0,50 4 9 4 9 4 9 4 9 0 0 1 1 2 2 3 3 90 95 00 ------an- an- an- j j j dec dec dec dec dec dec dec dec Datum

5.4.4 Rivierverruiming knelpunt Gorinchem Waal Op de Waal zijn geen effecten waarneembaar ten gevolge van de maatregelen. De effecten op bodemligging verplaatsen zich dus niet in bovenstroomse richting. Dit kwam ook naar voren uit de theoretische analyse.

Boven Merwede Uit figuur 5.12 volgt dat er een behoorlijke sedimentatie optreedt ten gevolge van de maatregelen. De sedimentatie voltrekt zich snel, na 10 jaar al is 2/3 van de uiteindelijke sedimentatie bereikt. De totale sedimentatie na 10 jaar bedraagt 1.200.000 m3. Tijdens de modelruns is

49 Morfologische gevolgen van rivierverruimende maatregelen

rekening gehouden met een toekomstige baggerwerkzaamheden van 2.500.000 m3 per 10 jaar, dat gebaseerd is op de baggerwerkzaamheden tijdens de periode 1990 – 2000. Om de bodem na toepassing van de voorgestelde maatregelen op peil te houden, zal dit baggerbeleid dus moeten worden geïntensiveerd tot 3.700.000 m3 per 10 jaar, een stijging met 50%.

...... Figuur 5.12 Boven Merwede Totale bodemverandering voor Boven Merwede na toepassen van de 2500 rivierverruimende maatregelen bij Gorinchem 2000

1500

1000 Nulsituatie Gorinchem 500

0 Bodemverandering [1000 m3] -500

-1000 na 10 jaar na 25 jaar na 50 jaar

Figuur 5.13 geeft een meer gedetailleerd inzicht in de effecten. Vrijwel de gehele sedimentatie doet zich voor tussen kmr 956 en 958, het tracé waar de zomerbedverbreding heeft plaatsgevonden. De effecten hiervan zijn zodanig groot, dat de effecten van de nevengeul nauwelijks zichtbaar zijn. Ter plekke van de nevengeul is na 50 jaar de bodemligging zo’n 10 cm hoger dan in de nulsituatie, terwijl ter plekke van de zomerbedverbreding de bodem uiteindelijk ruim 2 meter hoger komt te liggen. De overeenkomst met de theoretische analyse is dus vrij goed. Alleen rond kmr 956 verloopt de bodemsprong veel geleidelijker dan uit figuur 5.6 volgt. Dit is te wijten aan het feit dat er door SOBEK gemiddelde waarden worden berekend op een rekenpunt (er wordt gemiddeld over het halve vak erna en het halve vak ervoor). Rond kmr 958 treedt dit probleem niet op, omdat zich hier een knooppunt bevindt in verband met de aansluiting van de nevengeul.

...... Figuur 5.13 Effecten maatregelen Gorinchem op Boven Merwede Effecten bodemligging Boven Merwede tgv 4,00 de maatregelen bij Gorinchem tov de na 50 jaar nulsituatie. De stippellijnen geven het gebied ] van de zomerbedverbreding aan. m 3,00 na 25 jaar e [ i

t na 10 jaar a u

t 2,00 i na 6 jaar s l na 3 jaar

nu 1,00 v o t l i

h 0,00 c s r e -1,00

demv -2,00 bo -3,00 4 5 6 7 8 9 0 1 2 95 95 95 95 95 95 96 96 96 kmr

50 Morfologische gevolgen van rivierverruimende maatregelen

In de figuur zijn ook de effecten na 3 en 6 jaar weergegeven. Hier zijn de korte termijn effecten zichtbaar. Deze effecten komen goed overeen met de verwachte effecten op basis van de theoretische analyse (bijlage X): een geconcentreerde sedimentatie aan het begin van het maatregelentraject en een geconcentreerde erosie (-1.5 m tov nulsituatie) aan het eind van het traject van zomerbedverbreding. Deze effecten verplaatsen zich vervolgens in benedenstroomse richting als bodemgolven. Na 10 jaar bereikt de erosiegolf het einde van de Boven Merwede. De vorm van de erosiegolf is dan veranderd, minder diep en breder. De sedimentatiegolf is minder goed zichtbaar. Dit is voornamelijk te wijten aan de korte afstand tussen de twee bodemvormen, waardoor de sedimentatiegolf eigenlijk zorgt voor de opvulling van de erosiegolf.

Nieuwe Merwede De effecten van de maatregelen in de Boven Merwede zijn ook merkbaar op de takken benedenstrooms ervan. De zich voortplantende erosiegolf op de Boven Merwede heeft na 10 jaar de Nieuwe en Beneden Merwede bereikt. Door de verdeling over de twee takken en door diffusie is de golf duidelijk verondiept (figuur 5.14). De bodemligging in het dal van de golf is zo’n 80 cm lager in vergelijking met de nulsituatie, terwijl dit op de Boven Merwede nog 1.5 m bedroeg.

...... Figuur 5.14 Effecten maatregelen Gorinchem op Nieuwe Merwede Effecten bodemligging Nieuwe Merwede tgv de maatregelen bij Gorinchem tov de 0,20 nulsituatie. ] m [ e i 0,00 uat t i s

nul -0,20

v o t l i h

c -0,40 s r e -0,60 na 50 jaar na 25 jaar bodemv na 10 jaar -0,80 961 963 965 967 969 971 973 975 977 979 981 kmr

Naast de erosie die in de nulsituatie optreedt, ontstaat door de golf een extra erosie van 500.000 m3 die na de volledige intrede van de golf vrij constant blijft (figuur 5.15). Een deel van de extra erosie zal waarschijnlijk verdwijnen nadat de golf de Nieuwe Merwede heeft verlaten, maar dit is na 50 jaar nog niet het geval. Het andere deel van de extra erosie blijft, bovenstrooms van de golf (na 50 jaar) is de bodemligging namelijk lager dan in de nulsituatie. Een verklaring hiervoor zou kunnen zijn dat in de nulsituatie de bodem van de Nieuwe Merwede ook erodeert. Na 50 jaar is de evenwichtssituatie nog niet bereikt en dus zal de bodem hierna nog verder dalen. Door de passage van de erosiegolf is dit erosieproces versneld en is de bodem dichter tot de evenwichtssituatie genaderd.

51 Morfologische gevolgen van rivierverruimende maatregelen

...... Figuur 5.15 Nieuwe Merwede Totale bodemverandering voor Nieuwe Merwede na toepassen van de 0 rivierverruimende maatregelen bij Gorinchem -1000

-2000

Nulsituatie -3000 Gorinchem

-4000

Bodemverandering [1000 m3] -5000

-6000 na 10 jaar na 25 jaar na 50 jaar

Beneden Merwede Door de passage van de bodemgolf zal de voorspelde sedimentatie die in de nulsituatie naar voren komt, veel lager uitvallen (figuur 5.16).

...... Figuur 5.16 Beneden Merwede Totale bodemverandering voor Beneden Merwede na toepassen van de 1800 rivierverruimende maatregelen bij 1600 Gorinchem ] 3 1400

1200

ng [1000 m 1000 Nulsituatie 800 Gorinchem

veranderi 600

odem 400 B 200

0 na 10 jaar na 25 jaar na 50 jaar

Hoe de morfologische ontwikkeling zich na de eerste 10 jaar voortzet is met het huidige model lastig te bepalen. In de nulsituatie blijft de bodem sedimenteren gedurende de gehele simulatieperiode. Na 50 jaar is ook de sedimentatie in de situatie met de maatregelen bij Gorinchem aanzienlijk toegenomen. Waarschijnlijk zal de na 25 jaar volgens het model opgetreden sedimentatie in werkelijkheid niet optreden, maar zal er juist erosie optreden door de passage van de erosiegolf. Na 50 jaar zal de bodem weer terugkeren naar de oude situatie en zal de sedimentatie of erosie bij voortzetting van het huidige baggerbeleid klein zijn.

Na 50 jaar heeft de erosiegolf het einde van de Beneden Merwede bereikt. De bodemligging is dan gemiddeld zo’n 20 cm lager dan in de nulsituatie. Naar alle waarschijnlijkheid zal de bodem hierna weer stijgen, aangezien de evenwichtstoestand van de bodem nog niet bereikt is. De snelheid waarmee de bodemgolven zich over de Beneden en Nieuwe Merwede voortplanten zijn vrijwel gelijk. Aan het eind van de

52 Morfologische gevolgen van rivierverruimende maatregelen

simulatie liggen de dalen van beide golven in beide takken net na kmr 975. De oorsprong van de golf ligt ter hoogte van kmr 958, dus de voortplantingssnelheid bedraagt zo’n 330 meter per jaar.

...... Figuur 5.17 Effecten maatregelen Gorinchem op Beneden Merwede Effecten bodemligging Beneden Merwede tgv de maatregelen bij Gorinchem tov de 0,20 nulsituatie. ] [m

e 0,00 i t a u t i s

l -0,20 nu v o

t -0,40

l i h c s -0,60 er

mv na 50 jaar

de -0,80 na 25 jaar bo na 10 jaar -1,00 1 3 5 7 9 1 3 5 7 96 96 96 96 96 97 97 97 97 kmr

Effecten benedenstroomse takken Uit de resultaten blijkt dat erosiegolf die veroorzaakt is door de zomerbedverbreding op de Boven Merwede na 50 jaar nog duidelijk aanwezig is. Deze golf zal zeker ook de takken benedenstrooms van de Nieuwe en Beneden Merwede bereiken. Het model is hier niet gecalibreerd en daarom zijn deze effecten niet in kaart gebracht.

53 Morfologische gevolgen van rivierverruimende maatregelen

54 Morfologische gevolgen van rivierverruimende maatregelen

6 Conclusies en aanbevelingen

......

6.1 Conclusies

De doelstelling van dit onderzoek was het opzetten en calibreren van het morfologische NDB model met behulp van SOBEK teneinde een tweetal maatregelen in het kader van het project ‘Ruimte voor de Rivier’, die worden toegepast in de Nieuwe en Boven Merwede, te onderzoeken op hun grootschalige morfologische effecten. Hieronder zijn de conclusies van dit onderzoek gegeven. Allereerst volgen enkele conclusies met betrekking tot de calibratie van het morfologische model: • De calibratie van het model is uitgevoerd voor de drie Merwede takken en de Waal. Het model is afgeregeld op de grootschalige morfologie, ofwel de totale bodemverandering (m3) per tak in de periode 1990–2000. In tegenstelling tot het voorgaande morfologische model is dit model robuuster; er zijn weinig problemen geconstateerd tijdens de simulaties over de 10-jarige periode. • Het model is afgeregeld met (calibratievariabelen) de keuze van de transportformule, de calibratiefactor van de transportformule, Chézy waarden en korreldiameters. De keuze van transportformule heeft een relatief grote invloed op de modeluitkomsten. De Chézy waarden en korreldiameters zijn geschikte variabelen om de morfologische ontwikkeling binnen een tak af te regelen. • Als calibratie-eis is gehanteerd dat de modelvoorspelling moeten overeenkomen met de metingen, rekening houdend met een marge van 30% van de hoeveelheid baggerwerk. • De calibratie voor de Beneden Merwede is niet geslaagd conform de vooraf gestelde eisen. Voor de overige onderzoekstakken is dit wel gelukt. De beste resultaten zijn behaald voor de Nieuwe Merwede. • De problemen met de Beneden Merwede uitte zich in een sedimentatie op de gehele tak, terwijl uit de metingen juist volgde dat de bodem onveranderd bleef. Een oorzaak zou kunnen zijn dat de gebruikte rekentijdstap te groot is gekozen voor een goede modellering van de getijbeweging. • Voor het slagen van de calibratie in de Nieuwe Merwede was een behoorlijke aanpassing in korreldiameter op de laatste 10 km nodig. Uit onderzoek is gebleken dat er zich op dit traject van de Nieuwe Merwede slib bevindt. • Tijdens de calibratie is gebleken dat de baggeractiviteiten in de onderzoekstakken een significante rol op de morfologische ontwikkeling hebben. Tevens is de relatieve onzekerheid in de baggerdata aanzienlijk. Hierdoor wordt de onzekerheid op de grootschalige morfologische ontwikkeling voornamelijk bepaald door de onzekerheid in de baggerdata. De calibratie-eis is dan ook opgesteld op basis van deze onzekerheid.

55 Morfologische gevolgen van rivierverruimende maatregelen

De conclusies met betrekking tot de effecten van de rivierverruimende maatregelen zijn: • Het model is in staat gebleken om een periode van 50 jaar door te rekenen. Hierbij hebben zich geen grote problemen of modelcrashes voorgedaan, waardoor het model een robuust hulpmiddel is bij het bepalen van de effecten van ingrepen in de rivier. • De morfologische effecten van twee maatregelen zijn met behulp van SOBEK bepaald. De eerste maatregel betreft het afvoeren van water vanaf de Nieuwe Merwede door het Steurgat naar de Biesbosch bij extreme afvoer. De andere maatregel bestaat uit twee ingrepen in de Boven Merwede bij Gorinchem. Het gaat hierbij om de optimalisatie van de bestaande nevengeul aan de noordoever en tevens een zomerbedverbreding. De effecten van de maatregelen zijn in relatieve zin beschreven. Dat wil zeggen dat de effecten altijd ten opzichte van de bodemontwikkelingen in de nulsituaties geanalyseerd worden. • SOBEK berekent in de nulsituatie na 50 jaar erosie op de Nieuwe Merwede, de Boven Merwede en de Waal. Op de Nieuwe Merwede is deze erosie het grootst. De bodem is na 50 jaar gemiddeld met zo’n 50 cm gedaald. Op de Beneden Merwede wordt sedimentatie berekend. Echter, de sedimentatie is na 10 jaar van dezelfde orde van grootte als de overschatting van sedimentatie tijdens de calibratie. Naar alle waarschijnlijkheid blijft de bodem dan ook vrij stabiel, bij de voortzetting van het huidige baggerbeleid. • De maatregel betreffende het Steurgat brengt in de onderzochte vorm geen significante morfologische effecten met zich mee op de takken Beneden Merwede, Boven Merwede en de Waal. • Op de Nieuwe Merwede zijn ten gevolge van deze maatregelen wel enige effecten zichtbaar. Ten opzichte van de nulsituatie is de bodem op een groot deel van de tak na 50 jaar met zo’n 4 cm gestegen. Verder komen er na perioden van extreme afvoer, waarbij de maatregel actief wordt, sedimentatiegolven voor, die zich voortplanten in benedenstroomse richting. In vergelijking met de bodemveranderingen die zich in de nulsituatie voordoen zijn bovengenoemde effecten echter gering. • De maatregelen in de Boven Merwede bij Gorinchem veroorzaken grote morfologische effecten. Dit is vooral te wijten aan de zomerbedverbreding. De nevengeul is slechts bij hoge afvoeren stroomvoerend en heeft onder normale omstandigheden geen invloed op de morfologische ontwikkeling van het zomerbed. • Het eerste effect is de bodemstijging ten opzichte van de nulsituatie met zo’n 2 meter ter plekke van het traject waar de verbreding van het zomerbed is uitgevoerd. Het tweede belangrijke effect is het ontstaan van een erosiegolf net benedenstrooms van de zomerbedverbreding. Het dal van de golf ligt bij het ontstaan ervan zo’n 2 meter onder het bodemniveau van de nulsituatie en verplaatst zich in benedenstroomse richting. Tijdens de voortplanting van de golf verandert de vorm ervan. De golf wordt minder diep, maar wel breder. Na ongeveer 10 jaar bereikt de golf

56 Morfologische gevolgen van rivierverruimende maatregelen

de Beneden en Nieuwe Merwede. De golven hebben na 50 jaar nog niet het eind van beide takken bereikt. De verwachting is wel dat de effecten zich ook op de takken benedenstrooms zullen voordoen. • Naast de SOBEK berekeningen is tevens een analytisch onderzoek gedaan naar de effecten van de maatregelen. De overeenkomsten tussen deze analyse en de SOBEK resultaten zijn goed.

6.2 Aanbevelingen

Enkele aanbevelingen met betrekking tot de morfologische calibratie van het model: • De problemen met de morfologische modellering aan het eind van de Nieuwe Merwede door de aanwezigheid van slib, zouden aangepakt kunnen worden door het inbouwen van de mogelijkheid van zand-slib-segregatie in SOBEK. Tevens maakt dit een verdere uitbreiding van de calibratie naar de meer slibrijke takken mogelijk. • Ook zou een aanpassing in SOBEK betreffende de sedimentvoerende breedte gewenst zijn. Deze is nu beperkt tot de breedte van de hoofdgeul, waardoor de morfologische veranderingen bij hoge afvoeren op een relatief kleine breedte ten opzichte van de totale stroomvoerende breedte plaatsvinden. Dit levert onstabiele modelruns op, omdat de bodemstijgingen te grof zijn. • Uitbreiden van de Baseline-schematisaties naar de benedenstroomse takken, waardoor de stabiliteit en nauwkeurigheid van de modelberekeningen hier toenemen. • Calibratie verder uitwerken voor overige takken van het NDB. Tezamen met de uitbreiding van de Baseline-schematisatie zou het model verder gecalibreerd kunnen worden. Zolang er geen mogelijkheden voor modellering van zand-slib-segregatie ingebouwd zijn, is verdere uitbreiding van de morfologische calibratie alleen mogelijk op de overige zandige takken. Het bovenstroomse deel zou eerst afgewerkt kunnen worden door bijvoorbeeld de Maas en de Lek te calibreren.

Hieronder volgen de aanbevelingen naar aanleiding van het onderzoek naar de effecten van de rivierverruimende maatregelen: • De maatregelen in de Boven Merwede bij Gorinchem, zorgen voor een dusdanige bodemstijging dat het primaire doel van deze maatregelen – waterstandverlaging bij extreme afvoeren – deels teniet gedaan zou kunnen worden, wanneer de baggerwerkzaamheden op dit traject niet worden geïntensiveerd. Aanbevolen wordt om maatregelen te treffen die geen effect hebben tijdens normale omstandigheden. De nevengeul is hier een voorbeeld van, maar ook de maatregel ‘Steurgat’. Door het inzetten van deze maatregelen alleen tijdens extreme afvoeren blijven de morfologische effecten die deze maatregelen met zich meebrengen zeer beperkt. • Gedetailleerd gevoeligheidsonderzoek uitvoeren: er zouden enkele variaties op de maatregelen doorgerekend kunnen worden om zo de gevoeligheid ervan op de morfologische effecten te achterhalen.

57 Morfologische gevolgen van rivierverruimende maatregelen

Hierbij valt te denken aan een andere drempelhoogte voor het Steurgat en voor de nevengeul bij Gorinchem, maar ook bijvoorbeeld variaties in de zomerbedverbreding. Er zijn al varianten beschreven op de in dit hoofdstuk beschreven maatregelen. De morfologische effecten van deze varianten zouden ook met het model beschouwd kunnen worden. • Naast de in dit rapport beschreven maatregelen liggen er ook nog andere plannen op tafel voor de Merwede takken in het kader van ‘ruimte voor de rivier’. Het model zou ook in deze gevallen gebruikt kunnen worden om de morfologische effecten ervan te onderzoeken. Bovendien kunnen bij uitbreiding van de morfologische calibratie maatregelen in andere takken onderzocht worden.

58 Morfologische gevolgen van rivierverruimende maatregelen

Referenties

......

Literatuur Allersma, E., 1988, Morfologisch Onderzoek Noordelijk Deltabekken, Morfologische Modellering, deel 1t/m 4. WL | Delft Hydraulics, Z71- 03.

Bouwdienst Rijkswaterstaat, in opdracht van de projectorganisatie Ruimte voor de Rivier, 2002, Startnotitie MER, in het kader van de PKB-procedure Ruimte voor de Rivier, Min. V&W, Min. VROM en Min Landb., Nat. en Vis.

Brinke, W.B.M. ten, L.J. Bolwidt, E. Snippen, L.W.J. van Hal, 2001, Sedimentbalans Rijntakken 2000, Een actualisatie van de Sedimentbalans voor Slib, Zand en Grind van de rijntakken in het Beheersgebied van de Directie Oost-Nederland, RIZA, rapport 2001.043.

Bureau Benedenrivieren, 2002, Maatwerk, Overzicht van Maatregelen, Ontwerpproces met de Regio – Benedenrivieren, uitgave in het kader van ‘Ruimte voor de Rivier’.

Dreumel, P.F. van, 1995, Slib- en Zandbeweging in het Noordelijk Deltabekken in de periode 1982 –1992, Rijkswaterstaat Directie Zuid- Holland, Afdeling Watersysteemkennis.

Fioole, A., B. Steenkamp, 2000, Effect Hoge Afvoeren op de Bodemontwikkeling van de Nieuwe Merwede, RIZA-WST, werkdocument 2000.200X.

Hartman, M.R., A.P.P. Termes, J. Udo, 2001, SOBEK-model Noordelijk Deltabekken, Nieuwe Merwede, Beneden Merwede en de Brabantse Biesbosch, HKV PR499.

Houwing, E.J., G. Blom, 2001, Sedimenttransport en Morfologische Ontwikkeling van het Noordelijk Deltabekken, RIZA-WST, rapport 2001.021.

Jesse, P., D.F. Kroekenstoel, 2001, 1-D Morfologische model Rijntakken, Verslag Calibratie en Validatie, RIZA, rapport 2001.040.

KNMI, 2003, De Toestand van het Klimaat in Nederland 2003.

Koomans, R.L., 2003, Verontreiniging van de Waterbodem, Amer, Merwede, Bergse Maas, Medusa Explorations, project 2003-P-027.

59 Morfologische gevolgen van rivierverruimende maatregelen

Ledden, M. van, 1999, Zand-Slib-Segregatie, Modellering van Zand- slibstructuren in de Nieuwe Merwede, RIZA, werkdocument 98.155X

Ledden, M. van, 2003, Sand-Mud Segragation in Estuaries and Tidal Basins, proefschrift, TUD, ISBN 90-9016786-2.

Linden, M. van der, J.W. van Zetten, 2002, Een SOBEK-model van het Noordelijk Deltabekken: Bouw, Calibratie en Verificatie, RIZA, rapport 2002.002

Linden, M van der, 2003, Memo RIZA Maatregelen Concept 27-8- 2003 (in voorbereiding), RIZA Memo.

Lodder, Q., Q. Gao, 2003, Concept Rapportage verdiepingsslag BER Gorinchem, RIZA Memo.

Ogink, H.J.M, J.G. Grijsen, 1986, Hydraulische ruwheid van de Merweden en de Maas, Verslag van Onderzoek, WL | Delft Hydraulics rapport R2393/Q388.

Ogink, H.J.M., 2002, Quick Scan Zomerbedverdieping Merweden, WL | Delft Hydraulics, notitie Q3304.

Ribberink, J.S., R.H. Buijsrogge, 2002, Transportverschijnselen en Morfologie, Civiele Techniek Universiteit Twente, dictaat nr. 881.

Ribberink, J.S., S.J.H.M. Hulscher, R.H. Buijsrogge, 2002, Ondiepwaterstromingen, Civiele Techniek Universiteit Twente, dictaat nr. 866.

Rijkswaterstaat, 2001, Beschrijving Functies Boven Merwede, Beneden Merwede en Nieuwe Merwede, BPN versie januari 2001.

SOBEK Help Desk, 2001, SOBEK River/Estuary User Manual, WL | Delft Hydraulics, afkomstig van Sobekversie 2.51.001 cd

SOBEK Help Desk, 2000, SOBEK River/Estuary Technical Reference Manual, WL | Delft Hydraulics.

Wang, Z.B., R.J. Fokkink, M. de Vries, A. Langerak, 1995, Stability of River Bifurcations in 1D Morphodynamic Models, J. of Hydr. Res., Vol. 33, No. 6.

Wang, Z.B., A.W. Weck, E.J. van der Houwing, 2000, Morfologie van het Noordelijk Deltabekken, geschiktheid van numerieke modellen voor voorspelling van morfologische ontwikkelingen, WL | Delft Hydraulics, rapport Z2662.

Wang, Z.B., 2000, Morfologie van het Noordelijk Deltabekken, Pilot toepassing van 1D netwerk modellen SOBEK-MOR en ESTMORF, WL | Delft Hydraulics, rapport Z2662.

60 Morfologische gevolgen van rivierverruimende maatregelen

Wang, Z.B., 2001, Een 1D Netwerk Morfodynamisch Model voor het Noordelijk Deltabekken, Opzet van het Model, WL | Delft Hydraulics, rapport Z3146.

Werf, J.J. van der, 2001, Graded Sediment at River Bifurcations, an Analytical and Numerical Study of the Rhine Bifurcations, Universiteit Twente, Document 2001.147x.

WL | Delft Hydraulics en HKVlijn in water, 2003, Rivierkundige analyses maatregelen Benedenrivierengebied ten behoeve van PKB Ruimte voor de Rivier, Rapport Q3244.29.

WL | Delft Hydraulics en HKVlijn in water, 2003, PKB Ruimte voor de Rivier – Hydraulische effectbepaling – Aanvullende werkzaamheden, Rapport Q3244.56.

Gegevensbronnen DONAR database Rijkswaterstaat.

Fugro Ingenieursbureau, 2002, Laboratoriumresultaten betreffende Morf. Modelleren (Nr. 610106004), opdrachtnummer H-4086.

RWS-DZH Rotterdam, 2002, Evaluatie Baggerwerk 1982-2001 voor Merwede.

61 Morfologische gevolgen van rivierverruimende maatregelen

62 Morfologische gevolgen van rivierverruimende maatregelen

Bijlagen

......

Bijlage I: modellering waterbeweging SOBEK Bijlage II: SOBEK schematisatie takken NDB Bijlage III: beschrijving Baseline Bijlage IV: hydraulische calibratie model Bijlage V: korreldiameters NDB Bijlage VI: baggergegevens onderzoekstakken Bijlage VII: bodemliggingen onderzoekstakken Bijlage VIII: aanpassingen korreldiameter Bijlage IX: controle hydraulische calibratie Bijlage X: korte termijn effecten maatregelen Bijlage XI: bodemontwikkeling nulsituatie

63 Morfologische gevolgen van rivierverruimende maatregelen

64 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage I Modellering waterbeweging SOBEK

......

De stroming wordt in SOBEK gemodelleerd door de continuïteitsvergelijking voor water en de bewegingsvergelijking voor water. De continuïteitsvergelijking luidt:

∂A ∂Q t +=q ∂∂txlat

De bewegingsvergelijking luidt: ∂∂QQ⎛⎞2 ∂hgQ Q τ g∂ρ ++αηgA +−W wi +gA ()+ξQ Q +A =0 ⎜⎟Bf2 ff lm ∂∂tx A ∂xCRA ρρ∂x ⎝⎠ffw w

In deze vergelijkingen zijn de volgende symbolen gebruikt: Q debiet [m3/s] t tijd [s] x afstand [m]

αB constante van Boussinesq [-] 2 Af stroomvoerende dwarsdoorsnede [m ] 2 qlat lateraal debiet [m /s] g zwaarteversnelling [m/s2] h waterstand (t.o.v. referentieniveau) [m] C Chézy coëfficiënt [m1/2/s] R hydraulische straal [m]

Wf stroomvoerende breedte [m] 2 τwi schuifspanning t.g.v. wind [N/m ] 3 ρw dichtheid water [kg/m ] η eerste additionele weerstandscoëfficiënt [-] ξ tweede additionele weerstandscoëfficiënt [-] 3 Alm eerste orde moment dwarsdoorsnede [m ]

65 Morfologische gevolgen van rivierverruimende maatregelen

66 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage II SOBEK schematisatie takken NDB

......

SOBEK vakken in het oorspronkelijke model.

Bovenstaande figuur geeft een goede indruk van de vakken schematisatie in het SOBEK NDB model. Dit model is echter in deze studie niet gebruikt. In het basismodel van deze studie zijn de takken van de meest bovenstroomse rivieren al volgens Baseline geschematiseerd, wat hier nog niet gebeurd is. Bovendien is tijdens deze studie de Baseline schematisatie uitgebreid naar het gebied rondom de Biesbosch (paragraaf 3.3).

Takken en knopen in het huidige model:

In bovenstaande figuur is het netwerk van takken en knopen weergegeven. Ook zijn de belangrijke kunstwerken in het gebied aangegeven, die als ‘structures’ zijn geschematiseerd in SOBEK.

67 Morfologische gevolgen van rivierverruimende maatregelen

68 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage III Beschrijving Baseline

......

Voor de aanmaak van de vernieuwde dwarsprofielen voor het SOBEK- model is Baseline gebruikt. Baseline is een softwarepakket waarmee schematisaties ten behoeve van de invoerbestanden voor SOBEK gestructureerd kunnen worden uitgevoerd. Baseline bevat een ARC/INFO-database waarmee geografische informatie van een rivierengebied wordt verwerkt tot een netwerkstructuur met doorstroomprofielen, bergend oppervlak, zomerkaden, enz. Om SOBEK-dwarsprofielen te genereren dienen de Baseline basisbestanden aanwezig te zijn. De aanmaak van deze bestanden is als volgt: 1. Bepalen vakgrenzen; een tak binnen het SOBEK-model is opgedeeld in SOBEK-vakken. Deze vakken hebben een lengte van 500 m. De eindvakken van een tak hebben een halve vaklengte (250 m). De topografische gegevens binnen een SOBEK-vak worden in Baseline gemiddeld tot één representatief dwarsprofiel. Het dwarsprofiel wordt in het midden van het vak gelegd. Dit is tevens de locatie van het rekenpunt. 2. Begrenzing zomerbed, oeversectie en winterbed; een SOBEK- vak bestaat uit maximaal drie secties: het zomerbed, de oeversectie en het winterbed. De secties worden begrensd door de normaallijn (grens tussen het zomerbed en de oeversectie), de oeverlijn (grens oeversectie en het winterbed) en de bandijk (buitenste begrenzing van het winterbed). Wanneer de oeversectie is bepaald, liggen de langsbegrenzingen van het zomerbed en winterbed vast. De oeversectie wordt begrensd door de kribkoppen (normaallijn) en de locaties waar de kribruggen in de bodem verdwijnen (oeverlijn). Dit kan bepaald worden met behulp van topografische kaarten. 3. Begrenzing stroomvoering en berging; de gebruikelijke werkwijze is dat met behulp van een WAQUA-model (2D- stromings model) de stromingspatronen in het zomer- en winterbed worden berekend voor een aantal stromingsomstandigheden. Daarmee kan de grens tussen het stroomvoerende deel en het bergende deel van het winterbed worden bepaald. Het onderscheid tussen beide delen wordt gedaan op basis van een stroomsnelheidscriterium. Gedeelten met zeer lage stroomsnelheden worden aan het bergende deel toegekend, de overige delen van het winterbed aan het stroomvoerende deel. Bij het ontbreken van een actueel WAQUA-model kan met behulp van topografische kaarten en ‘expert judgement’ de grens tussen het stroomvoerende en bergende deel van het winterbed langs de waterlopen worden bepaald. 4. Kades en plassen; binnen SOBEK wordt een zomerkade gebruikt om aan te geven vanaf welke waterstand een bergend dan wel een stroomvoerend deel van het dwarsprofiel in de

69 Morfologische gevolgen van rivierverruimende maatregelen

bewegingsvergelijkingen wordt meegenomen. Het kadeniveau kan de hoogte zijn wanneer een omkaad gebied inundeert of een plas gaat meestromen. De zomerkade in Baseline is een standaardniveau binnen een dwarsprofiel. De ligging van zomerkades langs de waterlopen en omkade gebieden in het winterbed wordt vanaf topografische kaarten bepaald. Aangegeven moet worden wat de instroomhoogte is van een omkaad gebied. Het laagst gelegen punt van het omkade gebied wordt als instroomhoogte genomen. Het Baseline- programma berekent vervolgens het bergende oppervlak achter deze kades per profiel met behulp van de zomerkadehoogte (instroomhoogte). Veerstoepen en wegen op kades die min of meer dwars op de stroming lopen worden in de schematisatie van de dwarsprofielen meegenomen. In het ene geval kunnen ze als overlaat worden geschematiseerd om daarmee een lokaal energieverlies te simuleren. In het andere geval worden ze als lokale verhoging in het terrein meegenomen. Per voorkomende situatie wordt het verloop van dit type kade bekeken en geschematiseerd. Plassen in het winterbed kunnen aangetakt of niet-aangetakt zijn. Aangetakte plassen liggen in het winterbed en staan in open verbinding met het zomerbed (bijvoorbeeld havens). Niet- aangetakte plassen liggen in het winterbed en staan niet in verbinding met het zomerbed. Beide type plassen kunnen (geheel of gedeeltelijk) in het stroomvoerende deel van het winterbed liggen of in het bergende deel. En beide typen plassen kunnen dus ook (gedeeltelijk) stroomvoerend zijn tijdens hoogwater mits het winterbed meestroomt.

Het converteren van de GIS-gegevens naar een SOBEK-dwarsprofiel met behulp van de Baseline gaat in twee stappen: 1. Het omzetten van de basisbestanden naar een grid. 2. Het omzetten van de grids naar dwarsprofielen.

70 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage IV Hydraulische calibratie model

......

Aangezien de Merwede takken aangepast zijn in het nieuwe model, moet het model voor deze takken opnieuw gecalibreerd worden. De calibratie is uitgevoerd is samenwerking met J.W. van Zetten (RIZA). De toegepaste calibratie-methodiek is gelijk aan die van voorgaande calibraties. De calibratie geschiedt aan de hand van metingen van 5 meetstations gedurende de hoogwatergolf in de periode 24 oktober tot 15 november 1998. De 5 meetstations zijn: Werkendam, Vuren, Deneplaat, Moerdijk en Dordrecht (tabel IV.1). Als calibratievariabele wordt de ruwheid gebruikt, die in SOBEK is gegeven door middel van een waterstandafhankelijke Chézy waarde voor de Nieuwe Merwede en een Manning coëfficiënt voor de Beneden Merwede. Doel van deze calibratie is dezelfde of betere resultaten te realiseren als tijdens de calibratie van het vorige hydraulische model [van der Linden en van Zetten, 2000]. Het vorige hydraulische model is gecalibreerd met de doelstellingen om de waterstanden die het model voorspelt gemiddeld niet meer te laten afwijken dan 5 cm van de gemeten waarden met een standaardafwijking van maximaal 5 cm. Deze doelstellingen zijn tijdens deze eerdere calibratie niet in alle gevallen gehaald.

Meetstation Riviertak Kilometerraai SOBEK locatie Vuren Boven 952 Waal_2_37944 Merwede Werkendam Nieuwe 962,5 Nime_1837 Merwede Deneplaat Nieuwe 977,2 Nime_16168,9 Merwede Moerdijk Hollandsch 985 Hodi055_722,5 Diep Dordrecht Beneden 976,5 Beme_15119 Merwede Tabel IV.1: gebruikte meetstations tijdens hydraulische calibratie.

De volgende stappen zijn gehanteerd tijdens het afregelen van het model: 1. Initiële run: zonder aanpassingen is het model gedraaid en zijn de metingen met de SOBEK-resultaten vergeleken. 2. Op basis van de resultaten hiervan, zijn een aantal alternatieven opgesteld, waarin op verschillende wijze de Chézy-waarden zijn aangepast. 3. De resultaten van de alternatieven zijn vergeleken en op basis hiervan is weer een nieuwe set alternatieven opgesteld. 4. Stap 2 en 3 zijn een aantal keer herhaald, totdat er geen verbeteringen meer zijn opgetreden. 5. Het alternatief met de beste overeenkomst met de gemeten waarden is uiteindelijk gekozen.

71 Morfologische gevolgen van rivierverruimende maatregelen

De gemiddelde waterstandafwijkingen en standaarddeviaties voor de verschillende locaties voldoen niet allen aan de eisen (tabel IV.2). De gemiddelde waterstandafwijkingen zijn op Vuren na wel allen kleiner dan 0,05. Een verklaring hiervoor zou kunnen zijn dat bij de calibratie aan de randen gewerkt wordt met debieten, die met behulp van zogenaamde Q-h-relaties zijn afgeleid uit gemeten waterstanden. Deze Q-h-relaties blijken vrij onnauwkeurig te zijn, waardoor de afregeling van het model wordt bemoeilijkt. De onnauwkeurigheid wordt onder andere veroorzaakt doordat de gebruikte relaties verouderd zijn. Ook wordt in de relaties nog geen rekening gehouden met het zogenaamde hysterese verschijnsel: tijdens toenemende afvoeren geldt een andere Q-h-relatie dan bij afnemende afvoeren. De standaardafwijkingen zijn allen groter dan 0,05. Bij de calibratie van het voorgaande model is dit ook geconstateerd. De oorzaak hiervan ligt in het optreden van een faseverschil van de getijgolf in SOBEK ten opzichte van de werkelijkheid; doordat de voortplantingssnelheid van de getijgolf niet exact overeenkomt met de werkelijke waarden zullen op locaties verder van de zeerand afwijkingen ontstaan in de waterstanden.

Locatie Werkendam Vuren Dordrecht Deneplaat Moerdijk Gemiddelde afwijking [m] -0,01 0,06 0,03 0,02 0,03 Standaarddeviatie [m] 0,06 0,07 0,06 0,06 0,07 Tabel IV.2: resultaten hydraulische calibratie

De effecten op morfologie van bovenbeschreven faseverschil zijn verwaarloosbaar. Een klein faseverschil in de getijgolf zal geen drastische verandering teweeg brengen in de door het getij veroorzaakte component van de stroomsnelheid. Hiermee zal het effect op sedimenttransport en morfologische ontwikkeling ook nihil zijn. Dat een klein faseverschil al wel behoorlijke standaarddeviaties veroorzaakt blijkt uit figuur IV.1. Hierin zijn twee identieke golven weergegeven met een klein faseverschil ten opzichte van elkaar. De gemiddelde waterstanden van beide golven zijn gelijk en dus is de gemiddelde afwijking gelijk aan nul, doch door de structurele afwijking in waterstand, worden relatief hoge standaarddeviaties veroorzaakt.

Golf 1 Golf 2 (faseverschil +10%) Verschil

Figuur IV.1: verschil tussen twee (sinus) golven met een faseverschil van 10%.

72 Morfologische gevolgen van rivierverruimende maatregelen

In vergelijking met het eerdere hydraulische model zijn de resultaten niet verslechterd. Het gebruik van dit model als basis voor het morfologische model is hiermee gerechtvaardigd en wenselijk, gezien de verbeterde schematisering van enkele belangrijke takken.

73 Morfologische gevolgen van rivierverruimende maatregelen

74 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage V Korreldiameters NDB

......

In onderstaande figuren zijn de korrelgegevens weergegeven voor de takken van het NDB die zijn ingevoerd in het model. De gegevens voor de onderzoekstakken zijn hier niet opgenomen, deze zijn reeds in hoofdstuk 3 getoond. De figuren zijn gemaakt op basis van de gegevens van Fugro [2002].

Maas - Bergsche Maas - Amer

1 ruw e data 5 km gemiddeld 0,8

) 0,6 m m

50 ( 0,4 D

0,2

0 200 210 220 230 240 250 260 270 kmr

Hollandsch Diep - Haringvliet

0,5 ruw e data 5 km gemiddeld 0,4

) 0,3 m m

50 ( 0,2 D

0,1

0 980 985 990 995 1000 1005 1010 1015 1020 1025 1030 kmr

75 Morfologische gevolgen van rivierverruimende maatregelen

Dor dts che Kil

0,6 ruw e data 0,5 5 km gemiddeld

0,4

mm) 0,3 ( 0 5 D 0,2

0,1

0 983 984 985 986 987 988 989 990 kmr

Noor d

0,5

0,4

) 0,3 m m

50 ( 0,2 D

0,1 ruw e data 5 km gemiddeld 0 976 977 978 979 980 981 982 983 984 985 kmr

Lek - Nieuwe Maas - Nieuwe Waterweg

1,4 ruw e data 1,2 5 km gemiddeld

1 )

m 0,8 m

0 ( 0,6 D5 0,4

0,2

0 950 960 970 980 990 1000 1010 1020 1030 kmr

76 Morfologische gevolgen van rivierverruimende maatregelen

Oude Maas

0,6 ruw e data 0,5 5 km gemiddeld

0,4 ) m

m 0,3 50 ( D 0,2

0,1

0 975 980 985 990 995 1000 1005 kmr

Spui 0,5 ruw e data 5 km gemiddeld 0,4

0,3 mm) ( 0 5 0,2 D

0,1

0 995 998 1001 1004 1007 1010 1013 kmr

77 Morfologische gevolgen van rivierverruimende maatregelen

78 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage VI Baggergegevens onderzoekstakken

......

Baggerwerkzaamheden op de Waal per vak (vak 1 = kmr 913-922, vak 2 = kmr 923-932, vak 3 = kmr 933-942, vak 4 = kmr 943-952):

WAAL evaluatie baggerwerk per jaar tijdvak 1990- 1996 600000

480000

360000

240000 gebaggerde hoeveelheden in m3 120000

0 1990 1991 1992 1993 1994 1995 1996 jaar

vak 1 vak 2 vak 3 vak 4 Totaal

Baggerwerkzaamheden op de Boven Merwede per vak (vak 1 = kmr 952.5-957.5, vak 2 = kmr 957.5 - 961.3):

BOVENMERWEDE evaluatie baggerwerk per jaar tijdvak 1990 - 1999 400000

350000

300000

250000

200000

150000

gebaggerde hoeveelheden in m3 100000

50000

0 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 jaar

vak 1 vak 2 totaal

79 Morfologische gevolgen van rivierverruimende maatregelen

Baggerwerkzaamheden op de Nieuwe Merwede per vak (vak 1 = kmr 961.3-964.2, vak 2 = kmr 964.2-968.5, vak 3 = kmr 968.5-979.8):

NIEUWE MERWEDE evaluatie baggerwerk per jaar tijdvak 1990 - 1999

200000

175000

150000

125000

100000

75000

50000 gebaggerde hoeveelheden in m3

25000

0 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 jaar

vak 1 vak 2 vak 3 totaal

Baggerwerkzaamheden op de Beneden Merwede per vak (vak 1 = kmr 961.3-968.3, vak 2 = kmr 968.3 – 976.0):

BENEDEN MERWEDE evaluatie baggerwerk per jaar 1990 - 1999

125000

100000

75000

50000 gebaggerde hoeveelheden in m3 25000

0 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 jaar

vak 1 vak 2 totaal

80 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage VII Bodemliggingen onderzoekstakken

......

De bodemliggingen die in deze bijlage zijn getoond hebben allen betrekking op de gemiddelde bodemhoogte van de hoofdgeul.

Gemeten bodemligging in 1990 en 2000 voor de Waal:

Waal

-1,00 Bodemligging 2000

) Bodemligging 1990 m -2,00 (

P A

N -3,00 ng tov -4,00 ggi i l

-5,00 bodem

-6,00 915 920 925 930 935 940 945 950 kmr

Gemeten bodemligging in 1990 en 2000 voor de Boven Merwede:

Bove n M e r w e de

-4,00

) -4,50 m (

P

A -5,00 N v o -5,50 ng t ggi i -6,00

Bodemligging 2000 bodeml -6,50 Bodemligging 1990 -7,00 952 954 956 958 960 kmr

81 Morfologische gevolgen van rivierverruimende maatregelen

Gemeten bodemligging in 1990 en 2000 voor de Beneden Merwede:

Beneden Merwede

-4,00 Bodemligging 2000 )

m Bodemligging 1990 ( -5,00 P A N

-6,00 ng tov ggi i

-7,00 bodeml

-8,00 961 963 965 967 969 971 973 975 977 kmr

Gemeten bodemligging in 1990 en 2000 voor de Nieuwe Merwede:

Nieuwe Merwede

-3,00 ) m ( -3,50 P A N v o -4,00 ng t ggi i l

-4,50 bodem Bodemligging 2000 Bodemligging 1990 -5,00 961 963 965 967 969 971 973 975 977 979 981 kmr

82 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage VIII Aanpassingen korreldiameter

......

In onderstaande figuren zijn de aanpassingen in korreldiameter weergegeven voor de verschillende onderzoekstakken.

Waal en Boven Merw ede 2,50 5 km gemiddeld aangepaste korrelset 2,00 ruw e data )

m 1,50 m ( 0 1,00 D5

0,50

0,00 0 5 10 15 20 25 30 35 40 45 50 Afstand v.a. begin (km)

Nieuwe Merwede 0,60 5 km gemiddeld

0,50 aangepaste korrelset ruw e data 0,40 ) m m

( 0,30 0 D5 0,20

0,10

0,00 0 2 4 6 8 101214161820 Afstand v.a. begin (km)

Beneden Merwede 1,00 5 km gemiddeld aangepast korrelset 0,80 ruw e data )

m 0,60 m ( 0 0,40 D5

0,20

0,00 0 2 4 6 8 10121416 Afstand v.a. begin (km)

83 Morfologische gevolgen van rivierverruimende maatregelen

84 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage IX Controle hydraulische calibratie

......

In de bovenste tabel staan de resultaten van de hydraulische calibratie, beschreven in bijlage IV. Het hydraulische model is opnieuw gerund waarin de Chézy waarden zijn aangepast zoals in hoofdstuk 5 beschreven is. De resultaten hiervan staan in de twee onderste tabellen. Gemiddeld gezien levert een verandering van Chézy waarden op alleen de Beneden Merwede geen verslechtering op in het hydraulische model. Wanneer ook nog de Chézy waarden op de Nieuwe Merweden wordt gewijzigd, verslechtert het model wel.

Oorspronkelijk Chézy waarden (resultaten hydraulische calibratie, bijlage IV)

Locatie Werkendam Vuren Dordrecht Deneplaat Moerdijk Gemiddelde afwijking [m] -0,01 0,06 0,03 0,02 0,03 Standaarddeviatie [m] 0,06 0,07 0,06 0,06 0,07

Chézy Beneden Merwede aangepast

Locatie Werkendam Vuren Dordrecht Deneplaat Moerdijk Gemiddelde afwijking [m] 0,00 0,07 0,02 0,03 0,03 Standaarddeviatie [m] 0,05 0,06 0,06 0,06 0,07

Chézy Beneden en Nieuwe Merwede aangepast

Locatie Werkendam Vuren Dordrecht Deneplaat Moerdijk Gemiddelde afwijking [m] 0,03 0,08 0,02 0,04 0,03 Standaarddeviatie [m] 0,05 0,07 0,06 0,07 0,07

85 Morfologische gevolgen van rivierverruimende maatregelen

86 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage X Korte termijn effecten maatregelen

......

Gebruikte formules: Bepalen evenwichtswaterstanden na ingreep: 2 ⎛⎞Q 3 h = ⎜⎟ (X.1) ⎝⎠BCi Bepalen verloop stroomsnelheid en sedimenttransport met: Q u = (X.2) Bh sm= u5 (X.3) Hieruit volgt de gradiënt in sedimenttransport en de initiële bodemveranderingen met: ∂z ∂S b ∼ − (X.4) ∂∂tx

In bovenstaande formules zijn de volgende symbolen gebruikt: h waterdiepte [m] u stroomsnelheid [m2/s] Q debiet [m3/s] i bodemverhang [-] B breedte rivier [m] C Chézy coëfficiënt [m1/2/s] S sedimenttransport [m3/s]

zb bodemligging [m]

Steurgat: Door de debietonttrekking ter plaatse van de inlaat naar het Steurgat ontstaat een discontinuïteit in het verloop van het debiet. Volgens vergelijking X.1 daalt de evenwichtswaterstand benedenstrooms van de inlaat hierdoor. Bovenstrooms van het Steurgat is de evenwichtswaterstand (stippellijn) onveranderd, het debiet is hier immers niet gewijzigd. Het verschil in waterstanden wordt overbrugd middels een verhanglijn bovenstrooms van de ingreep (figuur X.1.a). Met behulp van vergelijking X.2 en X.3 is het nieuwe verloop van u en s bepaald (figuur X.1.b). Vervolgens is de transportgradiënt (figuur X.1.c) bepaald en de initiële bodemontwikkeling (figuur X.1.d). Uit figuur X.1 blijkt dat door de debietonttrekking de bodem bovenstrooms van de onttrekking erodeert. Ongetwijfeld zal deze erosie ook invloed hebben op de Boven Merwede (het Steurgat ligt niet ver benedenstrooms van de splitsing Merweden) en op de Beneden Merwede. Bovendien verandert de waterverdeling tussen de Nieuwe en Beneden Merwede, waardoor er ook morfologische veranderingen in de Beneden Merwede plaatsvinden. Ter plekke van de onttrekking zelf ontstaat echter een geconcentreerde sedimentatieplek. Aangezien het hier om een initiële verandering gaat zal de sedimentatie op dit punt niet eeuwig doorgaan. Deze

87 Morfologische gevolgen van rivierverruimende maatregelen

sedimentatieplek zal zich als een bodemgolf in benedenstroomse richting gaan verplaatsen. Overigens is bij deze redenering geen rekening gehouden met het feit dat er naast een wateronttrekking ook een sedimentonttrekking plaatsvindt, waardoor de voorspelde bodemstijging ter hoogte van de inlaat lager uitvalt. Door de aanwezigheid van een drempel ter plekke van de inlaat naar het Steurgat zal echter het bodemtransport nauwelijks wijzigen.

a: waterstandsverloop h

Steurgat: -∆Q b: stroomsnelheid en sedimenttransport (u en s)

c: gradiënt sedimenttransport (ds/dx)

d: initiële bodemontwikkeling (dz/dt)

Figuur X.1: bepaling initiële bodemveranderingen na toepassing maatregel Steurgat.

88 Morfologische gevolgen van rivierverruimende maatregelen

Rivierverruiming knelpunt Gorinchem Vergelijking X.1 leert dat de evenwichtswaterstand op het traject waar de maatregelen gepland zijn daalt (figuur X.2.a), want door de nevengeul daalt het debiet door de hoofdgeul en verder zal de breedte toenemen. Hierdoor ontstaan er verhanglijnen in de waterspiegel. Met behulp van vergelijking X.2 en X.3 kunnen de bijbehorende stroomsnelheden bepaald worden evenals de sedimenttransporten (figuur X.2.b). Vervolgens worden de transportgradiënt (figuur X.2.c) en de initiële bodemveranderingen (figuur X.2.d) bepaald.

a: waterstandsverloop h

Traject maatregelen: -∆Q en +∆B b: stroomsnelheid en sedimenttransport (u en s)

c: gradiënt sedimenttransport (ds/dx)

d: initiële bodemontwikkeling (dz/dt)

Figuur X.2: bepaling initiële bodemveranderingen na toepassing maatregelen Gorinchem.

89 Morfologische gevolgen van rivierverruimende maatregelen

Figuur X.2.d laat zien dat benedenstrooms van het maatregelentraject geen initiële veranderingen in bodemligging optreden. Ter plaatse van de rivierverruimende maatregelen sedimenteert de hoofdgeul. Bovenstrooms van dit traject erodeert de bodem. Verder vindt er een geconcentreerde sedimentatie plaats juist bovenstrooms van het maatregeltraject en een geconcentreerde erosie op de benedenstroomse rand van het betreffende traject. Net als bij de maatregel Steurgat, planten deze geconcentreerde bodemveranderingen zich als bodemgolven voort in benedenstroomse richting.

90 Morfologische gevolgen van rivierverruimende maatregelen

Bijlage XI Bodemontwikkeling nulsituatie

......

De volgende figuren laten de bodemontwikkelingen zien in de nulsituatie volgens SOBEK. De veranderingen van de bodem na 10 jaar, na 25 jaar en na 50 jaar zijn hierin weergegeven. De uitgangssituatie is 1990.

bodemontwikkeling nulsituatie Waal

0,80

] 0,40 m g [ in 0,00 nder a r e

v -0,40

em na 50 jaar

bod -0,80 na 25 jaar na 10 jaar -1,20 913 918 923 928 933 938 943 948 953 kmr

bodemontwikkeling nulsituatie NIME

0,20

0,00 ] m

ng [ -0,20 i r

ande -0,40 r e v -0,60 dem na 50 jaar bo -0,80 na 25 jaar na 10 jaar -1,00 1 3 5 7 9 1 3 5 7 9 1 96 96 96 96 96 97 97 97 97 97 98 kmr

bodemontwikkeling nulsituatie BOME

0,80 na 50 jaar

] na 25 jaar m 0,40 na 10 jaar ng [ i

ander 0,00 r e v

dem -0,40 bo

-0,80 954 955 956 957 958 959 960 961 962 kmr

91 Morfologische gevolgen van rivierverruimende maatregelen