All priority areas under the three targets in scenario ‘a’ (PAs must always be included) are more compact with less edge and have smaller mean boundary lengths (Table S3), whereas all scenarios where PAs are not obligatory to include (scenario 'b') have the lowest mean cost (Table S3). The total size of the priority areas is larger under scenario 'a,' but when only those from outside the current PA are considered, the solutions from 'b' variants are always larger in size (Table 1). For example, scenario 17% 'a' has a much larger priority area (20,924 km2) than scenario 'b' (4,906 km2), but when we considered only the areas outside the current PA network, scenario ‘b’ has a larger priority area (3,368 km2) than scenario ‘a’ (2,322 km2) (Table 1; Fig. 3).
Table 1
Size of the priority area (for planning units with a selection frequency >75) in scenarios ‘a’ and ‘b’.
Scenario target
|
Total priority area including PA (km2)
|
Priority area outside current PA (km2)
|
a
|
b
|
a
|
b
|
17%
|
20,924
|
4,906
|
2,322
|
3,368
|
30%
|
31,582
|
27,771
|
12,980
|
21,296
|
50%
|
67,338
|
63,685
|
48,736
|
51,218
|
Under scenario ‘b’, the size of the priority area within those already designated as a PA is always smaller for all three targets. This is demonstrated by scenario 17% ‘b’ which met this target by selecting a significantly smaller area within the PAs (1,538 km2) as a priority, as opposed to scenario 17% ‘a’, which select an area significantly larger (18,602 km2) out of the current PA network of 18,822 km2 (Fig. 3). Because of the small size of the total priority area within PA in scenario ‘b’ several of the established PAs were only partially or not at all selected.
The network of 79 PAs in Sulawesi includes 16.3% of the entire forest area coverage in 2018 (14,475 km2) and 8.9% of all karst ecosystems (1,681 km2), from the total available in the study area (Table S4). Of the island's 46,156 km2 area with the highest quarter (top quantile) of soil carbon organic, only 17.6% (8,109 km2) is protected (located within current PA). The remaining area with top quantile soil carbon organic is not protected and is in clusters 2 and 3. Within these PAs, the sub-montane forest has the most coverage (20.1%; 3,805 km2), while the low elevation hills forest has the least (11.7%; 2,945 km2) (Table S4).
As demonstrated by the representativeness within the priority area, larger does not always imply better. Apart from scenarios with a 17% target, the representativeness of the target conservation features for the 30% and 50% target are higher in scenario 'b'. Scenario 'b' increased the representative of area with the top quantile for carbon storage to 17,734 km2 (30% target) and 31,982 km2 (50% target) (Table S4; Fig. 4), an increase from 17.6–38.4% and 69.3%, respectively. In scenario ‘b’, the protection of low elevation hill forest, as the forest type with the smallest area under protection, increased from 11.7–15.5% (3,893 km2, 30% target) and 48% (12,024 km2, 50% target).
Based on their irreplaceability, the analysis identified several priority areas (Fig. 5). The weighted Kappa statistic indicated the degree of agreement between the two scenarios was lowest for the 17% target (κ = 0.41). This differed most between scenario ‘b’, where PA inclusion is not mandatory in the selection process, that selected a smaller area inside the PA (1,538 km2), compared to its equivalent under scenario ‘a’ that required 18,602 km2 to meet the same target (Figs. 3 & 5). There was moderate concurrence between scenarios ‘a’ and ‘b’ for the 30% target (κ = 0.77) and high concurrence between these scenarios for the 50% target (κ = 0.92). Under scenarios ‘a' and ‘b,' the representativeness of priority areas was higher outside the PA network than inside. Focusing on the areas outside of the PA that we recommended be considered for inclusion in current PA network, the results showed that scenario 'b' has a larger area (Table 2).
Table 2
Recommended order for assigning protection to priority areas outside of Sulawesi's current protected area network.
Order
|
Scenario ‘a’ (km2)
|
Scenario ‘b’ (km2)
|
First
|
2,263
|
3,368
|
Second
|
10,776
|
17,926
|
Third
|
35,697
|
29,926
|
Total area
|
48,736
|
51,220
|
The total size of priority areas identified outside of the current PA network is 48,736 km2 in scenario ‘a’ and 51,220 km2 in scenario ‘b’ (Table 2). These vast areas would require nearly three times the current land area being allocated for PAs and in order to achieve this we propose a sequential three-stage implementation plan. The priority areas that are aligned with Aichi Target 11 are represented in the first order (17%). Scenario ‘a' identifies 2,262.9 km2 for inclusion in the PA network, primarily through expanding the current PA network boundaries (Fig. 5).
The second order in scenario ‘a’ aligns with Target 3 of the first draft of the post-2020 GBF by: i) adding a buffer area to the existing PAs; ii) creating new stand-alone PAs; and iii) creating corridors that connect pre-existing PAs. Scenario ‘a’ focused on existing PA buffer zones, whereas scenario ‘b’ identified large stand-alone areas, such as Mekongga (1,576 km2) in cluster 5, Banggai (1,004 km2) in cluster 4, and Popayato-Paguat (1,299 km2) in cluster 2 (Fig. 5). These three large areas were also identified in scenario ‘a’ but as a second order selection priority with a similar spatial configuration for Mekongga and Banggai, but with Popayato-Paguat forming a corridor connecting two nearby PAs and two KBAs.
Outside of the current PA network, the following conservation features were chosen for each cluster, for which we consider possible impacts from different land management sectors.
North Sulawesi (Cluster 1)
The main first order forest type selected was sub-montane forest, with scenario ‘b’ selecting the inner part of Karakelang WR and Bogani Nani Wartabone NP (BNWNP). The remaining first order areas in scenario ‘a’ were all PAs. The buffer of Karakelang WR and BNWNP, mostly consisting of low and medium elevation hill forest, was in the second order for scenarios ‘a’ and ‘b’. In the third order, the notable additions were most of Sangihe island and three patches of forest outside Gunung Ambang NR with high carbon stock (Fig. S2). Karst ecosystem was not a main feature in this cluster. An identified impact in this cluster is from the operations of the PT Tambang Mas Sangihe concession that covers half of Sangihe Island. The priority area identified along the northern boundary of BNWNP is located within the PT Huma Sulut Natural Forest Management Permit area, meaning that it may be logged in the future.
Gorontalo and Central Sulawesi (Cluster 2)
This cluster contained a large portion of the top quantile for carbon stock that was located outside of the current PA network. Overall, the priority areas cover all the high carbon areas. The first priority under scenario ‘b’ includes a patch of lowland elevation hill forest (360 km2) as the primary habitat type (Fig. 5), while scenario ‘a’ mostly includes current PAs (Table S1: ID 14, 15, 19, 20, and 21). The remaining priority areas are a network of corridors that connect PAs and KBAs. In scenario ‘a’, a second order constitutes a forest area of 2,742 km2. It is dominated by low and medium elevation hills (1,497 km2) and includes a 1,299 km2 corridor connecting Nantu WR and Panua NR, as well as the Papayato-Paguat KBA and Gunung Ile-Ile KBA (Fig. S2). The second order area in scenario ‘b’ (4,276 km2) extends from Gunung Dako NR to Tinombala NR and includes the Buol-Tolitoli KBA. It is dominated by medium elevation hills (1,621 km2). This cluster includes the logging concession areas (HPH) of PT Sentral Pitulempa and PT Taman Hutan Asri and the oil palm concessions of PT Inti Global Laksana and PT Cipta Cakra Murdaya. Additionally, the PT Gorontalo Sejahtera mining concession area overlaps with the southern part of Panua NR.
Central Sulawesi (Cluster 3)
Outside of the PA network, the second order priority areas from scenario 'a' are a 3.223 km2 area dominated by lower montane forest (1,114 km2). These areas are in the buffer zones of Gandang Dewata NP (ID 37) and Faruhumpenai NR (ID 40), with standalone areas in the central Pegunungan Tokalekaju KBA and Pambuang KBA, the southern extent of Lore Lindu NP, and the northwest of Morowali NR. Under scenario 'b,' these last two areas were also chosen as first-order priorities. Potential impacts in this cluster include the overlap of Pegunungan Tokalekaju KBA with the Natural Forest Management Permit area of PT Satya Sena Indratama, Pambuang KBA with the Natural Forest Management Permit area of PT Inhutani I Mamuju, and Sulawesi Tengah Grand Forest Park with the mining concession PT. Citra Palu Mineral. Another area of this company’s concession is in Luwu Utara, to the South of the Lore Lindu NP.
Eastern Sulawesi (Cluster 4)
The 1,004 km2 Banggai landscape is a second-order priority area located in the western part of Lombuyan I and II WR is the most notable area in this cluster (ID 35). The most common habitat types are sub-montane forest (734 km2), karst ecosystem (870 km2), and a carbon-rich area. Around their respective borders, Bangkiriang WR and Morowali NR have first order priority areas. The Bangkiriang priority area is linked to a second order priority area (209 km2) that includes a portion of the Bangkiriang KBA. The presence of several plantation concessions with Natural Forest Management Permits, includes PT. Bina Balantak Raya, PT. Dahatama Adi Karya, PT. Palopo Timber Company, and PT. Satyaguna Sulajaya, which straddles the large patch of Banggai, as well as the mining concession of PT. Sinar Makmur Cemerlang inside Bangkiriang KBA.
Southeast Sulawesi (Cluster 5)
The 1,576 km2 Mekongga landscape is the most prominent feature in this cluster. It is a second order priority area that fully encompasses the Mekongga KBA’s northern portion. Sub-montane forest (1,041 km2) and karst ecosystem (1,158 km2) are the most common here. There are mining concessions, notably on Kabaena Island that has multiple mining concessions that cover the entire island.
South Sulawesi (Cluster 6)
This cluster has the smallest priority areas outside of the PA network of any of the clusters. The third order is the largest type of priority area of scenario ‘a’ in this cluster. This area of 799 km2 is dominated by montane forest and has the largest continuous forest patch (542 km2) that overlaps with the Pegunungan Latimojong KBA, which contains Sulawesi's highest mountain (Mount Latimojong). Scenario 'b' has second order priority areas (178 km2) outside of the PA network in the interior of the Pegunungan Latimojong KBA, which is dominated by tropical upper mountain and subalpine forest.