Evolving graphs: dynamical models, inverse problems and propagation

Applications such as neuroscience, telecommunication, online social networking, transport and retail trading give rise to connectivity patterns that change over time. In this work, we address the resulting need for network models and computational algorithms that deal with dynamic links. We introduce a new class of evolving range-dependent random graphs that gives a tractable framework for modelling and simulation. We develop a spectral algorithm for calibrating a set of edge ranges from a sequence of network snapshots and give a proof of principle illustration on some neuroscience data. We also show how the model can be used computationally and analytically to investigate the scenario where an evolutionary process, such as an epidemic, takes place on an evolving network. This allows us to study the cumulative effect of two distinct types of dynamics.

Complexation of copper by sewage sludge-derived dissolved organic matter: Effects on soil sorption behaviour and plant uptake

The complexation of Cu by sewage sludge-derived dissolved organic matter (SSDOM) is a process by which the environmental significance of the element may become enhanced due to reduced soil sorption and, hence, increased mobility. The work described in this paper used an ion selective electrode procedure to show that SSDOM complexation of Cu was greatest at intermediate pH values because competition between hydrogen ions and Cu for SSDOM binding sites, and between hydroxyl ions and SSDOM as Cu ligands, was lowest at such values. Batch sorption experiments further showed that the process of Cu complexation by SSDOM provided an explanation for enhanced desorption of Cu from the solid phase of a contaminated, organic matter-rich, clay loam soil, and reduced adsorption of Cu onto the solid phase of a sandy loam soil. Complexation of Cu by SSDOM did not affect uptake of Cu by spring barley plants, when compared to free ionic Cu, in a sand-culture pot experiment. However, it did appear to lead to greater biomass yields of the plant; perhaps indicating that the Cu-SSDOM complex had a lower toxicity towards the plant than the free Cu ion.

Effects of carbon dioxide on the searching behaviour of the root-feeding clover weevil Sitona lepidus (Coleoptera : Curculionidae)

The respiratory emission of CO2 from roots is frequently proposed as an attractant that allows soil-dwelling insects to locate host plant roots, but this role has recently become less certain. CO2 is emitted from many sources other than roots, so does not necessarily indicate the presence of host plants, and because of the high density of roots in the upper soil layers, spatial gradients may not always be perceptible by soil-dwelling insects. The role of CO2 in host location was investigated using the clover root weevil Sitona lepidus Gyllenhall and its host plant white clover (Trifolium repens L.) as a model system. Rhizochamber experiments showed that CO2 concentrations were approximately 1000 ppm around the roots of white clover, but significantly decreased with increasing distance from roots. In behavioural experiments, no evidence was found for any attraction by S. lepidus larvae to point emissions of CO2, regardless of emission rates. Fewer than 15% of larvae were attracted to point emissions of CO2, compared with a control response of 17%. However, fractal analysis of movement paths in constant CO2 concentrations demonstrated that searching by S. lepidus larvae significantly intensified when they experienced CO2 concentrations similar to those found around the roots of white clover (i.e. 1000 ppm). It is suggested that respiratory emissions of CO2 may act as a 'search trigger' for S. lepidus, whereby it induces larvae to search a smaller area more intensively, in order to detect location cues that are more specific to their host plant.

Survival, Pb-uptake and behaviour of three species of earthworm in Pb treated soils determined using an OECD-style toxicity test and a soil avoidance test

Mature (clitellate) Eisenia andrei Bouche (ultra epigeic), Lumbricus rubellus Hoffmeister (epigeic), and Aporrectodea caliginosa (Savigny) (endogeic) earthworms were placed in soils treated with Pb(NO3)(2) to have concentrations in the range 1000 to 10 000 mg Pb kg(-1). After 28 days LC50(-95%confidence limit) (+95%confidence limit) values were E. andrei 5824(-361)(+898) mg Pb kg(-1), L. rubellus 2867(-193)(+145) mg Pb kg(-1) and A. caliginosa 2747(-304)(+239) mg Pb kg(-1) and EC50s for weight change were E. andrei 2841(-68)(+150) Pb kg(-1), L. rubellus 1303(-201)(+204) mg Pb kg(-1) and A. caliginosa 1208(-206)(+212) Mg Pb kg(-1). At any given soil Pb concentration, Pb tissue concentrations after 28 days were the same for all three earthworm species. In a soil avoidance test there was no difference between the behaviour of the different species. The lower sensitivity to Pb exhibited by E. andrei is most likely due to physiological adaptations associated with the modes of life of the earthworms, and could have serious implications for the use of this earthworm as the species of choice in standard toxicological testing. (c) 2005 Elsevier Ltd. All rights reserved.

An Assessment of the Recent Past and Future Climate Change, Glacier Retreat, and Runoff in the Caucasus Region Using Dynamical and Statistical Downscaling and HBV-ETH Hydrological Model

The paper discusses the observed and projected warming in the Caucasus region and its implications for glacier melt and runoff. A strong positive trend in summer air temperatures of 0.05 degrees C a(-1) is observed in the high-altitude areas providing for a strong glacier melt and continuous decline in glacier mass balance. A warming of 4-7 degrees C and 3-5 degrees C is projected for the summer months in 2071-2100 under the A2 and B2 emission scenarios respectively, suggesting that enhanced glacier melt can be expected. The expected changes in winter precipitation will not compensate for the summer melt and glacier retreat is likely to continue. However, a projected small increase in both winter and summer precipitation combined with the enhanced glacier melt will result in increased summer runoff in the currently glaciated region of the Caucasus (independent of whether the region is glaciated at the end of the twenty-first century) by more than 50% compared with the baseline period.

The Dubawnt Lake palaeo-ice stream: evidence for dynamic ics sheet behaviour on the Canadian Shield and insights regarding the controls on ice-stream location and vigour

We report evidence for a major ice stream that operated over the northwestern Canadian Shield in the Keewatin Sector of the Laurentide Ice Sheet during the last deglaciation 9000-8200 (uncalibrated) yr BP. It is reconstructed at 450 km in length, 140 km in width, and had an estimated catchment area of 190000 km. Mapping from satellite imagery reveals a suite of bedforms ('flow-set') characterized by a highly convergent onset zone, abrupt lateral margins, and where flow was presumed to have been fastest, a remarkably coherent pattern of mega-scale glacial lineations with lengths approaching 13 km and elongation ratios in excess of 40:1. Spatial variations in bedform elongation within the flow-set match the expected velocity field of a terrestrial ice stream. The flow pattern does not appear to be steered by topography and its location on the hard bedrock of the Canadian Shield is surprising. A soft sedimentary basin may have influenced ice-stream activity by lubricating the bed over the downstream crystalline bedrock, but it is unlikely that it operated over a pervasively deforming till layer. The location of the ice stream challenges the view that they only arise in deep bedrock troughs or over thick deposits of 'soft' fine-grained sediments. We speculate that fast ice flow may have been triggered when a steep ice sheet surface gradient with high driving stresses contacted a proglacial lake. An increase in velocity through calving could have propagated fast ice flow upstream (in the vicinity of the Keewatin Ice Divide) through a series of thermomechanical feedback mechanisms. It exerted a considerable impact on the Laurentide Ice Sheet, forcing the demise of one of the last major ice centres.

Chemical behaviour of the Wheal Jane bioremediation system

The effectiveness of remediation of the highly acidic and transition metal polluted mine water discharge from the Wheal Jane Mine by the Wheal Jane Passive Treatment Plant is described. The success of the remediation required that all the system components work as predicted. The study shows considerable success in the removal of key toxic metals and clearly demonstrates the potential for natural attenuation of acid mine drainage, particularly iron oxidation, by microbial populations. The Wheal Jane Passive Treatment Plant provides the only experimental facility of its kind. (C) 2004 Elsevier B.V. All rights reserved.

Convectively coupled equatorial waves in high resolution Hadley centre climate models

Current global atmospheric models fail to simulate well organised tropical phenomena in which convection interacts with dynamics and physics. A new methodology to identify convectively coupled equatorial waves, developed by NCAS-Climate, has been applied to output from the two latest models of the Met Office/Hadley Centre which have fundamental differences in dynamical formulation. Variability, horizontal and vertical structures, and propagation characteristics of tropical convection and equatorial waves, along with their coupled behaviour in the models are examined and evaluated against a previous comprehensive study of observations. It is shown that, in general, the models perform well for equatorial waves coupled with off-equatorial convection. However they perform poorly for waves coupled with equatorial convection. The vertical structure of the simulated wave is not conducive to energy conversion/growth and does not support the correct physical-dynamical coupling that occurs in the real world. The following figure shows an example of the Kelvin wave coupled with equatorial convection. It shows that the models fail to simulate a key feature of convectively coupled Kelvin wave in observations, namely near surface anomalous equatorial zonal winds together with intensified equatorial convection and westerly winds in phase with the convection. The models are also not able to capture the observed vertical tilt structure and the vertical propagation of the Kelvin wave into the lower stratosphere as well as the secondary peak in the mid-troposphere, particularly in HadAM3. These results can be used to provide a test-bed for experimentation to improve the coupling of physics and dynamics in climate and weather models.

Annular variability and eddy-zonal flow interactions in a simplified atmospheric GCM. Part 1: Characterization of high and low frequency behaviour

Experiments have been performed using a simplified, Newtonian forced, global circulation model to investigate how variability of the tropospheric jet can be characterized by examining the combined fluctuations of the two leading modes of annular variability. Eddy forcing of this variability is analyzed in the phase space of the leading modes using the vertically integrated momentum budget. The nature of the annular variability and eddy forcing depends on the time scale. At low frequencies the zonal flow and baroclinic eddies are in quasi equilibrium and anomalies propagate poleward. The eddies are shown primarily to reinforce the anomalous state and are closely balanced by the linear damping, leaving slow evolution as a residual. At high frequencies the flow is strongly evolving and anomalies are initiated on the poleward side of the tropospheric jet and propagate equatorward. The eddies are shown to drive this evolution strongly: eddy location and amplitude reflect the past baroclinicity, while eddy feedback on the zonal flow may be interpreted in terms of wave breaking associated with baroclinic life cycles in lateral shear.

Modelling the inorganic nitrogen behaviour in a small Mediterranean forested catchment, Fuirosos (Catalonia)

The aim of this work was to couple a nitrogen (N) sub-model to already existent hydrological lumped (LU4-N) and semi-distributed (LU4-R-N and SD4-R-N) conceptual models, to improve our understanding of the factors and processes controlling nitrogen cycling and losses in Mediterranean catchments. The N model adopted provides a simplified conceptualization of the soil nitrogen cycle considering mineralization, nitrification, immobilization, denitrification, plant uptake, and ammonium adsorption/desorption. It also includes nitrification and denitrification in the shallow perched aquifer. We included a soil moisture threshold for all the considered soil biological processes. The results suggested that all the nitrogen processes were highly influenced by the rain episodes and that soil microbial processes occurred in pulses stimulated by soil moisture increasing after rain. Our simulation highlighted the riparian zone as a possible source of nitrate, especially after the summer drought period, but it can also act as an important sink of nitrate due to denitrification, in particular during the wettest period of the year. The riparian zone was a key element to simulate the catchment nitrate behaviour. The lumped LU4-N model (which does not include the riparian zone) could not be validated, while both the semi-distributed LU4-R-N and SD4-R-N model (which include the riparian zone) gave satisfactory results for the calibration process and acceptable results for the temporal validation process.

A contour-advective semi-lagrangian numerical algorithm for simulating fine-scale conservative dynamical fields

This paper describes a novel numerical algorithm for simulating the evolution of fine-scale conservative fields in layer-wise two-dimensional flows, the most important examples of which are the earth's atmosphere and oceans. the algorithm combines two radically different algorithms, one Lagrangian and the other Eulerian, to achieve an unexpected gain in computational efficiency. The algorithm is demonstrated for multi-layer quasi-geostrophic flow, and results are presented for a simulation of a tilted stratospheric polar vortex and of nearly-inviscid quasi-geostrophic turbulence. the turbulence results contradict previous arguments and simulation results that have suggested an ultimate two-dimensional, vertically-coherent character of the flow. Ongoing extensions of the algorithm to the generally ageostrophic flows characteristic of planetary fluid dynamics are outlined.

Dynamical evolution of the 2003 southern hemisphere stratospheric winter using Envisat trace-gas observations

Data from the MIPAS instrument on Envisat, supplemented by meteorological analyses from ECMWF and the Met Office, are used to study the meteorological and trace-gas evolution of the stratosphere in the southern hemisphere during winter and spring 2003. A pole-centred approach is used to interpret the data in the physically meaningful context of the evolving stratospheric polar vortex. The following salient dynamical and transport features are documented and analysed: the merger of anticyclones in the stratosphere; the development of an intense, quasi-stationary anticyclone in spring; the associated top-down breakdown of the polar vortex; the systematic descent of air into the polar vortex; and the formation of a three-dimensional structure of a tracer filament on a planetary scale. The paper confirms and extends existing paradigms of the southern hemisphere vortex evolution. The quality of the MIPAS observations is seen to be generally good. though the water vapour retrievals are unrealistic above 10 hPa in the high-latitude winter.