At over 40 years, the Landsat satellites provide the longest temporal record of space-based land surface observations, and the successful 2013 launch of the Landsat-8 is continuing this legacy. Ideally, the Landsat data record should be consistent over the Landsat sensor series. The Landsat-8 Operational Land Imager (OLI) has improved calibration, signal to noise characteristics, higher 12-bit radiometric resolution, and spectrally narrower wavebands than the previous Landsat-7 Enhanced Thematic Mapper (ETM+). Reflective wavelength differences between the two Landsat sensors depend also on the surface reflectance and atmospheric state which are difficult to model comprehensively. The orbit and sensing geometries of the Landsat-8 OLI and Landsat-7 ETM+ provide swath edge overlapping paths sensed only one day apart. The overlap regions are sensed in alternating backscatter and forward scattering orientations so Landsat bi-directional reflectance effects are evident but approximately balanced between the two sensors when large amounts of time series data are considered. Taking advantage of this configuration a total of 59 million 30m corresponding sensor observations extracted from 6,317 Landsat-7 ETM+ and Landsat-8 OLI images acquired over three winter and three summer months for all the conterminous United States (CONUS) are compared. Results considering different stages of cloud and saturation filtering, and filtering to reduce one day surface state differences, demonstrate the importance of appropriate per-pixel data screening. Top of atmosphere (TOA) and atmospherically corrected surface reflectance for the spectrally corresponding visible, near infrared and shortwave infrared bands, and derived normalized difference vegetation index (NDVI), are compared and their differences quantified. On average the OLI TOA reflectance is greater than the ETM+ TOA reflectance for all bands, with greatest differences in the near-infrared (NIR) and the shortwave infrared bands due to the quite different spectral response functions between the sensors. The atmospheric correction reduces the mean difference in the NIR and shortwave infrared but increases the mean difference in the visible bands. Regardless of whether TOA or surface reflectance are used to generate NDVI, on average, for vegetated soil and vegetation surfaces (0 ≤ NDVI ≤ 1), the OLI NDVI is greater than the ETM+ NDVI. Statistical functions to transform between the comparable sensor bands and sensor NDVI values are presented so that the user community may apply them in their own research to improve temporal continuity between the Landsat-7 ETM+ and Landsat-8 OLI sensor data. The transformation functions were developed using ordinary least squares (OLS) regression and were fit quite reliably (r2 values >0.7 for the reflectance data and >0.9 for the NDVI data, p-values <0.0001).
Hydrologic soil groups (HSGs) are a fundamental component of the USDA curve-number (CN) method for estimation of rainfall runoff; yet these data are not readily available in a format or spatial-resolution suitable for regional- and global-scale modeling applications. We developed a globally consistent, gridded dataset defining HSGs from soil texture, bedrock depth, and groundwater. The resulting data product—HYSOGs250m—represents runoff potential at 250 m spatial resolution. Our analysis indicates that the global distribution of soil is dominated by moderately high runoff potential, followed by moderately low, high, and low runoff potential. Low runoff potential, sandy soils are found primarily in parts of the Sahara and Arabian Deserts. High runoff potential soils occur predominantly within tropical and sub-tropical regions. No clear pattern could be discerned for moderately low runoff potential soils, as they occur in arid and humid environments and at both high and low elevations. Potential applications of this data include CN-based runoff modeling, flood risk assessment, and as a covariate for biogeographical analysis of vegetation distributions.
Biomass burning is a global phenomenon and systematic burned area mapping is of increasing importance for science and applications. With high spatial resolution and novelty in band design, the recently launched Sentinel-2A satellite provides a new opportunity for moderate spatial resolution burned area mapping. This study examines the performance of the Sentinel-2A Multi Spectral Instrument (MSI) bands and derived spectral indices to differentiate between unburned and burned areas. For this purpose, five pairs of pre-fire and post-fire top of atmosphere (TOA reflectance) and atmospherically corrected (surface reflectance) images were studied. The pixel values of locations that were unburned in the first image and burned in the second image, as well as the values of locations that were unburned in both images which served as a control, were compared and the discrimination of individual bands and spectral indices were evaluated using parametric (transformed divergence) and non-parametric (decision tree) approaches. Based on the results, the most suitable MSI bands to detect burned areas are the 20 m near-infrared, short wave infrared and red-edge bands, while the performance of the spectral indices varied with location. The atmospheric correction only significantly influenced the separability of the visible wavelength bands. The results provide insights that are useful for developing Sentinel-2 burned area mapping algorithms.
12 April 2017 in International Journal of Digital Earth
A global operational land imager (GOLI) Landsat-8 daytime active fire detection algorithm is presented. It utilizes established contextual active fire detection approaches but takes advantage of the significant increase in fire reflectance in Landsat-8 band 7 (2.20 μm) relative to band 4 (0.66 μm). The detection thresholds are fixed and based on a statistical examination of 39 million non-burning Landsat-8 pixels. Multi-temporal tests based on band 7 reflectance and relative changes in normalized difference vegetation index in the previous six months are used to reduce commissions errors. The probabilities of active fire detection for the GOLI and two recent Landsat-8 active fire detection algorithms are simulated to provide insights into their performance with respect to the fire size and temperature. The algorithms are applied to 11 Landsat-8 images that encompass a range of burning conditions and environments. Commission and omission errors are assessed by visual interpretation of detected active fire locations and by examination of the Landsat-8 images and higher spatial resolution Google Earth imagery. The GOLI algorithm has lower omission and comparable commission errors than the recent Landsat-8 active fire detection algorithms. The GOLI algorithm has demonstrable potential for global application and is suitable for implementation with other Landsat-like reflective wavelength sensors.
Satellite based fire radiant energy retrievals are widely applied to assess biomass consumed and emissions at regional to global scales. A known potential source of uncertainty in biomass burning estimates arises from fuel moisture but this impact has not been quantified in previous studies. Controlled fire laboratory experiments are used in this study to examine the biomass consumed and the radiant energy release (Fire Radiative Energy, FRE, (MJ)) for western white pine needle fuels burned with water content (WC, unitless) from 0.01 to 0.14. Results indicate a significant relationship: FRE per kilogram of fuel consumed = −5.32 WC + 3.025 (r2 = 0.83, n = 24, P < 0.001) and imply that not taking into account fuel moisture variations in the assumed relationship between FRE and fuel consumed can lead to systematic biases. A methodological framework to derive a revised formula that enables the estimation of biomass consumed from FRE, which explicitly takes into account fuel water content, is presented.