Predicting long-term streamflow variability in moist eucalypt forests using forest growth models and a sapwood area index

A major challenge in surface hydrology involves predicting streamflow in ungauged catchments with heterogeneous vegetation and spatiotemporally varying evapotranspiration (ET) rates. We present a top-down approach for quantifying the influence of broad-scale changes in forest structure on ET and hence streamflow. Across three catchments between 18 and 100 km² in size and with regenerating <i>Eucalyptus regnans</i> and <i>E. delegatensis</i> forest, we demonstrate how variation in ET can be mapped in space and over time using LiDAR data and commonly available forest inventory data. The model scales plot-level sapwood area (SA) to the catchment-level using basal area (BA) and tree stocking density (N) estimates in forest growth models. The SA estimates over a 69 year regeneration period are used in a relationship between SA and vegetation induced streamflow loss (L) to predict annual streamflow (Q) with annual rainfall (P) estimates. Without calibrating <i>P, BA, N, SA,</i> and <i>L</i> to <i>Q</i> data, we predict annual Q with R2 between 0.68 and 0.75 and Nash Sutcliffe efficiency (NSE) between 0.44 and 0.48. To remove bias, the model was extended to allow for runoff carry-over into the following year as well as minor correction to rainfall bias, which produced R2 values between 0.72 and 0.79, and NSE between 0.70 and 0.79. The model under-predicts streamflow during drought periods as it lacks representation of ecohydrological processes that reduce L with either reduced growth rates or rainfall interception during drought. Refining the relationship between sapwood thickness and forest inventory variables is likely to further improve results.