Open Collections

UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

Mixing processes from CTD profiles using a lake-specific equation of state : Quesnel Lake James, Christina 2004

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


831-ubc_2005-0058.pdf [ 4.74MB ]
JSON: 831-1.0063438.json
JSON-LD: 831-1.0063438-ld.json
RDF/XML (Pretty): 831-1.0063438-rdf.xml
RDF/JSON: 831-1.0063438-rdf.json
Turtle: 831-1.0063438-turtle.txt
N-Triples: 831-1.0063438-rdf-ntriples.txt
Original Record: 831-1.0063438-source.json
Full Text

Full Text

MIXING PROCESSES FROM CTD PROFILES USING A LAKE-SPECIFIC EQUATION OF STATE: QUESNEL L A K E by  CHRISTINA JAMES B . S c , University of British Columbia, 2001  A THESIS' S U B M I T T E D IN P A R T I A L F U L F I L M E N T O F T H E REQUIREMENTS FOR T H E D E G R E E O F M A S T E R O F APPLIED SCIENCE  in T H E F A C U L T Y O F G R A D U A T E STUDIES  (Department of Civil- Engineering; Environmental Fluid Mechanics)  U N I V E R S I T Y O F BRITISH C O L U M B I A  November 2004  © C h r i s t i n a James, 2004  Abstract Quesnel Lake, is a deep (511m maximum depth) fjord-type lake i n northeast British Columbia, Canada. M i x i n g processes in the lake exchange deep-water with surface water and contribute to the renewal of surface-water nutrients and oxygenated deep-water. These processes are of great consequence to the lake's trophic dynamics and understanding them will enable better management of the large salmon resources in Quesnel Lake. To better understand large-scale convective processes, a lake-specific equation of state was developed. Water samples were collected at locations around Quesnel Lake and analysed for ionic and non-ionic composition as well as other quantities that are integral to determining the lake's equation of state including p H , alkalinity and specific conductance.  A relationship was developed to  find lake water salinity from C T D data. Salinity was i n t u r n related to density using a modified form of a general limnological equation of state. T h e equation of state developed for Quesnel Lake gives densities accurate to ± 0.0018 k g / m  3  whereas the general equation of state (based on seawater composition) is only accurate to ± 0.0158kg/m for Quesnel Lake water samples. 3  The lake-specific equation of state was used to identify gravitational instability i n density profiles estimated from C T D data. In order to compare water parcel density within a profile, the hydrostatic pressure effect must be removed. The three quantities that are used for this purpose, potential density, quasidensity and standard density, were compared. Quasi-density was found to be most appropriate for Quesnel Lake's deep water which is near the temperature of maximum density. Quesnel lake water column stability was quantified using the Brunt-Vaisala frequency calculated using quasi-density.  ii  Contents Abstract  ii  Table of Contents  .'til,  List of Figures  vii  List of Tables  viii  List of Notations  ix  Acknowledgements  xi  1  Introduction  1  2  Literature Review  4  2.1 2.2  2.3  2.4  3  Other Works related to Equations of State State Variables 2.2.1 Temperature 2.2.2 Salinity 2.2.3 Pressure Relationships 2.3.1 Relating Conductivity to Specific Conductance 2.3.2 Relating Specific Conductance to Salinity 2.3.3 Relating Salinity to Density Profile Stability 2.4.1 Potential Density 2.4.2 Quasi-Density 2.4.3 Standard Density 2.4.4 Brunt-Vaisala Frequency  Methods 3.1 Water Samples 3.2 Calculating Conductivity 3.3 Finding Salinity 3.4 Determining the Coefficient of Haline Contraction 3.5 Measuring Density 3.6 C T D Profiles 3.7 Applying Quesnel Lake's Equation of State 3.7.1 Quasi-Density  iii  4 5 5 9 12 13 13 15 16 17 18 18 20 21 23 23 26 30 31 33 34 34 34  3.7.2 3.7.3 4  Standard Density Brunt Vaisala Frequency  36 36  Quesnel Lake's Equation of State  38  4.1 Equation for f • 4.2 Equation for Salinity 4.2.1 Finding the Constant of Proportionality, A 4.2.2 Salinity Equation 4.3 Equation of State for Quesnel Lake 4.3.1 Coefficient of Haline Contraction 4.3.2 Quesnel Lake's Equation of State 4.3.3 Error T  5  ;  . . . . .  38 41 41 42 43 43 45 45  Local Stability  50  5.1 5.2 5.3 5.4  50 55 57 66  Quasi-Density Standard Density Instability Identification Density Currents  ,  6  Conclusion  70  7  Future Work  73  7.1 Thorpe Scale 7.2 M9 Thermistor Chain 7.3 Further Investigation of Density Currents  73 74 75  References  76  A Polynomial Coefficients  B  81  A',1 Limiting Equivalent Conductivity A.2 Reduction Factor  81 82  The Carbonate System  83  C Matlab Equation of State Toolbox  C.l C.2 C.3 C.4 C.5 C.6  densitycalc.m relationships.m C25.m carbonates.m equivalents.m LEC.m  86  '.  iv  88 89 91 93 94 95  C.7 reductioncoeff.m C.8 keff.m . . C.9 salinity.m  ••  v  96 98 99  List of Figures 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35  Map of Quesnel Lake Density Dependence on Temperature Illustration of a Stable Profile that crosses TMD TMD d Tf i for seawater Density dependence on Salinity TMD dependence on Pressure Quasi-Density Illustration Schematic of Brunt Vaisala Frequency Map of water sample stations, July 2003 Composition of Quesnel Lake water samples Map of water sample stations, June 2004. Sample fits for Aj and f Concentrations of Si vs. C a Temperature factor fx function for Quesnel Lake Measured and Calculated specific conductance for Quesnel Lake . . . Measured and Calculated Specific Conductance for several B C lakes. . C25 vs Salinity for determining A Comparing density with and without Quesnel Lake salinity contribution Measured and calculated densities for Quesnel Lake Measured and calculated densities for Quesnel Lake Compared with  2 6 8 9 11 13 19 22 23 26 27 29 38 39 40 41 42 44 48  EOSCM  49  a n  reez  ng  • • •  Comparison of Potential and Quasi-Density Profiles Quasi-Density found using E O S C M and E O S Q L Quasi-Density Profiles for 2001, 2002, and 2003 Quasi-Density Profiles for for 2003 throughout Quesnel Lake Standard density calculated varying T : Density as a function of Temperature and Pressure (2D) B i n Average Quasi-Density Power spectrum for 2002 Quasi-density profile Brunt Vaisala Frequency found using Various B i n Sizes Brunt Vaisala Frequency found using E O S C M and EOSQL Brunt Vaisala Frequency found using Potential Density and QuasiDensity Profiles of the Brunt-Vaisala Frequency for 2001, 2002, and 2003 . . . Profiles of the Brunt-Vaisala Frequency for 2003 throughout Quesnel Lake Illustration of density current Map showing location of Mooring 9 •. Q  vi  50 52 53 54 55 56 58 59 60 61 62 64 65 69 74  36  The Carbonate System  83  37  Matlab Toolbox Flowchart  87  vii  List of Tables 1 2  Depth of Water Samples, July 2003 Summary of Tests on Quesnel Lake water samples  24 25  3 4  Depths of Water Samples, June, 2004 Stoichiometric combination of ions.  27 32  5 6  Anton Paar Densitometer calibration values Depths at which thermisors are attached to Mooring 9 in Quesnel Lake's East A r m  46  viii  75  List of Notations List of Symbols A  Ci C C .... D fi fp.... fx . . . . m .... N .... p r R .... S S .... T T .... TMD • • x z z' z  the coefficient of proportionality between specific conductance and salinity percentage offset of Quesnel Lake coefficient of haline contraction to coefficient of haline contraction in equation of state developed by Chen and Millero (1986) measured concentration of ion i conductivity conductivity measured at temperature T and pressure P standard density reduction coefficient of ion i pressure factor temperature factor concentration of a dissolved solid in Quesnel Lake Brunt Vaisala frequency or stability frequency pressure conductivity ratio coefficient of determination salinity reference salinity used in standard density temperature reference temperature used in standard density temperature of maximum density square root of concentration of all dissolved solids in Quesnel Lake depth integration variable for depth used in quasi^density reference depth used in quasi-density and potential density  a 3 7 Aj . . . . 'A°°.... A ^ p.....  coefficient of thermal expansion coefficient of haline contraction temperature-dependent coefficient for compressibility equivalent conductivity of ion i limiting equivalent conductivity of ion i equivalent conductivity of a salt adiabatic lapse rate of density density  B  T  2  0  0  r  ix  p .... PCM • •  mean density density calculated with the equation of state developed by Chen and Millero (Chen and Millero, 1986) Pp ... potential density Pquasi • • quasi-density PQuesnei density calculated with the equation of state developed specifically for Quesnel Lake cr p-1000, [kg/m ] 6 potential temperature 0  0t  3  List of Abbreviations EOSCM  C T D .. DDW . PESC . PVT .. EOSQX  Equation Of State developed by Chen and Millero (Chen and Millero, 1986) ' Conductivity - Temperature - Depth probe De-ionised Distilled Water Pacific Environmental Sciences Centre Pressure Volume Temperature properties Equation Of State developed specifically for Quesnel Lake  x  Acknowledgements I would like to express my gratitude to the following people for their assistance in the completion of this thesis. I would like to thank Dr. Bernard Laval, my supervisor, who has been a great source of ideas and guidance through this process. I am very grateful to the researchers from the Institute of Ocean-Sciences, Fisheries and Oceans Canada including Dr: Eddy Carmack, John Morrison and Dr. Svein Vagle who have helped me carry out my field work and develop some of the ideas presented here. Fellow students including Daniel Potts and Ryan North also helped with field work. Roger Pieters from the Department of Earth and Ocean Sciences, University of British Columbia spent much time discussing the finer points of this project with me. Thank you also to Dr. Svein Vagle of the Institute of Ocean Sciences, the Provincial Ministry of Water Land and A i r Protection and Fisheries and Oceans Canada's Cultus Lake Laboratory for providing ship time. I was supported by Fisheries and Oceans Canada under the Science Subvention Research grant. Dr. Greg Lawrence from the Department of C i v i l Engineering, University of British Columbia also provided generous assistance in this project. I am also very thankful to Brock Wilson and Sabina and Douglas James, for their love through this process.  xi  1  Introduction  Quesnel Lake is a part of the Fraser River system and is an important habitat for salmon. Historically, Quesnel Lake received up to one third of the entire Fraser River run of salmon as they spawn in the lake or pass through it on their way to their natal streams. This makes Quesnel Lake a gateway to the catchment which contributes the largest salmon population to the Fraser River watershed (Royal, 1966). After the rock slide at Hell's Gate in 1913, the salmon population in Quesnel Lake was reduced to the point where there were too few salmon to make observations (Thompson, 1945). In recent years, salmon escapement has returned to near record numbers increasing to 3.5 million in 2001 (Fisheries and Oceans Canada, 2001). Thus, a critical management issue concerns the lake's primary productivity, and the role lake mixing and convection plays in returning nutrients to the surface (euphotic) layer of the lake.  The deep  water of Quesnel Lake has oxygen levels near or at saturation suggesting that deep water is being ventilated by oxygenated surface water. Also, nutrients from bottom waters are, in some way, being recycled to the surface where these nutrients then support primary production. Therefore, a motivation for studying Quesnel Lake is to determine what mechanisms contribute to vertical mixing and deep water renewal and how these mechanisms affect the Quesnel Lake salmon resource. Quesnel Lake, located on the western edge of the Cariboo mountains in British Columbia (BC), Canada, is a steep-sided fjord-type lake consisting of three long, narrow arms (Fig 1). With a maximum depth of 511 m, it is the deepest lake in B C and one of the deepest lakes in the world. It has a volume of approximately 41.5 k m  3  and a surface area of 266.3 k m . Quesnel Lake has three major in-flowing rivers, the 2  Horsefly River, Mitchell River and Niagara Creek and one major outflow, the Quesnel River. Conductivity - temperature - depth (CTD) profilers are routinely used for oceanographic and limnological research. These instruments measure the conductivity, temperature and pressure as they descend through the water column. In order to study  1  Figure 1: Map of Quesnel Lake. Data from Campbell (2001). Place names are used in text.  vertical mixing in Quesnel Lake, density profiles, calculated from in-situ C T D measurements, were used. To find density profiles from C T D data, a lake specific equation of state must be employed. The equation of state relates a fluid's density, p, to its temperature, T , salinity, S, and pressure, p:  p = p(T,S,p)  (1)  While T and p are both directly measured by a C T D , lake-specific relations are required to estimate salinity from in-situ measurements.  Once a salinity relation  has been established, a modified general equation of state (e.g. Chen and Millero, 1986) can be used to determine water column density. Regions of mixing can often be identified from density profiles by isolating sections where density inversions, and therefore instabilities, occur (Hohmann et al., 1997). Identifying where Quesnel Lake's  2  water columns become unstable is an integral part of understanding how this deep (> 250m) lake mixes vertically. This paper will describe the processes by which water samples were taken from Quesnel Lake and analysed in order to determine the lake-specific equation of state for Quesnel Lake, hereafter referred to as EOSQZ,. Chapter 2 will outline work done by previous researchers on this topic.  Chapter 3 will give detail about the field  work that was carried out at Quesnel Lake along with the method by which the data was analysed.  Chapter 4 will show the results of this analysis, namely the  developed E O S Q L . Chapter 5 will investigate local stability using the E O S Q L and the Brunt-Vaisala Frequency in conjuction with C T D profiles and describe a proposed mechanism that is contributing to the deep water ventilation of Quesnel Lake.  3  2  Literature Review  This section presents a summary of other works that have addressed equations of state and local stability. In doing so, proceedures that are applied to Quesnel Lake water samples and C T D profiles are identified.  Specifically, this section addresses  other works related to equations of state (Section 2.1), the introduction of state variables (Section 2.2) and the specific issues related to each one, the developement of the relationships required to determine an equation of state (Section 2.3), and the strategies other researchers have used to investigate profile stability (Section 2.4).  2.1 Other Works related to Equations of State In the ocean, the relative proportions of dissolved salts are nearly constant throughout the world (Lewis, 1980). Therefore, only one equation of state is required to characterise stability throughout the world's oceans. In contrast, the salinity in lakes may vary from lake to lake or even between different water masses within a lake (i.e., differences in salt content between the hypo-, meta- and epilimnion) (Wuest et al., 1996). Varying local geology and river inflow can contribute to differing ion ratios (Millero, 2000a; Gray et al., 1979) which can affect the relationship between measured salinity and calculated salinity. For this reason, there is no analogous universal equation of state for all lakes. Compared to other deep lakes in the world, mixing processes in Quesnel Lake have received little scientific attention to date. With the similar motivation of understanding mixing processes, lake-specific salinity and density relationships have been developed by numerous researchers studying other deep lakes: • Lake Baikal in Siberia (Hohmann et al., 1997), at 1632m, which is the deepest lake in the world and is located in the Baikal Rift Zone. • Lake Malawi (Wuest et al., 1996), which is 700m deep and lies on the borders of Zambia, Tanzania Mozambique and Malawi, is also a rift lake and is the third  deepest lake i n the world. • Lake Issyk-kul (Vollmer et a l , 2002b), which is 668m deep and is located i n the T i e n Shan mountains of K y r g y z s t a n is the fifth deepest lake i n the world. • Crater Lake ( M c M a n u s et al., 1992) i n Oregon, which is 590m deep and is i n the collapsed caldera of M o u n t M a z a m a . It has hydrothermal vents that introduce thermally and chemically enriched water. It is the deepest lake i n the U n i t e d States and the seventh deepest lake i n the world. • Island Copper M i n e P i t (Fisher, 2002) on Vancouver Island, B C , C a n a d a , which is 340m deep. It was initially filled w i t h salt water but, due to the inflow of acid rock drainage and disassociation of ions from the pit walls, the water's chemical composition changed w i t h time such that the p r a c t i c a l salinity scale, developed for seawater, no longer applies. A time dependent salinity correction was developed to be used w i t h the equation of state for seawater. A l t h o u g h the situation for this lake is fundamentally different from Quesnel Lake, Island Copper M i n e P i t is mentioned here because similar principles were employed i n developing its equation of state. In these examples, the lake-specific equations are p r i m a r i l y determined from chemical analysis of water samples taken from the given lakes.  2.2  State Variables  A s described by E q . 1, the equation of state depends on three variables: T, S and p. T h i s section w i l l outline how density is affected by each, as well as describe some of the complications associated w i t h them. 2.2.1  Temperature  F i g . 2 shows the dependence of density on temperature. A b o v e 4 ° C , density generally increases w i t h decreasing temperature. A t zero salinity and at atmospheric pressure,  5  Temperature fC)  Figure 2: Dependence of density on temperature at atmospheric pressure and zero salinity. The vertical dashed line marks TMD- This plot is based on the seawater equation of state given in Millero et aL (1980).  density reaches a maximum of 999.972 kg-m~ at a temperature of 3.98°C. This is 3  known as the temperature of maximum density, TMD, and plays an important role in turnover in shallow (< 150 m) lakes. In shallow lakes, vertical mixing occurs when the stratification within the lake deteriorates due to temperature changes resulting from downward mixing of cooler surface water in spring and autumn. When the entire water column approaches 4°C, the stratification is reduced to the point of near uniform density. This loss of the water column stability results in seasonal turnover. The temperature of maximum density, TMD, has been a point of interest for other researchers in lake studies. In the early 1960s there was a short debate in Science between Eklund (1963, 1965) and Johnson (1964) relating to this subject. The focus was the Strom line (identified by Strom, 1945) which was incorrectly identified as TMD-  Eklund showed that TMD is actually approximately 4 ° C at the surface and has  6  a temperature gradient of — 0 . 0 2 1 ° C / b a r (Eklund, 1963) whereas the Strom line is the temperature gradient at which a profile has maximum stability (Eklund, 1965). Another point of contention was whether or not temperature profiles could cross over the temperature of maximum density. Imboden and Wuest (1995) address this topic and show that profiles may pass over the temperature of maximum density if they have a positive slope when profile temperatures are greater than TMD and a negative slope when profile temperature are equal or less than TMD, as illustrated in Fig. 3. Fig.  3 shows the change in TMD with depth along with a theoretical winter  temperature profile when reverse stratification occurs.  Although the temperature  profile crosses TMD, it is stable. If, at any place in the profile, a water parcel was moved upwards, it would be heavier than the surrounding water (i.e. closer to TMD) and would fall back to its original position. Conversely, if a water parcel was pushed downwards, it would be less dense as it would be further away from the TMD line, and would return to its original position due to buoyancy forces. This is the condition for stability. Other researchers have investigated the effects OIT D M  o lake stability problems as n  it plays an integral role in several interesting mixing mechanisms. For example, TMD plays a large role in the dynamics of thermobaric instabilities which have been studied by Sherstyankin and Kuimova (2002); Killworth et al. (1996) and Weiss et al. (1991). Thermobaric instabilities take place when wind or other external forcing pushes the cold water in the upper part of a reversely stratified water column down to where it is warmer than the local TMD- If this occurs, the overlying water becomes denser than the underlying water, resulting in episodic, local mixing. See Weiss et al. (1991) for a detailed explanation. Cabelling, another phenomena that relies on the anomaly of TMD, occurs when two water masses, one cooler than TMD and one warmer, mix together to make heavier water with a temperature closer to TMD- This heavier water then sinks and creates a plume of down-welling surface water called a thermal bar. Thermal bars occur when incoming river water or water near a lake's shore heats up while the bulk of the lake remains cold. A thermal bar can- move across a lake  7  Temperature  Figure 3: Illustration of a stable profile that crosses TMD- Based on the seawater equation of state in Millero et al. (1980).  as the season changes. This phenomenon has been studied by researchers including Shimaraev et al. (1993), Carmack (1979b) and McDougall (1987). TMD effects do not cause complications in oceanic waters. Fig. 4, which shows how seawater TMD (calculated using the equation of state for seawater extended to T = — 5.5°C for illustration purposes only) and seawater Tj  change with  reezing  pressure. At the surface, the extrapolated TMD is cooler than Tf  .  reezing  is increased, these values diverge. In waters with high (> 24.7 ° /  0 0  As pressure  ) salinity content  (like in the bulk of the ocean) water will always freeze before it reaches its TMDFor this reason, TMD effects need not be considered. However, in fresh water like in Quesnel Lake, TMD > Tf i  reez ng  making TMD effects relevant.  8  Temperature CC) Figure 4: T D M  2.2.2  and Tf  reezing  for seawater at S=35 as a function of pressure.  Salinity  Often fresh water lakes are assumed to have zero salinity. Making this assumption can lead to errors when calculating the physical and chemical properties of lakes. This was first pointed out by Chen and Millero (1977a) in relation to mistakes made in estimating Lake Ontario properties under the assumption that the water had zero salinity. The average salinity in Lake Ontario is 0.209 ° /  0 0  which, according to the  equation of state developed by Chen and Millero (1976), is enough to depress TMD by more than 0.04°C Chen and Millero (1977a). Lake water salinity, the composition of dissolved solids, cannot be found by making one measurement alone. Only a fully detailed chemical analysis on a water sample can determine its salinity (Lewis, 1980). This process is too time consuming for routine use and consequently a simpler method must be developed in order to measure salinity throughout a water column. This is done by relating salinity to easily measured C T D  9  parameters. For ocean water, C T D parameters are related to salinity using the practical salinity scale of 1978 (PSS78) (Dauphinee, 1980; Lewis, 1980; Hill et a l , 1986).  This  relationship is based on extensive measurements conducted at five seperate laboratories (Dauphinee, 1980) and is the standard equation used to determine salinity in seawater. PSS78 determines a "practical salinity" based on measurements of the conductivity ratio, r, defined in Eq. 2 as the ratio of the conductivity of the sample to the conductivity of a standard potassium chloride solution (KCI) at S=35 and T=15 °C at atmospheric pressure (Lewis, 1980):  (7(35,15,0)  V  '  PSS78 assumes that all waters with the same conductivity ratio, r, have the same salinity. The original PSS78 equation is only valid for salinities in the range of 2 to 42. Hill et al. (1986) extended the validity of PSS78 to salinities between 0 and 2 by adding a correction term. However, this extended equation still only applies to waters with the same compositions as seawater. Fig. 5 shows isopycnals of seawater with temperature and salinity changes to help illustrate density's dependence on salinity. At most temperatures the density of a lake's water is heavily dependent on the temperature of the water. However, in the temperature range near TMD, f r approaches 0. In this range, changes in density are less sensitive to changes in temperature, and salinity effects may become important. The salt content of freshwater lakes must be considered when calculating density in deep lakes such as Quesnel Lake where much of the water has temperature near TMD throughout the year (Fig. .6 a). Millero (2000b) proposes that using the equation of state for seawater at low salinities for lake water with the same salinity is adequate. He argues that the changes in the properties of low salinity solutions are independent of the specific salt added because, when diluted, the ion-ion interactions of salts in solution are reduced. He compares relative density (pi - u n  s  u  — pmean) calculated based on the composition of  10  Salinity (°/ ) Figure 5: Plot showing the dependence of density on temperature and salinity using seawater at low salinities for illustration. The horizontal dashed line marks TMDContour interval is lkg/m?. This plot is based on the equation of state.developed for the limnological range in Chen and Millero (1986) at atmospheric pressure.  Lakes Baikal (S=0.0963 % ) , Malawi (S=0.205 ° / ) , and Tanganyika (S=0.568 ° / 0  0 0  0 0  )  with the relative density for seawater diluted to the same total salinity and finds that the values agree to within the error inherent in the seawater equations.  He  concludes that lake salinity can be taken to have the same composition as seawater diluted to the same concentration of total dissolved solids as the lake water being studied.  This approach for brackish Lake Issyk-kul ( S = 6 . 0 6 ° / ) was disputed by o o  Vollmer (2002) who showed that, at this concentration, the proportion of the ions dissolved is important and so the composition of the lake must be considered. In order to use any equation of state, a water body, ocean or lake, must have a composition with a constant proportion of concentrations. As the total concentration  11  of ions change, the individual ions dissolved in the water must change at the same rate in order to maintain constant proportionality of ions and apply an equation of state. Changes in conductivity must reflect changes in both ionic and non-ionic constituent concentrations. McManus et al. (1992) discuss the issue of constant relative composition in their study of Crater Lake. The concentration of silicic acid, a non-ionic constituent introduced to Crater Lake via submerged thermal vents, was compared to the concentration of a strong electrolyte (Na ) meant to represent the lake's ionic +  salinity. McManus et al. (1992) showed that in Crater Lake silicic acid varies directly with ionic constituents and, therefore, variations in specific conductance will reflect the variations in overall concentration of dissolved solids.  The proportional varia-  tion of ionic and non-ionic constituents is assumed in studies of other lakes (Vollmer et a l , 2002b; Hohmann et al., 1997; Wuest et al., 1996). If this were not the case, the non-ionic contributors to salinity could not be parameterised in terms of easily measured quantities such as those measured by a C T D . The only alternative would be to measure non-ionic constituent concentrations directly. 2.2.3  Pressure  Pressure plays a role by mechanically increasing density due to compression. Because cold water is more compressible than warm water, TMD decreases with increasing pressure. Fig. 6b shows the relationship between pressure and TMD- In deep lakes, renewal processes may take longer than in shallow lakes because bottom water is constrained by the pressure dependence of TMD-  For example, as surface water cools  below 4 ° C , it will not mix downwards because it is less dense than water at TMD below. Cool bottom water (i.e. T < 4°C) will not mix upwards because it is closer to TMD when under pressure. In this way, interaction between surface and bottom waters can be suppressed. Small variations in the temperature profiles of Fig. 6b may or may not represent instabilities. The fluctuation in temperature may be compensated by a fluctuation in salinity resulting in a stable profile. Only with the use of an accurate equation of state  12  Figure 6: (a) Temperature profiles in East Arm of Quesnel Lake from 2001 to 2003. (b) Blow up of temperature profiles over temperature range of 3 to 4 °C.  can instabilities be identified. Once density profiles are established, quasi-density and the Brunt-Vaisala Frequency will be used to locate instabilities. In this way, mixing in Quesnel Lake can be analysed.  2.3 2.3.1  Relationships R e l a t i n g C o n d u c t i v i t y to Specific C o n d u c t a n c e  Conductivity, the ability of an aqueous solution to carry an electric charge, is dependent on the concentration of ions present, and on the temperature and pressure of the solution. Often, measured conductivity is converted to specific conductance at a reference temperature and pressure. Doing this removes the temperature and pressure signals and leaves only a conductivity signal that reflects the concentration  13  of ionic constituents. In North America, the reference temperature and pressure for specific conductance are T = 25°C and p = 0 bar, respectively (Physical and Aggregate Properties, 2000). A conversion between the measured in-situ conductivity, C£, to the specific conductivity, C ^ , can be approximated by:  Cr  — fp • fr • C ° 5  (3)  where fp is the temperature-dependent pressure correction, which accounts for the increase in conductivity due to compression of water, and fr is the temperature correction, which accounts for the dependence of conductivity on temperature for a given composition (Wuest et al., 1996). This conversion is analogous to the conductivity ratio, r, used by oceanographers (Perkin and Lewis, 1980) (See E q . 2). Bradshaw and Schleicher (1965) experimentally determined an empirical relationship between pressure and the electrical conductance of seawater. Other researchers (e.g. Wuest et al., 1996; Hohmann et al., 1997) who have studied deep lakes have extended Bradshaw and Schleicher's formula to low salinities through linear extrapolation and applied it to their respective lake waters. In addition, Hohmann et al. (1997) experimentally verified the extrapolated values. The temperature and pressure sensors of a C T D were wrapped in plastic before it was lowered to keep the same water parcel throughout the depth of the profile. The change in electrical conductivity with pressure was measured directly. The experiment yielded pressure correction values that are remarkably close to the extrapolated values for the same parameter (Hohmann et a l , 1997). In his study of Lake Issyk-Kul, Vollmer (2002) uses a pressure dependence that is developed specifically for this lake in order to eliminate small errors that are introduced by using equations developed for seawater. In the present study, for lack of a study that is more appropriate for describing pressure effects in lakes, the formula developed by Bradshaw and Schleicher (1965) extrapolated to low salinities was used to determine the pressure correction factor fp, in the equation of state for Quesnel Lake. The temperature correction depends on a lake's composition and so must be found  14  specifically for each lake. The temperature correction is determined by finding the ratio of conductivity at zero pressure, C j . , to C$ over a range of temperatures, at 5  a fixed pressure (P = P tmos heric a  P  = Odbar) and then fitting this data to a function.  Researchers have found the temperature correction both experimentally and analytically. Gray et al. (1979) and Johnson (1989) did this experimentally by measuring the conductivity of a sample of lake water over a range of temperatures in a laboratory and then fitting the collected data to a function to find the form of f . T  Analytically,  the temperature correction is determined by analysing a water sample for ionic composition, pH and alkalinity and using the collected data to calculate conductivity over a range of temperatures. This procedure was outlined in detail by Wuest et al. (1996) and used by Hohmann et al. (1997). There are a variety of functions that can be used to fit the conductivity and temperature data. Standard methods suggest the following form for this relationship:  f = % T  = \+.*{T-2S)  (4)  °25  where a is a temperature coefficient, calculated for individual water samples. Other forms that can be used are an exponential function (Gray, 1979) or, more commonly, a polynomial form (Wuest et al., 1996). The  analytical method will be applied to Quesnel Lake water samples overO to  25 ° C , the natural temperature range found in Quesnel Lake and will be described in more detail in Section 3.  2.3.2  R e l a t i n g Specific C o n d u c t a n c e t o Salinity  Rather than convert conductivity data to specific conductance, McManus et al. (1992) directly related salinity to temperature and in-situ conductivity using a least squares regression.  A more common method is to relate salinity to specific conductance  through one of two ways. For a dilute solution, a relationship between specific conductance and salinity is approximated to be linear (Sorensen and Glass, 1987) as was done for Arrow Lakes by Wuest (1999). Alternatively, if more precision is required,  15  a temperature dependent third order polynomial can be used as was done by Wuest et al. (1996) for Lake Malawi and Hohmann et al. (1997) for Lake Baikal. Quesnel Lake, which is ultra-oligotrophic, has ion concentrations that are low enough to make ion-ion interactions negligible and for this reason, salinity is assumed to be linearly related to conductivity. 2.3.3  R e l a t i n g Salinity to D e n s i t y  To calculate density in seawater, U N E S C O (1981) recommends using the equation of state developed in 1980 (EOS80) for temperatures between -4 and 40 ° C , pressures between 0 to 1000 bar and salinities between 0 to 40 ° /  0 0  . Developed separately, Chen  and Millero (1977b) initially presented an equation of state suitable for limnological systems based on temperature and pressure effects of fresh water combined with seawater effects on salinity. This limnological equation of state assumes that a lake's salt composition is identical to that of seawater, but diluted to the same concentration (Chen and Millero, 1977b). This equation was modified to have a claimed accuracy better than 2 x 10 g- c m _6  - 3  and re-presented, along with relationships for other P V T  properties such as thermal expansibility, isothermal compressibility and temperature of maximum density by Chen and Millero (1986). The equation of state presented by Chen and Millero (1986) specifically for the limnological range will hereafter be referred to as E O S C M - E O S C M , which was used by McManus et al. (1992) for Crater Lake, is valid for limnological ranges of temperatures between 0 and 30 ° C , pressures between 0 and 1700 bars and salinities between 0 and 0.6 ° /  0 0  .  E O S C M does not consider the unique composition of lakes. Millero (2000a) later suggested that one possible way to deal with the differing composition of water is to use an existing equation of state for seawater composition with a correction in density for added solutes. Reseachers including Wuest et al. (1996) and Hohmann et al. (1997) have found corrections for their respective lakes by determining an adjusted coefficient of haline contraction and using it in E O S C M -  16  The coefficient of haline contraction,  8, defined by: (5) describes the influence of salinity on density, i.e.: p  = (l Po  +  p-S)  (6)  where p is the average density. 0  Wuest et al. (1996) found the coefficient of haline contraction for Lake Malawi at p=0 bar and T = 2 5 ° C was 1.028 times greater than that found by E O S M due to the C  differing composition from dilute seawater. In calculations for density structure and stability, Wuest et al. (1996) use E O S M with B C  of  pChen+Muiero_  M a l a w i  = 1.028 •  pChen+Muiero  i  n  p  l  a  c  e  Hohmann et al. (1997) use the same method. The method outlined  by Wuest et al. (1996) will also be applied to Quesnel Lake. Density can also be measured directly. Fisher (2002) used density, measured with a densitometer, to back calculate the salinity of Island Copper Mine pit water samples. Grafe et al. (2002) showed that density can also be measured in-situ to an accuracy of pm 1 k g / m with an ultrasonic pulse transmitted through a quartz glass cylinder 3  and received by a piezo ceramic. Direct density measurements of Quesnel Lake water were done with a densitometer to verify the lake-specific relationships for salinity and density developed in this paper.  2.4  Profile Stability  Local vertical stability is described by Imboden and Wuest (1995) as a condition in which, if a water parcel is displaced isentropically by a small vertical distance, it returns to its original depth due to buoyancy forces. Because pressure plays a large role in determining in-situ density, profiles of density often increase monotonically simply because of the compressibility of water. In situ-density is not a good representation of stability because, if a water parcel experiences forcing (i.e. internal waves), it could be displaced vertically and its density would change due to pressure.. This change in density in the displaced parcel could result in a non-monotonically increasing, or  17  unstable, density profile. In order to address water column stability, in-situ density is often converted into a quantity that allows for easy comparison of densities of water parcels experiencing different pressures. Three such quantities are described below.  2.4.1 Potential Density Often, potential density, p tentiai, the density a water parcel would have it if were po  moved isentropically from its original pressure to a reference pressure, is used to remove the compression effect on density and allow comparison of water parcels within a C T D profile (Stansfield et al., 2001; Wuest et a l , 1996; McManus et a l , 1996). For potential density, the integrated adiabatic lapse rate of the parcel between the original position and reference pressure is subtracted from the in-situ density.  Ppotential(z, Z ) = p(z) - J r  <&(z,z')dz'  (7)  where z is the depth (positive upwards), z is a reference depth, and \& is the adiabatic r  lapse rate for density of a water parcel undergoing isentropic transport, defined by: . / T/  , N  R  *\. li 2) z  z  dp[ ( u 2),S(z ),p(z )} d  z  z  r  1  — —\  2  1  jisentropic-  {°J  az  2  where 6 is potential temperature, the temperature a water parcel would have it if were moved isentropically from its original pressure to a reference pressure. However, because potential density is a path-independent quantity, the second term reduces to two simple terms and the expression for p entiai simplifies to: pot  Ppotentializ, Z ) = (z) + ( (S, 0, 0) - (z)) T  2.4.2  P  P  P  = (S, 6, 0) P  (9)  Quasi-Density  Quasi-density, p si, introduce by Peeters et al. (1996), also removes the compression qua  effect on density. Peeters et al. (1996) define quasi-density as: Pquasi(z, Z ) = {z) - J " r  P  18  z')dz'  (10)  Temperature  Figure 7: (a) Illustration of important variables for potential and quasi densities, (b) Illustration of water parcel-being moved through TMD when calculating potential density. Quasi-density differs from potential density by comparing in-situ density with the integrated adiabatic lapse rate of the background water through which the parcel would pass if moved between the original pressure and a reference pressure. Fig. 7 a illustrates the difference between potential density (i.e. if parcel 1 was moved to the surface without interacting with its environment) and quasi-density (i.e. if parcel 1 was moved to the surface while interacting with its environment such as parcel 2 at Z=Z'). Considering a parcel's interaction with its environment makes quasi-density path-dependent.  A vertical path is assumed in its calculation.  Quasi-density has  no physical representation making it more difficult to conceptualise than potential density.  19  The reference depth often used in limnology and oceanography is the surface (i.e. z = 0). For potential density, a water parcel's density is determined at the surface r  which is equivalent to calculating the parcel's thermal expansion coefficient, a =  y^f,  at the surface: p = p (l + BS-aT)  (11)  0  However, because a changes sign at TMD, calculating a deep water parcel's a at the surface can cause sign errors to occur (Fig. 7b). Quasi-density takes into account the water surrounding the water parcel, and it's stability is determined locally so the sign errors, introduced in potential density, are avoided. Near TMD, potential density is not adequate for determining stability (Ekman, 1934), whereas quasi-density gives accurate vertical structure in the water column (Peeters et al., 1996). Peeters et al. (1996) argues that quasi-density, is the quantity most appropriate for the comparison of density from various depths in deep lakes that reach temperatures near TMDQuasi-density was used by Hohmann et al. (1997) for Baikal and will be applied to Quesnel Lake C T D profiles in the present work.  2.4.3  Standard Density  Foster (1995) also introduced a quantity that removes the compressibility signal from in-situ density. Standard density, D, defined in E q . 12, assumes that the measured water parcel's density's pressure dependence is linearly related to the pressure dependence of a reference water parcel. This is described by Foster's E q . 1 replicated here: D(S, T, P) = p(S, T, P)  •4 | ^ 4  (12)  p{S ,T ,P) 0  where D is the depressurised standard density, S  0  0  and T  0  are the salinity and tem-  perature of the reference water parcel respectively. Foster (1995) chose S  Q  — 35 and  T = 0°C which are close to the measured temperature and salinity of the water in his 0  study. He also suggests that standard density should be used only for small ranges of temperature and uses standard density to investigate deep water formation in the Antarctic in the temperature range of -2 to +2 ° C . 20  2.4.4  Briint-Vaisala Frequency  The Briint-Vaisala frequency, N, also known as the stability frequency, is the frequency at which a parcel will naturally oscillate if displaced from a position of neutral buoyancy. The stability, E, is proportional to N and so is important in identifying 2  where instabilities occur in the water column (Millard et al., 1990). N  2  is defined  according to E q . 13 (Imboden and Wuest, 1995).  N> = - - . / p oz 9  When N  2  .  9  > 0, the water column is stable, when N  2  (13)  < 0 the water column is unstable  (Gill, 1982). Millard et al. (1990) describes various common ways of calculating N  2  from C T D measurements and compares the results. Peeters et al. (1996) calculated the Briint-Vaisala frequency according to its definition (Imboden and Wuest, 1995) using quasi-density in the place of potential density (Eq. 14). =  r L _ . dynast  Pquasi  ^  9z  Because quasi-density provides only local stability information, the Briint Vaisala frequency calculated with quasi-density also only provides information on local stability. By taking the derivative of p i quas  with z, the following expression is found (Eq.  24 from Peeters et al. (1996)):  dpquasi _ ,dp dz ~ dz {  )isen  _ dp dz  . 1  . '  Term 1 on the right hand side is the change in density with depth, or adiabatic lapse rate of density, of a parcel of water moved a finite distance isentropically through the water column. It is compared with term 2 on the right hand side which is the adiabatic lapse rate of density of the surrounding water moved the same distance. Fig. 8 helps to illustrate the stability criteria. If the change in density of the water parcel (term 1 in Eq. 15) over the distance Az due to a change in pressure is smaller than the change in density of the surrounding water (term 2 in E q . 15), the parcel  21  p  3  (PJ)  p <p, 3  N <  Parcel Surrounding water  Figure 8: Schematic of B r i i n t V a i s a l a frequency, N . If a water parcel is moved, its change i n density due to a change i n pressure must be smaller than the change i n density of the surrounding water for the parcel to be stable and have a positive N value. 2  2  w i l l be forced back to its original position due to buoyancy forces. T h i s situation is considered stable and leads to a positive value for N . Note, 2  is a negative quantity  since density increases as depth decreases. Peeters et a l . (1996) calculated the B r i i n t V a i s a l a frequency from quasi-density values which had been averaged into 10 dbar bins. T h e y discuss the merits of choosing a moderately sized bin; a b i n size that balances resolving small instabilities, w i t h having a noisy signal for N w i t h large errors. 2  22  3  Methods  3.1  Water Samples  Between July 23rd and 25th 2003, water samples were collected at locations around Quesnel Lake in order to determine the quantities necessary to estimate the lake's equation of state. Stations were selected with an attempt to sample every limnion. Water composition throughout the lake was sampled in order to test if the proportions of salts is constant as required by the theory presented in Section 2.2.2. Samples were taken at each of the incoming rivers, Mitchell River, Horsefly River and Niagara Creek, as well as at stations in each arm that would be representative of that arm (Fig. 9). At each of these lake stations, samples were taken at the surface, at depth (near the bottom at that station) and at one or two depths in between. This was done to acquire samples representative of each strata at each station. In the East Arm, the deepest station, additional water samples (a total of 4) were taken in order to gain increased vertical resolution. Table 1 summarises the water samples taken in July, 2003.  Figure 9: Map of water sample stations, Jule 2003.  23  Table 1: Depth (m) of Quesnel Lake water samples taken in July 2003. See Fig 9 for station locations. Station 1  Station 2  Station 3  Station 4  10  10 110 160  10  50 100 125  125 350 485  85  Rivers surface samples only  Water samples were collected at 3 depths concurrently using cascading 1.71 L Niskin bottles.  From each Niskin a 250mL acid washed polyethelene bottle was  filled and the sample was acidified to a p H of less than 1.5 with Nitric acid in order to prevent metals from precipitating or absorbing into the container.  Also, a IL  polyethelyne bottle was filled as full as possible with remaining water. The IL bottles were on average 2/3 full. After being collected, water samples were sent to the Pacific Environmental Sciences Centre (PESC), an analytical chemistry lab in North Vancouver, British Columbia P E S C was requested to perform a variety of tests on the collected water samples. These test fall into two categories:  physical properties and concentration of con-  stituents, summarized in Table 2. In the development of a specific conductance relationship for in-situ C T D measurements, the conductivity of Quesnel Lake water over a range of temperatures is required. The 10 ions present in Quesnel Lake water in greatest concentration are used in this analysis. The contribution to conductivity of ions present in very small amounts (< 5% of total concentration of ionic constituents) is negligible compared to the contribution of the major ions. The results obtained from P E S C are summarized in Fig. 10 which displays the distribution of concentration of the 10 dissolved ions with greatest concentration in Quesnel Lake based on the 15 water samples collected in July 2003. In July 2003, surface water samples were also collected in 4L polyethelyne bottles  24  Table 2: Summary of analytical chemistry test performed on Quesnel Lake water samples collected in July 2003. Method  Quantity pH Alkalinity  Samples  Liquid Junction Potential electrometric titration  all samples  conductivity meter  specific conductance Anions  Ion Chromatography (ICA)  Cations  Inductively Coupled' Plasma Spectrometry (ICP)  all samples  Nitrate/Nitrite  Colorimetry  Ammonium  Colorimetry  Orthophosphate  Colorimetry  Total dissolved phosphate  Colorimetry  surface and deep  Ascorbic Acid Reduction and Colorimetry  samples in East A r m  Total phosphate  at each of the stations in Quesnel Lake.  all samples  Additional water samples were collected  on June 22nd and 23rd, 2004 for the purpose of verifying the density estimates of E O S Q £ . These water samples were taken in the East and North Arms (Fig. 11) with a 1.71L Niskin bottle primarily at a depth of 30m although some adjustments were made due to the shallower bottom depth at some stations. Table 3 summarises the water samples taken in June, 2004.  25  Na  4  Ca  Mg  CI  N0  3  S0  4  C0  3  HCO3  Figure 10: Average concentrations and standard deviations of the 10 ions with greatest concentration in Quesnel Lake in July 2003.  3.2  Calculating Conductivity  By summing the conductivities of individual ions, the conductivity for a water sample can be calculated. For ionic constituents, the sample conductivity is given by:  C  T  =  £  C  V  ( 1 6 )  i=0 where d is the conductivity of ion, i. Equivalent conductivity, A;, is the conductivity per ion equivalent at finite dilution and must be multiplied by the ion's concentration to determine the conductivity contributed by that ion. The relationship between an ion's conductivity, Cj, the ion's equivalent conductivity (at finite dilution), A,, and the ion's concentration, c», is given by:  $ =V * 26  t  ( 1 7 )  Figure 11: Map of water sample stations, June 2004.  Table 3: Depths of Quesnel Lake water samples taken in June, 2004 Station 5  Station 6  Station 7  Station 8  Station 9  Station 10  20 20  30 30  30 30  30 30  30 30  30 30  Niagara Creek surface surface  where conductivity is in units of fiS • cm^ . 1  Onsager's equation describes the concentration dependence of A^ in terms of the limiting equivalent conductivity (an ion's equivalent conductivity at infinite dilution), Xf, and sample concentration (Noggle, 1989). ^ = A°° + k •  (18)  At infinite dilution, the distance between nearest neighbouring ions is large, and the effect of the electric field due to surrounding ions is not felt. The limiting equivalent conductivity is a temperature and ion dependent quantity. Its dependence has been determined for ions that contribute significantly to Quesnel Lake's conductivity by  27  fitting empirical data (e.g. Robinson, 1959) to a third order polynomial. Fig. 12 (a) shows the fit of limiting equivalent conductivity for calcium as an example. Appendix A.2 shows the coefficients for the polynomial fit of the limiting equivalent conductivity for the top 10 ions dissolved in Quesnel Lake. Introducing a new variable /* — a reduction coefficient that takes into account the finite concentration of dissolved solids - Onsager's equation can be rewritten as follows: x  l  where the A, has units of /j,S • cm  -1  =  • eq"  xr-fi 1  (19)  and  = l +  (20)  As the distance between neighbouring ions decreases, the electric field applied on an ion due to its proximity to other charges retards the ion's motion. The reduction coefficient is a factor that quantifies the proportion by which the conducting capacity of an individual ion is reduced at finite dilution. Although the functional form of fi is unique for each ion, it depends on the equivalent concentration of the entire solution. This assumes the conducting capacity of a given ion is influenced by the presence of all other charges in the solution irrespective of its origin. It is assumed that' the proportion of ion concentrations in a lake's water is constant within that lake in both time and space. Therefore, the concentration dependence needs only to be found once and incorporated into the temperature correction which will be used in developing the equation of state that describes density throughout Quesnel Lake. Kolhrausch's Law of independent migration of ions says that, at low concentrations, the equivalent conductivity of a solution, A, is given by the sum of the equivalent conductivities of its cationic and anionic contributors: A = A+ + A_  (21)  The equivalent conductivity contributed by each ion in a salt is proportional to the ion's molecular mass. One can find the equivalent conductivity of an ion over a range  28  I  ,  .  .  ,  1  0.7'  •  •  •  0  10  20  30  40  50  0  0.1  0.2  0.3  Temperature fC)  sqrt(concentration) (mmole/L)  (a)  1/2  (b)  Figure 12: Polynomial fits for calcium of (a) limiting equivalent conductivity over a range of temperatures, (b) the reduction coefficient over a range of concentrations. of concentrations using published data for salts (Robinson, 1959; or Dean, 1979) at a reference temperature. Then, by dividing each ion's equivalent conductivities with its limiting equivalent conductivity at that temperature, one can develop a data set over which to fit a function in terms of the square root of concentration to obtain a concentration dependent form This method has been implemented using a third order polynomial for all the ions of interest in Quesnel Lake. Fig. 12 (b) shows the fit for calcium as an example. Appendix A.2 shows the coefficients for the polynomial fit of the reduction coefficient for the top 10 ions dissolved in Quesnel Lake. By substituting E q . 19 into Eq. 17, one arrives at: C = 4  ^ • Af •/<  (22)  Eq. 22 shows that the conductivity contributed by each ion is the product of the ion's equivalent concentration, its limiting equivalent conductivity and the reduction coefficient. The ion's concentration is measured directly by P E S C , and the tempertaure dependence of the limiting equivalent conductivity and concentration dependence of the reduction coefficient are determined based on published data allowing conductivity to be calculated for Quesnel Lake water over a range of temperatures. These  29  conductivities are used to find a relationship for f  T  which is described in Section 2.3.1  and used in E q . 3 to find a lake-specific function for in-situ C® in terms of C T D 5  parameters.  3.3  Finding Salinity  By assuming that the ratio of specific conductance, C% , to S is constant for a given 5  lake (Section 2.3.2), the ratio of  to S found from Quesnel Lake samples can be  related to the ratio of C ° to S measured in-situ with a C T D at Quesnel Lake in the 5  future. This is illustrated with the following equation:  S  S  1  sample °25,  CTD  1r\n\  A  2S TD  U o  m  p  (  e  C  where S  = total dissolved solids in the Quesnel Lake samples, [mg • L  sampU  ^25  s a m p  i  = e  s a m  pl ' e  s  - 1  ],  measured specific conductance, [fiS • c m ] , -1  = in-situ specific conductance calculated from C T D output, [/zS • c m ] , - 1  S  CTD  — in-situ salinity - the quantity for which an expression will be derived,  [rng.-L ] -1  A = lake-specific constant of proportionality found empirically from sample data. To use E q . 23 to find salinity in terms of C T D measurements, the empirical constant, A, is required for use with the lake-specific function for in-situ specific conductance described in Section 3.2. The value of the empirical constant for Quesnel Lake was determined from the samples collected over the summer of 2003. The salinity is calculated by summing the mass contributed by each dissolved solid in a sample. This includes acids such as carbonic acid and its equilibrium components, bicarbonate and carbonate, which  30  were not measured directly but can be calculated indirectly from measurements of pH and alkalinity (See Appendix B). Also, non-ionic constituents must be included in the calculation of salinity since, although they do not contribute to conductivity, they do contribute to density. In order to include non-ionic constituents, one must assume that the ionic and non-ionic components vary proportionately. This assumption will be further discussed in Section 4. The salinity is the sum of the mass of each dissolved constituent (including both ionic and non-ionic), and is given by: S = Y, i i m  (24)  where m; is the concentration of the dissolved constituent i and the units for salinity are mg • L  _ 1  .  The constant A is determined by finding the ratio of salinity to specific conductance. The value of A is the slope of the best fit line that passes through the origin found using a least squares fit.  3.4  Determining the Coefficient of Haline Contraction  The coefficient of haline contraction, (3, which describes how density changes with salinity was found using the procedure described by Wuest et al. (1996). Eq.  Wuest's  12(b) for B, derived from the definitions of density and the haline contraction  coefficient, is replicated here: PHViCi  B = l-~f-—  22 Mid  (25)  i  where d is the ion's concentration, M , is the ion's molar mass, and vt is the ion's specific volume, the change in volume due to an addition of a salt. Specific volume describes the expansion and contraction of water due to dissolved salts. Although specific volume is a property of individual ions, specific volume data is only available for concentrations of salts (i.e.  paired ions) (Millero et al., 1977). Therefore, ions  dissolved in Quesnel Lake must be combined into salts stoichiometrically to use the 31  available data. Since specific volume is a property of the ion, it makes no difference which salts are made as long as all ions present in the sample are accounted for through this process. The coefficient of haline contraction for Quesnel Lake was calculated based on the average composition of the 15 samples collected in July 2003. Table 4 shows how ions were combined stoichiometrically. The sum of charges, shown in Table 4, is used as a check to make sure anion and cation concentrations balance. Here, the balance is less than 3% of the total concentration of anions indicating that nearly all ions are accounted for in this method. Wuest et al. (1996) Eq. 12(b) is applied on the concentration of the assembled salts present in Quesnel Lake water to determine uesnel •  Table 4: Stoichiometric combination of ions needed for determining the coefficient of haline contraction. Ion  HCO3 0.990  Cone. [meq/L] Salts Na S0 NaCl KCI MgCOa Mg(HC0 ) Ca(HC0 ) MgSC-4 2  c o  3  0.002  3  Balance  4  Cl 0.015  Ca 0.884  Mg  Na  K  0.165  0.048  0.013  0.045 0.003  0.045  4  3  S0 0.133  0.003 0.013  0.013 0.002  0.002 2  2  0.106  0.106 0.884  0.884 0.057  0.057 0  0  0.032  0  0  0  0  0  Cations = 1.141 meq/L Anions = 1.109 meq/L Balance = 0.032 meq/L  The coefficient of haline contraction effect on density can be expressed according to Eq. 26.  p = p -(l + 0S-aT)p  = p(T)-{l  o  32  + 3-S) r  (26)  The second term in Eq. 26 adds a salinity correction to the original density. E O S C M was used as a starting point for Quesnel Lake's equation of state.  EOSCM  has the same structure as Eq. 26: PCM  where pa and  Ochen+Muiero  = P{T) • (1 + Pchen+Millero a  r  " S)  (27)  functions of temperature and pressure and S is a  e  function of temperature, pressure, and conductivity. It is assumed that BQ  i  uesne  proportionally to  varies  (Eq. 28).  Bchen+Miiiero  ^Quesnel  Once related to Pchen+MMero,  =  B  BQ  i  • Bchen+Millero  (28)  can be used in substituted into E O S C M  uense  (illustrated in Eq. 29) to get E O S Q . L  PQuesnel p (l 0  3.5  + B  • Bchen+Millero  ' S)  (29)  Measuring Density  Density of water samples taken in 2003 and 2004 were measured directly using an Anton Paar D M A 5000 densitometer. The densitometer readings were used to quantify the error present in the E O S Q L -  The Anton Paar densitometer determines the  resonant frequency of the water by oscillating the sample. The resonant frequency is inversely proportional to the root of water's mass. From the mass and the known volume (1.5mL) into which the sample is injected, the densitometer calculates the sample's density. Water samples were injected into the densitometer by syringe. Before using the syringe, the syringe's tip was rinsed with the sample water to be measured and dried. Then, the syringe body was rinsed with the sample water 4 times. The densitometer cell was also rinsed with 5.5mL of the sample. As the sample was being injected into the densitometer, the cell was watched in order to assure no air bubbles were pushed into the cell. The densitometer then automatically warmed or cooled the sample to  33  2 0 . 0 0 0 ± 0 . 0 0 1 ° C . Densities for Quesnel Lake water samples were found at atmospheric pressure. Deionized, distilled water (DDW) was run through the machine at the beginning and end of each day as well as periodically throughout the time the measurements were being made in order to assure the offset of measured values was low and that drift in values was maintained at a minimum. When drift in baseline measurements exceeded 5 ppm, the Anton Paar densitometer was cleaned using its automated cleaning process and D D W was again measured to assure baseline values within 5ppm were attained.  3.6  C T D Profiles  Conductivity, temperature and depth (CTD) profiles presented here were taken during the summers of 2001 to 2003. Each summer, a different profiler was used. Data from 2001 and 2002 was collected by a team from the Institute of Ocean Sciences made up of Dr. Fiona McLaughlin, Dr. Mike Foreman, and John Morrison. In July 2003, C T D profiles were taken at each of the stations where water samples were collected (Fig. 9), as well as one in the junction (Fig. 1) using a Sea-Bird Electronics Inc. S B E 19 C T D profiler. The S B E 19 used in 2003, had a sample period of 0.5s, and a resolution of ± 0 . 0 0 0 8 ° C in temperature, ± 0 . 1 7 5 pS/cm in conductivity and ± 0 . 0 9 dbar in pressure. The C T D was used in narrow conductivity mode to work in low salinity water. Only the C T D down-casts were used in quasi-density and N  2  calculations.  3.7  Applying Quesnel Lake's Equation of State  3.7.1  Quasi-Density  Analysis of vertical mixing requires that in-situ density profiles, which are dominated by hydrostatic pressure, be converted into a quantity that allows water parcels to be compared within a profile. Quasi-density profiles were found from C T D profiles, based on E O S Q L , using the procedure described by Peeters et al. (1996).  34  • C T D data were sorted so that only the downcast was used and pressure was always increasing. Data collected in error was removed from the data set. Where two or more measurements were taken at the same depth, data was averaged so that each pressure had a unique set of C and T measurements. Quasi-density was then calculated from the following algorithm which was executed for each measurement (i) in the C T D cast. Each calculation of quasi-density requires 3 data points (i —  + 1).  Line 1: ft = 6(Si- T _ P^, PA t  U  u  Potential temperature was determined over the whole water column using data point, i — 1, the initial pressure and the measurement directly above it, i, as the reference pressure.  Line 2: a  =EOS (S - ,Ti- ,P . )  u  QL  i  1  1  i 1  In-situ density was determined using E O S Q L for each measurement, i — 1, in the C T D cast.  Line 3:  a  12  =EOS (S . ,6 ,PA QL  i l  1  Density was also determined using E O S Q L for the measurement directly above i — 1, at i, as though the parcel at i — 1 had been lifted isentropically (without change in heat or salt) to the pressure at i. For this reason, the density is found using the salinity at i — 1, and the potential temperature, 6, as if the parcel was lifted to Pi with no heat loss.  Line 4:  (j^-i =  In order to find the local adiabatic lapse rate for density (defined in Eq. 8), the change in density with depth was found for two adjacent depth intervals and averaged together to find the adiabatic lapse rate at i.  Line 5:  6 = 6(3,, T P P ) 2  it  u  i+1  35  Line 6:  a  21  —EOSQ (Si,Ti, Pi)  Line 7:  a  22  =EOS {Si, 9 ,  Line 8-  C^V =^ " ^ i  L  QL  P ) i+1  2  Lines 6 through 9 are a repetition of line 2 through 5 with measurements at i and i + 1 to find the adiabatic lapse rate for density at the second depth.  Line 9:  integrand(i)  = ^ {  )  1  + { 2  ^  ) i  In each iteration, one value for ^dz(i)  • *±IZ£«=I is determined and added onto the  end of the vector. In this way, computation time is saved because throughout the whole water column is not calculated in each iteration. Line 10:  X =  ^{integrand)  All entries in the integrand vector were summed to integrate the adiabatic lapse rate of density using the trapezoidal rule for integration (See Press et al., 1987).  Line 11:  p u  =EOS (Si, T , Pi) - X  quas  Finally, p i quas  w  QL  a  s  t  calculated by subtracting the integrated adiabatic  lapse rate from the in-situ density. This was repeated until p i quas  was  . calculated for each point throughout the profile. 3.7.2  Standard Density  Throughout this paper, standard density is calculated according to the definition presented in the paper written by Foster (1995) and replicated here in E q . 12. 3.7.3  B r i i n t V a i s a l a Frequency  The Briint-Vaisala frequency is used to investigate water column stability from C T D profiles. It is calculated according to its definition (Imboden and Wuest, 1995) using quasi-density in place of potential density (Eq. 14) as recommended by Peeters et al.  36  (1996).  Quasi-density was calculated using E O S Q L and the algorith presented in  Section 3.7.1. Without bin averaging C T D data, the N profile contains a lot of information and 2  is very noisy. In order to remove the noise, quasi-density data was smoothed through bin averaging. Bin averaging was performed before N was calculated because, oth2  erwise, the quasi-density gradient found from data points taken closer together (ie with minimal spacial resolution) are disproportionately large. The size of bin chosen for averaging is important. When averaging is done, information about N is lost, 2  and, simultaneously, the error in N is reduced. The bin size must be chosen in a way 2  that balances resolving structure of interest in the profile and maintaining a small enough error so that the signal that is resolved is meaningful. This compromise is also discussed by Peeters et al. (1996). The error introduced through bin averaged  p si qua  is estimated by finding the sum  of the squares of the error in each bin as follows; Jv  1 ^Pquasi  —  !j~ ^ .(.Pquasi;  Pquasi)  (30)  i  where  N  is the number of values of  the bin, and  p  all values of  p  qua  is the mean value of p  i  quaS  q u a s  in the bin,  p si  i  q u a s  P q u a a  i.  is each value of  p  q u a s  i  .in  in the bin. This approach assumes that  i  are without error and the error introduced is done so because the  information contained in the bin averaged data is lost. After N is calculated, its error, 5N , is estimated according to 2  2  6N = 2  9—- • V2 • 6p  (31)  quasi  Pquasi  ' ^  z  This approach assumes that the error contributed to S N from Az is negligible com2  pared to the error contributed by  Ap . quasi  37  4  Quesnel Lake's Equation of State  4.1  E q u a t i o n for f  T  In order to have an equation that describes the whole lake, one must assume that the conductivity reflects the total concentration of dissolved solids. Non-ionic constituents, such as silicic acid, must be factored into this equation. This is done by assuming concentrations of non-ionic constituents is proportional to concentrations of ionic-constituents.  To test this assumption for Quesnel Lake, the concentration  of silicon (Si) and a strong electrolyte present in Quesnel Lake, calcium (Ca) are compared (Fig. 13). 18.4  -i  r  -i  1  1  1  1  L_  T~  18.2  18h  E  17.8  5 a  17.6 \  8  17h  O •fi 17.4 h  c  o  O  16.8 16.6 16.4  0.8  _l  0.9  1.1  1.2  I  I  I  1.3  1.4  1.5  1.7  1.8  Concentration of Silicon (mg/L)  Figure 13: Concentration of silicon (Si) versus calcium (Ca) for 15 samples taken from Quesnel Lake July, 2003.  If Si and C a varied proportionately, the data would fall on a line that passed through the origin. However, some variability is present in this data. The alternative to making this assumption would be to make direct measurements of non-ionic constituents with each C T D cast. This is prohibitively expensive and time consuming.  38  Temperature f C) Figure 14: Temperature factor f functions of interest.  T  function for Quesnel Lake compared to other  f  T  The assumption that the conductivity readings reflect all dissolved solids (both ionic and non-ionic) is one source of error in the equation of state developed specifically for Quesnel Lake. This error will be assessed later in this section. The relationship for the temperature factor, fx, was determined according to the method described in Section 2.3.1: fx = 0.5369 + 0.0156 • T + 1.339xl0  -4  •T  2  - 7.552xl0~ • T 7  3  (32)  Out of interest, this curve is compared to other curves generated for different lakes (Fig. 14). Notice the seawater and freshwater curves which bracket the Quesnel Lake  39  Measured Specific Conductance (u, S/ cm) Figure 15: Measured and calculated specific conductance for 15 samples taken from Quesnel Lake July, 2003.  curve. Notice also that the Arrow Lake (which have a similar composition to Quesnel Lake) has a very similar fx curve to Quesnel Lake. To verify that the fx curve for Quesnel Lake can accurately predict specific conductance values, calculated specific conductances are compared to specific conductance measured in the lab (Fig. 15). This data does not lie perfectly on the one to one line. However, the R value for this data is R = 0.85 indicating that more than 50%.of 2  2  the data is described by the 1 to 1 line. Error introduce by making this assumption will be discussed further in Section 4.3.3. Quesnel Lake data were compared to calculated specific conductances from other B.C. lakes including Arrow Lake, Slocan Lake and Kootenay Lake, for comparison over a larger range of specific conductances (Fig. 16). The plot shows that the Quesnel Lake data falls in a cluster over a small range of specific conductance compared to the  40  Figure 16: Comparison of Measured and Calculated Specific Conductance for A r row, Slocan, Kootenay and Quesnel Lakes. Arrow, Slocan and Kootenay Lakes data courtesy Roger Pieters unless otherwise marked.  range of specific conductances found in nearby lakes. Although, Fig. 15 shows that the calculated specific conductances vary from the measured specific conductances, Fig. 16 suggests that this variability may be a symptom of the samples having been taken over a small range of salinity. Over a larger range of salinity, this variability appears less drastic. A further discussion of estimates of error can be found in Section 4.3.3.  4.2 4.2.1  Equation for Salinity Finding the Constant of Proportionality, A  Salinity and measured specific conductance were plotted for all 15 samples collected from Quesnel Lake in July, 2003 (Fig. 17). A is the slope of the best fit line which  41  100  105  110 115 120 Measured Specific Conductance (u,S/cm)  125  130  Figure 17: Measured specific conductance plotted as a function of salinity in order to determine the value for A for Quesnel Lake. Data from all 15 Quesnel Lake samples are plotted with a solid line showing the best fit value for A .  passes through the origin. The value of A was determined empirically to be: A = 0.836m# • cm • L"  • fiS'  1  (33)  1  Fig. 17 shows some variability in the value of A determined from individual samples. Again, this variability could cause error that will be discussed further in Section 4.3.3.  4.2.2  Salinity Equation  Eq. 3 can be substituted into Eq. 23 and rewritten as: Q  CTD  U  — A .  f  (34)  f  IT - JF  In this way, once A (Eq. 33) and the functional form of f  T  (Eq. 32) are determined  for a specific lake, the salinity can be found throughout the water column with C T D  42  measurements.  This was done for the samples collected from Quesnel Lake in the  summer of 2003 to give a lake specific expression for salinity:  C T D  ~ (0.5369 + 0.0156 • T + 1.339X10" • T  2  4  - 7.552a;10- • T ) • f (P, 7  3  P  T)  K  1  where the pressure function taken from Wuest et al. (1996) is given by: /P(P,T) = 1+ - P  (36)  7  and 7 = 10~ • (1.856 - 0.05601 • T + 0.0007 • T ) 5  2  In E q . 35, C , T , and P are all measured with a C T D .  (37)  After a relationship for  salinity is determined in terms of easily measurable quantities, E O S Q L can be applied to calculate density as described in the next section.  4.3 4.3.1 The  E q u a t i o n o f S t a t e for Q u e s n e l L a k e Coefficient of H a l i n e  Contraction  analysis of the chemical constituents in Quesnel Lake water samples described  in Section 3.4 predicts the value for the coefficient of haline contraction, 3, to be: ^Quesnel = 0.827mg • I T  1  (38)  This Quesnel Lake specific coefficient of haline contraction can now be integrated into E O S C M to obtain EOSQL(1996) is used. BQ  i  To do this, the procedure described by Wuest et al.  is related to the Chen and Millero coefficient of haline con-  uesne  traction, Pchen+Millero0Quesnel = 1-096 • Pchen+Millero  (39)  and this new coefficient of haline contraction will be used to modify E O S C M according to Eq. 40.  PQuesnel Po (1 + 1.096-Pchen+Millero-S) 43  (40)  7  (a)  —==•  (c)  /• t  100  100  100  200  200  300  300  400  400  /' i  /  =§ 200  •c • 1 1  ;1  3 UJ  1  :i  £ 300  • t  •/'  Q_  .7 .7  •/  400  1 500 .' 80 3  s  T T  MD  ( j I  1 ' '. "  85 Salinity (mg/L) 4  P -P =oS  500 90 999 1000 1001 1002 „3> Density (Kg/m J  S  500 0.069 0.07 Density difference (Kg/m ) 3  5  Temperature f C) Figure 18: Comparing density with and without Quesnel Lake salinity contribution: (a) T and S profiles from the deep basin of Quesnel Lake, 2001. (b) Profiles of density, with and without the contribution from S. (c) Difference betwen density estimated using S calculated using E q . 35 and density estimated with S=0.  Fig.  18a shows T and S profiles from the deep basin of Quesnel Lake.  Below  approximately 100m, the temperature roughly follows TMD- At temperatures near TMD, the thermal expansion coefficient is near-zero, and so small salinity effects play an important role in determining water column stability. To illustrate this, Fig. 18b shows the corresponding profiles of density, both with and without the contribution from salinity and Fig. 18c shows the difference between these two curves. This difference highlights the importance of salinity's contribution to density in understanding' vertical mixing processes in the deep waters of Quesnel Lake.  44  4.3.2  Quesnel Lake's E q u a t i o n of State  E O S C M is modified for Quesnel Lake salinity according to Eq. 29 to produce E O S Q / , : PQuesnel = 0.9998395 + 6.7914 • 10~ • T - 9.0894 • 10" • T 5  +1.0171 • 10~ • T 7  +1.592 • 1 0 ~  U  •T  b  6  3  2  - 1.2846 • 10~ • T  - 5.0125 • 1 0 "  9  14  •T  +1.096 • (8.181 • 10~ - 3.85 • 10~ • T + 4.96 • 10~ • T ) 4  6  8  2  •, •  4  •• •  6  •• • • S.  (41)  where S is determined from C T D data using a lake-specific equation (eq. 35). 4.3.3  Error  Fifteen samples taken throughout Quesnel Lake have been analysed to find fr. Since this function is dependent on composition and since they have been taken from different parts of the lake, the question remains: to what accuracy can density be estimated from C T D data using a single equation of state? Error is introduced into .the developed E O S Q L by many avenues.  They include  the following: • Error in the original measurements by P E S C of concentrations of the major ions dissolved in Quesnel Lake; • Error is introduced by using marine equations to find concentrations of components of the carbonate system for Quesnel Lake water; • Error is introduced in the fitting process of polynomials to limiting equivalent conductivity and reduction coefficients data; • E O S Q X has been found based on the composition of 15 samples taken from various places around Quesnel Lake, because some variability exists between these samples, E O S Q L , which is based on the average, will not be specific to different lake regions;  45  • Variability exists in the plot of Si vs. C a (Fig. 13) showing that the assumption that they are regularly proportional is inaccurate, and; • The charge balance for the average composition of the 15 samples is not zero, which may have an impact on the coefficient of haline contraction, 3, determined in E O S  Q L  .  To assess error in E O S Q L , calculated densities were compared to densities measured directly with an Antor Paar densitometer. To calibrate the densitometer, the density of D D W was measured. Table 5 shows the expected and measured results of these measurements.  Based on the root mean square of the differences between  the expected and measured values, an offset of 2.6 x 1 0  - 6  g/cm  3  was assumed in  subsequent density measurements.  Table 5: Measured Density values of distilled, deionized water (DDW) used to calibrate the Anton Paar Densitometer , Density of D D W Measured  Expected  g/cm  g/cm  3  3  0.998204  0.998203  0.998205 0.998205 0.998206 0.998207  The  results of density measurements of all Quesnel Lake water samples collected  in July 2003 and June 2004 are shown in Figures 19 and 20. Fig.  19 shows measured densities and densities calculated with E O S Q L along with  a 1:1 line. The outlier in the upper right hand corner of the plot was taken at the Niagara Creek inflow to Quesnel Lake's East Arm.  This water had high silt content  at the time of collection which may account for its uncharacteristically high specific conductivity.  The linear correlation coefficient of the data in Fig.  46  19, excluding  the outlier, is B? = 0.929 indicating that 93% of the variation in density can be explained with the 1:1 line.  This correlation coefficient is close to 1 suggesting a  linear relationship between measured and calculated densities. Also, to assess the need of developing a lake specific equation of state, densities calculated E O S C M were compared to densities calculated with E O S Q L and densities measured with the densitometer (Fig. 20). Although the densities calculated using E O S Q L are often on the low end of the measured values, they are within the estimated error of the measured values. Also, the density values calculated with E O S Q £ are closer to the measured valued than those calculated with E O S C M - The standard deviation between measured values of density and those calculated with EOSQZ, is 0.0018 k g / m . The standard deviation between measured values of density and those 3  calculated with E O S C M is 0.0157 k g / m  3  which is signficantly larger than that cal-  culated with EOS<2£. E O S C M describes Quesnel Lake water density less accurately than reported by Chen and Millero (1986). The outlier samples (samples 7 and 8 in Fig.  20) are again apparent. Neither E O S Q L nor E O S C M accurately predict density  for water in that area of the lake.  47  -1.65 -1.66  x measured vs. calculated - - 1:1 line R =0.929 (ommitting outlier) 2  -1,67 -1.68 -1.69 n)  - 1 7  -1.71 -1.72 -1.73 -1.73  X  -1.725  X  -1.72 -1.715 Calculated o (kg/m  -1.71  -1.705  3  Figure 19: Measured and calculated densities for water samples collected from Quesnel Lake in July 2003 and June 2004. The outlier in the upper right hand corner of the plot was taken at the Niagara Creek inflow to Quesnel Lake's East A r m .  48  IT  -1.65 -1.66  * o x  Chen and Millero EOS Quesnel Lake EOS Measured  -1.67 -1.68 \  -1.69 ~  -1.7 o  o  -1.71 -1.72 -1.73 -1.74  11  * * **  * 8  *  * *• 5 J 5 ?  10 12 Sample Number  * * ** 14  16  18  20  22  Figure 20: Comparison of densities at T=20 °C and P=0 bar in Quesnel Lake: density measured directly with the Anton Paar densitometer; density calculated with E O S Q L , and density calculated with E O S C M -  The outlier samples (samples 7 and 8) were  taken at Niagara Creek inflow to Quesnel Lake's East A r m . Neither E O S Q I , nor E O S C M accurately predict density for water in that area of the lake.  49  5  Local Stability  Developing EOSQL  was done in the interest of investigating the deep water ventilation  in Quesnel Lake. A preliminary look at vertical mixing gives an indication of where and how this ventilation is taking place. For this reason, profiles of quasi-density and the Briint-Vaisala frequency have been found based on C T D profiles.  Figure 21: A comparison of potential and quasi densities for C T D profile taken in 2001: a) plot of in-situ density, b) potential and quasi-density, c) zoom-in of potential and quasi-density.  5.1  Quasi-Density  Quasi-density profiles were found according the procedure explained in Section 3.7.1. Fig. 21 illustrates the importance of using quasi-density for investigating local stability. Fig. 21a shows a plot of in-situ density. The hydrostatic pressure is the dominant  50  signal in this density plot. In order to get a plot of a quantity that can be used to assess column stability and compare densities from within a profile, this pressure signal must be removed. Both potential and quasi-density serve this purpose and plots of both for a C T D cast taken in the East A r m in 2001 are shown in Fig. 21b. Both of these return a virtually vertical profile - one with no hydrostatic pressure signal. Once zoomed in (Fig. 21c), one can see that the potential and quasi-density plots are not perfectly vertical. The potential density plot shows that below the thermocline, density gradually decreases thereby predicting a mostly unstable profile. Quasi-density predicts a profile which contains some instability but has nearly monotonically increasing densities and is stable.  Quasi-density was shown to be more appropriate  than potential density by arguments made in Section 2.4.2 and is supported by the validity of the pressure-free profiles predicted by each. In Quesnel Lake, quasi-density should be used to study local stability rather than potential density. Quasi-density was found from a C T D profile using E O S Q L and E O S C M to calculate in-situ density. Fig.22 shows this comparison. Quasi-density calculated with  EOSQL  is different than that calculated with E O S C M - However, the difference is primarily a simple offset in values; the structure of the profiles are still very similar. This will be explored further in Section 5.3 where the Briint Vaisala Frequency is discussed. Fig. Arm  23 displays quasi-density profiles from the C T D data collected in the East  over three years (2001, 2002 and 2003) in order to gain insight into how the  vertical structure in the water column varies over a period of time.  Quasi-density  profiles in the East A r m over the three year period all demonstrate very similar characteristics. Profiles have generally increasing slopes with density flucturations of simlar magnitude. Fig.  24 displays quasi-density profiles from the C T D data collected in July 2003  at each of the stations presented in Fig. 9 as well at the junction in order to gain insight on how vertical structure varies spatially throughout the lake. Quasi-density profiles from 2003, taken around Quesnel Lake, show different characteristics. The profile taken at the junction echoes characteristics of profiles taken in the North and  51  Figure 22: E O S M and E O S for C T D profile taken in 2001. Although density values are different, the water column stability is roughly the same in each profile. C  Q L  West Arms. For example, three small spikes in the West A r m profile at 65, 145 and 200 dbar are also present in the Junction profile at the same depths. A small step in quasi-density at 105 dbar in the North A r m is also present in the Junction profile at the same depth.  52  53  400'  ' 0.044  : 0.046 1  1  1  0.048  0.05  1  1  0.052  '  1  0.054  0.056  Standard Density, D (kg/m ) 3  Figure 25: Standard density calculated by varying the temperature, T , of the reference parcel to temperatures bracketing TMD in Quesnel Lake (3°C to 4 ° C ) . 0  5.2  Standard Density  The value calculated for standard density is very sensitive to which temperature is chosen as the temperature of the reference parcel and less sensitive to the reference parcel's salinity. Fig. 25 shows the standard density profiles calculated for the C T D cast collected at Quesnel Lake in 2001 with varying temperatures which bracket the variation in deep temepratures in Quesnel Lake chosen as the reference. Fig. 25 illustrates that the choice of reference temperature is critical to calculating a quantity that is meaningful for finding local stability in the water column. Standard density calculated from the 2001 Quesnel Lake C T D casts with T = 3 ° C shows an 0  almost entirely unstable profile while the same C T D cast calculated with T  A  shows an almost entirely stable profile.  55  = 4°C  1.002 Density using QL-EOS, S=0 (g/cm ) 3  Figure 26: Density as a function of Temperature and Pressure (2D). The inset shows a blow up of the indicated region with bolded line illustrating the difference in slope between density gradients at 0 ° C and 7°C.  The reason for this difference lies in the assumption that the change in density due to pressure at one temperature is equal to the change in density due to pressure at another temperature. Figures 26 helps to illustrate this point. Fig. 26 is a plot of density as a function of pressure with multiple lines showing this dependence for varying temperatures. The bold lines highlight the density gradients at two temperatures, one higher than TMD at 7°C, and one lower at 0°C. These lines are not parallel. The gradients of density with pressure are different at adjacent temperatures. This make it inappropriate to try to predict the change in density with pressure based on a parcel at a reference temperature different than the in-situ temperature with water neat TMDSignificant error would be introduced if standard density were calculated for a profile near 0°C with a reference parcel assumed to be at 7 ° C . This example illus-  56  trates an extreme case and gradients do become more similar at closer temperatures. However, this difference in gradients is large enough near TMD to cause errors when using standard density to assess stability. Standard density performed well for the work Foster (1995) did in the Antarctic. He suggests that standard density is only valid for small temperature ranges. More accurately, standard density is only valid for small ranges of temperature away from TMD where the changes of density with pressure are more uniform for neighbouring temperatures.  Because Foster (1995) is applying the standard density concept in  salt water, water that will freeze before approaching its temperature of maximum density, the anomalies which occur near TMD in fresh water are inconsequential to his work. However, standard density cannot be applied to the investigation of deep water ventilation of Quesnel Lake water, whose water temperature mostly lies near TMD-  5.3  Instability Identification  Quasi-density is used to calculate the Briint Vaisala frequency in order to identify where density inversions occur. Quasi-density is noisy when calculated directly from C T D profiler data. Consequently, the Briint Vaisala frequency calculated from quasidensity is also noisy — too noisy to identify instabilities. To filter out noise and to gain a meaningful signal, the quasi-density profile is smoothed through bin averaging. The total change in density over the instability is less than the fluctuations in unfiltered quasi-density, illustrating the need to bin average. By identifying the trends in the quasi-density data over a resolution of interest, instabilities of that resolution size are identified. Fig. 27 shows a blown up section of the quasi-density profile with an inset of the entire profile. The blown up section of the profile also features the profile after bin averaging has been performed with various bin sizes. A power spectrum for this quasi-density profile is shown in Fig. 28 on a log-log scale. Because so many frequencies are present in the quasi-density profile in Fig. 27, the power spectrum plot was generated in order to determine if there is some obvious  57  250 0.056 — 0.0585  i  0.059 0.0595 Quasi Density, p '  0.06  no bin averaging binsize=2 dbar binsize=8 dbar binsize=20 dbar  J  0.058  0.058  0.06  0.0605  "quasi  Figure 27: Quasi-density of C T D profile taken in 2003 in Quesnel Lake's Junction bin averaged to remove noise from signal into bin of sizes 2 dbar, 8 dbar and 20 dbar. The insert shows the entire profile.  feature that could impartially help select a cutoff frequency for which signals with greater frequency could be considered noise. Fig. 28 displays frequency power data with a relatively constant negative slope and no such feature is present. Briint Vaisala frequency plots are examined with various frequencies filtered out using bin averaging. The Briint Vaisala frequency, N , and its error have been found 2  for three bin sizes: 2, 8 and 20 dbar (Fig. 29). The N  2  profiles where found using  E O S Q L - For progressively smaller bin sizes N becomes less noisy and its error also 2  becomes smaller. The mean error in N dbar and 20 dbar are 5.5 x 1 0  - 7  2  for profiles found with bin sizes 2 dbar, 8  s , 1.6 x 1 0 -2  - 7  s~ , and 7.4 x 1 0 2  - 8  s  - 2  respectively  (see Section 3.7.3) . However, vertical resolution is lost in this process. It is crucial to know what size instability is of interest to know what size bin to choose. Choosing  58  io-  3  Wave Number, k ( m ' ) Figure 28: Power Spectrum for the unfiltered quasi-density profile taken in 2003 at the Junction in Quesnel Lake.  too large a bin may reduce resolution to a point where instabilities of interest are no longer visible as was done for this profile with a bin size of 20 dbar. For the purpose of this study, a bin size of 8 dbar was used to compare N  2  profiles calculated from  C T D profiles taken in Quesnel Lake's East A r m in 2001, 2002, and 2003. Also, N  2  profiles taken in 2003 from Quesnel Lake's three Arms, Junction and West Basin are also compared. From this series of profiles, one is able to deduce the resolution at which water column is completely stable. For example, the centre plot in Fig. 29 shows that when looking at a resolution of 8 dbar, there are some places in the profile that are slightly unstable. However, at a resolution of 20 dbar, the entire water column appears to be neutrally stable.  59  LO  f  o  k< .Q  T3 o  OJ II d> N  '(/J C  Figure 29: Briint Vaisala frequency of C T D profile taken in 2003 in Quesnel Lake's Junction found using Quasi-density averaged in to bins of sizes 2 dbar, 8 dbar and 20 dbar. The error in the Briint Vaisala frequency is also plotted for each of these bin sizes.  60  N (p  .(QL-EOS))-N (p  2  2  .(CM-EOS))  Figure 30: a)Briint Vaisala Frequency found using E O S C M and E O S Q / , for C T D profile taken in 2001, b) the difference between N calculated with / ^ u a s j ( E O S Q ^ ) and 2  Pquasi ( E O S C M ) .  Fig.  30 shows N  2  calculated for the profile taken in 2001 (for comparison with  profiles featured in Fig. 22) with quasi-density using both E O S Q X and E O S C M for comparison as well as the difference in the two quantities.  The two equations of  state predict very similar profiles. The difference indicates that E O S Q L shows greater stratification through the water column and so estimates that instabilities are more unstable and that stabilities are more stable than what E O S C M estimates.  How-  ever, the structure (i.e. gradient) predicted by the two equations of state are nearly indistinguishable. Fig. 31 shows N calculated with both potential and quasi-density for comparison 2  based on densities calculated with EOSQ/_,.  In Fig. 21, potential density portrays a  profile with an unstable region (i.e. positive slope) below approximately 250 m. This is echoed in N profile of the potential density profile. 2  61  0 50 100  bar)  150  CD CO W CD CL  200  5=  *  250 300 350 400 450 500 -5  4  5  N (s- ) 2  2  x 10~  7  Figure 31: Briint Vaisala Frequency found using taken in 2001.  62  p t tiai po  en  and.p  i-  quaS  for C T D profile •  Fig. 32 displays N profiles from the C T D data collected in the East A r m over 2  three years (2001, 2002 and 2003) in order to gain some perspective how the vertical structure in the water column varies over a period of time. Fig. 33 displays N profiles 2  from the C T D data collected at each of the stations presented in Fig. 9 as well at the junction in order to gain some perspective on how vertical structure varies spatially throughout the lake. Again, N profiles in the East Arm.over the three year period 2  all demonstrate very similar characteristics. Profiles indicate that the majority of the water column below approximately 50 m is close to being neutrally stable, save some locations in 2001 and 2002 where the profile becomes slightly unstable. N  2  profiles  from 2001 and 2002 both exhibit these unstable regions near 400m depth. Quasi-density profiles taken in 2003, from the different areas around Quesnel Lake, show differing characteristics.  The East Arm profile has N  2  ~ 0 s  - 2  below 50 m.  Compared to other areas in Quesnel Lake, the East Arm is relatively neutrally stable. This may be an indicator as to how deep water ventilation is taking place in this region of the lake.  63  to hrj O O CTQ H  - 1  to O o to  Pi l-S  « CO to • •  0  O- O  to o o  B  2002  2001  2003 0  CO  00  50  -  _ J ~  CD  5' Cd  \  J  i CD  >-i P=  00  100 150  <  200  CD  CL P;  fj  o00 P-  '> •-  CD  CD  >  O  r*  *  -  \  250  -  300  -  --'  \  <^  350  V  400  CD  ex o  450  O H O N 2  (° - ) s  2  o B  CD"  CO  cT 3  X10-  7  N  2° < -2) S  5  x10-  7  -  -  500 -5  N  2  (s -2'  x10  E <  I  1  '  -  .'  •  \  •  OJ <? o  r -  •c  o z  1  -  ^—|  '  l  ,  -r-  T-  1  o X  -  (/> 03 *  O LO  O O  4  O LO i-  O  O OJ  O  LO OJ  o o  CO  o  LO CO  O O  I  (jeqp) sjnssejd Figure 33: Profiles of the Briint Vaisala Frequency based on C T D profiles taken in 2003 at stations in the East A r m , North A r m , Junction, West A r m , and West Basin bin-averaged to 8 dbar.  65  5.4  Density Currents  The formation of a density current at the m o u t h of the N i a g a r a Creek is one mechnism that could contribute to the deep water ventilation i n Quesnel Lake. N i a g a r a Creek water originates i n a glaciated catchment and is heavily loaded w i t h silt. T h i s silt may increase the density of incoming river water enough to overcome the barrier caused by the temperature of m a x i m u m density previously described. To illustrate the mechanism, F i g . 34 shows a schematic of the density current. In F i g . 34a, where incoming riverine water is near the deep (P « 500 dbar) T D M  (2.93  ° C ) , cold, silt loaded river water passes through surface waters (which, i n the absence of silt, would be lighter than water at the surface TMD), and, v i a the density current, is driven to the b o t t o m of the lake. Entrainment of the surrounding water would increase the volume of the density current as it plunges.  T h i s , together w i t h the  gradually smaller density difference between the density current a n d its surroundings as the denisty current reaches greater depths, would cause the density current to lose momentum. T h e density current water could potentially reach the lake b o t t o m and ventilate Quesnel Lake's deep water. In F i g . 34b, where incoming riverine water is near the surface TMD, river water is again driven to the lake b o t t o m v i a the density current, but as the water loses momentum, silt w i l l settle out of the river water.  T h e clarified river water then  becomes buoyant a n d begins to move upward to a depth where the it has the same temperature as the surrounding water's TMD and is neutrally stable. In the latter case, b o t t o m water w i l l be entrained and the mixed river/lake water w i l l move upward together contributing to the flux of nutrients to the surface. Due to seasonal variation i n river temperature a n d flow, w h i c h entrains varying amounts of silt, the density difference between incoming river water a n d lake water w i l l vary. Because of this variation i n the density current's density, as it loses m o m e n t u m , varying degrees of silt w i l l be unloaded before the water becomes buoyant. Together, these factors lead to cases where the process described above is modified. Above, two  66  idealised cases are described, whereas in reality, Niagara Creek plunges to varying depth, and surfaces with varying clarity. The re-surfacing density current waters, which have been cycled up from deep parts of the lake, will be cold and have some suspended particulates. Evidence for the density driven mechanism comes from direct density measurements on samples collected from Niagara Creek on June 23, 2004. These density measurements show that Niagara Creek water density is significanly larger than rest of lake.  At atmospheric pressure and a temperature of 20.000 ± 0.001 ° C ,  Pmagara = 998.3480 ± 0.0026%/m and p nei 3  QueS  = 998.2818 ± 0.0026%/m . 3  Nia-  gara Creek water would remain denser than Quesnel Lake water through mechanical compression as it plunges.  Samples 7 and 8 in Fig. 20 were collected within the  turbid plume at the Niagara Creek inflow. Samples 9 and 10 were collected just outside of that plume. Density calculated using on E O S C M and E O S Q L based on water sample conductivity and temperature greatly underestimate the measured density. The measured density includes the density from non-ionic constituents in the sample water which contribute significantly in Niagara Creek water density. This additional non-ionic material is likely the silt mentioned above. The N profiles found for the East Arm in 2001, 2002, and 2003 (Fig. 32) indicate 2  that the whole water column is close to neutrally stable. This may result from the homogenizing effect the plume generates as the density current loses silt and convects upwards to its depth of neutral stability.  Density current water then remains at  a depth where it is neutrally buoyant until it experiences any subsequent external forcing. Indirect evidence also supports this mechanism. Visual observations suggest that water in the East Arm has lower transparency compared to the rest of the lake (personal communication from Daniel Potts, June, 2004; also observed by Christina James in July 2003, and September 2004). This is likely due to unsettled silt moved to the surface by the buoyant density current water. Note, the surface water in the East Arm was found to be colder than surface  67  temperatures in the rest of the lake in the summer months in two studies. In the fall of 1981, surface temperatures in the lake's East Arm were 4.5 °C colder than the West Basin (Stockner and Costella, 1980). In 1982, surface temperatures in the East Arm were again observed to be colder than the junction by 8.5 °C in June and by 4.5 °C in August (Morton and Williams, 1990).  68  Ventilation of deep water  (a) Ventilation o f surface water  T= 4 °C  -  Settling out o f sedimentation  (b) Figure 34: Illustration of density current: (a) water with T ~ 3°C sinks through surface water due to addition of suspended sediments (b) water with T ~ 4°C sinks through bottom water due to sediment load and, after sediments settle out, water rises to surface.  69  6  Conclusion  Deep water ventilation in Quesnel Lake, which controls the oxygenation of deep water and recycling of nutrients to surface waters, is of importance to fish within the lake and impacts the salmon population in the Fraser River. The objective of this research was to develop an accurate equation of state for Quesnel Lake for use with CTD profiler data in order to be able to investigate deep water ventilation. Water samples and CTD profiles were collected in July 2003. Water samples were analysed for concentrations of constituents, pH, alkalinity and specific conductance. From these quantities, an equation relating in-situ salinity to CTD parameters was developed: C T D  0.836 • C£ ~~ (0.5369 + 0.0156 • T + 1.339xl0- • T - 7.552xl0- • T ) • f (P, T) 4  2  7  3  P  where the pressure function taken from Wuest et al. (1996) is given by: /p(p,r) = i +  7  -P  and 7 = 10~ • (1.856 - 0.05601 [ d b a r ^ ] • T + 0.0007[d6ar- C- ] • T ) 5  -1  -1  lo  The coefficient of haline contraction was determined to be PQuesnel  2  2  0.827 m g / L  =  - 1  or PQuesnel — 1 -096 • Pchen+Muiero for the composition of Quesnel Lake. The coefficient of haline contraction was used to modify the equation of state developed for waters in the limnological range by Chen and Millero (1986) as follows: PQuesnel =  0.9998395 + 6.7914 • 1 0 • T - 9.0894 • 10~ • T • • • -5  6  2  +1.0171 • 10~ • T - 1.2846 • 10" • T • • • 7  3  9  4  +1.592 • 10" • T - 5.0125 • 10~ • T • • • 11  5  14  6  +1.096 • (8.181 • 10~ - 3.85 • 10~ • T + 4.96 • 1 0 • T ) • S . 4  6  -8  2  where S is determined from CTD data using the above lake-specific equation. 70  Based on the values shown in Fig. 20, one can conclude that E O S Q L accurately estimates density for Quesnel Lake to within 0.0018 k g / m ; significantly better than 3  the alternative, E O S C M , which is accurate to 0.0158 k g / m . 3  The equation of state was developed in order to obtain accurate estimates of density to investigate vertical stability in a water column. Generally, when investigating vertical stability, potential density is used to remove the hydrostatic compressibility in order to compare water parcel densities within a profile. Also, standard density has been used for similar reasons. However, neither standard density nor potential density can be applied to the investigation of deep water ventilation of Quesnel Lake, where water temperature lies near TMD-  Quasi-density profiles show generally increas-  ing density through the water that provide a more likely representation of stability than either potential or standard density profiles. In Quesnel Lake, and other deep temperate lakes, quasi-density should be used to study local stability. Quasi-density profiles were generated from C T D data in the East Arm from 2001, 2002 and 2003 as well as the North Arm, West Arm, West Basin and Junction for 2003. East A r m profiles from all three years show very smooth quasi-density profiles. In 2003, features present in quasi-density profiles from the North Arm and West A r m are echoed in the profile from the Junction. Quasi-density profiles were used to compute profiles of the Briint Vaisala frequency, N , which identify where instabilities occur within the water column. To 2  find N , p uasi data should filtered by bin-averaging using a bin size that is appro2  q  priate for the size of density inversion that is of interest. N profiles calculated with 2  EOSQZ, and E O S C M were compared. There was not a significant difference in the stability predicted by each. N profiles calculated with potential density and quasi2  density were also compared. Potential density predicted a mostly unstable profile while quasi-desnity predicted a more stable and, therefore, more realistic profile. Briint Vaisala frequency profiles indicate that the water column in Quesnel Lake's East A r m is close to neutrally stable in the three years in this study, with weak instabilities occurring near 400m depth. However, Briint Vaisala frequency profiles  71  for the other areas in the lake in 2003 show larger regions of instability between 50m and 200m depth indicating that mixing may be taking place in those areas. Briint Vaisala frequency profiles showing neutral stability in Quesnel Lake's East Arm is evidence of a density current contributing to deep water ventilation. The limited transparency, as well as cooler surface temperatures through the summer months also support this theory.  72  7  Future Work  The natural progression of the research carried out on Quesnel Lake in future will include determining a Thorpe Scale for overturns identified in C T D profiles, analysing data available from a thermistor mooring deployed in the lake's deepest basin and determining the extent of deep water ventilation by the density current initiated at Niagara Falls in Quesnel Lake's East Arm.  7.1  Thorpe Scale  Where the water column is unstable, a Thorpe scale can be used to quantify overturn scales in individual profiles. Instabilities can be identified from density profiles by isolating regions where the density inversions occur (Hohmann et al., 1997). The Thorpe Length, named after Steven A . Thorpe who introduced the concept in 1977 (Thorpe, 1977), can be used to estimate the average eddie size in a region of density inversion. The Thorpe Length is the root mean square of the displacements required to reorder a density profile to be monotonically increasing (Thorpe, 1977). Dillon (1982) compared the Thorpe Scale to the Ozmidov scale, a scale for eddie size predicted using the local turbulent energy budget, and found that, in regions away from the surface, they were in good agreement. Galbraith and Kelley (1996) outline a procedure for identifying overturn size from C T D data using the Thorpe Scale.  They address issues of spatial resolution and  outline two tests that are used to isolate real overturns due to density inversions from systematic or instrumental error.  These tests are the run length test, which sets  out a minimum number of reordered points differentiating real overturns from noise, and the water mass test, which differentiates temperature-salinity spikes from real overturns. This method was illustrated with two profiles, one from the Saint Lawrence Estuary in Canada which is actively mixed and highly stratified and one from near Svalbard in Norway which is less actively mixed and is more weakly stratified. The method showed that with C T D data that is able to resolve density difference on the  73  Figure 35: Map showing location of Mooring 9 order of 1 0 k g / m , overturns thinner than 2m cannot be detected if N 0.003s -3  3  -1  (i.e.  weakly stratified). Stansfield et al. (2001) applied the Galbraith and Kelley (1996) method to C T D profiles from the Juan de Fuca Strait which were already reordered and profiles with random noise added to them to evaluate if this method would reject false overturns. Stansfield et al. (2001) concluded that Galbraith and Kelley's method is robust. This method should be applied to Quesnel Lake C T D data to assess if it can be used to identify overturns in Quesnel Lake's weakly stratified water.  7.2  M 9 Thermistor Chain  A thermistor chain, called M9, was deployed in the East Arm (see Fig. 35) on September 24th, 2003. Five Brancker TR-1000s and 4 Vemco thermistors were attached to a kevlar line starting from the bottom and set to take temperature readings every hour over the winter months. Table 6 outlines the structure of the mooring. Note distances are measured from the bottom. The year-long time series of the thermal structure in Quesnel Lake's East A r m produced by this data will has not been analysed to date.  74  Its analysis will allow  Table 6: Depths at which thermisors are attached to Mooring 9 in Quesnel Lake's East Arm. Distance from Bottom  Im 100m 200m 300m 400m 425m 450m 480m 500m  one to answer questions of whether or not the deep water renewal is an episodic or continual process. Also, presumably the deep water temperature doesn't change significantly throughout the year and data from this mooring will allow us to verify this hypothesis. If the ventilation is indeed episodic, future field work could include the deployment of two thermistor chains in the East Arm. This will allow one to determine from which direction deep water is being ventilated.  7.3 Further Investigation of Density Currents Transmissivity data, collected in the summer of 2003, could give insight as to weather the limited transparency of the water in Quesnel Lake's East A r m is due to density plume. Also, unloading of density current water could be investigated with settling experiments using Niagara Creek river water.  The Department of Civil Engineer-  ing at U B C has recently acquired an autonomous underwater vehicle (AUV) capable of simultaneously measuring in situ conductivity, tempeature, pressure, optical backscatter, chlorophyll A , and water velocity at depths of up to 500m. Further research could include use of this A U V , coupled with C T D profiles to determined the extent of the density current and potentially an Acoustic Doppler Current Profiler to determine the velocity of the density current.  75  References Bradshaw, A . and K . Schleicher (1965). The effect of pressure on the electrical conductance of sea water. Deep-Sea Research 12, 151-162. Campbell, I. J. (2001). Drawing  No. CP01247.  Sidney, B.C.: Coast Pilot, Ltd. ph.  no. 250-656-2109. Carmack, E . C . (1979b).  Combined influence of inflow and lake temperatures on  spring circulation in a riverine lake. Journal of Physical  Oceanography  9, 422-434.  Chen, C . and F . Millero (1976). The specific volume of seawater at high pressures. Deep-Sea Research 23, 595-612. Chen, C . and F . Millero (1977a). Effect of salt content on the temperature of maximum density and on static stability in lake Ontario. Limnology  and  Oceanogra-  phy 22(11), 158-159. Chen, C . and F . Millero (1977b). The use and misuse of pure water pvt properties for lake waters. Nature 266, 707-708. Chen, C . and F . Millero (1986). Precise thermodynamic properties for natural waters covering only the limnological range. Limnology  and Oceanography  31, 657-662.  Dauphinee, T . M . (1980). Introduction to the special issue on the practical salinity scale 1978. IEEE  Journal  of Oceanic Engineering  OE-5(l),  1-3.  Dillon, T . M . (1982). Wertical overturns: A comparison of thorpe and ozmidov length scales. Journal  of Geophysical Research <?7(C12), 9601-9613.  Eklund, H . (1963). Fresh water: Temperature of maximum density calculated from compressibility. Science 142, 1457-1458. Eklund, H . (1965).  Stability of lakes near the temperature of maximum density.  Science 149, 632-633. Ekman, V . W . (1934). Review of das bodenwasser und die gliederung del atlantischen tiefsee by georg wust. ICES  Journal  Of Marine  Fisher, T . S. R. (2002, April). Determination  Science 9, 102-104.  of Salinity  and Density in a  Meromictic  Mine Pit Lake. Ph. D. thesis, Department of Civil Engineering and Department of Earth and Ocean Science, University of British Columbia, Vancouver, B . C .  76  Fisheries and Oceans Canada (2001). Fraser River Sockeye Escapement Near Final  Estimates,  BC Interior  Area.  Estimates;  Fisheries and Oceans Canada.  Data  Report. Foster, T . D. (1995). Abyssal water mass formation of the eastern wilkes land coast of antarctica. Deep-Sea Research I 42(A), 501-522. Galbraith, P. S. and D. E . Kelley (1996). Identifying Overturns in C T D Profiles. Journal  of Atmospheric  and Oceanic Technology 13, 688-702.  Gill, A . E . (1982). Atmospheric-Ocean  Dynamics.  Academic Press.  Grafe, H . et al. (2002). Ultrasonic in situ measurements of density, adiabatic compressibility and stability frequency. Limnology and Oceanography ^7(4), 1255-1260. Gray, C . B. J . , R. C. Wiegand, and R. A . Kirkland (1979). Data report on the limnology of kootnay lake, british Columbia, conductivity temperature salinity relationships. Technical Report 6, National Water Research Institute, West Vancouver, B.C. Hill, K. D., T . M . Dauphinee, and D. J . Woods (1986). The extension of the practical salinity scale 1978 to low salinities. IEEE  Journal of Oceanic Engineering  OE-11,  109-111. Hohmann, R., R. Kipfer, F . Peeters, G . Piepke, , and D. M . Imboden (1997). Processes of deep-water renewal in lake baikal. Limnology and Oceanography ^2(5), 841-855. Imboden, D. M . and A . Wuest (1995). Mixing mechanisms in lakes. In A . Lerman, D. Imboden, and J. Gat (Eds.), Physics  and Chemistry  of Lakes (Second ed.).  Springer. Johnson, D. (1989). Conductivity temperature coefficients for the great lakes. Technical report, Environment Canada. Honours Agriculture Semester Five. Johnson, L . (1964). Temperature regime of deep lakes. Science 144, 1336-1337. Killworth, P. D., E . C. Carmack, P. F . Weiss, and R. Matear (1996).  Modeling  deep-water renewal in lake baikal. Limnology and Oceanography 41 (7), 1521-1538. Lewis, E . L . (1980). The practical salinity scale 1978 and its antecedents. Journal of Oceanic Engineering  OE-5(l),  3-8. 77  McDougall, T . J . (1987). Thermobaricity, cabelling and water-mass conversion. Journal of Geophysical Research 92, 5448-5464. McManus, J . , R. Collier, J. Dymond, C. G . Wheat, and G . L . Larson (1992). Physical properties of crate lake, Oregon: A method for the determination of a conductivityand temperature-dependent expression for salinity.  Limnology  and  Oceanogra-  phy 37(1), 41-53. McManus, J . , R. Collier, J. Dymond, C . G . Wheat, and G . L . Larson (1996). Spatial and temporal distribution of dissolved oxygen in crater lake, Oregon.  Limnology  and Oceanography ^i(4), 722-731. Millard, R. C , W . B. Owens, , and N-. P. Fofonoff (1990). On the calculation of the brunt- vaisala frequency. Deep-Sea Research 37(1), 167-181. Millero, F . J . (2000a). Effect of changes in the composition of seawater on the densitysalinity relationship. Deep Sea Research I J±7, 1583-1590. Millero, F . J. (2000b). The equation of state of lakes. Aquatic  Geochemistry  Millero, F . J . , C . T . Chen, A . Bradshaw, and K . Schleicher (1980).  6, 1-16.  A new high  pressure equation of state for seawater. Deep-Sea Research 27A, 255-264. Millero, F . J . , A . Laferriere, and P. Chetirkin (1977). The partial molal volumes of electrolytes in 0.725m sodium chloride solutions at 25°c.  Journal  of  Physical  Chemisty 81, 1737-1745. Morton, K . F . and I. V . Williams (1990). Sockeye salmon oncorhynchus nerka utilization of quesnel lake, british Columbia. Technical report, Canadian technical report of fisheries and aquatic sciences. Peeters, F . , G . Piepke, R. Hohmann, and D. M . Imboden (1996). Description of stability and neutrally buoyant transport in freshwater lakes. Limnology  and Oceanog-  raphy 41(8), 1711-1724. Perkin, R. G . and E . L . Lewis (1980). The practical salinity scale 1978: Fitting the data. IEEE  Journal  of Oceanic Engineering  OE-5(l),  9-16.  Press, W . H . , B. P. Flannery, S. A . Teukolsky, and W . T . Vetterling (1987). Recipes, The Art of Scientific  Computing.  78  Cambridge University Press.  Numerical  Robinson, R. A . and R. H . Stokes (1959). Electrolyte  Solutions (Second ed.). London:  Butterworths. Royal, L . A . (1966). Problems in rehabilitating the quesnel sockeye run and their possible solution. Technical report, Internation Pacific Salmon Fisheries Commission, New Westminster, B . C . , Canada. Sherstyankin, P. P. and L. N. Kuimova (2002). Thermobaric stability and instability of deep natural waters in lake baikal. Doklady Earth Sciences 385(2), 247-251. Shimaraev, M . N . , N . G . Granin, and A . A . Zhdanov (1993). Deep ventilation of lake baikal waters due to spring thermal bars. Limnology and Oceanography  38(5),  1068-1072. Sorensen, J. A . and G . E . Glass (1987). Ion and temperature dependence of electrical conductance for natural waters. Analytical  Chemistry  59, 1594-1597.  Stansfield, K . L . , C. J. R. Garrett, and R. K . Dewey (2001). The probability distribution of the Thorpe displacement within overturns in Juan de Fuca Strait. of Physical  Journal  Oceanography 12, 3421-3434.  Stockner, J . G . and A . C. Costella (1980). The paleolimnology of eight sockeye salmon (Oncorhynchus nerka) nursery lakes in British Columbia, Canada. Technical report, Canadian technical report of fisheries and aquatic sciences. Strom, K . M . (1945). The temperature of maximum density for fresh waters. Geofysiske Publikasjoner  16(8).  Thompson, W . F . (1945). Effect of the obstruction at hell's gate on the sockeye salmon of the fraser river.  Technical Report 1, International Pacific Salmon Fisheries  Commission. Thorpe, S. A . (1977). Turbulence and mixing in a Scottish loch. Philisophical actions of the Royal Society of London. Series A, Mathematical  and Physical  TransSci-  ence ^ ( 1 3 3 4 ) , 125-181. U N E S C O (1981, 1-5 September 1980). Tenth report of the joint panel on oceanographic tables and standards. In UNECA ume 36, Sidney, B C , Canada.  79  technical papers in marine science, Vol-  Vollmer, M . K . , R. F . Weiss, R. T . Williams, K . K . Falkner, X . Qiu, E . A . Ralph, and V . V . Romanovsky (2002b). Physical and chemical properties of the waters of saline lakes and their importance for deep-water renewal: Geochimica  et Cosmochimica  Lake issyk-kul, kyrgyzstan.  Acta 66(24), 4235-4246.  Weiss, R. F . , E . C . Carmack, and V . M . Koropalov (1991). Deep-water renewal and biological production in lake baikal. Nature 349, 665-669. Wuest, A . (1999). Equation of state for arrow reservoir water. In R. Pieters (Ed.), Arrow Reservoir Limnology and Trophic Status Year 2 (1998/99) Report, Volume 5. Province of British Columbia Ministry of Environment, Lands and Parks. Wuest, A . , G . Piepke, and J . Halfman (1996). oclimatology  of the East African  Limnology,  Climatology  and  Pale-  Lakes, Chapter Combined Effects of Dissolved  Solids and Temperature on the Density Stratification of Lake Malawi, pp. 183-202. Toronto: Gordon and Breach.  80  A A.l  Polynomial Coefficients Limiting Equivalent Conductivity  Coefficients for the third order polynomials for limiting equivalent conductivity for ions in large concentration in Quesnel Lake based on data published in Robinson and Stokes (1959). Where limited data is available (i.e. less that 4 temperatures), a second order polynomial is used and a "-" is used to mark the absense of the forth coefficient. A°° = A + B-T Ion Sodium Chloride Magnesium Calcium Nitrate Sulphate Potassium Carbonate Bicarbonate  A 26.4118 40.9701 28.9000 31.1911 40.0000 41.0000 40.6780 36.00 24.63  + C-T  2  B 0.7707 1.2891 0.6957 0.9299 1.2175 1.4251 1.1817 1.3804 0.70136  81  + D-T  (42)  3  C 0.0082 0.0054 0.0107 0.0087 3.7871 x 10-5 0.0054 0.0061 -0.001073 0.0053739  D -4.5001 x 10-5 -1.0670 x 10-5  -2.8035 x 10-5 6.3940 x 10-5  -3.1577 x 10-5  -0.00006485  A.2  Reduction Factor  Coefficients for the third order polynomials for the reduction coefficient for-ions in large concentration in Quesnel Lake.'  fi = A + B-x + C-x  2  where x =  + D-x  (43)  3  Vconcentration. Ion Sodium Chloride Magnesium Calcium Nitrate Sulphate Potassium Carbonate Bicarbonate  A 0.9999 1.0000 1.0008 1.0011 0.9991 1.0022 1.0003 1.0202 0.9999  B -0.8248 -0.6271 -2.516 -2.3785 -0.6579 -2.0353 -0.674 -2.5638 -0.8871  82  C 1.586 0.6855 5.8629 5.4562 -0.1461 3.2068 0.9295 2.4655 1.2214  D -0.8723 -0.2695 -6.3597 -5.6325 0.7767 -2.0426 -0.658 2.8852 -4.2751  B  The Carbonate System  T h e carbonate system, which plays an important role i n buffering acids and bases, contains most of the carbon found i n natural waters and controls the transport of carbon through air-surface interfaces.  Gaseous carbon dioxide i n the atmosphere  equilibrates w i t h aqueous carbon dioxide i n the water and i n doing so enters the carbonate system. T h e carbonate system is represented i n the diagram below as a series of equilibria between the various components of the carbonate system.  co ( ) 2  C0  2  s  — (aq) X  *  H C0 2  HCO;  3  1  C O "  2  Figure 36: D i a g r a m depicting the equilibrium of components of the Carbonate System.  A l t h o u g h a l l three components contribute to the density of the water i n which they are dissociated, two of three components of the Carbonate System are ionic. For this reason, it is important to find the concentrations of the components and identify their contribution to the conductivity of the water so that the water's density can be determined through conductivity measurements. T h e dissociation constant i n the reaction i n which [H COs] is ionised, ki, is defined 2  by:  k l  _ ~  [HCO?][H«] [H C0 ] 2  ( 4 4 j  3  Similarly, for the ionisation of [HCO3 }, the dissociation constant, k , is defined 2  k 2  ~  [HCOt ] 1  S3  ( 4 5 )  Water is also ionised to a small degree and its dissociation constant, k , is denned w  in a similar way:  Equations 44 and 45 can be rearrange to give formulas for finding respectively [H C0 ] 2  3  and [CO3 } in terms of [HCO3 *] as follows: 2  [H C03] = [HCO-, ] • ^  (47)  1  2  [CO3- ] = [HCO?\ 2  •- H _  ( 4 8  )  These two equations give the concentration of two of the three components of the Carbonate System in terms of pk\ and pk which are found empirically and pH which 2  can be measured from a water sample. The third, [HCO^ ],  will have to be found  1  indirectly through the measurement of alkalinity. The definition of total alkalinity, TA,  is the concentration of all bases that can accept TA = [HCO3 } 1  H: +  + 2 • [COf] + [OH' ] - [H ] 1  (49)  +1  Note that when converting alkalinity concentrations from mg C a C O s / L to equivalents/L, one must remember that both anion and cation have a valence of two and therefore this factor must be included in the conversion. By substituting Eq. 48 for [CO^ ] and solving for [HCO^ ], 2  [ H C  °*  1  ~  J  E q . 49 becomes:  [H+1] + 2 • k  ( 5 0 ) 2  Or, by substituting Eq. 46 with [H 0] = 1, E q . 50 becomes: 2  _ [  H  C  U  ! _ TA[H"] -k + [H^f12 ~ [H^rigkt]+2-k w  3  (51)  J  2  Or, ,  [HC0 \ 3  TA-10-P -10-^+ H  =  1 0  (10- ) pH  ,  2  - P * + 2•  1  Now, all three parts of Carbonate System are given in terms of pk , pk , pk x  2  w  .  ' and  pH. The temperature dependence of these variables at low salinities have been found  84  empirically and are summarised by Millero (2000). The following empirical formulae assume that the water is so dilute that it can be considered fresh water (i.e. 5 = 0): - 45.0575 • ln(T)))  pfci = - log(exp(290.9097 pk  = - log(exp(207.6548  pk  = -log(exp(148.9802  2  w  11843.79  3  3  6  4  g  5  _  l  n  (  r  )  )  )  (53) (  5  4  )  T 13847.26  • 23.6521 • ln(T)))  (55)  where T is in Kelvin. With the above equations coupled with measurements of pH and total alkalinity, the temperature dependence of the concentrations of the ionic components of the Carbonate System can be found.  85  C  Matlab Equation of State Toolbox  The Equation of State Toolbox for Matlab has been developed to assist in the delopment of equations of state with lake water chemistry data. The toolbox produces coefficients for the expression of a lake's temperature factor, f , T  and the constant of  proportionality, A, required to produce an equation relating in-situ salinity to C T D parameters. In order to relate salinity to density, the coefficient of haline contraction must be determined by the procedure described in Section 3.4 of this text. The .m files printed below require chemistry data files for input that takes the specific form of chemfile=[pH alk Na NH4 K C a Mg CI F l N 0 3 S04 C25 Si] where all concentrations are in units of m g / L . The flow chart in Fig. 37 illustrates the sequence of steps carried out by the toolbox. A n initial script would contain the names of the chemistry files and C T D profiles. Here, densitycalc.m is used as an example, densitycalc.m calls relationships.m which in turn calls C25.m. C25.m, which calculates conductivity of the solution over a range of temperatures, calls a series of operations including: carbonates.m (which calculated the concentration of the components of the carbonate system), equivalents.m (which finds charge difference and equivalent concentration of sample), L E C . m ( which finds limiting equivalent conductivity (LEC) for ions over the given temperature range), reductioncoeff.m (which finds reduction coefficient for ions over equivalent concentration range) and keff.m (which finds the effective conductivities for each ion). The conductivity over a range of temperatures outputted to relationship.m which then uses it to calculate the coefficients of fr-  salinity.m is called to calculate the total  salinity and the ratio between salinity and specific conductivity — required for determining A. The coefficients of f  T  and A are outputted to densitycalc.m which then  uses the C T D measurements to calculate the in-situ salinity and density profiles.  86  densitycalc.m Input chemistry and CTD profile data. relationships.m C25.m equivalents.m carbonates.m  T  LEC.m  T reductioncoeff.m  1 1  1  keff.m  T  Average f coefficients for all samples, Fit A values to line. 1 salinity.m  -T  Calculates salinity based on CTD parameters.  * Calculates density with equation of state modified to use lake-specific coefficient of haline contraction. Figure 37: Flowchart of Matlab Equation of State Toolbox.  87  Cl  density calc.m  % This s c r i p t i s used t o c a l c u l a t e an i n - s i t u density p r o f i l e °/from CTD p r o f i l e s - and chemistry data. 0  clear a l l °/name relevant chemistry data f i l e s and CTD f i l e s 0  chemf iles=-[' f ilename'}; ctdfiles={'filename'>; b i n s i z e = 8 ; % s i z e of b i n to smooth quasi-density p r o f i l e Zr=0;  .  "/.reference depth from with t o s t a r t c a l c u l a t i n g quasi density °/find c o e f f i c i e n t s f o r f i n d i n g s a l i n i t y from CTD data 0  relationships(chemfiles) load coeff.mat for 1=1:length(ctdfiles) °/load c t d data and get data s t a r t i n g at reference depth 0  XX=load(ctdfiles{l}); C=XX(:,1); T=XX(:,3); P=XX(:,2); K=find(P>Zr);  %throw out data before reference depth  P=[Zr; P(K)]; C=[C(K(1));C(K)]; T=[T(K(1));T(K)]; [PP,I] =sort(P); °/,sort data so that pressure i s i n c r e a s i n g CC=C(I);  88  TT=T(I); for  i=l:length(P)-l  "/.pick out indices with non-reoccuring pressures p(i)=PP(i+l)-PP(i); end PI=find(p~=0); new_P=PP(PI)  "/.pick out data with non-reoccuring pressures  new_C=CC(PI) new_T=TT(PI) 7 calculate s a l i n i t y 0  for  j=l:length(new_P) g=(10~(-5))*(1.856-0.05601*new_T(j)+0.0007*new_T(j)"(2))  7,to get f p from  fit+T  fp=l+g*new_P(j); S (j) =A*new_C (j) / ((ft4+ft3*new_T(j) +f t2*new_T(j) " (2) +'f 11. *new_T(j)-(3))*(fp)); end S=S'; SS=S/1000;  "/.convert to rhocm u n i t s  "/.calculate i n - s i t u desnity from E0S$_{QL}$ for  j=l:length(new_P) rhoQL(j)=rhocm_Quesnel(SS(j),new_T(j),new_P(j));  end rhoQL=rhoQL*1000; end  C.2  relationships.m  "/.RELATIONSHIPS Develop l a k e - s p e c i f i c r e l a t i o n s h i p s . '/. "/.  [ ] r e l a t i o n s h i p s (chemf i l e s )  89  y. "/,  Develops lake s p e c i f i c r e l a t i o n s h i p s f o r C25=C25(C,T,P)  "/,  and determines a value f o r the empirical constant  "/.  chemfiles with format of filename-[pH a l k Na NH4 K Ca. . .  "/.  A using  Mg CI F l N03 S04 C25 S i ] .  •/. "/. "/,  Function c a l l e d by  densitycalc.m.  '/. %  Input: name of v a r i a b l e containing chemistry  "/,  Creates  "/,  Saves c a l c u l a t e d c o e f f i c i e n t s as coeff.mat  f i l e names.  .mat f i l e which contains A and f t c o e f f i c i e n t s .  '/. '/,  Author: C h r i s t i n a James, 2004.  function [ ]-relationships(chemfiles) % Calculates C25 and makes polynomial  f i t f o r temperature f a c t o r  a=zeros(size(chemfiles)); "/.vector - measured s p e c i f i c conductance b=zeros(size(chemfiles));  '/.vector - c a l c u l a t e d s p e c i f i c conductance  c=zeros(length(chemfiles) ,4); '/.matrix f o r c o e f f i c i e n t s of f t for  i=l:length(chemfiles) [C25exp,C25calc,ftfit]=C25(chemfiles{i»; a(i)=a(i)+C25exp; b(i)=b(i)+C25calc; c(i,:)=ftfit;  end ftl=mean(c(:,1)); ft2=mean(c(:,2)); ft3=mean(c(:,3)); ft4=mean(c(:,4)); d=zeros(size(chemf i l e s ) ) ; "/.vector of s a l i n i t y  90  e=zeros(size(chemf i l e s ) ) ; for  "/vector of A  i=l:length(chemfiles) [S,A]=salinity(chemfiles{i}); d(i)=d(i)+S; e(i)=e(i)+A;  end % to make s a l i n i t y - s p e c i f i c % through  conductivity l i n e go  (0,0)  b=[0 b'] ; d=[0  d'];  °/ fit data to f i n d A 0  [fitA]=polyfit(b,d,1); Afit=polyval(fitA.b); % save c o e f f i c i e n t s % a .mat  c a l c u l a t e d i n t h i s function as  file  A=fitA(l); save coeff  C.3  A f t l ft2 f t 3  ft4  C25.m  °/ C25 Find r e l a t i o n s h i p f o r s p e c i f i c 0  conductance i n term of C , T , P .  % %  [C25exp,C25calc,ftfit]=C25(fname)  7. %  Finds r e l a t i o n s h i p between s p e c i f i c  °/  measured by a CTD ( i . e .  %  over the temperature range of 0 to 25 degC.  0  conductance and parameters  c o n d u c t i v i t y , temperature and pressure)  7. %  C a l l e d by r e l a t i o n s h i p s . m  •/. %  Inputs: filename  (pH alk Na NH4 K Ca Mg C l F l N03 S04 C25 S i )  '/,  Outputs: C25 c a l c , C25 exp, and f t f i t 91  7.  7o  Author:  C h r i s t i n a James, 2004.  f u n c t i o n [C25exp,C25calc,ftfit]=C25(fname) 7o load and declare v a r i a b l e s X=load(fname); t=0:25;  % temperature range of i n t e r e s t : 0 to'25 oC unless 7o otherwise s p e c i f i e d  X(find(isnan(X)))=0;  °/,sets nan=0  C25exp=X(13); 7 . sample's measured s p e c i f i c 7. f o r Na, NH4,  K, Ca, Mg, C l , F l , N03,  mmass=[22.99 18.0 39.1 40.08 24.31 charge=[l 1 1 2  2 2 1112  conductance  S04, HC03, C03  35.45 19.0 62.0 96.0 61.01  60.01]  12];  7 f i n d s the concentration of i o n i c components of carbonate system 0  [H2C03,HC03,C03]=carbonates(t,X(l),X(2)); 7o f i n d m o l a r i t y of each i o n (besides HC03 and C03) from measured 7,  concentrations  for  i=3:ll molarity(i-2)=X(i)/mmass(i-2);  end 7o f i n d s charge d i f f e r e n c e and equivalent concentration of sample [diff,equivConc]=equivalents(molarity,charge,HC03,C03); 7. f i n d s l i m i t i n g equivalent c o n d u c t i v i t y (LEC) f o r ions over 7. T range [l]=LEC(t); 7. f i n d s reduction c o e f f i c i e n t f o r ions over equivalent 7o concentration range [f]=reductioncoeff(equivConc,t);  92  "/, f i n d s the e f f e c t i v e c o n d u c t i v i t i e s f o r each i o n from LEC, "/, f i and concentration [k]=keff(molarity,charge, 1, f,HC03,C03,t); "/.finds c o n d u c t i v i t y of sample over range of temperatures for i=l:length(t) Ksample(i)=k(i,l)+k(i,2)+k(i,3)+k(i,4)+k(i,5)+k(i,6)+k(i,7)+... k(i,8)+k(i,9)+k(i,10)+k(i,ll); end "/.calculated s p e c i f i c conductance f o r comparison with measured C25calc=Ksample(26); "/.find r a t i o of CT/C25 ( i . e . f T f o r C25=f(T,C,P) r e l a t i o n s h i p ) for i=l:length(t) ft(i)=Ksample(i)/Ksample(26); end "/, f i t t h i s data f o r f T to a polynomial  and p l o t  [ftfit,s]=polyfit(t,ft,3);  C.4  carbonates.m  "/.CARBONATES Calculate concentration of components of carbonate system.  7. "/„  [H2C03, HC03, C03] -carbonates (t, pH, alk)  7. 7,  Using equations  "/,  Equation f o r c a l c u l a t i n g concentration of HC03 i s derived on  "/,  d e f i n i t i o n of a l k a l i n i t y , equation f o r c a l c u l a t i n g C03 derived  "/,  from d e f i n i t i o n of e q u i l i b r i u m constant.  f o r pKl, pK2, and pKw from M i l l e r o 1995.  y. "/,  A f u n c t i o n c a l l e d by C25.m.  7. "/,  Input: temperature, pH and a l k a l i n i t y  */,  Output: concentration of H2C03, HC03 and C03 i n mmol/L  93  7,  %  Author:  function  C h r i s t i n a James, 2004. [H2C03,HC03,C03]=carbonates(t,pH,alk)  alk=alk*2/100.078;  % sample's a l k a l i n i t y i n mmol/L a v a i l a b l e f o r % accepting protons (charge of 2)  "/.finding the constants of d i s s o c i a t i o n for i=l:length(t) pKl(i)=6320.81/(t(i)+273.15)-126.3405+19.568*log(t(i)+273.15); pK2(i)=5143.69/(t(i)+273.15)-90.1833+14.613*log(t(i)+273.15); Kw(i)=-logl0(exp(148.9802-13847.26/(t(i)+273.15)-23.6521*... P  log(t(i)+273.15))); end °/ finding the concentration 0  of i o n i c  constituents  for j=l:length(t) HC03(j)=(10-(-pH)*alk/1000-10-(-pKw(j))+(10"(-pH))-2)*... 1000/(10"(-pH)+2*10"(-pK2(j))); C03(j)=HC03(j)*10~(-pK2(j))/10~(-pH); H2C03(j)=HC03(j)*10~(-pH)/10~(-pKl(j)); end H2C03=H2C03'; HC03=HC03'; C03=C03';  C.5  equivalents, m  "/EQUIVALENTS Calculate the equivalent  concentration  of a sample.  "/. °/ 0  [ d i f f , equivConc] =equivalents (molarity, charge, HC03, C03)  7.  7o  A f u n c t i o n c a l l e d by C25.m.  7.  94  7.  Inputs: molarity, charge, concentration of HC03, C03.  7.  Outputs: charge d i f f e r e n c e on ions, and equivalent concentration  %  of sample.  7.  %  A u t h o r : C h r i s t i n a James, 2004.  f u n c t i o n [diff,equivConc]=equivalents(molarity,charge,HC03,C03) °/find sum of p o s i t i v e 0  charges  f o r i=l:5; P(i)=[molarity(i)*charge(i)]; end pos=sum(P); °/find sum of negative charges 0  f o r i=6:10; N(i)=[molarity(i)*charge(i)]; end neg=sum(N)+HC03(l)+2*C03(l); % f i n d s d i f f e r e n c e between p o s i t i v e and negative charges diff=pos-neg; '/converts from meq/L to eq/L and f i n d sqrt of equiv concentration aveCharge=(pos+neg)*0.001/2; equivConc=sqrt(aveCharge);  C.6  LEC.m  %LEC Find l i m i t i n g equivalent conductivity over a range of 7,  temperatures.  7.  7,  [l]=LEC(t)  7.  %  L i m i t i n g equivalent conductivity c a l u l a t e d over a temperature  '/,  range appropriate t o the lake i n question ( i . e . T=0 t o T=25  7o  deg C) .  7o  data (Robinson, 1959); 2nd order polynomial where l i m i t e d data  7o  available  Uses 3rd order polynomial f i t s based on published ( i . e . l e s s than 4 temperatures).  95  % %  A f u n c t i o n c a l l e d by C25.m.  '/. %  Input: temperature vector  7,  Output: an n x m matrix where n i s the number of ions i n the  %  sample and m i s the length of the temperature  vector.  %  %  Author:  function  C h r i s t i n a James, 2004.  [l]=LEC(t)  for i=l:length(t) lNa(i)=26.4118+0.7707*t(i)+0.0082*t(i)"2-0.000045001*t(i)"3; lNH4(i)=40.2+1.2976*t(i)-0.001437*(i)"2+0.000113*(i)"3; lK(i)=40.678+1.1817*t(i)+0.0061*t(i)~2+0.000031577*t(i)~3; lCa(i)=31.1911+0.9299*t(i)+0.0087*t(i)"2-0.000028035*t(i)~3; lMg(i)=28.9+0.6957*t(i)+0.0107*t(i)"2; lCl(i)=40.9701+1.2891*t(i)+0.0054*t(i)"2-0.00001067*t(i)"3; 1N03(i)=40+1.2175*t(i)-0.000037871*t(i)"2+0.00006394*t(i)"3; 1S04(i)=41+1.4251*t(i)+0.0054*t(i)"2; lHC03(i)=24.63+0.70136*t(i)+0.0053739*t(i)"2-0.00006485*t(i)"3; 1C03(i)=36+1.3804*t(i)-0.001073*t(i)"2+0*t(i)"3; end l=[lNa' 1NH4' IK' lCa''lMg' 1C1' 1N03' 1S04' 1HC03' 1C03'];  C.7  reductioncoeff.m  7 REDUCT10NC0EFF Calculates the reduction c o e f f i c i e n t as a 0  '/,  f u n c t i o n of equivalent concentration.  7.  [f ]=reductioncoef f (equivConc ,t)  I  %  96  %  Uses polynomial f i t s to published data (Robinson, 1959; Dean,  %  1979) to give reduction c o e f f i c i e n t  '/„  equivalent concentration of sample.  as a function of sqrt of  7. 7  0  A function c a l l e d by C25.m.  7. %  Input: equivlent concentration (produced by equivalents.m),  %  temp vector.  7o  Output: m x n matrix of reduction c o e f f i c i e n t s where m i s  %  number of ions and n i s the length of the temperature v e c t o r .  the  7. %  A u t h o r : C h r i s t i n a James, 2004.  function  [f]=reductioncoeff(equivConc,t)  for i=l:length(t) fNad)=0.9999-0.8248*equivConc+l.586*equivConc~2-0.8723*equivConc~3; fNH4(i)=l-0.674*equivConc+0.6658*equivConc~2-0.3118*equivConc~3; fK(i)=l.0003-0.674*equivConc+0.9295*equivConc~2-0.658*equivConc~3; fCa(i)=1.. 0011-2.3785*equivConc+5.4562*equivConc~2-5.6325*equivConc~3; fMg(i)=1.0008-2.516*equivConc+5.8629*equivConc~2-6.3597*equivConc~3; fCI(i)=l-0.6271*equivConc+0.6855*equivConc"2-0.2695*equivConc~3; fN03(i)=0.9991-0.6579*equivConc-0.1461*equivConc~2+0.7767*equivConc"3; fS04(i)=l.0022-2.0353*equivConc+3.2068*equivConc"2-2.0426*equivConc"3; fHC03(i)=0.9999-0.8871*equivConc+l.2214*equivConc~2-4.2751*equivCoric"3; fC03(i)=l.0202-2.5638*equivConc+2.4655*equivConc"2+2.8852*equivConc~3; end f=[fNa'  fNH4' fK' f C a ' fMg' f C l ' fN03' fS04' fHC03'  97  fC03'];  C.8  keff.m  "/.EQUIVALENTS Calculate the equivalent concentration of a sample. "/  [diff,equivConc]=equivalents(molarity,charge,HC03,C03)  "/ "/  A f u n c t i o n c a l l e d by C25.m.  '/. ' "/  Inputs: molarity, charge, concentration of HC03, C03.  '/,  Outputs: charge d i f f e r e n c e on ions, and equivalent concentration of sample.  '/. "/  A u t h o r : C h r i s t i n a James, 2004.  f u n c t i o n [diff,equivConc]=equivalents(molarity,charge,HC03,C03) °/find sum of p o s i t i v e 0  charges  f o r i=l:5; P(i) = [molarity(i)*charge(i)] ; end pos=sum(P); "/find sum of negative charges f o r i=6:10; N(i)=[molarity(i)*charge(i)] ; end neg=sum(N)+HC03(l)+2*C03(l); °/finds d i f f e r e n c e between p o s i t i v e and negative charges 0  diff=pos-neg; "/converts from meq/L t o eq/L and f i n d sqrt of equiv "/concentration aveCharge=(pos+neg)*0.001/2; equivConc=sqrt(aveCharge);  98  C.9  salinity, m  "/SALINITY Calculates s a l i n i t y of samples. '/. 7.  . ' [S,A]=salnity(fname)  '/. "/  Calculates s a l i n i t y by summing concentrations  "/  componentes ( i o n i c  "/  between s a l i n i t y and measured s p e c i f i c  and non-ionic).  "/  C a l l e d by relationships.m.  of a l l d i s s o l v e d  Also determines the r a t i o conductance.  7. "/  Input: filename  "/  Output: t o t a l  (ie.chemistry f i l e entered  i n E0S_CJ.m)  s a l i n i t y i n mg/L and A i n mg*cm/(L*microS).  7. %  Author: C h r i s t i n a James, 2004.  function  [S,A]=salinity(fname)  X=load(fname); t=25; "/calculate s a l i n i t y at a reference  temperature  [H2C03,HC03,C03]=carbonates(t,X(l),X(2)); "/calcuates components of carbonate system H2C03=H2C03*62.01;, HC03=HC03*61.01;, C03=C03*60.01; "/convert back t o mg/L x=[X(3) X(4) X(5) X(6) X(7) X(8) X(9) X(10) X ( l l )  ...  X(13) H2C03 HC03 C03] ; S=sum(x); "/sum of a l l components d i s s o l v e d i n sample A=S/X(12);  99  (A)  


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



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