Search Documents

By I. C. Spangenberg

The status of sampling practices in the gold mining industry in Africa was determined as an initial step in a process to standardize sampling practices in the mining industry. Twenty-one gold mines, twenty metallurgical plants, and thirteen laboratories were rated for the potential influence of the relevant sampling errors on each component of the particular sampling system. The findings of the study on the status of equipment, standards, procedures, and management principles were presented in a dissertation (Spangenberg, 20121) and are now used by AngloGold Ashanti?s mines in the quest for correct sampling practices. The information is also used to compile a guideline on leading practice procedures for sampling methods in gold mining. This article is an overview of sampling best practice found in African mining Keywords theory of sampling, best practise, sampling errors, sampling standards.

Jan 1, 2014

By B. P. Watson

The Bushveld platinum group metal deposits are two distinct, shallow-dipping stratiform tabular orebodies and strike for many hundreds of kilometres. Mining is extensive, with depths ranging from close-to-surface to 2 000 m. The mining method is a variation of planar open stoping. Pillars are widely used to support the openstopes. Little work has been done in the past to determine pillar strength and pillars have been designed using experience and formulae developed for other hard-rock mines. This has led to oversize pillars with consequential loss of ore or undersize pillars resulting in potential collapses or pillar bursts. This paper describes a back-analysis of pillar failures at Impala Platinum, just north of Rustenburg, and provides insights into pillar behaviour. Strength formulae are presented for Merensky Reef pillars at Impala. These formulae may be used cautiously on other platinum mines with similar geotechnical and geomechanical conditions.

Jan 1, 2008

By Q. G. Reynolds

The presence of dust and fume suspended in the freeboard region is a common feature of the operation of direct current (DC) plasma smelting furnaces. This occurs primarily as a result of the use of fine feed materials together with the open-arc, open-bath operation of such smelters, and is exacerbated by the high velocities and turbulent mixing of the gas in the vicinity of the arc jet. Dust and fume losses into the furnace off-gas system can be significant in some cases and may have economic, operational, and environmental impacts on the process. A computational modelling study is presented in which the concentration of dust material was considered as a continuous field subject to a governing partial differential equation. Settling behaviour was calculated as a function of particle size, local gas/plasma temperature, and other physical properties. Development of the coupling between the concentration field and a magnetohydrodynamic description of the arc is shown, and the resulting models were used to compute various aspects of the behaviour of the concentration field in the arc region for a variety of furnace conditions. Time-averaged as well as transient models of the arc were used to generate the results presented. Qualitative case studies produced several practical suggestions for furnace operation, including increased dust capture by the bath when feed ports are located closer to the electrode, and the possible effects of feed segregation in the furnace freeboard based on dust particle size and density.

Jan 1, 2015

Cokriging uses a sparsely sampled, but accurate and precise primary dataset, together with a more abundant secondary data-set, for example grades in a polymetallic orebody, containing both error and bias, to provide improved results compared to estimation with the primary data alone, as well as filtering the error and mitigating the effects of conditional bias. The method described here may also be applied in polymetallic orebodies and in other cases where the primary and secondary data could be collocated, and one of the data-sets need not be biased, unreliable, etc. An artificially created reference data-set of 512 lognormally distributed precious metal grades sampled at 25×25 m intervals constitutes the primary data-set. A secondary data-set on a 10×10 m grid comprising 3200 samples drawn from the reference data-set includes 30 per cent error and 1.5 multiplicative bias on each measurement. The primary and secondary non-collocated data-sets are statistically described and compared to the reference data-set. Variograms based on the primary data-set are modelled and used in the kriging of 10×10 m blocks using the 25×25 m and 50×50 m data grids for comparison against the results of the cokriged estimation. A linear model of coregionalization (LMC) is established using the primary and secondary data-sets and cokriging using both data-sets is shown to be a significant improvement over kriging with the primary data-set alone. The effects of the error and bias are filtered and removed during the cokriging estimation procedure. Thus cokriging using the more abundant secondary data, even though it contains error and bias, significantly improves the estimation of recoverable reserves.

Jan 1, 2014

By F. Lombard

We consider the Grubbs estimation of the precision of an online gauge. Typically, this type of estimation involves independent results from two or more reference instruments (sampling and laboratory analysis operations). The properties of the estimator are then independent of the product variability. However, the use of more than one reference instrument entails significant additional costs. The two-instrument Grubbs estimator, which is based on results from a gauge and a single reference instrument, has the disadvantage that its standard error is heavily dependent on the degree of product variability. We propose a new estimator that has a variance that is more or less independent of product variability. In fact, the variance is typically less than that of the Grubbs estimator based on the use of two reference instruments. In order to function successfully, our methodology requires some prior knowledge of the extent of product variability and of gauge precision. In practice, such prior knowledge is often available and it is a weakness of the traditional two - and three-instrument Grubbs methods that no use is made of such knowledge. The efficacy and robustness of the new method is illustrated by Monte Carlo simulation.

Jan 1, 2012

By P. Venter, N. Naude

"Coal mining is a thriving industry and 53% of the coal mined in South Africa is used for electricity generation. Mechanization has made coal mining more efficient, but fines generation has subsequently increased. Up to 6% of the run of mine material can report to the -200 µm fraction. Common problems associated with fines handling include dust formation, storage problems, and high moisture levels. A method to turn this material into a saleable product instead of stockpiling it can add value to a company.Briquetting is a pressure agglomeration method where loose material is compacted into a dense mass (FEECO International, 2014). The briquettes must be able to withstand rigorous handling and transport operations without disintegrating. This study aims to investigate the optimum binder and moisture conditions required to produce a mechanically strong briquette using two different binders - a PVA powder (binder A) and a starch powder (binder B).It was found that for binder A the optimum moisture level was 12% to 14%. At this moisture level the greatest compression strength gains were observed, and low amounts of fines produced in impact and abrasion tests. The minimum amount of binder added while still obtaining a strong briquette was 0.5% binder A. For binder B the optimum moisture level was also 12% and the minimum amount of Binder B to be added was found to be 1%. Briquettes that were dried outside reached their peak strength after about four days, whereas the briquettes that dried inside took about 20 days to reach their strength plateau. Hardly any degradation took place on the surface of the binder A film after exposure of 300 hours of artificial weathering. Thermogravimetric analysis confirmed that neither binder A nor binder B will add to the ash content of the coal fines, as both binders totally decompose above 530°C.Binder B yielded stronger briquettes after 15 days and also generated less fines. It is therefore superior to binder A and would be recommended for further use.Paper written on final year project work carried out in partial fulfilment of B.Eng degree in Coal Beneficiation"

Jan 1, 2015

Many industries base their decisions on the assaying of tiny analytical sub-samples. The problem is that most of the time several sampling and sub-sampling stages are required before the laboratory provides its ultimate assays using advanced chemical and physical methods of analysis. As long as each sampling and sub-sampling stage is the object of due diligence using the theory of sampling it is likely that the integrity of the sought after information has not been altered and the generated database is still capable to fulfil its informative mission. Unfortunately, more often than not, unawareness of the basic properties of heterogeneous materials combined with the unawareness of stringent requirements listed in the theory of sampling, lead to the conclusion that massive discrepancies may be observed between the expensive outcome of along chain of sampling and analytical custody, and reality. There are no areas that are more vulnerable to such misfortune than sampling and assaying for trace amounts of constituents of interest in the environment, in high purity materials, in precious metals exploration, food chain, chemicals, and pharmaceutical products. Without the preventive suggestions of the theory of sampling serious difficulties may arise when making Gaussian approximations or even lognormal manipulations in the subsequent interpretations. A complementary understanding of Poisson processes injected in the theory of sampling may greatly help the practitioner understand structural sampling problems and prevent unfortunate mistakes from being repeated over and over until a crisis is reached. This paper presents an overview of the theoretical, practical and economic difficulties often vastly underestimated in the search for quantifying trace amounts of valuable or unwelcome components.

Jan 1, 2010

By A. v. E. Conradie, I. Nieuwoudt

Comminution and flotation circuits are frequently perturbed by disturbances, such as feed particle size variations, ore hardness and varying mineral liberation, which can have a significant impact on the recoveries and grades. Although advanced control strategies can neutralize many of these disturbances, few such systems have been implemented at present. One reason is that it is only recently that the requisite sensor systems have been developed and used to control froth flotation systems. Another is owed to the complexity of these non-linear interactive systems. In this paper, a plant-wide control strategy is considered, wherein both unit operations have a single overall control objective. The automated design of the controller is based on an advanced evolutionary algorithm viz., Symbiotic, Memetic Neuro-Evolution (SMNE), for developing neural network controllers in a reinforcement learning framework. A rigorous phenomenological comminution and flotation model was used to develop suitable neurocontrollers, which used all possible sensor inputs and available manipulated variables. The neurocontroller managed to identify optimal or near-optimal operating conditions consistently and showed robust performance in the presence of significant disturbances. The results suggest that evolutionary reinforcement learning offers a very competitive approach towards the development of integrated non-linear controllers for mineral processing plants. Keywords: neural networks, control, flotation, comminution, evolutionary reinforcement learning

Jan 1, 2003

By J. H. Venter

This paper contributes towards the solution of the problem of characterizing the accuracy of estimation of ore quality based on a random sample from the ore bulk.

By M. Stephton, L. A. Cornish, R. Süss, N. B. Maledi, J. H. Potgieter, L. Chown

In today?s gas turbines, the emphasis is on saving energy and reducing the amounts of pollutants emitted. This can be attained only by designing alloys with higher melting points and the capability to retain mechanical integrity at increased temperatures. The development of nickel-based superalloys (NBSAs) has reached its limit, as the alloys are operating at critical temperatures close to their melting points. Further increase in the operating temperatures will result in dissolution of the strengthening phases and even melting. Recent trends are towards developing materials that are adequately strong but have much higher melting points and corrosion resistance than the NBSAs that are currently in use. Studies conducted on Pt alloys with similar structure to the NBSAs have presented mechanical results that are promising. This has led to further development of these alloys as a replacement for NBSAs. The operations of gas turbines are usually limited by hot corrosion, which could lead to catastrophic failures. Increased temperatures have an adverse effect on the corrosion of the alloys by making the environment more aggressive. Hence, during a material selection process, the corrosion behaviour of any selected material for a highly specialized application such as gas turbines needs to be assessed. Generally, coatings are applied to combat degradation by salt deposits, hence coated and uncoated NBSAs were selected as benchmark alloys. Corrosion mechanisms and kinetic studies of the Pt alloys were conducted. High temperature hot corrosion (HTHC) tests were conducted at 950°C for540 hours. Analysis of the corrosion products was conducted using a Raman spectrometer, scanning electron microscope with energy dispersive X-ray spectroscopy (EDS) optical microscopy, and X-ray diffraction study of the phases that formed. Both coated and uncoated NBSAs failed under the test conditions. Samples were covered with a porous non-adherent greenish layer which is characteristic of NiO. Pt alloys displayed superior corrosion properties compared with the benchmark alloys.

Jan 1, 2006

By R. J. Davidson

The chemistry related to the adsorption of gold, uranium, and transition-metal ions onto algal biomass (activated carbon) is related to the genesis of the Witwatersrand. Detrital gold being washed into the early basin, together with cyanide generated in the atmosphere, would have solubilized as the stable aurocyanide complex, the necessary oxygen being produced from the photosynthesis of carbon dioxide by primitive algae. Pebble-type conglomerates, while providing for the free movement of solutions containing the vital nutrients for prolific algal growth, would have screened of the sunlight, which would otherwise have degraded any free cyanide in solution. With the subsequent decomposition of the algal deposits, it is surmised that carbon-rich layers having adsorptive properties similar to those of activated carbon formed in the conglomerates. Under these conditions, gold (silver) in solution would have been adsorbed selectively as the cyanide complex, together with uranium as the carbonate complex. Any transition-metal ions present would also have been adsorbed onto such a substrate. The subsequent burial and compression of the gold-rich conglomerate with temperatures rising to about 400DC would then have reduced the adsorbed gold to the metal in a single segregated gold-silver metal phase (95 to 99 per cent of the contained gold), while leaving behind a sulphide phase rich in iron, base metals, and arsenic and containing some gold and silver. The recrystalization of this later phase may explain the presence of the rather consistent contents of secondary pyrite in the Witwatersrand conglomerates. Similarly, uranium minerals, some rich in titanium and carbon, may also have formed. An adsorption model would explain the very consistent trends in the gold-to-silver ratios of individual reefs in the Witwatersrand, which suggest an extensive hydrothermal system approaching isothermal equilibrium. Also, as gold grades increase, so silver grades generally decrease, indicating the sequential displacement of silver by gold as classically obtained with activated carbon.

Jan 1, 1990

By Ishwinder Grewal, Michael Fullam

Many factors affect the overall performance of a gravity circuit beyond that of the primary gravity recovery device. It is critical to examine the gravity circuit in its entirety in order to determine the dominant operating variables that affect overall performance. Knelson has developed a mathematical model which has been used to assist in the design of gravity circuits, and in particular, in making accurate predictions of plant scale gravity recovery. More recently, inclusion of recovery by particle size class and additional data from various mine sites and research studies has made the model more robust. The model is now able to predict GRG (gravity-recoverable-gold) recoveries using various devices such as centrifugal concentrators, shaking tables, jigs and flash-flotation within the grinding circuit. The benefits derived by these recent advancements in modeling capability are two-fold; the first is as a tool to assist in circuit design and correct equipment selection and the second is as an analytical tool to help to evaluate and enhance the performance of existing gravity circuits.

Jan 1, 2006

By C. Steyn, J. D. le Roux

A dynamic model is used to represent the single-stage closed grinding mill circuit at Mogalakwena North Concentrator. The circuit consists of a ball mill closed by a hydrocyclone. The circuit is modelled according to a semi-mechanistic non-linear dynamic model available in literature. An algebraic routine is used to fit model parameters to steady-state process data. Preliminary validation results against dynamic data indicate that the model is capable of capturing the main dynamic trends of the plant. Since classical system identification techniques proved insufficient to produce accurate models; the non-linear model can be used as a reliable source for model-predictive controller design. Keywords: Grinding, modelling, process control, parameter estimation

Jan 1, 2020

By A. M. Neale

Kriel Colliery has experienced some sub-standard blast results on the overburden. Blast results are considered poor when the fragmentation is too big, back break into the new highwall, capping is experienced and/or excess overburden material is thrown into the void. Poor blasting leads to lower productivity, equipment breakdowns, and poor drilling and blasting results on the adjacent blast, contamination of coal or loss of coal. A simulation was developed in Microsoft Excel to model blasting results while changing various parameters. From the simulation results, it could be seen that a change of blast design could be beneficial. As an example, smaller drilling patterns were simulated (10 m x10 m to 8 m x8 m, and 8 m x8 m to 7 m x7 m in Pit 5 and Pit 6 respectively) in combination with changes in explosives (P700 emulsion vs. ANFO) and drill bit diameter (311 mm vs. 251 mm). The results show that each of the options has a specific scenario that will be most beneficial to the mine. These results are quite specific and subject to variable inputs. The purpose of the example is to demonstrate how the simulation can be used as a tool to asses blast design variations.

Jan 1, 2010

By M. A. Hermanus

This paper examines the occupational health and safety performance (OHS) of the South African mining sector against the backdrop of changes in the composition of the sector, international trends in OHS performance, and the agreement on OHS milestones and targets by mining stakeholders at the Mine Health and Safety Summit of 2003. Although OHS performance has improved, progress is slow and there is need for significant rather than incremental change if the targets are to be realized. Since 2003, fatalities and injuries are 20?25% short of annual milestones and not all commodities show consistent improvement. Statistics on occupational health impacts are unavailable, but noise and respirable dust levels are known to be high in mines, with much work lying ahead to address exposures at source. Trends in regulating and addressing OHS hazards are discussed with reference to the notion of ?systems-thinking? and the targets. Efforts to improve OHS and to respond to changes in the sector are constrained by a lack of training and consistency in risk management, guidance for junior, small and artisanal miners, and holistic approaches to risk. For example, treating risk holistically would involve taking proper account of contractors and women in the workplace, and attending to human factors and ergonomics.

Jan 1, 2007

By A. Tosun

"The size distribution of the muckpile formed as a result of open pit blasting operations has a considerable effect on the efficiency of loading, hauling, and crushing. Various researchers have developed specialized computer software that uses image analysis methods for determining the size distribution of the muckpile. However, these methods have some limitations. One of the most important of these limitations is that the very fine fragments in the muckpile cannot be used in the size distribution calculation. In this study, 18 test blasts in total were carried out in two limestone quarries belonging to Batiçim Corp. in Izmir, Turkey. A new model was developed in order to ensure that very fine fragments are used in the size distribution calculation. The size distributions of the test blasts were calculated by both the Wipfrag computer program and the new model. Correlations were established between the muckpile size distributions determined by both methods and the parameters determining the efficiency of the loader. IntroductionThe size distribution of the muckpile formed as a result of open pit blasting operations has a considerable effect on the efficiency of loading, hauling, and crushing. Various researchers have established certain correlations for predicting the efficiency of the operations according to the muckpile size distribution so that production can be carried out economically (Tunstall and Bearman, 1997; Nielsen and Kristiansen, 1996; Workman and Eloranta, 2004; Michaud and Blanchet, 1996; Frimpong, Kabongo, and Davies, 1996; Osanloo and Hekmat, 2005; Molotilov et al., 2010). In these studies, the 50% passing size (X50) is generally used. These methods rely on the correct measurement of the size distribution of the muckpile formed by the blast. The most reliable measurement method would be to subject the entire muckpile to a sieve analysis, but this is obviously impracticable. Therefore, computational method based on image analysis have been developed. These include IPACS (Dahlhielm, 1996), Tucips (Havermann and Vogt, 1996), Fragscan (Schleifer and Tessier, 1996), Cias (Downs and Kettunen, 1996), GoldSize (Kleine annd Cameron, 1996), Wipfrag (Maerz, Palangio, and Franklin, 1996), Split Desktop (Kemeny, 1994), PowerSieve (Chung and Noy, 1996), and Fragalyst (Raina et al., 2002). In these methods, the size limits of the fragments forming the material are determined by image analysis on photographs of the muckpile. However, the computer software has some limitations. It does not take into account the third dimension of the muckpile, nor the size distribution of very fine fragments. In particular, the fragment size ranges below 2.5 to 3 cm cannot be determined, and therefore these cannot be included in the muckpile size distribution. However, the muckpile formed by a bench blast contains a wide range of particle sizes (Tosun et al., 2015)."

Oct 1, 2018

By D. Minson, G. Lindner, S. Miskovic

Thickening is an essential process in the mining industry and currently one of the most costeffective technologies for recovering water from tailings streams. Improved performance of thickeners can lead to significant energy savings, and increased recovery of water close to the mill. Its performance is directly affected by the rheological properties of slurry at the raking zone. To further understand the dynamics of slurry inside thickeners, a multiphase particle-in-cell method (MP-PIC), is used to investigate its behaviour on a sliding channel. Parameters like particle-size distribution (PSD), rheological properties, and domain discretisation were benchmarked to study its influence on highly dense flows. The slurry concentration by weight of the samples ranged between 60 to 70% of solids by mass, under which conditions it is known to behave like a non-newtonian fluid. An experimental flume, 0.2 m wide and 2 m long, was used for validation under different angles and slurry composition. Visual analysis of the slurry height and shape, as well as velocity measurements, were compared with numerical data. The preliminary results show the ability of MP-PIC method to represent different types of highly dense flows, while still being computationally affordable and practical for industrial use. Keywords: High density thickener, CFD, MP-PIC, multiphase-flow

Jan 1, 2020

By M. Bester, R. Carey, T. Russell, J. van Heerden

"In order to reach the goal of optimal and safe extraction of the mineral reserve, the mining industry is increasingly focusing on reconciliation across the mining value chain to ensure that the mining process occurs in a progressively more predictable manner. This will consequently improve the understanding of mine plans and their associated implementation risks, ultimately leading to increased investor confidence. In terms of geotechnical engineering, mine to design performance assessment is a critical aspect of mine to plan adherence focusing on the complex aspects influencing slope performance. A mine to design reconciliation process therefore needs to systematically track implementation and geotechnical aspects in order to review or update slope design criteria and pit layouts according to the mine planning cycle as the mine is developed. This paper describes the development and implementation of a standardized mine to design reconciliation system for operations utilizing laser scanning technology, analysis tools, as well as associated key performance indicators. Performance is evaluated on final pit boundaries in terms of geometry achieved as well as blast performance evaluation of face conditions. These aspects as well as the resulting catch berm capacity should be captured in an information system, enabling longer term trend analysis and detailed root cause analysis where unsatisfactory outcomes occur. Adverse outcomes should subsequently be incorporated into the mine’s risk management process informing risk mitigation measures including blast design adjustments, slope design changes, or rockfall protection measures. IntroductionIn order to reach the goal of optimal and safe extraction of the mineral reserve, the mining industry is increasingly focusing on reconciliation across the mining value chain to ensure that the mining process occurs in a progressively more predictable manner. This will consequently improve the understanding of mine plans and their associated implementation risks, ultimately leading to increased investor confidence. In terms of geotechnical engineering, mine to design performance assessment is a critical aspect of mine to plan adherence focusing on the complex aspects influencing slope performance. A mine to design reconciliation process therefore needs to systematically track implementation and geotechnical aspects in order to review or update slope design criteria and pit layouts according to the mine planning cycle as the mine is developed."

Jan 1, 2016

By A. Marek

The safe and optimal installation of roof support remains one of the biggest challenges in underground mining today. Traditional hand-held, roofbolt installation in conventional mining has been replaced, in a lot of mining operations, by mechanized equipment. This was a first step to moving the operators out of the unsafe working environment and a start on the trend toward a changed support regime in mechanized mining. However, the challenge of handling of the drill steel itself, as well as inserting resin cartridges and bolts, remains a major hazard and the cause of many lost-time injuries. The paper describes and discusses developments in this area, of mechanized bolting and on-board automation, presenting solutions to known problems, achievements to date, and future developments. The solutions presented will identify with addressing customer needs while touching on possibilities for further technology changes, keeping focus on the safety aspect for operators and the need to bring the operators even further away from the unsafe working environment. The-state-of-the-art technology in mechanized bolting equipment is presented and a brief outlook on the past, current, and future developments is given.

Jan 1, 2012

By K. Basaure, S. Palma, R. L. Castro

"Mud rushes are sudden inflows of fine wet ore within underground mines; they can harm people and equipment as well as cause production delays. Block and panel caving mines are prone to mud rushes in their production levels, with mud flowing from the drawpoints. Different mines have developed methods to control drawpoint conditions and thus perform safe ore extraction. The aim of this study is to develop a decision matrix and apply it to samples from drawpoints for the Diablo Regimiento mine, Codelco, El Teniente Division (Chile). The Abrams’s Cone was used to perform slump tests on mud samples and characterize their consistency by setting the saturation degree and packing. The unconfined strength for different saturation degrees and packing was then evaluated. The results show that conditions for fluid response in slump tests correspond to conditions for low strength in unconfined compression tests. We also found the water content value at which the flowing consistency changes from plastic to soft behaviour according to the slump test classification for each sample. This value varied between 12.2% and 16.9%, depending on the sample tested. Finally, a fluid security factor was defined in order to design an extraction pattern for the samples tested, based on the water content. We conclude that a very densely packed ore is not prone to flowing. The flow properties depend strongly on the specific properties of the ore tested, and consequently a specific criterion should be developed for different types of ore at a particular mine. The criterion developed for Diablo Regimiento ore provides critical water content values from 11% to 15%, depending on the ore type. IntroductionBlock/panel caving has the lowest cost and the highest production rates of all underground methods (Heslop, 2000; Hubert et al., 2000). Good management of the risks is important to ensure the security of personnel and infrastructure. Operational risks for the method include rockbursts, air blasts, and mud rushes (Hubert et al. 2000). A mud rush is defined as a sudden inflow of mud from a drawpoint or other underground opening into the mine (Butcher, Joughin, and Stacey, 2000; Butcher, Stacey, and Joughin, 2005). Mud rushes can cause dilution, production delays, and damage to equipment as well as injuries and fatalities. The elements necessary for mud rushes are mud-forming material, disturbances, a discharge point, and water (Butcher, Joughin, and Stacey, 2000). All these elements are present in block/panel caving mines. Mud-forming material comes from ore milled by the caving process; water comes from underground and surface sources. Drawpoints in the production level are discharge points for the mud. Finally, several types of disturbance could trigger the mud rush phenomena: blasting, earthquakes, operation of equipment (Call & Nicholas, Freeport McMoRan, and Hydrologic Consultants, 1998; Jakubec, Clayton, and Guest, 2012), arch collapse above drawbells (Butcher, Joughin, and Stacey, 2000; Butcher, Stacey, and Joughin, 2005; Jakubec, Clayton, and Guest, 2012) and ore drawing (Butcher, Joughin, and Stacey, 2000; Butcher, Stacey, and Joughin, 2005; Call & Nicholas, Freeport McMoRan, and Hydrologic Consultants, 1998; Jakubec, Clayton, and Guest, 2012)."

Jan 1, 2017