> For the complete documentation index, see [llms.txt](https://swatplus.gitbook.io/io-docs/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://swatplus.gitbook.io/io-docs/theoretical-documentation/section-6-management-practices/septic-systems/biozone-algorithm/soil-moisture-and-percolation.md).

# Soil Moisture and Percolation

Soil porosity (or saturated moisture content) is generally constant in natural soil; however, the porosity of biozone changes with time. The actual porosity of biozone decreases as the suspended solids from STE accumulate in the pore space and the mineralized biomass (dead body) increases in the biozone.

                                                    $$\theta\_s=\theta\_{si}-\frac{plaque}{\rho\_{bm}}$$                                                                   (8)

where $$\theta\_{si}$$ is initial soil porosity with zero plaque (mm). The moisture content at each time step is estimated using the mass balance of water within the biozone.

                                                 $$\theta^t=\theta^{t-1}+\frac{Q\_{STE}}{10\*A\_d}-I\_p-ET-Q\_{lat}$$                                    (9)

where $$ET$$ is evaportranspiration from biozone (mm/day) and $$Q\_{lat}$$ is lateral flow (mm/day). Percolation to a subsoil layer is triggered if moisture content exceeds the field capacity in the biozone layer. Potential percolation is the maximum amount of water that can percolate during the time interval.    

                                                  $$I\_{p,pot}=K\_{bz}\*\Delta t$$                                                                   (10)

where $$I\_{p,pot}$$ is the potential amount of percolation (mm/day). The amount of water percolating to the sub-soil layer is calculated using storage routing methodology (Neitsch et al., 2005).

           $$I\_{p,excess}=(\theta -\theta\_f)(1-exp\[\frac{-\Delta t}{TT\_{perc}}])$$         if      $$\theta > \theta\_f$$

            $$I\_{p,excess}=0$$                                                  if      $$\theta \le \theta\_f$$                                            (11)

where  $$I\_{p,excess}$$ is the minimum amount of percolation (mm/day) and $$TT\_{perc}$$ is travel time for percolation $$(TT\_{perc}=(\theta\_s-\theta\_f)/K\_{bz})$$ in hour. The actual percolation is the smaller of the potential and minimum percolation.

                                                    $$I\_p=min(I\_{p,pot},I\_{p,excess})$$                                                    (12)