COUPLED ENGINEERED AND NATURAL DRAINAGE NETWORKS: DATA-MODEL SYNTHESIS IN URBANIZED RIVER BASINS

Abstract

In urbanized river basins, sanitary wastewater and urban runoff (non-sanitary water) from urban agglomerations drain to complex engineered networks, are treated at centralized wastewater treatment plants (WWTPs) and discharged to river networks. Discharge from multiple WWTPs distributed in urbanized river basins contributes to impairments of river water-quality and aquatic ecosystem integrity. The size and location of WWTPs are determined by spatial patterns of population in urban agglomerations within a river basin. Economic and engineering constraints determine the combination of wastewater treatment technologies used to meet required environmental regulatory standards for treated wastewater discharged to river networks. Thus, it is necessary to understand the natural-human-engineered networks as coupled systems, to characterize their interrelations, and to understand emergent spatiotemporal patterns and scaling of geochemical and ecological responses. My PhD research involved data-model synthesis, using publicly available data and application of well-established network analysis/modeling synthesis approaches. I present the scope and specific subjects of my PhD project by employing the Drivers-Pressures-Status-Impacts-Responses (DPSIR) framework. The defined research scope is organized as three main themes: (1) River network and urban drainage networks (Foundation-Pathway of Pressures); (2) River network, human population, and WWTPs (Foundation-Drivers-Pathway of Pressures); and (3) Nutrient loads and their impacts at reach- and basin-scales (Pressures-Impacts).Three inter-related research topics are: (1) the similarities and differences in scaling and topology of engineered urban drainage networks (UDNs) in two cities, and UDN evolution over decades; (2) the scaling and spatial organization of three attributes: human population (POP), population equivalents (PE; the aggregated population served by each WWTP), and the number/sizes of WWTPs using geo-referenced data for WWTPs in three large urbanized basins in Germany; and (3) the scaling of nutrient loads (P and N) discharged from ~845 WWTPs (five class-sizes) in urbanized Weser River basin in Germany, and likely water-quality impacts from point- and diffuse- nutrient sources. I investigate the UDN scaling using two power-law scaling characteristics widely employed for river networks: (1) Hack’s law (length-area power-law relationship), and (2) exceedance probability distribution of upstream contributing area. For the smallest UDNs, length-area scales linearly, but power-law scaling emerges as the UDNs grow. While area-exceedance plots for river networks are abruptly truncated, those for UDNs display exponential tempering. The tempering parameter decreases as the UDNs grow, implying that the distribution evolves in time to resemble those for river networks. However, the power-law exponent for mature UDNs tends to be larger than the range reported for river networks. Differences in generative processes and engineering design constraints contribute to observed differences in the evolution of UDNs and river networks, including subnet heterogeneity and non-random branching.In this study, I also examine the spatial patterns of POP, PE, and WWTPs from two perspectives by employing fractal river networks as structural platforms: spatial hierarchy (stream order) and patterns along longitudinal flow paths (width function). I propose three dimensionless scaling indices to quantify: (1) human settlement preferences by stream order, (2) non-sanitary flow contribution to total wastewater treated at WWTPs, and (3) degree of centralization in WWTPs locations. I select as case studies three large urbanized river basins (Weser, Elbe, and Rhine), home to about 70% of the population in Germany. Across the three river basins, the study shows scale-invariant distributions for each of the three attributes with stream order, quantified using extended Horton scaling ratios; a weak downstream clustering of POP in the three basins. Variations in PE clustering among different class-sizes of WWTPs reflect the size, number, and locations of urban agglomerations in these catchments. WWTP effluents have impacts on hydrologic attributes and water quality of receiving river bodies at the reach- and basin-scales. I analyze the adverse impacts of WWTP discharges for the Weser River basin (Germany), at two steady river discharge conditions (median flow; low-flow). This study shows that significant variability in treated wastewater discharge within and among different five class-sizes WWTPs, and variability of river discharge within the stream order