Forest plantations have a large potential for carbon sequestration, playing an important role in the global carbon cycle. However, despite the large amount of research carried out worldwide, the absolute contribution of forest plantations is still incomplete for some parts of the world. To help bridge this gap, we calculated the amount of C stock in three fast growing forest species in Chile. Carbon pools in above-ground and below-ground biomass, forest floor, and soil were considered for this analysis. Across the plantation forests of Chile, carbon accumulated in the above-ground biomass was 181–212 Mg · ha−1 for Pinus radiata, 147–180 Mg · ha−1 for Eucalyptus nitens, and 95–117 Mg · ha−1 for Eucalyptus globulus (age 20–24 years for P. radiata and 10–14 years for Eucalyptus). Total C stocks were for 343 Mg · ha−1 for P. radiata, 352 Mg · ha−1 for E. nitens, and 254 Mg · ha−1 for E. globulus, also at the end of a typical rotation. The carbon pool in the forest floor was found to be significantly lower (less than 4% of the total) when compared to the other pools and showed large spatial variability. Our results agree with other studies showing that 30–50% of the total C stock is stored in the soil. The baseline data will be valuable for modelling C storage changes under different management regimes (changes in species, rotation length and stocking) and for different future climates. Given the contribution of soils to total carbon stocks, special attention should be paid to forest management activities that affect the soil organic carbon pool.

Site Index has been widely used as an age normalised metric in order to account for variation in forest height at a range of spatial scales. Although previous research has used a range of modelling methods to describe the regional variation in Site Index, little research has examined gains that can be achieved through the use of regression kriging or spatial ensemble methods. In this study, an extensive set of environmental surfaces were used as covariates to predict Site Index measurements covering the environmental range of Pinus radiata D. Don plantations in Chile. Using this dataset, the objectives of this research were to (i) compare predictive precision of a range of geostatistical, parametric, and non-parametric models, (ii) determine whether significant gains in precision can be attained through use of regression kriging, (iii) evaluate the precision of a spatial ensemble model that utilises predictions from the five most precise models, through using the model prediction with lowest error for a given pixel, and (iv) produce a map of Site Index across the study area. The five most precise models were all geostatistical and they included ordinary kriging and four regression kriging models that were based on partial least squares or random forests. A spatial ensemble model that was constructed from these five models was the most precise of those developed (RMSE = 1.851 m, RMSE% = 6.38%) and it had relatively little bias. Climatic and edaphic variables were the strongest determinants of Site Index and, in particular, variables that are related to soil water balance were well represented within the most precise predictive models. These results highlight the utility of predicting Site Index using a range of approaches, as these can be used to construct a spatial ensemble that may be more precise than predictions from the constituent models.

Site Index has been widely used as an age normalised metric to account for variation in forest height at a range of spatial scales. Although previous research has used a range of modelling methods to describe regional variation in Site Index little research has examined gains that can be achieved through use of regression kriging or spatial ensemble methods. In this study an extensive set of environmental surfaces were used as covariates to predict Site Index measurements covering the environmental range of \textit{Pinus radiata} D. Don plantations in Chile. Using this dataset, the objectives of this research were to (i) compare predictive precision of a range of geostatistical, parametric and non-parametric models, (ii) determine if significant gains in precision can be attained through use of regression kriging, (iii) evaluate the precision of a spatial ensemble model that utilises predictions from the five most precise models, through using the model prediction with lowest error for a given pixel and (iv) produce a map of Site Index across the study area. The five most precise models were all geostatistical and included ordinary kriging and four regression kriging models that were based on partial least squares or random forests. A spatial ensemble model constructed from these five models was the most precise of those developed (RMSE = 1.851 m, RMSE% = 6.38%) and had relatively little bias. Climatic and edaphic variables were the strongest determinants of Site Index and in particular, variables related to soil water balance were well represented within the most precise predictive models. These results highlight the utility of predicting Site Index using a range of approaches, as these can be used to construct a spatial ensemble that may be more precise than predictions from the constituent models.

Representing the spatial distribution of trees and competition interactions in growth models improves growth prediction and provides insights into spatially explicit forecasts for precise silvicultural interventions. However, this information is rarely taken into account over large areas because obtaining the spatial distribution of individual trees and estimating their competition is both expensive and time consuming. Airborne laser scanning enables rapid estimation of tree height and other attributes over large areas. In this study, we implemented an individual tree detection approach to first extract tree attributes of Pinus radiata D. Don plantations, and second to use this spatially explicit information on tree location and competition to forecast potential tree height, defined as a maximum projected tree height at rotation age. To do so, using a chronosequence of tree heights, we developed a tree height growth model using a Chapman–Richards function, utilizing the effect of inter-tree competition and stand-level top height (TH) on the tree height growth. The results showed that using chronosequence of heights, competition, and TH resulted in accurate predictions of potential tree height (root mean square error = 2.9 m; mean absolute percentage error = 0.154%). We concluded that individual tree height growth is significantly influenced by competition, with increased competition values associated with reductions in potential height growth by 22.2% at 30 years.