Published online 14 July 2009
Published in Soil Sci Soc Am J 73:1504-1509 (2009)
DOI: 10.2136/sssaj2008.0106
© 2009 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
PEDOLOGY
Soil Carbon and Nitrogen Accumulation in a Forested Debris Flow Chronosequence, California
Judith K. Turk* and
Robert C. Graham
Soil and Water Sciences Program, Dep. of Environmental Sciences, Univ. of California, Riverside, CA 92521
* Corresponding author (judith.turk{at}email.ucr.edu).
 |
ABSTRACT
|
|---|
The role of forest soils in the biogeochemical cycling of C and N is most dynamic during the early stages of soil development. To define C and N trends that occur with soil development in a mixed coniferous forest, a chronosequence formed by debris flows was studied. The accumulation rates of total organic C (TOC) and total N (TN) were evaluated in soils on 10 debris flow deposits, ranging from <1 to 244 yr old. Analysis of the mineral soils was restricted to the 30-cm depth, since this was the depth of the shallowest debris flows. Carbon was found to accumulate in the organic horizons at a rate of 26.5 g m–2 yr–1 throughout the time span of the chronosequence. Total organic C accumulation in the mineral horizons (0–30 cm) occurred from 0 to 82 yr at a rate of 13 g m–2 yr–1, and was nearly stable from 82 to 244 yr. Total N accumulated at a rate of 0.57 g m–2 yr–1 in the organic horizons and a rate of 0.17 g m–2 yr–1 in the mineral horizons (0–30 cm) throughout the 244 yr chronosequence. This study suggests that C accumulation in the upper mineral horizons of young forest soils occurs for <100 yr, while N accumulation is a slower process that occurs for >250 yr. Carbon and N accumulation in the organic horizons, however, both follow a linear trend over the 244-yr period. The rates of C accumulation suggest a rapid recovery of the soil organic C pool following disturbance.
Abbreviations: SDEF, San Dimas Experimental Forest SMRNA, Shasta Mudflow Research Natural Area TN, total N TOC, total organic C
 |
INTRODUCTION
|
|---|
Soil organic matter is a unique soil property in that it reaches an approximate steady state in <5000 to 20,000 yr, which is rapid compared with other soil properties (Birkeland, 1999). In a review of chronosequence studies in a variety of ecosystems, Schlesinger (1990) found a much higher and more variable rate of C accumulation in soils <3 ka in age (5.7 ± 1.8 m–2 yr–1) compared with soils ranging from 3 to 10 ka (2.4 ± 0.7 m–2 yr–1). Thus, when considering the development of soil organic matter pools, it is important to include studies focusing on the early stages of soil development because this is when the most rapid accumulation takes place.
The organic C reservoir in soils is estimated to be 1580 Pg, which is greater than the C pool in land plant biomass and the atmosphere combined (estimated to be 1360 Pg) (Schimel et al., 1995). In particular, forest soils are considered to be an important sink for organic C. Coniferous forests cover 3.4% of the Earth's land area (Whittaker and Likens, 1975) and have an average soil organic C content of 15.1 kg m–2 (Kern, 1994). Using these values, it can be estimated that there are 75.5 Pg of C in coniferous forest soils, which is 5% of the global C pool in soils.
In the United States, 40% of forests occur in zones designated as landslide hazard areas (Fig. 1
). Thus, the ecosystem recovery and development of soils following disturbance by mass wasting is a natural component of the biogeochemical cycling of forests. In southern California, debris flows are a common type of mass wasting, occurring in several mountainous regions (Morton et al., 2008).
In the current study, debris flow deposits along a particularly active channel were studied as a high-resolution chronosequence, including 10 debris flows ranging in age from <1 to 244 yr. This approach was useful for focusing on soil formation during the phase in which C and N accumulate most rapidly. The study was conducted to resolve trends in TOC and TN which may not be apparent when longer time frames of soil development are considered.
 |
MATERIALS AND METHODS
|
|---|
Site Description
The study site is near the town of Forest Falls in the San Bernardino Mountains of southern California (Fig. 1). The elevation of the site is approximately 1675 m and an adjacent ridge (Yucaipa Ridge) rises approximately 500 m above the study area. Yucaipa Ridge is composed mainly of highly fractured gneiss, which produces debris flows during short-duration, high-intensity rainfall events (Morton et al., 2008). The debris flows leave discontinuous levee deposits and lobate deposits with steep-fronted snouts in the forested areas surrounding the debris flow channels. Debris flows on Yucaipa Ridge occur approximately once every 5 yr in three primary debris flow channels (Morton et al., 2008). The channels coalesce at the foot of the ridge to form an approximately 30-ha fan.
The deposits used for this study are on the western bank of one of the debris flow channels, Rattlesnake Creek (N 34.1° W 116.9°, Datum WGS84). The ages of the deposits range from <1 to 244 yr old, as determined using dendrochronological techniques (Turk et al., 2008). The deposits have a north-facing aspect, 9 to 22% slope gradient, and a convex/convex slope shape. Individual deposits range from 70 to 900 m2 in area and from 0.3 to 2 m in thickness. The deposits are composed of loamy coarse sand and unsorted rocks (36–65% by volume), which range from gravel to stone-size (2–600 mm). Some woody debris (e.g., branches) is also incorporated into the initial deposits. Very little textural differentiation occurred between the horizons with time in the fine-earth fraction, though there was a decrease in rock fragments of the surface horizons over the course of the chronosequence (Turk et al., 2008). Soils at the site were classified as either sandy-skeletal, mixed, mesic Typic Xerorthents or sandy-skeletal, mixed, mesic Typic Xerofluvents (Turk et al., 2008).
The forest is composed predominantly of incense cedar (Calocedrus decurrens Torr. Florin), with interspersed ponderosa pine (Pinus ponderosa P. Laws. & C. Laws.), Coulter pine (Pinus coulteri D. Don), California black oak (Quercus kelloggii Newb.), white fir (Abies concolor [Gord. and Glend.] Lindl. ex Hildebr.), and canyon live oak (Quercus chrysolepis Liebm.).
The mean annual precipitation at Forest Falls is 855 mm, with most of the precipitation occurring between November and March, as rain and snow (San Bernardino County Flood Control District, 2008). The mean annual temperature is 12.3°C, as recorded at a nearby weather station (Camp Angelus, 34.15°N 116.98°W, 1762 m elevation) (San Bernardino County Flood Control District, 2008).
Soil Sampling
One soil pit was excavated to the 0.7-m depth on each debris flow. Given the small area of the individual debris flows (70–900 m2), analysis of more than one soil pit per debris flow deposit was not justified, as it would have lead to pseudoreplication and excessive disturbance to the site. The location of the soil pit was chosen in areas that were least influenced by rocky debris flow snouts and erosion processes. The soil pits were described and the <75-mm fraction of each horizon was sampled for laboratory analysis. Site A was sampled in more than 1 yr, using separate pits (Table 1
). Each organic horizon was sampled by removing all material from the horizon within the 189-cm2 area outlined by a 15.5-cm diameter metal cylinder placed on the soil surface. This was done in triplicate and the three samples were combined. Most organic horizon samples were taken during the summer, so that calculations would be representative of the yearly minimum forest litter. However, the Oi-horizon sample of the 244-yr-old soil was taken in the fall of 2003, before removal of the 244-yr-old cohort of trees by logging.
Rock fragments >75 mm in length were separated from soil material during excavation of the soil pits and weighed in the field. To determine volume of the soil pits, plastic bags were stuffed with known volumes of Styrofoam peanuts (10–100 dm3) and molded into the pit until it was filled. The volume of bags that could fit in the pit was used to approximate the volume of the pit. Extra room was left around the peanuts in the bags so that they could be moved within the bags and the bags could conform to the edges of the pit and the surrounding bags. Small holes were poked in the bags to prevent trapped air from affecting the volume.
Laboratory Methods
Organic horizon samples were fractionated by hand sorting. The fractions included needles and leaves, woody structures (twigs, bark, and roots), reproductive structures (pine cones, catkins, and acorns), and decomposed unknowns (those components whose botanical origins were not recognizable). Each of the fractions was dried at 65°C, weighed, and ground with a Wiley mill for C and N analysis.
Mineral soil samples were air-dried and sieved to remove gravels (2–75 mm). The fine-earth (<2-mm) and gravel fractions were weighed. The lab weights and field measurements of rock weight and pit volume (described in the previous section) were combined to determine volume fraction of rock fragment in each horizon (assuming an even distribution of >75 mm rock fragments throughout the profile) as follows:
Where, Xv-cfi is the volume fraction of coarse fragments in horizon i, Wcb+ is the weight of the >75-mm fragments in the pit (in kg), Vp is the volume of the pit (in m3),
p is the particle density (in kg m–3), Wgri is the weight of gravels in the sample from horizon i (in kg), Wfei is the weight of the fine earth fraction in the sample from horizon i, and
bi is the bulk density of the fine-earth fraction of horizon i.
The particle density of the parent material was 2.7 x 103 kg m–3, as measured by the gas pycnometer method (Blake and Hartge, 1986b), using an AccuPyc 1330 pycnometer (Micromeritics, Norcross, GA). Bulk density of the organic horizons was calculated based on the dry weight and core volume of the field samples. Bulk density of the mineral soil horizons was measured using the clod method (Blake and Hartge, 1986a) with corrections for gravel (Hirmas and Furquim, 2006).
The fine-earth fraction was ground to pass through a 100-mesh (149 µm) sieve for total C and N analysis and carbonate analysis. Total C and N were determined by dry combustion using a Thermo-Finnigan Flash EA112 soil combustion N and C analyzer system (CE Elantech, Lakewood, NJ). Carbonate C, which was present in the <1- and 3-yr-old soils, was measured by the manometer method (Williams, 1948), using a Dwyer series 477 digital-manometer (Dwyer Instruments, Inc., Michigan City, IN). Total organic C (TOC) was determined by subtracting carbonate C from total C.
Concentrations were converted to mass per area (in g m–2) as follows:
(Boone et al., 1999). Where M is mass (in g), A is area (in m2), n is the number of horizons in the soil, Ti is the thickness of horizon i (in m),
bi is the bulk density of horizon i (in kg m–3), Ci is the concentration of C or N in horizon i (in g kg–1), and Xv-cfi is the volume fraction of coarse fragments in horizon i. Mass of TOC and TN per area were calculated for the organic horizons and for the upper 30 cm of the mineral soil. The depth of 30 cm was chosen because this was the shallowest depth at which buried A horizons were encountered. We wanted to assess TOC and TN accumulation in fresh parent material, so buried horizons needed to be excluded.
 |
RESULTS AND DISCUSSION
|
|---|
Organic Carbon Accumulation
Most of the C accumulation in the soils formed at the Forest Falls debris flows occurred in the organic horizons. The rate of C accumulation in the organic horizons was determined to be 26.5 g m–2 yr–1 (Fig. 2
). Accumulation of C in the organic horizons occurred predominantly in two fractions: the decomposed fraction (decomposed beyond recognition of botanical origins) and the woody materials (twigs, bark, and roots). The rates of C accumulation in these fractions (18.8 g m–2 yr for decomposed fraction and 6.1 g m–2 yr for woody materials) nearly account for the total C accumulation rate in the organic horizons. Other organic horizon fractions (needles and leaves and reproductive structures) were variable in their contribution to the C pool, but did not display a significant trend in with time. Carbon originating from needles and leaves was apparently accumulated in the decomposed fractions, reflecting the rapid decomposition of these structures (Swift et al., 1979).

View larger version (18K):
[in this window]
[in a new window]
|
Fig. 2. Total organic carbon (TOC) in the organic horizons (all fractions) and separate components of the organic horizons (highly decomposed and woody fraction).
|
|
Organic C accumulation in the mineral soil is slower and shorter-lived than in the organic horizons. In the mineral soil, organic C accumulation occurs from <1 to 82 yr and the rate of accumulation during this stage is 13 g m–2 yr–1 (Fig. 3
). From 82 to 244 yr the C pool in the mineral soil is nearly stable, there is a slight decline in the size of the TOC pool, but the trend is not significant (Fig. 3). These trends are similar to those observed in the chronosequence at Shasta Mudflow Research Natural Area (SMRNA), located in the Cascade Range of northern California. Both sites show the most rapid C accumulation in soils <100 yr old (Dickson and Crocker, 1953b).
To interpret the C accumulation rates calculated from this chronosequence in a broader context, it is important to recognize the limitations and assumptions of the calculated chronofunctions. In most of the debris flows (Sites A, B, D, F, G, and I) buried A horizons were observed within 100 cm of the soil surface. Although not observed, the other sites are assumed to have buried A horizons within 2 m of the surface, based on the maximum thickness of individual deposits observed in cross-sections exposed in the debris flow channels. To calculate the amount of C accumulated since the most recent debris flow, it was necessary to exclude the C in the buried soil, which predates the debris flows at the surface. Therefore, the C pool reported here only includes the depth to the shallowest buried A horizons, which were at 30 cm. This calculation leaves out a significant portion of the C pool that accumulates deeper in the soils, through 100-cm depth, due to root growth (Turk et al., 2008). At sites where a buried A horizon was not encountered within 100 cm of the mineral soil surface (Sites E and J) the upper 30 cm of the soil accounted for only 50% of the TOC stored in the mineral soil. Thus, it cannot reasonably be assumed that most C from root deposition and decay accumulates in the upper 30 cm, but by necessity, the chronofunctions presented here exclude this deeper C.
The shallow buried A horizons and the small area of the debris flows (70–900 m2) probably themselves affect ecosystem development, further limiting the applicability of the chronofunctions to soil development on larger geologic deposits. Since the developing forest grows both in the fresh debris flow deposit and in what may be a well-developed A horizon as little as 30 cm beneath the surface, the debris flows are somewhat intermediate between a primary and secondary successional chronsequence. There are also two important effects of the small area of the individual debris flows (i) colonization from surrounding vegetated areas is almost immediate and (ii) much of the litter deposited on a fresh debris flow surface originates from vegetation on older debris flows adjacent to the young surface. These effects were illustrated on the 3-yr-old debris flow deposit, which had a dense cover of incense cedar seedlings (150 m–2) (Turk et al., 2008) and a nearly continuous, 1-cm-thick layer of litter from mature trees on the surrounding debris flows.
In spite of these limitations, some comparisons can be made with other studies of the early stages of soil formation under coniferous vegetation in California. Related studies include the lysimeter biosequence at the San Dimas Experimental Forest (SDEF) in the San Gabriel Mountains in southern California (Ulery et al., 1995) and the chronosequence at SMRNA (Dickson and Crocker, 1953a) (Fig. 1). Comparison between Forest Falls, SMRNA, and SDEF reveals differences in the profile distributions and rates of C accumulation (Table 2
). Forest Falls has a distinctly different profile distribution of C compared with the other sites. At Forest Falls, C accumulates rapidly in the organic horizons and at a lower rate in the mineral soil, however, the inverse is true at SMRNA and SDEF. This trend may be caused by trees growing on older debris flows which contribute leaf litter to the younger debris flows, but do not contribute root biomass. Another factor may be the influence of poor air quality in the San Bernardino Mountains on the tree growth at Forest Falls. Ozone-sensitive conifers exposed to high levels of ozone and atmospheric N deposition have high shoot/root ratio and high rates of foliar turnover, leading to a thick forest floor (Takemoto et al., 2001).
View this table:
[in this window]
[in a new window]
|
Table 2. Climate conditions, vegetation, and rates of carbon and nitrogen accumulation calculated for short-term (<300 yr) soil formation over under coniferous vegetation in California.
|
|
Comparing the sum of C accumulation rate in the organic and mineral horizons (Table 2), the data show lower rates of accumulation at the warmer, drier sites (Forest Falls and SDEF) compared to the cooler, wetter SMRNA chronosequence. These data support the general observation of greater soil C pools as temperature decreases and precipitation increases (Post et al., 1982). Overall, the rates of C accumulation at all three young, forested sites (39.3–139.8 g m–2 yr–1) (Table 2) are high in comparison with those observed in chronosequence studies involving older landscapes and other ecosystems in California. Rates of 0.26 to 5.4 g m–2 yr–1 for soil organic C through the 30-cm depth can be calculated from the earliest stages of soil development considered in chronosequence studies of grasslands and coastal scrub (Harden, 1987, Harden et al., 1987). This difference in C accumulation rates can be attributed partly to differences in the C balance of these ecosystems. Generally higher C storage is observed in forest soils compared with grasslands and scrublands (Kern, 1994). However, it is also important to consider the time span of the chronosequence. In the grassland and coastal shrub chronosequences the difference in age between the youngest two surfaces is >1000 yr (Harden, 1987, Harden et al., 1987). These early years of soil development may include trend reversals and changes in the C accumulation trend, as observed in the mineral soil horizons at Forest Falls and Mt. Shasta (Dickson and Crocker, 1953b) and an overall trend toward a slower accumulation rate with time, as described by Schlesinger (1990).
Given a longer time span at Forest Falls, it is expected that the linear trend in C accumulation of the organic horizons would approach an approximate steady state level. The oldest site observed in the study (244 yr old) stored about 9 kg m–2 of C through the 1-m depth (starting from the surface of the organic horizons), which is lower than the average C content for soils in temperate coniferous forest (15.1 kg m–2) (Kern, 1994). A prediction of when steady state will be reached can be made by extrapolating the trends from this chronosequence (26.5 g m–2 accumulation in the organic horizons and no further accumulation in the mineral soil) and assuming that the steady state level of C is around 15.1 kg m–2. If the assumptions of this prediction hold, it is estimated that the C pool stabilizes 230 yr beyond the last age represented in the chronosequence, or within about 500 yr of soil development. This rapid achievement of steady state was also observed in the SMRNA chronosequence, in which the soil C pool was observed to be nearly stabile after 600 yr (Dickson and Crocker, 1953b, Lilienfein et al., 2003). This represents a much more rapid stabilization of the soil C pool than the general observation that soils reach steady-state organic C levels in <5000 to 20,000 yr (Birkeland, 1999). These results indicate that the C pool of soils of the mixed coniferous forests in California recover rapidly from disturbance by mass wasting.
Nitrogen Accumulation
Nitrogen accumulation in the organic horizons occurred at a rate of 0.63 g m–2 yr–1 (Fig. 4
). As observed for C storage, the N accumulation occurred primarily in two fractions of the organic horizon: the decomposed fraction and the woody materials. Nitrogen accumulated in the decomposed fraction at a rate of 0.49 g m–2 yr–1 and in the woody materials at a rate of 0.08 g m–2 yr–1, thus accounting for most of the N accumulation rate in the organic horizons. The rate of N accumulation in the mineral horizons was 0.17 g m–2 yr, but exhibits a high level of variability (Fig. 5
). The upper 30 cm accounts for 52% of the N storage through the 100-cm depth in the 244-yr-old debris flow, revealing that the upper 30 cm provides an underestimate of total N stored in the mineral horizons of the soils.

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 4. Total nitrogen (TN) in the organic horizons (all fractions) and separate components of the organic horizons (highly decomposed and woody fraction).
|
|
The total rate of N accumulation in the Forest Falls chronosequence (0.80 g m–2 yr–1) is lower than the rates of N accumulation determined for both the SMRNA (2.33 g m–2 yr–1) and the SDEF (1.9 g m–2 yr–1) (all calculated through 30-cm depth) (Table 2). The higher N accumulation rate at SMRNA could be related to the abundance of plant species that are symbiotic N-fixers. Nitrogen-fixing species at SMRNA (Ceanothus spp. and Purshia tridentata Pursh.) were documented on nearly every mudflow surface (Dickson and Crocker, 1953a). However, at the Forest Falls chronosequence there was only one symbiotic N-fixing plant species (Cercocarpus ledifolius Nutt.), which was observed on the 20 yr-old surface only (Turk, 2007). Likewise, at SDEF, the Coulter pine lysimeter was not planted with any N-fixing species.
At Forest Falls and SDEF, atmospheric deposition is likely to be an important source of N. The N accumulation rate at Forest Falls (0.8 g m–2 yr–1) (Table 2) is within the range of dry deposition rates of N at Barton Flats (0.5–0.9 g m–2 yr–1) (Bytnerowicz and Fenn, 1996), which is five miles north of Forest Falls. Likewise, the N accumulation rate observed at the SDEF lysimeters (1.9 g m–2 yr–1) (Table 2) could be entirely accounted for by the high the N deposition rate at the site (2.0–3.5 g m–2 yr–1) (Bytnerowicz and Fenn, 1996). The distribution of N between organic and mineral horizons is subject to the same influences discussed for C accumulation, including surface litter originating from adjacent older surfaces and increased litterfall due to ozone damage (Takemoto et al., 2001).
 |
SUMMARY AND CONCLUSIONS
|
|---|
Accumulation of C and N in the soils at Forest Falls occurs mainly in the organic horizons. Carbon accumulation in the organic horizons can be described by a linear trend with a slope of 26.5 g m–2 yr–1, which is uniform throughout the time frame of the chronosequence. Nitrogen accumulation in the organic horizons also follows a linear trend, with a slope of 0.63 g m–2 yr–1. Most of the C and N accumulation in the organic horizons occurred in fractions that were highly decomposed or composed of woody materials.
In the mineral horizons, C accumulation occurred from <1 to 82 yr, at a rate of 13 g m–2 yr–1 and was nearly stable from 82 to 244 yr. The rate of N accumulation was best described by a linear trend with a slope of 0.17 g m–2 yr–1 and there was no change in trend over the course of the chronosequence. These rates reflect accumulation in the upper 30 cm of the mineral soils only. Limiting the analysis to this depth was necessary to include the shallow debris flow deposits so that enough data was available to resolve trends, such as the change in C accumulation rate at 82 yr. However, another implication is that C and N stored deeper in the soil have not been included.
By the end of this 244-yr chronosequence, the C pool, including the organic horizons and the mineral soil through the 1-m depth, was 9 kg m–2. This approaches the size expected in a mature coniferous forest. This study, therefore, suggests that there is a rapid recovery of the soil organic C pool following disturbance by debris flows. The N, however, accumulates very slowly in the soils and shows no sign of stabilization within the recurrence interval between debris flows at the site.
 |
ACKNOWLEDGMENTS
|
|---|
This research was supported by the Kearney Foundation for Soil Science. Doug Morton introduced us to the field area. Susan Akers, Jacob Berkowitz, Sheila Furquim, Brett Goforth, Daniel Hirmas, Katherine Kendrick, and Rodrigo Vargas were involved in the preliminary research at the site. Joshua Barraza and Shahriar N. Uddin helped with field work.
 |
NOTES
|
|---|
All rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher. Permission for printing and for reprinting the material contained herein has been obtained by the publisher.
Received for publication March 29, 2008.
 |
REFERENCES
|
|---|
- Birkeland, P.W. 1999. Soils and geomorphology. 3rd ed. Oxford Univ. Press, New York.
- Blake, G.R., and K.H. Hartge. 1986a. Bulk density. p. 363–376. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Blake, G.R., and K.H. Hartge. 1986b. Particle density. p. 377–383. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Boone, R.D., D.F. Grigal, P. Sollins, R.J. Ahrens, and D.E. Armstrong. 1999. Soil sampling, preparation, archiving, and quality control. p. 3–28. In G. Philip (ed.) Standard soil methods for long term ecological research. Oxford Univ. Press, New York.
- Bytnerowicz, A., and M.E. Fenn. 1996. Nitrogen deposition in California forests: A review. Environ. Pollut. 92:127–146.[CrossRef][Medline]
- Dickson, B.A., and R.L. Crocker. 1953a. A chronosequence of soils and vegetation near Mt. Shasta, California. 1. Definition of the ecosystem investigated and features of the plant succession. J. Soil Sci. 4:123–141.[CrossRef]
- Dickson, B.A., and R.L. Crocker. 1953b. A chronosequence of soils and vegetation near Mt. Shasta, California. 2. The development of the forest floors and the carbon and nitrogen profiles of the soil. J. Soil Sci. 4:143–153.
- Dickson, B.A., and R.L. Crocker. 1954. A chronosequence of soils and vegetation near Mt. Shasta, California. 3. Some properties of the mineral soils. J. Soil Sci. 5:173–191.[CrossRef]
- Godt, J.W. 2001. Landslide incidence and susceptibility in the conterminous United States. Available at http://www.nationalatlas.gov/mld/lsoverp.html (Verified 18 June 2009). USGS, Reston, VA.
- Harden, J.W. 1987. Soil development in granitic alluvium near Merced, California. U.S. Geol. Surv. Bull. 1590-A. U.S. Gov. Print. Office, Washington, DC.
- Harden, J.W., A.M. Sarna-Wojcicki, and G.R. Dembroff. 1987. Soils developed on coastal and fluvial terraces near Ventura, California: U.S. Geol. Surv. Bull. 1590-B. U.S. Gov. Print. Office, Washington, DC.
- Hirmas, D.R., and S.A.C. Furquim. 2006. Simple modification of the clod method for determining bulk density of very gravelly soils. Commun. Soil Sci. Plant Anal. 37:899–906.[CrossRef][Web of Science]
- Kern, J.S. 1994. Spatial patterns of soil organic-carbon in the contiguous United-States. Soil Sci. Soc. Am. J. 58:439–455.[Abstract/Free Full Text]
- Lilienfein, J., R.G. Qualls, S.M. Uselman, and S.D. Bridgham. 2003. Soil formation and organic matter accretion in a young andesitic chronosequence at Mt. Shasta, California. Geoderma 116:249–264.
- Morton, D.M., R.M. Alvarez, K.R. Ruppert, and B. Goforth. 2008. Contrasting rainfall generated debris flows from adjacent watersheds at Forest Falls, southern California, USA. Geomorphology 96:322–338.[CrossRef][Web of Science]
- Post, W.M., W.R. Emanuel, P.J. Zinke, and A.G. Stangenberger. 1982. Soil carbon pools and world life zones. Nature 298:156–159.[CrossRef]
- Quideau, S.A., R.C. Graham, O.A. Chadwick, and H.B. Wood. 1998. Organic carbon sequestration under chaparral and pine after four decades of soil development. Geoderma 83:227–242.[CrossRef][Web of Science]
- San Bernardino Flood Control District. 2008. San Bernardino County Flood Warning Systme. Available at: http://www.sbcounty.gov/trnsprtn/pwg/alert/reports.html. (verified 1 June 2009)
- Schimel, D., I.G. Enting, M. Heimann, T.M.L. Wigley, D. Raynaud, D. Alves, and U. Siegenthaler. 1995. CO2 and the carbon cycle. p. 35–71. In J.H. Houghton et al (ed.) Climate change 1994. Cambridge Univ. Press, Cambridge.
- Schlesinger, W.H. 1990. Evidence from chronosequence studies for a low carbon-storage potential of Soils. Nature 348:232–234.[CrossRef]
- Swift, M.J., O.W. Heal, and J.M. Anderson. 1979. Decomposition in terrestrial ecosystems University of California Press, Berkeley.
- Takemoto, B.K., A. Bytnerowicz, and M.E. Fenn. 2001. Current and future effects of ozone and atmospheric nitrogen deposition on California's mixed conifer forests. For. Ecol. Manage. 144:159–173.[CrossRef]
- Turk, J.K. 2007. Soil development in a forested debris flow chronosequence in southern California. M.S. thesis. Univ. of California, Riverside.
- Turk, J.K., B.R. Goforth, R.C. Graham, and K.J. Kendrick. 2008. Soil morphology of a debris flow chronosequence in a coniferous forest, southern California. Geoderma 146:157–165.[CrossRef][Web of Science]
- Ulery, A.L., R.C. Graham, O.A. Chadwick, and H.B. Wood. 1995. Decade-scale changes of soil carbon, nitrogen and exchangeable cations under chaparral and pine. Geoderma 65:121–134.[CrossRef][Web of Science]
- USDA-Forest Service. 2002. Forest cover types. Available at http://nationalatlas.gov/mld/foresti.html (Verified 1 June 2009). National Atlas of the United States, Reston, VA.
- Whittaker, R.H., and G.E. Likens. 1975. The biosphere and man. p. 305–328. In H. Lieth and R.H. Whittaker (ed.) Primary productivity of the biosphere. Ecological studies 14. Springer-Verlag, New York.
- Williams, D.E. 1948. A rapid manometric method for the determination of carbonate in soils. Soil Sci. Soc. Am. Proc. 13:127–129.