Abstract

Succession following deforestation in Neotropical forests has been investigated extensively, yet rarely have studies connected nutrient dynamics with vegetation. This study was conducted in lowland wet forests of Maquenque, Costa Rica. The objectives were (1) to compare carbon (C), nitrogen (N), and phosphorus (P) characteristics and understory vegetation diversity between regenerating forests and primary forests; and (2) to use these variables to evaluate P. macroloba’s successional role. Four 300 m2 plots were established in primary and secondary forests where P. macroloba was the dominant N-fixing tree. Soil and vegetation data were collected from 2008 to 2010. Values of indicators of C, N, and P cycle activity were generally greater in primary than in secondary forest soils. Efficiency of organic C use and the relative contribution of respiration and organic C to soil biomass were also greater in the primary forest. These trends corresponded with greater richness, biomass, and cover of total and leguminous plant species, greater volume of P. macroloba in primary stands, and greater density of P. macroloba in secondary stands. As cleared regions of former primary forest regenerate, P. macroloba is the important dominant N-fixing tree and a critical driver of C, N, and P recuperation and ecosystem recovery.

1. Introduction

Tropical forests originally covered up to 99% of the land in Costa Rica [1]. Approximately 90% of the original forests in Costa Rica have been destroyed in recent decades, and approximately 46% of the total area has been converted into cattle pasture [1, 2]. In the Northern Zone of Costa Rica, four decades of deforestation have resulted in the loss of about 70% of the lowland forests [27]. The resulting fragmentation has yielded concern among regional scientists regarding whether the remaining primary forests will be able to regenerate at a rate matching the deforestation [7].

Today, tropical forest land is characterized by an expanding proportion of secondary forests [8]. As the area of the secondary forests grows at the cost of primary forests, the secondary forests will have to be managed and used in the future [9]. Rates of recovery for abandoned farmland in tropical areas are accelerated if prior land-use intensity was low, areas recovering are small in size, soils are fertile, and there are remnant forests nearby [10]. However, many areas cleared and used as pasture in Costa Rica have suffered such a severe deterioration of soils that forest successional processes are hindered. Instead of forests returning to the area, researchers are finding colonization of scrub growth and invasive grasses and ferns [1, 11]. Therefore, Costa Rica has a great need for reforestation and assisted natural regeneration [12]. Indeed, attempts to remediate the effects of agricultural activities have resulted in implementation of a variety of restoration strategies, including development of an extensive array of secondary forests [5, 6, 1315].

Soils in secondary forests in the tropics are thought to follow a successional process of development during forest regeneration that parallels the increase in complexity of the vegetation community [1619]. This process results in increases in the abundance and complexity of soil organic matter [2026]. Greater amounts of nitrogen N fixation and ammonium oxidation are also thought to occur in these soils, as the forest develops greater density, diversity, and volume of leguminous vegetation during natural succession or following disturbance [22, 2731].

After tropical forests are cleared, plants with N-fixing root nodule symbionts are considered to be the principal pathway by which these secondary forests recuperate N. These N-fixing species restore soil fertility, stimulating the growth of other plant species and soil biomass development [30, 32, 33]. The increase in abundance and diversity of understory vegetation in forests also results in enhanced leaf litter quantity and quality in forest floor debris and production of more organic carbon (C), N, and phosphorus (P) through increased rates of decomposition and enhancement of the soil organic matter (e.g., [19, 3444]). These activities can stimulate rhizodeposition and increased production of labile root-derived carbohydrates, further enhancing the microbial community [4547]. It is also believed that as forest soils undergo succession, the fungal community becomes more complex [48, 49], including the shift from nitrophobic to nitrophilic fungal groups that conduct nitrification in acidic forest soils [50, 51].

Succession in lowland Neotropical secondary forests has been investigated for almost a century, but more work is needed which connects the below-ground nutrient dynamics with above-ground vegetation characteristics [52]. Research has shown an increase in above-ground biomass, and species diversity generally characterizes older forests, yet much less work has been done to examine below-ground processes that occur during tropical secondary succession [8]. For example, Chazdon et al. [53] discussed the rates of changes in vegetation occurring in secondary forests in Costa Rica following the original disturbance, but there were no connections made to below-ground nutrient cycle dynamics. Eaton et al. [44] showed that 25- and 30-year-old secondary forests in the region had developed a more complex soil microbial community and nutrient structure than the younger 15-year-old forest soils, but there was no connection with above-ground forest vegetation community composition and structure.

In 2001, the Costa Rican government established the 1,204,812 ha San Juan-La Selva Biological Corridor (SJLSBC) to help protect the Northern Zone ecosystems (http://www.cct.or.cr/publicaciones/Estrategia-Exitosa-CBSS.pdf), the core conservation unit of which is the Maquenque National Wildlife Refuge (MNWLR), which conserves the portion of the corridor with the highest percentage of forest cover. The MNWLR contains a patchwork of old growth (primary) forests and harvested forests with natural regrowth secondary forests. Pentaclethra macroloba (Willd.) Kuntze (Fabaceae) is the dominant N-fixing tree in these Costa Rican hardwood forests and is thought to be important in N and C cycle dynamics, biomass enhancement, and forest succession [54, 55] and is presumed to be a fundamental early colonizing tree in these secondary forests. However, to date, there have been no studies that examine the role of this tree in succession of vegetation and nutrient dynamics in these Costa Rican secondary forests. It has recently been discovered that the root nodules of P. macroloba are associated with populations of Frankia spp., while Rhizobium spp. appear to colonize rhizosphere soils somewhat further from the root nodules [56]. The symbiotic association between these N-fixing microorganisms and P. macroloba represents an essential mechanism for N availability in the acidic and infertile soils of the wet tropics of Costa Rica [56, 57].

This study was conducted in primary and secondary lowland wet forests of Costa Rica’s MNWLR. The main objectives were (1) to compare the C, P, and N-cycle characteristics and the diversity and abundance of the leguminous vegetation in these forest stands and (2) to use these variables to begin to evaluate the role of P. macroloba in primary and secondary forest soil and successional processes. We hypothesized that we would find evidence to suggest that P. macroloba plays a greater role in secondary stands, having a greater effect on nutrient cycle dynamics, as these soils recover from disturbance, as compared to the primary stands, which have already undergone succession, and, consequently, would have greater complexity of vegetation and soil ecosystem conditions.

2. Methods

The study sites for this project were within the MNWLR located in the Northeast region of Costa Rica, about 15 km south of the Nicaraguan border (10.7151, −84.1697). During June of 2008, we established four 20 m × 15 m (300 m2) plots both in primary and secondary forests of MNWLR in areas where P. macroloba was present as the dominant N-fixing tree for a total of eight 300 m2 plots. All habitats were once part of a single primary forest of the same terrain, soil type, and topography but have been managed differently in the past ~30 years. The primary forests are upland forests typical of the region, which have not been disturbed in the known history of the area (covering at least 100 years). The secondary forests used in this study were part of this same upland forest but were cleared and have now been regenerating into secondary forests for 15 years.

2.1. Soil Investigation

In June of 2008 and 2009, we collected 20 randomly located 2 cm wide ×15 cm deep soil cores within each plot and composited for analysis over two consecutive days and sieved at field moist conditions through 8 mm mesh. Soils were held at ambient conditions to avoid temperature shock, and analysis was performed on field moist soils at the field station within 4 days of collection. The percent of the water by weight in the soil, pH, and bulk density was determined at 10 randomly located sites within each plot. All nutrient and microbial activity data are combined results from the two years of sampling and have been adjusted for dry weight and bulk density of the soil and are presented as such.

The amount of ammonium (NH4-N), as an indicator of possible N fixation and organic N decomposition, nitrate (NO3-N), as an indicator of ammonium oxidation, and total mineral N (TMN) were determined following 2 M KCl extraction of 10 g of soil by the ammonium salicylate and cadmium reduction spectrophotometric methods, using a HACH DR 3200 spectrophotometer (Hach Company, Loveland, Colorado, 80539-0389). The TMN was calculated as the amounts of NH4-N and NO3-N. Phosphate (PO4) was measured after Bray 1 extraction (typically used in this region due to the low pH of these soils) from 2 g of soil using the molybdate reduction method and the HACH DR 3200 spectrophotometer.

Microbial biomass C (Cmic) was determined by the fumigation-extraction method as the difference between K2SO4 extracted dissolved organic carbon (DOC) levels in ethanol-free chloroform-fumigated and unfumigated 10 g soil subsamples. The DOC levels were determined by dry combustion analysis at the CATIE labs in Turrialba, Costa Rica, using the methods of Anderson and Ingram [58] and an autoanalyzer (Alliance Instruments). The rate of respiration was determined as amount of CO2 released over a 1 hour period using a Qubit SR1LP Respiration system (Kingston, ON, Canada). From these, the microbial metabolic quotients (qCO2, as a ratio of CO2 from respiration/Cmic) and the ratio of Cmic/DOC were determined to estimate the efficiency of utilization of organic C [48, 59]. The lower the qCO2 and the greater the Cmic/DOC, the more efficient the microbial community is using organic C and is presumably converting more C into biomass, rather than releasing it as CO2 [48].

The recent approach of Blagodatskaya et al. [60] and Kuzyakov [61] was used to determine ratios of primary to secondary forest DOC, respiration, and biomass to estimate the relative contribution (RC) of these individual sources of C by habitat. Also following these methods, we combined the RC ratios of the biomass to respiration and biomass to DOC to provide a comparison of the level of contribution to biomass in these soils by habitat as a further indicator of C-use efficiency. Blagodatskaya et al. [60] and Kuzyakov [61] used these methods to estimate rate of C turnover in soils.

2.2. Vegetation Investigation

In June, 2010, vegetation studies were conducted within the same plots as had previously been used for soil collection. By conducting the vegetation analyses at a different time than the soil studies, the understory vegetation was not disturbed by what typically occurs during soil collection. Within each of these 20 m × 15 m plots, we recorded the density of adult P. macroloba trees (>2 m in height) and seedlings of young trees (in classes of <1 m and >1 m). The diameter at breast height (DBH) was recorded for all trees over 2 m in height. The height of each tree was calculated using the principle of triangulation with a clinometer, and each seedling or young tree was measured using a tape measure. Basal area of P. macroloba per hectare was calculated using the mean DBH of P. macroloba in each 300 m2 plot. Volume of P. macroloba trees over 2 m in height was calculated using the formula from Tilki and Fisher [57]

Tilki and Fisher [57] applied this calculation on P. macroloba in Costa Rica, yet the volume equation uses a general form factor of 0.4 regardless of species.

Within each plot, five nested 1 m2 quadrats were randomly placed to sample understory vegetation. Within each quadrat, percent cover and density were calculated for all plant species, and the total percent cover of the understory vegetation and overstory was estimated. Cover values were determined by layer, resulting in cover exceeding 100% in individual quadrats. Vegetative biomass for 1 m2 was calculated by clipping all above-ground understory vegetation in the quadrats, bringing it back to the lab, drying in an oven at 100°C for 24 hours and weighing the dried material. Samples of plants unidentifiable in the field were pressed and brought back to be identified by local experts using Manual de Plantas des Costa Rica (Hammel et al. manuscript in preparation). Currently accepted Latin names for species according to the Angiosperm Phylogeny Group III. Classification were found using http://www.theplantlist.org/.

The species richness and abundance of vegetation other than P. macroloba were recorded with the number of stems rooted within the 1 m2 quadrats as density. The density of each species in the five 1 m2 quadrats within a plot was combined to get a representative density for the larger plot area. We conducted Shannon-Weiner () diversity analysis on both total plant species and separately for just leguminous species to compare diversity between the sites.

2.3. Data Analysis

Statistical analyses on these data were performed using the statistical software SPSS to determine if there were meaningful differences in the mean values of the parameters measured from the two different habitats. A weight of evidence approach was used involving the percent difference (PD), -test values (), and Hedge’s effect size statistic (>0.7 is considered a large effect size difference). These measures were used as indicators of biologically meaningful differences between mean values of parameters measured, as recommended for analysis of small sample sizes by Di Stefano et al. [62]. Bivariate correlation analyses were also conducted to determine the strength of the relationships between vegetation and soil variables. For vegetative data, dependent variables compared between primary and secondary forests on the 300 m2 plot scale include number of P. macroloba seedlings, number of P. macroloba adults, P. macroloba seedling height (in classes of <1 m and >1 m and average), P. macroloba tree DBH, and heights of the trees. Dependent variables pooled among the five 1 m2 quadrats in each plot included percent cover of overstory, understory, percent cover for each species, and density of each woody species and density of legumes (including P. macroloba) rooted in the quadrats, species richness and percent cover of legumes, and biomass.

3. Results

Indicators of C, N, and P cycle activity were generally greater in the primary forest than in the secondary forest. Soils of the primary forests had greater levels of P, NO3, and TMN. The DOC levels and respiration rates were greater in the secondary forest soils, but the soil C biomass levels were greater in the primary forest soils (Table 1). Consistent with these findings, the indicators of C-use efficiency showed that the soils in the primary forest were likely converting more organic C into biomass, as the qCO2 values were less and the C biomass to DOC ratio was greater in the primary forest soils. The RC of primary forest soils DOC and CO2 (as respiration) in both habitats was far less than that of the secondary forest, but the RC of primary forest soils biomass was much greater than that of the secondary forest. A comparison of the combined ratios of the RC of the biomass to respiration and biomass to DOC between the primary and secondary forest soils showed greater than a 5-fold and 3-fold contribution, respectively, in the primary forest soils (Table 2).

At the 300 m2 plot level, the number of P. macroloba young seedlings was not significantly different between the primary and secondary forests. The number of seedlings above 1 m in height was significantly higher in the secondary forest, as was the total density of adult P. macroloba in the study area. In the primary forest plots, however, adult P. macroloba was greater in height, DBH, and volume than in the secondary forest (Table 1). The primary stands were characterized by a mean P. macroloba basal area of 2.09 m2/ha, whereas the mean basal area of P. macroloba in the secondary stands was 0.148 m2/ha.

With lower seedling establishment and density of P. macroloba in the primary forest, the role of N fixation may be filled by other species. When investigating understory plant diversity in the 1 m2 quadrats, we found greater richness, biomass, and percent cover of total and leguminous plant species as well as biomass in the primary forest. There were nine legume species in the primary forests, and all three of the species that occurred in the secondary forests (including P. macroloba) were also found in the primary forests (Table 3). There was a significant difference not only in the number of leguminous species in the primary versus secondary forest stands, but also in the percent cover of the quadrats that these species occupied (Table 1). It was not confirmed through this study if these legumes had N-fixing nodules, and since some were not identified to the species level, it is not possible to find existing literature regarding the identity of N-fixing symbionts for each species.

The Shannon-Weiner index confirmed the higher diversity of legume species in the primary forest () when compared to the secondary forest (). Total plant species richness was also higher in the primary forest stands (total of 52 species over the surveyed 1,200 m2 of primary forest) than the secondary stands (total of 43 species over the surveyed 1,200 m2 of secondary forest). For all species in the sampled area of the primary forest, and in the secondary forest .

Correlation analysis showed that a greater volume of P. macroloba and percent cover of leguminous vegetation in the primary forests was correlated with a general enhancement of the indicators of microbial activity in the C, N, and P cycles. The number of young P. macroloba > 1 m and the overall height of the seedlings decreased with an increase in P. macroloba DBH (; ), height (; ), and volume (; ). Percent of the leguminous cover increased with an increase in P. macroloba DBH (; ), height (; ), volume (; ), a decrease in the total number of P. macroloba trees (; ), and to a lesser degree, an increase in young tree height (; ) and a decrease in the number of P. macroloba seedlings >1 m (; ). Also, the number of legumes present on the forest floor followed the same pattern of increasing with a decrease in the number of P. macroloba (; ) and an increase in the P. macroloba DBH (; ) and volume (; ). Increases in P. macroloba DBH, height, and volume occurred with increases in soil nitrate (; ), total mineral N (; ), the percent of N as nitrate (; ), the DOC (; ), the C biomass (; ), indicators of a more efficient use of soil C (qCO2: ; ; C biomass per unit of DOC: 0.693–0.834; ), and with decreases in the respiration rate (; ).

4. Discussion

The importance of secondary forests is evident in Costa Rica, where nearly all primary lowland forests outside of protected areas are disturbed or converted by deforestation, and the amount of abandoned agricultural land continues to increase [1, 9]. Trends found in this research can be used as indicators of success in tropical wet lowland forest management projects. In areas where agricultural land use has been abandoned and soils left degraded, Chazdon [12] points out that rehabilitation through planting of carefully selected trees can improve soil fertility. Additionally, Yang et al. [63] and Jin et al. [64] showed that as degraded tropical forests undergo succession naturally, the quality and quantity of the forest litter, root zone, and below-ground indicators of soil ecosystem conditions were more enhanced when compared to adjacent tree plantation soils. The current research emphasizes the role of P. macroloba as a pioneer tree species in the region, as we observed that P. macroloba is a key driver of succession in both above- and below-ground communities. Similar to the findings of Yang et al. [63] and Jin et al. [64], we found that the increase in total and leguminous vegetation complexity was positively correlated with both the mass and volume of P. macroloba during forest succession, and the metrics were used to indicate enhanced soil ecosystem activities.

The presence of P. macroloba at high densities in secondary forests hastens the recovery of soil nutrients such as N and organic C and perhaps facilitates the development of a more complex N-fixing vegetation community following disturbance. While the P. macroloba seedlings established at similar rates in the secondary and primary forests, we saw a change in leguminous species assemblages between successional stages, with P. macroloba as one of the only leguminous species in the secondary stands. It complements the work of Peña-Claros [65] who found that canopies of tropical secondary forests tend to be dominated by a small number of pioneer species, and as these pioneer tree species begin to die off, an increase in species diversity is expected because other tree species reach the canopy. Similarly, our primary forests were higher in plant diversity and had P. macroloba as one of several, albeit the dominant, N-fixing trees that occurred with an increased leguminous herbaceous population to enhance N and C cycle dynamics. For example, three species of the Inga genus were found in the sampled areas of the primary forest, and only one in the secondary forest sampling areas. This genus of tree may also be playing an important role in these primary forests, as they are believed to be critical ecosystem components within their range, serving as important sources of inorganic N [66] and enhancing the quantity and quality of the soil organic C [6770], N, and P [69, 70]. In addition, Inga species have been found to stimulate greater above-ground biomass development than other non-Inga species in their habitats [68, 71] and are thought to stimulate C sequestration and long-term N accumulation [67, 68]. It appears that both Inga spp. and P. macroloba along with the increased understory leguminous vegetation in the primary forest are playing the important role of enhancing soil ecosystem nutrient quantity and quality one would expect in these forests [30, 32, 33], while P. macroloba dominates this role in the secondary forests.

The primary forest soils had greater amounts of nitrate and a greater percent of the inorganic N as nitrate but similar levels of ammonium as the secondary forest soils, suggesting either greater or more efficient ammonium oxidizing activity in the primary forests. The enhanced microbial decomposition of organic matter associated with increases in successional stage of forests can also provide amine and other reduced N groups used for biomass development and also be used as an electron source for ammonium oxidizers [51]. These suggestions should be confirmed in future studies; however, they would be expected to occur in concert with the increasing vegetation complexity that occurs during secondary forest succession [22, 2731].

The primary forest soils’ greater amount of soil C biomass, contribution of CO2 and DOC into soil C biomass, rate of C turnover, and efficiency of organic C use further illustrate the complexity of the soil community in these older successional stage forests. These findings can be indicators of a more fungal dominant population decomposing the organic matter more efficiently for incorporation into the greater vegetation biomass [48, 59, 72, 73].

Greater levels of inorganic P, such as those found in the primary forest soils, have been associated with increased decomposition of the more complex and greater amounts of woody debris and litter [26, 36, 37, 41, 42, 74]. It is of interest to note that the role of inorganic P in soil microbial activity in the tropical soils of this region remains unclear, as differences in observations have been reported from soils in similar forests of the Northern Zone of Costa Rica since 2002. Some researchers have shown that increases in inorganic P have been associated with increases in decomposition activity, soil ammonium oxidizing activity, rates of incorporation of inorganic N into cells, increased efficiency of C use and soil biomass development, and stimulation of the fungal community activity [19, 3840, 73, 75, 76]. Others have found the opposite relationship [44, 77, 78]. A pattern similar to the first scenario was observed in the primary forests in the current study, as the increased soil biomass development and increased C-use efficiency were associated with much greater levels of inorganic C and a more complex vegetation community in the primary forests. More work is clearly needed to understand the role of P in microbial biomass development in these soils.

The correlation analyses conducted on the data from this study showed that as the forests get older, there are fewer P. macroloba trees, which, not surprisingly, are larger and make up a greater amount of the volume in the forests. Along with this increase in size of these trees, there is a corresponding increase in the number and percent cover of legume species in the understory. These differences in vegetation structure corresponded with increased values of the indicators of microbial activity associated with the C and N cycles, and an enhanced efficiency of the use of the organic C by the soil biota [48, 52, 79]. Thus, it appears that while the greater density of younger P. macroloba is important in the C, N, and P recuperation of the earlier stages of secondary forest regeneration, perhaps this role shifts to the more diverse leguminous understory vegetation in the older forests. It would be valuable to study the different ages of secondary forests in the region to confirm that this shift does occur and determine at what age of secondary forest it takes place.

The patterns observed in the current study are consistent with what has been found elsewhere when comparing younger forests recovering from disturbance with more established primary or old growth forests. The increased amount of woody debris, litter, and leguminous vegetation that develops in the older forests is associated with an enhanced microbial community, which, typically, produces an accumulation of inorganic N and P nutrients and organic matter. These processes ultimately result in an increased efficiency in the use of the organic matter and biomass development [36, 37, 48, 59, 7881].

5. Conclusions

Despite the limited dataset in this preliminary study, we did find that these primary and secondary lowland tropical wet forest stands differ in both their above-ground diversity and below-ground activities. We found support of our hypothesis that P. macroloba is critical for soil nutrient recovery in secondary forest stands as the primary N-fixing species. Future research conducted by the authors will investigate competitive relationships between plant species in this region and study whether the replacement of P. macroloba by other leguminous species in the older growth forests is by virtue of its inability to outcompete these species in areas in later successional stages. We have begun to expand this research to include additional series from earlier and midsuccession.

More studies such as this one are needed to understand the connections between below-ground activities and vegetative biomass and diversity in tropical forests. Such research demonstrates how pioneer species and other key groups affect the development of secondary forests into older stage forests and can be used to develop management strategies for recuperating soil nutrients in such degraded tropical lands. It appears that P. macroloba is an important pioneer species in forest recovery that is playing this role but is also a part of a more complex system of N-fixing vegetation that develops in the primary forests that is associated with increased soil microbial activity, production of inorganic P, and driving the conversion of both the N and C into biomass development that is ultimately resulting in a more efficient use of the organic material generated in these older forest soils.

Increasingly, researchers and land managers are recognizing plantations and managed forests as important tools for restoration of abandoned tropical pasturelands that can support diverse plant communities in addition to accomplishing economic objectives [8183]. Plantations or managed forests with P. macroloba can be used to catalyze successional processes and have been found to sustain high plant diversity and organic carbon [83]. Our findings illustrate the fundamental role that P. macroloba plays in nutrient cycling in early successional forests. Restoration that supports either existing P. macroloba recruitment or adds this species to a site from which it has been extirpated can provide a method of minimum intervention that removes successional impediments. Our research suggests thata high density of P. macroloba in degraded lowland tropical forests may initiate a self-repair process that will lead toward a properly functioning ecosystem in which nutrient cycling is efficiently performed by a greater diversity of leguminous vegetation in the later successional stages.

Conflict of Interests

There is no conflict of interests between the authors of this paper and companies used for analysis (Hach Company, Loveland, Colorado, and IBM SPSS statistical software, Armonk, NY, USA).

Acknowledgments

This work was supported by the National Science Foundation Research Experiences for Undergraduates Program (Grant no. 1034896) and the Kean University Office of Research and Sponsored Programs. The authors are very grateful to the staff of Laguna Del Lagarto Lodge for hosting them when they conducted this research and to Orlando Vargas (Organization of Tropical Studies, La Selva Station) and Nelson Zamora (Instituto Nacional de Biodiversidad) for assisting with plant identification. They also thank Kate Niemiera, Katie McGee, Ashley Casimir, Caitlin Looby, Brandon Nwobu, and Dana D’Angelo for their assistance in the field and laboratory.