Search Documents

By D. Shishin

Continuous integrated experimental and thermodynamic modelling work is conducted in Pyrometallurgy Innovation Center (PYROSEARCH). The target is to develop a tool, which is capable to predict the compositions of phases, phase equilibria and enthalpy balance for pyrometallugical processes with higher accuracy. All these can be achieved by thermodynamic calculations based on Gibbs energy minimization algorithm. Recent developments made it possible to predict important factors, such as the chemical solubility of Cu in slag as a function of matte grade, temperature, P(SO2), Fe/SiO2 ratio in slag, wt% Al2O3, CaO and MgO in slag. In industrial practice, it is not only chemical solubility, but the overall amount of slag, which affects the recovery of copper and other valuable metals. It will also impact the energy balance. The presence of solids in the slag matrix influences the physical entrainment of matte in slag. Oxygen enrichment which is sometimes used to control the temperature of the furnace, may also affect the volatilization of minor elements into the gas phase. Same is true if additional fuel is used to compensate for heat balance. In other words, all key process parameters mentioned above are interrelated. Present study demonstrates the ways thermodynamic modelling can be used to link those process parameters on the example of copper smelting. INTRODUCTION Research methodology involving integrated thermodynamic modelling and experimental approach has been recently reported by the authors (Shishin, Hayes, & Jak, 2018 Shishin, Hidayat, Jak, Decterov, & Belov, 2016). Present paper reviews the updates and discusses challenges in the development of a large self- consistent thermodynamic database. The database covers the Al-Ca-Cu-Fe-Mg-O-Pb-S-Si-Zn-(Ag, As, Au, Bi, Ni, Sb, Sn) chemical system. The database works in the environment of FactSage software Bale et al., 2016 and related family of products: ChemApp (Petersen & Hack, 2007), ChemSheet and SimuSage. A list of phases, and thermodynamic models is shown in Table 1. For the copper industry, most important phases are Gas, Slag, Spinel and Liquid matte/metal/speiss. Black-colored elements represent the base system, elements in green can be considered as slag modifiers (slagging elements), and elements in red are minor. Ni is highlighted because it was included recently, in 2019, and currently under development. In addition to copper, the database covers the conditions typical for lead and zinc pyrometallurgical extraction. It is developed in consistency with the experimental results for high concentrations of Pb and Zn. Examples of applications in the present study include the thermodynamic analysis of fluxing for an agitated copper reactor to show the effect of Fe/SiO2, and calculations of copper slag cleaning using reductant.

Jan 1, 2019

By C. Prévost, V. Leclerc, R. Pinto, M. Ruel

"A buffer’s ability to dampen production variations is directly related to its size. Small buffers cannot handle large variations in production, resulting in production losses when their inventory is depleted or saturated. This paper presents a predictive control algorithm used to regulate the inventory in a buffer of limited size, subject to large production variations caused by an upstream batch process. The algorithm was designed to manage the acid flow rates in an acid regeneration process, servicing batch leach reactors at Rio Tinto Fer et Titane, in Sorel-Tracy, Québec. The objective was to increase production by reducing equipment downtime caused by imbalances in the acid tanks servicing the leach reactors. The controller attempts to maintain acceptable acid tank volumes by optimizing the leach reactors acid/slag ratio and reaction time as a function of anticipated acid needs, all while maintaining product quality.INTRODUCTION Some beneficiation processes in extractive metallurgy require processing bulk material in batches. Typical examples of metallurgical beneficiation batch processes are smelting, leaching, electrowinning, elution and carbon-in-pulp. These processes are often situated midstream in an ore beneficiation flow sheet and either receive or transfer ore to continuous processes. Batch processes, by their very nature, introduce large variability in the flow sheet, because they process large finite quantities of material at nearly periodic time intervals. In order to dampen variations in production, various buffers–tanks and stockpiles–are utilized to accumulate ore upstream and downstream of batch processes. A buffer’s ability to dampen production variations is directly related to its size (Phillis & Gumede, 2011). Small buffers cannot handle large production variations, resulting in equipment downtime when their inventory is depleted or saturated (Everett, 2010, Berton et al., 2013). Consequently, the batch equipment production rates must be adequately managed to minimize production losses, attributed to depleted or saturated inventories, while maintaining product quality. Algorithms that anticipate or predict future production variations are often utilized to regulate buffer inventories, subject to large flow rate variations (Salvo et al., 2002; Bonvin, 2013). This paper presents a predictive control algorithm used to maintain an acceptable volume in a spent acid tank servicing several batch leach reactors. The algorithm is implemented at the acid leach plant of Rio Tinto Fer et Titane, in Sorel-Tracy, Québec."

Jan 1, 2017

By Hökmark. H., M. Lönnqvist

"The Swedish Nuclear Fuel and Waste Management Company are planning to construct a permanent storage facility for high-level nuclear waste at about 500 m depth in fractured crystalline rock. One potential concern for the long-term performance of the repository relates to large fracture shear displacements, since these may damage intersected waste containers. Although details in fracture surface geometries are known to influence slip magnitudes and distributions, fractures are commonly approximated as planar features with lab-scale properties due to lack of relevant field-scale data and to limitations in available modelling tools. The purpose of the work described in this paper is to explore the possibilities to use the distinct element code 3DEC to generate and analyze models with isolated, non-planar fractures and, based on the results of such models, get perspectives on the influence of basic disturbances such as ridges, steps and waves. Given the model set-up, fracture segments with orientations deviating from the critical one will contribute to stabilize the fracture. As a result, all analyzed non-planar fractures slip less than the critically oriented planar fracture. We demonstrate that the angle of deviation is the most influential parameter determining the slip magnitude for comparable surface geometries. Similarly inclined steps and ridges have almost equal slip reductions despite differences in surface area restricting the movement. Although it may be possible to translate the influence of such irregularities to an effective coefficient of friction that could be used as input to models with idealized planar fractures, additional analyses with randomly shaped and distributed asperities, different fracture sizes and with alternative stress loads, need to be carried out to devise a defensible strategy for this type of parameter variation. INTRODUCTIONThe Swedish Nuclear Fuel and Waste Management Company’s (SKB) concept for deep geological storage of high-level nuclear waste is the KBS-3 method according to which the waste will be encapsulated in copper canisters and deposited, surrounded by compacted bentonite clay, in vertical deposition holes at about 500 m depth in fractured crystalline rock (SKB 2011). A site has been selected in the Forsmark area, which is located some 120 km north of the Swedish capital Stockholm. One potential concern for the long-term performance of the repository that has been identified by SKB (SKB 2011) is related to large shear displacements along repository host rock fractures induced by a nearby earthquake, since these may cause mechanical damage to intersected waste containers. Due to lack of relevant data at the appropriate field scale and to limitations in available modelling tools, fractures are commonly approximated as planar features with lab-scale properties. The standard numerical tool, 3DEC (Itasca 2013), used in SKB’s rock mechanics analyses, for example, generates planar fractures by default. Results from dynamic 3DEC earthquake-models (Fälth et al. 2010) currently provide input to layout decisions, i.e., the diameter of the smallest planar fractures that must be avoided when deciding on waste container positions at different distances from potential earthquake faults. It is, however, well known that fracture surfaces are rough or undulating on a wide range of scales (e.g., Candela et al. 2012) and that such deviations from planarity influence the mechanical behaviour of the fractures (e.g., Dieterich and Smith 2009, Marshall and Morris 2012, Ritz and Pollard 2012). It will probably never be a realistic option to attempt to represent the geometrical complexity of a population of fractures of different size, shape and surface structure in any large-scale rock mechanical analysis relevant for the performance assessments of the nuclear waste repository. It may however be possible to find approximate ways to represent different styles and degrees of non-planarity by assigning mod"

Jan 1, 2015

By M. Hamamoto, M. Furuta, S. Tanaka

The Tamano smelter and refinery has expanded its production from 228kt/y to 260kt/y since 2006. It was not a big jump, but efficient and economical improvements were executed. These were not only de-bottlenecking and modernization of the smelting plant, including the acid plant, but also additional processes were commissioned, e.g. slag cleaning furnace and steam dryer. On the other hand, a conversion to permanent cathode process in the #2 tank house made the biggest impact on production efficiency for the refinery plant. In order to overcome the drastic change of copper business, although each scale is not so big, continuous improvements are very important to survive. This paper mainly describes the recent improvements of the smelting plant.

Jan 1, 2007

By K. Svens, M. Lahtinen, L. Lehtinen

"Outotec's direct leaching process is an atmospheric leaching process for sulphidiczinc concentrates. The direct leaching process, based on special leaching reactors developed by Outokumpu and leaching know-how, is the most economical way to produce electrolytic zinc. The newest industrial-scale application of Outotec's direct leaching process will be started up in 2008 in China at Zhuzhou smelter. The process isintegrated to the existing production plant, replacing part of the old production facilities by modern technology and increasing the production capacity significantly. The process equipment supplied by Outotec includes also tailor-made reactors designed on the basisof the work and experience of Outotec Research Centre."

Jan 1, 2008

By M. Vaillé, B. Bossé, N. Durupt, A. Michaut

Imouraren is a sedimentary uranium deposit (> 150,000 tU @ 800 ppm), located approximately 100 km North-Northwest of Agadez (Niger); its depth reaches 110 to 170 meters. For low-grades ores, heap leaching is commonly considered as the appropriate technology for metal recovery. More precisely, the ore is crushed (≤ 300 mm) and agglomerated before being stacked. In this case, sulfuric acid can be added during agglomeration and/or irrigation (for three months to a year) to improve uranium recovery. It is known that mineralogical, chemical and physical properties of the uranium ore impact recovery, yet many challenges remain to increase recovery using this process. For example, ore heterogeneity and spatial variability (e.g., variation in uranium, carbonate and/or clay contents) can affect heap leaching performances (i.e., metal recovery, acid consumption or neutralization, percolation rate). However, laboratory column tests using composite ore sample that do not consider ore heterogeneity are widely used to anticipate recovery. Consequently, a new approach that consists in using Imouraren drill core ore samples was recently investigated. This latter also enabled a better assessment of the main parameters that influence uranium recovery and proposed a new predictive model to optimize mining operations.

Jan 1, 2020

By Eugene Barsky

"At present, theoretical concepts of gravitational separation processes are mainly limited by the study of the interaction between a flow and a solitary particle. Then the main conclusions of this analysis are applied to a mass process. However, this analysis is absolutely insufficient, because some elementary phenomena such as the ingress of coarse particles into the fine product output and vice versa cannot be interpreted from its standpoint. It is also impossible to predict the results of the classification of a material with a specific granulometric composition on a specific apparatus on its basis.In the present paper, an extended model of the process is suggested. The model is based on the behavior of all the particles of a specific narrow size class, and not on that of an individual particle. Similarity criteria of the process, which were empirically derived previously, are obtained from the analysis of this model.INTRODUCTIONIn case of gravitational separation, a part of solid phase moves, as a rule, vertically upwards with the flow, while another part settles against the flow. This flow is extremely complicated and chaotic due to various random factors. Among these factors are mass interactions of particles with each other and with the apparatus walls, non-uniform character of concentration and velocity fields, turbulent vortices and local fluctuations, real particles shape difference from spherical one, etc.Particles rotation results in a transverse force (Magnus' effect). At the motion of a non-rotating particle in a non-uniform velocity field, a transverse force also arises (Saffman's effect).A great number of publications have been devoted to the account for all these factors in two-phase flows in recent decades. Very complicated models have been developed, but the problem has not been completely solved as yet [1, 3].At the same time, simplified models of the behavior of a solitary particle in a flow are being developed. Despite their schematic nature, such models often prove to be useful for establishing general laws of the process and its similarity criteria [1]."

Jan 1, 2016

By Irvine J. Berscheid

This paper describes a rock drill bit which, although conventional in the processes of manufacture and use on the job, has incorporated into the outer wall or periphery of the insert, where wear is always the greatest, a variable amount of extra metal, distributed in units of height of insert to withstand this excessive wear. Explained are variations of construction, how the extra metal enables the bit to drill more footage, prevent antitaper and the need to gauge grind, increase penetration rate and the length of run per regrind, reduce shank breakage and ma-chine costs, and permit increased production per man-shift, without requiring any change in equipment already on the job. F OR purposes of clarity, reference to bits here means integral rods or multi-winged bits which have their inserts sharpened to a conventional shape. When dull, the insert has worn, on its face, a flat which increases in width as we progress from near the center of the bit toward the outer edge. Very

Jan 1, 1965

By Robert Hamilton

Aur Resources is Canada's leading copper mining company with producing mines and development projects in North, Central, and South America. Copper production from Louvicourt Mine in Quebec, Quebrada Blanca, and Andacollo Mines in Chile had a record of 245 million pounds in 2001. Located 600 km North of Montreal, Louvicourt Mine is a joint venture between Aur Resources (30%), Teck Cominco (45%), and Novicourt (45%). Aur Resources is the mine operator with 280 employees working to produce a copper and zinc concentrate. In operation since 1994, the mine has produced a total of 12,606,000 metric tonnes of ore grading 3.6 % Cu, 1.5 % Zn, 26gr./m.t Ag, 0.9gr.m.t Au. This highly mechanized mine reached full production in August 1995 at 4,000 tpd. From 1996 to 2001, the average mine production rate was more than 4,300 tpd. The Louvicourt Mine continued its record of solid performance in 2002 and is expected to continue to do so for the balance of its 2.5 years mine life.

May 1, 2003

Thermal barrier coating (TBC) systems provide thermal insulation to high temperature for aero¬propulsion. A fracture mechanical model developed earlier for multilayer systems is applied on TBC to realize superior performance by enhancement of fracture resistance. Functionally gradient materials (FGM) solution concept is the basis of approach followed here. A few realistic FGM TBCs are considered in the analysis in addition to a conventional TBC as the base case and a hypothetical worst case with YSZ-YSZ interface. Parameters like number, composition and thickness of graded layers and elastic modulus are considered for the analysis and quantitative improvement of fracture resistance of TBCs. The analysis is confined to deep crack with crack size to width ratio as 0.5.

Jan 1, 2010

By Forrest A. Kerr

The subject of the Coast Range batholith and its eastern contact zone needs no introduction to British Columbia mining men. Because of its reputed potentialities, this zone has aroused much interest. Because it has touched the pockets of the people by actually producing dividend-paying mines, it has been able, more particularly in certain isolated sections, to retain and accentuate that interest. However, the Stikine River section of the zone has never had the good fortune to have a mine developed, which, through profits or lack of profits, could create such interest. The Cassiar district, beyond the Stikine River area, has been a lure to the prospector for over half a century and was still strong enough in the last few years to attract close to a million dollars in capital and labour. All of this has passed through the Stikine river area, has looked at it, and ignored it. Sufficient interest has never been aroused in the area to bring forth an effort to see what is there, in spite of the fact that it belongs to the same zone in which some of our richest mining camps are located. Because the contact zone of the Coast Range batholith has long been recognized as an ore producer of considerable potentiality, the Geological Survey some years ago formulated a plan to map as much of this area as possible in order to assist in promoting its development. This plan has been furthered during the last fifteen years whenever opportunity presented itself, until today a fairly large part of the area has been covered. The work which has been done on the Stikine during the last two years is the beginning on what is the largest remaining gap in the general plan.

Jan 1, 1928

By K. Robertson, K. A. Scalise, K. C. Kocsis, M. B. Teixeira

"Many deep and highly mechanized underground metal mines require cooling. In many cases, appropriate working conditions might be fulfilled by ventilation alone, but some require the implementation of cooling systems. The most common cooling systems are bulk air cooling and localized cooling. These systems are effective, but they use considerable amounts of electricity. Ground source heat pumps (GSHP) can supplement existing cooling and heating systems, while geothermal systems can cool the mines in their entirety if there are significant sources of geothermal energy are present. A major advantage of geothermal and GSHP systems are their renewability, which reduces the carbon footprint of the mine.INTRODUCTIONThe western United States has abundant alternative energy sources including a high potential for geothermal energy. Geothermal utilization is a viable energy alternative that many areas throughout the world have access to because high geothermal potentials are located near areas of volcanic activity such as along plate boundaries, mid-oceanic ridges, rift valleys and hot spots (IRENA, 2017). There is significant potential to use geothermal energy for many applications. With constantly evolving and innovative technologies being developed, an increasing number of geothermal applications are becoming feasible. The western United States is an ideal candidate for alternative energy applications due to the geothermal availability and the emergence of modern technology. This paper investigates the incorporation of renewable methods and the alternative energy possibilities of Barrick Goldstrike, located in northern Eureka County, NV and Rio Tinto/BHP Billiton’s Resolution Copper mine located in Superior, AZ.RENEWABLE SYSTEMSGeothermal SystemsGeothermal temperatures are typically defined as either high or low temperature systems(Shengjun, Huaixian, & Tao, 2011). The high-temperature reservoirs (>220°C / 428°F) are feasible forlarge, commercial production of electricity with dry stream and flash stream systems (Shengjun, Huaixian, & Tao, 2011). Low temperature systems, below 220°C (428°F) (Shengjun, Huaixian, & Tao, 2011), aremuch more difficult to use for energy production. Low-temperature reservoir heat can be harnessedthrough the implementation of a binary cycle plant. The working fluid for a binary cycle plant is eitherorganic (Organic Rankine Cycle (ORC)) or an ammonia-water mixture (Kalina cycle). Organic fluids andammonia-water mixtures have low boiling points that allow for the recovery of heat from low-temperaturesources (Eymel Campos Rodriguez, et al., 2013). The efficiency of the binary cycle plant varies depending upon the fluid, thus the proper selection of the fluid is critical for maximum energy output"

Jan 1, 2019

By G. C. Nowlan, A. Illis, H. J. Koehler

At its iron ore recovery plant at Copper Cliff, INCO produces about 17 million pounds per year of a pure, acid-soluble, dense nickel oxide from nickeliferous pyrrhotite. The process used is also applicable to a broad range of ammoniacal nickel-bearing streams. The product has been widely accepted by the chemical, ceramic and electronic industries.

Jan 1, 1970

By R. W. Lawrence

"The management of acid rock drainage( ARD) to prevent impact when tailings are disposed of in the environment requires understanding of the phenomena giving rise to acidic drainage, a rigorous evaluation of the weathering characteristics of the tailings, and the successful implementation of control measures that will be effective in the long term. This paper provides a brief theoretical background to the complex reactions taking place in mill tailings containing residual sulfide minerals which must be appreciated to enable waste characterization program to be effective in predicting the future behavior of the tailings and in designing control methods which will work. The principals and practice of prediction and control as part of an ARD management strategy will then be discussed.Paper Not Available at Time of Publication"

Jan 1, 1997

By I. C. Brown, J. F. Henderson

"IntroductionThe Yellowknife greenstone belt lies along the west side of Yellowknife bay on the north shore of Great Slave lake, some 600 miles north of Edmonton, Alberta. The producing gold mines along the belt include the Con-Rycon and Negus which came into production in 1938 and 1939, and the Giant Yellowknife, which followed in 1948. In addition to these, the Akaitcho and Crestaurum mines are potential producers and there are several promising prospects.Early geological work along the belt was done by A. W. Jolliffe who first recognized a number of great faults, the largest of which is the West Bay fault with a horizontal movement of more than 3 miles. In 1943, A. S. Dadson, following an examination of the Giant-Yellowknife property, advanced the hypothesis that Baker Creek valley was underlain by a series of major gold-bearing shear zones. Subsequent diamond drilling in 1944 bore out the correctness of this view and outlined the great ore-bearing shear zones that are now being developed and mined by Giant Yellowknife Gold Mines. Near the south boundary of the Giant property, the shear zones are cut off by the West Bay fault. This prompted Neil Campbell, District Geologist for the Consolidated Mining and Smelting Company, to make a detailed study of the geology along the West Bay fault in the hope of finding the faulted extension of the Giant shear-zone system on the west side of this fault. In this he was entirely successful. By careful geological mapping he matched several faulted diabase dykes on both sides of the fault and, knowing their dip, determined the horizontal and vertical components of movement along the fault. Subsequent •diamond drilling based on Campbell's work located the presumed correlative of the Giant shear-zone system on the Con-Rycon and Negus properties, where it has been named the Campbell zone or system. At the present time, all production at the Negus Mines property is coming from orebodies along the Campbell zone."

Jan 1, 1950

By G. W. H. Norman

"During the past few years emphasis has been placed on the importance of major faults, fault zones, and associated structures in localizing ore deposition. This interest has developed partly through experience in mining camps in other parts of the world, and partly by the discovery of such structures associated with ore deposits in the Abitibi region, northwestern Quebec. The recognition of these structures has been brought about chiefly by detailed mapping that has been gradually extended along most of the southern gold belt of Quebec and across adjacent parts of Ontario. In the Chibougamau region of northern Abitibi, progress has been made in mapping major faults by using aeroplane photographs with a scale of about 1,000 feet to 1 inch for field mapping. These studies have shown that many large faults occur in the Abitibi region; also that there is a relationship, direct or indirect, between faults and mineral deposition. The faults are, however, of more than one type and age, and this should be clearly appreciated in any planned programme for prospecting this great region adequately.Geology of the Abitibi RegionGeneral GeologyFaults of more than one age and type occur in the Abitibi region. In order to appreciate some of the age relationships of the faults, a brief summary of the historical aspects and major structural features of the geology (fig. 1) is necessary. These historical aspects may be briefly summarized in chronological order as follows:(1) Accumulation of Archæan lavas and associated pyroclastic rocks, either alternating with, or essentially followed by, deposition of sedimentary rocks, and interrupted at least once by local and perhaps regional folding and erosion."

Jan 1, 1954

By Shoma Sinha, Nirmal K. Sinha

Deformation in a material at any temperature is known to consist of two major components - a reversible (elastic) part and an irreversible part (often called plastic). What is not universally known is the fact that at high temperatures greater than about 0.4Tm, where Tm, is the melting point in Kelvin, the reversible part includes not only elastic but also a delayed elastic (recoverable but time dependent) component. A very simple computational model is presented here for a constant-strain stress relaxation test (SRT) on the basis of a three-component constant-stress rheological equation derived from strain relaxation and recovery test (SRRT) data. It is shown that the rapid stress reduction during the early periods of relaxation when the stresses are high is governed primarily by delayed elasticity irrespective of the imposed strain level. This is in contradiction to the popular concept prevailing today. The asymptotic approach to a quasi-stable stress rate at longer times and lower stresses is, however, controlled by viscous flow (power-law creep). The validity of the theory is examined by comparing the predictions with experimental results obtained for a titanium-base alloy Ti-6246 at 0.45Tm.

Jan 1, 2004

By L. E. Djingheuzian

"ForewordTHE FOLLOWING excerpt is taken from Fine Grinding Investigations at Lake Shore Mines (C.I.M., Trans., Vol. XLIII, 1940, p. 307):""MacCalman, who has carried out exhaustive testing on the accuracy of screen analysis, makes the following statement in his report, The Accuracy of Sieving Tests, Industrial Chemistry (London, England), Vol. XIV, p. 231:"" 'It will be shown that, as no specific relationship between the sieving value and the various dimensional features of the sieve is evident, though there is an apparent tendency for the sieving value to follow the average aperture size, it is apparently impossible to define the sieving value of a screen from dimensions alone, and the calibration of sieves for important tests on a specific material, whereby the sieving values of these sieves are related to that of an arbitrarily chosen master sieve by correction values, will be recommended."" 'It will be emphasized that these correction values are probably not applicable to other materials, or the same material of different size composition. Sieve users, for important tests or in cases of dispute, will be recommended to calibrate their sieves against a mutually acceptable standard with specified conditions of test, and to recalibrate at intervals' ""."

Jan 1, 1953

By Hani S. Mitri

Overhand cut-and-fill mining of narrow vein orebodies inevitably leads to the creation of horizontal pillars. As the thickness of the overlying crown pillar diminishes, stresses build up in the pillar. More often than not, this situation causes an increased level of seismicity and possible pillar burst. One well-known strategy to reduce the hazard of crown and sill pillar burst in cut-and-fill vein mining is that of destress blasting. The paper first reviews the methods used for the monitoring and assessment of potential pillar burst in narrow vein cut-and-fill mining. The influence of horizontal stress, mine depth and pillar thickness on stress build-up in the crown pillar is examined through a characteristic narrow vein numerical model. The paper then presents, with reference to three mine case studies, the benefits of destress blasting and the factors affecting its success as a ground control technique.

Oct 1, 2001

By Marc Symbaluk

Core Programs Current ? Aboriginal Involvement Program ? Adaptive Forest Management / Forest History Program ? Grizzly Bear Research Program ? Natural Disturbance Program ? Foothills Growth and Yield Association ? Foothills Stream Crossing Program ? Fish and Watershed Program ? Caribou Landscape Management Association ? Social Science Program ? Local Level Indicators Program ? Communications and Extension Program ? Geographic Information Systems Program ? Mountain Pine Beetle Ecology Program In Development ? Climate Change Program ? Forest Community Program ? Water, Quantity and Quality Program (Water for Life) ? Circumboreal Initiative ? Yellowhead Ecosystem Management Group

May 1, 2008