One embodiment provides a method comprising receiving metadata comprising univariate time series data for each variable of a multivariate time series. The method comprises, for each variable of the multivariate time series, applying a hybrid and hierarchical model selection process to select an anomaly detection model suitable for the variable based on corresponding univariate time series data for the variable and covariations and interactions between the variable and at least one other variable of the multivariate time series, and detecting an anomaly on the variable utilizing the anomaly detection model selected for the variable. Based on each anomaly detection model selected for each variable of the multivariate time series, the method further comprises performing ensemble learning to determine whether the multivariate time series is anomalous at a particular time point.