Search

This study presents auxiliary support techniques for tree selection strategies based on the spatial structure indices and three competition indices in secondary forests, and discusses the importance of tree competition in forest management. The spatial structure parameter in the structured management is used as a quantitative index—the uniform angle index and three competition indices are used in the design of the algorithm for selective thinning for secondary forest. Based on the target tree-based management principles, simulation of selective thinning was carried out using GIS and C# programming languages. Data for this study were collected from experimental sample plots at Jilin Wangqing Forestry Bureau in China. The simulation results strongly support the use of auxiliary technology for scientifically selecting trees for thinning, avoiding the subjectivity of the traditional manual selection. Selection is largely based on the uniform angle index and competition index. Hegyi's competition index and its improved version used in the algorithm provided almost identical simulation results, i.e., thinning intensities suggested by these indices for the first sample plot are 21.8% and 21.5%, respectively, and for the second plot are 21.3% and 21.1%, respectively. Thus, one of these competition indices can be used to select trees for thinning. The comprehensive competition index (CCI, a combination of an improved version of Hegyi's competition index with tree species mingling) can avoid the selection of individual trees with high mingling and help maintain the tree species diversity. CCI suggests thinning intensities of 18.3% and 18.4% for the first and second sample plots, respectively. Presented methods and results may provide auxiliary supports for scientific thinning and help promote the application of information technology in forest management. [ABSTRACT FROM AUTHOR]

Northern China has been long threatened by aeolian desertification. In recent years, all levels of the Chinese government have performed a series of ecological protection and sand control projects. To grasp the implementation effects of these projects and adjust policies in time, it is necessary to understand the process of aeolian desertification quickly and accurately. Remote sensing technologies play an irreplaceable role in aeolian desertification monitoring. In this study, the Zhenglan Banner, which is in the hinterland of the Hunshandake Sandy Land, was considered as the research area. Based on unmanned aerial vehicle (UAV) images, ground survey data, and Landsat images called in Google Earth Engine (GEE), the aeolian desertified land (ADL) in 2000, 2004, 2010, 2015, and 2019 was extracted using spectral mixture analysis. A desertification index (DI) was constructed to evaluate the spatial and temporal dynamics of the ADL in the Zhenglan Banner. Finally, a residual analysis explored the driving forces of aeolian desertification. The results showed that (1) the ADL area in the Zhenglan Banner has been trending downwards over the past 20 years but rebounded from 2004 to 2010; (2) over the past 20 years, the area of slightly, moderately, and severely desertified land has decreased at annual rates of 0.4%, 2.7%, and 3.4%, respectively; (3) human activities had significantly positive and negative impacts on the aeolian desertification trend for 20.0% and 21.0% of the study area, respectively, but not for the rest. This paper explored new techniques for rapid aeolian desertification monitoring and is of great significance for controlling and managing aeolian desertification in this region. [ABSTRACT FROM AUTHOR]

Forest growing stem volume (GSV) reflects the richness of forest resources as well as the quality of forest ecosystems. Remote sensing technology enables robust and efficient GSV estimation as it greatly reduces the survey time and cost while facilitating periodic monitoring. Given its red edge bands and a short revisit time period, Sentinel-2 images were selected for the GSV estimation in Wangyedian forest farm, Inner Mongolia, China. The variable combination was shown to significantly affect the accuracy of the estimation model. After extracting spectral variables, texture features, and topographic factors, a stepwise random forest (SRF) method was proposed to select variable combinations and establish random forest regressions (RFR) for GSV estimation. The linear stepwise regression (LSR), Boruta, Variable Selection Using Random Forests (VSURF), and random forest (RF) methods were then used as references for comparison with the proposed SRF for selection of predictors and GSV estimation. Combined with the observed GSV data and the Sentinel-2 images, the distributions of GSV were generated by the RFR models with the variable combinations determined by the LSR, RF, Boruta, VSURF, and SRF. The results show that the texture features of Sentinel-2's red edge bands can significantly improve the accuracy of GSV estimation. The SRF method can effectively select the optimal variable combination, and the SRF-based model results in the highest estimation accuracy with the decreases of relative root mean square error by 16.4%, 14.4%, 16.3%, and 10.6% compared with those from the LSR-, RF-, Boruta-, and VSURF-based models, respectively. The GSV distribution generated by the SRF-based model matched that of the field observations well. The results of this study are expected to provide a reference for GSV estimation of coniferous plantations. [ABSTRACT FROM AUTHOR]

Forest canopy height is one of the most important spatial characteristics for forest resource inventories and forest ecosystem modeling. Light detection and ranging (LiDAR) can be used to accurately detect canopy surface and terrain information from the backscattering signals of laser pulses, while photogrammetry tends to accurately depict the canopy surface envelope. The spatial differences between the canopy surfaces estimated by LiDAR and photogrammetry have not been investigated in depth. Thus, this study aims to assess LiDAR and photogrammetry point clouds and analyze the spatial differences in canopy heights. The study site is located in the Jigongshan National Nature Reserve of Henan Province, Central China. Six data sets, including one LiDAR data set and five photogrammetry data sets captured from an unmanned aerial vehicle (UAV), were used to estimate the forest canopy heights. Three spatial distribution descriptors, namely, the effective cell ratio (ECR), point cloud homogeneity (PCH) and point cloud redundancy (PCR), were developed to assess the LiDAR and photogrammetry point clouds in the grid. The ordinary neighbor (ON) and constrained neighbor (CN) interpolation algorithms were used to fill void cells in digital surface models (DSMs) and canopy height models (CHMs). The CN algorithm could be used to distinguish small and large holes in the CHMs. The optimal spatial resolution was analyzed according to the ECR changes of DSMs or CHMs resulting from the CN algorithms. Large negative and positive variations were observed between the LiDAR and photogrammetry canopy heights. The stratified mean difference in canopy heights increased gradually from negative to positive when the canopy heights were greater than 3 m, which means that photogrammetry tends to overestimate low canopy heights and underestimate high canopy heights. The CN interpolation algorithm achieved smaller relative root mean square errors than the ON interpolation algorithm. This article provides an operational method for the spatial assessment of point clouds and suggests that the variations between LiDAR and photogrammetry CHMs should be considered when modeling forest parameters. [ABSTRACT FROM AUTHOR]

The forest growth and yield models, which are used as important decision-support tools in forest management, are commonly based on the individual tree characteristics, such as diameter at breast height (DBH), crown ratio, and height to crown base (HCB). Taking direct measurements for DBH and HCB through the ground-based methods is cumbersome and costly. The indirect method of getting such information is possible from remote sensing databases, which can be used to build DBH and HCB prediction models. The DBH and HCB of the same trees are significantly correlated, and so their inherent correlations need to be appropriately accounted for in the DBH and HCB models. However, all the existing DBH and HCB models, including models based on light detection and ranging (LiDAR) have ignored such correlations and thus failed to account for the compatibility of DBH and HCB estimates, in addition to disregarding measurement errors. To address these problems, we developed a compatible simultaneous equation system of DBH and HCB error-in-variable (EIV) models using LiDAR-derived data and ground-measurements for 510 Picea crassifolia Kom trees in northwest China. Four versatile algorithms, such as nonlinear seemingly unrelated regression (NSUR), two-stage least square (2SLS) regression, three-stage least square (3SLS) regression, and full information maximum likelihood (FIML) were evaluated for their estimating efficiencies and precisions for a simultaneous equation system of DBH and HCB EIV models. In addition, two other model structures, namely, nonlinear least squares with HCB estimation not based on the DBH (NLS and NBD) and nonlinear least squares with HCB estimation based on the DBH (NLS and BD) were also developed, and their fitting precisions with a simultaneous equation system compared. The leave-one-out cross-validation method was applied to evaluate all estimating algorithms and their resulting models. We found that only the simultaneous equation system could illustrate the effect of errors associated with the regressors on the response variables (DBH and HCB) and guaranteed the compatibility between the DBH and HCB models at an individual level. In addition, such an established system also effectively accounted for the inherent correlations between DBH with HCB. However, both the NLS and BD model and the NLS and NBD model did not show these properties. The precision of a simultaneous equation system developed using NSUR appeared the best among all the evaluated algorithms. Our equation system does not require the stand-level information as input, but it does require the information of tree height, crown width, and crown projection area, all of which can be readily derived from LiDAR imagery using the delineation algorithms and ground-based DBH measurements. Our results indicate that NSUR is a more reliable and quicker algorithm for developing DBH and HCB models using large scale LiDAR-based datasets. The novelty of this study is that the compatibility problem of the DBH model and the HCB EIV model was properly addressed, and the potential algorithms were compared to choose the most suitable one (NSUR). The presented method and algorithm will be useful for establishing similar compatible equation systems of tree DBH and HCB EIV models for other tree species. [ABSTRACT FROM AUTHOR]

Increasing the area of planted forests is rather important for compensation the loss of natural forests and slowing down the global warming. Forest growing stem volume (GSV) is a key indicator for monitoring and evaluating the quality of planted forest. To improve the accuracy of planted forest GSV located in south China, four L-band ALOS PALSAR-2 quad-polarimetric synthetic aperture radar (SAR) images were acquired from June to September with short intervals. Polarimetric characteristics (un-fused and fused) derived by the Yamaguchi decomposition from time series SAR images with different intervals were considered as independent variables for the GSV estimation. Then, the general linear model (GLM) obeyed the exponential distribution were proposed to retrieve the stand-level GSV in plantation. The results show that the un-fused power of double bounce scatters and four fused variables derived from single SAR image is highly sensitive to the GSV, and these polarimeric characteristics derived from the time series images more significantly contribute to improved estimation of GSV. Moreover, compared with the estimated GSV using the semi-exponential model, the employed GLM model with less limitations and simple algorithm has a higher saturation level (nearly to 300 m3/ha) and higher sensitivity to high forest GSV values than the semi-exponential model. Furthermore, by reducing the external disturbance with the help of time average, the accuracy of estimated GSV is improved using fused polarimeric characteristics, and the estimation accuracy of forest GSV was improved as the images increase. Using the fused polarimetric characteristics (Dbl×Vol/Odd) and the GLM, the minimum RRMSE was reduced from 33.87% from single SAR image to 24.42% from the time series SAR images. It is implied that the GLM is more suitable for polarimetric characteristics derived from the time series SAR images and has more potential to improve the planted forest GSV. [ABSTRACT FROM AUTHOR]

Rapidly advancing airborne laser scanning technology has become greatly useful to estimate tree- and stand-level variables at a large scale using high spatial resolution data. Compared with that of ground measurements, the accuracy of the inferred information of diameter at breast height (DBH) from a remotely sensed database and the models developed with traditional regression approaches (e.g., ordinary least square regression) may not be sufficient. Thus, this regression approach is no longer appropriate to develop accurate models and predict DBH from remotely sensed-related variables because DBH is subject to the random effects of forest stands. This study developed a generalized nonlinear mixed-effects DBH estimation model from remotely sensed imagery data. The light detection and ranging (LiDAR)-derived stand canopy density, crown projection area, and tree height were used as predictors in the DBH estimation model. These variables can be more readily measured over an extensive forest area with higher accuracy compared to the conventional field-based methods. The airborne LiDAR data for a total of 402 Picea crassifolia Kom trees on a sample plot that were divided into 16 sub-sample plots and located in the most important distribution region of western China were used. The leave-one sub-sample plot-out cross-validation method was applied to evaluate the model's prediction accuracy. The results indicated that the random effects of the sub-sample plot on the prediction of DBH were large and their inclusion into the DBH model significantly improved the prediction accuracy. The prediction accuracy of the proposed model at the mean (M) response was also substantially improved relative to the accuracy obtained from the base model. Among several tree selection alternatives evaluated, a sample size of the two largest trees per sub-sample plot used in estimating the random effects showed a significantly higher accuracy compared to other sampling alternatives. This sample size would balance both the measurement cost and potential prediction errors. The nonlinear mixed-effects DBH estimation model at the M response can also be applied if obtaining the estimates of individual tree DBH with a relatively lower cost, and a lower prediction accuracy was the purpose of the study. [ABSTRACT FROM AUTHOR]

China has been facing serious land degradation and desertification in its north and northwest arid and semi-arid areas. Monitoring the dynamics of percentage vegetation cover (PVC) using remote sensing imagery in these areas has become critical. However, because these areas are large, remote, and sparsely populated, and also because of the existence of mixed pixels, there have been no accurate and cost-effective methods available for this purpose. Spectral unmixing methods are a good alternative as they do not need field data and are low cost. However, traditional linear spectral unmixing (LSU) methods lack the ability to capture the characteristics of spectral reflectance and scattering from endmembers and their interactions within mixed pixels. Moreover, existing nonlinear spectral unmixing methods, such as random forest (RF) and radial basis function neural network (RBFNN), are often costly because they require field measurements of PVC from a large number of training samples. In this study, a cost-effective approach to mapping PVC in arid and semi-arid areas was proposed. A method for selection and purification of endmembers mainly based on Landsat imagery was first presented. A probability-based spectral unmixing analysis (PBSUA) and a probability-based optimized k nearest-neighbors (PBOkNN) approach were then developed to improve the mapping of PVC in Duolun County in Inner Mongolia, China, using Landsat 8 images and field data from 920 sample plots. The proposed PBSUA and PBOkNN methods were further validated in terms of accuracy and cost-effectiveness by comparison with two LSU methods, with and without purification of endmembers, and two nonlinear approaches, RF and RBFNN. The cost-effectiveness was defined as the reciprocal of cost timing relative root mean square error (RRMSE). The results showed that (1) Probability-based spectral unmixing analysis (PBSUA) was most cost-effective and increased the cost-effectiveness by 29.3% 29.3%, 33.5%, 50.8%, and 53.0% compared with two LSU methods, PBOkNN, RF, and RBFNN, respectively; (2) PBSUA, RF, and RBFNN gave RRMSE values of 22.9%, 21.8%, and 22.8%, respectively, which were not significantly different from each other at the significance level of 0.05. Compatibly, PBOkNN and LSU methods with and without purification of endmembers resulted in significantly greater RRMSE values of 27.5%, 32.4%, and 43.3%, respectively; (3) the average estimates of the sample plots and predicted maps from PBSUA, PBOkNN, RF, and RBFNN fell in the confidence interval of the test plot data, but those from two LSU methods did not, although the LSU with purification of endmembers improved the PVC estimation accuracy by 25.2% compared with the LSU without purification of endmembers. Thus, this study indicated that the proposed PBSUA had great potential for cost-effectively mapping PVC in arid and semi-arid areas. [ABSTRACT FROM AUTHOR]

Uncertainties in forest aboveground biomass (AGB) estimates resulting from over- and underestimations using remote sensing data have been widely studied. The uncertainties may occur due to the spatial effects of the plot data. In this study, we collected AGB data from a total of 147 Pinus densata forest sample plots in Yunnan of southwestern China and analyzed the spatial effects on the estimation of AGB. An ordinary least squares (OLS) and four spatial regression methods were compared for the estimation using Landsat 8-OLI images. Through the spatial analysis of AGB and residuals of model predictions, it was found that the spatial autocorrelation and heterogeneity of the plot data could not be ignored. Compared with the OLS, the impact of the spatial effects on AGB estimation could be reduced slightly by the spatial lag model (SLM) and the spatial error model (SEM) and greatly reduced by the linear mixed effects model (LMM) and geographically weighted regression (GWR) based on the distributions of prediction residuals, global Moran's I, and Z score. The spatial regression models had better performance for model fitting and prediction because of the reduction in overestimations and underestimations for the forests with small and large AGB values, respectively. However, the reductions in the overestimations and underestimations varied depending on the spatial regression models. The GWR provided the most accurate predictions with the largest R2 (0.665), the smallest root mean square error (34.507), and mean relative error (−9.070%) by greatly reducing the AGB interval for overestimations occurring and significantly increasing the threshold of AGB from 150 Mg/ha to 200 Mg/ha for underestimations. Thus, GWR offered the greatest potential of improving the estimation of Pinus densata forest AGB in Yunnan of southwestern China. [ABSTRACT FROM AUTHOR]

Rapid and efficient assessment of cultivated land quality (CLQ) using remote sensing technology is of great significance for protecting cultivated land. However, it is difficult to obtain accurate CLQ estimates using the current satellite-driven approaches in the pressure-state-response (PSR) framework, owing to the limitations of linear models and CLQ spectral indices. In order to improve the estimation accuracy of CLQ, this study used four evaluation models (the traditional linear model; partial least squares regression, PLSR; back propagation neural network, BPNN; and BPNN with genetic algorithm optimization, GA-BPNN) to evaluate CLQ for determining the accurate evaluation model. In addition, the optimal satellite-derived indicator in the land state index was selected among five vegetation indices (the normalized vegetation index, NDVI; enhanced vegetation index, EVI; modified soil-adjusted vegetation index, MSAVI; perpendicular vegetation index, PVI; and soil-adjusted vegetation index, SAVI) to improve the prediction accuracy of CLQ. This study was conducted in Conghua District of Guangzhou, Guangdong Province, China, based on Gaofen-1 (GF-1) data. The prediction accuracies from the traditional linear model, PLSR, BPNN, and GA-BPNN were compared using observations. The results demonstrated that (1) compared with other models (the traditional linear model: R2 = 0.14 and RMSE = 91.53; PLSR: R2 = 0.33 and RMSE = 74.58; BPNN: R2 = 0.50 and RMSE = 61.75), the GA-BPNN model based on EVI in the land state index provided the most accurate estimates of CLQ, with the R2 of 0.59 and root mean square error (RMSE) of 56.87, indicating a nonlinear relationship between CLQ and the prediction indicator; and (2) the GA-BPNN-based evaluation approach of CLQ in the PSR framework was driven to map CLQ of the study area using the GF-1 data, leading to an RMSE of 61.44 at the regional scale, implying that the GA-BPNN-based evaluation approach has the potential to map CLQ over large areas. This study provides an important reference for the high-accuracy prediction of CLQ based on remote sensing technology. [ABSTRACT FROM AUTHOR]

Spectral reflectance distortions caused by terrain and solar illumination seriously reduce the accuracy of mapping forest carbon density, especially in mountainous regions. Many models have been developed for mitigating or eliminating the terrain effects on the quality of remote sensing images in hilly and mountainous areas. However, these models usually use global parameters, which may lead to overcorrections for regions with poor illumination and steep slopes. In this study, we present a local parameter estimation (LPE) method based on a pixel-moving window for topographic correction (TC), which can be considered as a general optimization framework for most semiempirical TC models. We set seven kernel sizes for the presented framework, which are 15 pixels, 25 pixels, 50 pixels, 100 pixels, 250 pixels, 500 pixels, and 1000 pixels, respectively. The proposed method was then applied to four traditional TC models, Minnaert (MIN), C Correction (CC), Sun Canopy Sensor + C (SCSC) and Statistical Empirical Correction (SEC), to form four new TC models. These new models were used to estimate forest carbon density of a mountainous area in Southern China using field plot data and a Landsat 8 image. Four evaluation methods, including correlation analysis, the stability of land covers, comparison of reflectance between sunlit and shaded slopes, and accuracy assessment of forest carbon density, were employed to evaluate the contributions of moving window sizes, and assess the performance of the TC models for forest carbon density estimation. The results show that the four TC models with LPE perform much better than the traditional TC models in reducing the topographic effects and improving the estimation accuracy of forest carbon density for the study area. Among the traditional TC models, SEC performs slightly better than SCSC, CC, and MIN. Therefore, the SEC-based model with LPE, that is, LPE-SEC, gets greater R2 and smaller relative RMSE values in estimating forest carbon density than other models. Moreover, all the means of the predicted forest carbon density values fall in the confidence interval of the validation data at a significant level of 0.05. Overall, this study implies that the proposed method with LPE provides great potential to improve the performance of TC and forest carbon density estimation for the study area. It is expected that the improved TC method can be applied to other mountainous areas to improve the quality of remotely sensed images. [ABSTRACT FROM AUTHOR]

This study proposes a method for determining the optimal image date to improve the evaluation of cultivated land quality (CLQ). Five vegetation indices: leaf area index (LAI), difference vegetation index (DVI), enhanced vegetation index (EVI), normalized difference vegetation index (NDVI), and ratio vegetation index (RVI) are first retrieved using the PROSAIL model and Gaofen-1 (GF-1) images. The indices are then introduced into four regression models at different growth stages for assessing CLQ. The optimal image date of CLQ evaluation is finally determined according to the root mean square error (RMSE). This method is tested and validated in a rice growth area of Southern China based on 115 sample plots and five GF-1 images acquired at the tillering, jointing, booting, heading to flowering, and milk ripe and maturity stage of rice in 2015, respectively. The results show that the RMSEs between the measured and estimated CLQ from four vegetation index-based regression models at the heading to flowering stage are smaller than those at the other growth stages, indicating that the image date corresponding with the heading to flowering stage is optimal for CLQ evaluation. Compared with other vegetation index-based models, the LAI-based logarithm model provides the most accurate estimates of CLQ. The optimal model is also driven using the GF-1 image at the heading to flowering stage to map CLQ of the study area, leading to a relative RMSE of 14.09% at the regional scale. This further implies that the heading to flowering stage is the optimal image time for evaluating CLQ. This study is the first effort to provide an applicable method of selecting the optimal image date to improve the estimation of CLQ and thus advanced the literature in this field. [ABSTRACT FROM AUTHOR]

For the planning and sustainable management of forest resources, well-managed plantations are of great significance to mitigate the decrease of forested areas. Monitoring these planted forests is essential for forest resource inventories. In this study, two ALOS PALSAR-2 quad-polarimetric synthetic aperture radar (SAR) images and ground measurements were employed to estimate growing stem volume (GSV) of Chinese fir plantations located in a hilly area of southern China. To investigate the relationships between forest GSV and polarization characteristics, single and fused variables were derived by the Yamaguchi decomposition and the saturation value of GSV was estimated using a semi-exponential empirical model as a base model. Based on the estimated saturation values and relative root mean square error (RRMSE), the single and fused characteristics and corresponding models were selected and integrated, which led to a novel saturation-based multivariate method used to improve the GSV estimation and mapping of Chinese fir plantations. The new findings included: (1) All the original polarimetric characteristics, statistically, were not significantly correlated with the forest GSV, and their logarithm and ratio transformation fused variables greatly improved the correlations, thus the estimation accuracy of the forest GSV. (2) The logarithm transformation of surface scattering resulted in the greatest saturation, value but the logarithm transformation of double-bounce scattering resulted in the smallest RRMSE of the GSV estimates. (3) Compared with the single transformations, the fused variables led to more reasonable saturation values and obviously reduced the values of RRMSE. (4) The saturation-based multivariate method led to more accurate estimates of the forest GSV than the univariate method, with the smallest value (29.64%) of RRMSE achieved using the set of six variables. This implied that the novel saturation-based multivariate method provided greater potential to improve the estimation and mapping of the forest GSV. [ABSTRACT FROM AUTHOR]

Quickly and efficiently monitoring soil heavy metal content is crucial for protecting the natural environment and for human health. Estimating heavy metal content in soils using hyperspectral data is a cost-efficient method but challenging due to the effects of complex landscapes and soil properties. One of the challenges is how to make a lab-derived model based on soil samples applicable to mapping the contents of heavy metals in soil using air-borne or space-borne hyperspectral imagery at a regional scale. For this purpose, our study proposed a novel method using hyperspectral data from soil samples and the HuanJing-1A (HJ-1A) HyperSpectral Imager (HSI). In this method, estimation models were first developed using optimal relevant spectral variables from dry soil spectral reflectance (DSSR) data and field observations of soil heavy metal content. The relationship of the ratio of DSSR to moisture soil spectral reflectance (MSSR) with soil moisture content was then derived, which built up the linkage of DSSR with MSSR and provided the potential of applying the models developed in the laboratory to map soil heavy metal content at a regional scale using hyperspectral imagery. The optimal relevant spectral variables were obtained by combining the Boruta algorithm with a stepwise regression and variance inflation factor. This method was developed, validated, and applied to estimate the content of heavy metals in soil (As, Cd, and Hg) in Guangdong, China, and the Conghua district of Guangzhou city. The results showed that based on the validation datasets, the content of Cd could be reliably estimated and mapped by the proposed method, with relative root mean square error (RMSE) values of 17.41% for the point measurements of soil samples from Guangdong province and 17.10% for the Conghua district at the regional scale, while the content of heavy metals As and Hg in soil were relatively difficult to predict with the relative RMSE values of 32.27% and 28.72% at the soil sample level and 51.55% and 36.34% at the regional scale. Moreover, the relationship of the DSSR/MSSR ratio with soil moisture content varied greatly before the wavelength of 1029 nm and became stable after that, which linked DSSR with MSSR and provided the possibility of applying the DSSR-based models to map the soil heavy metal content at the regional scale using the HJ-1A images. In addition, it was found that overall there were only a few soil samples with the content of heavy metals exceeding the health standards in Guangdong province, while in Conghua the seriously polluted areas were mainly distributed in the cities and croplands. This study implies that the new approach provides the potential to map the content of heavy metals in soil, but the estimation model of Cd was more accurate than those of As and Hg. [ABSTRACT FROM AUTHOR]