ALBERT

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: 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: A nested‐grid ocean circulation modeling system is used to study the response of Lunenburg Bay in Nova Scotia, Canada, to local wind‐forcing, tides, remotely generated waves, and buoyancy forcing in the summer and fall of 2003. Quantitative comparisons between observations and model results demonstrate that the modeling system reproduces reasonably well the observed sea level, temperature, salinity, and currents in the bay. Numerical results reveal that the spatial and temporal variability of temperature and salinity in the bay during the study period is mainly forced by the local wind stress and surface heat/freshwater fluxes, with some contribution from tidal circulation. In particular, the local heat balance on the monthly timescale is dominated by cooling due to vertical advection and warming due to horizontal advection and net surface heat flux, while high‐frequency variations (timescales of 1–30 days) are mainly associated with vertical advection, i.e., wind‐induced upwelling and downwelling. There is also a strong baroclinic throughflow over the deep water region outside Lunenburg Bay that is strongly influenced by wind‐forcing. The vertically integrated momentum balance analysis indicates a modified geostrophic balance on the monthly timescale and longer, and is dominated by the pressure term and wind minus bottom stress in the high‐frequency band.

Description: Northern Europe and the UK experienced an exceptionally warm and wet winter in 2019/20, driven by an anomalously positive North Atlantic Oscillation (NAO). This positive NAO was well forecast by several seasonal forecast systems, suggesting that this winter the NAO was highly predictable at seasonal lead times. A very strong positive Indian Ocean dipole (IOD) event was also observed at the start of winter. Here we use composite analysis and model experiments, to show that the IOD was a key driver of the observed positive NAO. Using model experiments that perturb the Indian Ocean initial conditions, two teleconnection pathways of the IOD to the north Atlantic emerge: a tropospheric teleconnection pathway via a Rossby wave train travelling from the Indian Ocean over the Pacific and Atlantic, and a stratospheric teleconnection pathway via the Aleutian region and the stratospheric polar vortex. These pathways are similar to those for the El Niño Southern Oscillation link to the north Atlantic which are already well documented. The anomalies in the north Atlantic jet stream location and strength, and the associated precipitation anomalies over the UK and northern Europe, as simulated by the model IOD experiments, show remarkable agreement with those forecast and observed.

Description: There is a zonally oriented teleconnection pattern over the high-latitude Eurasian continent, which is maintained through baroclinic energy conversion. In this study, we investigate the unique features of the maintenance mechanism of this teleconnection. It is found that the baroclinic energy conversion is most efficient in both the mid-troposphere and the lower troposphere, and that the baroclinic energy conversion in the lower troposphere is comparable to that in the mid-troposphere. Further results indicate that the basic state plays a crucial role in the baroclinic energy conversion. For both the mid and lower troposphere, the atmospheric stability is low and the Coriolis parameter is large over high-latitude Eurasia, favoring strong baroclinic energy conversion. Particularly, in the lower troposphere, the atmospheric stability exhibits a clear land-sea contrast, favoring baroclinic energy conversion over the continents rather than the oceans. Furthermore, in the lower troposphere, the in-phase configuration of the meridional wind and temperature anomalies, which results from the strong meridional gradient of mean temperature around the north edge of the Eurasian continent, also significantly contributes to baroclinic energy conversion. This study highlights the role of the basic state of temperature rather than zonal wind in maintaining the high-latitude teleconnection through baroclinic energy conversion.

Description: A method using a linear shallow water model is presented for decomposing the temporal variability of the barotropic streamfunction in a high‐resolution ocean model. The method is based in the vertically‐averaged momentum equations and is applied to the time series of annual mean streamfunction from the model configuration VIKING20 for the northern North Atlantic. An important result is the role played by the nonlinear advection terms in VIKING20 for driving transport. The method is illustrated by examining how the Gulf Stream transport in the recirculation region responds to the winter North Atlantic Oscillation (NAO). While no statistically significant response is found in the year overlapping with the winter NAO index, there is a tendency for the Gulf Stream transport to increase as the NAO becomes more positive. This becomes significant in lead years 1 and 2 when the mean flow advection (MFA) and eddy momentum flux (EMF) contributions, associated with nonlinear momentum advection, dominate. Only after 2 years, does the potential energy (PE) term, associated with the density field, start to play a role and it is only after 5 years that the transport dependence on the NAO ceases to be significant. It is also shown that the PE contribution to the transport streamfunction has significant memory of up to 5 years in the Labrador and Irminger Seas. However, it is only around the northern rim of these seas that VIKING20 and the transport reconstruction exhibit similar memory. This is due to masking by the MFA and EMF contributions.

Description: Equatorial deep jets (EDJ) are vertically stacked, downward propagating zonal currents that alternate in direction with depth. In the tropical Atlantic, they have been shown to influence both surface conditions and tracer variability. Despite their importance, the EDJ are absent in most ocean models. Here we show that EDJ can be generated in an idealized ocean model when the model is driven only by the convergence of the meridional flux of intraseasonal zonal momentum diagnosed from a companion model run driven by steady wind forcing, corroborating the recent theory that intraseasonal momentum flux convergence maintains the EDJ. Additionally, the EDJ in our model nonlinearly generate mean zonal currents at intermediate depths that show similarities in structure to the observed circulation in the deep equatorial Atlantic, indicating their importance for simulating the tropical ocean mean state.

Description: Mesoscale eddies can be strengthened by the absorption of submesoscale eddies resulting from mixed-layer baroclinic instabilities. This is shown for mesoscale eddies in the Agulhas Current system by investigating the kinetic energy cascade with a spectral and a coarse-graining approach in two model simulations of the Agulhas region. One simulation resolves mixed-layer baroclinic instabilities and one does not. When mixed-layer baroclinic instabilities are included, the largest submesoscale near-surface fluxes occur in winter-time in regions of strong mesoscale activity for upscale as well as downscale directions. The forward cascade at the smallest resolved scales occurs mainly in frontogenetic regions in the upper 30 m of the water column. In the Agulhas ring path, the forward cascade changes to an inverse cascade at a typical scale of mixed-layer eddies (15 km). At the same scale, the largest sources of the upscale flux occur. After the winter, the maximum of the upscale flux shifts to larger scales. Depending on the region, the kinetic energy reaches the mesoscales in spring or early summer aligned with the maximum of mesoscale kinetic energy. This indicates the importance of submesoscale flows for the mesoscale seasonal cycle. A case study shows that the underlying process is the mesoscale absorption of mixed-layer eddies. When mixed-layer baroclinic instabilities are not included in the simulation, the open-ocean upscale cascade in the Agulhas ring path is almost absent. This contributes to a 20 %-reduction of surface kinetic energy at mesoscales larger than 100 km when submesoscale dynamics are not resolved by the model.