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 routing weighting¶
Sediment routing probability for a given cell to neighbor cell is computed according to:
where is the local routing direction and is a unit vector pointing to neighbor from cell , and is the cellular distance to neighbor ( for cells in main compass directions and for corner cells. is a resistance estimated as an inverse function of local water depth ():
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:
where is the volume of sediment to be eroded or deposited from the bed at a given cell along the parcel walk.
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;
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.
Model stability depends on…
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.
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
In the following example,
N_crossdiff takes the default value.
The impact of topographic diffusion is minor compared to the bed elevation change driven by parcel erosion or deposition (sand and mud routing effects).
The reference volume () impacts model stability. This volume characterizes the volume on one inlet-channel cell, from the channel bed to the water surface:
where is the inlet channel depth (meters) and is the cell length (meters).