A modified Green-Ampt model was developed to simulate water infiltration in fine-textured soil with a coarse interlayer. Because under such a soil structure, the two soils may not be fully saturated during infiltration, the model introduced two parameters-that is, the saturation coefficients a and b, to reflect the incomplete saturation condition and their influence on infiltration processes. In order to analyze the variation pattern of the two parameters in the above proposed model, scenarios were set for soil column infiltration in fine-textured soil with a coarse interlayer under different buried depths. A Richards equation-based model (RE-Model) was built for simulating the above scenarios and to obtain the evolution of soil water content along the soil profiles. Simulation results show that the infiltration rate decreased to a constant value when the wetting front crossed the upper interface between the fine and coarse soil layer. The soil matrix suction (ψ 2 ) at the upper interface remained unchanged after the wetting front advanced into the coarse layer, and the steady value of ψ 2 showed a linear relationship with the buried depth of the coarse layer. Based on the simulation results of the RE-Model, a method was proposed to determine the saturation coefficients related to the relative hydraulic conductivity and water content at ψ 2 in the modified Green-Ampt model. Then, the modified model was tested under various infiltration conditions with different soil layered structures, and the results showed good agreement with the experimental data.Water 2019, 11, 787 2 of 17 Wang et al. [15] found that the constant infiltration rate decreased with the buried depth of sand interlayers in loess soils after the wetting front reached the upper interface of the sand layer.Infiltration models have been developed to understand the laws embedded in experiments and evaluate the soil water dynamics in different infiltration cases [16,17]. Empirical models, such as the Kostiakov equation (1932), Horton equation (1940), and Philip equation (1957) are widely used because of their simplicity and capability for infiltration computation. However, these models are usually used in homogeneous soils, and may be questionable for simulating water infiltration in layered soils [9]. The Green and Ampt equation (1911) is a half-empirical and half-theoretical model describing the soil water infiltration process, originally in homogeneous soil and later extended to layered soils [18,19]. Wang et al. [13] developed a modified Green-Ampt model for two-layer soil infiltration, and found that the top soil layer's thickness greatly affects layered soil water infiltration. To estimate unstable infiltration in layered soil with non-uniform initial soil water content, Liu et al. [19] derived an infiltration model based on the Green-Ampt method. These modified Green-Ampt models were claimed to reasonably describe infiltration in layered soils. However, the water content and hydraulic conductivity of the wetted zone were usually assumed to be th...