Search Documents

By V Gau, L Overmeyer, M Kesting

Economically efficient mining processes can be significantly assisted by the implementation of sensor technologies for the detection of materials and boundary layers between materials inside the stratum. These sensors are meant to facilitate the machinery operator and illustrate furthermore the prerequisites in assisting tools towards fully automated mining equipment. Therefore, project æSensory Identification of Material Type and Detection of the InterfacesÆ (SIMT) aims at the integration of ground penetrating radar (GPR) and geoelectric sensor technology by endowing with them a bucket of a bucket wheel excavator. For this reason, a sensor bucket was developed and tested in the Inden Mine of RWE Power AG. The feasibility of integration and deployment of such sensor technology has been proven effective for the conditions of open cast mining. The main focus of current research lies on developing algorithms to analyse the acquired measurement data from the sensors to automatically detect layers and materials. A computer network receives the GPR and geoelectric measurement data as well as the positioning data of the sensor bucket. The positioning data mainly consists of GPS and encoder information, which is merged with the sensor data by using the UTC time stamp. With the aid of consecutive processes the relevant information will be extracted and provided to the user. The aim of SIMT is to replace material identification, traditionally carried out by the operator, by implementing a continuous operating interpretation system in mining equipment based on GPR and geoelectric sensor technology, resulting in more detailed information regarding the working face during an excavation. This detailed information provides an opportunity for efficient planning and controlling of selective mining of natural deposits, such as the process of lignite mining.

Sep 26, 2011

By J Foo, T Kleinschmidt, B Reynolds, K Kennewell

The consideration of complex supply chains often highlights multiple system bottlenecks that are localised and dynamic – they shift depending on what variability last influenced the system. Using the understanding of a dynamic system model, synergies between items that tackle localised bottlenecks can be understood to generate a suite of solutions. These can be simple solutions, such as a couple of new conveyors placed in series, or more complex ones, including buffer dynamics (capacity, location, behaviour, etc).Each solution in the suite will have an associated uplift in system performance and an associated cost of implementation. Dynamic simulations can address each situation in turn and estimate the projected improvement it would offer. However, these solutions will have different benefits when implemented together, and each solution will give different uplift depending on what other solutions have already been assumed. If the number of solutions is numerous, traditional dynamic simulation analysis becomes the bottleneck, limiting options analysis to the processing power available.This paper will illustrate the use of a new technique applied as part of a capital growth optimisation project with a major iron ore producer. Metamodelling is the creation of a model of a model; in this case study, the focus is on developing a response surface representing the suite of solutions. This can be achieved through a combination of dynamic modelling and an appropriate mathematical process called design of experiments. The result is an approximation surface of the outcomes of the dynamic model, encapsulating interactions between key elements of the value chain.This modelling approach provides a number of advantages compared to traditional dynamic modelling approaches, including but not limited to:early and rapid evaluation of a large range of potential solutionsaccommodation of multiple initial states of the value chainoptimisation of the overall suite of upgrades (down to route-level upgrades) based on business objectiveseasy re-evaluation of outcomes when key inputs such as costs change.CITATION:Kleinschmidt, T, Foo, J, Reynolds, B and Kennewell, K, 2015. Case study of iron ore value chain optimisation, in Proceedings Iron Ore 2015, pp 605–612 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 13, 2015

By J Yamatomi, Y Tamanishi, S Miwa

The four-dimensional (4D) network relaxation method, originally developed by Akaike in 1999, uses Lagrange relaxation and is known for its ability to include the price change and possible better outputs, but the computation time is huge. However, the continuous advancement in computer technologies allows the use of higher-speed computers with more memory capacity. Therefore, in this paper the authors have tried to implement the method with adjustments and improvements in its algorithm to suit the current computing system. As a result, the computation time has been drastically reduced and more constraints can be included, which makes the scheduling a more realistic one. The case study of an open pit gold mine with 129 500 blocks can be solved considering the capacity constraints in 71 seconds. In previous case studies, the mined ore had to be processed or disposed of immediately. Now, stockpiling can be included as another option. Since the details of the algorithm are still being improved, the computation time is a little longer than the previous cases without the stockpiling option. However, it is found that scheduling with the stockpiling option is possible in a reasonable time. In conclusion, the 4D network relaxation method has become more useful with advances in computer technologies and provides better results.

Sep 26, 2011

By T Roesch, M Neurohr, S Stratton

"Dwindling commodity prices are forcing today’s mining operations to increase output. Sometimes the solution is to upgrade the existing infrastructure; however, in most cases it has to be achieved with the same or reduced workforce and with the same equipment capacity through improved process efficiency and higher uptime.Measure, Manage and Improveis the mantra for this years’ Mill Operators’ Conference and while this in most cases relates to the above process efficiency, Siemens has set its focus on the availability and uptime aspect since unscheduled downtime is one of the biggest threats to any grinding operation. Pinion misalignments, resonance, unbalance and bearing and gear meshing faults in geared trains, as well as long-term changes in the condition of the mechanical components or structures and auxiliary systems in gearless systems are the main reasons for unscheduled downtimes in today’s mill drive trains. Our paper outlines the efforts that have gone into adapting conventional condition monitoring methods in conventional, pinion driven mill drive trains as well as gearless systems to capture the gradual deterioration of wearing parts or structural changes in the drive train not only on fixed speed application but also on variable speed applications with multiple operating points. The paper further explores the use of condition monitoring systems as a tool to assess spare thermal motor capacity to potentially push equipment beyond its nameplate rating, along with the process that was followed to correlate measured frequency spectra and envelopes to ambient and load data to understand the kinematics of the drive train.Citation:Neurohr, M, Stratton, S and Roesch, T, 2016. Advancements in Condition Monitoring on Mill Drive Trains for Proactive Maintenance and Early Fault Detection Purposes, in Proceedings 13th AusIMM Mill Operators’ Conference 2016, pp 263–272 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Oct 10, 2016

By M C. Schodlok

Banded iron-formation (BIF) and iron ores were measured using the CSIRO developed TIRHyLogger ™ to assess its ability to map quartz, carbonates and riebeckite. The TIR-HyLogger™ measures continuously spectral refl ectance data from diamond drill core or drill chips with a high spatial and spectral resolution. Samples from twelve BIF, bedded iron deposits (BID) and channel iron deposits (CID) have been analysed. The spectral data were processed with two different approaches. The first approach, focused on the identifi cation of the major minerals, was performed using principal component analysis (PCA). It identified four minerals (quartz, ankerite, riebeckite and goethite) and provided pure spectra called end-members. The second approach included a mineral diagnostic spectral feature extraction involving wavelength position, relative peak heights and relative depths of absorptions. Quartz has a diagnostic absorption at 8.625 µm and ankerite shows a refl ectance peak at 11.275 µm. In addition to the identified end-members, two carbonate populations (calcite and siderite) and two silica species (quartz and amorphous silica) were identified.

Jul 11, 2011

By J Rapkoch, M Cerna

As long as the tendency of Chilean copper ores continues to feature high levels of impurities, the smelters will tend to produce high-impurity mattes. The most notable impurities include As, Sb, Pb and Zn. Producing high-quality copper with such materials is a major challenge to metallurgical engineers and operators. More investigations into process optimisations and/or alternatives are needed, using both theoretical knowledge and operational skill. Converting of matte in Chile depends on the matte grade using established technologies. Mattes up to 66 per cent Cu, such as those from Outokumpu Flash furnaces, are converted in two stages, while typical Teniente Converter mattes, ranging from 70 - 74 per cent Cu are converted in one blow, directly to copper. In this paper thermodynamic data are used to evaluate possible variations in the use of traditional technologies applied to smelt copper concentrates with high impurity levels. The results of this analysis help to predict the impurity distributions between matte-slag-gas during the slag blow, and between copper-slag-gas during the copper blow. For converting high-impurity white metal, the effectiveness of volatilisation in the slag blow is analysed. In the copper blow, the type and the amount of flux (SiO2, CaO and others) are altered with the aim of optimising impurity distribution to the slag and gas.

Jan 1, 2005

By D R. L Vedoy, L Botha, S P. Gumfekar, J B. P Soares

Most of the known in-place reserves of oil sands are found in Canada. These oil sands ores are mined using various processes to extract bitumen. The by-product slurry resulting from the extraction process, known as tailings, is discharged into tailings ponds. Fine kaolin clays are an important ingredient of the tailings because they remain suspended in water for many years. Current polymeric flocculation/dewatering technologies use acrylamide-based polymers such as polyacrylamide (PAM), but sediments produced using PAM-based flocculants hold considerable amounts of water due to the formation of loosely-packed flocs. Anionic polymers are better flocculants due to bridging and the charge-neutralisation mechanism. In this paper, we will compare the synthesis of two flocculants: thermosensitive, anionic flocculants and thermosensitive, hydrophobically modified anionic flocculants. The first flocculant is a copolymer of acrylic acid (AA) and N-isopropyl acrylamide (NIPAM), in which AA is the anionic monomer. The second flocculant is a copolymer of NIPAM, AA and the hydrophobic monomer N-tert-butyl acrylamide. We evaluated the dewatering efficiency of 2 wt per cent kaolin using both flocculants by comparing settling rates at different flocculant dosages, charge neutralisation ability and supernatant clarity. To better understand floc formation, we measured the floc size formed by both flocculants at different flocculant dosages. We propose that the hydrophobic monomer in the terpolymer flocculant helps expel water out of flocs, ultimately decreasing the water content in the sediments.CITATION:Gumfekar, S P, Vedoy, D R L, Botha, L and Soares, J B P, 2015. Thermosensitive, hydrophobically modified anionic flocculants for dewatering kaolin clays, in Proceedings Tailings and Mine Waste Management for the 21st Century , pp 45–50 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 27, 2015

By G McKay, M Erskine, S Macklin

The City of Sydney is currently investing in a number of infrastructure projects that require tunnel excavations within the Triassic Hawkesbury Sandstone. It is generally accepted that this formation has high residual stresses of ‘tectonic’ origin, so much so that these stresses have been implicated in a number of well-known failures and problems for tunnelling. These include failure of ‘roof beams’ comprising discrete beds in the crown of tunnel excavations; stress-induced spalling and ‘slabbing’ of tunnel walls; shear displacement between beds causing the failure of rock bolts; strain in adjacent infrastructure and, the necessity for re-grouting of high conductivity fractures due to stress relief dilation after excavation. Whilst there is an extensive published database of in situ stress measurements comprising ‘hydrofrac’, ‘over-coring’ and ‘slotter’ stress measurements, these data indicate a considerable spread of results. The traditional approach to interpreting this data has been to attempt to fit linear correlations to stress versus depth data, based on personal experience and ‘back-analysis’ of monitoring data during construction. However, there are a number of such correlations in the literature and there is divergence of opinion regarding the assessment of ‘background’ in situ stress and appropriate upper and lower bound values for design. This paper provides a brief review of the published data; suggestions on an alternative ‘risk-based approach’ to assigning in situ stresses for design; and, proposals on how in situ stress tests may be conducted in future to improve the quality of the database available to designers.CITATION:Macklin, S, McKay, G and Erskine, M, 2014. A review of the in situ stress database for tunnel design in the Hawkesbury Sandstone, Sydney, New South Wales, in Proceedings 15th Australasian Tunnelling Conference 2014 , pp 171–178 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Sep 17, 2014

By S E. Box, J L. Mauk

The Middle Proterozoic metasedimentary rocks of the Belt Supergroup contain important metallic mineral deposits like the diagenetic sediment-hosted stratiform copper mineralization, including the Troy (Spar Lake) deposit and the epigenetic hydrothermal veins in the lower Belt and Ravalli Group in the Coeur d'Alene district. Magnetite (Fe3O4) is a common and widespread accessory mineral in these rocks and can have multiple origins including ore stage in deposits, mobilization stage in source rocks, metamorphism of red beds and alteration mineral in thin halos enveloping Ag-rich veins in the Coeur dÆAlene district. As shown in several classic studies, the composition of magnetite from igneous rocks varies significantly in response to controlling factors such as melt composition, temperature, silica activity, and oxygen and sulphur fugacity. Therefore, the distinctively different geologic settings within the Belt Supergroup strongly suggest that there may be differences in the magnetite compositions. However, initial microprobe analyses detected no variation in magnetite composition.   Three possibilities may give an adequate explanation for these initial results: 1.      the magnetite formed primarily as pure Fe3O4, 2.      one or more tectonothermal events led to regional scale re-equilibration that affected all rocks in the field area, and 3.      the variations in the composition are below the detection limit of the microprobe.

Jan 1, 2007

By L Y. Yang, Z Wang

This work presents the fundamental properties of CITIC magnetite concentrate and its performance in pelletising and sintering. CITIC concentrate has been widely used in pelletising to produce pellets for blast furnaces. It is also used in sintering to benefit from its sintering performance. CITIC concentrate has excellent balling properties owing to its size distribution, particle surface properties, large specific surface area and balling index, leading to strong green pellets and product pellets. When used in sintering at appropriate levels, sinter bed permeability improves, productivity is increased while coke consumption reduces, and sinter strength improves.

Nov 8, 2021

Most of the Ta-Nb ores in China occur in altered granites, and Ta-Nb minerals occur in these granites as accessory minerals. This type of Ta-Nb ores are characterised by whole mineralised granite mass, large reserves, low grade of tantalum and niobium, less economic minerals, complex mineralogy, difficulty in separation and concentration of Ta-Nb minerals and high operating cost of ore dressing. The Ta-Nb ore presented in this paper is a low-grade Ta-Nb ore occurring in altered granite in Jiangxi province. It has large reserves, with the grade of Ta2O5 0.013 per cent and Nb2O5 0.023 per cent. The process mineralogy study carried out by the authors determined ore mineralogy and ascertained occurrence of Nb and Ta and their mineral processing characteristics. The results are used as the guide for developing the optimum mineral processing flow sheet. The results also indicated that this type of Ta-Nb ore contains a high content of Li, K and Na, a low content of deleterious impurity Fe, and the tailings can be used as feed for glass and ceramic materials. This study provided important mineralogy information for exploiting and utilising low-grade Ta-Nb ore resources, developing a mineral processing flow sheet and further utilising the tailings.

Jan 1, 2008

By H Arvidson, V Chamberlain, R Marinho, B Parsons, T Rowl, M Noppé, M Mattera

With the introduction of CRIRSCO (2019) based mineral asset disclosure, Regulation S-K part 1300, 2019 (S-K1300) by the United States (US) Securities and Exchange Commission (SEC), applicable to companies listed on the New York Stock Exchange from 2021, Mineral Resource is now required to be reported Exclusive of Mineral Reserve, referred to here as ‘Exclusive Mineral Resource’ (EMR). EMR is equivalent to reporting Mineral Resource additional to Mineral Reserve, which has been permitted in jurisdictions outside the US for many years. Despite this, there are no practical guidelines in the literature or relevant reporting codes and rules for Competent/Qualified Persons to refer to, which creates a risk that disclosure of EMR estimates across the industry will be inconsistent and potentially ambiguous for analysts, investors, potential investors, and other stakeholders.With the introduction of CRIRSCO (2019) based mineral asset disclosure, Regulation S-K part 1300, 2019 (S-K1300) by the United States (US) Securities and Exchange Commission (SEC), applicable to companies listed on the New York Stock Exchange from 2021, Mineral Resource is now required to be reported Exclusive of Mineral Reserve, referred to here as ‘Exclusive Mineral Resource’ (EMR). EMR is equivalent to reporting Mineral Resource additional to Mineral Reserve, which has been permitted in jurisdictions outside the US for many years. Despite this, there are no practical guidelines in the literature or relevant reporting codes and rules for Competent/Qualified Persons to refer to, which creates a risk that disclosure of EMR estimates across the industry will be inconsistent and potentially ambiguous for analysts, investors, potential investors, and other stakeholders.

May 24, 2023

By P T. Stoker, A Allman, D P. Carville

Classification of mineral resource estimates is one of the most important responsibilities of the Competent Person. In recent years, classification decisions have been driven more by the detailed block-by-block attributes generated by the now widely applied geostatistical estimation methods, and less by a general geological overview. This is increasingly resulting in æspotted dogÆ outputs, in which blocks of Inferred Resources or unclassified material separate blocks of Measured and/or Indicated Resources, or individual drill holes are surrounded by annuli of Measured and Indicated Resource blocks. Not only are such outputs potentially misleading, since they ignore fundamentals such as continuity of geology and mineralisation between drill holes and the imprecise nature of resource estimation, they can also cause substantial problems for engineers undertaking mine designs and estimating ore reserves, particularly for underground mines. Competent Persons for mineral resources must keep in mind the purpose of their work, and should use their experience and judgement to avoid or smooth out æspotted dogÆ classifications, providing a result commensurate with the level of geological and resource estimation confidence. For their part, ore reserve Competent Persons must understand the resource estimation and classification process, and should question resource classifications that are not consistent with the level of geological and resource estimation confidence.

Jan 1, 2006

By Edmonds W, Giles C. W, Webb R. J, Gum J. C

Gold production from The Granites - Tanami Goldfield during 1993 is expected to exceed 8 t, eclipsing that of the traditional Northern Territory gold mining districts of Tennant Creek and Pine Creek. This is the result of dedicated and successful exploration, and the consequent discovery of new orebodies such as those at Dead Bullock Soak.Following granting of access to regional exploration areas in 1988 and identification of favourable host rocks in the Schist Hills area, North Flinders Mines (NFM) began a methodical exploration program of detailed geological mapping, rock chip sampling, costeaning and various types of drilling. The first resource of 600 000 t of 3.3 g/t gold was announced in mid-1989. This has grown progressively with systematic exploration to the currently published resource inventory of 6.6 mT @ 5.4 g/t, (as at 31.12.92) contained in six separate orebodies within the Dead Bullock Soak Mineral Lease at Schist Hills.

Jan 1, 1993

By A Tawadrous, B Wheeler, D S. Preece, S A. Silling

The two key results of rock blasting are fragmentation and movement of geomaterials. Movement/flow of the blasted rock is often referred to as heave. Moving rock with explosives rather than mechanical equipment can improve mining efficiency. Such is the case with coalmine cast blasting. Controlling and directing rock movement through explosive deck delay timing can greatly improve ore recovery by segregating ore and waste. In both cases, understanding, predicting and controlling heave is an important element of successful blasting.Computer programs for predicting 2D blasting-induced heave have been used in mining for decades, helping to drive the innovation of new and more efficient mining methods, ie cast blasting. The past few years have seen the advent of 3D distinct element modelling of blast-induced heave. This is an important advancement because most blasting situations are three-dimensional in nature, especially when controlling the direction and intensity of rock movement through explosive deck delay timing.This paper documents the first full-scale blasting simulations in 3D with distinct element models that include hundreds of blastholes and in excess of a million particles, made possible with parallel processing. Surface mineral blasting is the focus of this study which treats the effect of blasthole delay timing on ore recovery. Ore waste and dilution resulting from typical row-on-row delay timing will be demonstrated. Also demonstrated will be the ore recovery possible with delay timing designs that significantly enhance ore/waste segregation. This type of high fidelity modelling and the use of electronic detonators makes possible the sculpting of rock blasting muck into separate piles of ore and waste.CITATION:Preece, D S, Tawadrous, A, Silling, S A and Wheeler, B, 2015. Modelling full-scale blast heave with three-dimensional distinct elements and parallel processing, in Proceedings 11th International Symposium on Rock Fragmentation by Blasting, pp 127–136 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 24, 2015

By Davies K. E

Several nuclear techniques have been assessed for application to online conveyor belt analysis of oil yield from Rundle oil shale. The most favourable technique is neutron inelastic scattering which determines the total carbon content of bulk shale samples. Correlations of chemical laboratory assays for 540 Rundle oil shale samples show that oil yield can be predicted from total carbon to within 11 litre/tonne.Laboratory experimental work on 17 bulk samples of Rundle oil shale indicated that carbon and oil yield can be determined to within 0.6 wt per cent and 9 litre/tonne, respectively, using the neutron inelastic scatter technique. A practical neutron inelastic scatter gauge for application on a 23 mm thick reinforced conveyor belt could consist of a 5 curie 241 Am-Be neutron source and a 150 mm diameter x 100 mm thick Nal(TI) detector. With this gauge, the analysis time to achieve a counting statistical error of 6 litre/tonne should be about five minutes.

Jan 1, 1985

By R Pax

Feed ore variability, the major disturbance into a mineral processing plant receives considerable attention operationally, and at research and strategic levels to ensure that the throughput of the ore are maintained at all times whilst encouraging a sustainable approach to recovering the valuable minerals from ore bodies. Effort and insights from mineralogical and processing characterisation (eg. ore hardness) of sub-ore body domains allows sequencing for processing to ensure consistency and throughput for the size sensitive grinding unit processes. Optimising the blast fragmentation has benefit as does pre-concentration using ore sorting. Blending of ore domains with different characteristics is common for operations. In recent years a number of technologies have become available for on-line and more rapid assessment of ore and have been discussed in a recent paper by Cameron et.al. (2016) as an opportunity for further utilisation of mine-to-mill (and mine to product) strategies. Current concepts have focused on extensive characterisation of ores so that current control practices can be realised and made more effective with relatively slow measurements and actuation. CITATION: Pax, R, 2018. Dynamic concepts for mine to mill operation and optimisation, in Proceedings Complex Orebodies Conference 2018, pp 124–125 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 21, 2018

By R J. Higgins, M Edraki

What is the meaning of the Life of a Mine (LOM) if that mine results in a requirement for water treatment in perpetuity? ‘In perpetuity’ is of course a duration with no relationship to a Life-of-Mine or what might rightly be regarded by external stakeholders as a fair and reasonable post-mine recovery period. It exceeds the life span of investors, mine managers and governments, indeed even of capital markets as we know them, and of corporate and individual memories.In mine and mine closure planning, decisions are usually based at least to a significant extent on traditional economic analysis underlain by discounted cash flows (DCF). For investments in income-producing facilities, such as mining development and equipment, and ore-processing plant and equipment, this analysis reflects commercial reality and allows for comparison of competing investments. However, for environmental and social infrastructure, the same analysis using the same input parameters, especially the discount rate applied to cash flows, can lead to a bias against the adoption of best practices for non-investor stakeholders. CITATION:Higgins, R J and Edraki, M, 2018. Water treatment in perpetuity – or not?, in Proceedings Life-of-Mine 2018, pp 131–132 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 25, 2018

Blasthole cone samples are used for grade control in most of the iron ore mines of the Pilbara of Western Australia. This method of sampling is straightforward when the target sampling interval matches the mining bench height but is problematic when the blasted bench is 'flitch' mined in two or more passes. Crucially, the sampler needs to decide which part of the blast cone represents the upper and lower flitches to be mined. This sampling selection decision is usually made by assuming a certain proportion represents the subdrill component of a sampling face in a blasthole cone and the two flitches are split equally from the remaining portion. In this study of a West Australian channel iron deposit, sampling experiments were carried out to determine the actual proportions of a sampling face in a blasthole cone by spray painting the blast cones at flitch and subdrill boundaries during the drilling process. Measurements were then taken to assess the differences between assumed sampling face proportions and actual proportions. A set of paired samples was also collected from a number of blastholes with one set using the assumed proportions of each flitch for each sample, while the other set sampled the actual proportions, which were identified from the cone painting process.

Mar 1, 2010

By P Dowd, A Adeli, C Xu, X Emery

Geometallurgical modelling is increasingly being incorporated into mineral resource modelling and estimation as a means of increasing efficiency, decreasing operating costs and reducing risk in mining operations. Traditional geostatistical workflows for geometallurgical modelling face several difficulties, in particular, related to the non-additivity of process response variables, to highly heterotopic data sets and/or to preferential sampling designs, all of which may explain why most attempts to estimate geometallurgical response variables from assay data and mineralogy are still based on statistics and machine learning. The aim of the work presented here is to demonstrate the applicability of geostatistics to estimate copper recovery from copper sulfide ores at the Prominent Hill Iron Oxide Copper-Gold (IOCG) deposit in Australia. As recovery is non-additive, the preferred workflow is the joint estimation of the masses of metal in the feed and in the concentrate, both of which are additive variables. While the mass of metal in the feed can be estimated from abundant online assay analyses, the mass of metal in the concentrate can only be estimated from a very limited number of laboratory-scale batch flotation tests. Traditional approaches to estimating the masses of metal in the feed and in the concentrate, such as simple and ordinary co-kriging, are described and discussed. A modified version of co-kriging is introduced, which incorporates linear relationships between the mean values of the input variables, in addition to a spatial correlation model. Co-kriging with related means is shown to outperform traditional simple and ordinary co-kriging in terms of precision, accuracy and consistency of the estimates, and offers a practical and efficient tool to deal with highly heterotopic and preferential sampling designs.

May 24, 2023