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.

Flood detection in urban areas using TerraSAR-X

Flooding is a major hazard in both rural and urban areas worldwide, but it is in urban areas that the impacts are most severe. An investigation of the ability of high resolution TerraSAR-X Synthetic Aperture Radar (SAR) data to detect flooded regions in urban areas is described. The study uses a TerraSAR-X image of a 1 in 150 year flood near Tewkesbury, UK, in 2007, for which contemporaneous aerial photography exists for validation. The DLR SAR End-To-End simulator (SETES) was used in conjunction with airborne scanning laser altimetry (LiDAR) data to estimate regions of the image in which water would not be visible due to shadow or layover caused by buildings and taller vegetation. A semi-automatic algorithm for the detection of floodwater in urban areas is described, together with its validation using the aerial photographs. 76% of the urban water pixels visible to TerraSAR-X were correctly detected, with an associated false positive rate of 25%. If all urban water pixels were considered, including those in shadow and layover regions, these figures fell to 58% and 19% respectively. The algorithm is aimed at producing urban flood extents with which to calibrate and validate urban flood inundation models, and these findings indicate that TerraSAR-X is capable of providing useful data for this purpose.

Construction of coordination polymers of cadmium(II) with mixed hexamethylenetetramine and terephthalate or thiocyanate ligands

Three new polymeric complexes [Cd(hmt)(SCN)(2)(H2O)(2)](n) (1), [Cd-3(mu(2)-hmt)(2)(SCN)(6)(H2O)(2)](n) (2), and [Cd-2(hmt)(2)(tP)(2)(H2O)(6)](n) (3) [hmt = hexamethylenetetramine, tp = terephthalate] have been synthesized and characterized by single crystal X-ray diffraction. Both the compounds 1 and 2 are 1-D polymers where Cd units are linked by double end-to-end thiocyanate bridges but in 2 the chain is wider containing three cadmium atoms instead of one as found in 1. In both compounds the coordination environment around cadmium atom is distorted octahedral. Compound 3 is a three-dimensional polymer having water-filled microporous channels. Both tp and brut are mu(2)-bridged. One of the acid groups of tp is coordinated in chelating bidentate and the other in monodentate fashion. Half of its Cd atoms are hexa-coordinated and the rest are hepta-coordinated. Thermogravimetric analysis and X-ray diffraction study of 3 show that its framework remains intact upon removal of water molecules. The flexibility of coordination number around cadmium atoms (six or seven) probably plays an important role in establishing the rigidity of the framework. (C) 2003 Elsevier Ltd. All rights reserved.

Facile synthesis of Cu(II) complexes of monocondensed N,N,N donor Schiff base ligands: crystal structure, spectroscopic and magnetic properties

Four new copper(II) complexes, [((CuLN3)-N-1)(2)](ClO4)(2) (1), [(CuL2 N-3)(2)](ClO4)(2) (2), [CuL3(N-3)ClO4)](n) (3) and [CuL4(mu-1,1-N-3)(mu-1,3-N-3)(ClO4)](n) (4) where L-1 = N-1-pyridin-2-yl-methylene-propane-1,3-diamine, L-2 = N-1-(1-pyridin-2-yl-ethylidene)propane-1,3-diamine, L-3 =N-1-(1-pyridin-2-yl-ethylidene)ethane-1,2-diamine and L-4=N-1-(1-pyridin-2-yl-ethylidene)propane-1,2-diamine are four tridentate N,N,N donor Schiff base ligands, have been derived and structurally characterized by X-ray crystallography. Compounds 1 and 2 consist of double basal-apical end-on (EO) azide bridged dinuclear Cu-II complexes with square-pyramidal geometry. In complex 3 the square planar mononuclear [CuL3 (N-3)] units are linked by weakly coordinated perchlorate ions in the axial positions of Cu-II to form a one-dimensional chain. Two such chains are connected by hydrogen bonds involving perchlorate ions and azide groups. Compound 4 consists of 1-D chains in which the Cu-II ions with a square-pyramidal geometry are alternately bridged by single EO and end-to-end (EE) azido ligands, both adopting a basal-apical disposition. Variable temperature (300-2 K) magnetic susceptibility measurements and magnetization measurements at 2 K have been performed. The results reveal that complexes 1 and 2 are antiferromagnetically coupled through azido bridges (J= -12.18 +/- 0.09 and -4.43 +/- 0.1 cm(-1) for 1 and 2, respectively). Complex 3 shows two different magnetic interactions through the two kinds of hydrogen bonds; one is antiferromagnetic (J(1) = - 9.69 +/- 0.03 cm(-1)) and the other is ferromagnetic (J(2) = 1.00 +/- 0.01 cm(-1)). From a magnetic point of view complex 4 is a ferromagnetic dinuclear complex (J= 1.91 +/- 0.01 cm(-1)) coupled through the EO bridge only. The coupling through the EE bridge is practically nil as the N(azido)-Cu-II (axial) distance (2.643 angstrom) is too long. (C) 2006 Elsevier Ltd. All rights reserved.

Syntheses, structural analyses, and magneto-structural correlations of three polymeric Fe(II) complexes with azide ligand

Three new metal-organic polymeric complexes, [Fe(N-3)(2)(bPP)(2)] (1), [Fe(N-3)(2)(bpe)] (2), and [Fe(N-3)(2)(phen)] (3) [bpp = (1,3-bis(4-pyridyl)-propane), bpe = (1,2-bis(4-pyridyl)-ethane), phen = 1,10-phenanthroline], have been synthesized and characterized by single-crystal X-ray diffraction studies and low-temperature magnetic measurements in the range 300-2 K. Complexes 1 and 2 crystallize in the monoclinic system, space group C2/c, with the following cell parameters: a = 19.355(4) Angstrom, b = 7.076(2) Angstrom, c = 22.549(4) Angstrom, beta = 119.50(3)degrees, Z = 4, and a = 10.007(14) Angstrom, b = 13.789(18) Angstrom, c = 10.377(14) Angstrom, beta = 103.50(1)degrees, Z = 4, respectively. Complex 3 crystallizes in the triclinic system, space group P (1) over bar, with a = 7.155(12) Angstrom, b = 10.066(14) Angstrom, c = 10.508(14) Angstrom, alpha = 109.57(1)degrees, beta = 104.57(1)degrees, gamma = 105.10(1)degrees, and Z = 2. All coordination polymers exhibit octahedral Fe(II) nodes. The structural determination of 1 reveals a parallel interpenetrated structure of 2D layers of (4,4) topology, formed by Fe(II) nodes linked through bpp ligands, while mono-coordinated azide anions are pendant from the corrugated sheet. Complex 2 has a 2D arrangement constructed through 1D double end-to-end azide bridged iron(11) chains interconnected through bpe ligands. Complex 3 shows a polymeric arrangement where the metal ions are interlinked through pairs of end-on and end-to-end azide ligands exhibiting a zigzag arrangement of metals (Fe-Fe-Fe angle of 111.18degrees) and an intermetallic separation of 3.347 Angstrom (through the EO azide) and of 5.229 Angstrom (EE azide). Variable-temperature magnetic susceptibility data suggest that there is no magnetic interaction between the metal centers in 1, whereas in 2 there is an antiferromagnetic interaction through the end-to-end azide bridge. Complex 3 shows ferro- as well as anti-ferromagnetic interactions between the metal centers generated through the alternating end-on and end-to-end azide bridges. Complex I has been modeled using the D parameter (considering distorted octahedral Fe(II) geometry and with any possible J value equal to zero) and complex 2 has been modeled as a one-dimensional system with classical and/or quantum spin where we have used two possible full diagonalization processes: without and with the D parameter, considering the important distortions of the Fe(II) ions. For complex 3, the alternating coupling model impedes a mathematical solution for the modeling as classical spins. With quantum spin, the modeling has been made as in 2.

Mono, di and polynuclear Cu(II)-azido complexes incorporating N,N,N reduced Schiff base: syntheses, structure and magnetic behavior

Three kinds of copper(II) azide complexes have been synthesised in excellent yields by reacting Cu(ClO4)(2) . 6H(2)O with N,N-bis(2-pyridylmethyl)amine (L-1); N-(2-pyridylmethyl)-N',N'-dimethylethylenediamine (L-2); and N-(2-pyridylmethyl)-N',N'-diethylethylenediamine (L-3), respectively, in the presence of slight excess of sodium azide. They are the monomeric Cu(L-1)(N-3)(ClO4) (1), the end-to-end diazido-bridged Cu-2(L-2)(2)(mu-1,3-N-3)(2)(ClO4)(2) (2) and the single azido-bridged (mu-1,3-) 1D chain [Cu(L-3)(mu-1,3-N-3)](n)(ClO4)(n) (3). The crystal and molecular structures of these complexes have been solved. The variable temperature magnetic moments of type 2 and type 3 complexes were studied. Temperature dependent susceptibility for 2 was fitted using the Bleaney-Bowers expression which led to the parameters J = -3.43 cm(-1) and R = 1 X 10(-5). The magnetic data for 3 were fitted to Baker's expression for S = 1/2 and the parameters obtained were J = 1.6 cm(-1) and R = 3.2 x 10(-4). Crystal data are as follows. Cu(L-1)(N-3)(ClO4): Chemical formula, C12H13ClN6O4Cu; crystal system, monoclinic; space group, P2(1)/c; a = 8.788(12), b = 13.045(15), c = 14.213(15) Angstrom; beta = 102.960(10)degrees; Z = 4. Cu(L-2)(mu-N-3)(ClO4): Chemical formula. C10H17ClN6O4Cu: crystal system, monoclinic; space group, P2(1)/c; a = 10.790(12), b = 8.568(9), c = 16.651(17) Angstrom; beta = 102.360(10)degrees; Z = 4. [Cu(L-3)(mu-N-3)](ClO4): Chemical formula, C12H21ClN6O4Cu; crystal system, monoclinic; space group, P2(1)/c; a = 12.331(14), b = 7.804(9), c = 18.64(2) Angstrom; beta = 103.405(10)degrees; Z = 4. (C) 2004 Elsevier B.V. All rights reserved.

A spectral performance model for the High Resolution Dynamics Limb Sounder (6-18µm)

A spectral performance model, designed to simulate the system spectral throughput for each of the 21 channels in the HIRDLS radiometer, is described. This model uses the measured spectral characteristics of each of the components in the optical train, appropriately corrected for their optical environment, to determine the end-to-end spectral throughput profile for each channel. This profile is then combined with the predicted thermal emission from the atmosphere, arising from the height of interest, to establish an in-band (wanted) to out-of-band (unwanted) radiance ratio. The results from the use of the model demonstrate that the instrument level radiometric requirements for the instrument will be achieved. The optical arrangement and spectral design requirements for filtering in the HIRDLS instrument are described together with a presentation of the performance achieved for the complete set of manufactured filters. Compliance of the predicted passband throughput model to the spectral positioning requi rements of the instrument is also demonstrated.

Clustered functional MRI of overt speech production

To investigate the neural network of overt speech production, eventrelated fMRI was performed in 9 young healthy adult volunteers. A clustered image acquisition technique was chosen to minimize speechrelated movement artifacts. Functional images were acquired during the production of oral movements and of speech of increasing complexity (isolated vowel as well as monosyllabic and trisyllabic utterances). This imaging technique and behavioral task enabled depiction of the articulo-phonologic network of speech production from the supplementary motor area at the cranial end to the red nucleus at the caudal end. Speaking a single vowel and performing simple oral movements involved very similar activation of the corticaland subcortical motor systems. More complex, polysyllabic utterances were associated with additional activation in the bilateral cerebellum,reflecting increased demand on speech motor control, and additional activation in the bilateral temporal cortex, reflecting the stronger involvement of phonologic processing.

GNSS occultation sounding for climate monitoring

Considerable efforts are currently invested into the setup of a Global Climate Observing System (GCOS) for monitoring climate change over the coming decades, which is of high relevance given concerns on increasing human influences. A promising potential contribution to the GCOS is a suite of spaceborne Global Navigation Satellite System (GNSS) occultation sensors for global long-term monitoring of atmospheric change in temperature and other variables with high vertical resolution and accuracy. Besides the great importance with respect to climate change, the provision of high quality data is essential for the improvement of numerical weather prediction and for reanalysis efforts. We review the significance of GNSS radio occultation sounding in the climate observations context. In order to investigate the climate change detection capability of GNSS occultation sensors, we are currently performing an end-to-end GNSS occultation observing system simulation experiment over the 25-year period 2001 to 2025. We report on this integrated analysis, which involves in a realistic manner all aspects from modeling the atmosphere via generating a significant set of stimulated measurements to an objective statistical analysis and assessment of 2001–2025 temporal trends.

Maximin SNR beamforming strategies for two-way relay channels

In this letter, we consider beamforming strategies in amplified-and-forward (AF) two-way relay channels, where two terminals and the relay are equipped with multiple antennas. Our aim is to optimize the worse end-to-end signal-to-noise ratio of the two links so that the reliability of both terminals can be guaranteed. We show that the optimization problem can be recast as a generalized fractional programing and be solved by using the Dinkelbach-type procedure combined with semidefinite programming. Simulation results confirm the efficiency of the proposed strategies.

Mobility based energy efficient and multi-sink algorithms for consumer home networks

With the fast development of the Internet, wireless communications and semiconductor devices, home networking has received significant attention. Consumer products can collect and transmit various types of data in the home environment. Typical consumer sensors are often equipped with tiny, irreplaceable batteries and it therefore of the utmost importance to design energy efficient algorithms to prolong the home network lifetime and reduce devices going to landfill. Sink mobility is an important technique to improve home network performance including energy consumption, lifetime and end-to-end delay. Also, it can largely mitigate the hot spots near the sink node. The selection of optimal moving trajectory for sink node(s) is an NP-hard problem jointly optimizing routing algorithms with the mobile sink moving strategy is a significant and challenging research issue. The influence of multiple static sink nodes on energy consumption under different scale networks is first studied and an Energy-efficient Multi-sink Clustering Algorithm (EMCA) is proposed and tested. Then, the influence of mobile sink velocity, position and number on network performance is studied and a Mobile-sink based Energy-efficient Clustering Algorithm (MECA) is proposed. Simulation results validate the performance of the proposed two algorithms which can be deployed in a consumer home network environment.

Two-way CSI-assisted AF relaying with HPA nonlinearity

In this paper, we investigate half-duplex two-way dual-hop channel state information (CSI)-assisted amplify-and-forward (AF) relaying in the presence of high-power amplifier (HPA) nonlinearity at relays. The expression for the end-to-end signal-to-noise ratio (SNR) is derived as per the modified system model by taking into account the interference caused by relaying scheme and HPA nonlinearity. The system performance of the considered relaying network is evaluated in terms of average symbol error probability (SEP) in Nakagami-$m$ fading channels, by making use of the moment-generating function (MGF) approach. Numerical results are provided and show the effects of several parameters, such as quadrature amplitude modulation (QAM) order, number of relays, HPA parameters, and Nakagami parameter, on performance.

Tracking performance evaluation on PETS 2015 Challenge datasets

This paper presents a quantitative evaluation of a tracking system on PETS 2015 Challenge datasets using well-established performance measures. Using the existing tools, the tracking system implements an end-to-end pipeline that include object detection, tracking and post- processing stages. The evaluation results are presented on the provided sequences of both ARENA and P5 datasets of PETS 2015 Challenge. The results show an encouraging performance of the tracker in terms of accuracy but a greater tendency of being prone to cardinality error and ID changes on both datasets. Moreover, the analysis show a better performance of the tracker on visible imagery than on thermal imagery.

Large deviations of Rouse polymer chain: first passage problem

The purpose of this paper is to investigate several analytical methods of solving first passage (FP) problem for the Rouse model, a simplest model of a polymer chain. We show that this problem has to be treated as a multi-dimensional Kramers' problem, which presents rich and unexpected behavior. We first perform direct and forward-flux sampling (FFS) simulations, and measure the mean first-passage time $\tau(z)$ for the free end to reach a certain distance $z$ away from the origin. The results show that the mean FP time is getting faster if the Rouse chain is represented by more beads. Two scaling regimes of $\tau(z)$ are observed, with transition between them varying as a function of chain length. We use these simulations results to test two theoretical approaches. One is a well known asymptotic theory valid in the limit of zero temperature. We show that this limit corresponds to fully extended chain when each chain segment is stretched, which is not particularly realistic. A new theory based on the well known Freidlin-Wentzell theory is proposed, where dynamics is projected onto the minimal action path. The new theory predicts both scaling regimes correctly, but fails to get the correct numerical prefactor in the first regime. Combining our theory with the FFS simulations lead us to a simple analytical expression valid for all extensions and chain lengths. One of the applications of polymer FP problem occurs in the context of branched polymer rheology. In this paper, we consider the arm-retraction mechanism in the tube model, which maps exactly on the model we have solved. The results are compared to the Milner-McLeish theory without constraint release, which is found to overestimate FP time by a factor of 10 or more.

Finality revived: powers and intentionality

Proponents of physical intentionality argue that the classic hallmarks of intentionality highlighted by Brentano are also found in purely physical powers. Critics worry that this idea is metaphysically obscure at best, and at worst leads to panpsychism or animism. I examine the debate in detail, finding both confusion and illumination in the physical intentionalist thesis. Analysing a number of the canonical features of intentionality, I show that they all point to one overarching phenomenon of which both the mental and the physical are kinds, namely finality. This is the finality of ‘final causes’, the long-discarded idea of universal action for an end to which recent proponents of physical intentionality are in fact pointing whether or not they realise it. I explain finality in terms of the concept of specific indifference, arguing that in the case of the mental, specific indifference is realised by the process of abstraction, which has no correlate in the case of physical powers. This analysis, I conclude, reveals both the strength and weakness of rational creatures such as us, as well as demystifying (albeit only partly) the way in which powers work.