Search USGS⌕ Search

SEARCH · Search USGS

Results for “Algorithms”

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 793 records · Page 44Linked to original sources

Estimating abiotic thresholds for sagebrush condition class in the western United States

Sagebrush ecosystems of the western United States can transition from extended periods of relatively stable conditions to rapid ecological change if acute disturbances occur. Areas dominated by native sagebrush can transition from species-rich native systems to altered states where non-native annual grasses dominate, if resistance to annual grasses is low. The non-native annual grasses provide relatively little value to wildlife, livestock, and humans and function as fuel that increases fire frequency. The more land area covered by annual grasses, the higher the potential for fire, thus reducing the potential for native vegetation to reestablish, even when applying restoration treatments. Mapping areas of stability and areas of change using machine-learning algorithms allows both the identification of dominant abiotic variables that drive ecosystem dynamics and the variables’ important thresholds. We develop a decision-tree model with rulesets that estimate three classes of sagebrush condition (i.e., sagebrush recovery, tipping point [ecosystem degradation], and stable). We find rulesets that primarily drive development of the sagebrush recovery class indicate areas of midelevations (1 602 m), warm 30-yr July temperature maximums (tmax) (30.62°C), and 30-yr March precipitation (ppt) averages equal to 26.26 mm, about 10% of the 30-yr annual ppt values. Tipping point and stable classes occur at elevations that are lower (1 505 m) and higher (1 939 m), respectively, more mesic during March and annually, and experience lower 30-yr July tmax averages. These defined variable averages can be used to understand current dynamics of sagebrush condition and to predict where future transitions may occur under novel conditions.

Arizona, California, Colorado, Idaho, Montana, Neb↗

A comparative analysis of the Global Land Cover 2000 and MODIS land cover data sets

Accurate and up-to-date global land cover data sets are necessary for various global change research studies including climate change, biodiversity conservation, ecosystem assessment, and environmental modeling. In recent years, substantial advancement has been achieved in generating such data products. Yet, we are far from producing geospatially consistent high-quality data at an operational level. We compared the recently available Global Land Cover 2000 (GLC-2000) and MODerate resolution Imaging Spectrometer (MODIS) global land cover data to evaluate the similarities and differences in methodologies and results, and to identify areas of spatial agreement and disagreement. These two global land cover data sets were prepared using different data sources, classification systems, and methodologies, but using the same spatial resolution (i.e., 1 km) satellite data. Our analysis shows a general agreement at the class aggregate level except for savannas/shrublands, and wetlands. The disagreement, however, increases when comparing detailed land cover classes. Similarly, percent agreement between the two data sets was found to be highly variable among biomes. The identified areas of spatial agreement and disagreement will be useful for both data producers and users. Data producers may use the areas of spatial agreement for training area selection and pay special attention to areas of disagreement for further improvement in future land cover characterization and mapping. Users can conveniently use the findings in the areas of agreement, whereas users might need to verify the informaiton in the areas of disagreement with the help of secondary information. Learning from past experience and building on the existing infrastructure (e.g., regional networks), further research is necessary to (1) reduce ambiguity in land cover definitions, (2) increase availability of improved spatial, spectral, radiometric, and geometric resolution satellite data, and (3) develop advanced classification algorithms.

Remote Sensing of Environment↗

Mapping an invasive plant, Phragmites australis, in coastal wetlands using the EO-1 Hyperion hyperspectral sensor

Mapping tools are needed to document the location and extent of Phragmites australis, a tall grass that invades coastal marshes throughout North America, displacing native plant species and degrading wetland habitat. Mapping Phragmites is particularly challenging in the freshwater Great Lakes coastal wetlands due to dynamic lake levels and vegetation diversity. We tested the applicability of Hyperion hyperspectral satellite imagery for mapping Phragmites in wetlands of the west coast of Green Bay in Wisconsin, U.S.A. A reference spectrum created using Hyperion data from several pure Phragmites stands within the image was used with a Spectral Correlation Mapper (SCM) algorithm to create a raster map with values ranging from 0 to 1, where 0 represented the greatest similarity between the reference spectrum and the image spectrum and 1 the least similarity. The final two-class thematic classification predicted monodominant Phragmites covering 3.4% of the study area. Most of this was concentrated in long linear features parallel to the Green Bay shoreline, particularly in areas that had been under water only six years earlier when lake levels were 66??cm higher. An error matrix using spring 2005 field validation points (n = 129) showed good overall accuracy-81.4%. The small size and linear arrangement of Phragmites stands was less than optimal relative to the sensor resolution, and Hyperion's 30??m resolution captured few if any pure pixels. Contemporary Phragmites maps prepared with Hyperion imagery would provide wetland managers with a tool that they currently lack, which could aid attempts to stem the spread of this invasive species. ?? 2006 Elsevier Inc. All rights reserved.

Remote Sensing of Environment↗

Using LiDAR and quickbird data to model plant production and quantify uncertainties associated with wetland detection and land cover generalizations

Spatiotemporal data from satellite remote sensing and surface meteorology networks have made it possible to continuously monitor global plant production, and to identify global trends associated with land cover/use and climate change. Gross primary production (GPP) and net primary production (NPP) are routinely derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard satellites Terra and Aqua, and estimates generally agree with independent measurements at validation sites across the globe. However, the accuracy of GPP and NPP estimates in some regions may be limited by the quality of model input variables and heterogeneity at fine spatial scales. We developed new methods for deriving model inputs (i.e., land cover, leaf area, and photosynthetically active radiation absorbed by plant canopies) from airborne laser altimetry (LiDAR) and Quickbird multispectral data at resolutions ranging from about 30??m to 1??km. In addition, LiDAR-derived biomass was used as a means for computing carbon-use efficiency. Spatial variables were used with temporal data from ground-based monitoring stations to compute a six-year GPP and NPP time series for a 3600??ha study site in the Great Lakes region of North America. Model results compared favorably with independent observations from a 400??m flux tower and a process-based ecosystem model (BIOME-BGC), but only after removing vapor pressure deficit as a constraint on photosynthesis from the MODIS global algorithm. Fine-resolution inputs captured more of the spatial variability, but estimates were similar to coarse-resolution data when integrated across the entire landscape. Failure to account for wetlands had little impact on landscape-scale estimates, because vegetation structure, composition, and conversion efficiencies were similar to upland plant communities. Plant productivity estimates were noticeably improved using LiDAR-derived variables, while uncertainties associated with land cover generalizations and wetlands in this largely forested landscape were considered less important.

Remote Sensing of Environment↗

Modeling fire severity in black spruce stands in the Alaskan boreal forest using spectral and non-spectral geospatial data

Biomass burning in the Alaskan interior is already a major disturbance and source of carbon emissions, and is likely to increase in response to the warming and drying predicted for the future climate. In addition to quantifying changes to the spatial and temporal patterns of burned areas, observing variations in severity is the key to studying the impact of changes to the fire regime on carbon cycling, energy budgets, and post-fire succession. Remote sensing indices of fire severity have not consistently been well-correlated with in situ observations of important severity characteristics in Alaskan black spruce stands, including depth of burning of the surface organic layer. The incorporation of ancillary data such as in situ observations and GIS layers with spectral data from Landsat TM/ETM+ greatly improved efforts to map the reduction of the organic layer in burned black spruce stands. Using a regression tree approach, the R2 of the organic layer depth reduction models was 0.60 and 0.55 (pb0.01) for relative and absolute depth reduction, respectively. All of the independent variables used by the regression tree to estimate burn depth can be obtained independently of field observations. Implementation of a gradient boosting algorithm improved the R2 to 0.80 and 0.79 (pb0.01) for absolute and relative organic layer depth reduction, respectively. Independent variables used in the regression tree model of burn depth included topographic position, remote sensing indices related to soil and vegetation characteristics, timing of the fire event, and meteorological data. Post-fire organic layer depth characteristics are determined for a large (N200,000 ha) fire to identify areas that are potentially vulnerable to a shift in post-fire succession. This application showed that 12% of this fire event experienced fire severe enough to support a change in post-fire succession. We conclude that non-parametric models and ancillary data are useful in the modeling of the surface organic layer fire depth. Because quantitative differences in post-fire surface characteristics do not directly influence spectral properties, these modeling techniques provide better information than the use of remote sensing data alone.

Remote Sensing of Environment↗

Spectral assessment of new ASTER SWIR surface reflectance data products for spectroscopic mapping of rocks and minerals

ASTER reflectance spectra from Cuprite, Nevada, and Mountain Pass, California, were compared to spectra of field samples and to ASTER-resampled AVIRIS reflectance data to determine spectral accuracy and spectroscopic mapping potential of two new ASTER SWIR reflectance datasets: RefL1b and AST_07XT. RefL1b is a new reflectance dataset produced for this study using ASTER Level 1B data, crosstalk correction, radiance correction factors, and concurrently acquired level 2 MODIS water vapor data. The AST_07XT data product, available from EDC and ERSDAC, incorporates crosstalk correction and non-concurrently acquired MODIS water vapor data for atmospheric correction. Spectral accuracy was determined using difference values which were compiled from ASTER band 5/6 and 9/8 ratios of AST_07XT or RefL1b data subtracted from similar ratios calculated for field sample and AVIRIS reflectance data. In addition, Spectral Analyst, a statistical program that utilizes a Spectral Feature Fitting algorithm, was used to quantitatively assess spectral accuracy of AST_07XT and RefL1b data.Spectral Analyst matched more minerals correctly and had higher scores for the RefL1b data than for AST_07XT data. The radiance correction factors used in the RefL1b data corrected a low band 5 reflectance anomaly observed in the AST_07XT and AST_07 data but also produced anomalously high band 5 reflectance in RefL1b spectra with strong band 5 absorption for minerals, such as alunite. Thus, the band 5 anomaly seen in the RefL1b data cannot be corrected using additional gain adjustments. In addition, the use of concurrent MODIS water vapor data in the atmospheric correction of the RefL1b data produced datasets that had lower band 9 reflectance anomalies than the AST_07XT data. Although assessment of spectral data suggests that RefL1b data are more consistent and spectrally more correct than AST_07XT data, the Spectral Analyst results indicate that spectral discrimination between some minerals, such as alunite and kaolinite, are still not possible unless additional spectral calibration using site specific spectral data are performed. ?? 2010.

Remote Sensing of Environment↗

A simple and effective method for filling gaps in Landsat ETM+ SLC-off images

The scan-line corrector (SLC) of the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) sensor failed in 2003, resulting in about 22% of the pixels per scene not being scanned. The SLC failure has seriously limited the scientific applications of ETM+ data. While there have been a number of methods developed to fill in the data gaps, each method has shortcomings, especially for heterogeneous landscapes. Based on the assumption that the same-class neighboring pixels around the un-scanned pixels have similar spectral characteristics, and that these neighboring and un-scanned pixels exhibit similar patterns of spectral differences between dates, we developed a simple and effective method to interpolate the values of the pixels within the gaps. We refer to this method as the Neighborhood Similar Pixel Interpolator (NSPI). Simulated and actual SLC-off ETM+ images were used to assess the performance of the NSPI. Results indicate that NSPI can restore the value of un-scanned pixels very accurately, and that it works especially well in heterogeneous regions. In addition, it can work well even if there is a relatively long time interval or significant spectral changes between the input and target image. The filled images appear reasonably spatially continuous without obvious striping patterns. Supervised classification using the maximum likelihood algorithm was done on both gap-filled simulated SLC-off data and the original "gap free" data set, and it was found that classification results, including accuracies, were very comparable. This indicates that gap-filled products generated by NSPI will have relevance to the user community for various land cover applications. In addition, the simple principle and high computational efficiency of NSPI will enable processing large volumes of SLC-off ETM+ data.

Remote Sensing of Environment↗

Modeling the height of young forests regenerating from recent disturbances in Mississippi using Landsat and ICESat data

Many forestry and earth science applications require spatially detailed forest height data sets. Among the various remote sensing technologies, lidar offers the most potential for obtaining reliable height measurement. However, existing and planned spaceborne lidar systems do not have the capability to produce spatially contiguous, fine resolution forest height maps over large areas. This paper describes a Landsat–lidar fusion approach for modeling the height of young forests by integrating historical Landsat observations with lidar data acquired by the Geoscience Laser Altimeter System (GLAS) instrument onboard the Ice, Cloud, and land Elevation (ICESat) satellite. In this approach, “young” forests refer to forests reestablished following recent disturbances mapped using Landsat time-series stacks (LTSS) and a vegetation change tracker (VCT) algorithm. The GLAS lidar data is used to retrieve forest height at sample locations represented by the footprints of the lidar data. These samples are used to establish relationships between lidar-based forest height measurements and LTSS–VCT disturbance products. The height of “young” forest is then mapped based on the derived relationships and the LTSS–VCT disturbance products. This approach was developed and tested over the state of Mississippi. Of the various models evaluated, a regression tree model predicting forest height from age since disturbance and three cumulative indices produced by the LTSS–VCT method yielded the lowest cross validation error. The R 2 and root mean square difference (RMSD) between predicted and GLAS-based height measurements were 0.91 and 1.97 m, respectively. Predictions of this model had much higher errors than indicated by cross validation analysis when evaluated using field plot data collected through the Forest Inventory and Analysis Program of USDA Forest Service. Much of these errors were due to a lack of separation between stand clearing and non-stand clearing disturbances in current LTSS–VCT products and difficulty in deriving reliable forest height measurements using GLAS samples when terrain relief was present within their footprints. In addition, a systematic underestimation of about 5 m by the developed model was also observed, half of which could be explained by forest growth that occurred between field measurement year and model target year. The remaining difference suggests that tree height measurements derived using waveform lidar data could be significantly underestimated, especially for young pine forests. Options for improving the height modeling approach developed in this study were discussed.

Remote Sensing of Environment↗

Mapping ground surface deformation using temporarily coherent point SAR interferometry: Application to Los Angeles Basin

Multi-temporal interferometric synthetic aperture radar (InSAR) is an effective tool to detect long-term seismotectonic motions by reducing the atmospheric artifacts, thereby providing more precise deformation signal. The commonly used approaches such as persistent scatterer InSAR (PSInSAR) and small baseline subset (SBAS) algorithms need to resolve the phase ambiguities in interferogram stacks either by searching a predefined solution space or by sparse phase unwrapping methods; however the efficiency and the success of phase unwrapping cannot be guaranteed. We present here an alternative approach – temporarily coherent point (TCP) InSAR (TCPInSAR) – to estimate the long term deformation rate without the need of phase unwrapping. The proposed approach has a series of innovations including TCP identification, TCP network and TCP least squares estimator. We apply the proposed method to the Los Angeles Basin in southern California where structurally active faults are believed capable of generating damaging earthquakes. The analysis is based on 55 interferograms from 32 ERS-1/2 images acquired during Oct. 1995 and Dec. 2000. To evaluate the performance of TCPInSAR on a small set of observations, a test with half of interferometric pairs is also performed. The retrieved TCPInSAR measurements have been validated by a comparison with GPS observations from Southern California Integrated GPS Network. Our result presents a similar deformation pattern as shown in past InSAR studies but with a smaller average standard deviation (4.6 mm) compared with GPS observations, indicating that TCPInSAR is a promising alternative for efficiently mapping ground deformation even from a relatively smaller set of interferograms.

Remote Sensing of Environment↗

A comprehensive change detection method for updating the National Land Cover Database to circa 2011

The importance of characterizing, quantifying, and monitoring land cover, land use, and their changes has been widely recognized by global and environmental change studies. Since the early 1990s, three U.S. National Land Cover Database (NLCD) products (circa 1992, 2001, and 2006) have been released as free downloads for users. The NLCD 2006 also provides land cover change products between 2001 and 2006. To continue providing updated national land cover and change datasets, a new initiative in developing NLCD 2011 is currently underway. We present a new Comprehensive Change Detection Method (CCDM) designed as a key component for the development of NLCD 2011 and the research results from two exemplar studies. The CCDM integrates spectral-based change detection algorithms including a Multi-Index Integrated Change Analysis (MIICA) model and a novel change model called Zone, which extracts change information from two Landsat image pairs. The MIICA model is the core module of the change detection strategy and uses four spectral indices (CV, RCVMAX, dNBR, and dNDVI) to obtain the changes that occurred between two image dates. The CCDM also includes a knowledge-based system, which uses critical information on historical and current land cover conditions and trends and the likelihood of land cover change, to combine the changes from MIICA and Zone. For NLCD 2011, the improved and enhanced change products obtained from the CCDM provide critical information on location, magnitude, and direction of potential change areas and serve as a basis for further characterizing land cover changes for the nation. An accuracy assessment from the two study areas show 100% agreement between CCDM mapped no-change class with reference dataset, and 18% and 82% disagreement for the change class for WRS path/row p22r39 and p33r33, respectively. The strength of the CCDM is that the method is simple, easy to operate, widely applicable, and capable of capturing a variety of natural and anthropogenic disturbances potentially associated with land cover changes on different landscapes.

Remote Sensing of Environment↗

A comprehensive evaluation of two MODIS evapotranspiration products over the conterminous United States: using point and gridded FLUXNET and water balance ET

Remote sensing datasets are increasingly being used to provide spatially explicit large scale evapotranspiration (ET) estimates. Extensive evaluation of such large scale estimates is necessary before they can be used in various applications. In this study, two monthly MODIS 1 km ET products, MODIS global ET (MOD16) and Operational Simplified Surface Energy Balance (SSEBop) ET, are validated over the conterminous United States at both point and basin scales. Point scale validation was performed using eddy covariance FLUXNET ET (FLET) data (2001–2007) aggregated by year, land cover, elevation and climate zone. Basin scale validation was performed using annual gridded FLUXNET ET (GFET) and annual basin water balance ET (WBET) data aggregated by various hydrologic unit code (HUC) levels. Point scale validation using monthly data aggregated by years revealed that the MOD16 ET and SSEBop ET products showed overall comparable annual accuracies. For most land cover types, both ET products showed comparable results. However, SSEBop showed higher performance for Grassland and Forest classes; MOD16 showed improved performance in the Woody Savanna class. Accuracy of both the ET products was also found to be comparable over different climate zones. However, SSEBop data showed higher skill score across the climate zones covering the western United States. Validation results at different HUC levels over 2000–2011 using GFET as a reference indicate higher accuracies for MOD16 ET data. MOD16, SSEBop and GFET data were validated against WBET (2000–2009), and results indicate that both MOD16 and SSEBop ET matched the accuracies of the global GFET dataset at different HUC levels. Our results indicate that both MODIS ET products effectively reproduced basin scale ET response (up to 25% uncertainty) compared to CONUS-wide point-based ET response (up to 50–60% uncertainty) illustrating the reliability of MODIS ET products for basin-scale ET estimation. Results from this research would guide the additional parameter refinement required for the MOD16 and SSEBop algorithms in order to further improve their accuracy and performance for agro-hydrologic applications.

Remote Sensing of Environment↗

Monitoring conterminous United States (CONUS) land cover change with Web-Enabled Landsat Data (WELD)

Forest cover loss and bare ground gain from 2006 to 2010 for the conterminous United States (CONUS) were quantified at a 30 m spatial resolution using Web-Enabled Landsat Data available from the USGS Center for Earth Resources Observation and Science (EROS) (http://landsat.usgs.gov/WELD.php). The approach related multi-temporal WELD metrics and expert-derived training data for forest cover loss and bare ground gain through a decision tree classification algorithm. Forest cover loss was reported at state and ecoregional scales, and the identification of core forests' absent of change was made and verified using LiDAR data from the GLAS (Geoscience Laser Altimetry System) instrument. Bare ground gain correlated with population change for large metropolitan statistical areas (MSAs) outside of desert or semi-desert environments. GoogleEarth™ time-series images were used to validate the products. Mapped forest cover loss totaled 53,084 km2 and was found to be depicted conservatively, with a user's accuracy of 78% and a producer's accuracy of 68%. Excluding errors of adjacency, user's and producer's accuracies rose to 93% and 89%, respectively. Mapped bare ground gain equaled 5974 km2 and nearly matched the estimated area from the reference (GoogleEarth™) classification; however, user's (42%) and producer's (49%) accuracies were much less than those of the forest cover loss product. Excluding errors of adjacency, user's and producer's accuracies rose to 62% and 75%, respectively. Compared to recent 2001–2006 USGS National Land Cover Database validation data for forest loss (82% and 30% for respective user's and producer's accuracies) and urban gain (72% and 18% for respective user's and producer's accuracies), results using a single CONUS-scale model with WELD data are promising and point to the potential for national-scale operational mapping of key land cover transitions. However, validation results highlighted limitations, some of which can be addressed by improving training data, creating a more robust image feature space, adding contemporaneous Landsat 5 data to the inputs, and modifying definition sets to account for differences in temporal and spatial observational scales. The presented land cover extent and change data are available via the official WELD website (ftp://weldftp.cr.usgs.gov/CONUS_5Y_LandCover/ftp://weldftp.cr.usgs.gov/CONUS_5Y_LandCover/).

Remote Sensing of Environment↗

Rapid dispersal of saltcedar (Tamarix spp.) biocontrol beetles (Diorhabda carinulata) on a desert river detected by phenocams, MODIS imagery and ground observations

We measured the rate of dispersal of saltcedar leaf beetles ( Diorhabda carinulata ), a defoliating insect released on western rivers to control saltcedar shrubs ( Tamarix spp.), on a 63 km reach of the Virgin River, U.S. Dispersal was measured by satellite imagery, ground surveys and phenocams. Pixels from the Moderate Resolution Imaging Spectrometer (MODIS) sensors on the Terra satellite showed a sharp drop in NDVI in midsummer followed by recovery, correlated with defoliation events as revealed in networked digital camera images and ground surveys. Ground surveys and MODIS imagery showed that beetle damage progressed downstream at a rate of about 25 km yr −1 in 2010 and 2011, producing a 50% reduction in saltcedar leaf area index and evapotranspiration by 2012, as estimated by algorithms based on MODIS Enhanced Vegetation Index values and local meteorological data for Mesquite, Nevada. This reduction is the equivalent of 10.4% of mean annual river flows on this river reach. Our results confirm other observations that saltcedar beetles are dispersing much faster than originally predicted in pre-release biological assessments, presenting new challenges and opportunities for land, water and wildlife managers on western rivers. Despite relatively coarse resolution (250 m) and gridding artifacts, single MODIS pixels can be useful in tracking the effects of defoliating insects in riparian corridors.

Arizona, Nevada↗

Characterizing regional soil mineral composition using spectroscopyand geostatistics

This work aims at improving the mapping of major mineral variability at regional scale using scale-dependent spatial variability observed in remote sensing data. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data and statistical methods were combined with laboratory-based mineral characterization of field samples to create maps of the distributions of clay, mica and carbonate minerals and their abundances. The Material Identification and Characterization Algorithm (MICA) was used to identify the spectrally-dominant minerals in field samples; these results were combined with ASTER data using multinomial logistic regression to map mineral distributions. X-ray diffraction (XRD)was used to quantify mineral composition in field samples. XRD results were combined with ASTER data using multiple linear regression to map mineral abundances. We testedwhether smoothing of the ASTER data to match the scale of variability of the target sample would improve model correlations. Smoothing was donewith Fixed Rank Kriging (FRK) to represent the mediumand long-range spatial variability in the ASTER data. Stronger correlations resulted using the smoothed data compared to results obtained with the original data. Highest model accuracies came from using both medium and long-range scaled ASTER data as input to the statistical models. High correlation coefficients were obtained for the abundances of calcite and mica (R2 = 0.71 and 0.70, respectively). Moderately-high correlation coefficients were found for smectite and kaolinite (R2 = 0.57 and 0.45, respectively). Maps of mineral distributions, obtained by relating ASTER data to MICA analysis of field samples, were found to characterize major soil mineral variability (overall accuracies for mica, smectite and kaolinite were 76%, 89% and 86% respectively). The results of this study suggest that the distributions of minerals and their abundances derived using FRK-smoothed ASTER data more closely match the spatial variability of soil and environmental properties at regional scale.

Remote Sensing of Environment↗

Mapping irrigated areas in Afghanistan over the past decade using MODIS NDVI

Agricultural production capacity contributes to food security in Afghanistan and is largely dependent on irrigated farming, mostly utilizing surface water fed by snowmelt. Because of the high contribution of irrigated crops (> 80%) to total agricultural production, knowing the spatial distribution and year-to-year variability in irrigated areas is imperative to monitoring food security for the country. We used 16-day composites of the Normalized Difference Vegetation Index (NDVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor to create 23-point time series for each year from 2000 through 2013. Seasonal peak values and time series were used in a threshold-dependent decision tree algorithm to map irrigated areas in Afghanistan for the last 14 years. In the absence of ground reference irrigated area information, we evaluated these maps with the irrigated areas classified from multiple snapshots of the landscape during the growing season from Landsat 5 optical and thermal sensor images. We were able to identify irrigated areas using Landsat imagery by selecting as irrigated those areas with Landsat-derived NDVI greater than 0.30–0.45, depending on the date of the Landsat image and surface temperature less than or equal to 310 Kelvin (36.9 ° C). Due to the availability of Landsat images, we were able to compare with the MODIS-derived maps for four years: 2000, 2009, 2010, and 2011. The irrigated areas derived from Landsat agreed well r 2 = 0.91 with the irrigated areas derived from MODIS, providing confidence in the MODIS NDVI threshold approach. The maps portrayed a highly dynamic irrigated agriculture practice in Afghanistan, where the amount of irrigated area was largely determined by the availability of surface water, especially snowmelt, and varied by as much as 30% between water surplus and water deficit years. During the past 14 years, 2001, 2004, and 2008 showed the lowest levels of irrigated area (~ 1.5 million hectares), attesting to the severe drought conditions in those years, whereas 2009, 2012 and 2013 registered the largest irrigated area (~ 2.5 million hectares) due to record snowpack and snowmelt in the region. The model holds promise the ability to provide near-real-time (by the end of the growing seasons) estimates of irrigated area, which are beneficial for food security monitoring as well as subsequent decision making for the country. While the model is developed for Afghanistan, it can be adopted with appropriate adjustments in the derived threshold values to map irrigated areas elsewhere.

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,↗

Validation of the USGS Landsat Burned Area Essential Climate Variable (BAECV) across the conterminous United States

The Landsat Burned Area Essential Climate Variable (BAECV), developed by the U.S. Geological Survey (USGS), capitalizes on the long temporal availability of Landsat imagery to identify burned areas across the conterminous United States (CONUS) (1984–2015). Adequate validation of such products is critical for their proper usage and interpretation. Validation of coarse-resolution products often relies on independent data derived from moderate-resolution sensors (e.g., Landsat). Validation of Landsat products, in turn, is challenging because there is no corresponding source of high-resolution, multispectral imagery that has been systematically collected in space and time over the entire temporal extent of the Landsat archive. Because of this, comparison between high-resolution images and Landsat science products can help increase user's confidence in the Landsat science products, but may not, alone, be adequate. In this paper, we demonstrate an approach to systematically validate the Landsat-derived BAECV product. Burned area extent was mapped for Landsat image pairs using a manually trained semi-automated algorithm that was manually edited across 28 path/rows and five different years (1988, 1993, 1998, 2003, 2008). Three datasets were independently developed by three analysts and the datasets were integrated on a pixel by pixel basis in which at least one to all three analysts were required to agree a pixel was burned. We found that errors within our Landsat reference dataset could be minimized by using the rendition of the dataset in which pixels were mapped as burned if at least two of the three analysts agreed. BAECV errors of omission and commission for the detection of burned pixels averaged 42% and 33%, respectively for CONUS across all five validation years. Errors of omission and commission were lowest across the western CONUS, for example in the shrub and scrublands of the Arid West (31% and 24%, respectively), and highest in the grasslands and agricultural lands of the Great Plains in central CONUS (62% and 57%, respectively). The BAECV product detected most (> 65%) fire events > 10 ha across the western CONUS (Arid and Mountain West ecoregions). Our approach and results demonstrate that a thorough validation of Landsat science products can be completed with independent Landsat-derived reference data, but could be strengthened by the use of complementary sources of high-resolution data.

Remote Sensing of Environment↗

Remote sensing of dryland ecosystem structure and function: Progress, challenges, and opportunities

Drylands make up roughly 40% of the Earth's land surface, and billions of people depend on services provided by these critically important ecosystems. Despite their relatively sparse vegetation, dryland ecosystems are structurally and functionally diverse, and emerging evidence suggests that these ecosystems play a dominant role in the trend and variability of the terrestrial carbon sink. More, drylands are highly sensitive to climate and are likely to have large, non-linear responses to hydroclimatic change. Monitoring the spatiotemporal dynamics of dryland ecosystem structure (e.g., leaf area index) and function (e.g., primary production and evapotranspiration) is therefore a high research priority. Yet, dryland remote sensing is defined by unique challenges not typically encountered in mesic or humid regions. Major challenges include low vegetation signal-to-noise ratios, high soil background reflectance, presence of photosynthetic soils (i.e., biological soil crusts), high spatial heterogeneity from plot to regional scales, and irregular growing seasons due to unpredictable seasonal rainfall and frequent periods of drought. Additionally, there is a relative paucity of continuous, long-term measurements in drylands, which impedes robust calibration and evaluation of remotely-sensed dryland data products. Due to these issues, remote sensing techniques developed in other ecosystems or for global application often result in inaccurate, poorly constrained estimates of dryland ecosystem structural and functional dynamics. Here, we review past achievements and current progress in remote sensing of dryland ecosystems, including a detailed discussion of the major challenges associated with remote sensing of key dryland structural and functional dynamics. We then identify strategies aimed at leveraging new and emerging opportunities in remote sensing to overcome previous challenges and more accurately contextualize drylands within the broader Earth system. Specifically, we recommend: 1) Exploring novel combinations of sensors and techniques (e.g., solar-induced fluorescence, thermal, microwave, hyperspectral, and LiDAR) across a range of spatiotemporal scales to gain new insights into dryland structural and functional dynamics; 2) utilizing near-continuous observations from new-and-improved geostationary satellites to capture the rapid responses of dryland ecosystems to diurnal variation in water stress; 3) expanding ground observational networks to better represent the heterogeneity of dryland systems and enable robust calibration and evaluation; 4) developing algorithms that are specifically tuned to dryland ecosystems by utilizing expanded ground observational network data; and 5) coupling remote sensing observations with process-based models using data assimilation to improve mechanistic understanding of dryland ecosystem dynamics and to better constrain ecological forecasts and long-term projections.

Remote Sensing of Environment↗