1. Introduction
Deserts cover nearly one third of the Earth’s land surface [
1,
2]. The biome supports surprisingly high biodiversity despite stressful environmental conditions. Desert environments are sensitive to environmental shifts, and reversing environmental impacts can take decades or even centuries [
3]. Many desert lands in the world undergo significant changes caused by various forces, including climate change [
4], exotic species invasion [
5], over grazing [
6], excessive recreational use [
7], and extensive development [
8]. In the United States, the desert southwest region has been experiencing considerable pressure from energy development, such as pipelines, transmission lines, wind farms, and solar energy facilities. Solar energy development in the United States arid desert areas has been a particular focus over the past decade, and guidance has been developed to mitigate impacts from solar development on public lands [
9]. Despite the benefits of solar energy for decreasing greenhouse gas emissions in comparison with fossil-fuel based energy sources [
8], some impacts of large-scale solar energy development on desert environments are of concern because of the large development footprint and potential to impact microclimates and land cover [
10,
11]. Degradation of desert ecosystems is a serious concern in other parts of the world. For example, exotic species invasion has led to significant reduction or extinction of native plant species in the Middle Eastern deserts [
5]. Plant diversity and productivity have severely declined due to sheep grazing in the desert steppe of northern China [
6]. To combat the widespread degradation of desert lands, develo** the ability to detect early signs of habitat disturbance and to monitor changes efficiently over time is critical to ensuring mitigation of the initial impact and prevention of further degradation [
12].
Researchers and land resource managers strategize to monitor desert-land changes effectively and sustainably. In the United States, the U.S. Department of Interior, Bureau of Land Management (BLM) has established standardized monitoring protocols based on field surveys, known as the Assessment, Inventory, and Monitoring (AIM) protocol [
13]. The monitoring protocols primarily consist of systematic field surveys of plant functional groups, canopy height, presence of non-native species, and canopy gaps. These field data are used to calculate fractional cover to indicate landscape conditions. A group of researchers in Europe also developed a standardized procedure to monitor 130 habitat types, known as General Habitat Categories, or GHCs, including the desert biome [
14]. This landscape-level monitoring protocol adopted the concept of plant life-form types that has been widely applied in biogeography since the 19th century. Visually estimated fractional cover of plant life-form types or species (when cover >30%) is recorded in areas larger than 400 m
2. This protocol aims to inform policy making and research objective decisions and has been rigorously tested at field sites across Europe.
Many researchers and land managers, including the group supporting AIM monitoring efforts, have been actively looking into applications of the use of remote sensing in conjunction with field survey data to expand areal coverage of desert-land monitoring [
12,
15,
16]. For example, researchers found a strong relationship between percent cover of shrubs, herbaceous vegetation, and bare ground that were collected in the field and the same metrics derived from remotely sensed imagery [
17]. The cover metrics derived using the two methods were most similar for wood and shrub cover types. There was less correspondence between the two methods for herbaceous vegetation cover due to difficulties in distinguishing herbaceous vegetation in the imagery [
17].
Using remote sensing, researchers have studied conditions of vegetated landscapes often based on plant species abundance [
18] and richness [
19], plant life-form fractional cover [
18,
19,
20,
21,
22], the amount of exposed soil [
18], the amount of aboveground biomass [
20], and the presence of invasive plant species [
22]. While those studies focus on the presence or abundance of vegetation as an indicator for community conditions, other studies characterize vegetation communities based on spatial structure of the landscape, such as surface cover patterns, spatial arrangement of plants, and heterogeneity of soil resources, and use the information to evaluate the conditions of vegetation communities [
23,
24,
25,
26]. The use of spatial structure is considered particularly promising for the desert environment because of its seemingly distinguishable vegetation types and distribution patterns across a landscape.
Researchers have quantitatively characterized spatial structure or arrangement of individual plants or patches within an arid landscape by applying a variogram analysis for remotely sensed imagery [
23,
27]. The variogram is a geostatistical tool that measures spatial continuity and quantitatively characterizes spatial autocorrelation in the data [
28]. For example, researchers examined spatial characteristics of five major desert shrubland and grassland communities in the state of New Mexico using the Normalized Difference Vegetation Index (NDVI), derived from 1 m resolution multispectral images [
23]. They determined a suitable sampling interval for vegetation map** using <1 m resolution videography. Using variograms generated from images with pixel size incrementally coarsened from 50 cm to 16 m, the authors found that community-specific characteristics became obscured when the pixel size exceeded the average canopy size of the communities. Another group of researchers studied spatial patterns of vegetation communities within the Mojave Desert in California based on the parameters of variograms [
27]. By applying the canopy width as a guide to compute the variograms, the authors successfully quantified distinct spatial patterns across the vegetation community types that were dominated by different shrub species (e.g., creosote bush and bush sweepweed). They also found that each variogram parameter could indicate distinct vegetation community properties, such as the average plant canopy size and its variability in the landscape.
More recent studies have characterized the spatial distribution of NDVI across the landscape using a variogram analysis of red and near infrared bands in fine spatial resolution images. To examine bias in remotely sensed estimates of land surface variables, a research group characterized the scale of spatial variation on which the processes of landscape structure operate using the variogram [
29]. Another group unveiled spatial and temporal trends in the change rate of fractional vegetation cover using a variogram analysis in conjunction with Landsat images and vegetation index products derived from Moderate Resolution Imaging Spectroradiometer or MODIS [
30]. Variograms have also been applied to land management activities. For example, variograms were used to examine how grazing treatments affected the scale of patchiness across multiple vegetation cover types, as well as how spatial autocorrelation varied by season, and for native and non-native species [
31]. These studies clearly demonstrated the utility of the variogram applied to remotely sensed imagery to examine the spatial structures of arid vegetation communities.
The goal of this study was to determine relationships between observable remotely sensed landscape structures using a variogram method by identifying the linkage between variogram parameters and characteristics of dryland vegetation communities. Objectives were to (1) gain insight into whether or not a variogram technique is capable of detecting subtle changes in arid landscape properties, which may indicate impending irreversible habitat conversion; (2) understand trade-offs between spatial resolution of images and detectable changes within each vegetation habitat type representative of desert landscapes; and (3) determine considerations for the application of a variogram method for monitoring arid environments. We conducted this study in the context of long-term monitoring of desert lands to support sustainable utility-scale solar energy development in desert regions of Southern California using a very high spatial resolution (VHSR) multispectral aerial image.
3. Results and Discussion
3.1. Variogram Parameters for Vegetation Alliances Across Lags
Experimental variograms for the four vegetation alliances—LATR–ENFA, CHRI–GECA, LATR–AMUD, and PAFL–OLTE—with selected lags—50 cm, 1 m, 2.5 m and 5 m—are shown in
Figure 2. Experimental variograms for all lags examined are presented in
Supplemental Material (Figure S1). All variograms for the four alliances show either spherical or exponential form across most of the lags. The spherical form is characterized by the semivariance increasing at a relatively constant rate and plateauing when reaching the sill. The exponential form is characterized by the semivariance gradually increasing and forming a smooth curve when reaching the sill. Alliances dominated by shrubs and forbs, such as LATR–ENFA, CHRI–GECA, and LATR–AMDU, share a spherical form, while PAFL–OLTE has an exponential form, which likely reflects varying canopy size within the habitat (
Figure S1). In fact, the standard deviation of canopy size for PAFL–OLTE is the greatest of all alliances (2.7 m,
Table 1). The shape and characteristics of variograms become obscure as the lag increases.
The sill remains nearly constant across the lags for all alliances, which indicates that the sill is relatively insensitive to averaging of semivariances through increasing the bin widths or lag sizes. Thus, intrinsic vegetation density or cover of each alliance can be captured with any lag sizes within the lags examined (
Figure 2). The partial sill becomes increasingly smaller as the lag size increases, which suggests that variation within the landscape associated with spatial structure becomes difficult to observe when a large lag size is applied. Spatial structure becomes indistinguishable across the shrub- or forb-dominant habitats (i.e., LATR–ENFA, CHRI–GECA, and LATR–AMUD) when a lag of 1 m or greater is used. Very sparsely vegetated alliances, such as CHRI–GECA (7%) and LATR–AMDU (3%), tend to have sills lower than 13 (
Table 1 and
Figure 2), which is considerably low compared with LATR–ENFA (18%) and PAFL-OLTE (20%), where sills are higher than 27 (
Table 1 and
Figure 2).
At the 50 cm lag, the range relates to the average canopy size across the alliances. LATR–ENFA, which is characterized by small shrubs (average canopy diameter = 1.5 m) has the smallest range (0.9 m), and PAFL–OLTE that consists of trees and large shrubs (average canopy diameter = 3.9 m) has the largest range (8.8 m) (
Table 1 and
Figure 2). However, the relationship begins to dissipate at the 1 m lag; thus, the range calculated with a lag of 1 m or greater would not be indicative of the average canopy size within the alliance.
LATR–ENFA, which is characterized by moderate cover of small shrubs (
Table 1), appears to lose its characteristic variogram parameters beyond the 50 cm lag, such as the range indicative of canopy size and the partial sill indicative of spatial structure (
Figure 2). This suggests that the range and partial sill of this alliance are sensitive to smoothing, and a lag of 50 cm or smaller would be needed in order to observe spatial structure of the alliance. However, the sill indicative of vegetation cover remains constant around 20, regardless of the lag.
PAFL–OLTE containing trees and large shrubs (
Table 1) maintains its distinct variogram parameters even at the 5 m lag, including range, sill, and partial sill, despite the noticeable changes in nugget and partial sill when transposing the lag from 50 cm to 1 m (
Figure 2). This means that those variogram parameters are relatively robust to smoothing, and the habitat characteristics, such as canopy size, vegetation cover, and spatial structure, can be examined using a relatively large lag compared to other alliances.
CHRI–GECA, characterized by very sparse herbaceous cover, and LATR–AMDU, characterized by very sparse shrub cover (
Table 1), appear to mostly maintain their characteristic variogram parameters with lags up to 1 m (
Figure 2). Thus, the habitat characteristics, such as canopy size, vegetation cover, and spatial structure, could be characterized using a lag up to 1 m for these two sparsely vegetated alliances.
Based on the observations above, the lag analysis confirmed that the published relationships between variogram parameters and landscape properties [
23,
45] hold true for the desert landscape examined in this study. When using a remotely-sensed vegetation index, variogram parameters are indicative of the spatial structure of vegetation in the landscape. The variogram form indicates variability of plant canopy sizes, the sill corresponds to overall vegetation cover or density within each alliance, and the range relates to the average canopy size in a relative sense.
Averaging effects or smoothing of semivariances associated with bin widths or lag sizes has varying impacts on the variogram parameters and on alliance types. While the sill appears relatively insensitive, the partial sill is sensitive to increased bin widths because the nugget (indicating random variation in the data) becomes increasingly high when averaging over a larger bin width. These observations suggest that overall vegetation cover would be characterized relatively well, regardless of the lag across the community types, but characterizing spatial structure of the communities would require a small lag. The range appears sensitive to smoothing, and the intrinsic range quickly dissipates when a larger bin width is applied. This suggests that the average canopy size within the communities needs to be considered when selecting an appropriate lag. When considering the average canopy size, the variogram parameters of the alliance having the largest average canopy size (PAFL–OLTE) appear insensitive, while those of the alliance having the smallest average canopy size (LATR–ENFA) exhibit extreme sensitivity to smoothing. This observation clearly suggests that a small lag (~50 cm) would be required to effectively characterize intrinsic landscape structure of habitats dominated by small plants. Overall, quantifying characteristic spatial structure for each of the four alliances in a single image would require a 50 cm or smaller lag.
3.2. Variogram Parameters for Vegetation Alliances Across Image Resolutions
Experimental variograms for the four vegetation alliances—LATR–ENFA, CHRI–GECA, LATR–AMUD, and PAFL–OLTE—examined at four spatial resolutions—30 cm, 50 cm, 1 m, and 1.5 m—are shown in
Figure 3. Characteristic variogram parameters of each alliance, such as the range indicative of average canopy size and the partial sill indicative of spatial structure, become increasingly obscure as spatial resolution coarsens. This suggests that characteristic landscape properties and structure are not observable when using images having coarse spatial resolution. The sill remains constant across the spatial resolutions examined, which indicates that the vegetation cover or density can be characterized even at 1.5 m resolution.
When comparing the spatial-resolution analysis (
Figure 3) and the lag analysis (
Figure 2 and
Figure S1), the variogram parameters at a comparable unit distance (e.g., 1 m spatial resolution and 1 m lag) are not comparable. For example, the partial sill is consistently higher with the 50 cm lag than that with 50 cm resolution. The range is smaller with 50 cm lag than with 50 cm resolution, except for PAFL–OLTE. These patterns indicate that landscape characteristics, such as the average canopy size and spatial structure, which are observed in a fine resolution image with 50 cm lag would not be observable in a 50-cm resolution image. The sill, on the other hand, is comparable when comparing the variogram parameters at a comparable unit distance. The partial sill generally decreases as the resolution becomes coarser across the alliances. This indicates that detection of landscape variation associated with spatial structure decreases as spatial resolution coarsening. Particularly, at spatial resolution of 1 m or coarser, variation within the landscape becomes indistinguishable among the shrub- or forb-dominant alliances (i.e., LATR–ENFA, CHRI–GECA, and LATR–AMDU). The same pattern is observed in the lag analysis (
Figure 2).
For LATR–ENFA characterized by moderate cover (18%,
Table 1) of small shrubs and forbs, the partial sill decreases from 24.3 to 6.5 when coarsening the spatial resolution from 30 cm to 50 cm (
Figure 3). This indicates that the habitat’s intrinsic variation associated with spatial structure is no longer captured at the spatial resolution beyond 30 cm, although its overall vegetation cover is characterized by the sill, regardless of the spatial resolution. For PAFL–OLTE consisting of trees and large shrubs, the partial sill remains between 50.6 and 69.9 across the spatial resolutions, which is noticeably greater than that of other alliances (
Figure 3). This suggests that for this vegetation alliance, landscape variation associated with spatial structure can be detected with 1.5 m or finer spatial resolution. The shape of the variogram appears exponential compared to that of other alliances, indicating greater canopy size variation (standard deviation = 2.7,
Table 1).
3.3. Synthesis
It is clear that the LATR–ENFA alliance characterized by a small canopy size would require 50 cm or finer resolution to capture most of its inherent landscape characteristics. The CHRI–GECA and LATR–AMDU alliances that share similarities in landscape properties produced very similar variogram parameters. Although this could present a limitation when classifying habitat types, this result would be an indication of robust relationships between variogram parameters and landscape properties, regardless of habitat types.
The PAFL–OLTE alliance has consistent variogram parameters across a range of image resolutions. Thus, a variogram technique applied using publicly available imagery, such as NAIP imagery (60 cm~1 m), could provide useful information for detecting and monitoring internal changes in the alliance, such as changes in overall vegetation cover and canopy (e.g., degradation of shrubs, reduction of tree canopy size). Contrasting variogram forms between the tree-dominant PAFL–OLTE alliance and other shrub- and herb-dominant alliances suggests that the variogram form of PAFL–OLTE is expected to shift from an exponential form to a spherical form when canopy size diversity declines. In the context of monitoring of the Riverside East solar energy zone, a possible scenario that could be detected would be a reduction in tree abundance and size and an increase in shrub dominance in a part of the major ephemeral stream due to a change in surface water flow.
One of the most promising landscape properties that can be monitored using a variogram technique is vegetation cover within an alliance. Despite smoothing effects induced by large lags and coarse spatial resolutions, the sill appears to maintain its intrinsic height across the vegetation alliances. Thus, changes in vegetation cover can likely be detected using a range of lags and spatial resolutions. However, it should be noted that the sill cannot be directly interpreted as absolute fractional vegetation cover. In this study, the LATR–ENFA and PAFL–OLTE alliances that have similar vegetation cover (18% and 20%, respectively,
Table 1) had very different sills (
Figure 2 and
Figure 3). Therefore, it is essential to establish a baseline sill for each alliance as an indicator of fractional vegetation cover at the initiation of monitoring. Very different sills for comparable vegetation cover may be partly attributed to differences in vegetation greenness since variograms were generated from a spectral vegetation index. If this is true, the sill could indicate habitat health, as well as vegetation cover within the habitat. This hypothesis warrants further study.
When monitoring the four alliances using a single image, LATR–ENFA—having the smallest average canopy size—would determine an appropriate image resolution, which is likely the least cost-effective image type due to its fine image resolution. For sustainable long-term monitoring, this consideration is important for prioritizing monitoring efforts, such as what or how many habitats, how often, and how extensively to monitor using a variogram technique.
The comparison of variogram parameters between the lag analysis and the image-resolution analysis revealed a potential problem with determining an appropriate image resolution for applying a variogram technique to monitoring landscape structure changes. A lag up to 50 cm with a 15 cm resolution image appeared to have good potential for detecting changes in characteristic spatial structure for each of the four alliances in a single image. However, it is unlikely that a 50 cm resolution image indicates the same landscape properties for the alliances, and thus an image having a finer spatial resolution would be required, particularly for habitats characterized by small canopy sizes and sparse vegetation. This discrepancy suggests that the smoothing effect induced by increased lag sizes and that induced by degraded image resolutions do not have equal impacts on the variogram. Images having 30 cm spatial resolution or finer may be useful for monitoring internal landscape changes for dryland habitats because at this resolution, variogram parameters were relatively stable across the vegetation alliances examined (
Figure 3).
It should be noted that the results in this study may have been affected by the 13-month time span between image collection and field survey. However, because of desert vegetation characterized by slow-growth and drought-tolerant nature, the impact of the time lag is considered minimal.
4. Conclusions
By utilizing a fine-resolution spectral vegetation index, we examined the variogram technique for quantifying spatial structure of dryland environments using the four most abundant vegetation alliances in the Riverside East solar energy zone—one habitat characterized by moderate cover of small plants (LATR–ENFA), two habitats with very sparse vegetation cover of small or small to medium size canopies (LATR–AMDU and CHRI–GECA), and one habitat with moderate vegetation cover consisting of a range of plant sizes (PAFL–OLTE). Our aims were to (1) gain insights into whether or not a variogram technique has the potential to detect subtle changes in properties of desert vegetation habitats; (2) understand trade-offs between image resolution and detectable changes within each desert habitat type; and (3) determine considerations for the application of a variogram method for dryland monitoring. The findings from this study conducted in the desert region of the southwestern United States could apply to other desert environments. Key findings likely applicable for many desert regions are summarized as follows:
Variogram parameters that are generated from a spectral vegetation index using fine resolution imagery can provide information about unique, intrinsic characteristics of dryland habitats. The sill, range, form, and partial sill of the variogram relate to overall vegetation cover or density, average canopy size, canopy size variation, and spatial structure of plants within a habitat, respectively, which are consistent with published studies [
23,
27,
45].
Establishing a baseline of variogram parameters for each habitat is an important first step of monitoring because variogram parameters do not indicate absolute values of landscape properties, and they may be similar or indistinguishable between habitat types that have a similar canopy size. Comparing variogram parameters of each habitat against their baselines could indicate internal change within the habitat (e.g., vegetation cover, average canopy size, and canopy size variation). However, understanding the magnitude of changes detectable using variogram parameters requires additional studies.
For the variogram generated from a spectral vegetation index, the sill may be influenced by vegetation greenness, which could indicate the health of vegetation communities. If we can successfully detect changes in habitat conditions using remote sensing, it could greatly contribute to monitoring of extensive desert lands over time in a financially sustainable manner by substantially reducing field-based monitoring costs. To determine the effectiveness of variograms for monitoring vegetation community health, further studies are warranted.
For monitoring landscape changes, field survey and traditional image processing methods would continue to play an important role in monitoring landscape changes. The variogram technique demonstrated in this study could provide complementary knowledge on the spatial structure of the vegetated landscape, which is difficult to obtain from traditional monitoring methods.
Once the baseline variogram is established, periodic variogram analysis using publicly available imagery (e.g., 60 cm resolution NAIP imagery), commercially available fine spatial resolution imagery (e.g., 30~50 cm resolution WorldView-3 imagery), or custom very fine resolution imagery using unmanned aerial systems, or UASs, could reduce the required frequency of field data collection. Because of the accessibility of these fine resolution image options, a variogram technique is now considered a viable option for cost-effective long-term environmental monitoring. A combination of land cover information obtained from existing monitoring methods and spatial structure information extracted from the variogram analysis could improve effectiveness and timeliness of detection and monitoring of desert land landscape changes.