Open Collections

UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

A hydrogeological and geochemical study of the origin and nature of the prairie flats uranium deposit,.. Rossel, Kathy 1999

You don't seem to have a PDF reader installed, try download the pdf

Item Metadata


ubc_1999-0289.pdf [ 15.2MB ]
JSON: 1.0052927.json
JSON-LD: 1.0052927+ld.json
RDF/XML (Pretty): 1.0052927.xml
RDF/JSON: 1.0052927+rdf.json
Turtle: 1.0052927+rdf-turtle.txt
N-Triples: 1.0052927+rdf-ntriples.txt
Original Record: 1.0052927 +original-record.json
Full Text

Full Text

A HYDROGEOLOGICAL AND GEOCHEMICAL STUDY OF THE ORIGIN AND NATURE OF THE PRAIRIE FLATS URANIUM DEPOSIT, SUMMERLAND, BC by Kathy Rossel BASc, Queen's Univeristy, 1996  A THESIS SUBMITTED IN P A R T I A L F U L F I L L M E N T OF THE REQUIREMENTS FOR T H E D E G R E E OF M A S T E R OF APPLIED SCIENCE in  THE F A C U L T Y OF G R A D U A T E STUDIES (Department of Geological Engineering)  We accept this thesis as conforming to the required standard:  THE UNIVERSITY OF BRITISH C O L U M B I A April, 1999 © Kathy Rossel, 1999  In  presenting  degree  at  this  the  thesis  in  partial  fulfilment  of  University  of  British  Columbia,  I agree  freely available for copying  of  department publication  this or of  reference  thesis by  this  for  his thesis  and study. scholarly  or for  her  Department  of  The University of British Vancouver, Canada  DE-6 (2/88)  Columbia  purposes  gain  shall  requirements that  agree  may  representatives.  financial  permission.  I further  the  It not  be is  that  the  Library  permission  granted  by  understood be  for  allowed  an  advanced  shall make for  the that  without  it  extensive  head  of  my  copying  or  my  written  ABSTRACT An investigation of the hydrogeology and groundwater geochemistry was carried out at the Prairie Flats surficial uranium deposit in Summerland, B.C. The deposit contains an estimated 230 tonnes of uranium, most of which is concentrated in the upper half-metre of soil. It is believed to have formed out of the discharge of uraniferous groundwaters into organic-rich sediments over a period of 10,000 years. Objectives of this study were 1) to trace the origins of the groundwaters transporting uranium into the site, 2) to measure current rates of groundwater discharge and uranium deposition, 3) to identify the mechanism(s) of uranium retention, such as adsorption, reductive precipitation or evaporative precipitation, and 4) to comment on" the likelihood of uranium remobilization. A literature review was first carried out to characterize the area's local and regional groundwater systems. Next, a network of 13 piezometers was installed across the site, with completion depths ranging from 1 to 3 metres. Using these, hydraulic conductivities of the hydrostratigraphic units were measured, and relative head values were monitored at four different times of year. Groundwater and surface waters were tested for pH, Eh, conductivity, and concentrations of U , Ca, Mg, Na, K, N 0 , H C 0 , S0 , and Cl. 3  3  4  Measurements of groundwater discharge into Prairie Creek, which crosses the site, were also carried out.  Results show that the flats are a discharge zone for locally recharged groundwaters, that is groundwaters that infiltrate within a few kilometers of the site and travel at depths of less than 100m within glacial deposits and shallow bedrock. These groundwaters are neutral in pH, relatively oxidizing, and enriched in calcium and bicarbonate. Discharge rates are on the order of 9450 m /year, most of which flows vertically upward from below the deposit. As uranium 3  concentrations in the incoming groundwaters are up to 100u.g/L, current uranium deposition rates are estimated to be around 1 kg/year. This is at least ten times lower than that calculated using the estimated size and age of the deposit, which suggests that uranium deposition rates were higher in the past than they are today. A major fraction of the uranium is held by adsorption to organics, however desorption by the formation of soluble complexes with bicarbonate is also evident. Two field observations show that soil aeration.or exposure to septic discharge may also remobilize uranium. Uranium which is not held by adsorption is precipitated as a reduced uranium mineral, probably U0 ( ). Relatively more reducing conditions near ground surface than 2  C  at depth may help to explain the high concentrations of uranium within the top half-metre of soil.  ii  TABLE OF CONTENTS Abstract List of Tables List of Figures List of Plates Acknowledgements  ii v vi x xi  1  INTRODUCTION 1.1 Site description 1.2 Depositional controls  1 2 7  2  REGIONAL GEOLOGY Introduction 2.1 Bedrock Geology 2.2 Sources of Uranium 2.3 Surficial Geology Summary  10 11 11 17 18 21  3  REGIONAL AND L O C A L GROUNDWATER SYSTEMS Introduction 3.1 Regional Groundwater Flow 3.2 Local Groundwater Flow Summary  22 22 22 28 33  4  RESULTS OF PRAIRIE FLATS FIELD INVESTIGATION Introduction 4.1 Outline of FieldProgram 4.2 Topography 4.3 Hydrostratigraphy 4.4 Hydraulic Conductivity 4.5 Flow Paths 4.6 Seasonal Variations 4.7 Stream Measurements 4.8 Geochemistry Summary  34 34 34 36 36 40 44 44 51 52 52  5  DISCUSSION O F G R O U N D W A T E R D I S C H A R G E AND DEPOSITIONAL HISTORY Introduction 5.1 Source of Prairie Flats Groundwaters 5.2 Quantifying Groundwater Discharge 5.3 New Insights on Depositional History Summary  57 56 56 59 63 66  in  6  URANIUM GEOCHEMISTRY Introduction 6.1 Mobile Uranium 6.2 Immobile Uranium Summary  67 67 67 68 74  7  DISCUSSION O F G R O U N D W A T E R G E O C H E M I S T R Y AND U R A N I U M DEPOSITIONAL CONTROLS Introduction 7.1 Controls on the Groundwater Geochemistry 7.2 Controls on Uranium Fixation 7.3 Controls on Soil Uranium Distribution 7.4 Potential for Uranium Remobilization Summary  75  8  75 75 89 102 103 106  G E N E R A L CONCLUSIONS AND SUGGESTIONS FOR FUTURE W O R K  107  REFERENCES  110  APPENDICES A H Y D R O G E O L O G I C A L FIELD METHODS Piezometer Installation Hydrogeological Tests Stream Measurements  116 116 124 127  B  130  PIEZOMETER LOGS  C H Y D R A U L I C CONDUCTIVITY TEST D A T A  143  D W A T E R S A M P L I N G FIELD M E T H O D S Piezometer Preparation Sample Collection Field Analyses  156 156 158 158  •  E L A B A N A L Y S E S OF C H E M I C A L P A R A M E T E R S Sample Preservation and Storage Lab Instrumentation and Analysis Quality Control Results F G E O C H E M I C A L M O D E L I N G USING PHREEQC A N D WATEQ4F D A T A B A S E  162 162 163 . 1 6 7 172  iv  LIST OF TABLES Table 2-1 Bedrock lithologies in the Summerland and White Lake basins  Table 4-1 Field and laboratory methods used for geochemical analyses Table 4-2 Summary chart of horizontal hydraulic conducitivities calculated using 3 methods Table 4-3 Vertical hydraulic gradients measured across Prairie Flats at different times of year Table 4-4 Hydraulic gradient and seepage results for Prairie Creek Table 4-5 Geochemical sampling results from September, '97 and March, '98  Table A - l Grab sample descriptions Table A-2 Head measurements taken between July '97 and March '98 Table A-3 Vertical gradient data in Prairie Creek Table A-4 Seepage meter data in Prairie Creek  Table E - l Quality control results Table E-2 Error bars assigned to each measured chemical parameter  Table F - l PHREEQC output for Prairie Flats groundwaters Table F-2 Output of PHREEQC sensitivity analysis on U mineral saturation  LIST OF FIGURES Figure 1-1 Location of the (a) Summerland and White Lake basins and (b) Prairie Flats uranium deposit Figure 1-2 Site map showing the known extent of the Prairie Flats uranium deposit Figure 1-3 Histogram of precipitation data measured at the Summerland Agricultural Research Station, 1995. Figure 1-4 (a) Plan view of uranium distribution in uraniferous layer in lbs UsdVm , modified from Culbert (1979), and (b) Uranium distribution in cross-section, modified from Culbert and Leighton (1988) Figure 1-5 Variations with depth of (a) uranium content and organic content (measured by Loss on Ignition, or LOI) and (b) U / Th activity ratio, for two peat core samples taken from Prairie Flats, from Levinson et al. (1984) Figure 2-1 Geology of the Summerland basin, from Jessop and Church (1991) Figure 2-2 Location of the Summerland and White Lake basins, showing approximate locations of cross-sections Figure 2-3 Diagrammatic cross-section of the Summerland basin, from Jessop and Church (1991) Figure 2-4 Cross-section of the White Lake basin, from Church (1973). Symbols: Z = pre-Tertiary basement rocks; 0= Springbrook Formation; X = Undivided Matron Formation; 1= Yellow Lake Member, Matron Formation; 3 = Kitley Lake Member, Matron Formation; 5 = Nimpit Lake Member, Matron Formation; 6 = Park Rill Member, Matron Formation; 7 = Marama Formation; 8 = volcanic conglomerate, sandstones, and shales, White Lake Formation; 0 = mostly volcanic breccias, White Lake Formation; 10 = Skaha Formation.  Figure 2-5 Surficial geology surrounding the Prairie Flats, modified from Nasmith (1962) and K v i l l (1976)  Figure 3-1 Regional and local groundwater flow paths in (a) aerial view and (b) geological cross-section, modified from Jessop and Church, 1991 Figure 3-2 Hydrogeological features in the Summerland basin Figure 3-3 Piper plot of regional groundwaters sampled in the Summerland and White Lake basins. Data from Piteau & Associates (1984) and Grant and Michel (1983)  Figure 3-4 Piper plot of local groundwaters in the Summerland and White Lake basins Data from Grant and Michel (1983) and Piteau & Associates (1984)  Figure 4-1 Piezometer and stream measurement locations on the Prairie Flats Figure 4-2 Topography of the Prairie Flats Figure 4-3 Isopach map showing thicknesses of peat & clay unit Figure 4-4 Frequency distribution of all calculated horizontal hydraulic conductivities Figure 4-5 Cooper-Bredehoeft plots of slug test data in a) shallow piezometers (peat & clay) and b) deep piezometers (till). Ho is the initial hydraulic head value, and H is the subsequent hydraulic head value after start of the test Figure 4-6 Plot of infiltration test data in peat & clay unit. Ho is the initial hydraulic head value, and H is the subsequent hydraulic head value after start of the test Figure 4-7 Contour plot of relative head values measured in shallow piezometers (peat & clay) on Aug 1, '97 Figure 4-8 Contour plot of relative head values measured in the deep piezometers (till) on Aug 3, '97 Figure 4-9 Contour plot of relative head values measured in the deep piezometers (till) on Sept 26, '97, the morning after a rainfall event Figure 4-10 Summary diagram of Prairie Flats hydrostratigraphy and groundwater flow patterns during summer months  Figure 5-1 Piper plot of Prairie Flats groundwaters and nearby springs and surface waters. Shaughnessy springs data is taken from Piteau & Associates (1984) Figure 5-2 Components of recharge zones and discharge zones in a hydrologic budget, modified from Freeze and Cherry (1979, p. 3)  Figure 6-1 A Langmuir adsorption isotherm  Figure 6-2 Eh-pH diagram for aqueous species and solids in the system U-O2-CO2-H2O at 25 °C and 1 bar total pressure. Solid/aqueous boundaries (stippled) are drawn for I U = 10" M . U D C and UTC are U 0 ( C 0 ) " and U 0 ( C 0 ) 3 " respectively. From Langmuir, 1997, p. 505. 5  2  2  3  2  4  2  3  Figure 7-1 Plot of Prairie Flats groundwaters in relation to waters undergoing open system dissolution of calcite with Pco2 of 10" atm 15  Figure 7-2  Relationship between HCO3 and Ca in Prairie Flats groundwaters, suggesting calcite as a source of HCO3  Figure 7-3 Inconsistency between measured Eh and Fe Flats groundwaters  concentrations in Prame  Figure 7-4 a) Sequences of important redox processes at pH 7 in natural systems, from Stumm and Morgan, 1981 and b) Berner's (1971) plot of changes in groundwater composition with depth, with depth scale replaced by an Eh scale according to a) Figure 7-5 Change in concentration of a) NO3, Fe , and SO4 with depth in Prairie Flats groundwaters and b) showing NO3 and F e only at a smaller vertical scale 2+  2+  Figure 7-6 Depth-Eh plot showing relatively reducing groundwater conditions in the peat unit and relatively oxidizing conditions in the underlying till unit Figure 7-7 Increase in concentrations of Ca, Na, Mg,K, surface in Prairie Flats groundwaters  Fe oT, T  and U towards ground  Figure 7-8 Increase in concentrations of C l , N 0 , SO4, and HCO3 towards ground surface in Prairie Flats groundwaters 3  Figure 7-9 Upward concentration of dissolved groundwater species in relation to C l Figure 7-10 Evidence of uranyl-bicarbonate complexation in Prairie Flats groundwaters in (a) March '98 and (b) September '97, inferred from the (c) H C 0 conductivity relationship found in March data 3  Figure 7-11 Evidence of uranyl-carbonate complexation in waters in and around Prairie Flats as opposed to complexation with dissolved organic carbon and PO4 Figure 7-12 Relationship between soil U and dissolved U in a) March '98 and b) September '97 Figure 7-13 Contour plot of dissolved U concentrations (in (J-g/L) in shallow groundwaters underlying Prairie Flats (March '98)  Figure 7-14 Contour plot of dissolved U concentrations (in u,g/L) in shallow groundwaters underlying Prairie Flats (September '97)  Figure 7-15 Eh-pH diagram showing calculated Eh and measured pH in shallow groundwaters underlying Prairie Flats. Also shown are the collapsed stability zones of U30g( ), LUOc^c) and U307( ) C  C  Figure 7-16 Plot of saturation indices with respect to U02( ) and piezometer location C  U^O^c)  at each  Figure 7-17 Effect of Eh, pH, HCO3 concentration and dissolved U concentration on saturation indices with respect to U02( ) and L^Oc^c) C  Figure 7-18 Change in UCfyc) saturation indices with depth in Prairie Flats groundwaters, showing more favourable precipitation conditions near-surface Figure 7-19 Location of septic field and excavation pit for sewage system installation  Figure A - l Diagram of mini-piezometer installation in a streambed Figure A-2 Diagram of seepage meter installation in a streambed  Figure D-l Alkalinity titration data  ix  LIST OF PLATES Plate 7-1 Oily film on puddled water on the flats, thought to be a mixture of ferrihydrite and organic matter  Plate A - l A 1.5m-long piezometers with other lm-long piezometer sections in background  Plate A-2 Rotating dutch auger used to dig boreholes, total length 3m Plate A-3 Gas-powered post-hole digger Plate A-4 Sand bailer with lm-long extensions Plate A-5 3-inch ID P V C casing used during deep piezometer installation. The apparatus used to remove the casing lies on top of it.  Plate D - l Set-up of peristaltic pump (geopump) and flow-through cell (centre) with pH meter attached  x  ACKNOWLEDGEMENTS I would like to state first and foremost that this project was made possible by the combined efforts of many different people. Among them are  •  my field assistants: Khaled Shahin, Craig Nichol, Rina Mackillop, Tracy Bellehumeur, and Melody Taylor  •  friends and acquaintances in Summerland: Coral, Dave, and Dawn Knight, Yvonne and Vince Stelzer, Allen Wiens, Mark Siemens, Dr Tom Northcoat, Bernie Zebarth, George Redlich, Peter Rodd, Ron Mraz, Barry Cowan, Darrol Colgar, Malcolm Lane, Jacob Schulte, Peter Peto, and Tito Canape  •  technicians: Roy Rodway, Doug Poison, Bryon Cranston, and Marc Baker, Bert Mueller, Susan Harper, Paula Parkinson, Maureen Soon, Richard Friedman and Herb Lans  •  U B C faculty and grad students: Dr Ken Hall, Dr Marc Bustin, Dr Lee Groat, Dr Jim Mortensen, Dr Kristin Orians, Dr Mike Church, and Dr Les Lavkulich, Dr K Fletcher, Cari Deyell, Kristina Laretei, Rich Harris, Dr Chris Daughney, Raphael Wust  •  other scientists outside the U B C community: Doug Leighton, David Morley, Dr Fred Michel, Dr Richard Culbert, Dr Robert Zielinski and Dr James Otton  •  my thesis committee members: Myles Parsons, Dr Leslie Smith, Dr Oldrich Hungr, and Dr Greg Dipple  I would like to especially thank my advisor, Dr Roger Beckie, for his encouragement, great ideas, and ever-reliable open-door policy, and Sebastien Tixier for his unfailing love and understanding. This research was funded both by Dr Beckie and the National Science and Engineering Research Council of Canada.  xi  1 INTRODUCTION The Prairie Flats uranium deposit in the Okanagan region of British Columbia is one of the largest and most studied of its kind. It is a secondary, surficial enrichment wherein organic-rich sediments have stripped uranium from discharging groundwaters.  While  studies of uranium concentrations in the host soils, crop cover, and surrounding surface waters and well waters have been carried out, no detailed investigation of the hydrogeology and groundwater geochemistry has ever been attempted.  Therefore,  questions regarding the origins of these uraniferous groundwaters, the mechanisms of uranium entrapment, and the history of the deposit remain unanswered.  This thesis addresses these questions using published hydrogeological and geochemical data from the surrounding region, as well as data collected in a field study at the Prairie Flats. The thesis is divided into eight chapters. Chapter 1 introduces the deposit and summarizes what is known about it to date. Chapters 2 and 3 place the deposit within a more regional geological and hydrogeological context. Chapter 2 maps out the bedrock geology and glacial deposits and highlights the most likely source rocks of uranium. Both the local and regional groundwater flow regimes are traced and characterized in Chapter 3. This information helps to form hypotheses on where the uranium came from, whether it be local or regional bedrock sources.  Chapter 4 presents the results of the hydrogeological and geochemical field study carried out by the author from July 1997 to May 1998. This began with the installation of a network of 13 piezometers across the flats, which provided a better picture of the hydrostratigraphy. Using the piezometers, hydraulic conductivities were measured and hydraulic heads were monitored at different times of year.  Groundwater and surface  water samples were also collected for chemical analysis.  Interpretation of this data is divided into two chapters. Chapter 5 identifies the source of Prairie Flats groundwaters based on their geochemical signature and presents a simplified 1  hydrologic budget for the site. It quantifies groundwater discharge into the deposit and uses this figure to comment on its formational history. Chapter 6 reviews aspects of uranium geochemistry that are relevant to the understanding of uranium mobility in natural environments.  Chapter 7 applies these principles to the geochemical data  collected on site, looking at adsorption, reductive precipitation, and evaporitive precipitation as likely uranium retention mechanisms.  New theories explaining the  spatial distribution of uranium are presented, as well as hypotheses regarding its potential remobilization.  1.1  SITE D E S C R I P T I O N  This section summarizes what was already known about the Prairie Flats uranium deposit before the start of this thesis.  The first section describes the site's general location,  climate, hydrology, stratigraphy and uranium distribution. The second section presents current theories about the nature of uranium fixation, including the conclusions drawn from isotope analyses and bench-scale leaching experiments done on soil cores taken from the flats.  Terrain The Prairie Flats deposit is located in the centre of the Summerland basin in the Okanagan region of B C . It underlies a hay field at the base of Prairie Valley, just south of the Summerland town centre (see Figure 1-1). It is surrounded by a number of small mountains: Cartwright Mountain to the northwest, Mount Conkle to the south and Giant's Head to the east. This topographic setting prompted Culbert and Leighton (1988) to identify it as a "collector basin" type of uranium deposit, where "upwelling groundwaters, surface waters, and runoff meet in a marshy bowl." The flats were formerly a marshland up until the end of this century when drainage ditches were dug across it to drain the land for cultivation (personal communication  2  C o n t o u r s in f t a b o v e s e a l e v e l C o n t o u r i n t e r v a l 5 0 0 ft  Figure 1-1 Location of the (a) Summerland and White Lake basins and (b) Prairie Flats uranium deposit  3  with the Summerland Museum, 1998). These drainage ditches direct Prairie Creek the flats, as is illustrated in Figure 1-2.  Figure 1-2  Site map showing the known extent o f the Prairie Flats uranium deposit  Climate The climate of the Summerland area is generally arid with hot summers and cold winters. Mean annual precipitation, including snow and rain, is around 300 mm (Piteau & Associates, 1984).  Monthly precipitation values measured  at the Summerland  Agricultural Research Station in 1995 were highest in late spring and mid winter, and lowest in late winter and mid fall, as shown in Figure 1-3. Average temperatures are lowest in January (-2.7 °C) and highest in July (20.9 °C) (Piteau & Associates, 1984). The combination of low rainfall rates, warm summer temperatures, and high sun exposure results in a high potential evapotranspiration rate of approximately 652 mm per year (Piteau & Associates, 1984). Potential evapotranspiration rate is the amount of water that would be removed from the land surface by evaporation and transpiration processes i f sufficient water were available in the soil to meet the demand.  70 60 ^  50  co  40  Jan  Feb  Mar  Apr  May Jun July  Aug Sept  Oct  Nov Dec  Figure 1-3 Histogram of precipitation data measured at the Summerland Agricultural Research Station, 1995.  5  Stratigraphy The near-surface stratigraphy of the flats consists of four principal lithologic units: peat, glaciolacustrine silty clay, glaciofluvial till, followed by bedrock.  The peat unit is a  highly-organic mixture of fine silts and decomposing plant material that accumulated before the marshland was drained for cultivation. The silty clay occurs either as lenses within the peat or as a continuous unit below the peat. It was deposited by a lake that covered most of Prairie valley during the late stages of glacial retreat (Nasmith, 1962; Kvill, 1976), which is explained further in section 2.3. The collective thickness of the peat and silty clay units ranges from roughly 0.5 m to 3 m across the flats. It is thickest on the south end of the deposit and thins out to the north, exposing the underlying till unit (see Figure 1-2).  The till unit is probably an extension of the "kame terrace and meltwater channel" deposit (Nasmith, 1962) bordering the north side of the deposit (see section 2.3). Culbert (1980) proposes that this deposit may have been the original barrier which impounded the marsh whose peat now underlies the flats. The till consists of sand and gravel and occasional pebbles, fining downwards. It is called "till" because of its variety of grain sizes, and should not be thought of as hydraulically impermeable, as is often the nature of till. It is estimated to be less than 2 metres thick since drilling records for a well bordering the flats log bedrock at a depth of seven feet (Ministry of Environment, Lands and Parks B.C., 1998).  Hydraulic conductivity measurements of fractured bedrock in the Summerland basin and Trapping Creek Basin (near Kelowna) suggest that the bedrock underlying the till unit is relatively less permeable. In-situ packer testing in a test drillhole for the Summerset Inn on Giant's Head (Golder Associates, 1980) gave a 10" to 10" m/s range in hydraulic 10  conductivity at depths of less than 10m.  7  Lawson's (1968) estimates of hydraulic  conductivity in the Trapping Creek basin were between 10" and 10" m/s at less than 30 9  5  metres' depth. Based on this information, Piteau & Associates assigned bulk hydraulic  6  conductivity values on the order of 10" m/s to bedrock in a regional groundwater flow model of the Summerland basin.  Uranium Distribution Culbert (1980) also states that "it is very unusual to find so thin a topsoil over apparently well-drained till which has concentrated uranium to the extent shown here." In 1980, D.G. Leighton & Associates made approximately 35 auger holes across Prairie Flats and sampled these at every half-meter depth. Soil uranium concentrations were reported as pounds of uranium oxide (as U 0 ) per square metre, and are contoured in plan view in 3  8  Figure l-4a. In total, approximately 250 000 kg of uranium as U 0 3  8  are estimated to  underlie the flats (Culbert, 1979). The formula U 0 is used as a stoichiometric quantity, 3  8  and does not describe the actual form of uranium. Church et al. (1990) also found the uranium to be most concentrated in the east-central part of the deposit.  In vertical profile (see Figure l-4b), the highest concentrations of uranium were found within the top half-metre of soil, except near the southern edge of the site where it is concentrated at a few metres depth (Culbert, 1980). Similarly, a more recent soil survey measured highest uranium concentrations in the 8 to 16-inch depth interval, and a few at greater depth (Ministry of Health B.C., 1981).  1.2  DEPOSITIONAL CONTROLS  The Prairie Flats uranium deposit is a post-glacial, secondary deposit resulting from the mixing of uraniferous groundwaters with organic-rich sediments. Levinson et al. (1984) measured  234  U/  230  T h ratios on soil cores taken from the flats and used these to estimate the  age of the deposit. This method assumes that the decay of  2 3 4  U,  therefore the higher the  234  u  /230  230  T h nuclide is created in-situ by the  T h ratio, the younger the deposit. Ages  were found to increase with depth (see Figure l-5a), and were everywhere less than 10 000 years, calculated from the  234  U / T h ratio and half life of 230  2 3 4  U . Because the uranium  was deposited recently enough so as not to attain equilibrium with it's shorter-lived 7  daughter products, it is not highly radioactive. Many other post-glacial uranium deposits exist in the Okanagan region of B C and the northwestern United States and are further described in Culbert et al. (1984) and Culbert and Leighton (1988).  Figure 1-4 (a) Plan view o f uranium distribution in uraniferous layer i n lbs U 3 0 g / m 2 modified ;  from Culbert (1979), and (b) Uranium distribution in cross-section, modified f r o m Culbert and Leighton (1988)  8  The patchy uranium distribution at the Prairie Flats suggests a complex combination of depositional controls, of which topography, evaporation, and adsorption onto organics are thought to dominate (Culbert et al, 1984).  Levinson et al. (1984) found uranium  concentrations to correlate very closely to the organic contents of the soils (as is measured by Loss on Ignition or LOI in Figure l-5b) which was reaffirmed later by Church et al. (1991).  Levinson et al. (1984) believe that "the uranium was initially  concentrated in the surface layers and, with subsequent burial and compaction, has not been removed from the organic matter." Therefore, they proposed that adsorption is the major process of uranium concentration, minimizing but not eliminating the role of evaporation.  In addition to adsorption, reductive precipitation of uranium is also likely occurring at the flats, however the actual mineral form has yet to be identified. Between 40 and 60% of the uranium is in an unstable, remobilizable form, based on bench-scale leaching experiments of two soil samples (Culbert, 1980). This fraction is retained by adsorption onto organics, while the remaining fraction is held by "precipitation as a micron-sized four-valency mineral, probably uraninite ( U 0  2(c)  ) " (Culbert, 1980).  Because such a significant portion of the uranium is loosely-held, it could be remobilized under the right geochemical conditions.  Culbert (1980) stated that re-release of the  uranium could be initiated by as little as one unit increase in pH of the waters in contact with the peat. Culbert (1980c) also states that " i f [the uranium] was combined with an early spread of nitrate fertilizer, then conditions could exist where this loosely-held uranium would be naturally eluted and would join the drainage runoff, and would report into the lake, or even into the local water supplies." Therefore although the uranium does not impose a radioactive health hazard, it may threaten water quality i f land-use practices in and around Prairie valley are not carefully managed.  9  !  i  DEPTH  i  Figure 1-5 Variations with depth of (a) uranium content and organic content (measured by Loss on Ignition, or LOI) and (b) 238rj/230xh activity ratio, for two peat core samples taken from Prairie Flats, from Levinson et al. (1984)  10  2 REGIONAL GEOLOGY This chapter examines the bedrock and surficial geology of the Summerland basin and identifies possible source rocks of uranium. The geology and hydrogeology of the White Lake basin is also reviewed because of its similarity to the Summerland basin. Due to a lack of well data in the Summerland basin, well data collected in the White Lake basin will be used in upcoming chapters to make inferences about the Summerland basin groundwater regime.  2.1  BEDROCK GEOLOGY  The Prairie Flats deposit lies in the middle of the Summerland basin, which covers approximately 28 square kilometers between Okanagan Lake to the east and the Summerland reservoir to the west (see Figure 2-1). According to Church (1973), crustal extension along the Okanagan valley 50 million years ago caused rift and graben development which tilted and folded these rocks into an elliptical 5 by 10 kilometre synclinal trough. They have also undergone significant fabric loosening by folding and faulting.  Down-faulted volcanic beds meet Giant's Head along the Summerland fault,  which defines the southern boundary of the Summerland basin. Other major faults run approximately north-south along Eneas and Darke Creek Valleys (Figure 1-1).  Both the Summerland and White Lake basins are made up of Tertiary volcanic and sedimentary rocks overlying pre-Tertiary basement granites. The Tertiary volcanic rocks belong to the "Penticton Tertiary Outlier," (Jessop and Church, 1991), which is believed to be a remnant of a continuous belt of volcanic rocks passing through central Washington and south-central B C (Lewis, 1984). Cross sections of both the Summerland basin and White Lake basin are given in Figures 2-3 and 2-4, and their approximate locations in Figure 2-2.  11  LEGEND EOCENE WHITE LAKE FORMATION MOSTLY PCLYMCTIC VOLCANIC CONGLOMERATE AND SANDSTONE WW MINOR CARC0NACE0US SEAMS  i  5502000N  MARAMA FORMATION " ' OACJTE LAVA AND BRECCIA f  (  MARRON FORMATION NIMPrT LAKE MEMBER;TAN TRACHYTE ANO TRACHrANOESTE LAVA AND ASH FLOWS  +  KmEY LAKE MEMBER.-'ucSUY fELO— SPATHC TRACHYANDESTTE LAVA YELLOW LAKE MEMBER: RHOUB-POR( W W  S L  KETTLE RIVER FORMATION MOSTLY GRANITE 80 ULCER CONCLOUERATE ANO BRECCU  PRE—TERTIARY ^] GRANITIC BASEMENT ROCKS  OKANAGAN LAKE  49*33+ 119* 44' Km Figure 2-1 Geology of the Summerland basin, from Jessop and Church (1991)  12  Figure 2-2 Location of the Summerland and White Lake basins, showing approximate locations of cross-sections  13  ca i T3 fi — • 1C ca O &,  c o  S II .S fa" cV2  U5  o  PH  C ca  —  fi O  C  fc o  o C3 o 'S c a o •  ca  >  ii N  O PH  C .fi .2  H J  C  II 00  o 13  CU  s  ca oo O  O  PH  - H  ca S  S  ca  c  oo E  II  • o f^HH Os o  1  tr  °  cu  5 ca .2 h-i PH  >  N  O tl o ra  U 6 o  •H  JO O  ca J 3 C ca h-l jo cu M o ca pi  .-2t—<II  HH £  ca  ^  CH"PH  C  .2  fi fi O fc S  o ^  o  PH  ca  r  E  o  8  o o o  6 o IT)  2$ o 3> « o 10  E  o  8  CD  Is  c 13  1^ ca cu  > ^2 cu JO co  6  Swo X scu•  ^  ca J3  16  The Penticton Tertiary Outlier consists of sedimentary rocks belonging to the White Lake Formation overlying an alternating series of lava flows and ash flows belonging to the Marama and Marron Formations.  At the base of these is the Kettle River Formation,  consisting of conglomerates and breccia. The total thickness of this assemblage is greater than 1000 metres (Jessop and Church, 1991). Underlying the Penticton Tertiary Outlier are the pre-Tertiary basement granites originating from two periods of intrusions. The oldest is Nelson intrusive series (Little, 1961), followed by the Okanagan Highlands Intrusive Complex (Boyle, 1982).  Detailed descriptions of the lithologies of the Tertiary and pre-Tertiary rocks are listed from oldest to youngest in Table 2-1.  This is a compilation of the works of Church  (1982), Little (1961), and Boyle (1982).  Table 2-1 Bedrock lithologies in the Summerland and White Lake basins GEOLOGIC PERIOD TERTIARY  NAME PENTICTON T E R T I A R Y OUTLIER (White Lake Basin) White Lake Formation Marama Formation Marron Formation: Yellow Lake member Kitley Lake member Nimpit Lake member Kettle River Formation  LITHOLOGIES  polymictic volcanic conglomerate and sandstone with carbonaceous seams, shale dome-forming dacite lava and breccias rhomb-porphyry sill feldspathic trachyandesite lavas tan trachyte and trachyandesite lava and ash flows granite boulder conglomerates and breccias  PRETERTIARY (MESOZOIC) Late CretaceousPaleocene  OKANAGAN HIGHLANDS INTRUSIVE COMPLEX  granitoids: quartz monzonite, porphyritic granite, pegmatite  Jurassic-Early Cretaceous  NELSON INTRUSIONS  granitoids: granodiorite, diorite, granite, monzonite  16  2.2  SOURCES OF URANIUM  It is generally believed that the uranium found in surficial uranium deposits dotting the Okanagan Valley is coming from diffuse sources, rather than a single weathering ore body, as stated by Culbert (1980b):  ".. originally we viewed [the deposits] at least to a major extent as so-called geochemical anomalies, and we tried to follow them up until we found that there was just so many of them, that it wasn't necessarily any primary deposit from which they were coming."  Therefore, knowledge of the weatherability and uranium-content of bedrock in this region helps to highlight which rocks are the most likely sources.  Pre-Tertiary Granites In their investigation of surficial uranium deposits in the Okanagan area, Culbert et al. (1984) state that "virtually all deposits always occur in areas underlain by intrusive rocks varying from intermediate to felsic in composition." The Okanagan Highlands Intrusive Complex has the highest labile (ie. leachable) uranium content of all the intrusive rocks in the Okanagan region (Boyle, 1982). These rocks have relatively high carbon contents (Boyle, 1982), which enhances uranium leaching by forming highly soluble uranylcarbonate complexes.  Tertiary Volcanics and Sedimentary Rocks Another potential source of uranium is the Tertiary Volcanics, particularly the Nimpit lake member of the Marron formation. This ash flow sequence has an average uranium content of 6 ppm, and consists of partially welded fragments of potassium-rich glass that can be easily devitrified (Church, 1979; Zielinski, 1981). Furthermore, the Nimpit Lake member is widely exposed and stratigraphically thick throughout the Summerland basin. A l l of the cuttings from a 715 metre deep borehole (EPB/GSC 495) just outside of  17  Summerland (see Figure 3-1) were identified as Nimpit Lake volcanics (Jessop and Church, 1991).  2.3  SURFICIAL G E O L O G Y  The surficial geology of the Summerland basin consists mainly of glaciofluvial deposits, glaciolacustrine deposits, and exposed bedrock. Those deposits underlying or adjacent to the Prairie Flats deposit are shown in Figure 2-5.  Glacial activity which shaped the landscape of the Summerland basin has been described in detail by Nasmith (1962) and Kvill (1976). Nasmith has divided the deglaciation period into 4 distinct stages. During the first stage, a thick ice sheet covered all of Prairie Valley, the town of Summerland, and the main valley that currently hosts Lake Okanagan. During the next stage, the west margin of this ice sheet started to retreat, leaving a lobe of ice extending up Prairie Valley as far as the Summerland reservoir. This ice lobe retreated further east in stage three, and meltwaters obstructed by the local topography ponded to form a proglacial lake covering Summerland and the base of Prairie Valley. During the last stage of glacial retreat, the proglacial lake disappeared and ice sheet in the main valley melted to yield the present-day Lake Okanagan.  In the vicinity of the Prairie Flats there are two classes of glacial deposits: glaciofluvial sands and gravels and glaciolacustrine silts and clays. Each is described further in the following paragraphs.  At the top of Prairie Valley is a series of steep-sided sediment ridges between which the Summerland reservoir currently sits. These consist primarily of coarse sand with pebblesized inclusions. Nasmith (1962) identifies the ridges as a moraine which was deposited directly by the ice lobe that occupied Prairie Valley.  Evidence of sorting and  stratification in this deposit led K v i l l (1976) to call it a kame deposit, which implies that the depositional medium was meltwater, not ice.  Therefore, this is a glaciofluvial  deposit.  18  Pre-Tertiary Granites Tertiary Volcanics and Sedimentary Rocks /}  Glaciolacustrine deposits Glaciofluvial deposits 1 kame deposit 2 kame terrace and meltwater channel 3 fan or outwash terrace (^Z^)  Prairie Flats Uranium Deposit Summerland Reservoir  2  km  3  Figure 2-5 Surficial geology surrounding the Prairie Flats, modified from Nasmith (1962) and Kvill (1976)  19  The rest of Prairie valley including the Prairie Flats is blanketed by a thin layer of glaciolacustrine deposits. Nearby bedrock outcrops and a slightly irregular surface morphology in this area suggest that bedrock is not far below ground surface (Kvill, 1976). These deposits consist of silt and clay that was deposited by the proglacial lake that occupied this area during stage three of deglaciation.  Extending north of the Prairie Flats to the outlet of Eneas Creek (Garnet Valley) is a broad, flat apron of well-rounded, moderately sorted and stratified sands and gravels (Kvill, 1976). K v i l l (1976) identifies it as a "fan," which is a glaciofluvial/alluvial deposit formed when meltwater streams draining upland areas emerge into a larger valley and drop their sediment load in response to the decrease in gradient and widening of the valley. This deposit formed during the last stage of deglaciation after the drainage of the proglacial lake covering Prairie Valley.  Nasmith (1962), however, refers to this "fan" as an "outwash terrace" and places a "kame terrace and meltwater channel" deposit on either side of it (see Figure 2-5). A kame terrace is material that was deposited along the floor of streams occupying the space between a stagnant ice lobe and a valley wall.  Therefore, these "kame terrace and  meltwater channel" deposits were created while blocks of ice were still present at base of Garnet Valley and Prairie Valley, and before the glaciolacustrine silt and clay deposits were complete. Of particular interest is the kame terrace and meltwater channel deposit adjacent to the north side of the Prairie Flats. This may be the same "till- or outwash barrier" referred to by Culbert (1980) (see Chapter 1).  In more recent times, a considerable thickness of peat has accumulated on top of the silt and clay deposits covering the Prairie Flats. Most of this peat accumulated when the Prairie flats was a marshland, receiving water from natural groundwater drainage and runoff coming down Prairie Valley. In the later part of this century, it was drained for cultivation by means of drainage ditches which currently relay Prairie Creek across it (personal communication with Summerland Museum, 1998). To the author's knowledge, no dating of this peat has ever been done, however Levinson (1984) analyzed the  20  U/  Th activity ratios of the uranium on two peat cores and calculated ages of the  uranium ranging from approximately 8000 years (0.5 m depth) to 13000 years (2.5 m) for one core and 1900 (1.0 m) to 18000 years (3.0 m depth) for the other.  Still, it is  impossible to tell the age of the peat since the accumulation of peat and uranium may not have been coincident.  The above history explains the origins of the uppermost stratigraphic units of the Prairie Flats, which are referred to as "peat," "clay," and "till" throughout this thesis. The till unit is glaciofluvial material which was deposited by the lobe of ice that once extended up Prairie Valley. It is probably an extension of Nasmith's "kame terrace and meltwater channel deposit" on the north side of the flats, which precedes the glaciolacustrine silts and clays. Its western extent up Prairie Valley or possible connection to Kvill's "kame" deposit around the Summerland reservoir remain unclear. Further borehole drilling in and around the flats would help to clarify the origin of the till. The "clay" unit comes from the pro-glacial lake and the "peat" unit from more recent detrital accumulation.  SUMMARY The Prairie Flats lie in the middle of the Summerland basin, which consists of faulted lava and ash flow sequences on top of granitic basement rocks. The volcanic rocks are a remnant of the Penticton Tertiary Outlier, which also underlies the White Lake basin. Both the Tertiary volcanic rocks and pre-Tertiary granites are potential sources of uranium to the Prairie Flats. During glacial retreat, coarse-grained glaciofluvial material was deposited at the top of Prairie Valley and just north of Prairie Flats, and silts and clays were deposited under a meltwater lake. The uraniferous peats sit on top of these silts and clays and the underlying till is connected to the glaciofluvial deposits.  21  3 REGIONAL AND LOCAL GROUNDWATER SYSTEMS The Prairie Flats deposit was created by discharging uraniferous groundwaters, however the origins of these groundwaters, whether they be regional or local, remains unclear. This chapter defines regional and local groundwaters in terms of their flow paths and aqueous geochemistry.  Evidence is taken from studies carried out in both the  Summerland and White Lake basins, which is reasonable given their geological similarity.  3.1  REGIONAL GROUNDWATER FLOW  For this thesis, the term "regional" is used for groundwaters recharged several kilometers outside of the Summerland topographic basin and which descend to depths of over 100m in the bedrock.  They either discharge within the basin or into Lake Okanagan (see  Figure 3-1 a). To discharge within the basin, they would have to migrate up through the thick sequence of Tertiary volcanics that overlie the basement granites (see Figure 3-lb).  Boyle (1982) believes that regional groundwaters percolating in the Okanagan Highlands Intrusive Complex (pre-Tertiary granites) are discharging along major fault zones. These basement granites are an ideal setting for a deep-seated, regional groundwater flow regime owing to numerous interconnected faults and fracture systems created by tectonic and intrusive activity (Boyle, 1982). Water samples taken from creeks and boreholes along these fault zones were highly uraniferous (10 to 50 u,g/L), enriched in carbonate (150 to 400 mg/L), and slightly oxidizing (dissolved oxygen 4 to 6 mg/L), which Boyle (1982) assumes to be representative of these regional groundwaters. Such high dissolved oxygen concentrations have been reported in other granitic environments (Leenheer et al., 1974), and can be sustained in the basement rocks of the Summerland basin because of the scarce presence of reductants such as organic matter and sulfide (Boyle, 1982). 22  Cartwright Mtn  Giant's Head  Boundary of Summerland Basin local groundwater flow ^  regional groundwater flow  mmm t  ^  4  4  mmm  t  +  |[  a  |i  *  Bifes  4  4  1  4  Hr „|i,  !B  4i  4F  I  4  «|i  4"V  - "*"«" ' -mmm -  :  16  t  mmm  „|  |„  •  b)  s|s  a|«  1  4 i«|  3:  Hr""""Hr ,|  K  <!  #<  >#  4* "I"  w;> |:  Figure 3-1 Regional and local groundwater f l o w paths in (a) aerial view and (b) geological cross-section, modified from Jessop and Church, 1991  23  Increasing uranium concentrations with distance along Eneas and Darke creeks, which overlie two of these fault zones, support Boyle's hypothesis. Water samples taken from Darke creek and nearby boreholes show increasing uranium concentrations and decreasing tritium counts in the downstream direction. This suggests that older (> 25 years) groundwaters flowing southward along the fault zone are discharging in the downstream sections of the creek. Sampling of Eneas Creek (Ministry of Health B.C. and Ministry of Energy, Mines and Petroleum Resources B . C . , 1981) also found a consistent increase in uranium concentration in the downstream direction, doubling between the Garnet Lake spillway and a site 10km downstream (Ministry of Health B.C. and Ministry of Energy, Mines and Petroleum Resources B.C., 1980 and 1981).  An  increase in concentration of this magnitude cannot be caused by evaporation alone, but is likely due to the discharge of increasingly uranium-enriched groundwaters in the downstream direction.  On the contrary, Piteau and Associates (1984) believed that shallow (10 to 50m deep) groundwater flow systems are recharged in the local upland areas and discharge into Darke, Eneas and Trout creeks; while deeper flow systems are recharged several to tens of kilometers away and discharge into Lake Okanagan or along the lake shore. If deep groundwaters are discharging into these creeks, Piteau and Associates (1984) claim that dilution by shallow groundwaters prevents their chemical identification. They estimate that only 5% of the total groundwater recharge in the Summerland basin goes into the deep, regional flow system, which is comparable to Lawson's (1968) estimate of 2% in the Trapping Creek Basin just south of Kelowna.  Deep boreholes in the Tertiary volcanics of the Summerland basin provide limited evidence of regional groundwater flow.  Two boreholes drilled in the White Lake  formation (depth on the order of 100m) are artesian and are still flowing today (personal communication with driller, Ron Mraz, 1997). These are labeled Mraz well and Mraz well #2 on Figure 3-2. However, the White Lake formation would make a poor regional 24  3 Mraz well #2 4 Faulder Meadow Valley wells  9 Cowan Springs 10 Indian springs  Figure 3-2 Hydrogeological features in the Summerland basin 25  aquifer since it is discontinuous and local in extent (Grant and Michel, 1983). In 1990, a 715m deep well (EPB-GSC 495) was drilled just outside of Summerland in hopes of locating a geothermal groundwater reservoir in the brecciated fault zones of the Kettle River formation (Jessop and Church, 1991). Although brief flows of groundwater and evidence of hydrothermal alteration were found at depth, no hot water reservoir was located.  As part of a geothermal exploration program in the White Lake Basin, Lewis (1984) tested five deep boreholes (on the order of 400m deep) completed in the Tertiary volcanics.  Among these are two artesian boreholes, the P-Well and 78-4, which are  390m and 450m deep, respectively. Measured temperature profiles indicate upwardflowing waters, probably recharged at White Lake and moving up-dip to the west (Lewis, 1984). No large reservoir of hot water was identified, and large differences in water flows and heat flows between boreholes was noted. Therefore these groundwaters appear to be isolated and structurally controlled, confined to fractures sets, faults, the tops of lava flows, and interbeds of conglomerates and sands (Church and Jessop, 1991, Lewis, 1984).  Geochemical sampling of groundwater in these and other deep boreholes (> 100m) in the White Lake Basin have been carried out by Grant and Michel (1983) and by Piteau & Associates (1984), and are plotted in Figure 3-3. Both parties found the groundwater to be older (greater than 30 yrs) and more chemically evolved, based on their isotope ratios, Na-Cl chemistries and high dissolved solids concentrations.  Water from 78-4 had  elevated temperatures (11-14 °C) and a measured Eh of -210 mV, which is typical of waters isolated from atmospheric influence (Piteau & Associates, 1984).  Grant and  Michel (1983) found that the P-Well, 78-4, and three other bedrock wells had very different chemistries, leading them to conclude that deep bedrock groundwaters reflect local rock-water interactions and flow very slowly or not at all. Because of geological similarity, the behaviour of deep-seated, regional groundwaters in the Summerland basin is inferred to be similar to that of the White Lake basin. 26  loo  A  Bedrock wells (depth > 100 m) in White Lake basin Diamond drill holes in White Lake basin (P-well, 78-4)  Figure 3-3 Piper plot of regional groundwaters sampled in the Summerland and White Lake basins. Data from Piteau & Associates (1984) and Grant and Michel (1983)  27  3.2  LOCAL GROUNDWATER FLOW  Local groundwaters are defined as those recharged within or slightly outside the Summerland basin and flowing at depths less than 100m through the bedrock and surficial sediments. The hydrogeological features referred to in this section are found in Figure 3-2.  Small mountains such as Mt Conkle, Giant's Head, Rattlesnake Ridge, and Cartwright mountain (see Figure 3-2) provide suitable recharge areas for shallow groundwater aquifers in the Summerland basin. These mountains are made up of exposed bedrock with well-developed fracture sets (Piteau & Associates, 1984).  Lawson's (1968)  measurements of hydraulic conductivities of shallow volcanic rocks (< 30m) in the Trapping Creek Basin just south of Kelowna were between 10" and 10" m/s. A similar 9  5  range of 10" m/s to 10" m/s was measured by Golder Associates (1980) at depths of up 10  7  to 84 m in a test hole for sewage effluent disposal on Giant's Head.  Most creek valleys in and around Summerland are filled with unconsolidated glacial sediments that comprise ideal aquifers. When drilling the EPB-GSC 495 well at the base of Eneas creek valley, 56 metres of gravel and glacial till were penetrated before encountering bedrock. The Faulder-Meadow Valley district, which runs northward from the junction of Trout and Darke Creeks (see Figure 3-2), sits on top of a gravel aquifer that fills the valley.  A report by Pacific Hydrology Consultants Ltd (1983) states that  "this buried valley acts as a drain with water flowing through the gravel southward and then eastward, probably under Summerland to Okanagan Lake."  The kame deposit  surrounding the Summerland reservoir (see Figure 2-5) may be a connecting aquifer which directs groundwater into Prairie Valley. In addition to natural groundwater flow, significant amounts of water are diverted from Trout Creek Valley to Prairie Valley via the flume connecting Trout Creek to the Summerland reservoir. In turn, this water is transferred to Prairie Valley through irrigation lines and reservoir seepage. 28  These aquifers contribute significant baseflow to local creeks within the basin. A yearlong sampling program of Eneas Creek and a nearby water well gave consistent groundwater to surface water uranium concentration ratios of between 0.8 and 1.30, thereby indicating significant groundwater baseflow contribution to the creek (Ministry of Health B.C. and Ministry of Energy, Mines and Petroleum Resources B.C., 1980). Both also followed parallel cycles in uranium concentration over the year (Ministry of Health B.C. and Ministry of Energy, Mines and Petroleum Resources B.C., 1980). Darke creek and adjacent well waters, as well as Trout creek waters show the same pattern (Ministry of Health B.C. and Ministry of Energy, Mines and Petroleum Resources B.C., 1980).  The following sections present geochemical evidence that the groundwater in these aquifers is locally recharged and can be distinguished from regional groundwaters described in the previous section. Water quality analyses from the following features are found on the Piper plot of Figure 3-4.  Shaughnessy Springs A large groundwater reservoir is believed to occupy the glacial overburden and shallow bedrock channels beneath the Summerland business district and to receive waters from Prairie and Eneas Creek valleys (Piteau & Associates, 1984, Ministry of Environment, Lands and Parks B . C . , 1986).  K v i l l (1976) states that discharge from Eneas Creek  infiltrates the coarse-textured fan material north of Summerland town centre and flows towards the lake as groundwater.  Shaughnessy springs (see Figure 3-2), whose outlet  lies at the base of the silt cliffs just off Highway 97, is believed to be a discharge point for this aquifer. The springs supply the Summerland Trout hatchery with a steady year-round flow rate of 67.5 L/s at a temperature of 11 °C (Piteau & Associates, 1984).  29  100  Wells in overburden: White Lake basin Summerland Landfill monitoring wells Spring waters: White Lake basin Shaughnessy Springs, Summerland Figure 3-4 Piper plot of local groundwaters in the Summerland and White Lake basins. Data from Grant and Michel (1983) and Piteau & Associates (1984)  30  The chemistry of Shaughnessy springs waters points to local sources of recharge, including rainfall, irrigation water, and local creeks. Tritium levels measured in 1986 dated the water to be less than 30 years old, and hydrogen and oxygen isotope ratios were found to match those of local creeks and irrigation waters. (Ministry of Environment, Lands and Parks B . C . , 1986).  Other analyses showed these waters to be relatively  oxidizing (Eh +180 mV), slightly alkaline (pH = 7.90), with moderate amounts of dissolved solids (roughly 300 mg/L), and a C a - H C 0  3  chemistry (Ministry of  Environment, Lands and Parks B . C . , 1997 and Piteau & Associates, 1984).  This  chemical signature is typical of locally recharged groundwaters circulating in shallow glacial sediments that have undergone carbonate-mineral dissolution under open or partially open C 0 conditions (Freeze and Cherry, 1979, p. 284). Since the early 1990's, 2  increases in concentrations of nitrate and coliform bacteria have been noted in the spring water, which probably comes from agricultural runoff and septic discharge in Eneas and Prairie valleys (Ministry of Environment, Lands and Parks, B.C., 1997).  Perpetual Landslide Another example of a surficial aquifer receiving significant recharge from agricultural activity is the "Perpetual Landslide," located on the north side of Trout Creek Canyon. Here unconsolidated Quaternary deposits are in a state of constant creep, made worse in recent years by the increased pore water pressures from diffuse groundwater discharge (Piteau Associates, 1984). Roughly a dozen spring pools were identified on the side of the landslide by the author in May '97. This is likely irrigation water from the nearby golf course and fruit orchards. Church (1980) measured a uranium concentration of 43 Lig/L in some of these springs.  Summerland reservoir A network of shallow monitoring wells was installed in the kame deposit at the top of Prairie Valley just above the town reservoir (Golder Associates Ltd., 1994) to monitor the 31  migration of landfill leachate. Groundwaters were found to be migrating eastward, and groundwater samples taken from a background well were neutral to alkaline (pH 7.5 to 8), with low dissolved solids (230 mg/L), significant Ca (50 mg/L), and minor M g and Na concentrations (9 mg/L each). Uranium concentrations of 1 |ug/L have been detected in the reservoir itself (Church, 1980).  Indian Springs A uranium concentration of 1 Lig/L was measured at Indian Springs, located just west of the town reservoir (see Figure 3-2) (Church, 1980). This spring discharges at a rate of 2 L/s, and is slightly acidic (pH 5.8) with low dissolved solids (conductivity of 170 LiS/cm). The dominant ions were H C 0 and Ca, with minor but equal concentrations of M g and 3  Na. This is believed to represent even shallower groundwaters recharged within a few kilometers of the springs.  Cowan Springs A number of springs surface just downgradient from the Summerland reservoir. One of these is Cowan Springs, located in the backyard of Barry Cowan. These were sampled by the author in March '98 and their analyses are plotted later in Figure 5-1. Compared to Indian springs waters, concentrations of C l were slightly elevated (possibly due to mixing with leachate from the landfill upgradient of the reservoir), but the overall chemistry was the same.  White Lake Basin Groundwaters sampled in surficial wells in the White Lake basin had chemistries typical of shallower, more locally recharged groundwaters than those sampled in bedrock (Grant and Michel, 1983). Tritium counts indicated recharge in the mid 1970's; and their ion composition was predominantly C a - H C 0 and Ca-Mg-Na-HC0 with low concentrations 3  3  32  of dissolved solids.  SUMMARY  For this thesis, regional groundwaters are defined as originating several kilometers west of the Summerland basin and flowing at depths of over 100m within the bedrock toward Lake Okanagan. They are generally older than 25 years and have high levels of dissolved solids.  Borehole testing in the White Lake basin further characterized regional  groundwaters as relatively reducing, enriched in Na and C l , and confined to poorlyconnected fracture sets and porous units deep underground.  However, Boyle (1982)  argues that regional groundwaters form a well-established flow regime rising along major fault zones in and around the Summerland basin, and are relatively oxygenated and enriched in H C 0 . 3  Creeks and springs in the Summerland basin are discharge points for local groundwaters. These are defined as being recharged within or slightly outside the basin and flowing at depths of less than 100m in the shallow bedrock and glacial deposits.  Evidence  supporting their recent recharge and short travel times include their C a - H C 0 chemistries, 3  young age (< 25 years) relatively oxidizing Eh, and low dissolved solids concentrations. If regional groundwaters are discharging within the Summerland basin, they are probably diluted beyond chemical recognition. Therefore, evidence of regional uranium transport into the Summerland basin is scarce.  33  4 RESULTS OF FIELD INVESTIGATION This chapter presents the results of a field program designed to further characterize the hydrogeology and groundwater geochemistry of the Prairie Flats. It begins with a general outline of the field program carried out from July '97 to May '98. This is followed by a more detailed characterization of the topography and hydrogeology of the flats, including measurements of hydraulic conductivities and hydraulic gradients in the peat, clay and till units. Streamflow measurements are also presented, as well as the geochemical sampling results. This data is used in Chapters 5 and 7 to construct hypotheses regarding the origins and nature of the uranium deposit.  4.1  OUTLINE OF FIELD PROGRAM  A network of 13 piezometers was installed across the Prairie Flats in July and August '97 (see Figure 4-1). Four deep piezometers (labelled -01) were completed in the till unit at depths of about 3 metres. Nine shallow piezometers (labelled -02) were completed in the peat and clay units at depths of about 1 metre, which is just below the water table. The deep and shallow piezometer pairs (101/102, 401/402, 601/602, 701/702) were each installed within 50 cm of each other to measure vertical hydraulic gradients and changes in chemistry with depth.  Further details on piezometer installation are given in  Appendix A , and piezometer logs are located in Appendix B .  Slug tests were carried out in the piezometers to measure horizontal hydraulic conductivities of the peat, clay and till units, and a borehole infiltration test was used to measure the vertical hydraulic conductivity of the peat unit.  Hydraulic heads were  measured across the site over four seasons (late July/early August '97, late September '97, early March '98, and mid May '98) to determine lateral and vertical components of the hydraulic gradient and the general direction of groundwater flow.  34  Figure 4-1 Piezometer and stream measurement locations on the Prairie Flats  Measurements of groundwater-surface water interaction were made at locations A , B , C, and D on Prairie Creek using mini piezometers and seepage meters. Field methods for these tests are described in detail in Appendix A , and the slug test data can be found in Appendix C.  Water samples were collected from the piezometers and from surrounding springs and creeks in September '97 and March '98. Parameters tested in September included pH, temperature, conductivity, and dissolved uranium. A more thorough analysis was carried out in March with the addition of ferrous iron, dissolved carbon (total, inorganic, and organic), major cations (Ca, Na, K , Mg, F e  T0T  ) , major anions ( H C 0 , Cl, S0 ) and nitrate. 3  4  Field and laboratory methods are explained further in Appendices D and E , and are summarized in Table 4-1.  4.2  TOPOGRAPHY  The topography of the Prairie Flats is shown in plan view in Figure 4-2, based on a level survey carried out in August, '97. The ground surface rises gently to the southwest with an average gradient of 0.5%.  At the east end of the flats the gradient flattens to about  0.1% and the topographic contours curve around an elevation low. West of the deposit the terrain continues to climb up Prairie Valley to the Summerland reservoir, and about half a kilometer to the east it encounters the steep cliffs of Giant's Head.  4.3  HYDROSTRATIGRAPHY  The near-surface lithology of the Prairie Flats consists of four principal units: peat, clay, till, and bedrock.  These can be grouped into three hydrostratigraphic units:  a low  permeability peat & clay unit, a higher permeability till unit, followed by low permeability, fractured bedrock. Further discussion of the field results will be presented from this hydrostratigraphic perspective. 36  •O  "O  CD  CD  Id  • S-l CD  CO  CD  pi o fl  >H  00 cd i-  15 tn  ? CD  cd 001  CO  CD  &  cd 00  CD  O fl  O fl  CD  CD  cd  a  CD  13 T3  CD  fl  JP  0  o  Id  Pi  ccj  PI  -a  T 3  CD  CD CM  fl T2  jo O  cd  CD  cd  T3  IH  CD +-»  CD  CD  w  S-H  CD  CD  -4—*  CD  -*-»  -+—»  fl  _3  •o"  CD  t-  CD  -4—»  PI  -o  CD  >H  <3  co  0)  M  Id  O  PI  Pi  CD  CD  CD  co"  CD  Pi  to  O fl  O fl  CO  XI  I CD  CD N  CD  CD  H  CD  PI  PI  CD  c  T 3  o  331  CD  CD  Id  fl  CD  CD  O  CD  CO  T 3  X> «+H  00  VO  CD  CD  w  CD  fl  A  •O O XI  -i—> CD  PI  _o o pi  .2 g & 'co O co  T 3  CD  CD  fl -*-»  o  CD  o  cd cd  K  CD  M CD  O  fl  CD  G O PH  PH  13 PI o  CD  1 + O  CH  PH  CN CD  CD  W) 00 fl w  cd  U  O z z cd  g  cd  fl o  6b fl o O  .2  o  g  C^  s  CO ' f l  •o  ^ -o fl  uI PI  U o  o U  s PH  u  u o o° PH  OO  rj  9 <-> 55  H  3}  Total thicknesses of the peat & clay unit encountered during this field program ranged from 0.5 to 2.5 m, and are contoured in Figure 4-3. The base of the till unit is below the deepest piezometer, therefore the total thickness of the till unit was not measured.  contours in m relative elevation  Figure 4-2 Topography of the Prairie Flats  38  Contours in m thickness of peat & clay  Figure 4-3 Isopach map showing thicknesses of peat & clay unit  31  4.4  HYDRAULIC CONDUCTIVITY  Horizontal Hydraulic Conductivity Slug tests were performed to measure the horizontal hydraulic conductivity, K , of the H  peat & clay and till units. This data can be found in Appendix C. The data was analyzed by the Hvorslev (1951), Cooper-Bredehoeft-Papadopulos (1967), and Bouwer-Rice (1976) methods (using AQTESOLV™ for the latter two).  A l l three approaches were  tried since the site doesn't perfectly fit any of their sets of geometrical conditions. These discrepancies are discussed below.  The Hvorslev method assumes a fully completed well in a homogeneous, isotropic, perfectly confined aquifer of infinite lateral extent and unlimited thickness. However, the peat & clay unit is unconfmed and the piezometers are not fully completed (ie. screened across the entire thickness of peat & clay).  Also, both the Bouwer-Rice (1976) and  Hvorslev methods neglect specific storage. Demir and Narasimhan (1994) pointed out that specific storage effects cause head data on a semi-log plot to deviate from a straight line (ie. bend concave upward), which was observed at a few piezometers (202, 402, and 902).  The Cooper- Bredehoeft-Papadopulos (1967) method assigns a specific storage  value to the data, but the assumptions of Hvorslev are also inherent in this method.  Despite the limitations of these three methods, calculated hydraulic conductivity values were within the same order of magnitude for both the peat & clay unit (10" m/s) and the 7  till unit (10~ m/s), as shown in Table 4-2. Furthermore, the frequency distribution of all 5  hydraulic conductivity values (see Figure 4-4) is bimodal, with peaks corresponding to the peat & clay unit (3xl0" m/s) and till unit (6xl0" m/s). 7  6  ^0  Table 4-2 Summary chart o f horizontal hydraulic conductivities calculated using 3 methods  101 101 601 601  Lithology at screen intake till till till till  102 202 302 402 602 602 602 802 902 902  peat, clay till clay peat clay clay clay peat clay clay  zometer  Hvorslev K(m/s) 7.0E-07 1.1E-05 2.6E-05 2.1E-05  Cooper & Bredehoeft K(m/s) 3.7E-07 6.9E-06 2.1E-05 2.2E-05  Bower & Rice K(m/s) 1.4E-06 3.3E-05 4.5E-05 6.9E-05  2.0E-07 2.5E-07 4.7E-07 6.6E-07 6.7E-07 1.2E-06 1.0E-06 6.0E-07 1.5E-07 3.5E-07  1.7E-08 n/a 1.0E-06 1.8E-06 1.1E-06 1.3E-06 1.3E-06 8.9E-07 1.4E-07 4.5E-07  3.7E-07 n/a 8.6E-07 1.5E-06 1.3E-06 2.2E-06 1.9E-06 1.1E-06 5.0E-07 5.9E-07  12  10  o c  0 3  CU  Ln K Figure 4-4 Frequency distribution o f all calculated horizontal hydraulic conductivities  4-1  The alignment of head recovery data in the Cooper-Bredehoeft plots provides insight on the homogeneity of the peat & clay unit and the till unit. Figure 4-5a suggests that hydraulic conductivities of the peat & clay unit are relatively uniform across the site. Deviation from this trend at piezometers 102 and 202 may be due to the persistence of fines in the piezometer after development. The hydraulic conductivities measured in till unit (see Figure 4-5b) are not as uniform, which is reasonable given the unsorted nature of till. a)  1.2 o PZ-102 • PZ-202 A  PZ-302  X  PZ-402  o PZ-602 X PZ-802 + PZ-902  10  100  1000  10000  100000  Time (sec)  b)  <>PZ-101 • PZ-101  o X X  100  1000  10000  A  PZ-601  x  PZ-601  100000  Time (sec) Figure 4-5 Cooper-Bredehoeft plots of slug test data in a) shallow piezometers (peat & clay) and b) deep piezometers (till). HQ is the initial hydraulic head value, and H is the subsequent hydraulic head value after start of the test  42  Representative hydraulic conductivity values of 6x10" m/s for the peat & clay unit, and 7  2xl0" m/s for the till unit were calculated from the curves of best fit to the Cooper5  Bredehoeft-Papadopulos plots. These values are considered to be the best estimates of hydraulic conductivity because they account for specific storage effects.  Hydraulic  conductivities calculated using the Bouwer-Rice method were consistently higher than those calculated using the other two methods and thus are likely to be less representative.  Some piezometers were slug-tested repeatedly in order to detect any dependence of head recovery on the initial displacement, H , or to detect the possible evolution of a low 0  permeability well-skin, as is recommended by Butler et al. (1996). On the contrary, hydraulic conductivity values tended to increase, which means that some piezometers may not have been optimally developed before slug testing.  Vertical Hydraulic Conductivity A borehole infiltration test was carried out to estimate the vertical hydraulic conductivity, K of the peat & clay unit. The test method is described in Appendix A and the resultant v  head data is plotted in Figure 4-6.  The Hvorslev method was used to interpret the data since it is only when measuring the horizontal hydraulic conductivity (ie. radial-flow across the piezometer screen) that the neglect of specific storage incurs significant error (Demir and Narasimhan, 1996). A shape factor appropriate to the flow geometry was chosen (case C: "flush bottom in uniform soil" Hvorslev, 1951) and the mean hydraulic conductivity was calculated using:  K ^ T I D / I I T  where K,,, = V K X V K H  v  (4.1)  is the mean hydraulic conductivity (cm/sec)  K and K are the horizontal and vertical hydraulic conductivities (cm/sec) H  v  D is the diameter of the casing (cm) T is the basic time lag (sec)  45  0  50000  100000  150000  200000  Time (sec) Figure 4-6 Plot of infiltration test data in peat & clay unit. Ho is the initial hydraulic head value, and H is the subsequent hydraulic head value after start of the test  To calculate K^, a representative K value of 6xl0~ cm/sec was used for the peat & clay 5  H  unit.  This gave a K of about 4.4 xlO" cm/sec or 4xl0" m/s, about two orders of 7  9  v  magnitude lower than the horizontal hydraulic conductivity. A lower K is reasonable v  given that this material is made up of peat intermixed with silt and clay that has been repeatedly compacted by the heavy farming machinery. However, 10" m/s is an order of 9  magnitude lower than the 10" assigned to the underlying bedrock by Piteau & Associates 8  (1984), and corresponds to the hydraulic conductivity of shale (10"" to 10" , Freeze and 7  Cheery, 1979, p. 29), which is more compact and impermeable than peat. Therefore, 10"  8  m/s was chosen as a more reasonable value for the K of the peat & clay unit. v  4.5  F L O W PATHS  Contour plots of measured head values enable the lateral component of the hydraulic gradient to be calculated and the general groundwater flow direction to be traced.  44  Lateral Flow in the Peat & Clay and Till Units A representative contour plot of the hydraulic heads in the peat & clay unit is shown in plan view in Figure 4-7. This head pattern changed minimally throughout the year. By comparing Figure 4-7 and Figure 4-2, it is evident that the head contours closely mimic the topographic contours. The lateral component of groundwater flow is generally west to east, and the gradient ranges from 0.009 on the west end of the flats to 0.001 on the east end.  Because of the different completion depths of the deep piezometers, hydraulic heads measured here offer only an approximate measure of the lateral component of groundwater flow through the till unit, which is illustrated in Figure 4-8. The lateral component of the hydraulic gradient in the till is about half that in the peat & clay because of the higher hydraulic conductivity of the till.  Measurements taken the morning after a rainfall event showed a surprising rearrangement of heads in the till unit (see Figure 4-9), which returned to normal the next day. No such perturbation was noted in the peat & clay unit. This quick response suggests that the till unit is connected to nearby sources of recharge, possibly storm runoff that collected on residential streets and infiltrated along the northern boundary of the flats where the peat & clay unit is thinner.  46  \  Approximate direction of lateral component of groundwater flow o  Shallow Piezometer  •  Deep Piezometer  0  100  200  300 m  Contours in m relative head  Figure 4-7 Contour plot o f relative head values measured in shallow piezometers (peat & clay) on A u g 1,  '97  4<<  Approximate direction of lateral component of groundwater flow o  Shallow Piezometer  •  Deep Piezometer  100  200  o  rv  -8.  \ o  300 m  Contours in m relative head  Figure 4 - 8 Contour plot o f relative head values measured in the deep piezometers (till) on 3,  '97  Contours in m relative head  Figure 4-9 Contour plot o f relative head values measured in the deep piezometers (till) on Sept 26, '97, the morning after a rainfall event  Vertical Flow in the Peat & Clay Unit Table 4-3 summarizes the measurements of the vertical component of the hydraulic gradient taken at the four nested pairs of piezometers. These gradients were calculated from the head difference between a point just below the water table in the peat & clay unit (shallow piezometer) and a point at the top of the till unit (deep piezometer), over a difference in elevation of at least 110 cm, (except for 28 cm at 601/602). Given the difference in hydraulic conductivity between the peat & clay unit and the till unit, most of the head losses occur across the peat & clay, therefore these measurements represent the vertical component of the hydraulic gradient across the peat & clay unit. Gradients are described as upward where the head at depth exceeds that above, thus where groundwater would tend to flow upwards across the peat & clay unit (and vice versa).  Unlike the lateral component of the hydraulic gradient through the peat & clay unit, temporal variations were noted in the vertical component. The piezometer pairs located in the middle of the deposit, 401/402 and 601/602, show both downward and upward hydraulic gradients. The downward hydraulic gradients are probably an expression of recent climatic events, such as rainfall or snowmelt, rather than seasonal trends. However, upward gradients persist throughout most of the year, and average 0.04. More frequent head measurements are necessary to confirm seasonal and shorter-term trends.  Piezometers 101/102 and 701/702 show consistent upward gradients at all sampling times, averaging 0.04 and 0.10 respectively. 101/102 are located near an impermeable clay liner that was put in place under the Giant's Head School playfield to prevent ponding during wet periods. Groundwaters that formerly ponded in the middle of the field are now forced to surface along its edge, creating these upward gradients. Piezometers 701/702 are located near the base of a hill, where springs have been located by local residents (personal communication, 1997, 1998) and by Piteau & Associates (1985). Therefore these upward gradients are probably topographically-driven.  0.02  o  cd  0.01  O  -0.07  0.11  -0.17  -0.02  0.08  -0.14  -0.04  n/a  -0.15  -0.02  0.05  CD o  CN O  CN O  d  d  d  VO o  ^o  ro O  ro O  fl  d  d  d  IS  o d  CN O  ro O  VO o  ro O  CN O  d  d  o  VO o  c3  oo CN I  CM  o  -0.06  0.02 0.03 0.02  ro O  0.08 0.08  vo  0.16  0.15  -0.03 -0.02 -0.05 -0.04  -0.01 -0.01 n/a n/a  0.04  5?  0.04  CN  CD 0 0  CN i  CM  CD CO  CN i  OH  CD 0 0  0.02 0.04  <  o  CN O  o o  o  CN ro CN  1-1  d  d  d d  n/a  0 0  d  VO o  n/a  0 0  d  n/a  ro O  i  n/a  g  n/a  bO  n/a  ro  CN  CN o vo  CN O  o  o  o o  VO  *—H  4.6  STREAM MEASUREMENTS  Prairie Creek flows through a series of open drainage and underground culverts from entry point (E) to the Prairie Flats to exit point (A) (see Figure 4-1). Mini piezometers and seepage metres were installed in the creek bed at points B , C, and D to measure the vertical component of the hydraulic gradient and the seepage rate between the creek and the underlying groundwater. These field methods are further described in Appendix A , and their results are listed in Table 4-4.  Table 4-4 Hydraulic gradient and seepage results for Prairie Creek  Location IB C D  Date July 19, 1997 July 22, 1997 July 20, 1997  Vertical Hydraulic Gradient 013 -0.06 -0.32  Seepage Rate (m yf'm" ) 25 -10 not measured 3  2  * positive gradients and seepage rates indicate upward discharge of groundwater into the stream  From the above data set, it is difficult to state in general terms whether groundwaters contribute to streamflow across the flats or vice versa. Vertical hydraulic gradients and seepage rates varied significantly in both magnitude and direction between locations B , C, and D, and possibly across the rest of the site. This may be attributed to natural variations in creek bed topography or hydraulic conductivity, or to changes in the natural hydrology due to culvert installation or nearby residential development. For example, the relatively large upward gradient and seepage rate measured at B may be explained by its proximity to the impermeable liner under Giant's Head School playing field.  5|  4.7  GEOCHEMISTRY  Results of geochemical sampling in and around the Prairie Flats are found in Table 4-5. The waters are generally neutral, with pH values between 6.5 and 8.0. They are enriched in Ca and H C 0 , with minor Na. Uranium concentrations span a large range, from 10 to 3  nearly 1000 ppb. N 0 and C l may be coming from fertilizer, septic discharge, or landfill 3  leachate (from the Summerland landfill behind the town reservoir) migrating eastward down Prairie Valley.  Further discussion of these results is given in the upcoming  chapters.  SUMMARY  Borehole logs, slug tests, hydraulic head measurements, and a water sampling program provide further insight into the groundwater regime underlying the Prairie Flats. The hydrostratigraphy consists of a low-permeability peat & clay unit 0.5 to 2.5 metres thick, followed by a higher permeability till unit of unknown thickness.  Average lateral  hydraulic conductivity measurements in the till unit were around 2xl0" m/s, and 6xl0" 5  7  m/s in the peat & clay unit. A vertical hydraulic conductivity of 10" m/s was assigned to 8  the peat & clay. While the lateral component of hydraulic gradient in the peat & clay unit was small (0.005) and relatively unchanging throughout the year, the vertical component of the hydraulic gradient was an order of magnitude greater (0.04) and more variable. This is attributable to the varying effects of evapotranspiration, rainfall, and snowmelt at different times of year.  A summary sketch of the hydrostratigraphy and summer  groundwater flow conditions across Prairie Flats is given in Figure 4-10.  5z  Table 4-5 Geochemical sampling results from September, 1997 and March, 1998  Location  Sampling Exposed Depth (m) unit  PZ-101  2.58  PZ-401  2.47  PZ-601  1.24  PZ-701  2.29  PZ-102  0.98  PZ-202 PZ-302 PZ-402  0.88 1.70 1.12  PZ-502  1.16  PZ-602  0.97  PZ-702  1.14  PZ-802  1.14  PZ-902 Mraz well  1.06 100  Prairie Creek at E Prairie Creek at A Shaughnessy Springs Cowan Springs Indian Springs  till  Sampling Date  Mar 05 1998 Sep 25 1997 till Mar 07 1998 Sep 25 1997 till Mar 05 1998 Sep 24 1997 till Mar 06 1998 Sep 26 1997 peat, clay Mar 05 1998 Sep 24 1997 till Mar 06 1998 clay Mar 07 1998 peat Mar 07 1998 Sep 25 1997 peat Mar 06 1998 Sep 26 1997 clay Mar 05 1998 Sep 24 1997 peat Mar 06 1998 Sep 26 1997 peat Mar 07 1998 Sep 26 1997 Mar 05 1998 clay bedrock Mar 05 1998 Sep 25 1997 Mar 04 1998 Sep 24 1997 Mar 04 1998 Sep 25 1997 Mar 04 1998  pH  7.84 7.49 7.72 6.66 6.7 7.19 7.68 6.8 6.83 7.02 7.18 7.19 6.82 6.58 7.3 6.8 6.84 7.05 7.53 7.8 6.96 6.85 6.04 8.2 7.83 8.09 7.89 7.57 7.25  Sep 24 1997 Mar 06 1998 7.68 Mar 06 1998 7.3  Total Fe Temp Cond Meas Fe mg/L mg/L C uS/cm Eh mV 153.6 0.0 0.0 8.6 350 14.6 250 7.8 670 233.9 0.0 0.1 15.8 1300 1.2 640 154.5 1.8 6.9 15.5 420 0.0 625 358.7 <0.1 6.6 16.5 530 5.2 3.2 7.0 880 191 18 420 3.6 7.5 5.9 3300 155.6 500 263.9 0.7 0.2 6.5 1440 231.4 3.8 12.0 6.5 6.7 16.6 1550 7.4 650 196.6 2.7 9.5 16 660 2.1 3.1 6 750 198.8 16.8 490 0.4 0.3 6.7 670 243.3 17.2 720 9.5 4.8 250 222.1 3.3 15.3 340 264 3.5 7.5 6.1 270 0.4 6.4 1010 92 0.8 17.2 940 400 377 0.0 7.8 2+  7.7 15.4 12.8  540 360 560  437  7 7.4  670 120  368.5 393  HC0 mg/L 3  239 604 549 464 525 1354 610 647 498 573 555 203 88 732 283  0.0  345  0.0  0.0  317  0.0 0.0  0.0 0.0  488 83  380  53  Table 4-5 Continued Geochemical sampling results from September, 1997 and March, 1998  Location  TC mg/L 31  IC mg/L 30  PZ-601  103 89 91  59 50 74  PZ-701  70  69  PZ-102  97  72  25  PZ-202 PZ-302 PZ-402 PZ-502  250 42 142 99 76  93 39 87 67 62  157 3 55 32 14  PZ-602  106  83  23  PZ-702  81  79  PZ-802  41  30  PZ-101 PZ-401  so  OC N0 P0 4 mg/L mg/L mg/L mg/L 1 0.3 0.9 28.9 0.1 44 0.4 0.8 222.9 39 0.6 17 0.4 1.7 42.9 0.4 1 6.7 0.8 50.1 3  14.5 1.4 43.2 3.6 5.0 1.9 14.5  4  K Na U mg/L mg/L ug/L 12 11.9 2.25 14 46.2 64.0 1.86 743 7 30.0 2.30 36 35.7 57 41.3 70.0 4.75 81 166 42.4 61.0 3.05 74 31 165.5 630.0 10.10 3961 16.1 3.03 120 17.7 99.0 126.0 2.86 566 40.0 47.0 0.37 967 27.0 34.8 2.10 154 338 39.4 37.2 2.45 208 63 48.8 69.0 4.30 75 151 6.4 6.7 2.00 14  Cl Ca Mg mg/L mg/L mg/L 19.0 46 10.7 18.7  116  31.9  104  60.6  43  48.9  116  1,1  229.5  0.7 1.7 0.9  2490.3 265.2 15.2 38.6 788.3 24.8  ND  29.5  21.1  381 85 201 346 95  1.1  82.8  51.0  113  2  21.5 0.6 13.0  0.7  20.6  72.7  56  11  6.0  0.9  8.9  74.1  47  1Q  PZ-902 Mraz well Prairie Creek at E Prairie Creek at A Shaughnessy Springs Cowan Springs Indian Springs  50  23  27  33.2 1.2 0.1 3.3  0.7 1.2  10.8 121.4  81.9 68.1  16 7  7.3 4.8  27.9 1.42 267.0 2.08  0.8  23.5  67.7  54  17.2  19.0  3.72  3.7 0.5 0.0  1.0  62.7  93.4  46  18.2  27.1  2.86  0.9  41.9  262.0  35  7.0  10.6  1.49  9 1 2 34 22 117 34 10  4.7 4.3 1.0  0.9 ND  40.6 5.1  368.7 23.8  84 16  34.5 3.6  37.5 4.2  5.60 1.33  24 106 1  Figure 4-10 Summary diagram of Prairie Flats hydrostratigraphy and groundwater flow patterns during summer months  55  5 DISCUSSION OF SOURCE AND QUANTITY OF GROUNDWATER DISCHARGE This section uses the hydrogeological and geochemical field data presented in Chapter 4 to formulate hypotheses on the origins and quantities of groundwater discharge into the Prairie Flats. The first section presents geochemical arguments that the flats are a discharge zone for locally-recharged groundwaters. Next, groundwater sources and sinks are quantified within the framework of a simplified hydrologic budget and used to estimate the current uranium deposition rate, thereby providing some insight into the history of the deposit.  5.1  SOURCE OF PRAIRIE FLATS GROUNDWATERS  Both the groundwater geochemistry and physical flow regime indicate that the Prairie Flats is a groundwater discharge zone. The accumulation of uranium in this area is an obvious proof, since it was leached from the surrounding rocks rather than coming from atmospheric sources.  In most wetland areas, the decomposition of organic matter  generates carbonic and organic acids which lower groundwater pH. In groundwater-fed areas, these acids are neutralized by mineral bases (eg. CaCOs) dissolved in the discharging groundwaters (Shotyk, 1988). This results in a groundwater p H of between 6 and 8 rather than between 4 and 6, as is the case in the Prairie Flats. The measurement of upward flow gradients across the peat unit, as well as the immediate change in head in the till unit after a rainstorm (Figure 4-9) confirm that groundwaters are discharging here. Infiltrating rainwater could not have caused such a fast response in the underlying till owing to the low hydraulic conductivity of the peat & clay unit.  Further geochemical analysis indicates that the current source of these discharging groundwaters is local rather than regional. A Piper diagram of Prairie Flats groundwaters is given in Figure 5-1, which closely resembles the Piper diagram of other local  56  100  ©  Prairie Flats g r o u n d w a t e r s  A  M r a z well  O  Indian springs Prairie creek  X  C o w a n springs  (J)  S h a u g h n e s s y springs  Figure 5-1 Piper plot of Prairie Flats groundwaters and nearby springs and surface waters. Shaughnessy springs data is taken from Piteau & Associates (1984)  57  groundwaters (Figure 3-4), as opposed to that of regional groundwaters (Figure 3-3). Prairie Flats groundwaters are relatively fresh (TDS < 1000 mg/L), neutral in pH, and are dominated by Ca and HCO3, which is characteristic  of shallowly circulating  groundwaters in glacial sediments (Freeze and Cherry, p. 284). Summertime irrigation activity further up Prairie Valley is a major source of recharge to this shallow groundwater flow system. If regional groundwaters are discharging here, they are diluted beyond recognition by local groundwaters.  Water samples from the Mraz well, Indian Springs, Shaughnessy springs, and Cowan springs add further argument to the local origin of Prairie Flats groundwaters, and have been added to Figure 5-1. Indian springs is the freshest of the springs sampled (TDS < 100 mg/L) and plots close to Prairie Flats groundwaters on the Piper plot. The Mraz well groundwater sample shows a predominance of Na and K , plotting close to samples taken from the deep bedrock wells in the White Lake basin (Figure 3-3). This, plus its low uranium content (< 2 u.g/L) and reducing redox conditions (Eh = -172 mV) suggest that deep bedrock groundwaters have a negligible role in the formation of this deposit.  High uranium concentrations measured at Cowan and Shaughnessy springs (106 u.g/L and 10-24 u.g/L, respectively) suggest that shallow groundwaters carry significant quantities of uranium. Mass balance calculations for a similar surficial uranium deposit in Nevada (Zephyr Cove) proved that modest losses of uranium from near-surface bedrock within the 1.75 k m drainage basin could supply the 40 000 kg of uranium 2  present in the deposit (Otton et. al, 1989). Therefore, the uranium underlying the Prairie Flats need not have migrated great distances, but rather may have come from local bedrock sources within the Summerland basin.  58  5.2  QUANTIFYING GROUNDWATER DISCHARGE  A Simplified Hydrologic Budget A hydrologic budget balances the inflow and outflow components of a given drainage system.  With respect to the groundwater flow regime, it is divided into areas of  recharge (usually at higher elevation) and discharge (usually at lower elevation), as is illustrated in Figure 5-2 (Freeze and Cherry, 1979, pp. 3, 205-207). In a recharge zone, precipitation that is not lost to surface runoff (ie. streamflow) or to evapotranspiration infiltrates the ground surface and recharges the groundwater regime:  R = P - E T  R  - Q  (5.1)  where R is the groundwater recharge (m /yr) P is the precipitation (m /yr) E T is the evapotranspiration rate over the recharge area (m /yr) 3  R  Q is the surface runoff (m /yr) 3  Some of the groundwater recharge goes into storage while the rest moves down the hydraulic gradient as groundwater flow until reaching a groundwater discharge zone. At the groundwater discharge zone, this water may go into storage, evaporate, or exit as stream baseflow (ie. the groundwater contribution to stream flow):  D =B+ET +S D  (5.2)  D  where D is the groundwater discharge (m /yr) B is the baseflow component of streamflow (m /yr) E T D is the groundwater lost to evapotranspiration over the discharge area (m /yr) SD is the groundwater that is stored (m /year) 3  59  Precipitation can be neglected from the above equation because discharge areas usually cover only a small fraction of the total area of the drainage basin, as is illustrated in Figure 5-2.  Similarly, the approximate area of Prairie Valley, from the reservoir to the  base of Giant's head, is 6 km , and that of the Prairie Flats is about 0.75 km . The 2  2  insignificance of precipitation is particularly true during the summer months when high temperatures cause most of the precipitation to evaporate before infiltration.  0<x>o<  Groundwater discharge zone  Figure 5-2 Components of recharge zones and discharge zones in a hydrologic budget, modified from Freeze and Cherry (1979, p. 3)  60  Piteau & Associates (1984) estimated that the groundwater recharge rate for shallow groundwater flow systems in the Trout Creek basin (which they defined as being recharged in local upland areas and having flow depths in the order of 10 to 50m) is approximately 4% of the precipitation. A n average annual shallow groundwater recharge rate of 20.2 mm was calculated by dividing the total annual baseflow to Trout Creek by the catchment area, then multiplying by a factor of 1.3 to account for flow into the creek alluvium and evapotranspiration losses. This is assumed to be equivalent to the average annual shallow groundwater discharge rate in the basin.  Estimating Groundwater Discharge into the Prairie Flats This section estimates the current rate of shallow groundwater discharge into the Prairie Flats by calculating the discharge rate into the peat & clay unit. •  The above approach used by Piteau & Associates (1984) could not be used to calculate the shallow groundwater discharge rate into the Prairie Flats because flow rate data for Prairie Creek during the winter months (when irrigation is effectively turned off) is unavailable.  Instead, the average annual groundwater discharge into the flats was  estimated using Darcy's law, which is written  (5.3)  Q=KiA •7  where Q is the flow rate (m /year) K is the hydraulic conductivity of the porous medium (m/sec) A is the cross-sectional area of flow (m ) 2  As discussed in Chapter 4 and illustrated in Figure 4-10, the horizontal component of the groundwater flux into the flats is negligible compared to the vertical component, therefore total discharge (D) can be approximated using  D = K iv A v  (5.4)  61  where D is the total discharge into the flats (m /year) 3  K is the vertical hydraulic conductivity of the peat and clay units (m/sec) v  i is the vertical hydraulic gradient across the peat and clay units v  A is the approximate surface area of the known extent of the uranium deposit (m )  From Chapter 4, K v is approximately 10" m/s, the average annual iv is approximately 8  0.04, and A is about 0.75 km . Substituting these values into equation 5.4 gives a total 2  discharge rate of 9450 m /year, which is equivalent to about 13 mm/year over the area of 3  the flats. This is of the same order of magnitude as Piteau & Associate's (1984) estimate of 20.2 mm/year for the Trout Creek basin.  During the summer months, most of the groundwater discharge exits the flats as evapotranspiration across the shallow water table (ED). Evapotranspiration can also be calculated from shallow water table fluctuations (White, 1932, Meyboom 1967) but this method requires estimates of groundwater inflow and specific yield, as well as a continuous record of water levels over several days. Using potential evapotranspiration to estimate E would yield too high a value since potential evapotranspiration includes D  water lost from all surfaces (surface water, soil, vegetation, etc.) and not just from the water table.  Relative to evapotranspiration, other sinks such as groundwater storage, So, and baseflow to Prairie Creek, B , are negligible during the summer months. Losses to storage can be ignored since water levels are at an annual low.  Groundwater exiting the flats as  baseflow to Prairie Creek is relatively small since the surface area of the streambed (1700 m excluding culverted sections) represents less than one percent of the total area of the 2  flats (750 000 m ). 2  This statement assumes that the creek bed is also peat ( K  v  approximately 10~ m/s) and that the vertical hydraulic gradient across it is the same as 8  that across the rest of the flats. Also, measurements of hydraulic gradient and seepage rate along Prairie Creek offer no discernable pattern of groundwater discharge or recharge to or from the creek.  62  During the winter months however, evapotranspiration across the water table is effectively turned off due to low temperatures and snow accumulation on ground surface. At this time, most of the discharging groundwater goes into storage, causing a rise in the water table.  For example, in March 1997, average water levels in the shallow  piezometers were 45 cm higher than in August, 1996. Flow rate measurements in Prairie Creek at this time of year, when irrigation is ceased, would yield better estimates of the fraction of groundwater discharge contributing to baseflow of Prairie Creek.  Given the limited amount of field data collected at the Prairie flats, the above calculations provide only rough estimates of the various components of the hydrologic budget. Seasonal climate changes and irrigation practices significantly affect this budget and therefore require further consideration. More stream and groundwater measurements at different locations and at different times of year are needed to improve these estimates.  5.3  N E W INSIGHTS O N DEPOSITIONAL HISTORY  Using the calculated discharge rate and measured uranium concentrations of Prairie Flats groundwaters, the current flux of uranium into Prairie Flats was calculated. Next, this was compared with former estimates to gain new insights on the depositional history.  It has been established by Levinson (1984) that the Prairie flats uranium deposit is less than 10,000 years old, based on U  /Th  isotope ratios. This is a maximum value, as  minimum uranium ages on two soil cores were 2,000 and 8,000 years (Levinson, 1984). Based on a soil survey done by Culbert (1979), there is an estimated 230 000 kg of uranium underlying the flats. Dividing this by the maximum age of the deposit gives a minimum average uranium deposition rate of about 23 kg/yr.  The current flux of uranium into Prairie Flats is calculated using the equation:  C = D[U]  (5.3)  63  where C is flux of uranium into the flats (kg/yr) D is the average groundwater discharge (m /yr) [U] is the average incoming uranium concentration (kg/m )  Uranium concentrations in water samples taken from the deep piezometers, which represent incoming groundwaters, were between 10 and 100 ug/L. This range of uranium concentrations is consistent with other stream and groundwater measurements in the Okanagan area (Boyle, 1982, Ministry of Health B.C., 1981). Uranium concentrations on the order of 1000 ug/L have been measured only at a number of hydrologically isolated, highly alkaline ponds that are subject to high rates of evaporation, which does not represent conditions at the flats. Therefore, 100 ug/L is assumed to be the maximum uranium concentration in groundwaters discharging into the Prairie Flats.  When  multiplied by the estimated discharge rate of 9450 m /yr (section 5.2), the estimated 3  current maximum flux of uranium into the Prairie Flats is 0.95 kg/yr.  C - (9450 m /year)(l x 10" kg/m ) 3  4  3  = 0.95 kg/yr  This calculation assumes that the amount of uranium entering or leaving the flats via Prairie Creek is negligible.  Instead, all of the uranium is assumed to come from  discharging groundwaters, most of which probably evaporates during the summer months and goes into storage during the winter months, according to the budget presented in section 5.2. Ignoring uranium losses or gains from groundwater-surface water interaction is reasonable because:  1. The drainage ditches relaying Prairie Creek across the flats were excavated in the later part of this century whereas the deposit is thousands of years old. Prior to agricultural development the flats were a marshland, 2. There is no visual correlation between the distribution of uranium and the courses of drainage ditches across the flats,  64  3. There is minimal topographic relief around the drainage ditches that would direct groundwater into them. Similarly, the hydraulic head distribution in the vicinity of the drainage ditches does not suggest significant groundwater-surface  water  interaction (see Figure 4-7).  Therefore, the current maximum uranium flux into the flats is at least an order of magnitude lower than the minimum rate estimated by Culbert and Leighton (1988). This discrepancy suggests that a deposit of this size could not have been created i f the uranium deposition rates were the same in the past as they are today. At some point they were much higher, either due to increased groundwater discharge rates and/or higher uranium concentrations in the groundwater. For example, rocks that were scoured, fractured, and granulated by glacial processes were probably stripped of most of their labile uranium thousands of years ago. Furthermore, the Okanagan Valley was probably wetter and cooler at this time and not the dry, hot, semi-desert that it is today. Other past climatic events such as flooding may have enhanced uranium leaching and deposition. This uranium may have initially been deposited in the silt and clay deposits left by the meltwater lake that covered Prairie Valley, and then subsequently redeposited in the overlying peat layers by upward flowing groundwaters.  Another hypothesis is that the Prairie Flats uranium deposit is a more recent phenomenon associated with agricultural development in Prairie Valley. Cultivation and construction activity may have prompted a release of uranium from rocks and sediments, which were flushed down Prairie Valley by increased volumes of groundwater and surface runoff from irrigation and reservoir leakage.  Upon arrival at the Prairie Flats, this water  maintained the marshy conditions and which are responsible for the accumulation of peat. Still, this theory can be disputed by the facts that:  1. Analyses of the U  /Th  activity ratios by Levinson (1984) of peat cores taken  from the flats conclude that the uranium has been in place for at least 2000 years.  65  2. Many other surficial uranium deposits analogous to the Prairie Flats deposit have been discovered in the Okanagan Valley, many in isolated areas unaffected by agricultural development.  Obviously, further research is required to better define the timing of uranium deposition at the Prairie Flats. Further isotopic analysis of the uranium and carbon dating of the peat layer would help to resolve whether they accumulated simultaneously or at different times. A more thorough historical analysis of glacial and anthropogenic activity in this valley would also be useful.  SUMMARY  This chapter provides evidence that Prairie Flats is a discharge zone for groundwaters recharged within or just outside the Summerland basin. These groundwaters encounter carbonate minerals along short, shallow flow paths to the site, which upon discharge neutralize the acids generated in the peat.  Piper plots show that Prairie Flats  groundwaters and other local groundwaters are similar in chemical composition.  A simple hydrologic budget for the Prairie Flats suggests that evapotransporation in the summer, and groundwater storage in the winter are the prevalent sinks for discharging groundwaters. Also, the rate of uranium deposition appears to be much smaller today than in the past. Today's estimated discharge rates are on the order of 9450 m /year, which when combined with maximum incoming uranium concentrations of 100 (J-g/L, give an incoming uranium flux of 0.95 kg/year, at least ten times less than the average rate estimated from the size and age of the deposit.  It is proposed that uranium  deposition rates were higher in the past than they are today, however further scientific investigation is needed to confirm this hypothesis.  66  6 URANIUM GEOCHEMISTRY This chapter summarizes the geochemical behaviour of uranium in natural groundwater environments.  This information is used to form hypotheses on the nature of uranium  fixation at the Prairie Flats based on the geochemical data collected in this thesis. The first section looks at how uranium mobility is enhanced by complexation with soluble species, especially carbonate. The next section presents processes that fix uranium, such as  adsorption,  reductive  precipitation,  evaporative  precipitation,  and microbial  assimilation. Most of this material is taken from laboratory studies with peat.  6.1  MOBILE URANIUM  Uranium exists in two valence states in nature: U  and U . U  exists in primary ore-  forming uranium minerals such as uraninite (UO2), pitchblende (amorphous UO2) and coffinite (USi04).  These minerals are highly insoluble in reducing groundwaters.  However, upon contact with relatively oxygenated waters, U  4 +  is oxidized to U  6 +  and is  2+  released into solution as the uranyl (UO2 ) ion. •  •  2  Uranyl readily forms soluble complexes with oxygen-bearing hgands (eg. CO3 PO4 ", F ", OH"), thereby enhancing its mobility in aqueous environments. 2  2  2  SO4 ",  In the pH  range of most groundwaters (5 to 7), uranyl carbonate complexes are the most important. Uranyl dicarbonate (U02(C03)2 ~) and uranyl bicarbonate (U02(C03)3 ~) form at pH 4.5 2  4  and well into the alkaline range. They are also highly stable, with stability constants on the order of 10 and 10 , respectively (Langmuir, 1997, p. 552). This stability and their 17  22  negative charge make them highly mobile and resistant to processes like adsorption. The stability constant (or formation constant) of a complex is the equilibrium constant, K , for the reaction of a metal with its ligands to form a metal-ligand complex. For example, the equilibrium constant for the reacton U 0 2  2+  + 2CO3 " = U02(C03)2 \ is equal to 2  2  Auo2(co3)2/(Auo2)(Aco3) , where A is the activity of each species. 2  67  Uranyl phosphate complexes are also highly stable in the pH range of 4 to 9, with stability constants on the order of 1 0 1997, p. 5 5 2 ) .  14  for U 0 P 0 " and 1 0 for U 0 H P 0 ° (Langmuir, 8  2  4  2  Phosphate can compete with carbonate  4  only when phosphate  concentrations are significantly enriched (ie. [PC>4 "]T/[C03 "]T > 0.1), which is unlikely in 3  2  most natural groundwater environments (Sandino and Bruno, 1 9 9 2 ) .  6.2  IMMOBILE URANIUM  In the absence of solubilizing ligands such as carbonate, uranyl can be removed from solution by a variety of mechanisms:  (1) adsorption onto organics (2) adsorption onto inorganics (3) reduction and precipitation of uranium minerals (4) precipitation by evaporation (5) microbial assimilation  Each of these is discussed in detail in the following sections:  Adsorption onto Organics Field and laboratory studies present many examples of uranium adsorption onto organicrich soils (Owen and Otton, 1 9 9 5 ; Shotyk, 1 9 8 8 ; Lopatkina, 1 9 6 7 ; Szalay, 1964). The superior ability of organic molecules to adsorb uranium is attributable to their large surface area and abundance of negatively charged ligands which attract UO2 . The term "adsorption" as it is used here implies a chemical complexation or ion exchange process as opposed to purely physical adsorption on molecular surfaces.  Organics are concentrated in the humus layer of soils, which is a dark material made up of decaying plant and animal matter. Humus has three main constituents: humic acid, 68  fulvic acid, and humin, which are defined behaviourally. Fulvic acid is soluble at all pH values, and humin is insoluble at all pH values. Humic acid is insoluble under acid conditions (pH < 6) but soluble under alkaline conditions.  During adsorption to organics, most of the uranium is retained by humic acid (Kochenov et al, 1965; Shanbhag and Choppin, 1981), which exchanges H ions on its carboxyl +  groups for U 0  ions (Szalay, 1964; Borovec et al., 1979; Idiz et. al, 1986). Published  2 + 2  5  8  stability constants for uranyl-humic acid complexes are on the order of 10 to 10° for 1:1 J  complexes and 10 to 10 for 1:2 complexes (Kribek and Podlaha, 1980; Shanbhag and 9  12  Choppin, 1981).  Conditions affecting Uranium Adsorption onto Organics The optimum pH for uranium adsorption onto humic acids is between 4 and 6 (Manskaya et. al 1956; Shanbhag and Choppin, 1981). In sorption experiments with peat, Lopatkina (1967) noted a sharp decline in uranium retention under more alkaline conditions. Between pH 6 and 7.2, nearly all the uranium was adsorbed; at pH 7.8 about one half was adsorbed, and at pH 8.3, none was absorbed. This loss of sorptive capacity at alkaline pH is attributable to the peptization or dissolution of humic acids (Titayeva, 1967) and to the appearance of other complexing agents that scavenge uranium (eg. hydroxides and carbonates).  Uranium sorption rates decrease with increasing dissolved uranium concentration according to a Langmuir adsorption isotherm (Borovec et al, 1979; Szalay, 1964). The Langmuir isotherm relates quantities sorbed to quantities in solution by means of the expression:  X/C = bX /(l + C ) m  (6.1)  b  where X is the quantity adsorbed (mol/g), C is the quantity in solution at equilibrium (mol/L), b is the binding constant of the substrate (L/mol), and X  m  is the sorption  69  capacity of the substrate (mol/g). Another constant, called the geochemical enrichment factor (G.E.F.), is the slope of the tangent of the isotherm (X/C) at low dissolved U concentrations. These are all shown graphically in Figure 6-1.  _  Xm  C (mol/L)  Figure 6-1 A Langmuir Adsorption Isotherm  The sorption capacity for uranium on peat is around 2 to 3 meq/g (Szalay, 1964), and "near perfect" G.E.F.'s on the order of 10,000:1 have been reported (Szalay,1964; Idiz et al.1986; Owen and Otton, 1995).  Adsorption onto Inorganics Studies of uranyl adsorption onto iron oxides and iron oxyhydroxides generally conclude that adsorption increases with increasing pH between pH 5 and 8.5 (Langmuir, 1997, p. 509). Typical G.E.F.'s for uranyl on amorphous Fe oxyhydroxides and on goethite are on the order of 10 and 10 respectively (Langmuir, 1978). Uranyl adsorption onto clays is 6  3  weak at best, for example montmorillonite and kaolinite have G.E.F.'s of only 6 and 2 (Langmuir, 1978).  Where organics and inorganics are present together, organics will dominate and adsorb most of the uranium. Lopatkina (1967) collected both organic and inorganic material from a uranium-rich peatland and found the sorptive capacity of the organic material to  70  be 10 to 100 times greater.  Idiz et al. (1986) found uranium to be preferentially sorbed  to the organic fraction in bog sediments despite the presence of Fe and M n oxides.  Reductive Precipitation of Uranium Minerals The oxidation of organic matter during peat diagenesis creates oxygen depleted, relatively reducing conditions that are ideal for the reductive precipitation of uranium to coffinite (USiCv), uraninite (UO2), or pitchblende (amorphous UO2). can be the organic matter itself, F e C H , or F e 4  2+  2+  Reducing agents  minerals such as pyrite or marcasite, or mobile H S , 2  in the peat (Langmuir, 1997, p.509). This reduction of U  step process. The reaction UC>2 - » U 0 2 2+  +  6 +  to U  4 +  is a two  is instantaneous, whereas UC>2 —» U  slower because a U - 0 bond is broken (DeVoto, 1978, p.4).  +  Therefore  4 +  is  reductive  precipitation of uranyl is a slower and less reversible process than adsorption.  In peatland environments it has been suggested that adsorption and reduction work together in fixing uranium (Andreqev & Chumanchenko, 1964), with adsorption serving to preconcentrate the uranium so that the slower process of reductive precipitation can proceed. Langmuir (1978) states that "whereas mineral solubilities limit only maximal uranyl concentrations, sorption limits uranyl concentrations at all levels below saturation with the least soluble mineral in a given water." Furthermore, i f reduction does not follow adsorption, uranyl can be desorbed by an increase in alkalinity at constant pH, or by raising the pH (Langmuir, 1997, p 509).  Most naturally-occurring uraninites have chemical formulas between  UO2.00  and  UO2.67  (Langmuir, 1997, p. 505) due to partial oxidation after precipitation (Langmuir, 1978). Figure 6-2 shows the stability fields of these species, which overlap with those conditions commonly found in groundwater.  71  2  4  6 pH  8  10  12  Figure 6-2 Eh-pH diagram for aqueous species and solids in the system U-O2-CO2-H2O at 25 °C and 1 bar total pressure. Solid/aqueous boundaries (stippled) are drawn for I U = 1(T M . UDC and UTC are U0 (C0 ) " and U0 (C0 ) " respectively. From Langmuir, 1997, p. 505. 2  2  3  2  4  2  3  3  72  Evaporative Precipitation Uranyl minerals (U ) can be formed by evaporative concentration of uranyl under arid conditions.  Some of the most common minerals, in order of decreasing solubility, are  schoepite ( p - U 0 • 2H 0), autinite (Ca(U0 )2(P0 )2), carnotite ( K ( U 0 ) ( V 0 ) ) and 3  2  2  tyuyamunite (Ca(U02)2(V04)2).  4  2  2  2  4  2  While carnotite and tyuyamunite precipitate at  concentrations in the ppb range, shoepite requires concentrations in the ppm range. Except for carnotite and tyuyanmunite, uranyl minerals are not known to form extensive deposits and are limited to very arid environments and l o w - C 0 waters (Langmuir, 1997, 2  p. 497).  Microbial Assimilation of Uranium Microbes play a part in retaining uranium both directly and indirectly. Indirectly, they help to degrade plant material to produce peat, whose high surface area and high organic content enhance uranium adsorption (Robbins et. al 1990). Humic material also provides carbon energy sources for sulfate-reducing bacteria. These organisms take S 0 dissolved 4  in groundwater and produce H S, which in turn reduces U 2  6 +  to U . 4 +  Further studies have shown that sulfate-reducing bacteria, such as the desulfuricans  Desulfovibrio  can also directly assimilate uranium. Mohagheghi et al. (1984) postulated  that this may be due to the ability of their cell walls to adsorb uranium, thereby holding it in place during the more sluggish reduction step. In work with marine sediments, Lovely et al. (1993) later demonstrated that the Desulfovibrio  desulfuricans  can precipitate  uraninite directly from solution through enzymatic mechanisms.  73  SUMMARY  The speciation and mobility of uranium is controlled by many factors: redox, pH, its concentration in solution, and the presence of complexing ligands. Uranium is most mobile in oxygen-rich, carbonate-bearing waters having a neutral to alkaline pH. Under these conditions uranium as U  6 +  is complexed with carbonate as uranyl dicarbonate and  uranyl tricarbonate.  In peaty groundwater environments where the pH is between 4 and 6, organics can overcome the uranyl-carbonate complex and thereby immobilize uranium.  Organics  contribute humic acid ligands which adsorb uranium, and also foster a reducing environment. Uranyl which is initially adsorbed may be reduced to a U  4 +  mineral by the  organic matter itself, by mobile reductants such as H2S and F e , or by microbial 2+  processes.  Other fixation mechanisms such as adsorption on colloidal iron oxides and  evaporite formation are less important in such organic-rich environments.  74  7 DISCUSSION OF GROUNDWATER GEOCHEMISTRY AND URANIUM DEPOSITIONAL CONTROLS This chapter begins with a discussion of key processes which control the geochemistry of Prairie Flats groundwaters, based on the aqueous geochemical data collected on site. These controls include carbonate mineral-solution equilibria, redox reactions, and evapotranspiration.  Next, their role in uranium retention via adsorption, reductive  precipitation, and evaporative precipitation is addressed. Redox conditions may further explain uranium concentration near ground surface, and the hydrogeology offers new insights on its lateral distribution. Finally, the issue of uranium remobilization due to cultural inputs is addressed using two field observations.  7.1  CONTROLS ON T H E GROUNDWATER GEOCHEMISTRY  Carbonate Mineral-Solution Equilibria Chemical analyses of Prairie Flats groundwaters show that carbonate mineral-solution equilibria exert  a major  control on the  groundwater  chemistry.  Bicarbonate  concentrations are relatively high, averaging between 500 and 650 mg/L. This maintains the pH in a neutral range of 6.0 to 7.8, and accounts for the high acid-neutralizing capacity. Bicarbonate likely comes from the incorporation of soil C 0 in groundwater 2  recharge zones and the subsequent dissolution of carbonate minerals during groundwater migration to the site.  Of these carbonate minerals, calcite is a likely source of H C 0 . In a plot of H C 0 vs pH 3  3  (see Figure 7-1), most samples plot on or slightly above the saturation line with respect to calcite.  Similarly, PHREEQC calculations give saturation indices with respect to  calcite of between -0.5 and +0.5 (see Appendix F). The linear relationship between H C 0 and Ca also suggests that calcite is the main source of H C 0 (see Figure 7-2). 3  3  75  Calcite is preferred over dolomite since the Ca to M g ratios are consistently 2 or more, whereas in equilibrium with dolomite they would be closer to unity (Langmuir, 1997, p. 209).  Figure 7-1 Plot of Prairie Flats groundwaters in relation to waters undergoing open system dissolution of calcite with Pc02 °f lO'l-^ atm  Figure 7-2 Relationship between H C O 3 and Ca in Prairie Flats groundwaters, suggesting calcite as a source of H C O 3  76  Redox Conditions  Three methods were used to interpret redox conditions in Prairie Flats groundwaters. First, Eh was measured directly using a redox probe. Second, concentrations of redoxsensitive species, namely N 0 , S 0 and Fe 3  4  2+  were compared with the redox reaction  sequences introduced by Froelich et al. (1978).  Lastly, measured pH's and Fe  2+  concentrations were used to calculate Eh assuming equilibrium with ferrihydrite, Fe(OH) , also referred to as hydrous ferric oxide, F e ( O H ) n H 0 or H F O . 3  3  2  Eh measurements Meaningful measurements of Eh using platinum electrode redox probes have been obtained for groundwaters enriched in iron (Langmuir, 1997 p. 411). However, since no correlation was found between measured Eh and dissolved Fe  2+  concentrations (see  Figure 7-3) the Eh measurements were discarded due to suspected fouling of the platinum electrode by organics.  100  200  300  400  Measured Eh (mV) Figure 7-3 Inconsistency between measured Eh and Fe2+ concentrations in Prairie Flats groundwaters  77  Redox Reaction Sequences Froelich et al. (1978) indroduced a series of redox reactions that occur during organic diagenesis in upper marine sediments. Assuming that organic matter is the principal electron-donor, 0 , followed by N 0 , F e 0 or FeOOH, M n 0 , S 0 are subsequently 2  3  2  3  reduced in order of greatest free energy yield.  2  4  Thus under equilibrium conditions, each  redox-sensitive species is expected to exist only within a particular Eh range. The same applies to microorganisms that act as mediators or catalysts in these redox reactions, including aerobic heterotrophs, denitrifiers, sulfate-reducers, and so on (Sturnm and Morgan, 1981, p. 458).  Based on Froelich's redox reaction series, Berner (1981, p. 359-365) illustrates how concentrations of various redox-sensitive species are expected to vary within the sediment column, from oxic conditions near the surface to anoxic conditions at depth. Stumm and Morgan (1981, p. 460) assigned an Eh scale to the same sequence of reactions (see Figure 7-4a). Thus, by replacing Berner's "depth" axis with Stumm and Morgan's "Eh" scale, we can assign an approximate Eh range to a water sample based on the subset of redox-sensitive species present. This is illustrated in Figure 7-4b.  At the Prairie Flats, N 0 , S 0 , and Fe concentrations were measured at each of the four 2+  3  4  nested piezometer pairs, and plotted versus depth in Figure 7-5. S 0 , N 0 , and Fe were 2+  4  3  found together just below the water table, which from Figure 7-4 between roughly -150 and +50 mV. About one metre below that, Fe  is possible only 2+  has disappeared  and only S 0 and N 0 remain, corresponding to an Eh of between +50 and 250 mV. 4  3  Dilution by incoming groundwaters can not account for the disappearance of iron in the till unit because measurable amounts of S 0 and N 0 are still present. 4  3  78  "p.  cu  m •a oo  0\  u  t>J0  HZ)  ^  M  c .S ctf  HZ)  3  «  e  s  O  W  S g £ >, £ ^  0* I  *  M> V <H  O j to o»  (0  (0 m  [ OO Kl  co  Kl O  O  —  CM  •  t-iisis Si S  E  £I  I •? 51 ii<H 5  °s  "7^  CU  fr  £  t  ii  o  in  i  1_  i  -j-i te ^  s  • i i  5 g t; .2  CX  co  o  £ o o u  s a ° .5  o  3 t  "J3  O ^3 cu "H  I  s  s  a •£ to co CU O  -2  •3• TJ  111  51 i l l  5  u  1  o S  •  in c o  -si  /2s  2 c  <5  cu  i3 * c S  a*o  51  T3  -o •o  CU  a* (U  oo ^  HZS  -  g o  <* feb  T.5 ^ CU  >-  3  t>0  CO  cu  bo  C  eg  s -§ 71  300  a)  250 Ui E  200  c  o  o N03  fl 150 ™  • Fe 2+  8 100  A  S04  c  o o  50  r-£H  |-5H  0 0.5  1  1.5  2  2.5  Depth (m)  16 b) 14  j  12  £10  o N03  c  o  • Fe2+  re 8 a) o c  o O  c  D  4 2 0 0.5  1  1.5  2  2.5  Depth (m)  Figure 7-5 Change in concentration of a) N O 3 , Fe?+, and SO4 with depth in Prairie Flats groundwaters and b) showing N O 3 and Fe^ " only at a smaller vertical scale 4  80  Still, it is unusual to find N 0 and Fe together since N 0 is a strong oxidant capable of 2+  3  3  oxidizing Fe according to the reaction: 2+  10Fe + 2N0 " + 14H 0 - » lOFeOOH + N + 18H 2+  3  2  +  2  (7.1)  One explanation for the persistence of N 0 is its introduction at ground surface in 3  amounts that exceed the immediate reducing capacity of the soil. Nitrate concentrations in the 10's of mg/L (as measured at piezometers 202, 502, 602 and 902) are much greater than the 0.1 mg/1 typical of peat environments (Shotyk, 1988), pointing to fertilizer and manure application up Prairie Valley as likely sources of N 0 . 3  The formation of  intermediates like N0 ", N O and N 0 before N may also slow this reduction step. At the 2  same time, Fe  2+  2  2  may be continuously regenerated by the dissolution of Fe  2+  bearing  minerals and the reduction of Fe-oxyhydroxides in the soil.  Fe / Ferrihydrite Equilibria 2+  Langmuir (1997, p. 436) states that in aquatic environments, the oxidation of Fe  2+  forms  ferric hydroxides and oxyhydroxides that are mixtures of amorphous material and goethite (a-FeOOH). The amorphous material is most commonly ferrihydrite - Fe(OH) . 3  Iron-oxidizing bacteria can catalyze its production (Steinmann and Shotyk, 1997), and soil horizons rich in carbon will preferentially precipitate ferrihydrite before goethite, hematite, and lepidocrocite (Schwertmann et al, 1986). Colloidal mixtures of ferrihydrite and organic matter produce a thin "oily" film on puddles of water at ground surface, as was observed by the author at the Prairie Flats (see Plate 7-1).  81  10cm  Therefore ferrihydrite is believed to be in equilibrium with the redox environment of the Prairie Flats peats according to the reaction:  Fe + 2 H 0 <-» Fe(OH) 2+  2  (7.2)  + 3H + e +  3(am)  Eh was calculated using:  Eh = E - 0.18pH - 0.0592 log a(Fe ) E = 0.975 V to 1.06 V 2+  0  0  (7.3)  (Langmuir, 1997, pp. 417,437,482)  The value of E is difficult to pinpoint since it depends on the mass ratio of ferrihydrite to 0  other species in the precipitate (eg. goethite), and on its degree of crystallization. The solubility product (K^) of a mineral is defined as the equilibrium constant for the dissolution of a solid salt to give its ions in solution. For the reaction M N m  n ( s )  omlVT* +  nN" ", the K is equal to ( A ) ( A ) where A is the activity of each species. The K for 1  m  s p  n  M  N  s p  ferrihydrite can range from 10" to 10" M (moles/L), higher for fresh precipitates and 37  lower for more crystallized forms. A K K  39  of 10" corresponds to an E of 1.06, and for a 37  s p  0  of IO , E is 0.975 (Langmuir, 1997, pp. 417, 437, 482). Since the water chemistry -39  s p  0  data offer no means of distinguishing the crystal nature of ferrihydrite, a mean E of 1.02 0  was used to calculate the Eh value of each water sample. These were assigned error bars of ± 0.0425 V , which correspond to the maximum and minimum E values. For those 0  showing no detectable iron, the default Eh of 240 mV was chosen, which corresponds to relatively oxygenated waters (Langmuir, 1997, p. 411).  Calculated Eh values are plotted against depth in Figure 7-6. Redox potentials in the peat unit just below the water table generally lie between -100 and 200 mV, and in the till unit underlying the peat, between 50 and 250 mV. This gives an overall Eh range of between -100 mV and 250 mV, which is the same as that inferred from the redox reaction sequences.  83  Figure 7-6 Depth-Eh plot showing relatively reducing groundwater conditions i n the peat unit and relatively oxidizing conditions in the underlying t i l l unit  In summary, both methods of determining redox conditions indicate conditions to be relatively more reducing near surface than at depth. This is due to a combination of two processes: 1) the intense degradation of organic matter in the root zone and just below the water table which sets up an oxygen-poor and relatively reducing environment, and 2) the upward discharge of relatively oxygenated groundwaters from the till unit below the peat. Consistent upward hydraulic gradients measured at the piezometers pairs together with the consistent absence of dissolved Fe  2+  in the deep piezometers confirm the second  process. Eh values probably continue to decrease with depth below the water table, but at some point undergo a redox reversal upon encountering the zone of influence of the discharging groundwaters.  84  Evapotranspiration At all piezometer pairs the major ions were found to generally increase in concentration towards ground surface, as shown in Figure 7-7 and Figure 7-8.  If evapotranspiration were the cause of this upward increase in ion concentration, these ions would show the same rate of enrichment as C l , which is a conservative species. Thus, a histogram of the percent increase in concentration over depth, normalized by that of C l was constructed (see Figure 7-9). From the histogram, it is evident that H C 0 , Ca, 3  M g , Na, and K are being concentrated by evapotranspiration. The fact that H C 0 shows 3  the same pattern as these other ions further proves that the Prairie Flats is a groundwater discharge zone. If this were a recharge zone, a more dramatic increase in H C 0 would be 3  expected due to the incorporation of C 0 during infiltration, which is produced by root 2  respiration and the oxidation of organics.  Upward increases in N 0 S 0 , and Fe concentrations are more acute than for C l due to 2+  3  4  the change in redox environment between the peat unit and the underlying till. Also, the fact that N 0 and Fe underwent more drastic declines than S 0 indicates that conditions 2+  3  4  just below the water table are not more reducing than -250 mV, since this is the point at which S 0 reduction is initiated. Greater changes in dissolved U concentration relative to 4  C l may be explained by desorption or redissolution of fixed U .  85  0  2  4  6  0  500  Fe2+ (mg/L) Figure 7-7 Increase in concentrations of Ca, Na, Mg,K, FejOT, in Prairie Flats groundwaters  1000  Uranium (ug/L) a n  d U towards ground surface  86  -I 0  ,  500  S 0 4 (mg/L)  1  1000  3 _| 0  ,  200  ,  400  ,  600  1  800  H C 0 3 (mg/L)  Figure 7-8 Increase in concentrations of Cl, NO3, SO4, and HCO3 towards ground surface in Prairie Flats groundwaters  87  CM O CD  CM O  CM O  CM O  ^—  T  T—  T—  o  o  O  o  5;  06  • • •  CD  — 1  1  frOS CON €00H O -vf  LO CO  O CO  LO CM  O CM  LO  <r-  eoejjns punojB S P J B M O J [10] u| 9 S B 9 J O U | /aoejjns punojB S P J B M O J M ui esBajoui  7.2  CONTROLS ON URANIUM FIXATION  Adsorption Uranium adsorption refers to the replacement of a proton or other cation on a humic acid molecule by the uranyl ion, U 0 permanent, as U 0  2 + 2  2 + 2  .  Adsorption is a relatively quick process, but not  can later be desorbed from the exchange site. Indirect proof that U  adsorption is occurring at this site comes from evidence of desorption. In brief, these observations include:  1) the detection of dissolved U in amounts that exceed background levels of « 10 p.g/L 2) increasing U levels in Prairie Creek from entrance (E) to exit (A) of the flats 3) a relationship between dissolved U and carbonate 4) a relationship between dissolved U and soil U  A review of the data presented in Table 4-5 confirms the first two points. The last three require further explanation.  Dissolved U and Bicarbonate Figure 7-10a plots dissolved U against H C 0 concentrations measured in Prairie Flats 3  groundwaters in March '98.  Here, U values exponentially increase with increasing  H C 0 , due to the high stability and poor sorptive behaviour of uranyl-carbonate 3  complexes. Points falling below the line of best fit may denote U that is not easily desorbed, but rather is precipitated as a U mineral. 4+  The same relationship between dissolved U and H C 0 can be inferred from the Sept '97 3  data set where U and conductivity were measured but H C 0  3  was not.  Figure 7-1 Ob  shows an exponential relationship between dissolved U and conductivity. The March 89  data set shows that conductivity is proportional to the H C 0 concentration (see Figure 3  10c).  E .o  1500  CO  1000  > U 3  •o c o  100  10000  Uranium (ug/L)  o  2000 1500 1000 . . . . IWL •  500 0 0  500  1000  HC03 (mg/L)  7-10 Evidence of uranyl-bicarbonate complexation in Prairie Flats groundwaters in (a) March '98 and (b) September '97, inferred from the (c) HCO3- conductivity relationship found in March data Figure  90  The relationship between dissolved U and H C 0 is even more apparent when all of the 3  field data, including spring and surface waters, is plotted (see Figure 7-11). Furthermore, speciation calculations using PHREEQC and a corrected version of the WATEQ4F database found dissolved U to be complexed predominantly with carbonate (see Appendix F). In these calculations, uranyl-carbonate complexes outnumbered uranylphosphate complexes by more than 10 . 3  Therefore other complexing agents, such as  phosphate and dissolved organic acids play a negligible role in enhancing U mobility.  1000 -,  100 23  10  TJ d) > O  (0 11  s  0.1 200  400  600  800  HC03 (mg/L)  1000 •- •I  100  I ~ I  _ _J  i—*—i  W §>  1000  i—T==i*-  I ^I  10  100  I— ^  10  0.1 0  20  40  Dissolved Organic Carbon (mg/L)  60  0.5  1  1.5  P04 (mg/L)  Figure 7-11 Evidence of uranyl-carbonate complexation in waters in and around Prairie Flats as opposed to complexation with dissolved organic carbon and PO4 91  Dissolved U and Soil U Carbonate-bearing groundwaters in prolonged contact with high concentrations of sorbed U would be expected to have higher concentrations of dissolved U . In this study, soil U concentrations nearest to the shallow piezometers were obtained from a soil survey of Prairie Flats by Culbert (1979). Graphs of soil U versus dissolved U for the Sept '97 and March '98 sampling rounds is shown in Figure 7-12. adsorption isotherms.  Both resemble Langmuir  The same relationship between soil U and dissolved U is  illustrated in the contour plots of dissolved U concentrations shown in Figures 7-13 and 7-14. Dissolved U concentrations in Sept '97 are higher than those measured in March '98, probably because evapotranspiration rates are higher in late summer than in early spring. Like the soil U distribution shown in Figure l-3a, the highest values are centred around piezometers 401/402 (Figure 4-1).  92  Figure 7-12 Relationship between soil U and dissolved U in a) March '98 and b) September '97  93  Contours in ug/L uranium  Figure 7-13 Contour plot of dissolved U concentrations (in u.g/L) in shallow groundwaters underlying Prairie Flats (March '98)  94  Contours in ug/L uranium  Figure 7-14 Contour plot of dissolved U concentrations (in u.g/L) in shallow groundwaters underlying Prairie Flats (September '97)  95  Reductive Precipitation  In a reducing environment, the U 0 (U0  2(c)  ), pitchblende ( U 0  adsorption.  2(am)  2 + 2  ion may be reduced to U  4 +  and precipitate uraninite  ) , or coffmite (USI0 ), which usually occurs after 4  Reoxidation and resolubilization of U is impeded by sluggish reaction  kinetics, thereby making it a more "permanent" means of U fixation than adsorption. Evidence of reductive precipitation can be drawn from an analysis of the aqueous geochemical data collected on site. These interpretive methods are detailed below.  Eh-pH Diagrams Figure 7-15 is an Eh-pH diagram for the U - 0 - C 0 - H 0 system that incorporates the 2  2  2  mean dissolved U and H C 0 concentrations of Prairie Flats groundwaters. Measured pH 3  and calculated Eh values are used to plot the shallow groundwater data. Compared to Figure 6-2, the U 0 ( C 0 ) 7 U 0 2  2  3  2  3  and U 0 ( C 0 ) 7 U 0 2  7(C)  2  3  2  4  9(c)  solid-solution boundaries  have been shifted down because of the higher levels of dissolved carbonate. This causes a collapse of the U 0 3  8(c)  , U 0 4  9 ( c )  and U 0 3  7 ( c )  stability zones, thereby making U 0  2 ( c )  the more  likely initial precipitate.  Uranium mineral saturation indices The aqueous geochemical data for each water sample were entered into PHREEQC to give saturation indices with respect to various U mineral phases (see Appendix F). Saturation indices for U 0  2 ( c )  and U 0 4  9 ( c )  are plotted in Figure 7-16. In all of the shallow  piezometers, conditions are saturated to supersaturated with respect to U 0 indices with respect to U 0 4  is being oxidized to U 0 4  9 ( c )  9 ( c )  2(c)  . Saturation  are even greater, although whether or not this means U 0 is unclear.  Saturation indices for U 0 3  consistently negative, and the database does not include U 0 3  probable initial U precipitate is uraninite, U 0  2(c)  7(c)  8 ( c )  and U 0  2 ( a m )  2 ( c )  are  . Therefore, the most  .  96  300 250  U02(C03) U 0 2 ( C 0 3 ) 2 2U 0 2 ( C 0 3 ) 3 4-  200 i/"k i rfl  150 100  >  50  .c  U 0 2 (c)  LU  0 -50 U409-UDC  -100  U307-UDC  -150 8  6  -200  PH Figure 7-15 Eh-pH diagram showing calculated Eh and measured pH in shallow groundwaters underlying Prairie Flats. Also shown are the collapsed stability zones of U30s( ), U409( ) and U3C>7( ) c  C  C  97  Effect of Eh, pH, and H C Q on U-Mineral Saturation 3  To test the sensitivity of reductive precipitation of U to changes in Eh, pH, H C 0 and 3  dissolved U , a number of PHREEQC simulations were made on a representative sample. PZ-702 was chosen for its proximity to the average U , H C 0 and redox conditions, and 3  for its good charge balance. A l l relevant analytical data was entered into the model and pH, Eh, and H C 0 indices of U 0  2 ( c )  3  values were varied in turn across the measured range.  and U 0 4  Saturation  are plotted against these parameters in Figure 7-17.  9 ( c )  Of these variables, Eh and pH were found to have the greatest effect on TJ mineral saturation. U 0  2 ( c )  and U 0 4  saturation is reached at Eh's below 0 mV and at pH's below  9 ( c )  7.5. Under these conditions, U can precipitate at dissolved U concentrations as low as 100 ug/L. Increasing H C 0 concentrations slightly increases U 0 3  2 ( c )  solubility because of  its tendency to form soluble uranyl-carbonate complexes.  Although PHREEQC simulations gave both U 0 phases, these plots show U 0 become more reducing, U 0  2 ( c )  2 ( c )  2 ( c )  and U 0 4  to be the more stable U-precipitate. As Eh conditions  is the first to attain supersaturation. At Eh values around 0  mV, it maintains supersaturation up to pH's of 7.7, whereas U 0 4  up to p H 7.5.  U0  2 ( c )  as possible U mineral  9 ( c )  saturation is less affected by H C 0 , precipitating at H C 0 3  concentrations up to 700 mg/L, unlike the 500 mg/L for U 0 4  precipitate at much lower dissolved U concentrations than U 0 4  partial oxidation to U 0 4  can precipitate only  9 ( c )  9 ( c )  9(c)  9(c)  .  .  Also, U 0  2 ( c )  3  can  Therefore, although  is possible, it is likely to be limited.  98  s  o  -"'  OJ  O ZD  o CT>  O ZD  CO CO  •  a1 as a° as 2I  c o cd o o  s E  o N  o  EL to  'cl o u ca o oC  O  °- m  D  as ai as as N  ID  U CN  O  o  o  4—1  o p, co OJ  co U  o +3  a  ,g c o  ca  N  Q_  ca CO (4-1  O  o OH ^3  01 ha  3  xepu| uo;;ejn;es  ex  o o T3  C  Cd •  'o" CN  o  o  o CU  a  % CO  CU  o  •3 c c o  3ed co C O  c o  '-*c-d » J3  oo co ^j- eg o  CN  C  TJ-  xapu| uoueimes  xepu| uojiejmes  cu o c o  o  *0 CU  O  CN  O D  §> II  j> "o  _l  CO CO  3  O I" O m JL, II  a.  o  i i  LO  E  S  Lu £ K 2.  /  •3 XI  C  cd  C .o cd  U C cu  o x— / o /  C o  o  CO  O U ffi  > JC  a  ''O  1  /  rd  W  o c£  1  C+H  W  CN  xepu| uojiBjnies  xapu| uorjBjrnBS  t--  cu 0£  Evidence from Analogous Deposits Studies of other surficial U deposits in the United States have made similar conclusions about the role of adsorption and reductive precipitation in retaining U . Those that have been researched in sufficient detail are the Zephyr Cove deposit on Lake Tahoe (Zielinski et al. 1988, Johnson et al. 1987, Otton et al. 1989, Owen and Otton 1995) and the Flodelle Creek deposit in Stevens County, N E Washington state (Zielinski et al. 1987, Zielinski et al. 1986, Owen and Otton 1995). Geochemical studies carried out at these two sites concluded that where Eh conditions in the saturated zone are above 100 mV, U is fixed by adsorption onto humic acids.  Where Eh conditions fall below 100 mV,  reductive precipitation occurs, either directly or indirectly following the burial and diagenesis of previously sorbed U .  Evaporative Precipitation Geochemical analyses of shallow and deep groundwaters suggest that evaporation may increase concentrations of dissolved constituents near the top of the peat unit, but not enough to precipitate U indices for U  6 +  6 +  minerals.  According to PHREEQC simulations, saturation  minerals such as schoepite ((3-U0 • 2H 0), autinite (Ca(U0 ) (P0 ) ), 3  2  2  2  4  2  carnotite ( K ( U 0 ) ( V 0 ) and tyuyamunite (Ca(U0 (V0 ) ) never rise above 0. This can 2  2  4  2  2  4  2  be explained by the high solubilities of these minerals (most in the mg/L range), which are enhanced in high-C0 waters due to uranyl-carbonate complexation (Langmuir, 1997, 2  p. 497).  101  7.3  C O N T R O L S O N SOIL U R A N I U M D I S T R I B U T I O N  Vertical Distribution One puzzling feature of the Prairie Flats deposit is that the highest concentrations of U exist near the ground surface (see Figure l-3b), which is also true of the organic content of the soils. Simply put, the organics provide 1) sorption sites for the initial entrapment of U , and 2) prolonged reducing conditions which facilitate U reduction and precipitation to a more stable mineral form.  PHREEQC simulations show that the near-surface environment is more conducive to the reductive precipitation of U 0  2 ( c )  than at depth, as shown in Figure 7-18. This is very  different from most groundwater environments where conditions tend to be more reducing at depth. The cause of this redox reversal is the relatively oxygenated waters discharging into the peat unit from the till unit below.  This causes a continual  refreshening of the base of the peat and resolubilization of U .  0 0.5  3 -12  -8  -4  0  4  SIU02 Figure 7-18 Change i n U 0 2 ( c ) saturation indices w i t h depth i n Prairie Flats groundwaters, showing more favourable precipitation conditions near-surface  102  Lateral Distribution With respect to the lateral distribution of U across Prairie Flats, it is interesting to note that the highest soil U concentrations (see Figure l-3a) coincide with the thickest sections of peat and clay (see Figure 4-3). Both sets of contours encircle PZ 401/402 and have the same kidney-bean shape. Groundwaters discharging from the till into the peat would be slowed down due to the drop in hydraulic conductivity, thereby allowing ample time for U entrapment. Therefore, the peat unit not only provides the right chemistry for U retention, but also the right flow conditions.  7.4  POTENTIAL FOR URANIUM REMOBILIZATION  Evidence presented in this study and by Culbert (1980) shows that a major proportion of the U in the Prairie Flats is loosely retained by adsorption and could potentially be remobilized under the right geochemical conditions.  Once remobilized, substantial  amounts may enter Prairie Creek and discharge into Lake Okanagan.  In 1981, concentrations of dissolved U in runoff waters around Summerland were at or above the drinking water standard of 20 L t g / L , possibly indicating some leaching of the Prairie Flats deposit (Ministry of Health, 1981). In the water sampling results from this thesis, two other possible pieces of evidence are highlighted, namely;  • remarkably anomalous dissolved U concentrations near a household septic field and • a marked increase in dissolved U concentration at the exit point of Prairie Creek during sewer system installation in the spring of 1998.  Septic Field Piezometer 202 is located adjacent to a private septic field, as shown in Figure 7-19. The groundwater sampled here in March '98 had a highly anomalous U value (nearly 4  103  mg/L!), as shown in Table 7-1. It also had high concentrations of N 0 , P 0 , S 0 , and 3  4  4  dissolved organic carbon, as is characteristic of septic effluent.  Table 7-1: Water chemistry data at piezometer 202 (March, 1998) pH  Eh(V) H C 0  3  OC NO  PO4  SO4  Cl  Ca  Mg  Na  K  3  7.1  -0.027  1354  157  43  U (Mg/L)  0.7  2490  26  38  5  1  166  630  10.1  3961  * all values given in mg/L unless otherwise specified  This water sample exemplifies the potential impact of increased septic discharge into the flats. The assortment of oxidants and complexing agents in septic effluent comprises a strong leach solution that can strip U from the peat.  Such high dissolved U  concentrations at PZ 202 may be due to reoxidation and desorption working together, wherein N 0 reoxidizes U 0 to uranyl, which then complexes with H C 0 . Therefore, it 3  2  3  is important that mixing of septic effluent with Prairie Flats groundwaters be avoided as residential development continues to expand up Prairie Valley and Cartwright mountain.  Prairie Creek Results from both water sampling rounds show an increase in TJ levels in Prairie Creek from entrance (E) to exit (A) of the flats. In Sept '97 concentrations increased from 22 to 34 ug/L, and in March '98 from 34 to 117 ug/L. One explanation for the dramatic increase in March is the expansion of the town's sewer system on the property adjacent to the flats (see Figure 7-19).  Here the excavation of ditches exposed the peat unit to air  and thereby caused oxidization of any U 0 present to uranyl, which subsequently be 2  washed into Prairie Creek by rainstorms and other runoff events. Further field evidence is required at this point to support this claim.  104  105  Since the Prairie Flats U deposit has been around for thousands of years already, the likelihood of natural remobilization is low. However, this section points out that cultural disturbances such as fertilizer use, septic discharge, and soil excavation may increase the chances of such an event.  SUMMARY This chapter identified the principle controls of the aqueous geochemistry of Prairie Flats groundwaters and demonstrated how these affect the mobility and distribution of U in the deposit. Calcite mineral-solution equilibria control the H C 0 concentrations, and redox 3  conditions reflect the near-surface oxidation of organic material and the discharge of fresh, relatively oxygenated groundwaters from below. After initial adsorption to the organics, U is likely reduced and precipitated as U 0  2 ( c )  in this near-surface reducing  environment. This is shown by Eh-pH diagrams and PHREEQC simulations, which also indicate that evaporative  precipitation of U to some U  6 +  mineral is unlikely.  Remobilization of U and subsequent migration to Okanagan Lake via Prairie Creek could occur i f agricultural and residential activities continue to intensify in the area, but at this point the likelihood and consequences of such an event require further consideration.  106  8 GENERAL CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK This thesis further characterizes the hydrogeology and aqueous geochemistry of the Prairie Flats uranium deposit. When combined with the soil and surface water studies already carried out at the site and surrounding region, it offers a more complete picture of how the uranium deposit was formed, and in what condition it lies today.  The hydrostratigraphy of the Prairie Flats consists of a low-permeability peat & clay unit overlying a well-drained till unit, followed by bedrock. The peat unit is up to 3m thick and has a horizontal hydraulic conductivity of about 6 x 10~ m/s and vertical hydraulic 7  conductivity on the order of 10" m/s. The till unit has a horizontal hydraulic conductivity 8  on the order of 1C/ m/s. Further augering or drilling in and around the Prairie Flats 5  would help to delineate the vertical and lateral extent of the till unit and clarify its connection to the neighbouring glaciofluvial deposits.  Physical measurements of the groundwater flow regime underlying the flats indicate that groundwaters are discharging upwards into the uranium-bearing peat unit. During the summer months, the average lateral component of the hydraulic gradient through the peat is around 0.005, whereas the average vertical component is 0.04.  This upward flow is  topographically-driven since the flats lie at the base of Prairie Valley and are encircled by a number of small mountains.  More frequent monitoring of water levels in the  piezometers is needed to identify the effects of seasonal climate change, rainfall events, irrigation practices and other anthropogenic activities on the groundwater hydrology of the flats.  Darcy's law and a simplified hydrologic budget were used to estimate the amount of groundwater entering and leaving the Prairie Flats, respectively. The estimated total groundwater discharge rate is approximately 9450 mVyear, assuming that recharge to  107  groundwater from Prairie Creek are negligible. During the summer months, most of this discharge leaves the flats as evapotranspiration, and in the winter months, it goes into storage. The limited understanding of groundwater-surface water interaction across the flats could be improved with more frequent and extensive measurements of seepage rates, hydraulic gradients, flow rates, etc. in the drainage ditches (Prairie Creek).  Combining the groundwater discharge rate and average incoming uranium concentration gives an estimate of the current uranium deposition rate in the Prairie Flats of 0.95 kg U per year.  This is an order of magnitude smaller than the rate predicted given the  estimated tonnage and age of the uranium in the deposit. Therefore, uranium deposition rates were probably much higher in the past than they are today. For example, recentlyscoured and fractured rock surfaces combined with a colder, wetter climate would result in greater amounts of uranium leaching and deposition in this region. Radiocarbon dating of the peat unit that hosts the uranium may help to determine whether the deposit was formed thousands of years ago closer to deglaciation, or more recently as a consequence of the onset of agricultural activity in Prairie Valley.  Water chemistry analyses confirm that the Prairie Flats are a groundwater discharge zone for local groundwaters.  B y definition, local groundwaters are recharged within the  Summerland basin and flow at shallow depths (< 100m) in the glacial sediments and shallow bedrock. Incoming groundwaters are relatively oxidizing and have a C a - H C 0  3  chemistry with low concentrations of dissolved solids. Piper plots show a match between Prairie Flats groundwaters and other local groundwaters.  If regional waters are  discharging at this site, dilution by local groundwaters prevents their identification. Assuming that groundwater flow paths are unchanging in time, it is possible that the uranium in this deposit came from source rocks within the Summerland basin.  The chemistry of Prairie Flats groundwaters indicates that both adsorption and reductive precipitation of uranium are dominant fixation mechanisms.  Both processes are  facilitated by organic molecules in the peat which initially adsorb uranium and whose  108  oxidation sets up an oxygen-poor, relatively reducing environment suited to the reductive precipitation of U 0  2(c)  . Given the high competition for uranium by dissolved carbonate  species, neither process is fully efficient, therefore both were probably involved in generating a deposit of this size. Further analysis of soil samples by Scanning Electron Microscopy or X-ray diffraction may better identify uranium minerals on these soils and distinguish how much of the uranium is adsorbed, and how much is reduced.  Combining the geochemistry with the groundwater flow patterns offers further insight on the uranium distribution. Highest uranium values correspond to a thick mound of peat in the centre of the deposit because flow is slowed and contact times are higher. Uranium is concentrated near the ground surface because of the reducing environment set up by the oxidation of organics.  Retention at depth is inhibited by the discharge of relatively  oxygenating waters from the till unit below the peat. A multilevel piezometer would offer more resolution on the changes in groundwater chemistry with depth.  Anomalous counts of uranium near a septic field and during excavation of the peat warn that remobilization could occur i f the chemistry of Prairie Flats groundwaters is significantly changed. Leaching tests with septic discharge, for example, would help to further evaluate this risk.  109  REFERENCES Ahonen, Lasse, Ervanne, H . , Ruskeeniemi, Jaakkola, T., and Blomqvist, R. 1993. Uranium mineral - groundwater equilibration at the Palmottu natural analogue study site, Finland. Mat. Res. Soc. symp. proc. 294: 497-504.. Andreyev and Chumachenko. 1964. Reduction of uranium by natural organic substances. Geochemistry International. 1: 3-7. Appelo, C.A.J, and Postma, D. 1994. Geochemistry, Groundwater and Pollution. Rotterdam, The Netherlands: A . A . Balkema. ISBN 90 5410 105 9 Bates, David V., Murray, James W., and Raudsepp, Valter. 1980. Royal Commission of Inquiry, Health and Environmental Protection, Uranium Mining, Comissioner's Report, Volume 1. Oct 30, 1980. Berner, R.A. 1971. Principles of Chemical Sedimentology. New York: McGraw-Hill Book Company. 209 p. Borovec, Z., Kribek, B . and Tolar, V . Chemical Geology. 27: 39-46.  1979.  Sorption of Uranyl by Humic Acids.  Bouwer, Herman and Rice, R.C. 1976. Slug Test for Determining Hydraulic Conductivity of Unconfined Aquifers With Completely or Partially Penetrating Wells. Water Resources Research. 12: 423-428. Boyle, Dan R. 1982. The Formation of Basal-Type Uranium Deposits in South Central British Columbia. Economic Geology. 7:176-1209. Butler, James J. Jr., McElwee, Carl D., and Liu, Wenzhi. 1996. Improving the Quality of Parameter Estimates Obtained from Slug Tests. Groundwater. May/June: 480-490. Church, B.N., Jessop, A . M . , Bell, R., and Pettipas, G.S.B. 1991. Tertiary Outlier Studies: Recent Investigations in the Summerland Basin, South Okanagan Area, B . C . Geological Fieldwork 1990, Paper 1991-1. pp.163-170. Church, B.N. 1982. Geology of the Penticton Tertiary Outlier, B.C. Ministry of Energy, Mines and Petroleum Resources B.C. Revised Preliminary Map 35. Church, B . N . 1980. Anomalous Uranium in the Summerland Caldera. Fieldwork 1979, Paper 1980-1. pp. 11-15.  Geological  Church, B . N . 1973. Geology of the White Lake Basin. B.C. Department of Mines & Petroleum Resource Bulletin No. 61. 120 p.  110  Cooper, H . H . Jr., Bredehoeft, J.D., and Papadopulos, I.S. 1967. Response of FiniteDiameter Well to an Instantaneous Charge of Water. Water Resources Research. 3: 263269. Culbert, R.R. and Leighton, D.G. 313-330.  1988. Young Uranium. Ore Geology Reviews. 3:  Culbert, R.R., Boyle, D.R., and Levinson, A . A . 1984. Surficial Uranium Deposits in Canada. Surficial Uranium Deposits. Vienna: International Atomic Energy Agency. 252 p. Culbert, R.R. 1980. Report on Post-Glacial Uranium Concentrations in the Southern Okanagan Valley. For Ministry of Health B.C. June 4, 1980. Culbert, R.R. 1980b. Proceedings of Royal Commission of Inquiry, Health and Environmental Protection, Uranium Mining, Volume 69. p. 12662. Culbert, R.R. 1980c. Proceedings of Royal Commission of Inquiry, Health and Environmental Protection, Uranium Mining, Volume 69. p. 12626. Culbert, R.R. 1979. Post-Glacial Uranium Concentration in South Central British Columbia. For Royal Commission of Inquiry, Health and Environmental Protection, Uranium Mining. October 24, 1979. Culbert, R.R. and Leighton, D.G. 1978. Uranium in Alkaline Waters - Okanagan Area, British Columbia. CIMBulletin. May: 103-110. DeVeto, Richard H . 1978. Uranium Geology and Exploration - Lecture Notes and References. Golden, Colorado: Colorado School of Mines. Demir, Z. and Narasimhan, T.N. 1994. Improved Interpretation of Hvorslev Tests. Journal of Hydraulic Engineering. 120: 477-495. Drever, James I. 1988. The Geochemistry of Natural Waters. Englewood Cliffs, New Jersey: Prentice-Hall Inc. 437 p. Freeze, R.A. and Cherry, John A. 1979. Groundwater. Englewood Cliffs, New Jersey: Prentice-Hall Inc. 589 p. Froelich, P.N., Klinkhammer, G.P., Bender, M . L . , Luedtke, N . A . , Heath, G.R., Cullen, Doug, Dauphin, Paul, Hammond, Doug, Hartman, Blayne, and Maynard, Val. 1978. Early oxidation of organic matter in pelagic sediments of the eastern equatorial Atlantic: suboxic diagenesis. Geochimica et Cosmochimica Acta. 43: 1075-1090.  Ill  Giridhar, J. and Langmuir, D. 1991. Determination of E° for the U 0 / U Couple from Measurement of the Equilibrium: U 0 + C u + 4 H = U + C u + 2 H 0 at 25° C and some Geochemical Implications. Radiochimica Acta. 54:133-138. 2 +  4 +  2  2 +  2  +  (s)  4 +  2+  2  Golder Associates Ltd. 1994. Groundwater Monitoring Program - Summerland Landfill, Summerland, B C . For The Corporation of the District of Summerland. March 20, 1995. Golder Associates Ltd. 1980. Geotechnical Investigation, proposed injection well, proposed hotel development, College Road Summerland, B.C. For 5005 Holdings Ltd., c/o P. Avard, Naramata, B.C. April 7, 1980. Grant, M . B . and F.A. Michel. 1983. Study of the Hydrogeology of the White Lake Basin, British Columbia. For B . C . Department of Energy, Mines, and Resources. March, 1983. Contract Serial No. OSB83-00223. Hunkin Enginers, Inc. 1979. Preliminary Metallurgical Test Results and Conceptual Production Processes for Young Uranium Deposits with Construction and Operating Cost Estimates. For D.G. Leighton & Associates, Ltd., Vancouver, B.C. Hvorslev, M.J. 1951. Time lag and soil permeability in groundwater observations. 17.5. Army Corps Engrs. Waterways Exp. Sta. Bull. 36. Vicksburg, Miss. Idiz, Erdem F., Carlisle, Donald and Kaplan, I.R. 1986. Interaction between organic matter and trace metals in a uranium rich bog, Kern County, California, U.S.A. Applied Geochemistry. 1: 573-590. Jessop, A . M . and Church, B.N. 1991. Geothermal Drilling in the Summerland Basin, British Columbia, 1990. Geological Survey of Canada Open File 2348. 14p. Johnson, S.Y., Otton, J.K., and Macke, D.L. 1987. Geology of the Holocene surficial uranium deposit of the north fork of Flodelle Creek, northeastern Washington. Geological Society of America Bulletin. 98: 77-85. Kochenov, A . V . , Zinev'yev, and Lovaleva, S.A. 1965. Some Features of the Accumulation of Uranium in Peat Bogs. Geochemical International. 2: 65-70. Kribek, Bohdan and Podlaha, Jaroslav. 1980. The Stability Constant of the U 0 humic acid complex. Organic Geochemistry. 2: 93-97. 2  2 + 2  -  Kvill, Donald Raymond. 1976. Glacial History of the Trout Creek Basin, Summerland, British Columbia. Master of Science Thesis, Dept of Geography, University of Alberta. Langmuir, D. 1997. Aqueous Environmental Geochemistry. Upper Saddle River, New Jersey: Prentice Hall Inc. 590 p.  112  Langmuir, D. 1978. Uranium Solution-Mineral Equilibria at Low Temperatures with Applications to Sedimentary Ore Deposits. Geochimica et Cosmochimica Acta. 42: 547569. Lawson, D.W. 1968. Groundwater flow systems in the crystalline rocks of the Okanagan Highland, British Columbia. Canadian Journal of Earth Sciences. 5: 813-824. Lee, David R. and Cherry, John A. 1978. A Field Exercise on Groundwater Flow Using Seepage Meters and Mini-Piezometers. Journal of Geological Education. 27: 6-10. Leenheer, J.A., Malcolm, R.L., and White, W.R. 1974. Investigations of the reactivity and fate of certain organic components of an industrial waste after deep-well injection. Environmental Science and Technology. 10: 445-451 Levinson, A . A . , Bland, C J . and Dean, J.R. 1984. Uranium Series Disequilibrium in Young Surficial Uranium Deposits in Souther British Columbia. Canadian Journal of Earth Sciences. 21:559-566. Lewis, Trevor. 1984. Geothermal energy from Penticton Tertiary outlier, British Columbia: an initial assessment. Canadian Journal of Earth Sciences. 21: 181-188. Little, 1961. Geology of the Kettle River (West Half) area, British Columbia. Canadian Geological Survey Map 15-1961. Lopatkina, A.P. 1967. Conditions of accumulation of uranium in peat. Geochemistry International. 4: 577-588. Lovely, Derek R., Roden, Eric El, Phillips, E.J.P. and Woodward, J.C. 1993. Enzymatic iron and uranium reduction by sulfate-reducing bacteria. Marine Geology. 113: 41-53. Manskaya, S.M, Drozdova, T.V. and Emelianova, M.P. 1956. Association of uranium with humic acids and melanoidins. Geokhimiya. 4: 339-347. Meyboom, P. 1967. Groundwater studies in the Assiniboine River Drainage Basin: II. Hydrologic characteristics of phreatophytic vegetation in south-central Saskatchewan. Geol. Surv. Can. Bull 139. 64 p. Ministry of Environment, Lands and Parks, B.C. 1998. Groundwater Section website: Ministry of Environment, Lands and Parks, B . C . 1997. Memo to Mark Siemens, Summerland Trout Hatchery, Summerland, B.C. May 8, 1997. Ministry of Environment, Lands and Parks, B.C. 1986. Assessment of Water Resources at the Summerland Trout Hatchery, Report and Recommendations.  113  Ministry of Health, B.C. Uranium Deposits.  1981. Interim Report by the working group on Post-Glacial  Ministry of Health, B.C. and Ministry of Energy,.Mines and Petroleum Resources, B . C . 1981. Variations in Uranium and Radioactivity Levels in Surface and Groundwater at Selected Sites in British Columbia, April 1980-March 1981. 45 p. Ministry of Health, B.C. and Ministry of Energy, Mines, and Petroleum Resources, B.C. 1980. Variations in the Uranium and Radioactivity Levels of Potable Surface and Groundwater in the Okanagan and West Kootenay Regions of British Columbia. 63 p. Mohagheghi, A l i , Updegraff, David M . , and Goldhaber, Martin B . 1984. The Role of Sulfate-Reducing Bacteria in the Deposition of Sedimentary Uranium Ores. Geomicrobiology Journal. 4: 153-173. Nasmith, Hugh. 1962. Glacial History and Surficial Deposits of the Okanagan Valley, B C . B.C. Department of Mines and Petroleum Resources Bulletin No. 46. Owen, D.E. and Otton, J.K. 1995. Mountain wetlands: efficient uranium filters potential impacts. Ecological Engineering. 5:77-93. Otton, J.K., Zielinski, R . A . and Been, J.M. 1989. Uranium in Holocene Valley-fill Sediments and Uranium, Radon, and Helium in Waters, Lake Tahoe-Carson Range Area, Nevada and California, U.S.A. Environ. Geol. Water Sci. 13(1): 15-28. Pacific Hydrology Consultants Ltd. 1983. Groundwater Resources in the FaulderMeadow Valley Area. For Regional District of Okanagan-Similkameen. Nov 4, 1983. Piteau & Associates. 1984. Summerland Basin Hydrogeological Study. For Physics Branch, Energy, Mines and Resources Canada. January, 1984. Open File #85-29. Robbins, E.I., Zielinski, R.A., Otton, J.K., Owen, D.E., Schumann, R.R., and McKee, J.P. 1990. Microbially Mediated Fixation of Uranium, Sulfur, and Iron in a Peat-Forming Montane Wetland, Larimer County, Colorado. U.S.G.S. Research on Energy Resources 1990, Programs and Abstracts: pp. 70-71. Sandino, Amai and Bruno, Jordi. 1992. The solubility of (U0 ) (P0 ) -4H 0 and the formation of U(VI) phosphate complexes: Their influence in uranium speciation in natural waters. Geochimica et Cosmochimica Acta. 56: 4135-4145. 2  3  4  2  2  (s)  Schwertmann, U . , Kodama, H . , and Fischer, W.R. 1986. Mutual Interactions Between Organics and Iron Oxides. Interactions of Soil Minerals with Natural Organics and Microbes. Soil Science Society of America Spec. Pub. no 17. pp. 223-250.  114  Shanbhag, P . M . and Choppin, G. -19.81. Binding of uranyl by humic acid. J. Inorg. Nucl. Chem. 43: 3369-3372. : Shotyk, William. 1988. Review of the Inorganic Geochemistry of Peats and Peatland Waters. Earth Science Reviews. 25: 95-176. Steinmann, Philip, and Shotyk, William. 1997. Chemical composition, pH, and redox state of sulfur and iron in complete vertical porewater profiles from two Sphagnum peat bogs, Jura Mountains, Switzerland. Geochimica et Cosmochimica Acta. 61(6): 11431163. Stumrn, Werner and Morgan, James J. 1981. Aquatic Chemistry - An Introduction Emphasizing Chemical Equilibria in Natural Waters, 2 d Edition. New York: John Wiley & Sons Inc. n  Szalay, A . 1964. Cation exchange properties of humic acids and their importance in the geochemical enrichment of U 0 and other cations. Geochimica et Cosmochimica Acta. 28: 1605-1614. 2 +  2  Titayeva, N . A . 1967. Association of radium and uranium with peat. International. 4: 1168-1174.  Geochemistry  White, W.N. 1932. A method of estimating groundwater supplies based on discharge by plants and evaporation from soil. U.S. Geol. Surv. Water-Supply Paper 659-A. Zielinski, R . A . , Otton, J.K. Wanty, R.B. and Pierson, C T . 1988. The Aqueous Geochemistry of Uranium in a Drainage Containing Uraniferous Organic-rich Sediments, Lake Tahoe area, Nevada, U.S.A. Uranium. 4(4): 281-305. Zielinski, Robert A . , Otton, James K . , Wanty, Richard B . and Pierson, Charles. 1987. The geochemistry of Water Near a Surficial Organic-Rich Uranium Deposit Northeastern Washington State, U S A . Chemical Geology. 62: 263-289. Zielinski, Robert A . , Bush, C A . and Rosholt, J.N. 1986. Uranium series disequilibrium in a young surficial uranium deposit, northeastern Washington, U.S.A. Applied Geochemistry. 1: 503-511. Zielinski, Robert A . 1981. Experimental Leaching of Volcanic Glass: Implications for Evaluation of Glassy Volcanic Rocks as Sources of Uranium, pp. 1-10  115  A  HYDROGEOLOGICAL FIELD METHODS  This appendix describes the installation of the Prairie Flats piezometer network and the hydrogeological tests carried out in July, August, and September '97. It is divided into three main sections: 1) piezometer installation, 2) hydrogeological tests, and 3) stream measurements. PIEZOMETER INSTALLATION Soil Sampling During borehole excavation for piezometer installation, soil samples were taken to characterize the peat, clay and till units underlying the Prairie Flats. Grab samples of borehole cuttings were put in sealed plastic bags at piezometer locations 101/102, 302, 401/402, and 701/702. A list of the samples collected and their descriptions is found in Table A - l . Table A - l Grab sample descriptions Unit  Sample Location  Description  Till  101/102  Coarse sand and gravel fining downward with  401/402  rounded, pebble-sized clasts Grey Free of fines except at peat interface  Peat  101/102  Fine, powdery near surface, more decaying  701/702  plant material at depth Dark brown to black 0.5 to 2.5m thick  Clay  302  clay with some silt  402  blue-grey  116  Piezometer Construction Piezometers were made from sections of m-inch diameter, schedule 40 P V C pipe. The shallow and deep piezometers were roughly 1.5m and 3m in length, respectively. The deep piexometers were constructed from 1-metre flush-threaded sections, thereby allowing flexibility in length.  A n alternating pattern of one-eighth-inch holes was drilled into the pipe and covered with landscape filter fabric to make screens for the piezometers. The screened interval was 15cm long. The filter fabric was secured using plastic ties, and a plastic slip cap was fitted to the bottom. For piezometers 302 and 401, a fine nylon mesh was used in place of the filter fabric. A photo of these piezometers is shown in Plate A - l .  Plate A - l  A 1.5m-long piezometer with other lm-long piezometer sections in background  117  Borehole Digging Various digging tools were used to make boreholes for the piezometers. Choices depended on the desired depth of borehole and the type of material encountered. A description of each method is given below:  Rotating Dutch Auger A hand-held rotating dutch auger with a cutter head 4 inches in diameter and 25 cm long was used to drill shallow boreholes into the peat and clay for piezometers 102, 202, 302 and 602 (see Plate A-2). With extensions, it could attain a depth of 3 metres. Pebblesized clasts prevented the auger from penetrating the till unit beneath the peat, and lifting such water-logged, cohesionless material without losing it on its way up was impossible.  Plate A-2 Rotating dutch auger used to dig boreholes, 3m long with extensions  118  Post-hole Digger A two-man, gas-powered, post-hole digger with 4-inch solid stem was used to drill shallow boreholes into the peat and clay for piezometers 402, 702 and 502 (see Plate A 3). A 1-foot extention was fitted to a 4 foot flighted section to attain the desired depth of one metre.  Plate A-3 Gas-powered post-hole digger  119  Sand bailer The sand bailer (see Plate A-4) was the only effective means of removing the cohesionless till material beneath the peat unit, and was used at piezometers 101, 401, 601 and 701. A 3-inch diameter P V C casing (see Plate A-5) kept the formation open. The sand bailer consisted of a hollow aluminum cylinder 80 cm long and 8 cm in diameter with a one-way flap at its base. At its top an aluminum rod extended to a length of 4 metres. When pounded up and down inside the borehole, material became lodged inside the cylinder. Once full, it was brought up to surface, dumped out, and reinserted.  Plate A-4 Sand bailer with lm-long extensions  120  Plate A - 5 3-inch internal diameter P V C casing used during deep piezometer installation. The apparatus used to remove the casing lies on top o f it.  121  Piezometer Insertion Appendix B contains construction and stratigraphic details for each piezometer. Shallow piezometers were inserted deep enough so that their screened interval was below the water table, and deep piezometers were inserted as deep as possible.  Some piezometers were installed using a 3-inch diameter P V C casing and others not due to time and resource constraints. A 3-inch casing was used for piezometers 101, 102, 202, 302, 401, 601, 602 and 701 to prevent collapse of the peat around the borehole and to serve as a chute for sand, bentonite, and cuttings. First, the piezometer was slid into the open casing and allowed to fill with water. Next, clean, sorted silica filter sand was added to just above the screen. Bentonite chips were added to the top to create a 10cm seal, thereby isolating this depth interval for head measurements. This was especially important for the deep piezometers, but not so much for the shallow since these were completed just below the water table. The rest of the hole was filled with cuttings in the approximate order that they were removed, and more bentonite was added at the ground surface to create a final seal. With each filling step, the casing was pulled up to meet the new base of the hole and finally to remove it.  Piezometers installed without the use of casing include 402, 502, 702, 802 and 902. These were dropped into the open borehole and pushed into the soft saturated zone. No filter sand or bentonite seal was added around the screen. The borehole was filled with cuttings and sealed with bentonite at the ground surface.  122  Capping The final step in piezometer installation was capping. Excess standpipe was cut off with a hack saw, leaving about 10 cm above ground surface. A mark was made with a black marker on the top edge of the standpipe to indicate where water level measurements would be read and where the elevation of the piezometer would be surveyed in. A slip cap was put on and a 1/8-inch hole was drilled through the pipe and cap on each side. A piece of plastic-coated wire was threaded through these two holes and the ends were fastened with an electrical connector. This allowed air to move in and out of the standpipe with fluctuating water levels.  Surveying Two level surveys (the second as a check) were conducted using a Wild Heerbrugg No. 1 dumpy level to measure the elevations of the tops of the piezometers. A local benchmark was chosen and assigned an elevation of 10.000m, and all elevations were assigned relative to it. Angles and distances between piezometers were also measured to place them accurately on a map. Short distances were chained using a 30-metre nylon tape and chaining pins and horizontal angles were measured using a Pentax GT-48 transit. Range poles were positioned at each piezometer for easy sighting.  123  H Y D R O G E O L O G I C A L TESTS Head Measurements Several rounds of head measurements were made to draw piezometric contour maps of the site and to measure vertical hydraulic gradients. This data is presented in Table A-2. Depths to water were read relative to a black mark on the side of the piezometer standpipe using a water level tape, and relative head values were calculated using elevations from the level survey.  Slug Tests Slug tests were used to measure horizontal hydraulic conductivities of the peat and till units at most piezometer locations. For this method, a known volume of water is either withdrawn or added to the piezometer to cause an instantaneous change in head. Head recovery with time is recorded until it returns to its initial value. Head recovery data for the slug tests can be found in Appendix C.  Two rounds of slug tests were carried out over two visits to the site. On the first visit in July/August '97, injection tests were done at piezometers 202, 302, 402, and 902 using either tap water or a 1-inch diameter, 50cm long slug made of solid nylon. It is unlikely that this water affected the water chemistry results since ample time was allowed for reequilibration before sampling (two and seven months), and some of the highest uranium concentrations were found at 402 only two months later. More slug tests were performed in September '97 on piezometers 102, 602, 802, 902, 101, and 601, this time as withdrawal tests because of high water levels in the piezometers. For these a nylon bailer made from a 50cm length of 1-inch diameter nylon pipe was used.  124  CN  oo uo  CN ON NO  Tf  CO NO  ON  00  oo  od  ON  oo od od  uo SO  CN  SO  ON  NO CO CO  00 CO 00  uo  Tf  CN  CO 00  Tf  uo oo od od  00 NO  ,—I  00  00  od  ON  CN  VO  00 O  NO  Tf  ON  uo  ON  ON  Tf  NO Tf  r-H >—  1  ON  ON  ON  o  ^—1  oo  UO CO  00  Tf  UO  CN ON  uo oo uo  i—i  r-H *-H  ON  ON  ON  CO CO  uo  CN  CO  >—1  >—  Tf  ON  Tf  ruo  ON  ON  ON  ON  ON  ON  CN  Tf Tf  CN  oo  ON  NO NO CN  r—1  ON  ON  oo  00  ON  ON  r-  CN CO  Tf  © O  Tf  o  oo CN I  r-  PH  I CO ho CN  CO CN  i  CO 1  Tf  oo  q o  00  Tf  uo oo  ON  ON  ;  ON  ON  ON  00  o  ON  CN UO ON  ON  Tf  uo  ON UO  ON  ON  *—1  ON UO  CN O NO  ON 00  ON  ON  ON  ON  ON  00  O  CN NO  ON ON  r-  00  CN NO CN  Tf  NO ON  ON  ON  ON  ON  00  CN  so ,_;  UO  Tf  oo  00  >—  1  ON  ON  ON  ON  ON  ON  ON  00  NO  NO  UO  O  ON  Tf  Tf r-H  NO 00  uo  o  ON ON  ON  ON  00  00  ON  o\  Tf  Tf  NO  ON ON  oo od  od  ON  CO 00 NO  CN CO  ON  CO  .  uo uo  CN  CN  uo C N oo 0 0 od od  CO  Tf  Tf  ON  ON  t>  ON CN  UO  ON  uo CO  ,_; •ON NO  t>  UO  oo uo  ON  ON  ON ON CO  Tf  o  CO  ON  ON  ON  UO CN  o  Tf  CN  ON CO  ON  ON  ON  uo  CN  CO  r-  NO  q  ON  ON  UO CO  oo 00  CO CO NO  ON  od  00  o  Tf  NO  CN CO ON  uo  r~-  ON  ON  ON  od  od  CN  CN  ON  Tf  Tf  r-H  o  O  CN  o  ON  ON  os  oo  SO  Tf  UO  Tf  1  CN UO  U  r-oo od  o  t"-  CN  ON  ON  uo  UO  o  CN  ON  ON  00 NO  CN CO 00  oo od od  CN Tf  VO  uo O N r-; N O oo od  CO CN  NO  00  ON  uo  CN  l< 31  oo  ON  uo o  ON <—1 C O  ON  uo  Tf  so  SO SO  CN ON  ON  CN  ON  Tf Tf  CN 00  . r-H r-H  Tf  OH  OX) ,  00  uo uo od  CN  ON  Tf  o  00  ON  UO ON  ON  r-H  ICO  r-H  ON  T-H  oo od  ,_, SO  ON  uo  <D  NO  Tf (--  ON  NO UO  ON  UO  uo  ,_;  ON  ON NO ON  00  uo uo SO  ON  NO  ON CO  CN  i  CN  r-  NO  oo  Tf  ON  00  CN  o  r-H  00  00  uo  Tf  ON  (D ICO  Tf  NO NO ON  OS od  r-H  o  OH  ll  00  1—I  CN  oo 00  Tf Tf  NO ON  60  o  ON  CO  O  UO NO ON  ON  ON  o  00 =3  CN  oo  r-H  r-H  t-Tf  oo 0 0 od od  CO O;  CN ON NO  CN  ON NO CN  uo  o  00  00  ON  ON  od od  uo  NO  SO oo t> od od  o  ON 00  r—1  SO 00  uo uo  uo  o  ON 00  ON  o  Tf  '  UO ON  SO 00  q  r—1 r-H  NO CN  r-H  r-  Tf  o o  UO  Tf  uo  r-H  UO  Tf  CO NO  uo  ©  q  ON  ON  ON  ON  oo od  ON  ON  UO  CN UO NO  r-  CN  ON  NO  ON  ON  od  od  UO CO  Tf  ON  ^H  ,_;  uo uo o o  r—1  ON  ON  CO  r—1  NO  r—1  r-H  NO  . od  uo  s  CO  uo ©  r—1 r-H  s  UO O co r ON  -fa  U  S3  00  od od  CO C3  <U  a  o  s r-H  o  r-H  IS  CN  CN  CN  r-H  CN  CO  O  N N PH PH  O  N PH  r-H  CN  O  o  N PH  N N PH PH  Tf  O  Tf  CN  O uo  N PH  r-H  o  CN  r-H  CN  O  O O  N N PH PH  N N PH PH  NO  NO  cd  CN  O oo  N PH  CN ON  ml  N PH  *  Recent work by Butler and others (1996) recommends that three or more slug tests be carried out at a well during a particular test period, and using at least two different initial head displacements. This is to identify any evolving "skin effects" or dependence on initial head displacement. Therefore, three withdrawal tests with two different initial displacements were done at piezometer 602.  Borehole Infiltration tests in peat The vertical hydraulic conductivity of the peat was measured by means of an infiltration test. A 30 cm deep hole was bored into the peat using the dutch auger, into which a section of 3-inch diameter P V C casing was inserted to the bottom of the hole. The open cross-section at the bottom of the casing permits vertical flow only. A n initial head displacement of 28cm was made by pouring water into the casing, and falling head measurements were made over a two-day period.  126  STREAM MEASUREMENTS Vertical hydraulic gradient The vertical component of the hydraulic gradient between Prairie Creek and the underlying groundwater at points B, C, and D was measured using a miniature piezometer inserted in the streambed (see Figure A - l ) . The piezometer was made of a 0.3 cm diameter polyethylene tube with a perforated end (screen) wrapped with fine nylon mesh to prevent sediment intrusion. It was 2.20 m long to the midpoint of the screen. To install it, a 2 cm diameter steel pipe was driven into the streambed with a sledge hammer, the base of which was plugged with a steel bolt. The piezometer was inserted into the steel pipe and then the pipe removed, leaving the piezometer in place. For further details on the use of mini-piezometers see Lee and Cherry (1978).  J  dh  I  Figure A - l Diagram of mini-piezometer installation in a streambed  127  The distance over which the vertical component of the hydraulic gradient was measured (dl) was calculated by subtracting the exposed length of piezometer above the streambed (/) from the total length of the piezometer (L). The head difference (dh) was the difference in elevation between the stream water level and that inside the piezometer. Hydraulic gradient (i) was calculated as dh/dl, as shown in Table A-3.  Table A-3 Measurement of the vertical component of the hydraulic gradient in Prairie Creek Location Date  Time  L(m)  /(m)  dl(m)  dh(m)  i  B  July 19/97  17:00  2.20  0.88  1.33  0.12  0.13  C  July 22/97  12:50  2.20  0.60  1.60  -0.04  -0.06  D  July 20/97  17:00  2.20  0.44  1.76  -0.14  -0.32  Seepage Rate A seepage meter measures the seepage rate between groundwater and a stream bed (see Figure A-2). The base of the seepage meter consists of a 57 cm diameter metal drum with a 3 cm diameter hole drilled on the top and side. The base is pushed into the sediment about 10 cm and allowed to fill with water. Next, a one-hole stopper fitted with a 40 cm-long polyethylene tube is inserted into the top hole to serve as a vent. Another stopper is inserted into the side hole to which a plastic bag is attached with electric tape, taking care that it is well-sealed and completely submerged (thereby maintaining the same head as the stream). Where hydraulic gradients indicated groundwater discharge into the stream, the plastic bag was inserted empty and allowed to fill over a set time period. Where the opposite flow direction was measured, the bag was filled with a known quantity of water and allowed to drain. This technique is described in Lee and Cherry (1978).  Seepage rate was calculated as dV/dtA, where dV is the volume of water gained or lost, dT is the total time elapsed and A is the cross-sectional area of the seepage meter, as is tabulated below: 128  V  i  Figure A-2 Diagram of seepage meter installation in a streambed  Table A-4 Seepage meter measurements in Prairie Creek Location Date Time of insertion Time of removal Elapsed time (hrs) (dT) Change in volume (mL) (dV) Flow rate (mL/hr) area of seepage meter (cm ) (A) Seepage rate (mL/hr/cm ) 2  3  B July 19/97 16:40 17:50 1.17 870 745.71 2570 0.29  C July 21/97 14:15 14:55 0.67 -190.00 -285.00 2570 -0.11  B PIEZOMETER LOGS  P I E Z O M E T E R 101  Depth(m)  completed in sandy till Elevation of standpipe (w/o cap) (m)  9.250  Total depth ot piez. (m) Depth to eff midscreen (m) Avg depth to water (m)  2.730 2.585 0.300  Inside diameter (cm) Cross-sectional area (cm2) Purge volume (mL)  3.4 9.1 7435.9  Screen length (cm) Effective screen length (cm) Length of piezometer (m)  15.0 15.0 2.950  Z = elevation head (m)  6.665  Completion Date Drilling method  Jul 22 97 auger & bailer  0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3  Lithology  peat  peat and clay  pebbly sand fining up  130  P I E Z O M E T E R 102 m)  completed in peat Elevation of standpipe (w/o cap) (m)  9.250  Total depth of piez. (m) Depth to eff. midscreen (m) Avg depth to water (m)  1.120 0.975 0.300  Inside diameter (cm) Cross-sectional area (cm2) Purge volume (mL)  3.4 9.1 3050.6  Screen length (cm) Effective screen length (cm) Total length of piezometer (m)  15.0 15.0 1.545  Z = elevation head (m)  8.275  Completion Date Drilling Method  Lithology  see 101 mostly peat some clay  screen  Jul 22 97 Hand Auger  131  PIEZOMETER 202 completed in sandy till w/ sand pack Elevation ot standppe (w/o cap) (m)  Depth(m) 9.745  total depth ot piez. (m)  1.022  Depth to etr. midscreen (m) Avg depth to water (m)  0.877  Inside diameter (cm) Cross-sectional area (cm2) Purge volume (mL)  3.4 9.1 2783.7  1.540  Z = elevation head (m)  8.868  Drilling method  peat  0.700  Screen length (cm) Effective screen length (cm) total length ot piezometer (m)  Completion Date  Lithology  pebbly sand  15.0 ?  Jul-23-97 auger,  sand pack clay and peat cuttings bentonite seal  bailer, posthole digger  132  PIEZOMETER 302 completed in clay w/ sand pack  Depth(m)  Elevation of standpipe (w/o cap) (m)  10.450  Total depth of piez. (m) Depth to efi. midscreen(m) Avg depth to water (m)  1.845 1.700 1.100  Inside diameter (cm) Cross-sectional area (cm2) Purge volume (mL)  3.4 9.1 5025.3  Screen length (cm) Effective screen length (cm) Total length of piezometer (m)  15.0 40.0 2.050  Z = elevation head (m)  8.750  Completion Date  31-M-97  Drilling method  auger only  0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2  Lithology  I peat  Iclay  sand pack  clay  clay and peat cuttings  sand  gradin  bentonite seal  133  P I E Z O M E T E R 401 completed in pebbly sand w/ sand pack  Depth(m)  Lithology  Elevation of standpipe (w/o cap) (m) Total depth of piez. (m) Depth to eft! midscreen (m) Avg depth to water (m) Inside diameter (cm) Cross-sectional area (cm2) Purge volume (mL) Screen length (cm) Effective screen length (cm) Length of piezometer (m) Z = elevation head (m)  Completion Date Drilling method  134  P I E Z O M E T E R 402 completed in peat  Depth(m)  Elevation o f standpipe (w/o cap) (m)  9.630  Total depth o f piez. (m)  1.262  Depth to erf midscreen (m)  1.117  Avg depth to water (m)  0.700  Lithology  see 401  Inside diameter (cm)  3.4  Cross-sectional area (cm2)  9.1  Purge volume (mL)  3437.4  Screen length (cm)  15.0  Effective screen length (cm)  ?  Total length o f piezometer (m)  1.540  Z = elevation head (m)  8.513  Completion Date Drilling method  Jul-23-97 posthole  sand pack clay and peat cuttings bentonite seal  digger  135  P I E Z O M E T E R 502 Elevation o f standpipe (w/o cap) (m)  Lithology  Depth(m)  completed in peat 11.900  0.0 0.1  Total depth o f piez. (m)  1.305  0.2  Depth to effi nidscreen (m)  1.160  0.3  77.000  0.4  Avg depth to water (m)  peat  0.5 Inside diameter (cm)  3.4  0.6  Cross-sectional area (cm2)  9.1  0.7  3554.5  0.8  Purge volume (mL)  0.9 Screen length (cm) Effective screen length (cm)  15.0  1.0  ?  Total length o f piezometer (m)  1.540  Length of cutoff  0.235 sand pack  Z = elevation head (m)  10.740 clay and peat cuttings  Completion Date  Jul-23-97 bentonite seal  Drilling method  post-hole digger  136  P I E Z O M E T E R 601 Depth(m)  completed in sand w/ sand pack Elevation of standpipe (w/o cap) (m)  9.720  Total depth of piez. (m) Depth to efT. midscreen (m) Avg depth to water (m)  1.385 1.240 0.690  Inside diameter (cm) Cross-sectional area (cm2) Purge volume (mL)  3.4 9.1 3772.4  Screen length (cm) Effective screen length (cm) Length of piezometer (m)  15.0 26.5 1.540  Z = elevation head (m)  8.480  Completion Date Drilling method  Lithology  July 30'97 auger, bailer  137  P I E Z O M E T E R 602 completed in sandy  till/clay w/ sand pack  Depth(m)  Elevation o f standpipe (w/o cap) (m)  9.725  Total depth o f piez. (m)  1.115  Depth to efE nidscreen (m)  0.970  Avg depth to water (m)  0.690  Inside diameter (cm) Cross-sectional area (cm2)  Lithology  see 601  3.4 9.1  Purge volume (mL)  3037.0  Screen length (cm)  15.0  Effective screen length (cm) Total length o f piezometer (m)  1.965  Length o f cutoff (m)  0.850  Z = elevation head (m)  8.755  Completion Date  M-97  sand pack clay and peat cuttings bentonite seal Drilling method  auger  138  P I E Z O M E T E R 701 completed in sandy till w/ sand pack  Depth(m)  Lithology  [Elevation of standpipe (w/o cap) (m) Total depth of piez. (m) Depth to effi irridscreen (m) Avg depth to water (m) Inside diameter (cm) Cross-sectional area (cm2) Purge volume (mL) Screen length (cm) Effective screen length (cm) Length of piezometer (m) Z = elevation head (ni)  Completion Date Drilling method  2.8  2.9 3  139  P I E Z O M E T E R 702 Elevation of standpipe (w/o cap) (m)  Lithology  Depth(m)  completed in peat 10.465  Total depth of piez. (m)  1.287  Depth to efl midscreen (m)  1.142  Avg depth to water (m)  1.000 see 701  Inside diameter (cm) Cross-sectional area (cm2)  3.4 9.1  Purge volume (mL)  3505.5  Screen length (cm)  15.0  Effective screen length (cm)  ?  Total length of piezometer (m)  1.540  Z = elevation head (m)  9.323  sand pack  Completion Date  Jul 23 97  clay and peat cuttings  Drilling method  post-hole  bentonite seal  digger  140  P I E Z O M E T E R 802 Completed in peat Elevation of standpipe (w/o cap) (m)  Lithology  Depth(m)  9.580  0.0 0.1  Total depth of piez. (m)  1.285  0.2  Depth to eff. midscreen (m)  1.140  0.3  Avg depth to water (m)  0.700  0.4 0.5  Inside diameter (cm) Cross-sectional area (cm2)  3.4  0.6  9.1  0.7  Purge volume (mL)  3500.0  0.8  Screen length (cm)  15.0  1.0  Effective screen length (cm)  15.0  1.1  Total length of piezometer (m)  1.540  1.2  0.9  1.3 Z = elevation head (m)  8.440  1.4  o o  0  0 o O  sand pack Completion Date  31-M-97 clay and peat cuttings  Drilling method  auger bentonite seal  141  P I E Z O M E T E R 902 epth(m)  completed in peat Elevation of standpipe (w/o cap) (m)  9.122  Total depth of piez. (m)  1.201  Depth to eft rnidscreen (m)  1.056  Avg depth to water (m)  0.650  Inside diameter (cm) Cross-sectional area (cm2)  3.4 9.1  Purge volume (mL)  3271.2  Screen length (cm)  15.0  Effective screen length (cm)  15.0  Total length of piezometer (m)  1.540  Z = elevation head (m)  8.066  Completion Date Drilling method  Lithology  31-M-97 auger  142  C  H Y D R A U L I C C O N D U C T I V I T Y TEST D A T A  BAIL TEST Date  Aug2/97  P I E Z O M E T E R 101  Actual Time 9:50:00 10:00:00 10:02:00 10:04:00 10:06:00 10:15:00 10:21:00 10:26:00 11:06:00 11:10:00 11:16:00 11:24:00 11:45:00 13:13:00 14:03:00 16:30:00  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m) 0:00 0 0:02:00 0:04:00 0:06:00 0:15:00 0:21:00 0:26:00 1:06:00 1:10:00 1:16:00 1:24:00 1:45:00 3:13:00 4:03:00 6:30:00  120 240 360 900 1260 1560 3960 4200 4560 5040 6300 11580 14580 23400  0.496 2.540 2.352 2.230 2.115 1.680 1.406 1.290 0.642 0.601 0.542 0.513 0.487 0.494 0.498 0.509  8.754 6.710 6.898 7.020 7.135 7.570 7.844 7.960 8.608 8.649 8.708 8.737 8.763 8.756 8.752 8.741  Cooper-Bredehoeft  Hvorslev Plot  Plot  Time(s)  • •  r  (H-h)) (H-Ho)  •  •  •  •  • •  •  ±  •  •  •  •  Time (s)  143  P I E Z O M E T E R 101  BAIL TEST Date  Sep27/97  Actual Time 1 50:00 1 51:30 1 52:00 1 53:00 1 54:00 1 54:30 1 55:00 1 55:30 1 56:00 1 57:00 1 58:00 2 00:00 2 10:00 2 20:00 3 07:00  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m) 0:00 0:00 0:30 0:01:30 0:02:30 0:03:00 0:03:30 0:04:00 0:04:30 0:05:30 0:06:30 0:08:30 0:18:30 0:28:30 1:15:30  Hvorslev  0:00 30 90 150 180 210 240 270 330 390 510 1110 1710 4530  0.258 0.606 0.415 0.333 0.282 0.271 0.266 0.264 0.263 0.263 0.262 0.260 0.259 0.259 0.266  Plot  8.992 8.644 8.835 8.917 8.968 8.979 8.984 8.986 8.987 8.988 8.990 8.991 8.991 8.984 9.250  Cooper-Bredehoeft Plot  •  o  X  o T  mlm  •  I c  Time (s)  Time (s)  144  P I E Z O M E T E R 102  BAIL TEST Date  Sep27/97  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m)  Actual Time 1 30:00 1 30:30 1 31:00 1 32:00 1 33:00 1 34:00 1 35:00 1 40:00 1 45:00 1 50:00 2 00:00 2 10:00 2 20:00 3 09:00 4 10:00  0:00 0:00:00 0:00:30 0:01:30 0:02:30 0:03:30 0:04:30 0:09:30 0:14:30 0:19:30 0:29:30 0:39:30 0:49:30 1:38:30 2:39:30  Hvorslev  0:00 30 90 150 210 270 570 870 1170 1770 2370 2970 5910 9570  0.274 0.622 0.528 0.526 0.522 0.520 0.517 0.506 0.498 0.490 0.474 0.459 0.445 0.399 0.363  Plot  8.976 8.628 8.722 8.724 8.728 8.730 8.733 . 8.744 8.752 8.776 8.791 8.805 8.851 8.887 9.250  Cooper-Bredehoeft Plot  ln(H-h)jf(H-Ho)  • •  •  • • • •  •  2  •  • •  •  •  >  Time (s)  Time (s)  145  PIEZOMETER 202  SLUG TEST Date  Aug3/97  Actual Time 19:38:00 19:38:00  20:03:00 20:13:00 20:15:00 20:37:00 21:01:00  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m) 0 0 15 25 40 50 60 70 80 100 125 160 180 195 225 255 300 360 480 600 1500 2100 2220 3540 4980  0:00 0:00:15 0:00:25 0:00:40 0:00:50 0:01:00 0:01:10 0:01:20 0:01:40 0:02:05 0:02:40 0:03:00 0:03:15 0:03:45 0:04:15 0:05:00 0:06:00 0:08:00 0:10:00 0:25:00 0:35:00 0:37:00 0:59:00 1:23:00  i i  •  •  •  •  Time (s)  8.990 9.460 9.425 9.347 9.247 9.275 9.259 9.251 9.245 9.242 9.235 9.230 9.226 9.225 9.221 9.217 9.212 9.207 9.195 9.186 9.129 9.102 9.096 9.063 9.041  Cooper-Bredehoeft Plot  Hvorslev Plot  2"  0.755 0.285 0.320 0.398 0.498 0.470 0.486 0.494 0.500 0.503 0.510 0.515 0.519 0.520 0.524 0.528 0.533 0.538 0.550 0.559 0.616 0.643 0.649 0.682 0.704  P I E Z O M E T E R 302  SLUG TEST Date  Aug3/97  Actual Time  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m)  0:00 18:50:00 18:54:00 0:00:00 0:00:30 0:00:40 0:00:50 0:01:05 0:01:20 0:01:45 0:02:00 0:02:30 0:03:00 0:03:15 0:03:30 0:04:00 0:04:30 0:05:00 0:05:30 0:06:00 0:06:45 0:08:00 0:09:15 0:11:00 0:13:00 0:14:30 0:17:30 19:14:00 0:20:00 0:24:15 19:28:00 0:34:00 19:50:00 0:56:00 19:56:00 1:02:00 19:58:00 1:04:00 20:05:00 1:11:00 20:23:00 1:29:00  0 0 30 40 50 65 80 105 120 150 180 195 210 240 270 300 330 360 405 480 555 660 780 870 1050 1200 1455 2040 3360 3720 3840 4260 5340  1.194 0.694 0.700 0.704 0.710 0.715 0.720 0.726 0.732 0.738 0.746 0.749 0.752 0.758 0.764 0.771 0.777 0.784 0.793 0.807 0.821 0.84 0.86 0.875 0.905 0.922 0.948 1.009 1.09 1.105 1.11 1.12 1.145  9.256 9.756 9.750 9.746 9.740 9.735 9.730 9.724 9.718 9.712 9.704 9.701 9.698 9.692 9.686 9.679 9.673 9.666 9.657 9.643 9.629 9.610 9.590 9.575 9.545 9.528 9.502 9.441 9.360 9.345 9.340 9.330 9.305  Hvorslev Plot  •  Time (s)  Cooper-Bredehoeft Plot  •  •  Time (s)  P I E Z O M E T E R 402  S L U G TEST Date  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m)  Actual Time 18:29:00 18:29:00  Aug3/97  0 0 30 40 45 75 105 135 165 195 240 270 300 330 360 390 420 450 495  0:00 0:00 0:00:30 0:00:40 0:00:45 0:01:15 0:01:45 0:02:15 0:02:45 0:03:15 0:04:00 0:04:30 0:05:00 0:05:30 0:06:00 0:06:30 0:07:00 0:07:30 0:08:15  0.813 0.524 0.525 0.531 0.540 0.547 0.557 0.564 0.572 0.578 0.585 0.594 0.599 0.602 0.608 0.613 0.618 0.622 0.627  8.817 9.106 9.105 9.099 9.090 9.083 9.073 9.066 9.058 9.052 9.045 9.036 9.031 9.028 9.022 9.017 9.012 9.008 9.003  Cooper-Bredehoeft Plot  Hvorslev Plot • • •  o X  1  X  1c"  • • • •  o X  • •  X c  •  Time (s)  •  •  •  T i m e (s)  149  P I E Z O M E T E R 601  B A I L TESTS  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m)  Actual Time  Date  Sep27/97  10:43:00 10:43:00 10:44:00 10:44:30 10:45:00 10:45:30 10:46:00 10:46:30 10:47:00 10:49:00 10:50:00 11:00:00 11:05:00 11:10:00 11:30:00  0:00 0:00:00 0:01:00 0:01:30 0:02:00 0:02:30 0:03:00 0:03:30 0:04:00 0:06:00 0:07:00 0:17:00 0:22:00 0:27:00 0:47:00  11 33:00 11 33:00 11 33:10 0:00:10 11 33:30 0:00:30 11 34:00 0:01:00 11 34:30 0:01:30 11 35:00 0:02:00 11 37:00 0:04:00 11 40:00 0:07:00 15 27:00 3:54:00  60 90 120 150 180 210 240 360 420 1020 1320 1620 2820  0.250 0.598 0.420 0.340 0.269 0.260 0.257 0.255 0.255 0.255 0.254 0.258 0.257 0.257 0.257  9.470 9.122 9.300 9.380 9.451 9.460 9.463 9.465 9.465 9.465 9.466 9.462 9.463 9.463 9.463  10 30 60 90 120 240 420 3240  0.257 0.493 0.354 0.305 0.273 0.266 0.259 0.258 0.258 0.285  9.463 9.227 9.366 9.415 9.447 9.454 9.461 9.462 9.462 9.435  Hvorslev Plots  ln(H-h)ir(H-Ho)  A  Cooper-Bredehoeft Plots  • • •  o  A  • •  A  T  •  • •  A A  Time (s)  •  s  • A  • A  Time (s)  150  P I E Z O M E T E R 602  B A I L TESTS Date  Sep26/97  Sep27/97  Actual Time 3:52:00 3:52:00 3:53:00 3:54:00 3:55:00 3:56:00 3:57:00 3:58:00 3:59:00 4:00:00 4:01:00 4:03:00 4:05:00 4:10:00 4:15:00 4:20:00 4:25:00 4:30:00 4:40:00 4:50:00 4:55:00 6:53:00 10 54:00 10 55:00 10 55:30 10 56:00 10 57:00 10 58:00 10 59:00 11 00:00 11 02:00 11 04:00 11 06:00 11 08:00 11 10:00 11 15:00 11 20:00 11 25:00  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m)  0:00:00 0:01:00 0:02:00 0:03:00 0:04:00 0:05:00 0:06:00 0:07:00 0:08:00 0:09:00 0:11:00 0:13:00 0:18:00 0:23:00 0:28:00 0:33:00 0:38:00 0:48:00 0:58:00 1:03:00 3:01:00  0:00:30 0:01:00 0:02:00 0:03:00 0:04:00 0:05:00 0:07:00 0:09:00 0:11:00 0:13:00 0:15:00 0:20:00 0:25:00 0:30:00  60 120 180 240 300 360 420 480 540 660 780 1080 1380 1680 1980 2280 2880 3480 3780 10860  0.198 0.839 0.805 0.778 0.751 0.728 0.704 0.683 0.664 0.641 0.625 0.590 0.560 0.493 0.441 0.400 0.365 0.337 0.291 0.267 0.258 0.225  9.527 8.886 8.920 8.947 8.974 8.997 9.021 9.042 9.061 9.084 9.100 9.135 9.165 9.232 9.284 9.325 9.360 9.388 9.434 9.458 9.467 9.500  30 60 120 180 240 300 420 540 660 780 900 1200 1500 1800  0.251 0.524 0.512 0.503 0.484 0.469 0.456 0.441 0.417 0.395 0.379 0.361 0.347 0.32 0.299 0.286  9.474 9.201 9.213 9.222 9.241 9.256 9.269 9.284 9.308 9.330 9.346 9.364 9.378 9.405 9.426 9.439  11:30:00 11:40:00 11:50:00 12:00:00 12:21:00 15:27:00  Sept27/97  0:35:00 0:45:00 0:55:00 1:05:00 1:26:00 4:32:00  15:27:00 15:32:00 15:32:30 15:33:00 15:34:00 15:35:30 15:36:00 15:37:30 15:38:00 15:39:00 15:40:00 15:45:00 15:50:00 16:02:00  0:00:30 0:01:00 0:02:00 0:03:30 0:04:00 0:05:30 0:06:00 0:07:00 0:08:00 0:13:00 0:18:00 0:30:00  2100 2700 3300 3900 5160 16320  0.276 0.262 0.259 0.257 0.257 0.286  9.449 9.463 9.466 9.468 9.468 9.439  30 60 120 210 240 330 360 420 480 780 1080 1800  0.286 0.522 0.514 0.505 0.489 0.47 0.465 0.449 0.444 0.436 0.427 0.392 0.366 0.328  9.439 9.203 9.211 9.220 9.236 9.255 9.260 9.276 9.281 9.289 9.298 9.333 9.359 9.397  Cooper-Bredehoeft Plots  Hvorslev Plots  V • •A  • • •  1c* X  •  • • •  •  Time (s)  H/Ho  X X  • -  %  *x .x  Time (s)  152  P I E Z O M E T E R 802  BAIL TEST Date  Sep27/97  Actual Time 14:27:00 14:27:00 14:27:30 14:28:00 14:28:30 14:29:00 14:30:00 14:31:00 14:32:00 14:34:00 14:36:00 14:38:00 14:40:00 14:45:00 14:50:00 15:00:00 15:15:00  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m)  0:00:30 0:01:00 0:01:30 0:02:00 0:03:00 0:04:00 0:05:00 0:07:00 0:09:00 0:11:00 0:13:00 0:18:00 0:23:00 0:33:00 0:48:00  30 60 90 120 180 240 300 420 540 660 780 1080 1380 1980 2880  0.520 0.756 0.753 0.745 0.740 0.735 0.726 0.721 0.713 0.696 0.683 0.673 0.662 0.640 0.622 0.595 0.565  9.060 8.824 8.827 8.835 8.840 8.845 8.854 8.859 8.867 8.884 8.897 8.907 8.918 8.940 8.958 8.985 9.015  Cooper-Bredehoeft Plot  •  • • •  •  • •  •  Time (s)  153  P I E Z O M E T E R 902  SLUG TEST Date  Aug2/97  Actual Time 0.4416667 0.4420139 0.4423611 0.4430556 0.44375 0.4451389 0.4465278 0.4479167 0.4493056 0.45625 0.4604167 0.4673611 0.4770833 0.4784722 0.4861111 0.5013889 0.5541667 0.7111111  Elapsed Time 0 0 0.0003472 0.0010417 0.0017361 0.003125 0.0045139 0.0059028 0.0072917 0.0142361 0.0184028 0.0253472 0.0350694 0.0364583 0.0440972 0.059375 0.1121528 0.2690972  Elapsed Depth to Total Head Time (s) Water (m) Water (m) 8.645 8.931 8.922 8.91 8.901 8.888 8.873 8.863 8.855 8.818 8.805 8.782 8.759 8.752 8.746 8.724 8.684 8.653  0.477 0.191 0.2 0.212 0.221 0.234 0.249 0.259 0.267 0.304 0.317 0.34 0.363 0.37 0.376 0.398 0.438 0.469  0 30 90 150 270 390 510 630 1230 1590 2190 3030 3150 3810 5130 9690 23250  Cooper-Bredehoeft Plot  Hvorslev Plot  A *  </>  A A  2  o -TLP.  =£  I  \ A  A A  A A  T i m e (s)  A  T i m e (s)  154  PIEZOMETER 902  BAIL TEST Date  Sep27/97  Actual Time 9:17:00 9:17:00 9:17:30 9:18:00 9:19:00 9:20:00 9:21:00 9:23:00 9:25:00 9:27:00 9:35:00 9:40:00 9:45:00 9:55:00 10:05:00 10:15:00 10:30:00 15:04:00  Elapsed Elapsed Depth to Total Head Time Time (s) Water (m) Water (m)  0:00:30 0:01:00 0:02:00 0:03:00 0:04:00 0:06:00 0:08:00 0:10:00 0:18:00 0:23:00 0:28:00 0:38:00 0:48:00 0:58:00 1:13:00 5:47:00  30 60 120 180 240 360 480 600 1080 1380 1680 2280 2880 3480 4380 20820  0.291 0.531 0.530 0.529 0.526 0.524 0.522 0.518 0.511 0.506 0.490 0.481 0.473 0.455 0.439 0.425 0.409 0.321  8.831 8.591 8.592 8.593 8.596 8.598 8.600 8.604 8.611 •8.616 8.632 8.641 8.649 8.667 8.683 8.697 8.713 8.801  Cooper-Bredehoeft Plot  Hvorslev Plot  •  o  o  X X  -- T L -.  • •  x c  •  T i m e (s)  T i m e (s)  155  D  WATER SAMPLING FIELD METHODS  This appendix describes how water samples were collected from the piezometers in September '97 and March '98, and which chemical parameters were measured in the field. It is divided into three sections: 1) piezometer preparation 2) sample collection and 3) field analyses.  PIEZOMETER PREPARATION Immediately after installation, a piezometer is developed to clear away sediment from the base of the standpipe and around the screen. Just before taking a groundwater sample, purging is done to remove stagnant or non-representative water in the piezometer casing, filter pack, and formation. This section explains how development, purging, and sample retrieval were carried out on site. Piezometer Development A l l piezometers were developed using a Waterra™ inertial pump, applying high flow rates to free as much material as possible. The shallow piezometers were pumped to dryness three times. Due to the low hydraulic conductivity of the peat, this process required up to two days. Water clarity improved significantly (from grey to colourless) at all piezometers except 202. It is possible that the filter cloth covering the screen was pushed up or torn when the piezometer was inserted, allowing silt and clay to enter the piezometer.  At least three well volumes (one well volume is the volume of water in the standpipe) were pumped from each of the deep piezometers, stopping when good water clarity was achieved.  To avoid cross-contamination by the Waterra™ pump, the foot valve and  attached hose were rinsed with a 15% nitric acid solution, then twice with tap water between sampling locations.  156  Piezometer Purging Before taking a water sample from the deep piezometers, at least five well volumes were purged using a peristaltic pump at flow rates of around 750 mL/min. The pump was powered by a car battery and connected to a flow-through cell (see Plate D - l ) . The cell consisted of a plastic erlenmyer vacuum flask capped with a rubber stopper having three holes:  one for incoming water, another for the pH or Eh probe, and another for the  temperature probe. Water left the cell via the vacuum outlet. Separate lengths of tubing were inserted into each piezometer to prevent cross-contamination.  Water clarity was  good at all three locations and conductivity, pH, and temperature stabilized within minutes. Due to the slow recovery rates of the shallow piezometers, purging just prior to sample collection was not feasible.  Plate D-l  Set-up o f peristaltic pump (geopump) and flow-through cell (centre) w i t h p H meter  attached  157  SAMPLE COLLECTION  Samples were collected from the deep piezometers using the peristaltic pump once pH, conductivity, and temperature readings in the flow-through cell had stabilized (see Plate D-l). At the shallow piezometers, samples were collected using either a Waterra™ pump (September '97) or nylon bailer (March '98) within 24 hours of development, or as soon as the water levels had returned to equilibrium.  Oxygenation of the water in the  standpipe was not a concern since relatively high levels of ferrous iron (over 3 mg/L) were detected. To avoid cross-contamination, both the Waterra™ pump and bailer were rinsed twice with tap water and once with distilled water between sampling locations.  FIELD ANALYSES  The geochemical parameters measured in the field included Eh, ferrous iron, pH, conductivity, temperature, and alkalinity.  At all deep piezometers, Eh, pH, and  temperature were measured in the flow-through cell, and conductivity and alkalinity were measured on bailed samples. At the shallow piezometers, all parameters were measured on bailed samples. Further description of the methods and equipment used is given below. Eh  Redox potential was measured using a Cole-Palmer ORP (Redox) Combination Electrode probe with built in reference electrode and platinum band. This was hooked up to the p H meter. Before departure to the field, the Eh probe was tested in the lab with Quinhydrone solutions made up in p H 7 and pH 4 buffers. In the field, the platinum band was cleaned regularly to prevent fouling by organics and other agents in the peat waters, and care was taken to minimize sample exposure to air.  158  Ferrous Iron Ferrous iron was measured using a H A C H Ferrous Iron test kit. This method uses a 1,10 Phenanthroline indicator which turns the sample orange upon reaction with ferrous iron. Concentrations of up to 10 mg/L, in increments of 0.2 mg/L were read using a colour wheel. A l l containers were rinsed twice with tap water and once with distilled between stations. Care was also taken to minimize exposure to air.  pH and Temperature A n Orion liquid-filled pH probe connected to a Hanna HI 9025C p H meter was used to measure pH. Temperature was measured with a temperature probe connected to the pH meter. The p H meter was calibrated before each use using buffers 4.01 and 7.01, and rinsed with distilled water between readings.  Care was taken to minimize sample  exposure to air, since oxidation, degassing, and temperature increase affect pH.  Conductivity A hand-held conductivity meter (Hanna Dist WP 3) was inserted into the sample to measure conductivity. It had been pre-calibrated in the lab using standardized KC1 solutions.  Bicarbonate As an approximate measure of bicarbonate (HC0 ) concentrations in these waters, the 3  A S T M standard test method for alkalinity of water (D 1067-92) by electrometric titration was followed.  This assumes that most of the acid-neutralizing capacity of the water  comes from bicarbonate, which is reasonable given the neutral p H of these groundwaters and surface waters. Fifty mililitres of sample were titrated with small additions of 0.02N hydrochloric acid to a pH of 4.0.  Next, p H was plotted against the volume of acid added, and the bicarbonate inflection point (around pH 4.3) was interpolated. A l l titration data is plotted in Figure D-l. Assuming that all of the acid-neutralizing capacity comes from dissolved carbonate, and  159  knowing that bicarbonate is the dominant species at these neutral pH's, the bicarbonate concentration was calculated using:  mg/L HCO3  = (mL of titrant)(0.02 meq/L)(61 mg/meq)(1000 mL/L) mL sample = (mL of titrant) (24.4)  Replicate alkalinity titrations were done at every fifth sampling location to test the reproducibility of the bicarbonate values.  160  a> E t— o. W >.  CO CD  c  a.  Q CN O  r O  CM  o  CM  CM  CO  CO  o  T  o  CO c CD  CL  CO  CM O  CN O  CN  cn  CN  o  o CN  o  LU  CO CO CD  CN O CO  CO  3  I  <  cu CD  c £  ca  <  i  CN O CD  t £ t S ^ co  p | O V "|UJ  ^  E LAB ANALYSES OF CHEMICAL PARAMETERS This appendix summarizes the protocols followed in the lab and in the field for water samples requiring lab analysis. It is divided into two sections: 1) sample preservation and storage and 2) lab instrumentation and analysis.  SAMPLE PRESERVATION AND STORAGE  Bottled water samples were collected for laboratory analysis of dissolved uranium, dissolved cations (Ca, Mg, Na, K , F e  X0T  ), dissolved anions ( N 0 , P 0 , Cl, and S0 ), and 3  4  4  dissolved carbon (total, organic, and inorganic). Two sets of samples were taken, one for uranium and cations, and the other for cations and dissolved carbon, since each group shares the same preservation and storage requirements. Uranium and Cations (U, Ca, M g , Na, K , F e  TOT  )  Before departure to the field, all sampling and storage material was acid-washed in a 20% nitric acid solution, rinsed with distilled water, air-dried in a fume hood, and sealed in zip-lock plastic freezer bags.  In the field, water samples were syringe-filtered using  sterilized 0.45 pjn Millex-Ha Millipore filters, and stored in 30mL or 60mL nalgene bottles. Before capping, they were acidified to a p H < 2 using 50% v/v nitric acid. Syringes that had to be used more than once were soaked for 3 hours in a 20% nitric acid solution, then in distilled water, then sealed in a clean zip-lock bag. Back at the lab, these samples were placed in the refrigerator at 4°C until analysis.  Anions N O „ SQ , C l , P Q and Dissolved Carbon 4  4  In September '97, nitrate concentrations were measured using a Hach test kit. Samples were filtered in the field using a hand pump and vacuum flask fitted with Whatman 1 Qualitative filter paper. They were stored in lOOmL nalgene bottles in a cooler at 4°C,  162  and analyzed that same day.  A l l sampling and storage equipment was rinsed with  distilled water before and after use, with care taken to avoid contact with nitric acid.  In March '98, nitrate concentrations were measured using an ion chromatograph, along with other anions. Before departure to the field, new 60 mL nalgene bottles were rinsed twice with deionized water, air-dried in a fume-hood, and sealed in clean zip-lock bags. At the time of collection, samples were syringe-filtered through sterilized 0.45 um Millex-Ha Millipore filters and then bottled.  They were frozen within 6 hours of  collection, and thawed only in time for analysis.  L A B INSTRUMENTATION AND ANALYSIS Analyses for uranium, major cations (Ca, Mg, Na, K ) , anions ( N 0 , S 0 , P 0 , Cl) and 3  4  4  dissolved carbon (TC, IC, OC) were done in the lab within six weeks of sample collection. Descriptions of the laboratory procedures and instrumentation are given in this section.  Uranium Aqueous uranium ( U) concentrations were measured by inductively-coupled mass 238  spectrometry (ICP-MS). ICP-MS is ideally suited to analysis of trace heavy metals since it is rapid and sensitive to concentrations in the parts per billion range.  Also,  interferences by element isotopes of similar mass are not a problem with uranium.  First, external and internal standard solutions, and blank solutions were prepared. External standards are known concentrations of uranium whose measured response signal is plotted to create a calibration curve. Solutions of 1 u.g/L, 10 |J.g/L, and 25 (J.g/L uranium in a 1% nitric acid matrix were prepared from a 1000 mg/L standard reference solution. One percent is roughly the same concentration of nitric acid in the field samples after acidification to pH < 2. A l l dilutions and additions were carried out using precalibrated Eppendorf micropipetters. 163  Internal standards are known amounts of some non-analyte that are added to all samples in order to detect and correct for fluctuations in matrix composition, instrumental drift or blockages, and other processes that may cause loss of sample during analysis.  Indium  was chosen as the internal standard because its natural occurrence in the environment is negligible, it lies in the same mass range as uranium, it is easily ionized in the ICP, and it does not exhibit any interfering chemical effects.  A 1 mg/L Indium solution was  prepared from a 1000 mg/L Indium reference standard, then added to external standards, samples, and lab blanks to a concentration of 20 ug/L. The laboratory blank consisted of the 1% nitric acid matrix plus internal standard.  During analysis, the ICP-MS was set in pulse-counting mode, the number of counts being proportional to the concentration of analyte in the sample.  First the external standards  were measured, followed by the samples in order of expected uranium concentration. Pulse counts were later converted to uranium concentration values using the calibration curve.  To ensure quality results, a number of measures were taken:  • sample replicates were done every 5 samples, and the lab blank and 10 ug/L standard were run every 10 samples. Repeating the 10 ug/L standard allowed any shift in the calibration curve through the course of analysis to be detected.  • a 30 second instrument rinse with 1% nitric acid was done between samples, but i f remnant U in the sample loop was suspect, the lab blank was run again as a rinse.  • uranium values were corrected for instrumental drift, based on fluctuations in internal standard read-outs. The correction formula used is:  Corrected U = Uncorrected U * Average In counts / In counts 164  Anions ( N O , , SO/, C l , P Q / ) Nitrate samples collected in September '97 were analyzed by Hach's Medium-Range Nitrate (0-4.5 mg/L N 0 - N ) and Low-Range Nitrate (0-0.15 mg/L N 0 - N ) test methods. 3  3  These used powder pillows and a D R 2000 Spectrophotometer, and are based on the Cadmium Reduction Method outlined in Standard Methods for Wastewater Analysis. A N a N 0 solution of known concentration was used as a control standard. Water samples 3  from 102 and 402 underwent some orange tinting due to iron-oxide precipitation, which may have biased the results high.  For the March '98 water samples, concentrations of nitrate, sulfate, chloride, and phosphate were obtained by ion chromatography within two weeks of sample collection. This work was carried out by a lab technician in the Oceanography department at U B C , using a Dionex DX-100 Ion Chromatograph. As an accuracy check, two reference solutions of known concentration were analyzed every 10 samples.  Cations (Ca , M g , N a , K , F e ^ ) 2+  2 +  2+  +  Samples were analyzed by Flame Atomic Absorption (AA) and Atomic Emission (AE) using a ThermoJarell Ash Video  22 A A / A E  spectrophotometer  in the  Civil  Environmental Lab at U B C . Flame A A was used for iron, magnesium, and calcium; and A E was used for sodium and potassium.  External standards were prepared for Ca, Mg, Na, K and Fe by making dilutions of 1000 mg/L standard reference solutions. Their concentrations were chosen to straddle the expected range of cation concentrations in the samples while remaining within the instrument's optimal measurement range.  Dilutions were made using glass pipettes and  volumetric flasks, and a 1% nitric acid matrix was maintained for all standards and diluted samples.  165  During the analysis, samples were ordered from least to most concentrated. The machine was re-zeroed and re-calibrated every 10 samples, and lab replicates were also done every 5 samples. Calcium, iron, and magnesium concentrations were calculated automatically by the instrument, whereas non-linearity in the calibration curves for sodium and potassium required that each be interpolated manually. This was not a concern since nonlinear calibration curves are often encountered when operating in A E mode.  Dissolved Carbon A Shimadzu TOC-500 Total Carbon / Inorganic Carbon Analyzer was used to find concentrations of total dissolved carbon (TC), dissolved inorganic carbon (IC) and dissolved organic carbon (OC) in the samples. It uses the Combustion Infrared Method, described in Standard Methods for Wastewater Analysis. individually, and OC was calculated by difference.  TC and IC were measured  This analysis was done by a  technician in the Civil Environmental Lab at U B C .  Four standards were used to calibrate the instrument: 10 and 100 mg/L TC, and 10 and 100 mg/L IC.  As a precision check, two samples were analyzed twice, and as an  accuracy check, two other samples were spiked with known amounts of analyte.  When the carbon analyses were carried out, the samples had previously been thawed and refrozen. On the second thawing, a reddish-brown precipitate was observed in many of the samples, guessed to be a mixed precipitate of iron oxide, organics and silica. If they do contain organics, the dissolved carbon concentrations should not be interpreted quantitatively.  166  QUALITY CONTROL RESULTS  Various measures were taken to ensure quality results at all stages of sampling and analysis. This included the use of blanks, field duplicates, lab replicates, and control standards.  The results of this work are summarized in Table E - l . The left side gives the percent relative error values calculated for the control standards, lab replicates, and field duplicates. Absolute concentrations of analytes measured in the blanks is given on the right. Areas of concern are highlighted and discussed in the upcoming sections. Control Standards Where possible, solutions with known amounts of analyte were used to check the accuracy of results, particularly that of the standards. A l l relative errors were less than 8%, which is acceptable. Lab Replicates Lab replicates are the same samples analyzed twice in the lab. Agreement between these was consistently above 95%, which is excellent and therefore eliminates concern about instrument precision. Replicates were not carried out during anion analysis because these were done by another technician. Field Duplicates Field duplicates are separate samples taken at the same time from the same location. Duplicates were taken at piezometers 101, 401, 502, and at the exit point of Prairie Creek from the flats. The agreement between duplicates at 401 and 502 is excellent, but at 101 calculated errors are a consistent 25%, and at the stream, 100%. The variation in the Prairie Creek samples is not surprising since they weren't taken at exactly the same time from exactly the same spot along the stream.  167  The seriousness of the 25% variability seen at 101 is lightened somewhat by the fact that concentrations of most analytes at 101 are among the lowest measured in these groundwaters, some even near detection limits.  Thus, the relative error between  duplicates at other locations is likely to be lower. Field and Laboratory Blanks Five kinds of blanks were prepared for both the acidified and non-acidified samples: the Development blank, Bailer & Bottle blank, Syringe blank, Trip blank, and Cleaning blank. A brief description of each is given below:  • The Development blank consisted of tap water rinsed through the Waterra™ tubing right after the development of piezometer 402.  Its purpose was to detect any cross-  contamination by the Waterra™ during development and sampling.  • The Bailer and Bottle blank consisted of distilled water washed through the nylon bailer and sampling bottle after prior rinses with tap water-. It would detect any contamination of the sample by these items.  • The Cleaning blank was made up of distilled water that had been used in the final soak of bottle cleaning. This made sure that negligible amounts of analyte remained after cleaning.  • The Syringe blank was made up of distilled water passed through the syringe and filter after an acid cleaning in the field.  It detects any carry-over of analyte after cleaning  (especially N 0 ) . 3  • The Trip blank is distilled water brought out to the field, opened, acidified if required, and stored with the rest of the samples. This measures contamination during the entire process of transport, sampling, storage, and analysis.  168  Of all of these blanks only the Development blank had significant amounts of analyte. Its concentrations represent less than 15% of the measured concentrations in the groundwater at 402, except for phosphate (which is already near the detection limit) and chloride. Chloride has a higher concentration in the blank than in the groundwater, due to either analytical error or to chloride already present in the tapwater used to rinse the Waterra™ tubing. Even i f cross-contamination were occurring, pumping piezometers dry during development leaves virtually no water to contaminate fresh inflowing water. In general, the blanks results show that contamination of the water samples during all phases of field work and analysis is negligible. Error Bars  The largest contribution to error in the water chemistry data is sampling irreproducibility, since lab replicate errors are negligible by comparison. Therefore, error bars are based on the field duplicates, except for temperature, pH, Eh, and Fe , which are based on the 2+  measuring instrument or method. Ranges of error are summarized in Table E-2.  169  o  =8  H3  O  O  O  O  O  O  CN  —  o  o  o  o  o  o  o  —  o  ~  r o -d- —  —  —  oo  o  o  o  o  o  o  o  o  o  o  o  o  ^cd  ,cd  O  O  o  o  o  o  t/3 o  ^ 2  —1 PQ  o  o  o  —  6J0  cd  CN  o  O  _cd  c s oo o  —c  o  o  o  o  cn  CN  —  O  O  ro O  !_!  23 £  iO  cn  »i  »i  (N  CN  CN  CN  CN  O  O  O  n  oi  O  1  'd-  => <=  _cd  ,cd  ^cd  o  -g  2  2  vi f~  £ R  VO  2  S i l l  - °° - a - a « -a  o  o  o 8 cd H  St****  ,cd  ^cd  a a a  V~i  t—  CN  2  Q io  ,cd  cn  1-1  z; o  -£-§-£  8 *> ">  od  -—  -  a a a a  ro  cn  N"  -I  ,cd  CN  CN  3  ,cd  a e a a a e  P  u  .cd  —•  S55o  o o  Table  E-2 Error bars assigned to each chemical parameter  Parameter  Error  Temperature Conductivity pH Eh (measured)  ± 0.5 °C ± 10 uS/cm  Fe x  ± 5%  N0 P0  ±10% ±10% ± 10% ± 10%  TO  so  3  4  4  Cl Ca Mg ~Na K U  ± 0 . 1 pH unit ±mV ± 0.2 mg/L  ±5% ±5% ±5% ±5%  HCO3  ± 10% ±5%  Total Carbon (TC) Inorganic Carbon (IC) Organic Carbon (OC)  ± 10% ± 10% ± 20%  (20% at PZ 101) (20% at PZ 101)  (20% (20% (20% (20%  at at at at  PZ PZ PZ PZ  101) 101) 101) 101)  F  GEOCHEMICAL MODELING USING PHREEQC AND  WATEQ4F  DATABASE  This appendix explains the revisions made to the WATEQ4F database used by the PHREEQC program to do the geochemical modeling for this thesis. Aqueous species concentrations and mineral saturation indices calculated for the groundwater samples collected at the Prairie flats are also presented.  Revisions made to W A T E Q 4 F database  Since the original compilation of the WATEQ4F database in the 1980's, several advances have been made in uranium chemistry. The U0 (HP0 )2 ~ species has been discarded 2  2  4  (Sandino and Bruno, 1992), and the E° value for the U 0  2 + 2  /U  4 +  couple has been re-  measured (Giridhar and Langmuir, 1991). This has caused significant changes in the Gibbs free energy values for U , U O H , U 0 , and U(OH) ° and some U 4 +  3 +  +  2  4  4 +  minerals,  which have been accounted for in the thermodynamic tables found in Langmuir (1997, p. 551-2). Since the logK and AH values for other U  4 +  and U  6 +  species in the WATEQ4F  database matched those given by Langmuir (1997), the standard state conditions were also assumed to be consistent.  Therefore, the WATEQ4F database was updated  according to Langmuir (1997) through the following revisions.  U0  2  U0  2  U  + 2P0 " + 2 H = U 0 ( H P 0 ) 3  +  4  2 +  2  + 4 H + 2e" = U +  + H 0=UOH  4 +  2 + 2  + e" = U 0  4 +  +H  3 +  2  U0 U  2 +  +  log K = -0.65 log K = 1.49  +  2  4  +  2  am  U0 c) 2(  + 4H = U +  4 +  eliminated log K = 8.89  2  2  U0 ( ) + 4H = U  2 2  + 2H 0  + 4 H 0 = U(OH) ° + 4 H  4 +  4  +  log K =-12.0  2  log K = 4.16  + 2H 0  logK = -4.99  4 +  + 2H 0 2  172  Prairie Flats Geochemical Modeling  Mineral saturation indices and concentrations of aqueous species at equilibrium were calculated for each groundwater sample in the March '98 data set using P H R E E Q C and the revised WATEQ4F database. Fe  2+  Input parameters included pH, calculated Eh (from  concentrations), temperature, major anion concentrations (N03, C l , PO4, and SO4),  major cation concentratons (Ca, M g , Na, K , Fe oT), and uranium concentrations (see T  Table 4-5). The output is given in Table F-l:  Table F-l PHREEQC output for Prairie Flats groundwaters Location  Eh (V)  pE  SI U0  SI 2  U02(am)  SI U 0 4  SI 9  uo 3  SI Calcite  U-CO3  U-PO4  M  M  8  PZ 101  0.236  4.0  -8.0  -16.4  -26.2  -14.8  0.2  3.9E-08  1.3E-11  PZ 102  0.041  0.7  1.5  -6.7  3.1  -4.0  -0.4  1.5E-07  3.3E-10  PZ 202  -0.024  -0.4  2.6  -5.6  6.0  -3.7  0.6  1.6E-05  4.2E-11  PZ 302  0.018  0.3  1.2  -7.1  1.9  -4.9  0.1  3.5E-07  1.8E-10  PZ 401  0.236  4.0  -7.2  -15.5  -23.3  -13.0  0.7  2.8E-06  3.5E-11  PZ 402  0.035  0.6  2.2  -6.1  5.6  -2.3  -0.2  1.5E-06  7.5E-10  PZ502  -0.041  -0.7  3.1  -5.2  7.9  -2.5  0.1  4.6E-07  1.2E-12  PZ 601  0.077  1.3  0.3  -8.0  -0.8  -5.6  -0.5  9.4E-08  3.7E-10  PZ 602  0.047  0.8  1.6  -6.6  3.7  -3.2  -0.3  4.3E-07  7.5E-10  PZ 701  0.236  4.0  -7.6  -15.8  -24.8  -14.4  0.2  2.9E-07  1.5E-11  PZ 702  -0.035  -0.6  1.8  -6.5  3.2  -5.3  0.2  2.7E-07  1.2E-11  -0.9  4.2E-08  4.7E-10  -3.1  2.4E-08  3.5E-09  PZ 802  0.018  0.3  2.3  -5.8  5.8  -2.7  PZ 902  0.177  3.0  -0.6  -8.8  -2.2  -4.1  * U - C 0 and U-PO4 refer to the total concentrations of all uranyl-carbonate and uranyl-phosphate species, respectively 3  ** SI is the Saturation Index with respect to the given mineral phase  Results of a sensitivity analysis testing the effects of pH, HCO3 concentration, Eh, and U concentration on uranium mineral saturation indices are tabulated below:  173  Table F-2  Output of PHREEQC sensitivity analysis on U mineral saturation  pH  Eh (V)  HC0 mg/L  U ug/L  SI U 0  7.53  -0.1 0 0.1 0.2  555  75.4  3.99 0.61 -2.77 -6.17  9.76 -0.38 -10.52 -20.72  6 6.5 7 7.53 8  0  555  75.4  4.81 3.79 2.25 0.61 -0.82  13.39 10.31 5.13 -0.38 -5.13  7.53  0  100 300 500 700  75.4  2.64 1.39 0.74 0.3  7.76 2.74 0.17 -1.61  7.53  0  555  10 50 100 500 1000  -0.27 0.43 0.73 1.43 1.73  -3.89 -1.09 0.12 2.91 4.12  3  2  SI U 4 O 9  


Citation Scheme:


Usage Statistics

Country Views Downloads
Canada 15 6
United States 7 0
Japan 4 0
Saudi Arabia 3 0
Swaziland 1 0
Italy 1 0
City Views Downloads
Unknown 6 1
Kelowna 5 3
Tokyo 4 0
Calgary 2 0
Burnaby 2 2
Wilmington 2 0
Summerland 2 0
Ashburn 2 0
Lecce 1 0
Mbabane 1 0
Chicago 1 0
Kenner 1 0
Seattle 1 0

{[{ mDataHeader[type] }]} {[{ month[type] }]} {[{ tData[type] }]}
Download Stats



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items