Document Type : Original Research Paper
Authors
- Marco Antonio Cotrina-Teatino 1, 2
- Jairo Jhonatan Marquina-Araujo 1, 2
- Jose Nestor Mamani-Quispe 3, 2
- Solio Marino Arango-Retamozo 1, 2
- Joe Alexis Gonzalez-Vasquez 4
1 Department of Mining Engineering, Faculty of Engineering, National University of Trujillo, Trujillo, Peru
2 Laboratory for Mine Design and Planning – DIPLAMIN, National University of Trujillo, Trujillo, Peru
3 Faculty of Chemical Engineering, National University of the Altiplano of Puno, Puno, Peru
4 Department of Industrial Engineering, Faculty of Engineering, National University of Trujillo, Trujillo, Peru
10.22044/jme.2025.16188.3127
Abstract
Traditional geostatistical methods such as kriging exhibit limitations by assuming linear and symmetric dependencies, which can lead to smoothed estimates and the loss of local variability. To address these issues, this study applies Archimedean copulas (Clayton, Gumbel, and Frank) for the estimation of copper ore grades in a deposit located in Peru. A total of 5,654 composites, each 15 meters in length, were obtained from 185 diamond drill holes. The data were transformed to a uniform scale to allow for copula fitting. Dependence structures were modeled by lag distance, with the dependence parameter fitted using fifth-degree polynomials, and three-dimensional conditional estimation was implemented. Results indicate that ordinary kriging yielded RMSE = 0.161, MAE = 0.104, R2 = 0.692, and a correlation of 0.861. The Clayton copula slightly improved these metrics (RMSE = 0.154, MAE = 0.101, R2 = 0.717, R = 0.871), while the Gumbel copula captured higher local variability (RMSE = 0.161, MAE = 0.116, R2 = 0.692, R = 0.855). The Frank copula achieved the best performance with RMSE = 0.137, MAE = 0.090, R2 = 0.778, and R = 0.905. In conclusion, Archimedean copulas significantly enhance geostatistical estimation by better capturing spatial dependence, offering a robust alternative to classical geostatistical methods.
Keywords
Main Subjects
[1]. Cressie, N. (1990). The origins of kriging. Mathematical Geology, 22(3), 239-252.
[2]. Lloyd, C., & Atkinson, P. (2001). Assessing uncertainty in estimates with ordinary and indicator kriging. Computers & Geosciences, 27(8), 929-937.
[3]. Gräler, B., & Pebesma, E. (2011). The pair-copula construction for spatial data: A new approach to model spatial dependency. Procedia Environmental Sciences, 7, 206-211.
[4]. Da Rocha, M., & Yamamoto, J. (2000). Comparison between kriging variance and interpolation variance as uncertainty measurements in the Capanema iron mine, State of Minas Gerais, Brazil. Natural Resources Research, 9(3), 233-235.
[5]. Li, D., Zhang, L., Tang, X., Zhou, W., Li, J., Zhou, C., & Phoon, K. (2015). Bivariate distribution of shear strength parameters using copulas and its impact on geotechnical system reliability. Computers and Geotechnics, 68, 184-195.
[6]. Guo, N., Wang, F., & Yang, J. (2017). Remarks on composite Bernstein copula and its application to credit risk analysis. Insurance: Mathematics and Economics, 77, 38-48.
[7]. Lourme, A., & Maurer, F. (2017). Testing the Gaussian and Student’s t copulas in a risk management framework. Economic Modelling, 67, 203-214.
[8]. Zhu, H., Zhang, L., Xiao, T., & Li, X. (2017). Generation of multivariate cross-correlated geotechnical random fields. Computers and Geotechnics, 86, 95-107.
[9]. Gong, Y., Chen, Q., & Liang, J. (2018). A mixed data sampling copula model for the return-liquidity dependence in stock index futures markets. Economic Modelling, 68, 586-598.
[10]. Frahm, G., Junker, M., & Szimayer, A. (2003). Elliptical copulas: Applicability and limitations. Statistics & Probability Letters, 63(3), 275-286.
[11]. Li, D., & Peng, L. (2009). Goodness-of-fit test for tail copulas modeled by elliptical copulas. Statistics & Probability Letters, 79(8), 1097-1104.
[12]. Hashorva, E., & Jaworski, P. (2012). Gaussian approximation of conditional elliptical copulas. Journal of Multivariate Analysis, 111, 397-407.
[13]. Li, C., Huang, Y., & Zhu, L. (2017). Color texture image retrieval based on Gaussian copula models of Gabor wavelets. Pattern Recognition, 64, 118-129.
[14]. Schweizer, B. (2007). Introduction to copulas. Journal of Hydraulic Engineering, 12(4), 346.
[15]. Ravens, B. (2000). An introduction to copulas. Technometrics, 42(3), 317.
[16]. Hougaard, P. (1986). A class of multivariate failure time distributions. Biometrika, 73(3), 671-678.
[17]. Clayton, D. G. (1978). A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence. Biometrika, 65(1), 141-151.
[18]. Ali, M. M., Mikhail, N. N., & Haq, M. S. (1978). A class of bivariate distributions including the bivariate logistic. Journal of Multivariate Analysis, 8(3), 405-412.
[19]. Frank, M. J. (1979). On the simultaneous associativity of F(x, y) and x+y−F(x, y). Aequationes Mathematicae, 19(1), 194-226.
[20]. Shaked, M., & Joe, H. (1998). Multivariate models and dependence concepts. Journal of the American Statistical Association, 93(443), 1237-1238.
[21]. Alzaid, A., & Alhadlaq, W. (2024). A new family of Archimedean copulas: The half-logistic family of copulas. Mathematics, 12(1), 101.
[22]. Marshall, A. W., & Olkin, I. (1988). Families of multivariate distributions. Journal of the American Statistical Association, 83(403), 834-841.
[23]. Genest, C., & Mackay, R. J. (1986). Copules archimédiennes et familles de lois bidimensionnelles dont les marges sont données. Canadian Journal of Statistics, 14(2), 145-159.
[24]. Bárdossy, A. (2006). Copula-based geostatistical models for groundwater quality parameters. Water Resources Research, 42(11), W11416.
[25]. Quessy, J. F., Rivest, L. P., & Toupin, M. H. (2019). Goodness-of-fit tests for the family of multivariate chi-square copulas. Computational Statistics & Data Analysis, 140, 21–40.
[26]. Bedford, T., & Cooke, R. M. (2002). Vines: A new graphical model for dependent random variables. Annals of Statistics, 30(4), 1031–1068.
[27]. Aas, K., Czado, C., Frigessi, A., & Bakken, H. (2009). Pair-copula constructions of multiple dependence. Insurance: Mathematics and Economics, 44(2), 182–198.
[28]. Addo, E., Chanda, E. K., & Metcalfe, A. V. (2019). Spatial pair-copula model of grade for an anisotropic gold deposit. Mathematical Geosciences, 51(5), 553–578.
[29]. Dinda, K., & Samanta, B. (2021). Non-Gaussian copula simulation for estimation of recoverable reserve in an Indian copper deposit. Natural Resources Research, 30(1), 57–76.
[30]. Gräler, B. (2014). Modelling skewed spatial random fields through the spatial vine copula. Spatial Statistics, 10, 87–102.
[31]. Sohrabian, B. (2021). Geostatistical prediction through convex combination of Archimedean copulas. Spatial Statistics, 41, 100488.
[32]. Kazianka, H., & Pilz, J. (2011). Bayesian spatial modeling and interpolation using copulas. Computers & Geosciences, 37(3), 310–319.
[33]. Journel, A. G., & Alabert, F. (1989). Non‐Gaussian data expansion in the Earth sciences. Terra Nova, 1(2), 123–134.
[34]. Haslauer, C. P., Li, J., & Bárdossy, A. (2010). Application of copulas in geostatistics. In: geoENV VII – Geostatistics for Environmental Applications (pp. 395–404). Springer.
[35]. Gnann, S., Allmendinger, M., Haslauer, C., & Bárdossy, A. (2018). Improving copula-based spatial interpolation with secondary data. Spatial Statistics, 28, 105–127.
[36]. Agarwal, G., Sun, Y., & Wang, H. J. (2021). Copula-based multiple indicator kriging for non-Gaussian random fields. Spatial Statistics, 44, 100524.
[37]. Kazianka, H., & Pilz, J. (2010). Copula-based geostatistical modeling of continuous and discrete data including covariates. Stochastic Environmental Research and Risk Assessment, 24(5), 661–673.
[38]. Addo, E., Metcalfe, A. V., Chanda, E. K., Sepulveda, E., Assibey-Bonsu, W., & Adeli, A. (2019). Prediction of copper recovery from geometallurgical data using D-vine copulas. Journal of the Southern African Institute of Mining and Metallurgy, 119(11), 891–898.
[39]. Dinda, K., Samanta, B., & Chakravarty, D. (2022). A v-transformed copula-based simulation model for lithological classification in an Indian copper deposit. Scientific Reports, 12, 13728.
[40]. Sohrabian, B., & Tercan, A. (2024). Copula-based data-driven multiple-point simulation method. Spatial Statistics, 59, 100802.
[41]. Sohrabian, B., & Tercan, A. (2025). Grade estimation through the Gaussian copulas: A case study. Journal of Mining and Environment, 16(1), 1–13.
[42]. Hernández, H., Díaz-Viera, M., Alberdi, E., Oyarbide-Zubillaga, A., & Goti, A. (2024). Metallurgical copper recovery prediction using conditional quantile regression based on a copula model. Minerals, 14(6), 691.
[43]. Krysa, Z., Pactwa, K., Wozniak, J., & Dudek, M. (2017). Using copulas in the estimation of the economic project value in the mining industry, including geological variability. IOP Conference Series: Earth and Environmental Science, 95(4), 042001.
[44]. Xu, D., & Zhu, Y. (2020). A copula–Hubbert model for Co (by)-product minerals. Natural Resources Research, 29(5), 3069–3078.
[45]. Sohrabian, B., Soltani-Mohammadi, S., Pourmirzaee, R., & Carranza, E. (2023). Geostatistical evaluation of a porphyry copper deposit using copulas. Minerals, 13(6), 732.
[46]. Akbari Gharalari, M., Abdollahi-Sharif, J., & Sohrabian, B. (2022). Classification of reserve in Sungun mine based on Archimedean copulas estimates. Arabian Journal of Geosciences, 15, 1695.
[47]. Sotoudeh, F., Ataei, M., Kakaie, R., & Pourrahimian, Y. (2020). Application of sequential Gaussian conditional simulation into underground mine design under grade uncertainty. Journal of Mining and Environment, 11(3), 695–709.
[48]. Parhizkar, A., Ataei, M., Moarefvand, P., & Rasouli, V. (2011). Grade uncertainty and its impact on ore grade reconciliation between the resource model and the mine. Archives of Mining Sciences, 56(1), 119–134.
[49]. Parhizkar, A., Ataei, M., Moarefvand, P., & Rasouli, V. (2012). A probabilistic model to improve reconciliation of estimated and actual grade in open pit mining. Arabian Journal of Geosciences, 5(6), 1279–1288.
[50]. Farhadi, S., Tatullo, S., Konari, M., & Afzal, P. (2024). Evaluating StackingC and ensemble models for enhanced lithological classification in geological mapping. Journal of Geochemical Exploration, 260, 107441.
[51]. Fathi, M., Alimoradi, A., & Hemati Ahooi, H. (2021). Optimizing extreme learning machine algorithm using particle swarm optimization to estimate iron ore grade. Journal of Mining and Environment, 12(2), 397–411.
[52]. Afzal, P., Gholami, H., Madani, N., Yasrebi, A., & Sadeghi, B. (2023). Mineral resource classification using geostatistical and fractal simulation in the Masjed Daghi Cu–Mo porphyry deposit, NW Iran. Minerals, 13(3), 370.
[53]. Salarian, S., Oskooi, B., Mostafaei, K., & Smirnov, M. (2024). Improving the resource modeling results using auxiliary variables in estimation and simulation methods. Earth Science Informatics, 17, 4161–4181.
[54]. Mostafaei, K., & Ramazi, H. (2019). Mineral resource estimation using a combination of drilling and IP-Rs data with statistical and cokriging methods. Bulletin of the Mineral Research and Exploration, 160(160), 177–195.
[55]. Monjezi, M., Rajabalizadeh Kashani, M., & Ataei, M. (2013). A comparative study between sequential Gaussian simulation and kriging method grade modeling in open-pit mining. Arabian Journal of Geosciences, 6(1), 123–128
[56]. Tahernejad, M., Khalokakaie, R., & Ataei, M. (2018). Analyzing the effect of ore grade uncertainty in open pit mine planning: A case study of Rezvan iron mine, Iran. International Journal of Mining and Geo-Engineering, 52(1), 53–60.
[57]. Ghasemitabar, H., Alimoradi, A., Hemati Ahooi, H., & Fathi, M. (2024). Intelligent borehole simulation with Python programming. Journal of Mining and Environment, 15(2), 707–730.
[58]. Marquina, J., Cotrina, M., Mamani, J., Noriega, E., & Vega, J., & Cruz, J. (2024). Copper ore grade prediction using machine learning techniques in a copper deposit. Journal of Mining and Environment, 15(3), 1011–1027.
[59]. Cotrina, M., Marquina, J., & Mamani, J. (2025). Application of artificial neural networks for the categorization of mineral resources in a copper deposit in Peru. World Journal of Engineering.
[60]. Cotrina, M., Marquina, J., Mamani, J., Arango, S., Ccatamayo, J., Gonzalez, J., Donaires, T., & Calla, M. (2025). Categorization of mineral resources using Random Forest model in a copper deposit in Peru. Journal of Mining and Environment, 16(4), 947–962.
[61]. Marquina-Araujo, J., Cotrina-Teatino, M., Cruz-Galvez, J., Noriega-Vidal, E., & Vega-Gonzalez, J. (2024). Application of autoencoders neural network and K-means clustering for the definition of geostatistical estimation domains. Mathematical Modelling of Engineering Problems, 11(5), 1207–1218.
[62]. Cotrina-Teatino, M., Marquina-Araujo, J., & Riquelme, Á. (2025). Comparison of machine learning techniques for mineral resource categorization in a copper deposit in Peru. Natural Resources Research.
[63]. Sinclair, A., & Blackwell, G. (2000). Resource/reserve classification and the qualified person. CIM Bulletin, 93(1038), 29–35.
[64]. Addo, E., Chanda, E., & Metcalfe, A. (2017). Estimation of direction of increase of gold mineralisation using pair-copulas. In: Proceedings of the 22nd International Congress on Modelling and Simulation (MODSIM 2017), Hobart, Australia, 1541–1547.
[65]. Journel, A. (1986). Geostatistics: Models and tools for the earth sciences. Mathematical Geology, 18(1), 119–140.
[66]. Bevilacqua, M., Alvarado, E., & Caamaño-Carrillo, C. (2024). A flexible Clayton-like spatial copula with application to bounded support data. Journal of Multivariate Analysis, 201, 105277.
[67]. Wackernagel, H. (2003). Ordinary kriging. In: Multivariate Geostatistics (pp. 79–88). Springer, Berlin, Heidelberg.
[68]. Daya, A., & Bejari, H. (2015). A comparative study between simple kriging and ordinary kriging for estimating and modeling the Cu concentration in Chehlkureh deposit, SE Iran. Arabian Journal of Geosciences, 8(9), 6003–6020.
[69]. Novia, N. (2022). Estimation of ordinary kriging method with Jackknife technique on rainfall data in Malang Raya. International Journal on Information and Communication Technology (IJoICT), 8(2), 97–106.
[70]. Lamamra, A., Neguritsa, D., & Mazari, M. (2019). Geostatistical modeling by the ordinary kriging in the estimation of mineral resources on the Kieselguhr mine, Algeria. IOP Conference Series: Earth and Environmental Science, 362(1), 012057.
[71]. Klugman, S. (2011). Copula regression. Variance, 5(1), 60–68.
[72]. Bárdossy, A., & Li, J. (2008). Geostatistical interpolation using copulas. Water Resources Research, 44(7), W07412.
[73]. Bárdossy, A., & Hörning, S. (2023). Definition of spatial copula-based dependence using a family of non-Gaussian spatial random fields. Water Resources Research, 59(7), e2022WR033709.
)