A system identifies multiple causal anomalies in a power plant having multiple system components. The system includes a processor. The processor constructs an invariant network model having (i) nodes, each representing a respective system component and (ii) invariant links, each representing a stable component interaction. The processor constructs a broken network model having (i) the invariant network model nodes and (ii) broken links, each representing an unstable component interaction. The processor ranks causal anomalies in node clusters in the invariant network model to obtain anomaly score results. The processor generates, using a joint optimization clustering process applied to the models, (i) a model clustering structure and (ii) broken cluster scores. The processor performs weighted fusion ranking on the anomaly score results and broken cluster scores, based on the clustering structure and implicated degrees of severity of any abnormal system components, to identify the multiple causal anomalies in the power plant.