Disproportioned Performances of Protected Areas in the Beijing-Tianjin-Hebei Region

Designing and managing protected areas (PAs) has always been value-based. While all PAs provide unique values in a region, it has not been documented if their functional contributions are proportional to their land area. A challenge also arises from the governance of PAs due to differences in supporting resources, legislations, conservation goals, and guidance. Using the Beijing-Tianjin-Hebei (BTH) region as our testbed, we designed a study to test the hypothesis that the functional contributions of 84 PAs, measured by various ecosystem service values, vulnerability and human footprint (Hf), are not proportional to land areas relative to the regional total. Disproportionate performances are partially due to the governance (i.e., national, provincial and local) and the dominant ecosystem type of each PA. We confirm our hypothesis that the functional contributions of PAs in the BTH region far exceed their total land area, with forested PAs performing better. The functional measures of PAs are highly correlated but not equally provided among PA types and under different governances. Both governance and PA type play significant roles in determining functional values and human influences, but no significant interactive influences were detected. Results from this study are updated shreds of evidence for promoting PAs as well as for establishing new PAs through the ongoing Ecological Red Lines program.


Introduction
The concept of placing certain areas under protection was developed from the realization that there is unique and irreplaceable value in some landscapes [1]. The first official proposal to establish a protected area, led by Frederick Law Olmsted in 1865, argued that the Yosemite Valley and the Mariposa Big Tree Grove should be treated as parks so "the most certain value and gratification afforded by a park" could be experienced by people [2]. The proposal emphasized the importance of contact with the wilderness for human well-being, the effect of beautiful scenery on human perception, and the moral responsibility of governments to preserve regions of extraordinary natural beauty for the benefit of all people. The Organic Act of the US congress in 1916 officially established the National Park Service to "conserve the scenery and the natural and historic objects and the wildlife therein and to provide for the enjoyment of the same in such manner and by such means as will leave them unimpaired for the enjoyment of future generations". The act was largely influenced by John Muir's enthusiastic publications in 1913 and the persuasion of the government [3], as well as by public engagement (e.g., the establishment of the Sierra Club). This concept and associated practices have received global support from people and governments. By 2018, 8.3 million km 2 (21%) of global lands were under various forms of protection, reflecting a continued increasing trend [4]. In Europe, 1.34 million km 2 are

Study Area
The BTH region is located in the transitional zone between the Inner Mongolia Plateau and the North China Plain (Figure 1). The northwest is dominated by several mountains, which transition to plains in the southeast. Four ecological zones fall within the BTH region [24]: typical grasslands in the central and eastern Inner Mongolia Plateau, the Yanshan-Taihangshan mountain deciduous broad-leaved forest zone, urban agglomeration and sub-urban agricultural zone, and the agricultural zone.

PAs within BTH
Geographical location and associated information about PAs were extracted from the protected area database of National Forestry Administration [25], with additional statistics drawn from Wu et al. (2011) [9] and IUCN [26]. Within the database, some PAs were not described by their boundaries but by their total area and coordinates; we use symbols (Circles) to present these sites ( Figure 1). These PAs are mostly distributed in mountainous areas in the northwest parts of the region. The 84 PAs include all the nature reserves and some of the forest parks, wetland parks and geological parks in BTH, covering 10,413 km 2 (4.83% of BTH) ( Table 1). These PAs are administrated by national (23), provincial (42), and local (19) governments. The mean (STD) of national PAs (93 ± 129 km 2 ) is larger than provincial PAs (74 ± 123 km 2 ) and local PAs (31 ± 41 km 2 ). In terms of dominant cover type within each PA, forest PAs are the most common, in 44 PAs, which cover 5989 km 2 (57.5% of the total). On average, grassland-dominated PAs are the largest (205 km 2 ), compared to forest (75 km 2 ) and others (40-52 km 2 ). Among the three provinces, 45 PAs are in Tianjin, 28 in Beijing and 11 in Hebei (Table 1).

PAs within BTH
Geographical location and associated information about PAs were extracted from the protected area database of National Forestry Administration [25], with additional statistics drawn from Wu et al. (2011) [9] and IUCN [26]. Within the database, some PAs were not described by their boundaries but by their total area and coordinates; we use symbols (Circles) to present these sites ( Figure 1). These PAs are mostly distributed in mountainous areas in the northwest parts of the region. The 84 PAs include all the nature reserves and some of the forest parks, wetland parks and geological parks in BTH, covering 10,413 km 2 (4.83% of BTH) ( Table 1). These PAs are administrated by national (23), provincial (42), and local (19) governments. The mean (STD) of national PAs (93 ± 129 km 2 ) is larger than provincial PAs (74 ± 123 km 2 ) and local PAs (31 ± 41 km 2 ). In terms of dominant cover type within each PA, forest PAs are the most common, in 44 PAs, which cover 5989 km 2 (57.5% of the total). On average, grassland-dominated PAs are the largest (205 km 2 ), compared to forest (75 km 2 ) and others (40-52 km 2 ). Among the three provinces, 45 PAs are in Tianjin, 28 in Beijing and 11 in Hebei (Table 1).

Validation of Performances
Six indices as conservation objectives for regional importance were chosen. To justify the most appropriate conservation objectives, we first collected all relevant indices after reviewing over a hundred studies in BTH [27][28][29][30]. These indices were further evaluated and limited to six based upon typical indicator selection criteria of specificity, portable, measurability, sensitivity [31][32][33]. We calculated six indices as conservation objectives (COs). Four are ecosystem services: climate regulation (Cr), water retention (Wr), soil retention (Sr), and wind prevention (Wp). Two are for ecosystem vulnerability: soil erosion (Se) and desertification (De). We also calculated the human footprint index (Hf), which reflects the level of human activities as a stressor.
Water retention (Wr) volume is calculated using the water balance equation in InVEST [34,35]; the amount of soil retention (Sr) and the sensitivity index of soil erosion (Se) are calculated by the RUSLE model [36,37]. Rainfall data is extracted from the 1-km monthly and annual precipitation data of the Human Geosystem Sedisystem Sedodatabase (http://www.data.ac.cn/). The yield ratio (K) comes from the literature [36,38]. Vegetation cover is quantified from the 1-km NDVI Monthly Synthesis (MOD09A1) of NASA's EOS data center for 2000-2009 (https://modis.gsfc.nasa.gov/). The vegetation type is from the 1-km global land cover of the University of Maryland (http://glcfapp.glcf.umd.edu:8080/esdi/index.jsp). Soil type and soil texture are derived from the 1:1 million soil space (94,000 plots) attribute data of the second National Land Survey by the Institute of Soil Science, CAS. The 90-m Digital Terrain Elevation Model (DEM) is from the International Scientific Data Mirror Network (http://www.gscloud.cn) of the Computer Network Information Center, CAS. The precipitation, runoff ratio, soil, DEM and MODIS data for 2000-2009 are scaled to 1-km as required input variables in our calculations.
The climate regulation service index (Cr) is extracted from the 1-km dataset for China, which was developed by the Institute of Geosciences and Resources, CAS, in 2015 [39]. The wind prevention (Wp) and Desertification (De) in 2015 are from the 250-m grid data of the Institute of Environmental Science Planning and Design [28]. The 2009 human ecological footprint index (Hf) is 1-km grid of Venter, which is based on population density, land cover type, accessibility, and electricity infrastructure to model human impacts [40]. All the data are synchronized to 1-km grid spatial data for BTH within the ArcGIS platform.

Statistical Analysis
To answer the above questions and achieve our study objectives, we performed statistical analyses to test two specific hypotheses: (1) there exist significant differences in the functional contributions of PAs, which are highly dependent on a PA's governance and type; and (2) the ratio between conservation objective values of the PAs (CO PA ) and the total CO values of the BTH region (CO total ) is disproportional to their land areas (A PA ), wherein CO PA :CO total A PA :A total . Additionally, these disproportions vary by PA governance and type. Two-way analysis of variances (ANOVA) is performed to quantify the independent and interactive influences of governances (three levels) and types (forest and non-forest) on each CO measure. When interactive influences are not significant (p > 0.05), we proceed with the independent effects of governance and PA type for each objective value. As local management practices are based on multiple functions instead of individual ones, the CO values at county and planning unit (PU) levels may correlate with each other because of similar governance and PA type. There are 175 counties and 53,661 PUs within BTH; to provide useful evidence for future management, we calculated the correlations of CO values at PU level (i.e., 2 × 2 km).

Landscape Composition of the PAs and BTH
Farmlands account for 45.5% of the BTH region, followed by forests (22.6%), grasslands (18.1%), and urban build-ups (10.0%) ( Figure 2). In comparison, PAs have a composition of 49.7% forests, 19.1% farmlands, and 18.1% grasslands. The proportion of water bodies within the PAs is much higher (5.8%) than that of the region (1.8%). Not surprisingly, the ratio of urban build-ups within the PAs (4.7%) is less than half of the regional value. It is worth noting that the portion of wetlands within the PAs is small (0.2%) but still several orders of the magnitude higher than that in the BTH ( Figure 2). The landscape composition of PAs varies substantially by governance. National and local PAs have higher proportions of forests than provincial PAs, whereas water and grasslands are more common in provincial PAs than in national and local PAs. Relatively high proportions of farmland and built-up areas are also apparent in local PAs. (COtotal) is disproportional to their land areas (APA), wherein COPA:COtotal ≠ APA:Atotal. Additionally, these disproportions vary by PA governance and type. Two-way analysis of variances (ANOVA) is performed to quantify the independent and interactive influences of governances (three levels) and types (forest and non-forest) on each CO measure. When interactive influences are not significant (p > 0.05), we proceed with the independent effects of governance and PA type for each objective value. As local management practices are based on multiple functions instead of individual ones, the CO values at county and planning unit (PU) levels may correlate with each other because of similar governance and PA type. There are 175 counties and 53,661 PUs within BTH; to provide useful evidence for future management, we calculated the correlations of CO values at PU level (i.e., 2 × 2 km).

Landscape Composition of the PAs and BTH
Farmlands account for 45.5% of the BTH region, followed by forests (22.6%), grasslands (18.1%), and urban build-ups (10.0%) ( Figure 2). In comparison, PAs have a composition of 49.7% forests, 19.1% farmlands, and 18.1% grasslands. The proportion of water bodies within the PAs is much higher (5.8%) than that of the region (1.8%). Not surprisingly, the ratio of urban build-ups within the PAs (4.7%) is less than half of the regional value. It is worth noting that the portion of wetlands within the PAs is small (0.2%) but still several orders of the magnitude higher than that in the BTH ( Figure  2). The landscape composition of PAs varies substantially by governance. National and local PAs have higher proportions of forests than provincial PAs, whereas water and grasslands are more common in provincial PAs than in national and local PAs. Relatively high proportions of farmland and built-up areas are also apparent in local PAs.

Conservation Objective Value of the PAs and BTH
The average unit value of five COs in the PAs is higher than those of BTH; the De value is the same as for BTH region. As expected, the human footprint (Hf) of the PAs (14.1) is lower than that of BTH (18.4) ( Table 2). Overall, PAs provide the highest relative soil retention (Sr) at 7.75% and the lowest relative contribution to water retention (Wr) at 5.04%. The human footprint (Hf) of all PAs is 3.70%, which is far below that of the 4.83% of BTH land dedicated to PAs (Table 1). When compared among governances, the CO values per unit of land area in national PAs show higher amounts than those of the region, whereas provincial and local PAs provide similar or even lower values. The Hf of local PAs (16.0), for example, is surprisingly higher than that of national and provincial PAs (13.8-14.0), as well as the regional mean. Lower CO values in provincial PAs are coupled with their overall large sizes (Table 1), while local PAs have a substantially higher proportion of build-ups ( Figure 2) and human footprint ( Table 2). It is worth mentioning that Wr and Sr values per square kilometer at national and local levels are relatively higher than they are on a provincial level, while Se and De, the

Conservation Objective Value of the PAs and BTH
The average unit value of five COs in the PAs is higher than those of BTH; the De value is the same as for BTH region. As expected, the human footprint (Hf) of the PAs (14.1) is lower than that of BTH (18.4) ( Table 2). Overall, PAs provide the highest relative soil retention (Sr) at 7.75% and the lowest relative contribution to water retention (Wr) at 5.04%. The human footprint (Hf) of all PAs is 3.70%, which is far below that of the 4.83% of BTH land dedicated to PAs (Table 1). When compared among governances, the CO values per unit of land area in national PAs show higher amounts than those of the region, whereas provincial and local PAs provide similar or even lower values. The Hf of local PAs (16.0), for example, is surprisingly higher than that of national and provincial PAs (13.8-14.0), as well as the regional mean. Lower CO values in provincial PAs are coupled with their overall large sizes (Table 1), while local PAs have a substantially higher proportion of build-ups ( Figure 2) and human footprint (Table 2). It is worth mentioning that Wr and Sr values per square kilometer at national and local levels are relatively higher than they are on a provincial level, while Se and De, the two vulnerability values, are higher on a provincial level than the other two levels. When PAs are divided into forest and non-forest groups, CO values (except De) are substantially higher in forest-dominated PAs, and Hf values are higher in non-forest PAs. In sum, differences in CO values and human footprint vary by governance and type; local PAs appear to be more influenced by human activities (e.g., build-ups). Table 2. Values of six key conservation objectives and human footprint index (Hf) of the protected areas (PAs) in BTH region. The contribution of each conservation objective (CO) value is calculated as the CO per square kilometer (the measurement unit of climate regulation (Cr) is RMB, Wr is m 3 , and soil retention (Sr) and wind prevention (Wp) is ton; soil erosion (Se), desertification (De) and Hf are unitless; Se and De range 0-5.5 and 0-12.6, respectively; and Hf score ranges 12-96 within BTH), and its proportional contribution (%) is calculated as the CO values of the PAs to the total CO values of BTH. Our linear ANOVA indicates there is no significant interactive effect (p < 0.05) of governance (G) and PA type (E) on the six CO values and the human footprint ( Table 3). The ANOVA models explained small amounts of variation, with R 2 varying from 0.112 for De to 0.338 for Sr. The eta squared values further indicated that PA type is the primary independent variable in comparisons with governance and the interactive effect (E*G). For Cr, Sr and Se, PA type is responsible for 23.5%, 33.1% and 21.5% of the total variance, respectively. For human footprint, both G and E are significant independent influences. For the main effect, the p-value for Wp is <0.05 (Table 3), with the partial least square indicating that 11% of the variation is from PA type and 8% from governance. For De, the type is only significant (~11%) for provincial and local PAs ( Table 4). As for Hf, both governance and type show significance; PA type contributes 31.8% of the variation, compared with 14.2% for governance. For eight significant contributions (bolded data in Table 3), there is a significant difference in Cr, Wr, Sr, Wp, and Se between forest and non-forest PAs with an overall variation of 654 RMB, 103 m 3 , 62 t, 0.052 t, and 0.788, respectively (Table 4). Whereas for Wp, there are significant differences between national and local PAs, with an average difference of 0.057 t, and for De between provincial and local PAs with a difference of 0.736. Table 3. Statistics of the two-way ANOVA on CO valuesper unit land area with three levels of governance and two levels of PA type (n = 84). "E": dominant ecosystem type (forest and non-forest). "G": Governance (national, provincial and local). "E*G": Interactive effect of E and G.  Table 4. Statistical tests of conservation objective (CO) differences (t-test) when governance or type has a significant independent influence on each CO value at p < 0.05 (shown in Table 3). "F-NF": test between forest and non-forest PAs. "N-L": test between national and local PAs; "P-L" test between provincial-local PAs. "MD": mean difference.

Functional and Structural Contributions of PAs
Among the 23 national PAs, 12 are classified as forest (F-N) and 11 as non-forest (NF-N). For provincial PAs, 20 are forest (F-P) and 22 are non-forest (NF-P). Local PAs include 12 forest (F-L) and 7 non-forest (NF-L). The ratios of CO values from the PAs to those of regional means (i.e., CO PA :CO total ) are higher than the ratios of land area (i.e., A PA :A total ). Although there exist exceptions and variations by governance and PA type, our hypothesis is accepted for all seven measures of PA functions (Figure 3). The majority of forest PAs (in green) had higher functional contributions than their land portions, with many showing near a 2:1 ratio (e.g., Cr, Sr and Se) (Figure 3a,c,d). For De and Hf, the functional contributions of forest PAs are about 1-2 times of their area contributions (Figure 3e,f). For Wr, all forest PAs show near equal functional contributions (Figure 3b). The only PA showing a low ratio of Wr and Sr is a Xiaoshan volcanic geological heritage provincial protected area. For non-forest PAs, there appear to be no clear-cut conclusions. Most ratios for Cr and Sr are lower than their area proportion (Figure 3a,c), and their ratios for Wr and Hf are similar to the area ratios (Figure 3b,d,f). However, the De ratios are 1-2 times their area ratios (Figure 3e).
Governance also plays an important role in the relationships between CO PA :CO total and A PA :A total . National forest PAs (green triangle) have ratios above

Correlations among ES Measures
There exist variable correlations in both magnitude and direction among the seven performances of the PAs at three levels: PA, county and planning unit (PU) ( Table 5). Five correlations were not significant: De vs. Sr and De vs. Hf at PA level, Hf vs. De at both PA and county level, and Sr vs. Wr at county level. Overall, the correlations at county level are higher than those at PA and PU level.

Values of Protected Areas
Designing and managing protected areas has always been value-based. This is done by reducing or restricting human activities. Through estimating ecosystem services, vulnerability and human influences of 84 PAs within the BTH region, we confirm our hypothesis that the functional contributions of the PAs in the BTH region far exceed their total land area ( Figure 3). Some ecosystem services are double of their area portion, though with exceptions. These estimates provide us promising evidence supporting the success of PAs in BTH, which were developed mostly within the past three decades [24], albeit the first PA in China was established in 1956.
Regardless of the overall disproportionate performance of PAs, there exist substantial differences among the seven measures we examined, as well as the PA types. Compared to the 4.83% area contribution, all six indicators of the PAs provide 5%-7.8% contribution, while the human impact of the PAs (3.7%) is lower than that of the region. Most PAs worldwide were initially designed to protect biological diversity amid an increasing recognition of their other functions such as enhancing ecosystem services and regulating regional climate [41,42]. Within the BTH, 72 of the 84 PAs were originally designed to conserve species. Through testing a hypothesis on the area-function ratio, we demonstrate that PAs provide more value than other land cover types within the BTH region.
The average values of all four ecosystem services are higher for PAs compared with the regional means. This is coupled with lower human influences within the PAs (Table 2). However, we found that non-forest PAs did not always provide expected values when compared with regional means (Figure 3). This is especially true for their low contributions to climate regulation, soil regulation, and wind prevention. A possible reason is that BTH is in a temperate deciduous forest zone [18]. Restored non-forest PAs from croplands within the BTH would function differently than forests.
Finally, it is worth reiterating the correlations among functional measures of PAs because ecosystem services, vulnerability and human influences are not provided equally by each PA. Four pairs of metrics (Sr-Se, Cr-Sr, Cr-Se and De-Se) are highly correlated. While the high correlation of Sr-Se is likely because they are predicted with the same model, Cr's correlations with all other indicators (e.g., Wp and De) are reasonable because both wind erosion and desertification are related to climatic changes. In sum, PAs provide multiple functions in addition to the preservation of biodiversity. Integrating ecosystem functions and ecosystem services for assessing existing PAs and establishing new PAs is necessary in order to promote collaborative conservations and enhance the efficiency of conservation and management.

The Role of Governance and Type
Our second hypothesis is also accepted, i.e., both governance and PA type play significant roles in determining functional values and human influences. However, our initial premise was that governance might dominate over PA type in explaining variations because national PAs receive more resources from the government and have more restricted rules for human activities. However, this is proven untrue by our results. For Cr, Sr and Hf, PA type explains 20-33% of the total variation, whereas governance explains 8%-11%. The more important role of PA type compared with governance, meanwhile, emphasizes that the initial decision in establishing PA type is very crucial for its subsequent functions. Another surprise in this study is the insignificant interactions between governance and PA type. Forest PAs are one of the most protected ecosystems on a global level because of their great contributions to climate mitigation [43,44]. Our results confirm these continued efforts on the part of forest PAs within the BTH region.
The proportion of forests and water bodies in the PAs within the BTH is more than twice the regional average, while the proportion of grassland is similar to the regional mean. Plausibly, forest PAs accounts for the largest proportion, and water bodies account for the largest proportion in provincial protection spaces. Water bodies, meanwhile, are less sensitive to the area of protection than the others. In general, the trend among forest groups is obvious compared with non-forest groups, and forests are more concentrated even among human influence factors. The points in the forest group are basically sensitive to area increases, and certain characteristics, especially soil preservation and climate regulation, increase even more than twice the proportion of the increase in area. Nevertheless, this trend varies with different objectives, especially within the vulnerable group.
The governance, nonetheless, plays a significant role in Wp and De, especially on a local level. The overall function of national PAs is higher than those of provincial and local PAs ( Table 2). The primary reason for this may be higher supporting resources and stricter regulations, such as access to and activities allowed inside the parks. On the other hand, the Wr and Sr of local PAs are higher than those of provincial and national PAs due to a stronger interest in reducing wind and soil erosions through management. This suggests that management practices (and regulations) applied to national PAs can potentially be applied to provincial and local PAs. Meanwhile, local PAs are often closer to human populations or accessed more by residents. This means that they provide spaces and other functions to people that may not be replaceable by provincial and national PAs [45]. The overall higher services of Wr and Sr from local PAs are direct evidence for promoting more local PAs for human well-being. More research is needed to integrate social and economic elements on smaller scales [46].
Localized protected areas are closer to human settlements and have higher human impact. Their contribution to key ecosystem functions on a regional scale is generally lower, but this does not mean it is not important [47]. In holistic national or regional spatial planning, these localized protected areas scattered around human settlements not only provide diversified value to meet human needs but also contribute a continuous flow and network for wholesome ecosystem functions in the region, such as the ecosystem services and functions in this study. In large-scale spatial planning, it is necessary to recognize their roles and combine their contribution. For these close-knit protected areas, more research is needed to confirm their significance, as well as modes and patterns that are more meaningful on local, regional and national levels. Future studies should further identify and optimize the design of such spaces in the study area.

Implications for Future Research
This study sheds light on the importance of existing PAs in supporting the six conservation objects in BTH and achieving multifunctionality of PAs, which is increasingly important in China [48]. Effective planning and governance of PAs are much more complicated than we revealed in this paper, and we would suggest that future research addresses the following issues further. First, biodiversity is the primary goal of most current PAs, and more studies are needed to understand whether biodiversity has been achieved and what the relationships between biodiversity and other ecological objectives are. Second, aesthetic and recreational values have not been identified as important in many studies of BTH but have been highlighted in regional studies from other countries [49,50]. Third, conflicts between PAs and other land uses are pressing issues to be addressed. Conservation areas in BTH have been and will be in strong competition with other land uses such as urbanization. Maintaining the quantity and the quality of PAs must be placed in the context of other land uses and developments.
Achieving multiple conservation objectives including biodiversity, ecosystem functions and services is expected when designing and developing alternative management plans for PAs. Here, we provide updated evidence for PAs' performances within BTH. Lessons learned from this study will be useful for governments who are gearing up to increase the number of PA (i.e., from current 18% to 20% in the future). The evolving philosophy exists to enhance cumulative values through integrating existing and future PAs [27]; for example, the ongoing Ecological Red Lines project in China calls for multi-objective, systematic ecological conservation planning [51]. With rapid changes to the global climate and escalating human activities, many unprecedented processes were not foreseen by pioneer conservationists and policymakers when managing and designing PAs [52]. As global ecosystem degradations continue [53][54][55], public support for promoting PAs and adopting new conservation strategies is mounting. Pressing issues include addressing the conflicts between PAs and other land uses, increasing the quantity and quality of PAs, sustaining ecosystem services, and offsetting the effects of ecosystem degradation [56]. Our findings shed a bright light on these issues, at least for PAs within the BTH region.

Conclusions
Through quantifying the performances of protected areas (PAs) in the BTH region, we confirmed the effectiveness of PAs that provide overall high ecosystem service values, low vulnerability and reduced human footprint. Theses functional measures are not proportional to their land areas relative to the regional total. The majority of forest PAs had higher functional contributions than their land portions, with many showing near a 2:1 ratio. Disproportionate performances are found to vary by governance (i.e., national, provincial and local) and the dominant ecosystem type of PAs. These PAs provide the highest relative soil retention (Sr) at 7.75% and the lowest relative contribution to water retention (Wr) at 5.04%, which is still higher than the 4.83% of land areas dedicated to PAs. Meanwhile, the human footprint of all PAs is 3.70%. Surprisingly, we found governance and type of PA determined PA performances independently (versus their interactive effects), with national forested PAs performing the best. This study provides updated evidence supporting the concept of and continued efforts in promoting PAs. Results from this study can be useful for ongoing efforts in increasing the number, size and quality of protected areas within the BTH region, as well as for developing sustainable PAs elsewhere, achieving multiple conservation objectives through the engagement of all stakeholders.