Framework and modeling structure of Exploration and Mining Cadastre system in Iran

In this paper we introduce a framework for developing the Exploration and Mining Cadastre System (EMCS), which involves all activities related to exploration areas, including registration and subsequent legal issues in a broader sense. It is built with two main purposes in mind: (1) enable necessary surveys, observe applicant’s reserved parcel, register the parcel online, and issue and follow up necessary permits/licenses in the field of exploration, (2) provide a mutual platform for data exchange among various service systems in the field of exploration and mining in Iran. The EMCS model is designed in Unified Modeling Language (UML) with three main sets of classes: (1) RegisterObject (2) ActivityLicense, and (3) Person. Other relevant features of the model are also specified at the class level. These include: conditional registered parcels, unregisterable areas such as prohibited or protected areas, legal borders around registered or serving parcels, reference documents such as legal survey results documents, responsibilities, and restrictions. The model also considers issues related to exploration cost, parcel code, tax, royalty, commodity, priorities, facilitations, evaluation, transfer of ownership, changing parameters, and temporal registration. The temporal aspects of the classes are also supported by the model to take care of the time periods that involve the general procedures in the registration of exploration and mining properties. We anticipate that the framework and model structure provide a platform for data storage, improve the supervision of exploration and mining process, change manual into online registration system, prevent administrative corruption, and reduce bureaucracy.     

Multiple data-type processing as applied to the selection of areas with high gold potential in French Guyana

As part of the Guyana Shield, the basement rocks of French Guyana, between latitudes 3° and 3°30′N, consist of greenstones belts of Proterozoic age (2.0–2.1 Ga) enclosed by anatectic granitoïds and migmatites, and of acid or basic intrusives.     

Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization

This paper proposes a novel history-matching method where reservoir structure is inverted from dynamic fluid flow response. The proposed workflow consists of searching for models that match production history from a large set of prior structural model realizations. This prior set represents the reservoir structural uncertainty because of interpretation uncertainty on seismic sections. To make such a search effective, we introduce a parameter space defined with a “similarity distance” for accommodating this large set of realizations. The inverse solutions are found using a stochastic search method. Realistic reservoir examples are presented to prove the applicability of the proposed method.     

A hybrid AHP-VIKOR approach for prospectivity modeling of porphyry Cu deposits in the Varzaghan District,NW Iran

Identifying highly favorable areas related to a particular mineralization type is the main objective of mineral prospectivity modeling (MPM). The northwestern portion of Ahar-Arasbaran porphyry copper belt (AAPCB) is situated within the Urumieh-Dokhtar magmatic belt (UDMB). Because of owning many worthwhile Cu-Mo and Cu-Au porphyry deposits, this area is entitled to incorporate diverse spatial evidence layers for the MPM. In this paper, a hybrid AHP-VIKOR, as an improved knowledge-driven MPM procedure has been proposed for integration of various exploration evidence layers. For this, the AHP is used to calculate important weights of spatial criteria while the VIKOR is applied to outline ultimate prospectivity model. Six effective spatial evidence layers pertaining to the Varzaghan District are selected: (1) multi-elemental geochemical layer of Cu-Mo-Bi-Au; (2) remotely sensed data of argillic, phyllic, and iron oxide alteration layers; and (3) geological and structural layers of Oligo-Miocene intrusions and fault. In addition, a fuzzy prospectivity model (γ?=?0.9) is implemented to assess the AHP-VIKOR approach. Two credible validation methods comprising normalized density index and success rate curve are adapted for quantitative evaluation of predictive models and enhancing the probability of exploration success. The achieved results proved the higher accuracy of the AHP-VIKOR model compared with the fuzzy model in delimiting the favorable areas.     

A comprehensive VIKOR method for integration of various exploratory data in mineral potential mapping

The central Iranian volcanic-sedimentary belt in Kerman province of Iran that is located within the Urumieh-Dokhtar magmatic arc zone is chosen to integrate diverse evidential layers for mineral potential mapping. The studied area has high potential of mineral occurrences especially porphyry copper, and the prepared potential maps aim to outline new prospect zones for further investigation. Two evidential layers including the downward continued map and the analytic signal of filtered magnetic data are generated to be used as geophysical plausible traces of porphyry copper occurrences. The low values of the resistivity layer acquired from airborne frequency domain electromagnetic data are also used as an electrical criterion in this study. Four remote sensing evidential layers including argillic, phyllic, propylitic, and hydroxyl alterations are extracted from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images in order to map the altered areas associated with porphyry copper deposits. The Enhanced Thematic Mapper Plus (ETM+) images are used as well to prepare iron oxide layer. Since potassic alteration is generally the mainstay of copper ore mineralization, the airborne potassium radiometry data is used to explore both phyllic and potassic alteration. Finally, the geochemical layers of Cu/B/Pb/Zn elements and the main geochemical component responsible for ore mineralization extracted from principal component analysis are included in the integration process to prepare final potential maps. The conventional and the extended version of VIKOR method (as a well-known algorithm in multi-criteria decision making problems) produced two mineral potential maps, and the results were compared with the ones acquired from prevalent methods of the index overlay and fuzzy logic operators of sum and gamma. The final mineral potential maps based upon desired geo-data set indicate adequately matching of high potential zones with previous working and active mines of copper deposits.     

The Anarak, Jandaq and Posht-e-Badam metamorphic complexes in central Iran: New geological data, relationships and tectonic implications

The Anarak, Jandaq and Posht-e-Badam metamorphic complexes occupy the NW part of the Central-East Iranian Microcontinent and are juxtaposed with the Great Kavir block and Sanandaj-Sirjan zone. Our recent findings redefine the origin of these complexes, so far attributed to the Precambrian–Early Paleozoic orogenic episodes, and now directly related to the tectonic evolution of the Paleo-Tethys Ocean. This tectonic evolution was initiated by Late Ordovician–Early Devonian rifting events and terminated in the Triassic by the Eocimmerian collision event due to the docking of the Cimmerian blocks with the Asiatic Turan block.

The “Variscan accretionary complex” is a new name we proposed for the most widely distributed metamorphic rocks connected to the Anarak and Jandaq complexes. This accretionary complex exposed from SW of Jandaq to the Anarak and Kabudan areas is a thick and fine grain siliciclastic sequence accompanied by marginal-sea ophiolitic remnants, including gabbro-basalts with a supra-subduction-geochemical signature. New ⁴⁰Ar/³⁹Ar ages are obtained as 333–320 Ma for the metamorphism of this sequence under greenschist to amphibolite facies. Moreover, the limy intercalations in the volcano-sedimentary part of this complex in Godar-e-Siah yielded Upper Devonian–Tournaisian conodonts. The northeastern part of this complex in the Jandaq area was intruded by 215 ± 15 Ma arc to collisional granite and pegmatites dated by ID-TIMS and its metamorphic rocks are characterized by some ⁴⁰Ar/³⁹Ar radiometric ages of 163–156 Ma.

The “Variscan” accretionary complex was northwardly accreted to the Airekan granitic terrane dated at 549 ± 15 Ma. Later, from the Late Carboniferous to Triassic, huge amounts of oceanic material were accreted to its southern side and penetrated by several seamounts such as the Anarak and Kabudan. This new period of accretion is supported by the 280–230 Ma ⁴⁰Ar/³⁹Ar ages for the Anarak mild high-pressure metamorphic rocks and a 262 Ma U–Pb age for the trondhjemite–rhyolite association of that area. The Triassic Bayazeh flysch filled the foreland basin during the final closure of the Paleo-Tethys Ocean and was partly deposited and/or thrusted onto the Cimmerian Yazd block.

The Paleo-Tethys magmatic arc products have been well-preserved in the Late Devonian–Carboniferous Godar-e-Siah intra-arc deposits and the Triassic Nakhlak fore-arc succession. On the passive margin of the Cimmerian block, in the Yazd region, the nearly continuous Upper Paleozoic platform-type deposition was totally interrupted during the Middle to Late Triassic. Local erosion, down to Lower Paleozoic levels, may be related to flexural bulge erosion. The platform was finally unconformably covered by Liassic continental molassic deposits of the Shemshak.

One of the extensional periods related to Neo-Tethyan back-arc rifting in Late Cretaceous time finally separated parts of the Eocimmerian collisional domain from the Eurasian Turan domain. The opening and closing of this new ocean, characterized by the Nain and Sabzevar ophiolitic mélanges, finally transported the Anarak–Jandaq composite terrane to Central Iran, accompanied by large scale rotation of the Central-East Iranian Microcontinent (CEIM). Due to many similarities between the Posht-e-Badam metamorphic complex and the Anarak–Jandaq composite terrane, the former could be part of the latter, if it was transported further south during Tertiary time.     

Factor analysis applied to regional geochemical data: problems and possibilities

A large regional geochemical data set of C-horizon podzol samples from a 188,000 km² area in the European Arctic, analysed for more than 50 elements, was used to test the influence of different variants of factor analysis on the results extracted. Due to the nature of regional geochemical data (neither normal nor log-normal, strongly skewed, often multi-modal data distributions), the simplest methods of factor analysis with the least statistical assumptions perform best. As a result of this test it can generally be suggested to use principal factor analysis with an orthogonal rotation for such data. Selecting the number of factors to extract is difficult, however, the scree plot provides some useful help. For the test data, a low number of extracted factors gave the most informative results. Deleting or adding just 1 element in the input matrix can drastically change the results of factor analysis. Given that selection of elements is often rather based on availability of analytical packages (or detection limits) than on geochemical reasoning this is a disturbing result. Factor analysis revealed the most interesting data structures when a low number of variables were entered. A graphical presentation of the loadings and a simple, automated mapping technique allows extraction of the most interesting results of different factor analyses in one glance. Results presented here underline the importance of careful univariate data analysis prior to entering factor analysis. Outliers should be removed from the dataset and different populations present in the data should be treated separately. Factor analysis can be used to explore a large data set for hidden multivariate data structures.     

Groundwater modeling applied to production well management in a contaminated aquifer

A study was conducted to evaluate production strategies for a well field system near a source of groundwater contamination. Numerical modeling of groundwater flow was employed to generate hydraulic head configurations for different production scenarios. For a given scenario, an evaluation of contamination susceptibility was made by comparing head distributions in two aquifer units to the positions of the contaminant source and discharging water supply wells. The results of this study suggest that groundwater flow modeling can be a useful technique for planning the production of water supply wells in aquifers at risk of contamination from anthropogenic pollution sources.     

Cluster analysis applied to regional geochemical data: Problems and possibilities

Cluster analysis can be used to group samples and to develop ideas about the multivariate geochemistry of the data set at hand. Due to the complex nature of regional geochemical data (neither normal nor log-normal, strongly skewed, often multi-modal data distributions, data closure), cluster analysis results often strongly depend on the preparation of the data (e.g. choice of the transformation) and on the clustering algorithm selected. Different variants of cluster analysis can lead to surprisingly different cluster centroids, cluster sizes and classifications even when using exactly the same input data. Cluster analysis should not be misused as a statistical “proof” of certain relationships in the data. The use of cluster analysis as an exploratory data analysis tool requires a powerful program system to test different data preparation, processing and clustering methods, including the ability to present the results in a number of easy to grasp graphics. Such a tool has been developed as a package for the R statistical software. Two example data sets from geochemistry are used to demonstrate how the results change with different data preparation and clustering methods. A data set from S-Norway with a known number of clusters and cluster membership is used to test the performance of different clustering and data preparation techniques. For a complex data set from the Kola Peninsula, cluster analysis is applied to explore regional data structures.     

Determination of sand dune characteristics through geomorphometry and wind data analysis in central Iran (Kashan Erg)

The combination of wind measurements and remotely sensed geomorphometry indices provides a valuable resource in the study of desert landforms, because arduous desert environments are difficult to access. In this research, we couple wind data and geomorphometry to separate and classify different sand dunes in Kashan Erg in central Iran. Additionally, the effect of sand-fixing projects on sand dune morphology was assessed using geomorphometry indices (roughness, curvature, surface area, dune spacing and dune height). Results showed that a Digital Elevation Model of the National Cartographic Center of Iran (NCC DEM) with 10-m resolution and accuracy of 54% could discriminate geomorphometry parameters better than the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data with 30-m resolution and Shuttle Radar Topography Mission (SRTM) data with 90-m resolution and 45.2 and 1.6% accuracy, respectively. Low classification of SRTM DEM was associated with too many non-value points found in the DEM. Accuracy assessment of comparison ground control points revealed that ASTER DEM (RMSE = 4.25) has higher accuracy than SRTM and NCC DEMs in this region. Study of curvature showed that transverse and linear sand dunes were formed in concave topography rather than convex. Reduced slopes in fixed sand dunes were established due to wind erosion control projects. Measurements of dune height and spacing show that there is significant correlation in compound dunes (R ² = 0.546), linear dunes (R ² = 0.228) and fixed dunes (R ² = 0.129). In general, the height of dunes in Kashan Erg increases from the margin of the field to the center of the field with a maximum height of 120 m in star dunes. Analysis of wind data showed that sand drift potential is in low-medium class in Kashan Erg. Linear sand dunes in Kashan Erg show that they are following a global trend in forming of these. Finally, established of geomorphometry method in dune classification will help researchers to identify priority of land management and performance assessment of sand dunes fixing projects in arid arduous environment.     

Physicochemical modeling as applied to study of sulfoarsenide complexes in hydrothermal solutions

System As–Na–S–Cl–H–O was studied. The research was carried out in three stages: (1) selection of the most likely complexes resulting from arsenic sulfide dissolution, (2) calculation of their thermodynamic constants, and (3) comparison of calculated data with thermodynamic database obtained in tests with the solution of inverse thermodynamic problems using the Selektor program complex. The system As–Na–S–Cl–H–O included more than 230 dependent components, which were divided into two groups, base and functional. The former group includes components of the solution (NaCl, NaOH, Na₂S, NaHS, HCl, H₂S, H₂SO₄, sulfates, H₂SO₃, sulfites, thiosulfates, Na⁺, Cl⁻,HS⁻, S²⁻), gas phase (43 components), and solid phase (orpiment, red arsenic, arsenolite, claudetite, arsenic, sulfur, sodium salts). Thermodynamic constants of the base components are contained in the Selektor database (they were borrowed from reference-books). The latter group includes 77 complexes labile in the solution but determining the solubility of arsenic and stability of its solid phases. Physicochemical modeling was performed in H₂S (≥0.01 m, pH = 1–10), Na₂S, and NaHS solutions at 25–250 °C and saturated-vapor pressure. It has been established that the dissolution of arsenic sulfide mineral phases in subneutral and alkaline solutions at low oxidation potential is favored by the formation of sulfoarsenides, which are more stable than arsenides and arsenates. Thermodynamic constants of functional complexes determining the orpiment solubility were calculated within the experimental error. It is shown that in hydrothermal iron-free systems with a low oxidation potential, the concentration of As in the solution decreases on cooling and with acidity increase.     

Interpretative matrices approach to ranking lake sub-basin pollution potential: an applied study in Brazil

In 2005 a Cylindrospermopsis raciborskii bloom occurred in the Rio Verde Lake Basin (Brazil). To address this concern, a field analysis was performed to measure physicochemical variables and flows in 14 sub-basins, between 2008 and 2009. Measurements of mean total P (0.039 mg/L ± 0.018 mean SD), mean total Kjeldahl N (0.260 mg/L ± 0.226 mean SD), and mean BOD (1.2 mg/L ± 0.4 mean SD) concentrations were low in most streams, while COD reached a high of 27.1 mg/L (±4.9 mean SD). One tributary was responsible for 85 % of TP load, 77.1 % of TKN load, 78 % of t-BOD load, and 79 % of t-COD load. These concentrations and loads were used to develop the pollution potential assessment matrix (2PAM), which considered three different perspectives: stream water quality, reservoir ecosystem equilibrium, and sub-basin management. Each factor (TP, TKN, BOD and COD) was weighted based on concentration, total load and unit-area load. Pollution potential differed depending on which perspective was considered. The matrix developed, 2PAM, provides a new way to analyze concentrations and loads, enabling basin managers to prioritize action plans according to desired use within the basin.     

Geostatistical modeling of salinity as a basis for irrigation management and crop selection—A case study in central Tunisia

Spatial distribution of the salinity in groundwater used for irrigation and irrigated soil was investigated on Gammouda Plain, an important agricultural area in central Tunisia. Samples from 188 wells and irrigated fields were obtained in November 1988. Groundwater salinity ranges from 2 to 30 dS/m, and 50 percent of the sampled wells have a salinity higher than 3.5 dS/m. In most cases the soil samples have a salinity higher than 4 dS/m. Maps of salinity and leaching were produced by a geostatistical interpolation method based on ordinary kriging. In a great part of the plain, the estimated salinity levels indicate that there is a considerable risk for reduced yield potential for several of the most important crops grown. A rough estimation of the leaching fraction was computed from the soil and water salinity. Computer-based models of salinity ought to facilitate the general planning of irrigation management and crop selection.     

A technique for modeling transport/conversion processes applied to smectite-to-illite conversion in HLW buffers

This paper describes an application of a technique developed for modeling chemical processes in buffer materials that are controlled by a reaction rate and by the transport of one component, which is essential for the process in question to occur. The application described here is the illitization of smectite by fixation of potassium ions in cation exchange positions, and with diffusion of dissolved potassium being the transport process. The technique is verified by comparison with analytical solutions. An overview, based on small models, is given which outlines under what constellations of assumptions the time scale for conversion of the buffer is controlled by reaction rate parameters and under which conditions transport controls this time scale. Examples are given of calculations performed for deposition holes, with potassium being supplied from the surroundings to the upper parts of the highly compacted bentonite buffer. It is concluded that restrictions in nearfield transport capacity have a very significant effect on the conversion time scale. Towards the end of the heating period about 98% of the smectite is found to remain, even for reaction rates and buffer transport conditions that would have left only 10% of the smectite unconverted without nearfield transport restrictions. It is also concluded that the modeling technique can be applied to other, similar, transport/conversion processes.     

Comparison of new individual and hybrid machine learning algorithms for modeling and mapping fire hazard: a supplementary analysis of fire hazard in different counties of Golestan Province in Iran

Natural Hazards - The forest fire hazard mapping using the accurate models in the fire-prone areas has particular importance to predict the future fire occurrence and allocate the resources for...     

Advances in surrogate modeling for storm surge prediction: storm selection and addressing characteristics related to climate change

This paper establishes various advancements for the application of surrogate modeling techniques for storm surge prediction utilizing an existing database of high-fidelity, synthetic storms (tropical cyclones). Kriging, also known as Gaussian process regression, is specifically chosen as the surrogate model in this study. Emphasis is first placed on the storm selection for developing the database of synthetic storms. An adaptive, sequential selection is examined here that iteratively identifies the storm (or multiple storms) that is expected to provide the greatest enhancement of the prediction accuracy when that storm is added into the already available database. Appropriate error statistics are discussed for assessing convergence of this iterative selection, and its performance is compared to the joint probability method with optimal sampling, utilizing the required number of synthetic storms to achieve the same level of accuracy as comparison metric. The impact on risk estimation is also examined. The discussion then moves to adjustments of the surrogate modeling framework to support two implementation issues that might become more relevant due to climate change considerations: future storm intensification and sea level rise (SLR). For storm intensification, the use of the surrogate model for prediction extrapolation is examined. Tuning of the surrogate model characteristics using cross-validation techniques and modification of the tuning to prioritize storms with specific characteristics are proposed, whereas an augmentation of the database with new/additional storms is also considered. With respect to SLR, the recently developed database for the US Army Corps of Engineers’ North Atlantic Comprehensive Coastal Study is exploited to demonstrate how surrogate modeling can support predictions that include SLR considerations.     

Evaluation of scraping treatments to restore initial infiltration capacity of three artificial recharge projects in central Iran

 A limiting factor in developing artificial recharge of groundwater is clogging of the soil surface and consequent reduction of infiltration rates. In order to evaluate the degree of improving infiltration rates by scraping away various amounts of the upper soil materials, a study was conducted at three artificial recharge sites (Kohrouyeh, Bagh-Sorkh, and Kachak) in Isfahan Province, central Iran. Five treatments (T1–T5) were considered. Infiltration was measured: T1, on deposited sediment layer; T2, after removing the sediments; T3, scraping of sediments and 5 cm of soil; T4, scraping of sediments and 10 cm of soil; and T5, removing sediments and 15 cm of soil. Initial soil-moisture content of the sites ranged from 1.0–2.87% for Kohrouyeh, 1.18–3.47% for Bagh-Sorkh, and 1.89–3.93% for Kachak. The main texture of the soils was sandy loam. Clay particles have penetrated to a depth of more than 40 cm in some of the recharge basins. A significant increase in final infiltration rate of T5 as compared to T1 treatment was observed for all recharge sites. The final infiltration rates of T1 and T5 treatments for Kohrouyeh, Bagh-Sorkh, and Kachak sites were 0.35, 7.9; 1.22, 12.3; and 0.93, 6.2 cm/h, respectively. The differences between infiltration rates of T2, T3, and T4 treatments were not statistically significant. It is concluded that on average, the infiltration capacity of the untreated recharge facilities have reached 20.3% of the original values, and that scraping the top sediment layer and 15 cm of topsoil could restore 68.3% of the initial infiltration capacity. Received, July 1998 / Revised, April 1999, May 1999 / Accepted, June 1999