Discrimination of Mineralized Rock Types in a Copper-Rich Volcanogenic Massive Sulfide Deposit Through Fast Independent Component and Factor Analysis

Natural Resources Research - Multivariate methods are useful for simplifying the interpretation of variables in geochemical data and are widely used to uncover relationships between elements that...     

Compounding effects of human activities and climatic changes on surface water availability in Iran

Climatic Change - By combining long-term ground-based data on water withdrawal with climate model projections, this study quantifies the compounding effects of human activities and climate change...     

Optimizing urban stormwater control strategies and assessing aquifer recharge through drywells in an urban watershed

Hydrogeology Journal - A coupled simulation-optimization model (SOM) is developed in this work that links the US Environmental Protection Agency’s Storm Water Management Model (SWMM) with a...     

Spatial and temporal variations in the incidence of dust events over Iran

Large near-field tsunamis pose a significant threat to the Canadian West Coast due to its proximity to the circum-Pacific belt where a significant tsunami-inducing earthquake event from the Cascadia subduction zone is expected. This study investigated the risks associated with such an event in terms of pedestrian evacuation needs and plans for the Town of Tofino, a small community located on the West Coast of Vancouver Island. The population-at-risk within the hazard zone and its ability to evacuate to safety is evaluated using anisotropic path-distance modelling. Mitigation measures, such as vertical evacuation buildings, are quantitatively evaluated. Site-specific inundation modelling was not performed as part of this study; tsunami hazard and safe zones were computed using a range of run-ups varying between 3 and 25 m. It was established that up to 80% of the population is within the maximum hazard zone considered. This evacuation modelling exercise indicates that a maximum of 13% of the population would have insufficient time to reach safety when using a mobility-impaired ambulatory speed. The use of three vertical evacuation buildings can reduce the risk of losing population in this category by 99%. Although some conservative assumptions were used (vertical datum at higher high water, reductions in safe zones by generalization process and mobility-impaired evacuation speeds), the evacuation potential is likely overestimated due to the coarseness of the topographic data used in the evacuation modelling and from an overestimated first wave arrival time. This is the first Canadian study which used anisotropic evacuation modelling to evaluate the vulnerability of a Canadian community to tsunami inundation.

     

Multiple-point geostatistical simulation of dykes: application at Sungun porphyry copper system,Iran

Variogram-based methods are not capable of capturing high (>2) order statistics since the variogram measures the relationship between two points at a time only. Multiple-point geostatistics (MPS) has brought new insights into many geological modeling problems. The application of MPS methods has been well documented in realizing complex geological patterns. These methods have often been used in reservoir characterization since their advent in recent decades. The frequent non-linear behaviors of geologic continuity are not limited to reservoirs, but mineral deposits bear complicated formations in many cases. Relying on the power of MPS methods and considering the complexity of geological scenarios in mineral deposits, we have applied MPS in the modeling of mineral deposits. A training image (TI) is produced using geological data from upper horizons of a porphyry copper ore deposit which have been mined out during the previous mining operations. In this study, the SNESIM algorithm has been used. A number of realizations are produced using this multiple-point geostatistical method. Extensive validation steps are performed considering the TI as the reference model. These validations first show that the TI is representative for the domain under study and also illustrates some degrees of similarity between the TI and the realizations. Despite simplifications made to the problem, the application of MPS in mineral deposit modeling still faces many challenges.     

Alkaline basalt from the Central Iran, a mark of previously subducted Paleo-Tethys oceanic crust

In western part of the CEIM (Central-East Iranian Microcontinent) (Bayazeh area, Isfahan province, Iran), a series of Paleozoic basaltic rocks, occur. Major minerals of these basalts are olivine, clinopyroxene (diopside, augite), plagioclase (albite), sanidine, amphibole (kaersutite), phlogopite, ilmenite and magnetite. Secondary minerals include epidote, pumpellyite, albite, calcite and chlorite. Olivine and clinopyroxene are as phenocryst, while feldspars are restricted to groundmass. Chemical composition of clinopyroxenes indicates crystallization during ascending of magma. Geochemical analysis of whole rock samples shows that these rocks are characterized by low SiO₂ (43.21–48.45 wt %), high TiO₂ (1.81–3.00 wt %) and P₂O₅ (0.18–0.34 wt %). Petrography, chemistry of clinopyroxenes and whole rock analyses reveal an alkaline nature of these basalts. They are enriched in alkalis (Na₂O + K₂O = 4.1–7.7 wt %), LILE, HFSE and LREE. The Bayazeh alkali-basalts present strong enrichment in LREE relative to HREE (La/Lu ratio = 77.6–119.6) and were dominantly derived from partial melting of a metasomatized asthenospheric garnet-amphibole lherzolite. Field relationships reveal that junction of faults in west of the Bayazeh prepared a suitable path for ascending of magma from deep regions to surface and intra-plate continental magmatism. The Paleo-Tethys subduction from lower to upper Paleozoic is too enough for mantle enrichment in volatiles and basaltic alkaline magmatisrn in upper Paleozoic of Bayazeh area.     

Nodal-based three-dimensional discontinuous deformation analysis (3-D DDA)

This paper presents a new numerical model that can add a finite element mesh into each block of the three-dimensional discontinuous deformation analysis (3-D DDA), originally developed by Gen-hua Shi. The main objectives of this research are to enhance DDA block’s deformability. Formulations of stiffness and force matrices in 3-D DDA with conventional Trilinear (8-node) and Serendipity (20-node) hexahedral isoparametric finite elements meshed block system due to elastic stress, initial stress, point load, body force, displacement constraints, inertia force, normal and shear contact forces are derived in detail for program coding. The program code for the Trilinear and Serendipity hexahedron elements have been developed, and it has been applied to some examples to show the advantages achieved when finite element is associated with 3-D DDA to handle problems under large displacements and deformations. Results calculated for the same models by use of the original 3-D DDA are far from the theoretical solutions while the results of new numerical model are quite good in agreement with theoretical solutions; however, for the Trilinear elements, more number of elements are needed.     

Spatial modelling of hazardous elements at waste dumps using geostatistical approach: a case study Sarcheshmeh copper mine,Iran

Mining activities and resulting wastes can be considered as one of the most important sources of hazardous elements in the environment. Knowledge of the spatial distribution of toxic elements in waste dump systems is necessary to assess environmental hazard and strategy. To achieve this goal, this paper investigates spatial distribution of toxic elements using statistical and geostatistical analysis. A total of 58 soil samples were collected, and the amount of As, Cd, Co, Cr, Cu, Mn, Mo, Ni, Pb and Zn was then determined at “Sarcheshmeh” copper mine waste dumps. In order to evaluate the presence of multivariate outliers, Mahalanobis distance technique (D ²) was applied and the multivariate outlier samples were removed. This resulted in an increase in correlation coefficient. To reduce dimension of data set, principal component analysis was applied and four principal components were determined which indicate 83.463% of the total variance of data set. Estimated PCs together with the toxic elements maps based on the ordinary kriging display aggregation of toxic elements in some parts, and validity of predictions was evaluated using the leave-one-out cross-validation method. The regression coefficients of estimated and observed values presented the reliability of the kriging estimates. Sequential Gaussian simulation method was applied for principal components due to similar results of estimated principal components and toxic elements. The results of simulation maps are almost identical to estimated outcomes.