Inconsistencies of interannual variability and trends in long-term satellite leaf area index products

Understanding the long-term performance of global satellite leaf area index (LAI) products is important for global change research. However, few effort has been devoted to evaluating the long-term time-series consistencies of LAI products. This study compared four long-term LAI products (GLASS, GLOBMAP, LAI3g, and TCDR) in terms of trends, interannual variabilities, and uncertainty variations from 1982 through 2011. This study also used four ancillary LAI products (GEOV1, MERIS, MODIS C5, and MODIS C6) from 2003 through 2011 to help clarify the performances of the four long-term LAI products. In general, there were marked discrepancies between the four long-term LAI products. During the pre-MODIS period (1982-1999), both linear trends and interannual variabilities of global mean LAI followed the order GLASS>LAI3g>TCDR>GLOBMAP. The GLASS linear trend and interannual variability were almost 4.5 times those of GLOBMAP. During the overlap period (2003-2011), GLASS and GLOBMAP exhibited a decreasing trend, TCDR no trend, and LAI3g an increasing trend. GEOV1, MERIS, and MODIS C6 also exhibited an increasing trend, but to a much smaller extent than that from LAI3g. During both periods, the R2 of detrended anomalies between the four long-term LAI products was smaller than 0.4 for most regions. Interannual variabilities of the four long-term LAI products were considerably different over the two periods, and the differences followed the order GLASS>LAI3g>TCDR>GLOBMAP. Uncertainty variations quantified by a collocation error model followed the same order. Our results indicate that the four long-term LAI products were neither intraconsistent over time nor interconsistent with each other. These inconsistencies may be due to NOAA satellite orbit changes and MODIS sensor degradation. Caution should be used in the interpretation of global changes derived from the four long-term LAI products

Implications of whole-disc DSCOVR EPIC spectral observations for estimating Earth's spectral reflectivity based on low-earth-orbiting and geostationary observations

Earth’s reflectivity is among the key parameters of climate research. National Aeronautics and Space Administration (NASA)’s Earth Polychromatic Imaging Camera (EPIC) onboard National Oceanic and Atmospheric Administration (NOAA)’s Deep Space Climate Observatory (DSCOVR) spacecraft provides spectral reflectance of the entire sunlit Earth in the near backscattering direction every 65 to 110 min. Unlike EPIC, sensors onboard the Earth Orbiting Satellites (EOS) sample reflectance over swaths at a specific local solar time (LST) or over a fixed area. Such intrinsic sampling limits result in an apparent Earth’s reflectivity. We generated spectral reflectance over sampling areas using EPIC data. The difference between the EPIC and EOS estimates is an uncertainty in Earth’s reflectivity. We developed an Earth Reflector Type Index (ERTI) to discriminate between major Earth atmosphere components: clouds, cloud-free ocean, bare and vegetated land. Temporal variations in Earth’s reflectivity are mostly determined by clouds. The sampling area of EOS sensors may not be sufficient to represent cloud variability, resulting in biased estimates. Taking EPIC reflectivity as a reference, low-earth-orbiting-measurements at the sensor-specific LST tend to overestimate EPIC values by 0.8% to 8%. Biases in geostationary orbiting approximations due to a limited sampling area are between −0.7% and 12%. Analyses of ERTI-based Earth component reflectivity indicate that the disagreement between EPIC and EOS estimates depends on the sampling area, observation time and vary between −10% and 23%.The NASA/GSFC DSCOVR project is funded by NASA Earth Science Division. W. Song, G. Yan, and X. Mu were also supported by the key program of National Natural Science Foundation of China (NSFC; Grant No. 41331171). This research was conducted and completed during a 13-month research stay of the lead author in the Department of Earth and Environment, Boston University as a joint Ph.D. student, which was supported by the Chinese Scholarship Council (201606040098). DSCOVR EPIC L1B data were obtained from the NASA Langley Research Center Atmospheric Science Data Center. The authors would like to thank the editor who handled this paper and the two anonymous reviewers for providing helpful and constructive comments and suggestions that significantly helped us improve the quality of this paper. (NASA Earth Science Division; 41331171 - key program of National Natural Science Foundation of China (NSFC); 201606040098 - Chinese Scholarship Council)Accepted manuscrip

Generating global products of LAI and FPAR from SNPP-VIIRS data: theoretical background and implementation

Leaf area index (LAI) and fraction of photosynthetically active radiation (FPAR) absorbed by vegetation have been successfully generated from the Moderate Resolution Imaging Spectroradiometer (MODIS) data since early 2000. As the Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard, the Suomi National Polar-orbiting Partnership (SNPP) has inherited the scientific role of MODIS, and the development of a continuous, consistent, and well-characterized VIIRS LAI/FPAR data set is critical to continue the MODIS time series. In this paper, we build the radiative transfer-based VIIRS-specific lookup tables by achieving minimal difference with the MODIS data set and maximal spatial coverage of retrievals from the main algorithm. The theory of spectral invariants provides the configurable physical parameters, i.e., single scattering albedos (SSAs) that are optimized for VIIRS-specific characteristics. The effort finds a set of smaller red-band SSA and larger near-infraredband SSA for VIIRS compared with the MODIS heritage. The VIIRS LAI/FPAR is evaluated through comparisons with one year of MODIS product in terms of both spatial and temporal patterns. Further validation efforts are still necessary to ensure the product quality. Current results, however, imbue confidence in the VIIRS data set and suggest that the efforts described here meet the goal of achieving the operationally consistent multisensor LAI/FPAR data sets. Moreover, the strategies of parametric adjustment and LAI/FPAR evaluation applied to SNPP-VIIRS can also be employed to the subsequent Joint Polar Satellite System VIIRS or other instruments.Accepted manuscrip

Evaluation of MODIS LAI/FPAR product Collection 6. Part 1: consistency and improvements

As the latest version of Moderate Resolution Imaging Spectroradiometer (MODIS) Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation (FPAR) products, Collection 6 (C6) has been distributed since August 2015. This collection is evaluated in this two-part series with the goal of assessing product accuracy, uncertainty and consistency with the previous version. In this first paper, we compare C6 (MOD15A2H) with Collection 5 (C5) to check for consistency and discuss the scale effects associated with changing spatial resolution between the two collections and benefits from improvements to algorithm inputs. Compared with C5, C6 benefits from two improved inputs: (1) L2G–lite surface reflectance at 500 m resolution in place of reflectance at 1 km resolution; and (2) new multi-year land-cover product at 500 m resolution in place of the 1 km static land-cover product. Global and seasonal comparison between C5 and C6 indicates good continuity and consistency for all biome types. Moreover, inter-annual LAI anomalies at the regional scale from C5 and C6 agree well. The proportion of main radiative transfer algorithm retrievals in C6 increased slightly in most biome types, notably including a 17% improvement in evergreen broadleaf forests. With same biome input, the mean RMSE of LAI and FPAR between C5 and C6 at global scale are 0.29 and 0.091, respectively, but biome type disagreement worsens the consistency (LAI: 0.39, FPAR: 0.102). By quantifying the impact of input changes, we find that the improvements of both land-cover and reflectance products improve LAI/FPAR products. Moreover, we find that spatial scale effects due to a resolution change from 1 km to 500 m do not cause any significant differences.Help from MODIS & VIIRS Science team members is gratefully acknowledged. This work is supported by the MODIS program of NASA and partially funded by the National Basic Research Program of China (Grant No. 2013CB733402), the key program of NSFC (Grant No. 41331171) and Chinese Scholarship Council. (MODIS program of NASA; 2013CB733402 - National Basic Research Program of China; 41331171 - NSFC; Chinese Scholarship Council

Prototyping of LAI and FPAR retrievals from MODIS multi-angle implementation of atmospheric correction (MAIAC) data

Leaf area index (LAI) and fraction of photosynthetically active radiation (FPAR) absorbed by vegetation are key variables in many global models of climate, hydrology, biogeochemistry, and ecology. These parameters are being operationally produced from Terra and Aqua MODIS bidirectional reflectance factor (BRF) data. The MODIS science team has developed, and plans to release, a new version of the BRF product using the multi-angle implementation of atmospheric correction (MAIAC) algorithm from Terra and Aqua MODIS observations. This paper presents analyses of LAI and FPAR retrievals generated with the MODIS LAI/FPAR operational algorithm using Terra MAIAC BRF data. Direct application of the operational algorithm to MAIAC BRF resulted in an underestimation of the MODIS Collection 6 (C6) LAI standard product by up to 10%. The difference was attributed to the disagreement between MAIAC and MODIS BRFs over the vegetation by −2% to +8% in the red spectral band, suggesting different accuracies in the BRF products. The operational LAI/FPAR algorithm was adjusted for uncertainties in the MAIAC BRF data. Its performance evaluated on a limited set of MAIAC BRF data from North and South America suggests an increase in spatial coverage of the best quality, high-precision LAI retrievals of up to 10%. Overall MAIAC LAI and FPAR are consistent with the standard C6 MODIS LAI/FPAR. The increase in spatial coverage of the best quality LAI retrievals resulted in a better agreement of MAIAC LAI with field data compared to the C6 LAI product, with the RMSE decreasing from 0.80 LAI units (C6) down to 0.67 (MAIAC) and the R2 increasing from 0.69 to 0.80. The slope (intercept) of the satellite-derived vs. field-measured LAI regression line has changed from 0.89 (0.39) to 0.97 (0.25).This work was funded by NASA Earth Science Division to MODIS (NNX14AI71G) and VIIRS (NNX14AP80A) programs through grants to Boston University (Ranga B. Myneni, PI), and HBO contract # 21205-14-036 to Yuri Knyazikhin. (NNX14AI71G - NASA; NNX14AP80A - NASA; 21205-14-036 - HBO contract)http://www.mdpi.com/2072-4292/9/4/370Published versio

New Approach for Temporal Stability Evaluation of Pseudo-Invariant Calibration Sites (PICS)

Pseudo-Invariant Calibration Sites (PICS) are one of the most popular methods for in-flight vicarious radiometric calibration of Earth remote sensing satellites. The fundamental question of PICS temporal stability has not been adequately addressed. However, the main purpose of this work is to evaluate the temporal stability of a few PICS using a new approach. The analysis was performed over six PICS (Libya 1, Libya 4, Niger 1, Niger 2, Egypt 1 and Sudan 1). The concept of a Virtual Constellation was developed to provide greater temporal coverage and also to overcome the dependence limitation of any specific characteristic derived from one particular sensor. TOA reflectance data from four sensors consistently demonstrating stable calibration to within 5%the Landsat 7 ETM+ (Enhanced Thematic Mapper Plus), Landsat 8 OLI (Operational Land Imager), Terra MODIS (Moderate Resolution Imaging Spectroradiometer) and Sentinel-2A MSI (Multispectral Instrument)were merged into a seamless dataset. Instead of using the traditional method of trend analysis (Students T test), a nonparametric Seasonal Mann-Kendall test was used for determining the PICS stability. The analysis results indicate that Libya 4 and Egypt 1 do not exhibit any monotonic trend in six reflective solar bands common to all of the studied sensors, indicating temporal stability. A decreasing monotonic trend was statistically detected in all bands, except SWIR 2, for Sudan 1 and the Green and Red bands for Niger 1. An increasing trend was detected in the Blue band for Niger 2 and the NIR band for Libya 1. These results do not suggest abandoning PICS as a viable calibration source. Rather, they indicate that PICS temporal stability cannot be assumed and should be regularly monitored as part of the sensor calibration process

Assessment of two aerosol optical thickness retrieval algorithms applied to MODIS aqua and terra measurements in Europe

© Author(s) 2012. This work is distributed under the Creative Commons Attribution 3.0 LicenseThe aim of the present study is to validate AOT (aerosol optical thickness) and A° ngström exponent (α), obtained from MODIS (MODerate resolution Imaging Spectroradiometer) Aqua and Terra calibrated level 1 data (1 km horizontal resolution at ground) with the SAER (Satellite AErosol Retrieval) algorithm and with MODIS Collection 5 (c005) standard product retrievals (10 km horizontal resolution), against AERONET (AErosol RObotic NETwork) sun photometer observations over land surfaces in Europe. An inter-comparison of AOT at 0.469 nm obtained with the two algorithms has also been performed. The time periods investigated were chosen to enable a validation of the findings of the two algorithms for a maximal possible variation in sun elevation. The satellite retrievals were also performed with a significant variation in the satellite-viewing geometry, since Aqua and Terra passed the investigation area twice a day for several of the cases analyzed. The validation with AERONET shows that the AOT at 0.469 and 0.555 nm obtained with MODIS c005 is within the expected uncertainty of one standard deviation of the MODIS c005 retrievals (1AOT =±0.05±0.15 ·AOT). The AOT at 0.443 nm retrieved with SAER, but with a much finer spatial resolution, also agreed reasonably well with AERONET measurements. The majority of the SAER AOT values are within the MODIS c005 expected uncertainty range, although somewhat larger average absolute deviation occurs compared to the results obtained with the MODIS c005 algorithm. The discrepancy between AOT from SAER and AERONET is, however, substantially larger for the wavelength 488 nm. This means that the values are, to a larger extent, outside of the expected MODIS uncertainty range. In addition, both satellite retrieval algorithms are unable to estimate accurately, although the MODIS c005 algorithm performs better. Based on the inter-comparison of the SAER and MODIS c005 algorithms, it was found that SAER on the whole is able to obtain results within the expected uncertainty range of MODIS Aqua and Terra observations.Peer reviewe

Earth observations from DSCOVR EPIC instrument

The National Oceanic and Atmospheric Administration (NOAA) Deep Space Climate Observatory (DSCOVR) spacecraft was launched on 11 February 2015 and in June 2015 achieved its orbit at the first Lagrange point (L1), 1.5 million km from Earth toward the sun. There are two National Aeronautics and Space Administration (NASA) Earth-observing instruments on board: the Earth Polychromatic Imaging Camera (EPIC) and the National Institute of Standards and Technology Advanced Radiometer (NISTAR). The purpose of this paper is to describe various capabilities of the DSCOVR EPIC instrument. EPIC views the entire sunlit Earth from sunrise to sunset at the backscattering direction (scattering angles between 168.5° and 175.5°) with 10 narrowband filters: 317, 325, 340, 388, 443, 552, 680, 688, 764, and 779 nm. We discuss a number of preprocessing steps necessary for EPIC calibration including the geolocation algorithm and the radiometric calibration for each wavelength channel in terms of EPIC counts per second for conversion to reflectance units. The principal EPIC products are total ozone (O3) amount, scene reflectivity, erythemal irradiance, ultraviolet (UV) aerosol properties, sulfur dioxide (SO2) for volcanic eruptions, surface spectral reflectance, vegetation properties, and cloud products including cloud height. Finally, we describe the observation of horizontally oriented ice crystals in clouds and the unexpected use of the O2 B-band absorption for vegetation properties.The NASA GSFC DSCOVR project is funded by NASA Earth Science Division. We gratefully acknowledge the work by S. Taylor and B. Fisher for help with the SO2 retrievals and Marshall Sutton, Carl Hostetter, and the EPIC NISTAR project for help with EPIC data. We also would like to thank the EPIC Cloud Algorithm team, especially Dr. Gala Wind, for the contribution to the EPIC cloud products. (NASA Earth Science Division)Accepted manuscrip

Trends in MODIS and AERONET derived aerosol optical thickness over Northern Europe

Long-term Aqua and Terra MODIS (MODerate resolution Imaging Spectroradiometer) Collections 5.1 and 6.1 (c051 and c061, respectively) aerosol data have been combined with AERONET (AERosol RObotic NETwork) ground-based sun photometer observations to examine trends in aerosol optical thickness (AOT, at 550 nm) over Northern Europe for the months April to September. For the 1927 and 1559 daily coincident measurements that were obtained for c051 and c061, respectively, MODIS AOT varied by 86 and 90%, respectively, within the predicted uncertainty of one standard deviation of the retrieval over land (ΔAOT = ±0.05 ± 0.15·AOT). For the coastal AERONET site Gustav Dalen Tower (GDT), Sweden, larger deviations were found for MODIS c051 and c061 (79% and 75%, respectively, within predicted uncertainty). The Baltic Sea provides substantially better statistical representation of AOT than the surrounding land areas and therefore favours the investigations of trends in AOT over the region. Negative trends of 1.5% and 1.2% per year in AOT, based on daily averaging, were found for the southwestern Baltic Sea from MODIS c051 and c061, respectively. This is in line with a decrease of 1.2% per year in AOT at the AERONET station Hamburg. For the western Gotland Basin area, Sweden, negative trends of 1.5%, 1.1% and 1.6% per year in AOT have been found for MODIS c051, MODIS c061 and AERONET GDT, respectively. The strongest trend of –1.8% per year in AOT was found for AERONET Belsk, Poland, which can be compared to –1.5% per day obtained from MODIS c051 over central Poland. The trends in MODIS and AERONET AOT are nearly all statistically significant at the 95% confidence level. The strongest aerosol sources are suggested to be located southwest, south and southeast of the investigation area, although the highest prevalence of pollution events is associated with air mass transport from southwest.Peer reviewe

Evaluation of MODIS LAI/FPAR product Collection 6. Part 2: Validation and intercomparison

The aim of this paper is to assess the latest version of the MODIS LAI/FPAR product (MOD15A2H), namely Collection 6 (C6). We comprehensively evaluate this product through three approaches: validation with field measurements, intercomparison with other LAI/FPAR products and comparison with climate variables. Comparisons between ground measurements and C6, as well as C5 LAI/FPAR indicate: (1) MODIS LAI is closer to true LAI than effective LAI; (2) the C6 product is considerably better than C5 with RMSE decreasing from 0.80 down to 0.66; (3) both C5 and C6 products overestimate FPAR over sparsely-vegetated areas. Intercomparisons with three existing global LAI/FPAR products (GLASS, CYCLOPES and GEOV1) are carried out at site, continental and global scales. MODIS and GLASS (CYCLOPES and GEOV1) agree better with each other. This is expected because the surface reflectances, from which these products were derived, were obtained from the same instrument. Considering all biome types, the RMSE of LAI (FPAR) derived from any two products ranges between 0.36 (0.05) and 0.56 (0.09). Temporal comparisons over seven sites for the 2001–2004 period indicate that all products properly capture the seasonality in different biomes, except evergreen broadleaf forests, where infrequent observations due to cloud contamination induce unrealistic variations. Thirteen years of C6 LAI, temperature and precipitation time series data are used to assess the degree of correspondence between their variations. The statistically-significant associations between C6 LAI and climate variables indicate that C6 LAI has the potential to provide reliable biophysical information about the land surface when diagnosing climate-driven vegetation responses.Help from MODIS and VIIRS Science team members is gratefully acknowledged. This work is supported by the MODIS program of NASA and partially funded by the National Basic Research Program of China (Grant No. 2013CB733402) and the key program of NSFC (Grant No. 41331171). Kai Yan gives thanks for the scholarship from the China Scholarship Council. (MODIS program of NASA; 2013CB733402 - National Basic Research Program of China; 41331171 - NSFC; China Scholarship Council