Model help:CEM

CEM

The CEM model, the Coastline Evolution Model, simulates the evolution of a shoreline due to gradients in breaking-wave-driven alongshore sediment transport. The original CEM has been componentized to consist of the longshore transport module (CEM) and a wave input module (WAVES).

Extended model introduction

The CEM model assumes that the coast consists of a high percentage of mobile sediment and its other assumptions are more applicable at shoreline lengths of km’s and larger. The model was initially designed to investigate an instability in the shape of the coast caused by waves approaching with ‘high’ angles (with the angle between offshore wave crests--i.e. before transformation over approximately shore-parallel contours--and the coast > 45 degrees).
Although a number of wave (and geometry) parameters can be entered, the most vital input control for CEM is the wave climate. The current version of the CEM is driven by simplified directional wave climate controlled by two main input parameters: the asymmetry of the incoming waves angle and the proportion of high-angle waves. This model is not designed to accurately simulate a specific geographic location in detail but rather to more generally represent how a shoreline with highly mobile sediment may respond to varying wave angles. The value in this model is in the breadth it offers in representing how different wave climates can result in different potentially interesting shoreline configurations. Ashton and Murray (2006b) present a more thorough description of the model parameters and theoretical underpinning.

Model parameters

CEM does not need input files from the user, its input is entirely specified in the CMT graphical user interface. To obtain output from this component make sure you toggle on the output files; as a default they are OFF.

Input Files and directories

Run Parameters

Grid

Output Grids

Parameter	Description	Unit
Site prefix	Part of the input and output file name e.g. the name of the location, or project	[-]
Case prefix	Part of the input and output file name that provides you the opportunity to do different scenario simulations for e.g. the same location, or project	[-]

Parameter	Description	Unit
Run duration	Number of simulation time steps	[days]
Shoreface slope	Longitudinal Slope of the Shoreface	[-]
Shoreface depth	Critical threshold depth defining the shoreface	[m]
Shelf slope	Longitudinal Slope of the Shelf	[-]

Parameter	Description	Unit
Number of grid rows	number of rows in the computational grid	[-]
Number of grid columns	number of columns in the computational grid	[-]
Spacing of grid rows	spaces of rows in the computational grid	[m]
Spacing of grid columns	spaces of columns in the computational grid	[-]

Parameter	Description	Unit
Output directory	path to output grid files	[-]
Interval between output files	Specify how often you want output to be written; Notice: files will be really large if you choose to save every timestep.	[-]
Depth file	Evolving land and ocean bottom grid	[m, grid-stacks in NetCDF format]
Percent_filled file	Records the fill of the individual cells for the entire grid.	[%, grid-stacks in NetCDF format]

Coupling parameters

Uses ports

CEM requires wave parameters as supplied by the WAVES component.(Model Help of WAVES).
CEM requires a sediment flux (bedload) originating from one or more river distributary channels. This is provided by the AVULSION component.

Provides ports

CEM provides an elevation grid of offshore geometry after erosion and deposition is done in every timestep.
The elevation is required for the AVULSION component.

Main equations

• Alongshore sediment transport

[math]\displaystyle{ Q_{s} = K_{2} H_{0} ^ \left ( {\frac{12}{5}} \right ) T ^ \left ( {\frac{12}{5}} \right ) cos ^ \left ( {\frac{6}{5}}\right ) \left ( \Phi _{0} - \theta \right ) sin \left ( \Psi _{0} - \theta \right ) }[/math]

(1)

[math]\displaystyle{ K_{2} = \left ( {\frac{\sqrt{g \gamma}}{2 \pi}}\right ) ^ \left ({\frac{1}{5}} \right ) K_{1} }[/math]

(2)

• Predictions of shoreline evolution

[math]\displaystyle{ {\frac{d \eta}{d t}} = - {\frac{K_{1}}{D_{sf}}} H_{b} ^ \left ({\frac{5}{2}} \right ) {\frac{d^2 \eta}{d x^2}} }[/math]

(3)

[math]\displaystyle{ cos^ \left ( {\frac{1}{3}}\right ) \left ( \Phi_{b} - \theta \right ) \approx 1 }[/math]

(4)

• Alongshore component of the radiation stress

[math]\displaystyle{ S_{xy} = H^2 sin \left ( \Phi - \theta \right ) cos \left ( \Phi - \theta \right ) }[/math]

(5)

• Shoreline instability

[math]\displaystyle{ {\frac{d \eta}{d t}} = - {\frac{K_{2}}{D}} H_{0} ^ \left ({\frac{12}{5}} \right ) T^ \left ( {\frac{1}{5}} \right ) \{cos^\left ( {\frac{1}{5}} \right ) \left ( \Phi_{0} - \theta \right ) [ cos^2 \left ( \Phi _{0} - \theta \right ) - \left ( {\frac{6}{5}} \right ) sin^2 \left ( \Phi_{0} - \theta \right )] \} {\frac{d^2 \eta}{d x^2}} }[/math]

(6)

• Shoreline adjustment

[math]\displaystyle{ \Delta F = {\frac{\left ( Q_{in} - Q_{out} \right )}{\left ( D_{sf} + B \right ) \Delta W^2}} }[/math]

(7)

[math]\displaystyle{ \Delta t \propto \left ( {\frac{K_{1}}{D_{sf}}} H_{0} ^ \left ({\frac{12}{5}} \right ) T^\left ( {\frac{1}{5}}\right ) \right ) \Delta x^2 }[/math]

(8)

[math]\displaystyle{ \Delta Y_{bb} = {\frac{\left ( D_{sf} + B \right ) }{\left ( D_{bb} + B \right )}}\Delta Y_{sl} }[/math]

(9)

[math]\displaystyle{ W_{c} = W_{0} + \Delta Y_{bb} - \Delta Y_{sl} }[/math]

(10)

[math]\displaystyle{ \Delta Y_{bb} = {\frac{\left ( W_{c} - W_{0} \right )}{\left ( 1 - {\frac{\left ( D_{bb} + B \right ) }{\left (D_{sf} + B \right )}} \right )}} }[/math]

(11)

[math]\displaystyle{ \Delta Y_{sl} = {\frac{\left ( W_{c} - W_{0}\right )}{\left ({\frac{\left ( D_{sf} + B \right ) }{\left ( D_{bb} + B \right ) } - 1} \right )}} }[/math]

(12)

Notes

The parameters: Shoreface Slope, Shoreface Depth and Shelf Slope set the initial geometry of the shoreface domain and the shelf domain. Simulations will use the shoreface depth as an effective erosion depth, but deposition can take place to deeper depths if the shoreface is accreting on a deeper shelf.

Examples

An example run with input parameters as well as a figure / movie of the output

Follow the next steps to include images / movies of simulations:

• Upload file: https://csdms.colorado.edu/wiki/Special:Upload
• Create link to the file on your page: [[Image:<file name>]].

See also: Help:Images or Help:Movies

Developer(s)

Andrew Ashton

References

• Ashton A., Murray B.A. Arnault O. Formation of Coastline Features by Large-Scale Instabilities Induced by High-Angle Waves. Nature Magazine. Volume 414. 15 November 2001
  
• Ashton A.D., Murray A.B. High-Angle Wave Instability and Emergent Shoreline Shapes: 1. Wave Climate Analysis and Comparisons to Nature. Journal of Geophysical Research. Volume 111. 15 December 2006.
  
• Ashton A.D., Murray A.B. High-Angle Wave Instability and Emergent Shoreline Shapes: 2. Wave Climate Analysis and Comparisons to Nature. Journal of Geophysical Research. Volume 111. 15 December 2006.
  

Links

Movies generated with the stand-alone CEM model are documented in the CSDMS movie gallery. https://csdms.colorado.edu/wiki/Coastal_animations