[1] This paper evaluates the use of field data on the spatial variability of snow water equivalent (SWE) to guide the design of distributed snow models. An extensive reanalysis of results from previous field studies in different snow environments around the world is presented, followed by an analysis of field data on spatial variability of snow collected in the headwaters of the Jollie River basin, a rugged mountain catchment in the Southern Alps of New Zealand. In addition, area-averaged simulations of SWE based on different types of spatial discretization are evaluated. Spatial variability of SWE is shaped by a range of different processes that occur across a hierarchy of spatial scales. Spatial variability at the watershed-scale is shaped by variability in near-surface meteorological fields (e.g., elevation gradients in temperature) and, provided suitable meteorological data is available, can be explicitly resolved by spatial interpolation/extrapolation. On the other hand, spatial variability of SWE at the hillslope-scale is governed by processes such as drifting, sloughing of snow off steep slopes, trapping of snow by shrubs, and the nonuniform unloading of snow by the forest canopy, which are more difficult to resolve explicitly. Subgrid probability distributions are often capable of representing the aggregate-impact of unresolved processes at the hillslope-scale, though they may not adequately capture the effects of elevation gradients. While the best modeling strategy is case-specific, the analysis in this paper provides guidance on both the suitability of several common snow modeling approaches and on the choice of parameter values in subgrid probability distributions.

Abstract. This paper describes the design and use of a recursive ensemble Kalman filter (REnKF) to assimilate streamflow data in an operational flow forecasting system of seven catchments in New Zealand. The REnKF iteratively updates past and present model states (soil water, aquifer and surface storages), with lags up to the concentration time of the catchment, to improve model initial conditions and hence flow forecasts. We found the REnKF overcame instabilities in the standard EnKF, which were associated with the natural lag time between upstream catchment wetness and flow at the gauging locations. The forecast system performance was correspondingly improved in terms of Nash-Sutcliffe score, persistence index and bounding of the measured flow by the model ensemble. We present descriptions of filter design parameters and explanations and examples of filter behaviour, as an information source for other groups wishing to assimilate discharge observations for operational forecasting.

No abstract

No abstract

No abstract