The freshwater budget and underice spreading ... - Wiley Online Library

2 downloads 0 Views 2MB Size Report
Jan 15, 1995 - construct a freshwater budget for the Mackenzie shelf, the nearshere within which ... Table 1. the Shelf. Inflow of the Mackenzie River and Area of ..... 0. 10. 20. 30. 40. Salinity (psu). Figure 4. Salinity versus 5180 for the water column data collected in ... Kugmallit Bay, 181-182 cm (n - 24) at Sachs Harbour,.
JOURNAL OF GEOPHYSICAL

RESEARCH,

VOL. 100, NO. C1, PAGES 895-919, JANUARY

15, 1995

The freshwater budget and under-ice spreading of Mackenzie

River

water

in the

based on salinity and and

Canadian

Beaufort

Sea

measurementsin water

ice

Robie W. Macdonald, David W. Paten, and Eddy C. Carmack Institute of Ocean Sciences, Sidney, British Columbia, Canada

Anders

Omstedt

Swedish Meteorological and Hydrological Institute, Norrk;Sping

Abstract. Observations of salinityand oxygenisotopecomposition (5180) were made for the Beaufort shelf-Mackenzie estuary waters in September 1990, just prior to ice formation, and for both the water column and ice in April-May 1991, at the end of winter. These measurementsare used to determine the apportioning of fresh water in the estuary between its two main souxces,runoff and sea ice xttex[. •nanges in dispositionof water between seasonsand amounts frozen into the growing ice sheet are also derived. Two domains are consideredin order to constructa freshwaterbudget for the Mackenzie shelf, the nearsherewithin which

landfast ice growsin winter and the outer shelf. Most of the winter inflow from the Mackenzie River appears to remain impounded as liquid under the ice within

the landfast zone at the end of winter, and about 15% of it is incorporatedinto

the landfastice. Oxygenisotopes (Slso) in ice corescollected from acrossthe shelf record the progressbeneath the ice of new Mackenzie inflow as it invades the

nearsherethroughoutwinter. Rates of spreadingare about 0.2 cm s-• away from the coastand 1.3 cm s-1 alongthe coast. As this inflow spreadsacrossthe shelf, it progressivelyshuts off convectiondriven by brine production at locations within

the landfastice. Salinityand 5•80 in the offshorewatercolumnsuggest that about 3 m of sea ice was formed in the outer shelf domain. Since both brine and newly formed sea ice can be advectedoff the shelf, a completebudget for brine or sea ice production cannot be establishedwithout first measuring the advection of one of these two components.

1. Introduction

posedlarge-scaleoil developmentboth in the Canadian

Recent oceanographicstudies of the Arctic tend to be motivated by issuesrelating to climate changeand contaminants. Within the Arctic Ocean the relatively dynamic margins, which make up about a third of the surfacearea, are arguably the most important regions, both for their physicaland chemicalprocessesand their biology. It is at thesemarginsthat changesin ice thicknessand duration of the ice seasonare expected to be most pronounced,should global warming occur. The shelvesare also the logxcalentry point for many contaminants becauserivers drain into them, becausethey provide large open areas with which to exchangewith the atmosphere,and becausethey are the sites of pro-

and Russian sectors of the Arctic. Processes that control the under-ice

flow of river wa-

ter in major Arctic estuariesstronglyinfluencethe wintertime oceanographyand transport pathwaysover the

shelves[Macdonaldand Carmack,1991].Fewmeasurements of the behavior of river water plumes under ice have been made owinglargely to instrument limitations

at the slowspreadingvelocitiesinvolvedlingyam,198!; Ingram and Zarouche,1987]. In the presentpaper we use observationsof oxygen isotope ratios and salinity to determine the evolution and spreading rates of the under-ice plume arising from the Mackenzie River inflow to the Beaufort shelf from October 1990 to May 1991.

Arctic shelves are important sources of both posiCopyright1995by the AmericanGeophysical Union.

tive (f•e•h w•t•)and negative(brine)buoyancyfluxes [Macdonaldand Carmack,1991].On the shelves,input

Papernumber94JC02700.

of fresh water by river inflow, ice melt, and net precipitation increasesthe averall stability of Arctic sur-

0148-0227/95/94JC-02700505.00 895

MACDONALDET AL.: FRESHWATERBUDGET OF MACKENZIERIVER

896

Table 1. the Shelf

Paramete• Total Mackenzie

only the fourth largest river entering the Arctic Ocean

Value

River flow 1990

MackenzieRiverinflow(October-May) River waterfrozen(winter) Net inflow(winter) Total

Arctic shelves[Macdonaldet al., 1987; Thomaset ai., 1986;Aa#aardand Carmack,1989]. The Mackenzieis

Inflow of the Mackenzie River and Area of

Area of shelf

Area of landfast ice Area of outer sheff

(Table 1) (seealsoAa#aardand Carmack[1989]),but

sincetheadjacent shelf.•srelatively small(64,000 km2,

300 km $ 73 km s 11 km s 62 km s

Figure 1) the annual freshwateryield for the shelfis about 5-6 m [Brunskill,1986;Macdonaldet al., 1987]. Becausethe Mackenzie ILlvet draws on large headwater

64,000km2 16,000km2 48,000km2

lakes,the winterdischarge remainshigh (3500m3 s-x) and delivers amounts of water comparable to the Ob

(4000m3 s-x) andYenisey(5500m3 s-x) (V.K. Parlor and S.L. Pfirmen, Hydrographicstructure and variability of the Kara Sea: Implications for pollutant distri-

bution, submittedto Deep Sea Research,1994). De-

facewaters,therebypromotingwinterice growthand suppressing deepconvection, both in the centralArctic and, possibly,by advectionout of the system,in the NorthAtlantic[Aagaard and Carmack,1989].At the sametime,cold,brine-enriched watersproduced by freezingseawater serveto coolanddeepenthe Arctic halocline,therebyfavoringfurther ice growthby sustaininga thick,cold,insulating layerto shieldthe surface,whereiceis formed,fromthe deeperwaterwhich contains heat[Aagaard et al., 1981;MellingandLewis, 1982].

spite this strong, year-round supply of buoyancy the Mackenzie shelf is known in certain winters to produce sufficient brine-enhanced water from sea ice formation

to ventilate the offshorehalocline[Melling and Lewis,

1982;MellingandMoore,•995]. Two circumstances appear to be critical for this to occuras follows:first, the upperlayersmust be preconditionedby removalof freshwatercomponentsprior to freeze-up[Macdonald et al., 1989; Melling and Moore, 1995; Omstedtet al., 1994];and second,the dilutinginfluenceof new winter inflow of river water must be blocked from areas fa-

Just as the Arctic Ocean, itself, acts as an estuary vorableto brine enrichment[Macdonaldand Carmack, for the world ocean, so do the shelvesact as estuarThe shelves,which behave as a short-term buffer, ies for the Arctic Ocean, and the Mackenzieshelf of the Canadian Beaufort Sea is the most estuarine of all modify fresh water before it is exported to the Arc-

71øN

ß 30 ß 29

200m

ß 50

ß28

ß 52 70 ø

ß:24 ß 23 ß 22 ß21

49

38 ß36 35

69 ø

140øW

135ø

130 ø

Figure1. Bathymetry andstation locations for(a) thefallof1990showing thelocation ofthe packiceat thetimeof sampling (patterned area)and(b) thespring of 1991showing thelarge flawleadthat startedto openin late Februaryof 1991(cross-hatched area).

MACDONALD

ET AL.:

FRESHWATER

BUDGET

OF MACKENZIE

RIVER

897

50 m

(b) 71øN 200 m

CD-4' CD-3'

C-6 ß

CB-7 CB-40

CD-2 ' L-3

50 rn

70 ø

'5 ' CM-4 ß

PI-•

L-1 ß

PI-2 ß PI-!

ß GI-2

AP-3. AP-2 ' AP-1

' CM-3

CM-2 ' CM-I'

ß GI-1

....... ktoyaktuk

140øW 130 ø

135 ø

Figure 1. (continue. d)

tic interior water

ocean.

on Arctic

Therefore

the residence

shelves and within

time of fresh

the Arctic

halocline

progresses,the shelf clears of ice so that by September

the edgeof the permanentpack is well offshore(Fig-

[(•stlund, 1982;M•cdon•ld et •l., 1989;Schlosser et •l., ures la and 2a). During the open-waterseasonthe 1994] and the effect of freshwaterrunoff on the heat Mackenzie plume dominates the nearshore. Autumn budget of the Arctic Ocean and global thermohaline (Figure 2b) is characterized by decliningriver inflow circulation have been• and will continue to be, impor- and intensewinds[Fisseland Melling, 1990]whichtotant topicsof investigation(see,for example,the new gether produce well-mixed, relatively saline nearshore program proposedby the Surface Heat Budget of the water [Carmacket al., 1989]. New ice usuallybegins

ArcticOcean(SHEBA) ScienceWorkingGroup[1994]). to form in October(Figure 2c) with relativelysmooth Models

of Arctic shelves tend to treat them as a homolandfast ice occupyingthe region inshore of the 20-m geneouscompartment into which runoff is added, ice isobath(Figure lb). The mostimportantattribute of formed or melted, and for which there is some charac- the landfast ice in the context of this paper is that it reteristic residencetime or flushingtime [cf. Macdonald mains in place throughout winter, although ridges may et al., 1989]. Recently,Omstedtet al. [1994]modeled form within it. At the beginningof the ice-growingseafreshwater cycling on the Mackenzie shelf by parame- son, ice forming in shallow regionsnear the mouths of terizing the fluxes of river inflow, ice advection, wind- the Mackenzie River incorporates fresh water from the driven transport, and vertical mixing. Because this plume which is spreadingbeneath it, while farther offmodel uses only a single, vertically resolvedreservoir shore, sea ice formation is accompaniedby brine rejeccoupled to the open ocean, it cannot describe the lat- tion (Figure2c). Duringwintera flawlead opensspoeral transport within the estuary. To construct a model radically at the outer edge of the landfast ice, usually that allows a shelf simultaneouslyto accept fresh water at about the 30-m isobath, in responseto wind-driven and producebrine requiresunderstandingthe distribu- movementof the offshorepackice (Figures2c and 2d)

tion of freshwaterwithin the system(ice,ice melt, and [Giovandoand Herlinveau•, 1981]. The flaw lead is a river runoff) and the spatialand temporaldisposition region of enhancedice and brine production in winter. of important storage compartmentson the shelf. The seasonalcycle for the Mackenzie shelf is shown conceptually in Figure 2. In summer the Mackenzie River dischargesabout 75% of its annual inflow be-

As winter advances,the seaward edge of the landfast ice becomesprogressivelymore heavily ridged, and ex-

tensiverubblefieldsare produced(Figure2d). These rubblefields(stamukhi)act as an inverteddam to trap

tweenJune and September(Figure 3). At the start lessdensefresh water in the nearshore. The landfast ice, of freshet the shelf is mostly ice covered. As summer

which consistsof large, fiat fields separated by relatively

898

MACDONALD

ET AL.:

FRESHWATER

BUDGET

(b}

OF MACKENZIE

RIVER

W•ND

(c)

is$ light

is$ heavy

........

Figure 2. A schematicdiagramof the runoff and ice cyclein the Mackenzieestuaryshowing

(a) summerfreshet,characterized by stronginflowfromthe Mackenzie River,plumesandfronts on the shelf,and the permanentice pack,locatedwell offshore;(b) autumn, characterized by declininginflow, and autumn stormswhich mix and redistributefresh water on a largely icefree shelf;(c) early winter, characterized by freezingtemperatures with the Mackenzieinflow beginningto spreadunderthe newlyforminglandfastice, leavingits imprint asisotopicallylight ice(hatchedarea). Brine(circleB) produced in theouterregions wherethewinterinflowfromthe MackenzieRiverhasnot penetrated,helpsto alestabilize and convectthe watercolumn;and (d) the end of winter just beforebreakup where the winter inflow is containedbehind the stamukhi

zone(SZ) that hasprogzessively grownthroughout winter,at the endof the landfastice. Beyond the stamukhi,enhancedbrineproductioncontinues in the flaw lead (FL). The growingsheetof icein Figures2c and 2d hasincorporatedvaryingquantitiesof river wateras evidencedby a light

oxygenisotopecomposition (•sO).

low-relief ridges, has the potential to record horizontal shelf more than 97% of the meteoric water comes from patterns from the underlying surface water, provided the Mackenzie River; annual local precipitation yields there is a conservativetracer of the pattern. The pat- only 0.10-0.13 m of fresh water [Bur•s, 1974]. Since tern of interest here is the spreadingof new inflow water ice forms with minor fractionation at the isotopic comas a thin plume directly beneath the ice; an appropriate position of the water from which it grows, landfast ice records the horizontal spreading with time of the isotraceris the oxygenisotopicratio expressed as •80. Oxygen isotope measurementstogether with salinity topically light plume over the isotopicallyheavier seaare used in polar oceansto distinguishbetween mete- water as illustrated schematicallyin Figures 2c and 2d. oric water, which includes river runoff and precipitaHere we useseasonaloxygenisotopeand salinity meation, and sea ice melt, the other major sourceof fresh surements in both ice and water to address two main water [Redfielda;zdFriedrna;,,1969; Ta;, a;zdStrai;,, themes. First, we distinguish between the processes

1980;Vetshteyn et al., 1974;(Sstlund andHut,1984]. that Meteoric water is isotopically light compared to seawater

and the ice formed

from

it.

For the Mackenzie

affect

the freshwater

balance

on the Mackenzie

shelf-river inflow and the formation and melting of sea ice. Second, we use the isotopic record in the ice at

MACDONALD

ET AL.:

FRESHWATER

BUDGET

OF MACKENZIE

RIVER

899 -15

30,000

Period of Core Growth

o,

20,000

.

"'-' " ' Winter Values (minimum flow) I :1 .,' \Summer Values -20

10,000

0

• 0

I 100



I 200

I

I 300

I

I 400

I

I 500

-25

Day(starting January 1, 1990) Figure 3. The hydrograph(solidline) for the MackenzieRiver from January1, 1990 to May

11, 1991(datasupplied by Environment Canada).The dashed lineshows the predicted 5180 compositionof the river below the confluencewith the Liard River as a function of time using KrouseandMackay's[1971]measurements, togetherwith flow data for the two riversabovethe confluence.

the end of winter to establish the pattern and rate of within the lead (station ISZ) and anotherwas taken river plume spreadingthat has occurredunder the ice from the ice on the northernedgeof the lead (station throughout winter. C-6). Ice coring, water samplingfor chemicaldeter-

minations,and conductivity-temperature-depth (CTD) profiling were carried out along sections that encom2.

Methods

Field Sampling During September 15-16, 1990, water sampleswere collectedby hydrocastsfrom the Canadian Coast Guard

Ship Henry Larsenalongtwo transectsacrossthe Mack-

enzieshelf (Figure la). At the time of samplingthe shelf was completely clear of ice as shown in Figure la. Sampling details and su•pporting oceanographic data

are givenby Macdonaldet al. [1991]. FromApril 26 to May 11, 1991,watersamplesand/or ice coreswere collectedat the locationsshownin Figure lb. In late February 1991 a lead 7-8 km wide opened intermittently about 60 km to the north of the Tuktoyaktuk Peninsula. By mid-April a continuouslead was open to the north of the peninsula up to the southwest side of Banks Island. At the time of samplingthe width

of thisleadexceeded 60 km in someplaces(Figurelb). We collectedmost of our samplesfrom the ice inshore of the lead, but one station was taken from a small boat

passedthe landfastice to beyondthe stamukhi(Figure lb). The general plan was to sample along designated transcots, but the actual site for each ice core was selected carefully from the air; we chose broad, flat regions of first year ice that were well away from ridges. The depth of snow was measured in the vicinity of the site, snow was cleared from the site, and a Sipre corer was used to obtain the ice in segmentsof about 70 cm in length. Water was collected from just under the ice through the hole made by the ice corer. Discontinuities in the ice were noted. The total length of the core was carefully reconciled with the depth of the hole, a procedure found to be particularly important for river ice which was brittle and tended to fracture easily during the coring process.The ice coreswere immediately cut with a saw into 10-era segmentswhich were placed in sealed plastic containers. The containers were returned to the laboratory in Tuktoyaktuk, where the contents were thawed, homogenized,and subsampledfor salinity and oxygen isotope determinations. Sampling details

900

MACDONALD

ET AL.:

FRESHWATER

BUDGET

OF MACKENZIE

RIVER

ters(ViennaStandardMeanOceanWater (V-SMOW), Universityof Washington,Seattle (-11.55% o), and ter column sampleswere obtained independently of the Antarcticstandards(-33.40%0)). Two in-housestanand ancillary measurementsare reported by Macdon-

ald et al. [1•1, 1•2] and Paton et al. [1•4]. Wa-

dards that were calibrated by the University of Wash-

ice-coringoperation by augeringa 25-cm hole and using either a pump or hydrowire/5-L Niskin bottle to collect water at various depths. Water was subsampled into glasssalinity bottles with care to prevent freezing during sampling or shipment back to the laboratory.

ingtonagainstV-SMOW (LTW3,-16.62% o and Tuk Snow, -27.17% o) were run daily. Equilibrationwas

carried out by stirring at 20øC for 15 hoursin batches of 12 sampleswith four standards. Moisture was frozen out of the gas usinga Pelltier coolerand the dried gas Analytical Methods fed to the inlet systemof a Nuchde ratio massspectromSalinity was analyzed at the Institute of Ocean Sci- eter for isotopic ratio determination. Errors were meaences(September1990 samples)and at the Depart- sured through replicate analysesand repeat measurement of Fisheries and Oceans laboratory in Tukyoy- ments of in-house standards. Results are expressedin standard 5180 notation 'with the V-SMOW

aktuk (May 1992) usinga GuildlincAutosal (model 8400A)instrument;data are reportedin practicalsalinity units [œewis aridPerkin, 1978].SampleswerestandardizedagainstStandardSeaWater (StandardSeawater Service,Institute of Oceanography, Wormley,U.K.).

the referencevalue. Agreement between replicate sam-

plesrun on the sameday was usuallywithin 4-0.1%o. Over the period of months required to analyze the sam-

ple suite, precisionwascloserto 4-0.2%o.

Repeat determinations •ndicate a precision of -t-0.003. Overall uncertainty based on bottle duplicateswas better than

standard as

3. Observations

4-0.02.

Oxygen isotopic compositionwas measuredby a me-

Water

Column

Data

thod similar to one developedby Whaite[1982].A de-

Salinity-5180dataforautumn(September 1990)and

tailed discussionof the procedure and error is presented

by Pator• et al. [1994]. Briefly, the analysiswas per- winter (April-May 1991)are summarized in Figure4. formed by equilibrating 5 mL of the water samplewith At salinities above 32 both the summer and winter data Matheson(instrumentgrade) CO2 of isotopiccompo- for salinityand/P80 covary.The samerelationshiphas sition -27.63% 0 establishedwith three referencewa- beenobservedin previousyears(September1986,Mac-

Sea Ice

•Oo

Late Summer 1990

o8 o

-5

[]

O

o

+

cO

+

Late Winter 1991

o RiverIce -15

RiverWater(-:•8.3) -2O

0

10

20

30

40

Salinity(psu) Figure 4. Salinityversus5180 for the watercolumndata collected in September 1990(solid squares),for the waterdata collectedin April-May 1991(crosses), andfor the ice cores(average of coresamplesis shown)collectedin April-May 1991(opencircles).

MACDONALD ET AL.' FRESHWATER BUDGET OF MACKENZIE RIVER

901

donaldet al. [1989];April, 1987andAugust1987,Mac- averagedabout 8 cm (as snow)which,again,is normal donaldand Carmack[1991]).At salinitieslessthan 32 for the regionat thistime of year lice Centre,1992].We the data pointsfall into two groups;the uppergroup found no apparent relationship between snowdepth and (Figure4, solidsquares) comprises the September 1990 amount of ice formed at a site. The 5x80 valuesshowcoherentpatternsboth within data,andthelowergroup(crosses) comprises the AprilMay 1991 data. The late September data have been eachicecoreandbetweenicecores(Figure6). Nearthe affected by the addition of sea ice melt. River water

top of a given core, which containsthe ice formed early

(bottomleft) hasalsoinfluenced manyof thesesamples in the season, 5180 tendsto be high(-5 to -1% 0). to a greater or lesser degree. The winter data extend The 5180 valuesusuallyincreaseslightlywith depth toward the bottom left cornershowinga stronginflu- near the top of the core, but for many coresthis trend enceof river water at somestations,but they are also is reversedabruptly at somedepth and 5180 decreases displacedto the right relative to the summerdata. This sharplywith depththereafter(for example,seecorePIdisplacementindicates the addition of salt to the water 2). The 5180 valuedropstowardbut not pasta limit column during sea ice production. of-15.7% 0 which is the value for ice growingfrom Ice

Core

Data

The ice cores collected from the zone of landfast

ice

pure Mackenzie River water. This sharp transition correspondedto depths where during sample collectionwe often noted the texture ½,fthe ice to changefrom plastic

were,with a few exceptions, of uniformlength (solid (seaice) to brittle and transparent(riverice). The ice coresrecord the sequenceof isotopicallylight portionsof the histogramsin Figure 5; rangeof 163 cm to 202 cm; average length 1.87-•- 0.11 m, number of water invading many but not all of the siteso For exsamplesn = 23). Theseice depthsare normalfor the ample,ice coresalongthe PI section(Figure 6) shift time of year; the Ice Centre[1992]givesaverageice to lighter5180 valuesat 20-30 cm in PI-1, 70-110cm depthsforApril 26 to May 11as203-210cm(n -- 4) at in PI-2, 120-150cm in PI-3, and 170-200 cm in PI-4;

KugmallitBay,181-182cm (n - 24) at SachsHarbour, on the seaward side of the stamukhi zone at PI-5 the and 180 (n - 27) at Cape Parry. Four of the coresin shiftto lighter5180doesnotoccurat all. Considering Figure 5 are much shorter than the others. Two of these that ice growsdownwar(!as winter progresses and that short coreswere collectedfrom beyondthe landfast ice MackenzieRiver wateris isotopicallylight comparedto zone(stationsCD-5 and AP-5), wherethe particular seawater, theicecoresections (PI, GI, CM, andAP) can patch of ice we sampledmay have startedgrowingin be seenas new winter inflow spreadingoutwardunder leads that opened late in the season. The other two the ice towardthe stamukhizone. Similarpatternsare shortcoreswerecollected fromthe river (April 25, HW- evidentin sectionsrunningparallel to the Tuktoyaktuk 13, 80 cm, and HW-14, 150 cm), wherethe conditions Peninsula (CM-1, L-2, AP-1, L-3, andCD-2). affectingicegrowth(e.g.,temperature,snowcover,and Salinitydata for our ice cores(not shown)tend to heatflow)clearlydifferfromthosein the estuaryandon have a high variance due to the nonconservativebethe shelf. For comparison,the norm for the Mackenzie haviorof salt duringand after freezing[Ma•lkut,1985]. Riverat Inuvikin April is 120q-27 cm (n - 29). Snow This behaviortends to make salinity a poor tracer for

cover(shaded portionsof the histograms in Figure5)

river ice that has formed beneath sea ice because brine

25O

200





"',

',

8

100

0

Location

Figure 5. Histograms showing the lengthsof the icecorescollected duringthe study(solid portion) andthedepthofsnow cover at eachsitewhenthecorewascollected (stippled portion).

902

MACDONALD

ET AL.'

FRESHWATER

BUDGET

OF MACKENZIE

ß ß

ß

ß ß

ß

8

ß ß ß ß ß ß ß

: ß

ß

RIVER

MACDONALD

ET AL.: FRESHWATER BUDGET OF MACKENZIE

draining downward from the sea ice will contaminate ice grown later. Under these circumstanceswe find that salinity can be used only as a very crude guide to dif-

chosen for the reasons

RIVER

outlined

below.

903 Sea ice melt

(SIM) hasa specialattributein that it can be positive (seaice is melted)or negative(seaice is formed).

œerentiate river ice from •ea ice [cf. Jeffrieset aL, 1988, Fractionation of' 5180 When Water Freezes 1989; Jeffries and Krouse, 1988]. T•ends within the cores,likethoseshownsoclearlyin the 5i80 data (FigUnder equilibrium conditions,ice forms with a small ure 6), are not observedin the salinitydata exceptin isotopicfractionationot' about 3% 0 heavier than the extreme caseswhere seawater has been completely re- liquid from which it was formed [O'Neil, 1968; Beck placed by fresh water at some time during ice growth andMuennich,1988;Le,•rnannand Siegenthaler, 1991]. (e.g., at PI-1 and PI-2). However,if the salinitiesand Our measurements of •lSO at the bottom of the ice 5180 valuesare averaged for the 15-20sections analyzed and in the water just under the ice showa fractionation within eachcore,then an increasing trend of 5180 with offsetof 2.66q-0.28%0 as the interceptof a very good salinityemerges(Figure4, opencircles). linearregression (Figure7). It shouldbe notedthat this observedfractionation applies to ice growing slowly at

4. Processes Affecting the Composition

the end of the seasonwhen it is 1.6-2 m thick. The slope

of Water

of the regression (0.999q-0.020) showsthat this offset doesnot dependon the 51so value nor, by inference,

and

Ice

on the salinity. Takingthe averagedifference for our

Water Types and Composition

paired data, we find ice to be heavier than the water

At any site in the Mackenzie estuary and adjacent

it growsfrom by 2.57 q- 0.10%0. This fractionationis

slightly higher than that determined by Melling and the followingthree primary types of water: a saline end Moore[19951(2.09%0) thr ice growingfrom water of a mem• (•o• m• •y• (P•L)), • •o• o• me- restrictedsalinity range of 29-31 on the Beaufort shelf. shelf the surface

waters

can be viewed

as a mixture

of

teoricfreshwater (MW) (essentiallyMackenzieRiver Polar mixed layer (PML), saline water. We runoffaugmentedby minor amountsof precipitation), have selected a seawater end member appropriate to anda source/sink of brackishwater(seaicemelt (SIM)

/formation)[Macdonald et al., 1989;Macdonald and the study period by referring to Figure 4. Below the

Carmack,1991; Melling and Moore, 1995]. To deter- point definedby 32.2, -2.5% 0, the 1990and 1991data

branch into two groups; above this the summer mine the compositionof a given sample with reference sets data narrow onto a well-defined curve which has been to these primary water types requires the measurement observedin previous years. We use this branch point of at leasttwo conservative properties(e.g.,salinityand 5180) and an appropriate assignment of properties to to define the underlying saline water mass commonto 1990 and April !991. This end-member repeachoftheprimary watertypes.6sflsnd andHut[1984] September resentssalinity of water at the depth to which winter outlined the principles of the calculationsas they apply mixing by brine-driven convectionand winter stirring to the Arctic Ocean am! gave massbalance equations has occurred. The properties of this saline end-member

to find the recipefor any samplewith salinity(S) and (PML) will differ from :yearto year dependingon the 5180(5)values water properties prior to freeze-up and the amount of

(1) ice grown over winter. (2) Mackenzie River meteoric water (MW). The (3) Mackenzie River has a well-established summer 5180 value of about -20.3% 0 [Krouseand Mackay, 1971; suggest where the subscriptsrefer to the three primary water Macdonaldet al., 1989],but our measurements types,F is the fractionof eachprimarywater type, and that a slightlyhigher51Sovalue (-18.3% o) prevails S and 5 refer to the salinity and 5180 compositions, in winter. Cyclic variationin the river's51so can be respectively, of the subscriptedwater types. The end- accountedfor by seasonalvariation in the flow of the FPMLq- FMWq-FSIM -- 1 FpMLSpML 4- FMWSMW4- FsIMSSIM= $ FpMLSPML q- FMWSMW q-FSIMSSIM= •

member compositions are listed in Table 2 and were _

isotopically lighter (21.3ø/•

Liard River. a major trib-

Table 2. Water Propertiesand ConversionFactorsfor Mass BalanceCalculations

WaterType

Salinity

MW (summer) MW (winter)

0.15 0.15

PML SIM River ice Fr actionation

32.2 6.2 0 -

0.07 0.07 0.1 0.4 0.1 -

--20.3 --18.3

-2.5 0.07 -15.7 2.57

0.66 0.09 0.1 0.1 0.06 0.12

Densityof ice is 900 kg m-a. Abbreviationsare s•, the standarderror of the mean; MW, meteoric water; PML, polar mixed water; and SIM, sea ice melt.

904

MACDONALD

ET AL.'

FRESHWATER

BUDGET

OF MACKENZIE

ß

ß

ß

RIVER

/

ß ß

//

/

ß ß // / / / / / / /

-5

/ / / / /

o

/ / / /

ß

-15

/ /

-20

' -20

i

I

-15



I

i

-10

I

i

-5

0

Water (•l•O

Figure ?. The •180 at the bottomof the ice coreplottedagainstthe •180 in the waterjust beneath the ice at the time of collection. The isotopic fractionation that occursduring freezing is evident as the offset of data from the diagonal which would be observedif there were no fractionation.

utary of the Mackenzie[cf. KrouseandMacka•t,197!]. (2.57%0). Similarly,for river ice the salinity is taken Above the confluence of the Liard the Mackenzie has to be 0 and 5180 to be -15.7% 0.

a 5180of-17.4%0. UsingKrouseandMackay's 5•80 valuesfor the two rivers abovetheir confluence,together

River ice melt (RIM).

The three massbalance

with flow data (Water Surveyof Canada),we predict equations allow us to solve only for the three primary

the 5180 belowthe confluenceto rise from -20.2%0 typesof water (PML, MW, and SIM). However,the efat peak flow in summerto a constant-18.2 -4-0.2%0 fect of a fourth water type, river ice me]t, must be considered. In previousstudiesthat have used5180 meathroughoutwinter(seedashedlineandright-handscale

in Figure3). Ourmeasurements onriverwaterandon surementsto distinguish between the Arctic Ocean's the two ice corestaken at HW-13 and HW-14 (Figure 6, 5180 - -15.7 4-0.06% 0, n- 23) confirmthat the river had a constant5180 value of about -18.3%o throughoutthe period of ice growth, providedwe allow for the 2.57%0 fractionation. The summerinflow of the Mackenzie River to the Beaufort Sea accounts for

freshwatersources;fractionationduring the freezingof meteoric water has been ignored. While this probably constitutes insignificant error for the Arctic Ocean interior, for the shelveswherethere is a substantialquantity of runoff, the effect of altering the compositionof freshwater by freezingsomeof it and then later melting

75% and the winter inflow, 25% of the total flow (Ta- it must be evaluated. For lack of a third conservative ble 1). For our massbalancecalculations wewill assume property we are left only with the option of determining that the summer5180 valueis appropriatefor Septem- what error the fractionationof MW into ice might produce. Using the end-memberassignmentsgiven in Taber 1990 and for stations outside the landfast zone in April-May 1991, where the winter inflow appears not ble 2 and applying the conservationequations,we find

to havepenetrated. The winter5180valuefortheriver that 1 kg of riverice (salinityof 0, 5180of-15.7% o)is

willbeusedwithinthelandfast zonefortheApril-May solvedalgebraicallyas an apparent mixture of 0.86 kg

of MW plus 0.17 kõ of SIM minus 0.03 kõ of PML. Therefore the conservationof properties requires that Sea ice melt ($IM). As seenin Figure4 (open the water out of which ] kg of river ice grew must, accircles),ice from the Mackenzieestuarydoesnot have cording to the algebra, contain an apparent or virtual uniform 5180 OI salinity.Here we treat ice as a two- extra 0.14 kg of MW and -0.17 kg of SIM. The algecomponentmixture of seaice and river ice. The salinity braic solutionbalancessalt by removinga smallamount of first-year sea ice is estimated as 6.2 from the six ice (0.03 kg) of PML. In terms of properties,the water corestaken from the seawardsideof the stamukhizone. would be slightly lighter isotopicallybut would witness The5180ofseaice,+0.07%0,isto bethe5180ofthe no changein salt content. If this ice is subsequently salineend-member (PML) plusthe fractionationfactor melted and mixed back into the water, then the net 1991 stations.

MACDONALD

ET AL.'

FRESHWATER

BUDGET

OF MACKENZIE

RIVER

905

position returns to 1 kg of MW. Provided we can independently estimate how much meteoric ice has been formed in the estuary, we can apply an appropriate correction to the algebraicsolutions.

estimated from the water column data should, in principle, provide independentaccountingof salt and fresh

5. The

and

Distribution

water

for the outer shelf.

Distributions Sea Ice

of Freshwater

Components in Ice and Water Ice

and

Water

Before discussingthe water column budgets, we must consider

River

of Meteoric Melt

Sea Ice

the

error

introduced

into

the

mass

balance

1000kg m-a). From Table 3, whichsummarizes the

equations by isotopic fractionation during the freezing of river water. We have estimated that an average of 0.72 m of MW was captured into the growing landfast ice; therefore consideringthe effect of fractionation during freezingas disscussedabove, this will produce an apparent 0.1 m of MW and an apparent -0.12 m of SIM in the water beneath the ice. That is• the equivalent heights of MW and SIM calculated from the mass balance equations will, on average, overestimate MW and underestimate SIM by these amounts, or, sincenegative values for SIM imply the formation of sea ice, the mass balance calculations will overestimate sea ice production by this amount. T},•.eseerrors are small relative to

equivalentheightof seaiceis 1.00m (Table4).

in Table 4).

An estimateof the equivalentheight (liquid) of river water

R and

seawater

H-

R that

has frozen

into

a

sample of ice from any location can be calculated from

61sO,the average61sOvaluewithin the ice core;

(61sO- 2.57- 61SOs•M)

(4) R-- •i-•-M•• •i•O-•I•IM) Hpice/Pwater where H is the measured depth of ice and Pice and

Pwater are,respectively, the densityof ice (900kg m-a) and water (taken here to be the freshwatervalue of

integrated componentsfor each station, we find the av- the estimatedequivalentheightsand their ranges(Taerage equivalent height of river ice for all ice corescol- ble 3). We haveusedthe estimatederror only to correct lected from the landfast ice is 0.72 m and the average the averagevaluesfor landfastice zonein winter (given The equivalent height calculations for river ice and

September 1990. 'The late summer water properties, shownschematicallyin Figures 2b and 4, are shown the relatively flat and stable landfast zone. However, in section in Figures 8a-Sd as isopleths of fractional beyond the landfast zone, ice can advect out of the re- composition. Autumn storms have smoothed out plume gion during winter as evidencedby the large polynya structure, and inflow from the Mackenzie has declined. sea ice at each station

are reasonable

estimates

within

which openedin late February 1991 (Figure lb) [see Sections of MW at this time show the river water to be Melling and Riedel, 1993]. Conversely,advectioncan distributed relatively uniformly through the water colbring offshoreice onto the shelves,and, particularly for umn among the stations. Only near the surface at the

sites that open later in the season,the observedthickness of new ice will not be a reliable

estimate

inner

stations

of the eastern

section

do we see a struc-

of the to-

ture identifiableas the MackenzieRiver plume (Figure tal amountof ice actuallygrownat that site (e.go,AP-5 8a, stations 21-22; also see Figure 4 where these two and CM-5). The ice on the outershelftendsto be more pointsare displacedtoward the bottom left). There is

heavily ridged which means that measurement of ice at flat locationsgood for landing aircraft will considerably underestimate the net ice production for the shelf even if the fraction of ice coveris correctly estimated. Using height and frequencydata for ice ridges, Prinsenberg

minor variation in the equivalent height of MW within each section, shown as bars at the top of Figure 8, but the westernmostsectioncontainssignificantlylarger

tributed an additional 25 cm of ice in Hudson Bay and as much as 58 cm of ice in Foxe Basin where extremely roughice is found; thesefigurestranslate into a net con-

In contrast to MW, SIM is distributed predominantly in the upper 10-15 m and it tends to be more concen-

quantitiesof MW (2.93 m versus2.37 m; probability P < 0.05). The averageequivalentheight of MW for [1988]estimatedthat ice accumulatedinto ridgescon- all of the Septemberstationsis 2.57 4-0.12 m (r, - 20).

tribution

of about 30% over what would be estimated

trated towardthe outer shelf(Figures8c and 8d). At

from ice thickness measurements alone. Melling and

this are negative valuesindicating a lens of water that Riedel[1993]suggestthat 30% is probablyan underes- contains remnant brine from sea ice production during

timate since the under-ice keels are much larger than the above ice sails used in the calculations.

The advan-

tage therefore of using the mass balance equations to integrate the water-type content for the water column is that we can determine the image of ice production

left behind as salt in the water (negativeSIM), that is, water will capture the imprint of sea ice produced whether

or not the ice is later

advected

off the shelf or

piled into ridges. Comparing the measurementsgiven in Table 3 for sea ice thickness and sea ice production

the previous winter. •his deeper• negative SIM tends to o•set •he positive S•M in the surface layer with the result that the equiw]ent heights c•culated for the top 40 m are less than whnt is contained just in the surface 10 m. •he average equivalent height for SIM in

Septemberis 0.82• 0.11m (•-

20).

April-May 1991. The distribution of MW under the ice at the end of the followingwinter is shownfor the

stationsin Figure lb as four sections(Figures9a-9d).

906

MACDONALD

ET AL.'

FRESHWATER

BUDGET

OF MACKENZIE

RIVER

Table 3. Equivalent Heights in Meters of Freshwater Components Water

Station

Latitude, øN

Longitude, øW

Column

Ice

MW, top 40 m

SIM, top 40 m

Rice

021

69 51.11

133 19.17

2.61

0.14

-

022

69 55.28

133 23•39

2•55

0.04

-

023 024

70 01.63 70 08.30

133 24.43 133 26.03

1.71 1.76

0.18 0.89

-

025 026

70 13.12

133 33.73

70 18.20

133 40.20

2.11 2.00

0.21 0.78

-

Sice

Tice

-

027

70 23.37

133 46.14

2.23

0.31

-

028

70 29.10

133 53.14

2.62

0.51

-

029 030

70 33.81 70 39.21

133 57.78 133 59.47

2.37 2.56

0.91 0.77

-

032 033

70 49.29 70 55.59

134 17.15 134 23.68

2.59 2.95

0.28 1.32

-

034 035

70 59.09 69 08.00

134 29.78 137 41.62

2.74 1.97

1.30 0.83

-

036

69 16.90

137 51.02

3.16

1.47

-

038

69 24.65

137 58.60

2.54

1.11

-

042 049

69 30.85 69 37.95

138 04.47 138 14.09

2.66 2.99

1.41 0.97

-

050

70 24.89

139 03.39

3.58

0.72

052

70 09.74

138 48.34

3.63

1.94

-

HI-1 HI-2 KP-2 GI-1 GI-2 PI-1 PI-2 PI-3 PI-4 PI-5 CM-1 CM-2 CM-3 CM-4 CM-5 ISZ C-6 AP-1 AP-2 AP-3 AP-4 AP-5 L-1

69 69 69 69 69 69 69 69 69 70 69 69 69 70 70 70 70 69 70 70 70 70 69

28.00 36.82 23.94 30.24 38.05 43.84 48.75 54.54 57.63 02.28 45.01 52.05 56.15 00.99 08.06 11.60 28.29 59.74 04.21 08.52 15.44 18.05 53.00

138 138 138 136 136 134 134 134 134 135 133 133 133 133 133 133 133 131 131 131 131 131 133

50.58 47.84 07.47 09.49 21.95 31.26 38.23 48.54 52.94 01.50

1.96 6.11 5.38 4.21 1.91

L-2 L-3 CD-2 CD-3 CD-4 CD-5

69 70 70 70 70 70

54.03 16.06 25.01 31.46 39.95 48.32

132 130 129 129 129 129

18.00 21.32 21.57 24.05 25.23 39.73 51.47 27.97 25.20 25.16 29.24 21.28 59.98 13.88 31.69 39.98 39.47 40.83 40.35

4.68 5.49 2.64 1.70 1.42 2.12 1.81 5.51 5.05 2.16 0.83 1.02 1.31 4.52 2.83 -

-0.05 -0.40 -1.36 -1.92 -2.15 -0.53 -1.26 -2.11 -2.09 -1.85 -2.52 -1.88 -0.60 -1.35 -2.17 -1.25 -1.43 -1.86

0.36 0.24 0.36 1.71 1.07 1.58 1.04 0.71 0.39 0.18 1.14 0.97 0.27 0.18 0.74 0.40 0.25 0.13 0.14 -

-0.94 -0.51 -

0.74 0.38 0.18 0.11 0.13 0.06

1.33

1.69

1.46

1.70

1.33

1.69

0.01

1.72

0.59

1.65

0.18

1.76

0.66

1.70

1.09

1.80

1.41

1.80

1.58

1.76

0.53

1.67

O.7O

1.67

1.22

1.49

1.56

1.74

-

1.71

-

0.99

0.82

1.56

1.11

1.51

1.5

1.75

1.47

1.6

1.24

1.38

-

1.80

1.08

1.82

1.30

1.68

1.47

1.65

1.36

1.47

1.64

1.77

0.82

0.88

Water column is the top 40 m. Abbreviationsare MW, meteoric water; SIM, sea ice melt; P•ice,river ice; Sice,seaice; Tice, total ice.

MW is contained predominantly withina 4- to 5-mdeep spreadingout to the stamukhi will be further discussed layer of fresh water under the ice at nearshore stations.

below. Despite the apparent containment of new winter

For each section there is a front across which stratifi-

inflow inside this stamukhi

cationall but vanishes(see,e.g., PI-4 to PI-5 or AP-3

tains a substantial quantity of MW at the end of winter

zone the offshore water con-

to AP-4). Thisboundaryof winterinflowspreading in (Table 3, averageequivalentheightof MW is 1.54m). late winter generally follows the stamukhi zone at the

Mixing by autumn storms,followedby brine production

edgeof the landfastice; the rate and extentof plume

in winter

at these outer stations

redistributes

the MW

MACDONALD

ET AL.:

FRESHWATER

BUDGET

OF MACKENZIE

RIVER

907

Table 4. Average Equivalent Heights of Meteoric Water, Sea Ice Melt, Sea Ice, and Meteoric

Ice for the Landfast

Zone and the Offshore

Landfast Region

m

Offshore Region

Period

WatexType

s•n)

m

s•n)

Sept. 1990

meteoricwater seaice melt

2.$1 0.73

0.14(10) 0.16(10)

2.82 0.90

May 1991

meteoricwater meteoricice seaice melt seaice

4.16 0.72 --0.94 1.00

0.39(13) 0.11(17) 0.19(13) 0.11(13)

1.54 0.15 --1.88 1.47

0.15(10) 0.16(10) 0.18(7) 0.02(8) 0.16(7) 0.05(7)

Abbreviations are m, mean; s•, standard error of the mean; n, the number of data points.

uniformly throughoutthe top 40 m where it contributes

On the cross-shelf sections(Figures9e-g) the mini3-5% of the water (Figures9a-9d). Furthermore, wa- mum (most negative)valuesof SIM tend to be found ter displaced from the landfast zone by winter inflow nearthe outsideedgeof the landfastzone(PI-4, CM-3, or exchangedbelow the depth of the ice keels has con- and AP-3) rather than in the polynyabeyondit where tributed

additional

MW

to the offshore zone.

Within

enhanced ice production occurs. The location of this

the landfast ice zone the equivalentheightsof MW are maximum is probably due to the capture of the rejected muchhigherthan they are offshore(Table 3 and solid brine in the buoyantspreadingplume itself, which then barsat the top of Figures9a-d; the averageis 4.16m). prevents the water column from conveeringto the botThe equivalent height of MW at some of the stations tom and thereby mixing the brine through a larger volnearest to shore is simply limited by depth, which ex- ume. This "snowplow"effect of the spreadingplume is plains why a maximum is reachedat the central stations more than of academic interest as it would undoubtedly (5-6 m at PI-2, PI-3, and CM-2). be an important processin the transport of contamiThe

total

amount

of winter

inflow

that

has frozen

nants released beneath

t.he landfast

ice. The nearshore

into the ice at eachsite (calculatedaccordingto (4)) stations, which obviously produced brine early in the decreases as oneproceedsawayfrom the river (Table 3 season,have apparently 'beenflushedby the winter river and cross-hatchedportion of the bars depicting equiva- inflow, and it is probably this brine that tends to aug-

lentwaterheightoficein Figures9a-d). At the stations ment the middle stations. The histogramsat the tops of beyond the stamukhi, MW makes only a small contri- Figures9e-h showthat the amountof seaice (hatched bution to the ice formed(0.06-0.19m). The MW in- part of the icehistograms) increases with distanceaway corporated into the offshoreice comesmost likely from runoff left over from the previous summer rather than

from the river mouth and that the equivalentheight of negative SIM contained in the water beneath the ice

from new river inflow in winter.

(solidbars)similarlytendsto increase.Within thelandIn contrastto the situationin autumn(Figures8cand fasticezonethe averageseaiceproduction(1.0 m) is al8d), all valuesfor the fractionof SIM in the water col- most identical to the equivalentheight of negative SIM umn are negativein late winter (Figures9e-h), clearly (-1.06 m) contained in the waterbelow.However,the

indicating that salt has been added. The shift from summer to late winter conditionscan be seen clearly in

similarity of these two numbers does not imply that the salt rejected by the sea ice exactly balances the brine the salinity-51soplot (Figure4); in summerall points added to the water column during the current winter; lie to the left of the MW-PML mixing line, indicating becausethere was an average+0.73-m equivalentheight the addition of SIM. In late winter all points are found of SIM in September, there must have been either a loss

salt which is seenin/he water columnas negativevalues or an inflow of negative.SIM from outer parts of the

for SIM. Near the branchon Figure 4 (32.2, -2.5%0) thereare a numberof points(crosses) displacedfar to the right. These points, arhichindicate that large quantities of brine have been injected, come from stations at

shelf(or a combinationof both processes) to balance the salt budget. Winter

Inflow

From

the

Mackenzie

River

the entranceto LiverpoolBay (Figure 9h), a location We calculated the total winter inflow of the Mackenshown to be particularly effective at enhancing brine in the water due to its shallow depth and semienclosed zie River from freeze-up(October16, 1990)until the

shape(E.G. Carmacket al., Water and ice relatedphe- day we collectedour last ice core (May 10, 1991) to

of Canada are nomenain the coastalregion of the Beaufort Sea: Some be73km3 (datafromtheWaterSurvey parallelsbetweennative experienceand westernscience, plottedin Figure3). To determinethe net amountof river water delivered to the estuary, we must subtract submittedto Arctic,1994).

908

MACDONALD

ET AL.' FRESHWATER

BUDGET

OF MACKENZIE

RIVER

33 34

II 2o

ß3o

1.

•o 10-

30. 0.050 4O

-4

-3

ß-• 'r' G) •

2!

0

.-

22

2:3

-

"

24 ß ß

21627

25 -,

ß

ß

2__829 ß 3o . ß

ß

ß

:32 i

33 34 I I ß

.2 1

ß

0

ß

lO

10



ß

.s:: 20

.. ß

30

..

4O

.•

2,

'r'

1-

• •o"

:

' --•.•-0.025

---

36

• K

10 \

:

ß ß L0.025 •••........••

I /

!

38

,

,

,

42

I

I

0.10•

50

49 I :

52

ß

0075•0.10_ '•

•oo•,,.o,o•

I :



40

0

50

1O0

150

Distance (kin)

Figure 8. Data for September 1990stations(Figurela) showingcontours of fractionalcompositionfor meteoricwater(MW) at (a) stations21-34and(b) stations35-52andfor seaicemelt (SIM) at (c) stations21-34and(d) stations35-52. Solidhistograms at the topsof the diagrams give the equivalentheight in metersof the freshwatertypes.

MACDONALD

ET AL.' FRESHWATER

7'

ß •'

Pl-2

PI-3

PI-4

5

Pl-i

PI-.S

ISZ

1 o

909

6

•' 5

:•

RIVER

7

6'

._o•4 3• 3

BUDGET OF MACKENZIE

4 3 1

o

o

o

0.1'

lO

ß

'

0.0S

10

. ß

2Q

40

10

0

20

30

40

50

CM-2 CM-I

CM-3 CM-4

ISZ

CM-S

C-6

0.1

2O

•. 2o-

ß

40

E•, 6

AP-I

50

55

40

80

90

AP-I

AP-2

CM-1

L.-2 L-3

AP-3

• •

AP-4

21

CD-2

AP-5

LB-I

CB-4 CBTb

0

10.

lO

10-

lO

20.

20

20

20

30-

-30

1•0

2'0

3•

30

.30

,40

40

0

50

100

150

200

2•0

Distance (km)

Figure 9. April-May 1991stations(Figurelb) showing contoursof fractionalMW composition at (a) •ectionPI, (b) 8eelionG•, (½)•ectionAP, and(d) section alongtheTuktoyaktuk ?eninsula andSIM at (e) sectionPI, (f) sectionCM, (g) sectionAP, and(It) sectionalongthe Tuktoyaktuk Peninsula. I-Iistog•amsor the equivalentliquid height or the variousr•,hwa•e• component•

dxown;•oliclba• •ere• •o integ•ateclwate• •olumn values,o•en/hatdteclba• •ere• to ice. t•atd•eaportionor the ice histogram•ere, to the •ame r•,hwate• componentgivenin the •olia •a• ro• eachaiag•am(i.•., MW ro• 0a-0a aria SI• ro• 0e-0h). eo• SI•, neeativewate•column e•uivalent height• a•e sitownas black ba•s below the line. Dashea hi•tog•am• a•e given

locations wi•ereonlythecle•thoricewa•rnea,•ea. Openwater(0aand0•; 0eariaOr)i• •hown

a•a •i•lea •u•race neaxly levelwititthetopor•the ice(iceeee•oa•a w. usually about•0

i 40

910

MACDONALD

ET AL.' FRESHWATER

BUDGET

OF MACKENZIE

RIVER

_3'1 PI-1PI-2.A-.B•/z,.,.••L' PI-3 PI-4 PI-S ISZ 13 ':•..:•':.•....,::........::::•:...•.'-,----..•.•,,'"'• I \ 'ß o 40' "'•:.! :'(":'"'":' ":"•l ':":'/" 0

10

20

30

, 50

40

40

,3O

ß40 90

0

10

L-2

CM-1



20

AP-1

:30

L-3

CD-2

40

LB-1

CB-4 CB7b

:

-0 075

--0.125

'

.:.:.•3• .......iiiiiiiii........... •:.....•...:•• •oo '7111

10

••'øø

2O

4O



5o

1oo

150

200

Distance (km)

Figure 9. (Continued)

2,•o

MACDONALD

ET AL.:

FRESHWATER

BUDGET

OF MACKENZIE

RIVER

91!

that the ice corescollectedat HW-13 and HW-14 (in the MackenzieRiver) and the very shortcoretaken at

the portion of this inflow that has been captured into the growinglandfast ice. The calculated averageequiv-

alentheightof riverice (0.72m) spreadoverthe region CM-5 will be treated separately. For all coreswe assign

withinthelandfast ice(16,000 km2) accounts forabout 12 km3 or about 15% of the winter inflow. The amount of Mackenzie River left to spread under the ice in the

the collection

date to the bottom

of the core.

Ice growth models have been reviewed by Maykut

water column from October 16, 1990 to May 10, 1991 is therefore about 61 km3.

[1985,1986],Orastedtand Wettlaufer[1992],and Leppdfanta [1993]. For simplicitywe will adopt an empirical modeldevisedby Anderson[1961]. For youngice,

6. Plume Spreading' The Record

Anderson found the relationship between the thickness of the ice H, and the cumulative freezing-degreedays

in the

(0 - fj(T! - T•) dt where T! isthefreezing pointof

Ice

water,Ta is the air temperature,andt is time) to be

Three stepsare requi,ed to use 5•80 profilesin ice cores to infer

the

estuarine

circulation

under

+

the ice.

The first and most important step is to convert depth in the ice core at eachsite into "time" by estimating the growth rate of the ice. The secondstep is to convert the

Accordingto Maykut [1986], good estimatesare obtained with Anderson's[1961]model without specific knowledge of winds, cloudiness,snowfall, or oceanic

5•sO valuein the coreto the corresponding 5•sO value heat flux. The cumulative freezing-degreeday record in the water from which it originated;this is simply the from Tuktoyaktukand Paulatuk(on the southerncoast fractionation correction shown in Figure 7. The third of AmundsenGulf about 300 km to the east) are alstep is to infer the water's salinity from the calculated mostidenticalfor the reriod in question(Figure 10). water 51sO.In this way,the 5•sO ice coredata can be Freezing-degreedays started to accumulateat the end usedto produce a seriesof maps showingthe evolution of September 1990, but presumably during the first 2 of surfacesalinity under the ice throughout the winter. weeksof October the seasurfacewas being cooledto the freezing point. Using the freezing-degreeday record at Assigning a Time Scale to the Cores

Tuktoyaktukstartingon October16, we find from (5) Weeklyice charts(EnvironmentCanada)from Oc- that by aboutthe endof December(Figure11) the first meter of ice had grown. During this time period only about 10 cm of snow had accumulated at Tuktoyaktuk.

tober 2, 1990, to November 27, 1990, show that new landfast

ice started

to form after October

9 when most

of the nearshorewater was open and before October 16 when most of the region was coveredby new ice. Comparison of the October 9 and October 16 charts shows that the landfast ice from which we gathered most of our cores became coveredby new ice almost synchronouslyø

The Anderson[1961]modelappliedfor the entirewinter predicts the ice thicknessto grow to a maximum of

1.70m (Figure 11, bottomline) whichis lessthan the

average 1.87 m we observedfor our ice cores. The inclusion in the model of a snow cover, which insulates We thereforeassignOctober16, 1990(4-1 week),asthe the ice, would tend to decreasethe predicted ice depth date at the top of all of the ice cores. Note, however, at the end of winter and thus cannot bring the pre5,000

4,000

3,000

Tuktoyaktuk 2,000

aulatuk 1,000



0 0

•'•

I 50

I

I 100

I

I

I

150

I 200

i

I 250

Day (startingOctober16, 1990)

Figure 10. Cumulatiw•.freezing-degree days starting October 16, 1990, and finishingMay 10,

1991(the last daya corewascollected),at Tuktoyaktuk(dashedline), and raulatuk (solidline).

912

MACDONALD

ET AL.: FRESHWATER

BUDGET

OF MACKENZIE

RIVER

Kugmallit Bay

2!f-

•,

!

Sachs Harbour

o ø•,•'"•

/

.. •Anderson

o

o

o•

0.5• • • 0

50

100

150

200

Day(startingOctober16, 1990) Figure 11. Ice thicknessas a function of time. The bottom line is producedfrom Anderson's [1961]formula usingthe freezing-degree day recordat Tuktoyaktukstarting October 16, 1990

(data suppliedby EnvironmentCanada).The opencirclesshowthe averagethickness of iceasa functionof time for KugmalhtBay duringwinter[/ce Centre,1992].The middle(SachsHarbour)

andtop (KugmaHitBay)linesshowthefit to "average" icecorescollected onday200(May 2) at the respectivesites using the method describedin the text. The first meter of the averagecore is assumedto conform t• the Anderson model with the 1990-1991 keezing-degreeday record at Tuktoyaktuk.

dieted

value closer to the observations.

Therefore

for

the deepersectionof the ice we have proceededby sim-

at Inuvik lice Centre,1992]and thereforemustbe considered

normal

for their

location.

From

the Ice Cen-

ply assumingthat after the first meterof growth,the tre data we could estimate depth as a function of time ice thicknessincreaseda• the squareroot of time. Each

for each river core simply by assumingthe minimum corewastreated separatelyaccordingto its total length growth scenario fo• the short co•e and the maximum ...

and th e' date it was collected by fitting an equation of

growthscenario for the longcore.Sincethe •zso signal

the form H - a + bt•/2 where a and b are fitted con- does not vary significantly within either of these two stants and t is time.

We have evaluated

the uncer-

cores(Figure6), thereis little benefitto be gainedby

tainty in this method of dating the coresby applying assigninga time frame except to note that the long core

it to the Ice Centre[1.c•92] averageice growthstatis- andprobablythe shortcorespanthe winterstartingin tics for KugmalhtBay (4 yearsof observations, Figure mid-October and therefore fairly representthe Macken11, opencircles)and SachsHarbour(28 yearsof ob- zie River •So composition throughoutthe periodof ice

servations, Figure11,closedcircles).We havefit these cover. For the anomMouslyshort core at CM-5 the Anaveragedata as if they werean actual corecollectedby derson model taken, together with the freezing-degree usonApr• 26;weused(5) with the 1990-1991freezing- day record(workingbackwardkom the time of collecdegreeday recordand assumedthat the corestarted to tion), suggests that thisicestartedto formaboutFebrugrow on October 16. The fit is shownas the middle ary 15, 1991. The •zso signalin the coreis consistent line for Sachs Harbour and the top line for Kugmalht

Bay in Figure 11. Kugmalht Bay givesthe worstfit in

with this picture; nevertheless, we do not assume any rigoroustime scalefor this core.

the regionwherethe ice thickness jumpsbetweenday 1!0 and day 120. Here our modellags the average Estimating Water Salinity from •180 CC

•,

data by about 2 weeks. For SachsHarbour, agreement

is usuallywithinabout5 days.Therefore weestimate Water salinity at the time of freezing is calculated

the uncertainty involvedin the datingprocess for our from the salinity-•lso correlationdiagram. However, ice coresto be -4-1weekat the beginning,perhaps-4-2 becausethe summer and winter correlationsdiffer, some adjustmentbetweenthe seasons is required. One month before freeze-up the water on the middle and outer The two river cores(HW-13 and HW-14) bracketthe shelveshad the •lSO-sa.linitypropertiesshownby the rangeof ice depthsreportedfor the Mackenzie River solidsquaresin Figure 4. Only two of the sampleswere weekstoward'the middle of the core, and within a few days toward the bottom of the core.

MACDONALD

ET AL.:

FRESHWATER

BUDGET

s•tonglyinfluencedby •he MackenzieRivet (surfacea• s•a•ions021 and 022). The MackenzieRivet•ssummer composition (-20.3) is wellknownfrom previousmeasurements in the riverandestuary[KrouseandMackarl,

OF MACKENZIE

RIVER

913

of changeof salinity, as evidencedby the slope of the salinity as it shifts from the top trend line toward the

bottom trend line, is surprisinglyconsistentamong the

cores(0.35-0.45d-•).

1971; Macdonald et al., 1989; Macdonald and Carmack,

The timing of river water invasion varies with dis1991]. Thereforewe are confidentthat the landfastice tance from the source; Figures 13a and 13b show the started to grow from surface water that was a mixture

time distribution

of the shelf surface water: and Mackenzie

trol buoyancyat a particularsite (i.e., salinitypasses the maximum)and the time at whichthe river plume

River water as

approximated by the to.p line in Figure 4. At the end of winter there is a closerelationshipbetween5•80 and

when river water

first starts

to con-

"front" as definedby the 15 isohalinearrives at a partic-

salinity (bottom curveon Figure 4), but it is different ular location. Once river water starts to control buoyfrom

the one established

in late

summer.

To convert

ancy at a site, it essentially reversesthe density en-

the 51sOwater valueto salinity,we haveusedthe two hancementtrend driven by brine rejectedfrom growing relationshipsshown in Figure 4, starting with the top

sea ice.

oneandmovinglaterallyto thebottomoneusinga t 1/2

Figures 13a and 13b show that the plume spreads faster along the coast of the Tuktyoyaktuk Peninsula

scale to mimic the rate of ice formation.

We estimate

the error in salinity determined this way to be about (1.3 cms-•) than it doesperpendicular to the coast +2.5 for the saltiestsurt;ace water early in the growing (0.2 cms-1). Coresat HI-1 and KP-2 yield a r•te of season but less than this for freshet water and for the spreadingalong the west coast of MackenzieBay that late winter

values.

is intermediate

to the •bove

two values.

Toward

the

end of winter, the under-ice river plume, as defined by

Spreading Rate of the Mackenzie Plume

the 15 isohaline,coversan area of about 16,000km2. Surface water salinities

Combining the growth rate of the ice with the cal-

measured

at the time the cores

werecollectedshowa gradientin salinity(at a salinity

culatedsalinityyieldsthe water'ssurfacesalinityas a

functionof time at eachice corelocation(Figure12). of 15)of about0.4km-• alongthecoastand1.9km-•

(perpendicularto the coast). Multiplying the spread-

The patterns are remarkably clear. Most coresshow a

salinitymaximumin time. They start to growin water ing rate by the respectivesalinity gradientwith distance with a surfacesalinity of 20-30. The steady increase suggeststhat the rate of changeof salinity at the various of surface salinity due to brine release and convection sitesshouldbe about 0.3-0.45 d-•; thesevaluesmatch of the water can be se,min the top cluster of trend the rates of change of salinity at the various sites in lines;for the outerlocations(PI-5, AP-5, AP-4, CD-3, Figure 12 calculated solely from the recordsin the ice CM-4, and HI-2) surfacesalinity eventuallyincreases cores. The thick dashedline in Figures 13a and to 30-32.5. At the inst•ore stations the surface salinshowsthe boundary of influence of winter inflow from ities depart from the top trend line and drop toward the Mackenzie River at the end of winter; the salinity zero as river water invadesa particular site. The rate at stations beyond this boundary was either constant

AP-4 AP-S CO-3

3O

CM-4 HI-2

AP-3

Pl-4

I

0

I

50

I

I

100



I

150



I

200

Day(startingOctober16, 1990)

Figure 12. Surfacewater salinity shownas a functionof time throughoutwinter at the ice core sites.Valuesweredeter]ninedusingthe 5•sO recordin the ice coresas describedin the text.

914

MACDONALD

ET AL.: FRESHWATER

BUDGET

OF MACKENZIE

RIVER

i

(.) 71øN

200 m

50 m

ß

70 ø

6g ø

140øW

130 ø

135 ø

(b) 71øN

200 m

70 ø

140øW 135 ø

130 ø

Figure 13. (a) Isochv•nsfor the day at whichsalinityreachedits maximumand beganto decrease within the landfastice and (b) isoehrons for the day at whichthe salinitydroppedto 15 within the landfast ice. The thick dashed line showsthe boundary beyond which salinities continued to increase up to the date the core was collected.

or continuingto increasewith time at the date of sam- made it to the outer siteIf through displacement, expling. The thick dashedline coincidesroughly with the changeof deeper water below the level of the ice keels, nearshore

side of the stamukhi

zone found

on the south

and leakagethroughthe ice boundary[Ornstedtet al.,

side of the polynya shown in Figure lb. The heavily 1994]. The advancingfront of the freshenedwater and its efridged ice along the southern boundary of the polynya is apparently an effective backstop for the winter inflow fect on surfacesalinity with time throughoutthe winter seasoncan also be visualizedby contouringsalinity on less• some of the water in the nearshore has probably a diagramof distanceversustime (Figures14a-14d).

MACDONALD

ET AL.: FRESHWATER

BUDGET

OF MACKENZIE

RIVER

915

Thesediagramsshowa ridge of high salinity (dashed line) which signifiesthe divide wherebysalt enhancementfromseaiceformation(belowthe dashedline) be-

2OO

comesoutstripped by freshwater supply from the river

1•0

inflow(abovethe dashecl line). If we visualizepassage of time along a sectionas the z axis movingupward at a constant rate in the various panels of Figure 14, we see the winter inflow progressivelyshutting off convection as the line crossesthe ridge at each of the locations. The advancingfreshwaterplume essentiallycaptures any brine rejected from growing ice, thus cutting off convectionto the bottom. As shown earlier in Figures 9e-9h, this advancing front acts as a "plow" to carry with it the brine or other substancesintroduced

100 •0

200

CM-1

CM-2

CM-3

I



I

iso (b) /

CM-4

at the inner

/ •

stations.

7. Budgets and Fluxes for MW, and

100

SIM,

Ice

The data presentedin Figures 8 and 9 and the average equivalent heights categorized by area in Table 4 may be used to summarizeregional budgetsfor MW

50

and SIM

for the shelf and to infer the fluxes that

must

have occurred between the two seasons. Considering

{o

the distributionof winterinflow(Figures13aand 13b), we divide the shelfinto two regions,the nearshorezone

AP-1

AP-2

AP-3

AP-4

(16,000kin2),whichincorporates theareawhereland-

AP-5

fast ice growsand into which the Mackenzie River continues to flow and spread during winter, and the outer

2oo] 1•0 t ' ' / ' /[•

• •oo t t o

shelfdeeper thanabout20m (48,000 kin2),whichdoes not receive much winter

inflow.

The choice of this divi-

sionis supportedby eartier work [Macdonaldand Cavmacl•,1991],and it is clearthat the processes occurring under

the landfast

ice out to its stamukhi

•.one bound-

ary are different than those on the outer shelf in winter. Meteoric

{o

Water

Figure 15a showsthe dispositionof MW betweenMay and September. The total volumesfor the sampling L-3

CD-2

'

'

CD-3

periodare shownin boldnumbers(in cubickilometers), while fluxes betweencompartmentsare given in italics within

the arrows.

The fluxes refer to volumes of the

watertype (cubickilometers)that haveenteredor left the compartmentsbetweenthe samplingperiods(i.e., fluxesin the late summer panel refer to the period from

late winterto late summ,•r).For example,weknowthat iv .LL •,,,

AS

.L.L L A.L

a

ß

During winter, 73 km of river water wasadded, 12 km

of river ice wasformed,and 67 km3 remainedin the water o

40

80

120

160

column

at the

end oœ winter.

To balance

this

budgetrequires thelossof31km3 fromtheinnerregion

during winter. Figure 15a shows that the nearshore tends to hold greater quantities of MW in winter than in summer, and the reverseis true for the outer shelf. Figure 14. Contours of salinity with distance versus of MW goinginto ice production accountfor timefor sections (a) PI, (b) CM, (c) AP, and(d) thesec- Amounts 15% of the inflow. In the nearshore the winter inflow tion running along the 3•ktoyaktuk coast. The dashed exactlybalanced by the estimated line showsthe top of the salinity ridge indicating the di- (73km3) is almost vide between salinity increasingwith time and salinity volume of MW in the water plus that which has gone into ice production. Nevertheless,this still requires a decreasingwith time at the variousstations. Distance(km)

916

MACDONALD

(a) MW

ET AL' FRESHWATER BUDGET OF MACKENZIE RIVER

Late Summer

6



23•> •37

3

' ,'

Late Winter

4

-•-

10

135

!

0.4

0.1

2.3

8

,



11

8

,

2.8

36

[

o.,

[

o.,•

l

1.0

I

3.0

I

•44

43

!

0.5

'•0.7

0.8

' ,

•0.9

'

•>•-0.9 ! •_ -1.9

Figure 15. Budgets(in cubickilometers) and fluxes(cubickilometers transferred between

seasons) of (a) MW and(b) SlM fortheshelf.Thenearshore zonecoincides with20 m depth, theapproximate positionof thelandfastice. Budgetsareestimated fromdatagivenin Table• for thetwosampling periods(September 1990andApril-May1991).Fluxeshavebeencalculated by

addingtheknownriverftowfortheperiod(ArcticY•edmyerstation,Wate•Surveyof Canada) andbalancingthe startingbudgetwith the endingbudget.The smallarrowsat the top of each panelfor late summerreferto ice;it hasbeenassumed herethat half the iceformedin eachzone

is melted within the gonein summer. Seetext for descriptionsof numbers.

substantial leakage of MW (31km•) to theoutershelf

Sea Ice

Melt

throughoutthe period. The volumesgivenin the top panelof Figure15ahavebeenconvertedinto equivalent SinceSIM (Figure15b) can haveboth positiveand heights(bottompanel)by usingthe respective surface negativevalues,fluxesc,fSIM in one directioncan alareasfor the regions(Table 1). For the fluxesbetween ternativelybe representedas fluxesof negativeSIM in the inner and outer boxes,however,we must give two the oppositedirectiono: by a combinationof both. In equivalentheights;one for the box being left and one Figure 15b we have adoptedthe conventionof giving for the box beingentered.The magnitudeof the fluxes the fluxes as positive SIM. As was the case for MW, comparedto the contentsof the boxessuggeststhat the nearshore SIM is utmost balanced in late winter, betempting to equate the 16km3 ofsea the flushingtime of MW on the shelfis longerin winter andit would

SIM in thewa(about6 months) thanin summer (about4 months) and iceformedwiththe 14km3 of negative that it is very short in the nearshorezone in summer ter. However,in Septemberthere remaineda significant

(1 month).

quantity ofSIM(11km3)sothatsome ofthismusthave

MACDONALD

ET AL.: FRESHWATER

BUDGET

OF MACKENZIE

RIVER

917

leakedto the outer shelfor negativeSIM must havebeen curs in winter but not in sufficientquantities to turn off imported as brine from 'the outer shelf. A combination convectionthere. Three processescombine to suppress oftheseprocesses is probably responsible for the9 km• the transport of fresh water to the outer shelfin winter: flux requiredto balancethe system,but we have no way Mackenzie River inflow is reduced, some of the inflow is of separating the two. frozen into ice, and the stamukhi zone acts as a backThe amount of sea ice produced on the outer shelf stop to impede progressof the plume in late winter. cannot be estimated with certainty from measurements Conversely,in summerthe impoundedand frozen fresh of ice cover and thickness alone. Even if one accounts for water is suddenlyreleasedand, togetherwith intensified the ice that has goneinto forming ridgesby using some inflow from the Mackenzie River in freshet, enters the appropriatefactor [e.g.,Pri•se•berg• 1988;Melli•g a•d ocean in a short time. Ice therefore amplifies the seaRiedel, 1993] and makes reasonableestimates of the sonality of the freshwater cycle on the Mackenzie shelf. amount of open water using remote sensing,a problem We can presently produce a reasonable accounting remainsin determiningexactly how much of the ice that of meteoricfreshwater (runoff) as it transitsthe shelf formed on the shelf within that year was subsequently because(1) we can distinguishthis water from seaice advected off the shelf. The difference in the outer shelf melt and estimate its quantity using isotopic composiwater column content of SIM between September and tion and (2) we havereliablemeasurements of MackenMay amounts to about143km3 allowing for the net zie inflow throughout the year. The other processthat influxof9 kms fromthenearshore (Figure15b).This affectsfreshwater balanceis the formation and melting

quantityof negativeSIM is equivalentto the growthof of ice (SIM). As wasthe casefor meteoricwater,wecan an averageof 3.0 m of seaice (equivalent liquidheight) distinguish SIM andestimateits quantityusingisotopic over the region or an ice thicknessof about 3.3 m. In

composition. However,sincewe presentlydo not have

Figure15bwe havebalancedthe SIM budgetsimplyby reliablemeasurements of the source(new ice producproducingsufficientice. While an averageof 3.3 m of tion) or sink (lossof brine from the shelf), we cannot new ice production over the outer shelf appears to be a reasonableestimate, it must be noted that any brine that has been formed and subsequentlyescapedfrom

as yet produce a complete budget. Acknowledgments.

We are extremelygratefulfor the

the region(bottomdashedarrow,Figure15b) would logistic support provided to us by officersand crew of the not leave a record in the water

and therefore

would be

missedin this budget. Melling and Lewis [1982]suggest that in regionsof .,;mallbottom slope, densewater may linger for many months and that drainage of salinewater from the shelf is of seasonalor longer duration.

Under

these circumstances

the water

column

CCGS Henry Larsen and the staff of the Polar Continental

StieffProjectw.ithout whomwe wouldnot havebeenableto collect the samples. In particular, we thank Captains Dave Johns, Steve Gomes, and Ivan C5t4 of the Canadian Coast Guard, Barry Hough and Claude Brunet of the Polar Continental Sheff Project, and Ron Sprang and Pieterre Paros., the aircraft pilots for our spring work. Adrian Abehennah and William Grieve assisted with the 5•sO determinations

record of ice production is probably a better estimate than one based on the ice content plus ice advection. at IOS. We are indebtedto the Frozensea ResearchGroup Unlessone carefully measuresthe total brine produced (IOS) for supportand loan of equipmentand to Gerry La-

andOceans, Winnipeg)forgenerously sharing (watercolumnbudgetchangeplusexport)or the total cho(Fisheries laboratory spaceat Tuktoyaktnk. Paul Squires(Water Surnew ice produced(budgetof new ice plus new ice ex- vey) providedflowdata for the Macken•.ieRiver system.The ported), one cannotproducea completebudgetas has work was funded by Indian and Northern Affairs Canada as

beenproposedfor Arctic climatestudies[World Climate part of the Northern Oil and Gas Action Program; we'thank David Stone and Russ Shearer for their continued interest in ResearchProgram,1994]. and support for our work. We acknowledgeRosalie Rutka

8. Conclusions

forthemanycarefulhoursspentonproducing diagrams and organizingthe final versionof the text. Finally, we greatly

appreciate colleagueswho took time to review e•lier drafts Measurements of salinityand $180 in the water col- and suggest numerous improvements; they are Knut Aa-

umn at the beginningand end of winter and in the ice gaard,KellyFalkner, GSte•stlund,Stephanie Pfirman,and grownthat winter showthat the Mackenzieshelfdivides Peter Scklosser. into two domainsduring winter. in the nearshore,invasion by Mackenzie River winter inflow progressively suppressesconvectionand producesa slowly advancing References

front (0.2-1 cms-1). In the offshore, seaice formation enhances the salt content of the water, thereby promoting convection. The border area between these two

domainscanbe narrow(lessthan a kilometer)andgenerally follows the stamukhi zone in late winter. The within-shelf partition in winter cannot be ignored if we

are to producemodelsthat correctlysimulatehowwater and associatedproperties transit the shelf to the interior Arctic Ocean. It •'isprobable that leakage of some fresh water from the nearshore to the offshore oc-

Aagaard, K., and E. C. Carmack, The role of sea ice and other fresh water in the Arctic circulation, J. Geophys. Res., 9.•, 14,485-14,498: 1989.

Aagaard, K., L. K. Coach:man,and E. C. Carmack, On the halocline of the Arctic Ocean, Deep Sea Res., Part A, •8• 529-545,

1981.

Anderson, D. L., Growth rate of seaice, J. Glacio!., 3• 11701172, 1961.

Beck, N., and Ko O. Muexmich, Freezing of water: Isotopic

fractionation(abstract) Chem. Geol., 70• 168, 1988.

918

MACDONALD

ET AL'

FRESHWATER

Brunskill, G. J., The Mackenzie River basin, in The Ecology oj• River Systems, edited by B. R. Davies and K. F. Walker, pp. 435-471, Dr. W. Junk, Dordrecht, The Netherlands, 1986. Burns, B. M., The climate of the Mackenzie Valley, vol. 1,

Publ. EN 57 7/24,227 I•P., Environ. Can., Toronto, 1974. Carmack, E. C., R. W. Macdonald, and J. E. Papadakis,

BUDGET

OF MACKENZIE

RIVER

Macdonald, R. W., R. Pearson,D. Sieberg,F. A. McLaughlin, M. C. O'Brien, D. W. Paton, E. C. Carmack, J. R. Forbes, and J. Barwell-Clarke, NOGAP B.6 Physical and chemical

data collected in the Beaufort

Sea and Mackenzie

River delta, April-May 1991, Can. Data Rep. Hydrogr. Ocean Sci. 10•, 154 pp., Inst. of Ocean Sci., Sidney, B.C., 1992.

Maykut, G. A., The ice e.avironment,in Marine BiologyPolar Regions, edited by R. Horner, pp. 21-82, CRC shelf estuary, J. Geoph•s. Res., 94, 18,043-18,055, 1989. Press, Boca Raton, Fla., 1985. Fissel, D. B., and H. Melling, Interannual variability of oceanographic conditions in the southeastern Beaufort Maykut, G. A., The surface heat and mass balance, in The Geophysicsof Sea Ice, edited by N. Untersteiner, pp. 395Sea, Can. Contract Rep. Hydrogr. Ocean Sci. 35, 102 463, Plenum, New York, 1986. pp., Inst. of Ocean Sci., Sidney, B.C., 1990. Giovando, L. F., and R. H. Herlinveaux, A discussionof Melling, H., and E. L. Lewis, Shelf drainage flows in the factorsinfluencingdispersionof pollutantsin the Beaufort Beaufort Sea and their effect on the Arctic Ocean pycnoSea, Pac. Mar. Sci. Rep. 81-4, 198 pp., Inst. of Ocean cline, Deep Sea Res., Part A, 29, 967-985, 1982. Sci., Sidney, B.C., 198!. Melling, H., and D. A. Riedel, Draft and movement of pack Ice Centre, Ice Thickness Climatology, 1961-1990 Normals, ice in the Beaufort Sea, April 1990-March 1991, Can. Publ. EN 57-•8•1961-•.990, 277 pp., Minist. of Supply Tech. Rep. Hydrogr. Ocean Sci. 151, 79 pp., Inst. of and Serv., Ottawa, Ont., 1992. Ocean Sci., Sidney, B.C., 1993. Ingram, R. G., Characteristics of the Great Whale River Melling, H., and R. M. Moore, Modification of halo cline plume, J. Geophys. Res., 86, 2017-2023, 1981. sonrce waters during freezing on the Beaufort Sea Shelf: Ingram, R. G., and P. Larouche, Variability of an under-ice Evidence from oxygen isotopes and dissolvednutrients, river plume in Hudson Bay, J. Geophys. Res., 92, 9541Cont. Shelj• Res., in press, 1995. 9547, 1987. Omstedt, A., and S. Wettlaufer, Ice growth and oceanicheat Jeffries, M. O., and H. R. ]•rouse, Salinity and isotope analflux: Models and measurements,J. Geophys. Res., 97, 9383-9390, 1992. ysis of some multi-year landfast sea-ice cores, northern Ellesmere Island, Canada, Ann. Glaciol., 10, 63-68, 1988. Omstedt, A., E. C. Carmack, and R. W. Macdonald, Mode•llingthe seasonalcycleof salinityin the Mackenzieshelf/ Jeffries, M. O., W. M. Sackinger, H. R. Krouse, and H. estuary, J. Geophys. Res., 99, 10,011-10,021, 1994. V. Setson, Water circulation and ice accretion beneath Ward Hunticeshelf(northernEllesmere Island,Canada), O'Neil, J. R., Hydrogen and oxygen fractionation between ice and water, J. Phys. Chem., 72, 3683-3684, 1968. deduced from salinity •nd isotope analysis of ice cores, Ann. Glaciol., 10, 68-72, 1988. •)stlund,H. G., The residence timeof the freshwater component in the Arctic Ocean, J. Geophys. Res., 87, 2035Jeffries,M. O., H. R. Krouse, W. M. Sackinger,and H. V. Water

mass structnre

a.nd boundaries

in the Mackenzie

Setson,Stable-isotope (LsO/•60)tracingof fresh,brackish, and sea ice in muir's-yearland-fast sea ice, Ellesmere Island, Canada, J. Glaciol., 35, 9-16, 1989.

Krouse,H. R., andJ. R. Mackay,ApplicationofH•sO/H•60 abundancesto the problem of lateral mixing in the LiardMackenzie River system, Can. J. Earth Sci., 8, 11071109, 1971.

3684, 1982.

•)stlund,H. G., and G. Hut, Arcticoceanwatermassbalance from isotope data, J. Geophys.Res., 89, 6373-6381, 1984.

Paton, D. W., A. Abehennah, W. Grieve, and R. W. Macdonald, NOGAP B.6 Oxygen isotopedata from water and ice cores from the Beaufort Sea, September 1990, May 1991 and September 1991, Can. Data. Rep. Hydrogr. OceanSci. 13•, 61 pp., Inst. of Ocean Sci., Sidney,B.C.,

Lehmann, M., and U. Siegenthaler,Equilibrium oxygenand hydrogen isotope fractionation between ice and water, J. Glaciol., 37, 23-26, 1991. 1994. Lepp&ranta, M., A review of analytical models of sea-ice Prinsenberg, S. J., Ice-cover and ice-ridge contributions to growth, Atmos. Ocean, 31, 123-138, 1993. the freshwater contents of Hudson Bay and Foxe Basin, Lewis, E. L., and R. G. P :rkin, The practical salinity scale Arctic, •1, 6-11, 1988. 1978: Conversionof existing data, Deep Sea Res., Part A, Redfield, A. C., and I. Friedman, The effect of meteoric 28, 307-328, 1978.

Macdonald, R. W., and E•. C. Carmack, The role of largescale underice topography in producing property flux di-

vergencein an Arctic shelf/estuary,Atmos. Ocean,29, 37-53, 1991.

Arctic

Macdonald, R. W., C. S. Wong, and P. Erickson, The distribution

of nutrients

in the southeastern

Beaufort

Sea: Im-

plications for water circulation and primary production, J. Geophys. Res., 9•, 2939-2952, 1987. Macdonald, R. W., E. C. Carmack, F. A. McLaughlin, K. Iseki, D. M. Macdonald, and M. C. O'Brien, Composition and modification

water, melt water and brine on the composition of Polar Sea water and of the deep waters of the ocean, Deep Sea Res., 16, suppl., 197-214, 1969. Scltlosser,P., D. Grabits, R. Fairbanks, and G. BSnisch,

of water

masses in the Mackenzie

shelf

estuary, J. Geophys. Res., 9•, 18,057-18,070, 1989. Macdonald, R. W., E. C. Carmack, F. A. McLaughlin, D. Sieberg, M. C. O'Brien, D. Paton, R. Pearson, Y. Liangfeng, and C. Gobeil, Oceanographicdata collected

from the Henry Larsen in the Beaufort Sea, AugustSeptember, 1990, Can. Data Rep. Hydrogr. Ocean Sci.

97, 142 pp., Inst. of OceanSci., Sidney,B.C., 1991.

river-runoff:

Mean

residence

time

on the shelves

and in the halocline, Deep Sea Res., 41, 1053-1068, 1994. SHEBA ScienceWorking Group, New program to research issuesof global climate in the Arctic, Eos Trans. A 75, 249-253, 1994.

Tan, F., and P.M. Strain, The distribution of sea ice meltwater in the eastern Canadian Arctic, J. Geophys. Res., 85, 1925-1932, 1980.

Thomas, D. J., R. W. Macdonald, and A. B. Cornford, Geochemical

mass-balance

calculations

for the coastal

Beau-

fort Sea, N.W.T., Canada. Rapp. P.-V. R•un. Cons. Explor. Met., 186, 165-184, 1986. Vetshteyn, V. Ye, G. A. Malyuk, and V. P. Rusanov, Oxygen -18 distribution in the central Arctic Basin, Oceanology, lg, 514-519, 1974.

MACDONALD

ET AL.:

FRESHWATER

Whaite, P., An automated sample line for the preparation

BUDGET

OF MACKENZIE

RIVER

Institute of Ocean Sciences,P.O. Box 6000, Sidney, B.C.,

of lSO/160 isotopeanalyses fromwatersamples, M.Sc. CanadaViL 4B2. (e-mail:robie•ios.bc.ca) thesis, 165 pp., Univ. of B.C., Vancouver, 1982. World Climate ResearchProgram, Arctic Climate System

A. Omstedt, SwedishMeteorologicaland HydrologicalInstitute, S-601 76 NorrkSping, Sweden.

Study (ACSYS), Initial implementationplan, WCRP-85, Geneva, 1994.

E. C. Carmack, R. W. Macdonald, and D. W. Paton,

(ReceivedApril 20, 1994;•evisedAugust 16, 1994; acceptedSeptember16, 1994.)