Spatial Modeling of Geometallurgical Properties: Techniques and a Case Study

High-resolution spatial numerical models of metallurgical properties constrained by geological controls and more extensively by measured grade and geomechanical properties constitute an important part of geometallurgy. Geostatistical and other numerical techniques are adapted and developed to construct these high-resolution models accounting for all available data. Important issues that must be addressed include unequal sampling of the metallurgical properties versus grade assays, measurements at different scale, and complex nonlinear averaging of many metallurgical parameters. This paper establishes techniques to address each of these issues with the required implementation details and also demonstrates geometallurgical mineral deposit characterization for a copper–molybdenum deposit in South America. High-resolution models of grades and comminution indices are constructed, checked, and are rigorously validated. The workflow demonstrated in this case study is applicable to many other deposit types.     

Depth classification of underwater targets based on complex acoustic intensity of normal modes

In order to solve the problem of depth classification of the underwater target in a very low frequency acoustic field, the active component of cross spectra of particle pressure and horizontal velocity (ACCSPPHV) is adopted to distinguish the surface vessel and the underwater target. According to the effective depth of a Pekeris waveguide, the placing depth forecasting equations of passive vertical double vector hydrophones are proposed. Numerical examples show that when the sum of depths of two hydrophones is the effective depth, the sign distribution of ACCSPPHV has nothing to do with horizontal distance; in addition, the sum of the first critical surface and the second critical surface is equal to the effective depth. By setting the first critical surface less than the difference between the effective water depth and the actual water depth, that is, the second critical surface is greater than the actual depth, the three positive and negative regions of the whole ocean volume are equivalent to two positive and negative regions and therefore the depth classification of the underwater target is obtained. Besides, when the 20 m water depth is taken as the first critical surface in the simulation of underwater targets (40 Hz, 50 Hz, and 60 Hz respectively), the effectiveness of the algorithm and the correctness of relevant conclusions are verified, and the analysis of the corresponding forecasting performance is conducted.     

Tropical Atlantic climate response to different freshwater input in high latitudes with an ocean-only general circulation model

Tropical Atlantic climate change is relevant to the variation of Atlantic meridional overturning circulation (AMOC) through different physical processes. Previous coupled climate model simulation suggested a dipole-like SST structure cooling over the North Atlantic and warming over the South Tropical Atlantic in response to the slowdown of the AMOC. Using an ocean-only global ocean model here, an attempt was made to separate the total influence of various AMOC change scenarios into an oceanic-induced component and an atmospheric-induced component. In contrast with previous freshwater-hosing experiments with coupled climate models, the ocean-only modeling presented here shows a surface warming in the whole tropical Atlantic region and the oceanic-induced processes may play an important role in the SST change in the equatorial south Atlantic. Our result shows that the warming is partly governed by oceanic process through the mechanism of oceanic gateway change, which operates in the regime where freshwater forcing is strong, exceeding 0.3 Sv. Strong AMOC change is required for the gateway mechanism to work in our model because only when the AMOC is sufficiently weak, the North Brazil Undercurrent can flow equatorward, carrying warm and salty north Atlantic subtropical gyre water into the equatorial zone. This threshold is likely to be model-dependent. An improved understanding of these issues may have help with abrupt climate change prediction later.     

Mechanisms of partly flooded loose sand deposits

Partly flooded loose sand deposits are left back from lignite mining, particularly in east Germany. After flooding, gas inclusions in macropores enhance the pore pressure increase by total stress cycles without drainage, and this can cause the transition into a kind of suspension. Thus, natural and technical changes of surface loads triggered flow slides in chain reactions of sometimes catastrophic extent. Such collapsible deposits were and are stabilized by blast-induced liquefaction and densification. Successive collapse fronts in flooded grain skeletons are seismogeneous, and this appears as wavering in the near-field and leads to power-law spectra in the far-field. Apart from this evidence of mild fractality, the wild fractality of chain reactions eludes mathematical treatment. Therefore, uncommon methods of geotechnical design, operations and monitoring are required.     

Intrinsic groundwater vulnerability determination at the aquifer scale: a methodology coupling travel time estimation and rating methods

Groundwater vulnerability has been subject of much research due to the valuable information it provides concerning groundwater protection and exploitation potential. Up to now, most groundwater vulnerability studies adopt subjective systems of rating the various factors and subsequently, their results are often ambiguous and contradicting. Within the present study a methodology for the estimation of intrinsic groundwater vulnerability at the aquifer scale is presented. The methodology is based on travel time estimation from specified sources of pollution to the aquifer. Besides the deterministic calculation of travel times, the methodology provides a rating system for each pollution source, based on its relative severity and the estimated threat that it poses to the aquifer. Therefore, it can be regarded as a hybrid method that couples the advantages provided by the physically based methods with those of the subjective rating systems. The methodology is applied to the Neon Sidirochorion aquifer, Northeastern Greece, an overexploited aquifer where river waters, sea waters and lake waters interact, causing groundwater quality deterioration to the aquifer. The results indicated that the proposed groundwater vulnerability assessment methodology is well capturing pollution related to saltwater intrusion and agricultural activities, while it is concluded that the conceptual model is significantly affecting the vulnerability assessment results and therefore has to be previously developed.     

Estimation of air-overpressure produced by blasting operation through a neuro-genetic technique

Blasting operations usually produce significant environmental problems which may cause severe damage to the nearby areas. Air-overpressure (AOp) is one of the most important environmental impacts of blasting operations which needs to be predicted and subsequently controlled to minimize the potential risk of damage. In order to solve AOp problem in Hulu Langat granite quarry site, Malaysia, three non-linear methods namely empirical, artificial neural network (ANN) and a hybrid model of genetic algorithm (GA)–ANN were developed in this study. To do this, 76 blasting operations were investigated and relevant blasting parameters were measured in the site. The most influential parameters on AOp namely maximum charge per delay and the distance from the blast-face were considered as model inputs or predictors. Using the five randomly selected datasets and considering the modeling procedure of each method, 15 models were constructed for all predictive techniques. Several performance indices including coefficient of determination (R ²), root mean square error and variance account for were utilized to check the performance capacity of the predictive methods. Considering these performance indices and using simple ranking method, the best models for AOp prediction were selected. It was found that the GA–ANN technique can provide higher performance capacity in predicting AOp compared to other predictive methods. This is due to the fact that the GA–ANN model can optimize the weights and biases of the network connection for training by ANN. In this study, GA–ANN is introduced as superior model for solving AOp problem in Hulu Langat site.     

Modeling the Population Effects of Hypoxia on Atlantic Croaker (<Emphasis Type="Italic">Micropogonias undulatus</Emphasis>) in the Northwestern Gulf of Mexico: Part 1—Model Description and Idealized Hypoxia

We developed a spatially explicit, individual-based model to analyze how hypoxia effects on reproduction, growth, and mortality of Atlantic croaker in the northwestern Gulf of Mexico lead to population-level responses. The model follows the hourly growth, mortality, reproduction, and movement of individuals on a 300 × 800 spatial grid of 1-km² cells for 140 years. Chlorophyll-a concentration, water temperature, and dissolved oxygen (DO) were specified daily for each grid cell and repeated for each year of the simulation. A bioenergetics model was used to represent growth, mortality was assumed stage- and age-dependent, and the movement behavior of juveniles and adults was modeled based on temperature and avoidance of low DO. Hypoxia effects were imposed using exposure effect submodels that converted time-varying exposures to low DO to reduced hourly growth, increased hourly mortality, and reduced annual fecundity. Results showed that 100 years of either mild or intermediate hypoxia produced small reductions in population abundance, while repeated severe hypoxia caused a 19% reduction in long-term population abundance. Relatively few individuals were exposed to low DO each hour, but many individuals experienced some exposure. The response was dominated by a 5% average reduction in annual fecundity of individuals. Under conditions of random sequences of mild, intermediate, and severe hypoxia years occurring in proportion to their historical frequency, the model predicted a 10% decrease in the long-term population abundance of croaker. A companion paper substitutes hourly DO values from a three-dimensional water quality model for the idealized hypoxia and results in a more realistic population reduction of about 25%.     

Formation by silicate–fluoride + phosphate melt immiscibility of REE-rich globular segregations within aplite dikes

Aplite dikes intruding the Proterozoic 1.42(±?3) Ga Longs Peak-St. Vrain Silver Plume-type peraluminous granite near Jamestown, Colorado, contain F, P, and rare earth element (REE)-rich globular segregations, with 40–46% REE, 3.7–4.8 wt% P₂O₅, and 5–8 wt% F. A combination of textural features and geochemical data suggest that the aplite and REE-rich globular segregations co-existed as two co-genetic liquids prior to their crystallization, and we propose that they are formed by silicate–fluoride?+?phosphate (+?S?+?CO₂) melt immiscibility following ascent, cooling, and decompression of what was initially a single homogeneous magma that intruded the granite. The REE distribution coefficients between the silica-rich aplites and REE-rich segregations are in good agreement with experimentally determined distribution coefficients for immiscible silicate–fluoride?+?phosphate melts. Although monazite-(Ce) and uraninite U–Th–Pb microprobe ages for the segregations yield 1.420(±?25) and 1.442(±?8) Ga, respectively, thus suggesting a co-genetic relationship with their host granite, ε_Nd1.42Ga values for the granites and related granitic pegmatites range from ??3.3 to ??4.7 (average ??3.9), and differ from the values for both the aplites and REE-rich segregations, which range from ??1.0 to ??2.2 (average ??1.6). Furthermore, the granites and pegmatites have (La/Yb)_N <50 with significant negative Eu anomalies, which contrast with higher (La/Yb)_N >100 and absence of an Eu anomaly in both the aplites and segregations. These data are consistent with the aplite dikes and the REE-rich segregations they contain being co-genetic, but derived from a source different from that of the granite. The higher ε_Nd1.42Ga values for the aplites and REE-rich segregations suggest that the magma from which they separated had a more mafic and deeper, dryer and hotter source in the lower crust or upper mantle compared to the quartzo-feldspathic upper crustal source proposed for the Longs Peak-St. Vrain granite.     

Carbon and nitrogen isotope,and mineral inclusion studies on the diamonds from the Pozanti–Karsanti chromitite,Turkey

The Pozanti–Karsanti ophiolite (PKO) is one of the largest oceanic remnants in the Tauride belt, Turkey. Micro-diamonds were recovered from the podiform chromitites, and these diamonds were investigated based on morphology, color, cathodoluminescence, nitrogen content, carbon and nitrogen isotopes, internal structure and inclusions. The diamonds recovered from the PKO are mainly mixed-habit diamonds with sectors of different brightness under the cathodoluminescence images. The total δ¹³C range of the PKO diamonds varies between ??18.8 and ??28.4‰, with a principle δ¹³C mode at ??25‰. Nitrogen contents of the diamonds range from 7 to 541 ppm with a mean value of 171 ppm, and the δ¹⁵N values range from ??19.1 to 16.6‰, with a δ¹⁵N mode of ??9‰. Stacking faults and partial dislocations are commonly observed in the Transmission Electron Microscopy foils whereas inclusions are rather rare. Combinations of (Ca_0.81Mn_0.19)SiO₃, NiMnCo-alloy and nano-sized, quenched fluid phases were observed as inclusions in the PKO diamonds. We believe that the ¹³C-depleted carbon signature of the PKO diamonds derived from previously subducted crustal matter. These diamonds may have crystallized from C-saturated fluids in the asthenospheric mantle at depth below 250 km which were subsequently carried rapidly upward by asthenospheric melts.