On the magnitude of upwelling fluxes in shelf-break canyons

A hydrodynamic model is employed to derive the magnitude of on-shelf fluxes through a shelf-break canyon for a wide range of canyon sizes and ambient oceanic conditions. Predicted canyon-upwelling fluxes are of the order of 0.05–0.1 Sv (1 Sv=1 million m³/s), being several orders of magnitude greater than upslope fluxes in the bottom Ekman layer on the ambient continental slope. On the basis of ∼150 simulations conducted, a bulk formula of upwelling flux in a submarine canyon is derived. For typical conditions, the upwelling flux varies quadratically with forcing strength (speed of incident flow), linearly with canyon depth, and is inversely proportional to the buoyancy frequency of the density stratification inside the canyon. Other parameters such as density stratification above shelf-break depth and bottom friction are found to have minor influences on the resultant canyon-upwelling flux.     

Disentangling the upwelling mechanisms of the South Brazil Bight

This article presents a suite of long-term numerical simulations that investigate the dynamical mechanisms controlling the circulation in the South Brazil Bight (SBB). The overarching goal of these simulations is to quantify the relative contributions of local wind forcing and the Brazil Current (BC) to the upwelling of nutrient-rich slope water onto the shelf. The model results indicate that the water mass structure of the SBB is controlled by the synergy between wind-driven, inner-shelf upwelling and geostrophic, shelf-break upwelling. The later extends yearlong but the former peaks during the austral summer and decreases towards the winter. The interaction between the poleward flow of the BC and the bottom topography greatly influences the shelf circulation, particularly in the bottom boundary layer. Changes of the SBB coastline direction and shelf width modulate the along-shore pressure gradient and the magnitude of the shelf-break upwelling and downwelling. Thus, although the summer upwelling winds extend over large part of the SBB surface temperatures are warmer in the south because of the cooling effect of the shelf-break upwelling in the northern region. At difference with previous studies of shelf-break dynamics the shelf-break upwelling in our model is not controlled by the uplifting associated with the presence of instabilities of the boundary current or nonlinear accelerations under a variable shelf width. The proposed mechanism is relatively simple. As the boundary current flows along the continental slope, changes in the coastline orientation and along-shore bottom topography modify the along-shore pressure gradient which through geostrophy leads to inshore bottom flow and hence shelf-break upwelling. Such a mechanism can provide insight into upwellings on other western boundary current regions where similar topographic variations exist.     

Numerical models of continental and submarine erosion: application to the northern Ligurian Margin (Southern Alps,France/Italy)

We present a new numerical surface process model allowing us to take into account submarine erosion processes due either to submarine landslides or to hyperpycnal currents. A first set of models show that the frequency of hyperpycnal flows influences the development of submarine canyons at the mouth of continental rivers. Further experiments show that an increase in submarine slope leads to faster regressive canyon erosion and a more dentritic canyon network, whereas increasing the height of the unstable sediment pile located on the shelf break leads to wider and less dendritic canyons. The models are then applied to the western segment of the north Ligurian margin (northwestern Mediterranean), which displays numerous submarine canyons with various sizes and morphologies. From west to east, canyon longitudinal profiles as well as margin‐perpendicular profiles progressively change from moderately steep, concave‐up shapes to steeper linear to convex‐up shapes suggesting increasing eastward margin uplift. Moreover, the foot of the margin is affected by a marked slope increase with evidences of mass transport due to landslides. Numerical models which reproduce well the North Ligurian margin morphologic features indicate that the western part of the margin is submitted to rather low (i.e. 0.4 mm yr^?1) uplift and intense submarine erosion due to frequent hyperpycnal currents, whereas the eastern part bears more rapid (i.e. 0.7 mm yr^?1) uplift and has little or no hyperpycnal currents. Copyright © 2014 John Wiley & Sons, Ltd.     

Temporal evolution of lead isotope ratios in sediments of the Central Portuguese Margin: A fingerprint of human activities

Stable Pb isotope ratios (²⁰⁶Pb/²⁰⁷Pb, ²⁰⁸Pb/²⁰⁶Pb), ²¹⁰Pb, Pb, Al, Ca, Fe, Mn and Si concentrations were measured in 7 sediment cores from the west coast of the Iberian Peninsula to assess the Pb contamination throughout the last 200 years. Independently of their locations, all cores are characterized by increasing Pb/Al rends not related to grain-size changes. Conversely, decreasing trends of ²⁰⁶Pb/²⁰⁷Pb were found towards the present. This tendency suggest a change in Pb sources reflecting an increased proportion derived from anthropogenic activities. The highest anthropogenic Pb inventories for sediments younger than 1950s were found in the two shallowest cores of Cascais and Lisboa submarine canyons, reflecting the proximity of the Tagus estuary. Lead isotope signatures also help demonstrate that sediments contaminated with Pb are not constrained to estuarine–coastal areas and upper parts of submarine canyons, but are also to transferred to a lesser extent to deeper parts of the Portuguese Margin.     

Internal tides and waves near the continental shelf edge

Abstract

The subject is reviewed from the viewpoints of theory, internal tide and wave structure and their implications.

A wider theoretical context suggests scope for further investigation of natural or nearly-trapped forms above the inertial frequency.

Although internal tides in many locations are observed to have first-mode vertical structure, higher modes are seen offshore from shallow shelf-break forcing and for particular Froude numbers, and may be expected locally near generation. Bottom intensification is often observed where the sea floor matches the characteristic slope. Solitons form from internal tides of large amplitude or at large changes of depth.

Internal tides and solitons are observed also at many sills and in straits, and to intensify in canyons.

Non-linear effects of the waves, especially solitons, include the conveyance of water, nutrients, ‘‘mixing potential'’ etc. away from their source to other locations, and the generation of mean currents. The waves transfer energy and possibly heat between the ocean and shelf, may be a source of medium frequency waves on the shelf (periods of minutes) and can contribute to interior mixing and overturning, bottom stirring and sediment movement.     

An analytical model for ice-edge upwelling

Abstract

This paper presents an analytical, two-dimensional model of the wind-induced homogeneous circulation near the edge of an ice pack floating on the ocean surface. It is shown that a vertical shear layer arises under the ice edge, by which the wind-driven geostrophic motion in the open ocean is matched to the flow region underneath the ice. As in coastal upwelling models, this shear layer consists of a thin E ^1/2-layer inside a thicker E ^1/4-layer (E being the Ekman number). Under certain conditions the shear layer produces a vertical mass flux from the bottom to the surface Ekman layer. Near the surface this upwelling flux is concentrated in the narrow E ^1/2-layer. Comparison with observations of upwelling at the edge of a polar ice pack shows good agreement.     

Eastern boundary drainage of the North Atlantic subtropical gyre

The eastern boundary of the North Atlantic subtropical gyre (NASG) is an upwelling favorable region characterized by a mean southward flow. The Canary Upwelling Current (CUC) feeds from the interior ocean and flows south along the continental slope off NW Africa, effectively providing the eastern boundary condition for the NASG. We follow a joint approach using slope and deep-ocean data together with process-oriented modeling to investigate the characteristics and seasonal variability of the interior–coastal ocean connection, focusing on how much NASG interior water drains along the continental slope. First, the compiled sets of data show that interior central waters flow permanently between Madeira and the Iberian Peninsula at a rate of 2.5?±?0.6 Sv (1 Sv = 10⁶ m³ s^-1 10⁹ km s^-1), with most of it reaching the slope and shelf regions north of the Canary Islands (1.5?±?0.7?Sv). Most of the water entering the African slope and shelf regions escapes south between the easternmost Canary Islands and the African coast: In 18 out of 22 monthly realizations, the flow was southward (?0.9?±?0.4?Sv) although an intense flow reversal occurred usually around November (1.7?±?0.9?Sv), probably as the result of a late fall intensification of the CUC north of the Canary Islands followed by instability and offshore flow diversion. Secondly, we explore how the eastern boundary drainage may be specified in a process-oriented one-layer quasigeostrophic numerical model. Non-zero normal flow and constant potential vorticity are alternative eastern boundary conditions, consistent with the idea of anticyclonic vorticity induced at the boundary by coastal jets. These boundary conditions cause interior water to exit the domain at the boundary, as if recirculating through the coastal ocean, and induce substantial modifications to the shape of the eastern NASG. The best model estimate for the annual mean eastward flow north of Madeira is 3.9?Sv and at the boundary is 3.3?Sv. The water exiting at the boundary splits with 1?Sv flowing into the Strait of Gibraltar and the remaining 2.3?Sv continuing south along the coastal ocean until the latitude of Cape Ghir. The model also displays significant wind-induced seasonal variability, with a maximum connection between the interior and coastal oceans taking place in autumn and winter, in qualitative agreement with the observations.     

Banda Sea surface-layer divergence

Sea-surface temperature (SST) within the Banda Sea varies from a low of 26.5 °C in August to a high of 29.5 °C in December and May. Ekman upwelling reaches a maximum in May and June of approximately 2.5 Sv (Sv=10⁶ m³ s^?1) with Ekman downwelling at a maximum in February of approximately 1.0 Sv. The Ekman pumping annual average is 0.75 Sv upwelling. During the upwelling period, from April through December the average Ekman upwelling velocity is 2.36 × 10^?6 m s^?1 (1.27 Sv). ENSO modulation is generally within 0.5 Sv of the mean Ekman curve, with weaker (stronger) July to October upwelling during El Niño (La Niña). Combined TOPEX/POSEIDON and ERS 1993–1999 altimeter data reveal a 33 cm maximum range of sea level. Steric effects are minor, with well over 80% of the sea level change due to mass divergence (some bias due to unresolved tidal aliasing may still be present). The annual and interannual sea level behavior follows the monsoonal and ENSO phenomena, respectively. Lower (higher) sea level occurs in the southeast (northwest) monsoon and during El Niño (La Niña) events. The surface-layer volume anomaly and the surface-layer divergence, assuming a two-layer ocean, are estimated. Maximum divergence is attained during the transitional monsoon months of October/November: 1.7 Sv gain (convergence), with matching loss (divergence) in the April/May. During the El Niño growth period of 1997 the surface layer is divergent, but in 1998 when the El Niño was on the wane, the average rate of change is convergent. Surface-layer divergence attains values as high as 4 Sv. Banda Sea surface-water divergence correlates reasonably well with the 3-month lagged export of surface (upper 100?m) water into the Indian Ocean as estimated by a shallow pressure gauge array. It is concluded that the Banda Sea surface-layer divergence influences the timing and transport profile of the Indonesian throughflow export into the Indian Ocean, as proposed by Wyrtki in 1958, and that satellite altimetry may serve as an effective means of monitoring this phenomena.     

High salinity winter outflow from a mega inverse-estuary—the Great Australian Bight

Conductivity-temperature-depth (CTD) observations taken in the Great Australian Bight (GAB) during ORV Franklin cruise Fr 07/94 in July 1994 indicated the presence of a dense bottom layer at the head of the GAB, which flowed along the sea floor towards the shelf-break as a gravity current The north central region of the GAB was stratified with a maximum salinity difference of between 0.4 and 0.5. The outflow was confined to the shelf and was directed in a south-easterly direction with little evidence of cross-shelf transport. The flow exhibited a well-defined bottom interface evident from the head of the GAB to near the mouth of Spencer Gulf (SG), where the surface-bottom salinity difference was about 0.3. The mean thickness of the outflow was about 15 m. An estimate of the speed of the outflow at the discharge over the shelf-break was made using the zero entrainment assumption. This yielded a speed of <16 cm s⁻¹, which remarkably was consistent with near bottom current meter measurements (16 cm s⁻¹) on the continental shelf edge, reported south of the Eyre Peninsula. A mass budget analysis indicated that the outflow, which probably is partially maintained by the gravity current and partly by a wind-driven circulation would exist over the period, July–December, with a peak transport of about 10⁶ m³s⁻¹ (1 Sverdrup) which is approximately twenty times that of the bottom outflow from the adjoining Spencer Gulf.     

Ocean circulation on the North Australian Shelf

The ocean circulation on Australia's Northern Shelf is dominated by the Monsoon and influenced by large-scale interannual variability. These driving forces exert an ocean circulation that influences the deep Timor Sea Passage of the Indonesian Throughflow, the circulation on the Timor and Arafura Shelves and, further downstream, the Leeuwin Current. Seasonal maxima of northeastward (southwestward) volume transports on the shelf are almost symmetric and exceed 10⁶ m³/s in February (June). The associated seasonal cycle of vertical upwelling from June to August south of 8.5°S and between 124°E and 137.5°E exceeds 1.5×10⁶ m³/s across 40 m depth. During El Niño events, combined anomalies from the seasonal means of high regional wind stresses and low inter-ocean pressure gradients double the northeastward volume transport on the North Australian Shelf to 1.5×10⁶ m³/s which accounts for 20% of the total depth-integrated transport across 124°E and reduce the total transport of the Indonesian Throughflow. Variability of heat content on the shelf is largely determined by Pacific and Indian Ocean equatorial wind stress anomalies with some contribution from local wind stress forcing.     

Trace elements and stable isotopic geochemistry of an Early Cambrian chert-phosphorite unit from the lower Yurtus Formation of the Sugetbrak section in the Tarim Basin

A chert-phosphorite unit from the Sugetbrak section in the Tarim Basin was analyzed for rare earth elements (Ce, Eu), redox sensitive proxies, and carbon isotopic compositions (δ¹³C_carb and δ ¹³C_org) in the lower Yurtus Formation of the Early Cambrian period. Redox sensitive element ratios (Th/U, V/Cr, Ni/Co, and V/Sc) were employed to determine the palaeoenvironmental conditions during this period. The ratios indicated that the depositional environment of the chert-phosphorite-black shale unit ranged from suboxic to anoxic. Negative Ce and positive Eu anomalies in the chert-phosphorite assemblages of the studied Yurtus Formation indicated the existence of a redox-stratified ocean, similar to that of South China. Overturn or upward expansion of the deep water-mass probably reached the shallow marine zone after the formation of the Yurtus phosphorites. The characteristics of the negative Ce anomaly may be due to phosphoritic inheritance from the Ce-depleted signature of the overlying water column. Subsequent hydrothermal inputs and reduced detrital supplies of the deep water caused by the upwelling affected certain redox sensitive elements in the sedimentary basin. δ ¹³C_carb and δ¹³C_org negative excursions in the Yurtus chert-phosphorite unit may be related to a transgression phase when episodic basinal upwelling moved ¹²C- and P-rich waters from the pelagic basin floor to the continental shelf. Although carbon isotopic compositions in the Yurtus chert-phosphorite assemblages may have suffered from diagenetic alteration, they can be used to probe diagenetic conditions. Multi-proxy geochemical studies indicated that the δ ¹³C_carb values of the Yurtus chert-phosphorite assemblages might be considered reflections of a predominantly suboxic environment that was subsequently affected by hydrothermal inputs due to the upwelling.     

The geochemical cycle of 9Be: a reconnaissance

An initial survey has been made of the geochemistry of ⁹Be in natural waters with a view to utilizing the isotope with the cosmogenic species ¹⁰Be (t_1/2 = 1.5 m.y.) in geochronological studies of sedimentary accumulation. The mobility of⁹Be in continental waters is a strong function of pH with acid streams (pH < 6) being strongly enriched as compared to alkaline carbonate rivers. Large-scale scavenging of beryllium occurs in the estuaries of rivers containing more than about 200 pM of the element. In alkaline systems of lower concentration the element appears to behave conservatively. The net flux of beryllium to the ocean is estimated as that equivalent to an average effective river concentration of about 150 pM. Beryllium is strongly enriched in ridge crest hydrothermal solutions over ambient seawater. The flux is approximately 70% of that from the continents. The deep water concentrations of ⁹Be vary by only 20% between the North Atlantic and North Pacific demonstrating that the “chemical saltation” produced by scavenging from the deep water column and release from the sediments predominates over simple advective concentration along flow lines. Isotopic homogenization of ¹⁰Be and ⁹Be is probably achieved in the ocean. The situation in continental waters appears much less promising.     

210Po and 210Pb distributions in ocean water profiles from the Eastern South Pacific

Two ocean profiles from the Peru Basin from regions with different surface productivities were analyzed for total²¹⁰Pb and²⁰¹Po to evaluate the influence of particulates in the water column on their distribution. Comparison with a published²²⁶Ra profile for the region was made. The profile closest to the coast, where upwelling and productivity are high, shows depletion of²¹⁰Pb relative to²²⁶Ra at all depths, with particularly marked excursions from radioactive equilibrium at the surface and in the bottom water.²¹⁰Po appears to be deficient relative to²¹⁰Pb at depth as well. Mean residence times in the deep water, relative to particulate removal from the water column to the sediments, of about 100 years for²¹⁰Pb and about two years for²¹⁰Po are indicated. The profile northwest of the upwelling region shows the²²⁶Ra²¹⁰Pb²¹⁰Po system close to equilibrium at all depths to 1500 m (except for the effect of atmospheric²¹⁰Pb input seen at the surface.     

Water and nutrient fluxes off Northwest Africa

A historical data set is used to describe the coastal transition zone off Northwest Africa during spring 1973 and fall 1975, from 17° to 26°N, with special emphasis on the interaction between subtropical (North Atlantic Central Waters) and tropical (South Atlantic Central Waters) gyres. The near-surface geostrophic circulation, relative to 300 m, is quite complex. Major features are a large cyclonic pattern north of Cape Blanc (21°N) and offshore flow at the Cape Verde front. The large cyclone occurs in the region of most intense winds, and resembles a large meander of the baroclinic southward upwelling jet. The Cape Verde frontal system displays substantial interleaving that may partly originate as mesoscale features at the coastal upwelling front. Property–property diagrams show that the front is an effective barrier to all properties except temperature. The analysis of the Turner angle suggests that the frontal system is characterized by large heat horizontal diffusion as a result of intense double diffusion, which results in the smoothing of the temperature horizontal gradients. Nine cross-shore sections are used to calculate along-shore geostrophic water-mass and nutrient transports and to infer exchanges between the coastal transition zone and the deep ocean (import: deep ocean to transition zone; export: transition zone to deep ocean). These exchanges compare well with mean wind-induced transports and actual geostrophic cross-shore transport estimates. The region is divided into three areas: southern (18–21°N), central (21–23.5°N), and northern (23.5–26°N). In the northern area geostrophic import is roughly compensated with wind-induced export during both seasons. In the central area geostrophic import is greater than wind-induced export during spring, resulting in net import of both water (0.8 Sv) and nitrate (14 kmol s⁻¹), but during fall both factors again roughly cancel. In the southern area geostrophy and wind join to export water and nutrients during both seasons, they increase from 0.6 Sv and 3 kmol s⁻¹ during fall to 2.9 Sv and 53 kmol s⁻¹ during spring.     

The half life of32Si determined from a varved Gulf of California sediment core

Excess²¹⁰Pb measurements and varve chronology were used to establish a sediment accumulation rate of 0.19 cm/yr in a 95-cm-long box core raised from the Gulf of California. Varve thickness is unchanged over the entire length of the core, indicating a constant rate of sediment accumulation. The³²Si specific activity of biogenic silica shows an exponential decrease with depth in this core. The half life of³²Si, calculated from these data and the 0.19-cm/yr sediment accumulation rate, is276 ± 32 years. As most of the silica and³²Si supplied to the Gulf of California is a result of upwelling of deep ocean water, this half life determination should be relatively insensitive to secular variations in the atmospheric supply of³²Si.     

Changes in trace metals in Thalassia testudinum after hurricane impacts

Major hurricanes Emily and Wilma hit the Mexican Caribbean in 2005. Changes in trace metals in the seagrass Thalassiatestudinum prior to (May 2004, 2005) and following passage of these hurricanes (May, June 2006) were determined at four locations along a ∼130 km long stretch of coast. Before the hurricanes, essential metals were likely limiting and concentrations of potentially toxic Pb were high in a contaminated lagoon (27.5 μg g⁻¹) and near submarine springs (6.10 μg g⁻¹); the likely sources were inland sewage disposal or excessive boat traffic. After the hurricanes, Pb decreased to 2.0 μg g⁻¹ in the contaminated lagoon probably through flushing. At the northern sites, essential Fe increased >2-fold (from 26.8 to 68.3 μg g⁻¹ on average), possibly from remobilization of anoxic sediments or upwelling of deep seawater during Wilma. Thus, hurricanes can be beneficial to seagrass beds in flushing toxic metals and replenishing essential elements.     

Lithium and lithium isotope profiles through the upper oceanic crust: a study of seawater-basalt exchange at ODP Sites 504B and 896A

Ocean Drilling Program (ODP) Hole 504B near the Costa Rica Rift is the deepest hole drilled in the ocean crust, penetrating a volcanic section, a transition zone and a sheeted dike complex. The distribution of Li and its isotopes through this 1.8-km section of oceanic crust reflects the varying conditions of seawater alteration with depth. The upper volcanic rocks, altered at low temperatures, are enriched in Li (5.6-27.3 ppm) and have heavier isotopic compositions (δ⁷Li=6.6-20.8‰) relative to fresh mid-ocean ridge basalt (MORB) due to uptake of seawater Li into alteration clays. The Li content and isotopic compositions of the deeper volcanic rocks are similar to MORB, reflecting restricted seawater circulation in this section. The transition zone is a region of mixing of seawater with upwelling hydrothermal fluids and sulfide mineralization. Li enrichment in this zone is accompanied by relatively light isotopic compositions (−0.8-2.1‰) which signify influence of basalt-derived Li during mineralization and alteration. Li decreases with depth to 0.6 ppm in the sheeted dike complex as a result of increasing hydrothermal extraction in the high-temperature reaction zone. Rocks in the dike complex have variable isotopic values that range from −1.7 to 7.9‰, depending on the extent of hydrothermal recrystallization and off-axis low-temperature alteration. Hydrothermally altered rocks are isotopically light because ⁶Li is preferentially retained in greenschist and amphibolite facies minerals. The δ⁷Li values of the highly altered rocks of the dike complex are complementary to those of high-temperature mid-ocean ridge vent fluids and compatible to equilibrium control by the alteration mineral assemblage. The inventory of Li in basement rocks permits a reevaluation of the role of oceanic crust in the budget of Li in the ocean. On balance, the upper 1.8 km of oceanic crusts remains a sink for oceanic Li. The observations at 504B and an estimated flux from the underlying 0.5 km of gabbro suggest that the global hydrothermal flux is at most 8×10⁹ mol/yr, compatible with geophysical thermal models. This work defines the distribution of Li and its isotopes in the upper ocean crust and provides a basis to interpret the contribution of subducted lithosphere to arc magmas and cycling of crustal material in the deep mantle.     

On the influence of adequate Weddell Sea characteristics in a large-scale global ocean circulation model

Global ocean circulation models usually lack an adequate consideration of high-latitude processes due to a limited model domain or insufficient resolution. Without the processes in key areas of the global thermohaline circulation, the characteristics and flow of deep and bottom waters cannot be modeled realistically. In this study, a high-resolution (~20 km) ocean model focused on the Weddell Sea sector of the Southern Ocean is combined with a low-resolution (2^° × 2^°) global ocean model applying the state estimation technique. Temperature, salinity, and velocity data on two Weddell Sea sections from the regional model are used as constraints for the large-scale model in addition to satellite altimetry and sea-surface temperatures. The differences between the model with additional constraints and without document that the Weddell Sea circulation exerts significant influence on the course of the Antarctic Circumpolar Current with consequences for Southern Ocean water mass characteristics and the spreading of deep and bottom waters in the South Atlantic. Furthermore, a warming trend in the period 1993–2001 was found in the Weddell Sea and adjacent basins in agreement with float measurements in the upper Southern Ocean. Teleconnections to the North Atlantic are suggested but need further studies to demonstrate their statistical significance.     

Submarine groundwater discharge into the coast revealed by water chemistry of man-made undersea liquefied petroleum gas cavern

The occurrence of submarine groundwater discharge (SGD) as well as its supply of many nutrients and metals to coastal seawaters is now generally known. However, previous studies have focused on the chemical and radiological analysis of groundwater, surface seawater, shallow marine sediments and their pore waters, as well as the measurement of upward flow through the marine sediments, as end members of the discharge process. In this study, chemical and isotopic analysis results of marine subsurface waters are reported. These were obtained from deep boreholes of an undersea liquefied petroleum gas (LPG) storage cavern, located about 8 km off the western coast of Korea. The cavern is about 130–150 m below the sea bottom, which is covered by a 4.8–19.5 m silty clay stratum. An isotopic composition (δ²H and δ¹⁸O) of the marine subsurface waters falls on a mixing line between terrestrial groundwater and seawater. Vertical EC profiling at the cavern boreholes revealed the existence of a fresh water zone. An increase in the contents of ferrous iron and manganese and a decrease in levels of nitrate, bicarbonate and cavern seepage were recorded in August 2006, indicating a decreased submarine groundwater flux originating from land, mainly caused by an elevated cavern gas pressure. It is suggested in this study that the main source of fresh waters in the man-made undersea cavern is the submarine groundwater discharge mainly originating from the land.     

Mesoscale dynamics and walleye pollock catches in the Navarin Canyon area of the Bering Sea

Large canyons incise the shelf break of the eastern Bering Sea to be preferred sites of the cross-shelf exchange. The mesoscale eddy activity is particularly strong near the shelf-break canyons. To study the mesoscale dynamics in the Navarin Canyon area of the Bering Sea, the time series of velocities derived from AVISO satellite altimetry between 1993 and 2015, drifters, Argo buoys, and ship-borne data are analyzed. We demonstrate that the strength of anticyclonic eddies along the shelf edge in spring and summer is determined by the wind stress in March–April. The increased southward wind stress in the central Bering Sea forced a supply of low-temperature and low-salinity outer shelf water to the deep basin and formation of the anticyclonic mesoscale circulation seaward of the Navarin Canyon. Enhanced northwestward advection of the Bering Slope Current water leads to increase in an ice-free area in March and April and increased bottom-layer temperature at the outer shelf. The strong (weak) northwestward advection of the eastern Bering Sea waters, determined by eastern winds in spring, creates favorable (unfavorable) conditions for the pollock abundance in the western Navarin Canyon area in summer.