Historical Patterns of Natural Forest Management in Costa Rica: the Good, the Bad and the Ugly

We present an in-depth analysis of natural forest management practices in Costa Rica based on a new historical forestry GIS database encompassing five conservation regions in the country where selective logging has been prevalent for the last 18 years. Natural forest management refers to the selective logging practices for the production of timber in natural forest ecosystems. The study considers natural forest management plans developed following a standardized format, based on Forestry Law 7575, established in 1996. Our results show a positive aspect of forest management where a large percentage of the managed forest is still standing, contrary to other tropical countries where selective logging promotes land use change. The negative finding is that managed forest areas occur in the last remaining forest fragments outside of protected areas, which challenges the continuity of selective logging in the near and mid future due to the potential impoverishment of the isolated forest stands. The most negative aspect is the recent establishment, by law, of a very short return logging cycle (10–15 years—minimum cutting cycle) as it contradicts ecological studies showing that tropical forests require over 60 years for the forest structure and composition to return to a state similar to pre-logging states. Our results should inform decision makers and managers in conservation areas to review current policies and establish new benchmarks for sustainable forest management in the country.


Introduction
Worldwide, the sustainable use of tropical forest areas outside of protected area networks (e.g., national parks, wildlife refuges, biosphere reserves) is the only option for preserving a wide range of species and ecosystems absent from protected areas [1].Therefore, tropical forests outside protected areas deserve special attention in terms of their management, anthropogenic dynamics and conservation.As compared to human activities (e.g., cash crops, mining, etc.) that eliminate forest ecosystems and their well-recognized services (scenic beauty, mitigation of greenhouse gasses, watershed protection, biodiversity), it has been argued that natural forest management (i.e., conventional logging, reduced impact logging) has the potential to be a land use that preserves forests in the long term [2,3] and has a low impact on species richness and abundance [4,5].At the same time, other studies have shown the negative impacts of forest management in terms of their immediate and mid-term effects at different spatial scales, such as changes in the hydrological patterns of streams, forest impoverishment and increases of forest edges [6][7][8].Nevertheless, the impact of natural forest management varies depending on the nature and extent of forest damage [9] and the policies that provide both a legal framework and its enforcement.The enforcement should allow forest stands to be preserved after the extraction of timber resources in the long term.
The estimated size of natural tropical Permanent Forest Estates (PFE: forests under state ownership and centralized control) worldwide is 761 million hectares, comprising 403 million hectares of production forest and 358 million hectares of protected forest [10].Between 2005 and 2010, the area of natural forest under management plans in International Tropical Timber Organization (ITTO) producer countries increased to 183 million hectares, which accounts for 24% of the PFE.Worldwide, certified natural-forest production in PFEs increased 63% and sustainable forest management (SFM) doubled during the same period [10].SFM is important for forest conservation purposes, since it is based on a series of principles, criteria and indicators that aim to maintain, over the long-term, the biological, physical and social characteristics of these ecosystems [11,12].
In Costa Rica, the management of natural forests has been regulated in different ways since the establishment of a forestry legislation framework.For instance, the first Forest Law (1969) did not include any restrictions on selective logging.Later, Laws 7032 and 7174 (1986 and 1990, respectively) required a registered professional to produce a management plan prior to forest exploitation [13].The 1996 Forestry Law 7575 established the Forest Manager with public's faith who is overseen by a Union of Professionals.The Forest Manager is in charge of designing and supervising forest operations [13] and most importantly, prohibited land use conversion from forest to non-forest.Public's faith means that the statements of fact made by a Forest Manager in a written instrument are considered as evidence of the facts recorded and that those statements are irrefutable unless a litigating party brings a special plea alleging fraud or forgery (Valverde, 2014 personal communication).
The 1996 Forestry Law 7575 enhanced forest protection and led to the development of technical procedures to facilitate selective logging in forest stands.These technical aspects, among others, include the collection of land tenure information, a forest inventory and census, technical specifications about logging roads, and topographic maps indicating stream networks and large tree locations (>60 cm DBH).All of this information is integrated into a technical document called the natural forest management plan produced by the Forest Manager and submitted, for approval, to a regional office of the Forestry Authority of the State.In this case a natural forest is defined by law as a native ecosystem that encompasses a minimum area of two hectares and is characterized by the presence of mature trees of different ages, species and forms, one or more canopy layers that cover more than 70% of the surface, and more than 60 trees per hectare with a DBH greater than 15 cm.In this paper, natural forest management refers to the selective logging practices for the production of timber in this natural ecosystem.Further definitions regarding natural forest management in Costa Rica can be found in [14].
Despite these policy efforts to regulate the forestry activity, not a single study has rigorously measured the sustainability of forest management in Costa Rica at the country level and from a historical perspective since the establishment of the law.For instance, information about historical natural management trends are almost unknown or only collected for certain periods of time or regions [15].This problem is not necessarily unique to Costa Rica.With the establishment of the UN REDD+ initiative as a major international policy to prevent forest degradation and deforestation in tropical countries as well as provide funding for sustainable management [16], many tropical countries need to provide detailed information about their forest resources (e.g., inventories) and strategies for their conservation.
In this study, we present an in-depth analysis of natural forest management practices in Costa Rica based on a new historical forestry GIS database encompassing five Conservation Areas in the country where selective logging has been prevalent for 18 years [14,15].A Conservation Area is a geopolitical division used in Costa Rica to regulate different activities related to the sustainable use of natural resources (timber, geothermal energy, water, etc.).Our main objective is to critically assess the overall impact of selective logging practices in different regions of the country.Our temporal assessment is based on an analysis of the historical trends of forest management, while our spatial analysis uses two approaches.The first approach aims to examine whether areas subject to forest management preserve the forest stand in the longer term (i.e., forest "fate") since, according to the Forestry Law 7575, land use change from forest to non-forest is prohibited.Therefore, the expectation is that forest areas subject to logging retain their forest cover.The second approach aims to determine the relationship between forest management areas, landscape structure and composition [17] in each region.
A legislation recently developed in Costa Rica defines the sustainability standards for natural forest management and the Code of Practices including the procedures for polycyclic forest management and the potential for a second harvest rotation of 15 years with a minimum of 10 years [18].We argue that this rotation time might not be adequate for forest recovery since it ignores the current research in forest growth under new climate change regimes [19,20], the changes in tree species composition after logging and the structural recovery of forest stands over time (depending on logging intensity) [21].We compare the modeled distribution of the numbers of tree harvested per hectare for each Conservation Area with the simulated distribution using low felling intensity (LFI) reduced-impact logging (RIL) principles (i.e., 3 ± 1 trees/ha).Based on these results, we provide a general discussion on the potential negative effects of such a short rotation time.

Experimental Section
Our analysis is based on a recently created historical forestry relational database (FRD) and Forestry Geographic Information System (FGIS) encompassing more than 600 selective logging records in five Conservation Areas (CAs) in Costa Rica, spanning more than 20 years [14].The CAs are: Huetar Norte Conservation Area (ACAHN), Coordillera Volcanica Central Conservation Area (ACCVC), Tortuguero Conservation Area (ACTO), La Amistad-Caribe Conservation Area (ACLA-C) and Osa Conservation Area (ACOSA).Conservation Areas in Costa Rica are under the jurisdiction of the Ministry of Environment and Energy (MINAE) through the National System of Conservation Areas (SINAC).Costa Rica is divided into 11 Conservation Areas, which are governed using the same administrative and development model concerning the management of natural resources in both public and private lands [13].Figure 1 shows the Conservation Areas included in this study as well as the number of management plans (n) used to calculate historical trends and spatial statistics.To analyze the historical trends of forest management in each conservation area we queried the FRD for each management plan that was presented by a Forest Manager in the regional office of the specific conservation area.We then calculated basic statistics (mean, median, standard deviation, and coefficient of variation) for total property area, and the area of the productive and protected zones of the forest stand (Table 1).Protected zones in the management unit are defined as areas near water bodies (e.g., 15 to 50 mt buffer around streams and 100 mts buffer around springs) and on steep slopes (>40%) where selective logging is prohibited, while forest areas outside the protected ones are considered productive zones.A more detailed description of natural forest management practices in Costa Rica can be found in [14].The total area of the property may also include other land uses (e.g., agriculture, pastures, etc.).
To assess the impact of forest management on the long-term permanence of forest cover (over the last 18 years), four datasets were used: the FGIS, the forest cover assessment for Costa Rica for the year 2005 [22], recent images in Google Earth (2012-2013), and an aerial color photo mosaic at 30 cm spatial resolution (2005).The reason for the use different sources was to ensure that the interpretation of forest cover was consistent between datasets.Moreover, for management units logged after 2005 the aerial photo mosaic and the forest cover assessment were inadequate.We assessed forest permanence for 422 management units, beginning in 1989.We only considered management plans where we were able to record the logging report in our FRD, therefore, we were certain that selective logging operations occurred.Since the FRD was based on the recovery of archival information in different conditions (e.g., incomplete documents), there may be management units that do not have a logging report but actually were harvested.Coordillera Volcanica Central Conservation Area (ACCVC), Tortuguero Conservation Area (ACTO), La Amistad-Caribe Conservation Area (ACLA-C) and Osa Conservation Area (ACOSA); the darker areas in the map represent the forest cover for the year 2005.In general, forest management plans are located outside protected areas such as national parks and wildlife refuges (not shown here to avoid map clutter) in fragmented forests; the source for the original data is the National System of Conservation Areas archives of the regional offices, however, data presented here is based on work developed by [14,15]; areas for ACTO were provided by the regional office and rectified when necessary.For each forest management unit in the FGIS we overlaid its polygon with the three datasets and visually analyzed the forest stand.Because the property units were digitized from the original forest management plans [15], geo-positional errors needed to be visually verified and manually corrected using ArcGIS 10.1 [23].To improve geographical accuracy we revised the digitized units against a geocorrected cartographic dataset (1:50,000) for the country in the Costa Rica Transverse Mercator 2005 coordinate system.For each polygon, a PDF document of the original management plan was also used to compare the location of the forest as illustrated by the forest manager, when the full description of the land use/land cover was given.In cases that this map was not available, tabular data about forest cover area (in hectares) was used.Our exhaustive analysis aimed to avoid overestimating or underestimating the original forest cover included in the productive and protected zones within the management unit (total property area).As a last step for the forest permanence analysis, three categories were created to summarize our findings: agreement with the original plan, partial agreement with the original plan and no agreement with the original plan.The category of agreement with the original plan signifies that the forest unit encompasses the full forest extent originally mapped by the management plan, regardless of the area of the forest.A partial agreement indicates a reduction of the original forest cover, which suggests a reduction of the initial forest extent.Lastly, the category of no agreement with the original management plan, specifies a land-use change (e.g., from forest to cash crops such pineapple, bananas, etc.) (Figure 2).Because land-use change was prohibited after 1996, we paid special attention to cases that resulted in a 'no match' in terms of geo-location in order to avoid false-positives.Geographic errors inherent in the coordinate transformation for Costa Rica and the on-screen digitization of archival documents were triple-checked with different datasets and distinct features on the ground when available (e.g., roads, rivers, etc.); this example is an agreement with the original management plan because the forest is still present in the management unit.
The second aspect of the spatial analysis was to assess the relationship between forest management location and landscape forest fragmentation.After aggressive deforestation that lasted until the mid-1980s [24], Costa Rican forest ecosystems outside of protected areas became highly threatened and fragmented [25,26].It is therefore key to understand not only where forest management occurred, but also its linkage to forest spatial and landscape features (e.g., patch size, connectivity).Selective logging in fragmented landscapes is of special consideration because little is known about its effect at Color ae rial pho to local and landscape levels for different biological groups [27,28].As a first step to understand the location of selective logging units and the fragmentation patterns within each Conservation Area, we calculated forest fragmentation statistics (mean patch size, standard deviation, median patch size) using Fragstats 4.2.Further, we carried out a cluster analysis using a Multi-Distance Spatial Cluster Analysis (Ripley's K Function) for the forest management units using ArcGIS 10.1.For the fragmentation analysis, we used the forest cover assessment for the year 2005 [22].Because the large fragments include protected areas, we calculated the forested area not under protection (e.g., national parks, wildlife refuges) to determine the potential effective area for forest management and the location of management plans.A first step for polycyclic management under the resolution is revising whether there is a previous management plan document (hard copy) in the regional office of the Forestry Authority of the State.Since there are inconsistencies in the historical records of the management plans in Costa Rica [15], the legislation provides a framework for management units with or without a previous management plan document.Under an existing record, the previous information serves as a baseline of comparison for the new plan and the development of an annual operational logging plan.Without a previous record, the forest manager has to follow the RIL guidelines stated in the forestry law resolution [18].
To conduct a preliminary analysis of the potential effect of forest logging and the 10-15 year second return established by the Code of Practices for polycyclic management, we calculated the probability distribution of the number of selectively logged trees per hectare extracted from the approved management plans.Then, we determined the best-fit probability density function, based on the AIC criterion, for each Conservation Area using Matlab 2014a.We compared the results of the density functions to that of a theoretical one modeled for low tree felling intensity (3-4 trees/hectare) assuming RIL operations [29].Reduced-impact logging encompasses a series of controlled timber harvesting practices that aim to reduce the well-recognized deleterious effects of selective logging in the forest stand [30].However, even under RIL with an average extraction of 6 trees/ha there still might be considerable damage to the canopy [29].Although RIL (low or medium intensity) is not a common practice in Costa Rica, it provides us with a baseline against which we can compare the more conventional logging practices used in the country.We use the number of logged trees per hectare as a proxy for forest disturbance, under the assumption that a larger number of trees extracted per hectare using conventional logging techniques would increase both the collateral damage to the forest stand (e.g., tree felling, roads) and the time for vegetative compositional and structural recovery [21,31].For instance, the minimum logging diameter in Costa Rica is greater than 60 cm in DBH, and the cutting intensity is a 60/40 ratio with 60% of the commercial volume available for harvesting and 40% left as remnants.We discuss the implication of current forest practices and the short time frame allocated for forest recovery from logging operations in Costa Rica.

Historical Forest Management Trends
The average managed extent (ha) for selective logging units ranged from 51 hectares in ACTO to 37.1 hectares in ACLA-C, and is similar between the five conservation areas (Table 1).In comparison to other tropical countries with large tracks of continuous forest [32,33], forest areas in Costa Rica are small and highly variable in terms of size and topographic features.The total managed area encompasses a productive area where selective logging operations occur and protected zones, where selective logging is prohibited based on legal regulations (Forestry Law 7575).Our results indicate that protected zones are low in comparison to the area undergoing logging.The protected zones inside the management unit are the areas with steep slopes and within the buffer zones around water bodies.For instance, in ACAHN, the average productive area is 35.6 ha and the average protected zone is 6.3 ha while, in ACCVC, the protected zone represents 8.7 ha of the total managed area 28.1 ha.In the case of ACOSA, the average productive area is 21.3 ha and the average protected zone is 19.2 ha (Table 1).In both ACAHN and ACOSA, however, the standard deviation of the productive area is notably larger (24.3 ha and 22.0 ha, respectively) (Table 1) than within the other Conservation Areas.In the case of ACOSA, most of the forest management is carried out in areas of moderate to steep slopes [34], which may explain this difference compared to other conservation areas, where forests located in productive areas is larger.
Our trend analysis indicates a concentration of forest management activity between 1995 and 1999 in all Conservation Areas with a decline toward the early 2000s (Figure 2).This trend is even more noticeable when the total numbers of management plans in all Conservation Areas are summed (Figure 2).This result is consistent with the findings of [15], indicating the potential effect of the Forest Certificate of Payment for Natural Forest Management (known as CAFMA) established in 1994.The CAFMA program was established as an economic subsidy to promote forest management by supporting the preparation of forest management plans and the implementation of silvicultural treatments in natural forests.Nonetheless, this incentive started decreasing in 1999, country-wide, due to the negative effects of the forest management model on biodiversity as perceived by different groups (e.g., environmentalists and biologists) and the lack of auditing and control for approved management plans.For example, NGOs pressured the government for increased protection and reforestation, as forest management was perceived as "hidden deforestation" [35].With a large number of management plans presented to the regional offices between 1996 to 1998, and a shortage of technical personnel to review the management plans in situ, many inconsistencies were found in the management plans that potentially negatively affected logging operations.For instance [36] found that in ACOSA 74.4% of the management plans had anomalies, such as incomplete documentation or a lack of rigor in the statistical analysis for the number of trees, the basal area, the volume and the sampling design of the inventory.In addition, poor mapping of tree locations was an issue.Parallel to the decline in CAFMA, the Payment for Environmental Services (PES) program was established in 1997 and financed private farmers' conservation efforts.The PES program was widely apply in the country between 1998 and 2005, accounting for 91% of the area covered since 1998, and for 95% of the enrolled area by the end of 2005 [37].As [38] states, subsidies to the forestry sector were politically unsustainable because of insufficient economic contributions from the forestry sector.Therefore, the PES program fulfilled the need in the country for protecting forest areas, while providing forest owners with an economic tool for conservation.

The Good: An Assessment of Forest Fate
Our assessment of forest fate indicates an agreement or partial agreement of the original forest extent in the management unit (i.e., forest still standing) for more than 90% of the forest units.The highest agreement in ACAHN where 90% of the forest management units retained their forest cover, while the lowest is in ACCVC and ACTO (Table 2).Three Conservation Areas (ACTO, ACLA-C, ACOSA) show a low percentage of forest conversion, 8% or less.Overall, our forest fate results indicate that most of the forest subject to selective logging maintains a forest cover after logging operations.In some instances, logging operations might have a negative impact on the ecosystem.For example, in a two million square kilometer area of the Brazilian Amazon [39] found that logging operations resulted in severe canopy damage and may have promoted deforestation (16% ± 1% of selectively logged areas were deforested within one year of logging and an annual deforestation rate of 5.4% for four years after timber harvests was found).Nevertheless, [40] found that natural forest concessions (parcels of natural forest leased out to companies to extract timber on a long term basis (>30 years)) in Kalimantan (Indonesian Borneo), have a similar impact on slowing forest cover loss than those provided by protected areas in comparison to expanding plantations.It is important to mention that in Costa Rica, natural forest management is small scale (Table 1) and is generally carried out by the owner of the property who contracts the services for timber extraction (from planning to harvest).

Table 2. Assessment of forest permanence after selective logging in five Conservation
Areas in Costa Rica; the category agreement with original plan corresponds to a management plan that kept its original forest cover (aerial view analysis) after selective logging; a partial agreement with original plan corresponds to units that might have lost some forest area after logging (qualitative assessment), while no agreement with original plan are areas with no forest cover, meaning that the land use/land cover changed from forest to non-forest.Our analysis does not take into account the forest stand quality after logging (e.g., changes in vegetation species richness and structure); values are given in percentage of the total number of management plans used for the analysis and n corresponds to the number of management plans used for the analysis.A notable aspect of our analysis is that it only assesses the presence of forest cover and does not evaluate the quality of the forest stand after logging and its recovery through time.For instance, depending on the disturbance created by felling gaps and extraction roads, a long time may be needed for the recovery of the forest structure and tree species composition [21,41].Moreover, studies show a high variability of responses to logging practices depending on the objective of the study (e.g., species chosen, spatial and/or temporal scales), the sampling design and the ecosystem [42].

The Bad: Relationship between Forest Management Location and Forest Fragmentation
Fragmentation statistics for the different conservation areas show a larger number of forest fragments in ACAHN (2241) and ACTO (1144), but these forest areas collectively represent a lower percentage of the land (24% and 47%, respectively) (Table 3).An indication of continuous forest areas is given by the large patch index (LPI) and the percentage of the total forest area this patch represents.In ACLA-C, large forest fragments encompass 69.7% of the total forest area, whereas in ACAHN, the large forest fragments only account for 5.4% of the total forest area (Table 3).When evaluating the average and variability of fragment size, the mean patch area and standard deviation follow the same pattern as the LPI, indicating a high variability in the landscape across Conservation Areas.In terms of the geometric characteristics of the fragments in the landscape, the shape index average indicates very irregular forms of the forest fragments, with values ranging from 1.81 (ACAHN) to 2.08 (ACCVC) (Table 3).These patterns provide a clear overview of the landscape characteristics in these Conservation Areas.Our results are in agreement with previous studies that indicate a high degree of forest fragmentation in Costa Rica as a result of the rapid advance of the agricultural frontier during the 1970s and 1980s [13] and continued deforestation occurring at smaller scales in the last 15 years [43].A more detailed depiction of the landscape in the different Conservation Areas is shown in Figure 3, which illustrates the distribution of the patch class areas as well as the percentage of total forest area of each class in the landscape.Although there are numerous small forest fragments in the landscape (3-10 ha), they represent less than 10% of the total forest area in all of the Conservation Areas.In addition, most of the large fragments (>5000 ha) contribute to a higher percentage of the forest area (>50%).Our results, not including the protected areas (e.g., national parks), show similar trends for all Conservation Areas, with a few large forest fragments accounting for a higher proportion (50%-80%) of the total forest area.Figure 3. Historical trends of forest management plans presented to the regional offices for five Conservation Areas in Costa Rica and the total number of records across Conservation Areas; the sources for the original data are archives of the regional offices of the National System of Conservation Areas, however, data presented here was is following processed and analyzed by [14,15].

Our cluster analysis results indicate a high degree of spatial clustering for management plans in all
Conservation Areas as illustrated in Figures 4 and 5.A visual analysis for each Conservation Area indicated that forest management plans are located and clustered around the larger patches in the landscape, although there are many fragments of less than 100 ha scattered throughout each Conservation Area.Given that our average forest management unit is between 37 and 51 hectares, we do not have a specific explanation for the lack of selective logging in smaller patches (<100 ha) across the landscape, but it may be that these fragments do not provide the necessary volume and species required for logging operations.There is a substantial amount of evidence that smaller fragments in tropical areas have low species richness and are more affected by edge effects [44].Moreover, in a study in southeastern Brazil, [45] found that fragment connectivity was highly correlated with the species richness of shade-tolerant species, with fragments of low connectivity containing a lower species richness of these shade-tolerant species.A second potential explanation is that many of these forest areas belong to small to medium landholders who might have chosen the Payments for Environmental Services program as a way to get an economic incentive for preserving their forest areas [46].Regardless of the variable(s) that explains this pattern, these small forest fragments represent the last forest remnants in the agricultural matrix and their conservation value is significant for many services as well as for landscape connectivity.

The Ugly: a Preliminary Discussion on Polycyclic Management in Costa Rica
Results indicate a higher probability of management plans extracting timber in the 5-7 trees/ha bins for the Conservation Areas, which is within the range of medium felling intensity conventional logging (Figure 6 A-E) [29].For example, for ACOSA the highest probability for felling intensity is 5 trees/ha with 56% of the management plans falling into the 5 to 11 trees/ha potentially logged range (Figure 6F).In ACTO, the highest probability for felling intensity is 4 trees/ha but 72% of the management plans fall into the 4 to 24 trees/ha potentially logged range.Figure 6F clearly identifies the differences between the fitted function of trees per hectare potentially logged for each Conservation Area in comparison to a low intensity reduced-impact logging scenario, centered on 3 ± 1 trees per hectare [31].These results also indicate that heavy intensity logging (≥6 trees/ha) may be present in all Conservation Areas to some extent (Figure 6).In the two Conservation Areas with the highest intensity felling, 38% (ACTO) and 42% (ACOSA) of the management plans potentially extracted six or more trees per hectare.Under a low felling intensity RIL scenario 48% of the management plans would extract between one and three trees per hectare.For the Conservation Areas only 5% (ACOSA) to 29% (ACAHN) of the potentially extracted trees fall under the low felling intensity scenario.It is important to note that RIL was not practiced in any of the Conservation Areas and we mention it here as a qualitative approximation of sustainable forest management [30].In this case, we use numbers of trees per hectare as a proxy for forest disturbance based on the canopy damage that felling trees could cause and the gaps they create in the forest stand [21,39].Given that forest logging operations focus on large adult trees (≥60 cm DBH) and are based on the extraction of 60% of the total commercial volume allowed by Costa Rican law in natural forest extractions (for all species) [47], it is likely that forest stands are highly affected by logging disturbance [21].In this study, we do not take into account damage created by the felling of non-target trees and any additional disturbances created by skid roads and timber patios, where heavy machinery (e.g., D5 tractors) is commonly used for forest extraction, primary and secondary road construction [30].As noticed by [29], most logging operations occurring in the tropical rainforest remain unplanned and are very destructive to the stand.Our preliminary assessment indicates Costa Rican forest extraction operations may not be and exception to this finding (Figure 6).Further analyses of our data will focus on analyzing information from the management plans to measure the total disturbance reported.Nevertheless, field evaluation and/or remote sensing might be necessary to have a more complete understanding of the damaging effects of logging practices in Costa Rican forests.The current Forestry Law and its Code of Practices for polycyclic management indicate a possible rotation cycle of 10-15 years after the first logging extraction.To the extent of our knowledge, there is no ecological, economic or social study that justifies this rotation cycle for the recovery of the forest stand, given the historical extraction practices and the potential damage caused by the conventional logging used in the country (Figure 6) (however, see [48]).As stated by Bonnell, et al. [41], the uncertainty in recovery times after selective logging would result in the establishment of inappropriate harvesting rates and/or a failure to achieve conservation targets.For example, in Uganda's Kibale National Park, Bonnell, et al. [41] found that the recovery of basal area will take between 97 years, for moderately logged compartments, and 112 years, for heavily logged compartments.A plot level study in the eastern Brazilian Amazon indicates that, under conventional logging (5.6 trees per hectare) forest recovered only 77% of its original above-ground biomass stock 16 years after logging, and the biomass of all size classes decreased when compared to pre-logging levels [49].More importantly, for large diameter individuals, trees >60 cm DBH only contained 48% of their previous biomass.At the regional level, a simulation study using geographical modelling found that, for a large area in the Brazilian Amazon, the loss of forest carbon due to selective logging lasts two to three decades following harvest, and the original live biomass would only recover after 100 years [50] without subsequent disturbances of the forest stand.As noticed by [51], in a meta-analysis of selective logging in tropical regions the decreases in timber yields after the first harvest from old-growth forests seem inevitable.For selectively logged forests to fully regain the volume of timber accumulated during the preceding centuries while using a 20-40 year harvest cycle, forest management practices would need substantial modification.This, we believe, is the case for Costa Rican forest management practices.A substantial modification to the new Code of Practices for polycyclic management needs to be implemented if forest management is to be a sustainable land use in future decades.
The Code of Practices for polycyclic management provides a framework with very specific guidelines for selective logging in areas with or without previous natural forest management plan records.For example, with the advancement of technology (e.g., GPS), there has been an effort to establish better geo-location of the management units, sampling plots, adult trees (>60 cm DBH) and protected zones.Previous management plans lack a consistent quality in their mapping efforts [15] and also showed bias towards higher biomass zones in the management unit [52].This combination of factors potentially had an effect in the approval of the management plans by the Forestry Authority of the State, especially during the period of 1996-1998 where a large number of logging permits were requested and approved by the regional offices (Figure 3).
As a way to better understand part of the ecological processes in forest regeneration [53], the Code of Practices encompasses a classification based on the tolerance of vegetation to light regimes (short lived-sciophytes, long live sciophytes and heliophytes).This classification is therefore used as a proxy for understanding the densities (#trees/ha) of timber species (light hardwoods to hardwoods) and the economic value of these species in the forest stand [18].This classification in combination with reference thresholds [54] adapted for each Conservation Area and their ecosystem characteristics would allow for specific management regimes in the forest stands and a better definition of logging intensities.An example of this threshold is the density of short-lived sciophytes, where they should not surpass 15% of the total number of trees >10 cm DBH in the forest stand, because it most likely indicates a less productive forest stand that is still under recovery.
For the most part, we believe the methodology presented in the new Code of Practices for polycyclic management uses ecological principles in combination with an already well-standardized methodology for collecting the biophysical characteristics of the forest unit [14].Still, we argue that assuming that 10-15 years is enough for a forest to fully recover its structure and tree composition given the results obtained here, in terms of the number of trees logged per hectare under conventional logging (Figure 6F) and the very fragmented landscapes, might be insufficient for the long-term conservation of these complex forest ecosystems.Finally, we hope our results are taken as an opportunity to further research in forest management recovery after selective logging [21] in the country and other tropical regions with similar characteristics, as well as a chance for decision makers and stakeholders (e.g., NGOs, scientists, etc.) to understand natural forest management issues at a broader landscape scale.

Conclusions
Our study takes advantage of a fully developed and novel forestry geographic information system and relational database based on historical natural forest management data from Costa Rica.We provide a first analysis of the historical patterns of forest management in five Conservation Areas where selective logging has been a common land use.Forest management has very likely been influenced by policies in the country such as the incentives for forest management during the mid-1990s and a shift to the conservation of forest stands under the umbrella of the environmental services payment program.This shows the importance of the implementation of conservation policies in the country to forest management.On a positive note, it could be argued that forest management practices in the mid-1990s period might have helped to preserve the forest stands around Conservation Areas where logging was carried out, especially during the transition period with the law banning landuse changes (forest to non-forest) around 1996.A less positive aspect of forest management is the pressure of forest management around protected areas; we found a highly clustered pattern of management plans in these regions.This might be a consequence of the highly fragmented landscape matrix outside of continuous forest and the understanding that extremely patchy forest fragments may not be adequate for selective logging practices.Finally, we open a discussion on the rather ad hoc logging cycle of 10-15 years proposed by the Forestry Law, despite significant evidence of much longer forest recovery times in tropical regional undergoing conventional logging.Regardless of this negative finding, we recommend that sustainable forest management regulations and the implementation of reduced impact logging should be considered for reducing forest degradation in the country's fragile and isolated ecosystems.This would allow the country to strengthen its REDD+ initiative and provide a better framework to understand forest practices at local and regional levels.Future research must focus on how to implement the UN-REDD+ program for establishing better sustainable forest management practices and develop a most needed research agenda to measure the carbon stocks and fluxes in selectively logged forests in Costa Rica.

Figure 1 .
Figure 1.Location of forest management plans in Costa Rica (black dots) for five Conservation Areas based on the Forestry GIS: Huetar Norte Conservation Area (ACAHN), Coordillera Volcanica Central Conservation Area (ACCVC), Tortuguero Conservation Area (ACTO), La Amistad-Caribe Conservation Area (ACLA-C) and Osa Conservation Area (ACOSA); the darker areas in the map represent the forest cover for the year 2005.In general, forest management plans are located outside protected areas such as national parks and wildlife refuges (not shown here to avoid map clutter) in fragmented forests; the source for the original data is the National System of Conservation Areas archives of the regional offices, however, data presented here is based on work developed by[14,15]; areas for ACTO were provided by the regional office and rectified when necessary.

Figure 2 .
Figure 2. Process for assessing the permanence of the forest stand after selective logging.Geographic errors inherent in the coordinate transformation for Costa Rica and the on-screen digitization of archival documents were triple-checked with different datasets and distinct features on the ground when available (e.g., roads, rivers, etc.); this example is an agreement with the original management plan because the forest is still present in the management unit.

Figure 4 .
Figure 4. Number of forest fragments (left y axis) and percentage of the total % forest area (right y axis) with and without protected areas (e.g., National Parks) for five Conservation Areas in Costa Rica; the class area (ha) represents intervals for different forest fragment sizes encompassing small forest fragments to very large continuous forest areas (5000+ has).

Figure 5 .
Figure 5. Multi-Distance Spatial Cluster Analysis (Ripley's K Function) of forest management units for ACAHN (Northern Costa Rica) and ACOSA (Southern Pacific); the other three conservation areas also indicate a high degree of clustering at all spatial scales (Data not shown).

Figure 6 .
Figure 6.Probability distributions based on trees per hectare proposed to be logged for management plans in five Conservation Areas in Costa Rica (A-E); light colors on the bars represent low logging densities (trees/ha); more intense gray color indicates conventional logging practices and dark gray bars very heavy logging; (F) shows a fitted probability density function for each Conservation Area compared to a low density logging fit (3 trees/ha).

Table 1 .
Descriptive statistics calculated from management plans for five Conservation Areas in Costa Rica; the total property area uses the cadastral map information, while the managed area (forest stand), productive area and protected zone is defined by the Forest Manager.

Table 3 .
Fragmentation statistics for five Conservation Areas in Costa Rica based on the Forest Cover Assessment 2005.