A NUMERICAL INVESTIGATION OF THE SHELF SLOPE IN THE SCALING FOR UPWELLING OVER SUBMARINE CANYONS by Tara Howatt A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF BACHELOR OF SCIENCE (HONOURS) in THE FACULTY OF SCIENCE (Oceanography and Geophysics) This thesis conforms to the required standard Supervisor UNIVERSITY OF BRITISH COLUMBIA (Vancouver) MARCH 2012 c Tara M. Howatt, 2012Abstract Submarine canyons are known regions of enhanced upwelling. Previous studies have looked into the dynamics of the upwelling and recently a scaling system by Allen and Hickey [2010] has been developed. A closer look into this scaling has shown that the continental shelf slope surrounding the canyon plays a role and needs to be included in the scaling. The slope plays a role in the depth of upwelling and the upwelling ux because the slope will create a back pressure restricting the ow from moving further up the shelf. This thesis will discuss the determination of a non-dimensional number to characterize the shelf slope in these upwelling quantities. A numerical model will be used to test the new scaling and to assist in nding a relation between the shelf slope and the depth of upwelling and the upwelling ux. A non-dimensional number to take into account the e ect of the shelf slope has been found; however, further research is required for a physical explanation of the form of this number. Looking at the results from the model, relationships between the slope, strati cation, Coriolis parameter and upwelling depth and ux can also be made, where a atter slope, low strati cation, and a low Coriolis parameter will have the greatest upwelling response. iiContents List of Figures vii List of Tables viii List of Appendices ix Acknowledgments x 1 Introduction 1 2 Background 3 2.1 Upwelling over Submarine Canyons . . . . . . . . . . . . . . . . . . . . . 3 2.1.1 Phases of Upwelling . . . . . . . . . . . . . . . . . . . . . . . . . 3 2.1.2 Components of Flow Over Canyons . . . . . . . . . . . . . . . . . 3 2.2 Initial Scaling of Upwelling Over Submarine Canyons . . . . . . . . . . . 4 2.2.1 Assumptions of the Scaling . . . . . . . . . . . . . . . . . . . . . . 4 2.2.2 Parameters in the Scaling . . . . . . . . . . . . . . . . . . . . . . 4 2.2.3 Scaling of the Upwelling Depth . . . . . . . . . . . . . . . . . . . 5 2.2.4 Scaling of the Upwelling Flux . . . . . . . . . . . . . . . . . . . . 6 3 Methodology 6 3.1 The Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 3.2 The Experiment Setup . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 3.3 Incorporating Shelf Slope in the Scaling . . . . . . . . . . . . . . . . . . . 10 4 Initial Results 10 4.1 Scaling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 4.2 Time Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 4.2.1 Depth of Upwelling . . . . . . . . . . . . . . . . . . . . . . . . . . 12 4.2.2 Volume of Water Upwelled . . . . . . . . . . . . . . . . . . . . . . 13 5 Inclusion of Shelf Slope in the Scaling 14 5.1 Slope E ect on the Pressure Field . . . . . . . . . . . . . . . . . . . . . . 14 5.2 Scaling for the Depth of Upwelling . . . . . . . . . . . . . . . . . . . . . 18 5.3 Scaling for the Upwelling Flux . . . . . . . . . . . . . . . . . . . . . . . . 19 iii6 Determination and Veri cation of 19 6.1 Early Attempts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 6.2 Later Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 6.3 Error Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 6.3.1 Error in the Depth of Upwelling . . . . . . . . . . . . . . . . . . . 20 6.3.2 Error in the Upwelling Flux . . . . . . . . . . . . . . . . . . . . . 29 6.3.3 Error in Topographic Burger Number and 1 . . . . . . . . . . . 29 7 Discussion 30 7.1 Assumptions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 7.2 Slope E ect on Upwelling Characteristics . . . . . . . . . . . . . . . . . . 30 7.2.1 Depth of Upwelling . . . . . . . . . . . . . . . . . . . . . . . . . . 30 7.2.2 Upwelling Flux . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31 7.3 Scales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32 7.3.1 Depth of Upwelling Scales . . . . . . . . . . . . . . . . . . . . . . 32 7.3.2 Upwelling Flux Scales . . . . . . . . . . . . . . . . . . . . . . . . 32 8 Conclusions 32 References 34 Appendices 36 A Initial Scaling 36 A.1 Depth of Upwelling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 A.2 Upwelling Flux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 ivList of Figures 1 The original scaling system of Allen and Hickey [2010] for a) the depth of upwelling and b) the upwelling ux. The di erent dots represent the di er- ent canyon cases which varied the geometry of the canyon, the strati cation,N ; the Coriolis parameter, f ; and the U velocity upstream of the canyon. The variations in slope are shown by the red and green stars and are o the line of best t. The x in the best t refers to the x-axis. The norm of residuals is normalized to the x-axis. . . . . . . . . . . . . . . . . . . . . 2 2 Schematic of ow over a submarine canyon (from Allen and Hickey [2010]). 4 3 The geometry of the canyon is shown. The white contours represent the 0.02 m shelf-break depth. Wsb is the shelf break width, W is the width of the canyon, L is the length of the canyon from the shelf break width to the head of the canyon, and R is the radius of curvature of the canyon. . 5 4 A contour plot of the Basic slope topography in a tank. The topography contours are in centimeters. The white contours show the 2 cm shelf break depth and in the center is the cylinder. . . . . . . . . . . . . . . . . . . . 7 5 Cross sections of the topography from the center to the wall of the tank as shown in Figure 4. These cross sections are taken away from the canyon and illustrates the di erent slopes used in the various canyon cases. . . . 8 6 Observed versus scaled plots for a) the depth of upwelling and b) the upwelling ux using the original scales from Allen and Hickey [2010]. The line of best t was found using all points but those corresponding to the double or at slope cases because those slopes were extreme; and in a) the high or low N cases were also unused because they were outliers. The norm of residuals is normalized to the x-axis. . . . . . . . . . . . . . . . . 11 7 Time series for the upwelling depth for the whole duration of the upwelling event. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 8 Time series for the volume of water upwelled onto the shelf for the whole duration of the upwelling event. . . . . . . . . . . . . . . . . . . . . . . . 14 v9 Isohalines over a) a at slope canyon under low f conditions, b) a at slope canyon under high f conditions, c) a double slope canyon under low f conditions, and d) a double slope canyon under high f conditions. The head depth, Hh; upwelling depth, Z; and head of the canyon are shown. The circled regions highlight the thickness of the upwelled water on the shelf. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16 10 A cross section of a double slope canyon with tilted isohalines. The key depth parameters are also shown: the head depth, Hh; the shelf break depth, Hs; and the depth of upwelling, Z (it should be noted that the Z was drawn a bit deeper to make it more visibly distinct from Hs). The length of the canyon, L, is also shown. The ow is coming out of the page as indicated by the . The gravitational force acting down the slope is balanced by the Coriolis force acting up slope. . . . . . . . . . . . . . . . 17 11 Observed versus scaled plots for a) a good linear relation between the scaled and observed upwelling ux using f=(sN) to account for the slope. b) shows the scaled and observed upwelling depth using the same non- dimensional number used to produce a linear ux plot. The lines of best t and the norm of residuals were found without using the at slope cases. The norm of residuals is normalized to the x-axis. . . . . . . . . . . . . . 22 12 Observed versus scaled plots for the upwelling depth using the non-dimensional number to account for the e ect of the shelf slope. a) Shows the result using a Rossy number based on the radius of curvature (Ro = U=fR) and b) shows the result using a Rossby number based on the canyon width (Ro = U=fW ). The norm of residuals is normalized to the x-axis. . . . . 23 13 Observed versus scaled plots for the upwelling ux using the non-dimensional number to account for the e ect of the shelf slope. a) Shows the result using a Rossy number based on the radius of curvature (Ro = U=fR) and b) shows the result using a Rossby number based on the canyon width (Ro = U=fW ). The norm of residuals is normalized to the x-axis. . . . . 24 vi14 Linearized observed and scaled upwelling depth plots where a) Shows the result using a Rossy number based on the radius of curvature (Ro = U=fR) and b) shows the result using a Rossby number based on the canyon width (Ro = U=fW ). The dotted line is a line of best t using all points weighted equally. The solid line is a weighted line of best t, where smaller error bars in Figure 12 are weighted more than those with larger error bars. The slopes of the weighted lines gives the values for a1 for each Rossby number. The norm of residuals is normalized to the y-axis. . . . . . . . . 25 15 Linearized observed and scaled upwelling ux plots where a) Shows the result using a Rossy number based on the radius of curvature (Ro = U=fR) and b) shows the result using a Rossby number based on the canyon width (Ro = U=fW ). The dotted line is a line of best t using all points weighted equally. The solid line is a weighted line of best t, where smaller error bars in Figure 13 are weighted more than those with larger error bars. The slopes of the weighted lines gives the values for a2 for each Rossby number. The norm of residuals is normalized to the y-axis. . . . . . . . . 26 viiList of Tables 1 Shows the values of the parameters varied. . . . . . . . . . . . . . . . . . 8 2 Values for the di erent cases with corresponding error. f is the Coriolis parameter, N is the Brunt-V ais al a frequency, s is the shelf slope, R is the radius of curvature of the canyon, W is the width of the canyon, U is the velocity upstream of the canyon, and Ro are Rossby numbers. Determination of the error will be explained later. . . . . . . . . . . . . . 9 3 These are the values for the non-dimensional number used in the scaling with corresponding error. . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 4 This table shows the values for the coe cients a1 and a2 with correspond- ing error for both the upwelling depth and upwelling ux, respectively, using scales based on the radius of curvature, R, or the canyon width, W . 27 viiiList of Appendices Initial Scaling 36 Depth of Upwelling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 Upwelling Flux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 ixAcknowledgments I would like to rst thank my supervisor Susan Allen, who provided much guidance along the way and made the time to meet with me every week. I would like to thank Jordan Dawe for showing me how to use the MITgcm model and Westgrid and also for his patience while I was learning. On that note I would also like to thank those from MIT who created and work on the MITgcm model which I used for my data; and to Westgrid which allowed me to run more cases than would have otherwise been possible on my own computer. And lastly, thank-you to the undergraduate thesis supervisors Erin Lane and Elspeth Barnes who supported each of us as we struggled to complete our theses. x1 Introduction Submarine canyons are found along the continental shelf and enhance the upwelling in that region. A modeling study by K ampf [2010] shows that the presence of submarine canyons in the eastern Great Australian Bight can increase the volume ux of water onto the shelf by a factor of 3.5. Upwelling is signi cant for bringing nutrient rich water from depth onto the continental shelf allowing these regions to be biologically active [Allen et al., 2001]. A scaling of upwelling over submarine canyons has been developed by Allen and Hickey [2010]. The quantities scaled included the depth of upwelling, the ux of water upwelled onto the shelf, and the deep water stretching. A number for determining the presence of a shelf break eddy has also been developed using a scaling method. A numerical investigation of a canyon in a laboratory tank, with a diameter of one meter, indicated that this scaling can represent the upwelling fairly accurately; however, there is still need for improvement. The deep water stretching scaling deviated most from the observed values and the model used in the numerical investigation produced no shelf break eddies; thus these scales will not be explored further in this paper. This paper will focus on investigating the in uence of the continental shelf slope on an upwelling event and investigate how to t it into the scaling of the upwelling depth and the ux of water onto the shelf. Motivation for examining the shelf slope came upon by looking into the scaling system developed by Allen and Hickey [2010] for the depth of upwelling and ux of water onto the shelf. Figure 1 shows that variations in the shelf slope produced results that deviated from the line of best t and thus does not follow the proposed scaling. It is important to include the shelf slope in the scaling because many continental regions have a slope and that will make this scaling more applicable to real canyon scenarios. In section 2 a background of the upwelling over submarine canyons will be discussed as well as the initial scaling system. Section 3 includes the methodology of the experiment and a brief description on the approach of adding the e ect of the shelf slope in the scaling. Section 4 will include the results of the numerical model and a comparison with the original scales by Allen and Hickey [2010]. Section 5 will include the new scaling approach in more detail. Section 6 will include a veri cation of the new scales using the data from MITgcm. Section 7 will be a discussion of the results and how the new scaling better represents the upwelling over submarine canyons. 10.2 0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 F1/2 RL 1/2 Upwelling Depth D hâˆ’ 1 a) Best Fit = 1.5 Â± 0.1 x + 0.04 Â± 0.06 Norm of Residuals = 0.5 Basic Small Big Big Width Small Width Big L Small L Big Wsb Small Wsb Big W Small W Big R Small R Double Slope Half Slope Low N High N Low f High f High U Low U Medium High U Medium Low U 0 2 4 6 8 10 12 14 16 18 20 0 20 40 60 80 100 120 140 160 180 U W sb Dh F 3/2 RL 1/2 (cm3sâˆ’1) Flux (c m3 sâˆ’ 1 ) b) Best Fit = 8 Â± 1 x + 11 Â± 8 Norm of Residuals = 17 Basic Small Big Big Width Small Width Big L Small L Big Wsb Small Wsb Big W Small W Big R Small R Double Slope Half Slope Low N High N Low f High f High U Low U Medium High U Medium Low U Figure 1: The original scaling system of Allen and Hickey [2010] for a) the depth of up- welling and b) the upwelling ux. The di erent dots represent the di erent canyon cases which varied the geometry of the canyon, the strati cation,N ; the Coriolis parameter, f ; and the U velocity upstream of the canyon. The variations in slope are shown by the red and green stars and are o the line of best t. The x in the best t refers to the x-axis. The norm of residuals is normalized to the x-axis. 22 Background 2.1 Upwelling over Submarine Canyons 2.1.1 Phases of Upwelling There are three phases to an upwelling event: an initial transient phase, a near steady advection dominated phase, and a relaxation phase [Allen and Durrieu de Madron, 2009]. The initial transient phase is a time dependent response as the ow along the shelf- break increases and can be explained by linear dynamics [Allen, 1996]. The advection dominated phase is a response to when the shelf-break ow is steady, this phase is non- linear and a system of scales will be examined further in this paper. The relaxation phase occurs when the shelf-break ow decreases. 2.1.2 Components of Flow Over Canyons Geostrophic ows have zero horizontal divergence and are forced to follow isobaths [Allen and Durrieu de Madron, 2009]. This being the case the ow must be ageostrophic in order to have upwelling. Some conditions for vertical water movement include large frictional forces, the importance of time dependence, or advection dominated ow [Allen, 2004]. Which of these actually contributes to the upwelling depends on the timing and the phase of upwelling. There are di erent layers to be considered when looking at canyon upwelling: the surface ow, the upwelling current, and the deep ow [Allen, 2004] (Figure 2). The surface ow goes right over the canyon because it is only weakly a ected by the canyon [Allen, 2004]. The rim depth eddy can be generated by the stretching of the water column as the ow crosses the canyon just above the canyon rim [Allen et al., 2001]. The upwelling current is water that ows into the canyon and is upwelled onto the shelf [Allen, 2004]. The surface ow crossing the canyon is in balance with a pressure gradient at the rim depth and drives the upwelling within the canyon [Freeland and Denman, 1982]. There is stronger upwelling closer to shelf-break depth and weaker upwelling lower in the canyon which causes a stretching of the water column at depth. The deep ow has a cyclonic vorticity due to this stretching of the water column [Allen, 2004]. 3Figure 2: Schematic of ow over a submarine canyon (from Allen and Hickey [2010]). 2.2 Initial Scaling of Upwelling Over Submarine Canyons 2.2.1 Assumptions of the Scaling A scaling system was developed by Allen and Hickey [2010]. Assumptions of the scaling include: weak to moderate, uniform incoming ow; uniform strati cation near the rim depth; a shallow shelf break depth; that the canyon is deep, narrow, has steep walls and does not cut too close to the coast; and that the canyon has a regular shape at the upstream corner. 2.2.2 Parameters in the Scaling The parameters used in the scaling for upwelling over submarine canyons include the along isobath velocity upstream of the canyon, U ; the Coriolis parameter, f ; the strati cation given by the Brunt-V ais al a buoyancy frequency, N ; the length of the canyon, L; the width of the canyon at half the canyon length, W ; the width of the canyon at the shelf break, Wsb; the radius of curvature of the canyon, R; and the slope of the continental shelf, s. Figure 3 shows the canyon geometry with some of these scales drawn on which are used in the scaling described below. 4Figure 3: The geometry of the canyon is shown. The white contours represent the 0.02 m shelf-break depth. Wsb is the shelf break width, W is the width of the canyon, L is the length of the canyon from the shelf break width to the head of the canyon, and R is the radius of curvature of the canyon. Other scales which are referred to are the depth to the head of the canyon, Hh, and the depth to the shelf break, Hs. Some non-dimensional numbers used in the scaling include two Rossby numbers, Ro = U=fR and U=fW , which determines whether the ow will follow the isobaths or cross over the canyon and drive upwelling; a function F which changes U velocities along isobaths to V velocities across isobaths, F = c1Ro=(c2 + Ro) where c1 and c2 are constants taken to be 1 and 0.9 respectively [Allen and Hickey, 2010]; another Rossby number, RL = U=fL; and a depth scale Dh = fL=N . 2.2.3 Scaling of the Upwelling Depth The depth of upwelling is the depth below the head at which the deepest water arrives on the shelf. The scaling for the upwelling depth was done by estimating the strength of the ow crossing the canyon and determining the resulting pressure gradient [Allen and Hickey, 2010]. From this the deformation of the density eld is found and by balancing the pressure gradients at the head and mouth of the canyon allows for the determination 5of the depth of upwelling [Allen and Hickey, 2010]. The scaling for the depth of upwelling is as follows: Z = fULF N2 = Dh(FRL) 1=2 (1) where Z is the upwelling depth, f is the Coriolis parameter, U is the velocity upstream of the canyon, L is the length of the canyon, F = c1Ro=(c2 + Ro) is a function used to relate velocities U and V (the radially outward velocity) as ow crosses the canyon where Ro is the Rossby number Ro = U=(fR), N is the Brunt-V ais al a frequency, Dh = fL=N is the depth scale, RL = U=(fL) is another Rossby number. The full derivation of this scaling by Allen and Hickey [2010] can be found in the appendix. 2.2.4 Scaling of the Upwelling Flux The upwelling ux is the ux of water upwelled onto the shelf. The scaling for the ux of water upwelled onto the shelf was accomplished by estimating the speed and vertical thickness of the upwelling stream [Allen and Hickey, 2010]. The scaling for the upwelling ux is as follows: = U WsbZ = (UWsbDh)F 3=2R1=2L (2) where is the upwelling ux, U is the velocity of the upwelling stream, and Wsb is the width of the canyon at the shelf break. The full derivation of this scaling by Allen and Hickey [2010] can also be found in the appendix. 3 Methodology 3.1 The Model Using the Massachusetts Institute of Technology general circulation model (MIT- gcm) [Marshall et al., 1997] numerous experiments of upwelling over canyons in a lab- oratory tank were simulated. This model was created to model the ocean, atmosphere and climate. It uses the Navier Stokes Equations and has the capability of being non- hydrostatic. The model can be used on both large and small scales and can be applied to this laboratory experiment. 63.2 The Experiment Setup The scenario modeled consists of a canyon in a laboratory tank as per Dawe and Allen [2010]. There is a canyon cutting into the continental shelf with the center consisting of a at abyss (Figure 4). A central cylinder was added to the topography to maintain numerical stability [Dawe and Allen, 2010]. x (cm) y (cm ) âˆ’50 âˆ’40 âˆ’30 âˆ’20 âˆ’10 0 10 20 30 40 50 âˆ’40 âˆ’30 âˆ’20 âˆ’10 0 10 20 30 40 âˆ’9 âˆ’8 âˆ’7 âˆ’6 âˆ’5 âˆ’4 âˆ’3 âˆ’2 âˆ’1 Figure 4: A contour plot of the Basic slope topography in a tank. The topography contours are in centimeters. The white contours show the 2 cm shelf break depth and in the center is the cylinder. The model was run using cylindrical coordinates. The model used a viscosity of 10 6 m2s 1 and a di usion coe cient of 1.5 10 9 m2s 1, both are molecular values for water. The model was run assuming that the system was non-hydrostatic to allow for internal waves, but the results were very similar to the hydrostatic case [Dawe and Allen, 2010]. Other speci ed values included a forcing velocity of 0.3 ms 1 which is owing clockwise in the tank and the topography of the tank. The di erent cases included in the experiment consisted of varying the strati cation in the tank by changing the Brunt-V ais al a frequency, N ; the rotation rate of the tank 7by varying the Coriolis parameter, f ; and by varying the slope of the continental shelf (Figure 5). A table consisting of the varied parameters are included in Table 1 and Table 2 includes the values of these parameters for the di erent cases along with the radius of curvature, R, and the canyon width, W , both of which are used for di erent calculations for the Rossby number, Ro; and the upstream velocity U . 0 5 10 15 20 25 30 35 40 45 50 âˆ’10 âˆ’8 âˆ’6 âˆ’4 âˆ’2 0 2 r (cm) z (cm ) Flat Slope Half Slope Small Slope Basic Slope Big Slope Double Slope Figure 5: Cross sections of the topography from the center to the wall of the tank as shown in Figure 4. These cross sections are taken away from the canyon and illustrates the di erent slopes used in the various canyon cases. Table 1: Shows the values of the parameters varied. Parameter Basic High Value Low Value N (s 1) 2.2 4.4 1.1 f (s 1) 0.52 0.62 0.42 Slope Basic Double Flat Slope Range Flat Half Small Basic Big Double 0.000 0.039 0.058 0.077 0.09 0.15 8T abl e 2: V alue s fo r th e di ere nt case s wit h corres p ondin g error . f is th e Corioli s par ameter , N is th e Bru nt-V ais al a frequenc y, s is th e shel f slo p e, R is th e radiu s of cur vatur e of th e ca ny on , W is th e widt h of th e ca ny on , U is th e vel oci ty upstrea m of th e ca ny on , an d R o ar e Ross by nu m b ers . Determinatio n of th e erro r wil l b e explaine d later . Cas e f (s 1 ) N (s 1 ) s R (cm ) W (c m ) U R o = U =f R R o = U =f W Basi c 0.5 2 2. 2 0.07 7 0.00 8 3. 8 0. 2 2. 69 0.0 3 1. 35 0.0 2 0.6 9 0.0 4 0.9 7 0.0 2 Basi c L ow N 0.5 2 1. 1 0.07 7 0.00 8 3. 8 0. 2 2. 69 0.0 3 1.2 9 0.0 2 0.6 6 0.0 4 0.9 2 0.0 2 Basi c Hig h N 0.5 2 4. 4 0.07 7 0.00 8 3. 8 0. 2 2. 69 0.0 3 1.4 3 0.0 2 0.7 3 0.0 4 1.0 2 0.0 2 Basi c L ow f 0.4 2 2. 2 0.0 77 0.00 8 3. 8 0. 2 2. 69 0.0 3 1.4 0 0.0 2 0.8 8 0.0 5 1.2 4 0.0 2 Basi c Hig h f 0.6 2 2. 2 0.07 7 0.00 8 3. 8 0. 2 2. 69 0.0 3 1.3 0 0.0 2 0.5 5 0.0 3 0.7 8 0.0 1 Doubl e Slo p e 0.5 2 2. 2 0.1 5 0.0 2 3. 8 0. 2 1.97 0 0.00 9 1.3 1 0.0 2 0.6 7 0.0 4 1.2 8 0.0 2 Doubl e Slo p e L ow N 0.5 2 1. 1 0.1 5 0.0 2 3. 8 0. 2 1.97 0 0.00 9 1.3 0 0.0 2 0.6 6 0.0 4 1.2 7 0.0 2 Doubl e Slo p e Hig h N 0.5 2 4. 4 0.1 5 0.0 2 3. 8 0. 2 1.97 0 0.00 9 1.3 3 0.0 2 0.6 8 0.0 4 1.3 0 0.0 2 Doubl e Slo p e L ow f 0.4 2 2. 2 0.1 5 0.0 2 3. 8 0. 2 1.97 0 0.00 9 1.3 4 0.0 2 0.8 4 0.0 5 1.6 2 0.0 3 Doubl e Slo p e Hig h f 0.6 2 2. 2 0.1 5 0.0 2 3. 8 0. 2 1.97 0 0.00 9 1.2 7 0.0 2 0.5 4 0.0 3 1.0 4 0.0 2 Hal f Slo p e 0.5 2 2. 2 0.03 9 0.00 4 4. 0 0. 2 3. 18 0.0 2 1.2 6 0.0 2 0.6 1 0.0 3 0.7 6 0.0 1 Hal f Slo p e L ow N 0.5 2 1. 1 0.03 9 0.00 4 4. 0 0. 2 3. 18 0.0 2 1.1 5 0.0 2 0.5 6 0.0 3 0.7 0 0.0 1 Hal f Slo p e Hig h N 0.5 2 4. 4 0.03 9 0.00 4 4. 0 0. 2 3. 18 0.0 2 1.3 5 0.0 2 0.6 6 0.0 3 0.8 2 0.0 1 Hal f Slo p e L ow f 0.4 2 2. 2 0.03 9 0.00 4 4. 0 0. 2 3. 18 0.0 2 1.3 3 0.0 2 0.7 8 0.0 4 1.0 0 0.0 2 Hal f Slo p e Hig h f 0.6 2 2. 2 0.03 9 0.00 4 4. 0 0. 2 3. 18 0.0 2 1.2 1 0.0 2 0.5 0 0.0 3 0.6 1 0.0 1 Fla t Slo p e 0.5 2 2. 2 0.00 0 0.00 2 2. 9 0. 2 2. 97 0.0 4 1.2 0 0.0 2 0.7 9 0.0 6 0.7 8 0.0 2 Fla t Slo p e L ow N 0.5 2 1. 1 0.00 0 0.00 2 2. 9 0. 2 2.9 7 0.0 4 1.0 9 0.0 2 0.7 2 0.0 5 0.7 1 0.0 2 Fla t Slo p e Hig h N 0. 52 4. 4 0.00 0 0.00 2 2. 9 0. 2 2.9 7 0.0 4 1.2 7 0.0 2 0.8 4 0.0 6 0.8 2 0.0 2 Fla t Slo p e L ow f 0.4 2 2. 2 0.00 0 0.00 2 2. 9 0. 2 2.9 7 0.0 4 1.2 8 0.0 2 1.0 4 0.0 7 1.0 3 0.0 2 Fla t Slo p e Hig h f 0.6 2 2. 2 0.0 00 0.00 2 2. 9 0. 2 2.9 7 0.0 4 1.1 3 0.0 2 0.6 2 0.0 4 0.6 1 0.0 1 Bi g Slo p e 0.5 2 2. 2 0.0 9 0.0 1 2. 9 0. 2 1.7 6 0.0 1 1.3 5 0.0 2 0.9 0 0.0 6 1.4 8 0.0 2 Bi g Slo p e L ow N 0.5 2 1. 1 0.0 9 0.0 1 2. 9 0. 2 1.7 6 0.0 1 1.3 1 0.0 2 0.8 7 0.0 6 1.4 3 0.0 2 Bi g Slo p e Hig h N 0.5 2 4. 4 0.0 9 0.0 1 2. 9 0. 2 1.7 6 0.0 1 1.4 3 0.0 2 0.9 5 0.0 7 1.5 6 0.0 3 Bi g Slo p e L ow f 0.4 2 2. 2 0.0 9 0.0 1 2. 9 0. 2 1.7 6 0.0 1 1.4 1 0.0 2 1.1 6 0.0 8 1.9 1 0.0 3 Bi g Slo p e Hig h f 0.6 2 2. 2 0.0 9 0.0 1 2. 9 0. 2 1.7 6 0.0 1 1.3 1 0.0 2 0.7 3 0.0 5 1.2 0 0.0 2 Smal l Slo p e 0.5 2 2. 2 0.05 8 0.00 6 3. 0 0. 2 1.9 9 0.0 1 1.3 1 0.0 2 0.8 5 0.0 6 1.2 7 0.0 2 Smal l Slo p e L ow N 0.5 2 1. 1 0.05 8 0.00 6 3. 0 0. 2 1.9 9 0.0 1 1.2 3 0.0 2 0.8 0 0.0 6 1.1 9 0.0 2 Smal l Slo p e Hig h N 0.5 2 4. 4 0.05 8 0.00 6 3. 0 0. 2 1.9 9 0.0 1 1. 41 0.0 2 0.9 2 0.0 6 1.3 6 0.0 2 Smal l Slo p e L ow f 0.4 2 2. 2 0.05 8 0.00 6 3. 0 0. 2 1.9 9 0.0 1 1. 37 0.0 2 1.1 1 0.0 8 1.6 4 0.0 3 Smal l Slo p e Hig h f 0.6 2 2. 2 0.05 8 0.00 6 3. 0 0. 2 1. 99 0.0 1 1. 26 0.0 2 0.6 9 0.0 5 1.0 2 0.0 2 93.3 Incorporating Shelf Slope in the Scaling As mentioned earlier the initial scaling did not include a variation of the shelf slope, but after investigation it appears to have an e ect on the upwelling. To account for the shelf slope in the scaling of the depth of upwelling and the ux of water upwelled onto the shelf a non-dimensional number will be introduced. The origin of this number will be determined by re-evaluating the derivation of the original scaling and trying to incorporate the shelf slope. It is likely that the non-dimensional number will be a function of the slope of the shelf, and other vertical and horizontal parameters such as the buoyancy frequency and Coriolis parameter, respectively. To test the scaling plots of observations from the MITgcm model versus scaled values for the depth of upwelling and upwelling ux will be made. The more linear the plot and less scatter means that the scaled quantities are close to the observations. The norm of residuals will be used to measure goodness of t of the scaled values to the observations. Error bars will be included to set a range of acceptable values; the method for determining error will be discussed later. 4 Initial Results 4.1 Scaling Simulations of upwelling events over submarine canyons in a one meter diameter tank were run using the MITgcm model. These simulations varied the shelf slope along with the strati cation and Coriolis parameter. Calculations for the depth of upwelling and the upwelling ux are made from the output data. The depth of upwelling is de ned as the di erence between the deepest water that arrives on the shelf at head depth and the head depth of the canyon. The upwelling ux is de ned as the volume of water that comes up through the canyon during the last period of upwelling, in this case, over the last ve seconds of the run time. Then using the canyon geometry and the velocity over the canyon scales for the upwelling depth and ux were found using the scales derived by Allen and Hickey [2010]. The observed and scaled values of the upwelling depth and ux are included in Figure 6. These plots are not as linear as those in Figure 1 and also have a larger norm of residuals indicating much scatter from the line of best t. Figure 6a shows that in general the normalized measured upwelling depth is greatest for the low Coriolis parameter cases followed by low strati cation cases, cases with only 100.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 F1/2 RL 1/2 Upwelling Depth D hâˆ’ 1 a) Best Fit = 1.0 Â± 0.1 x + 0.28 Â± 0.05 Norm of Residuals = 1.3 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f 0 1 2 3 4 5 6 7 8 9 10 0 50 100 150 200 250 300 U W sb Dh F 3/2 RL 1/2 (cm3sâˆ’1) Flux (c m3 sâˆ’ 1 ) b) Best Fit = 19 Â± 5 x âˆ’ 1e+01 Â± 3e+01 Norm of Residuals = 65 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f Figure 6: Observed versus scaled plots for a) the depth of upwelling and b) the upwelling ux using the original scales from Allen and Hickey [2010]. The line of best t was found using all points but those corresponding to the double or at slope cases because those slopes were extreme; and in a) the high or low N cases were also unused because they were outliers. The norm of residuals is normalized to the x-axis. 11variation in slope, high Coriolis parameter cases, and with the high strati cation cases having the shallowest normalized upwelling depth. Within each of these groupings the cases with smaller slopes tend to have the deepest upwelling depths and steeper slope cases tend to have lower upwelling depths. The low strati cation cases have very little variation even with di erent shelf slopes. The measured depth of upwelling for the high strati cation cases have much variation compared to the basic and low strati cation, or Coriolis parameter variation cases. The variations in the Coriolis parameter cases follow the line of best t better than the strati cation cases. Figure 6b shows that in general the measured ux is largest for low strati cation cases, lowest for high strati cation cases and everything else lies in between. Within each strati cation or Coriolis parameter grouping the atter slope cases have the largest ux and the cases with steeper slopes have the smaller uxes. Looking at the plot it can be seen that there is much scatter about the line of best t which is also indicated by a large norm of residuals. 4.2 Time Series 4.2.1 Depth of Upwelling Time series for the depth of upwelling were created to better understand what is happening during the upwelling event (Figure 7). The high strati cation cases have the lowest depth of upwelling and the low strati cation cases have the deepest depth of upwelling over the full 30 s. The other cases have less variation, but the low Coriolis parameter has a slightly smaller upwelling depth than the high Coriolis parameter case over the full 30 s. 120 5 10 15 20 25 30 0 0.5 1 1.5 2 2.5 3 Depth of upelling (cm ) Time (s) Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f Figure 7: Time series for the upwelling depth for the whole duration of the upwelling event. The variation between the di erent slopes for the low strati cation cases is much smaller than the variation between the di erent slopes for the high strati cation case. For the low strati cation cases the at slope had the deepest upwelling depth followed by the big slope, half slope, small slope, double slope, and with the basic slope having the shallowest upwelling depth at 30 s. For the high strati cation cases the atter slope cases had the deepest upwelling depths and steeper slope cases have shallower upwelling depths at 30 s. 4.2.2 Volume of Water Upwelled Time series for the volume of water onto the shelf were created to better understand what is happening during the upwelling event (Figure 8). The high strati cation cases have the smallest volume of water upwelled onto the shelf after 25 s and the low strati - cation cases have the largest amount of water upwelled onto the shelf after 10 s. Before these times the volume is very close and nearly indistinguishable between the other cases. The variations in the Coriolis parameter and the cases in which only the slope was varied cover about the same range of upwelled water volume at all times. 13The slope of the lines for the low strati cation cases is much greater than the slope of the line for the high strati cation cases. The slope of these lines between 25 s and 30 s gives the upwelling ux. 0 5 10 15 20 25 30 0 1000 2000 3000 4000 5000 6000 7000 8000 9000 Volume (c m3 ) Time (s) Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f Figure 8: Time series for the volume of water upwelled onto the shelf for the whole duration of the upwelling event. 5 Inclusion of Shelf Slope in the Scaling 5.1 Slope E ect on the Pressure Field As mentioned earlier submarine canyons enhance the upwelling onto the shelf. As a larger volume of denser water arrives on the shelf this has the potential to create a back pressure making further upwelling more di cult. If the slope of the shelf were at then the water upwelled would likely continue to ow up the shelf or be carried away. If the shelf slope were to have a steep slope then the head of the canyon would be at a shallower depth and thus shallower water would make it onto the shelf, but since the de nition of the upwelling depth is the di erence of the depth of the deepest water which reached the head and the head depth there is another factor which will cause this di erence to 14shorten. The isopycnals which reach the shelf will follow the slope of the shelf. If the slope is not at then the water that get upwelled will have a more di cult time moving up the shelf because there is dense water already there. This denser water will create a back pressure and restrict the movement of water up the shelf. To nd the baroclinic pressure gradient force at the upwelling depth, Z, use the hydrostatic pressure gradient: @p @z = g (3) where p is pressure, is the density, and g is the gravitational acceleration. Rearranging (3) to nd the change in pressure at the upwelling depth between the mouth and head of the canyon and adding the e ect of the slope to the pressure: p = gZ (s;N; f) (4) where p is the pressure di erence along depth Z +Hh where Z is the upwelling depth and Hh is the head depth; is the density di erence along depth Z + Hh and given by = @ =@zZ (s;N; f), where is the undisturbed horizontal density; and g is the gravitational acceleration. The function adds some back pressure due to the slope, where N is the Brunt-V ais al a frequency, s is the shelf slope, and f is the Coriolis parameter. The function in equation 4 acts to add some back pressure due to the slope and restrict the water from moving further up the shelf. It should be a function of the slope, s, because a greater slope should create a greater back pressure due to gravity acting on the dense water on the shelf; also a function of N because there is a stronger density gradient and with more dense water on the shelf creates more back pressure; and a function of f 1 because the bottom Ekman layer is de ned by: = p 2 f (5) where is the thickness of the Ekman layer and is the viscosity. Thus upwelled water on the shelf with a lower f will have a thicker layer and creates more back pressure. Figure 9 shows how water movement is restricted for steeper slopes and low f . The isohalines (which are taken to re ect the density) follow the slope of the shelf, so if the shelf slope were at there should be no back pressure and if the shelf slope were steep 15r (cm) z (cm ) a) 10 20 30 40 âˆ’10 âˆ’8 âˆ’6 âˆ’4 âˆ’2 0 10 20 30 40 50 60 r (cm) z (cm ) b) 10 20 30 40 âˆ’10 âˆ’8 âˆ’6 âˆ’4 âˆ’2 0 10 20 30 40 50 60 r (cm) z (cm ) c) 10 20 30 40 âˆ’10 âˆ’8 âˆ’6 âˆ’4 âˆ’2 0 10 20 30 40 50 60 r (cm) z (cm ) d) 10 20 30 40 âˆ’10 âˆ’8 âˆ’6 âˆ’4 âˆ’2 0 10 20 30 40 50 60 Hh Z Hh Hh Hh Z Z Z Surface SurfaceSurface Surface Head Head HeadHead Figure 9: Isohalines over a) a at slope canyon under low f conditions, b) a at slope canyon under high f conditions, c) a double slope canyon under low f conditions, and d) a double slope canyon under high f conditions. The head depth, Hh; upwelling depth, Z; and head of the canyon are shown. The circled regions highlight the thickness of the upwelled water on the shelf. then ow up the shelf will have the weight of the water further up the shelf acting against it and creating a larger pressure. For lower f the width of the upwelled water for a given isohaline is thicker compared to a higher f which is shown in the circles by an arrow between the topography and an isohaline. A thicker water mass that reaches the same horizontal extent from the head will have a larger mass and a greater weight associated with it, so this will create more back pressure. If s, N , and f 1 are all included then all components of the topographic Burger num- ber will be present. The topographic Burger number, sN=f , is found in many examples of ow along slopes [Brink and Lentz, 2010; Lentz and Chapman, 2004; Jacox and Edwards, 2011]. The topographic Burger number measures the importance of the bottom slope 16and so a large Burger number means that the buoyancy transport e ects are important to the problem [Brink and Lentz, 2010]. Another explanation for the origin of the Burger number, and back pressure in this case, can be related to a study by Nof [1982]. As the ow comes out of the page in Figure 10 there is a gravitational force acting to pull the water down the slope which is balanced by the Coriolis force pulling the water further up slope [Nof, 1982]. Writing this balance as: g0s = fu (6) Where g0 is the reduced gravity de ned by g0 = ( = )g, and u is the along-isobath velocity. Rearranging for u: u = g0s f (7) Since N = ( g= o)@ =@z and is a function of g there is a relation to the velocity of the ow and the topographic Burger number. Figure 10: A cross section of a double slope canyon with tilted isohalines. The key depth parameters are also shown: the head depth, Hh; the shelf break depth, Hs; and the depth of upwelling, Z (it should be noted that the Z was drawn a bit deeper to make it more visibly distinct from Hs). The length of the canyon, L, is also shown. The ow is coming out of the page as indicated by the . The gravitational force acting down the slope is balanced by the Coriolis force acting up slope. 17Here the gravitational force is also balanced by the Coriolis force which should have the u velocity to the left, but in this case the ow is to the right of the Coriolis force. The dense water on the shelf will want to ow the other way keeping the shallower water to the right of the ow [Nof, 1982]. This could result in the potential for the water to move down slope, back towards the canyon, and thus is a potential source of the back pressure due to the water mass on the shelf. Now that the parameters signi cant to creating the added back pressure characterized by has been established the change in pressure at Z +Hh becomes: p = g @ @z Z2 (s;N; f) (8) This is the change in pressure between the mouth and the head of the canyon. 5.2 Scaling for the Depth of Upwelling To determine the scale for the depth of upwelling, Z, follow the same steps by Allen and Hickey [2010] which are found in the appendix. The Baroclinic pressure gradient force is found by dividing (8) by L. Equating this to the Barotropic pressure gradient force above the canyon which is de ned as ofUF [Allen and Hickey, 2010]. Now rearranging the same way as Allen and Hickey [2010] the depth of upwelling scale becomes: Z = (s;N; f) 1=2Dh(FRL) 1=2 (9) where Z is the upwelling depth, f is the Coriolis parameter, U is the velocity upstream of the canyon, L is the length of the canyon, F is a function used to relate velocities U and V as ow crosses the canyon, N is the Brunt-V ais al a frequency, Dh = fL=N is the depth scale, RL = U=(fL) is a Rossby number, and is a non-dimensional number which accounts for the shelf slope. From experimentation with data from the MITgcm model to t a linear line for the observed and scaled values of the depth of upwelling it can be found that: 1 = exp asN f (10) where a is a coe cient which has included the 1/2 exponent and determines the steepness of the exponential, s is the slope of the shelf, N is the Brunt-V ais al a frequency, and f 18is the Coriolis parameter. Veri cation of this number will be shown later using the data from the MITgcm model. Thus, the depth of upwelling scale can be expressed as: Z = exp a1sN f Dh(FRL) 1=2 (11) 5.3 Scaling for the Upwelling Flux Using the same methodology as Allen and Hickey [2010] as outlined in the appendix, the scaling for the upwelling ux can be expressed as: = exp a2sN f (UWsbDh)F 3=2R1=2L (12) where is the upwelling ux, U is the velocity of the upwelling stream, and Wsb is the width of the canyon at the shelf break. Note that the coe cients a1 and a2 in the depth of upwelling and upwelling ux are di erent. Their determination will be explained in the next section. 6 Determination and Veri cation of 6.1 Early Attempts As shown above in Figures 1 and 6, the shelf slope clearly a ects both the depth of upwelling and the amount of water that upwells onto the shelf. Steeper slopes have a shallower upwelling depth and a smaller ux. These plots also do not have a very linear relation between the observed and the scaled values as there is much scatter about the line of best t. To improve the scaling and create a more linear plot with less scatter, a non-dimensional number will be determined and incorporated into the original scaling from Allen and Hickey [2010]. To start examining the in uence of the shelf slope in the upwelling depth and the upwelling ux a hypothesis was made that the role of the shelf slope would have to do with the shelf slope, the buoyancy frequency and the Coriolis parameter. To test this hypothesis variations of the above were included into the original scaling and plotted against the measured values. The more linear the plot and the smaller the norm of residuals the better the scaling is at representing the modelled results. 19The most linear plot was given by the ux with the non-dimensional number repre- senting the slope as f=(sN) where s is the shelf slope (Figure 11a). However, did not work with the depth of upwelling(Figure 11b). Also because there was a division by the shelf slope the at slope cases went to in nity. 6.2 Later Results To correct for the at slope cases going to in nity, an exponential was used, making the new non-dimensional number 1 = exp(aNs=f) (Figures 12 and 13). The values for these non-dimensional numbers can be found in Table 3. Previous work has used the Rossby number as a function of the radius of curvature [Allen and Hickey, 2010], but here both a Rossby number as a function of the radius of curvature, R and as a function of the canyon width, W are used to see which scales better. The signi cance of this is whether W or R determines if the ow will follow the isobaths or cross the canyon and drive upwelling [Allen and Hickey, 2010]. These gures (12 and 13) have a coe cient a1 for the depth of upwelling and a2 for the upwelling ux, which is determined by nding the slope of the weighted least-squares t to the linearized plots found in Figures 14 and 15 where those points with the smallest error associated with them have the highest weighting, and in some cases masking out a highly weighted outlier. Values for these coe cients are found in Table 4. 6.3 Error Analysis Error bars are added to the plots to get a better understanding of how well the scaling ts the observations. If the line of best t goes through the error bars of most of points then the scaling is successful. The error bars are found by using standard fractional error techniques [e.g. Bevington, 1969]. 6.3.1 Error in the Depth of Upwelling The fractional error in the y-axis on the depth of upwelling plots in Figures 1a, 6a, and 12 is: y y = s 2Z Z2 + 2Dh D2h (13) 20T abl e 3: Thes e ar e th e value s fo r th e non-dimensiona l nu m b er use d in th e scalin g wit h corres p ondin g error . Cas e sN =f 1 Z (R ) 1 Z (W ) 1 f lu x( R ) 1 f lu x( W ) Basi c 0.3 3 0.0 4 0.8 5 0.0 2 1.2 4 0.0 4 0.4 2 0.0 5 0.3 6 0.0 5 Basi c L ow N 0.1 6 0.0 2 0.9 2 0.0 1 1.1 1 0.0 2 0.6 5 0.0 4 0.6 0 0.0 4 Basi c Hig h N 0.6 5 0.0 7 0.7 2 0.0 4 1.5 4 0.0 9 0.1 7 0.0 4 0.1 3 0.0 3 Basi c L ow f 0.4 0 0.0 4 0.8 2 0.0 3 1.3 1 0.0 5 0.3 4 0.0 5 0.2 8 0.0 5 Basi c Hig h f 0.2 7 0.0 3 0.8 7 0.0 2 1.2 0 0.0 3 0.4 8 0.0 5 0.4 2 0.0 5 Doubl e Slo p e 0.6 5 0.0 7 0.7 2 0.0 4 1.5 4 0.0 9 0.1 7 0.0 4 0.1 3 0.0 3 Doubl e Slo p e L ow N 0. 33 0.0 3 0.8 5 0.0 2 1.2 4 0.0 4 0.4 2 0.0 5 0.3 6 0.0 5 Doubl e Slo p e Hig h N 1. 3 0. 1 0.5 2 0.0 6 2. 4 0. 3 0.0 3 0.0 1 0.01 6 0.00 8 Doubl e Slo p e L ow f 0.8 1 0.0 9 0.6 7 0.0 5 1. 7 0. 1 0.1 2 0.0 3 0.07 8 0.0 3 Doubl e Slo p e Hig h f 0.5 5 0.0 6 0.7 6 0.0 4 1.4 3 0.0 7 0.2 3 0.0 5 0.1 8 0.0 4 Hal f Slo p e 0.1 6 0.0 2 0.9 2 0.0 1 1.1 1 0.0 2 0.6 5 0.0 4 0.6 0 0.0 4 Hal f Slo p e L ow N 0.08 1 0.00 9 0.96 0 0.00 7 1.0 6 0.00 8 0.8 0 0.0 2 0.7 7 0.0 3 Hal f Slo p e Hig h N 0.3 3 0.0 4 0.8 5 0.0 2 1.2 4 0.0 4 0.4 2 0.0 5 0.3 5 0.0 5 Hal f Slo p e L ow f 0.2 0 0.0 2 0.9 0 0.0 2 1.1 4 0.0 2 0.5 8 0.0 4 0.5 3 0.0 5 Hal f Slo p e Hig h f 0.1 4 0.0 2 0.9 3 0.0 1 1.0 9 0.0 1 0.6 9 0.0 4 0.6 5 0.0 4 Fla t Slo p e 0.00 0 0.00 7 1.00 0 0.00 3 1.00 0 0.00 4 1.0 0 0.0 2 1.0 0 0.0 2 Fla t Slo p e L ow N 0.00 0 0.00 3 1.00 0 0.00 2 1.00 0 0.00 2 1.00 0 0.00 9 1.0 0 0.0 1 Fla t Slo p e Hig h N 0.0 0 0.0 1 1.00 0 0.00 7 1.00 0 0.00 9 1.0 0 0.0 4 1.0 0 0.0 4 Fla t Slo p e L ow f 0.00 0 0.00 8 1.00 0 0.00 4 1.00 0 0.00 6 1.0 0 0.0 2 1.0 0 0.0 3 Fla t Slo p e Hig h f 0.00 0 0.00 6 1.00 0 0.00 3 1.00 0 0.00 4 1.0 0 0.0 2 1.0 0 0.0 2 Bi g Slo p e 0.3 8 0.0 4 0.8 2 0.0 3 1.2 9 0.0 4 0.3 6 0.0 5 0.3 0 0.0 5 Bi g Slo p e L ow N 0.1 9 0.0 2 0.9 1 0.0 1 1.1 3 0.0 2 0.6 0 0.0 4 0.5 4 0.0 4 Bi g Slo p e Hig h N 0. 77 0.0 8 0.6 8 0.0 4 1. 7 0. 1 0.1 3 0.0 4 0.0 9 0.0 3 Bi g Slo p e L ow f 0.4 7 0.0 5 0.7 9 0.0 3 1.3 7 0.0 6 0.2 8 0.0 5 0.2 2 0.0 4 Bi g Slo p e Hig h f 0.3 2 0.0 3 0.8 5 0.0 2 1.2 4 0.0 4 0.4 2 0.0 5 0.3 6 0.0 5 Smal l Slo p e 0.2 4 0.0 3 0.8 8 0.0 2 1.1 8 0.0 3 0.5 2 0.0 5 0.4 6 0.0 5 Smal l Slo p e L ow N 0.1 2 0.0 1 0.9 4 0.0 1 1.0 8 0.0 1 0.7 2 0.0 3 0.6 8 0.0 3 Smal l Slo p e Hig h N 0.4 9 0.0 5 0.7 8 0.0 3 1.3 8 0.0 6 0.2 7 0.0 5 0.2 1 0.0 4 Smal l Slo p e L ow f 0.3 0 0.0 3 0.8 6 0.0 2 1.2 2 0.0 3 0.4 4 0.0 5 0.3 8 0.0 5 Smal l Slo p e Hig h f 0.2 0 0.0 2 0.9 0 0.0 2 1.1 4 0.0 2 0.5 8 0.0 4 0.5 2 0.0 4 210 10 20 30 40 50 60 70 80 0 50 100 150 200 250 300 (f sâˆ’1 Nâˆ’1)U W sb Dh F 3/2 RL 1/2 (cm3sâˆ’1) Flux (c m3 sâˆ’ 1 ) a) Best Fit = 3.0 Â± 0.2 x + 9 Â± 4 Norm of Residuals = 4 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 (f sâˆ’1 Nâˆ’1) F1/2 RL 1/2 Upwelling Depth D hâˆ’ 1 b) Best Fit = 0.06 Â± 0.02 x + 0.57 Â± 0.03 Norm of Residuals = 0.28 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f Figure 11: Observed versus scaled plots for a) a good linear relation between the scaled and observed upwelling ux using f=(sN) to account for the slope. b) shows the scaled and observed upwelling depth using the same non-dimensional number used to produce a linear ux plot. The lines of best t and the norm of residuals were found without using the at slope cases. The norm of residuals is normalized to the x-axis. 220.15 0.2 0.25 0.3 0.35 0.4 0.45 0.2 0.3 0.4 0.5 0.6 0.7 0.8 exp(a1 s N f âˆ’1)(F RL) 1/2 Upwelling Depth D hâˆ’ 1 a) Best Fit = 1.3 Â± 0.2 x + 0.24 Â± 0.06 Norm of Residuals = 0.81 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.2 0.3 0.4 0.5 0.6 0.7 0.8 exp(a1 s N f âˆ’1)(F RL) 1/2 Upwelling Depth D hâˆ’ 1 b) Best Fit = 1.4 Â± 0.1 x + 0.19 Â± 0.05 Norm of Residuals = 0.64 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f Figure 12: Observed versus scaled plots for the upwelling depth using the non-dimensional number to account for the e ect of the shelf slope. a) Shows the result using a Rossy number based on the radius of curvature (Ro = U=fR) and b) shows the result using a Rossby number based on the canyon width (Ro = U=fW ). The norm of residuals is normalized to the x-axis. 230 1 2 3 4 5 6 0 50 100 150 200 250 300 exp(a2 s N f âˆ’1) U W sb Dh F 3/2 RL 1/2 (cm3 sâˆ’1) Flux (c m3 sâˆ’ 1 ) a) Best Fit = 32 Â± 3 x âˆ’ 1 Â± 9 Norm of Residuals = 60 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f 0 1 2 3 4 5 6 7 8 0 50 100 150 200 250 300 exp(a2 s N f âˆ’1) U W sb Dh F 3/2 RL 1/2 (cm3 sâˆ’1) Flux (c m3 sâˆ’ 1 ) b) Best Fit = 25 Â± 3 x + 7 Â± 9 Norm of Residuals = 62 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f Figure 13: Observed versus scaled plots for the upwelling ux using the non-dimensional number to account for the e ect of the shelf slope. a) Shows the result using a Rossy number based on the radius of curvature (Ro = U=fR) and b) shows the result using a Rossby number based on the canyon width (Ro = U=fW ). The norm of residuals is normalized to the x-axis. 240 0.5 1 1.5 âˆ’0.2 âˆ’0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 s N fâˆ’1 log(Upwelling Depth ( D h F R L)âˆ’ 1 ) a) Best Fit = âˆ’0.50 Â± 0.07 x + 0.69 Â± 0.02 Norm of Residuals = 0.91 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f 0 0.5 1 1.5 âˆ’0.4 âˆ’0.2 0 0.2 0.4 0.6 0.8 s N fâˆ’1 log(Upwelling Depth ( D h F R L)âˆ’ 1 ) b) Best Fit = âˆ’0.66 x Â± 0.05 + 0.65 Â± 0.02 Norm of Residuals = 0.83 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f Figure 14: Linearized observed and scaled upwelling depth plots where a) Shows the result using a Rossy number based on the radius of curvature (Ro = U=fR) and b) shows the result using a Rossby number based on the canyon width (Ro = U=fW ). The dotted line is a line of best t using all points weighted equally. The solid line is a weighted line of best t, where smaller error bars in Figure 12 are weighted more than those with larger error bars. The slopes of the weighted lines gives the values for a1 for each Rossby number. The norm of residuals is normalized to the y-axis. 250 0.5 1 1.5 1 1.5 2 2.5 3 3.5 4 s N fâˆ’1 log(Flux (U W s b D h F3/ 2 R L1/ 2 )âˆ’ 1 ) a) Best Fit = âˆ’2.7 Â± 0.2 x + 3.41 Â± 0.08 Norm of Residuals = 0.8 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f 0 0.5 1 1.5 0.5 1 1.5 2 2.5 3 3.5 4 s N fâˆ’1 log(Flux (U W s b D h F3/ 2 R L1/ 2 )âˆ’ 1 ) b) Best Fit = âˆ’3.2 x Â± 0.2 + 3.31 Â± 0.08 Norm of Residuals = 1.0 Basic Basic Low N Basic High N Basic Low f Basic High f Double Slope Double Slope Low N Double Slope High N Double Slope Low f Double Slope High f Half Slope Half Slope Low N Half Slope High N Half Slope Low f Half Slope High f Flat Slope Flat Slope Low N Flat Slope High N Flat Slope Low f Flat Slope High f Big Slope Big Slope Low N Big Slope High N Big Slope Low f Big Slope High f Small Slope Small Slope Low N Small Slope High N Small Slope Low f Small Slope High f Figure 15: Linearized observed and scaled upwelling ux plots where a) Shows the result using a Rossy number based on the radius of curvature (Ro = U=fR) and b) shows the result using a Rossby number based on the canyon width (Ro = U=fW ). The dotted line is a line of best t using all points weighted equally. The solid line is a weighted line of best t, where smaller error bars in Figure 13 are weighted more than those with larger error bars. The slopes of the weighted lines gives the values for a2 for each Rossby number. The norm of residuals is normalized to the y-axis. 26Table 4: This table shows the values for the coe cients a1 and a2 with corresponding error for both the upwelling depth and upwelling ux, respectively, using scales based on the radius of curvature, R, or the canyon width, W . a1 a2 Ro(R) -0.50 0.07 -2.7 0.2 Ro(W ) -0.66 0.05 -3.2 0.2 Where y = Z=Dh, Z is the depth of upwelling and Dh is a depth scale de ned earlier. is referring to the error in each of its subscripts. Since Z was found by taking the di erence between the deepest water on the shelf at head depth (Z + Hh)and the head depth (Hh) the error in the measured depth of upwelling was determined by: Z = q 2Hh + 2 Hh+Z (14) Where Hh and Hh+Z are the errors in the z direction. These are found by taking half of the z grid spacing de ned in the model. Error in the depth scale Dh was calculated as: Dh = f N L (15) Where L is the error in the length of the canyon. This was determined by taking half of the di erence in the allowable canyon lengths determined from the minimum and maximum acceptable head depths. The fractional error in the scaled normalized upwelling depth on the x-axes (FRL)1=2 for Figures 1a and 6a was determined to be: x x = s 2F 4F 2 + 2RL 4R2L (16) Where x = (FRL)1=2, F is a function of Ro (c1Ro=(c2+Ro)), and RL is a Rossby number. But for Figure 12 where the x-axis is exp(a1sN=f)(FRL)1=2 the fractional error was calculated as: x x = s 2F 4F 2 + 2RL 4R2L + 2a1N 2s2 f 2 + 2sN2a 2 1 f 2 (17) Where x = exp(asN=f)(FRL)1=2, s is the shelf slope, and a1 is a coe cient described 27earlier. The fractional error in F was calculated as: F F = Ro Ro c2 +Ro 1 c2 +Ro (18) Where c2 = 0:9 and is a constant [Allen and Hickey, 2010]. The fractional error in Ro was calculated as: Ro Ro = r 2U U2 + 2R R2 (19) Where U is the velocity upstream of the canyon and R is the radius of curvature of the canyon. The error in U was found by taking half of the largest di erence of the velocity at three locations upstream of the canyon at a depth of 1 cm; this value was found to be 0.02 cm/s. Since R was determined by using the Levenberg-Marquandt algorithm for non-linear least-squares curve tting, the error in R was determined by the di erence in the solutions by using di erent lengths for which the solution was found. The di erence between the solutions was approximately 0.2 cm. In the cases where the Rossby number was found using the width of the canyon rather than the radius of curvature, W would be used instead of R. The error in W was determined by taking half the di erence in the widths for the maximum and minimum allowable canyon lengths. Error in the shelf slope, s, was calculated as: s = s 2 z 1 r 2 + 2 r z r 2 (20) Where z is the di erence between the depth at the wall of the tank and the depth near the shelf break at 1.7 cm where the shelf has a linear slope within error, r is the di erence in radius between the same points as z. The error z and r are found by p 2 2z and p 2 2z respectively where the error in z and r is half of the grid spacing. Error in a1 was returned in the least-squares solution to the slope of the weighted line of best t in Figure 14. 286.3.2 Error in the Upwelling Flux The error for the upwelling ux calculated from the model on the y-axis of Fig- ures 1b, 6b, and 13 is found by calculating the ux between 29 s and 30 s and taking the di erence from the ux value used over 25 s to 30 s. The fractional error in the scaled ux used in Figures 1b and 6b was calculated as: = s 2U U2 + 2Wsb W 2sb + 2Dh D2h + 9 2F 4F 2 + 2RL 4RL 2 (21) Where is the upwelling ux and Wsb is the width of the canyon at the shelf-break. The error in Wsb is found by taking half of the di erence between the maximum and minimum acceptable values. The fractional error for the upwelling ux using = exp(a2Ns=f)(UWsbDh)F 3=2R 1=2 L was calculated as: = s 2U U2 + 2Wsb W 2sb + 2Dh D2h + 9 2F 4F 2 + 2RL 4RL 2 + 2a2N 2s2 f 2 + 2sN2a 2 2 f 2 (22) Where a2 is the analogous to the a1 coe cient but for the upwelling ux rather than the depth of upwelling. 6.3.3 Error in Topographic Burger Number and 1 Error in the topographic Burger number (sN=f) was calculated as: sN=f = s N f (23) The error in the non-dimensional number 1 = exp (asN=f) included in the scaling was calculated as: 1 = s 2aN2s2 f 2 exp 2asN f + 2sa2N2 f 2 exp 2asN f (24) Where the a and itâ€™s corresponding error can be either a1 or a2. 297 Discussion 7.1 Assumptions An assumption for the initial scaling by Allen and Hickey [2010] is that the slope is slightly sloped. This thesis tried to remove this assumption and allow the scaling of the upwelling depth and ux to be applicable to variations in shelf slope. A non-dimensional number, , was found that characterized the back pressure which becomes present with increasing slope. This number when incorporated into the scaling and rearranged acts to shallow the depth of upwelling and lower the upwelling ux with increasing slope by a factor of exp (asN=f). There are some bene ts to using an exponential: this also allows the slope to be used in such a way that if the slope were at (equal to zero) then the term would be equal to one and there would be no added back pressure to the upwelling ow. And if there were any slope would increase representing the added back pressure to upwelling ow and restricting its movement further up the shelf. An assumption here is that the slope of the shelf is positive, further tests would have to be done to see what would happen if the slope declined after the head of the canyon towards the coast. The dense water will likely ow down the slope which may pull the water behind it and acting to reduce the acting pressure. This new scaling will still follow other assumptions by Allen and Hickey [2010] de- scribed earlier. 7.2 Slope E ect on Upwelling Characteristics 7.2.1 Depth of Upwelling The depth of upwelling is greatly a ected by the slope which can be seen by the separation of the points along the observed axis in Figure 6a for cases with the same slope, Coriolis, and most strati cation cases. For example: with a Coriolis parameter of f=0.52 s 1 and a buoyancy frequency of N=2.2 s 1 a at slope will have a greater upwelling depth and ux than a double slope case. It can be seen that the strati cation can be low enough that the slope has no or little e ect on the depth of upwelling as all low strati cation cases have approximately the same depth of upwelling. The low strati cation will have a large upwelling depth that is independent of the shelf slope. This can be thought of the less dense water on the shelf is creating little 30back pressure and thus there is little restriction of the upwelling. However, this is not observed in the scaled values (Figure 12) where the slope e ect on strati cation is clearly observed. The high strati cation cases on the other hand have a steeper slope than the line of best t in Figure 6a which shows that the slope e ect is greater than the strati cation e ect for denser water. The combination of a steep slope and very dense water will produce a larger back pressure and further restrict the water from being upwelled onto the shelf. 7.2.2 Upwelling Flux In comparison to the depth of upwelling the upwelling ux seems to have the opposite e ect where the low strati cation cases have a greater variation in ux with slope and the high strati cation cases have a smaller variation with changing slope. This could be due to the back pressure of the water already on the shelf. The slopes of the lines in Figure 8 for volume of water on the shelf are much steeper for the low strati cation cases than the high strati cation cases and so the ux will also be higher. Within the slope cases the slope of lines vary more for the low strati cation and thus there will be a greater variation in ux. The slopes of the lines for the high strati cation cases are more similar and thus the ux for each of these will be more similar. Less dense water will result in a lower back pressure because the weight of the water on the shelf would be less. This makes it easier for the water push past itself with a steeper slope allowing the volume of water making it onto the shelf in the last ve seconds of the model run time higher than what would have been present at earlier times. In this case the weight of the water has little e ect on restricting movement further up the slope and so it is the shelf slope that is the key factor in limiting the rate at which the water arrives on the shelf. The more dense the water the larger back pressure because the weight of the water on the shelf is larger. With a stronger strati cation the heavier water has a greater momentum to push past the less dense water already on the shelf; therefore, changes in the shelf slope have little e ect on the rate that the dense water can make it onto the shelf and the uxes would be similar. 317.3 Scales 7.3.1 Depth of Upwelling Scales The scaling for the depth of upwelling with the included non-dimensional term to account for the shelf slope shows a good relationship with the observed upwelling depth. This can be seen by the low norm of residuals (Figure 12). These gures also show that the scale works best using a Rossby number as a function of the canyon width, W , because it has a lower norm of residuals compared to the scale using a Rossby number as a function of the radius of curvature, R. This means the width of the canyon determines whether ow crosses the canyon and drives upwelling rather than following the isobaths. The two results have a percent di erence of about 23% in the norm of residuals. This is a di erent relation than that posed by Allen and Hickey [2010], who say the radius of curvature, R, will determine whether there is upwelling. 7.3.2 Upwelling Flux Scales The scaling for the upwelling ux with the included non-dimensional term to account for the shelf slope shows a good relationship with the observed upwelling ux; however, the low strati cation cases are still very scattered. This can be seen by the low norm of residuals (Figure 13). These gures also show that the scale works best using a Rossby number as a function of the radius of curvature, R, because it has a lower norm of residuals compared to the scale using a Rossby number as a function of the canyon width, W . This means that the radius of curvature, R, determines whether the ow will cross the canyon rather than follow the isobaths. The two results have a percent di erence of about 3% in the norm of residuals. This agrees with Allen and Hickey [2010] but not with the depth of upwelling in this study. Since the percent di erence is very small in this case more work would have to be done either by re-analyzing the statistics or looking at more canyons to con rm whether it is R or W that determines if the ow will cross the canyon. Since the two results are very close it is possible that the canyon width could be that determining factor. 8 Conclusions A previous investigation into the accuracy of the scaling developed by Allen and Hickey [2010] showed that the shelf slope has an impact on the upwelling depth and ux 32for ow over submarine canyons. This thesis investigated how the shelf slope plays a role in the upwelling over submarine canyons by creating a non-dimensional number which takes into account the shelf slope. This number is incorporated into the scaling by acting as a back pressure to the ow reaching the shelf. As water piles up on the shelf this creates a back pressure restricting water from moving up the shelf. The e ect of the back pressure on the depth of upwelling and upwelling ux was found to be exp(asN=f). A comparison of Rossby numbers used in the scaling show that the radius of curvature, R, has a better t with the data for the depth of upwelling, but the canyon width, W has a better t with the data for the upwelling ux. Future research will include looking for a more physical explanation of why an expo- nential term works for the scaling or whether it is some other function of sN=f . Also further research should look into whether the canyon width or radius of curvature de- termines whether the ow will follow the isobaths or cross the canyon and drive the upwelling. This was only one experiment and the two theories of 1) exp(asN=f) char- acterizes the back pressure and 2) that it is likely the width of the canyon determines whether the ow follows the isobaths or crosses them to ow over the canyon and drive upwelling should be tested with other canyon geometries and with real-world canyons. 33References Allen, S. E., 1996: Topographically generated, subinertial ows within a nite length canyon. J. Phys. Oceanogr., 26, 1608{1632. Allen, S. E., 2004: Restrictions on deep ow across the shelf-break and the role of submarine canyons in facilitating such ow. Surv. Geophys., 25, 221{247. Allen, S. E. and X. Durrieu de Madron, 2009: A review of the role of submarine canyons in deep-ocean exchange with the shelf. Ocean Sci., 5, 607{620. Allen, S. E. and B. M. Hickey, 2010: Dynamics of advection-driven upwelling over a shelf break submarine canyon. J. Geophys. Res., 115, C08 018. Allen, S. E., C. Vindeirinho, R. E. Thomson, M. G. G. Foreman, and D. L. Mackas, 2001: Physical and biological processes over a submarine canyon during an upwelling event. Can. J. Fish. Aquat. Sci., 58, 671{684. Bevington, P. R., 1969: Data Reduction and Error Analysis for the Physical Sciences. McGraw-Hill Book Company, 56-65 pp. Brink, K. H. and S. J. Lentz, 2010: Buoyancy arrest and bottom ekman transport. part I: Steady ow. J. Phys. Oceanogr., 40, 621{635. Dawe, J. T. and S. E. Allen, 2010: Solution convergence of ow over steep topography in a numerical model of canyon upwelling. J. Geophys. Res., 115, C05 008. Freeland, H. J. and K. L. Denman, 1982: A topographically controlled upwelling center o southern vancouver island. J. Mar. Res., 40, 1069{1093. Holton, J. R., 1992: An Introduction to Dynamic Meteorology. 3d ed., Academic Press, 511 pp. Jacox, M. G. and C. A. Edwards, 2011: E ects of strati cation and shelf slope on nutrient supply in coastal upwelling regions. J. Geophys. Res., 116, C03 019. K ampf, J., 2010: On preconditioning of coastal upwelling in the eastern great australian bight. J. Geophys. Res., 115, C12 071. Lentz, S. J. and D. C. Chapman, 2004: The importance of nonlinear cross-shelf momen- tum ux during wind-driven coastal upwelling. J. Phys. Oceanogr., 34, 2444{2457. 34Marshall, J., A. Adcroft, C. Hill, L. Perelman, and C. Heisey, 1997: A nite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers. J. Geophys. Res., 102, 5753{5766. Nof, D., 1982: The translation of isolated cold eddies on a sloping bottom. Deep-Sea Research, 30, 171{182. 35Appendices A Initial Scaling This is the scaling for upwelling over submarine canyons developed by Allen and Hickey [2010]. The scaling is based on the advection-driven upwelling and is assumed to be steady. The shallow water equations for steady, inviscid, strati ed ow with Boussinesq, incom- pressible and non-di usive approximations are: !u r !u fk^ !u = 1 o rp (25a) r !u = 0 (25b) !u r = 0 (25c) @p @z = g (25d) where !u is the horizontal velocity, f is the Coriolis parameter, k^ is the vertical unit vector, o is a constant reference density, p is the pressure, and g is the gravitational acceleration. A.1 Depth of Upwelling Using equation (25a) in polar coordinates it is possible to nd the pressure change along the canyon that drives the upwelling. ur @u @r + u r @u @ + u ur r + fur = 1 or @p @ (26) The rst term of this equation is small compared to the third term and can be neglected. By conservation of volume the second term can be written as: u 1 r @u @ = u @w @z @ur @r (27) Where @w=@z is small in comparison to @ur=@r and will also be neglected. The second term of equation (26) can then be scaled as UV=R, where R is the radius of curvature of 36the canyon (Figure 3). The third term can be scaled as UV=R, fur can be scaled as fV . A function F is de ned as a function of the Rossby number (Ro = U=fR) which turns along-isobath ow into cross-isobath ow, so V = F (Ro)U and F (Ro) = c1Ro c2 +Ro (28) Where c1 and c2 are constants. Now rearranging the scaled equation (26) becomes: ofUF = 1 r @p @ (29) This is the pressure gradient along the upstream rim edge of the canyon. The pressure di erence between the head and the mouth of the canyon (equation 30) is the pressure gradient multiplied by a length scale, L. Here the length scale is de ned as the distance between the mouth at the shelf-break depth and the head of the canyon which is determined by the last dimpled isobath of the canyon (L in Figure 3). ofUFL = L @p @ (30) Now to determine the deformation of the density eld, which will allow the depth of upwelling to be determined, equations (25) are rewritten using the natural coordinate system [e.g. Holton, 1992]: u @u @s = 1 o @p @s (31a) u2 R0 + fu = 1 o @p @n (31b) @u @s + @w @z = 0 (31c) u @p @s + w @p @z = 0 (31d) @p @z = g (31e) Where u is the velocity following the horizontal direction of ow, s; w is the vertical velocity; is the density eld; and R0 is the radius of curvature of the streamlines. The deepest streamline to cross the rim of the canyon is found to be at a depth Z+Hs where Z is the upwelling depth and Hs is the depth of the shelf break. Along this 37streamline the density is conserved, which is given in equation (31d). Scaling equation (31d) gives: U L oN2 g = 0 (32) Where U is the velocity of the upwelling stream, is the scale of horizontal density perturbations, L is a length scale from the shelf-break to the head of the canyon, and is the vertical velocity scale. Scaling the conservation of mass equation (31c) gives: U L Z = 0 (33) Where Z is the depth of upwelling. Since the pressure change along the canyon that drives upwelling is in balance with the density gradient due to the perturbation density equation (31e) will be scaled as: ofUFL Z = g (34) To nd the scale for the depth of upwelling equations (32), (33), and (34) will be combined: Z = Dh(FRL) 1=2 (35) Where RL is a Rossby number given by RL = U=fL and Dh is a depth scale given by Dh = fL=N . A.2 Upwelling Flux The scale for upwelling ux can be found by scaling equation (31b). This ow is found directly below the shelf-break ow crossing the canyon and has a scaled pressure gradient of fUF . At the lowest part of the upwelling stream at Z there is no pressure gradient. Between the shelf-break and Z the pressure gradient is half of that at the shelf-break, so the scaled pressure gradient of the upwelling stream is given by (1=2)fUF . U2 Wsb + fU = 1 2 fUF (36) Where Wsb is the shelf-break width of the canyon and is the turning radius of the ow into the canyon. The rst term is small and can be neglected [Allen and Hickey, 2010], 38doing this equation (36) becomes: U = 1 2 UF (37) The ux is found by, = U multiplied by the area of upwelling stream. The area is found by the multiplication of the depth of upwelling, Z and the shelf-break width, Wsb. = U ZWsb = (UWsbDH)F 3=2R1=2L (38) 39
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Undergraduate Research /
- A Numerical Investigation of the Shelf Slope in the...
Open Collections
UBC Undergraduate Research
A Numerical Investigation of the Shelf Slope in the Scaling for Upwelling Over Submarine Canyons. Howatt, Tara 2012-12-31
pdf
Page Metadata
Item Metadata
Title | A Numerical Investigation of the Shelf Slope in the Scaling for Upwelling Over Submarine Canyons. |
Creator |
Howatt, Tara |
Date Issued | 2012 |
Description | Submarine canyons are known regions of enhanced upwelling. Previous studies have looked into the dynamics of the upwelling and recently a scaling system by Allen and Hickey [2010] has been developed. A closer look into this scaling has shown that the continental shelf slope surrounding the canyon plays a role and needs to be included in the scaling. The slope plays a role in the depth of upwelling and the upwelling ux because the slope will create a back pressure restricting the ow from moving further up the shelf. This thesis will discuss the determination of a non-dimensional number to characterize the shelf slope in these upwelling quantities. A numerical model will be used to test the new scaling and to assist in nding a relation between the shelf slope and the depth of upwelling and the upwelling ux. A non-dimensional number to take into account the e ect of the shelf slope has been found; however, further research is required for a physical explanation of the form of this number. Looking at the results from the model, relationships between the slope, strati cation, Coriolis parameter and upwelling depth and ux can also be made, where a atter slope, low strati cation, and a low Coriolis parameter will have the greatest upwelling response. |
Type |
Text |
Language | eng |
Series | University of British Columbia. EOSC 449 |
Date Available | 2012-04-19 |
Provider | Vancouver : University of British Columbia Library |
DOI | 10.14288/1.0053610 |
URI | http://hdl.handle.net/2429/42084 |
Affiliation |
Science, Faculty of Earth and Ocean Sciences, Department of |
Citation | Howatt, Tara. 2012. A Numerical Investigation of the Shelf Slope in the Scaling for Upwelling Over Submarine Canyons. Undergraduate Honours Thesis. Department of Earth and Ocean Sciences. University of British Columbia. |
Peer Review Status | Unreviewed |
Scholarly Level | Undergraduate |
Copyright Holder | Howatt, Tara |
Aggregated Source Repository | DSpace |
Download
- Media
- 52966-Howatt_Tara_UBC_2012_EOSC_Honours_Thesis.pdf [ 885.56kB ]
- Metadata
- JSON: 52966-1.0053610.json
- JSON-LD: 52966-1.0053610-ld.json
- RDF/XML (Pretty): 52966-1.0053610-rdf.xml
- RDF/JSON: 52966-1.0053610-rdf.json
- Turtle: 52966-1.0053610-turtle.txt
- N-Triples: 52966-1.0053610-rdf-ntriples.txt
- Original Record: 52966-1.0053610-source.json
- Full Text
- 52966-1.0053610-fulltext.txt
- Citation
- 52966-1.0053610.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.52966.1-0053610/manifest