Next Article in Journal
Innovative Seismic Imaging of the Platinum Deposits, Maseve Mine: Surface and In-Mine
Previous Article in Journal
Impact of Humidity and Freeze–Thaw Cycles on the Disintegration Rate of Coal Gangue in Cold and Arid Regions: A Case Study from Inner Mongolia, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Enhanced Weighted Drainage Catchment Basin Method for Extracting Geochemical Anomalies

1
School of Earth Sciences and Resources, China University of Geosciences (Beijing), Beijing 100083, China
2
China Railway SIYUAN Survey and Design Group Co., Ltd., Wuhan 430063, China
3
Institute of Applied Ecology, Chinese Academy of Sciences, Shenyang 110016, China
*
Author to whom correspondence should be addressed.
Minerals 2024, 14(9), 912; https://doi.org/10.3390/min14090912
Submission received: 23 July 2024 / Revised: 26 August 2024 / Accepted: 5 September 2024 / Published: 5 September 2024
(This article belongs to the Section Mineral Exploration Methods and Applications)

Abstract

:
Geochemical measurements of stream sediments are practical for small-scale mineral exploration. However, traditional grid interpolation methods cause element concentrations to diffuse and smooth out anomalies, particularly in complex terrains, making it challenging to reflect the actual distribution of elements accurately. We applied the Dynamic Enhanced Weighted Drainage Catchment Basin (DE-WDCB) method to enhance the retention and identification of local anomalies by limiting the scope of analysis to specific drainage units. This method reduces interference from varying background values across different watersheds, effectively enhancing geochemical element anomalies and aligning better with geomorphic conditions. The DE-WDCB method was tested in the Duobaoshan–Heihe area, a significant copper polymetallic mineral district in northeastern China. Compared with traditional grid interpolation methods, the DE-WDCB method retained and strengthened low and weak abnormal information of favorable mineralization elements, particularly in the Luotuowaizi area. The method demonstrated a higher spatial coverage rate with mineral points and a more vital ore-indicating ability. Specifically, the DE-WDCB method identified anomalies with a mean accuracy of 63.57% (p < 0.05, 95% CI: 47.64%–79.50%), compared to 50.53% for traditional methods. In conclusion, in regions with a complex topography and watershed differences, the DE-WDCB method effectively reduces local geochemical background interference, accurately identifies low and weak geochemical anomalies, and better reflects the actual distribution of elements. This makes it a significantly advantageous method for geochemical anomaly extraction, delineating higher-confidence exploration targets in the Sandaowan–Luotuowaizi area in the east and the triangular area between Duobaoshan, Yubaoshan, Sankuanggou, and the midstream highlands of the Guanbird River in the west.

1. Introduction

Traditional mineral exploration methods primarily rely on surface mapping to trace mineralized outcrops. However, as mineral exploration in China has advanced, the focus has gradually shifted toward locating concealed, semi-concealed, and difficult-to-identify deposits [1,2,3]. The decreasing availability of ‘high, large, and complete’ anomaly information necessitates new methods for identifying weak mineralization anomalies, which are critical for future exploration efforts [4,5,6].
With the continuous development of geological big data, the re-analysis of exploration geochemical data has become increasingly important. The extraction of geochemical anomalies from stream sediments, a vital geochemical exploration method, plays a crucial role in understanding the mineralization patterns of concealed deposits, delineating mineralization-induced anomalies, optimizing exploration targets, and evaluating resource potential [7,8,9,10]. In traditional geological work, grid-based continuous interpolation methods do not account for local variations in element background values caused by topographical differences within the area. This often results in false high-intensity anomaly zones and the smoothing or masking of weak anomaly information. Consequently, there is a discrepancy with the actual anomaly distribution, leading to the loss of significant weak anomaly information indicative of mineralization [11,12,13,14].
To address the above issues, researchers have proposed the “catchment-controlled area” method [15], which treats catchment basins as control units for stream sediment sampling points. As early as 1987, Bonham-Carter et al. [16] proposed Sample Catchment Basin Analysis (SCB) for predicting lead–zinc deposits in the Cobequid Highlands of Nova Scotia, Canada. This method demonstrated its effectiveness in improving the tracing and extraction capabilities of favorable element anomalies in mountainous and hilly regions. In 2006, Spadoni M. et al. [17] proposed the Extended Sample Catchment Basins (ESCB) method, which refines the catchment basin model by calculating sediment source concentration cut points through interpolation curves. They applied this method to extract anomalies in the southern part of the Mignone Basin in Italy. In 2013, Yousefi M. [18] proposed the Weighted Drainage Catchment Basin (WDCB) method, assigning fuzzy weights to each basin and predicting the mineral prospect of porphyry copper deposits in the Kerman Province of southeastern Iran. In 2019, Farahbakhsh E. [19] proposed the Catchment Basin Score (CBS) method, which describes the suitability of sampling points within catchment basins and improves the anomaly extraction capability of the WDCB method by adjusting watershed weights. This method successfully predicted porphyry copper–gold deposits in the Macquarie River Basin in eastern New South Wales, Australia. Over nearly four decades, methods for identifying and extracting geochemical anomalies based on catchment basins have developed into a systematic workflow that effectively avoids mathematical interference between samples and extracts local anomalies, particularly in mountainous and hilly regions [17,20,21].
The Duobaoshan–Heihe ore cluster lies in Heilongjiang Province, in northeastern China’s Greater Khingan Range. It is one of China’s three significant porphyry copper polymetallic mining areas. Situated at the Xingmeng orogenic belt intersection and the active continental margin of the western Pacific [22], the Duobaoshan area has a rich tectonic evolution history, providing favorable conditions for mineralization. The region hosts a variety of discovered deposits with complex mineralization characteristics, including super-large deposits such as the Duobaoshan porphyry copper deposit, the Tongshan porphyry copper deposit, and the Zhengguang gold deposit. Consequently, the Duobaoshan area has long been a focal point for geological exploration and mineralization theory research in Northeast China [23]. Previous geochemical surveys conducted at scales of 1:50,000 and 1:200,000 in the Duobaoshan ore cluster area have employed traditional grid-based interpolation methods to create contour maps of element anomalies [24,25]. However, these methods are unsuitable for expressing stream sediment measurements’ results, resulting in less reliable geological body reflections and highlighting the need for more appropriate mapping methods.
This study aims to evaluate the effectiveness of the Dynamic Enhanced Weighted Drainage Catchment Basin (DE-WDCB) method in accurately identifying geochemical anomalies in complex terrains, specifically within the Duobaoshan–Heihe ore cluster area. We hypothesize that the DE-WDCB method will outperform traditional grid interpolation methods by reducing background noise and enhancing anomaly detection. This study employs inverse distance weighting (IDW), Kriging, and the DE-WDCB method to extract and analyze geochemical anomalies from 1:200,000 scale stream sediment geochemical data in the Duobaoshan area. The findings are expected to provide new insights into geochemical mineralization environments and offer a scientific basis for future regional copper exploration.

2. Methodology

2.1. Multivariate Statistical Analysis

Multivariate statistical analysis is an effective method for analyzing geochemical data. Standard techniques include correlation analysis, R-type cluster analysis, and principal component analysis (PCA). In this study, we employed these methods to extract element combination characteristics, reduce the dimensionality of the measurement data, and explore the internal structure among multiple variables [26]. We used ArcGIS versiom 10.8 for spatial analysis and SPSS version 28 for statistical computations, including multivariate statistical analysis and trend surface fitting.
This approach not only investigates the co-association patterns of various elements but also characterizes the spatial distribution of primary ore-forming element combinations based on principal component scores. It can indicate the mineralization environment, discern the sources of ore-forming materials, simulate and reconstruct the element migration and evolution processes, and reveal geochemical patterns [27]. These insights provide a scientific basis for mineralization prediction and resource evaluation.

2.2. Dynamic Enhanced Weighted Drainage Catchment Basin Method

The Dynamic Enhanced Weighted Drainage Catchment Basin (DE-WDCB) method, an optimized algorithm based on the Weighted Drainage Catchment Basin (WDCB) method, was then applied [28]. Both methods focus on discretely modeling the geochemistry of the catchment basin where the sampling points are located and assigning these characteristics to their respective basins. Unlike the WDCB method, the DE-WDCB method considers variations in the geochemical background field at different locations. It revises the WDCB method’s approach by employing trend surface analysis theory to fit different geochemical background values for each basin unit and grades anomalies for the basins of different sampling points to facilitate anomaly extraction.
Due to the influence of various activities, such as the transportation, migration, and mixing of surface water and groundwater in upstream basins, the concentration of indicative elements tends to decrease with increasing distance from the anomaly source downstream. Thus, strong anomalies often exhibit sliding displacement along the stream distribution in practical geochemical anomaly extraction. However, the complexity and uncertainty of downstream influences mean that the relationship between indicative element concentrations and distance downstream is nonlinear. Therefore, the DE-WDCB method, building on the theoretical foundation of the WDCB method, uses the “suitability degree allocation” weighting method to fit this nonlinear relationship. It introduces “accumulated flow volume” as a parameter to represent the suitability of each sampling point’s basin. The larger the accumulated flow volume, the more likely the location is downstream and prone to depositing elements transported from upstream. Therefore, it should be assigned a smaller weight in anomaly calculations to reduce false anomalies. Conversely, elements in relatively upstream locations are more likely to be carried downstream by water flow and should be assigned a more significant weight to restore the upstream anomaly quantity. Ultimately, this enhances the ability to identify and extract indicative element anomaly sources and accurately reconstruct the anomaly levels in the catchment basin.
The DE-WDCB method for extracting geochemical anomalies mainly involves the following two parts:
Using trend surface analysis, the DE-WDCB method determines the basin background value and anomaly threshold. The following equation expresses the mathematical formula:
Zi = Ti + Mi
where Zi is the observed value, Ti is the regional trend value, and Mi is the residual component, consisting of local variations and random values. During geochemical anomaly extraction, areas where the residual component exceeds the trend (Zi-Ti > 0) are significant as they often indicate regions of element enrichment.
The DE-WDCB method assigns anomaly scores to each drainage basin. The mathematical formula is given by:
C x = i = 1 n a F A i i = 1 n a F A i · n i + 1 M i N i N t · 100
where Cx is the DE-WDCB score, n is the total number of anomaly grades for sampling point elements, FAi is the average value of each type of flow accumulation, Mi is the median concentration of elements classified as the highest anomaly grade (i-grade), Ni is the number of sediment samples in the i-grade, and Nt is the total number of sediment samples in each drainage basin unit. a is an introduced constant with a value between 0 and 1.
The DE-WDCB method can employ various grading rules for anomaly classification. In this study, the cumulative frequency method was used for grading to compare with the anomaly extraction results of the control group. Thresholds for low, medium, and high anomalies were set at the 85th, 95th, and 98th percentiles, respectively, i.e., n = 3.

3. Workflow

Figure 1 illustrates the detailed workflow of this study, outlining the steps from data collection to anomaly extraction and analysis. The workflow is structured as follows:
  • Data collection and preprocessing using ArcGIS.
  • Statistical analysis and extraction of favorable metallogenic elements using SPSS.
  • Establishment of the Digital Watershed Model.
  • Application of the DE-WDCB method for anomaly extraction.
  • Interpretation of results with a focus on spatial distribution and geochemical patterns.
By integrating these advanced statistical techniques and geospatial analysis tools, we aimed to improve the accuracy and reliability of geochemical anomaly detection in complex terrains.

4. Application Examples

4.1. Geological Setting

4.1.1. Regional Geological Setting

The Duobaoshan–Heihe ore cluster area is located in the Greater Khingan Range arc-basin system, on the west side of the significant fault separating the Xing’an block and the Songnen block (Figure 2a), and belongs to the eastern end of the Central Asian–Xingmeng geosyncline fold belt [29,30]. The region has undergone several geological events, including the formation of the Duobaoshan island arc in the Early Paleozoic, the development of the Handaqi intra-arc rift in the Late Paleozoic, and intense Mesozoic basin-and-range tectonic and magmatic activities (Figure 2b) [31]. These long-term tectonic and magmatic activities have resulted in the formation of the Duobaoshan–Kuanhe rhombohedral structure, controlled by NE- and NW-trending significant faults [32], which significantly influence the distribution of strata, magmatic rocks, and mineral deposits within the area. The sedimentary strata in the area primarily include the Ordovician, Silurian, and small amounts of the Upper Proterozoic, Cambrian, Devonian, Carboniferous, Permian, Jurassic, Cretaceous, and Cenozoic.
The ore cluster area mainly hosts several super-large to large Early Paleozoic epithermal low-temperature hydrothermal copper–molybdenum and copper–gold deposits, such as the Duobaoshan copper deposit, the Tongshan copper deposit, and the Zhengguang gold deposit [33]. The periphery of the ore cluster area predominantly develops Mesozoic magmatic–hydrothermal copper–gold, copper–molybdenum, and copper–lead–zinc deposits (or occurrences), such as the Sankuangou copper–molybdenum deposit and the Xiaoduobaoshan copper deposit. The Duobaoshan-type copper polymetallic deposits primarily formed due to the combined effects of early Caledonian tectonic and magmatic activities. During the Yanshanian period, hydrothermal activities were most intense, significantly impacting mineral enrichment and the destruction of ore deposits in the area [34].
Figure 2. (a) Schematic diagram of the Eurasian continent (adapted from [35]). (b) Tectonic location of the Duobaoshan–Heihe metallogenic belt.
Figure 2. (a) Schematic diagram of the Eurasian continent (adapted from [35]). (b) Tectonic location of the Duobaoshan–Heihe metallogenic belt.
Minerals 14 00912 g002

4.1.2. Overview of Typical Duobaoshan-Type Copper Deposits in the Study Area

The strata in the study area are well developed, ranging from the Ordovician to the Quaternary. The Upper Silurian Wotu Formation (S3w) unconformably overlies the Lower Devonian Heitai Formation (D2h), whereas the other strata maintain conformable contacts (Figure 3). Typical super-large copper deposits in the area, such as Duobaoshan and Tongshan, mainly occur in the Lower Ordovician Tongshan Formation (O1t) and the Middle Ordovician Duobaoshan Formation (O2d), consisting of andesite, andesitic tuff, tuffaceous sandstone, quartz sandstone, siltstone, and carbonaceous slate [36].
The area has experienced frequent magmatic activities, with multiple magmatic intrusions during the Jinning, Caledonian, Variscan, Late Indosinian, and Yanshanian periods. Granite intrusions spread widely across the region. The Early Caledonian tectonic layer, comprising Neoarchean to Early Cambrian volcanic and carbonate clastic rocks, forms the foundation of the Duobaoshan mineralization area [37]. The granitic porphyry and granodiorite of the Duobaoshan volcanic arc are the central ore-bearing rock bodies in the area.
The structural features in the area are well developed, with NE- and NW-trending faults jointly controlling the Yanshanian mineralization belts. The NW-trending Duobaoshan–Taxi fault and the Huasuizhai–Kesu Gulan River fault, formed in the Late Paleozoic, control the distribution of strata in the area [31]. Copper mineralization is primarily controlled by faults formed by submarine volcanic rocks of the Duobaoshan Formation and associated NE- and NW-trending structures.
The alteration associated with the porphyry copper deposits and granitic porphyry in the area exhibits distinct zoning characteristics, forming a planar, progressive porphyry–copper alteration zone. This zone extends outward from the central potassic silicification, followed by zones of potassic feldspar porphyry, quartz–sericite, and chlorination, with localized occurrences of monetization, skarn, and hornfels alterations.

4.2. Basic Data Characteristics and Multivariate Statistical Analysis

This study utilizes 1:200,000 scale stream sediment measurement data from the China National Geochemical Mapping Project to statistically analyze 14 elements (Au, As, Ag, Cd, Co, Cr, Cu, Nb, Ni, Pb, Sb, Ti, W, Zn) in the study area.

4.2.1. Basic Characteristics of Single-Element Distribution

After undergoing a single or multiple natural geochemical processes, the element content in the geological body tends to be normally distributed [38]. The frequency histogram of element concentration (Figure 4) illustrates the distribution pattern of each element. Through the iteration method, we confirmed that the data of all 14 elements passed the regular distribution test. Then, we calculated these elements’ mean, standard deviation, anomaly classification threshold, and other related parameters to understand the spatial distribution and enrichment or depletion.
The statistical analysis in Table 1 reveals that, except Sb and Cd, the average concentrations of other elements in this region exceed the background levels of the Greater Khingan metallogenic belt, indicating significant enrichment of Au, As, Ag, Co, Cr, Cu, Nb, Ni, Pb, Ti, W, and Zn. Notably, Cu, W, Ag, and Au exhibit high variation coefficients (>1.50) and concentration coefficients (>1.3), suggesting local solid differentiation and marking them as key ore-forming elements. Mineralization points for Cu, Ag, and Au have been confirmed in the study area, reinforcing its potential for further exploration.

4.2.2. Characteristics of Element Combinations

This study employs correlation analysis, cluster analysis, and principal component analysis (PCA) to investigate the correlations between elements and identify indicative element combinations for copper polymetallic deposits in the study area.
1.
Cu as the Primary Ore-Forming Element and its Indicative Correlations
To address compositional data’s “closure effect” [40], we standardized the data using a log transformation before calculating Pearson correlation coefficients. As shown in Figure 5, Cu exhibits a significant positive correlation with Co, Cr, and Ni and a significant negative correlation with Ti, all with absolute coefficients greater than 0.5. These findings suggest that Cu is the primary ore-forming element in the study area, while Co, Cr, Ni, and Ti are likely indicative elements for mineral exploration.
2.
Element Grouping Revealing Copper Polymetallic Mineralization Patterns
We performed an R-type clustering analysis on 14 elements and divided them into four groups at a Euclidean distance of 18 (Figure 6). The first and third groups present the combinations closely related to copper polymetallic mineralization.
  • Group 1: Cu, Cr, Ni, Co, Ti
This group is closely related to copper polymetallic mineralization. As typical chalcophile elements, Cu, Ni, and Co are strongly associated with intermediate to acidic magmatic–hydrothermal activities, often enriching and precipitating in high-temperature hydrothermal environments. Cr, primarily a lithophile element, is typically found in the outer alteration zones formed during the quartz–sericite stage of mineralization. Ti, mainly present in ilmenite (FeTiO3) and rutile (TiO2), is active in high-temperature, high-pressure, and acidic hydrothermal fluids but remains stable and precipitates in lower-temperature environments [41]. It contrasts with the accumulation and deposition environments of Cu, Co, and Cr. This group’s distribution makes it a critical indicator for locating and assessing copper polymetallic deposits.
  • Group 3: Cd, Zn, Pb
These elements, also chalcophile in nature, typically coexist with copper polymetallic minerals and are associated with medium to high-temperature hydrothermal mineralization. They often deposit in lower-temperature environments, typically on the peripheries of copper ore bodies, exhibiting distinct primary halo zoning. This distribution pattern indicates that these elements have undergone hydrothermal mineralization under varying temperature conditions, suggesting the area’s multi-staged or complex geological processes.
  • Group 4: As, Sb, Au, Ag
This group indicates low-temperature hydrothermal mineralization, notably different from the high-temperature magmatic activities indicated by the first group. We can view these variations in temperature and pressure conditions as markers of multiple geological processes in the Duobaoshan–Heihe area.
3.
High-Temperature Element Combination Reflecting Copper Mineralization Potential
According to the significance test (Table 2), the KMO value is 0.733, the degree of freedom, df, is much larger than 0, and the significance is less than 0.01, which aligns with the requirement of using principal component analysis.
The principal component analysis extracted four component factors with eigenvalues greater than 1 (Table 3).
  • Fac 1: This factor is dominated by high-temperature elements, including Co, Cr, Cu, Ni, and Ti. It primarily reflects copper mineralization information, indicating the substantial potential for copper exploration in the area.
  • Fac 3: This factor is composed of Cd and Zn, elements typically enriched in medium to low-temperature metallogenic zones. It reflects the primary halo zoning phenomena exhibited by wall rock alteration during various stages of copper mineralization.

4.3. Extraction of Stream Sediment Geochemical Anomalies Using DE-WDCB

This section describes the extraction of watershed characteristics in the study area using Aster 30m resolution DEM data to establish the geochemical landscape of catchment basins (Figure 7 and Figure 8). Using the DE-WDCB method, we determined dynamic background values and anomaly thresholds for different sampling points, effectively addressing geochemical background differences between basin landscapes.

4.3.1. Construction of the Catchment Basin Landscape Model

Before constructing the catchment basin model, the original data were preprocessed with sink-filling analysis to restore the correct outlet positions of the basins, ensuring the continuity of the generated water systems and the accuracy of the model (Figure 8a). Next, using the D8 algorithm, the flow direction was set to eight possible directions, represented by 2n (where n ranges from 0 to 7). ArcGIS calculates the slope relationships between each grid cell and its eight neighboring cells to determine the steepest downhill neighbor as the flow direction (Figure 8b). We refer to the flow accumulation of all downhill cells flowing into a grid cell as the flow accumulation of that cell (Figure 8c). Higher flow accumulation indicates more concentrated flow and a greater likelihood of runoff formation, reflecting the transport trend of weathering products like rocks across the surface.
We used the Strahler classification method to classify the water systems in the area into four levels (Figure 8d). Rivers originating from the headwaters are classified as first-order rivers, which, as shown in the figure, have a broad influence area, indicating minimal pollution from natural and human factors in the survey area. Finally, we set the outlets of the various water system levels as the spill points for the catchment basins. By capturing these spill points, a catchment basin model matching the water system was generated (Figure 8e,f).

4.3.2. Extraction and Mapping of Geochemical Anomalies for Favorable Ore-Forming Elements and Element Combinations

We used dynamic trend surface analysis to calculate trend values for each catchment basin based on the spatial distribution of high indicator element contents and the catchment basin divisions in the study area. This approach allowed us to obtain the distribution of geochemical background fields in the study area.
The results of the background field extraction (Figure 9) show significant differences in the geochemical background fields of Co, Cr, Cu, Ni, Ti, and their combination elements across different basin landscapes in the study area. The study area’s northwest to west side shows a high trend value concentration for favorable ore-forming single elements. In contrast, low trend values are concentrated in the central to eastern side, with the overall distribution trend weakening from northwest to southeast. The western and northeastern edges of the study area exhibit high trend values for element combinations, with small high-value areas showing a weak distribution trend from southwest to northeast.
The DE-WDCB algorithm’s inversion results reveal that single-element and element combination anomalies predominantly occur in the Duobaoshan–Sankuangou and Luotuowaizi–Sandaowan areas (Figure 10). The anomalies form closely connected block-point patterns in these regions, with primary ore-forming Cu and element combinations showing similar distribution characteristics.
In the Duobaoshan–Sankuangou area on the study area’s west side, high-value Cu and element combination anomalies extend in an NWW direction, forming an elliptical distribution over a large area. The anomaly area exposes a significant amount of favorable ore-forming strata, such as the Duobaoshan Formation (O2d) and the Tongshan Formation (O1-2t), and includes multiple complex NNW-trending faults. These faults generally align with the extension direction of the Sankuangou–Duobaoshan–Luohe ore-controlling fold belt (Figure 2). The anomaly area covers two super-large copper deposits, Duobaoshan and Tongshan, as well as 11 medium to small copper–molybdenum and copper–lead–zinc deposits and occurrences such as Xiaoduobaoshan and Sankuangou.
High anomalies cluster in the Sandaowan area on the northeastern side of the study area, forming three dispersed high-concentration centers of Cu elements. The overall distribution area is large, with numerous NNE-trending secondary faults of the Nenjiang-Xinkailing lithospheric fault zone extending within the anomaly area.

4.4. Comprehensive Comparison and Evaluation of Anomaly Extraction Results

The study area has a topography characterized by high central and low eastern and western regions, with the northern ridge of the Xiaoxing’an Mountains cutting through the center, dividing the study area into eastern and western regions (Figure 11). Data and geological analyses were conducted separately for these regions (Figure 12).
This study used inverse distance weighting (IDW) and Kriging methods in the control group to process element content data from geochemical sampling points. Table 4 displays the optimal interpolation model parameters.
The cumulative frequency method was used for grading in the IDW, Kriging, and DE-WDCB methods, with the 85th, 95th, and 98th percentiles as thresholds for low, medium, and high anomalies, respectively, to extract single elements and element combinations related to copper mineralization (Table 5).

4.4.1. Analysis of Traditional Anomaly Extraction Results Based on Grid Interpolation—IDW and Kriging

1.
Favorable Ore-Forming Single Elements
The control group, employing IDW and Kriging methods (Table 6), identified Cu, Co, Cr, and Ni anomalies primarily concentrated in the Duobaoshan–Sankuangou area (Region A), with a spatial distribution characterized by higher values in the west and lower values in the east. These anomalies encompassed 14.21% of the study area and overlapped with 48.21% (p < 0.05, 95% CI: 31.61%–64.81%) of known copper mineralization points, with an average Cu anomaly coverage of 48.57% (p < 0.05, 95% CI: 31.98%–65.16%). Significant deposits, such as Duobaoshan, Xiaoduobaoshan, and Tongshan, were mainly associated with medium to low anomaly levels.
Ti anomalies, which negatively correlate with copper mineralization, were predominantly located in Region A’s northern and southern parts. These anomalies corresponded closely with the shale, gneiss, and grey mudstone distribution, exhibiting a distinct negative correlation with Cu, Co, Cr, and Ni anomalies.
Kriging outperformed IDW slightly, providing smoother anomaly boundaries and a closer alignment with NNW-trending faults. The anomaly area identified by Kriging was 0.98 times that of IDW, while its mineral point coverage was 1.14 times higher.
2.
Favorable Ore-Forming Element Combinations
The element combination anomalies of Co-Cr-Cu-Ni-Ti extracted by both methods in the control group show a similar spatial distribution to the Cu anomalies. Anomalies are only identified in the Duobaoshan–Sankuangou area (Region A), mainly exhibiting an elliptical distribution along the NW direction, controlled by the Duobaoshan–Luohe structural belt. The average anomaly area covers 12.83% of the total area, with an average mineral point coverage rate of 52.85% (p < 0.05, 95% CI: 36.23%–69.47%).

4.4.2. Analysis of Anomaly Extraction Results Based on the Discrete Interpolation of the Catchment Basin Model—DE-WDCB

1.
Favorable Ore-Forming Single Elements
Compared to the control group results, the DE-WDCB method extracts high to medium anomalies of Cu, Co, Ni, and Ti over a similar proportion of the total study area while significantly increasing the area of low anomalies. The average three-level anomalies cover 28.34% of the study area, 1.99 times higher than the control group, with an average mineral point coverage rate of 55.71% (p < 0.05, 95% CI: 39.25%–72.17%), which is 1.16 times higher than the control group. Specifically, the coverage rate of Cu anomalies over mineral points is 57.14% (p < 0.05, 95% CI: 50.75%–73.53%), an increase of 1.18 times, with typical copper deposits like Duobaoshan, Tongshan, and Xiaoduobaoshan located in high anomaly areas.
The DE-WDCB-extracted Cu, Co, Cr, and Ni anomalies show a general trend of low values in the central area and high values in the eastern and western regions, consistent with the geological body distribution. In the Duobaoshan–Sankuangou area (Region A), an NWW-trending high anomaly zone aligns with the Duobaoshan–Luohe structural belt, indicating strong structural control over copper mineralization. In the Sandaowan–Songshugangzi area (Region B), NEE-trending anomalies around the Luotuowaizi copper–molybdenum deposit and low anomalies along NEE-trending secondary faults suggest robust structural control on the eastern side.
Analyzing Co, Cr, and Ni anomalies reveals a shift in high anomaly concentration in Region A along an NE-SW axis, reflecting the upstream–downstream water system distribution. In the northern part of Region A (Independence Mountain), where the control group identified no anomalies, the DE-WDCB method shows low anomalies within multiple first-order basins. Ti anomalies, negatively correlated with copper mineralization, mainly appear as widespread low anomalies that do not overlap with copper-related mineral points, highlighting their distinct negative correlation.
2.
Favorable Ore-Forming Element Combinations
Compared to the traditional methods, the DE-WDCB method extracts high and medium anomalies of favorable ore-forming element combinations over a similar proportion of the study area but significantly increases the area of low anomalies. The total anomalies cover 27.48% of the study area, with a mineral point coverage rate of 71.42% (p < 0.05, 95% CI: 56.45%–86.39%), 2.29 times the anomaly area and 1.4 times the mineral point coverage rate of the IDW method. The overall distribution of element combination anomalies shows higher values in the east and west and lower values in the center, similar to the distribution of single-element anomalies. Unlike the control group, the DE-WDCB method identifies anomalies in Regions A and B, with significant changes in anomaly locations.
In Region A, high and medium anomaly information is mainly concentrated in the Duobaoshan–Sankuangou area. Compared to the control group, the DE-WDCB extracted medium anomalies cover large copper deposits such as Duobaoshan, Tongshan, Xiaoduobaoshan, and Nangou copper–molybdenum deposits more extensively. It also covers four copper–lead–zinc mineral points on the southwestern side of Duobaoshan, indicating more substantial mineralization potential. The low anomaly concentration centers around the Zhan copper point northeast of Xiaoduobaoshan have shifted along an NE-SW axis, and researchers identified many low anomalies in the first-order basin areas around Independence Mountain.
Region B has fewer independently distributed basins with high and medium anomalies, while basins with low anomalies are more numerous and distinctly zoned. In the Sandaowan–Luotuowaizi–Songshugangzi area, significant low anomalies extend in an NE direction, consistent with the extension direction of secondary faults of the Xinkailing lithospheric fault zone, indicating robust structural control features.

4.4.3. Comparative Analysis and Advantages of the DE-WDCB Method over Traditional Approaches

Traditional methods rely on uniform background reference values, which can effectively extract distinct third-level anomalies in regions with prominent high anomaly clusters, such as Area A. However, the pronounced topographic differences between the eastern and western ends of the Duobaoshan–Heihe study area present significant challenges. Using excessively high background reference values in these methods often leads to the diffusion and smoothing of relatively weak anomaly signals in Area B. This can result in misrepresentation, where areas in Area B that should display anomalies instead appear as “non-anomalous”, thereby hindering an accurate depiction of the actual elemental distribution.
In contrast, the DE-WDCB method addresses these challenges by accounting for geochemical background variations across different basin units. By integrating trend surface analysis with the “suitability allocation” weighted scoring method, the DE-WDCB approach can accurately align background values to basin units according to their specific geomorphic characteristics.
The comparative analysis demonstrates that the DE-WDCB method significantly outperforms traditional approaches. The average anomaly mineral point coverage using the DE-WDCB method was 63.57% (p < 0.05, 95% CI: 47.64%–79.50%), compared to 50.53% (p < 0.05, 95% CI: 33.97%–67.09%) for traditional methods. This improvement underscores the DE-WDCB method’s enhanced ability to preserve and amplify low and weak anomalies, particularly in Area B. Consequently, this approach achieves a more significant spatial overlap between the extracted anomalies and known mineral points, highlighting its superiority in accurately reflecting elemental distribution, especially in regions with complex topography and variable geochemical backgrounds.

5. Discussion

5.1. Interpretation of Results

Our experiments revealed that the traditional IDW and Kriging methods yield results comparable to the DE-WDCB method in identifying anomalies in the western part of the study area (Region A). All three methods successfully detected large-scale high-to-moderate anomalies of single elements and element combinations in the Duobaoshan–Sankuangou region (Figure 13). These anomalies trend in an NNW direction, consistent with the extension of the Duobaoshan–Luohe fault zone, and align well with the geological formations such as the Upper Carboniferous granite porphyry, Lower Ordovician Tongshan Formation (O1t) andesite, and Middle Ordovician Duobaoshan Formation (O2d) granite porphyry. This finding validated the effectiveness of these methods in environments with high local element concentrations.
However, in the eastern part of the study area (Region B)(Figure 14), traditional methods like IDW and Kriging identified only minor, patchy low to weak Cu anomalies, which did not spatially coincide with known mineral deposits. This discrepancy indicates the susceptibility of traditional grid interpolation methods to local background value differences in areas with lower element concentrations, thereby limiting their effectiveness in uniform anomaly extraction.
Conversely, the DE-WDCB method demonstrated superior performance in Region B by effectively reducing the interference of complex topography through trend surface analysis. This method minimized the smoothing effect of large non-anomalous areas, successfully identifying low to weak anomalies in the Sandaowan–Luotuowaizi-Songshugangzi region. These NE-trending anomalies, aligned with secondary faults in the Xinkailing lithosphere fault zone, exhibited strong spatial congruence with known copper–molybdenum deposits, underscoring the DE-WDCB method’s enhanced sensitivity and accuracy in detecting low to weak anomalies.

5.2. Comparison with Previous Studies

The geological understanding of the Duobaoshan area obtained from this study is consistent with the results of Liu B.S. et al. (2020) [30]. At the same time, the results of this study also confirmed the research of Zhu R.W. (2022) [28] on the ability of the DE-WDCB method to extract weak anomalies. They also emphasized the advantages of basin-based anomaly detection methods. The DE-WDCB method can retain and enhance low to weak anomalies, especially in areas with complex terrain, which confirms its superiority over traditional grid-based methods. Unlike traditional methods that tend to smooth and distort the distribution of anomalies, the DE-WDCB method can describe more accurately and in detail element anomalies, especially in challenging geological environments.

5.3. Study Limitations

Despite the promising results, the DE-WDCB method has certain limitations. Its effectiveness depends on the division of watershed units and the fitting of background values for basin units. While advantageous in regions with complex topography and significant watershed differences, the method’s benefits may diminish in areas with flat terrain or minimal watershed variation. In such cases, the results might not differ significantly from those obtained through traditional methods. Additionally, the success of the DE-WDCB method relies heavily on the density and uniformity of sampling points. Sparse or unevenly distributed sampling can lead to excessive smoothing or deviations from the actual geological background, potentially obscuring critical low and weak anomaly information.

5.4. Future Research Directions

Future research should focus on refining the DE-WDCB method to address its limitations. Specifically, improving the method’s applicability in flat terrain and areas with minimal watershed variation is crucial. Enhancing the trend surface analysis component to perform reliably even with sparse or unevenly distributed sampling points could further increase the method’s accuracy. Integrating more comprehensive geological data and exploring advanced statistical techniques could contribute to more reliable anomaly detection. Expanding the application of the DE-WDCB method to various geological settings and comparing its performance with emerging anomaly detection techniques will further validate and refine its use in mineral exploration.

6. Conclusions

  • The DE-WDCB method significantly enhances the detection of geochemical anomalies in the complex terrain of the Duobaoshan–Heihe area. Our study identifies two distinct geochemical element groups indicative of different mineralization environments: one associated with medium-acid, high-temperature hydrothermal activity (Cu, Cr, Ni, Co, Ti) and the other with low-temperature hydrothermal processes (As, Sb, Au, Ag). Among these, Cu exhibits strong localized enrichment and spatial differentiation, making it a key indicator for copper polymetallic deposits.
  • Traditional methods like IDW and Kriging effectively detected high to moderate anomalies in the western Duobaoshan–Sankuangou area but failed to identify anomalies in the eastern regions. Traditional methods’ average anomaly mineral point coverage was 50.53% (p < 0.05, 95% CI: 33.97%–67.09%). In contrast, the DE-WDCB method, by accounting for topographic complexities and applying trend surface analysis, successfully detected and amplified low to weak anomalies, particularly in the Luotuowaizi area. The DE-WDCB method showed an average anomaly mineral point coverage of 63.57% (p < 0.05, 95% CI: 47.64%–79.50%). This method demonstrated superior spatial coverage and a stronger correlation with known mineralization points, underscoring its effectiveness in regions with complex topographic features.
  • Future exploration efforts should prioritize the DE-WDCB method, especially in areas like the triangular region between Duobaoshan, Yubaoshan, and Sankuangou (Region A) and the Sandaowan–Luotuowaizi area (Region B). In Region B, the exploration should focus on avoiding the Ti anomaly-rich Songshugangzi area and instead on the higher-confidence zones identified by the DE-WDCB method. Moreover, we believe that future research should focus on expanding the application of the DE-WDCB method to diverse geological settings, particularly in areas with varied geomorphological and geochemical characteristics, to validate its effectiveness further. Integrating advanced geospatial data, such as remote sensing images, could enhance anomaly detection accuracy. Optimizing sampling strategies by exploring the impact of different sampling densities and distributions will also improve geochemical anomaly extraction.

Author Contributions

Z.C. and R.Z. designed the project; Z.C. conducted the original literature reviews; funding acquisition and project management by J.C.; Z.C. and Q.Z. carried out the experiment; Z.C. wrote and organized the paper with a careful discussion and revision by Q.Z., G.Z., Z.J. and C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Key Research and Development Program of China project “Quantitative Assessment of Resource Potential in Key Metallogenic Zones of Strategic Minerals” (No. 2023YFC2906404).

Data Availability Statement

The data presented in this study are not publicly available because they are confidential. We want to ensure the confidentiality and integrity of sensitive data related to proprietary exploration techniques and geographic information that may compromise ongoing and future research projects. However, we are committed to providing legitimate researchers access to the data upon reasonable request, ensuring that the data can be independently verified and utilized for further scientific exploration.

Acknowledgments

The authors acknowledge Xiao Keyan for providing geochemical sampling point data.

Conflicts of Interest

The co-author, Renwei Zhu, is affiliated with the Company China Railway Fourth Survey and Design Institute Group Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Wang, Q.Y.; Hu, Y.P. Considerations on the Scarcity of Metal Resources and the Exploration of Concealed Deposits. Geol. Prospect. 2004, 75–79. [Google Scholar] [CrossRef]
  2. Zhao, P.D. Digital Prospecting and Quantitative Evaluation in the Era of Big Data. Geol. Bull. China 2015, 34, 1255–1259. [Google Scholar]
  3. Chen, J.P.; Cui, N.; Zhu, X.T.; Zhang, Y.; Xiang, J. Assessment of Copper Deposit Potential in China; Geological Publishing House: Beijing, China, 2017; pp. 217–226. [Google Scholar]
  4. Zuo, R.G. Geochemical Data Mining and Weak Anomaly Identification in Exploration. Earth Sci. Front. 2019, 26, 67–75. [Google Scholar]
  5. Zhao, P.D. Quantitative prediction and evaluation of solid mineral resources. In Digital Geology; Science Press: Beijing, China, 2024; pp. 380–384. [Google Scholar]
  6. Liu, B.; Cui, X.; Wang, X. The Delineation of Copper Geochemical Blocks and the Identification of Ore-Related Anomalies Using Singularity Analysis of Stream Sediment Geochemical Data in the Middle and Lower Reaches of the Yangtze River and Its Adjacent Areas, China. Minerals 2023, 13, 1397. [Google Scholar] [CrossRef]
  7. Wang, R.T.; Mao, J.W.; Ren, X.H.; Wang, J.Y.; Ouyang, J.P.; Yuan, B.Q. Current Status and Problems in the Evaluation of Regional Geochemical Anomalies. Geol. China 2005, 168–175. [Google Scholar] [CrossRef]
  8. Weng, W.F.; Wang, D.E.; Wang, B.M.; Ding, Y.; Wang, Y.J. Geochemical Characteristics and Prospecting Direction of Stream Sediments in the Qimen-Yixian Area, Anhui Province. Geophys. Geochem. Explor. 2020, 44, 1–12. [Google Scholar]
  9. Huang, W.B.; Luo, X.R.; Liu, P.F.; Zheng, C.J.; He, W.; Yang, X.X.; Xiao, X.Q.; Wang, S.L. Geochemical Characteristics and Prospecting Prediction of Stream Sediment Measurements in the Shihuigou Area, Qinghai Province. Bull. Geol. Sci. Technol. 2020, 39, 150–159. [Google Scholar]
  10. Zeng, K.; Liu, H.; Huang, D.J.; Guo, W.; Qi, S.L.; Si, X.H.; Yang, Y.Z. An Analysis of the Anomalous Characteristics and Prospecting Effects of 1:50,000 Stream Sediment Measurements in the Mengweng Area, Yunnan Province. Mod. Geol. 2021, 35, 270–280. [Google Scholar]
  11. Rose, A.W.; Dahlberg, E.C.; Keith, M.L. A Multiple Regression Technique for Adjusting Background Values in Stream Sediment Geochemistry. Econ. Geol. 1970, 65, 156–165. [Google Scholar] [CrossRef]
  12. Hawkes, H.E. The Downstream Dilution of Stream Sediment Anomalies. J. Geochem. Explor. 1976, 6, 345–358. [Google Scholar] [CrossRef]
  13. Zou, R.; Wang, J.; Chen, G.; Yang, M. Identification of Weak Anomalies: A Multifractal Perspective. J. Geochem. Explor. 2015, 148, 12–24. [Google Scholar]
  14. Huang, X.K.; Wei, J.H.; Shi, W.J.; Zhang, X.M.; Gao, Q.; Wang, S. Identification and Evaluation of Geochemical Anomalies Based on Catchment Basins: A Case Study of 1:50,000 Stream Sediment Geochemical Measurements in the Wulastai Area, East Kunlun. Bull. Geol. Sci. Technol. 2023, 42, 324–338. [Google Scholar]
  15. Zhang, W.L. Study on the Extraction Methods of Stream Sediment Information in the Duolong Ore Cluster Area, Tibet. Master’s Thesis, Chengdu University of Technology, Chengdu, China, 2018. [Google Scholar]
  16. Bonham-Carter, G.F.; Rogers, P.J.; Ellwood, D.J. Catchment Basin Analysis Applied to Surficial Geochemical Data, Cobequid Highlands, Nova Scotia. J. Geochem. Explor. 1987, 29, 259–278. [Google Scholar] [CrossRef]
  17. Spadoni, M. Geochemical Mapping Using a Geomorphologic Approach Based on Catchments. J. Geochem. Explor. 2006, 90, 183–196. [Google Scholar] [CrossRef]
  18. Yousefi, M.; Carranza, E.J.M.; Kamkar-Rouhani, A. Weighted Drainage Catchment Basin Mapping of Geochemical Anomalies Using Stream Sediment Data for Mineral Potential Modeling. J. Geochem. Explor. 2013, 128, 88–96. [Google Scholar] [CrossRef]
  19. Farahbakhsh, E.; Chandra, R.; Eslamkish, T.; Müller, R.D. Modeling Geochemical Anomalies of Stream Sediment Data through a Weighted Drainage Catchment Basin Method for Detecting Porphyry Cu-Au Mineralization. J. Geochem. Explor. 2019, 204, 12–32. [Google Scholar] [CrossRef]
  20. Wu, J.Y. Comparative Study on the Extraction and Method Effectiveness of Geochemical Anomalies in the Lanping-Simao Area. Master’s Thesis, China University of Geosciences (Beijing), Beijing, China, 2020. [Google Scholar]
  21. Kong, Y.H. Modeling Geochemical Element Migration Based on Catchment Basins: A Case Study of the Jiama Mining Area in Tibet. Master’s Thesis, Chengdu University of Technology, Chengdu, China, 2021. [Google Scholar]
  22. Wang, L.; Qin, K.Z.; Pang, X.Y.; Song, G.X.; Jin, L.Y.; Li, G.M.; Zhao, C. Geological characteristics and alteration zoning of the Tongshan porphyry copper deposit in Duobaoshan ore field: Implications for hydrothermal-mineralization centers and deep exploration. Miner. Depos. 2017, 36, 1143–1168. [Google Scholar]
  23. Li, C.L.; Fu, A.Z.; Xu, W.X.; Yuan, M.W.; Liu, B.S.; Yang, W.P.; Zhao, R.J.; Zhao, Z.H. Characteristics of Structural Overlap Halos and Deep Prospecting Prediction of the Yongxin Gold Deposit in the Duobaoshan Area, Heilongjiang Province. Mod. Geol. 2023, 37, 674–689. [Google Scholar]
  24. Xu, D.H.; Wan, T.P.; Shi, G.M. Geochemical Characteristics and Metallogenic Prospect Zoning of Stream Sediments in the Duobaoshan Area, Heilongjiang Province. Gold 2019, 40, 18–22. [Google Scholar]
  25. Zhang, L. Comparative Study on Geochemical Exploration Methods of 1:50000 Scale in the Duobaoshan Area, Heilongjiang Province. Geophys. Geochem. Explor. Comput. Technol. 2022, 44, 525–532. [Google Scholar]
  26. Singer, D.A.; Kouda, R. Some Simple Guides to Finding Useful Information in Exploration Geochemical Data. Nat. Resour. Res. 2001, 10, 137–147. [Google Scholar] [CrossRef]
  27. Cheng, Q.; Bonham-Carter, G.; Wang, W.; Zhang, S.; Li, W.; Qinglin, X. A Spatially Weighted Principal Component Analysis for Multi-Element Geochemical Data for Mapping Locations of Felsic Intrusions in the Gejiu Mineral District of Yunnan, China. Comput. Geosci. 2011, 37, 662–669. [Google Scholar] [CrossRef]
  28. Zhu, R.W. Research on Enhanced Extraction Methods for Geochemical Mineralization Anomalies in Stream Sediments. Master’s Thesis, China University of Geosciences (Beijing), Beijing, China, 2022. [Google Scholar]
  29. Ge, W.C.; Wu, F.Y.; Zhou, C.Y.; Zhang, J.H. Metallogenic Epoch and Geodynamic Significance of Porphyry Cu and Mo Deposits in the Eastern Segment of the Xingmeng Orogenic Belt. Chin. Sci. Bull. 2007, 52, 2407–2417. [Google Scholar] [CrossRef]
  30. Liu, B.S. Multiphase Mineralization and Superimposed Transformation of the Duobaoshan Porphyry Copper Deposit in Heilongjiang. Geol. Rev. 2020, 66, 29–32. [Google Scholar]
  31. Yang, X.P.; Ma, J.S.; Pang, X.J.; Yang, Y.J.; Jiang, B.; Fu, J.Y. Reconstruction of the Early Paleozoic trench-arc-basin system in Duobaoshan, Heilongjiang Province. Acta Petrol. Sin. 2022, 38, 2269–2291. [Google Scholar]
  32. Liu, Y.; Cheng, X.Z.; Wang, X.C.; Liu, J.Y.; Wang, L.; Wang, X.L. Copper metal source and enrichment law of the Duobaoshan porphyry copper deposit in Heilongjiang Province. Geol. Sci. 2008, 43, 671–684. [Google Scholar]
  33. Gao, R.; Xue, C.; Lü, X.; Zhao, X.; Yang, Y.; Li, C. Genesis of the Zhengguang Gold Deposit in the Duobaoshan Ore Field, Heilongjiang Province, NE China: Constraints from Geology, Geochronology and S-Pb Isotopic Compositions. Ore Geol. Rev. 2017, 84, 202–217. [Google Scholar] [CrossRef]
  34. Wang, X.C.; Wang, X.L.; Wang, L.; Liu, J.Y.; Xia, B.; Deng, J.; Xu, X.M. Mineralization and later transformation of the Duobaoshan super-large porphyry copper deposit in Heilongjiang Province. Geol. Sci. 2007, 124–133. [Google Scholar] [CrossRef]
  35. Zhou, J.; Han, J.; Zhou, G.; Zhang, X.; Cao, J.; Wang, B.; Pei, S. The Emplacement Time of the Hegenshan Ophiolite: Constraints from the Unconformably Overlying Paleozoic Strata. Tectonophysics 2015, 662, 398–415. [Google Scholar] [CrossRef]
  36. Bai, C.L.; Xie, G.Q.; Zhao, J.K.; Li, W.; Zhu, Q.Q. Discussion on the Metallogenic Characteristics and Deposit Model of the Porphyry Copper and Epithermal Low-Temperature Gold System in the Duobaoshan Ore Field, Eastern Central Asian Orogenic Belt. Earth Sci. Front. 2024, 31, 1081–1103. [Google Scholar]
  37. Li, X.Y.; Cui, J.; Hu, W.S.; Li, C.L. Application of machine learning methods based on multi-source geophysical data in geological body classification: A case study of the Duobaoshan mining area in Heilongjiang. Chin. J. Geophys. 2022, 65, 3634–3649. [Google Scholar]
  38. Yang, Z.Y.; Tang, J.X.; Ren, D.X.; Deng, A.; Wang, Y.; Wu, X. Progress in Geophysical and Geochemical Exploration of the Sinongduo Silver-Polymetallic Deposit in Tibet. Earth Sci. 2024, 49, 1081–1103. [Google Scholar]
  39. Shi, C.Y.; Liang, M.; Feng, B. Background Values of 39 Elements in Stream Sediments in China. Earth Sci. 2016, 41, 234–251. [Google Scholar]
  40. Parsa, M.; Maghsoudi, A.; Carranza, E.J.M.; Yousefi, M. Enhancement and Mapping of Weak Multivariate Stream Sediment Geochemical Anomalies in Ahar Area, NW Iran. Nat. Resour. Res. 2017, 26, 443–455. [Google Scholar] [CrossRef]
  41. Sun, W.; Zheng, Y.; Wang, W.; Feng, X.; Zhu, X.; Zhang, Z.; Hou, H.; Ge, L.; Lv, H. Geochemical Characteristics of Primary Halos and Prospecting Significance of the Qulong Porphyry Copper-Molybdenum Deposit in Tibet. Minerals 2023, 13, 333. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the DE-WDCB Method. This diagram outlines the steps involved in applying the DE-WDCB method, including data collection, trend surface analysis, and anomaly detection.
Figure 1. Flowchart of the DE-WDCB Method. This diagram outlines the steps involved in applying the DE-WDCB method, including data collection, trend surface analysis, and anomaly detection.
Minerals 14 00912 g001
Figure 3. Geological map of Duobaoshan area.
Figure 3. Geological map of Duobaoshan area.
Minerals 14 00912 g003
Figure 4. Frequency distribution histograms of element contents in stream sediments. This figure illustrates the distribution pattern of 14 elements in the study area, highlighting their spatial distribution after natural geochemical processes. (Units for Au, Ag, and Cd are 10−9 (ppb); units for other elements are 10−6 (ppm)).
Figure 4. Frequency distribution histograms of element contents in stream sediments. This figure illustrates the distribution pattern of 14 elements in the study area, highlighting their spatial distribution after natural geochemical processes. (Units for Au, Ag, and Cd are 10−9 (ppb); units for other elements are 10−6 (ppm)).
Minerals 14 00912 g004
Figure 5. Correlation coefficient heat map for key elements.
Figure 5. Correlation coefficient heat map for key elements.
Minerals 14 00912 g005
Figure 6. R-Type clustering of elements indicative of mineralization patterns.
Figure 6. R-Type clustering of elements indicative of mineralization patterns.
Minerals 14 00912 g006
Figure 7. Catchment basin modeling workflow. This figure illustrates the steps involved in constructing the catchment basin landscape model, including sink-filling analysis, flow direction determination, flow analysis, Strahler classification, and final catchment basin delineation.
Figure 7. Catchment basin modeling workflow. This figure illustrates the steps involved in constructing the catchment basin landscape model, including sink-filling analysis, flow direction determination, flow analysis, Strahler classification, and final catchment basin delineation.
Minerals 14 00912 g007
Figure 8. Modeling process of catchment basins (a) DEM; (b) Flow Direction Analysis Map; (c) Flow Accumulation Map; (d) Stream Network Extraction Map; (e) Distribution Map of Sampling Points and Spill Points; (f) Catchment Basin Extraction Map.
Figure 8. Modeling process of catchment basins (a) DEM; (b) Flow Direction Analysis Map; (c) Flow Accumulation Map; (d) Stream Network Extraction Map; (e) Distribution Map of Sampling Points and Spill Points; (f) Catchment Basin Extraction Map.
Minerals 14 00912 g008
Figure 9. Extraction of single element and element combination background values for catchment basins using the Dynamic Trend Surface method.
Figure 9. Extraction of single element and element combination background values for catchment basins using the Dynamic Trend Surface method.
Minerals 14 00912 g009
Figure 10. Visualization of single element and element combination anomalies extracted using the DE-WDCB method (Unit: ppm).
Figure 10. Visualization of single element and element combination anomalies extracted using the DE-WDCB method (Unit: ppm).
Minerals 14 00912 g010
Figure 11. Three-level anomaly mapping using DE-WDCB and IDW methods (Unit: 10−6) (Region A: west of the northern ridge of Xiaoxing’an Mountains; Region B: east of the northern ridge of Xiaoxing’an Mountains).
Figure 11. Three-level anomaly mapping using DE-WDCB and IDW methods (Unit: 10−6) (Region A: west of the northern ridge of Xiaoxing’an Mountains; Region B: east of the northern ridge of Xiaoxing’an Mountains).
Minerals 14 00912 g011
Figure 12. DE-WDCB tertiary anomaly distribution mapping (unit: ppm).
Figure 12. DE-WDCB tertiary anomaly distribution mapping (unit: ppm).
Minerals 14 00912 g012
Figure 13. Local anomaly analysis in the Duobaoshan–Sankuangou area (Region A). (a) Local geological map of Region A; (b) anomaly mapping results using IDW method in Region A (Cu, Fac, Ti); (c) anomaly mapping results using Kriging method in Region A (Cu, Fac, Ti); (df) anomalies of Cu element, element combination, and Ti element extracted using DE-WDCB method.
Figure 13. Local anomaly analysis in the Duobaoshan–Sankuangou area (Region A). (a) Local geological map of Region A; (b) anomaly mapping results using IDW method in Region A (Cu, Fac, Ti); (c) anomaly mapping results using Kriging method in Region A (Cu, Fac, Ti); (df) anomalies of Cu element, element combination, and Ti element extracted using DE-WDCB method.
Minerals 14 00912 g013
Figure 14. Local anomaly analysis in the Sandaowan–Songshugangzi area (Region B). (a) Local geological map of Region B; (b) anomaly mapping results using IDW method in Region B (Cu, Fac, Ti); (c) anomaly mapping results using Kriging method in Region B (Cu, Fac, Ti); (df) anomalies of Cu element, element combination, and Ti element extracted using DE-WDCB method.
Figure 14. Local anomaly analysis in the Sandaowan–Songshugangzi area (Region B). (a) Local geological map of Region B; (b) anomaly mapping results using IDW method in Region B (Cu, Fac, Ti); (c) anomaly mapping results using Kriging method in Region B (Cu, Fac, Ti); (df) anomalies of Cu element, element combination, and Ti element extracted using DE-WDCB method.
Minerals 14 00912 g014
Table 1. Statistics of geochemical parameters of stream sediment elements in the Duobaoshan area.
Table 1. Statistics of geochemical parameters of stream sediment elements in the Duobaoshan area.
ElementIterative
Data
Exclusion
MeanStandard DeviationBackground Value of Greater Khingan RangeRange
Coefficient (%)
Kurtosis
Coefficient
Skewness
Coefficient
Concentration
Coefficient
Variation
Coefficient
Abnormal Outer BandAbnormal Middle BandAbnormal Inner Band
Au41.41.91.0126.50.71.21.41.41.72.22.2
As410.07.57.7104.90.80.91.30.711.813.414.6
Ag14135.9210.573.71093.522.04.81.21.597.7104.1105.7
Cd875.172.080.559.46.1−1.00.91.079.882.284.8
Co314.35.39.034.90.11.21.60.415.215.515.8
Cr055.016.027.540.2−0.60.02.00.359.861.963.0
Cu1222.2131.716.478.40.08.41.45.925.729.931.0
Nb013.72.212.540.7−0.70.01.10.215.115.816.0
Ni224.26.611.577.8−0.40.32.10.328.230.130.6
Pb522.117.322.319.2−1.90.11.00.823.123.523.7
Sb00.40.10.597.3−0.10.50.90.30.50.60.6
Ti44854.33558.22991.324.81.60.11.60.74986.05098.15161.3
W22.03.91.552.73.10.41.31.92.12.22.2
Zn773.583.258.328.60.71.71.31.176.879.581.7
Units: This sentence is not intended to explain the asterisk but rather to clarify the units for the element contents in the table. The asterisk can be removed. for Au, Ag, and Cd are 10−9 (ppb); units for other elements are 10−6 (ppm); background values for elements in the Greater Khingan Range [39].
Table 2. KMO and Bartlett’s test of sphericity.
Table 2. KMO and Bartlett’s test of sphericity.
KMO0.733
Bartlett’s test of sphericityApproximate chi-square3579.201
Degrees of freedom (df)91
Significance0.000
Table 3. The load matrix of the original component.
Table 3. The load matrix of the original component.
ElementComponent Factors
Fac 1Fac 2Fac 3Fac 4
Au−0.3200.6740.3810.153
As−0.1020.884−0.053−0.078
Ag−0.3340.191−0.269−0.540
Cd0.324−0.1100.6950.424
Co0.7850.4510.193−0.236
Cr0.879−0.0100.120−0.218
Cu0.7710.378−0.0470.188
Nb0.404−0.286−0.4390.561
Ni0.888−0.1720.111−0.077
Pb−0.245−0.6960.313−0.099
Sb−0.1810.9380.1010.103
Ti0.804−0.0270.097−0.329
W0.1540.192−0.3380.698
Zn−0.339−0.0810.8010.152
Table 4. Optimal model parameters for control group using IDW and Kriging methods.
Table 4. Optimal model parameters for control group using IDW and Kriging methods.
MethodOptimal ParametersQualitative Parameters
IDWOptimal Distance Index (OOD)1 ≤ α ≤ 2
Search Step Number12
KrigingKernel Function ModelGaussian Model
Search DirectionFour Directions
Search Step Number12
Table 5. Comparison of anomaly areas and mineralization indicators delineated by IDW, Kriging, and DE-WDCB methods.
Table 5. Comparison of anomaly areas and mineralization indicators delineated by IDW, Kriging, and DE-WDCB methods.
MethodsElementsProportion of Anomalies in Total Area (%)Number of Mineral Points within AnomaliesTotal Mineral Point Coverage (%)
High AnomalyMedium AnomalyLow
Anomaly
Total AnomalyA
Region
B
Region
Total
Region
IDWCo0.92 3.03 9.14 13.09 1501542.8
Cr0.95 2.60 14.92 18.48 1501542.8
Cu0.95 2.87 8.68 12.49 1601645.7
Ni0.32 2.65 10.30 13.28 1701748.5
Ti0.72 3.67 7.71 12.10 0225.7
Fac1.00 2.96 8.02 11.99 1701748.5
KrigingCo1.10 3.87 8.88 13.84 1501542.8
Cr1.45 2.46 9.60 13.51 2102160.0
Cu2.15 3.49 8.59 14.23 1711851.4
Ni1.31 3.70 9.75 14.76 1801851.4
Ti1.64 3.65 8.43 13.73 1012.8
Fac0.08 3.13 10.48 13.69 2002057.1
DE-WDCBCo0.81 2.92 26.48 30.21 1721954.3
Cr1.09 2.24 27.47 30.81 1721954.3
Cu1.05 2.16 23.29 26.50 1822057.1
Ni0.53 2.21 23.09 25.84 1822057.1
Ti1.16 2.68 22.88 26.72 0112.9
Fac1.36 2.13 23.99 27.48 1962571.4
Table 6. Exception level 3 threshold table for IDW and Kriging methods (Unit: ppm).
Table 6. Exception level 3 threshold table for IDW and Kriging methods (Unit: ppm).
Thresholds for Abnormal Grading (ppm)Low ValueMedium ValueHigh Value
Co15.215.515.7
Cr59.861.863.0
Cu25.729.831.0
Ni28.130.130.5
Ti4985.95098.05161.2
Fac1.11.61.9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cui, Z.; Chen, J.; Zhu, R.; Zhang, Q.; Zhou, G.; Jia, Z.; Liu, C. Dynamic Enhanced Weighted Drainage Catchment Basin Method for Extracting Geochemical Anomalies. Minerals 2024, 14, 912. https://doi.org/10.3390/min14090912

AMA Style

Cui Z, Chen J, Zhu R, Zhang Q, Zhou G, Jia Z, Liu C. Dynamic Enhanced Weighted Drainage Catchment Basin Method for Extracting Geochemical Anomalies. Minerals. 2024; 14(9):912. https://doi.org/10.3390/min14090912

Chicago/Turabian Style

Cui, Zijia, Jianping Chen, Renwei Zhu, Quanping Zhang, Guanyun Zhou, Zhen Jia, and Chang Liu. 2024. "Dynamic Enhanced Weighted Drainage Catchment Basin Method for Extracting Geochemical Anomalies" Minerals 14, no. 9: 912. https://doi.org/10.3390/min14090912

APA Style

Cui, Z., Chen, J., Zhu, R., Zhang, Q., Zhou, G., Jia, Z., & Liu, C. (2024). Dynamic Enhanced Weighted Drainage Catchment Basin Method for Extracting Geochemical Anomalies. Minerals, 14(9), 912. https://doi.org/10.3390/min14090912

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop