Predicting future dynamics of landscape structure within protected areas using CA-Markov model (Case study: Dizmar protected area)

Document Type : Full length article


1 Ph.D. Candidate in Environmental Planning, Faculty of Environment, University of Tehran, Iran

2 Ph.D. in Environmental Planning, Faculty of Environment, University of Tehran, Iran


Land use and land cover change (LUCC) is a complex issue resulted from biophysical, socio-economic, cultural, organizational and technological factors in different spatial and temporal scales. LUCCs have both direct and indirect effects on the environment not only globally but also locally. LUCC is considered as an important threat to biodiversity as causing the fragmentation, natural vegetation destruction and natural areas isolation. The regions which managed by environmental protection organizations all over the world are established based on a common goal to maintain biodiversity. Current insufficient preservative and management actions in the protected areas (especially in Iran) are unable to guarantee the areas protection. Therefore, analyzing previous and current land use and land cover (LULC) status and predicting the future pattern within and surrounding protected areas are likely to provide more efficient information for proper natural resources management. RS data is cost effective means to detect changes on the Earth's surface and provide up to date information. Over the last decades, several methods and models are developed for extracting LULC maps, detecting LUCCs and modelling the future pattern using remotely sensed data. The objective of this study is to analyze spatiotemporal patterns of LUCC from the past to the future within Dizmar protected area in Iran. Firstly, LULC maps of 1989, 2000 and 2013 were extracted and then future LULC was predicted using CA- Markov models from 2013 to 2037.
Study area
Dizmar protected area is a mountainous-forested region located in the north of Eastern Azerbaijan province, Iran. It lies between the 41◦38' to 57◦38' N and 18◦40' to 46◦40' E with total area about 68576 ha. Its connections to Kiamaky nature reserve in the west, Arasbaran protected area in the east, the national park of Zagatay (in the Republic of Azerbaijan and Armenia) in the north, makes an important wildlife corridor in local, national and international levels.  It is home to 849 plant species (76 of them is endemic) and about 320 species of wildlife (such as Tetrao mlokosiewiczi placed in the IUCN list of globally threatened species). This protected area was faced with extra pressure on natural sources causes LUCC.
Materials and Methods
This study used Landsat satellite images (1989, 2000 and 2013) to extract LULC maps. After preprocessing step (such as image enhancement using Histogram Equalization) unsupervised classification was done. Then, the supervised classification was performed using the Maximum Likelihood Classifier (MLC) based on signatures file (generated from ground reference data that gathered in the field survey) for each of the images separately. Three LULC categories were extracted from TM, ETM+, and OLI images. Stratified random method in ERDAS Imagine 2013 is used to assess the accuracy of each obtained maps. CA­-Markov model was applied to project LULC in the study area for 2037. Validating the LULC prediction model is carried out using KIA (Kappa Agreement Index). LUCCs during studied timespans were calculated using the cross tabulation technique in Idirisi Selva environment.
Results and Discussion
The distribution, coverage and percentage of major LULC types (classified as agricultural land, barren/range land and forestland) for 1989, 2000 and 2013 are shown in Figure 1and Table 1.
Fig. 1. LULC maps of actual (1989, 2000 and 2013) and simulated (2013 and 2037)
The overall classification accuracy of each map for 1989, 2000 and 2013 are estimated to be 89%, 90% and 91%, respectively. The Kappa values also yield 0.81, 0.84 and 0.88, respectively. The main types of LULC was forestland (with 62.20% and 54.30% of the total area) from 1989 to 2000 but it changed to barren/range land in 2013 (52.53% of the total area). 
Results show a reduction in forestland between 1989 and 2013. Subsequently, agricultural land is increased from 0.72% in 1989 to 2.14% in 2013 due to the fact that the traditional livelihood remains farming.
Table 1. Distribution of LULC type in Dizamr protected area (in ha)






LULC type













Agricultural land






Barren/range land







The projected land use map by the CA-Markov model indicates that if the current management continues, barren/range land and agricultural land  reach to 62.03% and 3.22% of the total area at the expense of decreasing forestland area to 34.73% by 2030 (Fig. 2). In order to validate the CA-Markov model outputs the VALIDATE module existing in the IDRISI Selva was used. This is done by comparing simulated land use maps of 2013 with the actual ones based on Kappa statistics. Resulting Kappa values (Kno= 0.9295, Kstandard= 0.918, KlocationStrata= 0.9273 and Klocation= 0.9273) were all greater than 0.9 showing a satisfactory level of accuracy.
Fig. 2. Area of LULC types (as percentage of the total area) over the studied period
The area of 701.46 ha has been deforested and changed into agricultural land during 1989 to 2013. The amount of deforestation will be 521.19 ha by 2013. On the contrary, only 5.22 ha will be forested by transformation of barren/range land during 2013-2037.
Table 2. LULC conversions types during studied time spans (areas in ha)

Times pan

LULC conversions









F to A





F to BR





A to BR





BR to F





BR to A






*A: Agricultural land, BR= Barren/range land, F= Forestland
LUCC within and surrounding the protected areas probably continue to be expanded and intensified. Monitoring and projecting these changes can play key roles in preventing negative consequences of the changes by providing up to date information to planners and managers. This study shows the important role of LUCCs analysis and modeling to provide proper information for the protected area management. We applied dynamic approach to analyze LUCCs by analyzing previous and current LULC maps and predict the future trends. The results indicate the high capability of CA-Markov model to predict future LUCC in the study area. Therefore, it can be useful in the protected area's land use policy and action design.  Indeed, between 1989 and 2013, there has been a notable reduction in forestland and it was predicted to continue the reduction over the next 24 years. Agricultural land has been steady in increment during 1989-2013 and this trend continues by 2037. Expansion of agricultural land and barren/range land in the study area has led to rapid changes in landscape dynamics. Thus, it is recommended to create and strengthen non-farm/off-farm income. Adoption of agricultural policy based on the agroecological condition of the Dizmar protected area is important. Analysis of other factors such as land capability, stakeholders and LUCC drivers along with  the obtained results can be useful in proper LULC planning and management. The strategies of land resources (especially forest resources) development, attempts to overcome the current deterioration and avoid further extinction of remnant forest within the study area.


Main Subjects

ادهمی، س. و خلاقی، س. (1386). مفاهیم پردازش تصویر با تأکید بر نرم‏افزار ERDAS IMAGINE، سبزوار: امید مهر.
امینی، پ.و. (1393). مدل‏سازی اثرات احتمالی تغییر کاربری اراضی پیرامونی بر مدیریت ذخیره‏گاه زیست‌کرة ارسباران، استاد راهنما: احمدرضا یاوری، رشتة برنامه‏ریزی، مدیریت و آموزش محیط‏ زیست، پایان‏نامة کارشناسی ارشد، دانشکدة محیط‏ زیست، دانشگاه تهران.
امینی، پ.و.؛ صالحی، ا.؛ عادلی، ش. و عزیزی، ع. (1394). شبیه‏سازی تغییرات پویای کاربری زمین با استفاده از مدل تلفیقی CA-MARKOV (بررسی موردی: شهرستان ملکان)، فصلنامة علوممحیطی، 13(3): 133 ـ 142.
بنفشه، م.؛ رستم‌زاده، ه. و فیضی‌زاده. ب. (1386). بررسی و ارزیابی روند تغییر سطوح جنگل با استفاده از سنجش از دور و GIS (مطالعة موردی جنگل‏های ارسباران 1997 ـ 2005)، پژوهش‏هایجغرافیایی، 39(62): 143 ـ 159.
سازمان حفاظت محیط‏زیست (1391). گزارش مطالعة توجیهی منطقة حفاظت‏شدة دیزمار، آذربایجان شرقی.
نژادی، ا. (1391). تدوین سامانة پشتیبان تصمیم‏گیری مدیریت مناطق حفاظت‏شده بر مبنای مدل‏سازی تغییرات کاربری اراضی، مطالعة موردی: منطقة حفاظت‏شدة لیسار، استاد راهنما: حمیدرضا جعفری، رشتة برنامه‏ریزی محیط‏زیست، رسالة دکتری، دانشکدة محیط‏ زیست، دانشگاه تهران.
هاشمی، ن. (1392). شناسایی نواحی اولویت‏دار حفاظتی با رهیافت مدل‏سازی تغییر کاربردی اراضی (مطالعة موردی: منطقة حفاظت‏شدة جاجرود)، استاد راهنما: یاوری ا، ر.، رشتة برنامه‏ریزی، مدیریت و آموزش محیط‏ زیست، پایان‏نامة کارشناسی ارشد، دانشکدة محیط‏ زیست، دانشگاه تهران.
Adhami, S. and Khalaghi, S. (2007). Concepts of image processing in Erdas Imagine, Omid Mehr Press, Sbazevar, Iran.
Alain, O.G. (2014). IDRISI Selva Tuterial, Clark Labs, Clark University, United States.
Alers, M.; Bovarnick, A.; Boyle, T.; Mackinnon, K. and Sobrevila, C. (2007). Reducing Threats to Protected Areas Lessons from the Field, United Nations Development Programme, UN Plaza, New York.
Almatar, M. (2011). Utilizing geographic information systems and remote sensing to investigate urbanization processes: in both the US and KUWAIT. Phd thesis, University of Florida.
Amini, P.V. (2014). Modeling Plausible Impacts of Land Use Changes of the Surrounding Buffer Zone on the Management of the Arasbaran Biosphere Reserve, MSc thesis, Supervisor, Yavari, A, R., Environmental Planning, manangement and education, Faculty of Environment, University of Tehran, Iran.
Amini, P.V.; Salehi, E.; Adeli, S. and Azizi. A. (2015). Simulation of Land Use Change Dynamics Based on the CA-Markov Model (Case Study: Malekan County, Iran), Journal of Environmental Sciences,  13(3): 133-144.
Balzter, H. (2000). Markov chain models for vegetation dynamics, Ecological Modelling, 126 (2-3): 139-154.
Banafshe, M.; Rostamzadeh, H. and Feyzi Zadeh, B. (2007). Studying and analyzing forest area change trends using RS and GIS (case study: Arasbaran forests, 1987-2005), Physical Geography Research Quarterly, 39(62): 143-159.
Barrett, C.B.; Brandon, K.; Gibson, C.C. and Gjertsen, H. (2001). Conserving tropical biodiversity amid weak institutions, Bioscience, 51: 497-502.
Bates, D. and Rudel, T.K. (2000). The political ecology of conserving tropical rain forests: A cross national analysis, Society & Natural Resources,  13(7): 619-634.
Department of Environment (2012). Baseline studies of Dizmar protected Area, East Azerbaijan, Iran.
Foody, G.M. (2002). Status of land cover classification accuracy assessment, Remote sensing of environment,  80(1): 185-201.
Gaston, K.J.; Charman, K.; Jackson, S.F.; Armsworth, P.R.; Bonn, A.; Briers, R.A.; Callaghan, C.S.Q.; Catchpole, R.; Hopkins, J.; Kunin, W.E.; Latham, J.; Opdam, P.; Stoneman, R.; Stroud, D.A. and Tratt, R. (2006). The ecological effectiveness of protected areas: the United Kingdom, Biological Conservation,  132(1): 76-87.
Geist, H.J. and Lambin, E.F. (2002). Proximate causes and underlying driving forces of tropical deforestation, BioScience,  52(2): 143-150.
Hashemi, N. (2014). Determining priority of protection areas using land use changes modeling approach, case study: Jajrood protected area, MSc thesis, Supervisor, Yavari, A, R., Environmental Planning, manangement and education, Faculty of Environment, University of Tehran, Iran.
Houet, T. and Hubert-Moy, L. (2006). Modelling and projecting landuse and land-cover changes with a cellular automaton in considering landscape trajectories: an improvement for simulation of plausible future states, EARSeL eProceedings,  5: 63-76.
Kharouba, H.M. and Kerr, J.T. (2010). Just passing through: Global change and the conservation of biodiversity in protected areas, Biological Conservation,  143(5): 1094-1101.
Koomen, E. and Beurden, J.B-V. (2011). Land-use modelling in planning practice, The GeoJournal Library, Springer.
Lambin, E.F. and Geist, H.J. (2006). Land-use and land-cover change Local Processes and Global Impacts, Springer-Verlag Berlin Heidelberg.
Lambin, E.F.; Geist, H.J. and Lepers, E. (2003). Dynamics of land-use and land-cover change in tropical regions, Annual Review of Environment and Resources,  28: 205-241.
Martínez, M.L.; Pérez-Maqueo, O.; Vázquez, G.; Castillo-Campos, G.; García-Franco, J.; Mehltreter, K.; Equihua, M. and Landgrave, R. (2009). Effects of Land Use Change on Biodiversity and Ecosystem Services in Tropical Montane Cloud Forests of Mexico, Forest Ecology and Management,  258(9): 1856-63.
Nejadi, A. (2012). Developing a Decision Support System to manage the protected areas based on land use change modeling Case study: Lisar protected area. Ph.D thesis, Supervisor, Jafari, H., Faculty of Environment, University of Tehran, Iran.
Nejadi, A.; Jafari, H.R.; Makhdoum, M.F. and Mahmoudi, M. (2012). Modeling plausible impacts of land use change on wildlife habitats, application and validation: Lisar protected area, Iran, International Jurnal of Environmental Resech,  6: 883-892.
Nelson, G.C.; Bennett, E.; Berhe, A.A.; Cassman, K.; DeFries, R.; Dietz, T.; Dobermann, A.; Dobson, A.; Janetos, A.; Levy, M.; Marco, D.; Nakicenovic, N.; O’Neill1 B.; Norgaard, R.; Petschel-Held, G.; Ojima, D.; Pingali, P.; Watson, R. and Zurek, M. (2006). Anthropogenic drivers of ecosystem change: an overview. Ecology and Society, 11 (2): 29.
Pouzols, M.F.; Toivonen, T.; Di Minin, E.; Kukkala, A.S.; Kullberg, P.; Kuusterä, J.; Lehtomäki, J.; Tenkanen, H.; Verburg, P.H. and Moilanen, A. (2014). Global protected area expansion is compromised by projected land-use and parochialism, Nature,  516: 383-386.
Sang, L.; Zhang, C.; Yang, J.; Zhu, D. and Yun, W. (2011). Simulation of land use spatial pattern of towns and villages based on CA-Markov model, Mathematical and Computer Modelling,  54: 938-943.
Subedi, P.; Subedi, K. and Thapa, B. (2013). Application of a Hybrid Cellular Automaton – Markov (CA-Markov) Model in Land-Use Change Prediction: A Case Study of Saddle Creek Drainage Basin, Florida, Applied Ecology and Environmental Sciences, 16: 126-132.
Turner, W.R.; Brandon, K.; Brooks, T.M.; Costanza, R.; da Fonseca, G.A.B. and Portela, R. (2007). Global conservation of biodiversity and ecosystem services, BioScience, 57: 868-873.
Verburg, P.H.; Overmars, K.P.; Huigen, M.G.A.; de Groot, W.T. and Veldkamp, A. (2006). Analysis of the effects of land use change on protected areas in the Philippines, Applied Geography, 26: 153-173.
Viera, A.J. and Garrett, J.M. (2005). Understanding interobserver agreement: The kappa statistic, Family Medicine,  37(5): 360-363.
Vliet, J. Van. (2009). Assessing the Accuracy of Changes in Spatial Explicit Land Use Change Models, 12th AGILE International Conference on Geographic Information Science, Leibniz Universität Hannover, Germany, PP. 1-9.
Wilson, T.S.; Sleeter, B.M. and Davis, A.W. (2015). Potential future land use threats to California’s protected areas, Regional Environmental Change,  15: 1051-1064.