Land use and land cover change research has been applied to landslides, erosion, land planning and global change. Based on the CA-Markov model, this study predicts the spatial patterns of land use in 2025 and 2036 based on the dynamic changes in land use patterns using remote sensing and geographic information system. CA-Markov integrates the advantages of cellular automata and Markov chain analysis to predict future land use trends based on studies of land use changes in the past. Based on Landsat 5 TM images from 1992 and 2003 and Landsat 8 OLI images from 2014, this study obtained a land use classification map for each year. Then, the genetic transition probability from 1992 to 2003 was obtained by IDRISI software. Based on the CA-Markov model, a predicted land use map for 2014 was obtained, and it was validated by the actual land use results of 2014 with a Kappa index of 0.8128. Finally, the land use patterns of 2025 and 2036 in Jiangle County were determined. This study can provide suggestions and a basis for urban development planning in Jiangle County.

Increasing numbers of explanatory variables tend to result in information redundancy and “dimensional disaster” in the quantitative remote sensing of forest aboveground biomass (AGB). Feature selection of model factors is an effective method for improving the accuracy of AGB estimates. Machine learning algorithms are also widely used in AGB estimation, although little research has addressed the use of the categorical boosting algorithm (CatBoost) for AGB estimation. Both feature selection and regression for AGB estimation models are typically performed with the same machine learning algorithm, but there is no evidence to suggest that this is the best method. Therefore, the present study focuses on evaluating the performance of the CatBoost algorithm for AGB estimation and comparing the performance of different combinations of feature selection methods and machine learning algorithms. AGB estimation models of four forest types were developed based on Landsat OLI data using three feature selection methods (recursive feature elimination (RFE), variable selection using random forests (VSURF), and least absolute shrinkage and selection operator (LASSO)) and three machine learning algorithms (random forest regression (RFR), extreme gradient boosting (XGBoost), and categorical boosting (CatBoost)). Feature selection had a significant influence on AGB estimation. RFE preserved the most informative features for AGB estimation and was superior to VSURF and LASSO. In addition, CatBoost improved the accuracy of the AGB estimation models compared with RFR and XGBoost. AGB estimation models using RFE for feature selection and CatBoost as the regression algorithm achieved the highest accuracy, with root mean square errors (RMSEs) of 26.54 Mg/ha for coniferous forest, 24.67 Mg/ha for broad-leaved forest, 22.62 Mg/ha for mixed forests, and 25.77 Mg/ha for all forests. The combination of RFE and CatBoost had better performance than the VSURF–RFR combination in which random forests were used for both feature selection and regression, indicating that feature selection and regression performed by a single machine learning algorithm may not always ensure optimal AGB estimation. It is promising to extending the application of new machine learning algorithms and feature selection methods to improve the accuracy of AGB estimates.

An individual-tree diameter growth model was developed for Cunninghamia lanceolata in Fujian province, southeast China. Data were obtained from 72 plantation-grown China-fir trees in 24 single-species plots. Ordinary non-linear least squares regression was used to choose the best base model from among 5 theoretical growth equations; selection criteria were the smallest absolute mean residual and root mean square error and the largest adjusted coefficient of determination. To account for autocorrelation in the repeated-measures data, we developed one-level and nested two-level nonlinear mixed-effects (NLME) models, constructed on the selected base model; the NLME models incorporated random effects of the tree and plot. The best random-effects combinations for the NLME models were identified by Akaike's information criterion, Bayesian information criterion and −2 logarithm likelihood. Heteroscedasticity was reduced with two residual variance functions, a power function and an exponential function. The autocorrelation was addressed with three residual autocorrelation structures: a first-order autoregressive structure [AR(1)], a combination of first-order autoregressive and moving average structures [ARMA(1,1)] and a compound symmetry structure (CS). The one-level (tree) NLME model performed best. Independent validation data were used to test the performance of the models and to demonstrate the advantage of calibrating the NLME models.

A multiple linear model was developed for individual tree crown width of Cunninghamia lanceolata (Lamb.) Hook in Fujian province, southeast China. Data were obtained from 55 sample plots of pure China-fir plantation stands. An Ordinary Linear Least Squares (OLS) regression was used to establish the crown width model. To adjust for correlations between observations from the same sample plots, we developed one level linear mixed-effects (LME) models based on the multiple linear model, which take into account the random effects of plots. The best random effects combinations for the LME models were determined by the Akaike’s information criterion, the Bayesian information criterion and the -2logarithm likelihood. Heteroscedasticity was reduced by three residual variance functions: the power function, the exponential function and the constant plus power function. The spatial correlation was modeled by three correlation structures: the first-order autoregressive structure [AR(1)], a combination of first-order autoregressive and moving average structures [ARMA(1,1)], and the compound symmetry structure (CS). Then, the LME model was compared to the multiple linear model using the absolute mean residual (AMR), the root mean square error (RMSE), and the adjusted coefficient of determination (adj-R²). For individual tree crown width models, the one level LME model showed the best performance. An independent dataset was used to test the performance of the models and to demonstrate the advantage of calibrating LME models.

The importance of the height-diameter (H-D) relationship in forest productivity is well known. The general nonlinear regression model, based on the mean regression technical, is not able to give a complete description of the H-D relationship. This study aims to evaluate the H-D relationship among relative competition levels and develop a quantile regression (QR) model to fully describe the H-D relationship. The dominance index was applied to determine the relative competition levels of trees for the Chinese fir. Based on the basic Weibull growth model, the mean regression for five relative competition levels and 11 QR models was constructed with 10-fold cross-validation. We have demonstrated that the H-D relationship for the Chinese fir strongly correlated with relative competition states, but the five curves from mean regression models did not show a notable difference between the trends of H-D relationship under different competition levels. Similar regression results were found in QR models of the specific quantiles; the average tree height of five competition levels varied between 5.78% and 17.65% (i.e., about 0.06 and 0.18 quantiles). In addition, some special curves of the H-D relationship such as the QR models of the 0.01 and 0.99 quantile showed the H-D relationship under certain conditions. These findings indicate that the QR models not only evaluated the rates of change of the H-D relationships in various competition levels, but also described their characteristics with more information, like the upper and lower boundary of the conditional distribution of responses. Although the flexible QR curves followed the distribution of the data and showed more information about the H-D relationships, the H-D curves may not intersect with each other, even when the trees reached their maximum height. Hence, the QR model requires further practice in assessing the growth trajectory of the tree’s diameter or tree height to gain better results.