Field and numerical experiments of subsurface drainage systems in saline and clay interlayered fields in arid regions

A reasonable layout of subsurface drainage systems is considered essential for maximizing its drainage and salt control effectiveness. In the saline-alkali farmland of arid regions in Northwest China, clay interlayers are common; however, the influence of clay interlayers on the layout of the subsurface drainage has not been extensively considered in the literature. This study investigated the process of subsurface drainage and salt discharge in salt-affected fields with clay layers using field experiments and numerical simulations. Four field experiments were conducted, considering three different relative positions between the drainage pipes and clay interlayers. The results showed that the clay interlayers hindered water infiltration; however, the distribution of soil salinity in the soil profile varied with the buried depth of drainage pipes at different positions relative to the clay layer. When the buried depth of drainage pipes increased, the amount of water and salt discharged from drainage pipes increased, and the increase rate in salt discharge was greater than water drainage. A numerical model was calibrated and validated using the field experiment data, and 25 orthogonal numerical experiments were conducted to investigate the soil desalination effects of buried depth of drainage pipes, spacing between the pipes, permeability of the interlayer, and position of the clay interlayer. The results showed that the drainage pipe buried depth, spacing, and permeability of the clay layer had significant effects on the desalination rate (P 0.05). Therefore, subsurface drainage pipes should be placed below the clay interlayer. The desalination rate linearly increased with the buried depth of drainage pipe and permeability of the interlayer, and it increased exponentially with decreased spacing. An empirical formula for soil desalination rate considering interlayer and subsurface drainage pipe layout parameters was fabricated and used to determine the appropriate layout parameters.