Search Documents

By Brian G. D. Smart

A scheme is presented whereby the complexities of design and operation of a hydraulically powered coalmine support are modelled using the finite element method. This enables predictions to be made of the interaction of the support with its environment in a variety of operating conditions. The model is a two dimensional representation. The modelling scheme was applied in two case studies. In the first, the performance of a four leg support was uprated by increasing the front leg force. The finite element model enabled changes in overall support performance to be verified. In the second case study, the change in lateral force developed by a two leg support was assessed with changes in support operating height and floor foundation stiffness. The modelling scheme is capable of further development with the possibility of extending the model into three dimensions and the addition of control elements that would enable an actual mining sequence to be represented.

Jan 1, 1992

By L. Mahony

This paper describes the design, construction and commissioning testing program for a laboratory facility at the School of Mining Engineering, UNSW, Australia, aimed at evaluating reinforcement tendon performance when subjected to shear loading. The project has been financially supported through the coal industry ACARP research initiative. The single failure plane design adopted in the test rig has been successful in allowing shear loading to be directly applied to fully installed (and where relevant, pre-tensioned) rock bolts. Initial results are extremely encouraging, and illustrate a significant axial component of load development as shearing takes place - particular in lower strength, softer material types. The paper also outlines future plans for use of this facility.

Jan 1, 2005

By Erik Westman, Benjamin Mirabile

"Quantitative evaluation of secondary support performance in longwall tailgates in U.S. coal mines was introduced with the U.S. Bureau of Mines (USBM) Wood Crib Performance Model. This model evaluated crib performance based on load capacity, stability, and convergence characteristics. The introduction of a large variety of alternative secondary supports in the 1990s prompted a full-scale load-displacement testing program conducted by the USBM and, subsequently, NIOSH. The application of the ground response curve concept to secondary supports provided a more complete estimation of secondary support loading in longwall abutment zones. Secondary support loading from floor heave, horizontal stress effects, and pillar yield are captured by the ground response curve. A point on the ground response curve can be measured by instrumenting a secondary support with load and convergence measurement devices. A segment of the ground response curve can be measured by instrumenting two secondary support types of varying stiffness. Numerical models have been used to construct complete ground response curves for longwall tailgates based on measured secondary support performance. The ground response curve concept can be used as a design methodology for secondary supports in the #2 entry, as well as the longwall tailgate. Extending published evaluation and design methodologies for secondary support in longwall tailgates, a ground response curve unique to the loading conditions and required functionality of the #2 entry can be developed. Numerical models can be used to construct full ground response curves based on measured data in the #2 entry.INTRODUCTIONA longwall ventilation system is designed to continuously dilute and move methane-air mixtures and other contaminants from the active face, to bleeder systems, and out of the mine. The #2 entry of the tailgate gateroad serves to transport contaminated air from the longwall face to the bleeder system. The #2 entry is required to function as an air course between two longwall gobs and is subject to significant ground movement and loading. Current ground support design and evaluation practices are based on empirical evidence and little to no quantitative data on support performance exists for these entries. Substantial volumes of research have been conducted on ground movement and support performance in longwall tailgate entries. Following the successful ground support evaluation and design practices established for longwall tailgate entries, a practical and efficient design methodology can be adapted for supporting the #2 entry of a three-entry gateroad system."

Jan 1, 2018

By Mingju Liu

In this research, we measured fractal characteristics of over ten typical outburst areas on different scales and related them to their outburst proneness, in attempts to provide new methods for outburst prediction. It was shown that both faults and folds follow the power law on coal basin scale, mine scale or mining district scale, whichever it is. Fractal dimensions are appropriate indices to describe geological structures quantitatively. Preliminary relationships have been established between outburst proneness and fractal dimension. Comparisons have been made among fractal characteristics and degrees of outburst proneness on all three scales. The results indicate that multi-scale assessment of outburst proneness can be made based on fractal measurements independently or supplemented by other assessment methods.

Jan 1, 2000

By Quanxi Wang

This paper presents a new stability mapping system which tightly integrates structural and geologic mapping with geo- mechanical stress analysis for use in support design and mine plan- ning. The structural and geologic input can include such local characteristics as: rock strength, layer thickness, faults, stream val- ley effects, etc. The geo-mechanical stress inputs can include overburden stress, multiple-seam stresses, abutment stresses, etc. The mapping system utilizes the popular AutoCAD/SurvCADD platform as a foundation for inputting the geologic characteristics, and then integrates the boundary-element program LaModel for determining the geo-mechanical stress influences. Additional func- tionality in the system allows the user to transform each individual geologic or stress factor into an individual rating and then the user can apply various weighting functions to the individual factors to generate an overall stability index map for the mine property. A sample case study is included in the paper to demonstrate the basic concepts, procedures and utility of the stability mapping system.

Jan 1, 2005

By William B. Beerbower

When full length resin grouted bolts were introduced into U. S. coal mines some 11 years ago, it was assumed that safety would be the primary benefit. There is ample testimony to substantiate the safety factor as an important reason for using resin roof bolts. More recently, it has been determined that many coal operators not only concur with the safety consideration but also benefit economically from the safer working conditions resulting from a reduction in roof falls. IMPROVED SAFETY Some typical roof fall reduction experiences with resin bolting follow: ? Western Kentucky Mine - Resin bolting reduced roof falls by 80%. ? Southern West Virginia Mine - Experienced 141 roof falls between 1971 and 1972. In 1973 with resin bolting, roof falls were reduced 85%. ? Western Pennsylvania Mine - In 1973 had 59 roof falls with mechanical bolts. During the same period of 1974 with resin, only 19 falls were experienced - a reduction of 66%. The overall experience with resin bolted roof support has been an average 50% reduction in roof falls. Long term, this can translate into improved worker morale. As a matter of fact, a general mine foreman at the same Southern West Virginia mine mentioned above commented that improved miner morale was noticeable following the introduction of resin bolting. HOW GREATER SAFETY CAN BE COMBINED WITH INCREASED PRODUCTIVITY IN COAL MINING OPERATIONS Experience has also shown that a reduction of roof falls with resin bolting has resulted in increased production for many coal operators. Some typical examples with mines using resin bolting follow: ? Western Kentucky Mine - Resin bolting contributed to increasing production from 7000 to 11,000 tons per day, due to a 80% reduction in roof falls. This was accomplished in part by a reduction of 50% in timbering requirements and the reassignment of three roof fall cleanup crews to coal production. ? Pennsylvania Mine - Following resin bolting introduction the section showed a 50% increase in production. ? Central West Virginia Mine - Continuous miner production averaged 6 cuts per unit shift prior to resin bolting. Reduction of cleanup time, resulting from improved roof support with resin bolting, increased production to an average of 9 cuts per unit shift. PROPERTIES OF RESIN BOLTING SYSTEMS Several of the important properties of resin bolting systems which contribute to its improved performance are listed below: ? Full length anchored resin bolts and many resin anchored, point anchored bolts provide effective resistance to changing rock strata stresses throughout their length, thereby improving stabilization of rock strata over long periods of time. ? Individual layers of rock are anchored into a single high strength beam providing excellent resistance to vertical and lateral strata movement, the most common cause of roof failure. ? Moisture cannot corrode the fully grouted bolt. ? Has twice the pull strength of mechanical bolts (See Figure 1)

Jan 1, 1982

By William Kendall

INTRODUCTION Stillwater Mining Company, with its two underground, hard rock mines in South Central Montana, is the largest primary producer of palladium and platinum in the Western Hemisphere. Platinum Group Metals (PGMs) are mined from a specific zone within the Stillwater Complex, a layered mafic/ultramafic igneous intrusion, named as the J-M Reef. Over much of the J-M Reef, the economically minable, mineralized horizon is contained at typical horizontal widths of 6.5 ft to 8 ft wide. This feature of narrow mining widths has made it challenging to safely and effectively mine these areas of the J-M Reef using mechanized cut and fill and sublevel sill development mining methods. Mechanization of most or all of the steps of the mining process is important to Stillwater in the continuous pursuit of high standards in safety and health. Mechanized bolting within the J-M Reef ore production headings is an important piece of reaching these objectives. In general, metal/non-metal mining requires varied mining techniques, as compared to more traditional bedded mineral/resource deposits. Narrow vein mining is a selective mining technique that generates a minimal amount of waste material while obtaining the more valuable ore deposits. Drilling and bolting in narrow vein mining has traditionally been performed with air-powered, hand-held ?jackleg? equipment because mechanized equipment is, typically, too large for the selective narrow vein techniques. In addition, traditional hand-held pneumatic equipment provides significant flexibility for small-opening, narrow-vein developments, reducing waste rock handling and improving ore recovery management. The hand-held equipment is also relatively inexpensive to purchase. Due to its simplicity, it is also relatively easy to train operators to use hand-held equipment. Hand-held equipment has been routinely used to safely install ground support and to develop production faces in hard rock mines for many years. Although it is effective, the hand-held equipment requires operator exposure to certain inherent hazardous conditions in the active face/development mine areas. These include falling rock, tripping hazards, bending/twisting/lifting hazards, high noise levels, and chemical fumes from the hammer exhaust. Based on health and safety statistics from the Stillwater mines, operation of the rotary-percussive hand-held equipment adversely affects the long-term health of the operator?continual jackleg operation frequently results in long term damage to hearing, shoulders, arm and wrist joints, and back injuries. This reduces the quality of life for the operator and also results in higher personnel turnover rates for the mining company.

Jan 1, 2014

By Simon H. Prassetyo

Violent failures of coal pillars, known in practice as coal mine bumps, have long been a subject of investigation. Many field investigations have considered geological conditions that create high stress in the pillar to be the main causative factor leading to bumps. In recent years, constraint at interfaces has been observed to be associated with the occurrence of coal bursts. This paper shows how interface friction and width-to-height (W/H) ratio affect the degree of violence of coal specimen failure. A series of UCS tests on 178 coal samples from bump-prone (Utah, called UT coal) and non-bump-prone (West Virginia, called WV coal) mines were conducted. Three violent failure parameters? peak sound pressure level, core zone failure, and ultimate stress? were used to assess the violence of failure. The degree of violence were investigated at varying interface frictions (high: µ = 0.40, medium: µ = 0.22, and low: µ = 0.13) and W/H ratios (W/H = 2, 3, 4, 5, 6, 7, 8, 10, and 12). A simulation series of modeled coal specimen under uniaxial loading was also conducted using the finite element method to observe the equivalent plastic strain distribution experienced by the modeled coal specimens. It was found that the violence of coal specimen failure depends on W/H ratio and interface friction. Interface friction has more significant influence on the violence of failure at high W/H ratios than it does at low W/H ratios.

Jan 1, 2011

By Peter Altounyan

Goedehoop Colliery in the Witbank Coalfield of South Africa produces some 8 million tonnes of coal from fully mechanised room and pillar mining methods, using resin anchored bolts as roof support. Investigation of the quality of the installed bolts resulted in the finding that it was possible to pull out fully encapsulated roofbolts. This led to a drive to improve the quality of installation by making it easier to install a roofbolt correctly and more difficult to install incorrectly, whilst at the same time achieving high bond strength. The method of installation which was investigated, developed and implemented is known as 'spin to stall.' The roofbolt is pushed through the resin capsule to the top of the hole and then spun to mix the resin, until the increasing resistance to rotation exceeds the breakout limit on the bolt torque nut. The nut then runs up the bolt thread and is tightened to the roof as the drill stalls. This procedure is simple and avoids human factors such as non adherence to 'spin' and 'wait' times which contribute to poor installation with conventional methods. However spin to stall installation has been avoided in the past because of concerns over damage to the resin bond. This paper reports on the efforts made to investigate and understand the factors involved, and the progress made in overcoming the problem of resin damage so that the advantages of this installation method can be realised.

Jan 1, 2003

By Bruce K. Hebblewhite

BHP Billiton Illawarra Coal operates Dendrobium Mine in an area 10-20km west-northwest of Wollongong in the Southern Coalfield of New South Wales, Australia. The mine recently completed longwall mining in Area 3A adjacent to an overhanging rock feature known as Sandy Creek Waterfall, which had been identified as a significant, and potentially sensitive natural surface feature of the environment. This waterfall lies within the Sydney drinking water catchment area, in a surface terrain of variable topography where non-conventional subsidence behaviour such as valley closure, and far-field horizontal movements are known to occur, over and above conventional subsidence responses to mining. As a part of the conditions of approval to mine, Illawarra Coal undertook to protect the waterfall and the section of Sandy Creek immediately upstream of the waterfall from the effects of longwall mining using an innovative management process and an array of very high resolution monitoring systems. This paper outlines the type of non-conventional subsidence effects that can occur; the range of different monitoring regimes adopted in response to such effects in order to develop an understanding of the valley closure phenomenon in proximity to the waterfall and provide suitable early warning data. The paper also describes the management processes that were adopted, incorporating the use of these monitoring regimes in order to successfully protect the waterfall from the valley closure effects of mining four adjacent longwall panels in close proximity to the waterfall while continuing to maximise recovery of the coal resource in the area. The management structure adopted involved a technical committee, a steering committee, and an external independent reviewer. The technical committee was responsible for design of the monitoring systems, interpretation of the results, and making timely recommendations to the company steering committee who were then able to make informed decisions on when to cease mining each adjacent panel. The steering committee comprised Illawarra Coal management and technical personnel. Although the steering committee took advice from the technical committee, all decisions relating to mining were made by the steering committee. An external reviewer was engaged by the steering committee at the end of each longwall panel to review the results, interpretation and management decisions. This management structure was effective in successfully protecting the very sensitive structure of Sandy Creek Waterfall from potential impacts of nearby mining in the midst of active ongoing natural erosion processes.

Jan 1, 2014

By Nicole E. Evanek, Anthony T. Iannacchione

"Over the past 36 years, a great deal of information has been presented and discussed at the International Conference on Ground Control in Mining (ICGCM). The topic of surface subsidence engineering has resonated throughout each of the 36 conferences, producing 155 related papers. These papers have come from many countries. The major focus has been underground coal, especially longwall mining, but hard rock mining has also received attention. This impressive collection of reports presents a unique opportunity to reflect on the progress made towards understanding subsidence basin characteristics and developing methods to predict and mitigate its impact.In the initial meetings, the focus was on measuring subsidence and developing prediction methods. As understanding increased, predictive methods were developed and tested, especially in North America, China, Australia, Europe, and India. Efforts were made to understand how structures, aquifers, and surface waters were impacted by the formation of subsidence basins. As time progressed, the number of papers dropped slightly and became more focused on site-specific conditions. Most recently, advanced monitoring techniques have revealed new information, spurring a resurgence in research activity. Another recent focus has been on anomalous horizontal deformations occurring in sites with significant relief. Several reports from Australia and one from the US have brought this issue to light. In closing, all of these surface subsidence engineering papers have chronicled the maturing of this most important topic and, in the process, highlighted the importance of the ICGCM in the evolution of our current state of practice.INTRODUCTIONThe aim of this paper is to examine the role of ICGCM in 1) shaping our understanding of the mining induced subsidence phenomena, 2) developing tools to help predict its impact on the surface, and 3) evaluating controls to mitigate these impacts.This analysis involved the consideration of 155 papers available within the ICGCM conference proceedings and listed under the subsidence topic. Three of the papers where only available as abstracts and were not included in this analysis. Four other reports were focused on topics other than underground coal mining-induced subsidence issues and were also not included. It should also be noted that papers dealing with overburden behavior were not included in this analysis. As a result, 148 papers were reviewed and analyzed (Figure 1)."

Jan 1, 2018

By John Kucish

Roof bolting in the Pittsburgh Coalbed takes many forms today: coupled partially grouted bolts, fully grouted rebar, passive cable bolts, tensioned cable bolts, and torque tension rebar. The Federal No. 2 Mine of Peabody Energy has recently made a significant advance in longwall headgate primary roof support, with the change fiom an 8 ft long coupled partially grouted bolt to a 6 ft long torque tension system. The coupled partially grouted bolt consists of a 4 ft long 0.804 in diameter rolled deformation bar thread coupled to a smooth 718 in headed bar. The coupled partially grouted bolt is anchored into a 1 318 in borehole with 4 ft of high speed resin. The 718 in diameter smooth bar is tensioned by rotation into the treaded coupler, upon setting of the resin, compressing the immediate roof. The change was initially driven by occasional failure of the coupled partially grouted bolt at the threaded coupler, due to lateral movement caused by horizontal stress. A fully grouted % in diameter torque tension bolt, in a 1 in borehole, was proposed as a cost effective primary bolting for the high horizontal stress conditions. The upper 2 112 ft of the rebar is anchored in fast set, medium viscosity resin and the bottom 3 K ft is anchored in slow set, special low insertion force resin. The lower end of the rebar is threaded with a nut that permits the counter clockwise spinning of the bolt to mix the resin, at maximum rotational torque. Upon setting of the high speed top resin, the nut is spun clockwise to tighten the nut against a metal dome plate. The plate is forced against the roof to compress the immediate roof, A complete headgate entry has been successfully supported with the torque tension system. The mining of the adjacent longwall panel was completed in September 2004. The torque tension system provided a 25% savings over the previously employed coupled partially grouted bolt system.

Jan 1, 2005

By Lorraine Kent

Excessive floor and side deformation are major problems in deep coal mine entries in Europe, and in particular in the UK where rockbolt support is used at depths of working of around 3000 ft. This has major economic and safety implications for the industry resulting in the potential for sidewall instability, disruption to mine infrastructure, (belt conveyors, floor mounted transport systems and powered supports), and the requirement to safely manage excavation of large quantities of floor material. European research aimed at improving deformation control is described. The outcomes include innovations in support systems, materials and installation practices, as well as better methods for predicting and monitoring deformation levels. As coal mines get deeper, some entry deformation is inevitable, and methods for mechanised removal of floor lift and repair of sidewalls have been developed to allow production to continue as mining clearances are re-established.

Jan 1, 2014

By Michael Jenkins

This paper reviews the history of rock burst research in the Coeur d'Alene Mining District of northern Idaho over the last 30 years. Rock bursting is a problem because of a combination of geological factors in the region--complex structure, high horizontal stresses, and strong, brittle rock that readily stores strain energy. The Bureau of Mines has cooperated with industry to study rock burst phenomenon for many years. Early work focused on methods to mitigate bursting problems within the context of existing mining methods. In the 1970's, researchers began investigating methods for altering or controlling conditions that lead to rock bursts. These control measures become part of the stope development that precedes mining or are changes in the mining process itself. Although, researchers continue to develop microseismic techniques, the major research effort focuses on developing a new mining method to reduce the number and severity of rock bursts in order to improve safety and increase productivity in the Coeur d'Alene Mining District.

Jan 1, 1987

By J. R. M. Hill

The U.S. Bureau of Mines is involved in research to improve geological sensing (or geo-sensing) techniques. This paper presents results from field trials of two research projects concerned with new technologies for geosensing. The goal of the first project was to develop a "smart" roof bolter drill so that it could probe into roof rock to evaluate rock conditions. The goal of the second project was to modify an existing mine-wide monitoring system in a continuous mining operation so that the system could interface with geotechnical instruments to provide near real-time geosensing. Integration of this mine-wide monitoring system with the smart roof drill to collect, reduce, and analyze geotechnical data and the potential applications of the resulting technology in mining, are described.

Jan 1, 1992

By Li-min Liu

The Semi-Analytical method can be used to solve a variety of problems in geomechanics and civil engineering. When it comes to layered rock strata and rock masses, it is superior to all other analytical methods. This paper presents the application of a 3-D semi- analytical method to 3-D surface subsidence and stress analysis caused by mining. In this semi-analytical method the layered elements and triangular elements ate used alternatively. Their displacement pattern, strain mahix, elastic matrix, stiffiness matrix, load matrix and stress matrix are presented. A numerical simulation program based on the semi-analytical method bas been developed using Fortran 90. This program is suitable for subsidence prediction and stress analysis of mom and pillar mining, strip longwall and mechanized longwall mining methods. It is applicable to all geological conditions and mining methods. Comparing with the finite method, this method reduces significantly the number of simultaneous equations to be solved and the amount of data preparation. It not only reduces the computer memory, but also runs faster and cost less. Key Words Semi-analytical method, surface subsidence, stress analysis, layered element, triangular prism element

Jan 1, 2002

By Ross Seedsman

There are many models for how roof bolts behave, one of which is prevention of shear along bedding. The basis of this model is that a roof beam will not delaminate if bedding parallel shear can be prevented. A key point in such a model is that bedding discontinuities are present naturally in the roof and that mining induced stresses do not generate new surfaces. The presence of these discontinuities means that the rock mass must be modelled as a transversely isotropic material if continuum codes are used. A limit-equilibrium type design method for specifying roof bolts is being developed using voussoir beam and excess shear force concepts. The method incorporates much of what operators know intrinsically about roof reinforcement, and allows mining engineers to focus on the key uncertainty - the frequency and location of bedding surfaces in the immediate roof. Areas where more research is required are highlighted.

Jan 1, 2000

By Amirhossein Molavi Tabrizi

In this paper we introduce a realistic model for mine openings that explicitly considers the effects of local irregular topography. Also, we consider the effect of entry orientation and rock anisotropy on the induced stress fields around mine openings under the given in-situ stress field. A triangular mesh formulation is derived for representing complicated irregular surface topography and underground openings. Each element contains three internal collocation nodes. Parts of the surface which are far from the mine openings are discretized by infinite elements. The rectangular infinite elements with three internal nodes are introduced specifically for this problem so they are compatible with the triangular elements. With the corresponding boundary element method (BEM) formulation, we then investigate the combined effect of the topography and entry orientation on the stress fields induced by openings in anisotropic rocks under in-situ stresses. A single rectangular opening in the anisotropic ground with topography and in-situ stresses is simulated using data from the Carroll Hollow mine in Ohio. The results of these simulations show that topography, rock anisotropy, orientation of the regional stress field and mine opening geometry all have appreciable effects on the stress fields around underground openings, yet these factors can all be accounted for in an efficient code with no need to discretize the model domain. Future applications using this new BEM formulation will provide useful guidance on selecting the optimal opening orientation in the mine in order to minimize undesirable stress concentrations around underground mine workings.

Jan 1, 2014

By Sheng Zhang

The compressive strength of coal in Yanlong mine area of China is less than 3 MPa. Basically, it is powdery. The roof and floor rocks are mudstone. Therefore this coal seam is a typical ?three soft? coal seam. Roof bolting cannot be used due to the very soft strata. When the conventional U-steel support was used, roadways deformation exceeded 30% of its original dimensions in some sections one half month after development. The required repair work was huge, seriously affecting the daily mine production. This paper describes how to support the soft coal roadways. The deformation characteristics of soft coal roadway were investigated. The conventional U-steel support is not subjected to uniform load, so the whole supporting structure must be strengthened. In addition, drilling holes in the coal seam can transfer high stress from near the ribs to the deeper parts of the surrounding strata, which is beneficial to the stability of the roadway. This new support method is to combine the strengthened U-steel support with pressure-relieving drill holes, providing an economic and efficient way to support the very soft coal roadways.

Jan 1, 2013

By Bob Coutts, Matt Martin, Dan Lynch, Dan Payne

"The Broadmeadow punch longwall coal mine in Central Queensland Australia has experienced significant highwall movement associated with the effect of longwall subsidence when the longwalls approach their final position close to the open-cut highwall. In response to this movement, Broadmeadow employed two types of broadscale highwall monitoring (radar and laser scanners) to provide full coverage measurement throughout three consecutive longwalls approaching the highwall. This was to attain a better understanding of the mechanism causing the movement and potentially enable prediction of instability. Results from the monitoring found the highwall is displaced to magnitudes unlike those typically measured in open-cut mining and in direct contrast with typical longwall subsidence behaviour.This paper discusses the ground movements measured, monitoring methods used, safety measures established, as well as theorising the failure mechanism. Recommendations are made for mine and pit designs for future punch longwall layouts. The paper shows how the movements measured are more aligned to some measurements made during stream valley closure studies previously presented at the ICGCM and challenge the mechanisms suggested by previous literature.BACKGROUNDThe mining of longwalls under or adjacent to large voids (e.g., stream valleys, escarpments, or cliffs) is commonly associated with heavily vegetated or steep surface areas. In areas of extreme topographic variance, access to traditional survey pegs and stations or even new radar or laser technologies is limited. In addition, seldom have the longwall layouts aligned themselves parallel or perpendicular to the surface feature, making interpretation of any available surface movement data more complicated.The punch longwall layout (Figure 1) is also quite uncommon (only undertaken at a small number of longwall mines in Australia) but creates the perfect configuration to enable a somewhat controlled study of the effect of longwall subsidence on a steeply dipping surface feature (an un-vegetated, evenly excavated open-cut highwall). The relatively recently developed radar and laser scanning technology has also enabled near continuous, real-time, sub-millimetre monitoring of a full 500m wide x 100m high highwall. Because the punch longwall enables access and a clear view, the technology could be easily deployed as compared to the highly vegetated and variable topography of stream valleys."

Jan 1, 2018