Interannual variability of the shelf-slope front position between 75° and 50° W

The shelf-slope front (SSF) is a continuous shelf-break front running from the Tail of the Grand Banks to Cape Hatteras, North Carolina, separating colder and less-saline continental shelf waters from warmer and more saline slope waters. Time series containing mean monthly SSF positions were produced along each of 26 longitude lines between 75° and 50°W by workers located at Bedford Institute of Oceanography by digitizing individual frontal charts and computing mean monthly latitudinal positions over a 29-year (1973–2001) period. After removing seasonal variability at each longitude, interannual variability (IAV) of the SSF position at each longitude was computed as the annual mean of all monthly SSF position anomalies for each year over the 29-year period. Despite some missing data, a longitude-time plot reveals alternating bands of offshore (late-1970s, late-1980s, late-1990s) and onshore (early-1980s, early-1990s, early-2000s) annual mean SSF anomaly values, exhibiting a period of approximately 10 years. Annual mean SSF anomaly amplitudes are largest in the east, with maxima of O (± 100 km) located east of 60° W for years when data are available. Empirical orthogonal function (EOF) modes 1–4 (accounting for > 90% of the variance) form a set of basis functions that describe the SSF anomaly data and allow reconstruction of the entire data set since missing data are relatively few (14%). A complex empirical orthogonal function (CEOF) analysis using the “reconstructed” data reveals a wavelength scale of approximately 20° of longitude, a distance nearly equal to the entire study domain, along with steady, westward phase propagation of SSF anomalies over approximately the same distance. Speed calculations for the westward-propagating features yield a value of approximately 1.2 to 2.4 cm s^− 1 (1 to 2 km d^− 1), with annual mean SSF anomalies thus requiring about 4 years to propagate from the Tail of the Grand Banks in the east to Cape Hatteras, North Carolina, in the west. This propagation speed and the timing of the SSF positional anomalies at the Tail of the Grand Banks for the 29-year study period are in agreement with speeds computed for the propagation of quasi-decadal salinity anomalies through the Labrador Sea and the time of their arrival at the Tail of the Grand Banks. The small westward SSF anomaly propagation speed is an order of magnitude smaller than the associated currents, in agreement with a highly damped flow-through system originating from both Davis Strait and the West Greenland Current as discussed by other workers. Observations from both southern and northern portions of the study domain, within both continental shelf and slope waters, show that interannual changes in the volume of shelf water along with shelf water bulk properties exhibit a strong relationship with IAV of the SSF position over long time periods.     

Validation of a hybrid optimal interpolation and Kalman filter scheme for sea surface temperature assimilation

A hybrid data assimilation scheme designed for operational assimilation of satellite sea surface temperatures (SST) into an ocean model has been developed and validated against in-situ observations. The scheme consists of an optimal interpolation (OI) part and a greatly simplified Kalman filter (KF) part.The OI is performed only in the longitudinal and latitudinal directions. A climatological field is used as a background field for the interpolation. It is constructed by fitting daily averages of satellite SST to the annual mean, annual, and semiannual harmonics in a 20 km by 20 km grid. The background error covariance is approximated by a spatially varying two-dimensional exponential covariance model. The parameters of the covariance model are fitted to the deviations of the satellite data from the background field using data from a full year.The simplified KF uses ocean model forecasts as a background field. It is based on the assumption that it is possible to neglect horizontal SST covariances in the filter and that the typical time scale for vertical mixing in the mixed layer is much shorter than the average time between observations. We therefore assume that the error variance in a column of water is evenly spread out throughout the mixed layer. The result of these simplifications is a computationally very efficient KF.A one year validation of the scheme is performed for year 2001 using an operational eddy resolving ocean model covering the North Sea and the Baltic Sea. It is found that assimilation of sea surface temperature data reduces the model root mean square error from 1.13 °C to 0.70 °C. The hybrid scheme is found to reduce the root mean square error slightly more than the simplified KF without OI to 0.66 °C. The inclusion of spatially varying satellite error variances does not improve the performance of the scheme significantly.     

Assessment of sea surface temperature observational networks in the Baltic Sea and North Sea

The satellite and in situ Sea Surface Temperature (SST) observational networks in the Baltic Sea and North Sea are evaluated based on the quality of the gridded SST products generated from the networks. A multi-indicator approach is applied in the assessment. It includes evaluation of data quality, effective data coverage, field reconstruction error and model nowcast error. The results show that the best available full-coverage SST product is generated by assimilating the SST observations to obtain a yearly mean model bias of 0.07 °C and RMSE of 0.64 °C. The effective data coverage rate is 31% by using AVHRR (Advanced Very High Resolution Radiometer) data from NOAA (National Ocean and Atmosphere Administration) satellites 12, 14 and 16. The data redundancy increases rapidly with the number of infrared sensors. Using either NOAA satellite 12 or all 3 satellites makes a small difference with regard to derived effective coverage and the ocean model nowcast error. The influence of using the in situ SST observations in the SST field reconstruction is negligibly small. Instead, the major role of in situ SST observations is in calibrating the satellite observations. To study the relative importance of data quality and data coverage, an assessment is done for two satellite products: one product is based entirely on NOAA 12 data and has larger coverage but lower quality. The other product is a subset of the SAF products (derived from NOAA 14 and 16) and has lower coverage but higher quality. Based on current monitoring, modelling and assimilation technology, the results suggest that the data quality is an important factor in further improving the quality of the gridded SST products. Recommendations are made for possible further improvements of the existing SST observational networks.     

The Gulf Stream, rings and North Atlantic eddy structures from remote sensing (Altimeter and SeaWiFS)

Seasonal SeaWiFS chlorophyll a concentrations cycles and annual changes of altimeter Sea Level Anomaly are derived for the subtropical North Atlantic near  35°N and along a Gulf Stream axis. Spatial structure of SeaWiFS, is defined in terms of deviations from a local seasonal cycle and examined in relation to altimeter eddy structure. In the subtropical region near 35°N, SeaWiFS structure is evident during the spring bloom period with a scale of  430 km, or about twice the eddy scale. A Gulf Stream axis has been selected as a region where the Sea Level Anomaly variance is a maximum. Eddy propagation speeds and scales are examined. Cold-core (cyclonic) rings correspond to areas of high SeaWiFS chlorophyll a. Warm-core (anticyclonic) rings relate to areas of low chlorophyll concentration. Both SeaWiFS structure and eddy structure have a spatial scale of  450 km or twice the ring scale along the Gulf Stream axis. SeaWiFS chlorophyll anomalies and Altimeter Sea Level Anomaly structure have an overall negative correlation coefficient of r = − 0.34. Swirl currents between eddies redistribute surface chlorophyll concentrations and can spatially bias maximum and minimum concentration levels off eddy centre.     

The Mercator global ocean operational analysis system: Assessment and validation of an 11-year reanalysis

This paper presents Prototype Système 2 Global (PSY2G), the first Mercator global Ocean General Circulation Model (OGCM) to assimilate along-track sea level anomaly (SLA) satellite data. Based on a coarse resolution ocean model, this system was developed mainly for climatic purposes and will provide, for example, initial oceanic states for coupled ocean-atmosphere seasonal predictions. It has been operational since 3 September 2003 and produces an analysis and a two-week forecast for the global ocean every week. The PSY2G system uses an incremental assimilation scheme based on the Cooper and Haines [Cooper, M., Haines, K., 1996. Data assimilation with water property conservation. J. Geophys. Res., 101, 1059-1077.] lifting–lowering of isopycnals. The SLA increment is obtained using an optimal interpolation method then the correction is partitioned into baroclinic and barotropic contributions. The baroclinic ocean state correction consists of temperature, salinity and geostrophic velocity increments and the barotropic correction is a barotropic velocity increment. A reanalysis (1993–2003) was carried out that enabled the PSY2G system to perform its first operational cycle. All available SLA data sets (TOPEX/Poséïdon, ERS2, Geosat-Follow-On, Jason1 and Envisat) were assimilated for the 1993–2003 period. The major objective of this study is to assess the reanalysis from both an assimilation and a thermodynamic point of view in order to evaluate its realism, especially in the tropical band which is a key region for climatic studies. Although the system is also able to deliver forecasts, we have mainly focused on analysis. These results are useful because they give an a priori estimation of the qualities and capabilities of the operational ocean analysis system that has been implemented. In particular, the reanalysis identifies some regional biases in sea level variability such as near the Antarctic Circumpolar Current, in the eastern Equatorial Pacific and in the Norwegian Sea (generally less than 1 cm) with a small seasonal cycle. This is attributed to changes in mean circulation and vertical stratification caused by the assimilation methodology. But the model's low resolution, inaccurate physical parameterisations (especially for ocean–ice interactions) and surface atmospheric forcing also contribute to the occurrence of the SLA biases. A detailed analysis of the thermohaline structure of the ocean reveals that the isopycnal lifting–lowering tends to diffuse vertically the main thermocline. The impact on temperature is that the surface layer (0–200 m) becomes cooler whereas in deeper waters (from 500 to 1500 m), the ocean becomes slightly warmer. This is particularly true in the tropics, between 30°N and 30°S. However it can be demonstrated that the assimilation improves the variability in both surface currents and sub-surface temperature in the Equatorial Pacific Ocean.     

Westward moving waves or eddies (Storms) on the Subtropical/Azores Front near 32.5°N? Interpretation of the Eulerian currents and temperature records at moorings 155 (35.5°W) and 156 (34.4°W)

Three Argos buoy-years of Lagrangian data in westward-moving cyclonic eddies, or Storms, near 32.5°N, together with hydrographic measurements, have shown that Storms move westward at nearly 3 km day⁻¹. Water in eddies can be trapped and moved westward by advection within the eddy or by phase propagation of the eddy pattern, so we cannot say that the flow field (or Eulerian mean) is 3 km day⁻¹ westward. Two moorings (155 and 156) deployed in the Storm Corridor have provided a further 8 instrument-years of Eulerian data. The temperature and current records confirmed that two Storms a year passed each mooring over the 2-year measurement period. As expected, there is a lag of 1.3 month at mooring 155 (which is 102 km to the west of mooring 156) with respect to conditions at mooring 156. The progressive vector diagrams (PVDs) derived from the current meter records exhibit fairly regular X (east or zonal) and Y (north or meridional) displacement scales that repeat with semi-annual periodicity (SAP). The dominant current signal is the north component of the SAP, which reaches an amplitude of 18 cm s⁻¹ for the upper layer or near surface current record (242-m depth). The geostrophic north component values derived from altimetry were in good agreement with the upper layer current meter measurements. The large north component amplitude was not interpreted as evidence for Rossby Waves but rather due to the passage of nine eddies (eight complete) of alternate sign (cyclonic, anticyclonic) passing the mooring rigs during the 2-year deployment period. The Y scale shows that the near surface characteristic or mean maximum azimuthal speed is about 35 cm s⁻¹ for cyclonic eddies or Storms, and that this value is reduced to 4 cm s⁻¹ at 1400-m depth. The residual or mean Eulerian currents range from 8 cm s⁻¹ for the upper layer currents to 1 cm s⁻¹ for the deeper currents at 1400-m depth and are predominantly westward. Simple theoretical considerations and idealised numerical simulations show that the mean westward Eulerian current depends markedly on whether the eddy centres pass to the north or south of the rigs. This raises the question as to what do we mean by Eulerian residual currents, even for relatively long records (2 years). It is shown that the strong near surface westward current (6 km day⁻¹) measured at mooring 155 is largely due to a westward-moving eddy field with variable centre offsets. The magnitude of the near surface east–west component of flow was estimated as eastward at 2 cm s⁻¹. The north–south component of mean flow was southward at 2 cm s⁻¹. The deeper records gave a weak westward flow of 1 cm s⁻¹ but did not show a significant southward component for the mean Eulerian flow field. 7.4 float-years of Lagrangian ALACE data in the Subtropical Front region near 740 dbar gave mean east–west flows that were <0.5 cm s⁻¹. Overall, it is shown that the eddy structures propagate westward mainly by phase propagation (i.e. a westward-moving pattern with no westward advection for the current meter to measure), though plane Rossby Wave dynamics appeared inappropriate. Theoretical and modeling considerations show that a speed of 3-km day⁻¹ westward is too large a value for the self-advection of eddies due to the beta effect.     

Recovery of North-East Atlantic temperature fields from profiling floats: Determination of the optimal float number from sampling and instrumental error analysis

Argo is an international project that is deploying an array of temperature and salinity profiling floats over the global ocean. Here we use the error formulation derived from Optimal Statistical Interpolation to estimate statistical errors associated with the recovery of the temperature field in the North-East Atlantic ocean. Results indicate that with the present distribution of floats (119 in the considered domain), scales of wavelength larger than 500 km can be recovered with a relative uncertainty (rms error relative to the standard deviation of the field) of about 7% at 50 m, 8% at 200 m and 10% at 1000 m. This corresponds to mean absolute errors of 0.111 °C at 50 m, 0.104 °C at 200 m and 0.073 °C at 1000 m.The splitting of total errors into instrumental and sampling contributions reveals that, in the present scenario, errors are more due to the small number of floats than to instrumental errors, especially at upper levels. For scales larger than 500 km this will hold true until 200–250 floats are deployed (less than 200 for deep levels). In such a simulated scenario, the number of observations and the technology become approximately equally limiting factors for the accuracy of the temperature field mapping, with total relative errors of less than 2% at upper levels and about 3% at 1000 m.     

Interannual variation of the Polar Front in the Japan/East Sea from summertime hydrography and sea level data

The Polar Front in the Japan/East Sea separates the southern warm water region from the northern cold water region. A merged TOPEX/POSEIDON and ERS-1/2 altimeter dataset and upper water temperature data were used to determine the frontal location and to examine the structure of its interannual variability from 1993 to 2001. The identified frontal location, where sea surface height gradient has a maximum about 10–20 cm over the horizontal distance of 100 km, corresponds well to the maximum subsurface horizontal temperature gradient. The front migrates more widely (36°N–41°N) in the western part of the sea than in the eastern part. The interannual migration induces large variability in upper water temperatures and sea surface height in the western region. Responsible physical mechanisms were studied using a reduced-gravity model. Differences between inflow and outflow change the total volume of warm water, and total warm water volume change in the warm water region uniformly pushes the front in the meridional direction across its mean position in the model simulation. Interannual variation of wind stress causes relatively wide migration of the modeled front in the western part.     

Eddies around Madagascar — The retroflection re-considered

The Agulhas Current with its retroflection and attendant eddy-shedding is the cause of some of the greatest mesoscale variability in the ocean. This paper considers the area to the south and east of Madagascar, which provides some of the source waters of the Agulhas Current, and examines the propagating sea surface height signals in altimetry and output from a numerical model, OCCAM. Both show bands of variability along the axis of the East Madagascar Current (EMC) and along a zonal band near 25°S. Sequences of images plus associated temperature data suggest that a number of westward-propagating eddies are present in this zonal band. The paper then focuses on the region to the south of the island, where ocean colour and infra-red imagery are evocative of an East Madagascar Retroflection. The synthesis of data analysed in this paper, however, shows that remotely observed features in this area can be explained by anticyclonic eddies moving westward through the region, and this explanation is consistent with numerical model output and the trajectories of drifting buoys.     

Assimilation of sea surface temperature in the POL Coastal Ocean Modelling System

A one-dimensional scheme is used to assimilate satellite Sea Surface Temperature data into the Proudman Oceanographic Laboratory Coastal Ocean Modelling System, set up in the Irish Sea with a fine resolution ( 1.8 km). The capabilities of the assimilation scheme are investigated using two different sets of satellite data, of lower and similar resolution to that of the model respectively. Comparison of results with independent data show that assimilation improves the modelled Sea Surface Temperature, but does not address model representation of the temperature vertical structure. It is concluded that for the Irish Sea and at the scales resolved by the model, the assimilation problem cannot be approached in a one-dimensional framework. It is also pointed out that forecast error needs to account explicitly for errors in the representation of the vertical structure of the thermal field.Three-dimensional methods that are suited for coastal systems are then suggested.     

Air–sea gas transfer velocity estimates from the Jason-1 and TOPEX altimeters: Prospects for a long-term global time series

Estimation of global and regional air–sea fluxes of climatically important gases is a key goal of current climate research programs. Gas transfer velocities needed to compute these fluxes can be estimated by combining altimeter-derived mean square slope with an empirical relation between transfer velocity and mean square slope derived from field measurements of gas fluxes and small-scale wave spectra [Frew, N.M., Bock, E.J., Schimpf, U., Hara, T., Hauβecker, H., Edson, J.B., McGillis, W.R., Nelson, R.K., McKenna, S.P., Uz, B.M., Jähne, B., 2004. Air–sea gas transfer: Its dependence on wind stress, small-scale roughness and surface films, J. Geophys. Res., 109, C08S17, doi: 10.1029/2003JC002131.]. We previously reported initial results from a dual-frequency (Ku- and C-band) altimeter algorithm [Glover, D.M., Frew, N.M., McCue, S.J., Bock, E.J., 2002. A Multi-year Time Series of Global Gas Transfer Velocity from the TOPEX Dual Frequency, Normalized Radar Backscatter Algorithm, In: Gas Transfer at Water Surfaces, editors: Donelan, M., Drennan, W., Saltzman, E., and Wanninkhof, R., Geophysical Monograph 127, American Geophysical Union, Washington, DC, 325–331.] for estimating the air–sea gas transfer velocity (k) from the mean square slope of short wind waves (40–100 rad/m) and derived a 6-year time series of global transfer velocities based on TOPEX observations. Since the launch of the follow-on altimeter Jason-1 in December 2001 and commencement of the TOPEX/Jason-1 Tandem Mission, we have extended this time series to 12 years, with improvements to the model parameters used in our algorithm and using the latest corrected data releases. The prospect of deriving multi-year and interdecadal time series of gas transfer velocity from TOPEX, Jason-1 and follow-on altimeter missions depends on precise intercalibration of the normalized backscatter. During the Tandem Mission collinear phase, both satellites followed identical orbits with a mere 73-s time separation. The resulting collocated, near-coincident normalized radar backscatter (σ°) data from both altimeters present a unique opportunity to intercalibrate the two instruments, compare derived fields of transfer velocity and estimate the precision of the algorithm. Initial results suggest that the monthly gas transfer velocity fields generated from the two altimeters are very similar. Comparison of along-track Ku-band and C-band σ° during the collinear phase indicates that observed discrepancies are due primarily to small offsets between TOPEX and Jason-1 σ°. The Jason-1 k values have an apparent bias of + 4% relative to TOPEX, while the precision estimated from the two observation sets is 5–7% and scales with k. The resultant long-term, global, mean k is 16 cm/h.     

Observations of the Fawn Trough Current over the Kerguelen Plateau from instrumented elephant seals

Due to its great meridional extent and relatively shallow depths, the Kerguelen Plateau constitutes a major barrier to the eastward flowing Antarctic Circumpolar Current in the Indian sector of the Southern Ocean. While most of the Antarctic Circumpolar Current transport is deflected north of the Kerguelen Islands, the remainder ( 50 Sv, 1 Sv = 10⁶ m³ s^− 1) must pass south of the islands, most probably through the Fawn and Princess Elizabeth Troughs. However, the paucity of finely resolved quasi-synoptic hydrographic data in this remote and infrequently sampled area has limited the progress in our knowledge of the regional circulation. Since 2004, a new approach using elephant seals from the Kerguelen Islands as autonomous oceanographic profilers has provided new information on the hydrography over the Kerguelen Plateau, covering the entire Antarctic Zone between the Polar Front and Antarctica, with a mean along-track resolution of about 25 km. These finely resolved bio-logged data revealed details of a strong northeastward current found across the Fawn Trough (sill depth: 2600 m; 56°S, 78°E). This so-called Fawn Trough Current transports cold Antarctic waters found mostly south of the Elan Bank, between the Ice Limit (58°S) and the Antarctic Divergence (64°S) in the eastern Enderby Basin, toward the Australian–Antarctic Basin. Our analysis also demonstrates that the Deep Western Boundary Current, which carries cold Antarctic water along the eastern flank of the southern Kerguelen Plateau collides with Fawn Trough Current at the outlet of the Fawn Trough sill. In other words, the Fawn Trough constitutes a veritable bottleneck, channelling the quasi-totality of the Antarctic Circumpolar flow found south of the Polar Front. Thanks to the unprecedented fine resolution of seal-borne data, a branch of flow centered at the Winter Water isotherm of 1 °C is also revealed along the northern escarpment of the Elan Bank, and then along the southern edge of Heard Island. Further analysis of different supplementary data reveals a complex circulation pattern in the entire Enderby Basin, with several distinctive branches of flow being strongly controlled by prominent topographic features such as the Southwest Indian Ridge, Conrad Rise, Elan Bank, and Kerguelen Plateau. This newly emerged frontal structure refines considerably previous large-scale circulation schematics of the area.     

Modelling the hydrodynamics and ecosystem of the North-West European continental shelf for operational oceanography

This paper outlines an approach to complex spatio-temporal marine ecosystem modelling as applied to the North Western European Continental Shelf. The model presented here combines an eddy-permitting (approximately 6 km horizontal resolution) baroclinic model, the Proudman Oceanographic Laboratory Coastal Ocean Modelling System (POLCOMS), with the European Regional Seas Ecosystem Model (ERSEM). This has been run within an operational framework using operationally available high resolution atmospheric and lateral boundary forcing, allowing hindcast and near-real time nowcast simulations to be performed. The modelled surface temperature and chlorophyll distributions are presented, and interannual variations discussed. Validation of both the physical and ecosystem submodels show the system to be effective, whilst highlighting areas where improvements in the system can be made. Distinct regional differences in predictive skill are shown. The system presented is ready for operational implementation to provide products and services for use both scientifically and in coastal zone and shelf seas management activities. A programme of work to update the system is already in place.     

Optimal interpolation of sea surface temperature for the North Sea and Baltic Sea

Sea surface temperature fields of the North Sea and Baltic Sea have been constructed for the year 2001 using a multiplatform Optimal Interpolation scheme. The analyzed fields are constructed every 12 h on a 10 km spatial grid. The product is based upon observations from the three NOAA satellites 12, 14 and 16 together with a large amount of in situ observations. Space dependent covariance functions are estimated from the satellite observations and account for spatial and temporal lags. Several independent methods have been used to assess the error on the sea surface temperature product. Compared against independent in situ observations, the mean RMS difference for the year 2001 is 0.78 °C. The spatial distribution of the errors reveals that the Baltic Sea in general show higher errors than the North Sea. The error statistics throughout the year show a temporal variation of the errors with maximum during summer and winter. Tests with a varying number of satellite observations show that the accuracy of the satellite observations is the most important parameter in terms of reducing the errors on the interpolated sea surface temperature product.     

Comparison between a reanalyzed product by 3-dimensional variational assimilation technique and observations in the Ulleung Basin of the East/Japan Sea

Reanalyzed products from a MOM3-based East Sea Regional Ocean Model with a 3-dimentional variational data assimilation module (DA-ESROM), have been compared with the observed hydrographic and current datasets in the Ulleung Basin (UB) of the East/Japan Sea (EJS). Satellite-borne sea surface temperature and sea surface height data, and in-situ temperature profiles have been assimilated into the DA-ESROM. The performance of the DA-ESROM appears to be efficient enough to be used in an operational ocean forecast system.Comparing with the results from Mitchell et al. [Mitchell, D. A., Watts, D. R., Wimbush, M., Teague, W.J., Tracey, K. L., Book, J. W., Chang, K.-I., Suk, M.-S., Yoon, J.-H., 2005a. Upper circulation patterns in the Ulleung Basin. Deep-Sea Res. II, 52, 1617-1638.], the DA-ESROM fairly well simulates the high variability of the Ulleung Warm Eddy and Dok Cold Eddy as well as the branching of the Tsushima Warm Current in the UB. The overall root-mean-square error between 100 m temperature field reproduced by the DA-ESROM and the observed 100-dbar temperature field is 2.1 °C, and the spatially averaged grid-to-grid correlation between the two temperature fields is high with a mean value of 0.79 for the inter-comparison period.The DA-ESROM reproduces the development of strong southward North Korean Cold Current (NKCC) in summer consistent with the observational results, which is thought to be an improvement of the previous numerical models in the EJS. The reanalyzed products show that the NKCC is about 35 km wide, and flows southward along the Korean coast from spring to summer with maximum monthly mean volume transport of about 0.8 Sv in August–September.     

Review on the seasonal variation of the surface circulation in the Japan/East Sea

The seasonal variation of the surface circulation in the Japan/East Sea (JES) and the Tsushima/Korea Straits (TKS) is reviewed and discussed, focusing on mesoscale and submesoscale variabilities.The monsoon modified by coastal geographical features near Vladivostok generates a dipole of vortex off Vladivostok which induces dramatic changes in the surface circulation in the northwest JES, splitting the Subpolar Gyre into two smaller gyres by generating the Vladivostok Dome. Between these two smaller gyres, the Northwest Thermal Front is formed and current reversal is induced along the North Korean coast. The winter monsoon also induces a current reversal along the Sakhalin coast. The volume transport of the surface Subpolar Gyre has two maxima in January and August. The maximum in August is induced by the summer intensification of the Liman-North Korean Cold Current and the shallow and narrow surface coastal jet generated by the sea ice and snow melting. The maximum in January is induced by the northwest monsoon and associated cooling.Salient features in the TKS are the submesoscale variabilities. In the western channel, submesoscale eddies with length scale of about 80 km and time scale of 5–6 days develop in the cold period. On the lee side of the Tsushima Islands, Karman-like vortex pairs are generated in the warm period. Anticyclonic vortices generated at the northern tip of the Tsushima Islands have a time scale of 5 to 8 days, length scale of 35 to 60 km, and propagate toward the JES with a phase speed of 8 cm/s. Cyclonic vortices south of the anticyclonic counter part of the vortex pairs are rather stationary with intermittent occasional propagation toward the east. The development of stratification seems to be necessary for the development of Karman-like vortex pairs.Summarizing the results above, a schematic surface circulation with seasonal change is proposed.     

Microscale variability in the distributions of large fluorescent particles observed in situ with a planar laser imaging fluorometer

It has been known for decades that particle-size and biomass spectra show regular patterns in the ocean, and that these patterns often show systematic variations with other properties such as total biomass, nutrient concentration, season, and distance (both vertical and horizontal). The recent finding of the ubiquitous nature of layers of phytoplankton < 1 m thick prompted us to explore the fine- and microscale vertical variations of size- and fluorescence-abundance spectra in the ocean. Using a two-dimensional planar laser imaging system mounted on a free-falling platform, we quantified the properties of large fluorescent particles ( 20 μm–2 cm) through the water column, obtaining images every 10–30 cm. These images showed systematic relationships of the spectral properties to total chlorophyll: increased proportions of the smallest particles at high chlorophyll concentrations, and a lengthening of the spectral size range at high total chlorophyll concentrations (more large particles at high chlorophyll concentrations). Further, we observed significant variations of the spectral properties over scales of 1 m and less, and recorded the frequent occurrence of unusual layers of large particles. Our new instrument, which is sensitive to thin layers of enhanced phytoplankton biomass, shows the planktonic community to be highly structured vertically on scales of 1–2 m, particularly within the DCM.     

Air–sea CO2 fluxes in the Caribbean Sea from 2002–2004

Air–sea fluxes in the Caribbean Sea are presented based on measurements of partial pressure of CO₂ in surface seawater, pCO_2sw, from an automated system onboard the cruise ship Explorer of the Seas for 2002 through 2004. The pCO_2sw values are used to develop algorithms of pCO_2sw based on sea surface temperature (SST) and position. The algorithms are applied to assimilated SST data and remotely sensed winds on a 1° by 1° grid to estimate the fluxes on weekly timescales in the region. The positive relationship between pCO_2sw and SST is lower than the isochemical trend suggesting counteracting effects from biological processes. The relationship varies systematically with location with a stronger dependence further south. Furthermore, the southern area shows significantly lower pCO_2sw in the fall compared to the spring at the same SST, which is attributed to differences in salinity. The annual algorithms for the entire region show a slight trend between 2002 and 2004 suggesting an increase of pCO_2sw over time. This is in accord with the increasing pCO_2sw due the invasion of anthropogenic CO₂. The annual fluxes of CO₂ yield a net invasion of CO₂ to the ocean that ranges from − 0.04 to − 1.2 mol m^− 2 year^− 1 over the 3 years. There is a seasonal reversal in the direction of the flux with CO₂ entering into the ocean during the winter and an evasion during the summer. Year-to-year differences in flux are primarily caused by temperature anomalies in the late winter and spring period resulting in changes in invasion during these seasons. An analysis of pCO_2sw before and after hurricane Frances (September 4–6, 2004), and wind records during the storm suggest a large local enhancement of the flux but minimal influence on annual fluxes in the region.     

Effect of wave boundary layer on sea-to-air dimethylsulfide transfer velocity during typhoon passage

A full-spectral third-generation ocean wind–wave model (Wavewatch-III) implemented in the South China Sea is used to investigate the effects of the wave boundary layer on the drag coefficient and the sea-to-air transfer velocity of dimethylsulfide (DMS) during passage of Typhoon Wukong (September 5–11, 2000) with a maximum sustained wind speed of 38 m s^− 1. The model is driven by the reanalyzed surface winds (1° × 1°, four times daily) from the National Centers for Environmental Prediction. It is found that the wave boundary layer evidently enhances (16.5%) the drag coefficient (in turn increases the momentum flux across the air–sea interface), and reduces (13.1%) the sea-to-air DMS transfer velocity (in turn decreases the sea-to-air DMS flux). This indicates the possibility of important roles of wave boundary layer in atmospheric DMS contents and global climate system.     

Atmospheric CO2 measurements and error analysis on seasonal air–sea CO2 fluxes in the Bay of Biscay

Atmospheric molar fraction of CO₂ (xCO₂^atm) measurements obtained on board of ships of opportunity are used to parameterize the seasonal cycle of atmospheric xCO₂ (xCO₂^atm) in three regions of the eastern North Atlantic (Galician and French offshore and Bay of Biscay). Three selection criteria are established to eliminate spurious values and identify xCO₂^atm data representative of atmospheric background values. The filtered data set is fitted to seasonal curve, consisting of an annual trend plus a seasonal cycle. Although the fitted curves are consistent with the seasonal evolution of xCO₂^atm data series from land meteorological stations, only ship-board measurements can report the presence of winter xCO₂^atm minimum on Bay of Biscay. Weekly air–sea CO₂ flux differences (mmol C·m^− 2 day^− 1) produced by the several options of xCO₂^atm usually used (ship-board measurements, data from land meteorological stations and annually averaged values) were calculated in Bay of Biscay throughout 2003. Flux error using fitted seasonal curve relative to on board measurements was minimal, whereas land stations and annual means yielded random (− 0.2 ± 0.3 mmol C·m^− 2·day^− 1) and systematic (− 0.1 ± 0.4 mmol C·m^− 2 day^− 1), respectively. The effect of different available sources of sea level pressure, wind speed and transfer velocity were also evaluated. Wind speed and transfer velocity parameters are found as the most critical choice in the estimate of CO₂ fluxes reaching a flux uncertainty of 7 mmol C·m^− 2·day^− 1 during springtime. The atmospheric pressure shows a notable relative effect during summertime although its influence is quantitatively slight on annual scale (0.3 ± 0.2 mmol C·m^− 2·day^− 1). All results confirms the role of the Bay of Biscay as CO₂ sink for the 2003 with an annual mean CO₂ flux around − 5 ± 5 mmol C m^− 2 day^− 1.