Published online 27 October 2006
Published in J Environ Qual 35:2202-2211 (2006)
DOI: 10.2134/jeq2006.0042
© 2006 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
TECHNICAL REPORTS
Surface Water Quality
Nitrate and Sediment Fluxes from a California Rangeland Watershed
David J. Lewisa,
Michael J. Singerb,*,
Randy A. Dahlgrenb and
Kenneth W. Tatec
a University of California Cooperative Extension, Santa Rosa, CA 95403
b Department of Land Air and Water Resources, University of California, Davis, CA 95616
c Department of Plant Sciences, University of California, Davis, CA 95616
* Corresponding author (mjsinger{at}ucdavis.edu)
Received for publication January 27, 2006.
 |
ABSTRACT
|
|---|
Long-term water quality records for assessing natural variability, impact of management, and that guide regulatory processes to safeguard water resources are rare for California oak woodland rangelands. This study presents a 20-yr record (19812000) of nitrate-nitrogen (NO3N) and suspended sediment export from a typical, grazed oak woodland watershed (103 ha) in the northern Sierra Nevada foothills of California. Mean annual precipitation over the 20-yr period was 734 mm yr1 (range 3661205 mm yr1). Mean annual stream flow was 353 mm y1 (range 87848 mm yr1). Average annual stream flow was 48.1 ± 16% of precipitation. Mean annual NO3N export was 1.6 kg ha1 yr1 (range 0.183.6 kg ha1 yr1). Annual NO3N export significantly (P < 0.05) increased with increasing annual stream flow and precipitation. Mean daily NO3N export was 0.004 kg ha1 d1 (range 105 to 0.55 kg ha1 d1). Mean annual suspended sediment export was 198 kg ha1 yr1 (range 23479 kg ha1 yr1). There was a positive relationship (P < 0.05) between annual suspended sediment export, annual stream flow and precipitation. Mean daily suspended sediment export was 0.54 kg ha1 d1 (range 104 to 155 kg ha1 d1). Virtually no sediment was exported during the dry season. The large variation in daily and annual fluxes highlights the necessity of using long-term records to establish quantitative water quality targets for rangelands and demonstrates the difficulty of designing a water quality monitoring program for these ecosystems.
Abbreviations: LOD, limit of detection TMDL, total maximum daily load
 |
INTRODUCTION
|
|---|
APPROXIMATELY 3 million hectares of oak woodlands and savanna in the interior valleys and foothills of California are used extensively for cattle (Bos taurus) grazing, firewood production, wildlife habitat, and other land uses (Bolsinger, 1987; Griffin, 1988; Standiford and Howitt, 1993). Urbanization is rapidly occurring across parts of the region with potentially negative consequences for wildlife habitat, water quality, and other ecosystem services (Wacker and Kelly, 2004). These rangelands contribute significant surface drinking water supplies, with two-thirds of all drinking water reservoirs in the State located within grazed oak woodland/annual grassland ecosystems (FRRAP, 1988). Concerns about drinking water quality degradation and declines in fisheries, such as coho salmon (Oncorhynchus kisutch) and steelhead trout (Oncorhynchus mykiss), have lead to the listing of associated water bodies (e.g., rivers, reservoirs) as impaired following Section 303d of the Federal Clean Water Act (Lewis et al., 2002). Impairment indicates that the water body is not achieving its designated beneficial use(s) (e.g., municipal drinking water, cold water fisheries habitat). During the listing process, water resources agency staff must identify the cause of the impairment (e.g., excessive N and sediment) and land use(s) generating the impairment. Once a water body is listed as impaired, water resources protection agencies must develop and implement watershedwater body specific maximum allowable loading targets known as total maximum daily loads (TMDLs) for the specific nonpoint source pollutant of concern.
Development of the allowable loading limit and a watershed-scale strategy to achieve it depends on accurate estimation, among other things, of: (i) loading rates from specific land uses as well as natural background levels within the watershed; (ii) load reductions required from each land use to achieve the beneficial use(s); and (iii) watershed management scenarios that achieve the required load reductions (USEPA, 1991, 1999a, 1999b). Load estimations can and have been made from watershed specific water quality data, from watershedwater quality simulation models developed and calibrated with data representative of local watershed conditions, and from remote techniques such as aerial photographic review of historic and current erosion features (CRWQCB, 2002). By definition, a TMDL is the total maximum daily load (kg d1) of specified constituents that a water body can receive and still meet its designated beneficial use. However, in application a TMDL can take the shape of an annual load estimation for sediment (CRWQCB, 2002), or a seasonal calculation for N and P derived from a theoretical or empirical watershed-scale water quality model (ADEQ, 2000).
Long-term, watershed-scale water quality data are invaluable in the TMDL development process, and in evaluation of progress toward achieving loading reductions and beneficial uses of listed water bodies. In general, there are insufficient long-term data available to confidently estimate loading rates for specific land-use types, or to develop and calibrate watershed-scale simulation models to estimate water body response to various load reduction scenarios. This is certainly the case for California's oak woodlands and annual grasslands, where significant annual, seasonal, and storm-event variability in watershed-scale flow and water quality have been demonstrated at multiple spatial scales (Tate et al., 1999; Lewis et al., 2000, 2002; Holloway and Dahlgren, 2001; Ahearn et al., 2004; Dahlgren et al., 2005). The TMDL development guidelines acknowledge the need to account for inherent variability (e.g., annual weather fluctuation) in the establishment of maximum allowable loading targets via the addition of a margin of safety to the maximum allowable loading estimate (USEPA, 1991). However, without long-term records of weather, streamflow, and water quality, it is extremely difficult to establish credible margins of safety that incorporate natural variability into estimates of daily, seasonal, or annual loading estimates, without creating undue economic hardship and managerial constraints for land managers.
The purpose of this paper is to report a unique, continuous 20-yr record (19812000) of NO3N and suspended sediment export from a typical, grazed annual rangeland watershed (103 ha) in the northern Sierra Nevada foothills of California. Impairment by nutrients and sediments, along with microbial pathogens, are the primary water quality concerns on California rangelands. We report water quality data as daily and annual loading rates with the goal that they have utility for: (i) development of maximum allowable loading rates from rangeland watersheds; (ii) incorporation of inherent temporal variability in development of margins of safety for maximum allowable loading rate estimations; and (iii) development of loading rate reduction targets for annual rangeland watersheds. In addition, we report statistical relationships identified between precipitation and stream flow, and annual suspended sediment and NO3N fluxes.
 |
METHODS
|
|---|
Study Area
The study site is located on the University of California's Sierra Foothill Research and Extension Center, approximately 96 km northeast of Sacramento (Fig. 1) (39°15' N, 121°17' W). Elevations in the 103 ha Schubert watershed range from 152 to 427 m. The average slope is 18% with a range from 2 to 50%. Climate is Mediterranean with more than 90% of annual precipitation falling from November to March as rainfall. Over the 20-yr monitoring period, mean annual rainfall (±1 SD) was 734 ± 253 mm yr1 (range 3661205 mm yr1) and median annual rainfall was 684 mm yr1. Mean annual stream flow (±1 SD) was 353 ± 225 mm yr1 (range 87848 mm yr1) and median annual stream flow was 330 mm yr1. The mean annual air temperature has been 15°C.
Vegetation at the site is typical of Sierran foothill oak woodlands that are dominated by blue oak (Q. douglasii Hook and Arn.) with minor interior live oak (Q. wislizenii A.DC.) and foothill pine (Pinus sabiniana Douglas) (Griffin, 1988). Vegetation patterns across the watershed are a mosaic of open grassland and woodland (Epifanio, 1989; Jansen, 1987). Blue oak density ranges from 90 to 200 trees per hectare with approximately 70% average canopy cover for the watershed. Annual grasses and legumes dominate the ground cover, with differing species diversity under the oak canopy compared to open grasslands (Jackson et al., 1990; Jansen, 1987).
Soils in the watershed are Ruptic-Lithic Xerochrepts on the steep side slopes and Mollic Haploxeralfs on the more level areas (Lytle, 1998). Soils formed in residuum from basic metavolcanic (greenstone) bedrock (Beiersdorfer, 1979) are shallow to moderately deep, medium textured, gravelly, and rocky. Surface horizons have a high organic C content (63 g kg1) and a relatively high infiltration rate (1.8 cm h1) (Huang, 1997). Soils are rich in Fe oxides, with a mean dithionite-extractable Fe content of 31 and 35 g kg1 in the A/AB and Bt horizons, respectively. They have very stony clay subsoils with an average clay content of 400 g kg1 and a saturated hydraulic conductivity of about 0.68 cm h1 (Huang, 1997). The horizon in contact with the weathered bedrock had the lowest saturated hydraulic conductivity (0.012 cm h1), which in conjunction with the shallow bedrock (11.5 m depth) limits deep percolation of water (Dahlgren and Singer, 1994).
The watershed is used for light to moderate seasonal beef cattle grazing primarily from January to March and August to October. Average stocking rate over the 20-yr period was 0.22 animal unit months ha1 (Lewis, 1998). Grazing levels on the watershed for the 20-yr record fell within recommended levels for moderate grazing, as defined by October residual dry matter levels (6001000 kg ha1), on annual rangelands given site specific land slope and annual rainfall levels (Clawson et al., 1982; Bartolome et al., 2002). In addition to normal grazing management, 1352 oaks were selectively harvested from 15 ha in the northern upland portion of the watershed between July 1984 and May 1986 (Epifanio, 1989). Trees were removed during the spring and summer months when soils were dry. Trees were not removed from drainage-ways or from slopes steeper than 30%.
Streamflow Measurement and Water Quality Sampling
Hydrology and water quality have been measured in the watershed since 1978 with dependable records starting in 1981. The watershed was instrumented with a 0.912-m (3-foot) Parshall flume for storm flow measurements and a tandem 0.304-m (1-foot) 90° V-notch weir for base flow measurements (Lewis et al., 2000). Stage height was continuously recorded with stage height recorders located within stilling wells connected to the flume and weir. Stage height was converted to flow using standard methods (Linsley et al., 1982). Rainfall amount and intensity were recorded with tipping bucket rain gauges and recorders at 198, 290, and 328 m elevations within the watershed. Using this continuous hydrologic record, annual rainfall and runoff were calculated as the sum of all rainfall events or discharge measured during each water year (1 October30 September). No storm hydrographs were missed during the 20-yr period, but water flow during base-flow periods was not measured at times due to equipment failure. Flows for these times were estimated as linear changes between measured flows.
Stream-water samples were collected during rain-free periods, as well as on a storm-event basis, from the center of the stream 10 m upstream from the flume and weir using an ISCO (ISCO, Lincoln, NE) automatic pump sampler. Sampling frequency over the 20-yr record was more intense during storms (about one sample every 14 h) than during rain-free periods (about onetwo samples per 7 d). This sampling protocol is consistent with the concept of flow proportional sampling. Samples were refrigerated at 3°C through completion of all analyses. Samples were filtered through a 0.45-µm membrane filter and analyzed for pH, electrical conductivity, total alkalinity, suspended solids, and Ca2+, K+, Na+, Mg2+, Cl, NH4+, NO3, ortho-PO43, and SO42 concentrations. Only suspended sediment and NO3N data are presented here. The weight of suspended sediments was determined by the difference in weight of the 0.45-µm filter before and after filtering a known volume of water (LOD = 0.1 mg L1). Nitrate concentrations were determined using a Dionex ion chromatograph (LOD = 0.010.04 mg N L1).
Constituent fluxes or loads were calculated using a period-weighting method (Moldan and Cerny, 1994). This method multiplied the average of the concentration at the beginning and end of a time interval by the discharge volume during that time interval. Any missing baseflow sample concentrations were estimated from the 20-yr average base flow concentration for each constituent. Daily loads were then summed to generate annual loads. This method was chosen based on recommendations from studies comparing various flux calculation methods (Semkin et al., 1994; Dann et al., 1986). The period-weighting method is commonly employed in watershed studies (e.g., Likens et al., 1995) and provides daily and annual load estimations.
Statistical Analysis
Standard probability plots were created using log transformed daily NO3N and sediment masses. Daily values were estimated using the period-weighting method described previously. Simple linear regression analysis was used to determine relationships among annual precipitation, stream flow, and suspended sediment and NO3N fluxes from the watershed. Annual suspended sediment and NO3N loads were log10 transformed to provide a normal distribution of the data.
 |
RESULTS AND DISCUSSION
|
|---|
Precipitation and Stream Flow
Precipitation and stream flow varied significantly over the 20-yr record (Fig. 2). Mean annual rainfall (±1 SD) was 734 ± 253 mm yr1 (range 3661205 mm yr1) and mean annual stream flow (±1 SD) was 353 ± 225 mm yr1 (range 87848 mm yr1). Average annual stream flow, as a percentage of precipitation was 48.1 ± 16% and ranged from 19 to 76%. Annual stream flow increased significantly with increased annual precipitation (Table 1). We found no measurable change in stream flow that could be attributed to selective removal of oak trees early in the 20-yr study period (data not shown).
View this table:
[in this window]
[in a new window]
|
Table 1. Linear regression equations depicting relationships between annual rainfall (P), stream flow (Q), total suspended solid load (TSS), and nitrate-nitrogen load (NO3N) discharged from the Schubert Watershed over a 20-yr period (19812000).
|
|
Precipitation amount and timing during the rainy season were important factors regulating the stream flow/rainfall ratio. Over the course of the rainy season, the first storms produced no appreciable change in stream flow (Lewis et al., 2000). A mean precipitation depth of 47 mm (range 2669 mm) was required to "prime" the watershed and generate a measurable stream response to precipitation after the dry summer. Once the watershed soils reached field capacity, each storm event produced a rapid response in the stream hydrograph. This is illustrated by the 1984 water year (October 1983September 1984), which had near average precipitation for the 20-yr record. Four early season precipitation events totaling 37.8 mm produced only a small increase in the stream flow. The next 8 mm produced a notable rise in the stream hydrograph and later precipitation events produced large and distinct increases in the hydrograph (Fig. 3). In the spring, the hydrologic response of the stream declined as storm size and intensity declined along with an increase in evapotranspiration due to increasing temperature and plant growth (data not shown). Thus, water yields are greatest for winter season precipitation events following soil water replenishment in the fall and before elevated evapotranspiration rates in the spring.

View larger version (15K):
[in this window]
[in a new window]
|
Fig. 3. Early season stream flow and precipitation for 1984 illustrating early season "priming" of the watershed soils.
|
|
Nitrate-Nitrogen
We found considerable variability in the year-to-year NO3N export from the watershed (Fig. 4). Mean annual NO3N export from the watershed was 1.6 kg ha1 yr1 (range 0.183.6 kg ha1 yr1) (Fig. 4). Median NO3N export was 1.2 kg ha1 yr1. Daily NO3N export was also highly variable, with the greatest loads occurring between November and March (Fig. 5b). Mean daily NO3N export was 0.004 kg ha1 d1 (range 105 to 0.55 kg ha1 d1) (Fig. 5b). Nitrate was the only significant form of mineral N exported by the stream from the watershed; average NH4N export was <0.1 kg ha1 yr1. While annual NO3N fluxes reached values as high as 3.6 kg N ha1 yr1, loss of NO3 from the watershed was substantially less than the combined inputs of mineral N (NH4+ + NO3) in the bulk precipitation (mean 9.4 and maximum 17.8 kg N ha1 yr1). As a result, there is an apparent net retention of mineral N within the watershed. We did not measure organic forms of N or gaseous losses, and have not accounted for N exported by animal consumption of the grass, so we were not able to calculate a complete mass balance for N in the watershed. We expect volatilization losses to be small. Herman et al. (2003) measured gas emissions of <1 ng N cm2 h1 and concluded that N volatilization losses from this watershed were not an important pathway for N loss.

View larger version (29K):
[in this window]
[in a new window]
|
Fig. 5. (A) Mean daily stream flow for the Schubert watershed 1981 to 2000. (B) Mean daily nitrate-nitrogen (NO3N) load for the Schubert watershed 1981 to 2000.
|
|
Annual NO3N load was positively and significantly (P < 0.05) related to both annual precipitation and stream flow (Table 1). The years with the highest export were a function of a combined very high flow volume and high NO3N concentrations during some storm events. Daily export was often high during the early portion of the rainy season when NO3N concentrations are generally highest in stream water (Fig. 5a, 5b, and Fig. 6). In California oak woodlands, stream water NO3 concentrations are highest during the early rainfall season due to a temporal asynchrony in the N cycle (Tate et al., 1999; Holloway and Dahlgren, 2001; Jackson et al., 2006). Once early season rainfall begins to produce a significant change in the stream hydrograph, it also flushes N from the system (Fig. 6). Subsequent storms have lower N concentrations as the season progresses (Fig. 6). The reason for this behavior is that NO3 accumulates in the soil during the summer and fall periods, when plant uptake and leaching processes are very low. Nutrient uptake by oak trees is limited during the summer due to the soil water deficit and they become dormant in the fall and remain inactive until late March. Similarly, annual grasses die in June, and seeds germinate with the first fall precipitation. However, their growth is limited by cold temperatures until about the beginning of March. Therefore, NO3 that accumulates in the soil over the dry period is highly susceptible to leaching. The accumulated NO3 is progressively flushed out of the soil, resulting in lower NO3N concentrations in stream water with each subsequent storm. The export of NO3N decreases throughout the spring as plant uptake removes NO3 from the soil zone. The decreased NO3N fluxes also result from a progressive decrease in baseflow from spring to summer (Fig. 5a). Thus, early to midwinter precipitation events result in the highest yields of NO3N from the watershed.

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 6. Early season stream flow and nitrate-nitrogen concentration for the 1984 water year, and early season precipitation events for the 1984 water year.
|
|
Suspended Sediment
Mean annual suspended sediment export from the watershed was 198 kg ha1 yr1 (range 23479 kg ha1 yr1) (Fig. 7). Annual suspended sediment export was significantly (P < 0.05) related to annual stream flow and precipitation, with increasing export as annual stream flow and precipitation increased (Table 1). Mean daily suspended sediment export was 0.54 kg ha1 d1 (range 104 to 155 kg ha1 d1) (Fig. 8b). Virtually no sediment was exported from the watershed during the dry season (Fig. 8b). Maximum daily suspended sediment export was linked to high sediment concentrations and high stream flow occurring during large storm events (Fig. 8a).

View larger version (28K):
[in this window]
[in a new window]
|
Fig. 8. (A) Mean daily stream flow for the Schubert watershed 1981 to 2000. (B) Mean daily suspended sediment load for the Schubert watershed 1981 to 2000.
|
|
Suspended sediment exports are small relative to several other annual rangeland watersheds with similar hydrologic characteristics, but differing soil types (Lewis et al., 2002). The high concentrations of secondary Fe oxides and organic matter promote high infiltration rates and strong aggregate stability, leading to low erosivity. The high canopy coverage (70%) and adequate herbaceous cover (6001000 kg ha1 October residual dry matter) further attenuate surface runoff and erosion by limiting raindrop impact on the soil surface. Ground cover by herbaceous plants exceeds 90% in the watershed during the rainy season, helping to ensure a low erosion rate.
Implications for the Total Maximum Daily Load Development Process
These data have several important implications for the TMDL process when applied to annual rangeland watersheds, and illustrate the overall value of investment in long-term, strategic monitoring by agencies involved in water quality management and protection. An important use of long-term hydrology and water quality data is to calibrate and validate watershed-scale water quality models to conditions typical of annual rangeland watersheds. The use of such models will be required to address issues of spatial scale and variability in source assessments, to bridge data gaps, and to evaluate the potential effectiveness of alternative watershed management scenarios to reduce pollutant loadings. To have confidence in loading estimates and watershed management plans developed during the TMDL process, the models employed to develop these estimates and plans must be evaluated against real watershed-scale data. For example, Lewis et al. (2000) were able to use a portion of this dataset to evaluate methods of applying the Curve Number rainfall-runoff model for annual rangeland watersheds providing credibility to an important hydrologic modeling component common among watershed-scale water quality models.
To the best of our knowledge, these data represent the most comprehensive, watershed scale assessment of long-term NO3N and suspended sediment loading rates for moderately grazed annual rangeland watersheds. Under the grazing levels and climatic conditions realized during the 20 yr of data collection the watershed actually served as a sink for mineral N deposited as dry and wet atmospheric deposition. While there was certainly export of NO3N from the watershed on a daily and annual basis, these data raise the possibility that annual rangeland watersheds sequester more N than they generate.
Approximately 90% of annual rangelands in the region are privately owned and 80% are grazed by domestic livestock (Wacker and Kelly, 2004), so an important question to address is how grazing management influences the capacity of these watersheds to retain N and suspended sediments. It is likely that many watersheds throughout the region may be grazed more heavily than the studied watershed, and certainly some watershed management proposals call for elimination of grazing by domestic livestock. These data demonstrate that moderate intensity grazing, such as that practiced within the experimental watershed, does not contribute to water quality degradation by N or sediment compared to many agricultural and range lands.
These data also demonstrate the temporal variability in loading rate at both the daily and annual scale. Fundamentally, these data illustrate that the establishment of a set daily or annual loading rate (either as a regulatory target, or as a management expectation) ignores basic hydrologic and biogeochemical processes that dictate background N and sediment flux from these watersheds. For instance, it is clear from these data that a maximum allowable loading rate and expected loading reductions should be scaled according to season. At a minimum, one could justify different loading rates for dry versus wet season (Fig. 5 and 8). These data and the existing literature (Tate et al., 1999; Holloway and Dahlgren, 2001; Jackson et al., 2006) justify application of a different loading rate expectation for early vs. late wet season NO3N. There is an inherent flushing of NO3N from annual rangelands due to asynchrony between N availability and plant demand, and hydrologic transport events early in the wet season. Even in the absence of grazing, California oak woodland/annual grasslands are susceptible to NO3N leaching due to this asynchrony in the N cycle. This highlights the importance of considering the natural background contributions of pollutants, such as nutrients and suspended sediments, before setting regulatory criteria.
The temporal variability in NO3 and suspended sediment concentrations documented by this study provides a valuable record to assist in designing a cost-effective and accurate monitoring program. For several constituents, the majority of the annual constituent flux is transported during a few large storm events each year. Because constituent concentrations may change appreciably during storm events, intensive sampling of the storm hydrograph is required to adequately characterize fluxes. In contrast, constituent concentrations during baseflow conditions typically display small changes and therefore do not need to be sampled at a high frequency.
The probability plots reported in Fig. 9 (n = 7300) represent a useful, risk-based approach to assigning a margin of safety to maximum allowable loadings, as well as selecting reasonable allowable loading estimates. For example, given the probability information, a watershed manager or regulator could identify a maximum allowable loading for NO3N based on probability of occurrence rather than based on a mean or median value. Thus, a defendable risk of exceedance-based margin of safety could be selected to accommodate watershed specific risk tolerance both from a water quality and stake holder's risk tolerance perspective.

View larger version (15K):
[in this window]
[in a new window]
|
Fig. 9. Probability plots of daily nitrate-nitrogen (NO3N) and suspended sediment loads in stream water for the Schubert watershed.
|
|
The 20-yr record of NO3 and sediment flux reported here illustrates the importance of long-term monitoring and data collection to facilitate management and regulatory decision-making. These data indicate that NO3 and suspended sediment fluxes in stream water from annual rangelands are highly variable at the daily and annual time scales. Year-to-year NO3N and suspended sediment fluxes may vary by over an order of magnitude, making long-term data records essential for evaluating watershed input-output budgets and establishment of allowable loading rates for nonpoint source constituents from rangelands.
 |
ACKNOWLEDGMENTS
|
|---|
Partial support for this work came from the Integrated Hardwood Range Management Program, International Agricultural Development Graduate Group and the Jastro-Shields fund. Thanks to Mike Conner, Dave Labadie, Don Springsteen, Tim Clark, and the Sierra Foothills Research and Extension Center staff for continued support in sampling and instrumentation. Their work combined with that of Peter Janitzky, Charlette Epifanio, and Xiaohong Huang, as researchers and custodians of the Schubert Watershed record, have generated the long-term data used in this study.
 |
REFERENCES
|
|---|
- Ahearn, D.S., R.W. Sheibley, R.A. Dahlgren, and K.E. Keller. 2004. Temporal dynamics of stream water chemistry in the last free-flowing river draining the western Sierra Nevada, California. J. Hydrol. 295:4763.[CrossRef]
- Arizona Department of Environmental Quality. 2000. Luna Lake TMDL. Prepared by Tetra Tech for the ADEQ, Phoenix, AZ.
- Bartolome, J.W., W.E. Frost, N.K. McDougald, and M. Connor. 2002. California guidelines for residual dry matter (RDM) management on coastal and foothill annual rangelands. Publ. 8092. Univ. of California, Div. Agric. and Nat. Resources, Oakland.
- Beiersdorfer, R.E. 1979. Metamorphic petrology of the Smartville Complex, Northern Sierra Nevada Foothills. M.S. thesis. Univ. of California, Davis.
- Bolsinger, C.L. 1987. Major findings of a statewide resource assessment in California. p. 291297. In T.R. Plumb and N.H. Pillsbury (tech. coord.) Proc. Symp. Multiple-use management of California's hardwood resources, San Luis Obispo, CA. 1214 Nov. 1986. Gen Tech. Rep. PSW-100. Pacific Southwest Forest Range Exp. Stn., USDA Forest Serv., Berkeley, CA.
- Clawson, W.J., N.K. McDougald, and D.A. Duncan. 1982. Guidelines for residue management on annual range. Leaflet 21327. Univ. of California, Div. of Agric. and Nat. Resources, Oakland.
- Coast Region Water Quality Control Board. 2002. Action plan for the Garcia River Watershed Sediment TMDL. North Coast Regional Water Quality Control Board, Santa Rosa, CA.
- Dahlgren, R.A., and M.J. Singer. 1994. Nutrient cycling in managed and non-managed oak woodlandgrass ecosystem. Final Report to Integrated Hardwood Range Management Program. Land, Air and Water Resources Pap. 100028. Univ. of California, Davis.
- Dahlgren, R.A., K.W. Tate, and D.S. Ahearn. 2005. Watershed scale, water quality monitoring: Water sample collection. p. 547564. In R.D. Down and J.H. Lehr (ed.) Handbook of environmental instrumentation. John Wiley & Sons, New York.
- Dann, M.S., J.A. Lynch, and E.S. Corbett. 1986. Comparison of methods for estimating sulfate export from a forested watershed. J. Environ. Qual. 15:140145.[Abstract/Free Full Text]
- Epifanio, C.R. 1989. Hydrologic impacts of blue oak harvesting and evaluation of the modified USLE in the northern Sierra Nevada. M.S. thesis. Univ. of California, Davis.
- Forest and Rangeland Resources Assessment Program. 1988. California's forests and rangelands: Growing conflict over changing uses. Forest and Rangeland Resources Assessment Program. California Dep. of Forestry and Fire Protection.
- Griffin, J. 1988. Oak woodland. p. 383415. In M. Barbour and J. Major (ed.) Terrestrial vegetation of California. John Wiley & Sons, New York.
- Herman, D.J., L.J. Halverson, and M.K. Firestone. 2003. Nitrogen dynamics in an annual grassland: Oak canopy, climate, and microbial population effects. Ecol. Appl. 13:593604.
- Holloway, J.M., and R.A. Dahlgren. 2001. Seasonal and event-scale variations in solute chemistry for four Sierra Nevada catchments. J. Hydrol. 250:106121.[CrossRef]
- Huang, X. 1997. Watershed hydrology, soil, and biogeochemistry in an oak woodland annual grassland ecosystem in the Sierra Foothills, California. Ph.D. diss. Univ. of California, Davis.
- Jackson, L.E., R.B. Strauss, M.K. Firestone, and J.W. Bartolome. 1990. Influence of tree canopies on grassland productivity and nitrogen dynamics in deciduous oak savanna. Agric. Ecosyst. Environ. 32:89105.[CrossRef]
- Jackson, R.D., B. Allen-Diaz, L.G. Oates, and K.W. Tate. 2006. Spring-water nitrate increased by grazing removal in a Californian oak savanna. Ecosystems 9:115.[CrossRef]
- Jansen, H. 1987. The effect of blue oak removal on herbaceous production on a foothill site in the northern Sierra Nevada. p. 343350. In T.R. Plumb and N.H. Pillsbury (tech. coord.) Proc. Symp. Multiple-Use Management of California's Hardwood Resources, San Luis Obispo, CA. 1214 Nov. 1986. Gen. Tech. Rep. PSW-100. Pacific Southwest Forest Range Exp. Stn., USDA Forest Serv., Berkeley, CA.
- Lewis, D.J. 1998. Water quality and yield in a California oak woodland watershed: A 17-year study. M.S. thesis. Univ. of California, Davis.
- Lewis, D.J., M.J. Singer, R.A. Dahlgren, and K.W. Tate. 2000. Hydrology in a California oak woodland watershed: A 17 year study. J. Hydrol. 240:106117.[CrossRef]
- Lewis, D.J., K.W. Tate, R.A. Dahlgren, and J. Newell. 2002. Turbidity and total suspended solid concentration dynamics in streamflow from California oak woodland watersheds. p. 107118. In R.B. Standiford et al. (tech. coord.) Proc. Symp. on Oak Woodlands in California's Changing Landscape, San Diego, CA. 2225 Oct. 2001. Gen. Tech. Rep. PSW-184. Pacific Southwest Forest and Range Exp. Stn., Berkeley, CA.
- Likens, G.E., F.H. Bormann, R.S. Pierce, J.S. Eaton, and N.M. Johnson. 1995. Biogeochemistry of a forested ecosystem. Springer-Verlag, New York.
- Linsley, R.K., Jr., M.A. Kohler, and J.J.H. Paulhus. 1982. Hydrology for engineers. 3rd ed. McGraw-Hill, New York.
- Lytle, D.J. 1998. Soil survey of Yuba County, California. USDA Nat. Res. Conserv. Serv., Washington, DC.
- Moldan, B., and J. Cerny. 1994. Small watershed research. p. 189206. In B. Moldan and J. Cerny (ed.) Biogeochemistry of small watersheds: A tool for environmental research. Scope 51. John Wiley & Sons, New York.
- Semkin, R.G., D.S. Jeffries, and T.A. Clair. 1994. Hydrochemical methods and relationships for study of stream output from small watersheds. p. 163188. In B. Moldan and J. Cerny (ed.) Biogeochemistry of small watersheds: A tool for environmental research. Scope 51. Wiley & Sons, New York.
- Standiford, R.B., and R.E. Howitt. 1993. Multiple use of California's hardwood rangelands. J. Range Manage. 46:176182.
- Tate, K.W., R.A. Dahlgren, M.J. Singer, B. Allen-Diaz, and E.R. Atwill. 1999. Timing frequency of sampling affect accuracy of water-quality monitoring. Calif. Agric. 53:4448.
- USEPA. 1991. Guidance for water quality-based decisions: The TMDL process. USEPA Rep. 440/4-91-001. USEPA, Washington, DC.
- USEPA. 1999a. Protocol for developing nutrient TMDLs. USEPA Rep. 841/B-99-007. USEPA, Washington, DC.
- USEPA. 1999b. Protocol for developing sediment TMDLs. USEPA Rep. 841/B-99-004. USEPA, Washington, DC.
- Wacker, M.J., and N.M. Kelly. 2004. Ranchers vs. ranchettes in California's oak rangelands. Rangelands 26:1722.