Morphodynamics

pyDeltaRCM approximates sediment dispersal through the use of a weighted random walk dictated by water flux. In turn, sediment dispersal drives bed elevation change in the model domain by mass conservation.

See 1 for a complete description of morphodynamic assumptions in the DeltaRCM model. In this documentation, we focus on the details of model implementation, rather than model design.

Sediment Transport

Note

Incomplete.

Sediment routing weighting

Note

Incomplete.

Sediment routing probability for a given cell j to neighbor cell i is computed according to:

w_i = \frac{\frac{1}{R_i} \max(0, \mathbf{F}\cdot\mathbf{d_i})}{\Delta i},

where \mathbf{F} is the local routing direction and \mathbf{d_i} is a unit vector pointing to neighbor i from cell j, and \Delta_i is the cellular distance to neighbor i (1 for cells in main compass directions and \sqrt{2} for corner cells. R_i is a resistance estimated as an inverse function of local water depth (h_i):

R_i = \frac{1}{{h_i}^\theta}.

Here, \theta takes the value of coeff_theta_sand for sand routing probabilities, and coeff_theta_mud for mud routing.

(png, hires.png)

../_images/sediment_weights_examples.png

Changes in the bed elevation

Along the walk of a sediment parcel, the sediment parcel volume is modulated on each step, according to the sediment transport rules described above in Sediment Transport. As the volume of the sediment parcel changes, the channel bed elevation at the current parcel location is updated to reflect this volume change (_update_fields), i.e., the bed is eroded or sediment is deposited on the bed. The vertical change in the channel bed is dictated by sediment mass conservation (i.e., Exner equation) and is equal to:

\Delta \eta = \Delta V / dx^2

where \Delta V is the volume of sediment to be eroded or deposited from the bed at a given cell along the parcel walk.

Note

Total sediment mass is preserved, but individual categories of sand and mud are not. I.e., it is assumed that there is an infinite supply of sand and/or mud to erode and entrain at any location in the model domain.

Following a change in the bed elevation, the local flow depth is updated and then local flow velocity is updated according to fluid mass conservation (i.e., uw = qw / h; _update_fields; 1).

Sediment parcels are routed through the model domain step-by-step and in serial, such that changes in the bed elevation caused by one sediment parcel will affect the weighted random walks of all subsequent sediment parcels (Sediment routing weighting), due to the updated flow field.

Sediment parcel routing is handled by first routing all sand parcels, applying a topographic diffusion (see below and topo_diffusion()), and then routing all mud parcels. The impact of routing all sand and mud parcels on bed elevation is shown in the table below.

initial bed

route_all_sand_parcels()

route_all_mud_parcels()

(png, hires.png)

../_images/_initial_bed_state.png

(png, hires.png)

../_images/route_all_sand_parcels.png

(png, hires.png)

../_images/route_all_mud_parcels.png

Model Stability

Model stability depends on…

Note

Incomplete.

Topographic diffusion

Abrupt change in bed elevation (i.e., steep local bed slope) may lead to numerical instability. To prevent this, a topographic diffusion is applied immediately following the routing of all sand parcels in the model sequence.

Hint

Topographic diffusion is applied between routing sand parcels and routing mud parcels.

In implementation, topographic smoothing convolves topography with 3x3 cell kernels configured to a diffusive behavior. The diffusion is repeated over the entire model domain N_crossdiff times. In the following example, N_crossdiff takes the default value.

(png, hires.png)

../_images/topo_diffusion.png

The impact of topographic diffusion is minor compared to the bed elevation change driven by parcel erosion or deposition (sand and mud routing effects).

Reference Volume

The reference volume (V_0) impacts model stability. This volume characterizes the volume on one inlet-channel cell, from the channel bed to the water surface:

V_0 = h_0 {\delta_c}^2

where h_0 is the inlet channel depth (meters) and \delta_c is the cell length (meters).

References

1(1,2)

A reduced-complexity model for river delta formation – Part 1: Modeling deltas with channel dynamics, M. Liang, V. R. Voller, and C. Paola, Earth Surf. Dynam., 3, 67–86, 2015. https://doi.org/10.5194/esurf-3-67-2015