Westward Intensification of the Mean Circulation on the Bering Sea Shelf

Kinder T. H., Chapman D. C. & Whitehead J. A.

Woods Hole Oceanographic Institution, Woods Hole, Massachusetts

Journal of Physical Oceanography, Volume 16, January 1986

Geophysical Fluid Dynamics – MAST 806 Presented by: Berit Rabe and Felipe Pimenta - Fall - 2005 Outline

1. Introduction 2. Oceanographic Background of Outflow 3. Laboratory Model 4. Numerical Model 5. Conclusions 5. Study guide questions 1. Introduction Study area

Artic Ocean

Bering Strait Alaska

Bering Sea

PACIFIC OCEAN Mean: 0.5 - 1.0 Sv northward

Thesis: Bulk of the transport supplying the Bering Strait outflow occurs in the west, crossing the shelfbreak near Cape Navarin and flowing northward along the Siberian coast.

The Bering Sea (depth in meters)

Consequence: significant mean cross-isobath flow must exist from the shelfbreak (~170m depth) to the Strait (~50m depth). Dynamical implications: draw Bering Strait outflow across a shoaling bottom on the rotating earth. This can explain the apparent Westward Intensification. Arctic Ocean Bounded by land on three sides

Isobaths Bering Strait: narrow (85km), shallow (50m)

River runoff onto the shelf 40% abyssal 40% continental plain shelf (<200m) (>3500m) 2. Oceanographic Background of Bering Strait Outflow • Classify waters on the shelf by salinity

Significant feature: High S tongue, extending from the deep basin near Cape Navarin through the Gulf of to Bering Strait

Implication: There is a considerable flow of water across the shelfbreak toward the Strait. Salinity distribution at the bottom of the shelf in summer 1962-64 Siberia Alaska • Higher S (“Anadyr water”) on the Siberian side, lower S on Alaskan side • Cooler water on Siberian side, warmer water on Alaskan side

• Coachman et al. (1975): three water masses based on S – Median S of 33 for Anadyr water – 32.5 for “Bering Shelf” water – 31.5 for “Alaskan coastal” water Mean S and T distribution for Bering Strait for early summer Eastern Bering Strait section

Diomede Diomede T Island S Alaska Island Alaska

July 2003 CTD salinity section across the eastern channel of Bering Strait from Diomede Preliminary CTD islands (left) to Alaska (right). Sections, Bering Strait Note the fresh (warm) ACC on [Woodgate, 2005, Laurier the Alaskan coast. Mooring report] [“Revising the Bering Strait Freshwater Flux into the Arctic Ocean”, R. A. Woodgate and K. Aagaard, accepted for publ. in Geophysical Research Letters, Dec 2004] • Available direct current measurements (in 1986!): inadequate to estimate transports and to infer locations of currents across the shelf

• Use of estimates of river runoff and S as a water mass tracer synthesis of the shelf water balance can be constructed

• Bering Strait mean flow: driven by sea level difference between Pacific and Atlantic. Estimated transports: 2.2 and 1.6 Sv northward.

Hypothesized western boundary current probably supplies the northward flow of higher S water in the western 2/3 of the Strait.

Velocity (cm/s) obtained by lowering a current meter from an anchored ship during Aug 1967. Shaded areas: southward flow

• High speed jet (low S coastal water) on the Alaskan side • Northward flow (higher S Anadyr water) on the Siberian side. • Bering Strait flow: annual cycle with mean: ~0.6Sv • Two paths to the Strait: • Eastern current appears to be weak in the south but persistent throughout the year. • River runoff only accounts for a small fraction of the Bering Strait outflow.

• This coastal flow plus river Schematic of currents runoff accounts for only 15% that supply the Bering Strait outflow. Numbers of the Bering Strait outflow. are transports in Sv, and • Roughly corresponds to the are intended to suggest low S Alaskan coastal water relative sizes of the flows only. found on the Alaskan side of the Strait. • Northwestward flowing Bering Slope Current: 5 Sv, but parallel to the slope until Cape Navarin • Strong S signal close to Siberian coast, but no measurements. • Most likely location of remaining 85% of the overflow supply: along the Siberian coast, in a western boundary current (coincident with S tongue). • High S Anadyr water does not contribute 85% to the overflow Schematic of currents transport in the Strait. that supply the Bering • Some of the water has to undergo Strait outflow. Numbers are transports in Sv, and significant modifications before it are intended to suggest reaches the Strait. relative sizes of the flows only.

• The western boundary current contribution corresponds to the Anadyr and Bering Shelf waters. 3. Laboratory model • Concept of western-intensified source- sink flow was applied to Bering Strait Shelf: – Placement of a local sink near the center of a coast. – In the (topographic) β-plane analogy, the sink is in the center of the zonal northern coast. – Bottom with a slope of 1:5 was placed in 2m-turntable. – Rotation period of 15s (fluid experienced topographic β-effect). • Walls corresponds to eastern, southern, western walls of topographic β-plane. • Tank covered with transparent plastic (little or no wind stress). • Moveable sink placed at northern shoreline. • Run commenced after fluid spin up, pump (sink) was started, evolution of the dye recorded by camera above the tank

Diagram of laboratory apparatus. Left: top view, right: side view (with the vertical exaggerated by 2). Arrows show individual locations of the sink, circle in the lower left corner is the gravity return flow through a hole in the false bottom. West sink East

10cm Top view of the evolution of a rectangular grid of dye in a rotating fluid with a sloping 10s 20s bottom. Dye has been in place for at least 60s before the first photograph was taken. Times after start of the sink are indicated in 30s 40s figure. Period of rotation: 15.05s. Volume flux of the sink: 24cm3s-1

50s 120s West sink East Current flowed up the western wall, along the shallow part of 10s 20s the shelf to the sink. Little or no current elsewhere, neither to the 30s 40s right (east) nor in the interior of the fluid.

50s 120s sink

5s 25s

Same type of experiment as in previous figure except sink is near the west coast. Period of rotation: 15.5s, volume flux: 16cm3s-1.

• Same experiment, except sink in the West. • Same result as in last experiment. • Most of the water moving up the shelf did so in the western boundary current. • Lack of currents elsewhere. sink

5s 25s

Same type of experiment as in previous figure except sink is near the east coast. Period of rotation: 15.4s, volume flux: 53cm3s-1.

• Same experiment, except sink in the East. • Same result as in last experiments. • Most of the water moving up the shelf did so in the western boundary current. • Lack of currents elsewhere. J.E.Overland et al., Direct evidence of northward flow on the northwestern Bering Sea shelf. J. Geophys. Res., 101, 8971-8976, 1996.

Positions of buoy 7210 from June 1 to August 13, 1994 • Two classical length scales for width of western boundary current (linear flow): – Stommel scale (due to bottom friction):

2 1/ 2 Ws = (" /2 f tan #) – Munk scale (due to lateral friction of interior flow):

1/ 3 %1/ 2 1/ 3 %1/ 2 Wm = 2"(# /$) 3 = 2"(#d / f tan&) 3 ! • f: Coriolis parameter • ν: kinematic viscosity of the fluid

• tanα: bottom slope with β0 = f*tan α/d ! • d: depth of fluid • Plugging in values: – d=10cm, f=4π/15s-1, tanα=0.2ν=0.01cm2s -1

– Ws = 0.4cm – Wm = 3.1cm

• Observed width of western boundary current: – At d=10cm: ~10 cm in exp.1 – ~7-8 cm in exp. 2 – ~15 cm in exp. 3

• Munk scale closer to observations, but still too small. • Third possible length scale: – Inertial boundary layer, nonlinear advection dominates – Degree of non-linearity: compare bottom friction time scale to advective time scale – ratio: ~0.77 advective effects are of same order as bottom friction effects – inertial boundary layer appears possible! • Steady-state width of inertial boundary current: 1/ 2 WI = (VI d / f tan")

• VI: interior velocity normal to boundary layer

! • For inertial boundary layer of width WI=10cm, an interior velocity -1 of VI=-1.7cms is required.

• Large westward displacements of north-south dye lines were not observed in experiments.

• Internal scale does not apply here, or steady state was not reached in experiment.

• Calculated volume flux of western boundary current, compare to that of the pump: western boundary current contains the volume flux of the pump, current carries all the water on the shelf that goes into the sink. • Conclusion of lab experiments:

– Laboratory experiment shows a western boundary current with a width somewhat greater than that expected by frictional theory.

– Boundary current may be inertial, then it must still evolve in time.

– All of the flow into the sink comes via the western boundary current, and from nowhere else. 4. Numerical model 4. Numerical model

Goal: Investigation of the dynamics of the circulation of the Bering Sea and of Artic Ocean the laboratory experimental results Alaska Siberia Neither friction nor nonlinear effects seemed to account for the western current width. Bering Sea

Simple linear numerical model to PACIFIC OCEAN illustrate the existence of the Western Boundary Current (WBC):

Governing equation (~advection - diffusion) #$ ru " fv = "g " $ f dh dx ' *# $ 2dh dx ' *# #x h " 2# + & ) +& ) = 0 #$ rv % r ( *y % h ( *x fu = "g " #y h Solved with an upwind differencing method of "(uh) "(vh) Fiadeiro and Veronis (1977) + = 0 ! ! "x "x

#$ #$ uh = " vh = #y #x ! transp. stream-function

! 4a. Simulation of Lab. Experiment

Model characteristics and parameters 4" "f #1 f = = 0.8378s#1 r = = 0.065cm.s T 2 h(x) cm outflow: 5 cm

! ! " =1 ! " = 0 h(x) = 0.2x

Return flow Lab. appears to develop " =1 ! a wider current (8-15 cm) ! ! ny=39, nx=21, dx=dy=5 cm #$ uh = " = 0 #y ! Scale given by: r/fhx ~ 0.39 cm, model does not Results resolved details of it

2 ! 1 Cross isobath flow analogous to WBC 1st simulation described by Stommel (1948)

1 2 Northern coast: current turns and spreads into an alongshelf boundary layer, analogous to ATW solution (Csanady, 1978) 4a. Simulation of Lab. Experiment

Bathymetric domain " = 0 " =1 cm h(x) 5 cm 2nd simulation

! ! ! h(x) = 0.2x " = 0 " =1 " =1 3 !

! flat bottom " 2# " "(vh) 2 = = 0 ! "x "x "x !

Could the location of the gap in the offshe!lf wall be responsible for the cross isobath flow?

3 Inflow spread over the entire offshelf boundary, turning sharply toward west

Has little effect on the interior flow! 4b. Mean Circulation of Bering Sea

80 km, topographic effects are larger 40 m f =1.27 "10#4 " #1$10%13 cm%1s%1 1Sv "1 f dh dx r = 0.1cm.s " # # 5 $10%12cm%1s%1 T h " = 0 ! 40 + 4 "10#4 x x < 400km ! h(x) 200 + 0.035(! x # 4 "105 ), 400km < x < 480km slope 1 ! " = 3000 x > 480km

" 2# ! 2 = 0 ny=51, nx=28, dx=dy=20 km ! "x

1 Predicts strong WBC as in lab. 3rd sim!ul ation experiment 2 Inflow over entire deep basin, turns to Siberian Coast 1 3 Stommel scale of bc: r / f dh dx = 20km (shelf ) 2 r / f dh dx = 0.22km (slope) 3 Very little transport occurs outside the boundary current ! 4b. Mean Circulation of Bering Sea y Velocity structure

x 3rd simulation

-1 ~10-20cm.s increase due to Siberian coast the decrease in bottom depth

Momentum balance (x=260 km) p " #f$ " #g%

“-u geostrophically balanced “Along-current pressure gradient nearly with cross-current sea slope…” balances the drag due to bottom friction…” ! X-momentum balance

#$ rv #$ ru fu = "g " " fv = "g " #y h #x h Y-momentum balance

! ! 4b. Mean Circulation of Bering Sea y

x 3rd simulation

~10-20cm.s-1 Siberian coast

“approaching the northern boundary the along-shelf p " #f$ " #g% Momentum balance (y=100 km) current v is nearly in geostrophic balance with the offshore sea slope …” “friction effects become large near the boundary, but dimin!is h outside of the geostrophic flow…” X-momentum balance

#$ ru " fv = "g " #$ rv #x h Y-momentum fu = "g " balance #y h

! ! 5b. Mean Circulation of Bering Sea

Observations suggest that there may be a substantial flow in the deep part of the Bering Sea (4-5 Sv) (Kinder, 1975)

The effect can be qualitatively modeled by an imposing transport trough Alaska Peninsula

4th simulation 3

“Suggest that Bering Strait outflow may draw most of its water from the entire 2 uniformly slope and little from 1 Imposed deep basin…” Seaward of 1000 m

1 Deep inflow fills deep basin and supplies the transport for the Bering Strait outflow 3 Shelf flow is virtually unaltered by deep 2 Current parallel to the shelf break is mean flow generated and held there by deep flow 5. Conclusions 5. Conclusions

Based on observations, major part of supply Aspects that were excluded from the to the mean Bering Strait outflow occurs as modeling work: a Western Boundary Current (WBC) -density variations: The existence of that WBC could not be Shelf is strongly stratified, this could partially confirmed by direct measurements, but mitigate the effect of the varying bottom were explained in dynamical terms Anadyr could modify pressure gradients on Siberia coast Analogy to Stommel’s (1948) explanation for WBC: meridional gradient of Coriolis -wind stress and atm. pressure: parameter is replaced by the gradient of depth relatively high values and large fluctuations (days to years). It can dominate the variability near the Strait.

Laboratory and numerical experiments -tides: demonstrated the topographic β effect: kinectic energy dominate on south but is less Critical shelf features were included: important over northern shelf (r can be improved) -topography: Gulf of Anadyr: induce curvature -shoaling bottom, earth rotation, friction and St. Lawrance Island: split flow forced outflow Siberian coast: not a vertical wall: could move current away off coast “envisaged mechanism is physically real for a sink at the coast and the numerical model Other oceanographic situations: Western showed it can be applicable to the Bering Mediterranean: flow forced across shoaling Sea” topography of Strait of Gibraltar forms a narrow boundary against African slope 6. Study guide questions Study guide questions 1a)The source/sink system drives the mean circulation

1) What drives the mean circulation in both laboratory 1b) No. In the large scale case, winds play and numerical model, i.e. , does a western boundary the role of the source/sink. But essentially, current always require a wind stress? the existence of the WBC requires the presence of bottom friction in the large scale models. In Bering Sea, the intensification can be seen as follows:

Scale of u width: Cross-shore intensification: "# $ fh ' $ 2h ' assume vh = $ 0 " 2# + & x )# *& x )# = 0 "x % r ( y % h ( x 2 $ fhx ' $ 2hx ' " $ "# ' $ fh ' "# " # + & )# y *& )# x & ) + & x ) = 0 % r ( % h ( "y % "y ( % r ( "y ! "# # $ fh ' ! assume vh = $ 0 " (uh) + & x )( "uh) = 0 "x #y % r ( ! ! "u # fhx & # fhx & h + % ( uh = 0 y "yy + % (( )uh) = 0 "y $ r ' ! $ r ' ! u'+( fhx r)u = 0 x $ r ' $ 1 ' ! # fh r y ( x ) The shelf break ! u"# & )* yy& ) u(y) " Uoe INFLOW drives the % fhx ( {% h(x)+( ! mean circulation 12 3 cte cte hand waving…

! ! Study guide questions

2) What baroclinic current speeds and directions would expect from the section shown across the Bering Strait?

density pressure X

X

[Woodgate, 2005, Laurier Mooring report]

"V #g % "$( #9.8 % 22.5 # 25( = ' * + #4 + #' * "z $o f & "x ) (1025)(1.3,10 ) & 30000 )

$ 22.5 # 25' "V = & ) 40 * 0.24m.s#1 ! % 30000 (

! Study guide questions

3) What physical process are considered in the 4) What is the governing equation that forms laboratory and numerical models and how do the essence of the model applied to the they relate to non-dimensional parameters we Bering Sea and how does it relate to what derived in class? we discussed in class?

Laboratory: The essence of the physics is given by a vorticity equation in terms of the -barotropical model, no wind stress at surface transport stream-function: nor density gradients -shoaling bottom, -tank rotation, 2 $ fhx ' $ 2hx ' -bottom friction " # + & )# y *& )# x = 0 -forced inflow/outflow % r ( % h ( -note: non-linearity is present and experiments evolve in time Any ideas on how to relate it to the Numerical Model ! potential vorticity conservation? We have “stronger control” so we can set:

Ro <<1 R <<1 " = 0 steady state solution oT "t

! simple formulation of friction ! ! Study guide questions #$ " fv = "g #x 5) What are the dominant force and vorticity $% rv balances in different part of the domain fu " #g # $y h ! y $% fu " #g $y ! x ru #$ = "g h #x !

!

2 $ fh ' $ 2h ' x x % 2hx ( " # + & )# y *& )# x = 0 " # 0 " $' *" = 0 % r ( % h ( y xx & h ) x $ fh ' " # 0 " + & x )" = 0 x yy % r ( y ! !

! Study guide questions

6) Does the basic geostrophic flow follow the bottom contours?

Not in the entire domain. Near the shelf break and on the north coast the flow follows the isobaths.

The western boundary current on the other hand, even demonstrating a near- geostrophic nature, explicitly flows against the slope Study guide questions

7) Which physical process breaks the geostrophic constraint, that is, what is the origin of the relative vorticity term in the main governing equation?

#$ ru "(uh) "(vh) # "u "v& )u dh " fv = "g " % + ( = + = 0 # $ "x "x' h dx #x h #$ #$ "x "x #$ rv uh = " vh = fu = "g " #y #x "v "u "v % 1 "(vh) v dh( #y h # $ = ' # * "x "y "x & h "x h dx ) ! 2 ! "v " #! r "u f = g + y x y h y ! " " " " 1 $ "vh "uh' $ fu v ' 2 ! # = h + 2 h "u " # % 1 "v "(1/h)( & ) & x x ) f = $g $ r' + v * h % "x "y ( % r h ( "x "x"y & h "x "x ) % % (( # "u "v & r "u # 1 "v v dh& 1 " "# " "# % fu v ( f 0 r ' $ ' $ ** = ' hx + 2 hx * % + ( = + ) % ) 2 ( h x x y y & r h ) $ "x "y ' h "y $ h "x h dx ' ! & " " " & " )) ! # u dh& r # )v )u& # rv dh& $ 2 2 ' f " # " # $ fhx hx ' " % ( = " % " ( + % 2 ( + + # * 2 # = 0 $ h dx ' h $ )x )y' $ h dx ' & 2 2 ) & y x ) ! ! % "x "y ( % r h ( r $ "v "u' dh $ fu rv ' & # ) = & + 2 ) h % "x "y( dx % h h ( ! !

! Study guide questions

8) Why does the Rossby radius of deformation not enter the dynamics that are certainly dominated by rotation and ambient vorticity gradients?

Specifically for the WBC, as seen its width it’s given by the Stommel scale of the boundary current

#( fhx r)y u(y) " Uoe

r For a decay of 1/e: y = fhx !

Specifically !fo r the WBC, as seen its width it’s given by the Stommel scale of the boundary current End