Development of Failure Cause–Impact–Duration (CID) Plots for Water Supply and Distribution System Management

: Understanding the impact and duration (consequences) of di ﬀ erent component failures (cause) in a water supply and distribution system (WSDS) is a critical task for water utilities to develop e ﬀ ective preparation and response plans. During the last three decades, few e ﬀ orts have been devoted to developing a visualization tool to display the relationship between the failure cause and its consequences. This study proposes two visualization methods to e ﬀ ectively show the relationship between the two failure entities: A failure cause–impact–duration (CID) plot, and a bubble plot. The former is drawn for an e ﬀ ective snapshot on the range (extent) of failure duration and the impact of di ﬀ erent failures, whereas the latter provides failure frequency information. A simple and practical failure classiﬁcation system is also introduced for producing the two proposed plots e ﬀ ectively. To verify the visualization schemes, we collected records of 331 WSDS component failures that occurred in South Korea between 1980 and 2018. Results showed that (1) the proposed CID plot can serve as a useful tool for identifying most minor and major WSDS failures, and (2) the proposed bubble plot is useful for determining signiﬁcant component failures with respect to their failure consequences and occurrence likelihoods. to medium-sized districts). Another form of the CID plot (i.e., ID plot) can be drawn to show the impact and duration of individual failure (cause) in a single plot. Therefore, a total of nine ID plots are produced with identical x - and y -axis scales for e ﬀ ective comparison.


Introduction
A water supply and distribution system (WSDS) consists of various components (e.g., water treatment plant (WTP), pump station (PS), reservoir, pipes) to deliver the required quantity of water from source to customers with acceptable pressure and water quality [1]. WSDS component failures can result in water service outages, while the impact area and duration differ for each failure event [2][3][4]. For example, the segment (i.e., a subsection of the system) delineated by the surrounding valves should be closed (segment isolation) during replacement of a broken pipe [5]. A prolonged suspension of a WTP's operation would cause a service interruption across an entire city area [6]. Generally, the former case is recovered within several hours, whereas the latter case requires several days. On the other hand, component failures of the same type tend to have similar consequences despite individual failures have different timing and location [7,8]. For example, it is difficult that a local distribution pipe failure causes a service interruption for more than two weeks. Because the ultimate goal of water utilities is classification system is also introduced to produce the two proposed plots effectively. To verify the visualization schemes, we collected the records of 331 WSDS failures (WTP, reservoir, PS, water main, distribution pipe, etc.) that occurred in South Korea between 1980 and 2018.

Methodologies
This section provides the details of the proposed failure classification system, failure CID plot, and CID bubble plot.

Failure Classification
The WSDS failure cause classification system introduced here first classifies "causes" into internal and external triggers ( Table 1). Note that WSDS failure here is indicating any incidents that lead to water service interruption. The ultimate cause of an internal trigger originates within the system, e.g., the malfunction of a pumping unit in a WTP or a pipe burst. On the other hand, the cause of an external trigger comes from outside of the system (e.g., natural disasters such as an earthquake, intentional contamination, etc.). Further divided into facility and pipe triggers, internal triggers comprise a total of seven types of WSDS failure triggers: Raw water pumping facility (RWPF), WTP, reservoir (RES), PS, water supply pipe (WSP), distribution main pipe (DMP), and distribution pipe (DP) ( Table 1). The first four triggers occur in a WSDS facility that generally occupies a large space in a building structure with associated subcomponents to support its function. Therefore, the failure triggers can be caused by deterioration and malfunction of subcomponents, etc [18,19]. However, the last three are WSDS failure triggered by the pipes that occur in a line-shaped component of a WSDS that delivers water from one location to another. Pipe burst is the most common cause of this kind of failure [18,19]. In this study, external failures are either classified as a water quality failure (WQF) or a natural disaster (ND). Red water (mainly caused by the release of corrosion products from pipes [20]) and turbid water (caused by particles suspended or dissolved in service water) are examples of WQF. Note that the aforementioned WSDS component failures specially by catastrophic disaster events (e.g., the ceasing of a submerged pump station due to typhoon flooding) are classified in the category of ND, and represent the minority of failure cases. This is because a component failure from earthquakes or typhoons should be separately considered with that occurred under normal operating conditions. While the frequency characteristics of the latter tend to follow natural randomness by system internal issues (e.g., natural aging), the former is the secondary failure whose characteristics are controlled by system external forces.

Visualization Methods
This study introduces two 2-dimensional plots (i.e., failure CID plot and bubble plot) to effectively display the relationship between the cause and consequences of WSDS component failures. Water utilities can use and refer to the proposed plots to prepare for and respond to failure conditions. The most significant advantage of the visualization tool is its easiness of communication with the non-professional stakeholders (i.e., policymakers, practitioners, lawyers, etc.) which helps them to understand the needs of the system to determine investment and financing options for the further planning and management of the WSDS. Successful visualization will provide a useful comparison among the failure triggers and also help to identify the most vulnerable or critical failure trigger that requires close scrutiny. It is worthwhile to highlight that the visualization tool can also support risk analysis related to failure probability, but it is beyond the scope of this study.
This subsection describes the characteristics and steps of the plots. First, WSDS component failure records should be collected from water utilities that include the failure consequence data, i.e., failure impact and duration. Data collection should be limited in the extent to the area of interest. Failure impact data is generally recorded either by the number of households affected by failure (the number of households out of service) or the size of the administrative districts out of service. In the latter case, the number of affected households should be estimated. The failure duration data represents the period of water service outage due to failure (e.g., hours, days). Then, the cause of a failure event is classified with respect to the categorization system in Table 1. Therefore, each failure is coded as one of the nine failures in Table 1.

Failure CID Plot
The proposed failure CID plot is a 2-dimensional plot with the failure duration along the x-axis and the failure impact along the y-axis ( Figure 1). Figure 1 illustrates an example of the CID plot, which is only for demonstration purpose and not drawn with the actual data. The cause of failure is displayed in the middle of a box subplot on the 2-dimensional plot bounded by the minimum and maximum values of the failure duration and impact of each failure classification in Table 1. Note that, while each failure event of the same type (cause) of failure can have different consequences, it is more effective to use such area plots to show the overall range of the consequences of failure. For example, as shown by the light gray box labeled DP in the lower-left corner of Figure 1, the range of failure duration for DP failure is from 1 to 12 h, whereas its impact varies from 10 to approximately 1000 households (small to medium-sized districts). Another form of the CID plot (i.e., ID plot) can be drawn to show the impact and duration of individual failure (cause) in a single plot. Therefore, a total of nine ID plots are produced with identical xand y-axis scales for effective comparison. them to understand the needs of the system to determine investment and financing options for the further planning and management of the WSDS. Successful visualization will provide a useful comparison among the failure triggers and also help to identify the most vulnerable or critical failure trigger that requires close scrutiny. It is worthwhile to highlight that the visualization tool can also support risk analysis related to failure probability, but it is beyond the scope of this study. This subsection describes the characteristics and steps of the plots. First, WSDS component failure records should be collected from water utilities that include the failure consequence data, i.e., failure impact and duration. Data collection should be limited in the extent to the area of interest. Failure impact data is generally recorded either by the number of households affected by failure (the number of households out of service) or the size of the administrative districts out of service. In the latter case, the number of affected households should be estimated. The failure duration data represents the period of water service outage due to failure (e.g., hours, days). Then, the cause of a failure event is classified with respect to the categorization system in Table 1. Therefore, each failure is coded as one of the nine failures in Table 1.

Failure CID Plot
The proposed failure CID plot is a 2-dimensional plot with the failure duration along the x-axis and the failure impact along the y-axis ( Figure 1). Figure 1 illustrates an example of the CID plot, which is only for demonstration purpose and not drawn with the actual data. The cause of failure is displayed in the middle of a box subplot on the 2-dimensional plot bounded by the minimum and maximum values of the failure duration and impact of each failure classification in Table 1. Note that, while each failure event of the same type (cause) of failure can have different consequences, it is more effective to use such area plots to show the overall range of the consequences of failure. For example, as shown by the light gray box labeled DP in the lower-left corner of Figure 1, the range of failure duration for DP failure is from 1 to 12 h, whereas its impact varies from 10 to approximately 1000 households (small to medium-sized districts). Another form of the CID plot (i.e., ID plot) can be drawn to show the impact and duration of individual failure (cause) in a single plot. Therefore, a total of nine ID plots are produced with identical x-and y-axis scales for effective comparison.  The proposed CID plot can be divided into four quadrants by setting a threshold for each of the failure durations and impacts beyond which water utilities should adopt system-wide emergency and recovery actions ( Figure 1). Quadrant 1 (Q1) in the upper right corner of the CID plot includes the most critical WSDS component failures that have a significant consequence on the system in general with more than one week of service interruption to a number of households greater than 10,000 (at least the majority of a small city is affected) (Figure 1). For example, RWPF failure would cause water service failure to an entire city having a single raw water intake location. In such failures, an emergency interconnection pipe should be operated to receive water from neighboring cities with sufficient resources [7,21]. While the most frequently occurring daily failures (e.g., DP) are located in quadrant 4 (Q4), quadrant 2 (Q2) shows failures with great impact, but that can be recovered from within days. Those positioned in quadrant 3 (Q3) are prolonged failures with a minor impact on the system. In 2010, the breakage of a 150 mm WSP installed under the sea to supply 70 customers living on an island in Korea resulted in a water service outage for 220 h (about 9 d). System-specific thresholds should be determined by the water utility considering the system conditions and management target.
Water utility with a high management standard will have the cross of the threshold lines at the lower-left corner of the CID plot with the largest quadrant area in Q1. On the other hands, a tolerant utility would position the cross at the upper-right corner in the plot with a large Q4 area.

CID Bubble Plot
Multiple plots should be used to capture all failure-characteristic information laid in the historical records. Although the proposed CID plot provides an effective snapshot of the ranges (extent) of failure duration and the impact of different failures, failure frequency information is difficult to incorporate and visualize in a plot. Therefore, we also propose a CID bubble plot in which each failure is drawn as a bubble on the 2-dimensional plot with the same axis as the CID plot, where each bubble has a radius proportional to the frequency of the corresponding failure ( Figure 2). For example, based on Figure 2, DMP and DP failures are the most frequently occurring, whereas there are few RWPF failures. The center of the bubble is positioned at either the average or maximum values for failure duration and impact from the corresponding failure in the historical failure records collected which can be expressed as follows: where (x bubble and y bubble ) j indicate the x and y coordinates of the jth bubble (failure cause of jth categorization in Table 1) in the CID bubble plot; FD ij and FI ij are the failure duration and impact of the jth failure's ith event (e.g., j = 1 for RWPF failure in Table 1); nj is the total number of events of the jth failure recorded. Note that nj is more likely to be different for different failures in Table 1 and ∨ is the "Or" sign. Therefore, compared with the CID plot, the CID bubble plot provides rough information on failure duration and impact but has the advantage of providing relative failure frequency information.

Data
In order to verify the proposed CID and bubble plots, we collected the records of 331 WSDS component failures that occurred in South Korea between 1980 and 2018. Note that nation-wide failure data was collected and used to produce the plots because of data scarcity. Thus, the detailed description of the WSDS characteristics were not available. This is only for demonstration purposes, and one can use a single utility failure data, if sufficient, to create the plots suitable for site-specific decision-making tools.
Most water utilities in Korea stipulate filling out a failure event response sheet in which the date, cause, response action is taken, impact and duration (consequences), lessons, and figures should be recorded in the event of a component failure. It is worthwhile to highlight that the administration in South Korea generally records impact on the household unit (e.g., the number of households out of service) for simplicity purpose. The two largest water providers in Korea, the Office of Waterworks of Seoul Metropolitan Government (about 10 million customers) and K-water, regularly publish a WSDS failure casebook, i.e., the collection of the response sheets. In addition, the aforementioned failure cause and consequences information are often found in news articles, especially when the failure consequence is relatively significant. Therefore, this study used the historical failure data collected from the casebooks (1980-2015) and news articles (2008-2018). Table 2 summarizes several statistics for the collected WSDS component failures. DP and DMP are the most common failures in WSDSs in Korea, covering 31.4 and 29.6% of the total 331 failures, respectively. The longest average and maximum failure durations (41 and 120 h, respectively) were RWPF failures. However, the maximum failure impact, which is the maximum number of households out of service, was observed in a WTP failure (0.623 million households). Note that the average number of persons in a household in Korea has remained steady at 2.5 since 2015. A more comprehensive and effective understanding of the failure data can be achieved with the proposed CID plots in the next section. Therefore, compared with the CID plot, the CID bubble plot provides rough information on failure duration and impact but has the advantage of providing relative failure frequency information.

Data
In order to verify the proposed CID and bubble plots, we collected the records of 331 WSDS component failures that occurred in South Korea between 1980 and 2018. Note that nation-wide failure data was collected and used to produce the plots because of data scarcity. Thus, the detailed description of the WSDS characteristics were not available. This is only for demonstration purposes, and one can use a single utility failure data, if sufficient, to create the plots suitable for site-specific decision-making tools.
Most water utilities in Korea stipulate filling out a failure event response sheet in which the date, cause, response action is taken, impact and duration (consequences), lessons, and figures should be recorded in the event of a component failure. It is worthwhile to highlight that the administration in South Korea generally records impact on the household unit (e.g., the number of households out of service) for simplicity purpose. The two largest water providers in Korea, the Office of Waterworks of Seoul Metropolitan Government (about 10 million customers) and K-water, regularly publish a WSDS failure casebook, i.e., the collection of the response sheets. In addition, the aforementioned failure cause and consequences information are often found in news articles, especially when the failure consequence is relatively significant. Therefore, this study used the historical failure data collected from the casebooks  and news articles (2008-2018). Table 2 summarizes several statistics for the collected WSDS component failures. DP and DMP are the most common failures in WSDSs in Korea, covering 31.4 and 29.6% of the total 331 failures, respectively. The longest average and maximum failure durations (41 and 120 h, respectively) were RWPF failures. However, the maximum failure impact, which is the maximum number of households out of service, was observed in a WTP failure (0.623 million households). Note that the average number of persons in a household in Korea has remained steady at 2.5 since 2015. A more comprehensive and effective understanding of the failure data can be achieved with the proposed CID plots in the next section.

CID Plot
Two CID and nine ID plots were drawn using the 331 component failure records (Figures 3 and 4, respectively). The former is the CID plot with box subplots bounded by the minimum and maximum (Figure 3a, CID-1) and the average and maximum (Figure 3b, CID-2) of the duration and impact of nine WSDS failure triggers as classified in Table 1. A single 2-dimensional ID plot is drawn for each of the nine triggers with identical xand y-axis (failure duration and impact, respectively) scales for effective comparison (Figure 4).

CID Plot
Two CID and nine ID plots were drawn using the 331 component failure records (Figures 3 and  4, respectively). The former is the CID plot with box subplots bounded by the minimum and maximum (Figure 3a, CID-1) and the average and maximum (Figure 3b, CID-2) of the duration and impact of nine WSDS failure triggers as classified in Table 1. A single 2-dimensional ID plot is drawn for each of the nine triggers with identical x-and y-axis (failure duration and impact, respectively) scales for effective comparison (Figure 4).    Table 1.   Table 1.  Table 1.  As shown in Figure 3a, the CID-1 plot has box subplots overlapping each other because each of the nine failure triggers has a wide range of failure duration and impact (Table 1). This is due to the use of nation-wide data in this study. Opposite to what would be expected, components located upstream of WSDS were not always positioned in the upper right corner of the CID plot (Q1) constructed from the 331 failures considered in this study. For example, the maximum failure duration of RES (approximately a half-day) was shorter than that of DMP (10 d), which resulted in a RES box with narrow width being positioned in the left corner of the CID-1 plot (Figure 3a), and a big DMP box covering every quadrant (Figure 3a). A better CID-1 plot would be obtained if a sufficient number of failure records collected in an individual city are used.
The CID-2 plot provides a better view of an individual failure's overall consequence range, as the coordinate of the left edge of each box is equal to the average failure duration, whereas that of the lower edge is the average number of households experiencing water service outage (Figure 3b). While the RES box shrunk toward Q2, the CID-2 plot clearly shows that the DP and ND failures were the most minor failures compared with other failures considered in this study, because their boxes are located near Q4. Generally, only a local district in the downstream of a broken DP would be isolated for recovery. Thirteen failure events classified as ND (external failure cause) were WSDS component failures directly from typhoons and flooding with insignificant failure consequences. It was clearly confirmed from CID-2 plot that the component failures by ND considered in this study were not critical ones. The operation of an emergency readiness team during the disasters would help limit any further increase in the failure consequences, but the box would expand toward Q1 if seismic failures were considered (i.e., multiple component failures) [12] (Figure 3b). Note that other boxes, e.g., RWPF, WTP, PS, WSP, WQF, also overlapped in the CID-2 plot.
While the two CID plots can serve as useful tools to identify the most minor and major component failures (with the historical records considered in this study), the individual ID plots can provide a firm platform for detailed failure consequence comparison between the nine failures ( Figure 4). The nine individual ID plots not only present a bounded box area, but also contain the consequence data points of each failure category by which the point distribution could be identified. In this study, we use two different xand y-axis scales to better display the distribution of the data points. The xand y-axes for RWPF and WTP are consistently set to range from 0 to 120 h (10 d) and from 0 to 0.8 M households (e.g., the entire city), respectively (Figure 4a,b). Note that the ranges are 0 to 80 h (6.7 d) and from 0 to 0.1 M households for other components (Figure 4c-i). Data points beyond these ranges were assumed to be the maximum x and y-axis values (e.g., the maximum failure impact of DMP was 0.505 M households, Table 2), which was represented as 0.1 M households (the maximum y-axis value in Figure 5f). While some of the failures have few data points (Figure 4a,h,i), most of the box area (its width and height) was determined under the bias of a few extreme data points.
The shape of some of the boxes (i.e., the relative width and height) in the individual ID plots (Figure 4c,g,i) is similar to those included in the CID-2 plot (e.g., RES, DP, and ND boxes in Figure 3b). Most data points were gathered in the lower left corner of the WSP, DMP, and DP (i.e., pipe failures) ID plots (Figure 4e-g, respectively), indicating that the majority of the failure events were responded to and managed well. Therefore, water utility should focus on developing special strategies to promptly and effectively respond to few extreme failure events which cause a relatively long failure duration and great impact compared to other failures in the same category (Figure 4).

CID Bubble Plot
By using the failure data, two CID bubble plots were drawn: One with the center of the bubble positioned at the maximum (Figure 5a, ID-1), and the other with the center at the average value of the failure duration and impact (Figure 5b, ID-2). Because the bubbles are shifted toward the lower left corner of the CID plot in the latter case (the maximum is bigger than the average), the two bubble plots have different xand y-axis scales for better visualization. Note again that the radius of each bubble is proportional to each failure's frequency and proportion (occurrence likelihood) in Table 2.
The most critical failure can be determined from the two ID plots based on the radius and location of the bubble: The bubble with a large radius (high occurrence likelihood) located in the upper right corner (Q1) indicates the most critical one ( Figure 5). It was confirmed from the ID-1 plot that DMP and DP are the most frequent types of failures while the DMP also has the longest maximum failure duration (Figure 5a). However, the ID-2 plot indicates that their expected (average) failure consequences are smaller than other failures (Figure 5b). Based on the ID-2 plot, RWPF failure is the most critical failure but has a low likelihood because a small bubble is located in Q1. Note that only four failure events are of the RWPF type among the total 331 failures (1.2%) considered in this study (Table 2). Considering both ID plots, RWPF, WTP, WSP, and DMP are the critical components in WSDSs in Korea for which failure can cause significant consequences. It was verified that the CID bubble plots are useful for comparing and selecting WSDS component failures with respect to their failure consequences and occurrence likelihood. The most critical failure can be determined from the two ID plots based on the radius and location of the bubble: The bubble with a large radius (high occurrence likelihood) located in the upper right corner (Q1) indicates the most critical one ( Figure 5). It was confirmed from the ID-1 plot that DMP and DP are the most frequent types of failures while the DMP also has the longest maximum failure duration (Figure 5a). However, the ID-2 plot indicates that their expected (average) failure consequences are smaller than other failures (Figure 5b). Based on the ID-2 plot, RWPF failure is the most critical failure but has a low likelihood because a small bubble is located in Q1. Note that only four failure events are of the RWPF type among the total 331 failures (1.2%) considered in this study (Table 2). Considering both ID plots, RWPF, WTP, WSP, and DMP are the critical components in WSDSs in Korea for which failure can cause significant consequences. It was verified that the CID bubble plots are useful for comparing and selecting WSDS component failures with respect to their failure consequences and occurrence likelihood.  Table 2.

Broader Impact
The goal of CID visualization is to provide the management strategy of the WSDS intuitively. For the past few decades, resilient and robust planning has become an important consideration of the WSDS planning and management purpose owing to the high uncertainty of disturbances and the consequences of component failure [10,21]. Thus, knowing the variation (or uncertainty) of the failure consequences and preparing for the variation is of paramount importance for the utility and all other stakeholders to which both CID-1 and CID-2 plots (Figure 3) can deliver an important message for the purpose. First, CID-1 is not only showing how severe the consequence was of each failure trigger but also which trigger induces the most uncertain consequence; while CID-2 shows a similar vision but focusing on more severe events. For example, DMP in CID-1 shows a wide range of failure duration and impact, which indicate the need for a better isolation strategy or the installation of isolation valves to reduce the variation of the consequence. Moreover, the prioritization of management planning can be sought based on the location of each trigger, depending on the preference of the utility. For the given example in Figure 3a, and assuming the utility is more aware of significant failure (e.g., right upper corner of each rectangular), DMP, RWPF, and WSP will be the high priority of management. In other words, close scrutiny is required for DMP, RWPF, and WSP. On the other hand, if the utility only cares about the failure impact, then WTP, DMP, and RWPF will be the high priority of management. Thus, both CID plots can serve as a reference for resilient and robust WSDS planning and management. Regarding CID bubble plots ( Figure 5), a similar conclusion is also derived by the size (radius) and the location of the bubble. Additionally, CID bubble plots are showing the frequency of the failure as well, which is not shown in CID plots.  Table 2.

Broader Impact
The goal of CID visualization is to provide the management strategy of the WSDS intuitively. For the past few decades, resilient and robust planning has become an important consideration of the WSDS planning and management purpose owing to the high uncertainty of disturbances and the consequences of component failure [10,21]. Thus, knowing the variation (or uncertainty) of the failure consequences and preparing for the variation is of paramount importance for the utility and all other stakeholders to which both CID-1 and CID-2 plots (Figure 3) can deliver an important message for the purpose. First, CID-1 is not only showing how severe the consequence was of each failure trigger but also which trigger induces the most uncertain consequence; while CID-2 shows a similar vision but focusing on more severe events. For example, DMP in CID-1 shows a wide range of failure duration and impact, which indicate the need for a better isolation strategy or the installation of isolation valves to reduce the variation of the consequence. Moreover, the prioritization of management planning can be sought based on the location of each trigger, depending on the preference of the utility. For the given example in Figure 3a, and assuming the utility is more aware of significant failure (e.g., right upper corner of each rectangular), DMP, RWPF, and WSP will be the high priority of management. In other words, close scrutiny is required for DMP, RWPF, and WSP. On the other hand, if the utility only cares about the failure impact, then WTP, DMP, and RWPF will be the high priority of management. Thus, both CID plots can serve as a reference for resilient and robust WSDS planning and management. Regarding CID bubble plots ( Figure 5), a similar conclusion is also derived by the size (radius) and the location of the bubble. Additionally, CID bubble plots are showing the frequency of the failure as well, which is not shown in CID plots.
The individual ID plots shown in Figure 4 can also be expressed as a form of the response surface, which is a heat map-like graph reflecting data with different coloring. It has been frequently used for adaptation planning by visualizing the impact of climate change for decision making purposes [22,23]. Suggesting a response surface is not exactly same with the one used for those studies, but we can adapt the heatmap format to re-create the individual ID plots. However, for this type of visualization, more data is better since the color variation will be better illustrated. For this reason, only WSP, DMP, and DP ID response surface is illustrated in Figure 6. The individual ID plots shown in Figure 4 can also be expressed as a form of the response surface, which is a heat map-like graph reflecting data with different coloring. It has been frequently used for adaptation planning by visualizing the impact of climate change for decision making purposes [22,23]. Suggesting a response surface is not exactly same with the one used for those studies, but we can adapt the heatmap format to re-create the individual ID plots. However, for this type of visualization, more data is better since the color variation will be better illustrated. For this reason, only WSP, DMP, and DP ID response surface is illustrated in Figure 6. As shown in Figure 6, the color changes as the probability of the given failure impact and duration changes. Note that darker color is a higher probability of occurrence. Figure 6 implies the same message as Figure 4 does, but enhances the visualization of Figure 4 if the sufficient amount of data can be acquired. In addition, since this type of visualization is frequently used for adaptation planning, this visualization approach can also provide the high probability of an event that utility should aware of or probability of high failure consequence. Especially providing probability of high failure consequence can help the utility to determine whether investment for the high failure consequence is worthwhile or not. Moreover, the quantified risk of each failure consequence can be provided if the probability can be fully investigated.

Summary and Conclusions
This study proposed CID and bubble plots to effectively display the relationship between the failure cause and consequences (duration and impact) of WSDSs and to help communicate with nonprofessional stakeholders (i.e., policymakers, practitioners, lawyers, etc.). The proposed CID plot provides an effective snapshot of the ranges (extents) of failure duration and the impact of different failures, whereas the CID bubble plot contains failure frequency information. This study also introduced a simple and practical failure classification system that categorized WSDS failure into nine component failures (RWPF, WTP, RES, PS, WSP, DMP, DP, WQF, and ND). To verify the proposed visualization methods, we produced the proposed plots using the historical records of 331 WSDS component failures, which were collected in South Korea from 1980 to 2018. DP and DMP were the most frequent failures in Korean WSDSs, comprising 31.4 and 29.6% of the total 331 failures, respectively.
First, two CID plots were drawn using the 331 failure records: (1) One with box subplots bounded by the maximum and minimum of the failure duration and impact, and (2) the other with subplots bounded by the average and maximum values. While most box subplots overlapped in the former, the latter plot clearly identified the most minor and major component failures. Then, nine individual ID plots were drawn, each of which showed the relationship between the failure duration and impact of an individual failure cause. It was confirmed from the ID plots that most of the box area (its width and height) was determined under the bias of a few extreme data points. The individual ID response surface is also suggested as an alternative form of the individual ID plot. Lastly, the CID bubble plots were produced with the center of the bubbles positioned at either the maximum or the average values of the failure duration and impact. The radius of the bubble indicated As shown in Figure 6, the color changes as the probability of the given failure impact and duration changes. Note that darker color is a higher probability of occurrence. Figure 6 implies the same message as Figure 4 does, but enhances the visualization of Figure 4 if the sufficient amount of data can be acquired. In addition, since this type of visualization is frequently used for adaptation planning, this visualization approach can also provide the high probability of an event that utility should aware of or probability of high failure consequence. Especially providing probability of high failure consequence can help the utility to determine whether investment for the high failure consequence is worthwhile or not. Moreover, the quantified risk of each failure consequence can be provided if the probability can be fully investigated.

Summary and Conclusions
This study proposed CID and bubble plots to effectively display the relationship between the failure cause and consequences (duration and impact) of WSDSs and to help communicate with non-professional stakeholders (i.e., policymakers, practitioners, lawyers, etc.). The proposed CID plot provides an effective snapshot of the ranges (extents) of failure duration and the impact of different failures, whereas the CID bubble plot contains failure frequency information. This study also introduced a simple and practical failure classification system that categorized WSDS failure into nine component failures (RWPF, WTP, RES, PS, WSP, DMP, DP, WQF, and ND). To verify the proposed visualization methods, we produced the proposed plots using the historical records of 331 WSDS component failures, which were collected in South Korea from 1980 to 2018. DP and DMP were the most frequent failures in Korean WSDSs, comprising 31.4 and 29.6% of the total 331 failures, respectively.
First, two CID plots were drawn using the 331 failure records: (1) One with box subplots bounded by the maximum and minimum of the failure duration and impact, and (2) the other with subplots bounded by the average and maximum values. While most box subplots overlapped in the former, the latter plot clearly identified the most minor and major component failures. Then, nine individual ID plots were drawn, each of which showed the relationship between the failure duration and impact of an individual failure cause. It was confirmed from the ID plots that most of the box area (its width and height) was determined under the bias of a few extreme data points. The individual ID response surface is also suggested as an alternative form of the individual ID plot. Lastly, the CID bubble plots were produced with the center of the bubbles positioned at either the maximum or the average values of the failure duration and impact. The radius of the bubble indicated the relative occurrence likelihood of corresponding failure. RWPF, WTP, WSP, and DMP were the most critical components in Korean WSDSs based on the bubble plots.
It was confirmed that (1) the CID plots can serve as a useful tool to identify the most minor and major WSDS failures, and (2) the bubble plots are useful in determining the most significant failures with respect to their failure consequences and occurrence likelihoods. Collecting a sufficient number of historical failure records for various components will help in constructing more reliable CID and bubble plots.
This study has several limitations that future research must address. First, we collected records for component failures that occurred across multiple WSDSs in Korea. A water utility may want to limit their data collection to failures that occurred in their WSDS to have a system-specific relationship between failure cause and consequences. To that end, a sufficient number of failure records should be collected, with each failure event classified following the proposed classification system, and the failure duration and impact data properly documented. In the worst case, an estimation of the information should be recorded. Second, the proposed CID and bubble plots can be developed and compared under various WSDS conditions and characteristics. For example, it is worthwhile to compare plots with and without emergency pipes between district metering areas. Third, a 3-dimensional version of the CID and bubble plots could be constructed by displaying the temporal changes of the duration and impact over months/years along the z-axis. By indicating the major changes in WSDSs (e.g., parallel pipe construction in WSP lines, development of industrial complex in part of the system) in the z-axis (with arrow and text), these plots will effectively show the transition of the failure consequences upon system condition changes. Lastly, probability can be assigned to the CID plots to guide the water utility to conduct a risk analysis. The quantified risk of failure impact and duration can also be calculated once the probability is fully established.