Climate change and heat-related mortality in six cities Part 1: model construction and validation

Heat waves are expected to increase in frequency and magnitude with climate change. The first part of a study to produce projections of the effect of future climate change on heat-related mortality is presented. Separate city-specific empirical statistical models that quantify significant relationships between summer daily maximum temperature (T max) and daily heat-related deaths are constructed from historical data for six cities: Boston, Budapest, Dallas, Lisbon, London, and Sydney. ‘Threshold temperatures’ above which heat-related deaths begin to occur are identified. The results demonstrate significantly lower thresholds in ‘cooler’ cities exhibiting lower mean summer temperatures than in ‘warmer’ cities exhibiting higher mean summer temperatures. Analysis of individual ‘heat waves’ illustrates that a greater proportion of mortality is due to mortality displacement in cities with less sensitive temperature–mortality relationships than in those with more sensitive relationships, and that mortality displacement is no longer a feature more than 12 days after the end of the heat wave. Validation techniques through residual and correlation analyses of modelled and observed values and comparisons with other studies indicate that the observed temperature–mortality relationships are represented well by each of the models. The models can therefore be used with confidence to examine future heat-related deaths under various climate change scenarios for the respective cities (presented in Part 2).

Flux comparison of Eulerian and Lagrangian estimates of Agulhas leakage: A case study using a numerical model

Estimating the magnitude of Agulhas leakage, the volume flux of water from the Indian to the Atlantic Ocean, is difficult because of the presence of other circulation systems in the Agulhas region. Indian Ocean water in the Atlantic Ocean is vigorously mixed and diluted in the Cape Basin. Eulerian integration methods, where the velocity field perpendicular to a section is integrated to yield a flux, have to be calibrated so that only the flux by Agulhas leakage is sampled. Two Eulerian methods for estimating the magnitude of Agulhas leakage are tested within a high-resolution two-way nested model with the goal to devise a mooring-based measurement strategy. At the GoodHope line, a section halfway through the Cape Basin, the integrated velocity perpendicular to that line is compared to the magnitude of Agulhas leakage as determined from the transport carried by numerical Lagrangian floats. In the first method, integration is limited to the flux of water warmer and more saline than specific threshold values. These threshold values are determined by maximizing the correlation with the float-determined time series. By using the threshold values, approximately half of the leakage can directly be measured. The total amount of Agulhas leakage can be estimated using a linear regression, within a 90% confidence band of 12 Sv. In the second method, a subregion of the GoodHope line is sought so that integration over that subregion yields an Eulerian flux as close to the float-determined leakage as possible. It appears that when integration is limited within the model to the upper 300 m of the water column within 900 km of the African coast the time series have the smallest root-mean-square difference. This method yields a root-mean-square error of only 5.2 Sv but the 90% confidence band of the estimate is 20 Sv. It is concluded that the optimum thermohaline threshold method leads to more accurate estimates even though the directly measured transport is a factor of two lower than the actual magnitude of Agulhas leakage in this model.

Design and optimisation of a large-area process-based model for annual crops

The formulation of a new process-based crop model, the general large-area model (GLAM) for annual crops is presented. The model has been designed to operate on spatial scales commensurate with those of global and regional climate models. It aims to simulate the impact of climate on crop yield. Procedures for model parameter determination and optimisation are described, and demonstrated for the prediction of groundnut (i.e. peanut; Arachis hypogaea L.) yields across India for the period 1966-1989. Optimal parameters (e.g. extinction coefficient, transpiration efficiency, rate of change of harvest index) were stable over space and time, provided the estimate of the yield technology trend was based on the full 24-year period. The model has two location-specific parameters, the planting date, and the yield gap parameter. The latter varies spatially and is determined by calibration. The optimal value varies slightly when different input data are used. The model was tested using a historical data set on a 2.5degrees x 2.5degrees grid to simulate yields. Three sites are examined in detail-grid cells from Gujarat in the west, Andhra Pradesh towards the south, and Uttar Pradesh in the north. Agreement between observed and modelled yield was variable, with correlation coefficients of 0.74, 0.42 and 0, respectively. Skill was highest where the climate signal was greatest, and correlations were comparable to or greater than correlations with seasonal mean rainfall. Yields from all 35 cells were aggregated to simulate all-India yield. The correlation coefficient between observed and simulated yields was 0.76, and the root mean square error was 8.4% of the mean yield. The model can be easily extended to any annual crop for the investigation of the impacts of climate variability (or change) on crop yield over large areas. (C) 2004 Elsevier B.V. All rights reserved.

Mixed Logit Model Performance and Distributional Assumptions: Preferences and GM foods

Using mixed logit models to analyse choice data is common but requires ex ante specification of the functional forms of preference distributions. We make the case for greater use of bounded functional forms and propose the use of the Marginal Likelihood, calculated using Bayesian techniques, as a single measure of model performance across non nested mixed logit specifications. Using this measure leads to very different rankings of model specifications compared to alternative rule of thumb measures. The approach is illustrated using data from a choice experiment regarding GM food types which provides insights regarding the recent WTO dispute between the EU and the US, Canada and Argentina and whether labelling and trade regimes should be based on the production process or product composition.

A model of seed dormancy in wild oats (Avena fatua) for investigating genotype x environment interactions

Dormancy is an adaptive trait in seed populations that helps ensure that seed germination is distributed over time and occurs in environmental conditions suitable for seedling growth. Several genes.. associated with seed dormancy in various plant species, have been integrated into a hypothetical dormancy model for Avena fatua L. (wild oats). Generally, the synthesis of, and sensitivity to, abscisic acid (ABA) during imbibition determines whether genes similar to those during maturation are expressed leading to a maintenance of dormancy during extended imbibition. Alternatively, there may be a shift towards expression of genes associated with gibberellins leading to germination. Environmental factors during maturation, after-ripening and imbibition are likely to interact with the genotype to affect gene expression and hence whether or not a seed germinates. In spite of the difficulties of working on a hexaploid species, A. fatua was selected for study because of its worldwide importance as a weed. Dormant and non-dormant genotypes of this species were also available. Gene expression studies are being carried out on three A.fatua genotypes produced tinder different environmental conditions to investigate the role of specific genes in dormancy and genotype X environment interactions in relation to dormancy.

Culture and food: a model of yoghurt consumption in the EU

The completion of the Single European Market was expected to create a large market that would enable firms to capture economies of scale that would in turn result in lower prices to European consumers. These benefits are only likely to be realised if consumers in the various countries of the EU wish to consume the same products and respond to similar marketing strategies (with respect to promotion, distribution etc). This study examines, through a model of yoghurt consumption, whether cultural differences continue to determine food-related behaviour in the EU. The model is derived from the marketing literature and views the consumption decision as the outcome of a multi-stage process in which yoghurt knowledge, attitudes to different yoghurt attributes (such as bio-bifidus, low-fat, organic) and overall attitude towards yoghurt as a product all feed into the frequency with which yoghurt is consumed at breakfast, as a snack and as a dessert. The model uses data collected from a consumer survey in I I European countries and is estimated using probit and ordinal probit methods. The results suggest that important cultural differences continue to determine food-related behaviour in the I I countries of the study. (c) 2004 Elsevier Ltd. All rights reserved.

The tissue-specific processing of pro-opiomelanocortin

It is just over 30 years since the definitive identification of the adrenocorticotrophin (ACTH) precursor, pro-opiomelanocotin (POMC). Although first characterised in the anterior and intermediate lobes of the pituitary, POMC is also expressed in a number of both central and peripheral tissues including the skin, central nervous tissue and placenta. Following synthesis, POMC undergoes extensive post-translational processing producing not only ACTH, but also a number of other biologically active peptides. The extent and pattern of this processing is tissue-specific, the end result being the tissue dependent production of different combinations of peptides from the same precursor. These peptides have a diverse range of biological roles ranging from pigmentation to adrenal function to the regulation of feeding. This level of complexity has resulted in POMC becoming the archetypal model for prohormone processing, illustrating how a single protein combined with post-translational modification can have a diverse number of roles.

Overview of mathematical approaches used to model bacterial chemotaxis I: the single cell

Mathematical modeling of bacterial chemotaxis systems has been influential and insightful in helping to understand experimental observations. We provide here a comprehensive overview of the range of mathematical approaches used for modeling, within a single bacterium, chemotactic processes caused by changes to external gradients in its environment. Specific areas of the bacterial system which have been studied and modeled are discussed in detail, including the modeling of adaptation in response to attractant gradients, the intracellular phosphorylation cascade, membrane receptor clustering, and spatial modeling of intracellular protein signal transduction. The importance of producing robust models that address adaptation, gain, and sensitivity are also discussed. This review highlights that while mathematical modeling has aided in understanding bacterial chemotaxis on the individual cell scale and guiding experimental design, no single model succeeds in robustly describing all of the basic elements of the cell. We conclude by discussing the importance of this and the future of modeling in this area.

Apical leaf necrosis and leaf nitrogen dynamics in diseased leaves: a model study

Apical leaf necrosis is a physiological process related to nitrogen (N) dynamics in the leaf. Pathogens use leaf nutrients and can thus accelerate this physiological apical necrosis. This process differs from necrosis occurring around pathogen lesions (lesion-induced necrosis), which is a direct result of the interaction between pathogen hyphae and leaf cells. This paper primarily concentrates on apical necrosis, only incorporating lesion-induced necrosis by necessity. The relationship between pathogen dynamics and physiological apical leaf necrosis is modelled through leaf nitrogen dynamics. The specific case of Puccinia triticina infections on Triticum aestivum flag leaves is studied. In the model, conversion of indirectly available N in the form of, for example, leaf cell proteins (N-2(t)) into directly available N (N-1(t), i.e. the form of N that can directly be used by either pathogen or plant sinks) results in apical necrosis. The model reproduces observed trends of disease severity, apical necrosis and green leaf area (GLA) and leaf N dynamics of uninfected and infected leaves. Decreasing the initial amount of directly available N results in earlier necrosis onset and longer necrosis duration. Decreasing the initial amount of indirectly available N, has no effect on necrosis onset and shortens necrosis duration. The model could be used to develop hypotheses on how the disease-GLA relation affects yield loss, which can be tested experimentally. Upon incorporation into crop simulation models, the model might provide a tool to more accurately estimate crop yield and effects of disease management strategies in crops sensitive to fungal pathogens.

Diagnosis of variability and trends in a global precipitation dataset using a physically motivated statistical model

A physically motivated statistical model is used to diagnose variability and trends in wintertime ( October - March) Global Precipitation Climatology Project (GPCP) pentad (5-day mean) precipitation. Quasi-geostrophic theory suggests that extratropical precipitation amounts should depend multiplicatively on the pressure gradient, saturation specific humidity, and the meridional temperature gradient. This physical insight has been used to guide the development of a suitable statistical model for precipitation using a mixture of generalized linear models: a logistic model for the binary occurrence of precipitation and a Gamma distribution model for the wet day precipitation amount. The statistical model allows for the investigation of the role of each factor in determining variations and long-term trends. Saturation specific humidity q(s) has a generally negative effect on global precipitation occurrence and with the tropical wet pentad precipitation amount, but has a positive relationship with the pentad precipitation amount at mid- and high latitudes. The North Atlantic Oscillation, a proxy for the meridional temperature gradient, is also found to have a statistically significant positive effect on precipitation over much of the Atlantic region. Residual time trends in wet pentad precipitation are extremely sensitive to the choice of the wet pentad threshold because of increasing trends in low-amplitude precipitation pentads; too low a choice of threshold can lead to a spurious decreasing trend in wet pentad precipitation amounts. However, for not too small thresholds, it is found that the meridional temperature gradient is an important factor for explaining part of the long-term trend in Atlantic precipitation.

Graphical tracking systems revisited: a practical approach to computer scheduling in horticulture

Graphical tracking is a technique for crop scheduling where the actual plant state is plotted against an ideal target curve which encapsulates all crop and environmental characteristics. Management decisions are made on the basis of the position of the actual crop against the ideal position. Due to the simplicity of the approach it is possible for graphical tracks to be developed on site without the requirement for controlled experimentation. Growth models and graphical tracks are discussed, and an implementation of the Richards curve for graphical tracking described. In many cases, the more intuitively desirable growth models perform sub-optimally due to problems with the specification of starting conditions, environmental factors outside the scope of the original model and the introduction of new cultivars. Accurate specification for a biological model requires detailed and usually costly study, and as such is not adaptable to a changing cultivar range and changing cultivation techniques. Fitting of a new graphical track for a new cultivar can be conducted on site and improved over subsequent seasons. Graphical tracking emphasises the current position relative to the objective, and as such does not require the time consuming or system specific input of an environmental history, although it does require detailed crop measurement. The approach is flexible and could be applied to a variety of specification metrics, with digital imaging providing a route for added value. For decision making regarding crop manipulation from the observed current state, there is a role for simple predictive modelling over the short term to indicate the short term consequences of crop manipulation.

A genetic linkage map of microsatellite, gene-specific and morphological markers in diploid Fragaria

Diploid Fragaria provide a potential model for genomic studies in the Rosaceae. To develop a genetic linkage map of diploid Fragaria, we scored 78 markers (68 microsatellites, one sequence-characterised amplified region, six gene-specific markers and three morphological traits) in an interspecific F2 population of 94 plants generated from a cross of F.vesca f. semperflorens × F. nubicola. Co-segregation analysis arranged 76 markers into seven discrete linkage groups covering 448 cM, with linkage group sizes ranging from 100.3 cM to 22.9 cM. Marker coverage was generally good; however some clustering of markers was observed on six of the seven linkage groups. Segregation distortion was observed at a high proportion of loci (54%), which could reflect the interspecific nature of the progeny and, in some cases, the self-incompatibility of F. nubicola. Such distortion may also account for some of the marker clustering observed in the map. One of the morphological markers, pale-green leaf (pg) has not previously been mapped in Fragaria and was located to the mid-point of linkage group VI. The transferable nature of the markers used in this study means that the map will be ideal for use as a framework for additional marker incorporation aimed at enhancing and resolving map coverage of the diploid Fragaria genome. The map also provides a sound basis for linkage map transfer to the cultivated octoploid strawberry.

Generation of anti-complement "prodrugs": cleavable reagents for specific delivery of complement regulators to disease sites

Expression of biologically active molecules as fusion proteins with antibody Fc can substantially extend the plasma half-life of the active agent but may also influence function. We have previously generated a number of fusion proteins comprising a complement regulator coupled to Fc and shown that the hybrid molecule has a long plasma half-life and retains biological activity. However, several of the fusion proteins generated had substantially reduced biological activity when compared with the native regulator or regulator released from the Fc following papain cleavage. We have taken advantage of this finding to engineer a prodrug with low complement regulatory activity that is cleaved at sites of inflammation to release active regulator. Two model prodrugs, comprising, respectively, the four short consensus repeats of human decay accelerating factor (CD55) linked to IgG4 Fc and the three NH2-terminal short consensus repeats of human decay accelerating factor linked to IgG2 Fc have been developed. In each, specific cleavage sites for matrix metalloproteinases and/or aggrecanases have been incorporated between the complement regulator and the Fc. These prodrugs have markedly decreased complement inhibitory activity when compared with the parent regulator in vitro. Exposure of the prodrugs to the relevant enzymes, either purified, or in supernatants of cytokine-stimulated chondrocytes or in synovial fluid, efficiently cleaved the prodrug, releasing active regulator. Such agents, having negligible systemic effects but active at sites of inflammation, represent a paradigm for the next generation of anti-C therapeutics.

Qualitative modelling of sustainable energy scenarios: An application and extension of the Bon Qualitative Input-Output model

Climate change is one of the major challenges facing economic systems at the start of the 21st century. Reducing greenhouse gas emissions will require both restructuring the energy supply system (production) and addressing the efficiency and sufficiency of the social uses of energy (consumption). The energy production system is a complicated supply network of interlinked sectors with 'knock-on' effects throughout the economy. End use energy consumption is governed by complex sets of interdependent cultural, social, psychological and economic variables driven by shifts in consumer preference and technological development trajectories. To date, few models have been developed for exploring alternative joint energy production-consumption systems. The aim of this work is to propose one such model. This is achieved in a methodologically coherent manner through integration of qualitative input-output models of production, with Bayesian belief network models of consumption, at point of final demand. The resulting integrated framework can be applied either (relatively) quickly and qualitatively to explore alternative energy scenarios, or as a fully developed quantitative model to derive or assess specific energy policy options. The qualitative applications are explored here.

Mathematical model for low density lipoprotein (LDL) endocytosis by hepatocytes

Individuals with elevated levels of plasma low density lipoprotein (LDL) cholesterol (LDL-C) are considered to be at risk of developing coronary heart disease. LDL particles are removed from the blood by a process known as receptor-mediated endocytosis, which occurs mainly in the liver. A series of classical experiments delineated the major steps in the endocytotic process; apolipoprotein B-100 present on LDL particles binds to a specific receptor (LDL receptor, LDL-R) in specialized areas of the cell surface called clathrin-coated pits. The pit comprising the LDL-LDL-R complex is internalized forming a cytoplasmic endosome. Fusion of the endosome with a lysosome leads to degradation of the LDL into its constituent parts (that is, cholesterol, fatty acids, and amino acids), which are released for reuse by the cell, or are excreted. In this paper, we formulate a mathematical model of LDL endocytosis, consisting of a system of ordinary differential equations. We validate our model against existing in vitro experimental data, and we use it to explore differences in system behavior when a single bolus of extracellular LDL is supplied to cells, compared to when a continuous supply of LDL particles is available. Whereas the former situation is common to in vitro experimental systems, the latter better reflects the in vivo situation. We use asymptotic analysis and numerical simulations to study the longtime behavior of model solutions. The implications of model-derived insights for experimental design are discussed.