A new statistical process control method to monitor and diagnose bivariate normal mean vectors and covariance matrices simultaneously