Evaluating Watershed Vulnerability in Bernalillo County, New Mexico Using Expert Testimony, Fuzzy Analytic Hierarchy Process, and GIS

Geographic Information System (GIS) software was used to create a watershed vulnerability model for Bernalillo County, New Mexico. Watershed vulnerability was investigated as a function of soil erosion and infiltration criteria: precipitation, land slope, soil erodibility (K-factor), vegetation cover (NDVI), land use, drainage density, saturated hydraulic conductivity, and hydrologic soil group. Respective criteria weights were derived using a Fuzzy Analytic Hierarchy Process (FAHP) supported by expert opinion. A survey of 10 experts, representing New Mexico Institute of Mining and Technology (NMT), the New Mexico Bureau of Geology and Mineral Resources (NMBGMR), and the United States Geologic Survey (USGS), provided model input data for an integrated pair-wise comparison matrix for soil erosion and for infiltration. Individual criteria weights were determined by decomposing the respective fuzzy synthetic extent matrix using the centroid method. GIS layers were then combined based on criteria weights to produce maps of soil erosion potential and infiltration potential. A composite watershed vulnerability map was generated by equal weighting of each input map. Model results were categorized into five vulnerability categories: not vulnerable (N), slightly vulnerable (SV), moderately vulnerable (MV), highly vulnerable (HV), and extremely vulnerable (EV). The resulting FAHP/GIS model was used to generate a watershed vulnerability map of discrete areas in Bernalillo County, which may be vulnerable to stormwater run-off events and soil erosion. Such high volume run-off events can cause erosion damage to property and infrastructure. AlternaHow to cite this paper: Sadler, T., Richardson, C.P. and Miller, P. (2019) Evaluating Watershed Vulnerability in Bernalillo County, New Mexico Using Expert Testimony, Fuzzy Analytic Hierarchy Process, and GIS. Journal of Water Resource and Protection, 11, 866-885. https://doi.org/10.4236/jwarp.2019.117053 Received: June 12, 2019 Accepted: July 27, 2019 Published: July 30, 2019 Copyright © 2019 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/


Introduction
During precipitation events, water which is not absorbed into the ground flows over the ground surface eventually contributing to surface water features such as rivers and lakes. This phenomenon is called run-off and it is a major contributor of non-point source pollution, which can contaminate surface water features.
Run-off can also cause erosion damage to property and infrastructure and in extreme cases, risk to human life. From May 1996 to December 2017, the National Climatic Data Center (NCDC) reported a total of 1292 flood/flash flood events caused by stormwater run-off with 32 deaths, 15 injuries, over $115 million in property damage, and over $5 million in crop damage [1]. Urban runoff is a major contributor to elevated levels of suspended sediments and of harmful pollutants in streams, lakes, rivers, and oceans. Common runoff pollutants include nutrients, organic materials, pathogens, hydrocarbons, pesticides, metals, chlorides, grease, trash, and toxic substances [2]. These pollutants pose a serious threat to surface water quality and risk contamination of surrounding areas. Additionally, arid and semi-arid watersheds are especially vulnerable to run-off, due to infrequent but high-intensity episodic rainfall events which precipitate high volumes of water in a short period of time.
The management and control of stormwater run-off are of principle concern to community planners and local governments, as a matter of hazard mitigation as well as prevention of contamination to scarce local water resources. Geographic information systems (GIS) based approaches are being developed as tools to address this concern, focus efforts, and minimize issues associated with stormwater run-off. Understanding influencing watershed factors such as topography, precipitation, soil properties, vegetation, and drainage area, etc. is key to Journal of Water Resource and Protection the identification of problem areas for stormwater run-off. These areas, if geographically identified, may be evaluated for stormwater control measure (SCM) development. Structural and nonstructural SCMs are techniques, methods, activities, and practices that have the primary objective of controlling run-off volume and preventing or reducing the quantity of pollutants present in runoff water before reaching natural water systems (Richardson et al. 2018). Implementation of SCMs in strategic locations has the potential to prevent catastrophic run-off events and present a solution for decreasing surface waterway contamination.
Successful implementation of SCM strategies should be based on the unique geographic characteristics, local requirements of each project, and watershed-specific needs. Incorporating all of these potential decision making criteria for planning optimal SCM sites necessitates a system which allows for the comparison and prioritization of diverse criteria on a large scale. An extensively used technique for improving SCM efficiency is incorporating multicriteria decision support (MCDS) systems. MCDS involves the use of scientific and mathematical tools to evaluate complex trade-offs among different environmental and socioeconomic factors. Analytic hierarchy process (AHP) is one such commonly adopted MCDS technique. AHP is a measurement method to derive a scale of relative importance of alternatives or criteria from pairwise comparisons [3]. AHP input can represent actual quantitative measurements or designations of subjective judgment. A more evolving concept is the integration of MCDS and geographic information system (GIS) processes into SCM suitability studies, which allows for the pinpointing of optimal locations for SCM development.
An evaluation of watershed vulnerability in Bernalillo County, NM was performed by Amankwatia [4]. The research herein is based on this work, but uses fuzzy AHP (FAHP) to incorporate additional subjective judgment input regarding Bernalillo watershed processes. These subjective inputs were obtained from surveyed experts in the fields of geology and hydrology. Input from local experts was used to obtain a watershed vulnerability assessment of Bernalillo County comparable to the model produced in 2015.
One potential benefit of GIS-based MCDS models is that they can be easily implemented in regions where precise spatial data are not readily available without having to rely on expensive and time-consuming field surveys [2]. A second potential benefit of these methods is that they use local expert input to customize solutions to regionally specific watershed phenomena. Information from MCDS models can prove useful for overall watershed management and resource planning. This study combines specific soil erosion and infiltration criteria within a GIS-based fuzzy AHP framework to develop a watershed vulnerability map for Bernalillo County, New Mexico. The output could be beneficial for exploring optimization of SCM approaches to watershed management for targeted vulnerable areas within the county. Journal of Water Resource and Protection

Stormwater Control Measures in the Desert Southwest
Despite widespread support of stormwater control measures (SCM) and projects in the United States, there has been little representation of such approaches in the desert Southwest [5]. The lack of implementation of SCMs in the Southwest is closely tied to the perception that these regions experience little or no rainfall and, hence, the expense and burden of SCMs seem unjustifiable. However, Gautam et al. [5] argue that one important reason for establishing SCM in desert environs is that limited, but high-intensity, episodic rainfall does occur. Such events can be a major contributing factor to the elevated levels of harmful pollutants found in natural watercourses. Precipitation in New Mexico is highly variable with some months exhibiting little to no precipitation and summer monsoon storms bringing torrential downpours. These storms typical supply 50 percent of New Mexico's annual precipitation during the months of July, August, and early September [6].
In the semi-arid environment of New Mexico, stormwater is both an invaluable and limited resource. Limited availability of water resources and sustainability problems in this region actually intensify the need for the adoption of SCMs.
Gautam et al. [5] suggested that given the rapid urbanization and severe water availability issues in the desert Southwest, the best strategy for implementing SCMs in such areas should be based on sustainable land use, water conservation, and water reuse techniques. The authors also stated that one critical SCM for long-term sustainability is through low-impact development by reducing impervious areas and limiting directly connected impervious surfaces.

MCDS Approach to Watershed Vulnerability and SCM Using AHP and FAHP
The growing interest in the use of analytical hierarchy process (AHP) in SCM suitability studies is associated with the ability of AHP to deal with imprecision in environmental data. As reported by Young et al. [7], use of AHP supported with input from GIS can provide an objective, mathematically-based alternative to the existing subjective SCM approaches. Fuzzy AHP (FAHP) provides the same benefits and incorporates the ability to crowdsource decision prioritization. AHP has been successfully applied in assessing regional scale urban flood hazard [8] [9]. Research has also been carried out on the effective use of AHP and GIS to manage surface-water and groundwater resources. For example, Agarwal et al. [10] suggested that AHP can be effective to evaluate groundwater potential risks. A more recent attempt suggested that AHP techniques provide an effective framework for reassessing vulnerability of groundwater resources and aquifer systems [11]. Additional AHP studies have focused on soil erosion risk [12] and have shown it to be an effective technique of predicting soil erosion risks at precisions equivalent to field observations [13]. Fuzzy AHP adds an additional element of flexibility in data integration. The objective of this study was to implement a composite FAHP model based on soil erosion and infiltration T. Sadler et al. Journal of Water Resource and Protection with prioritization of decision criteria determined by a fuzzy analysis of expert knowledge of the study region. The purpose of using a weighted decision system is to provide a geographic framework for watershed vulnerability assessment, which is uniquely suited to the study watershed. This technique results in a vulnerability map, which can inform SCM implementation predicated by expert knowledge of the interest area.

Factors Influencing Watershed Vulnerability
Watershed vulnerability, defined as high soil erosion and low infiltration potential, is a complex process to model. Several watershed factors, such as slope gradient, soil texture, soil erodibility, soil permeability, land use and land cover, vegetation cover, and drainage density, influence the development of the watershed vulnerability model. The characteristics and effect of these factors on soil erosion and infiltration are outlined in this section.

Slope Gradient
Slope gradient (S), usually expressed as percent slope or degree slope, is used to represent the degree of variability in elevation and is found to strongly influence

Soil Texture
Soil texture, defined as ratios of sand and clay, is known to have a strong influ- Hydraulic soil class dataset was downloaded from the USDA database and accessed using Soil Data Viewer.

Soil Erodibility
The soil erodibility (K) factor is a quantitative soil property that plays a major role in determining the extent to which soil particles are eroded or transported by runoff water. Parameters such as soil moisture, soil permeability, and organic matter content influence the degree of soil erodibility and stability. Soils with Journal of Water Resource and Protection low K values are highly stable compared with those with higher K values.
K-factor values were obtained through SSURGO using Soil Data Viewer.

Vegetation Cover
The normalized difference vegetation index (NDVI) is a widely used standardization parameter to represent vegetation cover. NDVI values typically range from 1.0 to −1.0 with non-vegetated areas such as streams and rivers producing small or slightly negative values, and vegetated areas producing values ranging from 0.4 to 1.0 [14]. Landsat 8 data was downloaded for the area covering Bernalillo County using Earth Explorer provided by the United States Geologic Survey (USGS). Data was collected from cloud-free fly-overs in June 2018. June represents the driest month of the year in New Mexico and conditions which would be present at the start of the annual monsoon season [6]. Landsat 8 data was converted to reflectance, corrected for sun angle, and NDVI was calculated from the red and near infrared data bands.

Soil Permeability
Soil permeability, or more correctly saturated hydraulic conductivity (Ksat), is an important physical soil property that refers to the ability of soils to conduct water. High Ksat values represent the portions of watersheds that are highly drained and are suitable for infiltration practices. Conversely, lower Ksat values correspond to watershed regions that experience high runoff and are prone to soil erosion. Ksat values were obtained through SSURGO using Soil Data Viewer.

Land Use
Land use for purposes herein refers to the use of land for human activities. Land use directly translates to the amount of water penetration into soil layers or being produced as runoff water. For example, rangeland or grassland will have different infiltration or erosion potential than land which is used for industrial purposes. Land use data was downloaded from the USDA Geospatial Data Gateway and classified into risk categories.

Precipitation
Area precipitation characteristics (type, duration, intensity, and amount) are highly important in drainage projects because of their direct influence on the volume of generated runoff. Precipitation characteristics play a significant role in soil erosion. In civil and environmental engineering, storm classification is used to determine appropriate stormwater control measures. Storm intensity and duration are important factors in determining the volume of water precipitated by a storm and the erosional or run-off risk that the storm will pose to a watershed. In semiarid watersheds, high-intensity precipitation storms are believed to be the main component for soil erosion. Precipitation frequency (PF) estimates are used as design criteria for a wide variety of civil and environmental T. Sadler et al. Journal of Water Resource and Protection engineering applications, such as the design of hydraulic structures such as culverts, roadway drainages, and bridges. PF estimate data were obtained from NOAA's Hydrometeorological Design Studies Center-Precipitation Frequency Data Server (PFDS). This dataset consists of spatially distributed point PF estimates formatted as a georaster of the Southwestern United States. A high intensity, high volume storm is likely to pose a higher risk for erosion and storm water run-off. For this project, a precipitation depth for a 10-year, 6-hour storm was used to approximate precipitation for a severe monsoon storm, but alternate storm design frequencies could easily be substituted in the FAHP-GIS method.

Drainage Density
Drainage density is a quantitative indicator that represents the ratio of the total length of streams to the total watershed area. Low drainage density values are favorable for infiltration, whereas larger values indicate high runoff potential. High drainage density values are an indication of high soil erosion vulnerability, whereas lower values signify minor risk. The linear drainage file was downloaded from USDA Geospatial Data Gateway and drainage density was determined using the ArcMap line density function.

Site Description
This study was conducted using maps and digital data from Bernalillo County, New Mexico. The county is located in North Central New Mexico, centered at approximately longitude 106˚41'56.40''W and latitude 35˚3'14.4''N. Bernalillo County has a total area of approximately 3023 km 2 (1167 mi 2 ). Bernalillo County is the most highly populated county in the State of New Mexico, and is home to the City of Albuquerque, its county seat. The region is covered by shrubs and scrubs (71%), deciduous, mixed, and evergreen forests (18%), associated urban (9%), and open water (1%). The two dominant soil orders are Entisols (47%) and Aridisols (32%) [2]. Albuquerque's urban sprawl consisting of residential, commercial, military, and industrial land perpetuates a high risk of surface run-off events. Urban surface conditions result in contamination of surface water run-off, which poses risk nearby waterways and nearby agricultural development. Limited water availability, unique soil characteristics, dominant shrub-land vegetation, and intricate climate, coupled with a forecasted 35% population growth over the next 20 years [15], makes Bernalillo County a suitable case study area for the current research.

Modeling Approach
The overall modeling approach used for this project is based on Fuzzy AHP multicriteria decision making utilizing expert opinions collected by survey ( Figure 1). The principal benefit of this method as an approach for evaluating watershed vulnerability is its ability to objectively and simultaneously consider an unlimited number of relative criteria. For this project, two submodels are used to evaluate watershed vulnerability. Journal of Water Resource and Protection One submodel assesses soil erosion potential and the other assesses infiltration potential; each model incorporates multiple criteria. Coupled, these processes influence watershed vulnerability for stormwater run-off. The premise herein is that areas with high soil erosion potential and low infiltration potential are most vulnerable.

Soil Erosion Submodel
Soil erosion poses an imminent threat to infrastructure and property. As recently as 2018, heavy monsoon rains in the area flooded homes, eroded dirt roads, blocked access to homes and businesses, and left thick blankets of sediment on agricultural crops. The City of Corrales, New Mexico in Sandoval County adjacent to Bernalillo County saw up to three feet of silt and sand deposited in areas and extreme erosion in areas during the July 2013 floods [1]. The soil erosion submodel combines precipitation, soil erodibility, land slope, land use, and vegetation cover.

Infiltration Submodel
Infiltration is the entry of water into the soil surface and its subsequent vertical motion through the soil profile. Areas which experience low soil infiltration will be more highly vulnerable for surface run-off events. The soil infiltration submodel uses drainage density, hydrologic soil classification (HSC), land slope, land use, and saturated hydraulic conductivity as criteria to establish run-off potential.

ArcGIS
ArcGIS is geographic information software system, which provides an extensive toolset to manage multiple layers of spatially-distributed geographic data.
Modeling approaches based on GIS have been shown to be successful in assessing regional-scale soil erosion risks [16]. GIS-based models have also been extensively used in evaluating groundwater vulnerability in order to aid policy formulation. ArcGIS provides an infrastructure for making maps and working with large data sets of geographic and remotely sensed data. In this study, ArcMap was used for a range of geoprocessing tasks including spatial referencing, resampling, mosaicking, clipping, reformation, and classifying geographic data. ArcGIS's raster calculator function was used to combine thematic data layers, creating erosion and infiltration models. Data analysis and decision support tools were critical in creating custom watershed vulnerability maps.

Soil Data Viewer
The NRCS Soil Data Viewer is an ArcGIS extension for analyzing soil data and building soil-based thematic maps. The program enables map makers to interpret soil attributes and properties without having to deal with the complexity associated with querying and processing soil databases. Soil Data Viewer was used to create and export thematic layers such as soil texture, hydrologic soil groups, soil depth, and saturated hydraulic conductivity maps for the study area.

Data Processing
All layers were projected in ArcGIS to North American Datum (NAD) 1983 UTM zone 13 projection prior to assimilation into the model. A grid size of 30 × 30 (m) pixels was chosen for all spatial layers to cover the study area extent.
All feature datasets were converted to raster format and all raster datasets were spatially resampled to 30 (m) resolution using the Data Management toolbox of ArcGIS. The Data Management toolbox was then used to clip the input thematic layers to conform to the spatial extent of the study area. In order to produce the soil data maps (including soil texture K-factor and Ksat), SSURGO data were loaded into Soil Data Viewer and individual layers were selected for export.

Expert Survey
A survey conducted of experts in the fields of geology and hydrology was con- slightly more, moderately more, strongly more, or absolutely more important than criterion B. Or conversely, participants indicated whether criterion B was equally as, slightly more, moderately more, strongly more, or absolutely more important than criterion A. In this distributed survey (soil erosion and infiltration, respectively), each of the five influencing criteria is compared against each other criteria for the same phenomenon. Table 1 shows the relative scale of importance used for the respective pair-wise comparisons.

Initial Comparison Matrix
Results from the criteria evaluation expert survey (10 participants) were used to construct an initial comparison matrix for each response as shown in Figure 2 and Applying a triangular distribution to the numbers results in a minimum, mean, and maximum value for each importance comparison. The composite initial matrix is transformed into an integrated fuzzified pair-wise comparison matrix (FPCM) using a geometric mean of fuzzy triangular functions [17]. An example of the integrated FPCM for soil erosion is provided as

Consistency Evaluation
The principal behind using FAHP as a multicriteria decision making process (MCDS) relies on the consistency of input for the pairwise comparison matrix.
In order to ensure that AHP judgments are accurate and reliable, a consistency ratio (CR) for each submodel was calculated according to method proposed by [18] ( ) ( )  where n is the number of criteria in the comparison, λ max is the maximum Eigenvalue, and RI is a random consistency index based on the number of criteria.
The more consistent a comparison matrix is, the closer the value of computed maximum Eigenvalue (λ max ) is to the number of criteria. Ideally, the consistency ratio should be less than 10% for each submodel [19]. For both erosion and infiltration the consistency ratio was found to be less than 0.02% or 2%. Because there are five comparison criteria for each submodel, the maximum Eigenvalue should approach a value of 5 for both. The calculated consistency ratios and maximum Eigenvalues for the erosion and infiltration sub models are CR erosion = 0.0161; λ max erosion = 5.07 and CR infiltration = 0.0146; λ max infiltration = 5.06, respectively.

Defuzzification of Fuzzy Pair-Wise Comparison Matrix
Determining the crisp weight of each criteria from the integrated fuzzy comparison matrix uses the fuzzy extent analysis (FEA) method by Chang [20] with the Wang et al. [21] modification of the normalization process applied to calculate the values in the n x n fuzzy synthetic extent matrix (FSEM). The centroid method was then used to calculate crisp criteria weights for each of the influencing criteria. The centroid method calculates the average for row i in the FSEM as (l i + m i + u i )/3, which are summed to give a total for n rows. Respective crisp criteria weights (w i ) are determined by dividing row i centroid average by the total.

Vulnerability Rating of Fuzzy AHP Model
A five-category rating system was used to classify vulnerability for this research.
Each category is used to define a condition of vulnerability from "not vulnerable" to "extremely vulnerable". Vulnerability categorizations are as follows: 1)

Vulnerability Assessment
The next step is to assign vulnerability ratings for each criterion. Ranges were established to classify each criterion into risk categories for vulnerability. Rating categories were modeled after the works of Richardson and Amankwatia [2].
Categories were defined to reflect the level of vulnerability a criterion poses for stormwater run-off. For example, a slope of 0 -2 degrees is considered not vulnerable, 2 -5 degrees is slightly vulnerable, 5 -10 degrees is moderately vulnerable, 10 -18 degrees is highly vulnerable. and any slope greater than 18 degrees is extremely vulnerable. Each class is attributed to a value 1 -5. All of the weight values for each submodel are scaled to add up to 1, such that after all the criteria are weighted and combined, the maximum vulnerability value remains on a scale from 1 to 5. Also note that land cover classifications "herbaceous", "hay and pasture", and "cultivated crops" were reclassified for this research as "cultivated". Criteria vulnerability ratings and the crisp weights derived from fuzzy AHP are presented in Table 3 and Table 4.

Evaluation of Watershed Vulnerability
The weighted overlay of thematic map layers created a map of Bernalillo County, which indicates erosion vulnerability and infiltration vulnerability due to a storm event. The map of soil erosion potential for Bernalillo County is shown in Figure 4. Infiltration potential for Bernalillo County is shown in Figure 5. These two maps were combined to produce a map of overall watershed vulnerability for run-off related events. The scaled watershed vulnerability of Bernalillo County is given in Figure 6. Warm or red and orange-colored areas indicate regions of high vulnerability and cool-colored areas or blue and green areas indicate areas of lower vulnerability. For this study, the relationship of erosion and infiltration processes is combined with equal weighting to produce a map of combined vulnerability. Weighting factors for each influencing parameter included in the study of erosion potential and infiltration potential are shown in Table 3 and Table 4, respectively. NDVI is a dominant criterion influencing soil erosion, followed by precipitation. Saturated hydraulic conductivity and soil texture are equally dominant as influencing infiltration.
Areas of high infiltration vulnerability (or low infiltration) and areas of high erosion vulnerability designate areas are highly or extremely vulnerable for the occurrence of run-off related events. These high vulnerability areas, shown as red in Figure 6, are areas were stormwater control measures (SCM) can be focused for the protection of property and infrastructure as well as prevention of surface waterway contamination. Journal of Water Resource and Protection     Figure 8. The mean vulnerability value of the model was determined to be 3.07 (on a scale of 1 to 5) with a standard deviation of 0.28. The majority of the pixels lie within two standard deviations of the mean. The distribution indicates that the "not vulnerable" and "extremely vulnerable" categories do not represent a significant portion of total land area within Bernalillo County.

Advantages of Fuzzy AHP Model Approach
Multicriteria decision support approaches such as FAHP provide a relatively easy and customizable way of determining watershed vulnerability distributions on a regional scale. These methods can provide several strategic benefits to decision makers including the ability to incorporate a variety of quantitative data, the experience of local experts, and project specific objectives. An advantage of FAHP is the method's ability to deal with data uncertainty and lack of sampled erosion, infiltration, or other watershed data. However, developing pairwise comparisons for AHP hinge on subjective judgment and can create a great deal of variability in model decision making priorities. FAHP provides the opportunity to combine the judgement of multiple information sources, to provide additional stability to the normalized priority vector values. A sensitivity analysis can be performed wherein the input data are slightly modified in order to observe Journal of Water Resource and Protection  the impact on the weighted priority results. If the rank order does not change, the results are said to be robust; otherwise, they are sensitive.
The FAHP approach presented herein offers several benefits in developing an understanding of regional-scale watershed vulnerability. The integration of FAHP and GIS techniques into watershed vulnerability studies for SCM suitability can be easily implemented in regions where precise spatial data are not readily available without having to rely on expensive and time-consuming field surveys. Information from MCDS models can prove useful for large scale watershed management and resource planning, narrowing down areas which can later be tested for SCM suitability.

Conclusions
This study combines specific soil erosion and infiltration criteria within a GIS-based fuzzy AHP (FAHP) framework to develop a watershed vulnerability map for Bernalillo County, New Mexico. Quantifying a regional-scale risk index map can provide spatially relevant information for engineering management and decision making. In view of this, a 30 × 30 (m) watershed vulnerability risk map Journal of Water Resource and Protection was generated for the Bernalillo County, New Mexico using an FAHP methodology based on judgements of experts on local watershed processes. The approach incorporated drainage density, slope, precipitation, vegetative cover, land use, soil erodibility, soil texture, and soil permeability as the primary factors contributing to watershed vulnerability. The final map was produced by merging two submodels representing soil erosion and infiltration potential using the ArcGIS software. Analysis of the vulnerability map identifies that the most problematic areas are present in the Eastern and Northwestern portions of the study area. However, based on the analysis less than 1% of the total area lies within the lowest and highest vulnerability categories with the majority of land area centered around moderate vulnerability.
Results were compared with those using AHP derived weights from the analysis is required to determine model sensitivity to each influencing criterion. Fuzzy AHP is, therefore, an easily applied numerical method for evaluating criteria weights for MDCS systems. It allows for the quantitative evaluation of data from multiple, sometimes incomparable, sources by creating a system of numeric weights for a quantitative comparison and data incorporation. The method used in this process, couples two key hydrologic phenomena to assess watershed vulnerability (erosion and infiltration). The purpose of the application of this process is to identify vulnerable areas in need of implementation of stormwater management measures. FAHP relies on the experience of local experts and allows for the inclusion of additional evaluation datasets (criteria) customizable to an area's regionally specific watershed phenomena. This method produces site-specific vulnerability estimations that can be used for local discussions about individual risk. Information from the watershed vulnerability model can be useful for managing new structural or nonstructural SCM projects or improving existing ones to mitigate impact to stormwater contamination and run-off volume for identified vulnerable areas. Better prediction of problems allows for better planning and management of hazards, which contribute to watershed vulnerability.

Recommendation for Future Work
Future work should include field investigations of regional susceptibility to stormwater run-off and erosion in order to validate the usefulness of the model and for use in stormwater and watershed management. In addition to photographic documentation from various field sites, and quantitative erosion and infiltration data points could be used to adjust the vulnerability model and test for accuracy using a confusion matrix approach [22]. Journal of Water Resource and Protection Impacts of variation in spatial resolution of datasets should also be explored.
While most datasets including SPADTM, SSURGO, Landsat 8, and USDA landuse data had a resolution of 30 (m) or less, the climate precipitation data was provided in an 800 (m) raster. This results in reduced spatial variation in the precipitation data when it is incorporated into the model and may impact interpreted vulnerability by area. Amankwatia [4] used a 100 yr 24 hr precipitation depth criterion versus a 10 yr 6 hr precipitation depth criterion herein. NDVI data appropriate to the critical premonsoonal period should be extracted for multiple years and averaged to better address soil erosion potential. Amankwatia