ALBERT

Beschreibung: The authors derive a string function that describes the propagation of large-scale, potentially large amplitude, baroclinic energy anomalies in a two-layer ocean with variable topography and rotation parameter. The generality of the two-layer results allows results for the 1-layer, 1.5-layer, inverted 1.5-layer, lens, and dome models to be produced as limiting-cases. The string function is a scalar field that acts as a streamfunction for the propagation velocity. In the linear case the string function is simply c2o/f, where co is the background baroclinic shallow water wave speed, and typically describes propagation poleward on the eastern boundaries, westward (with some topographic steering) over the middle ocean, and equatorward on the western boundaries. In the more general nonlinear case, the string function is locally distorted by the anomaly. In the fully nonlinear examples of a lens or dome, there is no rest or background string function; the string function is generated entirely by the disturbance and propagation is due to asymmetric distribution of the anomalous mass over the string function contours. It is shown that conventional beta/topographic propagation results (e.g., beta drift of eddies, the Nof speed of cold domes) can be obtained as limiting cases of the string function. The string function provides, however, more general propagation velocities that are also usually simpler to derive. The first baroclinic mode string function for the global oceans is calculated from hydrographic data. The westward propagation speeds in the ocean basins as derived from the meridional gradient of the string function are typically two to five times faster than those expected from standard theory and agree well with the propagation speeds observed for long baroclinic Rossby waves in the TOPEX/Poseidon data.

Beschreibung: Empirical orthogonal function (EOF) analyses (rotated or not) are widely used in climate research. In recent years there have been several studies in which EOF analyses were used to highlight potential physical mechanisms associated with climate variability. For example, several SST modes were identified such as the “Tropical Atlantic Dipole,” the “Tropical Indian Ocean Dipole,” and different SLP modes in the Northern Hemisphere winter. In this note it is emphasized that caution should be used when trying to interpret these statistically derived modes and their significance. Indeed, from a synthetic example it is shown that patterns derived from EOF analyses can be misleading at times and associated with very little climate physics.

Beschreibung: Comparisons are made between a time series of meteorological surface layer observational data taken on board the R/V Knorr, and model analysis data from the European Centre for Medium-Range Weather Forecasting (ECMWF) and the National Centers for Environmental Prediction (NCEP). The observational data were gathered during a winter cruise of the R/V Knorr, from 6 February to 13 March 1997, as part of the Labrador Sea Deep Convection Experiment. The surface layer observations generally compare well with both model representations of the wintertime atmosphere. The biases that exist are mainly related to discrepancies in the sea surface temperature or the relative humidity of the analyses. The surface layer observations are used to generate bulk estimates of the surface momentum flux, and the surface sensible and latent heat fluxes. These are then compared with the model-generated turbulent surface fluxes. The ECMWF surface sensible and latent heat flux time series compare reasonably well, with overestimates of only 13% and 10%, respectively. In contrast, the NCEP model overestimates the bulk fluxes by 51% and 27%, respectively. The differences between the bulk estimates and those of the two models are due to different surface heat flux algorithms. It is shown that the roughness length formula used in the NCEP reanalysis project is inappropriate for moderate to high wind speeds. Its failings are acute for situations of large air–sea temperature difference and high wind speed, that is, for areas of high sensible heat fluxes such as the Labrador Sea, the Norwegian Sea, the Gulf Stream, and the Kuroshio. The new operational NCEP bulk algorithm is found to be more appropriate for such areas. It is concluded that surface turbulent flux fields from the ECMWF are within the bounds of observational uncertainty and therefore suitable for driving ocean models. This is in contrast to the surface flux fields from the NCEP reanalysis project, where the application of a more suitable algorithm to the model surface-layer meteorological data is recommended

Beschreibung: A simple stochastic atmosphere model is coupled to a realistic model of the North Atlantic Ocean. A north–south SST dipole, with its zero line centered along the subpolar front, influences the atmosphere model, which in turn forces the ocean model by surface fluxes related to the North Atlantic Oscillation. The coupled system exhibits a damped decadal oscillation associated with the adjustment through the ocean model to the changing surface forcing. The oscillation consists of a fast wind-driven, positive feedback of the ocean and a delayed negative feedback orchestrated by overturning circulation anomalies. The positive feedback turns out to be necessary to distinguish the coupled oscillation from that in a model without any influence from the ocean to the atmosphere. Using a novel diagnosing technique, it is possible to rule out the importance of baroclinic wave processes for determining the period of the oscillation, and to show the important role played by anomalous geostrophic advection in sustaining the oscillation.

Beschreibung: In this paper a decadal climate cycle in the North Atlantic that was derived from an extended-range integration with a coupled ocean–atmosphere general circulation model is described. The decadal mode shares many features with the observed decadal variability in the North Atlantic. The period of the simulated oscillation, however, is somewhat longer than that estimated from observations. While the observations indicate a period of about 12 yr, the coupled model simulation yields a period of about 17 yr. The cyclic nature of the decadal variability implies some inherent predictability at these timescales. The decadal mode is based on unstable air–sea interactions and must be therefore regarded as an inherently coupled mode. It involves the subtropical gyre and the North Atlantic oscillation. The memory of the coupled system, however, resides in the ocean and is related to horizontal advection and to the oceanic adjustment to low-frequency wind stress curl variations. In particular, it is found that variations in the intensity of the Gulf Stream and its extension are crucial to the oscillation. Although differing in details, the North Atlantic decadal mode and the North Pacific mode described by M. Latif and T. P. Barnett are based on the same fundamental mechanism: a feedback loop between the wind driven subtropical gyre and the extratropical atmospheric circulation.

Beschreibung: Output from an eddy-resolving model of the North Atlantic Ocean is used to estimate values for the thickness diffusivity κ appropriate to the Gent and McWilliams parameterization. The effect of different choices of rotational eddy fluxes on the estimated κ is discussed. Using the raw fluxes (no rotational flux removed), large negative values (exceeding −5000 m2 s−1) of κ are diagnosed locally, particularly in the Gulf Stream region and in the equatorial Atlantic. Removing a rotational flux based either on the suggestion of Marshall and Shutts or the more general theory of Medvedev and Greatbatch leads to a reduction of the negative values, but they are still present. The regions where κ 〈 0 correspond to regions where eddies are acting to increase, rather than decrease (as in baroclinic instability) the mean available potential energy. In the subtropical gyre, κ ranges between 500 and 2000 m2 s−1, rapidly decreasing to zero below the thermocline in all cases. Rotational fluxes and κ are also estimated using an optimization technique. In this case, |κ| can be reduced or increased by construction, but the regions where κ 〈 0 are still present and the optimized rotational fluxes also remain similar to a priori values given by the theoretical considerations. A previously neglected component (ν) of the bolus velocity is associated with the horizontal flux of buoyancy along, rather than across, the mean buoyancy contours. The ν component of the bolus velocity is interpreted as a streamfunction for eddy-induced advection, rather than diffusion, of mean isopycnal layer thickness, showing up when the lateral eddy fluxes cannot be described by isotropic diffusion only. All estimates show a similar large-scale pattern for ν, implying westward advection of isopycnal thickness over much of the subtropical gyre. Comparing ν with a mean streamfunction shows that it is about 10% of the mean in midlatitudes and even larger than the mean in the Tropics.

Beschreibung: Classical theory concerning theEliassen–Palmrelation is extended in this study to allowfor a unified treatment of midlatitude inertia–gravity waves (MIGWs), midlatitude Rossby waves (MRWs), and equatorial waves (EQWs). A conservation equation for what the authors call the impulse-bolus (IB) pseudomomentum is useful, because it is applicable to ageostrophic waves, and the associated three-dimensional flux is parallel to the direction of the group velocity of MRWs. The equation has previously been derived in an isentropic coordinate system or a shallow-water model. The authors make an explicit comparison of prognostic equations for the IB pseudomomentum vector and the classical energy-based (CE) pseudomomentum vector, assuming inviscid linear waves in a sufficiently weak mean flow, to provide a basis for the former quantity to be used in an Eulerian time-mean (EM) framework. The authors investigate what makes the three-dimensional fluxes in the IB and CE pseudomomentum equations look in different directions. It is found that the two fluxes are linked by a gauge transformation, previously unmentioned, associated with a divergence-form wave-induced pressure L. The quantity L vanishes for MIGWs and becomes nonzero for MRWs and EQWs, and it may be estimated using the virial theorem. Concerning the effect of waves on the mean flow, L represents an additional effect in the pressure gradient term of both (the three-dimensional versions of) the transformed EM momentum equations and the merged form of the EMmomentumequations, the latter of which is associated with the nonacceleration theorem.

Beschreibung: The turbulent dissipation rate ɛ is a key parameter to many oceanographic processes. Recently gliders have been increasingly used as a carrier for microstructure sensors. Compared to conventional ship-based methods, glider-based microstructure observations allow for long duration measurements under adverse weather conditions, and at lower costs. The incident water velocity U is an input parameter for the calculation of the dissipation rate. Since U can not be measured using the standard glider sensor setup, the parameter is normally computed from a steady-state glider flight model. As ɛ scales with U2 or U4, depending whether it is computed from temperature or shear microstructure, flight model errors can introduce a significant bias. This study is the first to use measurements of in-situ glider flight, obtained with a profiling Doppler velocity log and an electromagnetic current meter, to test and calibrate a flight model, extended to include inertial terms. Compared to a previously suggested flight model, the calibrated model removes a bias of approximately 1 cm s−1 in the incident water velocity, which translates to roughly a factor of 1.2 in estimates of the dissipation rate. The results further indicate that 90% of the estimates of the dissipation rate from the calibrated model are within a factor of 1.1 and 1.2 for measurements derived from microstructure temperature sensors and shear probes, respectively. We further outline the range of applicability of the flight model.

Beschreibung: Simple idealized layered models and primitive equation models show that the meridional gradient of the zonally averaged pressure has no direct relation with the meridional flow. This demonstrates a contradiction in an often-used parameterization in zonally averaged models. The failure of this parameterization reflects the inconsistency between the model of Stommel and Arons and the box model of Stommel, as previously pointed out by Straub. A new closure is proposed. The ocean is divided in two dynamically different regimes: a narrow western boundary layer and an interior ocean; zonally averaged quantities over these regions are considered. In the averaged equations three unknowns appear: the interior zonal pressure difference Delta p(i), the zonal pressure difference Delta p(b) of the boundary layer, and the zonal velocity us at the interface between the two regions. Here Delta p(i) is parameterized using a frictionless vorticity balance, Delta p(b), by the difference of the mean pressure in the interior and western boundary, and u(delta) by the mean zonal velocity of the western boundary layer. Zonally resolved models, a layer model, and a primitive equation model validate the new parameterization by comparing with the respective zonally averaged counterparts. It turns out that the zonally averaged models reproduce well the buoyancy distribution and the meridional flow in the zonally resolved model versions with respect to the mean and time changes.

Beschreibung: A simple point-vortex “heton” model is used to study localized ocean convection. In particular, the statistically steady state that is established when lateral buoyancy transfer, effected by baroclinic instability, offsets the localized surface buoyancy loss is investigated. Properties of the steady state, such as the statistically steady density anomaly of the convection region, are predicted using the hypothesis of a balance between baroclinic eddy transfer and the localized surface buoyancy loss. These predictions compare favorably with the values obtained through numerical integration of the heton model. The steady state of the heron model can be related to that in other convection scenarios considered in several recent studies by means of a generalized description of the localized convection. This leads to predictions of the equilibrium density anomalies in these scenarios, which concur with those obtained by other authors. Advantages of the heton model include its inviscid nature, emphasizing the independence of the fluxes affected by the baroclinic eddies from molecular processes, and its extreme economy, allowing a very large parameter space to be covered. This economy allows us to examine more complicated forcing scenarios: for example, forcing regions of varying shape. By increasing the ellipticity of the forcing region, the instability is modified by the shape and, as a result, no increase in lateral fluxes occurs despite the increased perimeter length. The parameterization of convective mixing by a redistribution of potential vorticity, implicit in the heton model, is corroborated; the heton model equilibrium state has analogous quantitative scaling behavior to that in models or laboratory experiments that resolve the vertical motions. The simplified dynamics of the heton model therefore allows the adiabatic advection resulting from baroclinic instability to be examined in isolation from vertical mixing and diffusive processes. These results demonstrate the importance of baroclinic instability in controlling the properties of a water mass generated by localized ocean convection. A complete parameterization of this process must therefore account for the fluxes induced by horizontal variations in surface buoyancy loss and affected by baroclinic instability.