Calibrating the self-thinning frontier

Calibration of the self-thinning frontier in even-aged monocultures is hampered by scarce data and by subjective decisions about the proximity of data to the frontier. We present a simple model that applies to observations of the full trajectory of stand mean diameter across a range of densities not necessarily close to the frontier. Development of the model is based on a consideration of the slope s = ln(N_t/N_t−1)/ln(D_t/D_t−1) of a log-transformed plot of stocking N_t and mean stem diameter D_t at time t. This avoids the need for subjective decisions about limiting density and allows the use of abundant data further from the self-thinning frontier. The model can be solved analytically and yields equations for the stocking and the stand basal area as an explicit function of stem diameter. It predicts that self-thinning may be regulated by the maximum basal area with a slope of −2. The significance of other predictor variables offers an effective test of competing self-thinning theories such Yoda's −3/2 power rule and Reineke's stand density index.     

Estimation of biomass and carbon stocks in plants,soil and forest floor in different tropical forests

An accurate characterization of tree carbon (TC), forest floor carbon (FFC) and soil organic carbon (SOC) in tropical forest plantations is important to estimate their contribution to global carbon stocks. This information, however, is poor and fragmented. Carbon contents were assessed in patula pine (Pinus patula) and teak (Tectona grandis) stands in tropical forest plantations of different development stages in combination with inventory assessments and soil survey information. Growth models were used to associate TOC to tree normal diameter (D) with average basal area and total tree height (H_T), with D and H_T parameters that can be used in 6–26 years old patula pine and teak in commercial tropical forests as indicators of carbon stocks. The information was obtained from individual trees in different development stages in 54 patula pine plots and 42 teak plots. The obtained TC was 99.6 Mg ha⁻¹ in patula pine and 85.7 Mg ha⁻¹ in teak forests. FFC was 2.3 and 1.2 Mg ha⁻¹, SOC in the surface layer (0–25 cm) was 92.6 and 35.8 Mg ha⁻¹, 76.1 and 19 Mg ha⁻¹ in deep layers (25–50 cm) in patula pine and teak, respectively. Carbon storage in trees was similar between patula pine and teak plantations, but patula pine had higher levels of forest floor carbon and soil organic carbon. Carbon storage in trees represents 37 and 60% of the total carbon content in patula pine and teak plantations, respectively. Even so, the remaining percentage corresponds to SOC, whereas FFC content is less than 1%. In summary, differences in carbon stocks between patula pine and teak trees were not significant, but the distribution of carbon differed between the plantation types. The low FFC does not explain the SOC stocks; however, current variability of SOC stocks could be related to variation in land use history.     

New dynamic site equation that fits best the Schwappach data for Scots pine (Pinus sylvestris L.) in Central Europe

Using historical growth series data of Scots pine (Pinus sylvestris L.) in Central Europe we examine all the dynamic site equations previously used for modeling the height growth of this species as well as a new dynamic site equation that has not been used previously in the context of this forestry data. The tested models included two groups of anamorphic and polymorphic dynamic site equations (three-dimensional site–height–age models, such as Y = f(t,t₀,y₀)). One group of the models is based on the algebraic difference approach (ADA) implementation of different, preexisting base equations (two-dimensional equations, such as Y = f(t)). The other group of models is based on newer generalized algebraic difference approach (GADA) formulations of new site–height–age relationships that may use older models only as a part of their structure. The models were selected because they were relevant to Scots pine height growth modeling in other studies. We compared all the models with each other in terms of the sum of square deviations associated with fitting them simultaneously to all sites represented by the Scots pine data. All the fits were based on base-age invariant stochastic regressions, in which the global model parameters that are common to all growth series are estimated simultaneously with the site-specific effects that are different for each of the site productivity series. Cieszewski's model [Cieszewski, C.J., 2005. A new flexible GADA based dynamic site equation with polymorphism and variable asymptotes. PMRC Technical Report 2005-2] best described the data.     

Above- and belowground biomass in a Brazilian Cerrado

Cerrado is a biome that occupies about 25% of the Brazilian territory and is characterized by a gradient of grassland to savanna and forest formations and by high species richness. It has been severely affected by degradation and deforestation and has been heavily fragmented over the past 4-5 decades. Despite the recognized overall ecological importance of the Cerrado, there are only few studies focusing on the quantification of biomass in this biome. We conducted such a case study in the South-East of Brazil in a cerrado sensu stricto (cerrado s.s.) with the goal to produce estimates of above- and belowground biomass and to develop allometric equations. A number of 120 trees from 18 species were destructively sampled and partitioned into the components: leaves, branches and bole. Five models with DBH (D), height (H), D²H and wood density (WD) as independent variables were tested for the development of allometric models for individual tree aboveground biomass (leaves + branches + bole). One model based on basal area (BA) as a stand parameter was also tested as an alternative approach for predicting aboveground biomass in the stand level. Belowground biomass was estimated by subsampling on 10 sample plots. Mean aboveground tree biomass (bole, branches and leaves) was estimated to be 62,965.5 kg ha⁻¹(SE = 14.6%) and belowground biomass accounted for 37,501.8 kg ha⁻¹ (SE = 23%). The best-fit equation for the estimation of individual tree aboveground biomass include DBH and wood density as explanatory variables (R² = 0.898; SEE = 0.371) and is applicable for the diameter range of this study (5.0-27.6 cm) and in environments with similar conditions of the cerrado s.s. sampled. In the stand level, the model tested presented a higher goodness of fit than the single tree models (R² = 0.934; SEE = 0.224). Our estimates of aboveground biomass are higher than reported by other studies developed in the same physiognomy, but the estimates of belowground biomass are within the range of values reported in other studies from sites in cerrado s.s. Both biomass estimates, however, exhibit relatively large standard errors. The root-to-shoot ratio of the sample trees is in the magnitude of reported values for savanna ecosystems, but smaller than estimated from other studies in the cerrado s.s.     

Soil–atmosphere CO2, CH4 and N2O fluxes in boreal forestry-drained peatlands

Greenhouse gas emissions from managed peatlands are annually reported to the UNFCCC. For the estimation of greenhouse gas (GHG) balances on a country-wide basis, it is necessary to know how soil–atmosphere fluxes are associated with variables that are available for spatial upscaling. We measured momentary soil–atmosphere CO₂ (heterotrophic and total soil respiration), CH₄ and N₂O fluxes at 68 forestry-drained peatland sites in Finland over two growing seasons. We estimated annual CO₂ effluxes for the sites using site-specific temperature regressions and simulations in half-hourly time steps. Annual CH₄ and N₂O fluxes were interpolated from the measurements. We then tested how well climate and site variables derived from forest inventory results and weather statistics could be used to explain between-site variation in the annual fluxes. The estimated annual CO₂ effluxes ranged from 1165 to 4437 g m⁻² year⁻¹ (total soil respiration) and from 534 to 2455 g m⁻² year⁻¹ (heterotrophic soil respiration). Means of 95% confidence intervals were ±12% of total and ±22% of heterotrophic soil respiration. Estimated annual CO₂ efflux was strongly correlated with soil respiration at the reference temperature (10 °C) and with summer mean air temperature. Temperature sensitivity had little effect on the estimated annual fluxes. Models with tree stand stem volume, site type and summer mean air temperature as independent variables explained 56% of total and 57% of heterotrophic annual CO₂ effluxes. Adding summer mean water table depth to the models raised the explanatory power to 66% and 64% respectively. Most of the sites were small CH₄ sinks and N₂O sources. The interpolated annual CH₄ flux (range: −0.97 to 12.50 g m⁻² year⁻¹) was best explained by summer mean water table depth (r² = 64%) and rather weakly by tree stand stem volume (r² = 22%) and mire vegetation cover (r² = 15%). N₂O flux (range: −0.03 to 0.92 g m⁻² year⁻¹) was best explained by peat CN ratio (r² = 35%). Site type explained 13% of annual N₂O flux. We suggest that water table depth should be measured in national land-use inventories for improving the estimation of country-level GHG fluxes for peatlands.     

An ecosystem approach to biodiversity effects: Carbon pools in a tropical tree plantation

This paper presents a synthesis of experiments conducted in a tropical tree plantation established in 2001 and consisting of 22 plots of 45 m × 45 m with either one, three or six native tree species. We examined the changes in carbon (C) pools (trees, herbaceous vegetation, litter, coarse woody debris (CWD), and mineral topsoil at 0-10 cm depth) and fluxes (decomposition of CWD and litter, as well as soil respiration) both through time and among diversity levels. Between 2001 and 2009 the aboveground C pools increased, driven by trees. Across diversity levels, the mean observed aboveground C pool was 7.9 ± 2.5 Mg ha⁻¹ in 2006 and 20.4 ± 7.4 Mg ha⁻¹ in 2009, a 158% increase. There was no significant diversity effect on the observed aboveground C pool, but we found a significant decrease in the topsoil C pool, with a mean value of 34.5 ± 2.4 Mg ha⁻¹ in 2001 and of 25.7 ± 5.7 Mg ha⁻¹ in 2009 (F_1,36 = 52.12, p < 0.001). Assuming that the biomass C pool in 2001 was negligible (<1 Mg ha⁻¹), then the plantation gained in C, on average, ∼20 and lost ∼9 Mg ha⁻¹ in biomass and soil respectively, for an overall gain of ∼11 Mg ha⁻¹ over 8 years. Across the entire data set, we uncovered significant effects of diversity on CWD decomposition (diversity: F_2,393 = 15.93, p < 0.001) and soil respiration (monocultures vs mixtures: t = 15.35, df = 11, p < 0.05) and a marginally significant time × diversity interaction on the loss of total C from the mineral topsoil pool (see above). Monthly CWD decomposition was significantly faster in monocultures (35.0 ± 24.1%) compared with triplets (31.3 ± 21.0%) and six-species mixtures (31.9 ± 26.8%), while soil respiration was higher in monocultures than in mixtures (t = 15.35, df = 11, p < 0.001). Path analyses showed that, as diversity increases, the links among the C pools and fluxes strengthen significantly. Our results demonstrate that tree diversity influences the processes governing the changes in C pools and fluxes following establishment of a tree plantation on a former pasture. We conclude that the choice of tree mixtures for afforestation in the tropics can have a marked influence on C pools and dynamics.     

Dynamics of wood recruitment in streams of the northeastern US

Wood is an important component of forested stream ecosystems, and stream restoration efforts often incorporate large wood. In most cases, however, stream restoration projects are implemented without information regarding the amount of wood that historically occurred or the natural rates of wood recruitment. This study uses a space-for-time analysis to quantify large wood loading to 28 streams in the northeastern US with a range of in-stream and riparian forest characteristics. We document the current volume and frequency of occurrence of large wood in streams with riparian forests varying in their stage of stand development as well as stream size and gradient. Linear models relating stream wood characteristics to stream geomorphic and forest characteristics were compared using Akaike's Information Criterion (AIC) model selection. The AIC analysis indicated that the volume and frequency of large wood and wood accumulations (wood jams) in streams was most closely associated with the age of the dominant canopy trees in the riparian forest (best models: log₁₀(large wood volume (m³ 100 m⁻¹)) = (0.0036 × stand age) − 0.2281, p < 0.001, r² = 0.80; and large wood frequency (number per 100 m) = (0.1326 × stand age) + 7.3952, p < 001, r² = 0.63). Bankfull width was an important factor accounting for wood volume per unit area (m³ ha⁻¹) but not the volume of wood per length of stream (100 m⁻¹). The empirical models developed in this study were unsuccessful in predicting wood loading in other regions, most likely due to difference in forest characteristics and the legacy of forest disturbance. However, these models may be applicable in other streams in the northeastern US or in streams with comparable riparian forests, underlying geology, and disturbance regimes—factors that could alter long-term wood loading dynamics. Our results highlight the importance of understanding region-specific processes when planning stream restoration and stream management projects.     

Tree biomass estimation in regenerating areas of tropical dry vegetation in northeast Brazil

Allometric equations have been developed for various different vegetation types but have rarely been validated in the field and never for dry tropical forest such as caatinga. In three areas of semi-arid Brazil, with regenerating caatinga vegetation, we measured and weighed twelve hundred individuals of four tree species and used the data to validate equations previously determined in mature caatinga. They and several other equations developed for tropical vegetations overestimate the biomass (B) of trees from the regeneration areas by more than 20%, possibly because these trees have reduced crowns, with lower branch masses. We then determined new allometric equations for them, validating equations for one site against data of the others and pooling the data if they were cross-validated. The best equations were power ones, based on diameter at breast height (D), with little improvement by including height, crown area and/or wood density (Caesalpinia pyramidalis, B = 0.3129D^1.8838; Croton sonderianus, B = 0.4171D^1.5601; Mimosa ophthalmocentra, B = 0.4369D^1.8493; and Mimosa tenuiflora, B = 0.3344D^1.9648 and 0.4138D^1.7718).     

Radiation use efficiency in adjacent hardwood and pine forests in the southern Appalachians

The efficiency with which trees convert photosynthetically active radiation (PAR) to biomass has been shown to be consistent within stands of an individual species, which is useful for estimating biomass production and carbon accumulation. However, radiation use efficiency (?) has rarely been measured in mixed-species forests, and it is unclear how species diversity may affect the consistency of ?, particularly across environmental gradients. We compared aboveground net primary productivity (ANPP), intercepted photosynthetically active solar radiation (IPAR), and radiation use efficiency (? = ANPP/IPAR) between a mixed deciduous forest and a 50-year-old white pine (Pinus strobus L.) plantation in the southern Appalachian Mountains. Average ANPP was similar in the deciduous forest (11.5 Mg ha⁻¹ y⁻¹) and pine plantation (10.2 Mg ha⁻¹ y⁻¹), while ? was significantly greater in the deciduous forest (1.25 g MJ⁻¹) than in the white pine plantation (0.63 g MJ⁻¹). Our results demonstrate that late-secondary hardwood forests can attain similar ANPP as mature P. strobus plantations in the southern Appalachians, despite substantially less annual IPAR and mineral-nitrogen availability, suggesting greater resource-use efficiency and potential for long-term carbon accumulation in biomass. Along a 260 m elevation gradient within each forest there was not significant variation in ?. Radiation use efficiency may be stable for specific forest types across a range of environmental conditions in the southern Appalachian Mountains, and thus useful for generating estimates of ANPP at the scale of individual watersheds.     

Transpiration along an age series of Eucalyptus globulus plantations in southeastern Australia

Many of the world's Eucalyptus plantations are grown on short rotations of 15 years or less, which often covers the most rapid phase of stand development and peaks in growth rates and leaf areas. Since transpiration is related to stand leaf area these short rotations that make use of rapid early growth rates, may also maximise plantation water use, which has implications for predicting their water requirements and impacts on catchment hydrology. This study examined the transpiration, leaf area and growth rates of Eucalyptus globulus Labill. plantations aged 2–8 years. Transpiration (E), estimated using the heat pulse technique, increased from 0.4 mm day⁻¹ at age 2 years to a peak of about 1.6–1.9 mm day⁻¹ in stands aged 5–7 years. This was associated with similar trends for stand leaf area index (LAI) and periodic annual increments of aboveground biomass, which both peaked at about age 4–6 years resulting in a linear relationship between E and LAI. While stand sapwood areas were continuing to increase at age 8 years, E was already declining due to reductions in sap velocity, from 13.5 cm h⁻¹ at age 2 years to 6.3 cm h⁻¹ at age 8 years and reduced sapwood area growth rates. Trees compensated for this reduction in sap velocity with declines in the leaf area (A_L) to sapwood area (A_S) relationship (A_L:A_S) with age. There was also a reduction in growth efficiency (aboveground biomass increment per LAI) with age. However, reductions in WUE were small after age 4 years, which explained the linear relationship between E and LAI. If E continues to decline successive short rotation lengths may not only make use of rapid early growth rates but could also increase plantation water use compared to longer rotations over the same period of time.     

Validation of plantation transpiration in south-eastern Australia estimated using the 3PG+ forest growth model

Forest plantations for wood production are an increasingly important land use in southern Australia, and there are potentially important hydrologic consequences of what is mostly a change in land use from agriculture to silviculture. An ability to predict, with some degree of accuracy, the impact of plantation expansion on surface water and groundwater resources is essential. A validated process-based modelling approach, integrating the many interacting environmental and management factors which may influence plantation growth and transpiration, can be used for this purpose. The 3PG forest growth model has been evaluated for a number of species from widely differing climate and site conditions. While growth predictions have been validated, little attention has been given to testing the accuracy of the transpiration predictions or the model's representation of the water balance. We enhanced the 3PG forest growth model (known as 3PG+) and then integrated it into the Catchment Analysis Tool (CAT), so that it now interfaces with a more detailed multi-layered, daily time step representation of the soil water balance. Simulated transpiration using 3PG+ in CAT was compared with field measurements in 30 plots (across 15 sites) representing 5 common plantation species (Eucalyptus globulus, E. nitens, E. grandis, E. regnans and Pinus radiata) across ages 2–31 years. Mean daily plot transpiration during the measurement periods ranged between 0.4 and 4.2 mm day⁻¹ (average 2.0 mm day⁻¹). Simulated mean daily plot transpiration using 3PG+ in CAT for Eucalyptus was good (coefficient of efficiency = 0.80; R² = 0.81). While the model tended to slightly under-predict transpiration at higher measured rates (>3.5 mm day⁻¹), predictions at monthly timescales had acceptable accuracy. The integration of 3PG+ into CAT resulted in an improvement in accuracy and applicability of CAT, and provides for the spatial application of 3PG+ across diverse and mixed land use catchments for investigation into carbon and water movement in forest systems.     

Formulating allometric equations for estimating biomass and carbon stock in small diameter trees

Biomass and carbon sequestration rate of a young (four year old) mixed plantation of Dalbergia sissoo Roxb., Acacia catechu Willd., and Albizia lebbeck Benth. growing in Terai region (a level area of superabundant water) of central Himalaya was estimated. The plantation is seed sown in the rainy season of year 2004 and spread over an area of 44 ha. Allometric equations for both above and below ground components were developed for three tree species. The density of trees in the plantation was 1322 trees ha⁻¹ The diameters of trees were below 10 cm. Five diameter classes were defined for D. sissoo and A. catechu and 3 for A. lebbeck. 5 trees were harvested in each diameter class. Individual tree allometry was exercised for developing the allometric equations relating tree component (low and above ground) biomass to d.b.h. Post analysis equations were highly significant (P > 0.001) for each component of all species. In the plantation Holoptelia integrifolia Roxb. (Family Ulmaceae) has been reduced to shrub form because of frost. Only the aboveground biomass of H. integrifolia and other shrubs were estimated by destructive harvesting method. Herbaceous forest floor biomass and leaf litter fall were also estimated. The total forest vegetation biomass was 10.86 Mg ha⁻¹ in 2008 which increased to 19.49 Mg ha⁻¹ in 2009. The forest is sequestering carbon at the rate of 4.32 Mg ha⁻¹ yr⁻¹.     

Chloroplast DNA diversity in Castanopsis hystrix populations in south China

Nineteen Chinese populations of Castanopsis hystrix were examined to quantify genetic diversity and genetic structure at chloroplast DNA. Microsatellites (SSR) were analyzed by PCR using conserved primers. The average within population gene diversity (H_S), the total gene diversity (H_T), and the differentiation for unordered alleles (G_ST) and for ordered alleles (N_ST) were measured. Fourteen different haplotypes were detected, two of them very common. The level of differentiation among populations (G_ST = 23.6%) indicates a highly efficient seed dispersal mechanism. In addition, the difference between G_ST and N_ST for the species is not significant, suggesting that the phylogeographic structure is weak or absent. The geographical pattern of C. hystrix haplotypes could be attributed to its migration from the numerous and scattered refugia, where the species confined during the last glacial period. These results provide an important insight into patterns of postglacial recolonization of this tree species.     

Carbon dioxide exchange of a larch forest after a typhoon disturbance

A typhoon event catastrophically destroyed a 45-year-old Japanese larch plantation in southern Hokkaido, northern Japan in September 2004, and about 90% of trees were blown down. Vegetation was measured to investigate its regeneration process and CO₂ flux, or net ecosystem production (NEP), was measured in 2006–2008 using an automated chamber system to investigate the effects of typhoon disturbance on the ecosystem carbon balance. Annual maximum aboveground biomass (AGB) increased from 2.7 Mg ha⁻¹ in 2006 to 4.0 Mg ha⁻¹ in 2007, whereas no change occurred in annual maximum leaf area index (LAI), which was 3.7 m² m⁻² in 2006 and 3.9 m² m⁻² in 2007. Red raspberry (Rubus idaeus) had become dominant within 2 years after the typhoon disturbance, and came to account for about 60% and 50% of AGB and LAI, respectively. In comparison with CO₂ fluxes measured by the eddy covariance technique in 2001–2003, for 4.5 months during the growing season, the sum of gross primary production (GPP) decreased on average by 739 gC m⁻² (64%) after the disturbance, whereas ecosystem respiration (RE) decreased by 501 gC m⁻² (51%). As a result, NEP decreased from 159 ± 57 gC m⁻² to −80 ± 30 gC m⁻², which shows that the ecosystem shifted from a carbon sink to a source. Seasonal variation in RE was strongly correlated to soil temperature. The interannual variation in the seasonal trend of RE was small. Light-saturated GPP (P_max) decreased from 30–45 μmol m⁻² s⁻¹ to 8–12 μmol m⁻² s⁻¹ during the summer season through the disturbance because of large reduction in LAI.     

Forest thinning and subsequent bark beetle-caused mortality in Northeastern California

The Warner Mountains of northeastern California on the Modoc National Forest experienced a high incidence of tree mortality (2001–2007) that was associated with drought and bark beetle (Coleoptera: Curculionidae, Scolytinae) attack. Various silvicultural thinning treatments were implemented prior to this period of tree mortality to reduce stand density and increase residual tree growth and vigor. Our study: (1) compared bark beetle-caused conifer mortality in forested areas thinned from 1985 to 1998 to similar, non-thinned areas and (2) identified site, stand and individual tree characteristics associated with conifer mortality. We sampled ponderosa pine (Pinus ponderosa var ponderosa Dougl. ex Laws.) and Jeffrey pine (Pinus jeffreyi Grev. and Balf.) trees in pre-commercially thinned and non-thinned plantations and ponderosa pine and white fir (Abies concolor var lowiana Gordon) in mixed conifer forests that were commercially thinned, salvage-thinned, and non-thinned. Clusters of five plots (1/50th ha) and four transects (20.1 × 100.6 m) were sampled to estimate stand, site and tree mortality characteristics. A total of 20 pre-commercially thinned and 13 non-thinned plantation plot clusters as well as 20 commercially thinned, 20 salvage-thinned and 20 non-thinned mixed conifer plot clusters were established. Plantation and mixed conifer data were analyzed separately. In ponderosa pine plantations, mountain pine beetle (Dendroctonus ponderosae Hopkins) (MPB) caused greater density of mortality (trees ha⁻¹ killed) in non-thinned (median 16.1 trees ha⁻¹) compared to the pre-commercially thinned (1.2 trees ha⁻¹) stands. Percent mortality (trees ha⁻¹ killed/trees ha⁻¹ host available) was less in the pre-commercially thinned (median 0.5%) compared to the non-thinned (5.0%) plantation stands. In mixed conifer areas, fir engraver beetles (Scolytus ventralis LeConte) (FEN) caused greater density of white fir mortality in non-thinned (least square mean 44.5 trees ha⁻¹) compared to the commercially thinned (23.8 trees ha⁻¹) and salvage-thinned stands (16.4 trees ha⁻¹). Percent mortality did not differ between commercially thinned (least square mean 12.6%), salvage-thinned (11.0%), and non-thinned (13.1%) mixed conifer stands. Thus, FEN-caused mortality occurred in direct proportion to the density of available white fir. In plantations, density of MPB-caused mortality was associated with treatment and tree density of all species. In mixed conifer areas, density of FEN-caused mortality had a positive association with white fir density and a curvilinear association with elevation.     

Simulating age-related changes in carbon storage and allocation in a Chinese fir plantation growing in southern China using the 3-PG model

Chinese fir [(Cunninghamia lanceolata (Lamb.) Hook (Taxodiaceae)] plantations are helping to meet China's increasing demands for timber, while, at the same time, sequestering carbon (C) above and belowground. The latter function is important as a means of slowing the rate that CO₂ is increasing in the atmosphere. Available data are limited, however, and even if extensive, would necessitate consideration of future changes in climatic conditions and management practices. To evaluate the contribution of Chinese fir plantations under a range of changing conditions a dynamic model is required. In this paper, we report successful outcome in parameterizing a process-based model (3-PG) and validating its predictions with recent and long-term field measurements acquired from different ages of Chinese fir plantations at the Huitong National Forest Ecosystem Research Station. Once parameterized, the model performed well when simulating leaf area index (LAI), net primary productivity (NPP), biomass of stems (W_S), foliage (W_F) and roots (W_R), litterfall, and shifts in allocation over a period of time. Although the model does not specifically include heterotrophic respiration, we made some attempts to estimate changes in root C storage and decomposition rates in the litterfall pool as well as in the total soil respiration. Total C stored in biomass increased rapidly, peaking at age 21 years in unthinned stands. The predicted averaged above and belowground NNP (13.81 t ha⁻¹ a⁻¹) of the Chinese fir plantations between the modeling period (from 4 to 21-year-old) is much higher than that of Chinese forests (4.8–6.22 t ha⁻¹ a⁻¹), indicating that Chinese fir is a suitable tree species to grow for timber while processing the potential to act as a C sequestration sink. Taking into account that maximum LAI occurs at the age of 15 years, intermediate thinning and nutrient supplements should, according to model predictions, further increase growth and C storage in Chinese fir stands. Predicted future increases (approximately 0–2 °C) in temperature due to global warming may increase plantation growth and reduce the time required to complete a rotation, but further increases (approximately 2–6 °C) may reduce the growth rate and prolong the rotational age.     

Effects of [CO2] and nitrogen on morphological and biomass traits of white birch (Betula papyrifera) seedlings

To investigate the interactive effects of CO₂ concentration ([CO₂]) and nitrogen supply on the growth and biomass of boreal trees, white birch seedlings (Betula papyrifera) were grown under ambient (360 μmol mol⁻¹) and elevated [CO₂] (720 μmol mol⁻¹) with five nitrogen supply regimes (10, 80, 150, 220, and 290 μmol mol⁻¹) in greenhouses. After 90 days of treatment, seedling height, root-collar diameter, biomass of different organs, leaf N concentration, and specific leaf area (SLA) were measured. Significant interactive effects of [CO₂] and N supply were found on height, root-collar diameter, leaf biomass, stem biomass and total biomass, stem mass ratio (SMR), and root mass ratio (RMR), but not on root mass, leaf mass ratio (LMR), leaf to root ratio (LRR), or leaf N concentration. The CO₂ elevation generally increased all the growth and biomass parameters and the increases were generally greater at higher levels of N supply or higher leaf N concentration. However, the CO₂ elevation significantly reduced SLA (13.4%) and mass-based leaf N concentration but did not affect area-based leaf N concentration. Increases in N supply generally increased the growth and biomass parameters, but the relationships were generally curvilinear. Based on a second order polynomial model, the optimal leaf N concentration was 1.33 g m⁻² for height growth under ambient [CO₂] and 1.52 g m⁻² under doubled [CO₂]; 1.48 g m⁻² for diameter under ambient [CO₂] and 1.64 g m⁻² under doubled [CO₂]; 1.29 g m⁻² for stem biomass under ambient [CO₂] and 1.43 g m⁻² under doubled [CO₂]. The general trend is that the optimal leaf N was higher at doubled than ambient [CO₂]. However, [CO₂] did not affect the optimal leaf N for leaf and total biomass. The CO₂ elevation significantly increased RMR and SMR but decreased LMR and LRR. LMR increased and RMR decreased with the increasing N supply. SMR increased with increase N supply up to 80 μmol mol⁻¹ and then leveled off (under elevated [CO₂]) or stated to decline (under ambient [CO₂]) with further increases in N supply. The results suggest that the CO₂ elevation increased biomass accumulation, particularly stem biomass and at higher N supply. The results also suggest that while modest N fertilization will increase seedling growth and biomass accumulation, excessive application of N may not stimulate further growth or even result in growth decline.     

Acid deposition effects on forest composition and growth on the Monongahela National Forest,West Virginia

The northern and central Appalachian forests are subject to high levels of atmospheric acid deposition (AD), which has been shown in some forests to negatively impact forest growth as well as predispose the forest system to damage from secondary stresses. The purpose of this study was to evaluate the possible contribution of AD to changes in composition and productivity of the Monongahela National Forest, and to evaluate soil-based indicators of acidification that might be useful for detecting AD-related forest changes. Soils adjacent to 30 Forest Inventory and Analysis (FIA) sites were sampled and analyzed for a suite of acidity indicators. These indicators were correlated with the periodic mean annual volume increment (PMAVI) of the forest stands on FIA plots for the 10-yr period 1989–2000. PMAVI ranged from −9.5 to 11.8 m³ ha⁻¹ yr⁻¹, with lower-than-expected growth (<3 m³ ha⁻¹ yr⁻¹) on two-thirds of the sites. In the surface horizon, effective base saturation, Ca²⁺ concentration, base saturation, K⁺ concentration, Ca/Al molar ratio, and Mg/Al molar ratio, were positively correlated with PMAVI and Fe concentration was negatively correlated with PMAVI (p ≤ 0.1). In the subsurface horizon pH_(w) and effective base saturation were positively correlated and Al³⁻ concentration and K⁺ concentration were negatively correlated with PMAVI. We hypothesized that NO₃-N/NH₄-N ratio would also be correlated with PMAVI, but it was not. Correlations between soil chemical indicators and PMAVI suggest that AD may contribute, in part, to the lower-than-expected forest growth on the Monongahela National Forest.     

Ecosystem carbon stocks and distribution under different land-uses in north central Alberta,Canada

Land-use and land cover strongly influence carbon (C) storage and distribution within ecosystems. We studied the effects of land-use on: (i) above- and belowground biomass C, (ii) soil organic C (SOC) in bulk soil, coarse- (250–2000 μm), medium- (53–250 μm) and fine-size fractions (<53 μm), and (iii) ¹³C and ¹⁵N abundance in plant litter, bulk soil, coarse-, and medium- and fine-size fractions in the 0–50 cm soil layer in Linaria AB, Canada between May and October of 2006. Five adjacent land-uses were sampled: (i) agriculture since 1930s, (ii) 2-year-old hybrid poplar (Populusdeltoides × Populus × petrowskyana var. Walker) plantation, (iii) 9-year-old Walker hybrid poplar plantation, (iv) grassland since 1997, and (v) an 80-year-old native aspen (Populus tremuloides Michx.) stand. Total ecosystem C stock in the native aspen stand (223 Mg C ha⁻¹) was similar to that of the 9-year-old hybrid poplar plantation (174 Mg C ha⁻¹) but was significantly greater than in the agriculture (132 Mg C ha⁻¹), 2-year-old hybrid poplar plantation (110 Mg C ha⁻¹), and grassland (121 Mg C ha⁻¹). Differences in ecosystem C stocks between the land-uses were primarily the result of different plant biomass as SOC in the 0–50 cm soil layer was unaffected by land-use change. The general trend for C stocks in soil particle-size fractions decreased in the order of: fine > medium > coarse for all land-uses, except in the native aspen stand where C was uniformly distributed among soil particle-size fractions. The C stock in the coarse-size fraction was most affected by land-use change whilst the fine fractions the least. Enrichment of the natural abundances of ¹³C and ¹⁵N across the land-uses since time of disturbance, i.e., from agriculture to 2- and then 9-year-old hybrid poplar plantations or to grassland, suggests shifts from more labile forms of C to more humified forms of C following those land-use changes.     

Does soil organic matter variation affect the retention of NH4 and NO3 in forest ecosystems?

Human activities have fundamentally changed global nitrogen (N) cycling, leading to elevated N deposition in most parts of the world. The fate of deposited N, whether being retained to sustain plant growth or causing ecosystem N saturation, is critical to the global carbon (C) cycling and local environment. In a short-term laboratory experiment, we used ¹⁵N-labeled NH₄⁺ and NO₃⁻ to study the fate of N inputs in forest soils and what regulates N retention. Soils with a wide range of organic matter content and other attributes were collected from a 70-year-old plantation containing monotypic stands of Norway spruce (Picea abies), red pine (Pinus resinosa), sugar maple (Acer saccharum), and red oak (Quercus rubra), and separated into 0-5 cm and 5-15 cm layers. Nitrogen added to the soil was either immediately extracted (Time 0: T0) with K₂SO₄ solution, or incubated for 7 d (T7) and then extracted. Retention of ¹⁵N into the non-extractable soil pool at T0 was limited; but after the 7-d incubation, between 20 and 70% of the ¹⁵NH₄⁺ was retained. Nitrification transformed on average 50% of the ¹⁵NH₄⁺ into ¹⁵NO₃⁻ during the incubation while retention of ¹⁵NO₃⁻ at T7 remained low (7.40 ± 1.08%). Retention of ¹⁵NH₄⁺ into non-extractable soil at T7 was positively correlated to the percentage of soil organic matter (SOM) (r² = 0.323, P < 0.001), and was significantly higher (P < 0.001) in the high-SOM 0-5 cm layer than in the low-SOM 5-15 cm layer. Conversion of ¹⁵NH₄⁺ to ¹⁵NO₃⁻ during incubation significantly reduced the ¹⁵NH₄⁺ retention (P < 0.001). Our results suggest that the variations of SOM and other soil attributes play strong roles in the retention of newly deposited inorganic N and could affect forest ecosystem responses to chronic N deposition.