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

A Near Four-Decade Time Series Shows the Hawaiian Islands Have Been Browning Since the 1980s

The Hawaiian Islands have been identified as a global biodiversity hotspot. We examine the Normalized Difference Vegetation Index (NDVI) using Climate Data Records products (0.05 × 0.05°) to identify significant differences in NDVI between neutral El Niño-Southern Oscillation years (1984, 2019) and significant long-term changes over the entire time series (1982-2019) for the Hawaiian Islands and six land cover classes. Overall, there has been a significant decline in NDVI (i.e., browning) across the Hawaiian Islands from 1982 to 2019 with the islands of Lāna'i and Hawai'i experiencing the greatest decreases in NDVI (≥44%). All land cover classes significantly decreased in NDVI for most months, especially during the wet season month of March. Native vegetation cover across all islands also experienced significant declines in NDVI, with the leeward, southwestern side of the island of Hawai'i experiencing the greatest declines. The long-term trends in the annual total precipitation and annual mean Palmer Drought Severity Index (PDSI) for 1982-2019 on the Hawaiian Islands show significant concurrent declines. Primarily positive correlations between the native ecosystem NDVI and precipitation imply that significant decreases in precipitation may exacerbate the decrease in NDVI of native ecosystems. NDVI-PDSI correlations were primarily negative on the windward side of the islands and positive on the leeward sides, suggesting a higher sensitivity to drought for leeward native ecosystems. Multi-decadal time series and spatially explicit data for native landscapes provide natural resource managers with long-term trends and monthly changes associated with vegetation health and stability.

Drivers and trends of global soil microbial carbon over two decades

Soil microorganisms are central to sustain soil functions and services, like carbon and nutrient cycling. Currently, we only have a limited understanding of the spatial-temporal dynamics of soil microorganisms, restricting our ability to assess long-term effects of climate and land-cover change on microbial roles in soil biogeochemistry. This study assesses the temporal trends in soil microbial biomass carbon and identifies the main drivers of biomass change regionally and globally to detect the areas sensitive to these environmental factors. Here, we combined a global soil microbial biomass carbon data set, random forest modelling, and environmental layers to predict spatial-temporal dynamics of microbial biomass carbon stocks from 1992 to 2013. Soil microbial biomass carbon stocks decreased globally by 3.4 ± 3.0% (mean ± 95% CI) between 1992 and 2013 for the predictable regions, equivalent to 149 Mt being lost over the period, or ~1‰ of soil C. Northern areas with high soil microbial carbon stocks experienced the strongest decrease, mostly driven by increasing temperatures. In contrast, land-cover change was a weaker global driver of change in microbial carbon, but had, in some cases, important regional effects.

Spatiotemporal Dynamics of NDVI, Soil Moisture and ENSO in Tropical South America

We evaluated the coupled dynamics of vegetation dynamics (NDVI) and soil moisture (SMOS) at monthly resolution over different regions of tropical South America and the effects of the Eastern Pacific (EP) and the Central Pacific (CP) El Niño–Southern Oscillation (ENSO) events. We used linear Pearson cross-correlation, wavelet and cross wavelet analysis (CWA) and three nonlinear causality methods: ParrCorr, GPDC and PCMCIplus. Results showed that NDVI peaks when SMOS is transitioning from maximum to minimum monthly values, which confirms the role of SMOS in the hydrological dynamics of the Amazonian greening up during the dry season. Linear correlations showed significant positive values when SMOS leads NDVI by 1–3 months. Wavelet analysis evidenced strong 12- and 64-month frequency bands throughout the entire record length, in particular for SMOS, whereas the CWA analyses indicated that both variables exhibit a strong coherency at a wide range of frequency bands from 2 to 32 months. Linear and nonlinear causality measures also showed that ENSO effects are greater on SMOS. Lagged cross-correlations displayed that western (eastern) regions are more associated with the CP (EP), and that the effects of ENSO manifest as a travelling wave over time, from northwest (earlier) to southeast (later) over tropical South America and the Amazon River basin. The ParrCorr and PCMCIplus methods produced the most coherent results, and allowed us to conclude that: (1) the nonlinear temporal persistence (memory) of soil moisture is stronger than that of NDVI; (2) the existence of two-way nonlinear causalities between NDVI and SMOS; (3) diverse causal links between both variables and the ENSO indices: CP (7/12 with ParrCorr; 6/12 with PCMCIplus), and less with EP (5/12 with ParrCorr; 3/12 with PCMCIplus).

Monitoring for Changes in Spring Phenology at Both Temporal and Spatial Scales Based on MODIS LST Data in South Korea

This study aims to monitor spatiotemporal changes of spring phenology using the green-up start dates based on the accumulated growing degree days (AGDD) and the enhanced vegetation index (EVI), which were deducted from moderate resolution imaging spectroradiometer (MODIS) land surface temperature (LST) data. The green-up start dates were extracted from the MODIS-derived AGDD and EVI for 30 Mongolian oak (Quercus mongolica Fisch.) stands throughout South Korea. The relationship between green-up day of year needed to reach the AGDD threshold (DoYAGDD) and air temperature was closely maintained in data in both MODIS image interpretation and from 93 meteorological stations. Leaf green-up dates of Mongolian oak based on the AGDD threshold obtained from the records measured at five meteorological stations during the last century showed the same trend as the result of cherry observed visibly. Extrapolating the results, the spring onset of Mongolian oak and cherry has become earlier (14.5 ± 4.3 and 10.7 ± 3.6 days, respectively) with the rise of air temperature over the last century. The temperature in urban areas was consistently higher than that in the forest and the rural areas and the result was reflected on the vegetation phenology. Our study expanded the scale of the study on spring vegetation phenology spatiotemporally by combining satellite images with meteorological data. We expect our findings could be used to predict long-term changes in ecosystems due to climate change.

PACO: Python-Based Atmospheric COrrection

The atmospheric correction of satellite images based on radiative transfer calculations is a prerequisite for many remote sensing applications. The software package ATCOR, developed at the German Aerospace Center (DLR), is a versatile atmospheric correction software, capable of processing data acquired by many different optical satellite sensors. Based on this well established algorithm, a new Python-based atmospheric correction software has been developed to generate L2A products of Sentinel-2, Landsat-8, and of new space-based hyperspectral sensors such as DESIS (DLR Earth Sensing Imaging Spectrometer) and EnMAP (Environmental Mapping and Analysis Program). This paper outlines the underlying algorithms of PACO, and presents the validation results by comparing L2A products generated from Sentinel-2 L1C images with in situ (AERONET and RadCalNet) data within VNIR-SWIR spectral wavelengths range.

Can We Use the QA4ECV Black-sky Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) using AVHRR Surface Reflectance to Assess Terrestrial Global Change?

NOAA platforms provide the longest period of terrestrial observation since the 1980s. The progress in calibration, atmospheric corrections and physically based land retrieval offers the opportunity to reprocess these data for extending terrestrial product time series. Within the Quality Assurance for Essential Climate Variables (QA4ECV) project, the black-sky Joint Research Centre (JRC)-fraction of absorbed photosynthetically active radiation (FAPAR) algorithm was developed for the AVHRR sensors on-board NOAA-07 to -16 using the Land Surface Reflectance Climate Data Record. The retrieval algorithm was based on the radiative transfer theory, and uncertainties were included in the products. We proposed a time and spatial composite for providing both 10-day and monthly products at 0.05º × 0.05º. Quality control and validation were achieved through benchmarking against third-party products, including Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) datasets produced with the same retrieval algorithm. Past ground-based measurements, providing a proxy of FAPAR, showed good agreement of seasonality values over short homogeneous canopies and mixed vegetation. The average difference between SeaWiFS and QA4ECV monthly products over 2002–2005 is about 0.075 with a standard deviation of 0.091. We proposed a monthly linear bias correction that reduced these statistics to 0.02 and 0.001. The complete harmonized long-term time series was then used to address its fitness for the purpose of analysis of global terrestrial change.

Validation of Earth Observation Time-Series: A Review for Large-Area and Temporally Dense Land Surface Products

Large-area remote sensing time-series offer unique features for the extensive investigation of our environment. Since various error sources in the acquisition chain of datasets exist, only properly validated results can be of value for research and downstream decision processes. This review presents an overview of validation approaches concerning temporally dense time-series of land surface geo-information products that cover the continental to global scale. Categorization according to utilized validation data revealed that product intercomparisons and comparison to reference data are the conventional validation methods. The reviewed studies are mainly based on optical sensors and orientated towards global coverage, with vegetation-related variables as the focus. Trends indicate an increase in remote sensing-based studies that feature long-term datasets of land surface variables. The hereby corresponding validation efforts show only minor methodological diversification in the past two decades. To sustain comprehensive and standardized validation efforts, the provision of spatiotemporally dense validation data in order to estimate actual differences between measurement and the true state has to be maintained. The promotion of novel approaches can, on the other hand, prove beneficial for various downstream applications, although typically only theoretical uncertainties are provided.

Improving the AVHRR Long Term Data Record BRDF Correction

The Long Term Data Record (LTDR) project has the goal of developing a quality and consistent surface reflectance product from coarse resolution optical sensors. This paper focuses on the Advanced Very High Resolution Radiometer (AVHRR) part of the record, using the Moderate Resolution Imaging Spectrometer (MODIS) instrument as a reference. When a surface reflectance time series is acquired from satellites with variable observation geometry, the directional variation generates an apparent noise which can be corrected by modeling the bidirectional reflectance distribution function (BRDF). The VJB (Vermote, Justice and Bréon, 2009) method estimates a target’s BRDF shape using 5 years of observation and corrects for directional effects maintaining the high temporal resolution of the measurement using the instantaneous Normalized Difference Vegetation Index (NDVI). The method was originally established on MODIS data but its viability and optimization for AVHRR data have not been fully explored. In this study we analyze different approaches to find the most robust way of applying the VJB correction to AVHRR data, considering that high noise in the red band (B1) caused by atmospheric effect makes the VJB method unstable. Firstly, our results show that for coarse spatial resolution, where the vegetation dynamics of the target don’t change significantly, deriving BRDF parameters from 15+ years of observations reduces the average noise by up to 7% in the Near Infrared (NIR) band and 6% in the NDVI, in comparison to using 3-year windows. Secondly, we find that the VJB method can be modified for AVHRR data to improve the robustness of the correction parameters and decrease the noise by an extra 8% and 9% in the red and NIR bands with respect to using the classical VJB inversion. We do this by using the Stable method, which obtains the volumetric BRDF parameter (V) based on its NDVI dependency, and then obtains the geometric BRDF parameter (R) through the inversion of just one parameter.

Prior Season Crop Type Masks for Winter Wheat Yield Forecasting: A US Case Study

Monitoring and forecasting crop yields is a critical component of understanding and better addressing global food security challenges. Detailed spatial information on crop-type distribution is fundamental for in-season crop condition monitoring and yields forecasting over large agricultural areas, as it enables the extraction of crop-specific signals. Yet, the availability of such data within the growing season is often limited. Within this context, this study seeks to develop a practical approach to extract a crop-specific signal for yield forecasting in cases where crop rotations are prevalent, and detailed in-season information on crop type distribution is not available. We investigated the possibility of accurately forecasting winter wheat yields by using a counter-intuitive approach, which coarsens the spatial resolution of out-of-date detailed winter wheat masks and uses them in combination with easily accessibly coarse spatial resolution remotely sensed time series data. The main idea is to explore an optimal spatial resolution at which crop type changes will be negligible due to crop rotation (so a previous seasons’ mask, which is more readily available can be used) and an informative signal can be extracted, so it can be correlated to crop yields. The study was carried out in the United States of America (USA) and utilized multiple years of NASA Moderate Resolution Imaging Spectroradiometer (MODIS) data, US Department of Agriculture (USDA) National Agricultural Statistics Service (NASS) detailed wheat masks, and a regression-based winter wheat yield model. The results indicate that, in places where crop rotations were prevalent, coarsening the spatial scale of a crop type mask from the previous season resulted in a constant per-pixel wheat proportion over multiple seasons. This enables the consistent extraction of a crop-specific vegetation index time series that can be used for in-season monitoring and yield estimation over multiple years using a single mask. In the case of the USA, using a moderate resolution crop type mask from a previous season aggregated to 5 km resolution, resulted in a 0.7% tradeoff in accuracy relative to the control case where annually-updated detailed crop-type masks were available. These findings suggest that when detailed in-season data is not available, winter wheat yield can be accurately forecasted (within 10%) prior to harvest using a single, prior season crop mask and coarse resolution Normalized Difference Vegetation Index (NDVI) time series data.

Contribution of Land Surface Temperature (TCI) to Vegetation Health Index: A Comparative Study Using Clear Sky and All-Weather Climate Data Records

The Vegetation Health Index (VHI) is widely used for monitoring drought using satellite data. VHI depends on vegetation state and thermal stress, respectively assessed via (i) the Vegetation Condition Index (VCI) that usually relies on information from the visible and near infra-red parts of the spectrum (in the form of Normalized Difference Vegetation Index, NDVI); and (ii) the Thermal Condition Index (TCI), based on top of atmosphere thermal infrared (TIR) brightness temperature or on TIR-derived Land Surface Temperature (LST). VHI is then estimated as a weighted average of VCI and TCI. However, the optimum weights of the two components are usually not known and VHI is usually estimated attributing a weight of 0.5 to both. Using a previously developed methodology for the Euro-Mediterranean region, we show that the multi-scalar drought index (SPEI) may be used to obtain optimal weights for VCI and TCI over the area covered by Meteosat satellites that includes Africa, Europe, and part of South America. The procedure is applied using clear-sky Meteosat Climate Data Records (CDRs) and all-sky LST derived by combining satellite and reanalysis data. Results obtained present a coherent spatial distribution of VCI and TCI weights when estimated using clear- and all-sky LST. This study paves the way for the development of a future VHI near-real time operational product for drought monitoring based on information from Meteosat satellites.

Quality Assurance Framework Development Based on Six New ECV Data Products to Enhance User Confidence for Climate Applications

Data from Earth observation (EO) satellites are increasingly used to monitor the environment, understand variability and change, inform evaluations of climate model forecasts, and manage natural resources. Policymakers are progressively relying on the information derived from these datasets to make decisions on mitigating and adapting to climate change. These decisions should be evidence based, which requires confidence in derived products, as well as the reference measurements used to calibrate, validate, or inform product development. In support of the European Union’s Earth Observation Programmes Copernicus Climate Change Service (C3S), the Quality Assurance for Essential Climate Variables (QA4ECV) project fulfilled a gap in the delivery of climate quality satellite-derived datasets, by prototyping a generic system for the implementation and evaluation of quality assurance (QA) measures for satellite-derived ECV climate data record products. The project demonstrated the QA system on six new long-term, climate quality ECV data records for surface albedo, leaf area index (LAI), fraction of absorbed photosynthetically active radiation (FAPAR), nitrogen dioxide (NO2), formaldehyde (HCHO), and carbon monoxide (CO). The provision of standardised QA information provides data users with evidence-based confidence in the products and enables judgement on the fitness-for-purpose of various ECV data products and their specific applications.

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.

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.