Soil Science Society of America Journal 67:225-234 (2003)
© 2003 Soil Science Society of America
DIVISION S-5PEDOLOGY
Hillslope Soils and Organic Matter Dynamics within a Native American Agroecosystem on the Colorado Plateau
J. B. Norton*,a,
J. A. Sandorb and
C. S. Whitec
a Dep. of Plants, Soils, and Biometeorology, Utah State Univ., Logan, UT 84322-4820, currently at: Cooperative Extension, Tuolumne County, Sonora, CA 95370
b Dep. of Agronomy, Iowa State Univ., Ames, IA 50011-1010
c Dep. of Biology, Univ. of New Mexico, Albuquerque, NM 87131
* Corresponding author (jnorton{at}co.tuolumne.ca.us)
 |
ABSTRACT
|
|---|
Zuni farmers of western New Mexico demonstrate knowledge about soil and hydrological processes that link upland watersheds to alluvium-derived soils that have crucial cultural, ecological, and hydrological functions. To define how hillslopes contribute to productivity of soils derived from alluvium, we studied soil-vegetation-landform parameters in three headwater drainages on the Zuni Indian Reservation. Analyses along eight summit to toeslope transects show that soil properties follow parabolic and linear trends with changes driven by elevation, lithology, and vegetation in the mesa-canyon pinyonjuniperGambel oak (Pinus edulis Engelm.Juniperus spp.Quercus gambelii Nutt.) woodlands. Total organic C, N, and P concentrations in surface horizons follow negative parabolic trends and are highest on wooded backslopes and lowest on summits and toeslopes. Inorganic N and available P concentrations and total organic to inorganic N and P ratios increase linearly from summit to toeslope. Taken together, soil, landform, and vegetation data suggest: (i) summit positions are relatively stable with immobilizing microbial environments; (ii) inorganic nutrients increase progressively down steep and erodible backslopes as inputs of forest litter are mixed with surface soil; (iii) influx of mixed sediment and organic materials from backslopes maintains concentrations of inorganic nutrients on footslopes and toeslopes. Entrenchment of drainage ways can circumvent these translocation processes. Without the influx of organic materials, footslopes and toeslopes may become nutrient-depleted as immobilization becomes the dominant microbial process. The results underscore the importance of functional connectivity between upland hillslopes and alluvial soils.
Abbreviations: SOM, soil organic matter
 |
INTRODUCTION
|
|---|
SOME OF THE oldest agricultural soils in North America lie within the Colorado Plateau on alluvial fans farmed by the Zuni of western New Mexico (Damp et al., 2002). Zuni farmers recognize the role of upland hillslopes in their strategy to produce corn (Zea mays L.) and other crops in the semiarid environment (Sandor et al., 2002). They actively seek fresh deposits of sediments and organic material for cultivation, and work to enhance fluvial processes that link hillslopes to their fields (Cushing, 1920; Norton et al., 2001; Norton et al., 2002). Sustained productivity of Zuni agricultural soils exemplifies the importance of hillslope processes that support crucial cultural, ecological, and hydrological functions of headwater alluvial fans (Bull, 1997; Peterson et al., 2001). This ancient agriculture, which is intimately linked to natural landscape processes, creates an excellent setting for investigating important relationships between upland hillslopes and nutrient cycling and retention in alluvium-derived soils of headwater ephemeral streams.
Peterson et al. (2001) create a convincing case for conservation that maintains or restores the storage capacity of headwater streams. Their data show that headwater streams play a key role in the physical and biological integrity of downstream waterways because of their great extent on the landscape and their proximity to uplands that are source areas for runoff, sediments, and nutrients. Properly functioning headwater systems retain and recycle a large portion of the nutrients, particularly N, from uplands. This retention and recycling creates the richest soils and highest biological productivity in many arid areas. Bull (1997) points out that many Native American farmers recognized the value of the system and their activities maintained these functions, while modern practices often have the opposite effect.
Although many authors describe hillslope models that link soil and nutrient-cycling processes from summits to alluvial toeslopes (Ruhe and Walker, 1968; Conacher and Dalrymple, 1977; Schimel et al., 1985; Aguilar and Heil, 1988), the unique terrain of the Colorado Plateau creates lithologically segmented hillslopes different from those described in accepted hillslope models. Our objective was to develop a hillslope model based on linked slope morphology, soils, vegetation, and nutrient dynamics on the Zuni Indian Reservation in New Mexico. We focused on hillslopes in headwater drainages whose sandy to loamy alluvial fans have been used by generations of Native American farmers. This type of farming continues to be an essential component of many Native Southwestern cultures. By storing runoff and sediments, mitigating destructive floods, and producing diverse, productive vegetation for livestock and wildlife, alluvium-derived soils of small upland watersheds also serve vital ecological and hydrological functions (Lagasse et al., 1990; Bull, 1997).
 |
MATERIALS AND METHODS
|
|---|
We analyzed soil, landform, and vegetation distribution in three small watersheds above long-term runoff agricultural fields on the Zuni Indian Reservation, New Mexico. The reservation lies on the southeastern part of the Colorado Plateau at 1800 to 2400 m elevation (Fig. 1)
and receives an average of 300 mm of precipitation annually, the majority of which often comes during thunderstorms in July, August, and September. Hillslope morphology on the Zuni Indian Reservation, as in much of the Colorado Plateau, is largely a function of horizontal layers of Cretaceous sandstones and shales. At our study sites, relatively indurate, cliff-forming sandstones act as caprocks over more erodible, slope-forming shales, shaping hillslope profiles much like that defined by Conacher and Dalrymple (1977) (Fig. 2)
. Headwater drainages typically form narrow, steep-walled canyons that drain onto broad alluvial valleys.

View larger version (196K):
[in this window]
[in a new window]
|
Fig. 1. Study area location and physiographic provinces of the Southwest (modified from Cordell, 1997).
|
|

View larger version (26K):
[in this window]
[in a new window]
|
Fig. 2. Typical hillslope profile and distribution of soils and vegetation on eight study hillslopes. All soils are in the Aridic subgroup.
|
|
Zuni is one of 19 pueblo tribes of Arizona and New Mexico well known for persistent, agriculturally based traditions (Ferguson and Hart, 1985). Zuni subsistence depended on corn grown in nonirrigated fields for at least 2000 yr (Kintigh, 1985; Damp et al., 2002). Socio-economic change and assimilation policies beginning with U.S. occupation of Zuni lands in the mid 19th century caused drastic declines in agriculture (Cleveland et al., 1995). We selected three study sites (the Sanchez, Laate, and Weekoty sites) that are representative of headwater drainages above long-term nonirrigated agricultural fields. We worked closely with a group of traditional Zuni farmers who assisted with study site location and provided invaluable insights about Zuni farming techniques.
The Sanchez watershed (68 ha) drains to the north with east- and west-facing hillslopes, the Laate watershed (7 ha) drains to the southwest with northwest- and south-facing hillslopes, and the Weekoty watershed (125 ha) drains to the east with north- and south-facing hillslopes (see Fig. 3)
. Steep-walled canyons cut into sandstone and shale members of the Gallup Sandstone formation characterize each of the watersheds (Anderson et al., 1989). The canyon-floor alluvium in each watershed is bisected by an entrenched ephemeral streambed (arroyo) ranging from a 1-m depth at the Laate site to >6 m at the Sanchez and Weekoty sites. In each case, the arroyo channel ends above a runoff agricultural field, near the canyon mouth.

View larger version (92K):
[in this window]
[in a new window]
|
Fig. 3. Topography, slope positions, transect locations, and representative soil profile locations (see Table 1) within the Sanchez (upper left), Weekoty (upper right), and Laate watershed study sites on the Zuni Indian Reservation.
|
|
To determine the aerial extent of each slope position, we mapped field observations in each watershed on fine-scale (1:5000; 50-cm contour interval) base maps (created by Koogle & Pouls Engineering, Inc., Albuquerque, NM, from 1988 aerial photography). Slope positions were entered as map units and area was calculated using Arc View 3.0 (Fig. 3).
We established two cross-watershed transects at the Sanchez watershed and one each at the Weekoty and Laate watersheds (Fig. 3) for a total of eight divide-to-drainage hillslope transects. To adequately sample the hillslopes, soil, vegetation, and landform samples were collected each 5 m along the transects at the Sanchez site. At the Weekoty and Laate sites, three samples were collected each 15 m along each transect: one on the transect and one 5 m each direction perpendicular from the transect. At each sample point we noted slope position and signs of erosion and deposition, classified and estimated soil cover, and collected 0- to 15-cm depth soil samples. Total vegetation cover by species, litter, bare soil, gravel, cobbles, stones, and boulders were estimated using a square frame constructed from PVC tubing. Plot size varied to reflect size of units sampled (Grieg-Smith, 1983), with 0.25 m2 for herbaceous plants and other soil cover parameters, 1 m2 for shrubs, and 4 m2 for trees. We estimated cover classes in the field and then converted to percentages (cover class midpoints) and log 10 transformed for analysis (Daubenmire, 1968). Soil samples consisted of five to seven 0- to 15-cm depth samples from within each 1-m2 plot mixed and subsampled to provide one sample per plot. We described and sampled soil profiles in each slope position.
Soil samples were air dried in the field and transported to labs for chemical and physical analyses. Particle-size distribution was determined using the sieve and pipette method (Gee and Bauder, 1986), with samples pretreated with 30% hydrogen peroxide for soil organic matter (SOM) digestion and a sodium hexametaphosphate solution for clay dispersion. Soil pH was measured electrometrically using a 1:1 suspension (weight basis) of soil in distilled water using a glass electrode (McLean, 1982). Total C and N concentrations were determined on subsamples ground to pass a 76-µm sieve using a Fissions EA1100 dry combustion CNSHO analyzer (Fissions Inst., Inc., Milan, Italy). Inorganic C concentration was determined with a coulombmeter on a subset of samples and found to be insignificant relative to total C concentration in these surface soils. Total C values were therefore interpreted as organic C. Phosphorus was analyzed in soil samples from the Weekoty and Laate hillslopes and one of the Sanchez hillslopes. Total P concentrations were determined by alkaline oxidation (Dick and Tabatabai, 1977). Available P concentrations were measured by the Olsen extraction method (Olsen and Sommers, 1982). Nitrate
and ammonium
concentrations for the four Sanchez hillslopes were determined in 2M potassium chloride extracts with a Lachat flow-injection procedure (Method 12-107-04-1-B, Lachat Instruments, Milwaukee, WI).
Slope, soil, and landform observations were analyzed as a summit to toeslope continuum and as slope positions as defined by Ruhe and Walker (1968). For continuum analyses, we pooled data from all eight hillslopes based on relative distance from interfluve (actual distance/total transect length). We grouped points by increments of 10 relative distance units and calculated means and standard errors for those groupings. We then ran regression analyses (trendline function in Microsoft Excel) on means calculated for each 10-unit increment along the composite hillslope. We also ran regressions on data from each hillslope, the entire pooled data set, and on unweighted slope position means. The 10-unit increment analyses are presented here for SOM components because they represent transport and transformation along a continuum through the hillslope system.
For analysis of parameters that reflect slope position rather than a hillslope continuum (e.g., soil texture, SOM content, and vegetation), we calculated means and standard errors by slope position for the pooled data set. Means of slope position data from each individual hillslope confirm the pooled means. The means were compared by calculating least significant differences with the GLM procedure (SAS Institute, 1999). To ensure normal distribution and equivalent variance, particle-size distribution data (presented as percentages) were arcsine transformed and soil cover data (also presented as percentages) were log 10 transformed before statistical analyses. Additional information is available in Norton (2000)(p. 180).
 |
RESULTS
|
|---|
Distribution of soils on the eight hillslopes corresponds to distribution of the lithologically controlled slope positions in the steep-walled, mesa/canyon topography (Fig. 2, 3; Table 1). Summit positions in our study areas are generally broad, nearly level mesa tops with shallow Entisols formed in sandstone residuum or sandy eolian materials. Shoulder positions are often dominated by bare sandstone with pockets of shallow, sandy material. Backslope soils typically have A horizons formed in loamy colluvium over siltstone residuum with well developed argillic horizons on north and east facing slopes and none on south and west facing slopes (Table 1). This sequencesandy or loamy surface materials over clayey subsoilsfollows a continuum from shallow A-2Cr horizon sequences on upper backslopes, gradually thickening to deep, well-developed A-Bt-2Cr Alfisols on footslopes. Soils on alluvial toeslopes feature well-developed Bt horizons that indicate long-term stability of the alluvium, overlain by stratified loamy sediments that generally lack soil structure and significant SOM accumulation. Backslopes are both the most extensive and steepest slope position at all three study watersheds (Fig. 3, Table 2). We noted evidence of dynamic erosion and deposition at each of the slope positions as rills, pedestaled plants and pebbles, and accumulations of slope wash materials, including sand, gravel, and organic material, behind downed logs and woody debris.
Vegetation patterns in the three watersheds are marked by distinct patchiness typical of pinyon-juniper woodlands and scrub-shrub grassland communities (Table 3). The patchiness creates high variability in cover plot data (standard deviations are generally higher than means) but consistent patterns emerge from composite analyses of the eight hillslopes. Herbaceous ground cover (grasses and forbs) and tree canopy cover (Table 3) follow opposite trends on the hillslopes. Tree cover is densest on backslopes where there is very little ground cover.
View this table:
[in this window]
[in a new window]
|
Table 3. Estimated percent areal cover averaged across all sites by slope position. Data from all eight hillslopes.
|
|
As with soils, plant species are divided into distinct associations by slope position (Table 3). Backslopes are characterized by pinyon, juniper, and Gambel oak growing singly or in clumps with thick litter layers beneath canopies and bare soil in interspaces. Summit and shoulder positions typically have the largest component of ponderosa pine (Pinus ponderosa P. Lawson & C. Lawson) along with pinyon and juniper, and some Douglas fir [Pseudotsuga menziesii (Mirb.) Franco] in protected areas. Tree stands are more open on summits than on backslopes. The wider interspaces are often covered by grass [typically Stipa spp., blue grama {Bouteloua gracilis (Kunth) Lag. ex Griffiths, nom. illeg. [syn. B. oligostachya (Nutt.) Torr. ex A. Gray]}, bottlebrush squirrel tail [Sitanion hystrix (Nutt.) J.G. Sm. (= Elymus elymoides subsp. elymoides)], and mutton grass [Poa fendleriana (Steud.) Vasey] shrubs {typically wavy-leaf oak (Q. undulata Torr.), mountain mahogany (Cercocarpus montanus Raf.), flowering ash (Fraxinus cuspidata Torr.), and antelope bitterbrush [Purshia tridentata (Pursh) DC.]}, and microbiotic crusts. Toeslopes are covered by big sagebrush (Artemisia tridentata Nutt.), western wheatgrass {Agropyron smithii Rydb. [= Pascopyrum smithii (Rydb.) A. Love]}, and blue grama with components of rabbitbrush [Chrysothamnus nauseosus (Pall.) Britton] and weedy herbaceous vegetation in wash areas.
Particle-size distribution follows distinct trends along the hillslope transects (Table 4) with finest soils derived from siltstone parent materials on backslopes and coarsest soils derived from sandstone parent materials on summits and loamy alluvium on toeslopes. This pattern is opposite of the broadly accepted hillslope models presented by Ruhe and Walker (1968) that describe slopes with constant parent material.
Contents of total organic C, N, and P, as well as C:N ratios, of surface soil closely fit negative parabolic models with highest values on backslopes (Fig. 4)
. Plant-available N and P contents as proportions of total N and P, as well as pH, follow linear models with lowest values on summits and highest on toeslopes. Ammonium-N concentrations on the Sanchez hillslopes increase linearly through the backslope and footslope sections but remain relatively constant along the hillslope continuum as a whole. Nitrate-N content is highly variable and does not follow a distinct trend along the hillslopes. These models fit both the composite hillslope continuum (Fig. 4) and the composite mean values by slope position (Table 4). While we did not differentiate between total- and organic-P fractions, the strong correlation between total P and organic C (P < 0.025) suggests that changes in total-P concentration along the hillslopes likely result from SOM dynamics rather than rock weathering processes.

View larger version (26K):
[in this window]
[in a new window]
|
Fig. 4. Soil organic matter components along summit to toeslope composite hillslope continuum. Backslopes extend from 27 to 55 and toeslopes from 89 to 100 relative distance units along the x-axis. Each point is average of 10 units. Error bars represent standard error; *, significant at the 0.05 level; **, significant at the 0.005 level.
|
|
 |
DISCUSSION
|
|---|
The horizontal stratigraphy of the Colorado Plateau combined with the high-elevation semiarid climate and intense convective summer thunderstorms, results in distinctive mesa-canyon topography with broad, level hilltops and steep, lithologically segmented hillslopes. This topography creates a unique landscape model with discrete landforms that have different but interconnected SOM processes. Fig. 2 is representative of our hillslope transects but the landscape as a whole is more complex, with the sequence of shoulder through footslope positions often repeating several times on mesa slopes (see Fig. 3). In general, the degree of soil development suggests dynamic erosion and deposition of surface horizon material over stable residuum and argillic subsurface horizons at each of the slope positions (Table 1). Similar distribution of clay content, vegetation, and SOM content among the five slope positions, with densest tree canopy and highest SOM content on the soils of backslopes, which have the highest clay contents, suggests that these parameters are controlled by underlying lithology (Table 3, 4, Fig. 4), and, more directly, its impacts on soil water holding capacity. The availability of N and P increases linearly and independently of slope position, however; even as total SOM content decreases (Fig. 4). This suggests that N and P availability are influenced by distance along the hillslope as well as plant community factors, possibly facilitated by mixing and physical disintegration as organic matter and sediments move downslope.
Summit Positions
Surface soils on the broad, gently sloping summit positions in this study are relatively stable with respect to erosion and deposition, as indicated by well-established grasses and microbiotic crusts in canopy interspaces and weathered bedrock residuum parent materials (Table 1). Litter from the mixed ponderosa pine, pinyon, juniper, and oak forests covers >30% of the soil surface (Table 3), but soil nutrients (organic C, total N, and total P) are low compared with those of downslope positions that have even less herbaceous ground cover. This may suggest a lack of mixing from erosion and deposition on the relatively level summits. Low soil C:N ratios suggest in situ decomposition rather than removal and replacement with fresh litter. Low mineral N and P concentrations, as well as low available P:total P and NH+4N:total N ratios suggest an immobilizing microbial environment typical of stable plant communities (Schimel, 1986). This type of stability, where disturbance (e.g., cultivation, fire, or dynamic erosion and deposition) is not a driving component, leads to well-established soil microbial communities that rapidly immobilize inorganic nutrients so SOM turnover rates may be relatively high but concentrations of inorganic nutrients are low (Stark and Hart, 1997).
Backslopes
Backslopes, with steeper gradients, slowly permeable subsoils, and lack of herbaceous understory vegetation, frequently generate runoff (Lagasse et al., 1990). The relatively large areal extent of backslopes suggests that the dynamic erosion and deposition are important factors in the landscape (Fig. 3, Table 2). Surface disturbance in the form of erosion and deposition (evident in numerous rills, pedestaled plants and pebbles, and accumulations of slopewash material) mixes, sorts, and transports soils and forest litter, possibly preventing establishment of an understory plant community. Relatively long slope lengths mean that these effects increase progressively downslope. Soil concentrations of organic C, total N, and total P are highest in the backslope positions, suggesting that slopewash processes mix forest litter with the loamy surface soils. Moving down the hillslope transects, C:N ratios increase markedly at the top of the backslopes where the density of pinyon, juniper, and oaks increases (Table 3). Then C:N ratios gradually decrease through the steep forested slopes. This suggests that organic materials are both moving downslope and decomposing. Rising NH+4N and available P concentrations, along with increasing available P:total P and NH+4N:total N ratios (Fig. 4), suggest a more mineralizing soil environment on backslopes than summit positions, which may also be a function of forest litter decomposing and being carried downslope. The clayey textures of the backslope soils may also contribute to a more mineralizing soil environment by enhancing soil moisture conditions, but the linear increases in mineral N and P concentrations through the foot- and toeslope positions, with progressively lower clay contents, may be better explained as progressive decomposition as material is moved downslope during frequent summer runoff events.
These trends may reflect increasing erosive power with increasing slope length. Rain splash and sheet flow on upper backslopes remove forest litter as it begins to decompose so that SOM is dominated by relatively fresh, high C:N ratio litter in an N-limited, immobilizing soil environment. Plant-available nutrients increase with decreasing elevation as summer runoff events mix partly decomposed material from upslope with fresh litter beneath the pinyon, juniper, and oak canopy, possibly stimulating mineralization. This is suggested by decreasing total concentrations of C, N, and P (overall loss of SOM) but increasing mineral N and P forms through lower backslope and footslope positions (Fig. 4). We observed fresh accumulations of mixed sediment and organic matter after precipitation events.
Footslopes and Alluvial Toeslopes
The soil C:N ratios continue to decrease as mineral N and P increase through footslope and toeslope positions (Fig. 4). Abundant herbaceous ground cover in lower slope positions may be a result of enhanced moisture and available-nutrient conditions. This could lead to a more stable surface soil environment and to microbial immobilization, but the influx of relatively highly mineralized slopewash material (colluvium), combined with occasional erosion and deposition, apparently outpaces and disrupts immobilization processes to maintain long-term productivity even where soils on alluvial valley floors have been farmed for many generations (Norton, 1996; Bull, 1997; Homburg, 2000). Schimel et al. (1985) noted the opposite in the less dynamic landscape (i.e., with denser ground cover, lower erosion rates, and gentler precipitation) of eastern Colorado's shortgrass prairie: decreased mineral N and higher C:N ratios on lower slopes. They attributed this to greater biomass production due to enhanced soil moisture, which caused microbial immobilization to outpace mineralization rates.
Classical catena studies show that quantity and quality of SOM, as well as soil texture, pH, and other properties change systematically down hillslopes (Conacher and Dalrymple, 1977; Jenny, 1980; Gerrard, 1992). Many studies document changes in soil properties with slope position, and define processes associated with sediment movement, sorting, and accumulation as it pertains to soil morphology (Ruhe and Walker, 1968; Kleiss, 1970; Honeycutt et al., 1990). Each of these studies defines increasing SOM concentrations with distance from summit along trends that mirror fine soil fractions. Aguilar and Heil (1988) noted increasing C, N, and P concentrations, as well as narrowing C:N ratios, with decreasing elevation on North Dakota rangeland hillslopes. They suggest accumulation of runoff-transported SOM on lower slopes and infer progressive downslope increases in inorganic N concentrations. Our results generally show opposite trends in both particle-size distribution and concentrations of SOM constituents from these models of hillslope processes. We attribute these differences to the distinct mesa-canyon topography that controls vegetation patterns of the Colorado Plateau landscape. Summit position soils form in sandy parent materials derived from the sandstone caprocks; backslope and footslope soils form in residuum and colluvium of the weathered, slope-forming shales; and toeslopes, which in this hillslope system double as alluvial fans or terraces, form in loamy colluvium and alluvium transported from adjacent hillslopes and deposited over sandy strata from the larger fluvial system (see Fig. 3, Table 1).
Studies of the effects of vegetation on SOM content and composition show that different species can have distinct impacts (Klemmedson, 1991; Klemmedson and Wienhold, 1991). Klemmedson (1991), for instance, found that Gambel oak had significant positive impacts on soil fertility in Arizona ponderosa pine stands; much more than actinomycetes-associated N-fixing shrubs such as mountain mahogany. The predominance of Gambel oak, pinyon, and juniper trees on the steep backslopes of our study areas may partly explain the peak in SOM contents on this slope position.
Similar parabolic trends followed by total organic C, N, and P content, clay content, and tree cover suggest that changes in SOM content are partly associated with particle-size distribution and/or vegetation. Each of these important factors in the soil microbial environment increases from summits to backslopes, possibly from enhanced moisture content with finer soil texture, and then decreases through footslopes and toeslopes as surface soil texture becomes sandier. Carbon:N ratios, however, are high on backslopes and lower on foot- and toeslopes, which is opposite what may be expected if soil texture and associated moisture holding capacity enhanced mineralization. Like SOM and clay content, NH+4N and available P concentrations in surface soils (especially relative to total N and P concentrations) also increase from summits to backslopes but then continue to increase independently of other soil, vegetation, and SOM parameters (Fig. 4). This suggests that mineralization increases as a function of distance through the hillslope system independently of the soil environment (e.g., SOM content and soil texture). We believe that increasing mineralization reflected in decreasing C:N ratios and increasing N and P availability results largely from physical breakdown and translocation of organic materials with frequent summer runoff events that alternate with hot dry conditions.
Most studies of SOM distribution in Southwestern pinyon-juniper woodlands and scrub-shrub grasslands focus on tree and shrub encroachment that increases soil heterogeneity, creating vegetated patches with depleted interspaces (Allen, 1991; Schlesinger et al., 1996). Development of these resource islands is both cause and effect of changing runoff and erosion regimes in a positive feedback relationship (Abrahams et al., 1995). Recent work has focused on effects of vegetation patches on runoff and erosion yields (Wilcox, 1994), redistribution of sediments among pinyon-juniper patches and intercanopy zones (Reid et al., 1999), and effects of shrub-patch-induced accelerated erosion on nutrient losses (Schlesinger et al., 1999). The backslopes of our study sites, with nearly mutually exclusive woody canopy and herbaceous ground cover, appear to be most affected by the resource island effect, which probably contributes to increased runoff, which transports and facilitates decomposition of forest litter.
 |
CONCLUSIONS
|
|---|
The results of this research provide a basis for understanding watershed hillslope characteristics that contribute to sustained productivity of semiarid alluvium- and colluvium-derived soils of Colorado Plateau-type landscapes. Relationships between soil, landform, and vegetation patterns point to the forested backslopes as a driving force behind movement and processing of runoff, sediment, and organic materials. The relatively large areal extent, steepness, slowly permeable subsoils, and combination of bare soil and tree species that produce large amounts of forest litter on backslopes contribute to favorable hydrological properties and fertility of traditional agricultural fields on soils derived from alluvium.
Traditional alluvial fan farming remains an important cultural activity among the Zuni and other Southwestern Native American Tribes, but the importance of hydrological connectivity between upland hillslopes and alluvial landforms extends beyond its value to traditional agriculture. Sediments, organic materials, and runoff from hillslopes are valuable resources for sustaining ecological diversity and productivity in flood plains, riparian areas, alluvial fans, and downstream aquatic systems. When such products of hillslope processes are deposited on functional flood plains and alluvial fans they store sediments, attenuate peak flows, and absorb runoff, which can maintain down stream perennial flows. Loss of functional hydrological connectivity (due to channel incision or channelization) means that products of hillslope erosion are lost to alluvial landforms and become environmental liabilities in a feedback spiral of soil degradation where sediments and constricted flows damage downstream aquatic systems. Channel incision, or arroyo cutting, is an important driver of environmental degradation across the southwestern USA and other arid and semiarid regions of the world (Cooke and Reeves, 1976; Elliot et al., 1999).
This research improves understanding of one component of the connection between upland slopes and alluvium-derived soils on the Zuni Indian Reservation and describes a hillslope model that applies to lithologically segmented landscapes such as the Colorado Plateau. The fluvial system that delivers hillslope materials to alluvial fans, as well as nutrient dynamics and agricultural use and conservation within alluvial fans are other important components warranting further study. Understanding these aspects of key landform processes could contribute to more effective conservation, protection, and restoration of degraded semiarid landscapes.
 |
ACKNOWLEDGMENTS
|
|---|
This work was funded by National Science Foundation grant no. DEB-9528458. We are indebted to the Zuni Sustainable Agriculture Project, the Zuni Conservation Project, and the Zuni Tribe. We thank Jeff Homburg, Todd Carlson, Marnie Criley, and Clara Wheeler for laboratory analysis and Troy Lucio, Lindsay Quam, Vanissa Laahty, and participants in Zuni's Job Training Partnerships Act (JTPA) program for field assistance. We are grateful to Tom DeLuca, Urszula Choromanska, Stephen Siebert, and two anonymous reviewers for technical advice and editorial review. Manuscript preparation was supported by the Utah Agriculture Experimental Station.
Received for publication January 31, 2002.
 |
REFERENCES
|
|---|
- Abrahams, A.D., A.J. Parsons, and J. Wainwright. 1995. Effects of vegetation change on interrill runoff and erosion, Walnut Gulch, southern Arizona. Geomorphology 13:3748.
- Aguilar, R., and R.D. Heil. 1988. Soil organic carbon, nitrogen, and phosphorus quantities in Northern Great Plains rangeland. Soil Sci. Soc. Am. J. 52:10761081.[Abstract/Free Full Text]
- Allen, E.B. 1991. Temporal and spatial organization of desert plant communities. p. 193208. In J. Skujins (ed.) Semiarid lands and deserts: Soil resource and reclamation. Marcel Dekker, New York.
- Anderson, O.J., S.G. Lucas, D.W. Love, and S.M. Cather. 1989. Southeastern Colorado Plateau: New Mexico. Geological Society Fortieth Annual Field Conf., Albuquerque, NM. 28 Sept.1 Oct. 1989. New Mexico Geological Society, Albuquerque, NM.
- Bull, W.B. 1997. Discontinuous ephemeral streams. Geomorphology 19:227276.
- Cleveland, D.A., F.J. Bowannie, D.F. Eriacho, A. Laahty, and E. Perramond. 1995. Zuni farming and United States Government policy: The politics of biological diversity in agriculture. Agric. Human Values 12:218.
- Conacher, A.J., and J.B. Dalrymple. 1977. The nine-unit landsurface model: An approach to pedogeomorphic research. Geoderma 18:1154.
- Cooke, R.U., and R.W. Reeves. 1976. Arroyos and environmental change in the American Southwest. Oxford University Press, New York.
- Cushing, F.H. 1920. Zuni breadstuff. Museum of the American Indian: Heye Foundation, New York.
- Damp, J.E., S.A. Hall, and S. Smith. 2002. Early irrigation on the Colorado Plateau near Zuni Pueblo, New Mexico. Am. Antiq. 67:665676.
- Daubenmire, R. 1968. Plant communities: A textbook of plant synecology. Harper & Row, New York.
- Dick, W.A., and M.A. Tabatabai. 1977. An alkaline oxidation method for determination of total phosphorus in soils. Soil Sci. Soc. Am. J. 41:511514.[Abstract/Free Full Text]
- Elliot, J.G., A.C. Gellis, and S.B. Aby. 1999. Evolution of arroyos: Incised channels of the southwestern United States. p. 153185. In S.E. Darby and A. Simon (ed.) Incised river channels. John Wiley & Sons, New York.
- Ferguson, T.J., and E.R. Hart. 1985. A Zuni atlas. Univ. of Oklahoma Press, Norman, OK.
- Gee, G.W., and J.W. Bauder. 1986. Particle-size analysis. p. 383411. In A. Klute (ed.) Methods of soil analysis. Physical and mineralogical methods. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Gerrard, J. 1992. Soil geomorphology: An integration of pedology and geomorphology. Chapman & Hall, London.
- Grieg-Smith, P. 1983. Quantitative plant ecology. Blackwell Scientific Publ., Oxford.
- Homburg, J.A. 2000. Anthropogenic influences on American Indian agricultural soils of the southwestern United States. Ph.D. diss. Iowa State Univ., Ames, IA. (Publ. No. AAT 9962820).
- Honeycutt, C.W., R.D. Heil, and C.V. Cole. 1990. Climatic and topographic relations on three Great Plains soils: II. Carbon, nitrogen, and phosphorus. Soil Sci. Soc. Am. J. 54:476483.[Abstract/Free Full Text]
- Jenny, H. 1980. The Soil Resource. Springer-Verlag, New York.
- Kintigh, K.W. 1985. Settlement, subsistence, and society in late Zuni prehistory. Anthropological papers of the University of Arizona. Univ. of Arizona Press, Tucson, AZ.
- Kleiss, H.J. 1970. Hillslope sedimentation and soil formation in Northeastern Iowa. Soil Sci. Soc. Amer. Proc. 34:287290.
- Klemmedson, J.O. 1991. Oak influence on nutrient availability in pine forests of central Arizona. Soil Sci. Soc. Am. J. 55:248253.[Abstract/Free Full Text]
- Klemmedson, J.O., and B.J. Wienhold. 1991. Aspect and species influences on nitrogen and phosphorus availability in Arizona Chaparral soils. Soil Sci. Soc. Am. J. 55:17351740.[Abstract/Free Full Text]
- Lagasse, P.F., J.D. Schall, and S.G. Wells. 1990. Runoff and erosional processes in an arid watershed. p. 614632. In Hydraulics/hydrology of arid lands: Proceedings of an international symposium. American Society of Civil Engineers, Fort Collins, CO.
- McLean, E.O. 1982. Soil pH and lime requirements. p. 199222. In A.L. Page et al. (ed.) Methods of soil analysis. Chemical and microbiological properties. Vol. 2. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Norton, J.B. 1996. Soil, geomorphic, and ecological factors in Zuni runoff agriculture. M.S. thesis. Iowa State Univ., Ames, IA.
- Norton, J.B. 2000. Agroecology, hydrology, and conservation of ephemeral streams and alluvial fans, Zuni Pueblo, New Mexico. Publ. No. AAT 9993970. Univ. of Montana, Missoula, MT.
- Norton, J.B., F.J. Bowannie, P. Peynetsa, W. Quandelacy, and S.F. Siebert. 2002. Native American methods for conservation and restoration of semiarid ephemeral streams. J. Soil Water Conserv. 57:250258.
- Norton, J.B., R.R. Pawluk, and J.A. Sandor. 2001. Farmer-Scientist collaboration for research and agricultural development on the Zuni Indian Reservation, New Mexico, USA. p. 107120. In W.A. Payne and D. Keeney (ed.) Sustainability in agricultural systems in transition. ASA Spec. Publ. 64. ASA, CSSA, and SSSA, Madison, WI.
- Olsen, S.R., and L.E. Sommers. 1982. Phosphorus. p. 403430. In A.L. Page et al. (ed.) Methods of soil analysis. Chemical and biological methods. Vol. 2. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Peterson, B.J., W.M. Wollheim, P.J. Mulholland, J.R. Webster, J.L. Meyer, J.L. Tank, E. Marti, W.B. Bowden, H.M. Valett, A.E. Hershey, W.H. McDowell, W.K. Dodds, S.K. Hamilton, S. Gregory, and D.D. Morrall. 2001. Control of nitrogen export from watersheds by headwater streams. Science (Washington, DC) 292:8690.[Abstract/Free Full Text]
- Reid, K.D., B.P. Wilcox, D.D. Breshears, and L. MacDonald. 1999. Runoff and erosion in a pinion-juniper woodland: Influence of vegetation patches. Soil Sci. Soc. Am. J. 63:18691879.[Abstract/Free Full Text]
- Ruhe, R.V., and P.H. Walker. 1968. Hillslope models and soil formation: I. Open systems. 9th Int. Congr. Soil Sci. 4:551560.
- Sandor, J.A., J.B. Norton, R.R. Pawluk, J.A. Homburg, D.A. Muenchrath, C.S. White, S.E. Williams, C.I. Havener, and P.D. Stahl. 2002. Soil knowledge embodied in a Native American runoff agroecosystem. Trans. 17th World Congr. Soil Sci., 2002, Bangkok, Thailand.
- SAS Institute. 1999. SAS user's guide. SAS Inst., Cary, NC.
- Schimel, D.S. 1986. Carbon and nitrogen turnover in adjacent grassland and cropland ecosystems. Biogeochemistry 2:345357.
- Schimel, D., M.A. Stillwell, and R.G. Woodmansee. 1985. Biogeochemistry of C, N, and P in a soil catena of the shortgrass steppe. Ecology 66:276282.
- Schlesinger, W.H., A.D. Abrahams, A.J. Parsons, and J. Wainwright. 1999. Nutrient losses in runoff from grassland and shrubland habitats in southern New Mexico: I. Rainfall simulation experiments. Biogeochemistry 45:2134.
- Schlesinger, W.H., J.A. Raikes, A.E. Hartley, and A.F. Cross. 1996. On the spatial pattern of soil nutrients in desert ecosystems. Ecology 77:364374.[ISI]
- Stark, J.M., and S.C. Hart. 1997. High rates of nitrification and nitrate turnover in undisturbed coniferous forests. Nature (London) 385:6164.
- Wilcox, B.P. 1994. Runoff and erosion in intercanopy zones of pinyon-juniper woodlands. J. Range Manage. 47:285295.
This article has been cited by other articles:

|
 |

|
 |
 
J. B. Norton, J. A. Sandor, C. S. White, and V. Laahty
Organic Matter Transformations through Arroyos and Alluvial Fan Soils within a Native American Agroecosystem
Soil Sci. Soc. Am. J.,
April 5, 2007;
71(3):
829 - 835.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
J. B. Norton, J. A. Sandor, and C. S. White
Runoff and Sediments from Hillslope Soils within a Native American Agroecosystem
Soil Sci. Soc. Am. J.,
March 12, 2007;
71(2):
476 - 483.
[Abstract]
[Full Text]
[PDF]
|
 |
|