Infiltration into Frozen Silty Clay Loam Soil with Different Soil Water Contents in the Red River of the North Basin in the USA

Predicting surface runoff and flooding in seasonally frozen areas such as the Red River of the North Basin (RRB) in USA is a challenging task. It depends on the knowledge of the complex process of infiltration in frozen soil, such as phase changes of water, ice content and distribution in the infiltration zone (the top 0−30 cm of the soil profile), soil pore size distribution, soil temperature and freeze−thaw cycles. In this study, the infiltration rates into frozen soil (Colvin silty clay loam according to the United States Department of Agriculture (USDA) Classification, and Chernozem according to Food and Agriculture Organization of the United Nations (FAO) international soil Classification) were measured at three different initial water contents: permanent wilting point (PWP), θ_pwp; field capacity (FC), θ_fc; and between FC and PWP, θ_mid. Laboratory infiltration experiments were conducted using a Cornell sprinkle infiltrometer with three replications for each initial water content. Volumetric soil water content (θ_v) and soil temperature at three depths were also continuously monitored using sensors. The average infiltration rates were 0.66, 0.38, and 0.59 cm/min for three initial water contents (θ_pwp, θ_mid, and θ_fc, respectively). Initial infiltration into frozen soil occurred quickly in the soil with θ_pwp because the soil was dry. Melted ice water contributed to the total soil water content over time, so it made the initial infiltration comparatively slower in the soil with θ_mid. Initial infiltration was also slower in the soil with θ_fc because the wet soil had very small pore space, so the soil rapidly reached its saturation after the infiltration started. The Horton infiltration equation was fitted with the observed infiltration rates for the soils with three initial water contents, and the goodness of fit was evaluated by using the coefficient of determination (R²) and the root-mean-square error (RMSE). The final infiltration rates from the fitted Horton equations were 0.060, 0.010, and 0.027 cm/min for the initial water contents (θ_pwp, θ_mid, and θ_fc, respectively). The soil water content along the soil profile changed with the amount of infiltrating water over time. However, the initial soil water content and melt water from ice resulting from soil temperature rise regulated the change in soil water content. The amount of ice melt water contribution to soil water content change varied among the soils with different initial water contents (θ_pwp, θ_mid, and θ_fc, respectively). The θ_v changed gradually in the θ_pwp soil, rapidly at 0 °C in the θ_mid soil, and less in the θ_fc soil. The change in pore distribution due to freeze−thaw cycles and soil packing altered the soil hydraulic properties and the infiltration into the soil. This study can provide critical information for flood forecasting model and subsurface drainage design in the RRB.