Quantifying Basalt Rock Outcrops in NRCS Sol Map Units Using Landsat-5 Data Carol A. Moore, Glenn A. Hoffmann, and Nancy F. Glenn abstract Rock outcrop places many limitations on land use. For this reason, identifying its presence in a soil map unit is as important as identifying major soils. Presently the most accurate way to de- termine the percentage of rock outcrop in any one soil map unit involves transecting, which is time-consuming and not practical for large areas. Remote sensing can provide more accurate data for the surveys and a tremendous time saving. Endmembers, which are spectrally unique materials selected within the scene, were chosen for rock outcrop by analyzing National Aerial Photography Program (NAPP) photography under a stereoscope. Spectral Angle Mapper (SAM) classifications were performed on Landsat data of the study area to classify vegetation and rock outcrop. Global Positioning Satellite (GPS) ground-truth data with 30-meter buffers were used to validate the rock classifications. Although it is likely the amount of rock outcrop is overestimated in the final image, additional methods could ultimately produce a successful methodology to accurately map and quantify rock outcrop. he objective of this work was to develop a remote sensing 250,000 acres. Clark County lies at the foot of the continental divide Tmethodology to quantify basalt rock outcrop for NRCS soil mapping and includes the cold, arid, northeastern Snake River plain. The purposes in the Clark Area Soil Survey. Soil surveys provide users with mountainous area is dissected with many canyons and creeks, all assessments of vegetation, soil, and rock outcrop in given areas. A draining south. The soils are formed dominantly in loess and eobian common member of any soil map unit description in this survey area is sands and have varying degrees of development depending on the age exposed basalt bedrock. The amount of rock outcrop present is highly of the lava flows and where they are found. variable because of the relief of the lava flows and the amount of soil The specific area of study for this remote sensing project is that has been deposited on them. The quantity of forage available for concentrated along the mountain ranges and valley surrounding the domestic animals and wildlife, as well as the placement of routes for cities of Dubois and Spencer, ID (Fig. 1). The elevations for these water pipelines and roads, are highly dependent on the amount and cities are 5148 and 5883 feet, respectively. Pleistocene basalt lava location of rock outcrops. Due to the limitations rock outcrop places covers much of the central and eastern valley floor. Quaternary alluvial on land use, identifying its presence in a map unit is as important as deposits cover the western portion of the valley floor, moving into identifying major soils. Presently the most accurate way to determine Pliocene and Upper Miocene felsic volcanic rocks and rhyolite flows the percentage of rock outcrop in any one soil map unit involves in the mountain range. Vegetation in conducting transects where the soil scientist walks across the landform this area includes antelope bitter brush, and counts the total number of paces. The number of paces that fall on sagebrush, bushy birds beak, buckwheat, rock outcrop divided by the total gives the percentage of rock. While this blue-bunch wheatgrass, and needle-and- is an accurate estimation method for a small area, it is not a practical thread grass (Fig. 2). The average annual air means of determining the amounts of rock outcrop on 250,000 acres of temperature is approximately 43° Fahrenheit land. Remote sensing can provide a valuable tool for these purposes, and the average annual precipitation is providing more accurate data for the surveys as well as a tremendous approximately 13 inches. The average time savings. number of days with at least one inch of The Clark Area Soil Survey covers approximately 770,000 acres snow on the ground is 118. Farming and in Clark and Butte counties. Of this area, lava fields comprise about ranching are the most important industries in the county. Most irrigated lands are C.A. Moore, Idaho State University, Dep. of Geosciences, Idaho Falls, ID situated near the 5200-ft elevation level (email: [email protected]); G.A. Hoffmann, MLRA Soil Survey Project and are a hay base for the livestock. These Leader, USDA-NRCS, Idaho Falls, ID; N.F. Glenn, Research Associate Professor, Idaho State University, Dep. of Geosciences, Boise, ID.Carol lands are irrigated from the Eastern Snake Moore is currently an Idaho State University Junior pursuing a Bachelor River Plain Aquifer. Wheat and other Fig. 1. General location of of Science degree in Earth and Environmental Systems. This research study area within the Clark was designed and conducted as part of a class project and independent small grains are grown on dry farms. Area Soil Survey (Idaho study at ISU. Work was completed during the spring and fall semesters Other nonirrigated areas are mainly of 2006 while the senior author was a student employee with NRCS. Geospatial Data Clearing- Published in Soil Surv. Horiz. 48:59-62 (2007). wildlife habitat. house, 2004). FALL 2007 Endrnember Collection Spectra Veg W2 X. Brush NE2 Brush Outcrop South Cliff Bands West i I - Band Number 4 Data and Methodology Landsat 30-rn resolution data acquired July 6, 2006 were used for this project, along with USDA National Agriculture Imagery Program (NAIP) 1-rn resolution imagery, a digital elevation model (DEM), and black and white National Aerial were created for rock outcrop: one representing basalt outcrop in the Photography Program (NAPP) 1998 photographs at 1:24,000 scale, a Eastern portion of the valley and one representing some cliff bands in geology map, and GPS ground-truth data. All imagery was processed the mountain range surrounding the Western portion of the valley. Using using ESRI ArcGIS 9.1 (Environmental Systems Research Institute, the false color Landsat image and field knowledge as a guide, eight 2005) and Environment for Visualizing Images (ENVI) 4.3 (ITT Visual additional regions of interest were created for various other components Information Solutions, 2006), a software package for analyzing all types such as irrigated agricultural land, brush, grasses, and other vegetation. of remote sensing data. These were selected from arbitrary areas in the image that were The first step in the process was to identify several areas of large believed to contain these components. Due to time constraints, no outcrop to be used as training data. Outcrops of 30 m or larger were ground-truth data were obtained to provide confirmation. These 10 desired due to the 30-rn resolution of the Landsat data. The NAIP image regions of interest were used to define the training data endmembers, or was visually inspected for this purpose, but even when draped over the spectrally unique materials selected within the scene (Fig. 3). Spectral DEM, it was a cumbersome task to differentiate rock outcrop from other profiles of healthy, green vegetation are generally characterized by features. For this reason, the NAPP photography was inspected under higher values in NIR band 4 than in SWIR band 5; however, this can vary a stereoscope. This made the relief much more apparent and seemed greatly depending on variables such as species of plant, environmental to make the task of locating the outcrops a bit easier since many of the conditions during growth, and background reflectance (mixed larger outcrops are located on pressure ridges or large mounds. Several components within the pixel). We can see in Fig. 3 the spectral profiles areas believed to be outcrop were chosen and digitized over the NAIP for Rock, Cliff Bands, and Irrigated Ag clearly fit their respective criteria. image, taking care to ensure the geographic coordinate system datum However, the spectral profiles for the five vegetation endmembers at the matched that of the Landsat datum. top of the plot resemble criteria for rock in SWIR band 5, and the Brush North spectral profile contains high values in both NIR band 4 and SWIR Landsat data often are used for geological mapping. Traditional band 5. Mixed components within the pixels are likely the reason for this, interpretation of the data couples false color composites, principal but we can only speculate due to the lack of ground-truth data. Although component analysis, and ratio images with image interpretation to visually similar, each endmember is spectrally unique and therefore determine the geology of a region (Riley et al., 2006). The traditional provides more data from which to differentiate the rock and cliff band near-infrared false color composites display healthy vegetation as red components during classification. in color, making readily apparent the areas that are not likely to be rock or bare soil. The Landsat data were displayed in false color, and Having now defined the training data, the next step was to perform the shapefile of digitized training polygons was then used as a vector a spectral classification of the image using the Spectral Angle Mapper overlay. Spectral profiles were interactively displayed for various pixels (SAM) classification method. SAM determines the spectral similarity within the digitized training polygons. The Landsat spectral profile of between two spectra by calculating the angle between them. The closer rock is characterized by higher values in the short wave infrared (SWIR) the angle, the more similarity exists. The SAM classification method was band 5 than in near-infrared (NIR) band 4. The spectral profiles for the run using three different methods: (1) using the original stacked Landsat selected pixels seem to fit the basic reflectance characteristics for rock. image, (2) creating a Principal Component Analysis (PCA) 2/1 ratio image and adding it to the Landsat stack to maximize and enhance the A spatial subset of the image was created to reduce the large spectral differences in the image, and (3) running the classification with Landsat image to a smaller size concentrated over the specific area of a spectral subset of the original Landsat stack by removing the thermal study.