The influence of stratospheric vortex displacements and splits on surface climate

A strong link exists between stratospheric variability and anomalous weather patterns at the earth’s surface. Specifically, during extreme variability of the Arctic polar vortex termed a “weak vortex event,” anomalies can descend from the upper stratosphere to the surface on time scales of weeks. Subsequently the outbreak of cold-air events have been noted in high northern latitudes, as well as a quadrupole pattern in surface temperature over the Atlantic and western European sectors, but it is currently not understood why certain events descend to the surface while others do not. This study compares a new classification technique of weak vortex events, based on the distribution of potential vorticity, with that of an existing technique and demonstrates that the subdivision of such events into vortex displacements and vortex splits has important implications for tropospheric weather patterns on weekly to monthly time scales. Using reanalysis data it is found that vortex splitting events are correlated with surface weather and lead to positive temperature anomalies over eastern North America of more than 1.5 K, and negative anomalies over Eurasia of up to −3 K. Associated with this is an increase in high-latitude blocking in both the Atlantic and Pacific sectors and a decrease in European blocking. The corresponding signals are weaker during displacement events, although ultimately they are shown to be related to cold-air outbreaks over North America. Because of the importance of stratosphere–troposphere coupling for seasonal climate predictability, identifying the type of stratospheric variability in order to capture the correct surface response will be necessary.

The role of anomalous SST and surface fluxes over the Southeastern North Atlantic in the explosive development of windstorm Xynthia

In late February 2010 the extraordinary windstorm Xynthia crossed over Southwestern and Central Europe and caused severe damage, affecting particularly the Spanish and French Atlantic coasts. The storm was embedded in uncommon large-scale atmospheric and boundary conditions prior to and during its development, namely enhanced sea surface temperatures (SST) within the low-level entrainment zone of air masses, an unusual southerly position of the polar jet stream, and a remarkable split jet structure in the upper troposphere. To analyse the processes that led to the rapid intensification of this exceptional storm originating close to the subtropics (30°N), the sensitivity of the cyclone intensification to latent heat release is determined using the regional climate model COSMO-CLM forced with ERA-Interim data. A control simulation with observed SST shows that moist and warm air masses originating from the subtropical North Atlantic were involved in the cyclogenesis process and led to the formation of a vertical tower with high values of potential vorticity (PV). Sensitivity studies with reduced SST or increased laminar boundary roughness for heat led to reduced surface latent heat fluxes. This induced both a weaker and partly retarded development of the cyclone and a weakening of the PV-tower together with reduced diabatic heating rates, particularly at lower and mid levels. We infer that diabatic processes played a crucial role during the phase of rapid deepening of Xynthia and thus to its intensity over the Southeastern North Atlantic. We suggest that windstorms like Xynthia may occur more frequently under future climate conditions due to the warming SSTs and potentially enhanced latent heat release, thus increasing the windstorm risk for Southwestern Europe.

The dichotomous structure of the warm conveyor belt

The warm conveyor belt (WCB) of an extratropical cyclone generally splits into two branches. One branch (WCB1) turns anticyclonically into the downstream upper-level tropospheric ridge, while the second branch (WCB2) wraps cyclonically around the cyclone centre. Here, the WCB split in a typical North Atlantic cold-season cyclone is analysed using two numerical models: the Met Office Unified Model and the COSMO model. The WCB flow is defined using off-line trajectory analysis. The two models represent the WCB split consistently. The split occurs early in the evolution of the WCB with WCB1 experiencing maximum ascent at lower latitudes and with higher moisture content than WCB2. WCB1 ascends abruptly along the cold front where the resolved ascent rates are greatest and there is also line convection. In contrast, WCB2 remains at lower levels for longer before undergoing saturated large-scale ascent over the system's warm front. The greater moisture in WCB1 inflow results in greater net potential temperature change from latent heat release, which determines the final isentropic level of each branch. WCB1 also exhibits lower outflow potential vorticity values than WCB2. Complementary diagnostics in the two models are utilised to study the influence of individual diabatic processes on the WCB. Total diabatic heating rates along the WCB branches are comparable in the two models with microphysical processes in the large-scale cloud schemes being the major contributor to this heating. However, the different convective parameterisation schemes used by the models cause significantly different contributions to the total heating. These results have implications for studies on the influence of the WCB outflow in Rossby wave evolution and breaking. Key aspects are the net potential temperature change and the isentropic level of the outflow which together will influence the relative mass going into each WCB branch and the associated negative PV anomalies at the tropopause-level flow.

The ability of CMIP5 models to simulate North Atlantic extratropical cyclones

The ability of the climate models participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5) to simulate North Atlantic extratropical cyclones in winter [December–February (DJF)] and summer [June–August (JJA)] is investigated in detail. Cyclones are identified as maxima in T42 vorticity at 850 hPa and their propagation is tracked using an objective feature-tracking algorithm. By comparing the historical CMIP5 simulations (1976–2005) and the ECMWF Interim Re-Analysis (ERA-Interim; 1979–2008), the authors find that systematic biases affect the number and intensity of North Atlantic cyclones in CMIP5 models. In DJF, the North Atlantic storm track tends to be either too zonal or displaced southward, thus leading to too few and weak cyclones over the Norwegian Sea and too many cyclones in central Europe. In JJA, the position of the North Atlantic storm track is generally well captured but some CMIP5 models underestimate the total number of cyclones. The dynamical intensity of cyclones, as measured by either T42 vorticity at 850 hPa or mean sea level pressure, is too weak in both DJF and JJA. The intensity bias has a hemispheric character, and it cannot be simply attributed to the representation of the North Atlantic large- scale atmospheric state. Despite these biases, the representation of Northern Hemisphere (NH) storm tracks has improved since CMIP3 and some CMIP5 models are able of representing well both the number and the intensity of North Atlantic cyclones. In particular, some of the higher-atmospheric-resolution models tend to have a better representation of the tilt of the North Atlantic storm track and of the intensity of cyclones in DJF.

Maintenance of the African easterly jet

The maintenance of the African easterly jet (AEJ) has been examined using a zonally symmetric general circulation model with simple parametrizations. It is shown that the AEJ is maintained in association with two diabatically forced meridional circulations: one associated with surface fluxes and dry convection in the Saharan heat-low region and one associated with deep moist convection in the intertropical convergence zone equatorward of this. the heat-low heating, which reaches the height of the AEJ around 700 mb, is particularly important in maintaining the AEJ and its associated meridional gradients in potential vorticity. It is concluded that the mean observed AEJ results from a combination of the diabatically forced meridional circulations which maintain it and easterly waves which weaken it.

Development of convection along the SPCZ within a Madden-Julian oscillation

A subtropical Rossby-wave propagation mechanism is proposed to account for the poleward and eastward progression of intraseasonal convective anomalies along the South Pacific convergence zone (SPCZ) that is observed in a significant proportion of Madden–Julian oscillations (MJOs). Large-scale convection, associated with an MJO, is assumed to be already established over the Indonesian region. The latent heating associated with this convection forces an equatorial Rossby-wave response with an upper-tropospheric anticyclone centred over, or slightly to the west of, the convection. Large potential-vorticity (PV) gradients, associated with the subtropical jet and the tropopause, lie just poleward of the anticyclone, and large magnitude PV air is advected equatorwards on the eastern side of the anticyclone. This ‘high’ PV air, or upper-tropospheric trough, is far enough off the equator that it has associated strong horizontal temperature gradients, and it induces deep ascent on its eastern side, at a latitude of about 15–30°. If this deep ascent is over a region susceptible to deep convection, such as the SPCZ, then convection may be forced or triggered. Hence convection develops along the SPCZ as a forced response to convection over Indonesia. The response mechanism is essentially one of subtropical Rossby-wave propagation. This hypothesis is based on a case study of a particularly strong MJO in early 1988, and is tested by idealized modelling studies. The mechanism may also be relevant to the existence of the mean SPCZ, as a forced response to mean Indonesian convection.

Efficient fully nonlinear data assimilation for geophysical fluid dynamics

A potential problem with Ensemble Kalman Filter is the implicit Gaussian assumption at analysis times. Here we explore the performance of a recently proposed fully nonlinear particle filter on a high-dimensional but simplified ocean model, in which the Gaussian assumption is not made. The model simulates the evolution of the vorticity field in time, described by the barotropic vorticity equation, in a highly nonlinear flow regime. While common knowledge is that particle filters are inefficient and need large numbers of model runs to avoid degeneracy, the newly developed particle filter needs only of the order of 10-100 particles on large scale problems. The crucial new ingredient is that the proposal density cannot only be used to ensure all particles end up in high-probability regions of state space as defined by the observations, but also to ensure that most of the particles have similar weights. Using identical twin experiments we found that the ensemble mean follows the truth reliably, and the difference from the truth is captured by the ensemble spread. A rank histogram is used to show that the truth run is indistinguishable from any of the particles, showing statistical consistency of the method.

Can polar lows be objectively identified and tracked in the ECMWF operational analysis and the ERA-Interim reanalysis?

Polar lows are maritime meso-cyclones associated with intense surface wind speeds and oceanic heat fluxes at high latitudes. The ability of the ERA-Interim (ERAI) reanalysis to represent polar lows in the North Atlantic is assessed by comparing ERAI and the ECMWF operational analysis for the period 2008-2011. First, the representation of a set of satellite observed polar lows over the Norwegian and Barents Seas in the operational analysis and ERAI is analysed. Then, the possibility of directly identifying and tracking the polar lows in the operational analysis and ERAI is explored using a tracking algorithm based on 850 hPa vorticity with objective identification criteria on cyclone dynamical intensity and atmospheric static stability. All but one of the satellite observed polar lows with a lifetime of at least 6 hours have an 850 hPa vorticity signature of a co-located mesocyclone in both the operational analysis and ERAI for most of their life cycles. However, the operational analysis has vorticity structures that better resemble the observed cloud patterns and stronger surface wind speed intensities compared to those in ERAI. By applying the objective identification criteria, about 55% of the satellite observed polar lows are identified and tracked in ERAI, while this fraction increases to about 70% in the operational analysis. Particularly in ERAI, the remaining observed polar lows are mainly not identified because they have too weak wind speed and vorticity intensity compared to the tested criteria. The implications of the tendency of ERAI to underestimate the polar low dynamical intensity for future studies of polar lows is discussed.

An inter-hemispheric comparison of the tropical storm response to global warming

Model studies do not agree on future changes in tropical cyclone (TC) activity on regional scales. We aim to shed further light on the distribution, frequency, intensity, and seasonality of TCs that society can expect at the end of the twenty-first century in the Southern hemisphere (SH). Therefore, we investigate TC changes simulated by the atmospheric model ECHAM5 with T213 (~60 km) horizontal resolution. We identify TCs in present-day (20C; 1969–1990) and future (21C; 2069–2100) time slice simulations, using a tracking algorithm based on vorticity at 850 hPa. In contrast to the Northern hemisphere (NH), where tropical storm numbers reduce by 6 %, there is a more dramatic 22 % reduction in the SH, mainly in the South Indian Ocean. While an increase of static stability in 21C may partly explain the reduction in tropical storm numbers, stabilization cannot alone explain the larger SH drop. Large-scale circulation changes associated with a weakening of the Tropical Walker Circulation are hypothesized to cause the strong decrease of cyclones in the South Indian Ocean. In contrast the decrease found over the South Pacific appears to be partly related to increased vertical wind shear, which is possibly associated with an enhanced meridional sea surface temperature gradient. We find the main difference between the hemispheres in changes of the tropical cyclones of intermediate strength with an increase in the NH and a decrease in the SH. In both hemispheres the frequency of the strongest storms increases and the frequency of the weakest storms decreases, although the increase in SH intense storms is marginal.

Power transfer in velocity-vortex acceleration

A theoretical model is presented of an electron acceleration-as-oscillator method derived from the work of Joseph Larmor unified with J. Clerk Maxwell’s theory of vorticity for the displacement of radiation into free-space at an antenna interface.

Curl free pressure gradients over orography in a solution of the fully compressible Euler equations with implicit treatment of acoustic and gravity waves

Steep orography can cause noisy solutions and instability in models of the atmosphere. A new technique for modelling flow over orography is introduced which guarantees curl free gradients on arbitrary grids, implying that the pressure gradient term is not a spurious source of vorticity. This mimetic property leads to better hydrostatic balance and better energy conservation on test cases using terrain following grids. Curl-free gradients are achieved by using the co-variant components of velocity over orography rather than the usual horizontal and vertical components. In addition, gravity and acoustic waves are treated implicitly without the need for mean and perturbation variables or a hydrostatic reference profile. This enables a straightforward description of the implicit treatment of gravity waves. Results are presented of a resting atmosphere over orography and the curl-free pressure gradient formulation is advantageous. Results of gravity waves over orography are insensitive to the placement of terrain-following layers. The model with implicit gravity waves is stable in strongly stratified conditions, with N∆t up to at least 10 (where N is the Brunt-V ̈ais ̈al ̈a frequency). A warm bubble rising over orography is simulated and the curl free pressure gradient formulation gives much more accurate results for this test case than a model without this mimetic property.

Climate at high-obliquity

The question of climate at high obliquity is raised in the context of both exoplanet studies (e.g. habitability) and paleoclimates studies (evidence for low-latitude glaciation during the Neoproterozoic and the ”Snowball Earth” hypothesis). States of high obliquity, φ, are distinctive in that, for φ ≥54◦, the poles receive more solar radiation in the annual mean than the Equator, opposite to the present day situation. In addition, the seasonal cycle of insolation is extreme, with the poles alternatively “facing” the sun and sheltering in the dark for months. The novelty of our approach is to consider the role of a dynamical ocean in controlling the surface climate at high obliquity, which in turn requires understanding of the surface winds patterns when temperature gradients are reversed. To address these questions, a coupled ocean-atmosphere-sea ice GCM configured on an aquaplanet is employed. Except for the absence of topography and modified obliquity, the set-up is Earth-like. Two large obliquities φ, 54◦ and 90◦, are compared to today’s Earth value, φ=23.5◦. Three key results emerge at high obliquity: 1) despite reversed temper- ature gradients, mid-latitudes surface winds are westerly and trade winds exist at the equator (as for φ=23.5◦) although the westerlies are confined to the summer hemisphere, 2) a habitable planet is possible with mid-latitude temperatures in the range 300-280 K and 3) a stable climate state with an ice cap limited to the equatorial region is unlikely. We clarify the dynamics behind these features (notably by an analysis of the potential vorticity structure and conditions for baroclinic instability of the atmosphere). Interestingly, we find that the absence of a stable partially glaciated state is critically linked to the absence of ocean heat transport during winter, a feature ultimately traced back to the high seasonality of baroclinic instability conditions in the atmosphere.

Simulation of the global ENSO–Tropical cyclone teleconnection by a high-resolution coupled general circulation model

This study assesses the influence of the El Niño–Southern Oscillation (ENSO) on global tropical cyclone activity using a 150-yr-long integration with a high-resolution coupled atmosphere–ocean general circulation model [High-Resolution Global Environmental Model (HiGEM); with N144 resolution: ~90 km in the atmosphere and ~40 km in the ocean]. Tropical cyclone activity is compared to an atmosphere-only simulation using the atmospheric component of HiGEM (HiGAM). Observations of tropical cyclones in the International Best Track Archive for Climate Stewardship (IBTrACS) and tropical cyclones identified in the Interim ECMWF Re-Analysis (ERA-Interim) are used to validate the models. Composite anomalies of tropical cyclone activity in El Niño and La Niña years are used. HiGEM is able to capture the shift in tropical cyclone locations to ENSO in the Pacific and Indian Oceans. However, HiGEM does not capture the expected ENSO–tropical cyclone teleconnection in the North Atlantic. HiGAM shows more skill in simulating the global ENSO–tropical cyclone teleconnection; however, variability in the Pacific is overpronounced. HiGAM is able to capture the ENSO–tropical cyclone teleconnection in the North Atlantic more accurately than HiGEM. An investigation into the large-scale environmental conditions, known to influence tropical cyclone activity, is used to further understand the response of tropical cyclone activity to ENSO in the North Atlantic and western North Pacific. The vertical wind shear response over the Caribbean is not captured in HiGEM compared to HiGAM and ERA-Interim. Biases in the mean ascent at 500 hPa in HiGEM remain in HiGAM over the western North Pacific; however, a more realistic low-level vorticity in HiGAM results in a more accurate ENSO–tropical cyclone teleconnection.

On the structure and dynamics of Indian monsoon depressions

ERA-Interim reanalysis data from the past 35 years have been used with a newly-developed feature tracking algorithm to identify Indian monsoon depressions originating in or near the Bay of Bengal. These were then rotated, centralised and combined to give a fully three-dimensional 106-depression composite structure – a considerably larger sample than any previous detailed study on monsoon depressions and their structure. Many known features of depression structure are confirmed, particularly the existence of a maximum to the southwest of the centre in rainfall and other fields, and a westward axial tilt in others. Additionally, the depressions are found to have significant asymmetry due to the presence of the Himalayas; a bimodal mid-tropospheric potential vorticity core; a separation into thermally cold- (~–1.5K) and neutral- (~0K) cores near the surface with distinct properties; and that the centre has very large CAPE and very small CIN. Variability as a function of background state has also been explored, with land/coast/sea, diurnal, ENSO, active/break and Indian Ocean Dipole contrasts considered. Depressions are found to be markedly stronger during the active phase of the monsoon, as well as during La Niña. Depressions on land are shown to be more intense and more tightly constrained to the central axis. A detailed schematic diagram of a vertical cross-section through a composite depression is also presented, showing its inherent asymmetric structure.

Friction in mid-latitude cyclones: an Ekman-PV mechanism

The mechanism by which the atmospheric boundary layer reduces the intensity of mid-latitude cyclones is investigated. It is demonstrated that two alternative theories, Ekman pumping and the baroclinic potential vorticity (PV) mechanism, in fact act in union to maximize the spin-down. Ekman pumping aids the ventilation of PV from the boundary layer, and shapes the resulting PV anomaly into one of increased static stability. PV inversion techniques are used to demonstrate how this anomaly reduces the coupling between the upper- and lower-levels within the cyclone, reducing the growth rate.