Stable reduced-order models for discrete-time systems

The Routh-stability method is employed to reduce the order of discrete-time system transfer functions. It is shown that the Routh approximant is well suited to reduce both the denominator and the numerator polynomials, although alternative methods, such as PadÃ�Â(c)-Markov approximation, are also used to fit the model numerator coefficients.

Periodic time series modeling of environmental proxy records with guaranteed positive growth rate estimation

Identifying a periodic time-series model from environmental records, without imposing the positivity of the growth rate, does not necessarily respect the time order of the data observations. Consequently, subsequent observations, sampled in the environmental archive, can be inversed on the time axis, resulting in a non-physical signal model. In this paper an optimization technique with linear constraints on the signal model parameters is proposed that prevents time inversions. The activation conditions for this constrained optimization are based upon the physical constraint of the growth rate, namely, that it cannot take values smaller than zero. The actual constraints are defined for polynomials and first-order splines as basis functions for the nonlinear contribution in the distance-time relationship. The method is compared with an existing method that eliminates the time inversions, and its noise sensitivity is tested by means of Monte Carlo simulations. Finally, the usefulness of the method is demonstrated on the measurements of the vessel density, in a mangrove tree, Rhizophora mucronata, and the measurement of Mg/Ca ratios, in a bivalve, Mytilus trossulus.

A stress-controlled mechanism for the intensity of very large magnitude explosive eruptions

Large magnitude explosive eruptions are the result of the rapid and large-scale transport of silicic magma stored in the Earth's crust, but the mechanics of erupting teratonnes of silicic magma remain poorly understood. Here, we demonstrate that the combined effect of local crustal extension and magma chamber overpressure can sustain linear dyke-fed explosive eruptions with mass fluxes in excess of 10^10 kg/s from shallow-seated (4–6 km depth) chambers during moderate extensional stresses. Early eruption column collapse is facilitated with eruption duration of the order of few days with an intensity of at least one order of magnitude greater than the largest eruptions in the 20th century. The conditions explored in this study are one way in which high mass eruption rates can be achieved to feed large explosive eruptions. Our results corroborate geological and volcanological evidences from volcano-tectonic complexes such as the Sierra Madre Occidental (Mexico) and the Taupo Volcanic Zone (New Zealand).

Testing for non-stationarity and cointegration allowing for the possibility of a structural break: an application to EuroSterling interest rates

In this paper we examine the order of integration of EuroSterling interest rates by employing techniques that can allow for a structural break under the null and/or alternative hypothesis of the unit-root tests. In light of these results, we investigate the cointegrating relationship implied by the single, linear expectations hypothesis of the term structure of interest rates employing two techniques, one of which allows for the possibility of a break in the mean of the cointegrating relationship. The aim of the paper is to investigate whether or not the interest rate series can be viewed as I(1) processes and furthermore, to consider whether there has been a structural break in the series. We also determine whether, if we allow for a break in the cointegration analysis, the results are consistent with those obtained when a break is not allowed for. The main results reported in this paper support the conjecture that the ‘short’ Euro-currency rates are characterised as I(1) series that exhibit a structural break on or near Black Wednesday, 16 September 1992, whereas the ‘long’ rates are I(1) series that do not support the presence of a structural break. The evidence from the cointegration analysis suggests that tests of the expectations hypothesis based on data sets that include the ERM crisis period, or a period that includes a structural break, might be problematic if the structural break is not explicitly taken into account in the testing framework.

The variability of tropical ice cloud properties as a function of the large-scale context from ground-based radar-lidar observations over Darwin, Australia

The high complexity of cloud parameterizations now held in models puts more pressure on observational studies to provide useful means to evaluate them. One approach to the problem put forth in the modelling community is to evaluate under what atmospheric conditions the parameterizations fail to simulate the cloud properties and under what conditions they do a good job. It is the ambition of this paper to characterize the variability of the statistical properties of tropical ice clouds in different tropical "regimes" recently identified in the literature to aid the development of better process-oriented parameterizations in models. For this purpose, the statistical properties of non-precipitating tropical ice clouds over Darwin, Australia are characterized using ground-based radar-lidar observations from the Atmospheric Radiation Measurement (ARM) Program. The ice cloud properties analysed are the frequency of ice cloud occurrence, the morphological properties (cloud top height and thickness), and the microphysical and radiative properties (ice water content, visible extinction, effective radius, and total concentration). The variability of these tropical ice cloud properties is then studied as a function of the large-scale cloud regimes derived from the International Satellite Cloud Climatology Project (ISCCP), the amplitude and phase of the Madden-Julian Oscillation (MJO), and the large-scale atmospheric regime as derived from a long-term record of radiosonde observations over Darwin. The vertical variability of ice cloud occurrence and microphysical properties is largest in all regimes (1.5 order of magnitude for ice water content and extinction, a factor 3 in effective radius, and three orders of magnitude in concentration, typically). 98 % of ice clouds in our dataset are characterized by either a small cloud fraction (smaller than 0.3) or a very large cloud fraction (larger than 0.9). In the ice part of the troposphere three distinct layers characterized by different statistically-dominant microphysical processes are identified. The variability of the ice cloud properties as a function of the large-scale atmospheric regime, cloud regime, and MJO phase is large, producing mean differences of up to a factor 8 in the frequency of ice cloud occurrence between large-scale atmospheric regimes and mean differences of a factor 2 typically in all microphysical properties. Finally, the diurnal cycle of the frequency of occurrence of ice clouds is also very different between regimes and MJO phases, with diurnal amplitudes of the vertically-integrated frequency of ice cloud occurrence ranging from as low as 0.2 (weak diurnal amplitude) to values in excess of 2.0 (very large diurnal amplitude). Modellers should now use these results to check if their model cloud parameterizations are capable of translating a given atmospheric forcing into the correct statistical ice cloud properties.

Amplifying vibrational circular dichroism by manipulation of the electronic manifold

Vibrational circular dichroism is a powerful technique to study the stereochemistry of chiral molecules, but often suffers from small signal intensities. Electrochemical modulation of the energies of the electronically excited state manifold is now demonstrated to lead to an order of magnitude enhancement of the differential absorption. Quantum-chemical calculations show that increased mixing between ground and excited states is at the origin of this amplification.

Evolution of network structure in polymer-stabilised liquid crystals

Experiments were performed to investigate the evolution of structure and morphology of the network in polymer-stabilised liquid crystals. In situ optical microscopy revealed that the morphology was significantly altered by extraction of the LC host, while scanning electron microscopy showed that the network morphology was also dependent on the polymerisation conditions and closely related to the depletion of monomer, as monitored by high performance liquid chromatography. Transmission electron microscopy allowed observation of internal structure, resolving microstructure on the order of 0. 1 μm.

Side-chain liquid crystalline elastomers: the relationship between the orientational ordering of the polymer backbone and the length of the coupling chain

The influence of cross-linking on the phase behaviour of a series of side-chain liquid crystalline elastomers has been studied. For samples cross-linked in the temperature range corresponding to the nematic phase, the phase transition was shifted compared to that observed when an identical sample was cross-linked in the isotropic phase. This shift represented a stabilisation of the nematic phase in the former case, in line with theoretical expectations. By utilising a novel, slow cross-linking method, which allows the polymer backbone to take up an equilibrium conformation prior to network formation, it proved possible to monitor the shifts in phase transition temperature as a function of the length of the methylene chain coupling the mesogenic units to the polymer backbone. The results obtained are related to the backbone anisotropy and indicate that the level of orientational order of the polymer in the nematic phase backbone increases with a reduction in the length of the coupling chain.

Evidence of dispersion relations for the nonlinear response of the Lorenz 63 system

Along the lines of the nonlinear response theory developed by Ruelle, in a previous paper we have proved under rather general conditions that Kramers-Kronig dispersion relations and sum rules apply for a class of susceptibilities describing at any order of perturbation the response of Axiom A non equilibrium steady state systems to weak monochromatic forcings. We present here the first evidence of the validity of these integral relations for the linear and the second harmonic response for the perturbed Lorenz 63 system, by showing that numerical simulations agree up to high degree of accuracy with the theoretical predictions. Some new theoretical results, showing how to derive asymptotic behaviors and how to obtain recursively harmonic generation susceptibilities for general observables, are also presented. Our findings confirm the conceptual validity of the nonlinear response theory, suggest that the theory can be extended for more general non equilibrium steady state systems, and shed new light on the applicability of very general tools, based only upon the principle of causality, for diagnosing the behavior of perturbed chaotic systems and reconstructing their output signals, in situations where the fluctuation-dissipation relation is not of great help.

New results on the thermodynamical properties of the climate system

In this paper the authors exploit two equivalent formulations of the average rate of material entropy production in the climate system to propose an approximate splitting between contributions due to vertical and eminently horizontal processes. This approach is based only on 2D radiative fields at the surface and at the top of atmosphere. Using 2D fields at the top of atmosphere alone, lower bounds to the rate of material entropy production and to the intensity of the Lorenz energy cycle are derived. By introducing a measure of the efficiency of the planetary system with respect to horizontal thermodynamic processes, it is possible to gain insight into a previous intuition on the possibility of defining a baroclinic heat engine extracting work from the meridional heat flux. The approximate formula of the material entropy production is verified and used for studying the global thermodynamic properties of climate models (CMs) included in the Program for Climate Model Diagnosis and Intercomparison (PCMDI)/phase 3 of the Coupled Model Intercomparison Project (CMIP3) dataset in preindustrial climate conditions. It is found that about 90% of the material entropy production is due to vertical processes such as convection, whereas the large-scale meridional heat transport contributes to only about 10% of the total. This suggests that the traditional two-box models used for providing a minimal representation of entropy production in planetary systems are not appropriate, whereas a basic—but conceptually correct—description can be framed in terms of a four-box model. The total material entropy production is typically 55 mW m−2 K−1, with discrepancies on the order of 5%, and CMs’ baroclinic efficiencies are clustered around 0.055. The lower bounds on the intensity of the Lorenz energy cycle featured by CMs are found to be around 1.0–1.5 W m−2, which implies that the derived inequality is rather stringent. When looking at the variability and covariability of the considered thermodynamic quantities, the agreement among CMs is worse, suggesting that the description of feedbacks is more uncertain. The contributions to material entropy production from vertical and horizontal processes are positively correlated, so that no compensation mechanism seems in place. Quite consistently among CMs, the variability of the efficiency of the system is a better proxy for variability of the entropy production due to horizontal processes than that of the large-scale heat flux. The possibility of providing constraints on the 3D dynamics of the fluid envelope based only on 2D observations of radiative fluxes seems promising for the observational study of planets and for testing numerical models.

Composite control of non-linear singularly perturbed systems: a geometric approach

Using a geometric approach, a composite control—the sum of a slow control and a fast control—is derived for a general class of non-linear singularly perturbed systems. A new and simpler method of composite control design is proposed whereby the fast control is completely designed at the outset. The slow control is then free to be chosen such that the slow integral manifold of the original system approximates a desired design manifold to within any specified order of ε accuracy.

Composite control of non-linear singularly perturbed systems: a geometric approach

Using a geometric approach, a composite control—the sum of a slow control and a fast control—is derived for a general class of non-linear singularly perturbed systems. A new and simpler method of composite control design is proposed whereby the fast control is completely designed at the outset. The slow control is then free to be chosen such that the slow integral manifold of the original system approximates a desired design manifold to within any specified order of ε accuracy.

Chromium nanostructures formed by dewetting of heteroepitaxial films on W(100)

In this paper, we report the surprising formation of square-based facetted islands with linear dimension of the order of 500 nm upon dewetting of a Cr multilayer onW(100).We show that these square islands are composed of inclined facets surrounding a depressed center such that the facet slopes inward with the outer edges of the islands thicker than the centers. The islands’ shapes do not represent traditional equilibrium crystal shapes as expected for a Wulf construction. In situ UV and x-ray photoelectron emission microscopy allied to spatially resolved spectroscopy throws considerable light on the nature of the dewetting and shows that the metal surface between the islands remains covered by a thin pseudomorphic wetting layer of ∼1 ML. Low-energy electron diffraction and scanning tunneling and atomic force microscopies allow quantification of facet slopes, and we identify a predominance of tilted Cr(100) facets ±5◦ off of the substrate normal bound by (210) planes at ∼26◦. The epitaxial Cr islands adopt the bulk Cr lattice constant but are tilted with respect to the surface normal.We suggest that the Cr crystallite tilting creates a vicinal-like interface structure that determines the island morphology

Influence of a non-ionic amphiphilic copolymer on the self-assembly of a peptide amphiphile that forms nanotapes

The influence of a non-ionic polymeric surfactant on the self-assembly of a peptide amphiphile (PA) that forms nanotapes is investigated using a combination of microscopic, scattering and spectroscopic techniques. Mixtures of Pluronic copolymer P123 with the PA C16-KTTKS in aqueous solution were studied at a fixed concentration of the PA at which it is known to self-assemble into extended nanotapes, but varying P123 concentration. We find that P123 can disrupt the formation of C16- KTTKS nanotapes, leading instead to cylindrical nanofibril structures. The spherical micelles formed by P123 at room temperature are disrupted in the presence of the PA. There is a loss of cloudiness in the solutions as the large nanotape aggregates formed by C16-KTTKS are broken up, by P123 solubilization. At least locally, b-sheet structure is retained, as confirmed by XRD and FTIR spectroscopy, even for solutions containing 20 wt% P123. This indicates, unexpectedly, that peptide secondary structure can be retained in solutions with high concentration of non-ionic surfactant. Selfassembly in this system exhibits slow kinetics towards equilibrium, the initial self-assembly being dependent on the order of mixing. Heating above the lipid chain melting temperature assists in disrupting trapped non-equilibrium states.

The influence of Langmuir turbulence on the scaling for the dissipation rate in the oceanic boundary layer

A model for estimating the turbulent kinetic energy dissipation rate in the oceanic boundary layer, based on insights from rapid-distortion theory, is presented and tested. This model provides a possible explanation for the very high dissipation levels found by numerous authors near the surface. It is conceived that turbulence, injected into the water by breaking waves, is subsequently amplified due to its distortion by the mean shear of the wind-induced current and straining by the Stokes drift of surface waves. The partition of the turbulent shear stress into a shear-induced part and a wave-induced part is taken into account. In this picture, dissipation enhancement results from the same mechanism responsible for Langmuir circulations. Apart from a dimensionless depth and an eddy turn-over time, the dimensionless dissipation rate depends on the wave slope and wave age, which may be encapsulated in the turbulent Langmuir number La_t. For large La_t, or any Lat but large depth, the dissipation rate tends to the usual surface layer scaling, whereas when Lat is small, it is strongly enhanced near the surface, growing asymptotically as ɛ ∝ La_t^{-2} when La_t → 0. Results from this model are compared with observations from the WAVES and SWADE data sets, assuming that this is the dominant dissipation mechanism acting in the ocean surface layer and statistical measures of the corresponding fit indicate a substantial improvement over previous theoretical models. Comparisons are also carried out against more recent measurements, showing good order-of-magnitude agreement, even when shallow-water effects are important.