Growth of the lava dome and extrusion rates at Soufrière Hills Volcano, Montserrat, West Indies: 2005–2008

The third episode of lava dome growth at Soufrière Hills Volcano began 1 August 2005 and ended 20 April 2007. Volumes of the dome and talus produced were measured using a photo-based method with a calibrated camera for increased accuracy. The total dense rock equivalent (DRE) volume of extruded andesite magma (306 ± 51 Mm3) was similar within error to that produced in the earlier episodes but the average extrusion rate was 5.6 ± 0.9 m3s−1 (DRE), higher than the previous episodes. Extrusion rates varied in a pulsatory manner from <0.5 m3s−1 to ∼20 m3s−1. On 18 May 2006, the lava dome had reached a volume of 85 Mm3 DRE and it was removed in its entirety during a massive dome collapse on 20 May 2006. Extrusion began again almost immediately and built a dome of 170 Mm3 DRE with a summit height 1047 m above sea level by 4 April 2007. There were few moderate-sized dome collapses (1–10 Mm3) during this extrusive episode in contrast to the first episode of dome growth in 1995–8 when they were numerous. The first and third episodes of dome growth showed a similar pattern of low (<0.5 m3s−1) but increasing magma flux during the early stages, with steady high flux after extrusion of ∼25 Mm3

Rapid model quality assessment for protein structure predictions using the comparison of multiple models without structural alignments

MOTIVATION: The accurate prediction of the quality of 3D models is a key component of successful protein tertiary structure prediction methods. Currently, clustering or consensus based Model Quality Assessment Programs (MQAPs) are the most accurate methods for predicting 3D model quality; however they are often CPU intensive as they carry out multiple structural alignments in order to compare numerous models. In this study, we describe ModFOLDclustQ - a novel MQAP that compares 3D models of proteins without the need for CPU intensive structural alignments by utilising the Q measure for model comparisons. The ModFOLDclustQ method is benchmarked against the top established methods in terms of both accuracy and speed. In addition, the ModFOLDclustQ scores are combined with those from our older ModFOLDclust method to form a new method, ModFOLDclust2, that aims to provide increased prediction accuracy with negligible computational overhead. RESULTS: The ModFOLDclustQ method is competitive with leading clustering based MQAPs for the prediction of global model quality, yet it is up to 150 times faster than the previous version of the ModFOLDclust method at comparing models of small proteins (<60 residues) and over 5 times faster at comparing models of large proteins (>800 residues). Furthermore, a significant improvement in accuracy can be gained over the previous clustering based MQAPs by combining the scores from ModFOLDclustQ and ModFOLDclust to form the new ModFOLDclust2 method, with little impact on the overall time taken for each prediction. AVAILABILITY: The ModFOLDclustQ and ModFOLDclust2 methods are available to download from: http://www.reading.ac.uk/bioinf/downloads/ CONTACT: l.j.mcguffin@reading.ac.uk.

Calculations of rovibrational energies and dipole transition intensities for polyatomic molecules using MULTIMODE

We report rigorous calculations of rovibrational energies and dipole transition intensities for three molecules using a new version of the code MULTIMODE. The key features of this code which permit, for the first time, such calculations for moderately sized but otherwise general polyatomic molecules are briefly described. Calculations for the triatomic molecule BF(2) are done to validate the code. New calculations for H(2)CO and H(2)CS are reported; these make use of semiempirical potentials but ab initio dipole moment surfaces. The new dipole surface for H(2)CO is a full-dimensional fit to the dipole moment obtained with the coupled-cluster with single and double excitations and a perturbative treatment of triple excitations method with the augmented correlation consistent triple zeta basis set. Detailed comparisons are made with experimental results from a fit to relative data for H(2)CS and absolute intensities from the HITRAN database for H(2)CO.

Rovibrational energy levels for the electronic ground state of A1OH

The vibrational-rotational energy levels of aluminum monohydroxide in its electronic ground state, (A) over tilde (1)A' AlOH, have been predicted using the variational method. The potential energy surface of the (X) over tilde (1)A' ground state of AlOH was determined employing the ab initio coupled cluster method with single, double, and perturbative triple excitations [CCSD(T)] and the correlation-consistent polarized valence quadruple zeta (cc-pVQZ) basis set. Low-lying J= 0 and J= 1 vibrational levels are reported. These are analyzed in terms of the quasilinearity of the molecule. Coriolis effects are shown to be significant. We hope that our predictions will be of value in the future when assigning rovibrational transitions in spectroscopic studies. (c) 2006 Elsevier B.V. All rights reserved.

UV absorption spectra of methyl-substituted hydroxy-cyclohexadienyl radicals in the gas phase

UV absorption spectra of five methyl-substituted hydroxy-cyclohexadienyl radicals, formed by the addition of the hydroxyl radical (OH) to toluene (methyl benzene), o-, m- and p-xylene (1,2-, 1,3- and 1,4-dimethyl benzene, respectively) and mesitylene (1,3,5-trimethylbenzene), have been determined at 298 K, 1 atm pressure (N-2 + O-2), and the corresponding absolute absorption cross-sections measured, using laser flash photolysis and time-resolved UV absorption detection. As observed for other cyclohexadienyl-type radicals, a strong absorption band is present in the 260-340 nm spectral region, with maximum cross-sections in the range (0.9-2.2) x 10(-17) cm(2) molecule(-1). The shape of the band varies significantly from one radical to the next for the series of aromatic precursors investigated. The nature and yields of hydroxylated ring-retaining oxidation products, identified in previous studies of the OH-initiated oxidation of aromatic hydrocarbons, and the results of theoretical density functional theory (DFT) calculations indicate that one or more possible isomers of the various OH-adducts may contribute to the observed spectra. Isomers where the OH-group is ortho- (or both ortho- and ipso-) to a substituent methyl-group are likely to be the most abundant but other isomers may also be formed to a significant extent. Nonetheless, the present study provides absorption spectra of the adduct radicals formed from the gas phase addition of OH to the aromatic hydrocarbons considered, near room temperature and I atm pressure. (c) 2005 Elsevier B.V. All rights reserved.

Redox chemistry and electronic properties of 2,3,5,6-tetrakis(2-pyridyl)pyrazine-bridged diruthenium complexes controlled by N,C,N '-biscyclometalated ligands

To investigate the consequences of cyclometalation for electronic communication in dinuclear ruthenium complexes, a series of 2,3,5,6-tetrakis(2-pyridyl)pyrazine (tppz) bridged diruthenium complexes was prepared and studied. These complexes have a central tppz ligand bridging via nitrogen-to-ruthenium coordination bonds, while each ruthenium atom also binds either a monoanionic, N,C,N'-terdentate 2,6-bis(2'-pyridyl)phenyl (R-N boolean AND C boolean AND N) ligand or a 2,2':6',2 ''-terpyridine (tpy) ligand. The N,C,N'-, that is, biscyclometalation, instead of the latter N,N', N ''-bonding motif significantly changes the electronic properties of the resulting complexes. Starting from well-known [{Ru(tpy)}(2)(mu-tppz)](4+) (tpy = 2,2':2 '',6-terpyridine) ([3](4+)) as a model compound, the complexes [{Ru(R-N boolean AND C boolean AND N)}(mu-tppz){Ru(tpy)}](3+) (R-N boolean AND C(H)boolean AND N = 4-R-1,3-dipyridylbenzene, R = H ([4a](3+)), CO2Me ([4b](3+))), and [{Ru(R-N boolean AND C boolean AND N)}(2)(mu-tppz)](2+), (R = H ([5a](2+)), CO2Me ([5b](2+))) were prepared with one or two N,C,N'-cyclometalated terminal ligands. The oxidation and reduction potentials of cyclometalated [4](3+) and [5](2+) are shifted negatively compared to non-cyclometalated [3](4+), the oxidation processes being affected more significantly. Compared to [3](4+), the electronic spectra of [5](2+) display large bathochromic shifts of the main MLCT transitions in the visible spectral region with low-energy absorptions tailing down to the NIR region. One-electron oxidation of [3](4+) and [5](2+) gives rise to low-energy absorption bands. The comproportionation constants and NIR band shape correspond to delocalized Robin-Day class III compounds. Complexes [4a](3+) (R = H) and [4b](3+) (R = CO2Me) also exhibit strong electronic communication, and notwithstanding the large redox-asymmetry the visible metal-to-ligand charge-transfer absorption is assigned to originate from both metal centers. The potential of the first, ruthenium-based, reversible oxidation process is strongly negatively shifted. On the contrary, the second oxidation is irreversible and cyclometalated ligand-based. Upon one-electron oxidation, a weak and low-energy absorption arises.

Characterization and classification of tumor lesions using computerized fractal-based texture analysis and support vector machines in digital mammograms

Objective: This paper presents a detailed study of fractal-based methods for texture characterization of mammographic mass lesions and architectural distortion. The purpose of this study is to explore the use of fractal and lacunarity analysis for the characterization and classification of both tumor lesions and normal breast parenchyma in mammography. Materials and methods: We conducted comparative evaluations of five popular fractal dimension estimation methods for the characterization of the texture of mass lesions and architectural distortion. We applied the concept of lacunarity to the description of the spatial distribution of the pixel intensities in mammographic images. These methods were tested with a set of 57 breast masses and 60 normal breast parenchyma (dataset1), and with another set of 19 architectural distortions and 41 normal breast parenchyma (dataset2). Support vector machines (SVM) were used as a pattern classification method for tumor classification. Results: Experimental results showed that the fractal dimension of region of interest (ROIs) depicting mass lesions and architectural distortion was statistically significantly lower than that of normal breast parenchyma for all five methods. Receiver operating characteristic (ROC) analysis showed that fractional Brownian motion (FBM) method generated the highest area under ROC curve (A z = 0.839 for dataset1, 0.828 for dataset2, respectively) among five methods for both datasets. Lacunarity analysis showed that the ROIs depicting mass lesions and architectural distortion had higher lacunarities than those of ROIs depicting normal breast parenchyma. The combination of FBM fractal dimension and lacunarity yielded the highest A z value (0.903 and 0.875, respectively) than those based on single feature alone for both given datasets. The application of the SVM improved the performance of the fractal-based features in differentiating tumor lesions from normal breast parenchyma by generating higher A z value. Conclusion: FBM texture model is the most appropriate model for characterizing mammographic images due to self-affinity assumption of the method being a better approximation. Lacunarity is an effective counterpart measure of the fractal dimension in texture feature extraction in mammographic images. The classification results obtained in this work suggest that the SVM is an effective method with great potential for classification in mammographic image analysis.

A forward-constrained regression algorithm for sparse kernel density estimation

Using the classical Parzen window (PW) estimate as the target function, the sparse kernel density estimator is constructed in a forward-constrained regression (FCR) manner. The proposed algorithm selects significant kernels one at a time, while the leave-one-out (LOO) test score is minimized subject to a simple positivity constraint in each forward stage. The model parameter estimation in each forward stage is simply the solution of jackknife parameter estimator for a single parameter, subject to the same positivity constraint check. For each selected kernels, the associated kernel width is updated via the Gauss-Newton method with the model parameter estimate fixed. The proposed approach is simple to implement and the associated computational cost is very low. Numerical examples are employed to demonstrate the efficacy of the proposed approach.

Numerical modelling of two-way propagation of non-linear dispersive waves

We describe and implement a fully discrete spectral method for the numerical solution of a class of non-linear, dispersive systems of Boussinesq type, modelling two-way propagation of long water waves of small amplitude in a channel. For three particular systems, we investigate properties of the numerically computed solutions; in particular we study the generation and interaction of approximate solitary waves.

High frequency scattering by convex curvilinear polygons

We consider the scattering of a time-harmonic acoustic incident plane wave by a sound soft convex curvilinear polygon with Lipschitz boundary. For standard boundary or finite element methods, with a piecewise polynomial approximation space, the number of degrees of freedom required to achieve a prescribed level of accuracy grows at least linearly with respect to the frequency of the incident wave. Here we propose a novel Galerkin boundary element method with a hybrid approximation space, consisting of the products of plane wave basis functions with piecewise polynomials supported on several overlapping meshes; a uniform mesh on illuminated sides, and graded meshes refined towards the corners of the polygon on illuminated and shadow sides. Numerical experiments suggest that the number of degrees of freedom required to achieve a prescribed level of accuracy need only grow logarithmically as the frequency of the incident wave increases.

Advancements in the estimation of ice particle fall speeds using laboratory and field measurements

Accurate estimates for the fall speed of natural hydrometeors are vital if their evolution in clouds is to be understood quantitatively. In this study, laboratory measurements of the terminal velocity vt for a variety of ice particle models settling in viscous fluids, along with wind-tunnel and field measurements of ice particles settling in air, have been analyzed and compared to common methods of computing vt from the literature. It is observed that while these methods work well for a number of particle types, they fail for particles with open geometries, specifically those particles for which the area ratio Ar is small (Ar is defined as the area of the particle projected normal to the flow divided by the area of a circumscribing disc). In particular, the fall speeds of stellar and dendritic crystals, needles, open bullet rosettes, and low-density aggregates are all overestimated. These particle types are important in many cloud types: aggregates in particular often dominate snow precipitation at the ground and vertically pointing Doppler radar measurements. Based on the laboratory data, a simple modification to previous computational methods is proposed, based on the area ratio. This new method collapses the available drag data onto an approximately universal curve, and the resulting errors in the computed fall speeds relative to the tank data are less than 25% in all cases. Comparison with the (much more scattered) measurements of ice particles falling in air show strong support for this new method, with the area ratio bias apparently eliminated.

Radiative efficiencies and global warming potentials using theoretically determined absorption cross-sections for several hydrofluoroethers (HFEs) and hydrofluoropolyethers (HFPEs)

Integrated infrared cross-sections and wavenumber positions for the vibrational modes of a range of hydrofluoroethers (HFEs) and hydrofluoropolyethers (HFPEs) have been calculated. Spectra were determined using a density functional method with an empirically derived correction for the wavenumbers of band positions. Radiative efficiencies (REs) were determined using the Pinnock et al. method and were used with atmospheric lifetimes from the literature to determine global warming potentials (GWPs). For the HFEs and the majority of the molecules in the HG series HFPEs, theoretically determined absorption cross-sections and REs lie within ca. 10% of those determined using measured spectra. For the larger molecules in the HG series and the HG′ series of HFPEs, agreement is less good, with theoretical values for the integrated cross-sections being up to 35% higher than the experimental values; REs are up to 45% higher. Our method gives better results than previous theoretical approaches, because of the level of theory chosen and, for REs, because an empirical wavenumber correction derived for perfluorocarbons is effective in predicting the positions of C–F stretching frequencies at around 1250 cm−1 for the molecules considered here.

A survey of gradual solar energetic particle events

We develop a database of 110 gradual solar energetic particle (SEP) events, over the period 1967–2006, providing estimates of event onset, duration, fluence, and peak flux for protons of energy E > 60 MeV. The database is established mainly from the energetic proton flux data distributed in the OMNI 2 data set; however, we also utilize the McMurdo neutron monitor and the energetic proton flux from GOES missions. To aid the development of the gradual SEP database, we establish a method with which the homogeneity of the energetic proton flux record is improved. A comparison between other SEP databases and the database developed here is presented which discusses the different algorithms used to define an event. Furthermore, we investigate the variation of gradual SEP occurrence and fluence with solar cycle phase, sunspot number (SSN), and interplanetary magnetic field intensity (Bmag) over solar cycles 20–23. We find that the occurrence and fluence of SEP events vary with the solar cycle phase. Correspondingly, we find a positive correlation between SEP occurrence and solar activity as determined by SSN and Bmag, while the mean fluence in individual events decreases with the same measures of solar activity. Therefore, although the number of events decreases when solar activity is low, the events that do occur at such times have higher fluence. Thus, large events such as the “Carrington flare” may be more likely at lower levels of solar activity. These results are discussed in the context of other similar investigations.

Benchmarking and valuation issues in measuring depreciation for European office markets

Purpose – The paper addresses the practical problems which emerge when attempting to apply longitudinal approaches to the assessment of property depreciation using valuation-based data. These problems relate to inconsistent valuation regimes and the difficulties in finding appropriate benchmarks. Design/methodology/approach – The paper adopts a case study of seven major office locations around Europe and attempts to determine ten-year rental value depreciation rates based on a longitudinal approach using IPD, CBRE and BNP Paribas datasets. Findings – The depreciation rates range from a 5 per cent PA depreciation rate in Frankfurt to a 2 per cent appreciation rate in Stockholm. The results are discussed in the context of the difficulties in applying this method with inconsistent data. Research limitations/implications – The paper has methodological implications for measuring property investment depreciation and provides an example of the problems in adopting theoretically sound approaches with inconsistent information. Practical implications – Valuations play an important role in performance measurement and cross border investment decision making and, therefore, knowledge of inconsistency of valuation practice aids decision making and informs any application of valuation-based data in the attainment of depreciation rates. Originality/value – The paper provides new insights into the use of property market valuation data in a cross-border context, insights that previously had been anecdotal and unproven in nature.

Assessing hepatic metabolic changes during progressive colonization of germ-free mouse by 1H NMR spectroscopy

It is well known that gut bacteria contribute significantly to the host homeostasis, providing a range of benefits such as immune protection and vitamin synthesis. They also supply the host with a considerable amount of nutrients, making this ecosystem an essential metabolic organ. In the context of increasing evidence of the link between the gut flora and the metabolic syndrome, understanding the metabolic interaction between the host and its gut microbiota is becoming an important challenge of modern biology.1-4 Colonization (also referred to as normalization process) designates the establishment of micro-organisms in a former germ-free animal. While it is a natural process occurring at birth, it is also used in adult germ-free animals to control the gut floral ecosystem and further determine its impact on the host metabolism. A common procedure to control the colonization process is to use the gavage method with a single or a mixture of micro-organisms. This method results in a very quick colonization and presents the disadvantage of being extremely stressful5. It is therefore useful to minimize the stress and to obtain a slower colonization process to observe gradually the impact of bacterial establishment on the host metabolism. In this manuscript, we describe a procedure to assess the modification of hepatic metabolism during a gradual colonization process using a non-destructive metabolic profiling technique. We propose to monitor gut microbial colonization by assessing the gut microbial metabolic activity reflected by the urinary excretion of microbial co-metabolites by 1H NMR-based metabolic profiling. This allows an appreciation of the stability of gut microbial activity beyond the stable establishment of the gut microbial ecosystem usually assessed by monitoring fecal bacteria by DGGE (denaturing gradient gel electrophoresis).6 The colonization takes place in a conventional open environment and is initiated by a dirty litter soiled by conventional animals, which will serve as controls. Rodents being coprophagous animals, this ensures a homogenous colonization as previously described.7 Hepatic metabolic profiling is measured directly from an intact liver biopsy using 1H High Resolution Magic Angle Spinning NMR spectroscopy. This semi-quantitative technique offers a quick way to assess, without damaging the cell structure, the major metabolites such as triglycerides, glucose and glycogen in order to further estimate the complex interaction between the colonization process and the hepatic metabolism7-10. This method can also be applied to any tissue biopsy11,12.