Revealing climate and vegetation indices interactions through Cross Recurrence Techniques

Author(s):  
Andrés Felipe Almeida Ñauñay ◽  
Rosa María Benito Zafrilla ◽  
Miguel Quemada Sáenz-Badillos ◽  
Juan Carlos Losada ◽  
Ana María Tarquis Alfonso

<p>Grasslands are one of the world's major ecosystems groups many of them are now being acknowledged as having a multifunctional role in producing food and rehabilitating croplands, in environmental management and cultural heritage. Multiple studies showed pasture grasslands as a complex agroecological system, depending on multiple variables with a nonlinear dynamic greatly affected by climate fluctuations over time. Remote sensing methods proved to be one of the most effective techniques for monitoring variations over wide areas. In this line, vegetation indices (VIs) demonstrated to be an appropriate indicator of vegetation cover condition. This study aims to perform a method to visualize and quantify the complexity between semiarid grasslands and climate. With this goal, VIs and climate time series are analysed simultaneously with non-linear techniques to reveal the dynamic behaviour of both series over time and their interaction.</p><p>A semi-arid grassland area characterized by a Mediterranean climate with a continental character and low precipitation throughout the year were chosen. VIs time series were constructed from MODIS TERRA (MOD09Q1.006) product. Multispectral images composed by 8-days were acquired from 2002 till 2018 and four pixels with a spatial resolution of 250 x 250 m<sup>2</sup> were chosen in the selected area. Normalized Difference Vegetation Index (NDVI) and Modified Soil-Adjusted Vegetation Index (MSAVI) were calculated based on these images. Temperature and precipitation series were acquired from a near meteorological station and adapted to 8-day time step.</p><p>Cross-Recurrence plots (CRP) and recurrence quantification analysis (RQA) were performed to analyse the climate and vegetation dynamics simultaneously. To achieve this goal, several measures of complexity were computed, such as Determinism (DET), average diagonal length (LT) and entropy (ENT).</p><p>Our results showed different CRPs depending on the climate variable and the utilized VIs. Temperature and VIs CRPs showed a periodical pattern, implying the temperature seasonality over time. In contrast, precipitation and VIs CRPs showed more chaotical behaviour, due to the occurrence of extreme events and seasonal shifts. These results are quantified by the DET and ENTR values.</p><p>Our results indicate that temperature and precipitation present a dynamical complexity that is revealed in VIs response. Both indices showed different results of complexity measures, implying that MSAVI has a higher complexity than NDVI. This fact is probably due to the addition of a variable soil adjustment factor. Consequently, MSAVI should be considered as a potential alternative to NDVI in semiarid areas.</p><p><strong>Reference</strong></p><p>Almeida-Ñauñay, A. F., Benito, R. M., Quemada, M., Losada, J. C., & Tarquis, A. M. Complexity of the Vegetation-Climate System Through Data Analysis. In International Conference on Complex Networks and Their Applications. Springer, Cham., 609-619, 2020</p><p><strong>Acknowledgements</strong></p><p>The authors acknowledge support from Project No. PGC2018-093854-B-I00 of the Spanish Ministerio de Ciencia Innovación y Universidades of Spain and the funding from the Comunidad de Madrid (Spain), Structural Funds 2014-2020 512 (ERDF and ESF), through project AGRISOST-CM S2018/BAA-4330 and the financial support from Boosting Agricultural Insurance based on Earth Observation data - BEACON project under agreement Nº 821964, funded under H2020_EU, DT-SPACE-01-EO-2018-2020.</p>

2020 ◽  
Author(s):  
David Rivas-Tabares ◽  
Juan J. Martín-Sotoca ◽  
Antonio Saa-Requejo ◽  
Ana María Tarquis

<p>Crop yields of rainfed cereal are highly dependent of the soil-plant-atmosphere system, especially referred to the weather conditions and soil properties. The study of this interaction is feasible through the earth observations of historical data. Remote sensing data and agricultural survey work together identifying and analyzing plots with monocrop cereal sequences. In this research, we investigate the relation of the Normalized Difference Vegetation Index (NDVI) residual time series behavior relative to soil classes from Self-Organizing Maps (SOM) and the precipitation residual time series.</p><p>The midlands of Eresma-Adaja watershed (Dueros’ River basin, Spain) is historically depicted to rainfed cereal agriculture, some evidence of monocropping sequences are worrisome the water availability in the area. Within this area, two contrasting soil properties sites were selected to assess plots with at least 20 years of rainfed monocropping sequences but under similar weather regime. This allows analyzing the effect and relationships of this practice by soil type in time. For this, we treat the NDVI and precipitation time residual series as signals. The use of the Generalized Structure Function applied to these residual time series and the Hurst exponent, serve to confirm the soil properties differences from SOM and to reinforce the scaling properties of soil-climate interaction in semiarid regions for cereals in monocrop. As a result, the NDVI and precipitation series present an antipersistence behavior supporting that precipitation regime is influencing as the same manner the NDVI residual time series among complimentary factors.</p><p><strong>ACKNOWLEDGEMENTS</strong></p><p>Finding for this work was partially provided by Boosting agricultural Insurance based on Earth Observation data - BEACON project under agreement Nº 821964, funded under H2020_EU, DT-SPACE-01-EO-2018-2020. The authors also acknowledge support from Project No. PGC2018-093854-B-I00 of the Spanish <em>Ministerio de Ciencia Innovación y Universidades</em> of Spain. The data provided by ITACyL and AEMET is greatly appreciated.</p><p> </p>


2021 ◽  
Author(s):  
David Rivas-Tabares ◽  
Ana María Tarquis Alfonso

<p>Rainfed crops as cereals in the semiarid are common and extensive land cover in which climate, soils and atmosphere interact trough nonlinear relationships. Earth Observations coupled to ground monitoring network allow to improve the understanding of these relationships during each cropping season. However, novel analysis is required to understand these relationships in larger periods to improve sustainability and suitability of the productive areas in the semiarid.</p><p>The aim of this work is to use a joint multifractal approach using vegetation indices, precipitation, and temperatures to analyze atmosphere-plant nonlinear relationships. For this, time series of 20 cropping seasons were used to characterize these relationships in central Spain. The Generalized Structure Function and the derived Generalized Hurst Exponent analysis were implemented to investigate precipitation, vegetation indices and temperature time series. For this, an exhaustive selection based on land use and a land cover change analysis was performed to detect plots in which cereal crop sequences are dedicated to barley and wheat over the period 2000 to 2020.</p><p>As a result, two agro zones were characterized by different multifractal properties. Precipitation series show antipersistent characteristics and fractal properties between zones while original vegetation indices show trending behavior but shifted between analyzed zones. Nonetheless, soils and rainfall events can vary interannual conditions in which the crop is developing. For vegetation indices long-term series the trend is persistent. Even so, the dynamics of vegetation indices also provide more information when annual patterns are extracted from the series, exhibiting fractal properties mainly from rainfall pattern of each zone. Finally, in this case, the joint multifractal analysis served to characterize agro zones using earth observation and climate data for extensive cereals in Central Spain.</p><p><strong>Reference</strong></p><p>Rivas-Tabares D., Tarquis A.M. (2021) Towards Understanding Complex Interactions of Normalized Difference Vegetation Index Measurements Network and Precipitation Gauges of Cereal Growth System. In: Benito R.M., Cherifi C., Cherifi H., Moro E., Rocha L.M., Sales-Pardo M. (eds) Complex Networks & Their Applications IX. COMPLEX NETWORKS 2020 2020. Studies in Computational Intelligence, vol 943. Springer, Cham. https://doi.org/10.1007/978-3-030-65347-7_51</p><p><strong>Acknowledgements</strong></p><p>The authors acknowledge support from Project No. PGC2018-093854-B-I00 of the Spanish Ministerio de Ciencia Innovación y Universidades of Spain and the funding from the Comunidad de Madrid (Spain), Structural Funds 2014-2020 512 (ERDF and ESF), through project AGRISOST-CM S2018/BAA-4330 and the financial support from Boosting Agricultural Insurance based on Earth Observation data - BEACON project under agreement Nº 821964, funded under H2020_EU, DT-SPACE-01-EO-2018-2020.</p>


2011 ◽  
Vol 4 (4) ◽  
pp. 1103-1114 ◽  
Author(s):  
F. Maignan ◽  
F.-M. Bréon ◽  
F. Chevallier ◽  
N. Viovy ◽  
P. Ciais ◽  
...  

Abstract. Atmospheric CO2 drives most of the greenhouse effect increase. One major uncertainty on the future rate of increase of CO2 in the atmosphere is the impact of the anticipated climate change on the vegetation. Dynamic Global Vegetation Models (DGVM) are used to address this question. ORCHIDEE is such a DGVM that has proven useful for climate change studies. However, there is no objective and methodological way to accurately assess each new available version on the global scale. In this paper, we submit a methodological evaluation of ORCHIDEE by correlating satellite-derived Vegetation Index time series against those of the modeled Fraction of absorbed Photosynthetically Active Radiation (FPAR). A perfect correlation between the two is not expected, however an improvement of the model should lead to an increase of the overall performance. We detail two case studies in which model improvements are demonstrated, using our methodology. In the first one, a new phenology version in ORCHIDEE is shown to bring a significant impact on the simulated annual cycles, in particular for C3 Grasses and C3 Crops. In the second case study, we compare the simulations when using two different weather fields to drive ORCHIDEE. The ERA-Interim forcing leads to a better description of the FPAR interannual anomalies than the simulation forced by a mixed CRU-NCEP dataset. This work shows that long time series of satellite observations, despite their uncertainties, can identify weaknesses in global vegetation models, a necessary first step to improving them.


2019 ◽  
Vol 19 (8) ◽  
pp. 1685-1702 ◽  
Author(s):  
Juan José Martín-Sotoca ◽  
Antonio Saa-Requejo ◽  
Rubén Moratiel ◽  
Nicolas Dalezios ◽  
Ioannis Faraslis ◽  
...  

Abstract. Vegetation indices based on satellite images, such as the normalized difference vegetation index (NDVI), have been used in countries like the USA, Canada and Spain for damaged pasture and forage insurance over the last few years. This type of agricultural insurance is called satellite-index-based insurance (SIBI). In SIBI, the occurrence of damage is defined as normal distributions. In this work a pasture area at the north of the Community of Madrid (Spain) has been delimited by means of Moderate Resolution Imaging Spectroradiometer (MODIS) images. A statistical analysis of NDVI histograms was applied to seek for alternative distributions using the maximum likelihood method and χ2 test. The results show that the normal distribution is not the optimal representation and the generalized extreme value (GEV) distribution presents a better fit through the year based on a quality estimator. A comparison between normal and GEV is shown with respect to the probability under a NDVI threshold value throughout the year. This suggests that an a priori distribution should not be selected and a percentile methodology should be used to define a NDVI damage threshold rather than the average and standard deviation, typically of normal distributions. Highlights. The GEV distribution provides better fit to the NDVI historical observations than the normal one. Differences between normal and GEV distributions are higher during spring and autumn, which are transition periods in the precipitation regimen. NDVI damage threshold shows evident differences using normal and GEV distributions both covering the same probability (24.20 %). NDVI damage threshold values based on percentile calculation are proposed as an improvement in the index-based insurance in damaged pasture.


2020 ◽  
Vol 12 (6) ◽  
pp. 907 ◽  
Author(s):  
Teodoro Semeraro ◽  
Andrea Luvisi ◽  
Antonio O. Lillo ◽  
Roberta Aretano ◽  
Riccardo Buccolieri ◽  
...  

Forests are important in sequestering CO2 and therefore play a significant role in climate change. However, the CO2 cycle is conditioned by drought events that alter the rate of photosynthesis, which is the principal physiological action of plants in transforming CO2 into biological energy. This study applied recurrence quantification analysis (RQA) to describe the evolution of photosynthesis-related indices to highlight disturbance alterations produced by the Atlantic Multidecadal Oscillation (AMO, years 2005 and 2010) and the El Niño-Southern Oscillation (ENSO, year 2015) in the Amazon forest. The analysis was carried out using Moderate Resolution Imaging Spectroradiometer (MODIS) images to build time series of the enhanced vegetation index (EVI), the normalized difference water index (NDWI), and the land surface temperature (LST) covering the period 2001–2018. The results did not show significant variations produced by AMO throughout the study area, while a disruption due to the global warming phase linked to the extreme ENSO event occurred, and the forest was able to recover. In addition, spatial differences in the response of the forest to the ENSO event were found. These findings show that the application of RQA to the time series of vegetation indices supports the evaluation of the forest ecosystem response to disruptive events. This approach provides information on the capacity of the forest to recover after a disruptive event and, therefore is useful to estimate the resilience of this particular ecosystem.


Author(s):  
H. Bendini ◽  
I. D. Sanches ◽  
T. S. Körting ◽  
L. M. G. Fonseca ◽  
A. J. B. Luiz ◽  
...  

The objective of this research is to classify agricultural land use in a region of the Cerrado (Brazilian Savanna) biome using a time series of Enhanced Vegetation Index (EVI) from Landsat 8 OLI. Phenological metrics extracted from EVI time series, a Random Forest algorithm and data mining techniques are used in the process of classification. The area of study is a region in the Cerrado in a region of the municipality of Casa Branca, São Paulo state, Brazil. The results are encouraging and demonstrate the potential of phenological parameters obtained from time series of OLI vegetation indices for agricultural land use classification.


Author(s):  
D. Ratha ◽  
D. Mandal ◽  
S. Dey ◽  
A. Bhattacharya ◽  
A. Frery ◽  
...  

Abstract. In this paper, we present two radar vegetation indices for full-pol and compact-pol SAR data, respectively. Both are derived using the notion of a geodesic distance between observation and well-known scattering models available in the literature. While the full-pol version depends on a generalized volume scattering model, the compact-pol version uses the ideal depolariser to model the randomness in the vegetation. We have utilized the RADARSAT Constellation Mission (RCM) time-series data from the SAMPVEX16-MB campaign in the Manitoba region of Canada for comparing and assessing the indices in terms of the change in the biophysical parameters as well. The compact-pol data for comparison is simulated from the full-pol RCM time series. Both the indices show better performance at correlating with biophysical parameters such as Plant Area Index (PAI) and Volumetric Water Content (VWC) for wheat and soybean crops, in comparison to the state-of-art Radar Vegetation Index (RVI) of Kim and van Zyl. These indices are timely for the upcoming release of the data from the RCM, which will provide data in both full and compact-pol modes, aimed at better crop monitoring from space.


Author(s):  
E. Çolak ◽  
M. Chandra ◽  
F. Sunar

Abstract. Recently, the demand for nuclear power plants has been increasing in developing countries in line with global energy demands. Turkey, one of the developing economies, is also making plans for nuclear power generation since 1970. The Sinop Nuclear Power Plant was a planned nuclear plant located in the Turkey's most northern point in an area where 99% of the land is forest, in Sinop Peninsula. If disputes are resolved and its construction continues, the plant is expected to be put into service in 2028. On the other hand, due to the construction of the nuclear power plant, the land cover in and around the plant site has changed, potentially causing major environmental changes. As an example, more than 650000 trees have been cut down so far for the construction of a nuclear power plant, which may have a negative impact on the region's ecological balances by endangering biodiversity and causing ecological damage. The aim of this study is to detect changes in forest areas from the start of nuclear power plant construction through December 2020 using Sentinel 1 SAR and Sentinel 2 optical time series images. For this purpose, different radar and optical vegetation indices such as Modified Radar Vegetation Index (mRVI), Modified Radar Forest Degradation Index (mRFDI), and Normalized Difference Vegetation Index (NDVI) were applied using Google Earth Engine (GEE) Sentinel 1/2 satellite time series for 2015–2020 period. As a result, the indices used were found to yield findings consistent with the reported negative land cover change. In addition, correlation analysis were made between the radar vegetation indices used and a very high negative correlation (−0.99) was found. The annual distributions of the values of the three indices used were statistically evaluated using boxplots.


2021 ◽  
Vol 13 (20) ◽  
pp. 4126
Author(s):  
Yang Li ◽  
Ziti Jiao ◽  
Kaiguang Zhao ◽  
Yadong Dong ◽  
Yuyu Zhou ◽  
...  

Vegetation indices are widely used to derive land surface phenology (LSP). However, due to inconsistent illumination geometries, reflectance varies with solar zenith angles (SZA), which in turn affects the vegetation indices, and thus the derived LSP. To examine the SZA effect on LSP, the MODIS bidirectional reflectance distribution function (BRDF) product and a BRDF model were employed to derive LSPs under several constant SZAs (i.e., 0°, 15°, 30°, 45°, and 60°) in the Harvard Forest, Massachusetts, USA. The LSPs derived under varying SZAs from the MODIS nadir BRDF-adjusted reflectance (NBAR) and MODIS vegetation index products were used as baselines. The results show that with increasing SZA, NDVI increases but EVI decreases. The magnitude of SZA-induced NDVI/EVI changes suggests that EVI is more sensitive to varying SZAs than NDVI. NDVI and EVI are comparable in deriving the start of season (SOS), but EVI is more accurate when deriving the end of season (EOS). Specifically, NDVI/EVI-derived SOSs are relatively close to those derived from ground measurements, with an absolute mean difference of 8.01 days for NDVI-derived SOSs and 9.07 days for EVI-derived SOSs over ten years. However, a considerable lag exists for EOSs derived from vegetation indices, especially from the NDVI time series, with an absolute mean difference of 14.67 days relative to that derived from ground measurements. The SOSs derived from NDVI time series are generally earlier, while those from EVI time series are delayed. In contrast, the EOSs derived from NDVI time series are delayed; those derived from the simulated EVI time series under a fixed illumination geometry are also delayed, but those derived from the products with varying illumination geometries (i.e., MODIS NBAR product and MODIS vegetation index product) are advanced. LSPs derived from varying illumination geometries could lead to a difference spanning from a few days to a month in this case study, which highlights the importance of normalizing the illumination geometry when deriving LSP from NDVI/EVI time series.


2018 ◽  
Author(s):  
Henning Lebrenz ◽  
Andras Bárdossy

Abstract. The widely applied geostatistical interpolation methods of Ordinary Kriging (OK) or External Drift Kriging (EDK) interpolate the variable of interest to the unknown location, providing a linear estimator and an estimation variance as measure of uncertainty. The methods implicitly pose the assumption of Gaussianity on the observations, which is not given for many variables. The resulting best linear and unbiased estimator from the subsequent interpolation optimizes the mean error over many realizations for the entire spatial domain and, therefore, allows a systematic under- (over-) estimation of the variable in regions of relatively high (low) observations. In case of a variable with observed time-series, the spatial marginal distributions are estimated separately for one time step after the other, and the errors from the interpolations might accumulate over time in regions of relatively extreme observations. Therefore, we propose the interpolation method of Quantile Kriging (QK) with a two step procedure prior to interpolation: we firstly estimate distributions of the variable over time at the observation locations and then estimate the marginal distributions over space for every given time step. For this purpose, a distribution function is selected and fitted to the observed time-series at every observation location, thus converting the variable into quantiles and defining parameters. At a given time step, the quantiles from all observation locations are then transformed into a Gaussian-distributed variable by a twofold quantile-quantile transformation with the Beta- and the Normal-distribution function. The spatio-temporal description of the proposed method accommodates skewed marginal distributions and resolves the spatial non-stationarity of the original variable. The Gaussian-distributed variable and the distribution parameters are now interpolated by OK and EDK. At the unknown location, the resulting outcomes are reconverted back into the estimator and the estimation variance of the original variable. As a summary, QK newly incorporates information from the temporal axis for its spatial marginal distribution and subsequent interpolation and, therefore, could be interpreted as a space-time version of Probability Kriging. In this study, QK is applied for the variable of observed monthly precipitation from raingauges in South Africa. The estimators and estimation variances from the interpolation are compared to the respective outcomes from OK and EDK. The cross-validations shows that QK improves the estimator and the estimation variance for most of the selected objective functions. QK further enables the reduction of the temporal bias at locations of extreme observations. The performance of QK, however, declines when many zero-value observations are present in the input data. It is further revealed that QK relates the magnitude of its estimator with the magnitude of the respective estimation variance as opposed to the traditional methods of OK and EDK, whose estimation variances do only depend on the spatial configuration of the observation locations and the model settings.


Sign in / Sign up

Export Citation Format

Share Document