Search Documents

By I Canbulat, A Mcquillan

This paper presents a parametric analysis showing the impact of virgin in situ stresses on a shallow (less than 100 to 150 m) high slope in a jointed sedimentary rock environment. Slope stability is assessed at varying ratios of slope orientation and horizontal stresses. The results of this study indicate stress conditions should be included in stability assessments of excavated slopes, particularly where structure is orientated perpendicular to the direction of horizontal stress. Some consideration of stress state is applied in deep open cut pits (eg those with planned excavations of 400–1200 m) (Myrvang et al, 1993; Stacey et al, 2003; de Bruyn et al, 2014; Kozyrev et al, 2015). Dodd and Anderson (1972), Kalkani and Piteau (1976), Lee (1978), Coulthard et al (1992), Noorani et al (2011) have studied the effect or horizontal stress and its impact on slope stability, showing that tensile stresses tend to develop at the crest of the slope, with larger tensile zones for higher horizonal stresses and steeper slope angles (Stacey, 1970, 1973; Sjoberg, 2013). However, the consideration of virgin stresses on shallow slope stability is less frequent (Lynch et al, 2005; Dight, 2006; Lucas, 2006; Sjoberg, 2013). This has led many geotechnical engineers to believe the impact of virgin in situ stresses on open cut slope stability is minimal and not worth time or consideration in modelling, or their effects on slope stability are still poorly understood (Stacey et al, 2003; Sjoberg, 2013).

Nov 29, 2022

Contemporary approaches to safety management appear to be failing short of meeting its mark in improving mine safety. This is evidenced by the high workers compensation, high incidence rates and fatalities. Evidence from high-risk and complex organisations points towards safety culture as being important in improving site safety. In more recent years resilience engineering has been touted as a new and innovative way of managing safety. This paper reviews and synthesises previous literature on safety culture and resilience engineering. It then highlights methods that can be used to measure safety culture and resilience engineering, and explores similarities and differences between complex organisations and gold mining to identify opportunities for more innovative approaches to improving safety in gold mining operations through safety culture and resilience engineering.

Aug 1, 2010

By B J. Wilde

Evaluating the resources and reserves in a mineral deposit requires many interdependent subjective decisions. The decision of stationarity is arguably the most important. This decision includes: the choice of domains that subset the entire deposit, modelling the boundaries of each domain, and understanding and considering the nature of the contacts between domains. Reasonable choices are made in consideration of the geological context, the available data, the goals of the study and the available resources for the study. This paper reviews aspects of the decision of stationarity in modern resource assessment with a focus on quantifying uncertainty and using numerical/statistical tools to help the geomodeller. The use of volume/distance functions to quantify uncertainty in domain limits is particularly important and is presented in detail. The use of contact plots provides help with the continuity of grades across domain boundaries. Realisations are not assembled as a branching combinatorial of all possibilities, but as a set of realisations that provide a measure of uncertainty. This paper amounts to a definition of stationarity for modern resource assessment.

Sep 26, 2011

By Mez W

the German ruining industry uses different residues for hackfilling and stowing. A hydraulic pipe conveying technology for these materials has shown advantages for the transportation into the mine. Pipeline conveying systems are able to bridge long distances, they save room within the underground workings and they save staff and costs. Absolutely necessary for successful use of the hydraulic conveying technology is the knowledge about the unusual non-Newtonian flow properties of solid/liquid mixtures like pastes and slurries. The flow resistance of these materials cannot be described easily by one viscosity. Owing to a yield point and a pseudo-plastic behaviour the viscosity changes with the shear rate and with the streaming velocity inside the pipe. Also these flow behaviours depend very sensitively on the liquid/solid-ratio, often so dramatically that gravimetric or volumetric dosing equipment is not in a position to control the production of backfill pastes. Quality changes of the solid raw materials (particle size distribution, density, mineralogical composition) lead to additional difficulties. The first experiences with hydraulic backfilling plants were troubled by plugged pipes. a bad backfill quality caused by unnecessarily high water content and high costs for chemical binders to solidify the backfill. The installation of viscosity controlled paste production leads to success with the hydraulic backfill technology. A development of the l)MT, a bypass pipe-viscometer, has been the key technology. The system is able to analyse on-line, difficult, sticky and abrasive substances and determine their non-Newtonian flow properties. Because the system measures real physical values it is in a position to pre-calculate process relevant parameters like the pressure drop of the backfilling paste in the conveying pipeline. Newer types of the viscometer system also measure the slurry density and calculate its composition. These systems have proven a success controlling backfill paste production plants. They offer it secure prevention of pipeline plugging and a production of high quality backfill pastes with a water content as low as possible.

Jan 1, 1998

By D Hunter

A production information management system (PIMS) was recently implemented at Queensland NickelÆs Yabulu Refinery. Although a strong business case was developed to justify the significant investment in PIMS some potential beneficiaries remained sceptical of the anticipated benefits. The decision to install PIMS was justified as the on-line availability of continuous production data enabled the optimisation of some crucial operations resulting in significant cost savings, improved throughput and final product quality. Furthermore the rapid access to extended historical data has allowed the cause of repeated failure of equipment to be determined. The investment is expected to pay for itself within two years and has already contributed towards significant savings. This paper focuses on some of the diverse production improvements gained and compares the anticipated benefits in the business case to those achieved. A brief outline of the phased implementation of the OSIsoft PI data historian and complementary client software at Yabulu is included, along with details of the lessons learned.

Jan 1, 2003

By N Davies, M J. Coward, P L. Hissey

Mount Isa Mines has four mining operations in Mt Isa with the X41 and Enterprise Mines producing copper and the Lead Mine and George Fisher producing lead zinc and silver. Each stream contributes ore to their respective concentrator at a feed rate and therefore metal ratio set out by the monthly schedule. Mine Planning and scheduling lead to a sales and operating plan (S&OP) and revenue projections which ultimately derives actual revenue /tonne and therefore cost/tonne. At Mount Isa, both the copper and lead/zinc streams produce ore from relatively new operations to be blended with material from more mature operations. It is critical to accurately predict the metal contribution from each operation in order to gauge the relative performance of each mine. Development mining ore contributions are estimated and reconciled, as are stope and bench sources. Performance is continually measured and assessed during extraction and grade adjustments are made to account for underbreak, overbreak and design modifications. It is critical to not understate or overstate the contribution from one operation relative to the other as this can produce a false picture of operational performance. Identifying process improvements are key to improving the business overall. Ore from the lead stream passes 17 measurement stations and two midstream reconciliation points before delivery to the heavy media plant (HMP), where the metallurgical circuit begins. From this point two concentrates are produced. A zinc concentrate for treatment off site and a lead concentrate for further on site processing. The lead smelter takes 100 per cent of Isa concentrates as well as a contribution from BHP-BillitonÆs Cannington Mine.

Jan 1, 2002

By J E. Everett

An essential requirement of product quality control in the mining industry is to be able to reliably predict key quality properties of finished product from the data available before the extraction of the ore. From a production viewpoint, the unit of data collection is generally the input and output data set for each shift of crusher production, but could be any period where mine precrusher data can be reliably matched with product data.This paper demonstrates that linear regression models can be used to predict crush grades from blast grades, even where the crush material is blended from multiple sources or pits for each of which differing regression models might apply. The advantage of using increasingly complex regressions models to match the available data to give increasingly more powerful models is examined. The best model for any application will be a balance between required predictability, available data and the business’ tolerance for complex models.The regression modelling approach has several advantages over the more classic method of run-of-mine (ROM) crusher trials, which will be discussed.The models can use any predictor variable such as grade, geotype and in situ density provided the preextraction data set can be reliably matched with post-crusher data and is significant as a predictor. The models have been used extensively in the generation of the daily crusher plan with the aim of maintaining finished product grade. This approach has also been used associated with exploration drilling and long-term planning.It is acknowledged that there are inherent problems in fitting lump and fines grade to a linear model. However, these problems are minor when they are used for interpolation within the window spanned by the shift blend records used to produce the model. This paper discusses some of the issues limiting linear regression models in this application, and suggests methods enabling consistent models to be formulated

Sep 5, 2011

By D Veneziano, J Van Dyck, R Snow, W Pytel, O Zanoli

Best practice standards generally recommend assessment of the environmental and societal risks posed by tailings storage facilities, but they rarely specify how the analysis should be made. As a case study, this paper describes a safety assessment of the Zelazny Most tailings pond operated by the Hydrotechnical Unit of KGHM Polska Miedz SA in Poland. The complete assessment includes a failure mode effect criticality analysis (FMECA) that considers all conceivable failure modes, but this paper focuses on one key mode, the loss of stability during the planned raising of the crest above its current level. Reliability against static dam failure considers uncertainty on the safety factor in a limit equilibrium approach. The threat of dynamic (earthquake induced) instability is assessed by developing a model of future seismicity as a function of the mining plan and considering uncertainty on the induced permanent displacements along critical sliding surfaces. We quantify the additional level of protection against outflow provided by the beach using an event tree with expert-assigned probabilities. Continuous monitoring of the dam for premonitory signs of instability and a set of loss-prevention plans in the case of a spill provide additional layers of protection; their effectiveness is quantified by estimating the fraction of the population that could be evacuated as a function of warning time. The probabilistic analysis indicates that the safety targets are met and produces practical recommendations on how to maintain and improve the safety level while the dam crest is raised above its current level.CITATION:Van Dyck, J, Pytel, W, Veneziano, D, Snow, R and Zanoli, O, 2015. Safety of a tailings storage facility against dam instability, in Proceedings Tailings and Mine Waste Management for the 21st Century , pp 253–262 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Jul 27, 2015

By H Wen

LongWangGou coalmine auxiliary inclined shaft heading face has a large section, fast heading speed, and large dust production. The problems of high concentration and wide distribution of dust in the roadway air are very prominent. To improve the efficiency of roadway ventilation and dustremoval, long-compression and short-suction (after this referred to as LCSS) are adopted. The dust concentration and airflow field distribution in the roadway are numerically simulated by FLUENT. The change of the position between the pressure duct and the exhaust duct and the different pressure-extraction ratios affect the distribution of airflow velocity, the formation of vortex zone, and the distribution of dust concentration in the roadway. When Ly = 25 m, Lc = 5.0 m, W = 1.2, a wind curtain is formed at 5 m in front of the heading face to prevent dust from spreading deeper into the roadway.

Oct 12, 2022

Since the early 1990s it has been evident that electricity demand has outstripped investment in generation and transmission. At the same time, concerns were being expressed about the risks associated with the run down of the Maui gas field. It is axiomatic that the necessary investments firstly require a stable and supportive regulatory and fiscal environment, and secondly that such investments have long lead times. A unwillingness by successive Governments to fund the necessary infrastructure, however, lead to the profound restructuring of the electricity with the ultimate objective of transferring the investment risk from Government to the private sector. Successive power crises since 1992 result from a failure to put in place a regulatory framework necessary to provide the environment to allow the private sector to make the needed investments within the necessary timeframe. The electricity sector has now stabilised to the extent that the regulatory risks associated with generation investment have been reduced. This does not mean that risk has been eliminated, however, and both economic and environmental issues will continue to constrain investment in some power proposals particularly those in the South Island and especially any South Island coal or gas thermal proposal. This is firstly because of several very simple economic issues:  1. while the South Island may have significant coal resources it lacks a significant electricity market; 2. except for a six to ten week period in a dry year, low cost hydro will price any large-scale South Island thermal out of the market which will restrict any large scale thermal station to a dry year and/or peaking load; 3. except for a six to ten week period in a dry year, export capacity on the HVDC link to the North Island is constrained to around 964 MWe (and 482 MWe from North to South in dry years); and 4. construction of any new thermal plant in the northern half of the North Island will price any South Island thermal station out of the market. There are also significant environmental issues associated with the establishment of a large scale thermally based power station in the South Island. Because the South Island lacks any large scale economically proven gas reserves, the environmental constraints considered are those associated with coal production. A coal mine dedicated to power generation will have significant on-site impacts associated, in particular those associated with disposal of opencast overburden and visual impacts. These will, however, be partially offset by the relative remoteness of such sites to population centres. What cannot be offset, however, is the off-site impacts, in particular, CO2, SOx and NOx emissions, water disposal and acid mine drainage.

Jan 1, 2003

By D K. Weatherley, A N. Mainza, I Govender, L M. Tavares, M S. Powell

It is understood that rocks do not instantaneously transport along a grinding mill and are consequently not perfectly mixed, but this is conveniently overlooked in mill modelling. A lack of understanding of the transport phenomenon and the difficulty of modelling realistic transport in a complex tumbling system, contribute to this. Visual observation of segregation in AG and SAG mills and measured lag time of mill response to changes in feed provide evidence of the influence of axial transport within mills. Incorporating this lag in dynamic models should underpin dynamic simulation of milling circuits for design and for model-based control. Limitations and errors in predicting the influence of increased power via mill length versus diameter, inability to predict pebble production, and a lack of the differential in residence time across rock sizes and competences, all contribute to inadequate design prediction that can be especially acute for large Autogenous mills. Examples of dramatically under – and overloaded pebble regrind circuits and the inability to predict pebble load in AG milling discouraging the application of this energy-efficient option, point to the need for predictive modelling of breakage by competence linked to transport through and out of mills. Incorporating progress in granular flow and mixing, coupled to computational advances enable models to be based on complex breakage and transport models at high spatial resolution. This allows mechanistic models to simulate thousands of times faster than real time, opening the door to more accurate mill design, including dynamics in process simulation and to live simulation for process control – all with advanced mechanistic models. The value of grasping this opportunity is not easily refuted.

Aug 24, 2022

By P J. Hawke, M F. Flis

The Hamersley Basin of Western Australia holds significant tonnages of detrital iron deposits. Formed by a cyclic weathering and eroding sequence of nominally barren banded iron formation, they may form rich proximal accumulations of iron cemented gravels distal from any obvious hard-rock iron mineralisation. While generally small, they nonetheless represent excellent economic targets as they are near surface and easy to mine. The relatively high density of iron detritals, when compared with their sedimentary hosts, makes the gravity method the primary method of exploring for these blind deposits. Identification of possible trap sites is used as a precursor to gravity surveying, with magnetics used as an ancillary method prior to drilling. Application of the gravity method in this terrain is not straightforward. Detrital trap sites are best developed next to the spectacular cliffs that form the range fronts to the Hamersley Ranges. Precise terrain corrections and use of first vertical derivatives are routinely applied to gravity data in an effort to map every possible gravity high that may reflect the presence of a detrital accumulation. Extensions to Hamersley Iron Pty LimitedÆs BS2D deposit were discovered by the routine application of gravity. This methodology has been used to create a model by which exploration for detrital deposits continue to the present day.

Jan 1, 2002

By M Zoorabadi, B Hebblewhite, S Saydam, W Timms

Based on the cubic law, the hydraulic aperture of rock discontinuities has the most significant effect on the hydraulic conductivity of jointed rocks. The hydraulic aperture of rock discontinuities is a stress-dependant parameter and reduces with increasing applied stress. Furthermore, the rock stress increases with increasing depth. Therefore, hydraulic aperture is a depth-dependant parameter. In this paper, new empirical equations are introduced to calculate the hydraulic aperture of rock discontinuities at different depths. These equations were developed on the basis of back analysis of compiled packer tests from different actual cases. They provide the capability to calculate range of variation and average hydraulic aperture when the orientation of discontinuity and field stress tensor are given. Furthermore, it is possible to calculate the hydraulic aperture of rock discontinuities when its value is given for another depth. A hydraulic aperture of 7.7?µm was calculated by this method, comparing well with a reported site value of 8.3?µm for rock discontinuities at the Stripa mine project at a depth of 385?m. These new equations can be used by analytical or numerical methods to estimate the hydraulic conductivity of jointed rocks.CITATION:Zoorabadi, M, Saydam, S, Timms, W and Hebblewhite, B, 2014. New empirical equations to estimate the hydraulic aperture of rock discontinuities at different depths, in Proceedings AusRock 2014: Third Australasian Ground Control in Mining Conference , pp 375–380 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Nov 5, 2014

By G W. Wilson

Net solar radiation is one of the key factors affecting evaporation from soil covers on waste rock dumps. While net solar radiation is normally taken as a uniform parameter over a site for modelling purposes, there can be significant variation in the radiation actually received on various portions of the dump. These variations are a function of the difference in the cover slope and direction of exposure. A predictive model has been developed to calculate the actual radiation received on the dump surfaces, based on the net radiation data measured at a weather station on a horizontal surface. The model has been designed for application at any mine site in the world, and can be used in conjunction with evaporative flux models to help estimate the effect of three-dimensional variations in the cover surface exposure. The waste rock dump at Equity Silver site (located in British Columbia, Canada) has a barrier-type soil cover system to limit oxygen entry and water infiltration. This cover has been monitored for the past nine years. The new net solar radiation model has been applied to this site and the predicted net radiation was compared to the values of actual net radiation measured on the site. The predicted variations in net radiation over the three-dimensional surface of the waste rock dump have been applied to calculate spatial variations in potential and actual evaporation (or evapotranspiration) from the site. Preliminary analyses indicate actual evaporation from a site at this latitude could vary by a factor of two for exposed versus shaded slopes. The implications of these variations for the site are considered.

Jan 1, 2003

By P J. Terbrugge

Many ore deposits have considerable vertical extent which is not always known at the time of commencement of mining. In many of these cases, initial mining is by opencast methods, and as ore reserves are proved to greater depths, the pits are often planned to go deeper than had been originally envisaged. In such cases, surface plants and critical underground facilities such as conveyer tunnels, access ramps, ore passes, hoisting shafts and ventilation shafts are often located much closer to the pit rim and the orebody than is desirable for the deeper mining situation. This usually raises the question of the stability of such underground excavations, which may be critical for the viable longer-term mass underground operation of the mine. In this paper the key considerations in the evaluation of the stability of the overall mining environment for ædeepened minesÆ are defined. Aspects which are dealt with are: open pit stability, shaft stability, dilution, mudrushes and airblasts, and selection of mining method. The necessary stability evaluations are put into the appropriate context by considering two case histories. Critical parameters and considerations which need to be addressed during the evaluation of excavation interaction are also discussed.

Jan 1, 2000

By M I. Pownceby, N A. S Webster, A J. Studer, I C. Madsen, J A. Kimpton, M J. Fisher-White, L Lu

In situ synchrotron X-ray and in situ neutron diffraction-based experimentation have been implemented to characterise the formation of the complex calcium ferrite iron ore sinter bonding phases silico-ferrite of calcium and aluminium (SFCA) and SFCA-I. Experiments were carried out using fine-grained (<20 ¦m) sinter mixtures containing ~77 wt per cent Fe2O3, 14 wt per cent CaO, 3.5 wt per cent SiO2 and 5 wt per cent Al2O3, with diffraction data collected during heating until 1300 - 1350¦C, which was enough to ensure complete melting of the SFCA phases. In the case of the in situ synchrotron experimentation, results showed that altering the nature of the starting sinter mixture (ie substitution of goethite for hematite; substitution of an amorphous Al2O3-based mineral for gibbsite) did not have a significant effect on the thermal stability range of the SFCA phases. It did, however, have a profound effect on the formation and rate of consumption of the calcium-rich ferrite phases C2F and CF, and on the formation mechanisms of SFCA and SFCA-I. The in situ neutron diffraction experimentation is the first described in the context of iron ore sintering, and the observation of both SFCA-I and SFCA formation in the in situ diffraction data represents the first step in achieving the goal of charactering iron ore sinter phase formation in large volumes of industrial sinter starting materials with a wide range of particle sizes (up to 6.3 mm).CITATION:Webster, N A S, Pownceby, M I, Madsen, I C, Studer, A J, Lu, L, Manuel, J R, Kimpton, J A and Fisher-White, M J, 2013. In situ x-ray and neutron diffraction studies of silico-ferrite of calcium and aluminium iron ore sinter phase formation, in Proceedings Iron Ore 2013 , pp 487-494 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 12, 2013

Project-finance risk assessment is mature and approached systematically. Mining industry due diligence remains fragmented with specialists addressing different aspects of the project such as reserves/JORC, environmental impact, commodity market, even financial analysis/ modelling. Each professional has a natural tendency to believe that the risk factors they have examined are the key ones necessary to be reviewed, as if their review will stand in isolation from others in the way a project materialises. By first laying out the project risk categories and risk matrices, this paper examines the due diligence process with a keen eye to the risk assessment and overlap/interaction. Some cross fertilisation of due diligence techniques from other industry sectors is placed into the mining context. Six case studies are then presented with regards to due diligence features in each projectÆs feasibility review or financial information/ offering memorandum. A power case first shows the differences between the due diligence expectations of banks versus bond financings. Then some Australian and Indonesian mining development cases are analysed to see the due diligence learning points. Some of these were project/due diligence failures. The conclusion focuses on the value of a systematic risk assessment for project evaluation, all from the project-finance perspective.

Jan 1, 2007

By D Pearce

The 15-year long debate between the industry, Government and Aboriginal people about the rights and wrongs of Land Rights, the right of Aboriginal people to control development on their land and the relationship between the parties, has been long and often bitter. It has so far been a war of words and often the truth has been the first casualty. This paper looks at the debate and presents the perspective of Aboriginal people, as well as looking at how the industry and Aboriginal people can work together for everyone's benefit. It addresses such matters as: ò why Aboriginal people believe they have the right to control development on their land; ò the perception that Aboriginal people are anti-development; ò the perception that Aboriginal land is 'locked-up'; and ò how companies can negotiate with Aboriginal people.

Jan 1, 1994

By J J. Eksteen

"In the platinum group metals (PGM) industry, smelters are often perceived as ‘black boxes’ that have to process concentrates of variable mineralogy and originating from multiple mines operating along strike or downdip from each other in the same geographical region, as is the case in the Western Bushveld of South Africa. However, PGM smelter complexes are much more sensitive to concentrate blend variation than most metallurgists or geologists would anticipate.This paper will review the impact of concentrate blends deriving from a number of concentrators processing ores from the Bushveld Igneous Complex (BIC) on smelting and converting behaviour and the constraints they impose upon smelter operation. Significant variation in behaviour is not only noted between the main mineralisation types, namely the Merensky and Upper Group 2 (UG2) reefs, but also within a given reef type. Due to the way PGMs are mineralised and finely disseminated, concentrators have to grind concentrates ultra-fine to liberate the PGMs and provide a sufficiently high-grade concentrate for smelting.The chemistry and mineralogy of the concentrates vary with the type of mine (open pit versus underground), the reef type and varies along strike; this variation impacts numerous key smelter parameters such as the base metals content, the matte fall, the converter aisle capacity, sulfur dioxide gas production, furnace matte and slag temperatures, and the accumulation of high smelting chromite-spinel minerals that has led to significant smelter upset conditions such as catastrophic furnace hearth and sidewall failures.Due to the production intensity of smelting, ie its high rate of processing per unit volume, smelters have been found to be far more sensitive to rapid variation in concentrate properties, highlighting the need for adequate blend control.CITATION:Eksteen, J J, 2016. Impact of Mineralogical Variation on the Smelting of Cu-Ni-PGM Concentrates, in Proceedings The Third AusIMM International Geometallurgy Conference (GeoMet) 2016, pp 321–326 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Jun 15, 2016

By Caress ME, Knight MD

The eruption 75,000 y ago of the Younger Toba Tuff ignimbrite and associated ash-fall ranks as the greatest known volcanic event of the past million years. The ignimbrite is known to cover 21,000 km2 (it may extend farther under alluvium of the coastal plain) and forms a low-angle depositional shield centered on the Lake Toba caldera. The ignimbrite is largely nonwelded and shows systematic outward decreases in lithic and crystal contents and clast sizes. Coarse lithic lag breccias that occur at places 85 km apart point to two or more vents on a fissure having this length. The extensive associated ash-fall traceable as far as Malaysia and India is evidently a co-ignimbrite ash but could have a phreatomag-matic component. The Lake Toba caldera measures95 x 35 km and is >1 kn deep. It is not a single caldera but comprises several contiguous depressions and grabens diversified by two re-surgent domes. Earlier parts of the caldera are related to at least two earlier ignimbrites including the Older Toba Tuff dated by Chesner at 0.5 Ma. The caldera is elongated parallel to, and is only 25 km from, the great Sumatra strike-slip fault, and it can be considered almost certain that the caldera is in part a kind of pull-apart structure. The ignimbrites are very crystal-rich. There is no reason to think basaltic magma was involved in their emplacement. Together with ash-falls they have a total volume of several thousand cubic kilometers. They can be categorized as originating from partially crystallized granitic magmas that erupted because tectonic pull-apart occurred several times across the pluton. Small rhyolitic lava domes on fault-lines comprise the youngest volcanic products.

Jan 1, 1990