A 24 h investigation of the hydrogeochemistry of baseflow and stormwater in an urban area impacted by mining: Butte,Montana

Changes in water quality during a storm event were continuously monitored over a 24 h period at a single location along an urban stormwater drain in Butte, Montana. The Butte Metro Storm Drain (MSD) collects groundwater baseflow and stormwater draining Butte Hill, a densely populated site that has been severely impacted by 130 years of mining, milling, and smelting of copper‐rich, polymetallic mineral deposits. On the afternoon of 26 June 2002, a heavy thunderstorm caused streamflow in the MSD to increase 100‐fold, from 0·2 ft³ s⁻¹ to more than 20 ft³ s⁻¹. Hourly discharge and water quality data were collected before, during, and following the storm. The most significant finding was that the calculated loads (grams per hour) of both dissolved and particulate copper passing down the MSD increased more than 100‐fold in the first hour following the storm, and remained elevated over baseline conditions for the remainder of the study period. Other metals, such as zinc, cadmium, and manganese, showed a decrease in load from pre‐storm to post‐storm conditions. In addition to the large flush of copper, loads of soluble phosphorus increased during the storm, whereas dissolved oxygen dropped to low levels (<2 mg l⁻¹). These results show that infrequent storm events in Butte have the potential to generate large volumes of runoff that exceed Montana water quality standards for acute exposure of aquatic life to copper, as well as depressed levels of dissolved oxygen. This study has important implications to ongoing reclamation activities in the upper Clark Fork Superfund site, particularly with respect to management of storm flow, and may be applicable to other watersheds impacted by mining activities. Copyright © 2005 John Wiley & Sons, Ltd.     

Floodplain sedimentation rates in an alpine watershed determined by radionuclide techniques

Vertical profiles of the activities of ¹³⁷Cs and ²¹⁰Pb were measured on floodplain sediment cores and upland soil cores along the Soda Butte Creek and the Yellowstone River to determine floodplain sedimentation rates. The position of mine tailings from a 1950 impoundment failure was used as a stratigraphic marker to estimate the sedimentation rates and to make comparisons with rates provided by radionuclide‐based methods. Mass accumulation (sedimentation) rates calculated from the position of the mine tailings ranged from 0·00 to 0·17 g cm^?2 yr^?1 and were in good agreement with sedimentation rates calculated from the inventories of ¹³⁷Cs and ²¹⁰Pb. Sedimentation rates calculated from the position of the ¹³⁷Cs peak generally overestimated the sedimentation rates, probably because of increased downward migration of ¹³⁷Cs caused by the low pH of water moving through the mine tailings or the high permeability of floodplain sediments relative to upland reference soils. This study demonstrates that the ¹³⁷Cs and ²¹⁰Pb inventory methods for determining sedimentation rates can be applied to an alpine floodplain where sedimentation events are episodic and where orographic effects on precipitation generate strong downstream gradients in the delivery of atmospheric radionuclides. Copyright © 2007 John Wiley & Sons, Ltd.     

Fracture patterns at lava–ice contacts on Kokostick Butte, OR, and Mazama Ridge, Mount Rainier, WA: Implications for flow emplacement and cooling histories

Cooling lava commonly develop polygonal joints that form equant hexagonal columns. Such fractures are formed by thermal contraction resulting in an isotropic tensional stress regime. However, certain linear cooling fracture patterns observed at some lava–ice contacts do not appear to fit the model for formation of cooling fractures and columns because of their preferred orientations. These fracture types include sheet-like (ladder-like rectangular fracture pattern), intermediate (pseudo-aligned individual column-bounding fractures), and pseudopillow (straight to arcuate fractures with perpendicular secondary fractures caused by water infiltration) fractures that form the edges of multiple columns along a single linear fracture. Despite the relatively common occurrence of these types of fractures at lava–ice contacts, their significance and mode of formation have not been fully explored. This study investigates the stress regimes responsible for producing these unique fractures and their significance for interpreting cooling histories at lava–ice contacts.Data was collected at Kokostick Butte dacite flow at South Sister, OR, and Mazama Ridge andesite flow at Mount Rainier, WA. Both of these lava flows have been interpreted as being emplaced into contact with ice and linear fracture types have been observed on their ice-contacted margins. Two different mechanisms are proposed for the formation of linear fracture networks. One possible mechanism for the formation of linear fracture patterns is marginal bulging. Melting of confining ice walls will create voids into which flowing lava can deform resulting in margin-parallel tension causing margin-perpendicular fractures. If viewed from the ice-wall, these fractures would be steeply dipping, linear fractures. Another possible mechanism for the formation of linear fracture types is gravitational settling. Pure shear during compression and settling can result in a tensional environment with similar consequences as marginal inflation. In addition to this, horizontally propagating cooling fractures will be directly influenced by viscous strain caused by the settling of the flow. This would cause preferential opening of fractures horizontally, resulting in vertically oriented fractures.It is important to note that the proposed model for the formation of linear fractures is dependent on contact with and confinement by glacial ice. The influence of flow or movement on cooling fracture patterns has not been extensively discussed in previous modeling of cooling fractures. Rapid cooling of lava by the interaction with water and ice will increase the ability to the capture and preserve perturbations in the stress regime.     

Compositions of magmatic hydrothermal fluids determined by LA-ICP-MS of fluid inclusions from the porphyry copper–molybdenum deposit at Butte, MT

We analyzed 85 fluid inclusions from seven samples from the porphyry Cu–Mo deposit in Butte, MT, using laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS). The Butte deposit formed at unusually great depth relative to most porphyry deposits, and fluid inclusions in deep veins trapped a low-salinity, CO₂-bearing, magmatically derived, supercritical fluid as a single aqueous phase. This fluid is interpreted to be the parent fluid that cooled, decompressed, unmixed, and reacted with wall rock to form the gigantic porphyry Cu deposit at Butte. Few previous analyses of such fluids exist.Low-salinity, aqueous fluids from the earliest veins at Butte are trapped in deep veins with biotite-rich alteration envelopes (EDM veins). These veins, and the Butte quartz monzonite surrounding them, host much of the Butte porphyry Cu mineralization. Twenty fluid inclusions in one EDM quartz vein are dominated by Na, K, Fe (from 0.1 to 1 wt.%) and contain up to 1.3 wt.% Cu. These inclusions contain only small amounts (tens of ppm) of Pb, Zn, and Mn, and typically contain Li, B, Ca, As, Mo, Ag, Sn, Sb, Ba, and W in less than detectable quantities. The abundance of Cu in early fluids indicates that a low-salinity, Cu-rich, aqueous ore fluid can be directly produced by aqueous fluid separation from a granitic magma. Similar inclusions (eight) in an early deep quartz–molybdenite vein with a K-feldspar selvage have similar compositions but contain significantly less Cu than most inclusions in the biotite-altered vein. Analyzed inclusions in both veins contain less than detectable concentrations of Mo even though one is molybdenite-bearing.Low-salinity, CO₂-bearing aqueous fluids are also trapped in pyrite–quartz veins with sericitic selvages. These veins cut both of the above vein types and contain inclusions that were trapped at lower pressure and temperature. Thirty-nine inclusions in two such veins have compositions similar to early fluids, but are enriched by up to a factor of 10 in Mn, Pb, and Zn relative to early fluids, and are slightly depleted in Fe. Many of these inclusions contain as much or more Cu than early fluids, although little chalcopyrite is found in or around pyrite–quartz veins.Eighteen halite-bearing inclusions from three veins from both chalcopyrite-bearing and barren veins with both K-silicate and sericitic selvages were analyzed as well. Halite-saturated inclusions are dominated by Na, K, Fe, and in some inclusions Ca. Whereas these inclusions are significantly enriched in Ca, Mn, Fe, Zn, and Pb, fluids in all three veins contain significantly less Cu than early, high temperature, low-salinity inclusions.Analyses of all inclusion types show that whereas bulk-salinity of the hydrothermal fluid must be largely controlled by the magma, fluid–rock interactions have a significant role in controlling fluid compositions and metal ratios. Cu concentrations range over an order of magnitude, more than any other element, in all four samples containing low-salinity inclusions. We infer that variations are the result of fluid trapping after different amounts of fluid–rock reaction and chalcopyrite precipitation. Enrichment, relative to early fluids, of Mn, Pb, and Zn in fluids related to sericitic alteration is also likely the result of fluid–rock reaction, whereby these elements are released from biotite and feldspars as they alter to sericite. In halite-bearing inclusions, concentrations of Sr, Ca, Pb, and Ba are elevated in inclusions from the pyrite–quartz vein with sericitic alteration relative to halite-bearing inclusions from unaltered and potassically altered samples. Such enrichment is likely caused by the breakdown of plagioclase and K-feldspar in the alteration envelope, releasing Sr, Ca, Pb, and Ba.     

Petrological and geochronological constraints on regional metamorphism along the northern border of the Bitterroot batholith

Quantitative thermobarometry in pelites and garnet amphibolites from the Bitterroot metamorphic core complex, combined with U–Pb dating of metamorphic monazite and zircon from footwall rocks, provide new constraints on the P – T  – t evolution of footwall rocks. The thermobarometric and geochronological results, when correlated with observations from other regions bordering the Bitterroot batholith, define a regional metamorphic history for the northern margin of the Bitterroot batholith consisting of three distinct events beginning with early prograde metamorphism (M1) coincident with arc-related magmatism and crustal shortening at c .  100–80 Ma. Magmatism and crustal thickening led to regional upper-amphibolite facies metamorphism (M2) and anatectic melting between 64 and 56 Ma. Mineral textures related to high-temperature isothermal decompression (M3), coincident with late stages of magmatism in the Bitterroot complex footwall (56–48 Ma), are only preserved in areas adjacent to extensional structures. The close temporal relationship between peak metamorphism and the onset of footwall decompression indicates that thermal weakening was an important factor in the initiation of Early Eocene regional extension and tectonic denudation of the Bitterroot complex and possibly the Boehls Butte metamorphic terrane. The morphology of the decompressional P – T  – t path derived for Bitterroot footwall rocks is similar to other trajectories reported for Cordilleran core complexes and may represent a transition in the deformational style of core-bunding detachments responsible for exhumation.