The Bardaskan area is in the geological map on scale 1:100,000 is one of the rectangular sheets of Kashmar map which is on scale 1:250,000. The study area is in geographical coordinates 57˚00' - 57˚30' eastern longitude and 35˚00' - 35˚30' northern latitude. The range of the Bardaskan sheet is among the cities of Khorasan Razavi province of Iran. Major geological subdivisions of Iran and geology map of Bardaskan area are shown in Figure 1. Khorasan province has a good position regarding the capacity and diversity of mineral reserves due to its special location, which is the location of the major geological phenomena. The vast reserves of oil and gas in the north of the province are abundant. Also, deposits of coal, copper, iron, manganese, fire clay, magnesite, white clay, various types of building stones such as marble, travertine, iron ore and gypsum are found in this province. The most important minerals in Khorasan Razavi province that have economic mines include:

- Magnesite.

- Fluorite.

- Bauxite.

- Chromite.

- Copper.

- Turquoise.

- Coal.

- Iron.

Considering the existence of different economic mines in Khorasan Razavi province, the Bardaskan region, which is one of the most susceptible areas in the province, was selected to investigate exploratory geochemistry by applying statistical methods.

- Geomorphology:

From the morphological viewpoint, the Bards can area can be studied in two separate parts which are distinguished from each other by the Daroone fault. The part of the area located above the fault is a mountain range with a wide variety. But the southern part of the mentioned fault, except for the southeast hills, is a continuation of the outcrops of the Uzbak mountain range, with an average elevation of 850 meters above sea level, which occurs at a vast surface of quaternary units such as alluvial terraces, alluvial fans, clay and salt formations [7] .

- Stratigraphy:

Precambrian: This section is located south of Takanar main fault. The wedge form of this section, which is located between the Daroone and Taknar faults. Based on the existence of outcrops from the Taknar-Precambrian formation, and the covering of the Paleozoic and Mesozoic rocks it has been designated as an erosional window indicating the uplift of Precambrian Paleozoic basement rock of the Iran central tectonic zone in the Tertiary Age [8] .

Taknar formation: The Taknar formation consists of a thick sequence of schists, tuffs, green schists, and quartzite sandstones that have undergone a mild metamorphism in the sub-greenschist facies. In this sequence there are massive metamorphic rhyolites and rhyodacites. One of the important features of this formation is its contact with two intrusive masses. One Intrusive Mass of Precambrian that includes granite and granitoid, and another mass of granite whose time of influence is Eocene-Oligocene [9] [10] .

The type of sampling is lithogeochemical and is performed according to a regular network. A total of 483 lithogeochemical samples were collected from the area. Samples have been analyzed using ICP-AES method. Sampling network location is shown in Figure 2.

Correlation is used to test relationships between quantitative variables or categorical variables. In other words, it’s a measure of how things are related. The study of how variables are correlated is called correlation analysis [11] .

Correlations are useful because if you can find out what relationship variables have, you can make predictions about future behavior [12] . Correlation analysis in geochemical data is very important. One of the valuable results of the investigation of the correlation between chemical elements is the identification of the mineralization and geochemical behavior of the region [13] .

In data mining and statistics, hierarchical clustering (also called hierarchical cluster analysis or HCA) is a method of cluster analysis which seeks to build a hierarchy of clusters. Strategies for hierarchical clustering generally fall into two types [14] :

• Agglomerative Method:

This is a “bottom up” approach: each observation starts in its own cluster, and pairs of clusters are merged as one moves up the hierarchy.

• Divisive Method:

This is a “top down” approach: all observations start in one cluster, and splits are performed recursively as one moves down the hierarchy (see Figure 3).

Discriminant analysis is used as a tool for the separation of communities [15] . This method is described by the number of categories that is possessed by the dependent variable [16] . Linear discriminant analysis works under the simplifying assumption that Σ k = Σ for each class k. In other words, the classes share a common within-class covariance matrix. Since x Σ x term is constant across classes, this simplifies the discriminant function to a linear classifier, δ k ( x ) are known as discriminant functions [17] [18] :

δ k ( x ) = − μ k Σ − 1 x + 1 2 μ k Σ − 1 μ k + log (πk)

you can see the steps in discriminant analysis method in Figure 4.

The following plot shows the linear classification boundaries that result when a sample data set of two variables is modelled using linear discriminant analysis (Figure 5).

The results of the samples analysis which were taken regularly from rock units were analyzed using single-variable statistics. Regarding the field limitations of the elements in the area as well as their measured value in the samples, Ti, Cu, Mo and Zn elements are known as anomalous elements in the region.The result of this study is presented in Table 1.

Using Spearman method, correlation coefficients among the elements in the region were investigated. Table 2 shows the correlation of the elements with each other.

Considering the need for clustering and creating enough visibility to understand mineral elements and also the separation of these elements from the elements that have created the area’s rocks, hierarchical cluster method was used (Figure 6). This clustering can be done by different methods. In this paper, the squared

Euclidean distance was used. Titanium, iron, and magnesium appear together, which reflects the mineralization in the mafic bed rock. In the correlation coefficient (Table 2), correlation of these elements was also observed.

Regarding the choice of the titanium element as target element, its histogram was drawn and the type of distribution of the statistical society of that, was detected as log-normal distribution (Figure 7). The cumulative probability curve for the titanium element is given below in Figure 8.

The separation of communities from the titanium probability chart is considered to be an important point in determining the limits of the society and thus determining the geochemical threshold. For this purpose, a line is fitted to the probability diagram, and the fracture points of the diagram are investigated. Figure 9 illustrates the implications of this fit and determination of communities. Important communities are starting at 1000, 3000, and 9000 ppm which are considered as background, weak anomaly and strong anomaly of titanium in the region.

In this section, LDA was performed using the limits obtained in the previous section. Figure 10 shows the separated communities. In Figure 11, the identified

communities are considered with consideration of all the elements and finally the accuracy of this separation is given in Table 3. In Figure 10, by applying the Linear Discriminant Analysis (LDA) method, three general communities are

separated. In Figure 11, these communities are shown with details and members. The elements in these communities are as follows:

Ÿ 1^st community: Ti.

Ÿ 2^nd community: Be, Fe, Ba, Hg, Mn.

Ÿ 3^rd community: Ni, Zn, W, Ag, Bi, Au, Sn, Cr, Pb, Sr, As, Mo, Cu, Co, Sb.

In Table 3, validation of each of the separated communities in the analysis of differentiation is shown in relation to the threshold values of the titanium element. According to this table, it has shown a high degree of accuracy for this grouping.

After performing various affairs and identifying the relationship between the elements and also recognizing the titanium element as important element in the region. Regarding the behavior of this element, which has been studied in different parts, in the form of anomalies, a map of the highest concentration limits in the region should be prepared. This map represents the best locations for detailed exploration and further exploration. The Kriging Interpolation method (with a variogram that was extracted from a radial survey) was used to prepare this map.

The map is shown in Figure 12. This map, from the blue (the lowest concentration of Ti) to the red (the highest concentration of Ti), is shown in the northern part of an anomaly and in the central part of a small anomaly. It should be noted that this map can be further investigated by changing the interpolation methods and analyzing the locations of the anomalies according to the geological map.