MODELIZACIÓN DEL FLUJO TRANSITORIO EN SUELOS FRANCO-LIMOSOS UTILIZANDO HYDRUS-1D

C. Rubio1, P. Llorens1 y M. Th. Van Genuchten2

1 Instituto de Ciencias de la Tierra Jaume Almera (CSIC), Lluís Solé i Sabarís s/n, 08028 Barcelona, España, crubio@ija.csic.es
2 George E. Brown Jr. Salinity Laboratory, USDA-ARS, 450 Big Springs Rd. Riverside, CA, USA.

RESUMEN. El modelo Hydrus-1D (Simunek et al. 1998a), basado en la ecuación de Richards, se ha utilizado para estudiar los flujos transitorios en profundidad, de una terraza de cultivo abandonada con una cubierta vegetal herbácea (cuencas experimentales de Vallcebre, NE España). Con dicho modelo, se ha realizado un ejercicio de calibración y validación utilizando potenciales mátricos y humedades volumétricas del perfil. Los resultados obtenidos han mostrado un buen ajuste a los datos de tensiometría medidos en campo cuando el perfil se humecta, para la profundidad de -20 cm (RMSE=1kPa en validación), mientras que para el resto del perfil se observa una reacción más lenta del modelo. Los contenidos hídricos simulados sobreestiman ligeramente los medidos. Cuando la demanda evapotranspirativa es elevada, el modelo sobreestima ligeramente los potenciales mátricos y los contenidos hídricos, sobre todo para los niveles superficiales del perfil. Estas diferencias pueden relacionarse con procesos de contracción-expansión de este tipo de suelos franco limosos.

ABSTRACT. The Hydrus-1D software package (Simunek et al. 1998a), based on the Richards equation, was used to study transient flow in an abandoned terraced area under grass. The field site was located in the Vallcebre research catchments (Eastern Pyrenees, North-Eastern Spain). A split test was carried out to calibrate and evaluate the model using observed matric potential and water content profiles. Results obtained showed a good agreement with the field data during wetting at 20cm depth (RMSE=1kPa for the validation), but a slow reaction of the simulated potentials with depths. Observed water contents and pressure heads were overestimated during periods of high evaporative demand. Some of the results may have been caused by the presence of shrinkage/swelling processes in the silt loam soil, which were only partially accounted for in the model.