Orography in a contour dynamics model of large-scale atmospheric flow

The influence of orography on the structure of stationary planetary Rossby waves is studied in the context of a contour dynamics model of the large-scale atmospheric flow. Orography of infinitesimal and finite amplitude is studied using analytical and numerical techniques. Three different types of orography are considered: idealized orography in the form of a global wave, idealized orography in the form of a local table mountain, and the earth's orography. The study confirms the importance of resonances, both in the infinitesimal orography and in the finite orography cases. With finite orography the stationary waves organize themselves into a one-dimensional set of solutions, which due to the resonances, is piecewise connected. It is pointed out that these stationary waves could be relevant for atmospheric regimes.

Use of a crop model ensemble to quantify CO2 stimulation of water-stressed and well-watered crops

Increased atmospheric concentrations of carbon dioxide (CO2) will benefit the yield of most crops. Two free air CO2 enrichment (FACE) meta-analyses have shown increases in yield of between 0 and 73% for C3 crops. Despite this large range, few crop modelling studies quantify the uncertainty inherent in the parameterisation of crop growth and development. We present a novel perturbed-parameter method of crop model simulation, which uses some constraints from observations, that does this. The model used is the groundnut (i.e. peanut; Arachis hypogaea L.) version of the general large-area model for annual crops (GLAM). The conclusions are of relevance to C3 crops in general. The increases in yield simulated by GLAM for doubled CO2 were between 16 and 62%. The difference in mean percentage increase between well-watered and water-stressed simulations was 6.8. These results were compared to FACE and controlled environment studies, and to sensitivity tests on two other crop models of differing levels of complexity: CROPGRO, and the groundnut model of Hammer et al. [Hammer, G.L., Sinclair, T.R., Boote, K.J., Wright, G.C., Meinke, H., Bell, M.J., 1995. A peanut simulation model. I. Model development and testing. Agron. J. 87, 1085-1093]. The relationship between CO2 and water stress in the experiments and in the models was examined. From a physiological perspective, water-stressed crops are expected to show greater CO2 stimulation than well-watered crops. This expectation has been cited in literature. However, this result is not seen consistently in either the FACE studies or in the crop models. In contrast, leaf-level models of assimilation do consistently show this result. An analysis of the evidence from these models and from the data suggests that scale (canopy versus leaf), model calibration, and model complexity are factors in determining the sign and magnitude of the interaction between CO2 and water stress. We conclude from our study that the statement that 'water-stressed crops show greater CO2 stimulation than well-watered crops' cannot be held to be universally true. We also conclude, preliminarily, that the relationship between water stress and assimilation varies with scale. Accordingly, we provide some suggestions on how studies of a similar nature, using crop models of a range of complexity, could contribute further to understanding the roles of model calibration, model complexity and scale. (C) 2008 Elsevier B.V. All rights reserved.

Difference in competence for in vitro proliferation and ex vitro growth of genetically identical mature and juvenile clones of apomictic Malus species

The apomictic system in Malus wits Used Is a model to examine rejuvenation by generating genetically identical tissue culture lines that had two entirely different developmental origins: either embryo-derived tissues (juvenile clones) or somatic tissue from the adult/mature tree (mature clones). These two lines were then subsequently used to examine in vitro difference between mature (M) and juvenile (J) tissues in potential for shoot, root proliferation and ex vitro (glasshouse) growth. The M clones of M. hupehensis and M. toringoides in vitro had significantly fewer total shoots and shoot more than 2 cm in length per proliferating explant than the J clones and also rooted less efficiently. Ex vitro (glasshouse) juvenile clones had shorter internodes, a greater number of leaves and more dry weight compared to their mature counterparts.

Seeds of growth? Agricultural productivity and the transitional dynamics of the Ramsey model

A two-sector Ramsey-type model of growth is developed to investigate the relationship between agricultural productivity and economy-wide growth. The framework takes into account the peculiarities of agriculture both in production ( reliance on a fixed natural resource base) and in consumption (life-sustaining role and low income elasticity of food demand). The transitional dynamics of the model establish that when preferences respect Engel's law, the level and growth rate of agricultural productivity influence the speed of capital accumulation. A calibration exercise shows that a small difference in agricultural productivity has drastic implications for the rate and pattern of growth of the economy. Hence, low agricultural productivity can form a bottleneck limiting growth, because high food prices result in a low saving rate.

Evaluation of a predictive model for Mycosphaerella graminicola for economic and environmental benefits

A method was developed to evaluate crop disease predictive models for their economic and environmental benefits. Benefits were quantified as the value of a prediction measured by costs saved and fungicide dose saved. The value of prediction was defined as the net gain made by using predictions, measured as the difference between a scenario where predictions are available and used and a scenario without prediction. Comparable 'with' and 'without' scenarios were created with the use of risk levels. These risk levels were derived from a probability distribution fitted to observed disease severities. These distributions were used to calculate the probability that a certain disease induced economic loss was incurred. The method was exemplified by using it to evaluate a model developed for Mycosphaerella graminicola risk prediction. Based on the value of prediction, the tested model may have economic and environmental benefits to growers if used to guide treatment decisions on resistant cultivars. It is shown that the value of prediction measured by fungicide dose saved and costs saved is constant with the risk level. The model could also be used to evaluate similar crop disease predictive models.

Self assembly of a model amphiphilic phenylalanine peptide/polyethylene glycol block copolymer in aqueous solution

There has been great interest recently in peptide amphiphiles and block copolymers containing biomimetic peptide sequences due to applications in bionanotechnology. We investigate the self-assembly of the peptide-PEG amphiphile FFFF-PEG5000 containing the hydrophobic sequence of four phenylalanine residues conjugated to PEG of molar mass 5000. This serves as a simple model peptide amphiphile. At very low concentration, association of hydrophobic aromatic phenylalanine residues occurs, as revealed by circular dichroism and UV/vis fluorescence experiments. A critical aggregation concentration associated with the formation of hydrophobic domains is determined through pyrene fluorescence assays. At higher concentration, defined beta-sheets develop as revealed by FTIR spectroscopy and X-ray diffraction. Transmission electron microscopy reveals self-assembled straight fibril structures. These are much shorter than those observed for amyloid peptides, the finite length may be set by the end cap energy due to the hydrophobicity of phenylalanine. The combination of these techniques points to different aggregation processes depending on concentration. Hydrophobic association into irregular aggregates occurs at low concentration, well-developed beta-sheets only developing at higher concentration. Drying of FFFF-PEG5000 solutions leads to crystallization of PEG, as confirmed by polarized optical microscopy (POM), FTIR and X-ray diffraction (XRD). PEG crystallization does not disrupt local beta-sheet structure (as indicated by FTIR and XRD). However on longer lengthscales the beta-sheet fibrillar structure is perturbed because spheruilites from PEG crystallization are observed by POM. (C) 2009 Elsevier B.V. All rights reserved.

Effect of fecal water on an in vitro model of colonic mucosal barrier function

Fecal water (FW) has been shown to exert, in cultured cells, cytotoxic and genotoxic effects that have implications for colorectal cancer (CRC) risk. We have investigated a further biological activity of FW, namely, the ability to affect gap junctions in CACO2 cell monolayers as an index of mucosal barrier function, which is known to be disrupted in cancer. FW samples fi-om healthy, free-living, European subjects that were divided into two broad age groups, adult (40 +/- 9.7 yr; n = 53) and elderly (76 +/- 7.5 yr; n = 55) were tested for effects on gap junction using the transepithelial resistance (TER) assay. Overall, treatment of CACO2 cells with FW samples fi-om adults increased TER (+ 4 %), whereas FW from elderly subjects decreased TER (-5%); the difference between the two groups was significant (P < 0.05). We also measured several components of FW potentially associated with modulation of TER, namely, short-chain fatty acid (SCFA) and ammonia. SCFAs (propionic, acetic, and n-butyric) were significantly lower in the elderly population (-30%, -35%, and -21%, respectively, all P pound 0.01). We consider that FW modulation of in vitro epithelial barrier function is a potentially useful noninvasive biomarker, but it requires further validation to establish its relationship to CRC risk.

A neurofuzzy network knowledge extraction and extended gram-Schmidt algorithm for model subspace decomposition

This paper introduces a new neurofuzzy model construction and parameter estimation algorithm from observed finite data sets, based on a Takagi and Sugeno (T-S) inference mechanism and a new extended Gram-Schmidt orthogonal decomposition algorithm, for the modeling of a priori unknown dynamical systems in the form of a set of fuzzy rules. The first contribution of the paper is the introduction of a one to one mapping between a fuzzy rule-base and a model matrix feature subspace using the T-S inference mechanism. This link enables the numerical properties associated with a rule-based matrix subspace, the relationships amongst these matrix subspaces, and the correlation between the output vector and a rule-base matrix subspace, to be investigated and extracted as rule-based knowledge to enhance model transparency. The matrix subspace spanned by a fuzzy rule is initially derived as the input regression matrix multiplied by a weighting matrix that consists of the corresponding fuzzy membership functions over the training data set. Model transparency is explored by the derivation of an equivalence between an A-optimality experimental design criterion of the weighting matrix and the average model output sensitivity to the fuzzy rule, so that rule-bases can be effectively measured by their identifiability via the A-optimality experimental design criterion. The A-optimality experimental design criterion of the weighting matrices of fuzzy rules is used to construct an initial model rule-base. An extended Gram-Schmidt algorithm is then developed to estimate the parameter vector for each rule. This new algorithm decomposes the model rule-bases via an orthogonal subspace decomposition approach, so as to enhance model transparency with the capability of interpreting the derived rule-base energy level. This new approach is computationally simpler than the conventional Gram-Schmidt algorithm for resolving high dimensional regression problems, whereby it is computationally desirable to decompose complex models into a few submodels rather than a single model with large number of input variables and the associated curse of dimensionality problem. Numerical examples are included to demonstrate the effectiveness of the proposed new algorithm.

Model selection approaches for non-linear system identification: a review

The identification of non-linear systems using only observed finite datasets has become a mature research area over the last two decades. A class of linear-in-the-parameter models with universal approximation capabilities have been intensively studied and widely used due to the availability of many linear-learning algorithms and their inherent convergence conditions. This article presents a systematic overview of basic research on model selection approaches for linear-in-the-parameter models. One of the fundamental problems in non-linear system identification is to find the minimal model with the best model generalisation performance from observational data only. The important concepts in achieving good model generalisation used in various non-linear system-identification algorithms are first reviewed, including Bayesian parameter regularisation and models selective criteria based on the cross validation and experimental design. A significant advance in machine learning has been the development of the support vector machine as a means for identifying kernel models based on the structural risk minimisation principle. The developments on the convex optimisation-based model construction algorithms including the support vector regression algorithms are outlined. Input selection algorithms and on-line system identification algorithms are also included in this review. Finally, some industrial applications of non-linear models are discussed.

Finite- N effects for ideal polymer chains near a flat impenetrable wall

This paper addresses the statistical mechanics of ideal polymer chains next to a hard wall. The principal quantity of interest, from which all monomer densities can be calculated, is the partition function, G N(z) , for a chain of N discrete monomers with one end fixed a distance z from the wall. It is well accepted that in the limit of infinite N , G N(z) satisfies the diffusion equation with the Dirichlet boundary condition, G N(0) = 0 , unless the wall possesses a sufficient attraction, in which case the Robin boundary condition, G N(0) = - x G N ′(0) , applies with a positive coefficient, x . Here we investigate the leading N -1/2 correction, D G N(z) . Prior to the adsorption threshold, D G N(z) is found to involve two distinct parts: a Gaussian correction (for z <~Unknown control sequence '\lesssim' aN 1/2 with a model-dependent amplitude, A , and a proximal-layer correction (for z <~Unknown control sequence '\lesssim' a described by a model-dependent function, B(z).

Evaluation of a Lagrangian box model using field measurements from EASE (Eastern Atlantic Summer Experiment) 1996

The Cambridge Tropospheric Trajectory model of Chemistry and Transport (CiTTyCAT), a Lagrangian chemistry model, has been evaluated using atmospheric chemical measurements collected during the East Atlantic Summer Experiment 1996 (EASE '96). This field campaign was part of the UK Natural Environment Research Council's (NERC) Atmospheric Chemistry Studies in the Oceanic Environment (ACSOE) programme, conducted at Mace Head, Republic of Ireland, during July and August 1996. The model includes a description of gas-phase tropospheric chemistry, and simple parameterisations for surface deposition, mixing from the free troposphere and emissions. The model generally compares well with the measurements and is used to study the production and loss of O3 under a variety of conditions. The mean difference between the hourly O3 concentrations calculated by the model and those measured is 0.6 ppbv with a standard deviation of 8.7 ppbv. Three specific air-flow regimes were identified during the campaign – westerly, anticyclonic (easterly) and south westerly. The westerly flow is typical of background conditions for Mace Head. However, on some occasions there was evidence of long-range transport of pollutants from North America. In periods of anticyclonic flow, air parcels had collected emissions of NOx and VOCs immediately before arriving at Mace Head, leading to O3 production. The level of calculated O3 depends critically on the precise details of the trajectory, and hence on the emissions into the air parcel. In several periods of south westerly flow, low concentrations of O3 were measured which were consistent with deposition and photochemical destruction inside the tropical marine boundary layer.

Model inter-comparison between statistical and dynamic model assessments of the long-term stability of blanket peat in Great Britain (1940-2099)

We compared output from 3 dynamic process-based models (DMs: ECOSSE, MILLENNIA and the Durham Carbon Model) and 9 bioclimatic envelope models (BCEMs; including BBOG ensemble and PEATSTASH) ranging from simple threshold to semi-process-based models. Model simulations were run at 4 British peatland sites using historical climate data and climate projections under a medium (A1B) emissions scenario from the 11-RCM (regional climate model) ensemble underpinning UKCP09. The models showed that blanket peatlands are vulnerable to projected climate change; however, predictions varied between models as well as between sites. All BCEMs predicted a shift from presence to absence of a climate associated with blanket peat, where the sites with the lowest total annual precipitation were closest to the presence/absence threshold. DMs showed a more variable response. ECOSSE predicted a decline in net C sink and shift to net C source by the end of this century. The Durham Carbon Model predicted a smaller decline in the net C sink strength, but no shift to net C source. MILLENNIA predicted a slight overall increase in the net C sink. In contrast to the BCEM projections, the DMs predicted that the sites with coolest temperatures and greatest total annual precipitation showed the largest change in carbon sinks. In this model inter-comparison, the greatest variation in model output in response to climate change projections was not between the BCEMs and DMs but between the DMs themselves, because of different approaches to modelling soil organic matter pools and decomposition amongst other processes. The difference in the sign of the response has major implications for future climate feedbacks, climate policy and peatland management. Enhanced data collection, in particular monitoring peatland response to current change, would significantly improve model development and projections of future change.

Finite departure from convective quasi-equilibrium: periodic cycle and discharge-recharge mechanism

A simple self–contained theory is proposed for describing life cycles of convective systems as a discharge–recharge process. A closed description is derived for the dynamics of an ensemble of convective plumes based on an energy cycle. The system consists of prognostic equations for the cloud work function and the convective kinetic energy. The system can be closed by intro ducing a functional relationship between the convective kinetic energy and the cloud–base mass flux. The behaviour of this system is considered under a bulk simplification. Previous cloud–resolving mo delling as well as bulk statistical theories for ensemble convective systems suggest that a plausible relationship would be to assume that the convective kinetic energy is linearly proportional to the cloud–base mass flux. As a result, the system reduces to a nonlinear dynamical system with two dependent variables, the cloud–base mass flux and the cloud work function. The fully nonlinear solution of this system always represents a periodic cycle regardless of the initial condition under constant large–scale forcing. Importantly, the inclusion of energy dissipation in this model does not in itself lead the system to an equilibrium.

Comparing the UK Met Office climate prediction system DePreSys with idealized predictability in the HadCM3 model

The initial condition effect on climate prediction skill over a 2-year hindcast time-scale has been assessed from ensemble HadCM3 climate model runs using anomaly initialization over the period 1990–2001, and making comparisons with runs without initialization (equivalent to climatological conditions), and to anomaly persistence. It is shown that the assimilation improves the prediction skill in the first year globally, and in a number of limited areas out into the second year. Skill in hindcasting surface air temperature anomalies is most marked over ocean areas, and is coincident with areas of high sea surface temperature and ocean heat content skill. Skill improvement over land areas is much more limited but is still detectable in some cases. We found little difference in the skill of hindcasts using three different sets of ocean initial conditions, and we obtained the best results by combining these to form a grand ensemble hindcast set. Results are also compared with the idealized predictability studies of Collins (Clim. Dynam. 2002; 19: 671–692), which used the same model. The maximum lead time for which initialization gives enhanced skill over runs without initialization varies in different regions but is very similar to lead times found in the idealized studies, therefore strongly supporting the process representation in the model as well as its use for operational predictions. The limited 12-year period of the study, however, means that the regional details of model skill should probably be further assessed under a wider range of observational conditions.

A detailed single chain model for molten poly(tetrafluoroethylene) from a novel structure refinement technique on neutron scattering data

We present a new methodology that couples neutron diffraction experiments over a wide Q range with single chain modelling in order to explore, in a quantitative manner, the intrachain organization of non-crystalline polymers. The technique is based on the assignment of parameters describing the chemical, geometric and conformational characteristics of the polymeric chain, and on the variation of these parameters to minimize the difference between the predicted and experimental diffraction patterns. The method is successfully applied to the study of molten poly(tetrafluoroethylene) at two different temperatures, and provides unambiguous information on the configuration of the chain and its degree of flexibility. From analysis of the experimental data a model is derived with CC and CF bond lengths of 1.58 and 1.36 Å, respectively, a backbone valence angle of 110° and a torsional angle distribution which is characterized by four isometric states, namely a split trans state at ± 18°, giving rise to a helical chain conformation, and two gauche states at ± 112°. The probability of trans conformers is 0.86 at T = 350°C, which decreases slightly to 0.84 at T = 400°C. Correspondingly, the chain segments are characterized by long all-trans sequences with random changes in sign, rather anisotropic in nature, which give rise to a rather stiff chain. We compare the results of this quantitative analysis of the experimental scattering data with the theoretical predictions of both force fields and molecular orbital conformation energy calculations.