Remote Sensing Applications and Decision Support

Development of an integrated multiplatform approach for assessing brush management conservation efforts in semiarid rangelands

[+] Author Affiliations
Chandra D. Holifield Collins, Guillermo Ponce-Campos

USDA ARS Southwest Watershed Research Center, 2000 E. Allen Road, Tucson, Arizona 85719, United States

Mark A. Kautz, Sapana Lohani, John Hottenstein

University of Arizona, 1230 E. Speedway Boulevard, P.O. Box 210104, Tucson, Arizona 85721, United States

Ronald Tiller

The Nature Conservancy, 1510 E. Fort Lowell Road, Tucson, Arizona 85719, United States

Loretta J. Metz

USDA NRCS Resource Assessment Division (CEAP-Grazing Lands), Grassland, Soil and Water Research Laboratory, 808 E. Blackland Road, Temple, Texas 76502, United States

J. Appl. Remote Sens. 9(1), 096057 (May 07, 2015). doi:10.1117/1.JRS.9.096057
History: Received October 6, 2014; Accepted April 7, 2015
Text Size: A A A

Open Access Open Access

Abstract.  Millions of dollars have been spent on brush management, or removal of unwanted woody vegetation, as a conservation practice to control the presence of woody species. Land managers need an inexpensive means of monitoring the effects of brush management conservation methods for decreasing degradation in rangeland systems. In this study, free, publically available, high-resolution (1 m) imagery from the National Agricultural Imagery Program (NAIP) and moderate-resolution (30 m) Landsat-5 Thematic Mapper (TM) imagery were combined to produce a large-scale technique for mapping woody cover. High-resolution imagery-based estimates of woody cover were found to be reasonable (RMSE=3.8%, MAE=2.9%) surrogates for ground-based woody cover. An equation for TM-derived woody cover was developed. TM scenes of woody cover (TMWC) were produced and validated using NAIP and ground-based data. Results showed that the developed relation produced viable (RMSE=8.5%, MAE=6.4%) maps of woody cover that could be used to successfully track the occurrence of brush removal, as well as monitor the presence or lack of subsequent reemergence. This work provides land managers with an operational means of determining where to allocate resources to implement brush management, as well as a cost-effective method of monitoring the effects of their efforts.

Figures in this Article

There has been an increase in trees and shrubs on rangeland systems in the past 150 years.13 This increase in woody species comes at the expense of perennial grasses whose losses can cause changes in primary production, nutrient cycling, and the accumulation of soil organic matter.4 In order to decrease the loss of inherent ecosystem function, vegetation mortality, and increased runoff associated with rangeland system degradation, it is necessary to utilize effective conservation practices. One such practice is brush management, which typically employs prescribed burning, mechanical methods (e.g., chaining, roller chopping, root plowing, shredding, and bulldozing), chemical methods (i.e., herbicide application), or a combination of these methods to control unwanted woody vegetation.5 Between 1997 and 2003, the United States Department of Agriculture’s (USDA) Natural Resources Conservation Service (NRCS) spent nearly $34 million dollars on conservation practices, over $19 million being for brush management alone, on 188 million ha of central and western rangelands and grazed forests.6 With such a large investment of resources being allocated toward the implementation of brush management, it is essential to have a mechanism for evaluating its effectiveness. The cost combination of labor, equipment, vehicles, travel, and supplies, make traditional ground-based methods of data collection on large scales prohibitively expensive. A pilot study to sample 400 sites (800 points) on 3.1 million ha of rangeland and pinyon-juniper woodlands, reported costs of sampling 448 points by 12 people working 10–12 h per day (data collection and travel) from early June to mid-October at approximately $400,000, with field data collection cost before analysis, averaging $893 per sample point.7 Evaluating nearly 200 million ha solely using on-the-ground resources is simply not feasible. A new large-area method is required.

Due to its ability to provide information faster and more cost-effectively than ground-based data collection, remote sensing has been used to map, monitor, and assess vegetation over large areas for over 30 years.817 Within the last 10 years, new government programs, policy changes, and technological advances have made remotely sensed data more accessible than ever before. The USDA Farm Service Agency’s (FSA) National Agriculture Imagery Program (NAIP) has been collecting and providing free, publically available, high-resolution aerial imagery for the continental United States since 2003. Further, the global, multidecadal Landsat satellite data record (1980s to present), together with Landsat’s free data policy, make county, state, regional, and national scale vegetation assessment possible. There was a time when cost made the use of remotely sensed products for long-term monitoring purposes unrealistic. However, with the availability of these no-cost, high spatial, and high temporal resolution datasets, remotely sensed assessment of the effectiveness of conservation practices that strongly influence vegetation cover (e.g., brush management) is now attainable.

Remote sensing has been used to assess woody vegetation in arid and semiarid rangeland systems around the world.8,1824 Assessments are commonly performed using spectral mixture modeling18,19,2529 and land cover maps2022,24,30 produced through supervised classification techniques. These methods often employ the use of high-resolution (<15m) satellite imagery, aerial photography, or hyperspectral data. However, the nature of these methods does not easily lend itself to operational use. High-resolution imagery is often expensive to obtain, not readily available on an extended temporal basis, and the requirements that attend its use (e.g., software and skilled labor for data preprocessing, development of image processing algorithms, and data management and storage) present additional challenges. Likewise, resource investments (people, time, money, and equipment) for collection of the data required for ground-truth land cover maps are often prohibitive for repetitive map production. To be of practical use by land managers, the quantification of large-scale woody vegetation cover must be possible using a fairly straightforward and easy to apply methodology that utilizes inexpensive (free) imagery and requires minimal image processing (due to time constraints of overburdened personnel). Thus, the objective of this study was to develop an operationally oriented method for assessing woody vegetation cover using publically available, no-cost remotely sensed data sources on semiarid rangelands in southeastern Arizona.

The goal of this work was to develop an operational remote sensing protocol for quantifying woody vegetation cover using no-cost imagery. This protocol is intended for use by land managers to aid in implementing, monitoring, and assessing the effects of brush management conservation practice applications. The woody species Prosopis L. (mesquite) was the focus of this particular study. Sections 2.12.3 include descriptions of the study areas, ground-based datasets, imagery used, its processing, the integration process used for subpixel discernment of woody vegetation, and the generation of the relation to quantify woody cover.

Study Areas

The study was conducted within major land resource area (MLRA)-41 Southeastern Arizona Basin and Range31 and on the Empire Ranch in southeastern Arizona (Fig. 1). The Arizona portion of MLRA-41 covers approximately 36,281km2. Elevation ranges from 975 to 1525 m and the average annual precipitation ranges from 304 to 406 mm with 60% of the rainfall occurring between July and September. The Empire Ranch, which is operated by the Bureau of Land Management (BLM), is located within the Las Cienegas National Conservation Area in the Sonoita Valley of southern Arizona. It is a high profile multiple use area which includes BLM, Pima County, State Trust, and private lands. The ranch is managed cooperatively with significant citizen input, is actively used for cattle ranching, and has experienced a number of conservation management practices (prescribed grazing, prescribed burning, and mechanical grubbing for brush removal). The dominant vegetation is native grass, including Bouteloua gracilis (Willd. ex Kunth) Lag. ex Griffiths (blue grama), Bouteloua eriopoda (Torr.) Torr. (black grama), Bouteloua curtipendula (Michaux) Torr. (sideoats grama), Bouteloua hirsuta (Lagasca) (hairy grama), and Aristida spp. Woody species are mainly comprised of Prosopis L. (mesquite) and Populus L. (cottonwood).

Graphic Jump Location
Fig. 1
F1 :

Major land resource area 41 overlaid with a Landsat-5 TM tile for path 35, row 38, and the Empire Ranch study area with National Agricultural Imagery Program (NAIP), National Resources Inventory (NRI), The Nature Conservancy (TNC), brush removal, and Empire Ranch sampling block (ERBLOCK) study sites. NAIP and NRI represent high-resolution imagery subset locations. TNC and ERBLOCK represent ground-based sampling locations.

Mechanical grubbing for brush removal was performed at several locations on the Empire Ranch. Brush removal treatments occurred from October 18, 2010 to January 26, 2011 on three sites covering approximately 180 ha (Fig. 1). These sites were used to conduct a pre- and post-treatment assessment of woody cover.

Data Collection and Processing

Five years of pre-existing ground-based data collected on the Empire Ranch and at National Resources Inventory (NRI) sites within the MLRA were used for the study (Table 1). The NRI and Empire Ranch sampling block (ERBLOCK) data were collected by NRCS and USDA Agricultural Research Service Southwest Watershed Research Center personnel, respectively. The Nature Conservancy (TNC) data were collected by TNC employees and volunteers as part of an annual monitoring program for the Empire Ranch. The line-point intercept32 method was used to measure foliar canopy cover by species at all locations. Only woody species cover data were used for analysis. Five years of remotely sensed data were also obtained for the study (Table 2). The datasets NRI, NAIP, and Landsat-5 thematic mapper (TM) imagery were used for training, validation, and equation development. Descriptions of image selection criteria and processing procedures are provided in subsections 2.2.12.2.3.

Table Grahic Jump Location
Table 1Ground-based dataset description. NRI=National Resources Inventory, TNC=The Nature Conservancy, ERBLOCK=Empire Ranch sampling block.
Table Grahic Jump Location
Table 2Imagery dataset description. NRI=National Resources Inventory, NAIP=National Agricultural Imagery Program, TM=Landsat-5Thematic Mapper.

It should be noted that dataset acquisition dates ranged from 2002 to 2013. However, barring the occurrence of major disturbance (e.g., mechanical removal, fire), mesquite cover changes very slowly (<1% per year) in the Southwest.3335 Given that no disturbance events occurred within the areas examined for this study, the assumption was made that comparisons made of datasets acquired within a 3-year window of one another were valid.

National Resources Inventory imagery

The NRCS rangeland NRI image scenes are ultra-high-resolution (0.305-m), three-band (red, green, and blue) aerial photographs that have been collected as part of the NRI36,37 data collection efforts since 2000. Five scenes were used for the study (Table 2). Each scene covered approximately 1609m2 and was subset to correspond with ground data collection sites. All subsets were classified as woody or nonwoody through supervised object-oriented classification performed with Overwatch Systems LTD’s Feature Analyst®. Object-oriented analysis utilizes both the spectral characteristics of pixels and their spatial arrangement.24

National Agriculture Imagery Program imagery

Three-band (red, green, and blue) mosaicked NAIP images for Pima (NAIPP) and Cochise (NAIPC) counties collected on August 07, 2013, and July 23, 2013, respectively, were obtained from the USDA NRCS Geospatial Data Gateway.38 The NAIP scenes were 1-m ground sample distance ortho-imagery rectified within ±6m to true ground at a 95% confidence level and formatted to a UTM NAD83 coordinate system. Subsets of the images (NAIPP1, NAIPP2, NAIPC1, and NAIPC2) were selected to represent vegetation found in the mesquite grassland portions of the TM tiles (Fig. 2). NAIPP1 and NAIPC1 are examples of grassland areas with low to moderate levels of mesquite. NAIPC1 was more sparsely vegetated with greater soil heterogeneity. NAIPP2 and NAIPC2 are examples of mesquite dominated areas, with NAIPP2 being largely composed of mesquite and the large bunchgrass Sporobolus airoides (Torr.) Torr. (alkali sacaton). All subsets were classified as woody or nonwoody through supervised object-oriented classification performed with Overwatch Systems LTD’s Feature Analyst®. Accuracy assessments of each scene were performed using error matrices with 100 randomly distributed sample points per class.39

Graphic Jump Location
Fig. 2
F2 :

Examples of NAIP subsets from Pima and Cochise counties collected on August 07, 2013 and July 23, 2013, respectively.

Thematic Mapper imagery

Moderate-resolution (30 m) TM image tiles of path 35, row 38, selected for the study were downloaded from the USGS Earth Explorer website40 and had been atmospherically corrected to surface reflectance using the Landsat Ecosystem Disturbance Adaptive Processing System.41,42 Scenes came processed to standard terrain correction (level 1T), which provides systematic radiometric and geometric accuracy by integrating ground control points while utilizing a digital elevation model for topographic accuracy. All scenes were converted to UTM North American Datum of 1983 (NAD83) zone 12.

The landscape of MLRA-41 is a mix of heterogeneous vegetative cover and soil types. Therefore, characterizing its vegetation canopies will require the use of a spectral vegetation index that can handle such variability. One of the greatest concerns for assessing the sites used in this study was the influence of soil background conditions, which can have considerable influence on calculated vegetation indices for partial canopy spectra.43,44 The normalized difference vegetation index (NDVI)45Display Formula

NDVI=ρNIRρredρNIR+ρred,(1)
where ρNIR and ρred are near-infrared (NIR) and red reflectance, is one of the most commonly used indices for assessing vegetation with remotely sensed data. However, the NDVI is sensitive to optical properties of soil background, particularly when vegetation is sparse, making it problematic for application to data spanning different soil types.43,45,46 The soil-adjusted vegetation index (SAVI) Display Formula
SAVI=ρNIRρredρNIR+ρred+L(1+L),(2)
is capable of minimizing the effects of variable soil backgrounds through the use of an L adjustment factor.43 However, to obtain the most precise measures of SAVI, knowledge of vegetative densities for the area of interest is needed to select the most appropriate value for L (e.g., L=1 for very low vegetation densities, L=0.5 for intermediate densities, and L=0.25 for higher densities).43,44 This can be an issue when looking at large areas containing varying densities of vegetation. Unlike SAVI, the modified soil-adjusted vegetation index (MSAVI) Display Formula
MSAVI=2ρNIR+1(2ρNIR+1)28(ρNIRρred)2,(3)
requires no prior knowledge about vegetation densities in order to select the best L value. Instead optimal L values are determined based on the present reflectance.47 Due to the heterogeneity of the study area, MSAVI was selected for use in this study.

Given that a vegetative index will only provide a general measure of greenness, it was necessary to use phenologic timing to discriminate between woody and nonwoody vegetation canopies. Initial selection of TM scenes for the study was based on a short preliminary study conducted to identify the time period that would provide the greatest contrast between mesquite and herbaceous canopies. Homogeneous areas of mesquite and herbaceous vegetative cover were sampled from cloud-free TMMSAVI images ranging from January 1, 2002, through December 31, 2011. Mean monthly MSAVI values of the mesquite and herbaceous areas were calculated and plotted for the time series (Fig. 3). Optimal timing for capturing the greatest contrast between the two canopy types was between June and July, with mid- or late- June being most favorable.

Graphic Jump Location
Fig. 3
F3 :

Comparison of MSAVI-based greenness for mesquite and herbaceous canopies over time.

Cloud-free scenes of June were then screened for precipitation. Occurrence of precipitation events was evaluated using data generated from Daymet,48 a collection of algorithms and computer software that interpolates and extrapolates from daily meteorological observations to produce gridded estimates of daily weather parameters over large regions at 1-km spatial resolution. Studies in semiarid environments have demonstrated that germination of annual forbs and grasses can begin at precipitation levels ranging from 10 to 25 mm.4951 This appearance of green vegetation sparked by winter rains is often followed by a subsequent senescent period.52 It is during this senescent period that deep-rooted woody vegetation will green-up, with perennial grasses following after the onset of summer rains (Fig. 3). However, an occurrence of subsequent rainfall events following germination can lead to the presence of a green herbaceous cover extending into the critical window needed for discriminating woody from herbaceous canopies. Moreover, it has been observed that 50.8 mm (2 in.) of winter precipitation can sustain herbaceous vegetation until summer precipitation arrives.53 Therefore, all scenes with a cumulative precipitation (January 1 to the image date) of 50.8 mm or greater were eliminated in an effort to decrease the likelihood that green herbaceous vegetation was present. Finally, TM scene selection was narrowed to years spanning the 2006 to 2011 timeframe to facilitate application with available coincident datasets (Table 2).

Subpixel Discernment of Woody Vegetation and Equation Development

Moderate-resolution satellite imagery is useful for examining vegetation cover over large areas. However, when trying to discriminate between vegetation types, the presence of mixed pixels can be problematic. In this study, a method was needed to “unmix” pixels without having to undertake the more complicated methods used for spectral unmixing.5457 This was accomplished by integrating moderate 30-m resolution imagery with higher 1-m resolution imagery.

Each classified NAIP subset was aggregated up to 30-m resolution using a script created in MATLAB® (Fig. 4). The following is a short description of operations performed within the script for each NAIP scene. The 2011 TMMSAVI and the classified NAIP images were cropped together so that their dimensions ensured each TM pixel region comprised exactly 900 (30×30) NAIP pixels. The cropped NAIP image was then reclassified using a binary system, with each pixel given a value of one for woody or zero for nonwoody. This binary image was aggregated to 30-m resolution by passing a 30×30 pixel moving window over the image and calculating the percent woody cover per iteration [Fig. 4(c)]. Percent woody cover was calculated using the following equation: Display Formula

Percent woody cover=i=1nvin×100,(4)
where vi is a reclassified binary pixel value, n is the total number of pixels in the moving window, and NAIP-based woody cover (NAIPWC) scenes were produced.

Graphic Jump Location
Fig. 4
F4 :

Process for aggregating 1-m resolution NAIP to 30-m resolution: (a) a section of unclassified 1-m NAIP image, (b) section of 1-m NAIP image classified as woody or nonwoody, (c) sample of binary aggregation calculation of percent woody cover using Eq. (2) for a 5×5m moving window where n=25, (d) section of classified NAIP aggregated to 30-m resolution.

The 2011 TMMSAVI and NAIPWC images were used to develop a relationship for producing TM images of woody cover (TMWC). A MATLAB script was created to perform randomly stratified sampling of the TMMSAVI and NAIPWC scenes. Stratification class bounds were based on the NAIPWC values (0 to 100%) and set at 10% increments. The number of samples was set to 10 samples per class without replacement. A 3×3 pixel (90×90m) window was used for sampling the scenes and the mean TMMSAVI and NAIPWC values for the window region were extracted.58 The resulting data were split to serve as training and validation datasets. Regression of the training dataset pairs of TMMSAVI and NAIPWC was used to develop a relation for TMMSAVI-derived woody cover. The resulting equation was applied to the three TMMSAVI scenes to create TMWC scenes. Validation of the TMWC scenes from 2006, 2008, and 2011 was conducted using the NAIPWC, TNC06, and TNC08 datasets.

The objective of this study was to develop an operationally oriented method for assessing woody vegetation cover using publically available, no-cost remotely sensed data sources on semiarid rangelands in southeastern Arizona. The following is a discussion of the results of the method’s development, including accuracy assessments, validation, and end product usage.

It is often necessary to use pre-existing datasets when conducting large-scale land cover assessments. Traditionally, ground-based measures have been the preferred source of reference data when performing accuracy assessments of remotely sensed products. Unfortunately, ground-based data are not only expensive and time consuming to obtain, but when collected are rarely at a scale appropriate for use with moderate-resolution scenes.7,36,39 Alternatively, studies5962 have used high-resolution airborne imagery as a surrogate for ground-based data. When using this method, a subset comparison of ground-collected and airborne data to verify reliability is recommended.39 Therefore, a comparison of high-resolution image-based woody cover and ground-based woody cover is shown (Fig. 5). Good agreement (RMSE=3.8%, MAE=2.9%) was found between the datasets. The error fell within bounds considered reasonable for accepting imagery-based woody cover estimates for use as ground-based surrogates.

Graphic Jump Location
Fig. 5
F5 :

Comparison of high-resolution (NAIP and NRI) imagery-based woody cover and ground-based woody cover. Root-mean-square error (RMSE) and mean absolute error (MAE) are presented.

NAIP imagery was used to produce ground-based woody cover in this study. Accurate classification of this imagery was imperative for generating reasonable estimates of woody cover. Error matrices of the classified NAIP subsets indicated overall accuracies for the subsets that ranged from 93% to 97% (Fig. 6). Misclassification of dark and/or shaded grasses as woody vegetation was a source of error for all subsets. Overclassification of tree clusters was the main source of error for the mesquite dominated scenes (NAIPP2 and NAIPC2).

Graphic Jump Location
Fig. 6
F6 :

Error matrices to assess the accuracy of classifications performed on four NAIP subsets of images from Pima (NAIPP1, NAIPP2) and Cochise (NAIPC1, NAIPC2) counties captured on August 07, 2013 and July 23, 2013, respectively.

The practice of relating spectral vegetation indices to vegetative properties (e.g., canopy cover, LAI, biomass, leaf water content, and chlorophyll) is well documented.6365 In this study, a strong (R2=0.92) linear relationship was found between TMMSAVI and NAIPWC (Fig. 7). The resulting equation: Display Formula

Percent Woody Cover=599.64(TMMSAVI)49.05,(5)
made it possible to quantify woody cover from TM imagery.

Graphic Jump Location
Fig. 7
F7 :

Relationship of 2013 NAIPWC to 2011 TMMSAVI. Regression is significant.

Validation is a critical step in the production of remotely sensed products.39,66,67 Assessment of the 2006, 2008, and 2011 TMWC scenes revealed a number of factors responsible for the resultant error (RMSE=8.5%, MAE=6.4%) (Fig. 8). First, saturation of TMWC corresponding to the NAIPWC(P2) data was likely due to the possible presence of broadleaved woody species mixed with mesquite canopies that were visually undetectable during initial classification efforts. Next, apparent underestimation of TMWC for NAIPWC(P2) data seen in the mid- to high (40% to 85%) range of observed cover was the result of an inflation of NAIPWC caused by overclassification of tree clusters and dark and/or shaded grasses. Further, TMWC was overestimated for NAIPWC(P2) data within the lower (5% to 25%) range of observed cover. This was largely due to a resulting disconnect due to the ability to visually detect the sacaton in the NAIPP2 scene when determining observed woody cover, and an inability of TMMSAVI to differentiate between mesquite and sacaton greenness. Given that green-up of sacaton canopies can start as early as April and last past October, discernment from mesquite greenness was simply unfeasible. Finally, uncertainty in cover estimates when TMWC was less than 10% was largely dependent on how well the woody signal overcame the underlying ground cover signal. Despite the varied nature of these localized errors, the RMSE and MAE of the TMWC scenes were within acceptable bounds for operational use to quantify large-scale woody cover for brush management decisions.

Graphic Jump Location
Fig. 8
F8 :

Comparison of observed (NAIPWC(P1), NAIPWC(P2), NAIPWC(C1), NAIPWC(C2), TNC06, and TNC08) and predicted (2006, 2008, and 2011 TMWC) woody cover. Root-mean-square error (RMSE) and mean absolute error (MAE) are presented.

These large-scale images of woody cover can provide a picture of the spatial distribution of woody cover that is integral in establishing where removal should occur, and in tracking its re-emergence for subsequent removal. For example, subsets of TMWC scenes of pre and postbrush removal on the Empire Ranch showed a decrease in woody cover from 30% to 60% in 2008 to 10% to 20% in 2011 within areas treated with brush removal (Fig. 9). Changes seen in the area southwest of the treatments where woody cover dropped from roughly 11% to 60% down to 0% to 5% cover were caused by a wildfire on June 01, 2011. The continuation of the Landsat time series through the launch of Landsat-8 will make tracking reemergence within these treatments post-2011 possible.

Graphic Jump Location
Fig. 9
F9 :

Map subsets of TMWC on the Empire Ranch for June 22, 2008, and June 15, 2011. Brush removal treatment areas are outlined in black. The years of 2008 and 2011 represent pre- and post-treatment, respectively.

An operationally oriented protocol for quantifying woody vegetation cover with no-cost, publically available imagery is needed to aid land managers in implementing, monitoring, and assessing brush management conservation practice applications. This work showed that it was possible to produce viable (RMSE=8.5%, MAE=6.4%) maps of woody cover to track the occurrence of brush removal and monitor the presence or lack of subsequent re-emergence. However, there are a number of conditions that must be upheld to achieve valid results. With regard to Landsat image selection, appropriate phenologic timing for capturing the greatest contrast between the woody species of interest and other vegetative cover is critical. Precipitation and cloud cover must also be considered. In this study, the woody species of interest was mesquite and optimal timing was mid- or late- June, with less than 50.8 mm of cumulative precipitation occurring before the cloud-free image date. Use of the protocol in other locations (e.g., Texas, New Mexico, Colorado, etc.) will require new determinations of the optimal timing to allow for greatest contrast between woody and nonwoody vegetation, as well as precipitation conditions that limit the occurrence of green forbs and grasses.

Generation and validation of TM scenes to estimate woody cover from TMMSAVI requires corresponding image-based woody cover values. Given that these image-based woody cover estimates are calculated from classifying high-resolution imagery, accurate classification of woody versus nonwoody vegetation is essential. Misclassification due to shadowing and clumping of vegetation will manifest as error in the validation of Landsat-based woody cover estimates. Early green-up of herbaceous vegetation will also contribute to misclassification error and should be avoided if possible. The relation for woody cover from MSAVI in this study (Eq. 5) was developed to target mesquite cover. It does detect other woody species; however, to achieve the best results, new relations should be developed for use with woody species like Juniperus L. (juniper) or Populus L. (cottonwood) that have reflective signatures that are very different from mesquite.

This study demonstrated that it was possible to produce a viable method for developing an operationally oriented method for assessing woody vegetation using publically available, no-cost remotely sensed data on semiarid rangelands in southeastern Arizona. Future work will include development of new relationships to better quantify juniper and cottonwood species, as well as test protocol applicability to Landsat-8 imagery and other rangeland locations in the central and western United States.

This research was supported by NRCS Interagency Agreement 67-3A75-14-115. The authors wish to express sincere thanks to Mary White, Jason Wong, David Goodrich, and the personnel at the USDA-ARS Southwest Watershed Research Center, Emilio Carrillo from the USDA-NRCS area office in Tucson, Arizona, Lee Norfleet from the USDA-NRCS Resource Assessment Division (CEAP Modeling Team) in Temple, Texas, Ken Laterza and Tony Garcia from the USDA-NRCS Central Remote Sensing Laboratory in Ft. Worth, Texas, Veronica Lessard from the Iowa State University’s Center for Survey Statistics and Methodology, and Stuart Marsh, Willem Van Leeuwen, Kyle Hartfield, and Andy Honaman from the University of Arizona, for their cooperation and assistance. Mention of proprietary product does not constitute a guarantee or warranty of the product by USDA or the author and does not imply its approval to the exclusion of other products that may also be suitable.

Archer  S. R.  et al., “Brush management as a rangeland conservation strategy: a critical evaluation,” in Conservation Benefits of Rangeland Practices: Assessment, Recommendations, and Knowledge Gaps. , Briske  D. D., Ed., pp. 107 ,  United States Department of Agriculture ,  Natural Resources Conservation Service  (2011).
Barger  N. N.  et al., “Woody plant proliferation in North American drylands: a synthesis of impacts on ecosystem carbon balance,” J. Geophys. Res.. 116, , 1 –17 (2011) G00K07, . 0148-0227 CrossRef
Knapp  A. K.  et al., “Shrub encroachment in North American grasslands: shifts in growth form dominance rapidly alters control of ecosystem carbon inputs,” Global Change Biol.. 14, , 615 –623 (2008). 1354-1013 CrossRef
Archer  S., , Boutton  T. W., and Hibbard  K. A., “Trees in grasslands: biogeochemical consequences of woody plant expansion,” in Global Biogeochemical Cycles in the Climate System. , Schulze  M., , Heimann  S., , Harrison  E., , Holland  J., and Lloyd  I., Eds., pp. 115 –138,  Academic Press ,  San Diego, California  (2001).
Welch  T. G., , Smith  R. P., and Rasmussen  G. A., “Brush management technologies,” in Integrated Brush Management Systems for South Texas: Development and Implementation. , Vol. 1493, pp. 71 ,  Agricultural Experimental Station Bulletin ,  College Station, Texas  (1985).
Briske  D. D.  et al., “Introduction to the conservation effects assessment project and the rangeland literature synthesis,” in Conservation Benefits of Rangeland Practices: Assessment, Recommendations, and Knowledge Gaps. , Briske  D. D., Ed., pp. 4 –8,  United States Department of Agriculture, Natural Resources Conservation Service ,  Lawrence, Kansas  (2011).
Pellant  M., , Shaver  P., and Spaeth  K., “Field test of a prototype rangeland inventory procedure in the western USA,” in People and Rangelands: Building the Future. , Eldridge  D., and Freudenberger  D., Eds.,  Proc. VI Int. Rangeland Congr.  Vol. 4814, pp. 766 –767,  The Congress ,  Townsville, Queensland  (1999).
Graetz  R. D.  et al., “The application of Landsat image data to rangeland assessment and monitoring: an example from South Australia,” Aust. Rangel. J.. 5, (2 ), 63 –73 (1983).CrossRef
Graetz  R. D., , Pech  R. P., and Davis  A. W., “The assessment and monitoring of sparsely vegetated rangelands using calibrated Landsat data,” Int. J. Remote Sens.. 9, (7 ), 1201 –1222 (1988). 0143-1161 CrossRef
Pickup  G., , Bastin  G. N., and Chewings  V. H., “Remote-sensing-based condition assessment for nonequilibrium rangelands under large-scale commercial grazing,” Ecol. Appl.. 4, (3 ), 497 –517 (1994). 1051-0761 CrossRef
Elmore  A. J.  et al., “Quantifying vegetation change in semiarid environments: precision and accuracy of spectral mixture analysis and the normalized difference vegetation index,” Remote Sens. Environ.. 73, , 87 –102 (2000). 0034-4257 CrossRef
Muchoney  D.  et al., “Application of the MODIS global supervised classification model to vegetation and land cover mapping of Central America,” Int. J. Remote Sens.. 21, (6–7 ), 1115 –1138 (2000). 0143-1161 CrossRef
Hunt  E.  Jr.  et al., “Applications and research using remote sensing for rangeland management,” Photogramm. Eng. Remote Sens.. 69, (6 ), 675 –693 (2003). 0099-1112 CrossRef
Geerken  R., and Ilaiwi  M., “Assessment of rangeland degradation and development of a strategy for rehabilitation,” Remote Sens. Environ.. 90, (4 ), 490 –504 (2004). 0034-4257 CrossRef
Marsett  R. C.  et al., “Remote sensing for grassland management in the arid southwest,” Rangeland Ecol. Manag.. 59, (5 ), 530 –540 (2006). 1550-7424 CrossRef
Verbesselt  J.  et al., “Detecting trend and seasonal changes in satellite image time series,” Remote Sens. Environ.. 114, , 106 –115 (2010). 0034-4257 CrossRef
Esch  T.  et al., “Combined use of multiseasonal high and medium resolution satellite imagery for parcel related mapping of cropland and grassland,” Int. J. Appl. Earth Obs. Geoinf.. 28, , 230 –237 (2014). 0303-2434 CrossRef
Asner  G. P.  et al., “Net changes in regional woody vegetation cover and carbon storage in Texas drylands, 1937–1999,” Global Change Biol.. 9, , 316 –335 (2003). 1354-1013 CrossRef
DeFries  R. S.  et al., “A new global 1-km dataset of percentage tree cover derived from remote sensing,” Global Change Biol.. 6, , 247 –254 (2000). 1354-1013 CrossRef
Mirik  M., and Ansley  R. J., “Utility of satellite and aerial images for quantification of canopy cover and infilling rates of the invasive woody species honey mesquite (Prosopis glandulosa) on rangeland,” Remote Sens.. 4, , 1947 –1962 (2012). 2072-4292 CrossRef
Moleele  N. M.  et al., “More woody plants? the status of bush encroachment in Botswana’s grazing areas,” J. Environ. Manage.. 64, , 3 –11 (2002). 0301-4797 CrossRef
Belay  T. A., , Totland  Ø., and Moe  S. R., “Woody vegetation dynamics in the rangelands of lower Omo region, southwestern Ethiopia,” J. Arid Environ.. 89, , 94 –102 (2013).CrossRef
Hansen  M. C.  et al., “Towards an operational MODIS continuous field of percent tree cover algorithm: examples using AVHRR and MODIS data,” Remote Sens. Environ.. 83, , 303 –319 (2002). 0034-4257 CrossRef
Laliberte  A. S.  et al., “Object-oriented image analysis for mapping shrub encroachment from 1937 to 2003 in southern New Mexico,” Remote Sens. Environ.. 93, , 198 –210 (2004). 0034-4257 CrossRef
Sohn  Y., and McCoy  R. M., “Mapping desert shrub rangeland using spectral unmixing and modeling spectral mixtures with TM data,” Photogramm. Eng. Remote Sens.. 63, (6 ), 707 –716 (1997). 0099-1112 
Hamada  Y.  et al., “Assessing and monitoring semi-arid shrublands using object-based image analysis and multiple endmember spectral mixture analysis,” Environ. Monit. Assess.. 185, , 3173 –3190 (2013). 0167-6369 CrossRef
Hamada  Y., , Stow  D. A., and Roberts  D. A., “Estimating life-form cover fractions in California sage scrub communities using multispectral remote sensing,” Remote Sens. Environ.. 115, , 3056 –3068 (2011). 0034-4257 CrossRef
Kuemmerle  T., , Röder  A., and Hill  J., “Separating grassland and shrub vegetation by multidate pixel-adaptive spectral mixture analysis,” Int. J. Remote Sens.. 27, (15 ), 3251 –3271 (2006). 0143-1161 CrossRef
Okin  G. S.  et al., “Practical limits on hyperspectral vegetation discrimination in arid and semiarid environments,” Remote Sens. Environ.. 77, , 212 –225 (2001). 0034-4257 CrossRef
Davies  K. W.  et al., “Estimating juniper cover from National Agriculture Imagery Program (NAIP) imagery and evaluating relationships between potential cover and environmental variables,” Rangeland Ecol. Manage.. 63, , 630 –637 (2010). 1550-7424 CrossRef
United States Department of Agriculture, Natural Resources Conservation Service, Land Resource Regions and Major Land Resource Areas of the United States, the Caribbean, and the Pacific Basin. , pp. 107 , U.S. Department of Agriculture Handbook 296 (2006).
Herrick  J. E.  et al., Monitoring Manual for Grassland, Shrubland and Savanna Ecosystems. Volume 1: Quick Start. , USDA-ARS Jornada Experimental Range,  The University of Arizona Press ,  Tucson  (2005).
Ansley  R. J., , Wu  X. B., and Kramp  B. A., “Observation: long-term increases in mesquite canopy cover in a North Texas savanna,” J. Range Manage.. 54, , 171 –176 (2001). 0022-409X CrossRef
Warren  A., , Holechek  J., and Cardenas  M., “Honey mesquite influences on Chihuahuan desert vegetation,” J. Range Manage.. 49, , 46 –52 (1996). 0022-409X CrossRef
Archer  S.  et al., “Autogenic succession in a subtropical savanna: conversion of grassland to thorn woodland,” Ecol. Monogr.. 58, (2 ), 111 –127 (1988). 0012-9615 CrossRef
Nusser  S. M., and Goebel  J. J., “The national resources Inventory: a long-term multi-resource monitoring programme,” Environ. Ecol. Stat.. 4, , 181 –204 (1997). 1352-8505 CrossRef
Spaeth  K. E.  et al., “New proposed national resources inventory protocols on nonfederal rangeland,” J. Soil Water Conserv.. 58, (1 ), 18A –21A (2003). 0022-4561 
Congalton  R. G., and Green  K., Assessing the Accuracy of Remotely Sensed Data. ,  CRC Press ,  Boca Raton, Florida  (2008).
Masek  J. G.  et al., “A Landsat surface reflectance data set for North America,” IEEE Geosci. Remote Sens. Lett.. 3, , 68 –72 (2006). 1545-598X CrossRef
Masek  J. G.  et al., LEDAPS Calibration, Reflectance, Atmospheric Correction Preprocessing Code, Version 2. Model Product.. http://daac.ornl.gov] from Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee (2013).
Huete  A. R., “A soil-adjusted vegetation index (SAVI),” Remote Sens. Environ.. 25, , 295 –309 (1988). 0034-4257 CrossRef
Rondeaux  G., , Steven  M., and Baret  F., “Optimization of soil-adjusted vegetation indices,” Remote Sens. Environ.. 55, , 95 –107 (1996). 0034-4257 CrossRef
Pettorelli  N.  et al., “Using the satellite-derived NDVI to assess ecological responses to environmental change,” Trends Ecol. Evol.. 20, (9 ), 503 –510 (2005). 0169-5347 CrossRef
Baret  F., and Guyot  G., “Potentials and limits of vegetation indices for LAI and APAR assessment,” Remote Sens. Environ.. 35, , 161 –173 (1991). 0034-4257 CrossRef
Qi  J.  et al., “A modified soil adjusted vegetation index,” Remote Sens. Environ.. 48, (2 ), 119 –126 (1994). 0034-4257 CrossRef
Thornton  P. E., , Running  S. W., and White  M. A., “Generating surfaces of daily meteorology variables over large regions of complex terrain,” J. Hydrol.. 190, , 214 –251 (1997). 0022-1694 CrossRef
Juhren  M., , Went  F. W., and Phillips  E., “Ecology of desert plants. IV. Combined field and laboratory work on germination of annuals in the Joshua Tree National Monument, California,” Ecol.. 37, (2 ), 318 –330 (1956). 0012-9658 CrossRef
Halvorson  W. L., and Patten  D. T., “Productivity and flowering of winter ephemerals in relation to Sonoran desert shrubs,” Am. Midland Nat.. 93, (2 ), 311 –319 (1975). 0003-0031 CrossRef
Bowers  J. E., “El Niño and displays of spring-flowering annuals in the Mojave and Sonoran deserts,” J. Torrey Bot. Soc.. 132, (1 ), 38 –49 (2005). 1095-5674 CrossRef
Svoray  T., and Karnieli  A., “Rainfall, topography and primary production relationships in a semiarid ecosystem,” Ecohydrology. 4, (1 ), 56 –66 (2011). 1936-0584 CrossRef
Emilio Carrillo, personal communication, USDA Natural Resources Conservation Service, Tucson, Arizona (2015).
Asner  G. P., and Heidebrecht  K. B., “Spectral unmixing of vegetation, soil and dry carbon cover in arid regions: comparing multispectral and hyperspectral observations,” Int. J. Remote Sens.. 23, (19 ), 3939 –3958 (2002). 0143-1161 CrossRef
Huang  C., and Townshend  J. R. G., “A stepwise regression tree for nonlinear approximation: applications to estimating subpixel land cover,” Int. J. Remote Sens.. 24, (1 ), 75 –90 (2003). 0143-1161 CrossRef
Foody  G. M.  et al., “Non-linear mixture modeling without end-members using an artificial neural network,” Int. J. Remote Sens.. 18, (4 ), 937 –953 (1997). 0143-1161 CrossRef
Schmitt-Harsh  M., , Sweeney  S. P., and Evans  T. P., “Classification of coffee-forest landscapes using Landsat TM imagery and spectral mixture analysis,” Photogramm. Eng. Remote Sens.. 79, (5 ), 457 –468 (2013). 0099-1112 CrossRef
Stehman  S. V., and Czaplewski  R. L., “Design and analysis for thematic map accuracy assessment: fundamental principles,” Remote Sens. Environ.. 64, (3 ), 331 –344 (1998). 0034-4257 CrossRef
Justice  C.  et al., “Developments in the ‘validation’ of satellite sensor products for the study of the land surface,” Int. J. Remote Sens.. 21, , 3383 –3390 (2000). 0143-1161 CrossRef
Xian  G., and Homer  C., “Updating the 2001 national land cover database impervious surface products to 2006 using Landsat imagery change detection methods,” Remote Sens. Environ.. 114, , 1676 –1686 (2010). 0034-4257 CrossRef
Xiao  J., and Moody  A., “A comparison of methods for estimating fractional green vegetation cover within a desert-to-upland transition zone in central New Mexico, USA,” Remote Sens. Environ.. 98, , 237 –250 (2005). 0034-4257 CrossRef
Wulder  M. A.  et al., “Monitoring Canada’s forests. Part 1: completion of the EOSD land cover project,” Can. J. Remote Sens.. 34, (6 ), 549 –562 (2008). 0703-8992 CrossRef
Adam  E., , Mutanga  O., and Rugege  D., “Multispectral and hyperspectral remote sensing for identification and mapping of wetland vegetation: a review,” Wetlands Ecol. Manage.. 18, , 281 –296 (2010).CrossRef
Purevdorj  T.  et al., “Relationships between percent vegetation cover and vegetation indices,” Int. J. Remote Sens.. 19, (18 ), 3519 –3535 (1998). 0143-1161 CrossRef
Hagen  S. C.  et al., “Mapping total vegetation cover across western rangelands with moderate-resolution imaging spectroradiometer data,” Rangeland Ecol. Manag.. 65, (5 ), 456 –467 (2012). 1550-7424 CrossRef
Foody  G. M., “Status of land cover classification accuracy assessment,” Remote Sens. Environ.. 80, , 185 –201 (2002). 0034-4257 CrossRef
Congalton  R. G., “A review of assessing the accuracy of classifications of remotely sensed data,” Remote Sens. Environ.. 37, , 35 –46 (1991). 0034-4257 CrossRef
© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Chandra D. Holifield Collins ; Mark A. Kautz ; Ronald Tiller ; Sapana Lohani ; Guillermo Ponce-Campos, et al.
"Development of an integrated multiplatform approach for assessing brush management conservation efforts in semiarid rangelands", J. Appl. Remote Sens. 9(1), 096057 (May 07, 2015). ; http://dx.doi.org/10.1117/1.JRS.9.096057


Figures

Graphic Jump Location
Fig. 1
F1 :

Major land resource area 41 overlaid with a Landsat-5 TM tile for path 35, row 38, and the Empire Ranch study area with National Agricultural Imagery Program (NAIP), National Resources Inventory (NRI), The Nature Conservancy (TNC), brush removal, and Empire Ranch sampling block (ERBLOCK) study sites. NAIP and NRI represent high-resolution imagery subset locations. TNC and ERBLOCK represent ground-based sampling locations.

Graphic Jump Location
Fig. 2
F2 :

Examples of NAIP subsets from Pima and Cochise counties collected on August 07, 2013 and July 23, 2013, respectively.

Graphic Jump Location
Fig. 3
F3 :

Comparison of MSAVI-based greenness for mesquite and herbaceous canopies over time.

Graphic Jump Location
Fig. 4
F4 :

Process for aggregating 1-m resolution NAIP to 30-m resolution: (a) a section of unclassified 1-m NAIP image, (b) section of 1-m NAIP image classified as woody or nonwoody, (c) sample of binary aggregation calculation of percent woody cover using Eq. (2) for a 5×5m moving window where n=25, (d) section of classified NAIP aggregated to 30-m resolution.

Graphic Jump Location
Fig. 5
F5 :

Comparison of high-resolution (NAIP and NRI) imagery-based woody cover and ground-based woody cover. Root-mean-square error (RMSE) and mean absolute error (MAE) are presented.

Graphic Jump Location
Fig. 9
F9 :

Map subsets of TMWC on the Empire Ranch for June 22, 2008, and June 15, 2011. Brush removal treatment areas are outlined in black. The years of 2008 and 2011 represent pre- and post-treatment, respectively.

Graphic Jump Location
Fig. 8
F8 :

Comparison of observed (NAIPWC(P1), NAIPWC(P2), NAIPWC(C1), NAIPWC(C2), TNC06, and TNC08) and predicted (2006, 2008, and 2011 TMWC) woody cover. Root-mean-square error (RMSE) and mean absolute error (MAE) are presented.

Graphic Jump Location
Fig. 7
F7 :

Relationship of 2013 NAIPWC to 2011 TMMSAVI. Regression is significant.

Graphic Jump Location
Fig. 6
F6 :

Error matrices to assess the accuracy of classifications performed on four NAIP subsets of images from Pima (NAIPP1, NAIPP2) and Cochise (NAIPC1, NAIPC2) counties captured on August 07, 2013 and July 23, 2013, respectively.

Tables

Table Grahic Jump Location
Table 1Ground-based dataset description. NRI=National Resources Inventory, TNC=The Nature Conservancy, ERBLOCK=Empire Ranch sampling block.
Table Grahic Jump Location
Table 2Imagery dataset description. NRI=National Resources Inventory, NAIP=National Agricultural Imagery Program, TM=Landsat-5Thematic Mapper.

References

Archer  S. R.  et al., “Brush management as a rangeland conservation strategy: a critical evaluation,” in Conservation Benefits of Rangeland Practices: Assessment, Recommendations, and Knowledge Gaps. , Briske  D. D., Ed., pp. 107 ,  United States Department of Agriculture ,  Natural Resources Conservation Service  (2011).
Barger  N. N.  et al., “Woody plant proliferation in North American drylands: a synthesis of impacts on ecosystem carbon balance,” J. Geophys. Res.. 116, , 1 –17 (2011) G00K07, . 0148-0227 CrossRef
Knapp  A. K.  et al., “Shrub encroachment in North American grasslands: shifts in growth form dominance rapidly alters control of ecosystem carbon inputs,” Global Change Biol.. 14, , 615 –623 (2008). 1354-1013 CrossRef
Archer  S., , Boutton  T. W., and Hibbard  K. A., “Trees in grasslands: biogeochemical consequences of woody plant expansion,” in Global Biogeochemical Cycles in the Climate System. , Schulze  M., , Heimann  S., , Harrison  E., , Holland  J., and Lloyd  I., Eds., pp. 115 –138,  Academic Press ,  San Diego, California  (2001).
Welch  T. G., , Smith  R. P., and Rasmussen  G. A., “Brush management technologies,” in Integrated Brush Management Systems for South Texas: Development and Implementation. , Vol. 1493, pp. 71 ,  Agricultural Experimental Station Bulletin ,  College Station, Texas  (1985).
Briske  D. D.  et al., “Introduction to the conservation effects assessment project and the rangeland literature synthesis,” in Conservation Benefits of Rangeland Practices: Assessment, Recommendations, and Knowledge Gaps. , Briske  D. D., Ed., pp. 4 –8,  United States Department of Agriculture, Natural Resources Conservation Service ,  Lawrence, Kansas  (2011).
Pellant  M., , Shaver  P., and Spaeth  K., “Field test of a prototype rangeland inventory procedure in the western USA,” in People and Rangelands: Building the Future. , Eldridge  D., and Freudenberger  D., Eds.,  Proc. VI Int. Rangeland Congr.  Vol. 4814, pp. 766 –767,  The Congress ,  Townsville, Queensland  (1999).
Graetz  R. D.  et al., “The application of Landsat image data to rangeland assessment and monitoring: an example from South Australia,” Aust. Rangel. J.. 5, (2 ), 63 –73 (1983).CrossRef
Graetz  R. D., , Pech  R. P., and Davis  A. W., “The assessment and monitoring of sparsely vegetated rangelands using calibrated Landsat data,” Int. J. Remote Sens.. 9, (7 ), 1201 –1222 (1988). 0143-1161 CrossRef
Pickup  G., , Bastin  G. N., and Chewings  V. H., “Remote-sensing-based condition assessment for nonequilibrium rangelands under large-scale commercial grazing,” Ecol. Appl.. 4, (3 ), 497 –517 (1994). 1051-0761 CrossRef
Elmore  A. J.  et al., “Quantifying vegetation change in semiarid environments: precision and accuracy of spectral mixture analysis and the normalized difference vegetation index,” Remote Sens. Environ.. 73, , 87 –102 (2000). 0034-4257 CrossRef
Muchoney  D.  et al., “Application of the MODIS global supervised classification model to vegetation and land cover mapping of Central America,” Int. J. Remote Sens.. 21, (6–7 ), 1115 –1138 (2000). 0143-1161 CrossRef
Hunt  E.  Jr.  et al., “Applications and research using remote sensing for rangeland management,” Photogramm. Eng. Remote Sens.. 69, (6 ), 675 –693 (2003). 0099-1112 CrossRef
Geerken  R., and Ilaiwi  M., “Assessment of rangeland degradation and development of a strategy for rehabilitation,” Remote Sens. Environ.. 90, (4 ), 490 –504 (2004). 0034-4257 CrossRef
Marsett  R. C.  et al., “Remote sensing for grassland management in the arid southwest,” Rangeland Ecol. Manag.. 59, (5 ), 530 –540 (2006). 1550-7424 CrossRef
Verbesselt  J.  et al., “Detecting trend and seasonal changes in satellite image time series,” Remote Sens. Environ.. 114, , 106 –115 (2010). 0034-4257 CrossRef
Esch  T.  et al., “Combined use of multiseasonal high and medium resolution satellite imagery for parcel related mapping of cropland and grassland,” Int. J. Appl. Earth Obs. Geoinf.. 28, , 230 –237 (2014). 0303-2434 CrossRef
Asner  G. P.  et al., “Net changes in regional woody vegetation cover and carbon storage in Texas drylands, 1937–1999,” Global Change Biol.. 9, , 316 –335 (2003). 1354-1013 CrossRef
DeFries  R. S.  et al., “A new global 1-km dataset of percentage tree cover derived from remote sensing,” Global Change Biol.. 6, , 247 –254 (2000). 1354-1013 CrossRef
Mirik  M., and Ansley  R. J., “Utility of satellite and aerial images for quantification of canopy cover and infilling rates of the invasive woody species honey mesquite (Prosopis glandulosa) on rangeland,” Remote Sens.. 4, , 1947 –1962 (2012). 2072-4292 CrossRef
Moleele  N. M.  et al., “More woody plants? the status of bush encroachment in Botswana’s grazing areas,” J. Environ. Manage.. 64, , 3 –11 (2002). 0301-4797 CrossRef
Belay  T. A., , Totland  Ø., and Moe  S. R., “Woody vegetation dynamics in the rangelands of lower Omo region, southwestern Ethiopia,” J. Arid Environ.. 89, , 94 –102 (2013).CrossRef
Hansen  M. C.  et al., “Towards an operational MODIS continuous field of percent tree cover algorithm: examples using AVHRR and MODIS data,” Remote Sens. Environ.. 83, , 303 –319 (2002). 0034-4257 CrossRef
Laliberte  A. S.  et al., “Object-oriented image analysis for mapping shrub encroachment from 1937 to 2003 in southern New Mexico,” Remote Sens. Environ.. 93, , 198 –210 (2004). 0034-4257 CrossRef
Sohn  Y., and McCoy  R. M., “Mapping desert shrub rangeland using spectral unmixing and modeling spectral mixtures with TM data,” Photogramm. Eng. Remote Sens.. 63, (6 ), 707 –716 (1997). 0099-1112 
Hamada  Y.  et al., “Assessing and monitoring semi-arid shrublands using object-based image analysis and multiple endmember spectral mixture analysis,” Environ. Monit. Assess.. 185, , 3173 –3190 (2013). 0167-6369 CrossRef
Hamada  Y., , Stow  D. A., and Roberts  D. A., “Estimating life-form cover fractions in California sage scrub communities using multispectral remote sensing,” Remote Sens. Environ.. 115, , 3056 –3068 (2011). 0034-4257 CrossRef
Kuemmerle  T., , Röder  A., and Hill  J., “Separating grassland and shrub vegetation by multidate pixel-adaptive spectral mixture analysis,” Int. J. Remote Sens.. 27, (15 ), 3251 –3271 (2006). 0143-1161 CrossRef
Okin  G. S.  et al., “Practical limits on hyperspectral vegetation discrimination in arid and semiarid environments,” Remote Sens. Environ.. 77, , 212 –225 (2001). 0034-4257 CrossRef
Davies  K. W.  et al., “Estimating juniper cover from National Agriculture Imagery Program (NAIP) imagery and evaluating relationships between potential cover and environmental variables,” Rangeland Ecol. Manage.. 63, , 630 –637 (2010). 1550-7424 CrossRef
United States Department of Agriculture, Natural Resources Conservation Service, Land Resource Regions and Major Land Resource Areas of the United States, the Caribbean, and the Pacific Basin. , pp. 107 , U.S. Department of Agriculture Handbook 296 (2006).
Herrick  J. E.  et al., Monitoring Manual for Grassland, Shrubland and Savanna Ecosystems. Volume 1: Quick Start. , USDA-ARS Jornada Experimental Range,  The University of Arizona Press ,  Tucson  (2005).
Ansley  R. J., , Wu  X. B., and Kramp  B. A., “Observation: long-term increases in mesquite canopy cover in a North Texas savanna,” J. Range Manage.. 54, , 171 –176 (2001). 0022-409X CrossRef
Warren  A., , Holechek  J., and Cardenas  M., “Honey mesquite influences on Chihuahuan desert vegetation,” J. Range Manage.. 49, , 46 –52 (1996). 0022-409X CrossRef
Archer  S.  et al., “Autogenic succession in a subtropical savanna: conversion of grassland to thorn woodland,” Ecol. Monogr.. 58, (2 ), 111 –127 (1988). 0012-9615 CrossRef
Nusser  S. M., and Goebel  J. J., “The national resources Inventory: a long-term multi-resource monitoring programme,” Environ. Ecol. Stat.. 4, , 181 –204 (1997). 1352-8505 CrossRef
Spaeth  K. E.  et al., “New proposed national resources inventory protocols on nonfederal rangeland,” J. Soil Water Conserv.. 58, (1 ), 18A –21A (2003). 0022-4561 
Congalton  R. G., and Green  K., Assessing the Accuracy of Remotely Sensed Data. ,  CRC Press ,  Boca Raton, Florida  (2008).
Masek  J. G.  et al., “A Landsat surface reflectance data set for North America,” IEEE Geosci. Remote Sens. Lett.. 3, , 68 –72 (2006). 1545-598X CrossRef
Masek  J. G.  et al., LEDAPS Calibration, Reflectance, Atmospheric Correction Preprocessing Code, Version 2. Model Product.. http://daac.ornl.gov] from Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee (2013).
Huete  A. R., “A soil-adjusted vegetation index (SAVI),” Remote Sens. Environ.. 25, , 295 –309 (1988). 0034-4257 CrossRef
Rondeaux  G., , Steven  M., and Baret  F., “Optimization of soil-adjusted vegetation indices,” Remote Sens. Environ.. 55, , 95 –107 (1996). 0034-4257 CrossRef
Pettorelli  N.  et al., “Using the satellite-derived NDVI to assess ecological responses to environmental change,” Trends Ecol. Evol.. 20, (9 ), 503 –510 (2005). 0169-5347 CrossRef
Baret  F., and Guyot  G., “Potentials and limits of vegetation indices for LAI and APAR assessment,” Remote Sens. Environ.. 35, , 161 –173 (1991). 0034-4257 CrossRef
Qi  J.  et al., “A modified soil adjusted vegetation index,” Remote Sens. Environ.. 48, (2 ), 119 –126 (1994). 0034-4257 CrossRef
Thornton  P. E., , Running  S. W., and White  M. A., “Generating surfaces of daily meteorology variables over large regions of complex terrain,” J. Hydrol.. 190, , 214 –251 (1997). 0022-1694 CrossRef
Juhren  M., , Went  F. W., and Phillips  E., “Ecology of desert plants. IV. Combined field and laboratory work on germination of annuals in the Joshua Tree National Monument, California,” Ecol.. 37, (2 ), 318 –330 (1956). 0012-9658 CrossRef
Halvorson  W. L., and Patten  D. T., “Productivity and flowering of winter ephemerals in relation to Sonoran desert shrubs,” Am. Midland Nat.. 93, (2 ), 311 –319 (1975). 0003-0031 CrossRef
Bowers  J. E., “El Niño and displays of spring-flowering annuals in the Mojave and Sonoran deserts,” J. Torrey Bot. Soc.. 132, (1 ), 38 –49 (2005). 1095-5674 CrossRef
Svoray  T., and Karnieli  A., “Rainfall, topography and primary production relationships in a semiarid ecosystem,” Ecohydrology. 4, (1 ), 56 –66 (2011). 1936-0584 CrossRef
Emilio Carrillo, personal communication, USDA Natural Resources Conservation Service, Tucson, Arizona (2015).
Asner  G. P., and Heidebrecht  K. B., “Spectral unmixing of vegetation, soil and dry carbon cover in arid regions: comparing multispectral and hyperspectral observations,” Int. J. Remote Sens.. 23, (19 ), 3939 –3958 (2002). 0143-1161 CrossRef
Huang  C., and Townshend  J. R. G., “A stepwise regression tree for nonlinear approximation: applications to estimating subpixel land cover,” Int. J. Remote Sens.. 24, (1 ), 75 –90 (2003). 0143-1161 CrossRef
Foody  G. M.  et al., “Non-linear mixture modeling without end-members using an artificial neural network,” Int. J. Remote Sens.. 18, (4 ), 937 –953 (1997). 0143-1161 CrossRef
Schmitt-Harsh  M., , Sweeney  S. P., and Evans  T. P., “Classification of coffee-forest landscapes using Landsat TM imagery and spectral mixture analysis,” Photogramm. Eng. Remote Sens.. 79, (5 ), 457 –468 (2013). 0099-1112 CrossRef
Stehman  S. V., and Czaplewski  R. L., “Design and analysis for thematic map accuracy assessment: fundamental principles,” Remote Sens. Environ.. 64, (3 ), 331 –344 (1998). 0034-4257 CrossRef
Justice  C.  et al., “Developments in the ‘validation’ of satellite sensor products for the study of the land surface,” Int. J. Remote Sens.. 21, , 3383 –3390 (2000). 0143-1161 CrossRef
Xian  G., and Homer  C., “Updating the 2001 national land cover database impervious surface products to 2006 using Landsat imagery change detection methods,” Remote Sens. Environ.. 114, , 1676 –1686 (2010). 0034-4257 CrossRef
Xiao  J., and Moody  A., “A comparison of methods for estimating fractional green vegetation cover within a desert-to-upland transition zone in central New Mexico, USA,” Remote Sens. Environ.. 98, , 237 –250 (2005). 0034-4257 CrossRef
Wulder  M. A.  et al., “Monitoring Canada’s forests. Part 1: completion of the EOSD land cover project,” Can. J. Remote Sens.. 34, (6 ), 549 –562 (2008). 0703-8992 CrossRef
Adam  E., , Mutanga  O., and Rugege  D., “Multispectral and hyperspectral remote sensing for identification and mapping of wetland vegetation: a review,” Wetlands Ecol. Manage.. 18, , 281 –296 (2010).CrossRef
Purevdorj  T.  et al., “Relationships between percent vegetation cover and vegetation indices,” Int. J. Remote Sens.. 19, (18 ), 3519 –3535 (1998). 0143-1161 CrossRef
Hagen  S. C.  et al., “Mapping total vegetation cover across western rangelands with moderate-resolution imaging spectroradiometer data,” Rangeland Ecol. Manag.. 65, (5 ), 456 –467 (2012). 1550-7424 CrossRef
Foody  G. M., “Status of land cover classification accuracy assessment,” Remote Sens. Environ.. 80, , 185 –201 (2002). 0034-4257 CrossRef
Congalton  R. G., “A review of assessing the accuracy of classifications of remotely sensed data,” Remote Sens. Environ.. 37, , 35 –46 (1991). 0034-4257 CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.