Sustainable Stormwater Management: Examining the Role of Local Planning Capacity in Mitigating Peak Surface Runoff

The Chesapeake Bay, the largest estuary in the United States, is rich in natural resources. Its watershed has been impacted by excessive and degraded stormwater runoff from rapid urbanization. We used an empirical approach to investigate how local planning capacity in the Chesapeake Bay watershed affected stream flow. A multiple regression analysis was employed to examine to what extent that the planning factors and other contextual variables were associated with peak runoff. Counterintuitively, we found that sub-basins included in the sample jurisdictions with a relatively high plan quality score tend to generate higher volumes of peak runoff. Results further indicate that specific geographical, basin characteristic, and biophysical factors affected mean annual peak runoff significantly. Overall, our findings highlight the importance of local planning capacity and sustainable stormwater management concepts in mitigating excessive runoff.


Introduction
The occurrence of excessive runoff and flooding events is increasing in the United States due to rapid urbanization and aging stormwater infrastructure.According to the most recent U.S. Census, from 1950 to 2010, urbanized areas expanded by almost 210 percent, and population in urban areas increased by more than 130 percent.Land consumption rate is outpacing the population shift from urban areas to suburban areas [1].The ability of nature to respond to change has decreased due to rapid urbanization and urban sprawl.Conventional low-density development patterns, which caused environmental degradation, has significantly enlarged the area taken up by impervious surfaces, and thus facilitated landscape fragmentation, habitat displacement, and flood risks [2,3].The influences of land use changes, such as urbanization and deforestation, led to the rising increment of stormwater runoff volume and pollution [4,5].Previous studies [6][7][8] have discovered that increased impervious surfaces caused by urbanization generate negative hydrologic consequences, including excessive overflow, lack of infiltration, and insufficient aquifer recharge.Downstream water pollution and flooding have been exacerbated because of the early stormwater runoff system design and aging pipeline infrastructure.Specifically, conventional stormwater management approaches have focused on removing stormwater as promptly as possible in order to mitigate impacts from flooding in a particular subdivision [9].Hence, old pipeline drainage systems have increased the volume and velocity of runoff as well as peak flows, which incur greater danger to downstream water bodies in the form of flooding [9].Maintenance and replacement costs for these pipelines are relatively expensive compared to other on-site management systems such as Best Management Practices (BMPs) and Low Impact Development (LID) techniques [10,11].Unfortunately, the majority of local jurisdictions have historically paid little attention to stormwater management related infrastructure, and funding has been limited by regional and state governments compared to other governmental infrastructure activities such as road and land construction, which are classified as mainstream works [12,13].
In sum, these two problems are significant issues resulting in excessive runoff and will become more problematic as they continue to disturb the hydrological cycle and increase flood damage.Effective control and regulation in the early phases of development can help forestall or resolve these issues.Planning includes diverse planning processes, incorporating the active participation of various stakeholders.The decision-making processes before development provides local governments an opportunity to more effectively and comprehensively address runoff issues by embracing a wide range of goals toward sustainable stormwater management.In addition, planning is a process directed by a plan document that must be a long-range blueprint for a community's future development [14].Thus, incorporating stormwater management policies while adopting a plan may play a critical role in establishing stormwater management strategies for implementation in the initial stage and help effectively minimize adverse impacts from flooding and overflow.Most importantly, since many factors relating to stormwater runoff-such as rapid urbanization, urban sprawl, and inadequate drainage systems-are at the local level, the role of local land use decision making is becoming more crucial in managing stormwater and reducing excessive runoff [9,15].
Local governments are responsible for land use planning; they guide and regulate various urban environments and developments that may directly affect the stormwater system.Therefore, stormwater management should be addressed in the local comprehensive plan to proactively prepare for future flood risks and manage stormwater in a manner incorporated into larger concepts such as hazard, environmental, and ecosystem planning.
For almost two decades, local planning instruments have been evaluated in relation to the aspects of resilience, natural hazards, climate change, sustainability, smart growth, urban sprawl, citizen participation, green infrastructure, ecosystem management, and environmental planning [1,[15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32].However, no studies have examined whether the concepts, policies, and strategies of sustainable stormwater management are incorporated into planning documents.In addition, although the relationships among various factors and surface runoff have been examined in the past [2,7,[33][34][35][36], few studies [6] have thoroughly explored the effects of planning capacity on flood mitigation.Given these gaps in the previous research, this study is intended to assess the impact of local planning capacity on mean annual peak runoff.Seventy-five sub-basins within the Chesapeake Bay watershed were selected for investigation.The results will provide valuable information to nearby local decision-makers and watershed planners on how to improve their stormwater management planning and plan documents that may effectively minimize the volume of runoff and enhance the overall health of the bay area in the long term.

Planning Capacity and Its Impact on Runoff
Internal planning capacity refers to features that can be controlled by local governments, such as the planning resource and process as well as the institutional capacity.It has been known to be a key criteria in examining whether the function of local governments work well or not [37].Local jurisdictions with high commitment to local planning will likely to have stronger awareness and ability to mitigate stormwater runoff even though limitations still exist to verifying the degree to which planning factors affect the implementation of plan parameters in practice.Planning capacity investigated in this study included plan quality score, plan adopted year, number of planning staff, and involvement of consultant because they are considered most related to creating planning documents.
Plans that incorporate stormwater management mechanisms, such as non-structural tools (e.g., regulations on land use, taxes, site design, building codes, and public participation and education programs) and structural tools (e.g., LID practices, BMPs, and green infrastructures), tend to have local governments with higher commitment to controlling runoff.Brody and Highfield [24] used 18 plan quality indicators to assess whether wetland permit clusters in Florida conform to the original designs of comprehensive plans.They identified that plan quality scores of specific environmental and implementation policies had significant correlations with the degree of plan implementation (e.g., wetland development).Nelson and French [38] discovered that seismic safety elements within local comprehensive plans may have a positive effect in minimizing earthquake damage.Kang [39] found that plan quality scores of flood mitigation policies were positively associated with insured flood losses, even though the coefficient was statistically insignificant.
When localities have more resources and expertise, higher-quality plans can be generated, and thus specific policies have better chances to be implemented [21,39,40].Several previous studies [40][41][42] underscored planning staffs because of their crucial role in mitigating hazards, especially regarding flood damage.Brody et al. [1] found that jurisdictions with more planning agency staff had stringent sprawl-mitigation measures in their local comprehensive plans.Tang and Brody [32] discovered that more planners in the staff would contribute to higher-quality local environmental plans.Furthermore, more recently updated plans are likely to include up-to-date information, naturaland built-environmental conditions, and techniques, and thus they may enable local governments to develop better plans and encourage their implementation.Hiring private consultants may also bring more technical and human resources to the table with which to improve plan quality and facilitate implementation.Some studies further examined the impacts of internal planning factors such as budget, collaborative effort, planners' commitment, participation, and leadership on plan outcome [24,38,39].While this study did not include surveys with planning staff or community leaders, acquiring this information may better represent the entire local planning efforts and capacities; and thus, help explain the variations of surface runoff generation.
In this study, we tested the following hypotheses: (1) Sub-basins with plans of higher quality have lower peak runoff; (2) sub-basins that are included in jurisdictions that have recently adopted the plan, are less likely to generate excessive runoff; (3) sub-basins that are included in jurisdictions with more planners while drafting a local plan, will generate less peak runoff; and (4) sub-basins that are included in jurisdictions that engage private consultants for drafting a local plan will generate less peak runoff.

Conceptual Model
In this study, a conceptual model was developed with two phases to examine the effects of four specific factors (planning capacity, geographical, basin characteristics, and biophysical variables) on mean annual peak runoff (Figure 1).In Phase 1, plan quality score was derived by evaluating whether local jurisdictions in the sample sufficiently integrate the key principles of sustainable stormwater management into local comprehensive plans.A coding protocol for the plan evaluation was developed through the review of the literature associated with stormwater management.In Phase 2, the variance of mean annual peak runoff was analyzed by conducting multivariate regression analysis with four specific variables.original designs of comprehensive plans.They identified that plan quality scores of specific environmental and implementation policies had significant correlations with the degree of plan implementation (e.g., wetland development).Nelson and French [38] discovered that seismic safety elements within local comprehensive plans may have a positive effect in minimizing earthquake damage.Kang [39] found that plan quality scores of flood mitigation policies were positively associated with insured flood losses, even though the coefficient was statistically insignificant.When localities have more resources and expertise, higher-quality plans can be generated, and thus specific policies have better chances to be implemented [21,39,40].Several previous studies [40][41][42] underscored planning staffs because of their crucial role in mitigating hazards, especially regarding flood damage.Brody et al. [1] found that jurisdictions with more planning agency staff had stringent sprawl-mitigation measures in their local comprehensive plans.Tang and Brody [32] discovered that more planners in the staff would contribute to higher-quality local environmental plans.Furthermore, more recently updated plans are likely to include up-to-date information, natural-and built-environmental conditions, and techniques, and thus they may enable local governments to develop better plans and encourage their implementation.Hiring private consultants may also bring more technical and human resources to the table with which to improve plan quality and facilitate implementation.Some studies further examined the impacts of internal planning factors such as budget, collaborative effort, planners' commitment, participation, and leadership on plan outcome [24,38,39].While this study did not include surveys with planning staff or community leaders, acquiring this information may better represent the entire local planning efforts and capacities; and thus, help explain the variations of surface runoff generation.
In this study, we tested the following hypotheses: (1) Sub-basins with plans of higher quality have lower peak runoff; (2) sub-basins that are included in jurisdictions that have recently adopted the plan, are less likely to generate excessive runoff; (3) sub-basins that are included in jurisdictions with more planners while drafting a local plan, will generate less peak runoff; and (4) sub-basins that are included in jurisdictions that engage private consultants for drafting a local plan will generate less peak runoff.

Conceptual Model
In this study, a conceptual model was developed with two phases to examine the effects of four specific factors (planning capacity, geographical, basin characteristics, and biophysical variables) on mean annual peak runoff (Figure 1).In Phase 1, plan quality score was derived by evaluating whether local jurisdictions in the sample sufficiently integrate the key principles of sustainable stormwater management into local comprehensive plans.A coding protocol for the plan evaluation was developed through the review of the literature associated with stormwater management.In Phase 2, the variance of mean annual peak runoff was analyzed by conducting multivariate regression analysis with four specific variables.

Study Sample
The target population of this study is sub-basins within the Chesapeake Bay watershed.The Chesapeake Bay is the largest and most biologically diverse estuary in North America located in the Mid-Atlantic region [43].The watershed covers approximately 166,000 km 2 and a total of 203 counties and independent cities lie within or adjoining the bay watershed (see Figure 2).The study area has historically been polluted by human developments and impervious surfaces accompanied by rapid population growth.The population in the watershed has doubled between 1950 and 2000 (from 8 to 16 million).This growth has contributed to an impaired bay ecosystem, including habitat loss and water quantity/quality degradation [44].Approximately 15 percent of the total nitrogen entering the bay originates from urban and suburban polluted runoff, which has been recently recognized as the greatest threat to bay water quality [45].
the Mid-Atlantic region [43].The watershed covers approximately 166,000 km and a total of 203 counties and independent cities lie within or adjoining the bay watershed (see Figure 2).The study area has historically been polluted by human developments and impervious surfaces accompanied by rapid population growth.The population in the watershed has doubled between 1950 and 2000 (from 8 to 16 million).This growth has contributed to an impaired bay ecosystem, including habitat loss and water quantity/quality degradation [44].Approximately 15 percent of the total nitrogen entering the bay originates from urban and suburban polluted runoff, which has been recently recognized as the greatest threat to bay water quality [45].
The study area for this research was chosen based on the following steps.First, local jurisdictions that overlap with the Chesapeake Bay watershed boundary by more than 50 percent were selected to avoid the jurisdictions that may not directly influence the entire watershed ecosystem; Second, the sample was limited to jurisdictions with populations greater than 10,000 to prevent skew toward small jurisdictions, where areas exert little influence on the bay and often lack the resources to initiate a sufficient planning effort [19].Third, jurisdictions that adopted comprehensive plans between 2000 and 2010 were selected to determine the implementation effect of planning capacity factors on mean annual peak runoff from 2011 to 2014.The interpretation of results, however, should be made carefully because the lag time between plan implementation and hydrologic change can be large and we cannot guarantee whether the plan has been implemented; Finally, sub-basins that overlap with the boundary of a specific jurisdiction by more than 80 percent were chosen for the final sample, in order to represent the planning factors where the unit of analysis is at the county level.Through the above selection process, a total of 42 local jurisdictions and 75 sub-basins were contained in the sample.

Unit of Analysis
The unit of analysis for this study is sub-basin.The sample sub-basins have been delineated based on stream gauge data from the United States Geological Survey (USGS) [46] by following three The study area for this research was chosen based on the following steps.First, local jurisdictions that overlap with the Chesapeake Bay watershed boundary by more than 50 percent were selected to avoid the jurisdictions that may not directly influence the entire watershed ecosystem; Second, the sample was limited to jurisdictions with populations greater than 10,000 to prevent skew toward small jurisdictions, where areas exert little influence on the bay and often lack the resources to initiate a sufficient planning effort [19].Third, jurisdictions that adopted comprehensive plans between 2000 and 2010 were selected to determine the implementation effect of planning capacity factors on mean annual peak runoff from 2011 to 2014.The interpretation of results, however, should be made carefully because the lag time between plan implementation and hydrologic change can be large and we cannot guarantee whether the plan has been implemented; Finally, sub-basins that overlap with the boundary of a specific jurisdiction by more than 80 percent were chosen for the final sample, in order to represent the planning factors where the unit of analysis is at the county level.Through the above selection process, a total of 42 local jurisdictions and 75 sub-basins were contained in the sample.

Unit of Analysis
The unit of analysis for this study is sub-basin.The sample sub-basins have been delineated based on stream gauge data from the United States Geological Survey (USGS) [46] by following three sampling processes.First, a gauge that has its outlet located within a reservoir or has a dam on the upstream was excluded from the sample since the data can be impacted by storage capacity.Second, only gauges that have streamflow records between 2011 and 2014 were chosen for the final study, in order to examine the implementation effects of local plans that were adopted from 2000 to 2010.Third, for data efficiency and accuracy only gauges that have at least 90 percent of streamflow records per year were selected [34].
By using StreamStats, a Web-based GIS application that was developed by the USGS and Environmental Systems Research Institute (ESRI) for water resources planning and management, a distinct sub-basin boundary from each gauge station was delineated.

Dependent Variable
Mean annual peak runoff of 75 sub-basins from water years of 2011 to 2014 were obtained from the USGS gauge stations.Because stations provided the annual peak discharge rates for each sub-basin with the unit of cubic meter per second (m 3 /s), this study converted the flows into total annual runoff depth (in millimeters).Specifically, the converting method that the USGS applied for its estimation was employed.First, 86,400 s per day was multiplied to convert the value into a total annual flow volume (m 3 ).Second, runoff volume expressed in depth was computed by dividing the total annual flow volume by the contributing drainage area measured by ArcGIS.Third, meter measurement has been converted into the millimeter measurement by multiplying 1000.The mean annual peak runoff from each gauge station was obtained by averaging four-year (2011-2014) annual peak runoff.To better approximate a normal distribution, mean annual peak runoff was log-transformed (see Table 1).

Independent Variables
For examining the plan implementation process, several studies have adopted plan quality score as a causal variable.Through employing the content analysis methodology that is widely used as an evaluation protocol in plan assessments [19,23,47], this study conceptualized local plan quality on sustainable stormwater management, based on five key plan components: (1) Factual basis; (2) goals and objectives; (3) inter-organizational coordination; (4) policies, tools, and strategies; and (5) implementation.Sixty-two indicators were developed in evaluating local comprehensive plan by referring the concepts of sustainable stormwater management that were established in previous research [48][49][50][51][52][53][54][55][56][57] and various federal, state, and local water resources and stormwater management and planning guidelines [58][59][60][61][62][63][64][65].Table 2 shows the detailed indicators, measurements, and descriptive statistics for each plan component.Total plan quality scores for each jurisdiction was measured using Equations ( 1) and (2) (see Brody [47] for more details of the calculation process) and they are shown in the Appendix A (Table A1).
where PCQ j refers to the quality of the j-th plan component; m j refers to the total number of indicators within the j-th plan component (scale: 0-10); and I i refers to the i-th indicator's scores (scale: 0-2; scale for the "goals and objectives" component: 0-1).
where TPQ refers to the total plan quality scores (scale: 0-50).The measurement of total plan quality score followed the procedure of previous plan evaluation studies and they are computed by four steps.

•
Step 1: Score each indicator (scale: 0-2) within a plan component and add them all to gain total plan component score.Indicators were coded on a 0-2 ordinal scale.Indicators for the "goals and objectives" component, however, have been scored on a 0-1 scale.Specifically, an indicator scored two points when it was fully identified and demonstrated within a plan.If an indicator was explained or identified without a detailed description, it received 1. Zero points were given to an indicator when it was never mentioned within a plan.For "policies, tools, and strategies" component, an indicator scored 2 when it was clearly mentioned with a firm commitment words, such as "require," "must," "shall," and "will."Score of 1 was received when an indicator was portrayed with vague commitment words (e.g., "encourage," "should," "may," and "consider").When an indicator was specified but was not described with detailed information (e.g., "what," "where," how," and "when"), it received one point.

•
Step 2: Each plan component was standardized by dividing the total indicator scores within a component by the total available scores of a component.

•
Step 3: Multiply each plan component score by ten in order to make a 0-10 scale.
To maintain an inter-coder reliability and reduce personal bias in judgment, two trained scorers have evaluated all 42 local comprehensive plans.The plan indicators were pre-tested by the first scorer and re-tested by the second scorer using the same plan coding protocol.The percent agreement score, which is a generally accepted technique to measure inter-coder reliability in past plan evaluation studies, was computed through "ReCal", a Web-based tool [66].The overall average percent agreement score calculated from the double-coded data was about 84 percent.Generally, past plan quality evaluation studies considered a score higher than 80 percent as acceptable [67,68].
To examine the level of inter-item consistency and reliability, Cronbach's Alpha test, which assesses the degree to which a set of indicators are correlated as a group, was conducted in this study.The Cronbach's alpha exceeded 70 percent for all five plan components.An α value in the range of 70 percent or above is typically considered as an adequate reliability by many researchers [69,70].
The plan adoption year data were computed by subtracting the year that a plan was adopted from the year 2010.Data on the number of planners and the participation of consultants while drafting a plan were obtained from each local jurisdiction's comprehensive plan.Individual contacts have been made with local planning department officials where sufficient information was not provided within a plan.
Land use/land cover dataset (2011) was obtained from the USGS National Land Cover Database (NLCD) [71] at a 30 m resolution.Developed areas were represented by grouping three land use/land cover classes (LULC Class: 22-24): Low-intensity, medium-intensity, and high-intensity developed areas.These intensities were classified based on the percentage of impervious cover, and each comprises 21-49 percent, 50-79 percent, and 80-100 percent of impervious surfaces, respectively.Land uses for low-and medium-intensity developed areas are typically single-family housing, whereas high-intensity developed areas generally contain apartment complexes and commercial/industrial facilities [72].Wetland areas were represented by two LULC classes (woody wetland and emergent herbaceous wetland; LULC Class: 90-95).The percentages of LULC distribution were calculated by ArcGIS [73] with the Geospatial Modelling Environment (GME) extension [74].
Both mean slope and basin shape were measured by using ArcGIS.Specifically, mean slope was calculated based on the 30 m resolution DEMs obtained from the National Hydrography Dataset (NHD) Plus Version 2 [75].From the several basin shape measurements, such as circularity ratio, length to width ratio, and elongation ratio, this study employed the elongation ratio approach, which is frequently used in recent hydrological research.The value of elongation ratio was attained through calculating Equation (3).
where A refers to the basin area; and L refers to the basin length from the gauge station to the farthest point within a basin boundary.Four biophysical factors that may directly/indirectly influence the quantity of stormwater runoff are included in this variable: average monthly precipitation, natural drainage density, percentage of 100-year floodplain, and soil characteristics.Average monthly precipitation data were acquired from the PRISM Climate Group for the period from 2011 to 2014.The PRISM Climate Group produced a continuous record of surface precipitation by using the Climatologically-Aided Interpolation (CAI) approach.Each basin's average monthly precipitation was summed over the water year (1 October to 30 September) and was measured for the study period using ArcGIS to calculate average weighted mean of raster data.Natural drainage density was measured using ArcGIS with the national hydrography dataset obtained from the USDA's GeoSpatial Data Gateway [76].The ratio of total stream length to basin area was calculated.The digital flood insurance rate map (DFIRM) and Q3 data were obtained from the FEMA Map Service Center [77] to calculate the percentage overlapping a FEMA-defined 100-year floodplain with the basin area.To obtain the saturated hydraulic conductivity (K sat ) value, which is often used in soil interpretation, the Soil Survey Geographic Database (SSURGO) was obtained from the USDA's Web Soil Survey [78] and run using the Soil Data Viewer 6.1.Average K sat value of each sub-basin was then created using ArcGIS GME extension to weight the value according to the proportional areas.

Data Analysis
An ordinary least squares (OLS) technique was used to test how the independent variables (planning capacity, geographical, basin characteristics, and biophysical factors) explain the variance of dependent variable.Due to the relatively small sample size (n = 75) compared to the number of independent variables, variables were analyzed by four block groups.This approach was frequently adopted by several studies that had small number of sample sizes in order to alleviate the impact of each variable on the validity of statistical conclusion [30,32,39,79,80].Five models have been analyzed in this analysis.Specifically, Model 1 (baseline model) included only the block group of planning capacity variables.Geographical, basin characteristics, and biophysical variables were then added one by one to create the next models (Models 2-4).Only statistically significant variables in each of four models were then chosen for the final fully specified model (Model 5).By following Equation (4), multiple regression analyses were conducted.
where MAPR refers to mean annual peak runoff; α refers to regression intercept; β x refers to partial regression coefficients; X 1 refers to planning capacity factors; X 2 refers to geographical factors; X 3 refers to basin characteristics factors; and X 4 refers to biophysical factors.
To ensure that OLS regression assumptions were not violated and to check whether the OLS would yield best, linear, and unbiased estimates, this study tested model specification, outliers, multicollinearity, heteroskedasticity, and spatial autocorrelation.No major violations were detected through the diagnostics.

Results
We examined the influence of planning capacity, geographical, basin characteristics, and biophysical variables employing multiple regression analysis.Table 3 reports both multivariate regression coefficients and standardized coefficients for mean annual peak runoff.The variance of the dependent variable was explained the most by Model 5 (65 percent), followed by Model 2 (46 percent), Model 4 (42 percent), Model 3 (35 percent) and Model 1 (25 percent).
With respect to the association between planning capacity variables and peak runoff, plan quality score was constantly positive and significant in all models except Model 2. Its degree of coefficients, however, was relatively weak to explain the variance of the dependent variable, and directions were opposite to our initial expectation that sub-basins of higher plan quality would generate less peak runoff.Involvement of consultants was negative and statistically significant only in Models 1 and 3.Both plan adopted year and number of planners were statistically insignificant in all models.
For geographical variables, impervious surface had a positive and statistically significant relationship with mean annual peak runoff in Models 2 and 5.For 1 percent increase in impervious surface, peak runoff can be increased by approximately 1.7 percent.Specifically, impervious surface was the most powerful predictor in explaining the variance of mean annual peak runoff (Beta = 0.4634).Wetland was negatively associated with mean annual peak runoff in Model 2, but did not show a significant result in the fully-specified model.Although no serious multicollinearity was detected in Model 5, relatively high correlation (r = 0.53) between wetland and soil might be a possible reason for a reduced statistical effect on mean annual peak runoff.
Among the two basin characteristics variables, only basin shape (elongation ratio) was positively and significantly associated with mean annual peak runoff in Models 3 and 5.One unit increase in elongation ratio increase more than double, in the peak runoff.Average basin slope, however, was not statistically significant while the direction followed the expected signs.
Biophysical factors were highly related to runoff generation.Average monthly precipitation, the percentage of floodplain, and saturated hydraulic conductivity (soil) had statistically significant effects on mean annual peak runoff in Models 4 and 5. Specifically, sub-basins with a high percentage of floodplain were more likely to generate less peak runoff.Saturated hydraulic conductivity (soil) displayed a negative relationship with peak runoff.This result supports that sub-basins containing a higher percentage of permeable soils are likely to generate less mean annual peak runoff.Natural drainage density had a positive but statistically insignificant relationship with the dependent variable.

Discussion and Policy Implications
The explanatory results of this study have revealed several facets that are worth further consideration in minimizing surface runoff.First, local plan quality score was positively associated with mean annual peak runoff.This result suggests that possessing a high-quality plan does not always result in minimizing surface runoff.Perhaps this relationship may stem from several reasons.Although a jurisdiction develops a thorough comprehensive plan incorporating various policies and action strategies associated with stormwater management, those policies may not be implemented in practice.Several indicators used for the plan evaluation in this study may also be difficult for local planners and administrators to measure whether the implementation has occurred.Thus, we recommend that state and local agencies develop a plan implementation evaluation system that assesses whether plan outcomes conform to the initial intent of a plan.For example, they might adopt the methodology that Laurian et al. [81,82] used to identify whether land development permitting processes followed a plan's development policies.Such plan implementation evaluation systems may also play an important role in regular plan updates by discerning how certain policies have been actually implemented.
In addition, jurisdictions that have frequently experienced damages due to flooding or excessive runoff may have already recognized their vulnerability to stormwater runoff, and thus have integrated diverse stormwater management policies and tools into their comprehensive plans beforehand.If this is the case, plan quality score will be a reactive measure to the previous flooding experiences.As a result, even though a sub-basin is included in a jurisdiction that has a high-quality plan, the sub-basin may generate more runoff compared to sub-basins that generate less stormwater runoff historically.
The scores obtained from the plan evaluation process might also be a reason for the unexpected result.Although we have used the five plan components approach that is most commonly employed in the previous plan evaluation research, considering other components, such as "monitoring", "public participation," and "organization and presentation" may allow more number of practical indicators to be utilized for the assessment and alter the overall scores [83].Moreover, while indicators here were developed based on concerning several literature, guidelines, and opinions from experts, weighting schemes were not applied.Assigning unequal weights and values to different indicators based on multiple practitioners and planning researchers' judgment may lead to the different result.
Lastly, the volume of peak runoff may be significantly influenced by upstream human disturbances.Although the sampled sub-basins were delineated based on the topography and flow direction and accumulation, upstream development pressures may considerably impact the quantity and quality of interconnected downstream flow.Therefore, degraded upstream sub-basins will impact downstream ones on peak runoff and flooding.For this reason, local planners should actively cooperate with the upstream and nearby jurisdictions in regulating peak runoff even though they have developed a thorough comprehensive plan towards mitigating floods.
Because of the above interpretations, further research in the relationship between plan quality and surface runoff generation is needed.Comprehensive plan is a long-range policy document that guides a community's future development.It takes efforts to implement a plan and time to observe the outcomes.Communities that consider only immediate and short-term concerns of mitigating stormwater runoff without long-term visions, goals, objectives, and action strategies will fail to manage stormwater sustainably.Thus, local planners and decision-makers should continuously monitor whether specific policies have been successfully implemented and then determine if those implemented policies were effective.
Second, despite the unexpected result, the directions of plan updated year versus peak runoff and the participation of consultants versus peak runoff support the findings of previous studies: (1) The latest information and circumstances should be regularly included in updated plans and, thus, encouraging local planners to apply up-to-date techniques within the action strategies; (2) For small jurisdictions, assistance from outside consultants may also have positive influence in controlling runoff by adding more planning resources, analyses, and GIS resources into a plan.
Third, the percentage of floodplain had a negative association with peak runoff.This may be because land developments within the 100-year floodplain were well regulated by local governments, and thus the amount of excessive runoff might be minimized.
Fourth, impervious surfaces, which accounted for an average of 21.6 percent of land cover in the study area, were highly associated with peak runoff: a one percent increase in the impervious surface resulted in approximately 1.7 percent increase in mean annual peak runoff, holding other variables constant.That is, sound land regulations and strategies should be made to reduce the impervious surfaces during the development process.Local/state planners and agencies are strongly recommended to manage and monitor the spatial distributions and configurations of impervious land cover in order to effectively control the excessive runoff that might be caused by indiscriminate land development [84].In addition, BMPs and LID practices should be installed in places where the percentage of impervious surface is too high to efficiently manage runoff and prevent flooding events.These proactive planning approaches may also lead local governments to save initial construction and maintenance costs because on-site source control practices are more cost-effective than conventional drainage systems [85].

Conclusions
Although the relationship between land use and stormwater quantity has been researched for a while, this study is unique by utilizing planning capacity linked with hydrologic measures.In contrast to large-scale flooding caused by hurricanes and extreme rainfalls, stormwater runoff and/or flooding can be more effectively controlled at the local level.This study found that local planning and other contextual factors may significantly influence mean annual peak runoff in the Chesapeake Bay watershed.By far, a point increase in plan quality score would increase the peak runoff, inferring that the majority of local governments may already recognize the significance of stormwater runoff and flood occurrence, and thus substantially incorporate concepts of sustainable stormwater management in their plans.The link between comprehensive planning and successful stormwater management, however, is still tenuous and can be influenced by other factors that are not well represented in this study.Thus, interpretation of plan implementation results should be made carefully.Findings also support the pattern of previous plan evaluation studies that continuous and regular plan updates as well as more human resources may contribute significantly on minimizing the occurrence of flooding.Although the results of our study supported efforts to better understand the relationship of local planning capacity and the generation of peak runoff, some methodological limitations exist that need further investigation.
First, due to our stringent selection criteria and processes of the study area, a relatively small number of samples, 75 sub-basins, were chosen for this study.This is a threat to sound statistical conclusions.Further studies should increase the sample size by employing an alternative way of representing jurisdictions by watersheds, such as the weighting approach that Brody et al. [15] applied in their study.Because of the weak statistical power, the findings should be generalized to other areas with care, especially where natural and built environments have dissimilar patterns.
Second, this study equally weighted the indicators while evaluating the plan quality in order to reduce personal bias.Although typical plan evaluation studies have avoided prioritizing specific indicators for the purpose of minimizing vague judgments, future studies should assign differential weighting schemes by considering more opinions from planning practitioners, which may enable this evaluation protocol to be more useful in the practice [32,83].
Third, temporal limitations exist for runoff and precipitation data.The data in this study were examined on an annual basis.However, the amount of runoff and precipitation vary significantly by each month or season.For example, this study ruled out hydrological fluctuations that may be caused by snowmelt.Further studies should address temporal impacts of surface runoff, precipitation, and other natural environmental attributes.
Finally, to account for these temporal dimension issues and explain the causal relationship between independent and dependent variables, longitudinal analysis or panel analysis should be performed in future research rather than cross-sectional analysis if data are available.Particularly, panel analysis may better explain whether planning capacities have implementation effects in surface runoff by looking at the percent change of two distinct periods.

Figure 2 .
Figure 2. Selected sub-basins in the Chesapeake Bay watershed.

Figure 2 .
Figure 2. Selected sub-basins in the Chesapeake Bay watershed.
Notes: Number of observations is 75 for all variables.

Table 2 .
Indicators in sustainable stormwater management plan coding protocol.

Table 3 .
Factors influencing mean annual peak runoff.