In many countries cadmium (Cd) and arsenic (As) commonly coexist in soils contaminated by mining activities, and can easily enter the human body via consumption of leafy vegetables, like the popularly consumed pakchoi (Brassica chinensis L.), causing major health concerns. In the present study, bioaccessibility and human exposure of Cd and As were assessed in twenty genotypes of pakchoi cultured at two different levels of co-contamination to identify low health risk genotypes. The bioaccessibilities of Cd and As represent a fraction of the total metals content could be bioaccessible for human, in the present study, significant differences in pakchoi Cd and As bioaccessibility were observed among all tested genotypes and co-contaminated levels. Cd and As bioaccessibility of pakchoi were in the ranges of 24.0–87.6% and 20.1–82.5%, respectively, for in the high level co-contaminated soils, which was significantly higher than for low level co-contaminated soils with 7.9–71.8% for Cd bioaccessibility and 16.1–59.0% for As bioaccessibility. The values of bioaccessible established daily intakes (BEDI) and the total bioaccessible target hazard quotients (TBTHQ) of Cd and As were also considerably higher in high level co-contaminated soils than in low level co-contaminated soils. Two genotypes (Meiguanqinggengcai and Zhenqing60F1) contained relatively low concentrations and bioaccessible Cd and As and, their BEDI and TBTHQ for Cd and As ranged below the tolerable limits set by the FAO/WHO (BEDI of Cd < 0.83 μg kg⁻¹ bw day⁻¹, BEDI of As < 3 μg kg⁻¹ bw day⁻¹) and United States Environmental Protection Agency (TBTHQ for Cd and As < 1), this applied for both levels of co-contaminated soils for adults and children. Consequently, these findings suggest identification of safe genotypes in leafy vegetable with low health risk via genotypic screening and breeding methods could be a useful strategy to ensure the safety of food crops grown in those Cd and As co-contaminated fields due to mining activities.

Estimating the aboveground biomass (AGB) at the plot level plays a major role in connecting accurate single-tree AGB measurements to relatively difficult regional AGB estimates. However, AGB estimates at the plot level suffer from many uncertainties. The goal of this study is to determine whether combining machine learning with spatial statistics reduces the uncertainty of plot-level AGB estimates. To illustrate this issue, this study evaluates and compares the performance of different models for estimating plot-level forest AGB. These models include three different machine learning models [support vector machine (SVM), random forest (RF), and a radial basis function artificial neural network (RBF-ANN)], one spatial statistic model (P-BSHADE), and three combinations thereof (SVM & P-BSHADE, RF & P-BSHADE, and RBF-ANN & P-BSHADE). The results show that the root mean square error, mean absolute error, and mean relative error of all combined models are substantially smaller than those of any individual model, with the RF & P-BSHADE combined method generating the smallest values. These results indicate that a combined approach using machine learning with spatial statistics, especially the RF & P-BSHADE model, improves the accuracy of plot-level AGB models. These research results contribute to the development of accurate large-forested-landscape AGB maps.

Sensitivity analysis and parameter optimization of stand models can improve their efficiency and accuracy, and increase their applicability. In this study, the sensitivity analysis, screening, and optimization of 63 model parameters of the Physiological Principles in Predicting Growth (3PG) model were performed by combining a sensitivity analysis method and the Markov chain Monte Carlo (MCMC) method of Bayesian posterior estimation theory. Additionally, a nine-year observational dataset of Chinese fir trees felled in the Shunchang Forest Farm, Nanping, was used to analyze, screen, and optimize the 63 model parameters of the 3PG model. The results showed the following: (1) The parameters that are most sensitive to stand stocking and diameter at breast height (DBH) are nWs(power in stem mass vs. diameter relationship), aWs(constant in stem mass vs. diameter relationship), alphaCx(maximum canopy quantum efficiency), k(extinction coefficient for PAR absorption by canopy), pRx(maximum fraction of NPP to roots), pRn(minimum fraction of NPP to roots), and CoeffCond(defines stomatal response to VPD); (2) MCMC can be used to optimize the parameters of the 3PG model, in which the posterior probability distributions of nWs, aWs, alphaCx, pRx, pRn, and CoeffCond conform to approximately normal or skewed distributions, and the peak value is prominent; and (3) compared with the accuracy before sensitivity analysis and a Bayesian method, the biomass simulation accuracy of the stand model was increased by 13.92%, and all indicators show that the accuracy of the improved model is superior. This method can be used to calibrate the parameters and analyze the uncertainty of multi-parameter complex stand growth models, which are important for the improvement of parameter estimation and simulation accuracy.

Disturbances to forests are getting worse with climate change and urbanization. Assessing the functionality of forest ecosystems is challenging because it requires not only a large amount of input data but also comprehensive estimation indicator methods. The object of the evaluation index of forest ecosystem restoration relies on the ecosystem function instead of the area. To develop the appropriate index with ecological implications, we built the hybrid assessment approach including ecosystem structure-function-habitat representatives. It was based on the Normalized Burn Ratio (NBR) spectral indicator and combined with the local forest management inventory (LFMI), Landsat, Light Detection and Ranging (LiDAR) data. The results of the visual interpretation of Google Earth’s historical imagery showed that the total accuracy of the hybrid approach was 0.94. The output of the hybrid model increased as the biodiversity index value increased. Furthermore, to solve the multi-source data availability problem, the random forest model (R² = 0.78, RMSE = 0.14) with 0.77 total accuracy was built to generate an annual recovery index. A random forest model based on tree age is provided to simplify the hybrid approach while extending the results on time series. The recovery index obtained by the random forest model could facilitate monitoring the forest recovery rate of cold spots. The regional ecological recovery time could be predicted. These two results could provide a scientific basis for forest managers to make more effective forest restoration plans. From the perspective of space, it could ensure that the areas with slow recovery would be allocated enough restoration resources. From the perspective of time, the implementation period of the closed forest policy could also be estimated.

Urbanization transforms landscapes from natural ecosystems to configurations of impervious surfaces and green spaces, leading to urban heat island effects that impact health and ecosystem sustainability. This study in Xiamen City, China, categorizes urban areas into functional zones, employs Random Forest and Stepwise Regression models to assess thermal differences, and proposes optimization measures for the building–green space landscape. The optimization involves altering the characterization of the building–green space landscape pattern. Results indicate: (1) due to the spatial heterogeneity of the building–green space landscape pattern in different functional zones, the surface temperature also shows strong spatial heterogeneity in different functional zones; (2) different optimization measures for the building–green space pattern are needed for different functional zones; taking the urban residential zone as an example, the Normalized Difference Vegetation Index (NDVI) in the hot spot area can be adjusted according to the value range of the cold spot area; (3) considering the solar radiation process, Sun View Factor (S_unVF) plays an important role in indicating the change in surface temperature in the commercial service area, and as S_unVF increases, the surface temperature of the functional zone tends to rise. This research offers insights into urban thermal environment improvement and landscape pattern optimization.