FACTORS AFFECTING HERBICIDE YIELDS IN THE CHESAPEAKE BAY WATERSHED,JUNE 19941

ABSTRACT: Median concentrations and instantaneous yields of alachlor, metolachlor, atrazine, cyanazine, and simazine were generally highest at sites in the Lower Susquehanna River Basin and in agricultural subbasins. Instantaneous herbicide yields are related to land use, hydrogeologic setting, streamflow yield, and agricultural row cropping practices. The significance of these relations may be affected by the interdependence of the factors. The percentage of basin area planted in corn is the most influential factor in the prediction of herbicide yield. Instantaneous yields of all five herbicides measured in June 1994 related poorly to averaged 1990–94 herbicide use. Annually averaged herbicide-use data are too general to use as a predictor for short-term herbicide yields. An evaluation of factors affecting herbicide yields could be refined with more-current land use and land cover information and a more accurate estimate of the percentage of basin area planted in corn. Factors related to herbicide yields can be used to predict herbicide yields in other basins within the Chesapeake Bay watershed and to develop an estimate of herbicide loads to Chesapeake Bay.     

ANALYSIS AND MODELING OF LONG-TERM STREAM TEMPERATURES ON THE STEAMBOAT CREEK BASIN,OREGON: IMPLICATIONS FOR LAND USE AND FISH HABITAT1

ABSTRACT: Steamboat Creek basin is an important source of timber and provides crucial spawning and rearing habitat for anadromous steelhead trout (Oncorhynchus mykiss). Because stream temperatures are near the upper limit of tolerance for the survival of juvenile steelhead, the possible long-term effect of clear-cut logging on stream temperatures was assessed. Twenty-year (1969–1989) records of summer stream temperature and flow from four tributaries and two reaches of Steamboat Creek and Boulder Creek (a nearby unlogged watershed) were analyzed. Logging records for the Steamboat Creek basin and air temperature records also were used in the analysis. A time-series model of the components of stream temperature (seasonal cycle of solar radiation, air temperature, streamflow, an autoregressive term of order 1, and a linear trend variable) was fitted to the water-temperature data. The linear trend variable was significant in all the fitted models except Bend Creek (a tributary fed by cool ground-water discharge) and Boulder Creek. Because no trends in either climate (i.e., air temperature) or streamflow were found in the data, the trend variable was associated with the pre-1969 loss and subsequent regrowth of riparian vegetation and shading canopies.     

Multi-scale analysis of oxygen demand trends in an urbanizing Oregon watershed, USA

Human alteration of the landscape has an extensive influence on the biogeochemical processes that drive oxygen cycling in streams. We estimated trends from the mid-1990s to 2003, using the seasonal Mann-Kendall's test, for percent saturation dissolved oxygen (DO), chemical oxygen demand (COD), total Kjeldahl nitrogen (TKN), and ammonia-nitrogen (NH(3)-N) for 12 sites in the Rock Creek watershed, northwest Oregon, USA. In order to understand the influence of landscape change, scale, and stormwater runoff management on dissolved oxygen trends, we calculated land cover change through aerial photo interpretation at full-basin, local (near sample point) basin, and 100m stream buffer scales, for the years 1994 and 2000. Significant (p < or = 0.05) trends occurred in DO (increasing at five sites), COD (decreasing at seven sites), TKN (decreasing at five sites, increasing at one site), and NH(3)-N (decreasing at one site, increasing at one site). Significant land cover change occurred in agricultural land cover (-8% for the entire basin area) and residential land cover (+10% for the entire basin area) (p < or = 0.05). Correlation results indicated that: (1) forest cover negatively influenced COD at the full basin scale and positively influences NH(3)-N at local scales, (2) residential land cover influenced oxygen demand variables at local scales, (3) agricultural land cover did not influence oxygen demand, (4) local topography negatively influenced TKN and NH(3)-N, and (5) stormwater runoff management infrastructure correlated positively with COD at the local scale. This study indicates that landscape factors influencing DO conditions for the study streams act at multiple scales, suggesting that better knowledge of scale-process interactions can guide watershed managers' decision making in order to maintain improving water quality conditions.     

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

SIMULATING THREE-DIMENSIONAL GROUND WATER RESPONSE IN A SMALL MOUNTAINOUS WATERSHED1

ABSTRACT: Snowmelt from deep mountainous snowpacks is seldom rapid enough to exceed infiltration rates; thus, the source of streamflow in many mountainous watersheds is snowmelt recharge through shallow ground water systems. The hydrologic response and interaction between surface and sub-surface flow processes in these watersheds, which is controlled by basin structure, the spatial distribution of snowmelt, and the hydrogeology of the subsurface, are not well understood. The purpose of this study was to test a three-dimensional ground water model using simulated snowmelt input to simulate ground water response to spatially distributed snowmelt on the Upper Sheep Creek Watershed located within the Reynolds Creek Experimental Watershed in Southwestern Idaho. The model was used to characterize the mountainous aquifer and to delineate the subsurface flow mechanisms. Difficulty in finding a reasonable combination of grid spacing and time stepping within the model was encountered due to convergence problems with the Picard solution to the non-linear variably saturated ground water flow equations. Simulation results indicated that flow may be either unconfined or confined depending on inflow rate and hydrogeologic conditions in the watershed. The flow mechanism had a much faster response time when confined flow occurred. Response to snowmelt from a snow drift approximately 90 m away took only a few hours when flow was confined. Simulated results showed good agreement with piezometer measurements both in magnitude and timing; however, convergence problems with the Picard solution limited applicability of the model.     

ANALYSIS AND PREDICTION OF SURFACE RUNOFF IN AN URBANIZING WATERSHED USING SATELLITE IMAGERY1

ABSTRACT: This paper demonstrates how satellite image data [e.g., from Landsat 5 Thematic Mapper (TM)], in conjunction with an urban growth model and simple runoff calculations, can be used to estimate future surface runoff and, by implication, water quality within a watershed. To illustrate the method, predictions of land use change and surface runoff are shown for Spring Creek Watershed, a medium sized urbanizing watershed in Central Pennsylvania. Land cover classifications for this watershed were created from images for summertime 1986 and 1996 and subsequently used as input to the Clarke urban growth model, called SLEUTH, to predict land use changes to the year 2025. Simulations with this model show a progressive growth in the percentage of urban pixels and in impervious surface area in the watershed but also an increase in woodland, primarily in previously clear‐cut areas. Given that woodland area will continue to increase in area, surface runoff into Spring Creek is predicted to remain only slightly above present level. However, should the woodland amount fail to increase, surface runoff is then predicted to increase more significantly during the next 25 years. Finally, the concept of urban sprawl is addressed within the context of predicted increases in urbanization by relating the implied increase in impervious surface area to population density within the watershed.     

Impact of Climate and Land Use Change on Streamflow and Sediment Yield in a Snow‐Dominated Semiarid Mountainous Watershed

This study investigates the impact of climate and land use change on the magnitude and timing of streamflow and sediment yield in a snow‐dominated mountainous watershed in Salt Lake County, Utah using a scenario approach and the Hydrological Simulation Program — FORTRAN model for the 2040s (year 2035–2044) and 2090s (year 2085–2094). The climate scenarios were statistically and dynamically downscaled from global climate models. Land use and land cover (LULC) changes were estimated in two ways — from a regional planning scenario and from a deterministic model. Results indicate the mean daily streamflow in the Jordan River watershed will increase by an amount ranging from 11.2% to 14.5% in the 2040s and from 6.8% to 15.3% in the 2090s. The respective increases in sediment load in the 2040s and 2090s is projected to be 6.7% and 39.7% in the canyons and about 7.4% to 14.2% in the Jordan valley. The historical 50th percentile timing of streamflow and sediment load is projected to be shifted earlier by three to four weeks by mid‐century and four to eight weeks by late‐century. The projected streamflow and sediment load results establish a nonlinear relationship with each other and are highly sensitive to projected climate change. The predicted changes in streamflow and sediment yield will have implications for water supply, flood control and stormwater management.     

Watershed Evapotranspiration Increased due to Changes in Vegetation Composition and Structure Under a Subtropical Climate1

Abstract: Natural forests in southern China have been severely logged due to high human demand for timber, food, and fuels during the past century, but are recovering in the past decade. The objective of this study was to investigate how vegetation cover changes in composition and structure affected the water budgets of a 9.6‐km² Dakeng watershed located in a humid subtropical mountainous region in southern China. We analyzed 27 years (i.e., 1967‐1993) of streamflow and climate data and associated vegetation cover change in the watershed. Land use/land cover census and Normalized Difference of Vegetation Index (NDVI) data derived from remote sensing were used to construct historic land cover change patterns. We found that over the period of record, annual streamflow (Q) and runoff/precipitation ratio did not change significantly, nor did the climatic variables, including air temperature, Hamon’s potential evapotranspiration (ET), pan evaporation, sunshine hours, and radiation. However, annual ET estimated as the differences between P and Q showed a statistically significant increasing trend. Overall, the NDVI of the watershed had a significant increasing trend in the peak spring growing season. This study concluded that watershed ecosystem ET increased as the vegetation cover shifted from low stock forests to shrub and grasslands that had higher ET rates. A conceptual model was developed for the study watershed to describe the vegetation cover‐streamflow relationships during a 50‐year time frame. This paper highlighted the importance of eco‐physiologically based studies in understanding transitory, nonstationary effects of deforestation or forestation on watershed water balances.     

APPLICATION OF THE RHESSys MODEL TO A CALIFORNIA SEMIARID SHRUBLAND WATERSHED1

ABSTRACT: Distributed hydrologic models which link seasonal streamflow and soil moisture patterns with spatial patterns of vegetation are important tools for understanding the sensitivity of Mediterranean type ecosystems to future climate and land use change. RHESSys (Regional Hydro‐Ecologic Simulation System) is a coupled spatially distributed hydroecological model that is designed to be able to represent these feedbacks between hydrologic and vegetation carbon and nutrient cycling processes. However, RHESSys has not previously been applied to semiarid shrubland watersheds. In this study, the hydrologic submodel of RHESSys is evaluated by comparing model predictions of monthly and annual streamflow to stream gage data and by comparing RHESSys behavior to that of another hydrologic model of similar complexity, MIKESHE, for a 34 km² watershed near Santa Barbara, California. In model intercomparison, the differences in predictions of temporal patterns in streamflow, sensitivity of model predictions to calibration parameters and landscape representation, and differences in model estimates of soil moisture patterns are explored. Results from this study show that both models adequately predict seasonal patterns of streamflow response relative to observed data, but differ significantly in terms of estimates of soil moisture patterns and sensitivity of those patterns to the scale of landscape tessellation used to derive spatially distributed elements. This sensitivity has implications for implementing RHESSys as a tool to investigate interactions between hydrology and ecosystem processes.     

STORM RUNOFF PREDICTION BASED ON A SPATIALLY DISTRIBUTED TRAVEL TIME METHOD UTILIZING REMOTE SENSING AND GIS1

ABSTRACT: In this study, remotely sensed data and geographic information system (GIS) tools were used to estimate storm runoff response for Simms Creek watershed in the Etonia basin in northeast Florida. Land cover information from digital orthophoto quarter quadrangles (DOQQ), and enhanced thematic mapper plus (ETM+) were analyzed for the years 1990, 1995, and 2000. The corresponding infiltration excess runoff response of the study area was estimated using the U.S. Department of Agriculture (USDA), Natural Resources Conservation Service Curve Number (NRCS‐CN) method. A digital elevation model (DEM)/GIS technique was developed to predict stream response to runoff events based on the travel time from each grid cell to the watershed outlet. A comparison of predicted to observed stream response shows that the model predicts the total runoff volume with an efficiency of 0.98, the peak flow rate at an efficiency of 0.85, and the full direct runoff hydrograph with an average efficiency of 0.65. The DEM/GIS travel time model can be used to predict the runoff response of ungaged watersheds and is useful for predicting runoff hydrographs resulting from proposed large scale changes in the land use.     

Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain1

Epps, Thomas H., Daniel R. Hitchcock, Anand D. Jayakaran, Drake R. Loflin, Thomas M. Williams, and Devendra M. Amatya, 2012. Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12000 Abstract: Hydrologic monitoring was conducted in two first‐order lower coastal plain watersheds in South Carolina, United States, a region with increasing growth and land use change. Storm events over a three‐year period were analyzed for direct runoff coefficients (ROC) and the total storm response (TSR) as percent rainfall. ROC calculations utilized an empirical hydrograph separation method that partitioned total streamflow into sustained base flow and direct runoff components. ROC ratios ranged from 0 to 0.32 on the Upper Debidue Creek (UDC) watershed and 0 to 0.57 on Watershed 80 (WS80); TSR results ranged from 0 to 0.93 at UDC and 0.01 to 0.74 at WS80. Variability in event runoff generation was attributed to seasonal trends in water table elevation fluctuation as regulated by evapotranspiration. Groundwater elevation breakpoints for each watershed were identified based on antecedent water table elevation, streamflow, ROCs, and TSRs. These thresholds represent the groundwater elevation above which event runoff generation increased sharply in response to rainfall. For effective coastal land use decision making, baseline watershed hydrology must be understood to serve as a benchmark for management goals, based on both seasonal and event‐based surface and groundwater interactions.     

FLOW ANALYSIS OF LANDSLIDE DAMMED LAKE WATERSHEDS: A CASE STUDY1

Abstract: The Chi-Chi earthquake, which occurred on September 21, 1999, and had a magnitude of 7.3 on the Richter scale, resulted in an extensive landslide that blocked the Ching-Shui Creek in Taiwan, forming a large lake with a storage volume of 40 million m³. This paper describes an analytical procedure used to perform flow analysis of the Tsao-Ling watershed, which includes the new landslide dammed lake. In this study, a digital elevation model was applied to obtain the watershed geomorphic factors and stage-area storage function of the landslide dammed lake. Satellite images were used to identify the landslide area and the land cover change that occurred as a result of the earthquake. Two topography-based runoff models were applied for long term and short term streamflow analyses of the watershed because the watershed upstream of the landslide dam was ungauged. The simulated daily flow and storm runoff were verified using limited available measured data in the watershed, and good agreement was obtained. The proposed analytical procedure for flow analysis is considered promising for application to other landslide dammed lake watersheds.     

MODELING ACTUAL EVAPOTRANSPIRATION FROM FORESTED WATERSHEDS ACROSS THE SOUTHEASTERN UNITED STATES1

ABSTRACT: About 50 to 80 percent of precipitation in the southeastern United States returns to the atmosphere by evapotranspiration. As evapotranspiration is a major component in the forest water balances, accurately quantifying it is critical to predicting the effects of forest management and global change on water, sediment, and nutrient yield from forested watersheds. However, direct measurement of forest evapotranspiration on a large basin or a regional scale is not possible. The objectives of this study were to develop an empirical model to estimate long‐term annual actual evapotranspiration (ART) for forested watersheds and to quantify spatial AET patterns across the southeast. A geographic information system (GIS) database including land cover, daily streamflow, and climate was developed using long term experimental and monitoring data from 39 forested watersheds across the region. Using the stepwise selection method implemented in a statistical modeling package, a long term annual AET model was constructed. The final multivariate linear model includes four independent variables—annual precipitation, watershed latitude, watershed elevation, and percentage of forest coverage. The model has an adjusted R² of 0.794 and is sufficient to predict long term annual ART for forested watersheds across the southeastern United States. The model developed by this study may be used to examine the spatial variability of water availability, estimate annual water loss from mesoscale watersheds, and project potential water yield change due to forest cover change.     

DEVELOPMENT AND TESTING OF A SNOWMELT-RUNOFF FORECASTING TECHNIQUE1

ABSTRACT: The snowmelt-runoff model (SRM) was used to produce accurate simulations of streamfiow during the snowmelt period (April-September) for ten years on the Rio Grande Basin (3419 km²) near Del Norte, Colorado, U.S.A. In order to use SRM in the forecast situation, it was necessary to develop a family of snow cover depletion curves for each elevation zone based on accumulated snow water equivalent on April 1. Selection of an appropriate curve for a particular year from snow course measurements allows input of the daily snow cover extent to SRM for forecast purposes. Data from three years (1980, 1981, and 1985) were used as a quasi-forecast test of the procedure. In these years forecasted snow cover extent data were input to SRM, but observed temperature and precipitation data were used. The resulting six-month hydrographs were very similar to the hydrographs in the ten simulation years previously tested based on comparisons of performance evaluation criteria. Based on this result, the Soil Conservation Service (SCS) requested SRM forecasts for 1987 on the Rio Grande. Using the same procedure but with SCS estimated temperature and precipi-tation data, SRM produced a forecast hydrograph that had a r²= 0.82 and difference in seasonal volume of 4.4 percent. To approximate actual operational conditions, SRM computed daily flows were updated every seven days with measured flows. The resulting forecast hydrograph had a R²= 0.90 and a difference in volume of 3.5 percent. The method developed needs to be refined and tested on additional years and basins, but the approach appears to be applicable to operational runoff forecasting using remote sensing data.     

Evaluation of Short-to-Medium Range Streamflow Forecasts Obtained Using an Enhanced Version of SRM1

Harshburger, Brian J., Karen S. Humes, Von P. Walden, Brandon C. Moore, Troy R. Blandford, and Albert Rango, 2010. Evaluation of Short-to-Medium Range Streamflow Forecasts Obtained Using an Enhanced Version of SRM. Journal of the American Water Resources Association (JAWRA) 46(3):603-617. DOI: 10.1111/j.1752-1688.2010.00437.x Abstract: As demand for water continues to escalate in the western United States, so does the need for accurate streamflow forecasts. Here, we describe a methodology for generating short-to-medium range (1 to 15 days) streamflow forecasts using an enhanced version of the Snowmelt Runoff Model (SRM), snow-covered area data derived from MODIS products, data from Snow Telemetry stations, and meteorological forecasts. The methodology was tested on three mid-elevation, snowmelt-dominated basins ranging in size from 1,600 to 3,500 km². To optimize the model performance and aid in its operational implementation, two enhancements have been made to SRM: (1) the use of an antecedent temperature index method to track snowpack cold content, and (2) the use of both maximum and minimum critical temperatures to partition precipitation into rain, snow, or a mixture of rain and snow. The comparison of retrospective model simulations with observed streamflow shows that the enhancements significantly improve the model performance. Streamflow forecasts generated using the enhanced version of the model compare well with the observed streamflow for the earlier leadtimes; forecast performance diminishes with leadtime due to errors in the meteorological forecasts. The three basins modeled in this research are typical of many mid-elevation basins throughout the American West, thus there is potential for this methodology to be applied successfully to other mountainous basins.     

Water Quality Signals from Rural Land Use and Exurbanization in a Mountain Landscape: What's Clear and What's Confounded?

In mountainous landscapes with high climatic and geomorphic variability, how do rural land uses and exurbanization alter hydrology and water quality? We evaluated effects of rural land use and exurbanization on streamflows, suspended sediment concentrations and loads, specific conductance, and summer water temperatures in 12 streams and rivers within the Upper Little Tennessee River basin in the southern Appalachian Mountains. Eleven streams featured low levels of development (>61% forest cover) but differed in land use patterning, basin size, annual precipitation, and watershed morphology. One urban stream, located within the largest town in the basin, provided the high development comparative endpoint. Even low levels of rural development and exurbanization were associated with substantial increases in suspended sediment concentrations, sediment loads, and summer stream temperature daily maxima and diurnal variation. Observed summer temperature increases were much larger than would be expected due to global climate change over the next century. Specific conductance was idiosyncratic among the smaller streams. These water quality changes were not accompanied by streamflow changes that were discernible amid the high natural variation in precipitation and geomorphology. The water quality findings suggest the need for applying the best management practices, including riparian buffers, to even low levels of rural development.     

Hydrologic Modeling Uncertainty Resulting From Land Cover Misclassification1

Abstract: A stochastic, spatially explicit method for assessing the impact of land cover classification error on distributed hydrologic modeling is presented. One‐hundred land cover realizations were created by systematically altering the North American Landscape Characterization land cover data according to the dataset’s misclassification matrix. The matrix indicates the probability of errors of omission in land cover classes and is used to assess the uncertainty in hydrologic runoff simulation resulting from parameter estimation based on land cover. These land cover realizations were used in the GIS‐based Automated Geospatial Watershed Assessment tool in conjunction with topography and soils data to generate input to the physically‐based Kinematic Runoff and Erosion model. Uncertainties in modeled runoff volumes resulting from these land cover realizations were evaluated in the Upper San Pedro River basin for 40 watersheds ranging in size from 10 to 100 km² under two rainfall events of differing magnitudes and intensities. Simulation results show that model sensitivity to classification error varies directly with respect to watershed scale, inversely to rainfall magnitude and are mitigated or magnified by landscape variability depending on landscape composition.     

Drivers of future streamflow changes in watersheds across the Northeastern United States

Accurate projections of streamflow, which have implications for flooding, water resources, hydropower, and ecosystems, are critical to climate change adaptation and require an understanding of streamflow sensitivity to climate drivers. The northeastern United States has experienced a dramatic increase in extreme precipitation over the past 25 years; however, the effects of these changes, as well as changes in other drivers of streamflow, remain unclear. Here, we use a random forest model forced with a regional climate model to examine historical and future streamflow dynamics of four watersheds across the Northeast. We find that streamflow in the cold season (November–May) is primarily driven by 3-day rainfall and antecedent wetness (Antecedent Precipitation Index) in three rainfall-dominant watersheds, and 30-day rainfall, antecedent wetness, and 30-day snowmelt in the fourth, more snowmelt-dominated watershed. In the warm season (June–October), streamflow is driven by antecedent wetness and rainfall in all watersheds. By the end of the century (2070–2099), cold season streamflow depends on the importance placed on snow in the machine learning model, with changes ranging from −7% (with snow) to +40% (without snow) in a single watershed. Simulated future warm season streamflow increases in two watersheds (56% and 193%) due to increased precipitation and antecedent soil wetness, but decreases in the other two watersheds (−6% and −27%) due to reduced precipitation.     

A Simple Process‐Based Snowmelt Routine to Model Spatially Distributed Snow Depth and Snowmelt in the SWAT Model1

Abstract: We present a method to integrate a process‐based (PB) snowmelt model that requires only daily temperature and elevation information into the Soil and Water Assessment Tool (SWAT) model. The model predicts the spatiotemporal snowpack distribution without adding additional complexity, and in fact reduces the number of calibrated parameters. To demonstrate the utility of the PB model, we calibrate the PB and temperature‐index (TI) SWAT models to optimize agreement with stream discharge on a 46‐km² watershed in northwestern Idaho, United States, for 10 individual years and use the calibrated parameters for the year with the best agreement to run the model for 15 remaining years. Stream discharge predictions by the PB and TI model were similar, although the PB model simulated snowmelt more accurately than the TI model for the remaining 15‐year period. Spatial snow distributions predicted by the PB model better matched observations from LandSat imagery and a SNOTEL station. Results for this watershed show that including PB snowmelt in watershed models is feasible, and calibration of TI‐based watershed models against discharge can incorrectly predict snow cover.     

INTEGRATING LANDSCAPE ASSESSMENT AND HYDROLOGIC MODELING FOR LAND COVER CHANGE ANALYSIS1

ABSTRACT: Significant land cover changes have occurred in the watersheds that contribute runoff to the upper San Pedro River in Sonora, Mexico, and southeast Arizona. These changes, observed using a series of remotely sensed images taken in the 1970s, 1980s, and 1990s, have been implicated in the alteration of the basin hydrologic response. The Cannonsville subwatershed, located in the Catskill/Delaware watershed complex that delivers water to New York City, provides a contrast in land cover change. In this region, the Cannonsville watershed condition has improved over a comparable time period. A landscape assessment tool using a geographic information system (GIS) has been developed that automates the parameterization of the Soil and Water Assessment Tool (SWAT) and KINEmatic Runoff and EROSion (KINEROS) hydrologic models. The Automated Geospatial Watershed Assessment (AGWA) tool was used to prepare parameter input files for the Upper San Pedro Basin, a subwatershed within the San Pedro undergoing significant changes, and the Cannonsville watershed using historical land cover data. Runoff and sediment yield were simulated using these models. In the Cannonsville watershed, land cover change had a beneficial impact on modeled watershed response due to the transition from agriculture to forest land cover. Simulation results for the San Pedro indicate that increasing urban and agricultural areas and the simultaneous invasion of woody plants and decline of grasslands resulted in increased annual and event runoff volumes, flashier flood response, and decreased water quality due to sediment loading. These results demonstrate the usefulness of integrating remote sensing and distributed hydrologic models through the use of GIS for assessing watershed condition and the relative impacts of land cover transitions on hydrologic response.