Compaction formula (Subsidence Overlay)

The amount of subsidence due to compaction is calculated by the following formula:

$S_{sub}=(a_{peat}\cdot p_{x,y}+a_{tl}\cdot tl_{x,y})\cdot log10(t_{days})+a_{h}\cdot \Delta {h}+b_{p}\cdot p_{x,y}+b_{tl}\cdot tl_{x,y}$

where

• $p_{x,y}$ is the peat fraction of the cell at coordinate (x,y)
• $tl_{x,y}$ is the top layer thickness of the cell at coordinate (x,y)
• $\Delta _{h}$ The height increase as a result of the actions taken during a session, such as the creation of levees.
• $t_{days}$ is the number of days in a year times the current year being calculated.
• $a_{peat}$ is the first regression constant (0.015853041) for the effect of the peat fraction and time in days.
• $a_{tl}$ is the first regression constant (0.006617643) for the top layer thickness and time in days.
• $a_{h}$ is a regression constant (0.200468677) for the change in surface height, for example caused by added materials.
• $b_{peat}$ is the second regression constant (0.02348519) for the effect of the peat fraction.
• $b_{tl}$ is the second regression constant (-0.010061616) for the top layer thickness.

This formula is based on provided expert data in the form of a reference table, indicating the amount of subsidence based on the parameters used in the formula above. The formula's results conform to the reference table to within an average of a tenth of the margin of error of the original table.