The Hollandsche IJssel plays an important role in the freshwater provision of the province Zuid-Holland. Consequently, for Rijkswaterstaat it is key that salt intrusion is minimal in the Hollandsche IJssel. Recent studies noted that salt intrusion in the Hollandsche IJssel is limited due to a phase difference between tidal velocities in the main channel, the Nieuwe Maas, and the side channel, the Hollandsche IJssel. Earlier research investigated the impact of phase differences between branches and found it can lead to increased dispersion in the main channel, through a process known as tidal trapping. At the same time, this phase difference can prevent the saltiest water from entering the side channel, as was found at the Hollandsche IJssel. Because of this role, it is relevant to find out how this phase difference may be influenced by sea level rise, more extreme river discharges and particularly how it depends on the geometry of the main and the side channels. Especially the latter could help Rijkswaterstaat to minimize salt intrusion at locations relevant to freshwater intake, such as the Hollandsche IJssel.

The main objective of this thesis is to investigate how the geometry of the side and main channel influences the tidal phase difference between these two channels, and how this may impact the salt dispersion in the side channel. For this, an analytical model is developed describing harmonic wave propagation in multi-branch systems and this is used next to results from a 3D numerical model for the Rhine Meuse Delta (RMM3D). First, the influence of changes in geometry and forcing is systematically investigated for a network containing a single junction. This shows that the length and depth of the side channel are the most significant variables. The depth is one of the main variables impacting friction, which governs the type of wave which can form in the system. A decrease in friction allows a wave to transform into a standing wave pattern as the return wave becomes more important, while increased friction transforms it into a propagating wave. The length also controls the type of wave which can form as it determines the distance along which the friction can work. Additionally, the length also governs potential resonance in the side channel.

Next, the phase differences of the M2, M4 and M6 tide are determined for the junction with the Hollandsche IJssel in the Rhine Meuse Delta (RMD) based on the RMM3D model. The main tidal constituent regarding tidal trapping was found to be M2. However, this does not fully represent the time difference between flow reversal at the Hollandsche IJssel and the Nieuwe Maas, which was found to be around 75 minutes. Additionally, the phase difference at the Lek was investigated. For the M2 tide at the Hollandsche IJssel and Lek, a phase difference of 55⁰ and 31⁰ was found, respectively. These phase differences prevent salt intrusion in the respective side channels. The inflow of the side channels starts while the main channel still flows to the sea during the ebb. At this moment, the salt concentrations in the main channel have already returned to background levels...

Due to climate change, low river discharges are expected more frequently. As a result, salt intrusion in the Rotterdam Waterway will be critical at least five times more frequently, negatively impacting freshwater availability. A possible new measure to reduce salt intrusion is providing additional mixing by internal wave generation, which reduces the estuarine circulation. The potential to reduce salt intrusion in the Rotterdam Waterway by internal wave generation is investigated. To determine this the irreversible mixing generated by internal waves is quantified and related to salt intrusion length. Numerical 2DV simulations of a decelerating flow over an undular bed were executed for an idealized section of the Rotterdam Waterway. The mixing was quantified by determining the background potential energy in a chosen control volume. This method can be applied if disturbances in the velocity or density profile do not reach the boundaries of the control volume. The mixing observed in the control volume is related to salt intrusion by the buoyancy frequency of the background state. The net salt flux imported by estuarine circulation was found to be proportional to the squared buoyancy frequency. Estuarine circulation is the main forcing mechanism of salt intrusion in the seaward part of the Rotterdam Waterway. For the small idealized section of the Rotterdam Waterway, internal waves generated over geometrical features are capable of reducing the salt intrusion caused by estuarine circulation by approximately 1.7% of the initial value, if bed waves are present over 27.5% of the control volume.

In this study the contribution of the low-frequency residuals to the sea-level variability has been examined. This is done using 106-year old sea-level record obtained at Hoek van Holland. Computing the mean sea-level per season each year and the corresponding standard deviation one finds an increase in both these features of the sea-level record. The rise of the mean sea-level implies the effect of climate change. Moreover, it is found that the standard deviation in the sea-level, thus the intensity of variation, is the highest during fall and winter. This implies that the sea-level variability has a seasonal dependency. Furthermore, one finds an increase in the standard deviation on the long term. However, since a big part of the mean sea-level is influenced by tidal events, the increase in standard deviation is possibly linked to climate change in meteorological factors as has been found is in the study carried out by Gerkema and Duran-Matute(2017). With the use of the computational algorithms such as the Fast Fourier Transformation and the Wavelet Transformation one can extract the low-frequency residuals from the sea-level record. Inspired by the study of Gerkema and Duran-Matute(2017) a correlation between the wind speed and the low-frequency residuals have been found. Using the Wavelet Transformation it is found that the low-frequency residuals obtain the most energy during fall and winter, this empowers the finding that these low-frequency residuals are seasonal dependent. Moreover, when studying the low-frequency residuals closely it is found that the frequencies below 0.60 1 day obtain the most energy during fall and winter, especially the frequencies near 0.10 1 day . With these findings one can say that the contribution of the low-frequency signals, thus the low-frequency residuals, is correlated to meteorological events, such as the wind. Moreover, it is found that the variation in the low-frequency residuals is much smaller during spring and summer compared to the case during winter and fall. This seems not to be the case for the sea-level variability due to the tidal constituents and the high-frequency waves. Therefore, one may conclude that these low-frequency residuals contribute to a great extent in the standard deviation of the sea-level.

In recent years, several logistic optimisation models have emerged as powerful tools to assess and optimise the planning, costs, and workability of marine operations. Nevertheless, these models often rely on two underlying assumptions: (1) the significant wave height and peak wave period adequately describe (directional) wave fields, and (2) these two parameters sufficiently describe the conditions causing weather downtime. However, little is known about how these underlying assumptions affect the reliability of logistic optimisation tools. This study, therefore, presents an alternative model that integrates response motions of vessels and turbine structures into a logistic optimisation tool to address and expose the implications that follow from using these assumptions. A case-study approach, on two recently realised offshore wind farm projects, showed that integrating response motions results in, up to 10%, less favourable workability conditions. Further analysis on the data showed that it is crucial to include the two-dimensional wave energy distribution to expose more complex sea states that induce weather downtime. Moreover, a failure analysis approach found that the conditions inducing downtime events are more accurately described by response motions instead of sea state parameters. Therefore, the findings of this study suggest that integrating response motions into logistic optimisation models improves the reliability of the model estimates. Besides, this study suggest that the industry’s approach potentially overestimates the true workability and, therefore, imposes unnecessary operational risks. Hence, the results of this study demonstrate the importance of integrating hydrodynamic engineering knowledge into the assessment and optimisation of project planning, costs, and workability.

Key points:- Pre-closure salinity intrusion into Haringvliet-Hollands Diep is known to have reached Biesbosch National Park at flood tide. Post-Delta21 salinity intrusion is projected to be less extensive than pre-closure. Maximum intrusion is estimated to reach the westernmost Moerdijk ports mainly due to diminished tidal flow at the estuary mouth.- Freshwater intake in the Haringvliet is projected to be compromised for the majority of the year whereas freshwater intake in Hollands Diep is compromised during prolonged drought.- Salinity outwash from the Haringvliet on ebb tide is projected to be poor under both drought and normal conditions due to widespread salinity diffusion in lateral and vertical directions on flood tide. The Delta21 framework, in which this study is positioned, aims at increased flood safety and ecological restoration of the Haringvliet, a former estuary in the Dutch southwestern delta that has been closed as part of the Delta works. To achieve this, Delta21 proposes to reopen the Haringvliet freshwater basin to tidal effects in an attempt to restore its estuarine character. Implementation of Delta21 introduces the risk of compromising agricultural and industrial activities around the Haringvliet-Hollands Diep by reintroducing salinity and tidal movement in the basin. The central problem lies in quantifying the extent and stability of the expected periodic salt intrusion post-Delta21. This has been done by projecting a geographical bandwidth of salinity intrusion patterns depending on the severity of SLR (2020-2100) and low-to-mean Rhine-Meuse discharges (in 2100 CE) using a numerical model. Mitigation strategies are then recommended based on these numerical results. The OSR-HV model (owner: Port of Rotterdam) is used to run predictive scenarios of salinity intrusion for the lower reaches of the Rhine-Meuse basin. OSR-HV runs in TRIWAQ (Rijkswaterstaat), which is 3D numerical modelling software that employs coupled hydrodynamics and constituent transport thereby resolving salinity transport. An upper-end critical scenario consists of a Rhine (Lobith) discharge averaging 1000 m3/s for 31 days combined with 85cm SLR in 2100. This resulted in an intrusion pattern reaching the westernmost port of Moerdijk in Hollands Diep. Results show that the basin geometry, possibly nudged by Coriolis deflection, initially causes a preferential path of salinity intrusion along the southern bank of the Haringvliet. Lateral and vertical mixing is extensive in the western part of Haringvliet which is thought to be a combined effect of weakened tidal flow at the estuary mouth and robust and erratic geometry of the basin. Further up-estuary, the historical flood-ebb tidal channel structures is the main transporter of salinity. The Haringvliet shows distinctly different estuarine behaviour compared to the neighbouring Rotterdam Waterways where stratification is more stable, causing less up-estuary diffusion. The relative robustness of the Haringvliet and mild freshwater forcing give rise to extensive 3D mixing which subsequently limits the maximum horizontal excursion of salinity. Salinity intrusion into the Old Meuse is observed to aggravate upon opening of the Haringvliet sluices due to flow reversal in Spui which connects Haringvliet with the Rotterdam Waterways. This effect can even result in salt intrusion from the Rotterdam Waterways via Old Meuse and Spui back into the Haringvliet, causing a secondary spike in salinity on ebb tide there. Finally, height-limitation of the Haringvliet sluice gates is somewhat effective in countering horizontal excursion of salinity but runaway diffusion in the Haringvliet results in similar salinity profiles compared to full opening of the sluice gates. A shipping channel that crosses the current Haringvliet front delta is included in Delta21. This deep feature promotes advection of salt through the Haringvliet sluices. Subsequent diffusion patterns cause for poor washout of salinity from the basin upon ebb tide. Limiting the depth of this channel is recommended if salt intrusion is to be reduced. Full opening of the Haringvliet sluices furthermore causes approx. 0.50m lowering of Mean Low Water (MLW) near Moerdijk which affects busy shipping routes between Rotterdam-Moerdijk-Scheldt. Partial reduction of the conveying area of the Haringvliet sluices may be used to suppress the tidal wave penetration into Haringvliet-Hollands Diep. Significant gain in ecological value is likely post-Delta21 due to addition of approx. 1900 ha of intertidal areas and a 40 km salinity gradient (excl. front delta). Ample recommendations on further research have been made in this exploratory study. It is recommended to further study the effects of Delta21 interventions on macro hydrodynamics of the Dutch coastal shelf. The interconnected nature of the region requires a larger modelling domain to prevent the occurrence of non-physical effects obtained from the current numerical schematization. Furthermore, hydrological relations and bathymetry were generated from 2020 data. It is therefore recommended to precede the assessment of salinity intrusion into the Haringvliet with numerical projections on change to these environmental factors. Lastly, significant gain in accuracy may be obtained from applying spatially varying temperature and wind to the domain to better replicate baroclinic flows and turbulent mixing.

The amount of suspended sediment spill released during dredging activities in the Fehmarnbelt strait is monitored along the excavated trench to prevent negative effects on local ecology. Stratification due to the interaction of saline water from the North Sea and fresh water from the Baltic Sea forces bi-directional flow of the layer above and below the density gradient, causing bi-directional spread of the dredging plume perpendicular to the trench as well. During such an event, monitoring on both sides of the trench is required. By mapping out the spatial and temporal behaviour of stratification, a prediction can be given on where and when measuring on both sides of the trench will be needed to include the baroclinic effect. Figures created with monitoring data show that the density profile is influenced by the salinity, rather than the temperature. Therefore, stratification predominantly depends on the in and outflow of water with respect to the Baltic Sea. The figures also show a larger density gradient during out than inflow. It can also be observed that where the water depth is restricted to 10 meters, wind and bottom friction mix the entire water column. Therefore, stratification occurs predominantly during outflow in sections deeper than 10 meters, indicating the need for monitoring the bi-directional plume spread during such circumstances. Whether stratification occurs during inflow in sections deeper than 10 meters likely depends on the duration and strength
of the wind forcing and the initial strength of the density gradient prior to the inflow event. Further analysis should be done to confirm this. Signs of Ekman transport, return flows and the deflection of the currents towards deeper water can also be observed in the measurement figures. Since these processes affect the plume spread direction, additional research can be done on the behaviour of the current direction in the Fehmarnbelt.

Two current assumptions of the Zeesluisformulering (ZSF) regarding a homogeneous lock density and instantaneous shipping are tested against refined descriptions according to literature: a stratified lock density and shipping according the continuity approach. These are implemented into a conceptual model constructed in Python, after which optimizations are quantified by means of a sensitivity analysis and a case study. It was found that by using the new assumption set, a maximum of approximately 10% less salt intrusion is obtained as a result of the influence of the return current on the lock exchange.

Due to climate change and human interventions, saltwater intrusion is becoming a topic of increasing concern worldwide. Salt water intrudes into the Rotterdam Waterway (RWW) by an exchange flow, where the denser sea water propagates landwards at the bottom. The main competing mechanism for this stratified exchange flow is vertical mixing, which can be realised by internal wave induced shear instabilities or wave breaking. The goal of this study is to investigate whether internal waves generated over undular bottom topography in the RWW can generate additional vertical mixing. The underlying assumption is that a decrease in stratification decreases salt intrusion.   The approach to answer the main research question is a combination of an analytical and a numerical analysis. The analytical study is based on frictionless linear theory. Internal wave behaviour is further analysed with FinLab, a finite element model which includes the non-hydrostatic processes and effects of density differences. FinLab is evaluated for the application of this study by means of a validation case.  In the analytical study, linear theory is applied to obtain a relation between the bed wave parameters and average internal wave energy density E for internal waves generated over sinusoidal bottom topography in a linearly stratified fluid. The derived expression describes that the bottom topography amplitude h0 and bed wave number kT both have a positive quadratic relation with the energy. Additionally, kTkinfluences the resonance conditions.  To validate FinLab for internal wave breaking and mixing an experiment in a wave tank, according to an example from literature, is simulated. The validation case reveals a shortcoming in the turbulent mixing parameterization. However, on scales relevant for the RWW the effect of this will not have the same significance. The validation case offers a suggestion for a subgrid closure of diffusion, where density effects are taken into account.  Numerical simulations of a 2D channel stretch with sinusoidal bottom topography, a linearly stratified fluid and a linearly varying background velocity, show generation of resonant trapped internal waves for the first two resonant modes. These occurrences correspond to the highest values of kinetic energy as function of vertical velocity averaged over the bed wave domain. The vertical buoyancy flux b is downward directed during occurrences of internal waves and becomes upward directed for increasing background flow. Vertical mixing is associated with an increase in average potential energy Ep, which is 17% higher for the base case (containing bed waves) than for a similar case without bed waves. This increase is larger when bottom shear stress increases. Richardson numbers below 0.25, associated with shear instabilities and mixing, are only observed near the bed, mainly when internal waves are present. The effect of variations in bottom topography wavelength LT and amplitude h0 on internal wave energy can be explained by the analytical formulation. The effect of bed wave parameter changes on b and relative increase in Ep can be related to the effect of the changed amount of bed friction rather than the difference in wave energy.   The first resonant mode is the most energetic, however, the average energy density found for these waves is only 0.4% to 6.7% of the potential energy anomaly (PEA); the energy required to fully mix a stratified water column. In the simulations the only mechanism that could transfer internal wave energy to turbulent kinetic energy are shear instabilities near the bed. Over the full simulation, the net vertical buoyancy transport is of negligible magnitude, where Ep shows significant increases between 6% and 99% compared to similar cases without bed waves and is enhanced during the presence of internal waves.  The main discussion point is that the quantification of vertical mixing requires improvement, particularly to determine the importance of mixing by internal wave-induced shear instabilities and by bed shear. Mixing by local shear instabilities (of which the relevant scales cannot be resolved with the current grid resolution) does not have an adequate parameterization, because density effects are not included in the turbulence closure. The bed friction parameter, which greatly influences the behaviour of the system, has to be validated. Furthermore, cases where internal waves might break in practice (e.g. at banks) were not considered. Finally, the observed internal wave energy is of small magnitude, however field measurements by Pietrzak(1991) shows that turbulence production by internal waves was significant.   

Phenomena like sea level rise, global warming and erosion together contribute to increasing flood risk in vulnerable coastal areas. As this flood risk increases, initiatives to mitigate the effects of climate change in the coastal zone are also increasingly sought. During recent years, the view that these measures should be circular has gained significant momentum. In line with this, the idea that sediment should not be treated as waste, but as a valuable product in a circular economy has been widely accepted. Within this theme, the sediment uses as resources in circular and in territorial economies (SURICATES) was created to develop and execute eco-innovative solutions for the reuse of sediments in Western Europe. As a part of SURICATES, over the course of a nine week pilot, a total of 500 tonnes of sediment was reallocated in a designated area in the Rotterdam Waterway, with the expectation that this sediment would be transported out of the Rhine-Meuse estuary into the North Sea. The SURICATES project is one of the first real large-scale efforts to reuse sediment for economic as well as ecosystem services. Effective implementation of the SURICATES project, however, requires a thorough understanding of the governing physical processes impacting the distribution of the deposited sediment locally. Furthermore, to assess the feasibility of similar future applications, efficient modelling of the pilot project is necessary. The work presented here aims to elucidate the hydrodynamic processes that are governing in the Rotterdam Waterway, within the framework of the SURICATES pilot project, and assess the reproducibility of these processes by (two) predictive models that are currently operational at the Port of Rotterdam. A special 6-hour monitoring survey was set up to measure salinities, velocities, temperature, and suspended particulate matter (SPM) along a transect crossing the reallocation location. In combination with a literature review, this dataset provides the basis for research into the predictive capabilities of two currently operational hydrodynamic models. Analysis of this dataset reveals the dominant terms in the momentum balance, the influence of Coriolis, the occurrence of internal waves, and the effect that all these mechanisms may have on the SPM distribution around the reallocation location. When the system dynamics are elucidated, the model performance of the two hydrodynamic models is assessed—both quantitatively and qualitatively. It is also investigated whether phase shifts are introduced in the models. It is found that the primary hydrodynamic processes in the Rotterdam Waterway are related to the barotropic tidal asymmetry imposed at the river mouth, the tidal excursion of the salt wedge, baroclinic exchange flow processes, and turbulence damping at the pycnocline. Turbulence damping at the pycnocline generally poses an upper limit to the (re)distribution of SPM over the water column, although field data suggests that this damping may not be sufficient to counteract diffusion processes locally. This effect occurs under certain forcing conditions and during low water slack. Furthermore, an internal Froude number analysis provides evidence for the possible generation of internal waves in one of the river’s bends. The evaluated models, however, are not capable of reproducing all of these hydrodynamic processes adequately. Although both models adequately reproduce water levels and the vertical velocity structure, they have difficulties predicting the pycnocline height. Additionally, it is found that both models introduce a small phase shift in the velocity and salinity prediction. The research presented here is a contribution to the understanding of the governing hydrodynamic processes in the Rotterdam Waterway, and the effect that (in)accurate modelling of these processes may have on future studies. Recommendations following from this research could improve future modelling practices.

The port of Rotterdam is located within the Rhine-Meuse estuary where a substantial amount of fine sediment transport takes place. Therefore, the port of Rotterdam is subject to significant siltation, requiring maintenance dredging to guarantee a sufficient nautical depth of fairways and harbour basins. To optimise the dredging strategy in the port of Rotterdam, a pilot study has been carried out wherein sediment is reallocated in the Rotterdam Waterway, during ebb, instead of offshore in the North Sea. Between May and November 2019, 210,000 tons of sediment has been reallocated. This pilot study has been carried out in the context of the larger EU-Interreg Sediment Uses as Resources in Circular and Territorial Economies (SURICATES) project. The main goals are to re-use the sediment as a resource and to reduce the sailing time of the dredging vessels. Both ideas comply with the Building with Nature philosophy; a concept gaining popularity over recent years in The Netherlands focusing on, amongst others the optimisation of dredging strategies. It is expected that the reallocated sediment is mainly transported offshore, while at the same time some of the sediment will nourish the river banks of the Rotterdam Waterway enhancing its flood resilience. This thesis focuses on understanding the fine sediment behaviour of the SURICATES pilot project on two different scales. This is done by analysing measurements and model hindcasts. The measurement campaign is set-up by Deltares and the Port of Rotterdam. For the model hindcasts an operational hydrodynamic and sediment model is used. On the small scale this is done by focusing on the behaviour of a single disposal over a tidal cycle. The large scale focuses on the cumulative long term behaviour of all sediment reallocations. For the small scale two measurement surveys are analysed. In both surveys it is found that a sediment reallocation executed by bow coupling is subject to mixing up to halfway the water column. Subsequently, the sediment plume is advected around and below the pycnocline. Further measurements in the mid field are lacking, but it is hypothesised that the majority of the sediment settles during subsequent low water slack. For the other execution method, drawing the bottom doors, which is used to reallocate the majority of the sediment, useful measurements are absent. It is hypothesised that the majority of this reallocated sediment is confined in the salt wedge and therefore mainly transported upstream over time. To assess the long term behaviour of the cumulative behaviour of all sediment reallocations, a different measurement campaign is set-up. In this measurement campaign, bed samples are collected prior to and during the pilot study to determine the change of the bed composition. In this campaign an indication for increased sedimentation related to the pilot study is found for nearly all the sample locations. The short term model study is set-up to enhance the understanding of the short term behaviour of a sediment plume, to derive an accurate source term for the sediment disposals and to carry out a sensitivity analysis. This sensitivity analysis is executed to derive the influence of differences in disposal method, timing of disposal, and uncertainties in the model. It is found that the execution method has the largest influence on the critical sediment fluxes on the short term, followed by the timing of disposal. From the long term model hindcast, in which the entire pilot study is hindcasted, it is found that 27% of the total amount of reallocated sediment flows downstream from the location of disposal and 73% upstream. These estimations are in line with the hypothesis and long term measurement results, but a thorough calibration of the results is lacking. To conclude, in this thesis a pilot study utilising a different sediment reallocation strategy in the port of Rotterdam has been investigated. This study shows that majority of the sediment disposed, in the current set-up of the pilot study, is estimated to flow upstream. In the sensitivity analysis, it is predicted that this might be caused by the timing of disposal or method of execution. It is also found that the initial behaviour of the sediment plume and the long term measurement contain a large amount of uncertainties. As most important recommendation for future work an expansion of the current measurement survey is proposed with at least two fixed locations: one downstream and one upstream of the location of disposal. In this way sediment fluxes can be established, which can also be used to verify and calibrate the sediment model.

TheAmazon-Orinoco river plume is a buoyant freshwater lens of 1.2 × 10⁶km², which has been traced over 2000 km from the Amazon river mouthinto the Caribbean Sea and along the Lesser Antilles. The river plume iswarmer than the surrounding open-ocean waters, with temperaturedifferences up to 1.5 ∘C caused by a stratification-induced barrier layerinhibiting vertical mixing and coloured matter increasing solar energyabsorption. Due to its magnitude, the river plume affects thehydrodynamics and the oceanic conditions in the Western Tropical NorthAtlantic (WTNA) substantially, but its variations on interannual time scales andthe corresponding relation to local sea-surface temperature (SST) are notwell understood. The Caribbean Sea is a region of high ecological value as itis home to extensive coral reefs, which are especially sensitive topersistent high SST. Therefore, this study investigates the interannualvariability of the Amazon-Orinoco river plume and its relationship to SSTsin the Caribbean Sea and the WTNA. It is hypothesised that fresh anomaliesof the river plume salinity pattern are indicative of a more extensivetransportation of the heat contained in the river plume. As a result, itis expected that interannual variations of dominant river plume pathwaysaffect the magnitude and location of anomalous SSTs. To test thishypothesis, model reanalysis fields of oceanic conditions from 1993 to 2017are used to conduct statistical analyses. In this context, the river plumevariability is determined using specific regions of freshwater influenceestablished using Empirical Orthogonal Function (EOF) analysis ofanomalous sea-surface salinity (SSS). Cross-correlations analysis relatingthese EOF modes of with atmospheric processes show that the interannualvariability of the river plume is dominated by wind-inducedadvective transport and -mixing. Strong winds along the Brazilian shelfare related locally increased SSS, while a weak southward component makesfor extensive spreading of the low-salinity plume waters. Additionally, weshow that high river discharge affects SSS east of the Lesser Antillesafter a lag of three months. Through its modulation of these atmosphericprocesses, there is a strong indication that the El Niño-SouthernOscillation affects SSS variability in the main along-shelf northwestplume pathway, with low SSS 1–9 months after a La Niña event. DecreasedSSS are found in phases 2 and 3 of the Madden-Julian Oscillation, whileincreased SSS was observed in phases 6 and 7. However, the evidence forthis relation is weak and should be investigated in further research.
The results show that, opposed to the hypothesis,a more extensive river plume is not associated with higher SSTs in theCaribbean Sea. However, strong correlations are found between river plumesurface area and SSTs at a lag of 1 year. Based upon results of previousstudies, we argue that the river plume has the ability to pre-heat themixed layer in the WTNA leading to extreme temperatures in the followingyear. It is wise to conduct a Lagrangian parcel back-tracking experimentto verify this mechanism.

The Rhine Region of Fresh Water Influence (ROFI) is a shallow frictional river plume in front of the Dutch coast. Each tidal cycle a new tidal plume front with fresh-water is released. Recently, internal gravity waves have been observed in this plume. Using a Froude number analysis, support for the internal wave generation mechanism by a tidal plume front is found. It is shown, that on averaged neap tides the density stratification is large. This results in a larger area where the internal waves can be released from the tidal plume fronts compared to spring tides. As the Rhine ROFI is located in shallow water, it is investigated what the effect of bathymetry variation is on internal waves. This is studied by extending the standard KdV model derivation to account for a variable bed. This resulted in a small correction on the propagation speed inside the KdV equation. Observations of internal waves have been used to validate the KdV model. By scaling analysis of the observed waves it is obtained that the relative wave height and relative depth balance for most of the observed events. For these events the wave period and velocity amplitude of the KdV model are well matched with the measurements. This showed that the KdV model may be used for a first estimate of internal waves in the Rhine ROFI. By developing a TGE fitting procedure, it was possible to obtain the pycnocline depth and the direction of wave propagation with only limited data available. Therefore, the parameters have been varied and the error between the TGE solution and the velocity potential obtained from the velocity measurements has been minimized.

This study aims to get an insight into the particle trajectories that microplastics follow, after having been released in the North Sea.

For the computations daily-mean values of the surface currents are used, retrieved from the Mercator global ocean model. 2D particles trajectories are simulated for a year, with a 3rd party Python toolbox for Lagrangian simulation of particles: OceanParcels. Particles released from any location in the North Sea eventually get trapped in the Norwegian Coastal Current (NCC). From here they are being further advected to the North, at different moments in time for the particles released at different locations. The coastal processes in the NCC are mainly linked to wind and stratification, hence variations in ow patterns near the coast are linked to the seasons. When these ow pattern include large scale eddies, the particles follow a meandering and erratic path. Floating plastic particles released in the North Sea will flow northwards along the coast of Norway. Eventually those particles will end up in the Arctic region or get
trapped in the Norwegian fjords, independently of the location of release. However, the time scale of the northward advection depends both on where the particle has been released and the environmental conditions.

Coastal bays in the Caribbean accommodate different marine ecosystems, including seagrass meadows and coral reefs, which provide important ecosystem services. However, these marine ecosystems are endangered. Seagrass ecosystems have a key role in coastal bays, but despite their alarming rates of loss, they receive little attention compared to other marine ecosystems. This study aims to provide a better understanding of the role of seagrass ecosystems in coastal bays by assessing their impact on the morphodynamic behavior. In order to achieve this, an existing hydrodynamic model has been extended to include morphodynamics. Two different bays have been investigated: Baie Orientale and Baie de l’Embouchure (St. Martin). These bays represent a partially exposed and fully sheltered bay, respectively. Both regular wave conditions, as well as extreme storm conditions, have been investigated to understand the consequences of seagrass loss on coastal erosion for different environmental climates. The coastal erosion has been found to be more prominent in the exposed region, whereas the more sheltered areas are more resilient. The erosion takes place in shallow waters but remains limited under regular swell conditions. The increased wave energy, during the extreme storm event, increases the erosion rates considerably. However, the regular swell conditions are normative in shaping the morphological development of these coastal bays when longer timescales are considered. The seagrass counters erosion most effectively in the foreshore and much less in deeper regions. Removal of seagrass in the foreshore (between 1-3m) halves the sediment stabilization. In contrast, the sediment stabilization increases if the meadows are able to spread especially towards the shore. In addition, the wave energy and waveform have a significant impact on the sediment stabilization. The stabilization typically decreases for higher and shorter (storm) waves, whereas it increases for smaller and longer (swell) waves. The long-distance interactions between seagrass meadows and coral reefs have led to their mutual coexistence, which is also found to be beneficial for the erosion control of the coastal bays. The dissipation of wave energy on top of the reefs fosters the sediment stabilization by seagrass. Moreover, the impact of the reefs on sediment stabilization also increases in the presence of seagrass meadows due to the additional drag exerted by the seagrass. The seagrass meadows are typically more effective under swell waves, whereas the reefs are more dominant in the dissipation of storm waves. Seagrass meadows and coral reefs form thus a synergy, in which the stabilization of sediments provided by the individual ecosystems is facilitated, reinforced, and complemented by the proximate presence of the other ecosystem. Finally, the role of seagrass on sediment stabilization has been explored under climate change. The coastal erosion increases for the considered sea-level rise scenarios. The impact of the seagrass on sediment stabilization decreases when the bays, seagrass and reef ecosystems are not able to keep up with sea-level rise. The synergy between seagrass and reefs emphasizes that the responses of both ecosystems are of importance for the future of tropical sheltered bays and should therefore not be managed in isolation.

This thesis introduces a new algorithm for optimising shipping routes within a dredging project. Highly dynamic and time-dependent hydrodynamic features influence shipping routes. Due to the complex interactions between the horizontal tide, vertical tide, stratifying forces, wind-driven forces, and limited water-depth, shipping routes were previously only optimised for large scale routes (order of 1000 km). This study presents an algorithm that can optimise shipping routes that are influenced by these small scales (order of 10 km) hydrodynamic features. This algorithm uses graph theory to solve for the time-dependent fastest path between start and destination. Graph theory searches for the optimal path through a set of nodes that are connected with edges. This study uses the time-dependent shortest path algorithm which accounts for the FIFO-criteria (Waiting criteria) and can solve the non-convex nature of the problem

The input of this algorithm is a hydrodynamic model. These models are Computational Fluid Dynamic (CFD) models that calculate currents and water levels in a specific domain. The domain is discretised into cells and nodes to calculate these hydrodynamic features. This study uses the nodes of this hydrodynamic model as the vertices of the graph. However, for some cases, the hydrodynamic model has too many nodes for the shortest path algorithm. This study presents a method for reducing the number of nodes without reducing the spatial resolution. The nodes are reduced based on a combination of the vorticity and the magnitude of the flow.

This algorithm is implemented in a python software package named Hydrodynamic Algorithm for Logistic enhancement Module (HALEM). HALEM can determine the optimal shipping route for a given hydrodynamic model. Defining different cost functions results in different optimisation purposes. This thesis presents cost functions for the fastest route, shortest route, cheapest route and least polluting route. This software is then implemented in the OpenCLSim software so that this combination of software can optimise routes of entire projects. A case study simulates a beach-nourishment at Schouwen Westkop Noord to demonstrate the practical use of HALEM and OpenCLSim. For this project, 425,500 m3 sand should be dredged offshore and pumped onto the beach. Due to the narrow gullies and tidal changes in hydrodynamic features, the routes were hard to predict. The simulation with HALEM and OpenCLSim shows an increase in the production with 21 % compared to the simulation with just OpenCLSim.

Climate models predict increased Arctic precipitation and subsequent Arctic freshening as a response to increased green house gas concentrations. Eulerian studies have shown that with increased Arctic precipitation AMOC (Atlantic Meridional Overturning Circulation) strength decreases. Decrease in AMOC strength comes with a decreased redistribution of heat from lower to higher latitudes which can have severe effects on our climate. Therefore, understanding the effects and mechanisms of Arctic precipitation change is a crucial building block for predicting and possibly preventing climate change. This study used a Lagrangian approach. The pathways of water at Fram Strait were investigated for present-day climate (control run) and two scenario runs with increased Arctic precipitation (+50% and +300% respectively). Importantly, it was found that Arctic water reaches the Labrador Sea through Denmark Strait for all three runs. Thus, the extra fresh water in the Arctic can possibly impact sinking and convection zones in the Labrador Sea. The total amount of Arctic water, passing Denmark Strait from Fram Strait, increases for the weak scenario and decreases for the strong scenario of this study. On the other side of Iceland, for the strong scenario of this study, Arctic water stops passing the Iceland-Faroe-Ridge through the Faroe Bank Channel. The amount of Arctic water going into and staying in the Nordic Seas remained almost unchanged with increased Arctic precipitation. The two routes passing from Fram Strait into the North Atlantic were analysed further with respect to depth changes and properties. On both routes particles were fresher compared to the control run when increasing Arctic precipitation. For the weak scenario particles were usually colder than the control run on both routes. For the strong scenario, particles were only colder at Fram Strait, but got warmer than the control run along the pathway.

Globally, at each river mouth freshwater is discharged to coastal waters. At each river mouth a freshwater plume is formed. These fresh water plumes are also known as the regions of freshwater influence (or ROFI’s). Furthermore, the freshwater plume will affect local coastal water levels. In order to improve global flood forecasting, research on this topic has been conducted. Through a combination of numerical (DFlow-fm) and probabilistic (vine copula) models predictions can be made for the size- and hydraulic effect of these river plumes.

This Master thesis presents an improved model for the transport of fine suspended sediment (SPM) in the Rhine region of fresh water influence (ROFI) and in the Port of Rotterdam. It is known that the transport of SPM in the Rotterdam harbor depends on marine as well as fluvial processes. SPM transport inside the harbor has been modelled in the past, but these studies did not take into account the dynamic nature of SPM concentrations at sea because they employed constant sediment boundary conditions and did not resolve wind waves. At the same time, numerous model studies of SPM transport in the Southern North Sea have been carried out, but none included the Rotterdam harbor as more than a discharge point. This MSc thesis is the first study to combine a North Sea SPM model with a detailed model of the Port of Rotterdam. By doing so, it incorporates the dynamic nature of sediment concentrations in the North Sea, and its effect on sediment transport in the Rotterdam harbor.

A new suspended sediment model of the Rotterdam harbor was set up with the
Delft3D-WAQ software package and the NSC-Course model grid already in use by Port of Rotterdam. A validated hydrodynamic model was available for this grid schematization, and has been used as model forcing. At the open sea boundaries, the validated ZUNO-DD model for SPM transport on the North Sea was used to compute suspended sediment concentrations, which were then applied as boundary conditions. In addition, a wave buoy assimilation technique was used to include the effect of wind
waves on the resuspension of sediment. A set of initial conditions was created by repeating a single spring-neap cycle with calm, virtually waveless conditions until sediment concentrations reached a dynamic equilibrium. A byproduct of this method was that it created a large sediment availability in the bed, which lead to higher than expected resuspension during storms. To investigate the response of the system to variable forcing conditions, three 14 day periods have been selected from measured environmental conditions of the year 2007, containing different combinations of storms and river discharges.

The model shows positive results, but does need improvement. Principally, the episodic nature of siltation in the harbor basins lining the mouth of the Rotterdam Waterway was reproduced. In line with observations, storms at sea correspond to high SPM concentrations at sea, large influxes through the harbor mouth, and increased siltation rates in the harbor basins. Furthermore, differential advection of the salinity structure and the trapping of SPM in the Rotterdam Waterway was observed. However, erosion was also observed, and harbor basins lining the Rotterdam Waterway showed a strong response to events at sea which contradicts observations. This is attributed to
the wave stress assimilation performed on NSC-Coarse grid which generates large wave stresses inside the harbor and causes unwanted resuspension. Three different forcing scenarios were chosen to highlight the differences between marine and fluvial processes, but a clear distinction could not be made. This is expected to be caused by the erroneous wave shear stresses which obfuscate the effects of different siltation mechanisms. It is recommended that this aspect of the model is improved and that the model is reassessed using measurements.

In this report a first assessment is made of the impact of hurricanes on Saint Martin. Two bays at the east coast of Saint Martin, Baie Orientale and Baie de L'Embouchure, are studied under the conditions of hurricane Irma and Maria, which both have passed over in September 2017. During these hurricanes the bays were exposed to large wind speeds up to 40 [m/s]. Also the bays were subject to runoff due to large rainfall. The runoff has been flowing into the bays as fresh water, but also as high saline waters. This due to the salt ponds which are located in the region.

The model has been run in Delft3D and as result large wind-setup driven circulations were obtained. These circulations caused a drop of waterlevel during Irma and an increase under Maria. The saline and fresh waters discharged during the hurricanes were mostly kept onshore and is distributed over the two basins.

Shallow bays in the Caribbean, like Baie Orientale and Baie de L'Embouchure in Saint Martin, are often sheltered by coral reefs and covered by seagrass meadows. They provide valuable services as tourism and coastal protection. The ecosystems are linked through biological, chemical and physical processes. But they are under pressure due to sea-level rise. The response of one of the ecosystems to climate change could impact the other ecosystem.

In order to predict the impact of sea-level rise on the biogeomorphology in Baie Orientale and Baie de L'Embouchure, the hydrodynamic model Delft3D Flexible Mesh is applied. The effect of seagrass meadows and coral reefs on both flow and waves are captured with this model. In this way, the long term change in average hydrodynamic conditions due to sea-level rise is determined depending on the response of the ecosystems.

A wave-driven circulation is found in both bays with flows of 0.5 m/s over the reefs and currents of 0.2 m/s inside the bays. The hydrodynamic conditions are mainly determined by the reef height. Depending on the response of coral reefs to climate change and the amount of sea-level rise, the wave height inside the bays and the wave-induced currents increase. Under the worst-case scenario, where coral reefs degrade and seagrass meadows die, flow velocities increase by more than 100% in Baie de L'Embouchure and by 200% in Baie Orientale under a sea-level rise of 0.87 m. The significant wave height rises to 300% in Baie Orientale and doubles in Baie de L'Embouchure. But this increase of hydrodynamic stresses is not expected to lead to devastating damage to coral reefs and seagrass meadows. Instead, the response of coral reefs will be determined by changing water temperatures and ocean acidification. A shift in seagrass occurrence due to the changed hydrodynamics is expected.

The long term impact of sea-level rise on the biogeomorphology of Baie de L'Embouchure and Baie Orientale seems to be limited. The ability to mitigate the impact of sea-level rise is shown and the resilience of the ecosystems proved, which is very promising for other shallow Caribbean bays that are threatened by sea-level rise.