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.

Can the dynamical impact of the stratosphere on the troposphere be described by large-scale adjustment to the stratospheric PV distribution?

Recent numerical experiments have demonstrated that the state of the stratosphere has a dynamical impact on the state of the troposphere. To account for such an effect, a number of mechanisms have been proposed in the literature, all of which amount to a large-scale adjustment of the troposphere to potential vorticity (PV) anomalies in the stratosphere. This paper analyses whether a simple PV adjustment suffices to explain the actual dynamical response of the troposphere to the state of the stratosphere, the actual response being determined by ensembles of numerical experiments run with an atmospheric general-circulation model. For this purpose, a new PV inverter is developed. It is shown that the simple PV adjustment hypothesis is inadequate. PV anomalies in the stratosphere induce, by inversion, flow anomalies in the troposphere that do not coincide spatially with the tropospheric changes determined by the numerical experiments. Moreover, the tropospheric anomalies induced by PV inversion are on a larger scale than the changes found in the numerical experiments, which are linked to the Atlantic and Pacific storm-tracks. These findings imply that the impact of the stratospheric state on the troposphere is manifested through the impact on individual synoptic-scale systems and their self-organization in the storm-tracks. Changes in these weather systems in the troposphere are not merely synoptic-scale noise on a larger scale tropospheric response, but an integral part of the mechanism by which the state of the stratosphere impacts that of the troposphere.

The impact of the stratosphere on the troposphere during the southern hemisphere stratospheric sudden warming, September 2002

Recent research has established that a small but statistically significant link exists between the stratosphere and the troposphere in the northern hemisphere extratropics. In this paper it is shown that a similar link exists between the stratosphere and troposphere during the unprecedented September 2002 sudden warming in the southern hemisphere. Two ensemble forecasts of the stratospheric sudden warming are run which have different stratospheric initial conditions and identical tropospheric initial conditions. Stratospheric initial conditions have an impact on the tropospheric flow at the peak of the major warming (5 days into the run) and on longer time-scales (18 days into the run). The character of this influence is a localized, equatorward shift of the tropospheric storm track. The averaged impact of the change in the position of the storm-track maps strongly onto the Southern Annular Mode structure, but does not have an annular character.

On the climate response of the low-latitude Pacific Ocean to changes in the global freshwater cycle

Under global warming, the predicted intensification of the global freshwater cycle will modify the net freshwater flux at the ocean surface. Since the freshwater flux maintains ocean salinity structures, changes to the density-driven ocean circulation are likely. A modified ocean circulation could further alter the climate, potentially allowing rapid changes, as seen in the past. The relevant feedback mechanisms and timescales are poorly understood in detail, however, especially at low latitudes where the effects of salinity are relatively subtle. In an attempt to resolve some of these outstanding issues, we present an investigation of the climate response of the low-latitude Pacific region to changes in freshwater forcing. Initiated from the present-day thermohaline structure, a control run of a coupled ocean-atmosphere general circulation model is compared with a perturbation run in which the net freshwater flux is prescribed to be zero over the ocean. Such an extreme experiment helps to elucidate the general adjustment mechanisms and their timescales. The atmospheric greenhouse gas concentrations are held constant, and we restrict our attention to the adjustment of the upper 1,000 m of the Pacific Ocean between 40°N and 40°S, over 100 years. In the perturbation run, changes to the surface buoyancy, near-surface vertical mixing and mixed-layer depth are established within 1 year. Subsequently, relative to the control run, the surface of the low-latitude Pacific Ocean in the perturbation run warms by an average of 0.6°C, and the interior cools by up to 1.1°C, after a few decades. This vertical re-arrangement of the ocean heat content is shown to be achieved by a gradual shutdown of the heat flux due to isopycnal (i.e. along surfaces of constant density) mixing, the vertical component of which is downwards at low latitudes. This heat transfer depends crucially upon the existence of density-compensating temperature and salinity gradients on isopycnal surfaces. The timescale of the thermal changes in the perturbation run is therefore set by the timescale for the decay of isopycnal salinity gradients in response to the eliminated freshwater forcing, which we demonstrate to be around 10-20 years. Such isopycnal heat flux changes may play a role in the response of the low-latitude climate to a future accelerated freshwater cycle. Specifically, the mechanism appears to represent a weak negative sea surface temperature feedback, which we speculate might partially shield from view the anthropogenically-forced global warming signal at low latitudes. Furthermore, since the surface freshwater flux is shown to play a role in determining the ocean's thermal structure, it follows that evaporation and/or precipitation biases in general circulation models are likely to cause sea surface temperature biases.

The influence of different artificial soil types on the acute toxicity of carbendazim to the earthworm Eisenia fetida in laboratory toxicity tests

Field populations of earthworms have shown a varied response in mortality to the fungicide carbendazim, the toxic reference substance used in agrochemical field trials. The aim of this study was to determine the influence of soil conditions as a potential cause of this variation. Laboratory acute toxicity tests were conducted using a range of artificial soils with varying soil components (organic matter, clay, pH and moisture). Batch adsorption/desorption studies were run to determine the influence of the soil properties on carbendazim behaviour. Adsorption was shown to be correlated with organic matter content and pH and this in turn could be linked to Eisenia fetida mortality, with lower mortality occurring with increased adsorption. Overall while E.fetida mortality did vary significantly between several of the soils the calculated LC50 values in the different soils did not cover a wide range (6.04-16.00 mg kg(-1)), showing that under these laboratory conditions soil components did not greatly influence carbendazim toxicity to E.fietida. (c) 2007 Elsevier Masson SAS. All rights reserved.

Impact of hydrographic data assimilation on the modelled Atlantic meridional overturning circulation

Here we make an initial step toward the development of an ocean assimilation system that can constrain the modelled Atlantic Meridional Overturning Circulation (AMOC) to support climate predictions. A detailed comparison is presented of 1° and 1/4° resolution global model simulations with and without sequential data assimilation, to the observations and transport estimates from the RAPID mooring array across 26.5° N in the Atlantic. Comparisons of modelled water properties with the observations from the merged RAPID boundary arrays demonstrate the ability of in situ data assimilation to accurately constrain the east-west density gradient between these mooring arrays. However, the presence of an unconstrained "western boundary wedge" between Abaco Island and the RAPID mooring site WB2 (16 km offshore) leads to the intensification of an erroneous southwards flow in this region when in situ data are assimilated. The result is an overly intense southward upper mid-ocean transport (0–1100 m) as compared to the estimates derived from the RAPID array. Correction of upper layer zonal density gradients is found to compensate mostly for a weak subtropical gyre circulation in the free model run (i.e. with no assimilation). Despite the important changes to the density structure and transports in the upper layer imposed by the assimilation, very little change is found in the amplitude and sub-seasonal variability of the AMOC. This shows that assimilation of upper layer density information projects mainly on the gyre circulation with little effect on the AMOC at 26° N due to the absence of corrections to density gradients below 2000 m (the maximum depth of Argo). The sensitivity to initial conditions was explored through two additional experiments using a climatological initial condition. These experiments showed that the weak bias in gyre intensity in the control simulation (without data assimilation) develops over a period of about 6 months, but does so independently from the overturning, with no change to the AMOC. However, differences in the properties and volume transport of North Atlantic Deep Water (NADW) persisted throughout the 3 year simulations resulting in a difference of 3 Sv in AMOC intensity. The persistence of these dense water anomalies and their influence on the AMOC is promising for the development of decadal forecasting capabilities. The results suggest that the deeper waters must be accurately reproduced in order to constrain the AMOC.

Impact of land cover changes on the South African climate

The Joint UK Land Environmental Simulator (JULES) was run offline to investigate the sensitivity of land surface type changes over South Africa. Sensitivity tests were made in idealised experiments where the actual land surface cover is replaced by a single homogeneous surface type. The vegetation surface types on which some of the experiments were made are static. Experimental tests were evaluated against the control. The model results show among others that the change of the surface cover results in changes of other variables such as soil moisture, albedo, net radiation and etc. These changes are also visible in the spin up process. The model shows different surfaces spinning up at different cycles. Because JULES is the land surface model of Unified Model, the results could be more physically meaningful if it is coupled to the Unified Model.

Rice versus fish revisited: On the integrated management of floodplain resources in Bangladesh

Disproportionately little attention has been paid to the dry season trade-off between rice and (inland capture) fish production on the floodplains of Bangladesh, compared to the same trade-off during the flood season. As the rural economy grows increasingly dominated by dry-season irrigated rice production, and floodplain land and water come under ever-increasing pressure during the dry winter months, there is an urgent need to focus attention on these dry months that are so critical to the survival and propagation of the floodplain resident fish, and to the poor people that depend on these fish for their livelihood. This article examines three important dry-season natural resource constraints to floodplain livelihoods in Bangladesh, and finds a common factor at the heart of all three: rice cultivation on lands at low and very low elevations. The article articulates the system interlinkages that bind these constraints and the long-run trend towards irrigated rice cropping on lower-lying lands, and suggests a management approach based on locally tailored strategies to arrest this trend. Apart from its direct relevance to the floodplains of Bangladesh, which support more than 100 million people, these lessons have relevance for river floodplain systems elsewhere in the developing world, notably the Mekong Delta.

Shrinking the posterior: A note on the Nerlovian model

Diebold and Lamb (1997) argue that since the long-run elasticity of supply derived from the Nerlovian model entails a ratio of random variables, it is without moments. They propose minimum expected loss estimation to correct this problem but in so-doing ignore the fact that a non white-noise-error is implicit in the model. We show that, as a consequence the estimator is biased and demonstrate that Bayesian estimation which fully accounts for the error structure is preferable.

The influence of larval density, food availability and habitat longevity on the life history and population growth rate of the midge Chironomus riparius

Despite long-standing interest in the forms and mechanisms of density dependence, these are still imperfectly understood. However, in a constant environment an increase in density must reduce per capita resource availability, which in turn leads to reduced survival, fecundity and somatic growth rate. Here we report two population experiments examining the density dependent responses under controlled conditions of an important indicator species, Chironomus riparius. The first experiment was run for 35 weeks and was started at low density with replicate populations being fed three different rations. Increased ration reduced generation time and increased population growth rate (pgr) but had no effect on survival, fecundity and female body weight in the first generation. In the second generation there was a six-fold increase in generation time, presumably due to the greatly reduced per capita resource availability as the estimated initial densities of the second generation were 300 times greater than the first. Juvenile survival to emergence, fecundity, adult body weight and pgr declined by 90%, 75%, 35% and 99%, respectively. These large between-generation effects may have obscured the effects of the threefold variation in ration, as only survival to emergence significantly increased with ration in the second generation. These results suggest that some chironomid larvae survive a reduction in resource availability by growing more slowly. In the ephemeral habitats sometimes occupied by C. riparius, the effects of population density may depend crucially on the longevity of the environment. A second experiment was therefore performed to measure pgr from six different starting densities over an eight-week period. The relationship between pgr and density was concave, viewed from above. At densities above 16 larvae per cm(2), less than 1% of the population emerged and no offspring were produced. Under the conditions of experiment 2 - an 8-week habitat lifespan carrying capacity was estimated as 8 larvae per cm(2).

In vitro examination of the effect of Orlistat on the ability of the faecal microbiota to utilize dietary lipids

Orlistat is an anti-obesity treatment with which several gastrointestinal (GI) side-effects are commonly associated in the initial stages of therapy. There is no physiological explanation as to why two-thirds of those who take the drug experience one or more side-effects. It has been hypothesized that the GI microbiota may protect from or contribute to these GI disturbances. Using in vitro batch culture and human gut model systems, studies were conducted to determine whether increased availability of dietary lipids and/or orlistat affect the composition and/or activity of the faecal microbiota. Results from 24-h batch culture fermentation experiments demonstrated no effect of orlistat in the presence or absence of a dietary lipid (olive oil) on the composition of bacterial communities [as determined by fluorescence in situ hybridization (FISH) and denaturing gradient gel electrophoresis (DGGE) analyses], but did show there was great variability in the lipolytic activities of the microbiotas of individuals, as determined by gas chromatography analysis of long-chain fatty acids in samples. Subsequent studies focused on the effect of orlistat in the presence and absence of lipid in in vitro human gut model systems. Systems were run for 14 days with gut model medium (GMM) only (to steady state, SS), then fed at 12-h intervals with 50 mg orlistat, 2 g olive oil or a mixture of both for 14 days. FISH and DGGE were used to monitor changes in bacterial populations. Bacteria were cultivated from the GMM only (control) systems at SS. All strains isolated were screened for lipolytic activity using tributyrin agar. FISH and DGGE demonstrated that none of the compounds (singly or in combination) added to the systems had any notable effect on microbial population dynamics for any of the donors, although Subdoligranulum populations appeared to be inhibited by orlistat in the presence or absence of lipid. Orlistat had little or no effect on the metabolism of indigenous and added lipids in the fermentation systems, but there was great variability in the way the faecal microbiotas of the donors were able to degrade added lipids. Variability in lipid degradation could be correlated with the number and activity of isolated lipolytic bacteria. The mechanism by which orlistat and the GI microbiota cause side-effects in individuals is unknown, but several hypotheses have been proposed to account for their manifestation. The demonstration of great variability in the lipolytic activity of microbiotas to degrade lipids led to a large-scale cultivation-based study of lipolytic/lipase-positive bacteria present in the human faecal microbiota. Of 4,000 colonies isolated from 15 donors using five different agars, 378 strains were identified that had lipase activity. Molecular identification of strains isolated from five donors demonstrated that lipase activity is more prevalent in the human GI microbiota than previously thought, with members of the phyla Firmicutes, Bacteroidetes and Actinobacteria identified. Molecular identification and characterization of the substrate specificities of the strains will be carried out as part of ongoing work.

The impact of the state of the troposphere on the response to stratospheric heating in a simplified GCM

Previous studies have made use of simplified general circulation models (sGCMs) to investigate the atmospheric response to various forcings. In particular, several studies have investigated the tropospheric response to changes in stratospheric temperature. This is potentially relevant for many climate forcings. Here the impact of changing the tropospheric climatology on the modeled response to perturbations in stratospheric temperature is investigated by the introduction of topography into the model and altering the tropospheric jet structure. The results highlight the need for very long integrations so as to determine accurately the magnitude of response. It is found that introducing topography into the model and thus removing the zonally symmetric nature of the model’s boundary conditions reduces the magnitude of response to stratospheric heating. However, this reduction is of comparable size to the variability in the magnitude of response between different ensemble members of the same 5000-day experiment. Investigations into the impact of varying tropospheric jet structure reveal a trend with lower-latitude/narrower jets having a much larger magnitude response to stratospheric heating than higher-latitude/wider jets. The jet structures that respond more strongly to stratospheric heating also exhibit longer time scale variability in their control run simulations, consistent with the idea that a feedback between the eddies and the mean flow is both responsible for the persistence of the control run variability and important in producing the tropospheric response to stratospheric temperature perturbations.

Lumbricus terrestris L. does not impact on the remediation efficiency of compost and biochar amendments

The aim of this study was to test the impact of compost and Biochar, with or without earthworms, on the mobility and availability of metals, and on the growth of grass to re-vegetate contaminated soil from the Parys Mountain mining site, Anglesey. We also determined if the addition of earthworms compromises remediation efforts. In a laboratory experiment, contaminated soil (1343 mg Cu kg−1, 2511 mg Pb kg−1 and 262 mg Zn kg−1) was remediated with compost and/or Biochar. After 77 days Lumbricus terrestris L. earthworms were added to the treatment remediated with both compost and Biochar, and left for 28 days. L. terrestris was not able to survive in the Biochar, compost or unamended treatments. A germination and growth bioassay, using Agrostis capillaris (Common Bent) was then run on all treatments for 28 days. The combination of Biochar and compost decreased water soluble Cu (from 5.6 to 0.2 mg kg−1), Pb (0.17 to less than 0.007 mg kg−1) and Zn (3.3 to 0.05 mg kg−1) in the contaminated soil and increased the pH from 2.7 to 6.6. The addition of L. terrestris to this treatment had no effect on the concentration of the water soluble metals in the remediated soil. The compost was the only treatment that resulted in germination and growth of A. capillaris suitable for re-vegetation purposes. However, the combination of compost, Biochar (with or without L. terrestris) produced the lowest concentrations of Cu (8 mg kg−1) and Zn (36 mg kg−1) in the aboveground biomass, lower than the compost treatment (15 mgCu kg−1 and 126 mgZn kg−1). The addition of Biochar and compost both separately and as co-amendments was effective in reducing the mobility and availability of metals. The addition of L. terrestris did not re-mobilise previously sequestered metals.

The underpricing of IPOs on the stock exchange of Mauritius

This paper investigates the underpricing of IPOs on the Stock Exchange of Mauritius (SEM). Taking into account the whole population of firms which went public since the inception of the SEM until 2010, the results show an average degree of underpricing within the range 10 to 20%. Using a regression approach, we demonstrate that the aftermarket risk level and auditor's reputation both have a significant positive impact on initial returns. We propose the use of the Z-score as a composite measure of a firm's ex ante financial strength, and find that it has a significant negative effect on the degree of short-run underpricing.

The effect of nonlinearity on the variational assimilation of satellite observations using a simple column model

Cloud imagery is not currently used in numerical weather prediction (NWP) to extract the type of dynamical information that experienced forecasters have extracted subjectively for many years. For example, rapidly developing mid-latitude cyclones have characteristic signatures in the cloud imagery that are most fully appreciated from a sequence of images rather than from a single image. The Met Office is currently developing a technique to extract dynamical development information from satellite imagery using their full incremental 4D-Var (four-dimensional variational data assimilation) system. We investigate a simplified form of this technique in a fully nonlinear framework. We convert information on the vertical wind field, w(z), and profiles of temperature, T(z, t), and total water content, qt (z, t), as functions of height, z, and time, t, to a single brightness temperature by defining a 2D (vertical and time) variational assimilation testbed. The profiles of w, T and qt are updated using a simple vertical advection scheme. We define a basic cloud scheme to obtain the fractional cloud amount and, when combined with the temperature field, we convert this information into a brightness temperature, having developed a simple radiative transfer scheme. With the exception of some matrix inversion routines, all our code is developed from scratch. Throughout the development process we test all aspects of our 2D assimilation system, and then run identical twin experiments to try and recover information on the vertical velocity, from a sequence of observations of brightness temperature. This thesis contains a comprehensive description of our nonlinear models and assimilation system, and the first experimental results.