| 
 
      The combination of dynamic 
      simulation with GIS for evaluation and prediction of ecological benefit of 
      shelter forest  
 Zhang Zhiyong 
      National Lab of REIS, Institute of Geography, Chinese Academy of 
      Sciences, Beijing, China
 
 Zhou Xintie
 National Remote 
      Sensing Centre of China
 
 Hsu Chochun Li Qinming
 Dept. of 
      Computer Science, Peking University, Beijing, 
      China
 
 Abstract
 This paper puts forward a 
      new method combining the mathematical models in social and economic 
      development, in which system dynamics model is used as a min body, with 
      spatial analysis in ecological environment based on GIS. This method has 
      been applied in Pingquan Country.
 
 Introduction
 It is 
      known that the regional socio-economic development greatly affects the 
      ecological environments, and the ecological environments constrain the 
      life of human being and economical development. Therefore, it is important 
      to make a harmonic coordination between the regional socio-economic 
      development and the changes of ecological environments. However, the 
      current researches on socio-economic development and prediction are often 
      divorced from those on ecological environmental changes. It is to find a 
      research method which models both of the development of society and 
      economics and the improvement. Through the project of remote sensing 
      investigation in Pingquan County, a series of thematic maps, such as land 
      use map, forest distribution map, and forest dynamic map, etc. have been 
      completed. Meanwhile, a new GIS software system-Geo-Union-has been 
      developed. On the basis of the research work mentioned above, the 
      following research work is addressed:
 
        General Situations of the study 
      areaBy inspecting the relations between forest coverage and soil and 
        water loss, the regional comprehensive development and prediction models 
        have been established.
 
Based on the spatial analysis function of GIS, the predictive data 
        of land use types, land construction and forest coverage, etc. have been 
        spatially planned and an series of thematic maps in the future time 
        sequence have been produced. 
 
According to the changes of land use and plant coverage, the 
        prediction of soil and water loss and ecological environment change have 
        been analysed.
 
 Pingquan County, located in the north-east of Hebei Province, 
      is a mid-low mountain area with a rough ratio of seven-tenth Mountain, 
      one-tenth water and two-tenth farmland. The unreasonable human activities 
      in the past, have resulted in a lot of social, economic and environmental 
      problems:
 
        Among the problems, the most serious one is the ecological 
      environmental problem, that is, the soil and water loss. It is the main 
      reason of the proposing construction of shelter forest system in the 
      region.Rapid expansion of population has intensified the seek for food and 
        energy, which have caused farming on unsuitable land exceovercutting of 
        forest, and over-grazing of grass land. Consequently the intensive soil 
        and water loss and deterioration of ecological environment have also 
        been caused. With the construction of shelter forest system in recent 
        years, the improvement of ecological environment is promising.
 
The ill structure of agriculture and land use makes the level of 
        economic benefit low.
 
There exists barren mountain about 40% of the total area, which 
        cases not only the waste of land resources, but also soil and water 
        loss.  
 The establishment of models for regional comprehensive 
      development and prediction
 The general structure of the analysis 
      system based on GIS is shown in figure 1.
 
 
 There are four models-the 
      population prediction models (PPM), the farming structure dynamic 
      optimization goal programming model (FSDO-GPM), the forest growth 
      prediction model (FGPM) and the regional comprehensive development system 
      dynamics model (RCDSDM). In RCDSDM, Eight Subsystems, Farming, Forestry, 
      animal-husbandry, living energy, population, industry, investment and soil 
      and water loss, are considered. The cause and effect relation in RCDSDM is 
      shown in figure 2.
 
 
 
 
 The results of the models 
      running
 Through the analysis of several simulation plans, it is 
      thought that the feasible scheme for regional comprehensive development 
      should be:
 
        The models expands 
      over 1987 to 2000. The prospects of the socio - economic development in 
      the country could be:To promote the proper investment of farming to raise the 
        productivity of farmland and make farming production steadily 
        grow.
 
To return the sloppy farmland (unsuitable for farming) to forest at 
        a prerequisite that the production of grain meets the demand of food for 
        the predicted population.
 
To plant trees and grass in barren mountain area, especially to 
        plant economic trees and fire woods, so as to promote the economic 
        benefit and meet the needs of living energy.  
        With the returning of unsuitable sloppy farmland to forest, the 
        farmland decreases, which causes the slow increase of farming output 
        value
 
Through planting trees and grass on Barren Mountain, especially the 
        economic trees and fire woods, the forestry and animal -husbandry output 
        value increases rapidly.
 
Concerning the total agricultural value, agriculture develops 
        steadily and the average value increases rapidly.  The productive data are shown in Table 1.  
        
        
          | Landuse AreaYear | Farmland | Sloppy farm land | Terraced farm land | Forest land returned from sloppy farm land | Economi trees land | Shelter and wood forest land | Fire wood forest land | Grassland | Barren mountalu | Out put value a griculture |  
          | 1987 | 911.6 | 363.4 | 171.1 | 0 | 21.8 | 1255 | 268.5 | 171.9 | 2080 | 151.5 |  
          | 1988 | 894.0 | 329.5 | 187.5 | 17.6 | 44.5 | 1315 | 298.0 | 246.0 | 1911 | 163.1 |  
          | 1989 | 877.3 | 297.9 | 202.3 | 16.6 | 67.2 | 1375 | 327.5 | 316.5 | 1755 | 174.9 |  
          | 1990 | 861.3 | 268.6 | 215.7 | 16.1 | 89.5 | 1432 | 356.7 | 380.8 | 1613 | 186.4 |  
          | 1991 | 846.0 | 241.2 | 227.8 | 15.2 | 111.0 | 1486 | 385.3 | 437.5 | 1485 | 197.5 |  
          | 1992 | 831.5 | 215.8 | 238.7 | 14.6 | 131.0 | 1536 | 413.0 | 487.1 | 1371 | 208.1 |  
          | 1993 | 817.6 | 192.3 | 248.4 | 13.8 | 151.0 | 1582 | 440.0 | 529.9 | 1269 | 218.2 |  
          | 1994 | 804.4 | 170.4 | 257.0 | 13.2 | 169.4 | 1624 | 466.0 | 566.6 | 1179 | 227.6 |  
          | 1995 | 791.8 | 150.1 | 264.7 | 12.6 | 186.5 | 1661 | 491.1 | 597.8 | 1099 | 236.4 |  
          | 1996 | 779.8 | 131.4 | 2715 | 12.1 | 202.5 | 1695 | 515.2 | 624.2 | 1028 | 244.6 |  
          | 1997 | 768.4 | 114.0 | 277.4 | 11.4 | 217.3 | 1725 | 538.5 | 646.5 | 965.2 | 252.1 |  
          | 1998 | 757.5 | 98.0 | 282.5 | 10.9 | 230.9 | 1752 | 561.0 | 665.4 | 909.2 | 259.1 |  
          | 1999 | 747.1 | 83.2 | 286.9 | 10.4 | 243.4 | 1777 | 582.6 | 681.3 | 859.4 | 265.5 |  
          | 2000 | 737.2 | 69.6 | 290.7 | 9.9 | 254.9 | 1798 | 603.5 | 694.7 | 815.0 | 271.3 | Spatial planning of 
      predictive data of land use types
 Table 1 gives various predictive 
      data of land use types and agricultural output value from 1987 to 2000. 
      However, these predictive data neither tell us where they should be 
      planned. nor answer how the ecological environment will become after 
      suitable measures have been take. Therefore, it is necessary to plan the 
      predictive data onto the suitable geographical space.
 
 Having 
      noticed that, spatial planning of predictive data depends on the regional 
      natural environmental conditions and socio-economic conditions, the 
      planning process is divided into two steps:
 
 
 
        Spatial Planning Constrained by Natural Environmental 
        Conditions
 
 
          Land Resource Evaluation: The fuzzy mathematical model is applied 
          to evaluate land resources in Pingqual. With this model, the resources 
          can be divided into four types farming favourable land, forestry 
          favourable land, animal-husbandry favourable land and unfavorable 
          land. In farming and forestry favourable land, three land quality 
          classes high, medium and low-are classified respectively. Each has 
          special land suitability and land use purpose. The land resource 
          evaluation map is shown in Figure 3.
 
Creation of the planning Goal Map: The planning goal map is 
          created by overlaying land use map and land evaluation map based on 
          Geo-Union, which provides a set of map operations for planning of 
          predictive data.
 
The Establishment of Land Use Suitability Table (Table 2) 
 
 
          Table2. The Land Suitability Table for Planning  
            
            
              | Land gult-ablity types | Land quality classes | Land use direction | Constraint conditions |  
              | Elevation (m) | Slope | Soil erosion | Soil orange matter content | Soil texture |  
              | Farming favourable land | High class | Rice and wheat crop | 335-450 | 00 -30 | 0 | 0-1 | 0-1 |  
              | Medium class | Course grain crop | 335-500 | 00- 70 | 0-1 | 0-1 | 0-1 |  
              | Low class | Coarse grain crop to be improvement | 335-700 | 00- 00 | 0-2 | 0-1 | 0-1 |  
              | Forestry favourable land | High class | Economic tree | 335-1000 | 00= 50 | 0-3 | 0-2 | 0-2 |  
              | Medium class | Wood forest | 335-1800 | 00= 00 | 0-3 | 0-2 | 0-3 |  
              | Low class | Firewood | 335-1800 | 00-50 | 0-4 | 0-2 | 0-4 |  
              | Animal Hubandary favourable Land | Grass | 335-1800 | 00-50 | 0-4 | 0-3 | 0-5 |  
              | Unfavourable land |  | 335-1800 | > 450 | 5 | 4 | 6 |  
Spatial planning the basic principle is to search for proper land 
          according to the land suitability table. Generally, higher quality 
          land has greater priority in planning. In order to plan the predictive 
          data reasonably, the fuzy score is calculated for each planning unit 
          under the constraint of natural environmental factors. The formula of 
          fuzzy comprehensive score model is 
 Vi =P1* Ui1 + P2 + P3* Ui3 
          + P4*Ui4=P5*Ui5………..(A)
 
 Uij =1/[ 1+AJ* ( Wij -Cj) (Wij -Cj)] 
          …………..( B)
 
 Where Vi is the fuzzy score of planning unit, Pj is 
          the weight of constraining factor. Formula (B) is the fuzzy 
          mathematical function of each factor.
 
 Taken DTM grid as 
          planning unit , the predictive data are planned in the planning goal 
          may by searching for proper area grid by grid according to the fuzzy 
          score of each grid . When the searched area is equal to the value of a 
          predictive datum, the searching work stops and the attribute of the 
          predictive datum is given to the searched area in planning goal map.
 
Spatial Planning Constrained by Socio-Economic Conditions
 In 
        the spatial planning, socio-economic conditions should also be 
        considered. For example, we should not plant economic trees in a small 
        portion of land where it is very difficult to access, even though it is 
        very suitable to economic trees. In this paper, the production and 
        management conditions to farming and economic trees concerned. The 
        processes are as follows:
 
 Overlaying road distribution map, resident distribution map with 
          the original planning maps.
 
Searching for those unsuitable farming land (or economic tree land 
          ) polygons that the area of each polygon is less than D (where D is a 
          constant ).
 
Beginning with the smallest polygon, calculating the core point of 
          each polygon.
 
Creating a changeable circle by taking the core point as its 
          centre and r as its radius (where r is a variable, it increases from 
          zero). The circle gradually expands with the increase of r. When the 
          are of the circle hits the nearest boundary of any road, resident or 
          farmland (or economic tree land) the circle stops expanding and the 
          value R of the radius is written down. If R > L (where L is a 
          constant) , it is thought that the polygon is not suitable for farming 
          (or economic trees) and it is merged to other type. 
        
 Similarly, polygon by polygon, all the unsuitable farmland 
        (or economic trees) is filtered through step I and step II, the 
        predictive data in Table 1 are planned in geographical space, and a 
        series of predictive thematic maps from 1987 to 2000 , such as land use 
        maps, forest distribution maps and plant coverage distribution maps , 
        etc. are created . The land use maps of two periods, 1987 and 2000, are 
        sleeted to the shown in figure 4 and Figure 5.
 
 The prediction 
        of soil and water loss changes
 According to the land use maps 
        and plant coverage maps, the changing tendency of soil and water loss 
        has been studied by applying the soil and water common equation. The 
        formula is
 
 
 A= R*K*L*S*C*P  Where A is the amount of soil 
        erosion, R is the rainfall factor, K is the soil property factor, L is 
        the length factor of sloping land . So is the gradient factor of sloping 
        land, C is biological factor and P is the protective measure factor.
 
 The soil erosion distribution maps in 1987 and 2000, are sleeted 
        to be shown in figure 6 and Figure 7. The tendency curve of soil erosion 
        from 1983 to 2000 is shown in Figure 8.
 
 Reference
 
 
 
          Jay W. Forrester, Industrial Dynamic, Mass : The MIT Press,1980 
          ZK)
 
Burrough, P.A. Principles of Geographical Information Systems for 
          Land Resource Assessment, Oxford University Press, 1986.ZK)
 
Walsh, S.J. Geographical Information Systems for Natural Resource 
          Management, Journal of Soil and Water Conservation, Vol.40, 202-205, 
          1985.  |