The effects of topographic depressions on multiscale overland flow connectivity: A high‐resolution spatiotemporal pattern analysis approach based on connectivity statistics

In watershed modelling, the traditional practice of arbitrarily filling topographic depressions in digital elevation models has raised concerns. Advanced high‐resolution remote sensing techniques, including airborne scanning laser altimetry, can identify naturally occurring depressions that impact overland flow. In this study, we used an ensemble physical and statistical modelling approach, including a 2D hydraulic model and two‐point connectivity statistics, to quantify the effects of depressions on high‐resolution overland flow patterns across spatial scales and their temporal variations in single storm events. Computations for both models were implemented using graphic processing unit‐accelerated computing. The changes in connectivity statistics for overland flow patterns between airborne scanning laser altimetry‐derived digital elevation models with (original) and without (filled) depressions were used to represent the shifts of overland flow response to depressions. The results show that depressions can either decrease or increase (to a lesser degree and shorter duration) the probability that any two points (grid locations) are hydraulically connected by overland flow pathways. We used macro‐connectivity states (Φ) as a watershed‐specific indicator to describe the spatiotemporal thresholds of connectivity variability caused by depressions. Four states of Φ are identified in a studied watershed, and each state represents different magnitudes of connectivity and connectivity changes (caused by depressions). The magnitude of connectivity variability corresponds to the states of Φ, which depend on the topological relationship between depressions, the rising/recession limb, and the total rainfall amount in a storm event. In addition, spatial distributions of connectivity variability correlate with the density of depression locations and their physical structures, which cause changes in streamflow discharge magnitude. Therefore, this study suggests that depressions are “nontrivial” in watershed modelling, and their impacts on overland flow should not be neglected. Connectivity statistics at different spatial scales and time points within a watershed provide new insights for characterizing the distributed and accumulated effects of depressions on overland flow.     

Surface and subsurface water contributions to streamflow from a mesoscale watershed in complex mountain terrain

An understanding of surface and subsurface water contributions to streamflow is essential for accurate predictions of water supply from mountain watersheds that often serve as water towers for downstream communities. As such, this study used the end‐member mixing analysis technique to investigate source water contributions and hydrologic flow paths of the 264 km² Boulder Creek Watershed, which drains the Colorado Front Range, USA. Four conservative hydrochemical tracers were used to describe this watershed as a 3 end‐member system, and tracer concentration reconstruction suggested that the application of end‐member mixing analysis was robust. On average from 2009 to 2011, snowmelt and rainwater from the subalpine zone and groundwater sampled from the upper montane zone contributed 54%, 22%, and 24% of the annual streamflow, respectively. These values demonstrate increased rainwater and decreased snow water contributions to streamflow relative to area‐weighted mean values derived from previous work at the headwater scale. Young water (2.3 ± 0.8 months) fractions of streamflow decreased from 18–22% in the alpine catchment to 8–10% in the lower elevation catchments and the watershed outlet with implications for subsurface storage and hydrological connectivity. These results contribute to a process‐based understanding of the seasonal source water composition of a mesoscale watershed that can be used to extrapolate headwater streamflow generation predictions to larger spatial scales.     

Factors controlling nitrogen and dissolved organic carbon exports across timescales in two watersheds with different land uses

Investigating factors controlling the temporal patterns of nitrogen (N) and dissolved organic carbon (DOC) exports on the basis of a comparative study of different land uses is beneficial for managing water resources, especially in agricultural watersheds. We focused our research on an agricultural watershed (AW) and a forested watershed (FW) located in the Shibetsu watershed of eastern Hokkaido, Japan, to investigate the temporal patterns of N and DOC exports and factors controlling those patterns at different timescales (inter‐annual, seasonal, and hydrological event scales). Results showed that the annual patterns of N and DOC exports significantly varied over time and were probably controlled by climate. Higher discharge volumes in 2003, a wet year, showed higher N and DOC loadings in both watersheds. However, this process was also regulated by land use associated with N inputs. Higher concentrations and loadings were shown in the agricultural watershed. At the seasonal scale, N and DOC exports in the AW and the FW were more likely controlled by sources associated with land use. The Total N (TN) and Nitrate‐N (NO₃^?‐N) had higher concentrations during snowmelt season in the AW, which may be attributed to manure application in late autumn or early winter in the agricultural watershed. Concentrations of TN, NO₃^?‐N, dissolved organic nitrogen (DON), and DOC showed higher values during the summer rainy season in the FW, related to higher litter decomposition during summer and autumn and the fertilizer application in the agricultural area during summer. Higher DOC concentrations and loadings were observed during the rainy season in the AW, which is probably attributed to higher DOC production related to temperature and microbial activity during summer and autumn in grasslands. Correlations between discharge and concentrations differed during different periods or in different watersheds, suggesting that weather discharge can adequately represent the fact that N export depends on N concentrations, discharge level, and other factors. The differing correlations between N/DOC concentrations and the Si concentration indicated that the N/DOC exports might occur along different flow paths during different periods. During baseflow, the high NO₃^?‐N exports were probably derived from deep groundwater and might have percolated from uplands during hydrological events. During hydrological events, NO₃^?‐N exports may occur along near‐surface flow paths and in deep groundwater, whereas DOC exports could be related to near‐surface flow paths. At the event scale, the relationships between discharge and concentrations of N and DOC were regulated by antecedent soil moisture (shallow groundwater condition) in each watershed. These results indicated that factors controlling N and DOC exports varied at different timescales in the Shibetsu area and that better management of manure application during winter in agricultural lands is urgently needed to control water pollution in streams. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling the interaction of climate,forest ecosystem,and hydrology to estimate catchment dissolved organic carbon export

The potential for increased loads of dissolved organic carbon (DOC) in streams and rivers is a concern for regulating the water quality in water supply watersheds. With increasing hydroclimatic variability related to global warming and shifts in forest ecosystem community and structure, understanding and predicting the magnitude and variability of watershed supply and transport of DOC over multiple time scales have become important research and management goals. In this study, we use a distributed process‐based ecohydrological model (Regional Hydro‐Ecological Simulation System [RHESSys]) to explore controls and predict streamflow DOC loads in Biscuit Brook. Biscuit Brook is a forested headwater catchment of the Neversink Reservoir, part of the New York City water supply system in the Catskill Mountains. Three different model structures of RHESSys were proposed to explore and evaluate hypotheses addressing how vegetation phenology and hydrologic connectivity between deep groundwater and riparian zones influence streamflow and DOC loads. Model results showed that incorporating dynamic phenology improved model agreement with measured streamflow in spring, summer, and fall and fall DOC concentration, compared with a static phenology. Additionally, the connectivity of deep groundwater flux through riparian zones with dynamic phenology improved streamflow and DOC flux in low flow conditions. Therefore, this study suggests the importance of inter‐annual vegetation phenology and the connectivity of deep groundwater drainage through riparian zones in the hydrology and stream DOC loading in this forested watershed and the ability of process‐based ecohydrological models to simulate these dynamics. The advantage of a process‐based modelling approach is specifically seen in the sensitivity to forest ecosystem dynamics and the interactions of hydroclimate variability with ecosystem processes controlling the supply and distribution of DOC. These models will be useful to evaluate different forest management approaches toward mitigating water quality concerns.     

Watershed structural influences on the distributions of stream network water and solute travel times under baseflow conditions

Watershed structure influences the timing, magnitude, and spatial location of water and solute entry to stream networks. In turn, stream reach transport velocities and stream network geometry (travel distances) further influence the timing of export from watersheds. Here, we examine how watershed and stream network organization can affect travel times of water from delivery to the stream network to arrival at the watershed outlet. We analysed watershed structure and network geometry and quantified the relationship between stream discharge and solute velocity across six study watersheds (11.4 to 62.8 km²) located in the Sawtooth Mountains of central Idaho, USA. Based on these analyses, we developed stream network travel time functions for each watershed. We found that watershed structure, stream network geometry, and the variable magnitude of inputs across the network can have a pronounced affect on water travel distances and velocities within a stream network. Accordingly, a sample taken at the watershed outlet is composed of water and solutes sourced from across the watershed that experienced a range of travel times in the stream network. We suggest that understanding and quantifying stream network travel time distributions are valuable for deconvolving signals observed at watershed outlets into their spatial and temporal sources, and separating terrestrial and in‐channel hydrological, biogeochemical, and ecological influences on in‐stream observations. Copyright © 2016 John Wiley & Sons, Ltd.     

An analytical solution for rapidly predicting post-fire peak streamflow for small watersheds in southern California

Following wildfires, the probability of flooding and debris flows increase, posing risks to human lives, downstream communities, infrastructure, and ecosystems. In southern California (USA), the Rowe, Countryman, and Storey (RCS) 1949 methodology is an empirical method that is used to rapidly estimate post-fire peak streamflow. We re-evaluated the accuracy of RCS for 33 watersheds under current conditions. Pre-fire peak streamflow prediction performance was low, where the average R² was 0.29 and average RMSE was 1.10 cms/km² for the 2- and 10-year recurrence interval events, respectively. Post-fire, RCS performance was also low, with an average R² of 0.26 and RMSE of 15.77 cms/km² for the 2- and 10-year events. We demonstrated that RCS overgeneralizes watershed processes and does not adequately represent the spatial and temporal variability in systems affected by wildfire and extreme weather events and often underpredicted peak streamflow without sediment bulking factors. A novel application of machine learning was used to identify critical watershed characteristics including local physiography, land cover, geology, slope, aspect, rainfall intensity, and soil burn severity, resulting in two random forest models with 45 and five parameters (RF-45 and RF-5, respectively) to predict post-fire peak streamflow. RF-45 and RF-5 performed better than the RCS method; however, they demonstrated the importance and reliance on data availability. The important parameters identified by the machine learning techniques were used to create a three-dimensional polynomial function to calculate post-fire peak streamflow in small catchments in southern California during the first year after fire (R² = 0.82; RMSE = 6.59 cms/km²) which can be used as an interim tool by post-fire risk assessment teams. We conclude that a significant increase in data collection of high temporal and spatial resolution rainfall intensity, streamflow, and sediment loading in channels will help to guide future model development to quantify post-fire flood risk.     

Fluvial transport of suspended sediment and organic carbon during flood events in a large agricultural catchment in southwest France

Water draining from a large agricultural catchment of 1 110 km² in southwest France was sampled over an 18‐month period to determine the temporal variability in suspended sediment (SS) and dissolved (DOC) and particulate organic carbon (POC) transport during flood events, with quantification of fluxes and controlling factors, and to analyze the relationships between discharge and SS, DOC and POC. A total of 15 flood events were analyzed, providing extensive data on SS, POC and DOC during floods. There was high variability in SS, POC and DOC transport during different seasonal floods, with SS varying by event from 513 to 41 750 t; POC from 12 to 748 t and DOC from 9 to 218 t. Overall, 76 and 62% of total fluxes of POC and DOC occurred within 22% of the study period. POC and DOC export from the Save catchment amounted to 3090 t and 1240 t, equivalent to 1·8 t km^?2 y^?1 and 0·7 t km^?2 y^?1, respectively. Statistical analyses showed that total precipitation, flood discharge and total water yield were the major factors controlling SS, POC and DOC transport from the catchment. The relationships between SS, POC and DOC and discharge over temporal flood events resulted in different hysteresis patterns, which were used to deduce dissolved and particulate origins. In both clockwise and anticlockwise hysteresis, POC mainly followed the same patterns as discharge and SS. The DOC‐discharge relationship was mainly characterized by alternating clockwise and anticlockwise hysteresis due to dilution effects of water originating from different sources in the whole catchment. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulated water budget of a small forested watershed in the continental/maritime hydroclimatic region of the United States

Annual streamflows have decreased across mountain watersheds in the Pacific Northwest of the United States over the last ~70 years; however, in some watersheds, observed annual flows have increased. Physically based models are useful tools to reveal the combined effects of climate and vegetation on long‐term water balances by explicitly simulating the internal watershed hydrological fluxes that affect discharge. We used the physically based Simultaneous Heat and Water (SHAW) model to simulate the inter‐annual hydrological dynamics of a 4 km² watershed in northern Idaho. The model simulates seasonal and annual water balance components including evaporation, transpiration, storage changes, deep drainage, and trends in streamflow. Independent measurements were used to parameterize the model, including forest transpiration, stomatal feedback to vapour pressure, forest properties (height, leaf area index, and biomass), soil properties, soil moisture, snow depth, and snow water equivalent. No calibrations were applied to fit the simulated streamflow to observations. The model reasonably simulated the annual runoff variations during the evaluation period from water year 2004 to 2009, which verified the ability of SHAW to simulate the water budget in this small watershed. The simulations indicated that inter‐annual variations in streamflow were driven by variations in precipitation and soil water storage. One key parameterization issue was leaf area index, which strongly influenced interception across the catchment. This approach appears promising to help elucidate the mechanisms responsible for hydrological trends and variations resulting from climate and vegetation changes on small watersheds in the region. Copyright © 2015 John Wiley & Sons, Ltd.     

Seasonal controls on DOC dynamics in nested upland catchments in NE Scotland

Spatial and temporal variability of hydrological responses affecting surface water dissolved organic carbon (DOC) concentrations are important for determining upscaling patterns of DOC export within larger catchments. Annual and intra‐annual variations in DOC concentrations and fluxes were assessed over 2 years at 12 sites (3·40–1837 km²) within the River Dee basin in NE Scotland. Mean annual DOC fluxes, primarily correlated with catchment soil coverage, ranged from 3·41 to 9·48 g m^?2 yr^?1. Periods of seasonal (summer–autumn and winter–spring) DOC concentrations (production) were delineated and related to discharge. Although antecedent temperature mainly determined the timing of switchover between periods of high DOC in the summer‐autumn and low DOC in winter‐spring, inter‐annual variability of export within the same season was largely dependent on its associated water flux. DOC fluxes ranged from 1·39 to 4·80 g m^?2 season^?1 during summer–autumn and 1·43 to 4·15 g m^?2 season^?1 in winter–spring.Relationships between DOC areal fluxes and catchment scale indicated that mainstem fluxes reflect the averaging of highly heterogeneous inputs from contrasting headwater catchments, leading to convergent DOC fluxes at catchment sizes of ca 100 km². However, during summer–autumn periods, in contrast to winter–spring, longitudinal mainstem DOC fluxes continue to decrease, most likely because of increasing biological processes. This highlights the importance of considering seasonal as well as annual changes in DOC fluxes with catchment scale. This study increases our understanding of the temporal variability of DOC upscaling patterns reflecting cumulative changes across different catchment scales and aids modelling of carbon budgets at different stages of riverine systems. Copyright © 2010 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Dissolving the mystery of subsurface controls on snowmelt–discharge dynamics in karst mountain watersheds using hydrologic timeseries

Streamflow generation in mountain watersheds is strongly influenced by snow accumulation and melt as well as groundwater connectivity. In mountainous regions with limestone and dolomite geology, bedrock formations can host karst aquifers, which play a significant role in snowmelt–discharge dynamics. However, mapping complex karst features and the resulting surface-groundwater exchanges at large scales remains infeasible. In this study, timeseries analysis of continuous discharge and specific conductance measurements were combined with gridded snowmelt predictions to characterize seasonal streamflow response and evaluate dominant watershed controls across 12 monitoring sites in a karstified 554 km² watershed in northern Utah, USA. Immense surface water hydrologic variability across subcatchments, years and seasons was linked to geologic controls on groundwater dynamics. Unlike many mountain watersheds, the variability between subcatchments could not be well described by typical watershed properties, including elevation or surficial geology. To fill this gap, a conceptual framework was proposed to characterize subsurface controls on snowmelt–discharge dynamics in karst mountain watersheds in terms of conduit flow direction, aquifer storage capacity and connectivity. This framework requires only readily measured surface water and climatic data from nested monitoring sites and was applied to the study watershed to demonstrate its applicability for evaluating dominant controls and climate sensitivity.     

Understanding hydrological processes in glacierized catchments: Evidence and implications of highly variable isotopic and electrical conductivity data

The spatial and temporal characterization of geochemical tracers over Alpine glacierized catchments is particularly difficult, but fundamental to quantify groundwater, glacier melt, and rain water contribution to stream runoff. In this study, we analysed the spatial and temporal variability of δ²H and electrical conductivity (EC) in various water sources during three ablation seasons in an 8.4‐km² glacierized catchment in the Italian Alps, in relation to snow cover and hydro‐meteorological conditions. Variations in the daily streamflow range due to melt‐induced runoff events were controlled by maximum daily air temperature and snow covered area in the catchment. Maximum daily streamflow decreased with increasing snow cover, and a threshold relation was found between maximum daily temperature and daily streamflow range. During melt‐induced runoff events, stream water EC decreased due to the contribution of glacier melt water to stream runoff. In this catchment, EC could be used to distinguish the contribution of subglacial flow (identified as an end member, enriched in EC) from glacier melt water to stream runoff, whereas spring water in the study area could not be considered as an end member. The isotopic composition of snow, glacier ice, and melt water was not significantly correlated with the sampling point elevation, and the spatial variability was more likely affected by postdepositional processes. The high spatial and temporal variability in the tracer signature of the end members (subglacial flow, rain water, glacier melt water, and residual winter snow), together with small daily variability in stream water δ²H dynamics, are problematic for the quantification of the contribution of the identified end members to stream runoff, and call for further research, possibly integrated with other natural or artificial tracers.     

Contribution of glacier meltwater to streamflow in the Wind River Range,Wyoming, inferred via a Bayesian mixing model applied to isotopic measurements

The Wind River Range (WRR) of Wyoming has the largest concentration of alpine glaciers in the American Rockies and contributes to several major river systems in the western United States. Declines in the areal extent and volume of these glaciers are well documented, and eventual loss of alpine glaciers will reduce the amount of water available for agricultural and domestic use. The contribution of glacial melt to streamflow remains largely unquantified in Wyoming. We used isotope measurements and Bayesian modeling to estimate the fractional contribution of glacier meltwater to Dinwoody Creek (DC) in the WRR on bi‐weekly and seasonal (spring, summer, and fall) time scales over 2 years. In 2007 and 2008, we made temporally intensive measurements of the stable isotope composition of water from the DC watershed. Samples of the primary sources of streamflow (snowmelt, glacier melt, rain, and baseflow) were collected during field campaigns, and automated collection of stream samples occurred over the melt season. Isotope data (D and ¹⁸O) were analyzed within a hierarchical Bayesian framework that incorporated temporal and spatial correlations. Glacial melt contributed a significant proportion (～53–59%) to streamflow in a low‐flow year (2007) or when streamflow was low during a high‐flow year (2008). In 2008, a large and persistent snowpack contributed significantly (～0·42–51%) to streamflow in mid‐summer. The large contribution of glacial melt to streamflow suggests that the loss of glaciers may impact riparian ecosystems and human water supplies in the late summer and in years with low snowpack. Copyright © 2011 John Wiley & Sons, Ltd.     

Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.     

Regional patterns of postwildfire streamflow response in the Western United States: The importance of scale‐specific connectivity

Wildfires can impact streamflow by modifying net precipitation, infiltration, evapotranspiration, snowmelt, and hillslope run‐off pathways. Regional differences in fire trends and postwildfire streamflow responses across the conterminous United States have spurred concerns about the impact on streamflow in forests that serve as water resource areas. This is notably the case for the Western United States, where fire activity and burn severity have increased in conjunction with climate change and increased forest density due to human fire suppression. In this review, we discuss the effects of wildfire on hydrological processes with a special focus on regional differences in postwildfire streamflow responses in forests. Postwildfire peak flows and annual water yields are generally higher in regions with a Mediterranean or semi‐arid climate (Southern California and the Southwest) compared to the highlands (Rocky Mountains and the Pacific Northwest), where fire‐induced changes in hydraulic connectivity along the hillslope results in the delivery of more water, more rapidly to streams. No clear streamflow response patterns have been identified in the humid subtropical Southeastern United States, where most fires are prescribed fires with a low burn severity, and more research is needed in that region. Improved assessment of postwildfire streamflow relies on quantitative spatial knowledge of landscape variables such as prestorm soil moisture, burn severity and correlations with soil surface sealing, water repellency, and ash deposition. The latest studies furthermore emphasize that understanding the effects of hydrological processes on postwildfire dynamic hydraulic connectivity, notably at the hillslope and watershed scales, and the relationship between overlapping disturbances including those other than wildfire is necessary for the development of risk assessment tools.     

Estimation of potential evapotranspiration in the mountainous Panama Canal watershed

Spatially distributed hydrometeorological and plant information within the mountainous tropical Panama Canal watershed is used to estimate parameters of the Penman–Monteith evapotranspiration formulation. Hydrometeorological data from a few surface climate stations located at low elevations in the watershed are complemented by (a) typical wet‐ and dry‐season fields of temperature, wind, water vapour and pressure produced by a mesoscale atmospheric model with a 3 × 3 km² spatial and hourly temporal resolution, and (b) leaf area index fields estimated over the watershed during a few years using satellite data with two different spatial and temporal resolutions. The mesoscale model estimates of spatially distributed surface hydrometeorological variables provide the basis for the extrapolation of the surface climate station data to produce input for the Penman–Monteith equation. The satellite information and existing digital spatial databases of land use and land cover form the basis for the estimation of Penman–Monteith spatially distributed parameter values. Spatially distributed 3 × 3 km² potential evapotranspiration estimates are obtained for the 3300 km² Panama Canal watershed. Estimates for Gatun Lake within the watershed are found to reproduce well the monthly and annual lake evaporation obtained from submerged pans. Sensitivity analysis results of potential evapotranspiration estimates with respect to cloud cover, dew formation, leaf area index distribution and mesoscale model estimates of surface climate are presented and discussed. The main conclusion is that even the limited spatially distributed hydrometeorological and plant information used in this study contributes significantly toward explaining the substantial spatial variability of potential evapotranspiration in the watershed. These results also allow the determination of key locations within the watershed where additional surface stations may be profitably placed. Copyright © 2006 John Wiley & Sons, Ltd.     

Relationship of topography to surface water chemistry with particular focus on nitrogen and organic carbon solutes within a forested watershed in Hokkaido,Japan

We studied the relationships between streamwater chemistry and the topography of subcatchments in the Dorokawa watershed in Hokkaido Island, northern Japan, to examine the use of topography as a predictor of streamwater chemistry in a watershed with relatively moderate terrain compared with other regions of Japan. Topographic characteristics of the Dorokawa watershed and its subcatchments were expressed as topographic index (TI) values, which ranged from 4·5 to 20·4 for individual grid cells (50 × 50 m²), but averaged from 6·4 to 7·4 for the 20 subcatchments. Streamwater samples for chemical analyses were collected four times between June and October 2002 from 20 locations in the watershed. The pH of water that passed through the watershed increased from ～5·0 to 7·0, with major increases in Na⁺ and Ca²⁺ and marked decreases in NO₃^? and SO . Distinctive spatial patterns were observed for dissolved organic carbon (DOC), dissolved organic nitrogen (DON), and NO₃^? concentrations of streamwater across the watershed. Statistical analyses indicated significant linear relationships between the average TI values of subcatchments and DOC, DON, and NO₃^? concentrations. Furthermore, the proportion of DOC in streamwaters in the wet season increased with TI values relative to other nitrogen species, whereas NO₃^? concentrations decreased with TI. The gradients of soil wetness and the presence of wetlands explained many of the observed spatial and temporal patterns of DOC, DON, and NO₃^? concentrations in the surface waters of the Dorokawa watershed. Our results suggest that the TI is especially useful for predicting the spatial distribution of DOC, DON and NO₃^? in the surface waters of Hokkaido, where topographical relief is moderate and wetlands more common than in other regions of Japan. Copyright © 2006 John Wiley & Sons, Ltd.     

Response characteristics of DOC flushing in an alpine catchment

The spatial distribution of source areas and associated residence times of water in the catchment are significant factors controlling the annual cycles of dissolved organic carbon (DOC) concentration in Deer Creek (Summit County, Colorado). During spring snowmelt (April–August 1992), stream DOC concentrations increased with the rising limb of the hydrograph, peaked before maximum discharge, then declined rapidly as melting continued. We investigated catchment sources of DOC to streamflow, measuring DOC in tension lysimeters, groundwater wells, snow and streamflow. Lysimeter data indicate that near-surface soil horizons are a primary contributor of DOC to streamflow during spring snowmelt. Concentrations of DOC in the lysimeters decrease rapidly during the melt period, supporting the hypothesis that hydrological flushing of catchment soils is the primary mechanism affecting the temporal variation of DOC in Deer Creek. Time constants of DOC flushing, characterizing the exponential decay of DOC concentration in the upper soil horizon, ranged from 10 to 30 days for the 10 lysimeter sites. Differences in the rate of flushing are influenced by topographical position, with near-stream riparian soils flushed more quickly than soils located further upslope. Variation in the amount of distribution of accumulated snow, and asynchronous melting of the snowpack across the landscape, staggered the onset of the spring flush throughout the catchment, prolonging the period of increased concentrations of DOC in the stream. Streamflow integrates the catchment-scale flushing responses, yielding a time constant associated with the recession of DOC in the stream channel (84 days) that is significantly longer than the time constants observed for particular locations in the upper soil. © 1997 John Wiley & Sons, Ltd.     

Patuxent Landscape Model. III. Model calibration

Using the LHEM/SME the Patuxent Landscape Model (PLM) was built to simulate fundamental ecological processes in the watershed scale driven by temporal (nutrient loadings, climatic conditions) and spatial (land use patterns) forcings. The model addresses the effects of both the magnitude and spatial patterns of land use change and agricultural practices on hydrology, plant productivity, and nutrient cycling in the landscape. The spatial resolution for the full Patuxent watershed is 1 km², while subwatersheds are analyzed at a 200 × 200 m resolution to allow adequate depiction of the pattern of ecosystems and human settlement on the landscape. The temporal resolution is different for various components of the model. We used a modular, multiscale approach to calibrate and test the model. Model results show good agreement with data.     

Spatial and Temporal Distribution of Dust‐Bound Trace Metal Elements around Kuhdasht Watershed Area in Iran

Dust, as a source of trace metal elements, affects the health of society. The spatial and temporal concentrations of dust‐bound trace metals (Cd, Pb, Ni, Zn, Cu, and Mn) in Kuhdasht watershed (456 km²), Lorestan Province, Iran, is investigated. Dust is collected using glass traps placed in ten research stations in the region. The spatial and temporal distribution of dust trace metals are plotted using ARC‐GIS. The highest and the lowest concentrations of Zn (9751150 mg kg^?1), Pb (46.352.9 mg kg^?1), and Cd (2.443.30 mg kg^?1) are obtained in winter, of Ni (98110 mg kg^?1) and Cu in autumn (16.053.5 mg kg^?1), and of Mn in summer (385505 mg kg^?1). The spatial concentrations of dust‐bound trace metals indicate all, except Cu, show a decreasing trend from the mountains toward the plains, similar to that of soil and of dust, except for Zn, which shows higher concentrations in dust than in soil. The potential sources of dust‐bound trace metals and their rate of contamination are also investigated using the enrichment and contamination factors. The major sources of Cd and Zn in the dust of watershed are due to anthropogenic activities or from activities outside the borders.