Vegetation change detection research of Dunhuang city based on GF-1 data

This study selected the Dunhuang city with the unique landscape pattern, its oasis surround by desert, for studying the vegetation change. Based on two period GaoFen-1 images, combined with field survey data, the dimidiate pixel model and fractional vegetation coverage estimation model were applied to detect the changes of local vegetation coverage from July 2013 to July 2015. Analysis showed that during the study years, the Grain for Green project had a remarkable effect on the changes. The vegetation types of intermediate high and very low grades increased by 0.66 and 6.78 km2 respectively. The low vegetation coverage decreased by 23.87 km2. The vegetation coverage types of intermediate grade increased by 88.97 km2 because of the planted forest, which accounted for 10.84% of the study area. The vegetation coverage types of high grade reduced by 72.47 km2. This change effectively prevents the spread of desert and lays a good foundation for the ecological construction of Dunhuang city.


Introduction
Oasis refers to a green land with weeds in the desert. Dunhuang oasis, is one of the main bases for human settlement and production in the arid and semi-arid regions of Northwest China. Fractional Vegetation Cover (FVC) is the basic index of ecological system, which is affected by the changes of natural environment and human activities. In order to understand the current situation, trends, driving forces, and environmental impacts of the vegetation cover in desert oasis, using remote sensing date, numerous studies have been performed (Yang et al. 2010;Liu et al. 2013;Zhang et al. 2015). However, there are few research based on the Chinese high resolution satellite constellation. Based on two period GaoFen-1 (GF-1) images, the dimidiate pixel model and vegetation coverage estimation model were applied to detect the changes of vegetation covers from 2013 to 2015 in Dunhuang oasis. The goal is of this study is to provide a reference for accurately understanding the situation of regional ecological change and evaluating the effectiveness of the ecological projects.

Study area and data
Dunhuang city of the northwest of Gansu province, is located at the ancient Silk Road. The territory is surrounded by mountains at three aspects. The Danghe river is the only water source for the survival of the city. GF-1 images (two for July 2013 and two for July 2015) covering Dunhuang city were obtained from China Center for Resources Satellite Data and Application. GF-1 is a kind of high resolution earth observation satellite with independent & Xin Tian tianxin@caf.ac.cn property rights in China (Fig. 1). Data pre-processing was conducted by the ENVI software (version 5.4.1).

Image data preprocessing
The RPC file of GF-1 images in 2013 and the corresponding DEM from Geospatial Data Cloud were used to make the geographic correction. After that, the GF-1 images in 2015 were processed to match with the GF-1 images from 2013 by means of an image-to-image matching method provided by the ENVI Imagine software. The root mean square (RMS) control of each control point is less than 1 pixels, and the total root mean square error (RMSE) is 0.49 (Fig. 2).

NDVI and vegetation coverage
The normalized difference vegetation index is known to be an index of vegetation conditions (Carlson and Ripley 1997). This index is expressed as: NIR and R represent the reflectivity of the near infrared and red bands, respectively (Ye et al. 2009). Combined with the content of the study, the NDVI gray scale images of different time images are obtained by using NDVI. It is necessary to transform the vegetation index into the degree of coverage. In a certain extent, FVC can make up for the disadvantage of NDVI, in which NDVI is difficult to distinguish the low coverage vegetation and the coverage is easy to be saturated. On the contrary, so as to effectively stretch the range of vegetation information. The vegetation cover can be divided into two types, pure and mixed pixels by the different structure characteristics. A mixed pixel includes two kinds of information of soil and vegetation. Assume that the vegetation index is NDVIveg, and the vegetation index is NDVIsoil, and the mixed pixel vegetation index is: The vegetation coverage model is: In the above formula, NDVIveg represents the normalized difference vegetation index of the whole vegetation cover, NDVIsoil represents the normalized difference vegetation index of bare soil.
The key of the mixed pixel decomposition method is the acquisition of NDVIveg and NDVIsoil. In this paper, NDVIsoil and NDVIveg parameters were generated by using the NDVI threshold of land cover type. The following expressions are used to calculate the NDVIsoil and NDVIveg: where b1, b2, b3, b4 and b5, are forest land cover document, agricultural land cover document, urban land use mask files, water Mask documents and other land use mask file respectively; n1 and n6, n2 and n7, n3 and n8, n4 and n9, n5 and n10, are the NDVI minimum and maximum value of corresponding each mask file respectively expressed. The Eqs. 2-3 were applied to conduct the transformation of gray scale map of vegetation cover (Zhang et al. 2017).

Classification of vegetation coverage
According to the criteria of classification in desertification survey of China, and the characteristics of Dunhuang oasis vegetations, the coverage vegetation in the study area were divided into five grades, (1) very low coverage vegetation,  Table 1. Normally, the transformation of the low level vegetation to the high level is called the vegetation growth, and the transformation of the higher vegetation to the lower one is named as the vegetation degradation. According to the above classification criteria, the vegetation coverage map of the different periods is obtained (Fig. 3). Through field survey and land use map can know vegetation classification results is right. Then, the vegetation area and its changes from 2013 to 2015 were obtained by difference operation.  5 Analysis on the factors of vegetation cover change

Natural factors
Natural factors were mainly analyzed from two aspects, namely, precipitation and vegetation. Dunhuang city is characterized by little precipitation but large evaporation. The average amount of evaporation is 60 times of precipitation, and there almost no conversion among precipitation, soil surface water and ground water. Precipitation is one of the most important factors affecting vegetation coverage. The precipitation from 2013 to 2015 in Dunhuang city is shown in Table 5, which showed that both the annual precipitation and the summer precipitation in 2013 and 2015 were very small. As a result, the water supply for agriculture and forestry mainly came from groundwater. Water stress is the mainly problem of Dunhuang city and because of the different farming conditions, the local areas change greatly from 2013 to 2015.  Grassland and desert are the main vegetation types of the Dunhuang city. The main types of plants are Populus euphratica, Tamarisk and Manaplant Alhagi, which can resistant to drought, wind and sand and adapt to extreme drought climate conditions. In particular, the oasis fringe area, its vegetation types are more barren, and the vegetation species structure is simple with extremely low FVC. Which is consistent with that shown in Fig. 3, the oasis surrounding is mainly the extremely low vegetation coverage.

Social factors
In combination with various factors, social factors were also main factors affecting the vegetation cover change in Dunhuang oasis. Social factors were mainly analyzed from these aspects, namely, population, economic, government industrial structure conduction and so on.
As the artificial forests increase, the effects of Grain for Green project are significant. According to the information from official website of the Forestry Bureau of the Dunhuang City. In recent years, the amount of afforestation in Dunhuang is 7.12 million hectares, 14.5 million hectares of afforestation from the project 0, and 160 million acres natural vegetation were protected and managed. Only in 2015, the afforestation area is beyond 453 hectares, 1180 hectares of economic forests being newly planted, and 91.2 km of farmland shelterbelts were built and improved. As a result of large amounts of farmlands were converted to artificial forests, the medium vegetation coverage increased largely and the high vegetation coverage dereased.
Dunhuang was named as national comprehensive demonstration area of desertification prevention and control in 2013. 190 hectares sands were got under control in Mingsha Mountain Montenegro, Mogao Grottoes and other *0-0.00001 is a special treatment for the classification, in order to remove background values and not allow them to participate in the classification With the rapid development of economy, the local industrial structure has been constantly adjusted. The Dunhuang City is an indispensable part for the ''Silk Road Economic Zone'' of ''the Belt and Road'' (B&R). As a result, ''the Belt and Road'' strategy not only brought more opportunities for the development, but also promoted its industrial restructuring, especially for the service industry. With the optimization and upgrading of industrial structure, the service industry, which is dominated by tourism, has became to be the main force for city's economic development. The income from the tourism and cultural industries is shown in Fig. 4.
During the period of the China 12th Five-Year Plan, which is scheduled from 2011 to 2015, the three industries of the city were fully developed, with the primary, second and service industries increasing by 3.7, 18.1 and 15.5%, respectively. These three industrial structures tend to be optimized from 21.  (2015). The proportion of the first and secondary industry decreased by 2.9% and 2.0%, the proportion of the tertiary industry increased by 4.9%. The city's three industrial structures were shown in Table 6.
With a slight increase of the local population, the employment structure was gradually optimized. The census data published by Dunhuang City Government Network showed that by the end of 2012 and the end of 2016, the number of household population is 141,969 and 143,800, respectively. However, with the gradual optimization of the employment structure, the labor force has been transferred from the first and second industries to the tertiary industry, and as a result, some of them moved from rural to urban areas. According to the 12th Five-Year Plan Economic and Social Development Review for Dunhuang city, published by the Dunhuang government network, the employment structure of the three industries in this city changed from 41. 8:15.4:42.8 in 2010 to 37.8:14.9:47.3 in 2015. The reduction of labor force in the primary industry has a direct impact on the reduction of cultivated land area. This is also

Conclusions
In this study, based on GF-1 NDVI information and dimidiate pixel model, changes of vegetation coverage between 2013 and 2015 in the Dunhuang oasis were analyzed. The results showed that the implement of Grain for Green project has largely affected the changes. According to the changing static information, the middle and high vegetation coverage, and extremely low vegetation coverage increased 0.66 and 6.78 km 2 , respectively, and low vegetation coverage reduced by 23.87 km 2 , and the intermediate vegetation coverage increased by 88.97 km 2 , and high vegetation coverage area decreased by 72.47 km 2 . In summary, the social factors are the most significant factors affecting the changes of the Dunhuang. In the light of the present vegetation coverage situation of the Dunhuang oasis, the protection work for the oasis should be continually conducted in order to keep the oasis stability. Meanwhile, the local water saving agriculture and Grain for Green should also be further developed.