A landscape connectivity index for assessing desertification

  • 格式:pdf
  • 大小:424.77 KB
  • 文档页数:13

Abstract As a global and regional environ-mental problem,desertification assessment is an instrumental component in developing global/regional actions plans aimed at preventing and/or eradicating desertification ing a landscape assessment approach allows for rela-tively quick assessments of desertification that can then be used in developing practicable application plans at the regional level in desert-ification prevention planning and decision-mak-ing.This study was conducted to determine whether a cost-distance connectivity index could both reveal evidence of,and act as an indicator for,desertification.Cost-distance,a simple GIS-calculated connectivity measure,was applied to a 1997land use map to indicate desertificationin Minqin County,China.The results showed that connectivity based on cost-distance follows basic landscape ecological principles including species–area relations and edge effects,and indicates desertification.Although grain size had a significant effect on the cost-distance index,especially at patch boundaries,log cost-distance closely corresponded with degree of desertification within grain size.Changing the extent of analysis had no significant effect on the cost-distance index and its relevance to degree of desertification.The Minqin County landscape had a high level of connectivity,although the area’s grasslands,oasis irrigated cultivated lands,alkali-saline lands and forest-lands played important roles in resisting desert-ification.Three areas require restoration of native vegetation or afforestion to cut the connectivity of desertified patches.The applica-tion of connectivity based on cost-distance provides a straightforward,easily visualized description of desertification.In addition,land use data is readily available in China,allowing for relatively easy and quick assessments of regional level desertification for planning and decision-making.Keywords Landscape connectivity ÆCost-distance ÆLand use types ÆDesertification assessment ÆChinaD.Sun ÆB.LiMinistry of Education Laboratory for Plant–Soil Interaction Processes,College of Natural Resources and Envir.Science,China Agricultural University,Beijing 100094,ChinaH.LiBeijing Academy of Agriculture and Forestry,Beijing 100089,ChinaR.Dawson ÆR.Wei (&)Laboratory for Earth Surface Processes,College of Environmental Science,Peking University,Beijing 100871,Chinae-mail:rwdawson@Landscape Ecol (2007)22:531–543DOI 10.1007/s10980-006-9046-6RESEARCH ARTICLEA landscape connectivity index for assessing desertification:a case study of Minqin County,ChinaDanfeng Sun ÆRichard Dawson ÆHong Li ÆRong Wei ÆBaoguo LiReceived:7December 2005/Accepted:4September 2006/Published online:19October 2006ÓSpringer Science+Business Media B.V.2007IntroductionDesertification—land degradation in arid,semi-arid and dry sub-humid areas resulting from various factors,including climatic variations and human activities—is becoming one of the world’s more serious environmental problems(Glenn et al.1998).In China,and particularly in west China,the over use and misuse of land has contributed to a number of serious environmental problems,especially desertification.In1999Chi-na had2.67million km2of land(27.9%of total land area)classified as experiencing some degree of desertification,an increase of52,000km2 compared to1995(State Forestry Administration of China2002).Because west China is the focus of many national economic development directives, the ability to monitor and model patterns and processes of desertification is important not only for the rehabilitation of existing lands,but also to mitigate the processes before the lands become desertified.The Minqin Oasis in China’s Gansu Province has become one of China’s more severely deser-tified ecological regions over the last several decades and represents a serious source of sand-storms that increasingly plague northeast China (Man et al.2001;Gu2002;Wang et al.2003;Liu 2004).Desertification in Minqin County is prin-cipally driven by intensification of agricultural activities,which caused the deficiency of surface water resources in the county(Ma and Wei2003; Yang2003;Sun et al.2006).Although a campaign to combat desertification in Minqin County was initiated in2002,the lack of proper data made it difficult to quantify the extent of desertification and its rate of change(Li and Chen2001).The role spatial pattern plays in ecological processes(Pickett and Cadenasso1995)has made possible the development of landscape assess-ment approaches for regional environmental quality assessment and for monitoring land use and land cover types(O’Neill et al.1997;Lee et al.1999;Kepner et al.2000).Desertification is a global environmental problem and its assess-ment has made increasing use of landscape ecology ndscape metrics derived from land use and land cover maps have been used to quantify environmental change in arid and semi-arid regions(Kepner et al.2000;Jia et al.2004).Spatial heterogeneity has been shown to be one of the most reliable indicators of desertification(Schlesinger et al.1990).In this research,landscape metrics and spatial heteroge-neity indicators were combined to model desert-ification changes in Minqin County,China(Sun et al.2005).Landscape connectivity can be defined by metrics to reflect spatial heterogeneity because connectivity represents‘‘the degree to which the landscape facilitates or impedes the movement of organisms among source patches’’(Taylor et al. 1993;Tischendorf and Fahirg2000).Recently, however,there has been a shift in the form of approach with a move away from structural connectivity to one centered more on functional connectivity.The approach has gained some momentum as research results better clarify how different landscape elements affect dispersal, which can be both species and process specific (Chardon et al.2003;Adriaensen et al.2003; Marulli and Mallarach2005).Recent studies have shown that most landscape metrics,including connectivity,are scale-sensitive(Wu et al.2002; Wu2004).From a functional connectivity point of view, we can think of desertification as a resistance component within the landscape matrix such that land use types can hinder or enhance movements between patches with bare soil or sand.This connectivity,therefore,relates to not only the pattern of desertified patches and desertification movement,but also to the patterns of land use and land cover types having different levels of resistance.By measuring connectivity,it may be possible to reveal the interaction of desertifica-tion process with landscape pattern for assessing both the risk of desertification as well as its stage of development within a region.Connectivity may also help in rehabilitating or combating existing desertified areas by identifying landscapes requir-ing restoration.The major objective of this study was to evaluate whether connectivity can be used to assess desertification and serve as an indicator of desertification.We used Minqin County as a case study to:(1)calculate connectivity using the cost-distance index and to analyze its relationship todegree of desertification,and (2)explore the relationships between connectivity and landscape patch structure.We also addressed how scale or spatial resolution (grain size)affects the cost-distance measure as an indicator of desertifica-tion.Materials and methods Study areaMinqin County is located in Gansu Province,northwest China (ranging from 101°49¢to 104°12¢E and from 38°03¢to 39°28¢N)with an estimated area of 15,870km 2.The Baidan Jilin Desert,China’s third largest,is situated to the west and north while the Tengger Desert,the fourth largest,is located to the east.Because the study area is situated between two major deserts (Fig.1),oases represent important ecological corridors as well as key areas for national desert-ification monitoring and prevention assessment.The county has a predominant oasis landscape that slopes downward from ~1,460m in the southwest to ~1,295m in the northeast.The area has an arid continental climate with an average annual temperature of 7.8°C (average maximum 23.2°C in July and average minimum of –9.6°C in January).Mean annual precipitation in the county is ~110mm,half of which comes in July and August.Annual average evaporation is 2,664mm,24times annual mean precipitation.The only surface water source for irrigationcomes from the Shiyang River.However,the expansion of agricultural activity in the upper reaches of Shiyang River since the 1990s has led to the over-consumption of water resources to the extent that annual surface water flow (1.3·108m 3–1.5·108m 3in mid-1990s)now satisfies only 20%of the total water requirements for agriculture in the oasis (Minqin Water Con-servancy Bureau 1996).Economically,Minqin County is best characterized as an oasis farming area in an arid desert zone with agriculture accounting for 80%of gross domestic production.Of this,plantation agriculture and stock hus-bandry account for 77%and 20%of agricultural incomes,respectively.The Minqin oasis has become desertified over the last decade as the result of expanded regional economic growth initiatives and continued popu-lation pressure.For lack of proper data to reflect the extent of desertification and its rate of change since 1980s,Sun et al.(2005)produced a series of land desertification maps for 1988,1992and 1997based on TM series remote sensing classification of land surface spectral information,albedo and NDVI images.Analysis of the mapping results indicated that wind erosion was the dominant cause of desertification in 1997,affecting more than half of the study area (53.34%).Although moderate desertification was found to be the dominant class of desertification (43.64%of total area),extreme/severe desertification had ex-panded to become the second largest class (26.15%of total area in 1997)since the late 1980s.Areas of non-desertified oasislandscapeFig.1The Minqin study areawere becoming more fragmented and isolated in addition to being less stable from1988–1997. Cost-distance methodologyThe cost-distance method,based on a least-cost analysis from graph theory,is a standard function in GIS software(such as ARC/INFO,ARC-VIEW and IDRISI).Two GIS grid-layers,a source layer and a friction/resistance layer form the input of the model.The source layer indicates the target patches from which connectivity is to be calculated,while the resistance layer(cost-grid)is a map of the landscape matrix identifying land use types.Every cell on this cost-grid has a resistance value(cost)depending on its land use type and land cover nd use types that hinder ecological processes get values according to the strength of hindrance.Ideally, these values are based on empirical data,but in those cases where data are lacking it is necessary to rely on expert judgment(Chardon et al.2003). For target patches,ecological processes have no resistance.Thus,a value of1indicates that one step completely transfers between the sources cells(Chardon et al.2003;Adriaensen et al. 2003).Therefore,model output represents a cost-distance grid in which every grid cell has an exact value indicating its least(accumulated)cost path(based on single paths)from a defined source layer.This allows the cost-distance value to be interpreted as a connectivity measure that is weighted by the intervening matrix.For example, a cost value of n indicates the cost of moving through n cells with a resistance value of1,or through one cell with a resistance value of n (Adriaensen et al.2003).This cost-distance, which can be converted to an equivalent distance measure by multiplying the cost by cell size,can also be interpreted as the effective distance for the ecological translation of the calculated cost-values(Ricketts2001;Chardon et al.2003). Setting the source layer and the resistancelayer in the land use mapThe Minqin County land use map(1:10,000scale) for1997,collected from the Minqin County Land Management Bureau,was based on the national standardized land classification scheme(Minqin Land Management Bureau1993)to identify land use types in Minqin County(Table1).To run the cost-distance model,the target patches had to be first identifind use types consisting of sand or bare soil were considered desert while rock and/or gravel surfaces were considered Gobi.Table1Minqin county land types following China’s national standardized land classification schemeLabel number Land use type Definition13Irrigated cultivated land The cultivated land with irrigation device21Orchards The land planted fruit tree with canopy cover over50%31Forest The natural or artificial trees with closure percent over30% 32Shrub The shrub cover over40%33Sparse forest Forest with trees closure between10%and30%41Natural pasture Natural herbaceous plants for herd51City&town52Rural settlement53Stand-alone industry/mining site71River73Reservoir74Pond75Bulrush land76Beach81Waste grassland Scattered grasses with less than10%ground cover82Alkali-saline land Land with halophyte only83Swamp Often covered by water and growing hygrophyte84Sand land Surface covered by sand with scarcely any vegetation85Bare soil Surface covered by soil with scarcely any vegetation86Rock and/or gravel Surface covered by over70%rock and/or gravelTherefore,patches of these land use types com-posed the source layer.Because many rivers have beach sections,areas capable of providing sand for wind erosion,beach patches were also chosen as target patches for measuring connectivity.The second input,the cost-grid,contains the estimated cost(resistance value)to move through each grid cell in the study area.The resistance values were estimated for each land use type based on its canopy cover,which could be viewed as an attribute that affects sand transport and desertification process.The resistance values in Table2are conservative estimates of stopping a given percent of sand from moving for each land use type as adopted from the FAO(1979)and research results of Li(1988)and Zhao et al. (2005).Built-up areas,including cities and towns, rural settlements and stand-alone industry/mining sites,also strongly impede sand movement. Resistance values for these land use types were based on research results from Zhang and Wang (1999).Rivers,reservoirs,ponds and bulrush lands were assigned high resistance value because they greatly hinder sand encroachment and pro-vide water for vegetation growth.Four land use types were assigned resistance values of1(Fig.2).A grid cell size of30by30m was adopted as the base spatial resolution of the land use map to coincide with the1997Desertification Grade Map of Minqin County produced by Sun et al.(2005), the size of which was5,107by6,702pixels.To investigate the effects of changing grain size,the spatial resolution of the1997land use map was systematically increased by2n(n from1to10)on the base size(Wu et al.2002).These land use maps were then used to derive the corresponding source and resistance layers for the cost-distance calculation.To investigate the effects of changing extent,we systematically extracted a series of land use sub-maps at1,000by1,000,2,000by 2,000,3,000by3,000,4,000by4,000to5,000by 5,000pixels with Minqin County as the maps center while keeping the grain size at30by30m. In addition,the cost-distance,as the measure of desertification connectivity,with the smaller cost value was taken to indicate higher desertification connectivity.Data analysisInputting the corresponding source patches layers and the resistance layers produced a series ofTable2The set of resistance values(r)for land use types of Minqin CountyLand use label1321313233415152537173747576818283848586r50707085603080808060606060130560111Fig.2Land use type resistance values and desertification source patches(the number in parentheses is the land use type code)cost-distance maps that were transformed to natural logarithms to compress the range for display and further analysis.A cost-distance grid map with base grain size was used to analyze the relationship between desertification connectivity and landscape patch ing each desert-ification level in1997as a zone,a zonal mean value of cost-distance was calculated and plotted (scattergram)to reflect the relationship with degree of desertification.FRAGSTASTS(McGa-rigal et al.2002)was used to calculate land use patch metrics including AREA(patch area), PARA(patch perimeter-area ratio),SHAPE, FRAC(fractal dimension index)and CONTIG (contiguity index)to measure patches size,shape complexity and spatial connectedness.Bivariate correlation analysis was use to calculate the coefficients between patch connectivity and these patch metrics.The slope grid of connectivity was derived from the cost-distance grid and used to explore gradient change particularly at patch boundaries.Soil and vegetation are two important factors influencing desertification movement and resis-tance anic matter and CaCO3are components important in relating soil texture to susceptibility to wind erosion.The fraction of soils subject to erosion(SEF)in Minqin County was calculated according to Fryrear et al.(1994).A Normalized Difference Vegetation Index (NDVI)reflecting vegetation status was also calculated from August1997TM remote sensing data.The correlation coefficients of cost-distance, SEF and NDVI were then calculated to validate the rationality of cost-distance based on resis-tance values of land use types.In order to investigate the effects of changes in grain size on the cost-distance calculation,a series of simulated cost-distance maps with grain sizes ranging from30m by2n(n from1to10)were produced by adding all the pixels values at the 30-m level that were covered by a higher level pixel.The series of simulated cost-distance maps were also transformed to natural logarithm.Abso-lute value maps of the difference between the calculated cost-distance maps and the correspond-ing simulated maps were adopted to reflect the effects of changing grain size on the cost-distance calculation;the larger the difference value the more influence of grain size had on calculation. The mean absolute difference values in desertifi-cation classes and the whole county area were then calculated and plotted(scattergram)to reflect the effects of changing grain size.In addition,the class mean values of connectivity with corresponding grain size were calculated and plotted(scatter-gram)to compare the effects on their relationship with the degree of desertification.In order to investigate the effect of changing extent on the cost-distance calculation,a series of cost-distance maps with corresponding extent were extracted from the cost-distance map calcu-lated from the whole map extent.The absolute value maps of the difference between the cost-distance maps calculated from the extracted land use map and the corresponding extracted cost-distance maps were adopted to reflect the effects on the cost-distance calculations from changes in extent,with larger difference values indicating the increased influence of extent on the calculation. The mean absolute difference values in desertifi-cation class and the whole area were then calculated and plotted(scattergram)to reflect the effects of changing extent.Class mean values of connectivity with corresponding extent were also calculated and plotted(scattergram)to com-pare the effects on their relationship with degree of desertification.ResultsThe1997land use landscape(base grain size) had high overall connectivity(low cost-distance) (Fig.3),although oasis and grassland(including natural pasture and waste grassland)patches had low connectivity levels(high cost-distance). Grasslands accounted for24.11%of landscape area in1997,but accounted for60.63%of the cost-distance connectivity of the total landscape in resisting desertification(Table3).While alkali-saline lands and forestlands(including forest,shrub and sparse forest)played an important role in resisting desertification,their island like distribution within the desert land-scape(Figs.2,3)was not sufficient to stop the two deserts from merging(connecting)(Fig.3; points A,B&C).The log cost-distance calculated using base grain size had a staircase-like response curve with degree of desertification (Fig.4),indicating that more severe desertification had higher connectiv-ity (smaller cost-distance).Mean log cost-distance differences were significantly large between mod-erate desertification and severe desertification,and between non-desertification and mild desert-ification.There was only a small difference between mild and moderate desertification.These findings reflect the difficulty of using the cost-distance alone in classifying desertification as mild or moderate,and in correctly classifying changes from non-desertified to mild desertifica-tion.When grain size was below 27times basesize,log cost-distance had a similar increasing response for all degrees of desertification,but when grain size changed to 29times base size,the cost-distance declined for mild and non-desertifi-cation classes.The cost-distance of patches had a significant positive correlation with patch area,shape com-plexity,fractal dimension and spatial connected-ness and a significant negative correlation with the ratio of patch perimeter to area (Table 4).The AREA,SHAPE,FRAC and CONTIG indices had significant positive correlations with each other and a significant negative correlation with PARA.The change gradient of cost-distance took place mainly near the interfacebetweenFig.3The 1997cost-distance grid after natural log transformation with base grain size.Thedivision into five contours is chosen for visibility reasons,but can beadjusted according to user requirements since every grid cell has its own cost-distanceTable 3Land use type area percents andlogarithmic cost-distance percent for the whole landscape matrixLand use label 41,81138231,32,33Sum Area (%)24.11 6.17 6.69 2.7739.74Cost-distance (%)60.6315.9913.676.796.99source patches and resistance patches,decreasing rapidly beyond the patch boundary (Fig.5).The soils in Minqin County are very suscepti-ble to wind erosion with a high erodible fraction.The relatively low values were distributed in the oasis areas in the middle of Minqin County (Fig.6).Vegetative coverage was also distributedaround the center of the county with the oases,although most other parts had negative NDVI values suggesting no vegetation (Fig.7).SEF and NDVI had similar spatial distributions on the basis of the cost-distance map.Differences that could be observed occurred in the northwestern part where the patches were the result of naturalTable 4Bivariate correlation coefficients of patch logarithmic cost-distances and the selected patches metricsCost-distanceAREA PARA SHAPE FRAC CONTIG Cost-distance 10.32**–0.15**0.36**0.14**0.15**AREA 0.32**1–0.14**0.35**0.15**0.14**PARA –0.15**–0.14**1–0.30**–0.13**–0.99**SHAPE 0.36**0.35**–0.30**10.80**0.31**FRAC 0.14**0.15**–0.13**0.80**10.15**CONTIG0.15**0.14**–0.99**0.31**0.15**1**Correlation is significant at the 0.01level (2-tailed),and patches number (n )is1,092Fig.5(a )the gradient contour of the log cost-distance grid with land use boundary,(b )the zooming gradient contour with land use boundarypastureland use,although this contradiction may be the result of error in identifying land use type.Table 5shows that cost-distance had a significantnegative correlation with SEF,and a significant positive correlation with NDVI indicating that the higher the erodible fraction of soils with no vegetation results in lower cost-distance and higher connectivity.As such,the relationship between cost-distance,SEF and NDVI reflects the rationality of cost-distance based on resis-tance values of land use type.Effects of changing grain size and extent Mean absolute difference of log cost-distance had an exponential relation to changes in grain size (Fig.8),which was similar for all degreesofFig.6Soil erodible fraction map of MinqinCountyFig.7Normalized Difference Vegetation Index map in August,1997.TM image of Minqin CountyTable 5Bivariate correlation coefficients between logarithmic cost-distances,soil erodible fraction and normalization difference vegetation index in Minqin CountyCost-distanceSEF NDVI Cost-distance 1–0.501**0.299**SEF –0.501**1–0.348**NDVI0.299**–0.348**1**Correlation is significant at the 0.01level (2-tailed),and pixels number (n )is 17,162,990desertification(Table6).In other words,as grain size increased above24times base size,there was a rapid increase in cost-distance differences (Fig.9).Mean absolute difference log cost-distance appeared to decrease as analysis extent expanded (Fig.10),but these differences were not signifi-cant.Log cost-distance had similar staircase-like response curves to degree of desertification for different analysis extents(Fig.11),becoming more linear as analysis extent decreased. Discussion and conclusionsOurfindings support our hypothesis that cost-distance connectivity for land use maps can serve as a useful indicator of desertification because it incorporates the connecting resistance for differ-ent land use types,hence,reflecting landscape function patterning.As a global and regional environmental problem,desertification assess-ment,including its status,change and trend,will be instrumental in developing global and regional action plans to prevent and eradication desertifi-cation.Overall,the Minqin landscape had a high level of connectivity(low cost-distance)with only grasslands,oasis irrigated cultivated lands,alkali-saline lands and forests showing signs of resis-tance.A landscape assessment approach based on land use data has several advantages for assessing desertification.First,land use is a familiar attri-bute in landscape descriptions and relatively easy to interpret either from ground surveys or from remote-sensing imagery.Second,historic land use data is relatively easy to collect for purposes of assessing change.In China,a large-scale national land survey took place from1984to1996,which was based on a standardized land classification scheme consisting of eight1-digit categories and forty-six2-digit categories.Since then,all land use change data are adjusted each year on31Octo-ber.This means that the updating of land use maps is routine in a manner similar to that of the national population census(Ma2000).The land use data now available allows for easy,quick assessments and applications at the regional level for desertification prevention planning and deci-sion-makings.Third,land use reflects the inter-action of biophysical factors and socio-economic conditions at specific points in time and,there-fore,represents an invaluable‘‘summary’’of the physical,socio-economic environment of the desertification region.It also can be validated from the relationship between cost-distance based on resistance values of land use types,SEF and NDVI.Finally,given the falling costs and increas-ing accessibility of remotely sensed imagery,it is possible that desertification monitoring based on landscape assessment with land use data could be undertaken cheaply and efficiently.Remote sens-ing may be used to identify target patches and to set resistance values in ways similar to how vegetation cover changes have been successfullyTable6The function of mean absolute difference log-distance with grain sizeNon-desertification Milddesertification ModeratedesertificationSeveredesertificationWholeCountyEquation y=0.0406e0.5285n y=0.0406e0.5185n y=0.0313e0.5347n y=0.0315e0.5206n y=0.0338e0.5284n R20.9810.9790.9870.9830.984y is the mean absolute different log cost-distance,n is the exponent of2n times of base grain size。