Search Documents

By C F. Prins, C LantuTjoul

Mineralisation in alluvial diamond or vein gold deposits is often hosted in trapping mechanisms. For instance, diamonds are nested in bedrock depressions such as gullies and potholes which were carved out by the sea long before the diamond deposition took place. Similarly, gold mineralisation can precipitate in the necks of a boudinage structure. In such deposits, the spatial distribution of the mineralisation is governed by the mutual arrangement of the trapsites and by the grade variations within each trapsite. Reliable estimates can be obtained only if they rest on a proper knowledge of both degrees of variability. This creates a challenging sampling problem, as during the exploration phase it is important to have a good global mean estimate to enable decision making towards the future. A similar problem occurs for local estimates at a later stage of development. To understand and manage the sampling for global means, specific sampling procedures are required based on an appropriate choice of sample support size and spacing. Estimates from such sampling campaigns should be unbiased, and the less variability they display, the better the sampling parameters. The performance and impact of sampling parameters are tested and compared on different simulated deposits.

Jan 1, 2005

By G Chitombo, A Perryman, A D. Campbell, J White, B Hollis

This paper describes the development and implementation of an optimisation tool developed as part of the Mass Mining Technology (MMT) Project. The tool, known as the SLC Schedule Optimiser, is capable of optimising the schedule of development and production activity for a sublevel cave, given a fixed draw strategy, for the life of the mine.This paper starts by describing the development and production scheduling problem for a sublevel cave, the geotechnical constraint that a schedule of activities must adhere to, and outlines the pre-existing manual planning procedures used to produce schedules.Following this, the paper then outlines the optimisation tool developed, describes the changes to the planning process it has enabled and finally quantifies the improvement to operations this tool has brought using Glencore’s Ernest Henry Mining Operations as a case study.CITATION:White, J, Perryman, A, Campbell, A D, Hollis, B and Chitombo, G, 2016. Life-of-mine development and production scheduling optimisation for sublevel cave operations, in Proceedings Seventh International Conference and Exhibition on Mass Mining (MassMin 2016), pp 133–138 (The Australasian Institute of Mining and Metallurgy: Melbourne).

May 9, 2016

By D E. G Connelly, K Y. Wong

Gold at the Harbour Lights Mine (HLM), WA occurs principally in pyrite (Py) and arsenopyrite (AsPy) minerals rendering it inaccessible to direct cyanidation. Froth flotation produces a concentrate with low mass recovery and with high recoveries of sulfides. After investigating several alternatives, a novel flotation process has been developed to separate the concentrate into pyrite and arsenopyrite rich fractions. The pyrite fraction is finely ground and cyanide leached, while the arsenopyrite is treated off-site. This paper describes the separation process from laboratory bench scale through to commercial operation at the HLM gold plant.

Jan 1, 2005

Conditional simulation can be used to determine the probabilities of block iron-grades being above selected cut-offs. However, there are much faster alternatives that can produce acceptable results. In this paper blasthole data from a selected area of Joffre Member hosted ore of the Brockman Iron Formation at the Mt Whaleback orebody is used as the basis of a case study to assess alternatives to simulation. Iron grade, in per cent, is interpolated from a blasthole dataset into a block model using ordinary kriging. Samples are then removed from this blasthole dataset to produce data subsets. A complete distribution of iron grades for points or blocks is obtained utilising the data subsets and the following techniques: conditional simulation, median indicator kriging, the probability from conditional expectation and confidence intervals attached to the estimation of blocks in the framework of the discrete Gaussian model. The complete blasthole dataset is considered to represent the true iron grades at sample point support while the ordinary kriging estimate derived from these samples is considered to represent the true block grades. The distributions obtained from point estimates, block estimates and simulations are compared to each other and to the true estimates of iron grades using graphs and statistics. The conditional expectation technique reproduced the results obtained from conditional simulation extremely well and it is a viable alternative to simulation at point support. At block support the confidence intervals technique reproduced the results of the conditional simulation more closely than the median indicator kriging and in the data subsets tested the confidence intervals technique actually reproduced the true grade distribution more closely than the conditional simulation.

Jan 1, 2005

By Y Ren

The algorithms employed by the commercial tools are either heuristic algorithm or iterative optimisation of a subsystem. This often results in a suboptimal schedule. In this paper, the author presents a new scheduling algorithm that incorporates with a novel stockpile configuration system, an improved algorithm on multiple processes and a mixed-integer linear programming (MILP) global optimisation model. It can simultaneously optimise stockpiles, blending, multiple processes as well as multiple mines. In addition, it allows the user to interact with the solution to adjust the solution to be practical and desirable.CITATION:Ren, Y, 2013. A simultaneous optimisation tool for open pit mine scheduling, in Proceedings World Gold 2013 , pp 427-432 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Sep 26, 2013

By M Ogle

"Booster fans are relatively rare in Australian coal mines. However, as mines become progressivelydeeper it is likely that such fans will be increasingly adopted as a method to sustain or improve the ventilation of underground workings. The reasons why booster fans may be an att ractive alternative include: increased cost of sinking additional ventilation shafts with greater seam depths; environmental objections to shaft sinking projects; and increasing urbanisation and surface constraints limiting availability of potential ventilation shaft sites. As Tahmoor Colliery movedinto its northern lease the ventilation constraints became the limiting factor on production. Thelimitations of the old workings meant that the existing ventilation system would be inadequateto sustain the required production. Extensive ventilation modelling was carried out to evaluatevarious potential solutions and the recommended option was the installation of a booster fanor fans, to gain underground access to an already approved shaft site in the northern lease.Proposals were sought from recognised fan manufacturers and various alternative arrangements were tendered. Ultimately the choice was between one large centrifugal fan and a set (three or four) of smaller æauxiliary-likeÆ fans in parallel. The successful tender was for a single centrifugal, an arrangement very similar to the booster fans that have been operating eff ectively at Westcliff Colliery for a number of years. This paper describes the current ventilation network, evaluation of the ventilation options, the fan selection process, and the project status to date."

Jun 1, 2012

By R J. Hardy

The aim of this study is to examine the volumetric and mass movement characteristics of load and haul operations using statistical modelling of relatively large samples of bucket loads that are transported as small samples in the form of truck payloads. Generally, the factors influencing bucket load variability are random, difficult to predict and to control. In the process of loading trucks the underlying variability of individual bucket loads will determine the variability of truck payloads. Bucket loads of material prepared by drilling and blasting for loading into trucks are fundamental quantities in open pit mine production. A study of payload data from a large surface mining operation indicates that bucket loads and truck payloads that are a small sample of bucket loads are normally distributed. Truck payloads have a reduced dispersion compared with bucket loads that is an inverse 0.5 power function of the number of passes per truck payload. Comparison of statistical modelling results with the æ10/10/20Æ policy promoted by Caterpillar Global Mining Inc, indicate that compliance with the æ10/10Æ part of the policy, applying to truck overloading is facilitated, as the coefficient of variation of truck payloads decreases, as a result of the increasing number of loading passes. Compliance with the more stringent æ20Æ part of the policy, relating to safe steering and braking is also facilitated, albeit at lower coefficient of variation for similar reasons.

Sep 26, 2011

By Abrahamsson S, Niklasson B

In a follow-up of 12 production blasts in a quarry at the Mrica project on Java, tightly spaced porous bands were found to govern fragmentation to the extent that the average fragment size k5o was nearly independent of the specific charge within the range 0.30-0.56 kghn3. This fording resulted in a reduction of the running costs for drilling and blasting by 20-40%. Our work has given considerable practical benefits and insights into the blasting process, even if the results were of no direct use to check out fragmentation predictions based on models, which was one of our main Goals.

Jan 1, 1990

By D Kentwell

IGO Ltd’s (IGO) Nova Operation (Nova) is 160 km east-north-east of the town of Norseman in Western Australia (WA). Discovered by Sirius Resources NL (Sirius) in 2012 and acquired by IGO in 2015, the Nova-Bollinger Deposit (Nova-Bollinger) occurs within the Fraser Zone of the AlbanyFraser Orogen. Nova-Bollinger is a complex nickel-copper-cobalt (Ni-Cu-Co) magmatic intrusive sulfide deposit, which is characterised by high-grade folded and undulating breccia lodes that host the bulk of the payable metals. Significant sulfide mineralisation is also found in several magmatic intrusive styles including net-textured, disseminated and blebby styles. A smaller but material volume of sulfides has also been remobilised into the deposit’s footwall rocks as ‘stringer’ bands in metasediments, or as discrete massive sulfide lenses. When interpolating Mineral Resource estimate (MRE) block grades in the presence of locally varying grade connectivity trends, such as folds or undulations, the sample search neighbourhood, which is used to identify the composites subset that is used for each block estimate, should be oriented to match the local trends. In theory, the most robust method of estimation where locally varying angles (LVAs) are applied to search neighbourhoods, are techniques where the composites and the blocks to be estimated are transformed to an unfolded coordinate system (UCS). However, UCS methods can be tedious to implement and require subjective decisions to be made that cannot always be fully defended. An alternative LVA method commonly used in the industry is dynamic anisotropy (DA), which is much simpler to implement and is often assumed (but not that often demonstrated) to result in similar results to UCS-style estimation. In 2018, Nova’s Mine Geology team (MGT) estimated the grades and density of the Nova-Bollinger MRE using the UCS approach. For the 2020 MRE update, Nova’s MGT tested and compared the USC grade interpolation method with the DA approach. As applied to the Nova-Bollinger’s 2020 MRE, the comparison, advantages and disadvantages of the two estimation methods are the subject of this paper

Mar 22, 2022

By P B. Hills

This paper discusses the computerisation of the drill hole database, some historical records and the creation of a three-dimensional model of the gold bearing structures for resource estimation and mine design at the Beaconsfield Gold Mine in northern Tasmania. The phase of the work reported here commenced in November 1997 and was completed in March 1998. The historical mining on the Tasmania Reef had reached a depth of 450 m below surface when work ceased in 1914. Production up to this time was 854 000 ounces gold. Recent activity has concentrated on testing depth extensions of the previously mined ore structure, to a depth of 1000 m below surface, and the establishment of an underground mining operation to exploit this resource. The lode structure at Beaconsfield is a typical mesothermal vein style deposit. Gold is contained within narrow quartz, carbonate and minor sulphide reefs or veins. Reef width is variable from a few centimetres to several metres with æsplitsÆ or bifurcations being common. The reef averages 2.7 m true width in the economic zones. The key feature of this style of mineralisation is that reef continuity can be confirmed by drilling to a high level of confidence, but the gold grades are somewhat more erratic in nature. The resource estimates were based entirely on drill data. The drill intercepts on the reef were coded for each split directly into a database and these allocations validated graphically. The distribution of grades and widths show an increasing gold grade with reef width for both the drill hole and historical data. The same histogram patterns, mixed populations and overall arithmetic average grade for the historical mine (at 21.52 g/t > 3 g/t Au) and new drill data (at 21.55 g/t > 3 g/t Au) would indicate that there has been no major change in the overall style or grade of the ore from the historically mined areas to the new resource. The reef, host geology, area of historical workings and resource zones at the Beaconsfield Mine were defined as a series of closed wireframes. The Tasmania Reef was modelled as five mutually exclusive wireframes, one each for the Main Reef, Footwall split, Second footwall split, Third footwall split and Hangingwall split. The data for this exercise was extracted directly from the drill hole database, stored in Microsoft Access format. The formation of the solids was largely automated with manual validation and refinements allowing for multiple runs using different interpolation methods and combinations of drill holes. A rotated 3D-block model in the plane of the Tasmania Reef was created to contain the grade model and allow for tonnage, grade and reef width estimates. An inverse distance squared (ID2) interpolation method was used with ranges of 25 m along strike and 125 m down plunge with blank areas infilled with 50 m by 250 m. Each reef split was filled separately using only drill intercepts from that split. The method was independently validated used MIK methods. Longitudinal sections have been prepared showing gold grade, reef width and width +grade. Estimates were also made for arsenic, sulphur and other minor elements. The long sections show three distinct mining zones. A method is presented for the conversion of the resources to reserves by the addition of mining parameters for minimum mining width, minimum pillar, cut-off grade and internal and external dilution. The selection is done on columns within the block model, and allows for all splits to be considered at once. Multiple runs can be made for varying combinations of parameters to test different mining methods and their impact on reserves, mine life and NPV.

Jan 1, 1999

Opening AddressThank you for inviting me here today to open what promises to be a significant conference on the technological aspects of extractive metallurgy. I would especially like. to welcome delegates who have travelled great distances in order to contribute to, and be participants in, this examination of the technical innovations which have always been a crucial part of our industry.

Jan 1, 1992

This study examines the application of Landsat TM data to the identification of hydro- thermal alteration in the arid terrain of the El Salvador region of north Chile. Numerical reflectance values were extracted from the digital Landsat TM data for a variety of rock surfaces, including four parts of the El Salvador gossan, for each of six spectral bands. These reflectance values were analysed statistically in order to select the three spectral bands, combined as a colour composite image, that are most efficient in discriminating different varieties of alteration and for general geolog- ical interpretation. The most cost effective com- posite image for this area is a combination of bands 1,4 and 7 as the blue, green and red components respectively, with simple contrast enhancement. This image is superior to some much more expensive enhancement techniques and allows unambiguous identification of areas of hydrothermal alteration larger than about 50m. The display includes a practical guide to the use of Landsat TM imagery for volcanic gold Exploration.

Jan 1, 1987

By E M. Gama, L Resende

The tailings are waste from mineral processing operations. The type of tailings of iron ore processing are classified into sandy tailings and mud. Already the overburden are showing significant levels of useful mineral. Both iron ore`s waste have in your composition phyllosilicates/clay that indicates the potential for use in niche markets as paving and cement. According to several studies, there are various technologies in bench and pilot levels that could be used on a large scale but there is not an economic environment, tax, regulatory and tax set. Many of these technologies are associated with the reuse of tailings and overburden in construction. All these factors influence for that so many initiatives remain in search formats and cannot advance to the new business format. There is no way to set goals of reuse of this waste, as already occurs in China, if there is not a fiscal, tax and economic structure to support such measures. The way to foster and ensure that reuse of tailings and overburden is an incentive for promoting a culture of circular economy that would allow the emergence of an ecosystem between the private, public and Government sectors with the objective of accelerating and subsidize the entire network of reuse of this waste. The closure of the mine could be the tool used to ensure and boost this ecosystem of reutilization to align it with the vocation of the area of the mining law. Following this model, the end of the mining activity could provide a financial return to the local population, a more effective environmental compensation and a reduction in costs related to the closure of the mine. In this study are shown how the Vale is working in the reuse of tailings and overburden.CITATION:Resende, L and Gama, E M, 2018. The challenges of the reuse of tailings and overburden, in Proceedings Mine Waste and Tailings Stewardship Conference 2018, pp 117–125 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 23, 2018

By C Li

Removal of gangue minerals from ilmenite concentrates can be accomplished by relatively low temperature roasting in controlled gas atmospheres to selectively enhance the magnetic separation characteristics of the ilmenite. At roasting temperatures below 800¦C the magnetising reactions are often not consistent with phase equilibria considerations. Metastable phases can nucleate and control the roasting products. In weathered ilmenites the phase assemblages and microstructures in the altered ilmenite grains can influence the roasting reactions. In this paper we give an overview of the phase chemistry of chemically weathered ilmenites and of metastable phases formed by low temperature roasting. We also report the results of roast experiments conducted in a neutral gas atmosphere. Major increases in magnetic properties after short roasting times at low temperatures were obtained for a range of ilmenite concentrates. For example heating a secondary ilmenite concentrate from Western Australia for 15 minutes in nitrogen at 600¦C gave an increase in the three kilogauss magnetic fraction from three to 92 wt per cent. The solid state phase changes accompanying the inert atmosphere roasts were characterised by combining elemental analyses with quantitative phase analyses of the ilmenite concentrates and their roast products. It was found that the phase reaction that is principally responsible for the increase in the magnetic properties of weathered ilmenites involves a transformation of disordered ferrous-containing pseudorutile to an ordered ilmenite-hematite solid solution containing typically ten to 35 mole per cent of the hematite component. Examples are given of the application of the inert atmosphere roast to effect separations of ilmenite from gangue minerals.

Jan 1, 2001

By N Belluz, L Johnson, N Wiwoho, P Manning

The Kucing Liar skarn deposit is located in the Ertsberg mining district, which also contains other skarn deposits and the Grasberg porphyry copper û gold deposit. The Ertsberg district is located in the province of Irian Jaya, IndonesiaÆs easternmost province. Mining development began in the Ertsberg district in 1972 with the mining of the Ertsberg deposit. Daily production in the district now averages about 230 000 tonnes of ore, with the bulk of the ore coming from the Grasberg open pit and the balance from underground mining of skarn deposits. During 1998, 45 218 additional metres of drilling were conducted to delineate and understand the Kucing Liar orebody. The deposit was first intersected in drill holes in 1994 and then was intersected in underground workings located about 1000 metres from the discovery drill hole in 1995. Initially, the two intersections were thought to be different orebodies but are now realised to be different parts of the same orebody. The Kucing Liar deposit consists of magnetite û copper û gold replacement and skarn mineralisation in favorable Tertiary and Cretaceous units. The deposit is approximately 2000 metres in length and averages about 200 to 250 metres in thickness. The width varies from 400 to 600 metres. Drilling in 1998 was directed towards delineating the deposit on its NW and SE ends and to determine its down dip extension. Reserves for the deposit at the end of 1998 stand at 320 million tonnes at 1.4 per cent copper, 1.4 g/t gold, and 5.3 g/t silver. Many important questions remain to be answered concerning the relationship of the Kucing Liar orebody to the major faults in the district and to the Grasberg intrusive complex which contains the Grasberg orebody. Continued exciting developments are anticipated as work to understand the deposit by drilling and geologic study continues.

Jan 1, 1999

By M Dusci, A Trueman, S Poos, R Baumgartner, M Brittan

The Canahuire epithermal Au-Cu-Ag deposit in southern Peru is a complex and challenging deposit for metallurgical studies. It has issues of gold locking in sulfides and non-sulfides, deleterious elements and acid rock drainage (ARD). From scoping study to feasibility study, there has been considerable improvement in the understanding of the important geometallurgical characteristics that need to be measured and modelled. This paper presents the current understanding and the requirements to continue to provide optimal models for mining studies.Significant effort has been put into understanding the genesis of the deposit. Improved spatial and temporal understanding of lithology, alteration, structure and mineralogy was used to build a geological and alteration model of the deposit.Early metallurgical programs quantified the effect of head grade on process plant recoveries, but the geometallurgical and statistical analyses did not identify other physical characteristics that influence these recoveries. The studies improved the understanding of the nature of locked gold which can occur as minute inclusions (< 5 ¦m) in sulfides (primarily pyrite with lesser arsenopyrite, pyrite being 100 times more abundant than arsenopyrite) and gangue minerals (quartz and carbonates) as well as solid solution in arsenopyrite and pyrite.The presence of pyrite and the limited buffering capacity of waste rocks increases the potential ARD, which will require mitigation costs before, during and after operations. Waste rock ARD potential has been classified and incorporated in the geometallurgical model. The current classification of waste rock at Canahuire, using international classification standards, is limited due to the complexity of the deposit.Complex deposit metallurgy necessitates a complex process flow sheet, including flotation and cyanidation. Economic recovery will be constrained by the relationships between economic variable grade, other elements and mineralogy. These constraints and inter-relationships are difficult to establish throughout the deposit with a single test. A follow-up geometallurgical program is under investigation and will attempt, with the help of additional test work, to map areas of higher and lower recoveries. It is hoped that metallurgical response can be linked to measurable parameters other than head grade. If this is possible, these factors can be incorporated into future geometallurgical models.CITATION:Baumgartner, R, Dusci, M, Trueman, A, Brittan, M and Poos, S, 2013. Building a geometallurgical model for the Canahuire Epithermal Au-Cu-Ag Deposit, Southern Peru - past, present and future, in Proceedings The Second AusIMM International Geometallurgy Conference (GeoMet) 2013 , pp 51-58 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Sep 29, 2013

By K G. Fleetwood

This paper proposes an approach to determine the influence of non-ideal detonation characteristics related to the charge diameter, explosive type and confinement conditions on the available shock energy of ANFO and a representative chemically-sensitised (gassed) emulsion. The proposed degrees of non-ideality of the two different blasting agents are indicated by the Gurney ratios, determined from existing cylinder expansion testing data. The influence of the non-ideal detonation properties on blast design and performance are also compared using a series of theoretical, equal powder factor open pit blasting patterns with hole diameters ranging from 102 mm to 270 mm. The non-ideal shock energy factor (NSEF), which is the available shock energy per tonne of blasted material (MJ/t), is used to compare the different blasting patterns as charged with either ANFO or the gassed emulsion of equal charge density (0.85 g/cm3). The intended purpose of this investigation is to provide explosive users with an approach to compare different explosive types or charge diameters in the absence of dedicated non-ideal detonation modelling, which is generally not available to explosives users.The results of the comparisons suggest that at low charge densities, the example gassed emulsion displays a higher velocity of detonation (VOD) than ANFO in smaller charge diameters but similar degrees of non-ideality (based on proposed Gurney ratios). The results also suggest that the confined, non-ideal gassed emulsion VOD at hole diameters greater than 270 mm can be less than ANFO due to the low density effects on the emulsion matrix. The differences in non-ideal detonation properties and Gurney ratios between the ANFO and gassed emulsion resulted in gassed emulsion NSEF values ranging from 15 per cent lower to 65 per cent higher than ANFO for equal powder-factor patterns in blasthole diameters from 102 mm to 270 mm.

Nov 8, 2011

By S Nie, I Bergqvist, F Ouchterlony

In blasting, a few or many cracks are driven from the borehole into the rock. But what causes the cracks? The most common theory of breakage consists of two stages; first the shock wave causes radial cracks to form around the hole then the gases penetrate into the cracks, and widen them and make them longer. Another theory presented by Brinkmann suggests that the back damage is primarily controlled by shock and that the gas penetration is the mechanism controlling breakout of the burden. He did his experimental work using blasthole liners. Recent research at SveBeFo has examined this matter further. In a quarry a number of benching holes have been blasted simultaneously. In some of these holes tubular Swellex bolts were inflated and decoupled charges put inside the tubes without stemming. Other holes were identically charged but without the lining. Finally some holes were also stemmed. After blasting the cracks in the remaining rock were studied. There was no difference in crack lengths between holes charged normally (no stemming) and holes where the charges were inside the bolts. On the other hand when stemming was used, the crack lengths increased for some explosives but remained the same for an emulsion explosive. In another set up blasted granite blocks were charged in the same way as above. Then we could also measure the bore hole pressure. The pressure gauge consists of a small carbon resistor inside a steel cylinder. It is called LHM (Location-fixed Hydrodynamic Measuring Cup) and is placed at the bottom of the hole. A smaller exit hole from the bottom is drilled for the cables. The paper presents the technique and the results obtained from both the quarry blasting and the blasting of the blocks.

Jan 1, 2001

By Z U Bajwah

During the Palaeoproterozoic, a large number of granitic intrusions were emplaced in the various parts of the Pine Creek Geosyncline (PCG). In the central region, the granitic magmatism is represented either by the coalesced mass, eg the Cullen Batholith, or individual plutons such as the Mount Bundey Granite which together with the Mount Goyder Syenite intrude the metasedimentary strata, causing wide spread contact aureoles. Some of these granitic intrusions are temporally and genetically related to gold (and less commonly base-metals) mineralisation present in the adjacent metasediments. These plutons are characterised by the presence of biotite and homblende which appear late in the paragenetic sequence. Electron microprobe analyses of biotite and homblende show that these minerals are poorly to moderately zoned. Biotites are characterised by a decrease in Fe/(Fe + Mg) ratios when plotted against AIVI. Similarly a modest decrease in Ti from centre to the margins of the biotite crystals has also been observed. Furthermore, an increase in Mg/(Mg + EFe) ratios from centre to the margins as well as in mafic to felsic rocks is also apparent. On the other hand, homblendes from the granitic bodies display less well marked trends and are poorly zoned. However, they show systematic increases in Mg/(Mg + EFe) ratios, and decreases in Fe, Al"' and Fe/(Fe + Mg) ratios, when plotted against Si. Magnetic susceptibility measurements on powdered samples of the plutons give results characterised by high values ranging from 0.8 x 10-4 to 2.0 x 10-3 emu/g which place the granites in magnetite- series field. This indicates that the oxidation state of the magma was high during crystallisation. The range of chemical compositions of the mafic phases indicates a rise in the oxidation state of the magma as it progressively became saturated in H20-volatiles and other metals, not partitioned into the earlier crystallised phases in the final stages of magma consolidation. These variations in chemical compositions also indicate that initially fOz was high and became progressively higher during fractional crystallisation of magmas; this subsequently led to the emanation of hydrothermal fluids responsible for the mineralisation in the adjacent metasediments.

Jan 1, 1994

By S Hapugoda, J R. Manuel, M J. Peterson, A Poliakov, E Donskoi

################

Aug 12, 2013