Forest Mapping Research Articles

Comparing harmonic regression and GLAD Phenology metrics for estimation of forest community types and aboveground live biomass within forest inventory and analysis plots

The Landsat multispectral time series is a valuable source of moderate spatial resolution data to support forest mapping and monitoring tasks. Using United States Department of Agriculture (USDA) Forest Service Forest Inventory and Analysis (FIA) plots within the states of Michigan, Oregon, and West Virginia, two methods to summarize time series observations, harmonic regression coefficients and Global Land Analysis & Discovery (GLAD) Phenology metrics, are compared for predicting forest community type, total aboveground live biomass (AGLBM), and species-specific AGLBM. Harmonic regression coefficients, which provided mean overall accuracies (OAs) between 62.8% and 73.1% and map image classification efficacies (MICEs) varying between 0.455 and 0.566 for the three studied states and calculated using multiple model replicates, generally provided better predictive performance for differentiating forest community types in comparison to GLAD phenology metrics. However, differences were not always statistically significant. DTM-derived terrain variables improved classification performance in some landscapes when using machine learning, random forest classification (for example, the highest obtained mean OA of 78.5% (MICE = 0.557) was obtained for West Virginia when using the combined harmonic regression and terrain variables). For the regression-based estimation of total and species-specific AGLBM, the spectral data used and the incorporation of terrain variables had less of an impact on model performance. Further, the regression models in Oregon provided larger mean R-squared values (0.452 to 0.490) in comparison to those in Michigan and West Virginia, where all R-squared values were lower than 0.200, suggesting that AGLBM prediction may be more or less challenging in different landscapes depending on forest characteristics, terrain, management practices, and disturbance histories.

Winter remote sensing images are more suitable for forest mapping in Jiangxi Province

ABSTRACT Jiangxi Province boasts the second-highest forest coverage in China. Its forests play a crucial role in providing essential ecosystem services and maintaining the ecological health of the region. High-resolution and high-precision forest mapping are significant in the timely and accurate monitoring of dynamic forest changes to support sustainable forest management. This study used Sentinel-2 images from four seasons in the Google Earth Engine (GEE) platform to map forest distribution. Moreover, the classification results were compared and analyzed using different classification algorithms and feature-variable combinations. Based on the overall accuracy, the optimal image seasonality, feature combinations and classification algorithms were selected, and the forest maps of Jiangxi Province were mapped from 2019 to 2021. The accuracy evaluation showed that the winter image classification results had the highest accuracy (above 0.88). The red edge bands carried by Sentinel-2 could effectively improve the classification accuracy. The Random Forest classifier is the optimal classification algorithm for forest mapping in Jiangxi Province. The forest mapping obtained can be used for ecological health assessment and ecosystem function. The study provides a scientific basis for accurate and timely extraction of forest cover and can serve as a valuable resource for forest management planning and future research.

Spatial database of planted forests in East Asia

Planted forests are critical to climate change mitigation and constitute a major supplier of timber/non-timber products and other ecosystem services. Globally, approximately 36% of planted forest area is located in East Asia. However, reliable records of the geographic distribution and tree species composition of these planted forests remain very limited. Here, based on extensive in situ and remote sensing data, as well as an ensemble modeling approach, we present the first spatial database of planted forests for East Asia, which consists of maps of the geographic distribution of planted forests and associated dominant tree genera. Of the predicted planted forest areas in East Asia (948,863 km2), China contributed 87%, most of which is located in the lowland tropical/subtropical regions, and Sichuan Basin. With 95% accuracy and an F1 score of 0.77, our spatially-continuous maps of planted forests enable accurate quantification of the role of planted forests in climate change mitigation. Our findings inform effective decision-making in forest conservation, management, and global restoration projects.

The development of a global LAI and FAPAR product using GCOM-C/SGLI data

The Japan Aerospace Exploration Agency (JAXA) launched the Global Change Observation Mission - Climate (GCOM-C) satellite on December 23rd, 2017. As a part of the standard products of JAXA’s GCOM-C satellite, we developed the Leaf Area Index (LAI) and the Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) product. In comparison to other global LAI and FAPAR products by the National Aeronautics and Space Administration (NASA) and the Copernicus Global Land Service (CGLS), our products have the following advantageous features: high spatial and temporal resolution (250 m and daily); the retrieval algorithm (look-up-table) is based on “the Forest Light Environmental Simulator” (FLiES) which is a full 3-D canopy radiative transfer model as well as a vegetation landscape model and a global forest landscape map created for this study. The algorithm and ancillary data sets are completely independent of the existing products that are more or less influenced by the MODIS LAI/FAPAR, providing an alternate and independent source of LAI and FAPAR. The dataset has been created since April 2018 and is freely available on JAXA’s G-Portal site (https://gportal.jaxa.jp/gpr). As a result of accuracy assessment with in-situ data, the Root Mean Square Errors (RMSEs) for total LAI and FAPAR were 0.776 and 0.112, respectively. We also compared our LAI and FAPAR estimates with three existing global datasets (PROBA-V/OLCI, MODIS and VIIRS product). The spatial distribution of LAI and FAPAR showed different characteristics among the datasets although their temporal trends were similar. The differences in estimation results among datasets due to differences in the retrieval algorithms, base map (or land cover data) and input data need to be evaluated in the future.

Spatial-temporal mapping of forest vegetation cover changes along highways in Brunei using deep learning techniques and Sentinel-2 images

Infrastructure development is a leading driver of forest cover loss in the tropics, resulting in a significant decrease in biodiversity. With recent advancements in digital image processing and satellite imagery techniques, this study presents an approach for tracking forest cover alterations along a recently constructed highway in Brunei. Specifically, this paper harnesses the power of cutting-edge deep learning (DL) models; U-Net and Deeplabv3+, to classify forest and non-forest features from Sentinel-2 satellite images captured along the Telisai–Lumut highway between 2018 and 2022. A computational performance analysis and network comparison revealed that the Deeplabv3+ had a remarkable semantic segmentation performance, with a mean intersection over union of 0.9247, followed by an ensemble of four networks, with a mean intersection over union of 0.9171. The best model consistently outperformed other models in segmenting forest and non-forest features, with an intersection over union (IOU) of 0.9463 and 0.9031, respectively. Moreover, deploying hyperparameter optimization allowed us to assess the model's sensitivity to different hyperparameters for better results. Furthermore, forest change statistics were computed using the best model's generated forest change maps, and the results indicated a 16.7% gain in the overall forest vegetation cover along the highway, demonstrating the validity of the results. These reproducible findings will provide forest managers and ecologists with near real-time information on the status, health, and extent of forest ecosystems. This study's results have critical implications for using DL to monitor forest cover changes and will be of great interest to the scientific community as they open new avenues for future research and underline the need for continued innovation in this important field.

Sentinel-1 Time Series for Predicting Growing Stock Volume of Boreal Forest: Multitemporal Analysis and Feature Selection

Copernicus Sentinel-1 images are widely used for forest mapping and predicting forest growing stock volume (GSV) due to their accessibility. However, certain important aspects related to the use of Sentinel-1 time series have not been thoroughly explored in the literature. These include the impact of image time series length on prediction accuracy, the optimal feature selection approaches, and the best prediction methods. In this study, we conduct an in-depth exploration of the potential of long time series of Sentinel-1 SAR data to predict forest GSV and evaluate the temporal dynamics of the predictions using extensive reference data. Our boreal coniferous forests study site is located near the Hyytiälä forest station in central Finland and covers an area of 2500 km2 with nearly 17,000 stands. We considered several prediction approaches and fine-tuned them to predict GSV in various evaluation scenarios. Our analyses used 96 Sentinel-1 images acquired over three years. Different approaches for aggregating SAR images and choosing feature (predictor) variables were evaluated. Our results demonstrate a considerable decrease in the root mean squared errors (RMSEs) of GSV predictions as the number of images increases. While prediction accuracy using individual Sentinel-1 images varied from 85 to 91 m3/ha RMSE, prediction accuracy with combined images decreased to 75.6 m3/ha. Feature extraction and dimension reduction techniques facilitated the achievement of near-optimal prediction accuracy using only 8–10 images. Examined methods included radiometric contrast, mutual information, improved k-Nearest Neighbors, random forests selection, Lasso, and Wrapper approaches. Lasso was the most optimal, with RMSE reaching 77.1 m3/ha. Finally, we found that using assemblages of eight consecutive images resulted in the greatest accuracy in predicting GSV when initial acquisitions started between September and January.

Forest Height Inversion via RVoG Model and Its Uncertainties Analysis via Bayesian Framework—Comparisons of Different Wavelengths and Baselines

Accurate estimation of forest height over a large area is beneficial to reduce the uncertainty of forest carbon sink estimation, which is of great significance to the terrestrial carbon cycle, global climate change, forest resource management, and forest-related scientific research. Forest height inversion using polarimetric interferometry synthetic aperture radar (PolInSAR) data through Random volume over ground (RVoG) models has demonstrated great potential for large-area forest height mapping. However, the wavelength and baseline length used for the PolInSAR data acquisition plays an important role during the forest height estimation procedure. In this paper, X–, C–, L–, and P–band PolInSAR datasets with four different baseline lengths were simulated and applied to explore the effects of wavelength and baseline length on forest height inversion using RVoG models. Hierarchical Bayesian models developed with a likelihood function of RVoG model were developed for estimated results uncertainty quantification and decrease. Then a similar procedure was applied in the L– and P–band airborne PolInSAR datasets with three different baselines for each band. The results showed that (1) Wavelength showed obvious effects on forest height inversion results with the RVoG model. For the simulated PolInSAR datasets, the L– and P–bands performed better than the X– and C–bands. The best performance was obtained at the P–band with a baseline combination of 10 × 4 m with an absolute error of 0.05 m and an accuracy of 97%. For the airborne PolInSAR datasets, an L–band with the longest baseline of 24 m in this study showed the best performance with R2 = 0.64, RMSE = 3.32 m, and Acc. = 77.78%. (2) It is crucial to select suitable baseline lengths to obtain accurate forest height estimation results. In the four baseline combinations of simulated PolInSAR datasets, the baseline combination of 10 × 4 m both at the L– and P–bands performed best than other baseline combinations. While for the airborne PolInSAR datasets, the longest baseline in three different baselines obtained the highest accuracy at both L– and P–bands. (3) Bayesian framework is useful for estimation results uncertainty quantification and decrease. The uncertainties related to wavelength and baseline length. The uncertainties were reduced obviously at longer wavelengths and suitable baselines.

A novel HHO-RSCDT ensemble learning approach for forest fire danger mapping using GIS

Accurate prediction models for spatial prediction of forest fire danger play a vital role in predicting forest fires, which can help prevent and mitigate the detrimental effects of such disasters. This research aims to develop a new ensemble learning model, HHO-RSCDT, capable of accurately predicting spatial patterns of forest fire danger. The HHO-RSCDT method combines three distinct components, namely Random Subspace (RS), Credal Decision Tree (CDT), and Harris Hawks Optimizer (HHO). Herein, RS generates a series of subspace datasets, which are subsequently utilized to produce individual CDT classifiers. Then, HHO optimizes the ensemble model, enabling the model to achieve higher predictive performance. The model was trained and validated using a Phu Yen province, Vietnam dataset. The dataset includes 306 forest fire locations and ten influencing factors from the study province. The results showed the capability of the HHO-RSCDT model in predicting forest fire danger, with an accuracy rate of 83.7%, a kappa statistic of 0.674, and an AUC of 0.911. A comparison between the HHO-RSCDT model and two state-of-the-art machine learning methods, i.e., support vector machine (SVM) and random forest (RF), indicated that the HHO-RSCDT model could perform better, making it a valuable tool for modeling forest fire danger. The forest fire danger map produced using this novel model could be a new tool for local authorities in the Phu Yen province, assisting them in managing and protecting the forest ecosystem. By providing a detailed overview of the are as most susceptible to forest fires, the map can help authorities to develop targeted and effective forest management strategies, such as focusing on areas with high fuel loads or implementing controlled burning programs.

A New Approach Based on TensorFlow Deep Neural Networks with ADAM Optimizer and GIS for Spatial Prediction of Forest Fire Danger in Tropical Areas

Frequent forest fires are causing severe harm to the natural environment, such as decreasing air quality and threatening different species; therefore, developing accurate prediction models for forest fire danger is vital to mitigate these impacts. This research proposes and evaluates a new modeling approach based on TensorFlow deep neural networks (TFDeepNN) and geographic information systems (GIS) for forest fire danger modeling. Herein, TFDeepNN was used to create a forest fire danger model, whereas the adaptive moment estimation (ADAM) optimization algorithm was used to optimize the model, and GIS with Python programming was used to process, classify, and code the input and output. The modeling focused on the tropical forests of the Phu Yen Province (Vietnam), which incorporates 306 historical forest fire locations from 2019 to 2023 and ten forest-fire-driving factors. Random forests (RF), support vector machines (SVM), and logistic regression (LR) were used as a baseline for the model comparison. Different statistical metrics, such as F-score, accuracy, and area under the ROC curve (AUC), were employed to evaluate the models’ predictive performance. According to the results, the TFDeepNN model (with F-score of 0.806, accuracy of 79.3%, and AUC of 0.873) exhibits high predictive performance and surpasses the performance of the three baseline models: RF, SVM, and LR; therefore, TFDeepNN represents a novel tool for spatially predicting forest fire danger. The forest fire danger map from this study can be helpful for policymakers and authorities in Phu Yen Province, aiding sustainable land-use planning and management.

High-resolution forest age mapping based on forest height maps derived from GEDI and ICESat-2 space-borne lidar data

Forest age is a key parameter for estimating forest growth and carbon uptake and for forest management. Remote sensing provides indirect but useful information for mapping forest age at large scales. However, existing regional and global forest age products were generated at low spatial resolutions (often 1000 m) and are not useful for most forest stands in China that are smaller than 1000 m. This study aims to map forest age at the 30 m resolution based on forest height maps mainly derived from the Global Ecosystem Dynamics Investigation (GEDI) and Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) data, and analyze the roles of auxiliary data including temperature, precipitation, slope, and aspect in forest age mapping. Forest age is defined as the average age of dominant tree species within a pixel. Five commonly-used stand growth equations and twelve machine learning methods were tested for their suitability for mapping forest age of different tree species. We found that the Logistic equation performed the best among the tested stand growth equations and the Random Forest (RF) was the best among the tested machine learning methods. According to RF, forest height contributed predominantly to the variance in forest age mapping, while temperature, precipitation, slope, and aspect also had an overall non-negligible and variable contribution among different tree species. By integrating the climate and topographical variables, RF was applicable for forest age mapping without classifying the tree species. These results show that forest height maps derived from space-borne lidar data such as GEDI and ICESat-2 data are highly useful for mapping forest stand age, and the methodology developed in this study highlights a perspective for generating national and global forest age products at a high spatial resolution.

Potential priority areas for forest-dwelling species in Spain based on the degree of forest fragmentation

ABSTRACT Habitat fragmentation is a process that may impair the “habitat provision” ecosystem service provided by forests. Thus, it is relevant to quantify the degree of forest fragmentation, since higher levels are expected to have greater effects on forest species. Our goal was to deploy the Forest Area Density (FAD) metric as a tool to derive maps localizing potential priority areas for species requiring large, relatively undisturbed blocks of forest. The Spanish Forest Map (1:50,000) provided comprehensive data on forestland in the country. We defined forest types considering the ten most abundant species in Spain and pure and mixed stands. Then, we calculated FAD by forest type and produced a set of maps showing the results by general management scenario, namely habitat conservation and restoration. To develop forest planning, specific actions can be implemented within these areas in later steps, fostering biodiversity at national or smaller scales.

Reconstructing Long‐Term Forest Age of China by Combining Forest Inventories, Satellite‐Based Forest Age and Forest Cover Data Sets

AbstractForest age is one of the most important ecosystem characters for accurately estimating the magnitude and potential of carbon sink in forest ecosystems. During the past 40 years, national ecological restoration projects have led to the near doubling of the forest cover area in China, which has also substantially affected the dynamics of forest age. Therefore, there is an urgent need to generate long‐term forest age maps for China. This study reconstructed China forest age datasets (CFAD) from 1980 to 2015 at five year intervals at a 1 km spatial resolution by merging a satellite‐based forest age map in 2010 and forest cover dynamic maps from 1980 to 2015. The random forest method was used to reconstruct the forest age where forest age could not be inferred from the forest age base map in 2010 directly. CFAD showed a good agreement with the province‐level mean forest age derived from the several national forest inventories (R2 ranged from 0.66 to 0.86). In general, the younger forests are mainly distributed in southern and eastern China. The older forests are mainly distributed in the mountain areas of northeast, northwest and southwest China. The average age of China's forests increased from 18.2 to 44.0 years old from 1980 to 2015. Based on the current forest age and future afforestation planning, the average forest age in China is predicted to reach 71.6 years old in 2060. The CFAD provides an alternative data set to obtain improved estimates of local and national forest carbon sinks in China.

WILDFIRE RISKS IN THE SOUTHWEST OF KY SON DISTRICT, NGHE AN PROVINCE – A MULTI-CRITICAL MODEL

Wildfires are considered one of the most common hazards in mountainous areas, posing a serious threat to flora, fauna, and humans. This article focuses on building a forest fire risk map based on GIS techniques in the southwest area of Ky Son district, Nghe An province. The selected factors to create the forest fire risk map include 9 indices: NDDI drought index, surface temperature, slope, elevation, distance to roads, distance to rivers and streams, distance to upland fields, distance to residential areas, and aspect. The results of determining the weights of these factors using the AHP hierarchical analysis technique show that drought and surface temperature are the two factors that have the greatest influence on the forest fire risk index. The limit values that lead to forest fire risks for the drought criterion are 0.2 - 0.4 mm and for the surface temperature criterion are 25 - 30°C.The area with high fire risk accounts for the largest proportion, 38%, covering an area of 24,519.78 ha, followed by the average fire risk with 21,012.48 ha, accounting for 33% of the entire study area. The area with extremely high fire risk occupies the smallest area, 1,283.76 ha (2%). While controlling wildfires is challenging, forest fire risk maps can detect areas with higher risk and mitigate them. Therefore, creating forest fire risk maps for the identified areas is necessary.

Monitoring direct drivers of small-scale tropical forest disturbance in near real-time with Sentinel-1 and -2 data

Advancements in satellite-based forest monitoring increasingly enable the near real-time detection of small-scale tropical forest disturbances. However, there is an urgent need to enhance such monitoring methods with automated direct driver attributions to detected disturbances. This would provide important additional information to make forest disturbance alerts more actionable and useful for uptake by different stakeholders. In this study, we demonstrate spatially explicit and near real-time methods to monitor direct drivers of small-scale tropical forest disturbance across a range of tropical forest conditions in Suriname, the Republic of the Congo and the Democratic Republic of the Congo. We trained a convolutional neural network with Sentinel-1 and Sentinel-2 data to continuously classify newly detected RAdar for Detecting Deforestation (RADD) alerts as smallholder agriculture, road development, selective logging, mining or other. Different monitoring scenarios were evaluated based on varying sensor combinations, post-disturbance time periods and confidence levels. In general, the use of Sentinel-2 data was found to be most accurate for driver classifications, especially with data composited over a period of 4 to 6 months after the disturbance detection. Sentinel-1 data showed to be valuable for more rapid classifications of specific drivers, especially in areas with persistent cloud cover. Throughout all monitoring scenarios, smallholder agriculture was classified most accurately, while road development, selective logging and mining were more challenging to distinguish. An accuracy assessment throughout the full extent of our study regions revealed a Macro-F1 score of 0.861 and an Overall Accuracy of 0.897 for the best performing model, based on the use of 6-month post-disturbance Sentinel-2 composites. Finally, we addressed three specific monitoring use cases that relate to rapid law enforcement against illegal activities, ecological impact assessments and timely carbon emission reporting, by optimizing the trade-off in classification timeliness and confidence to reach required accuracies. Our findings demonstrate the strong capacities of high spatiotemporal resolution satellite data for monitoring direct drivers of small-scale forest disturbance, considering different user interests. The produced forest disturbance driver maps can be accessed via: https://bartslagter94.users.earthengine.app/view/forest-disturbance-drivers.

Advancing wildlife connectivity in land use planning: a case study with four‐toed salamanders

AbstractStable habitat connections that wildlife can safely traverse are essential to biodiversity conservation and healthy ecosystems. We developed high‐resolution landscape connectivity models to predict resistance to movement by a threatened wetland‐obligate amphibian, the four‐toed salamander (Hemidactylium scutatum), and identified priority management areas on the 13,000‐ha Department of Energy Oak Ridge Reservation (ORR) from 2019 to 2022. We developed a resistance surface based on aerial light detection and ranging data (LiDAR), >30 years of field‐based mapping of forest, hydrologic, and geologic features, and contemporary population surveys, alongside derived predictors at <1‐m resolution. We then modeled predicted movement corridors using a circuit theory‐based modeling approach. We worked closely with land management and natural resources personnel to integrate ecological modeling with broader land use priorities, monetary costs, and feasibility. We identified important terrestrial and aquatic areas on ORR and simulated management scenarios to promote stable connections for four‐toed salamanders. This approach allowed us to narrow down a list of 438 potential habitat manipulation sites to 10 sites where open‐bottomed culverts and buffers could be implemented. This smaller‐scale restoration approach produced a similar increase in landscape connectivity while costing <20% of a larger‐scale approach based on barrier removal. We successfully identified feasible, cost‐effective management strategies that integrated knowledge from a variety of sources. We offer a strategy that permitted integration of wildlife management goals into infrastructure upgrades wherein wildlife was not an initial consideration.

A Rapid and Easy Way for National Forest Heights Retrieval in China Using ICESat-2/ATL08 in 2019

Continuous and extensive monitoring of forest height is essential for estimating forest above-ground biomass and predicting the ability of forests to absorb CO2. In particular, forest height at the national scale is an important indicator reflecting the national forestry economic construction, environmental governance, and ecological balance. However, the lack of inventory data restricts large-scale monitoring of forest height to some extent. Conducting manual surveys of forest height for large-scale areas would be labor-intensive and time-consuming. The successful launch of the new generation of spaceborne light detection and ranging (LiDAR) (The Ice, Cloud, and Land Elevation Satellite-2/the Advanced Topographic Laser Altimeter System, ICESat-2/ATLAS) has brought new opportunities for national-scale forestry resource surveys. This paper explores a method to survey national forest canopy height from the new generation of ICESat-2/ATLAS data. In view of the sparse sampling and little overlap between repeated spaceborne LiDAR data, a strategy for assessing the overall change of canopy height for large scales is provided. Some spatially continuous ancillary data were used to assist ICESat-2/ATLAS data to generate a wall-to-wall (spatially continuous) forest canopy height map in China by using the machine learning approach and then quantifying the analysis of forest canopy height in various provinces. The results show that there is a good correlation between the model forest height and the verification data, with a root mean squared error (RMSE) of 3.30 m and a coefficient of determination (R2) of 0.87. This indicates that the method for retrieving national forest canopy height is reliable. There are some limitations in areas with lower vegetation coverage or complex topography which need additional filtering or terrain correction to achieve higher accuracy in measuring forest canopy height. Our analysis suggests that ICESat-2/ATLAS data can achieve the retrieval of national forest height at an overall level, and it would be feasible to use ICESAT-2/ATLAS products to estimate forest canopy height change for large-scale areas.

Source analysis and health risk assessment of heavy metals in agricultural land of multi-mineral mining and smelting area in the Karst region – a case study of Jichangpo Town, Southwest China

In the Karst region of Southwest China, the content of soil heavy metals is generally high because of the geological background. Moreover, Southwest China is rich in mineral resources. A large number of mining and smelting activities discharge heavy metals into surrounding soil and cause superimposed pollution, which has drawn widespread concern. Due to the large variation coefficients of soil heavy metals in the Karst region, it is particularly essential to select appropriate analysis methods. In this paper, Jichangpo in Puding County, a Karst area with multi-mineral mining and smelting, is selected as the research object. A total of 368 pieces of agricultural topsoil in the study area are collected. The pollution level of heavy metals in agricultural soil is evaluated by the geological accumulation index (Igeo) and enrichment factor (EF). Absolute Factor Score/Multiple Linear Regression (APCS/MLR), geographic information system (GIS), self-organizing mapping (SOM), and random forest (RF) are used for the source allocation of soil heavy metals. Finally, the combination of APCS/MLR and health risk assessment model is adopted to evaluate the risks of heavy metal sources and determine the priority-control source. The results show that the average values of soil heavy metals in the study area (Cd, Hg, As, Pb, Cr, Cu, Zn, and Ni) exceed the background values of corresponding elements in Guizhou Province. Three sources of heavy metals are identified by combining APCS/MLR, GIS, SOM, and RF. Zn (63.47%), Pb (55.77%), Cd (58.98%), Hg (32.17%), Cu (14.41%), and As (5.99%) are related to lead-zinc mining and smelting; Cr (98.14%), Ni (90.64%), Cu (76.93%), Pb (43.02%), Zn (35.22%), Cd (28.97%), Hg (22.44%), and As (5.84%) are mixed sources (natural and agricultural sources); As (88.17%), Hg (45.39%), Cd (12.04%), Cu (8.66%), and Ni (6.72%) are related to the mining and smelting of coal and iron. The results of health risk assessment show that only As poses a non-carcinogenic risk to human health. 3.31% of the sampling points of As have non-carcinogenic risks to adults and 10.22% to children. In terms of carcinogenic risks, As, Pb, and Cr pose carcinogenic risks to adults and children. Combined with APCS/MLR and the health risk assessment model, the mining and smelting of coal and iron is the priority-control pollution source. This paper provides a comprehensive method for studying the distribution of heavy metal sources in areas with large variation coefficients of soil heavy metals in the Karst region. Furthermore, it offers a theoretical basis for the management and assessment of heavy metal pollution in agricultural land in the study area, which is helpful for researchers to make strategic decisions on food security when selecting agricultural land.

Identification of macrophage-related genes in sepsis-induced ARDS using bioinformatics and machine learning

Sepsis-induced acute respiratory distress syndrome (ARDS) is one of the leading causes of death in critically ill patients, and macrophages play very important roles in the pathogenesis and treatment of sepsis-induced ARDS. The aim of this study was to screen macrophage-related biomarkers for the diagnosis and treatment of sepsis-induced ARDS by bioinformatics and machine learning algorithms. A dataset including gene expression profiles of sepsis-induced ARDS patients and healthy controls was downloaded from the gene expression omnibus database. The limma package was used to screen 325 differentially expressed genes, and enrichment analysis suggested enrichment mainly in immune-related pathways and reactive oxygen metabolism pathways. The level of immune cell infiltration was analysed using the ssGSEA method, and then 506 macrophage-related genes were screened using WGCNA; 48 showed differential expression. PPI analysis was also performed. SVM-RFE and random forest map analysis were used to screen 10 genes. Three key genes, SGK1, DYSF and MSRB1, were obtained after validation with external datasets. ROC curves suggested that all three genes had good diagnostic efficacy. The nomogram model consisting of the three genes also had good diagnostic efficacy. This study provides new targets for the early diagnosis of sepsis-induced ARDS.

Topographic effects amplify forest disturbances detected by yearly wide-time-window Landsat time series

ABSTRACT Widespread topographic effects in remote sensing images of mountainous regions with rugged terrain severely hinder satellite-based global land surface monitoring and change detection. However, topographic effects in Landsat time series (LTSs) disturbance detection remain highly debated and unquantified. In this study, we proposed a novel postprocessing approach to quantify the impact of topographic effects on long time-series (1989 − 2020) forest disturbance detection, taking the Hengduan Mountains Region (HDMR) of southwest China as an example. This approach applied a pixel-by-pixel simulation based on a semiempirical sun-canopy-sensor with a C-correction (SCS+C) model to identify and remove topography-induced disturbances from the LandTrendr-detected forest disturbances. Filtering of the detected forest disturbances with different patch sizes was conducted, and its effectiveness in removing topography-induced disturbances was quantitatively evaluated. The results showed that 10.43% of the total area of LandTrendr-detected forest disturbances was topography-induced. Removing topography-induced disturbances increased the traditional and area-adjusted overall accuracies (OAs) of the forest disturbance map by 3.50% and 0.62%. Topography-induced disturbances occurred mainly on northwest-facing (300°−360°) slopes with a gradient of 30°−55° and higher latitude and between the day of year (DOY) pairs (such as DOYs 90, 330, and 360) with a considerable disparity in solar azimuth and zenith angle. The pixel-by-pixel simulation approach outperformed spatial noise filtering with a minimum mapping unit (MMU), which considerably decreased the area percentage of topography-induced disturbances at the expense of decreasing the producer’s accuracy (PA) in detecting disturbed forests. This study highlighted the importance of removing topographic effects from the detected forest disturbances and shed new light on the promising potential of the postprocessing topographic correction method to replace the preprocessing topographic correction of the LTSs to map forest disturbances in mountainous regions globally.

Zoning Prediction and Mapping of Three-Dimensional Forest Soil Organic Carbon: A Case Study of Subtropical Forests in Southern China

Accurate soil organic carbon (SOC) maps are helpful for guiding forestry production and management. Different ecological landscape areas within a large region may have different soil–landscape relationships, so models specifically for these areas may capture these relationships more accurately than the global model for the entire study area. The aim of this study was to investigate the role of zonal modelling in predicting forest SOC and to produce highly accurate forest SOC distribution maps. The prediction objects were SOC at five soil depths (0–20, 20–40, 40–60, 60–80, and 80–100 cm). First, the forest type map and soil texture class map were used to divide the relative homogeneous regions in Shaoguan City, Guangdong Province, China. Second, seven terrain variables derived from a 12.5-m digital elevation model (DEM) and five vegetation variables generated from 10-m Sentinel-2 remote sensing images were used as predictors to develop regional artificial neural network (ANN) models for each homogeneous region, as well as a global ANN model for the entire study area (1000 sample points). Finally, 10-fold cross-validation was used to assess the ANN prediction model performance, and independent validation was used to evaluate the produced forest SOC prediction maps (194 additional samples). The cross-validation results showed that the accuracies of the regional models were better than that of the global model. Independent validation results also showed that the precision (R2) of 0- to 100-cm forest SOC maps generated by forest type modelling had an improvement of 0.05–0.15, and that by soil texture class modelling had an improvement of 0.07–0.13 compared to the map generated by the global model. In conclusion, delineating relatively homogeneous regions via simple methods can improve prediction accuracy when undertaking soil predictions over large areas, especially with complex forest landscapes. In addition, SOC in the study area is generally more abundant in broadleaf forest and clay areas, with overall levels decreasing with soil depth. Accurate SOC distribution information can provide references for fertilization and planting. Plants with particularly high soil fertility requirements may perhaps be planted in broadleaf forests or clay areas, and plants with particularly developed roots may require furrow application of a small amount of SOC.