Evaluation of Landsat-8 and Sentinel-2A aerosol optical depth retrievals across Chinese cities and implications for medium spatial resolution urban aerosol monitoring.

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.

Author Correction: Global land change from 1982 to 2016.

In this Letter, errors in Supplementary Table 1 have been corrected.

Global land change from 1982 to 2016.

Land change is a cause and consequence of global environmental change1,2. 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 services1-4. 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 globally5-tree cover has increased by 2.24 million km2 (+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 km2 (-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 change1-4,6.

Transitioning from MODIS to VIIRS: an analysis of inter-consistency of NDVI data sets for agricultural monitoring.

The Visible/Infrared Imager/Radiometer Suite (VIIRS) aboard the Suomi National Polar-orbiting Partnership (S-NPP) satellite was launched in 2011, in part to provide continuity with the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument aboard National Aeronautics and Space Administration's (NASA) Terra and Aqua remote sensing satellites. The VIIRS will eventually replace MODIS for both land science and applications and add to the coarse-resolution, long term data record. It is, therefore, important to provide the user community with an assessment of the consistency of equivalent products from the two sensors. For this study, we do this in the context of example agricultural monitoring applications. Surface reflectance that is routinely delivered within the M{O,Y}D09 and VNP09 series of products provide critical input for generating downstream products. Given the range of applications utilizing the normalized difference vegetation index (NDVI) generated from M{O,Y}D09 and VNP09 products and the inherent differences between MODIS and VIIRS sensors in calibration, spatial sampling, and spectral bands, the main objective of this study is to quantify uncertainties related the transitioning from using MODIS to VIIRS-based NDVI's. In particular, we compare NDVI's derived from two sets of Level 3 MYD09 and VNP09 products with various spatial-temporal characteristics, namely 8-day composites at 500 m spatial resolution and daily Climate Modelling Grid (CMG) images at 0.05° spatial resolution. Spectral adjustment of VIIRS I1 (red) and I2 (near infra-red - NIR) bands to match MODIS/Aqua b1 (red) and b2 (NIR) bands is performed to remove a bias between MODIS and VIIRS-based red, NIR, and NDVI estimates. Overall, red reflectance, NIR reflectance, NDVI uncertainties were 0.014, 0.029 and 0.056 respectively for the 500 m product and 0.013, 0.016 and 0.032 for the 0.05° product. The study shows that MODIS and VIIRS NDVI data can be used interchangeably for applications with an uncertainty of less than 0.02 to 0.05, depending on the scale of spatial aggregation, which is typically the uncertainty of the individual dataset.

Atmospheric Correction Inter-comparison eXercise.

The Atmospheric Correction Inter-comparison eXercise (ACIX) is an international initiative with the aim to analyse the Surface Reflectance (SR) products of various state-of-the-art atmospheric correction (AC) processors. The Aerosol Optical Thickness (AOT) and Water Vapour (WV) are also examined in ACIX as additional outputs of an AC processing. In this paper, the general ACIX framework is discussed; special mention is made of the motivation to initiate this challenge, the inter-comparison protocol and the principal results. ACIX is free and open and every developer was welcome to participate. Eventually, 12 participants applied their approaches to various Landsat-8 and Sentinel-2 image datasets acquired over sites around the world. The current results diverge depending on the sensors, products and sites, indicating their strengths and weaknesses. Indeed, this first implementation of processor inter-comparison was proven to be a good lesson for the developers to learn the advantages and limitations of their approaches. Various algorithm improvements are expected, if not already implemented, and the enhanced performances are yet to be investigated in future ACIX experiments.

Combined use of Landsat-8 and Sentinel-2A images for winter crop mapping and winter wheat yield assessment at regional scale.

Timely and accurate information on crop yield is critical to many applications within agriculture monitoring. Thanks to its coverage and temporal resolution, coarse spatial resolution satellite imagery has always been a source of valuable information for yield forecasting and assessment at national and regional scales. With availability of free images acquired by Landsat-8 and Sentinel-2 remote sensing satellites, it becomes possible to enable temporal resolution of an image every 3-5 days, and therefore, to develop next generation agriculture products at higher spatial resolution (30 m). This paper explores the combined use of Landsat-8 and Sentinel-2A for winter crop mapping and winter wheat assessment at regional scale. For the former, we adapt a previously developed approach for Moderate Resolution Imaging Spectroradiometer (MODIS) at 250 m resolution that allows automatic mapping of winter crops taking into account knowledge on crop calendar and without ground truth data. For the latter, we use a generalized winter wheat yield model that is based on NDVI-peak estimation and MODIS data, and further downscaled to be applicable at 30 m resolution. We show that integration of Landsat-8 and Sentinel-2A has a positive impact both for winter crop mapping and winter wheat yield assessment. In particular, the error of winter wheat yield estimates can be reduced up to 1.8 times comparing to the single satellite usage.

Multi-spectral misregistration of Sentinel-2A images: analysis and implications for potential applications.

This study aims at analyzing sub-pixel misregistration between multi-spectral images acquired by the Multi-Spectral Instrument (MSI) aboard Sentinel-2A remote sensing satellite, and exploring its potential for moving target and cloud detection. By virtue of its hardware design, MSI's detectors exhibit a parallax angle that leads to sub-pixel shifts that are corrected with special pre-processing routines. However, these routines do not correct shifts for moving and/or high altitude objects. In this letter, we apply a phase correlation approach to detect sub-pixel shifts between B2 (blue), B3 (green) and B4 (red) Sentinel-2A/MSI images. We show that shifts of more than 1.1 pixels can be observed for moving targets, such as airplanes and clouds, and can be used for cloud detection. We demonstrate that the proposed approach can detect clouds that are not identified in the built-in cloud mask provided within the Sentinel-2A Level-1C (L1C) product.

Automatic sub-pixel co-registration of Landsat-8 OLI and Sentinel-2A MSI images using phase correlation and machine learning based mapping.

This study investigates misregistration issues between Landsat-8/OLI and Sentinel-2A/MSI at 30 m resolution, and between multi-temporal Sentinel-2A images at 10 m resolution using a phase correlation approach and multiple transformation functions. Co-registration of 45 Landsat-8 to Sentinel-2A pairs and 37 Sentinel-2A to Sentinel-2A pairs were analyzed. Phase correlation proved to be a robust approach that allowed us to identify hundreds and thousands of control points on images acquired more than 100 days apart. Overall, misregistration of up to 1.6 pixels at 30 m resolution between Landsat-8 and Sentinel-2A images, and 1.2 pixels and 2.8 pixels at 10 m resolution between multi-temporal Sentinel-2A images from the same and different orbits, respectively, were observed. The non-linear Random Forest regression used for constructing the mapping function showed best results in terms of root mean square error (RMSE), yielding an average RMSE error of 0.07±0.02 pixels at 30 m resolution, and 0.09±0.05 and 0.15±0.06 pixels at 10 m resolution for the same and adjacent Sentinel-2A orbits, respectively, for multiple tiles and multiple conditions. A simpler 1st order polynomial function (affine transformation) yielded RMSE of 0.08±0.02 pixels at 30 m resolution and 0.12±0.06 (same Sentinel-2A orbits) and 0.20±0.09 (adjacent orbits) pixels at 10 m resolution.

A 30+ year AVHRR Land Surface Reflectance Climate Data Record and its application to wheat yield monitoring.

The Advanced Very High Resolution Radiometer (AVHRR) sensor provides a unique global remote sensing dataset that ranges from the 1980's to the present. Over the years, several efforts have been made on the calibration of the different instruments to establish a consistent land surface reflectance time-series and to augment the AVHRR data record with data from other sensors such as the Moderate Resolution Imaging Spectroradiometer (MODIS). In this paper, we present a summary of all the corrections applied to the AVHRR Surface Reflectance and NDVI Version 4 Product, developed in the framework of the National Oceanic and Atmospheric Administration (NOAA) Climate Data Record (CDR) program. These corrections result from assessment of the geo-location, improvement of the cloud masking and calibration monitoring. Additionally, we evaluate the performance of the surface reflectance over the AERONET sites by a cross-comparison with MODIS, which is an already validated product, and evaluation of a downstream Leaf Area Index (LAI) product. We demonstrate the utility of this long time-series by estimating the winter wheat yield over the USA. The methods developed by [1] and [2] are applied to both the MODIS and AVHRR data. Comparison of the results from both sensors during the MODIS-era shows the consistency of the dataset with similar errors of 10%. When applying the methods to AVHRR historical data from the 1980's, the results have errors equivalent to those derived from MODIS.

A Workflow for Automated Satellite Image Processing: from Raw VHSR Data to Object-Based Spectral Information for Smallholder Agriculture.

Earth Observation has become a progressively important source of information for land use and land cover services over the past decades. At the same time, an increasing number of reconnaissance satellites have been set in orbit with ever increasing spatial, temporal, spectral, and radiometric resolutions. The available bulk of data, fostered by open access policies adopted by several agencies, is setting a new landscape in remote sensing in which timeliness and efficiency are important aspects of data processing. This study presents a fully automated workflow able to process a large collection of very high spatial resolution satellite images to produce actionable information in the application framework of smallholder farming. The workflow applies sequential image processing, extracts meaningful statistical information from agricultural parcels, and stores them in a crop spectrotemporal signature library. An important objective is to follow crop development through the season by analyzing multi-temporal and multi-sensor images. The workflow is based on free and open-source software, namely R, Python, Linux shell scripts, the Geospatial Data Abstraction Library, custom FORTRAN, C++, and the GNU Make utilities. We tested and applied this workflow on a multi-sensor image archive of over 270 VHSR WorldView-2, -3, QuickBird, GeoEye, and RapidEye images acquired over five different study areas where smallholder agriculture prevails.

Preliminary analysis of the performance of the Landsat 8/OLI land surface reflectance product.

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.

Widespread decline of Congo rainforest greenness in the past decade.

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.

Amazon forests maintain consistent canopy structure and greenness during the dry season.

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.

Land and cryosphere products from Suomi NPP VIIRS: Overview and status.

[1] The Visible Infrared Imaging Radiometer Suite (VIIRS) instrument was launched in October 2011 as part of the Suomi National Polar-Orbiting Partnership (S-NPP). The VIIRS instrument was designed to improve upon the capabilities of the operational Advanced Very High Resolution Radiometer and provide observation continuity with NASA's Earth Observing System's Moderate Resolution Imaging Spectroradiometer (MODIS). Since the VIIRS first-light images were received in November 2011, NASA- and NOAA-funded scientists have been working to evaluate the instrument performance and generate land and cryosphere products to meet the needs of the NOAA operational users and the NASA science community. NOAA's focus has been on refining a suite of operational products known as Environmental Data Records (EDRs), which were developed according to project specifications under the National Polar-Orbiting Environmental Satellite System. The NASA S-NPP Science Team has focused on evaluating the EDRs for science use, developing and testing additional products to meet science data needs, and providing MODIS data product continuity. This paper presents to-date findings of the NASA Science Team's evaluation of the VIIRS land and cryosphere EDRs, specifically Surface Reflectance, Land Surface Temperature, Surface Albedo, Vegetation Indices, Surface Type, Active Fires, Snow Cover, Ice Surface Temperature, and Sea Ice Characterization. The study concludes that, for MODIS data product continuity and earth system science, an enhanced suite of land and cryosphere products and associated data system capabilities are needed beyond the EDRs currently available from the VIIRS.

MODIS derived fire characteristics and aerosol optical depth variations during the agricultural residue burning season, north India.

Agricultural residue burning is one of the major causes of greenhouse gas emissions and aerosols in the Indo-Ganges region. In this study, we characterize the fire intensity, seasonality, variability, fire radiative energy (FRE) and aerosol optical depth (AOD) variations during the agricultural residue burning season using MODIS data. Fire counts exhibited significant bi-modal activity, with peak occurrences during April-May and October-November corresponding to wheat and rice residue burning episodes. The FRE variations coincided with the amount of residues burnt. The mean AOD (2003-2008) was 0.60 with 0.87 (+1σ) and 0.32 (-1σ). The increased AOD during the winter coincided well with the fire counts during rice residue burning season. In contrast, the AOD-fire signal was weak during the summer wheat residue burning and attributed to dust and fossil fuel combustion. Our results highlight the need for 'full accounting of GHG's and aerosols', for addressing the air quality in the study area.

Radiative transfer codes for atmospheric correction and aerosol retrieval: intercomparison study.

Results are summarized for a scientific project devoted to the comparison of four atmospheric radiative transfer codes incorporated into different satellite data processing algorithms, namely, 6SV1.1 (second simulation of a satellite signal in the solar spectrum, vector, version 1.1), RT3 (radiative transfer), MODTRAN (moderate resolution atmospheric transmittance and radiance code), and SHARM (spherical harmonics). The performance of the codes is tested against well-known benchmarks, such as Coulson's tabulated values and a Monte Carlo code. The influence of revealed differences on aerosol optical thickness and surface reflectance retrieval is estimated theoretically by using a simple mathematical approach. All information about the project can be found at http://rtcodes.ltdri.org.

Validation of a vector version of the 6S radiative transfer code for atmospheric correction of satellite data. Part II. Homogeneous Lambertian and anisotropic surfaces.

This is the second part of the validation effort of the recently developed vector version of the 6S (Second Simulation of a Satellite Signal in the Solar Spectrum) radiative transfer code (6SV1), primarily used for the calculation of look-up tables in the Moderate Resolution Imaging Spectroradiometer (MODIS) atmospheric correction algorithm. The 6SV1 code was tested against a Monte Carlo code and Coulson's tabulated values for molecular and aerosol atmospheres bounded by different Lambertian and anisotropic surfaces. The code was also tested in scalar mode against the scalar code SHARM to resolve the previous 6S accuracy issues in the case of an anisotropic surface. All test cases were characterized by good agreement between the 6SV1 and the other codes: The overall relative error did not exceed 0.8%. The study also showed that ignoring the effects of radiation polarization in the atmosphere led to large errors in the simulated top-of-atmosphere reflectances: The maximum observed error was approximately 7.2% for both Lambertian and anisotropic surfaces.

Large seasonal swings in leaf area of Amazon rainforests.

Despite early speculation to the contrary, all tropical forests studied to date display seasonal variations in the presence of new leaves, flowers, and fruits. Past studies were focused on the timing of phenological events and their cues but not on the accompanying changes in leaf area that regulate vegetation-atmosphere exchanges of energy, momentum, and mass. Here we report, from analysis of 5 years of recent satellite data, seasonal swings in green leaf area of approximately 25% in a majority of the Amazon rainforests. This seasonal cycle is timed to the seasonality of solar radiation in a manner that is suggestive of anticipatory and opportunistic patterns of net leaf flushing during the early to mid part of the light-rich dry season and net leaf abscission during the cloudy wet season. These seasonal swings in leaf area may be critical to initiation of the transition from dry to wet season, seasonal carbon balance between photosynthetic gains and respiratory losses, and litterfall nutrient cycling in moist tropical forests.