Recent contamination events in water distribution networks (WDNs) suggest the need for resilience strategies to address the uncertainty of contamination intrusion. Resilience-based decision-making requires a resilience evaluation process. However, the suggestions on resilience options can vary with the definition of the system's functionality considered in the resilience evaluation. Thus, this study characterized resi-lience and its attributes (i.e., robustness (RS), loss rate (LR), recovery rate (RR), failure duration (FD), and recovery completeness (RC)) in multiple functional dimensions for a WDN to contamination events. The resilience evaluation was performed using a resilience measure based on a time-dependent functionality variation during contamination events. This study considered the functionalities in four dimensions: the contaminated node (CN), compromised demand (CD), biodegradable dissolved organic carbon (BDOC) concentration (BC), and mass con-sumed (MC). The hydraulic and water quality models were simulated using EPANET-MSX to evaluate the functionality variation under bacterial intrusion events. The results noted that the resilience levels significantly varied with the functional dimensions, relatedness, and contamination event conditions. The results also identified that different profiles of resilience attributes could characterize similar levels of multi-dimensional resilience. The findings suggest insights on incorporating the concept of multi-dimensional resilience in decision-making processes, therefore achieving the goal of improving the overall resilience of the system in diverse aspects of functionality.