ALBERT

Description: In the present study, the influence of some major tropical modes of variability on northern hemisphere regional blocking frequency variability during boreal winter is investigated. Reanalysis data and an ensemble experiment with the ECMWF model using relaxation towards the ERA-Interim reanalysis data inside the tropics are used. The tropical modes under investigation are El Niño Southern Oscillation (ENSO), the Madden-Julian Oscillation (MJO) and the upper tropospheric equatorial zonal-mean zonal wind . An early (late) MJO phase refers to the part of the MJO cycle when enhanced (suppressed) precipitation occurs over the western Indian Ocean and suppressed (enhanced) precipitation occurs over the Maritime Continent and the western tropical Pacific. Over the North Pacific sector, it is found that enhanced (suppressed) high latitude blocking occurs in association with El Niño (La Niña) events, late (early) MJO phases and westerly (easterly) . Over central to southern Europe and the east Atlantic, it is found that late MJO phases, as well as a suppressed MJO are leading to enhanced blocking frequency. Furthermore, early (late) MJO phases are followed by blocking anomalies over the western North Atlantic region, similar to those associated with a positive (negative) North Atlantic Oscillation. Over northern Europe, the easterly (westerly) phase of is associated with enhanced (suppressed) blocking. These results are largely confirmed by both the reanalysis and the model experiment.

Description: The Eddy Kinetic Energy (EKE) associated with the Subtropical Countercurrent (STCC) in the western subtropical South Pacific is known to exhibit substantial seasonal and decadal variability. Using an eddy-permitting ocean general circulation model, which is able to reproduce the observed, salient features of the seasonal cycles of shear, stratification, baroclinic production and the associated EKE, we investigate the decadal changes of EKE. We show that the STCC region exhibits, uniquely among the subtropical gyres of the world’s oceans, significant, atmospherically forced, decadal EKE variability. The decadal variations are driven by changing vertical shear between the STCC in the upper 300 m and the South Equatorial Current below, predominantly caused by variations in STCC strength associated with a changing meridional density gradient. In the 1970s, an increased meridional density gradient results in EKE twice as large as in later decades in the model. Utilizing sensitivity experiments, decadal variations in the wind field are shown to be the essential driver. Local wind stress curl anomalies associated with the Interdecadal Pacific Oscillation (IPO) lead to up- and downwelling of the thermocline, inducing strengthening or weakening of the STCC and the associated EKE. Additionally, remote wind stress curl anomalies in the eastern subtropical South Pacific, which are not related to the IPO, generate density anomalies that propagate westward as Rossby waves and can account for up to 30–40 % of the density anomalies in the investigated region.

Description: The mechanisms involved in setting the annual cycle of the Florida Current transport are revisited using an adjoint model approach. Adjoint sensitivities of the Florida Current transport to wind stress reproduce a realistic seasonal cycle with an amplitude of ~1.2 Sv (1 Sv ≡ 106 m3 s−1). The annual cycle is predominantly determined by wind stress forcing and related coastal upwelling (downwelling) north of the Florida Strait along the shelf off the North American coast. Fast barotropic waves propagate these anomalies southward and reach the Florida Strait within a month, causing an amplitude of ~1 Sv. Long baroclinic planetary Rossby waves originating from the interior are responsible for an amplitude of ~0.8 Sv but have a different phase. The sensitivities corresponding to the first baroclinic mode propagate westward and are highly influenced by topography. Considerable sensitivities are only found west of the Mid-Atlantic Ridge, with maximum values at the western shelf edge. The second baroclinic mode also has an impact on the Florida Current variability, but only when a mean flow is present. A second-mode wave train propagates southwestward from the ocean bottom on the western side of the Mid-Atlantic Ridge between ~36° and 46°N and at Flemish Cap, where the mean flow interacts with topography, to the surface. Other processes such as baroclinic waves along the shelf and local forcing within the Florida Strait are of minor importance.

Description: The residual effect of surface gravity waves on mean flows in the upper ocean is investigated using thickness weighted mean (TWM) theory applied in a vertically Lagrangian and horizontally Eulerian coordinate system. Depth-dependent equations for the conservation of volume, momentum, and energy are derived. These equations allow for (i) finite amplitude fluid motions, (ii) the horizontal divergence of currents and (iii) a concise treatment of both the kinematic and viscous boundary conditions at the sea surface. Under the assumptions of steady and monochromatic waves and a uniform turbulent viscosity, the TWM momentum equations are used to illustrate the pressure- and viscosity-induced momentum fluxes through the surface that are implicit in previous studies of the wave-induced modification of the classical Ekman spiral problem. The TWM approach clarifies, in particular, the surface momentum flux associated with the so-called virtual wave stress of Longuet-Higgins. Overall the TWM framework can be regarded as an alternative to the three-dimensional Lagrangian mean framework of Pierson. Moreover the TWM framework can be used to include the residual effect of surface waves in large-scale circulation models. In specific models that carry the TWM velocity appropriate for advecting tracers as their velocity variable, the turbulent viscosity term should be modified so that the viscosity acts only on the Eulerian mean velocity.

Description: This study focuses on an important aspect of air–sea interaction in models, namely, large-scale, spurious heat fluxes due to false pathways of the Gulf Stream and North Atlantic Current (NAC) in the “storm formation region” south and east of Newfoundland. Although high-resolution eddy-resolving models show some improvement in this respect, results are sensitive to poorly understood, subgrid-scale processes for which there is currently no complete, physically based parameterization. A simple method to correct an ocean general circulation model (OGCM), acting as a practical substitute for a physically based parameterization, is explored: the recently proposed “semiprognostic method,” a technique for adiabatically adjusting flow properties of a hydrostatic OGCM. The authors show that application of the method to an eddy-permitting model of the North Atlantic Ocean yields more realistic flow patterns and watermass characteristics in the Gulf Stream and NAC regions; in particular, spurious surface heat fluxes are reduced. Four simple modifications to the method are proposed, and their benefits are demonstrated. The modifications successfully account for three drawbacks of the original method: reduced geostrophic wave speeds, damped mesoscale eddy activity, and spurious interaction with topography. It is argued that use of a corrected (eddy permitting) OGCM in a coupled modeling system for simulating present climate (as now becomes possible because of increasing computer power) should lead to a more realistic simulation in regions of strong air–sea interaction as compared with that obtained with an uncorrected model. The method is also well suited for the simulation of the uptake and transport of passive tracers, such as anthropogenic carbon dioxide or components of ecosystem models.

Description: We perform eddy-resolving and high-vertical-resolution numerical simulations of the circulation in an idealized equatorial Atlantic Ocean in order to explore the formation of the deep equatorial circulation (DEC) in this basin. Unlike in previous studies, the deep equatorial intraseasonal variability (DEIV) that is believed to be the source of the DEC is generated internally by instabilities of the upper ocean currents. Two main simulations are discussed: Solution 1, configured with a rectangular basin and with wind forcing that is zonally and temporally uniform; and Solution 2, with realistic coastlines and with an annual cycle of wind forcing varying zonally. Somewhat surprisingly, Solution 1 produces the more realistic DEC: The large-vertical-scale currents (Equatorial Intermediate Currents or EICs) are found over a large zonal portion of the basin, and the small-vertical-scale equatorial currents (Equatorial Deep Jets or EDJs) form low-frequency, quasi-resonant, baroclinic equatorial basin modes with phase propagating mostly downward, consistent with observations. We demonstrate that both types of currents arise from the rectification of DEIV, consistent with previous theories. We also find that the EDJs contribute to maintaining the EICs, suggesting that the nonlinear energy transfer is more complex than previously thought. In Solution 2, the DEC is unrealistically weak and less spatially coherent than in the first simulation probably because of its weaker DEIV. Using intermediate solutions, we find that the main reason for this weaker DEIV is the use of realistic coastlines in Solution 2. It remains to be determined, what needs to be modified or included to obtain a realistic DEC in the more realistic configuration.

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

Description: A generalization of the transformed Eulerian and temporal residual means is presented. The new formulation uses rotational fluxes of buoyancy, and the full hierarchy of statistical density moments, to reduce the cross-isopycnal eddy flux to the physically relevant component associated with the averaged water mass properties. The resulting eddy-induced diapycnal diffusivity vanishes for adiabatic, statistically steady flow, and is related to either the growth or decay of mesoscale density variance and/or the covariance between small-scale forcing (mixing) and density fluctuations, such as that associated with the irreversible removal of density variance by dissipation. The relationship between the new formulation and previous approaches is described and is illustrated using results from an eddying channel model. The formalism is quite general and applies to all kinds of averaging and to any tracer (not just density).

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.