Solutions of the linearized Richardson–Richards Equation (RRE) for one‐dimensional soil water flow in layered soils with sinusoidal flux in the frequency domain are derived. We evaluate the accuracy of our analytical and other analytical solutions by comparing them with results from a standard numerical model. Our analytical solution agrees with the numerical solution under multi‐layered heterogeneous soil, while others disagree. We also demonstrate the capability of the proposed solution to simulate soil moisture dynamics under a realistic, multi‐frequency flux case. The procedure described in the paper is valid for any series of arbitrary periodic flux superpositions for layered heterogeneous . Moreover, our solution is efficient in the calculation compared with numerical solutions, especially when dealing with long‐time series soil moisture, which can provide a validation of numerical models.