Barotropic wind-driven circulation patterns in a closed rectangular basin of variable depth influenced by a peninsula or an island

We study how a coastal obstruction (peninsula or coastal island) affects the three-dimensional barotropic currents in an oblong rectangular basin with variable bathymetry across the basin width. The transverse depth profile is asymmetric and the peninsula or island lies in the middle of the long side of the rectangle. A semi-spectral model for the Boussinesq-approximated shallow water equations, developed in Haidvogel et al. and altered for semi-implicit numerical integration in time in Wang and Hutter, is used to find the steady barotropic state circulation pattern to external winds. The structural (qualitative) rearrangements and quanti2tative features of the current pattern are studied under four principal wind directions and different lengths of the peninsula and its inclination relative to the shore. The essentially non-linear relationships of the water flux between the two sub-basins (formed by the obstructing peninsula) and the corresponding cross-sectional area left open are found and analysed. It is further analysed whether the depth-integrated model, usually adopted by others, is meaningful when applied to the water exchange problems. The flow through the channel narrowing is quantitatively estimated and compared with the three-dimensional results. The dynamics of the vortex structure and the identification of the up-welling/down-welling zones around the obstruction are discussed in detail. The influence of the transformation of the peninsula into a coastal island on the global basin circulation is considered as are the currents in the channel. The geometric and physical reasons for the anisotropy of the current structure which prevail through all obtained solutions are also discussed.Dedicated to Professor L.A. Mysak on the occasion of his sixtieth birthday.     

Internal waves and temperature fronts on slopes

Time series measurements from an array of temperature miniloggers in a line at constant depth along the sloping boundary of a lake are used to describe the ‘internal surf zone” where internal waves interact with the sloping boundary. More small positive temperature time derivatives are recorded than negative, but there are more large negative values than positive, giving the overall distribution of temperature time derivatives a small negative skewness. This is consistent with the internal wave dynamics; fronts form during the up-slope phase of the motion, bringing cold water up the slope, and the return flow may become unstable, leading to small advecting billows and weak warm fronts. The data are analysed to detect ‘events’, periods in which the temperature derivatives exceed a set threshold. The speed and distance travelled by ‘events’ are described. The motion along the slope may be a consequence of (a) instabilities advected by the flow (b) internal waves propagating along-slope or (c) internal waves approaching the slope from oblique directions. The propagation of several of the observed ‘events’ can only be explained by (c), evidence that the internal surf zone has some, but possibly not all, the characteristics of the conventional ‘surface wave’ surf zone, with waves steepening as they approach the slope at oblique angles.     

Combining T/P altimetric data with hydrographic data to estimate the mean dynamic topography of the North Atlantic and improve the geoid

The mean dynamic topography of the surface of the North Atlantic is estimated using an inverse model of the ocean circulation constrained by hydro-graphic and altimetric observations. In the North Atlantic, altimetric observations have no significant impact on the topography estimate because of the limited precision of available geoid height models. They have a significant impact, however, when uncertainties in the density field are increased to simulate interpolation errors in regions where hydrographic data are scarce. This result, which moderates the conclusion drawn by Ganachaud and co-workers of no significant contribution of altimetric observations to the determination of the large-scale steady circulation, reflects the simple idea that altimetric data are most useful near the surface of the ocean and in areas where the hydrography is poorly determined. One application of the present inverse estimate of the mean dynamic topography is to compute a geoid height correction over the North Atlantic which reduces the uncertainty in the geoid height expanded to spherical harmonic 40 down to a level of about 5 cm.     

Wind-induced baroclinic response of Lake Constance

We present results of various circulation scenarios for the wind-induced three-dimensional currents in Lake Constance, obtained with the aid of a semi-spectral semi-implicit finite difference code developed in Haidvogel et al. and Wang and Hutter. Internal Kelvin and Poincaré-type oscillations are demonstrated in the numerical results, whose periods depend upon the stratification and the geometry of the basin and agree well with measured data. By solving the eigenvalue problem of the linearized shallow water equations in the two-layered stratified Lake Constance, the interpretation of the oscillations as Kelvin and Poincaré-type waves is corroborated.     

Reflection of equatorial Kelvin waves at eastern ocean boundaries Part I: hypothetical boundaries

A baroclinic shallow-water model is developed to investigate the effect of the orientation of the eastern ocean boundary on the behavior of equatorial Kelvin waves. The model is formulated in a spherical polar coordinate system and includes dissipation and non-linear terms, effects which have not been previously included in analytical approaches to the problem. Both equatorial and middle latitude response are considered given the large latitudinal extent used in the model. Baroclinic equatorial Kelvin waves of intraseasonal, seasonal and annual periods are introduced into the domain as pulses of finite width. Their subsequent reflection, transmission and dissipation are investigated. It is found that dissipation is very important for the transmission of wave energy along the boundary and for reflections from the boundary. The dissipation was found to be dependent not only on the presence of the coastal Kelvin waves in the domain, but also on the period of these coastal waves. In particular the dissipation increases with wave period. It is also shown that the equatorial β-plane approximation can allow an anomalous generation of Rossby waves at higher latitudes. Nonlinearities generally have a small effect on the solutions, within the confines of this model.     

Water mass distribution in Fram Strait and over the Yermak Plateau in summer 1997

The water mass distribution in northern Fram Strait and over the Yermak Plateau in summer 1997 is described using CTD data from two cruises in the area. The West Spitsbergen Current was found to split, one part recirculated towards the west, while the other part, on entering the Arctic Ocean separated into two branches. The main inflow of Atlantic Water followed the Svalbard continental slope eastward, while a second, narrower, branch stayed west and north of the Yermak Plateau. The water column above the southeastern flank of the Yermak Plateau was distinctly colder and less saline than the two inflow branches. Immediately west of the outer inflow branch comparatively high temperatures in the Atlantic Layer suggested that a part of the extraordinarily warm Atlantic Water, observed in the boundary current in the Eurasian Basin in the early 1990s, was now returning, within the Eurasian Basin, toward Fram Strait. The upper layer west of the Yermak Plateau was cold, deep and comparably saline, similar to what has recently been observed in the interior Eurasian Basin. Closer to the Greenland continental slope the salinity of the upper layer became much lower, and the temperature maximum of the Atlantic Layer was occasionally below 0.5 °C, indicating water masses mainly derived from the Canadian Basin. This implies that the warm pulse of Atlantic Water had not yet made a complete circuit around the Arctic Ocean. The Atlantic Water of the West Spitsbergen Current recirculating within the strait did not extend as far towards Greenland as in the 1980s, leaving a broader passage for waters from the Atlantic and intermediate layers, exiting the Arctic Ocean. A possible interpretation is that the circulation pattern alternates between a strong recirculation of the West Spitsbergen Current in the strait, and a larger exchange of Atlantic Water between the Nordic Seas and the inner parts of the Arctic Ocean.     

A dynamically consistent analysis of circulation and transports in the southwestern Weddell Sea

An inverse model is applied for the analysis of hydrographic and current meter data collected on the repeat WOCE section SR4 in the Weddell Sea in 1989–1992. The section crosses the Weddell Sea cyclonic gyre from Kapp Norvegia to the northern end of the Antarctic Peninsula. The concepts of geostrophy, conservation of planetary vorticity and hydrostatics are combined with advective balances of active and passive properties to provide a dynamically consistent circulation pattern. Our variational assimilation scheme allows the calculation of three-dimensional velocities in the section plane. Current speeds are small except along the coasts where they reach up to 12 cm/s. We diagnose a gyre transport of 34 Sverdrup which is associated with a poleward heat transport of 28 × 10¹² W corresponding to an average heat flux of 15 Wm^–2 in the Weddell Sea south of the transect. This exceeds the estimated local flux on the transect of 2 Wm^–2. As the transect is located mostly in the open ocean, we conclude that the shelf areas contribute significantly to the ocean-atmosphere exchange and are consequently key areas for the contribution of the Weddell Sea to global ocean ventilation. Conversion of water masses occuring south of the section transform 6.6 ± 1.1 Sv of the inflowing warm deep water into approximately equal amounts of Weddell Sea deep water and Weddell Sea bottom water. The volume transport of surface water equals in the in-and outflow. This means that almost all newly formed surface water is involved in the deep and bottom water formation. Comparison with the results obtained by pure velocity interpolation combined with a hydrographic data subset indicates major differences in the derived salt transports and the water mass conversion of the surface water. The differences can be explained by deviations in the structure of the upper ocean currents to which shelf areas contribute significantly. Additionally a rigorous variance analysis is performed. When only hydrographic data are used for the inversion both the gyre transport and the poleward heat transport are substantially lower. They amount to less than 40% of our best estimate while the standard deviations of both quantities are 6.5 Sv and 37 × 10¹² W, respectively. With the help of long-term current meter measurements these errors can be reduced to 2 Sv and 8 × 10¹² W. Our result underlines the importance of velocity data or equivalent information that helps to estimate the absolute velocities.     

Turbulent viscosity optimized by data assimilation

As an alternative approach to classical turbulence modelling using a first or second order closure, the data assimilation method of optimal control is applied to estimate a time and space-dependent turbulent viscosity in a three-dimensional oceanic circulation model. The optimal control method, described for a 3-D primitive equation model, involves the minimization of a cost function that quantifies the discrepancies between the simulations and the observations. An iterative algorithm is obtained via the adjoint model resolution. In a first experiment, a k ± L model is used to simulate the one-dimensional development of inertial oscillations resulting from a wind stress at the sea surface and with the presence of a halocline. These results are used as synthetic observations to be assimilated. The turbulent viscosity is then recovered without the k + L closure, even with sparse and noisy observations. The problems of controllability and of the dimensions of the control are then discussed. A second experiment consists of a two-dimensional schematic simulation. A 2-D turbulent viscosity field is estimated from data on the initial and final states of a coastal upwelling event.     

Reflection of equatorial Kelvin waves at eastern ocean boundaries Part II: Pacific and Atlantic Oceans

The effect of viscosity, non linearities, incident wave period and realistic eastern coastline geometry on energy fluxes are investigated using a shallow water model with a spatial resolution of 1/4 degree in both meridional and zonal directions. Equatorial and mid-latitude responses are considered. It is found that (1) the influence of the coastline geometry and the incident wave period is more important for the westward energy flux than for the poleward flux, and (2) the effect of the inclination of the eastern ocean boundary on the poleward energy flux, for the Pacific and Atlantic Oceans, decline as the period of the incident wave increases. Furthermore, the model simulations suggest that the poleward energy fluxes from meridional boundaries give plausible results for motions of seasonal and annual periods. For comparatively shorter periods, a realistic coastline geometry has to be included for more accurate results. It is recommended that any numerical model involving the reflection of baroclinic Rossby waves (of intraseasonal, seasonal or annual periods) on the eastern Pacific or Atlantic Oceans, should consider the effect of the coastline geometry in order to improve the accuracy of the results.     

The evolution of an internal bore at the Malin shelf break

Observations of internal waves were made at the Malin shelf edge during SESAME (Shelf Edge Studies Acoustic Measurement Experiment), a part of the NERC LOIS-SES experiment, in August-September 1996. These measurements provide a high resolution dataset demonstrating internal wave generation and propagation. This note presents observations of the evolution of an internal bore. The process is shown clearly in a sequence of thermistor chain tows across the shelf break covering a complete tidal cycle, as the double-sided bore transforms into a group of undulations and eventually into more distinct solitary waveforms. Current structures associated with the bore and waves were also observed by ship-mounted ADCP. Analysis of the waveforms in terms of the linear modes and empirical orthogonal functions (EOFs) indicate the dominance of the first mode, which is typical of a shallow water seasonal thermocline environment. Determination of the phase speed of the waves from the consecutive ship surveys enabled the Doppler shift in the towed data to be removed, allowing analysis of the real length scales of the waves. The bore evolution has been modelled using a first order non-linear KdV model for the first mode, initialised with the waveform in the first survey. Comparison of the model and the observations show close agreement in the amplitudes, length scales, phase speeds and separations of the leading internal waves as they evolve. Finally, analysis of the observed internal wave shapes indicates that, within the uncertainties of measurement, the wave-lengths lie between those predicted by first and second order soliton theory.     

Climatic features of the Mediterranean Sea detected by the analysis of the longwave radiative bulk formulae

Some important climatic features of the Mediterranean Sea stand out from an analysis of the systematic discrepancies between direct measurements of longwave radiation budget and predictions obtained by the most widely used bulk formulae. In particular, under clear-sky conditions the results show that the surface values of both air temperature and humidity over the Mediterranean Sea are larger than those expected over an open ocean with the same amount of net longwave radiation. Furthermore, the twofold climatic regime of the Mediterranean region strongly affects the downwelling clear-sky radiation. This study suggests that a single bulk formula with constant numerical coefficients is unable to reproduce the fluxes at the surface for all the seasons.     

Temperature anomalies in high northerly latitudes and their link with the El Niño/Southern Oscillation

I report the discovery of a low frequency temperature oscillation in the eastern North Atlantic (NA), which was significantly correlated with the Southern Oscillation Index (SOI) in the tropical Pacific, but led the latter index by a number of months. This discovery is significant, because it demonstrates a link between the tropical Pacific and the high northerly latitudes which cannot readily be explained in terms of El Niño/Southern Oscillation (ENSO) feedbacks from the tropics, and opens up the possibility that ENSO and temperature anomalies in northerly climes, may actually have a common origin within, or even external to, the global climate system.     

Drifter observations of the Hebrides slope current and nearby circulation patterns

The mean flow at and around the Hebrides and Shetland Shelf slope is measured with ARGOS tracked drifters. Forty-two drifters drogued at 50 m were deployed in three circles over the Hebrides slope at 56.15°N in two releases, one on 5th December, 1995 and the second on 5–9th May, 1996. The circles span a distance of some 20 km from water depths of 200 m to 1200 m. Drifters are initially advected poleward along-slope by the Hebrides slope current at between 0.05 and 0.70 m s^–1 in a laterally constrained (25–50 km wide) jet-like flow. Drifters released in winter remained in the slope current for over 2000 km whilst summer drifters were lost from the slope current beyond the Wyville-Thomson Ridge, a major topographic feature at 60°N. Dispersion from the slope region into deeper waters occurs at bathymetric irregularities, particularly at the Anton Dohrn Seamount close to which the slope current is found to bifurcate, both in summer and winter, and at the Wyville-Thomson Ridge where drifters move into the Faeroe Shetland Channel. Dispersion onto the continental shelf occurs sporadically along the Hebrides slope. The initial dispersion around the Hebrides slope is remarkably sensitive to initial position, most of the drifters released in shallower water moving onto the shelf, whilst those in 1000 m or more are mostly carried away from the slope into deeper water near the Anton Dohrn Seamount. The dispersion coefficients estimated in directions parallel and normal to the local direction of the 500 m contour, approximately the position of the slope current core, are approximately 8.8 × 10³ m² s^–1 and 0.36 × 10³ m² s^–1, respectively, during winter, and 11.4 × 10³ m² s^–1 and 0.36 x 10³ m² s^–1, respectively, during summer. At the slope there is a minimum in across-slope mean velocity, Reynolds stress, and across-slope eddy correlations. The mean across-slope velocity associated with mass flux is about 4 × 10^–3 m s^–1 shelfward across the shelf break during winter and 2 × 10^–3 m s^–1 during summer. The drifters also sampled local patterns of circulation, and indicate that the source of water for the seasonal Fair Isle and East Shetland currents are the same, and drawn from Atlantic overflows at the Hebrides shelf.     

Intercomparison of oceanic and atmospheric forced and coupled mesoscale simulations. Part I: Surface fluxes

A mesoscale non-hydrostatic atmospheric model has been coupled with a mesoscale oceanic model. The case study is a four-day simulation of a strong storm event observed during the SEMAPHORE experiment over a 500 × 500 km² domain. This domain encompasses a thermohaline front associated with the Azores current. In order to analyze the effect of mesoscale coupling, three simulations are compared: the first one with the atmospheric model forced by realistic sea surface temperature analyses; the second one with the ocean model forced by atmospheric fields, derived from weather forecast re-analyses; the third one with the models being coupled. For these three simulations the surface fluxes were computed with the same bulk parametrization. All three simulations succeed well in representing the main oceanic or atmospheric features observed during the storm. Comparison of surface fields with in situ observations reveals that the winds of the fine mesh atmospheric model are more realistic than those of the weather forecast re-analyses. The low-level winds simulated with the atmospheric model in the forced and coupled simulations are appreciably stronger than the re-analyzed winds. They also generate stronger fluxes. The coupled simulation has the strongest surface heat fluxes: the difference in the net heat budget with the oceanic forced simulation reaches on average 50 Wm⁻² over the simulation period. Sea surface-temperature cooling is too weak in both simulations, but is improved in the coupled run and matches better the cooling observed with drifters. The spatial distributions of sea surface-temperature cooling and surface fluxes are strongly inhomogeneous over the simulation domain. The amplitude of the flux variation is maximum in the coupled run. Moreover the weak correlation between the cooling and heat flux patterns indicates that the surface fluxes are not responsible for the whole cooling and suggests that the response of the ocean mixed layer to the atmosphere is highly non-local and enhanced in the coupled simulation.     

A three-dimensional model study of warm core ring interaction with continental shelf and slope

The interaction of warm core rings with a western wall and shelf/slope is examined with a three-dimensional primitive-equation model. The model ring is initialized with an axisymmetric Gaussian-type anticyclonic eddy placing far from the coastal boundary to allow the ring to freely propagate towards the wall and shelf/slope. The ring initially propagates steadily to the southwest at about 3 km/day under the combined planetary β and nonlinear effects. When colliding with a wall, the ring adjusts into a ‘D’ shape and moves poleward under primarily the image effect. When colliding with a shelf and slope, the ring however becomes stalled and bounces on and off the shelf/slope with little net movement. Small cyclones marked by strong upwelling are generated near the shelfbreak. Cyclones and anticyclones also are spawned at the periphery of the ring. Satellite SST images and concurrent ADCP transects are used to illustrate the strong interaction of a Gulf Stream warm ring (99B) with the Middle Atlantic Bight.     

Changes in sediment and organic carbon accumulation in a highly-disturbed ecosystem: The Sacramento-San Joaquin River Delta (California, USA)

We used the Sacramento-San Joaquin River Delta CA (Delta, hereafter) as a model system for understanding how human activities influence the delivery of sediment and total organic carbon (TOC) over the past 50–60 years. Sediment cores were collected from sites within the Delta representing the Sacramento River (SAC), the San Joaquin River (SJR), and Franks Tract (FT), a flooded agricultural tract. A variety of anthropogenic tracers including ¹³⁷Cs, total DDE (∑DDE) and brominated diphenyl ether (BDE) congeners were used to quantify sediment accumulation rates. This information was combined with total organic carbon (TOC) profiles to quantify rates of TOC accumulation. Across the three sites, sediment and TOC accumulation rates were four to eight-fold higher prior to 1972. Changes in sediment and TOC accumulation were coincident with completion of several large reservoirs and increased agriculture and urbanization in the Delta watershed. Radiocarbon content of TOC indicated that much of the carbon delivered to the Delta is “pre-aged” reflecting processing in the Delta watershed or during transport to the sites rather than an input of predominantly contemporary carbon (e.g., 900–1400 years BP in surface sediments and 2200 yrs BP and 3610 yrs BP at the base of the SJR and FT cores, respectively). Together, these data suggest that human activities have altered the amount and age of TOC accumulating in the Delta since the 1940s.     

Stratospheric background aerosol and polar cloud observations by laser backscattersonde within the framework of the European project “Stratospheric Regular Sounding”

The Stratospheric Regular Sounding project was planned to measure regularly the vertical profiles of several tracers like ozone, water vapor, NO_x, ClO_x and BrO_x radicals, aerosol, pressure and temperature, at three latitudes, to discriminate between the transport and photochemical terms which control their distribution. As part of this project, the “Istituto di Fisica dell’Atmosfera” launched nine laser backscattersondes (LABS) on board stratospheric balloons to make observations of background aerosol and PSCs. LABS was launched with an optical particle counter operated by the University of Wyoming. Observations have been performed in the arctic, mid-latitudes and tropical regions in different seasons. Polar stratospheric clouds have been observed in areas inside and outside the polar vortex edge. A background aerosol was observed both in mid-latitudes and in arctic regions with a backscattering ratio of 1.2 at 692 nm. Very stratified aerosol layers, possibly transported into the lower stratosphere by deep convective systems, have been observed in the lower stratosphere between 20 and 29 km in the tropics in the Southern Hemisphere.     

Stationary magnetospheric convection on November 24, 1981. 1. A case study of “pressure gradient/minimum-B” auroral arc generation

We present two case studies in the night and evening sides of the auroral oval, based on plasma and field measurements made at low altitudes by the AUREOL-3 satellite, during a long period of stationary magnetospheric convection (SMC) on November 24, 1981. The basic feature of both oval crossings was an evident double oval pattern, including (1) a weak arc-type structure at the equatorial edge of the oval/polar edge of the diffuse auroral band, collocated with an upward field-aligned current (FAC) sheet of ≈1.0 μA m⁻², (2) an intermediate region of weaker precipitation within the oval, (3) a more intense auroral band at the polar oval boundary, and (4) polar diffuse auroral zone near the polar cap boundary. These measurements are compared with the published magnetospheric data during this SMC period, accumulated by Yahnin et al. and Sergeev et al., including a semi-empirical radial magnetic field profile B_Z in the near-Earth neutral sheet, with a minimum at about 10–14 R_E. Such a radial B_Z profile appears to be very similar to that assumed in the “minimum B/cross-tail line current” model by Galperin et al. (GVZ92) as the “root of the arc”, or the arc generic region. This model considers a FAC generator mechanism by Grad-Vasyliunas-Boström-Tverskoy operating in the region of a narrow magnetic field minimum in the near-Earth neutral sheet, together with the concept of ion non-adiabatic scattering in the “wall region”. The generated upward FAC branch of the double sheet current structure feeds the steady auroral arc/inverted-V at the equatorial border of the oval. When the semi-empirical B_Z profile is introduced in the GVZ92 model, a good agreement is found between the modelled current and the measured characteristics of the FACs associated with the equatorial arc. Thus the main predictions of the GVZ92 model concerning the “minimum-B” region are consistent with these data, while some small-scale features are not reproduced. Implications of the GVZ92 model are discussed, particularly concerning the necessary conditions for a substorm onset that were not fulfilled during the SMC period.     

Basaltic volcanism in the northern sector of the main Ethiopian rift

The different basalt types related to rift structure development have been investigated, starting from the pre-rift stage in the northern Ethiopian rift and its eastern escarpment and plateau.The basic volcanic rocks are represented mainly by transitional basalts, both in the pre-rift (plateau) and rift (escarpment and rift floor) stages. A striking feature is that although the plateau basalts show clear tholeiitic affinity and the rift basalts reveal a somewhat pronounced “alkaline” character, the REE and LILE element abundances, however, progressively decrease from the “tholeiitic” basalts of the plateau to the “alkaline” basalts of the rift.All data support the view that such contrasting features may be attributed to a continuous depletion of hygromagmatophile (REE, LILE) elements in the mantle source material, related to the large volumes of magmas produced in the early phase of rift structure development. The transition from “tholeiitic” (plateau) to “alkaline” (rift) transitional basalts is related to decreasing intensity of extensional movements.