Hydrologic Landscape Classification to Estimate Bristol Bay,Alaska Watershed Hydrology

Hydrologic landscapes (HLs) have proven to be a useful tool for broad scale assessment and classification of landscapes across the United States as they help organize larger geographical areas into areas of similar hydrologic characteristics. We developed a HL classification for the Bristol Bay watershed of southwest Alaska that incorporates indices of annual climate and seasonality, terrain, geology, and the influences of large lakes and glaciers. A HL classification is particularly useful in this large watershed because of its hydrologic and landscape variability, important salmon fishery, variety of environmental and potential anthropogenic stressors, and lack of widespread hydrologic data. Following creation of Bristol Bay basin‐wide HL classes, we compared the HL distributions within watersheds grouped by two calculated runoff parameters derived from available long‐term streamflow records and found HL distributions within these groups provided predictive insight on hydrologic behavior. Using these developed runoff groups, we estimated expected hydrologic behavior in watersheds across the larger Bristol Bay watershed that lacked gauged streamflow records. The HL approach provides a scientific basis for estimating the first‐order hydrologic behavior of watersheds and landscapes that lack detailed hydrologic information.     

Hydrologic Landscape Characterization for the Pacific Northwest,USA

We update the Wigington et al. (2013) hydrologic landscape (HL) approach to make it more broadly applicable and apply the revised approach to the Pacific Northwest (PNW; i.e., Oregon, Washington, and Idaho). Specific changes incorporated are the use of assessment units based on National Hydrography Dataset Plus V2 catchments, a modified snowmelt model validated over a broader area, an aquifer permeability index that does not require preexisting aquifer permeability maps, and aquifer and soil permeability classes based on uniform criteria. Comparison of Oregon results for the revised and original approaches found fewer and larger assessment units, loss of summer seasonality, and changes in rankings and proportions of aquifer and soil permeability classes. Differences could be explained by three factors: an increased assessment unit size, a reduced number of permeability classes, and use of smaller cutoff values for the permeability classes. The distributions of the revised HLs in five groups of Oregon rivers were similar to the original HLs but less variable. The improvements reported here should allow the revised HL approach to be applied more often in situations requiring hydrologic classification and allow greater confidence in results. We also apply the map results to the development of hydrologic landscape regions.     

Revealing the Diversity of Natural Hydrologic Regimes in California with Relevance for Environmental Flows Applications

Alterations to flow regimes for water management objectives have degraded river ecosystems worldwide. These alterations are particularly profound in Mediterranean climate regions such as California with strong climatic variability and riverine species highly adapted to the resulting flooding and drought disturbances. However, defining environmental flow targets for Mediterranean rivers is complicated by extreme hydrologic variability and often intensive water management legacies. Improved understanding of the diversity of natural streamflow patterns and their spatial arrangement across Mediterranean regions is needed to support the future development of effective flow targets at appropriate scales for management applications with minimal resource and data requirements. Our study addresses this need through the development of a spatially explicit reach‐scale hydrologic classification for California. Dominant hydrologic regimes and their physio‐climatic controls are revealed, using available unimpaired and naturalized streamflow time‐series and generally publicly available geospatial datasets. This methodology identifies eight natural flow classes representing distinct flow sources, hydrologic characteristics, and catchment controls over rainfall‐runoff response. The study provides a broad‐scale hydrologic framework upon which flow‐ecology relationships could subsequently be established towards reach‐scale environmental flows applications in a complex, highly altered Mediterranean region.     

Physical and Temporal Isolation of Mountain Headwater Streams in the Western Mojave Desert,Southern California1

Abstract: Streams draining mountain headwater areas of the western Mojave Desert are commonly physically isolated from downstream hydrologic systems such as springs, playa lakes, wetlands, or larger streams and rivers by stream reaches that are dry much of the time. The physical isolation of surface flow in these streams may be broken for brief periods after rainfall or snowmelt when runoff is sufficient to allow flow along the entire stream reach. Despite the physical isolation of surface flow in these streams, they are an integral part of the hydrologic cycle. Water infiltrated from headwater streams moves through the unsaturated zone to recharge the underlying ground‐water system and eventually discharges to support springs, streamflow, isolated wetlands, or native vegetation. Water movement through thick unsaturated zones may require several hundred years and subsequent movement through the underlying ground‐water systems may require many thousands of years – contributing to the temporal isolation of mountain headwater streams.     

INCREASED BASEFLOW IN IOWA OVER THE SECOND HALF OF THE 20TH CENTURY1

ABSTRACT: Historical trends in annual discharge characteristics were evaluated for 11 gauging stations located throughout Iowa. Discharge records from nine eight‐digit hydrologic unit code (HUC‐8) watersheds were examined for the period 1940 to 2000, whereas data for two larger river systems (Cedar and Des Moines Rivers) were examined for a longer period of record (1903 to 2000). In nearly all watersheds evaluated, annual base flow, annual minimum flow, and the annual base flow percentage significantly increased over time. Some rivers also exhibited increasing trends in total annual discharge, whereas only the Maquoketa River had significantly decreased annual maximum flows. Regression of stream discharge versus precipitation indicated that more precipitation is being routed into streams as base flow than as storm flow in the second half of the 20th Century. Reasons for the observed stream flow trends are hypothesized to include improved conservation practices, greater artificial drainage, increasing row crop production, and channel incision. Each of these reasons is consistent with the observed trends, and all are likely responsible to some degree in most watersheds.     

Evaluating Flow Metric‐Based Stream Classification Systems to Support the Determination of Ecological Flows in North Carolina

Hydroecological classification systems are typically based on an assemblage of streamflow metrics and seek to divide streams into ecologically relevant classes. Assignment of streams to classes is suggested as an initial step in the process of establishing ecological flow standards. We used two distinct hydroecological river classification systems available within North Carolina to evaluate the ability of a hydrologic model to assign the same classes as those determined by observed streamflows and to assess the transferability of such systems to ungaged streams. Class assignments were examined by rate of overall matches, rate of class matches, spatial variability in matches, and time period used in class assignment. The findings of this study indicate assignments of stream class: (1) are inconsistent among different classification systems; (2) differ between observed and modeled data; and (3) are sensitive to the period of record within observed data. One clear source of inconsistency/sensitivity in class assignments lies with the use of threshold values for metrics that distinguish stream classes, such that even small changes in metric values can result in different class assignments. Because these two hydroecological classification systems are representative of other classification systems that rely on quantitative decision thresholds, it can be surmised that the use of such systems based on stream flow metrics is not a reliable approach for guiding ecological flow determinations.     

Temporal Changes in Streamflow and Attribution of Changes to Climate and Landuse in Wisconsin Watersheds

Previous historic trends analyses on 21st Century hydrologic data in the United States generally focus on annual flow statistics and have continued to use USGS hydro‐climatic data network (HCDN) stations, although post‐1988 diversions and runoff regulations are not reflected in the HCDN. Using a more recent dataset, Geospatial Attributes of Gages for Evaluating Streamflow, version II (GAGES II), compiled by Falcone (2012), which includes more watersheds with reference conditions, a comprehensive analysis of changes in seasonal, and annual streamflow in Wisconsin watersheds is demonstrated. Given the pronounced influence of seasonal hydrology in Wisconsin watersheds, the objective of this study is to elucidate the nature of temporal (annual, seasonal, and monthly) changes in runoff. Considerable temporal and regional variability was found in annual and seasonal streamflow changes between the two historic periods 1951‐1980 and 1981‐2010 considered in the study. For example, the northern watersheds show relatively small changes in streamflow discharge ranging from ?6.0 to 4.2%, while the southern watersheds show relatively large increases in streamflow discharge ranging from 13.1 to 18.2%. To apportion streamflow changes to climate and nonclimatic factors, a method based on potential evapotranspiration changes is demonstrated. Results show that nonclimatic factors account for more than 60% of changes in annual runoff in Wisconsin watersheds considered in the study.     

Hydrologic Connectivity and the Contribution of Stream Headwaters to Ecological Integrity at Regional Scales1

Abstract: Cumulatively, headwater streams contribute to maintaining hydrologic connectivity and ecosystem integrity at regional scales. Hydrologic connectivity is the water‐mediated transport of matter, energy and organisms within or between elements of the hydrologic cycle. Headwater streams compose over two‐thirds of total stream length in a typical river drainage and directly connect the upland and riparian landscape to the rest of the stream ecosystem. Altering headwater streams, e.g., by channelization, diversion through pipes, impoundment and burial, modifies fluxes between uplands and downstream river segments and eliminates distinctive habitats. The large‐scale ecological effects of altering headwaters are amplified by land uses that alter runoff and nutrient loads to streams, and by widespread dam construction on larger rivers (which frequently leaves free‐flowing upstream portions of river systems essential to sustaining aquatic biodiversity). We discuss three examples of large‐scale consequences of cumulative headwater alteration. Downstream eutrophication and coastal hypoxia result, in part, from agricultural practices that alter headwaters and wetlands while increasing nutrient runoff. Extensive headwater alteration is also expected to lower secondary productivity of river systems by reducing stream‐system length and trophic subsidies to downstream river segments, affecting aquatic communities and terrestrial wildlife that utilize aquatic resources. Reduced viability of freshwater biota may occur with cumulative headwater alteration, including for species that occupy a range of stream sizes but for which headwater streams diversify the network of interconnected populations or enhance survival for particular life stages. Developing a more predictive understanding of ecological patterns that may emerge on regional scales as a result of headwater alterations will require studies focused on components and pathways that connect headwaters to river, coastal and terrestrial ecosystems. Linkages between headwaters and downstream ecosystems cannot be discounted when addressing large‐scale issues such as hypoxia in the Gulf of Mexico and global losses of biodiversity.     

POTENTIAL IMPACTS OF CLIMATE CHANGE ON CALIFORNIA HYDROLOGY1

ABSTRACT: Previous reports based on climate change scenarios have suggested that California will be subjected to increased wintertime and decreased summertime streamflow. Due to the uncertainty of projections in future climate, a new range of potential climatological future temperature shifts and precipitation ratios is applied to the Sacramento Soil Moisture Accounting Model and Anderson Snow Model in order to determine hydrologic sensitivities. Two general circulation models (GCMs) were used in this analysis: one that is warm and wet (HadCM2 run 1) and one that is cool and dry (PCM run B06.06), relative to the GCM projections for California that were part of the Third Assessment Report of the Intergovernmental Panel on Climate Change. A set of specified incremental temperature shifts from 1.5°C to 5.0°C and precipitation ratios from 0.70 to 1.30 were also used as input to the snow and soil moisture accounting models, providing for additional scenarios (e.g., warm/dry, cool/wet). Hydrologic calculations were performed for a set of California river basins that extend from the coastal mountains and Sierra Nevada northern region to the southern Sierra Nevada region; these were applied to a water allocation analysis in a companion paper. Results indicate that for all snow‐producing cases, a larger proportion of the streamflow volume will occur earlier in the year. The amount and timing is dependent on the characteristics of each basin, particularly the elevation. Increased temperatures lead to a higher freezing line, therefore less snow accumulation and increased melting below the freezing height. The hydrologic response varies for each scenario, and the resulting solution set provides bounds to the range of possible change in streamflow, snowmelt, snow water equivalent, and the change in the magnitude of annual high flows. An important result that appears for all snowmelt driven runoff basins, is that late winter snow accumulation decreases by 50 percent toward the end of this century.     

Mid‐Range Streamflow Forecasts Based on Climate Modeling – Statistical Correction and Evaluation1

Abstract: Mid‐range streamflow predictions are extremely important for managing water resources. The ability to provide mid‐range (three to six months) streamflow forecasts enables considerable improvements in water resources system operations. The skill and economic value of such forecasts are of great interest. In this research, output from a general circulation model (GCM) is used to generate hydrologic input for mid‐range streamflow forecasts. Statistical procedures including: (1) transformation, (2) correction, (3) observation of ensemble average, (4) improvement of forecast, and (5) forecast skill test are conducted to minimize the error associated with different spatial resolution between the large‐scale GCM and the finer‐scale hydrologic model and to improve forecast skills. The accuracy of a streamflow forecast generated using a hydrologic model forced with GCM output for the basin was evaluated by forecast skill scores associated with the set of streamflow forecast values in a categorical forecast. Despite the generally low forecast skill score exhibited by the climate forecasting approach, precipitation forecast skill clearly improves when a conditional forecast is performed during the East Asia summer monsoon, June through August.     

Use of Hydrologic Landscape Classification to Diagnose Streamflow Predictability in Oregon

We implement a spatially lumped hydrologic model to predict daily streamflow at 88 catchments within the state of Oregon and analyze its performance using the Oregon Hydrologic Landscape (OHL) classification. OHL is used to identify the physio‐climatic conditions that favor high (or low) streamflow predictability. High prediction catchments (Nash‐Sutcliffe efficiency of (NS) > 0.75) are mainly classified as rain dominated with very wet climate, low aquifer permeability, and low to medium soil permeability. Most of them are located west of the Cascade Mountain Range. Conversely, most low prediction catchments (NS < 0.6) are classified as snow‐dominated with high aquifer permeability and medium to high soil permeability. They are mainly located in the volcano‐influenced High Cascades region. Using a subset of 36 catchments, we further test if class‐specific model parameters can be developed to predict at ungauged catchments. In most catchments, OHL class‐specific parameters provide predictions that are on par with individually calibrated parameters (NS decline < 10%). However, large NS declines are observed in OHL classes where predictability is not high enough. Results suggest higher uncertainty in rain‐to‐snow transition of precipitation phase and external gains/losses of deep groundwater are major factors for low prediction in Oregon. Moreover, regionalized estimation of model parameters is more useful in regions where conditions favor good streamflow predictability.     

Spatial Calibration and Temporal Validation of Flow for Regional Scale Hydrologic Modeling1

Abstract: Physically based regional scale hydrologic modeling is gaining importance for planning and management of water resources. Calibration and validation of such regional scale model is necessary before applying it for scenario assessment. However, in most regional scale hydrologic modeling, flow validation is performed at the river basin outlet without accounting for spatial variations in hydrological parameters within the subunits. In this study, we calibrated the model to capture the spatial variations in runoff at subwatershed level to assure local water balance, and validated the streamflow at key gaging stations along the river to assure temporal variability. Ohio and Arkansas‐White‐Red River Basins of the United States were modeled using Soil and Water Assessment Tool (SWAT) for the period from 1961 to 1990. R² values of average annual runoff at subwatersheds were 0.78 and 0.99 for the Ohio and Arkansas Basins. Observed and simulated annual and monthly streamflow from 1961 to 1990 is used for temporal validation at the gages. R² values estimated were greater than 0.6. In summary, spatially distributed calibration at subwatersheds and temporal validation at the stream gages accounted for the spatial and temporal hydrological patterns reasonably well in the two river basins. This study highlights the importance of spatially distributed calibration and validation in large river basins.     

Adaptive Water Resource Planning in the South Saskatchewan River Basin: Use of Scenarios of Hydroclimatic Variability and Extremes

The South Saskatchewan River Basin is one of Canada's most threatened watersheds, with water supplies in most subbasins over‐allocated. In 2013, stakeholders representing irrigation districts, the environment, and municipalities collaborated with researchers and consultants to explore opportunities to improve the resiliency of the management of the Oldman and South Saskatchewan River subbasins. Streamflow scenarios for 2025‐2054 were constructed by the novel approach of regressing historical river flows against indices of large‐scale ocean‐atmosphere climate oscillations to derive statistical streamflow models, which were then run using projected climate indices from global climate models. The impacts of some of the most extreme scenarios were simulated using the hydrologic mass‐balance model Operational Analysis and Simulation of Integrated Systems (OASIS). Based on stakeholder observations, the project participants proposed and evaluated potential risk management and adaption strategies, e.g., modifying existing infrastructure, building new infrastructure, changing operations to supplement environmental flows, reducing demand, and sharing supply. The OASIS model was applied interactively at live modeling sessions with stakeholders to explore practical adaptation strategies. Our results, which serve as recommendations for policy makers, showed that forecast‐based rationing together with new expanded storage could dramatically reduce water shortages.     

Modeling Hydrology in a Small Rocky Mountain Watershed Serving Large Urban Populations1

Abstract: This article describes the development of a calibrated hydrologic model for the Blue River watershed (867 km²) in Summit County, Colorado. This watershed provides drinking water to over a third of Colorado’s population. However, more research on model calibration and development for small mountain watersheds is needed. This work required integration of subsurface and surface hydrology using GIS data, and included aspects unique to mountain watersheds such as snow hydrology, high ground‐water gradients, and large differences in climate between the headwaters and outlet. Given the importance of this particular watershed as a major urban drinking‐water source, the rapid development occurring in small mountain watersheds, and the importance of Rocky Mountain water in the arid and semiarid West, it is useful to describe calibrated watershed modeling efforts in this watershed. The model used was Soil and Water Assessment Tool (SWAT). An accurate model of the hydrologic cycle required incorporation of mountain hydrology‐specific processes. Snowmelt and snow formation parameters, as well as several ground‐water parameters, were the most important calibration factors. Comparison of simulated and observed streamflow hydrographs at two U.S. Geological Survey gaging stations resulted in good fits to average monthly values (0.71 Nash‐Sutcliffe coefficient). With this capability, future assessments of point‐source and nonpoint‐source pollutant transport are possible.     

Heterogeneity of Hydrologic Response in Urban Watersheds1

Meierdiercks, Katherine L., James A. Smith, Mary Lynn Baeck, and Andrew J. Miller, 2010. Heterogeneity of Hydrologic Response in Urban Watersheds. Journal of the American Water Resources Association (JAWRA) 46(6):1221–1237. DOI: 10.1111/j.1752-1688.2010.00487.x Abstract: The changing patterns of streamflow associated with urbanization are examined through analyses of discharge and rainfall records for the study watersheds of the Baltimore Ecosystem Study (BES). Analyses utilize a decade (1999-2008) of observations from a dense U.S. Geological Survey stream gaging network and Hydro-NEXRAD radar rainfall fields. The principal study watershed of the BES is Gwynns Falls (171 km²). Focus is given to two Gwynns Falls basins with contrasting patterns and histories of development, Dead Run and Upper Gwynns Falls. The sharp contrasts in streamflow properties between the basins reflect the differences in urban development prior to implementation of stormwater management regulations (much of Dead Run) and development for which stormwater management is an integral part of the hydrologic system (Upper Gwynns Falls). The mean annual runoff in Dead Run (558 mm) is 35% larger than that of Upper Gwynns Falls; larger contrasts in runoff properties typify the “warm season” and are linked to storm event hydrologic response. Spatial heterogeneities in storm event response are reflected in seasonal and diurnal properties of streamflow. Analyses of storm event response are presented for June 2006, during which monthly rainfall over the BES region ranged from less than 150 to more than 500 mm. Baisman Run, the BES forest reference watershed, and Moores Run, a highly urbanized watershed in Baltimore City, provide “end-member” representations of urban impacts on streamflow.     

How Hydrologic Connectivity Regulates Water Quality in River Corridors

Downstream flow in rivers is repeatedly delayed by hydrologic exchange with off‐channel storage zones where biogeochemical processing occurs. We present a dimensionless metric that quantifies river connectivity as the balance between downstream flow and the exchange of water with the bed, banks, and floodplains. The degree of connectivity directly influences downstream water quality — too little connectivity limits the amount of river water exchanged and leads to biogeochemically inactive water storage, while too much connectivity limits the contact time with sediments for reactions to proceed. Using a metric of reaction significance based on river connectivity, we provide evidence that intermediate levels of connectivity, rather than the highest or lowest levels, are the most efficient in removing nitrogen from Northeastern United States’ rivers. Intermediate connectivity balances the frequency, residence time, and contact volume with reactive sediments, which can maximize the reactive processing of dissolved contaminants and the protection of downstream water quality. Our simulations suggest denitrification dominantly occurs in riverbed hyporheic zones of streams and small rivers, whereas vertical turbulent mixing in contact with sediments dominates in mid‐size to large rivers. The metrics of connectivity and reaction significance presented here can facilitate scientifically based prioritizations of river management strategies to protect the values and functions of river corridors.     

SIMULATING THE EFFECTS OF REDUCED PRECIPITATION ON GROUND WATER AND STREAMFLOW IN THE NEBRASKA SAND HILLS1

ABSTRACT: The Nebraska Sand Hills have a unique hydrologic system with very little runoff and thick aquifers that constantly supply water to rivers, lakes, and wetlands. A ground water flow model was developed to determine the interactions between ground water and streamflow and to simulate the changes in ground water systems by reduced precipitation. The numerical modeling method includes a water balance model for the vadose zone and MOD‐FLOW for the saturated zone. The modeling results indicated that, between 1979 and 1990, 13 percent of the annual precipitation recharged to the aquifer and annual ground water loss by evapotranspiration (ET) was only about one‐fourth of this recharge. Ground water discharge to rivers accounts for about 96 percent of the streamflow in the Dismal and Middle Loup rivers. When precipitation decreased by half the average amount of the 1979 to 1990 period, the average decline of water table over the study area was 0.89 m, and the streamflow was about 87 percent of the present rate. This decline of the water table results in significant reductions in ET directly from ground water and so a significant portion of the streamflow is maintained by capture of the salvaged ET.     

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.     

Estimation of Evapotranspiration Across the Conterminous United States Using a Regression With Climate and Land‐Cover Data1

Sanford, Ward E. and David L. Selnick, 2012. Estimation of Evapotranspiration Across the Conterminous United States Using a Regression with Climate and Land‐Cover Data. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12010 Abstract: Evapotranspiration (ET) is an important quantity for water resource managers to know because it often represents the largest sink for precipitation (P) arriving at the land surface. In order to estimate actual ET across the conterminous United States (U.S.) in this study, a water‐balance method was combined with a climate and land‐cover regression equation. Precipitation and streamflow records were compiled for 838 watersheds for 1971‐2000 across the U.S. to obtain long‐term estimates of actual ET. A regression equation was developed that related the ratio ET/P to climate and land‐cover variables within those watersheds. Precipitation and temperatures were used from the PRISM climate dataset, and land‐cover data were used from the USGS National Land Cover Dataset. Results indicate that ET can be predicted relatively well at a watershed or county scale with readily available climate variables alone, and that land‐cover data can also improve those predictions. Using the climate and land‐cover data at an 800‐m scale and then averaging to the county scale, maps were produced showing estimates of ET and ET/P for the entire conterminous U.S. Using the regression equation, such maps could also be made for more detailed state coverages, or for other areas of the world where climate and land‐cover data are plentiful.     

Impact of 20 Years of Land-Cover Change on the Hydrology of Streams in the Southeastern United States1

Cruise, James F., Charles A. Laymon, and Osama Z. Al-Hamdan, 2010. Impact of 20 Years of Land-Cover Change on the Hydrology of Streams in the Southeastern United States. Journal of the American Water Resources Association (JAWRA) 46(6):1159–1170. DOI: 10.1111/j.1752-1688.2010.00483.x Abstract: Land-cover changes for portions of Alabama, Georgia, and Tennessee were estimated for the years 1980, 1990, and 2000 using classified Landsat images, and associated with hydrologic indices for 12 watersheds in the region. Rainfall-adjusted mean annual streamflow, an ET proxy (precipitation minus runoff), frequency of inundation above thresholds, and duration of inundation were used to characterize the hydrologic response of the test basins over the two-decade study period. Results indicate that several of the watersheds had undergone significant (>20%) reductions in agricultural land cover with a coincident increase in forested land. Attempts to correlate the hydrologic results with the land-cover changes were only partially successful. Watersheds with the largest land conversion from agriculture to forest (20% or more) did show significant trends in hydrologic indices indicating decreasing streamflow; however, other basins evidenced ambiguous results. The net conclusions of the study are that land-cover effects on hydrologic variables may be nuanced and can sometimes be only indirectly evident, and that a rigorous and detailed land-cover classification effort along with a battery of statistical tests with the same objective may be necessary to uncover these effects.