Effects of spatial resolution in the simulation of daily and subdaily precipitation in the southwestern US

We evaluate the effects of spatial resolution on the ability of a regional climate model to reproduce observed extreme precipitation for a region in the Southwestern United States. A total of 73 National Climate Data Center observational sites spread throughout Arizona and New Mexico are compared with regional climate simulations at the spatial resolutions of 50 km and 10 km for a 31 year period from 1980 to 2010. We analyze mean, 3-hourly and 24-hourly extreme precipitation events using WRF regional model simulations driven by NCEP-2 reanalysis. The mean climatological spatial structure of precipitation in the Southwest is well represented by the 10 km resolution but missing in the coarse (50 km resolution) simulation. However, the fine grid has a larger positive bias in mean summer precipitation than the coarse-resolution grid. The large overestimation in the simulation is in part due to scale-dependent deficiencies in the Kain-Fritsch convective parameterization scheme that generate excessive precipitation and induce a slow eastward propagation of the moist convective summer systems in the high-resolution simulation. Despite this overestimation in the mean, the 10 km simulation captures individual extreme summer precipitation events better than the 50 km simulation. In winter, however, the two simulations appear to perform equally in simulating extremes.

The role of air–sea coupling in the simulation of the Madden–Julian oscillation in the Hadley Centre model

The role of air–sea coupling in the simulation of the Madden–Julian oscillation (MJO) is explored using two configurations of the Hadley Centre atmospheric model (AGCM), GA3.0, which differ only in F, a parameter controlling convective entrainment and detrainment. Increasing F considerably improves deficient MJO-like variability in the Indian and Pacific Oceans, but variability in and propagation through the Maritime Continent remains weak. By coupling GA3.0 in the tropical Indo-Pacific to a boundary-layer ocean model, KPP, and employing climatological temperature corrections, well resolved air–sea interactions are simulated with limited alterations to the mean state. At default F, when GA3.0 has a poor MJO, coupling produces a stronger MJO with some eastward propagation, although both aspects remain deficient. These results agree with previous sensitivity studies using AGCMs with poor variability. At higher F, coupling does not affect MJO amplitude but enhances propagation through the Maritime Continent, resulting in an MJO that resembles observations. A sensitivity experiment with coupling in only the Indian Ocean reverses these improvements, suggesting coupling in the Maritime Continent and West Pacific is critical for propagation. We hypothesise that for AGCMs with a poor MJO, coupling provides a “crutch” to artificially augment MJO-like activity through high-frequency SST anomalies. In related experiments, we employ the KPP framework to analyse the impact of air–sea interactions in the fully coupled GA3.0, which at default F shows a similar MJO to uncoupled GA3.0. This is due to compensating effects: an improvement from coupling and a degradation from mean-state errors. Future studies on the role of coupling should carefully separate these effects.

Boundary condition effects in the simulation study of equilibrium properties of magnetic dipolar fluids

In this paper we investigate the equilibrium properties of magnetic dipolar (ferro-) fluids and discuss finite-size effects originating from the use of different boundary conditions in computer simulations. Both periodic boundary conditions and a finite spherical box are studied. We demonstrate that periodic boundary conditions and subsequent use of Ewald sum to account for the long-range dipolar interactions lead to a much faster convergence (in terms of the number of investigated dipolar particles) of the magnetization curve and the initial susceptibility to their thermodynamic limits. Another unwanted effect of the simulations in a finite spherical box geometry is a considerable sensitivity to the container size. We further investigate the influence of the surface term in the Ewald sum-that is, due to the surrounding continuum with magnetic permeability mu(BC)-on the convergence properties of our observables and on the final results. The two different ways of evaluating the initial susceptibility, i.e., (1) by the magnetization response of the system to an applied field and (2) by the zero-field fluctuation of the mean-square dipole moment of the system, are compared in terms of speed and accuracy.

Understanding advances in the simulation of the Madden-Julian Oscillation in a numerical weather prediction model.

The Madden-Julian Oscillation (MJO) is the dominant mode of intraseasonal variability in the Trop- ics. It can be characterised as a planetary-scale coupling between the atmospheric circulation and organised deep convection that propagates east through the equatorial Indo-Pacific region. The MJO interacts with weather and climate systems on a near-global scale and is a crucial source of predictability for weather forecasts on medium to seasonal timescales. Despite its global signifi- cance, accurately representing the MJO in numerical weather prediction (NWP) and climate models remains a challenge. This thesis focuses on the representation of the MJO in the Integrated Forecasting System (IFS) at the European Centre for Medium-Range Weather Forecasting (ECMWF), a state-of-the-art NWP model. Recent modifications to the model physics in Cycle 32r3 (Cy32r3) of the IFS led to ad- vances in the simulation of the MJO; for the first time the observed amplitude of the MJO was maintained throughout the integration period. A set of hindcast experiments, which differ only in their formulation of convection, have been performed between May 2008 and April 2009 to asses the sensitivity of MJO simulation in the IFS to the Cy32r3 convective parameterization. Unique to this thesis is the attribution of the advances in MJO simulation in Cy32r3 to the mod- ified convective parameterization, specifically, the relative-humidity-dependent formulation for or- ganised deep entrainment. Increasing the sensitivity of the deep convection scheme to environmen- tal moisture is shown to modify the relationship between precipitation and moisture in the model. Through dry-air entrainment, convective plumes ascending in low-humidity environments terminate lower in the atmosphere. As a result, there is an increase in the occurrence of cumulus congestus, which acts to moisten the mid-troposphere. Due to the modified precipitation-moisture relationship more moisture is able to build up which effectively preconditions the tropical atmosphere for the transition to deep convection. Results from this thesis suggest that a tropospheric moisture control on convection is key to simulating the interaction between the physics and large-scale circulation associated with the MJO.

Effects of dendritic morphology on CA3 pyramidal cell electrophysiology: a simulation study

We investigated the effect of morphological differences on neuronal firing behavior within the hippocampal CA3 pyramidal cell family by using three-dimensional reconstructions of dendritic morphology in computational simulations of electrophysiology. In this paper, we report for the first time that differences in dendritic structure within the same morphological class can have a dramatic influence on the firing rate and firing mode (spiking versus bursting and type of bursting). Our method consisted of converting morphological measurements from three-dimensional neuroanatomical data of CA3 pyramidal cells into a computational simulator format. In the simulation, active channels were distributed evenly across the cells so that the electrophysiological differences observed in the neurons would only be due to morphological differences. We found that differences in the size of the dendritic tree of CA3 pyramidal cells had a significant qualitative and quantitative effect on the electrophysiological response. Cells with larger dendritic trees: (1) had a lower burst rate, but a higher spike rate within a burst, (2) had higher thresholds for transitions from quiescent to bursting and from bursting to regular spiking and (3) tended to burst with a plateau. Dendritic tree size alone did not account for all the differences in electrophysiological responses. Differences in apical branching, such as the distribution of branch points and terminations per branch order, appear to effect the duration of a burst. These results highlight the importance of considering the contribution of morphology in electrophysiological and simulation studies.

Spatial and temporal effects of free-air CO2 enrichment (POPFACE) on leaf growth, cell expansion and cell production in a closed canopy of poplar

Leaf expansion in the fast-growing tree,Populus × euramericana was stimulated by elevated [CO2] in a closed-canopy forest plantation, exposed using a free air CO2 enrichment technique enabling long-term experimentation in field conditions. The effects of elevated [CO2] over time were characterized and related to the leaf plastochron index (LPI), and showed that leaf expansion was stimulated at very early (LPI, 0–3) and late (LPI, 6–8) stages in development. Early and late effects of elevated [CO2] were largely the result of increased cell expansion and increased cell production, respectively. Spatial effects of elevated [CO2] were also marked and increased final leaf size resulted from an effect on leaf area, but not leaf length, demonstrating changed leaf shape in response to [CO2]. Leaves exhibited a basipetal gradient of leaf development, investigated by defining seven interveinal areas, with growth ceasing first at the leaf tip. Interestingly, and in contrast to other reports, no spatial differences in epidermal cell size were apparent across the lamina, whereas a clear basipetal gradient in cell production rate was found. These data suggest that the rate and timing of cell production was more important in determining leaf shape, given the constant cell size across the leaf lamina. The effect of elevated [CO2] imposed on this developmental gradient suggested that leaf cell production continued longer in elevated [CO2] and that basal increases in cell production rate were also more important than altered cell expansion for increased final leaf size and altered leaf shape in elevated [CO2].

Bioactive films produced from self-assembling peptide amphiphiles as versatile substrates for tuning cell adhesion and tissue architecture in serum-free conditions

The development of versatile bioactive surfaces able to emulate in vivo conditions is of enormous importance to the future of cell and tissue therapy. Tuning cell behaviour on two-dimensional surfaces so that the cells perform as if they were in a natural three-dimensional tissue represents a significant challenge, but one that must be met if the early promise of cell and tissue therapy is to be fully realised. Due to the inherent complexities involved in the manufacture of biomimetic three-dimensional substrates, the scaling up of engineered tissue-based therapies may be simpler if based upon proven two-dimensional culture systems. In this work, we developed new coating materials composed of the self-assembling peptide amphiphiles (PAs) C16G3RGD (RGD) and C16G3RGDS (RGDS) shown to control cell adhesion and tissue architecture while avoiding the use of serum. When mixed with the C16ETTES diluent PA at 13 : 87 (mol mol-1) ratio at 1.25 times 10-3 M, the bioactive {PAs} were shown to support optimal adhesion, maximal proliferation, and prolonged viability of human corneal stromal fibroblasts ({hCSFs)}, while improving the cell phenotype. These {PAs} also provided stable adhesive coatings on highly-hydrophobic surfaces composed of striated polytetrafluoroethylene ({PTFE)}, significantly enhancing proliferation of aligned cells and increasing the complexity of the produced tissue. The thickness and structure of this highly-organised tissue were similar to those observed in vivo, comprising aligned newly-deposited extracellular matrix. As such, the developed coatings can constitute a versatile biomaterial for applications in cell biology, tissue engineering, and regenerative medicine requiring serum-free conditions.

Comparative in situ analyses of cell wall matrix polysaccharide dynamics in developing rice and wheat grain

Cell wall polysaccharides of wheat and rice endosperm are an important source of dietary fibre. Monoclonal antibodies specific to cell wall polysaccharides were used to determine polysaccharide dynamics during the development of both wheat and rice grain. Wheat and rice grain present near synchronous developmental processes and significantly different endosperm cell wall compositions, allowing the localisation of these polysaccharides to be related to developmental changes. Arabinoxylan (AX) and mixed-linkage glucan (MLG) have analogous cellular locations in both species, with deposition of AX and MLG coinciding with the start of grain filling. A glucuronoxylan (GUX) epitope was detected in rice, but not wheat endosperm cell walls. Callose has been reported to be associated with the formation of cell wall outgrowths during endosperm cellularisation and xyloglucan is here shown to be a component of these anticlinal extensions, occurring transiently in both species. Pectic homogalacturonan (HG) was abundant in cell walls of maternal tissues of wheat and rice grain, but only detected in endosperm cell walls of rice in an unesterified HG form. A rhamnogalacturonan-I (RG-I) backbone epitope was observed to be temporally regulated in both species, detected in endosperm cell walls from 12 DAA in rice and 20 DAA in wheat grain. Detection of the LM5 galactan epitope showed a clear distinction between wheat and rice, being detected at the earliest stages of development in rice endosperm cell walls, but not detected in wheat endosperm cell walls, only in maternal tissues. In contrast, the LM6 arabinan epitope was detected in both species around 8 DAA and was transient in wheat grain, but persisted in rice until maturity.

Regulation of cell migration in cancer: investigation into the regulation of cancer cell blebbing in the extracellular matrix

Cell migration is a highly coordinated process and any aberration in the regulatory mechanisms could result in pathological conditions such as cancer. The ability of cancer cells to disseminate to distant sites within the body has made it difficult to treat. Cancer cells also exhibit plasticity that makes them able to interconvert from an elongated, mesenchymal morphology to an amoeboid blebbing form under different physiological conditions. Blebs are spherical membrane protrusions formed by actomyosin-mediated contractility of cortical actin resulting in increased hydrostatic pressure and subsequent detachment of the membrane from the cortex. Tumour cells use blebbing as an alternative mode of migration by squeezing through preexisting gaps in the ECM, and bleb formation is believed to be mediated by the Rho-ROCK signaling pathway. However, the involvement of transmembrane water and ion channels in cell blebbing has not been examined. In the present study, the role of the transmembrane water channels, aquaporins, transmembrane ion transporters and lipid signaling enzymes in the regulation of blebbing was investigated. Using 3D matrigel matrix as an in vitro model to mimic normal extracellular matrix, and a combination of confocal and time-lapse microscopy, it was found that AQP1 knockdown by siRNA ablated blebbing of HT1080 and ACHN cells, and overexpression of AQP1-GFP not only significantly increased bleb size with a corresponding decrease in bleb numbers, but also induced bleb formation in non-blebbing cell lines. Importantly, AQP1 overexpression reduces bleb lifespan due to faster bleb retraction. This novel finding of AQP1-facilitated bleb retraction requires the activity of the Na+/H+ pump as inhibition of the ion transporter, which was found localized to intracellular vesicles, blocked bleb retraction in both cell lines. This study also demonstrated that a differential regulation of cell blebbing by AQP isoforms exists as knockdown of AQP5 had no effect on bleb formation. Data from this study also demonstrates that the lipid signaling PLD2 signals through PA in the LPA-LPAR-Rho-ROCK axis to positively regulate bleb formation in both cell lines. Taken together, this work provides a novel role of AQP1 and Na+/H+ pump in regulation of cell blebbing, and this could be exploited in the development of new therapy to treat cancer.

Crop yield reduction in the tropics under climate change: Processes and uncertainties

Many modelling studies examine the impacts of climate change on crop yield, but few explore either the underlying bio-physical processes, or the uncertainty inherent in the parameterisation of crop growth and development. We used a perturbed-parameter crop modelling method together with a regional climate model (PRECIS) driven by the 2071-2100 SRES A2 emissions scenario in order to examine processes and uncertainties in yield simulation. Crop simulations used the groundnut (i.e. peanut; Arachis hypogaea L.) version of the General Large-Area Model for annual crops (GLAM). Two sets of GLAM simulations were carried out: control simulations and fixed-duration simulations, where the impact of mean temperature on crop development rate was removed. Model results were compared to sensitivity tests using 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., and Bell, M.J., 1995, A peanut simulation model: I. Model development and testing. Agron. J. 87, 1085-1093]. GLAM simulations were particularly sensitive to two processes. First, elevated vapour pressure deficit (VPD) consistently reduced yield. The same result was seen in some simulations using both other crop models. Second, GLAM crop duration was longer, and yield greater, when the optimal temperature for the rate of development was exceeded. Yield increases were also seen in one other crop model. Overall, the models differed in their response to super-optimal temperatures, and that difference increased with mean temperature; percentage changes in yield between current and future climates were as diverse as -50% and over +30% for the same input data. The first process has been observed in many crop experiments, whilst the second has not. Thus, we conclude that there is a need for: (i) more process-based modelling studies of the impact of VPD on assimilation, and (ii) more experimental studies at super-optimal temperatures. Using the GLAM results, central values and uncertainty ranges were projected for mean 2071-2100 crop yields in India. In the fixed-duration simulations, ensemble mean yields mostly rose by 10-30%. The full ensemble range was greater than this mean change (20-60% over most of India). In the control simulations, yield stimulation by elevated CO2 was more than offset by other processes-principally accelerated crop development rates at elevated, but sub-optimal, mean temperatures. Hence, the quantification of uncertainty can facilitate relatively robust indications of the likely sign of crop yield changes in future climates. (C) 2007 Elsevier B.V. All rights reserved.

Elimination of the Greenland Ice Sheet in a High CO2 Climate

Projections of future global sea level depend on reliable estimates of changes in the size of polar ice sheets. Calculating this directly from global general circulation models (GCMs) is unreliable because the coarse resolution of 100 km or more is unable to capture narrow ablation zones, and ice dynamics is not usually taken into account in GCMs. To overcome these problems a high-resolution (20 km) dynamic ice sheet model has been coupled to the third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3). A novel feature is the use of two-way coupling, so that climate changes in the GCM drive ice mass changes in the ice sheet model that, in turn, can alter the future climate through changes in orography, surface albedo, and freshwater input to the model ocean. At the start of the main experiment the atmospheric carbon dioxide concentration was increased to 4 times the preindustrial level and held constant for 3000 yr. By the end of this period the Greenland ice sheet is almost completely ablated and has made a direct contribution of approximately 7 m to global average sea level, causing a peak rate of sea level rise of 5 mm yr-1 early in the simulation. The effect of ice sheet depletion on global and regional climate has been examined and it was found that apart from the sea level rise, the long-term effect on global climate is small. However, there are some significant regional climate changes that appear to have reduced the rate at which the ice sheet ablates.

The physical properties of the atmosphere in the new Hadley Centre Global Environmental Model (HadGEM1). Part II: Aspects of variability and regional climate

The performance of the atmospheric component of the new Hadley Centre Global Environmental Model (HadGEM1) is assessed in terms of its ability to represent a selection of key aspects of variability in the Tropics and extratropics. These include midlatitude storm tracks and blocking activity, synoptic variability over Europe, and the North Atlantic Oscillation together with tropical convection, the Madden-Julian oscillation, and the Asian summer monsoon. Comparisons with the previous model, the Third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3), demonstrate that there has been a considerable increase in the transient eddy kinetic energy (EKE), bringing HadGEM1 into closer agreement with current reanalyses. This increase in EKE results from the increased horizontal resolution and, in combination with the improved physical parameterizations, leads to improvements in the representation of Northern Hemisphere storm tracks and blocking. The simulation of synoptic weather regimes over Europe is also greatly improved compared to HadCM3, again due to both increased resolution and other model developments. The variability of convection in the equatorial region is generally stronger and closer to observations than in HadCM3. There is, however, still limited convective variance coincident with several of the observed equatorial wave modes. Simulation of the Madden-Julian oscillation is improved in HadGEM1: both the activity and interannual variability are increased and the eastward propagation, although slower than observed, is much better simulated. While some aspects of the climatology of the Asian summer monsoon are improved in HadGEM1, the upper-level winds are too weak and the simulation of precipitation deteriorates. The dominant modes of monsoon interannual variability are similar in the two models, although in HadCM3 this is linked to SST forcing, while in HadGEM1 internal variability dominates. Overall, analysis of the phenomena considered here indicates that HadGEM1 performs well and, in many important respects, improves upon HadCM3. Together with the improved representation of the mean climate, this improvement in the simulation of atmospheric variability suggests that HadGEM1 provides a sound basis for future studies of climate and climate change.

A Nyström method for a boundary value problem arising in unsteady water wave problems

This paper is concerned with solving numerically the Dirichlet boundary value problem for Laplace’s equation in a nonlocally perturbed half-plane. This problem arises in the simulation of classical unsteady water wave problems. The starting point for the numerical scheme is the boundary integral equation reformulation of this problem as an integral equation of the second kind on the real line in Preston et al. (2008, J. Int. Equ. Appl., 20, 121–152). We present a Nystr¨om method for numerical solution of this integral equation and show stability and convergence, and we present and analyse a numerical scheme for computing the Dirichlet-to-Neumann map, i.e., for deducing the instantaneous fluid surface velocity from the velocity potential on the surface, a key computational step in unsteady water wave simulations. In particular, we show that our numerical schemes are superalgebraically convergent if the fluid surface is infinitely smooth. The theoretical results are illustrated by numerical experiments.

Sequence analysis of the cupin gene family in Synechocystis PCC6803

The recently described cupin superfamily of proteins includes the germin and germinlike proteins, of which the cereal oxalate oxidase is the best characterized. This superfamily also includes seed storage proteins, in addition to several microbial enzymes and proteins with unknown function. All these proteins are characterized by the conservation of two central motifs, usually containing two or three histidine residues presumed to be involved with metal binding in the catalytic active site. The present study on the coding regions of Synechocystis PCC6803 identifies a previously unknown group of 12 related cupins, each containing the characteristic two-motif signature. This group comprises 11 single-domain proteins, ranging in length from 104 to 289 residues, and includes two phosphomannose isomerases and two epimerases involved in cell wall synthesis, a member of the pirin group of nuclear proteins, a possible transcriptional regulator, and a close relative-of a cytochrome c551 from Rhodococcus. Additionally, there is a duplicated, two-domain protein that has close similarity to an oxalate decarboxylase from the fungus Collybia velutipes and that is a putative progenitor of the storage proteins of land plants.

Use of global gene expression patterns in mechanistic studies of oestrogen action in MCF7 human breast cancer cells

Over the years, the MCF7 human breast cancer cell line has provided a model system for the study of cellular and molecular mechanisms in oestrogen regulation of cell proliferation and in progression to oestrogen and antioestrogen independent growth. Global gene expression profiling has shown that oestrogen action in MCF7 cells involves the coordinated regulation of hundreds of genes across a wide range of functional groupings and that more genes are down regulated than upregulated. Adaptation to long-term oestrogen deprivation, which results in loss of oestrogen-responsive growth, involves alterations to gene patterns not only at early time points (0-4 weeks) but continuing through to later times (20-55 weeks), and even involves alterations to patterns of oestrogen-regulated gene expression. Only 48% of the genes which were regulated >= 2-fold by oestradiol in oestrogen-responsive cells retained this responsiveness after long-term oestrogen deprivation but other genes developed de novo oestrogen regulation. Long-term exposure to fulvestrant, which resulted in loss of growth inhibition by the antioestrogen, resulted in some very large fold changes in gene expression up to 10,000-fold. Comparison of gene profiles produced by environmental chemicals with oestrogenic properties showed that each ligand gave its own unique expression profile which suggests that environmental oestrogens entering the human breast may give rise to a more complex web of interference in cell function than simply mimicking oestrogen action at inappropriate times. (C) 2009 Elsevier Ltd. All rights reserved.