Home Application of multivariate regression on magnetic data to determine further drilling site for iron exploration
Article Open Access

Application of multivariate regression on magnetic data to determine further drilling site for iron exploration

  • Faranak Feizi EMAIL logo , Amir Abbas Karbalaei-Ramezanali and Sasan Farhadi
Published/Copyright: February 17, 2021
Become an author with De Gruyter Brill

Abstract

In this study, a new approach of the multivariate regression model has been applied to make a precise mathematical model to determine further drilling for the detailed iron exploration in the Koohbaba area, Northwest of Iran. Furthermore, to figure out the additional drilling locations, the ore length to the total core ratio for the drilled boreholes has been used based on the geophysical exploration dataset. Hence, different regression analyses including linear, cubic, and quadratic models have been applied. In this study, the ore length to the total core ratio of the chosen drilled boreholes has been considered as a dependent variable; besides, the outputs of the magnetic data using the UP10 (10m upward-continuation), RTP (reduction to the pole), and A.S. (analytic signal) techniques have been designated as independent variables. Based on probability value (p-value), coefficients of determination (R 2 and R adj 2 ), and efficiency formula (EF), the fourth regression model has revealed the best results. The accuracy of the model has been confirmed by the defined ratio of boreholes and demonstrated by four additional drilled boreholes in the study area. Therefore, the results of the regression analysis are reasonable and can be used to determine the additional drilling for the detailed exploration.

1 Introduction

The modern mineral exploration is the definitive aim for the geophysical examination. Hence, several geophysical maps should generate to examine the underground mineral perception [1]. To determine the best zone, drilling of some boreholes is fundamental. Although the most reliable examination for the deposit potential is the drilling [2], it is the most expensive procedure of the mineral exploration [3]. Therefore, using proper methods is essential to decrease the drilling risks and improve the accuracy of the drilling sites [4,5]. Statistical methods play an important role to enhance the success rate and overcome the cost of mineral exploration [6,7,8].

In the past few years, the quantitative study of geoscientific data has increased rapidly [9]. There are several probabilistic, statistical, and mining models proposed for mineral exploration [4,5], such as logistic regression [2,10], ridge regression [11], multitemporal nonlinear regression [12], weight of evidence [13,14], artificial neural networks [15], Bayesian network classifiers [16], and multiple regression analyses [17,18]. All of these methods and techniques have shown promising results and have applied successfully to mineral resource appraisal [1].

Multivariate regression analysis is used successfully to model subsurface mineralization based on the geochemical dataset [19,20], iron mineral mapping [5], and rock properties predictions [21] and to improve prediction model to estimate the sediment yield [22]. In this study, a novel application of the multivariate regression method is proposed to determine additional drilling based on the geophysical exploration dataset. To this end, six different types of multivariate regression have been employed for the iron exploration by using several techniques including RTP (reduction to the pole), A.S. (analytic signal), and UP10 (Upward Continuation) to designate magnetic susceptibility concentrations. The outcomes of these methods have been compared with the log reports from eight different boreholes in the Koohbaba (Qoja-Kandi) area, and the results have demonstrated the proper accuracy of the technique. The log reports of the boreholes have been used to perform further drilling for iron exploration. To achieve this goal, multivariate regression analysis of the ground magnetic data layers has been performed.

1.1 Geological study

This study has been conducted at the Koohbaba area within the Urumieh-Dokhtar Magmatic Arc (UDMA), located in the Northwest of Iran. The study location has a surface area of 1.44 km2, positioned between 46°59′40″ and 47°1′50″ east longitude and 36°51′30″ and 36°53′40″ north latitude, East Azerbaijan (Figure 1). Magmatic activity in UDMA was originated in the Eocene and continued to Pliocene with its climax in the Middle Eocene [23]. Moreover, concerning some recent research works [24,25], UDMA was dominated by the Eocene magmatic rocks, and this fact was confirmed by the geochemical analysis [26].

Figure 1 
                  The Koohbaba position in the physiographic-tectonic zoning map [27].
Figure 1

The Koohbaba position in the physiographic-tectonic zoning map [27].

The UDMA has been dominated by plutonic rocks together with felsic volcanic rocks [23] and forms intrusive–exclusive complex with over 4 km thickness [28,29] and [30]. There is a wide range of composition in the study area, such as schist and shale (Kahar formation), dolomite and limestone (Elika formation), shale, sandstone and limestone (Shemshak formation), limestone, marl, sand stone, conglomerate, and andesite. Omrani et al. [26] have explained that UDMA volcanic rocks form a wide range in composition, which include andesite, minor basalts, and dacites. Magnetite associated with andesite units caused iron mineralization in the study area. The UDMA hosts large metal deposits such as iron and copper [32,33], as shown in Figure 2. Moreover, Mansouri et al. [31] have introduced some iron ore deposits close to the research area.

Figure 2 
                  (a) Cu and Fe outcrops of the Urumieh-Dokhtar Magmatic Arc (UDMA) and (b) magnetite outcrops in the study area.
Figure 2

(a) Cu and Fe outcrops of the Urumieh-Dokhtar Magmatic Arc (UDMA) and (b) magnetite outcrops in the study area.

1.2 Multivariate regression

The regression analysis shows the relation between one or more responses (dependent variables) and one or more predictors (independent variables) and also predicts the values of the responses for a given set of predictors. Usually, these variables are quantitative, i.e., interval or ratio. The following mathematical formula can express the simple relationship between those variables:

(1) Y = F ( X i ) ,

where Y represents dependent variables and X expresses the independent variables. It would be a linear regression if Y and X i have a linear relation; otherwise, it would be a nonlinear regression. The main aim of the regression analysis is to represent the dependent variables as an independent variable function [34]. In this section, the mathematics approach of this technique is presented, and an in-depth presentation of this subject is available in the literature [35].

Multivariate regression is a beneficial statistical method to evaluate the linear relationships between several independent and dependent variables [5]. Therefore, it is the multiple regression expansion with an equal number of equations as the number of response variables. One advantage of using multivariate analysis is that the type 1 error can be determined, and it does not cover the number of variables [34]. Consequently, this method is applied in this study. The linear regression model can be expressed as follows [36,37]:

(2) Y i = β ( X i ) + ε i i = 1 , 2 , , n

where β is a p-dimensional column vector of unknown regression coefficients, Y i is the ith component of the n-dimensional column vector Y, X i is the ith row of the n × p design matrix X, and ε is the random error. The random error value indicates the amount of dispersion in the estimation of the Y value [19]. With n independent observation on Y and the associated values of Z i , the complete model can be expressed as follows [38,39]:

(3) Y 1 = β 0 + β 1 X 11 + β 2 X 12 + + β k X 1 k + ε 1 Y 2 = β 0 + β 1 X 21 + β 2 X 22 + + β k X 2 k + ε 2 Y n = β 0 + β 1 X n 1 + β 2 X n 2 + + β k X n k + ε n ,

In the matrix form, it becomes:

(4) Y 1 Y 2 Y n = 1 X 11 X 12 X 1 k 1 X 21 X 22 X 2 k 1 X n 1 X n 2 X n k β 1 β 2 β n + ε 1 ε 2 ε n .

The error assumptions are as the following:

  1. E(ε) = 0

  2. Cov(ε) = σ 2 I

where I is the n × n identity matrix. The ordinary least-square method has been used to estimate the β matrix [40,41], which is the regression coefficient matrix (Formula 5).

(5) [ β ] = ( [ X ] T [ X ] ) 1 [ X ] T [ Y ]

In multivariate regression, the relationship between the dependent variable Y and the independent variables X i is measured by the coefficients of determination [35]. This is the most frequent statistical approach to estimate the fitness of the model [42]. The mathematical expressions are as follows:

(6) R 2 = SSR SST = 1 SSE SST .

where SSR shows the regression sum of squares, SSE represents the residual (error) sum of squares, and SST is the total sum of squares total. The model can fit the data more reliable as the value of R 2 increases [43]. The R 2 value may increase by adding a new independent variable to the function, although its presence can be required or not. Therefore, it looks complicated whether an increase in R 2 value is significant. In that case, the adjusted determination of coefficients is expressed as follows [19]:

(7) R adj 2 = 1 ( 1 R 2 ) ( n 1 ) ( n p )

where n is the number of samples (data) and k represents the number of coefficients. The adjusted coefficient of determination R adj 2 has been submitted to decrease the bias in R 2 [42]. By adding a different independent variable to the regression model, if the mean square residual (MSR) reduces, the R adj 2 will increase.

2 Materials and methods

2.1 Data analysis

Ground magnetic data were collected by following the same method used by Mansouri et al. [31] in the Koohbaba area. The required geophysical data have provided by magnetometer GSM-19T in the research region (±0.2 nT absolute accuracy). As shown in Figure 3, the total magnetic intensity (TMI) represents the magnetic anomalies in the E–W direction in the north and center of the site. In general, there are three dipolar magnetic anomalies (one magnetic anomaly in the north and two magnetic anomalies in the west and east of the center). These three dipolar magnetic anomalies are 130 related to magnetite dikes in andesite units [5]. Accordingly, the multivariate regression method has been used to determine further drilling sites for iron exploration.

Figure 3 
                  TMI map of the study area with the location of eight drilled boreholes.
Figure 3

TMI map of the study area with the location of eight drilled boreholes.

2.1.1 Preparing dependent variables

The output of the drilled boreholes has been considered as the dependent variable because this variable reveals the accuracy for defining the drilling points. The value for the ratio Ore Length Total Core is between 0 and 1, and for the best boreholes, it is closer to one. The log reports have been collected, and the magnetite thickness of each has been measured respect to the total core (Table 1). The accepted boundary for the ore length is 20% of the total Fe. Besides, Table 2 presents the statistical factors of the ratio using the regression analysis.

Table 1

Log report of boreholes with RTP classification based on fractal method [31]

Borehole id Total core (m) Magnetite range (m) (grades greater than 20% Fe total) Magnetite thickness (m) in total core (grades greater than 20% Fe total) Ore length/total core
From To
BH1 136.5 19.3 25.2 52.4 0.38
60.7 85.2
109.4 131.4
BH2 171.2 4 12.2 47.2 0.27
50.2 53.5
130.6 166.3
BH3 151.2 80 102 32 0.21
112 122
BH4 106 44 48 12.5 0.11
81 89.5
BH5 58.9 0 0
BH6 136.5 69 72 3 0.02
BH7 172 44 47 14 0.08
61.5 63.5
156 164
BH8 157 70 90 29 0.18
133 142
Table 2

The statistical parameters of dependent (ore length/total core) and unique independent (UP10, RTP, and A.S.) variables used in the regression modeling

Variables N Mean St. Dev Skewness Kurtosis
Dependent Sample Y ore length/total core 8 0.156 0.13 0.522 −0.436
Unique independent Raster maps (pixel) X 1 UP10 8 52261.8 3464.8 0.065 −1.273
X 2 RTP 8 56116.5 7494.7 0.170 −1.793
X 3 A.S. 8 279.7 276.04 0.356 −2.035

2.1.2 Preparing independent variables

In the regression analysis, selecting the independent variable is essentially important as these variables must be relevant to the models. Therefore, to make the model, three geophysical raster maps such as upward continuation (UP10), reduce to pole (RTP), and analytic signals (A.S.) have been generated by using Oasis Montaj V.8.4. The upward continuation method is a proper method for deep and semi-deep iron and porphyry-Cu deposit exploration [44]. This method distinguishes the magnetic field far away from the source, and it can decrease the effect of shallow magnetic frequency to create a better map. Figure 4a shows the UP10 map in the Koohhbaba area.

Figure 4 
                     Raster maps of the Koohbaba area: (a) upward continued to 10 m (UP10), (b) reduction to the pole (RTP), and (c) analytic signal (A.S.).
Figure 4

Raster maps of the Koohbaba area: (a) upward continued to 10 m (UP10), (b) reduction to the pole (RTP), and (c) analytic signal (A.S.).

The RTP approach converts magnetic anomalies to a symmetrical pattern, and it can make the magnetic anomalies shape with higher accuracy to the spatial site. Therefore, the magnetic anomalies can interpret much easier [44,45]. In this study, the TMI map has been converted to RTP by using a magnetic declination (4.93) and inclination (55.43). Figure 4b shows the RTP raster map of the study location. The A.S. technique is a wildly known filter to enhance the magnetic field and for locating the magnetic anomalies edges. The basic concept of this method has been discussed in the literature [29,46,47]. The A.S. raster map is shown in Figure 4c.

After generating maps, the values of UP10, RTP, and A.S. raster layers have been extracted in the exact location of eight boreholes (dependent variable); consequently, the values of three unique independent variables (UIVs) have been obtained. The statistical parameters of UIVs (UP10, RTP, and A.S.) have been used in the regression modeling (Table 2).

2.2 Regression analysis

In the regression analysis such as linear regression, to have the best model, it is essential to create many models [19]. In this study, six types of multivariate regression analyses, including linear, quadratic, and cubic models, have been employed (Table 3) to find out the best drilling points for iron exploration.

Table 3

Multivariate regression formula used for detailed iron exploration

Types of regression Number of coefficients Formula
First degree 4 Y 1 = a 0 + a 1 x 1 + a 2 x 2 + a 3 x 3 + ε
First degree 7 Y 2 = Y 1 + a 4 x 1 x 2 + a 5 x 1 x 3 + a 6 x 2 x 3
Second degree 10 Y 3 = Y 2 + a 7 x 1 2 + a 8 x 2 2 + a 9 x 3 2
Second degree 13 Y 4 = Y 3 + a 10 x 1 2 x 2 2 + a 11 x 1 2 x 3 2 + a 12 x 2 2 x 3 2
Third degree 16 Y 5 = Y 4 + a 13 x 1 3 + a 14 x 2 3 + a 15 x 3 3
Third degree 19 Y 6 = Y 5 + a 16 x 1 3 x 2 3 + a 17 x 1 3 x 3 3 + a 18 x 2 3 x 3 3

The models became more complex, from Y 1 to Y 6 as the coefficients increased too. To select the best model, there are some criteria considered for the regression analysis. The values for R 2, R adj 2 , and p-value (ANOVA test) for different regression models is presented in Table 4. Besides, the unknown regression coefficients values are presented in Tables 5 and 6. Other independent variables, which are not presented in tables, are excluded, and they do not affect the statistical models.

Table 4

Values for multivariate regression models

Models R 2 R adj 2 p-value (ANOVA)
Y 1 0.910 0.887 0
Y 2 0.93 0.912 0.003
Y 3 0.937 0.917 0.001
Y 4 0.986 0.966 0.001
Y 5 0.986 0.955 0.001
Y 6 0.986 0.948 0.005
Table 5

The calculated coefficients of regression models (1–3)

Model 1 Model 2 Model 3
Variables Coefficients ( a i ) Variables Coefficients ( a i ) Variables Coefficients ( a i )
CST −0.866 CST −0.59 CST −6.333
x 1 5.504 x 1 −3.351 x 2 −0.014
x 2 4.016 x 2 1.628 x 3 −0.022
x 3 4.016 x 3 −3.354 x 2 2 −1.469
x 1 x 3 2.051 x 3 2 2.292
x 1 2 −7.266
Table 6

The calculated coefficients of regression models (4–6)

Model 4 Model 5 Model 6
Variables Coefficients ( a i ) Variables Coefficients ( a i ) Variables Coefficients ( a i )
CST −4.481 CST −2.210 CST −1.413
x 1 −1.540 x 2 6.206 x 1 −2.839
x 2 0.018 x 3 0.016 x 2 6.188
x 3 −0.024 x 1 2 x 3 2 −8.798 x 3 −2.539
x 1 2 x 2 2 −2.012 x 1 3 −5.654 x 1 2 x 3 2 −7.294
x 1 2 x 3 2 6,215 x 2 3 −9.229 x 1 3 x 2 3 −1.299
x 3 3 4.354 x 1 3 x 3 3 −1.974

3 Results and discussions

To select the best model, some criteria have to been considered. First, the computed variance and random error mean values confirmed the acceptable value for all regression models. Furthermore, based on Table 3, the p-value (ANOVA test) of the models is acceptable (≪0.05). Therefore, all these values demonstrated the accuracy of the regression models.

As Table 4 represents, the lowest value for R 2 has been obtained by the first model (Y 1), and the highest one has been achieved by the last three models (Y 4, Y 5, and Y 6). The fourth model is a second-degree function, and it has lower complexity than other models (Table 3). Even though R 2 is a proper parameter for examining the model with the same number of independent variables, it cannot be adequate for comparing the models with different numbers of independent variables as increasing the number of independent variables will increase the R 2 values. Hence, R adj 2 has been computed, and model number 4 has indicated a fitted model with the highest value (0.966). Therefore, based on the results obtained by coefficients of determination, model Y 4 is the most appropriate model in this study, and it can be applied to determine further drilling exploration sites in the study area.

To approve this fact, initially, the efficiency formula (Formula 8) has been computed as follows [48]:

(8) EF = i = 1 n ( P i O i ) 2 i = 1 n ( O ¯ i O i ) 2 ,

where n is the number of observations, O is the observed values, O ¯ mean of absolute value, and P i is the estimated value. The mean value of the observation becomes more reliable than the estimated values if the model efficiency (EF) approaches zero; therefore, to avoid the model limitations, EF values should be closer to 1 to have a workable model [49].

According to the EF values (Table 4), Y 4 is considered as the best model, followed by Y 5 and Y 6. This result confirmed the output result from the regression analysis by considering the coefficients of determination and p-values (ANOVA test). To determine the further drilling sites in the study area (Koohbaba), the raster map is obtained from Y 4. To generate the intended figures, ArcGIS V.10.1 has been employed (with using the raster calculator toolbox). Figure 5a represents the final raster map of the study by considering the fourth regression model.

Figure 5 
               The final raster map of the study area based on the Y
                  4 regression model with drilled boreholes.
Figure 5

The final raster map of the study area based on the Y 4 regression model with drilled boreholes.

Moreover, four new boreholes (borehole No. 9–12) have been drilled in the study area, belonging to class 5. The result is presented in Table 7 and Figure 5. The accepted boundary for the magnetic thickness is 20% of the total Fe. Concerning this information, the final results are very promising and appropriate.

Table 7

Log report of new drilled boreholes

Borehole ID Total core (m) Magnetite range (m) Magnetite thickness (m) Ore length/total core
From To
BH9 170.8 39.3 41.8 69.9 0.41
43.5 81.2
82.6 92.3
136.5 138.4
142.4 160.5
BH10 169.2 9.4 12.1 107.7 0.63
18.2 52.5
59.3 62.1
68.3 102.4
129.2 163
BH11 171.3 19.8 23.2 65.5 0.38
59.8 82.3
108.2 128.4
143.2 162.6
BH12 170.7 3.3 11.1 77.1 0.45
51.1 53.4
59.2 93.4
130.4 163.2

4 Conclusion

The conclusions are presented as follows:

  1. Regression analysis is a proper and direct statistical method to identify the potential favorable drilling exploration sites with high accuracy in Koohbaba, Northwest of Iran.

  2. The application of the multivariate regression analysis has been confirmed in this area. In this research study, multivariate regression has been developed to create a mathematical model (with reasonable accuracy) for iron mineral exploration by using geophysical data as a new approach.

  3. Six different types of multivariate regression models such as two linear, two quadratic, and two cubic equations have been employed to identify the additional drilling area. According to the results of the coefficients of determination (R 2 and R adj 2 ), p-value, and EF, the fourth regression model (quadratic equation) has been the best response, and it has been confirmed by the ratio (ore length/total core) values of the former drilled boreholes and the further drilled boreholes.

  4. The accuracy of the model has been approved by drilling four new boreholes. This additional field investigation has shown promising results.

  5. The results confirm that the regression analysis model using the geophysical data is an effective approach as it can reduce the time and cost of exploration.

  6. Conclusively, 17092.11 m2 of the study area has been considered as a suitable candidate zone for the detailed studying and determining additional drilling for iron exploration in the area of interest.

References

[1] Oh HJ, Lee S. Regional probabilistic and statistical mineral potential mapping of gold–silver deposits using GIS in the Gangreung area, Korea. Resour Geol. 2008 Jun;58(2):171–87.10.1111/j.1751-3928.2008.00050.xSearch in Google Scholar

[2] Harris DV, Pan G. Mineral favorability mapping: a comparison of artificial neural networks, logistic regression, and discriminant analysis. Nat Resour Res. 1998 Aug;8(2):17.Search in Google Scholar

[3] Marjoribanks RW. Geological methods in mineral exploration and mining. 2nd edn. Berlin, New York: Springer; 2009.10.1007/978-3-540-74375-0Search in Google Scholar

[4] Xiong Y, Zuo R, Carranza EJM. Mapping mineral prospectivity through big data analytics and a deep learning algorithm. Ore Geol Rev. 2018 Nov;102:811–7.10.1016/j.oregeorev.2018.10.006Search in Google Scholar

[5] Mansouri E, Feizi F, Rad AJ, Arian M. Remote-sensing data processing with the multivariate regression analysis method for iron mineral re- source potential mapping: a case study in the Sarvian area, central Iran. Solid Earth. 2018 Mar;9(2):373–84.10.5194/se-9-373-2018Search in Google Scholar

[6] Chen Y, Wu W. Mapping mineral prospectivity using an extreme learning machine regression. Ore Geol Rev. 2017 Jan;80:200–13.10.1016/j.oregeorev.2016.06.033Search in Google Scholar

[7] Ramezanali AK, Feizi F, Jafarirad A, Lotfi M. Application of best-worst method and additive ratio assessment in mineral prospectivity mapping: a case study of vein-type copper mineralization in the Kuhsiah- e-Urmak area, Iran. Ore Geol Rev. 2020 Feb;117:103268.10.1016/j.oregeorev.2019.103268Search in Google Scholar

[8] Feizi F, KarbalaeiRamezanali A, Mansouri E. Calcic iron skarn prospectivity mapping based on fuzzy AHP method, a case study in Varan area, Markazi province, Iran. Geosci J. 2017 Feb;21(1):123–36.10.1007/s12303-016-0042-9Search in Google Scholar

[9] Chen C, Dai H, Liu Y, He B. Mineral prospectivity mapping integrating multi-source geology spatial data sets and logistic regression modelling. Proceedings 2011 IEEE international conference on spatial data mining and geo- graphical knowledge services. Fuzhou, China: IEEE; 2011 Jun. p. 214–710.1109/ICSDM.2011.5969034Search in Google Scholar

[10] Xiong Y, Zuo R. GIS-based rare events logistic regression for mineral prospectivity mapping. Comput Geosci. 2018 Feb;111:18–25.10.1016/j.cageo.2017.10.005Search in Google Scholar

[11] Hang R, Liu Q, Song H, Sun Y, Zhu F, Pei H. Graph regularized nonlinear ridge regression for remote sensing data analysis. IEEE J Sel Top Appl Earth Observ Remote Sens. 2017 Jan;10(1):277–85.10.1109/JSTARS.2016.2574802Search in Google Scholar

[12] Kim HJ, Seo DK, Eo YD, Jeon MC, Park WY. Multi- temporal nonlinear regression method for landsat image simulation. KSCE J Civ Eng. 2019 Feb;23(2):777–87.10.1007/s12205-018-1157-5Search in Google Scholar

[13] Tangestani MH, Moore F. Porphyry copper potential mapping using the weights- of- evidence model in a GIS, northern Shahr-e-Babak, Iran. Aust J Earth Sci. 2001 Oct;48(5):695–701.10.1046/j.1440-0952.2001.485889.xSearch in Google Scholar

[14] Agterberg FP, Bonham-Carter GF. Measuring the performance of mineral-potential maps. Nat Resour Res. 2005 Mar;14(1):1–17.10.1007/s11053-005-4674-0Search in Google Scholar

[15] Singer DA, Kouda R. Application of a feedforward neural network in the search for Kuroko deposits in the Hokuroku district, Japan. Math Geol. 1996 Nov;28(8):1017–23.10.1007/BF02068587Search in Google Scholar

[16] Porwal A, Carranza EJM, Hale M. Bayesian network classifiers for mineral potential mapping. Comput Geosci. 2006 Feb;32(1):1–16.10.1016/j.cageo.2005.03.018Search in Google Scholar

[17] Carranza EJM. Catchment basin modelling of stream sediment anomalies revisited: incorporation of EDA and fractal analysis. Geochem Explor Environ Anal. 2010 Nov;10(4):365–81.10.1144/1467-7873/09-224Search in Google Scholar

[18] Carranza EJM. Mapping of anomalies in continuous and discrete fields of stream sediment geochemical landscapes. Geochem Explor Enviro Anal. 2010 May;10(2):171–87.10.1144/1467-7873/09-223Search in Google Scholar

[19] Granian H, Tabatabaei SH, Asadi HH, Carranza EJM. Multivariate regression analysis of lithogeochemical data to model subsurface mineralization: a case study from the Sari Gunay epithermal gold deposit, NW Iran. J Geochem Explor. 2015 Jan;148:249–58.10.1016/j.gexplo.2014.10.009Search in Google Scholar

[20] Ramezanali AK, Feizi F, Jafarirad A, Lotfi M. Geochemical anomaly and mineral prospectivity mapping for vein-type copper mineralization, Kuhsiah-e-Urmak area, Iran: application of sequential gaussian simulation and multivariate regression analysis. Nat Resour Res. 2020 Feb;29(1):41–70.10.1007/s11053-019-09565-7Search in Google Scholar

[21] Ma YZ. Pitfalls in predictions of rock properties using multivariate analysis and regression methods. J Appl Geophys. 2011;75:390–400.10.1016/j.jappgeo.2011.06.022Search in Google Scholar

[22] Grauso S, Pasanisi F, Tebano C, Grillini M, Peloso A. Investigating the sediment yield predictability in some Italian rivers by means of hydro-geomorphometric variables. Geosciences. 2018 May;8:249.10.3390/geosciences8070249Search in Google Scholar

[23] Kananian A, Sarjoughian F, Nadimi A, Ahmadian J, Ling W. Geochemical characteristics of the Kuh-e Dom intrusion, Urumieh–Dokhtar magmatic arc (Iran): implications for source regions and magmatic evolution. J Asian Earth Sci. 2014 Aug;90:137–48.10.1016/j.jseaes.2014.04.026Search in Google Scholar

[24] Yeganehfar H, Ghorbani MR, Shinjo R, Ghaderi M. Mag- matic and geodynamic evolution of Urumieh–Dokhtar basic volcanism, Central Iran: major, trace element, isotopic, and geochronologic implications. Int Geol Rev. 2013 April;55(6):767–86.10.1080/00206814.2012.752554Search in Google Scholar

[25] Babazadeh S, Ghorbani MR, Cottle JM, Bröcker M. Multistage tectono-magmatic evolution of the central Urumieh-Dokhtar magmatic arc, south Ardestan, Iran: insights from zircon geochronology and geochemistry. Geol J. 2019 July;54((4):2447–71.10.1002/gj.3306Search in Google Scholar

[26] Omrani J, Agard P, Whitechurch H, Benoit M, Prouteau G, Jolivet L. Arc-magmatism and subduction history beneath the Zagros mountains, Iran: a new report of adakites and geodynamic consequences. Lithos. 2008 Dec;106(3–4):380–98.10.1016/j.lithos.2008.09.008Search in Google Scholar

[27] Arian M. Physiographic-tectonic zoning of Iran’s sedimentary basins. Open J Geol. 2013;3(3):169–77.10.4236/ojg.2013.33020Search in Google Scholar

[28] Alavi M. Tectonics of the Zagros erogenic belt of Iran: new data and interpretations. Tectonophysics. 1994;229:211–38.10.1016/0040-1951(94)90030-2Search in Google Scholar

[29] Ramezanali AK, Mansouri E, Faranak F. Integration of aeromagnetic geophysical data with other exploration data layers based on fuzzy AHP and C-A fractal model for Cu-porphyry potential mapping: a case study in the Fordo area, central Iran. Boll di Geofisica Teorica ed Applicata. 2017;58(1):55–73.Search in Google Scholar

[30] Feizi F, Mansouri E, Ramezanali AK. Prospecting of Au by remote sensing and geochemical data processing using fractal modelling in Shishe-Botagh, area (NW Iran). J Indian Soc Remote Sens. 2016 Aug;44(4):539–52.10.1007/s12524-015-0510-0Search in Google Scholar

[31] Mansouri E, Feizi F, Karbalaei-Ramezanali A. Identification of magnetic anomalies based on ground magnetic data analysis using multifractal modelling: a case study in Qoja-Kandi, East Azerbaijan province, Iran. Nonlinear Process Geo-Phys. 2015 Oct;22(5):579–87.10.5194/npg-22-579-2015Search in Google Scholar

[32] Feizi F, Karbalaei-Ramezanali A, Tusi H. Mineral potential mapping via TOPSIS with hybrid AHP–shannon entropy weighting of evidence: a case study for porphyry-Cu, Farmahin area, Markazi province, Iran. Nat Resour Res. 2017 Oct;26(4):553–70.10.1007/s11053-017-9338-3Search in Google Scholar

[33] Feizi F, Mansuri E. Separation of alteration zones on ASTER data and integration with drainage geochemical maps in Soltanieh, Northern Iran. Open J Geol. 2013;03(02):134–42.10.4236/ojg.2013.32017Search in Google Scholar

[34] Duleba A, Olive D. Regression analysis and multivariate analysis. Sem Reprod Med. 1996 May;14(2):139–53.10.1055/s-2007-1016322Search in Google Scholar PubMed

[35] Rencher AC. Methods of multivariate analysis. Wiley series in probability and mathematical statistics. 2nd edn. New York: Wiley; 2002.10.1002/0471271357Search in Google Scholar

[36] Chung CF, Agterberg FP. Regression models for estimating mineral resources from geological map data. J Int Assoc Math Geol. oct 1980;12(5):473–88.10.1007/BF01028881Search in Google Scholar

[37] Zhang D. A coefficient of determination for generalized linear models. Am Stat. 2017 Oct;71(4):310–6.10.1080/00031305.2016.1256839Search in Google Scholar

[38] Johnson RA, Wichern DW. Applied multivariate statistical analysis. 6th edn. New Jersey: Peardon; 2007.10.1002/0471667196.ess6094Search in Google Scholar

[39] Chung CJF, Van Westen CJ. Multivariate regression analysis for landslide hazard zonation. Geographical information systems in assessing natural hazards, vol. 5. Netherlands, Dordrecht: Springer; 1995. p. 107–33.10.1007/978-94-015-8404-3_7Search in Google Scholar

[40] Scott AJ, Holt D. The effect of two-stage sampling on ordinary least squares methods. J Am Stat Assoc. 1982;77(380):7.10.1080/01621459.1982.10477897Search in Google Scholar

[41] Chang J, Olive DJ. OLS for 1D regression models. Commun Stat Theory Methods. 2010 May;39(10):1869–82.10.1080/03610920902923494Search in Google Scholar

[42] Akossou AYJ, Palm R. Impact of data structure on the estimators R-square and adjusted R-square in linear regression. Int J Math Comput. 2013;20:10.Search in Google Scholar

[43] Cornell JA. Factors that influence the value of the coefficient of determination in simple linear and nonlinear regression models. Phytopathology. 1987;77(1):63.10.1094/Phyto-77-63Search in Google Scholar

[44] Abedi M, Norouzi GH. Integration of various geophysical data with geological and geochemical data to determine additional drilling for copper explo- ration. J Appl Geophys. 2012 Aug;83:35–45.10.1016/j.jappgeo.2012.05.003Search in Google Scholar

[45] Golshadi Z. Interpretation of magnetic data in the Chenar-e Olya area of Asad- abad, Hamedan, Iran, using analytic signal, euler deconvolution, horizontal gradient and tilt-derivative methods. Boll di Geofisica Teorica ed Applicata. 2016 Dec;57(4):329–42.Search in Google Scholar

[46] Nabighian MN. The analytic signal of two-dimensional magnetic bodies with polygonal cross-section: its properties and use for automated anomaly interpretation. Geophysics. 1972 Jan;37(3):507–17.10.1190/1.1440276Search in Google Scholar

[47] Nabighian MN. Additional comments on the analytic sig- nal of two-dimensional magnetic bodies with polygonal cross- section. Geophysics. 1974 Feb;39(1):85–92.10.1190/1.1440416Search in Google Scholar

[48] Li J, Heap AD. A review of spatial interpolation methods for environ- mental scientists. Canberra: Geoscience Australia; 2008. p. 154.Search in Google Scholar

[49] Vicente-Serrano SM, Saz-Sanchez MA, Cuadrat JM. Comparative analysis of inter- polation methods in the middle Ebro valley (Spain): application to annual precipitation and temperature. Clim Res. 2003;24:161–80.10.3354/cr024161Search in Google Scholar

Received: 2020-05-17
Revised: 2021-01-13
Accepted: 2021-01-16
Published Online: 2021-02-17

© 2021 Faranak Feizi et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Lithopetrographic and geochemical features of the Saalian tills in the Szczerców outcrop (Poland) in various deformation settings
  3. Spatiotemporal change of land use for deceased in Beijing since the mid-twentieth century
  4. Geomorphological immaturity as a factor conditioning the dynamics of channel processes in Rządza River
  5. Modeling of dense well block point bar architecture based on geological vector information: A case study of the third member of Quantou Formation in Songliao Basin
  6. Predicting the gas resource potential in reservoir C-sand interval of Lower Goru Formation, Middle Indus Basin, Pakistan
  7. Study on the viscoelastic–viscoplastic model of layered siltstone using creep test and RBF neural network
  8. Assessment of Chlorophyll-a concentration from Sentinel-3 satellite images at the Mediterranean Sea using CMEMS open source in situ data
  9. Spatiotemporal evolution of single sandbodies controlled by allocyclicity and autocyclicity in the shallow-water braided river delta front of an open lacustrine basin
  10. Research and application of seismic porosity inversion method for carbonate reservoir based on Gassmann’s equation
  11. Impulse noise treatment in magnetotelluric inversion
  12. Application of multivariate regression on magnetic data to determine further drilling site for iron exploration
  13. Comparative application of photogrammetry, handmapping and android smartphone for geotechnical mapping and slope stability analysis
  14. Geochemistry of the black rock series of lower Cambrian Qiongzhusi Formation, SW Yangtze Block, China: Reconstruction of sedimentary and tectonic environments
  15. The timing of Barleik Formation and its implication for the Devonian tectonic evolution of Western Junggar, NW China
  16. Risk assessment of geological disasters in Nyingchi, Tibet
  17. Effect of microbial combination with organic fertilizer on Elymus dahuricus
  18. An OGC web service geospatial data semantic similarity model for improving geospatial service discovery
  19. Subsurface structure investigation of the United Arab Emirates using gravity data
  20. Shallow geophysical and hydrological investigations to identify groundwater contamination in Wadi Bani Malik dam area Jeddah, Saudi Arabia
  21. Consideration of hyperspectral data in intraspecific variation (spectrotaxonomy) in Prosopis juliflora (Sw.) DC, Saudi Arabia
  22. Characteristics and evaluation of the Upper Paleozoic source rocks in the Southern North China Basin
  23. Geospatial assessment of wetland soils for rice production in Ajibode using geospatial techniques
  24. Input/output inconsistencies of daily evapotranspiration conducted empirically using remote sensing data in arid environments
  25. Geotechnical profiling of a surface mine waste dump using 2D Wenner–Schlumberger configuration
  26. Forest cover assessment using remote-sensing techniques in Crete Island, Greece
  27. Stability of an abandoned siderite mine: A case study in northern Spain
  28. Assessment of the SWAT model in simulating watersheds in arid regions: Case study of the Yarmouk River Basin (Jordan)
  29. The spatial distribution characteristics of Nb–Ta of mafic rocks in subduction zones
  30. Comparison of hydrological model ensemble forecasting based on multiple members and ensemble methods
  31. Extraction of fractional vegetation cover in arid desert area based on Chinese GF-6 satellite
  32. Detection and modeling of soil salinity variations in arid lands using remote sensing data
  33. Monitoring and simulating the distribution of phytoplankton in constructed wetlands based on SPOT 6 images
  34. Is there an equality in the spatial distribution of urban vitality: A case study of Wuhan in China
  35. Considering the geological significance in data preprocessing and improving the prediction accuracy of hot springs by deep learning
  36. Comparing LiDAR and SfM digital surface models for three land cover types
  37. East Asian monsoon during the past 10,000 years recorded by grain size of Yangtze River delta
  38. Influence of diagenetic features on petrophysical properties of fine-grained rocks of Oligocene strata in the Lower Indus Basin, Pakistan
  39. Impact of wall movements on the location of passive Earth thrust
  40. Ecological risk assessment of toxic metal pollution in the industrial zone on the northern slope of the East Tianshan Mountains in Xinjiang, NW China
  41. Seasonal color matching method of ornamental plants in urban landscape construction
  42. Influence of interbedded rock association and fracture characteristics on gas accumulation in the lower Silurian Shiniulan formation, Northern Guizhou Province
  43. Spatiotemporal variation in groundwater level within the Manas River Basin, Northwest China: Relative impacts of natural and human factors
  44. GIS and geographical analysis of the main harbors in the world
  45. Laboratory test and numerical simulation of composite geomembrane leakage in plain reservoir
  46. Structural deformation characteristics of the Lower Yangtze area in South China and its structural physical simulation experiments
  47. Analysis on vegetation cover changes and the driving factors in the mid-lower reaches of Hanjiang River Basin between 2001 and 2015
  48. Extraction of road boundary from MLS data using laser scanner ground trajectory
  49. Research on the improvement of single tree segmentation algorithm based on airborne LiDAR point cloud
  50. Research on the conservation and sustainable development strategies of modern historical heritage in the Dabie Mountains based on GIS
  51. Cenozoic paleostress field of tectonic evolution in Qaidam Basin, northern Tibet
  52. Sedimentary facies, stratigraphy, and depositional environments of the Ecca Group, Karoo Supergroup in the Eastern Cape Province of South Africa
  53. Water deep mapping from HJ-1B satellite data by a deep network model in the sea area of Pearl River Estuary, China
  54. Identifying the density of grassland fire points with kernel density estimation based on spatial distribution characteristics
  55. A machine learning-driven stochastic simulation of underground sulfide distribution with multiple constraints
  56. Origin of the low-medium temperature hot springs around Nanjing, China
  57. LCBRG: A lane-level road cluster mining algorithm with bidirectional region growing
  58. Constructing 3D geological models based on large-scale geological maps
  59. Crops planting structure and karst rocky desertification analysis by Sentinel-1 data
  60. Physical, geochemical, and clay mineralogical properties of unstable soil slopes in the Cameron Highlands
  61. Estimation of total groundwater reserves and delineation of weathered/fault zones for aquifer potential: A case study from the Federal District of Brazil
  62. Characteristic and paleoenvironment significance of microbially induced sedimentary structures (MISS) in terrestrial facies across P-T boundary in Western Henan Province, North China
  63. Experimental study on the behavior of MSE wall having full-height rigid facing and segmental panel-type wall facing
  64. Prediction of total landslide volume in watershed scale under rainfall events using a probability model
  65. Toward rainfall prediction by machine learning in Perfume River Basin, Thua Thien Hue Province, Vietnam
  66. A PLSR model to predict soil salinity using Sentinel-2 MSI data
  67. Compressive strength and thermal properties of sand–bentonite mixture
  68. Age of the lower Cambrian Vanadium deposit, East Guizhou, South China: Evidences from age of tuff and carbon isotope analysis along the Bagong section
  69. Identification and logging evaluation of poor reservoirs in X Oilfield
  70. Geothermal resource potential assessment of Erdaobaihe, Changbaishan volcanic field: Constraints from geophysics
  71. Geochemical and petrographic characteristics of sediments along the transboundary (Kenya–Tanzania) Umba River as indicators of provenance and weathering
  72. Production of a homogeneous seismic catalog based on machine learning for northeast Egypt
  73. Analysis of transport path and source distribution of winter air pollution in Shenyang
  74. Triaxial creep tests of glacitectonically disturbed stiff clay – structural, strength, and slope stability aspects
  75. Effect of groundwater fluctuation, construction, and retaining system on slope stability of Avas Hill in Hungary
  76. Spatial modeling of ground subsidence susceptibility along Al-Shamal train pathway in Saudi Arabia
  77. Pore throat characteristics of tight reservoirs by a combined mercury method: A case study of the member 2 of Xujiahe Formation in Yingshan gasfield, North Sichuan Basin
  78. Geochemistry of the mudrocks and sandstones from the Bredasdorp Basin, offshore South Africa: Implications for tectonic provenance and paleoweathering
  79. Apriori association rule and K-means clustering algorithms for interpretation of pre-event landslide areas and landslide inventory mapping
  80. Lithology classification of volcanic rocks based on conventional logging data of machine learning: A case study of the eastern depression of Liaohe oil field
  81. Sequence stratigraphy and coal accumulation model of the Taiyuan Formation in the Tashan Mine, Datong Basin, China
  82. Influence of thick soft superficial layers of seabed on ground motion and its treatment suggestions for site response analysis
  83. Monitoring the spatiotemporal dynamics of surface water body of the Xiaolangdi Reservoir using Landsat-5/7/8 imagery and Google Earth Engine
  84. Research on the traditional zoning, evolution, and integrated conservation of village cultural landscapes based on “production-living-ecology spaces” – A case study of villages in Meicheng, Guangdong, China
  85. A prediction method for water enrichment in aquifer based on GIS and coupled AHP–entropy model
  86. Earthflow reactivation assessment by multichannel analysis of surface waves and electrical resistivity tomography: A case study
  87. Geologic structures associated with gold mineralization in the Kirk Range area in Southern Malawi
  88. Research on the impact of expressway on its peripheral land use in Hunan Province, China
  89. Concentrations of heavy metals in PM2.5 and health risk assessment around Chinese New Year in Dalian, China
  90. Origin of carbonate cements in deep sandstone reservoirs and its significance for hydrocarbon indication: A case of Shahejie Formation in Dongying Sag
  91. Coupling the K-nearest neighbors and locally weighted linear regression with ensemble Kalman filter for data-driven data assimilation
  92. Multihazard susceptibility assessment: A case study – Municipality of Štrpce (Southern Serbia)
  93. A full-view scenario model for urban waterlogging response in a big data environment
  94. Elemental geochemistry of the Middle Jurassic shales in the northern Qaidam Basin, northwestern China: Constraints for tectonics and paleoclimate
  95. Geometric similarity of the twin collapsed glaciers in the west Tibet
  96. Improved gas sand facies classification and enhanced reservoir description based on calibrated rock physics modelling: A case study
  97. Utilization of dolerite waste powder for improving geotechnical parameters of compacted clay soil
  98. Geochemical characterization of the source rock intervals, Beni-Suef Basin, West Nile Valley, Egypt
  99. Satellite-based evaluation of temporal change in cultivated land in Southern Punjab (Multan region) through dynamics of vegetation and land surface temperature
  100. Ground motion of the Ms7.0 Jiuzhaigou earthquake
  101. Shale types and sedimentary environments of the Upper Ordovician Wufeng Formation-Member 1 of the Lower Silurian Longmaxi Formation in western Hubei Province, China
  102. An era of Sentinels in flood management: Potential of Sentinel-1, -2, and -3 satellites for effective flood management
  103. Water quality assessment and spatial–temporal variation analysis in Erhai lake, southwest China
  104. Dynamic analysis of particulate pollution in haze in Harbin city, Northeast China
  105. Comparison of statistical and analytical hierarchy process methods on flood susceptibility mapping: In a case study of the Lake Tana sub-basin in northwestern Ethiopia
  106. Performance comparison of the wavenumber and spatial domain techniques for mapping basement reliefs from gravity data
  107. Spatiotemporal evolution of ecological environment quality in arid areas based on the remote sensing ecological distance index: A case study of Yuyang district in Yulin city, China
  108. Petrogenesis and tectonic significance of the Mengjiaping beschtauite in the southern Taihang mountains
  109. Review Articles
  110. The significance of scanning electron microscopy (SEM) analysis on the microstructure of improved clay: An overview
  111. A review of some nonexplosive alternative methods to conventional rock blasting
  112. Retrieval of digital elevation models from Sentinel-1 radar data – open applications, techniques, and limitations
  113. A review of genetic classification and characteristics of soil cracks
  114. Potential CO2 forcing and Asian summer monsoon precipitation trends during the last 2,000 years
  115. Erratum
  116. Erratum to “Calibration of the depth invariant algorithm to monitor the tidal action of Rabigh City at the Red Sea Coast, Saudi Arabia”
  117. Rapid Communication
  118. Individual tree detection using UAV-lidar and UAV-SfM data: A tutorial for beginners
  119. Technical Note
  120. Construction and application of the 3D geo-hazard monitoring and early warning platform
  121. Enhancing the success of new dams implantation under semi-arid climate, based on a multicriteria analysis approach: Case of Marrakech region (Central Morocco)
  122. TRANSFORMATION OF TRADITIONAL CULTURAL LANDSCAPES - Koper 2019
  123. The “changing actor” and the transformation of landscapes
Downloaded on 19.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/geo-2020-0165/html?lang=en
Scroll to top button