Wave-activity conservation laws and stability theorems for semi-geostrophic dynamics: part 2. pseudoenergy-based theory

This paper represents the second part of a study of semi-geostrophic (SG) geophysical fluid dynamics. SG dynamics shares certain attractive properties with the better known and more widely used quasi-geostrophic (QG) model, but is also a good prototype for balanced models that are more accurate than QG dynamics. The development of such balanced models is an area of great current interest. The goal of the present work is to extend a central body of QG theory, concerning the evolution of disturbances to prescribed basic states, to SG dynamics. Part 1 was based on the pseudomomentum; Part 2 is based on the pseudoenergy. A pseudoenergy invariant is a conserved quantity, of second order in disturbance amplitude relative to a prescribed steady basic state, which is related to the time symmetry of the system. We derive such an invariant for the semi-geostrophic equations, and use it to obtain: (i) a linear stability theorem analogous to Arnol'd's ‘first theorem’; and (ii) a small-amplitude local conservation law for the invariant, obeying the group-velocity property in the WKB limit. The results are analogous to their quasi-geostrophic forms, and reduce to those forms in the limit of small Rossby number. The results are derived for both the f-plane Boussinesq form of semi-geostrophic dynamics, and its extension to β-plane compressible flow by Magnusdottir & Schubert. Novel features particular to semi-geostrophic dynamics include apparently unnoticed lateral boundary stability criteria. Unlike the boundary stability criteria found in the first part of this study, however, these boundary criteria do not necessarily preclude the construction of provably stable basic states. The interior semi-geostrophic dynamics has an underlying Hamiltonian structure, which guarantees that symmetries in the system correspond naturally to the system's invariants. This is an important motivation for the theoretical approach used in this study. The connection between symmetries and conservation laws is made explicit using Noether's theorem applied to the Eulerian form of the Hamiltonian description of the interior dynamics.

Wave-activity conservation laws and stability theorems for semi-geostrophic dynamics: part 1. pseudomomentum-based theory

There exists a well-developed body of theory based on quasi-geostrophic (QG) dynamics that is central to our present understanding of large-scale atmospheric and oceanic dynamics. An important question is the extent to which this body of theory may generalize to more accurate dynamical models. As a first step in this process, we here generalize a set of theoretical results, concerning the evolution of disturbances to prescribed basic states, to semi-geostrophic (SG) dynamics. SG dynamics, like QG dynamics, is a Hamiltonian balanced model whose evolution is described by the material conservation of potential vorticity, together with an invertibility principle relating the potential vorticity to the advecting fields. SG dynamics has features that make it a good prototype for balanced models that are more accurate than QG dynamics. In the first part of this two-part study, we derive a pseudomomentum invariant for the SG equations, and use it to obtain: (i) linear and nonlinear generalized Charney–Stern theorems for disturbances to parallel flows; (ii) a finite-amplitude local conservation law for the invariant, obeying the group-velocity property in the WKB limit; and (iii) a wave-mean-flow interaction theorem consisting of generalized Eliassen–Palm flux diagnostics, an elliptic equation for the stream-function tendency, and a non-acceleration theorem. All these results are analogous to their QG forms. The pseudomomentum invariant – a conserved second-order disturbance quantity that is associated with zonal symmetry – is constructed using a variational principle in a similar manner to the QG calculations. Such an approach is possible when the equations of motion under the geostrophic momentum approximation are transformed to isentropic and geostrophic coordinates, in which the ageostrophic advection terms are no longer explicit. Symmetry-related wave-activity invariants such as the pseudomomentum then arise naturally from the Hamiltonian structure of the SG equations. We avoid use of the so-called ‘massless layer’ approach to the modelling of isentropic gradients at the lower boundary, preferring instead to incorporate explicitly those boundary contributions into the wave-activity and stability results. This makes the analogy with QG dynamics most transparent. This paper treats the f-plane Boussinesq form of SG dynamics, and its recent extension to β-plane, compressible flow by Magnusdottir & Schubert. In the limit of small Rossby number, the results reduce to their respective QG forms. Novel features particular to SG dynamics include apparently unnoticed lateral boundary stability criteria in (i), and the necessity of including additional zonal-mean eddy correlation terms besides the zonal-mean potential vorticity fluxes in the wave-mean-flow balance in (iii). In the companion paper, wave-activity conservation laws and stability theorems based on the SG form of the pseudoenergy are presented.

Multiple stationary solutions of an irradiated slab

A mathematical model describing the heat budget of an irradiated medium is introduced. The one-dimensional form of the equations and boundary conditions are presented and analysed. Heat transport at one face of the slab occurs by absorption (and reflection) of an incoming beam of short-wave radiation with a fraction of this radiation penetrating into the body of the slab, a diffusive heat flux in the slab and a prescribed incoming heat flux term. The other face of the slab is immersed in its own melt and is considered to be a free surface. Here, temperature continuity is prescribed and evolution of the surface is determined by a Stefan condition. These boundary conditions are flexible enough to describe a range of situations such as a laser shining on an opaque medium, or the natural environment of polar sea ice or lake ice. A two-stream radiation model is used which replaces the simple Beer’s law of radiation attenuation frequently used for semi-infinite domains. The stationary solutions of the governing equations are sought and it is found that there exists two possible stationary solutions for a given set of boundary conditions and a range of parameter choices. It is found that the existence of two stationary solutions is a direct result of the model of radiation absorption, due to its effect on the albedo of the medium. A linear stability analysis and numerical calculations indicate that where two stationary solutions exist, the solution corresponding to a larger thickness is always stable and the solution corresponding to a smaller thickness is unstable. Numerical simulations reveal that when there are two solutions, if the slab is thinner than the smaller stationary thickness it will melt completely, whereas if the slab is thicker than the smaller stationary thickness it will evolve toward the larger stationary thickness. These results indicate that other mechanisms (e.g. wave-induced agglomeration of crystals) are necessary to grow a slab from zero initial thickness in the parameter regime that yields two stationary solutions.

Integrable evolution equations in time-dependent domains

We discuss the implementation of a method of solving initial boundary value problems in the case of integrable evolution equations in a time-dependent domain. This method is applied to a dispersive linear evolution equation with spatial derivatives of arbitrary order and to the defocusing nonlinear Schrödinger equation, in the domain l(t)

Enhancement of the aerosol direct radiative effect by semi-volatile aerosol components: airborne measurements in North-Western Europe

A case study of atmospheric aerosol measurements exploring the impact of the vertical distribution of aerosol chemical composition upon the radiative budget in North-Western Europe is presented. Sub-micron aerosol chemical composition was measured by an Aerodyne Aerosol Mass Spectrometer (AMS) on both an airborne platform and a ground-based site at Cabauw in the Netherlands. The examined period in May 2008 was characterised by enhanced pollution loadings in North-Western Europe and was dominated by ammonium nitrate and Organic Matter (OM). Both ammonium nitrate and OM were observed to increase with altitude in the atmospheric boundary layer. This is primarily attributed to partitioning of semi-volatile gas phase species to the particle phase at reduced temperature and enhanced relative humidity. Increased ammonium nitrate concentrations in particular were found to strongly increase the ambient scattering potential of the aerosol burden, which was a consequence of the large amount of associated water as well as the enhanced mass. During particularly polluted conditions, increases in aerosol optical depth of 50–100% were estimated to occur due to the observed increase in secondary aerosol mass and associated water uptake. Furthermore, the single scattering albedo was also shown to increase with height in the boundary layer. These enhancements combined to increase the negative direct aerosol radiative forcing by close to a factor of two at the median percentile level. Such increases have major ramifications for regional climate predictions as semi-volatile components are often not included in aerosol models. The results presented here provide an ideal opportunity to test regional and global representations of both the aerosol vertical distribution and subsequent impacts in North-Western Europe. North-Western Europe can be viewed as an analogue for the possible future air quality over other polluted regions of the Northern Hemisphere, where substantial reductions in sulphur dioxide emissions have yet to occur. Anticipated reductions in sulphur dioxide in polluted regions will result in an increase in the availability of ammonia to form ammonium nitrate as opposed to ammonium sulphate. This will be most important where intensive agricultural practises occur. Our observations over North-Western Europe, a region where sulphur dioxide emissions have already been reduced, indicate that failure to include the semi-volatile behaviour of ammonium nitrate will result in significant errors in predicted aerosol direct radiative forcing. Such errors will be particularly significant on regional scales.

Predicting the effects of nitrogen and planting density on maize water use in semi-arid Kenya

Models are important tools to assess the scope of management effects on crop productivity under different climatic and soil regimes. Accordingly, this study developed and used a simple model to assess the effects of nitrogen fertiliser and planting density on the water use efficiency (q) of maize in semi-arid Kenya. Field experiments were undertaken at Sonning, Berkshire, UK, in 1996 (one sowing) and 1997 (two sowings). The results from the field experiments plus soil and weather data for Machakos, Kenya (1 degree 33'S, 37 degree 14'E and 1560 m above sea level), were then used to predict the effects that N application and planting density may have on water use by a maize crop grown in semi-arid Kenya. The increase in q due to N application was greater under irrigated (15%-19%) than rainfed (7%-8%) conditions. Also, high planting density increased q (by 13%) under irrigation but decreased q (by 17%) under rainfed conditions. The current study has shown the significance of crop modelling techniques in assessing the influence of N and planting density on maize production in one region of semi-arid Kenya where there is high variability of rainfall.

Simple linear asymmetrical complexes of silver(I): NC-Ag-NH3 and Br-Ag-NH3

The compounds Ag(CN)(NH3) and Ag(Br)(NH3) are remarkable in that they form solids containing the simple molecular units NC-Ag-NH3 and Br-Ag-NH3, rather than extended solids, and are the first examples of simple linear asymmetric complexes of silver(I).

Comparing tube models for predicting the linear rheology of branched polymer melts

The hierarchical and "bob" (or branch-on-branch) models are tube-based computational models recently developed for predicting the linear rheology of general mixtures of polydisperse branched polymers. These two models are based on a similar tube-theory framework but differ in their numerical implementation and details of relaxation mechanisms. We present a detailed overview of the similarities and differences of these models and examine the effects of these differences on the predictions of the linear viscoelastic properties of a set of representative branched polymer samples in order to give a general picture of the performance of these models. Our analysis confirms that the hierarchical and bob models quantitatively predict the linear rheology of a wide range of branched polymer melts but also indicate that there is still no unique solution to cover all types of branched polymers without case-by-case adjustment of parameters such as the dilution exponent alpha and the factor p(2) which defines the hopping distance of a branch point relative to the tube diameter. An updated version of the hierarchical model, which shows improved computational efficiency and refined relaxation mechanisms, is introduced and used in these analyses.

Identification of linear systems in the presence of piecewise polynomial disturbances

The paper proposes a method of performing system identification of a linear system in the presence of bounded disturbances. The disturbances may be piecewise parabolic or periodic functions. The method is demonstrated effectively on two example systems with a range of disturbances.

Neurofuzzy control using Kalman filtering state feedback with coloured noise for unknown non-linear processes

This paper presents a controller design scheme for a priori unknown non-linear dynamical processes that are identified via an operating point neurofuzzy system from process data. Based on a neurofuzzy design and model construction algorithm (NeuDec) for a non-linear dynamical process, a neurofuzzy state-space model of controllable form is initially constructed. The control scheme based on closed-loop pole assignment is then utilized to ensure the time invariance and linearization of the state equations so that the system stability can be guaranteed under some mild assumptions, even in the presence of modelling error. The proposed approach requires a known state vector for the application of pole assignment state feedback. For this purpose, a generalized Kalman filtering algorithm with coloured noise is developed on the basis of the neurofuzzy state-space model to obtain an optimal state vector estimation. The derived controller is applied in typical output tracking problems by minimizing the tracking error. Simulation examples are included to demonstrate the operation and effectiveness of the new approach.

Structure adaptation of linear MMSE adaptive filters

The tap-length, or the number of the taps, is an important structural parameter of the linear MMSE adaptive filter. Although the optimum tap-length that balances performance and complexity varies with scenarios, most current adaptive filters fix the tap-length at some compromise value, making them inefficient to implement especially in time-varying scenarios. A novel gradient search based variable tap-length algorithm is proposed, using the concept of the pseudo-fractional tap-length, and it is shown that the new algorithm can converge to the optimum tap-length in the mean. Results of computer simulations are also provided to verify the analysis.

Linear-wavelet models applied to the identification of a two-link manipulator

This paper shows that a wavelet network and a linear term can be advantageously combined for the purpose of non linear system identification. The theoretical foundation of this approach is laid by proving that radial wavelets are orthogonal to linear functions. A constructive procedure for building such nonlinear regression structures, termed linear-wavelet models, is described. For illustration, sim ulation data are used to identify a model for a two-link robotic manipulator. The results show that the introduction of wavelets does improve the prediction ability of a linear model.

Enhanced second order non-linear optical properties of silica sol-gel thin films through photo-induced poling

Sol-gel derived inorganic materials are of interest as hosts for non-linear optically active guest molecules and they offer particular advantages in the field of non-linear optics. Orientationally ordered glasses have been prepared using a sol-gel system based on tetramethoxysilane, methyltrimethoxysilane and a non-linear optical chromophore Disperse Red 1. The novel technique of photo-induced poling was used to generate enhanced levels of polar order. The level of enhancement is strongly dependent on the extent of gelation and an optimum preparation time of ∼100 h led to an enhancement factor of ∼5. Films prepared in this manner exhibited a high stability of the polar order.

Relevance of sampling schemes in light of Ruelle's linear response theory

We reconsider the theory of the linear response of non-equilibrium steady states to perturbations. We �rst show that by using a general functional decomposition for space-time dependent forcings, we can de�ne elementary susceptibilities that allow to construct the response of the system to general perturbations. Starting from the de�nition of SRB measure, we then study the consequence of taking di�erent sampling schemes for analysing the response of the system. We show that only a speci�c choice of the time horizon for evaluating the response of the system to a general time-dependent perturbation allows to obtain the formula �rst presented by Ruelle. We also discuss the special case of periodic perturbations, showing that when they are taken into consideration the sampling can be �ne-tuned to make the de�nition of the correct time horizon immaterial. Finally, we discuss the implications of our results in terms of strategies for analyzing the outputs of numerical experiments by providing a critical review of a formula proposed by Reick.

Regularization of Descriptor Systems by Derivative and Proportional State Feedback

For linear multivariable time-invariant continuous or discrete-time singular systems it is customary to use a proportional feedback control in order to achieve a desired closed loop behaviour. Derivative feedback is rarely considered. This paper examines how derivative feedback in descriptor systems can be used to alter the structure of the system pencil under various controllability conditions. It is shown that derivative and proportional feedback controls can be constructed such that the closed loop system has a given form and is also regular and has index at most 1. This property ensures the solvability of the resulting system of dynamic-algebraic equations. The construction procedures used to establish the theory are based only on orthogonal matrix decompositions and can therefore be implemented in a numerically stable way. The problem of pole placement with derivative feedback alone and in combination with proportional state feedback is also investigated. A computational algorithm for improving the “conditioning” of the regularized closed loop system is derived.