United States Department of Agriculture
Natural Resources Conservation Service
Go to Accessibility Information
Skip to Page Content
picture of a mesa





Characterization of the Ecological Integrity of Commercially Grazed Rangelands Using Remote Sensing-based Ecological Indicators

Updated 02/27/2009

Robert A. Washington-Allen1,2 , Neil E. West1, R. Douglas Ramsey3, and Carolyn T. Hunsaker4

1Department of Rangeland Resources, Utah State University, Logan, Utah 84322-5230
2 Environmental Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6407
3 Department of Geography and Earth Resources, Utah State University, Logan, Utah 84322-5240
4 USDA Forest Service, Pacific Southwest Research Laboratory, 2081 E. Sierra Avenue, Fresno, California 93710.

Summary of Research:

We are using 27 years of wet and dry season Landsat satellite imagery from (1972 to 1998); a Geographic Information System (GIS) digital database of site biological, physical, and administrative characteristics and analysis tools; historical and current ranch management records; and a multiple-time by nested multiple-scale experimental design to establish causal links between possible threshold response of ecological indicators diagnostic of land degradation and human management interventions in order to assess the ecological integrity of ecosystems within a semi-arid landscape subject to commercial livestock grazing (Fig. 1).

Range scientists define rangeland degradation as: (1) a change in plant species composition, e.g., an increase in low quality forage species compared with high quality forage; (2) a decrease in plant productivity; (3) a reduction in soil quality; (4) accelerated soil erosion (Pickup 1989, Behnke and Scoones 1993); and a change in landscape composition and pattern that affect ecosystem functions. We view ecological terminology such as landscape and ecosystem as spatio-temporal scale-independent criteria (sensu Allen and Hoeskstra 1992). Ecological integrity is defined as an equilibrial or non-equilibrial ecosystem or landscape which possesses a full set of natural components and processes in "good working order" (West et al. 1994). Composition and pattern refer to the spatial organization of ecosystem characteristics when viewed as patches within a mosaic, e.g., a shrub vegetation patch within a bare soil and grassland mosaic (Noy-Meir 1973, Ludwig and Tongway 1995). Measurement indicators derived from these five assessment endpoints (an explicit statement of the actual environmental value that is to be protected, Suter 1993) will provide a diagnostic basis for an estimate of the degree of land degradation. We are using indices of vegetation cover and soil integrity, the Soil Adjusted Vegetation Index (SAVI, Huete 1988) and the Soil Stability Index (SSI, Pickup and Nelson 1984), and metrics of landscape structure (Turner 1989) on categorical maps of plant lifeforms and bare soil. Figure 2 is an overview of how the analyses of the measurement idicators are related to their derivation from the same landscape and its changes in space and time (Figure 2A). The analyses concentrate on possible thresholds associated with state changes of ecological sites (Figure 2B) and landscape pattern (Figure 2C), changes at concentration areas, e.g., waterpoints, and the use of change vector analyses to measure ecological resilience (Figure 2D). Discrete or rapid change in an indicator is indicative of a threshold. A transition threshold is defined as the boundary in space and time between two seral states, e.g., the observed changes of grassland to predominantly woody vegetation (Archer 1989, Figure 3). Land managers have a different perception of thresholds that is primarily informed by ecological and economic considerations. From their perspective, a transition threshold is the point where the initial shift across the boundary is not reversible on a practical time scale without substantial intervention of cultural inputs, e.g., heavy machinery, or fire (Archer 1989, Friedel 1991, Figure 3 (the right-hand Y-axis)).

The study site selected for this research is the 88,000 ha Deseret Land and Livestock Ranch (McMurrin 1989) located in the northeast corner of the Utah panhandle. The first stage of our project has been GIS and remote sensing data acquisition and development (Table 1). Extension activities with Deseret's land managers has yielded eighteen years of livestock (1980 to 1997) and sixteen years of wild ungulate (1982 to 1997) management data (Figures 4 and 5). Twenty-six years of wet and dry season imagery (49 scenes) and two years of annual imagery (1985 and 1986, 9 scenes/yr) are being image-to-image rectified. Image normalization of the data set is occurring via transformation to exo-atmospheric reflectance values and use of an empirical radiometric normalization technique to relatively atmospherically-correct imagery to a reference image. Landsat MSS data collection ceased after 1992, thus wet and dry season Landsat Thematic Mapper (TM) imagery from 1986 to 1997 (Table 1.) were acquired and multiple regression techniques are being used to allow TM sensor response to simulate MSS sensors. The image processing procedures used in this study are labor intensive, but we are automating this process in collaboration with scientists from a DoD Strategic Environmental Research and Development Program (SERDP) conservation project that is concerned with developing remote sensing-based environmental monitoring protocols for military bases.

Results

Seasonal images from 1972 to 1992 were converted to the soil-adjusted vegetation index (SAVI) images and analyzed (Fig. 6). SAVIseasonal state space mean/variance plots for Deseret at the landscape scale reveal that vegetation cover dynamics clearly discriminate into two domains of high mean/low variance in the 1970s and low mean/high-low variance in the 1980's and 90s. The two domains corresponded to wet and drought periods indicated by the Palmer Drought Severity Index (PDSI). The mean/variance plots temporal trajectories were aperiodic with somewhat periodic orbit pattern which identify the two domains as "strange attractors". This pattern is reminiscent of Lorenz's Butterfly and thus indicative of a chaotic system. In preliminary studies, seasonal images from 1972 to 1992 were converted to SAVI images and there trend analyzed Three hypotheses were tested at the landscape and administrative scales: which were 1) at the landscape scale, the El Niño wet year of 1983-84 would be manifest when compared to mean conditions; 2) The period from 1972 to 1983 would have lower vegetational cover than the period 1984 to 1992 because of better ranch management; and at the administrative scale, 3) vegetation cover within Bureau of Land Management (BLM) islands would not differ from that of the privately-owned paddocks they were within because Deseret Land and Livestock does not separately manage on the basis of land ownership. We found that 1) the El Niño wet year of 1983-84 allowed most of the landscape to recover from an apparent drought in 1982, except for riparian areas which had lower than mean vegetational cover 2) Contrary to our hypothesis, the period from 1972 to 1983 had higher vegetation cover than the period 1984 to 1992; and 3) as expected, vegetation cover between BLM islands and the paddocks which contained them did not differ with correlation coefficients ranging from 0.86 to 0.99 (Figure 7). The overall SAVI trend suggests a decline in vegetational cover and examination of the trend at the scale of individual vegetation types revealed a similar response compared to the landscape scale. This response suggests a climatic constraint. However, comparison at the landscape-scale of the lagged water year PDSI and the SAVI had a low correlation of 0.17. There are anomalous years: 1978 (high SAVI and dry year), 1981 (high SAVI and dry year) , and 1982 (low SAVI and wet year) which are not consistent with the wet and dry precipitation periods. However, the relationship was fairly consistent for the dry period 1987 to 1992. A possible explanation for the behavior in 1982 may be related to the livestock and wildlife grazing record and other land management practices which have occurred, e.g., prescribed burning. Later analyses will consider these factors.

Future Activities:

We will complete the database, including production of soil attribute layers and digital elevation models, acquisition and scanning of aerial photos for accuracy assessment, and then we will proceed to concentrate on analyses. We will analyze the temporal and spatial behaviour of remote sensing-based ecological indicators derived from various definitions of rangeland degradation. The assessment of these indicators will occur at multiple scales including: landscape, watershed, administrative (i.e., public versus private land), individual paddock, ecological site, and piosphere (waterpoints). Continued extension activities will occur with Deseret's natural resource managers, including transfer of Deseret's GIS database to CD-ROM for use by Deseret.

Literature Cited

Allen, T. F. H. and T. W. Hoekstra. 1992. Toward a Unified Ecology, Columbia University Press, New York, 384 p.

Archer, S. 1989. Have southern Texas savannas been converted to woodlands in recent history? American Naturalist 134:545-561.

Behnke, R.H. and I. Scoones. 1993. Rethinking Range Ecology: Implications for rangeland management in Africa. pp. 1-30, In Behnke, R.H., I. Scoones and, C. Kerven (Eds.) Range Ecology at Disequilibrium. Overseas Development Institute, International Institute for Environment and Development, and Commonwealth Secretariat, London.

Friedel, M.H. 1991. Range condition assessment and the concept of thresholds: A Viewpoint. Journal of Range Management 44:422-426.

Huete, A.R. 1988. A soil-adjusted vegetation index (SAVI). Remote Sensing of Environment 25:295-309.

Ludwig, J.A. and D.J. Tongway. 1995. Spatial organization of landscapes and its function in semi-arid woodlands, Australia. Landscape Ecology 10:51-63.

McMurrin, J. 1989. The Deseret Live Stock Company : the first fifty years, 1890-1940. MS Thesis, University of Utah, Salt Lake City, Utah.

Noy-Meir, I. 1973. Desert ecosystems: environment and producers. Annual Review of Ecology and Systematics 4:25-51.

Pickup, G. 1989. New land degradation survey techniques for arid Australia: problems and prospects. Australian Rangeland Journal 11:74-82.

Pickup, G. and D.J. Nelson. 1984. Use of Landsat radiance parameters to distinguish soil erosion, stability and deposition in arid central Australia. Remote Sensing of Environment 16:195-209.

Schlesinger, W. H., J.F. Reynolds, G. L. Cunningham, L.F. Huenneke, W. M. Jarrell, R.A. Virginia, and W.G. Whitford. 1990. Biological feedbacks in Global Desertification. Science 247:1043-1048

Suter II, G.W. 1993. Ecological Risk Assessment. Lewis Publishers, Boca Raton. 538 p.

TGUCT (Task Group on Unity in Concepts and Terminology). 1995. New Concepts for assessment of rangeland condition. Journal of Range Management 48:271-282.

Turner, M.G. 1989. Landscape Ecology: The effect of pattern on process. Annual Review of Ecology and Systematics 20:171-197.

West, N.E., K. McDaniel, E.L. Smith, P.T. Tueller, and S. Leonard. 1994. Monitoring and interpreting ecological integrity on arid and semi-arid lands of the western United States. Western Regional Research Coordinating Committee-40. New Mexico Range Improvement Task Force, Report No.37, Las Cruces, New Mexico.


Table 1. Summary of GIS data available to characterize the natural resources of Deseret Land and Livestock Ranch.
Theme Attributes Dates Spatial Coverage Resolution/Spatial Error Estimate Data Type Reference
Vegetation Cover 31 Vegetation Cover Types and 5 Land Use Types based on UNESCO Hiearchichal classification system June -July 1984 - 1993 Utah Entire State 100 ha minimal mapping unit (MMU)

Riparian and Wetland have 40 ha MMU

Raster Map derived from mosaiced Landsat Thematic Mapper Scences (acquired from Utah Gap Analysis Project), UTM Projection, NAD 27 Edwards, Jr., T.C., C.G. Homer, S.D. Bassett, A. Falconer, R.D. Ramsey, and D.W. Wight. 1995. Utah Gap Analysis: An Environmental Information System. Final Project Report 95-1, Utah Cooperative Fish and Wildlife Research Unit, Utah State University, Logan, Utah 84322-5210.
Land Thematic Mapper Imagery 7-Band sensed image (Infrared x 3, thermal red, green, and blue), 8 bit Raster July 1986 to September 1998 180 km x 180 km scene

Path 38 Row 

30 m resampled to 57 m

 

Raster Landsat Thematic Mapper (TM) 7 x 8-bit bands raster image, full scenes acquired from EROS Data Center, UTM Projection, NAD 27 Receiving Station: USGS EROS Data Center, South Dakota, USA

Earth Observation Satellite Company, 1985, User's guide for Landsat thematic mapper computer-compatible tapes: Lanham, Maryland, Earth Observation Satellite Company [variously paged].

Earth Observation Satellite Company, 1994, Landsat system status report--September 1994: Lanham, Maryland, Earth Observation Satellite Company, p. 1-11.

Digital Elevation Model (DEM) topographic map hypsography overlay 1992 7.5 Minute Quadrangle  30 m 

A vertical RMSE of 7 m -

15 m

Level-1 DEM scanned from NHAP/NAPP aerial photography or from digitized contour lines from cartographic maps, Single band, 16-bit Digital Elevation Model (DEM), UTM Projection, NAD 27 U.S. Geological Survey, 1993, Digital elevation models -- data users

guide 5: Reston, Virginia, U.S. Geological Survey, 48 p.

U.S. Geological Survey, 1993, US GeoData digital elevation model,

factsheet: [Reston, Virginia], U.S. Geological Survey, 2 p.

Landsat Multi spectral Imagery 4-Band sensed image (Infrared x 2, red and green), 6 and 7 bit Raster September 1972 to September 1992 180 km x 180 km scene 25m raster Landsat Multi spectral Scanner 4 x 6 and 7-bit bands raster image, full scenes acquired from EROS Data Center. UTM Projection, NAD 27 <

NALC: a Landsat MSS triplicate that is comprised of three dates of co-registered imagery and a DEM . Imagery were acquired in 1973, 1986, and 1991 plus or minus one year. the DEM is 16-bit (INTEGER*2) data.

Receiving Station: USGS EROS Data Center, SD, USA

Thomas, V.L., 1977, Generation and physical characteristics of the Landsat-1, -2 and -3 MSS computer compatible tapes: Greenbelt, Maryland, National Aeronautics and Space Administration, Goddard Space Flight Center, Technical Memorandum 78018 [variously paged]

U.S. Geological Survey and National Oceanic and Atmospheric Administration, 1984, Landsat 4 data users handbook: [Washington, D.C.], U.S. Geological Survey and National Oceanic and Atmospheric Administration [variously paged].

Lunetta, R.S., and Sturdevant, J.A., 1993, The North American Landscape Characterization Landsat Pathfinder Project, in Pettinger, L.R., ed., Pecora 12 Symposium, Land Information from Space-Based Systems, Proceedings: American Society of Photogrammetry and Remote Sensing, Bethesda, Maryland, pp. 363-371.

Geology formation 1980 Utah, Entire State 25m in raster vector/polygon, UTM Projection, NAD 27, Hand digitized from 1:500,000 map Hintze, L.F. (Compiler)1980. Geological Map of Utah. 1/500,000 colored with cross sections and stratigraphic columns. Published by compiler, Provo, Utah.

Hand digitized by Utah State University R/S GIS Research Laboratory in 1992.

Soils soil series and related physical attributes from MUIR NRCS database. 1982 Rich County, Utah 1 ha vector/polygon, UTM Projection, NAD 27 Digitized by Robert A. Washington-Allen from

Rich County Soil Survey 1982

Soil Conservation Service. 1982. Soil Survey of Rich County, Utah. USDA. U.S. Government Printing Office, Washington, D.C. \

Joined to National Map Unit Interpretation Record (MUIR) database from NRCS for Rich County, Utah

USDA-NRCS.1994. National Map Unit Interpretation Record (MUIR) Database Publication Information. USDA-NRCS, Fort Worth, Texas.

Soils Soil associations and related physical attributes from STATSGO database 1982 United States of America 625 ha vector/polygon, UTM Projection, NAD 27 USDA-NRCS State Soil Geographic Database (STATSGO) U.S. Department of Agriculture, 1994, State soil geographic (STATSGO) data base  -  data use information, miscellaneous publication number 1492(rev. ed.): Fort Worth, Texas, Natural Resources Conservation Service [variously paged].
Roads primary

secondary

tertiary

trails

jeep roads

1,994 Utah, Entire State within ± 80 m vector/line, UTM Projection, NAD 27 SGID Utah State Geographic Information Database by The Automated Geographic Reference Center (AGRC) a division of the State of Utah, Department of Administrative Services, Information Technology Services, Salt Lake City, Utah 84114, (801) 538-3165
Grazing Paddocks Paddock name joined to digital eighteen year livestock and sixteen year wildlife grazing management data. 1,998 88,000 ha 90% of data within ± 40 ft vector/line and polygon, UTM Projection, NAD 27 Digitized by Robert A. Washington-Allen from

USGS 7.5 Min Topographic Maps 1976 and 1996

Hydrology streams, ponds, and tributaries 1,994 Utah Entire State within ± 80 m vector/line and polygon, UTM Projection, NAD 27 SGID Utah State Geographic Information Database by The Automated Geographic Reference Center (AGRC) a division of the State of Utah, Department of Administrative Services, Information Technology Services, Salt Lake City, Utah 84114, (801) 538-3165
Administrative Boundary Desert ownership boundary 1994 88,000 ha 90% of data within ± 40 ft vector/line and polygon, UTM Projection, NAD 27 Digitized by Robert A. Washington-Allen from

USGS 7.5 Min Topographic Maps 1976 and 1996

Shrub Cover Density Shrub density 1988 Rich County 30 m raster raster, UTM Projection, NAD 27 Homer, C.G., T.C. Edwards, Jr, R.D. Ramsey, and K.P. Price. 1993. Use of Remote Sensing methods in modeling sage grouse winter habitat. Journal of Wildlife Management 57:78-84.
Grass and Shrub Density Cover Map Vegetation Density 1988 Rich County 30 m raster raster, UTM Projection, NAD 27 Hunnicutt, M. 1992. Use of Landsat Imagery and Geographical Information Systems in the assessment of Rangeland Cover and Wildlife Habitat. M.S. Thesis, Utah State University, Logan. USA.
Waterpoints location 1997 88,000 ha  90% of data within ± 40 ft vector/point, UTM Projection, NAD 27 Digitized by Robert A. Washington-Allen from

USGS 7.5 Min Topographic Maps 1976 and 1996

Utah_state _bnd State Boundary 1994 Utah Entire State 90% of data within ± 40 ft vector/line and polygon, UTM Projection, NAD 27 SGID Utah State Geographic Information Database by The Automated Geographic Reference Center (AGRC) a division of the State of Utah, Department of Administrative Services, Information Technology Services, Salt Lake City, Utah 84114, (801) 538-3165
Rich_Bnd County Boundary 1994 Entire County 90% of data within ± 40 ft vector/line and polygon, UTM Projection, NAD 27 SGID Utah State Geographic Information Database by The Automated Geographic Reference Center (AGRC) a division of the State of Utah, Department of Administrative Services, Information Technology Services, Salt Lake City, Utah 84114, (801) 538-3165
aerial photographs B&W and Colour IR 1972 to 1997 indexed on 1:100,000 - 250,000 map reference sheets 1 to 2 meter object resolution analog B&W

analog CIR aerial photo

analog Orthophotoquadrangle, UTM Projection, NAD 27

U.S. Geological Survey, 1992, The National Aerial Photography Program (NAPP), fact sheet: Reston, Virginia, U.S. Geological Survey, 1 p. 

U.S. Geological Survey, 1992, NHAP and NAPP photographic enlargements, fact sheet: Reston, Virginia, U.S. Geological Survey, 1 p.

Figure 1. An overview of the proposed GIS/Remote Sensing-based framework for assessing the degree of land degradation on a landscape managed for commercial livestock grazing.

Figure 2. Conceptual basis for the development of a multi-temporal Remote Sensing/GIS-based assessment of rangeland integrity.

Figure 3. State Changes and Transition Thresholds

Figure 4. The Stocking Rate (Animal Unit (AU) per hectare for a year) on Deseret Land & Livestock Ranch from 1980 to 1997.

Figure 5. Population of elk (A), pronghorn (B), and deer (C) on Deseret Land & Livestock Ranch from 1982 to 1997. The acronym AD means animal days, but in this case is equal to animal numbers.

Figure 6. Mean /Variance plots of the Wet and Dry Season Soil Adjusted Vegetation Index (SAVI). SAVI indicates the amount of vegetation cover on the ranch where high variance and low mean years are indicative of drought and increased susceptibility to accelerated soil erosion.

Figure 7. Comparison from 1975 to 1992 of the soil adjusted vegetation index (SAVI), an index of vegetation cover, between Bureau of Land Management land islands and the private land which contains them.



< Back to Agenda