On the use of MODIS EVI to assess gross primary productivity of North American ecosystems

[1] Carbon flux models based on light use efficiency (LUE), such as the MOD17 algorithm, have proved difficult to parameterize because of uncertainties in the LUE term, which is usually estimated from meteorological variables available only at large spatial scales. In search of simpler models based entirely on remote‐sensing data, we examined direct relationships between the enhanced vegetation index (EVI) and gross primary productivity (GPP) measured at nine eddy covariance flux tower sites across North America. When data from the winter period of inactive photosynthesis were excluded, the overall relationship between EVI and tower GPP was better than that between MOD17 GPP and tower GPP. However, the EVI/GPP relationships vary between sites. Correlations between EVI and GPP were generally greater for deciduous than for evergreen sites. However, this correlation declined substantially only for sites with the smallest seasonal variation in EVI, suggesting that this relationship can be used for all but the most evergreen sites. Within sites dominated by either evergreen or deciduous species, seasonal variation in EVI was best explained by the severity of summer drought. Our results demonstrate that EVI alone can provide estimates of GPP that are as good as, if not better than, current versions of the MOD17 algorithm for many sites during the active period of photosynthesis. Preliminary data suggest that inclusion of other remote‐sensing products in addition to EVI, such as the MODIS land surface temperature (LST), may result in more robust models of carbon balance based entirely on remote‐sensing data

Vegetation Dynamics Revealed by Remote Sensing and Its Feedback to Regional and Global Climate

This book focuses on some significant progress in vegetation dynamics and their response to climate change revealed by remote sensing data. The development of satellite remote sensing and its derived products offer fantastic opportunities to investigate vegetation changes and their feedback to regional and global climate systems. Special attention is given in the book to vegetation changes and their drivers, the effects of extreme climate events on vegetation, land surface albedo associated with vegetation changes, plant fingerprints, and vegetation dynamics in climate modeling

Soil drought anomalies in MODIS GPP of a Mediterranean broadleaved evergreen forest

The Moderate Resolution Imaging Spectroradiometer (MODIS) yields global operational estimates of terrestrial gross primary production (GPP). In this study, we compared MOD17A2 GPP with tower eddy flux-based estimates of GPP from 2001 to 2010 over an evergreen broad-leaf Mediterranean forest in Southern France with a significant summer drought period. The MOD17A2 GPP shows seasonal variations that are inconsistent with the tower GPP, with close-to-accurate winter estimates and significant discrepancies for summer estimates which are the least accurate. The analysis indicated that the MOD17A2 GPP has high bias relative to tower GPP during severe summer drought which we hypothesized caused by soil water limitation. Our investigation showed that there was a significant correlation (R-2 = 0.77, p < 0.0001) between the relative soil water content and the relative error of MOD17A2 GPP. Therefore, the relationship between the error and the measured relative soil water content could explain anomalies in MOD17A2 GPP. The results of this study indicate that careful consideration of the water conditions input to the MOD17A2 GPP algorithm on remote sensing is required in order to provide accurate predictions of GPP. Still, continued efforts are necessary to ascertain the most appropriate index, which characterizes soil water limitation in water-limited environments using remote sensing

Estimating causal networks in biosphere–atmosphere interaction with the PCMCI approach

Local meteorological conditions and biospheric activity are tightly coupled. Understanding these links is an essential prerequisite for predicting the Earth system under climate change conditions. However, many empirical studies on the interaction between the biosphere and the atmosphere are based on correlative approaches that are not able to deduce causal paths, and only very few studies apply causal discovery methods. Here, we use a recently proposed causal graph discovery algorithm, which aims to reconstruct the causal dependency structure underlying a set of time series. We explore the potential of this method to infer temporal dependencies in biosphere-atmosphere interactions. Specifically we address the following questions: How do periodicity and heteroscedasticity influence causal detection rates, i.e. the detection of existing and non-existing links? How consistent are results for noise-contaminated data? Do results exhibit an increased information content that justifies the use of this causal-inference method? We explore the first question using artificial time series with well known dependencies that mimic real-world biosphere-atmosphere interactions. The two remaining questions are addressed jointly in two case studies utilizing observational data. Firstly, we analyse three replicated eddy covariance datasets from a Mediterranean ecosystem at half hourly time resolution allowing us to understand the impact of measurement uncertainties. Secondly, we analyse global NDVI time series (GIMMS 3g) along with gridded climate data to study large-scale climatic drivers of vegetation greenness. Overall, the results confirm the capacity of the causal discovery method to extract time-lagged linear dependencies under realistic settings. The violation of the method's assumptions increases the likelihood to detect false links. Nevertheless, we consistently identify interaction patterns in observational data. Our findings suggest that estimating a directed biosphere-atmosphere network at the ecosystem level can offer novel possibilities to unravel complex multi-directional interactions. Other than classical correlative approaches, our findings are constrained to a few meaningful set of relations which can be powerful insights for the evaluation of terrestrial ecosystem models

An analysis of long-term effects of climate change and socioeconomic activities on grassland productivity of inner Mongolia

In recent years, researchers have recognized the complexity of the interactions between the ecological system and the economic development of human society. However, the complicated relationships overwhelm traditional statistical procedures and require an innovative approach to investigate their dynamics. We proposed this study to provide a unique perspective in analyzing the long-term causal relationships between the grassland productivity, climate change, and socioeconomic development of Inner Mongolia Autonomous Region (IMAR) of China. Our attempt began with acquiring remotely sensed satellite imagery, climatic variations, and aggregated annual reports of the socio-economy of the IMAR in vegetation growing seasons for 15 years. The spatial and temporal dissimilarities of the raw observations prevented us from exploiting the potential of this valuable dataset; thus, we interpolated and extrapolated the data to generate a panel dataset with consistent spatial and temporal resolutions. Then, we took another step to preprocess the panel data by applying a signal filter to isolate the long-term trend of change from the inter- and intra-annual cyclic patterns and used the trends as the input for a panel data model. The results from our statistical analysis indicated that the independent variables explained the variations in the dependent variable extremely well, while the polynomial terms of climatic variables were significant with limited marginal effect and most of the climatic variables showed negative linear impact on the grassland productivity. In the meantime, we found not all socioeconomic variables we attempted to include into the model significantly affected grassland productivity, especially the variables describing the financial status of the IMAR residents

Impacts of Climate Extremes on Terrestrial Productivity

Terrestrial biosphere absorbs approximately 28% of anthropogenic CO2 emissions. This terrestrial carbon sink might become saturated in a future climate regime. To explore the issues associated with this topic, an accurate estimate of gross primary production (GPP) of global terrestrial ecosystems is needed. A major uncertainty in modeling global terrestrial GPP is the parameter of light use efficiency (LUE). Most LUE estimates in global models are satellite-based and coarsely measured with emphasis on environmental variables. Others are from eddy covariance towers with much greater spatial and temporal data quality and emphasis on mechanistic processes, but in a limited number of sites. In this study, we conducted a comprehensive global study of tower-based LUE from 237 FLUXNET towers, and scaled up LUEs from in-situ tower level to global biome level. We integrated the tower-based LUE estimates with key environmental and biological variables at 0.5º × 0.5º grid-cell resolutions, using a random forest regression (RFR) approach. Then we developed a RFR-LUE-GPP model using the grid-cell LUE data. In order to calibrate the LUE model, we developed a data-driven RFR-GPP model using random forest regression method only. Our results showed LUE varies largely with latitude. We estimated a global area-weighted average of LUE at 1.23±0.03 gC m-2 MJ-1 APAR, which led to an estimate of global gross primary production (GPP) of 107.5±2.5 Gt C /year from 2001 to 2005. Large uncertainties existed in GPP estimations over sparsely vegetated areas covered by savannas and woody savannas at middle to low latitude (i.e. 20ºS to 40ºS and 5ºN to 40ºN) due to the lack of available data. Model results were improved by incorporating Köppen climate types to represent climate/meteorological information in machine learning modeling. This brought a new understanding to the recognized problem of climate-dependence of spring onset of photosynthesis and the challenges in accurately modeling the biome GPP of evergreen broad leaf forests (EBF). The divergent responses of GPP to temperature and precipitation at mid-high latitudes and at mid-low latitudes echo the necessity of modeling GPP separately by latitudes. We also used a perfect-deficit approach to identify forest canopy photosynthetic capacity (CPC) deficits and analyze how they correlate to climate extremes, based on observational data measured by the eddy covariance method at 27 forest sites over 146 site-years. We found that droughts severely affect the carbon assimilation capacities of evergreen broadleaf forest and deciduous broadleaf forest. The carbon assimilation capacities of Mediterranean forests were highly sensitive to climate extremes, while marine forest climates tended to be insensitive to climate extremes. Our estimates suggest an average global reduction of forest canopy photosynthetic capacity due to unfavorable climate extremes of 6.3 Pg C (~5.2% of global gross primary production) per growing season over 2001-2010, with evergreen broadleaf forests contributing 52% of the total reduction. At biome-scale, terrestrial carbon uptake is controlled mainly by weather variability. Observational data from a global monitoring network indicate that the sensitivity of terrestrial carbon sequestration to mean annual temperature (T) breaks down at a threshold value of 16oC, above which terrestrial CO2 fluxes are controlled by dryness rather than temperature. Here we show that since 1948 warming climate has moved the 16oC T latitudinal belt poleward. Land surface area with T \u3e16oC and now subject to dryness control rather than temperature as the regulator of carbon uptake has increased by 6% and is expected to increase by at least another 8% by 2050

Intermediate snowpack melt-out dates guarantee the highest seasonal grasslands greening in the Pyrenees

In mountain areas, the phenology and productivity of grassland are closely related to snow dynamics. However, the influence that snow melt timing has on grassland growing still needs further attention for a full understanding, particularly at high spatial resolution. Aiming to reduce this knowledge gap, this work exploits 1 m resolution snow depth and Normalized Difference Vegetation Index observations acquired with an Unmanned Aerial Vehicle at a sub-alpine site in the Pyrenees. During two snow seasons (2019–2020 and 2020–2021), 14 NDVI and 17 snow depth distributions were acquired over 48 ha. Despite the snow dynamics being different in the two seasons, the response of grasslands greening to snow melt-out exhibited a very similar pattern in both. The NDVI temporal evolution in areas with distinct melt-out dates reveals that sectors where the melt-out date occurs in late April or early May (optimum melt-out) reach the maximum vegetation productivity. Zones with an earlier or a later melt-out rarely reach peak NDVI values. The results obtained in this study area, suggest that knowledge about snow depth distribution is not needed to understand NDVI grassland dynamics. The analysis did not reveal a clear link between the spatial variability in snow duration and the diversity and richness of grassland communities within the study area

The uncertainty analysis of the MODIS GPP product in global maize croplands

Gross primary productivity (GPP) is very important in the global carbon cycle. Currently, the newly released estimates of 8-day GPP at 500 m spatial resolution (Collection 6) are provided by the Moderate Resolution Imaging Spectroradiometer (MODIS) Land Science Team for the global land surface via the improved light use efficiency (LUE) model. However, few studies have evaluated its performance. In this study, the MODIS GPP products (GPPMOD) were compared with the observed GPP (GPPEC) values from site-level eddy covariance measurements over seven maize flux sites in different areas around the world. The results indicate that the annual GPPMOD was underestimated by 6%‒58% across sites. Nevertheless, after incorporating the parameters of the calibrated LUE, the measurements of meteorological variables and the reconstructed Fractional Photosynthetic Active Radiation (FPAR) into the GPPMOD algorithm in steps, the accuracies of GPPMOD estimates were improved greatly, albeit to varying degrees. The differences between the GPPMOD and the GPPEC were primarily due to the magnitude of LUE and FPAR. The underestimate of maize cropland LUE was a widespread problem which exerted the largest impact on the GPPMOD algorithm. In American and European sites, the performance of the FPAR exhibited distinct differences in capturing vegetation GPP during the growing season due to the canopy heterogeneity. In addition, at the DE-Kli site, the GPPMOD abruptly produced extreme low values during the growing season because of the contaminated FPAR from a continuous rainy season. After correcting the noise of the FPAR, the accuracy of the GPPMOD was improved by approximately 14%. Therefore, it is crucial to further improve the accuracy of global GPPMOD, especially for the maize crop ecosystem, to maintain food security and better understand global carbon cycle