| ... | @@ -9,13 +9,13 @@ The DEMO model has a very simple formulation for the hydrodynamics (exchange bet |
... | @@ -9,13 +9,13 @@ The DEMO model has a very simple formulation for the hydrodynamics (exchange bet |
|
|
flowchart LR;
|
|
flowchart LR;
|
|
|
Comp_i <==>|Ex_i,j| Comp_j;
|
|
Comp_i <==>|Ex_i,j| Comp_j;
|
|
|
```
|
|
```
|
|
|
The same formulation is used for the exchange between the North Sea and the compartment neighboring the North Sea. The exchange parameters ($`Ex_{i,j}`$) are estimated by fitting a tracer model to data, artificially generated by a tracer run with the 3-D hydrodynamic model of Jiang et al. (2019)[^3049].
|
|
The same formulation is used for the exchange between the North Sea and the compartment neighboring the North Sea ($`Ex_{i,NS}`$). The exchange parameters ($`Ex_{i,j}`$ and $`Ex_{i,NS}`$) are estimated by fitting a tracer model to data, artificially generated by a tracer simulation with the 3-D hydrodynamic model of Jiang et al. (2019)[^3049].
|
|
|
|
|
|
|
|
[Top](Hydrodynamics)
|
|
[Top](Hydrodynamics)
|
|
|
# 3D Hydrodynamic model
|
|
# 3D Hydrodynamic model
|
|
|
A 3D hydrodynamic model has been developed by Jiang et al. (2019) [^3049]. The model was developed in the open source General Estuarine Transport Model [GETM](https://getm.eu) in combination with the General Open Turbulence Model [GOTM](http://gotm.net). The model has a 300x300 meter rectangular grid in the Oosterschelde and a part of the Voordelta in the North Sea.
|
|
A 3D hydrodynamic model has been developed by Jiang et al. (2019) [^3049]. The model was developed in the open source General Estuarine Transport Model [GETM](https://getm.eu) in combination with the General Open Turbulence Model [GOTM](http://gotm.net). The model has a 300x300 meter rectangular grid in the Oosterschelde and a part of the Voordelta in the North Sea.
|
|
|
|
|
|
|
|
With the model a tracer simulation was run. At the start of the simulation the fraction of North Sea water in the North Sea was 1 (100%) and in he Oosterschelde the fraction of North Sea water was 0 (0%). The model was run for one year from 1 January 2009 to 31 December 2009. For each hour, the depth-averaged tracer concentration and the water level in each (300x300 m) grid-cell was exported. The intrusion of North Sea water into the Oosterschelde with the tidal currents is shown in the movie below.
|
|
With the 3D model a tracer simulation was done. At the start of the simulation the fraction of North Sea water in the North Sea was 1 (100%) and in he Oosterschelde the fraction of North Sea water was 0 (0%). The model was run for one year from 1 January 2009 to 31 December 2009. For each hour, the depth-averaged tracer concentration and the water level in each (300x300 m) grid-cell was exported. The intrusion of North Sea water into the Oosterschelde with the tidal currents that was calculated by the model is shown in the movie below .
|
|
|
|
|
|
|
|

|
|

|
|
|
|
|
|
| ... | @@ -24,7 +24,7 @@ The results of the tracer simulation with the 3D hydrodynamic model were used to |
... | @@ -24,7 +24,7 @@ The results of the tracer simulation with the 3D hydrodynamic model were used to |
|
|
|
|
|
|
|
[Top](Hydrodynamics)
|
|
[Top](Hydrodynamics)
|
|
|
# Calibration exchange rates
|
|
# Calibration exchange rates
|
|
|
The exchange rates between the neighboring compartments are estimated by fitting the output of the DEMO tracer model to the artificial data generated by the GetM/GotM model. Initial conditions for each compartment were calculated from the artificial data as the average fraction of North Sea water in each compartment during the first four tidal cycles. The models are fitted using the modFit function in R (Soetaert and Petzold, 2010)[^3228]. The 1D model and the models with the 4 and 28 compartments are fitted separately. The parameters were estimated using the Levenberg-Marquardt algorithm where the objective was to minimize the sum of squared residuals (SSR) between model and "observations".
|
|
The exchange rates between the neighboring compartments are estimated by fitting the output of the DEMO tracer model to the artificial data generated by the 3D GetM/GotM model. Initial conditions for each compartment were calculated from the artificial data as the average fraction of North Sea water in each compartment during the first four tidal cycles. The models are fitted using the modFit function in R (Soetaert and Petzold, 2010)[^3228]. The 0D model and the models with the 4 and 28 compartments are fitted separately. The parameters were estimated using the Levenberg-Marquardt algorithm where the objective was to minimize the sum of squared residuals (SSR) between model and "observations".
|
|
|
|
|
|
|
|
[Top](Hydrodynamics)
|
|
[Top](Hydrodynamics)
|
|
|
## 0D model
|
|
## 0D model
|
| ... | | ... | |