Anomalies and extremes are often detected using univariate peak-over-threshold approaches in the geoscience community. The Earth system is highly multivariate. We compare eight multivariate anomaly detection algorithms and combinations of data preprocessing. We identify three anomaly detection algorithms that outperform univariate extreme event detection approaches. The workflows have the potential to reveal novelties in data. Remarks on their application to real Earth observations are provided.