Knight Inlet Research Articles

Synallactesmcdanieli sp. nov., a new species of sea cucumber from British Columbia, Canada and the Gulf of Alaska, USA (Holothuroidea, Synallactida).

The family Synallactidae comprises mostly deep-sea forms and is the least-studied large taxon amongst deep-sea cucumbers. They are part of the abyssal megafauna and play an important role in modifying the sediment landscape and structuring the communities that live within it. The family embraces the genus Synallactes, which contains approximately twenty-five species from the Pacific, Atlantic (six species), Indian (seven species) and Antarctic Oceans (one species). Synallactesmcdanieli sp. nov. is described from the Northeast Pacific, Knight Inlet, British Columbia, Canada to Kodiak Island, Gulf of Alaska, USA, at depths from 21 to 438 m. This new species is unique amongst the species of the genus Synallactes because of the number and arrangement of dorsal papillae, number of Polian vesicles, together with the entire ossicle arrangement. In addition, this species has the shallowest bathymetric distribution ever recorded for this genus.

Paired environmental DNA and dive surveys provide distinct but complementary snapshots of marine biodiversity in a temperate fjord

AbstractMarine biodiversity is a key indicator of ecosystem health and can be assessed using a variety of methods, including environmental DNA (eDNA) sampling. However, the ecology of eDNA in physically dynamic nearshore environments remains uncertain, particularly with regards to how eDNA stratifies with depth. Here, we paired eDNA sampling with dive surveys at six sites in Knight Inlet, British Columbia, Canada. eDNA samples were collected from the surface, midwater column and bottom (8–25 m depth) at each site, while dive surveys focused on the bottom (benthic) habitat. Amplicon sequencing using the mitochondrial 12S rRNA gene (targeting fish) and the COI gene (targeting marine invertebrates and algae) resolved significant differences in community composition in surface waters compared with midwater and bottom. Differences by depth were greater than differences across sites, with surface waters dominated by salmon (Oncorhynchus spp.) and rotifer DNA, and midwater and bottom samples largely dominated by Pacific herring, copepods, and mussels. eDNA samples collected at the surface, therefore, may not accurately capture benthic communities, particularly in systems with high levels of freshwater input such as coastal temperate fjords. Over small spatial scales, particularly in systems with strong stratification, adding samples from different depths may be more effective at maximizing inferred diversity rather than sampling more sites. In general, there was low overlap in species detection between dive and eDNA surveys (less than 10% for each taxonomic group – fish, invertebrates, and algae). However, we observed clear strengths for each method – dive surveys provided better taxonomic resolution, while eDNA resolved greater total diversity. These results suggest that the two survey methods can be used in tandem to provide distinct and complementary snapshots of marine biodiversity in the nearshore environment.

Deep Waters in British Columbia Mainland Fjords Show Rapid Warming and Deoxygenation From 1951 to 2020

AbstractMany complex fjord systems cross British Columbia's coastline. A 70 year (1951–2020) time series analysis of temperature, salinity, and oxygen in four such fjords between ∼54 and 50oN (Douglas Channel, Rivers Inlet, Knight Inlet and Bute Inlet) shows that changes were greatest in deep waters between the sill and the bottom. In Rivers, Knight and Bute Inlet, the deep water temperature increased by 1.2–1.3°C over 70 years, up to two times the global average for open ocean waters at corresponding depths, while salinity increased by 0.1–0.2, and oxygen decreased by 0.4–0.7 mLL−1. The most northern inlet, Douglas Channel, showed a temperature increase of 0.8°C from 1951 to 2016, while trends in oxygen and salinity were not statistically significant. An analysis of Apparent Oxygen Utilization suggests that the deep waters in Douglas Channel are more readily exchanged with the outer coast than the three other fjords.

U-Th-Pb Geochronology and Lu-Hf Isotope Geochemistry of Detrital Zircons in Metasedimentary Rocks of the Southern Coast Mountains Batholith

Abstract Most metasedimentary rocks in the southern Coast Mountains batholith are of uncertain tectonic affinity because they occur in discontinuous pendants surrounded by large intrusive bodies, and many protolith features are obscured by regional deformation and metamorphism. This study uses U-Th-Pb ages and Lu-Hf isotope signatures of detrital zircons in metasedimentary rocks in Bute, Loughborough, and Knight Inlets in an effort to test possible correlations with the adjacent Wrangellia, Alexander, Taku, Yukon-Tanana, and Stikine terranes. Detrital zircons from metasedimentary samples yield ages that belong to age groups of 590-528 Ma (peak age of 560 Ma), 485-432 Ma (peak age of 452 Ma), 356-286 Ma (peak age of 307 Ma), and 228-185 Ma (peak ages of 215 and 198 Ma). A small number of ~1.1-1.9 Ga grains are also present. εHft values of the 590-185 Ma grains yield a progression from intermediate (0 to +5) values to more juvenile (mostly +4 to +15) values from Neoproterozoic through early Mesozoic time. The Comparison of these results with similar data sets from adjacent terranes demonstrates that primary connections with the Yukon-Tanana and Taku terranes are unlikely but are consistent with primary connections with the Wrangellia, Stikine, and/or Alexander terranes. Unfortunately, the available constraints are not sufficient to eliminate any of these options or the possibility that the pendants are a unique tectonic fragment. Zircons from the metasedimentary samples also yield U-Th-Pb ages of 165-128 Ma (peak age of 152 Ma) and 114-88 Ma (peak age of 102 Ma). εHft values of these zircon domains are mostly juvenile (+7 to +13). Comparison of U concentrations, U/Th values, and CL textures of zircons from the metasedimentary samples, leucocratic sills that intrude the pendants, and surrounding plutonic bodies suggests that most of the young grains, as well as widespread younger rims on older grains, grew during metamorphism associated with emplacement of the adjacent plutonic bodies. Some young grains were derived from thin felsic sills or veins that were unintentionally included in the sampled material.

Internal solitary wave generation by tidal flow over topography

Oceanic internal solitary waves are typically generated by barotropic tidal flow over localised topography. Wave generation can be characterised by the Froude number $F=U/c_{0}$, where $U$ is the tidal flow amplitude and $c_{0}$ is the intrinsic linear long wave phase speed, that is the speed in the absence of the tidal current. For steady tidal flow in the resonant regime, $\unicode[STIX]{x1D6E5}_{m}<F-1<\unicode[STIX]{x1D6E5}_{M}$, a theory based on the forced Korteweg–de Vries equation shows that upstream and downstream propagating undular bores are produced. The bandwidth limits $\unicode[STIX]{x1D6E5}_{m,M}$ depend on the height (or depth) of the topographic forcing term, which can be either positive or negative depending on whether the topography is equivalent to a hole or a sill. Here the wave generation process is studied numerically using a forced Korteweg–de Vries equation model with time-dependent Froude number, $F(t)$, representative of realistic tidal flow. The response depends on $\unicode[STIX]{x1D6E5}_{max}=F_{max}-1$, where $F_{max}$ is the maximum of $F(t)$ over half of a tidal cycle. When $\unicode[STIX]{x1D6E5}_{max}<\unicode[STIX]{x1D6E5}_{m}$ the flow is always subcritical and internal solitary waves appear after release of the downstream disturbance. When $\unicode[STIX]{x1D6E5}_{m}<\unicode[STIX]{x1D6E5}_{max}<\unicode[STIX]{x1D6E5}_{M}$ the flow reaches criticality at its peak, producing upstream and downstream undular bores that are released as the tide slackens. When $\unicode[STIX]{x1D6E5}_{max}>\unicode[STIX]{x1D6E5}_{M}$ the tidal flow goes through the resonant regime twice, producing undular bores with each passage. The numerical simulations are for both symmetrical topography, and for asymmetric topography representative of Stellwagen Bank and Knight Inlet.

Passing on the excitement of experimental oceanography

Having graduated from the University of Victoria in 1995 with David Farmer as my advisor, I learned how exciting it is to work at the cutting-edge of experimental oceanography. During my Ph.D., I studied bubble mediated air-sea gas exchange which involved chasing winter storms in the NE Pacific by ship. That work inspired our study of the role of tropical cyclones in the global carbon cycle using air deployed gas sensing floats equipped with ambient noise recorders to evaluate wave breaking effects. I was also inspired when David enthusiastically showed me on the echosounder internal waves passing under the boat in Knight Inlet. I will present our recent AUV measurements made in prominent estuarine features (e.g., salt wedge and internal hydraulic jump) and discuss the impact of these features on suspended sediment distributions and underwater acoustic communications.

Generation of internal solitary waves over a large sill: From Knight Inlet to Luzon Strait

A fully nonlinear, nonhydrostatic numerical model is utilized to investigate the generation of Internal Solitary Waves (ISWs) upstream of the Knight Inlet sill. While an upstream hydraulic jump initiates the ISW generation and both hydraulic jump and upstream influence contribute to the generation, it is found that upstream influence is dominant and the hydraulic jump is not necessary for the ultimate generation of ISWs. Decreasing the tidal forcing or upstream sill width may render the flow subcritical (i.e., the hydraulic jump disappears) and ISWs can be generated by nonlinear steepening of long wave disturbances induced by upstream influence. Increasing the tidal forcing or upstream sill width may generate a hydraulic jump blocking strong upstream propagating disturbances. The jump subsequently becomes a turbulent bore and later disperses into a train of ISWs as the tide relaxes. Further increase in the tidal forcing may sweep the turbulent bore downstream and a train of ISWs is emitted upstream toward the end of waning tide. By reducing the stratification strength by 1 order of magnitude, the near-sill flow is in the transcritical regime and ISWs are resonantly generated over the lee side slope. Connections to the internal tide release mechanism at Luzon Strait and to the unsteady lee wave model are also discussed. The present work provides some more insights into the ISW generation process at Knight Inlet and the connection between the generation mechanism at Knight Inlet and that at Luzon Strait is identified.

Structure, stability, and tsunami hazard associated with a rock slope in Knight Inlet, British Columbia

Abstract. Rockfalls and rockslides during the past 12 000 years have deposited bouldery debris cones on the seafloor beneath massive rock slopes throughout the inner part of Knight Inlet. The 885 m high rock slope, located across from a former First Nations village destroyed in the late 1500s by a slide-induced wave, exposes the contact between a Late Cretaceous dioritic pluton and metamorphic rocks of the Upper Triassic Karmutsen Formation. The pluton margin is strongly foliated parallel to primary and secondary fabrics in the metamorphic rocks, resulting in highly persistent brittle structures. Other important structures include a set of sheeting joints and highly persistent mafic dykes and faults. Stability analysis indicates that planar and wedge rock slope failures up to about 500 000 m3 in volume could occur. We suspect that failures of this size in this setting would have the potential to generate locally hazardous waves. As several similar rock slopes fronted by large submarine debris cones exist in the inner part of Knight Inlet, it is clear that tsunami hazards should be considered in coastal infrastructure development and land-use planning in this area.

Space–time cluster analysis of sea lice infestation (Caligus clemensi and Lepeophtheirus salmonis) on wild juvenile Pacific salmon in the Broughton Archipelago of Canada

Sea lice infestation levels on wild chum and pink salmon in the Broughton Archipelago region are known to vary spatially and temporally; however, the locations of areas associated with a high infestation level had not been investigated yet. In the present study, the multivariate spatial scan statistic based on a Poisson model was used to assess spatial clustering of elevated sea lice (Caligus clemensi and Lepeophtheirus salmonis) infestation levels on wild chum and pink salmon sampled between March and July of 2004 to 2012 in the Broughton Archipelago and Knight Inlet regions of British Columbia, Canada. Three covariates, seine type (beach and purse seining), fish size, and year effect, were used to provide adjustment within the analyses. The analyses were carried out across the five months/datasets and between two fish species to assess the consistency of the identified clusters. Sea lice stages were explored separately for the early life stages (non-motile) and the late life stages of sea lice (motile). Spatial patterns in fish migration were also explored using monthly plots showing the average number of each fish species captured per sampling site. The results revealed three clusters for non-motile C. clemensi, two clusters for non-motile L. salmonis, and one cluster for the motile stage in each of the sea lice species. In general, the location and timing of clusters detected for both fish species were similar. Early in the season, the clusters of elevated sea lice infestation levels on wild fish are detected in areas closer to the rivers, with decreasing relative risks as the season progresses. Clusters were detected further from the estuaries later in the season, accompanied by increasing relative risks. In addition, the plots for fish migration exhibit similar patterns for both fish species in that, as expected, the juveniles move from the rivers toward the open ocean as the season progresses The identification of space–time clustering of infestation on wild fish from this study can help in targeting investigations of factors associated with these infestations and thereby support the development of more effective sea lice control measures.

Paleogeography of the Insular and Intermontane terranes reconsidered: Evidence from the southern Coast Mountains Batholith, British Columbia

New geologic and paleomagnetic data from Knight Inlet in the southwestern Coast Mountains Batholith, British Columbia, support signifi cant revision to the paleogeography of the Insular and Intermontane terranes. Recompilation of radiometric ages confi rms that after 100 Ma, a magmatic arc migrated northeastward across the Coast Mountains Batholith at ~2 km/m.y. Magmatic age patterns suggest that plutons older than 100 Ma intruded the Intermontane terrane, not the expected Insular terrane. The distribution of brittle faults along Knight Inlet defi nes a structurally intact central domain, ~45 km wide, fl anked to the SW and NE by faulted domains, with no evidence of the widespread Tertiary extension affecting the batholith farther north. Al-in-hornblende geobarometry yields emplacement depths of ~2.5‐4 kbar and does not reveal systematic postemplacement tilting. Plutons in the central structural domain yield a consistently oriented paleomagnetic remanence presumably acquired as the Late Cretaceous arc cooled from ca. 110 to 85 Ma. In the absence of recognizable tilting, this result indicates ~1700 km of northward translation since ca. 85 Ma, which is signifi cantly less than predicted for the Insular terrane in the “Baja British Columbia” model but similar to results from the Intermontane terrane. The pluton ages and the paleomagnetic results suggest that the Intermontane terrane, not the Insular terrane, underlies the southwestern fl ank of the Coast Mountains Batholith. This conclusion is compatible with a paleogeographic model in which the Vancouver Island fragment of Wrangellia was juxtaposed against the Intermontane terrane prior to ca. 120‐100 Ma and emplaced in southern British Columbia after ca. 75 Ma.

Sea lice infestations on juvenile chum and pink salmon in the Broughton Archipelago, Canada, from 2003 to 2012

Juvenile pink salmon Oncorhynchus gorbuscha and chum salmon O. keta were sampled by beach or purse seine to assess levels of sea lice infestation in the Knight Inlet and Broughton Archipelago regions of coastal British Columbia, Canada, during the months of March to July from 2003 to 2012. Beach seine data were analyzed for sea lice infestation that was described in terms of prevalence, abundance, intensity, and intensity per unit length. The median annual prevalence for chum was 30%, ranging from 14% (in 2008 and 2009) to 73% (in 2004), while for pink salmon, the median was 27% and ranged from 10% (in 2011) to 68% (in 2004). Annual abundance varied from 0.2 to 5 sea lice per fish with a median of 0.47 for chum and from 0.1 to 3 lice (median 0.42) for pink salmon. Annual infestation followed broadly similar trends for both chum and pink salmon. However, the abundance and intensity of Lepeophtheirus salmonis and Caligus clemensi, the 2 main sea lice species of interest, were significantly greater on chum than on pink salmon in around half of the years studied. Logistic regression with random effect was used to model prevalence of sea lice infestation for the combined beach and purse seine data. The model suggested inter-annual variation as well as a spatial clustering effect on the prevalence of sea lice infestation in both chum and pink salmon. Fish length had an effect on prevalence, although the nature of this effect differed according to host species.

Submarine channel evolution: active channels in fjords, British Columbia, Canada

Seafloor channel systems are well developed in several British Columbia (Canada) fjords, including Knight, Bute and Toba inlets. The channels are active conduits for turbidity currents during episodic failures of delta fronts at the fjord heads. Recently acquired swath multibeam bathymetric data enabled a complete and detailed assessment of the form of these channel systems. Repeated multibeam surveys in Bute and Toba inlets (Mar. 2008 and Nov. 2010) indicate that the deeply incised channels may undergo accretion or erosion over a period as short as 2.5 years. Locally significant (>15 m) bathymetric changes were documented in the channel systems, with bathymetric differences exceeding 5 m occurring along approx. 20–25% of the resurveyed channel lengths. Slope failure involving the upper delta front observed in Toba Inlet is a mechanism for the rapid initiation of new channel segments. Sediment accretion on the inner side of some channel bends indicates that point bars may form as lateral channel migration proceeds. Channel sidewall slumps are found at some outer channel bends and mounded deposits may occur at these sites. Lateral migration is accompanied by outer bend erosion, overall deepening of the channel and the development of terraces.

Response of Euphausia pacifica to small-scale shear in turbulent flow over a sill in a fjord

Zooplankton in the ocean respond to visual and hydro-mechanical cues such as small-scale shear in turbulent flow. In addition, they form strong aggregations where currents intersect sloping bottoms. Strong and predictable tidal currents over a sill in Knight Inlet, Canada, make it an ideal location to investigate biological behaviour in turbulent cross-isobath flow. We examine acoustic data (38, 120 and 200 kHz) collected there during the daylight hours, when the dominant zooplankters, Euphausia pacifica have descended into low light levels at ∼90 m. As expected, these data reveal strong aggregations at the sill. However, they occur consistently 10–20 m below the preferred light depth of the animals. We have constructed a simple model of the flow to investigate this phenomenon. Tracks of individual animals are traced in the flow and a variety of zooplankton behaviours tested. Our results indicate that the euphausiids must actively swim downward when they encounter the bottom boundary layer (bbl) to reproduce the observed downward shift in aggregation patterns. We suggest that this behaviour is cued by the small-scale shear in the bbl. Furthermore, this behaviour is likely to enhance aggregations found in strong flows at sills and on continental shelves.

Turbulence closure: turbulence, waves and the wave-turbulence transition – Part 1: Vanishing mean shear

Abstract. This paper extends a turbulence closure-like model for stably stratified flows into a new dynamic domain in which turbulence is generated by internal gravity waves rather than mean shear. The model turbulent kinetic energy (TKE, K) balance, its first equation, incorporates a term for the energy transfer from internal waves to turbulence. This energy source is in addition to the traditional shear production. The second variable of the new two-equation model is the turbulent enstrophy (Ω). Compared to the traditional shear-only case, the Ω-equation is modified to account for the effect of the waves on the turbulence time and space scales. This modification is based on the assumption of a non-zero constant flux Richardson number in the limit of vanishing mean shear when turbulence is produced exclusively by internal waves. This paper is part 1 of a continuing theoretical development. It accounts for mean shear- and internal wave-driven mixing only in the two limits of mean shear and no waves and waves but no mean shear, respectively. The new model reproduces the wave-turbulence transition analyzed by D'Asaro and Lien (2000b). At small energy density E of the internal wave field, the turbulent dissipation rate (ε) scales like ε~E2. This is what is observed in the deep sea. With increasing E, after the wave-turbulence transition has been passed, the scaling changes to ε~E1. This is observed, for example, in the highly energetic tidal flow near a sill in Knight Inlet. The new model further exhibits a turbulent length scale proportional to the Ozmidov scale, as observed in the ocean, and predicts the ratio between the turbulent Thorpe and Ozmidov length scales well within the range observed in the ocean.

Numerical modelling of stratified tidal flow over a fjord sill

Field observations of tidally driven stratified flow in the sill area of Knight Inlet (British Columbia) revealed a very complicated structure, which includes solitary waves, upstream bifurcation, hydraulic jump and mixing processes. Recent observations suggest that the flow instabilities on the plunging pycnocline at the lee side of the sill may contribute to solitary wave generation through a subharmonic interaction. The present study reports on a series of numerical experiments of stratified tidal flow in Knight Inlet performed with the help of a fine resolution fully non-linear non-hydrostatic numerical model. The model reproduces all important stages of the baroclinic tidal dynamics observed in Knight Inlet. Results demonstrate that solitary waves are generated apart from the area of hydrodynamic instability. Accelerating tidal flux forms a baroclinic hydraulic jump just above the top of the sill, whereas the bifurcations and zones of shear instabilities are formed downstream of the sill. The first baroclinic mode having the largest velocity escapes from the generation area and propagates upstream, disintegrating further into a packet of solitary waves reviling the classical “non-subharmonic” mechanism of generation. The remaining part of the disturbance (slow baroclinic modes) is arrested by tidal flow and carried away to the lee side of the obstacle, where shear instability, billows and mixing processes are developed. Some sensitivity runs were performed for different value of tidal velocity.

Destruction of the first nations village of Kwalate by a rock avalanche‐generated tsunami

The First Nations (Da ‘naxda ‘xw) village of Kwalate, Knight Inlet, British Columbia was located along the shore of a funnel‐shaped bay. Archaeological investigations show that this was a major village that stretched 90 m along the shoreline and was home to possibly 100 or more inhabitants. Oral stories indicate that the village was completely swept away by a tsunami that formed when an 840‐m high rock avalanche descended into the water on the opposite side of the fjord. Shipboard geological mapping, combined with empirical tsunami modelling, indicate that the tsunami was likely 2 to 6 m high prior to run‐up into the village. Radiocarbon dates reveal that the village was occupied from the late 1300s CE until the late 1500s CE when it was destroyed by the tsunami.

Mapping fish schools with a high-resolution multibeam sonar

In late 2002 an evaluation of a high-resolution multibeam sonar (RESON SeaBat 8125) for mapping mid-water fish schools was performed in Knight Inlet, BC. The ship-mounted multibeam sonar operated at 455 kHz and provided 240 one-half-degree beams in the across-track plane with one-degree along-track beam width. Simultaneous measurements were made with a conventional 200-kHz echo-sounder. Acoustic calibrations of this multibeam sonar will be discussed, showing some range- and angle-dependent variations in sensitivity. The multibeam sonar was able to provide well-resolved mid-water images of fish schools (likely juvenile herring), sampling many fish schools missed by the echo-sounder and revealing cross-track fish school shape. Comparison of near-vertical volumetric scattering strengths between the multi-beam sonar and echo-sounder show reasonable agreement. Discussions on the combined use of conventional and multibeam sonars for fish surveys will be given.

Amplification and Suppression of Internal Waves by Tides over Variable Bottom Topography

Abstract The energy exchange between internal waves and barotropic currents over inclined bottom topography is studied theoretically and in the framework of the numerical model. The energy balance equation derived for a continuously stratified fluid predicts that energy can either be transferred toward or away from the internal wave depending on the direction of propagation of both the wave and current. Four scenarios of wave–flow interaction over the inclined bottom were identified. An internal wave extracts energy from the background tidal flow during its propagation upslope–upstream or downslope–downstream and its amplitude grows. The wave loses energy propagating downslope–upstream or upslope–downstream and reduces in amplitude. This mechanism of suppression or amplification of internal waves by a current over an inclined bottom is verified numerically. When applied to the area of the Knight Inlet sill, a high-resolution fully nonlinear, nonhydrostatic model reproduces the packets of internal waves generated by supercritical tidal flow over the sill. Careful inspection of the wave fields revealed the presence of an irregular wave structure within wave packets—namely, internal waves are not arranged by amplitude. This phenomenon, obtained numerically and observed in situ, is treated in terms of the mechanism of wave–flow interaction: the energy exchange between the tidal current and generated internal waves over the inclined bottom topography is the reason for the absence of traditional rank-ordered waves in the packet.

Estuarine and Tidal Currents in the Broughton Archipelago

Current meter observations, Conductivity‐Temperature‐Depth (CTD) profiles, and river discharges are combined with two numerical models to understand better the estuarine and tidal circulation in the Broughton Archipelago, a complex region of islands, channels, and fiords that has become a primary location for salmon farms in British Columbia. Though tidal currents are strong in many sub‐regions, the primary transport mechanisms in the archipelago are the estuarine flows resulting from river and glacial runoff, and the near‐surface currents that arise from strong winds. The harmonic finite element model, TIDE3D, is shown to reproduce the barotropic tidal currents with reasonable accuracy, but unlike other regions of the British Columbia coast where an extensive archive of historical observations has permitted a diagnostic calculation of average seasonal flows, sparse and noisy CTD observations did not allow a similar computation here. In order to simulate these background flows, the prognostic finite volume model, ELCIRC, was initialized with a smoothed version of these same historical temperature and salinity fields and forced with tides and river discharge. Though the near‐surface flows were reproduced with acceptable accuracy, the estuarine and tidal flows at depth were found to be much too weak as a result of numerical dissipation arising from the Eulerian‐Lagrangian time stepping. Nevertheless, ELCIRC did confirm current observations suggesting that the bottom estuarine flow in Knight Inlet actually comes from Queen Charlotte Strait via the “back‐door” of Fife Sound and Tribune Passage. Consistent with anecdotal evidence, ELCIRC also showed that the surface estuarine flow coming down Knight Inlet bifurcates with part going down Tribune Channel and Fife Sound and part continuing down Knight Inlet. The relevance of these background flows for aquaculture issues, such as oxygen renewal and the transport of sea lice and viruses, is discussed.

Energy Fluxes due to the Surface and Internal Tides in Knight Inlet, British Columbia

Abstract A laterally integrated (two dimensional) nonlinear numerical model is used to examine the flux of M2 tidal energy in Knight Inlet. The simulated flux of tidal energy into the inlet is somewhat smaller than that estimated using the change in phase of the M2 tidal height along the inlet, a method that does not account for the effect of the internal tide on the surface elevation. The simulated energy flux into the inlet is close to the energy flux of the internal tide away from the sill determined from observations using an acoustic Doppler current profiler (ADCP). The net flux due to the internal tide is significantly less than (<1/2 of) the rate at which energy is removed from the surface tide. Earlier linear models of the internal tide produced energy fluxes that agreed with those estimated from the phase change of the tidal height but were larger than the fluxes that could be found in the observations. The reason for this discrepancy is not that these simple models neglected nonlinear effects, but rather that they did not take reflections of the internal tide into account. Also, the simulated flux of energy into the inlet less the net flux of internal tidal energy away from the sill is about equal to the simulated dissipation within 2 km on either side of the sill. The simulated net flux of internal tidal energy away from the sill is in agreement with the flux determined from the ADCP observations on the downinlet side of the sill, but not on the upinlet side of the sill. A possible explanation is that only the first internal mode (which is surface intensified) was important on the downinlet side but the first three internal modes were important on the upinlet side. The flux calculation using the ADCP observations took variations in the inlet width into account but did not take depth variations into account; thus, the reflection coefficients of the second and third modes may have been underestimated.