Very large inverse problems in atmosphere and ocean modelling

For the very large nonlinear dynamical systems that arise in a wide range of physical, biological and environmental problems, the data needed to initialize a numerical forecasting model are seldom available. To generate accurate estimates of the expected states of the system, both current and future, the technique of ‘data assimilation’ is used to combine the numerical model predictions with observations of the system measured over time. Assimilation of data is an inverse problem that for very large-scale systems is generally ill-posed. In four-dimensional variational assimilation schemes, the dynamical model equations provide constraints that act to spread information into data sparse regions, enabling the state of the system to be reconstructed accurately. The mechanism for this is not well understood. Singular value decomposition techniques are applied here to the observability matrix of the system in order to analyse the critical features in this process. Simplified models are used to demonstrate how information is propagated from observed regions into unobserved areas. The impact of the size of the observational noise and the temporal position of the observations is examined. The best signal-to-noise ratio needed to extract the most information from the observations is estimated using Tikhonov regularization theory. Copyright © 2005 John Wiley & Sons, Ltd.

Empirical mode decomposition: a novel technique for the study of tremor time series

Tremor is a clinical feature characterized by oscillations of a part of the body. The detection and study of tremor is an important step in investigations seeking to explain underlying control strategies of the central nervous system under natural (or physiological) and pathological conditions. It is well established that tremorous activity is composed of deterministic and stochastic components. For this reason, the use of digital signal processing techniques (DSP) which take into account the nonlinearity and nonstationarity of such signals may bring new information into the signal analysis which is often obscured by traditional linear techniques (e.g. Fourier analysis). In this context, this paper introduces the application of the empirical mode decomposition (EMD) and Hilbert spectrum (HS), which are relatively new DSP techniques for the analysis of nonlinear and nonstationary time-series, for the study of tremor. Our results, obtained from the analysis of experimental signals collected from 31 patients with different neurological conditions, showed that the EMD could automatically decompose acquired signals into basic components, called intrinsic mode functions (IMFs), representing tremorous and voluntary activity. The identification of a physical meaning for IMFs in the context of tremor analysis suggests an alternative and new way of detecting tremorous activity. These results may be relevant for those applications requiring automatic detection of tremor. Furthermore, the energy of IMFs was visualized as a function of time and frequency by means of the HS. This analysis showed that the variation of energy of tremorous and voluntary activity could be distinguished and characterized on the HS. Such results may be relevant for those applications aiming to identify neurological disorders. In general, both the HS and EMD demonstrated to be very useful to perform objective analysis of any kind of tremor and can therefore be potentially used to perform functional assessment.

A step toward the wet surface chemistry of glycine and alanine on Cu{110}: destabilization and decomposition in the presence of near-ambient water vapor

The coadsorption of water with organic molecules under near-ambient pressure and temperature conditions opens up new reaction pathways on model catalyst surfaces that are not accessible in conventional ultrahigh-vacuum surfacescience experiments. The surface chemistry of glycine and alanine at the water-exposed Cu{110} interface was studied in situ using ambient-pressure photoemission and X-ray absorption spectroscopy techniques. At water pressures above 10-5 Torr a significant pressure-dependent decrease in the temperature for dissociative desorption was observed for both amino acids, accompanied by the appearance of a newCN intermediate, which is not observed for lower pressures. The most likely reaction mechanisms involve dehydrogenation induced by O and/or OH surface species resulting from the dissociative adsorption of water. The linear relationship between the inverse decomposition temperature and the logarithm of water pressure enables determination of the activation energy for the surface reaction, between 213 and 232 kJ/mol, and a prediction of the decomposition temperature at the solidliquid interface by extrapolating toward the equilibrium vapor pressure. Such experiments near the equilibrium vapor pressure provide important information about elementary surface processes at the solidliquid interface, which can be retrieved neither under ultrahigh vacuum conditions nor from interfaces immersed in a solution.

The study of chiral adsorption systems using synchrotron- based structural and spectroscopic techniques: stereospecific adsorption of serine on Au-modified chiral Cu{531} surfaces

We apply modern synchrotron-based structural techniques to the study of serine adsorbed on the pure andAumodified intrinsically chiral Cu{531} surface. XPS and NEXAFS data in combination with DFT show that on the pure surface both enantiomers adsorb in l4 geometries (with de-protonated b-OH groups) at low coverage and in l3 geometries at saturation coverage. Significantly larger enantiomeric differences are seen for the l4 geometries, which involve substrate bonds of three side groups of the chiral center, i.e. a three-point interaction. The l3 adsorption geometry, where only the carboxylate and amino groups form substrate bonds, leads to smaller but still significant enantiomeric differences, both in geometry and the decomposition behavior. When Cu{531} is modified by the deposition of 1 and 2ML Au the orientations of serine at saturation coverage are significantly different from those on the clean surface. In all cases, however, a l3 bond coordination is found at saturation involving different numbers of Au atoms, which leads to relatively small enantiomeric differences.

FAMOUS, faster: using parallel computing techniques to accelerate the FAMOUS/HadCM3 climate model with a focus on the radiative transfer algorithm

We have optimised the atmospheric radiation algorithm of the FAMOUS climate model on several hardware platforms. The optimisation involved translating the Fortran code to C and restructuring the algorithm around the computation of a single air column. Instead of the existing MPI-based domain decomposition, we used a task queue and a thread pool to schedule the computation of individual columns on the available processors. Finally, four air columns are packed together in a single data structure and computed simultaneously using Single Instruction Multiple Data operations. The modified algorithm runs more than 50 times faster on the CELL’s Synergistic Processing Elements than on its main PowerPC processing element. On Intel-compatible processors, the new radiation code runs 4 times faster. On the tested graphics processor, using OpenCL, we find a speed-up of more than 2.5 times as compared to the original code on the main CPU. Because the radiation code takes more than 60% of the total CPU time, FAMOUS executes more than twice as fast. Our version of the algorithm returns bit-wise identical results, which demonstrates the robustness of our approach. We estimate that this project required around two and a half man-years of work.

Decomposition of surface electromyographic signal using Hidden Markov Model

The detection of physiological signals from the motor system (electromyographic signals) is being utilized in the practice clinic to guide the therapist in a more precise and accurate diagnosis of motor disorders. In this context, the process of decomposition of EMG (electromyographic) signals that includes the identification and classification of MUAP (Motor Unit Action Potential) of a EMG signal, is very important to help the therapist in the evaluation of motor disorders. The EMG decomposition is a complex task due to EMG features depend on the electrode type (needle or surface), its placement related to the muscle, the contraction level and the health of the Neuromuscular System. To date, the majority of researches on EMG decomposition utilize EMG signals acquired by needle electrodes, due to their advantages in processing this type of signal. However, relatively few researches have been conducted using surface EMG signals. Thus, this article aims to contribute to the clinical practice by presenting a technique that permit the decomposition of surface EMG signal via the use of Hidden Markov Models. This process is supported by the use of differential evolution and spectral clustering techniques. The developed system presented coherent results in: (1) identification of the number of Motor Units actives in the EMG signal; (2) presentation of the morphological patterns of MUAPs in the EMG signal; (3) identification of the firing sequence of the Motor Units. The model proposed in this work is an advance in the research area of decomposition of surface EMG signals.

Adsorption, orientation and thermal decomposition of copper(II) hexafluoroacetylacetonate on rutile TiO2(110)

We have investigated the adsorption and thermal decomposition of copper hexafluoroacetylacetonate (Cu-11(hfaC)(2)) on single crystal rutile TiO2(110). Low energy electron diffraction shows that room temperature saturation coverage of the Cu-II(hfac)(2) adsorbate forms an ordered (2 x 1) over-layer. X-ray and ultra-violet photoemission spectroscopy of the saturated surface were recorded as the sample was annealed in a sequential manner to reveal decomposition pathways. The results show that the molecule dissociatively adsorbs by detachment of one of the two ligands to form hfac and Cu-1(hfac) which chemisorb to the substrate at 298 K. These ligands only begin to decompose once the surface temperature exceeds 473 K where Cu core level shifts indicate metallisation. This reduction from Cu(I) to Cu(0) takes place in the absence of an external reducing agent and without disproportionation and is accompanied by the onset of decomposition of the hfac ligands. Finally, C K-edge near edge X-ray absorption fine structure experiments indicate that both the ligands adsorb aligned in the < 001 > direction and we propose a model in which the hfac ligands adsorb on the 5-fold coordinated Ti atoms and the Cu-1(hfac) moiety attaches to the bridging O atoms in a square planar geometry. The calculated tilt angle for these combined geometries is approximately 10 degrees to the surface normal.

Decomposition in soil and chemical characteristics of pollen

The input to soils made by pollen and its subsequent mineralization has rarely been investigated from a soil microbiological point of view even though the small but significant quantities of C and N in pollen may make an important contribution to nutrient cycling. The relative resistance to decomposition of pollen exines (outer layers) has led to much of the focus of pollen in soil being on its preservation for archaeological and palaeo-ecological purposes. We have examined aspects of the chemical composition and decomposition of pollen from birch (Betula alba) and maize (Zea mays) in soil. The relatively large N contents, small C-to-N ratios and large water-soluble contents of pollen from both species indicated that they would be readily mineralized in soil. When added to soil and incubated at 16 degrees C an amount of C equivalent to 22-26% of the added pollen C was lost as CO2 within 22 days, with the Z. mays pollen decomposing faster. For B. alba pollen, the water-soluble fraction decomposed faster than the whole pollen and the insoluble fraction decomposed more slowly over 22 days. By contrast, there were no significant differences in the decomposition rates of the different fractions from Z. mays pollen. Solid-state C-13 nuclear magnetic resonance (NMR) revealed no gross chemical differences between the pollen of these two species, with strong resonances in the alkyl- and methyl-C region (0-45 p.p.m.) indicative of aliphatic compounds, the O-alkyl-C (60-90 p.p.m.) and the acetal- and ketal-C region (90-110 p.p.m.) indicative of polysaccharides, and the carbonyl-C region indicative of peptides and carboxylic acids. In addition, both pollens gave a small but distinct resonance at 55 p.p.m. attributed to N-alkyl-C. The resonances attributed to polysaccharides were lost completely or substantially reduced after decomposition.

Decomposition and phosphorus release from organic residues as affected by residue quality and added inorganic phosphorus

The combined use of organic residue and inorganic fertiliser-phosphorus (P) is appropriate in meeting both the short and long-term P requirement of crops. To assess the influence of added inorganic fertiliser-P on the processes of decomposition and P release from the residue and the relationships with quality, prunings of Gliricidia sepium, Leucaena leucocephela, Senna siamea, Acacia mangium and Paraserienthus falcataria were incubated without and with added inorganic fertiliser-P for 56 days. Soil was added only as inoculum. Decomposition rate and amounts of acid extractable-P (P release) were in the same order: G. sepium > S. siamea > L. leucocepheta > P falcataria > A. mangium. Unlike the other residues, A. mangium released no P despite the loss of half its mass during the 8 weeks of incubation. The residue P content correlated with P release. However, decomposition rate did not correlate with residue P content but with the lignin, polyphenol and cellulose content, and ratios to P. These ratios were negatively correlated with P release suggesting that lignin and polyphenol contents influence P release more when the residue-P content is low. Results suggest that rate of decomposition influences the release of P. The critical residue P content for P release was estimated to be 0.12% < P < 0.19%. Added P had no effect on decomposition and P release from the residues.

A review of dissolved oxygen modelling techniques for lowland rivers

This review introduces the methods used to simulate the processes affecting dissolved oxygen (DO) in lowland rivers. The important processes are described and this provides a modelling framework to describe those processes in the context of a mass-balance model. The process equations that are introduced all require (reaction) rate parameters and a variety of common procedures for identifying those parameters are reviewed. This is important because there is a wide range of estimation techniques for many of the parameters. These different techniques elicit different estimates of the parameter value and so there is the potential for a significant uncertainty in the model's inputs and therefore in the output too. Finally, the data requirements for modelling DO in lowland rivers are summarised on the basis of modelling the processes described in this review using a mass-balance model. This is reviewed with regard to what data are available and from where they might be obtained. (C) 2003 Elsevier Science B.V. All rights reserved.

The use of earth observation techniques to improve catchment-scale pollution predictions

Remote sensing can potentially provide information useful in improving pollution transport modelling in agricultural catchments. Realisation of this potential will depend on the availability of the raw data, development of information extraction techniques, and the impact of the assimilation of the derived information into models. High spatial resolution hyperspectral imagery of a farm near Hereford, UK is analysed. A technique is described to automatically identify the soil and vegetation endmembers within a field, enabling vegetation fractional cover estimation. The aerially-acquired laser altimetry is used to produce digital elevation models of the site. At the subfield scale the hypothesis that higher resolution topography will make a substantial difference to contaminant transport is tested using the AGricultural Non-Point Source (AGNPS) model. Slope aspect and direction information are extracted from the topography at different resolutions to study the effects on soil erosion, deposition, runoff and nutrient losses. Field-scale models are often used to model drainage water, nitrate and runoff/sediment loss, but the demanding input data requirements make scaling up to catchment level difficult. By determining the input range of spatial variables gathered from EO data, and comparing the response of models to the range of variation measured, the critical model inputs can be identified. Response surfaces to variation in these inputs constrain uncertainty in model predictions and are presented. Although optical earth observation analysis can provide fractional vegetation cover, cloud cover and semi-random weather patterns can hinder data acquisition in Northern Europe. A Spring and Autumn cloud cover analysis is carried out over seven UK sites close to agricultural districts, using historic satellite image metadata, climate modelling and historic ground weather observations. Results are assessed in terms of probability of acquisition probability and implications for future earth observation missions. (C) 2003 Elsevier Ltd. All rights reserved.