Dynamical control of the mesosphere by orographic and non-orographic gravity wave drag during the extended northern winters of 2006 and 2009

A version of the Canadian Middle Atmosphere Model (CMAM) that is nudged toward reanalysis data up to 1 hPa is used to examine the impacts of parameterized orographic and non-orographic gravity wave drag (OGWD and NGWD) on the zonal-mean circulation of the mesosphere during the extended northern winters of 2006 and 2009 when there were two large stratospheric sudden warmings. The simulations are compared to Aura Microwave Limb Sounder (MLS) observations of mesospheric temperature, carbon monoxide (CO) and derived zonal winds. The control simulation, which uses both OGWD and NGWD, is shown to be in good agreement with MLS. The impacts of OGWD and NGWD are assessed using simulations in which those sources of wave drag are removed. In the absence of OGWD the mesospheric zonal winds in the months preceding the warmings are too strong, causing increased mesospheric NGWD, which drives excessive downwelling, resulting in overly large lower mesospheric values of CO prior to the warming. NGWD is found to be most important following the warmings when the underlying westerlies are too weak to allow much vertical propagation of the orographic gravity waves to the mesosphere. NGWD is primarily responsible for driving the circulation that results in the descent of CO from the thermosphere following the warmings. Zonal mean mesospheric winds and temperatures in all simulations are shown to be strongly constrained by (i.e. slaved to) the stratosphere. Finally, it is demonstrated that the responses to OGWD and NGWD are non-additive due to their dependence and influence on the background winds and temperatures.

Idealised simulations of sting-jet cyclones

An idealised modelling study of sting-jet cyclones is presented. Sting jets are descending mesoscale jets that occur in some extratropical cyclones and produce localised regions of strong low-level winds in the frontal fracture region. Moist baroclinic lifecycle (LC1) simulations are performed with modifications to produce cyclones resembling observed sting-jet cyclones. A sting jet exists in the idealised control cyclone with similar characteristics to the sting jet in a simulation of windstorm Gudrun (a confirmed sting-jet case). Unlike in windstorm Gudrun, a low-level layer of strong moist static stability prohibits the descent of the strong winds from above the boundary layer to the surface in the idealised case. Conditional symmetric instability (CSI) exists in the cloud head and dissipates as the sting jet leaves the cloud head and descends. The descending, initially moist, sting-jet trajectories consistently have negative or near-zero saturated moist potential vorticity but moist static stability and inertial stability, consistent with CSI release; the moist static stability becomes negative during the period of most rapid descent, by which time the air is relatively dry implying conditional instability release is unlikely. Sensitivity experiments show that the existence of the sting jet is robust to changes in the initial state, and that the initial tropospheric static stability significantly impacts the descent rate of the sting jet. Inertial and conditional instability are probably being released in the experiment with the weakest initial static stability. This suggests that sting jets can arise through the release of all three instabilities associated with negative saturated moist potential vorticity. While evaporative cooling occurs along the sting-jet trajectories, a sensitivity experiment with evaporation effects turned off shows no significant change to the wind strength or descent rate of the sting jet implying that instability release is the dominant sting-jet driving mechanism.

Satellite observations and modeling of transport in the upper troposphere through the lower mesosphere during the 2006 major stratospheric sudden warming

An unusually strong and prolonged stratospheric sudden warming (SSW) in January 2006 was the first major SSW for which globally distributed long-lived trace gas data are available covering the upper troposphere through the lower mesosphere. We use Aura Microwave Limb Sounder (MLS), Atmospheric Chemistry Experiment-Fourier Transform Spectrometer (ACE-FTS) data, the SLIMCAT Chemistry Transport Model (CTM), and assimilated meteorological analyses to provide a comprehensive picture of transport during this event. The upper tropospheric ridge that triggered the SSW was associated with an elevated tropopause and layering in trace gas profiles in conjunction with stratospheric and tropospheric intrusions. Anomalous poleward transport (with corresponding quasi-isentropic troposphere-to-stratosphere exchange at the lowest levels studied) in the region over the ridge extended well into the lower stratosphere. In the middle and upper stratosphere, the breakdown of the polar vortex transport barrier was seen in a signature of rapid, widespread mixing in trace gases, including CO, H2O, CH4 and N2O. The vortex broke down slightly later and more slowly in the lower than in the middle stratosphere. In the middle and lower stratosphere, small remnants with trace gas values characteristic of the pre-SSW vortex lingered through the weak and slow recovery of the vortex. The upper stratospheric vortex quickly reformed, and, as enhanced diabatic descent set in, CO descended into this strong vortex, echoing the fall vortex development. Trace gas evolution in the SLIMCAT CTM agrees well with that in the satellite trace gas data from the upper troposphere through the middle stratosphere. In the upper stratosphere and lower mesosphere, the SLIMCAT simulation does not capture the strong descent of mesospheric CO and H2O values into the reformed vortex; this poor CTM performance in the upper stratosphere and lower mesosphere results primarily from biases in the diabatic descent in assimilated analyses.

Determination of eddy-diffusivity in the lowermost stratosphere

We present a 2D-advection-diffusion model that simulates the main transport pathways influencing tracer distributions in the lowermost stratosphere (LMS). The model describes slow diabatic descent of aged stratospheric air, vertical (cross-isentropic) and horizontal (along isentropes) diffusion within the LMS and across the tropopause using equivalent latitude and potential temperature coordinates. Eddy diffusion coefficients parameterize the integral effect of dynamical processes leading to small scale turbulence and mixing. They were specified by matching model simulations to observed CO distributions. Interestingly, the model suggests mixing across isentropes to be more important than horizontal mixing across surfaces of constant equivalent latitude, shining new light on the interplay between various transport mechanisms in the LMS. The model achieves a good description of the small scale tracer features at the tropopause with squared correlation coefficients R2 = 0.72…0.94.

Interdecadal enhancement of the Walker circulation over the Tropical Pacific in the late 1990s

The Walker circulation is one of the major components of the large-scale tropical atmospheric circulation and variations in its strength are critical to equatorial Pacific Ocean circulation. It has been argued in the literature that during the 20th century the Walker circulation weakened, and that this weakening was attributable to anthropogenic climate change. By using updated observations, we show that there has been a rapid interdecadal enhancement of the Walker circulation since the late 1990s. Associated with this enhancement is enhanced precipitation in the tropical western Pacific, anomalous westerlies in the upper troposphere, descent in the central and eastern tropical Pacific, and anomalous surface easterlies in the western and central tropical Pacific. The characteristics of associated oceanic changes are a strengthened thermocline slope and an enhanced zonal SST gradient across the tropical Pacific. Many characteristics of these changes are similar to those associated with the mid-1970s climate shift with an opposite sign. We also show that the interdecadal variability of the Walker circulation in the tropical Pacific is inversely correlated to the interdecadal variability of the zonal circulation in the tropical Atlantic. An enhancement of the Walker circulation in the tropical Pacific is associated with a weakening zonal circulation in the tropical Atlantic and vise versa, implying an inter-Atlantic-Pacific connection of the zonal overturning circulation variation. Whether these recent changes will be sustained is not yet clear, but our research highlights the importance of understanding the interdecadal variability, as well as the long-term trends, that influence tropical circulation.

Generation of a buoyancy-driven coastal current by an Antarctic Polynya

Descent and spreading of high salinity water generated by salt rejection during sea ice formation in an Antarctic coastal polynya is studied using a hydrostatic, primitive equation three-dimensional ocean model called the Proudman Oceanographic Laboratory Coastal Ocean Modeling System (POLCOMS). The shape of the polynya is assumed to be a rectangle 100 km long and 30 km wide, and the salinity flux into the polynya at its surface is constant. The model has been run at high horizontal spatial resolution (500 m), and numerical simulations reveal a buoyancy-driven coastal current. The coastal current is a robust feature and appears in a range of simulations designed to investigate the influence of a sloping bottom, variable bottom drag, variable vertical turbulent diffusivities, higher salinity flux, and an offshore position of the polynya. It is shown that bottom drag is the main factor determining the current width. This coastal current has not been produced with other numerical models of polynyas, which may be because these models were run at coarser resolutions. The coastal current becomes unstable upstream of its front when the polynya is adjacent to the coast. When the polynya is situated offshore, an unstable current is produced from its outset owing to the capture of cyclonic eddies. The effect of a coastal protrusion and a canyon on the current motion is investigated. In particular, due to the convex shape of the coastal protrusion, the current sheds a dipolar eddy.

Distinguishing the cold conveyor belt and sting jet air streams in an intense extratropical cyclone

Strong winds equatorwards and rearwards of a cyclone core have often been associated with two phenomena, the cold conveyor belt (CCB) jet and sting jets. Here, detailed observations of the mesoscale structure in this region of an intense cyclone are analysed. The {\it in-situ} and dropsonde observations were obtained during two research flights through the cyclone during the DIAMET (DIAbatic influences on Mesoscale structures in ExTratropical storms) field campaign. A numerical weather prediction model is used to link the strong wind regions with three types of ``air streams'', or coherent ensembles of trajectories: two types are identified with the CCB, hooking around the cyclone center, while the third is identified with a sting jet, descending from the cloud head to the west of the cyclone. Chemical tracer observations show for the first time that the CCB and sting jet air streams are distinct air masses even when the associated low-level wind maxima are not spatially distinct. In the model, the CCB experiences slow latent heating through weak resolved ascent and convection, while the sting jet experiences weak cooling associated with microphysics during its subsaturated descent. Diagnosis of mesoscale instabilities in the model shows that the CCB passes through largely stable regions, while the sting jet spends relatively long periods in locations characterized by conditional symmetric instability (CSI). The relation of CSI to the observed mesoscale structure of the bent-back front and its possible role in the cloud banding is discussed.

Fast identification algorithms for Gaussian process model

A class identification algorithms is introduced for Gaussian process(GP)models.The fundamental approach is to propose a new kernel function which leads to a covariance matrix with low rank,a property that is consequently exploited for computational efficiency for both model parameter estimation and model predictions.The objective of either maximizing the marginal likelihood or the Kullback–Leibler (K–L) divergence between the estimated output probability density function(pdf)and the true pdf has been used as respective cost functions.For each cost function,an efficient coordinate descent algorithm is proposed to estimate the kernel parameters using a one dimensional derivative free search, and noise variance using a fast gradient descent algorithm. Numerical examples are included to demonstrate the effectiveness of the new identification approaches.

Genetic analysis reveals a novel locus in Vicia faba decoupling pigmentation in the flower from that in the extra-floral nectaries

The genetics of the stipule spot pigmentation (SSP) in faba bean (Vicia faba L.) was studied using four inbred lines, of which Disco/2 was zero-tannin (zt2) with colourless stipule spots, ILB938/2 was normal-tannin (ZT2) with colourless stipule spots, and both Aurora/2 and Mélodie/2 were ZT2 with coloured stipule spots. Crosses Mélodie/2 × ILB 938/2, Mélodie/2 × Disco/2, ILB 938/2 × Aurora/2 and ILB 938/2 × Disco/2 (A, B, C and D, respectively) were prepared, along with reciprocals and backcrosses, and advanced through single-seed descent. All F1 hybrid plants had pigmented stipule spots, and in the F2 generation, the segregation ratio fit 3 coloured:1 colourless in crosses A, B and C and 9:7 in cross D. In the F3 generation, the ratio fit 5:3 in crosses A and C and 25:39 in cross D, and in the F4 generation, 9:7 in cross A. SSP was linked to the zero-tannin characteristics (white flower) only in cross B. The results show that coloured stipule spot is dominant to colourless and that colouration is determined by two unlinked complementary recessive genes. We propose the symbols ssp2 for the gene associated with zt2 in Disco/2 and ssp1 for the gene not associated with tannin content in ILB938/2. The novel ssp1 locus was mapped at F5 in cross ‘A’ using Medicago truncatula-derived single-nucleotide polymorphism and was on chromosome 1 of faba bean, in a well-conserved region of M. truncatula chromosome 5 containing some candidate Myb and basic helix–loop–helix transcription factor genes.

Financial services for the poor within the formal/informal nexus in sub-Saharan Africa

This paper draws upon fieldwork undertaken across Kenya, Zambia, Mozambique and South Africa to present a reflective overview of the use of financial services amongst the poorest members of society. It considers the role that access to a portfolio of financial products and services may have as a contributory factor in poverty alleviation, but also how inappropriate use of these mechanisms may exacerbate a descent into poverty. This work draws upon the notions of poverty pools and the rise of fall of low income households in and out of poverty, alongside the contributory nature of vicious cycles of economic and political poverty. Drawing on fieldwork experiences it presents a synopsis of the types of financial mechanisms commonly in use on the African continent, as well as examples of public, private and civil society partnerships that are producing services specifically tailored for those in extreme and absolute poverty.

Estimation of Gaussian process regression model using probability distance measures

A new class of parameter estimation algorithms is introduced for Gaussian process regression (GPR) models. It is shown that the integration of the GPR model with probability distance measures of (i) the integrated square error and (ii) Kullback–Leibler (K–L) divergence are analytically tractable. An efficient coordinate descent algorithm is proposed to iteratively estimate the kernel width using golden section search which includes a fast gradient descent algorithm as an inner loop to estimate the noise variance. Numerical examples are included to demonstrate the effectiveness of the new identification approaches.

Design of delay-tolerant space–time codes with limited feedback

In cooperative communication networks, owing to the nodes' arbitrary geographical locations and individual oscillators, the system is fundamentally asynchronous. Such a timing mismatch may cause rank deficiency of the conventional space-time codes and, thus, performance degradation. One efficient way to overcome such an issue is the delay-tolerant space-time codes (DT-STCs). The existing DT-STCs are designed assuming that the transmitter has no knowledge about the channels. In this paper, we show how the performance of DT-STCs can be improved by utilizing some feedback information. A general framework for designing DT-STC with limited feedback is first proposed, allowing for flexible system parameters such as the number of transmit/receive antennas, modulated symbols, and the length of codewords. Then, a new design method is proposed by combining Lloyd's algorithm and the stochastic gradient-descent algorithm to obtain optimal codebook of STCs, particularly for systems with linear minimum-mean-square-error receiver. Finally, simulation results confirm the performance of the newly designed DT-STCs with limited feedback.

The effects of springtime mid-latitude storms on trace gas composition determined from the MACC reanalysis

The relationship between springtime air pollution transport of ozone (O3) and carbon monoxide (CO) and mid-latitude cyclones is explored for the first time using the Monitoring Atmospheric Composition and Climate (MACC) reanalysis for the period 2003–2012. In this study, the most intense spring storms (95th percentile) are selected for two regions, the North Pacific (NP) and the North Atlantic (NA). These storms (∼60 storms over each region) often track over the major emission sources of East Asia and eastern North America. By compositing the storms, the distributions of O3 and CO within a "typical" intense storm are examined. We compare the storm-centered composite to background composites of "average conditions" created by sampling the reanalysis data of the previous year to the storm locations. Mid-latitude storms are found to redistribute concentrations of O3 and CO horizontally and vertically throughout the storm. This is clearly shown to occur through two main mechanisms: (1) vertical lifting of CO-rich and O3-poor air isentropically, from near the surface to the mid- to upper-troposphere in the region of the warm conveyor belt; and (2) descent of O3-rich and CO-poor air isentropically in the vicinity of the dry intrusion, from the stratosphere toward the mid-troposphere. This can be seen in the composite storm's life cycle as the storm intensifies, with area-averaged O3 (CO) increasing (decreasing) between 200 and 500 hPa. The influence of the storm dynamics compared to the background environment on the composition within an area around the storm center at the time of maximum intensity is as follows. Area-averaged O3 at 300 hPa is enhanced by 50 and 36% and by 11 and 7.6% at 500 hPa for the NP and NA regions, respectively. In contrast, area-averaged CO at 300 hPa decreases by 12% for NP and 5.5% for NA, and area-averaged CO at 500 hPa decreases by 2.4% for NP while there is little change over the NA region. From the mid-troposphere, O3-rich air is clearly seen to be transported toward the surface, but the downward transport of CO-poor air is not discernible due to the high levels of CO in the lower troposphere. Area-averaged O3 is slightly higher at 1000 hPa (3.5 and 1.8% for the NP and NA regions, respectively). There is an increase of CO at 1000 hPa for the NP region (3.3%) relative to the background composite and a~slight decrease in area-averaged CO for the NA region at 1000 hPa (-2.7%).

Mechanisms underlying temperature extremes in Iberia: a Lagrangian perspective

The mechanisms underlying the occurrence of temperature extremes in Iberia are analysed considering a Lagrangian perspective of the atmospheric flow, using 6-hourly ERA-Interim reanalysis data for the years 1979–2012. Daily 2-m minimum temperatures below the 1st percentile and 2-m maximum temperatures above the 99th percentile at each grid point over Iberia are selected separately for winter and summer. Four categories of extremes are analysed using 10-d backward trajectories initialized at the extreme temperature grid points close to the surface: winter cold (WCE) and warm extremes (WWE), and summer cold (SCE) and warm extremes (SWE). Air masses leading to temperature extremes are first transported from the North Atlantic towards Europe for all categories. While there is a clear relation to large-scale circulation patterns in winter, the Iberian thermal low is important in summer. Along the trajectories, air mass characteristics are significantly modified through adiabatic warming (air parcel descent), upper-air radiative cooling and near-surface warming (surface heat fluxes and radiation). High residence times over continental areas, such as over northern-central Europe for WCE and, to a lesser extent, over Iberia for SWE, significantly enhance these air mass modifications. Near-surface diabatic warming is particularly striking for SWE. WCE and SWE are responsible for the most extreme conditions in a given year. For WWE and SCE, strong temperature advection associated with important meridional air mass transports are the main driving mechanisms, accompanied by comparatively minor changes in the air mass properties. These results permit a better understanding of mechanisms leading to temperature extremes in Iberia.

An Ensemble Smoother with Error Estimates

A smoother introduced earlier by van Leeuwen and Evensen is applied to a problem in which real obser vations are used in an area with strongly nonlinear dynamics. The derivation is new , but it resembles an earlier derivation by van Leeuwen and Evensen. Again a Bayesian view is taken in which the prior probability density of the model and the probability density of the obser vations are combined to for m a posterior density . The mean and the covariance of this density give the variance-minimizing model evolution and its errors. The assumption is made that the prior probability density is a Gaussian, leading to a linear update equation. Critical evaluation shows when the assumption is justified. This also sheds light on why Kalman filters, in which the same ap- proximation is made, work for nonlinear models. By reference to the derivation, the impact of model and obser vational biases on the equations is discussed, and it is shown that Bayes’ s for mulation can still be used. A practical advantage of the ensemble smoother is that no adjoint equations have to be integrated and that error estimates are easily obtained. The present application shows that for process studies a smoother will give superior results compared to a filter , not only owing to the smooth transitions at obser vation points, but also because the origin of features can be followed back in time. Also its preference over a strong-constraint method is highlighted. Further more, it is argued that the proposed smoother is more efficient than gradient descent methods or than the representer method when error estimates are taken into account