The variability of the sediment plume and ocean circulation features of the Nass River Estuary, British Columbia

The Nass River discharges into Nass Bay and Iceberg Bay, which are adjoining tidal inlets located within the northern inland waters of British Columbia, Canada. After the Skeena River, the Nass River is the second longest river within northern British Columbia, which discharges directly into Canadian waters of the Pacific Ocean. It is also supports one of the most productive salmon fisheries in northern British Columbia. The Nass River discharges into the eastern end of Nass Bay. Nass Bay, in turn feeds into Portland Canal and the fresh surface waters then flows westward to the Pacific Ocean via Dixon Entrance. The tides in Northern British Columbia are very large with a tidal height range of just over 7 m. Nass Bay is a shallow inlet of less than 10 km in length with typical water depths of than 10 m or less. The existing knowledge of oceanographic processes in Nass and Iceberg Bays was rudimentary until three years ago, when the first modern oceanographic measurements were obtained. In this study, the seasonal and tidal variability of the lateral extent of the Nass River surface plume is mapped from analyses of Landsat satellite data spanning the period from 2008 to 2015. A high resolution coupled three dimensional (3D) hydrodynamic model was developed and implemented, within the widely used and accepted Delft3D modeling framework, which was forced and validated using recent 2013-2016 in-situ oceanographic measurements. The combined satellite and numerical modeling methods are used to study the physical oceanographic and sediment transport regime of Nass and Iceberg Bays and the adjoining waters of Portland Inlet and Observatory Inlet. The ocean circulation of Nass and Iceberg Bays was found to be dominated by tidal currents, and by the highly seasonal and variable Nass River freshwater discharges. Complex lateral spatial patterns in the tidal currents occur due to the opening of the southwestern side of Nass Bay onto the deeper adjoining waters of Iceberg Bay. Surface winds are limited to a secondary role in the circulation variability. The sediment dynamics of the Nass Bay system features a very prominent surface sediment plume present from the time of freshet in mid-spring through to large rainfall runoff events in the fall. The time-varying turbidity distribution and transport paths of the Nass River sediment discharges in the study area were characterized using the model results combined with an analysis of several high-resolution multi-year Landsat satellite data sets

Circulation in the vicinity of Mackenzie Canyon from a year-long mooring array

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Lin, P., Pickart, R. S., Fissel, D., Ross, E., Kasper, J., Bahr, F., Torres, D. J., O'Brien, J., Borg, K., Melling, H., & Wiese, F. K. Circulation in the vicinity of Mackenzie Canyon from a year-long mooring array. Progress in Oceanography, 187, (2020): 102396, doi:10.1016/j.pocean.2020.102396.Data from a five-mooring array extending from the inner shelf to the continental slope in the vicinity of Mackenzie Canyon, Beaufort Sea are analyzed to elucidate the components of the boundary current system and their variability. The array, part of the Marine Arctic Ecosystem Study (MARES), was deployed from October 2016 to September 2017. Four distinct currents were identified: an eastward-directed flow adjacent to the coast; a westward-flowing, surface-intensified current centered on the outer-shelf; a bottom-intensified shelfbreak jet flowing to the east; and a recirculation at the base of the continental slope within the canyon. The shelf current transports −0.120.03 Sv in the mean and is primarily wind-driven. The response is modulated by the presence of ice, with little-to-no signal during periods of nearly-immobile ice cover and maximum response when there is partial ice cover. The shelfbreak jet transports 0.030.02 Sv in the mean, compared to 0.080.02 Sv measured upstream in the Alaskan Beaufort Sea over the same time period. The loss of transport is consistent with a previous energetics analysis and the lack of Pacific-origin summer water downstream. The recirculation in the canyon appears to be the result of local dynamics whereby a portion of the westward-flowing southern limb of the Beaufort Gyre is diverted up the canyon across isobaths. This interpretation is supported by the fact that the low-frequency variability of the recirculation is correlated with the wind-stress curl in the Canada Basin, which drives the Beaufort gyre.The authors are indebted to Fisheries and Oceans Canada for building the logistics for MARES into the at-sea missions of the Integrated Beaufort Observatory. We are grateful to the captain and crew of the CCGS Sir Wilfred Laurier for ably deploying and recovering the MARES array. Marshall Swartz assisted with the cruise preparation logistics. We thank the two anonymous reviewers for their input which helped improve the paper. This project was funded by the US Bureau of Ocean Energy Management (BOEM), on behalf of the National Ocean Partnership Program. The Canadian contribution was supported by the Environmental Studies Research Fund (ESRF Project 2014-02N). MARES publication 003

Comparison of prestellar core elongations and large-scale molecular cloud structures in the Lupus 1 region

Turbulence and magnetic fields are expected to be important for regulating molecular cloud formation and evolution. However, their effects on sub-parsec to 100 parsec scales, leading to the formation of starless cores, are not well understood. We investigate the prestellar core structure morphologies obtained from analysis of the Herschel-SPIRE 350 mum maps of the Lupus I cloud. This distribution is first compared on a statistical basis to the large-scale shape of the main filament. We find the distribution of the elongation position angle of the cores to be consistent with a random distribution, which means no specific orientation of the morphology of the cores is observed with respect to the mean orientation of the large-scale filament in Lupus I, nor relative to a large-scale bent filament model. This distribution is also compared to the mean orientation of the large-scale magnetic fields probed at 350 mum with the Balloon-borne Large Aperture Telescope for Polarimetry during its 2010 campaign. Here again we do not find any correlation between the core morphology distribution and the average orientation of the magnetic fields on parsec scales. Our main conclusion is that the local filament dynamics---including secondary filaments that often run orthogonally to the primary filament---and possibly small-scale variations in the local magnetic field direction, could be the dominant factors for explaining the final orientation of each core

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

We present and analyze observations of polarized dust emission at 850 μm toward the central 1 7 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields in Star-forming Region Observations survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis–Chandrasekhar–Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from Herschel data and the C18O (J = 3 - 2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 \ub1 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 \ub1 0.02. Additionally, the mean Alfv\ue9n Mach number is 0.35 \ub1 0.01. This suggests that, in Mon R2, the magnetic fields provide resistance against large-scale gravitational collapse, and the magnetic pressure exceeds the turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically subcritical

HAWC+/SOFIA Multiwavelength Polarimetric Observations of OMC-1

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy. We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 μm at angular resolutions of 5'', 8'', 14'', and 19'' for the four bands, respectively. The photometric maps enable the computation of improved spectral energy distributions for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the "hourglass" configuration seen in previous studies, indicating magnetically regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin–Neugebauer/Kleinman–Low outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions

A Numerical Study of Long-Return Period Near-Bottom Ocean Currents in Lower Cook Inlet, Alaska

Lower Cook Inlet (LCI) is an important waterway with very large tides and high marine productivity. Oceanographic forcing in LCI is complex due to a combination of tides, seasonal winds, and large freshwater discharges, as well as inflow from the Alaska Coastal Current. From an analysis of historical current meter data sets, deeper ocean currents of LCI were found to have large differences resulting from the dominance of large tides in the northeast portion of LCI while subtidal contributions to the deeper currents are more important relative to the reduced tidal currents in central and western parts of LCI. To compute the largest values of the near-bottom currents of LCI, a 3D hydrodynamical model was developed over a large model domain extending over the full 300 km length of Cook Inlet as well as a large portion of the adjoining Alaska continental shelf region. At the open model boundaries, nine major tidal height constituents were specified based on National Oceanic and Atmospheric Administration (NOAA) tidal gauge data. The model was forced by the spatially varying winds and freshwater discharges for the six gauged rivers in Cook Inlet. The model was verified using available current meter data in the study area. Model runs were carried out for 21 case studies to derive the near-bottom currents for return periods of 1, 10, and 100 years. Within LCI, the extremal values for near bottom currents arise from quite different forcing regimes. Tidal currents are completely dominant in the northeast portion of LCI while for central and western portions, remote wind forcing over the Alaskan continental shelf current, which generates the Alaska Coastal Current, becomes more important

High Resolution 3-D Finite-Volume Coastal Ocean Modeling in Lower Campbell River and Discovery Passage, British Columbia, Canada

The 3-D unstructured-grid, Finite-Volume Coastal Ocean Model (FVCOM) was used to simulate the flows in Discovery Passage including the adjoining Lower Campbell River, British Columbia, Canada. Challenges in the studies include the strong tidal currents (e.g., up to 7.8 m/s in Seymour Narrows) and tailrace discharges, small-scale topographic features and steep bottom slopes, and stratification affected by the Campbell River freshwater discharges. Two applications of high resolution 3-D FVCOM modeling were conducted. One is for the Lower Campbell River extending upstream as far as the John Hart Hydroelectric dam. The horizontal resolution varies from 0.27 m to 32 m in the unstructured triangular mesh to resolve the tailrace flow. The bottom elevation decreases ~14 m within the distance of ~1.4 km along the river. This pioneering FVCOM river modeling demonstrated a very good performance in simulating the river flow structures. The second application is to compute ocean currents immediately above the seabed along the present underwater electrical cable crossing routes across Discovery Passage. Higher resolution was used near the bottom with inter-layer spacing ranging from 0.125 to 0.0005 of total water depth. The model behaves very well in simulating the strong tidal currents in the area at high resolution in both the horizontal and vertical. One year maximum near bottom tidal current along the routes was then analyzed using the model results

A Numerical Study of Long-Return Period Near-Bottom Ocean Currents in Lower Cook Inlet, Alaska

Lower Cook Inlet (LCI) is an important waterway with very large tides and high marine productivity. Oceanographic forcing in LCI is complex due to a combination of tides, seasonal winds, and large freshwater discharges, as well as inflow from the Alaska Coastal Current. From an analysis of historical current meter data sets, deeper ocean currents of LCI were found to have large differences resulting from the dominance of large tides in the northeast portion of LCI while subtidal contributions to the deeper currents are more important relative to the reduced tidal currents in central and western parts of LCI. To compute the largest values of the near-bottom currents of LCI, a 3D hydrodynamical model was developed over a large model domain extending over the full 300 km length of Cook Inlet as well as a large portion of the adjoining Alaska continental shelf region. At the open model boundaries, nine major tidal height constituents were specified based on National Oceanic and Atmospheric Administration (NOAA) tidal gauge data. The model was forced by the spatially varying winds and freshwater discharges for the six gauged rivers in Cook Inlet. The model was verified using available current meter data in the study area. Model runs were carried out for 21 case studies to derive the near-bottom currents for return periods of 1, 10, and 100 years. Within LCI, the extremal values for near bottom currents arise from quite different forcing regimes. Tidal currents are completely dominant in the northeast portion of LCI while for central and western portions, remote wind forcing over the Alaskan continental shelf current, which generates the Alaska Coastal Current, becomes more important

Modeling the Transport and Fate of Sediments Released from Marine Construction Projects in the Coastal Waters of British Columbia, Canada

Major marine construction projects, resulting in the release of sediments, are subject to environmental assessment and other regulatory approval processes. An important tool used for this is the development of specialized numerical methods for these marine activities. An integrated set of numerical methods addresses four distinct topics: (1) The near-field release and mixing of suspended sediments into the water column (i.e., the initial dilution zone); (2) the transport of the suspended sediments under the influence of complex ocean currents in the far-field; (3) the settling of the transported suspended sediments onto the seabed; and (4) the potential for resuspension of the deposited sediments due to sporadic occurrences of unusually large near-bottom currents. A review of projects subjected to environmental assessment in the coastal waters of British Columbia, from the year 2006 to 2017, is presented to illustrate the numerical models being used and their ongoing development. Improvements include higher resolution model grids to better represent the near-field, the depiction of particle size dependent vertical settling rates and the computation of resuspension of initially deposited sediments, especially in relation to temporary subsea piles of sediments arising from trenching for marine pipelines. The ongoing challenges for this numerical modeling application area are also identified