1. Introduction
Regulated stormwater management aims to improve the ecological function of urbanized watersheds. To protect threatened watersheds in Wisconsin, the Wisconsin Pollutant Elimination System (WPDES) determines standards for stormwater management through a permitting process that sets total maximum daily loads (TMDL) for targeted pollutants such as total suspended solids (TSS) and total phosphorus (TP). It was recently reported that as of 2022, 1,526 bodies of water in the state were considered impaired [
1]. Nestled in the heart of Wisconsin, the University of Wisconsin-Madison campus borders more than a fifth of Lake Mendota, its largest water body in the county that provides critical ecosystem services to the surrounding watershed. Thus, any pollutant discharge from the campus compromises the integrity of the lake and surrounding terrestrial environments (
e.g., Lakeshore Nature Preserve) [
2]. To mitigate these potential effects, the university works to meet pollutant reduction goals as dictated by the YaharaWins project, a collaborative effort to meet WPDES permit regulations [
3,
4]. Specifically, the university seeks to reduce TSS pollution by 40% and TP pollution by 27% partially through the use of stormwater control measures [
3]. In 2021, progress toward meeting these goals was estimated at 32.76% and 25.3% reduction in TSS and TP, respectively [
5].
As a critical part of this effort, UW-Madison has incorporated green and gray infrastructure in its landscape design to help mitigate the impacts of stormwater [
3]. One such method includes the use of porous surfaces such as porous concrete, asphalt, or interlocking impervious pavers to promote infiltration and reduce runoff that would otherwise carry pollutants into surrounding water bodies [
6,
7,
8]. The Wisconsin Department of Natural Resources (WDNR) sets performance standards for porous pavements to infiltrate a minimum of 25 cm hr
-1 [
9]. This standard is tested using the Source Loading and Management Model for Windows (WinSLAMM) that estimates pollutant concentrations and runoff volumes from user-defined landscape parameters; the model can incorporate green and gray infrastructure, a particularly useful feature for stormwater monitoring over complex landscapes [
10,
11]. While the model provides a useful tool to evaluate stormwater pollution, it is rarely validated empirically as it is not required by the state government in many cases (
e.g. Wisc. Admin. Code Ch. NR 152) [
11,
12]. Moreover, the model estimates are only as good as the validity of the input parameters (e.g., existing landcover), which are especially challenging to quantify in heterogeneous urbanized landscapes exemplified by Madison, Wisconsin [
13,
14]. Thus, the long-term fidelity of these mitigation practices is at best unsystematically monitored and, at worst, unknown, a discrepancy that can compromise estimates of stormwater impact to surface waters and the surrounding terrestrial environments.
To test whether pollutant reduction values are being accurately modeled on the UW-Madison campus, we conducted a case study with two primary objectives. Our first objective was to empirically quantify total suspended solids (TSS), total phosphorus (TP), total Kjeldahl nitrogen (TKN), and chloride (Cl
-) concentrations in stormwater runoff from the UW-Madison campus to compare these observations to modeled values from WinSLAMM and international averages. These pollutants were chosen for their well-documented adverse effects on surface water health [
15,
16,
17,
18]. Excess TSS can diminish water quality via increased turbidity and increase the transport of pollutants such as phosphorus that can bind to the suspended sediment [
15]. Similarly, phosphorus and nitrogen loading can result in eutrophication and harmful algal blooms of surface waters [
16]. In addition, excess chlorides derived from road salts can impact aquatic organisms by disrupting their ability to osmotically regulate water [
18]. Our second objective was to quantify the infiltration capacity of four differently engineered porous surfaces used on campus. We expected the modeled results to deviate from field estimates for each pollutant measured due to temporal limitations of the model. We also expected that the infiltration capacity of all porous surfaces (asphalt, concrete, pavers with aggregates, and pavers without aggregates) would not meet WDNR infiltration standards and thus be overestimated by the WinSLAMM model.
3. Results
Stormwater samples were collected during eight rainfall events from August to September. The average precipitation event during that time averaged 1.02 cm and ranged from 0.15 cm to 3.33 cm during the time of sampling (
Figure 2). Over the entire sampling period, the concentration of TSS did not exceed international averages as reported by Shahzad et al. (2022) (
Figure 3a) [
25]. In contrast, TP concentrations did exceed international averages in four locations, most notably ME-11, and ME-26, which doubled this benchmark (
Figure 3c). Because international averages for TKN do not exist (to the best of our knowledge), the international average for TN as reported by Shahzad et al. (2022) was used for comparison; no TKN concentrations sampled in this study exceeded this benchmark (
Figure 3b) [
25]. Likewise, none of the outlets sampled for chlorides over the course of this study exceeded the benchmark concentrations for surface water toxicity as reported by the EPA and WDNR (
Figure 3d) [
26,
27].
The WinSLAMM-modeled averages of TSS, TP, and TKN were greater than the observed concentrations from the stormwater outlets sampled in this study. Observed TSS concentrations for every catchment area or group of catchment areas were less than modeled predictions. On average, observed TSS concentrations were 148 mg L
-1 lower than modeled concentrations with a standard deviation of 24.9 mg L
-1. Observed and modeled TP values were on average 0.0280 mg L
-1 apart with high variability; all modeled values exceeded the TP international average (
Figure 4c). The measured TP concentrations for WC-8 and the EC-9 catchment group were over 0.200 mg L
-1 less than their modeled concentrations. Similarly, the average measured TKN concentration for WC-8 was 0.700 mg L
-1 of its modeled concentration. Otherwise, observed TKN values were on average 0.394 mg L
-1 below modeled values (
Figure 4b).
Infiltration rates varied between the four porous pavement and paver sites measured in this study (
Figure 5). The infiltration rate for both porous pavers exceeded WDNR requirements for a porous pavement (0.42 cm min
-1) [
9]. The Carson Gulley porous pavers with aggregates exceeded this requirement with an average infiltration rate of 4.57 cm min
-1, as did the Dejope porous pavers without aggregates, but by a lesser amount (0.57 cm min
-1.) There was higher variability in the Dejope site with a standard deviation that exceeded the mean (0.92 cm min
-1). Neither porous pavement (porous asphalt or porous concrete) met WDNR requirements, although the porous asphalt did have a slightly higher infiltration rate than the porous concrete (0.25 and 0.13 cm min
-1, respectively).
4. Discussion
Our results suggest that TSS concentrations in stormwater runoff may be less of a concern compared to TP for the university within the area and for the period sampled here. The TP of several catchments exceeded the international benchmark, which is noteworthy in an already phosphorus-impaired watershed [
4]. However, the TKN and Cl
- concentrations were similar to other cold-climate municipalities and, although potentially detrimental to the watershed, were not at levels of concern over the period assessed in this study [
33]. The WinSLAMM model yielded similar TP and TKN concentrations compared to empirical values in this study, however it yielded TSS concentrations almost fourfold higher than empirical data. The porous pavers measured in this study were more effective at infiltrating stormwater compared to the porous pavements, which did not meet WDNR infiltration requirements.
The TSS concentrations observed here were well below the international average [
25] and nearly half of those from similar catchments around Madison, Wisconsin [
34,
35]. These differences may be explained from the lower proportion of impervious surfaces that characterize the university grounds sampled in this study compared to the city of Madison [
34,
35]. It was demonstrated with the use of geotextiles that urban areas with greater extent of pervious ground showed reductions in TSS concentrations specifically when that pervious ground is managed to prevent erosion [
36]. Because most pervious areas in the sampled catchments were comprised of well-maintained lawns or forested areas which can stabilize soil macroaggregates, we hypothesize that erosion was likely minimized overall reducing the potential TSS concentrations in these catchments [
37]. Alternately, we may have underestimated mean annual concentrations of TSS by the limited two-month sampling period, which may not capture the seasonal variability of this pollutant [
38,
39]. For example, concentrations of TSS have been shown to increase significantly at the beginning of the stormwater season (early May in our region), appearing as a “first-flush” effect, followed by a decrease as the summer progresses [
38,
39]. Adding to this potential error, samples taken in this study immediately after the rainfall events may not represent the mean TSS concentrations of the entire event. Research has demonstrated how rainfall characteristics can lead to differing TSS concentrations over the course of a stormwater event [
40]. Whether this disparity is the result of differences in land cover, the sampling methods used in this study, or both, the implication these data have for TSS loading warrants further investigation.
Both TP concentrations and its variability between catchments were greater for the area studied compared to international averages, exceeding those averages during many rainfall events. This variability may be attributed to differences in land use. Typically, higher phosphorus concentrations are associated with effluent water sourced from pervious areas receiving phosphorus fertilizer (e.g., lawns) or areas where organic litter can be eroded in stormwater events (e.g., disturbed construction areas) [
16,
41,
42]. Consistent with these drivers, the outlets in this study that drained more pervious catchment areas supporting lawn and urban trees displayed elevated levels of TP compared to outlets draining more impervious areas. Similarly, our measurements of TP were within the range of low-density residential areas encompassing lawns and gardens and commercial areas that incorporate green and gray infrastructure [
35,
41,
43]. Despite the observed elevated TP concentrations from the more pervious catchments, overall, TP concentrations from all land uses were similar to other studied urban areas [
35,
41,
43].
Both TKN and chloride concentrations measured in this study corresponded with concentrations in similar landscapes. The concentrations of TKN approached TN values of the more urbanized area of Madison, Wisconsin and an urbanized catchment of Columbus, Ohio (e.g. Smith et al. 2020). Research in Australia demonstrated that TKN represents a high percentage of TN, nearly 84% in some cases, which would make the TKN levels in this study comparable to the TN international benchmark and other research (
e.g. Lucke et al., 2018) reasonable [
44,
45]. As demonstrated in an urban catchment in Florida, a likely source of organic nitrogen was leaf litter and grass clippings, detritus was likewise observed in our samples [
44]. These comparisons suggest that campus management of vegetative detritus may be critical in mitigating stormwater TKN. Similarly, the chloride concentrations measured in this study were consistent with levels found in watersheds within cold climates requiring road salts [
39,
46]. As with TSS, a seasonal bias in sampling may explain these results, as chloride concentrations in runoff typically fall to their lowest levels in late summer and fall months, as documented in a study of urban stormwater in Helinski, Finland [
39]. Thus, we hypothesize that chloride concentrations would increase dramatically following snowmelt, a hypothesis worth testing specifically on large university campuses within temperate climates that feed large water bodies (such as the one here) that require ice abatement [
47,
48].
In this study, observed TSS concentrations were four times lower than the modeled values generated by WinSLAMM. These differences may be explained by three possible sources of error. First the aforementioned bias in sampling may have contributed to the discrepancy between modeled and observed values. Second, WinSLAMM modeled an average TSS concentrations over five years using parametrized precipitation values from a rainfile spanning 1980-1985, which likely varies from the two-month sampling average used in this study. Third, the overestimations of the WinSLAMM modeled TSS concentrations may be due to the model’s lack of appropriate parameters to represent an increasingly complex landscape that integrates stormwater management infrastructures [
34,
49]. For example, a study reviewing the role of particle size distribution in stormwater modeling found that particle size distribution can have a significant effect on modeled outputs of loading [
49,
50,
51], a physical aspect of the soil influenced by stormwater mitigation and not quantified in our study. Uncalibrated particle sizes coupled with inaccurate landcover data and parameterization for the modeled area may have synergistically combined to produce overestimated concentrations of solids present in effluent water observed in this study.
In contrast with TSS, the WinSLAMM model more closely characterized empirical estimates of stormwater TP and TKN concentrations in the university catchments. This alignment is somewhat counterintuitive given runoff with phosphorus and nitrogen is often sourced from suspended particles, which the WinSLAMM model utilizes in part to estimate TP and TKN concentrations [
16,
43,
49]. If TP and TKN were indeed accurately modeled, this suggests that TSS was likewise correctly modeled. Alternately, we hypothesize that the discrepancy between TSS and TKN and TP could be due to higher contributions of non-particulate bound sources of TP and TKN from litter or from the dissolved sources of these two pollutants. Studies have found that the non-particulate and dissolved fractions of nitrogen and phosphorus were significant contributors to their overall concentration in stormwater [
44,
52,
53]. This disparity supports investigation into sources of phosphorus and nitrogen and their relative contribution to stormwater generated from this setting.
As referenced earlier in the discussion, the differences between the modeled and field-measured results in this study may solely be explained by a temporal discontinuity [
38,
39] between the sampling period of this study relative to the five-year averages based on the 1980-1985 rainfile for Madison, WI used in the current WinSLAMM model. Because it was recommended that the WinSLAMM model is only suited for periods over a year [
54], we decided to test if modeled values of TSS would more closely align with our empirical estimates if the sampling period was constrained to the two-month sampling period used in this study (August – September). This comparison yielded TSS concentrations that were similar to the values generated by the five-year model; seasonally-modeled values for TSS exceeded the five-year model by an average of 16.1 mg L
-1 with a standard deviation of 7.59 mg L
-1. The seasonal model also exceeded empirical measurements of TSS concentrations by an average of 164 mg L
-1 with a standard deviation of 14.1 mg L
-1.
The infiltration rate of the porous pavers with aggregates outperformed the WDNR standard by over four cm min
-1, as well as the other three porous surfaces measured in this study, which did not meet this standard. We hypothesize that this difference may be explained by an interaction between the landcover type producing the influent stormwater, the age of the system, and the type of porous material used, all of which affect clogging between aggregates. In a study testing the change in infiltration and flow rate on a lab-simulated porous pavement, variable clogging was observed with changes in influent suspended particles changed used as a proxy for different landcover [
55]. This effect of the surrounding landcover may be enhanced as the porous system ages. Case studies have also shown porous pavers, porous concrete, and porous asphalt to clog at significantly different rates in as little as three years post instillation, although regular maintenance of these systems can minimize these differences [
55,
56]. An age difference may explain the higher infiltration rate of the porous pavers, the youngest systems by two years, however it is not sufficient to explain the differences between the other two, similar-aged paver systems or the higher infiltration rate of the older porous asphalt relative to the porous concrete. We speculate that the fourfold increase in infiltration rate between the two porous pavers of the same age may either be attributed to the inclusion of the aggregates in the interlocking gaps of the paver system. Assuming minimal differences in their immediate environments, the aggregates could have contributed to decreased rates of clogging leading to sustained infiltration rates relative to the other systems [
57]. Relative to our finding of porous asphalt, a recent study of porous pavements demonstrated that porous asphalt retains its infiltration rate more effectively than porous concrete [
56]. Although this study only assessed data over three years, we surmise that the porous asphalt either retained its pore connectivity since installation or its surrounding environment did not produce the solids necessary to clog its pores at a comparable rate to the porous concrete. Although strong inference cannot be generated from our case study, these results do suggest that the matching of the porous surface type to its site is essential to the long-term maintenance plan and mitigating effectiveness of these systems.