Simulation modelling applied to runoff modelling using mapscript

Dynamic spatial analysis addresses computational aspects of space–time processing. This paper describes the development of a spatial analysis tool and modelling framework that together offer a solution for simulating landscape processes. A better approach to integrating landscape spatial analysis with Geographical Information Systems is advocated in this paper. Enhancements include special spatial operators and map algebra language constructs to handle dispersal and advective flows over landscape surfaces. These functional components to landscape modelling are developed in a modular way and are linked together in a modelling framework that performs dynamic simulation. The concepts and modelling framework are demonstrated using a hydrological modelling example. The approach provides a modelling environment for scientists and land resource managers to write and to visualize spatial process models with ease.

SimuMap: a computational system for spatial modelling

The paper presents a computational system based upon formal principles to run spatial models for environmental processes. The simulator is named SimuMap because it is typically used to simulate spatial processes over a mapped representation of terrain. A model is formally represented in SimuMap as a set of coupled sub-models. The paper considers the situation where spatial processes operate at different time levels, but are still integrated. An example of such a situation commonly occurs in watershed hydrology where overland flow and stream channel flow have very different flow rates but are highly related as they are subject to the same terrain runoff processes. SimuMap is able to run a network of sub-models that express different time-space derivatives for water flow processes. Sub-models may be coded generically with a map algebra programming language that uses a surface data model. To address the problem of differing time levels in simulation, the paper: (i) reviews general approaches for numerical solvers, (ii) considers the constraints that need to be enforced to use more adaptive time steps in discrete time specified simulations, and (iii) scaling transfer rates in equations that use different time bases for time-space derivatives. A multistep scheme is proposed for SimuMap. This is presented along with a description of its visual programming interface, its modelling formalisms and future plans. (C) 2003 Elsevier Ltd. All rights reserved.

Incorporating level set methods in Geographical Information Systems (GIS) for land-surface process modelling

Land-surface processes include a broad class of models that operate at a landscape scale. Current modelling approaches tend to be specialised towards one type of process, yet it is the interaction of processes that is increasing seen as important to obtain a more integrated approach to land management. This paper presents a technique and a tool that may be applied generically to landscape processes. The technique tracks moving interfaces across landscapes for processes such as water flow, biochemical diffusion, and plant dispersal. Its theoretical development applies a Lagrangian approach to motion over a Eulerian grid space by tracking quantities across a landscape as an evolving front. An algorithm for this technique, called level set method, is implemented in a geographical information system (GIS). It fits with a field data model in GIS and is implemented as operators in map algebra. The paper describes an implementation of the level set methods in a map algebra programming language, called MapScript, and gives example program scripts for applications in ecology and hydrology.

A relational algebra approach to ETL modeling

The MAP-i Doctoral Programme in Informatics, of the Universities of Minho, Aveiro and Porto

Local dimension and finite time prediction in coupled map lattices

Forecasting, for obvious reasons, often become the most important goal to be achieved. For spatially extended systems (e.g. atmospheric system) where the local nonlinearities lead to the most unpredictable chaotic evolution, it is highly desirable to have a simple diagnostic tool to identify regions of predictable behaviour. In this paper, we discuss the use of the bred vector (BV) dimension, a recently introduced statistics, to identify the regimes where a finite time forecast is feasible. Using the tools from dynamical systems theory and Bayesian modelling, we show the finite time predictability in two-dimensional coupled map lattices in the regions of low BV dimension. © Indian Academy of Sciences.

Differential algebra techniques for space applications

This project investigates the utility of differential algebra (DA) techniques applied to the problem of orbital dynamics with initial uncertainties in the orbital determination of the involved bodies. The use of DA theory allows the splitting of a common Monte Carlo simulation in two parts: the generation of a Taylor map of the final states with regard to the perturbation in the initial coordinates, and the evaluation of the map for many points. A propagator is implemented exploiting DA techniques, and tested in the field of asteroid impact risk monitoring with the potentially hazardous 2011 AG5 and 2007 VK184 as test cases. Results show that the new method is able to simulate 2.5 million trajectories with a precision good enough for the impact probability to be accurately reproduced, while running much faster than a traditional Monte Carlo approach (in 1 and 2 days, respectively).

The Hitchin map for one-nodal base curves of compact type

Studying moduli spaces of semistable Higgs bundles (E, \phi) of rank n on a smooth curve C, a key role is played by the spectral curve X (Hitchin), because an important result by Beauville-Narasimhan-Ramanan allows us to study isomorphism classes of such Higgs bundles in terms of isomorphism classes of rank-1 torsion-free sheaves on X. This way, the generic fibre of the Hitchin map, which associates to any semistable Higgs bundle the coefficients of the characteristic polynomial of \phi, is isomorphic to the Jacobian of X. Focusing on rank-2 Higgs data, this construction was extended by Barik to the case in which the curve C is reducible, one-nodal, having two smooth components. Such curve is called of compact type because its Picard group is compact. In this work, we describe and clarify the main points of the construction by Barik and we give examples, especially concerning generic fibres of the Hitchin map. Referring to Hausel-Pauly, we consider the case of SL(2,C)-Higgs bundles on a smooth base curve, which are such that the generic fibre of the Hitchin map is a subvariety of the Jacobian of X, the Prym variety. We recall the description of special loci, called endoscopic loci, such that the associated Prym variety is not connected. Then, letting G be an affine reductive group having underlying Lie algebra so(4,C), we consider G-Higgs bundles on a smooth base curve. Starting from the construction by Bradlow-Schaposnik, we discuss the associated endoscopic loci. By adapting these studies to a one-nodal base curve of compact type, we describe the fibre of the SL(2,C)-Hitchin map and of the G-Hitchin map, together with endoscopic loci. In the Appendix, we give an interpretation of generic spectral curves in terms of families of double covers.

Pseudodifferential operators with C*-algebra-valued symbols: Abstract characterizations

Given a separable unital C*-algebra C with norm parallel to center dot parallel to, let E-n denote the Banach-space completion of the C-valued Schwartz space on R-n with norm parallel to f parallel to(2)=parallel to < f, f >parallel to(1/2), < f, g >=integral f(x)* g(x)dx. The assignment of the pseudodifferential operator A=a(x,D) with C-valued symbol a(x,xi) to each smooth function with bounded derivatives a is an element of B-C(R-2n) defines an injective mapping O, from B-C(R-2n) to the set H of all operators with smooth orbit under the canonical action of the Heisenberg group on the algebra of all adjointable operators on the Hilbert module E-n. In this paper, we construct a left-inverse S for O and prove that S is injective if C is commutative. This generalizes Cordes' description of H in the scalar case. Combined with previous results of the second author, our main theorem implies that, given a skew-symmetric n x n matrix J and if C is commutative, then any A is an element of H which commutes with every pseudodifferential operator with symbol F(x+J xi), F is an element of B-C(R-n), is a pseudodifferential operator with symbol G(x - J xi), for some G is an element of B-C(R-n). That was conjectured by Rieffel.

Low dimensional behavior in three-dimensional coupled map lattices

The analysis of one-, two-, and three-dimensional coupled map lattices is here developed under a statistical and dynamical perspective. We show that the three-dimensional CML exhibits low dimensional behavior with long range correlation and the power spectrum follows 1/f noise. This approach leads to an integrated understanding of the most important properties of these universal models of spatiotemporal chaos. We perform a complete time series analysis of the model and investigate the dependence of the signal properties by change of dimension. (c) 2008 Elsevier Ltd. All rights reserved.

Tick saliva induces regulatory dendritic cells: MAP-kinases and Toll-like receptor-2 expression as potential targets

Ticks (Acari: Ixodidae) are bloodsucking ectoparasitic arthropods of human and veterinary medical importance. Tick saliva has been shown to contain a wide range of bioactive molecules with vasodilatory, antihemostatic, and immunomodulatory activities. We have previously demonstrated that saliva from Rhipicephalus sanguineus ticks inhibits the maturation of dendritic cells (DCs) stimulated with LPS. Here we examined the mechanism of this immune subversion, evaluating the effect of tick saliva on Toll-like receptor (TLR)-4 signalling pathway in bone marrow-derived DCs. We demonstrated that R. sanguineus tick saliva impairs maturation of DCs stimulated with LIPS, a TLR-4 ligand, leading to increased production of interleukin (IL)-10 and reduced synthesis of IL-12p70 and TNF-alpha. The immunomodulatory effect of the tick saliva on the production of pro-inflammatory cytokines by DCs stimulated with LPS was associated with the observation that tick saliva inhibits the activation of the ERK 1/2 and p38 MAP kinases. These effects were independent of the expression of TLR-4 on the surface of DCs. Additionally, saliva-treated DCs also presented a similar pattern of cytokine modulation in response to other TLR ligands. Since the recent literature reports that several parasites evade immune responses through TLR-2-mediated production of IL-10, we evaluated the effect of tick saliva on the percentage of TLR-2(+) DCs stimulated with the TLR-2 ligand lipoteicoic acid (LTA). The data showed that the population of DCs expressing TLR-2 was significantly increased in DCs treated with LTA plus saliva. In addition, tick saliva alone increased the expression of TLR-2 in a dose- and time-dependent manner. Our data suggest that tick saliva induces regulatory DCs, which secrete IL-10 and low levels of IL-12 and TNF-alpha when stimulated by TLR ligands. Such regulatory DCs are associated with expression of TLR-2 and inhibition of ERK and p38, which promotes the production of IL-10 and thus down-modulates the host`s immune response, possibly favouring susceptibility to tick infestations. (C) 2009 Elsevier B.V. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part I - A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops

An updated flow pattern map was developed for CO2 on the basis of the previous Cheng-Ribatski-Wojtan-Thome CO2 flow pattern map [1,2] to extend the flow pattern map to a wider range of conditions. A new annular flow to dryout transition (A-D) and a new dryout to mist flow transition (D-M) were proposed here. In addition, a bubbly flow region which generally occurs at high mass velocities and low vapor qualities was added to the updated flow pattern map. The updated flow pattern map is applicable to a much wider range of conditions: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to +25 degrees C (reduced pressures from 0.21 to 0.87). The updated flow pattern map was compared to independent experimental data of flow patterns for CO2 in the literature and it predicts the flow patterns well. Then, a database of CO2 two-phase flow pressure drop results from the literature was set up and the database was compared to the leading empirical pressure drop models: the correlations by Chisholm [3], Friedel [4], Gronnerud [5] and Muller-Steinhagen and Heck [6], a modified Chisholm correlation by Yoon et al. [7] and the flow pattern based model of Moreno Quiben and Thome [8-10]. None of these models was able to predict the CO2 pressure drop data well. Therefore, a new flow pattern based phenomenological model of two-phase flow frictional pressure drop for CO2 was developed by modifying the model of Moreno Quiben and Thome using the updated flow pattern map in this study and it predicts the CO2 pressure drop database quite well overall. (C) 2007 Elsevier Ltd. All rights reserved.

Spectral properties of chaotic signals generated by the skew tent map

Chaotic signals have been considered potentially attractive in many signal processing applications ranging from wideband communication systems to cryptography and watermarking. Besides, some devices as nonlinear adaptive filters and phase-locked loops can present chaotic behavior. In this paper, we derive analytical expressions for the autocorrelation sequence, power spectral density and essential bandwidth of chaotic signals generated by the skew tent map. From these results, we suggest possible applications in communication systems. (C) 2009 Elsevier B.V. All rights reserved.

Extension of the core map of common bean with EST-SSR, RGA, AFLP, and putative functional markers

Microsatellites and gene-derived markers are still underrepresented in the core molecular linkage map of common bean compared to other types of markers. In order to increase the density of the core map, a set of new markers were developed and mapped onto the RIL population derived from the `BAT93` x `Jalo EEP558` cross. The EST-SSR markers were first characterized using a set of 24 bean inbred lines. On average, the polymorphism information content was 0.40 and the mean number of alleles per locus was 2.7. In addition, AFLP and RGA markers based on the NBS-profiling method were developed and a subset of the mapped RGA was sequenced. With the integration of 282 new markers into the common bean core map, we were able to place markers with putative known function in some existing gaps including regions with QTL for resistance to anthracnose and rust. The distribution of the markers over 11 linkage groups is discussed and a newer version of the common bean core linkage map is proposed.

An integrated molecular map of yellow passion fruit based on simultaneous maximum-likelihood estimation of linkage and linkage phases

The development of genetic maps for auto-incompatible species, such as the yellow passion fruit (Passiflora edulis Sims f.flavicarpa Deg.) is restricted due to the unfeasibility of obtaining traditional mapping populations based on inbred lines. For this reason, yellow passion fruit linkage maps were generally constructed using a strategy known as two-way pseudo-testeross, based on monoparental dominant markers segregating in a 1:1 fashion. Due to the lack of information from these markers in one of the parents, two individual (parental) maps were obtained. However, integration of these maps is essential, and biparental markers can be used for such an operation. The objective of our study was to construct an integrated molecular map for a full-sib population of yellow passion fruit combining different loci configuration generated from amplified fragment length polymorphisms (AFLPs) and microsatellite markers and using a novel approach based on simultaneous maximum-likelihood estimation of linkage and linkage phases, specially designed for outcrossing species. Of the total number of loci, approximate to 76%, 21%, 0.7%, and 2.3% did segregate in 1:1, 3:1, 1:2:1, and 1:1:1:1 ratios, respectively. Ten linkage groups (LGs) were established with a logarithm of the odds (LOD) score >= 5.0 assuming a recombination fraction : <= 0.35. On average, 24 markers were assigned per LG, representing a total map length of 1687 cM, with a marker density of 6.9 cM. No markers were placed as accessories on the map as was done with previously constructed individual maps.

Development of a physical and genetic map of the virulent Wolbachia strain wMelPop

We report here the construction of a physical and genetic map of the virulent Wolbachia strain, wMelPop. This map was determined by ordering 28 chromosome fragments that resulted from digestion with the restriction endonucleases FseI, ApaI, SmaI, and AscI and were resolved by pulsed-field gel electrophoresis. Southern hybridization was done with 53 Wolbachia-specific genes as probes in order to determine the relative positions of these restriction fragments and use them to serve as markers. Comparison of the resulting map with the whole genome sequence of the closely related benign Wolbachia strain, wMel, shows that the two genomes are largely conserved in gene organization with the exception of a single inversion in the chromosome.