Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Remote Sensing of Environment”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10Linked to original sources

Marsh dieback, loss, and recovery mapped with satellite optical, airborne polarimetric radar, and field data

Landsat Thematic Mapper and Satellite Pour l'Observation de la Terre (SPOT) satellite based optical sensors, NASA Uninhabited Aerial Vehicle synthetic aperture radar (UAVSAR) polarimetric SAR (PolSAR), and field data captured the occurrence and the recovery of an undetected dieback that occurred between the summers of 2010, 2011, and 2012 in the Spartina alterniflora marshes of coastal Louisiana. Field measurements recorded the dramatic biomass decrease from 2010 to 2011 and a biomass recovery in 2012 dominated by a decrease of live biomass, and the loss of marsh as part of the dieback event. Based on an established relationship, the near-infrared/red vegetation index (VI) and site-specific measurements delineated a contiguous expanse of marsh dieback encompassing 6649.9 ha of 18,292.3 ha of S. alterniflora marshes within the study region. PolSAR data were transformed to variables used in biophysical mapping, and of this variable suite, the cross-polarization HV (horizontal send and vertical receive) backscatter was the best single indicator of marsh dieback and recovery. HV backscatter exhibited substantial and significant changes over the dieback and recovery period, tracked measured biomass changes, and significantly correlated with the live/dead biomass ratio. Within the context of regional trends, both HV and VI indicators started higher in pre-dieback marshes and exhibited substantially and statistically higher variability from year to year than that exhibited in the non-dieback marshes. That distinct difference allowed the capturing of the S. alterniflora marsh dieback and recovery; however, these changes were incorporated in a regional trend exhibiting similar but more subtle biomass composition changes.

Louisiana↗

Modeling a historical mountain pine beetle outbreak using Landsat MSS and multiple lines of evidence

Mountain pine beetles are significant forest disturbance agents, capable of inducing widespread mortality in coniferous forests in western North America. Various remote sensing approaches have assessed the impacts of beetle outbreaks over the last two decades. However, few studies have addressed the impacts of historical mountain pine beetle outbreaks, including the 1970s event that impacted Glacier National Park. The lack of spatially explicit data on this disturbance represents both a major data gap and a critical research challenge in that wildfire has removed some of the evidence from the landscape. We utilized multiple lines of evidence to model forest canopy mortality as a proxy for outbreak severity. We incorporate historical aerial and landscape photos, aerial detection survey data, a nine-year collection of satellite imagery and abiotic data. This study presents a remote sensing based framework to (1) relate measurements of canopy mortality from fine-scale aerial photography to coarse-scale multispectral imagery and (2) classify the severity of mountain pine beetle affected areas using a temporal sequence of Landsat data and other landscape variables. We sampled canopy mortality in 261 plots from aerial photos and found that insect effects on mortality were evident in changes to the Normalized Difference Vegetation Index (NDVI) over time. We tested multiple spectral indices and found that a combination of NDVI and the green band resulted in the strongest model. We report a two-step process where we utilize a generalized least squares model to account for the large-scale variability in the data and a binary regression tree to describe the small-scale variability. The final model had a root mean square error estimate of 9.8% canopy mortality, a mean absolute error of 7.6% and an R 2 of 0.82. The results demonstrate that a model of percent canopy mortality as a continuous variable can be developed to identify a gradient of mountain pine beetle severity on the landscape.

Montana↗

The evolution of mapping habitat for northern spotted owls (Strix occidentalis caurina): A comparison of photo-interpreted, Landsat-based, and lidar-based habitat maps

Wildlife habitat mapping has evolved at a rapid pace over the last few decades. Beginning with simple, often subjective, hand-drawn maps, habitat mapping now involves complex species distribution models (SDMs) using mapped predictor variables derived from remotely sensed data. For species that inhabit large geographic areas, remote sensing technology is often essential for producing range wide maps. Habitat monitoring for northern spotted owls ( Strix occidentalis caurina ), whose geographic covers about 23 million ha, is based on SDMs that use Landsat Thematic Mapper imagery to create forest vegetation data layers using gradient nearest neighbor (GNN) methods. Vegetation data layers derived from GNN are modeled relationships between forest inventory plot data, climate and topographic data, and the spectral signatures acquired by the satellite. When used as predictor variables for SDMs, there is some transference of the GNN modeling error to the final habitat map. Recent increases in the use of light detection and ranging (lidar) data, coupled with the need to produce spatially accurate and detailed forest vegetation maps have spurred interest in its use for SDMs and habitat mapping. Instead of modeling predictor variables from remotely sensed spectral data, lidar provides direct measurements of vegetation height for use in SDMs. We expect a SDM habitat map produced from directly measured predictor variables to be more accurate than one produced from modeled predictors. We used maximum entropy (Maxent) SDM modeling software to compare predictive performance and estimates of habitat area between Landsat-based and lidar-based northern spotted owl SDMs and habitat maps. We explored the differences and similarities between these maps, and to a pre-existing aerial photo-interpreted habitat map produced by local wildlife biologists. The lidar-based map had the highest predictive performance based on 10 bootstrapped replicate models (AUC = 0.809 ± 0.011), but the performance of the Landsat-based map was within acceptable limits (AUC = 0.717 ± 0.021). As is common with photo-interpreted maps, there was no accuracy assessment available for comparison. The photo-interpreted map produced the highest and lowest estimates of habitat area, depending on which habitat classes were included (nesting, roosting, and foraging habitat = 9962 ha, nesting habitat only = 6036 ha). The Landsat-based map produced an estimate of habitat area that was within this range (95% CI: 6679–9592 ha), while the lidar-based map produced an area estimate similar to what was interpreted by local wildlife biologists as nesting (i.e., high quality) habitat using aerial imagery (95% CI: 5453–7216). Confidence intervals of habitat area estimates from the SDMs based on Landsat and lidar overlapped. We concluded that both Landsat- and lidar-based SDMs produced reasonable maps and area estimates for northern spotted owl habitat within the study area. The lidar-based map was more precise and spatially similar to what local wildlife biologists considered spotted owl nesting habitat. The Landsat-based map provided a less precise spatial representation of habitat within the relatively small geographic confines of the study area, but habitat area estimates were similar to both the photo-interpreted and lidar-based maps. Photo-interpreted maps are time consuming to produce, subjective in nature, and difficult to replicate. SDMs provide a framework for efficiently producing habitat maps that can be replicated as habitat conditions change over time, provided that comparable remotely sensed data are available. When the SDM uses predictor variables extracted from lidar data, it can produce a habitat map that is both accurate and useful at large and small spatial scales. In comparison, SDMs using Landsat-based data are more appropriate for large scale analyses of amounts and general spatial patterns of habitat at regional scales.

Remote Sensing of Environment↗

Examining the utility of satellite-based wind sheltering estimates for lake hydrodynamic modeling

Satellite-based measurements of vegetation canopy structure have been in common use for the last decade but have never been used to estimate canopy's impact on wind sheltering of individual lakes. Wind sheltering is caused by slower winds in the wake of topography and shoreline obstacles (e.g. forest canopy) and influences heat loss and the flux of wind-driven mixing energy into lakes, which control lake temperatures and indirectly structure lake ecosystem processes, including carbon cycling and thermal habitat partitioning. Lakeshore wind sheltering has often been parameterized by lake surface area but such empirical relationships are only based on forested lakeshores and overlook the contributions of local land cover and terrain to wind sheltering. This study is the first to examine the utility of satellite imagery-derived broad-scale estimates of wind sheltering across a diversity of land covers. Using 30 m spatial resolution ASTER GDEM2 elevation data, the mean sheltering height, h s , being the combination of local topographic rise and canopy height above the lake surface, is calculated within 100 m-wide buffers surrounding 76,000 lakes in the U.S. state of Wisconsin. Uncertainty of GDEM2-derived h s was compared to SRTM-, high-resolution G-LiHT lidar-, and ICESat-derived estimates of h s , respective influences of land cover type and buffer width on h s are examined; and the effect of including satellite-based h s on the accuracy of a statewide lake hydrodynamic model was discussed. Though GDEM2 h s uncertainty was comparable to or better than other satellite-based measures of h s , its higher spatial resolution and broader spatial coverage allowed more lakes to be included in modeling efforts. GDEM2 was shown to offer superior utility for estimating h s compared to other satellite-derived data, but was limited by its consistent underestimation of h s , inability to detect within-buffer h s variability, and differing accuracy across land cover types. Nonetheless, considering a GDEM2 h s -derived wind sheltering potential improved the modeled lake temperature root mean square error for non-forested lakes by 0.72 °C compared to a commonly used wind sheltering model based on lake area alone. While results from this study show promise, the limitations of near-global GDEM2 data in timeliness, temporal and spatial resolution, and vertical accuracy were apparent. As hydrodynamic modeling and high-resolution topographic mapping efforts both expand, future remote sensing-derived vegetation structure data must be improved to meet wind sheltering accuracy requirements to expand our understanding of lake processes.

Remote Sensing of Environment↗

Oil detection in the coastal marshes of Louisiana using MESMA applied to band subsets of AVIRIS data

We mapped oil presence in the marshes of Barataria Bay, Louisiana following the Deepwater Horizon oil spill using Airborne Visible InfraRed Imaging Spectrometer (AVIRIS) data. Oil and non-photosynthetic vegetation (NPV) have very similar spectra, differing only in two narrow hydrocarbon absorption regions around 1700 and 2300 nm. Confusion between NPV and oil is expressed as an increase in oil fraction error with increasing NPV, as shown by Multiple Endmember Spectral Mixture Analysis (MESMA) applied to synthetic spectra generated with known endmember fractions. Significantly, the magnitude of error varied depending upon the type of NPV in the mixture. To reduce error, we used stable zone unmixing to identify a nine band subset that emphasized the hydrocarbon absorption regions, allowing for more accurate detection of oil presence using MESMA. When this band subset was applied to post-spill AVIRIS data acquired over Barataria Bay on several dates following the 2010 oil spill, accuracies ranged from 87.5% to 93.3%. Oil presence extended 10.5 m into the marsh for oiled shorelines, showing a reduced oil fraction with increasing distance from the shoreline.

Louisiana↗

A global reference database from very high resolution commercial satellite data and methodology for application to Landsat derived 30 m continuous field tree cover data

The methodology for selection, creation, and application of a global remote sensing validation dataset using high resolution commercial satellite data is presented. High resolution data are obtained for a stratified random sample of 500 primary sampling units (5 km × 5 km sample blocks), where the stratification based on Köppen climate classes is used to distribute the sample globally among biomes. The high resolution data are classified to categorical land cover maps using an analyst mediated classification workflow. Our initial application of these data is to evaluate a global 30 m Landsat-derived, continuous field tree cover product. For this application, the categorical reference classification produced at 2 m resolution is converted to percent tree cover per 30 m pixel (secondary sampling unit)for comparison to Landsat-derived estimates of tree cover. We provide example results (based on a subsample of 25 sample blocks in South America) illustrating basic analyses of agreement that can be produced from these reference data. Commercial high resolution data availability and data quality are shown to provide a viable means of validating continuous field tree cover. When completed, the reference classifications for the full sample of 500 blocks will be released for public use.

Remote Sensing of Environment↗

Characterization of shrubland ecosystem components as continuous fields in the northwest United States

Accurate and consistent estimates of shrubland ecosystem components are crucial to a better understanding of ecosystem conditions in arid and semiarid lands. An innovative approach was developed by integrating multiple sources of information to quantify shrubland components as continuous field products within the National Land Cover Database (NLCD). The approach consists of several procedures including field sample collections, high-resolution mapping of shrubland components using WorldView-2 imagery and regression tree models, Landsat 8 radiometric balancing and phenological mosaicking, medium resolution estimates of shrubland components following different climate zones using Landsat 8 phenological mosaics and regression tree models, and product validation. Fractional covers of nine shrubland components were estimated: annual herbaceous, bare ground, big sagebrush, herbaceous, litter, sagebrush, shrub, sagebrush height, and shrub height. Our study area included the footprint of six Landsat 8 scenes in the northwestern United States. Results show that most components have relatively significant correlations with validation data, have small normalized root mean square errors, and correspond well with expected ecological gradients. While some uncertainties remain with height estimates, the model formulated in this study provides a cross-validated, unbiased, and cost effective approach to quantify shrubland components at a regional scale and advances knowledge of horizontal and vertical variability of these components.

California, Idaho, Nevada, Oregon↗

Distribution of near-surface permafrost in Alaska: estimates of present and future conditions

High-latitude regions are experiencing rapid and extensive changes in ecosystem composition and function as the result of increases in average air temperature. Increasing air temperatures have led to widespread thawing and degradation of permafrost, which in turn has affected ecosystems, socioeconomics, and the carbon cycle of high latitudes. Here we overcome complex interactions among surface and subsurface conditions to map nearsurface permafrost through decision and regression tree approaches that statistically and spatially extend field observations using remotely sensed imagery, climatic data, and thematic maps of a wide range of surface and subsurface biophysical characteristics. The data fusion approach generated medium-resolution (30-m pixels) maps of near-surface (within 1 m) permafrost, active-layer thickness, and associated uncertainty estimates throughout mainland Alaska. Our calibrated models (overall test accuracy of ~85%) were used to quantify changes in permafrost distribution under varying future climate scenarios assuming no other changes in biophysical factors. Models indicate that near-surface permafrost underlies 38% of mainland Alaska and that near-surface permafrost will disappear on 16 to 24% of the landscape by the end of the 21st Century. Simulations suggest that near-surface permafrost degradation is more probable in central regions of Alaska than more northerly regions. Taken together, these results have obvious implications for potential remobilization of frozen soil carbon pools under warmer temperatures. Additionally, warmer and drier conditions may increase fire activity and severity, which may exacerbate rates of permafrost thaw and carbon remobilization relative to climate alone. The mapping of permafrost distribution across Alaska is important for land-use planning, environmental assessments, and a wide-array of geophysical studies.

Alaska↗

LiDAR based prediction of forest biomass using hierarchical models with spatially varying coefficients

Many studies and production inventory systems have shown the utility of coupling covariates derived from Light Detection and Ranging (LiDAR) data with forest variables measured on georeferenced inventory plots through regression models. The objective of this study was to propose and assess the use of a Bayesian hierarchical modeling framework that accommodates both residual spatial dependence and non-stationarity of model covariates through the introduction of spatial random effects. We explored this objective using four forest inventory datasets that are part of the North American Carbon Program, each comprising point-referenced measures of above-ground forest biomass and discrete LiDAR. For each dataset, we considered at least five regression model specifications of varying complexity. Models were assessed based on goodness of fit criteria and predictive performance using a 10-fold cross-validation procedure. Results showed that the addition of spatial random effects to the regression model intercept improved fit and predictive performance in the presence of substantial residual spatial dependence. Additionally, in some cases, allowing either some or all regression slope parameters to vary spatially, via the addition of spatial random effects, further improved model fit and predictive performance. In other instances, models showed improved fit but decreased predictive performance—indicating over-fitting and underscoring the need for cross-validation to assess predictive ability. The proposed Bayesian modeling framework provided access to pixel-level posterior predictive distributions that were useful for uncertainty mapping, diagnosing spatial extrapolation issues, revealing missing model covariates, and discovering locally significant parameters.

Colorado, Minnesota↗

Combining NLCD and MODIS to create a land cover-albedo database for the continental United States

Land surface albedo is an essential climate variable that is tightly linked to land cover, such that specific land cover classes (e.g., deciduous broadleaf forest, cropland) have characteristic albedos. Despite the normative of land-cover class specific albedos, there is considerable variability in albedo within a land cover class. The National Land Cover Database (NLCD) and the Moderate Resolution Imaging Spectroradiometer (MODIS) albedo product were combined to produce a long-term (14 years) integrated land cover-albedo database for the continental United States that can be used to examine the temporal behavior of albedo as a function of land cover. The integration identifies areas of homogeneous land cover at the nominal spatial resolution of the MODIS (MCD43A) albedo product (500 m × 500 m) from the NLCD product (30 m × 30 m), and provides an albedo data record per 500 m × 500 m pixel for 14 of the 16 NLCD land cover classes. Individual homogeneous land cover pixels have up to 605 albedo observations, and 75% of the pixels have at least 319 MODIS albedo observations (≥ 50% of the maximum possible number of observations) for the study period (2000–2013). We demonstrated the utility of the database by conducting a multivariate analysis of variance of albedo for each NLCD land cover class, showing that locational (pixel-to-pixel) and inter-annual variability were significant factors in addition to expected seasonal (intra-annual) and geographic (latitudinal) effects.

Remote Sensing of Environment↗

Quantifying soil carbon loss and uncertainty from a peatland wildfire using multi-temporal LiDAR

Peatlands are a major reservoir of global soil carbon, yet account for just 3% of global land cover. Human impacts like draining can hinder the ability of peatlands to sequester carbon and expose their soils to fire under dry conditions. Estimating soil carbon loss from peat fires can be challenging due to uncertainty about pre-fire surface elevations. This study uses multi-temporal LiDAR to obtain pre- and post-fire elevations and estimate soil carbon loss caused by the 2011 Lateral West fire in the Great Dismal Swamp National Wildlife Refuge, VA, USA. We also determine how LiDAR elevation error affects uncertainty in our carbon loss estimate by randomly perturbing the LiDAR point elevations and recalculating elevation change and carbon loss, iterating this process 1000 times. We calculated a total loss using LiDAR of 1.10 Tg C across the 25 km 2 burned area. The fire burned an average of 47 cm deep, equivalent to 44 kg C/m 2 , a value larger than the 1997 Indonesian peat fires (29 kg C/m 2 ). Carbon loss via the First-Order Fire Effects Model (FOFEM) was estimated to be 0.06 Tg C. Propagating the LiDAR elevation error to the carbon loss estimates, we calculated a standard deviation of 0.00009 Tg C, equivalent to 0.008% of total carbon loss. We conclude that LiDAR elevation error is not a significant contributor to uncertainty in soil carbon loss under severe fire conditions with substantial peat consumption. However, uncertainties may be more substantial when soil elevation loss is of a similar or smaller magnitude than the reported LiDAR error.

Virginia↗

Developing a 30-m grassland productivity estimation map for central Nebraska using 250-m MODIS and 30-m Landsat-8 observations

Accurately estimating aboveground vegetation biomass productivity is essential for local ecosystem assessment and best land management practice. Satellite-derived growing season time-integrated Normalized Difference Vegetation Index (GSN) has been used as a proxy for vegetation biomass productivity. A 250-m grassland biomass productivity map for the Greater Platte River Basin had been developed based on the relationship between Moderate Resolution Imaging Spectroradiometer (MODIS) GSN and Soil Survey Geographic (SSURGO) annual grassland productivity. However, the 250-m MODIS grassland biomass productivity map does not capture detailed ecological features (or patterns) and may result in only generalized estimation of the regional total productivity. Developing a high or moderate spatial resolution (e.g., 30-m) productivity map to better understand the regional detailed vegetation condition and ecosystem services is preferred. The 30-m Landsat data provide spatial detail for characterizing human-scale processes and have been successfully used for land cover and land change studies. The main goal of this study is to develop a 30-m grassland biomass productivity estimation map for central Nebraska, leveraging 250-m MODIS GSN and 30-m Landsat data. A rule-based piecewise regression GSN model based on MODIS and Landsat (r = 0.91) was developed, and a 30-m MODIS equivalent GSN map was generated. Finally, a 30-m grassland biomass productivity estimation map, which provides spatially detailed ecological features and conditions for central Nebraska, was produced. The resulting 30-m grassland productivity map was generally supported by the SSURGO biomass production map and will be useful for regional ecosystem study and local land management practices.

Nebraska↗

The global Landsat archive: Status, consolidation, and direction

New and previously unimaginable Landsat applications have been fostered by a policy change in 2008 that made analysis-ready Landsat data free and open access. Since 1972, Landsat has been collecting images of the Earth, with the early years of the program constrained by onboard satellite and ground systems, as well as limitations across the range of required computing, networking, and storage capabilities. Rather than robust on-satellite storage for transmission via high bandwidth downlink to a centralized storage and distribution facility as with Landsat-8, a network of receiving stations, one operated by the U.S. government, the other operated by a community of International Cooperators (ICs), were utilized. ICs paid a fee for the right to receive and distribute Landsat data and over time, more Landsat data was held outside the archive of the United State Geological Survey (USGS) than was held inside, much of it unique. Recognizing the critical value of these data, the USGS began a Landsat Global Archive Consolidation (LGAC) initiative in 2010 to bring these data into a single, universally accessible, centralized global archive, housed at the Earth Resources Observation and Science (EROS) Center in Sioux Falls, South Dakota. The primary LGAC goals are to inventory the data held by ICs, acquire the data, and ingest and apply standard ground station processing to generate an L1T analysis-ready product. As of January 1, 2015 there were 5,532,454 images in the USGS archive. LGAC has contributed approximately 3.2 million of those images, more than doubling the original USGS archive holdings. Moreover, an additional 2.3 million images have been identified to date through the LGAC initiative and are in the process of being added to the archive. The impact of LGAC is significant and, in terms of images in the collection, analogous to that of having had two additional Landsat-5 missions. As a result of LGAC, there are regions of the globe that now have markedly improved Landsat data coverage, resulting in an enhanced capacity for mapping, monitoring change, and capturing historic conditions. Although future missions can be planned and implemented, the past cannot be revisited, underscoring the value and enhanced significance of historical Landsat data and the LGAC initiative. The aim of this paper is to report the current status of the global USGS Landsat archive, document the existing and anticipated contributions of LGAC to the archive, and characterize the current acquisitions of Landsat-7 and Landsat-8. Landsat-8 is adding data to the archive at an unprecedented rate as nearly all terrestrial images are now collected. We also offer key lessons learned so far from the LGAC initiative, plus insights regarding other critical elements of the Landsat program looking forward, such as acquisition, continuity, temporal revisit, and the importance of continuing to operationalize the Landsat program.

Remote Sensing of Environment↗

Evaluating Landsat 8 evapotranspiration for water use mapping in the Colorado River Basin

Evapotranspiration (ET) mapping at the Landsat spatial resolution (100 m) is essential to fully understand water use and water availability at the field scale. Water use estimates in the Colorado River Basin (CRB), which has diverse ecosystems and complex hydro-climatic regions, will be helpful to water planners and managers. Availability of Landsat 8 images, starting in 2013, provides the opportunity to map ET in the CRB to assess spatial distribution and patterns of water use. The Operational Simplified Surface Energy Balance (SSEBop) model was used with 528 Landsat 8 images to create seamless monthly and annual ET estimates at the inherent 100 m thermal band resolution. Annual ET values were summarized by land use/land cover classes. Croplands were the largest consumer of “blue” water while shrublands consumed the most “green” water. Validation using eddy covariance (EC) flux towers and water balance approaches showed good accuracy levels with R 2 ranging from 0.74 to 0.95 and the Nash–Sutcliffe model efficiency coefficient ranging from 0.66 to 0.91. The root mean square error (and percent bias) ranged from 0.48 mm (13%) to 0.60 mm (22%) for daily (days of satellite overpass) ET and from 7.75 mm (2%) to 13.04 mm (35%) for monthly ET. The spatial and temporal distribution of ET indicates the utility of Landsat 8 for providing important information about ET dynamics across the landscape. Annual crop water use was estimated for five selected irrigation districts in the Lower CRB where annual ET per district ranged between 681 mm to 772 mm. Annual ET by crop type over the Maricopa Stanfield irrigation district ranged from a low of 384 mm for durum wheat to a high of 990 mm for alfalfa fields. A rainfall analysis over the five districts suggested that, on average, 69% of the annual ET was met by irrigation. Although the enhanced cloud-masking capability of Landsat 8 based on the cirrus band and utilization of the Fmask algorithm improved the removal of contaminated pixels, the ability to reliably estimate ET over clouded areas remains an important challenge. Overall, the performance of Landsat 8 based ET compared to available EC datasets and water balance estimates for a complex basin such as the CRB demonstrates the potential of using Landsat 8 for annual water use estimation at a national scale. Future efforts will focus on (a) use of consistent methodology across years, (b) integration of multiple sensors to maximize images used, and (c) employing cloud-computing platforms for large scale processing capabilities.

Arizona, California, Colorado, Nevada, New Mexico,↗

Perspectives on monitoring gradual change across the continuity of Landsat sensors using time-series data

There are many types of changes occurring over the Earth's landscapes that can be detected and monitored using Landsat data. Here we focus on monitoring “within-state,” gradual changes in vegetation in contrast with traditional monitoring of “abrupt” land-cover conversions. Gradual changes result from a variety of processes, such as vegetation growth and succession, damage from insects and disease, responses to shifts in climate, and other factors. Despite the prevalence of gradual changes across the landscape, they are largely ignored by the remote sensing community. Gradual changes are best characterized and monitored using time-series analysis, and with the successful launch of Landsat 8 we now have appreciable data continuity that extends the Landsat legacy across the previous 43 years. In this study, we conducted three related analyses: (1) comparison of spectral values acquired by Landsats 7 and 8, separated by eight days, to ensure compatibility for time-series evaluation; (2) tracking of multitemporal signatures for different change processes across Landsat 5, 7, and 8 sensors using anniversary-date imagery; and (3) tracking the same type of processes using all available acquisitions. In this investigation, we found that data representing natural vegetation from Landsats 5, 7, and 8 were comparable and did not indicate a need for major modification prior to use for long-term monitoring. Analyses using anniversary-date imagery can be very effective for assessing long term patterns and trends occurring across the landscape, and are especially good for providing insights regarding trends related to long-term and continuous trends of growth or decline. We found that use of all available data provided a much more comprehensive level of understanding of the trends occurring, providing information about rate, duration, and intra- and inter-annual variability that could not be readily gleaned from the anniversary date analyses. We observed that using all available clear Landsat 5–8 observations with the new Continuous Change Detection and Classification (CCDC) algorithm was very effective for illuminating vegetation trends. There are a number of potential challenges for assessing gradual changes, including atmospheric impacts, algorithm development and visualization of the changes. One of the biggest challenges for studying gradual change will be the lack of appropriate data for validating results and products.

Remote Sensing of Environment↗

Including land cover change in analysis of greenness trends using all available Landsat 5, 7, and 8 images: A case study from Guangzhou, China (2000–2014)

Remote sensing has proven a useful way of evaluating long-term trends in vegetation “greenness” through the use of vegetation indices like Normalized Differences Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI). In particular, analyses of greenness trends have been performed for large areas (continents, for example) in an attempt to understand vegetation response to climate. These studies have been most often used coarse resolution sensors like Moderate Resolution Image Spectroradiometer (MODIS) and Advanced Very High Resolution Radiometer (AVHRR). However, trends in greenness are also important at more local scales, particularly in and around cities as vegetation offers a variety of valuable ecosystem services ranging from minimizing air pollution to mitigating urban heat island effects. To explore the ability to monitor greenness trends in and around cities, this paper presents a new way for analyzing greenness trends based on all available Landsat 5, 7, and 8 images and applies it to Guangzhou, China. This method is capable of including the effects of land cover change in the evaluation of greenness trends by separating the effects of abrupt and gradual changes, and providing information on the timing of greenness trends. An assessment of the consistency of surface reflectance from Landsat 8 with past Landsat sensors indicates biases in the visible bands of Landsat 8, especially the blue band. Landsat 8 NDVI values were found to have a larger bias than the EVI values; therefore, EVI was used in the analysis of greenness trends for Guangzhou. In spite of massive amounts of development in Guangzhou from 2000 to 2014, greenness was found to increase, mostly as a result of gradual change. Comparison of the greening magnitudes estimated from the approach presented here and a Simple Linear Trend (SLT) method indicated large differences for certain time intervals as the SLT method does not include consideration for abrupt land cover changes. Overall, this analysis demonstrates the importance of considering land cover change when analyzing trends in greenness from satellite time series in areas where land cover change is common.

Guangzhou↗

Mapping changing distributions of dominant species in oil-contaminated salt marshes of Louisiana using imaging spectroscopy

The April 2010 Deepwater Horizon (DWH) oil spill was the largest coastal spill in U.S. history. Monitoring subsequent change in marsh plant community distributions is critical to assess ecosystem impacts and to establish future coastal management priorities. Strategically deployed airborne imaging spectrometers, like the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS), offer the spectral and spatial resolution needed to differentiate plant species. However, obtaining satisfactory and consistent classification accuracies over time is a major challenge, particularly in dynamic intertidal landscapes. Here, we develop and evaluate an image classification system for a time series of AVIRIS data for mapping dominant species in a heavily oiled salt marsh ecosystem. Using field-referenced image endmembers and canonical discriminant analysis (CDA), we classified 21 AVIRIS images acquired during the fall of 2010, 2011 and 2012. Classification results were evaluated using ground surveys that were conducted contemporaneously to AVIRIS collection dates. We analyzed changes in dominant species cover from 2010 to 2012 for oiled and non-oiled shorelines. CDA discriminated dominant species with a high level of accuracy (overall accuracy = 82%, kappa = 0.78) and consistency over three imaging dates (overall 2010 = 82%, overall 2011 = 82%, overall 2012 = 88%). Marshes dominated by Spartina alterniflora were the most spatially abundant in shoreline zones (≤ 28 m from shore) for all three dates (2010 = 79%, 2011 = 61%, 2012 = 63%), followed by Juncus roemerianus (2010 = 11%, 2011 = 19%, 2012 = 17%) and Distichlis spicata (2010 = 4%, 2011 = 10%, 2012 = 7%). Marshes that were heavily contaminated with oil exhibited variable responses from 2010 to 2012. Marsh vegetation classes converted to a subtidal, open water class along oiled and non-oiled shorelines that were similarly situated in the landscape. However, marsh loss along oil-contaminated shorelines doubled that of non-oiled shorelines. Only S. alterniflora dominated marshes were extensively degraded, losing 15% (354,604 m 2 ) cover in oiled shoreline zones, suggesting that S. alterniflora marshes may be more vulnerable to shoreline erosion following hydrocarbon stress, due to their landscape position.

Louisiana↗

Multi-index time series monitoring of drought and fire effects on desert grasslands

The Western United States is expected to undergo both extended periods of drought and longer wildfire seasons under forecasted global climate change and it is important to understand how these disturbances will interact and affect recovery and composition of plant communities in the future. In this research paper we describe the temporal response of grassland communities to drought and fire in southern Arizona, where land managers are using repeated, prescribed fire as a habitat restoration tool. Using a 25-year atlas of fire locations, we paired sites with multiple fires to unburned control areas and compare satellite and field-based estimates of vegetation cover over time. Two hundred and fifty Landsat TM images, dating from 1985–2011, were used to derive estimates of Total Vegetation Fractional Cover (TVFC) of live and senescent grass using the Soil-Adjusted Total Vegetation Index (SATVI) and post-fire vegetation greenness using the Normalized Difference Vegetation Index (NDVI). We also implemented a Greenness to Cover Index that is the difference of time-standardized SATVI-TVFC and NDVI values at a given time and location to identify post-fire shifts in native, non-native, and annual plant cover. The results highlight anomalous greening and browning during drought periods related to amounts of annual and non-native plant cover present. Results suggest that aggressive application of prescribed fire may encourage spread of non-native perennial grasses and annual plants, particularly during droughts.

Remote Sensing of Environment↗