Find out and document what `waveHybridCA` does

[NiklasHohmann]

Currently, usage of both the simpleWave and waveHybridCA result in dominance of the light green facies, and 100 % stratigraphic completeness. This effect (1) is present in both attached and isolated carbonate platforms and (2) vanishes as soon as sea level fluctuations are present. It might not be a boundary effect (bc it is present in both isolated and attached platforms) and not depend on the CA (bc it is present in simpleWave, which has no CA)

[NiklasHohmann]

The effect is absent under productionCA (produces a mound geometry) and waveCA (drowns) when starting from the simple isolated platform bathymetry (isolated_platform.txt)

[Burgesski]

Can you pop an image of the output in here for me to have a look at? Also, can you check and tell me the grid size for this model e.g. 200 x10 or perhaps 200x20 would be a typical size. Also, perhaps pop me the paramter input files in an e-mail and I can run it on my version of CarboCAT lite and see what happens...

[NiklasHohmann]

Image sent via email. For the record: @xyl96 and I found that under constant SL, wave energy across the whole platform top is zero. In this case, the simpleWave procedure then always selects the facies that can best deal with low wave energy, which will then immediately fill the accommodation space generated by subsidence.

[NiklasHohmann]

We changed the wave energy tolerance parameters to be 0 (minimum) and 1 (maximum) for all factories, making them effectively independent of wave energy (which is effectively removing any environmental niche). This lead to the following results:

• for simpleWave, the factory dominating the platform changes to factory 1. This is because if all factory have the same wave energy tolerance, the algorithm always picks the first one, making the results dependent on the enumeration of factory
• for waveHybridCA, the platform drowns. Before drowning, the factory with highest production rate dominates, and there is spatial heterogeneity visible.

[xyl96]

wave hybrid: first calculate wavesimple: if multiple facies fit in, just choose the one has the smallest WaveEnergyImbalance. If no facies has been found in this process, i.e., wave energy selection does not work, run CA and production. This is why wavehybridCA looks so similar to simplewave: the wave (niche) is calculated first.

[Burgesski]

If you check Burgess (2013) it explains how the CA rotates the facies checking order to avoid the first choice bias. This snippet of code from the basic CA function could be popped into waveHybridCA too, of course, to cover for the sitations where the wave energy and water depth ranges do not determine which facies to choose. Or perhaps just turn wave energy calculation off and use the basic CA? Depends on the philosophy behind the triangle experiment - what EXACTLY are we trying to test with the triangle? I know in general, but perhaps some more specific thought is required to decide on the final CarboCAT formulation to use; one of my PhD examiners said to me "Model less, think more" and they definitely had a point ;-)