Modelling the conversion of Colombian lowland ecosystems since 1940: Drivers, patterns and rates

In biologically mega-diverse countries that are undergoing rapid human landscape transformation, it is important to understand and model the patterns of land cover change. This problem is particularly acute in Colombia, where lowland forests are being rapidly cleared for cropping and ranching. We apply a conceptual model with a nested set of a priori predictions to analyse the spatial and temporal patterns of land cover change for six 50-100 km(2) case study areas in lowland ecosystems of Colombia. Our analysis included soil fertility, a cost-distance function, and neighbourhood of forest and secondary vegetation cover as independent variables. Deforestation and forest regrowth are tested using logistic regression analysis and an information criterion approach to rank the models and predictor variables. The results show that: (a) overall the process of deforestation is better predicted by the full model containing all variables, while for regrowth the model containing only the auto-correlated neighbourhood terms is a better predictor; (b) overall consistent patterns emerge, although there are variations across regions and time; and (c) during the transformation process, both the order of importance and significance of the drivers change. Forest cover follows a consistent logistic decline pattern across regions, with introduced pastures being the major replacement land cover type. Forest stabilizes at 2-10% of the original cover, with an average patch size of 15.4 (+/- 9.2) ha. We discuss the implications of the observed patterns and rates of land cover change for conservation planning in countries with high rates of deforestation. (c) 2005 Elsevier Ltd. All rights reserved.

Acid-mediated conversion of methylene-interrupted bisepoxides to tetrahydrofurans: A biomimetic transformation

The acid-mediated transformation of syn and anti methylene interrupted cis,cis and cis,trans bisepoxides to tetrahydrofurans is high yielding, and demonstrates both regioselectivity and stereoselectivity. Trans,trans methylene interrupted bisepoxides do not yield tetrahydrofurans under the same conditions.

A model for adsorption of condensable vapors on nonporous materials

A new isotherm is proposed here for adsorption of condensable vapors and gases on nonporous materials having type II isotherms according to the Brunauer-Deming-Deming-Teller (BDDH) classification. The isotherm combines the recent molecular-continuum model in the multilayer region, with other widely used models for sub-monolayer coverage, some of which satisfy the requirement of a Henry's law asymptote. The model is successfully tested using isotherm data for nitrogen adsorption on nonporous silica, carbon and alumina, as well as benzene and hexane adsorption on nonporous carbon. Based on the data fits, out of several different alternative choices of model for the monolayer region, the Freundlich and the Unilan models are found to be the most successful when combined with the multilayer model to predict the whole isotherm. The hybrid model is consequently applicable over a wide pressure range. (C) 2000 Elsevier Science B.V. All rights reserved.

Conversion of glucosamine to galactosamine and allosamine derivatives: Control of inversions of stereochemistry at C-3 and C-4

The reactions of sodium benzoate with a series of trimesylates derived from glucosamine have been examined in an attempt to gain facile access to galactosamine analogues. Trimesylate 17, in which the amino group was protected as a phthalimide, underwent double displacement at positions 4 and 6 to give the dibenzoate 18 with the desired galactosamine configuration. In contrast, trimesylates 21 and 27, in which the amino groups were protected as acetamides, unexpectedly underwent double displacement at positions 3 and 6, giving products 22 and 28, respectively, with allosamine configurations.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Model-based analysis of anaerobic acetate uptake by a mixed culture of polyphosphate-accumulating and glycogen-accumulating organisms

An increasing number of studies shows that the glycogen-accumulating organisms (GAOs) can survive and may indeed proliferate under the alternating anaerobic/aerobic conditions found in EBPR systems, thus forming a strong competitor of the polyphosphate-accumulating organisms (PAOs). Understanding their behaviors in a mixed PAO and GAO culture under various operational conditions is essential for developing operating strategies that disadvantage the growth of this group of unwanted organisms. A model-based data analysis method is developed in this paper for the study of the anaerobic PAO and GAO activities in a mixed PAO and GAO culture. The method primarily makes use of the hydrogen ion production rate and the carbon dioxide transfer rate resulting from the acetate uptake processes by PAOs and GAOs, measured with a recently developed titration and off-gas analysis (TOGA) sensor. The method is demonstrated using the data from a laboratory-scale sequencing batch reactor (SBR) operated under alternating anaerobic and aerobic conditions. The data analysis using the proposed method strongly indicates a coexistence of PAOs and GAOs in the system, which was independently confirmed by fluorescent in situ hybridization (FISH) measurement. The model-based analysis also allowed the identification of the respective acetate uptake rates by PAOs and GAOs, along with a number of kinetic and stoichiometric parameters involved in the PAO and GAO models. The excellent fit between the model predictions and the experimental data not involved in parameter identification shows that the parameter values found are reliable and accurate. It also demonstrates that the current anaerobic PAO and GAO models are able to accurately characterize the PAO/GAO mixed culture obtained in this study. This is of major importance as no pure culture of either PAOs or GAOs has been reported to date, and hence the current PAO and GAO models were developed for the interpretation of experimental results of mixed cultures. The proposed method is readily applicable for detailed investigations of the competition between PAOs and GAOs in enriched cultures. However, the fermentation of organic substrates carried out by ordinary heterotrophs needs to be accounted for when the method is applied to the study of PAO and GAO competition in full-scale sludges. (C) 2003 Wiley Periodicals, Inc.

On a resampling approach for tests on the number of clusters with mixture model-based clustering of tissue samples

We consider the problem of assessing the number of clusters in a limited number of tissue samples containing gene expressions for possibly several thousands of genes. It is proposed to use a normal mixture model-based approach to the clustering of the tissue samples. One advantage of this approach is that the question on the number of clusters in the data can be formulated in terms of a test on the smallest number of components in the mixture model compatible with the data. This test can be carried out on the basis of the likelihood ratio test statistic, using resampling to assess its null distribution. The effectiveness of this approach is demonstrated on simulated data and on some microarray datasets, as considered previously in the bioinformatics literature. (C) 2004 Elsevier Inc. All rights reserved.

Speeding up the EM algorithm for mixture model-based segmentation of magnetic resonance images

Mixture models implemented via the expectation-maximization (EM) algorithm are being increasingly used in a wide range of problems in pattern recognition such as image segmentation. However, the EM algorithm requires considerable computational time in its application to huge data sets such as a three-dimensional magnetic resonance (MR) image of over 10 million voxels. Recently, it was shown that a sparse, incremental version of the EM algorithm could improve its rate of convergence. In this paper, we show how this modified EM algorithm can be speeded up further by adopting a multiresolution kd-tree structure in performing the E-step. The proposed algorithm outperforms some other variants of the EM algorithm for segmenting MR images of the human brain. (C) 2004 Pattern Recognition Society. Published by Elsevier Ltd. All rights reserved.

Angiotensin II facilitates autoregulation in the perfused mouse kidney: An optimized in vitro model for assessment of renal vascular and tubular function

Background and Aims: We have optimized the isolated perfused mouse kidney (IPMK) model for studying renal vascular and tubular function in vitro using 24-28 g C57BL6J mice; the wild type controls for many transgenic mice. Methods and Results: Buffer composition was optimized for bovine serum albumin concentration (BSA). The effect of adding erythrocytes on renal function and morphology was assessed. Autoregulation was investigated during stepped increases in perfusion pressure. Perfusion for 60 min at 90-110 mmHg with Krebs bicarbonate buffer containing 5.5% BSA, and amino acids produced functional parameters within the in vivo range. Erythrocytes increased renal vascular resistance (3.8 +/- 0.2 vs 2.4 +/- 0.1 mL/min.mmHg, P < 0.05), enhanced sodium reabsorption (FENa = 0.3 +/- 0.08 vs 1.5 +/- 0.7%, P < 0.05), produced equivalent glomerular filtration rates (GFR; 364 +/- 38 vs 400 +/- 9 muL/min per gkw) and reduced distal tubular cell injury in the inner stripe (5.8 +/- 1.7 vs 23.7 +/- 3.1%, P < 0.001) compared to cell free perfusion. The IPMK was responsive to vasoconstrictor (angiotensin II, EC50 100 pM) and vasodilator (methacholine, EC50 75 nM) mediators and showed partial autoregulation of perfusate flow under control conditions over 65-85 mmHg; autoregulatory index (ARI) of 0.66 +/- 0.11. Angiotensin II (100 pM) extended this range (to 65-120 mmHg) and enhanced efficiency (ARI 0.21 +/- 0.02, P < 0.05). Angiotensin II facilitation was antagonized by methacholine (ARI 0.76 +/- 0.08) and papaverine (ARI 0.91 +/- 0.13). Conclusion: The IPMK model is useful for studying renal physiology and pathophysiology without systemic neurohormonal influences.

An FDTD model for calculation of gradient-induced eddy currents in MRI system

In most magnetic resonance imaging (MRI) systems, pulsed magnetic gradient fields induce eddy currents in the conducting structures of the superconducting magnet. The eddy currents induced in structures within the cryostat are particularly problematic as they are characterized by long time constants by virtue of the low resistivity of the conductors. This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) scheme in cylindrical coordinates for eddy-current calculation in conductors. This model is intended to be part of a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The singularity apparent in the governing equations is removed by using a series expansion method and the conductor-air boundary condition is handled using a variant of the surface impedance concept. The numerical difficulty due to the asymmetry of Maxwell equations for low-frequency eddy-current problems is circumvented by taking advantage of the known penetration behavior of the eddy-current fields. A perfectly matched layer absorbing boundary condition in 3-D cylindrical coordinates is also incorporated. The numerical method has been verified against analytical solutions for simple cases. Finally, the algorithm is illustrated by modeling a pulsed field gradient coil system within an MRI magnet system. The results demonstrate that the proposed FDTD scheme can be used to calculate large-scale eddy-current problems in materials with high conductivity at low frequencies.

The use of simultaneous H-1 & P-31 magic angle spinning nuclear magnetic resonance measurements to characterize energy metabolism during the conversion of muscle to meat

In order to investigate the potential of magic angle spinning nuclear magnetic resonance (MAS NMR) in the elucidation of post-mortem metabolism in muscle biopsies, simultaneous H-1 and (31)p MAS NMR measurements were made continuously on postmortem (20 min to 24 h) muscle longissimus samples from rabbits. The animals had either been or not been given adrenaline (0.5 mg kg(-1) 4 h pre-slaughter) to deplete stores of muscle glycogen. The intracellular pH was calculated from H-1 spectra, and the post-mortem rate of formation of lactate was followed and quantified. Comparison of measurements made on muscle samples from rabbits treated with adrenaline with measurements made on muscle samples from untreated' rabbits revealed significant effects of adrenaline treatment on both pH (pH24 h = 6.42 vs. pH24 It = 5.60) and formation of lactate (16 mmol g(-1) vs. 65 mmol g(-1)). The P-31 NMR spectra were used to follow the rate of degradation of ATP and phosphocreatine. The present study clearly shows that MAS NMR has potential for the study of post-mortem energy metabolism.

A Bethe ansatz solvable model for superpositions of Cooper pairs and condensed molecular bosons

We introduce a general Hamiltonian describing coherent superpositions of Cooper pairs and condensed molecular bosons. For particular choices of the coupling parameters, the model is integrable. One integrable manifold, as well as the Bethe ansatz solution, was found by Dukelsky et al. [J. Dukelsky, G.G. Dussel, C. Esebbag, S. Pittel, Phys. Rev. Lett. 93 (2004) 050403]. Here we show that there is a second integrable manifold, established using the boundary quantum inverse scattering method. In this manner we obtain the exact solution by means of the algebraic Bethe ansatz. In the case where the Cooper pair energies are degenerate we examine the relationship between the spectrum of these integrable Hamiltonians and the quasi-exactly solvable spectrum of particular Schrodinger operators. For the solution we derive here the potential of the Schrodinger operator is given in terms of hyperbolic functions. For the solution derived by Dukelsky et al., loc. cit. the potential is sextic and the wavefunctions obey PT-symmetric boundary conditions. This latter case provides a novel example of an integrable Hermitian Hamiltonian acting on a Fock space whose states map into a Hilbert space of PE-symmetric wavefunctions defined on a contour in the complex plane. (c) 2006 Elsevier B.V. All rights reserved.