Mode-1 internal tides in the western North Atlantic Ocean

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

The impact of the M2 internal tide on data-assimilative model estimates of the surface tide

We use a synthetic data experiment to assess the accuracy of ocean tides estimated from satellite altimetry data, with emphasis on the impact of the phase-locked internal tide, which has a surface expression of several centimeters near its sites of genesis. Previous tidal estimates have regarded this signal as a random measurement error; however, it is deterministic and not scale-separated from the barotropic (surface) tide around complex bathymetric features. The synthetic data experiments show that the internal tide has a negligible impact on the barotropic tidal fields inferred under these circumstances, and the barotropic dissipation (a quadratic functional of the tidal fields) is in good agreement with the energetics of the three-dimensional regional primitive equations model which is the source of the synthetic data.     

Tidal data inversion: interpolation and inference

Initial efforts in applying inverse methods to studies of ocean tides have focused on making the best use of a small number of observations to map tidal fields in a large area. As such, inversion can be viewed as an objective analysis scheme which uses a dynamically appropriate spatial covariance, derived from the shallow water equations, to interpolate and smooth a sparse data set. Data from recent altimetry missions are not sparsely distributed relative to tidal wavelengths in the open ocean, apparently reducing the need for complicated dynamically based interpolation schemes. Altimetric data sets are also quite large, making application of rigorous inversion methods to global tidal modeling a challenging computational problem. We describe here a new iterative solution scheme which allows us for the first time to fit the full set of TOPEX/Poseidon cross-over differences. The resulting solution (TPXO.3) fits validation tide gauges significantly better than previous inverse solutions. TPXO.3 also reduces residual cross-over variances relative to other recent inverse and empirical solutions, particularly in shallow water where improvements are dramatic. With the new solution approach very significant improvements in global tidal models should be possible in shallow areas and in the vicinity of complex bathymetry, where high-accuracy tidal modeling remains a challenging problem. With the recent improvements in the definition of tidal elevations in the open ocean it should now also be possible to resolve some long unanswered questions about tidal energetics and dynamics. Inverse methods provide a natural framework for addressing these issues, and making inferences about tidal dynamics. In particular, by bringing data and dynamics together in a single solution, we can rigorously test the consistency of the two. We present results of global and local inversions which suggest that over elongated bathymetric features oriented perpendicular to tidal flows, energy dissipation in the open ocean is significantly enhanced, presumably due to conversion of barotropic tidal motions into baroclinic modes. For the M₂ tide our preliminary results suggest that perhaps as much as 0.5 TW of energy is dissipated in this manner. However, due to the simplified linear dynamics and limited spatial resolution used for our inversion, there are significant uncertainties associated with these results. A more careful application of inverse methods to make more rigorous inferences about tidal energetics, including use of more reasonable prior dynamics, and the highest possible spatial resolution, should allow for closure of the tidal energy budget within the next few years.     

Performance comparison between vertical and horizontal arrays for geoacoustic inversion

Most of the research on model-based geoacoustic inversion techniques has concentrated on data collected using moored vertical receiver arrays. However, there are many advantages to considering geoacoustic inversion using a towed horizontal array. Towed arrays are easily deployed from a moving platform; this mobility makes them well suited for surveying large areas for sea-bed properties. Further, if a model-based geoacoustic inversion scheme uses both a towed source and array, the separation between the two can be kept short, which reduces the requirement for range-dependent modeling. Range-independent modeling is used for inverting all the horizontal array data considered in this paper. Using the Inversion Techniques Workshop Benchmark Test Cases, the performance of a horizontal (simulated towed) and vertical arrays are compared and found to be very similar. However, it will be shown that, for Benchmark Test Case 3, where the bathymetry is flat and a hidden bottom intrusion exists, a towed horizontal array is ideal for determining the range-dependent sea-bed properties. The practical advantages of using a towed array are clear and the purpose of this paper is to show that the performance is similar (and in some cases better) than using moored vertical arrays.     

Geoacoustic Information Content of Horizontal Line Array Data

This paper examines the effectiveness of horizontal line arrays (HLAs) for matched-field inversion (MFI) by quantifying geoacoustic information content for a variety of experiment and array factors, including array length and number of sensors, source range and bearing, source-frequency content, and signal-to-noise ratio (SNR). Emphasis is on bottom-moored arrays, while towed arrays are also considered, and a comparison with vertical line array (VLA) performance is made. The geoacoustic information content is quantified in terms of marginal posterior probability distributions (PPDs) for model parameters estimated using a fast Gibbs sampler approach to Bayesian inversion. This produces an absolute, quantitative estimate of the geoacoustic parameter uncertainties which can be directly compared for various experiment and array factors.     

Inverse reconstructions of ecosystem flows in investigating regime shifts: impact of the choice of objective function

Inverse analysis is increasingly used in ecosystem modelling to objectively reconstruct a large number of unknown flows or interactions from a small number of observations. This type of analysis may be useful in relating observed regime shifts in ecosystem structure to underlying processes. Inversions of ecosystem flow networks currently use a constrained least-squares solution which at the same time minimizes the squared norm (the sum of squares) of the reconstructed flows. This minimum norm (MN) inversion is thought to be a parsimonious solution to the ecosystem flow inverse problem, but it may well not reflect how ecosystems are organised. It has been proposed instead that ecosystems evolve to maximize energy/mass flows or that they maximize the information content of the network weighted by ecosystem flows (ascendancy). We used simulated inverse experiments, where inverse analyses are applied to simulations of flow networks, to explore objective functions different than the MN generally used. We could not compute inverse solutions that maximize ascendancy because the objective function is unbounded. We could calculate inversions that maximize flows; however, these generally overestimated the simulated flows, even though the simulations were designed to maximize flows. It appears that the ecosystem flow inverse problem is too under-determined (too few data relative to the number of unknowns) to allow the use of these maximizing goal functions. We introduce a new minimization that simultaneously minimizes the squared flows and the squared differences between flows. This smoothing minimization makes the inverse flows as even as possible and it helps with some technical issues with MN inversions. The simulated inverse experiments indicated that this smoothed norm (SM) is the most robust in comparative analyses of contrasting ecosystem states, such as those that can be associated with regime shifts. Like the MN inversion, the SM inversion has no ecological basis. However, it is a conservative norm that is less likely to produce false differences between the dynamics of regimes.     

The influence of lateral mixing on a phytoplankton bloom: Distribution in the Kerguelen Plateau region

A unique phytoplankton bloom appears every year during the austral spring/summer in the northern Kerguelen Plateau region. The Kerguelen Ocean and Plateau compared Study (KEOPS) showed that an increase in subsurface iron coming up from the seafloor through vertical mixing was responsible for the observed increase in chlorophyll-a above the plateau. We demonstrate that the bloom pattern is not a simple increase of biomass over shallow water: it is strongly influenced by the bathymetry and its spatial extent controlled by strong currents around the plateau. Here we focus on lateral mixing process to explain the particular shape of the bloom. We use the Smagorinsky [1963. General circulation experiments with the primitive equations. I. The basic experiment. Monthly Weather Review 91 (3), 99–164] formula to estimate and map fields of lateral mixing time scales (τ) due to barotropic tidal currents, barotropic atmospheric forced currents, Ekman velocities and geostrophic velocities. Results show that short time scale mixing is strongly influenced by the tides while the other processes have minor influences. Comparisons of τ and satellite chlorophyll-a images show that the spatial pattern of the bloom seems to be delimited by a barrier of high lateral mixing that is essentially due to tides. This emphasises the role played by the tides over the Kerguelen Plateau in supplying iron to the phytoplankton and containing the horizontal shape of the bloom. This is one of the first times such a link has been demonstrated, which has implications for the study of iron advection in the ocean.     

A barotropic model of the linear semidiurnal tide over a continental slope

Semidiurnal tides, and especially the lunar tide M₂, are dominant dynamics in the Bay of Biscay. Strong tidal currents are associated with the presence of a significant continental slope. By combining Newton's gravitation laws and Euler's equations, Laplace's equations contain the astronomical forcing responsible for the observed semidiurnal tides. In shallow waters, this direct forcing is often neglected. We study here its influence on the tidal dynamics over the continental slope through the development of a simple model describing the barotropic semidiurnal dynamics on a transect perpendicular to the slope. This new model results from the combination of two different models, i.e. the one developed by Rosenfeld and Beardsley (1987), which takes into account the tide-generating force, and that of Battisti and Clarke (1982), which neglects it. A first model is developed by neglecting the direct astronomical forcing in equations: it consists in solving a second-order homogeneous propagation equation for the barotropic semidiurnal tide and needs only coastal conditions as well as the knowledge of the along-slope wave number of the solution. For a mean slope typical of the South Brittany area, this non-forced model provides results in accordance with those of Battisti and Clarke and Le Cann (1990): in particular, in the upper part of the slope, it shows a polarization inversion of tidal ellipses characteristic of the tidal dynamics observed in this area. Then, the direct astronomical forcing is kept in equations. The simple model developed without this forcing is fitted in order to solve the resulting forced propagation equation for the barotropic tide. The solution of this second model is the sum of a forced wave responding to the direct astronomical forcing and of a free wave generated at the coastal boundary. Under the same boundary conditions, the results obtained with the influence of the tide-generating force are then compared with those obtained without it. This comparison allows one to apprehend the importance of the direct astronomical forcing on tidal dynamics across the slope: in particular, the main difference appears in deep waters where this forcing induces a phase-lag between the plain and the shelf for the sea-surface slope.     

Non-linear interaction between tidal and subinertial barotropic flows in the Strait of Gibraltar

The present work establishes clear relationships of the amplitude and phase variation of the barotropic M₂ signal in the velocity of the current with the barotropic subinertial flow in the Strait of Gibraltar. The analytical procedure is applied on data from Gibraltar Experiment in order to obtain barotropic subinertial series and the amplitude and phase variation of the M₂ signal series involves harmonic analysis, empirical orthogonal function analysis and complex demodulation. In addition, cross-spectral analysis has been applied to study these relationships, concerning which non-linear interaction between M₂ and the subinertial oscillation is proposed as the responsible physical mechanism. An analytical solution characterizing this type of non-linear interaction is offered in explanation of the experimental results.     

Tidal observations on the West Coast,South Island,New Zealand

Harmonic tidal constants, calculated from sea surface elevation observations at Jackson Bay on the West Coast of the South Island, are consistent with available semi‐diurnal and diurnal tidal phase distributions. Current observations taken over a 111 day period at mid‐depth in 1505 m of water on the southern flank of the Challenger Plateau and over a 240 day period in 1430 m of water on the South Island western coast continental slope, are subject to tidal analysis. At both sites there is a component of energy flux directed across the isobaths and only at the northern site for the M2 tide is the phase consistent with a dominant progressive barotropic tide. The successive 30 day harmonic constants at the southern continental slope site exhibit a trend in the M2 tidal ellipse speed and ellipticity suggesting the presence of a regular internal tide. Superposition of ‘internal tidal’ and barotropic tidal flows, as prescribed from progressive‐ and standing‐wave elevations, to fit the observations indicates that the ‘internal tide’ is probably associated with the first baroclinic mode. At the current‐meter depths the speeds of the ‘internal tide’ for the M2 tide are about the same as the barotropic speeds, whereas, the S2 ‘internal tide’ speeds are larger than those of the barotropic tide. The consistency of the trend in ellipse parameters lends support to the theoretical progressive trapped barotropic tidal flows being a good approximation to the actual barotropic tide. Some support for the hypothesis that the S2 tide on the West Coast of New Zealand has a substantial standing wave contribution is given by the northern observations, where the ratio of the S2: M2 internal tidal ellipse current amplitudes are substantially larger than the ratio of the elevations, the internal tide being generated by across‐isobath flows.     

Internal tide generation by seamounts

The generation of internal tidal wave fields by barotropic tidal flow past a representative seamount is computed by modelling the seamount as a pillbox, and linearising the equations for internal wave dynamics. This is justifiable for mid-ocean seamounts, which constitute steep topography for internal waves of tidal frequency. For linearly polarised barotropic tidal flow, the resulting flow field consists of conical beams radiating from the region above the seamount, with largest velocities aligned with the barotropic flow. These beams vary with azimuthal angle but resemble the corresponding beams from two-dimensional steep topography, particularly in the barotropic flow direction. They are primarily forced by the barotropic flow over the seamount, which is amplified by the topography and is independent of the stratification if the radius of the seamount is sufficiently large. In a barotropic tidal flow of 1 cm/s amplitude, energy fluxes from individual seamounts are of order 10⁶ W. Summing this over all seamounts higher than 1 km gives baroclinic energy generation of order 5.10⁹ W, a number that is less than estimates of baroclinic energy flux from the continental slopes and the Hawaiian ridge, but is comparable with them.     

Modeling the tide and wind-induced circulation in Buzzards Bay

Hydrodynamic model application to Buzzards Bay is performed using a three-dimensional Boundary-fitted Hydrodynamic model in this study. The model is forced with observed tidal harmonic constants along the open boundaries and winds on the surface. The main focus of the present study is to model the detailed wind and tide-induced circulation in Buzzards Bay. The observed surface elevations and currents given in [Butman, B., Signell, R., Shoukimas, P., Beardsley, R.C., 1988. Current Observations in Buzzards Bay, 1982–1986. Open File Report 88-5. United States Geological Survey] and the tide and current harmonics given in [Signell, R.P., 1987. Tide- and Wind-forced Currents in Buzzards Bay, Massachusetts. Technical Report WH-87-15. Woods Hole Oceanographic Institution, Woods Hole, Massachusetts] are used to validate the model predictions. The calibrated model is then used to study the relative contributions of tidal and wind forcing on the instantaneous and residual circulation in Buzzards Bay. The amplitudes and phases of the principal tidal constituents at 10 tidal stations in Buzzards Bay obtained from a harmonic analysis of a 60-day simulation compare well with the observed data. The predicted amplitude and phase of the M₂ tidal constituent of surface elevations at these stations are, respectively, within 4 cm and 5° of the observed data. The errors in the model-predicted M₂ harmonic principal current speeds are less than 6 cm/s, and the principal current directions and phases are within 14° of the observations. The observed surface elevations and currents given in [Butman, B., Signell, R., Shoukimas, P., Beardsley, R.C., 1988. Current Observations in Buzzards Bay, 1982–1986. Open File Report 88-5. United States Geological Survey] are used to validate the model-predicted low-frequency surface elevations and currents. The model predictions in low-frequency surface elevations at Woods Hole closely follow the trends seen in the observations with a correlation coefficient of 0.735, but fail to capture some of the peak surges seen in the observations. The model-predicted low-frequency currents in the east–west direction at stations in Buzzards Bay compare well with the observations with the correlation coefficient exceeding 0.811 and the model capturing the trends seen in the observations, for the most part. However, the model-predicted north–south velocities does not compare well with the observations. The model-predictions agree with the observations that the tidal currents in Vineyard Sound lagged the currents in Buzzards Bay by more than 3 h. The interaction of wind stress with large bathymetric gradients was shown to cause many vortices in Buzzards Bay, as seen from the model predictions. Model simulations show that the winds play a more dominant role than the tides in the generation of the barotropic residual currents in Buzzards Bay, while the model-predicted tide-induced residual current was seen to be small.     

Centimeter-Level Positioning of Seafloor Acoustic Transponders from a Deeply-Towed Interrogator

An array of three seafloor transponders was acoustically surveyed to centimeter precision with a deeply-towed interrogator. Measurements of two-way acoustic travel time and hydrostatic pressure made as the interrogator was towed above the array were combined in a least-squares adjustment to estimate the interrogator and transponder positions in two surveys spanning two years. No transponder displacements were expected at this site in the interior of the Juan de Fuca Plate (48^?11′ N, 127^?12′ W) due to the lack of active faults. This was confirmed to a precision of ±2 cm by least-squares adjustment. Marginally detectable blunders in the observations were shown to affect the transponder position estimates by no more than 3 mm, demonstrating the geometric strength of the data set. The accumulation of many hundreds of observations resulted in a significant computational burden on the least-squares inversion procedure. The sparseness of the normal matrix was exploited to reduce by a factor of 1000 the number of calculations. The acoustic survey results suggested that the near-bottom sound speed fields during the two surveys were in better agreement than inferred from yearly single-profile conductivity, temperature, and pressure (CTD) measurements.     

A three-dimensional numerical study of river plume mixing processes in Otsuchi Bay,Japan

The three-dimensional numerical model SUNTANS is applied to investigate river plume mixing in Otsuchi Bay, an estuary located along the Sanriku Coast of Iwate, Japan. Results from numerical simulations with different idealized forcing scenarios (barotropic tide, baroclinic tide, and diurnal wind) are compared with field observations to diagnose dominant mixing mechanisms. Under the influence of combined barotropic, baroclinic and wind forcing, the model reproduces observed salinity profiles well and achieves a skill score of 0.94. In addition, the model forced by baroclinic internal tides reproduces observed cold-water intrusions in the bay, and barotropic tidal forcing reproduces observed salt wedge dynamics near the river mouths. Near these river mouths, vertically sheared flows are generated due to the interaction of river discharge and tidal elevations. River plume mixing is quantified using vertical salt flux and reveals that mixing near the vicinity of the river mouth, is primarily generated by the barotropic tidal forcing. A 10 ms^?1 strong diurnal breeze compared to a 5 ms^?1 weak breeze generates higher mixing in the bay. In contrast to the barotropic forcing, internal tidal (baroclinic) effects are the dominant mixing mechanisms away from the river mouths, particularly in the middle of the bay, where a narrow channel strengthens the flow speed. The mixing structure is horizontally asymmetric, with the middle and northern parts exhibiting stronger mixing than the southern part of the bay. This study identifies several mixing hot-spots within the bay and is of great importance for the coastal aquaculture system.     

Influence of the White Sea on tides in adjacent marginal seas of the North European Basin

An investigation into the interaction of surface M ₂ tides in the system of marginal seas of the North European Basin is carried out using the three-dimensional finite-element hydrostatic model QUODDY-4. Three numerical experiments are performed for this purpose. In the first (control), the model equations are solved in the system of the Norwegian, Greenland, Barents, and White seas; thereby the interaction of the tides in these seas is explicitly taken into account. In the second experiment, the White Sea is excluded from consideration and the no-flux condition is posed at the entrance to the sea. The third experiment uses an approach in which the observed tidal elevations that determine the existence of a finite horizontal transport of barotropic energy to the White Sea are specified at the open boundary of the White Sea. It is shown that changes in tidal dynamics represented by changes in the amplitudes and phases of tidal elevations and in the barotropic tidal velocity ellipse parameters are within the model noise in experiments 2 and 3 when compared with the control experiment. On the contrary, changes in energy characteristics (the horizontal wave transport, density, and dissipation rate of barotropic tidal energy) are equal to or greater (in order of magnitude) than the energy characteristics themselves.     

Particle-Velocity-Field Difference Smoothing for Coherent Source Localization in Spatially Nonuniform Noise

This communication considers the problem of estimating 2-D directions of arrival (DOAs) of multiple coherent signals under spatially nonuniform noise (spatially inhomogeneous temporary white noise) using an array of vector hydrophones. A novel preprocessing method called particle-velocity-field difference smoothing (PVFDS) is proposed. The key idea underlying the PVFDS is to remove the spatially nonuniform noise by using the matrix difference of pairs of particle-velocity data correlation matrices, and to decorrelate the coherent signals by summing these difference correlation matrices. Unlike most of other existing preprocessing techniques, such as spatial smoothing and forward–backward averaging, the PVFDS processing does not decrease the array aperture. For arbitrary array geometries, the PVFDS can resolve up to four coherent signals, and for centro–symmetric arrays, forward–backward averaging can double this number to eight. Monte Carlo simulations illustrate that the PVFDS-based eigenstructure algorithms can offer better performance than the particle-velocity-field smoothing (PVFS)-based counterparts.      

Quantifying data information content in geoacoustic inversion

This paper examines the information content in matched-field geoacoustic inverse problems as a function of a variety of experiment factors, with the aim of guiding data collection and processing to achieve the best possible inversion results. The information content of the unknown geoacoustic parameters is quantified in terms of their marginal posterior probability distributions, which define the accuracy expected in inversion. Marginal distributions are estimated using a fast Gibbs sampler approach to Bayesian inversion, which provides an efficient, unbiased sampling of the multi-dimensional posterior probability density. When sampled to convergence, the marginal distributions are found to have simple, smooth shapes that facilitate straightforward comparisons. The approach is general; the specific examples considered here include factors such as the number of sensors in the receiving array, array length, source-receiver range, source frequency, number of frequencies, source bandwidth, and signal-to-noise ratio     

Mass, Heat and Salt Balances in the Eastern Barents Sea Obtained by Inversion of Hydrographic Section Data

Standard hydrological section data, collected in the eastern Barents Sea in September 1997, have been analyzed using a variational data assimilation technique. This method allows us to obtain temperature, salinity and velocity fields that are consistent with observations and dynamically balanced within the framework of a steady-state model describing large-scale nearly geostrophic circulation. Error bars of the optimized fields are computed by explicit inversion of the Hessian matrix. The optimized velocity field is in agreement with independent velocity observations derived from surface drifter trajectories in the southwestern part of the Barents Sea. Optimized fields provide the following estimates of integral characteristics of the circulation in the region: i) the North Cape current transport is 2.12 ± 0.25 Sv; ii) the Karskie Vorota Strait throughflow is 0.7 ± 0.06 Sv; iii) heat flux with Atlantic water is 4.7 ± 0.16⋅10¹¹ W; iv) salt import from the Atlantic Ocean is 7.41 ± 0.46⋅10³ kg/s. The imbalance of the heat budget in the eastern part of the Barents Sea indicates the presence of statistically insignificant surface heat fluxes which are less than 1 W/m². This revised version was published online in July 2006 with corrections to the Cover Date.     

Wave model predictions in the Black Sea: Sensitivity to wind fields

This paper evaluates the impact of using different wind field products on the performance of the third generation wave model SWAN in the Black Sea and its capability for predicting both normal and extreme wave conditions during 1996. Wind data were obtained from NCEP CFSR, NASA MERRA, JRA-25, ECMWF Operational, ECMWF ERA40, and ECMWF ERA-Interim. Wave data were obtained in 1996 at three locations in the Black Sea within the NATO TU-WAVES project. The quality of wind fields was assessed by comparing them with satellite data. These wind data were used as forcing fields for the generation of wind waves. Time series of predicted significant wave height (H_mo), mean wave period (T_m02), and mean wave direction (DIR) were compared with observations at three offshore buoys in the Black Sea and its performance was quantified in terms of statistical parameters. In addition, wave model performance in terms of significant wave height was also assessed by comparing them against satellite data.The main scope of this work is the impact of the different available wind field products on the wave hindcast performance. In addition, the sensitivity of wave model forecasts due to variations in spatial and temporal resolutions of the wind field products was investigated. Finally, the impact of using various wind field products on predicting extreme wave events was analyzed by focussing on storm peaks and on an individual storm event in October 1996. The numerical results revealed that the CFSR winds are more suitable in comparison with the others for modelling both normal and extreme events in the Black Sea. The results also show that wave model output is critically sensitive to the choice of the wind field product, such that the quality of the wind fields is reflected in the quality of the wave predictions. A finer wind spatial resolution leads to an improvement of the wave model predictions, while a finer temporal resolution in the wind fields generally does not significantly improve agreement between observed and simulated wave data.     

Using satellite altimetry to correct mean temperature and salinity fields derived from Argo floats in the ocean regions around Australia

We present results from a suite of methods using in situ temperature and salinity data, and satellite altimetric observations to obtain an enhanced set of mean fields of temperature, salinity (down to 2000-m depth) and steric height (0/2000 m) for a time-specific period (1992–2007). Firstly, the improved global sampling resulting from the introduction of the Argo program, enables a representative determination of the large-scale mean oceanic structure. However, shortcomings in the coverage remain. High variability western boundary current eddy fields, continental slope and shelf boundaries may all be below their optimal sampling requirements. We describe a simple method to supplement and improve standard spatial interpolation schemes and apply them to the available data within the waters surrounding Australia (100°E–180°W; 50°S–10°N). This region includes a major current system, the East Australian Current (EAC), complex topography, unique boundary currents such as the Leeuwin Current, and large ENSO related interannual variability in the southwest Pacific. We use satellite altimetry sea level anomalies (SLA) to directly correct sampling errors in in situ derived mean surface steric height and subsurface temperature and salinity fields. The surface correction is projected through the water column (using an empirical model) to modify the mean subsurface temperature and salinity fields. The errors inherent in all these calculations are examined. The spatial distribution of the barotropic–baroclinic balance is obtained for the region and a ‘baroclinic factor’ to convert the altimetry SLA into an equivalent in situ height is determined. The mean fields in the EAC region are compared with independent estimates on repeated XBT sections, a mooring array and full-depth CTD transects.