Next Article in Journal
Using Macronutrient Distributions within Trees to Define a Branch Diameter Threshold for Biomass Harvest in Sugar Maple-Dominated Stands
Next Article in Special Issue
Effect of Soil Moisture on the Response of Soil Respiration to Open-Field Experimental Warming and Precipitation Manipulation
Previous Article in Journal
Assessing Precision in Conventional Field Measurements of Individual Tree Attributes
Previous Article in Special Issue
Spatial Upscaling of Soil Respiration under a Complex Canopy Structure in an Old‐Growth Deciduous Forest, Central Japan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temporal Variability of Soil Respiration in Experimental Tree Plantations in Lowland Costa Rica

Department of Ecology, Evolution & Organismal Biology, Iowa State University, Ames, IA 50011, USA
Forests 2017, 8(2), 40; https://doi.org/10.3390/f8020040
Submission received: 10 January 2017 / Revised: 1 February 2017 / Accepted: 2 February 2017 / Published: 8 February 2017
(This article belongs to the Special Issue Forest Soil Respiration under Climate Changing)

Abstract

:
The principal objective of this study was to determine if there is consistent temporal variability in soil respiration from different forest plantations in a lowland tropical rainforest environment. Soil respiration was measured regularly over 2004 to 2010 in replicated plantations of 15- to 20-year-old evergreen tropical trees in lowland Costa Rica. Statistically significant but small differences in soil respiration were observed among hours of the day; daytime measurements were suitable for determining mean fluxes in this study. Fluxes varied more substantially among months, with the highest average emissions (5.9 μmol·m−2·s−1) occurring in September and low emissions (3.7 μmol·m−2·s−1) occurring in January. Three of the six tree species had significantly increasing rates of soil respiration across 2004–2010, with fluxes increasing at an average of 0.09 μmol·m−2·s−1 per year: the three other species had no long-term trends. It was hypothesized that there would be a tradeoff between carbon allocation aboveground, to produce new leaves, and belowground, to sustain roots and mycorrhizae, but the relationship between canopy leaf fall—a surrogate for canopy leaf flushing—and soil respiration was significantly positive. The similarities observed among temporal trends across plantation types, and significant relationships between soil respiration, soil water content and soil temperature, suggest that the physical environment largely controlled the temporal variability of soil respiration, but differences in flux magnitude among tree species were substantial and consistent across years.

1. Introduction

The tropical forest biome covers only about 12% of Earth’s land surface [1,2], yet it harbors much of the world’s biological diversity [3,4,5,6] and disproportionally influences global biogeochemistry. Tropical forests are responsible for one-third of the world’s terrestrial plant productivity [1,2,7], one-third of the world’s terrestrial respiration [8] and one-third of terrestrial evapotranspiration [9]. Although they vary greatly in composition and structure, tropical rainforests often contain substantial living biomass, and large stocks of detrital carbon and organic nitrogen [10,11,12,13] that sustain high rates of microbial activity. In brief, the warm, moist conditions that prevail year-round in tropical forests support abundant life and rapid land-atmosphere exchanges of energy, water and carbon. Both warmer temperatures and an accelerated hydrologic cycle are expected in tropical forest regions in future years [14]. Those changes will alter tropical forests and soils and, thus, land-atmosphere exchanges of carbon, energy and water. We do not fully envision what the resulting impacts will be.
Large stocks of detritus, warm temperatures, and humid conditions together create a strong potential for warming (without drying) to stimulate decomposer activity and the mineralization of detritus to CO2 and other trace gases. That could generate a positive feedback to global warming, because CO2 is an important greenhouse gas [15,16,17]. The likelihood that this particular feedback will emerge has been extensively discussed [18,19,20,21,22,23] but has not been ruled out [24,25,26,27,28,29]. As a result, despite decades of discussion and debate, the question persists. Although greater warming is expected in northern latitudes than in the tropics [14], CO2 fluxes from humid tropical forests are already large, and small changes in large fluxes can be important. For example, a 3% increase in a tropical soil respiration of 1500 g·C·m−2·year−1 is greater than a 20% increase in a tundra soil-respiration rate of 200 g·C·m−2·year−1. Based on studies of moist tropical forests, it has been proposed that decomposition rates increase faster at warmer temperatures than does net primary productivity (NPP) [30]. That would be consistent with a positive-feedback scenario for tropical forest soil respiration. A variety of evidence suggests that the tropics already influence interannual variations in atmospheric CO2 concentrations [31]. Measurements of soil respiration are insufficient, by themselves, to address this debate, but tropical forests play an important role in the earth’s carbon cycle; they will be impacted by climatic changes; they are likely to influence future atmospheric changes; and they remain understudied in relation to their potential to cause change.

Objectives

The over-riding objective of this study was to determine if there is consistent temporal variability in soil respiration (Rsoil) from evergreen forest plantations in a lowland tropical rainforest environment. To do so, Rsoil was measured over diel, monthly and multi-annual time frames. Second, I sought to identify the environmental variables underlying observed temporal variations in Rsoil, to inform the process-based understanding of tropical forest carbon dynamics that is needed to improve predictive modeling [32]. Despite an increasing body of empirical data from tropical forests, there is a lack of agreement on several issues. For instance, diel variations in Rsoil have been observed in some evergreen tropical forests, but not in others [33,34,35,36,37,38]. That issue underlies an important concern: most measurements of tropical forest Rsoil are made only in the daytime, especially in remote locations. If fluxes are significantly different at night than during daylight hours, then we have a poor grasp of what average tropical forest Rsoil truly is [35]. Additionally, daytime photosynthesis by a productive tree canopy may drive diel variations in Rsoil, via temporally variable canopy-to-root carbon translocation [39,40,41]. Over longer time frames, most studies indicate that temporal variability in tropical forest Rsoil is attributable to variations in soil moisture content or temperature [37,38,42], indicating that climatic variability is important. However, plant phenology may independently influence the temporal dynamics of Rsoil, via seasonally variable rates of aboveground litter production or root growth or turnover [43]. Furthermore, if Rsoil is limited by carbon (C) fluxes to roots, then trade-offs between aboveground and belowground C allocation may also cause seasonal variability in Rsoil. For instance, belowground fluxes may be lower at times when forest canopies are producing a new crop of leaves. If so, that might suggest that root-shoot allocation patterns vary seasonally within a forest, and that carbon limitation influences the magnitude of Rsoil, at least during some months.
I tested the following hypotheses, based on the measurements collected. First, Hypothesis 1: soil respiration varies on a diel basis. Such temporal variations in Rsoil have important implications for interpretation of daytime-only measurements, and might result from a variety of processes. For instance, lowland, humid tropical forests have a unique climatic feature: diel variations in temperature often exceed the variations in mean monthly temperatures. Warmer daytime temperatures may drive faster daytime CO2 production in, and diffusion from, soils. Also, observed night-time warming is greater than observed daytime warming in the tropics, probably due to increased cloud cover [44,45,46,47], yet most data available on soil CO2 emissions from tropical rainforests are derived from daytime measurements. In tropical forests, soil temperatures typically exceed air temperatures at night, and that could promote movement of CO2-rich soil air into the atmosphere. I also tested Hypothesis 2, that Rsoil varies among months of the year. Even within relatively aseasonal wet tropical forests, there is climatic variability and there are phenological differences among tree species [48] that likely influence the magnitudes of aboveground and belowground C fluxes. However, seasonal patterns in belowground carbon fluxes in evergreen tropical forests are very poorly understood. To further evaluate seasonal variability, I also tested Hypothesis 3, that Rsoil and canopy leaf flushing are negatively correlated with one another. This hypothesis is based on the supposition that there is a trade-off between plant carbon use aboveground, to produce a new crop of canopy leaves, and belowground, to support root systems. Leaf flushing by tropical evergreen broad-leaved trees represents an important and seasonally variable investment of available carbon in a new crop of leaves. As new leaves expand out of their buds, old leaves are dropped. Thus, leaf flushing appears as seasonality of leaf fall by canopy trees. Finally, I tested Hypothesis 4, that Rsoil increases with forest age in developing forest stands, as suggested by [49]. This hypothesis is put forth to determine if there is a multi-annual trajectory in Rsoil, either positive or negative, in the 15- to 20-year old experimental tree plantations that were the focus of this research. It is widely accepted that forest productivity increases with age in young forests, to a point, but does not continue to increase continuously, and very often decreases as forests mature [50,51,52]. Such non-linear trends may be difficult to disprove, but tropical tree plantations can grow rapidly [53]; quantifying belowground dynamics in a variety of different forest types as they mature and age, within a single rainforest location, can provide valuable insights into the temporal dynamics of belowground processes, about which we know very little.

2. Materials and Methods

2.1. Study Site

This research was undertaken at the La Selva Biological Station of the Organization for Tropical Studies, in the Caribbean lowlands of northeastern Costa Rica (10°26′ N, 84°03′ W). Over 1997–2009, annual precipitation averaged 4537 mm and mean air temperature was 25.1 °C [54]. The driest season at La Selva generally extends from February through April, but rainfall averages > 150 mm·mol−1 in every month. Forest evapotranspiration averages > 2000 mm·year−1 [55]. The native vegetation at this site is species-rich broad-leaved evergreen rainforest that has a high abundance of Pentaclethra macroloba (Willd.) Kuntze (Fabaceae, Mimosoideae) trees, and many subcanopy palms. In primary forest on rolling terrain, there are an average of 528 trees·ha−1 having a basal area of 23 m2·ha−1 [56]. Long-term studies of old-growth forest at La Selva show that litterfall averages 9.1 Mg·ha−1·year−1, 78% of which is leaves [54]. Tree growth averages half that and their sum, aboveground net primary productivity (ANPP), averages 13.7 Mg·ha−1·year−1 [54]. More information about the physical and biological features of La Selva is provided in [57].
Field studies for this research were conducted in replicated experimental tree plantations that were established on formerly grazed pastures, to test the influence of tree species identity on soil properties after reforestation [58,59]. The pastures were established from mature forest in the mid-1950s, and trees were planted, after removal of the cattle, in the winter of 1988–1989 [53,58,60]. This study began 15 years later. The study plots were at 10°26′ N, 83°59′ W, on hilly terrain with elevations of 44 to 89 m. They were approximately 3.3 km from the biological station and weather station. The soils were Oxisols derived from volcanic parent materials, and belong to the Matabuey consociation [61]. In brief, they are deep, acidic, highly permeable clays with low base saturation and relatively high soil organic matter contents [61,62,63]. Surface-soil pH (in water) was <4.5 in all plots [63]. The available data indicate that total soil C and nitrogen (N) stocks were lower in the plantations than in nearby, undisturbed forest, on average; and that soil organic C mineralization rates (gCO2-C gSoil−1·day−1) were lower in plantation soils than in mature forest soils [8,53,63]. Nevertheless, rates of N cycling in the plantations were very high [64]. Fine-root biomass was concentrated in the surface 15 cm of mineral soil, and declined rapidly with increasing soil depth to 1 m [65]. The forest floors were comprised primarily of leaves and particulate debris that decomposed rapidly: forest floor turnover rates among the plantations we studied ranged from 1.5 to 2.3 year−1 [66]. The diversity and abundance of litter-dwelling organisms were high.
The experimental design included four randomized blocks composed of twelve 50 m × 50 m plots, each of which originally contained a single tree species planted at 3 m × 3 m spacing, except that the control plots were not planted or otherwise managed in any way [58,67]. Each block was centered across a single hill and thus encompassed a range of slope and aspect positions. Plantations of the fastest growing species were thinned three years after planting, and additional thinning was done three years later [67,68]. Each of four replicate plots of six plantation types, a total 24 plots, were included in this study. These included plantations of Pinus patula subsp. tecunumanii (Eguiluz and J. P. Perry) Styles; Hieronyma alchorneoides Allemao; Pentaclethra macroloba (Willd.) Kuntze; Virola koschnyi Warb.; Vochysia ferruginea Mart.; and Vochysia guatemalensis Donn. Sm. Pentachlethra is a Mimosoid legume with occasionally nodulated roots; Pinus patula was the only conifer and the only non-native species included in this study. Site management included understory clearing over the first four years, and thinning, generally to 50% of the standing basal area (i.e., every other tree), but the timing of thinning varied among species [67]. Hieronyma and Vochysia spp. coppiced after being thinned, and Pentaclethra apparently was never thinned. In 2004, when this study began, these plantations had an average of 405 trees·ha−1 with a total basal area of 23 m2·ha−1 and average tree heights that ranged from 14 to 26 m (Table 1). Over the course of this study, there was a decline in tree stem densities whereas overstory basal areas and tree heights increased, as is typical in growing forests. The exceptions were plantations of Pinus patula and Vochysia ferruginea, which lost basal area as a result, apparently, of fungal attacks. The understories of the experimental plantations developed dense and diverse understories [67,68] that continued to grow throughout this study.

2.2. Measurements of Soil Respiration

Measurements of soil respiration were initiated in September 2004 in each of the four replicate plots of six plantation types and continued until February 2010 in four of the plantation types. All studies were based on LI-COR© automated soil gas flux instruments (LI-COR Biosciences, NE, USA) that were returned to the factory for cleaning and recalibration annually. The specific measurement chambers utilized varied, but they provided fully comparable measurements of in situ soil respiration, Rsoil. I herein report the results of two principal studies: diel study and survey measurements. The diel study characterized variations in Rsoil each hour across 2-day and longer periods, based on automated sampling with a single instrument that was moved to sample each of the 24 study plots at least twice. The Survey study involved sampling each of the experimental plots at one time, during daytime hours, so that the effects of the experimental treatments could be cleanly compared. It sometimes took two days to sample all the chambers in all plots. Measurements were conducted through time so that seasonal variations in Rsoil could be quantified. Summary results of the first two years of this study were reported earlier [53,65], but data from 2008 to 2010 have not been previously reported, nor have data from plantations of Vochysia ferruginea or from the diel study. Based on the data collected during both studies, I also tested for longer-term trends in fluxes in the maturing plantations (Hypothesis 4), which were 16 years old in 2004 and 21 years old in 2009. Thus, hourly, monthly, and inter-annual variations in Rsoil were assessed for this article.

2.2.1. Diel Measurements

Soil CO2 emissions were measured every hour over a total of 52 continuous periods of >2 days between November 2004 and February 2010, with a single chamber in one plot at a time, to investigate diel variations in Rsoil. Measurements were made with a LI-COR 8100 soil CO2 flux system attached to an 8100-101 long-term chamber (LI-COR© Environmental, Lincoln, NE, USA) and powered by a deep-charge marine battery, all of which were back-packed into an experimental plot and then left in place for the duration of measurements. The battery allowed for measurements to be collected every hour over 2–3 day periods. A day before measurements began, a single 20-cm diameter, 12-cm tall polyvinyl chloride collar was inserted approximately 2 cm deep into the surface soil of a plot by carefully cutting through the forest floor with a sharp knife. All measurements were at randomly located positions without regard to tree positions. On the measurement date, the automatic chamber was carefully situated such that it closed firmly on top of the collar, and lifted freely after each measurement. During measurements, the chamber remained tightly closed for 2 min, and within-chamber CO2 concentrations were monitored as the chamber headspace was circulated through the IRGA (Infrared Gas Analyzer, total system volume averaged 7 L). The first 30 s of measurements were ignored to allow for internal mixing to complete. Over the remaining 90 s, the concentration of CO2 within the chamber was measured every second, and the rate of change in headspace CO2 concentrations was used to quantify Rsoil. The resulting time-series data were analyzed using LI-COR’s embedded FV8100 file viewer software. In most cases, exponential fits were utilized because the rate at which within-chamber CO2 concentrations increased typically declined through time. The R2 values of curve fits averaged 0.99 (N = 2650) and the minimum R2 was 0.75. The first hourly measurement collected at a new location was discarded because of consistently high ambient CO2 concentrations that were attributable to human respiration. A summary of the 52 measurement periods is provided in Table S1. Due to time commitments to other studies, no diel measurements were made from May through September.
During diel measurements, near-surface air temperatures, relative humidity, atmospheric pressure, and CO2 concentrations were monitored continuously as the measurement chamber closed and then within the closed chamber. Mineral-soil temperature at 10 cm depth was monitored with a soil thermistor probe (LI-COR 8150-203). Surface-soil volumetric moisture content was assessed at four locations surrounding the chamber, at the beginning and end of each measurement period, using a CS620 Hydrosense system with 12-cmrods (Campbell Scientific, Inc., Logan, UT, USA). Following successful completion of one measurement period (i.e., one continuous sequence of hourly measurements at one plot, typically lasting >50 h), the equipment was disassembled and relocated to a different plot. All 24 plots were measured twice between January 2008 and February 2010, to capture both dry and wet periods (Supplementary Table S1). Additional measurements collected in 2004, before sensor malfunctions halted measurements, are included in the analyses (Table S1).
Weather data encompassing 2008 through 2010 were obtained from the La Selva meteorological station, which was approximately 3.3 km from the experimental plots. These data included half-hourly air temperature (°C), relative humidity (RH, %), rainfall (mm), mean wind speed (m·s−1), solar radiation (W·m−2) and photosynthetic photon flux density (PPFD, μmol·m−2·s−1). To meld data into a single file for analyses of flux-×-environment relationships, measurements were matched to the nearest half-hourly meteorological data. I also assigned each measurement to its nearest hour so that ‘hour-of-day’ could be used as a categorical variable for statistical analyses. Detailed weather data were not available for 2004.

2.2.2. Across-Site Surveys of Soil Respiration

The study site was on hilly terrain, and measurements were conducted with backpacked instruments and were limited to daytime hours. Altogether, 99.5% of the measurements were collected between 07:30 and 16:30 sunrise and sunset were at approximately 06:00 and 18:00 daily. The survey measurements were undertaken to test for significant influences of plantation type (tree species) on the magnitude of Rsoil. These data are utilized herein for three purposes: (a) to quantify monthly variability in Rsoil (Hypothesis 2); (b) to test Hypothesis 3, which posits a temporal offset between aboveground and belowground C fluxes; and (c) to identify longer-term (2004–2010) trends in Rsoil (Hypothesis 4). The survey measurements undertaken for this study are summarized in Table S2.
Soil respiration within the plantations was measured 66 times between August 2004 and March 2010 using a LI-COR® 8100 soil CO2 flux system. From 2004 through 2008, measurements were made with an 8100-102 chamber placed on top of a 10-cm diameter, 5-cm tall plastic collar that was carefully inserted by cutting with a sharp knife through the forest floor and about 1 cm deep into the mineral soil. There were 3–4 collars per plot (12–16 collars per species each date), and they were moved to new, randomly selected locations every year or whenever they were disturbed. At each collar on each measurement date, soil CO2 fluxes were monitored every second for 90 s. Fluxes were calculated over the final 70 s of chamber closure, using the embedded LI-COR file viewer software, with the first 20 s being excluded to allow for full mixing of the within-chamber atmosphere. The mean flux from the final 70 s was based on an exponential fit of the CO2 concentration over time relationship. Initial measurements included six species, in one plot per block of each. Measurements in Vochysia ferruginea were discontinued after May 2006 because entire trees began falling as a result of butt rot, which seemed to spread among adjacent trees via root grafts. Measurements were discontinued in Pinus palustris after March 2008 because trees dropped all their needles and disintegrated. From February 2009 through February 2010 an 8100-103 (20-cm diameter) chamber was utilized, with four 20-cm diameter collars in each of the four plots of each of the remaining four tree species (Supplementary Table S2). Within-chamber CO2 concentrations were monitored every second for four minutes, with the first 40 s being allowed for full mixing of the air within the chamber, and the subsequent 200 s being used to calculate CO2 emissions. It was sometimes possible to collect data from all soil collars within all four blocks in one day, but it often took two. Each measurement cycle began at a randomly selected block to minimize the potential for temporal bias to influence comparisons among blocks and species.
At the time of soil respiration measurements, a thermistor probe (LI-COR 8100-201) was inserted 5 cm into the mineral soil approximately 20 cm distant from the measurement collar at three locations, to measure soil temperature. In 2009, we started measuring soil temperature (Tsoil) at 10 cm depth. After chamber closure, we measured surface-soil water content with a hand-held electronic soil moisture sensor that was inserted into three locations around but >20 cm from the chamber. Early measurements were made with an ECH2O Dielectric Aquameter (Model EC-10, Decagon Devices, Inc., Pullman, WA, USA; LI-COR part 8100-202). In 2005, that instrument was replaced with a Campbell 620 Hydrosense system with 12-cm long probes (Campbell Scientific, Inc., Logan, UT, USA), and it proved more durable under the site conditions. We cross-calibrated the two water-content sensors such that they provided the similar readings across the breadth of observed soil moisture contents. Local weather data were obtained from the La Selva Biological Station, which was 3.3 km from the closest study plots.

2.3. Canopy Leaf Fall

Total fine litterfall was measured with four traps per plot that were emptied at the middle and end of each month. Traps had wooden frames that had 1.3 × 0.4 m internal dimensions and 2-mm mesh screen bottoms. They were supported about 30 cm above the soil on steel legs. The collected materials, which included branches ≤ 1 cm in diameter and all non-woody materials, were combined within months to generate a single sample per plot that was sorted into four fractions: branches and bark; canopy leaves; other leaves; and miscellaneous materials [65,66]. Traps were removed, repaired and repositioned annually. More than four years of litterfall measurements were made from October 2003 to December 2009. They are used herein to test for temporal offsets between leaf fall and soil respiration (i.e., Hypothesis 3).

2.4. Statistical Analyses

Statistical analyses were conducted using JMP® 11.2.0 (SAS Institute, Inc., Cary, NC, USA). The underlying experimental design was a randomized complete-block with four blocks and six plantation types, which are referred to by the single tree species that was planted into each plot. A combination of parametric and non-parametric statistical analyses was applied to test for significant differences among the test variable, typically time (hour or month or year); tree species; and their interaction. Tests of hourly and monthly data were based on binned values, 24 h of day and 12 months of year, rather than the continuous time variables that were collected. Several features characterize the datasets. First, the hypotheses were designed to test for landscape-level patterns but, commonly, results varied among the plantation types. Also, sample sizes differed among treatments. Therefore, hypotheses typically were tested at both the all-plot (landscape) and treatment (tree species) levels. Homogeneity of variances was specifically tested with the Brown-Forsythe and Levene tests and was very frequently rejected. In those cases, Welch’s test was applied. In cases where data were not normally distributed, transformations were applied prior to statistical analysis. For instance, measurements of daytime Rsoil (N = 4962) were not normally distributed and included occasional very high fluxes (maximum = 46 μmol·m−2·s−1). Natural-log transformation of the data reduced kurtosis from 54.7 to 0.8, reduced skewness from 3.6 to −0.13, and changed the median-to-mean ratio from 0.94 to 1.01. Statistical tests of Rsoil therefore were based on ln(Rsoil), but results typically are presented as actual values to facilitate understanding. In all cases of comparisons among groups, a non-parametric Kruskal–Wallis test was applied to validate the results of parametric analyses of variance; reported results represent the more conservative result (i.e., lowest P value). Differences among factors were considered significant at P = 0.05.
The diel study comprised 52 independent observational studies that were designed to determine the effect of time-of-day on Rsoil. Each study included automated measurements of Rsoil that were made every hour over more than two continuous days, at one location in one plantation type. Measurements were made a minimum of two times in each plantation type in each block, at different times of year (Table S1). Thus, the measurements incorporated variability attributable to tree species, hour-of-day, and season-of-year, all of which were found to significantly affect fluxes, at least sometimes. Therefore, to quantify the hour-of-day effect, each of the 39 to 70 hourly measurements made during each of the 52 measurement periods was normalized to a measurement-period mean value of precisely 1.0. Within an individual measurement period, then, fluxes varied through time but had an overall average normalized flux rate of one. As a result, all measurement periods had fully comparable fluxes that varied only with time of day. The normalized data then were combined for statistical analyses. The resulting dataset had a sample size of 2650, a median = 0.99, a range of 0.025–1.798, a CV (coefficient of variation) of 0.14, a skewness of −0.12, and a kurtosis of 6.5. To test for diel variations in Rsoil (Hypothesis 1), data were assigned to hourly bins (categorical variable, hour 14 = 13:30–14:29).
The survey study similarly comprised individual observational comparative studies during which Rsoil was measured within plantations within all four blocks. This study included 66 such measurement periods: the species being measured varied through time (Table S2). To test for seasonal differences (Hypothesis 2), measurements were assigned to monthly bins. Block was treated as a random effect. Tree species and their interaction significantly influenced Rsoil (least squares ANOVA, N = 4858), so comparisons among months were done on a per-species basis. In most cases, variances among months were not homogeneous, so the results of Welch’s test are reported. In cases where variables were not normally distributed, the non-parametric ranked-sum Kruskal–Wallis Test was applied.
To compare the seasonality of leaf flushing with that of Rsoil, to test Hypothesis 3, mean monthly soil respiration data were paired with mean monthly canopy leaf fall data from the same species and time frame. Pearson product–moment correlations were calculated on a single dataset including all months and species, i.e., mean monthly canopy leaf fall and mean monthly Rsoil for each month and species (N = 72).
To assess whether there were long-term trends in Rsoil over the duration of the study (Hypothesis 4), data from both the diel and survey studies were combined into a single dataset, and the hypothesis was tested independently for each species. Least-squares linear regression modeling was applied with ln(Rsoil) as the dependent variable. Results from simple correlations, non-linear models, or multiple linear models having additional time-related variables (e.g., time2) added no meaningful information, and are not reported.

3. Results

This study ran from October 2003 through February 2010, during which largely normal weather patterns prevailed at La Selva. Rainfall averaged 4420 mm·year−1 and temperatures averaged 25.1 °C.; both are similar to longer-term means (Table S3). Dry-season precipitation (January through April) was below average in 2007 and 2008, but exceeded 65 mm in every month. Overall, both the diel and survey measurements included a broad range of the weather conditions that typify the climate of La Selva, including brief rain-free periods.

3.1. Diel Variability in Rsoil and the Environment

Soil respiration varied widely within and among days, without obvious regularity (Figure 1). Nevertheless, normalized fluxes did vary significantly among hours of the day (Kruskal–Wallis test, df = 23, χ2 < 0.0001). This finding supports Hypothesis 1, that Rsoil varies on a diel basis. However, the R2 was low (0.02) and variability within hours was high. Overall, Rsoil was lowest between 06:30 and 10:30, when it averaged 4.03 μmol·m−2·s−1 (Figure 2a). There was no equivalent identifiable period of maximum emissions; all other hours of the day had statistically equivalent fluxes.
Several meteorological variables that potentially influenced emissions also varied through time (Figure 2) and, thus, potentially provide insight into environmental controls over within-soil CO2 production and its rate of escape into the atmosphere. Note that data in Figure 2 represent averages calculated from hours during which Rsoil was measured. Across all measuremement dates, soil temperature (Tsoil) at 10-cm depth averaged <24.2 °C from 07:30–12:30 but >24.5 between 16:30 and 23:30 (Figure 2c). Air temperatures increased from an average daybreak minimum of 21.6 °C to a post-noon maximum of 27.3 °C (Figure 2d). The temperature difference between the soil and atmosphere varied more widely (Figure 2b), ranging from −3.1 °C at 13:00 to 2.7 °C at 06:00. Atmospheric relative humidities were quite high, on average (Figure 2e), but they too varied across the day: they averaged >93% between 02:30 and 08:30 but <83% between 11:30 and 16:30. Among those variables, Rsoil most closely paralleled changes in soil temperature, and increased with Tsoil at 10 cm depth with a Q10 of 3.3 (Least-squares linear regression, N = 2432, R2 = 0.07, P < 0.0001, temperature exponent = 0.1204 ± 0.0086 per °C).
Ambient CO2 concentrations above the forest floor varied significantly with hour-of-day in plantations of each of the study species (GLM (General Linear Modeling procedure), N = 398 to 539, χ2 < 0.0001) and across all species combined (N = 2650, χ2 < 0.0001). Carbon dioxide, presumably from the soil, accumulated on-site when wind speeds (measured in the open) averaged <1 m·s−1, and atmospheric CO2 typically exceeded 500 μmol·mol−1 during the nighttime (Figure 2g). Ambient CO2 and wind speed (Figure 2h) were strongly negatively correlated (N = 2492, r = −0.61, P < 0.0001), based on wind-speed measurements at the La Selva Biological Station, about 3.5 km from the field plots. Also measured at that weather station, over 2008 through 2010, were solar irradiance and photosynthetic photon flux density (PPFD). During diel measurements of Rsoil, irradiance averaged 551 W·m−2 at noon and had a slightly skewed distribution, with slightly more sunlight in the morning than afternoon (Figure 2f). Based on parallel measurements at the same location, photosynthetic photon flux density (PPFD, not shown) was linearly correlated with irradiance according to
PPFD (μmol·m−2·s−1) = 7.88 + 2.49 × Irradiance (W·m−2)
(Least-squares linear regression, N = 1889, R2 = 0.99, P < 0.0001). This suggests that either variable would be equally suitable for predictive modeling. In the current study, stepwise linear regression indicated that the normalized soil CO2 efflux was negatively related to incident sunlight received two days earlier, and positively related to sunlight received 5.5 days previously:
Normalized[ln(Rsoil)] = 1.00 − 0.0000204 × PPFD(−2.0 days) + 0.0000209 × PPFD(−5.5 days)
(stepwise linear regression, selected for minimum BIC (Bayesian Information Criterion), N = 1889, R2 = 0.02, P < 0.0001). The three parameters in this equation were significant at P < 0.0001, 0.0015 and 0.0025, respectively. The predictive utility of this relationship is virtually nil, but a 5.5-day lag between photosynthesis and soil-CO2 emissions is consistent with a 4 to 5 day lag for forest trees [41].

3.2. Seasonal Variability in Rsoil and Environment

Soil respiration during across-site surveys ranged from 0.85 to 46.1 μmol·m−2·s−1 (N = 4858). Fluxes varied significantly among species (Kruskal–Wallis χ2 < 0.0001, Welch’s P < 0.0001): they were lowest in plantations of Pinus, at 3.8 μmol·m−2·s−1; and highest in plantations of Vochysia spp., at 5.3 μmol·m−2·s−1 (Table 2). Fluxes varied significantly among months in all six plantation types (N = 274–1040, Welch’s P < 0.0001), consistent with Hypothesis 2. In all six plantation types, Rsoil was highest in either August or September, and lowest in January or February (Table 3, Table S4). That is, there were remarkably consistent seasonal patterns of soil respiration across the six different plantation types.
Soils were wettest, on average, in November, when 53% of the total soil volume was water, and driest (28%) in April, at the end of La Selva’s drier season (Table 3). Soil moisture content varied significantly among months within all six plantation types (Welch’s tests, P < 0.0001), and the correlation between mean monthly Rsoil and mean monthly soil moisture content was significantly negative (N = 72, Pearson’s r = −0.44, P < 0.0001). Soil temperature also varied significantly among months, based on measurements at both 5-cm and 10-cm depth (Welch’s tests, P < 0.0001). Soil temperatures were coolest in January or February, and were warmest in August or September, depending upon the depth of measurement. Thus, the seasonality observed in Tsoil matched that observed in Rsoil, and soil temperature and Rsoil correlated positively for most species (Table S6). Incident solar radiation is bimodally distributed at the study site, peaking in April and September each year, with a wintertime minimum (Table 3). Such seasonal variations in climate and microclimate potentially drive seasonal variability in Rsoil.

3.3. Litterfall and Leaf Fall

Total fine litterfall over the course of this study ranged from an average of 2.13 ± 0.09 g·m−2·day−1 in plantations of Virola to 3.12 ± 0.17 g·m−2·day−1 in plantations of V. ferruginea. Across all plantation types and years, both total and canopy leaf litter production varied by species, by month, and with a significant species × month interaction (GLM, N = 1072, χ2 < 0.0001). Canopy leaf fall varied significantly through the year in all species except Vochysia ferruginea. Canopy leaves made up an average of 53% of the total fine litterfall and largely defined the seasonality of leaf fall; other litterfall components were not strongly seasonal. Across all species and plots, canopy leaf fall had a bimodal distribution that enveloped the vernal and autumnal equinoxes (Table 3, Table S7). Contrary to Hypothesis 3, which posited that canopy leaf fall and Rsoil would be negatively correlated, mean monthly canopy leaf fall and mean monthly Rsoil were positively correlated (Figure 3, N = 72, Pearson’s r = 0.42, P < 0.001).

3.4. Longer-Term Trends in Rsoil

Soil respiration in plantations of Heironyma (N = 1486, χ2 < 0.0001), Pinus (N = 846, χ2 < 0.001) and Virola (N = 1493, χ2 < 0.0004) increased significantly over the five-year duration of this study, with slopes (± SE) of 0.13 ± 0.026, 0.079 ± 0.028 and 0.054 ± 0.031 μmol·m−2·s−1·year−1, respectively (Figure 4). There were no discernible influences of time, across years, on Rsoil in any of the other plantation types.

4. Discussion

The spatial variability of Rsoil at the study site was high: coefficients of variation (CV) among fluxes gathered over individual across-site surveys averaged 37% and, in some cases, exceeded 100% (Table S2). Previous studies at this site demonstrated differences in Rsoil among tree species. Those differences have persisted through time. For instance, over 2004–2010, Rsoil in plantations of Vochysia spp. averaged 30% greater than emissions in Virola plantations (Table 2), as they did after years one [65] and two [53]. Thus, some of the spatial variability encountered during this study resulted directly from the experimental treatment. Temporal variability was significant but somewhat less important quantitatively: CVs of flux measurements from individual chambers that were sequentially monitored every hour at a single location averaged 18% (Table S1). It was apparent in the field, while measurements were being collected, that the activities of fauna, particularly ants, also generated substantial spatial flux variability, as demonstrated elsewhere [69]. Despite this variability, demonstrable differences were apparent.
The principal objective of this study was to characterize temporal variability in Rsoil within evergreen tree plantations in a lowland tropical rainforest environment. Soil respiration varied significantly on diel, monthly and, in three of six species, multi-annual time frames (Figure 2, Table 3 and Figure 4, respectively). Temporal variations in soil respiration have been documented at other sites [38,49,70,71,72]. They commonly are attributed to temporal variations in plant productivity, or to environmental characteristics such as prevailing temperatures and water availability.

4.1. Diel Variability

The diel variability of Rsoil in this study was not pronounced (Figure 1), but significant differences among hours of the day were observed (χ2 < 0.0001). Specifically, lower fluxes were observed between 06:30 and 10:30, with sunrise occurring at about 06:00. One important reason for measuring hourly variations in Rsoil is to determine if daytime-only measurements cause bias in the estimation of annual fluxes. In this study, more than 99% of the across-site survey measurements of Rsoil were made between 07:30 and 16:30, over which Rsoil averaged 4.12 μmol·m−2·s−1. That value is 2% greater than the morning minima, and 2% lower than the 24-h average flux of 4.21 μmol·m−2·s−1, based on hourly measurements. In short, there was no evidence that daytime measurements of soil respiration generated biased estimates of monthly or seasonal soil respiration from the plantations studied, as a result of diel flux variability.
The notable differences between daytime and nighttime ambient conditions in tropical rainforest environments (e.g., Figure 2) provide the opportunity to apply diel measurements to identify environmental controls over Rsoil. In this study, hourly variations in Rsoil correlated with hourly variations in soil temperatures (Figure 5a, N = 24, Pearson’s r = 0.82, P < 0.0001), with high probability (χ2 < 0.0001) but low predictability (R2 = 0.02). This also was true at the monthly scale (Figure 5b, N = 72, r = 0.50, P < 0.0001). Warmer soils had greater Rsoil. This is consistent with our understanding that rates of metabolism and gas diffusion both increase with temperature. Soil respiration also varied with soil water content (Table S5). However, soil moisture content was not monitored continuously during the diel studies, so sample sizes were too low to evaluate meaningfully. Other considered factors did not correlate with Rsoil, including concurrent atmospheric humidity, ambient CO2 concentrations, the temperature difference between the air and soil, or solar radiation, all of which varied consistently across days (Figure 2).
In contrast to the relatively predictable diel dynamics of many physical factors (Figure 2), precipitation at La Selva occurs at all hours of the day and night (Figure 6). That simple observation may explain the great within-hour variability in CO2 emissions that was observed (Figure 1). Rainfall during this study averaged 4420 mm·year−1. Annual canopy interception of rain may be as high as 710 mm [55], and transpiration is perhaps 1500 mm. Every time that 1 cm of rain falls, 442 times a year, about 8 mm of water enters the soil, where it advectively displaces 8 liters per m−2 of CO2-rich soil atmosphere. The soil atmosphere, in other words, is awash with rainwater, more frequently than it is connected by air-filled pores through which soil gases might continuously diffuse. Every time that 1 cm of soil water is taken up by vegetation and transpired, 150 times each year, 10 liters of atmosphere are pulled into each m2 of soil, refreshing it with oxygen. Water movemet into and out of the soil at La Selva occurs all all hours of the day and, I submit, obfuscates the temperature- and CO2-concentration gradient-driven processes that typically control gas diffusion from soils. Carbon dioxide still escapes from the soil, at a very high rate, but likely through advection as frequently as diffusion, and when it can rather than when it is produced. Animal tunnels offer routes of escape.
Data from the diel study contained evidence of a time lag between PPFD and Rsoil (Equation (2)). Such a lag is expected [73], if PPFD is a meaningful surrogate for canopy photosynthesis. The estimated time lag between canopy fixation of CO2 and CO2 release from soils is 4–5 days for forest trees [41], but there are few data from tropical forests. Additional studies are required to verify that this result was not a Type 1 statistical error: the model included 30 different time lags for selection; two were found to be highly significant, and of opposite direction.
The time series that emerged from each diel measurement period further provide insight into processes. Across days without rain (Figure 7a,b) there are distinct diel variations in soil temperature that are not necessarily associated with parallel changes in Rsoil. Figure 7a shows the best example of what predictable diel variability in Rsoil would look like: peak emissions were offset from peak soil temperatures by about 13 h, and were either 4 h before or 20 h after mid-day. This pattern was not again observed, however. Moderate diel variations in soil temperature in Pentaclethra (Figure 7b) did not apparently influence Rsoil on 28 March. A sharp decline in emissions at 14:00 that coincided with a rainstorm was observed in the time series retrieved from Virola (Figure 7c). The time series from Vochysia guatemalensis (Figure 7d) shows a typically rainy day. Rain seems to temporarily seal the soil surface with a skin of wet soil, impeding the outward diffusion of CO2, but when the rain stops, emissions slowly increase, as the water penetrates the soil. A rainforest is defined by rain and, it seems, so are its soil-CO2 emission rates.

4.2. Monthly Variability

Monthly variability in Rsoil was pronounced: mean soil respiration in September was 60% greater than it was in January (Table 3). Mean monthly fluxes correlated negatively with mean monthly soil moisture content (N = 71, Pearson’s r = −0.45, P < 0.0001), demonstrating that Rsoil was greater in drier months. Emissions correlated positively with both mean monthly soil temperature (Figure 5b N = 72, r = 0.50, P < 0.0001) and irradiance (W·m−2) (N = 72, r = 0.29, P < 0.02). Those findings indicate that monthly weather data could be usefully applied to estimate monthly Rsoil at this site. They do not, however, include any direct information about plant-to-soil carbon fluxes, which fundamentally control the amount of carbon available for respiration.
A specific objective of this study was to seek a plant-growth variable that would provide insight into aboveground–belowground interactions. Hypothesis 3 posited that Rsoil and canopy leaf flushing would be negatively correlated, if carbon allocation aboveground, for the production of a new crop of leaves, was temporally offset from carbon allocation belowground, to produce fine roots and mycorrhizae. The fall of canopy tree leaves is an indirect measure of the timing of leaf flushing. Across all six plantation types, maximum canopy leaf fall occurred in March, in the dry season, before the coming of more dependable rains in May (Table 3 and Table S7). They also shared a secondary peak of canopy leaf fall in October, before the colder, darker, and wetter winter months set in. Over this study, including all species and blocks, canopy leaf fall and soil respiration were not offset from one another, they paralleled one another (N = 72, Pearson’s r = 0.42, P = 0.0002). Fortunately, the phenological characteristics of La Selva’s forest trees have been monitored. The number of tree species flushing many new leaves is maximum in February and March, with a secondary peak in September; and the number of species producing mature fruits is maximum in September [74]. Those observations support the inference that March and September are periods of high allocation of carbon to aboveground processes. They are also months having high soil respiration. Thus, the evidence gathered in this study strongly suggests that above- and belowground C use covaried temporally at this site, opposite of the a priori conjecture of Hypothesis 3. More evidence addressing this topic is warranted, but the alternative hypothesis, i.e., that roots and leaves are not competing with one another for photosynthate but, rather, together depend upon recently produced photosynthates, is reasonable and consistent with the data. Covariance between litterfall and soil respiration was previously found in three forest types in subtropical China [49].

4.3. Physical and Biological Factors

Evidence for a significant time lag between irradiance and soil-CO2 emissions was found (Equation 2), which suggests the possibility of a time lag between canopy photosynthesis and Rsoil. This finding is uncertain because of the large number of independent variables included in the model, and so deserves to be investigated further. The hypothesis that there would be a trade-off in C use between aboveground and belowground C use was tested and rejected. Positive relationships between soil temperatures and Rsoil were observed at diel and monthly time steps (Figure 5), and are evident in the entire dataset derived from across-site surveys (Figure 8a). A significantly negative relationship between soil moisture content and Rsoil was observed at the monthly scale (r = −0.44, P < 0.0001), which indicates that Rsoil was greater in drier months. Across all measurements, Rsoil (μmol·m−2·s−1) varied non-linearly with surface-soil moisture content (SoilH2O, cm·m−1) based on a quadratic polynomial fit:
ln(Rsoil) = 2.036 − 0.01075 × SoilH2O − 0.00037758 × (SoilH2O − 44.12)2
(Figure 8b, N = 4486, F2,4483 = 415, R2 = 0.16) with all three parameters being significant at P < 0.0001. Relationships between Rsoil and soil water content were observed for each of the six species, and were non-linear in all cases except for V. ferruginea (Table S5). At La Selva, soil respiration was limited by a lack of drier soils. Overall, these data support the conclusion that physical variables—soil moisture content, temperatures, solar radiation—can serve as useful inputs into models that are designed to evaluate the likely directional changes in Rsoil that would accompany changing climatic conditions. The absolute magnitude of fluxes was species-dependent, however.

5. Conclusions

In this study, the different tree plantations varied from 3.8 to 5.3 μmol·m−2·s−1 in mean annual emissions, i.e., by 40% (Table 2) at one site, with all plots exposed to the same weather and on the same soils. It is frightening to consider the global ramifications. We need to improve our capacity to translate plant (and animal) species effects to broader spatial scales in a way that could inform physical climate models. This study did not provide new insights into that issue.
It is sometimes put forth that modeling Rsoil is not particularly worthwhile (but see [75]), because Rsoil derives from different processes that respond differently to environmental conditions [76], and often is considered to be the sum of soil heterotrophic and autotrophic respiration [77]. There is merit in modeling those components, and perhaps others, such as roots as distinct from mycorrhizae, and soil arthropods and annelids as distinct from microbes, but we require all the information we can muster, to broadly advance understanding. The value of soil respiration measurements and models is that they reflect an empirically measurable and important soil-to-atmosphere CO2 flux. At a minimum, they provide useful information for testing models of subcomponents. The problem with Rsoil measurements is that they reflect the composite emissions of multiple within-soil CO2-producing sources that are individually influenced by multiple environmental factors.

Supplementary Materials

The following are available online at www.mdpi.com/1999-4907/8/2/40/s1, Table S1: Summary of diel measurements included in this study, Table S2: Summary of across-site survey measurements included in this study, Table S3: Annual climate data for the study site during this study, Table S4: Average soil respiration differed among months within species but in parallel among species, Table S5: Statistical relationships between observed Rsoil and soil water content based on data from the across-site survey study, Table S6: Linear relationships between observed Rsoil and soil temperature at 5 cm depth, Table S7: Mean monthly litterfall across 24 monodominant plantations of six tree species in the Atlantic lowlands of Costa Rica [78].

Acknowledgments

This work was funded by the U.S. National Science Foundation grants DEB-0343766 and DEB-0703561. The work benefitted substantially by the contributions of many, most notably Ricardo Bedoya Arrieta, Oscar Valverde-Barrantes, Dennes Chavarría, and Marlón Hernández, and from the staff and facilities of the La Selva Biological Station.

Conflicts of Interest

The author declares no conflict of interest. The funding agency had no role in the design of the study; in the collection, analyses, or interpretation of data; or in the preparation of the manuscript.

Abbreviations

The following abbreviations are used in this manuscript:
ANOVAanalysis of variance
BICBayesian Information Criterion
CVcoefficient of variation
GLMgeneral linear modeling statistical procedure
IRGAinfrared gas analyzer
PPFDphotosynthetic photon flux density
Rsoilsoil respiration
Tsoiltemperature of the soil

References

  1. Lieth, H. Primary production: Terrestrial ecosystems. Hum. Ecol. 1973, 1, 303–332. [Google Scholar] [CrossRef]
  2. Whittaker, R.H. Communities and Ecosystems, 2nd ed.; Macmillan Publishing Co., Inc.: New York, NY, USA, 1975; p. 385. [Google Scholar]
  3. Erwin, T.L. Tropical forests: Their richness in Coleoptera and other arthropod species. Coleopt. Bull. 1982, 36, 74–75. [Google Scholar]
  4. Dirzo, R.; Raven, P.H. Global state of biodiversity and loss. Annu. Rev. Environ. Resour. 2003, 28, 137–167. [Google Scholar] [CrossRef]
  5. Willig, M.R.; Kaufman, D.M.; Stevens, R.D. Latitudinal lgradients of biodiversity: Pattern, process, scale, and synthesis. Annu. Rev. Ecol. Evol. Syst. 2003, 34, 273–309. [Google Scholar] [CrossRef]
  6. Brown, J.H. Why are there so many species in the tropics. J. Biogeogr. 2014, 41, 8–22. [Google Scholar] [CrossRef] [PubMed]
  7. Saugier, B.; Roy, J.; Mooney, H.A. Estimations of global terrestrial productivity: Converging towards a single number? In Terrestrial Global Productivity; Roy, J., Saugier, B., Mooney, H.A., Eds.; Academic Press: San Diego, CA, USA, 2001; pp. 543–557. [Google Scholar]
  8. Raich, J.W.; Lambers, H.; Oliver, D.J. Respiration in Terrestrial Ecosystems. In Treatise on Geochemistry, 2nd ed.; Holland, H.D., Turekian, K.K., Eds.; Elsevier: Oxford, UK, 2014; pp. 613–649. [Google Scholar]
  9. Schlesinger, W.H.; Jasechko, S. Transpiration in the global water cycle. Agric. For. Meteorol. 2014, 189–190, 115–117. [Google Scholar] [CrossRef]
  10. Dixon, R.K.; Brown, S.; Houghton, R.A.; Solomon, A.M.; Trexler, M.C.; Wisniewski, J. Carbon pools and flux of global forest ecosystems. Science 1994, 263, 185–190. [Google Scholar] [CrossRef] [PubMed]
  11. Jobbágy, E.G.; Jackson, R.B. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol. Appl. 2000, 10, 423–436. [Google Scholar] [CrossRef]
  12. Raich, J.W.; Russell, A.E.; Kitayama, K.; Parton, W.J.; Vitousek, P.M. Temperature influences carbon accumulation in moist tropical forests. Ecology 2006, 87, 76–87. [Google Scholar] [CrossRef] [PubMed]
  13. Hansen, M.C.; Potapov, P.V.; Moore, R.; Hancher, M.; Turubanova, S.A.; Tyukavina, A.; Thau, D.; Stehman, S.V.; Goetz, S.J.; Loveland, T.R.; et al. High-resolution global maps of 21st-century forest cover change. Science 2013, 342, 850–853. [Google Scholar] [CrossRef] [PubMed]
  14. Intergovernmental Panel on Climate Change (IPCC). Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; p. 1535. [Google Scholar]
  15. Arrhenius, S. On the Influence of Carbonic Acid in the Air upon the Temperature of the Ground. Philos. Mag. 1896, 41, 237–276. [Google Scholar] [CrossRef]
  16. Schleser, G.H. The response of CO2 evolution from soils to global temperature changes. Z. Naturforsch. 1982, 37, 287–291. [Google Scholar] [CrossRef]
  17. Jenkinson, D.S.; Adams, D.E.; Wild, A. Model estimates of CO2 emissions from soil in response to global warming. Nature 1991, 351, 304–306. [Google Scholar] [CrossRef]
  18. Kirschbaum, M.U.F. The temperature dependence of soil organic matter decomposition, and the effect of global warming on soil organic C storage. Soil Biol. Biochem. 1995, 27, 753–760. [Google Scholar] [CrossRef]
  19. Kirschbaum, M.U.F. Soil respiration under prolonged soil warming: Are rate reductions caused by acclimation or substrate loss? Glob. Chang. Biol. 2004, 10, 1870–1877. [Google Scholar] [CrossRef]
  20. Davidson, E.A.; Janssens, I.A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 2006, 440, 165–173. [Google Scholar] [CrossRef] [PubMed]
  21. Hartley, I.P.; Heinemeyer, A.; Ineson, P. Effects of three years of soil warming and shading on the rate of soil respiration: Substrate availability and not thermal acclimation mediates observed response. Glob. Chang. Biol. 2007, 13, 1761–1770. [Google Scholar] [CrossRef]
  22. Hopkins, F.M.; Torn, M.S.; Trumbore, S.E. Warming accelerates decomposition of decades-old carbon in forest soils. Proc. Natl. Acad. Sci. USA 2012, 109, E1753–E1761. [Google Scholar] [CrossRef] [PubMed]
  23. Lu, M.; Zhou, X.; Yang, Q.; Li, H.; Luo, Y.; Fang, C.; Chen, J.; Yang, X.; Li, B. Responses of ecosystem carbon cycle to experimental warming: A meta-analysis. Ecology 2013, 94, 726–738. [Google Scholar] [CrossRef] [PubMed]
  24. Knorr, W.; Prentice, I.C.; House, J.I.; Holland, E.A. Long-term sensitivity of soil carbon turnover to warming. Nature 2005, 433, 298–301. [Google Scholar] [CrossRef] [PubMed]
  25. Schindlbacher, A.; Wunderlich, S.; Borken, W.; Kitzler, B.; Zechmeister-Boltenstern, S.; Jandl, R. Soil respiration under climate change: Prolonged summer drought offsets soil warming effects. Glob. Chang. Biol. 2012, 18, 2270–2279. [Google Scholar] [CrossRef]
  26. Cox, P.M.; Pearson, D.; Booth, B.B.; Friedlingstein, P.; Huntingford, C.; Jones, C.D.; Luke, C.M. Sensitivity of tropical carbon to climate change constrained by carbon dioxide variability. Nature 2013, 494, 341–344. [Google Scholar] [CrossRef] [PubMed]
  27. Huntingford, C.; Zelazowski, P.; Galbraith, D.; Mercado, L.M.; Sitch, S.; Fisher, R.; Lomas, M. Simulated resilience of tropical rainforests to CO2-induced climate change. Nat. Geosci. 2013, 6, 268–273. [Google Scholar] [CrossRef]
  28. Bradford, M.A.; Wieder, W.R.; Bonan, G.B.; Fierer, N.; Raymond, P.A.; Crowther, T.W. Managing uncertainty in soil carbon feedbacks to climate change. Nat. Clim. Chang. 2016, 6, 751–758. [Google Scholar] [CrossRef]
  29. Carey, J.C.; Tang, J.; Templer, P.H.; Kroeger, K.D.; Crowther, T.W.; Burton, A.J.; Dukes, J.S.; Emmett, B.; Frey, S.D.; Heskel, M.A.; et al. Temperature response of soil respiration largely unaltered with experimental warming. Proc. Natl. Acad. Sci. USA 2016, 113, 13797–13802. [Google Scholar] [CrossRef] [PubMed]
  30. Vitousek, P.M.; Aplet, G.H.; Raich, J.W.; Lockwood, J.P. Biological perspectives on Mauna Loa Volcano: A model system for ecological research. In Mauna Loa Revealed: Structure, Composition, History, and Hazards; Rhodes, J.M., Lockwood, J.P., Eds.; American Geophysical Union: Washington, DC, USA, 1995; pp. 117–126. [Google Scholar]
  31. Wang, X.; Piao, S.; Ciais, P.; Friedlingstein, P.; Myneni, R.B.; Cox, P.; Heimann, M.; Miller, J.; Peng, S.; Wang, T.; et al. A two-fold increase of carbon cycle sensitivity to tropical temperature variations. Nature 2014, 506, 212–215. [Google Scholar] [CrossRef] [PubMed]
  32. Bahn, M.; Janssens, I.A.; Reichstein, M.; Smith, P.; Trumbore, S.E. Soil respiration across scales: Towards an integration of patterns and processes. New Phytol. 2010, 186, 292–296. [Google Scholar] [CrossRef] [PubMed]
  33. Vargas, R.; Allen, M.F. Environmental controls and the influence of vegetation type, fine roots and rhizomorphs on diel and seasonal variation in soil respiration. New Phytol. 2008, 179, 460–471. [Google Scholar] [CrossRef] [PubMed]
  34. Zimmermann, M.; Meir, P.; Bird, M.; Malhi, Y.; Cahuana, A. Litter contribution to diurnal and annual soil respiration in a tropical montane cloud forest. Soil Biol. Biochem. 2009, 41, 1338–1340. [Google Scholar] [CrossRef]
  35. Zimmermann, M.; Meir, P.; Bird, M.I.; Malhi, Y.; Ccahuana, A.J.Q. Temporal variation and climate dependence of soil respiration and its components along a 3000 m altitudinal tropical forest gradient. Glob. Biogeochem. Cycles 2010, 24, GB4012. [Google Scholar] [CrossRef]
  36. Tan, Z.-H.; Zhang, Y.-P.; Liang, N.; Song, Q.-H.; Liu, Y.-H.; You, G.-Y.; Li, L.-H.; Yu, L.; Wu, C.-S.; Lu, Z.-Y.; et al. Soil respiration in an old-growth subtropical forest: Patterns, components, and controls. J. Geophys. Res. Atmos. 2013, 118, 2981–2990. [Google Scholar] [CrossRef]
  37. Wood, T.E.; Detto, M.; Silver, W.L. Sensitivity of soil respiration to variability in soil moisture and temperature in a humid tropical forest. PLoS ONE 2013, 8, e80965. [Google Scholar] [CrossRef] [PubMed]
  38. Hanpattanakit, P.; Leclerc, M.Y.; Mcmillan, A.M.S.; Limtong, P.; Maeght, J.-L.; Panuthai, S.; Inubushi, K.; Chidthaisong, A. Multiple timescale variations and controls of soil respiration in a tropical dry dipterocarp forest, western Thailand. Plant Soil 2015, 390, 167–181. [Google Scholar] [CrossRef]
  39. Kuzyakov, Y.; Cheng, W. Photosynthesis controls of rhizosphere respiration and organic matter decomposition. Soil Biol. Biochem. 2001, 33, 1915–1925. [Google Scholar] [CrossRef]
  40. Moyano, F.E.; Kutsch, W.L.; Rebmann, C. Soil respiration fluxes in relation to photosynthetic activity in broad-leaf and needle-leaf forest stands. Agric. For. Meteorol. 2008, 148, 135–143. [Google Scholar] [CrossRef]
  41. Kuzyakov, Y.; Gavrichkova, O. Time lag between photosynthesis and carbon dioxide efflux from soil: A review of mechanisms and controls. Glob. Chang. Biol. 2010, 16, 3386–3406. [Google Scholar] [CrossRef]
  42. Ohashi, M.; Kumagai, T.; Kume, T.; Gyokusen, K.; Saitoh, T.M.; Suzuki, M. Characteristics of soil CO2 efflux variability in an aseasonal tropical rainforest in Borneo Island. Biogeochemistry 2008, 90, 275–289. [Google Scholar] [CrossRef]
  43. Curiel Yuste, J.; Janssens, I.A.; Carrara, A.; Ceulemans, R. Annual Q10 of soil respiration reflects plant phenological patterns as well as temperature sensitivity. Glob. Chang. Biol. 2004, 10, 161–169. [Google Scholar] [CrossRef]
  44. Easterling, D.R.; Horton, B.; Jones, P.D.; Peterson, T.C.; Karl, T.R.; Parker, D.E.; Salinger, M.J.; Razuvayev, V.; Plummer, N.; Jamason, P.; et al. Maximum and minimum temperature trends for the globe. Science 1997, 277, 364–367. [Google Scholar] [CrossRef]
  45. Clark, D.A.; Piper, S.C.; Keeling, C.D.; Clark, D.B. Tropical rain forest tree growth and atmospheric carbon dynamics linked to interannual temperature variation during 1984–2000. Proc. Natl. Acad. Sci. USA 2003, 100, 5852–5857. [Google Scholar] [CrossRef] [PubMed]
  46. Vincent, L.A.; Peterson, T.C.; Barros, V.R.; Marino, M.B.; Rusticucci, M.; Carrasco, G.; Ramirez, E.; Alves, L.M.; Ambrizzi, T.; Berlato, M.A.; et al. Observed trends in indices of daily temperature extremes in South America 1960–2000. J. Clim. 2005, 18, 5011–5023. [Google Scholar] [CrossRef]
  47. Donat, M.G.; Alexander, L.V. The shifting probability distribution of global daytime and night-time temperatures. Geophys. Res. Lett. 2012, 39, L14707. [Google Scholar] [CrossRef]
  48. Newstrom, L.E.; Frankie, G.W.; Baker, H.G. A new classification for plant phenology based on flowering patterns in lowland tropical rain forest trees at La Selva, CostaRica. Biotropica 1994, 26, 141–159. [Google Scholar] [CrossRef]
  49. Yan, J.; Wang, Y.; Zhou, G.; Zhang, D. Estimates of soil respiration and net primary production of three forests at different succession stages in South China. Glob. Chang. Biol. 2006, 12, 810–821. [Google Scholar] [CrossRef]
  50. Kira, T.; Shidei, T. Primary production and turnover of organic matter in different forest ecosystems of the western Pacific. Jpn. J. Ecol. 1967, 17, 70–87. [Google Scholar]
  51. Ryan, M.G.; Binkley, D.; Fownes, J.H. Age-related decline in forest productivity: Pattern and process. Adv. Ecol. Res. 1997, 27, 213–262. [Google Scholar]
  52. Ryan, M.G.; Binkley, D.; Fownes, J.H.; Giardina, C.P.; Senock, R.S. An experimental test of the causes of forest growth decline with stand age. Ecol. Monogr. 2004, 74, 393–414. [Google Scholar] [CrossRef]
  53. Russell, A.E.; Raich, J.W.; Bedoya, R.; Valverde-Barrantes, O.; González, E. Impacts of individual tree species on carbon dynamics in a moist tropical forest environment. Ecol. Appl. 2010, 20, 1087–1100. [Google Scholar] [CrossRef] [PubMed]
  54. Clark, D.A.; Clark, D.B.; Oberbauer, S.F. Field-quantified responses of tropical rainforest aboveground productivity to increasing CO2 and climatic stress, 1997–2009. J. Geophys. Res. Biogeosci. 2013, 118, 783–794. [Google Scholar] [CrossRef]
  55. Loescher, H.W.; Gholz, H.L.; Jacobs, J.M.; Oberbauer, S.F. Energy dynamics and modeled evapotranspiration from a wet tropical forest in Costa Rica. J. Hydrol. 2005, 315, 274–294. [Google Scholar] [CrossRef]
  56. Hartshorn, G.S.; Hammel, B.E. Vegetation types and floristic patterns. In La Selva Ecology and Natural History of a Neotropical Rain Forest; McDade, L.A., Bawa, K.S., Hespenheide, H.A., Hartshorn, G.S., Eds.; University of Chicago Press: Chicago, IL, USA, 1994; pp. 73–89. [Google Scholar]
  57. McDade, L.A.; Bawa, K.S.; Hespenheide, H.A.; Hartshorn, G.S. La Selva Ecology and Natural History of a Neotropical Rain Forest; The University of Chicago Press: Chicago, IL, USA, 1994; p. 486. [Google Scholar]
  58. Fisher, R.F. Amelioration of degraded rain forest soils by plantations of native trees. Soil Sci. Soc. Am. J. 1995, 59, 544–549. [Google Scholar] [CrossRef]
  59. Haggar, J.P.; Briscoe, C.B.; Butterfield, R.P. Native species: A resource for the diversification of forestry production in the lowland humid tropics. For. Ecol. Manag. 1998, 106, 195–203. [Google Scholar] [CrossRef]
  60. González, J.E.; Fisher, R.F. Growth of native forest species planted on abandoned pasture land in Costa Rica. For. Ecol. Manag. 1994, 70, 159–167. [Google Scholar] [CrossRef]
  61. Sollins, P.; Sancho, F.M.; Mata, R.C.; Sanford, R.L., Jr. Soils and soil process research. In La Selva Ecology and Natural History of a Neotropical Rain Forest; McDade, L.A., Bawa, K.S., Hespenheide, H.A., Hartshorn, G.S., Eds.; University of Chicago Press: Chicago, IL, USA, 1994; pp. 34–53. [Google Scholar]
  62. Kleber, M.; Schwendenmann, L.; Veldkamp, E.; Rößner, J.; Jahn, R. Halloysite versus gibbsite: Silicon cycling as a pedogenetic process in two lowland Neotropical rain forest soils of La Selva, Costa Rica. Geoderma 2007, 138, 1–11. [Google Scholar] [CrossRef]
  63. Russell, A.E.; Raich, J.W.; Valverde-Barrantes, O.J.; Fisher, R.F. Tree species effects on soil properties in experimental plantations in tropical moist forest. Soil Sci. Soc. Am. J. 2007, 71, 1389–1397. [Google Scholar] [CrossRef]
  64. Russell, A.E.; Raich, J.W. Rapidly growing tropical trees mobilize remarkable amounts of nitrogen, in ways that differ surprisingly among species. Proc. Natl. Acad. Sci. USA 2012, 109, 10398–10402. [Google Scholar] [CrossRef] [PubMed]
  65. Valverde-Barrantes, O.J. Relationships among litterfall, fine root growth, and soil respiration for five tropical tree species. Can. J. For. Res. 2007, 37, 1954–1965. [Google Scholar] [CrossRef]
  66. Raich, J.W.; Russell, A.E.; Bedoya-Arrieta, R. Lignin and enhanced litter turnover in tree plantations of lowland Costa Rica. For. Ecol. Manag. 2007, 239, 128–135. [Google Scholar] [CrossRef]
  67. Haggar, J.P.; Wightman, K.; Fisher, R.F. The potential of plantations to foster woody regeneration within a deforested landscape in lowland Costa Rica. For. Ecol. Manag. 1997, 99, 55–64. [Google Scholar] [CrossRef]
  68. Powers, J.S.; Haggar, J.P.; Fisher, R.F. The effect of overstory composition on understory woody regeneration and species richness in 7-year-old plantations in Costa Rica. For. Ecol. Manag. 1997, 99, 43–54. [Google Scholar] [CrossRef]
  69. Ohashi, M.; Kume, T.; Yamane, S.; Suzuki, M. Hot spots of soil respiration in an Asian tropical rainforest. Geophys. Res. Lett. 2007, 34, L08705. [Google Scholar] [CrossRef]
  70. Oishi, A.C.; Palmroth, S.; Butnor, J.R.; Johsen, K.H.; Oren, R. Spatial and temporal variability of soil CO2 efflux in three proximate temperate forest ecosystems. Agric. For. Meteorol. 2013, 171, 256–269. [Google Scholar] [CrossRef]
  71. Wang, H.; Liu, S.; Wang, J.; Shi, Z.; Lu, L.; Zeng, J.; Ming, A.; Tang, J.; Yu, H. Effects of tree species mixture on soil organic carbon stocks and greenhouse gas fluxes in subtropical plantations in China. For. Ecol. Manag. 2013, 300, 4–13. [Google Scholar] [CrossRef]
  72. Huang, Z.; Yu, Z.; Wang, M. Environmental controls and the influence of tree species on temporal variation in soil respiration in subtropical China. Plant Soil 2014, 382, 75–87. [Google Scholar] [CrossRef]
  73. Högberg, P.; Read, D.J. Towards a more plant physiological perspective on soil ecology. Trends Ecol. Evol. 2006, 21, 548–554. [Google Scholar] [CrossRef] [PubMed]
  74. Frankie, G.W.; Baker, H.G.; Opler, P.A. Comparative phenological studies of trees in tropical wet and dry forests in the lowlands of Costa Rica. J. Ecol. 1974, 62, 881–919. [Google Scholar] [CrossRef]
  75. Phillips, C.L.; Bond-Lamberty, B.; Desai, A.R.; Lavoie, M.; Risk, D.; Tang, J.; Todd-Brown, K.; Vargas, R. The value of soil respiration measurements for interpreting and modeling terrestrial carbon cycling. Plant Soil 2016. [Google Scholar] [CrossRef]
  76. Raich, J.W.; Mora, G. Estimating root plus rhizosphere contributions to soil respiration in annual croplands. Soil Sci. Soc. Am. J. 2005, 69, 634–639. [Google Scholar] [CrossRef]
  77. Savage, K.; Davidson, E.A.; Tang, J. Diel patterns of autotrophic and heterotrophic respiration among phenological stages. Glob. Chang. Biol. 2013, 19, 1151–1159. [Google Scholar] [CrossRef] [PubMed]
  78. Raich, J.W.; Valverde-Barrantes, O.J. Soil CO2 Flux, Moisture, Temperature, and Litterfall at La Selva, Costa Rica, 2003–2010. Available online: http://dx.doi.org/10.3334/ORNLDAAC/1373 (accessed on 7 February 2017).
Figure 1. Soil respiration was highly variable at all hours of the day and night, within all plantation types. This figure includes all data collected during diel measurement periods (N = 2650).
Figure 1. Soil respiration was highly variable at all hours of the day and night, within all plantation types. This figure includes all data collected during diel measurement periods (N = 2650).
Forests 08 00040 g001
Figure 2. Observed diel variability in (a) soil respiration (Rsoil) was significant, but inconsistent (R2 = 0.02), whereas the variability in potentially controlling environmental variables was pronounced. The shaded hours highlight the period between 06:30 and 10:30 when Rsoil was lowest on average; (b) Mean surface air temperatures; (c) soil temperatures at 10-cm depth; (d) the difference between air and soil temperatures; (e) atmospheric relative humidity and (g) ambient atmospheric CO2 concentration were all measured in situ, with Rsoil. Shown are means (±S.E.) of each variable, based on all diel measurement dates and times; (f) Solar irradiance and (h) mean wind speeds were measured at the La Selva weather station, about 3.5 km from the field sites; the data included here were temporally synchronized with the diel measurements.
Figure 2. Observed diel variability in (a) soil respiration (Rsoil) was significant, but inconsistent (R2 = 0.02), whereas the variability in potentially controlling environmental variables was pronounced. The shaded hours highlight the period between 06:30 and 10:30 when Rsoil was lowest on average; (b) Mean surface air temperatures; (c) soil temperatures at 10-cm depth; (d) the difference between air and soil temperatures; (e) atmospheric relative humidity and (g) ambient atmospheric CO2 concentration were all measured in situ, with Rsoil. Shown are means (±S.E.) of each variable, based on all diel measurement dates and times; (f) Solar irradiance and (h) mean wind speeds were measured at the La Selva weather station, about 3.5 km from the field sites; the data included here were temporally synchronized with the diel measurements.
Forests 08 00040 g002
Figure 3. Mean monthly canopy leaf fall and soil respiration (Rsoil) correlated positively across all plantations studied (N = 72, P = 0.0002, r = 0.42). Symbol colors are as in Figure 1. Each symbol is the mean of all measurements from a particular plantation type and month, from 2004 to 2010.
Figure 3. Mean monthly canopy leaf fall and soil respiration (Rsoil) correlated positively across all plantations studied (N = 72, P = 0.0002, r = 0.42). Symbol colors are as in Figure 1. Each symbol is the mean of all measurements from a particular plantation type and month, from 2004 to 2010.
Forests 08 00040 g003
Figure 4. From late 2004 to early 2010, soil respiration (Rsoil) increased significantly in plantations of Hieronyma, Pinus and Virola, but not in other plantation types. Lines show the least-square regressions of significant relationships; among those three species the rate of increase averaged 0.09 μmol·m−2·s−1·year−1.
Figure 4. From late 2004 to early 2010, soil respiration (Rsoil) increased significantly in plantations of Hieronyma, Pinus and Virola, but not in other plantation types. Lines show the least-square regressions of significant relationships; among those three species the rate of increase averaged 0.09 μmol·m−2·s−1·year−1.
Forests 08 00040 g004
Figure 5. (a) Mean hourly soil respiration (Rsoil) in relation to mean hourly soil temperature over the course of the diel study (N = 24, Pearson r = 0.82). (b) Mean monthly soil respiration (Rsoil) in relation to mean monthly daytime soil temperature for all six species over the course of the across-site survey study (N = 72, Pearson’s r = 0.50). Symbol colors follow Figure 1; error bars show ±1 S.E.
Figure 5. (a) Mean hourly soil respiration (Rsoil) in relation to mean hourly soil temperature over the course of the diel study (N = 24, Pearson r = 0.82). (b) Mean monthly soil respiration (Rsoil) in relation to mean monthly daytime soil temperature for all six species over the course of the across-site survey study (N = 72, Pearson’s r = 0.50). Symbol colors follow Figure 1; error bars show ±1 S.E.
Forests 08 00040 g005
Figure 6. (a) The likelihood that rain will occur within any given hour at La Selva is quite uniform across the day. (b) The amount of rain that falls each hour of the day at La Selva is variable: afternoon showers bring more rain. These figures are based on 2009, a relatively typical year.
Figure 6. (a) The likelihood that rain will occur within any given hour at La Selva is quite uniform across the day. (b) The amount of rain that falls each hour of the day at La Selva is variable: afternoon showers bring more rain. These figures are based on 2009, a relatively typical year.
Forests 08 00040 g006
Figure 7. Soil temperatures at 10-cm depth (open symbols), soil respiration (Rsoil, dark red symbols and lines) and hourly rainfall (blue bars) over 52 continuous hours on different dates in plantations of four different tree species at La Selva. Dates refer to the beginning of the measurement period. (a) Hieronyma alchorneoides, in this case, exhibited increasing rates of Rsoil to morning peaks that were distinctly offset from maximum soil temperatures, on two consecutive rain-free days; (b) Pentaclethra macroloba showed no obvious trends in Rsoil across two consecutive rain-free days, despite obvious variations in soil temperatures; (c) A rainstorm in the early afternoon on day two resulted in a dramatic decline in Rsoil in a plantation of Virola koschnyi. Emissions slowly increased following the rain; (d) A typical day in a plantation of Vochysia guatemalensis: rainfall impedes CO2 escape from the soil, and emissions increase after the rain stops.
Figure 7. Soil temperatures at 10-cm depth (open symbols), soil respiration (Rsoil, dark red symbols and lines) and hourly rainfall (blue bars) over 52 continuous hours on different dates in plantations of four different tree species at La Selva. Dates refer to the beginning of the measurement period. (a) Hieronyma alchorneoides, in this case, exhibited increasing rates of Rsoil to morning peaks that were distinctly offset from maximum soil temperatures, on two consecutive rain-free days; (b) Pentaclethra macroloba showed no obvious trends in Rsoil across two consecutive rain-free days, despite obvious variations in soil temperatures; (c) A rainstorm in the early afternoon on day two resulted in a dramatic decline in Rsoil in a plantation of Virola koschnyi. Emissions slowly increased following the rain; (d) A typical day in a plantation of Vochysia guatemalensis: rainfall impedes CO2 escape from the soil, and emissions increase after the rain stops.
Forests 08 00040 g007
Figure 8. (a) The relationship between log (soil respiration), i.e., ln(Rsoil), and soil temperature at 5 cm depth was highly significant, signifying an exponential relationship between CO2 emissions and temperature. Colors follow Figure 1; (b) The relationship between soil moisture content and soil respiration was non-linear (Equation (3)). Both figures include data from all across-site surveys; every other datum is shown to improve legibility.
Figure 8. (a) The relationship between log (soil respiration), i.e., ln(Rsoil), and soil temperature at 5 cm depth was highly significant, signifying an exponential relationship between CO2 emissions and temperature. Colors follow Figure 1; (b) The relationship between soil moisture content and soil respiration was non-linear (Equation (3)). Both figures include data from all across-site surveys; every other datum is shown to improve legibility.
Forests 08 00040 g008
Table 1. Characteristics of the experimental tree plantations in 2004, fifteen years after planting, and in 2010, at the end of this study. Values are means ± 1 S.E. based on measurements of all trees in four replicates of each plantation type.
Table 1. Characteristics of the experimental tree plantations in 2004, fifteen years after planting, and in 2010, at the end of this study. Values are means ± 1 S.E. based on measurements of all trees in four replicates of each plantation type.
Tree SpeciesYear of MeasurementDensity (stems·ha−1)Basal Area (m2·ha−1)Height 1 (m)Diameter 2 (cm)
Hieronyma alchorneoides2004349 ± 1315.2 ± 0.931.1 ± 1.123.6 ± 1.0
2010307 ± 817.4 ±1.032.0 ± 0.626.9 ± 1.1
Pentaclethra macroloba2004611 ± 8721.2 ± 4.421.6 ± 1.721.2 ± 2.9
2010548 ± 5025.3 ± 5.923.3 ± 1.724.0 ± 3.3
Pinus patula2004289 ± 919.8 ± 1.028.9 ± 0.829.6 ± 1.2
2010155 ± 1814.0 ± 2.230.9 ± 2.033.8 ± 1.2
Virola koschnyi2004455 ± 4819.8 ± 3.626.7 ± 1.423.1 ± 1.7
2010432 ± 5022.2 ± 4.227.1 ± 1.425.1 ± 1.8
Vochysia ferruginea2004206 ± 4420.3 ± 2.331.3 ± 0.736.4 ± 2.2
2010156 ± 4519.1 ± 3.731.7 ± 0.741.0 ± 3.4
Vochysia guatemalensis2004519 ± 3040.2 ± 3.633.9 ± 1.731.4 ± 1.3
2010476 ± 3143.6 ± 1.736.2 ± 1.334.3 ± 1.5
1 The mean of the tallest tree in each plot (N = 4); 2 the diameter at breast height of an average-sized tree (i.e., the quadratic mean diameter).
Table 2. Annual average (±1 S.E.) rates of daytime soil respiration (Rsoil) and litterfall in experimental plantations of evergreen trees at La Selva, Costa Rica during this study. Means are weighted by days per month; and standard errors reflect variation among months. Both fluxes varied significantly among the plantation types.
Table 2. Annual average (±1 S.E.) rates of daytime soil respiration (Rsoil) and litterfall in experimental plantations of evergreen trees at La Selva, Costa Rica during this study. Means are weighted by days per month; and standard errors reflect variation among months. Both fluxes varied significantly among the plantation types.
Tree SpeciesRsoil (2004–2010) (μmol·m−2·s−1)Litterfall (2003–2009) (g·m−2·day−1)
Hieronyma alchorneoides4.91 ± 0.642.88 ± 0.21
Pentaclethra macroloba4.15 ± 0.422.65 ± 0.17
Pinus patula3.81 ± 0.842.63 ± 0.14
Virola koschnyi4.08 ± 0.402.13 ± 0.09
Vochysia ferruginea5.26 ± 0.973.12 ± 0.17
Vochysia guatemalensis5.26 ± 0.542.69 ± 0.09
Table 3. Mean monthly soil respiration (Rsoil), soil temperatures (Tsoil), soil water contents, canopy leaf fall and incident solar radiation across 24 monodominant plantations of six tree species in the Atlantic lowlands of Costa Rica. Not all species were measured over the duration of the study; the mean values for each month were determined from the mean monthly values determined for each species. Soil temperatures refer to 2008–2010, when they were measured at 10 cm depth. Error intervals are ±1 S.E.
Table 3. Mean monthly soil respiration (Rsoil), soil temperatures (Tsoil), soil water contents, canopy leaf fall and incident solar radiation across 24 monodominant plantations of six tree species in the Atlantic lowlands of Costa Rica. Not all species were measured over the duration of the study; the mean values for each month were determined from the mean monthly values determined for each species. Soil temperatures refer to 2008–2010, when they were measured at 10 cm depth. Error intervals are ±1 S.E.
MonthRsoil (μmol·m−2·s−1)Tsoil 1 (°C)Soil H2O (cm·m−1)Canopy Leaf Fall (g·m−2·day−1)Mean Daily Irradiance 2 (W·m−2)
January3.67 ± 0.1821.6 ± 0.0451.0 ± 1.81.41 ± 0.14177 ± 43
February3.80 ± 0.2724.0 ± 0.0952.0 ± 1.71.61 ± 0.17154 ± 24
March4.50 ± 0.3022.8 ± 0.0438.1 ± 1.02.01 ± 0.39214 ± 24
April4.73 ± 0.2924.7 ± 0.0431.2 ± 2.61.24 ± 0.19272 ± 23
May4.84 ± 0.3025.2 ± 0.0436.3 ± 1.81.24 ± 0.15232 ± 15
June4.68 ± 0.3125.4 ± 0.0538.8 ± 1.61.27 ± 0.15178 ± 16
July4.31 ± 0.3225.4 ± 0.0344.8 ± 2.01.29 ± 0.17187 ± 22
August4.58 ± 0.4925.5 ± 0.0347.6 ± 1.61.21 ± 0.15236 ± 29
September5.85 ± 0.3326.1 ± 0.0338.6 ± 0.91.60 ± 0.23214 ± 26
October4.97 ± 0.3825.3 ± 0.0350.6± 1.41.78 ± 0.25193 ± 22
November4.12 ± 0.3324.3 ± 0.0353.5 ± 1.01.34 ± 0.22176 ± 43
December4.25 ± 0.3124.0 ± 0.0552.4 ± 2.31.02 ± 0.12159 ± 22
1 January 2008 through February 2010; 2 on dates of Rsoil measurements; data courtesy of La Selva Biological Station.

Share and Cite

MDPI and ACS Style

Raich, J.W. Temporal Variability of Soil Respiration in Experimental Tree Plantations in Lowland Costa Rica. Forests 2017, 8, 40. https://doi.org/10.3390/f8020040

AMA Style

Raich JW. Temporal Variability of Soil Respiration in Experimental Tree Plantations in Lowland Costa Rica. Forests. 2017; 8(2):40. https://doi.org/10.3390/f8020040

Chicago/Turabian Style

Raich, James W. 2017. "Temporal Variability of Soil Respiration in Experimental Tree Plantations in Lowland Costa Rica" Forests 8, no. 2: 40. https://doi.org/10.3390/f8020040

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop