Modeling water flow and volumetric water content in a degraded peat comparing unimodal with bimodal porosity and flux with pressure head boundary condition

Degraded peatlands release large amounts of greenhouse gases. The development of effective mitigation and management measures requires an understanding of relevant site‐specific biogeochemical and hydraulic processes. However, the simulation of water fluxes and vadose zone state variables of degrading peatlands relies on proper process description, parameterization of hydraulic functions, and representation of the boundary conditions. The objective of this study was to analyze the effects of unimodal versus bimodal soil hydraulic functions and pressure head versus flux‐type lower boundary conditions (LBCs) on the calculated hydraulic characteristics of a degraded peat profile. HYDRUS‐1D was used to study the hydraulic flow dynamics parameterized with data from a weighable groundwater lysimeter for the period from May 1 to December 31, 2019. Simulations comparing uni‐ and bimodal hydraulic functions showed only minor differences. Simulations of soil water pressure at a depth of 30 cm using a flux‐type LBC (RMSE: 27 cm, where RMSE is root mean square error) performed better than simulations using a pressure head LBC (RMSE: 48 cm). The pressure head LBC performed better at simulating volumetric water contents in 30‐cm depth than the flux LBC variant (RMSE: 0.05 vs. 0.09 cm3 cm−3). For specific site conditions with a shallow, fluctuating groundwater table and temporary air entrapment, the choice of LBC was important for a more accurate simulation of soil water fluxes and volumetric water content.