Are Valuable and Representative Natural Habitats Su ﬃ ciently Protected? Application of Marxan model in the Czech Republic

: The joint impact of human activities and climate change on natural resources lead to biodiversity loss. Therefore, it is important to select protected areas through systematic conservation planning. The present study assessed how representative natural habitats are protected under the nature conservation network, and to identify new—but so far insu ﬃ ciently—protected areas containing these habitats for sustainable management. We used the Marxan model to select the most valuable insu ﬃ ciently protected natural habitats in the Czech Republic as a representative example for a conservation strategy for Central–Eastern European environments. We set three conservation targets (25%, 50%, and 75%), deﬁning how much percent area of valuable representative natural habitats should be added to the area of the habitats already included in the Nature Protection Network. To implement these conservation targets it is necessary to preserve 22,932 ha, 72,429, ha and 124,363 ha respectively of the conservation targets occurring in the insu ﬃ ciently protected areas, and 17,255 ha, 51,620 ha, and 84,993 ha respectively of the conservation features in the areas without protection status. Marxan was revealed to be an appropriate tool to select the most valuable and insu ﬃ ciently protected natural habitats for sustainable management. national with nature and landscape conservation


Introduction
The increasing human population and the rising demand of natural resources as world domestic product, global tree cover loss, fossil fuel consumption, and increasing of carbon dioxide (CO2) emissions and increasing the temperature together with the impact of global climate change clearly and unequivocally show that the planet Earth is facing a climate emergency. It is essential to prevent the collapse of ecosystem functioning by protection and restoration of ecosystems and biodiversity [1,2]. Therefore, it is important to design protected areas to maintain minimally a present status of biodiversity. One way of promoting biodiversity conservation is to implement systematic conservation planning by effectively defining protected areas in both natural and agricultural landscapes. In addition, the newly designated protected areas can provide ecosystem services and biodiversity protection [3]. Climate i.e., 23% of country land. The Natura 2000 sites cover about 21% of the land and marine waters territory of the European Environment Agency's member countries and collaborating countries [27]. The European Environmental Agency (EEA) is responsible for developing The European Environment Information and Observation Network (Eionet) and coordinates its activities together with National Focal Points (NFPs) in the countries. The Eionet constellation is supplemented by seven European Topic Centres [28].
In our study, we divided the various categories of protected areas into sufficiently and insufficiently protected according to the level of habitat biodiversity protection (see Table 1). The distribution of protected areas in the Czech Republic is shown in Figure 1. The aim of our study was to assess how the valuable representative natural and near-natural habitats, identified by the national Habitat Mapping of the Czech Republic project, are protected under the nature conservation network, and to identify the new areas containing valuable natural habitats for sustainable management, not yet sufficiently protected. Table 1. The status of Planning Units (PUs) based on the degree of protection under the Nature Protection Acts of the Czech Republic (Conserved-sufficiently protected areas, Available-insufficiently protected areas, Excluded-excluded from protection).

Conserved
Available Excluded

A Description of the Marxan Model
The Marxan model was designed by Ball et al. [14] to plan new representative and spatially compact protection sites for landscape planning. The model uses a combination of Planning Units (PUs) that include all defined conservation features (CFs), minimising the total cost of the selected PUs. The Marxan model divides the area of interest into PUs that have either hexagonal or square meshes [29]. The PUs are divided into (i) conserved units-situated in sufficiently protected areas; (ii) available units-which can be selected if they include suitable conservation features; (iii) excluded units-excluded from the selection (e.g., built-up areas). Areas with strict protection of habitats located in larger protected areas or having a larger protection zone around them are considered as sufficiently protected areas. This category includes: 1st and 2nd zones of NPs, 1st zone of LPAs, National Natural Reserve and small-scale plots of National Natural Monuments, Natural Reserves, and Nature Monuments, situated in the LPAs or NPs. The Marxan model uses the "object function" to calculate the optimal selection of PUs by comparing the selected combinations of PUs. For each combination of PUs, the software adds up the total costs needed to create them. It performs ongoing combination tests to achieve the lowest possible object function with the lowest costs.
The optimisation algorithm within Marxan attempts to find good systems of sites through simulated annealing, whereby different sets of potential conservation areas are compared with user-defined targets and costs, and the set of areas that achieves its objective most efficiently is determined. With the use of stochastic optimisation routines (simulated annealing) it generates spatial reserve systems that achieve particular biodiversity representation goals with reasonable optimality. Computationally, Marxan provides solutions to a conservation version of the 0-1 knapsack problem, where the objects of interest are potential reserve sites with given biological

A Description of the Marxan Model
The Marxan model was designed by Ball et al. [14] to plan new representative and spatially compact protection sites for landscape planning. The model uses a combination of Planning Units (PUs) that include all defined conservation features (CFs), minimising the total cost of the selected PUs. The Marxan model divides the area of interest into PUs that have either hexagonal or square meshes [29]. The PUs are divided into (i) conserved units-situated in sufficiently protected areas; (ii) available units-which can be selected if they include suitable conservation features; (iii) excluded units-excluded from the selection (e.g., built-up areas). Areas with strict protection of habitats located in larger protected areas or having a larger protection zone around them are considered as sufficiently protected areas. This category includes: 1st and 2nd zones of NPs, 1st zone of LPAs, National Natural Reserve and small-scale plots of National Natural Monuments, Natural Reserves, and Nature Monuments, situated in the LPAs or NPs. The Marxan model uses the "object function" to calculate the optimal selection of PUs by comparing the selected combinations of PUs. For each combination of PUs, the software adds up the total costs needed to create them. It performs ongoing combination tests to achieve the lowest possible object function with the lowest costs.
The optimisation algorithm within Marxan attempts to find good systems of sites through simulated annealing, whereby different sets of potential conservation areas are compared with user-defined targets and costs, and the set of areas that achieves its objective most efficiently is determined. With the use of stochastic optimisation routines (simulated annealing) it generates spatial reserve systems that achieve particular biodiversity representation goals with reasonable optimality. Computationally, Marxan provides solutions to a conservation version of the 0-1 knapsack problem, where the objects of interest are potential reserve sites with given biological attributes. The simulated annealing algorithm attempts to minimise the total cost of the reserve system, while achieving a set of conservation goals (typically that a certain percentage of each geographical/biological feature is represented by the reserve system). The object function of the total cost of the solution is calculated according to the following formula: C -total cost of the final design solution, C PU -the costs of all planning units involved in the solution, SPF -total penalty for non-compliant conservation features, C BL -the length of the border surrounding the final protected area, BLM -a coefficient accounting for the complexity of the final protected area.

Planning Units (PUs)
Our analysis was carried out considering the whole of the Czech Republic as the study area. Country land was divided into PUs with a hexagonal net (edge of 310 m, area of one hexagon was 25 ha) that were divided into three categories: conserved, available, and excluded. At the beginning of the analysis, the nature protection effectiveness of different types of protected areas in the Czech Republic was critically assessed [30] (Table 1). To describe PUs in the conserved category, we used data from the Nature Conservation Agency in the Czech Republic (NCA CR). To delineate PUs in the available category, ZABAGED and OpenStreetMap and NCA CR data were used, while for identifying PUs in the excluded category, ZABAGED and OpenStreetMap only were used. All data were in the scale 1:10,000 with a projection of EPSG 5514.
The PUs in the Nature Protection Network (NPN) of the Czech Republic with sufficient protection have been categorized as conserved units, see Table 1. The excluded units included built-up areas, industrial areas, roads and their immediate surroundings. The Marxan model did not select any PUs from these categories. The rest of the Czech Republic was included in the available units ( Table 1).
The cost of individual PUs was calculated as the weighted average of the point value of biodiversity determined by the Habitat Valuation Method (HVM) (see Table 1) for all natural and non-natural habitats [31][32][33]. The Boundary Length Modifier (BLM) parameter determines the length of PU borders. For the interfaces Conserved vs. Excluded and Available vs. Excluded, the length of the PUs was increased three times.

Conservation Features
The Conservation Feature (CF) is a measurable, spatially defined element of biodiversity that is present in the final solution of the protected network; the CF is also included in the PU. CFs can be individual species and habitats, but also larger units such as biomes [34]. Each CF has a target which represents the minimum limit of CFs to be included in the final solution of the Marxan model. If the target is not met, the object function of the final solution is penalized.
We identified CFs as the natural and near-natural habitats of the Habitat Mapping from NCA CR [35]. The near-natural habitats represented habitats established by humans, but having species and structural diversity similar to nature habitats (e.g., cultural orchids meadows). Unit classification was defined by the Habitat Catalogue of the Czech Republic [36]. We used only those CFs which complied with the strict criteria of two quantitative parameters. The first parameter included (i) higher levels of representativeness and the well-protected status, and it was divided into the classes AA, AB, BA, BB, for the habitats valuated according to the original Methods of Habitat Mapping of Natura 2000 and the Emerald Network in 2000-2009 [37], and (ii) representativeness and degradation, according to the methodology for updated valuation of habitats [38]. The second parameter for each habitat segment was the achievement of the minimum area necessary for ecosystem self-regulation, derived from the HVM method for each natural and near-natural habitat [31]. A total of 139 habitat types were identified as CFs, expressed in 149,906 segments and with a total area of 470,939 ha.
We created three scenarios representing the percent area of CFs included in the final proposal of the Marxan model to meet different protection targets (25%, 50%, and 75%), considering valuable habitats from all selected CFs in the Czech Republic. If part of the CF was already protected in the conserved PU, the protected area was deducted from the target area. The Marxan model selected the conservation features in hexagons with the highest biodiversity point values determined by the HVM method. If two hexagons had the same point value, the Marxan model chose a hexagon with a shorter boundary length than the edge of the PU.

Technical Description of Marxan Calculations
Our analysis was done in QGIS 3.6 and in the Marxan version 2.43 (64b). QGIS 3.6 was needed to pre-process spatial data and visualize the results. Each analysis was run in 20 steps with 1,000,000 iterations. To create input files with control parameters for Marxan, the Conservation Land-Use Zoning software (CLUZ) in QGIS was used. The CLUZ is a QGIS plugin that assists in designing protected area networks and other conservation landscapes [39]. After the simulation, we analysed how the PUs conform to the scenarios, and compared them with the rest of the existing network of protected areas in the Czech Republic. The compliance analysis was performed using GIS overlay operations.

PUs and Protection Scenarios
The area of the Czech Republic for PUs was divided into 318,006 PUs with a total area of 7,950,150 ha ( Table 2). This area is larger than the original area of the Czech Republic (7,886,879 ha) because the model includes the hexagons that cross the national borders. Marxan unifies PUs (hexagons) into patches. The total number of patches was 897, and the median patch area was 100 ha ( Table 3). The smallest patch is an area of one PU containing at least one segment of a CF. The PUs in the available category covered the largest area; however, the PUs in the conserved category had the highest cost per ha according to HVM [31]. The most valuable habitats, defined in terms of their CFs, occupied an area of 470,939 ha, and made up 6% of the area of the Czech Republic. The spatial distribution of CFs and PUs is shown in Figure 2.  PUs containing valuable natural and near-natural habitats (CFs) not yet conserved in protected areas were selected using the Marxan model (Table A1, Appendix A). To obtain a better overview of proposed new CFs, we aggregated the selected valuable habitats into habitat formation groups with similar environmental conditions, e.g., all types of meadows to the grasslands group and all wetlands and peatbogs to the wetlands group [36] (Table 4). The largest area of natural and nearnatural habitats was covered by beech and oak, and oak-hornbeam forests, and mesophilic meadows with an area in the range of 90,260 ha to 108,969 ha, respectively. The smallest area with rocks and rubble, natural Pinus mugo scrubs and alpine grasslands was found in the interval of 1242 ha to 1465 ha. The results show that the groups of habitats with the largest areas are not sufficiently protected over 50% of their area. A different situation was found for the following habitats: natural alpine grasslands, Pinus mugo scrubs, bog forests, spruce forests, peatbogs, and springs. These habitats are mostly protected over 50% of their area and for sufficient protection it is enough to add 392 ha. Groups of habitats with less than 10% of the total area sufficiently protected featured oak and oak-hornbeam forests, macrophytic vegetation of standing waters, alluvial forests, alluvial meadows and mesophilic meadows. Table 5 shows the required target areas for individual natural habitats classified as CFs, to be added to the proposal of new protected areas according to the three scenarios. PUs containing valuable natural and near-natural habitats (CFs) not yet conserved in protected areas were selected using the Marxan model (Table A1, Appendix A). To obtain a better overview of proposed new CFs, we aggregated the selected valuable habitats into habitat formation groups with similar environmental conditions, e.g., all types of meadows to the grasslands group and all wetlands and peatbogs to the wetlands group [36] (Table 4). The largest area of natural and near-natural habitats was covered by beech and oak, and oak-hornbeam forests, and mesophilic meadows with an area in the range of 90,260 ha to 108,969 ha, respectively. The smallest area with rocks and rubble, natural Pinus mugo scrubs and alpine grasslands was found in the interval of 1242 ha to 1465 ha. The results show that the groups of habitats with the largest areas are not sufficiently protected over 50% of their area. A different situation was found for the following habitats: natural alpine grasslands, Pinus mugo scrubs, bog forests, spruce forests, peatbogs, and springs. These habitats are mostly protected over 50% of their area and for sufficient protection it is enough to add 392 ha. Groups of habitats with less than 10% of the total area sufficiently protected featured oak and oak-hornbeam forests, macrophytic vegetation of standing waters, alluvial forests, alluvial meadows and mesophilic meadows. Table 5 shows the required target areas for individual natural habitats classified as CFs, to be added to the proposal of new protected areas according to the three scenarios. Table 4. The representation of the most valuable habitats selected as CFs in the total area of PUs and Conserved PUs. The all-natural and near-natural habitat codes are shown in Table 5, and abbreviations of these habitats according to the Habitat Catalogue of the Czech Republic [36] are reported in Table 3A in Appendix A.  Rocks, Rubble A5, A6, S1.1, S1.2, S1.3, S1.4, S1.5, S2 8.3 269.39 543.53

Assessment and Proposed Extension of the Protected Network
The existing NPN in the Czech Republic is considered sufficient in the representation and distribution of individual types of protected areas. Nevertheless, even if a CF area of 73,805 ha, corresponding to 15.6% of all CFs, is included in the sufficiently protected part of the NPN, it is necessary to provide increased biodiversity protection for the remaining valuable natural habitats classified as CFs. According to our proposal including three scenarios to protect 25%, 50%, or even 75% of CFs, it is necessary to conserve them in 22,932 ha, 72,429 ha, and 124,363 ha respectively in the insufficiently protected areas, and in 17,255 ha, 51,620 ha, and 84,993 ha respectively in the non-protected areas ( Table 6). For the Czech Republic, the Marxan model selected 3.66%, 10.96%, and 18.05% of the most valuable habitats without any protection according to our three scenarios. The largest non-protected land patches are situated mostly in the mountain forest areas along the borders of the Czech Republic and in military areas, e.g., in the Doupovské Mountains in Northwest Bohemia. To achieve the abovementioned scenarios, especially the scenario of adding valuable habitats to 75% of the area, Marxan selected and added isolated segments throughout the country (Figure 3).
protected land patches are situated mostly in the mountain forest areas along the borders of the Czech Republic and in military areas, e.g., in the Doupovské Mountains in Northwest Bohemia. To achieve the abovementioned scenarios, especially the scenario of adding valuable habitats to 75% of the area, Marxan selected and added isolated segments throughout the country (Figure 3).

Discussion
Our results show that the most valuable habitats selected as CFs are included in the current NPN of the Czech Republic (Table 4). However, only some of the most valuable habitats are sufficiently protected in these areas, such as alpine grasslands (only A1.1-A8.2), bog forests (only

Discussion
Our results show that the most valuable habitats selected as CFs are included in the current NPN of the Czech Republic (Table 4). However, only some of the most valuable habitats are sufficiently protected in these areas, such as alpine grasslands (only A1.1-A8.2), bog forests (only L10.3, L10.4), wetlands (only M2.2), peat bogs and springs (only R1.5, R3.1), and Isoëtes vegetation (V6) ( Table 5,  and Table A2 in the Appendix A). According to the Habitat Catalogue of the Czech Republic [36], these selected habitats cover only a small area and have specific natural conditions. All these habitats (except for the habitats R1.4, R1.5 included in peat bogs and springs) belong to the network Natura 2000 sites, and they are included as SACs in the national law of the Czech Republic [19]. This network consists of important habitats that are protected by the member states of the European Union. These habitats are important for biodiversity conservation but, because of their small size, they have little impact on the Territorial System of Ecological Stability of the Landscape (TSES). The objective of planning and developing the TSES is to halt the unfavourable trend in the development of ecological stability and to hamper biodiversity loss in the culture according to Czech law No. 114/1992 [19]. However, our three scenarios show that the most of the remaining habitats are not sufficiently protected. Each scenario represents the percentage of CFs that should be included in the final proposal of Marxan to achieve the targets of 25%, 50%, and 75% of all selected CFs. The most frequent habitats, covering large areas but with only 10% of the area sufficiently protected in NPN, belong mostly to the following habitats: Hercynian oak-hornbeam forests (L3.1), ash-alder alluvial forests (L2.2), willow-poplar forests of lowland rivers (L2.4), wet Filipendula grasslands (T1.6), intermittently wet Molinia meadows (T1.9), and mesophilic Arrhenatherum meadows (T1.1). According to Bastian [40], natural habitats occurring in insufficiently or non-protected areas, such as intermittently wet Molinia meadows (T1.9) and ash-alder alluvial forests (L2.2), are important for supplying ecosystem services, e.g., provisioning, regulating, and cultural services. The national Nature Conservation Agency evaluated the state and importance of the habitats included in the Natura 2000 sites in the years 2007 and 2013 in the Czech Republic. The Agency reported a higher number of important European habitats in 2013 compared to 2007, while admitting a possible estimation bias due to changes in assessment methods during this period. In addition, the number of habitats in an unfavourable state, translated into classification according to the Habitat Catalogue of the Czech Republic [36], e.g., macrophyte vegetation of water streams (V4), intermittently wet Molinia meadows (T1.9), and wet acidophilous oak forests (L7.2), has increased in the year 2017. All alpine grasslands were the only habitats with a better status (A1.1-A8. 2) in the year 2017 [41].
The division of protected areas into sufficiently and insufficiently protected areas was carried out by our expert assessment, based on management of habitat biodiversity in protected areas and on the definition of protected areas in nature and landscape conservation law [19]. Sufficiently protected areas are large-scale sites where the main objective is to protect biodiversity and ecosystem processes and human intervention is excluded or strongly limited. Furthermore, some small-scale protected areas that are situated in the territory of large-scale protected areas and have a sufficiently large protection zone containing human disturbance, are also considered as sufficiently protected areas. Most sites in the Natura 2000 network tend to overlap with large-scale protected areas, such as NPs and LPAs. However, most of these overlapped areas include zones with low protection such as the 3rd zone of NPs and protection zone of NP, and 2nd to 4th zones of LPAs. In the case of significant habitats or species occurrence, a management contract can be signed with the landowner or tenant [19].
There are issues with the inefficient use of subsidies from the EU and insufficient control of nature protection by the central EU authorities [42], including problems with the non-functional compensation mechanism for private forests included in the Natura 2000 network [43]. According to Křenová and Kindlmann [42], the negative effects on Natura 2000 habitats and species are stronger in the states with the weak position of the state representatives, such as the Czech Republic or, more generally, Eastern European countries. Most of all, post-communist countries, e.g., Slovakia, Hungary and Romania, have been relatively unsuccessful in implementing Natura 2000 in agricultural landscapes. The main causes can be found in the discrepancy between local communities and authorities and the formal protection of sites, and the lack of financial support for management [44]. In addition, some habitats of European importance are fully or partly dependent on the continuation of traditional agricultural management [45].
The only important areas for biodiversity protection not included in the NPN are military training areas, including habitats with manmade irregular mechanical disturbances. These interventions result in fine-scaled mosaics of ecological conditions important for many species of plants and insects, above all butterflies [46]. Warren and Buttner [47] confirmed the importance of active military training areas as refuges for disturbance-dependent endangered insects in Germany. After the army abandoned these areas, their conservation value is threatened. To conserve the threatened species, Cizek et al. [46] recommended implementing alternative projects mimicking army activities, or sustainable management systems.
There are several ways to establish new protected areas in the Czech Republic: (1) by creating new protected areas in the currently most valuable unprotected habitats identified as CFs by Habitat Mapping of the Czech Republic (e.g., in previous military areas); (2) by connecting selected CFs with existing protected areas or their buffer zones, mostly in the forest mountainous regions along the state borders; (3) by improving the protection of insufficiently protected areas through effective financial investments [38] (e.g., supporting traditional mowing management for keeping valuable meadows [45,48] or using the regional seed mixtures [49]); (4) by supporting the natural species composition in near-natural forests by close-to-nature forest management (e.g., oak forests). This type of management leads to an emphasis on stability, productivity, diversity, and continuity of forest conditions [50].
The Marxan model has been used in a number of regional studies to design protected sites, including a wide range of habitats as anthropogenic land, terrestrial lowland forests, mangrove forests, mountain areas, lakes, the coral reef and open sea with a size of one PU of 10 km 2 [6,7,16]. In our study, we have reduced the scale of habitats such as grasslands, forests, wetlands, and water bodies to only 0.25 km 2 for each PU. Our results show that it is possible to use Marxan for large areas with relatively detailed PUs. Further research should determine which scenario provides an optimal representation of valuable natural habitats, that is sufficient to maintain sustainable state of biodiversity in the cultural forest-agricultural landscape.
For better protection of biodiversity, it is also necessary to improve the connectivity among protected areas with different degrees of protection. The green infrastructures as a terrestrial ecosystem may play a key role in improving the coherence and resilience of protected areas [27]. For example, contractual agreements can be used to support the protection and connectivity of valuable habitats in insufficiently protected areas or in the areas without protection status at all. Three military territories with valuable habitats have been included in the Natura 2000 network as Special Protection Areas (SPAs) and Special Areas of Conservation (SACs), and at the same time they had the status of Contract Protected Area under Czech legislation. Another example is from the Landscape Protected Areas (LPA) Beskydy, where the forests of the Czech Republic manage the zone II of the LPA, overlapping with the SACs as a zone I of the LPA. In this case, forest owners may be supported by a financial contribution from landscape programs of the Ministry of the Environment.
Preliminary investigations on Natura 2000-protected areas suggest that there is now a relatively good spatial and functional connectivity among sites across national borders. The Alpine Network of Protected Areas (ALPARC) and the Carpathian Convention framework for protected areas are good examples of over-boundary cooperation. Both areas contain the large protected areas with another type of protection and both areas almost touch. Part of the Carpathian protected areas passes through the Czech Republic [27].
Results from the Marxan model can also serve to verify some categories of the Nature Protection Network (NPN). Not all types of protected areas have been proposed by the Ministry of the Environment of the Czech Republic on the basis of sufficient background data. At the same time, stakeholders such as foresters and farm holders, as well as fishing and hunting associations, should be involved in a sustainable management of near-natural agricultural and forestry management in natural valuable habitats [51][52][53]. Sustainable management is a conservation type of management containing appropriate procedures for conserving nature, e.g., extensive mowing and grazing, selective forest management, supported by subsidy funds from the Ministry of the Environment. Non-governmental organizations should also be involved in this issue, in particular by disseminating information on the importance of promoting valuable habitats in the cultural landscape to government officials and the general public.

Conclusions
The existing Nature Protection Network in the Czech Republic was found to be only partly adequate in terms of effective protection of representative valuable natural habitats mapped by the Nature Conservation Agency of the Czech Republic. A total area of 73,805 ha (15.6%), of the most valuable natural habitats classified as CFs, is sufficiently protected at present from the remaining CF area. Therefore, according to three conservation targets aimed at preserving 25%, 50%, and 75% of the CFs, it is still necessary to preserve 22,932 ha, 72,429 ha, and 124,363 ha of the CFs included in the insufficiently protected area, and 17,255 ha, 51,620 ha, and 84,993 ha of the CFs in non-protected areas. The group of most valuable habitats that sufficiently protect more than 50% of their total area included small habitats with specific environmental conditions, e.g., alpine grasslands and bog forests. The sufficiently protected area of large-area habitats, such as oak and alluvial forests, alluvial and Mesic meadows, covered less than 10% of their total area. Nevertheless, these significantly represented habitats in the CR affected by human activity are still important for the provision of ecosystem services and require protection in the form of appropriate land management.
We proved that Marxan model is applicable over relatively large regions with relatively small size of planning units. The Marxan model is a candidate tool for verification of currently protected areas by the Nature Conservation Agency of the Czech Republic. According to our results, only large segments of valuable habitats without any protection have been found in military areas that are considered particularly important sites for biodiversity in Eastern Europe. Stakeholders managing the landscape and the non-governmental organizations should be involved in the protection of selected valuable habitats selected by the Marxan model in the insufficiently or even completely unprotected areas. A sustainable land management strategy promoting less intensive agriculture and close-to-nature forestry should be specifically designed for (and applied to) these areas. The Marxan model revealed itself to be a suitable tool for supporting large-scale strategic planning in larger administrative units up to the national scale in the Czech Republic.
Funding: This work was supported by the Ministry of Education, Youth and Sports of the Czech Republic within the CzeCOS, grant number LM2015/061.

Acknowledgments:
We would like to thank Lenka Štěrbová for preparation of the nature and near-nature habitat databases with assignment of the value and price according to Seják et al. [31].

Conflicts of Interest:
The authors declare no conflict of interest.