Modeling Pit Lake Water Column Stability Using Ce-Qual-W2
Total Page:16
File Type:pdf, Size:1020Kb
Modeling pit lake water column stability using Ce-Qual-W2 Lee A. Colarusso & John A. Chermak MFG/Shepherd Miller, Fort Collins, Colorado, USA J. C. Priscu Department of Land Resources & Environmental Sciences Montana State University, Bozeman, MT, USA F. K. Miller Glenn Springs Holdings, Inc., Lexington, Kentucky, USA ABSTRACT: Field measurements of temperature, conductivity, and TDS in an acidic pit lake in southeast Tennessee show that the water column of the lake is currently density stratified. The relative depth of the pit is an intermediate value suggesting that while potentially stable, based on bathymetry alone, the pit’s water column could mix. The potential for the pit’s water column to turn over is a concern because mixing of the shallow low TDS water and the deep high TDS water could potentially degrade the water quality down- stream. Water quality, temperature and bathymetric data were used with meteorological observations to de- velop a limnologic model of the pit using the Ce-Qual-W2 (v.3.0). A series of model simulations were de- signed to test the long-term stability of the pit’s chemical stratification under changing temperatures, storm events, and water quality (TDS) conditions. 1 INRODUCTION known as the hypolimnion averaged 572 mg/L iron, 2966 mg/L sulfate, 14 ug/L copper and 37 mg/L An extensive field investigation occurred at a his- manganese. toric copper mine to characterize the physical and The shallow, low-conductivity water layer dis- chemical characteristics of a 550 million gallon played more variability in temperature and other (2,082,000 m3) stratified pit lake located in southeast water quality parameters than the hypolimnion. This Tennessee. The pit has a maximum depth of ap- type of shallow lacustrine layer is termed a “mixed proximately 60 meters and a surface area of ap- zone” because the water in this zone is affected by proximately 20 acres. A primary tributary enters at meteorological and seasonal fluctuations, changes in the north end of the pit and discharges from the influent water quality and mechanical mixing from south end of the pit. In cooperation with regulatory wind-driven and fluvial currents. The concentra- agencies, this stream was diverted into the pit as a tions of total metals in the pit’s mixed zone were sediment control measure. The boundary for this generally lower than those in the deep water layer field study encompassed approximately 220 acres averaging 3 mg/L iron, 261 mg/L sulfate, 105 mg/L located near the bottom of the tributary’s 1500-acre copper and 3 mg/L manganese. The water quality in drainage basin. the tributary both upstream and downstream of the Field studies conducted at the site between 1998 pit was comparable to the water quality in the mixed and 2002 suggest that the pit is chemically stratified zone of the pit lake. and that the water column is resistant to large-scale mixing known as turnover. Specific conductivity re- sults from depth profiles taken using a Hydrolab® 7 multi-parameter probe demonstrated the presence of 6 a deep high-conductivity water layer underlying a 5 shallow lower conductivity layer (Fig. 1). A distinct 4 chemocline observed at a depth of 7.5 to 8.5 meters 3 separates the low TDS and high TDS layers. 2 1 The results from bi-monthly sampling events Conductivity (ms/cm) conducted from August 2001 through August 2002 0 0 102030405060 yielded relatively constant temperature, pH, dis- Depth (m) solved oxygen, Eh, total dissolved solids (TDS) and conductivity profiles below the top 8.5 meters of the Figure 1. Pit conductivities (measured bi-monthly (8/01-8/02)) pit. The water in the deep, high-conductivity layer Both naturally occurring lakes and pit lakes gen- for the pit are described below. The Ce-Qual-W2 erally fall into one of three stratification categories: model performed volume balance, energy balance, well mixed lakes, seasonally mixed lakes, and lakes mass balance, wind, inflow, outflow, and heat ex- that do not mix (meromictic lakes) (Doyle & Run- change calculations. The pit lake simulations in- nells 1997, Wetzel 1983). While the pit lake ap- voked the CeQual-W2 default equations and used peared to display the chemical stratification of a values for variables recommend by the Ce-Qual-W2 meromictic lake throughout the 2001/2002 moni- User Guide for the evaporation model, transport so- toring period, the pit’s bathymetry suggests that the lution and hydraulic coefficients. The water quality pit’s water column could have the potential to turn constituent calculations were implemented for TDS over. Calculating the surface area to depth ratio in all runs. Ice cover calculations were not used be- known as the relative depth of water body can indi- cause the pit lake, located at temperate latitude and cate the lake or pit’s mixing status (Doyle & Run- at an elevation of 442 meters, does not develop sea- nells 1997, Anderson et al. 1985, Wetzel 1983). sonal ice coverage. Output results included hori- Lakes with a relative depth of greater than 20% gen- zontal and vertical water velocity, TDS, temperature erally do not mix. The classifications of lakes with and density. lower relative depths vary and depend on other fac- Ce-Qual-W2 calibration simulations lasted 1 year tors affecting the stability of a lake’s water column beginning September 2, 2001 and ending August 31, such as density gradients, climatic conditions, and 2002 thereby including the fall and spring transi- wind direction and intensity. While the pit lake is tional seasons as well as the coldest months of the currently stratified, the pit’s relative depth is 18.5%. year when the pit is the most susceptible to turnover. Based on a bathymetry analysis alone, the pit’s wa- The minimum timestep for the 12-month period was ter column has the potential to circulate and turn specified at 10 seconds and Ce-Qual-W2 averaged a over. The long-term stability of the stratification of timestep of 9 seconds over the duration of the simu- the pit lake was further evaluated using a limnologic lations. model. In order to accurately simulate conditions in the Water quality, temperature and previously col- pit lake, a series of initial conditions and time vary- lected bathymetric data were used with meteorologi- ing inputs were specified to reflect field observa- cal data to develop a limnologic model of the pit tions. lake using Ce-Qual-W2 version 3.0. Ce-Qual-W2 is This pit lake is a simply shaped oblong water a hydrodynamic model which provides a two dimen- body approximately 550 by 145 meters. It has sional representation of the mixing dynamics in res- stepped walls and no significant branches. To ervoirs and lakes. Ce-Qual-W2 was originally de- simulate this bathymetry, the pit was divided into 9 veloped by the Army Corps of Engineers, longitudinal segments and 47 horizontal depth layers Waterways Experiments Station, Vicksburg, MS and (Fig. 2). Shallowest layers were 0.5 meters thick. is currently supported by Portland State University Layers in the observed thermocline/chemocline were in conjunction with the Corps of Engineers. 0.1 meters thick. Layers in the deep section of the To calibrate the model for the pit lake, simulated pit were initially 5 meters thick but several 2.5- conditions predicted by Ce-Qual-W2 were compared meter layers were created near the base of the pit for to the temperature and TDS data collected on site better geometrical resolution. between August 2001 and August 2002. Once the simulated conditions were in good agreement with 1 2 3 4 5 6 7 8 9 End View the observed data, a series of model simulations 1 tested the long-term stability of the pit’s chemical 11 stratification under changing water quality (TDS) conditions and during large (high wind) storm Side View events. A discussion of the model’s final calibration Bathymetry 9 longitudinal segments simulation for the pit lake is described followed by a 6 47 depth layers discussion of a series of wind and water quality sen- 5 sitivity simulations. 4 Top View Figure 2. Pit Lake Bathymetry 2 CALIBRATION The initial vertical temperature and TDS profiles for the water column were specified using data from 2.1 Parameters and Initial Conditions field and laboratory measurements. Temperatures were measured using the Hydrolab® probe on The parameters, calculations, initial conditions and 8/14/2001 and ranged from approximately 24°C near time-varying inputs used in the simulation that best the surface to 11°C at depth (Fig. 3). TDS values predicted the observed temperature and TDS profiles were calculated by summing measured ion concen- other Hydrolab probe located at the upgradient trations from discrete water samples taken between weir gauge, which continuously monitored the water August 2001 and January 2002 (Fig. 4). An average temperature of the tributary. TDS concentration was calculated for each depth The concentration of TDS in the pit’s influent interval. The resulting profile ranged from 333 was modeled both at 220 mg/L and 350 mg/L for mg/L TDS in the shallow waters to close to 5000 two separate calibration simulations. A TDS con- mg/L TDS at depth. An initial density profile is not centration of 220 mg/L represents the average meas- specified in Ce-Qual-W2; the modeling routines cal- ured inflow concentration between August 2001 and culate density using temperature and the concentra- July 2002. tion of TDS for freshwater simulations. 30 7-Jan Field 7-Jan Model 25 8-April Field 8-April Model 22-April Field 20 22-April Model 14-Aug Field 15 12-Aug Model Temperature (C) Temperature 10 5 0 0 102030405060 Depth (m) Figure 3.