Nonpoint source pollution of urban stormwater runoff: a methodology for source analysis

Vertical distribution and migration of heavy metals in soil of green stormwater infrastructure receiving roof runoff

In order to achieve effective management of urban stormwater runoff, green stormwater infrastructure (GSI) has been widely used worldwide. However, the problem of heavy metal contamination in GSI soils has gradually become a limiting factor for their development. In this paper, concentrations of 6 heavy metals were detected in soils from 0 to 80 cm depth in the GSI receiving roof runoff. The total concentration of the 6 heavy metals in topsoil of the GSI was 395 mg·kg-1, which was significantly higher than that of the control site (p < 0.01). Vertically, heavy metal concentrations in GSI soils decreased with depth. The concentrations of As, Cr, Cu, Pb and Zn in the soils at different depths (0-80 cm) of GSI were higher than those in the soils at the corresponding depths at the control site. This suggests that heavy metals may have migrated below 80 cm. The decrease in Zn/Cu ratios with depth indicated that Zn has a higher initial concentration and more of Zn accumulated in upper layer soil than Cu. Finally, we simulated the vertical migration of dissolved heavy metals in GSI soils by the Hydrus-1D model. Heavy metal concentrations in runoff and soil particle size distributions were examined, while adsorption isotherm experiments were conducted to determine input parameters for the model. The simulation results showed that only dissolved Zn would accumulate in the top soil layer. A short time of infiltration (<60 min) may have little effect on the profile shape of dissolved heavy metal levels in the soil. Substantial increases in infiltration may result in the migration of dissolved heavy metals to deeper depths in the soil.

Review: The application of source analysis methods in tracing urban non-point source pollution: categorization, hotspots, and future prospects

The contribution of urban non-point source (NPS) pollution to surface water pollution has gradually increased, analyzing the sources of urban NPS pollution is of great significance for precisely controlling surface water pollution. A bibliometric analysis of relevant research literature from 2000 to 2021 reveals that the main methods used in the source analysis research of urban NPS pollution include the emission inventory approach, entry-exit mass balance approach, principal component analysis (PCA), positive matrix factorization (PMF) model, etc. These methods are primarily applied in three aspects: source analysis of rainfall-runoff pollution, source analysis of wet weather flow (WWF) pollution in combined sewers, and analysis of the contribution of urban NPS to the surface water pollution load. The application of source analysis methods in urban NPS pollution research has demonstrated an evolution from qualitative to quantitative, and further towards precise quantification. This progression has transitioned from predominantly relying on on-site monitoring to incorporating model simulations and employing mathematical statistical analyses for traceability. This paper reviews the principles, advantages, disadvantages, and the scope of application of these methods. It also aims to address existing problems and analyze potential future development directions, providing valuable references for subsequent related research.

Sediment recontamination potential and biological impacts of hydrophobic organics from stormwater in a mixed-use watershed

This study focused on understanding the influence of stormwater on sediment recontamination and biota bioaccumulation. Stormwater and sediment sampling was conducted at Paleta Creek San Diego, CA, a mixed-use urban watershed. Stormwater samples were size fractioned into clay, silt, sand, and filtered fractions to better understand the association of PAHs and PCBs with different particle sizes. Settling traps were deployed during storm seasons to capture solid particles deposited by stormwater. Bioaccumulation studies of the receiving water sediments were conducted on bent-nosed clams (Macoma nasuta) both in-situ and ex-situ. The results indicated that the majority of PAHs and PCBs were discharged during the early stages of the storm and were primarily associated with coarser particles (>20 μm). Sediment near the discharge exhibited higher levels of contaminants, consistent with the findings from the settling traps and stormwater discharges. Ex-situ bioaccumulation studies demonstrated that sediment contamination did not correlate with bioaccumulation in Macoma nasuta, whereas porewater accurately reflected bioaccumulation trends. This study highlights the importance of considering stormwater discharge patterns and particle sizes when assessing sediment recontamination but also that the resulting bulk sediment contamination may not reflect bioavailability as measured by organism bioaccumulation. These insights contribute to a better understanding of the impacts of stormwater runoff on sediment and biota in southern California coastal watersheds, aiding in the development of effective management strategies.

Accumulation and contamination of gully pot sediments from varied land-use types: metal loads, concentrations and speciation

Urban stormwater typically enters sewer networks through gully pots, which allow a primary sedimentation of solids upstream of the piped network. The regular removal and disposal of retained sediment are necessary, costly and can involve environmental risks due to the contamination of sediments with substances from the urban environment such as metals. The concentrations and speciation of Cd, Cr, Cu, Ni, Pb and Zn were analysed in sediments from 26 gully pots located in different land use areas in Stockholm, Sweden. In addition, accumulation rates of both sediment and metal masses were evaluated, providing a basis for optimising maintenance practices and better understanding of impacts of characteristic urban land use types. Metal concentrations varied by at most a factor of eight between samples and were always below Swedish polluted site guidelines for less sensitive land use, with only eight samples exceeding the guideline values for Cu and Zn for sensitive land use. Sequential extraction showed Pb and Zn to be the most mobile metals. Sediment accumulation rates varied from 0.003 to 0.197 kg/m2 impermeable surface/year. Metal accumulation rates were much more variable than metal concentrations, with a factor of up to 172 between the highest and lowest rates and the highest metal accumulation rates corresponding to the lower range of mass loads in road runoff. Differences in metal concentrations, sediment or metal mass accumulations could not be solely attributed to either traffic or catchment land use. In contrast, traction grit used for winter road maintenance, which has low (but detectable) metal concentrations, is identified as a major component of gully pot sediments, with a combined effect of both moderating metal concentrations and contributing to total mass.

Molecularly imprinted polymers (MIPs) as efficient catalytic tools for the oxidative degradation of 4-nonylphenol and its by-products

Water pollution is a current global concern caused by emerging pollutants like nonylphenol (NP). This endocrine disruptor cannot be efficiently removed with traditional wastewater treatment plants (WTPs). Therefore, this work aimed to evaluate the adsorption influence of molecularly imprinted polymers (MIPs) on the oxidative degradation (ozone and ultraviolet irradiations) of 4-nonylphenol (4-NP) and its by-products as a coadjuvant in WTPs. MIPs were synthesized and characterized; the effect of the degradation rate under system operating conditions was studied by Box-Behnken response surface design of experiments. The variables evaluated were 4-NP concentration, ozone exposure time, pH, and MIP amount. Results show that the MIPs synthesized by co-precipitation and bulk polymerizations obtained the highest retention rates (> 90%). The maximum adsorption capacities for 4-NP were 201.1 mg L-1 and 500 mg L-1, respectively. The degradation percentages under O3 and UV conditions reached 98-100% at 120 s of exposure at different pHs. The degradation products of 4-NP were compounds with carboxylic and ketonic acids, and the MIP adsorption was between 50 and 60%. Our results present the first application of MIPs in oxidation processes for 4-NP, representing starting points for the use of highly selective materials to identify and remove emerging pollutants and their degradation by-products in environmental matrices.

The fate of three typical persistent organic pollutants in bioretention columns as revealed by stable carbon isotopes

There is a lack of simple and effective methods to quantify the fate processes of persistent organic pollutants (POPs) in bioretention systems. In this study, the fate and elimination processes of three typical ¹³C-labeled POPs in regularly added bioretention columns were quantified using stable carbon isotope analysis techniques. The results showed that the modified media bioretention column removed more than 90% of Pyrene, PCB169 and p,p'-DDT. Media adsorption was the dominant removal mechanism for the reduction of the three exogenous organic compounds (59.1-71.8% of the input) although plant uptake (5.9-18.0%) was also important. Mineralization was effective in degrading pyrene (13.1%) but had a very limited effect on p,p'-DDT and PCB169 removal (<2.0%), the reason for which may be related to the aerobic conditions of the filter column. Volatilization was relatively weak and negligible (<1.5%). The presence of heavy metals inhibited the removal of POPs to some extent: media adsorption, mineralization and plant uptake were reduced by 4.3-6.4%, 1.8-8.3% and 1.5-3.6% respectively. This study suggests that bioretention systems are an effective measure for the sustainable removal of POPs from stormwater and that heavy metals can inhibit the overall performance of the system. Stable carbon isotope analysis techniques can help to investigate the migration and transformation of POPs in bioretention systems.

Spatiotemporal trends and impact of Covid-19 lockdown on eight sewage contaminants in Brisbane, Australia, from 2012 to 2020

This study aims to investigate the spatiotemporal trends and impact of COVID-19 lockdowns to the profile of physiochemical parameters in the influent of wastewater treatment plants (WWTPs) around Brisbane, Australia. One 24-hr composite influent sample was collected from 10 WWTPs and analyzed for a range of physiochemical parameters per week (i.e., chemical oxygen demand (COD), total organic carbon (TOC), total nitrogen (TN), total phosphorus (TP), ammonia, volatile suspended solid (VSS)) and per month (i.e., Ni and Cr) from 2012 to 2020, including the period of COVID-19 lockdowns in the region. The catchments studied were urban, with a mix of domestic and industrial activities contributing towards the contaminant profile. Statistical analysis identified that industrial and commercial land use, as well as population size had a large impact to the parameter loads and profile. Per capita mass loads of Cr in one catchment were 100 times higher than in others from one industrial point source. TP demonstrated a potential monotonic decrease over time due to practical reduction policies that have been implemented for phosphorous content in household detergents, except for one catchment where trade waste from food manufacturing industries contributed to an overall increase of 6.9%/year TP. The COVID-19 lockdown (March-April 2020) posed different impact on different catchments, either decrease (7-61%) or increase (2-40%) of most parameter loads (e.g., COD, TOC, TN, TP, VSS, Ammonia), which was likely driven by catchment characteristics (i.e., the proportion of residential, commercial, and industrial land uses). This study enhances our understanding of spatiotemporal trend of contaminants in the catchments for further effective source control.

Spatiotemporal occurrence of phthalate esters in stormwater drains of Hong Kong, China: Mass loading and source identification

Urban stormwater is an important pathway for transporting anthropogenic pollutants to water bodies. Phthalate esters (PAEs) are endocrine disruptors owing to their estrogenic activity and potential carcinogenicity and their ubiquitous presence has garnered global interest. However, their transportation by urban stormwater has been largely overlooked. This study, for the first time, investigated 15 PAEs in stormwater from six major stormwater drains in the highly urbanized Hong Kong, a major metropolitan city in China. The results showed that PAEs were ubiquitous in the stormwater of Hong Kong, with total concentrations (∑₁₅PAEs) spanning from 195 to 80,500 ng/L. Bis(2-n-butoxyethyl) phthalate (DBEP), diisopentyl phthalate (DiPP), dicyclohexyl phthalate (DCHP) and di-n-pentyl phthalate (DnPP) were detected in stormwater for the first time. Spatial variations in PAEs were observed among different stormwater drains, possibly due to the different land use patterns and intensities of human activities in their respective catchments. The highest and lowest levels of ∑₁₅PAEs were found in Kwai Chung (3860 ± 1960 ng/L) and the Ng Tung River (672 ± 557 ng/L), respectively. Additionally, significantly higher concentrations of ∑₁₅PAEs in stormwater were found in the wet season (2520 ± 2050 ng/L) than in the dry season (947 ± 904 ng/L). Principal component analysis classified domestic and industrial origins as two important sources of PAEs in the stormwater of Hong Kong. Stormwater played a crucial role in transporting PAEs, with an estimated annual flux of 0.705-29.4 kg. Thus, possible stormwater management measures were proposed to protect the receiving environment and local ecosystems from stormwater.

Coupling with high-resolution remote sensing data to evaluate urban non-point source pollution in Tongzhou, China

Urban non-point source (NPS) pollution has gradually become one of the important factors affecting the urban water environment. The quantitative evaluation of urban NPS pollution is the priority to identify key control area of urban NPS pollution. Current model applied in China is mainly focused on small-scale area, large-scale spatial continuous simulation is lacking. In this study A spatial continuous evaluation method coupled with high-resolution remote sensing data has been established and the method was applied to Tongzhou, China. With the spatial distribution of land-use type and built-up area which were been obtained by remote sensing technology, the accumulative and wash-off load of urban NPS nitrogen and phosphorus were estimated for the prominent problems of nitrogen and phosphorus nutrient pollution in the rivers in the study area. The main sources of urban NPS Nitrogen and phosphorus pollution are roof and road rainfall runoff respectively. Compared to other urban NPS pollution models, the method developed in this study can quickly realize spatial visualization assessment of urban NPS pollution and provide a means to estimate urban NPS loads in entire city or urban agglomeration, it is applicable for common urban NPS pollutants and also has advantages in areas without data.

Addressing the social barriers to green stormwater infrastructure in residential areas from a socio-ecological perspective

Despite the well-recognized financial limitations, social aspects can also impact the adoption of green stormwater infrastructure (GSI) as it is ingrained in the socio-ecological system we live in. Thus, this work focuses on gaining an understanding of the public's perceptions of GSI by considering cognitive biases that hinder its adoption. This work is composed of two forms of human-subject studies, including an online-based survey and a series of semi-structured interviews. The survey (n = 510) was conducted to gauge public opinions toward GSI, whereas the interviews with representatives of major local regulatory agencies were to learn about the logistics for GSI implementation in Mecklenburg County, NC. The results were interpreted using the theory of planned behavior of rational actors. Statistical results showed a weak interpretation through this theory to explain the survey participants' intention to adopt GSI measures. This could suggest that the incorporation of irrationality, such as cognitive biases, could further enhance the predictability of the theory. At the same time, an inconsistency between the findings from the survey and the interviews was identified: most survey participants showed an overall uniform positive attitude, intention, and behavior regarding GSI practice adoption, whereas the interviewed experts all suggested a wide diversity on such terms. Suggestions were made based on the findings for better policy-making on public engagement for local regulatory agencies. This research aims to help local stormwater management authorities explore shortcomings in current stakeholder engagement plans to gain sustainable support for GSI implementation in urbanized areas.

Analysis of micropollutants in a marine outfall using network analysis and decision tree

The presence of micropollutants (MPs), including pharmaceutical, industrial, and pesticidal compounds, threatens both human health and the aquatic ecosystem. The development and extensive use of new chemicals have also inevitably led to the accumulation of MPs in aquatic environments. Recreational beaches are especially vulnerable to contamination, affecting humans and aquatic animals via the absorption of MPs in water during marine activities (e.g., swimming, sailing, and windsurfing). Additionally, marine outfalls in an urbanized coastal city can cause serious chemical and microbial pollution on recreational beaches, leading to an increase in adverse effects on public health and the ecological system. Therefore, the aim of this study was to, with the use of network and decision tree analyses, identify the features and factors that influence the change in MP concentrations in a marine outfall. These analyses were conducted to inspect the relationship between each MP and its hierarchical structure as well as hydrometeorological variables. Additionally, a risk analysis was conducted in this study in which the MPs were prioritized based on their optimized risk quotient values. During our monitoring of MP concentrations over time at the marine outfall, high concentrations of pharmaceutical and industrial compounds were detected when the tide level was low after rainfall. Furthermore, results of the risk analysis and the prioritization revealed that a total of 18 substances identified in our study posed a risk to the ecosystem; these include major ecotoxicologically hazardous substances such as telmisartan, mevinphos, and methiocarb. Results of the network analysis demonstrated distinct trends for pharmaceutical and industrial substances, whilst those for pesticide compounds were irregular. Additionally, the hierarchical structures for most MPs consisted of rainfall, tide level, and antecedent dry hours; this implies that these factors influence MP dynamics. These findings will be helpful for establishing chemical contamination management plans for recreational beaches in the future.

Lead Inhibits Postecdysial Exoskeletal Calcification in the Blue Crab (Callinectes sapidus)

Postecdysial mineralization in crustaceans involves the deposition of carbonate salts, such as calcium carbonate, to the organic matrix. Because of the resemblance between Pb²⁺ and Ca²⁺ , the present study was carried out to investigate whether Pb is incorporated into the new shell during postecdysial mineralization using the blue crab (Callinectes sapidus) as the model crustacean. It was hypothesized that injected Pb would be deposited in the shell via calcium transporters in the epidermis during the mineralization process. Postecdysial blue crabs were injected with two doses of 5 µg Pb/g wet weight each in lead acetate, and then Pb, Ca, and Mg contents were analyzed in the exoskeleton, while only Pb bioaccumulation was quantified for the hepatopancreas, gills, muscles, and hemolymph. The results showed a statistically nonsignificant increase in exoskeletal Pb content in Pb-treated crabs compared to control, suggesting that exoskeletal Pb may not be a sensitive proxy for aquatic Pb pollution. There was a significant decrease in Ca content in Pb-treated crabs, suggesting that Pb hindered the deposition of Ca to crab exoskeleton, thereby obstructing calcification. A trend of a decrease in exoskeletal Mg was also observed in Pb-treated crabs. There was a significant increase in Pb content found in the gills, hepatopancreas, muscle, and hemolymph in Pb-treated crabs. The rank of the Pb level among three soft tissues in a decreasing order is hepatopancreas > gill > muscle. This is the first study to present evidence that Pb disrupts postecdysial exoskeletal calcification in a crustacean. Environ Toxicol Chem 2022;41:474-482. © 2021 SETAC.

Micropollutants in Urban Runoff from Traffic Areas: Target and Non-Target Screening on Four Contrasted Sites

Although runoff from trafficked urban areas is recognized as a potentially significant pathway of micropollutants, runoff pollution remains poorly documented, except for relatively few historical pollutants such as some metals and hydrocarbons. Therefore, in this work, road and parking lot runoff from four sites with contrasting traffic levels were analyzed for a very broad spectrum of molecules and elements. A total of 128 pollutants and micropollutants were monitored, including inorganic (n = 41) and organic (n = 87) pollutants. Both the dissolved and particulate phases were considered. For a reduced number of samples, non-targeted screening by high-resolution mass spectrometry (HRMS) was carried out. For targeted screening, the contamination profiles were quite homogeneous, but the concentrations significantly differed between the different sites. Sites with the highest traffic density exhibited the highest concentrations for polycyclic aromatic hydrocarbons (PAHs), some traffic-related metals, alkylphenols and phthalates. Overall, for most micropollutants, the parking lot runoff exhibited the lowest concentrations, and the specificity of this site was confirmed by its HRMS fingerprint. Non-target screening allowed the sites to be discriminated based on the occurrence of specific compounds. Unlike the results of targeted screening, the HRMS intra-site variability was lower than its inter-site variability. Unknown substances were tentatively identified, either characteristic of each site or ubiquitous of all samples.

Biological nitrogen removal from stormwater in bioretention cells: a critical review

Excess nitrogen in stormwater degrades surface water quality via eutrophication and related processes. Bioretention has been recognized as a highly effective low-impact development (LID) technology for the management of high runoff volumes and reduction of nitrogen (N) pollutants through various mechanisms. This paper provides a comprehensive and critical review of recent developments on the biological N removal processes occurring in bioretention systems. The key plant- and microbe-mediated N transformation processes include assimilation (N uptake by plants and microbes), nitrification, denitrification, and anammox (anaerobic ammonia oxidation), but denitrification is the major pathway of permanent N removal. Overall, both laboratory- and field-scale bioretention systems have demonstrated promising N removal performance (TN: >70%). The phyla Bacteroidetes and Proteobacteria are the most abundant microbial communities found to be enriched in biofilter media. Furthermore, the denitrifying communities contain several functional genes (e.g., nirK/nirS, and nosZ), and their concentrations increase near the surface of media depth. The N removal effectiveness of bioretention systems is largely impacted by the hydraulics and environmental factors. When a bioretention system operates at: low hydraulic/N loading rate, containing a saturation zone, vegetated with native plants, having deeper and multilayer biofilter media with warm climate temperature and wet storm events periods, the N removal efficiency can be high. This review highlights shortcomings and current knowledge gaps in the area of total nitrogen removal using bioretention systems, as well as identifies future research directions on this topic.

A study on water quality monitoring of Buddha Nullah, Ludhiana, Punjab (India)

Diverse genotoxic agents, entering the aquatic environment through natural and anthropogenic events, pose serious threats to its biotic components. The present study involves the monitoring of water quality by assessing the genotoxic effects and physico-chemical parameters including heavy metals of 10 surface water samples collected from different locations of Buddha Nullah, a tributary of Sutlej flowing through Ludhiana, Punjab (India). Genotoxicity was evaluated following Allium cepa root chromosomal aberration assay and DNA nicking assay using plasmid (pBR322) whilst the metal (cadmium, chromium, cobalt, copper, lead, nickel and zinc) analysis was conducted using atomic absorption spectrophotometer. All water samples collected from the study area had cobalt and lead content more than the permissible limits (0.04 and 0.01, respectively) recommended by the Bureau of Indian Standards and the World Health Organization. The samples also induced genotoxicity following both bioassays. The water samples collected from Gaunspur (GP), a site approx. 75.53 km upstream of the Sutlej-Buddha Nullah joining point, has shown the maximum genotoxic effect, i.e. 38.62% in terms of per cent total aberrant cells during A. cepa assay and 100% DNA damage during DNA nicking assay. The Pearson correlation indicated that genotoxicity had a significant positive correlation with the content of cobalt (at p ≤ 0.5). During cluster analysis, the samples from 10 sites formed four statistically significant clusters based on the level of pollution that was dependent on two factors like similarity in physico-chemical characteristics and source of pollution at a specific site.

Contamination of Urban Stormwater Pond Sediments: A Study of 259 Legacy and Contemporary Organic Substances

Stormwater ponds improve water quality by facilitating the sedimentation of particles and particulate contaminants from urban runoff. Over time, this function entails the accumulation of contaminated sediments, which must be removed periodically to maintain a pond's hydraulic and treatment capacity. In this study, sediments from 17 stormwater sedimentation facilities from four Swedish municipalities were analyzed for 259 organic substances likely to be found in the urban environment. A total of 92 substances were detected in at least one sample, while as many as 52 substances were detected in a single sample. A typical profile of urban contamination was identified, including polychlorinated biphenyls, polycyclic aromatic hydrocarbons, organotins, aliphatic hydrocarbons, phthalates, aldehydes, polybrominated diphenyl ethers, perfluorinated substances, and alkylphenols. However, levels of contamination varied greatly between ponds, influenced heavily by the dilution of urban pollutants and wear particles from other sources of particles such as eroded soil, sand, or natural organic matter. For 22 of 32 samples, the observed concentrations of at least one organic substance exceeded the regulatory threshold values derived from toxicity data for both sediment and soil.

Nature-Based Solutions and Real-Time Control: Challenges and Opportunities

Nature-based solutions (NBS) as green infrastructures to urban drainage are an effective mitigation strategy both in terms of quantity and quality of runoff. Real-time control (RTC) can complement both flood mitigation and improvement of water quality by controlling elements of the drainage and sewage system. This study assessed the improvement opportunities with RTC of three NBS-related techniques commonly applied in urban drainage with different spatial scales: green roof, bioretention and detention basin and the remaining challenges to integrate both methods. Additionally, our investigations showed that the main difficulties reported involve the planning and monitoring stages of the RTC system. All of the studied devices can benefit from RTC. It is possible to observe that, despite the good results reported in the literature, the application of RTC to NBS studies on urban drainage are very recent. There are several opportunities that can be explored to optimize the performance.

Performance of two contrasting pilot swale designs for treating zinc, polycyclic aromatic hydrocarbons and glyphosate from stormwater runoff

Swales are a widespread stormwater management solution to reduce pollutant concentrations in runoff. An innovative pilot facility was constructed to evaluate the treatment efficiency of the two main types of water-quality swales, i.e. standard swales and filtering swales. Using stormwater roof runoff, without any additions or spiked with organic micropollutants, 12 runoff simulation runs mimicking frequent storm events were discharged longitudinally or laterally over the pilot swales. The performance of each swale was assessed for 4 micropollutants, i.e. zinc (Zn), glyphosate, pyrene and phenanthrene. These substances were mainly found in the dissolved phase of the stormwater runoff used to supply the pilot swales. The standard swale, constructed from a silt loam soil, partially managed stormwater runoff by infiltration. Micropollutant concentration reductions were higher in the infiltrated water (35-85%) than in the overflow (-13-66%). The filtering swale, made of a sandy central part bordered by silt loam embankments, completely managed stormwater runoff by infiltration, providing high micropollutant concentration reductions (65-100%). Mass load reductions were higher for the filtering swale (67-90% for Zn and ≥89% for organic micropollutants) than for the standard swale (33-73% for Zn, 19-67% for glyphosate and ≥50% for both pyrene and phenanthrene). For both swales, lateral inflow was often associated with significantly higher concentration and mass reductions than longitudinal inflow. Consequently, when designing swales for the treatment of micropollutants, practitioners should preferentially promote filtering swales and installations providing lateral diffuse inflow over the facility.

The pollution conveyed by urban runoff: A review of sources

Urban stormwater and snowmelt pollution contributes significantly to the deterioration of surface waters quality in many locations. Consequently, the sources of such pollution have been studied for the past 50 years, with the vehicular transportation sector and the atmospheric deposition identified early as the major pollution sources. In search for mitigation of this pollution, source controls, besides other measures, were recognised as effective pollution mitigation tools, whose successful implementation requires a good knowledge of pollution sources. Even though great research efforts have been exerted to document specific sources of urban runoff pollution, or specific groups of pollutants present in urban runoff, a comprehensive overview of all known contributing sources is still missing. This review contributes to closing this gap by compiling findings of previous research and critically synthesizing the current knowledge of various stormwater pollution sources. As the emphasis is placed on the sources, the related issues of implications for urban surface water quality and possible source controls for individual sources are touched upon just briefly, where required. The review showed that the atmospheric deposition, vehicular transportation-related activities and metallic building envelopes continue to be among the major pollution sources, which have been studied in a far greater detail than other sources. Furthermore, it was noted that because of the rapid advances in clean manufacturing and pollution control technologies, a large part of the body of data on stormwater quality available in the literature should be considered as historical data, which may no longer describe well the current conditions. Progressing historical data obsolescence, combined with continuing releases of new materials and chemicals, and, in some cases of new substances of potential concern, into the environment, suggests that the identification of important stormwater runoff/snowmelt pollution sources, and the associated pollutants, has been and will remain to be a work in progress.

A preliminary pan-European assessment of pollution loads from urban runoff

Acknowledging the difficulty of modelling pollution conveyed by urban runoff, this contribution presents a first pan-European quantification of loads from this diffuse source. We estimate annual loads of 5-days biochemical oxygen demand (BOD₅), nitrogen (N), phosphorus (P) and total suspended solids (TSS) using a simple event mean concentration (EMC) model initially proposed by Heaney et al., 1976. On a European scale, this yields discharges corresponding to untreated wastewater of about 31 million population equivalents (PE) for BOD₅, about 18.5 million PE for N and P and about 280 million for TSS. These represent 51% of the pollution coming from treated wastewater for BOD₅, 15% for N and P and 461% for TSS. Although the model applied for the assessment was developed more than 40 years ago, the results are consistent with those obtained using more recent parameterizations, except for a tendency to underestimate P loads. Although lack of data on pollution from urban runoff makes model verification impossible, and the uncertainty on EMC models is known to be very high, urban runoff emerges as a significant source of pollution, and should be properly addressed as such. Reducing runoff volume from urban areas through improved water retention is not only key to pollution control, but also a no-regret option thanks to its co-benefits, especially when incorporated at early stages of planning and design.

Tracking faecal microorganisms using the qPCR method in a typical urban catchment in China

Faecal microorganisms represent a key threat to human health. Potential origins of faecal microbial contamination in a typical urban-representative micro-scale were evaluated. The quantitative polymerase chain reaction (qPCR) method was used in this study. The Bacteroidetes is selected as the indicative microorganism in runoff samples that are collected during four representative stormwater events in north China. The principal component analysis (PCA) method indicated the distribution feature of the environmental factors. The largest contributor is dog, followed by bird and human to the faecal pollution in stormwater runoff. The output of human and dog faecal pollutants in response to the first flush effect of nonpoint source pollution while the transmit time of bird faecal pollutant is relatively longer. In addition, the number of antecedent drying days represents the key factor for dog faecal pollution, while human faecal pollution is impacted by more factors. The results of this study will provide sound evidence for the tracking and management of nonpoint source faecal pollution in urban catchment areas.

Understanding spatiotemporal variability of in-stream water quality in urban environments - A case study of Melbourne, Australia

To support sustainable urban planning and the design of water pollution mitigation strategies, the spatial and temporal trends of water quality in urban streams needs to be further understood. This study analyses over ten years of surface water quality data from 53 upstream catchments (20 of them predominated by a single type of land use) and two lowland sites across Greater Melbourne, Australia. We evaluated the impact of various catchment characteristics, especially urban land uses, on spatial and temporal urban water quality trends. Here, we focused on common urban pollutants: total suspended solids (TSS), total phosphorous (TP), total nitrogen (TN), zinc (Zn), copper (Cu) and nickel (Ni). Site median nutrient and heavy metal concentrations were negatively correlated with the catchment's elevation and its average annual rainfall. Further analysis shows that such trends were driven by the geographical pattern of Melbourne - i.e. low-laying sites tend to have less rainfall and be more urbanised. Only median concentrations of heavy metals (Zn and Cu) were correlated to catchment imperviousness. Further characterising of the urban environment was done into specific land uses (residential, industrial and commercial), yet median concentrations of all pollutants were not significantly correlated with land uses. This is because simple metrics, such as land use proportions, do not adequately reflect the significant variability in pollution sources that can exist even within the same land use type. Indeed, our temporal analysis found that the water quality difference between catchments with similar land uses is likely caused by their site-specific pollutant sources (construction and illegal discharge) and environmental management actions (wastewater management actions) regardless of similarities in land use. A 3-stage urbanisation cycle (development, operation and renewal) is suggested to further explain the urban water quality variance, but more data from small areas of an urban catchment is required to directly understand the unique impact of each urbanisation stage on water quality.

Processes improving urban stormwater quality in grass swales and filter strips: A review of research findings

Increasing interest in urban drainage green infrastructure brings attention to grass swales and filter strips (GS&GFS) and their role in stormwater management. While the understanding of the hydrology and hydraulics of these stormwater control measures is adequate for current needs, there are knowledge gaps in understanding the water quality processes in GS&GFS and such a finding motivated preparation of the review paper that follows. The review revealed that most of the empirical studies of GS&GFS flow quality focused on the removal of pollutants associated with road runoff, and particularly solids, with relatively few studies addressing nutrients, traffic associated hydrocarbons, oxygen demanding substances, chloride, and faecal indicator bacteria. The reported results suffer from limitations caused by experimental conditions often representing a steady flow used to irrigate GS&GFS and generate runoff, non-submerged flows, no lateral inflows along swale side slopes, constant dosing of solids, emphasis on larger-than-typical solids, incomplete descriptions of experimental conditions, and limited attention to experimental uncertainties. Besides settling, other treatment processes, like adsorption/desorption, plant uptake, chemical precipitation and microbial degradation are often acknowledged, but without attempting to quantify their effects on flow quality. The modelling of GS&GFS flow quality would be beneficial for an improved understanding of green urban drainage infrastructure, but currently it is infeasible without a better knowledge of stormwater quality processes in GS&GFS facilities.

Comparison of two receptor models PCA-MLR and PMF for source identification and apportionment of pollution carried by runoff from catchment and sub-watershed areas with mixed land cover in South Korea

The application and comparison of receptor modeling techniques based on ambient air quality and particulate matter increasingly being studied. However, less is known about the comparison of receptor modeling techniques using spatial runoff quality data to identify and quantify the stormwater runoff pollution. This study compared the performance of principal component analysis-multiple linear regressions (PCA-MLR) and positive matrix factorization (PMF) models on stormwater runoff data collected from a small catchment (Site 1) with urban development activity and a sub-watershed outlet (Site 2). In both sites, the PCA-MLR model identified three pollution sources, whereas PMF identified five with a detailed source mechanism including two additional sources. Furthermore, the spatial land-use land-cover (LULC) analysis results indicate that the Site 1 exhibited a rapid conversion of the native area into a built-up area over the monitoring period compared to Site 2. Based on the modeling results, domestic wastewater and soil erosion were the major source of pollution at Site 1 and Site 2, respectively. The performance evaluation statistics including Nash coefficient (0.86-0.99), % error (<-14 to 2), and coefficient of determination (R² ≤ 0.99) showed better performance for the PMF model than the PCA-MLR model. Overall, the PMF receptor modeling approach was found to be more robust for the current study sites with different land use types. The findings of this study could provide a basis for further application of these receptor models and their comparison using spatial-temporal ionic and sediment related runoff monitoring data. Also, the models from this research could be combined with other receptor models on runoff quality data (e.g. CMB or UNMIX) to explore and inter-compare the outcomes, and to determine how the model results are affected by modifications to input data and model parameters. Therefore, further research is required to precisely assess the accuracy of both receptor models.

A Semi-distributed Model for Predicting Faecal Coliform in Urban Stormwater by Integrating SWMM and MOPUS

The microbial contamination of urban stormwater has an important impact on human health and stormwater reuse. This study develops an exploratory semi-distributed model, MOPUS_S, which can simulate faecal coliform levels in separate sewer systems in urban catchments. The MOPUS_S was built by coupling the SWMM model and the microbial MOPUS model. The parameters associated with the deposition and wash-off of microorganisms were more influential than those related to microorganism survival processes. Compared to other existing bacterial models, MOPUS_S showed comparable performance in predicting faecal coliform concentrations. The performance varied largely between rainfall events, with Nash-Sutcliffe efficiency (NSE) values ranging from −5.03 to 0.39 and R² ranging from −0.02 to 0.83, respectively. The model simulation results for low and medium concentrations were better than those for the peak concentrations. Poor simulation results of peak concentrations obviously affect the overall model performance. In general, MOPUS_S could be capable of predicting the faecal coliform load in urban catchments and be a useful tool for urban stormwater management planning.

Factors affecting microbial and physico-chemical pollutants in stormwater in a typical Chinese urban catchment

An understanding of microbial pollution characteristics is needed for stormwater reuse and development of microorganism simulations in urban stormwater. This study investigated the discharge characteristics of faecal indicator bacteria (faecal coliforms) in urban runoff by field sampling both the underlying surfaces and the stormwater pipe outlet. Faecal coliform contamination in urban runoff was found to be frequent, and the highest instantaneous concentration reached 2.42 × 106 MPN/100 ml. Faecal coliforms did not show a consistent first flush effect amongst the different surfaces sampled, and this was exacerbated under rainfall events with high intensity. PROMETHEE (Preference Ranking Organization Method for Enrichment Evaluation) and GAIA (Geometrical Analysis for Interactive Aid) analyses were further applied to explore the ranking of pollutants, the relationship among the pollutants, and the factors affecting the contamination in cases of multiple underlying surfaces, multiple pollutants and rainfall events. For the pollutants of suspended solids (SS), total phosphorus (TP) and chemical oxygen demand (COD), the road sample contamination was significantly higher than on the roof surfaces. No such trend in ranking of faecal coliforms was observed. Rainfall depth and intensity were found to have a significant influence on stormwater contamination by physico-chemical pollutants, while having a somewhat smaller influence on faecal coliform contamination. Faecal coliform contamination is closely associated with the index related to the antecedent dry period. The average temperature and average relative humidity also showed a positive relationship with faecal coliform contamination. The effects of antecedent dry period duration on contamination of physico-chemical pollutants and faecal coliform are completely opposite. Antecedent dry period duration was positively related to the contamination of physico-chemical pollutants, but negatively related to faecal coliform contamination. Therefore, three variables, i.e., antecedent dry period duration, average temperature and average relative humidity, might be used to model the survival/die-off of faecal coliform during the antecedent dry period.

Tire wear particles in the aquatic environment - A review on generation, analysis, occurrence, fate and effects

Tire wear particles (TWP), generated from tire material during use on roads have gained increasing attention as part of organic particulate contaminants, such as microplastic, in aquatic environments. The available information on properties and generation of TWP, analytical techniques to determine TWP, emissions, occurrence and behavior and ecotoxicological effects of TWP are reviewed with a focus on surface water as a potential receptor. TWP emissions are traffic related and contribute 5-30% to non-exhaust emissions from traffic. The mass of TWP generated is estimated at 1,327,000 t/a for the European Union, 1,120,000 t/a for the United States and 133,000 t/a for Germany. For Germany, this is equivalent to four times the amount of pesticides used. The mass of TWP ultimately entering the aquatic environment strongly depends on the extent of collection and treatment of road runoff, which is highly variable. For the German highways it is estimated that up to 11,000 t/a of TWP reach surface waters. Data on TWP concentrations in the environment, including surface waters are fragmentary, which is also due to the lack of suitable analytical methods for their determination. Information on TWP properties such as density and size distribution are missing; this hampers assessing the fate of TWP in the aquatic environment. Effects in the aquatic environment may stem from TWP itself or from compounds released from TWP. It is concluded that reliable knowledge on transport mechanism to surface waters, concentrations in surface waters and sediments, effects of aging, environmental half-lives of TWP as well as effects on aquatic organisms are missing. These aspects need to be addressed to allow for the assessment of risk of TWP in an aquatic environment.

Technical and natural conditions and operating efficiency of a municipal stormwater treatment plant

A decade of observations provided grounds for assessing the operation of one of the few stormwater treatment plants in Poland (system: screens-grit chambers-settler-retention pond) which collects effluents from 471 ha of the city. Among other aspects, the following were evaluated: treatment efficiency, relationship between the quality of treated stormwater and that of waters in the receiving body (the ox-bow lake of the Vistula river), operating stability of key units, significance of the facility for nature. During the assessment, the plant had a positive effect on the quality of stormwater effluents-the content of the analysed pollutants was reduced (more than 80% average efficiency for mineral forms of nitrogen and suspension matter) and oxygen ratios improved (23% increase in the average concentration of dissolved oxygen and more than 50% decrease in 5-day biochemical oxygen demand and chemical oxygen demand). Although the overall assessment of the facility's operation was good, some omissions and operating errors were noted (method of removing retained pollutants, stormwater flow control). Eliminating them is a prerequisite for maintaining the expected reliability of the system. An effect of stormwater ponds on the increase in biodiversity in the poor urbanised landscape has also been observed. The structures, forming a uniform system along with urban green areas, constitute specific enclaves which attract living organisms.

Total and settling velocity-fractionated pollution potential of sewer sediments in Jiaxing, China

Sewer sediments and their associated contaminant released along with wet-weather discharges pose potential pollution risks to environment. This paper presents total characteristics of sediments collected from Jiaxing, China. Size distribution and concentrations of volatile solids (VS) and four metals (Pb, Cu, Zn, Cr) of sediment samples from seven land use categories were analyzed. Then, the sediment samples were graded five fractions according to its settling velocity through the custom-built settling velocity-grading device. Sediment mass and pollution load distribution based on settling velocity were also assessed. The results show that there are relatively high level of heavy metal load in the sediment of separated storm drainage systems in Jiaxing, especially for the catchment of residential area (RA), road of developed area (RDA), and industrial area (IA). Although grain size follows a trend of increasing along with settling velocity, the methods of settling velocity grading are meaningful for stormwater treatment facilities with precipitation. For all land use categories, the pollution concentrations of the three lower settling velocity-fractionated sediment are relatively consistent and higher than others. Combined with mass distribution, the pollution percentage of fraction with different velocities for seven land use categories were also evaluated. Based on it, the statistical conclusion of design target settling velocity to different pollution load removal rates are drawn, which is helpful to guide design of on-site precipitation separation facilities.

Mineral and Anthropogenic Indicator Inorganics in Urban Stormwater and Snowmelt Runoff: Sources and Mobility Patterns

Inorganic chemicals in urban stormwater and snowmelt runoff originate from catchment geology and anthropogenic activities. The occurrence, partitioning and mobility of six minerals and six trace metal (TM) indicators of anthropogenic activities were studied in stormwater, snowmelt and baseflow in four urban catchments, and the sampling of inorganics was supplemented by measurements of electrical conductivity (EC), pH and total suspended solids (TSSs). Minerals occurred at concentrations several orders of magnitude higher (1-10² mg/L) than those of TMs (10^-2-10² μg/L) and reflected the composition of local groundwater seeping into sewers. Concentrations of Ca, K, Mg and Na were enhanced by baseflow contributions and followed closely the electrical conductivity. Al and Fe minerals occurred in insoluble forms, and their pollutographs were similar to those of TMs, whose concentrations mimicked, to some extent, the flux of TSS. The TMs with the highest and lowest particulate fractions were Cr&Pb and Cu&Zn, respectively. The concentrations of total TMs in snowmelt were two to four times higher than those in stormwater, and both sources likely exceeded some of the stormwater effluent limits (for Cd, Cu and Zn) proposed in Sweden. Where such concentrations depended on water hardness, the risk of toxicity might be reduced by elevated hardness of the monitored snowmelt and stormwater. Recognizing the good ecological status of the study area receiving water, Lake Storsjön, some protection against polluted runoff and snowmelt may be needed and could be achieved by implementing stormwater management measures controlling TSS and TMs.

Variability of zinc, copper and lead contents in sludge of the municipal stormwater treatment plant

Several years of observations enabled us to assess the extent and variability of heavy metal contamination (Zn, Cu and Pb) of the sludge of the municipal stormwater treatment plant (the town of Puławy, Poland). In terms of size (high capacity) and innovation, it is the only facility of this kind in the country. It collects rainwater run-offs from two catchments (separate inlets) with a total area of about 500 ha. The concentration of the analysed metals is characterised by a large spatial and temporal diversity. The reason for this may be differences in the manner and intensity of use of drained surfaces and different hydraulic conditions (of sedimentation) prevailing in the particular treatment devices. The highest pollution was found in sediments in the grit chamber and in the part of the settler from the side of the interceptor supplying sewage, i.e. from the main traffic route of the town (heavy traffic and developed technical infrastructure). The best-quality sludge was retained in the pond for treated wastewater. In the pool of analysed components, the largest share is Zn, which amounts to about 85%. The content of heavy metals limits the possibility of the natural use of sludge from the municipal stormwater treatment plant. In chemical terms, they should be seen as a potentially dangerous waste and undergo remediation.

Spatial distribution of heavy metals in the surface soil of source-control stormwater infiltration devices - Inter-site comparison

Stormwater runoff infiltration brings about some concerns regarding its potential impact on both soil and groundwater quality; besides, the fate of contaminants in source-control devices somewhat suffers from a lack of documentation. The present study was dedicated to assessing the spatial distribution of three heavy metals (copper, lead, zinc) in the surface soil of ten small-scale infiltration facilities, along with several physical parameters (soil moisture, volatile matter, variable thickness of the upper horizon). High-resolution samplings and in-situ measurements were undertaken, followed by X-ray fluorescence analyses and spatial interpolation. Highest metal accumulation was found in a relatively narrow area near the water inflow zone, from which concentrations markedly decreased with increasing distance. Maximum enrichment ratios amounted to >20 in the most contaminated sites. Heavy metal patterns give a time-integrated vision of the non-uniform infiltration fluxes, sedimentation processes and surface flow pathways within the devices. This element indicates that the lateral extent of contamination is mainly controlled by hydraulics. The evidenced spatial structure of soil concentrations restricts the area where remediation measures would be necessary in these systems, and suggests possible optimization of their hydraulic functioning towards an easier maintenance. Heterogeneous upper boundary conditions should be taken into account when studying the fate of micropollutants in infiltration facilities with either mathematical modeling or soil coring field surveys.

Model-based screening for critical wet-weather discharges related to micropollutants from urban areas

Wet-weather discharges contribute to anthropogenic micropollutant loads entering the aquatic environment. Thousands of wet-weather discharges exist in Swiss sewer systems, and we do not have the capacity to monitor them all. We consequently propose a model-based approach designed to identify critical discharge points in order to support effective monitoring. We applied a dynamic substance flow model to four substances representing different entry routes: indoor (Triclosan, Mecoprop, Copper) as well as rainfall-mobilized (Glyphosate, Mecoprop, Copper) inputs. The accumulation on different urban land-use surfaces in dry weather and subsequent substance-specific wash-off is taken into account. For evaluation, we use a conservative screening approach to detect critical discharge points. This approach considers only local dilution generated onsite from natural, unpolluted areas, i.e. excluding upstream dilution. Despite our conservative assumptions, we find that the environmental quality standards for Glyphosate and Mecoprop are not exceeded during any 10-min time interval over a representative one-year simulation period for all 2500 Swiss municipalities. In contrast, the environmental quality standard is exceeded during at least 20% of the discharge time at 83% of all modelled discharge points for Copper and at 71% for Triclosan. For Copper, this corresponds to a total median duration of approximately 19 days per year. For Triclosan, discharged only via combined sewer overflows, this means a median duration of approximately 10 days per year. In general, stormwater outlets contribute more to the calculated effect than combined sewer overflows for rainfall-mobilized substances. We further evaluate the Urban Index (A_{urban,impervious}/A_natural) as a proxy for critical discharge points: catchments where Triclosan and Copper exceed the corresponding environmental quality standard often have an Urban Index >0.03. A dynamic substance flow analysis allows us to identify the most critical discharge points to be prioritized for more detailed analyses and monitoring. This forms a basis for the efficient mitigation of pollution.

Source-Based Modeling Of Urban Stormwater Quality Response to the Selected Scenarios Combining Future Changes in Climate and Socio-Economic Factors

The assessment of future trends in urban stormwater quality should be most helpful for ensuring the effectiveness of the existing stormwater quality infrastructure in the future and mitigating the associated impacts on receiving waters. Combined effects of expected changes in climate and socio-economic factors on stormwater quality were examined in two urban test catchments by applying a source-based computer model (WinSLAMM) for TSS and three heavy metals (copper, lead, and zinc) for various future scenarios. Generally, both catchments showed similar responses to the future scenarios and pollutant loads were generally more sensitive to changes in socio-economic factors (i.e., increasing traffic intensities, growth and intensification of the individual land-uses) than in the climate. Specifically, for the selected Intermediate socio-economic scenario and two climate change scenarios (RSP = 2.6 and 8.5), the TSS loads from both catchments increased by about 10 % on average, but when applying the Intermediate climate change scenario (RCP = 4.5) for two SSPs, the Sustainability and Security scenarios (SSP1 and SSP3), the TSS loads increased on average by 70 %. Furthermore, it was observed that well-designed and maintained stormwater treatment facilities targeting local pollution hotspots exhibited the potential to significantly improve stormwater quality, however, at potentially high costs. In fact, it was possible to reduce pollutant loads from both catchments under the future Sustainability scenario (on average, e.g., TSS were reduced by 20 %), compared to the current conditions. The methodology developed in this study was found useful for planning climate change adaptation strategies in the context of local conditions.

Comparative analysis of water quality and toxicity assessment methods for urban highway runoff

In this study, comparative analyses of highway runoff samples obtained from seventeen storm events have been conducted between the traditional water quality assessment method and biotoxicity tests, using zebrafish (Danio rerio) embryos and luminous bacteria (Vibrio qinghaiensis. Q67) to provide useful information for ecotoxicity assessment of urban highway runoff. The study results showed that the Nemerow pollution index based on US EPA recommended Criteria Maximum Concentrations (CMC) (as traditional water quality assessment method) had no significant correlation with luminous bacteria acute toxicity test results, while significant correlation has been observed with two indicators of 72 hpf (hours post fertilization) hour hatching rate and 96 hpf abnormality rate from the toxicity test with zebrafish embryos. It is therefore concluded that the level of mixture toxicity of highway runoff could not be adequately measured by the Nemerow assessment method. Moreover, the key pollutants identified from the water quality assessment and from the biotoxicity evaluation were not consistent. For biotoxic effect evaluation of highway runoff, three indexes were found to be sensitive, i.e. 24 hpf lethality and 96 hpf abnormality of zebrafish embryos, as well as the inhibition rate for luminous bacteria Q67. It is therefore recommended that these indexes could be incorporated into the traditional Nemerow method to provide a more reasonable evaluation of the highway runoff quality and ecotoxicity.

Characterization and source identification of pollutants in runoff from a mixed land use watershed using ordination analyses

While identification of critical pollutant sources is the key initial step for cost-effective runoff management, it is challenging due to the highly uncertain nature of runoff pollution, especially during a storm event. To identify critical sources and their quantitative contributions to runoff pollution (especially focusing on phosphorous), two ordination methods were used in this study: principal component analysis (PCA) and positive matrix factorization (PMF). For the ordination analyses, we used runoff quality data for 14 storm events, including data for phosphorus, 11 heavy metal species, and eight ionic species measured at the outlets of subcatchments with different land use compositions in a mixed land use watershed. Five factors as sources of runoff pollutants were identified by PCA: agrochemicals, groundwater, native soils, domestic sewage, and urban sources (building materials and automotive activities). PMF identified similar factors to those identified by PCA, with more detailed source mechanisms for groundwater (i.e., nitrate leaching and cation exchange) and urban sources (vehicle components/motor oils/building materials and vehicle exhausts), confirming the sources identified by PCA. PMF was further used to quantify contributions of the identified sources to the water quality. Based on the results, agrochemicals and automotive activities were the two dominant and ubiquitous phosphorus sources (39-61 and 16-47 %, respectively) in the study area, regardless of land use types.

Water and nonpoint source pollution estimation in the watershed with limited data availability based on hydrological simulation and regression model

Nonpoint source (NPS) pollution is considered as the main reason for water quality deterioration; thus, to quantify the NPS loads reliably is the key to implement watershed management practices. In this study, water quality and NPS loads from a watershed with limited data availability were studied in a mountainous area in China. Instantaneous water discharge was measured through the velocity-area method, and samples were taken for water quality analysis in both flood and nonflood days in 2010. The streamflow simulated by Hydrological Simulation Program-Fortran (HSPF) from 1995 to 2013 and a regression model were used to estimate total annual loads of various water quality parameters. The concentrations of total phosphorus (TP) and total nitrogen (TN) were much higher during the flood seasons, but the concentrations of ammonia nitrogen (NH3-N) and nitrate nitrogen (NO3-N) were lower during the flood seasons. Nevertheless, only TP concentration was positively correlated with the flow rate. The fluctuation of annual load from this watershed was significant. Statistical results indicated the significant contribution of pollutant fluxes during flood seasons to annual fluxes. The loads of TP, TN, NH3-N, and NO3-N in the flood seasons were accounted for 58-85, 60-82, 63-88, 64-81% of the total annual loads, respectively. This study presented a new method for estimation of the water and NPS loads in the watershed with limited data availability, which simplified data collection to watershed model and overcame the scale problem of field experiment method.