Detecting dominant changes in irregularly sampled multivariate water quality data sets