Indicator favorability theory for mineral potential mapping

Favorability methods produce a unique measure for mineral potential mapping and quantitative estimation of mineral resources. Indicator favorability theory is developed in this study to account for spatial (auto and cross) correlations of regionalized geological, geochemical, and geophysical fields based on the indicator concept. Target and explanatory indicators are introduced to describe, respectively, direct and indirect evidence of the mineralization of interest. Mineralization is represented by a combination () of a set of target indicators. Indicator favorability theory estimates a regionalized favorability function in two stages: (1) estimate a linear combination of target indicators by maximizing var() and (2) estimate favorability functionF by minimizing estimation variance var[F–]. The model is established on the basis of a conceptual model of target. The favorability estimates can be justified by correlation analysis and cross validation in control areas. The indicator favorability theory is demonstrated on a case study for gold-silver mineral potential mapping based on geophysical, structural, and geochemical fields.     

Application of the Geo-Anomaly Unit Concept in Quantitative Delineation and Assessment of Gold Ore Targets in Western Shandong Uplift Terrain,Eastern China

A number of large and giant ore deposits have been discovered within the relatively small areas of lithospheric structure anomalies, including various boundary zones of tectonic plates. The regions have become the well-known intercontinental ore-forming belts, such as the circum-Pacific gold–copper, copper–molybdenum, and tungsten–tin metallogenic belts. These belts are typical geological anomalous areas. An investigation into the hydrothermal ore deposits in different regions in the former Soviet Union illustrated that the geologic structures of ore fields of almost all major commercial deposits have distinct features compared with the neighboring areas. These areas with distinct features are defined as geo-anomalies. A geo-anomaly refers to such a geologic body or a combination of bodies that their composition, texture–structure, and genesis are significantly different from those of their surroundings. A geo-anomaly unit (GU) is an area containing distinct features that can be delineated with integrated ore-forming information using computer techniques on the basis of the geo-anomaly concept. Herein, the GU concept is illustrated by a case study of delineating the gold ore targets in the western Shandong uplift terrain, eastern China. It includes: (1) analyses of gold ore-forming factors; (2) compilation of normalized regional geochemical map and extraction of geochemical anomalies; (3) compilation of gravitational and aeromagnetic tectonic skeleton map and extraction of gravitational and aeromagnetic anomalies; (4) extraction of circular and linear anomalies from remote-sensing Landsat TM images; (5) establishment of a geo-anomaly conceptual model associated with known gold mineralization; (6) establishment of gold ore-forming favorability by computing techniques; and (7) delineation and assessment of ore-forming units. The units with high favorability are suggested as ore targets.     

Lithodiversity and Its Spatial Association with Metallic Mineral Sites,Great Basin of Nevada

Geographical information system (GIS) techniques were used to investigate the spatial association between metallic mineral sites and lithodiversity in Nevada. Mineral site data sets include various size and type subsets of about 5,500 metal-bearing occurrences and deposits. Lithodiversity was calculated by counting the number of unique geological map units within four sizes of square-shaped sample neighborhoods (2.5-by-2.5, 5-by-5, 10-by-10, and 20-by-20 km) on three different scales of geological maps (national, 1:2,500,000; state, 1:500,000; county, 1:250,000). The spatial association between mineral sites and lithodiversity was observed to increase with increasing lithodiversity. This relationship is consistent for (1) both basin-range and range-only regions, (2) four sizes of sample neighborhoods, (3) various mineral site subsets, (4) the three scales of geological maps, and (5) areas not covered by large-scale maps. A map scale of 1:500,000 and lithodiversity sampling neighborhood of 5-by-5 km was determined to best describe the association. Positive associations occurred for areas having >3 geological map units per neighborhood, with the strongest observed at approximately >7 units. Areas in Nevada with more than three geological map units per 5-by-5 km neighborhood contain more mineral sites than would be expected resulting from chance. High lithodiversity likely reflects the occurrence of complex structural, stratigraphic, and intrusive relationships that are thought to control, focus, localize, or expose mineralization. The application of lithodiversity measurements to areas that are not well explored may help delineate regional-scale exploration targets and provide GIS-supported mineral resource assessment and exploration activity another method that makes use of widely available geological map data.     

Spatial Analysis Techniques as Successful Mineral-Potential Mapping Tools for Orogenic Gold Deposits in the Northern Fennoscandian Shield,Finland

Geoscientific Information Systems (GIS) provide tools to quantitatively analyze and integrate spatially referenced information from geological, geophysical, and geochemical surveys for decision-making processes. Excellent coverage of well-documented, precise and good quality data enables testing of variable exploration models in an efficient and cost effective way with GIS tools. Digital geoscientific data from the Geological Survey of Finland (GTK) are being used widely as spatial evidence in exploration targeting, that is ranking areas based on their exploration importance. In the last few years, spatial analysis techniques including weights-of-evidence, logistic regression, and fuzzy logic, have been increasingly used in GTK’s mineral exploration and geological mapping projects. Special emphasis has been put into the exploration for gold because of the excellent data coverage within the prospective volcanic belts and because of the increased activity in gold exploration in Finland during recent years. In this paper, we describe some successful case histories of using the weights-of-evidence method for the Au-potential mapping. These projects have shown that, by using spatial modeling techniques, exploration targets can be generated by quantitatively analyzing extensive amounts of data from various sources and to rank these target areas based on their exploration potential.     

An Integrated Mapping Approach to Petroleum Exploration on the Pratt Anticline,Pratt County,Kansas

A quantitative map comparison/integration technique to aid in petroleum exploration was applied to an area in south-central Kansas. The visual comparison and integration of maps has become increasingly difficult with the large number and different types of maps necessary to interpret the geology and assess the petroleum potential of an area; therefore, it is desirable to quantify these relationships. The algebraic algorithm used in this application is based on a point-by-point comparison of any number and type of spatial data represented in map form. Ten geological and geophysical maps were compared and integrated, utilizing data from 900 wells located in a nine-township area on the Pratt Anticline in Pratt County, Kansas. Five structure maps, including top of the Lansing Group (Pennsylvanian), Mississippian chert, Mississippian limestone, Viola Limestone (Ordovician), and Arbuckle Group (Cambro-Ordovician), two isopachous maps from top of Mississippian chert to Viola and Lansing to Arbuckle, a Mississippian chert porosity map, Bouguer gravity map, and an aeromagnetic map were processed and interpreted. Before processing, each map was standardized and assigned a relative degree of importance, depending on knowledge of the geology of the area. Once a combination of weights was obtained that most closely resembled the pattern of proved oil fields (target map), a favorability map was constructed based on a coincidence of similarity values and of geological properties of petroleum reservoirs. The resulting favorability maps for the study area indicate location of likely Mississippian chert and lower Paleozoic production.     

Large-scale mineral potential estimation for blind precious metal ore bodies

Conventional evaluation of quantitative mineral potential has focused on target selection at small scales. Mapping at small scales usually results in large-area targets, which may be suitable for grass-roots exploration or regional evaluation of potential. Unfortunately, the estimates in small-scale exploration are commonly associated with large uncertainties. Large-scale estimation is used for optimal in-fill drilling design and step-out drilling target selection. In-fill drilling helps to confirm ore-grade continuities and translate a portion of geological resources into minable reserves, whereas step-out target estimation is useful for finding new orebodies in the vicinity of known ore deposits. Both of these processes are necessary for mine development and production planning. A comprehensive methodology is proposed here, particularly for large-scale mineral exploration. The central information synthesizer is canonical or indicator favorability analysis. A case study is presented to demonstrate the methodology for large-scale target selection. The study involves a gold-mining district where step-out drilling targets are being sought to expand the resource base. Several drilling targets were delineated in the study region. Two of them were tested through surface sampling with positive results.     

Combining AHP with GIS for Predictive Cu Porphyry Potential Mapping: A Case Study in Ahar Area (NW,Iran)

Using the analytic hierarchy process (AHP) method for multi-index evaluation has special advantages, while the use of geographic information systems (GIS) is suitable for spatial analysis. Combining AHP with GIS provides an effective approach for studies of mineral potential mapping evaluation. Selection of potential areas for exploration is a complex process in which many diverse criteria are to be considered. In this article, AHP and GIS are used for providing potential maps for Cu porphyry mineralization on the basis of criteria derived from geologic, geochemical, and geophysical, and remote sensing data including alteration and faults. Each criterion was evaluated with the aid of AHP and the result mapped by GIS. This approach allows the use of a mixture of quantitative and qualitative information for decision-making. The results of application in this article provide acceptable outcomes for copper porphyry exploration.     

Spatial Modeling for Base-Metal Mineral Exploration Through Integration of Geological Data Sets

This study involves the integration of information interpreted from data sets such as LandsatTM, Airborne magnetic, geochemical, geological, and ground-based data of Rajpura—Dariba,Rajasthan, India through GIS with the help of (1) Bayesian statistics based on the weights ofevidence method and (2) a fuzzy logic algorithm to derive spatial models to target potentialbase-metal mineralized areas for future exploration. Of the 24 layers considered, five layers(graphite mica schist (GMS), calc-silicate marble (CALC), NE-SW lineament 0–2000 mcorridor (L4-NESW), Cu 200–250 ppm, and Pb 200–250 ppm) have been identified from theBayesian approach on the basis of contrast. Thus, unique conditions were formed based onthe presence and absence of these five map patterns, which are converted to estimate posteriorprobabilities. The final map, based on the same data used to determine the relationships, showsfour classes of potential zones of sulfide mineralization on the basis of posterior probability.In the fuzzy set approach, membership functions of the layers such as CALC, GMS, NE-SWlineament corridor maps, Pb, and Cu geochemical maps have been integrated to obtain thefinal potential map showing four classes of favorability index.     

GIS-Based Slope Stability Analysis,Chuquicamata Open Pit Copper Mine,Chile

The risk of slope failure in the Chuquicamata open-pit mine was analyzed using Geographic Information System (GIS) software and modeling techniques. Models incorporated various component layers at a relatively large map scale (1:5000): alteration, geotechnical unit, proximity to major faults (VIF), GSI (geological strength index), slope (from digital elevation model), proximity to watertable (difference grid between topography and modeled watertable), and composite structural density grid (VIF, smaller faults, and fracture frequency); not all layers were used in all models. Three modeling techniques were used: fuzzy logic, in which parameters in each component layer were ranked by mine geotechnical experts according to their influence in promoting slope failure, and two data-driven techniques, weights-of-evidence and logistic regression, in which statistical correlation of training points (known failures) with parameters were used to derive a relative probability of failure. Because most slope failures are controlled by structure, VIF and smaller faults were divided by orientation into subsets with dip direction parallel, opposite, and normal to slope aspect; these orientations promote circular and planar, toppling, and wedge-type failures, respectively. Density grids of these subsets show high-risk areas for individual failure types. The models demonstrate sensitivity of the analysis to (1) selection of component layers, (2) selection of training points, (3) classification and ranking of categorical parameters, and (4) data problems in certain layers. Predicted high-risk zones in the final models show a high degree of correspondence with recent, post-model failures. Such models can be used to anticipate future pit design concerns. The results presented here illustrate how vast amounts of data, in multiple geo-referenced layers, can be analyzed and modeled using GIS techniques for predictive studies at relatively large map scales. Such modeling techniques could provide a powerful tool for predictive modeling in a vast array of large-map-scale applications requiring similar data integration and evaluation.     

Use of GIS for a spatial and temporal analysis of Kenyan wildlife with generalised linear modelling

This paper applies generalised linear statistical techniques in a GIS to analyse wildlife data from a Kenyan wildlife reserve and its surrounding areas. Attention focuses on the spatial distribution of elephant during nine successive surveys, analysing their temporal and spatial relationship to 12 environmental covariates. A principal component analysis identifies five major determining factors, thereby reducing dimensionality in the data, while a simple spatial analysis procedure, suitable for wildlife data obtained from airborne surveys, quantfies clustering for different animal species. The number of explanatory variables appearing in abundance models is found to be subject to large variations during successive surveys with a minimum and maximum of four and eight variables, respectively. Species from highly clustered populations are found to have over 20 times more observations within short distances compared to the rest. The study concludes that a combination of generalised linear modelling and GIS gives deeper insight into the dynamics of wildlife species in and around well-defined nature reserves.     

Mineral-Potential Mapping for MVT Deposits with Limited Data Sets Using Landsat Data and Geological Evidence in the Borden Basin, Northern Baffin Island, Nunavut, Canada

Posterior probabilities of occurrence for Zn-Pb Mississippi Valley Type (MVT) mineralization were calculated based on evidence maps derived from regional geology, Landsat-TM, RADARSAT-1, a digital elevation model and aeromagnetic data sets in the Borden Basin of northern Baffin Island, Canada. The vector representation of geological contacts and fault traces were refined according to their characteristics identified in Landsat-TM, RADARSAT-1, DEM, slope, aspect, and shaded relief data layers. Within the study area, there is an association between the occurrence of MVT mineralization and proximity to the contact of platformal carbonates and shale units of the adjacent geological formation. A spatial association also tends to exist between mineralization and proximity to E-W and NW-SE trending faults. The relationships of known MVT occurrences with the geological features were investigated by spatial statistical techniques to generate evidence maps. Supervised classification and filtering were applied to Landsat-TM data to divide the Society Cliffs Formation into major stratigraphic subunits. Because iron oxides have been observed at some of the MVT occurrences within the Borden Basin, Landsat-TM data band ratio (3/1) was calculated to highlight the potential presence of iron-oxides as another evidence map. Processed Landsat-TM data and other derived geological evidence maps provided useful indicators for identifying areas of potential MVT mineralization. Weights of evidence and logistic regression were used independently to integrate and generate posterior probability maps showing areas of potential mineralization based on all derived evidence maps. Results indicate that in spite of the lack of important data sets such as stream or lake sediment geochemistry, Landsat-TM data and regional geological data can be useful for MVT mineral-potential mapping.     

Measuring the Performance of Mineral-Potential Maps

D. P. Harris and others have proposed a new method for comparative analysis of favorability mappings. In their approach, Weights-of-Evidence (WofE) consistently shows poorer results than other more flexible methods. Information loss because of discretization would be a second drawback of WofE. In this paper, we discuss that the random cell selection method proposed by Harris and others necessarily results in higher success ratios for more flexible methods but this does not necessarily indicate that these methods provide better mineral-potential maps. For example, a good point density contouring method that does not use any geoscience background information also would score high in the random cell selection approach. Additionally, we show that discretization usually is advantageous because it prevents occurrences of overly high posterior probabilities. For more detailed comparison, we have conducted a number of experiments on 90 gold deposits in the Gowganda Area of the Canadian Shield comparing WofE with the more flexible weighted logistic regression method. Mineral occurrences should be modeled as discoveries at points instead of randomly sampling them together with their surrounding environments in small cells.     

Exploring Spatially Non-stationary Relationships in the Determinants of Mineralization in 3D Geological Space

Exploring the spatial relationships between various geological features and mineralization is not only conducive to understanding the genesis of ore deposits but can also help to guide mineral exploration by providing predictive mineral maps. However, most current methods assume spatially constant determinants of mineralization and therefore have limited applicability to detecting possible spatially non-stationary relationships between the geological features and the mineralization. In this paper, the spatial variation between the distribution of mineralization and its determining factors is described for a case study in the Dingjiashan Pb–Zn deposit, China. A local regression modeling technique, geological weighted regression (GWR), was leveraged to study the spatial non-stationarity in the 3D geological space. First, ordinary least-squares (OLS) regression was applied, the redundancy and significance of the controlling factors were tested, and the spatial dependency in Zn and Pb ore grade measurements was confirmed. Second, GWR models with different kernel functions in 3D space were applied, and their results were compared to the OLS model. The results show a superior performance of GWR compared with OLS and a significant spatial non-stationarity in the determinants of ore grade. Third, a non-stationarity test was performed. The stationarity index and the Monte Carlo stationarity test demonstrate the non-stationarity of all the variables throughout the area. Finally, the influences of the degree of non-stationary of all controlling factors on mineralization are discussed. The existence of significant non-stationarity of mineral ore determinants in 3D space opens up an exciting avenue for research into the prediction of underground ore bodies.

     

Structural Analysis for Gold Mineralization Using Remote Sensing and Geochemical Techniques in a GIS Environment: Island of Lesvos,Hellas

Exploration for epithermal Au has been active lately in the Aegean Sea of the eastern Mediterranean Basin, both in the islands of the Quaternary arc and in those of the back-arc region. The purpose of this study was the structural mapping and analysis for a preliminary investigation of possible epithermal gold mineralization, using remotely sensed data and techniques, structural and field data, and geochemical information, for a specific area on the Island of Lesvos. Therefore, Landsat-TM and SPOT-Pan satellite images and the Digital Elevation Model (DEM) of the study area were processed digitally using spatial filtering techniques for the enhancement and recognition of the geologically significant lineaments, as well as algebraic operations with band ratios and Principal Component Analysis (PCA), for the identification of alteration zones. Statistical rose diagrams and a SCHMIDT projection Stereo Net were generated from the lineament maps and the collected field data (dip and strike measurements of faults, joints, and veins), respectively. The derived lineament map and the band ratio images were manipulated in a GIS environment, in order to study the relation of the tectonic pattern to both the alteration zoning and the geomorphology of the volcanic field of the study area. Target areas of high interest for possible mineralization also were specified using geochemical techniques, such as X-Ray Diffraction (XRD) analysis, trace-element, and fluid-inclusion analysis. Finally, preliminary conclusions were derived about possible mineralization, the type (high or low sulfidation), and the extent of mineralization, by combining the structural information with geochemical information.     

Mineral Favorability Mapping: A Comparison of Artificial Neural Networks,Logistic Regression,and Discriminant Analysis

A Probabilistic Neural Network (PNN) was trained to classify mineralized and nonmineralized cells using eight geological, geochemical, and geophysical variables. When applied to a second (validation) set of well-explored cells that had been excluded from the training set, the trained PNN generalized well, giving true positive percentages of 86.7 and 93.8 for the mineralized and nonmineralized cells, respectively. All artifical neural networks and statistical models were analyzed and compared by the percentages of mineralized cells and barren cells that would be retained and rejected correctly respectively, for specified cutoff probabilities for mineralization. For example, a cutoff probability for mineralization of 0.5 applied to the PNN probabilities would have retained correctly 87.66% of the mineralized cells and correctly rejected 93.25% of the barren cells of the validation set. Nonparametric discriminant analysis, based upon the Epanechnikov Kernel performed better than logistic regression or parametric discriminant analysis. Moreover, it generalized well to the validation set of well-explored cells, particularly to those cells that were mineralized. However, PNN performed better overall than nonparametric discriminant analysis in that it achieved higher percentages of correct retention and correct rejection of mineralized and barren cells, respectively. Although the generalized regression neural network (GRNN) is not designed for a binary—presence or absence of mineralization— dependent variable, it also performed well in mapping favorability by an index valued on the interval [0, 1]. However, PNN outperformed GRNN in correctly retaining mineralized cells and rejecting barren cells of the validation set.     

Application of GIS Processing Techniques for Producing Mineral Prospectivity Maps—A Case Study: Mesothermal Au in the Swayze Greenstone Belt,Ontario, Canada

A Geographic Information System (GIS) is used to prepare and process digital geoscience data in a variety of ways for producing gold prospectivity maps of the Swayze greenstone belt, Ontario, Canada. Data used to produce these maps include geologic, geochemical, geophysical, and remotely sensed (Landsat). A number of modeling methods are used and are grouped into data-driven (weights of evidence, logistic regression) and knowledge-driven (index and Boolean overlay) methods. The weights of evidence (WofE) technique compares the spatial association of known gold prospects with various indicators (evidence maps) of gold mineralization, to derive a set of weights used to produce the final gold prospectivity map. Logistic regression derives statistical information from evidence maps over each known gold prospect and the coefficients derived from regression analysis are used to weight each evidence map. The gold prospectivity map produced from the index overlay process uses a weighting scheme that is derived from input by the geologist, whereas the Boolean method uses equally weighted binary evidence maps.The resultant gold prospectivity maps are somewhat different in this study as the data comprising the evidence maps were processed purposely differently for each modeling method. Several areas of high gold potential, some of which are coincident with known gold prospects, are evident on the gold prospectivity maps produced using all modeling methods. The majority of these occur in mafic rocks within high strain zones, which is typical of many Archean greenstone belts.     

Evaluation of vegetative criteria for near-surface groundwater detection using multispectral mapping and GIS techniques in semi-arid Botswana

The critical need to consider all options in the search for groundwater in semi-arid areas has promoted work on the possible association of near-surface groundwater and vegetation characteristics using a combination of remote-sensing data and geographic information systems (GIS) techniques. Two vegetative criteria (dense woody cover and abundance of deep-rooting species) are identified as being indicative of near-surface groundwater, and their spatial distribution is tested against the location of aquifers in southeast Botswana. Vegetative criteria classes were combined in a GIS environment with the distribution of geomorphic units and bedrock geology to determine the degree of coincidence with assumed or known aquifers. Results indicate that the distribution of dense woody vegetation as mapped from Thematic Mapper imagery has some potential in identifying especially surficial but also bedrock near-surface groundwater sources in mostly naturally vegetated semi-arid areas. Dense woody cover classes tend to select aquifers in topographically higher areas while classes comprising some deep-rooting species tend to select low-lying aquifers such as those occurring in fossil valleys. Deep-rooting species, however, are less successful as a vegetative criterion. Although various technical refinements are suggested, this work shows that vegetative criteria mapping can however be used in conjunction with conventional geological/geophysical techniques to enhance the prospects for groundwater location in relatively undisturbed semi-arid areas.     

Developing a geographic expert system for regional mapping of volcanogenic massive sulfide (VMS) deposit potential

A personal computer-based geographic information system (GIS) is used to develop a geographic expert system (GES) for mapping and evaluating volcanogenic massive sulfide (VMS) deposit potential. The GES consists of an inference network to represent expert knowledge, and a GIS to handle the spatial analysis and mapping. Evidence from input maps is propagated through the inference network, combining information by means of fuzzy logic and Bayesian updating to yield new maps showing evaluation of hypotheses. Maps of evidence and hypotheses are defined on a probability scale between 0 and 1. Evaluation of the final hypothesis results in a mineral potential map, and the various intermediate hypotheses can also be shown in map form.The inference net, with associated parameters for weighting evidence, is based on a VMS deposit model for the Chisel Lake deposit, a producing mine in the Early Protoerzoic Snow Lake greenstone belt of northwest Manitoba. The model is applied to a small area mapped at a scale of 1:15,840. The geological map, showing lithological and alteration units, provides the basic input to the model. Spatial proximity to contacts of various kinds are particularly important. Three types of evidence are considered: stratigraphic, heat source, and alteration. The final product is a map showing the relative favorability for VMS deposits. The model is implemented as aFortran program, interfaced with the GIS. The sensitivity of the model to changes in the parameters is evaluated by comparing predicted areas of elevated potential with the spatial distribution of known VMS occurrences.     

Use of Fuzzy Membership Input Layers to Combine Subjective Geological Knowledge and Empirical Data in a Neural Network Method for Mineral-Potential Mapping

Use of GIS layers, in which the cell values represent fuzzy membership variables, is an effective method of combining subjective geological knowledge with empirical data in a neural network approach to mineral-prospectivity mapping. In this study, multilayer perceptron (MLP), neural networks are used to combine up to 17 regional exploration variables to predict the potential for orogenic gold deposits in the form of prospectivity maps in the Archean Kalgoorlie Terrane of Western Australia. Two types of fuzzy membership layers are used. In the first type of layer, the statistical relationships between known gold deposits and variables in the GIS thematic layer are used to determine fuzzy membership values. For example, GIS layers depicting solid geology and rock-type combinations of categorical data at the nearest lithological boundary for each cell are converted to fuzzy membership layers representing favorable lithologies and favorable lithological boundaries, respectively. This type of fuzzy-membership input is a useful alternative to the 1-of-N coding used for categorical inputs, particularly if there are a large number of classes. Rheological contrast at lithological boundaries is modeled using a second type of fuzzy membership layer, in which the assignment of fuzzy membership value, although based on geological field data, is subjective. The methods used here could be applied to a large range of subjective data (e.g., favorability of tectonic environment, host stratigraphy, or reactivation along major faults) currently used in regional exploration programs, but which normally would not be included as inputs in an empirical neural network approach.