Water distribution networks with dynamically adaptive connectivity offer greater operational flexibility. While this strategy has demonstrated improvements in pressure management and network resiliency, further research is needed to better understand its impact on water quality dynamics. This paper investigates the short-term variability of disinfectant residuals in a real-world distribution network operated with dynamic connectivity. We simulate water quality dynamics under two control configurations with pressure control and automatic flushing valve operations. Our simulation results inform the development of flow variability metrics to reveal the relationship between changing hydraulic conditions and increased water quality dynamics. These metrics can (i) improve observability by supporting the placement of additional water quality monitoring locations and (ii) enhance controllability by enabling the formulation of optimization problems that incorporate hydraulic surrogates for modelling water quality. Furthermore, we validate the identified regions of increased water quality dynamics using continuous disinfectant data from a large-scale experimental programme. Our findings emphasize the need for further analytical and experimental research to manage water quality in distribution networks with dynamically adaptive connectivity and hydraulic control.

  • Investigated the impact of dynamically adaptive connectivity on water quality in a real-world distribution network.

  • Proposed flow variability metrics to identify areas of the network with increased disinfectant variability.

  • Validated model simulations using time series analysis of continuous water quality data from a real-world experiment.

  • Highlighted the need for further research to improve the observability and controllability of water quality in dynamically adaptive networks.

District metered areas (DMAs) have proven highly effective for managing leakage and pressure in water distribution networks. However, traditional single-feed DMAs can present significant operational challenges such as reduced network resilience and water quality deterioration. In response to these challenges, Wright et al. (2014, 2015) developed a novel concept of adaptive DMAs, which can dynamically adapt network connectivity (topology) and hydraulic control to accommodate a range of operational conditions and objectives. Although this operational strategy has been shown to improve pressure management and network resiliency (Wright et al. 2015; Giudicianni et al. 2020; Bui et al. 2022; Ulusoy et al. 2023), its impact on water quality dynamics has not been studied. Further research is therefore needed to investigate the impacts of advanced hydraulic control on the temporal variability of disinfectant residuals.

Traditional single-feed DMAs often impose artificial dead ends at their boundaries, resulting in hydraulic conditions characterized by low-flow velocities and long residence times. Consequently, dead-end pipes are generally associated with poor water quality conditions (US Environmental Protection Agency 2002; Machell & Boxall 2014; Blokker et al. 2016; Dias et al. 2017; Armand et al. 2018). Excessive disinfectant decay at dead ends can be attributed to long residence times and the influence of dispersive transport mechanisms (Abokifa et al. 2016; Abokifa et al. 2019). Moreover, the presence of low-flow conditions is associated with greater sediment accumulation, biofilm growth, and iron uptake from corrosion scaling (in metallic pipes), increasing the risk of discolouration incidents (Armand et al. 2017, 2018). While adaptive DMAs are designed to reduce dead ends, the increased changes in hydraulic conditions, specifically the variability in flow paths, could potentially introduce new water quality challenges. The main water quality challenge in adaptive DMAs is the management of source water mixing from multiple inlets, particularly when inlets have markedly different source water characteristics (e.g. disinfectant concentrations and total organic carbon). Although source water mixing has been modelled at the network level to assess the effects of new/different water supplies (Montiel et al. 2002; Ostfeld & Salomons 2004; Mala-Jetmarova et al. 2015; Szuster–Janiaczyk & Bylka 2019), its magnitude and frequency caused by the operation of adaptive DMAs at the distribution level remains unexplored. Additionally, flow path variability can create flow reversals and/or a wider range of pipe flow velocities, potentially leading to higher water age in looped sections of the network. We note that, while not explicitly considered in this work, changes in flow velocities are also linked to discolouration risk in distribution networks (Horn et al. 2003; Ahn et al. 2011; Derlon et al. 2008; Furnass et al. 2013; Boxall et al. 2023), which can in turn affect disinfectant decay behaviour.

In this paper, we investigate the relationship between changing hydraulic conditions and short-term temporal disinfectant variability within dynamically adaptive DMAs. We focus on hydraulic conditions resulting from two control configurations: (i) minimizing average zone pressure (AZP) (Wright et al. 2015) and (ii) integrating AZP with a self-cleaning capacity (SCC) control mode (Jenks et al. 2023a, 2023b). The AZP objective aims to reduce background leakage through pressure reduction, a common operational practice in UK water networks. On the other hand, the SCC objective is a relatively new strategy which aims to proactively decrease the risk of discolouuration in distribution networks. These control configurations are evaluated using a multi-feed operational network located in Bristol, UK, where ongoing experimental and analytical research is carried out by Imperial College London, Bristol Water, and Badger Meter/Analytical Technology (ATi). For each configuration, we identify areas of increased hydraulic dynamics using a set of proposed flow variability metrics computed over a 24-h simulation. These metrics aim to assess the degree of source water mixing between DMA inlets and identify other hydrodynamic conditions known to influence water quality deterioration. The flow variability metrics can also serve as hydraulic surrogates for modelling water quality dynamics in optimization problems formulated for control. Additionally, we simulate water quality dynamics in the operational network with varying inlet disinfectant concentrations to study the potential effects of source water mixing. We then validate our simulation results for the AZP control configuration, which is currently implemented in the field, by performing a time series analysis of data collected from an ongoing water quality sensing programme. The time series analysis focuses on highlighting short-term temporal (i.e. hourly) fluctuations in disinfectant data across the network.

This study represents a significant advancement in the field as it comprises an experimental programme with continuous water quality sensing at high spatial and temporal resolution. Previous water quality field studies have mainly consisted of pilot-scale programmes relying on manual collection techniques (Vasconcelos et al. 1997; Liu et al. 2015; Maleki et al. 2023) or short-term deployment of continuous reagent-based chlorine sensors (US Environmental Protection Agency 2008; Lee et al. 2010). In this work, we analyse continuous disinfectant data collected in a real-world operational network to validate our simulation results, namely the modelling of source mixing (i.e. flow paths). Our simulations and exploratory analysis of sensor data inform the design of subsequent experiments to enhance the observability and controllability of water quality dynamics in adaptive DMAs.

The remainder of this paper is structured as follows. Section 2 introduces the operational case study network, which we use to evaluate different control configurations as well as collect continuous water quality data. Section 3 outlines the paper's methodology, including model development, the proposed flow variability metrics, and time series data processing procedures. In Section 4, we present and discuss our simulation and time series analysis results. Section 5 concludes with suggestions for increasing the observability and controllability of water quality in adaptive DMAs.

Overview

The Bristol Water Field Lab (Field Lab) is an operational distribution network supplying approximately 8,000 customer connections in Bristol, UK. The network was originally comprised of two single-feed DMAs hydraulically isolated via permanently closed isolation valves (IVs). Wright et al. (2014) introduced the Field Lab as an experimental programme for dynamically aggregating the isolated DMAs through the installation of dynamic boundary valves (DBVs), which are equipped with bidirectional control capabilities. Flow modulation profiles were also derived for advanced pressure control at DMA inlets. The Field Lab has since been a case network for numerous water distribution analysis and optimization studies and a testbed for advances in pressure and water quality sensing technologies. The subsequent sections describe the different control configurations evaluated in the Field Lab (Section 2.2) and the current water quality sensing programme conducted by our research group (Section 2.3).

Control configurations

In addition to the original single-feed DMA topology, we consider two control configurations to explore water quality dynamics in the Field Lab: (i) minimizing AZP and (ii) integrating AZP and SCC operations. Figure 1 illustrates the network topology and control valve locations for each configuration.
Figure 1

Control configurations applied to the Field Lab network: (a) single-feed DMAs, (b) adaptive AZP control, (c) adaptive AZP-SCC control.

Figure 1

Control configurations applied to the Field Lab network: (a) single-feed DMAs, (b) adaptive AZP control, (c) adaptive AZP-SCC control.

Close modal

The original single-feed DMA topology is illustrated in Figure 1(a). In this setup, several IVs sectorize different areas of the network, three of which separate DMAs 1 and 2, and pressure-reducing valves (PRVs) reduce DMA inlet pressures. The adaptive AZP control configuration in Figure 1(b) replaces two of the IVs with remote-control DBVs. These DBVs permit bidirectional control between DMAs, enabling the network to dynamically adapt between an isolated topology at night and an open topology during the day (Wright et al. 2014, 2015). Additionally, existing PRVs are retrofitted with flow modulated control to minimize the network's AZP (and thereby leakage) (Wright et al. 2015). This adaptive AZP configuration is currently implemented in the Field Lab.

We also consider a new water quality control strategy aimed at increasing the number of pipes in the Field Lab experiencing self-cleaning velocities (Jenks et al. 2023a, 2023b). The SCC of the network, defined as its ability to periodically resuspend and flush accumulated material, is improved by jointly optimizing the location and settings of PRV and automatic flushing valve (AFV) actuators, the latter of which discard water from the system in a controlled manner (Figure 1(c)) (Jenks et al. 2023a). In this configuration, PRV and AFV controls switch between the existing AZP mode and a 1-h self-cleaning mode where flushing is activated to mobilize material which has potentially accumulated in pipes. Such self-cleaning controls have been shown to significantly alter flow paths in the Field Lab (Jenks et al. 2023b), particularly due to the removal of existing IVs and DBVs. This operational mode creates a more looped network configuration, potentially causing large shifts in the hydraulic balance point. In this work, the self-cleaning configuration represents advanced control enabled by adaptive DMAs to reduce the risk of discolouration. We simulate the resulting water quality dynamics using a mathematical model of the Field Lab. Planned experimental work, which leverages the analytical methods and findings of this study, aims to deploy the adaptive AZP-SCC control configuration in the Field Lab.

Water quality sensing

Water quality in the Field Lab is continuously monitored at nine locations using MetriNet Q52 multi-parameter, reagent-free water quality sensors (Analytical Technology Inc. 2022). These sensors collect data on five parameters – chlorine, pH, turbidity, conductivity, and temperature – at frequencies ranging from one sample per second to every 15 min. The data are transmitted to our research database at Imperial College London for near real-time processing and analysis. Note that we perform monthly calibrations of the MetriNet sensors using manual DPD (N,N-diethyl-p-phenylenediamine) grab samples to maintain suitable measurement accuracy for modelling the transport and decay of disinfectant across the network. Figure 2 illustrates the locations of MetriNet water quality sensors deployed in the Field Lab.
Figure 2

Field Lab water quality sensor deployment.

Figure 2

Field Lab water quality sensor deployment.

Close modal
Figure 3 shows the three MetriNet sensor installation types in the Field Lab. Six devices are placed in telemetry kiosks (Figure 3(a)), corresponding to locations with pressure control such as DMA inlets and boundary valves; three devices are installed in fire hydrant chambers (Figure 3(b)); and one device is housed within a specially constructed telemetry bollard (Figure 3(c)). The telemetry bollard aims to avoid adverse conditions associated with below-grade installations, as shown in Figure 3(b).
Figure 3

MetriNet sensor installations: (a) telemetry kiosk, (b) fire hydrant chamber, (c) telemetry bollard.

Figure 3

MetriNet sensor installations: (a) telemetry kiosk, (b) fire hydrant chamber, (c) telemetry bollard.

Close modal

The time series data used in this study were downsampled to 15-min intervals to extend battery life and ensure the longevity of data collection from the devices, which collected data from April to July 2023.

Figure 4 illustrates the research methodology of this paper. The following subsections provide further detail for each of the applied methods.
Figure 4

Map of research methodology.

Figure 4

Map of research methodology.

Close modal

Simulation

Model development

The hydraulic model of the Field Lab comprises junction nodes, variable head nodes (e.g. DMA inlets), links, and discrete 15-min time steps. Temporal patterns for variable head and demand nodes were acquired from pressure and flow sensors deployed as part of the Field Lab's ongoing experimental programme. For each time step , we solve the steady-state flow () and hydraulic head () conditions using state-of-the-art modelling software EPANET (Rossman et al. 2020), accessed via the WNTR Python package (Klise et al. 2017). As discussed in Section 2.2, we introduce control through the placement and operation of PRV, DBV, and AFV actuators. For each time step , we set (i) pressure setpoints at the downstream node of PRV links, (ii) dynamic local loss coefficients at DBV links, and (iii) flushing demands at AFV nodes. Such controls are derived from previous work which focused on optimizing hydraulic control within adaptive WDNs for AZP (Wright et al. 2015; Jenks et al. 2023b) and SCC (Jenks et al. 2023a, 2023b) objectives. Details on PRV, DBV, and AFV actuators as well as other model inputs are available in the project's GitHub repository (Jenks et al. 2024).

We use EPANET to simulate the transport and fate of water quality constituents (e.g. disinfectant and source trace). For modelling disinfectant residuals, we assume a first-order reaction relationship with bulk decay coefficient and wall decay coefficients ranging from to . These decay coefficients were calibrated to water quality conditions in the Field Lab. We set initial conditions (at ) across all links and simulated water quality dynamics over a 7-day period to ensure periodic behaviour is reached across all nodes. At each DMA inlet , we assign source concentration . Here, corresponds to the mean or worst-case (10th and 90th percentiles) of disinfectant variability conditions observed at inlet s over the period April to July 2023 (Table 1). We assume a normal distribution with standard deviation mg L−1 for modelling inlet variability from upstream mixing/storage dynamics. Since this work focuses on short-term disinfectant variability resulting from source water mixing between DMAs, we do not explicitly consider any seasonal (climatic) variability in the assigned inlet concentration patterns. Additionally, no clear seasonal trend could be discerned in the time series data collected at DMA inlets mainly due to the relatively short sensing period. We note, however, that seasonal variability is certainly an important consideration in managing disinfectant residuals, and that future work needs to consider data collected over longer time periods.

Table 1

Mean inlet source concentrations for water quality simulations

InletMean scenario (mg L−1)Worst-case scenario (mg L−1)
DMA 1 0.75 0.95 
DMA 2 0.65 0.40 
InletMean scenario (mg L−1)Worst-case scenario (mg L−1)
DMA 1 0.75 0.95 
DMA 2 0.65 0.40 

For source trace simulations, the contribution of source water is modelled as a conserved constituent (i.e. no reaction kinetics).

Flow variability

We propose three flow variability metrics to identify areas in the network with increased changes in hydraulic conditions: (i) count of pipe flow reversals; (ii) coefficient of variation (CV) of diurnal pipe flow velocities; and (iii) percentage of source water mixing. These metrics are computed using the steady-state hydraulic model of the Field Lab over a 24-h simulation.

The number of pipe flow reversals in link j over the set of hydraulic time steps is defined as
(1)
where is the flow in link j computed at time step k. Pipe flow reversals can deteriorate water quality directly by increasing the water age in parts of the network, especially in looped configurations.
The flow velocity CV in link j over the set of hydraulic time steps K is defined as
(2)
where and are the standard deviation and mean pipe flow velocity in link j, respectively. The CV metric quantifies the variability in pipe flow velocity relative to the mean across a 24-h demand period. Similar to flow reversals, flow velocity variability directly influences water age and thereby disinfectant decay.
The mean source trace at node i over the set of hydraulic time steps K is defined as
(3)
where and represent the percentage of water supplied from DMA inlets 1 and 2, respectively, at node i and time step k. We take the maximum mean source trace over the simulation period K as the governing measure of source water mixing. The metric highlights the extent of mixing in the Field Lab, which is directly related to control at DMA inlets, flow across DBVs, and changes to the hydraulic gradient caused by relatively large flushing demands at AFVs.

Disinfectant variability

We define water quality variability as the standard deviation of disinfectant residuals during the last 24-h period of the 7-day water quality simulation:
(4)
where is the simulated disinfectant concentration at node i and time step k; and is the mean concentration at node i over the set of time steps , in which corresponds to the last 96 discrete 15-min time steps of the 7-day simulation. The first 6 days of the simulation are discarded to ensure water quality reaches periodic behaviour across all nodes in the network.

Sensor data

We analyse disinfectant time series from the deployed MetriNet sensors to validate the correlation between hydraulic variability and water quality dynamics in the Field Lab. The sensor data correspond to the adaptive AZP control configuration currently operated in the Field Lab. Our analysis first removes periods with missing data from the sensor time series. A data differencing step is then applied to accentuate short-term fluctuations and eliminate any trends from the inlet concentrations. The first-order differencing
(5)
creates a new time series for all sensors using the previous 15-min data point . Note that we exclude fluctuations in resulting from the removal of missing data entries. The 15-min interval between measurements effectively captures dynamics in the time series data; however, this parameter might vary with different network conditions. Future work should explore additional methods to highlight patterns of short-term variability. For instance, analysing the fractal dimension of the time series could reveal self-similarity across time scales corresponding to hydraulic control in the network (Higuchi 1988).

Simulation

Flow variability

Figure 5 plots the flow variability metrics described in Section 3.2, where results for the single-feed DMA, adaptive AZP, and adaptive AZP-SCC control configurations are presented from left to right. Subplots represent the number of flow reversals (Figure 5(a)–5(c)) and flow velocity coefficient of variation (Figure 5(d)–5(f)) across all links , and mean source trace at all nodes (Figure 5(g)–5(i)).
Figure 5

Flow variability metrics over a 24-h simulation. Subplots (a–c): R; (d–f): ; (g–i) . Left to right: single-feed DMA, adaptive AZP control, adaptive AZP-SCC control.

Figure 5

Flow variability metrics over a 24-h simulation. Subplots (a–c): R; (d–f): ; (g–i) . Left to right: single-feed DMA, adaptive AZP control, adaptive AZP-SCC control.

Close modal

The single-feed DMA configuration yielded minimal hydraulic dynamics due to the isolation of DMAs 1 and 2, characterized by the closed IVs, a single inlet connection, and predominantly branched topology. Apart from a few looped sections of the network (e.g. near Sensor 7), flow paths in the Field Lab remained unchanged throughout the simulation period. In the adaptive AZP configuration, however, we observed an increase in flow reversals (Figure 5(b)) and flow velocity CV (Figure 5(e)) at links in the vicinity of Sensors 2 and 5. This is a direct consequence of dynamic DBV operations switching between partially opened (daytime) and closed (nighttime) states. Moreover, the dynamically adaptive network connectivity resulted in a shift of the hydraulic balance point from the original DMA boundary to an area just upstream in DMA 2. Consequently, this part of the network conveyed a much larger range in flow and experienced several flow reversals. Such dynamics are corroborated by the mean source trace results in Figure 5(h).

Flow variability was exacerbated in the adaptive AZP-SCC configuration. Most notably, activated PRV and AFV settings during the 1-h self-cleaning period shifted the hydraulic balance point even further towards the inlet of DMA 2. This shift is best depicted by the mean source trace results in Figure 5(i). Additionally, increased flow reversals and pipe flow velocity CV in Figure 5(c) and 5(f), respectively, were observed in looped or dead-end links connecting AFV actuators, highlighting the elevated flow velocities generated by flushing demands. The increase in changing hydraulic conditions from the self-cleaning mode emphasize the potential impact of advanced hydraulic control on water quality dynamics.

Water quality dynamics

In this section, we present the water quality simulation results and investigate the correlation between disinfectant variability and the proposed flow variability metrics. Figure 6 illustrates the standard deviation of disinfectant residuals across the network. Subplots depict the mean (Figure 6(a)–6(c)) and worst-case (Figure 6(d)–6(f)) inlet concentration scenarios detailed in Table 1. Lastly, results for the single-feed DMA, adaptive AZP, and adaptive AZP-SCC control configurations are presented from left to right.
Figure 6

Standard deviation (SD) of disinfectant residuals over the last 24-h period of a 7-day simulation. Subplots (a–c): mean source conditions; (d–f): worst-case source conditions. Left to right: single-feed DMA, adaptive AZP control, adaptive AZP-SCC control.

Figure 6

Standard deviation (SD) of disinfectant residuals over the last 24-h period of a 7-day simulation. Subplots (a–c): mean source conditions; (d–f): worst-case source conditions. Left to right: single-feed DMA, adaptive AZP control, adaptive AZP-SCC control.

Close modal

The single-feed DMA configuration exhibited minimal disinfectant variability during the 24-h simulation, which was expected due to the relatively static hydraulic conditions. In contrast, the open network configurations corresponding to the adaptive AZP and AZP-SCC control objectives presented moderately higher disinfectant variability. Specifically, Figure 6(b) and 6(c) displays areas with higher standard deviation within the areas of increased hydraulic variability depicted in Figure 5. On the basis of these results, the proposed flow variability metrics can serve as a good surrogate measure for identifying increased disinfectant variability in the absence of a calibrated water quality model.

The observed disinfectant variability is mainly caused by the mixing of source water with different inlet concentrations. A larger disparity in inlet concentrations exacerbates this variability, as illustrated by the worst-case scenario in Figure 6(e) (AZP control) and Figure 6(f) (AZP-SCC control). Furthermore, the worst-case scenario exhibits a clear increase in the spatial extent of nodes with higher disinfectant variability. These simulations demonstrate the importance of controlling inlet concentrations for managing water quality in adaptive DMAs.

Time series analysis

We analysed differenced time series for all sensors deployed in the Field Lab to validate our simulation results for the adaptive AZP control configuration. Figure 7 presents a frequency heatmap of differenced data using bin sizes of 0.1 mg L−1, where each cell indicates the percentage of data points (rounded to the nearest 0.1%) within the specified interval. It should be noted that the average discrepancy between DPD measurements at all sensor locations recorded over the monthly calibration visits was on the order of 0.05 mg L−1, and therefore within the differenced data bin sizes presented in Figure 7.
Figure 7

Frequency heatmap of differenced disinfectant data for AZP control configuration.

Figure 7

Frequency heatmap of differenced disinfectant data for AZP control configuration.

Close modal

In Figure 7, the relatively high variability observed at Sensors 2 and 5 correlates with increased hydraulic and water quality dynamics from our model simulations, supporting our hypothesis that increases in flow variability from networks operated with dynamic connectivity can lead to greater disinfectant variability. However, as discussed in Section 4.1.2, we emphasize that this correlation depends on the disparity in inlet concentrations and the decay of disinfectant materializing between the inlet and node where mixing occurs; in this case, near the DBVs. Consequently, the relatively modest disinfectant variability observed at Sensors 2 and 5 is likely due to comparable inlet concentrations recorded at Sensors 1 (DMA 2) and 4 (DMA 1). Furthermore, changes in flow velocity for the AZP configuration were relatively minor compared to the AZP-SCC controls, which have a 1-h flushing period. This suggests that additional disinfectant decay materializing from hydraulic dynamics is limited in the Field Lab. Subsequent experiments are planned to observe disinfectant decay behaviour for the more dynamic AZP-SCC control configuration. At Sensor 7, our simulations revealed several pipes with flow reversals and higher flow velocity variability in this region (see Figure 5(b) and 5(e)). The hydraulic balance point in this section of looped pipes appears to be continuously shifting due to changing demands. As a result, flow comprising parcels of water with much higher age is likely contributing to the observed disinfectant variability at Sensor 7.

At Sensors 2 and 5, we observed a pattern of sharp drops in disinfectant residuals once or twice over each 24-h period. Figure 8 illustrates this variability at Sensor 2 over 4 days. The first and more consistent drop (shaded in grey) occurs when DBVs transition from their closed (nighttime) to open (daytime) positions, causing mixing of flows with different disinfectant concentrations across the DMA boundaries. Such variability is partly due to lower disinfectant levels in dead-end sections created when DBVs are closed. However, as discussed previously, it is also largely dictated by the disparity in DMA inlet concentrations. The second drop appears in all sensors supplied by DMA 2, indicating that it is likely caused by upstream storage and/or mixing dynamics. The impacts of such upstream dynamics are an equally important consideration for managing water quality dynamics. In any case, we highlight the repeatability of this temporal variability at Sensors 2 and 5 by comparing boxplots of differenced time series data across a 24-h period at all sensor locations (see Figure S1 in Supplementary Material).
Figure 8

Example of temporal disinfectant variability at Sensor 2 (installed at a DBV site). Shaded grey regions indicate a sharp drop in disinfectant coinciding with DBV operations.

Figure 8

Example of temporal disinfectant variability at Sensor 2 (installed at a DBV site). Shaded grey regions indicate a sharp drop in disinfectant coinciding with DBV operations.

Close modal

The impacts of dynamically adaptive connectivity and hydraulic control on water quality dynamics were investigated in a real-world water distribution network. We simulated the variability of disinfectant residuals resulting from two adaptive control configurations: (i) minimization of AZP and (ii) integration of AZP with a SCC mode. Model simulations were supported by continuous disinfectant data from a large-scale water quality experimental programme. The main conclusions of this study are as follows:

  • Flow variability metrics. The proposed metrics effectively identified areas of the network with greater short-term disinfectant variability due to hydraulic control. These metrics can enhance water quality observability by informing additional monitoring locations and can be used to model water quality dynamics in optimization problems formulated for control.

  • Validation of model simulations. Continuous disinfectant data collected from April to July 2023 validated our simulation results for the adaptive AZP control configuration, revealing greater disinfectant variability near the hydraulic balance point of the aggregated (opened) DMA configuration. This variability was more pronounced during the period when DBVs transitioned from their closed to open position, which caused mixing between DMAs.

  • Influence of inlet concentrations. The extent of disinfectant variability caused by DMAs with dynamic connectivity is influenced by the disparity in disinfectant concentrations at the inlets. Controlling inlet concentrations is therefore critical for managing water quality in dynamically adaptive networks.

These conclusions support further analytical and experimental research needed to enhance water quality management in real-world systems. Future work will focus on the controllability of disinfectant residuals in adaptive networks by incorporating the proposed flow variability metrics into the formulation of optimal design and control problems. These control problems should also explicitly model water quality dynamics as constraints and consider the impacts of uncertainties in disinfectant decay parameters. Concurrently, ongoing experimental research is being conducted to observe water quality conditions over a broader range of hydraulic conditions, including the self-cleaning control configuration presented in this study. This experimental work also aims to further characterize water quality dynamics using field data collected across different seasonal, short-term, and long-term conditions, which is critical for improving the accuracy of water quality models used in control applications.

Bradley Jenks is supported by Imperial College London, Bristol Water Plc, Analytical Technology, and the Natural Sciences and Engineering Research Council of Canada (PGSD-577767-2023). Angeliki Aisopou is a Daphne Jackson fellow, supported by the Royal Society and the Engineering and Physical Sciences Research Council (EPSRC). Ivan Stoianov is an Anglian Water Services and CLA-VAL UK/Royal Academy of Engineering Senior Research Fellow in Dynamically Adaptive Water Supply Networks (RCSRF2324-17-41).

All relevant data are available at: https://github.com/bradleywjenks/wq_dynamics.

The authors declare there is no conflict.

Abokifa
A. A.
,
Yang
J. Y.
,
Lo
C. S.
&
Biswas
P.
(
2016
)
Water quality modeling in the dead end sections of drinking water distribution networks
,
Water Research
,
89
,
107
117
.
doi:10.1016/j.watres.2015.11.025
.
Abokifa
A. A.
,
Maheshwari
A.
,
Gudi
R. D.
&
Biswas
P.
(
2019
)
Influence of dead-end sections of drinking water distribution networks on optimization of booster chlorination systems
,
Journal of Water Resources Planning and Management
,
145
(
12
),
04019053
.
doi:10.1061/(ASCE)WR.1943-5452.0001125
.
Ahn
J. C.
,
Lee
S. W.
,
Choi
K. Y.
,
Koo
J. Y.
&
Jang
H. J.
(
2011
)
Application of unidirectional flushing in water distribution pipes
,
Journal of Water Supply: Research and Technology – Aqua
,
60
(
1
),
40
50
.
doi:10.2166/aqua.2011.013
.
Analytical Technology Inc.
. (
2022
)
MetriNet Q52 Literature
. .
Armand
H.
,
Stoianov
I.
&
Graham
N. J. D.
(
2017
)
A holistic assessment of discolouration processes in water distribution networks
,
Urban Water Journal
,
14
(
3
),
263
277
.
doi:10.1080/1573062X.2015.1111912
.
Armand
H.
,
Stoianov
I.
&
Graham
N.
(
2018
)
Impact of network sectorisation on water quality management
,
Journal of Hydroinformatics
,
20
(
2
),
424
439
.
doi:10.2166/hydro.2017.072
.
Blokker
E. J. M.
,
Furnass
W. R.
,
Machell
J.
,
Mounce
S. R.
,
Schaap P
G.
&
Boxall J
B.
(
2016
)
Relating water quality and age in drinking water distribution systems using self-organising maps
,
Environments
,
3
(
2
),
10
.
doi:10.3390/environments3020010
.
Boxall
J.
,
Blokker
M.
,
Schaap
P.
,
Speight
V.
&
Husband
S.
(
2023
)
Managing discolouration in drinking water distribution systems by integrating understanding of material behaviour
,
Water Research
,
243
,
120416
.
doi:10.1016/j.watres.2023.120416
.
Bui
X. K.
,
Jeong
G.
&
Kang
D.
(
2022
)
Adaptive DMA design and operation under multiscenarios in water distribution networks
,
Sustainability
,
14
(
6
),
3692
.
doi:10.3390/su14063692
.
Derlon
N.
,
Massé
A.
,
Escudié
R.
,
Bernet
N.
&
Paula
E.
(
2008
)
Stratification in the cohesion of biofilms grown under various environmental conditions
,
Water Research
,
42
,
2102
2110
.
doi:10.1016/j.watres.2007.11.016
.
Dias
V. C. F.
,
Besner
M.-C.
&
Prévost
M.
(
2017
)
Predicting water quality impact after district metered area implementation in a full-scale drinking water distribution system
,
Journal AWWA
,
109
(
9
),
363
380
.
doi:10.5942/jawwa.2017.109.0099
.
Furnass
W. R.
,
Mounce
S. R.
&
Boxall
J. B.
(
2013
)
Linking distribution system water quality issues to possible causes via hydraulic pathways
,
Environmental Modelling & Software
,
40
,
78
87
.
doi:10.1016/j.envsoft.2012.07.012
.
Giudicianni
C.
,
Herrera
M.
,
di Nardo
A.
&
Adeyeye
K.
(
2020
)
Automatic multiscale approach for water networks partitioning into dynamic district metered areas
,
Water Resources Management
,
34
(
2
),
835
848
.
doi:10.1007/s11269-019-02471-w
.
Higuchi
T.
(
1988
)
Approach to an irregular time series on the basis of the fractal theory
,
Physica D: Nonlinear Phenomena
,
31
(
2
),
277
288
.
doi:10.1016/0167-2789(88)90081-4
.
Horn
H.
,
Reiff
H.
&
Morgenroth
E.
(
2003
)
Simulation of growth and detachment in biofilm systems under defined hydrodynamic conditions
,
Biotechnology and Bioengineering
,
81
(
5
),
607
617
.
doi:10.1002/bit.10503
.
Jenks
B.
,
Pecci
F.
&
Stoianov
I.
(
2023a
)
Optimal design-for-control of self-cleaning water distribution networks using a convex multi-start algorithm
,
Water Research
,
231
,
119602
.
doi:10.1016/j.watres.2023.119602
.
Jenks
B.
,
Ulusoy
A.-J.
,
Pecci
F.
&
Stoianov
I.
(
2023b
)
Dynamically adaptive networks for integrating optimal pressure management and self-cleaning controls
,
Annual Reviews in Control
,
55
,
489
497
.
doi:10.1016/j.arcontrol.2023.03.014
.
Jenks
B.
,
Aisopou
A.
,
Ulusoy
A.-J.
&
Stoianov
I.
(
2024
)
Investigating Water Quality Dynamics in Distribution Networks with Dynamic Topology and Hydraulic Control
. Available at: https://github.com/bradleywjenks/wq_dynamics
(Accessed 15 February 2024)
.
Klise
K. A.
,
Bynum
M.
,
Moriarty
D.
&
Murray
R.
(
2017
)
A software framework for assessing the resilience of drinking water systems to disasters with an example earthquake case study
,
Environmental Modelling & Software
,
95
,
420
431
.
doi:10.1016/j.envsoft.2017.06.022
.
Lee
Y.
,
Sinha
E.
,
Piao
H.
,
Stillman
J.
,
Hartman
D.
&
Bush
C.
(
2010
) ‘
Chlorine wall decay coefficients to calibrate the GCWW all-pipes distribution system model
’,
Proceedings of World Environmental and Water Resources Congress 2010: Challenges of Change
,
16–20 May
.
Providence, USA
, pp.
4405
4413
doi:10.1061/41114(371)448
.
Liu
M.
,
Craik
S.
&
Zhu
D.
(
2015
)
Determination of cast iron pipe wall decay coefficient for combined chlorine in a municipal water distribution system
,
Canadian Journal of Civil Engineering
,
42
(
4
),
250
258
.
doi:10.1139/cjce-2014-0449
.
Machell
J.
&
Boxall
J.
(
2014
)
Modeling and field work to investigate the relationship between age and quality of tap water
,
Journal of Water Resources Planning and Management
,
140
(
9
),
04014020
.
doi:10.1061/(ASCE)WR.1943-5452.0000383
.
Mala-Jetmarova
H.
,
Barton
A.
&
Bagirov
A.
(
2015
)
Exploration of the trade-offs between water quality and pumping costs in optimal operation of regional multiquality water distribution systems
,
Journal of Water Resources Planning and Management
,
141
(
6
),
04014077
.
doi:10.1061/(ASCE)WR.1943-5452.0000472
.
Maleki
M.
,
Ardila
A.
,
Argaud P
-O.
,
Pelletiera
G.
&
Rodriguez
M.
(
2023
)
Full-scale determination of pipe wall and bulk chlorine degradation coefficients for different pipe categories
,
Water Supply
,
23
(
2
),
657
670
.
doi:10.2166/ws.2023.020
.
Montiel
F.
,
Coutelan
J.
&
Nguyen
B.
(
2002
)
A new real-time drinking water network mixing model for the control and command centre of the city of Paris
,
Water Science and Technology
,
45
(
4–5
),
383
388
.
doi:10.2166/wst.2002.0631
.
Ostfeld
A.
&
Salomons
E.
(
2004
)
Optimal operation of multiquality water distribution systems: Unsteady conditions
,
Engineering Optimization
,
36
(
3
),
337
359
.
doi:10.1080/0305215042000207054
.
Rossman
L.
,
Woo
H.
,
Tryby
M.
,
Shang
F.
,
Janke
R.
&
Haxton
T.
(
2020
)
EPANET 2.2 User Manual
.
Cincinnati, OH
:
US Environmental Protection Agency
.
Ulusoy
A.-J.
,
Nerantzis
D.
&
Stoianov
I.
(
2023
)
Adaptive MPC for burst incident management in water distribution networks
,
IEEE Transactions on Control of Network Systems
,
10
(
4
),
2060
2071
.
doi:10.1109/TCNS.2023.3259103
.
US Environmental Protection Agency
. (
2002
)
Effects of Water Age on Distribution System Water Quality
.
EPA
,
Washington, DC, USA
.
US Environmental Protection Agency
. (
2008
)
Water Security Initiative: Cincinnati Pilot Post-Implementation System Status
.
EPA
,
Cincinnati, OH, USA
.
Vasconcelos
J.
,
Rossman
L.
,
Grayman
W.
,
Boulos
P.
&
Clark
R.
(
1997
)
Kinetics of chlorine decay
,
Journal AWWA
,
89
(
7
),
54
65
.
doi:10.1002/j.1551-8833.1997.tb08259.x
.
Wright
R.
,
Stoianov
I.
,
Parpas
P.
,
Henderson
K.
&
King
J.
(
2014
)
Adaptive water distribution networks with dynamically reconfigurable topology
,
Journal of Hydroinformatics
,
16
(
6
),
1280
1301
.
doi:10.2166/hydro.2014.086
.
Wright
R.
,
Abraham
E.
,
Parpas
P.
&
Stoianov
I.
(
2015
)
Control of water distribution networks with dynamic DMA topology using strictly feasible sequential convex programming
,
Water Resources Research
,
51
(
12
),
9925
9941
.
doi:10.1002/2015WR017466
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).

Supplementary data