# Comparison of Marine Spatial Planning Methods in Madagascar Demonstrates Value of Alternative Approaches

Thomas F. Allnutt, Timothy R. Mcclanahan, Serge Andrefouet, Merrill Baker, Erwann Lagabrielle, Caleb Mcclennen, Andry J. M. Rakotomanjaka, Tantely F. Tianarisoa, Reg Watson, Claire Kremen

<figure><img src="image_1.png" /></figure>

## To cite this version:

Thomas F. Allnutt, Timothy R. Mcclanahan, Serge Andrefouet, Merrill Baker, Erwann Lagabrielle, et al.. Comparison of Marine Spatial Planning Methods in Madagascar Demonstrates Value of Alternative Approaches. PLoS ONE, 2012, 7 (2), pp.e28969. 10.1371/journal.pone.0028969. hal-01239639

HAL Id: hal-01239639

https://univ-reunion.hal.science/hal-01239639v1

Submitted on 14 Jun 2018

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
# Comparison of Marine Spatial Planning Methods in Madagascar Demonstrates Value of Alternative Approaches

Thomas F. Allnutt¹*, Timothy R. McClanahan², Serge Andréfouët³, Merrill Baker¹, Erwann Lagabrielle⁴,⁵, Caleb McClennen², Andry J. M. Rakotomanjaka⁶, Tantely F. Tianarisoa⁶, Reg Watson⁷, Claire Kremen¹

1 Department of Environmental Science, Policy and Management, University of California, Berkeley, California, United States of America, 2 Marine Programs, Wildlife Conservation Society, Bronx, New York, United States of America, 3 Institut de Recherche pour le Développement, Nouméa Cedex, New Caledonia, 4 Institut de Recherche pour le Développement, Sainte-Clotilde, France, 5 Botany Department, Nelson Mandela Metropolitan University, George, South Africa, 6 Madagascar Program, Wildlife Conservation Society, Antananarivo, Madagascar, 7 Fisheries Center, University of British Columbia, Vancouver, Canada

## Abstract

The Government of Madagascar plans to increase marine protected area coverage by over one million hectares. To assist this process, we compare four methods for marine spatial planning of Madagascar's west coast. Input data for each method was drawn from the same variables: fishing pressure, exposure to climate change, and biodiversity (habitats, species distributions, biological richness, and biodiversity value). The first method compares visual color classifications of primary variables, the second uses binary combinations of these variables to produce a categorical classification of management actions, the third is a target-based optimization using Marxan, and the fourth is conservation ranking with Zonation. We present results from each method, and compare the latter three approaches for spatial coverage, biodiversity representation, fishing cost and persistence probability. All results included large areas in the north, central, and southern parts of western Madagascar. Achieving 30% representation targets with Marxan required twice the fish catch loss than the categorical method. The categorical classification and Zonation do not consider targets for conservation features. However, when we reduced Marxan targets to 16.3%, matching the representation level of the "strict protection" class of the categorical result, the methods show similar catch losses. The management category portfolio has complete coverage, and presents several management recommendations including strict protection. Zonation produces rapid conservation rankings across large, diverse datasets. Marxan is useful for identifying strict protected areas that meet representation targets, and minimize exposure probabilities for conservation features at low economic cost. We show that methods based on Zonation and a simple combination of variables can produce results comparable to Marxan for species representation and catch losses, demonstrating the value of comparing alternative approaches during initial stages of the planning process. Choosing an appropriate approach ultimately depends on scientific and political factors including representation targets, likelihood of adoption, and persistence goals.

**Citation:** Allnutt TF, McClanahan TR, Andréfouët S, Baker M, Lagabrielle E, et al. (2012) Comparison of Marine Spatial Planning Methods in Madagascar Demonstrates Value of Alternative Approaches. PLoS ONE 7(2): e28969. doi:10.1371/journal.pone.0028969

**Editor:** Richard K. F. Unsworth, Swansea University, United Kingdom

**Received** September 23, 2011; **Accepted** November 18, 2011; **Published** February 16, 2012

**Copyright:** © 2012 Allnutt et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by MacArthur Foundation grants 07-89636-000-GSS and 07-89630-001-GSS. The fish catch model was developed from the "Sea Around Us" global catch database. This database was developed through a scientific collaboration between the University of British Columbia and the Pew Environment Group. The Indian Ocean Commission, Fond Français pour l'Environnement Mondial and World Wildlife Fund funded the Indian Ocean Commission Protected Areas Network conservation planning workshop. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: tom.allnutt@berkeley.edu

## Introduction

While climate change, overfishing and land use change negatively impact the biodiversity and ecological function of marine ecosystems worldwide [1,2], increasing evidence shows that effective conservation and management can recover the resource base, conserve biodiversity, and increase fisher's incomes [3,4,5,6,7]. Yet targeting interventions to "maximize returns" across large regions remains a major challenge: in most tropical seas, for example, basic patterns and interactions among key biological, environmental and social variables remain poorly understood [8]. In addition, while appropriate targeting methods and data sources may exist, they can be underutilized due to a lack

of real-world examples and comparative evaluations [9]. In this paper we compare four alternative methods of successive technical complexity for identifying conservation and management priorities across Madagascar's west coast, a regionally and globally important tropical coral reef ecosystem [10,11].

The waters of the West Coast of Madagascar are home to 90% of Madagascar's coral reefs, large-scale export fisheries for shrimp, octopus, sea cucumbers, and tuna, and important artisanal fisheries. Nonetheless, the formal management of marine resources in the region is in its initial stages. There is only one fully decreed marine protected area (MPA) (Sahamalaza-Isles Radama), with several others in various stages of designation. Less than 1% of Madagascar's reefs are included in no-take areas, the lowest rate

<figure><img src="image_3.png" /></figure>
for five West Indian Ocean countries [12]. Few additional areas are under formal marine management across a region spanning thousands of kilometers of coastline and home to nearly two million people, many of whom are dependent on marine resources as an important protein source for local consumption and as a source of cash from export or sale [13,14]. Due to inadequate marine resource management across the majority of the region, fishers target lower trophic levels [15], and the use of illegal and destructive gear such as small-meshed beach seine nets is widespread [16]. As a result, large areas of the region's coral reef ecosystems are chronically stressed [17].

Several governmental, non-governmental and community organizations in Madagascar are interested in increasing the scale of marine resource management on Madagascar's west coast through MPAs and other strategies [18]. In general, increasing the scope and effectiveness of marine resource management across large areas requires a large-scale spatial synthesis of biological, socio-economic and environmental patterns [19]; such a synthesis is lacking across Madagascar's west coast, since the region is large, relatively inaccessible and insufficiently surveyed [11].

Many methods now exist for mapping conservation priorities and management actions in marine systems. Collectively termed “Marine Spatial Planning” [20], this approach may involve traditional techniques such as biogeographic classification [21], gap analysis [22], and “systematic” target-based planning methods [23,24,25]. When formulating marine spatial plans, several authors emphasize the need to include socioeconomic as well as biological data [26,27] and to consider exposure to climate change in the context of persistence of managed areas [28,29].

Existing examples and reviews of available planning approaches are biased towards marine systems occurring in more economically developed, data rich countries [9,21,30], and often focus exclusively on quantitative efficiency measures to identify an “optimal” result [31]. In the planning context within Madagascar, however, uncertainties surrounding potential social and ecological responses to management and threats could lead to lost opportunities, such as sustainable management of fisheries, if there is heavy reliance on a single optimization approach in early planning stages. Therefore, we apply multiple conceptual and analytic techniques to the same data sources to compare results, and provide alternative starting points for developing regional marine conservation and management plans.

# Methods

## Study area

Our study area is the coastal region of western Madagascar, from Cap Vohimena (Cap Sainte Marie) in the south, to Cap d'Ambre in the north. The region includes extensive fringing, patch and barrier reefs, many small islands, large areas of coastal shelf, mangroves, seagrass beds, and other typical tropical marine ecosystems. Measured on a 1:250,000 scale [32], the region has 7,000 km of coastline across a 14° latitudinal north-south gradient, and exhibits considerable variation in environmental conditions [33,34]. To define the study region, we first combined the boundaries of the three coastal bioregions that occur along Madagascar's west coast [18]. These bioregions cover the entire neritic zone from the coastline to a depth of roughly 200 m. As a final step defining the study area, we extended the outer boundary by 0.25 degrees to ensure complete coverage of coastal and neritic habitats, producing a total area of 201,057 km² (Figure 1).

Unlike Madagascar's east coast, most of the west coast has a wide, shallow, gently sloping continental shelf. At the widest point, between Besalampy and Juan de Nova Island, the region is

approximately 175 km wide. The narrowest zones are offshore from Toliara province, where faulting has produced a steeper, narrower shelf typically only 40–50 km wide. Ninety-two percent of the study area is within Madagascar's Exclusive Economic Zone (EEZ), adjacent to Mahajanga and parts of Toliara, and Antsiranana provinces. Small portions also fall within the Juan de Nova and Glorieuses EEZ, both territories of France (Figure 1).

According to a recent study of reef geomorphology [35], 90% of Madagascar's coral reefs are found within this study region. These reefs support the highest coral [36] and invertebrate macrofauna [37] species richness in the central and western Indian Ocean, and also have high levels of mollusk and fish diversity [38]. The study area is also important for marine megafauna, including significant populations of five species of marine turtles [39], 27 species of marine mammals, including several threatened species such as dugong (*Dugong dugon*) [40,41], coelacanth (*Latimeria chalumnae*) [42], as well as a number of regionally important seabird colonies [43], especially for Sooty (*Onychoprion fuscata*) and Crested Terns (*Thalasseus bergii*) [44]. Overall, however, the Indian Ocean is the least known of the tropical seas [45], and the biodiversity of the study area remains poorly surveyed [11]. The large majority of readily available biological data comes from the extensive and now largely degraded reef systems off of Toliara, in the south [46]. The next best-studied area is centered on the island of Nosy Be, 1,200 km to the north. Relatively little is known about the vast area between these two widely separated locations, although recent surveys have begun to specifically target these gaps [47].

Human and environmental impacts on the region are also not well studied but are likely to be highly variable, dependent on the intensity of fishing for subsistence or export, sedimentation from deforestation, and the strength of relatively recent severe sea temperature anomalies [33,34,48]. Threats to marine resources are changing with the rapid increase of coastal populations. Two of the west coast's three provinces, Tulear and Mahajanga, now have the highest birth rates in all of Madagascar: 6.3 and 6.1 children per women, respectively [14]. In some areas along the west coast, numbers of boats and fishermen have increased fivefold since the 1980's, contributing to resource decline, particularly in waters near urban centers such as Toliara [17,46]. Export-driven demand for products such as octopus, sea cucumbers and shark fins is another important driver of ecological impacts. Many of these fisheries show signs of overexploitation, decline, or transition from high- to low-value species [49,50,51].

Along the West Coast, Sahamalaza-Isles Radama is the only park established primarily for marine protection that is also fully decreed as an “existing” protected area as stipulated under the Madagascar Protected Areas Law [52,53]. In addition to this, there are four fully decreed terrestrial protected areas that include some marine or coastal habitats (Kirindy Mitea, Baie de Baly, Lokobe and Nosy Ve). Madagascar protected areas law also includes “temporary” and “new” designations. “Temporary” protected areas have completed an initial administrative process and await final legislative decree to become full protected areas. For “new” protected areas, the administrative procedure for implementation is ongoing [53]. Nosy Hara is the most significant marine area under “temporary” protected status. Since 2009, ten other areas have been designated as “new” protected areas. Along the west coast in total, there are 5,097 km² (9.2% marine, 90.8% terrestrial) of “existing” protected areas, 14,615 km² (10.6% marine, 89.4 km² terrestrial) of “temporary” protected areas, and 4,579 km² (74.8% marine, 25.2% terrestrial) of “new” protected areas. For our study area, this amounts to 0.3% designated as “existing”, 0.9% as “temporary”, and 2.0% as

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Map of study area on Madagascar's West Coast, and locations mentioned in the text. Study area is shown in grey with black outline. Most of the study area is in Madagascar's Exclusive Economic Zone with the exception of small areas that fall in Glorieuses and Juan de Nova. doi:10.1371/journal.pone.0028969.g001</figcaption></figure>

"new" protected areas. The terrestrial portions of reserves are excluded from these percent coverage statistics.

## Spatial data

We collected spatial data for three variables relevant to conservation decisions across the study area: fishing pressure, exposure to thermal stress, and biodiversity (Table 1).

**Fishing pressure.** Although there are many human activities in the region, such as tourism and marine transport, the human use with the most significant and direct impact on marine resources in the West Indian Ocean is fishing [45]. Here, we mapped fishing pressure as a combination of motorized and non-motorized coastal

fishing pressure (Figure 2a) using a recent model of anthropogenic drivers of marine change for the West Indian Ocean [54]. This model combines spatial data on global fisheries catches [55,56], tuna purse seine catch data supplied by the Indian Ocean Tuna Commission (IOTC) and data on coastal fisheries derived from national fisheries statistics and population density data.

**Exposure to thermal stress.** Of particular interest is the relative susceptibility of shallow tropical marine regions to heat stress that can cause coral bleaching, particularly under future climate scenarios. For this measure, we used results from Maina et al. [33] to show relative probability of exposure to thermal stress across the study area (Figure 2b). Here, "environmental

<figure><img src="image_6.png" /></figure>
<table><caption>**Table 1.** Data used in the study (I), sources, methods (M), and an overview of the comparisons (C).</caption><thead><tr><th></th><th></th></tr></thead><tbody><tr><td>Input variables (I)</td><td>Fishing pressure (catch) (F) [54]<br>Exposure to thermal stress (E) [33]<br>Biodiversity features (B)<br>B1. Bioregions [18]<br>B2. Coral reef geomorphology [35]<br>B3. Mangroves [58,59]<br>B4. Fish species distributions<br>B5. Biodiversity value calculated with Zonation algorithm<br>B6. Biological value calculated as species and habitat richness</td></tr><tr><td>Analysis methods (M)</td><td>M1a. Visual gradient overlay of inputs E, F and B5<br>M1b. Visual gradient overlay of inputs E, F, and B6<br>M2. Categorical classification of conservation and management actions from inputs E, F, and B5<br>M3a. Marxan with 30% feature targets from inputs E, F, and B1-B4<br>M3b. Marxan with 16% feature targets from inputs E, F, and B1-B4<br>M4. Zonation with inputs B1-B4 as positively weighted features, and E, F as negatively weighted opposing features</td></tr><tr><td>Comparisons (C)</td><td>C1. Visual comparisons between M1a, M1b, M2 and M3a<br>C2. Quantitative comparisons between M2 (strictest protection category only), M3a, M3b, M4 (threshold at rank with 16.3% average representation):<br>- Total overlap<br>- Average proportion of distributions included (B1-B4)<br>- Fishing pressure in terms of total catch (F)<br>- Exposure to thermal stress (E): Average value, Probability of result missing feature targets</td></tr></tbody></table>

doi:10.1371/journal.pone.0028969.t001

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Non-biological inputs to the analysis. A: Fishing pressure from fish catch model [54], units are metric tons/km²/yr.; B: Environmental exposure probability [33].
doi:10.1371/journal.pone.0028969.g002</figcaption></figure>

<figure><img src="image_8.png" /></figure>
exposure", a predictor of the degree to which coral communities
are susceptible to climate induced thermal stress, consists of a
weighted combination of nine satellite-derived environmental
variables [33].

**Biodiversity.** We collected spatial distribution data on four biodiversity features: bioregions, coral reefs, mangroves and fish species distributions (Figure 3, Figures 4a,b, Table 2). Bioregions are broad units of relatively similar biological and environmental conditions, frequently used to ensure representation within reserve networks [57]. We utilized a map of bioregions developed from species data, environmental data and expert opinion from a recent West Indian Ocean conservation planning exercise [18]. This map recognizes three coastal bioregions across the entire neritic zone of western Madagascar's coast: a north-western unit from Cap d'Ambre to 15 km north of the city of Mahajanga; a large central unit extending from north of Mahajanga to 150 km south of Toliara; and finally a southern unit continuing from south of Toliara around Cap Vohimena. We created a map of mangroves from a simple union of two sources: Harper et al. [58] and Moat & Smith [59]. Both authors used Landsat TM imagery from 1998 to 2005 to map mangroves at approximately 30 m resolution. We took reef geomorphological types from a West Indian Ocean map of reef morphology produced by Andréfouët et al. [35] from Landsat imagery acquired between 1999 and 2003. The reef data includes 16 reef geomorphological types within the study area.

To map fish distributions, we started with a list of 530 focal fish
species identified by Malagasy and international marine conser-
vation experts as conservation priorities in Madagascar. From this
list we collected available locality data from private and public
sources (Table 2), including the Global Biodiversity Information
Facility (http://data.gbif.org) and the Ocean Biogeographic
Information System (http://www.iobis.org). We then used Max-
Ent software [60] to model fish species distributions across the
planning region for the 274 species (Table S1) with at least eight
non-duplicate records (Methods in Materials S1), a common
minimum threshold [61]. Environmental data in MaxEnt models
consisted of the following nine environmental datasets with an
original resolution between 1 km² and 1 degree: chlorophyll [33],
current velocity [33], photosynthetically active radiation [33], sea-
surface temperature [33], ultraviolet irradiance [33], wind speed
[33], salinity [62], depth [63], and percent reef [35].

In choosing species and habitats for the analysis, we selected those
that were most important biologically (e.g. rare, threatened),
important as providers of ecosystem services (e.g. keystone habitats),
and important culturally (e.g. for fisheries). Unfortunately, this
selection was necessarily limited by data availability. We could not
include maps of sea-grass distributions, for example, despite the
critical importance of this habitat for regional ecosystem services,
such as fisheries and nodes for connectivity [64].

We resampled all grids to 25 km² for use in MaxEnt and
subsequent analyses. Although it is possible to resample to a
smaller grid, we felt this is the minimum resolution that is likely to
reliably reflect patterns of biodiversity, exposure and resource use
across the study area given the current understanding and
resolution of the data.

We used the biodiversity features maps to measure biodiversity
value in two distinct ways. One, we calculated biological richness
as the sum total of fish species in each 25 km² cell. To determine
fish species presence, we first applied a threshold to the fish species
distribution models to convert them to binary presence-absence
maps. We used a fixed logistic threshold value of 0.75, which
conservatively selects environmental conditions where fish are
highly likely to be present. Biological richness is a common
measurement of biodiversity value that is widely used for

<figure><img src="image_9.png" /><figcaption>**Figure 3.** Mangroves, reef geomorphology and bioregions.
doi:10.1371/journal.pone.0028969.g003</figcaption></figure>

prioritizing conservation areas [65]. Two, we used the benefit-
function conservation priority-setting software Zonation [66] to
map biodiversity value of each cell (Methods in Materials S1).
Zonation uses principles of complementarity to produce contin-
uous ranking of biodiversity value from input species distribution
models, habitats and bioregions, so that selection of any subset of
top-ranked grid cells (e.g. the top quartile) maximizes species,
habitat, and bioregion representation within the subset.

**Mapping conservation and management priorities**

We conducted four separate analyses of the same primary
datasets (e.g. fishing pressure, exposure to thermal stress and

<figure><img src="image_10.png" /></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 4.** Additional biological inputs into the analysis. A: Biological richness showing number of fish species per 25 km² grid cell; B: Biodiversity value measured using Zonation algorithm. doi:10.1371/journal.pone.0028969.g004</figcaption></figure>

biodiversity) (Table 1): 1) Visual gradient overlay of these three primary datasets layers in red, green, and blue (RGB) color space; 2) categorical classification of proposed conservation and management action zones; 3) target-based site optimization algorithm using Marxan, the most widely used target-based marine conservation planning optimization method [24,67]; and finally, 4) conservation priority ranking with Zonation, an algorithm that produces a continuous measure of conservation value [68].

**Visual gradient overlay in RGB color space.** To map each of the three variables in continuous RGB space, we first converted the original data values to a 0–255 scale, so as to render them as 8-bit TIFF files in the ENVI image processing software [69]. Next, we displayed the three variables together as single bands in a 3-band RGB composite image, with biodiversity value assigned to the red band, fishing pressure assigned to green, and environmental exposure assigned to blue. We produced two maps, the first with biodiversity value equal to richness, and the second with biodiversity value mapped by Zonation. By default, ENVI applies a 2% linear stretch to images, which facilitates interpretation of color patterns. The value or intensity of each primary color indicates the relative strength of that variable in that location. In addition, the additive combinations of the primary colors show how the three variables interact across the study area. Areas dominated by blue, for example, show where exposure is the dominant variable; areas dominated by red show where biodiversity is stronger than fishing pressure or exposure, and so on.

**Categorical classification of proposed conservation and management action zones.** We built on conservation action frameworks applied recently to coral reefs and fishing landing sites in five countries [12,70] to classify the study region based on binary splits and combinations of the three variables. First, we converted the continuous 25 km² measures of biodiversity value measured by Zonation, fishing pressure and exposure to binary “low” and “high” values, defining “high” as the top quartile of values, and “low” as the bottom three quartiles for each input variable. Next, we combined the three binary grids to produce a single map with eight classes representing all unique combinations of the “high” and “low” values. Finally, we proposed a conservation and management action category for each class based on these unique attributes [12,70]. For example, in areas characterized by high biodiversity value, low fishing pressure and low exposure, we recommended a management regime that emphasizes large closures and strict protection, as these waters are likely to contain biodiversity of global significance, along with a high probability of persistence due to low climate impacts and conflicts with fishing pressure.

**Marxan target-based site optimization algorithm.** Marxan maps areas that meet quantitative representation targets for conservation features (in our case the distribution maps for fish, habitats, and bioregions), while minimizing an overall objective function score that can include different types of costs [24]. Here, we configured the objective function to penalize failure to meet feature targets, and minimize boundary length and total costs of

<figure><img src="image_12.png" /></figure>
<table><caption>**Table 2.** Biodiversity data used in the analysis.</caption><thead><tr><th>Data type</th><th>Dataset</th><th>Source</th></tr></thead><tbody><tr><td>Bioregions</td><td>3 West Coast Madagascar Bioregions</td><td>Indian Ocean Commission [18]</td></tr><tr><td>Habitats</td><td>Mangroves</td><td>Harper et al. [58]; Moat & Smith [59]</td></tr><tr><td></td><td>Reef geomorphology: Atoll rim</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Bank barrier</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Bank lagoon</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Coastal barrier reef complex</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Coastal/fringing patch</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Continental lagoon</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Diffuse fringing</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Fringing of coastal barrier complex</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Intra-lagoon patch-reef complex</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Intra-seas exposed fringing</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Intra-seas patch-reef complex</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Lagoon exposed fringing</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Ocean exposed fringing</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Outer barrier reef complex</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Shelf patch-reef complex</td><td>Andréfouët et al. [35]</td></tr><tr><td></td><td>Reef geomorphology: Shelf slope</td><td>Andréfouët et al. [35]</td></tr><tr><td>Species</td><td>251 fish species models generated using MaxEnt</td><td>REBIOMA, unpubl. data</td></tr></tbody></table>

doi:10.1371/journal.pone.0028969.t002

the selected areas (Methods in Materials S1). We selected a feature representation target of 30% for all conservation features. 30% is a common representation level for marine conservation planning [71], and matches the representation target recommended at a recent Indian Ocean Commission Protected Area Network (IOCPAN) conservation planning workshop [18]. Since Marxan is generally run with only binary (presence-absence) species distribution models as an input, for fish species, we used the previously described presence-absence species distribution models where thresholds were established. We used the fish catch data to assign an economic value or selection “cost” to each 25 km² planning unit. We also used a feature of Marxan that allows the inclusion of an additional probability-based cost for each planning unit, typically representing the likelihood of a threatening process, such as catastrophic bleaching [28]. We used the environmental exposure data [33] directly as this measure of future threat probability. By minimizing this set of terms, Marxan attempts to identify reserve systems that meet 30% targets for conservation features, minimize total boundary length and costs, and maximize the probability that conservation features will persist in the face of future threats. We ran Marxan 100 times, selecting the “best” run (i.e. the one with the lowest overall objective function score) as our result.

**Zonation conservation priority ranking.** The most recent Zonation release (version 3.0.2) allows input features to be positively or negatively weighted. This enables the inclusion of both positively weighted features representing biodiversity, as well as negatively weighted features to be avoided, such as exposure or other threatening conditions. For this analysis, we input continuous fish species distribution models, habitat maps, and bioregions as positively weighted features, and fishing pressure and exposure to thermal stress as negatively weighted features (Methods in Materials S1). In all other respects, this Zonation

analysis (hereafter “weighted Zonation”) was similar to the one described previously.

**Comparisons.** We compared the Marxan result with the strictest protection category of the categorical classification result (hereafter “large closures” result) and the weighted Zonation, as all three primarily represent areas of high biodiversity value, low cost, and low exposure probability. We made four comparisons of these three results. First, we calculated the total overlap between the three results. Second, we compared the average proportion of conservation features included in each result (Methods in Materials S1). Third, we compared total fish catch included in each result. This is the opportunity cost of withdrawing these areas from production, and serves as a measure of marine reserve efficiency [72]. Fourth, we compared the predicted persistence of each result in terms of exposure probability, in two ways: first by looking at exposure values per grid cell across each result, and second using the value reported by Marxan that shows the probability of a solution missing its conservation feature targets [73]. Although this is a Marxan-specific measure, we locked the other results into runs of Marxan to be able to compare this value (Methods in Materials S1). Finally, we ran Marxan a second time. Here, instead of setting feature representation targets to 30% for conservation features, we used a target equal to the average proportion of species and habitat distributions included in the large closures map, and again compared overlap, average proportion of species ranges included, opportunity cost, and exposure. Similarly, because Zonation produces a continuous ranking of biodiversity value, we selected top ranked cells of the weighted Zonation result until the average proportion of species and habitats represented matched that of the large closures map, and set this as a threshold to produce a binary Zonation result. This thresholded Zonation result (hereafter “Zonation 16%”) and the second Marxan analysis facilitate comparison of these results to

<figure><img src="image_13.png" /></figure>
the large closures map by ensuring that each achieves approximately the same representation targets.

## Results

We produced two maps of biodiversity value, one based on the estimated numbers of fish species (Figure 4a) and one with biodiversity value calculated by the first run of the Zonation algorithm on biodiversity features alone (Figure 4b). The fish species richness map shows that areas of high species richness are concentrated on the fringing reefs of the southwest, the barrier reefs off Maintarano and Besalampy, Juan de Nova, and the fringing reefs, islands and atolls of the Northwest (Figure 4a). The map produced by Zonation highlights many of these same areas for high biodiversity value (Figure 4b), with some key differences, notably higher estimates for the Banc de Leven some 40 km west-northwest of the Nosy Mitsio group, and the large shallow banks southwest of Cap Sainte Marie.

The visual overlay of the use, exposure, and biodiversity in RGB color space provides a visual interpretation of the relative strength of each of the three primary variables – use represented by green, biodiversity by red, and exposure by blue (Figures 5a,b). For example, in both results, large areas of the near-coastal southwest are yellow, indicating both high biodiversity and fishing, but low

exposure, while the central coastal areas are light-blue indicating high exposure and fishing, but low biodiversity value. The far west is dark blue where exposure is the dominant variable; and the large, near-coastal areas of the northwest trend towards white, indicating a saturation of all three variables. These two maps also facilitate the comparison of the two alternative biodiversity measures, by showing how they interact with the non-biological variables.

The result of the categorical management action classification extends this simple visual classification to assign a conservation management priority along a gradient of protection, use and exposure risk corresponding to the eight possible combinations of high and low biodiversity, fishing and exposure values (Figure 6a; Table 3). By far the largest area is within the sustainable use and maintain category (Class 5: 87,682 km²), which is made up of relatively deep, offshore habitats, where all three variables are low; hence conditions are suitable for managing use up to sustainable limits. This class includes 49.7% of the total study area, and highlights the fact that the majority of the study area (81.0%) maps to the lowest quartile of environmental exposure. In contrast, only 23.7% of the area maps to high biodiversity value (classes 1–4). All of the smaller size classes have high exposure, with a mix of biodiversity and fishing pressure. The smallest class 8 (412 km²), suggests providing relief to local communities, is orders of

<figure><img src="image_14.png" /></figure>

**Figure 5.** Two results of RGB visual overlay of primary variables (biodiversity, fishing pressure, exposure). A: Biodiversity value expressed as fish species richness; B: Biodiversity value measured using the Zonation algorithm. Key shows classification in 3-dimensional RGB color cube, with biodiversity (letter B in the key) assigned to Red (z-axis), fishing (F) assigned to Green (y-axis), and exposure (E) assigned to Blue (x-axis). Only the colors formed on the inner and outer planes of the cube are visible. On the inner planes, one variable is always 0. On the outer planes, one variable is always 255. The inner corner (black) has 0 values for all variables. The outer corner (white) has values of 255 for all variables.
doi:10.1371/journal.pone.0028969.g005

<figure><img src="image_15.png" /></figure>
<figure><img src="image_16.png" /><figcaption>**Figure 6.** Two views of conservation and management priorities. A: results of the categorical classification (see Table 3 for class descriptions); B: target-based priority-setting with Marxan. doi:10.1371/journal.pone.0028969.g006</figcaption></figure>

magnitude smaller than the largest, and is characterized by high exposure and fishing pressure, and low biodiversity value.

The Marxan 30% target result shows areas that meet quantitative conservation targets at minimum cost and high probability of persistence (Figure 6b). Results are distinctly clustered within the south, central and northern parts of the study area, with a scattering of individual planning units selected in between. By far the least selected areas are in the northwest between Cap d'Ambre and Analalava, where exposure probability and, therefore, the probability of persistence is low, and fishing pressure is also relatively high.

The weighted Zonation result shows all areas ranked in order of their conservation value (Figure 7). Because this run of Zonation tends to avoid negatively weighted features by giving them a low ranking (here fishing pressure and exposure), the highest ranked areas are generally those furthest away from where the negatively weighted variables have their strongest influence. As a consequence, the dominant concentrations of high ranked cells are found in three places: offshore southwest, offshore central in the vicinity of Juan de Nova, the Barren Islands and associated barrier reefs, and offshore northwest, on the Banc de Leven.

Comparing the management actions categorization and Marxan results, we find that the large closures result covers 20,528 km², and the Marxan 30% result covers 39,585 km². The large closures result shares 61.1% of its selected area with the Marxan result; 32.0% of the Marxan result overlaps with the large closures result. The Marxan result includes an average of 28.9% of

the total distribution of each species or habitat. The large closures result includes 16.3% of the total distribution of each species on average, which is significantly less than the Marxan result (two-sample paired t-test, t(590) = 18.4655, p<0.0001). Finally, assuming fishing is excluded entirely from the selected areas and there are no spillover benefits, the impact of the large closures categorization is 6,340.2 tons of fish/year, compared to 14,410.9 tons/year for Marxan. After reducing feature targets in Marxan to 16.3% (hereafter Marxan 16% result) to match the average proportion of the range captured by the large closures result, we find the Marxan 16% result covers a comparable area (21,266 km²), and there is no longer a significant difference between species proportions represented (Marxan 16% result: 15.4% of each species range, two-sample paired t-test, t(590) = -1.2875, p = 0.19) (Table 4). The large closures result shares 38.8% of its selected area with the Marxan 16% result, while 37.4% of Marxan 16% result overlaps with the large closures result. We find the fisheries impact of the Marxan 16% result is 6,393.6 tons/year, 56.0% less than included in the Marxan 30% target result (Table 4), and similar to the large closures result.

Comparing the Zonation 16% result to the other results, we find that this thresholded result covers 33,516 km², and has a predicted fisheries impact of 6,580.8 tons/year, and is therefore larger in area and fisheries impact than the other two results that represent about 16% of species' ranges on average, but smaller in both respects than the Marxan 30% result (Table 4). The mean proportion of ranges represented (16.3%) is not significantly

<figure><img src="image_17.png" /></figure>
<table><caption>**Table 3.** Management categories for management actions assigned to combinations of high and low variables based on the conceptual model of McClanahan et al. [68].</caption><thead><tr><th>Class</th><th>Area (km²)</th><th>Exposure</th><th>Fishing pressure</th><th>Biodiversity Value</th><th>Category</th><th>Description</th></tr></thead><tbody><tr><td>1</td><td>20,528</td><td>Low</td><td>Low</td><td>High</td><td>Large closures</td><td>Maintain high biodiversity values through strategies including large closures, strict protection and minimal use where possible</td></tr><tr><td>2</td><td>14,622</td><td>Low</td><td>High</td><td>High</td><td>Protect and restore</td><td>Maintain and restore high biodiversity values through closures, strict protection, reduced use and active restoration</td></tr><tr><td>3</td><td>5,314</td><td>High</td><td>Low</td><td>High</td><td>Protect and recover</td><td>Protection, recovery and restoration necessary to maintain high biodiversity values with a degree of resignation to ecological degradation: management involves risk of failure due to high exposure</td></tr><tr><td>4</td><td>1,357</td><td>High</td><td>High</td><td>High</td><td>Protect and manage</td><td>Protection, recovery and restoration necessary to maintain high biodiversity values. High risk of failure due to high exposure and use</td></tr><tr><td>5</td><td>87,682</td><td>Low</td><td>Low</td><td>Low</td><td>Use sustainably and maintain</td><td>Manage towards limits of sustainable use. Lower biodiversity values, but lower risk due to lower exposure and use; overuse to extinction unlikely</td></tr><tr><td>6</td><td>19,426</td><td>Low</td><td>High</td><td>Low</td><td>Use sustainably and recover</td><td>Manage towards limits of sustainable use with recovery and restoration, reduce use towards the limits of sustainable yield</td></tr><tr><td>7</td><td>27,018</td><td>High</td><td>Low</td><td>Low</td><td>Use sustainably and manage</td><td>Manage towards limits with some resignation to ecological degradation; management interventions carry risk of failure due to high exposure, and low returns due to low biodiversity value.</td></tr><tr><td>8</td><td>412</td><td>High</td><td>High</td><td>Low</td><td>Provide relief</td><td>Communities here may require relief because of high existing fishing pressure and uncertain future resilience and sustainability</td></tr></tbody></table>

Total area sums to less than the study area (86%) because of the extension of the study area by 0.25 degrees inland.
doi:10.1371/journal.pone.0028969.t002

different from the other results (Marxan 16% result: mean 15.4% of each species range, two-sample paired t-test, t(590) = -1.2582, p = 0.2088; Large closures result: mean 16.3% of each species range, t(590) = -6e-04, p = .9995). The Zonation 16% result shares 21.4% and 31.7% of its selected area with Marxan 16% and the large closures result, respectively. The Marxan 16% result and the large closures result, on the other hand, share 33.7% and 51.8% of their selected area with the Zonation 16% result (Table 4, Figures 8–9).

Finally, we find that the distribution of exposure values in the large closures result is significantly lower than in both Marxan results (Wilcoxon rank-sum test, p<0.0001, two-tailed). The Zonation 16% result has the lowest mean exposure values, and this distribution is also significantly different from all other results (Wilcoxon rank-sum test, p<0.0001, two-tailed). The probability of missing conservation feature targets due to exposure, however, is lower in both Marxan results and highest in the Zonation 16%

result (Zonation 16%: 0.29, Large closures: 0.26, Marxan 16%: 0.14, Marxan 30%: 0.07, Table 4).

Discussion

We had two primary aims with this work. First, we applied several different methods for mapping conservation management options across a relatively vast, poorly known region, to see if broad trends emerge irrespective of approach. Second, we evaluated the conditions where alternative approaches might provide results comparable to methods using optimization. Overall, despite the different methods employed, the results show substantial similarities (Figures 5, 6, 7, 8, 9). For example, many areas dominated by red and purple in the continuous RBG color space results (areas of high biodiversity value) were also identified in the Marxan, large closures, and weighted Zonation results. The most prominent examples are the Barren Islands and associated

<figure><img src="image_18.png" /></figure>
<figure><img src="image_19.png" /><figcaption>**Figure 7. Weighted Zonation result.** This map shows a continuous ranking of conservation value by the Zonation algorithm. Higher ranked cells are more important for species representation, and tend to have lower fishing pressure and exposure values.
doi:10.1371/journal.pone.0028969.g007</figcaption></figure>

barrier reef systems north and east of Maintirano, and the Juan de
Nova Island reefs. Two additional large areas, the Banc de Leven
to the northwest, and the large shallow banks southwest of Cap
Sainte Marie, are highlighted as conservation priorities in all
results with the exception of the RGB map, which used species
richness as the measure of biodiversity value (Fig. 4a). Species
richness is the only biodiversity measure used here that does not
incorporate notions of complementarity, which likely explains why
these areas do not turn up in this one result. While the Barren

Islands have received conservation attention in recent years, to our
knowledge the other areas have not, and this convergence in
prioritization among methods warrants additional field investiga-
tions.

Comparing results: threats, opportunities and benefits

Each approach comes with a distinct set of strengths and
weaknesses. Although not strictly a planning approach per se, the
RGB maps are valuable for providing a rapid visual overview of
the distribution of input variables and their interactions. The
primary intent of this technique was descriptive, not prescriptive.
Therefore, translating these simple RGB maps into specific
management actions is challenging. The Marxan result, weighted
Zonation result, and the categorical management action classifi-
cation, on the other hand, can be used to recommend specific
management and conservation actions, but based on different
assumptions. As configured here, the Marxan solution shows areas
that meet species representation targets set by analysts, maximize
persistence in the face of environmental exposure, and reduce
impacts on fisheries by minimizing the opportunity costs from
foregone fishing revenues. The weighted Zonation result also
shows areas that maximize species representation while avoiding
areas of high exposure and fishing pressure. Unlike Marxan,
however, Zonation produces a continuous ranking of biodiversity.
In Zonation, therefore, it is up to users to decide which rank, or
solution area, to choose as a threshold for a final result. In
Marxan, by contrast, it is up to users to set targets for
representation. Typically, Marxan and Zonation results such as
these are used to guide decisions concerning the placement of
strictly protected reserves [25,74]. Because both the Marxan and
the weighted Zonation results provide little information on the
value, management and use of areas outside of those selected as
top conservation priorities, however, they have less utility as
comprehensive zoning tools.

The value of the categorical management action classification is that it can guide the placement of strict marine reserves and also zone the entire planning region into a variety of management categories that include various actions aimed at achieving more sustainable use of the seascape [70]. Zoning the entire area allows flexibility to consider management of “suboptimal” areas that would not be selected under optimization approaches, including degraded areas where effective management may lead to recovery. For example, maintaining herbivores through reduced fishing may help reefs recover from damage due to climate disturbances or destructive fishing [75]. The disadvantage of this method relative to optimization methods, however, is that they are likely to pick regions for strict conservation less efficiently, by not meeting biodiversity representation targets at minimum cost [76].

Fishing is a factor that has often undermined the success of protected area management [77,78]. Despite this, it is frequently treated only as a factor to avoid in marine conservation planning. For example, Marxan is routinely configured to minimize fisheries conflicts and opportunity costs by avoiding areas of high fishing value [23]. However, there can be equally valid reasons for prioritizing protection of heavily fished areas because they may be areas of high productivity and source populations [79] that could stimulate spillover to fisheries [80,81], and because closures can have their greatest fisheries benefits when fishing is beyond a maximum sustained yield [7,82,83]. Reserve networks that exclude fishing may not provide the economic benefits necessary for social adoption nor improve management more than standard fisheries management tools [84].

Many other ecosystem goods and services can be threatened by
unsustainable management, and improved spatial planning tools

<figure><img src="image_20.png" /></figure>
<table><caption>**Table 4.** Comparison of results.</caption><thead><tr><th>Result</th><th>Area (km²)</th><th>Average total distribution represented</th><th>Fishing cost (tons/year)</th><th>Mean exposure probability</th><th>Probability result misses targets</th></tr></thead><tbody><tr><td>Categorical classification "large closures"</td><td>20,528</td><td>16.3%</td><td>6,340.2</td><td>0.47</td><td>0.26</td></tr><tr><td>Marxan with 30 percent targets</td><td>39,585</td><td>28.9%</td><td>14,410.9</td><td>0.52</td><td>0.07</td></tr><tr><td>Marxan with 16 percent targets</td><td>21,266</td><td>15.4%</td><td>6,393.6</td><td>0.52</td><td>0.14</td></tr><tr><td>Zonation top ranked cells</td><td>33,516</td><td>16.3%</td><td>6,580.8</td><td>0.38</td><td>0.29</td></tr></tbody></table>

doi:10.1371/journal.pone.0028969.t003

are needed to handle this complexity. Emerging software, such as Marxan with Zones [67] and Marine InVEST [85], are promising. We did not compare either of these programs to the categorical management action classification because of our primary interest in evaluating Marxan and Zonation, the most widely used conservation planning tools in the region. Future efforts could consider a broader array of available and emerging approaches, however.

## Comparing results: impacts on fisheries and persistence

Marine spatial plans are frequently evaluated for their efficiency in meeting conservation objectives while minimizing economic impacts [72]. Given the assumption of no spillover and that strict protection completely displaces fishing activities, the impact of the large closures result on fisheries was ~55% less than the Marxan 30% result. This difference in economic impact is explained by differences in the area of the two results. The large closures result is half the size of the Marxan area and represents significantly less of the total distribution of each species, on average. Importantly, there is little difference in economic impact between methods when Marxan targets are reduced to match the species proportions covered by the large closures result (Table 4). The Zonation 16% result, on the other hand, is about 3–4% more costly in terms of fisheries impacts than the other two results that also achieve about 16% representation. Giving additional negative weight to this factor in Zonation could potentially reduce this difference.

Spatial overlap between the three 16% representation results is relatively low (Figure 9). Despite similarities in cost and species proportions represented, the three results only share 9% of their total selected area. Two factors may explain these differences. One, according to our models, species beta-diversity is relatively low across the study area. When beta-diversity is low, many areas have similar species compositions, and as a consequence, many different potential reserve configurations can achieve roughly equal representation. Two, the method Marxan uses to incorporate the exposure variable is quite different from the other approaches. Marxan minimizes exposure probabilities across conservation features (species, habitats, bioregions), in balance with other costs and penalties, whereas the large closure and weighted Zonation results tend to avoid areas of high exposure entirely. Comparing raw exposure values, the Zonation 16% result consequently has the lowest mean exposure, followed by the large closure result, then the two Marxan results (Table 4). The exposure distributions from both the Zonation 16% result and the large closure result are significantly different from each other and the Marxan results. On the surface, the fact that exposure probability is lower in the large closures and Zonation 16% results than the Marxan optimizations is surprising. Marxan, however, seeks solutions that minimize exposure probability in order to meet representation targets for each individual conservation feature [28]. As a result, the probability of the large closures result missing

feature targets is higher than both Marxan results, and highest in the Zonation 16% result (Table 4).

## Targets

By reducing Marxan representation targets to 16.3% to match the target level implied by the large closures result, we show that a classification-based method performs comparably to Marxan in terms of species representation and efficiency, though not as well in terms of the persistence of conservation features (Table 4). Similarly, we show that a 16% representation threshold on the weighted Zonation result produces comparable results to the other methods. Although necessary to compare approaches, this neglects an important question: what are adequate representation targets for marine spatial planning in Madagascar? While iterative, target-based approaches ensure species representation at the specified level [86], little is known about what constitutes an adequate target area for marine organisms [87]. For example, The Convention on Biological Diversity recommends “at least 10% of each of the world’s marine and coastal ecological regions be effectively conserved by 2010” [88], while other organizations call for two to three times this level [89]. Targets used here fall within the broad range of these recommendations, but neither are adequately informed by species requirements. For example, many large migratory species such as sea turtles and cetaceans have specific breeding and nesting requirements. Protecting a portion of their local distribution may fail to protect these species over the long term unless augmented by ecologically meaningful guidelines across their distributional range [90]. Consequently, future national level marine conservation planning in Madagascar will benefit from a more rigorous elaboration of conservation targets, based on a comprehensive review and understanding of individual species life history traits, area requirements, and sources of mortality.

## Improvements and next steps

We expect marine spatial planning in Madagascar to improve with the availability of additional biological and environmental data. For example, datasets currently unavailable at sufficient resolution across the planning region include data on seagrass communities, information on habitat condition and intactness, site connectivity, and potential climate-induced acidification. Further, predicted responses, willingness to adopt management recommendations, and adaptive capacity of people to climate change is needed [70].

The Government of Madagascar committed to a plan to increase MPA coverage by at least one million hectares in 2003 [91], and this work is expected to continue into 2012 despite political changes resulting from political unrest in 2009. Consequently, the work presented here provides a foundation for the national level analyses needed to fulfill that target. There is considerable overlap in the areas identified for the highest priority

<figure><img src="image_21.png" /></figure>
<figure><img src="image_22.png" /><figcaption>**Figure 8.** Comparison of results. A: Strict protection class of the categorical classification; B: Marxan 30% solution; C: Marxan 16% solution; D: Zonation 16% solution.
doi:10.1371/journal.pone.0028969.g008</figcaption></figure>

reserves and, as the national discussion and management evolves, these methods and outputs can further help guide these decisions. Because implementation is both a political and a scientific process, however [92,93], ultimately, the choice of approach should consider both goals for representation and persistence as well as the likelihood of adoption and compliance [78,94,95,96], and the timing of implementation [97]. With this study, we show that methods based on Zonation and a simple combination of variables can produce result for strict protection that are similarly representative and have similar economic impacts as ones based on optimization. More broadly, we demonstrate the utility of

<figure><img src="image_23.png" /><figcaption>**Figure 9.** Overlap between results. The number (one, two or three) indicates the number of 16% solutions represented; in other words, the number of times a planning unit has been selected by either the strict protection, Marxan 16% or Zonation 16% result.
doi:10.1371/journal.pone.0028969.g009</figcaption></figure>

comparing alternative methods early in the planning process for understanding patterns and interactions of key biodiversity and conservation variables.

**Supporting Information**

**Materials S1** This text provides additional information on methods including modeling fish species distributions with MaxEnt, details of Marxan and Zonation

<figure><img src="image_24.png" /></figure>
## runs, measuring representation of conservation features across results, and measuring exposure.
(DOC)

**Table S1 Focal fish species list.** This list shows the 274 fish species with distributions modeled in MaxEnt for use in the analysis and data sources. OBIS refers to the Ocean Biogeo- graphic Information System (http://www.iobis.org/). WCS, WWF, and CI refer to data provided by Wildlife Conservation Society, World Wildlife Fund, and Conservation International, respectively. COOKE refers to data provided by Andrew Cooke, and SAMOILY refers to data provided by Dr. Melita Samoily. Species marked with an asterix (*) were not included in the Marxan analyses because their continuous distribution value was below threshold.

## Acknowledgments

We thank the John D. and Catherine T. MacArthur Foundation for supporting this work. The fish catch model was developed from the “Sea

Around Us” global catch database. This database was developed through a scientific collaboration between the University of British Columbia and the Pew Environment Group. The Indian Ocean Commission, Fond Français pour l'Environnement Mondial (FFEM) and World Wildlife Fund supported the Indian Ocean Commission Protected Areas Network conservation planning workshop. We greatly appreciate data contributions from A. Cooke and M. Samoily. E. Lagabrielle acknowledges the African Coelacanth Ecosystem Programme and the National Research Foundation of South Africa for co-funding its contribution to this research. Finally, we thank the editor, reviewers, and G. Powell and J. Watson for constructive comments on earlier drafts of the paper.

## Author Contributions

Conceived and designed the experiments: TFA TRM CM AR TT CK. Performed the experiments: TFA AR TT. Analyzed the data: TFA TRM AR TT CK SA. Contributed reagents/materials/analysis tools: AR TT EL RW. Wrote the paper: TFA TRM CM AR TT CK RW EL MB SA. Provided reef geomorphology dataset: SA.

## References

1. Worm B, Barbier EB, Beaumont N, Duffy JE, Folke C, et al. (2006) Impacts of biodiversity loss on ocean ecosystem services. Science 314: 787-790.

2. Jackson JBC (2008) Ecological extinction and evolution in the brave new ocean. Proceedings of the National Academy of Sciences 105: 11458-11465.

3. Botsford LW, Castilla JC, Peterson CH (1997) The Management of fisheries and marine ecosystems. Science 277: 509-515.

4. Worm B, Hilborn R, Baum JK, Branch TA, Collie JS, et al. (2009) Rebuilding global fisheries. Science 325: 578-585.

5. McClanahan TR, Graham NAJ, Calnan JM, MacNeil MA (2007) Toward pristine biomass: reef fish recovery in coral reef marine protected areas in Kenya. Ecological Applications 17: 1055-1067.

6. Lester SE, Halpern BS, Grorud-Colvert K, Lubchenco J, Ruttenberg BI, et al. (2009) Biological effects within no-take marine reserves: a global synthesis. Marine Ecology Progress Series 384: 33-46.

7. McClanahan TR (2010) Effects of fisheries closures and gear restrictions on fishing income in a Kenyan coral reef. Conservation Biology 24: 1519-1528.

8. Knowlton N, Jackson JBC (2008) Shifting baselines, local impacts, and global change on coral reefs. PLoS Biology 6.

9. Leslie HM (2005) A synthesis of marine conservation planning approaches. Conservation Biology 19: 1701-1713.

10. Olson DM, Dinerstein E (1998) The Global 200: A representation approach to conserving the Earth's most biologically valuable ecoregions. Conservation Biology 12: 502-515.

11. Wafar M, Venkataraman K, Ingole B, Ajmal Khan S, LokaBharathi P (2011) State of knowledge of coastal and marine biodiversity of Indian Ocean countries. PLoS ONE 6: e14613.

12. McClanahan TR, Cinner JE, Maina J, Graham NAJ, Daw TM, et al. (2008) Conservation action in a changing climate. Conservation Letters 1: 53-59.

13. World Bank (2007) Commune-level census data for Madagascar. World Bank, FID-Madagascar.

14. INSTAT, ORC Macro (2005) Enquête démographique et de santé, Madagascar 2003-2004: Rapport de synthèse. Calverton, Maryland, USA. 34 p.

15. Laroche J, Razanoelisoa J, Fauroux E, Rabenevanana MW (1997) The reef fisheries surrounding the south-west coastal cities of Madagascar. Fisheries Management and Ecology 4: 285-299.

16. Rakotoson LR, Tanner K (2006) Community-based governance of coastal zone and marine resources in Madagascar. Ocean and Coastal Management 49: 855-872.

17. Harris A, Manahira G, Sheppard A, Gough C, Sheppard C (2010) Demise of Madagascar's once great barrier reef - change in coral reef condition over 40 years. Atoll Research Bulletin 574: 1-18.

18. Indian Ocean Commission (2010) A regional strategy for conserving marine ecosystems and fisheries of the Western Indian Ocean Islands Marine Ecoregion (WIOMER). Draft Report. Quatre Bornes, Mauritius and Antananarivo, Madagascar: Indian Ocean Commission and World Wildlife Fund 2010.

19. Foley MM, Halpern BS, Micheli F, Armsby MH, Caldwell MR, et al. (2010) Guiding ecological principles for marine spatial planning. Marine Policy 34: 955-966.

20. Douvere F (2008) The importance of marine spatial planning in advancing ecosystem-based sea use management. Marine Policy 32: 762-771.

21. Lourie SA, Vincent ACJ (2004) Using biogeography to help set priorities in marine conservation. Conservation Biology 18: 1004-1020.

22. Scott JM, David F, Csuti B, Noss R, Butterfield B, et al. (1993) Gap analysis: a geographic approach to protection of biological diversity. Wildlife Monographs 123: 1-41.

23. Sala E, Aburto-Oropeza O, Paredes G, Parra I, Barrera JC, et al. (2002) A general model for designing networks of marine reserves. Science 298: 1991-1993.

24. Ball IR, Possingham HP, Watts ME (2009) Marxan and relatives: software for spatial conservation prioritization. In: Moilanen A, Possingham HP, Wilson KA, eds. Spatial conservation prioritization: Quantitative methods and computational tools. Oxford: Oxford University Press. pp 185-195.

25. Leathwick J, Moilanen A, Francis M, Elith J, Taylor P, et al. (2008) Novel methods for the design and evaluation of marine protected areas in offshore waters. Conservation Letters 1: 91-102.

26. Ban NC, Klein CJ (2009) Spatial socioeconomic data as a cost in systematic marine conservation planning. Conservation Letters 2: 206-215.

27. Weeks R, Russ GR, Bucol AA, Alcala AC (2010) Shortcuts for marine conservation planning: The effectiveness of socioeconomic data surrogates. Biological Conservation 143: 1236-1244.

28. Game ET, Watts ME, Wooldridge S, Possingham HP (2008) Planning for persistence in marine reserves: A question of catastrophic importance. Ecological Applications 18: 670-680.

29. McLeod E, Salm R, Green A, Almany J (2009) Designing marine protected area networks to address the impacts of climate change. Frontiers in Ecology and the Environment 7: 362-370.

30. Beck MW (2003) The sea around: Conservation planning in marine regions. In: Groves C, ed. Drafting a conservation blueprint: A practitioner's guide to planning biodiversity. Washington, DC: Island Press. pp 319-344.

31. Klein C, Steinback C, Scholz A, Possingham H (2008) Effectiveness of marine reserve networks in representing biodiversity and minimizing impact to fishermen: a comparison of two approaches used in California. Conservation Letters 1: 44-51.

32. Wessel P, Smith WHF (1996) A global self-consistent, hierarchical, high-resolution shoreline database. Journal of Geophysical Research 101: 8741-8743.

33. Maina J, Venus V, McClanahan TR, Ateweberhan M (2008) Modelling susceptibility of coral reefs to environmental stress using remote sensing data and GIS models. Ecological Modelling 212: 180-199.

34. McClanahan TR, Ateweberhan M, Omukoto J, Pearson L (2009) Recent seawater temperature histories, status, and predictions for Madagascar's coral reefs. Marine Ecology Progress Series 380: 117-128.

35. Andréfouët S, Chagnaud N, Kranenburg CJ (2009) Atlas of Western Indian Ocean coral reefs. Centre IRD de Nouméa, New-Caledonia. 157 p.

36. Veron JEN, Turak E (2005) Zooxanthellate scleractinia of Madagascar. In: McKenna SA, Allen GR, eds. A rapid marine biodiversity assessment of the coral reefs of northwest Madagascar. Washington, DC: Conservation International. pp 23-25.

37. Cooke AJ, Lutjeharms J, Vasseur P (2003) Marine and coastal ecosystems of Madagascar. In: Goodman S, Benstead J, eds. Natural history of Madagascar. Chicago: Chicago University Press. pp 179-208.

38. McKenna SA, Allen GR, eds. A rapid marine biodiversity assessment of the coral reefs of Northwest Madagascar. Washington, DC: Conservation International.

39. Rakotonirina B, Cooke A (1994) Sea turtles of Madagascar - their status, exploitation and conservation. Oryx 28: 51-61.

<figure><img src="image_25.png" /></figure>
40. Kiszka J, Muir C, Poonian C, Cox TM, Amir OA, et al. (2009) Marine mammal bycatch in the Southwest Indian Ocean: review and need for a comprehensive status assessment. Western Indian Ocean Journal of Marine Sciences 7: 119-136.

41. Rosenbaum HC (2003) Marine mammals of Madagascar. In: Goodman S, Benstead J, eds. Natural history of Madagascar. Chicago, Illinois. pp 217-221.

42. Scharl M, Hornung U, Hismann K, Schauer J, Fricke H (2005) Genetics: Relatedness among east African coelacanths. Nature 435: 901-901.

43. Fishpool LDC, Evans MI (2001) Important bird areas in Africa and associated islands. Priority sites for conservation. Cambridge, UK: Pisces Publications and Birdlife International.

44. Le Corre M, Bemanaja E (2009) Discovery of two major seabird colonies in Madagascar. Marine Ornithology 37: 153-158.

45. Sheppard CRC (2000) Coral Reefs of the Western Indian Ocean: An Overview. In: McClanahan T, Sheppard CRC, Obura DO, eds. Coral reefs of the Indian Ocean. New York: Oxford University Press. pp 3-34.

46. Gabrié C, Vasseur P, Randriamiarana H, Maharavo J, Mara E (2000) The coral reefs of Madagascar. In: McClanahan T, Sheppard CRC, Obura DO, eds. Coral reefs of the Indian Ocean. New York: Oxford University Press. pp 411-444.

47. McClanahan TR, Maina JM, Muthiga NA (2011) Associations between climate stress and coral reef diversity in the Western Indian Ocean. Global Change Biology. pp no-no.

48. Davies TE, Beanjara N, Tregenza T (2009) A socio-economic perspective on gear-based management in an artisanal fishery in south-west Madagascar. Fisheries Management and Ecology.

49. Anderson SC, Flemming JM, Watson R, Lotze HK (In press) Serial exploitation of global sea cucumber fisheries. Fish and Fisheries.

50. Barnes DKA, Rawlinson KA (2009) Traditional coastal invertebrate fisheries in south-western Madagascar. Journal of the Marine Biological Association of the United Kingdom 89: 1589-1596.

51. McVean AR, Walker RCJ, Fanning E (2006) The traditional shark fisheries of southwest Madagascar: A study in the Toliara digital Fisheries Research 82: 280-289.

52. SAPM (2010) Protected areas of Madagascar digital GIS shapefiles. Système d'Aires Protégées de Madagascar, Antananarivo, Madagascar.

53. Ministre de L'Environnement et des Forêts, Ministre de L'Élevage, Ministre de L'Agriculture, Ministre de La Pêche et des Resources Halieutiques, Ministre des Mines et Des Hydrocarbures, et al. (2010) Arrêté interministériel n°52005/2010 du 20 décembre 2010 modifiant l'arrêté interministériel n°18633.

54. Lagabrielle E, Maina J, Massé L, Séré M, Watson R (2009) Mapping anthropogenic drivers of change for biodiversity conservation planning in the Western Indian Ocean. AntananarivoMadagascar: Indian Ocean Commission & World Wildlife Fund Madagascar. 53 p.

55. Watson R, Revenga C, Kura Y (2006) Fishing gear associated with global marine catches: I. Database development. Fisheries Research 79: 97-102.

56. Watson R, Revenga C, Kura Y (2006) Fishing gear associated with global marine catches: II. Trends in trawling and dredging. Fisheries Research 79: 103-111.

57. Pressey RL, Cowling RM, Rouget M (2003) Formulating conservation targets for biodiversity pattern and process in the Cape Floristic Region, South Africa. Biological Conservation 112: 99-127.

58. Harper GJ, Steininger MK, Tucker CJ, Juhn D, Hawkins F (2007) Fifty years of deforestation and forest fragmentation in Madagascar. Environmental Conservation 34: 325-333.

59. Moat J, Smith P (2007) Atlas of the vegetation of Madagascar. Kew: Kew Publishing, Royal Botanic Gardens.

60. Phillips SJ, Anderson RP, Schapire RE (2006) Maximum entropy modeling of species geographic distributions. Ecological Modelling 190: 231-259.

61. Pearson RG, Raxworthy CJ, Nakamura M, Townsend Peterson A (2007) Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar. Journal of Biogeography 34: 102-117.

62. Boyer TP, Antonov JI, Garcia HE, Johnson DR, Locarnini RA, et al. (2006) World Ocean Database 2005. Washington, D.C.: U.S. Government Printing Office.

63. Becker JJ, Sandwell DT, Smith WHF, Braud J, Binder B, Depner J, Fabre D, Factor J, Ingalls S, Kim S-H, Ladner R, Marks K, Nelson S, Pharaoh A, Trimmer R, Von Rosenberg J, Wallace G, Weatherall P (2009) Global Bathymetry and Elevation Data at 30 Arc Seconds Resolution: SRTM30_PLUS. Marine Geodesy 32: 355-371.

64. Unsworth RKF, Cullen LC (2010) Recognising the necessity for Indo-Pacific seagrass conservation. Conservation Letters 3: 63-73.

65. Fleishman E, Noss RF, Noon BR (2006) Utility and limitations of species richness metrics for conservation planning. Ecological Indicators 6: 543-553.

66. Moilanen A (2007) Landscape Zonation, benefit functions and target-based planning: Unifying reserve selection strategies. Biological Conservation 134.

67. Watts ME, Ball IR, Stewart RS, Klein CJ, Wilson K, et al. (2009) Marxan with Zones: Software for optimal conservation based land- and sea-use zoning. Environmental Modeling and Software 24: 1513-1521.

68. Moilanen A, Anderson BJ, Eigenbrod F, Heinemeyer A, Roy DB, et al. (2011) Balancing alternative land uses in conservation prioritization. Ecological Applications 21: 1419-1426.

69. Research Systems Incorporated (2004) Environment for Visualizing Images (ENVI) User's Guide. 4.7 ed. Boulder, CO: Research Systems Inc.

70. McClanahan TR, Cinner JE, Graham NAJ, Daw TM, Maina J, et al. (2009) Identifying reefs of hope and hopeful actions: Contextualizing environmental, ecological, and social parameters to respond effectively to climate change. Conservation Biology 23: 662-671.

71. Leslie H, Ruckelshaus R, Ball IR, Andelman S, Possingham HP (2003) Using siting algorithms in the design of marine reserve networks. Ecological Applications 13: S185-S198.

72. Stewart RR, Possingham HP (2005) Efficiency, costs and trade-offs in marine reserve system design. Environmental Modeling and Assessment 10: 203-213.

73. Watts M, Kircher L (2009) Advanced Marxan course manual. BrisbaneQueensland: The Ecology Centre, University of Queensland. pp 39.

74. Airame S, Dugan JE, Lafferty KD, Leslie H, McArdle DA, et al. (2003) Applying ecological criteria to marine reserve design: A case study from the California Channel Islands. Ecological Applications 13: 170-184.

75. Game ET, Bode M, McDonald-Madden E, Grantham HS, Possingham HP (2009) Dynamic marine protected areas can improve the resilience of coral reef systems. Ecology Letters 12: 1336-1345.

76. Pressey RL, Nicholls AO (1989) Efficiency in conservation evaluation: Scoring versus iterative approaches. Biological Conservation 50: 199-218.

77. McClanahan TR (1999) Is there a future for coral reef parks in poor tropical countries? Coral Reefs 18: 321-325.

78. McClanahan TR, Marnane MJ, Cinner JE, Kiene W (2006) A comparison of Marine Protected Areas and alternative approaches to coral reef management. Current Biology 16: 1408-1413.

79. Gaines SD, White C, Carr MH, Palumbi SR (2010) Designing marine reserve networks for both conservation and fisheries management. Proceedings of the National Academy of Sciences 107: 18286-18293.

80. Crowder LB, Lyman SJ, Figueira WF, Priddy J (2000) Source-sink population dynamics and the problems of siting marine reserves. Bulletin of Marine Science. pp 799-820.

81. Kellner JBT, Gaines I, Nisbet SD, R M (2007) Fishing the line near marine reserves in single and multispecies fisheries. Ecological Applications 17: 1039-1054.

82. Nowlis JS, Roberts CM (1999) Fisheries benefits and optimal design of marine reserves. Fishery Bulletin 97: 604-616.

83. Rodwell LD, Barbier EB, Roberts CM, McClanahan TR (2002) A model of tropical marine reserve-fisheries linkages. Natural Resource Modeling 15: 453-486.

84. Hilborn R, Stokes K, Maguire J-J, Smith T, Botsford LW, et al. (2004) When can marine reserves improve fisheries management? Ocean & Coastal Management 47: 197-205.

85. Ruckelshaus M, Guerry AD (2009) Valuing marine ecosystems? The Marine Scientist 26: 26-29.

86. Kirkpatrick JB (1983) An iterative method for establishing priorities for the selection of nature reserves: An example from Tasmania. Biological Conservation 25: 127-134.

87. Wabnitz CCC, Andréfouët S, Muller-Karger FE (2010) Measuring progress toward global marine conservation targets. Frontiers in Ecology and the Environment 8: 124-129.

88. Convention on Biological Diversity (2006) COP 8 - Eighth Ordinary Meeting of the Conference of the Parties to the Convention on Biological Diversity. Curitiba, 20-31 March 2006. Decision VIII/15: Framework for monitoring implementation of the achievement of the 2010 target and integration of targets into the thematic programmes of work. Annex 4. United Nations Environment Programme.

89. International Union for Conservation of Nature (2005) Benefits beyond boundaries. Proceedings of the Vth IUCN World Parks Congress. Gland, Switzerland and Cambridge, UK: IUCN.

90. Linquist S (2008) But is it progress? On the alleged advances of conservation biology over ecology. Biology and Philosophy 23: 529-544.

91. Norris S (2006) Madagascar defiant. BioScience 56: 960-965.

92. Knight AT, Driver A, Cowling RM, Maze K, Desmet PG, et al. (2006) Designing systematic conservation assessments that promote effective implementation: best practice from South Africa. Conservation Biology 20: 739-750.

93. Lagabrielle E, Botta A, Daré W, David D, Aubert S, et al. (2010) Modelling with stakeholders to integrate biodiversity into land-use planning - Lessons learned in Réunion Island (Western Indian Ocean). Environmental Modelling Software 25: 1413-1427.

94. Christie P (2004) Marine protected areas as biological successes and social failures in southeast Asia. American Fisheries Society Symposium 42: 155-164.

95. Johannes RE (2002) The renaissance of community-based marine resource management in Oceania. Annual Review of Ecology and Systematics 33: 317-340.

96. McClanahan TR, Mwaguni S, Muthiga NA (2005) Management of the Kenyan coast. Ocean & Coastal Management 48: 901-931.

97. Meir E, Andelman S, Possingham HP (2004) Does conservation planning matter in a dynamic and uncertain world? Ecology Letters 7: 615-622.

<figure><img src="image_26.png" /></figure>