Influence of Riparian Seepage Zones on Nitrate Variability in Two Agricultural Headwater Streams

Riparian seeps have been recognized for their contributions to stream flow in headwater catchments, but there is limited data on how seeps affect stream water quality. The objective of this study was to examine the effect of seeps on the variability of stream NO₃‐N concentrations in FD36 and RS, two agricultural catchments in Pennsylvania. Stream samples were collected at 10‐m intervals over reaches of 550 (FD36) and 490 m (RS) on 21 occasions between April 2009 and January 2012. Semi‐variogram analysis was used to quantify longitudinal patterns in stream NO₃‐N concentration. Seep water was collected at 14 sites in FD36 and 7 in RS, but the number of flowing seeps depended on antecedent conditions. Seep NO₃‐N concentrations were variable (0.1‐29.5 mg/l) and were often greater downslope of cropped fields compared to other land uses. During base flow, longitudinal variability in stream NO₃‐N concentrations increased as the number of flowing seeps increased. The influence of seeps on the variability of stream NO₃‐N concentrations was less during storm flow compared to the variability of base flow NO₃‐N concentrations. However, 24 h after a storm in FD36, an increase in the number of flowing seeps and decreasing streamflow resulted in the greatest longitudinal variability in stream NO₃‐N concentrations recorded. Results indicate seeps are important areas of NO₃‐N delivery to streams where targeted adoption of mitigation measures may substantially improve stream water quality.     

Long‐Term High‐Resolution Radar Rainfall Fields for Urban Hydrology

Accurate records of high‐resolution rainfall fields are essential in urban hydrology, and are lacking in many areas. We develop a high‐resolution (15 min, 1 km²) radar rainfall data set for Charlotte, North Carolina during the 2001‐2010 period using the Hydro‐NEXRAD system with radar reflectivity from the National Weather Service Weather Surveillance Radar 1988 Doppler weather radar located in Greer, South Carolina. A dense network of 71 rain gages is used for estimating and correcting radar rainfall biases. Radar rainfall estimates with daily mean field bias (MFB) correction accurately capture the spatial and temporal structure of extreme rainfall, but bias correction at finer timescales can improve cold‐season and tropical cyclone rainfall estimates. Approximately 25 rain gages are sufficient to estimate daily MFB over an area of at least 2,500 km², suggesting that robust bias correction is feasible in many urban areas. Conditional (rain‐rate dependent) bias can be removed, but at the expense of other performance criteria such as mean square error. Hydro‐NEXRAD radar rainfall estimates are also compared with the coarser resolution (hourly, 16 km²) Stage IV operational rainfall product. Stage IV is adequate for flood water balance studies but is insufficient for applications such as urban flood modeling, in which the temporal and spatial scales of relevant hydrologic processes are short. We recommend the increased use of high‐resolution radar rainfall fields in urban hydrology.     

Analysis of Daily Peaking and Run‐of‐River Operations with Flow Variability Metrics,Considering Subdaily to Seasonal Time Scales

Environmental flows are an important consideration in licensing hydropower projects as operational flow releases can result in adverse conditions for downstream ecological communities. Flow variability assessments have typically focused on pre‐ and post‐dam conditions using metrics based on daily averaged flow values. This study used subdaily and daily flow data to assess environmental flow response to changes in hydropower operations from daily peaking to run‐of‐river. An analysis tool was developed to quantify flow variability metrics and was applied to four hydropower projects. Significant differences were observed between operations at the 99% confidence level in the median flow values using hourly averaged flow datasets. Median daily rise and fall rates decreased on average 34.5 and 27.9%, respectively, whereas median hourly rise and fall rates decreased on average 50.1 and 50.6%, respectively. Differences in operational flow regimes were more pronounced in the hourly averaged flow datasets and less pronounced or nonexistent in the daily averaged flow datasets. These outcomes have implications for the development of ecology‐flow relationships that quantify effects of flow on processes such as fish stranding and displacement, along with habitat stability. Results indicate that flow variability statistics should be quantified using subdaily datasets to accurately represent the nature of hydropower operations, especially for daily peaking facilities.     

Effects of Land Use and Sample Location on Nitrate‐Stream Flow Hysteresis Descriptors during Storm Events

The U.S. Geological Survey's New Jersey and Iowa Water Science Centers deployed ultraviolet‐visible spectrophotometric sensors at water‐quality monitoring sites on the Passaic and Pompton Rivers at Two Bridges, New Jersey, on Toms River at Toms River, New Jersey, and on the North Raccoon River near Jefferson, Iowa to continuously measure in‐stream nitrate plus nitrite as nitrogen (NO₃ + NO₂) concentrations in conjunction with continuous stream flow measurements. Statistical analysis of NO₃ + NO₂ vs. stream discharge during storm events found statistically significant links between land use types and sampling site with the normalized area and rotational direction of NO₃ + NO₂‐stream discharge (N‐Q) hysteresis patterns. Statistically significant relations were also found between the normalized area of a hysteresis pattern and several flow parameters as well as the normalized area adjusted for rotational direction and minimum NO₃ + NO₂ concentrations. The mean normalized hysteresis area for forested land use was smaller than that of urban and agricultural land uses. The hysteresis rotational direction of the agricultural land use was opposite of that of the urban and undeveloped land uses. An r² of 0.81 for the relation between the minimum normalized NO₃ + NO₂ concentration during a storm vs. the normalized NO₃ + NO₂ concentration at peak flow suggested that dilution was the dominant process controlling NO₃ + NO₂ concentrations over the course of most storm events.     

Incorporating Variable Source Area Hydrology into a Spatially Distributed Direct Runoff Model1

Buchanan, Brian, Zachary M. Easton, Rebecca Schneider, and M. Todd Walter, 2011. Incorporating Variable Source Area Hydrology Into a Spatially Distributed Direct Runoff Model. Journal of the American Water Resources Association (JAWRA) 48(1): 43‐60. DOI: 10.1111/j.1752‐1688.2011.00594.x Abstract: Few hydrologic models simulate both variable source area (VSA) hydrology, and runoff‐routing at high enough spatial resolutions to capture fine‐scale hydrologic pathways connecting VSA to the stream network. This paper describes a geographic information system‐based operational model that simulates the spatio‐temporal dynamics of VSA runoff generation and distributed runoff‐routing, including through complex artificial drainage networks. The model combines the Natural Resource Conservation Service’s Curve Number (CN) equation for estimating storm runoff with the topographic index concept for predicting the locations of VSA and a runoff‐routing algorithm into a new spatially distributed direct hydrograph (SDDH) model (SDDH‐VSA). Using a small agricultural watershed in central New York, SDDH‐VSA results were compared to those from a SDDH model using the traditional land use assumptions for the CN (SDDH‐CN). The SDDH‐VSA model generally agreed better with observed discharge than the SDDH‐CN model (average, Nash‐Sutcliffe efficiency of 0.69 vs. 0.58, respectively) and resulted in more realistic spatial patterns of runoff‐generating areas. The SDDH approach did not correctly capture the timing of runoff from small storms in dry periods. Despite this type of limitation, SDDH‐VSA extends the applicability of the SDDH technique to VSA conditions, providing a basis for new tools to help identify critical management areas and assess water quality risks due to landscape alterations.     

Hydrologic Controls on Nitrogen and Phosphorous Dynamics in Relict Oxbow Wetlands Adjacent to an Urban Restored Stream

Although wetlands are known to be sinks for nitrogen (N) and phosphorus (P), their function in urban watersheds remains unclear. We analyzed water and nitrate (NO₃^?) and phosphate (PO₄^3?) dynamics during precipitation events in two oxbow wetlands that were created during geomorphic stream restoration in Baltimore County, Maryland that varied in the nature and extent of connectivity to the adjacent stream. Oxbow 1 (Ox1) received 1.6‐4.2% and Oxbow 2 (Ox2) received 4.2‐7.4% of cumulative streamflow during storm events from subsurface seepage (Ox1) and surface flow (Ox2). The retention time of incoming stormwater ranged from 0.2 to 6.7 days in Ox1 and 1.8 to 4.3 days in Ox2. Retention rates in the wetlands ranged from 0.25 to 2.74 g N/m²/day in Ox1 and 0.29 to 1.94 g N/m²/day in Ox2. Percent retention of the NO₃^?‐N load that entered the wetlands during the storm events ranged from 64 to 87% and 23 to 26%, in Ox1 and Ox2, respectively. During all four storm events, Ox1 and Ox2 were a small net source of dissolved PO₄^3? to the adjacent stream (i.e., more P exited than entered the wetland), releasing P at a rate of 0.23‐20.83 mg P/m²/day and 3.43‐24.84 mg P/m²/day, respectively. N and P removal efficiency of the oxbows were regulated by hydrologic connectivity, hydraulic loading, and retention time. Incidental oxbow wetlands have potential to receive urban stream and storm flow and to be significant N sinks, but they may be sources of P in urban watersheds.     

ESCHERICHIA COLI SURVIVAL IN MANTLED KARST SPRINGS AND STREAMS,NORTHWEST ARKANSAS OZARKS,USA1

Recent studies indicate fecal coliform bacterial concentrations, including Escherichia coli (E. coli), characteristically vary by several orders of magnitude, depending on the hydrology of storm recharge and discharge. E. coli concentrations in spring water increase rapidly during the rising limb of a storm hydrograph, peak prior to or coincident with the peak of the storm pulse, and decline rapidly, well before the recession of the storm hydrograph. This suggests E. coli are associated with resuspension of sediment during the onset of turbulent flow, and indicates viable bacteria reside within the spring and stream sediments. E. coli inoculated chambers were placed in spring and stream environments within the mantled karst of northwest Arkansas to assess long term (> 75 days) E. coli viability. During the 75‐day study, a 4‐log die‐off of E. coli was observed for chambers placed in the Illinois River, and a 5‐log die‐off for chambers placed in Copperhead Spring. Extrapolation of the regression line for each environment indicates E. coli concentration would reach 1 most probable number (MPN)/100 g sediment at Copperhead Spring in about 105 days, and about 135 days in the Illinois River, based on a starting inoculation of 2.5 × 107 MPN E. coli/100 g of sediment. These in situ observations indicate it is possible for E. coli to survive in these environments for at least four months with no fresh external inputs.     

A Daily Time Series Analysis of StreamWater Phosphorus Concentrations Along an Urban to Forest Gradient

During a 1-year period, we sampled stream water total phosphorus (TP) concentrations daily and soluble reactive phosphorus (SRP) concentrations weekly in four Seattle area streams spanning a gradient of forested to urban-dominated land cover. The objective of this study was to develop time series models describing stream water phosphorus concentration dependence on seasonal variation in stream base flows, short-term flow fluctuations, antecedent flow conditions, and rainfall. Stream water SRP concentrations varied on average by ±18% or ±5.7 μg/L from one week to another, whereas TP varied ±48% or ±32.5 μg/L from one week to another. On average, SRP constituted about 47% of TP. Stream water SRP concentrations followed a simple sine-wave annual cycle with high concentrations during the low-flow summer period and low concentrations during the high-flow winter period in three of the four study sites. These trends are probably due to seasonal variation in the relative contributions of groundwater and subsurface flows to stream flow. In forested Issaquah Creek, SRP concentrations were relatively constant throughout the year except during the fall, when a major salmon spawning run occurred in the stream and SRP concentrations increased markedly. Stream water SRP concentrations were statistically unrelated to short-term flow fluctuations, antecedent flow conditions, or rainfall in each of the study streams. Stream water TP concentrations are highly variable and strongly influenced by short-term flow fluctuations. Each of the processes assessed had statistically significant correlations with TP concentrations, with seasonal base flow being the strongest, followed by antecedent flow conditions, short-term flow fluctuations, and rainfall. Times series models for each individual stream were able to predict ∼70% of the variability in the SRP annual cycle in three of the four streams (r² = 0.57–0.81), whereas individual TP models explained ∼50% of the annual cycle in all streams (r² = 0.39–0.59). Overall, time series models for SRP and TP dynamics explained 82% and 76% of the variability for these variables, respectively. Our results indicate that SRP, the most biologically available and therefore most important phosphorus fraction, has simpler and easier-to-predict seasonal and weekly dynamics.     

RUNOFF MODELING OF A MOUNTAINOUS CATCHMENT USING TOPMODEL: A CASE STUDY1

ABSTRACT: The rainfall‐runoff response of the Tygarts Creek Catchment in eastern Kentucky is studied using TOPMODEL, a hydrologic model that simulates runoff at the catchment outlet based on the concepts of saturation excess overland flow and subsurface flow. Unlike the traditional application of this model to continuous rainfall‐runoff data, the use of TOPMOEL in single event runoff modeling, specifically floods, is explored here. TOPMODEL utilizes a topographic index as an indicator of the likely spatial distribution of rainfall excess generation in the catchment. The topographic index values within the catchment are determined using the digital terrain analysis procedures in conjunction with digital elevation model (DEM) data. Select parameters in TOPMODEL are calibrated using an iterative procedure to obtain the best‐fit runoff hydrograph. The calibrated parameters are the surface transmissivity, T_O, the transmissivity decay parameter, m, and the initial moisture deficit in the root zone, S_r0. These parameters are calibrated using three storm events and verified using three additional storm events. Overall, the calibration results obtained in this study are in general agreement with the results documented from previous studies using TOPMODEL. However, the parameter values did not perform well during the verification phase of this study.     

PREDICTING PEAK STREAM FLOW FROM AN UNDISTURBED WATERSHED IN THE CENTRAL APPALACHIANS1

ABSTRACT: Data from a small forested catchment were used to model peak stream flow as a function of basic hydrologic variables associated with 112 rain storms. Rainfall depth and initial stream flow rate accounted for 87.1 percent of peak flow variability. Forty expressions of rainfall intensity (describing both the temporal sequence of intensity for 20 equal storm intervals, and maximum intensity for 20 separate interval lengths) were used in an attempt to improve the predictability of basic models. None of the intensity parameters improved predictability by as much as 2 percent, apparently because the most intense rainfall bursts generally occurred near the beginning of storm periods. Mean rainfall intensity for entire storms was generally as effective as any of the shorter interval intensities, and its use helped to linearize the relationship between peak flow and rainfall depth and duration.     

Impacts of land cover on stream hydrology in the West Georgia Piedmont, USA

The southeastern United States is experiencing rapid urban development. Consequently, Georgia's streams are experiencing hydrologic alterations from extensive development and from other land use activities such as livestock grazing and silviculture. A study was performed to assess stream hydrology within 18 watersheds ranging from 500 to 2500 ha. Study streams were first, second, or third order and hydrology was continuously monitored from 29 July 2003 to 23 September 2004 using InSitu pressure transducers. Rating curves between stream stage (i.e., water depth) and discharge were developed for each stream by correlating biweekly discharge measurements and stage data. Dependent variables were calculated from discharge data and placed into 4 categories: flow frequency (i.e., the number of times a predetermined discharge threshold is exceeded), flow magnitude (i.e., maximum and minimum flows), flow duration (i.e., the amount of time discharge was above or below a predetermined threshold), and flow predictability and flashiness. Fine resolution data (i.e., 15-min interval) were also compared to daily discharge data to determine if resolution affected how streams were classified hydrologically. Urban watersheds experienced flashy discharges during storm events, whereas pastoral and forested watersheds showed less flashy hydrographs. Also, in comparison to all other flow variables, flow frequency measures were most strongly correlated to land cover. Furthermore, the stream hydrology was explained similarly with both the 15-min and daily data resolutions.     

The Water‐Year Water Balance of the Colorado River Basin

Model‐estimated monthly water balance components (i.e., potential evapotranspiration, actual evapotranspiration, and runoff (R)) for 146 United States (U.S.) Geological Survey 8‐digit hydrologic units located in the Colorado River Basin (CRB) are used to examine the temporal and spatial variability of the CRB water balance for water years 1901 through 2014 (a water year is the period from October 1 of one year through September 30 of the following year). Results indicate that the CRB can be divided into six subregions with similar temporal variability in monthly R. The water balance analyses indicated that approximately 75% of total water‐year R is generated by just one CRB subregion and that most of the R in the basin is derived from surplus (S) water generated during the months of October through April. Furthermore, the analyses show that temporal variability in S is largely controlled by the occurrence of negative atmospheric pressure anomalies over the northwestern conterminous U.S. (CONUS) and positive atmospheric pressure anomalies over the southeastern CONUS. This combination of atmospheric pressure anomalies results in an anomalous flow of moist air from the North Pacific Ocean into the CRB, particularly the Upper CRB. Additionally, the occurrence of extreme dry and wet periods in the CRB appears to be related to variability of the Atlantic Multidecadal Oscillation and the Pacific Decadal Oscillation.     

Effect of Soils on Water Quantity and Quality in Piedmont Forested Headwater Watersheds of North Carolina1

Boggs, Johnny, Ge Sun, David Jones, and Steven G. McNulty, 2012. Effect of Soils on Water Quantity and Quality in Piedmont Forested Headwater Watersheds of North Carolina. Journal of the American Water Resources Association (JAWRA) 1‐19. DOI: 10.1111/jawr.12001 Abstract: Water quantity and quality data were compared from six headwater watersheds on two distinct soil formations, Carolina Slate Belt (CSB) and Triassic Basins (TB). CSB soils are generally thicker, less erodible, and contain less clay content than soils found in TB. TB generated significantly more discharge/precipitation ratio than CSB (0.33 vs. 0.24) in the 2009 dormant season. In the 2009 growing season, TB generated significantly less discharge/precipitation ratio than CSB (0.02 vs. 0.07). Over the entire monitoring period, differences in discharge/precipitation ratios between CSB and TB were not significantly different (0.17 vs. 0.20, respectively). Storm‐flow rates were significantly higher in TB than CSB in both dormant and growing season. Benthic macroinvertebrate biotic index scores were excellent for all streams. Nutrient concentrations and exports in CSB and TB were within background levels for forests. Low‐stream nitrate and ammonium concentrations and exports suggested that both CSB and TB were nitrogen limited. Soils appear to have had a significant influence on seasonal and storm‐flow generation, but not on long‐term total water yield and water quality under forested conditions. This study indicated that watersheds on TB soils might be more prone to storm‐flow generation than on CSB soils when converted from forest to urban. Future urban growth in the area should consider differences in baseline hydrology and effects of landuse change on water quantity and quality.     

IMPACTS OF CALIFORNIA'S CLIMATIC REGIMES AND COASTAL LAND USE CHANGE ON STREAMFLOW CHARACTERISTICS1

ABSTRACT: To investigate the impacts of urbanization and climatic fluctuations on stream flow magnitude and variability in a Mediterranean climate, the HEC‐HMS rainfall/runoff model is used to simulate stream flow for a 14‐year period (October 1, 1988, to September 30, 2002) in the Atascadero Creek watershed located along the southern coast of California for 1929, 1998, and 2050 (estimated) land use conditions (8, 38 and 52 percent urban, respectively). The 14‐year period experienced a range of climatic conditions caused mainly by El Nino‐Southern Oscillation variations. A geographic information system is used to delineate the watershed and parameterize the model, which is calibrated using data from two stream flow and eight rainfall gauges. Urbanization is shown to increase peak discharges and runoff volume while decreasing stream flow variability. In all cases, the annual and 14‐year distributions of stream flow are shown to be highly skewed, with the annual maximum 24 hours of discharge accounting for 22 to 52 percent of the annual runoff and the maximum ten days of discharge from an average El Nino year producing 10 to 15 percent of the total 14‐year discharge. For the entire period of urbanization (1929 to 2050), the average increase in annual maximum discharges and runoff was 45 m³/s (300 percent) and 15 cm (350 percent), respectively. Additionally, the projected increase in urbanization from 1998 to 2050 is half the increase from 1929 to 1998; however, increases in runoff (22 m³/s and 7 cm) are similar for both scenarios because of the region's spatial development pattern.     

The Potential Use of Synthetic Aperture Radar for Estimating Methane Ebullition From Arctic Lakes1

Abstract: Arctic lakes are significant emitters of methane (CH₄), a potent greenhouse gas, to the atmosphere; yet no rigorous quantification of the magnitude and variability of pan‐Arctic lake emissions exists. In this study, we demonstrate the potential for a new method using synthetic aperture radar (SAR) imagery to detect methane bubbles in lake ice to scale up whole‐lake measurements of CH₄ ebullition (bubbling) to regional scales. We estimated ebullition from lakes, which is often the dominant mode of lake emissions, by mapping the distribution of bubble clusters frozen in early winter ice across surfaces of seven tundra lakes and one boreal forest lake in Alaska. Applying previously measured ebullition rates associated with four distinct classes of bubble clusters found in lake ice, we estimated whole‐lake emissions from individual lakes. The percent surface area of lake ice covered with bubbles (R² = 0.68) and CH₄ ebullition rates from lakes (R² = 0.59) and were correlated with radar return values from RADARSAT‐1 Standard Beam mode 3 for the tundra lakes, suggesting that with appropriate scaling and consideration for variability in lake‐ice conditions, this technique has the potential to be used for estimating broader‐scale regional and pan‐Arctic lake methane emissions.     

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.     

Empirical Estimation of Stream Discharge Using Channel Geometry in Low‐Gradient,Sand‐Bed Streams of the Southeastern Plains

Manning's equation is used widely to predict stream discharge (Q) from hydraulic variables when logistics constrain empirical measurements of in‐bank flow events. Uncertainty in Manning's roughness (n_M) is the major source of error in natural channels, and sand‐bed streams pose difficulties because flow resistance is affected by flow‐dependent bed configuration. Our study was designed to develop and validate models for estimating Q from channel geometry easily derived from cross‐sectional surveys and available GIS data. A database was compiled consisting of 484 Q measurements from 75 sand‐bed streams in Alabama, Georgia, South Carolina, North Carolina (Southeastern Plains), and Florida (Southern Coastal Plain), with six New Zealand streams included to develop statistical models to predict Q from hydraulic variables. Model error characteristics were estimated with leave‐one‐site‐out jackknifing. Independent data of 317 Q measurements from 55 Southeastern Plains streams indicated the model (Q = A_cR_H^0.6906S^0.1216; where A_c is the channel area, R_H is the hydraulic radius, and S is the bed slope) best predicted Q, based on Akaike's information criterion and root mean square error. Models also were developed from smaller Q range subsets to explore if subsets increased predictive ability, but error fit statistics suggested that these were not reasonable alternatives to the above equation. Thus, we recommend the above equation for predicting in‐bank Q of unbraided, sandy streams of the Southeastern Plains.     

Effect of Flow Depth and Velocity on Nitrate Loss Rates in Natural Channels1

Carleton, James N. and Yusuf M. Mohamoud, 2012. Effect of Flow Depth and Velocity on Nitrate Loss Rates in Natural Channels. Journal of the American Water Resources Association (JAWRA) 1‐12. DOI: 10.1111/jawr.12007 Abstract: Loss rates of nitrate from streams and rivers are governed by movement of the ion from water column to anoxic bed sediments. Quantitative representations of nitrate in streams and rivers have often treated such losses as governed by first‐order mechanisms that are invariant with respect to potential modulating factors other than temperature. Results of studies in recent years, however, suggest that rates of water column‐sediment mass transfer are influenced by stream geometry and associated hydraulics. We develop expressions for the instream nitrate loss rate coefficient, k, as a function of water velocity and depth, using hydraulic geometry to empirically relate velocity to depth for two cases: (1) variability in mean conditions among reaches; and (2) temporal variability in conditions at a single reach, under changing flow. The result is expressions for k as functions of water column depth. Measured stream k values reported in the literature are shown to be well represented by expressions developed for the first case, and the potential for application to probabilistic analysis is briefly examined. We explore the latter case using the Hydrologic Simulation Program – FORTRAN (HSPF) model, modified to incorporate the dependence of k on instantaneous stream depth. In example simulations of two nitrate‐exporting watersheds, the incorporation of depth‐dependence of k produces improvement in the model’s ability to match observed stream nitrate concentrations.     

Decoupling Streamflow Responses to Climate Variability and Land Use/Cover Changes in a Watershed in Northern China

Restored annual streamflow (Q_r) and measured daily streamflow of the Chaohe watershed located in northern China and associated long‐term climate and land use/cover data were used to explore the effects of land use/cover change and climate variability on the streamflow during 1961‐2009. There were no significant changes in annual precipitation (P) and potential evapotranspiration, whereas Q_r decreased significantly by 0.81 mm/yr (p < 0.001) over the study period with a change point in 1999. We used 1961‐1998 as the baseline period (BP) and 1999‐2009 the change period (CP). The mean Q_r during the CP decreased by 39.4 mm compared with that in the BP. From 1979 to 2009, the grassland area declined by 69.6%, and the forest and shrublands increased by 105.4 and 73.1%, respectively. The land use/cover change and climate variability contributed for 58.4 and 41.6% reduction in mean annual Q_r, respectively. Compared with the BP, median and high flows in the CP decreased by 38.8 and up to 75.5%, respectively. The study concludes that large‐scale ecological restoration and watershed management in northern China has greatly decreased water yield and reduced high flows due to the improved land cover by afforestation leading to higher water loss through evapotranspiration. At a large watershed scale, land use/cover change could play as much of an important role as climate variability on water resources.