Ocean current estimation using a Multi-Model Ensemble Kalman Filter during the Grand Lagrangian Deployment experiment (GLAD)

In the summer and fall of 2012, during the GLAD experiment in the Gulf of Mexico, the Consortium for Advanced Research on Transport of Hydrocarbon in the Environment (CARTHE) used several ocean models to assist the deployment of more than 300 surface drifters. The Navy Coastal Ocean Model (NCOM) at 1 km and 3 km resolutions, the US Navy operational NCOM at 3 km resolution (AMSEAS), and two versions of the Hybrid Coordinates Ocean Model (HYCOM) set at 4 km were running daily and delivering 72-h range forecasts. They all assimilated remote sensing and local profile data but they were not assimilating the drifter’s observations. This work presents a non-intrusive methodology named Multi-Model Ensemble Kalman Filter that allows assimilating the local drifter data into such a set of models, to produce improved ocean currents forecasts. The filter is to be used when several modeling systems or ensembles are available and/or observations are not entirely handled by the operational data assimilation process. It allows using generic in situ measurements over short time windows to improve the predictability of local ocean dynamics and associated high-resolution parameters of interest for which a forward model exists (e.g. oil spill plumes). Results can be used for operational applications or to derive enhanced background fields for other data assimilation systems, thus providing an expedite method to non-intrusively assimilate local observations of variables with complex operators. Results for the GLAD experiment show the method can improve water velocity predictions along the observed drifter trajectories, hence enhancing the skills of the models to predict individual trajectories.     

Dispersion of surface drifters and model-simulated trajectories

From a data set encompassing the years 1990–2008 pairs of surface drifters with maximum initial separations of 5, 10 and 25 km have been identified. Model trajectories have been calculated using the same initial positions and times as the selected pairs of surface drifters. The model trajectories are based on the TRACMASS trajectory code and driven by the ocean general circulation model NEMO. The trajectories are calculated off-line, i.e. with the stored velocity fields from the circulation model. The sensitivity of the trajectory simulations to the frequency of the stored velocity fields was tested for periods of 3 and 6 h as well as 5 days. The relative dispersion of the surface-drifter and model trajectories has been compared, where the latter was found to be too low compared to the relative dispersion of the drifters.Two low-order trajectory sub-grid parameterisations were tested and successfully tuned so that the total amplitude of the relative dispersion of the model trajectories is similar to that associated with the drifter trajectories. These parameterisations are, however, too simple for a correct simulation of Lagrangian properties such as the correlation time scales and the variance of the eddy kinetic energy.The importance of model-grid resolution is quantified by comparing the relative dispersion from an eddy-permitting and a coarse-resolution model, respectively. The dispersion rate is halved with the coarse grid. The consequences of the two-dimensionality of the trajectories is evaluated by comparing the results obtained with the 2D and the Lagrangian 3D trajectories. This shows that the relative dispersion is 15% stronger when the trajectories are freely advected with the 3D velocity field.     

How much data is enough? The importance of morphological sampling interval and duration for calibration of empirical shoreline models

The ability to robustly predict future shoreline position under the influence of changing waves and sea-level rise is a key challenge to scientists and engineers alike. While extrapolating a linear trend out in time is a common baseline approach, the recent development of a number of empirical shoreline models allows the prediction of storm and annual-scale variability as well. The largest constraint in applying these models is the availability of high quality, adequate duration data sets in order to calibrate model free parameters. This contribution outlines several such models and discusses the monitoring programs required to calibrate and hindcast shoreline change from 1 to 10 years at two distinct beach types: a storm-dominated site and the second exhibiting a large seasonal variability. The seasonally-dominated site required longer data sets but was less sensitive to sampling interval, while the storm-dominated site converged on shorter, more frequently sampled data sets. In general, calibration based on a single year of observed shorelines resulted in a large range of model skill and was not considered robust. Monitoring programs of at least two years, with shorelines sampled at dt ≤ 30 days were sufficient to determine initial estimates of calibration coefficients and hindcast short-term (1–5 years) shoreline variability. In the presence of unresolved model processes and noise, hindcasting longer (5 + years) data sets required longer (5 + years) calibration data sets, particularly when sampling intervals exceeded 60 days.     

Concept analysis and laboratory observations on a water piercing missile launcher

The water piercing missile launcher (WPML) is a new concept for launching missiles from submerged platforms. The WPML employs a high speed gas jet, using rocket exhaust as the gas source, to create a dry path underwater through which a missile may pass without contacting water. The gas jet is deflected due to cross flow and the gas jet trajectory is computed through a semi-empirical relationship commonly used to describe single phase jet deflection. This relationship, which requires an experimentally determined constant to predict jet deflection, is computed using experimental data. Uncoupled simulations of rocket exhaust and missile dynamics are shown to demonstrate how such a launcher could be utilized in launching a generic artillery missile. Although the results indicate the optimal launch depth, in terms of maximizing the launch depth while minimizing missile restraint time, is 14 m given an 8 m/s submarine speed it may be possible to launch a missile from a moving submarine at a speed of 5 m/s at a maximum depth of 20 m.     

Agulhas leakage into the Atlantic estimated with subsurface floats and surface drifters

Surface drifters and subsurface floats drifting at depths near 800 m were used to study the pathways of warm, salty Indian Ocean water leaking into the South Atlantic that is a component of the upper limb of the Atlantic meridional overturning circulation (MOC). Four drifters and 5 floats drifted from the Agulhas Current directly into the Benguela Current. Others looped for various amounts of time in Agulhas rings and cyclones, which translated westward into the Atlantic, contributing a large part of Indian Ocean leakage. Agulhas rings translated into the Benguela Current, where they slowly decayed. Some large, blob-like Agulhas rings with irregular shapes were found in the southeastern Cape Basin. Drifter trajectories suggest these rings become more circular with time, eventually evolving into the circular rings observed west of the Walvis Ridge. Agulhas cyclones, which form on the north side of the Agulhas Current south of Africa, translated southwestward (to 6°E) and contributed water to the southern Cape Basin. A new discovery is a westward extension from the mean Agulhas retroflection measured by westward drifting floats near 41°S out to at least 5°W, with some floats as far west as 25°W. The Agulhas extension appears to split the South Atlantic Current (SAC) into two branches and to transport Agulhas water westward, where it is mixed and blended with eastward-flowing water from the western Atlantic. The blended mixture flows northeastward in the northern branch of the SAC and into the Benguela Current. Agulhas leakage transport was estimated from drifters and floats to be at least 15 Sv in the upper 1000 m, which is equivalent to the transport of the upper layer MOC. It is suggested that the major component of the upper layer overturning circulation in the Atlantic is Agulhas leakage in the form of Agulhas rings.     

On the detectability of internal tides in Drake Passage

The Southern Ocean hosts significant topographic mixing that might be associated with internal tides. Tidal signals are evident in bottom temperature at 1000 m in Drake Passage, suggesting that internal tides with an amplitude of between ∼20 and 200 m may be present. Various necessary conditions for internal tide generation show that the steep topography in and around Drake Passage can initiate internal tides, and recent global tide models have suggested this region to generate very large interface displacements. Here, we present an attempt to detect internal tides in Drake Passage. During the last 10 years, combinations of bottom pressure recorders and inverted echo sounders have been deployed in the region. The bottom pressure recorders measure predominantly the barotropic tide; the inverted echo sounders measure travel time from sea bed to sea surface and therefore are influenced both by sea level (barotropic tide) and internal sound speed (internal tide). By subtracting one from the other, the internal tide should be detectable. Although the technique works successfully around Hawaii, it does not prove the existence of large internal tides in Drake Passage. The detectability of the internal tidal signal in Drake Passage is investigated using a six-layer one-dimensional model to simulate the bottom pressure and travel time signals of a semi-diurnal tide. The temperature and salinity stratification in Drake Passage is sufficiently weak that large vertical excursions are necessary to produce a signal in travel time detectable above the noise in Drake Passage. An internal tide of at least 70 or 20 m in northern and southern Drake Passage, respectively, would be detected. The fact that these are, perhaps surprisingly, not detected by the combination of bottom pressure and travel time, constrains the internal tides in Drake Passage to be ∼20 m in southern Drake Passage, and between 20 and 70 m in northern Drake Passage. The model also predicts that satellite altimetry would not be able to detect internal tides in Drake Passage, but would in the Brazil Basin and Hawaii regions.     

Directed drifter launch strategies for Lagrangian data assimilation using hyperbolic trajectories

The dependence of the fidelity of a Lagrangian data assimilation scheme on the initial launch locations of the observed drifters is studied in the context of a reduced gravity, primitive equation model of mid-latitude circulations. A directed launch strategy, based on tracking the Lagrangian manifolds emanating from strongly hyperbolic regions in the flow field, is developed. In a series of twin assimilation experiments, the rate of convergence of the data assimilating scheme is shown to be consistently higher in such directed launches compared to those produced from randomly selected initial drifter positions. By directing initial drifter positions along the out-flowing branch of identifiable Lagrangian boundaries, the relative dispersion of the drifters, the overall data coverage and the sampling of high kinetic energy features in the flow are optimized. In general, the performance of the assimilation procedure is shown to depend strongly on the independence of the observed drifter trajectories and the temporal persistence of the corrections provided by the data.     

The relative contributions of unicellular and filamentous diazotrophs to N2 fixation in the South China Sea and the upstream Kuroshio

We studied the seasonal, diel, and vertical distribution of phytoplankton N₂ fixation to understand the relative contributions of unicellular and filamentous nitrogen fixers (diazotrophs) to N₂ fixation and nitrogen recycling in the northern South China Sea (SCS) and the neighboring upstream Kuroshio. N₂-fixation rates were measured by the ¹⁵N₂ tracer technique (addition by bubble) on unicellular (<10 or 20 µm) and the filamentous diazotrophs (>10 or 20 µm, mostly Trichodesmium and Richelia) fractionated by 10- or 20-µm mesh sizes. The mean depth-integrated total (unicellular+filamentous) N₂-fixation rates in the SCS (51.7±6.2 µmol N m⁻² d⁻¹) averaged 1/3 of that in the Kuroshio (142.7±29.6 µmol N m⁻² d⁻¹), with higher rates in the winter than in other seasons in the SCS and the opposite seasonal pattern in the Kuroshio. Unicellular diazotrophs contributed 65% of the total N₂ fixation in the SCS, which were negatively correlated with surface temperature and, as for total N₂ fixation, were higher in the winter when Trichodesmium spp. were scarce. In comparison, the unicellular diazotrophs contributed 50% of total N₂ fixation in the Kuroshio, and their contributions were not significantly correlated with surface temperature. In both the SCS and the Kuroshio, the unicellular N₂ fixation was more important during the night than during the day, and in the deep euphotic layer than in the surface layer, even in the daytime. Our results show that the unicellular diazotrophs were important N₂ fixers and contributed significantly to N₂ fixation in the tropical marginal seas, more so in the SCS than the Kuroshio.     

A probabilistic method for constructing wave time-series at inshore locations using model scenarios

Continuous time-series of wave characteristics (height, period, and direction) are constructed using a base set of model scenarios and simple probabilistic methods. This approach utilizes an archive of computationally intensive, highly spatially resolved numerical wave model output to develop time-series of historical or future wave conditions without performing additional, continuous numerical simulations. The archive of model output contains wave simulations from a set of model scenarios derived from an offshore wave climatology. Time-series of wave height, period, direction, and associated uncertainties are constructed at locations included in the numerical model domain. The confidence limits are derived using statistical variability of oceanographic parameters contained in the wave model scenarios. The method was applied to a region in the northern Gulf of Mexico and assessed using wave observations at 12 m and 30 m water depths. Prediction skill for significant wave height is 0.58 and 0.67 at the 12 m and 30 m locations, respectively, with similar performance for wave period and direction. The skill of this simplified, probabilistic time-series construction method is comparable to existing large-scale, high-fidelity operational wave models but provides higher spatial resolution output at low computational expense. The constructed time-series can be developed to support a variety of applications including climate studies and other situations where a comprehensive survey of wave impacts on the coastal area is of interest.     

Near-surface circulation of the northeast Pacific Ocean derived from WOCE-SVP satellite-tracked drifters

Statistics of the near-surface circulation in the northeast Pacific Ocean were derived from the trajectories of nearly 100 surface drifters tracked between August 1990 and December 1995 as part of the World Ocean Circulation Experiment's (WOCE) Surface Velocity Program (SVP). Drifters were drogued within the mixed layer (15 m drogue depth) or near the top of the permanent halocline (120 m). All branches of the Alaskan Gyre were well-sampled at both depths, revealing a weak Subarctic Current, a bifurcation of the Subarctic Current near 48°N, 130°W at 15 m depth, and strong, variable flow in the Alaska Current and Alaskan Stream. At 120 m depth, northward flow in the Alaska Current occurred much farther offshore than within the mixed layer. The drifter trajectories revealed interannual variability, with evidence of an intensified Alaskan Gyre during the winters of 1991–92 and 1992–93 and more southerly transport during winter 1994–95. A minimum in eddy kinetic energy was found at both depths within the northern branch of the Subtropical Gyre. Eddy kinetic energies were nearly twice as high in the mixed layer compared to below, and were 2–3 times larger in winter than in summer throughout most of the near-surface Alaskan Gyre. High eddy energies observed near the eastern perimeter of the Alaskan Gyre may be due to the offshore intrusion of eddies formed by coastal current instabilities.Taylor's theory of single-particle dispersion was applied to the drifter ensembles to estimate Lagrangian decorrelation scales and eddy diffusivities. Both the initial dispersion and random walk regimes were identified in the dispersion time series computed for several regions of both ensembles. The integral time scales and eddy diffusivities computed from the dispersion scale linearly with r.m.s. velocity, which is consistent with drifter studies from the Atlantic. An exception is the meridional integral time scales, which were nearly constant throughout the study area and at both drogue depths. The magnitudes of the derived eddy statistics are comparable to those derived from surface drifters in other parts of the world ocean. These are the first Lagrangian estimates of particle dispersion over a broad region of the near-surface North Pacific, and the consistency of the results with previous studies from the Atlantic lends credence to the idea that the simplifying assumptions of Taylor (1921) (Proceedings of the London Mathematical Society Series A 20, 196–221) are reasonably valid throughout the upper ocean. This bodes well for the effective parameterization of near-surface diffusivities in general circulation models. Finally, the drifter-derived velocity statistics were used to speculate on the source regions of waters of possible coastal origin observed at offshore stations during the field studies of the Canadian Joint Global Ocean Flux Study.     

Silicon flux and distribution of biogenic silica in deep-sea sediments in the western North Pacific Ocean

We investigated biogenic silica, several biological components, and silicate in pore-water in the abyssal sediment to determine silicon flux of western North Pacific during several cruises. The surficial sediment biogenic silica content was high at high latitudes with the boundary running along the Kuroshio Extension, and maximum values (exceeding 20%) were found in the Oyashio region. In the subtropical region to the south, most stations showed less than 5% biogenic silica content. This distribution pattern reflected primary production and ocean currents in the surface layer very well. Pore-water samples were collected from 4 stations along the east coast of Japan. The highest asymptotic silicic acid concentration (670 μmol L^?1) in pore-water was observed at the junction of Kuroshio and Oyashio, followed by samples from the Oyashio region. It is at the southern station that the lowest value (450 μmol L^?1) was observed, and the primary production is low under the influence of Kuroshio there. The diffusive flux followed the same geographic trend as the asymptotic silicic acid concentrations did, ranging 77–389 mmol m^?2 yr ^?1. Multiple sampling of pore-water was conducted throughout the year at one station at high latitude. The average annual biogenic silica rain flux observed using sediment traps was 373 mmol m^?2 yr^?1; the diffusive flux and burial flux at the sediment–water interface were 305 and 9 mmol m^?2 yr^?1, respectively. We concluded that most of the settling silica particles dissolved and diffused at the sediment–water interface and approximately 3% only were preserved in this area. In addition, the obvious time lag observed between the peak rain flux and the maximum diffusive flux suggested that primary production in the surface layer has a great influence on the sedimentation environment of abyssal western North Pacific. These transitions of Si flux at the sediment–water interface were considerably greater in northwestern North Pacific than in southwestern North Pacific. In addition, a station in the Philippine Sea indicated high biogenic silica content because of Ethmodiscus ooze, which are scattered randomly on the sea floor in the subtropical region.     

The swarm dynamics of northern krill (Meganyctiphanes norvegica) and pteropods (Cavolinia inflexa) during vertical migration in the Ligurian Sea observed by an acoustic Doppler current profiler

A ship-mounted 153 kHz narrow-band ADCP and 1 m² MOCNESS were deployed between 16 and 24 Sept. 1997 in the Ligurian central zone (∼43°20′N 7°48′E). Results from both instruments showed that the zooplankton community performed vertical migrations that conformed to the classical pattern of ascent at dusk (∼18:30 h) and descent at dawn (∼06:30 h). Depth-discrete net samples between 0 and 500 m showed that the community was dominated by two species, the euphausiid Meganyctiphanes norvegica (Northern krill) and the pteropod Cavolinia inflexa, which migrated in separate discrete bands that were detectable by the ADCP. Information from the ADCP was used to estimate vertical migration speed in two ways: (i) from the trajectory of the back-scattering bands over time and (ii) from the Doppler-shift vertical velocity measured within depth zones at the corresponding time and depth of these bands. Estimates of the migration speed of C. inflexa were between 2 and 7 cm s⁻¹ upwards and between 4 and 7 cm s⁻¹ downwards. M. norvegica was estimated to migrate at speeds between 7 and 8 cm s⁻¹ upwards and over 11 cm s⁻¹ downwards. The consistently lower migration speeds estimated from Doppler measurements as compared with estimates obtained from measuring trajectories of back-scattering bands over time was believed to result from a methodological artefact. The Doppler measurements were nevertheless useful in a relative sense in revealing the relative speed of individuals within swarms. It was shown that individuals at the front of the upwardly migrating band of M. norvegica moved more slowly than those at the rear. These results illustrate the extra biological information that can be obtained by ADCPs compared with conventional echo-sounders.     

The importance of tides for sediment dynamics in the deep sea—Evidence from the particulate-matter tracer 234Th in deep-sea environments with different tidal forcing

Key aspects of deep-ocean fluid dynamics such as basin-scale (residual) and tidal flow are believed to have changed over glacial/interglacial cycles, with potential relevance for climatic change. To constrain the mechanistic links, magnitudes and temporal succession of events analyses of sedimentary paleo-records are of great importance. Efforts have been underway for some time to reconstruct residual-flow patterns from sedimentary records. Attempts to reconstruct tidal flow characteristics from deep-sea sediment deposits, however, are at a very early stage and first require a better understanding of the reflection of modern tides in sediment dynamics. In this context internal (baroclinic) tides, which are formed by the surface (barotropic) tide interacting with seafloor obstacles, are believed to play a particularly important role. Here we compare two modern deep-sea environments with respect to the effect of tides on sediment dynamics. Both environments are influenced by kilometre-scale topographic features but with vastly different tidal forcing: (1) two sites in the Northeast Atlantic (NEA) being surrounded by, or located downstream of, fields of short seamounts (maximum barotropic tidal current velocities ～5 cm s^?1); and (2) a site next to the Anaximenes seamount in the Eastern Mediterranean (EMed) (maximum barotropic tidal current velocities ～0.5 cm s^?1). With respect to other key fluid-dynamical parameters both environments are very similar. Signals of sedimentary particle dynamics, as influenced by processes taking place in the bottom boundary layer, were traced by the vertical water-column distribution of radioactive disequilibria (daughter/parent activity ratios≠1) between the naturally occurring, short-lived (half-life: 24.1 d) particulate-matter tracer ²³⁴Th relative to its very long-lived and non-particle-reactive parent nuclide ²³⁸U. Activity ratios of ²³⁴Th/²³⁸U<1 in water samples collected near the seafloor indicate fast ²³⁴Th scavenging onto particles followed by fast settling of these particles from the sampled parcel of water and, therefore, imply active sediment resuspension and dynamics on time scales of up to several weeks. In the Northeast Atlantic study region tides (in particular internal tides) are very likely to locally push total current velocities near the seafloor across the critical current velocity threshold for sediment erosion or resuspension whereas in the Eastern Mediterranean the tides are much too weak for this to happen. This difference in tidal forcing is reflected in a difference of the frequency of the occurrence of radioactive disequilibria <1 between total ²³⁴Th and ²³⁸U: In the near-bottom water column of the Northeast Atlantic region 59% of samples had detectable ²³⁴Th/²³⁸U disequilibria whereas at the Eastern Mediterranean site this fraction was only 7% (including a few disequilibria >1). The results of this study, therefore, add to the evidence suggesting that tides in the deep sea of the open oceans are more important for sediment dynamics than previously thought. It is hypothesised that (a) tide/seamount interactions in the deep open ocean control the local distribution of erosivity proxies (e.g., distributions of sediment grain sizes, heavy minerals and particle-reactive radionuclides) in sedimentary deposits and (b) the aforementioned topographically controlled sedimentary imprints of (internal) tides are useful in the reconstruction of past changes of tidal forcing in the deep sea.     

Mode-1 internal tides in the western North Atlantic Ocean

Mode-1 internal tides were observed the western North Atlantic using an ocean acoustic tomography array deployed in 1991–1992 centered on 25°N, 66°W. The pentagonal array, 700-km across, acted as an antenna for mode-1 internal-tides. Coherent internal-tide waves with O(1 m) displacements were observed traveling in several directions. Although the internal tides of the region were relatively quiescent, they were essentially phase locked over the 200–300 day data record lengths. Both semidiurnal and diurnal internal waves were detected, with wavenumbers consistent with those calculated from hydrographic data. The M₂ internal-tide energy flux was estimated to be about 70 W m⁻¹, suggesting that mode-1 waves radiate 0.2 GW of energy, with large uncertainty, from the Caribbean island chain at this frequency. A global tidal model (TPXO 5) suggested that 1–2 GW is lost from the M₂ barotropic tide over this region, but the precise value was uncertain because the complicated topography makes the calculation problematic. In any case, significant conversion of barotropic to baroclinic tidal energy does not occur in the western North Atlantic basin. It is apparent, however, that mode-1 internal tides have very weak decay and retain their coherence over great distances, so that ocean basins may be filled up with such waves. Observed diurnal amplitudes were an order of magnitude larger than expected. The amplitude and phase variations of the K₁ and O₁ constituents observed over the tomography array were consistent with the theoretical solutions for standing internal waves near their turning latitude. The energy densities of the resonant diurnal internal waves were roughly twice those of the barotropic tide at those frequencies.     

Eddies and a mesoscale deflection of the slope current in the Faroe–Shetland Channel

The mesoscale dynamics of the Scottish side of the Faroe–Shetland Channel have been investigated using synoptic in situ and remote sensing observations. A cold core cyclonic eddy, identified from an AVHRR image, had a diameter of about 50 km and surface current speeds of up to 50 cm s^-1; it appeared to be attached to the 800 m isobath as it moved north-eastward along the edge of the channel at about 8 cm s^-1. Speeds in the slope current were about 50 cm s^-1 but increased to 70 cm s^-1 where the current was compressed by the eddy. Offshore, over the 1000 m isobath in the cooler water, speeds in the current were slower (ca. 20 cm s^-1). North-west of the Shetlands the offshore edge of the slope current was deflected across the channel for a distance of about 70 km from the shelf edge. The speed of drifters in the slope current increased to over 60 cm s^-1 as they moved anti-cyclonically around this deflection. CTD profiles suggest that the movement of the surface waters was mirrored in the deep water of the channel. The deflection carried a very large quantity of North Atlantic Water into the central part of the channel; its cause and ultimate fate are not known, although it is likely to have had a significant impact on the dynamics of the channel.     

Biogeochemical patterns in the Atlantic Inflow through the Strait of Gibraltar

The effects of tidal forcing on the biogeochemical patterns of surface water masses flowing through the Strait of Gibraltar are studied by monitoring the Atlantic Inflow (AI) during both spring and neap tides. Three main phenomena are defined depending on the strength of the outflowing phase predicted over the Camarinal Sill: non-wave events (a very frequent phenomenon during the whole year); type I Internal wave events (a very energetic event, occurring during spring tides); and type II Internal wave events (less intense, occurring during neap tides).During neap tides, a non-wave event comprising oligotrophic open-ocean water from the Gulf of Cádiz is the most frequent and clearly dominant flow through the Strait. In this tidal condition, the inflow of North Atlantic Central Water (NACW) provides the main nutrient input to the surface layer of the Alboran Sea, supplying almost 70% of total annual nitrate transport to the Mediterranean basin. A low percentage of active and large phytoplankton cells and low average concentrations of chlorophyll (0.3–0.4 mg m⁻³) are found in this tidal phase. Around 50% of total annual phytoplankton biomass transport into the Mediterranean Sea through the Strait presents these oligotrophic characteristics.In contrast, during spring tides, patches of water with high chlorophyll levels (0.7–1 mg m⁻³) arrive intermittently, and these are recorded concurrently with the passage of internal waves coming from the Camarinal Sill (type I internal wave events). When large internal waves are arrested over the Camarinal Sill this implies strong interfacial mixing and the probable concurrent injection of coastal waters into the main channel of the Strait. These processes result in a mixed water column in the AI and can account for around 30% of total annual nitrate transport into the Mediterranean basin. Associated with type I internal wave events there is a regular inflow of large and active phytoplankton cells, transported in waters with relatively high nutrient concentrations, which constitutes a significant supply of planktonic resources to the pelagic ecosystem of the Alboran Sea (almost 30% of total annual phytoplankton biomass transport).     

Commercially-cultured oysters (Crassostrea gigas) exert top-down control on intertidal pelagic resources in Willapa Bay,Washington, USA

The capacity of filter feeders to reduce seston and phytoplankton concentrations in the water column has important implications for restoration and management of coastal ecosystems. We directly measured changes in chlorophyll a concentration on commercially stocked intertidal oyster beds (Crassostrea gigas) in Willapa Bay, Washington, USA by recording water properties near small drifters as they tracked parcels of water across tide flats. Chlorophyll declined 9.6% per half hour in water passing on-bottom adult oysters and 41% for longline adult oysters, whereas chlorophyll concentrations increased as water flowed across tide flats without adult oysters. Field filtration rates, which were fit to exponential declines in chlorophyll and accounted for oyster density and water depth, averaged 0.35 L g^− 1 h^− 1 (shucked dry weight) for on-bottom aquaculture and 0.73 L g^− 1 h^− 1 for longline culture, compared to values of 2.5–12 L g^− 1 h^− 1 reported from laboratory studies of C. gigas. Field filtration rates may be lower than laboratory rates due to unfavorable field conditions (e.g., low initial chlorophyll concentrations) or masked by resuspension of benthic microalgae. In addition to distinctions among on-bottom, longline, and no-oyster habitats, Akaike's Information Criterion analysis showed temperature, initial chlorophyll concentration, and depth related to chlorophyll decline. This research corroborates mathematical models suggesting that benthic suspension feeders are exerting top-down control of pelagic production in this estuary, with strong patterns in chlorophyll emerging across extensive tideflats populated by C. gigas despite low field filtration rates.     

Evolution of the Kuroshio Tropical Water from the Luzon Strait to the east of Taiwan

This study examines the evolution of the Kuroshio Tropical Water (KTW) from the Luzon Strait to the I-Lan Ridge northeast of Taiwan. Historical conductivity temperature depth (CTD) profiles are analyzed using a method based on the calculation of the root mean square (rms) difference of the salinity along isopycnals. In combination with analysis of the distribution of the salinity maximum, this method enables water masses in the Kuroshio and the vicinity, to be tracked and distinguished as well as the detection of the areas where water masses are modified. Vertical and horizontal eddy diffusivities are then calculated from hydrographic and current velocity data to elucidate the dynamics underlying the KTW interactions with the surrounding water masses. Changes in KTW properties mainly occur in the southern half of the Luzon Strait, while moderate variations are observed east of Taiwan on the right flank of the Kuroshio. In spite of a front dividing the KTW from the South China Sea Tropical Water (SCSTW) on Kuroshio׳s western side, mixing between these two water masses seemingly occurs in the Luzon Strait. These water masses׳ interaction is not evident east of Taiwan. The estimation of eddy diffusivities yields high horizontal diffusivities (K_h~10² m² s⁻¹) all along the Kuroshio path, due to the high current shear along the Kuroshio׳s flanks. The vertical diffusivity approaches 10⁻³ m² s⁻¹, with the highest values in the southern Luzon Strait. Instabilities generated when the Kuroshio encounters the rough topography of this region may enhance both vertical and horizontal diffusivities there.     

Wave climate variability in the North-East Atlantic Ocean over the last six decades

Ocean surface gravity waves play a major role in many engineering and environmental problems, both in the open ocean and in coastal zones. Therefore, it is essential to improve our knowledge on spatial and temporal variability of wave climate. This study aims at investigating this variability in the North-East Atlantic Ocean (25°W–0°W and 30°N–60° N), using a 57-year hindcast (1953–2009) obtained with a spectral wave model forced with reanalysis wind fields. The hindcast analysis reveals firstly strong seasonal fluctuations of wave climate, with winters characterized by large and long-period waves of mean direction spreading from south-west to north-west, and summers characterized by smaller and shorter-period waves originating from norther directions. From northern (55°N) to southern (35°N) latitudes, the significant wave height (Hs) decreases by roughly 40%, the mean wave direction (Mwd) rotates clockwise by about 25% while the peak period (Tp) only grows by 5%. These three parameters also exhibit a strong inter-annual variability, particularly when winter-means (from 1st of December to 1st of April) are considered. Linear trend analysis over the studied period shows spatially variable long-term trends, with a significant increase of Hs (up to 0.02 m yr^?1) and a counterclockwise shift of Mwd (up to ?0.1° yr^?1) at northern latitude, contrasting with a fairly constant trend for Hs and a clockwise shift of Mwd (up to +0.15° yr^?1) at southern latitudes. Long-term trends of Tp are less significant, with still a slight increase in the north-eastern part of the study area (up to +0.01 s yr^?1). Eventually, a comparison between the inter-annual variability of the winter-means of the three selected wave parameters and the North Atlantic Oscillation (NAO) reveals: (1) a strong positive correlation between Hs and the NAO index at northern latitudes (correlation coefficient up to R = 0.91) and a significant negative correlation at southern latitudes (up to R = ?0.6); (2) no significant correlation for Mwd north of 40°N and a clear positive correlation southward of 40°N (up to R = 0.8) and (3) a northward increasing correlation for Tp (up to R = 0.8). Long-term trends for Hs, Mwd and Tp are finally explained by a significant increase in the NAO index over the studied period.     

The Kuroshio Extension and its recirculation gyres

This paper reports on the strength and structure of the Kuroshio Extension and its recirculation gyres. In the time average, quasi-permanent recirculation gyres are found to the north and south of the Kuroshio Extension jet. The characteristics of these recirculations gyres are determined from the combined observations from the Kuroshio Extension System Study (KESS) field program (June 2004–June 2006) and include current meters, pressure and current recording inverted echo sounders, and subsurface floats. The position and strength of the recirculation gyres simulated by a high-resolution numerical model are found to be consistent with the observations. The circulation pattern that is revealed is of a complex system of multiple recirculation gyres that are embedded in the crests and troughs of the quasi-permanent meanders of the Kuroshio Extension. At the location of the KESS array, the Kuroshio Extension jet and its recirculation gyres transport of about 114 Sv. This represents a 2.7-fold increase in the transport of the current compared to the Kuroshio's transport at Cape Ashizuri before it separates from the coast and flows eastward into the open ocean. This enhancement in the current's transport comes from the development of the flanking recirculation gyres. Estimates from an array of inverted echo sounders and a high-resolution ocean general circulation model are of similar magnitude.