doi:10.1016/j.jenvman.2006.02.005
Copyright © 2006 Elsevier Ltd All rights reserved.
A quantitative method for zoning of protected areas and its spatial ecological implications
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
María del Carmen Sabatinia,
, Adriana Verdiellb,
,
, Ricardo M. Rodríguez Iglesiasa and Marta Vidalb, c
aCONICET and Departamento de Agronomía. Universidad Nacional del Sur. San Andrés 800, (8000) Bahía Blanca, Argentina
bDepartamento de Matemática. Universidad Nacional del Sur. Av. Alem 1200 (8000), Bahía Blanca, Argentina
cUniversidad Tecnológica Nacional, Facultad Regional Bahía Blanca. 11 de Abril 461 (8000), Bahía Blanca, Argentina
Received 4 December 2004;
revised 2 January 2006;
accepted 16 February 2006.
Available online 11 May 2006.
Abstract
Zoning is a key prescriptive tool for administration and management of protected areas. However, the lack of zoning is common for most protected areas in developing countries and, as a consequence, many protected areas are not effective in achieving the goals for which they were created. In this work, we introduce a quantitative method to expeditiously zone protected areas and we evaluate its ecological implications on hypothetical zoning cases. A real-world application is reported for the Talampaya National Park, a UNESCO World Heritage Site located in Argentina.
Our method is a modification of the zoning forest model developed by Bos [Bos, J., 1993. Zoning in forest management: a quadratic assignment problem solved by simulated annealing. Journal of Environmental Management 37, 127–145.]. Main innovations involve a quadratic function of distance between land units, non-reciprocal weights for adjacent land uses (mathematically represented by a non-symmetric matrix), and the possibility of imposing a connectivity constraint. Due to its intrinsic spatial dimension, the zoning problem belongs to the NP-hard class, i.e. a solution can only be obtained in non-polynomial time [Nemhausser, G., Wolsey, L., 1988. Integer and Combinatorial Optimization. John Wiley, New York.]. For that purpose, we applied a simulated annealing heuristic implemented as a FORTRAN language routine.
Our innovations were effective in achieving zoning designs more compatible with biological diversity protection. The quadratic distance term facilitated the delineation of core zones for elements of significance; the connectivity constraint minimized fragmentation; non-reciprocal land use weightings contributed to better representing management decisions, and influenced mainly the edge and shape of zones.
This quantitative method can assist the zoning process within protected areas by offering many zonation scheme alternatives with minimum cost, time and effort. This ability provides a new tool to improve zoning within protected areas in developing countries.
Keywords: Zoning; Protected areas; Simulated annealing
Fig. 1. Hypothetical case of land use aptitude and its surrounding area. Land aptitude is depicted using a 0 (no aptitude, white) or 1 (with aptitude, black) for extensive recreation and strict preservation.
Fig. 2. Talampaya land use aptitude. Land aptitude scale is from 0 (no aptitude, white) to 9 (maximum aptitude, black).
Fig. 3. Alternative zonings obtained with different assumptions for key parameters: distance between cells, connectivity constraint, and symmetry of compatibility coefficients; see text for further explanation. (a) Squared Euclidean distance, connectivity constraint for the strict preservation use and non reciprocal compatibility between uses (extensive recreation−strict preservation =+1, strict preservation−extensive recreation =−1); other parameters as in Table 2. (b) Without connectivity constraint for strict preservation; other parameters as in a). (c) Linear Euclidean distance; other parameters as in (a). (d) Reciprocal compatibility (−1) between extensive recreation and strict preservation; other parameters as in (a). (e) Reciprocal compatibility (+1) between extensive recreation and strict preservation; other parameters as in (a).
Fig. 4. Qualitative (right), and quantitative (left) zoning of Talampaya National Park. The oval encloses an area of marked disagreement between zonings; see text for further explanation. (a) Quantitative zoning. (b) Qualitative zoning by Dellafiore and Sylvester (2000).
Table 1.
Selected metrics for assessing zoning design

See McGarigal et al. (2002) for detailed descriptions and metrics formulae.
Table 2.
Parameters for the hypothetical zoning scenarios

Table 3.
Parameters for Talampaya National Park

Table 4.
Spatial metrics computed for comparing results from alternative hypothetical scenarios using FRAGSTATS
a At whole zoning level.
b At strict preservation use level.
This research was sponsored by UNS Project 24/L057 and CONICET PIP 02156.

Corresponding author. Tel.: +54 291 4595162; fax: +54 291 4595163.