Medium-sized Watersheds Research Articles

Soil and Water Assessment Tool model predictions of annual maximum pesticide concentrations in high vulnerability watersheds.

Recent national regulatory assessments of potential pesticide exposure of threatened and endangered species in aquatic habitats have led to increased need for watershed-scale predictions of pesticide concentrations in flowing water bodies. This study was conducted to assess the ability of the uncalibrated Soil and Water Assessment Tool (SWAT) to predict annual maximum pesticide concentrations in the flowing water bodies of highly vulnerable small- to medium-sized watersheds. The SWAT was applied to 27 watersheds, largely within the midwest corn belt of the United States, ranging from 20 to 386 km2 , and evaluated using consistent input data sets and an uncalibrated parameterization approach. The watersheds were selected from the Atrazine Ecological Exposure Monitoring Program and the Heidelberg Tributary Loading Program, both of which contain high temporal resolution atrazine sampling data from watersheds with exceptionally high vulnerability to atrazine exposure. The model performance was assessed based upon predictions of annual maximum atrazine concentrations in 1-d and 60-d durations, predictions critical in pesticide-threatened and endangered species risk assessments when evaluating potential acute and chronic exposure to aquatic organisms. The simulation results showed that for nearly half of the watersheds simulated, the uncalibrated SWAT model was able to predict annual maximum pesticide concentrations within a narrow range of uncertainty resulting from atrazine application timing patterns. An uncalibrated model's predictive performance is essential for the assessment of pesticide exposure in flowing water bodies, the majority of which have insufficient monitoring data for direct calibration, even in data-rich countries. In situations in which SWAT over- or underpredicted the annual maximum concentrations, the magnitude of the over- or underprediction was commonly less than a factor of 2, indicating that the model and uncalibrated parameterization approach provide a capable method for predicting the aquatic exposure required to support pesticide regulatory decision making. Integr Environ Assess Manag 2018;14:358-368. © 2017 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals, Inc. on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Performance of a Distributed Hydrological Model Based on Soil and Moisture Zone Maps

ABSTRACT The scarcity of field data to develop soil maps through a pedological survey is one of the main limitations to using distributed hydrological models, especially in small and medium-sized watersheds. The aim of this study was to compare a distributed hydrological model prediction to a soil map based on a pedological survey and a moisture zone map obtained using the Height Above the Nearest Drainage – HAND model. The Distributed Hydrology Soil Vegetation Model – DHSVM, which is a [...]

Comparison of mercury mass loading in streams to atmospheric deposition in watersheds of Western North America: Evidence for non-atmospheric mercury sources

Annual stream loads of mercury (Hg) and inputs of wet and dry atmospheric Hg deposition to the landscape were investigated in watersheds of the Western United States and the Canadian-Alaskan Arctic. Mercury concentration and discharge data from flow gauging stations were used to compute annual mass loads with regression models. Measured wet and modeled dry deposition were compared to annual stream loads to compute ratios of Hg stream load to total Hg atmospheric deposition. Watershed land uses or cover included mining, undeveloped, urbanized, and mixed. Of 27 watersheds that were investigated, 15 had some degree of mining, either of Hg or precious metals (gold or silver), where Hg was used in the amalgamation process. Stream loads in excess of annual Hg atmospheric deposition (ratio>1) were observed in watersheds containing Hg mines and in relatively small and medium-sized watersheds with gold or silver mines, however, larger watersheds containing gold or silver mines, some of which also contain large dams that trap sediment, were sometimes associated with lower load ratios (<0.2). In the non-Arctic regions, watersheds with natural vegetation tended to have low ratios of stream load to Hg deposition (<0.1), whereas urbanized areas had higher ratios (0.34–1.0) because of impervious surfaces. This indicated that, in ecosystems with natural vegetation, Hg is retained in the soil and may be transported subsequently to streams as a result of erosion or in association with dissolved organic carbon. Arctic watersheds (Mackenzie and Yukon Rivers) had a relatively elevated ratio of stream load to atmospheric deposition (0.27 and 0.74), possibly because of melting glaciers or permafrost releasing previously stored Hg to the streams. Overall, our research highlights the important role of watershed characteristics in determining whether a landscape is a net source of Hg or a net sink of atmospheric Hg.

Modeling Surface Runoff and Evapotranspiration using SWAT and BEACH for a Tropical Watershed in North Vietnam, Compared to MODIS Products

Accurate estimation of surface runoff (Q) and evapotranspiration (ET) is a challenging task but an important research topic because both Q and ET play vital roles in the study of the hydrological cycle, of climate change, water resources, flood management and so on. In this paper we will present the modeling method to estimate the daily Q and ET for a medium-sized watershed in the tropical region of the North of Vietnam using the Soil and Water Assessment Tool (SWAT) and Bridging Event and Continuous Hydrological (BEACH) models. The models were calibrated and validated for the river discharge for SWAT and evaporation (E) for BEACH in a 12 year period from 2001 to 2012. The simulated ETs by the models were compared with the satellite-based ET of MODIS products. Our simulation results show that the SWAT and BEACH models are capable of satisfactorily reproducing (with the NSE > 0.62 and R 2 > 0.78) the stream-gauged river discharge and the observed E, respectively. Daily ET varied from 0.3 to 14 mm day -1 and was highest from May to August and lowest from December to March. Although the monthly and yearly MODIS ETs were slightly higher than those of SWAT and BEACH, a strong relationship between them was found with a standard deviation ranging from three to 40 mm. A light decrease of ET values in the 12 years can be seen in the result analyses; however a longer simulation time might be needed to ensure this trend.

A generalized concentration curve (GCC) method for storm flow hydrograph prediction in a conceptual linear reservoir-channel cascade

The unit hydrograph (UH) is a hydrological tool that represents the unit response of a watershed to a unit input of rainfall. UH models based on lumped reservoir and channel conceptual cascade assume that rainfall is evenly distributed, thus limiting the use of UHs to relatively small watersheds of less than around 500 km2 in area. In this paper, a new hydrograph prediction method, named the generalized concentration curve (GCC), was derived that can be applied to large heterogeneous watersheds. The GCC method divides the watershed into subareas by isochrones. In each subarea, an independent linear reservoir-channel cascade model that considers both attenuation and translation is established. Comparative application of the GCC and the traditional Nash instantaneous unit hydrograph to 18 storm events from three medium-sized watersheds (727, 1,800 and 5,253 km2 in area) revealed superior performance of the GCC, with the average Nash–Sutcliffe efficiency coefficient higher by 7.66%, and the average peak discharge error lower by 4.14%. This study advances the theory of UH and expands the scope of application of UH to larger watersheds.

Assessing the effects of land cover and future climate conditions on the provision of hydrological services in a medium-sized watershed of Portugal

The separated and combined effects of land-cover scenarios and future climate on the provision of hydrological services were evaluated in Vez watershed, northern Portugal. Soil and Water Assessment Tool was calibrated against daily discharge, sediments and nitrates, with good agreements between model predictions and field observations. Four hypothetical land-cover scenarios were applied under current climate conditions (eucalyptus/pine, oak, agriculture/vine and low vegetation). A statistical downscaling of four General Circulation Models, bias-corrected with ground observations, was carried out for 2021-2040 and 2041-2060, using representative concentration pathway 4.5 scenario. Also, the combined effects of future climate conditions were evaluated under eucalyptus/pine and agriculture/vine scenario. Results for land cover revealed that eucalyptus/pine scenario reduced by 7% the annual water quantity and up to 17% in the summer period. Although climate change has only a modest effect on the reduction of the total annual discharge (-7%), the effect on the water levels during summer was more pronounced, between -15% and -38%. This study shows that climate change can affect the provision of hydrological services by reducing dry season flows and by increasing flood risks during the wet months. Regarding the combined effects, future climate may reduce the low flows, which can be aggravated with eucalyptus/pine scenario. In turn, peak flows and soil erosion can be offset. Future climate may increase soil erosion and nitrate concentration, which can be aggravated with agriculture scenario. Results moreover emphasize the need to consider both climate and land-cover impacts in adaptation and land management options at the watershed scale.

Impacts of climate change on precipitation and discharge extremes through the use of statistical downscaling approaches in a Mediterranean basin

Mediterranean region is characterized by high precipitation variability often enhanced by orography, with strong seasonality and large inter-annual fluctuations, and by high heterogeneity of terrain and land surface properties. As a consequence, catchments in this area are often prone to the occurrence of hydrometeorological extremes, including storms, floods and flash-floods. A number of climate studies focused in the Mediterranean region predict that extreme events will occur with higher intensity and frequency, thus requiring further analyses to assess their effect at the land surface, particularly in small- and medium-sized watersheds. In this study, climate and hydrologic simulations produced within the Climate Induced Changes on the Hydrology of Mediterranean Basins (CLIMB) EU FP7 research project were used to analyze how precipitation extremes propagate into discharge extremes in the Rio Mannu basin (472.5km2), located in Sardinia, Italy. The basin hydrologic response to climate forcings in a reference (1971–2000) and a future (2041–2070) period was simulated through the combined use of a set of global and regional climate models, statistical downscaling techniques, and a process based distributed hydrologic model. We analyzed and compared the distribution of annual maxima extracted from hourly and daily precipitation and peak discharge time series, simulated by the hydrologic model under climate forcing. For this aim, yearly maxima were fit by the Generalized Extreme Value (GEV) distribution using a regional approach. Next, we discussed commonality and contrasting behaviors of precipitation and discharge maxima distributions to better understand how hydrological transformations impact propagation of extremes. Finally, we show how rainfall statistical downscaling algorithms produce more reliable forcings for hydrological models than coarse climate model outputs.

Effectiveness of alternative management scenarios on the sediment load in a Mediterranean agricultural watershed

The Annualised Agricultural Non-point Source model was used to evaluate the effectiveness of different management practices to control the soil erosion and sediment load in the Carapelle watershed, a Mediterranean medium-size watershed (506 km&lt;sup&gt;2&lt;/sup&gt;) located in Apulia, Southern Italy. The model was previously calibrated and validated using five years of runoff and sediment load data measured at a monitoring station located at Ordona - Ponte dei Sauri Bridge. A total of 36 events were used to estimate the performance of the model during the period 2007-2011. The model performed well in predicting runoff, as the high values of the coefficients of efficiency and determination during the validation process showed. The peak flows predictions were satisfactory especially for the high flow events; the prediction capability of sediment load was good, even if a slight over-estimation was observed. Simulations of alternative management practices show that converting the most eroding cropland cells (13.5% of the catchment area) to no tillage would reduce soil erosion by 30%, while converting them to grass or forest would reduce soil erosion by 36.5% in both cases. A crop rotation of wheat and a forage crop can also provide an effective way for soil erosion control as it reduces erosion by 69%. Those results can provide a good comparative analysis for conservation planners to choose the best scenarios to be adopted in the watershed to achieve goals in terms of soil conservation and water quality.

Distributed hydrologic modeling of a sparsely monitored basin in Sardinia, Italy, through hydrometeorological downscaling

Abstract. The water resources and hydrologic extremes in Mediterranean basins are heavily influenced by climate variability. Modeling these watersheds is difficult due to the complex nature of the hydrologic response as well as the sparseness of hydrometeorological observations. In this work, we present a strategy to calibrate a distributed hydrologic model, known as TIN-based Real-time Integrated Basin Simulator (tRIBS), in the Rio Mannu basin (RMB), a medium-sized watershed (472.5 km2) located in an agricultural area in Sardinia, Italy. In the RMB, precipitation, streamflow and meteorological data were collected within different historical periods and at diverse temporal resolutions. We designed two statistical tools for downscaling precipitation and potential evapotranspiration data to create the hourly, high-resolution forcing for the hydrologic model from daily records. Despite the presence of several sources of uncertainty in the observations and model parameterization, the use of the disaggregated forcing led to good calibration and validation performances for the tRIBS model, when daily discharge observations were available. The methodology proposed here can be also used to disaggregate outputs of climate models and conduct high-resolution hydrologic simulations with the goal of quantifying the impacts of climate change on water resources and the frequency of hydrologic extremes within medium-sized basins.

Runoff and sediment yield modeling in a medium-size mediterranean watershed

The AnnAGNPS model was used to estimate runoff, peak discharge and sediment yield at the event scale in the Carapelle watershed, a Mediterranean medium-size watershed (506 km2) located in Apulia, Southern Italy. The model was calibrated and validated using five years of runoff and sediment yield data measured at a monitoring station located at Ordona – Ponte dei Sauri Bridge. A total of 36 events was used to estimate the output of the model during the period 2007-2011, in comparison to the corresponding observations at the watershed outlet. The model performed well in predicting runoff, as was testified by the high values of the coefficients of efficiency and determination during the validation process. The peak flows predictions were satisfactory especially for the high flow events; the prediction capability of sediment yield was good, even if a slight over-estimation was observed. Finally, the model was used to evaluate the effectiveness of different Management practices (MPs) on the watershed (converting wheat to forest, using vegetated streams, crop rotation corn/soybean, no tillage). While the maximum reduction in sediment yield was achieved converting wheat to forest, the best compromises between soil conservation and agriculture resulted to be crop rotations.

Runoff and sediment yield modeling in a medium-size mediterranean watershed

The AnnAGNPS model was used to estimate runoff, peak discharge and sediment yield at the event scale in the Carapelle watershed, a Mediterranean medium-size watershed (506 km2) located in Apulia, Southern Italy. The model was calibrated and validated using five years of runoff and sediment yield data measured at a monitoring station located at Ordona – Ponte dei Sauri Bridge. A total of 36 events was used to estimate the output of the model during the period 2007-2011, in comparison to the corresponding observations at the watershed outlet. The model performed well in predicting runoff, as was testified by the high values of the coefficients of efficiency and determination during the validation process. The peak flows predictions were satisfactory especially for the high flow events; the prediction capability of sediment yield was good, even if a slight over-estimation was observed. Finally, the model was used to evaluate the effectiveness of different Management practices (MPs) on the watershed (converting wheat to forest, using vegetated streams, crop rotation corn/soybean, no tillage). While the maximum reduction in sediment yield was achieved converting wheat to forest, the best compromises between soil conservation and agriculture resulted to be crop rotations.

Estimation of Runoff, Peak Discharge and Sediment Load at the Event Scale in a Medium‐Size Mediterranean Watershed Using the Annagnps Model

AbstractIn this paper, the Annualized Agricultural Non‐Point Source (AnnAGNPS) model has been used to estimate runoff, peak discharge and sediment load at the event scale in a Mediterranean watershed. The study area is the Carapelle torrent, Southern Italy (area = 506 km2), where continuous rainfall, streamflow and sediment load data are available. Nineteen flood events have been registered in the period 2007–2009 and were used for the application of the model.The aim of the paper is to evaluate the predictive accuracy of the model at the event scale, in a medium‐size watershed, given the specific conditions of the semi‐arid environments. A sensitivity analysis has been carried out to assign the correct parameterization: the mean normalized output variation of the most meaningful input parameters pointed out the influence of the curve number on runoff, peak discharge and sediment load predictions (values greater than 1); the MN Manning's roughness coefficient and K, C and P factors of the universal soil loss equation showed a moderate influence on sediment load simulations (values between 0·5 and 1).The selection of the Soil Conservation Service synthetic storm types has been based on the observed storm events analysis to improve the peak discharge simulations. The model prediction has proved to be good for runoff (R2 = 0·74, NSE = 0·75, W = 0·92) and peak discharge (R2 = 0·85, NSE = 0·70, W = 0·94), and satisfactory for sediment yield (R2 = 0·70, NSE = 0·63, W = 0·91). The relative error is lower for high events; this result is quite interesting in semi‐arid environments, where most of the annual sediment yield is concentrated in a few, severe events. Copyright © 2013 John Wiley &amp; Sons, Ltd.

Geospatially based distributed rainfall-runoff modelling for simulation of internal and outlet responses in a semi-forested lower Himalayan watershed

This study presents a Geographic Information System (GIS)-based distributed rainfall-runoff model for simulating surface flows in small to large watersheds during isolated storm events. The model takes into account the amount of interception storage to be filled using a modified Merriam (1960) approach before estimating infiltration by the Smith and Parlange (1978) method. The mechanics of overland and channel flow are modelled by the kinematic wave approximation of the Saint Venant equations which are then numerically solved by the weighted four-point implicit finite difference method. In this modelling the watershed was discretized into overland planes and channels using the algorithms proposed by Garbrecht and Martz (1999). The model code was first validated by comparing the model output with an analytical solution for a hypothetical plane. Then the model was tested in a medium-sized semi-forested watershed of Pathri Rao located in the Shivalik ranges of the Garhwal Himalayas, India. Initially, a local sensitivity analysis was performed to identify the parameters to which the model outputs like runoff volume, peak flow and time to peak flow are sensitive. Before going for model validation, calibration was performed using the Ordered-Physics-based Parameter Adjustment (OPPA) method. The proposed Physically Based Distributed (PBD) model was then evaluated both at the watershed outlet as well as at the internal gauging station, making this study a first of its kind in Indian watersheds. The results of performance evaluation indicate that the model has simulated the runoff hydrographs reasonably well within the watershed as well as at the watershed outlet with the same set of calibrated parameters. The model also simulates, realistically, the temporal variation of the spatial distribution of runoff over the watershed and the same has been illustrated graphically. Copyright © 2011 John Wiley & Sons, Ltd.

Validation of SCS Method for Runoff Estimation

The Soil Conservation Service (SCS) method is widely used to estimate runoff from small- to medium-sized watersheds. The most critical assumption of the SCS method is that the ratio of actual retention to potential retention is the same as the ratio of actual runoff to potential runoff; however, this assumption has not been empirically validated. The paper develops a framework to test this proportionality assumption that underpins the SCS method. Using data on rainfall intensity and storm runoff from 210 site events from Australia and Southeast Asian countries, this paper shows that a strong relationship exists between maximum retention on the basis of the SCS equation and the product of effective storm duration and spatially averaged maximum infiltration rate, and empirical support exists for the proportionality assumption for runoff estimation. Relating the maximum retention to the effective storm duration and maximum infiltration rate provides additional avenues for prediction of storm runoff amount and peak runoff rate, which are the key design parameters for storm water control and management.

Anthropogenic impacts on suspended sediment load in the Upper Yangtze river

Climate change and human disturbance drive catchment erosion and increase riverine sediment load sensitively in small and medium-sized watersheds. This is not always true in large basins, where aggregation and buffering effects have dampen the ability to determine the driving forces of sedimentation. Even though there are significant responses to sedimentation in large river basins, it is difficult to get a precise quantitative assessment of specific drivers. This paper develops a methodology to identify driving forces that change suspended sediment load in the Upper Yangtze river. Annual runoff and sediment load data from 1954 to 2005 at the Yichang gauging station in the Upper Yangtze basin, daily precipitation data from 60 meteorological stations, and survey data on reservoir sediment were collected for the study. Sediment load/rainfall erosivity (S/R), is a new proxy indicator introduced to reflect human activities. Since the mid-1980s, S/R in the Upper Yangtze has dramatically declined from 0.21 to 0.03 (x10(10) t ha h MJ(-1) mm(-1)), indicating that human activity has played a key role in the decline of the suspended sediment load. Before the mid-1980s, a higher average S/R is attributed to large-scale deforestation and land reclamation. A significant sediment decrease occurred from 1959 to 1961 during an extreme drought condition, and an increase in sedimentation in 1998 coincided with an extreme flood event, which was well recorded in the S/R curve. This indicates that the S/R proxy is able to distinguish anthropogenic from climate impacts on suspended sediment load, but is not necessarily indicatory in extreme climate events. In addition, typical drivers of riverine sediment load variation including soil conservation projects, reservoirs construction, and land use/cover change are discussed.

Performance of a distributed semi-conceptual hydrological model under tropical watershed conditions

Abstract Conceptual rainfall-runoff models at the watershed scale are useful tools for assisting in managing and planning water resources, making it possible to estimate streamflow and to predict hydrologic impacts due to land-use changes. The objective of this study was to calibrate and to validate the LAvras Simulation of Hydrology (LASH Model) (Beskow, 2009) in a Brazilian Tropical Watershed for daily streamflow. LASH is a continuous, distributed, semi-physically based model for simulation of different hydrologic components on a daily basis. The Shuffled Complex Evolution (SCE-UA) global search method was used with the LASH model in order to optimize model parameters that were found to be the most sensitive or not directly measurable. The LASH model was calibrated over a 2-year period, thereafter, the parameters obtained through the calibration were kept constant for the validation step using a different period of time from that analyzed during the calibration. The Nash–Sutcliffe coefficient (CNS) values found were 0.820 and 0.764 during calibration and validation, respectively, whereas, the CNS (log Q) values equal to 0.821 and 0.770 were obtained for the same periods. The SCE-UA method was found to be an efficient algorithm for searching ‘optimal’ model parameter values. It was possible to conclude that the model has a great potential for being applied in generating the long-term streamflow as well as flow-duration curves. Therefore, the model can reliably be applied under tropical conditions of this medium-sized watershed or other similar watersheds, thus making it useful to plan the sustainable development of similar tropical and subtropical watersheds.

A process-based, distributed hydrologic model based on a large-scale method for surface–subsurface coupling

Process-based watershed models are useful tools for understanding the impacts of natural and anthropogenic influences on water resources and for predicting water and solute fluxes exported from watersheds to receiving water bodies. The applicability of process-based hydrologic models has been previously limited to small catchments and short time frames. Computational demands, especially the solution to the three-dimensional subsurface flow domain, continue to pose significant constraints. This paper documents the mathematical development, numerical testing and the initial application of a new distributed hydrologic model PAWS (Process-based Adaptive Watershed Simulator). The model solves the governing equations for the major hydrologic processes efficiently so that large scale applications become relevant. PAWS evaluates the integrated hydrologic response of the surface–subsurface system using a novel non-iterative method that couples runoff and groundwater flow to vadose zone processes approximating the 3D Richards equation. The method is computationally efficient and produces physically consistent solutions. All flow components have been independently verified using analytical solutions and experimental data where applicable. The model is applied to a medium-sized watershed in Michigan (1169 km 2) achieving high performance metrics in terms of streamflow prediction at two gages during the calibration and verification periods. PAWS uses public databases as input and possesses full capability to interact with GIS datasets. Future papers will describe applications to other watersheds and the development and application of fate and transport modules.

Analysis of stream water temperature changes during rainfall events in forested watersheds

Despite continued interest in stream water temperature (Tw) analysis, there are few studies of Tw response to rainfall events at forested watersheds. We examined 61 sets of data on Tw for 21 rainfall events at 16 watersheds with various slope gradients (from 0.08 to 0.56) in four regions of Japan from June 2004 to December 2005. The investigation focused on the changes of specific discharge (ΔQs) and ΔTw at medium-sized watersheds (0.5–100 ha). The results clearly demonstrated different flow patterns expressed by Qs vs. Tw hysteretic loops. Those were clockwise in Period I (April–September) and counterclockwise in Period II (October–March), except for lower slope gradient at Aichi, where counterclockwise loops were observed in both periods. These differences in hysteretic loops could be explained by the differences in Tw and in response times to rainfall between surface/subsurface and groundwater flows. The response times were probably determined by the slope gradient and the vertical level of groundwater. We also found that the changes in air temperature (ΔTa) influenced ΔTw to a lesser degree than Qs. The average rainfall intensities in Period I and Period II (9.3 ± 1.7 and 5.4 ± 0.2 mm/h, respectively) affected the average values of ΔQs and ΔTw (6.62 ± 4.08 mm/h and 1.7 ± 0.4°C; 0.85 ± 0.68 mm/h and 0.9 ± 0.3°C, respectively). This indicates that slope gradient and Qs influenced ΔTw by changing the relative proportions of flow paths. In addition, the water table changes influenced the percentage of groundwater flow to the stream.

Comparative study of two models to simulate diffuse nitrogen and phosphorus pollution in a medium-sized watershed, southeast China

The aim of this study was to compare and assess two models to calculate diffuse nitrogen and phosphorus emissions in a selected watershed. The GIS-based empirical model and the physically-based AnnAGNPS model were evaluated for comparative purposes. The methodologies were applied for the Jiulong River watershed, covering 14,700 km 2, located in southeast China, with intensive agricultural activities. The calculated loadings by AnnAGNPS model was checked by the measured values at the watershed outlet, whereas the calculated nitrogen and phosphorus emission by GIS-based empirical model spatially provided the potential values in terms of sub-watersheds, districts/counties, and land use type. Both models gave similar levels of diffuse total nitrogen emissions, which also fit well with previous estimates made in the Jiulong River watershed. Comparatively, the GIS-based empirical model gave sound results of source apportionment of non-point source pollution (NPS) from the available input data and critical source areas identification of diffuse nitrogen and phosphorus pollution. The AnnAGNPS model predicted reasonable nitrogen loading at the watershed outlet and simulated well for NPS management alternatives under changing land use conditions. The study indicated that the GIS-based empirical model has its advantage in extensive studies as a decisions support tool for preliminary design since it is easily applied to large watersheds with fewer data requirements, while AnnAGNPS has its advantage in detailed emission assessment and scenario development.

Influences of forested watershed conditions on fluctuations in stream water temperature with special reference to watershed area and forest type

In order to gain insight into the effect of watershed conditions on fluctuations in stream water temperature, we statistically analyzed water temperature data for 1 year, using root mean square (Rms) and harmonic (A Amplitude, φ delay time) methods. The average values of delay time (days) between air and water temperatures (Ta and Tw) of small ( 100 ha) watersheds were 4.53 ± 0.82 days, 11.83 ± 3.88 days and 4.45 ± 1.52 days, respectively. Fluctuations in stream water temperature expressed by Rms (Rms Tw/Rms Ta) and harmonic methods (A −Tw/A −Ta) in the medium-sized watersheds with moderate slope gradients were 0.37 ± 0.09 and 0.56 ± 0.14, respectively. These values increased in the larger watersheds with low slope gradients, including five large rivers covered by various landscapes, with their averages of 0.53 ± 0.09 and 0.78 ± 0.09, respectively, indicating the influences of solar radiation and heat transfer processes. In the smaller watersheds with high slope gradients, these values were 0.73 ± 0.02 and 0.87 ± 0.03, respectively, suggesting that shorter passage time affected water temperatures. With respect to forest type, these values at badly managed hinoki forest watersheds (0.45 ± 0.04 and 0.73 ± 0.07) were larger than those at broadleaf forest (0.34 ± 0.04 and 0.51 ± 0.12) and well-managed hinoki forest (0.33 ± 0.04 and 0.51 ± 0.07) watersheds, indicating different proportions of flow paths.