Loading...
HomeMy WebLinkAboutEvaluation of ice problems Hydro in AK 1 of 3 1982HYD 050 EVALUATION OF ICE PROBLEMS ASSOLI.ATED Wllh H DROELECTRIC POKER GEN,ERATIOH IN ALASKA: FiRal Repori to the Sure of Aiaska D?partment of Commares aal 4onomiC Develcpiprt Ccl•lt. sm 08-73-7-913 0801-E-11h or Con+rac! AEC81005-3 FINAL REPORT TO THE STATE OF ALASKA DEPARTMENT OF COMMERCE AND ECONOMIC DEVELOPMENT Hydroelectric Icing 08-73-7-958/08-71-6-114 or Contract AEC81005-3 EVALUATION OF ICE PROBLEMS ASSOCIATED WITH HYDROELECTRIC POWER GENERATION IN ALASKA: FINAL REPORT TO THE STATE OF ALASKA, DEPARTMENT OF COMMERCE AND ECONOMIC DEVELOPMENT by J. P. Gosink and T. E. Osterkamp Contents Summary Appendices 1. A Survey of Ice Problems at Hydroelectric Facilities 2. Results of Questionnaire Study of 26 British Columbia Hydroelectric Stations 3. Ice Problems at Swedish Hydroelectric Plants 4. Length of Open Water Reach Below a Dam or Reservoir 5. List of Published Articles SUMMARY Information regarding ice, its adverse effects upon hydroelectric facilities, and reliable methods to minimize these effects is sparse, located in obscure references or proprietary. As a result, the engineering information and expertise necessary to deal with ice problems is not normally found in U.S. engineering firms. Since there has not been any hydroelectric development in Interior Alaska or along its colder coasts, there is no core of engineering experience to draw upon. The primary objective of this project has been to acquire, document and develop the necessary engineering information base to be used by hydroelectric power planners, designers and operators to eliminate, avoid or reduce ice problems associated with hydroelectric power production in Alaska's cold winter climate. We proposed to accomplish this objective by compilation of state-of-the-art engineering information, applied research where appropriate and publication of reports summarizing current world-wide engineering practice and research information. During the first two years of this project, reports were published which include: I. A survey of manufacturers, available equipment, applicability (head, discharge, KW) and experience with northern climates; 2. A bibliography listing sources of information on small hydropower with critical annotations regarding the usefulness of each; 3. A brief survey of ice problems and mitigating procedures in hydroelectric facilities in Canada, Switzerland, and Scandinavian countries. This final report completes the limited objectives for this project as set out in the revised work plan May 28, 1982. These objectives, including the complete survey of ice problems and mitigation procedures at hydroelectric sites in Sweden and British Columbia, and the development and documentation of a water temperature model for downstream thermal predictions, are addressed in appendices 1 through 5 of this report. These appendices represent both compilations of existing international engineering experience and methodology, and original applied research. The knowledge gained regarding ice problems should be made available to Alaskan hydropower engineers and planners. Correct site selection procedure, knowledge of potential problems and the means to alleviate those problems is of great benefit to Alaska. This information will allow rational management decisions to be made both in the planning and operational stage of hydroelectric development. Appendix 1 A Survey of Ice Problems at Hydroelectric Facilities: Report to the State of Alaska Department of Commerce and Economic Development by Greg Penn T. E. Osterkamp J. P. Gosink Geophysical Institute University of Alaska Fairbanks, Alaska 99701 There has been some interest in establishing small hydropower plants in northern regions, particularly Alaska. However, cold climates pose some special problems. These problems can be severe, especially for small hydro or run -of -river plants. The problems are often solvable and can be dealt with, but they need to be addressed before beginning con- struction of a hydropower plant in a very cold climate. A survey of difficulties experienced by existing power plants in northern regions can help identify what problems are likely to occur and how to deal with them. The survey would be useful for planning and for existing facilities. This is a survey of ice problems that hydropower plants have had or are having. A brief description of a problem and possible solutions is followed by a list of hydropower plants known to have experienced the problem. The intent of the paper is to show that difficulties do exist and to indicate what those difficulties are and their severity and incidence. Information for the survey was obtained from brief questionnaires sent to hydropower plants and utility companies, letters received from people dealing with hydropower, telephone conversations with power plant personnel, and personal visits to hydropower facilities. Intake Blocked With Ice A common problem is the build up of frazil ice on intake trash racks. Frazil ice is produced in the absence of ice cover in turbulent, supercooled water. The small ice crystals are carried downstream where they are drawn into intakes and cling to the trash racks. Build up can be rapid resulting in reduced or no flow to the turbines. Intakes can also be damaged or blocked by surface ice on a reservoir as the water level drops to or below the intake. The problem of ice on the trash rcks has been dealt with successfully in many cases, but it still causes difficulties under certain conditions. Some methods of preventing ice blocking the intake are heating trash racks, back flushing, creating an ice cover to minimize frazil production, lowering the intake water velocity to decrease drawdown and ensuring that intakes are in deep water. Gold Creek Juneau, Alaska 1.6 MW maximum Dewey Lakes Skagway, Alaska 30-375 kW Manitoba Hydro Newfoundland and Labrador Hydro Forces Motrices de Mauvois in Switzerland 3 plants 835 million kWh Soderfors DalalvenRiver Sweden The problem is not severe in this case since the plant is usually shut down due to low flow when it is cold. Here the water level in the reservoir must be watched to ensure it is deep enough to minimize frazil ice and to keep surface ice above the intakes. Frazil ice on the trash racks. Trash racks have been heated to minimize the problem. Spillimacheen Columbia River British Columbia, Canada 5MW Western Mica British Columbia, Canada Bennett Dam Peace River British Columbia, Canada Trash racks have been heated to minimize the problem. Air bubblers are used to prevent icing at the intake. Rapid ice formation in temporary diversion tunnels threatened flooding during construction. Flooding Caused by Ice Hanging dams, ice jams, anchor and frazil ice can restrict the flow in the normal river channel causing floods. The floods are sometimes severe and damage property including the power station. Careful planning is needed so that equipment and structures are above possible flood levels. Blasting with dynamite to break up ice jams has been frequently tried. Ice booms have been used to help establish a stable ice cover. Often, water that normally goes through the turbines must be either stored in the reservoir to float ice jams free or used to float ice masses over a spillway. This can reduce power output significantly. Dikes have been built to prevent property damage from floods. Manitoba Hydro Winnipeg, Saskatchewan, and Nelson Rivers Irve Tolles Real Data, Inc. Manchester, N.H. Town of Peace River British Columbia, Canada Near the Bennett Dam Town of Taylor British Columbia, Canada Near the Bennett Dam Ice damming "has become very critical on several occasions necessitating the mobilization of forces to prevent the topping of coffer dams". A plant under construction would have had water level controls damaged if they had been in place. Floods in 1973 and 1974. Dikes have prevented the problem since. A cold snap in 1979 caused increased electrical usage which necessitated greater discharge from the turbines. This affected the ice front upstream of Taylor producing a large ice jam that caused a flood with an 18' rise in 48 hours. Carefully controlling the discharge minimizes this problem. Icing of Structures from Spray Falling water at a dam or falls creates spray which wets nearby structures. When the spray freezes, it can damage those structures, often due to the weight of the ice on them. An icy coating can also be hazardous to people who must work in the area. Usually the ice is manually chipped away when it becomes a problem. Heating structures to melt the ice is also a possibility. Pidgeon River Plant Vanderbilt, Michigan 11-100 kW Forces Matrices de Mauvoisin Switzerland 3 plants 835 million kWh Manitoba Hydro The problem is minor. Pelican Creek A leaky woodstove penstock ices a Pelican, Alaska walkway making its use hazardous less than 500 kW and difficult. Shut Down Due to Low Flow in Winter The water supply to a power plant usually decreases in the winter in cold climates, often drastically. Decreased output or complete shut down result. Low flow in a pipe increases the danger of freezing in the pipe. Water storage in a reservoir or piping water in from other drainages can help make up for periods of low flow. Pelican Creek Pelican, Alaska less than 500 kW Gold Creek Juneau Alaska 1.6 MW maximum Dewey Lakes Low flow increases the danger of Skagway, Alaska pipes freezing. Pipes have split 30-375 kW and the plant has been close to shutting down on several occasions. Significant energy is used to heat pipes. Open Water Downstream In the winter, when warm water (4°C) is discharged downstream of a power plant, it can flow several hundred kilometers before it cools enough to freeze. This creates a long stretch, or reach, of open water. The upstream edge of the ice cover may be thin or unstable due to fluctuations in discharge. The frozen river may be a significant roadway for men and animals, and river crossings may be frequent. An open reach can disrupt migratory routes, river crossings, and winter travel on the river in general. Open water exposed to very cold air produces ice fog which can blanket a large area and create hazardous driving conditions. An open reach might encourage frazil ice production causing problems for other facilities downstream. It may be necessary to build bridges to allow river crossings over open reaches or to put up signs warning of thin ice. There are no known effective methods of controlling ice fog. However, the length of the open reach can be decreased with the use of ice booms or by changing the topography of the river bed to decrease water velocity. Controlled discharge can be used to stabilize the upstream edge of the ice cover. Manitoba Hydro Columbia River Ice fog increases due to open water British Columbia, Canada reaches on the Columbia River are expected to be about 5%. Equipment or Structures Damaged by Ice in Any Way The hydropower plants listed here are those that responded affirmatively to the statement "Equipment or structures damaged by ice in any way" on a questionnaire. The nature of the damage is unknown to us at this time. Annex Creek Alaska 2.8 MW Pidgeon River Plant The problem is minor. Vanderbilt, Michigan '11-100 kW Forces Motrices de Mauvoisin Switzerland 3 plants 835 million kWh Hnchor Ice Anchor ice forms on fixed objects such as the stream bed or manmade structures. It can block intakes or restrict flow in a channel. An insulating ice cover or heated structures can decrease anchor ice formation. Sometimes it is manually chopped out. Sheldon -Jackson Junior College Ice in a 2000 foot flume must be Sitka, Alaska manually chopped out. 50 kW Bruce P. Sloat Anchor ice has blocked intake Lancaster, N.H. structures. 15 kW Reservoir Ice Problems Changing water levels in a reservoir cause the ice cover to move up and down. This can damage structures, particularly dams. Also, as the water level drops, the center of the ice cover may sag giving the ice a slope that is dangerous for men and animals. Structures exposed to the ice must be made sturdy enough to withstand its abuse. Warning signs and fences may be needed to keep people and animals away from dangerous ice cover. Upper Salmon Creek Alaska 2.8 MW Crystal Lake Petersburg, Alaska Ice several feet thick rubbing against the upper face of the dam. Minor damage to the dam. The problem was solved with aluminum facing. Other Problems with Ice Annex Creek Icing of transmission line conductors. Alaska 2.8 MW Dewey Lakes Water seeps under the earth dam creating Skagway, Alaska a glacier on the creek bed below the 30-375 kW dam and threatening penstocks that come out of the dam. Snettisham Icing on transmission line destroyed Juneau, Alaska the line APPENDIX 2 RESULTS OF QUESTIONNAIRE STUDY OF 28 BRITISH COLUMBIA HYDROELECTRIC STATIONS Report to the State of Alaska, Department of Commerce and Economic Development compiled by T. E. Osterkamp Greg Penn J. P. Gosink Geophysical Institute University of Alaska Fairbanks, Alaska 99701 BRITISH COLUMBIA HYDRO AND POWER AUTHORITY BOX ,212, $55 WEST HASTINGS STREET, VANCOUVER, B-C. V68 4T6 CABLE ADDRESS "INTERPOW" TELEX 04-54456 December 14, 1981 File: 1206.10 Mr. T. Osterkamp, Professor of Physics, Geophysical Institute, C.T. Elvey Bldg., University of Alaska, Fairbanks, Alaska, 99701, U.S.A. Dear Sir, This is in reply to your letter dated nil to our Peace Canyon Project, which was referred to me. Enclosed are a reference map and completed ice questionaires for 28 B.C. Hydro hydroelectric generation stations. All but two of the plants responded and, at those two, ice problems are not usually experienced. The capacity of both unreported plants exceeds 1000KW. We would appreciate receiving copies of your survey results. We request, if possible, six copies for distribution to our region managers who completed the questionaires. Yours very truly, G.M. Salmon TPK/rt Manager, Encls. Development Department Please complete this questionnaire and return it in.the envolope. provided. Write the name and complete mailing address of your facility in the space below. Name of facility: C4 C- Ll l c.-» F� Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW_ 11-100 kW_ 101-1000 kW V� more than 1000 kW Have you had any difficulties with ice? Yes ✓ No If not, mail this questionnaire without continuing. _ Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice................i/ Flooding caused by ice jams, hanging dams, etc. . Icing of structures from spray.... .......... -- Equipment or structures damaged by iceinany way._ Other. Please specify. f._ I �• y r r� l7-n n` 1 n 11i J-71 i t C 1 _ Thank you for your cooperation. 702 W Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. . liame of facility: 8Vt3'f'.:&l5r� ..Put a mark next to the appropriate response. -The capacity of our facility is: 0-10 kW 11-100 kW�J 101-1000 kW more than 1000 kW Have you had any difficulties with ice? Yes No 1/ If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. . Icing of structures from spray.................. ... Equipment or structures damaged by ice in any way._ Other. Please specify. Thank you for your cooperation. s,— 0._1 2 'L I oo M �J Geophysical iratiww. C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 h10NE-907479.7282 TELEX: 35414GEOPH INSTSBK L Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. Yams of facility: ALOC'rrIl Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW 11-100 kW 101-1000 kW more than 1000 kW-/— $loco k-VJ Have you had any difficulties with ice? Yes No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter..._ Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. _ Icing of structures from spray .................... - - Equipment or structures damaged by ice in any way._ Other. Please specify. Thank you for your cooperation. G:-OP')':z-` I:-sunrtr, C.L Elvey Building, University of Alaska, Fairbanks, Alaska 99701 . PH0r:E_ PQ7-479-72&2 TELEX: 2514 GEOPH MST SBK e Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. . Mame of facility: G.M. Shrum Generating Station Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW 11-100 kW 101-1000 kW more than 1000 kW XX Have you had any difficulties with ice? Yes X. No_ If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................— Flooding caused by ice jams, hanging dams, etc. .. X Icing of structures from spray .................... Equipment or structures damaged by ice in any way. Other. Please specify. Thank you for your cooperation. 2.416CDO 6N! sus Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. Name. Of facility, Peace Canyon Generating Station 1 _ Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW .11-100 kW 101-1000 kW�A more than 1000 kW_✓ `M00Cp k 4 Have you had any difficulties with ice? Yes X. No If not, mail this questionnaire without continuing. — Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................ — Flooding caused by ice jams, hanging dams, etc._.:-X' Icing of structures from spray........ — - Equipment or structures damaged by ice in any. Other. Please specify. — Thank you for your cooperation. Please complete this questionnaire and return it in the envolope provided. - Write the name and complete mailing address of your facility in the space below. Name of facility: 10e1Jr 1'AJ e5• �• Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW 11-100 kW— 101-1000 kW more than 1000 kW X �o5�boo k,u} Have you had any difficulti2s with ice? Yes No X If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. _ Icing of structures.from spray .................... — — Equipment or structures damaged by ice in any way. — Other. Please specify. Thank you for your cooperation. G ps ysio; Institute, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 PHONE: 907A747262 TELEX: 3iS16 GEOPH INST SSK .1 Please complete this questionnaire and return it in the envolope provided. _ Write the name and complete mailing address of your facility in the space below. Xame of facility; vtE rf}"S C7• S. Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW 11-100 kW_ 101-1000 kW more than 1000 kW_,>Z 5215eD LW Have you had any difficulties with ice? Yes No. If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................ — Flooding caused by ice jams, hanging dams, etc... Icing of structures.from spray.:..................® Equipment or structures damaged by ice in any way. Other. Please specify. Thank you for your cooperation. G :, .i•a' Insttuti, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 . PHOUL, 907-074778:2 TELEX- IS414 GEOPH INKT SRK Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. dame of facility: �/��//�/7 �• Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW 11-100 kW_ 101-1000 kW more than 1000 kW_Z Have you had any difficulties with i.ce? Yes No X If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................®_ Flooding caused by ice jams, hanging dams, etc... Icingof structures.from spray .................... Equipment or structures damaged by ice in any way._ Other. Please specify. Thank you for your cooperation. 601000 kul Geopt:ysia1 Institu-te, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 . PHONF- 4:7<7471F TFI FY• 19 14 nFrIPN IllC fMe M Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. Name of facility: Put a mark next to the appropriate response. The capacity of our facility is: / 0-10 kW 11-100 kW101-1000 kW more than 1000 kW V Have you had any difficulties with ice? Yes—V No_ If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter..._ Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. . Icing of structures from spray .................... Equipment or structures damaged by ice in any way. te- Other. Please specify. - T V � 'u..e \9- d up 414^&4 I: � ..post c Lo_� OLJLl, Q C-d:(, w.Qw1 ia��a ._ Carp—" � n ba, G{A J- 4t, i c Thank you for yoUr cooperation. Geophysical Institute, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 P Mllu c-any wee ��m ru cv. w... ncnnu we-r c.... S Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. / Name of facility: ,L0- Toi< d'a„I„/Z4-r.' c S747T,ae/ Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW_ 11-100 kW 101-1000 kW • more than 1000 kW ✓ 22fOf,0 �c� Have you had any difficulties with ice? Yes ✓ No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter...® Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. Icing of structures from spray .................... Equipment or structures damaged by ice in any way. Other. Please specify. 1 Thank you for your cooperation. Geophysical Institute, C.T. Eivey Building, University of Alaska, Fairbanks, Alaska 99701 PHONE: 907-479"7282 TELEX: 35414 GEOPH INST SSK Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below, Name of facility; i3,Q7a6c 21 ✓cZ CZ14x447-,,µ%G ST,a%oq/ % Put,a mark next_to the appropriate response. The capacity of our facility is: B # 0-10 kW_ 11-100 kW_ 101-1000 kW more than 1000 kW 2. = �w Have you had any difficulties with ice? Yes / No 928,E If not, mail this questionnaire without continuing. _ Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................ — Flooding caused by ice jams, hanging dams, etc. — Icing of structures from spray .................... Equipment or structures damaged by ice in any way. Other. Please specify. Thank you for your cooperation. Geophysical Institute, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 PHONE: 907A79-7282 TELEX:35414 GEOPH INST SSK M Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. Name of facility:,c= -el ✓ic%L G'c4/G-,-A%/,,J 5%A71Osc/ 0'L Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW 11-100 kW 101-1000 kW more than 1000 kW Have you had any difficulties with ice? Yes No nu If not, mail this questionnaire without contiing. _ Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. . Icing of structures from spray.. ... ........ Equipment or structures damaged by ice in any way. Other. Please specify. -- of-ah, U. Y Thank you for your cooperation. 'B1+Z = 42,8/000lc0 Geophysical Institute, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 PHONE: 907479-7282 TELEX: 35414 GEOPH INST SSK Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. . Mame of facility: Put.a mark next to the appropriate response. The capacity of our facility is: 0-10 kW 11-100 kW_ 101-1000 kW more than 1000 kW '� 42,000 Have you had any difficulties with ice? Yes No_✓. If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open °:rater downstream of the facility in winter..._ Intake structures blocked with ice ................ - - Flooding caused by ice jams, hanging dams, etc... _ Icing of structures from spray.....................o Equipment or structures damaged by ice in any way. Other. Please specify. Thank you for your cooperation. Geo-tlysical lrst;tvte, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 PHOt:E907479-72S2 TELEX: 35414 GEOPH INSTSBK Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. lame of facility: Put a mark next to the appropriate response. The capacity of our facility is: _ 0-10 kW 11-100 kW 101-1000 kW_ more than 1000 kwo o �•v Have you had any difficulties with ice? Yes No ./ If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. _ Open water downstream of the facility in winter... Intake structures blocked with ice ...... ......... Flooding caused by ice jams, hanging dams, etc. Icing of structures from spray:,:. Equipment or structures damaged by ice in any Other. Please specify. Thank you for your cooperation. M Please complete this questionnaire and return it in -the envolope provided. Write the name and complete mailing address of your facility in the space below. . Elam® of facility: 962rr/o1.'e 67Ze w� Put a mark next to the appropriate response. The capacity of our facility is: - 0-10 kW 11-100 kW 101-1000 kW_ more than 1000 kW_, Have you had any difficulties with ice? Yes No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter..._ Intake structures blocked with ice... .............. Flooding caused by ice jams, hanging dams, etc. _ Icing of structures from spray.................—• Equipment or structures damaged by ice in any way,® _ Other. Please specify. Thank you for your cooperation. G�ophysid inst njtt, C.T, Elvey Building, University of Alaska, Fairbanks, Alaska 99701v A— P"0NE907479-72S2 TELEX: 35414GEOPH INSTSSK Please complete this quesnare anhe tionid returnAt in:t provided:- '= • _ _ - _ - >J!•-r. s . •-_�� �T.d�\:-°n••' �i'JC _•�`••'•"p�i .. � •. -•: a-=^_ � •TY'"TT+t' _ :.eTf+! a: r _ Write the name and complete mailing address of your facility in the space below. = _- Name.of ficlIIty.- . • Put -a mark next. to the appropriate_ response.-`_ The capacity'of our facility isi.r` ^'- -= _ -: =- ° z:.;r:: ' .'•.� ; - 0-10 kW_ 101-10001V ' more than IOOO kW Have -you had any difficulties* with ice? Yes Ho_'•- If not, mail this questionnaire without continuing: - Put a mark next to problems associated with ice'and cold climate that have occurred at your facility.- _ Open water downstream of the facility in winter... Intake structures blocked with ice.>..._:�..:.:.<.— Flooding caused by ice jams, hanging dams, etc:'.— _- Icing of structures from .- Equipment or structures damaged by_ ice in any way._: - Other. Please specify. Thank you for your cooperation. ' `"= • - ".��' •a;',•i • r ::. =:..: �S -. 1.1,E ._�"� �� : :��a�r_ r _T-::_y �� :i . _ S - _ •'` _ 1 • �t f:�. �. ••%E'er •%L. 1 ,•.• • _r.i �+•�iT• •• /T ..' t..:.. -mow. ':�: «.. -... -_ ="; �••7 •y:._ C.T. Elrry Huldrg, Unirersitof Alaska, Fairbank%Alaska 99707 _ =• `.:'�-,. _-� r�• `��.s^; �-.�•,;,J� ft*NLi07.47P7= TELEJ:2U1/GEOPHINSTS7K5CS'�-�-:iiM�*-:i , V. �� +�`• +i7r�.r•.�:.33.�rs�i •s !!'R_•il-'a:L'LX'f�isf�.ri y+�.z'-.yam:.-:-,2� �:.2.�.:.1':ix. .:t-_z-�-..,, .-_ � _ .� Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. liame 'of facility: Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW_ 11-I00 kW_ 301-1000 kW_ more than 1000 kW 007)$"00 40 Have you had any difficulties with ice? Yes No­� If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter..._ Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. . Icing of structures from spray .................... Equipment or structures damaged by ice in any way._ Other. Please specify'. Thank you for your cooperation. Geofiysial 4nstituty, CT. Elvey Building, University of Alaska, Fairbanks, Alaska'99701 PHONE °07A79-7262 TELEX: 35414 GEOPH INSTSSK _ _, M Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. Name of facility: C L /:V ej .a .S7s r: :- +. Put a mark next to the appropriate response. The capacity of our facility is: - 0-10 kW 11-100 kW_ 101-1000 kW_ more than 1000 kW Have you had any difficulties'with ice? Yes No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter,.. Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc... v Icing of structures from spray._..:.. Equipment or structures damaged by ice in any way.® Other. Please specify. Thank you for your cooperation. Cxophysicsl itutit:.v, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 PHONE 907<79-77R7 TFI cv• vu.e renew �.�� _ .- E G Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. Hame of facility- IAQ5 W 64Q Put a mark next to the appropriate response. The capacity of our facility is: - 0-10 kW_' 11-100 kW 101-1000 kW more than 1000 kW Have'you had any difficulties with ice? Yes ,/ No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. - Open water downstream of the facility in winter... Intake structures blocked with ice ................. Flooding caused by ice jams, hanging dams, etc. . Icing of structures from spray.. t....... .=/ Equipment or structures damaged by ice in�any�way. Other. Please specify. Thank you for your cooperation. Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. lame of facility: k,.)F�ATS l�� CI�NEi�ATINC� S�L1TtZ,/J l Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW .11-100 kW_!d- 101-1000 kW more than 1000 kW_� Have you had any difficulties with ice? Yes No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter..._ Intake structures blocked with ice ....... Flooding caused by ice jams, hanging dams, etc. . --- Icing of structures .from spray.... ....::.....:.... Equipment or structures damaged by ice in any way. . Other. Please specify. — Thank you for your cooperation. Geophysical In -juju. C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 _ _ . PHONE-907479-7282 TELEX. 35414 GEOPH INST SSK -�1���74Y y a ii MR, k AM Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. lame of facility: Wo�`�rer ��c �v.a�. Cle-,, — Cg000 }Zw� Put a mark next to the appropriate response. The capacity of our facility is: - 0-10 kW 11-100 kW 101-1000 kW more than 1000 kW Have you had any difficulties with ice? Yes / No If not, mail this questionnaire without continuing. _ Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ............... Flooding caused by ice jams, hanging dams, etc. Icing of structures from spray .....a.... .. Equipment or structures damaged by ice in. way.® Other. Please specify. Thank you for your cooperation. Please complete this questionnaire provided, d return it in the envolope Write the name and complete mailing address of your facility in the space below. iw-m- of facility: J 0 06 f1 N G Eft' !«V 6 is Put a mark next to the appropriate response. ` The capacity of our facility is _ / 0-10 kW 11-100 kW 101-1000 kW more than 1000 kW ✓ Have you had any difficulties with ice? Us No If not, mail this questionnaire without contTnuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. 150d000 k—w Open water downstream of the facility in winter... Intake structures blocked with ice ................ Flooding caused by ice Jams, hanging dams, etc. _ Icing of structures from spray.; ........... ­ - * pray.:.............:. Equipment or structures damaged by. ice in any way.=---�-14o-t�t,t�vy CotVc U A Vr=-- Other. Please specify.' Thank you for your cooperation. ft Please complete this questionna ' and return it in the envolope provided. Write the name and complete calling address of your facility in the space below. i pane or teciiitys- = A`, 94 4r? " Put a mark next to the appropriate -response. = Thc-capacfty of our fact ity�is• - 0-20 kW 11 100 kW ` 101 1000 kW�more than 1000akW :.' Have' -you had any difficulties with ice? Yes No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ...............:_ Flooding caused by ice jams, hanging dams, etc. . _ Icing of structures from spray ...................._ Equipment or structures damaged by ice in any way, Other. Please specify._ Thank you for your cooperation. r 2512-00 kv\ 6yL."S s,4o,.-7 Ga0%n S'0 //-0/%Yy r0/s9f<'✓gas a/S 0 6i i--L".leC-d 9K"/1� /o rQ i` t < O ry R�/ •` ✓ ¢ w 0I'177s L _ g t af3 Q l- �1fr c,__�iyt_Co�4c YaI✓4 G*op yicil instibrce C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 _ M40NE907-479-7282 TELEX:35414GE0PHIN57SBK Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. ` Same of facility: 1,cz�. Put a mark next to the appropriate response. -The capacity of our faci1' y is: 0-10 kW 11-100 kW 101-1000 kW_ more than 1000 kW V/. 27"'o k'J Have you had any difficulties with ice? Yes No 1/00 If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. - Open water downstream of the facility in winter... Intake structures blocked with ice ................—` Flooding caused by ice jams, hanging dams, etc. Icing of structures from spray ......... ....... .... — Equipment or structures damaged by ice in any way. Other. Please specify. Thank you for your cooperation. ' __" —. • may. - - -_ _ _ ... • n = � .J . 1 �. 3' -`' _ _ -` _ lam' _TGwphys dl 1nstiLrta, C.T. Elvey Building, University of Alaska, Fairbanks, Alaska 99701 _ • - _ • -' _ • _t rHONEc D07479-7282 TELEX: 3S414 GEOPH INST SBK Please complete this questionnaire and return it in the envolope provided. _ Write the name and complete mailing address of your facility in the space below. p Mame of facility: _ Put'a mark next to the appropriate response. The capacity of our faciljy is: 0-10 kW .11-100 kW V4 _ 101-1000 kW more than 1000 kW_�. (01 SCO �-) Have you had any difficulties with ice? Yes No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have Occurred at your facility. - Open water downstream of the facility in winter..._ Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. Icing of structures from spray .................... Equipment or structures damaged by ice in any way. Other. Please specify. Thank you for your cooperation. Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. ` iame of facility: Put a mark next to the appropriate response. The capacity of our fa c M is: 0-10 kW_ .11-100 kG1101-1000 kW_ more than 1000 kW_ Have you had any difficulties with ice? Yes No r� If not, mail this questionnaire without continuing. ' Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................® Flooding caused by ice jams, hanging dams, etc. �. Icing of structures from spray.. ................. . Equipment or structures damaged by ice in any way.® Other. Please specify. Thank you for your cooperation. 54)000 V-W Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. �f ia:ne of facility: - Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW .11-100 kW 101-1000 kWmore than 1000 kW ✓ 120� 000 LW Have you had any difficulties with ice? Yes No If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter..._ Intake structures blocked with ice ................ Flooding caused by ice jams, hanging_dams, etc... Icing of structures from spray .................. _ Equipment or structures damaged by ice in any way.® Other. Please specify. Thank you for your cooperation. Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below.; lame of facility: Put a mark next to the appropriate response. The capac�y of our facility is: 0-10 kW_ 11-100 kW_ 101-10.00 kW_ more than 1000 kW_I/ Have you had any difficulties with ice? Yes ✓ No If not, mail this questionnaire without continuing. _ Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Cj) (goo kVJ Open water downstream of the facility in winter... Intake structures blocked with ice ................ F-rnZ l sc�- Flooding caused by ice jams, hanging dams, etc..._ Icing of structures.from spray ...:...........:.... _ y Equipment or structures damaged by ice in any way. — Other. Please specify. Please complete this questionnaire and return it in the envolope provided. Write the name and complete mailing address of your facility in the space below. Name of facility:/ 5 i, w.. n i } (1 Put a mark next to the appropriate response. The capacity of our facility is: 0-10 kW 11-100 kW 101-1000 kW more than 1000 kW ✓ . l32 Have you had any difficulties with ice? Yes No-C If not, mail this questionnaire without continuing. Put a mark next to problems associated with ice and cold climate that have occurred at your facility. Open water downstream of the facility in winter... Intake structures blocked with ice ................ Flooding caused by ice jams, hanging dams, etc. Icing ofstructuresfrom spray.... Equipment or structures damaged by ice in any way. -- Other. Please specify., Thank you for your -cooperation. APPENDIX 3 ICE PROBLEMS AT SWEDISH HYDROELECTRIC POWER PLANTS by Lennart Billfalk an informal translation by the Geophysical Institute Report to the State of Alaska, Department of Commerce and Economic Development ICE PROBLEMS AT SWEDISH HYDROELECTRIC POWER PLANTS by Lennart Billfalk Hydraulics Laboratory Swedish State Power Board Alvkarleby, Sweden ABSTRACT In the following report, a summary of ice problems in hydropower plants in Swedish rivers is presented. The rivers from Lagan in the south, up to Lule river in the north, are included. The information has been collected by the Hydraulics Laboratory, by sending out questionnaires on two different occasions (1977 and 1973), both to the managers of the power plants, and to private companies. Mainly, four types of ice problems can be distinguished: (1) Ice pressure against dams and gates, and freezing -up of gates and blades, in the construction, because of ice build-up, (2) ice formation on intake gates, blades, or turbines, (3) ice floe and slush against the intakes, and (4) ice dams. As far as gates are concerned, the Power Board, in the beginning of the 40's, decided to have effective heating methods to be drawn up to guarantee the maneuverability. There has been further development in improving the safety of running the movable parts of the structures. In order to avoid ice pressure against the dams and gates, a small opening in the ice cover has been made in front of these constructions. The use of current generators, foot -lights, and, in certain sites, air bubblers has shown that nowadays it should not be any technical impossibility to eliminate ice pressure risks that affect sensitive structures. As far as ice forming is concerned it seems to be that the heating of intake gates usually results in lessening the problems, but ice forming can start in spite of the installed gate warmers. The effect that has been brought about by heating the iron of the gates is self- ` evident, but it seems to be difficult to eliminate entirely the risk of ice formation at certain stations. The only way to surely avoid icing problems is to bring about a fast and lasting freeze-up upstream of the station. Problems with ice floes and slush affecting the intake and the origin of ice dams are partly connected to how the power stations are regulated. Short term regulating, when there is a great difference between the daytime and nighttime use of water, increases the risk for this type of ice problem (even icing -becomes greater since the freeze-up has been made more difficult). In the following report (part 1), the description of ice problems, and the measures to prevent troubles are summed up and commented upon. The report has been based on the information submitted by administrators and power companies. The part 2 of the report is made up of a combination of the letters that were obtained in answer to the questionnaire that was the first one to be sent out in 1977. Finally it is noted that information about ice problems (information was obtained from the first questionnaire in 1977) was summed up as a contribution to IAHR's ice symposium in Lulea in 1978 [5]. 1 INTRODUCTION In the year 1937, the general director of the State Power Board gave C. E. Soderbaum an assignment to work systematically at questions that were connected with ice difficulties of hydraulic power plants. In the presence of a threat of a war, there were questions concerning the maneuverability of the gates, and after a study trip during the month of March, Soderbaum and Witalis compiled a report: "Ice problems at hydro power plants(1)". To sum up, it is stated in (1): "There is a need for greater safety against ice difficulties that the power stations have been facing as the technology has been developing. Thus, the gates, where icing can be feared to occur, are fitted with effective warming systems. When it it necessary, there are many different ways to protect buildings against ice pressure. In some waterfalls, the results have been satisfactory, but at this point a further development seems to be desirable. As far as mobile constructions are concerned, methods with fully proven effectiveness have been used, and as such, make it possible to have the gates constantly operating even during winter." Finally it is noted in the report (1) that "Proposal has been submitted that current (generators) and other underwater pumps with piping should be installed at several different power stations. By means of underwater pumps, the water is made to circulate and thereby ice build up is hindered. In the year 1959, the director of the Hydraulic Laboratory, Stig Angelin, urged the administrators at the hydro sites to compile various experiences that they have gained as far as many kinds of ice problems are concerned. Part of these experiences compiled, because of the request, were presented in an ice meeting held in connection with the IAHR's 8th Congress in Montreal, 1959 [(2),(3),(4)]. In January 1977, in order to have a general survey of current ice problems in our Swedish hydropower plants, the hydraulic laboratory sent forward a questionnaire partly to the administrators of powerplants and partly via VAST to a number of large power companies. As a result of this questionnaire a lot of good information about ice problems was obtained from certain power stations in sections of some rivers; while information from other regions was scarce or failed to come. In order to complete this inquiry paper a detailed questionnaire was distributed in March 1978 to power companies that administer power stations in Lagar, Kolbacks, Dal River, Ume River, and Lule River. The purpose of this other questionnaire was to get -as detailed a description as possible of the ice problems in Swedish power stations in a number of rivers from Lagar to Lule River in the north. Importance is placed this time also on getting information about the design of the power plants and other matters that are important in understanding the causes of the ice problems. Definition of Ice Terminology Terminology within the realm of ice is not entirely unambiguous. It also appears to vary within the country. An attempt is made below to briefly describe some phenomena. Ice build-up (slush, swell) Frazil ice particles are formed on open stretches of rivers and lakes as the water supercools. There the frazil ice particles grow in size and form small flo? and gradually they "fasten" into numerous round ice floes. In rivers, ice covers build up partially from the shores (surface ice) and partially from collections of ice particles which are 2 in their different stages of development. Frazil ice particles, as they evolve, fasten on all kinds of things that come their way. That way anchor ice builds up and blocks the intake of a power plant by freezing the gates. Ice dams Ice floes that drift with the current along an open stretch of a river can jam against the surface of the ice edge. They can either be stopped and build up an ice cover of loosely packed ice floes or be sucked under the ice barrier. The maximum flow velocity for the ice floes to build up an ice cover is about 0.6 m/sec. For collection of slush (and thin ice floes) the maximum speed is lower and depends on its thickness and strength. For loosely packed ice floes the maximum flow velocity can be 0.2-0.3 m/sec. If upstream of an ice edge, the stream velocity'is so great that ice floes get sucked under the ice cover and get deposited under the ice, the water level rises and the speed of the river flow is reduced. As a result ice can then accumulate against the edge of the ice sheet causing the ice cover to grow upstream. Ice damming can also start as a result of ice which has run aground. It grows then on stones and shores in such a fashion that damming occurs. Inquiry 1 (paper 1) The questionnaire with the requested informtion is given in Appendix 1. This information is concerned in part with ice barriers. In the appendices 2.1-2.7 a survey of our existing Swedish powerplants over 10 MW is shown. A number of stations less than 10 MW are also reported. 3 Answers to the questionnaire concerning ice problems for each power plant and river are reported in the appendices. The answers that have come in concerning ice problems and ice barriers are accounted for in a separate report (Ice problems in Swedish hydropowerplants, part 2). The most essential problems concerning existing ice problems are reported, besides what appears in appendices 2.1-2.7 in the section under "comments" in the following report. Inquiry 2 (paper 2) In March 1978 a detailed questionnaire was distributed to the owners and administrators of hydropowerplants of Lagan, Kolbacksan, Dalalven, Indalsalven, Umealven, and Lulealven. After a certain amount of persuasion, information has been obtained from the aforementioned power plants which are on selected rivers. The respondants did not answer all the questions, which, moreover, were often not too skillfully worded. The unanswered questions are marked with --- in the charts. Any left out answer that applies to questions about ice problems is, nevertheless, interpreted as a problem that does not occur. Certain questions could be misunderstood. Moreover, it appears that people have been less inclined to answer all the questions about power plants if there were no ice problems occurring. A summary of the reported ice problems is compiled in a chart form in the tables 4.1-4.15. Powerplants from six rivers are compiled in the report. In these tables information that is judged to•have direct interest for purposes of comparison between different power stations, is also given. The remaining information that was asked for in the questionnaire is not included. 2 (It can be obtained by request to the Hydropowerplant Laboratory). It ought.to be pointed out here that time has not been granted to work with the left out information. Below are comments about the meanings and thoughts behind the headings on the tables. Power plants Mean rate of flow and corresponding velocity at the upstream intake. All power plants upstream from the mouth of the river are considered. Information on the mean volume of water transported is obtained from the mean value for a winter day and winter night. The available information comes from hours of different discharge in the nearby power stations without the regulating possibilities. Such information cannot, therefore, be totally accurate. Outward bound stream speed has been calculated from the water mass transported and the width and depth of the power plant. Especially where there are no canals, it can be assumed that the speed of the river can be lower than calculated (max. depth stated). EVALUATION OF ICE PROBLEMS ASSOCIATED WITH HYDROELECTRIC POWER GENERATION IN ALASKA: Final Report to the State of Alaska Department of Commerce and Economic Development Contract 08-73-7-958/08-71-6-114 or Contract AEC81005-3 The size of the areas not covered by ice immediately upstream of the stations. The heating of gates Icing If one supposes that freezing occurs if the speed of the stream falls below 0.6 m/sec., it appears that there is an inconsistency between the calculated speed of the stream and data about ice -free areas. The inconsistency can depend upon the incorrectly calculated speed of the stream or that the data about the size of ice -free surfaces is valid during the freeze- up period. It is not clear from the questions in the questionnaire if the warmed part of the total gate area is concentrated in one (or several) intake while gates in front of the remaining intakes lack heating. Information about icing occurrence on turbine.blades (ledskenor) and gates and whether ice gatewarmers are switched on automatically or manually and at what temperature this occurs. There is information also about production drop as a result of icing. 6 Ice dams Presence of ice dams (description of location and cause goes under "other") as well as information about production drop as a result of reduced height of the fall. Ice floe Problems with ice floe against the intake. Information about causes, point of time, necessary work, and drop in production. Other Added information above all about icing problems and ice dams. Comments Besides the six rivers included in the paper 2, a good picture of the state of affairs in many power stations in Gbta River, Motala River, Mar River, and Skellefte River is obtained from the paper 2. Below are the comments on the information received about these ten rivers. Lagan According to Sydkraft's description, ice problems arise mainly in power plants where upstream is a river bed, inlet canal, or a tunnel. Ice problems.almost never start in plants that are directly connected with a reservoir. According to appendix 4.1, however, two out of three power plants which have intakes in a reservoir, indicated in paper 2, icing problems appear. The reason for this probably is that the reservoirs in question are relatively shallow. 7 Icing is found yearly in six (seven?) out of the eight power plants in spite.of the installed gate warmers (exception is Aby power station). It is also stated that problems occur even when the gate warmers were turned on before ice formation started. In three (four) of the power plants with icing problems icing occurs -not only on gates but even on turbine blades (ledskenorma). The gate warmers are turned on already when the water temperature is between 0.3 and 0.5 C. This is unusually early compared with the practices of other stations. If problems with operating the gates occur at +0.3 to 0.5 C, it will not probably be because of icing on the iron of the gates. In such high water temperatures it can be assumed that the ice problems occur when great quantities of slush block the gate openings. In certain canals strong anchor icing occurs at the bottom (1-1.5 m). When this anchor ice gets loose, it can contain stones (of 10-13 kg) which together with the ice can cause problems to the gates and turbines. In control dams the ice pressure causes damage on vertical level gates. It should be investigated if laying out generators could help this problem. Sydkraft has taken different steps to reduce the ice problems. In power stations, where the water upstream is constant, stream flow (velocity) is reduced during the freezing period to increase the rapidity of the ice formation upstream. Even laying out ice booms is used for this purpose. With variable upstream water levels, the canals must ordinarily have open water surfaces. The great difficulty is to keep the canals ice q free at the start. Ice is driven out according to a certain system of starting, stopping, and driving out the ice through the ice outlet many times and using a certain amount of acquired experience of opening ledskenor (blades?, guides?). Clean up work in the waterways has also brought about as a result a retarding effect on icing in some installations. In four plants they have installed automatic gate cleaners for general gate cleaning work (such as leaves, litter, ice slush), but the cleaners have even proved to have a good effect 'against ice formation. It is thought that the gates are given a certain vibration and that the cleaners break the ice film on the upstream side of the gate iron and, therefore, the handling of ice has become considerably easier. Sydkraft maintains that without taking proper steps, it would be scarcely possible to have power production at all times of ice difficulties. This is especially so since nowadays and there is a lack of personnel to solve difficult ice situations. Gota River To a certain degree, ice problems in G8ta River are special because of ship traffic which has its demands for passability. On the power plant's side, the people hope to have a fast and as complete a freeze up as possible to avoid icing problems that otherwise appear moderately often, especially at Lilla Edets power plant. The shipping office on their part hopes as a principle to keep the river ice -free for shipping. During recent years the shipping office and the power plant office have found certain solutions that both sides can accept. The objective is to reduce ice production in the river and this can be done by letting a great part of the river form an ice cover. Freeze-up can be facilitated 9 by placing ice booms in strategic places. Further the power plant's side aims at keeping a uniform and low discharge, when conditions for freeze-up occur. This, however, is difficult since power stations in Gbta River usually must do short term regulating. In the river even ice dams are sometimes built up, especially on the stretch of Lilla Edethavet. Ice dams appear, partially because of strong drift ice formation, partially as a result of anchor ice, and during certain years they have given cause to troubling floods. Motala Stream Icing has occurred only a few times in Motala and Malfors stations. Icing has never occurred in Nykvarn. In"Malfon and Nykvarn problems with ice floes appear occasionally. Ice floes build up in the intake canals at nighttime if the stations are not operating. Ice floes then break loose with the starting of operation in the morning and travel down to the intake gates. Loss of head appears as a result. The most difficult freeze-up appears during Saturday and Sunday when the stations do not operate from Saturday afternoon until Monday morning. Klar River Of Uddelholm's nine power stations in Klar River, there are troublesome ice problems in Munkfors and Forshaga. In Edsforsen, Skoga and Deje there are certain troubles with icing before freeze-up. Icing or ice problems seldom or never occur in the remaining power stations. In Nunkfors and Deje stations, there are trash racks that work with icing and it seems that they make it possible to keep the operation going on. At Deje, during part of the year, timber is laid down to protect against icing. 10 Downstream of H51jes power station there are about two ten kilometer stretches of undeveloped (rapids). In spite of that, -the water temperature at the station can be about +1°C in February; ice slush and drift ice are built up along the rapids stretch. This ice gets stored up where the water flows slowly and gives rise to damming. Since short term regulating has been started.at Hbljes, the risk of the problems has increased and things have been closely watched. According to the information from the management, the ice forming is worst in very cold weather and when discharge is low. Kolbacksgn Kolbacksgn is included in (Inquiry 2) Paper 2 in order to present detailed information about a small waterway with small power stations. Icing problems are found in five out of the eight stations that are characterized as river power plants or that have intake canals to bring in the water. (Gatewarmers are lacking in all power stations in Kolbackson.) Reports on icing appear from this waterfall even before the freeze-up, but in Trangfors icing can occur during the whole winter day and night. Trangfors power station has a 400 m long canal that brings in water. It is open the whole winter. On four out of the five power stations which have icing problems there is ice growth on both gates and ledskeneapparaten (blade apparatus?). Smaller problems with ice floe against intake gates are reported from two stations. The problems occur with greater discharge changes respectively in milder weather. Dal River In all thirteen power stations between the confluence of the Vaster ii and Osterdal Rivers (Lindbyn included) and the ocean, icing has occurred 11 in the intake gates. Ice growth at both gates and ledskenor (blades?) appears in five stations in a row from Avesta-Lillfors to Lanforsen, at the lower part of the river. Occasionally in Domnarvet power station icing has occurred on the turbine blades in one of the two Kaplan turbines, this caused an imbalance in the unit and there was some damage as a result. Ice formation in turbine blades (in Kaplan turbines) has also been discovered at the Lanforsen power station. In Vaster River, icing occurs in Eldeforsen each autumn, and only seldom in Hummelforsen. Icing in Hummelforsen happens even though there are gatewarmers. In bsterdal River, icing occurs more seldom the further upstream the power stations in the river are. Ice dams occur regularly downstream of Trangslet and about 6 km upstream of Vasa power station. Flooding caused by these ice dams causes damage to the surrounding settlements. It is to be assumed that in Tra�ngslet the problems are reduced after the installation of air bubblers in the reservoir by which warm bottom water is lifted towards the intake. During the winters 1977 and 1978, ice damming occurred in one of the forks of the river, upstream of Untran power station. Ice damming during these two winters was attributed to the increased water flow of the fork of the river. It happened as a result of diversion of water in connection of building a power station in Soderfors. The ice dams caused the water level to rise above the dam, damaging it. In doing so approximately 75 m3/sec water flowed on the site of Untran power station and lfo . intdrained river basin. Upstream of �lvkarleby power station, bottom anchor ice can occur in the shallow part of the river which results in the rise of the downstream 12 water level in Lanforsens power station. Problems with the ice floes against the intake gates occur in several power stations mainly downstream of where Vaster and Osterdal Rivers flow together. In order to prevent ice pressure against dams and gates, current generators are used with good results in many stations. In Avesta-Storfors they seem to get good results by using gate cleaning machines, among other things, to remove ice that cannot pass through the ice outlets. The power companies agree that the most effective measures to prevent icing are to keep low and even discharge of water during the freeze-up period so that an ice cover is formed. Indals River Ice build up on intake gates occurs in Svarthalsforsen power station (3-6 days a year). In the Stuguns and M6rsils power stations, ice build up has occurred once since the plant started operating since 1976, and in Hammarsforsen power station once during the last 30 years. In the rest of the power stations there are no problems with ice build up. Downstream of Bergeforsen the river has ice every year on the shores and bottom, decreasing the area of the river. This causes a 0.7-0.8 m drop in the pressure head at the station. During February -March 1972, there were serious problems with ice slush (not ice build up) at Svarthdlsforsen. Portal cranes with ice scoops were working in three shifts a day and night to remove the slush. Svarthglsforsen company points out that at this time a lot of water was lost as a result of repairs at the Krangede and Gammeldnge stations and they speculate that ice problems in Svarthdlsforsen appeared because of this spill. No problems, however, appeared at the Hammarforsen power 13 station that lies between Krangede and Gamnelange and Svathalforsen. An explanation could be that the water temperature at the outlet from Hammarforsen was lower than normal as a result of cooling in the tailrace in Krangede and Gammelange. Ice slush build up was lighter than normal on the relatively silty river stretch between Hammarforsen and Svathalsforsen. The dam in Stadsforsen power station is not considered to withstand from any ice pressure, and therefore, a canal towards the dam is kept open. It is reported that chain saws, ice brakers and scoops of portal cranes and also snowmobiles are used for this work. If no special circumstances prevail in Stadsforsen, an.opening along the dam can be maintained with current generators or floodlights which is done in other locations. Breaking up the ice cover together with changes in stage and ice floe movement towards the intake occur in Midskogs, Mtlrsil, and Svarthalsforsens power stations. Ice floe movements towards the gates occur in some other stations in connection with freeze-up or break up. However, these ice floes do not usually constitute any more serious problems. As an exception, however, considerable difficulties can arise. In M6rsils power station, thaw and strong wind in January 1973 caused a 10-15 cm thick ice on Liten Lake to break up and obstruct all of the 3 km long intake canal. The ice masses were 5-6 m deep in the dam and caused the power production to drop for almost 24 hours. Ume River Icing on intake gates occurs at a few years intervals at Storrnorrfors and Bdlforsen power stations. In Hdllforsen and Betsele, icing has even appeared on ledskenorna (blades?). In B81forsen, Betsele and Hallforsen power stations, icing caused a total stop of power production in 1971, 1973, and 1975. 14 At Tuggen power station ice dams appear yearly 5-10 km downstream of the station. The downstream water level at the power station is about 2 m more than the surface in ice -free conditions. During the day, because of water discharge, the speed of the stream in the outflow canal is so high that the canal does not freeze. In cold weather, great masses of ice slush build up along the canal and along the ice -free river stretch downstream, and come to the damming area. There the speed of the river becomes lower and ice slush gets accumulated under the ice. Karteringar has shown that large sections of the river are obstructed by slush in this manner. As a result of the increased resistance, the water level rises upstream of such an ice dam and the stream velocity decreases and conditions for freeze-up occur. The ice cover that gradually covers a greater part of the canal is made up of slush and ice floes, hence, gets a moderately uneven surface. Ice floe accumulation along the edges indicates that the stream velocity has been very close to the limit of where the ice floes get sucked under the ice cover. Apart from extensive excavation work in Tuggen outlet canal, the only possibility to reduce the risk of ice damming is to restrict short term regulating during the time when good conditions for freeze- up prevail. Skellefte River Information has been only received from Skellfte power plant. Descriptions of ice conditions in Skellfte power station are given by an administrator as follows. 1. Icing in Finnfors, Granfors, Krangfors, and Selfors power stations 1973-11-10 and 1973-11-11. 15 Because of low water temperature, strong winds, and low air temperature, strong cooling of the river occured without ice formation. We measured water temperatures to an accuracy of one -thousandth of a degree C. Because of freezing there, and the lack of ice cover, the pressure head dropped rapidly at the gates. So, gradually, production had to be reduced very strongly at the power stations. Some turbines had to be stopped because the intake gates froze totally together. In Selfors power station, ledskenorna (blades?) in the turbines froze together so that it became impossible to use the gear shifts. After about 40 hours, the strong winds decreased and freeze-up on the dams of power stations became possible (an ice cover formed). The water temperature started to climb' up a few tenths of a degree and so the freeze-up began. After about 44 hours, the power production was resumed to full extent. Production drop during this period was about 630 MWH. It should be pointed out that the above mentioned problem is very uncommon. No information like that has ever come up before. Therefore, our actions for above mentioned type of problem are limited to trying to keep constant power. The most useful approach is to keep constant water flow during the freeze-up period. 2. Ice dams downstream of Granfors power station during the winters of 1975/76 and 1976/77. During these winters we have had ice dams above all in downstream Granfors power station. Because of that we got flooding in G1 turbine pits at Granfors on January 7, 1976. The loss of pressure head on this occasion was about 2.7 m (the height of the water head). The immediate 16 procedure was to reduce the water volume going through the station, and to install a warning signal for high downstream water level. (The station is far from our management center in Skelleftea.) Attempts to break away ice dams were performed with certain success, but the problem still remained partially during the winter 75/76 for about a month. The drop in water height (pressure head) varied between 1.7-2.3 meters, and the discharge was about 180-190 m3/s. We believe the best manner to avoid the aforementioned problems, even here, is to try to keep the discharge constant (and possibly low) during the freeze-up period, that is to say, to keep short term regulating as low as possible. 3. Outflow stretch downstream of Kvistforsen power station. About 1 km downstream from the central part of Skelleftea city, ice dams often appear. In the winters of 74/75�and 75/76 the dams were particularly big. Damming is caused by ice slush and piled up ice blocks. The cross-section of the river is considerably smaller where the dams start than in the upstream part of the river stretch. In January 1975 in the city center the maximum damming was measured at 0.9 m, and in January 1976 about 1.2 m. That the ice dams became so big in these years can, among other things, be caused by the unusually high discharge by the power plants during these winters. During the winter 1975/76, damming was concentrated to quite a short stretch, and that is why it was judged to be possible to break away the ice. On January 8, 1976, an explosive bursting of ice was carried out, and it resulted in an immediate decrease of the damming by 0.8 m and gradually decreased further. 17 Casualties of damming occur foremost through embankment overflows and water seeping into the community drainage system. These require increased pumping cost and the risks of overflowing wells. Besides, crowding together (narrowed channels?) has caused the water velocity to become so great increased bottom erosion- can occur. Lule River Apart from certain problems with ice floes against intake gates at Letsi and Akkats power station, and a moderate risk of icing at Boden and Laxede, serious ice problems have appeared only at Vittjdw power station. During the first winter after the present administration took over at Vittjaw power station they were forced to spill all of the water through spillway which is equipped with a so called "ski jump" in order to reduce the velocity of water downstream. This spill together with a couple of stretches of stream open downstream of the station produced a great amount of slush. This slush got accumulated in the rapids and caused damming that reduced the water pressure head (fall height) down to 2 m. After major dredging work downstream, it appears that the risk of ice damming is considerably smaller than earlier. Icing at the intake gates has been a yearly recurring problem ever since the present administration took over. As a test, gate warmers were installed on the intake in one or two sets. No reduction of icing or ice build up problems on the warmed intake could be established, however. Therefore, the gatewarmers were disconnected. Neither has the extensive cleaning work performed upstream in order to hasten freeze- up, have substantially reduced the icing problem. M. INVENTORY OF ICE PROBLEMS The nature of problems - Ice formation on gates, ledskenor,(turbine blades) and so on - Ice dams, anchor ice (reduced water head, flooding) - Ice pressure against dams and locks - Other The causes of problems and the measures which have been taken - Descriptions of how often and, if possible, why the special problems occurred - Which measures have been taken (or should/can be taken), examples - Reduced water flow (the volume of water) during the freeze- up period - Laying out of ice boors - Changing of water ways - Heaters on gates and other constructions - Laying out current generators, releasing the warmed water or similar measures - Other - What effect have the measures brought about (the degree of difficulty of the problems and frequency before and after the steps had been taken) Inquiry about ice booms. The purpose of laying out ice booms. - Accumulate floating ice and thereby initiate the forming of firm ice cover 19 Prevent drift -ice from reaching the intake constructions, where the risk of the ice being sucked under or blocking occurs Direct ice to ice outlet Other Placing of booms (give reason for the choice of the placing of booms please, enclose a sketch or direction) - The surface velocity in the selected section is so low that drift ice is gathering toward the boom - Other Necessary maintenance - Recording of booms in summertime - Exchange of the non-functioning parts Length of life - How long have the booms been in use - If a boom is exchanged; give its life span and the reason to the exchange Expenses - The cost of manufacturing - The cost of placing (laying out) them, including anchoring 20 LO • y � i 1 C � O aJ a1 C O t. O i > t N •r r ro a)YL L _� A L 4J Y O 7 0 0 N O LO ro R N U N U N C O cm n N C C E L ro aJ r o E ro d . LO O •r •r aJ ' L v C r v� 3 n N O (aO E 4- >1 N U •r aJ a) .Q o O N nE N aj � a1 CJ L L Y L coi yt c w =LYY 1 {� •r •r N .G L) 4J E N O N N R 1 0.a1 L •N O 'a N L Y Vl a) .. .O• S- E a De • ro i •O eaO)v ro m "•• E O R 4J L N L 4- L- C y O. a) S- Y C 4- N •r W O O 4- r N C1 U ro ro N a1 U U Y 4• o a +? E a9 ' r a7 0 a) . Y r 1N O N ro N E N rC N toN R i i Y O O C.- m. w N 7 N> O) I a1 r >r r () R cu N 0 Y Oi 4- 4- F- U 3 a) W"a L u O + •r O L L O R i Y a1 O tm a) W U R v'Y F••• _ro cc U .-Urrv. •.U- r-a X a L` a1 ro x x x x xx x v z c L O O U + Y ro i Y Y R L• N Y 9 U � 4Y••. Y N • i L N a� C QI Y = - Vro9 i N N O [n R N N NN l!7 N SH Y i a1 • 3 > � 3 erp cn C 3 S ko N Y 3 �� � N 3 3 f • f E E E Y r=Q�Y 3 �r � N 00 CD CDO CDCDi 3 33 0 O Yt0 Y f N C I I O CD CDI 1 I 1 O rl (n c%j I 1 1• W N r r r r L C r ro L .n Lr r r r r r r R r b r I r 0 N R r 3 YYY 4. Y Y Y Y Y Y Y Y Y .— Oo O L L r R 4- iT> Y N N N N N N N N N N N N Y r i YY Ln r tp W ^ LO r r lD l0 C• J F- v ^� ZN N CD L C F ro10 3 C r L >L. Y y N C � N E ro c C r c cn L c! C R 7 L a1 an O R y C R ro R R Y N ro N Y a) R i Y ro ro b a7 E Y •.- �6 9 ro N Y C J s S W Z G 9 to cm j- I E W r L N U m O 4- N > d m 3 L U y n + m > > r a) Lv N E L L C •1 L +•+ a) a) a) N L m a) U c� N U U v o a ,p •.- aasd .••. Q) lY Cr CL' 4- L O U • U t 4� U •d L d O U N m aJ� •r 7 3� E i.� E C C • c d J m i� N ev m N m > N m N m > E > c a) I La LU ar m 0 L an d E= E= armor L U N•.- CM E o N+L E O Ea E d- au) 1 n c� am),—� d O O 0 0 0 >> 4- L U O C 4- 1 o O O 4- 4- +� O CM +.I N 3 LU L. d4- O a O N c LL 1] O d4_ d4- L L LU . (v d d U O a L L Q7 L Y N a) L Y • E N O �••' L a) a) a) a) U (V N > O1 O N d O L U. y u dN U Y u 4• -,i - a) C4-•r 7 N 7 '. y'a r - r m-, m m L C m a) L I a1 O L L a)d 61 Q) 33 E3 a) ta) m a' N C N mLmOENL= ^ L C E 'o0ao >A >i>rLca a) E aCOa) L c E m aO I 7QU ZUZ3z-3)3Lm V �4- L m I CL W X v a •^'v a L ++ U X X >C X X X X X >C X 2 c'•, ) c � r-. L 0 0 -W o L? m U L 4J •C m U d-> U - 4- L m m E L N • O N a CY � G Y E U c m 0 1 \ O 00 N d ]C N N L L C L L n N a1 N m N = _ _ -- _ - _ _ m C T 4v, O L�R O = F- O O 3 f 3 O 7- 3 3 3 3 ^ 3 O S �. ON 23^ CIJ rm � E^ f Nr r O 3 N N r O E +� N N a) m � � o� C O N V N 7 r r L L N LT C CD 4 a 1 ^ M 4- cm a) =0 C C L Q1 O a 00 4-4- L O m O O '- .• ^ rEE • 4- •n m E Y .0 N cm Y a) 4- L N N m d 0 L m .y 00 (v - >,L L O c'n N L L O m r E E U .mi N S F-- W h Y L.L_ Ln g p U Y G+ U W Lv C O = =10 cu a m > > C N m m Y 'c L ?p v m U L. a za a) >- m L O .0 4� G Y > L O SO Ln G Y } Y y � L d CJ \ Y R al � L 1� C O Q • i Y n•'- E L of y E C c O 'a Y Io ter- r-- Y C Y a) L — V- CL N t0 T C) Ya E L O Y O U N N N N O • + o C N > > CJ Io C) p N N N y N N N y N 41 ra N N y C1 C) Y Y a) w Y Y CJ a) Q) Y Y Y C) a (m C) a) C) C C a N ,O R i0 (a R M fa C1 lu e0 m OI tm C OI 01 Cn C OIL 0) 01 to C C O E c c c c c c c c c Y c c • c N • O O O o O O o 00 N 000 O al p y E 'a ° a Im of of c m a a� a� OI v o c c Y c •ccc c c la c� ca c as U U oo .r .' .� O O L .U, M r-e ar U U r U U U U U U^ U O U a o-e r-i (--i ... �..• .-. �..� rr I -a Y rr e-Uw M >< x X x x X X X X X X Y OI U C-)O _ O U? e L y C N U L) H O O L y R 4- •C �a L) C) N Q _ • CL CL - - - - - - - M - - 'I"I C O O E O CL CL rt7 i 4- c U y Y o L a y �6 Y to at R L L Oo OL Y Baj L ne C _ - - _ - - - L T y L > cm Ilz N -- 41= a e= L>} E C) V7 G S tr> > fli cn C1 CZ C.B Y 3 3 �3 3icc un C0 o 3 3 C D N 3 3 y L M ® M: 0 to C 3 M i fco E c E 3 3 4- h et N S M 3 3 Lid Ln r� C d ^ f L 3 co �1 3 a1 y 'II i tD Y sr CO , C- Y 3 C O O N E C (3) tb Lo O) O N > C >4- > OIM^ N O C C1 Y O y� E Y C y�� N a) N •- L E > >I y L L fp •� p ro r— S- _ t.[)'-- M L L L 00 +� C CJ Y N L .a C tp O� N� C y N L > '1• Y V O i O L Y M O Y Y to L �O �O •p a1 y >> •r f N L C L r- Mr C o1'a N to E C a7 C 4- N N L O OILY M O 4- Y �i F- Cl L b, Q? CO N C7 L. > 7 aJ y Y C > O c W m O C L 0 -j Y C] D J N M: G Z O J Q C J N cz m J N O C C) > 4- O C rEo C O a) L CJ y L Q1 Q1 �— U tEp N A 'O cu C C al C N C)M 7 Y O > > i d QJ O Y O C Z N .N R Z •^ N O O N 47 L 0 4- 4..•N N 0)^•r•r 4- IO v- E O •C r c ai E r n .... N •.- U 41 ro v a w L + L r O C +1 Of 41 L C O C y o O L• n•�- • • N 4• C N Na •.r N c �N 10 o ro•-ro E EL E • c L i= L L v N ro y N N •^ ro N O y N r; r U Eva nv �•'• Cd a- T L O U O O +j U O r En d N 0 C L n M. N O n O 0 7 C O 0! a N N S- O •^ O.r 4� L i r O L n•r i. 4� a L d O r 4J C O O C _ _ - _ (0• ^ C Z a)nN _ a)Ln }T 4J t (D f N •-Ui f X d L tA O •^ V L°^ R XX X X X XX XX X d 7 U X X 3 a-•^ M i O • 0= Q U O O O U - U O U UU O= O ro . qy c0 44-- 4J C.) C c 4- U _ _ O L.)U O c R L 4- R c 3 4-� O N roUC- Y L N 0 YE0 ro4- V- O L L ro- U _ C L. 4J '0 44- - i L Y L S- CO Y y L.r i O N C N V- b L 4- C Y N O1 N vE f ro R N >C C of y Y N • W r- O 0 ro C O L w =0O L C d > ro Z L L N ro N 4J L •n N _ r 4- N ro O 4- Y O)r L w O N= L C L C s X L ro r L R C L y 4 J 4- C R L is Y r cy GJ 3 O o= d Y ®O7 O O .�C : = a) c ro > L O N R CO 4Y w» > Z N z> C Y H q cn Y Z O Y m v] O N N m 3 i > O CD O 3 E CD 3 (aCiE '- Ig CD co3� _ � 3 � ^^�� O N Gr r N �� f 0 L i�rCli M_ i Ecn 4J 3N v Lor3rM C N Z n� M M N N N 4 E r 4J l A l D E co N L n O ro N N L tp O N L L O w r 4- ro 0 l0 E G) E d• r L. C. 0 ro r d Q7 L N O O w 0) N i q N'aatip E-0 20 a• 0 7 N4-r N C O 4- a ^ v 4- r= O 4- L r zo ro L. O7 cy GJ N> W y C =0 4- M=V a C ro .0 a Y r •n ro 4� •ry 4- O O L 4� ro C N ro 4• o a tm C C C> 0 ro 01 0, N ^ w EN C N C •w 4- -]e N ro 4� C ir 0- C 7 4- N R• C 4j E C nr n r w 4� W C L to y 0 a d) •^J R R ro 4� L L. L L 3 L ro L L r 044 -0d ro w 'M L N J O p] S (z •n J J J C+ L =0 L+ F- d' O =0 ro F- 7 d v •n S J O -C �-•. cn, 7 b •r ro 4J > gV) S Y 0 Ln = v- =0 Z Z y C N L C dJ N C r0 ` E •n °J L o d v N 4 L N L L N 4� S-C O O N .. •r m V) C C N O �. .Y 0 4O.j D 0 0 L 'o r0 YO3 ) O a). - a) o- �4' AES- aj �� i•r O AU 3 ar •r N ••- O IV N E In 3 0) N N C r6 r- � 4- a O) U E E •r 4J i'' O .n (W cu c o o'er s c 4- r- T L C Ln • I 00 N •r O L L •r D. O O O CJ C O, N m E O) 0) Ol to Y O]'c O- r0 > w , 41 U 4� O c0 •r O) 4�r v (1) Ad rro Z Z Fu O r-r U U C V U U O i--r •r •r o--r f-r ,-+ J X ti � N � Y 'r r0 ccv X X X X X X X X X e o-•r n c a O U 44 L 4- O) O O O i eo O U O Oo U U y 4• C] U U U U 0 4- 2 4� 4- (U C S- 4� C 10 N 4- V- 4- N N L L O Y 4- a) L O re' e6 i L- C C 0= ]z O N L r0 Y Y P L L C X le N Y U > S- a: L O O Y • L C C N Y N S- 4) U U- N d N N O W d 4/ r° N N m > 0) C O) w N > N L d C U y o_ > \ C L L N c N O)'c L N tm C r y N U N a) c) Z• O O Y m r0 9 'y C C w d l C r6 d= S- L: , c y n_ L C N E= r 07 0) 0 r = N O zo O O) c b r• a c c c c N 4-. c 4- c } L r0 > C r0 t0 N T r- B L L r m= ^= 03 N Y m N }}}}} N N N V) N •r r J m 0 2 L c7 L L M Y Y m L O O r0 > R O L CD Z Y m v) m cm C7 3 � 3 3 2: a 3N z,:3: nc LnO Ln 3 E i OCl r 'd' LO O C 3 E r 3On C L ^ ^ N N i 3 3N r-CO f ilNO �. Ct v O co3 C C E C > i �� �� 3 CAD N a)N 4- 0 0) y N O N L 3 i L^g 0 CN 3 M F N C C R G) L N N N N lD O f") t,f) V-. r Lo O r CD C) Ol N l•n O/ C O r L L N r L N. 4- r- t0 C ll'! r• 0) N Y L N 6) 00 L L O C O) L L 40. O r 4- a) O 4- N C rd � L N N cn C O) 4� `� +� r L r- '.-i O O _ N (EE0 4� N N > L 0) > •r� i� 'O C r 4! 4- 0) r- y rd 4- N +•` r0 O O N M A N 4- 4- +� � C ^ C N N 6) 9 L o Q) Ql N V�i ^ = R O m E E rd r0 L r- S- C to > E N L ¢0 r- d m N C) r r. 'O W C S- r0 r0 > 4� 9 34 •r 0) c r0 •o r0 -v O Mn O .n L L L =0 � r- r• O) d Ln N ]C c,7 2 H V) 2 -7 rn m -� m -D m J Ln C W Li S Ln qt 0) O O O •r) •r a0 6 0) ma O m Y m m LL (n) S cm G lid C C N > N � C � N E L c L N O N > � G n N C O ea ro V) N a O c 0 N C7 ro w d s N X 9CX Y L O y U > O41 F 4- U 4- O ro O L L U Y 4- Y ro c ro L N y y 24 r N 41 N L 4•. N N O OJ _ Y = 4� r C 414� a N >y d C m C7 N N (n N 4. V- N ■ X O U v 01 O) c } }}} S. }} In Ln N N Y N vn m 3 3Lo CD N = r., 0 O1 N r r t� C Cr LO CD CLO z r 00 C r r • f d i L 0 O C a O E O 0, N 4-4-. O Nr C N r E 4- L ro ro ro •P v0 O O 4- J J Y Z LL N n � � ro v N L N 4� n 3 4J rn � r N _ � E ro c ro ro � N N N N W M v E E E ro d 1 '00) Y r r r Or as a O O c i s- o, E C O. d Ofl m C r N 0 0 0 'mac Eccc Ln X X X X X X X O U 4- ro i Y N C v N L w > } > y > > o L/ 1 N Ln N N N N m A r O II C LO C'O N C C •� .� y O ^ Ln X O C d N E LO N H N O C t0 C drsE atv Y R a Y�v YY O Y U TLC LO • e0 ,O N t r O O Y Y C tl S- W W 4- E ^E c QJ Y Y L f6 1Q O. i0 o C L OG d 4j 4- SOT ci L LT R Y to G1 7 r >O N C OIr N Of O) N C r N L -n c LD N C C 3 r O Y N .- C I� GJ •� •r- O N N O L Y O U U•.-\ L U Uri Y Y O_•'- V �"'� In G•� •r N •'- C N A L d W Ln O c Ln Lnt Er u O O •0 E °�- O 7 Y r0 u •r d C) 00 t0 O O 0 0 N LO LV .^° C Y O Y Cu., O. L L O a S- L L> y C X U >> U Y L V U.G7 C L (0, C F-• ) N O N� n n> N N W r-. U �••• Y J �•-� X N T Y N � Y � Y `h 6 L Y •7 ?C YC X X G1 U I C � C r+ C y i L O N Y > L W Y LO N L Y m E a � Y G W L d C N 3 Y 7V)7?>>>7>V)>V)9 O N V; N to L/i fn V7 to Ln ti) to N v] t✓7 3 3 z 3 3 O f O 3 o3 • LLn . C i Y g C) i N3 M Ln d LD O Ln N M O M d C° O N O) O 1 0 Ln M w L LO rLOM N C tD N11r > V O L O RN V rt1r Lr Lh O O 4- O LF > N H to t 7 d:m W M4•• Y 6J •r Y •e- LO 7 O° rO a •r• •O °n C +� C C N N Y Y �Q N Y •ry N O) 0 L N Y R7 LO L6 L6 r e- L Y Y N L L CM N L X Y -O a'f C7 Y V) Y h C Q J C 2 J f C -j tl d Y > W d L CU y > > N Y C N J f FOLLOW-UP QUESTIONNAIRE 1. The name of the power station and of the river. 2. (a) The owner of the station. (b) The company responsible for the management. 3. The year when the present management took over. 4. The water head m. 5. The volume of water that goes through. (a) Max --- m3/s (b) The average water volume on a winter day m3/s. (c) The average water volume on a winter night m3/s. 6. Turbines (a)- Type (b) How many 7. Give the type of regulating (day and night, week regulating, etc.) and the variations of water level upstream of the power plant. 8. Type of outlet. 9. Are there current generators, air vents or similar structures installed in front of the dams or at the outlet; describe. 10. Type of water intake. (a) Open canal (b) River power plant (special intake canal missing) (c) Intake in connection of reservoir (directly or via a tunnel) 11. The dimensions of the intake canal. (a) length-m (b) breadth-m (c) depth 12. The dimensions of the river upstream of the power plant (river power plant). (a) breadth-m (b) depth-m 13. Water intake (a) Now many intake openings (b) The breadth and height of the intake openings 14. Gates: built in (a) Yes (b) No 15. Gates leaning out from vertical plane (0* for vertical gates). 16. The dimensions of vertical racks (a) diameter--mm (b) approximate separation--mm 17. Are there mechanical gate cleaners. (a) Yes (b) No 18. The heating of gates. (a) Lacking (b) The fraction of total gate area which is heated (example 2/3) (c) The entire gate area is warmed 19. The type of warming the gates have (a) Induction (b) Circulation of warm water (c) Other 20. The power (electric) on the gates (a) Total --kw (b) Per facing surface 21. Temperature and observing the ice formation. Reading slush term (0 meter) (kvicksilver term) ------------ times a day. 22. Re" a� temperatures recorded automatically. (a) Yes (b) No 23. If the account of recorded temperatures is missing, how often is the temperature observed when there is a risk of icing. -------- times a day. 24. Is there installed a meter over the gates for measuring the loss of the fall. (a) yes (b)- no 25. Are cables (chains?), cords etc., used for detecting the beginning of icing. (a) yes (b) no 26. At what temperature are the gate warmers switched on. ------ °C 27. Are the gatewarmers switched on manually or automatically. (a) manually (b) automatically 28. During winter, the areas which are not covered by ice, directly upstream of the station (including the intake canal), cracks (rifts?) (a) length ----- m (b) breadth ---- m - 29. Stretches of rapids (streams?) upstream. (a) distance from the station -------- m (b) the cracks length ---------------- m (c) the cracks breadth---------------- m 30. Stretches of rapids (stream?) downstream (incoming) minor rivers upstream of the station. 31. Minor river ------- m3/s 32. Temperature in relation to the main river. (a) same (b) colder (c) warmer 33. Measures to hasten freeze-up upstream of the station. Laying out 4 ice booms. (a) yes (b) no 34. Reduction of water flow during the freeze-up. (a) yes (b) no 35. Is the surface of water constant upstream during freeze-up. (a) yes (b) no 36. Other 37. Occurrence of icing on gates. (a) yes (b) no 38. Tracks or turbine blades (a) yes (b) no 39. Give the time (morning, daytime, evening, night) and the type of weather (temperature, precipitation, the direction of wind and the wind velocity, etc.) when ice forming usually happens. Also inform in what direction the intake canal is (example North -South). 40. How often does the icing occur. 41. Estimate the average production drop per year during the last 10 years ----------- kWh/a year 42. Has icing occurred in spite of that the gate warming was switched on before ice forming started. (a) yes (b) no 43. Occurrence of ice dams. Where do the ice dams originate (give likely causes). 44. Consequences of ice dams (reduced height of fall, flooding, etc.). 45. How often do ice dams appear. 46. Estimate the average drop in production during the last ten years ----=-------- kWh/a year 47. Ice floes against intake canals. Interferences in running (extra work input) ----------- man hours/a year 48. When do the ice dams appear. 49. Is there an ice outlet. 50. Does the ice outlet (isutskov) work. (a) yes (b) no (c) partly 51. How is ice removed if the ice outlet does.not function. 52. Estimate the average drop in production per year during the last 10 years (because of the decreased waterflow)-------- kWh/a year. 53. Other ice problems (ice pressure against dams, ice on gates and gate folds, etc.). 54. Give the effective methods of fighting against ice. 55. More information (the rest) (for example: details of the form of the station; details which have importance on the occurring ice problems). EXPLANATION FOR TABLE COLUMNS 1. Type of intake. 2. Average water flow (m3/s) and corresponding velocity (m/s) upstream of the intake (winter day). 3. Average water flow (m3/s) and corresponding velocity (m/s) upstream of the intake (winter night). 4. Size of the areas which are not covered by ice and are directly upstream of the station (m2). 5. The heating of gates. The total warmed gate area, (m2). 6. The heating of gates. Total power (Kw). 7. (°C) temeperature at which gate warming is switched on (a) auto- matically or (m) manually. 8. Icing. Tracks (?) 9. Icing. Gates. 10. Icing. Drop in production (MWh/a year). 11. Ice dams. 12. Ice dams. Production drop (MWh/a year). 13. Ice floes against the gates. Reason, point of time. 14. Ice floes against the gates. Work in man hours. 15. Ice floes against the gates. Production drop (MWh/a year). 16. Other. MA e• .O ^, 1 I 1 1 I N 10 I la 4J >I Y 1 1 I I tT N y L r L i > Z Y I Y I I O I 1.Vy a••I L N N I I I i T� CJ I 1 I 1 O 1 1 1 d 1� r cn I I I I 0 0 O I 1 1 1 y C X � 1 I 1 1 > 1a I z n e c 'o o CD E o 0 o a•-� N Z Z O Z I 1 1 1 b 1 1 i I N 'r °r• ° p t I t 1 O 1a N O r In N In 1 1 1 i 1 .0 'a 4] 41 > O r 1 I 1a w i 4J 1.1 Y L > r 41 O j I 1 1 4] v i •r- QJ C 41 X t vE x 4. O N = N Y Y 4] C 4J I{! o a E i 1 1 CD1 +t 1a M .. L N Y a v z i Z i i o ' U a c v- v- r''aa I ^-� N >l 0 v1 ra 0 0 i t41 1a • ° In I 1 1 I O 1a w V N Ln i I 1 I 1 d q 41 O I t B N t 41 Y Cm> E > > X r p j 1 1 1 v L •r• v C 4] .-� 1 I I >a L d L •r- L N E O X X O a..• O N N 1 C N O a) o.N O 4] 4) 1 I 1 O I 1a (] Z t"7 f cn y v Y > 1 Z 1 Q (U U c It.. m 4•. M1] 1 1 N >l O •.- o 3 0 O r CO I 1 1 I R N I O °_ 1 1 t LC) I 1 1 1 1 d b N a O 4] X 1 1 1a N .G QJ y p1 O 1 1 O fa 1 I 1 cm E L dL •r- L N O Y CD w 1 1 CDb R O 4) I 1 I I O 1 a.l 10 V L f" E F rO 4-' N v Z Y I 1 1 r 4] V C '1a 4-o 1 1 Cl I ?� o N •�- 0 3 o GI 0 r N O I rt1 • I 1 1 1 1 w O N 1 I 1 1 I I S r R E ? X r 1 1 t 1 I 1 L W C 1 I I I I tT L dL •.� IT In X In Y I N N I 1 I 1 t I O O N a.a E n B n O 1 Q7 1 O 1 I I 1 r U L v C' i-I I Y Y I Z 1 1 i i V U Cr 4-• ro L I 1 1 rr O 0 3 +•t �rl m In C r I I 1 1 1 1 I 1 i V • 1 I I I 1 I I 1 I 1 t I b 17 4J c O I I I 1 1 1 1 I 1 1 L 4J 44 (V O1 3 1 r 1 I 1 I 1 I 1 I 1 C •�- aJ � CDI 1 a I I 1 I I I I I 01 41 d 1a •ate In X 1 N 1 1 w N 1 1 I I C d N t T E 1 O I I 4) N 1 I 1 1 I I •V d L � 1 a.a I 1 Y Y I 1 1 1 1 1a C 4- ra . 1 I 1 1 �--� L •� 0 3 O co C L ^ p 1 t I 1 I 1 I � ^ CD I 1 1 I � tJ1 I I 1 1 1 A N x } ti I 1 1 1 1 1 1 I 1 1 1 1 1 a� N X O N 1 I 1 I C 41 16 to tVa 1 I t I I 1 I 1 1 co 0) O N M In lc r- L N N N 0. r r I 1 1 r3 C 1 1 1 L Y O O X T 1 I I L N G > 10 1 1 1 I O U 1 1 1 1 I 1 r CC C n^ to q 1 1 1 r r r r Z Z Z 1 1 1 1 N L N tm 3 1 O oC O 1 1 1 T L G Y N N- M tT• 1 1 I 4eT� V L Y O O X ••- 1 1 1 TU G O 1 1 1 1 G i> >_ try > m 1 1 O U N N 1 1 1 I N O N tQ O 10 N N 1 t 1 i•1 C Gr r r a r 1 z 1 I 1 o.a a R 44 r L � Im ^ 1 1 1 I 1 1 1 I U 1 I I e L `y T9(D 3VL- O lL L G Y O E W A ' 1 t 1 1 1 1 I I IT It N C +• _ .�— O I 1 I 1 1 I I I 1 3 � T O ^t0 O 1 1 I 1 I 1 I 1 1 1 TN O Nro1 L '� c a' r Ct 1 1 1 1 1 rU.1 0 3 d Y 10 G 01 w�� G r p 1 1 1 1 1 N • c •, � r r I tT - 1 1 1 t 1 1 1 1 o u 1 c 1 1 1 I e _ Y E OV OV 1 Y 1 1 1 1 t 1 I •, 1 I 1 U 1 1 1 I 1 Y N 1 1 1 10 O O 1 1 1 1 1 . to M N N 1 r Z Z 1 1 1 1 1 N t s 10 - d L Y U >tU C i.I N N N I c 1 L r o L Y 0 0 O L 7 a T Y I I T 3 1 3 a)N CD i N G I 1 4+ N O Y cn > rt3 1 1 3.G U N N O 1 C a M M m. G� r (LOj i 03 C q V 1-V• 3 L N N N C CD O I I CD cm1 1 1 L i> O O I 1 1 N G 1 1 X Y I I 1 1 > O 0 1 1 I 1 r N rt7 L L - C O / 1 1 1 a- 1 i a U •Ct ct p1 I 1 1 U N a X C 1 I I I O d 1 1 O Y 1 I I 1 1 O c N O 1 1 CD U y 1 1 I 1 C N N Cr r O N 1 O 1 1 T Y N 1 I 1 1 I 'F L 1 rU-i O 0. > 3 O O 7 0 1 1 a c r. ^ 1 I - 1 1 1 U N Y G 1 I 1 O d 1 O O X T 1 1 1 1 I QL1 Y Y 1 1 c N > I 00 H N 1 I 1 1 1 a N N ;.1J ^^ ^ 1 1 1 T Y N C to r4- •"' N M to to I� m p) O r N M i r E ro C p n� ro - � E O O X r 1 I Y 1 1 O V N w - iO ( co co co r Y Y N 2 O co ^ O i O O C ' X i ro y m 0 � CD 1 U i Y Y Y 1 >f r 1 C Y _ ro V Y C ro V tp Lo CO ro O Q O CD ' E 0 0 X m CJ Y Cl 1 1 V N O I 1 p m p jJ iO O O O O a) cm C) w N Y N O 0 0 S- L t to \ t.C) Cl 61 •.• N i m N r Y 'H Rf 1 e t t iro C p v O O X E o 0 00 0 17 Cl) N N i I C � U co m O O O E I X On I N OCD 0 1 O IJ O cu N L r C C) I- L N L CI r a Y ^ U 41 Y Y Y o '1- = C — 3 O ro O r I C r W a) 1 a) a) al r C L CJ N w t ro +- C 3 C al 01 7 L Y w CJ 3 •p ro 1 L c O Y �>IN I C) ro L Ih C� C Y 1 L L CJ L O M ear •r ro Y U J 1 L s- ea 41 WCw ar^a) 1 C e0 3N C a!Y (a C� 1 a% L E w C O r 1 L w C C O V C) ro +• U 1 3 'r r •r It'7 r . ro Y 3 V Y {O•i ro 1 A OI 1 1 1 aJ N >t nO L• 1 1 1 1 CDro T 1 1 i 1 3 L ro v L 1 1 1 1 IM w a) a) aJ i It 1 1 C C) L L Y C U-0 1 I 1 1 •ry E ea Cl N C C r I t 1 1 U •.- () L cu •r C? 0 I 1 1 I r•-e Y>, Y 3 3 1 V a) d O w N E m L+ w ro Y Y Z L U Y C t C 0.7 0 � r L N C r O CJ aJ Y C m Ar to ro O to O o 1 1 1 N Y m. 1 I 1 i 1 1 1_ 4e O O w •t- O 1 i 1 1 O i C 0 1 1 + C Y 1 ro 1 1 r ro V •P 3 V 6 d 4 O X t > C T i i F- QI O CJ •t•- � 1 1 G � i i tT us O 1 i C N Ln O E aaj awi 1 aNi i 7 Y Y >- I Y i O 4-. O O N O ^ o ro o w Y aJ L 0 0 a m w at L N C > t > ro O+-r 1 O V Y 1 U O E eo L 1 a i al 1 aJ O L Y nY>ro 1 .�. 7 N > � C/ 1 3 y w ro L r c o to Y C C C as roa ro al > C U 3 N•.� C� OfY ro r C ro C o C s-; o > > 7 a O C Y > I ro T N L L IT L 0 a/ 0 +cYE EH CD N r LO tiY4J y 3 w 1 � w EY.0 c c •.- ro 'O C1 L 3 Y ro L� Y C Y N I-> C N O C Gr, O c T E V Y a 3 a d r roc Y drnro T T L d Z V C ro V V 7 L N •.•• ro r ro 1 C C U a ar me A 4-ro C N r •r Cl •r- w Q) 4! O C GY 0 r W U ro w= ro 4- cn L ro V O y ro a U a ro C1 L o O 61 IL c O U C N O C to r--t ro ro V- ro 1 6a Y { N N L 91 O O 1 I 3 1 I O r r 1 1 O. OI 1 1 0 0 C 1 1 L 1.1 O Y 1 1 C1 C 1 I Y 1 1 > R 0 0 X V N Y r r r R O N O 1 1 N N L 0 0 3 1 t 1 1 I O r r OCL 1 cm 0 0 C I 1 1 1 L 4•1 X 1 I e a > > R 00 O U VI 1 t I j i C d N N r Z } 1 I I 1 1 L ' v 3 . O_ N N 1 I 10 O 1 B i 4J > R O O O O R O 0 O 1 1 C d N N r Z } O Z O 1 N L 3 of to 1 t G N r 1 ! + E O O) I 7 0 0 O N C 1 1 1 i 41 I I I 1 1 1J c •c X Y >R oo u w i O ^ O d 1 O 1 1 1 C d N i N z 1 i 1 1 i N O Cl) I N r EO rn i O O N 1 1 1 L Y E T Ol C 1 1 Y X > R O O C V 1 1 1 I C CD N 1 O I 1 I C' (n N i N = r 1 Z 1 1 a o RLOC u 0 0 w m I 1 Ol Y r N O W U N 1 p > 6 r r C 1 L N CJ CM cm� I 1 1 1 Y N W CDtf'J In O V N 1 I R +.J R i O N N L A R O Q J 1 1 C C L R t i V 3 4••• +� C�� I E O 1 dd L O1 C a O I d c L I R V r O^ 1 i R G E O O I I E X 'WO.r- O C I i S U N u Sri I 1 1 O O E 6 O O L +^•• Y c M r O L R U ++ R O V I GJ CD 1 t i C 01 N N ct 4•- 1 3 R Y i r N M V () 0 n co O O N M C Y 1 1 C1 1 1 U 1 1 % 1 1 O N 1 1 L 1 I O) R 1 1 C Cl 1 1 U 1 I � M r v C r C C Y S- R R i ; •� CJ d O v R I 01lO 3 Y B O i O 1 W> r 6 C 1 C L I v Y U 1 Y d /v I U> Y 00 i 1 O •ti O R U M O. G vQ) O E L L C + o B O (AY R N >1 H GY QIO O d Y C r r r C U Y r r QI C O O •V b U 3 3 6 cm 1 c r c i+ U Y L v w Y IA O C Y N i C Oal R 3 M c E OJ L Rw OY >,>, c(n L + + C r R O C L V 3 O O 1✓) N M C7 ...m Z O) 1 C t C UY U 1 1 Y N O C •r 1 1 N L C• 1 i cEai;�ow R\ 1 1 N N M 7 W-U.•Z R d u C O O C L Y R V d o-•r >I rn In QJ Y V U 0J O Q C O C R N i Y rn v i d O v R O O u C � .N•- L d ' O r r 1 1 d • Of I 1 OV OV L 4J t 1 O N C t I Y It 1 E > R 1 1 X U 1 1 J N C d N O r 2 Z 1 Z O i O O N N C 3r O R C O N N C +J L f. O N 3 0C R a) CD O 1 1- 1 1 U O a) O O> iJ L C O i O r r O 1 I O N ro 0 7, C N O d • e tT 1 1 00 1 1 R C N R R N eti R L J 4 1 1 X T +- L L 'n N C 1 1 .Y 1 / / CI R U v CJ • N O s- > R O_ O O U N 1 1 C C/ y o dr r r N L O Y •'•' r O O CJ / O 1 i • r U> I O E r R r i R N N � ero Z >- 1 Z O 1 O U C r tL N tT R R 7 LH O ... O ch I +J 3 .... CI T E C V 0, R C/ O t > r O i % R C N C 4J eJ U N Q r N .0E ea s L v , N 1 L 3 1 1 1 1 I i W i•f R OdL0 1 1 1 1 �••1 R yJ Iv CD C; 1 j 1 1 1 1 1 d C 'U L 1 Y 1 1 1 1 1 i FJ Y R i I i R 1 1 1 OCo r O > 1 3 -C '^ C Ln tLi N lV 1 1 O iJ L® r E O Z d Y 1 1 1 1 V= di Cl 1 1 I O 1 f--f .0 0 iJ G a C O O O R +J N L O L R N ea iJ .0Y Q N 1J C 0 1 L d V C y r R a '^ E 4� C V R 7 L L +J ^ 09 O 7 tT > C C > r 1 R C O •r R R I 1 U I I E 1 L U d d m I 1 1 1 r a•J I 1 U •p d O 6 V N N 1 1 1 I 1 1 1 1 ro a i v 1 E C) O 1 c •r 1 1 N 1 1 0 r+ 0 0 ci I 1 L I 1 J C'4..1 N Ln ta'! LC1 ( R O Q/ N 1 1 Y 1 W N N N 1 r = Y r 1 1 Q 1 1 rr d•J •r }J tH d U Y C O R r0 O ci a L S- N E R V1 R N O I + O N .0 N 41 V L 0 iJ L v R L 0.4- 6/ iJ R O d r d N L CL O Gl 11- E R L O W ro N Qj N O (v •� L O r r 1 1 N C d • O 1 1 1 I T O U N O O C I I 1 1 1 1 N 1 1 1 C r jJ dC L O 1 1 1 L 1 1 1 N C)7 0 U L L I UO O > I 1 i t X Y U I 1 t Y I 1 •.- 4 J E y 0 Z O 1 1 t LJ I 1 1 1 y N Y R O L - C 10 O r IY/f N L � L O Y O O O 1 i 1 I 1 U C d r r IT CD f f 1 I i-t U to 1 I 1 1 I N O L Y V V _ iJ 1 1 1 1 I 1 R p O C f I Y 1 1 f 1 f 1 1 .— Or C - > > ro I 1 CD CD'a x U o N aJ O a I 1 1 1 1 1 N v m f C:: O• lt] LD O I t t f r C a c Y O'r E r O Ll7 r O O Y L p r ro m a> N V C 0)X rr L. Co Lo O r i i T N Ewa=+ c 4- U O N ro .4) 3 0 O C L Y 3 a I az+110 = a+ t 4. C t0 y ro 47 T %- r d M1 L r C 7 OJ N 0/ (J Y c >o NaL.1 enY c IV 4- O O ro d T 4J N L 07 0 N i L >1 a%NLY alro i, L07++roci I _ 1 N r g m d 1 f d > 3 4..1 Y d 4J N ro 1 1 4J O E 10 L 1 I L w O L 01r i 1 O r Y. 4-1 N y dr O O 1 ' Ta-1T�E3ET z Z O 1 O ' 1 O L= ro 0 0 O v a+ O C 0. y Q d+.> N 1 1 1 1 1 1 1 1 1 1 1 t 1 / I e 1 1 B 1 I 1 Z Z O Z 1 1 O 1 N L >� ro a � ar u`+ E 4J. s a+ c ro > c > a ro n E c a+ L oj Y ro O 3 L ' 1 N L O rr 4J N 1 V d e . Ia d U 4- O O C 1 O V L L O 1 O Ol O > f i X u i Ol T O O CO tc'1 O ro r O Z 31 Y Y O 1 1 T ro U U L t0 l0 O Z O ! L OJ 1 1 3 1 1 I 1 d 1r cm 1 1 I L 41 py a p 1 G/ C I Y 1 I > m i 1 X ! T r O O ro O O O 1 C' M. m to O r Z 2 O Z O 1 O I Si N 3 I 1 d O 1 1 1 1 i a.1 C Ov ' p c 1 i w C 1 1 Y 1 1 1 > ro 1 1 X U ' I 0 0 ^ O O O 1 1 C d ccO i N O N N 1 I i d I � i. 1 I N ro GI N C 1 I 1 Y 1 I OI T > ro O O 1 U N 1 I C z n r ! O N ^ O z r z O I O 1 rr M r N c•'i c7 m Ql O 04 f•7 -Izr ti7 to L 00O (2) 71b � O p > N O i y Ci L i C d 3 R m 41 . U r d O C p O C CJ L to C T N cn: i 1 CJ .�- O R L 3 C;i 1 = C CJ o 1 1 +' 4. d 1 I L r r ^ Im C R ea O R L >> .r 1 1 0 7LM O X Y I 1 C 09 R O 4JR r X X 2 C Lo CO Z Z Z t 1 U C f.7 ® O F4 T T R m C ci LO 0 i N I O L lA 1 i C) 4— O I LOW L 4 1 r 3 O 1 E N i I L 4% M R d r I IT 1 i.l T C n R T 1 N >I .0 L N 41 O >> 1 1 C 1 Cl U N E 1 I TC OI C T ; C R R > R X x L V 44 1 C C N N I d R R 1 R O C) 1 ) N 1 T L >t CD�-r C R d� N i O i-Ve E 4J C X 1 1 C 1 41 R 1 L is CI T CJ C N Y L O 1 1 I t I 1 1 1 Cl B L 41 d C C O N 7 O C1 R R CJ 4.t N t,() O X 1 1 t 1 I 1 L T N y i-1 C C7 1 1 1 1 1 O 0, 4 a 4 a E N rs L mot' N O t 1 1 1 1 O Z Z O 1 T Q) R r d C L 1 Z I ea 1•-e Qi d R to Y A C C/ N C L 1 1 1 1 t 4J -0 R C CI > 1 1 1 1 1 cm R N O C T 1 1 t 1 r C N r r R CA 4J 0O -Nd1 1 1 1 I T T O X V t I L N i ti i o LO d'tn CDr 1 1 Z Z 1 1 1 U 7 C d d 41 1 Y N O 1 1 1 1 A N d tT 1 1 1 t 1 N Q) 1 C 1 1 1 1 1 C) S- S- 44 Y C t Y 1 1 1 1 1 L i > > R 0 0 1 V 1 I 1 1 1 4. ec 1 1 1 Rm C N 1 1 I 1 1 4Jd I I 1 1 ¢ 1 O N m O C L O I 1 1 4) N >E r. O O1 i C 1 I t i i 1 i L >> R C! CDR r 0 I U 1 1 C C R R R 1 I I I 1 1 I t 1 I 1 1 I 1 1 C W 4- 4J 1 N C O 0, t E O= E N L O Gt O O O C co a N d CLJ C I j Y y E -0 R C CD N > R 0 X v _ cu r 1 T CL LOLL O O N Cl O E L Cl i C/ r t\ 7 r O r Z Z O O d Y 4- O r N M V t0 10 n O Ot p .r- r N L E Cf O y w II WMA y O r OI a > o YL > r eaN x y i•1 N X V C iv R O R O L- E r Z Z Z i X O ci r N N O 1 v t2y 1 w 0 0 L y-a I 00 CI C 1 1 > 'fl 0 0 O I t BI r co r 0 0 1 O CC O. r p7 1 Z Z Z O O E ro v T L E R O L w C C 0 0 R N r LC) d >O U O. N C 1 1 C1 w Y L 1 O O y i t I41 O _ 1 +.o V C 1 UIV R 41 t79 1 I X y 1 C Y 1 1 1 = O C1 Y— C is w ro C p 1 U a)QI 1 1 V L C O Y O O R O N 1 Cl w 1 1 U C O ro CY LY� 10 t0 r Z � 1 > � 1 U 00 G N p Gl 1 1 0 C L O 100 1 1 y C 4•- Y Y a . t 1 ro 0 0 m O .° J 1 1 CJ-/Ja C O> OQj N y i> >> a R Gl > C iJ N Y R X Y B 1 O R 4Jis C v i r 0 1 1 U ro a w 3 nw i r. a 1 r Z Z 1 1 J 7 C w O r C O O Y a LO O Y 7 C i I 1 I 1 1 1 U I r °P O N 1 1 1 1 0 1 lU y r r 1 I 1 I 1 1 1 C R CJ > O Cd O 1 1 1 1 1 1 C N y i >> 1 1 I 1 1 1 O L 4J N X X X 1 1 1 I 1 1 I U.0 ro R O I 1 1 O O f 1 1 1 1 C w p N _ + c C i R o i0O ® a �' y CC •r ^"' i ° ° N 1 1 I 1 1 1 1 I Ny V N O O m a r 1 1 1 I 1 I 1 1 1 C ro C L O 1 I 1 I 1 C N O r i• }J Y L, L. ro Cl L. d >> X I 1 i y V^^ o Y 3 X X 1 1 1 1 1 1 i 1 U.0 C U W a S-U' d. ro R O 1 I 1 I 1 1 f 1 1 1 1 1 1 1 1 1 Y^ G O R V •+•VV}{11 i Y 'C i .%. E .E ct 1 1 1 a w v .••. � L i R •.L- N w O N 3 C> t O> y R O 0, L r r 1 I L ro i y> a.... d . I 1 d E C a R o w CIOt a 0 LC i i- 1 Y I Y dam. ro 7 C Re &- cia)> C L 1 t 1 y 1 N i a 3 Y •r CD > R O O 1 V I 1 D. N R r C CD t0 1 R Cl 47 0) N LR C aY •� Irl 1 r Z Z Z CD L L L G 1 L Y CD CDI- ++.r ++ a++ N� L 61 0 O ^• ^ I 1 a 1 m I Oj CD j I L Y •rC- 1 C 1 1 1 y t > ro > 1 r O to 1 R O O O 1 d' a CCl) 1 r Z Z Z O 1 O O M yr Ln b t\ tb m O r r N r M r Cr r Ln r l0 r- O 11Y O C 4� M r 1 N R R N G)O i i 1 i I i 1 E 0. wNr T i i 1 yt 4) o O N S r r I LT I i I 1 4J � 4- y }" >> i I 1 1 e I y T Z R Q I y i i 1 1 1�� M C) 4J c NN R A 1 i 1 , C A R Y Y 0 ; 1 m i 1 r I 1 1 1 rUs N N r c 1 r C t c +^ L 4a 3 41 R C L O4� YY O G) r L eEy R �'f N T Y R Ql I O S. r Y _ e N 3 p Ia O L 1 1 1 r R N C i N—= C) L C %1 N I 1 1 1 (D C O Cl E V .1L 1 p Y N I ? 6) R N L L R E R rn 4R) w I x Y t e p i U 0.4» R Q) eam M m o O be C) 1 I N a)=vavLE M r O r Z Z W °O CD M i 1 _a, C t.T R U C U 4J 1 S Y R U Y R Y N N L cu a cL N + Y r i = O 1 07 1 1 1 1 t a) s p I e 1 1 1 1 e L C) 7 L Y N 1 U 1 1 1 1 1 O� CO ko 1 10 1 1 1 ne-a i co S- (V 3 o c L N r Y Y j I 1 1 1 1 i� O O 1 IT I I 1 1 1 (1)L 1 1 1 4-GI I 1 Y I I I 1 1 L L R C) I 1 Y 1 1 1 I 1 7'o 4) N 1 1 t 1 1 .'y d c N CO 10 1 R 1 1 L N c0 1 r = Z 9 1 1 1 1 1 1 CD U O C Ot O 1 IT I 1 1 1 1 R� N 67 ' I 1 1 I I L U i N y i 1 1 Y 1 1 I 1 I L R L R N 1 / Y I 1 1 1 1 tT L y L J-1 N C' 1 U 1 1 c= ?. _ 1 1 I ^• L N CO1 r = Z I I 1 1 1 V U Q 1 1 1 S O N e i L O . 1 1 1 1 O O 1 tT 1 IH Y L i •.e- 1 I I 1 1 1 1 1 C 4.1 N /v G Iv !9 10 1 Ip 1 9 1 1 1 1 f1 •... L N W 1 r Z Z 1 1 I 1 1 1 N N S- S- r 1 • R C 0 N N O L C R 67 R tU C .r A L V t[] M R i N V t O R ro N- U I O C Y R tD L 0 0 C i i i M V tL O 1 C) c L O I I Y 1 1 L o Ia -VL- L O 0 0 0 x U N I O r V ^ r N L)'Y N N N r O R r O Y 1 1 1 41 0 v N Z Y 1 1 I =Y .U-. CJ Cn O.0 r ^ N M er CO Cf 0 04 ^ C ^ b q Y > i T C G > r T L q L y C 6 N q N y U L G G Y O O^ N G y r• 1 0 O tT y E N Y E; Q T 41 4j 4J O O O U y O LO G 1r /O •) N a r• Z Z 2 tN c y '••- 1 1 N T Y O 1 O I 1 y N y> V 1 C 1 E O T 1 T QI G m w 1 1 .NC 1 +•> c 0 1 X U i E i 4J N N n•• i 1 O 2 = Z 1 H N r- y N 3 N y i r T 1 1 1 1 1 T fT O O O tT 1 1 1 1 y O y L I O1 m c q y I O y 1 1 1 1 1 Y T T 7 i 4Y N 00 X U 1 1 1 1 1 1 y EEy O c y i P^ q O O 1 1 1 1 i E q I"1 N O' r• Z Z I 1 I 1 1 U) y G G N L T y y S- u > 31 L c L d M i 4- q _ I O O O N U E Y > g M V E O q i t c 0* V y y + y V L N i q L 1 y 1 I P• q q Y S. I N 1 1 1 y y L N t y 1 C 4J 4• B G 1 y�•rq Y • 1 1 1 a U 0 3 3 ci O O O Y I CD C i c w- Oy 4- q 1 Y _ 1 1 .N.. 1 1 1 1 1 r• y O Y Y N d•1 O O X U 1 1 I 1 •- N L O C G ^ O q O O 1 0 1 1 1 I 1 y 7 U g N H O ^ 1 1 1 3 g c 1 4> 0 O - N O y c G L M O y N O i 1 G I I 1 I 1 1 M 4. O Y >> I T 1 I I 1 1 1 y Y q y 1 Y 1 1 1 I 1 1 Of y y. 1•a q $� � 4j0 X X 1 V 1 1 1 1 I 1 1 G U C L i` q 0] C S- q q I ^ O O 1 1 1 1 •r C •+ L g g 1 Z Z 1 1 I 1 4J CL C 1 I •-Vi N f1 N O O r N O y M w TL y Y . • 1 I 1 1 1 1 1 i4E• -1 y C p y> O Cl I 1 c 1 I 1 1 1 1 1 i t p 4J -0CM F> t Y 1 1 y Y Y IC tT- y 4j N y X X I U N N I I 1 1 1 1 C V C i i © i q q f 1 q y y 1 1 Y q q IV 1 r• Y Y 1 I 1 i C I •--� N CL N O n co 01 O N M tPl t0 i 1 O j i 1 N 1 N L 1 1 C j T r r 1 d O O I C Y V V 1 T ro 1 i 1 Y I 1 1 1 V C M M 1 ro N N 1 � I I 1 I 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 i 1 I 1 t Z 1 1 1 w O tU N ro V GJ a L d U p T L y_• N O N 1 1 L 1 I M +1 C I L ro T 7 m �v Ov 1 I N m IA`. c I r•- T m 1 U t r T +1 O O X v b C ^� ^ O O O 1 LO t0 CO 01 CD I L T O L QI N Y i C r r T p p 4! Y B I O T J+ O O X U C In In R O O N C R C L 3 O f+9 en CL tT O O C L 1 1 p T 1 1 Y > > 00 % w T I[1 t0 (a R L 61 > O O Lfl Ll•) 1 1 1 1 1 1 I 1 1 1 1 1 1 O 1 Z I O s N Q1 Y L Y O > L I 1 Q7 I 1 N 1 i L N N 1 I 1 O O 1 1 C 1 1 1 1 T r r c 1 1 I 1 d op 1 1 Y 1 1 O T 1 1 R 1 I Y 1 1 O O X U N l 1 a I2Y 11'i R O Ql r O I 1 r N (`7 V lY9 <D n W 01 O N cn 01 G Cl N C O C LT L TO 4J 4J Lt R N H L U C) R U Q1 CIo Oea N O. O N T O T r O N C I Y cmcu L 9/ Cl » Lr cc � A •W R as o r o > r N (al r tT .Q L tp R C =O . CI N 01 4- T L •n C C r 3 5 C1 4• T R O E 41 .12 rs O.T N V R> w E a d U L r 1 T CL Cl +J 0 O N C8 V.. L C1 R E 01 OO R. L R C C U C) W •r- T 01 L O U N i. T4- L-- W O N N O Y C C L C R 0 09� U O 4 6 R N L O ae > L Y R O O ITT C) R Cl T C N C Y L C N E 4J c T u a R X O C R 11 C -0 C) U O C1 7 C1 N U C N C R C T r O 1 I N dT 1 1 id 1 1 N N R 1 1 C L i 1 i E i� Q1 1 1 O 1 1 1 C O 1 1 Q 1 I N L C 1 1 3 C1 I 1 Q1 4-1 I C r L I n R 1 1 U01 Y i 1 tr r N � Ifl t0 Y I(7 O C > ro i U a w OI a C L O C r a In N a O Y n Y C N CC O lt1 X CY.1 w N G T > to (n In ro O O O a T O X M r r r Z Z Z p a r c L E ro T Y G ® ro O a> vG O Y L N I O Y C ro a tT 1 L.1 N ro C p �( U w41 0 C a Y O O a Y In p r = Z p a l� Y T a > a. N � y O Y L r \\V. a ro ' G ro a x I y N iJ N O G 0 0 N c) w w . C a bl ro to ^ O a ro a '"' O) U r CV r Z Z N Z p a F-` Y 1- a c is r •"' a L 10 ro ro ro a G G G 1 1 Cn Ow ro G ON On 1 i r d U 7 3 ro V 1 1 a 0 a 1 1 G L r X Y e0 G ITY 1 E I E L) Y G U 1 1 a CD Cl N 0 0 O O C a .-.CD U a i 1 1 C U L O M O Yi d ^L ++ O 1 1 IT c eq Nd° M O N�++> w Z Z I 1 N C a O O Y r ei- L }roi dy Y N U > a L a N a L 1 In ro Y G 1 1 1 1 I 1 N U t a s a Y 1 1 1 O Y a OV 1 1 1 I i C° CDa ro 1 Y 1 1 1 a 01 d ys O X V 1 1 1 r iJ d ti M ro 1 1 1 N O M O In to > L T O . > r i ro a C d w ro a U 1 L 1 Y 1 01 1 ro 1 C Y O N 1 L i r ly IT 1 a 1 a Y E r 1 0 I 1 Y 1 +� O O M 01 to l0 t\ co 0) O N M U to t0 C N N � Y N ro C is C S= Y d i C •r C O � L C/ O N C! Y i N C T ro Y Y � X Y C CD to + mo L C Y cm V V G Y Y Y X Y C Co �r O m i Y O i C! N 41 L c v o0 0 Y i 1 1 i X Y Jro O [1 N m � r N x X CJ Y O C Vf d = 2 Z p N N T x as Y O C N 2 2 O Z 0 i Y O > L N 4� N y N L X C d Y Y C N X Y N co Z 2 Z p N m :r co Qt O N (") <f LO �O APPENDIX 4 LENGTH OF THE OPEN -WATER REACH BELOW A DAM OR RESERVOIR: Report to the State of Alaska, Department of Commerce Economic Development by J. P. Gosink Geophysical Institute University of Alaska Fairbanks, Alaska 99701 May 7, 1984 ABSTRACT The prediction of the open water length downstream from a dam is an essential safety concern for hydroelectric development in Alaska. This information provides the position of the ice front and determines the stability of that ice front during changes in atmospheric conditions and/or changes in discharge from the dam. During very cold weather (lair < -25°C) the open water reach will be the site of severe ice fog, causing icing on structures, visibility problems, and adversely affecting nearby residents. In addition, the open water reach may eliminate a traditional winter crossing route for man and animals. In this report we examine different approaches for the prediction of the open water length; they are compared for simplicity, for generality and for accuracy. Formulae for direct application of certain of the models are given in tabular and/or graphical format. Several simple -to - use analytic formulae are given for steady-state and transient boundary conditions. The impact of various complications, including lateral temperature gradients, effects of side streams, water clarity and braided channels, which characterize realistic conditions in Alaskan rivers but which unfortunately are not included in the simpler formulae, are discussed and methods are suggested for the quantitative analyses of these problems. Finally, a finite difference computer program of the transient river temperature distribution for the single channel, constant discharge case is given. INTRODUCTION Hydroelectric development in Alaska is proceeding at an accelerating pace. The recently completed hydroelectric projects at Solomon Gulch, Green Lake and Tyee Lake will provide 48.5 megawatts to cities in southeastern Alaska. Other projects either under construction or recommended for construction, including the Susitna site, can provide some 1800 megawatts to the state. The creation of new reservoirs or the deepening of existing lakes and reservoirs can drastically alter the thermal regime in the lake basin and in the downstream river. Water released from the dam will be warmer in the winter and colder in the summer than under pre -construction conditions. In Arctic and sub -arctic regions the temperature of the outfall water during winter is a critical parameter controlling the length of the open water reach downstream from the dam and the position of the leading ice edge. There is a great deal of concern regarding the length of the open water reach since the released water vapor will cause icing on nearby structures and equipment, and will produce thick ice fog during periods of extremely cold temperatures. In addition traditional winter crossing routes for man and animals would be eliminated by the open water reach. Several different methods have been used to determine the length of the open water reach. In general these methods could be classified as either statistical or semi -empirical. The first class uses data acquired for many years at a particular site to establish a curve or set of limits for the length of the open water reach as a function of meteorological and discharge parameters. Two examples of this procedure are the analyses of Goryunov and Perzhinskiy (1967) and of Gotlib and Gorina (1974). Only e: one of the examples of the statistical method shows actual comparisons of the predictions with measured open water lengths. The statistical method is useful only at sites where there exists a long data base for analysis. In addition, the predictions are no longer valid when the hydrology of the reservoir basin and river system are appreciably altered, as for example, by deepening of the reservoir. The other class of semi -empirical methods for finding the open water length is analytical in the sense that some attempt is made to model the basic physics of the problem. These models vary in the assumptions made, but in general, they utilize a semi -empirical heat balance for the open water reach. A major shortcoming of all the models considered in this report is that none takes the dynamics of the ice cover into consideration; that is, all of the models are primarily thermodynamic. This approach is suitable as long as the ambient conditions (discharge and meteorology) are relatively stable, so that changes in the ice conditions occur rela- tively slowly. These models are not applicable for example, during a sudden thaw or with a sudden drastic change in discharge. For stable winter conditions, the analytic models yield reasonably accurate predictions of open water length. Both steady state and transient models are available, and the steady state assumption allows a particularly simple closed form solution to be written for the open water length. In the present report we introduce a closed form solution to the transient problem which is exact whenever the air temperature and/or short wave radiation can be expressed as a sum of sinusoidal terms of arbitrary frequencies - a fairly common case. All closed form solutions are based on the assumption of uniform river hydrology, i.e., constant width, depth, velocity, discharge and specifically, no braided channels or 3 stream inflow. If variations in these parameters are to be included, a finite difference or finite element solution of the governing equations is necessary. An example of this type of finite difference model is given by Ashton (1979). His model allows arbitrary variations in air temperature, and changes in river width and depth and may be modified to improve the surface heat transfer expression or to include the thermal effects of inflowing streams. The purpose of this report is to summarize and assess models for the prediction of the open water length downstream from a dam in arctic and sub -arctic conditions. We include two statistical models to demonstrate their use and the required data. The primary emphasis is on analytical models which are of general applicability. We explain the derivation of the governing equations and differences in the surface heat transfer expression. Using comparisons with data from sub -arctic rivers, we demonstrate that the Dingman and Assur (1969) version of the "Russian winter equation" for linearized heat exchange provides the best estimate of surface flux. The simple closed form solutions of the heat balance equations are presented for both steady state conditions and for sinusoi- dally varying air temperatures. These closed form solutions are useful estimates of the open water length when there are no side streams entering the river, and little variation in river width and depth. Finally, for more general applications, we present a finite difference model based on the Ashton (1979) model, which may easily be extended to include heat flux from side streams and heat exchange by the Dingman, Weeks and Yen (1967) formulae. Other complications including water clarity and transverse mixing are discussed quantitatively, and recommendations are made for Alaskan applications. 0 Classification of Models to Determine the Lead Area Downstream of a Power Station There are two distinct types of models which predict the open lead area downstream of a power station. The first is the totally statistical technique suggested by Goryunov and Perzhinskiy (1967) and by Gotlib and Gorina (1974). All the remaining models discussed in this paper may be classified as semi -empirical models. The models to be discussed are listed for reference in Table 1. Statistical Techniques 1) Gotlib and Gorina (1974) Gotlib and Gorina present a graph (Figure 1) which represents the length of the open lead downstream from the Bratsk hydroelectric plant under cold winter conditions (air temp.: Dec. and Jan., - 290C). D is outflow discharge in m3/sec, and L is open water length in km. These curves represent the minimum length of the lead; a maximum length of 30-48 km is suggested for warm -winter conditions. Each curve is associated with a specific water outflow temperature at the dam ranging from 1.0°C to 3.0°C with increments of O.M. Apparently winter discharge tempera- tures at the Bratsk hydroelectric site always vary between 1.0°C and 3.0°C. From Figure 1, it is evident that the length of the open lead varies directly with reservoir discharge and with the temperature of the outflow water (Tw). No details are given by Gotlib and Gorina (1974) regarding their analysis; furthermore, no comparison with data is given. Although the length of the open water reach increases with the magnitude of the warm discharge and with the temperature of the discharge, neither increase G O M b N N N N N (gas/Ew) Q a6.�ey�s�a v i Y C q �o L GJ O d L 4 C Y Y E w Vi 3 O O L � s W4-u Y 0] Y N C 'O 00 d L r U Y 4. w C p Q 1 tA Y Y w Y �. ref cc C d N Y 4. L L d C O U G1 O 0 C Y r w+ q De Q7 rU N 09 L D w e Y Y o O U C R 0 w M O C B7 fl) O L O 6++ y J CY O� Y L w L U p r m L Y C14- O CJ QI Y C E ZI o E T 4- S- v 4- L .4 N L 7 \ O1 M LL L•" e s ! � F•. I 1 ; v I df t ® ' V 6NI V I"j V7 V ® o a v o 0 0 0 o v IT" � ro N L- it O L®�ro�cn O w ro fU w � .••. � � y N . e W \ C O i� as ro•� w�e70 4J L p Cd C O � � �• O.B ro O Q •P M o •.• � • N ea .••. L C ® Q2 w 0 as41 ro c cl 10 U O Y�'i t= •. roa.c�.av�E C1 L Q9� O O ro C i.a �i T _ wm 4fl .r 3: w V- ro U ro ro d L N c L ws+aq\ •.� cs a u w e" L L C O E O O 7 y C7 L 1�I9 M O ro O Ce L. 5 C. ro U ,--+ +j L U \ r•� c U1P5 O Ql •r \ C7 L ro M CD CA N E = � CD N a 00 O M CD � G ++r{N ro N m is linear. This implies that extrapolation beyond the range of the curves is impossible, since the systematic variation of the open water length with these parameters is not provided, nor is any indication of open water length for varying winter air temperatures. Gotlib and Gorina (1974) also provide graphical information con- cerning the date of the initiation of ice edge recession (i.e., date of onset of break up) as a function of positive degree days, reservoir discharge, and existing ice thickness (see Figure 2). Parameters for the six curves, discharge and initial ice thickness, are given in the caption. Figure 2 demonstrates that for constant ice thickness, break-up of the ice edge begins 3-4 days earlier when the discharge D = 3500 m3/sec than when D = 3000 m3/sec. Furthermore, for constant discharge, say 0 = 3500 m3/sec, ice edge recession is delayed 4 days for every 0.2 meter of ambient ice thickness above 0.6 meter. 2) Goryunov and Perzhinskiy Goryunov and Perzhinskiy (1967) present an empirical curve (see Figure 3) to represent the relation between the length of the open water lead, L, and the sum of the negative degree days. The formula suggested by Goryunov and Perzhinskiy (1967) is: L (in km) = 5.5 • 106 (E-Tair)-2 (1) This formula is applicable to the Lower Volga downstream of the Volgograd Reservoir. Discharge rates and meteorological conditions are not given; thus direct comparison with analytic techniques cannot be made. The data points represent the open lead length for a particular winter as a function of total negative degrees. It is not clear whether this length is a minimum or a winter average. It would be interesting to see whether discharge rates from the Volgograd Reservoir varied during 0 r f0c �- 1960-61 •- 1961 -B2 • - 1962 - 63 0- 1963-64 a- 1964 - 6S 0 e A � o � a a ZU° 300 400 SOO km v a „ � �• G O u = n � a u z r y z a Figure 3. (from Goryunov and Perzhinskiy, 1967) Open water length downstream of the Volgograd Reservoir in km vs. sum of negative degree days. Note that ice covers Yenotayevka Pond with less than 100 decree days. the years of observation, since the analytical models discussed in the following sections all predict a linear increase of open water length with discharge. The data from the Lower Volga suggest that the least variabi- lity in open water length occurs during the warmer periods (smaller degree days). Furthermore, local hydrologic effects would appear to play a major role in establishing ice coverage at Yenotayevka pond; although never stated explicitly, it is reasonable to assume that the pond is a wide river reach with slow water velocities. Finally, it should be noted that the general appearance of the L vs. -lair curve found by statistical methods by the Russian investigators is similar to the theoretical curve predicted by the semi -empirical models. The latter models suggest a relationship of the form L a In [1+Tw/-Tairl for steady-state conditions, and this logarithmic function approaches L - 1/-Tair when Tw << -Tair- Statistical models can provide useful guidelines at existing sites where a good data base already exists. They have no predictive value at the site for any major alterations in the reservoir -river system or for weather extremes. They are not useful as predictive tools for the planning of new projects or expansions. These statistical models yield the following qualitative information on open water length: open water length decreases with negative degree days and with decreasing dam water outflow temperature, and increases with reservoir discharge. With respect to the timing of ice cover break-up, the statistical models suggest accelerated break-up with increased discharge and thinner initial ice thickness. 10 Analytical Techniques Another approach to finding the open water length involves analysis of the basic physics or thermodynamics of the flow. Consider the thermal balance on a slab of fluid: UT + a(UT)/ax ox] Rate of change of = Heat - Heat + heat flux + Other heat in the slab convected in convected out through top heat sources (1) pcp a(hboxT)/at = pcphbUT-pcp(hbUT+a(hbUT)/ax ax) + Q box + E S or (2) pcp [a(hbT)/at + a(UhbT)/ax] = Qb + ES' where p is water density, cp is specific heat, h is river depth, b is river width, T is water temperature, t is time, U is average streamwise velocity, x is streamwise distance, Q is net surface heat exchange [W/m2], and E S' represents the sum of other heat sources including side stream inflow and longitudinal heat diffusion. Initially we will consider only rivers with constant discharge (Uhb = D = constant), constant width, depth and velocity, and zero stream inflow. Then the governing equation simplifies to: (3) pcph[aT/at + U aT/ax] = Q 11 All of the analytical models use a simplified form of equation (1), usually equation (3). Q, the surface heat transfer expression, is determined from semi -empirical models for radiative, turbulent, latent and bottom heat exchange; the formulae for Q vary substantially between the different models, and the complexity of the formulae for Q determines whether or not a closed form solution is available for T(x,t). The expressions for Q take three forms: Q is a function of atmospheric parameters only, Q is linearly proportional to the air and water temperature a difference, and Q is a complex function of water temperature and atmospheric parameters. Details of the second and third types of expression will be given in subsequent sections. In the following section we discuss models for open water length based on all three types of expressions for Q. We have listed the models in the order of increasing mathematical complexity, i.e., increasing complexity of the functional form of Q. In a subsequent section predictions of the models and field data will be compared to assess their realiability; finally guidelines will be offered for the selection of an appropriate model for a given application. Asvall (1972) Asvall (1972) greatly simplifies equation (1) by assuming steady state conditions, constant discharge, river depth and width, and a surface heat transfer expression for Q which depends only on atmospheric conditions. Asvall suggests using the net surface heat loss expression for Q from Devik (1964); this is given in equation 35 of this report and will be discussed subsequently. Since Q is assumed to be a constant (a known function of air temperature and wind velocity), equation (1) may be integrated directly to become, 12 (4) PcpUhTo/IQI = L where L is the open water length and To is the outflow temperature at the dam. Another way of determining this simple formula for L is by equating the net heat into the river at the dam (= PcpUhbTo) and the net heat lost over the open water area (=QLb). However equation 4 fails to take into consideration the fact that the surface exchange Q is a function of water temperature, time and river location. In order to incorporate in a simple way the variation of Q with water temperature, various linearized expressions have been determined for the surface heat exchange. A particularly useful linearization formula expresses Q as a linear function of water temperature, (5) Q=A+BT When an expression of the form of equation (5) is assumed for Q, simple steady state solutions of equation (1) exist, and yield more reliable estimates of the open water length. Consequently, there is no real advantage in using constant values of Q, and a real physical advantage in including water temperature dependency in the formulation of Q. Dingman and Assur (1969) Dingman and Assur (1969) introduce a simple steady state analysis for open water length. The major simplification comes from the linearization of the surface heat transfer expression Q, given previously in equation (5) and written in a more general form below as, (6) Q = -Qo - K(T-Tair) 13 The linear coefficients, Qo and K, were obtained by analysis of the empirical expressions of Dingman at al. (1967) for net long wave sensible and latent heat flux; net measured short wave radiation should be added to the expression. The coefficients are given as functions of wind velocity and cloud cover in Table 2. With this linearization, the steady state heat balance from equation (3) becomes (7) pcphU dT/dx = -Qo -K(T-Tair) and this equation has the closed form solution, (8) T = To -CQo/K + To - Tair][1-exp(-Kx/pcpUh)] where To = T(x=o) is the average well -mixed temperature at the outfall. The temperature of the water decreases exponentially with distance from the outfall, and approaches a theoretical equilibrium temperature (Te = Tair - Qo/K) at x = In actuality the temperature decreases to DOC at the leading edge of the ice; beyond this distance the expression for surface heat exchange is no longer valid, the water temperature remains DOC, and heat loss through the ice cover implies ice growth. The position of the zero isotherm, L, can be found by setting T = 0°C in equation (8): (9) L = (pcpUh/K) In [1 + KTo/(Qo-K Tair)] Note that the coefficients for atmospheric heat transfer, Qo and K, are simple functions of wind velocity and air temperature; clearly equation (9) is a very easy -to -use formula for open water length. However, the assumptions required for this derivation should be kept in mind. These include: 1) uniform and constant river discharge, width and depth, 2) constant air tempratures and wind velocity, 3) no inflowing streams, 14 4) no heat flux from or to the river bottom, 5) applicability of the linearized surface heat transfer expressions. Where these assunptions are violated, an appropriate strategy might be the use of equation (9) as a first estimate of the open water depth, with subsequent analysis of the effects of other parameters. Quantitative discussion of some of these parameters follows in a subsequent section. Paily, Magagno and Kennedy (1974) Paily et al. (1974) solve the following version of equation (2), (10) aT/at + U aT/ax = Q/pc ph + E 32T/3x2 with a linearized heat exchange expression for Q similar to the expression used by Dingman and Assur (1969), but involving different values of the linearization coefficients K and Qo. In Paily et al. (1974) the coeffi- cients are given in tabular form rather than as functions of wind speed and cloud cover; the coefficients are presented in Table 3 of the present report. The additional term, Ea2T/ax2, represents streamwise diffusion of heat. Obviously diffusion is a much less effective mechanism for heat transport in a river than is convection. Nevertheless it is included in this model for completeness and to demonstrate the relative effect of longitudinal diffusion of heat. For steady state cases, Paily et al. (1974) give a closed form solution of equation (10); (11) T = To - [(Qo/K) + To - Tair]C1-exp{(Ux/2E)(1- 1+4KE/pcphU,2)' 11 This solution was devised earlier by Daily and Harleman (1966). It is important to note that in the limit as E approaches zero, the argument of the exponential term goes to-Kx/pcpUh, exactly as predicted by the Dingman-Assur (1969) model (see equation 8); this can be seen either by 15 Table 3. Values of Qo and K from Paily et al. (1974) i I I I Heat ex - Base I change co - Wind I I heat I efficient, K I velocity I 1 exchange I in Watts I in miles I I rate, Qo I per square Air tem- I per hour I Relative I in Watts I meter perature, I (meters I humidity, I per square per degree in degrees I per I as a per- I meter I Celsius Celsius I second) I centage I (1) i (2) i (3) i (4) i (5) -1.0 11.0 70.0 16.25 31.40 (4.95) -3.0 11.0 70.0 65.35 32.50 (4.95) -5.0 1170.0 70.0 114.67 33.58 (4.95) -10.0 11.0 70.0 239.39 36.22 (4.95) -15.0 11.0 70.0 366.96 38.77 (4.95) -18.0 11.0 70.0 445.27 40.28 (4.95) -5.0 0.0 70.0 23.04 16.67 (0.0) -5.0 3.7 70.0 53.59 23.30 (1.65) -5.0 7.4 70.0 84.13 27.94 (3.30) -5.0 11.0 70.0 114.67 33.58 (4.95) -5.0 14.7 70.0 145.22 39.21 (6.60) -5.0 18.4 70.0 175.76 44.85 (8.25) -5.0 11.0 10.0 171.79 34.25 (4.95) -5.0 11.0 30.0 152.75 34.02 (4.95) -5.0 11.0 50.0 133.71 69.80 (4.95) -5.0 11.0 70.0 114.67 33.58 (4.95) -5.0 11.0 90.0 95.64 33.35 (4.95) -5.0 11.0 100.0 86.12 33.24 (4.95) -for a ues valid range of water temperature between 0 C and C; va ues o other meteorological variables are: barometric pressure = 996.0 mb; cloud height = 3,275 ft (1,000 m); cloud cover = .6; and visibility = 1.87 miles (3 km). applying L'Hopital's rule or by expanding the square root with the binomial expansion. The latter procedure yields the following series for the argument, (12) (Ux/2E) 1-2(KE/pcphU2) + 2 (KE/pcphU2)2 - 4(KE/pcphU2)3 + ... } This series shows that the diffusion term lessens the longitudinal temperature decrease, producing a slightly longer open water length. The tempering effect of the diffusion term can also be seen directly from ` equation (10) when it is noticed that the second derivative term is Positive definite in these problems. A closed form expression for the open water length can be written as follows, (13) L = (pc pUh/K)(1/2 + 1/4 + KE/pcphU12)ln[l + KTo/(Qo-K Tair)] Here the effect of the diffusion term on the open water length is immediately apparent. Clearly when E << 4pcphU2/K there is very little increase in open water length. Numerical estimates of this increment for typical Alaskan conditions will be given in a subsequent section. Note that if E is small, then differences between values of L calculated by the Paily et al. (1974) formulae and the Dingman and Assur (1969) formulae will depend primarily upon the linearization coefficients, K and Qo, in the respective formulae. Paily et al. (1974) also provide a closed form solution for the transient case of equation (10) for linearized surface heat exchange and particular initial boundary conditions. However the specific initial and boundary conditions assumed by these authors are not appropriate for the temperature regime for water released from a dam. Paily et al. (1974) 17 are interested in the temperature regime in a flowing river with a heat source at x = o in which the entire river including upstream (x < o) is subject to atmospheric heat transfer. That is to say, the water arriving at x = o from upstream is changing temperature due to atmospheric forcing. The application of the Paily et al. (1974) model is to temperature prediction in a river with a thermal effluent injected at x = o. Therefore they assume that the boundary temprature T(x = o, t) is not constant, but instead, equals the sum of the inflow temperature To plus the transient river response to uniform atmospheric heat transfer. This boundary condition is not appropriate for the water released from a dam. Water released from the dam is at a constant temperature since this water comes from depth below the ice cover, and reservoir water under the ice cover has very little if any diurnal temperature variation. During breakup or during intense wind mixing, or when alternative outlets from the dam are used, the released water temperature will vary, but the released water temperature cannot be predicted from a simple river temperature model. It is essential that a reliable reservoir temperature model be used to define the outflow temperature. In the present analysis we consider the outflow temperature as given either through measurements or by prediction from a reservoir model. Other analytic solutions An interesting and useful analytic solution can be found for the problem of the transient response of the river to periodically varying meteorological conditions. The meteorological condition may represent diurnal variation in air temperature and/or short wave radiation, or alternately, seasonal climatic variation. The formal statement of the problem consists of the governing equation (equation 3) with the atmospheric heat transfer expression as follows, I (14) Q = - Qo - K(T-Tair) + Q sin wt with the initial conditions, (15) T(x,o) = To - [(Qo/K) + To - Tair][1-exp(-Kx/pcpUh)] Note that this initial condition has the expected behavior at x = o, i.e., T(x = o, t) = To, the constant outflow temperature. Furthermore, the initial condition is actually the steady state temperature for the case when Q = 0. The solution then defines the transient river response to sinusoidal atmospheric forcing. The analytic solution to this problem is, (16) T = To - [(Qo/K) + To - Tair]11-exp(-kx/Uh)1 + AT {sin(wt-a) - exp(-kx/Uh) sin (wt-wx/U-g)} where k = K/pcp AT = Q/pcp k2+w2h2 and a = sin-1 [wh/ k2+w2h2] The form of the solution highlights the role played by the periodic air boundary conditions. If Q =_ 0, or a constant air temperature is assumed, the solution reduces to the steady-state case. When Q + 0, the periodic nature of the temperature distribution in the river becomes evident. The river temperature lags the air temperature by the phase angle a. This Phase lag is directly proportional to river depth and inversely proportional to the surface heat loss coefficient, matching the intuitive expectation for river temperature adjustment to air temperature variation. That is, shallow rivers (h + 0) cool faster than deeper rivers with the same discharge, and rapid heat transfer (k >> 0), which occurs for example with high winds, is characterized by rapid temperature adjustment. An 19 estimate of the typical diurnal adjustment time lag for winter conditions may be found by assuming reasonable values for K, w and h: w = 27r/(24.3600) s-1, h = 3 m, and k = 7 • 10-6 m/sec. This value of k corresponds to a coefficient of Tair in the Dingman-Assur formula (Table 2) equal to 30 W/m2. These parameters suggest a daily time lag between air and water of about 5.8 hours during winter conditions. A similar estimate can be made for a seasonal time lag when an annual period is assumed for the air tem- perature; this estimate suggests a time lag of about 5 days. A still more general transient solution may be found for the case where the atmospheric heat transfer can be represented by a sum of periodically varying terms of arbitrary frequency and magnitude. This boundary condition may represent the combination of diurnal and seasonal variation in air temperature, and in short wave radiation or other parameters, or it may represent a complex transient surface heat flux determined from measured values by harmonic analysis. For this general case the heat transfer expression is, N (17) Q = -Qo - K (T-Tair) + E Qi sin (wit+oi) i=1 and the initial conditions are given by equation (15). The analytic solution is, (18) T = To - C(Qo/K) + To - Tair][1 - exp(-kx/Uh)] N + E oTi{sin(wit+oi-si)-exp(-kx/Uh)sin(wit+oi-wix/U-si)} i=1 where oTi = Qi/Pco k2+wih2 and Bi = sin-1 Cwih/42+w?h2] 20 Each phase si can be calculated independently, and each phase lag is directly related to the period of the respective heat flux fluctuation. The amplitude of the periodic temperature waves in the river, oTi is inversely proportional to the forcing frequency; i.e., short period fluctu- ations in air temperature are hardly felt in the river and longer period fluctuations are strongly impressed upon river temperature. In all cases, the amplitude of periodic temperature waves in the river is inversely related to river depth, and if river depth is very small, that amplitude approaches the amplitude of air temperature variation (Qi/K). Finally, a slightly more general transient solution may be found for the case where the atmospheric heat transfer varies in a known way as a Polynomial function of river distance. This boundary condition may represent a spatially varying air temprature because of lapse rate, weather pattern, or systematic change in radiative heating. For this general case the heat transfer expression is, N — M (19) Q = -Qo - K(T - Tair) + E Qi sin(wit+(Di) + E gixi i=1 i=1 where qi represent the known longitudinal variation, and the initial conditions are given by equation (15). The analytic solution is given by equation (18) plus a linear summation from the longitudinal variation: (20) T = To - C(Qo/K) + To - Tair][1-exp (-kx/Uh)] N + E oTi{sin(wit+Oi-Si)-exp(-kx/Uh)sin(wit+oi-wix/U-si)} i=1 M + E gixi+l/(i+1)PcpUh i=1 21 This solution is the most general closed form expression for the temperature regime downstream from a dam when the surface heat transfer has been linearized. Allowable functional forms for the surface heat transfer (equation 19) can be quite general including differing periodicity of air temperature and radiation, as well as combinations of diurnal, seasonal and episodic events and arbitrary persistent longitudinal variation. Hone of the transient analytic solutions for temperature (equations 16, 18 or 20) can be directly inverted to determine open water length since the equations are transcendental. However the temperature regime can be easily calculated as a function of x and t, and, for a particular time, the open water length determined directly. It is important to remember the limitations of all the analytic models. First, none of the analytic thermal models include latent heat exchange with an ice cover and are therefore only useful for river temperatures above or equal to 0°C. They can be directly applied only in uniform river stretches, i.e., with no variation in river width, depth and velocity and no inflowing streams. The allowable heat transfer functions, although reasonably general, are based on linearized analysis of higher order surface heat transfer expressions, and the appropriateness of the linearizations must be considered. In the following section we shall consider semi -empirical formulae for surface heat expressions, and discuss some assumptions involved in the linearization of these formulae. Measurements of open water length in typical Alaskan conditions will be compared with predictions from the different linearization expressions. 22 Heat Transfer Expressions Dingman, Weeks and Yen (1967) provide a very extensive analysis of the mechanisms of heat transfer to a flowing stream. These authors consider the following eight heat transfer terms: (21) Q=QR+Qg+QE+QH+QS+QG+QGW+QF where QR short wave radiative flux QB net long wave exchange with the atmosphere QE evaporative heat exchange QH sensible or turbulent heat flux QS heat lost by influx of snow QG heat added by geothermal transfer QGW heat added by ground water QF heat added by friction from stream bottom The expressions for each of these terms are given in Table 2. Dingman et al. (1967) were particularly interested in the selection of appropriate expressions for QE and QH in arctic and sub -arctic conditions. They compared the formulae of Kohler (1954) and of Rimsha and Donchenko (1957) to cold region data and determined that the "Russian winter equation" as given by Rimsha and Donchenko (1957) was the more accurate of the two formulae. We have included both the Kohler (1954) and the Rimsha and Donchenko (1957) formulae for QE and QH in Table 2 for comparison. More recent formulae for water -atmospheric heat transfer have been given by the Tennessee Valley Authority (1972), Hicks (1972), Pond et al. (1974), and Holmgren and Weller (1968); however, the first three of these were devised primarily for temperate regions, and all four were devised for deep water. McFadden (1974) presented a comprehensive comparison of heat 23 transfer formula with measurements for arctic conditions; the reader is referred to that report for details of the comparisons. In this report we shall not attempt to compare in detail the formulae for heat transfer mechanisms given by each author. Instead we shall make recommendations for both the full empirical formulae and for the linearized versions of these formulae based on our calculations and those of Dingman et al. (1967) and of McFadden (1974). In all these discussions the units of heat flux, Q, are W/m2. QR: Short wave radiation Short wave radiation is always positive and represents a relatively small component of the heat budget of Alaskan rivers in winter. McFadden (1974) cites several references which report the daily flux of short wave radiation near 65° latitude to be less than 5 W/m2 in December. Wendler (1980) gives the average measured short wave radiative flux as less than 5 W/m2 during November, December and January. This contrasts with lower latitudes where the short wave radiation is often the dominant mode of heat transfer to a water surface (e.g., see Fischer et al. 1979). Because of the reliability and simplicity of short wave radiometer systems, it is recommended that short wave radiation be measured directly at the site, and the measured values used in the calculations for open water length. In the linearization formulae, QR can be added directly to the heat flux terms. If short wave radiation measurements are not available, then the following estimation procedure modified from Dingman et al. (1967) is recommended: (21) QR = 0.892 QRI + 1.397 • 10-4 QR12 [W/m2] and 24 QRI = QCL (0.96 - 0.61C) [W/m2] where C is cloud cover in tenths (e.g., complete cloud cover implies C = 1.0) alnd QCL is incoming short wave radiation for a cloudless sky. QCL may be found for various latitudes as functions of season in tablular and graphical form (see TVA, 1972 and Bolsenga, 1964). A distinctive feature of the short wave radiation fl ux is the fact that it is not completely absorbed at the water surface; it penetrates to some depth depending upon the water clarity and turbidity. The short wave flux available at a depth y is usually assumed to follow Bouger's Law for absorption: (22) QR (Y) = QR (Y=o) exp (-nY) where n is an extinction coefficient ranging from about 0.2 m-1 for very clear water,to 4.0 m-1 for turbid water. This implies that in a very clear shallow stream with depth equal to 1 meter, only 20% of the short wave radiation is absorbed by the water column, and the remaining 80% penetrates into the river bottom. At night some of this stored heat flux is released into the water column, implying an increase in geothermal heat flux QG which lags the short wave flux. In sediment laden streams Ti may be even larger than 4.0 m-1, and therefore virtually all short wave radiation is absorbed in the topmost meter of the water column. In order to reliably model the bottom flux it would be necessary to couple the river temperature model to a ground thermal model. However, due to the fact that QR in late fall is only a minor component in the thermal budgets for high latitude rivers, it is uaualiy possible to ignore heat absorption in the river bottom. In any case the main effect of bottom heat absorption on river temperatures would be a lag in the diurnal temperature maximum of the river 25 or a slower decrease in river temperature in the evening. Quantitatively QR will represent less than 5% of the overall river heat budget from late fall through early spring; therefore the lagged release of heat from bottom sediments may equal 4% of the river heat budget in clear streams. In this report we will not propose a mathematical model which couples predictions of the ground thermal regime to predictions of water temperatures. For rivers deeper than 2 meters and in sediment laden streams, we recommend assuming that QR is entirely absorbed by the river, unless it is critical at the particular site to determine the diurnal variation in water tempera- ture. For clear shallower streams we recommend that an experimental study be undertaken to determine the diurnal lag in river temperature due to gradual release of stored radiative heat in the river bottom. QB: Net long wave radiation exchange with the atmosphere Net long wave exchange with the atmosphere consists of the outgoing long wave radiation emitted from the water surface QW plus the net incoming radiation from the atmosphere Qp: (23) QB = -QW + Qq The net long wave exchange may be measured directly at the site. If these measurements are not available, then the long wave exchange may be estimated by semi -empirical formulae relating QB to water and air temperatures. The radiation from the water surface is modeled by the Stefan formula, (24) QW = eWa(T + 273)4 where eW is the emissivity of water (=0.97), a is the Stefan-Boltzman constant (5.67 • 10-8 W/m2K4) and T is the surface water temperature in °C. This formula is widely accepted in the literature and is recommended here. 26 There has been some speculation that at the time of ice formation a thin supercooled layer of water may exist on the river surface. While this assumption may be valid in quiescent ponds, it has been shown to be unfounded in turbulent rivers (Osterkamp et al., 1983). For rivers with mean velocity greater than about 0.6 m/sec the surface water temperature may be assumed to be equal to the mean river temperature. Usually the atmospheric radiation can be modeled by a Stefan formula (25) QA = f(e,C,H,a) v (Tair + 273)4 where f is a function of air vapor pressure (e), cloud cover (C), cloud height (H), absorptivity of the water surface (a) and Tair is the air temperature in °C at a specific height, usually 10 meters. McFadden (1974) has discussed various expressions for f in some detail, and proposed a complex formula especially for cooling ponds which includes an additional dependency on the cooling pond shape factor. McFadden (1974) also compared long wave radiation data at a site with ice fog with the predictions of long wave radiation determined by the formulas of Brunt (1944), Angstrom (1920), Elsasser (1942) and Andersen (1952), and used a correlation technique to modify these formula and thereby improve the agreement with the data. We recommend McFadden's (1974) modified version of the Andersen (1952) formula. The Anderson equation (1952), both in the original format and in the modified version, exhibited the minimum standard error of all those investigated. The Andersen equation (1952) was also adopted in the Dingman, Weeks and Yen (1967) river thermal model. The modified version of the Andersen equation is: (26) QA = E.814 + .11C exp (-.19H) + ea (.0054 - .000594 C exp (-.197H))] a (Tair + 273)4 27 where H is cloud height in km and ea is vapor pressure of the air in mb. QE: Evaporative heat flux QH: Turbulent heat flux Semi -empirical formulae for turbulent heat flux are usually written in the form, (27) QH = (A + Bw)(T-lair) where w is wind speed, and A and B are empirically determined parameters. There is an extensive core of literature related to the determination of A and B (e.g., Friehe and Schmitts, 1976; Kohler, 1954; Rimsha and Donchenko, 1957; TVA, 1972; Hicks, 1972; Kays, 1966). The form of the equation models the intensification of convective or turbulent heat transfer by strong winds and increased temperature difference between the air and water. In addition, the parameter A assures upward heat transfer from a water surface which is warmer than the air even when the wind velocity is small. This situation frequently occurs in interior Alaska where air temperatures 300 below water temperatures may exist with no wind. Under these conditions the air is buoyantly unstable, and strong vertical motion in the form of thermal plumes or buoyant convective cells may develop, facilitating surface heat transfer. Evaporative heat loss QE occurs when there isa net upward transport of vapor from the water surface; the heat loss is the product of the specific heat of the vapor and the evaporation rate. There is extensive literature on evaporative heat loss (e.g., Hicks, 1972 and 1975; TVA, 1972; Friehe and Schmitt, 1976; Anderson, 1954; Pasquill, 1949; Rimsha and Donchenko, 1957; Devik, 1964). It is usually assumed to be linearly proportional to the air - water specific humidity difference and is modeled by equations of the form, (28) QE = (C + Dw) (e - eair) a3 where w is wind speed, C and D are empirically determined parameters, a is the saturated vapor pressure of air at the temperature of the water, and eair is the vapor pressure of the air at a specific height, usually 10 meters. It should be noticed that the transfer of water vapor or any gas across the water surface is a complex problem and the subject of intense recent research (e.g., see Brutsaert and Jirka, 1984). Dingman et al. (1967) and McFadden (1974) reviewed several models for QH and QE, and compared the predictions of these models with data from arctic conditions. Both concluded that the Rimsha-Donchenko (1957) formulae for QH and QE more accurately predicted turbulent and evaporative heat exchange in arctic conditions than did other models under consideration. The Rimsha-Donchenko formulae are given here: (29) QH = [3.87 + 0.17 (T - Tair) + 1.89 w ] (T - Tair) (30) QE = [6.04 + 0.264 (T - Tair) + 2.94 w ] (e - Bair) where QH and QE are in W/m2, w is wind speed in m/sec, T is water temperature, e is saturated vapor pressure at T, Tair is air temperature at 2 meters, and eair is vapor pressure at 2 meters. QS, QG, QGW and QF For the four types of heat transfer, QS, QE, QGW and QF we follow the recommendations of Dingman et al. (1967). Latent heat exchange from snow QS is proportional to the snow accumulation rate A: (31) Qs = cA [G + Ci (T - Tair)] where A is given in g/cm day, a is the latent heat of fusion of ice in cal/g, Ci is the heat capacity of ice in cal/g °C, and c is a dimension 29 conversion constant c = 0.484 [W/m2 : cal/cm2 day]. If snow accumulation rate is not available, then A may be estimated as a function of visibility by an expression of the form (Mellor, 1964): (32) A = 7.85 v"2.375 where v is visibility in km. For consistency the river discharge should be increased by A times the river width, although the net change in discharge would be very small. QG is the geothermal heat flux below the river plus heat released from bottom sediments and must be determined from local data. The geothermal flux is expected to be small except possibly in areas of high geothermal flux (Osterkamp, Kawasaki and Gosink, 1983). As discussed earlier some of the daily short wave radiation QR may penetrate through the river and be absorbed into the river bottom. This stored heat may then be released later in the day, thus delaying the diurnal river temperature decrease. Accurate knowledge of this effect can only be established by analysis which couples temperature distribution in the river with temperature distribution in the bottom sediments. The effect will not be significant (< 4% of total heat flux) for rivers deeper than 1 meter with extinction coefficient greater than about 0.2 m-1. If it is essential to determine the diurnal temperature regime in a very shallow and clear stream, then a more complex coupled analysis of river and sediment temperature is necessary. As a general rule, where there are no indications of high geothermal heat flux, where the river is deeper than about 2 m, and where the short wave extinction coefficient is greater than 0.2 m-1, the total geothermal flux QG may be considered negligible. IN QGW is the heat added by flow of ground water and smaller streams into the river. In order to model this heat flux, information is needed on both the ground water recharge or stream discharge and the temperature of the inflowing water. Note that QGW affects both the right hand and left hand sides of equation (3), by changing the heat input and the river discharge respectively. If stream inflow and temperature measurements are known, these may be incorporated into the model by relatively small changes in the finite difference form of equation (3). QF is the heat added to the river due to friction of the water flowing over the river bottom. It is generally assumed that the decrease in potential energy in the river as it flows downhill is compensated for by the frictional drag at the bottom; subsequently, the drag creates turbulent eddies which, through the turbulent energy cascade, ultimately cause viscous heating. The major problem with this assumption is the neglect of the wall (river bottom) temperature, since if the river bottom is colder than the bulk river temperature, frictional heating will be directed downward into the sediment (Schlichting, 1968). Therefore, the model for frictional heating suggested here and in Dingman et al. (1967) or Starosolszky (1970) should be considered an upper limit to heat flux by frictional heating of the river. The relation between bottom shear stress and the change in potential energy of a volume of water is given by standard hydraulic theory (Henderson, 1966). The shear stress at the river bottom is, (33) Tw = PwghS [kg/m-secZ] where pw is water density in kg/m3, g is the gravitational constant in m/secZ, h is river depth in m, and S is the slope of the water surface. Then the heat flux generated by this stress is (Ince and Ashe, 1964), 31 (34) QF = U Tw = PwgUhS For steep rivers both U and S may be high, suggesting that frictional heating may be a significant fraction of the total heat transfer Q. For example, for h = 3m, U = 2m/sec and S = 10-3, QF = 60 W/m2, and for water at 0°C, the long wave radiative flux QW = 305 W/m2. Even if it is assumed that half the frictional heating is directed upward into the water, QF represents at least 10% of the long wave radiation and therefore should be included in the total budget. It should be noted that Dingman et al. (1967) suggest that QF is insignificant while Starosolszky (1970) recommends that QF be included in the heat budget. QF is relatively easy to estimate for a given river reach, and its magnitude may be included in the governing equation (equation 3) as an additive constant, posing no real complication to the solution of the governing equation. We suggest including QF when the river slope is greater than about 10-4 Linearization formulae The long wave radiation from the water surface (equation 24) and the turbulent heat flux (equation 29) depend nonlinearly upon water temperature; due to this fact an analytic solution of equation 3 is generally not available. However, it is possible to solve equation 3 when all heat fluxes are expressed as linear combinations of water temperature and other parameters, as demonstrated by the solutions given in equations 4, 8, 11, 16, 18 and 20. Therefore, several authors have determined linearized forms of several terms in the heat budget, specifically Q8, QH and QE. It is assumed that since the remaining heat flux terms, QR, QS, QG. QGW and QF are not dependent upon water temperature, their cumulative effect is equivalent to an additive constant in Q, i.e., they are simply added to the linearization constant Qo in equation 6: 32 (6) Q = QR + QS + QG + QGW + Qp - Qo - K (T - Tair) Using regression techniques, Dingman and Assur (1969) determined the following expressions for Qo and K: Qo 50.93 + 11.21 w (clear sky) (35) -35.28 + 4.40 w (cloudy sky) K = 16.99 + 2.05 w (clear sky) 17.97 + 2.22 w (cloudy sky) where w is wind speed in m/sec (the height of the anemometer was not given), and the units of Qo are W/m2 and of K, W/m2-°C. These expressions are linearizations of the Dingman et al. (1967) heat flux formulae for long wave radiation, and turbulent and latent heat flux (see equations 24, 26, 29 and 30). Other linearization expressions include formulae derived specifically for a reach of the St. Lawrence River by Pruden et al. (1954): QB + QH + QE = - 88.91 - 7.5 Tair - 20.87 (T - Tair) and the formulae given by Asvall (1972) and adapted from Devik (1964): (36) 136.05 + 2.09 w C = 0.0 Qo = 77.38 + 2.09 w C = 0.5 23.00 + 2.09 w C = 1.0 12.59 + 1.63 w C = 0.0 K' = 9.44 + 2.41 w C = 0.5 10.92 + 2.05 w C = 1.0 where C is cloud cover and K' multiplies - Tair (°C) instead of T - Tair (°C) as in equation (6). 33 Paily et al. (1974) also determined empirical fits to the Dingman et al. (1967) formulae for QB, QH and QE by a least squares polynomial approxi- mation technique. Values of Qo and K determined by Paily et al..(1974) are given in tabular form in Table 3. These coefficients differ from the set given by Dingman and Assur (1969) and there are two reasons for the differ- ences: 1) values of Qo and K from equations (35) are not dependent on air temperature while the Paily et al. (1974) coefficients are, and 2) the Paily et al. (1974) coefficients were selected as best fits over the range of air temperatures -180C < Tair < DOC, while the coefficients in equation (35) were selected as best fits over the range of air temperatures -50°C < Tair < 0°C. This latter effect becomes critical for application to Alaskan rivers. Although the Paily et al. (1974) expressions for Qo and K are reliable within their range of applicability, they deviate from the complex Dingman et al. (1967) formulae when air temperatures are substan- tially below -18°C. Since only discrete values of Qo and K are given in the Paily et al. (1974) report, we have determined the following interpolation formulae which agree with their tabular values with a maximum deviation of 1.7% and an average deviation less than 0.5%. (37) K = 14.795 + 3.45 w - 1.11 • 10-2 q + .540ITairl - 1.12 • 10-31Tair12 Qo=-32.796 + 18.513 w -.952q + (24.290-K)ITairl + 4.016 • 10-21Tair12 + 2.696 • 10-51Tair14 where w is wind speed in m/sec, q is humidity in % (100. is saturated), and Tair is air temperature in *C. Note that K must be calculated first, since it is used in the evaluation of Qo. 34 The linearized heat transfer from equations 35, 36 or 37 represents QS the sum of long wave radiative exchange plus evaporative and turbulent heat flux, as given by the approximate expression: QS = Qw - QA + QH + QE = Qo + K (Tw - Tair) Both left and right hand sides are functions of Tw the water temperature. The agreement between the different linearization formulae and the "exact" formulae may be tested for an appropriate range of river temperature and atmospheric conditions. We shall plot both sides of the expression for the range of values, 1.0 z Tw < 4.0, with the terms Qw, QA, QH and QE calculated from equations 24, 26, 29 and 30 respectively and Qo and K from equations 35, 36 and 37. In Figure 4, we assume zero wind velocity and clear sky or zero relative humidity. The "exact" values of QS (as given by equations 24, 26, 29 and 30) are shown for air temperatures Tair = {0, -10, -20, -30, -40 } by the vertical braces. (Note that in some cases the vertical braces have been shifted slightly left or right for clarity). Since the Devik (1964) formula (equation 36) is not a function of water temperature, only a single value of QS may be plotted at each air temperature. The Dingman and Assur (1969) expressions or equation 35 yield the range of QS denoted on Figure 4 by the solid bar; and the Paily et al. (1974) expressions or equation 37 yield the range of QS denoted by the open bar. There is a clear tendency for equation 37 to diverge from the exact solution, becoming less accurate as the air temperature decreases below -20°C. Equation 36 (from Devik (1964)) also diverges from the exact solution with decreasing air temperature. The Dingman and Assur (1969) expression or equation 35 provides the best overall estimate of the exact solution. 35 1200 1000 LEE 600 400 mil] I "exact" solution Dingman and Assur (1969) p Paily et al. (1974) x Devik (1964) 5 Tair x Figure 4. Comparison of complete heat flux equations with linearized approximations for zero wind velocity. W. In Figure 5, we assume a wind velocity of 10 m/sec and clear sky or zero relative humidity. Again the Dingman and Assur (1969) expression or equation 35 provides the best overall agreement with the exact solution, keeping pace with the intense heat transfer associated with high wind -low air temperature. The Devik (1964) expression consistently underestimates the heat transfer rate, and the Paily et al. (1974) expression diverges from the exact solution beyond about -15°C. The Dingman and Assur (1969) ex- pressions or equation 35 are significantly more accurate than the others at low air temperatures. Comparison with Data: Example 1 Studies of ice -free reaches downstream from a warm discharge seldom contain complete meteorological and hydrological conditions. For example, Carlson et al. (1978) do not report air temperature, wind velocity, cloud cover or radiation data. However, this information is sometimes available from local weather records. The information should be acquired from weather stations as close as possible to the study site to minimize errors in the determination of heat loss and whenever possible, at the study site. Carlson et al. (1978) specify that the data were recorded at the MUS Power Plant in Fairbanks during December 1971. Thus, referring to Fairbanks meteorological reports for this period, it is possible to calculate heat loss with the different linearization models, and then to compare calculated and measured open water areas. In particular we wish to compare the linearization formula of Dingman and Assur (1969) (equation 35), Asvall (1972) (equation 36), and Paily et al. (1972) (equation 37) and the analytic solutions for river temperature as given by equations 4, 9 and 13. For the month of December 1971, discharge rates for the Chena River and the MUS Power Plant are approximately 80o ft3/sec and 25 ft3/sec respectively. 37 1600 1600 1400 1000 L d t- x U- u.. AIR 400 200 ;act" solution igman and Assur (196 ly et al. (1974) ,ik (1964) Figure 5. Comparison of complete heat flux equations with linearized approximations for wind velocity = 10 m/s. If complete mixing near the discharge is assumed, the effective temperature rise becomes: To = 10° (25/800) = 0.31°C when the effluent temperature is 10°C. The MUS Power Plant uses two different types of discharge. The first and more conventional mode of discharge is the subsurface diffuser. When this technique is employed, there is considerable turbulent mixing near the diffuser. Consequently, mixing may be assumed to be complete, and the one- dimensional assumption implicit in the models is appropriate. Measurements of the open water length in the Chena when the subsurface diffuser was in use in December 1971 indicated an ice -free area of 15 acres (Carlson et al., 1978). Measurements of the open water length were also made when the second type of discharge, the surface dispersion field, was in use. In this mode, the effluent enters the stream at the surface through a series of pipes with little turbulent mixing. Hence, the dispersion field operates as a surface spreading scheme. Heat transfer is rapid, since heat loss is pro- portional to the temperature differences between the water and the air. As expected, the surface dispersion scheme produces smaller ice -free area; in December 1971; average areas of 8 acres were measured. The surface dispersion field is characterized by strong vertical and lateral temprature gradients. The existence of steep temperature gradients invalidates the assumptions implicit in the one-dimensional models, indicating that comparisons of prediction schemes with existing data are appropriate for only the subsurface diffuser. Fairbanks weather data for the month of December 1971 was compiled by the Environmental Data Service, National Oceanic and Atmospheric Administration. 39 Average temperature for the month was -21°C; wind speed, 3.7 mph = 1.65 m/sec; cloud cover, 0.7. Air temperature was about normal for December, and cloud cover, heavier than normal. We have used the above values of the mean air temperature, wind speed, and cloud cover in each of the linearization formulae (equations 35, 36 and 37) and have determined the parameters Qo and K listed in Table 4. For equations 35 and 36, interpolation between cloud covers of 5 and 10 was required. Equation 37 contains no functional dependence on cloud cover but is dependent upon humidity, which, for this test case was assumed to be 10%. There is a surprising lack of agreement of the calculated values of Qo and K between the different models, particularly between the Dingman and Assur (1969) and the Paily et al. (1974) formulae which both represent linearizations of the Dingman et al. (1967) formulae. However it is encouraging to note that the open water areas predicted by these two expressions are in good agreement and bracket the measured open water length of 15 acres. The Paily formula is somewhat sensitive to the selection of humidity, and when a humidity of 90% is assumed, the predicted open water area is 16.6 acres or identical with the Dingman and Assur (1969) prediciton. We have determined the open water area for the Paily (1974) model using E = 4.51 m2/sec which is the value recommended by Paily et al. (1974) and E = 0 to test the sensitivity of the longitudinal diffusion term; as expected the difference is negligible. It should be noted that 4.51 m2/sec is about twice the value of the longitudinal dispersion coefficient calculated by using the Fischer et al. (1978) expression for dispersion coefficient. As previously suggested, longitudinal diffusion of heat becomes important only for slow rivers, in particular for conditions in which the ratio KE/pcpA is about 0.1 or greater (see equation 13). For this example when E = 4.51, the ratio is 2 • 10-5. Whenever the ratio M In r N -ct V to to N N L U .� Q N tr Q N E .1 M p tp rl O CO tp In tci U N N Ly Ln N i I . p E ct W U 0 N n S n N N N 3 M N M M Y mr co L 3 O N M M .y to Ln MM) O• .-1 .M-1 t I _ Lo + F + f- + F- N C N lD Ln PN nN \ O N d° N N N N N 3 N O NM NM C M ICf 47 C f Y O c6 l0 R N- M M M M A Q C O .- 4C m A II7 W r N n C O\ r 01 C tp m to v C v y y r av as 4% m CM C > C N •r 01 , Q1 O N d fa .--I 'p ,..� 41 KE/PcpA is less than 0.1, we recommend the simpler Dingman and Assur (1969) formula for open water length (equation 9) over the Paily et al. (1974) formula (equation 13). Comparison with Data: Example 2 Data from W.A.C. Bennett Dam on the Peace River in British Columbia can also be used to compare the accuracy of the various models. Measurements of open water length downstream from the dam are available for the winters of 73/74, 74/75, 75/76, 76/77 and 77/78 (British Columbia Hydro and Power Authority, personal communication), with the length varying between 60 and 203 miles during these years. However, it is difficult to apply the theoretical models for open water length directly to the Peace River data for several reasons related to the assumptions implicit in the models: 1) the models (equations 4, 9 and 13) are all steady state cases, implying both steady discharge and meteological parameters; 2) the Peace River meanders in the region of interest and a typical river width is difficult to determine; 3) the closest meteorological data come from Fort St. John about 15 miles downstream from the dam, and meteorological data from this location often disagree substantially with data from the next downstream source, Peace River some 60 miles from the dam. Nevertheless it is useful to determine rough estimates of the open water length for the five winters by using "mean" meteorological and hydrau- lic parameters at the site. Discharge and average outflow temperature are known (British Columbia Hydro and Power Authority, personal communication). We assume a constant river width of 200 meters. A mean air temperature is probably the most subjective choice since it is not clear whether the period of averaging should include the entire winter or a specific period proceeding the time of the minimum open water length. We have chosen a degree day method to determine the mean air temperature. Using the measured air 42 temperatures from Fort St. John assembled by British Columbia Hydro we divide the maximum accumulated degree days by the number of degree days; these "mean" air temperatures are listed in column 4 of Table 5. The minimum measured open water length and the measured average winter discharge (Nov -Feb) were provided by British Columbia Hydro. We have determined mean winter wind velocity and cloud cover from the Meteorological Data for Canada. Humidity was not available, but we have assumed a constant 10% throughout the winter which may be slightly high considering the cold air temperatures 4 at the site. Using this combination of averaged meteorological and hydrological data, we determined Qo and K according to equations 35, 36, and 37 and applied the heat loss coefficients to the appropriate models of open water length, i.e., we used equations 36 and 4 to determine open water length according to Asvall (1972); equations 35 and 9 according to Dingman and Assur (1969); and equations 37 and 13 according to Paily et al. (1974) with E set equal to 0.0. These calculated open water lengths appear in columns 8, 9 and 10 of Table 5. Clearly the Asvall (1972) formulae consistently overpredicts open water length. The Dingman and Assur (1969) and the Paily et al. (1974) formulae are in substantial agreement, as should be expected considering that both heat loss formulae are linearizations of the earlier Dingman et al. (1967) equations. It appears that equation 35 is in better agreement with the data than equation 37, with the former yielding an aver- age deviation from measured open water length of 13 miles and the latter, an average deviation of 15 miles. However, considering the assumptions employed in determining an "average" air temperature, wind speed, discharge etc., the difference is not significant. It is worth noting that the Asvall (1972) formulae overpredicts open water length both in this example and in the earlier example, and that equations 35 and 37 predict the same trend in open water length as is found in the measured open water length. 43 N E T O to to N co C7 O r •" M M ^ M to O Uto CO N E a ,Ny ""1 r-1 v c O 4J N0 L� OI N r c IL l0 Ol Q1 O to d i > J O.M O N N ct 1d x Q ON c W ,�.1 cli ti O O. N O N L O 7 � J N N J M a r C1 to .0 5 2 M 1� M p � Lo U c (a cli t-I N C1 N U 3 W G 14, E 3 s � L 3 Cd O N O M N O > R t9 n lD tp r O U U C Y a• 3 UO N O N Cf E v LO <t to rn N v � i U/ i L p 7 Q 1J n P N O — c L U <t c C� L� IO C/ O •.i 1 .N-1 1 I 1 .E Ems" h C1 O L U M to � Cl) N L ep y O L N eti .-f M N to ON co 1\ N �O NM E .-1 •-1 .N-i ,may .^-1 f O c Gl L O C C C 0 N 4 J Cl L r N MO 00 M O N CD E Gl i0 f 3 L n ^ co 10 \ \ \ \ \ to Y n n n n 4 a Based on the two examples for the Chena River and the Peace River, we recommend either the Dingman and Assur (1969) heat loss expressions (equation 35) or the Paily et al. (1974) expression (equation 37) when a simplified version of surface heat transfer is to be used and when the air temperature is warmer than -19°C. Since the Paily et al. (1974) formulae have not been tested below about -21°C, and since they were derived explicitly for temperatures greater than -19°C, we suggest using equation 35 exclusively whenever air temperatures below -190C are possible. Finite difference methods In the foregoing sections of this report, we have primarily examined steady state and analytic models for the temperature regime in a river. These are important tools for environmental assessment for known meteorological forcing. That is, for design purposes when only the large scale hydrologic conditions and climatic variability are known, the analytic models provide useful estimates of the expected open water length. However for operational purposes on a day to day basis, a finite difference on finite element model is needed to simulate the site specific variations in river hydrology, and variations in discharge and meteorology. General finite difference models for arbitrary surface heat loss have been given by Dingman, Weeks and Yen (1967), by Ashton (1979), and more complicated models for coupled hydrodynamic and thermodynamic analysis have been given by Chaudry et al. (1983) and Bowles et al. (1977). The first model is for steady state conditions and therefore, except for allowing the non -linear surface heat transfer expressions (see equations 26 and 29), offers no real advantage over the analytic models when reliable linearization formulae are used (e.g., equations 35 and 37). The coupled hydrologic thermal models of Chaudry et al. (1983) and Bowles et al. (1977) represent 45 a very sophisticated approach to the analysis of river temperature. However, at this time we do not recommend this level of modeling for application to Alaskan rivers due to the scarcity of the necessary hydrological data. It should be noted that in order to simulate continually varying discharge it is necessary to use the coupled hydrologic thermal models; for gradually changing discharge the thermal models described in this report should give reasonable estimates of the water length. Furthermore, if a coupled hydrologic thermal model were to be used it is essential that the surface heat transfer expressions be based on formulae appropriate for arctic conditions as discussed earlier. It would be necessary to change the thermal portion of the model to follow the suggestions given earlier for surface heat transfer. The finite difference model from Ashton (1979) allows daily variations in meteorology and local variations in river width and mean velocity. Variations in discharge, both from changes at the dam and from stream inflow downstream are not included in the Ashton (1979) mo- del; ice dynamics are also not included. However, the Ashton (1979) model provides a useful framework for the study of transient effects, and is easily modified to include a variety of site specific adaptations. A copy of the Ashton (1979) computer model is included in Appendix A. In the present section, we shall briefly describe the model, its limitations and assumptions and discuss refinements which could be included. The Ashton (1979) model is a numerical solution of equation 2 in which the width b, depth h and mean velocity U are allowed to vary with downstream distance x; Q the surface heat exchange is calculated according to a simplified air temperature -wind velocity formula, and no other heat transfer terms S' are included. The river discharge (D = Uhb) is assumed to remain constant over the calculation period. The simulations are done in a M Lagrangian reference frame, following a fluid parcel downstream; therefore downstream distance steps are set internally depending on local velocity. The inclusion of local river geometry in terms of variable U, h and b is clearly an improvement over analytic models in which these terms are held constant. This feature is important for the heat balance since the net heat loss is directly proportional to river width b. presently, Ashton (1979) assumes that the river hydrology h(x), b(x) and D are known, and calculates U(x) locally assuming a rectangular basin. In principle, any measured river width and depth distribution, including the total river width in a braided section of the river, may be used as data. The extension of the model for alternate basin geometry (e.g., trapezoidal or multi -channel) is straightforward requiring only the inclusion of a flag variable to define basin geometry at each subreach (alter statement 40 in the model to define area discharge relation and statements 10 and 20 to define the basin geometry flag). In its present form, the model features a simplified expression for Q which is calculated daily based on mean air temperature and wind velocity. For application to Alaskan rivers we recommend using the linearized expressions from Dingman and Assur (1969) (see equation 35). These expressions are only slightly more complicated than those in the Ashton model, and programming changes to the model would be minimal (alter statement 87). If small streams enter the main channel, they will increase river discharge and alter the thermal balance. In principle, this effect can be handled by solution of equation (1) in which the other heat sources are the known stream input in terms of stream water and ice discharge and water temperature. The solution procedure will "step downstream", and a new 47 increased discharge calculated for the next calculation reach. As a practical matter, information on small stream discharge, temperature and particularly ice content, is usually not available. In addition, if the stream inflow is at a different temperature from the river, it will also be at a different density, and subsequently will not mix instantaneously with the main flow. However, a reasonably good literature exists regarding theoretical and field examinations of transverse mixing in rivers, and a rough estimate can be made of the distance required for complete transverse mixing. If this distance is substantially°less than the estimated open water length (from equation 9), then the thermal effects of inflowing streams can be simulated by adding discrete heat and mass sources to the governing equations at the appropriate locations. If the mixing distance is of the same order as the open water distance, then a two-dimensional model involving downstream convection of heat and cross -stream diffusion of heat must be used. An example of such a two-dimensional model is given by Ashton (1979) and is listed in Appendix B. The determination of whether a two-dimensional model is required hinges on the estimate for transverse mixing length Lt. Transverse mixing for open channel flow is determined by the transverse mixing coefficient et where et = chU* and h is depth, U* is friction velocity and c is a scale constant (Fischer et al., 1978). Ashton (1979) assumes c - 0.2, but a more recent compilation of typical values suggests c = 0.6 is more appropriate for the winding rivers characteristic of Alaska. Diffusion theory predicts that a passive tracer will diffuse as (time)1/2: /N where a2 is the variance of the diffusion and a is the appropriate diffusion coefficient. Fischer et al. (1979) suggest that a reasonable criterion for substantially complete transverse mixing is when the tracer is diffused to within 5% of its mean value everywhere on the cross-section. Assuming a Gaussian distribution for the tracer, this occurs when o = 0.5b where b is river width. The time required for this to occur following a fluid parcel is Tt m a2/et = 0.25 b2/et, and the downstream distance travelled is Lt = TtU = .25b2U/et = 0.25 b2U/(.6U*h). Since a reasonable approximate value for U* - 0.1U, we have (38) Lt = 4b2/h As a rule of thumb, the river and side stream inflow are well mixed at the distance Lt. If this distance is the same order of magnitude as L the estimated open water length from equation 9, then a two-dimensional model is necessary. If on the other hand Lt < 0.1 L, a one-dimensional model is acceptable. Assuming that Lt << L for all small streams entering the main river, a procedure could be devised to alter the one-dimensional Ashton model (1979) to include these additional thermal sources. The simplest way to do this appears to be: first, make_ discharge a variable (alter statements 8, 12, 15, 30, 40) in particular defining the subreach velocity by the reach characteristics (statement 40 becomes U(J) = DISCH(I)/(SB(I)* SD(I))), and second, define a new variable giving the temperature increment from the small stream and insert it where the subreach characteristics are defined, say after statement 46. It would be of the form TINC(J) = TINFL(I)*(DISCH(I) -DISCH(I-1))/DISCH(I), and TINFL is temperature of the stream water. This would have the effect of adding the additional heat only where the stream enters and weighting it proportional to the stream discharge. Finally, the third step would require that TINC(J) be added to the local temperature, by altering statement 80 to read TWOUT(J) = TWOUT(J) + DELTW + TINC(J). A few additional alterations would be required to change format statements, and to zero unaffected TINC(J), etc. A discussion of the Ashton (1979) finite difference model is incomplete without reference to the two modes of thermal equilibrium used in the model. The first mode states that length of the open water reach is coincidental with the position of the zero degree (°C) water temperature. This is the mode that has been assumed throughout this report, and is implicit in the analytic solutions (see equations 9 and 13). The Ashton model uses this definition (Tw(L) = 0°C) to define L when the ice cover is newly forming or melting. When an ice cover is already present, an alternative criterion for ice edge position is adopted in the Ashton (1979) model which is referred to as the equilibrium criterion. The equilibrium criterion is derived from the heat balance equation through the ice cover: (39) (Tm-Tair)/(n/ki + 1/hia) - hiw (Tw - Tm) = Pi a do/dt where n is ice thickness, Tm is the melting point (Tm = 0°C), Tair is air temperature, Tw is water temperature, ki is thermal conductivity of the ice, Pi is ice density, a is the heat of fusion and hia and hiw are the ice/air and ice/water heat transfer coefficients respectively. This equation in turn is derived from the energy balance at the water/ice interface: (40) �i - �wi = Pi a do/dt 50 where �i is the heat flux by conduction through the ice and Owi is the heat flux from the water to the ice. It is assumed that �i = pia where pia is the heat flux from the top surface of the ice to the atmosphere. The equilibrium criterion for the leading edge of the ice is determined from the condition that n = do/dt = 0 in equation (39). This condition then defines the equilibrium temperature of the water at the leading edge: (41) Twe = -hia/hiw Tair Clearly TWe is not in general equal to 0°C; in fact, TWe < 0 whenever Tair > 0 and TWe > 0 whenever Tair < 0. The first condition is clearly meaningless and therefore in the model the equilibrium criterion is inoperative whenever Tair > 0; the more standard zero isotherm criterion is adopted for the position of the ice edge if Tair > 0. The equilibrium criterion is used in the model only when a presently existly ice edge is growing or decreasing in length and the air temperature is less than zero; under all other conditions including the first formation of the ice, the zero isotherm criterion is used. There are several basic problems associated with the use of the equilibrium criterion. This criterion is determined from equation 40 with the additional assumption that the conductive heat transfer through the ice exactly balances an expression for ice/atmosphere heat transfer. It should be noted that: 1) equation 40 neglects the possibility of surface melt, defining all melting on the water/ice interface; 2) the expression used for conductive heat transfer across the ice is the steady state linear formula (�i =-kiTs/n where Ts is the top surface temperature of the ice) which is not realistic during a period of ice growth or decay; 3) the expression used for the ice/atmosphere heat transfer pia is a simple linearization 51 formula (Oia = hia (Ts - Tair)) and thus Oia effectively ignores effects of melt puddles and short wave radiative exchange; 4) equating pia to conductive heat transfer 01 is a questionable assumption, particularly when the ice is wet and Ts is close to DOC while Tair « 0 °C; and finally, 5) there are no data available which would indicate that the equilibrium criterion is actually an improvement on the zero isotherm criterion. The zero isotherm criterion may be implemented as the only criterion by the following program modification. Between statements 105 and 106 add the statement, IF (ETA(J).GT.O. .AND. TWOUT(J).GT.O.) ETA(J) = 0. Finally, note a correction to the Ashton model; statements 103 and 104 should be reversed. 52 Conclusions This report has reviewed several approaches to the problem of the determination of the length of open water downstream from a dam or thermal source in winter. Open water lengths have been predicted by several Russian studies by statistical approaches based on local data, and are appropriate only for particular locations. More general types of analyses for rivers were introduced by Asvall (1972), Dingman, Weeks and Yen (1967), Paily et al. (1974) and Harleman (1972). These analyses are based on semiempirical formulae for the rate of heat transfer from an open water surface to the atmosphere by evaporation, radiation and sensible heat transfer, and possibly including infiltration of ground water and frictional effects. The heat transfer expressions are applied to the one-dimensional equation for conservation of thermal energy in the river, yielding solutions which predict temperature in the river. Since the formulae for radiative heat transfer are non -linear functions of water temperature, in general numerical methods must be used to determine temperature distributions in the river. However, there are several "linearized" versions of the surface heat transfer expressions, including those by Paily et al. (1974) and by Dingman and Assur (1969). The application of these linear heat transfer expressions greatly simplifies the mathematics involved in the determination of river temperatures, and in fact, allows closed form analytic solutions to be found for a limited number of boundary conditions. The most obvious of these analytic solutions is the steady-state case, given by Dingman and Assur (1969) and defined in this report by equation 9. 53 We compared the steady-state solution with measured ice -free area in the Chena River and in the Peace River for three linearizations of the surface heat transfer expressions: Dingman and Assur (1969); Paily et al. (1974); and Asvall (1972). The linearizations given by Dingman and Assur (1969) and by Paily et al. (1974) were based on the "Russian winter equation" of Rimsha-Donchenko (1957) and produced the best agreement with the data. However, since the Paily et al. (1974) linearization formulae were derived primarily for air temperatures greater than -19°C, the Dingman and Assur (1969) formulae given by equation 35 are recommended. Paily et al. (1974) found an additional analytic solution, the transient response of an intially uniform river temperature distribution to a given temperature increment at x = o, with constant air temperature and solar radiation. In equation 20 of this report we introduce a new analytic solution, the transient response of river temperature to periodic air temperature and/or solar radiation. The latter analytic solution provides information on the phase lag between atmospheric forcing and river response, indicating that this lag increases with increasing river depth and decreases with surface heat loss rate and is independent of river width. This closed form solution also includes the effects of spatially varying air temperature and therefore, provides a general model for temperature prediction in rivers with uniform flow and uniform cross - sectional area. Another important use of this transient analytic solution is for comparison with numerical models. Since the analytic solution is exact, it provides a reliable gauge for the accuracy of finite difference or finite element models, thus providing confidence in the applicability of these models. 54 General finite difference models for arbitrary surface heat loss and changing river basin geometry have been given by Dingman et al. (1967), Ashton (1979), Chaudry et al. (1983) and Bowles et al. (1977). For general applicability in Alaskan rivers where the hydrological data base is sparse, we recommend the Ashton (1979) finite difference model for river temperature analysis. The model predicts the transient response of water temperature for constant discharge, spatially varying cross - sectional area, and temporally varying air temperature and discharge temperature. We have discussed several refinements to Ashton model including arbitrary (non -rectangular) cross -sectional area, the implementation of the Dingman and Assur (1969) heat transfer expressions, heat flux from small streams and an alternative criterion for leading ice edge position. Any of these modifications may be rather simply applied to the Ashton (1979) model. None of the models discussed in this report are applicable when river discharge is changing drastically. In this case, ice movement and ice front position is a mechanical -hydrodynamic problem, only slightly affected by thermal changes.' At this time there is no reliable theore- tical or numerical model available for ice front behavior with rapid changes in discharge. For clear strategic reasons, field measurements of these events are rare. In this report we have reviewed the thermo- dynamic models which are appropriate only for gradually changing discharge when the ice conditions and water temperatures are controlled by the discharge temperature and the local meteorology. Comparisons with data indicate that under these conditions, the appropriate thermodynamic models yield realistic estimates of open water length. 55 REFERENCES Anderson, E. R., 1952. Energy Budget Studies from Water Loss Investigations, Vol. 1, Lake Hefner Studies Technical Report, U.S. Geological Survey Circular #229, pp. 91-119, Washington, D.C. Angstrom, A., 1920. Application of Heat Radiation Measurements to the Problems of the Evaporation from Lakes and the Heat Convection at Their Surfaces, Geografiska Annales, Vol. 2, pp. 237-252. Ashton, G., 1979. Suppression of River Ice by Thermal Effluents, CRREL Report 79-30. Asvall, 1972. Proceedings of the Banff Symposia on The Role of Snow and Ice in Hydrology. Bolsenga, S. J., 1964. Daily Sums of Global Radiation for Cloudless Skies, U.S. Army CRREL Report 160, Hanover, N.H. Bowles, D. S., 1977. Coupled Dynamic-Streamflow-Temperature models, ASCE, Journal of the Hydraulics Division, Vol. 103, No. HYS, pp. 515-530. Brunt, D., 1944. Physical and Dynamical Meteorology, Cambridge University Press. Brutsaert, W. and G. Jirka (editors), 1984. Gas Transport at Water Surfaces, D. Reidel Publishing Co., Dordrecht, Holland. Carlson, R. F., T. Tilsworth and C. Hok, 1978. Effects of Thermal Discharge Upon a Subarctic Stream, Institute of Water Resources-49, University of Alaska, Fairbanks, Alaska. Chaudry, M. H., 1983. Modeling of Unsteady Flow Water Temperatures, ASCE, Journal of Hydraulic Engineering, Vol. 10.9, No. 5, pp. 657-669. Daily, J. W. and D. Harleman, 1966. Fluid Dynamics, Addison-Wesley Publishing Co., Reading Mass. Devik, 0., 1964. Present Experience on Ice Problems Connected with the Utilization of Water Power in Norway, Journal of the International Association for Hydraulic Research, Vol. 2, No. 1. Dingman, S. and A. Assur, 1969. Cold Regions Research and Engineering Laboratory Research Report 206, Part II, Hanover, N.H. Dingman, S., W. Weeks and Y. C. Yen, 1967. The Effects of Thermal Pollution in River Ice Conditions -I: A General Method of Calculation, U.S. Army Cold Regions Research and Engineering Laboratory Research Report 206, Pt. 1. Elsasser, W. M., 1942. Heat Transfer by Infrared Radiation in the Atmosphere, Harvard Meteorological Study No. 6, Harvard Univ. Blue Hill Met. Observatory. F: Fischer, H., J. Imberger, E. J. List, R. Koh and N. Brooks, 1979. Mixing in Inland and Coastal Waters, Academic Press, N.Y. Friehe, C. A. and K. F. Schmitt, 1976. Parameterization of air-sea inter- face fluxes of sensible heat and moisture by the bulk aerodynamic formulas, J. of Physical Oceanography, 61 801-809. Goryunov and Perzhinskiy, 1967. Soviet Hydrology, Selected Papers, Issue No.- 4. Gotlib and Gorina, 1974. Gidrotekhnicheske, Stroitel`stvo, No. 11. Harleman, D.R.F., 1972. Longitudinal Temperature Distribution in Rivers and Estuaries -One Dimensional Mathematical Models, Engineering As ects of Heat Disposal from Power Generation -MIT summer session, M. . ambri ge, ass. Henderson, F. M., 1966, Open Channel Flow, McMillan Publishing Co., N.Y. Hicks, P. B., 1975. A Procedure for the Formulation of Bulk Transfer Coefficients Over Water, Boundary Layer Meteorology, 8, pp. 515-524. Hicks, P. B., 1972. Some Evaluations of Drag and Bulk Transfer Coefficients Over Water Bodies of Different Sizes, Boundary Layer Meteorology, 3, pp. 201-213. — Holmgren, B. and G. Weller, 1974. Local Radiation Fluxes Over Open and Freezing Leads in the Polar Ice Pack, AIDJEX Bulletin, 27, pp. 149-166. Ince, S. and G. Ashe, 1964. Observations on the Winter Temperature Structure of the St Lawrence River, Proceedings of the Eastern Snow Conference, p. 1-13. Kays, W. M., 1966. Convective Heat and Mass Transfer, McGraw Hill, New York. Kohler, M. A., 1954. Lake and Pan Evaporation, in Water -Loss Investigations: Lake Hefner Studies, U.S. Geological Survey Professional Paper 269, PP• - McFadden, T., 1974. Suppression of Ice Fog from Power Plant Cooling Ponds, Ph.D. Dissertation, University of Alaska, Fairbanks, Alaska. Mellor, M., 1964. Properties of Snow, CRREL Monographs, Part III, Section Al. Osterkamp, T. E., R. E. Gilfilian, J. P. Gosink and C. S. Benson, 1983. Water Temperature Measurements in Turbulent Streams During Periods of Frazil Ice Formation in Rivers, Annals of Glaciology, 4, pp. 209-215. — Osterkamp, T. E., K. Kawasaki and J. P. Gosink, 1983. Shallow Magnetic Induction Meausrements for Delineating Near -Surface Hot Groundwater Sources in Alaskan Geothermal Areas, J. Energy Res. Tech., 105, pp. 156-161. 57 Paily, P., E. Macagno and J. F. Kennedy, 1974. Winter -Regime Thermal Response of Heated Streams, ASCE Hydraulics Division, HY4, pp. 531-550. Pasquill, F., 1948. Eddy Diffusion of Water and Heat Near the Ground, Proceedings of the Royal Society (London), Ser. A., Vol. 198, No. 1052, pp. 116-140. Pond, S., D. B. Fissel and C. A. Paulson, 1974. A Note on Bulk Aerodynamic Coefficients for Sensible Heat and Moisture Fluxes, Boundary -Layer Meteorology, 6, pp. 333-339. Pruden, F. W., R. T. Wardlaw, D. C. Baxter and I. T. Orr, 1954. A Study of Wintertime Heat Losses From a Water Surface and of Heat Conservation and Heat Addition to Combat Ice Formation in the St. Lawrence River, Report No. MD42, National Research Council of Canada, Ottawa, Canada. 4 Rimsha, V. A. and Donchenko, R. V., 1957. The Investigation of Heat Loss From Free Water Surfaces in Winter Time. Tr. Leningrad Gas. Gidrol. Inst. 65 (in Russian). Schlichting, H., 1968. Boundary Layer Theory, McGraw Hill, N.Y. Starosolszky, 0., 1970. Ice in Hydraulic Engineering, Report No. 70-1, Division of Hydraulic Engineering, The Norwegian Institute of Technology, The University of Trondheim. Tennessee Valley Authority, 1972. Heat and Mass Transfer Between a Water Surface and the Atmosphere, Division of Water Resources Research, TVA Report No. 14. Wendler, G., 1980. Solar Radiation Data for Fairbanks, Solar Energy Meteorological Research and Training Site, Geophysical Institute, University of Alaska, Fairbanks, Alaska. 0 Table 1 - Models of open water area 1) Statistical Models Investigator Relationship for open water area Gotlib and Gorina (1974) 1) Graphical relationship between length Godrotekhnicheskoe of open water (L), reservoir discharge Stroitel'stvo, No. 11 (D) and temperature of discharge (Tw) L = f (D, Tw) for cold winters 2) For warm -winter conditions, one specific L = f(D) Goryunov and Perzkinskiy (1967) Soviet Hydrology; Selected Papers, Issue No. 4 2) Semi -Empirical Models Investigator Asvall (1972) Proc. of Banff Symposia on the role of snow and ice in hydrology 3) Graphs for transient response of location of ice edge under warming conditions L (t)/L (t=o) = f(D,ni,E-lair) where ni = ice thickness, Tair = air temp. (Note: This graph is not well labeled. The location of the ice edge is not de- fined clearly. Only qualitative infor- mation regarding the effects of variation in D and ni may be discerned). L = 5.5 • 106 (E - Tair)-2 L gives the transient location of ice edge since t - Tair sums over all negative degree days Surface Heat Loss Definition Tabular values heat loss (Q), as a function of cloud cover (C), wind velocity (w) and air temp. (Tair); Q = f (C, w, Tair)• Q can be formulated, Q = ao + al Tair where ao, al = f (C,w) Relationship for Open Water Area Area (Lb) times sur- face loss (Q) equals heat input from reser- voir LbQ = Ro L = Ro/Qb E'�1 Surface Heat Relationship for Investigator Loss Definition Open Water Area Paily, Magagno, Kennedy Q = -K (T-TE) Solves the one- (1974), Jrnl. of where K is a surface exchange dimensional partial Hydraulics Div., ASCE coefficient diff. equation for T is local water temperature conservation of TE is "equilibrium temp"., thermal energy. It water temp. at which there is is assumed that no exchange of heat across T = T(x,t) only, and the water surface with the the equation is atmosphere. integrated over a cross -sectional area. Since a linear re- ation is assumed for tie equation may be solved analy- tically, yielding an expression T = T (x,t) Dingman, Weeks, Yen Q = non -linear function of Numerically integrates (1967), CRREL Res. T, Tair, w, e, D, S where the one-dimensional, Rpt. 206 new variables are steady-state partial e = evaporation pressures differential equation for D = discharge conservation of thermal S = river slope or energy assuming negli- Q = non -linear function of gible longitudinal T, Tair. w, D, S diffusion. Dingman, Assur (1969) Q = Q'o + K (T - Tair) Closed form solution of CRREL Res. Rpt. 206 Part II Q'o and K from regression one-dimensional steady - analysis of the non -linear state ordinary differ - function in Dingman et al. ential equation (linear) (1967). for conservation of thermal energy. Table 2 - Surface Heat Transfer Definitions Investigator Surface Heat Transfer Expressions in W/m2 Asvall (1972) Cloud cover = 0.0, Q = 136.05 + 2.09 w +(12.59 + 1.63 w)ITa = 0.5, Q = 77.38 + 2.09 w +(9.44 + 2.41 w)IT r 1 0, Q = 23.00 + 2.09 w + (10.92 + 2.05 w)T}a where Q = [W/m�], w = [m/sec] Tair = ['C] Paily, Macagno, Graphs of a and n for Kennedy (1974) Q = eT+n where e = e (Tair, w, R.H) n = n (Tair, w, R.H) and R.H. = relative humidity It is assumed that barometric pressure = 99.6 mb cloud height = 3,275 ft cloud cover = 6 visibility = 1.87 miles Dingman, Weeks, Yen and QR - QB - QE - QH - QS + QG + QGW + QF QR is heat from short wave radiation QB net loss of heat by exchange of long -wave rad. w. atmos. QE heat loss due to evaporation QH sensible heat loss QS heat lost by influx of snow QG heat added by flow of geothermal heat QGW heat added by flow of ground water QF heat added by friction on stream bottom QR = QRI QRR = incoming -reflected short wave radiation and QRI = QCL 1.17 + .30 (1-C)] QCL is incoming short wave radiation C is cloudiness in to the 2 QRR = .052 QRI - 3.28 • 10-� QRI 'QB = Qa - Qar - Qbs Qa = long wave radiation from atmosphere Qar = .03 Qa = reflected incoming long wave radiation Qbs = long, wave radiation from water surface 61 Qa = (a + bea) a Tair4 a = .36 + .12 C exp [-1.92 • 10-4 Z] b = 2.8 • 10-3 - 26.1 • 10-4 C exp [-1.97 10-4Z] ea = vapor pressure of air (mb) Z = cloud height (m) Qar = .03 Qa o Qbs = •97 a Tw QE and QH were estimated by two approaches. Kohler formulae: QE = (1.52 + 3.55 w) (eSW - ea) QH = (.92 + 2.16w) (Tw - Tair) w wind velocity at 2 m. Rimsha and Donchenko formula QE = (1.56 kN + 2.94w) (eSW - ea) QH = (kN + 1.89w)(Tw - Tair) kN = 3.87 + .17 (Tw - Tair) eSW = saturation vapor pressure (mb) Qs = A [a + Ci (Tw - Tair)] A is snow ar�cpglation rate A . 7.85 V- V is visibility in km Ci heat capacity of ice a latent heat of ice QG by local measurements of geothermal gradient QGW by local measurements of ground water flows QF = OyS/Jb D is river discharge [m3/sec] y is weight density of water = pwg = [kg/m2 sec21 S is water surface slope b is river width Dingman, Assur (1969) Clear (C=O.) Q = 50.93 + 11.21w + (16.99 + 2.05w) (Tw-Tair) Cloudy (C=1.0) Q = -35.28 + 4.40w + (17.97 + 2.22w) (Tw-Tair w = wind velocity Tw = local water temperature, Tw (x) Tair = ambient air temperature 62 APPENDIX 5 REPORTS Osterkamp, T. E., J. P. Gosink, K. Kawasaki and G. Penn, "Annotated biblio- graphy listing sources of information for small hydropower users in cold climates", An Interim Report to the Alaska Energy Center under Contract No. AEC81005-3, December 1981. Osterkamp, T. E., J. P. Gosink, K. Kawasaki and G. Penn, "Survey of manu- facturers of hydroelectric equipment", An Interim Report to the Alaska Energy Center under Contract No. AEC81-005-3, December, 1981. Osterkamp, T. E. and J. P. Gosink, Ice cover development on interior Alaska Streams, Report No. UAGR-293, Geophysical Institute, University of Alaska, Fairbanks, Alaska, December, 1982. PAPERS PRESENTED AT MEETINGS Gosink, J. P. and T. E. Osterkamp, Hydraulic resistance generated by frazil ice formation, Paper presented at the Workshop on Hydraulic Resistance of River Ice, Burlington, Ontario, Canada, September, 1980. Osterkamp, T. E., R. E. Gilfilian, J. P. Gosink and C. S. Benson, Water temperature measurements in turbulent streams during periods of frazil ice formation. Paper presented at the Second Symposium on Applied Glaciology, Hanover, NH, August 23-27, 1982. Gosink, J. P. and T. E. Osterkamp, Measurements and analyses of velocity profiles and frazil ice crystal rise velocities during periods of frazil ice formation on rivers. Paper presented at the Second Symposium on Applied Glaciology, Hanover, NH, August 23-27, 1982. Gosink, J. P. and T. E. Osterkamp, Preliminary evaluation of hydroelectric power generation in cold climates, Paper presented at the 33rd Alaska Science Conference, AAAS, September 1982, Fairbanks, Alaska. Gosink, J. P., T. E. Osterkamp and P. A. Hoffman, Modeling of ice covered lakes. Poster session presented at Frontiers in Hydrology Speciality Conference of ASCE, M.I.T., Cambridge, Massachusetts, August 1983.