• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Characterization of landscape pattern based on land economic niche change: A case study in Ganzhou,Gansu Province, China

    2018-07-04 10:05:36HuaLiTongPeiJiShiXueBinZhangZaiYanLi
    Sciences in Cold and Arid Regions 2018年3期

    HuaLi Tong, PeiJi Shi , XueBin Zhang, ZaiYan Li

    College of Geography and Environmental Science, Northwest Normal University, Lanzhou, Gansu 730070, China

    1 Introduction

    Landscape ecology is an interdisciplinary between geography and ecology, which was developed in the 1980s. As the core research contents, spatial statistical characteristics comparison (Li, 2000; Qiet al.,2013), landscape metrics (Wu, 2000; Wuet al., 2002)and Markov transfer matrix (Fenget al., 2010) have become the main methods and indicators for landscape spatial analysis. Landscape metrics are widely used in landscape pattern analysis and dynamic research (Fan and Ding, 2016; Robinson and Weckworth, 2016), which contain highly concentrated landscape information, and simply reflect the landscape structure and spatial allocation (Zhanget al., 2003).With the development of GIS and the upgrading of Fragstats (Chenet al., 2002), research on landscape metrics continues to develop, such as patch shape index, landscape diversity, fractal dimension, fragmentation, grain degree and other quantitative research landscape indices (Pennocket al., 1994; Wuet al.,2002; Fuet al., 2011). In terms of landscape representation, landscape metrics can reveal landscape patterns such as the evolution and alteration of land-scapes (Verburget al., 1999; Freitaset al., 2013), but does not reveal the driving and influencing factors (Liet al., 2004). Therefore, landscape metrics has become a key issue in understanding the driving and influencing factors that cause landscape changes at a regional scale (Chaplotet al., 2005; Diggelenet al.,2005; Bajoccoet al., 2016).

    A niche refers to the position and function of a species or population in the ecological environment,that is, the position and interdependence temporally and spatially (Zhang and Xie, 1997; Silvertown,2004). With the development of niche theory, some scholars have proposed natural, social and economic complex ecosystems to apply niche theory to a wider range of fields (Basilleet al., 2008; Nascimentoet al.,2010; Cardosoet al., 2017). In recent years, the theory has a close connection with the study of land resources (Haoet al., 2010; Ceccarelliet al., 2014; Huet al., 2016). Some scholars have constructed the land niche on the basis of the ecological niche theory and model (Zhanget al., 2002; Dan and Seifert, 2011).The land niche can reflect space, function and position of each land type in the land use system, and quantitatively express the ability of each land type to occupy or utilize a new habitat. At the spatial level,the land niche is a measurement of the interaction between different geographical units (Qiet al., 2007;Wu and Li, 2014; Huet al., 2016). The larger the land niche, the stronger the ability of the regional spatial unit to attract population and material circulation.

    The change of land use type driven by land niche has a direct or indirect effect on regional material circulation and energy flow, which is one of the important factors of regional land landscape pattern change(Yuet al., 2015). Therefore, it is of great significance to characterize the regional spatial landscape pattern using the land niche theory. Based on the ecological significance of landscape pattern on a county scale,this paper used Ganzhou District as the study area.We selected six landscape indices that characterize landscape quantity, structure and morphological characteristics and human disturbance, which clarifies the trend of land use change and explores the characterization of landscape pattern based on land economic niche.

    2 Study area and material

    2.1 Study area

    Ganzhou District is a county level district, administered by Zhangye, a prefecture-level city, and Ganzhou is the city proper of Zhangye. This county is located in the middle of the Hexi Corridor, between 100°04'E-100°52'E, 38°32'N-39°24'N. It is adjacent to Shandan and Minle counties to the east, Sunan Yugur Autonomous County to the south, Linze County to the west, and Inner Mongolia Autonomous Region to the north. The county administers a national economic and technological development zone, 18 townships, 5 street offices and 245 administrative villages. Heihe River, which traverses the county, is the second largest inland river in China, and forms a wellknown oasis known as "Sai Shang Jiang Nan", meaning Lush southern-type fields on the frontier. Ganzhou is 65 km long from east to west and 98 km from north to south. The total area is 3,657 km2.

    The climate in Ganzhou District is typical continental arid. There is scarce precipitation and strong evaporation and large daily temperature difference.The north side are Heli Mountain and Longshou Mountain, and south sides is Qilian Mountain, and the terrain tilts from southwest to northeast. Because of soil and geographical differences, Ganzhou District contains different types of landscape, including mountain forest, plain forests, desert shrub grassland,meadow, marsh and wetlands. The Sai Shang Jiang Nan oasis is formed from five major rivers, the Heihe,Suyoukou, Dayekou, Daciyao and Shandan. Annual runoff is 2.46 billion m3.

    Ganzhou District is a commodity grain base in Hexi Corridor and a tourist city along the Silk Road(Tonget al., 2015; Tonget al., 2018). It has a total population of over 513.6 thousand (2015), 48% of whom are urban dwellers. The total GDP of Zhangye City reached 15.675 billion Yuan (approximately 2.58 billion US Dollars) in 2015. Ganzhou District has a balanced industrial structure, with a ratio of primary:secondary: tertiary industry of 22.7:22.2:53.1 in 2015.

    2.2 Data and methodology

    2.2.1 Data process

    Data for this study include satellite remote sensing images of Landsat TM and ETM, and Ganzhou District Statistical Yearbook for 1995, 2000, 2005,2010 and 2015 (Ganzhou Statistics Bureau,2005-2015). Remote sensing data was obtained from the USGS website (http://glovis.usgs.gov), the images were mainly obtained in August and September.Vegetation characteristics are more robust during this period, which is conducive to classification and interpretation. Auxiliary data is from the second national land survey data of Ganzhou District. Based on the landscape classification system, the remote sensing image (30 m) of the study area was supervised by Envi5.0 software. Land use classification was obtained by geometric correction, cloud removal and image enhancement. On the basis of this classification, the land use map, Google Earth high resolution images, and field investigation results of different periods were used to correct interpretation results,and finally the land use data of Ganzhou District were obtained. Through the verification of these results,overall accuracy reached more than 90%, which meet the accuracy requirements of this study. Secondly,interpretation of the land use map was converted into a grid format in Arcmap10.2, using the Fragstasts 4.2 software. Thus, the landscape metrics can be calculated.

    Figure 1 Location of study area

    According to the classification of land use status(GB/T 21010-2007) and the classification method of land use/cover classification system of remote sensing in China, combined with the landscape characteristics of the study area, Ganzhou District is divided into six landscape types: farm land, forest land, grassland, construction land, watershed and bare land. The watershed includes rivers, lakes and wetlands. The construction land includes urban, rural residential and mining.

    2.2.2 Methodology

    (1) Land economic niche

    Niche is the relative position and function of the biological unit in the process of interacting with the environment in a particular ecosystem (Zhu, 1997).Mainly, it includes two aspects: first, the state of biological units (biomass, the number of individuals and biological units, energy, adaptability, resource share,level of economic development), which is the result of past growth and development, socio-economic development and the interaction with the environment, that is, the niche state. Second is the influence or dominance of biological units, such as the rate of energy,growth, economic growth and material transformation, the ability to occupy a new habitat, and productivity, that is, niche potential. The land ecology niche is a natural, social and economic multidimensional concept, but under the market economy condition, its main function is the economic niche (Qin and Gao, 2004; Wu and Li, 2014). The land economic niche is mainly composed of the comparative economic benefits of land use, and the status, effect and function of the land use ecological element in the process of regional economic sustainable development(Chenet al., 2010). The formula for the land economic niche is as follows (Xiaoet al., 1997):where,Niis the land economic niche of different land use types;Siis the state of different land use types andPiis the potential of different land use types;Aiis the dimension conversion coefficient;MiXiis the product value of land area and the unit area income, which is total revenue or the GDP of theithtype land use,Yiis the annual growth rate of the economic income or GDP of theithtype land;Riis the population attached to theithtype land.

    The economic benefits of the land type used in this paper are obtained from the statistical yearbook.The farmland output value is the benefit of agricultural in the year. The forest production value is the forestry benefit and the grassland output value mean the grassland benefit. The watershed output value is the benefit of the aquatic products. The construction land output value is the total benefit of second and tertiary industries.

    (2) Landscape metrics

    Landscape metrics is the most widely used method for landscape spatial pattern analysis. Based on previous studies, this paper selected the more commonly used and meaningful indicators (Qiet al.,2013), including total area of patches (CA), number of patches (NA), the mean patch area (AREA_MN),the aggregation index (AI), the fragmentation index(FN) and the fractal dimension index (FD). The calculation of landscape metrics is calculated by Fragstats 4.2 software. FN is calculated by the formula based on the data calculation of Fragstats 4.2 software. The aforementioned index can represent the characteristics and change of spatial pattern in Ganzhou District from four aspects.

    1) Quantity characteristics. CA determines landscape extent and the maximum scale of research and analysis, which is the basis for calculating other indicators. NP is the number of land use type patches that can be used to describe the heterogeneity of the entire landscape. AREA_MN is the average area of a patch type, which represents an average condition. The value change can reflect the rich landscape ecological information, which is the key to reflect landscape heterogeneity.

    2) Morphological characteristics. FD is a measurement of complex irregularity, which is between 1 and 2, and can quantitatively describe core area size of a certain type and the meandering/sinuosity of its boundary line. If FD approaches 1, patch geometry tends to be more regular; conversely, patch geometry tends to be irregular. This is mainly to reflect the validity of the complex form of space, which also reflects the impact of human activity on the landscape pattern (Krummelet al., 1987).

    3) Structure. AI indicates the adjacent relationship between similar types of patch pixels from the landscape level. When there is no common boundary between all the pixels in a certain type, the aggregation index is the lowest. Conversely, the larger the common boundary, the greater the aggregation index(maximum value is 100).

    4) Human disturbance. FN is used to describe landscape fragmentation. This refers to the process of a plaque mosaic from single, homogeneous and continuous to complex, heterogeneous and discontinuous individual under natural or human disturbance. The formula is as follows (Li, 1989):

    where, MPS is the mean patch area of all kinds of patches in the landscape,nis the total number of patches of the landscape type, andAis the landscape area.

    3 Results and analysis

    3.1 The change of land economic niche in Ganzhou

    The land economic niche reflects the function of land cover in Ganzhou from the economic output value. The landscape pattern can express the spatial distribution of land cover, and the landscape metrics reflects the spatial structural of land cover. The land economic niche combined with landscape pattern and metrics indicate the ecological possession in function,space and structure of land cover. According to the Formula(1), the land economic niche of Ganzhou District from 1995 to 2015 is calculated (Table 1).

    As presented in Table 1, during 1995-2015,niches of farmland, watershed and construction land in Ganzhou show a rising trend. The economic niche of cultivated land rose from 0.11 to 0.14, with an average annual growth rate of 1.2%. The economic niche of watershed increased from 0.07 in 1995 to 0.14 in 2015, with an average annual increase of 3.5%. In these five types of land cover, the annual growth rate of watershed economic niche is the highest. The study area is located in the northwestern arid region whose ecology is very fragile, with most of the output value of water lies in its indirect creation of economic value and protection of the regional ecological elements to guarantee normal operation and survival. Construction land niche rose from 0.34 in 1995 to 0.65 in 2015, with an average annual growth of 3.2%. The economic growth of forest land first dropped then rose to a peak of 0.06 in 2005, then decline again from 2005 to 2015. The overall trend of forest was declining, with an annual decline rate of 4.5%. The decline of grassland was the largest, from 0.43 in 1995 to 0.05 in 2015, with an annual decrease of 10.2%. This indicates that the grassland economic niche in the study area accounted for the most important position before 2000. With rapid urbanization,construction land gradually replaced grassland. The construction economic niche has become the largest since 2000. The construction land output value has the greatest contribution to the economic development of Ganzhou District. This is because China's rapid economic development led to rapid urbanization over the past 20 years.

    Table 1 Land economic ecological niche of landscape from 1995 to 2015 in Ganzhou District

    The general change of land economy ecological niche is as follows. The economic niche of construction land and watershed are basically increasing while grassland is decreasing. The economic niche of farmland shows a slight upward trend while grassland shows a dramatic decreasing trend. Combined with the trend of local land use change, the regional production value and urbanization level, the output value of the first industry in Ganzhou District show a downward trend, the second and third industries developed rapidly, the output value became high, resulting in change in the development model. Second and tertiary industries became the main source of income in Ganzhou District.

    3.2 Spatial analysis of landscape pattern change in Ganzhou District

    According to the local circumstances and the national "land use classification" (GB/T 21010-2007),the study area is divided into six categories which are farmland, forest, grassland, construction land, watershed, and bare land (Figure 2) through ENVI 5.0.Aiming at the acquired five phases of remote sensing image data, classification was carried out by a combination of visual interpretation and supervised classification. The total accuracy of the five phases of data after image classification is 90.6%, 94.2%, 93.8%,94.6% and 95.4%, respectively. The Kappa coefficients are 0.8709, 0.9326, 0.9209, 0.9338 and 0.9423,respectively, which meet the requirements of research and provides reliable data guarantee for later comprehensive analysis.

    As presented in Figure 2, during the study period,there was an obvious increase of farmland and construction land, while grassland and bare land continuously decreased. The degradation of grassland and conversion to farmland are obvious, which is the embodiment of large-scale agriculture and the result of increased human activity.

    3.3 The change of landscape metrics

    3.3.1 The change of landscape type

    As presented in Table 2, farmland is the main land use type in the study area, followed by grassland, watershed and construction land. Farmland, construction land and watershed areas increased continuously from 1995 to 2015. Among which, the CA of farmland increased from 895.625 km2to 1,120.500 km2, with an increase of 224.875 km2. NP of farmland decreased,AREA_MN and AI increased, and FN became smaller. CA of construction land increased from 99.94 km2to 146 km2, with an increase of 46.46 km2. NP of construction land decreased, AREA_MN and AI increased, and FN became smaller. CA of forest land changed slightly, with a reduction of 2.81 km2. CA of the watershed decreased from 1995 to 2005, but increased between 2005 and 2015. This increase is due to local government's implementation of the wetland protection policy to restore wetlands and the building of wetland parks and three artificial lakes in urban areas after 2006. Overall, the water area shows a decreasing trend generally from 135.63 km2to 123.38 km2, with a decrease of 12.25 km2. The main source of landscape change in Ganzhou District is grassland and bare land. From 1995 to 2015, the grassland decreased from 794.94 km2to 680.19 km2,a reduction of 114.75 km2. The bare land decreased from 1,647.5 km2to 1,506.1 km2, a reduction of 141.4 km2. NP of grassland and bare land shows an increasing trend while AREA_MN shows a decreasing trend.

    Figure 2 The spatial pattern of land use in Ganzhou District from 1995-2015

    3.3.2 The change of landscape morphology

    The FD absolute value change during 1995-2015(Table 2) for the watershed is the highest, followed by grassland, farmland, forest and construction land.According to the landscape ecology theory, the complex patch morphology indicates a high degree of naturalization and low self-similarity. The FD of grassland and construction land are relatively low, indicating that the patch of grassland and construction land have high self-similarity, and the patch morphology is regular, which is strongly influenced by human activity. Due to the dynamic change of the landscape's FD, the FD of farmland and construction land decreased from 1.562 to 1.544 and from 1.531 to 1.511,respectively. This decrease demonstrates how human activity consolidated the scope and shape of farmland,and construction land was programmatically expanded during the process of urbanization. The FD of the watershed increased from 1.702 to 1.719 due to local governmental implementation of the wetland restoration and protection policy and the Heihe River Basin program. The FD of forest and grassland shows an increasing trend. The FD of forest and grassland increased from 1.491 to 1.549 and from 1.571 to 1.584, respectively. This increase is due to overgrazing and land degradation, with portions of forest and grassland converted to farmland and construction land.

    3.3.3 The change of landscape structure

    Landscape structure is the spatial arrangement and composition of landscape components and elements.In landscape ecology, AI indicates the adjacency between pixels of the same type patch at a landscape level. As presented in Table 2, AI of farmland patch is the highest except for bare land, followed by grassland and forest, while watershed and construction land is relatively low. During 1995-2015, the AI of farmland, watershed and construction land show an increasing trend while forest, grassland and bare land show a decreasing trend. The AI of farmland increased from 79.75 to 82.68 (4% increase), watershed increased from 34.22 to 42.42 (23.94% increase), and construction land increased from 25.95 to 36.55 (41%increase). The AI of forest decreased from 60.55 to 55.89 (7.7% decrease), and grassland decreased from 64.59 to 57.99 (10% decrease). Combined with Figure 2, we can see that the AI of farmland, watershed and construction land increased because their total area increased. The reduction of forest and grassland AI are closely related to the reduction of the corresponding total area. Portions of forest land and grassland area were converted into cultivated land and construction land. Also, overgrazing and degradation are the main reasons for the reduction of forest and grassland areas.

    Table 2 Landscape indices of different land use types of Ganzhou District in 1995, 2000, 2005, 2010 and 2015

    3.3.4 The impact of human activity

    From the absolute value of landscape FN for Ganzhou District during 1995-2015, construction land is the largest; watershed is second, followed by forest,grassland, farmland and bare land. In landscape ecology, landscape structure of construction land is simple, the degree of homogenization is high, and has the highest land economic niche. Construction land is the most affected by human activity and has the highest FN in the study area. The study area is in an arid region, located in a fragile ecological zone in Northwest China. The total amount of water resources is unstable, and coupled with human interference, wetland degradation is serious, and the FN of watershed is high. The few natural forest is mostly located near water sources. In recent years, sheltered forest area has increased slightly due to China's national policy, but the distribution is heterogeneous and dependent on human activity. The spatial landscape of grassland is complex, heterogeneity is high, and land economic ecology is less than construction land and farmland. The FN of farmland and bare land are relatively low. Farmland is affected by terrain, water, soil quality and human activity. In order to facilitate cultivation, farmland in the study area is mostly large and continuous. The FN of bare land is low because most of the landscape in Ganzhou District is bare rock and desert, where conversion to other land types is difficult, and human activity has little influence on bare land. However, in recent years, numerous improved land and large-scale renovation projects have been started in Ganzhou, such as "Wa Qu Pai Jian"which means "Digging canals and removing alkali Project", and "Farming Plan in Gebi Desert". In order to improve the regional ecological environment, a new urban plan will be implemented, expanding west into barren lands.

    3.4 The relationship between land economic niche and landscape metrics

    In order to study the relationship between the land economic niche and corresponding landscape metrics,the Pearson correlation coefficient was obtained by using SPSS (Table 3).R1jis the correlation coefficient between the farmland economic niche and landscape metrics.R2jis the correlation coefficient between the forest economic niche and landscape metrics.R3jis the correlation coefficient between the grassland economic niche and landscape metrics.R4jis correlation coefficient between the watershed economic niche and landscape metrics.R5jis the correlation coefficient between the construction land economic niche and landscape metrics.

    As presented in Table 3, CA, NP, AREA_MN, AI,FN and FD of farmland and grassland are all correlated with the land economic niche. The watershed and forest have no significant correlation between the corresponding landscape metrics and the land economic niche. It can be found that NP and FN are negatively correlated with the land economic niche, and AI is positively correlated with the land economic niche. However, the NP of construction land is positively correlated with the land economic niche. Because the study area belongs to an oasis area,and some urban area is also an important ecological protection area.

    Table 3 Correlation coefficients between land economic ecological niche and landscape indices

    4 Discussion

    The theory of land economic niche combines natural and human factors, focusing on land use economic benefits, which describe the contribution of different land use to the local economy. Also, regional economic development is an important driving force for the change of landscape pattern. Results of the correlation between the land economic niche and corresponding landscape metrics show that partial landscape metrics has a significant correlation with the land with large direct economic benefits, such as farmland, grassland and construction land. The land economic niche is negatively correlated with NP, FN and FD, and positively correlated with AI. Correlation analysis of the land economic niche and corresponding landscape metrics shows that the land economic niche is only relevant to NP, FN, FD and AI of farmland, grassland and construction land, which are affected by change of the land economic niche. The economic niche of forest and watershed has no significant correlation with the corresponding landscape metrics. On the one hand, watershed and forest in the arid area does not increase more direct economic benefits. Furthermore, watershed and forest are determined by natural properties. Water is limited by weather and human factors, and the growth period of trees is long. In general, the change of land economic niche affects the change in regional landscape patterns.

    In agriculture, farmland area increased year by year, NP, FN and FD decreased with the increase of farmland economic niche, and AI increased, which is consistent with the correlation analysis. In the second and third industry, the construction land economic niche increased year by year, CA, NP and AI increased, FN and FD decreased with niche increase,which is consistent with the correlation analysis except for CA. In animal husbandry, the grassland economic niche show a decreasing trend, CA decreased,NP show an increasing trend, and AI show a decreasing trend which is consistent with the conclusions of the correlation analysis. Farmland in the study area has been developed over a large area, and construction land is planned to expand from the original oasis area. Some of the grassland has converted to cultivated land and construction land. Although there is a small amount of farmland and bare land converted to grassland, the grassland area is still shrinking. Thus,there is a serious ecological security threat of the Ganzhou oasis area. The watershed area has decreased,NP, FN and FD increased with the land economic niche decrease, which is consistent with the correlation analysis. Since Ganzhou District is in an arid ecological fragile area in northwest China, the wetland and Heihe River play an important role in maintaining regional ecological balance as the largest water body in the study area.

    In this study, although the economic niche of forest and watershed has no significant correlation with the corresponding landscape metrics, the reduction of forest and watershed area will inevitably lead to a decrease of regional water conservation function.Also, the increase and intensive use of farmland will aggravate agricultural water use. The shortage of water resources leads to the inability of anthropogenic development of bare land in the study area, and limited human activity leads to limited economic development. Even cause a series of problems such as ecological deterioration. The local government has made some effort to solve these problems where the watershed area has increased since 2005 and watershed economic niche also increased. Other measures need to be implemented such as converting low-quality farmland to grassland and forest which can resist salt intrusion and drought. Based on the regional farmland red line barren land should be converted into grassland and forest. At the same time, high-quality farmland is strictly restricted for other purposes. For bare lands around the oasis, we can improve the soil and plan grass and forest.

    Comparing the economic niche of grassland and construction land in 2010 and 2015, grassland was reduced by 58.3%, and construction land increased by 14%, which indicate the economic benefits of rapid development of second and tertiary industries accounted for a larger proportion. Also, production output of second and tertiary industries will have a certain impact on the agricultural economy. Under the current situation, how to adjust the regional industrial structure is still an important issue.

    5 Conclusions

    This paper introduced the theory and method of land economic niche into landscape ecology, which provide a new method for spatial characterization of urban and rural spatial landscape patterns. Our case study contributes to the literature and practice of characterization of landscape pattern by calculating land economic niche. Change of land economy niche can mark the direction of economic development of Ganzhou District. The general characteristic of the land economic niche in Ganzhou District is as follows,economic niche of construction land, watershed and farmland increased during 1995-2015, while grassland declined significantly. Combined with local land use change trend, regional GDP and urbanization level, the output value of the primary industry in Ganzhou District shows a decreasing trend. Also, second and third industries developed rapidly and the output value was high, which resulted in the development model change. The area of local construction land continued to grow, and the second and third industries became the major economic income in Ganzhou District.

    The change of land economic niche has a driving effect on the change of landscape pattern on a county scale. Economic development and the acceleration of urbanization process results in a change of local landscape type under the influence of human activity.Economic development led to changes in land cover.During 1995-2015, land cover changed greatly in Ganzhou District, where the conversion of grassland and bare land to farmland was obvious. The economic niche of farmland, construction land, watershed and grassland show a negative correlation with NP, FN and FD, and had a positive correlation with AI. There is no significant correlation between the forest land economic niche and landscape metrics. Also, the landscape pattern of forest land changed due to the conversion of cultivated land, construction land, grassland and land degradation. The landscape metrics also changed accordingly.

    The land economic niche can characterize the direction of regional human activity. The pattern of landscape changed greatly from 1995 to 2015 in Ganzhou District, and landscape heterogeneity show an increasing trend. Also, land use is developing in the direction of diversification and homogenization. In particular, the landscape structure in Ganzhou District in 2015 was complex. The landscape pattern in the study area is becoming more and more complicated, and the human disturbance has increased. The increasing areas of farmland and construction land are due to human activity. The farmland is developing in a large continuous patch while construction land is planned to expand from the original oasis area programmatically.The watershed economic niche increased year by year while the forest economic niche is small and has decreased. At the same time, the economic niche and area of the grassland decreased, and the degree of fragmentation is high.

    Acknowledgments:

    This work is funded by National Natural Science Foundation of China (Grant No. 41771130 and Grant No. 41661035). We highly appreciate the editor and the two reviewers for their valuable comments,suggestions and editing, which have greatly helped to improve this manuscript.

    :

    Bajocco S, Ceccarelli T, Smiraglia D,et al., 2016. Modeling the ecological niche of long-term land use changes: the role of biophysical factors. Ecological Indicators, 60: 231-236. DOI: 10.1016/j.ecolind.2015.06.034.

    Basille M, Calenge C, Marboutin E,et al., 2008. Assessing habitat selection using multivariate statistics: some refinements of the ecological-niche factor analysis. Ecological Modelling, 211(1-2):233-240. DOI: 10.1016/j.ecolmodel.2007.09.006.

    Ceccarelli T, Bajocco S, Perini L,et al., 2014. Urbanization and land take of high quality agricultural soils—exploring long-term land use changes and land capability in northern Italy. International Journal of Environmental Research, 8(1): 181-192.

    Chaplot V, Giboire G, Marchand P,et al., 2005. Dynamic modelling for linear erosion initiation and development under climate and land-use changes in northern Laos. Catena, 63(2-3): 318-328.DOI: 10.1016/j.catena.2005.06.008.

    Chen W, Xiao D, Li X, 2002. Classification, application, and creation of landscape indices. The journal of applied ecology, 13(1):121-125.

    Chen Y, Zhang RZ, Chen B,et al., 2010. Construction and application on the niche theory of land use. Arid Land Geography, 33(5):791-801. DOI: 10.13826/j.cnki.cn65-1103/x.2010.05.022.

    Dan LW, Seifert SN, 2011. Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria. Ecological Applications: A Publication of the Ecological Society of America, 21(2): 335-342.

    Diggelen RV, Sijtsma FJ, Strijker D,et al., 2005. Relating land-use intensity and biodiversity at the regional scale. Basic & Applied Ecology, 6(2): 145-159. DOI: 10.1016/j.baae.2005.01.004.

    Fan Q, Ding S, 2016. Landscape pattern changes at a county scale: a case study in Fengqiu, Henan Province, China from 1990 to 2013.Catena, 137: 152-160. DOI: 10.1016/j.catena.2015.09.012.

    Feng Y, Luo G, Zhou D,et al., 2010. Effects of land use change on landscape pattern of a typical arid watershed in the recent 50 years:A case study on Manas River Watershed in Xinjiang. Acta Ecologica Sinica, 30(16): 4295-4305.

    Freitas MWDD, Santos JRD, Alves DS, 2013. Land-use and land-cover change processes in the Upper Uruguay Basin: linking environmental and socioeconomic variables. Landscape Ecology, 28(2):311-327. DOI: 10.1007/s10980-012-9838-9.

    Fu B, Liang D, Lu N, 2011. Landscape ecology: Coupling of pattern,process, and scale. Chinese Geographical Science, 21: 385-391.DOI: 10.1007/s11769-011-0480-2.

    Ganzhou Statistic Bureau, 2015. Ganzhou District Statistical Yearbook.

    Hao SL, Cao LH, Li BC, 2010. Application of niche theory to land utilization: A case study on Shanghuang in Guyuan. Soil and Water Conservation in China, 3(3): 35-37.

    Hu M, Zhang J, Liang E,et al., 2016. Analysis of land utilization ecological niche change in Manas River Watershed over the past 37 years. Journal of Shihezi University, 34(5): 571-575.

    Krummel JR, Gardner RH, Sugihara G,et al., 1987. Landscape patterns in a disturbed environment. Oikos, 48(3): 321-324.

    Li BL, 2000. Fractal geometry applications in description and analysis of patch patterns and patch dynamics. Ecological Modelling,132(1-2): 33-50. DOI: 10.1016/S0304-3800(00)00303-3.

    Li H, 1989. Spatial Temporal Pattern Analysis of Managed Forest Landscapes. Ph.D. Thesis. Oregon: The Oregon State University.

    Li X, Bu R, Chang Y,et al., 2004. The response of landscape metrics against pattern scenarios. Acta Ecologica Sinica, 24: 123-134.

    Cardoso MB, Ladio AH, Lozada M, 2017. Niche breadth and redundancy: Useful indices to analyse fuelwood use in rural communities.Journal of Arid Environments, 145: 52-59. DOI: 10.1016/j.jaridenv.2017.05.007.

    Nascimento APBD, Ferreira ML, Molina SMG, 2010. Ecological niche theory: non-traditional urban and rural human populations. Journal of Human Ecology, 32(3): 175-182.

    Pennock DJ, Anderson DW, Jong ED, 1994. Landscape-scale changes in indicators of soil quality due to cultivation in Saskatchewan,Canada. Geoderma, 64(1-2): 1-19. DOI: 10.1016/0016-7061(94)90086-8.

    Qi X, Mao J, Cheng Y, 2007. Study on optimizing strategies of city spaces basing on niche theory: A case in Zengcheng City. Planners,23(2): 89-91.

    Qi Y, Wu JG, Li JL,et al., 2013. Landscape dynamics of medium- and small-sized cities in eastern and western China: a comparative study of pattern and driving forces. Acta Ecologica Sinica, 33(1):275-285.

    Qin JC, Gao M, 2004. Economic niche of land use and sustainable development in mountains: A case study of Chongqing. Journal of Mountain Research, 21(6): 702-706.

    Robinson H, Weckworth B, 2016. Landscape ecology: linking landscape metrics to ecological processes. Snow Leopards. Elsevier Inc. pp. 395-405.

    Silvertown J, 2004. The ghost of competition past in the phylogeny of island endemic plants. Journal of Ecology, 92(1): 168-173. DOI:10.1111/j.1365-2745.2004.00853.

    Tong HL, Shi PJ, Zhu GF,et al., 2015. The characteristics of oasis urban expansion and drive mechanism analysis: a case study on Ganzhou District in Hexi Corridor, China. Sciences in Cold and Arid Regions, 7(3): 282-292. DOI: 10.3724/SP.J.1226.2015.00282.

    Tong H, Shi P, Bao S,et al., 2018. Optimization of urban land development spatial allocation based on ecology-economy comparative advantage perspective. Journal of Urban Planning and Development, 144(2): 05018006, 1-14. DOI: 10.1061/(ASCE)UP.1943-5444.0000444.

    Verburg PH, Koning GHJD, Kok K,et al., 1999. A spatial explicit allocation procedure for modelling the pattern of land use change based upon actual land use. Ecological Modelling, 116(1): 45-61.DOI: 10.1016/S0304-3800(98)00156-2.

    Wu JG, 2000. Landscape Ecology: Pattern, Process, Scale and Hierarchy. Beijing: Higher Education Press.

    Wu J, Shen W, Sun W,et al., 2002. Empirical patterns of the effects of changing scale on landscape metrics. Landscape Ecology, 17(8):761-782. DOI: 10.1023/A:1022995922992.

    Wu Q, Li Y, 2014. The representation of landscape pattern to land economic niche change: A case study of Huizhou City in Guangdong Province. Scientia Geographica Sinica, 34(6): 705-710. DOI:10.13249/j.cnki.sgs.2014.06.705.

    Xiao DN, Bu RC, Li XZ, 1997. Spatial ecology and landscape heterogeneity. Acta Ecologica Sinica, 17(5): 453-461.

    Yu HY, Zhang F, Wang J,et al., 2015. Application of land economic ecological niche in landscape pattern analysis at county level: A case study of Jinghe County in Xinjiang, China. Chinese Journal of Applied Ecology, 26(12): 3849-3857.

    Zhang GM, Xie SC, 1997. Development of niche concept and its perspectives: A review. Chinese Journal of Ecology, 16(6): 46-51.

    Zhang QJ, Fu BJ, Chen LD, 2003. Several problems about landscape pattern change research. Scientia Geographica Sinica, 23(3):264-270.

    Zhang X, Ge XD, Pu LJ,et al., 2002. A study on the economic niche of land use theory and mechanisms of cropland conservation.Journal of Natural Resources, 17(6): 677-683. DOI:10.11849/zrzyxb.2002.06.004.

    Zhu CQ, 1997. The niche ecostate-ecorole theory and expansion hypothesis. Acta Ecologica Sinica, 17(3): 324-332.

    久久久久精品人妻al黑| 在线观看免费视频日本深夜| 国产精品一区二区在线不卡| 最新的欧美精品一区二区| 成人免费观看视频高清| 精品一区二区三卡| 国产亚洲精品第一综合不卡| 美女视频免费永久观看网站| 欧美日韩av久久| 欧美日韩av久久| 亚洲一区二区三区不卡视频| 大码成人一级视频| 在线观看日韩欧美| 成人黄色视频免费在线看| 日韩欧美免费精品| 国产av精品麻豆| 亚洲中文av在线| 久久精品aⅴ一区二区三区四区| 成人国产一区最新在线观看| 国产高清视频在线播放一区| 精品一区二区三区视频在线观看免费 | 久99久视频精品免费| 最近最新免费中文字幕在线| 午夜亚洲福利在线播放| videosex国产| 日韩欧美三级三区| 亚洲精品久久成人aⅴ小说| 免费少妇av软件| 丁香欧美五月| 村上凉子中文字幕在线| 极品少妇高潮喷水抽搐| 丰满饥渴人妻一区二区三| 黄网站色视频无遮挡免费观看| 欧美成人午夜精品| 色老头精品视频在线观看| 国产在视频线精品| 国产成人精品久久二区二区免费| 国产成人免费无遮挡视频| 国产激情欧美一区二区| 窝窝影院91人妻| 久久影院123| 最新美女视频免费是黄的| 黄色a级毛片大全视频| 成人影院久久| 男女高潮啪啪啪动态图| 婷婷丁香在线五月| 中出人妻视频一区二区| 国产人伦9x9x在线观看| 黄片播放在线免费| 久久国产精品人妻蜜桃| 最近最新中文字幕大全免费视频| 51午夜福利影视在线观看| 亚洲av日韩在线播放| 亚洲国产精品一区二区三区在线| 亚洲aⅴ乱码一区二区在线播放 | 亚洲成人手机| 精品高清国产在线一区| 一本综合久久免费| 又黄又爽又免费观看的视频| 一边摸一边做爽爽视频免费| 动漫黄色视频在线观看| 成人av一区二区三区在线看| www.自偷自拍.com| 国产野战对白在线观看| 日韩 欧美 亚洲 中文字幕| av视频免费观看在线观看| 日韩欧美三级三区| 国产一区二区激情短视频| 国产精品一区二区在线观看99| 少妇裸体淫交视频免费看高清 | 99热网站在线观看| 久久久久久亚洲精品国产蜜桃av| 欧美另类亚洲清纯唯美| 亚洲国产中文字幕在线视频| 中文字幕色久视频| 欧美另类亚洲清纯唯美| 国产片内射在线| 黄色怎么调成土黄色| 精品高清国产在线一区| 免费在线观看视频国产中文字幕亚洲| 午夜福利影视在线免费观看| 国产精品乱码一区二三区的特点 | 久久久久久久精品吃奶| 中亚洲国语对白在线视频| 久久九九热精品免费| 欧美人与性动交α欧美精品济南到| 99久久精品国产亚洲精品| 中出人妻视频一区二区| 久久ye,这里只有精品| 91大片在线观看| 怎么达到女性高潮| 亚洲av片天天在线观看| 老司机深夜福利视频在线观看| 免费日韩欧美在线观看| 国产精品一区二区精品视频观看| 亚洲欧美日韩另类电影网站| 12—13女人毛片做爰片一| www.自偷自拍.com| 亚洲精品一二三| 免费高清在线观看日韩| 国产av又大| 午夜老司机福利片| 亚洲成a人片在线一区二区| 熟女少妇亚洲综合色aaa.| 成人18禁在线播放| 久久久久久久久免费视频了| 黑人欧美特级aaaaaa片| 国产aⅴ精品一区二区三区波| 黄片播放在线免费| 91精品国产国语对白视频| 久久久久国产精品人妻aⅴ院 | 国内毛片毛片毛片毛片毛片| 一夜夜www| 亚洲五月婷婷丁香| av福利片在线| 大片电影免费在线观看免费| 成年女人毛片免费观看观看9 | 黄色成人免费大全| 黄网站色视频无遮挡免费观看| 亚洲av成人不卡在线观看播放网| av线在线观看网站| 精品国产一区二区三区四区第35| 国产免费av片在线观看野外av| 亚洲精品美女久久av网站| videosex国产| 久久人妻福利社区极品人妻图片| 夜夜爽天天搞| 欧美乱妇无乱码| 操美女的视频在线观看| 国产片内射在线| 99精品在免费线老司机午夜| 成年版毛片免费区| 999精品在线视频| 精品亚洲成a人片在线观看| 91麻豆精品激情在线观看国产 | 在线免费观看的www视频| 久久久水蜜桃国产精品网| 欧美激情高清一区二区三区| a在线观看视频网站| 中文字幕人妻熟女乱码| 大码成人一级视频| 亚洲在线自拍视频| 亚洲av片天天在线观看| 精品少妇一区二区三区视频日本电影| 动漫黄色视频在线观看| 国产一区有黄有色的免费视频| 人人妻,人人澡人人爽秒播| 丝袜在线中文字幕| 丝瓜视频免费看黄片| 精品人妻在线不人妻| 国产亚洲欧美在线一区二区| 91av网站免费观看| 国产精品偷伦视频观看了| 久久精品国产99精品国产亚洲性色 | 欧美精品人与动牲交sv欧美| 久久九九热精品免费| 亚洲精品乱久久久久久| 欧美精品av麻豆av| 男女高潮啪啪啪动态图| 亚洲精品美女久久久久99蜜臀| 一级片免费观看大全| 欧美色视频一区免费| 人人妻人人添人人爽欧美一区卜| 一二三四社区在线视频社区8| 深夜精品福利| 国产伦人伦偷精品视频| a在线观看视频网站| 人人妻人人澡人人爽人人夜夜| 高清黄色对白视频在线免费看| 久99久视频精品免费| 国产欧美日韩综合在线一区二区| tocl精华| 亚洲精品自拍成人| 欧美成狂野欧美在线观看| 久久草成人影院| 欧美最黄视频在线播放免费 | 又大又爽又粗| 国产熟女午夜一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一区二区免费欧美| 91九色精品人成在线观看| 手机成人av网站| 亚洲精品中文字幕在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品二区激情视频| 天天操日日干夜夜撸| 亚洲人成电影免费在线| 女人被狂操c到高潮| 窝窝影院91人妻| 老司机在亚洲福利影院| 亚洲 国产 在线| 国产91精品成人一区二区三区| 久久精品国产亚洲av高清一级| 欧美人与性动交α欧美精品济南到| 国产精品欧美亚洲77777| 一级毛片精品| 亚洲精品成人av观看孕妇| 亚洲少妇的诱惑av| 国产精品二区激情视频| 最近最新中文字幕大全免费视频| 欧美成狂野欧美在线观看| 国产激情久久老熟女| 成人永久免费在线观看视频| 精品国产亚洲在线| 啦啦啦免费观看视频1| 久久精品人人爽人人爽视色| 十八禁人妻一区二区| 狂野欧美激情性xxxx| 男男h啪啪无遮挡| 中文亚洲av片在线观看爽 | 国产精品美女特级片免费视频播放器 | 狠狠狠狠99中文字幕| 成人手机av| av网站在线播放免费| 在线国产一区二区在线| 麻豆乱淫一区二区| 免费看a级黄色片| 亚洲熟妇中文字幕五十中出 | 亚洲九九香蕉| 亚洲色图 男人天堂 中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 精品卡一卡二卡四卡免费| 婷婷成人精品国产| 80岁老熟妇乱子伦牲交| 国产精品av久久久久免费| 美女 人体艺术 gogo| 久久影院123| 一级黄色大片毛片| 亚洲精品久久成人aⅴ小说| 一夜夜www| 国产精品亚洲av一区麻豆| 精品国产国语对白av| 国产精品影院久久| 亚洲精品国产精品久久久不卡| 91国产中文字幕| 丰满人妻熟妇乱又伦精品不卡| 下体分泌物呈黄色| 国产真人三级小视频在线观看| 中文字幕人妻丝袜制服| 精品国产一区二区久久| 看免费av毛片| 精品国产亚洲在线| 日韩欧美一区二区三区在线观看 | 最新在线观看一区二区三区| av国产精品久久久久影院| 久久久久久免费高清国产稀缺| 精品人妻在线不人妻| 中文字幕另类日韩欧美亚洲嫩草| 男人操女人黄网站| √禁漫天堂资源中文www| 国产免费男女视频| av中文乱码字幕在线| 美女高潮到喷水免费观看| 免费在线观看黄色视频的| 50天的宝宝边吃奶边哭怎么回事| 建设人人有责人人尽责人人享有的| 久久精品国产亚洲av香蕉五月 | 色播在线永久视频| 精品电影一区二区在线| 亚洲av成人av| 欧美成狂野欧美在线观看| 久久热在线av| 亚洲精品在线美女| 午夜免费观看网址| 精品福利永久在线观看| 欧洲精品卡2卡3卡4卡5卡区| 色老头精品视频在线观看| av欧美777| 亚洲精品中文字幕一二三四区| 久久久久久亚洲精品国产蜜桃av| 一级毛片女人18水好多| 咕卡用的链子| 亚洲国产精品sss在线观看 | 国产xxxxx性猛交| 亚洲精品中文字幕一二三四区| 亚洲成国产人片在线观看| 午夜福利影视在线免费观看| 国产亚洲精品第一综合不卡| 日韩视频一区二区在线观看| 两人在一起打扑克的视频| 成人永久免费在线观看视频| 午夜福利视频在线观看免费| 亚洲熟妇中文字幕五十中出 | 99热国产这里只有精品6| 99精品在免费线老司机午夜| 亚洲专区中文字幕在线| 日本黄色日本黄色录像| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲 欧美一区二区三区| 日韩免费av在线播放| 欧美老熟妇乱子伦牲交| a在线观看视频网站| 两个人看的免费小视频| 日本黄色视频三级网站网址 | 久久人妻av系列| 婷婷成人精品国产| 黄色视频不卡| 99久久国产精品久久久| 久久 成人 亚洲| 亚洲av成人不卡在线观看播放网| 日韩熟女老妇一区二区性免费视频| 精品一区二区三区四区五区乱码| 在线观看舔阴道视频| 国产一区有黄有色的免费视频| 久久久久久久精品吃奶| 精品久久蜜臀av无| 国产精品国产高清国产av | 首页视频小说图片口味搜索| 两性午夜刺激爽爽歪歪视频在线观看 | a在线观看视频网站| 高潮久久久久久久久久久不卡| 午夜精品国产一区二区电影| 国产成人精品在线电影| 女人精品久久久久毛片| 欧美 亚洲 国产 日韩一| av在线播放免费不卡| 国产精品 国内视频| 人人妻人人添人人爽欧美一区卜| 男人舔女人的私密视频| 亚洲中文av在线| 精品欧美一区二区三区在线| 色老头精品视频在线观看| 热99久久久久精品小说推荐| 成年人黄色毛片网站| av国产精品久久久久影院| 久久人人97超碰香蕉20202| 极品教师在线免费播放| 大码成人一级视频| 黑人巨大精品欧美一区二区mp4| 99re6热这里在线精品视频| av福利片在线| 在线观看舔阴道视频| 国产不卡一卡二| 欧美亚洲 丝袜 人妻 在线| 黑人欧美特级aaaaaa片| 欧美日韩亚洲高清精品| 满18在线观看网站| 久久国产精品人妻蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品成人在线| 美女福利国产在线| 精品国产一区二区三区四区第35| 免费看十八禁软件| 50天的宝宝边吃奶边哭怎么回事| 精品电影一区二区在线| 精品久久久久久久毛片微露脸| 国产亚洲欧美98| 亚洲欧美精品综合一区二区三区| 午夜福利,免费看| 国产精品美女特级片免费视频播放器 | 久久久久视频综合| 97人妻天天添夜夜摸| 久久中文看片网| 亚洲欧美激情综合另类| www日本在线高清视频| 日韩欧美一区视频在线观看| 中文欧美无线码| 99精品在免费线老司机午夜| 91成人精品电影| 满18在线观看网站| 亚洲欧美日韩高清在线视频| www.自偷自拍.com| 人成视频在线观看免费观看| 久久影院123| 18禁观看日本| 欧美在线一区亚洲| 黄色a级毛片大全视频| 国产在线一区二区三区精| 69av精品久久久久久| 久久 成人 亚洲| 久久狼人影院| 亚洲人成电影观看| 国产一区二区三区视频了| 美女高潮到喷水免费观看| 成熟少妇高潮喷水视频| 日本黄色视频三级网站网址 | 亚洲欧美日韩高清在线视频| 日韩欧美国产一区二区入口| 69av精品久久久久久| 成人影院久久| 人人妻,人人澡人人爽秒播| 一级,二级,三级黄色视频| 日日爽夜夜爽网站| 国产精品久久久av美女十八| 大香蕉久久成人网| 美国免费a级毛片| 欧美黄色片欧美黄色片| 老司机午夜十八禁免费视频| 变态另类成人亚洲欧美熟女 | 色婷婷久久久亚洲欧美| 一二三四社区在线视频社区8| 欧美丝袜亚洲另类 | 淫妇啪啪啪对白视频| 国产男靠女视频免费网站| 91国产中文字幕| 国产精品久久久久成人av| 深夜精品福利| 777米奇影视久久| 身体一侧抽搐| 男女之事视频高清在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲成人免费av在线播放| 1024视频免费在线观看| 日日夜夜操网爽| 国产主播在线观看一区二区| 日本一区二区免费在线视频| 一边摸一边抽搐一进一出视频| 少妇被粗大的猛进出69影院| 欧美色视频一区免费| 成年人午夜在线观看视频| 黄片小视频在线播放| 首页视频小说图片口味搜索| 精品国产超薄肉色丝袜足j| 久久人妻熟女aⅴ| 人人澡人人妻人| 日韩 欧美 亚洲 中文字幕| 欧美成狂野欧美在线观看| 久久精品亚洲熟妇少妇任你| 日本黄色日本黄色录像| 午夜久久久在线观看| 极品少妇高潮喷水抽搐| 十八禁人妻一区二区| 咕卡用的链子| 久久久久精品人妻al黑| 老司机福利观看| www.自偷自拍.com| 十八禁人妻一区二区| 国产亚洲精品一区二区www | 国产精品 欧美亚洲| 新久久久久国产一级毛片| 1024视频免费在线观看| 女人爽到高潮嗷嗷叫在线视频| 人妻一区二区av| 亚洲综合色网址| 亚洲全国av大片| 亚洲国产欧美一区二区综合| 亚洲欧美一区二区三区久久| 美女高潮喷水抽搐中文字幕| 麻豆av在线久日| 宅男免费午夜| 夜夜夜夜夜久久久久| 成人特级黄色片久久久久久久| 国产亚洲精品久久久久5区| 我的亚洲天堂| 国产淫语在线视频| 视频区图区小说| 亚洲欧美精品综合一区二区三区| 露出奶头的视频| 男女高潮啪啪啪动态图| 国产在线一区二区三区精| 国产成人影院久久av| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成av片中文字幕在线观看| 老司机福利观看| 99久久人妻综合| 丝袜人妻中文字幕| 香蕉久久夜色| 国产午夜精品久久久久久| 久久影院123| 国产又色又爽无遮挡免费看| 国产欧美日韩综合在线一区二区| 侵犯人妻中文字幕一二三四区| 香蕉丝袜av| 亚洲精品av麻豆狂野| 亚洲成人国产一区在线观看| 免费不卡黄色视频| 黄色 视频免费看| 中亚洲国语对白在线视频| 色婷婷久久久亚洲欧美| 岛国毛片在线播放| 国产成人啪精品午夜网站| 精品久久久久久,| 免费观看a级毛片全部| 一级黄色大片毛片| 最近最新中文字幕大全免费视频| xxxhd国产人妻xxx| 他把我摸到了高潮在线观看| 免费在线观看亚洲国产| 老汉色∧v一级毛片| 免费看十八禁软件| av不卡在线播放| 女性生殖器流出的白浆| 午夜成年电影在线免费观看| 女人高潮潮喷娇喘18禁视频| 国产熟女午夜一区二区三区| 中文亚洲av片在线观看爽 | 国产成+人综合+亚洲专区| 亚洲中文av在线| 亚洲在线自拍视频| 另类亚洲欧美激情| 一级片'在线观看视频| 色在线成人网| 国产一区二区激情短视频| av一本久久久久| 日本欧美视频一区| 精品无人区乱码1区二区| 国产无遮挡羞羞视频在线观看| 亚洲av第一区精品v没综合| 亚洲在线自拍视频| 欧美老熟妇乱子伦牲交| av一本久久久久| 搡老岳熟女国产| 极品人妻少妇av视频| 久久久久久人人人人人| 黄片大片在线免费观看| 亚洲欧美精品综合一区二区三区| 精品国产乱码久久久久久男人| 少妇猛男粗大的猛烈进出视频| 天堂√8在线中文| 久久热在线av| 国产欧美日韩一区二区三| 国产精品欧美亚洲77777| 美女午夜性视频免费| 女人被狂操c到高潮| 91成人精品电影| 99热网站在线观看| av线在线观看网站| 免费观看a级毛片全部| 午夜两性在线视频| 国产成人影院久久av| 成人特级黄色片久久久久久久| 91国产中文字幕| 成人av一区二区三区在线看| 一进一出抽搐动态| 亚洲三区欧美一区| 日韩大码丰满熟妇| 国产精品99久久99久久久不卡| 交换朋友夫妻互换小说| 国产精品乱码一区二三区的特点 | a级毛片黄视频| 国产有黄有色有爽视频| 热re99久久国产66热| videosex国产| 真人做人爱边吃奶动态| 成人免费观看视频高清| 日日摸夜夜添夜夜添小说| 国产欧美亚洲国产| 免费在线观看日本一区| 国产欧美亚洲国产| 亚洲精品av麻豆狂野| 岛国毛片在线播放| 国精品久久久久久国模美| 男女床上黄色一级片免费看| 我的亚洲天堂| 在线观看一区二区三区激情| 夜夜夜夜夜久久久久| 国产男靠女视频免费网站| 操美女的视频在线观看| 女人高潮潮喷娇喘18禁视频| 极品人妻少妇av视频| 极品少妇高潮喷水抽搐| e午夜精品久久久久久久| 亚洲精品在线观看二区| 性少妇av在线| 男女午夜视频在线观看| x7x7x7水蜜桃| 女人高潮潮喷娇喘18禁视频| 免费不卡黄色视频| 国产在线观看jvid| 熟女少妇亚洲综合色aaa.| 精品视频人人做人人爽| 高清av免费在线| 久久精品国产亚洲av香蕉五月 | 欧美av亚洲av综合av国产av| 可以免费在线观看a视频的电影网站| av天堂在线播放| 国产区一区二久久| 波多野结衣av一区二区av| 精品欧美一区二区三区在线| 国产成人系列免费观看| 亚洲第一欧美日韩一区二区三区| 欧美av亚洲av综合av国产av| 欧美一级毛片孕妇| 国产精品一区二区在线不卡| 国产精品影院久久| 亚洲成a人片在线一区二区| 国产在线观看jvid| 亚洲avbb在线观看| 丰满的人妻完整版| 精品一区二区三区av网在线观看| 成年人免费黄色播放视频| 午夜日韩欧美国产| 精品国产乱码久久久久久男人| 国产亚洲精品一区二区www | 亚洲av片天天在线观看| 免费观看人在逋| 老司机亚洲免费影院| a级毛片在线看网站| 亚洲三区欧美一区| 如日韩欧美国产精品一区二区三区| 国产99白浆流出| 操出白浆在线播放| 涩涩av久久男人的天堂| 久久精品国产99精品国产亚洲性色 | 成年人免费黄色播放视频| 精品国产乱子伦一区二区三区| 王馨瑶露胸无遮挡在线观看| 丰满的人妻完整版| 99久久综合精品五月天人人| 久久久久久人人人人人| 夜夜爽天天搞| 亚洲专区字幕在线| 一个人免费在线观看的高清视频| 亚洲国产欧美网| 久久中文字幕一级| 91在线观看av| 人人妻人人澡人人爽人人夜夜| 免费观看精品视频网站| 日韩欧美三级三区| 亚洲欧美日韩高清在线视频| 天堂√8在线中文| 欧美最黄视频在线播放免费 | 久久久精品区二区三区| 久久九九热精品免费| 欧美亚洲日本最大视频资源|