Carbon storage estimation in a secondary tropical forest at CIEE Sustainability Center, Monteverde, Costa Rica.

Secondary growth tropical rainforests have the potential to sequester large amounts of atmospheric carbon dioxide and as such are an important carbon sink. To evaluate a local forest, a Carbon Neutrality Program was initiated at the Council on International Educational Exchange, San Luis Campus, Monteverde, Costa Rica. The study was conducted on 50 hectares of forest classified as Premontane Wet Forest. The forest, part of the Arenal-Monteverde Protected Zone, is estimated to be aproximately 50 years old and is in the upper regions of the San Luis valley at 1100 m elevation. Assessment of the carbon stock in trees was carried out in two permanent, 1 hectare plots, 100 m by 100 m, Camino Real and Zapote. The plots were divided into 25 subplots, 20 m by 20 m totaling 400 m2 per subplot. Ten subplots in each area were studied which represented 1.6% the total surface area of the forest. All of the trees were measured within the subplots that had a diameter at breast height ≥ 10 cm and the height of 10% of the trees measured. The estimated total CO2 sequestered by the campus forest was 18,210 ton (in 2019).

Introduction

Carbon dioxide is recognized as the primary greenhouse gas of the four gases that contribute the most to global warming. They include, carbon dioxide, CO2 (81%); methane, CH4 (10%); nitrous oxide, N2O (7%); and, halogen containing gases, CFCs (chlorofluorocarbons) (3%)[1-3]. Carbon dioxide results primarily from the combustion of fossil fuels which are derived from plants. Coal originates from decayed plant material occurring in swampy environments and is buried over geologic time. Oil originates from the phytoplankton and other organic matter that accumulates in oceanic sediments and after burial migrates underground to reservoirs. Coming full circle, Dyson in 1977 was the first to propose removing CO2 from the atmosphere using plant photosynthesis[4]. The present study further explores the importance of carbon capture and sequestration in a neotropical rainforest.

The concentration of carbon dioxide in the atmosphere has increased from preindustrial concentrations of approximately 280 ppm (parts per million by volume)[5] to the present 419 ppm worldwide[3,6-8]. As atmospheric CO2 has become a central focus in climate change science the study of all aspects of the biogeochemical carbon cycle is the focus of much research[1]. One of the most important aspects of the carbon cyle, is the estimation of the halflife of CO2 in the atmosphere[9-13]. However, because of the complexity of the carbon cycle, the estimates range from 30 to 95 years[14]. If one considers the long CO2 concentration tail in some models, it has been suggested that the halflife could be as long as 300 years[15]. Better constraining these estimates is a very active area of research and it is possible that the Climate Model Intercomparison Project now in its sixth iteration (CMIP6) will lead to better estimates of atmospheric CO2 lifetimes[16,17].

Limiting climate change will require that the atomspheric concentration of CO2 is stabilized[18,19]. The IPCC (2018) report entitled “Global Warming of 1.5 °C” emphasizes that new or renewed efforts are needed to reach the goal of stabilizing atmospheric concentrations of CO2[20]. Studies have been published addressing stabilizing atmospheric CO2 concentrations[21-24]. However, the concentration of CO2 continues to rise a little over 2 ppm per year and forests are likely to play an important role in the future[5].

The importance of forests was reviewed in a report prepared by the World Resources Institute[25], and more recently by the IPCC[26,27]. Forests (or planting forests) can be thought of as reverse engineering human deforestation and providing a natural sink for CO2[28-30]. It is estimated world-wide, that of the 4.06 billion hectares of forest, 45% (1.8 billion hectares) are tropical forests[31,32]. Tropical forests provide the highest potential for carbon capture and sequestration because they have the highest carbon density of all forests[33].

Many institutions are studying their environmental impact by initiating carbon neutral programs. To reduce their carbon emissions and offset those emissions which they can not reduce, a potential first step is the preservation of forests on-premises. In this approach, growing trees sequesters carbon. The objective of this study was to calculate the carbon stored by trees in 50 hectares of forests by studying in detail two one-hectare plots at the Council on International Educational Exchange (CIEE) San Luis campus, Monteverde, Costa Rica.

Results

The focus of this study was the above ground carbon captured. This was designed to define the carbon capture and sequestration of a secondary forest at one point in time. Our data, summarized in Fig. 1, shows that both plots have predominantly young tree growth with the diameter at breast height (DBH) between 10 and 30 cm and the distribution of tree sizes is approximatley the same. This confirms their young age and secondary forest status. However, Zapote (the top line) appears to be somewhat older and had a higher number of trees, when compared to Camino Real, at every DBH.

Figure 1

The distribution of trees based on diameter at breast height (DBH) for the Camino Real and Zapote plots.

A frequency distribution suggested that the four data points of DBH > 100 cm did not fit the data and suggested that these four trees were older that the rest of the plots and were considered outliers. Therefore, they were not included in the estimation of basal area for either of the plots.

To find the carbon storage within the plots, it was necessary to calculate the basal area of the subplots, which referred to the total surface area covered by tree trunks in relation to the surface area of a plot. A statistical analysis for “tree basal area” was performed to give the average basal areas per plot. The standard deviation, standard error, coefficient of variation, and, sampling mean error in basal areas were no more than 20%, at a confidence interval of 95% using Student´s Distribution Table (n − 1). After calculating the carbon stored in each subplot the value for each hectare was determined and then the result was multiplied by the number of hectares of each forest.

The basal area results for Zapote and Camino Real plots are summarized in Fig. 2, for each subplot map the basal average is given. It is noted that the Zapote (average area 49.3 m2 ha−1) plot had a much greater tree cover than Camino Real (average area 23.7 m2 ha−1). The results of the basal areas of each plot are consistent with the plot assignments given based on the “Google map” tree densities.

Figure 2

Basal areas (m2 ha−1) for the Zapote and Camino Real plots showing subplot locations.

In total there were 60 diferent species of trees in both plots. In the Zapote plots the most common species were Eugenia guatemalensis (21), Dendropanax arboreus (18), Cecropia obtusifolia (18), Inga punctata (28), and Nectandra salicina (38). In Camino Real, there was a substantial number of Inga punctata, in addition to the four above.

As seen in Table 1, based on the National Forest Inventory[34], a secondary forest has approximately 204 tonnes ha−1 of CO2 while a mature forest has 447 tonnes ha−1. The Camino Real plot with 243 tonnes ha−1 was close to that of a secondary forest while the Zapote plot with 518 tonnes ha−1 was closer to the mature forest of 447 tonnes ha−1. Zapote’s high basal area and carbon storage indicates that it is an older secondary forest with many young trees competing with the larger trees.

Table 1

Calculated values for the two permanent plots and the average values reported for a secondary forest and mature forest.

Variable	Zapote plot	Camino real plot	Typical secondary forest (National Forest Inventory of Costa Rica)[34]	Typical mature forest (National Forest Inventory of Costa Rica)[34]
Numbers of trees (ha−1)	737	475	415	469
Average basal area (m2 ha−1)	49.3	23.7	19.3	32.3
Volume (m3 ha−1)	560	269	159	338
CO2 (ton ha−1)	518	243	204	447

Table 2 summarizes the data and shows that the Zapote plot collected slightly over twice the CO2 of the Camino Real plot. This is reasonable because of the differences in the basal areas of the two plots. These two plots are much different from each other which suggests that they were good representatives of the total 50 hectare campus forest area. Taking the average CO2 per hectare of each plot and extrapolating to the campus forests provides an estimate of the total carbon sequestration.The estimation of the total CO2 collected by the campus forest to date is aproximately 18,210 tonnes.

Table 2

Estimation of the carbon storage in the CIEE forest, Monteverde, Costa Rica.

Plot	Average ton of CO2 plot−1	Average ton of CO2 ha−1	Corresponding hectares
Zapote	11,400	518	22
Camino real	6810	243	28
Total tonnes CO2 in the forest	18,210

Limitations

The tags of some trees (< 20 individuals) had fallen off; therefore, they were identified based on past measurements and locations recorded, if they could not be identified, they were marked as new trees. This should not affect the final results of our analysis since we are looking at the total basal area rather than individual trees.

In the Zapote plot, there were no marks on trees which would indicate previous measurements, so the trees were measured according to the Protocol, DBH at 1.3 m from the base unless a tree presented irregularities that required a different measurement location. However, 2018 measurements were taken 10 cm below tags. This likely led to a systematic underestimation of the forest growth, since many trees had tags which were located at 1.3 m.

The accuracy of height measurements was restricted by the tool and our ability to see the top of a tree canopy and the length of the measurement stick. Trees were generally chosen at random, but heights of trees above 40 ft (12.1 m) had to be estimated by eye and were subject to greater error. This could potentially lead to an underestimation of average tree height.

Other problems in assessing different factors; for example, form factor or basic specific weight of each tree, derived from the lack of data on non-commercial trees led to higher uncertainty.

Discussion

The numbers calculated for total carbon stored are typical for secondary forests. The Camino Real area historically has been more interrupted by human activity than the Zapote plot which is why it has less trees, less basal area and ultimately less carbon volume. The Zapote plot, on the other hand, has been less interrupted, so the trees in this permanent plot have had more time to grow and the tree density is lager at each DBH. This is additionally supported by the Google Earth images of both plots, where one can clearly see that the sections labeled as Zapote are denser while Camino Real is more scarce and divided by pastures. There are approximately 518 tonnes ha−1 of CO2 stored in the Zapote plot and 243 tonnes ha−1 of CO2 in Camino Real plot (Table 1).

Secondary forests, like the permanent plots of Zapote and Camino Real, have great potential to sequester carbon through tree growth, as they gain biomass more rapidly than older primary forests. Now that there is a baseline for the volume of carbon stored in these permanent plots, we can calculate the carbon sequestered in subsequent years. This can be done through the same measurement procedure as described in this study, comparing the initial and final carbon storage. This work could be used in the future calculations of carbon sequestration for CIEE’s carbon neutrality program and as a reference for other institutions which also wish to understand carbon sequestration of secondary forests, potentially for their own carbon neutrality calculations.

In the future the sample size should be increased. This study represented 1.6% of the total forested area at the lower end of 1–5% of the total area recommneded for studies of this nature[50].

Methods

Study area

This study was conducted on the CIEE San Luis campus, in the monitoring plots which were first established in 1997 and 2014 by Diana and Milton Lieberman on the grounds of Ecolodge San Luis and Biological Station (University of Georgia, Costa Rica Campus) and purchased by CIEE in 2019. The site has 50 hectares of forest in the upper regions of the San Luis valley at 1100 m above sea level[35]. It is classified as Premontane Wet Forest according to Holdridge and is part of the Arenal-Monteverde Protected Zone[36,37]. The property borders the Monteverde Cloud Forest Reserve and the Children’s Eternal Rainforest and hosts a unique assemblage of flora and fauna from the cloud forest above and dry forest below. In Costa Rica, the Premontane Wet Forest life zone faced rapid deforestation in the mid to late 1900s as vast areas were converted to pasture. As a result, what remains of the forests is second only to the tropical dry forest as the most heavily altered and reduced life zone in the country[38].

The study area is based on two pre-established plots on the campus premises, Camino Real and Zapote. Both plots are 100 m × 100 m (1 ha) divided into 25 subplots each. Camino Real is a plot located alongside the Camino Real trail with its southwest point located at 10° 16.800′ N, 84° 47.952′ W. It was established on a section of secondary forest south of main campus with DBH measurements beginning in 2014. The Zapote plot is on the north side of campus, by the Zapote trail, and has its northwest corner located 10° 17.321′ N, 84° 47.800′ W. Historically, much of the property was used for coffee, guava, banana, and pasture land until the 1970s. However, the area within the Zapote plot was less intensively used and more closely resembles primary forest when compared to the Camino Real plot on the lower slopes of the property near Rio Alondra as shown in[39] Fig. 3.

Figure 3

Map of the campus and location of Zapote and Camino Real Plots. Map of the campus outlined in red, the forest areas corresponding to Zapote in blue and Camino Real in green. (2018 Google, Image ©2019 Maxar Technologies).

Methods

Assessments of biomass, carbon capture and sequestration, is an active area of study and several methods have been published[30,39-49]. In this study the estimation of carbon stock was derived directly by dasometric variables of trees through measuring sample tree attributes in the field, such as diameters (DBH), tree heights and using an allometric equation. All the measurements in this study were conducted in accordance with the Costa Rican National Standard INTE/DN 03:2016 Methodology for Quantification and Reporting Greenhouse Gas Removals Resulting from Forestry Activities (National Standard INTE/DN 03:(2016)[50]. This is the official standard for developing inventories for the Carbon Neutrality Program, and the Protocol for the Establishment and Measurement of Permanent Sampling Plots in Natural Forest[39].

Experimental design

Each one-hectare plot was divided into 25 subplots and then 10 subplots, 20 m by 20 m designated for sampling (Fig. 4). This design represented 40% of the two one-hectare plots studied and 1.6% of the total forest area (20 plots each 0.04 ha/50 ha−1 = 1.6%).

Figure 4

Diagramatic representation of the subplots chosen for this study are the shaded subplots.

Plot reconnaissance

The plots, Zapote and Camino Real, were first set up in 2014. Therefore, prior to going to the field it was necessry to examine previous data collected. This included a list of the trees with DBH ≥ 10 cm, their assigned number, species, previous recorded diameter measurements and tree conditions as well as their coordinates. During the reconnaissance if a tree needed a new tag; for example, if the tag was missing, broken, fallen, or covered by tree growth, a new tag was made and afixed to the tree for future reference.

The northwest corner of the Zapote plot was identified with GPS, then a compass was used to scout subsequent subplots and checked with the permanent plot map. Where possible the visible markers at the corners of the subplots were located. In each subplot, tree tag numbers and referenced the data sheet were used to double check that the tree was in the assigned subplot. Since most tree measurement locations were not marked the tree was numbered, DBH measured (1.3 m from the base of the tree) and marked with paint for future reference.

The southwest corner of the Camino Real plot was identified with GPS and then a compass was used to scout subsequent subplots. In addition, the tree tags in Camino Real indicated the given subplot number (1–25) on the tag which helped to find the correct location of the plot as well. The tree measurement locations in Camino Real were well marked with paint and the DBH was verified.

Tree measurement procedure

The measurements for this study were conducted during the months of June to August 2019. All of the trees were measured that had a DBH ≥ 10 cm and the height of at least 10% of the trees in the plot.

Tree height was determined using a stick measure (telescopic cylinder for heights) to get an overall average tree height in the given subplot. If a tree was over the top of the pole (40 feet, 12.1 m) either the height was estimated or the average height within the subplot was used. For each tree, its condition was recorded according to the conditions specified in Sánchez-Monge[39].

If a tree looked like a new recruit, its DBH was measured to see if it was actually ≥ 10 cm. A new recruit was marked with colored tape with a new tree number, which was assigned as the next number in the subplot based on the datasheet. Then the tree number, diameter, condition, and location based on its proximity to other trees was recorded.

If a tree could not be found based on the coordinates on the datasheet the condition code of 16 (tree not found) would be written down. If a previously measured tree was dead, it would be recorded as dead and not measured.

Calculation of CO2

The following allometric equation, Eq. (1), was used to calculate the above ground carbon stored within the plots:

To estimate the amount of carbon captured, CO2e, a conversion factor of one ton of forest biomass has approximately 0.5 tonnes of carbon, was used.

The basal area is DBH2 * (π/4), and h is the tree height.

However, a tree is not a perfect cylinder so multiple factors have to be applied according to IPCC Guidelines the “fitting factors” are[51]: Finally, 3.67 is the ratio of the weight of carbon in carbon dioxide. This means that one tonne of carbon is equivalent to 3.67 tonnes of CO2 (obtained based on the molecular weights of CO2 and carbon, of 44/12).

FF is the “form fitting” factor which takes into consideration the change in diameter in relation to height. The FF varies according to the tree species but, if the species could not be identified, the IPCC allows a FF of 0.7 in forest trees, meaning the average diameter is 70% of the DBH.

CF is the carbon fraction which represents the estimated volume of carbon stored in the trunk and is considered to be approxcimately 50%, or 0.5.

BSW is the basic specific weight which refers to the weight of the dry wood in a tree. It varies by species and is related to the density of the wood. If the species of the tree is unknown, the IPCC recommends using 0.5[51].

BEFa and BEFg refer to the volume of the branches and leaves (BEFa) and roots (BEFg). For above-ground biomass the estimate for both was 1.2 of the volume of the trunk, calculated by the basal area times height.

Extrapolating to the 50 hectares

To estimate the total amount of carbon captured, a Google Earth image of CIEE Campus was used to estimate how much of the forest was similar to the Zapote, a mature plot, and how much was like Camino Real, an apparently more disturbed plot. Any dense tree areas were defined as more mature, while the areas with more pastureland and scarce growth were defined as younger, corresponding to Zapote and Camino Real, respectively. Figure 3 shows a map of the campus along with mature areas in blue and less mature areas in green. It was concluded that the landscape similar to Zapote, was 22 ha and that similar to Camino Real was 28 ha.