Next Article in Journal
Developing an Exploratory Framework Linking Australian Aboriginal Peoples’ Connection to Country and Concepts of Wellbeing
Previous Article in Journal
The Influence of Fish Length on Tissue Mercury Dynamics: Implications for Natural Resource Management and Human Health Risk
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Long-Term Relationship between Population Growth and Vegetation Cover: An Empirical Analysis Based on the Panel Data of 21 Cities in Guangdong Province, China

1
Sustainable Development Research Center, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou, Guangdong 510640, China
2
Key Laboratory of Marginal Sea Geology, Chinese Academy of Sciences, Guangzhou, Guangdong 510640, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
Guangdong Social Sciences Association, Guangzhou, Guangdong 510050, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2013, 10(2), 660-677; https://doi.org/10.3390/ijerph10020660
Submission received: 3 December 2012 / Revised: 25 January 2013 / Accepted: 28 January 2013 / Published: 7 February 2013

Abstract

:
It is generally believed that there is an inverse relationship between population growth and vegetation cover. However, reports about vegetation protection and reforestation around the World have been continuously increasing in recent decades, which seems to indicate that this relationship may not be true. In this paper, we have taken 21 cities in Guangdong Province, China as the study area to test the long-term relationship between population growth and vegetation cover, using an AVHRR NDVI data set and the panel cointegrated regression method. The results show that there is a long-term inverted N-shaped curve relationship between population growth and vegetation cover in the region where there are frequent human activities and the influence of climate change on vegetation cover changes is relatively small. The two turning points of the inverted N-shaped curve for the case of Guangdong Province correspond to 2,200 persons·km−2 and 3,820 persons·km−2, and they can provide a reference range for similar regions of the World. It also states that the population urbanization may have a negative impact on the vegetation cover at the early stage, but have a positive impact at the later stage. In addition, the Panel Error Correction Model (PECM) is used to investigate the causality direction between population growth and vegetation cover. The results show that not only will the consuming destruction effect and planting construction effect induced by the population growth have a great impact on vegetation cover changes, but vegetation cover changes in turn will also affect the population growth in the long term.

1. Introduction

Vegetation plays an important role in energy conversion and material circulation of the Earth, and vegetation cover is a visual sign of ecosystem health [1]. Although the dynamic evolution of vegetation cover is affected by both natural and anthropogenic factors [2,3,4], the impact of anthropogenic factors on the vegetation cover change is more significant than that of natural factors in the regions with frequent human activities [5,6,7,8].
The available studies about population and vegetation cover have mainly focused on the destructive effects of human activities on vegetation [9,10,11,12]. Green vegetation was excessively logged and used by human beings to support their productive and living activities [13]. Industrial growth in urban areas, including highway and railway construction, water use, mining and power generation, have directly consumed a huge amount of vegetation and caused chemical pollution [14,15]. The development of agriculture has destroyed a large part of the original vegetation, such as forests [16,17]. Tourism activities do not only lead to the death of some plants, but also subsequently exert a great negative effect on the growth of plants [18,19,20]. The destructive effects of human activities on vegetation have exerted an increasingly severe impact on the ecological environment, accordingly resulting in grassland degradation, deforestation, desertification, farmland reduction, soil erosion and other serious consequences [21].
During recent years, studies concerning how human activities were improving vegetation cover began to increase with more vegetation protection and reforestation [22,23,24,25]. Population growth does not necessarily mean vegetation destruction, because economic, social, political, technological and other developments can actually promote reforestation and improve vegetation cover [26,27,28]. A study has suggested that the forest covers of Europe, North America and Asia were increasing, and more and more countries and regions were shifting from deforestation to reforestation [29]. Some data from France indicated that although the population had been growing, there had been an upward trend in forest cover due to ongoing reforestation activities in France since 1830 [30]. From 1960 to 2006, the French forest area expanded more than a quarter, while the total population increased from 42 million to 61 million [29].
In summary, human activities can either increase or reduce the vegetation cover, which in the long term may therefore not just show an inverse relationship between population growth and vegetation cover as many scholars have stated [31,32,33,34], and a U-shaped or inverted N-shaped curve relationship might be more possible, yet, empirical research on this issue are quite few. Some studies have already investigated a positive relationship between population growth and vegetation cover in recent decades [35,36], but they all used statistical forest cover indexes which may cause inaccurate and non-comprehensive results. Firstly, the forest cover cannot reflect the cover of other vegetation types in an area, such as arable lands, grasslands and gardens. Secondly, the forest cover also cannot reflect the vegetation cover of urban areas with frequent human activities. Thirdly, the diverse statistical standards of forest cover index and other human factors will result in errors. NDVI is an indicator of vegetation’s growth state and spatial distribution, and has a positive correlation with vegetation cover. It reflects all the vegetation cover on land, including forests, arable lands, grasslands, gardens and parks. A greater NDVI value indicates higher vegetation cover and better vegetation growth [37,38]. Therefore, this paper uses NDVI instead of forest cover indexes to reflect vegetation cover and takes 21 cities in Guangdong Province, China as the study area to test the long-term relationship between population growth and vegetation cover, using the panel cointegrated regression method.

2. Materials and Methods

2.1. Study Area

The studied area covers the entire Guangdong Province in China (20°13′N–25°31′N, 109°39′E–117°19′E) (Figure 1). The province, with a land area of 179,757 km2, is located in southeast China and ranked fifteenth in size among all the Chinese provinces. It is divided into 21 administrative units called cities, and each city is used as the basic unit to analyze the vegetation cover and population characteristics.
Figure 1. Map showing the geographic location of the study area.
Figure 1. Map showing the geographic location of the study area.
Ijerph 10 00660 g001
The surface of the province is dominated by forests, accounting for 58.9% in 2006. Meanwhile, arable lands, grasslands, urban and built-up lands, and bodies of water account for 22.4%, 0.2%, 9.8% and 4.7%, respectively [39].
The climate of the province is of a humid subtropical monsoon type [40,41]. Its annual precipitation is generally over 1,300 mm and annual average temperature is about 22 °C. Due to the adequate rainfall, mild weather and forest-dominated landscape, the influence of climate change on the vegetation cover changes in this region is small [42].
As the forefront area of China’s reform and opening up, Guangdong Province has maintained rapid population growth since the 1980s, with annually average population growth rate above 2.27%. Due to the diversity of economic development and geographical environment, there is a significant difference in population growth among 21 cities of the province. The population grew rapidly in the developed Pearl River Delta Region, and, for example, Shenzhen’s permanent residents grew from 0.45 million to 8.46 million in 1982–2006, increasing by 1.783% in just 25 years; the population grew steadily in western Guangdong, and Zhanjiang’s permanent residents have grown from 4.49 million to 6.90 million in 1982–2006, increasing by 53.7% in 25 years; finally, the population grew slowly in northern Guangdong, and Shaoguan’s permanent residents have grown from 2.43 million to 2.95 million in 1982–2006, increasing by 21.2% in 25 years. This population growth diversity among the 21 cities made Guangdong Province form different population distribution patterns at different periods. This phenomenon provides us with an ideal panel data set to carry out empirical research on the long-term relationship between population growth and vegetation cover.

2.2. Data Sources and Processing

NDVI is a ratio of the near-infrared (NIR) and visible (VIS) radiances. It is calculated as: NDVI = (NIR − VIS)/(NIR + VIS). The NDVI data set used in this study was provided by the Environmental and Ecological Science Data Center for West China, National Natural Science Foundation of China (http://westdc.westgis.ac.cn). It is generated from the Global Inventory Monitoring and Modeling Studies (GIMMS) group, derived from the NOAA/AVHRR land data set, at a spatial resolution of 8 km × 8 km and taken at half a month intervals, for the period from January 1982 to December 2006. This data set has been processed with radiometric correction, geometric correction, replacing bad lines and cloud clearing. In order to explore inter-annual variations in vegetation cover, we use the Maximum Value Compositing (MVC) technique to reconstruct the half-monthly NDVI data set into the annual maximum NDVI data set. Data of population density is calculated by city according to the Guangdong Statistical Yearbooks from 1982 to 2006 (Table A1). The annual maximum NDVI values were averaged on the surface of the 21 cities to get yearly time series by cities (Table A2).
The trends in population density and NDVI among the 21 cities in the study area are diverse during the period from 1982 to 2006 (Figure 2, Figure 3). Some cities with fast population growth have experienced different NDVI trends from the cities with slow population growth. This phenomenon provides us an ideal panel data set to carry out the empirical research on the long-term relationship between population growth and vegetation cover. The statistical characteristics of the panel data are shown in Table 1.
Figure 2. Trends in population density of 21 cities during the period from 1982 to 2006. Note: A 5-year smoothing interval is used to clearly show the trends.
Figure 2. Trends in population density of 21 cities during the period from 1982 to 2006. Note: A 5-year smoothing interval is used to clearly show the trends.
Ijerph 10 00660 g002
Figure 3. Trends in NDVI of 21 cities during the period from 1982 to 2006. Note: A 5-year smoothing interval is used to clearly show the trends.
Figure 3. Trends in NDVI of 21 cities during the period from 1982 to 2006. Note: A 5-year smoothing interval is used to clearly show the trends.
Ijerph 10 00660 g003
Table 1. Descriptive statistics of population density and annual maximum NDVI data of 21cities in Guangdong Province during the period from 1982 to 2006.
Table 1. Descriptive statistics of population density and annual maximum NDVI data of 21cities in Guangdong Province during the period from 1982 to 2006.
P (Population Density)V (NDVI)
Mean0.6610060.588819
Median0.4250.62016
Maximum4.2870.75056
Minimum0.1340.350488
Std. Dev.0.6737030.097254
Observations525525
Cross sections2121

2.3. Methods

Based on a conceptual model, we constructed a simplified cubic polynomial regression equation to test the long-term relationship between population growth and vegetation cover, using the panel cointegrated regression method. In order to ensure the effectiveness of the fitting model, the two most commonly used panel unit root test methods, including the LLC test [43] and Fisher-ADF test [44], are applied to test the stationarity of the panel data (the LLC test is applicable for homogeneous panels, while the Fisher-ADF test is applicable for heterogeneous panels). Three kinds of panel cointegrated test methods, including Pedroni [45], Kao [46] and Johansen-Fisher [47], are applied to determine whether there is a long-term cointegrated relationship between population growth and vegetation cover or not. The random effects regression model is used to estimate the regression equation. The Panel Error Correction Model (PECM) is used to investigate the causality direction between population growth and vegetation cover.

3. Conceptual Model

The influences of population growth (the population growth here includes the natural population growth and immigration growth) on vegetation cover can be considered as two effects. One is the consuming destruction effect. Population growth can inevitably result in increasing demands for life necessities. To meet these demands, large areas with good vegetation cover will be exploited for construction of houses, roads, factories and shops, and many vegetation resources will be plundered, resulting in a vegetation cover decrease [5]. Another is the planting construction effect. Vegetation is an essential element for human development, and it can help improve the living environment as well as providing productive materials and a source of energy for humans. With the population growth, the demands for the ecological functions provided by vegetation will increase. To fulfill these demands, some constructive actions, such as vegetation protection, reforestation and so on, will be carried out, resulting in a vegetation cover increase [48].
Of course, vegetation cover can be influenced by other factors, especially the long-term inter-annual climate changes, which affect vegetation cover mainly through fluctuation of rainfall. However, in the studied area, the annual rainfall was more than 1,300 mm during the period from January 1982 to December 2006. It is considered where there are frequent human activities and the influence of climate change is relatively small on vegetation cover changes (anthropogenic factors play a dominant role in the vegetation cover changes in this region, while the natural factors have a small impact). Based on the above-stated two effects, we try to construct a conceptual model as follows: the long-term relationship between population growth and vegetation cover can be separated into three stages in a region where there are frequent human activities and the influence of climate change on vegetation cover changes is small. At the first stage, there is an inverse relationship between population growth and vegetation cover. As the vegetation cover is relatively high and the public facilities are imperfect at the early stage of population growth, the vegetation cover decreases fast with the population growth when the consuming destruction effect is much stronger than the planting construction effect. Subsequently, due to the improvement of public facilities, the consuming destruction effect induced by population growth tends to become small. Besides, the ecological deterioration resulting from vegetation cover decrease makes people aware of the importance of the ecological functions of vegetation. To improve the vegetation cover, the government implements some measures, such as more intensive residence areas, vegetation protection, reforestation and so on. Consequently, the planting construction effect gradually offsets the consuming destruction effect, and the inverse relationship between population growth and vegetation cover weakens at the latter part of the first stage.
At the second stage, there is a positive relationship between population growth and vegetation cover. With the development of society and the economy, the ecological functions of vegetation arouse more people’s attention. Then, some new green technologies are developed, which makes many places without vegetation cover (such as roofs, walls, fences, bridges, etc.) be planted and the permanent vegetation cover protected. At the same time, as the population continues to grow and the government strengthens the vegetation protection, to meet the need for vegetation resources, the food and resource sectors tend to increase imports, resulting in a further decline of the consuming destruction effect. These activities cause an increase of vegetation cover. At this stage, the planting construction effect directly or indirectly induced by population growth has already exceeded the consuming destruction effect, while the relationship between population growth and vegetation cover becomes positive, and the vegetation cover tends to increase slightly with the population growth.
At the third stage, there is an inverse relationship between population growth and vegetation cover again. With the population continuing to expand beyond a certain limit, the consuming destruction effect will again surpass the planting construction effect, and then vegetation cover will tend to decrease with population growth.
In the long term, we will observe an inverted N-shaped curve relationship between population growth and vegetation cover (Figure 4). What needs to be explained is that, with the economic development and birth rate changes, the third stage may not appear in some regions because of the zero or negative population growth. In these regions, the vegetation cover may be maintained at the level of the second stage or even improved with the population decrease.
To test the long-term inverted N-shaped curve relationship between population growth and vegetation cover, we constructed a simplified cubic polynomial regression equation after the Environmental Kuznets Curve (EKC) [49] as follows:
V = c + β1P + β2 P2 + β3 P3 + u
where V (NDVI) is the index reflecting the vegetation cover, P (population density) is the index reflecting the population, c is a constant, u is the random error, the parameters β1, β2 and β3 are the coefficients of the first, second and third term of P, respectively. If β3 > 0, β2 < 0, β1 > 0, there is a N-shaped curve relationship between population growth and vegetation cover, which means that the vegetation cover tends to increase first, then decrease and rise again with the population growth. If β3 < 0, β2 > 0, β1 < 0, there is an inverted N-shaped curve relationship between population growth and vegetation cover, which means that the vegetation cover tends to decrease first, then increase and decline again with the population growth. If β3 = 0, β2 > 0, β1 < 0, there is an U-shaped curve relationship between population growth and vegetation cover, which means that vegetation cover tends to decrease first and then increase with the population growth. If β3 = 0, β2 < 0, β1 > 0, there is an inverted U-shaped curve relationship between population growth and vegetation cover, which means that the vegetation cover tends to increase first and then decrease with the population growth. If β3 = β2 = 0, β1≠ 0, there is a linear relationship between population growth and vegetation cover.
Figure 4. The long-term inverted N-shaped curve relationship between population growth and vegetation cover.
Figure 4. The long-term inverted N-shaped curve relationship between population growth and vegetation cover.
Ijerph 10 00660 g004

4. Results and Discussion

4.1. Cointegration Relationship

In order to ensure the effectiveness of the fitting model, we test the stationarity of the panel data. The results (Table 2) show that all variables in the panel series are unable to reject the null hypotheses of non-stationarity at the specified significance levels, while all the 1st differential panel series reject the null hypotheses of non-stationarity at a significance level of 10% at least. This indicates that all variables can be considered as non-stationary panel series with first order integration (denoted as I(1)).
Table 2. Test statistics of the panel unit root test for the variables involved in Equation (1).
Table 2. Test statistics of the panel unit root test for the variables involved in Equation (1).
Test methodLevin, Lin & Chu tADF - Fisher Chi-square
Panel seriesPrimary1st differencePrimary1st difference
V1.04913−30.3697 ***36.4358595.897 ***
P−2.89544 ***−3.74377 ***27.521161.1483 **
P2−0.52219−3.02647 ***23.64756.5126 *
P31.56007−2.35949 ***23.049957.9320 *
The primary panel series are tested with intercept, while the 1st differential panel series are tested without intercept; *, ** and *** indicate that the null hypotheses of a unit-root are rejected at significance level of 10%, 5% and 1% respectively.
Since all variables are non-stationary panel series with first order integration (I(1)), we proceed with the panel cointegration test to determine whether there is a long-term relationship between population density and NDVI or not. The results (Table 3) show that all tests reject the null hypotheses of no cointegration among the four variables at significance level of 5%, which suggests that there is a long-term cointegration relationship between NDVI and population density.
Table 3. Panel cointegration tests.
Table 3. Panel cointegration tests.
PanelGroup
v-Strho-StPP-StADF-Strho-StPP-StADF-St
Pedroni−2.677−3.7234 ***−10.0978 ***−3.9372 ***−2.9945 ***−17.3462 ***−7.0680 ***
Kaot-Statistic: 1.7008 **
FisherFisher Stat. *: 474.2 ***
** and *** indicate that the null hypotheses of no cointegration are rejected at 5% and 1% levels respectively.

4.2. A Long-term Inverted N-shaped Curve

We use the Hausman test to determine whether a fixed or random effects model is appropriate to estimate the Equation (1), and the result shows that the null hypothesis of a random effects model isn’t rejected at 10% significance level. So we use the random effects model to estimate the equation, and the regression results (Table 4) show that there is a long-term inverted N-shaped curve relationship between population growth and vegetation cover, because the null hypotheses that the coefficients of variables are equal to zero are all rejected at a significance level of 1%, and the coefficients meet the condition β3 < 0, β2 > 0 and β1 < 0. The panel regression curve of population density and NDVI is shown in Figure 5. The first turning point, where the relationship between population density and NDVI changes from negative to positive correlation, occurs at the population density of about 2,200 persons·km−2. The second turning point, where the relationship between population density and NDVI changes from positive to negative correlation, occurs at the population density of about 3,820 persons·km−2. Most of the 21 cities were still at the first stage when there was an inverse relationship between population growth and vegetation cover. However, as the government continued to strengthen environmental and vegetation protection, some cities with a higher population density, such as Guangzhou, Zhuhai, Foshan, Zhongshan and so on, had entered the latter part of the first stage when the devegetation process had slowed down. Dongguan and Shantou had entered the second stage when the relationship between population growth and vegetation cover became positive. Shenzhen, whose population density was 4,287 persons·km−2 and vegetation cover was higher than that of Dongguan and Shantou, had just entered the third stage when the inverse relationship between population growth and vegetation cover appeared again.
Table 4. Estimation results from the panel random effects regression.
Table 4. Estimation results from the panel random effects regression.
VariablePP2P3c
Coefficient−0.228456 ***0.082598 ***−0.009194 ***0.684206 ***
(−8.819210)−5.592608(−3.940781)−41.59666
R-squared0.334048
F-statistic87.11297 ***
*** indicates that the null hypothesis that the coefficient is zero is rejected at 1% level.
Figure 5. The panel regression curve of population density and NDVI.
Figure 5. The panel regression curve of population density and NDVI.
Ijerph 10 00660 g005
In contrast to other studies’ inverse relationship [31,32,33,34], our empirical result shows that there is a long-term inverted N-shaped curve relationship between population growth and vegetation cover. It should be noted that 2,200 persons·km−2 and 3,820 persons·km−2 are the two turning points of the inverted N-shaped curve based on the experiences of Guangdong Province’s 21 cities, where there are similar natural conditions and social systems. However, they can just provide a reference range for similar regions of the World, where there are also frequent human activities and the influence of climate change on vegetation cover changes is also small, because different social systems and natural conditions will affect the location of the two turning points in different regions. For example, there will be faster vegetation restoration and the turning point will come earlier in more democratic regions, where people’s ecological needs can be better reflected.

4.3. Relationships in Different Regions

Guangdong is a province with unbalanced regional development, consisting of the developed Pearl River Delta region which has experienced a fast urbanization and the undeveloped peripheral region which has experienced a relatively slow urbanization. As there is big difference among 21 cities in different regions, we classify the cities into two types. One type is dominated by high-urbanization cities located in the Pearl River Delta region, including Guangzhou, Shenzhen, Zhuhai, Foshan, Huizhou, Dongguan, Zhongshan, Jiangmen and Zhaoqing. The other type is dominated by low-urbanization cities located in the peripheral region, including Shantou, Shaoguan, Heyuan, Meizhou, Shanwei, Yangjiang, Zhanjiang, Maoming, Qingyuan, Chaozhou, Jieyang and Yunfu. Then we draw the scatterplots of population density and NDVI with fitted curve, and find that there is an inverted N-shaped curve relationship between population growth and vegetation cover in the high-urbanization cities from 1982 to 2006 (Figure 6), while there is a negative linear relationship in the low-urbanization cities (Figure 7). Due to the free migration of population among cities in Guangdong province, the cities with higher urbanization would have larger population. Therefore, according to the conceptual model, we can explain that the cities with higher urbanization would enter the second and the third stage earlier. It also states that the population urbanization may have a negative impact on the vegetation cover at the early stage [50,51,52], but have a positive impact at the later stage.
Figure 6. The scatterplot of population density and NDVI with fitted curve in high-urbanization cities from 1982 to 2006.
Figure 6. The scatterplot of population density and NDVI with fitted curve in high-urbanization cities from 1982 to 2006.
Ijerph 10 00660 g006
Figure 7. The scatterplot of population density and NDVI with fitted curve in low-urbanization cities from 1982 to 2006.
Figure 7. The scatterplot of population density and NDVI with fitted curve in low-urbanization cities from 1982 to 2006.
Ijerph 10 00660 g007

4.4. Granger Causality

According to the Granger Representation Theorem [53], if the variables are cointegrated, there must be at least one direction of causal relationship between them. The results (Table 5) show that the null hypotheses that the coefficients of panel error correction terms ecm(v) and ecm(p) are equal to zero are rejected at significance levels of 1% and 5%, respectively, which indicates that V is Granger causal for P and P is Granger causal for V in the long term.
Table 5. The estimated results of the Panel Error Correction Model (PECM).
Table 5. The estimated results of the Panel Error Correction Model (PECM).
D(V)D(P)
VariableCoefficientVariableCoefficient
c−0.005207 *** (−3.950969)c0.005245 *** (5.450232)
ecm(v)−0.118505 *** (−2.715071)ecm(p)−0.005538 ** (−2.152945)
D[V(−1)]−0.392603 *** (−6.946931)D[V(−1)]−0.038606 (−1.123531)
D[V(−2)]−0.182085 *** (−3.956092)D[V(−2)]−0.000284 (−0.009368)
D[P(−1)]−0.006660 (−0.111667)D[P(−1)]1.165810 *** (26.33352)
D[P(−2)]−0.028675 (0.486901)D[P(−2)]−0.383236 *** (−8.155351)
ecm(v) and ecm(p) are the panel error correction terms; ecm(v) is the lag-1 panel series of residuals in fixed effects model regression of V on P, and ecm(v) is the lag-1 panel series of residuals in random effects model regression of V on P; We select the lag length in 2; D[V(−1)] and D[V(−2)] are the lag-1 and lag-2 1st differential panel series of V respectively; D[P(-1)]and D[P(−2)]are the lag-1 and lag-2 1st differential panel series of P respectively; ** and *** indicate that the null hypotheses that the coefficients of variables are equal to zero are rejected at 5% and 10% significance levels respectively.
We conclude that there is a long-term interactive relationship between population growth and vegetation cover. In other words, the impact of population growth on the vegetation cover change is a long-term process, and the long-term changes of vegetation cover will restrict the population growth. The null hypotheses that the regression coefficients of D(V) on D[P(−1)] and D[P(−2)], and D(P) on D[V(−1)] and D[V(−2)] are equal to zero are not rejected at significance level of 10% respectively, which indicates that V is not Granger causal for P and P is not Granger causal for V in the short term. We conclude that there is not a short-term interactive relationship between population growth and vegetation cover. In other words, the short-term changes of population may affect the vegetation cover changes little, while the short-term changes of vegetation cover may also affect the population little.

5. Conclusions

The purpose of this paper was to prove that there is a long-term inverted N-shaped curve relationship between population growth and vegetation cover in a region where there are frequent human activities and the influence of climate change is relatively small on vegetation cover changes [42]. In summary, we have obtained the following conclusions:
(1) There can be a positive relationship between population growth and vegetation cover. According to the empirical results of 21 cities in Guangdong Province, China, the vegetation cover will increase slightly with the population growth when the population density is in the range of 2,200–3,820 persons·km−2. Although the two turning points of population density from Guangdong Province’s experiences may not be entirely applicable to other parts of the World, it can provide a reference range for similar regions. At present, the global population is still growing, but some regions’ vegetation cover has been improved. In these regions, a lot of new technology (such as roof greening) is being developed to increase the vegetation cover, and many new measures (such as increasing imports) are implemented to reduce the consumption of vegetation resources. With the economic and social development, people are asking for greener environments and healthier living. As long as the population does not expand beyond a certain limit, the vegetation cover can remain unchanged or even improved in many regions.
(2) Unlimited population growth will eventually lead to a significant reduction of vegetation cover. According to the empirical results of 21 cities in Guangdong Province, China, the vegetation cover will decrease with the population growth when the population density is below 2,200 persons·km−2 or above 3,820 persons·km−2. For some regions in the World, the population growths are still at the first stage and are leading to a significant reduction of vegetation cover on land. For other regions, although the vegetation covers are experiencing a recovery with more people paying close attention to the ecological functions of vegetation, the degrees of vegetation restoration are low. If the population growth continues to expand beyond a certain limit, the consuming destruction effect will surpass the planting construction effect again, and the inverse relationship between population growth and vegetation cover will appear eventually.
(3) Population urbanization may free up space for vegetation ecological construction at some stage. According to the empirical results, there is an inverted N-shaped curve relationship between population growth and vegetation cover in the high-urbanization cities while there is only a negative linear relationship in the low-urbanization cities. This result shows that the cities with higher urbanization may have chances to improve their vegetation cover at the later stage, when the land saving from the urban agglomeration effect is larger than the land occupation from the urban expansionary effect, and more people promote more green building.
(4) Population growth and vegetation cover affect each other in the long term. According to the results of the Granger causality test, there is an interactive relationship between population growth and vegetation cover in the long term, but not in the short term. This result shows that not only will the consuming destruction effect and planting construction effect induced by the population growth have a great impact on vegetation cover changes in the long term, but vegetation cover changes in turn will also affect the population growth. For example, the suburban areas with better vegetation cover and ecological environment are attracting more urban population to move into them.
Finally, this paper has just tested a long-term inverted N-shaped curve relationships between population growth and vegetation cover through the panel data of 21 cities in Guangdong Province, China, where there are frequent human activities and the climate is of a humid subtropical monsoon type. However, maybe there will be the same relationships in most of the regions where there are frequent human activities, including humid, sub-humid, semiarid and arid climate regions, because anthropogenic factors may also play a dominant role in the vegetation cover changes in these regions. Therefore, further study is necessary to provide validation using data from more regions.

Appendix

Table A1. The average annual permanent population density of 21 cities in Guangdong Province from 1982 to 2006 (person·km−2).
Table A1. The average annual permanent population density of 21 cities in Guangdong Province from 1982 to 2006 (person·km−2).
YearGZSZZHSTFSSGHYMZHZSWDGZSJMYJZJMMZQQYCZJYYF
19826922092221,492630134155239177339468569339253363383197156614714247
19837032672261,522638135158244180347474573341256369388200159623727250
19847143422311,549647137160247183356480579344259374393202161632738253
19857274152391,574656140162250186364486585347262380397205163640748256
19867404652681,596667142164253188372495593349265387403207165646755258
19877545103021,616679144166256191380503601353268394410210167653762262
19887685783221,636692146165261194389510611357273402419212169660770265
19897826703441,657704148164265197398518621361277409428215171668779268
19908217923681,732749152162262203410617660364279426433216172681791271
19918731,0103931,819809156160257211420762717368279441434216172694808272
19929131,2674191,861857155158255218424868764373279445436217172702827273
19939551,5464481,904907154156253225429988814377279449439218171710845273
19949991,9174781,9499601531542512334341,125868382278454441219170718864274
19951,0452,2075101,9941,0161531532492404381,282925387278458443221169727884274
19961,0932,3865452,0401,0751521512472494431,460986392278462446222168735904275
19971,1432,5885822,0881,1481511492452574481,6621,051396278467448223167744924275
19981,1962,8376222,1361,2141511472442654531,8931,120401278471451224166753945276
19991,2513,1056642,1861,2741501462422744582,1561,194406278476453225166761967276
20001,3083,4157092,2391,3491491442402844632,4561,273412278481456227165771989277
20011,3393,6517472,2911,4071501472422944752,6351,3194182814914652301687831,013281
20021,3333,7677702,3251,4371521542473044892,6561,3294232865044772351737941,030287
20031,3173,9047892,3451,4561541602513125012,6571,3394252895134862391777991,041291
20041,3044,0438112,3711,4801561662553215142,6591,3464282935244972431828071,054295
20051,2884,1708312,3911,5011581732583295262,6611,3504302965335062461868121,065299
20061,2954,2878492,4021,5151591762593355322,7001,3694302985375142481888151,072301
Data of average annual permanent population density are calculated according to Guangdong Statistical Yearbook from1981 to 2006; GZ, SZ, ZH, ST, FO, SG, HY, MZ, HZ, SW, DG, ZS, JM, YJ, ZJ, MM, ZQ, QY, CZ, JY and YF are short for Guangzhou, Shenzhen, Zhuhai, Shantou, Foshan, Shaoguan, Heyuan, Meizhou, Huizhou, Shanwei, Dongguan, Zhongshan, Jiangmen, Yangjiang, Zhanjiang, Maoming, Zhaoqing, Qingyuan, Chaozhou, Jieyang and Yunfu respectively.
Table A2. The annual maximum NDVI of 21 cities in Guangdong Province from 1982 to 2006.
Table A2. The annual maximum NDVI of 21 cities in Guangdong Province from 1982 to 2006.
YearGZ SZ ZH STFSSG HY MZ HZ SW DG ZS JM YJ ZJ MM ZQ QY CZ JY YF
1982 0.5140.4240.3640.3640.4510.6610.6500.6580.6050.4480.4430.4190.5390.6090.4850.6380.6630.6410.5740.5460.671
1983 0.5600.5040.4190.3920.5130.7280.7130.7110.6620.5070.4760.5020.5790.6230.5160.6670.6890.6870.6140.5980.699
1984 0.5600.4700.4020.4130.4880.7450.6860.6790.6420.4730.4500.4800.5640.6070.5060.6420.6970.7210.5870.5690.673
1985 0.5570.4760.3860.3830.4890.7260.6760.6680.6310.4800.4600.4460.5270.5760.5020.6330.6670.6950.5600.5520.652
1986 0.5490.4830.3980.4000.5070.7020.6800.6790.6350.4940.4560.4670.5500.6120.5300.6480.6740.6800.5730.5540.672
1987 0.5480.4670.4220.3920.4580.7300.7100.7100.6530.4890.4390.4430.5600.5960.5220.6300.7020.7010.5970.5790.669
1988 0.5880.4820.4110.3950.4970.7310.7000.6990.6580.4970.4750.4710.5730.6500.5650.7080.7020.7150.5850.5760.716
1989 0.5360.4450.3790.3830.4500.6840.6540.6650.6210.4850.4440.4430.5360.5930.5330.6420.6750.6710.5940.5600.660
1990 0.5430.4570.3800.3860.4590.7210.6730.6810.6270.4830.4700.4460.5310.6070.5240.6280.6540.6920.6000.5750.663
1991 0.5650.4400.3850.3890.4580.7000.6800.6760.6300.4900.4420.4730.5430.5800.5030.6200.6590.6810.6070.5580.649
1992 0.5360.4300.3790.3690.4260.7160.6950.6830.6430.4770.4440.4180.5520.6030.5370.6510.6730.6800.5780.5590.668
1993 0.5220.4210.3690.3940.3840.7060.6770.6750.6410.4810.4000.4070.5420.6040.5540.6510.6900.6910.5710.5680.693
1994 0.5260.4220.3350.3630.3670.7220.6820.6560.6370.4720.3940.3720.5220.6040.5380.6660.6850.7190.5710.5550.674
1995 0.5340.4360.3680.3680.3930.6940.6750.6640.6290.4780.4050.3930.5120.5980.5310.6360.6930.6860.5570.5440.680
1996 0.4860.4030.3390.3530.3790.6600.6470.6560.5930.4610.3860.3690.5100.5920.5520.6480.6470.6510.5590.5320.650
1997 0.4760.4280.3390.3530.3680.6970.6660.6680.5930.4450.3890.3780.4990.5780.5200.6250.6210.6460.5510.5240.618
1998 0.5070.4550.3920.3770.4060.6690.6440.6550.6020.5010.4070.4040.5250.5770.5320.6330.6590.6450.5890.5650.655
1999 0.4730.4170.3460.3690.3730.6580.6470.6540.5790.4480.3830.4030.5190.5830.5430.6220.6290.6280.5590.5310.634
2000 0.4840.4320.3650.3620.3710.6620.6320.6240.5830.4580.3940.4060.5180.5920.5250.6220.6230.6470.5520.5140.631
2001 0.4720.4240.3440.3800.3770.7040.6650.6500.5890.4650.3840.3820.5030.5620.5090.5890.6390.6780.5600.5370.636
2002 0.4800.4000.3590.3510.3810.6850.6530.6430.5820.4590.3840.3870.5100.5670.5380.6020.6220.6500.5490.5190.617
2003 0.4780.4330.3800.3610.3830.6600.6410.6490.5860.4750.3750.3890.5350.5850.5260.6100.6440.6480.5660.5340.638
2004 0.4680.4020.3470.3510.3550.6340.6250.6300.5680.4590.3700.3750.4890.5550.5370.6080.5970.5980.5480.5140.610
2005 0.4380.3950.3500.3410.3470.5760.5890.6000.5440.4480.3470.3510.4980.5540.5200.6080.5860.5640.5370.5070.603
2006 0.4650.4360.3520.3440.3720.6020.6000.6050.5660.4680.3870.3920.5010.5480.5540.6110.5890.5740.5330.5070.610
The data set is provided by Environmental and Ecological Science Data Center for West China, National Natural Science Foundation of China (http://westdc.westgis.ac.cn); GZ, SZ, ZH, ST, FO, SG, HY, MZ, HZ, SW, DG, ZS, JM, YJ, ZJ, MM, ZQ, QY, CZ, JY and YF are short for Guangzhou, Shenzhen, Zhuhai, Shantou, Foshan, Shaoguan, Heyuan, Meizhou, Huizhou, Shanwei, Dongguan, Zhongshan, Jiangmen, Yangjiang, Zhanjiang, Maoming, Zhaoqing, Qingyuan, Chaozhou, Jieyang and Yunfu respectively.

Acknowledgments

The authors are indebted to Caiming Tang and Jingjing Wu for their assistance in improving the readability of this manuscript.

References

  1. Cihlar, J.; Stlaurent, L.; Dyer, J.A. Relation between the normalized difference vegetation index and ecological variables. Rem. Sens. Environ. 1991, 35, 279–298. [Google Scholar] [CrossRef]
  2. Bradshaw, R.; Hannon, G. Climatic-change, human influence and disturbance regime in the control of vegetation dynamics within Fiby Forest, Sweden. J. Ecol. 1992, 80, 625–632. [Google Scholar] [CrossRef]
  3. Vicente-Serrano, S.M.; Lasanta, T.; Romo, A. Analysis of spatial and temporal evolution of vegetation cover in the Spanish Central Pyrenees: Role of human management. Environ. Manag. 2004, 34, 802–818. [Google Scholar] [CrossRef]
  4. Kulakowski, D.; Bebi, P.; Rixen, C. The interacting effects of land use change, climate change and suppression of natural disturbances on landscape forest structure in the Swiss Alps. Oikos 2011, 120, 216–225. [Google Scholar] [CrossRef]
  5. Meyer, W.B.; Turner, B.L. Human-population growth and global land-use cover change. Ann. Rev. Ecol. Syst. 1992, 23, 39–61. [Google Scholar]
  6. Vanderknaap, W.O.; Vanleeuwen, J.F.N. Holocene vegetation succession and degradation as responses to climatic-change and human activity in the Serra de Estrela, Portugal. Rev. Palaeobot. Palyno. 1995, 89, 153–211. [Google Scholar] [CrossRef]
  7. Wang, Y.B.; Feng, Q.; Si, J.H.; Su, Y.H.; Chang, Z.Q.; Xi, H.Y. The changes of vegetation cover in Ejina Oasis based on water resources redistribution in Heihe River. Environ. Earth Sci. 2011, 64, 1965–1973. [Google Scholar] [CrossRef]
  8. Peng, J.; Xu, Y.Q.; Cai, Y.L.; Xiao, H.L. Climatic and anthropogenic drivers of land use/cover change in fragile karst areas of southwest China since the early 1970s: A case study on the Maotiaohe watershed. Environ. Earth Sci. 2011, 64, 2107–2118. [Google Scholar] [CrossRef]
  9. Foster, D.R.; Zebryk, T.; Schoonmaker, P.; Lezberg, A. Post-settlement history of human land-use and vegetation dynamics of a Tsuga canadensis (hemlock) woodlot in central new England. J. Ecol. 1992, 80, 773–786. [Google Scholar] [CrossRef]
  10. Andersen, U.V. Resistance of danish coastal vegetation types to human trampling. Biol. Conserv. 1995, 71, 223–230. [Google Scholar] [CrossRef]
  11. Kissling, M.; Hegetschweiler, K.T.; Rusterholz, H.P.; Baur, B. Short-term and long-term effects of human trampling on above-ground vegetation, soil density, soil organic matter and soil microbial processes in suburban beech forests. Appl. Soil. Ecol. 2009, 42, 303–314. [Google Scholar] [CrossRef]
  12. Shipigina, E.; Rees, W.G. Analysis of human impact on boreal vegetation around Monchegorsk, Kola Peninsula, using automated remote sensing technique. Polar Rec. 2012, 48, 94–106. [Google Scholar] [CrossRef]
  13. Watson, A. Soil-erosion and vegetation damage near ski lifts at cairn-gorm, Scotland. Biol. Conserv. 1985, 33, 363–381. [Google Scholar] [CrossRef]
  14. Ratcliffe, D.A. Post-medieval and recent changes in British vegetation—The culmination of human influence. New Phytol. 1984, 98, 73–100. [Google Scholar]
  15. Sasaki, N.; Takahara, H. Late holocene human impact on the vegetation around Mizorogaike Pond in northern Kyoto basin, Japan: A comparison of pollen and charcoal records with archaeological and historical data. J. Archaeol. Sci. 2011, 38, 1199–1208. [Google Scholar] [CrossRef]
  16. Delcourt, H.R. The impact of prehistoric agriculture and land occupation on natural vegetation. Trends Ecol. Evol. 1987, 2, 39–44. [Google Scholar] [CrossRef]
  17. Mann, D.; James, E.; Chase, J.; Beck, W.; Reanier, R.; Mass, M.; Finney, B.; Loret, J. Drought, vegetation change, and human history on Apa Nui (Isla de Pascua, easter Island). Quat. Res. 2008, 69, 16–28. [Google Scholar] [CrossRef]
  18. Burden, R.F.; Randerson, P.F. Quantitative studies of the effects of human trampling on vegetation as an aid to the management of semi-natural areas. J. Appl. Ecol. 1972, 9, 439–457. [Google Scholar] [CrossRef]
  19. Kerbiriou, C.; Leviol, I.; Jiguet, F.; Julliard, R. The impact of human frequentation on coastal vegetation in a biosphere reserve. J. Environ. Manag. 2008, 88, 715–728. [Google Scholar] [CrossRef]
  20. Lucas-Borja, M.E.; Bastida, F.; Moreno, J.L.; Nicolas, C.; Andres, M.; Lopez, F.R.; del Cerro, A. The effects of human trampling on the microbiological properties of soil and vegetation in Mediterranean mountain areas. Land Degrad. Dev. 2011, 22, 383–394. [Google Scholar]
  21. Bai, Z.; Dent, D. Recent land degradation and improvement in China. AMBIO 2009, 38, 150–156. [Google Scholar] [CrossRef]
  22. Rudel, T.K. Is there a forest transition? Deforestation, reforestation, and development. Rural Soc. 1998, 63, 533–552. [Google Scholar] [CrossRef]
  23. Nagendra, H. Drivers of Reforestation in Human-Dominated Forests. Proc. Natl. Acad. Sci. USA 2007, 104, 15218–15223. [Google Scholar] [CrossRef]
  24. Cortina, J.; Amat, B.; Castillo, V.; Fuentes, D.; Maestre, F.T.; Padilla, F.M.; Rojo, L. The restoration of vegetation cover in the semi-arid Iberian southeast. J. Arid. Environ. 2011, 75, 1377–1384. [Google Scholar] [CrossRef]
  25. Plieninger, T.; Schleyer, C.; Mantel, M.; Hostert, P. Is there a forest transition outside forests? Trajectories of farm trees and effects on ecosystem services in an agricultural landscape in eastern Germany. Land Use Pol. 2012, 29, 233–243. [Google Scholar] [CrossRef]
  26. Mather, A.S.; Needle, C.L. The forest transition: A theoretical basis. Area 1998, 30, 117–124. [Google Scholar]
  27. Mather, A.S.; Fairbairn, J. From floods to reforestation: The forest transition in Switzerland. Environ. Hist. 2000, 6, 399–421. [Google Scholar] [CrossRef]
  28. Nagendra, H. Reforestation and regrowth in the human dominated landscapes of South Asia. In Reforesting Landscapes: Linking Pattern and Process; Nagendra, H., Southworth, J., Eds.; Springer: Dordrecht, The Netherlands, 2010; Volume 10, p. 149. [Google Scholar]
  29. Kauppi, P.E.; Ausubel, J.H.; Fang, J.; Mather, A.S.; Sedjo, R.A.; Waggoner, P.E. Returning Forests Analyzed with the Forest Identity. Proc. Natl. Acad. Sci. USA 2006, 103, 17574–17579. [Google Scholar]
  30. Mather, A.S.; Fairbairn, J.; Needle, C.L. The course and drivers of the forest transition: The case of France. J. Rural Stud. 1999, 15, 65–90. [Google Scholar] [CrossRef]
  31. Allen, J.C.; Barnes, D.F. The causes of deforestation in developing countries. Ann. Ass. Am. Geogr. 1985, 75, 163–184. [Google Scholar] [CrossRef]
  32. Harrison, P. The Third Revolution: Environment, Population and a Sustainable World; Penguin Books: London, UK, 1993. [Google Scholar]
  33. DeFries, R.S.; Rudel, T.; Uriarte, M.; Hansen, M. Deforestation driven by urban population growth and agricultural trade in the twenty-first century. Nat. Geosci. 2010, 3, 178–181. [Google Scholar] [CrossRef]
  34. Hazra, S.; Sen, R.K. Population Growth and Forest: Degradation in North East India. In Population and Development in North East India; Deb, B.J., Ed.; Concept Publishing Company: New Delhi, India, 2010; pp. 99–110. [Google Scholar]
  35. Mather, A.; Needle, C. The relationships of population and forest trends. Geogra. J. 2000, 166, 2–13. [Google Scholar] [CrossRef]
  36. Nagendra, H. Reforestation and regrowth in the human dominated landscapes of South Asia. Reforest. Landsc. 2010, 149–174. [Google Scholar]
  37. Nordberg, M.L.; Evertson, J. Monitoring change in mountainous dry-heath vegetation at a regional scale using multitemporal LandsatTM data. AMBIO 2003, 32, 502–509. [Google Scholar]
  38. Purevdorj, T.; Tateishi, R.; Ishiyama, T.; Honda, Y. Relationships between percent vegetation cover and vegetation indices. Int. J. Rem. Sens. 1998, 19, 3519–3535. [Google Scholar] [CrossRef]
  39. Land-Use Change Survey Statistics of Guangdong Province. Available online: http://www.gdlr.gov.cn/newsAction.do?method=viewNews&classId=020019980000000488&newsId=020010040000014133 (accessed on 12 December 2012).
  40. Zheng, J.; Yin, Y.; Li, B. A new scheme for climate regionalization in China. Acta Meteorolog. Sinic. 2010, 65, 3–12. [Google Scholar]
  41. Zhao, S.-Q. A new scheme for comprehensive physical regionalization in China. Acta Geogr. Sinic. 1983, 38, 1–10. [Google Scholar]
  42. Li, X.-B.; Shi, P.-J. Sensitivity analysis of variation in NDVI, temperature and precipitation in typical vegetation types across China. Acta Phytoecolog. Sinic. 2000, 24, 379–382. [Google Scholar]
  43. Levin, A.; Lin, C.F.; Chu, C.S.J. Unit root tests in panel data: Asymptotic and finite-sample properties. J. Econometrics 2002, 108, 1–24. [Google Scholar] [CrossRef]
  44. Maddala, G.S.; Wu, S.W. A comparative study of unit root tests with panel data and a new simple test. Oxford B Econ. Stat. 1999, 61, 631–652. [Google Scholar]
  45. Pedroni, P. Panel cointegration: Asymptotic and finite sample properties of pooled time series tests with an application to the PPP hypothesis. Economet. Theor. 2004, 20, 597–625. [Google Scholar]
  46. Kao, C. Spurious regression and residual-based tests for cointegration in panel data. J. Econometrics 1999, 90, 1–44. [Google Scholar]
  47. Larsson, R.; Lyhagen, J.; Löthgren, M. Likelihood-based cointegration tests in heterogeneous panels. Econometrics J. 2001, 4, 109–142. [Google Scholar] [CrossRef]
  48. Fuller, R.A.; Gaston, K.J. The scaling of green space coverage in European cities. Biol. Lett. 2009, 5, 352–355. [Google Scholar]
  49. Grossman, G.M.; Krueger, A.B. Environmental Impacts of a North American Free Trade Agreement. Available online: http://www.nber.org/papers/w3914 (accessed on 12 December 2012).
  50. Mohan, M.; Pathan, S.K.; Narendrareddy, K.; Kandya, A.; Pandey, S. Dynamics of urbanization and its impact on land-use/land-cover: A case study of megacity delhi. J. Environ. Protec. 2011, 2, 1274–1283. [Google Scholar] [CrossRef]
  51. Piao, S.; Fang, J.; Zhou, L.; Guo, Q.; Henderson, M.; Ji, W.; Li, Y.; Tao, S. Interannual variations of monthly and seasonal normalized difference vegetation index (NDVI) in China from 1982 to 1999. J. Geophys. Res. 2003, 108. [Google Scholar] [CrossRef]
  52. Dewan, A.M.; Yamaguchi, Y. Land use and land cover change in Greater Dhaka, Bangladesh: Using remote sensing to promote sustainable urbanization. Appl. Geogr. 2009, 29, 390–401. [Google Scholar] [CrossRef]
  53. Granger, C.W.J. Ome recent development in a concept of causality. J. Econometrics 1988, 39, 199–211. [Google Scholar]

Share and Cite

MDPI and ACS Style

Li, C.; Kuang, Y.; Huang, N.; Zhang, C. The Long-Term Relationship between Population Growth and Vegetation Cover: An Empirical Analysis Based on the Panel Data of 21 Cities in Guangdong Province, China. Int. J. Environ. Res. Public Health 2013, 10, 660-677. https://doi.org/10.3390/ijerph10020660

AMA Style

Li C, Kuang Y, Huang N, Zhang C. The Long-Term Relationship between Population Growth and Vegetation Cover: An Empirical Analysis Based on the Panel Data of 21 Cities in Guangdong Province, China. International Journal of Environmental Research and Public Health. 2013; 10(2):660-677. https://doi.org/10.3390/ijerph10020660

Chicago/Turabian Style

Li, Chao, Yaoqiu Kuang, Ningsheng Huang, and Chao Zhang. 2013. "The Long-Term Relationship between Population Growth and Vegetation Cover: An Empirical Analysis Based on the Panel Data of 21 Cities in Guangdong Province, China" International Journal of Environmental Research and Public Health 10, no. 2: 660-677. https://doi.org/10.3390/ijerph10020660

APA Style

Li, C., Kuang, Y., Huang, N., & Zhang, C. (2013). The Long-Term Relationship between Population Growth and Vegetation Cover: An Empirical Analysis Based on the Panel Data of 21 Cities in Guangdong Province, China. International Journal of Environmental Research and Public Health, 10(2), 660-677. https://doi.org/10.3390/ijerph10020660

Article Metrics

Back to TopTop