Landsat Time Series Analysis of Vegetation Changes in Solar Energy Development Areas of the Lower Colorado Desert , Southern California

Land cover change from renewable energy development in southern California is receiving increasing attention due to potential impacts on protected area conservation, endangered species, and greenhouse gas emissions. This study was designed to quantify and map, for the first time, variations desert vegetation canopy density and related growth rates using 30 consecutive years of Landsat satellite image data across the Lower Colorado Desert. The time-series for mean normalized difference vegetation index (NDVI) values sampled from each of the three major land cover types (shrubland, barren sand dune, and developed urban) showed no significant positive or negative trend in vegetation canopy density. Three periods of significant decrease in NDVI were detected during the drought periods of 1989-1990, 2002-2003, and 2013-2015, indicating that annual precipitation has been the main controller of shrubland canopy growth and green cover. Shrubland canopy cover has been relatively stable in renewable energy development zones since the mid-2000s. NDVI change in the period after nearly all southern California solar energy developments were initiated (post-2010) could be attributed largely to topographic water flow pathways through canyons and desert washes, both in and around all solar energy development zones.


Introduction
Damage to native shrubland vegetation in southern California resulting from solar energy development projects has been speculated to cause negative impacts on biodiversity conservation and greenhouse gas emissions from the disturbance or destruction of desert plant biomass and surface soils [1]- [3].Renewable energy development projects sited near protected areas (e.g., state and national parks) have been suspected of disrupting dispersion capabilities of native species and facilitating invasions of non-native species through habitat loss, fragmentation, and isolation [2].These same assessment reports [2] [3] have called for, but have not included, rigorous statistical analyses that could quantify both internal and external threats to protected southern California desert ecosystems, as required for assessing vulnerability of landscape sustainability.
The objective of this remote sensing data analysis for southern California deserts was to quantify and characterize, for the first time, variations desert vegetation canopy density and related growth rates using 30 consecutive years of Landsat satellite image data across the Lower Colorado Desert region (Figure 1).The normalized difference vegetation index (NDVI) from Landsat has been closely correlated with percent cover measurements of green vegetation canopies in arid eco-systems, as reported in numerous published studies [4]- [10].More specifically, reference [11] found that NDVI from Landsat was tightly correlated (r = 0.88; p < 0.05) with total percent cover of live vegetation in a semi-arid sagebrush ecosystem in south-central Utah.Reference [12] reported similar results for NDVI correlations with fractional cover in creosote-dominated shrubland sites of New Mexico.
The Desert Renewable Energy Conservation Plan [13] (DRECP) has become a major component of the State of California Renewables Portfolio Standard.The DRECP is intended to provide effective environmental protections and natural resource conservation of desert ecosystems as appropriate development of renewable energy projects advances in southern California.The DRECP covers parts of seven California counties: Imperial, Inyo, Kern, Los Angeles, Riverside, San Bernardino, and San Diego.Approximately 91,000 km 2 of federal and nonfederal California desert land are part of the DRECP Area.By the end of 2016, the Final DRECP Environment Impact Statement (EIS) will be completed for public lands under agency jurisdictions within the DRECP area, including numerous zones designated as Development Focus Areas (DFAs), as shown in Figure 1.

Study Area
The area of interest for this Landsat data study of vegetation cover growth within the southern DRECP region was the Lower Colorado Desert of California (Figure 1).This area is bounded on the west by the Laguna, Santa Rosa, and San Jacinto mountain ranges, on the east by the California-Arizona state line, on the north by the gradual transition to the Mojave Desert, and on the south by the California-Mexico border [14].Low annual rainfall (50 -300 mm) and high temperatures (reaching 45˚C in the summer) make this area one of the most arid in North America.Vegetation communities in the Colorado Desert have been classified into seven basic types: creosote bush scrub, cactus scrub, saltbush scrub, alkali sink, microphyll woodland, psammaophytic scrub, and palm oasis [15].

Methods
For this study, near cloud-free imagery from the Landsat sensor was selected from the United States Geological Survey (USGS) Earth Explorer data portal (available online at earthexplorer.usgs.gov)for every year from 1985 to 2015.Landsat scenes from path/row 39/37 were acquired between April 1 and May 30 each year, around the peak flowering period in the Lower Colorado Desert growing season [15].All images used in this study were geometrically registered (UTM Zone 10) using terrain correction algorithms (Level 1T) applied by the USGS EROS Data Center.
For the Landsat 4 -5 Thematic Mapper (TM) images acquired between 1985 and 2011, 30-m resolution surface reflectance data were generated from the Landsat Ecosystem Disturbance Adaptive Processing System [16].Moderate Resolution Imaging Spectroradiometer (MODIS) atmospheric correction routines were applied to Level-1 TM data products.Water vapor, ozone, geopotential height, aerosol optical thickness, and digital elevation are input with Landsat data to the Second Simulation of a Satellite Signal in the Solar Spectrum (6S) radiative transfer models to generate top of atmosphere (TOA) reflectance, surface reflectance, brightness temperature, and masks for clouds, cloud shadows, adjacent clouds, land, snow, ice, and water.Landsat 8 (after 2012) surface reflectance products were generated from the L8SR algorithm, a method that uses the scene center for the sun angle calculation and then hard-codes the view zenith angle to 0. The solar zenith and view zenith angles are used for calculations as part of the atmospheric correction.
NDVI was computed for this study using Landsat bands as: where Red is the reflectance from 0.63 to 0.69 μm and NIR is the near-infrared reflectance from 0.77 to 0.90 μm [17].NDVI was scaled from 0 to 10,000, with low values of (near 0) indicating barren land cover and high values of NDVI (above 7000) indicating dense canopy vegetation cover [18].Land cover types for the study area were determined based on the National Agricultural Statistics Service (NASS), California Cropland Data Layer (CDL) from 2014 (available online at nassgeodata.gmu.edu/CropScape).The CDL is a raster, geo-referenced land cover data layer with a ground resolution of 30 m produced using satellite imagery from the Landsat sensor collected during the summer season.Elevation, slope, and aspect were determined at 30 m spatial resolution from the USGS National Elevation Dataset (NED).
Calculations of statistical significance of trends over time were computed to test the hypothesis of chance deviations from a flat trend, following the same analysis approach applied by reference [19].A 95% confidence interval for a yearly sample mean difference was reported as greater than (plus or minus) 2 standard errors (SE) of the mean (at p < 0.05) [20].Within each of the three major land cover types of the region (shrublands, barren sand dunes, and developed urban areas), 2000 point locations were randomly selected for statistical summaries and comparisons of yearly mean NDVI values.

Results and Discussion
Resulting plots of that the complete 30-year time-series for mean NDVI sampled from each of the three major land cover types showed no significant positive or negative trends in vegetation canopy density since 1985 (Figure 2).There were, nonetheless, three notable periods of significant (p < 0.05) decrease in the mean annual NDVI from the 30-year NDVI averages (by vegetation type), namely during the droughts of 1989-1990, 2002-2003, and 2013-2015.Across all 30 years, 2 SE of the mean yearly NDVI sampled for shrubland cover varied between 20 and 40 NDVI units, whereas the three largest negative departures from the 30-year shrubland NDVI mean were between 300 and 500 NDVI units during these extreme drought periods in southern California (as also documented by reference [21]).
Comparisons of the three different cover types showed that developed urban and industrial locations had vegetation green cover at levels consistently 1000 NDVI higher than desert shrubland cover locations, which receive no watering by residential or commercial irrigation systems (Figure 2).Since the mid-1980s, mean yearly NDVI levels for shrubland cover within DFAs have been 200 and 500 NDVI units below shrubland cover outside of the DFAs, which indicates that the most productive shrubland cover in the Lower Colorado Desert has not yet been included within DFAs for solar energy development.As expected, barren sand dune covered locations showed the lowest mean NDVI levels and the lowest within-year variations across all years.
The spatial pattern of change in NDVI (Figure 3), starting from the period before solar energy developments  were initiated in the region to the present, showed that DFAs of East Riverside County were largely unchanged in vegetation canopy cover, or saw increases of between 200 and 400 NDVI units following river channels and desert washes, such as Ford Dry Lake, Lower Palen Valley, McCoy Wash, and Cinnabar Wash.The most conspicuous patterns of recent NDVI increase followed these well-delineated hydrologic drainage basins, all at lower than 200 meters elevation [22].The only extensive areas near DFAs where NDVI recently decreased by between 200 and 400 units were located on steep slopes above 400 meters elevation in the McCoy and Palen Mountains where solar energy developments are not currently planned.

Conclusions
It can be concluded that inter-annual variations in precipitation accounted for nearly all the periodic changes in NDVI in shrubland communities observed since 1985 across the Lower Colorado Desert.There have been statistically significant decreases in shrubland cover to a stable baseline level of about 1000 NDVI units during each extreme drought period recorded over the past 30 years, followed by recovery to levels between 1500 and 2000 NDVI units during more favorable rainfall years.Hernandez et al. [2] reported that nearly 50% of utility scale solar energy development (greater than 1 megawatt capacity) in California is being sited in shrubland and scrubland cover types, potentially leading to increased isolation of endangered species and non-native species invasions within these ecosystems.However, results from analysis of Landsat NDVI 30-year trends and spatial patterns of change over the past decade do not support contentions that detectable disturbance of shrubland vegetation cover has yet resulted in adverse ecological impacts to communities of the Lower Colorado Desert region, either within DFAs of East Riverside County or in shrublands outside of the DFAs.Instead, shrubland canopy cover has been relatively stable since the mid-2000s.Subtle changes in NDVI in the period after nearly all southern California solar energy developments were initiated (post-2010, [2]) could be attributed largely to topographic water flow pathways through canyons and desert washes, both in and around the solar energy DFAs.

Figure 1 .
Figure 1.Map of the lower Colorado Desert (brown line) and DRECP (red line) study area, with development focus areas (DFAs) in other shaded polygon outlines.Smaller map in the right panel shows vegetation cover types for the entire DRECP area, with shrublands in light green, barren sand dunes in tan, developed urban areas in grey, croplands in magenta, and open water in blue.

Figure 2 .
Figure 2. Plot of the 30-year trend in Landsat NDVI for lower Colorado desert vegetation cover types.Error bars represent 2 SE of the yearly mean NDVI computed for 2000 randomly selected points within each vegetation cover type.

Figure 3 .
Figure 3. Map image of the change in Landsat NDVI between 2003 (prior to solar energy developments) and the present (2015) for East Riverside DFAs outlined in red.Positive changes indicate an increase in NDVI since 2003.Drainage basin names are included from Seaber et al. [22].