Validation of a Locally Revised Topographic Index in Central New Jersey, Usa

Saturation excess is a major runoff process in humid regions such as the US Northeast. Topographic index (TI) is used to simulate the pattern of runoff-contributing areas following a saturation excess runoff process. Although TI is useful to delineate saturated areas, i.e., hydrologically sensitive areas, for taking spatially distinctive actions in watersheds for improving water quality, local resource management practitioners often question its applicability to local conditions. This study introduces two methods to validate a locally revised TI in humid central New Jersey, USA. The revised TI uses soil moisture deficit instead of water table depth as the state variable in simulating the saturation excess runoff process. First, the calculated TI values were compared to the soil moisture measurements sampled in two sites in Tewksbury Township in the region to evaluate their correlations. Second, a watershed model Variable Source Load Function (VSLF) that incorporates TI was applied to the Neshanic River Watershed in the region and compared to another watershed model Soil and Water Assessment Tool (SWAT) to evaluate its capability in predicting the streamflow and its runoff and baseflow components. The positive correlations between soil moisture measurements and TI suggested TI is a good indicator of runoff-generating potential. VSLF achieves a modeling efficiency comparable to SWAT in simulating watershed hydrology. Such validation gives practitioners confidence to incorporate TI pattern into watershed management practices for improving their efficiency. The results are applicable to shallow, interflow-driven watersheds in humid regions.


Introduction
Saturation excess is a major runoff process in humid regions such as the US Northeast and can be better captured by a physical hydrology concept called Variable Source Area (VSA) hydrology [1][2][3].The VSA hydrology applies a saturation-excess hydrological process to explain that runoff-carrying pollutants contribute to the stream hydrograph are primarily generated in relatively small but predictable hydrologically sensitive areas in a watershed where soils are saturated [4][5][6].The VSA hydrology presents interesting insights regarding the complex hydrological connections between streams and their upland contributing areas and offers practical guidance on managing the upland contributing areas to achieve desired water quality in streams.The VSA hydrology has been closely linked to the delineation and management of critical source areas (CSAs) for achieving long-term water quality goals and the sustainability of managing water resources more efficiently and effectively [7][8][9][10][11][12][13].
The topographic index (TI) was first introduced as a basis for a rainfall-runoff model, the TOPography-based hydrological MODEL (TOPMODEL), which simulates the pattern of surface runoff-contributing areas following the VSA hydrology [13].The index measures the propensity of a given point in a watershed to become saturated and act as a source area for surface runoff [14].The index is often called a wetness index and is used to predict soil moisture distribution in the landscape [15,16].Because of its simplicity and ability to visualize the spatial predictions, TI has been widely adopted to simulate hydrologic processes with various modifications and extensions that reflect different modeling assumptions and take account of additional natural resource conditions [17][18][19].TI and its variants have been incorporated into rainfall-runoff models to predict streamflow [20][21][22][23].
TI and its variants had been subject to intensive investigation on its ability to predict the spatial pattern of runoff-contributing areas in hydrology literature.The predicted spatial patterns were compared to the measured soil moisture pattern [24][25][26], the shallow groundwater table [27,28], and the observed wetland occurrence [29,30].However, the validation results were mixed.While reasonable correlation between observed soil moisture patterns and TI was found [31,32], many found a poor correlation between predicted and observed spatial patterns, especially the point-to-point correlation [33][34][35].The poor correlation was not entirely unexpected because of measurement errors and varying local resource conditions [35].
One TI variant was to re-conceptualize TI using the soil moisture deficit instead of the water table depth as the state variable to make it applicable to shallow, interflow-driven watersheds in the US Northeast [36].The revised TI has been extensively evaluated in terms of its ability to predict the spatial pattern of runoff-contributing areas, and to predict streamflow in the New York Catskill Mountains.The patterns of runoff potential as predicted by the revised TI were compared to the field soil monitored moisture pattern and to the runoff pattern predicted by a fully distributed Soil Moisture Routing (SMR) model in the Town Brook Watershed in the same region [11,37].SMR is a mechanistic physical-based, spatially distributed model that captures the VSA hydrology of humid, temperate regions such as the US Northeast that have highly permeable, sloping soils over a restricting layer [38].They concluded the revised TI was a good indicator of runoff risk, especially in US Northeast.A normalized difference water index derived from medium resolution Landsat 7 Enhanced Thematic Mapper Plus (ETM+) images using an unsupervised classification method was used to delineate the saturated areas, and the delineated pattern of saturated areas matched the predictions made using the revised TI in the Town Brook watershed well [39].A variable source loading function (VSLF) model was developed by incorporating the revised TI to modify the Soil Conservation Service curve number (SCS-CN) rainfall-runoff method in the popular General Watershed Loading Function (GWLF) model [40,41].The VSLF model provides a much more valid spatial distribution of soil moisture and runoff generation as compared to the GWLF and other current SCS-CN-based water quality models.Consequently, using such a model could improve the simulation accuracy of pollutants transported by runoff [23].
The objective of this study is to validate this locally revised TI in central New Jersey, which is located in the humid US Northeast.Two methods were used to evaluate how well the revised TI would perform in terms of its ability to predict VSA patterns and watershed hydrology in central New Jersey.First, the soil moistures were measured using a time domain reflectometry (TDR) soil moisture meter in selected fields and the TDR measurements were compared to the TI values to evaluate their correlations.Second, the VSLF model was applied to a suburban watershed with mixed land uses located in central New Jersey and the modeling results will be compared to a distributed biophysical model called the Soil and Water Assessment Tool (SWAT), which was applied to the same watershed to evaluate its capability in predicting streamflow and its runoff and baseflow components.Located at the outskirts of the New York metropolitan area, central New Jersey has been experiencing rapid urban development as a part of suburbanization and exurbanization processes, which creates substantial stresses to natural resources.Local communities have been active in protecting their natural resources for achieving the long-term sustainability of water resources in the region.The TI validation results would help the local communities to adopt spatially targeting strategies to direct their protection efforts to CSAs and therefore achieve a greater efficiency in protection [12,13].

Study Area
The study was conducted in Hunterdon County, which is located in central New Jersey (Figure 1).Soil moisture measurement was conducted on two sites in Tewksbury Township to validate the derived TI.The Cold Brook Preserve is a 116 hectare (ha) preserved parkland managed by the Hunterdon County Parks and Recreation Division, much of which was contracted to a local farm for active cropping activities.The Christy Hoffman Farm Park is a 69 ha community park managed by Tewksbury Township for recreation.Both sites consist of gently rolling hills, have diverse topographic variation, and are easily accessible from county roads.They were selected because of those characteristics and the close contact between the research team and the public agencies that manage these sites.
Watershed modeling was conducted in the 7897 ha Neshanic River Watershed, a suburban watershed with mixed land uses located in Hunterdon County.The watershed is comprised of 36% agriculture, 31% urban, 21% forest, and 10% wetlands, with the remaining lands both barren and in water.The climate in the region is humid with typical hot and humid summers and cold winters.The average high temperature in summer is 27-30 °C and the average low temperature in winter is −7 to −5 °C.The mean annual precipitation in the watershed is 1218 mm.Snowfall in the winter season is 120-770 mm, but varies yearly.Urban development and agricultural operation have degraded water quality in the Neshanic River, leading to impairment from dissolved oxygen, phosphorus, total suspended solids, and nitrogen [42].Various best management practices were taken to improve water quality in the watershed [43].

Topographic Index
The most common form of TI (λ) is the wetness index: where α is the upslope contributing area per unit contour length in meters to a given point in a watershed and β is the local surface slope angle in radians.The revised TI that predicts the spatial distribution of runoff-contributing areas is approximated as follows [36]: where D is the local soil depth in meters to the fragipan, bedrock, or other type of restrictive layer; Ks is the mean saturated hydraulic conductivity of the soil profile in meters per day above the restrictive layer; and the other variables are defined as before.Equation (2) has two components: the wetness index, ln( ( ) ), and the soil water storage index,

ln(
).In addition to the topography, the revised TI considers the impacts of the soil water storage capacity (KsD) above the restrictive layer.In general, the deeper the topsoil (D) and the higher the saturated hydraulic conductivity (Ks), the lower the λ value, which implies a lower likelihood of generating surface runoff.In reality, there may be several topsoil layers with different Ks values above a restrictive layer or bedrock.In this case, Ks is defined as: where d is total thickness of the soil above the restrictive layer; di is the thickness of layer i; ki is the saturated hydraulic conductivity of layer i [44].
The revised TI presented in Equation ( 2) was derived using two spatial datasets: a 10 m digital elevation model (DEM) obtained from the New Jersey Department of Environmental Protection (NJDEP) and a Soil Survey Geographic (SSURGO) soil database obtained from the US Department of Agriculture Natural Resources Conservation Service.DEM was first processed using the open source ArcGIS extension Terrain Analysis Using Digital Elevation Models (TauDEM) Version 3.1 [45] and ArcGIS Raster Calculator to obtain a wetness index grid.The soil depth (D) and the mean saturated hydraulic conductivity (Ks) were extracted for each soil type from the SSURGO soil database.The product (Ks*D) for each soil type was calculated and linked to the soil spatial distribution GIS layer.The soil layer was converted into a soil water storage index raster layer based on the value of ln( ).The converted raster layer had the same resolution as the wetness index raster layer and was also in alignment with the wetness index raster layer.The ArcGIS Raster Calculator was used to manipulate the raster layers for the wetness index derived from the DEM and the soil water storage index from the SSURGO soil data and to calculate the TI based on Equation (2).

Soil Moisture Measurement
The FieldScout TDR 300 soil moisture meter with a 12 cm rod was used to measure soil moisture at the depth of 0-12 cm in the Cold Brook Reserve and Christy Hoffman Farm Park.The TDR meter converts a measured electrical signal into soil moisture content using an equation that was validated over a wide range of mineral soils.The calibration is not required for general use, and is only recommended for maximum accuracy [46].The TDR meter has two volumetric water content (VWC) modes for measuring soil moisture: one for standard soils and the other for highly clay soils.Based on the soil conditions in the fields, the mode for standard soils was used in this study and no calibration was conducted.Locations where soil moisture was measured were recorded using a Global Positioning System (GPS) system, Trimble GeoXH handheld.The GeoXH handheld provides real-time subfoot (<30 cm) accuracy with the internal antenna, or decimeter (10 cm/4 inch) accuracy after post-processing.NJDEP GPS Data Collection Standards for GIS Data Development were followed when collecting the location data [47].
The Cold Brook Reserve was visited three times on 8 June 2009, 1 December 2009, and 30 April 2010, respectively.Christy Hoffman Farm Park was visited two times on 1 December 2009 and 17 December 2009, respectively.Locations for taking soil moisture measurements using TDR were determined on site to ensure these locations were spread across variety of topographic changes such as flat and sloping areas.Different routes in these fields were taken for those visits.The adjacent sampling locations were roughly 10 meters apart, considering the restriction imposed by the DEM resolution.Once a location was selected, three soil moisture measurements were taken at three points around the location using the TDR meter.The three points were generally one meter apart, forming a triangle shape.The average of the three soil moisture measurements was used for further analysis.Table 1 summarizes the soil moisture measurement activities at the two sites.The recorded sampling locations and the measured soil moisture data were imported into ArcGIS for data processing and interpolation.The value of the derived TI for each sampling location was extracted from the calculated TI raster layer.The spatial patterns of the measured soil moistures were analyzed using various interpolation techniques in the Spatial Analyst to virtually evaluate whether they match the spatial patterns of the TI values.Since the sample points were selected on site as the TDR soil moisture was measured, there was no control over the corresponding TI values of the selected locations, which resulted in a uneven distribution of sample locations corresponding to different TI values.To offset the impacts of such unevenness in data distribution, the statistical analysis of the relationship between the TDR measurements and TI values was based on binned data.The pairs of data at all sampling locations were first sorted by TI values.Based on the ranges of TI values, the data was separated into several TI value bins with incremental intervals of TI values.The number of bins and the corresponding ranges of the TI values vary by each sampling event.The average TDR measurement was calculated for all sampling locations falling within a specific bin.The regression analysis and plotting were based on the lower bound of the TI value in each bin and the corresponding average TDR measurement.

VSLF Modeling
The VSLF model was applied to the 6493 ha upper portion of the Neshanic River Watershed which drains in to the U.S. Geological Survey (USGS) Reaville Gage Station at the intersection of Reaville Road and the Neshanic River as shown in Figure 2. The standard SCS-CN method used in the GWLF model was reconceptualized to estimate the runoff from VSAs in the VSLF model [23].The VSLF model simulates the watershed runoff response to rainfall using the standard SCS-CN method, and spatially distributes the runoff response according to a soil saturation index.While the runoff estimation from the impervious surfaces was kept the same as in GWLF, the runoff estimation in the pervious portion of the watershed was revised to capture the saturation excess runoff process.The watershed was divided into 10 saturation classes with equal areas based on the TI values, as shown in Figure 2. The watershed runoff depth Q is the aerially weighted sum of runoff depths qi for all discrete contributing areas: where i is the index of saturation class, ΔAs,i is the fraction of the watershed in the saturation class i, and qi is the per unit runoff depth for saturation class i predicted by the traditional SCS-CN rainfall-runoff method based on land use and soil conditions within the saturation class.
The 2002 land use/cover data developed from the aerial photographs taken in the spring of 2001 by NJDEP was overlapped with the derived saturation class layer to identify the land use distribution in each saturation class.The NJDEP land use/cover data used a modified Anderson classification system [48].
The land uses were classified into six broad land use categories including agriculture, barren, forest, urban, water, and wetlands, and 50 subcategories using four digital land use classification codes.The NJDEP land use/cover data contained the impervious surface rate for each land use polygon.The impervious surface rate was used to separate the total area of each land use category in each saturation class into two parts: pervious and impervious areas.The NJDEP four digital land use classification codes for the pervious areas in each saturation class were reclassified into the following land use categories required by the VSLF model: forest deciduous; coniferous forest; mixed forest; brushland; cropland and pasture; residential; commercial and industrial, roads; wetlands; and water.The impervious areas in all saturation classes were lumped into the following land use categories for the watershed: barnyard, residential impervious, commercial and industrial impervious, and roads.For example, the impervious areas from all agricultural lands and residential land uses were lumped into barnyard and residential impervious, respectively.The NJDEP land use classification did not distinguish the specific uses of agricultural lands.In the application, the total agricultural lands in each saturation class were further divided into pasture, permanent hay, and cropland based on the following percentages: 30%, 20%, and 50%, respectively.Those percentages were based on the field observation data obtained through two rounds of agricultural land use inventories throughout the watershed during the period 2007-2008.
The precipitation, maximum and minimum temperature during 1960-2008, solar radiation, and relative humidity during 1960-2008 in the Flemington weather station, which was just outside of the watershed, were collected for the modeling.Daily streamflow was obtained from the Reaville Gage Station, observed over the years 1960-2008, and was compiled and used as a model input.
The VSLF model was calibrated for the period 1997-2002.A 20-step program to calibrate each of the calibration parameters was used to calibrate the model.In each step the value of a parameter was adjusted until an optimization goal was reached.The optimization goal is either to minimize the bias (absolute value of difference between mean observed and mean simulated), or to maximize the Nash-Sutcliffe efficiency (NSE) [49] in comparing the simulated to the observed output variable affected.The calibrated model was then used to validate its capacity to simulate the streamflow during the period 2003-2008.To calibrate and validate the VSLF model, the observed and VSLF-simulated streamflow was further separated into two components: surface runoff and baseflow.Specifically, runoff was calculated by the baseflow separation method using beta = 0.925 [50].Eventswitch divides the time series into successive storm event periods using a minimum hydrograph rise threshold of 0.001 cm [4].The eventswitch time series was used in the model for defining events.The means and standard deviations of the VSLF-simulated streamflow and its two components were compared to the observed values.The goodness of fit, i.e., agreement between simulated and observed values, was assessed using NSE on the [−∞, 1] interval.The closer the coefficient is to 1, the better the fit.

SWAT Modeling
A similar watershed modeling effort was conducted in the Neshanic River Watershed using SWAT.SWAT is one of many distributed watershed models that has been widely used to understand the hydrological cycle and processes, simulate the hydrological and water quality impacts of land use and management practices, and evaluate alternative management strategies to improve water quality and ecosystem functions in watersheds [51,52].Unlike the VSLF model that distinguishes a watershed into 10 saturation classes to evaluate its hydrological responses under different types of land uses, SWAT delineates a watershed into smaller subbasins based on hydrological connections and further into hydrological response units (HRUs) base on land use, topography, and soil conditions [53,54].HRUs are unique combinations of land use, soil, and slope, which allow the model to reflect differences in evapotranspiration and other hydrologic conditions across subbasins and the watershed.
Based on DEM and the stream network, the watershed was delineated into 25 subbasins as shown in The streamflow observed at the Reaville Gage Station during 1960-2008 was similarly processed into surface runoff and baseflow as discussed before to calibrate and validate the SWAT model.NSE was used to evaluate the goodness of fit between the SWAT-simulated and the observed streamflow and their two components.The SWAT model was calibrated for the period 1997-2002 and validated for the period 2003-2008 to be compared to the VSLF model.

Comparison between DTR Measurements and TI Values
The interpolated soil moisture and the predicted TI patterns did not show good correlation, which further confirmed the difficulty of point-to-point comparison [35].Site-specific soil moisture is affected by many different factors such as vegetation, small variations in topography, and soil compact.The statistical approach based on binned data was used to eliminate the impacts of site-specific factors to assess the relationship between TI and TDR measurements.
Figure 3 presents the soil moisture sampling locations and the relationship between TDR measurements versus TI values at Cold Brook Farm during three field visits on 8 June 2009, 1 December 2009, and 30 April 2010.Three graphs showed a positive relationship between TDR measurements and TI values, i.e., a higher TI value corresponded to a higher TDR measurement.The relationship was especially strong when the sampling took place across a larger variation of topographic change indicated by a wider range of TI values and a wetter condition indicated by higher TDR measurements.The R 2 values were 0.94 and 0.88 on 8 June 2009 and 30 April 2010, respectively.On the other hand, the R 2 value was only 0.45 despite a positive relationship between the TDR measurement and TI value on 1 December 2009 when the sampling took place in a gradient with a small topographic variation in a dry weather condition.At Christy Hoffman Farm Park, as shown in Figure 4, the TDR measurements were also perfectly and positively correlated with the TI values with R 2 = 0.96 on 17 December 2009, but the data on 1 December 2009 showed a positive but very weak correlation between the TDR measurements and the TI values with R 2 = 0.07.The different results were due to the weather conditions around the sampling dates.On 1 December 2009, the ground was very dry as indicated by low TDR measurements with very little variation, while on 17 December 2009, the soils were saturated as indicated by the wide range of TDR measurements.

Comparison between VSLF and SWAT
Table 2 presents the seven VSLF calibration parameters, their affected hydrological variables, optimization goals being used, default values, original values calibrated by the model developer for a watershed in New York, and calibrated values from the Neshanic River Watershed.The calibrated values of the precipitation correction, SMIN and SMAX factors, were close to the original values.The calibrated values for the precipitation correction factor, runoff recess co-efficient, recess coefficient, and bypass coefficient were close to the default values.The only exception was the melt factor that determines the rate of snowpack melting.The calibrated value of 0.99 was much larger than the default value of 0.45 and the original value of 0.41 for the New York watershed.The difference in the melt factor value suggested the quite different snowpack melting processes among those watersheds due to their geographic location.The Neshanic River Watershed is located in the south of the New York watershed and much closer to the Atlantic Coast, and therefore snowpack melting would take place much faster.The SWAT calibration process and parameters were found in [55].The comparison between the SWAT-and VSLF-simulated and the observed daily, monthly, and annual runoff, baseflow, and streamflow at the USGS Reaville Gage Station during the calibration period 1997-2002 is presented in Table 3.The NSEs for the calibrated SWAT model were 0.58, 0.37, and 0.57 for the daily surface runoff, baseflow, and streamflow; 0.65, 0.61, and 0.68 for the monthly surface runoff, baseflow, and streamflow; and 0.48, −1.95, and 0.40 for the annual surface runoff, baseflow, and streamflow during the calibration period, respectively.The NSEs of the calibrated VSLF model were 0.54, 0.55, and 0.54 for the daily surface runoff, baseflow, and streamflow; 0.68, 0.75, and 0.74 for the monthly surface runoff, baseflow, and streamflow; and 0.85, 0.64, and 0.86 for the annual surface runoff, baseflow, and streamflow during the calibration period, respectively.For a monthly time step, the model performance is considered to be good for NSE values that fall in the range between 0.65 and 0.75 and to be very good for values between 0.75 and 1 [56].The reported NSE values of 0.68 and 0.74 for monthly streamflow indicated that both SWAT and VSLF models were well calibrated in the watershed.SWAT was better calibrated for daily streamflow while VSLF was better for monthly and annual streamflow.Both models underestimated the variability of streamflow as evidenced by the smaller standard deviations for the simulated streamflow.The calibrated models were then used to predict the surface runoff, baseflow, and streamflow for the validation period 2003-2008, which were then compared to the observed streamflow and the derived baseflow and runoff at the Reaville Gage Station in the period (Table 4).SWAT tended to underestimate and VSLF overestimate surface runoff and streamflow, while both were consistent in simulating baseflow showing by the means of simulated and observed flows.Both SWAT and VSLF tended to underestimate the variability of the flows showing lower standard deviations, but the variability captured by VSLF is closer to the observed variability than SWAT.Such findings were further confirmed by the hydrographs of the simulated and observed streamflow at the Reaville Gage Station (Figure 5).The NSEs for the SWAT model during the validation period were 0.33, 0.50, and 0.37 for the daily surface runoff, baseflow, and streamflow; 0.63, 0.71, and 0.68 for the monthly surface runoff, baseflow, and streamflow; and 0.52, 0.72, and 0.75 for the annual surface runoff, baseflow, and streamflow during the calibration period, respectively.The NSEs for the VSLF model were 0.30, 0.55, and 0.38 for the daily surface runoff, baseflow, and streamflow; 0.81, 0.59, and 0.79 for the monthly surface runoff, baseflow, and streamflow; and 0.38, 0.14, and 0.35 for the annual surface runoff, baseflow, and streamflow during the validation period, respectively.Following the same performance ratings provided by [56], these NSE values for the monthly streamflow indicated that the validation is very good for the VSLF model and good for the SWAT model.The differences in NSEs for the two models between the calibration and validation periods might be partially due to the changes in land uses in the watershed during the validation period while the model was still based on the land use conditions in 2002.Those measurements indicated both SWAT and VSLF achieved comparable predictions in surface runoff, baseflow, and streamflow in the watershed.

Conclusions
Natural hydrological processes coupled with historical land use management practices form various hydrological hotspots and/or CSAs in watersheds.The VSA hydrology and TI consistent with such a hydrological concept have been used to characterize hydrological sensitivity, and identify and delineate HSAs and/or CSAs in watersheds.Such characterization allows resource management practitioners to take spatially discrete land use management actions to protect water quality and aquatic ecosystem health in watersheds.It is essential to validate such a useful tool so that a suite of effective watershed management strategies can be built upon and implemented to improve management efficiency.
The study shows there are some practical difficulties to virtually validate the TI spatial pattern through field soil moisture sampling due to the significant impacts of micro-variations in soil, topography, and vegetation.The statistical analyses indicate a positive and high correlation between the TDR soil moisture measurements and the TI values, as shown on the binned graphs.Such correlation indicates that TI reasonably approximates the soil moisture condition, and therefore reasonably represents variability in runoff-generating potentials in landscapes.The lower correlation cases, as shown in Figures 3b and 4b, occurred because of the dry ground conditions on 1 December 2009 and the poor selection of the sampling locations in the Cold Brook Preserve.Future field studies on TI validation should select sampling locations with diverse TI values and should be conducted in wet seasons.For example, instead of randomly selecting sampling points on site, future validation studies might pre-select those sampling locations to cover a wider TI range and have a significant amount of locations in very different parts of a watershed to have a fuller representation of a watershed.
VSLF watershed modeling results indicate that incorporating TI into a lumped watershed model simulates the watershed hydrology reasonably well, as indicated by the goodness-of-fit in predicting surface runoff, baseflow, and streamflow in the Neshanic River Watershed during both the calibration and validation period.Such goodness-of-fit is comparable to the SWAT model in the same watershed using the same datasets.The simulation precision in watershed modeling is often achieved by dividing a watershed into smaller but different units with each unit representing a unique hydrological condition.Despite the similar modeling performance, the complexity of representing the watershed was quite different between the two models.In SWAT modeling, the same watershed was divided into 21 subbasins and 523 HRUs.Comparatively, VSLF used a much simpler watershed representation, i.e., intersections between land uses and 10 saturation class-based TI values.The watershed modeling incorporating TI can significantly improve the modeling efficiency.
In conclusion, TI is a good indicator for hydrological sensitivity in a watershed, provides a meaningful understanding on hydrological connections in a watershed and watershed hydrology, and therefore is a useful tool for watershed and water resource management.Such validation gives resource management practitioners the confidence to use TI and implement TI-based, spatially discrete land use management actions in watershed and water resource management practices.

Figure 1 .
Figure 1.The location of the study sites.

Figure 2 .
Figure 2. The spatial distribution of TI in the Neshanic River Watershed.

Figure 2 .
The outlet of Subbasin 12 coincides with the Reaville Gage Station.Similar to the VSLF modeling, the land use for the SWAT model was based on the 2002 NJDEP land use/cover data and the detailed land use activities in the watershed obtained by two rounds of agricultural land use inventories throughout the watershed during the period 2007-2008.The area percentage thresholds for defining HRUs were set at 1%, 20%, and 25% for land use, soil, and slope, respectively.As such, 625 individual HRUs were defined for the 25 subbasins delineated in the watershed.The drainage area of the Reaville Gage Station (i.e., the outlet of Subbasin 12) contains 21 subbasins and 523 HRUs.Similarly, the precipitation, maximum and minimum temperature during 1960-2008, solar radiation, and relative humidity during 1960-2008 in the Flemington weather station were used for the SWAT modeling.

Figure 5 .
Figure 5.Comparison of VSLF-and SWAT-simulated and USGS-observed streamflow at the Reaville Gage Station in the Neshanic River Watershed.

Table 1 .
Soil moisture measurements in two fields in Tewksbury Township.

Table 2 .
The default, original, and calibrated values for seven calibration parameters for the Variable Source Loading Function (VSLF) model.

Table 3 .
Comparison of daily, monthly, and annual streamflow and its two components, surface runoff and baseflow, during the calibration period 1997-2002.

Table 4 .
Comparison of daily, monthly, and annual streamflow and its two components, surface runoff and baseflow, during the validation period2003-2008.