ALBERT

Description: On 19 October 2000, Hurricane Michael merged with an approaching baroclinic trough over the western North Atlantic Ocean south of Nova Scotia. As the hurricane moved over cooler sea surface temperatures (SSTs; less than 25°C), it intensified to category-2 intensity on the Saffir–Simpson hurricane scale [maximum sustained wind speeds of 44 m s−1 (85 kt)] while tapping energy from the baroclinic environment. The large “hybrid” storm made landfall on the south coast of Newfoundland with maximum sustained winds of 39 m s−1 (75 kt) causing moderate damage to coastal communities east of landfall. Hurricane Michael presented significant challenges to weather forecasters. The fundamental issue was determining which of two cyclones (a newly formed baroclinic low south of Nova Scotia or the hurricane) would become the dominant circulation center during the early stages of the extratropical transition (ET) process. Second, it was difficult to predict the intensity of the storm at landfall owing to competing factors: 1) decreasing SSTs conducive to weakening and 2) the approaching negatively tilted upper-level trough, favoring intensification. Numerical hindcast simulations using the limited-area Mesoscale Compressible Community model with synthetic vortex insertion (cyclone bogus) prior to the ET of Hurricane Michael led to a more realistic evolution of wind and pressure compared to running the model without vortex insertion. Specifically, the mesoscale model correctly simulates the hurricane as the dominant circulation center early in the transition process, versus the baroclinic low to its north, which was the favored development in the runs not employing vortex insertion. A suite of experiments is conducted to establish the sensitivity of the ET to various initial conditions, lateral driving fields, domain sizes, and model parameters. The resulting storm tracks and intensities fall within the range of the operational guidance, lending support to the possibility of improving numerical forecasts using synthetic vortex insertion prior to ET in such a model.

Description: Previous attempts to derive the depth-dependent expression of the radiation stress have lead to a debate concerning (i) the applicability of Mellor’s approach to a sloping bottom, (ii) the introduction of the delta function at the mean sea surface in the later papers by Mellor, and (iii) a wave-induced pressure term derived in several recent studies. The authors use an equation system in vertically Lagrangian and horizontally Eulerian (VL) coordinates suitable for a concise treatment of the surface boundary, and obtain an expression for the depth-dependent radiation stress that is consistent with the vertically-integrated expression given by Longuet-Higgins and Stewart. Concerning (i)-(iii) in the above, the difficulty of handling a sloping bottom disappears when wave-averaged momentum equations in the VL coordinates are written for the development of (not the Lagrangian mean velocity but) the Eulerian mean velocity. There is also no delta function at the sea surface in the expression for the depth-dependent radiation stress. The connection between the wave-induced pressure term in the recent studies and the depth-dependent radiation stress term is easily shown by rewriting the pressure-based form stress term in the thickness-weighted-mean (TWM) momentum equations as a velocity-based term which contains the time derivative of the pseudomomentum in the TWM framework.

Description: The use of a coupled ocean/atmosphere/sea-ice model to hindcast (i.e. historical forecast) recent climate variability is described and illustrated for the cases of the 1976/77 and 1998/99 climate shift events in the Pacific. The initialization is achieved by running the coupled model in partially coupled mode whereby global observed wind stress anomalies are used to drive the ocean/sea-ice component of the coupled model while maintaining the thermodynamic coupling between the ocean/sea-ice and atmosphere components. Here we show that hindcast experiments can successfully capture many features associated with the 1976/77 and 1998/99 climate shifts. For instance, hindcast experiments started from the beginning of 1976 can capture sea surface temperature (SST) warming in the central-eastern equatorial Pacific and the positive phase of the Pacific Decadal Oscillation (PDO) throughout the 9 years following the 1976/77 climate shift, including the deepening of the Aleutian low pressure system. Hindcast experiments started from the beginning of 1998 can also capture part of the anomalous conditions during the 4 years after the 1998/99 climate. We argue that the dynamical adjustment of heat content anomalies that are present in the initial conditions in the tropics is important for the successful hindcast of the two climate shifts.

Description: The Equatorial Deep Jets (EDJs) are an ubiquitous feature of the equatorial oceans; in the Atlantic Ocean, they are the dominant mode of interannual variability of the zonal flow at intermediate depth. On the basis of more than 10 years of moored observations of zonal velocity at 23°W, the vertically propagating EDJs are best described as superimposed oscillations of the 13th to the 23th baroclinic modes with a dominant oscillation period for all modes of 1650 days. This period is close to the resonance period of the respective gravest equatorial basin mode for the dominant vertical modes 16 and 17. It is argued that since the equatorial basin mode is composed of linear equatorial waves, a linear reduced gravity model can be employed for each baroclinic mode, driven by spatially homogeneous zonal forcing oscillating with the EDJ period. The fit of the model solutions to observations at 23°W yields a basin wide reconstruction of the EDJs and the associated vertical structure of their forcing. From the resulting vertical profile of mean power input and vertical energy flux on the equator, it follows that the EDJs are locally maintained over a considerable depth range, from 500-2500 m, with the maximum power input and vertical energy flux at 1300 m. The strong dissipation closely ties the apparent vertical propagation of energy to the vertical distribution of power input and, together with the EDJs’ prevailing downward phase propagation, require the phase of the forcing of the EDJs to propagate downward.

Description: Seasonal variability of the tropical Atlantic circulation is dominated by the annual cycle, but semi-annual variability is also pronounced, despite weak forcing at that period. Here we use multi-year, full-depth velocity measurements from the central equatorial Atlantic to analyze the vertical structure of annual and semi-annual variations of zonal velocity. A baroclinic modal decomposition finds that the annual cycle is dominated by the 4th mode and the semi-annual cycle by the 2nd mode. Similar local behavior is found in a high-resolution general circulation model. This simulation reveals that the annual and semi-annual cycles of the respective dominant baroclinic modes are associated with characteristic basin-wide structures. Using an idealized linear reduced-gravity model to simulate the dynamics of individual baroclinic modes, it is shown that the observed circulation variability can be explained by resonant equatorial basin modes. Corollary simulations of the reduced-gravity model with varying basin geometry (i.e. square basin versus realistic coastlines) or forcing (i.e. spatially uniform versus spatially variable wind) show a structural robustness of the simulated basin modes. A main focus of this study is the seasonal variability of the Equatorial Undercurrent (EUC) as identified in recent observational studies. Main characteristics of the observed EUC including seasonal variability of transport, core depth, and maximum core velocity can be explained by the linear superposition of the dominant equatorial basin modes as obtained from the reduced-gravity model.

Description: By performing two sets of high-resolution atmospheric general circulation model (AGCM) experiments, we find that the atmospheric response to a sea surface temperature (SST) anomaly in the extratropical North Pacific is sensitive to decadal variations of the background SST on which the SST anomaly is superimposed. The response in the first set of experiments, in which the SST anomaly is superimposed on the observed daily SST of 1981-1990, strongly differs from the response in the second experiment, in which the same SST anomaly is superimposed on the observed daily SST of 1991-2000. The atmospheric response over the North Pacific during 1981-1990 is eddy-mediated, equivalent barotropic and concentrated in the east. In contrast, the atmospheric response during 1991-2000 is weaker and strongest in the west. The results are discussed in terms of Rossby wave dynamics, with the proposed primary wave source switching from baroclinic eddy vorticity forcing over the eastern North Pacific in 1981-1990 to mean flow divergence over the western North Pacific in 1991-2000. The wave source changes are linked to the decadal reduction of daily SST variability over the eastern North Pacific and strengthening of the Oyashio Extension front over the western North Pacific. Thus, both daily and frontal aspects of the background SST variability in determining the atmospheric response to extratropical North Pacific SST anomalies are emphasized by our AGCM experiments.

Description: 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.

Description: Variations in the global tropospheric zonal mean zonal wind ([U]) during boreal winter are investigated using Rotated Empirical Orthogonal Functions applied to monthly means. The first two modes correspond to the Northern and Southern Annular Mode and modes 3 and 4 represent variability in the tropics. One is related to El Niño Southern Oscillation and the other has variability that is highly correlated with the time series of [U] at 150 hPa between 5°N and 5°S ([U150]E) and is related to activity of the Madden-Julian Oscillation. The extratropical response to [U150]E is investigated using linear regressions of 500 hPa geopotential height onto the [U150]E time series. We make use of reanalysis data and of the ensemble mean output from a relaxation experiment using the European Center for Medium Range Weather Forecasts model in which the tropical atmosphere is relaxed towards reanalysis data. The regression analysis reveals that a shift of the Aleutian low and a wave train across the North Atlantic are associated with [U150]E. We find that the subtropical waveguides and the link between the North Pacific and North Atlantic are stronger during the easterly phase of [U150]E. The wave train over the North Atlantic is associated with Rossby wave sources over the subtropical North Pacific and North America. Finally, we show that a linear combination of both [U150]E and the Quasi Biennial Oscillation in the lower stratosphere can explain the circulation anomalies of the anomalously cold European winter of 1962/63 when both were in an extreme easterly phase.

Description: There is an ongoing discussion in the community concerning the wave-averaged momentum equations in the hybrid vertically Lagrangian and horizontally Eulerian (VL) framework and, in particular, the form stress term (representing the residual effect of pressure perturbations) which is thought to restrict the handling of higher order waves in terms of a perturbation expansion. The present study shows that the traditional pressure-based form stress term can be transformed into a set of terms that do not contain any pressure quantities but do contain the time derivative of a wave-induced velocity. This wave-induced velocity is referred to as the pseudomomentum in the VL framework, as it is analogous to the generalized pseudomomentum in Andrews and McIntyre. This enables the second expression for the wave-averaged momentum equations in the VL framework (this time for the development of the total transport velocity minus the VL pseudomomentum) to be derived together with the vortex force. The velocity-based expression of the form stress term also contains the residual effect of the turbulent viscosity, which is useful for understanding the dissipation of wave energy leading to transfer of momentum from waves to circulation. It is found that the concept of the virtual wave stress of Longuet-Higgins is applicable to quite general situations: it does not matter whether there is wind forcing or not, the waves can have slow variations, and the viscosity coefficient can vary in the vertical. These results provide a basis for revisiting the surface boundary condition used in numerical circulation models.

Description: The dynamical origin of the spectral and autocorrelation structure of annular variability in the troposphere is investigated by a deductive approach. Specifically, the structure of the power spectrum and autocorrelation function of the zonal-mean geopotential is analyzed for the case of a quasigeostrophic spherical atmosphere subject to a white noise mechanical forcing applied in a single Hough mode and concentrated at a particular level in the vertical, with vertically uniform Newtonian cooling and Rayleigh drag concentrated at a rigid lower boundary. Analytic expressions for the power spectrum are presented together with expressions for an approximate red noise (i.e., a Lorentzian-shaped) power spectrum. It is found that for an infinitely deep atmosphere the power spectrum can be well approximated by a red noise process for the first few Hough modes (associated with large Rossby heights), provided the distance from the forcing is not larger than about one Rossby height. When a frictional rigid lower boundary is included, however, the approximation is generally bad. The high-frequency part of the power spectrum exhibits near-exponential behavior and the autocorrelation function shows a transition from a rapid decay at short lags to a much slower decay at longer lags, if the thermal and mechanical damping time scales are sufficiently well separated. Since observed annular variability exhibits the same characteristics, the above results lead to the hypothesis that these characteristics may, to some extent, be intrinsic to the linear zonal-mean response problem—although the need for an additional contribution from eddy feedbacks is also implied by the results.