Running Climate Models On The NERC Cluster Grid Using G-Rex

Compute grids are used widely in many areas of environmental science, but there has been limited uptake of grid computing by the climate modelling community, partly because the characteristics of many climate models make them difficult to use with popular grid middleware systems. In particular, climate models usually produce large volumes of output data, and running them usually involves complicated workflows implemented as shell scripts. For example, NEMO (Smith et al. 2008) is a state-of-the-art ocean model that is used currently for operational ocean forecasting in France, and will soon be used in the UK for both ocean forecasting and climate modelling. On a typical modern cluster, a particular one year global ocean simulation at 1-degree resolution takes about three hours when running on 40 processors, and produces roughly 20 GB of output as 50000 separate files. 50-year simulations are common, during which the model is resubmitted as a new job after each year. Running NEMO relies on a set of complicated shell scripts and command utilities for data pre-processing and post-processing prior to job resubmission. Grid Remote Execution (G-Rex) is a pure Java grid middleware system that allows scientific applications to be deployed as Web services on remote computer systems, and then launched and controlled as if they are running on the user's own computer. Although G-Rex is general purpose middleware it has two key features that make it particularly suitable for remote execution of climate models: (1) Output from the model is transferred back to the user while the run is in progress to prevent it from accumulating on the remote system and to allow the user to monitor the model; (2) The client component is a command-line program that can easily be incorporated into existing model work-flow scripts. G-Rex has a REST (Fielding, 2000) architectural style, which allows client programs to be very simple and lightweight and allows users to interact with model runs using only a basic HTTP client (such as a Web browser or the curl utility) if they wish. This design also allows for new client interfaces to be developed in other programming languages with relatively little effort. The G-Rex server is a standard Web application that runs inside a servlet container such as Apache Tomcat and is therefore easy to install and maintain by system administrators. G-Rex is employed as the middleware for the NERC1 Cluster Grid, a small grid of HPC2 clusters belonging to collaborating NERC research institutes. Currently the NEMO (Smith et al. 2008) and POLCOMS (Holt et al, 2008) ocean models are installed, and there are plans to install the Hadley Centre’s HadCM3 model for use in the decadal climate prediction project GCEP (Haines et al., 2008). The science projects involving NEMO on the Grid have a particular focus on data assimilation (Smith et al. 2008), a technique that involves constraining model simulations with observations. The POLCOMS model will play an important part in the GCOMS project (Holt et al, 2008), which aims to simulate the world’s coastal oceans. A typical use of G-Rex by a scientist to run a climate model on the NERC Cluster Grid proceeds as follows :(1) The scientist prepares input files on his or her local machine. (2) Using information provided by the Grid’s Ganglia3 monitoring system, the scientist selects an appropriate compute resource. (3) The scientist runs the relevant workflow script on his or her local machine. This is unmodified except that calls to run the model (e.g. with “mpirun”) are simply replaced with calls to "GRexRun" (4) The G-Rex middleware automatically handles the uploading of input files to the remote resource, and the downloading of output files back to the user, including their deletion from the remote system, during the run. (5) The scientist monitors the output files, using familiar analysis and visualization tools on his or her own local machine. G-Rex is well suited to climate modelling because it addresses many of the middleware usability issues that have led to limited uptake of grid computing by climate scientists. It is a lightweight, low-impact and easy-to-install solution that is currently designed for use in relatively small grids such as the NERC Cluster Grid. A current topic of research is the use of G-Rex as an easy-to-use front-end to larger-scale Grid resources such as the UK National Grid service.

Coriolis perturbations in the infrared spectrum of ethylene

Rotational structure has been resolved and analyzed in two of the infrared‐active perpendicular bands of C2H4 vapor: the Type b fundamental band, ν10, at 826 cm—1, and the Type c fundamental band, ν7, at 949 cm—1. Many of the individual PP and RR branch lines have been observed. The analysis has been confined to values of the quantum number K≥3, for which energy levels ethylene shows no detectable deviations from a symmetric‐top rotational structure. The analysis reveals a Coriolis interaction between ν7 and ν10, and between ν4 and ν10, and values of the Coriolis constants ζ7,10z and ζ4,10y are obtained; these are related to normal coordinate calculations for the appropriate symmetry species, and force constants are derived to fit the observed zeta constants. The band center of ν10 has been revised from the original figure of 810 cm—1 to the new value, 826 cm—1, and the inactive frequency ν4 is estimated to lie at 1023±3 cm—1, in good agreement with the previous estimate of 1027 cm—1. The change in the value of ν10 leads to a suggested change in the value of the Raman‐active fundamental ν6 from 1236 to 1222 cm—1. New combination bands have been observed at 2174 cm—1, assigned as ν3+ν10; and at 2252 cm—1, assigned as ν4+ν6; also rotational structure has been resolved and analyzed in the ν6+ν10 band at 2048 cm—1. The new data obtained for the C2H4 molecule are summarized in Table XII, with all of the other data presently available on the vibrational and rotational constants.

Analytical potentials for triatomic molecules from spectroscopic data V. Application to HOX (X=F, Cl, Br, I)

Analytical potential energy functions are reported for HOX (X=F, Cl, Br, I). The surface for HOF predicts two metastable minima as well as the equilibrium configuration. These correspond to HFO (bent) and OHF (linear). Ab initio calculations performed for the HOF surface confirm these predictions. Comparisons are drawn between the two sets of results, and a vibrational analysis is undertaken for the hydrogen bonded OHF species. For HOCl, one further minimum is predicted, corresponding to HClO (bent), the parameters for which compare favourably with those reported from ab initio studies. In contrast, only the equilibrium configurations are predicted to be stable for HOBr and HOI.

A global climatology of wind–wave interaction

Generally, ocean waves are thought to act as a drag on the surface wind so that momentum is transferred downwards, from the atmosphere into the waves. Recent observations have suggested that when long wavelength waves, characteristic of remotely generated swell, propagate faster than the surface wind momentum can also be transferred upwards. This upward momentum transfer acts to accelerate the near-surface wind, resulting in a low-level wave-driven wind jet. Previous studies have suggested that the sign reversal of the momentum flux is well predicted by the inverse wave age, the ratio of the surface wind speed to the speed of the waves at the peak of the spectrum. ECMWF ERA-40 data has been used here to calculate the global distribution of the inverse wave age to determine whether there are regions of the ocean that are usually in the wind-driven wave regime and others that are generally in the wave-driven wind regime. The wind-driven wave regime is found to occur most often in the mid-latitude storm tracks where wind speeds are generally high. The wave-driven wind regime is found to be prevalent in the tropics where wind speeds are generally light and swell can propagate from storms at higher latitudes. The inverse wave age is also a useful indicator of the degree of coupling between the local wind and wave fields. The climatologies presented emphasise the non-equilibrium that exists between the local wind and wave fields and highlight the importance of swell in the global oceans.

Combinatorial peptide ligand libraries and plant proteomics: a winning strategy at a price

The mechanism of action and properties of a solid-phase ligand library made of hexapeptides (combinatorial peptide ligand libraries or CPLL, for capturing the "hidden proteome", i.e. the low- and very low-abundance proteins Constituting the vast majority of species in any proteome. as applied to plant tissues, are reviewed here. Plant tissues are notoriously recalcitrant to protein extraction and to proteome analysis, Firstly, rigid plant cell walls need to be mechanically disrupted to release the cell content and, in addition to their poor protein yield, plant tissues are rich in proteases and oxidative enzymes, contain phenolic Compounds, starches, oils, pigments and secondary metabolites that massively contaminate protein extracts. In addition, complex matrices of polysaccharides, including large amount of anionic pectins, are present. All these species compete with the binding of proteins to the CPLL beads, impeding proper capture and identification I detection of low-abundance species. When properly pre-treated, plant tissue extracts are amenable to capture by the CPLL beads revealing thus many new species among them low-abundance proteins. Examples are given on the treatment of leaf proteins, of corn seed extracts and of exudate proteins (latex from Hevea brasiliensis). In all cases, the detection of unique gene products via CPLL Capture is at least twice that of control, untreated sample. (c) 2008 Elsevier B.V. All rights reserved.

Designing Hospital Bedside Systems for use by Patients with Visual Impairments

This paper describes the development of an interface to a hospital portal system for information, communication and entertainment such that it can be used easily and effectively by all patients regardless of their age, disability, computer experience or native language. Specifically, this paper reports on the work conducted to ensure that the interface design took into account the needs of visually impaired users.

H∞ robust control design for unsupported paraplegic standing: experimental evaluation

This paper is concerned with the design of robust feedback H~-control systems for the control of the upright posture of paraplegic persons standing. While the subject stands in a special apparatus, stabilising torque at the ankle joint is generated by electrical stimulation of the paralyzed calf muscles. Since the muscles acting as actuators in this setup show a significant degree of nonlinearity, a robust H~-control design is used. The design approach is implemented in experiments with a paraplegic subject. The results demonstrate good performance and closed loop stability over the whole range of operation.

Evidence for involvement of both IKCa and SKCa channels in hyperpolarizing responses of the rat middle cerebral artery

Endothelium-derived hyperpolarizing factor responses in the rat middle cerebral artery are blocked by inhibiting IKCa channels alone, contrasting with peripheral vessels where block of both IKCa and SKCa is required. As the contribution of IKCa and SKCa to endothelium-dependent hyperpolarization differs in peripheral arteries, depending on the level of arterial constriction, we investigated the possibility that SKCa might contribute to equivalent hyperpolarization in cerebral arteries under certain conditions. METHODS: Rat middle cerebral arteries (approximately 175 microm) were mounted in a wire myograph. The effect of KCa channel blockers on endothelium-dependent responses to the protease-activated receptor 2 agonist, SLIGRL (20 micromol/L), were then assessed as simultaneous changes in tension and membrane potential. These data were correlated with the distribution of arterial KCa channels revealed with immunohistochemistry. RESULTS: SLIGRL hyperpolarized and relaxed cerebral arteries undergoing variable levels of stretch-induced tone. The relaxation was unaffected by specific inhibitors of IKCa (TRAM-34, 1 micromol/L) or SKCa (apamin, 50 nmol/L) alone or in combination. In contrast, the associated smooth-muscle hyperpolarization was inhibited, but only with these blockers in combination. Blocking nitric oxide synthase (NOS) or guanylyl cyclase evoked smooth-muscle depolarization and constriction, with both hyperpolarization and relaxation to SLIGRL being abolished by TRAM-34 alone, whereas apamin had no effect. Immunolabeling showed SKCa and IKCa within the endothelium. CONCLUSIONS: In the absence of NO, IKCa underpins endothelium-dependent hyperpolarization and relaxation in cerebral arteries. However, when NOS is active SKCa contributes to hyperpolarization, whatever the extent of background contraction. These changes may have relevance in vascular disease states where NO release is compromised and when the levels of SKCa expression may be altered.

Diabatic processes modifying potential vorticity in a North Atlantic cyclone

Potential vorticity (PV) succinctly describes the evolution of large-scale atmospheric flow because of its material conservation and invertibility properties. However, diabatic processes in extratropical cyclones can modify PV and influence both mesoscale weather and the evolution of the synoptic-scale wave pattern. In this investigation, modification of PV by diabatic processes is diagnosed in a Met Office Unified Model (MetUM) simulation of a North Atlantic cyclone using a set of PV tracers. The structure of diabatic PV within the extratropical cyclone is investigated and linked to the processes responsible for it. On the mesoscale, a tripole of diabatic PV is generated across the tropopause fold extending down to the cold front. The structure results from a dipole in heating across the frontal interface due to condensation in the warm conveyor belt flanking the upper side of the fold and evaporation of precipitation in the dry intrusion and below. On isentropic surfaces intersecting the tropopause, positive diabatic PV is generated on the stratospheric side, while negative diabatic PV is generated on the tropospheric side. The stratospheric diabatic PV is generated primarily by long-wave cooling which peaks at the tropopause itself due to the sharp gradient in humidity there. The tropospheric diabatic PV originates locally from the long-wave radiation and non-locally by advection out of the top of heating associated with the large-scale cloud, convection and boundary layer schemes. In most locations there is no diabatic modification of PV at the tropopause itself but diabatic PV anomalies would influence the tropopause indirectly through the winds they induce and subsequent advection. The consequences of this diabatic PV dipole for the evolution of synoptic-scale wave patterns are discussed.

Rights, symbolic violence and the micro-politics of the rural: the case of the Parish Paths Partnership scheme

In this paper microlevel politics and conflict associated with social and economic change in the countryside and linked changes in rural governance are explored with a focus upon research carried out on a recent rural policy initiative aimed at local 'empowerment'. This acts as a touchstone for a wider theoretical discussion. The paper is theorised within a conceptual framework derived and extended from the work of Pierre Bourdieu and others in order to explore case studies of the English Countryside Commission's Parish Paths Partnership scheme. The micropolitics involved with this scheme are examined and used to highlight more general issues raised by increased 'parish empowerment' in the 'postrural'.

Contributions of roots and rootstocks to sustainable, intensified crop production

Sustainable intensification is seen as the main route for meeting the world’s increasing demands for food and fibre. As demands mount for greater efficiency in the use of resources to achieve this goal, so the focus on roots and rootstocks and their role in acquiring water and nutrients, and overcoming pests and pathogens, is increasing. The purpose of this review is to explore some of the ways in which understanding root systems and their interactions with soils could contribute to the development of more sustainable systems of intensive production. Physical interactions with soil particles limit root growth if soils are dense, but root–soil contact is essential for optimal growth and uptake of water and nutrients. X-ray microtomography demonstrated that maize roots elongated more rapidly with increasing root–soil contact, as long as mechanical impedance was not limiting root elongation, while lupin was less sensitive to changes in root–soil contact. In addition to selecting for root architecture and rhizosphere properties, the growth of many plants in cultivated systems is profoundly affected by selection of an appropriate rootstock. Several mechanisms for scion control by rootstocks have been suggested, but the causal signals are still uncertain and may differ between crop species. Linkage map locations for quantitative trait loci for disease resistance and other traits of interest in rootstock breeding are becoming available. Designing root systems and rootstocks for specific environments is becoming a feasible target.

Thr54 allele carriers of the Ala54Thr variant of FABP2 gene have associations with metabolic syndrome and hypertriglyceridemia in urban South Indians

The intestinal fatty acid-binding protein gene is proposed as a candidate gene for diabetes because the protein it codes is involved in fatty acid absorption and metabolism. This study investigates the association of the Ala54Thr variant of the intestinal fatty acid-binding protein gene on type 2 diabetes mellitus and other related metabolic traits in Asian Indians. Ala54Thr polymorphism was genotyped by using polymerase chain reaction-restriction fragment length polymorphism in unrelated 773 type 2 diabetic and 899 normal glucose-tolerant (NGT) subjects, randomly chosen from the Chennai Urban Rural Epidemiology Study, an ongoing population-based study in South India. The Ala54Thr polymorphism was not associated with type 2 diabetes mellitus or obesity. However, genotype-phenotype study revealed that the NGT subjects carrying the Thr54 allele had significantly higher 2-hour plasma glucose (P = .007), glycated hemoglobin (P = .004), 2-hour insulin (P = .027), and fasting low-density lipoprotein cholesterol (P = .032) levels compared with those with the Ala54 allele. Normal glucose-tolerant subjects with Ala54Thr and Thr54Thr genotypes had significantly higher fasting serum triglyceride levels (P = .003) compared with those with Ala54Ala. The subjects were stratified into those with hypertriglyceridemia (serum triglyceride levels >or=150 mg/dL) and those without. The odds ratio for hypertriglyceridemia for the individuals carrying the Ala54Thr genotype was 1.491 (95% confidence interval [CI], 1.22-1.83, P < .0001), and for those carrying the Thr54Thr genotype, it was 1.888 (95% CI, 1.34-2.67; P < .0001). Subjects were also stratified into those with metabolic syndrome (MS) and those without, according to modified Adult Treatment Panel III guidelines. The odds ratio (adjusted for age and sex) for MS for the individuals carrying the Ala54Thr genotype was 1.240 (95% CI, 1.02-1.51; P = .03), whereas for those carrying the Thr54Thr genotype, it was 1.812 (95% CI, 1.28-2.57; P = .001). Carriers of the Thr54 allele have associations with MS and hypertriglyceridemia in this urban South Indian population.

A continuum anisotropic model of sea-ice dynamics

We develop the essential ingredients of a new, continuum and anisotropic model of sea-ice dynamics designed for eventual use in climate simulation. These ingredients are a constitutive law for sea-ice stress, relating stress to the material properties of sea ice and to internal variables describing the sea-ice state, and equations describing the evolution of these variables. The sea-ice cover is treated as a densely flawed two-dimensional continuum consisting of a uniform field of thick ice that is uniformly permeated with narrow linear regions of thinner ice called leads. Lead orientation, thickness and width distributions are described by second-rank tensor internal variables: the structure, thickness and width tensors, whose dynamics are governed by corresponding evolution equations accounting for processes such as new lead generation and rotation as the ice cover deforms. These evolution equations contain contractions of higher-order tensor expressions that require closures. We develop a sea-ice stress constitutive law that relates sea-ice stress to the structure tensor, thickness tensor and strain rate. For the special case of empty leads (containing no ice), linear closures are adopted and we present calculations for simple shear, convergence and divergence.