Modelling heteroscedasticty of streamflow times series