1. Introduction
1.1 The relationship between normalized difference vegetation index and carbon dioxide
The use of atmospheric carbon dioxide (CO2) seasonal index measures its sources and sinks with high precision. Understanding the seasonal dynamics of these emissions is necessary to enable an effective global mitigation. However, it is difficult to estimate carbon-driven emissions using only direct measurements, in contrast to regional estimation (Le Quere et al., 2009). Komhyr et al. (1989) established that CO2 measurements made at the Mauna Loa Observatory are conducted with great precision. The observatory, which measures the mole fraction of CO2 using infrared absorption, gives quantitative information about its global trend. Kaufmann et al. (2008) used autoregression to test their approach and conclude that measurements based on distance are adequate. This information contradicts the results of Le Quere et al. (2009), who developed a model that effectively replicated the activity of land-ocean CO2 sinks and its impact on the increase of CO2 levels from 2000-2008.
The increase of CO2 sources has been attributed to fossil fuel, and carbon intensity depends on the regional scale (Weber et al., 2008). The interplay between carbon-driven emissions and the resultant climate change effect is already well-known (Plass, 1956; Kaplan, 1960; Tans et al., 1989; Siegenthaler and Sarmiento, 1993; Houghton and Hackler, 1999; Seekell, 2016). An estimation using different model methods has shown that, over the decades, a significant proportion of CO2 released into the atmosphere has been absorbed by the biosphere and oceans (Rayner et al., 1999; Battle et al., 2000; Le Quéré et al., 2003; Manning and Keeling, 2006). The variability in relative carbon sequestration by the two sinks causes variability in the overall CO2 atmospheric loading and, by extension, it produces climate variability and change. The idea of using greening in carbon capture and storage aims to limit the corresponding concentrations and stabilize the Earth’s radiative forcing. Prior studies have demonstrated the applications of greening (Bartlett et al., 1990; Hope et al., 1993; Sharkey, 2006; Rigby et al., 2008; Turnbull et al., 2011; McKain et al., 2012; Rogers et al., 2017), therein validating the option as applicable to broader efforts in the mitigation of climate change. Calculating dynamical trends in CO2, temperature and normalized index across the regions of Africa, Europe and Asia would be a worthwhile, relevant exercise. The results of the dynamical trend will provide a possible mitigation strategy. The idea of using greening to regulate surface CO2 storage aims to limit the corresponding increasing frequency of heat waves and stabilize the Earth’s radiative forcing. The normalized difference vegetation index (NDVI) uses the quotient of near infrared radiation and visible wavelength radiation. Its usage is in global climate change and environmental issues. NDVI compensates for reflectivity changes at certain frequencies and tracks greening. The degree of greenness ranges from -1 to +1, where the positive value is indicative of greening. Using the seasonal index to validate the relationship between CO2 and NDVI is not an approach which has been well investigated. There is a robust review of literature from multiple prior analyses, but no one has investigated the interplay between NDVI and CO2 with latitude, seasons and the solar cycle. This paper discusses the current state of understanding of this interplay and attempts to bridge the knowledge gap between seasonal trends of carbon and their possible mitigation strategy.
Dagg and Lafleur (2014) estimated the NDVI and CO2 relationship over tundra vegetation regions. Their study used a portable chamber to measure CO2 fluxes and a handheld infrared camera for the estimation of NDVI and CO2 variations and their contribution to the understanding of landscape characteristics. They found a positive correlation of CO2 exchange over NDVI, and concluded that NDVI is a good predictor of CO2 exchange. In an attempt to estimate global CO2 concentration using MODIS, Guo et al. (2012) used a model to investigate the carbon cycle between biosphere and atmosphere combined with an observation from eddy covariance. Although the study covered a short period (2009-2011), it aimed to map monthly CO2 concentrations over Africa, Eurasia, and North and South America using vegetation cover. The results indicate that seasonal trends of CO2 have local variations and their maxima and minima occur during summer and winter, respectively, over Eurasia, while the fluctuations over Africa were quite small. The study, which also features other parameters, found that land surface temperature has a better potential for validating CO2 concentration trends. La-Puma et al. (2007) found a good relationship between NDVI and CO2 during the summer, when biomass reaches its peak. Their results suggest that NDVI should not be included in model calibration over CO2 fluxes. These results conflict with Cihlar et al. (1992), who measured NDVI and CO2 interactions from an aircraft over agricultural lands in Kansas, finding that their relationship was linearly correlated over a single day, but nonlinearly over a longer period. Stow et al. (1998) found similar results and suggested that NDVI is a good parameter to predict variations in CO2. The same results were corroborated in other studies (McMichael, 1999; Markon and Peterson, 2002; Boelman et al., 2003). However, some authors (Huete et al., 2002; Xiao et al., 2004) gave more credit to an enhanced vegetation index (EVI). Lim et al. (2004) suggested that NVDI and CO2 data correlation should not be done directly, for it does not give a better result. Their results show that CO2 and NDVI have a strong seasonal covariance.
2. Methodology and data source
In this study, validation of the seasonal variations of CO2 under the influence of NDVI using a seasonal index method of analysis and trend with temperature were validated. The data used in this study were retrieved from the National Aeronautics and Space Administration (NASA GESDISC) via https://giovanni.gsfc.nasa.gov. CO2 concentration was retrieved in parts per million (ppm) from an atmospheric infrared sounder (AIRS) with monthly temporal resolution and 2º × 2.5º spatial resolution (ID AIRX3C2M v. 005), and temperature with monthly temporal resolution and 1º spatial resolution (ID AIRX3STM v. 006). NDVI data with monthly temporal resolution and 0.05º spatial resolution was retrieved from MODIS-Terra for the period 2003-2008 over the study area (Fig. 1), which comprises Africa (Nigeria, Mauritania, Congo, Sudan), Europe (France, Finland, Turkey, Ukraine), and Asia (China, Mongolia, India, Afghanistan). Since the AIRS collection of CO2 data was discontinued in 2012, the corresponding active solar years were used in the analysis. AIRS is an essential tool for the study of CO2 distribution. The seasonal index was obtained by dividing the yearly mean for each monthly period (Eq. ENT#091;1ENT#093;); however, a difference correlation of NVDI-CO2 in percentage was placed to investigate the underlying covariance of both parameters and validate the seasonal change, maintaining the solar activity cycle. The method described in Ferreira et al. (2003) was adopted to examine the spectral seasonal index of CO2 over the selected regions. The regions in study vary regarding their meteorological cycles; therefore, the idea of using seasonality is to characterize its temperature variation, which gives a better insight for making new advances.
The aim is to validate the CO2 variation as a function of vegetation index by grouping the years based on maximum (2003-2004), intermediate (2005-2006) and very low (2007-2008) solar activity. Table I presents the coordinates of the study locations. The study regions have different meteorological characteristics and different CO2 concentrations, resulting from local biomass burning. Therefore, for a broader view, a larger study region will provide adequate results in characterizing the seasonal trend of CO2. Correlation analysis will reveal if there is any distinctive characteristic in the spatial heterogeneities of either variable, as well as temporal evolution from seasons to years as seen in Eq (2). Chen et al. (2001) used a similar approach to validate seasonal changes in the increase/decrease of CO2 concentrations.
where n is the number of variables, ∑xy is the sum of the product of the paired NDVI and CO2 variables, ∑x is the sum of NDVI variables, ∑y is the sum of CO2 variables, ∑x 2 is the sum of squared x variables, ∑y 2 is the sum of y variables, and r stands for correlation.
Station | Latitude | Longitude |
Nigeria | 9.0820º N | 8.6753º E |
Mauritania | 21.0079º N | 10.9408º W |
Congo | 4.0383º S | 21.7587º E |
Sudan | 12.8628º N | 30.2176º E |
France | 46.2276º N | 2.2137º E |
Finland | 61.9241º N | 25.7482º E |
Turkey | 38.9637º N | 35.2433º E |
Ukraine | 48.3794º N | 31.1656º E |
China | 35.8617º N | 104.1954º E |
Mongolia | 46.8625º N | 103.8467º E |
India | 20.5937º N | 78.9629º E |
Afghanistan | 33.9391º N | 67.7100º E |
3. Results and discussion
An adequate validation of the seasonal CO2 trend requires a broader regional sampling technique. This section discusses the scientific approach based on the region and presents new and deeper understanding of the NDVI and CO2 trend analysis as suggested by Lim et al. (2004). However, to obtain a significant result, the direct use of NDVI and CO2 data is not appropriate. The significance of this result shows clearly that seasonal trend analysis is potentially useful to study NDVI and CO2 variations and report the interactions on a regional basis.
3.1 Africa
Validations of the NDVI and CO2 relationship in different seasons over Africa are presented. Figure 2a shows a CO2 peak during winter, indicating high carbon concentrations; conversely, the incident photosynthetic uptake as the result of carbon storage via greening leads to a decline. Figure 2b, c shows the mass uptake of CO2 during the summer, which does not compensate the emission increase of this compound, leading as a result to a net accumulation of carbon. The inverse correlation of CO2 with NDVI indicates that this a good mitigation strategy to limit the concentration of carbon. Curtis and Wang (1998) obtained similar results and suggested an appropriate mitigation of CO2. In our results, NDVI and CO2 are positively correlated during the ascending solar activity cycle (Table II), indicating that NDVI records moderate values towards low emission of carbon. The growth in carbon concentration implies larger anthropogenic emission activities in Nigeria, corresponding to low greening. The results of this study also fill the knowledge gap discussed in Le Quere et al. (2009) regarding a year to year anthropogenic-driven CO2 emission rate. There is a substantial degree of variation in the degree of CO2 effect on temperature; however, the photosynthetic capacity during the growing season closely tracked the pace of warming and the temperature record shows convincingly that Africa has experienced the largest warming rate. Towards the intermediate solar activity, Figure 3a, b shows a decreasing trend of NDVI and an increase of CO2 during summer. This is an indicator of weak radiant energy, as displayed in Table II, which indicates a negative weak correlation of NDVI with CO2 and a strong correlation of CO2 as function of NDVI. However, towards the vernal and autumnal equinoxes, CO2 decreases with maximum NDVI at the corresponding months, manifesting a positive correlation. A similar result is seen in Barchivich et al. (2013), and it could therefore be concluded that anthropogenic carbon emission is seasonally dependent. In the descending solar activity cycle as seen in Table II and Figure 4a, b, CO2 reaches peaks during winter and summer, corresponding to a declining trend of NDVI at particular seasons. In Figure 5a, Mauritania displays a moderate CO2 concentration trend, while NDVI peaks during the vernal and autumnal seasons, indicating a reduced value of CO2 as seen in Figure 5a, b. However, in Figure 6a, b, the CO2 concentration reaches its maximum during the winter and summer, as seen during the intermediate solar activity year. Figure 7a, b shows similar results to Figure 5a, b, but varies in regard to the carbon storage rate. This shows that NDVI and CO2 variations strongly depend on seasonality.
Station | NDVI and CO2 |
Nigeria | |
2003-2004 | 0.03 |
2005-2006 | 0.73 |
2007-2008 | -0.36 |
Mauritania | |
2003-2004 | -0.44 |
2005-2006 | -0.54 |
2007-2008 | -0.67 |
Congo Kinshasa | |
2003-2004 | 0.18 |
2005-2006 | 0.02 |
2007-2008 | 0.22 |
Sudan | |
2003-2004 | -0.08 |
2005-2006 | -0.13 |
2007-2008 | -0.28 |
Seasonal signatures in the atmospheric concentration of CO2 reflect the changing balance of photosynthetic rate; however, biomass burning accounts for the annual cycle of CO2 increase. This result agrees with those of Markon and Peterson (2002). The CO2 and NDVI relationship rate indicates an inverse correlation (Table II). This indicates that the CO2 concentration rate is a landscape dynamic scale and depends on the seasons, which is shown in Figures 6a, b, 7a, b, 8a, b, 9a, b and 10a, b. The seasonal responses of NDVI depend on the atmospheric conditions, as depicted in Figures 11a, 12 and 13b. In Figure 11a, Sudan shows low concentrations of CO2 due to a high vegetation index during the ascending solar activity cycle, as depicted in Figure 11b, but also shows a reciprocal relation in Figure 13a, b.
3.2 Europe
The increasing or decreasing trend of CO2 is the major determinant of climate change. However, to validate this relationship, Figure 14a, b indicates a moderate CO2 concentration. This clearly shows that during the ascending solar activity cycle, the emission-driven CO2 was regulated. A similar result is depicted in Figure 15a, b, which shows carbon capture during both winter and summer. In Figure 16a there is a CO2 peak and a drastic increase in the vegetation index. Table II displays a positive weak correlation between NDVI and CO2. The decline in CO2 resulting from NDVI increases is shown in Figures 16a, 17a, 18a to 19a. In Figures 20, 21 and 22a, during the ascending, intermediate and weak solar activity, Turkey experiences maximum emissions of CO2 during winter and summer, which results from heating of the region, as confirmed in Figure 20b. This leads to the conclusion that high concentrations of CO2 caused the region’s heat wave, thereby resulting in climate change. NDVI and CO2 interactions are regional, and the range of change in climate also differs. This result is well supported by La-Puma et al. (2007). In Table III, NDVI and CO2 show a weak correlation towards the first and second years, but the moderate correlation in Turkey over the period confirms a high emission rate. In Figures 23a, 24 and 25b there is a moderate CO2 trend due to high concentrations of NDVI.
Station | NDVI and CO2 |
France | |
2003-2004 | 0.34 |
2005-2006 | 0.11 |
2007-2008 | 0.29 |
Finland | |
2003-2004 | 0.32 |
2005-2006 | 0.27 |
2007-2008 | 0.06 |
Turkey | |
2003-2004 | 0.69 |
2005-2006 | 0.62 |
2007-2008 | 0.58 |
Ukraine | |
2003-2004 | 0.47 |
2005-2006 | 0.28 |
2007-2008 | 0.48 |
3.3 Asia
TableIV shows a negative, weak correlation between NDVI and CO2 over China and Mongolia during the intermediate and weak solar activity years. This result is in agreement with Le Quere et al. (2009), as shown in Figure 27a, b, which suggests that over 30% of the CO2 concentrations in China are generated by industrial activities. However, it can be seen that the northern hemisphere is not sensitive to the dynamics of CO2 variation. Therefore, surface-based observation is recommended in this region. Additionally, the results of Murayama et al. (2004) confirm that direct measurement in the northern hemispheric zone is valuable. The results of this study contribute to this knowledge gap. Figures 28, 29, 30, 31, 32, 33, 34 and 35 show a moderate trend of CO2 as a result of an increase in NDVI during winter and autumn. Towards the intermediate solar activity year over India and Afghanistan (Figs. 33, 34, 35, 36 and 37) and the weak solar activity year over this last country, there is a seasonal trend of CO2 increase indicating high atmospheric accumulation leading to climate variability.
Station | NDVI and CO2 |
China | |
2003-2004 | 0.12 |
2005-2006 | -0.07 |
2007-2008 | -0.22 |
Mongolia | |
2003-2004 | 0.60 |
2005-2006 | 0.40 |
2007-2008 | 0.51 |
India | |
2003-2004 | -0.59 |
2005-2006 | -0.65 |
2007-2008 | -0.69 |
Afghanistan | |
2003-2004 | 0.83 |
2005-2006 | 0.52 |
2007-2008 | 0.61 |
4. Conclusion and policy recommendation
This study presents a combined assessment of continental CO2 dominance and seasonal changes. The results show that the sequestration of CO2 using greening is essential to maintain moderate atmospheric conditions. However, the study has revealed that the continental biomass burning rate is seasonally dependent and varies regionally. The variability of seasonally integrated NDVI presents a consistent picture of the green index and carbon mitigation rate across a region. The seasonal mean is valuable to enable the regional investigation of carbon concentration using vegetation index, in order to document seasonal responses to high carbon atmospheric content (Heiman et al. 1998). However, our results were also compared to the corresponding temperature variations. The seasonal changes of NDVI correspond to a terrestrial sink of regional CO2, mostly occurring during equinoctial months (Piao et al., 2008; Vesala et al., 2010). Additionally, the dynamical responses of rapid carbon concentration resulting from biomass burning have broad spatial and temporal extent. The regression analysis shows that Africa and Asia are the largest CO2 contributors, with an inverse relationship which indicates coupling of CO2 and NDVI. An extended validation including more regions is needed in future works, to better characterize the drivers of variability in CO2 time series.
The seasonal investigation of CO2 renders the direct atmospheric composition, which is similar to the results of Conway et al. (1994).
The seasonal variation of CO2 from the results of the present study analysis depends on geographic latitude.
The seasonal index analysis is better to analyze the trends in atmospheric CO2 variation.