SPIN-UP OVER STEEP TOPOGRAPHY S U B M A R I N E A N D T H EEFFECTS CANYON by Ramzi Mirshak B.Sc. (honours) M c G i l l University, Montreal, Canada, 1998 A THESIS SUBMITTED THE I N PARTIAL F U L F I L L M E N T OF REQUIREMENTS MASTER F O R THE D E G R E E OF OF SCIENCE in T H E FACULTY OF G R A D U A T E STUDIES Department of EARTH A N D OCEAN SCIENCES We accept this thesis as conforming to the required standard T H E UNIVERSITY O F BRITISH C O L U M B I A May 2001 © Ramzi Mirshak, 2001 OF A In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Earth and Ocean Sciences The University of British Columbia 2075 Wesbrook Place Vancouver, Canada V 6 T 1Z1 Date: Abstract Submarine canyons are common bathymetric features that cut into the continental shelf from the continental slope. During upwelling favourable conditions over the continental shelf, submarine canyons display enhanced upwelling, heavily impacting shelf-slope mass exchange. In this thesis, laboratory experiments were conducted to quantify how velocity, stratification and rotation affect canyon upwelling. Currents were forced by changing the rotation rate of an already rotating tank. Over time, the forced currents evolve, or spin-up, until the fluid within the tank is once again rotating at the same rate as the tank itself. A difference in the spin-up behaviour is observed when a canyon is not present in the laboratory topography. The dynamics that govern flow evolution when a canyon is not present are Ekman suction and diffusion, both of which can be affected by a sloped bottom boundary layer. A numerical model is developed that replicates the spin-up of a stratified fluid over the changing slopes of the laboratory topography in the absence of a canyon, incorporating both the decay of Ekman suction and the change in the diffusion of momentum which occur as the boundary layer flow is arrested by buoyancy effects. Observed spin-up with a canyon present is compared to predicted spin-up without a canyon. The difference measures a force imposed by the canyon, which is found to be proportional to U ' f°' /N, 2 5 5 where U is the shelfbreak velocity, / is the Coriolis frequency and N is the buoyancy frequency. This force relates to the flux of water crossing the shelfbreak through the canyon. The results are applied to Astoria canyon, suggesting that in a strong upwelling event, the flux of water through the canyon is 1.3 x 1 0 m s , nearly 5 3 -1 20 times larger than wind-induced upwelling over a similar length of the shelfbreak. ii Table of Contents Abstract ii L i s t of Tables vi L i s t of F i g u r e s vii Acknowledgements 1 Introduction 1 1.1 Background 1 1.1.1 Flow Characteristics in and around Submarine Canyons 4 1.1.2 The Rossby Number 8 1.1.3 Upwelling and Submarine Canyons 9 1.1.4 Modeling of Canyon Flow in the Laboratory 1.2 2 ix 11 Method of Solution 12 1.2.1 The Effects of Rotation and a Canyon on Flow Evolution 12 1.2.2 Organization of This Thesis 14 L a b o r a t o r y M e t h o d s a n d D a t a A n a l y s i s Techniques 15 2.1 Scaling Analysis 15 2.2 Experimental Trials 20 2.3 P I V techniques 22 2.3.1 24 Error Analysis and Propagation iii 3 Laboratory Spin-up W i t h o u t a Canyon 28 3.1 28 3.1.1 The Spin-up of a Homogeneous Fluid Over a Flat Bottom . . . . 31 3.1.2 The Spin-up of a Stratified Fluids Over a Flat Bottom 32 3.1.3 Sloping Ekman Layers 33 3.2 Objectives and Outline of This Chapter 34 3.3 Spin-up of a Homogeneous Fluid Over a Steep Slope 35 3.3.1 Boundary Layer Flow 36 3.3.2 Effect of a Slope on Ekman Suction 43 3.3.3 Effect of Ekman suction on Spin-up 45 3.3.4 Solving the Homogeneous Problem Numerically 46 3.3.5 Comparing Theoretical Predictions To Laboratory Results 3.3.6 Fitting Curves to the Velocity Data 49 3.3.7 Flow Comparisons For the Homogeneous Spin-up Model 52 3.4 3.5 4 Introduction to the Spin-up Problem . . . . 49 Spin-up of a Stratified Fluid Over a Steep Slope 52 3.4.1 Solving the Stratified Problem Numerically 57 3.4.2 Flow Comparison For the Stratified Spin-up Model 59 Summary 60 Effects of a S u b m a r i n e C a n y o n o n the L a b o r a t o r y F l o w 66 4.1 Introduction 66 4.2 Qualitative Evidence of Canyon Effects 68 4.3 Drag 70 4.3.1 4.4 Determining Canyon Drag Contributions To The Evolving Flow . 71 Analysis of Results 75 4.4.1 78 Shelfbreak Velocity iv 5 4.4.2 Stratification 78 4.4.3 Rate of Rotation 80 4.4.4 Combination of Dimensional Parameters 80 4.4.5 Non-Dimensional Parameterization of Cd 81 4.5 Case Study: Predicting the Upwelling Flux Through Astoria Canyon . . 83 4.6 Summary 83 Discussion a n d Conclusions 85 5.1 Laboratory Methods 85 5.2 Laboratory Spin-up With and Without a Canyon 86 5.3 Scatter in Data Measurements 87 5.4 Implications of Findings 88 5.5 Future Work 90 Bibliography Appendix A ' Canyon F l u x Measurements v 91 97 List of Tables 2.1 Scaling of parameters used in laboratory modelling 19 2.2 Trials performed in the laboratory 19 A.l Drag and upwelling flux predictions for trial 1 98 A.2 Drag and upwelling flux predictions for trial 4 98 A.3 Drag and upwelling flux predictions for trial 5 99 A.4 Drag and upwelling flux predictions for trial 6 99 A.5 Drag and upwelling flux predictions for trial 7 A.6 Drag and upwelling flux predictions for trial 8 vi . . . 100 100 L i s t of F i g u r e s 1.1 Cartoon of Generic Shelf Topography 2 1.2 Submarine canyons of the west coast of North America 3 1.3 Flow in and around a submarine canyon 5 1.4 Vortex stretching as a current travels over topography 6 1.5 Wind-driven coastal upwelling 10 2.1 Laboratory set-up 16 2.2 Bathymetry of Astoria canyon 18 2.3 Creation of a density gradient in the laboratory 21 2.4 Flow chart of P I V analysis 23 2.5 Plots used to determine particle image diameter 26 3.1 Spin-up in the laboratory 29 3.2 Secondary circulation of a homogeneous fluid in a cylinder 30 3.3 Secondary circulation of a homogeneous fluid in the laboratory 35 3.4 Rotation of coordinates 37 3.5 Contribution of shear to Ekman suction over a steep slope 3.6 Vertical velocity contribution of Ekman suction on a slope 3.7 Vertical Ekman pumping into a sloped boundary layer as a function of shear 44 3.8 Example of a velocity profile early in laboratory spin-up 46 3.9 Topography in the laboratory and in the numerical model 48 . 42 43 3.10 Laboratory measurements showing nonlinear spin-up in early time . . . . 50 3.11 Spin-up of a homogeneous fluid 53 vii 3.12 Nomenclature used in stratified spin-up coordinates 58 3.13 Spin-up for trial 1 without a canyon 61 3.14 Spin-up for trial 4 without a canyon 62 3.15 Spin-up for trial 6 without a canyon 63 3.16 Spin-up for trial 7 without a canyon 64 3.17 Spin-up for trial 8 without a canyon 65 4.1 Flow evolution with and without a canyon in the laboratory 67 4.2 Pathlines of flow past a canyon 69 4.3 Freebody diagram of forces on the 4.4 Calculation of drag by the canyon from model output 74 4.5 Method for determining error in 76 4.6 Drag forces as a function of the shelfbreak velocity squared 4.7 Relationship between stratification and drag and rotation and drag 4.8 Relationship between canyon drag, velocity, stratification and rotation . . 82 5.1 Comparison of upwelling predictions 89 fluid 72 fluxes viii 77 . . . 79 Acknowledgements I would like to give a heartfelt thank-you to my supervisor, Susan Allen. Her insight and enthusiasm have been a wonderful inspiration. I could not have hoped for a better supervisor. Rich Pawlowicz has been a tremendous help. From my many hours in the first row of his lectures, to the time in his office working on theoretical ideas and on to his suggestions on the text, his teachings and guidance have been superb. Noboru Yonemitsu gave tremendously useful and enthusiastic contributions to the developments in the laboratory. His selfless contribution of time and energy, his curious nature and his cycling advice have all been extremely helpful. David Jones and Doug Latornell kept the table turning in the laboratory. Tawnya Peterson is personally responsible for the fact that my head is still firmly attached to my shoulders. I dedicate this project to my parents and sisters who have given me their care and support throughout my time in Vancouver. ix Chapter 1 Introduction Submarine canyons are bathymetric features that cut into the continental slope, often extending past the shelfbreak and incising the continental shelf. During upwelling favourable conditions, the upwelling signal is enhanced in these canyons. The upwelled water is rich in nutrients, and as a result canyons are regions of high biological productivity. Upwelling events last for several days. A n upwelling episode contains an early transient phase and a longer lasting, steady advection phase. The transient phase is well described by linear dynamics (Allen, 1996) but the second phase is highly nonlinear and difficult to characterize. This thesis examines how the variation of current velocities, density stratification and Coriolis frequency control the amount of water that is upwelled through a submarine canyon. 1.1 Background The mesoscale current field on the continental shelf is forced by winds, gravity and topography. Gravity is a body force while wind and topography are felt through the boundaries of the ocean. Aside from near seamounts and islands, the effects of wind and topography are combined only at or near the edges of the ocean basins in a region known as the continental margin. This margin is comprised of the continental shelf, the shelfbreak, the continental slope, and the continental rise (Figure 1.1). The continental shelf extends seaward from the coast for a distance that varies from 3 km to 500 km. The edge of the continental shelf is marked by a sudden and considerable 1 Chapter 1. 2 Introduction XL Shelf-break Continental Shelf \ Continental Slope Continental Rise \ / Abyssal Plain Figure 1.1: Cartoon of generic shelf topography described in the text. Submarine canyons cut into the continental slope increase in the bottom slope. The location of this change in slope is known as the shelfbreak and is believed mark the edge of the sialic continental mass. The steep region downslope of the shelfbreak is the continental slope, where ocean depth increases quickly to over 1 km. The continental rise marks a decrease in the slope and connects the continental slope to the abyssal plains of the deep ocean. To lowest order, the ocean is dominated by a geostrophic balance at low frequencies, which effectively implies that horizontal pressure gradients in the water column are balanced by the rotational influences of the planet. One of the effects of this balance on the coastal ocean is that currents generally follow isobaths. Due to the steepness of the continental slope, exchange across the shelfbreak is limited to secondary circulation. A r eas of complex shelfbreak geometry, however, will have stronger cross-shelfbreak currents since geostrophic balances tend to break down in these regions. A n example of such a feature is a submarine canyon. Submarine canyons are common along ocean edges (Figure 1.2). They are conduits for sediment deposition to the deep ocean, and are often sites of genesis for turbidity currents. Chapter 1. Introduction 3 Figure 1.2: Submarine canyons of the west coast of North America. Reproduced with permission from Hickey 1995. Chapter 1. 4 Introduction Submarine canyons are only present in areas where the slope is greater than 3° (Thurman, 1994) and are believed to be formed by erosional action of these turbidity currents. Where these features reach past the shelfbreak and cut into the continental shelf, they are sites of enhanced upwelling and downwelling onto and off of the continental shelf and as such are important in determining watermass distribution over the continental shelf and slope as well as mass exchange across the shelfbreak. 1.1.1 Flow Characteristics i n and around Submarine Canyons Flow patterns within a submarine canyon are driven by the geostrophic pressure gradient (Freeland and Denman, 1982). The interaction of the density stratified ocean with the steep and complex topography results in a strongly non-linear, baroclinic, threedimensional flow-field (Figure 1.3), the details of which still remain unclear. In order for flow characteristics to be affected by topography, deviations from geostrophy within the canyon must be strong. Canyon effects are strong over more narrow canyons (those whose width is less than half of the Rossby deformation radius ; Klinck, 1988), although 1 recent measurements over Barkley Canyon (Vindeirinho, 1998; see Figure 1.2) showed that upwelling effects can be measured even over wider canyons. It also seems that the width of a canyon and its radius of curvature where it meets the shelfbreak play an important role in the volume of upwelling through a canyon (Allen and Hickey, in prep.). One of the best field datasets examining flow through a canyon was gathered over Astoria Canyon (Hickey, 1997), which cuts into the continental shelf west of the mouth to the Columbia River. Results of this field work demonstrated that while currents in the upper layers of the ocean (100 metres above the canyon-rim) are barely affected, The Rossby deformation radius is equal to -dgD/f where g is gravity, D is the depth and / is the Coriolis frequency. 1 Chapter 1. 5 Introduction Canyon Head Figure 1.3: Schematic diagram showing the flow features in and around a submarine canyon. if at all, by the presence of a canyon below, the water below this depth is stretched as it passes over the canyon rim (Figure 1.4). This stretching results in an addition of cyclonic vorticity to the flow field. Submarine canyons have an eddy at rim-depth (Allen et ai., submitted). Laboratory results suggest that barotropic controls of the rimdepth eddy include the strength of the velocity field (Perenne et al, 1997), while scaling arguments show that the shape of the canyon (triangular- versus square-shaped rim) is also important in baroclinic environments (Allen and Hickey,in prep.). Geometry also affects the nature of the upwelling environment. Some canyons, such as Juan de Fuca canyon, cut far into the continental shelf. (In the case of Juan de Fuca Canyon, it cuts all the way into the shelf, eventually becoming Juan de Fuca Strait.) Chapter 1. Introduction 6 Pressure Gradient Figure 1.4: Vortex stretching experienced by a fluid as it travels over topography. The topography in this case is meant to be the upstream lip of a canyon. The cross-section is oriented such that you are looking towards shore. The wavy line represents the ocean surface. Lines of constant depth (dotted lines) and constant density (dashed lines) are also shown. Arrows are also present giving directions of prevailing winds and currents. The pressure gradient coming out of the page is representative of the northern hemisphere (i.e. it would be into the page in the southern hemisphere). The fluid column (box filled by diagonal cross-hatching) increases its vertical length as it travels over the steep topography, resulting in an addition of local vorticity to the flow field. Chapter 1. Introduction 7 This canyon exhibits stronger, more consistent upwelling (Allen, 2000) while a tapered, more normal canyon such as Astoria Canyon exhibits pulses of upwelling (Hickey, 1997). Normal canyons appear to only become important when the along-shore currents are strong (Allen, pers. comm.; Allen and Hickey, in prep.; Hickey, 1997). The Astoria dataset has been used to calibrate numerical models. These numerical models revealed that the advection of deep water onto the shelf was related to the canyon dimensions and stratification (Klinck, 1996) and that upwelling through a submarine canyon is an order of magnitude stronger than the on-shore flux at the shelfbreak in the absence of a canyon (Allen, 1996). Submarine canyons not only affect vertical currents (upwelling and downwelling) but horizontal ones as well. Bottom currents downstream of the Hudson Shelf Valley (Han et al, 1980) and Astoria canyon (Hickey, 1997) were observed to be in phase with flow variations within the canyon during upwelling-favour able conditions. In a downwellingfavourable scenario, (one where the currents on the shelf are travelling in the same direction as the shelf waves,) the flow patterns change. When downwelling currents travel past the Grand-Rhone Canyon in the Mediterranean Sea, the flow does not extend onto the shelf downstream (Durrieu de Madron, 1994). Since upwelling takes place when the flow is travelling in the direction opposite to the propagation of shelf waves, a train of standing Rossby shelf waves forms downstream of the canyon. These waves do not form when the flow is downwelling favourable. The formation of these waves and their effect in a homogeneous fluid experiencing oscillatory forcing has been studied numerically and in the laboratory (Perenne et ah, 1997). The standing waves impose a drag on the flow during upwelling events, but there is no equivalent set of waves formed during the downwelling phase of the oscillatory forcing. The unbalanced drag results in a rectified flow travelling in the downwelling-favourable direction along the shelf. Chapter 1. 8 Introduction Tidal oscillatory forcing can result in a rectified flow travelling in the direction of shelf waves propagation over canyonous bathymetric regimes as may be inferred by the current regime near Grand-Rhone Canyon (Durrieu de Madron, 1994). Off the west coast of North America, however, the upwelling and downwelling patterns are seasonal, and while strong tidal oscillations are present (Foreman and Thomson, 1997), they do not affect the mean currents strongly enough to eliminate the upwelling fluxes during the summer months. It is therefore possible that, in this region of the world, standing waves remain present over the continental shelf during the summer months, while upwelling persists. 1.1.2 The Rossby Number The degree to which rotation is important in a flow is expressed in a non-dimensional number known as the Rossby number. The Rossby number is defined as R = I J (") where U is a characteristic velocity scale, L is a characteristic length scale and / is the Coriolis frequency, / = 2fisino:, where VL is the rotation rate of Earth and a is the latitude. The Rossby number results from a scaling analysis of the momentum equation, advection — + uV • u rotation + fk x u = ^2 Forces, (1.2) where u is the velocity and k is the unit vector in the vertical direction. When examining ocean currents at low frequency, it can be assumed that du/dt is small. This leaves only the advection and rotation terms, which are labeled above, on the left hand side of the equation. By scaling (1.2) with appropriate length and velocity scales, L and U, respectively, the ratio of the two terms above will be the Rossby number. When R <C 1, Chapter 1. 9 Introduction a reasonable assumption is that advective terms are unimportant and the geostrophic balance will govern flow behaviour. 1.1.3 Upwelling and Submarine Canyons Currents on the continental shelf are governed by geostrophy to lowest order. (A flow velocity of lCT m - s , a Coriolis frequency of 1 0 s , and a lengthscale of 10 m, the 1 -1 _ 4 _ 1 5 Rossby number is 10~ , indicating that rotation is far more important than advection 2 for these currents.) Along the shelfbreak a front is present that marks the shelf-slope interface and is an effect of the geostrophic balance. Upwelling currents are generated primarily as a secondary current caused by the divergence in surface currents away from the coast due to the southward wind stress (Figure 1.5). Off the west coast of North America, winds blowing to the south result in westward transport in the upper Ekman layer. By continuity, the depleted water at the coast is replaced from below, resulting in upwelling. The wind-driven offshore transport also creates an offshore pressure gradient, resulting in a southward flow over the continental shelf. Southward currents in this region of the world are typical in the summer upwelling season. When this southward travelling current encounters a submarine canyon, due to the curvature of canyon bathymetry (typically less than or on the order of 10 km), geostrophy begins to break down. Advective terms increase and the current directions deviate from the regional isobaths. The region of continental shelf extending along the west coast of British Columbia, Canada and Washington, USA is cut by numerous submarine canyons (Figure 1.2). In upwelling-favourable scenarios, nutrient-rich water (Freeland and Denman, 1982) is advected onto the continental slope across the entire shelfbreak, with enhanced advection present in canyons (Hickey, 1997; Allen, 1996; Klinck, 1996) resulting in enhanced primary productivity in these regions (Allen et al, 2001; Vindeirinho, 1998; Mackas et al, 1997). Analysis of chemical properties suggests that the waters Chapter 1. Introduction 10 Figure 1.5: Wind-driven coastal upwelling as described in the text. Reproduced with permission from Thomson (1981). Chapter 1. Introduction 11 being upwelled through Juan de Fuca canyon have a source depth of 400 m (Freeland and Denman, 1982), although the geometry of Juan de Fuca canyon makes it a special case (Allen, 2000) . Scaling analysis predicts source waters for upwelling in other canyons to come from a depth of about 300 m (Allen and Hickey, in prep.). Throughout most of the coastal ocean, primary productivity is dependent primarily on nutrients . Nutrient levels in the ocean peak at approximately 500 m depth, tapering to 2 a near-zero level at the surface. One of the sources of nutrients into the surface waters is upwelling, which brings deep, cold, nutrient-rich water to the euphotic zone. Downwelling results in the opposite, arresting fertilization of surface waters from below. Upwelling and downwelling can be forced by wind-driven currents and by topographically-driven currents. Topographically forced upwelling is strong over the continental shelf and is one of the principle reasons for high biological productivity in coastal regions around the world. 1.1.4 M o d e l i n g of C a n y o n F l o w i n t h e L a b o r a t o r y Due to the steep nature of canyon topography and the baroclinic controls on the threedimensional flow field associated with these features, the flow in and around submarine canyons is complex. In addition to field studies (Hickey et al, 1986; Hickey, 1997), analytical models (Chen and Allen, 1996) and numerical models (Allen, 1996;Klinck, 1996; Klinck, 1989), laboratory experiments (Allen et al, in prep; Boyer and Davies 2000; Perenne et al., 1997) have been used to realize flow around submarine canyons. Due to the limitations of laboratory scaling, however, the flow field produced in the laboratory is similar but not identical to oceanographic flows (Allen et al, submitted). Field observations over canyons give sparse data that are both costly and challenging to collect. As a result, advancements in understanding of currents in these regions relies 2 Regions of the coast are, at times, limited by the micronutrient iron (Hutchins and Bruland, 1998) Chapter 1. 12 Introduction heavily on both laboratory and numerical predictions. Developing a better understanding of flow through submarine canyons has been difficult due to the nonlinear nature of the flow and the steep topography of these features. Comparisons between laboratory and numerical model outputs have revealed that the two do not give identical results (Allen et al., in prep.). While laboratory models can be expected to truthfully obey the physics of a problem, quantifying the results is not always possible. Numerical models, on the other hand, have the opposite relation. The results that they give can be easily quantified, but some of the physics must be removed. In this study the goal has been to develop a method of quantifiably measuring the upwelling in a submarine canyon with a laboratory model. M e t h o d of S o l u t i o n 1.2 The laboratory experiments performed here as well as those studying canyons in the past (Allen et ai., submitted; Perenne et al, 1997) are performed in a cylindrical rotating tank where velocities in the tank are forced by altering the rotation rate. In order to create an upwelling favourable environment, the rotation rate of the tank is increased (causing the currents in the tank to travel in the direction opposite to shelf waves). Quantifying the upwelling flux directly is extremely challenging if not impossible. In this thesis, the upwelling flux through the canyon is measured by examining how the surface currents of the tank evolve. 1.2.1 T h e Effects of R o t a t i o n a n d a C a n y o n o n F l o w E v o l u t i o n A non-rotating current travelling over a flat bottom boundary will have a governing equation of ™ = ^ u, 2 (1.3) Chapter 1. 13 Introduction where u is the velocity and v is the kinematic viscosity. Given a no-slip boundary condition along the bottom boundary, the timescale for velocity evolution is H v~ , 2 l where H is the fluid depth. For a fluid that is 0.1 m deep and has a kinematic viscosity of 10 -6 m s , this timescale is 10 seconds, (or a little over 2 and a half hours). 2 -1 4 In this problem, however, the currents are rotating and the forcing mechanism itself is dependent on increasing the rotation rate. A byproduct of this forcing is the addition of anticyclonic vorticity to the flow field. As a result, an Ekman layer will form along the bottom boundary, which will advect the water radially outwards. This process causes vortex stretching in the interior, which will accelerate the flow evolution. Analysis of the boundary layer response to the induced current makes possible predictions of how the boundary layer will cause the interiorflowto evolve, or spin-up. For aflat-bottomtank, the timescale for flow evolution is iiTQ~ 1//2 i/ / -1 (Greenspan, 1968). For a rotation rate 2 of 0.25 s , and the same depth and viscosity as above, this will result in a timescale of -1 200 seconds. In this study, the model topography replicates the continental shelf and slope. The evolution of currents are observed in the presence and absence of a canyon. Experiments show that the timescale of flow evolution with a canyon present is significantly less than for the canyon-free scenario. Differences in flow evolution can be used to quantifiably measure the force imposed by the canyon on the flow. This force, F, which is used to determine the upwelling through a canyon, is defined as F = pC L L U , (1.4) 2 d 1 2 where p is the density, Cd is a drag constant, Li and L are appropriate lengthscales, 2 and U is the velocity. The appropriate lengthscales for this problem are not examined in this thesis, although suggestions by Allen and Hickey (in prep.) suggest that they are the width of the canyon at the shelfbreak and the radius of curvature of the shelfbreak Chapter 1. 14 Introduction as it turns into the upstream flank of the canyon. It is believed that C is dependent d on the velocity, stratification and rotation rate of the flow approaching the canyon as well as canyon geometry. Although effects of geometry are not studied here, the exact relationship between Cd and the other properties is examined in this thesis. 1.2.2 O r g a n i z a t i o n of T h i s Thesis Chapter 2 will discuss work done in the laboratory. The choice of scales used in the model, the laboratory methods performed for data collection and velocity measurement, and an error analysis are presented. Chapter 3 examines the spin-up behaviour in the absence of a canyon. It begins by giving a brief history of the spin-up problem and introduces the limitations of pre-existing theory to the this study. Extensions of current theory are developed and tested against laboratory data. Chapter 4 examines the spin-up problem when a canyon is introduced to the model. The effects of a canyon are measured and the forces imposed by the canyon in the various trials are determined. Analysis finds a relationship between the force and the flow properties. The thesis closes with a discussion of results and conclusions in Chapter 5. Chapter 2 L a b o r a t o r y M e t h o d s a n d D a t a A n a l y s i s Techniques The goal of laboratory modeling, also known as physical modeling, is to better understand a phenomenon found in nature by replicating it in a controlled laboratory setting. Quite often one is not able to reproduce the field scenario exactly, and the important non-dimensional quantities are matched by scaling the appropriate dimensional quantities instead. Clearly, the flow patterns being modeled in this study must be scaled, as those in the ocean have lengthscales that are O(km). The tank used for these experiments is one meter in diameter and the topography within mimics that of a generic coastal ocean, albeit at a 10-fold vertical exaggeration: an abyssal plain in the tank centre adjoins the continental slope, leading to a continental shelf and ending in a coastline at the tank edge (recall Figure 1.1). The vertical scaling is necessary in order to correctly replicate viscous effects in the laboratory. To study the effects of a submarine canyon, a 22° slice of the continental slope is removed and replaced by a topographic section that mimics the bathymetry of a submarine canyon (Figure 2.1). For some experiments, the topography was radially symmetric, i.e. there was no canyon present in the laboratory set-up. 2.1 Scaling A n a l y s i s The non-dimensional parameters that are important to this problem are the Ekman number, the Rossby and the Burger number. The Ekman number, Ek, is the ratio of the bottom boundary layer thickness, 8, to the water depth. Following Hewett (1998), 15 17 Chapter 2. Laboratory Methods and Data Analysis Techniques the chosen depth scale was the shelfbreak depth, H , giving an Ekman number of s Ek = ~ . tl (2.1) s The Rossby number, Ro, relates the importance of advection to rotation Ro=JL, (2.2) where U is the velocity, 1Z is the radius of curvature of the upstream flank of the canyon and / is the Coriolis frequency. The Burger number, Br, relates the importance of buoyancy to rotation and is equal to where H is the shelfbreak depth, L is the length of the canyon normal to the shelfbreak s c from its head to its mouth, N is the Brunt-Vaisala frequency or buoyancy frequency, g is gravity p is a reference density, p is the density and z is the vertical direction. It 0 was discovered during analysis that the Froude number is also important. The Froude number, Fr, is defined as =m - Fr (2 5) and gives the relative importance of velocity to buoyancy effects. The laboratory model was originally designed to mimic Astoria canyon (Figure 2.2), and the laboratory scaling is therefore matched to the oceanographic conditions found there (Table 2.1). The values N — 2.2s , / = 0.52s -1 -1 and U = 1.16cm s _ 1 were taken as a set of standard values. To see how these values affect the flow, they were changed separately in order to better determine how the parameters affect flow evolution. The trials performed are shown in Table 2.2. 18 Chapter 2. Laboratory Methods and Data Analysis Techniques 46°25 1 46°20' 46°1? 46 10' o 124°40' 124°35' 124°30' 124°25' 124°20' 124-15" 124°10' 124°05' 124°00' 123°55' Figure 2.2: Bathymetry of Astoria canyon. Reproduced with permission from Hickey (1997) 19 Chapter 2. Laboratory Methods and Data Analysis Techniques Parameter Astoria Lab L H 22 km* 150 m 4.5 km* 50 m^ 10-2 s- 20 cm s0.44 8 cm 2 cm 4 cm 2 mm 2 s0.5 s" 1.2 cm s 0.60 0.3 0.13 0.13 0.1 0.3 0.25 c n s f s N f U R 10-V 1 xt 1 1 1 Ek Bu Fr - 1 Table 2.1: Scaling of parameters used in laboratory modelling. The blank line marks a separation between dimensional and non-dimensional values t denotes values are from Allen and Hickey (in prep.); J denotes values are from Hewett (1998) Trial 1 2t 3* 4 5 6 7 8 / N 2.2 SB 0.40 0.70 A/ 0.12 0.12 0.12 0.12 0.12 0.12 0.52 0.52 0.09 0.06 0.52 0.52 0.52 0.52 0.0 1.1 4.4 2.2 2.2 2.2 2.2 Table 2.2: Trials performed in the laboratory. Each trial was performed with and without a canyon present in the tank topography. U = r A / is equal to the initial geostrophic current at the shelf-break, where r is the radial location of the shelf-break in the tank. Units for each of the dimensional parameters are N = [s ], / = [s ], Af = [s~ ] and U = [m-s ]. Values of i V , / and Af that have been changed from the standard values are accented in bold, t Trial 2 uses a homogeneous fluid; * Results for trial 3 were lost and could not be analyzed. SB SB SB -1 SB -1 SB S B -1 : Chapter 2. Laboratory Methods and Data Analysis Techniques 2.2 20 E x p e r i m e n t a l Trials In order to avoid mixing and to achieve solid body rotation more quickly, the tank was filled with stratified water while rotating over a period of approximately 90 minutes, and was left to spin-up for an additional two hours, ensuring that it had achieved solid-body rotation. Stratification was achieved by adding salt, hence density diffusion is assumed negligible over the course of an experiment. Stratification was created using the Osier method (Osier, 1965), where two containers are filled with equal volumes of fresh and salt water, respectively. These containers are connected at the bottom by a piece of tubing, and there is a mixer in the salt water container. Another piece of tubing attaches the container housing the saline water to the tank. If the Olser method is used to fill a tank containing only vertical walls, then it will result in a linear stratification. Due to the sloping nature of the topography in our tank, however, the stratification is stronger near the surface (Figure 2.3). Velocity forcing in the tank was achieved by varying the rotation rate of the tank, which was accelerated from an initial rotation rate to a final rotation rate over one rotation period of the tank. This impulse is sufficiently long to allow development of an Ekman layer along the bottom, but short enough to leave the interior rotating at approximately its original rate. Pliolite V T - L particles were used as tracers and were placed on the fluid surface. (These particles are produced by Goodyear Chemical Company.) In order to reduce surface tension, the water which was used to fill the tank was treated with small amounts of dish soap and Photo-Flo (produced by Kodak). A videocamera mounted in the rotating reference frame recorded the flow and velocities were calculated using MatPIV, a particle image velocimetry (PIV) package for Matlab (Sveen, 2000). This technique provides a continuous time-series of velocity data, allowing us to observe the evolution of the flow field at the tank surface. Pathlines were created by overlaying Chapter 2. Laboratory Methods and Data Analysis Techniques 21 Mixer 5Z7- Fresh Water Salt Water To Tank 1 m Figure 2.3: Creation of a density gradient in the laboratory a) Schematic showing method for creating a density stratified fluid to fill the tank, b) Cross-sectional view of the tank bathymetry. The location of the shelfbreak is highlighted by the cross-hairs. The tank was filled from bottom to top, with the stratified water coming from the buckets shown in a). This results in the density profile shown in c) and the buoyancy frequency, N, shown in d). Values shown in c) and d) are theoretical expectations giving a stratification representative of Trials 1, 7 and 8. The dotted line represents the location of the shelfbreak Chapter 2. Laboratory Methods and Data Analysis Techniques 22 images from the video camera, providing a tool for qualitative flow visualization. 2.3 P I V techniques P I V is a very powerful correlation tool for flow-field analysis as it does not suffer from complications arising from tracking individual particles. P I V traces particle distributions in a flow field, rather than the particles themselves. The fundamental concept behind this method is that a digital image can be considered as a signal field analagous to a digital time-series. In order to obtain velocity measurements, two sequential images are subsampled in small regions known as "interrogation windows", and the displacement of particles between two video frames within a given interrogation window is determined by cross-correlation methods (Figure 2.4). The cross-correlation function, Cr, is defined as Cr(x,y)= M/2 N/2 Y, J2 i=-M/2 where Im\ and lm 2 Im (i,j)Im (i 1 + x,j+ y) 2 (2.6) j=-N/2 are the interrogation windows being compared and, usually N — M. Calculating Cr using the method in (2.6) will require 0(N ) computations. This number A can be reduced by applying the correlation theorem: correlation in the spatial domain is equivalent to multiplication in the frequency domain. This theorem allows (2.6) to be rewritten as Cr = &{ifft [fft(Jmi) • fft(/m )* ]} r 2 (2.7) where fft is the Fast Fourier Transform, ifft is the Inverse Fast Fourier Transform and I* r is the reverse conjugate matrix of / . (For a matrix A m,n m,n-l a A* = r nm,n m x n , a:m—l,n m—l.n—1 a 1,71-1 n m—l,n a l.n 0 , 1 (2.8) Chapter 2. Laboratory Methods and Data Analysis Techniques 24 where the * denotes the complex conjugate.) Using this method decreases the computational expense from 0(N ) 4 to 0(N logN) computations. 2 If there is a correlation between the two images, this will be expressed by a noticeable peak in the function Cr. The position of this peak relative to the center of the first interrogation window gives the displacement between the two images. It is assumed that the particle displacement between two frames is homogeneous within a given interrogation window. 2.3.1 Error Analysis and Propagation A balance between efficiency of experimental procedures and minimization of error must be found. The accuracy of measurements attained by P I V techniques is a combination of several aspects in the procedure beginning with recording techniques and ending with evaluation methods. In order to determine the error in P I V measurements, Raffel et al. (1998) create an artificial set of images with known particle displacements numerically. They then analyze the artificial images using P I V in order to empirically determine the error in measurements. This section presents a breakdown of their findings, as relevant to this study. • Particle Image Diameter (ej ) D Raffel et al. (1998) show that for cross-correlation between two images, particle sizes slightly larger than 2 pixels is ideal. For the window size and resolution used in the experiments, an ideal particle size would be approximately 1.5 mm. While pliolite particles are definitely smaller than this, they have a tendency to coagulate. Examination of sections of three random images suggests that the average particle size is near 4 pixels in diameter (Fig. 2.5). The interrogation windows used were 32x32, giving the random error due to particle image diameter as approximately 25 Chapter 2. Laboratory Methods and Data Analysis Techniques 0.015 pixels. A bias error can be introduced when particles are small. Sub-pixel locations are determined using a three-point estimator. For particles smaller than one pixel, this bias can be significant although this error quickly diminishes for larger particles. For particles greater than two pixels, the error is less than 0(0.001 px) and can be ignored. • Particle Image Shift (e/5) PIV uses the Fourier transform which inherently assumes that windows are periodic. If displacements are larger than half the window size, aliasing will occur. Introducing a displacement between the two interrogation windows improves accuracy of P I V measurements. A shift was assumed to compensate for average velocities in the tank (0.5 cm-s ), but this will still allow errors in velocity calculations. Veloc-1 ities in the tank during P I V evaluations had a maximum value of 1.8 c m - s . This -1 is equal to less than 4.5 pixels per frame. For a Particle Image Diameter of 4 pixels, and a 32x32 interrogation window, a random error of 0.05 pixels is predicted. By shifting the second interrogation window, this error is reduced to 0.045 pixels. Late in time, when velocities are near zero, the shift described above gives a predicted error of 0.04 pixels. It will be assumed hereafter that the error is consistently 0.045 pixels. • Image Quantization (CQL) Video data were stored in 8-bit resolution. Willert (1996) and RafTei et al. (1998) both found that increasing quantization levels has very little effect on the uncertainty of measurements. In fact, decreasing the quantization levels from 8 bits per 27 Chapter 2. Laboratory Methods and Data Analysis Techniques pixel down to 4 has practically no influence on the anticipated error. For a particle density of 10.2 particles in a 32x32 interrogation window, which is significantly lower than the values from these experiments, the uncertainty is expected to be 0.025 pixels. • Background Noise For image quantization of 4 bits per pixel or greater, RMS-error is essentially unaffected by noise (Raffel et ah, 1998). • Displacement Gradients (CDG) Experimental results show that the steepest gradients in displacement will occur across the shelf-break late in experiments with strong stratification. Displacement gradients in these trials are found to have a maximum value of 0.07 (pixels per pixel). Raffel et al. (1998) find that for experiments similar to the ones presented here, this will result in an uncertainty of 0.15 pixels. It is crucial to note that all estimates for error above err on the side of caution; the error predicted above is representative of the flow field at its most troublesome velocities. The total error in calculations can be expected to be: (2.9) « V0.015 + 0.045 + 0.025 + 0.15 px 2 2 2 2 0.16px Clearly, (.DG dominates the R M S error over the shelf-break and continental slope. Over the continental shelf, CDG ~ 0.08 px, which results in e t ~ 0.10 px. This leads to to velocity uncertainties of about 5.6 m m - s the continental shelf. -1 on the continental slope and 3.5 m.m-s -1 over Chapter 3 Laboratory Spin-up Without a Canyon This thesis examines subtle differences in flow evolution with and without a submarine canyon present in the. laboratory topography to determine quantitative expectations of the upwelling behaviour of water passing through a submarine canyon. Understanding the governing mechanisms that describe flow evolution in the tank without a canyon is critical in being able to measure the effect of a canyon on the evolution of the flow. The flow in the laboratory is forced by increasing the rotation rate of the tank. The process by which the flow evolves, which ends when the fluid is uniformly rotating at the increased rotation rate of the tank, is called spin-up. A n example of the spin-up process in the laboratory is shown in Figure 3.1. This chapter develops the necessary theory to be able to successfully the predict spin-up processes when no canyon is present in the laboratory topography. 3.1 Introduction to the Spin-up Problem The upwelling problem described in Chapter 1 and the spin-up problem are very similar in nature as they both describe a system in which ageostrophic events create a secondary current that strongly affects an otherwise geostrophic flow. In the classic spinup experiment, a cylindrical, rigid-walled tank containing a viscous fluid in solid body rotation has its angular velocity increased impulsively by a small amount. Initially the fluid interior, which is away from any boundaries, is unaffected. Any walls not parallel to the axis of rotation, however, will develop Ekman boundary layers. In these layers, the 28 30 Chapter 3. Laboratory Spin-up Without a Canyon /=2Q A o 4 5 Figure 3.2: Cartoon showing the secondary circulation of a homogeneous fluid in a cylinder. The three-sided box (thick line) represents a vertical cross-section through a cylindrical tank holding water. The lowest order circulation of the flow, indicated by the vectors above the water, is anticyclonic, causing a divergence in the bottom Ekman layer (8, below the dotted line). The Ekman divergence flow is fed by Ekman suction out of the interior, which adds cyclonic vorticity to the interior through vortex tube stretching. Water exits the Ekman layer at the tank edge where it enters the side-wall boundary layers (K, outside the dashed lines). The side boundaries act as a replenishing source for water within the interior, which has a weak, radially inward flow. was used to explain the rapid spin-up of rotating flows both theoretically and in the laboratory (Greenspan and Howard, 1963; Greenspan, 1964; Greenspan, 1965; Greenspan, 1968). Other works went on to examine spin-up and spin-down in a homogeneous fluid where the Rossby number varied from -1 to 1 (e.g. Wedemeyer, 1964; Benton, 1973). These works were mathematical models, in that the flow was laminar and did not include instabilities characteristic of large Rossby numbers. The spin-up behaviour for larger Chapter 3. Laboratory Spin-up Without a Canyon 31 Rossby numbers is more complicated and hindered the development of a more general theory until more recently when more complex mathematical models (O'Donnell and Linden, 1991) were developed and non-linearities were further examined (Konijnenberg et ai, 1995). Other studies (e.g. Linden and van Heijst, 1984; Konijnenberg et al., 1996) have looked at spin-up in non-axisymmetric containers, but the question of how variable topography - axisymmetric or non - might affect spin-up has received little attention. 3.1.1 The Spin-up of a Homogeneous Fluid Over a Flat Bottom In homogeneous fluids, spin-up proceeds as described above (Figure 3.2). Greenspan (1963; 1968) examined the classic spin-up problem for cases where the increase in rotation rate, Afl, is considerably less than the initial rotation rate, fl, i.e. AQ/Q < 1. A homogeneous fluid contained by a circularly cylindrical, flat-bottomed container has three important parameters that will determine its spin-up rate: a temporal Rossby number, Ro = AQ/T2, which gives the ratio of the change in angular velocity to the final rate of t rotation; a Froude number, F = 4Q R /(gH), 2 2 where R is the radius of the tank, g is the gravitational acceleration, and H is the depth of the tank , which measures the relative 1 importance of rotation to gravity; and an Ekman number, E = v/(QH ), 2 where v is the liquid's kinematic viscosity, which represents the importance of viscosity to rotation (Greenspan, 1968). The spin-up of a fluid with a velocity described by u^,(r) = AQr, where AQ <C 0 , will obey = AQrexp(-t/V), where equal to (3.1) is the tangential velocity, r is the radial coordinate and r is a decay timescale 1/(QE~2). Adding a non-homogeneous fluid or a more complicated tank shape causes the spin-up This is a curious definition for the Froude number as it seems to resemble the Burger number. It would make more sense for it to be defined either as the Burger number or as 4AJ7 B?/(gH). 1 2 32 Chapter 3. Laboratory Spin-up Without a Canyon problem to become considerably more complex. 3.1.2 T h e S p i n - u p of a Stratified F l u i d s O v e r a F l a t B o t t o m Convection-style spin-up fails for non-homogeneous fluids as vertical motions are weakened, and eventually decay, due to the vertical structure of the fluid. In the late 1960's and early 1970's this problem began to receive much attention (Benton and Clark, 1974). In addition to the temporal Rossby number, the Ekman number and the Froude number, the Prandtl number and the Burger number become important parameters for stratified spin-up. The Prandtl number measures the relative importance of the diffusion of momentum against the diffusion of density and is defined as P = u/x, where v is the kinematic viscosity and x is the density diffusivity. The Burger number, S = (N/f) , 2 where N is the buoyancy frequency, p is density, z is the vertical coordinate, and / is the Coriolis parameter, measures the relative importance of rotational and gravitational effects on the fluid. Laboratory spin-up experiments of a thermally-stratified fluid in a cylinder showed that Ekman layers of thickness 0(E l ) l 2 still formed, but that Ekman-type effects die out before solid- body rotation is achieved, leaving the fluid in a quasi-steady state with the following characteristics (Holton, 1965): • close to the horizontal boundaries, the interior rotates at the same rate as the tank that contains it; • the angular velocity increases on a vertical length scale L/S / 1 2 3> S, where L is the horizontal lengthscale of the tank and 5 is the Ekman layer depth. 33 Chapter 3. Laboratory Spin-up Without a Canyon As there is no viscous stress within the Ekman layer, no Ekman pumping occurs and convection-like spin-up is halted, leaving diffusion of momentum as the mechanism for restoring solid-body rotation. Further mathematical examination of the problem brought the relevance of these results into question. The result presented by Holton (1965) was determined with noninsulating side-boundaries. For an insulated fluid, buoyancy forces result in side walls being unable to accept the mass flux out of the Ekman layer due to gravitational effects of the accumulation of dense water along bottom edge of the tank (Pedlosky, 1967). Pedlosky (1967) concluded that as a result, spin-up by non-diffusive means would be halted prematurely and that the quasi-steady solution could not be achieved. These conclusions turned out to be incorrect. Despite this behaviour by the side walls, a secondary circulation does form (Sakurai, 1969; Waiin, 1969) and the time required for the realization of the quasi-steady state described by Holton (1965) is 0(L/(VluS) ) 1/2 3.1.3 (Waiin, 1969). S l o p i n g E k m a n Layers Ekman layers on a shallow slope have also received considerable attention in the more recent past, mostly driven by the desire to better understand mixing processes invoked by the bottom boundary layer of the ocean. The behaviour of an Ekman boundary on a slope in a homogeneous fluid generates the familiar Ekman spiral (Pedlosky, 1987). The addition of stratification to flow over a slope makes the problem considerably more difficult. Initial attempts to solve the problem specified density at the boundary (Holton, 1967; Hsueh, 1969), something that is more appropriate for atmospheric flows. Follow-up studies provided an insulated boundary (Weatherly and Marin, 1978; Thorpe, 1987) but the solutions gave curious results. For example, infinite transport was predicted for flat slopes. When the problem was adjusted and solved for an oceanographically relevant case, Chapter 3. Laboratory Spin-up Without a Canyon 34 results showed that Ekman transport upslope causes the advection of dense water into a region of less dense water, imposing a downslope gravitational force on the fluid. Eventually the gravitational force is as strong as the rotational forces and the boundary layer arrives in a thermal-wind balance such that there is no cross-slope transport. "Fast spinup" processes with an Ekman-related solution are halted prematurely. Furthermore, the stagnant, or shutdown, boundary layer, acts as a relatively slippery boundary between the interior flow and its more quickly rotating container, thereby decreasing the rate of momentum diffusion to the interior (MacReady and Rhines, 1991; 1993). 3.2 Objectives a n d O u t l i n e of T h i s C h a p t e r The aim of this chapter is to determine the important spin-up processes in the laboratory when a canyon is not present in the tank topography. Since the laboratory work in this thesis involves an evolving flow over a steep slope, the spin-up behaviour of a fluid over a steep slope must be understood. Although the behaviour of an Ekman layer on a slope has become an area of increasing interest, to date studies have only examined Ekman suction on shallow slopes. This chapter examines the limitations of linear theory for homogeneous Ekman layers over a steep slope, then applies this theory to determine how suction by a sloping Ekman layer will affect laboratory spin-up. Also, spin-up theory for a stratified fluid over a shallow slope is tested to see if its predictions can be extended to steeper slopes. Findings are compared with laboratory results to determine if the theory accurately predicts the spin-up behaviour witnessed in the tank. Results from this chapter will be applied in Chapter 4 to measure the effects of a canyon on flow evolution. 35 Chapter 3. Laboratory Spin-up Without a Canyon 3.3 Spin-up of a Homogeneous Fluid Over a Steep Slope In the absence of a canyon and stratification, the spin-up process is governed almost purely by Ekman layer dynamics. Radially outward transport is generated in the bottom boundary layer. Continuity dictates that this flow must be balanced by a flux from the interior into the boundary layer and a weak, radially inward flux in the fluid interior (Figure 3.3). As this secondary flow is driven by the boundary layer response to bot- O interior flow Ekman transport Ekman suction Return flow top of boundary layer Y 8 H L f — slope angle — Ekman layer depth — depth scale length scale — Coriolis frequency — H /, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \^ \ \ \ \ \ \ \ \ \ \ \ \Y \\\\\\\\\\\\\\\\\N X Figure 3.3: Cartoon showing the secondary circulation of a homogeneous fluid in the laboratory. torn friction, (in this case, over a slope,) analysis begins by examining the momentum equations for a rotating bottom boundary layer on a slope. 36 Chapter 3. Laboratory Spin-up Without a Canyon 3.3.1 B o u n d a r y Layer F l o w If the fluid is assumed to be incompressible then the equations of motion and mass conservation in a rotating reference frame are ^ Dt + fkxu = - ^ + vV u, p (3.2) 2 V - u = 0. (3.3) In these equations u is the velocity vector (u,v,w) in Cartesian coordinates (x,y,z), f is the Coriolis frequency which is equal to 20, p is the perturbation pressure equal to the total pressure less the undisturbed hydrostatic pressure pgz, p is the fluid density, g is gravity, V is the gradient vector (d/dx,d/dy,d/dz), D/Dt is the total derivative d/dt + u • V , k is the unit vector (0,0,1) and v is the kinematic viscosity of water. Although the problem here involves spin-up over a slope that expresses radial symmetry, an assumption is made that the ratio of the boundary layer height to the radius of curvature is small 0 ( 1 O ) , allowing the boundary layer behaviour to be examined in a -2 Cartesian coordinate system. It is assumed that, at least at lowest order, the flow in the tank, like the topography, is radially symmetric. This allows us to assume d/dy terms to be small, if not zero. We also assume that in the boundary layer d/dx <C d/dz. To study the boundary layer on a slope, it is useful to introduce a rotation of coordinates about the y-axis through the slope angle (Figure 3.4), 7, such that x' = x cos 7 + z sin 7, x = x' cos 7 z' = z cos 7 — x sin 7, z = z' cos 7 + x' sin 7, u' = u cos 7 + w sin 7, u = v! cos 7 — — z' sin 7, w' sin 7, w = w cos 7 — u sin 7> w = w'cos^y + u'sin^. 1 38 Chapter 3. Laboratory Spin-up Without a Canyon Dw , . 1 dp „ — + / w s i n 7 = — ^f + u^w. ,„ „. (3.8 9 The viscous terms in the above equations act as a singular perturbation that is negligible in the fluid interior, but gives rise to the Ekman boundary layers where it is of the same importance as rotation. In the boundary layer, one can neglect the time derivative of the above equations by assuming that the boundary flow remains in a quasi-steady state with the interior flow, which evolves over a larger timescale. Neglecting the time-derivatives results in du ™--/,cos 1 dp du = - - _ + „ _ , 2 7 dv oz (3.9) dv v——, dz 2 w——h / (ifCOS7 ldw - — 2 — iusin7) = 1 I dp dw + /,sm = - - - + , — . (3.10) 2 (3.11) 7 The continuity equation is du dx dv dy dw _ dz . ^\ It is useful to introduce the scaling, z = 5(*, w = 6(V/L)w*, u = Vu*, v = Vv* and x = Lx*, where V = A£IL is the geostrophically balanced interior flow, L is the length scale, and 5 = y^2v/(f cos 7) is the Ekman layer depth. The asterisks denote nondimensional variables. The momentum and continuity equations are further treated with the non-dimensional parameters, R = V/(fL cos 7) and tan 9 — 5/L, where R is a Rossby number and (tanf?) is a ratio of the length-to-depth scales used in the problem. This allows us to scale the pressure field as p = pVLf cos'jp*. Incorporating these scalings, the momentum equations become (dropping the asterisks) du Rw— -v a; D _ dv dp 1 1 du = - — + — , dx 2 cos 7 a ; 2 . 1 1 Rw77- -\- u — VJ tan 7 tan 0 = — d( (3.13) dv /o 1 ^ 2 —-, 2 cos 7 dC 2 (^-14) 39 Chapter 3. Laboratory Spin-up Without a Canyon R tan 9 dw .dp ,d w — o ^r + vtawy = cot 6— + tan 9——. 2 d( o( aC 2 2 (3.15) Further analysis is aided by introducing the asymptotic expansion, oo (u,v,w) = R {u ,v ,w ). n n For the laboratory setup, L = 0.5m, / = 0.52s tanf? = 0(R ), 4 -1 n n and v — 1 0 m s . These values give - 6 2 - 1 allowing us to ignore contributions of this order. The results of (3.15) allow us to determine the vertical variation of the pressure field, which is not of interest for this problem and the equation is dropped from the analysis. Assuming that dp/dx is balanced by the geostrophic velocity in the interior, and implementing the asymptotic expansion gives, to first order, 1- v 1 2 u (3.16) 2 2 u 1 d cos 7 dC v which has a classic Ekman spiral solution u = A(x)e~ ^ V2 i — sin (^/2 cos 7 C ) (3.17) 1 — cos ^ 2 cos 7<^ v where A(x) is derived from the geostrophic velocity of the interior and dA/dx is small. Continuity is used to determine the flux out of the boundary layer: °° du dC, _ — f°° r w Jo dx 1 2y/2cos7 dA dx' ^ , 1 8 ^ If one examines (3.18), it appears that w ~ oo as 7 —> 90°. This finding is a result of the scaling of w, which also is also dependent on (COS7) / , - 1 2 and signals that at some stage the choice of scales for this problem breaks down. When this breakdown occurs, the results from this solution may no longer be correct. It is necessary, then, to determine 1. whether or not the solution presented is valid for the boundary layer present in the laboratory; and 40 Chapter 3. Laboratory Spin-up Without a Canyon 2. when the solution begins to break down. As the slope increases, so do the boundary layer thickness (S = ^2v/(f cos 7)) and the Rossby number (R = £ / / ( £ / c o s 7)), which affect the size of tan 9 and the validity of the asymptotic expansion, respectively. On the continental slope, (7 = 56° and L = 0.25 m. For V = 0.01 m s , / = 0.52 _1 s _1 and v = 10~6 m s ,) tan0 = 0.0026 and R = 0.3. The Rossby number is becoming 2 -1 large and the danger of nonlinearities due to the importance of advective terms looms. Numerical experiments of spin-up over a flat bottom have shown that linear theory holds for Rossby numbers smaller than about 0.5 (Hart, 2000; van de Konijnenberg and van Heist, 1995). With this information it is reasonable to assume that the behaviour of an Ekman layer on a slope should follow linear theory when R < 0.5, so long as tanf? remains small (I propose an arbitrary limit of tanf? < 0.05). For the conditions present in the laboratory, the above solution should hold until the slope reaches about 80° at which point the Rossby number approaches 0.45 and tan 0 approaches 0.02. This suggests that not only is the solution useful for solving the problem relevant to the laboratory set-up, but that it is well within the domain where its result can be trusted. With this information one can proceed, having faith that the flux into the boundary layer has been properly characterized. It is now necessary to determine the relation between the right side of (3.18), dA/dx, and the interior flow. This will be done by developing a physical explanation of A, determining its value in the laboratory, and implementing this value into (3.18). In the derivation that has brought us to this point, it has been assumed that although the tank is cylindrical, the boundary layer thickness is small enough that we can ignore any curvature effects and solve the problem in Cartesian coordinates. While defining A, Cartesian coordinates will remain in use even though the following explanation takes 41 Chapter 3. Laboratory Spin-up Without a Canyon place in the interior. Once a physical description of A is made, its value with respect to the cylindrical coordinate system is easily defined. To understand the physical meaning of A, begin by considering an infinitely long, sloping boundary of angle 7 below a half-space of homogeneous fluid on an /-plane. The fluid has an along slope velocity, v(x), with a horizontal shear (Figure 3.5a). If the velocity had a cross slope component to it as well, the curl of this velocity field normal to the boundary would be dv/dx' — du/dy but as there is no u-component here, it is simply dv/dx'. If the fluid is assumed to be inviscid except at the bottom where it interacts with the sloping boundary, then the boundary layer will respond to the component of the curl whose direction is normal to its plane, i.e. (0, 0, dv/dx'), but not those which lie in its plane, i.e. (—dv/dz',0,0), (Figure 3.5b). This relationship implies that as 7 increases, the amount of shear that contributes to the boundary layer evolution will decrease as cos 7. Returning to the laboratory setting, imagine that v(r) is the along-slope velocity in the interior, v can be non-dimensionalized as VA(r), where V is a characteristic velocity of the interior, and A(r) is the non-dimensional velocity present in (3.18). This definition lets us rewrite (3.18) as w = —. cos 7V x A(r), 2^2 cos 7 1 (3.19) or, in dimensional coordinates, w = \J^V xv { (3.20) where v is the velocity in the interior. l This result states that the slope angle has no effect on the flux out of the boundary layer, which may be surprising. However, if one considers that the total frictional drag imposed by the boundary layer over a given horizontal section will be increased by a factor of sec 7 but that the effective value of the Coriolis frequency in the boundary layer 43 Chapter 3. Laboratory Spin-up Without a Canyon will be decreased by a factor of cos 7, then the result is less worrisome. 3.3.2 Effect of a Slope on E k m a n S u c t i o n For small 7, it is often assumed that sin. y ~ 0, cos 7 « 1, and the vertical component / of Ekman suction is equal the total Ekman suction calculated over a given region. For steeper slopes, these assumptions do not hold. As the slope increases, the vertical flux into the boundary layer decreases as cos 7 (Figure 3.6). When sin 7 is large, Ekman suction will invoke a radially outward flow in the interior. For these laboratory experiments, it is assumed that horizontal flows due to Ekman suction are small and can be ignored for the desired accuracy of solution. The expected impact of a sloping boundary layer on the the vertical flux from the interior is shown in Figure 3.7. This flux, which is a function w sin Y X cosY Figure 3.6: Vertical velocity contribution to Ekman suction on a slope. of the interior flow, is responsible for the non-canyon and non-diffusive spin-up processes. 44 Chapter 3. Laboratory Spin-up Without a Canyon 0.1 F 11 i i i i i i i i 1 0 10 20 30 40 50 60 70 80 90 Slope A n g l e Figure 3.7: Vertical component of Ekman pumping into a sloped boundary layer (i.e. vertical flow into the Ekman layer, not the flow normal to the boundary layer) is a function of slope angle and shear. Vertical velocities are exaggerated by a factor of 1 0 . The angle to which this relationship remains valid is dependent on the geometry of the problem. For the laboratory set-up, the solution is expected to hold until 7 reaches an angle of about 80°. -3 45 Chapter 3. Laboratory Spin-up Without a Canyon 3.3.3 Effect of E k m a n suction o n S p i n - u p Once the flux into the boundary layer is determined, the next step is to understand how this flux will affect the evolution of flow in the tank. In the interior of the fluid, friction is negligible and at first order potential vorticity is conserved, i.e., where H is the water depth, and u# corresponds to the horizontal velocity field. The superscript, i denotes that these are values for the interior. (3.21) can be expanded to yield ^ (/ + + (/ + * „ ) £ £ = 0. (3.22) Although it may seem that H is a constant, it should be remembered that if there is a vertical flux out of the water column, this is equivalent to a change in depth. Hence, dH/dt — w . Assuming that / is a constant and that nonlinear terms can be ignored, l then (3.22) becomes x if H = ^ ( / + V x v? ) , (3.23) H where / v cos 7 ( V x 4) 4/ 2 w =w cos 7 , (3.24) V x v* = - | - M r or (3.25) H and v is the azimuthal (along slope) geostrophic velocity. (Figure 3.6). In the classic spin-up experiments with a flat bottom, V x t ? H decays exponentially in time, with no added spatial variation (i.e. V x v} = [V x t? ] exp(—t/r), where r is the H H 0 appropriate timescale). Due to the depth variation present in the laboratory set-up of these experiments, it is difficult give a classical spin-up timescale. During the time over which the ramping of the table takes place alone, the strong depth changes within the 46 Chapter 3. Laboratory Spin-up Without a Canyon tank cause the velocity structure to deviate from a linear dependence on the radius like that seen in the classic spin-up experiments (Figure 3.8). Spin-up over a slope cannot be easily solved by analytical techniques due to spatial variation in topography, and the problem is solved numerically. T 20 25 ; 1 30 1 1 r 35 40 45 50 Radius [cm] Figure 3.8: Example of the surface velocity in the tank sampled immediately after ramping of the rotation rate of the table. The dotted line at r = 28.2 cm represents the location of the shelfbreak. The tank edge is at r = 50 cm. In the classic spin-up problem, the azimuthal velocity follows u = ar, where a is dependent only on time. 3.3.4 S o l v i n g the Homogeneous P r o b l e m N u m e r i c a l l y Laboratory results suggest that (3.23) is not separable, and the equations are solved numerically. Determining how the velocity field will evolve over a given time requires three steps that must be iterated through the appropriate number of time steps: 1. the horizontal structure of the vertical velocity field, defined by (3.24), must be solved; 2. how the vertical velocity (and diffusion for stratified trials) affects the evolution of vorticity, defined by (3.23), must be determined; and 47 Chapter 3. Laboratory Spin-up Without a Canyon 3. the velocity must be found by integrating its curl. (3.24) and (3.23) solve the homogeneous system. The solution to the above problem is solved over the radial domain r = [25 cm, 45 cm], with a spatial resolution A r = 2 mm and a temporal resolution At = 0.1 s. The initial condition required for the problem is the velocity in the interior, which is set based on laboratory observations. The boundary conditions are no-slip and free-slip at the innerand outer-edges of the radial domain, respectively. Topography in the numerical model is relatively simple: the continental slope is inclined at 56° and the continental shelf has a grade of approximately 4.5° (Figure 3.9). The method for finding the value of w is straight forward. Discretizing (3.24) in cylindrical coordinates gives w , {n - cos t) - + _ j , (3.26) where the terms here are the same as those defined in (3.24) and the subscripts n represent the discretized locations in space over the r-domain (ri = 25 cm, r = 25 cm + A r , 2 etc.), the subscripts t represent discretized locations in time over the ^-domain (ti = 0 s, t = 0 s 4- At, etc.), at t and 2 t n = 1,2,3,A,...,N = v(r ) at t . V(n,t) n (3.26) gives values for w for t - 1. Once the iu( ,t) is determined, v is calculated a time-step At forward in time. (3.23) n is discretized to give A ( „,\ _ r (n,t) ( 2w H n+l V(n+l,t) f +H \ r n n + + V(n,t) V(n+l,t) i +r ~ V(n,t)\ Ar n T n+1 J + r n 2 (3.27) where A(rv)( n>t) is equal to A(ru)(„ The expression defining A(rv) n i4) = r n + l V( n + h t ) - r V(n, ) after the time-step, A i . n t can be rearranged to give the solution for v V(n,t) = («( +M) n+l - A(ur)(„ )) r n )t ^- (3.28) 48 Chapter 3. Laboratory Spin-up Without a Canyon 30 35 Radius [cm] 40 Figure 3.9: Topography present in the laboratory (solid line) and that used to test the theoretical model (dashed line). 45 Chapter 3. Laboratory Spin-up Without a Canyon 49 To match the boundary condition, the integrations are done from high values of n to low ones, i.e. n = N — 1, N — 2,..., 1, using the boundary condition = 0. Theo- retical expectations developed are compared to laboratory results, but the method for determining initial conditions must first be discussed. 3.3.5 Comparing Theoretical Predictions To Laboratory Results In order to compare the theoretical predictions to measurements from the laboratory, it is necessary to match the initial velocity. For the laboratory experiments, time is set to zero at the end of the forcing period. Flow evolution early in the ramping period is difficult to characterize as can be seen from laboratory data (Figure 3.10). Rather than attempting to solve for the velocity evolution over the ramping period, I chose to match the initial conditions to the measured velocity at the end of the forcing period. Unfortunately, the tank was not completely level, and this lead to a periodic oscillation in the flow field. In order to be consistent in measurements, I sampled velocities in the bottom of the first wave-trough after the end of forcing as it passed through the measurement area. Other obvious alternatives include sampling in the crest of the wave or at the node of the wave (when the velocity anomaly was equal to zero). The first alternative was rejected primarily because cross-shelf velocities were present near the edge of the tank during wave crest. The second was rejected because the rate of wave change is fastest at that time, and consistent identification of this location in time was challenging. 3.3.6 F i t t i n g C u r v e s t o the V e l o c i t y D a t a Once a location in time is selected, the velocity curve is fit to all data within the measurement window with the restriction of r > 20 cm. This was done because the numerical solutions to the spin-up equations described in sections §3.3 and §3.4 extended from 50 Chapter 3. Laboratory Spin-up Without a Canyon 1 cm/s v. \ AV; : ~; ^ /• " ^ , \ i Vt y s s ^ - ' - - " - / / / / / / / - - / / / ^ttt- '/////A'.' (d) V / • ' ' : I .•• ' s s / / rr/ - ' -/////// (f) I S ' Figure 3.10: Laboratory measurements showing nonlinear spin-up in early time. Velocity fields shown are (a) 1 seconds, (b) 3 seconds, (c) 5 seconds, (d) 7 seconds, (e) 9 seconds and (f) 11 seconds after ramping begins. The ramping trime is 27.3 seconds. The large dot represents the middle of the tank. The beginning of the continental slope (dotted line, r = 22 cm), the shelfbreak (dashed line, r = 28.2 cm) and r = 45 cm (dashed-dotted line) are shown to aid in spatial orientation of the image. Ideally, the velocity would be near 0 in figure (a). The disturbances in the flow field are likely due to the oscillation present in the tank, deviation from initial solid body rotation or non-linear effects in spin-up. 51 Chapter 3. Laboratory Spin-up Without a Canyon r = 45 cm to r = 25 cm. The data were fit to a fourth order polynomial that was forced to match the boundary condition, v(r = 50cm) = 0, i.e. the curve of best fit is equal to V(r) = d(r - 50) + C ( r - 50) + C ( r - 50) + C ( r - 50) , 2 2 3 3 4 4 (3.29) where the constants C\, C , Cz and C4 are determined by a least-squares approximation. 2 (3.29) can be rewritten in matrix form as VI JV = X R4XNC\ 4X The solution for Ci 4 X is found by matrix manipulation: C = (3.30) (R'R)- R'V L where R' is the transpose of R. Once C is found, a new matrix R, which represents the radial points in the finite difference numerical schemes, is used to determine the initial condition of the numerical solution, V. Although the value of C is found, it leaves an important question unanswered, namely, "what is the error in the approximation?" This question can be answered in several ways. One option is to follow the the error propagation from measurements of V to C in (3.30). But this method gives the errors in the individual constants, rather than our exact uncertainty in the line the data have produced. As a result, an alternative method was used to determine the error. Following a suggestion in §14.5 of Press et al. (1986) it is assumed that the values collected in V are close enough to the real values that the error in V has a normal distribution. For this reason, a bootstrap-like method was used in determining error estimates. A random selection of 100 points were taken from the laboratory data (whose number ranged from 110 to 172 points) and each had a (normally distributed) random error added to it. The value of C is recalculated in 100 iterations as C = i (R'R)- R'{V* L + QW), (3.31) 52 Chapter 3. Laboratory Spin-up Without a Canyon where V is a random set of V and Q\ 1 xN is & normally distributed set of random numbers ranging from -1 to 1 that is reevaluated for each iteration. The values for the constants, C\ are then each used to individually determine a set of vectors V . A value V is determined as 1 -i v=— 100 yv i 100 t with the error being twice the standard deviation from the mean. x 3.3.7 F l o w C o m p a r i s o n s F o r the Homogeneous S p i n - u p M o d e l A comparison of actual flow evolution to its theoretical expectation is shown in Figure 3.11. Laboratory and theoretical results are shown at 5 second intervals. The data represented in the plot are from an angular section of the tank (the same region illustrated in Figure 3.10). 3.4 S p i n - u p of a S t r a t i f i e d F l u i d O v e r a Steep Slope In this section, the theory developed in §3.3 is extended to include stratification. The previous analysis showed that Ekman suction in a homogeneous fluid over a steep slope can be calculated with relative ease. This finding is now employed as an assumption for the stratified fluid case. When stratification is introduced, however, the upslope transport that is present in the sloped Ekman layer redistributes the density field, advecting denser fluid into a regime of less dense fluid. This introduces buoyancy forces, which eventually become as large as Coriolis forces, effectively shutting down the boundary layer and insulating the interior from the bottom. Although an analytical model for this problem cannot be easily attained due to buoyancy effects in the equations of motion on a steep slope, laboratory results are predicted 53 Chapter 3. Laboratory Spin-up Without a Canyon 0.3 0 seconds 5 seconds 10 seconds 15 seconds 20 seconds 25 seconds 30 seconds 35 seconds 40 seconds 45 seconds 0.4 0.3 0.4 Figure 3.11: Spin-up of a homogeneous fluid (/ = 0.52 s—1, Af — 0.12 s—1) without a canyon, comparison of laboratory observations (dots) to theoretical predictions (lines). Error on data and line are only slightly larger than points and line thickness. The location of the shelfbreak is shown as a dotted line at r = 22 cm. Time slices are shown every 5 seconds. One rotation period is 24.2 seconds. 54 Chapter 3. Laboratory Spin-up Without a Canyon by a diffusion model that is based on the work of MacCready and Rhines (1991) and the theory introduced in the previous section. MacCready and Rhines (1991) predict the shutdown time, r, for a bottom Ekman boundary layer over a gentle slope to be T = cos7], where S/[(l + S)S f 2 S = (Nsin7//cos7) is a Burger number and N is the 2 Brunt-Vaisala frequency. When shutdown of the boundary layer is complete, it insulates the interior from the friction effects of the bottom. In this limit there is no Ekman suction and flow evolution will evolve as (after MacReady and Rhines, 1991): dv dv ( . S —+ u— + fu = v [sm 7 + 7 7 r - £ 9 9 c o s dv ^) ^ 2 \ ( + U S [T+~S . S where z is vertical (i.e. not normal to the slope) and V M 2 H , 9 7 + C 0 S 7 \ J _ r V o ' ' /o > represents the horizontal components of the Laplacian operator, and u is the radial velocity. Horizontal diffusion is not included in the model of MacCready and Rhines (1991) but must be considered for this problem. Near the shelfbreak, the depth change is rapid enough that a horizontal shear will evolve in the flow, requiring the consideration of horizontal diffusion for the problem. For a no-slip bottom boundary and only vertical diffusion, this suggests a spinup timescale of 15 minutes for water only 1 cm seaward of the shelfbreak, compared with about 8 minutes if horizontal diffusion is included. For laboratory simulations, the shutdown timescale over the continental shelf ranged from 2 seconds (/ = 0.52s" , N = 4.4s" ) to 20 seconds (/ = 0.52s" , N = 2.2s" ). 1 1 1 1 (For / = 0.52S , one inertial period is equal to 12 seconds.) This shutdown process -1 is a gradual one suggesting that Ekman suction into the boundary layer, and therefore its corresponding spin-up influence, will be important early in time but will decay on the shutdown timescale. Flux out of the boundary layer will then be a slightly altered version of (3.24): WB = ( V x ff ) exp (3.33) B where WB is the vertical Ekman flux out of the boundary layer and v7 is the horizontal B 55 Chapter 3. Laboratory Spin-up Without a Canyon velocity of the interior along the tank bottom. It is assumed that the Ekman suction of the stratified boundary layer is the same as that for the non-stratified Ekman layer, except that it decays exponentially in time on the shutdown timescale. Following (3.23), it is anticipated that the flow evolution due to Ekman effects is (3.34) A similar approach is taken to characterize the evolution of the diffusive behaviour in the fluid. Early in time, the familiar "isotropic diffusion" is present: (3.35) uV v/ . 2 H By the end of boundary layer shutdown, however, the diffusion will follow (3.32). The formation of this shutdown boundary layer is due to the gradual increase in gravitational force imposed on the fluid as it is advected upslope. As such, it can be expected that the diffusion will change from (3.35) to (3.32) over the shutdown timescale: ^ • 2 + S ^ ^ x ^ = ,{exp(-^)v + 2 \ From conservation of mass, d 2 V'H-Uh — 5 . 2 + [l-exp(-^) \ fl d sin 7 + cos J [- g2 7 r r d' dr2 2 + 2 uH (3.36) = —dw/dz which, given uniform stretching through the water column due to Ekman suction is —w (t)/H. Therefore V x (3.36) can be written l as 56 Chapter 3. Laboratory Spin-up Without a Canyon {I-" B) - I - (-;)] exp + (_|_ y2 v 1 ) (11 + sin27 + cosS exp *jjj d 2 [( V x f l sin27+ TTs ^) dz . ^ = cos ( / + V x 2 tf,) . (3.37) The theory of MacCready and Rhines (1991) being applied to model spin-up is derived to work only on gentle slopes. Due to the lack of a theory describing arrested or shutdown boundary layers on a steep slope, it is assumed here that the general theory of MacCready and Rhines (1991) adequately characterizes the behaviour of the shutdown Ekman layer over a steep slope. Applying this theory to the steep continental slope gives a shutdown timescales of O ( 1 0 s ) . This short timescale suggests that even if shutdown on steep -2 slopes occurs considerably more slowly than on gentle slopes, then the boundary layer will have shut down by the time the ramping of the tank velocity is complete. As in the homogeneous case, the initial conditions are set from the laboratory flow. This poses a slight problem: the flow below the surface is unknown as its vertical structure varies in time. As a result certain assumptions must be made to set the boundary conditions. The first is that effects from Ekman suction are vertically homogeneous. While stratification does have an effect on varying the Ekman flux, it is small over the depths being considered in this study. The second assumption is that over the time of the experiments, diffusive effects from the shutdown boundary layer barely alter the flow. These assumptions were tested by examining the vertical evolution of the flow with time in the numerical model. (Originally vertical variation in Ekman suction was considered, but its impact on flow evolution was negligible.) When model predictions are compared with actual laboratory output (see §3.3.5), the agreement suggests that these assumptions are acceptable. Equation (3.37) is solved numerically with boundary conditions of no-slip at the outer 57 Chapter 3. Laboratory Spin-up Without a Canyon radial edge of the modeling domain and along the bottom, and free-slip along the top and along the inner radial edge of the modeling domain. 3.4.1 S o l v i n g t h e Stratified P r o b l e m N u m e r i c a l l y The method for determining the result to (3.37) is indeed similar to that used for (3.23) in §3.3 except that for the stratified solution vertical variation and diffusion are present. This alters the method of solution to 1. determine the vertical velocity; 2. use Ekman suction and diffusion to evaluate the time-evolution of the velocity curl; and 3. calculate the velocity from the velocity curl. The problem is solved over the radial domain r = [25 cm, 49 cm], with a spacing, Ar = (49 — 25)/99 cm, the vertical domain z = [0 cm, 6 cm], with a spacing, Az = Ar, where z = 0 is at the surface, and the temporal domain t = [0 s, T], with a spacing A* = 0.0625s. Vertical velocity is determined by (3.33) using the same discretized form as (3.26): _ „ W [ n ' t ] ^ n C Q O / (V(H n+l,t) v S + V(H ,n,t) n+u 7 l 4 / l r n+1 V(H ,n+l,t) n + r n ~ n+1 + V{H ,n,t) \ ( 2t n Ar ) * [ r ^ (3.38) where the subscripts n and t are the same as §3.3.4, the subscript H represents the n bottom-most vertical layer within the model, i.e. the layer below-which the coordinates have left the fluid and are beyond the boundary (Fig. 3.12), and r i n+ timescale at r i, where r i = (r n+ n+ n+1 is the shutdown +r )/2. From hence forth, values at the discretized n step n+ \ are equal to the average of the two bounding values at n and n + l. The Ekman layer is not explicitly present in the model. 59 Chapter 3. Laboratory Spin-up Without a Canyon where Vq,n+l,t+l n+l ~ r (V x u)9,n+5,t+l r (3.42) Ar r +i + r n q,n,t+l n v ri Note that the subscript q is being used in the place of H in these diffusion terms. n The subscript q represents the discretized locations within the z-domain (i.e. z\ — 0 cm, Z2 = 0 cm + Az, etc.). H marks a special value of q at each radial location r . H n n n holds the bottom element in the ^-coordinate that is still within the fluid. Using these abbreviations allows the discretized form of (3.37) to be written as (V x t + 1 - (V x u ) 9 i n + i t w [f + (V xu) ± } qn+ H At 1 — exp sin 2 7 n + i + T•n+i „ x l S i +' cos w i + t 2 7 n + +k n 1 — exp + < v exp —- , —TT— 'Ar(V X + \vexp t i ] ^ A (V x « ) ZZ - 2 sin , 7_ r r i i 0 2 ,l + c o s + A ( V x «)(,,„+!,*) ( g > n + 7,1 (3.43) The velocity can then be back-calculated out of (V x u) t+\niq> ( v x u) , +u+i q n ( n+i r + r ) Ar n + ^n+l,g,t+l?"n+l n,q,t+1 V (3.44) As the boundary condition states that VN, ,t = 0, v is calculated by means of (3.44) q starting at n = N — 1 and traveling radially inward. 3.4.2 F l o w C o m p a r i s o n F o r the Stratified S p i n - u p M o d e l Using the same methods for determining initial conditions as described in §3.3.6, the model for stratified spin-up can be tested as well. The following figures (Figure 3.13 through Figure 3.17) show laboratory and theoretical results at 5 second intervals. The Chapter 3. Laboratory Spin-up Without a Canyon 60 data represented in these plots are from an angular section of the tank (the same region illustrated in Figure 3.10). As there is a wave traveling through the tank, this can sometimes cause local variations in the velocity field. A good example of this can be seen in the 10th panel (45 second-panel) of Figure 3.14. Another effect of this wave is directly linked to the sampling strategy. As mentioned in §3.3.6, the choice of temporal location for capturing the initial conditions was done so as to be at the trough of each wave. As a result, one can see in almost all comparisons that the prediction leads the laboratory initially, but that after one rotation period, the two match together again relatively well. It is the comparison between theory and lab after one rotation period that is important. Although there is a tendency for the velocity near the edge of the tank to approach zero faster in the theoretical expectations than in the laboratory results, the final implications of this are not severe, as is explained in greater detail in the following chapter. 3.5 Summary This chapter examines spin-up of a both homogeneous and stratified fluids in labo- ratory when a canyon is absent. The theory of how the boundary layer will evolve over a steep slope is examined, and the Ekman suction over the slopes present in laboratory topography is found. Analysis of the solution suggests that, so long as the Rossby number stays below 0.5 and the ratio of Ekman layer depth to the characteristic lengthscale of the interior is small, the flux into the Ekman layer will follow linear theory. As the slope of the boundary increases, the proportion of this flux that results in vortex tube stretching of the interior decreases as the cosine of the slope angle. Introduction of stratification to the fluid in the laboratory results in a transport of dense water upslope, which leads to the development of buoyancy forces in the boundary layer. Eventually the buoyancy effects match rotation and the bottom boundary layer 61 Chapter 3. Laboratory Spin-up Without a Canyon -0.01 * 5 seconds 0 seconds -0.03 I 1 -0.01 —at". • •*> • ' L~*—-~^ ,» « r -0.01 ** ' JV-»——^ 1 -0.01 — . • • _—^~-«—^^T- • " * * * *— . »^ •• • —'T • * 1, " *•" * ' •^^J > 45 seconds 40 seconds 0.3 • 35 seconds .— -0.03 • _y 30 seconds -0.03 • 25 seconds ' -0.01 • _ i S - * - i — * * T . 20 seconds -0.03 — 15 seconds 10 seconds ' I -c- * ii-—•— -0.03 I 0.4 i • 0.3 0.4 Figure 3.13: Spin-up for trial 1 (N = 2.2 s" , / = 0.52 s - l , A / = 0.12 s - l ) without a canyon, comparison of laboratory observations (dots) to theoretical predictions (lines). Error on data and line are only slightly larger than points and line thickness. The location of the shelfbreak is shown as a dotted line at r = 22 cm. Time slices are shown every 5 seconds. One rotation period is 24.2 seconds. 1 62 Chapter 3. Laboratory Spin-up Without a Canyon -0.01 • 5 seconds 0 seconds -0.03 • 1 1 • I 1 -0.01 15 seconds 10 seconds -0.03 i I i -0.01 • 1 -0.01 -0.03 i i t • !_«—tii- —^^^^ 1 35 seconds 30 seconds • • 1 -0.01 -0.03 25 seconds 20 seconds -0.03 4 " • • i 1 | »• • L X - - - - . .-Si*-—-*"*""""""" • " 45 seconds 40 seconds • • 0.3 0.3 0.4 0.4 Figure 3.14: Spin-up for trial 4 (N = 4.4 s" , / = 0.52 s - 1 , A / = 0.12 s-1) without a canyon, Plots are the same as Figure 3.13. One rotation period is 24.2 seconds. 1 shuts down, arresting all Ekman suction and, therefore, vortex tube stretching. The boundary layer acts as an insulator between the interior and the bottom boundary. Spinup continues to take place by diffusion, although at a decreased rate due to the insulated boundary. Laboratory results match well with theoretical expectations of flow evolution. Despite the presence of a wave traveling in the tank with a period of 4 7 r / , theoretical - 1 prediction of the flow still match with laboratory data when the two are compared one wave period after modeling begins. Chapter 3. Laboratory Spin-up Without a Canyon 63 Figure 3.15: Spin-up for trial 6 (TY = 2.2 s" , / = 0.7 s - 1 , A / = 0.12 s-1) without a canyon, Plots are the same as Figure 3.13. One rotation period is 18.0 seconds. 1 Chapter 3. Laboratory Spin-up Without a Canyon 64 Figure 3.16: Spin-up for trial 7 (N = 2.2 s" , / = 0.52 s - l , Af = 0.09 s - l ) without a canyon, Plots are the same as Figure 3.13. One rotation period is 24.2 seconds. 1 65 Chapter 3. Laboratory Spin-up Without a Canyon •••»"••• *«— -0.01 5 seconds 0 seconds -0.03 , -J, | * • » • • 1 * -0.01 -0.03 -0.01 -0.03 • • 25 seconds 20 seconds .• • •-i * 5 35 seconds 30 seconds • 1 •— —i- „ „ A* -0.01 -0.03 • • -0.01 -0.03 15 seconds 10 seconds • • _—* 45 seconds 40 seconds • i i I 1 0.3 0.4 I i i 0.3 0.4 I 1— Figure 3.17: Spin-up for trial 8 (N = 2.2 s" , / = 0.52 s - l , Af = 0.06 s - l ) without a canyon, Plots are the same as Figure 3.13. One rotation period is 24.2 seconds. 1 Chapter 4 Effects o f a S u b m a r i n e C a n y o n o n t h e L a b o r a t o r y F l o w 4.1 Introduction In the previous chapter, the problem of how stratified water in the tank with canyon- free topography would respond to a forcing event was studied. The results of that work will now act as the foundations for examining the effect of a canyon on flow evolution. Intuitive expectations are that the presence of any angular-dependence to the tank topography will decrease the spin-up time of the fluid, although to my knowledge, studies examining spin-up over variable topography in a cylindrical tank do not exist, with the exception of one study that examined spin-up over a flat bottom with a radial channel (Long, unpublished data). Studies of flow over topography on the continental shelf have shown that oscillatory or random forcing will result in a mean current (Holloway, 1992). This mean current results from unequal forcing effects of the topography in a rotating environment and can be quite strong. As is discussed in Chapter 1, such effects are seen in the flow field after it has passed a submarine canyon. Using the model presented in Chapter 3, the evolution of the flow in the tank with a canyon present can be compared to how the flow would evolve in the absence of a canyon. Comparing flow evolution of a laboratory trial with and without a canyon present show that the difference in flow evolution is considerable (Figure 4.1). The rate of flow evolution over the continental shelf with a canyon present is much quicker. 66 Careful Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 68 examination of the two trials shown in Figure 4.1 shows that their initial conditions (i.e. the velocity distribution across the shelf at the end of forcing) are different due to variations in flow evolution during the forcing period. Experiments performed for this thesis examined only the surface currents in the tank. Quantitatively determining the upwelling flux through a submarine canyon is one of the key questions driving research on currents near submarine canyons. The challenge, then, has been to determine how currents traveling near or through the submarine canyon will affect the flow evolution, and how one can examine the flow evolution to determine the amount of water upwelled through the canyon. 4.2 Q u a l i t a t i v e E v i d e n c e of C a n y o n Effects Before beginning quantitative analysis, it is useful to develop an understanding of how a canyon might affect the flow on the continental shelf. The flow behaviour is easily visualized for a homogeneous fluid in the laboratory. Pathlines can be generated by overlaying images in the flow (Figure 4.2), showing that the canyon effect is large. At the end of the forcing period (Figure 4.2a), adjustment to the presence of a canyon has begun; a towards-shore deviation in the flow can be seen downstream of the canyon. Over the next 30 seconds the flow evolves further (Figure 4.2b). The location where the flow begins to deviate in the up-canyon direction moves upstream - likely due to the decrease in flow rate past the canyon. It also seems to be getting stronger compared to the flow on the shelf. This could be because the water travelling over the slope spinsup more slowly than the water over the shelf. As a result, when it is advected over the shelfbreak through the canyon, its velocity is larger than the autochthonous water present over the shelf. Hence, the effect of water being carried onto the continental shelf from the continental slope is becoming larger compared to the flow travelling over Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow the shelf. 70 The apparent absence of particles downstream of the canyon may be due to the advective upwelling causing a surface divergence, although with no visualisation of the currents upstream of the canyon, one cannot be sure. These effects are further emphasized another 30 seconds later (Figure 4.2c). At this stage, the flow on the shelf has halted almost completely, as can be seen downstream of the canyon near the edge of the tank. The trials that receive quantitative analysis in this chapter are all stratified. As such it is emphasized that the effects seen here are qualitative and meant for visual or conceptual aid. The strong canyon effects on surface flow in a homogeneous fluid will not propagate to the surface in the same manner for a stratified fluid. The effects seen in Figures 4.1 and 4.2 are due to the upwelling component of the flow through the canyon, but they are also due to the drag effects imposed by the canyon on the flow as the two are correlated. 4.3 Drag When a fluid travels past any solid body, the body will impose a set of forces on the fluid. There will be a tangential stress, which is due to the inertial viscosity of the fluid, and there will be a normal stress. The normal stress is comprised of • a lift force, which is the component of the total drag force that is normal to the direction of motion; • an induced drag, which is due to the transfer of kinetic energy between the fluid and the body, and manifests itself as a trail of vortices; and • a form drag, which is the resultant pressure force parallel to the velocity (but in the opposite direction) after the induced drag has been subtracted (Batchelor, 1967). Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 71 Upwelling through a canyon is due to inviscid dynamics in the fluid-topography interaction. It follows that the slowing effects seen in Figure 4.1 are almost entirely due non-frictional drag (referred to hereafter simply as drag). It is the stresses normal to the canyon topography but parallel to the flow direction (induced drag and form drag) that will be used to determine upwelling here. The strong, almost swirling effects seen in the homogeneous case, supply evidence that as the currents travel past the canyon they cross bathymetry contours. This behaviour results in the flow direction not being parallel to the topography, creating drag. As is illustrated in §1.1.1, upwelling in submarine canyons occurs because the sudden changes in topography are strong enough that geostrophy breaks down and the water column experiences vortex tube stretching as it is advected across the bathymetry contours of the canyon. A n alternative way to visualize the problem is that as the fluid passes over the canyon, it feels a drag force imposed by the topography which results in upwelling (Figure 4.3). The presence of this drag force, therefore presents a method for considering canyon upwelling in an evolving flow. While this method does not permit one to understand the provinces of flow within a canyon, it does provide a technique for quantitatively determining upwelling through a submarine canyon in the laboratory. 4.3.1 Determining Canyon D r a g Contributions To The E v o l v i n g Flow Chapter 3 examines the important properties governing stratified spin-up in the laboratory. The results demonstrate that the important variables are the shear flow in the fluid, which determines the strength of flux into the bottom Ekman layer, as well as the stratification and bottom slope, which determine the rate at which the bottom boundary layer shuts down, eventually invoking diffusive spin-up over a slippery bottom boundary layer. Using velocity measurements for trials where a canyon is present, it is possible to use Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 73 the spin-up model derived and tested in Chapter 3 to determine the effect of a canyon on flow evolution. Velocity measurements were taken from trials with a canyon present and entered as initial conditions into the theoretical "no canyon" model. After one rotation period (4.7r/), the model velocities are compared to the measured velocities. If the effect of a canyon were non-existent, then one would expect the model output to agree with the laboratory results. It is already known that the canyon does affect flow evolution (Figure 4.1), and a difference between the model predictions and laboratory measurements is anticipated. This difference is used to calculate the difference in along-shore flux over the shelf, A $ , as A$(r) = Av(r)h(r), (4.1) where Av is the difference between the velocity in the laboratory and the model velocity. The difference in flux can be used to calculate the drag imposed by the canyon using Newton's second law of motion, Force = Mass x Acceleration: F = (pV)(Av/At), (4.2) which can be rewritten as (4.3) where FAt is the impulse imposed on the flow by the canyon in the time At, and F is the average force applied by the canyon over that time. Note that the force is calculated for the entire water mass within the radial limits of integration (through the entire angle 2TT), rather than simply over the canyon. This domain of integration is chosen because the measured flow evolution is assumed to be representative of the effect of the canyon on the entire shelf, so the force that is measured is one that is effectively applied to the entire mass of water over the shelf. An example of the results obtained by this method is shown in Figure 4.4, except that in the figure, (4.3) has not yet been integrated. Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 74 Figure 4.4: Calculation of drag by the canyon from model output, a) Initial conditions of flow (solid line) from laboratory measurements (dots); b) canyon-free expectation of flow evolution (solid line), and actual flow evolution (dashed line) fit to laboratory data (dots); c) points shown in a) with uncertainty and the line fit through the data in a). Dashed lines represent the uncertainty of the fit; d) same as c) except for b); e) difference in velocities given in b) (solid line) and calculated drag from the canyon (dashed line) dr is equal to the radial spacing in the theoretical model. Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 75 The uncertainty in the curves shown in Figure 4.4c and d are calculated using the same bootstrap-like method described in §3.3.6. To determine the uncertainty in predicted flow evolution, the three curves in Figure 4.4c (the average curve of best fit and the two "error curves") are then each run through the theoretical model. The outputs of the error bars represent what the horizontal current profile would look like if it started as the error bar. As a result, the uncertainty given by the outputs will not necessarily have a normal distribution and was calculated using the method shown in Figure 4.5. 4.4 A n a l y s i s of R e s u l t s One can make out a measurable effect of the canyon on spin-up, and this can be expressed as a drag force. Traditionally, a drag force in a non-rotating fluid is be described as F = where p is the density, d (4.4) pCdL^U 2 is the drag coefficient, L\ and L are the appropriate length2 scales, and U is the velocity. For the purposes here, p is assumed to have the constant value, 1 g c m L 2 - 3 (all units used in this analysis are c.g.i. units). The lengthscales L i and are not determined in this thesis and are combined into a constant, allowing (4.4) to be rewritten as break. The relationship between the drag force and the square of the velocity suggests that there is some variation in Cd between the various trials (Figure 4.6). (The results presented in this section are presented graphically. Raw data for U,N,f and F are compiled into tables and presented in Appendix A.) The stratification (TV), shelfbreak velocity Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 76 Radius Figure 4.5: Schematic showing the method used for determining errors in flux measurements. The thick solid line represents the initial conditions taken from the laboratorytrial. The thick dashed line represents the model predictions for how the flow field shown by the solid line would have evolved in the absence of a canyon. The thick dotted line represents how the flow in the laboratory with a canyon present actually evolved. Thin lines that flank the thicker lines of the same format represent the uncertainty in those measurements. The value for Av is equal to the distance expressed by line 2. Line 1 and line 3 represent minimum and maximum values, respectively. 77 Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 10' :l: : : : : : : : I: I : : : : : :l: : : : : 0 :::::.^::::::<T|L, oT 10 -1 in 1 E o TO LL * : AT • • A- 10 : v:::::::::;: V 10" j 0.4 0.5 i_ mmmmmmm J 0.7 , JED. o Trial 1 ::: A Trial 4 o Trial 5 0 Trial 6 I •': V 6 \• 0.6 O ; i i i 0.8 0.9 1 Trial 7 Trial 8 • '•; U [cms- ] 1 sb 10 10 tn E o LL 10 10 0.4 0.5 0.6 0.7 0.8 0.9 1 U [cms- ] 1 sb Figure 4.6: Drag forces as a function of the shelfbreak velocity squared are presented for all trials. The data in the top figure are identical to the bottom figure except that error bars are absent in the top figure for clarity. The bottom figure also shows the line of best fit through the data, giving a relationship of F oc U . 2 0 2 78 Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow (U b), and rotation rate (/) appear to be important in determining drag effects by the S canyon and it is expected that Cd = Cd{U b, N, / ) . The relationship between each of S these dimensional parameters and F/(pU ) 2 b is now compared. Although relationships can be found between these variables independently, it will be shown that they must be compared in combination to determine a final relationship. 4.4.1 Shelfbreak V e l o c i t y The shelfbreak velocity is expected to be one of the key parameters governing the drag imposed by the canyon. The faster the fluid travels, the more it will deviate from a geostrophic balance, resulting in a stronger deviation of the current paths from the isobaths. Perhaps it is surprising, then, that the results shown in Figure 4.6 give no real dependence upon velocity past U . 2 b 4.4.2 Stratification Stratification will affect the canyon spin-up because only water that is lifted out of the canyon and onto the shelf will contribute to the drag felt by the water over the continental shelf. If the stratification is increased, then a higher amount potential energy must be gained by a fluid parcel travelling up through the canyon from a given depth. If the stratification continues to increase, a time will arrive when fluid particles originating from a given depth will not reach the top of the canyon and therefore will not affect the measured drag at the surface. Examination of the relationship between stratification (N) and drag suggests that the drag is related to the stratification as N~ 0 83 (Figure 4.7). 79 Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow Figure 4.7: Left: the relationship between N and drag, showing data from all trials being compared. A line of best fit to the data provides the relationship F cx AT-°- . Right: the relationship between / and drag, showing data from all trials being compared. The line of best fit gives the relationship F oc / ° . 83 1 4 80 Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 4.4.3 R a t e of R o t a t i o n Not including its contribution to the Rossby number, / is expected to play a minor role, if any, in the drag force imposed by the canyon. This is because rotational effects should not affect the magnitude of the drag against the fluid once it has left the bathymetry contours and entered the canyon. There apparently is some other dependence on / as its exponent correlates positively with F/(pUg ). (Figure 4.7) b 4.4.4 C o m b i n a t i o n of D i m e n s i o n a l P a r a m e t e r s The result determined from the independent relations of U b, N and / to F/(pU ) 2 s does not produce a dimensionally feasible solution. b Using these results would give F/(pU A) = /°- /V-o-83r/o.o2 which gives C the approximate dimensions of s 2 14 b d - 0 7 . This is clearly incorrect, as Cd should be dimensionless. It is likely that the variations of each of these individual parameters is large enough that they cannot be solved individually and must be solved together as a set, i.e.: J_ P sb = Af^mrjn ( 4 6 ) U or log(F) - 2 \og{U ) = log(4.) + k log / + m log N + n log U 8b where A r = Ap. 0 01 b (4.7) This equation can be solved using a least squares fit which gives F/(pU ) = Af- - N-°- U° . 2 ab 96 36 b Unfortunately, this solution does not match the re- quired dimensions of the problem either. Given the difficulty in scaling F exactly (or nearly exactly) with dimensional coordinates, it makes sense to attempt to force a nondimensional nature upon Cd by scaling in non-dimensional space. 81 Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 4.4.5 N o n - D i m e n s i o n a l P a r a m e t e r i z a t i o n of Cd The dimensionless parameters that can be used to scale the flow are the Rossby number, the Froude number and the Burger number. If it is assumed that Cd = Ro Fr Br k m n then (4.5) can be rewritten as log(F) - 21og(i7 ) = \og(A ) + k\og(Ro) + m l o g ( F r ) s6 r (4.8) +n\og(Br). The solution of this problem by a least squares estimate inverts a near-singular matrix suggesting that the choice of scaling is probably poor. The solution to (4.8) is found by rewriting it in matrix form as F = PQ, where P are the parameters (1, \og(Ro), log(Fr) and log(-Br)) and Q are the exponents. The values for Q are determined as Q = (P'P)- P'F. To test the validity of this solution for (4.8), the columns of the P were rearranged. This rearrangement resulted in a different answer for Q, suggesting l that the choice of scaling was bad, indeed. Examination of the non-dimensional parameters identifies the problem in the choice of scaling to be Br — S (Ro/Fr), c where S is a constant that is dependent on the spatial c scalings in these parameters. Regardless of which of these three parameters is dropped, the final outcome is the same. If the Burger number is dropped then the least-squares problem solves (4.9) log(F) - 21og(C/, ) = log(A.) + k\og{Ro) + m l o g ( F r ) . 6 Error in the approximation is made using the bootstrap method and the predicted values are log A = 28±14, k = - 0 . 4 7 ± 0 . 2 7 , and m = 0.86±0.24. Using these values for A,m and n, (4.9) predicts the behaviour of the drag (Figure 4.8). Recalling that Ro = U t,/(flZ) s and Fr = U b/(NH), Cd can be written as s 7^0.47±0.25 rj0.49±0.35 J 0 . 4 7 ± 0 . 2 5 C = d ^0.86±0.25 _/y0.86±0.25 (4.10) 82 Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow 10' o Trial 1 A Trial 4 O Trial 5 0 Trial 6 V Trial 7 Trial 8 <N W 3 Figure 4.8: Relationship between drag, velocity, Rossby number and Froude number, based on the least squares solution to (4.8). Data are presented on a logarithmic scale. If the data fit the scaling proposed perfectly, they would lie on a straight line with a slope of 1. The line placed present in the figure has a slope of 1. It was not fit to the data but is meant to guide the eye along the correct aspect ratio. The drag induced by the canyon is therefore £0.47±0.25 F = p(28cm ± l W ) £j2.49±0.35 x0.47±0.25 J^ekii 2 f f l 8 M K > C- ) 4 11 or F^(30cm ) 2 r"0.5 0.5 7T2.5 7-/2.5 <f0.5 -c — ( 4 . 1 2 ) The upwelling velocity can be determined from this force, as is expressed schematically in Figure 4.3: - F (4.13) ~Jv~ Vup P where v up is the average cross-shelfbreak velocity imposed by the canyon and V is the volume of water affected by the force. In the laboratory, the upwelling flux can be calculated as $ up = 2nr h v , sb sb up (4.14) 83 Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow where <J> is the cross-shelfbreak flux due to the canyon and 2nr h up sb is the area of the sb cylindrical wall of water directly above the shelfbreak in the tank. Based on laboratory data, the upwelling flux in the tank is equal to $ 4.5 — K up ? 3 F . (4.15) V g / Case S t u d y : P r e d i c t i n g the U p w e l l i n g F l u x T h r o u g h A s t o r i a C a n y o n Taking the findings of this study, it is possible to make predictions for upwelling in coastal submarine canyons. Since most of the relevant information necessary to make this prediction is already in Table 2.1 for Astoria canyon, it will be used as the example. Typical scales for Astoria canyon are U = 20cm s , N = 7.5 x 1 0 s , / = 1 0 s , -1 H = 150m, L = 22km, Z = 60m, W sb _ 3 - 1 - 4 - 1 = 16km and IZ = 4.5km, where Z is the depth below the canyon head from which the deepest upwelled water originates as defined by Allen and Hickey (in prep.). If it is assumed that the scales of A increase by O(10 ) 6 (corresponding to a change in spatial scales from cm to k m - see (4.5)), then the drag 2 2 induced by Astoria canyon will be approximately 3 x 10 N . If the volume of Astoria 8 canyon is approximately LW Z sb and the area through which the flux enters is ZW , sb then the upwelling flux for Astoria canyon is *- = JT ' (4 16) P which is 1.3 x 1 0 m s . 5 4.6 3 _1 Summary The flow evolution of laboratory spin-up experiments with a canyon present were compared to theoretical expectations for the spin-up in the absence of a canyon. The theoretical model used velocity measurements from the laboratory as initial conditions 84 Chapter 4. Effects of a Submarine Canyon on the Laboratory Flow and laboratory data were compared to model output one tank rotation period later. The differences in velocity were used to determine the drag force imposed on the fluid by the canyon. Examining the relationships held between the form drag and the appropriate non-dimensional parameters suggests that the drag felt by the canyon will adhere to the relationship C/5/2 f1/2 The upwelling flux onto the shelfbreak due to the canyon is determined as <J> oc F/f, u p or TJ5/2 Analysis of these results findings for Astoria canyon suggests that the upwelling flux through Astoria upwelling favourable conditions is O(10 m s ). 5 3 _1 Chapter 5 Discussion a n d C o n c l u s i o n s The goal of this thesis was to determine how water and current properties upstream of a submarine canyon will affect the amount of upwelling through it. The problem was solved by studying how the presence of a submarine canyon affects the spin-up rate of a fluid in the laboratory, and using the flow evolution to determine drag induced by the canyon. This study is the first where the drag is measured through analysis of the spin-up problem. The measured drag can be used to estimate upwelling through the submarine canyon. 5.1 Laboratory Methods Originally, eight trials were intended to be run and compared. In Trial 2, the fluid was homogeneous and the results were necessary to verift the Ekman suction model presented in Chapter 3. While results for Trial 2 with a canyon present were used for qualitative visualisation of canyon effects on flow over the shelf, predictions of the ensuing form drag were not made. Results for Trial 3 with a canyon present were unfortunately lost and could not be analyzed. Laboratory measurements made in this thesis were contaminated by an oscillation in the velocity due to the fact that the axis of rotation was not perfectly normal to the vertical. Although attempts were made to level the tank, it seems that the accuracy required was beyond the methods used and made it difficult to examine the flow late in time. Nevertheless, quantitative predictions of how the canyon affects the flow were still 85 86 Chapter 5. Discussion and Conclusions found and match well with intuitive expectations of the flow behaviour. The use of P I V methods worked as an excellent tool for flow measurements. This project reintroduced the methods to the laboratory, and although the error in measurements still allowed meaningful results, the uncertainties could likely be reduced further. 5.2 Laboratory Spin-up W i t h and W i t h o u t a Canyon Analysis of the momentum equations for the boundary layer of a homogeneous fluid on a slope provides three important observations: • the thickness of the boundary layer will increase as (cos7) / ; -1 2 • the strength of vortex stretching in the interior due to Ekman suction will remain a constant of the flow vorticity; • Ekman suction will have a horizontal component proportional sin 7. These expectations are consistent with laboratory observations. The behaviour of the bottom boundary layer in a stratified, rotating fluid over a shallow slope is described by MacCready and Rhines(1991). They predict the boundary layer to shutdown on a timescale T = (l/cr + S)/[(l + S)S f cos7], after which diffusion 2 over a slippery bottom boundary layer will govern the spin-up process. In order to consider the effects of horizontal diffusion in regions of high horizontal shear, the diffusion solution proposed by MacCready and Rhines (1991) is modified slightly. This theory is used to successfully replicate the spin-up results for a stratified fluid in the laboratory when a canyon is not present. The introduction of a submarine canyon to the laboratory topography changes the spin-up regime in the laboratory. Using the theory mentioned above, comparisons of the flow evolution with and without a canyon provided a measurable effect of its presence, 87 Chapter 5. Discussion and Conclusions from which the force (or drag) imposed on the fluid by the canyon was approximated. The experiment was repeated for various values of / , U and ./V and an empirical relationship between the drag force and these parameters was found. Analyzing each of these parameters individually would give a possible scaling of F oc U N~ , 2 1 but this scaling is dimensionally incorrect. Scaling by U and the non-dimensional numbers Fr and Ro 2 gives a dimensionally correct scaling of F oc U ^ / / ./V . 5 2 _1 2 _1 The relationship between the drag force and the upwelling velocity is quite simple, which makes it possible to predict the relationship of the upwelling flux through the canyon to the above parameters as p 5.3 N f'' 1 2 Scatter i n D a t a M e a s u r e m e n t s The data presented in Chapter 4 displays a large amount of scatter, which is most likely because the velocity used for drag calculations is derived from the shelfbreak region, which falsely assumes that the flow over the length of the canyon is constant. While this assumption seems correct in the real ocean, it does not hold in the laboratory. By the end of the forcing period, the velocity at the head of the canyon is about 30% of what it is at the shelfbreak (Figure 4.1b). While a large portion of the drag experienced by the flow will be from water whose source is near the shelfbreak, some drag is likely from water closer to the head of the canyon. Also, high frequency changes in the velocity (seen in Figure 4.1) could cause values of the predicted drag to be higher or lower, depending on if the data were sampled in a high or low point in these bumps (either in the initial conditions or in the final measurements). This occurrence would alter the expected drag value considerably. These high-frequency oscillations are radially and azimuthally continuous, suggesting that it is not an artifact 88 Chapter 5. Discussion and Conclusions of the sampling method, such as windowing or biasing in the PIV process. Instead, it is likely that they are due to inconsistencies in table rotation, as they appear to be of the same magnitude as the precision of the instruments that control the rotation rate of the table. Despite these uncertainties, the signal of the upwelled water is strong and consistent enough that a believable relationship is still found. A scaling analysis of the flux through a laboratory canyon being developed by Allen and Hickey (in prep.) quantify the amount of water travelling through a canyon, </>, as <f> = W ZU FH, sb where Z = (y/fU L?)/N, ab H = ^Jl/(ARl) (5.1) sb T = 2(Ro + G- RoG)/3, G = -(l./Ro + 1/Ru, - l/(2R ) w and R w = U T / {fW ). sb sb W sb - ^4 + l/Ro )/2, 2 is the distance across the canyon at the shelfbreak, Z is the maximum depth from which water can enter the canyon and exit onto the shelf and U J-H sb is the amount of water entering the canyon at its head. The findings presented in Chapter 4 support this theory (Figure 5.1) and a linear relationship between them appears to be present. 5.4 I m p l i c a t i o n s of F i n d i n g s In §4.5, the findings of this study were applied to Astoria canyon to determine its upwelling flux as 1.3 x 10 m s . 5 3 _1 To put this value in context, it is compared to a prediction of wind-forced upwelling at the shelfbreak. Assuming an Ekman balance between the wind stress and the ocean currents, the wind stress, r, will create a surface flux U = T, y Pf v = - - (5.2) f PJ (5-3) 89 Chapter 5. Discussion and Conclusions "TJT O A m c • 0 v 3 tn CC CD Trial 1 Trial 4 Trial 5 Trial 6 Trial 7 Trial 8 4 O o 0 Air o o O. A JD CO V ^7 J4_ <y 3 4 Theory Figure 5.1: Comparing upwelling measurements from the laboratory with those predicted by the theory described by Allen and Hickey(in prep.) (<j>). Both are in units of c m s . A line of best fit through the origin gives a slope of 0.3 ± 0.02. 3 where the axes define a right-handed system relative to the coast with y increasing in the offshore direction, / is the Coriolis frequency, p is the density, u and v are the velocities in the x and y direction, respectively, and r and r are the wind stress in the x and y x y directions, respectively. Lentz (1992) found that the mean wind stress over the Oregon shelf from July 5 to August 28, 1973 was parallel to the coast with r = T = -0.45 dyn c m , where 1 dyn - 2 X = 1 0 N . For a fluid density of 1000 kg m _5 the wind will be 0.45 m s 3 _ 1 - 3 the offshore flux of surface water due to for each meter along the shelfbreak. By continuity, there must be a towards-shore return flow of equal magnitude. Astoria canyon has width of 1.6 x 10 m. A region of shelfbreak equal to the shelfbreak width of Astoria canyon will 4 therefore have a towards shore flux of 7.2 x 1 0 m s , suggesting that the flux through 3 3 -1 Astoria canyon is almost 20 times larger than that over a section of shelfbreak of equal width. - 1 Chapter 5. Discussion and Conclusions 5.5 90 Future W o r k The spatial and temporal features of a given canyon are unique. The magnitude of difference in upwelling determined between Astoria canyon and an appropriate length of shelfbreak is not identical for other canyons. A project gathering the dimensional parameters of all the canyons off the west coast of North America has already begun. These data will be used with the necessary flow and water property information to determine the upwelling flux of water through each of these submarine canyons. The laboratory work is to be followed by a similar set of experiments as those performed in this thesis except that their goal will be to determine the role of a canyon in down-welling favourable flow, which is characteristic of the winter season off of the coast of British Columbia. This may prove to be quite challenging due to the behaviour of boundary layers on a slope in a stratified fluid. In order to force downwelling currents in the laboratory, one must decrease the rotation rate. This will add cyclonic vorticity to the water in the tank and the water in the boundary layer will travel radially inward. As it travels downslope, however, the water in the boundary layer will be advected a region of higher density, causing it to become unstable. Instabilities in boundary layers travelling downslope have will have mixing efficiencies an order of magnitude higher than those travelling upslope (Condie, 1999). A second proposed laboratory study, which may help to explain the scatter in the drag-velocity relationship, would develop a three-dimensional map of submarine canyon flow. Bibliography Allen, S. (1996). Topographically generated, subinertial flows within a finite length canyon. J. Phys. Oceanogr., 26(8): 1608-1632. Allen, S. (2000). On subinertial flow in submarine canyons: Effect of geometry. J. Geophys. Res., 105:1285-1297. Allen, S., Dinniman, M . , Gorby, D., Klinck, J., Hewett, T., and She, J. (submitted). Numerical and physical modelling: A comparison over steep topography. J. Geophys. Res. Allen, S. and Hickey, B . (in prep.). Dynamics of advection-driven upwelling over a submarine canyon. Allen, S. E., Vindeirinho, C., Thomson, R. E., Foreman, M . G . , and Mackas, D. L . (2001). Physical and biological processes over a submarine canyon during an upwelling event. Can. J. Fish. Sci, 58(4)-.671-684. Batchelor, G . (1967). An Introduction to Fluid Dynamics. Cambridge University Press, Cambridge, 615 pp. Benton, E. (1973). Nonlinear hydrodynamic and hydromagnetic spin-up driven by ekmanhartmann boundary layers. J. Fluid Mech., 57:337-360. Benton, E. and Clark, Jr., A . (1974). Spin-up. Ann. Rev. Fluid Mech., 6:257-280. Boyer, D. L . and Davies, P. A . (2000). Laboratory studies of orographic effects in rotating and stratified flows. Ann. Rev. Fluid Mech., 32:165-202. 91 92 BIBLIOGRAPHY Charney, J. and Eliassen, A . (1949). A numerical method for predicting the perturbations of the middle latitude westerlies. Tellus, 1:38-54. Chen, X . and Allen, S. E . (1996). The influence of canyons on shelf currents: A theoretical study. J. Geophys. Res., 101:18043-18059. Condie, S. A . (1999). Ocean boundary mixing during Ekman layer arrest. J. Phys. Oceanogr., 29:2993-3001. Durrieau de Madron, X . (1994). Hydrography and nephaloid structure in the GrandRhone canyon. Contin. Shelf Res., 14:457-477. Ekman, V . W . (1905). On the influence of the earth's rotation on ocean currents. Arkiv. Matem., Ast. Fysik, 2:1-52. Foreman, M . G . G . and Thomson, R. E. (1997). Three-dimensional model simulations of tides and buoyancy currents along the West Coast of Vancouver Island. J. Phys. Oceanogr., 27:1300-1325. Freeland, H . and Denman, K . (1982). A topographically controlled upwelling center off Vancouver Island. J. Marine Res., 40:1069-1093. Greenspan, H . (1964). On the transient motion of a contained rotating fluid. J. Fluid Mech., 21:673-696. Greenspan, H . (1965). On the general theory of contained rotating fluid motions. J. Fluid Mech., 22:449-462. Greenspan, H . (1968). 327pp., Cambridge. The Theory of Rotating Flows. Cambridge University Press, 93 BIBLIOGRAPHY Greenspan, H . and Howard, L . (1963). On a time-dependent motion of a rotating fluid. J. Fluid Mech., 17:385-404. Han, G . , Hansan, D., and Gait, J. (1980). Steady-state diagnostic model of the New York Bight. J. Phys. Oceanogr., 10:1998-2020. Hart, J . (2000). A note on nonlinear corrections to the Ekman layer pumping velocity. Physics of Fluids, 12(1):131-135. Hewett, A . (1998). Laboratory simulation of a geophysical flow past a submarine canyon. Honour's thesis at the University of British Columbia. Hickey, B . (1997). The response of a narrow canyon to strong wind forcing. J. Phys. Oceanogr., 27(5):697-726. Hickey, B . M . (1995). Coastal submarine canyons. In Muller, P. and Henderson, D., editors, 'Aha Huliko'a Hawaiian Winter Workshop on topographic effec ts in the ocean, pages 95-110, University of Hawaii at Manoa. S O E S T Special Publication. Hickey, B . M . , Baker, E . , and Kachel, N . (1986). Suspended particle movement in and around Quinault submarine canyon. Marine Geol, 71:35-83. Holloway, G . (1992). Representing topographic stress for large-scale ocean models. J. Phys. Oceanogr., 22:1033-1046. Holton, J . R. (1965). The influence of viscous boundary layers on transient motions in a stratified fluid: part I. J. Atm. Sci., 22:402-411. Holton, J. R. (1967). The diurnal boundary layer wind oscillation above sloping terrain. Tellus, 19:199-205. Hsueh, Y . (1969). Buoyant Ekman layer. Phys. Fluids, 12:1757-1762. 94 BIBLIOGRAPHY Hutchins, D. and Bruland, K . (1998). Iron-limited diatom growth and Si:N uptake ratios in a coastal region. Nature, 393:561-564. Klinck, J. (1988). The influence of a narrow transverse canyon on an initially geostrophic flow. J. Geophys. Res., 93:509-515. Klinck, J. (1989). Geostrophic adjustment over submarine canyons. J. Geophys. Res., 94:6133-6144. Klinck, J. (1995). Circulation near submarine canyons: a modelling study. J. Geophys. Res., 101:1211-1223. Lentz, S. J. (1992). The surface boundary layer in coastal upwelling regions. J. Phys. Oceanogr., 122:1517-1539. Linden, P. and van Heifst, G . (1984). Two-layer spin-up and frontogenesis. J. Fluid Mech., 143:69-94. MacCready, P. and Rhines, P. B . (1991). Buoyant inhibition of Ekman transport on a slope and its effect on stratified spin-up. J. Fluid Mech., 223:631-661. MacCready, P. and Rhines, P. B . (1993). Slippery bottom boundary layers on a slope. J. Phys. Oceanogr., 23:5-22. Mackas, D. L . , Kieser, R., Saunders, M . , Yelland, D. R., Brown, R. M . , and Moore, D. F. (1997). Aggregation of euphasics and Pacific Hake (merluccius productus) along the outer continental shelf off Vancouver Island. Can. J. Aquat. Sci., 57:2080-2096. O'Donnell, J. and Linden, P. (1991). Free-surface effects on the spin-up of fluid in a rotating cylinder. J. Fluid Mech., 232:439-453. Osier, G . (1965). Density gradients. Scientific American, 213:70-76. 95 BIBLIOGRAPHY Pedlosky, J. (1987). Geophysical Fluid Dynamics, 2 na Ed. Springer-Verlag, Ann Arbor, USA, 710pp. Perrene, N . , Verron, J., Renouard, D., Boyer, D., and Zhang, X . (1997). Rectified barotropic flow over a submarine canyon. J. Phys. Oceanogr., 27(9):1868-T893. Press, W., Flannery, B . , Teukolsky, S., and Vetterling (1986). Numerical Receipes: the Art of Scientific Computing. Cambridge University Press, Cambridge, 818pp. Raffel, M . , Willert, C , and Kompenhans, J. (1998). Particle Image Velocimetry: A Practical Guide. Springer-Verlag, Berlin, Germany. Sveen, J. (2000). An introduction to MatPIV v.1.4- University of Oslo, Norway, pp. 18. Thomson, R. (1981). Oceanography of the British Columbia Coast. Department of Fisheries and Oceans. Thorpe, S. (1987). Current and temperature variability on the continental slope. Phil. Trans. Roy. Soc. Lon., A323:471-517. Thurman, H . (1994). Introductory Oceanography. MacMillan Publishing Company, New York, 550pp. van de Konijnenberg, J. and van Heijst, G . (1995a). Nonlinear spin-up of a circular cylinder. Phys. Fluids, 7:2989-2999. van de Konijnenberg, J., Wessels, T., and ven Heijst, G . (1996). Spin-up in a circular tank with a radial barrier. Phys. Fluids, 8(8):2048-2059. van de Konijnenberg, J. A . and van Heijst, G . (1995b). Nonlinear spin-up in a circular cylinder. Phys. Fluids, 7(12):2989-2999. BIBLIOGRAPHY 96 Vindeirinho, C. (1998). Water properties, currents and zooplankton distribution over a submarine canyon under upwelling favourable conditions. Master's thesis, University of British Columbia, Vancouver, B . C , Canada. Weatherly, G . and Marin, P. (1978). On the structure and dynamics of the ocean bottom boundary layer. J. Phys. Oceanogr., 8:557-570. Wedemeyer, E. (1964). The unsteady flow within a cylinder. J. Fluid Mech., 20:383-99. Willert, C. (1996). The fully digital evaluation of photographic P I V recordings. Appl. Sci. Res., 56:79-102. Appendix A Canyon F l u x Measurements The tables presented in this appendix give the raw data that are used in for the production of figures in Chapter 4. The value for the shelfbreak velocity, U b is calculated as the s mean surface velocity between 27 and 30 cm. 97 Appendix A Time s 15.8 25.6 26.4 27.0 30.6 31.2 31.8 35.4 41.6 Canyon Flux Measurements U. 'sb cm s" -1.20 +0.044/-0.044 -1.10 +0.031/-0.031 -1.50 +0.052/-0.052 -1.09 +0.028/-0.028 -1.09 +0.035/-0.035 -1.50 +0.053/-0.053 -1.08 +0.036/-0.036 -1.50 +0.049/-0.049 -1.50 +0.052/-0.052 Drag [g cm s 3.81 11.04 14.73 10.93 6.44 15.61 5.53 15.64 15.34 +4.032/-3.893 +2.515/-2.469 +5.075/-4.829 +2.463/-2.416 +2.669/-2.616 +5.162/-4.919 +2.814/-2.761 +5.074/-4.822 +5.487/-5.216 cm s" 3 2 up 0.52 +0.549/-0.530 1.50 +0.342/-0.336 2.00 +0.691/-0.657 1.49 +0.335/-0.329 0.88 +0.363/-0.356 2.12 +0.702/-0.669 0.75 +0.383/-0.376 2.13 +0.690/-0.656 2.09 +0.747/-0.710 Table A . l : Drag and upwelling flux predictions for trial 1 (/ = 0.52 s \ N = 2.2 s e [s U 10.8 11.6 13.2 13.4 13.6 14.4 14.6 16.2 16.4 17.0 34.8 35.0 35.2 37.0 38.0 59.4 -1 08 +0 -1 09 +0 -1 08 +0 -1 50 +0 -1 10 +0 -1 58 +0 -1 50 +0 -1 10 +0 -1 52 +0 -1 50 +0 -1 50 +0 -1 20 +0 -1 55 +0 -0 91 +0 -1 09 +0 -1 22 +0 sb [cm s- ] 1 Drag 5 73 036/- 0.036 028/- 0.028 10 93 5 53 036/- 0.036 053/- 0.053 15 61 031/- 0.031 11 04 020/- 0.020 17 95 052/- 0.052 14 73 033/- 0.033 11 00 058/- 0.058 10 11 049/- 0.049 15 64 052/- 0.052 15 34 044/- 0.044 3 81 044/- 0.044 17 82 022/- 0.022 12 08 6 44 035/- 0.035 044/--0.044 6 82 [g cm s ] +2 737/- 2 685 +2 463/- 2 416 +2 814/- 2 761 +5 162/- 4 919 +2 515/-•2 469 +3 302/- 3 141 +5 075/- 4 829 +2 615/- 2 570 +5 545/- 5 281 +5 074/--4 822 +5 487/--5 216 +4 032/- •3 893 +5 955/--5 762 +1 934/-•1 867 +2 669/--2 616 +3 668/--3 503 2 $ up cm s ] 3 0 78 +0.372/--0 365 1 49 +0.335/--0 329 0 75 +0.383/- 0 376 2 12 +0.702/--0 669 1 50 +0.342/- 0 336 2 44 +0.449/- 0 427 2 00 +0.691/- 0 657 1 50 +0.356/- 0 350 1 38 +0.754/--0 719 2 13 +0.690/--0 656 2 09 +0.747/-•0 710 0 52 +0.549/--0 530 1 80 +0.602/--0 582 2 14 +0.342/-•0 330 0 88 +0.363/--0 356 0 93 +0.499/--0 477 Table A.2: Drag and upwelling flux predictions for trial 4 (/ = 0.52 s \ N = 4.4 s Appendix A Time [s 4.8 7.4 15.8 16.8 17.6 23.4 24.0 25.6 34.6 35.8 39.2 46.8 49.2 Canyon Flux Measurements [cm s ] l U sb -1.10 -1.20 -1.09 -1.50 -1.10 -1.50 -1.50 -1.08 -1.09 -1.08 -1.22 -1.52 -1.50 +0.031/- 0.031 +0.044/- 0.044 +0.028/- 0.028 +0.049/- 0.049 +0.033/- 0.033 +0.053/--0.053 +0.052/--0.052 +0.036/--0.036 +0.035/--0.035 +0.036/--0.036 +0.044/--0.044 +0.058/--0.058 +0.052/--0.052 Drag g cm s ] $ up c m s ] 11.04 3.81 10.93 15.64 11.00 15.61 14.73 5.53 6.44 5.73 6.82 10.11 15.34 +2.515/-2.469 +4.032/-3.893 +2.463/-2.416 +5.074/-4.822 +2.615/-2.570 +5.162/-4.919 +5.075/-4.829 +2.814/-2.761 +2.669/-2.616 +2.737/-2.685 +3.668/-3.503 +5.545/-5.281 +5.487/-5.216 1 50 0 52 1 49 2 13 1 50 2 12 2 00 0 75 0 88 0 78 0 93 1 38 2 09 2 3 l +0.342/-0.336 +0.549/-0.530 +0.335/-0.329 +0.690/-0.656 +0.356/-0.350 +0.702/-0.669 +0.691/-0.657 +0.383/-0.376 +0.363/-0.356 +0.372/-0.365 +0.499/-0.477 +0.754/-0.719 +0.747/-0.710 Table A.3: Drag and upwelling flux predictions for trial 5 ( / = 0.40s , TV = 2.2 s x Time 10.4 13.8 14.4 15.6 21.8 34.4 35.4 50.6 51.8 Drag U.sb cm s" -1.50 +0.052/-0.052 14.73 -1.10 +0.031/-0.031 11.04 -1.50 +0.052/-0.052 15.34 -1.09 +0.035/-0.035 6.44 -1.50 +0.053/-0.053 15.61 -1.09 +0.028/-0.028 10.93 -1.20 +0.044/-0.044 3.81 -1.08 +0.036/-0.036 5.53 -1.50 +0.049/-0.049 JL5.64 [g cm s - 2 1 $ +5.075/-4.829 +2.515/-2.469 +5.487/-5.216 +2.669/-2.616 +5.162/-4.919 +2.463/-2.416 +4.032/-3.893 +2.814/-2.761 +5.074/-4.822 2.00 1.50 2.09 0.88 2.12 1.49 0.52 0.75 2.13 [cm s 3 M p 1 +0.691/-0.657 +0.342/-0.336 +0.747/-0.710 +0.363/-0.356 +0.702/-0.669 +0.335/-0.329 +0.549/-0.530 +0.383/-0.376 +0.690/-0.656 Table A.4: Drag and upwelling flux predictions for trial 6 (/ = 0.70 s , N = 2.2 s 1 Appendix A Time -12.6 -7.2 -1.0 0.8 6.2 21.2 21.8 35.2 65.6 67.4 Canyon Flux Measurements U,sb cm s -r -1 50 +0.049/- 0.049 -1 52 +0.058/- 0.058 -1 10 +0.031/- 0.031 -1 50 +0.053/- 0.053 -1 20 +0.044/- 0.044 -1 09 +0.028/- 0.028 -1 50 +0.052/- 0.052 -1 08 +0.036/- 0.036 -1 50 +0.052/--0.052 -1 09 +0.035/--0.035 Drag [g cm s - 2 1 15.64 10.11 11.04 15.61 3.81 10.93 14.73 5.53 15.34 6.44 +5.074/--4.822 +5.545/--5.281 +2.515/--2.469 +5.162/--4.919 +4.032/--3.893 +2.463/--2.416 +5.075/--4.829 +2.814/--2.761 +5.487/--5.216 +2.669/--2.616 cm s 3 up 2.13 1.38 1.50 2.12 0.52 1.49 2.00 0.75 2.09 0.88 +0.690/-0.656 +0.754/-0.719 +0.342/-0.336 +0.702/-0.669 +0.549/-0.530 +0.335/-0.329 +0.691/-0.657 +0.383/-0.376 +0.747/-0.710 +0.363/-0.356 Table A.5: Drag and upwelling flux predictions for trial 7 (/ = 0.52 s , N = 2.2 s 1 Time 1.2 3.6 12.4 19.2 22.2 cm s -1.09 +0.028/-0.028 -1.10 +0.031/-0.031 -1.09 +0.035/-0.035 -1.50 +0.053/-0.053 -1.08 +0.036/-0.036 sb Drag [g cm s $ 10.93 11.04 6.44 15.61 5.53 +2.463/-2.416 +2.515/-2.469 +2.669/-2.616 +5.162/-4.919 +2.814/-2.761 1.49 1.50 0.88 2.12 0.75 [cm s 3 u p r +0.335/-0.329 +0.342/-0.336 +0.363/-0.356 +0.702/-0.669 +0.383/-0.376 Table A.6: Drag and upwelling flux predictions for trial 8 (/ = 0.52 s , N = 2.2 s 1
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Spin-up over steep topography and the effects of a...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Spin-up over steep topography and the effects of a submarine canyon Mirshak, Ramzi 2001
pdf
Page Metadata
Item Metadata
Title | Spin-up over steep topography and the effects of a submarine canyon |
Creator |
Mirshak, Ramzi |
Date Issued | 2001 |
Description | Submarine canyons are common bathymetric features that cut into the continental shelf from the continental slope. During upwelling favourable conditions over the continental shelf, submarine canyons display enhanced upwelling, heavily impacting shelf-slope mass exchange. In this thesis, laboratory experiments were conducted to quantify how velocity, stratification and rotation affect canyon upwelling. Currents were forced by changing the rotation rate of an already rotating tank. Over time, the forced currents evolve, or spin-up, until the fluid within the tank is once again rotating at the same rate as the tank itself. A difference in the spin-up behaviour is observed when a canyon is not present in the laboratory topography. The dynamics that govern flow evolution when a canyon is not present are Ekman suction and diffusion, both of which can be affected by a sloped bottom boundary layer. A numerical model is developed that replicates the spin-up of a stratified fluid over the changing slopes of the laboratory topography in the absence of a canyon, incorporating both the decay of Ekman suction and the change in the diffusion of momentum which occur as the boundary layer flow is arrested by buoyancy effects. Observed spin-up with a canyon present is compared to predicted spin-up without a canyon. The difference measures a force imposed by the canyon, which is found to be proportional to U[sup 2.5]f[sup 0.5]/N, where U is the shelfbreak velocity, f is the Coriolis frequency and N is the buoyancy frequency. This force relates to the flux of water crossing the shelfbreak through the canyon. The results are applied to Astoria canyon, suggesting that in a strong upwelling event, the flux of water through the canyon is 1.3 x 10 ⁵m³ s⁻¹ , nearly 20 times larger than wind-induced upwelling over a similar length of the shelfbreak. |
Extent | 4197554 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-08-06 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0052635 |
URI | http://hdl.handle.net/2429/11840 |
Degree |
Master of Science - MSc |
Program |
Earth and Ocean Sciences |
Affiliation |
Science, Faculty of Earth, Ocean and Atmospheric Sciences, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 2001-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_2001-0472.pdf [ 4MB ]
- Metadata
- JSON: 831-1.0052635.json
- JSON-LD: 831-1.0052635-ld.json
- RDF/XML (Pretty): 831-1.0052635-rdf.xml
- RDF/JSON: 831-1.0052635-rdf.json
- Turtle: 831-1.0052635-turtle.txt
- N-Triples: 831-1.0052635-rdf-ntriples.txt
- Original Record: 831-1.0052635-source.json
- Full Text
- 831-1.0052635-fulltext.txt
- Citation
- 831-1.0052635.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0052635/manifest