Comparison of Gross Primary Productivity Derived from GIMMS NDVI3g, GIMMS, and MODIS in Southeast Asia

Gross primary production (GPP) plays an important role in the net ecosystem exchange of CO2 between the atmosphere and terrestrial ecosystems. It is particularly important to monitor GPP in Southeast Asia because of increasing rates of tropical forest degradation and deforestation in the region in recent decades. The newly available, improved, third generation Normalized Difference Vegetation Index (NDVI3g) from the Global Inventory Modelling and Mapping Studies (GIMMS) group provides a long temporal dataset, from July 1981 to December 2011, for terrestrial carbon cycle and climate response research. However, GIMMS NDVI3g-based GPP estimates are not yet available. We applied the GLOPEM-CEVSA model, which integrates an ecosystem process model and a production efficiency model, to estimate GPP in Southeast Asia based on three independent results of the fraction of photosynthetically active radiation absorbed by vegetation (FPAR) from GIMMS NDVI3g (GPPNDVI3g), GIMMS NDVI1g (GPPNDVI1g), and the Moderate Resolution Imaging Spectroradiometer (MODIS) MOD15A2 FPAR product (GPPMOD15). The GPP results were validated using ground data from eddy flux towers located in different forest biomes, and comparisons were made among the three GPPs as well as the MOD17A2 GPP products (GPPMOD17). Based on validation with flux tower derived GPP estimates the results show that GPPNDVI3g is more accurate than GPPNDVI1g and is comparable in accuracy with GPPMOD15. In addition, GPPNDVI3g and GPPMOD15 have good spatial-temporal consistency. Our results indicate that GIMMS NDVI3g is an effective dataset for regional GPP simulation in Southeast Asia, capable of accurately tracking the variation and trends in long-term terrestrial ecosystem GPP dynamics

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

Study on Regional Responses of Pan-Arctic Terrestrial Ecosystems to Recent Climate Variability Using Satellite Remote Sensing

I applied a satellite remote sensing based production efficiency model (PEM) using an integrated AVHRR and MODIS FPAR/LAI time series with a regionally corrected NCEP/NCAR reanalysis surface meteorology and NASA/GEWEX shortwave solar radiation inputs to assess annual terrestrial net primary productivity (NPP) for the pan-Arctic basin and Alaska from 1983 to 2005. I developed a satellite remote sensing based evapotranspiration (ET) algorithm using GIMMS NDVI with the above meteorology inputs to assess spatial patterns and temporal trends in ET over the pan-Arctic region. I then analyzed associated changes in the regional water balance defined as the difference between precipitation (P) and ET. I finally analyzed the effects of regional climate oscillations on vegetation productivity and the regional water balance. The results show that low temperature constraints on Boreal-Arctic NPP are decreasing by 0.43% per year ( P \u3c 0.001), whereas a positive trend in vegetation moisture constraints of 0.49% per year ( P = 0.04) are offsetting the potential benefits of longer growing seasons and contributing to recent drought related disturbances in NPP. The PEM simulations of NPP seasonality, annual anomalies and trends are similar to stand inventory network measurements of boreal aspen stem growth ( r = 0.56; P = 0.007) and atmospheric CO2 measurement based estimates of the timing of growing season onset (r = 0.78; P \u3c 0.001). The simulated monthly ET results agree well (RMSE = 8.3 mm month-1; R2 = 0.89) with tower measurements for regionally dominant land cover types. Generally positive trends in ET, precipitation and available river discharge measurements imply that the pan-Arctic terrestrial water cycle is intensifying. Increasing water deficits occurred in some boreal and temperate grassland regions, which agree with regional drought records and recent satellite observations of vegetation browning and productivity decreases. Climate oscillations including Arctic Oscillation and Pacific Decadal Oscillation influence NPP by regulating seasonal patterns of low temperature and moisture constraints to photosynthesis. The pan-Arctic water balance is changing in complex ways in response to climate change and variability, with direct linkages to terrestrial carbon and energy cycles. Consequently, drought induced NPP decreases may become more frequent and widespread, though the occurrence and severity of drought events will depend on future water cycle patterns

QUANTIFICATION OF ERROR IN AVHRR NDVI DATA

Several influential Earth system science studies in the last three decades were based on Normalized Difference Vegetation Index (NDVI) data from Advanced Very High Resolution Radiometer (AVHRR) series of instruments. Although AVHRR NDVI data are known to have significant uncertainties resulting from incomplete atmospheric correction, orbital drift, sensor degradation, etc., none of these studies account for them. This is primarily because of unavailability of comprehensive and location-specific quantitative uncertainty estimates. The first part of this dissertation investigated the extent of uncertainty due to inadequate atmospheric correction in the widely used AVHRR NDVI datasets. This was accomplished by comparison with atmospherically corrected AVHRR data at AErosol RObotic NETwork (AERONET) sunphotometer sites in 1999. Of the datasets included in this study, Long Term Data Record (LTDR) was found to have least errors (precision=0.02 to 0.037 for clear and average atmospheric conditions) followed by Pathfinder AVHRR Land (PAL) (precision=0.0606 to 0.0418), and Top of Atmosphere (TOA) (precision=0.0613 to 0.0684). ` Although the use of field data is the most direct type of validation and is used extensively by the remote sensing community, it results in a single uncertainty estimate and does not account for spatial heterogeneity and the impact of spatial and temporal aggregation. These shortcomings were addressed by using Moderate Resolution Imaging Spectrometer (MODIS) data to estimate uncertainty in AVHRR NDVI data. However, before AVHRR data could be compared with MODIS data, the nonstationarity introduced by inter-annual variations in AVHRR NDVI data due to orbital drift had to be removed. This was accomplished by using a Bidirectional Reflectance Distribution Function (BRDF) correction technique originally developed for MODIS data. The results from the evaluation of AVHRR data using MODIS showed that in many regions minimal spatial aggregation will improve the precision of AVHRR NDVI data significantly. However temporal aggregation improved the precision of the data to a limited extent only. The research presented in this dissertation indicated that the NDVI change of ~0.03 to ~0.08 NDVI units in 10 to 20 years, frequently reported in recent literature, can be significant in some cases. However, unless spatially explicit uncertainty metrics are quantified for the specific spatiotemporal aggregation schemes used by these studies, the significance of observed differences between sites and temporal trends in NDVI will remain unknown

The 2010 Russian Drought Impact on Satellite Measurements of Solar-Induced Chlorophyll Fluorescence: Insights from Modeling and Comparisons with the Normalized Differential Vegetation Index (NDVI)

We examine satellite-based measurements of chlorophyll solar-induced fluorescence (SIF) over the region impacted by the Russian drought and heat wave of 2010. Like the popular Normalized Difference Vegetation Index (NDVI) that has been used for decades to measure photosynthetic capacity, SIF measurements are sensitive to the fraction of absorbed photosynthetically-active radiation (fPAR). However, in addition, SIF is sensitive to the fluorescence yield that is related to the photosynthetic yield. Both SIF and NDVI from satellite data show drought-related declines early in the growing season in 2010 as compared to other years between 2007 and 2013 for areas dominated by crops and grasslands. This suggests an early manifestation of the dry conditions on fPAR. We also simulated SIF using a global land surface model driven by observation-based meteorological fields. The model provides a reasonable simulation of the drought and heat impacts on SIF in terms of the timing and spatial extents of anomalies, but there are some differences between modeled and observed SIF. The model may potentially be improved through data assimilation or parameter estimation using satellite observations of SIF (as well as NDVI). The model simulations also offer the opportunity to examine separately the different components of the SIF signal and relationships with Gross Primary Productivity (GPP)

Hydroclimatic Controls on the Means and Variability of Vegetation Phenology and Carbon Uptake

Long-term, global offline (land-only) simulations with a dynamic vegetation phenology model are used to examine the control of hydroclimate over vegetation-related quantities. First, with a control simulation, the model is shown to capture successfully (though with some bias) key observed relationships between hydroclimate and the spatial and temporal variations of phenological expression. In subsequent simulations, the model shows that: (i) the global spatial variation of seasonal phenological maxima is controlled mostly by hydroclimate, irrespective of distributions in vegetation type, (ii) the occurrence of high interannual moisture-related phenological variability in grassland areas is determined by hydroclimate rather than by the specific properties of grassland, and (iii) hydroclimatic means and variability have a corresponding impact on the spatial and temporal distributions of gross primary productivity (GPP)

Evaluating and Quantifying the Climate-Driven Interannual Variability in Global Inventory Modeling and Mapping Studies (GIMMS) Normalized Difference Vegetation Index (NDVI3g) at Global Scales

Satellite observations of surface reflected solar radiation contain informationabout variability in the absorption of solar radiation by vegetation. Understanding thecauses of variability is important for models that use these data to drive land surface fluxesor for benchmarking prognostic vegetation models. Here we evaluated the interannualvariability in the new 30.5-year long global satellite-derived surface reflectance index data,Global Inventory Modeling and Mapping Studies normalized difference vegetation index(GIMMS NDVI3g). Pearsons correlation and multiple linear stepwise regression analyseswere applied to quantify the NDVI interannual variability driven by climate anomalies, andto evaluate the effects of potential interference (snow, aerosols and clouds) on the NDVIsignal. We found ecologically plausible strong controls on NDVI variability by antecedent precipitation and current monthly temperature with distinct spatial patterns. Precipitation correlations were strongest for temperate to tropical water limited herbaceous systemswhere in some regions and seasons 40 of the NDVI variance could be explained byprecipitation anomalies. Temperature correlations were strongest in northern mid- to-high-latitudes in the spring and early summer where up to 70 of the NDVI variance was explained by temperature anomalies. We find that, in western and central North America,winter-spring precipitation determines early summer growth while more recent precipitation controls NDVI variability in late summer. In contrast, current or prior wetseason precipitation anomalies were correlated with all months of NDVI in sub-tropical herbaceous vegetation. Snow, aerosols and clouds as well as unexplained phenomena still account for part of the NDVI variance despite corrections. Nevertheless, this study demonstrates that GIMMS NDVI3g represents real responses of vegetation to climate variability that are useful for global models

Toward a better understanding of changes in Northern vegetation using long-term remote sensing data

Cascading consequences of recent changes in the physical environment of northern lands associated with rapid warming have affected a broad range of ecosystem processes, particularly, changes in structure, composition, and functioning of vegetation. Incomplete understanding of underlying processes driving such changes is the primary motivation for this research. We report here the results of three studies that use long-term remote sensing data to advance our knowledge of spatiotemporal changes in growing season, greenness and productivity of northern vegetation. First, we improve the remote sensing-based detection of growing season by fusing vegetation greenness, snow and soil freeze/thaw condition. The satellite record reveals extensive lengthening trends of growing season and enhanced annual total greenness during the last three decades. Regionally varying seasonal responses are linked to local climate constraints and their relaxation. Second, we incorporate available land surface histories including disturbances and human land management practices to understand changes in remotely sensed vegetation greenness. This investigation indicates that multiple drivers including natural (wildfire) and anthropogenic (harvesting) disturbances, changing climate and agricultural activities govern the large-scale greening trends in northern lands. The timing and type of disturbances are important to fully comprehend spatially uneven vegetation changes in the boreal and temperate regions. In the final part, we question how photosynthetic seasonality evolved into its current state, and what role climatic constraints and their variability played in this process and ultimately in the carbon cycle. We take the ‘laws of minimum’ as a basis and introduce a new framework where the timing of peak photosynthetic activity (DOYPmax) acts as a proxy for plants adaptive state to climatic constraints on their growth. The result shows a widespread warming-induced advance in DOYPmax with an increase of total gross primary productivity across northern lands, which leads to an earlier phase shift in land-atmosphere carbon fluxes and an increase in their amplitude. The research presented in this dissertation suggests that understanding past, present and likely future changes in northern vegetation requires a multitude of approaches that consider linked climatic, social and ecological drivers and processes

Analysis of monotonic greening and browning trends from global NDVI time-series

Remotely sensed vegetation indices are widely used to detect greening and browning trends; especially the global coverage of time-series normalized difference vegetation index (NDVI) data which are available from 1981. Seasonality and serial auto-correlation in the data have previously been dealt with by integrating the data to annual values; as an alternative to reducing the temporal resolution, we apply harmonic analyses and non-parametric trend tests to the GIMMS NDVI dataset (1981-2006). Using the complete dataset, greening and browning trends were analyzed using a linear model corrected for seasonality by subtracting the seasonal component, and a seasonal non-parametric model. In a third approach, phenological shift and variation in length of growing season were accounted for by analyzing the time-series using vegetation development stages rather than calendar days. Results differed substantially between the models, even though the input data were the same. Prominent regional greening trends identified by several other studies were confirmed but the models were inconsistent in areas with weak trends. The linear model using data corrected for seasonality showed similar trend slopes to those described in previous work using linear models on yearly mean values. The non-parametric models demonstrated the significant influence of variations in phenology; accounting for these variations should yield more robust trend analyses and better understanding of vegetation trends

Reconstruction of a Long-term spatially Contiguous Solar-Induced Fluorescence (LCSIF) over 1982-2022

Satellite-observed solar-induced chlorophyll fluorescence (SIF) is a powerful proxy for diagnosing the photosynthetic characteristics of terrestrial ecosystems. Despite the increasing spatial and temporal resolutions of these satellite retrievals, records of SIF are primarily limited to the recent decade, impeding their application in detecting long-term dynamics of ecosystem function and structure. In this study, we leverage the two surface reflectance bands (red and near-infrared) available both from Advanced Very High-Resolution Radiometer (AVHRR, 1982-2022) and MODerate-resolution Imaging Spectroradiometer (MODIS, 2001-2022). Importantly, we calibrate and orbit-correct the AVHRR bands against their MODIS counterparts during their overlapping period. Using the long-term bias-corrected reflectance data, a neural network is then built to reproduce the Orbiting Carbon Observatory-2 SIF using AVHRR and MODIS, and used to map SIF globally over the entire 1982-2022 period. Compared with the previous MODIS-based CSIF product relying on four reflectance bands, our two-band-based product has similar skill but can be advantageously extended to the bias-corrected AVHRR period. Further comparison with three widely used vegetation indices (NDVI, kNDVI, NIRv; all based empirically on red and near-infrared bands) shows a higher or comparable correlation of LCSIF with satellite SIF and site-level GPP estimates across vegetation types, ensuring a greater capacity of LCSIF for representing terrestrial photosynthesis. Globally, LCSIF-AVHRR shows an accelerating upward trend since 1982, with an average rate of 0.0025 mW m-2 nm-1 sr-1 per decade during 1982-2000 and 0.0038 mW m-2 nm-1 sr-1 per decade during 2001-2022. Our LCSIF data provide opportunities to better understand the long-term dynamics of ecosystem photosynthesis and their underlying driving processes