Finite departure from convective quasi-equilibrium: periodic cycle and discharge-recharge mechanism

A simple self–contained theory is proposed for describing life cycles of convective systems as a discharge–recharge process. A closed description is derived for the dynamics of an ensemble of convective plumes based on an energy cycle. The system consists of prognostic equations for the cloud work function and the convective kinetic energy. The system can be closed by intro ducing a functional relationship between the convective kinetic energy and the cloud–base mass flux. The behaviour of this system is considered under a bulk simplification. Previous cloud–resolving mo delling as well as bulk statistical theories for ensemble convective systems suggest that a plausible relationship would be to assume that the convective kinetic energy is linearly proportional to the cloud–base mass flux. As a result, the system reduces to a nonlinear dynamical system with two dependent variables, the cloud–base mass flux and the cloud work function. The fully nonlinear solution of this system always represents a periodic cycle regardless of the initial condition under constant large–scale forcing. Importantly, the inclusion of energy dissipation in this model does not in itself lead the system to an equilibrium.

Iterative Methods for Equality-Constrained Least Squares Problems

We consider the linear equality-constrained least squares problem (LSE) of minimizing ${\|c - Gx\|}_2 $, subject to the constraint $Ex = p$. A preconditioned conjugate gradient method is applied to the Kuhn–Tucker equations associated with the LSE problem. We show that our method is well suited for structural optimization problems in reliability analysis and optimal design. Numerical tests are performed on an Alliant FX/8 multiprocessor and a Cray-X-MP using some practical structural analysis data.

Model studies of the interannual variability of the northern-hemisphere stratospheric winter circulation: the role of the quasi-biennial oscillation

A series of experiments are described that examine the sensitivity of the northern-hemisphere winter evolution to the equatorial quasi-biennial oscillation (QBO). The prime tool for the experiments is a stratosphere-mesosphere model. The model is integrated over many years with the modelled equatorial winds relaxed towards observed values in order to simulate a realistic QBO. In experiment A the equatorial winds are relaxed towards Singapore radiosonde observations in the height region 16-32 km. In contrast to previous modelling studies, the Holton-Tan relationship (warm/cold winters associated with easterly/westerly QBO winds in the lower stratosphere) is absent. However, in a second experiment (run B) in which the equatorial winds are relaxed towards rocketsonde data over the extended height range 16-58 km, a realistic Holton-Tan relationship is reproduced. A series of further studies are described that explore in more detail the sensitivity to various equatorial height regions and to the bottom-boundary forcing. The experiments suggest that the evolution of the northern-hemisphere winter circulation is sensitive to equatorial winds throughout the whole depth of the stratosphere and not just to the lower-stratospheric wind direction as previously assumed.

A variational method to retrieve the extinction profile in liquid clouds using multiple field-of-view lidar

Liquid clouds play a profound role in the global radiation budget but it is difficult to remotely retrieve their vertical profile. Ordinary narrow field-of-view (FOV) lidars receive a strong return from such clouds but the information is limited to the first few optical depths. Wideangle multiple-FOV lidars can isolate radiation scattered multiple times before returning to the instrument, often penetrating much deeper into the cloud than the singly-scattered signal. These returns potentially contain information on the vertical profile of extinction coefficient, but are challenging to interpret due to the lack of a fast radiative transfer model for simulating them. This paper describes a variational algorithm that incorporates a fast forward model based on the time-dependent two-stream approximation, and its adjoint. Application of the algorithm to simulated data from a hypothetical airborne three-FOV lidar with a maximum footprint width of 600m suggests that this approach should be able to retrieve the extinction structure down to an optical depth of around 6, and total opticaldepth up to at least 35, depending on the maximum lidar FOV. The convergence behavior of Gauss-Newton and quasi-Newton optimization schemes are compared. We then present results from an application of the algorithm to observations of stratocumulus by the 8-FOV airborne “THOR” lidar. It is demonstrated how the averaging kernel can be used to diagnose the effective vertical resolution of the retrieved profile, and therefore the depth to which information on the vertical structure can be recovered. This work enables exploitation of returns from spaceborne lidar and radar subject to multiple scattering more rigorously than previously possible.

Climate variability-observations, reconstructions, and model simulations for the Atlantic-European and Alpine region from 1500-2100 AD

A detailed analysis is undertaken of the Atlantic-European climate using data from 500-year-long proxy-based climate reconstructions, a long climate simulation with perpetual 1990 forcing, as well as two global and one regional climate change scenarios. The observed and simulated interannual variability and teleconnectivity are compared and interpreted in order to improve the understanding of natural climate variability on interannual to decadal time scales for the late Holocene. The focus is set on the Atlantic-European and Alpine regions during the winter and summer seasons, using temperature, precipitation, and 500 hPa geopotential height fields. The climate reconstruction shows pronounced interdecadal variations that appear to “lock” the atmospheric circulation in quasi-steady long-term patterns over multi-decadal periods controlling at least part of the temperature and precipitation variability. Different circulation patterns are persistent over several decades for the period 1500 to 1900. The 500-year-long simulation with perpetual 1990 forcing shows some substantial differences, with a more unsteady teleconnectivity behaviour. Two global scenario simulations indicate a transition towards more stable teleconnectivity for the next 100 years. Time series of reconstructed and simulated temperature and precipitation over the Alpine region show comparatively small changes in interannual variability within the time frame considered, with the exception of the summer season, where a substantial increase in interannual variability is simulated by regional climate models.

Decadal to multidecadal variability and the climate change background

Three prominent quasi-global patterns of variability and change are observed using the Met Office's sea surface temperature (SST) analysis and almost independent night marine air temperature analysis. The first is a global warming signal that is very highly correlated with global mean SST. The second is a decadal to multidecadal fluctuation with some geographical similarity to the El Niño–Southern Oscillation (ENSO). It is associated with the Pacific Decadal Oscillation (PDO), and its Pacific-wide manifestation has been termed the Interdecadal Pacific Oscillation (IPO). We present model investigations of the relationship between the IPO and ENSO. The third mode is an interhemispheric variation on multidecadal timescales which, in view of climate model experiments, is likely to be at least partly due to natural variations in the thermohaline circulation. Observed climatic impacts of this mode also appear in model simulations. Smaller-scale, regional atmospheric phenomena also affect climate on decadal to interdecadal timescales. We concentrate on one such mode, the winter North Atlantic Oscillation (NAO). This shows strong decadal to interdecadal variability and a correspondingly strong influence on surface climate variability which is largely additional to the effects of recent regional anthropogenic climate change. The winter NAO is likely influenced by both SST forcing and stratospheric variability. A full understanding of decadal changes in the NAO and European winter climate may require a detailed representation of the stratosphere that is hitherto missing in the major climate models used to study climate change.

The structure of the optimal income tax in the quasi-linear model

Existing numerical characterizations of the optimal income tax have been based on a limited number of model specifications. As a result, they do not reveal which properties are general. We determine the optimal tax in the quasi-linear model under weaker assumptions than have previously been used; in particular, we remove the assumption of a lower bound on the utility of zero consumption and the need to permit negative labor incomes. A Monte Carlo analysis is then conducted in which economies are selected at random and the optimal tax function constructed. The results show that in a significant proportion of economies the marginal tax rate rises at low skills and falls at high. The average tax rate is equally likely to rise or fall with skill at low skill levels, rises in the majority of cases in the centre of the skill range, and falls at high skills. These results are consistent across all the specifications we test. We then extend the analysis to show that these results also hold for Cobb-Douglas utility.

A ‘Boscastle-type’ quasi-stationary convective system over the UK Southwest Peninsula

An investigation is presented of a quasi-stationary convective system (QSCS) which occurred over the UK Southwest Peninsula on 21 July 2010. This system was remarkably similar in its location and structure to one which caused devastating flash flooding in the coastal village of Boscastle, Cornwall on 16 August 2004. However, in the 2010 case rainfall accumulations were around four times smaller and no flooding was recorded. The more extreme nature of the Boscastle case is shown to be related to three factors: (1) higher rain rates, associated with a warmer and moister tropospheric column and deeper convective clouds; (2) a more stationary system, due to slower evolution of the large-scale flow; and (3) distribution of the heaviest precipitation over fewer river catchments. Overall, however, the synoptic setting of the two events was broadly similar, suggesting that such conditions favour the development of QSCSs over the Southwest Peninsula. A numerical simulation of the July 2010 event was performed using a 1.5-km grid length configuration of the Met Office Unified Model. This reveals that convection was repeatedly initiated through lifting of low-level air parcels along a quasi-stationary coastal convergence line. Sensitivity tests are used to show that this convergence line was a sea breeze front which temporarily stalled along the coastline due to the retarding influence of an offshore-directed background wind component. Several deficiencies are noted in the 1.5-km model’s representation of the storm system, including delayed convective initiation; however, significant improvements are observed when the grid length is reduced to 500 m. These result in part from an improved representation of the convergence line, which enhances the associated low-level ascent allowing air parcels to more readily reach their level of free convection. The implications of this finding for forecasting convective precipitation are discussed.

'Quasi'-norm of an arithmetical convolution operator and the order of the Riemann zeta function

In this paper we study Dirichlet convolution with a given arithmetical function f as a linear mapping 'f that sends a sequence (an) to (bn) where bn = Pdjn f(d)an=d. We investigate when this is a bounded operator on l2 and ¯nd the operator norm. Of particular interest is the case f(n) = n¡® for its connection to the Riemann zeta function on the line 1, 'f is bounded with k'f k = ³(®). For the unbounded case, we show that 'f : M2 ! M2 where M2 is the subset of l2 of multiplicative sequences, for many f 2 M2. Consequently, we study the `quasi'-norm sup kak = T a 2M2 k'fak kak for large T, which measures the `size' of 'f on M2. For the f(n) = n¡® case, we show this quasi-norm has a striking resemblance to the conjectured maximal order of j³(® + iT )j for ® > 12 .