Research connects soil hydrology and stream water chemistry in California oak woodlands
The UC Sierra Foothill Research and Extension Center (SFREC) is located in the heart of typical California blue oak and live oak woodlands within metavolcanic terrain of the Sierra Nevada foothills. These types of woodlands often exist at the interface between urban, wild and agricultural lands and are used extensively for livestock grazing, wildlife habitat and surface water supply. Soil surveys for this region and within SFREC depict relatively few soil types compared to areas that support more-intensive agricultural land uses. Despite this inferred homogeneity, our study showed that the biogeochemical and physical properties of soils vary sharply over short distances of less than 10 feet and also experience changes by season and as a result of storm events. An understanding of soil variability in this setting is important to assess rangeland productivity, perennial grass and oak restoration potential, carbon sequestration, stream flow generation and stream water chemistry.
UC Davis undergraduate Tony Orozco installs a sensor in the Lewis-1 watershed at the UC Sierra Foothill Research and Extension Center.
The bracket shows an abrupt clay increase associated with the upper boundary of the claypan in a typical soil profile at SFREC.
Oak woodlands are a patchy array of open grassland and oak canopy. The presence or absence of trees is one factor that imposes spatial variability in soil characteristics. Oak trees change the properties of soils beneath their canopies through a variety of nutrient-cycling processes, such as annual inputs of litter fall (Dahlgren et al. 2003). As leaves drop and decompose at the soil surface, organic matter and nutrients are resup-plied to the upper soil horizons, creating islands of fertility and enhanced soil quality. Physical (bulk density, water infiltration rates), chemical (nutrient enrichment, pH, cation exchange capacity) and biological (microbial biomass, soil respiration) properties are enhanced in soils under oak canopies (Dahlgren et al. 1997). In contrast, soils forming under annual grasses in the absence of oak canopy are less fertile, have higher bulk density and are more susceptible to surface-water runoff.
Soil processes exhibit unique patterns in Mediterranean ecosystems such as California, due to the timing of precipitation. Summers are warm and dry, and most of the rain falls in winter months when vegetation is dormant and the temperature is low. Most biogeochemical processes require hospitable temperatures and an adequate water supply. Interactions between the ecosystem services that soils provide, such as nutrient cycling and the regulation of water quality and quantity, are strongly influenced by the seasonality of climate and individual storm events.
Together, (1) the Mediterranean climate, (2) the inferred spatial homogeneity of soils and (3) oak tree-induced differences in soil properties, pose an interesting challenge for understanding the ecosystem services that soils provide and subsequent management implications.
Our primary objective was to examine the linkages between nitrogen cycling and soil hydrology in a manner that considers how the temporal and spatial variability of soil properties govern these processes and ultimately the soil's ability to regulate stream water quantity and quality in oak woodlands. We present results from multiple long-term monitoring projects at the UC Sierra Foothill Research and Extension Center (SFREC) that integrate soil, ecology and water-resources data from the area's California blue oak (Quercus douglasii) and live oak (Q. wislizeni) woodlands. The ability to co-locate these projects is rare and would not be possible without the support of the SFREC infrastructure and staff.
A natural laboratory
SFREC has served as a natural laboratory to study soils, ecology, hydrological processes and land management effects for many years. There is a nearly 30-year record of water flow and water quality for its Schubert watershed (Lewis et al. 2000, 2006). Schubert watershed is 255 acres in size, with tree coverage of approximately 50%, predominantly blue oak. A similar experimental 90-acre watershed at SFREC, called Lewis-1, was instrumented in 2002. These paired watersheds are within 3 miles of each other, and are two of four watersheds at SFREC that we have been monitoring for stream flow and water quality. They were chosen because they are similar in elevation, topography, geology, soils and vegetation (type and canopy coverage). The main difference between Lewis-1 and Schubert is size, Lewis-1 being 90 acres and Schubert, 255 acres. As part of a project beginning in 2002, Lewis-1 has been ungrazed for about 10 years, and Schubert has been heavily grazed.
In the experimental Lewis-1 watershed, we installed a weather station, 400 soil-moisture sensors, a perched water monitoring infrastructure and a stream-flow monitoring network. A perched water table occurs during the rainy season above a relatively impermeable clay-rich horizon situated approximately 10 inches below the soil surface. This infrastructure has allowed us to document the relationships between rainfall, soil moisture storage, lateral flow through soils and stream flow at 15-minute intervals since summer 2006.
Above, the Lewis-1 watershed at SFREC was equipped with instruments in 2002 to sample and monitor the interactions between soil and stream water. Top left, a perched water monitoring trench uses trays to collect water samples from different soil horizons (A, AB, Bt1 and BC). Top right, tipping buckets are mounted on the blue containers to measure lateral water flow.
To understand the temporal and spatial patterns in water content, sensors were placed in 100 soil profiles, typically at depths of 4, 12, 20 and 40 inches. Perched water flow was monitored continuously from a trench that contained three soil profiles fitted with water-collection trays inserted at four depths: 4 inches (bottom of A horizon), 12 inches (bottom of AB horizon), 24 inches (upper boundary of the clay pan) and 43 inches (below the clay pan). Water flow rates were recorded from each layer using tipping buckets, and water chemistry samples were collected in plastic drums. Stream flow was monitored continuously using water-height sensors connected to a combination V-notch weir for measurement at low-flow conditions, and a Parshall flume was used for high-flow measurements. Automated samplers were used to collect water samples from streams on an hourly basis during storm events, and base flow was collected periodically between storms.
The soil and ecological data that we collected reflects a variety of experiments that have been conducted at SFREC over the last two decades, including at Schubert watershed. Soil profiles under the oak canopy and open grasslands were sampled by genetic soil horizon to determine the effects of oak nutrient cycling on soil characteristics. Organic matter concentrations (total carbon and nitrogen) were determined for each soil horizon, and these values were converted to an area basis using the soil bulk density, horizon thickness and rock content. Two representative blue oaks (70 and 92 years old with 11.5-and 19-inch-diameter breast heights, respectively) were destructively sampled, and all aboveground components were dried and weighed to determine the size of the carbon and nutrient pools. Root biomass and nutrient concentrations were calculated from root quantification studies conducted for these same sites (Millikin and Bledsoe 1999). To document temporal changes in soil solution chemistry, 1.5-to-1 water extracts were prepared for each soil horizon monthly to determine the concentrations of soluble nutrients in soil horizons beneath the oak canopy and in adjacent grasslands.
Fig. 1. Nitrogen cycling with major pools of nitrogen (pounds/acre) for an oak woodland-grassland ecosystem in the Schubert watershed at SFREC.
Soils at SFREC and throughout the metavolcanic terrain within the foothill region are moderately deep to very deep, ranging from 20 inches to more than 60 inches. Soils possess a distinctive red color due to the high iron content of the parent material that is released through weathering. A typical soil profile has the following genetic horizon sequence: A-AB-Bt1-Bt2-BC-Cr or R.
The A and AB horizons have a clay-loam texture and appreciable soil organic-matter enrichment. Under oak trees, A horizons are thicker (4 versus 2 inches), with higher organic matter and lower bulk density compared to those in open grasslands (Dahlgren et al. 1997). The A and AB horizons are porous and extensively burrowed by pocket gophers. Soils with Bt1 and Bt2 horizons represent an accumulation of clay (about 35% to 50%, maximum) transported from overlying horizons by percolating water. Bt horizon textures range from clay loam to clay depending on the landscape position. Claypans (Bt2) are a subset of Bt horizons that have an abrupt clay increase (more than 20% relative to the overlying horizon) over a short vertical distance (less than 1 inch). Soil profiles containing claypans occupy less than half of the landscape area in the experimental watershed and tend to occur on level or gently sloping hillsides and land forms that accumulate water. Clay content decreases below the Bt horizons, resulting in a transition horizon (BC). The soil interface with bedrock is described by either Cr (soft bedrock) or R (hard bedrock) horizons.
Nitrogen cycling within the landscape is directly linked with the oak canopy distribution and hydrologic cycle. The majority of nitrogen is present as soil organic matter stored in the A and AB horizons (approximately the upper 15 inches of soil). We found that the major nitrogen pool is contained in soil profiles under oak trees at about 8,373 pounds per acre (fig. 1). The nitrogen pool contained within the oak trees was approximately 14% of the total soil nitrogen pool. Nitrogen stored in blue oak was about 1,138 pounds per acre (931 aboveground plus 207 pounds per acre in roots belowground), while understory components (e.g., grasses and herbs) added 39 and 53 pounds per acre under the oak canopy and in open grasslands, respectively. The soil nitrogen pool was nearly 1.6 times greater (8,373 versus 5,344 pounds per acre) in soils beneath the oak canopy compared to open grasslands (fig. 1).
Fig. 2. Stream flow and nitrate-nitrogen (NO3-N) load in the 1992-1993 water year, Dec. 1 through April 30, at the Schubert watershed.
The ability of oaks to alter their soil environment occurs primarily through the addition of organic matter and nutrient cycling. During a 3-year (1990-1992) study period, oaks returned about 4,008 and 75 pounds per acre annually of organic carbon and nitrogen, respectively, to the soil surface through aboveground litter fall. Belowground carbon inputs through oak and annual grass root turnover also contributed an appreciable amount of organic matter (Jackson et al. 1990; Millikin and Bledsoe 1999). An additional 5.3 pounds per acre of nitrogen was added annually to the soil surface by canopy throughfall (canopy drip) and stem flow (water flowing down the tree trunk). Nutrient flux in the canopy throughfall originates from the capture of atmospheric aerosols, gases and particulate matter, as well as from root uptake with subsequent foliar leaching.
Fig. 3. Nitrate-nitrogen (NO3-N, parts per million) concentrations in soil pore water under oak canopy and open grassland.
In addition, oaks are rooted considerably deeper than annual grasses, which reduces deep leaching and the loss of nutrients beneath the oak canopy. About 70% of the oak root biomass occurs within the upper 20 inches of the soil profile (Millikin and Bledsoe 1999). Oak roots also extend well beyond the edge of the canopy, where they can sequester nutrients from grassland soils. Over time, they concentrate these nutrients beneath the oak canopy in the form of litter fall and canopy throughfall. Given the organic-matter enrichment in soils and oak vegetation, there is considerable potential for carbon sequestration (about 89 tons per acre) associated with the conversion of grasslands to oak woodlands. To ensure long-term carbon sequestration, oak trees must be protected. Organic carbon and many of the enhanced soil properties are lost within 20 to 30 years following an oak's removal (Dahlgren et al. 2003).
Climate and seasons.
The Mediterranean climate contributes to strong spatial and seasonal patterns in soil-solution and stream nutrient concentrations, as illustrated by nitrate concentrations during the winter/spring rainy season (figs. 2 and 3). Peak soil-solution nitrate concentrations occur in fall and late spring and are considerably higher in soils under oak canopy during these times (fig. 3).
In contrast, nitrate concentrations are very low and similar between oak and open-grassland soils during the winter months when soil moisture storage and leaching are at their peaks, and soil temperatures and microbial activity are low.
In addition, annual grass growth and oak foliage production results in strong demand for biological nutrients in spring, further depleting the soil nitrate pool in spite of favorable moisture and temperature conditions for microbial mineralization. The increase in soil nitrate in late spring coincides with the senescence (die off) of annual grasses, which results in a decrease in nitrogen uptake demands coupled with the release of nitrogen from dead annual grass roots by microbial decomposition. Dry conditions in the upper soil profile throughout summer result in low nitrogen uptake by oak roots.
Instead of continuous nitrogen feedback among senescing plants, their soils and new growth, nitrogen is mineralized and accumulates in soils during the dry summer and fall months (Hart et al. 1993). Throughout summer, organic and inorganic nitrogen accumulates in the soil, and there is little plant growth to sequester these nutrients. Early fall rains (September to October) cause the leaching of ammonium, nitrate and dissolved organic nitrogen (DON) from annual grass residues (more than 50% of total DON) while senescent oak foliage accounts for more than 70% of total DON after the first storm (Chow et al. 2009). High concentrations of DON combined with water availability and warm residual soil temperatures lead to an explosion of microbial activity that results in rapid nitrogen mineralization and the accumulation of nitrate in the upper soil horizons.
An asynchrony within nutrient cycles occurs in the Mediterranean climate, causing a marked nitrate spike in stream water during the onset of the rainy season (December to January) (Holloway and Dahlgren 2001). Nitrate leaching from oak woodland watersheds at SFREC is strongly linked to the hydrologic cycle at time scales ranging from storm events (hours) to years. Pulses in stream nitrate concentrations are observed that match those observed in soil pore water in the fall. Each storm progressively flushes this nitrogen pool, so that by March there is little nitrate found in soils and stream water (figs. 2 and 3).
Hydrologic flow paths
Stream flow and water quality are related to soil water storage and hydrologic flow paths through soils. Our conceptual model of the dominant hydrologic flow path in soils of this area is lateral, perched water flow through the upper soil layers (A and AB horizons) following watershed priming. Approximately 7 to 10 inches of precipitation are required to bring the watershed soils to their maximum water-holding capacity before stream flow generation begins (Lewis et al. 2000). With the onset of winter storms, water infiltrates and percolates through the highly permeable A and AB horizons. The vertical percolation of water is impeded by the presence of the claypan approximately 24 inches below the soil surface. The claypan is slowly permeable to water, and as a result saturated conditions develop at the upper boundary of the claypan during storms. A saturation zone (perched water table) forms at the claypan boundary and extends upward into the A and AB horizons as precipitation continues. Once they are saturated, water moves with the force of gravity (down slope) through the highly permeable A and AB horizons, rapidly delivering water as lateral flow to streams.
The hydrologic monitoring infrastructure within the experimental watershed demonstrated the connectivity of lateral flow in soils with stream flow during a typical storm (fig. 4). The timing of stream and perched water flow during the course of a storm was similar. Peak stream flow was sustained for approximately 3 hours from 7 p.m. to 10 p.m. during a storm in January 2007 and corresponded with the initiation of lateral flow through the AB horizon located above the claypan. Perched water flow through the AB horizon was the dominant flow path. During this storm event the AB horizon delivered 60% of the total water volume, measured from the perched water monitoring infrastructure. The other horizons contributed relatively equal water volumes ranging from 10% to 17% (fig. 5). An additional line of evidence is apparent from the slope characteristics associated with the rising limb of the hydrographs, where the stream and AB horizon have nearly identical slopes, suggesting that the rapid increase in stream flow is a result of water transport through the AB horizon (fig. 4).
Fig. 4. (A) Stream and (B) perched water hydrographs during the first stream flow event of 2007.
The AB horizon supplied more water to the stream for a variety of reasons. The AB horizon is more permeable than the Bt and BC horizons, because it contains less clay, better soil structure and more large pores (e.g., gopher burrows and decayed root channels) to conduct water rapidly. While the permeability of the AB horizon is similar to that of the A horizon, the AB horizon supplied more water because it was saturated longer due to its proximity to the underlying water-restrictive clay layer. The delay in flow of the A relative to the AB horizon represents the time it took for the perched water table to extend upward into the A horizon. The AB horizon is also thicker than the A horizon (8 versus 4 inches), resulting in a greater cross-sectional area to deliver water.
Fig. 5. Contribution of hydrologic flow paths in soils at SFREC. Percentages refer to the total lateral flow contributed by each soil layer during a typical storm event.
The hydrologic flow paths through soils at SFREC directly affect water quality. With their permeable nature, surface horizons convey water rapidly to streams in large volumes. Infiltration rates range from 1 to 4 inches per hour, and as a result, surface runoff is rare when soils are not saturated (Lewis et al. 2000). Because plants are not actively transpiring during the rainy season, the soil water-storage capacity is frequently exceeded, which results in the transport of water via lateral flow through the nutrient-rich A and AB horizons, directly to the stream. The nitrate accumulated in the A and AB horizons is susceptible to leaching because little plant uptake occurs during the winter months. The shift to laterally dominated flow paths during storms results in elevated nitrogen concentrations in stream water as nitrate, which is leached directly to the stream via lateral subsurface flow through the A and AB horizons (figs. 4, 5 and 6). In contrast, elements that are associated with groundwater and chemical weathering in the lower Bt and BC horizons (e.g., calcium and sodium) display a distinct decrease in concentration during storms as the hydrologic flow path is short-circuited through the upper soil horizons (fig. 6).
Fig. 6. Changes in stream water chemistry during an early-season (Jan. 2007) storm event.
Fig. 7. Relationship of stream nitrate fluxes and stream flow volume (water flux) for a 20-year period (1981-2000) in the Schubert watershed at SFREC; r2 = 0.52; P < 0.001.
There is considerable year-to-year variability in nitrate leaching from oak woodland watersheds (fig. 7). The mean annual nitrate-nitrogen (NO3-N) export from the Schubert watershed was 1.4 pounds per acre annually (ranging from 0.16 to 3.2 pounds per acre annually) (Lewis et al. 2006). while annual nitrate-nitrogen flux at Schubert reached values as high as 3.2 pounds per acre annually, the loss of nitrate-nitrogen from the watershed was substantially less than atmospheric inputs of mineral nitrogen (ammonium [NH4] plus nitrate [NO3]) in the bulk precipitation (mean 8.4 pounds per acre annually). As a result, there was a net retention of mineral nitrogen within the Schubert watershed, assuming no loss of nitrogen through denitrification or volatilization. The annual nitrate-nitrogen load was positively and significantly (P < 0.05) related to both annual precipitation and stream flow (Lewis et al. 2006). Nitrate fluxes were a function of total stream flow volume and the timing of storm-flow events throughout the season. Years with high rainfall and runoff, combined with high runoff during the early portion of the rainy season — when nitrate concentrations are highest in the soil profile — resulted in the largest fluxes of nitrate-nitrogen from the watershed.
Differences in stream nitrate flux (figs. 2 and 6) can also be explained by the extent of the perched water table within the watershed and its connectivity to the stream. The observed variation in stream nitrate flux is a result of flushing of the A and AB horizons by the lateral flow path, because these horizons contain most of the nitrogen. The degree of flushing is a function of the amount and intensity of rainfall, the seasonal distribution of rainfall, the spatial extent of the claypan, and the soil moisture content prior to a storm event.
Perched water table.
Figure 8 shows the extent and connectivity of the perched water table within the AB horizon as indicated by saturated conditions (water content greater than 35%) during peak stream flow for two storms: (A) Jan. 4, with a large stream nitrate concentration and (B) Feb. 27, with a low stream nitrate concentration. The spatial extent of the perched water table was 59% of the watershed area during the first storm (fig. 8A) compared to 36% during the second storm (fig. 8B). The perched water table was connected to the stream, as indicated by the overlap of regions shaded in green with blue flow lines, which were generated from a digital elevation model to predict the routing of water by topography, using the slope and curvature of the landscape (fig. 8).
Fig. 8. Extent of saturated conditions in the experimental watershed within the AB horizon during peak flow of two storm events: (A) Jan. 4 and (B) Feb. 27, 2008. Saturated conditions are depicted in green where water content is 35% ± 3%.
The perched water table during the second storm occupied isolated areas with less connectivity to the stream, as indicated by areas of yellow (unsaturated) and green (saturated) intersecting the blue flow lines (fig. 8B). This suggests that the dominant source area of stream flow was limited to the center of the watershed during the Feb. 24 storm, when saturated conditions coincided with the main flow lines. Thus, variation in nitrate flux observed over a season (figs. 2 and 3) and during a storm (fig. 6) is a result of the sequential flushing of nitrate from the soil over time and the area of the watershed that is leached, determined by the spatial extent of the perched water table and its connectivity to streams.
Soils and rangeland management
More than two-thirds of California's drinking-water supply passes through or is stored in oak woodlands. Despite the seasonal peaks in nitrate in streams at SFREC, the concentrations are low relative to drinking-water standards (10 parts per million [ppm] nitrate-nitrogen). However, the temporal and spatial patterns between soils and streams observed in this study are relevant to understanding how soils affect water supply and its quality, and can help to guide rangeland management practices such as the stocking rate, season and location of livestock congregation activities (supplemental feeding stations).
Two-thirds of California's drinking water passes through or is stored in oak woodlands. Detailed knowledge about soil properties can help rangeland managers to assess the impacts of grazing on water supply and quality.
In some instances, the variability of soil properties at SFREC and in the Sierra Foothill region exceeds our ability to document them at scales relevant to rangeland management. Key soil properties, however, such as the distribution of claypans, should be documented by soil surveys in order to understand the spatial connectivity of hydrologic-transport processes (surface runoff, subsurface flow, stream flow) within watersheds. Currently, California's soil surveys in rangelands are not mapped at scales fine enough to portray these soil features across landforms and their connectivity to streams. It is important that stakeholders communicate this concern to the U.S. National Cooperative Soil Survey to focus soil survey updates in these important landscapes. This level of detail is warranted to better understand patterns in oak regeneration, perennial-grass restoration, forage productivity, rainfall-to-runoff characteristics and water-quality dynamics.