Comparison of Methods for Determining the Thresholds of Geochemical Anomalies and the Prospecting Direction—A Case of Gold Deposits in the Gouli Exploration Area, Qinghai Province

Geochemical anomalies play an important role in mineral exploration, because they are closely related to the formation of ore deposits. In this paper, the Au in the geochemical samples of 1:50,000 stream sediments in the Gouli exploration area, Qinghai Province, is selected as an example. The cumulative frequency method, local singularity analysis, and exploration data analysis method are used to process the Au data and determine the abnormal thresholds. Robust principal component analysis is used to explore the elemental associations. Among them, the effect of local singularity analysis is the best and it delineates the weak geochemical anomalies not delineated by cumulative frequency method and exploration data analysis method, which shows its superiority. Robust principal component analysis shows that Sb, W, and As are closely related to Au mineralization. Three anomalous areas are delineated in the Gouli exploration area. Gold ore-body is found to occur within anomaly area I. The other anomalous areas are highly consistent with abnormal elemental composition of Sb, W, and As. Moreover, these two areas are extremely conformable with the high-value areas determined by the binary state method. The optimality of the three methods are compared comprehensively, it is pointed out that it is more reasonable to use the local singularity analysis to determine the abnormal threshold. The anomaly areas II and III of Au determined by local singularity analysis have shown great potential for prospecting.


Introduction
At the International Conference in Sydney in 1976, the abnormal threshold was defined as a value determined by the geochemical makers according to some analysis and test results of samples. It can identify the anomalies associated with the mineralization. A general definition of abnormal threshold was put forward as the value that can distinguish abnormal and non-abnormal data. The determination of the threshold of geochemical anomalies is directly related to the identification of mineralization information. Therefore, the determination of the threshold of geochemical anomalies is very significant for mineral exploration.
Since the 1970s, many countries in the world have carried out regional geochemical surveys successively and obtained primary and secondary samples, such as rocks, stream sediments, water, and

Geological Settings
The Gouli exploration area is covered by grassland and is located in the South East part of the Eastern Kunlun tectonic belt (Eastern Kunlun Central Orogenic Belt) and between latitudes 35 • Figure 1C). Since the Phanerozoic eons, there have been two major tectonic movements in the region, namely, the early Paleozoic Caledonian movement (the original Tethys evolution), Late Paleozoic, and the early Mesozoic Paleotethys evolution.
The Gouli area is located between the Kunzhong fault and Kunbei fault, which is mainly restricted by the Kunzhong and Kunbei fault zones ( Figure 1A,B). The structural deformation is intense and fault structures have developed. The structures in the mining area are mainly NE-extending faults of different grades and sequences, followed by NW-extending and nearly EW-extending faults with genetic connection, and a few SN-extending and NNW-extending faults ( Figure 1C). The NE-extending faults are mainly compressive or compressive-torsional faults, which constitute the main faults and are characterized by episodic activities. These faults control the distribution direction of strata and magmatic rocks in the area. Results of LA-ICP-MS (Laser Ablation Inductively Coupled Plasma Mass Spectrometry) zircon U-Pb isotope dating of the study area conducted at Wuhan Sample Solution Analytical Technology Co. Ltd., gives 230 ± 1.2 Ma for monzonite granite, and 210 ± 2 Ma for syenite granite. The formation ages of the rock masses are basically the same as that of the gold deposits in Eastern Kunlun (250-205 Ma) [10,[19][20][21]. Hence, it is presumed that the fault structure may be genetically related to gold and polymetallic mineralization.   In addition to the Holocene alluvial deposits and residual slope deposits, the strata exposed in the mining area are mainly schist, gneiss, amphibolite, and marble of the Jinshuikou Group. Magmatic rocks are well developed in the mining area. The intrusive rocks are mainly syenite granite and monzonite granite, as well as some granodiorite and plagioclase granite. Ultrabasic dikes, diorite porphyrites dikes, and quartz veins are occasionally exposed. Most of the rock masses and dikes in this area are developed in the gneiss of Jinshuikou Group, and have intrusive contact and fault contact with the gneiss ( Figure 1C).
The Yelonggou and the Gazhima gold deposits are both located in the strata of the Jinshuikou Group and in a fault structural system ( Figure 1C). According to mining data in the area, the contents of Au in all rock types of Jinshuikou Group are 2-3 times higher than the average crustal value, and the abundance of Au in granite is higher, which provides fluid and some material sources for gold mineralization [22,23].

Data Sources and Eigenvalues Statistics
The 764 samples collected in the regional geochemical survey were 1:50,000 stream sediment geochemical analysis samples from an area of 125 km 2 ( Figure 2). Au, Ag, Cu, Cr, As, Sb, Ni, Zn, and other elements were analyzed in Hebei Institute of regional geology and Mineral Resources Research Laboratory (Hebei, China). The concentrations of Au were tested by graphite furnace-atomic absorption spectrometry (GF-AAS) [24,25]. The concentrations of Ag were tested by emission spectrometry (ES). The concentrations of Cu, Cr, Ni, and Zn were tested by X-ray fluorescence (XRF) [26]. The concentrations of As and Sb were tested by hydride generation-atomic fluorescence spectrometry (HG-AFS). The Yelonggou and the Gazhima gold deposits are both located in the strata of the Jinshuikou Group and in a fault structural system ( Figure 1C). According to mining data in the area, the contents of Au in all rock types of Jinshuikou Group are 2-3 times higher than the average crustal value, and the abundance of Au in granite is higher, which provides fluid and some material sources for gold mineralization [22,23].

Data Sources and Eigenvalues Statistics
The 764 samples collected in the regional geochemical survey were 1:50,000 stream sediment geochemical analysis samples from an area of 125 km 2 ( Figure 2). Au, Ag, Cu, Cr, As, Sb, Ni, Zn, and other elements were analyzed in Hebei Institute of regional geology and Mineral Resources Research Laboratory (Hebei, China). The concentrations of Au were tested by graphite furnace-atomic absorption spectrometry (GF-AAS) [24,25]. The concentrations of Ag were tested by emission spectrometry (ES). The concentrations of Cu, Cr, Ni, and Zn were tested by X-ray fluorescence (XRF) [26]. The concentrations of As and Sb were tested by hydride generation-atomic fluorescence spectrometry (HG-AFS). The estimated range of the concentrations of Au in the study area is shown in Figure 3. The majority of the samples yielded low concentration of Au, while smaller portion are characterized by relatively high concentration. In general, the maximum concentration of Au recorded in the study area is 8.5 ppb, and the minimum concentration is 0.5 ppb.
According to the above geological characteristics of the mining area, the eigenvalues of stratigraphic units and intrusive rocks of different ages in the study area are statistically analyzed The estimated range of the concentrations of Au in the study area is shown in Figure 3. The majority of the samples yielded low concentration of Au, while smaller portion are characterized by relatively high concentration. In general, the maximum concentration of Au recorded in the study area is 8.5 ppb, and the minimum concentration is 0.5 ppb. 2019, 9, x FOR PEER REVIEW Figure 3. The histogram of concentration of Au in the study area.  According to the above geological characteristics of the mining area, the eigenvalues of stratigraphic units and intrusive rocks of different ages in the study area are statistically analyzed respectively (Table 1). We know that the closer the skewness is to 0, the closer the kurtosis is to 3, then the closer the distribution of the data is to normal distribution. It can be seen from Table 1 that the distribution of Au in the study area does not conform to normal distribution. The characteristic distributions of the whole area are closer to those of Triassic intrusive rocks and Jinshuikou Group strata. Therefore, it is reasonable to use the uniform abnormal threshold in the study area.

Methods
This paper mainly presents statistics on the spatial distribution of Au anomalies in the Gouli exploration area. We know that the spatial distribution of Au must be influenced by strata, structures, magmatic rocks, and other ore-controlling factors. From the above, we can see that the enrichment of Au in strata of the Jinshuikou Group is beneficial to the enrichment of Au. The role of structures and magmatic rocks are also significant to Au enrichment. In other words, the more developed the strata of the Jinshuikou Group, the structures, and the magmatic rocks are, the more enriched the distribution of Au is. Au enrichment must be accompanied by associated elements of enrichment. In this paper, the cumulative frequency method, local singularity analysis, and exploration data analysis method are used to analyze the spatial distribution of Au, and robust principal component analysis is used to analyze the Au and its associated elements.

Cumulative Frequency Method
The method of cumulative frequency sorts a group of data in order from small to large, and takes the part of the data at the top of the element as the abnormal part to research. The determination of the top data ratio is based on the geological characteristics of the study area and the results of previous comprehensive research on the characteristics of the study area, such as geophysics, geochemistry, remote sensing, and other characteristics of the study area; the top data ratio is generally 85% [27][28][29]. Geochemical exploration practice shows that the concentrations of geochemical elements in the crust or in a geological body are randomly distributed. To accurately describe this set of random variables, according to the relevant theories of probability theory and mathematical statistics, the probability distribution function can be obtained; that is, the relationship between the element concentration and its corresponding probability is: In Equation (1), f (x) is the probability density function, and p(c ≤ x) is the probability that the concentration is less than or equal to a certain value x. One of the premises of using the statistical integral method is that F(x) exists for any element of any geological body, that is, f (x) is integrable for any element of any geological body. Because the distribution of geochemical elements is random, the expression of f (c) is determined by the specific distribution of elements.
If the distribution of geochemical data obeys normal distribution, its probability density function is: Therefore, the probability distribution function of this element can be expressed as: In Equations (2) and (3), µ is the expected value of the element (i.e., the average) and σ is the standard deviation of the element. The element concentration corresponding to the cumulative frequency of the element can be obtained by integral calculus with Equation (3), that is, the abnormal threshold value of the element is determined by the cumulative frequency method. The equations for calculating µ and σ are: In Equations (4) and (5), n represents the number of elements, and x represents the average number of elements.
It is convenient to use probability theory and mathematical statistics for geochemical data with normal distribution, but it is not suitable for geochemical data with general distribution. This is mainly because the integral function form is not well defined, and the integration process is also very complicated. According to Table 1, the distribution of Au in this study area does not obey normal distribution, so it is not suitable for integral calculus. In this paper, we will use the C program to determine the abnormal threshold of the Au; the code can be found in the Appendix A.

Local Singularity Analysis
The essences of fractal geometry are self-similarity and scale-invariance. Self-similarity refers to the statistical similarity between parts and the whole in shape, function, information, structure, and other aspects; Scale-invariance refers to enlarging or reducing the size appropriately without altering the structure of the object [30]. In the field of geology, fractal application has been very extensive, for example, Sanderson et al. [31] studied the fractal relationship between the thickness and grade of the La Codosera gold deposit vein in Spain. In china, Ding et al. [32] studied the fractal distribution of gold deposits, while Shi et al. [33] studied the continuous multi-fractal distribution model of geochemical fields. In this study, the authors will take the local singularity analysis in a multi-fractal analysis as an example.
The local singularity analysis is an important progress for fractal or multi-fractal modeling of geochemical data [4]. It is defined as the characterization of the anomalous behaviors of singular physical processes that often result in anomalous amounts of energy release or material accumulation within a narrow spatio-temporal interval [1,[34][35][36]. The singularity can be estimated from observed element concentration within small neighborhoods based on the following equation: where X represents element concentrations, c is a constant value, α is the singularity, r is a normalized distance measure, such as block cell edge, and E is the Euclidian dimension [37]. The window-based method can be used to estimate the singularity index α [1,34-36] based on either raw point geochemical data or raster maps using the following steps: (1) define a set of sliding squares (or other shapes) windows A(r) with variable window sizes, r min = < r 1 < r 2 < · · · < r n = r max , and calculate the average concentration C[A(r i )] for each window size for a given location on the map, where C[A(r i )] is equal to the sum of all the cell concentrations divided by the total number of cells within the window; (2) plot the pair data set of C[A(r i )] (i = 1, 2, · · · , n), and r i in a log-log graph, and obtain the following linear relationship: The value of α − 2 can be estimated from the slope of the fitted straight line; in step (3), perform a similar procedure to all other locations on the geochemical map. For a geochemical map, most of the areas are linked with a singularity value α close to 2, representing a normal background, whereas the areas with α > 2 or α < 2 represent a depleted pattern and enrichment pattern of element concentrations, respectively [34]. In other words, if α > 2, it means that the element is depleted, if α < 2, it means that the element is highly enriched. The value of singularity index reflects the degree of local mineralization singularity and mineralization enrichment intensity. In the case of α < 2, when α is infinitely close to 2, the concentration at the same window as α is the abnormal threshold of the element.

Exploration Data Analysis Method
The exploration data analysis method [38] uses boxplot to interpret the distribution of data and determine the abnormal threshold [39]; this method does not require that the data (raw data, after processing of positive cosine, inverse cosine, or logarithm) obey normal distribution. As shown in Figure 4, five parameters (maximum, upper quartile, median, lower quartile, and minimum) are defined in the boxplot. These parameters can describe the important characteristics, structure, and distribution of the data [40]. The detailed method sorts the data from large to small in order to determine the maximum, minimum, lower quartile (Q 1 ), median (Q 2 ), and upper quartile (Q 3 ). The internal divergence (hs), the upper limit of the abnormal points (F l ), and the threshold of the abnormal points (F u ) can be determined with the following Equations (8)- (10).
All data less than F l or greater than F u are outliers, in which above F u are called upper outliers (equivalent to positive anomaly), below F l are lower outliers (equivalent to negative anomaly); it can be seen from Equations (9) and (10)

Robust Principal Component Analysis
Robust Principal Component Analysis (RPCA) is a very effective method for identifying geochemical element assemblages. However, as typical compositional data, the closure effect of geochemical data may cause virtual correlation of the original data in the process of principal component analysis [41][42][43][44]. In order to weaken the influence of the closure effect, geochemical data need to be "opened" in pretreatment. At present, the main methods of "opening" data are the additive log-ratio (alr) [42], centered log-ratio (clr) [42], and isometric log-ratio (ilr) [45,46]. Carranza [46] compares in detail the application effects of the above three methods in the delineation of geochemical anomalies in stream sediments. Because the denominator selection of the alr-transformation is artificial, the sum of all variables after clr-transformation is zero [46][47][48] and the resulting data are collinear, so alr-transformation and clr-transformation are not suitable for robust statistics but ilr-transformation is. The ilr-transformation solves the problem of data collinearity resulting from the clr-transformation, while preserving all its advantages, such as isometry between the simplex and real space [45]. Even so, the dimensions of data are reduced by one dimension after ilr-transformation, so it is difficult to interpret the results of ilr-transformation. The problem of data interpretation is solved by using ilr-transformation to transform the compositional data into the European space, combining with the loadings and scores, and then transforming the loadings and scores into the clr space to get the connection with the original data [49]. The results of RPCA are presented in a compositional biplot to present loadings (representing information of variables) and scores (representing information of observations). The scores of RPCA are interpreted to reflect the presence of mineralization.

The Results Processed by Cumulative Frequency Method, Local Singularity Analysis and Exploration Data Analysis Method
The three methods mentioned above show that the abnormal threshold of Au determined by local singularity analysis is the smallest, followed by cumulative frequency method, as shown in Table 2. From Table 2, we can see that the gold deposits of the Gouli exploration area in Qinghai Province mainly occur in the Triassic intrusive rocks and Jinshuikou Group strata. According to the abnormal thresholds of Au determined by the above three methods, abnormal isoline graphs are drawn by Surfer software and superimposed on mineral geological maps, as shown in Figure 5.

The Results Processed by RPCA
The results of RPCA for clr-transformed regional geochemical data were used to explore the elemental associations. The resulting loadings were visualized graphically in a biplot ( Figure 6). The first group from the loadings of PC1 consists of Zn, Cu, Co, Cr, Ba, Pb, Fe, Hg, Mo, Mn, Bi, Sn, La, Ag, and Ni, and the second group from the loadings of PC2 consists of As, Sb, Au, and W. As can be seen from Figure 6, the angle between Au and Sb is the smallest, followed by W and As, indicating that Au and Sb, W, and As have strong correlation.
RPCA is used to explore the elemental associations. Sb, W, and As are closely related to Au mineralization; in other words, Au, Sb, W, and As may have the same behavior in magmatic fluids. So Sb, W, and As are good geochemical indicators for tracing the Au mineralization in the Gouli exploration area, Qinghai Province. To identify the anomalies related to Au, Sb, W, and As, the singularity indices and the concentrations of Au, Sb, W, and As in the same window as the singularity index α are calculated by Geodas and Arcgis software (Table 3).
Note: The dimension of abnormal threshold value is ppb, the dimension of abnormal area is km 2 , and the rest is dimensionless.
From Table 2, we can see that the gold deposits of the Gouli exploration area in Qinghai Province mainly occur in the Triassic intrusive rocks and Jinshuikou Group strata. According to the abnormal thresholds of Au determined by the above three methods, abnormal isoline graphs are drawn by Surfer software and superimposed on mineral geological maps, as shown in Figure 5.

The Results Processed by RPCA
The results of RPCA for clr-transformed regional geochemical data were used to explore the elemental associations. The resulting loadings were visualized graphically in a biplot ( Figure 6). The first group from the loadings of PC1 consists of Zn, Cu, Co, Cr, Ba, Pb, Fe, Hg, Mo, Mn, Bi, Sn, La, Ag, and Ni, and the second group from the loadings of PC2 consists of As, Sb, Au, and W. As can be seen from Figure 6, the angle between Au and Sb is the smallest, followed by W and As, indicating that Au and Sb, W, and As have strong correlation. RPCA is used to explore the elemental associations. Sb, W, and As are closely related to Au mineralization; in other words, Au, Sb, W, and As may have the same behavior in magmatic fluids. So Sb, W, and As are good geochemical indicators for tracing the Au mineralization in the Gouli exploration area, Qinghai Province. To identify the anomalies related to Au, Sb, W, and As, the singularity indices and the concentrations of Au, Sb, W, and As in the same window as the singularity index α are calculated by Geodas and Arcgis software (Table 3).

Relationship between Fracture and Mineralization
It can be seen that the distribution of Au anomalies is closely related to faults when the spatial distribution of Au anomalies are superimposed on the geological map. Distribution of the Au anomaly areas in Figure 5b-d are beaded in the NE direction, consistent with the NE-extending main faults; the anomaly area of Yelonggou is consistent with the faults developed in the Middle Triassic monzonite granite, and the extension direction of the anomaly area of Gazhima is also consistent with the nearly EW-extending faults. The distribution directions of other minor anomaly areas are roughly the same as the extending of sub-faults. It can be concluded that the fault structures not only control the distribution direction of gold deposits, but also provide channels for the migration and enrichment of Au. The relationship between fault structure and gold deposits in the study area is indeed very close.
According to the data of the study area, during the Late Variscan-Indosinian orogenic period, successive activities (with the characteristics of strike-slip fault) took place along the fault zone ( Figure 1A,B), accompanied by massive intrusion of rock-magma ( Figure 1C), which made metallogenic fluids migrate along the fault zone system in large amounts and extract gold and other metallogenic materials of wall rocks along the way. When the gold-bearing fluid entered the ore-controlling and ore-hosting structure, the gold-bearing hydrothermal fluid and wall rocks underwent metasomatic alteration due to the regional thermodynamic action, which eventually led to the formation of the Yelonggou gold deposit. The Deep fault in the Kunzhong region controlled the regional mineralization ( Figure 1B), the Kunzhong fault is the northern slope structural belt in the Caledonian period, while the Wahongshan-Hot Spring fault zone ( Figure 1A) underwent multiple transformations and superimposed transformation of deformation at different tectonic levels later, accompanied by massive intrusion of rock-magma, which made metallogenic fluids migrate along the fault zone system in large amounts and extract gold and other metallogenic materials of wall rocks along the way. The ore-bearing fluid was transported to the ore-hosting site and precipitated to form the Gazhima gold deposit. Therefore, it can be concluded that the faults in the study area are mainly ore-hosting structures.

The Choice of the Optimal Solution Method
As can be seen from Figure 5b-d, the anomaly areas are mainly distributed in the Middle Triassic monzonite granite, Late Triassic syenite granite, and Jinshuikou Group strata. However, the anomaly areas delineated by the three methods are distinctly different in many aspects.
The abnormal threshold calculated by cumulative frequency method is 1.9 ppb. The abnormal thresholds of elements delineated by accumulative frequency method have certain artificiality and empiricism. This mainly relies on the content of sampling elements to be sorted, which is more practical in the case of simple geological conditions, but less practical in areas with very complex geological conditions. The area delineated by exploration data analysis method is the largest (the abnormal threshold is 2.55 ppb), and many abnormal areas have been left out. The anomaly shapes are scattered and not continuous, and the shapes are small-scale block distributions. The concentrations of most elements in the study area are regarded as background values. This is because exploration data analysis method does not require data to obey normal distribution or lognormal distribution from the structure of data itself. It can resist the influence of ultra-high and ultra-low values. However, when geochemical sampling is carried out in mining areas or in areas with high background values, such as enrichment areas of individual elements, the range of data varies greatly, and there are more middle and high values, which makes the quartile spacing smaller and causes more missing anomalies.
As we can see, the anomaly areas calculated by the above two methods mainly develop along NE-extending faults, especially at the intersection of faults (Figure 5b,d). However, the above two methods only consider the frequency distribution and not the spatial variations of geochemical data from low background patterns. In other words, they do not consider the compositional nature of geochemical data [47,48,[50][51][52][53][54][55].
The formation of ore deposits is the result of the coupling of various geological ore-controlling factors (such as strata, structure, rock mass, alteration, etc.) and various metallogenic processes (physics, chemistry, and biology) in the specific geological tectonic environment [56]. The ore deposits are distributed in time and space in multi-stage zonation, and are regarded as geological anomalies. Local singularity analysis emphasizes that mineralization and the formation of ore deposits are the result of geological anomalies and other ore-controlling factors [34], which are not the characteristics of cumulative frequency method and exploration data analysis method. Moreover, the local singularity analysis delineates the weak anomaly areas in the Gouli exploration area (Figure 5c); in this respect, the local singularity analysis is not only a powerful multifractal tool to identify weak anomalies, but also a kind of local neighborhood statistical analysis that can reduce the effects of regional background and provide useful statistical information by involving the data within a small singularity around a specific spatial location [1,4,34,37,49,[57][58][59][60][61][62][63][64][65][66]. The coverage is still relatively thick in Qinghai Province with high altitude, which means that the geochemical anomalies in stream sediments are very weak and easy to be ignored. With the discovery of a large number of ore deposits, concealed ore-bodies have become the trend for prospecting in the future. For deep-buried concealed ore-bodies, the geochemical anomalies exposed on the surface are very weak. Comparatively speaking, the cumulative frequency method and the exploration data analysis method have better effects for areas with obvious anomalies, but for areas with weak anomalies, the effects of the two methods are not as significant as those of local singularity analysis. Therefore, the authors believe that it is optimal to use local singularity analysis to determine the abnormal threshold of Au in the Gouli exploration area, Qinghai Province.

Mineralization-Related Associated Elements
The spatial distributions of singularity indexes of Au, Sb, W, and As (Figure 7a-d) are generally consistent and mainly spread in the NE-extending fault, according to the main fault direction in the study area, suggesting that the mineralization of Au is similar with that of Sb, W, or As. A singularity index α < 2 represents elements that are enriched, while a singularity index α > 2 represents elements that are depleted. The areas of enriched elements in the study area are small and dispersed (Figure 7a-d). In the case of spatial scattered distribution, the Yelonggou gold deposit and Gazhima gold deposit occur in or near to the anomalous areas with low singularity index, suggesting that the local singularity analysis can well identify the geochemical anomalies. The spatial distributions of Sb, W, and As, and the compositional elements of Au, Sb, W, and As contour maps obtained by local singularity analysis are shown in Figure 8. Positions of high concentrations of Sb, W, and As located in Yelonggou area, and Gazhima area results were the same (Figure 8a-c), especially at the intersection of NNE-extending and NE-extending faults. The values of compositional elements of Au, Sb, W, and As (Figure 8d) are calculated based on the scores from the ilr-transformation, which is a comprehensive reflection of Au mineralization in the Gouli exploration area, Qinghai Province. The spatial distribution of the compositional elements not only highlights the known areas of Au mineralization, but also identifies new anomaly areas, which means the analysis by the compositional elements is more reasonable than those by single elements. Meanwhile, the abnormal spatial distribution of the compositional elements is highly consistent with that of Au. Therefore, the compositional elements of Au, Sb, W, and As can be good geochemical comprehensive indicators for tracing the Au mineralization. The spatial distributions of Sb, W, and As, and the compositional elements of Au, Sb, W, and As contour maps obtained by local singularity analysis are shown in Figure 8. Positions of high concentrations of Sb, W, and As located in Yelonggou area, and Gazhima area results were the same (Figure 8a-c), especially at the intersection of NNE-extending and NE-extending faults. The values of compositional elements of Au, Sb, W, and As (Figure 8d) are calculated based on the scores from the ilr-transformation, which is a comprehensive reflection of Au mineralization in the Gouli exploration area, Qinghai Province. The spatial distribution of the compositional elements not only highlights the known areas of Au mineralization, but also identifies new anomaly areas, which means the analysis by the compositional elements is more reasonable than those by single elements. Meanwhile, the abnormal spatial distribution of the compositional elements is highly consistent with that of Au. Therefore, the compositional elements of Au, Sb, W, and As can be good geochemical comprehensive indicators for tracing the Au mineralization.

Delineation of Prospecting Targets and Anomaly Verification
At present, three gold ore-bodies have been delineated in Yelonggou area, including one industrial ore-body, one low-grade ore-body, and one gold mineralized body. The ore-bodies were mostly veins, stratiform, and lenticular. The widths of gold ore-bodies are generally 0.21-2.57 m. The ore-bodies were controlled by single engineering, so it was speculated that the length of the longest was 80 m and the maximum extension was 133 m. The Au average grade ranged generally from 1.21 ppm to 8.85 ppm. The ore-bodies were strictly controlled by the mineralized alteration zone and produced in the fault gouges in the structural alteration zone. The ore minerals in gold ore-bodies are mainly pyrite, limonite, and a little arsenopyrite, while gangue minerals were mainly quartz, followed by feldspar, sericite, calcite, chlorite, epidote, amphibole, pyroxene, and graphite.
One copper-gold ore-body, one silver-gold ore-body, one gold ore-body, and three mineralized bodies have been delineated in Gazhima area [67]. The No. 1 copper-gold ore-body was lenticular and stratiform, controlled by trench 01, trench 02, and trench 11. The length of the ore-body is about 100 m and the thickness is 13.2 m, and the Au average grade is 0.73%. The ore-body occurred in quartzite of altered fracture zone FS1, and it could be seen that altered mineralizations were pyritization, ferritization, silicification, and partial malachitization. The wall rocks of the upper and lower wall were weakly altered biotite schist. The No. 2 gold ore-body is lenticular and stratiform, controlled by trench 04 and trench 12. The length of the ore-body was not determined and the thickness is about 1m, with an average Au grade of 1.2 ppm. The ore-body occurred in quartzite of altered fracture zone FS1. It could be seen that altered mineralizations were pyritization, ferritization, kaolinization, silicification, and partial malachitization. The wall rocks of the upper and lower wall were mica schist. The No. 3 ore-body was a silver-gold ore-body and the main ore-body in the working area. It was found in FS8 quartzite of altered zone, controlled by trench 21 and trench 32. The length of the ore-body was about 250 m and the thickness was 3 m, and the Au average grade was 3.04 ppm. The ore-body occurred in quartzite of the altered fractured zone. There were many impurities in the quartzite, which appears gray-white and smoky-gray. The main

Delineation of Prospecting Targets and Anomaly Verification
At present, three gold ore-bodies have been delineated in Yelonggou area, including one industrial ore-body, one low-grade ore-body, and one gold mineralized body. The ore-bodies were mostly veins, stratiform, and lenticular. The widths of gold ore-bodies are generally 0.21-2.57 m. The ore-bodies were controlled by single engineering, so it was speculated that the length of the longest was 80 m and the maximum extension was 133 m. The Au average grade ranged generally from 1.21 ppm to 8.85 ppm. The ore-bodies were strictly controlled by the mineralized alteration zone and produced in the fault gouges in the structural alteration zone. The ore minerals in gold ore-bodies are mainly pyrite, limonite, and a little arsenopyrite, while gangue minerals were mainly quartz, followed by feldspar, sericite, calcite, chlorite, epidote, amphibole, pyroxene, and graphite.
One copper-gold ore-body, one silver-gold ore-body, one gold ore-body, and three mineralized bodies have been delineated in Gazhima area [67]. The No. 1 copper-gold ore-body was lenticular and stratiform, controlled by trench 01, trench 02, and trench 11. The length of the ore-body is about 100 m and the thickness is 13.2 m, and the Au average grade is 0.73%. The ore-body occurred in quartzite of altered fracture zone FS1, and it could be seen that altered mineralizations were pyritization, ferritization, silicification, and partial malachitization. The wall rocks of the upper and lower wall were weakly altered biotite schist. The No. 2 gold ore-body is lenticular and stratiform, controlled by trench 04 and trench 12. The length of the ore-body was not determined and the thickness is about 1 m, with an average Au grade of 1.2 ppm. The ore-body occurred in quartzite of altered fracture zone FS1. It could be seen that altered mineralizations were pyritization, ferritization, kaolinization, silicification, and partial malachitization. The wall rocks of the upper and lower wall were mica schist. The No. 3 ore-body was a silver-gold ore-body and the main ore-body in the working area. It was found in FS8 quartzite of altered zone, controlled by trench 21 and trench 32. The length of the ore-body was about 250 m and the thickness was 3 m, and the Au average grade was 3.04 ppm. The ore-body occurred in quartzite of the altered fractured zone. There were many impurities in the quartzite, which appears gray-white and smoky-gray. The main mineralizations were galena mineralization and ferritization, while the wall rocks were relatively fragmented and strongly altered, and the main wall rock alterations were chloritization, kaolinization, and sericitization. The No. 4 gold mineralized body was lenticular and stratiform, controlled by trench 01, trench 02, trench 03, trench 04, trench 11, and trench 12. The length of the mineralized body was about 350 m and the thickness was 1-3 m, and the Au average grade was 0.24-0.9 ppm. The mineralized body occurred in the altered fracture zone. It could be seen that altered mineralizations were pyritization, ferritization, kaolinization, silicification, and partial malachitization. The wall rocks of the upper and lower wall were weakly altered biotite schist. The No. 5 silver-gold mineralized body was lenticular, controlled by trench 21, trench 23, trench 32, trench 24, trench 26, and trench 19. The length of the mineralized body was about 1600 m and the thickness was 1-15 m. The main mineralizations were pyritization and galena mineralization. The No. 6 gold mineralized body was lenticular and stratiform, controlled by trench 08 and trench 09. The length of the ore-body was unmeasured and the thickness was 1-4.5 m, and the Au average grade was 0.13-0.18 ppm. The mineralized body occurred in the altered fractured zone. It could be seen that altered mineralizations were pyritization, ferritization, kaolinization, and silicification. The wall rocks of the upper and lower walls were mica schist.
According to geological maps for the field and the Yelonggou and Gazhima gold deposits, we can conclude the output characteristics of gold deposits in the Gouli exploration area, Qinghai Province: (1) There are faults or structure alteration zones; (2) the geological unit is Triassic intrusive rock or Jinshuikou Group strata; (3) mineralization alteration is mainly pyritization, ferritization, and galena mineralization; (4) there are strong anomalies for Au; (5) there are the abnormal compositional elements of Sb, As, and W. According to the binary state thinking, the following scoring rules are stipulated (Table 4). The higher the cumulative sum of these scores, the higher the prospecting potential. In this paper, the whole area is divided into 100 small areas based on the density of 1.25 km × 1 km. According to the binary state scoring rules, the scores are shown in Figure 9.
The highest score is 5, followed by 4 ( Figure 9). The areas according to the two scores may have great prospecting potentials. We can see that Yelonggou and Gazhima gold deposits are both in the highest score areas. Compared with the high-score areas in Figure 9, the anomaly areas I, II, and III can be taken as prospecting targets for gold ( Figure 10) in the Gouli exploration area, which may have great prospecting potential for gold. about 20 m, silicification and ferritization are strong, and soil anomalies show that gold and silver are abnormally high, with a gold grade of 1.34 ppm. Anomaly areas I, II, and III are in a NE-extending fault, which indicates that the metallogenic conditions of anomaly areas II and III are similar to that of anomaly area I. Therefore, anomaly areas II and III have high prospecting potential for gold.   Trench 01 (TC01) (Figure 11) was used to verify the anomaly area I (Figure 10) in the Gouli exploration area. The Au mineralization alteration zone (Figure 11a) hosting one gold ore-body was found successively; the gold ore-body is strictly controlled by the NE-extending fault. There are honeycomb-like erosion holes in the ore-body, the average grade of gold ore-body is 63 ppm, while maximum grade is 202 ppm. Mineralization alteration is very intense and the width of mineralization alteration zone is about 13 m. Mineralization characteristics show strong ferritization (Figure 11b,f), malachitization (Figure 11c,d), pyritization in the quartz vein (Figure 11e), and silicification. In the northern part of the gold ore-body, the width of the fractured alteration zone is about 20 m, silicification and ferritization are strong, and soil anomalies show that gold and silver are abnormally high, with a gold grade of 1.34 ppm.   Anomaly areas I, II, and III are in a NE-extending fault, which indicates that the metallogenic conditions of anomaly areas II and III are similar to that of anomaly area I. Therefore, anomaly areas II and III have high prospecting potential for gold.

Conclusions
Based on the detailed comparison of the geochemical anomaly processing methods of Au and the study of the anomaly threshold of Au and its associated elements in the Gouli exploration area, Qinghai Province, the following conclusions are drawn in this paper: (1) Comparing with the cumulative frequency method and the exploration data analysis method, the authors consider that local singularity analysis is the optimum method to be applied in the study area. The local singularity analysis not only highlights the strong Au anomaly areas, but also delineates the weak gold anomaly areas. Under complex geological backgrounds, it is a powerful tool for identifying weak anomalies.
(2) Anomaly areas II and III have similar anomaly intensity to anomaly area I and are in the same tectonic system. Moreover, anomaly areas II and III are located in the anomaly areas of compositional elements of Sb, W, and As, which are most closely related to Au. Therefore, anomaly areas II and t=float(i+1); num=float(n); if(t/num>=0.85) {//Determine the location of elements at a cumulative frequency of 85%.