An agricultural watershed, Cole Farm, was established as the newest of the three subcatchments in the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO) in 2017. The catchment contains mostly pasture and crops, with a small portion of deciduous forest. The observations in Cole Farm afford an opportunity to test the spatially distributed land surface hydrologic model, Flux-PIHM, in farmland for the first time. In this study, we calibrated the model to only discharge and groundwater level observations at Cole Farm, but it’s able to capture the variations and magnitudes of soil moisture, latent heat (LE) and sensible heat (H) fluxes. Modeled soil moisture on the ridge top matched the observations well, but modeled soil moisture in the mid-slope differed from observations likely due to the existence of fragipan in the soil column. Flux-PIHM reproduced the seasonality and diurnal variations of watershed-average evapotranspiration (ET), sensible heat flux (H), though modeled ET in summer is about 25% greater than tower ET. To study the impact of land cover on hydrology, we imposed two different LAI forcings to the model: spatially distributed versus uniform LAI. Spatially distributed LAI produced higher ET and lower soil moisture in the forested part of the watershed due to higher LAI of deciduous forest in comparison to crops and pasture. But the impact of different LAI forcings on discharge was small. We further compared the water budget simulated by Flux-PIHM in the agricultural watershed (Cole Farm) to a forested watershed (Shale Hills). Flux-PIHM simulated less discharge and higher transpiration and bare soil evaporation in the Cole Farm watershed relative to Shale Hills watershed. Our work shows that with a few key observations, Flux-PIHM can be calibrated to simulate agricultural watershed hydrology, but spatially distributed LAI and soils data are needed to capture the spatial variations in soil moisture and ET.