Using a random displacement model to simulate turbulent particle motion in a baroclinic frontal zone: A new implementation scheme and model performance tests

Numerical artifacts can limit accurate simulation of turbulent particle motion when Lagrangian particle-tracking models are implemented in hydrodynamic models with stratified conditions like fronts. Yet, modeling of individual particle motion in frontal regions is critical for understanding sediment dynamics as well as the transport and retention of planktonic organisms. The objective of this research was to develop a numerical technique to accurately simulate turbulent particle motions in a particle-tracking model embedded within a hydrodynamic model of a frontal zone. A new interpolation scheme, the ‘water column profile’ scheme, was developed and used to implement a random displacement model for turbulent particle motions. A new interpolation scheme was necessary because linear interpolation schemes caused artificial aggregation of particles where abrupt changes in vertical diffusivity occurred. The new ‘water column profile’ scheme was used to fit a continuous function (a tension spline) to a smoothed profile of vertical diffusivities at the x–y particle location. The new implementation scheme was checked for artifacts and compared with a standard random walk model using (1) Well Mixed Condition tests, and (2) dye-release experiments. The Well Mixed Condition tests confirmed that the use of the ‘water column profile’ interpolation scheme for implementing the random displacement model significantly reduced numerical artifacts. In dye-release experiments, high concentrations of Eulerian tracer and Lagrangian particles were released at the same location up-estuary of the salt front and tracked for 4 days. After small differences in initial dispersal rates, tracer and particle distributions remained highly correlated (r = 0.84 to 0.99) when a random displacement model was implemented in the particle-tracking model. In contrast, correlation coefficients were substantially lower (r = 0.07 to 0.58) when a random walk model was implemented. In general, model performance tests indicated that the ‘water column interpolation’ scheme was an effective technique for implementing a random displacement model within a hydrodynamic model, and both could be used to accurately simulate diffusion in a highly baroclinic frontal region. The new implementation scheme has the potential to be a useful tool for investigating the influence of hydrodynamic variability on the transport of sediment particles and planktonic organisms in frontal zones.     

Inertial waves in the Ibiza channel

Currentmeter data taken in the Ibiza Channel show the almost permanent presence of near-inertial motions below the mixed layer. They correspond to downwards progressing waves with a vertical group velocity of some m/day. The presence of the Balearic Front sensibly affects the propagation of these inertial waves. Although situations exist in which the passage of atmospheric fronts along the Channel is clearly the generating force of these near-inertial motions we find others in which the energy density in mid-depths is higher than in any other depth. These last situations are closely related to the arrival to the Channel of some of the different water masses which flow around there. A clear correspondence between the presence of relatively strong inertial waves and a noticeable vertical shear of the subinertial flow, evidenced by an averaged Richardson number, is also observed. In some circumstances, the vertical shear of the whole flow (inertial plus subinertial) is higher than the stability limit, that is, Ri < 0.25, favoring the braking of the internal waves. This could be a plausible cause of their decay and a reason to explain why they do not penetrate further than certain depths.     

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.     

Mesoscale subduction at the Almeria–Oran front: Part 1: Ageostrophic flow

This paper presents a detailed diagnostic analysis of hydrographic and current meter data from three, rapidly repeated, fine-scale surveys of the Almeria–Oran front. Instability of the frontal boundary, between surface waters of Atlantic and Mediterranean origin, is shown to provide a mechanism for significant heat transfer from the surface layers to the deep ocean in winter. The data were collected during the second observational phase of the EU funded OMEGA project on RRS Discovery cruise 224 during December 1996. High resolution hydrographic measurements using the towed undulating CTD vehicle, SeaSoar, traced the subduction of Mediterranean Surface Water across the Almeria–Oran front. This subduction is shown to result from a significant baroclinic component to the instability of the frontal jet. The Q-vector formulation of the omega equation is combined with a scale analysis to quantitatively diagnose vertical transport resulting from mesoscale ageostrophic circulation. The analyses are presented and discussed in the presence of satellite and airborne remotely sensed data; which provide the basis for a thorough and novel approach to the determination of observational error.     

A two-layer model of the Rhine plume

The outflow of Rhine water into the shallow Southern Bight of the North Sea leads to almost discontinuous vertical density distributions and sharp frontal structures around the river mouth. Strong tidal motion, wind and baroclinic effects have large influence on the dynamics and dispersion of river water. A three-dimensional tidal model, including advective and diffusive transport of salinity, is used in the two-layer mode for simulation of Rhine water outflow to quantify the interaction of the different processes and the effect on dispersion and mixing of river water. Layer depths are adjusted in a way that no advective transports between upper and lower layer take place in case of sufficiently stable stratification. In case of weak or no stratification the upper layer depth is fixed, and advective transports between layers are computed. Model results show frontal eddy development and (limited) growing internal waves due to baroclinic instability. Comparisons with observational data are presented.     

Distributions of oxygen and carbon stable isotopes and CFC-12 in the water masses of the Southern Ocean at 30°E from South Africa to Antarctica: results of the CIVA1 cruise

This study presents oceanic distributions of stable isotopes (δ¹⁸O of water and δ¹³C of ΣCO₂) and CFC-12 from samples collected during the CIVA1 cruise (February/March 1993), across the Southern Ocean, along a meridian section at 30°E, from South Africa (44°S) to Antarctica (70°S). The isotopic measurements show important variations between the subantarctic surface waters with low δ¹⁸O–high δ¹³C values and the antarctic surface waters with very low δ¹⁸O–low δ¹³C values. The surface distributions of δ¹³C values follow the major frontal oceanic structures; the vertical distribution shows the progressive upwelling from the subantarctic zone to the antarctic divergence of ¹³C-depleted CO₂ derived from remineralization of organic matter. Along the Antarctic continental shelf, between 2500 and 4000 m, a core of water with δ¹⁸O values close to −0.1‰ is associated with a relative maximum in CFC-12 concentration, although this core is not detected by its temperature and salinity parameters. This water mass, which corresponds to recently formed deep water, may originate from the eastward extension of the Weddell gyre or from bottom waters coming from the East and formed near Prydz Bay.     

Mesozooplankton species distribution in the NW and N Iberian shelf during spring 2004: Relationship with frontal structures

We have analysed the mesozooplankton community structure in the southern Bay of Biscay shelf and its relationship with the hydrographic conditions during spring 2004. According to thermohaline characteristics, we observed two frontal zones of distinct origin along the shelf (around 7° and 3°W), that allowed us to differentiate three different hydrographic domains. The westernmost part of the shelf (WC), defined by the presence of relatively warm and salty water related to the presence of the Iberian Poleward Current (IPC), the easternmost region (EC), characterized by colder and fresher water and subject to the influence of freshwater inputs from the Adour river in the French coast, and a region in the Central Cantabrian Sea (CC), where thermohaline characteristics were intermediate between these two extremes. The mixing layer depth (MLD) regime in these areas was also different: the WC region was characterized by a mixed water column, whereas in the EC region the river discharges produces stratification of the upper meters of the water column (< 10 m); in the CC region, we found a distinct vertical mixing regime that separated coastal (stratification) from shelf (mixed water column) stations, giving rise to a notorious across-shelf front. We found a good match between the aforesaid hydrographic regions and the distribution of mesozooplankton species composition and community assemblages: the Mantel correlation between physical variables and mesozooplankton distribution was highly significant (n = 63, r = 0.70, α < 0.001). In the WC region, the community was dominated by Paracalanus parvus, Oithona helgolandica, Acartia clausi and Clausocalanus pergens, while in the EC region the most dominant species were Noctiluca scintillans, Oncaea media and Temora longicornis. The CC region showed similar composition of copepods than the WC region, but larvaceans (Oikopleura spp. and Fritillaria spp.) were more abundant in the CC region than in the WC region. Within each zone, the relative abundances of the dominant species differed between coastal and shelf locations.     

About process times in some deepwater formation events in the Mediterranean Sea

Experimental data gathered on some Mediterranean deepwater formation sites have been analysed in order to check the Chapman (Ch), Visbeck and Maxworthy scale laws, their final mixing water depth and process times estimates. All the analysed sites have a high probability to be “internally constrained” so that the baroclinic Rossby radius is the dominant horizontal length scale and the final chimney depth is independent of the rate of rotation as indicated by Visbeck and Maxworthy criteria; the Eady instability times are generally very fast (3/f) but the least final mixing process times are too long compared with the meteorological forcing ones. It is attempted here to identify some new criteria to give some insight on the processes time scale, which is still an open problem. Conclusions about the Rhodes gyre appear questionable, because the experimental data appear inconsistent with any theoretical criterion.     

The vernal bloom in heterogeneous convection: a numerical study of Baltic restratification

The onset of the vernal phytoplankton bloom in heterogeneous convection at a temperature below the temperature of maximum density (θ_ρ) is studied by means of an idealized numerical model. The modeled setting resembles the well-documented conditions in the Baltic Sea. For the low sea surface temperature found in early spring, the convection in this inland sea is driven by surface heating, in marked contrast to oceanic conditions, where cooling serves the purpose. The heterogeneity of the physical environment is enhanced by a lateral freshwater flux, which stratifies parts of the domain and allows the surface to warm above θ_ρ in the stratified areas. Consequently, the same surface heat flux may cause convection in some parts and stabilisation in other parts of the domain. The development of a layer of anomalously cold water with θ<θ_ρ is demonstrated numerically.The physical heterogeneity is strongly reflected in the development of the vernal bloom, with significant phytoplankton growth possible only in the stratified areas. The non-local organised structures in the mixing layer play a dominant role in the development of the bloom. An analysis of the model results shows the presence of biological intrusions during the bloom, highly reminiscent of their atmospheric counterparts.     

One-dimensional modelling of the plankton ecosystem of the north-western Corsican coastal area in relation to meteorological constraints

In order to study the influence of wind mixing on the spring variability of the plankton production of the north western Corsican coastal area, a one-dimensional (1D), vertical, coupled hydrodynamic/biological model (ECOHYDROMV) is used. A hydrodynamic 1D model of the water column with a k–l turbulent closure is applied. The biological model comprises six state variables, representing the plankton ecosystem in the spring period: phytoplankton, copepods, nitrate, ammonium, particulate organic matter of phytoplanktonic origin and particulate organic matter of zooplanktonic origin. The system is influenced by turbulence (expressed by the vertical eddy diffusivity), temperature and irradiance. The model takes into account momentum and heat surface fluxes computed from meteorological data in order to simulate a typical spring atmospheric forcing for the considered area. Results show that primary production vertical structure is characterised by a subsurface maximum which deepens with time and is regulated by the opposite gradients of nitrate concentration and irradiance. Surface plankton productivity is mainly controlled by turbulent vertical transport of nutrients into the mixed layer. The short time scale variability of turbulent mixing generated by the wind appears to be responsible for the plurimodal shape of plankton blooms, observed in the considered area. Furthermore, the model is applied to the study of the spring evolution of the plankton communities off the bay of Calvi (Corsica) for the years 1986 and 1988. In order to initiate and validate the model, time series of hydrological, chemical and biological data have been used. The model reproduces accurately the spring evolution of the phytoplankton biomass measured in situ and illustrates that its strong variability in those years was in close relation to the variability of the wind intensity.     

The structure of dissipation in the western Irish Sea front

We report on an intensive campaign in the summer of 2006 to observe turbulent energy dissipation in the vicinity of a tidal mixing front which separates well mixed and seasonally stratified regimes in the western Irish Sea. The rate of turbulent dissipation ε was observed on a section across the front by a combination of vertical profiles with the FLY dissipation profiler and horizontal profiles by shear sensors mounted on an AUV (Autosub). Mean flow conditions and stratification were obtained from a bed mounted ADCP and a vertical chain of thermistors on a mooring. During an Autosub mission of 60 h, the vehicle, moving at a speed of ~ 1.2 m s^− 1, completed 10 useable frontal crossings between end points which were allowed to move with the mean flow. The results were combined with parallel measurements of the vertical profile of ε which were made using FLY for periods of up to 13 h at positions along the Autosub track. The two data sets, which show a satisfactory degree of consistency, were combined to elucidate the space–time variation of dissipation in the frontal zone. Using harmonic analysis, the spatial structure of dissipation was separated from the strong time dependent signal at the M₄ tidal frequency to yield a picture of the cross-frontal distribution of energy dissipation. A complementary picture of the frontal velocity field was obtained from a moored ADCP and estimates of the mean velocity derived from the thermal wind using the observed density distribution. which indicated the presence of a strong (0.2 m s^− 1) jet-like flow in the high gradient region of the front. Under neap tidal conditions, mean dissipation varied across the section by 3 orders of magnitude exceeding 10^− 2 W m^− 3 near the seabed in the mixed regime and decreasing to 10^− 5 W m^− 3. in the strongly stratified interior regime. The spatial pattern of dissipation is consistent in general form with the predictions of models of tidal mixing and does not reflect any strong influence by the frontal jet.     

On the background level of Tc, Sr and Cs in the North Atlantic

Mean concentrations of the anthropogenic radioactive oceanographic tracers ⁹⁹Tc, ⁹⁰Sr and ¹³⁷Cs have been measured as 0.005, 1.6 and 2.5 Bq m⁻³ in oceanic Northeast Atlantic surface water, east and northeast of the Azores, in 1992. This is, apparently, the first published value for fallout “background” ⁹⁹Tc in oceanic Atlantic water.Comparison with older data indicates an observed half life for ⁹⁰Sr and ¹³⁷Cs in the northeast Atlantic surface water of 20 yr corresponding to a mean residence time of 80–100 yr for the stable elements.The observed ⁹⁹Tc/⁹⁰Sr ratio (3 × 10⁻³) in the Azores samples is 10 times higher than the theoretical fission yield decay corrected to 1992. This is in agreement with published data on rain water samples and may be characteristic for 1960's global fallout. Furthermore, the measured ¹³⁷Cs/⁹⁰Sr ratio is not significantly different from that observed for global fallout. There do not appear to be any additional significant sources of artificial radionuclides in this region.     

Variability of mesozooplankton spatial distribution in the North Aegean Sea, as influenced by the Black Sea waters outflow

The North Aegean Sea constitutes an important region of the Mediterranean Sea since in its eastern part the mesotrophic, low salinity and relatively cold water from the Black Sea (outflowing from the Dardanelles strait) meets the oligotrophic, warm and very saline water of Levantine origin, thus forming a thermohaline front. Mesozooplankton samples were collected at discrete layers according to the hydrology of the upper 100 m, during May 1997 and September 1998. In May highest biomass and abundance values (up to 66.82 mg m^− 3 and 14,157 ind m^− 3) were detected in the 10–20 m layer (within the halocline) of the stations positioned close to the Dardanelles strait. The front moved slightly southwards in September, characterized by high biomass and abundance values within the halocline layer. The areas moderately or non influenced by Black Sea water revealed lower standing stock values than the frontal area in both cruises and maxima were detected in the uppermost low salinity layer. Samples collected at the stations and/or layers more influenced by Black Sea water were distinguished from those collected at layers and/or stations more affected by Levantine waters in both periods. In May the former samples were characterized by the copepods Acartia clausi, Centropages typicus, Paracalanus parvus. The abundance of the above species decreased gradually with increasing salinity, in the horizontal and/or in the vertical dimension, with a parallel increase of the copepods Oithona plumifera, Oithona copepodites, Oncaea media, Ctenocalanus vanus, Farranula rostrata. During September the frontal area as well as that covered by the modified Black Sea water, were highly dominated by the cladoceran Penilia avirostris and doliolids. For both seasons, MDS plots, issued from the combination of mesozooplankton and water-type data, revealed the gradual differentiation of zooplankton composition from the frontal area towards the area covered by Levantine water, following the spreading and mixing of the Black sea water. The observed temporal and spatial variability in the distribution pattern of mesozooplankton standing stock and species composition seems to depend considerably on the variability of circulation and frontal flows.     

An inverse model for calculation of global volume transport from wind and hydrographic data

The P-vector inverse method has been successfully used to invert the absolute velocity from hydrographic data for the extra-equatorial hemispheres, but not for the equatorial region since it is based on the geostrophic balance. A smooth interpolation scheme across the equator is developed in this study to bring together the two already known solutions (P-vectors) for the extra-equatorial hemispheres. This model contains four major components: (a) the P-vector inverse model to obtain the solutions for the two extra-equatorial hemispheres, (b) the objective method to determine the Ψ-values at individual islands, (c) the Poisson equation-solver to obtain the Π-values over the equatorial region from the volume transport vorticity equation, and (d) the Poisson equation-solver to obtain the Ψ and depth-integrated velocity field (U, V) over the globe from the Poisson Ψ-equation. The Poisson equation-solver is similar to the box model developed by Wunsch. Thus, this method combines the strength from both box and P-vector models. The calculated depth-integrated velocity and Ψ-field agree well with earlier studies.     

Seasonal variability of fractionated phytoplankton, biomass and primary production in the Straits of Magellan

The results on the distribution of phytoplankton biomass (expressed as Chla) and primary production (¹⁴C assimilation), during three oceanographic cruises carried out during Austral spring and at the end of the summer and the autumn in the Straits of Magellan, suggest a strong variability of trophic levels for this ecosystem.Seasonal evolution of the biomass concentration goes from the spring maximum of 2.33 μg/l through a sharp decrease, 0.49 μg/l, observed at the end of summer, until the minimum of 0.24 μg/l measured during the autumn.The trophic conditions are dependent on hydrographic, meteo-climatic and geo-morphological characteristics: at the Atlantic entrance and between the two Angosturas the strong mixing of water column limit the development of phytoplankton; at the Western opening and along the Pacific arm the complex exchange mechanisms with the ocean, the glacio-fluvial contribution and the presence of a thermohaline front near the Isla Carlos III influence both biomass and primary production distributions. The maximum values are reached in the Central Zone (Paso Ancho) characterized by high stability of the water column.Primary production ranged from a minimum of 12.3 to a maximum of 125.9 mgC m⁻² h⁻¹. The overall trend seems to be a progressive and simultaneous increase from the Pacific and Atlantic openings to the Central Zone of Paso Ancho where the maximum value was reached. In general, biomass and primary production distributions correspond quite well except for the area of Isla Carlos III where biological and chemico-physical causes tend to limit ¹⁴C assimilation.Contribution of pico-phytoplankton (< 2 μm) to total biomass appears to be time dependent: in the blooms observed during spring a very modest incidence (< 6%) was observed whereas became more (> 50%) during the summer-autumn seasons when total biomass was decreasing.Within the Straits, at the end of summer, the contribution of pico-phytoplankton primary production is 59%, whereas nano and microplankton contribute 39% and 2%, respectively. At the oceanic external stations the photosynthetic activity of the bigger size-fraction (> 2 μm) is predominant (> 50%).These findings support the hypothesis that the pico-phytoplankton ( < 2 μm) is substantially constant, whereas temporal variations are due to the larger (> 10 μm) cells only.     

Experimental investigation on the drag reducing efficiency of the outer-layer vertical blades

An experimental assessment has been made of the drag reducing efficiency of the outer-layer vertical blades, which were first devised by Hutchins and Choi (Proceedings of ASME FEDSM’02 2002 ASME Fluids Engineering Division Summer Meeting, Montreal) and Hutchins (An investigation of larger-scale coherent structures in fully developed turbulent boundary layers, Hutchins N (2003), PhD thesis, University of Nottingham). The reported drag reduction efficiency, which was as much as 30%, was quantified only in terms of the reduction in the local skin-friction coefficient. The assessment of the drag reducing efficiency did not take the side effects of the inclusion of the blades into considerations. Those effects are the increase in the wetted surface area and the flow disturbances due to the presence of the blades. In the present study, a series of drag force measurements in towing tank has been performed toward the assessments of the total drag reduction efficiency of the outer-layer vertical blades. It was found that for the case of h 4.0 × z 4.0 (h/δ = 1.04), the outer-layer vertical blades array achieved about 9.6% drag reduction without considering the increase in the wetted surface area. A proper scaling method to give collapsed plot of drag reduction efficiency C _F/C _F0 was attempted, but the correlation remained limited. Of the two scaling methods, the outer scaling is found to be relevant one.     

A novel method for estimating vertical eddy diffusivities using diurnal signals with application to western Long Island Sound

We present an approach that allows the estimation of vertical eddy diffusivity coefficients from buoy measurements made at two or more depths. By measuring the attenuation and phase lag of a scalar signal generated periodically at the surface as it propagates downwards, the vertical eddy diffusivity coefficients can be calculated as K_v = ωΔz²/2ln²(α₂/α₁), where α₂/α₁ is the ratio of the real amplitudes at frequency ω at the two depths separated by Δz = z₂ − z₁; as K_V = ωΔz²/2φ², where φ is the phase lag at the frequency ω; or as K_v = iωΔz²/ln²(U₂/U₁), where U₂/U₁ is the ratio of the complex signal amplitudes at the two depths. The method requires that horizontal fluxes be small at the ω frequency and that the signal-to-noise ratios at the two depths allow the determination of the amplitude and phase of ω.Application of this method to summertime 2004 western Long Island Sound oxygen and temperature buoy measurements at two depths provides a time-series of two-day average vertical eddy diffusivity estimates. Using these eddy diffusivities in conjunction with measured vertical concentration gradients, we obtain a time-series of vertical transport rates for oxygen and heat and estimate mean downward fluxes for June and July as 150–260 mMol m^− 2 day^− 1 and 100–400 W m^− 2 respectively. These estimates are of a similar magnitude to sub-pycnocline O₂ and heat demands of 240 ± 200 mMol m^− 2 day^− 1 and 180 ± 60 W m^− 2 that we infer from simple budgets, implying that vertical transport is significant to both budgets.The eddy coefficients obtained from the independent O₂ and temperature measurements have a 68% correlation, and the O₂ flux estimates show a correlation of 41% to measured rates of change in bottom dissolved oxygen levels. Our results indicate that extended time-series of eddy diffusivity coefficients can be obtained from in situ buoy measurements and the method shows promise as a way to constrain the vertical transport variability in budgets of dissolved materials in estuaries.     

A 3D, finite element model for baroclinic circulation on the Vancouver Island continental shelf

This paper describes the development and application of a 3-dimensional model of the barotropic and baroclinic circulation on the continental shelf west of Vancouver Island, Canada. A previous study with a 2D barotropic model and field data revealed that several tidal constituents have a significant baroclinic component (the K1 in particular). Thus we embarked on another study with a 3D model to study the baroclinic effects on the residual and several selected tidal constituents.The 3D model uses a harmonic expansion in time and a finite element discretization in space. All nonlinear terms are retained, including quadratic bottom stress, advection and wave transport (continuity nonlinearity). The equations are solved as a global and a local problem, where the global problem is the solution of the wave equation formulation of the shallow water equations, and the local problem is the solution of the momentum equation for the vertical velocity profile. These equations are coupled to the advection-diffusion equation for density so that density gradient forcing is included in the momentum equations. However, the study presented here describes diagnostic calculations for the baroclinic residual circulation only.The model is sufficiently efficient that it encourages sensitivity testing with a large number of model runs. In this sense, the model is akin to an extension of analytical solutions to the domain of irregular geometry and bottom topography where this parameter space can be explored in some detail.In particular, the consequences of the sigma coordinate system used by the model are explored. Test cases using an idealized representation of the continental shelf, shelf break and shelf slope, lead to an estimation of the velocity errors caused by interpolation errors inherent in the sigma coordinate system. On the basis of these estimates, the computational grid used in the 2D model is found to have inadequate resolution. Thus a new grid is generated with increased accuracy in the region of the shelf break. However, even with increased resolution, spurious baroclinic circulation seaward of the shelf break and in the vicinity of Juan de Fuca canyon remained a significant problem when the pressure gradient terms were evaluated using the σ coordinate system and using a realistic density profile.With the new grid, diagnostic calculations of the barotropic and baroclinic residual circulation are performed using forcing from the observed σ_t (density) field and from the gradient of this field.