Rock fragment and spatial variation of soil hydraulic parameters are necessary on soil water simulation on the stony-soil hillslope

Whether the spatial variation and influence of rock fragment contents (RFCs) on soil hydraulic parameters (SHPs) should be jointly considered in the soil water simulation have been less investigated in previous studies. In this study, we tested whether considering these factors were necessary in the soil water simulation on a representative stony-soil hillslope located in Taihu Lake Basin, China. Five schemes of SHPs in HYDRUS-3D model were established. They were (i) MultiRFHet: spatially varied SHPs extracted from Durner's multimodal retention function based on the observed soil water retention data; (ii) RosHom: spatially uniformed SHPs derived from ROSETTA; (iii) RosHet: spatially varied SHPs derived from ROSETTA; (iv) RosRFCHom: spatially uniformed SHPs derived from ROSETTA and adjusted by RFCs; (v) RosRFCHet: spatially varied SHPs derived from ROSETTA and adjusted by RFCs. Results indicated when the spatial variation and influence of RFCs on SHPs were both considered (MultiRFHet and RosRFCHet), acceptable accuracies (Nash-Sutcliffe efficiency or NSE ≥ 0.58 for MultiRFHet and > 0.17 for RosRFCHet) were achieved in simulating the soil water storage (SWS) variation. Since the MultiRFHet required the calibration by the observed soil water retention data, RosRFCHet was a good alternative in the simulation. On the contrary, the SHPs acquired without considering neither the RFCs nor the spatial variation yielded unacceptable simulation results (NSE general < 0). This demonstrated that spatial variation and influence of RFCs on SHPs should be included in the SWS simulation on this stony-soil hillslope, although it had limited effect in subsurface flow simulation.