The Harmonized Landsat and Sentinel-2 (HLS) project is a NASA initiative aiming to produce a Virtual Constellation (VC) of surface reflectance (SR) data acquired by the Operational Land Imager (OLI) ...and Multi-Spectral Instrument (MSI) aboard Landsat 8 and Sentinel-2 remote sensing satellites, respectively. The HLS products are based on a set of algorithms to obtain seamless products from both sensors (OLI and MSI): atmospheric correction, cloud and cloud-shadow masking, spatial co-registration and common gridding, bidirectional reflectance distribution function normalization and spectral bandpass adjustment. Three products are derived from the HLS processing chain: (i) S10: full resolution MSI SR at 10 m, 20 m and 60 m spatial resolutions; (ii) S30: a 30 m MSI Nadir BRDF (Bidirectional Reflectance Distribution Function)-Adjusted Reflectance (NBAR); (iii) L30: a 30 m OLI NBAR. All three products are processed for every Level-1 input products from Landsat 8/OLI (L1T) and Sentinel-2/MSI (L1C). As of version 1.3, the HLS data set covers 10.35 million km2 and spans from first Landsat 8 data (2013); Sentinel-2 data spans from October 2015.
The L30 and S30 show a good consistency with coarse spatial resolution products, in particular MODIS Collection 6 MCD09CMG products (overall deviations do not exceed 11%) that are used as a reference for quality assurance. The spatial co-registration of the HLS is improved compared to original Landsat 8 L1T and Sentinel-2A L1C products, for which misregistration issues between multi-temporal data are known. In particular, the resulting computed circular errors at 90% for the HLS product are 6.2 m and 18.8 m, for S10 and L30 products, respectively. The main known issue of the current data set remains the Sentinel-2 cloud mask with many cloud detection omissions. The cross-comparison with MODIS was used to flag products with most evident non-detected clouds. A time series outlier filtering approach is suggested to detect remaining clouds. Finally, several time series are presented to highlight the high potential of the HLS data set for crop monitoring.
•HLS is a surface reflectance record from the Landsat/Sentinel-2 Virtual Constellation.•HLS products cover 10 million km2, and correct for spectral and angular sensors differences.•HLS v1.3 products are consistent with MODIS products with a deviation lower than 11%.•The HLS relative spatial co-registration is around 0.62 pixels.•Sentinel-2 cloud mask is not optimal and is the main HLS v1.3 known issue.
The surface reflectance, i.e., satellite derived top of atmosphere (TOA) reflectance corrected for the temporally, spatially and spectrally varying scattering and absorbing effects of atmospheric ...gases and aerosols, is needed to monitor the land surface reliably. For this reason, the surface reflectance, and not TOA reflectance, is used to generate the greater majority of global land products, for example, from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Visible Infrared Imaging Radiometer Suite (VIIRS) sensors. Even if atmospheric effects are minimized by sensor design, atmospheric effects are still challenging to correct. In particular, the strong impact of aerosols in the visible and near infrared spectral range can be difficult to correct, because they can be highly discrete in space and time (e.g., smoke plumes) and because of the complex scattering and absorbing properties of aerosols that vary spectrally and with aerosol size, shape, chemistry and density.
This paper presents the Landsat 8 Operational Land Imager (OLI) atmospheric correction algorithm that has been developed using the Second Simulation of the Satellite Signal in the Solar Spectrum Vectorial (6SV) model, refined to take advantage of the narrow OLI spectral bands (compared to Thematic Mapper/Enhanced Thematic Mapper (TM/ETM+)), improved radiometric resolution and signal-to-noise. In addition, the algorithm uses the new OLI Coastal aerosol band (0.433–0.450μm), which is particularly helpful for retrieving aerosol properties, as it covers shorter wavelengths than the conventional Landsat, TM and ETM+ blue bands. A cloud and cloud shadow mask has also been developed using the “cirrus” band (1.360–1.390μm) available on OLI, and the thermal infrared bands from the Thermal Infrared Sensor (TIRS) instrument. The performance of the surface reflectance product from OLI is analyzed over the Aerosol Robotic Network (AERONET) sites using accurate atmospheric correction (based on in situ measurements of the atmospheric properties), by comparison with the MODIS Bidirectional Reflectance Distribution Function (BRDF) adjusted surface reflectance product and by comparison of OLI derived broadband albedo from United States Surface Radiation Budget Network (US SURFRAD) measurements. The results presented clearly show an improvement of Landsat 8 surface reflectance product over the ad-hoc Landsat 5/7 LEDAPS product.
•Landsat 8 OLI Atmospheric correction algorithm•Algorithm takes advantage of improved sensor characteristics.•A step forward compared to the ad-hoc Landsat 5/7 LEDAPS surface reflectance product
Maintaining consistent datasets of Surface Reflectance (SR) is an important challenge to ensure long-term quality of Climate Data Records. The Landsat 5 and 7 archives offer a unique data source to ...monitor globally the land surface at high spatial resolution. The Landsat-5 TM and Landsat-7 ETM+ SR products, derived from the on-demand processing Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS), require periodic evaluation to check the data consistency. Two evaluation approaches are presented in this paper. The first approach used the Aerosol Robotic Network (AERONET) data set for the period 2000 to 2013 over 489 sites with 3600 Landsat-5 TM and Landsat-7 ETM+ scenes selected. For each scene, 10×10km subsets of LEDAPS-derived Landsat SR and AERONET-derived SR are compared. The latter are computed using Landsat top of atmosphere reflectance, AERONET measurements of atmospheric parameters, and the 6S radiative transfer model. Second, we introduce a methodology to cross-compare Landsat data and MODIS data acquired on the same day. The analysis is based on 4000 random Landsat scenes globally distributed from 2000 to 2013. This method includes: (i) a surface anisotropy adjustment, based on the VJB Bidirectional Reflectance Distribution Function (BRDF) method, to adjust Terra and Aqua MODIS data to Landsat 5 and 7 sun-view geometry, (ii) a spectral adjustment based on an artificial neural network trained with the PROSAIL vegetation radiative transfer model, to adjust MODIS data to TM and ETM+ spectral responses.
The overall results of both approaches show a good match in over 80% of the scenes, i.e. the TM and ETM+ SR uncertainty remained within the SR specification, defined as 0.05×SR+0.005. The worst results are found in the blue band used in LEDAPS to adjust the Aerosol Optical Thickness (AOT). The MODIS-Landsat SR cross-comparison confirms the utility of a BRDF adjustment method to decrease the scattering between Landsat sensors and MODIS sensors (Terra and Aqua). The spectral adjustment removes part of the biases related to spectral response differences. Global analysis is used to identify AOT retrieval issues over specific scenes, mostly over bright surfaces. From 2000 to 2013, no significant temporal variation of the performance is detected, which enhanced the consistency of LEDAPS-derived surface reflectance data set.
•We evaluate LEDAPS-derived Landsat-5 TM and −7 ETM+ surface reflectance (SR).•We evaluate the theoretical performances of LEDAPS over 489 AERONET sites.•We build a methodology to compare MODIS and Landsat SR with BRDF and spectral adjustments.•The Landsat SR uncertainty is below the specification (0.05ρ+0.005) except for blue band.•The Landsat SR uncertainty is stable through time.
Global land change from 1982 to 2016 Song, Xiao-Peng; Hansen, Matthew C; Stehman, Stephen V ...
Nature (London),
08/2018, Volume:
560, Issue:
7720
Journal Article
Peer reviewed
Open access
Land change is a cause and consequence of global environmental change
. Changes in land use and land cover considerably alter the Earth's energy balance and biogeochemical cycles, which contributes ...to climate change and-in turn-affects land surface properties and the provision of ecosystem services
. However, quantification of global land change is lacking. Here we analyse 35 years' worth of satellite data and provide a comprehensive record of global land-change dynamics during the period 1982-2016. We show that-contrary to the prevailing view that forest area has declined globally
-tree cover has increased by 2.24 million km
(+7.1% relative to the 1982 level). This overall net gain is the result of a net loss in the tropics being outweighed by a net gain in the extratropics. Global bare ground cover has decreased by 1.16 million km
(-3.1%), most notably in agricultural regions in Asia. Of all land changes, 60% are associated with direct human activities and 40% with indirect drivers such as climate change. Land-use change exhibits regional dominance, including tropical deforestation and agricultural expansion, temperate reforestation or afforestation, cropland intensification and urbanization. Consistently across all climate domains, montane systems have gained tree cover and many arid and semi-arid ecosystems have lost vegetation cover. The mapped land changes and the driver attributions reflect a human-dominated Earth system. The dataset we developed may be used to improve the modelling of land-use changes, biogeochemical cycles and vegetation-climate interactions to advance our understanding of global environmental change
.
This paper briefly describes the land surface reflectance product (MOD09), the current Moderate Resolution Imaging Spectroradiometer (MODIS) atmospheric correction (AC) algorithm and its recent ...updates, and provides the evaluation of the algorithm performance and product quality. The accuracy of the AC algorithm has been significantly improved owing to the use of the accurate Second Simulation of a Satellite Signal in the Solar Spectrum, Vector (6SV) radiative transfer code and a better retrieval of aerosol properties by a refined internal aerosol inversion algorithm. The Collection 5 MOD09 surface reflectance product computed by the improved AC algorithm was analyzed for the year of 2003 through the comparison with a reference data set created with the help of Aerosol Robotic Network (AERONET) measurements and the 6SV code simulations. In general, the MOD09 product demonstrated satisfactory quality in all used MODIS bands except for band 3 (470 nm), which is used for aerosol inversion. The impact of uncertainties in MOD09 upon the downstream product, such as vegetation indices and albedo, was also evaluated.
Tropical forests are global epicentres of biodiversity and important modulators of climate change, and are mainly constrained by rainfall patterns. The severe short-term droughts that occurred ...recently in Amazonia have drawn attention to the vulnerability of tropical forests to climatic disturbances. The central African rainforests, the second-largest on Earth, have experienced a long-term drying trend whose impacts on vegetation dynamics remain mostly unknown because in situ observations are very limited. The Congolese forest, with its drier conditions and higher percentage of semi-evergreen trees, may be more tolerant to short-term rainfall reduction than are wetter tropical forests, but for a long-term drought there may be critical thresholds of water availability below which higher-biomass, closed-canopy forests transition to more open, lower-biomass forests. Here we present observational evidence for a widespread decline in forest greenness over the past decade based on analyses of satellite data (optical, thermal, microwave and gravity) from several independent sensors over the Congo basin. This decline in vegetation greenness, particularly in the northern Congolese forest, is generally consistent with decreases in rainfall, terrestrial water storage, water content in aboveground woody and leaf biomass, and the canopy backscatter anomaly caused by changes in structure and moisture in upper forest layers. It is also consistent with increases in photosynthetically active radiation and land surface temperature. These multiple lines of evidence indicate that this large-scale vegetation browning, or loss of photosynthetic capacity, may be partially attributable to the long-term drying trend. Our results suggest that a continued gradual decline of photosynthetic capacity and moisture content driven by the persistent drying trend could alter the composition and structure of the Congolese forest to favour the spread of drought-tolerant species.
The seasonality of sunlight and rainfall regulates net primary production in tropical forests. Previous studies have suggested that light is more limiting than water for tropical forest productivity, ...consistent with greening of Amazon forests during the dry season in satellite data. We evaluated four potential mechanisms for the seasonal green-up phenomenon, including increases in leaf area or leaf reflectance, using a sophisticated radiative transfer model and independent satellite observations from lidar and optical sensors. Here we show that the apparent green up of Amazon forests in optical remote sensing data resulted from seasonal changes in near-infrared reflectance, an artefact of variations in sun-sensor geometry. Correcting this bidirectional reflectance effect eliminated seasonal changes in surface reflectance, consistent with independent lidar observations and model simulations with unchanging canopy properties. The stability of Amazon forest structure and reflectance over seasonal timescales challenges the paradigm of light-limited net primary production in Amazon forests and enhanced forest growth during drought conditions. Correcting optical remote sensing data for artefacts of sun-sensor geometry is essential to isolate the response of global vegetation to seasonal and interannual climate variability.
Since first light in early 2000, operational global quantitative retrievals of aerosol properties over land have been made from Moderate Resolution Imaging Spectroradiometer (MODIS) observed spectral ...reflectance. These products have been continuously evaluated and validated, and opportunities for improvements have been noted. We have replaced the surface reflectance assumptions, the set of aerosol model optical properties, and the aerosol lookup table (LUT). This second‐generation operational algorithm performs a simultaneous inversion of two visible (0.47 and 0.66 μm) and one shortwave‐IR (2.12 μm) channel, making use of the coarse aerosol information content contained in the 2.12 μm channel. Inversion of the three channels yields three nearly independent parameters, the aerosol optical depth (τ) at 0.55 μm, the nondust or fine weighting (η), and the surface reflectance at 2.12 μm. Retrievals of small‐magnitude negative τ values (down to −0.05) are considered valid, thus balancing the statistics of τ in near zero τ conditions. Preliminary validation of this algorithm shows much improved retrievals of τ, where the MODIS/Aerosol Robotic Network τ (at 0.55 μm) regression has an equation of: y = 1.01x + 0.03, R = 0.90. Global mean τ for the test bed is reduced from ∼0.28 to ∼0.21.
In- land surface models, which are used to evaluate the role of vegetation in the context of global climate change and variability, LAI and FAPAR play a key role, specifically with respect to the ...carbon and water cycles. The AVHRR-based LAI/FAPAR dataset offers daily temporal resolution, an improvement over previous products. This climate data record is based on a carefully calibrated and corrected land surface reflectance dataset to provide a high-quality, consistent time-series suitable for climate studies. It spans from mid-1981 to the present. Further, this operational dataset is available in near real-time allowing use for monitoring purposes. The algorithm relies on artificial neural networks calibrated using the MODIS LAI/FAPAR dataset. Evaluation based on cross-comparison with MODIS products and in situ data show the dataset is consistent and reliable with overall uncertainties of 1.03 and 0.15 for LAI and FAPAR, respectively. However, a clear saturation effect is observed in the broadleaf forest biomes with high LAI (>4.5) and FAPAR (>0.8) values.
In urban environments, aerosol distributions may change rapidly due to building and transport infrastructure and human population density variations. The recent availability of medium resolution ...Landsat-8 and Sentinel-2 satellite data provide the opportunity for aerosol optical depth (AOD) estimation at higher spatial resolution than provided by other satellites. A year of 30 m Landsat-8 and 10 m Sentinel-2A AOD data retrieved using the Land Surface Reflectance Code (LaSRC) were compared with coincident ground-based Aerosol Robotic Network (AERONET) Version 3 AOD data for 20 Chinese cities. Stringent selection criteria were used to select contemporaneous data - only satellite and AERONET data acquired within 10 minutes were considered. The average satellite retrieved AOD over a 1470 m × 1470 m window centered on each AERONET site was derived to capture fine scale urban AOD variations. AERONET Level 1.5 (cloud-screened) and Level 2.0 (cloud-screened and also quality assured) data were considered. For the 20 urban AERONET sites in 2016 there were 106 (Level 1.5) and 67 (Level 2.0) Landsat-8 AERONET AOD contemporaneous data pairs, and 118 (Level 1.5) and 89 (Level 2.0) Sentinel-2A AOD data pairs. The greatest AOD values (>1.5) occurred in Beijing, suggesting that the Chinese capital was one of the most polluted cities in China in 2016. The LaSRC Landsat-8 and Sentinel-2A AOD retrievals agreed well with the AERONET AOD data (linear regression slopes > 0.96; coefficient of determination r2 > 0.90; root mean square deviation < 0.175) and demonstrate that the LaSRC is an effective and applicable medium resolution AOD retrieval algorithm over urban environments. The Sentinel-2A AOD retrievals had better accuracy than the Landsat-8 AOD retrievals, which is consistent with previously published research. The implications of the research and the potential for urban aerosol monitoring by combining the freely available Landsat-8 and Sentinel-2 satellite data are discussed.