The relation between Cu/Au ratio and formation depth of porphyry-style Cu–Au ± Mo deposits

Constraints on gold and copper ore grades in porphyry-style Cu–Au ± Mo deposits are re-examined, with particular emphasis on published fluid pressure and formation depth as indicated by fluid inclusion data and geological reconstruction. Defining an arbitrary subdivision at a molar Cu/Au ratio of 4.0 × 10⁴, copper–gold deposits have a shallower average depth of formation (2.1 km) compared with the average depth of copper–molybdenum deposits (3.7 km), based on assumed lithostatic fluid pressure from microthermometry. The correlation of Cu/Au ratio with depth is primarily influenced by the variations of total Au grade. Despite local mineralogical controls within some ore deposits, the overall Cu/Au ratio of the deposits does not show a significant correlation with the predominant type of Cu–Fe sulfide, i.e., chalcopyrite or bornite. Primary magma source probably contributes to metal endowment on the province scale and in some individual deposits, but does not explain the broad correlation of metal ratios with the pressure of ore formation. By comparison with published experimental and fluid analytical data, the observed correlation of the Cu/Au ratio with fluid pressure can be explained by dominant transport of Cu and Au in a buoyant S-rich vapor, coexisting with minor brine in two-phase magmatic hydrothermal systems. At relatively shallow depth (approximately <3 km), the solubility of both metals decreases rapidly with decreasing density of the ascending vapor plume, forcing both Cu and Au to be coprecipitated. In contrast, magmatic vapor cooling at deeper levels (approximately >3 km) and greater confining pressure is likely to precipitate copper ± molybdenum only, while sulfur-complexed gold remains dissolved in the relatively dense vapor. Upon cooling, this vapor may ultimately contract to a low-salinity epithermal liquid, which can contribute to the formation of epithermal gold deposits several kilometers above the Au-poor porphyry Cu–(Mo) deposit. These findings and interpretations imply that petrographic inspection of fluid inclusion density may be used as an exploration indicator. Low-pressure brine + vapor systems are favorable for coprecipitation of both metals, leading to Au-rich porphyry–copper–gold deposits. Epithermal gold deposits may be associated with such shallow systems, but are likely to derive their ore-forming components from a deeper source, which may include a deeply hidden porphyry–copper ± molybdenum deposit. Exposed high-pressure brine + vapor systems, or stockwork veins containing a single type of intermediate-density inclusions, are more likely to be prospective for porphyry–copper ± molybdenum deposits.     

Nonlinear estimation of spatial distribution of rainfall — An indicator cokriging approach

Indicator cokriging (Journel 1983) is examined as a tool for real-time estimation of rainfall from rain gage measurements. The approach proposed in this work obviates real-time estimation of real-time statistics of rainfall by using ensemble or climatological statistics exclusively, and reduces computational requirements attendant to indicator cokriging by employing only a few auxiliary cutoffs in estimation of conditional probabilities. Due to unavailability of suitable rain gage measurements, hourly radar rain fall data were used for both indicator covariance estimation and a comparative evaluation. Preliminary results suggest that the indicator cokriging approach is clearly superior to its ordinary kriging counterpart, whereas the indicator kriging approach is not. The improvement is most significant in estimation of light rainfall, but drops off significantly for heavy rainfall. The lack of predictability in spatial estimation of heavy rainfall is borne out in the integral scale of indicator correlation: peaking to its maximum for cutoffs near the median, indicator correlation scale becomes increasingly smaller for larger cutoffs of rainfall depth. A derived-distribution analysis, based on the assumption that radar rainfall is a linear sum of ground-truth and a random error, suggests that, at low cutoffs, indicator correlation scale of ground-truth can significantly differ from that of radar rainfall, and points toward inclusion of rainfall intermittency, for example, within the framework proposed in this work.     

Spatial rainfall estimation by linear and non-linear co-kriging of radar-rainfall and raingage data

The feasibility of linear and nonlinear geostatistical estimation techniques for optimal merging of rainfall data from raingage and radar observations is investigated in this study by use of controlled numerical experiments. Synthetic radar and raingage data are generated with their hypothetical error structures that explicitly account for sampling characteristics of the two sensors. Numerically simulated rainfall fields considered to be ground-truth fields on 4×4 km grids are used in the generation of radar and raingage observations. Ground-truth rainfall fields consist of generated rainfall fields with various climatic characteristics that preserve the space-time covariance function of rainfall events in extratropical cyclonic storms. Optimal mean areal precipitation estimates are obtained based on the minimum variance, unbiased property of kriging techniques under the second order homogeneity assumption of rainfall fields. The evaluation of estimated rainfall fields is done based on the refinement of spatial predictability over what would be provided from each sensor individually. Attention is mainly given to removal of measurement error and bias that are synthetically introduced to radar measurements. The influence of raingage network density on estimated rainfall fields is also examined.     

Microscopic observation and contact stress analysis of granite under compression

To study the damage process of microscale to macroscale in coarse-grained granite specimen under uniaxial compressive stress, we have observed micro-damage localization and propagation by using a newly developed experimental system that allows us to observe the damaging process continuously.

The results showed that pre-existing microcracks lead to macroscopic shear fracture through the damage development process. The mechanism of micro-damage initiation in a granite specimen under uniaxial compressive stress may be considered for two cases. One is that two grains such as quartz and feldspar contact each other in the same direction as the axial stress, and the other is that a biotite grain inclined to the axial stress direction is surrounded by feldspar grains. The homogenization theory was applied to verify numerically the micromechanics of stress-induced damage in the mineral contacts. Local stress distribution in the periodic-micro structure was also calculated by the homogenization theory. It is shown that this analysis, which takes into account the initial state of the specimen, is well adapted to the behavior of two grains for which microcracking is the fundamental mechanism of damage.     

Hydrological modelling uncertainty analysis for different flow quantiles: a case study in two hydro-geographically different watersheds

This study focuses on analysis of hydrological model parameter uncertainty at varying sub-basin spatial scales. It was found that the variation in sub-basin spatial scale had little influence on the entire flow simulations. However, the different sub-basin spatial scales had a significant impact on the reproduction of the flow quantiles. The coarser sub-basin spatial scale provided a better coverage of most prediction uncertainty in observations. However, the finer sub-basin spatial scale produced the best single simulation output closer to the observations. In general, the optimal sub-basin spatial scales (ratio to the entire watershed size) in the two test watersheds were found to be in the ranges 14–19% and 2–4% for good simulation of high and low flows, respectively. It is therefore worthwhile to put more effort into reproducing different flow quantiles by investigating an appropriate sub-basin spatial scale.     

Wave prediction in a port using a fully nonlinear Boussinesq wave model

A wave forecasting system using FUNWAVE-TVD which is based on the fully nonlinear Boussinesq equations by Chen(2006) was developed to provide an accurate wave prediction in the Port of Busan, South Korea. This system is linked to the Korea Operational Oceanographic System(KOOS) developed by Park et al.(2015). The computational domain covers a region of 9.6 km×7.0 km with a grid size of 2 m in both directions, which is sufficient to resolve short waves and dominant sea states. The total number of grid points exceeds 16 millions,making the model computational expensive. To provide real-time forecasting, an interpolation method, which is based on pre-calculated results of FUNWAVE-TVD and SWAN forecasting results at the FUNWAVE-TVD offshore boundary, was used. A total of 45 cases were pre-calculated, which took 71 days on 924 computational cores of a Linux cluster system. Wind wave generation and propagation from the deep water were computed using the SWAN in KOOS. SWAN results provided a boundary condition for the FUNWAVE-TVD forecasting system. To verify the model, wave observations were conducted at three locations inside the port in a time period of more than 7 months. A model/model comparison between FUNWAVE-TVD and SWAN was also carried out. It is found that, FUNWAVE-TVD improves the forecasting results significantly compared to SWAN which underestimates wave heights in sheltered areas due to incorrect physical mechanism of wave diffraction, as well as large wave heights caused by wave reflections inside the port.     

A comparative study of toxicity identification using Daphnia magna and Tigriopus japonicus: implications of establishing effluent discharge limits in Korea

In Korea, the new permission criteria for industrial effluents based on Daphnia magna acute toxicity tests will be gradually implemented starting from 2011. Thus, in this study, toxicity assessment and identification using a marine species (Tigriopus japonicus) and the freshwater species (D. magna) was comparatively investigated. Effluent from an acid mine drainage treatment plant showed acute toxicity toward both organisms due to low pH, which was removed by neutralization of the effluent. Additionally, evaluation of the effluent of an electronics company revealed that Cu was attributable to the observed toxicity, and the effluent was more toxic toward T. japonicus than D. magna. Moreover, effluents from a metal plating factory were acutely toxic toward D. magna (6.50 TU), while they were not toxic against T. japonicus. Toxicity identification revealed that the high level of Cl- (12,841 mg L(-1)) was the cause of toxicity. Thus, the effluents had no effect on the marine species, T. japonicus. These findings suggest that a marine species rather than a freshwater species is more desirable for toxicity assessment of industrial effluent discharged into the saltwater, and thus should be considered in the legislation of toxicity-based discharge limits in Korea.