NMDC - ROOT of content - Wageningen UR

Daarnaast is veel tijd gaan zitten in de pre- en post processing, vooral omdat door de (zeer) fijne ruimtelijkeresolutie (zeer) veel data uit het modelinstrumentarium moesten worden verwerkt.ad 1 Koppeling MetaSWAP en WOFOST aan AMIGOVoor het NMDC-innovatie project bestond AMIGO uit het grondwatermodel MODFLOW voor de simulatie vangrondwater in de verzadigde zone en het hydrologisch model CAPSIM voor de onverzadigde zone. In decasestudie is CAPSIM vervangen door MetaSWAP. MetaSWAP berekent de grondwateraanvulling die alsbovenrandvoorwaarde wordt ingevoerd in MODFLOW, terwijl MODFLOW de grondwaterstand enkwel/wegzijgingsflux berekent die als onderrandvoorwaarde wordt ingevoerd in MetaSWAP. WOFOST berekentde gewasgroei onder invloed van straling, CO2-gehalte, temperatuur en vochtgehalten in de wortelzone. Ookwordt hierbij de diepte van de beworteling berekend. De gewasgroei wordt in de vorm van Leaf Area Index(LAI) doorgegeven aan MetaSWAP, die de evapotranspiratie vervolgens afhankelijk stelt van deze LAI. Decapillaire nalevering door zuigspanning in de wortelzone wordt in MetaSWAP afhankelijk gesteld van de metWOFOST berekende diepte van de wortels. Deze koppelingen zijn schematisch weergegeven in figuur 12.Figuur 12Uitwisseling van informatie tussen het verzadigde zone model MODFLOW, onverzadigde zone model MetaSWAP en hetgewasproductiemodel WOFOST.De koppelingen tussen MODFLOW-MetaSWAP-WOFOST zijn volledig geïntegreerd in het rekenproces; degekoppelde modellen wisselen geautomatiseerd informatie aan elkaar uit in iedere rekenstap of tijdstap. Er zijnook geen verschalingen, omdat ieder model in de keten met dezelfde ruimtelijke en temporele resolutie rekent.Voor de koppelingen met WOFOST gelden dezelfde (kritische) aandachtspunten als in de modelketen vanWalcheren. Een verschil hiermee is wel dat in MetaSWAP-WOFOST geen invloed van zout wordt meegenomen.Een belangrijke stap voorwaarts die binnen dit project gezet is, is het mogelijk maken van de MetaSWAPuitvoerin het zogenaamde IDF-formaat. Dit formaat is de standaard binnen de iMOD-omgeving en wordt bijDeltares gebruikt voor het opslaan en visualiseren van Modflow in- en uitvoer. IDF's zijn in feite niets andersdan binaire ascii-gridbestanden. Dat nu ook de MetaSWAP uitvoer in dit formaat beschikbaar is vergemakkelijktde visualisatie, verificatie en postprocessing van de modelresultaten. Een nog grotere, nog te zetten, stapvoorwaarts zou zijn als ook zoveel mogelijk MetaSWAP invoer in IDF-formaat aan de modelcode aangebodenzou kunnen worden. Dat maakt visualisatie van de invoer gemakkelijker en doordat vertaalslagen tussen deSVAT-units van MetaSWAP en de gridcellen van Modflow overbodig worden, de modelopzet minderfoutengevoelig.Zoals gezegd is het stroomgebied van de Baakse Beek voor de modellering onderverdeeld in vijftiendeelmodellen. Dit was noodzakelijk omdat het gehele model niet op de gebruikte rekenresolutie (te weten25x25 meter) gedraaid kan worden, zowel qua geheugengebruik als de doorlooptijd van de berekening. Devijftien deelmodellen kunnen parallel (tegelijkertijd) gedraaid worden, zodat er een enorme winst qua30 Alterra-rapport 2443