An application of the GLUE methodology for estimating the parameters of the INCA-N model

The conceptual and parameter uncertainty of the semi-distributed INCA-N (Integrated Nutrients in Catchments-Nitrogen) model was studied using the GLUE (Generalized Likelihood Uncertainty Estimation) methodology combined with quantitative experimental knowledge, the concept known as 'soft data'. Cumulative inorganic N leaching, annual plant N uptake and annual mineralization proved to be useful soft data to constrain the parameter space. The INCA-N model was able to simulate the seasonal and inter-annual variations in the stream-water nitrate concentrations, although the lowest concentrations during the growing season were not reproduced. This suggested that there were some retention processes or losses either in peatland/wetland areas or in the river which were not included in the INCA-N model. The results of the study suggested that soft data was a way to reduce parameter equifinality, and that the calibration and testing of distributed hydrological and nutrient leaching models should be based both on runoff and/or nutrient concentration data and the qualitative knowledge of experimentalist. (c) 2006 Elsevier B.V. All rights reserved.

Impacts of forestry on nitrogen in upland and lowland catchments: a comparison of the River Severn at Plynlimon in mid-Wales and the Bedford Ouse in south-east England using the INCA Model

The impacts of afforestation at Plynlimon in the Severn catchment, mid-Wales. and in the Bedford Ouse catchment in south-east England are evaluated using the INCA model to simulate Nitrogen (N) fluxes and concentrations. The INCA model represents the key hydrological and N processes operating in catchments and simulates the daily dynamic behaviour as well as the annual fluxes. INCA has been applied to five years of data front the Hafren and Hore headwater sub-catchments (6.8 km(2) area in total) of the River Severn at Plytilimon and the model was calibrated and validated against field data. Simulation of afforestation is achieved by altering the uptake rate parameters in the model. INCA simulates the daily N behaviour in the catchments with good accuracy as well as reconstructing the annual budgets for N release following clearfelling a four-fold increase in N fluxes was followed by a slow recovery after re-afforestation. For comparison, INCA has been applied to the large (8380 km(2)) Bedford Ouse catchment to investigate the impact of replacing 20% arable land with forestry. The reduction in fertiliser inputs from arable farming and the N uptake by the forest are predicted to reduce the N flux reaching the main river system, leading to a 33% reduction in N-Nitrate concentrations in the river water.

Impacts of climate change on in-stream nitrogen in a lowland chalk stream: An appraisal of adaptation strategies

The impacts of climate change on nitrogen (N) in a lowland chalk stream are investigated using a dynamic modelling approach. The INCA-N model is used to simulate transient daily hydrology and water quality in the River Kennet using temperature and precipitation scenarios downscaled from the General Circulation Model (GCM) output for the period 1961-2100. The three GCMs (CGCM2, CSIRO and HadCM3) yield very different river flow regimes with the latter projecting significant periods of drought in the second half of the 21st century. Stream-water N concentrations increase over time as higher temperatures enhance N release from the soil, and lower river flows reduce the dilution capacity of the river. Particular problems are shown to occur following severe droughts when N mineralization is high and the subsequent breaking of the drought releases high nitrate loads into the river system. Possible strategies for reducing climate-driven N loads are explored using INCA-N. The measures include land use change or fertiliser reduction, reduction in atmospheric nitrate and ammonium deposition, and the introduction of water meadows or connected wetlands adjacent to the river. The most effective strategy is to change land use or reduce fertiliser use, followed by water meadow creation, and atmospheric pollution controls. Finally, a combined approach involving all three strategies is investigated and shown to reduce in-stream nitrate concentrations to those pre-1950s even under climate change. (c) 2006 Elsevier B.V. All rights reserved.

LiDAR mapping of tidal marshes for ecogeomorphological modelling in the TIDE project

The European research project TIDE (Tidal Inlets Dynamics and Environment) is developing and validating coupled models describing the morphological, biological and ecological evolution of tidal environments. The interactions between the physical and biological processes occurring in these regions requires that the system be studied as a whole rather than as separate parts. Extensive use of remote sensing including LiDAR is being made to provide validation data for the modelling. This paper describes the different uses of LiDAR within the project and their relevance to the TIDE science objectives. LiDAR data have been acquired from three different environments, the Venice Lagoon in Italy, Morecambe Bay in England, and the Eden estuary in Scotland. LiDAR accuracy at each site has been evaluated using ground reference data acquired with differential GPS. A semi-automatic technique has been developed to extract tidal channel networks from LiDAR data either used alone or fused with aerial photography. While the resulting networks may require some correction, the procedure does allow network extraction over large areas using objective criteria and reduces fieldwork requirements. The networks extracted may subsequently be used in geomorphological analyses, for example to describe the drainage patterns induced by networks and to examine the rate of change of networks. Estimation of the heights of the low and sparse vegetation on marshes is being investigated by analysis of the statistical distribution of the measured LiDAR heights. Species having different mean heights may be separated using the first-order moments of the height distribution.

Forward modelling to determine the observational signatures of white-light imaging and interplanetary scintillation for the propagation of an interplanetary shock in the ecliptic plane

Recent coordinated observations of interplanetary scintillation (IPS) from the EISCAT, MERLIN, and STELab, and stereoscopic white-light imaging from the two heliospheric imagers (HIs) onboard the twin STEREO spacecraft are significant to continuously track the propagation and evolution of solar eruptions throughout interplanetary space. In order to obtain a better understanding of the observational signatures in these two remote-sensing techniques, the magnetohydrodynamics of the macro-scale interplanetary disturbance and the radio-wave scattering of the micro-scale electron-density fluctuation are coupled and investigated using a newly constructed multi-scale numerical model. This model is then applied to a case of an interplanetary shock propagation within the ecliptic plane. The shock could be nearly invisible to an HI, once entering the Thomson-scattering sphere of the HI. The asymmetry in the optical images between the western and eastern HIs suggests the shock propagation off the Sun–Earth line. Meanwhile, an IPS signal, strongly dependent on the local electron density, is insensitive to the density cavity far downstream of the shock front. When this cavity (or the shock nose) is cut through by an IPS ray-path, a single speed component at the flank (or the nose) of the shock can be recorded; when an IPS ray-path penetrates the sheath between the shock nose and this cavity, two speed components at the sheath and flank can be detected. Moreover, once a shock front touches an IPS ray-path, the derived position and speed at the irregularity source of this IPS signal, together with an assumption of a radial and constant propagation of the shock, can be used to estimate the later appearance of the shock front in the elongation of the HI field of view. The results of synthetic measurements from forward modelling are helpful in inferring the in-situ properties of coronal mass ejection from real observational data via an inverse approach.

Mathematical modelling of decision making processes in construction projects

Many different individuals, who have their own expertise and criteria for decision making, are involved in making decisions on construction projects. Decision-making processes are thus significantly affected by communication, in which a dynamic performance of human intentions leads to unpredictable outcomes. In order to theorise the decision making processes including communication, it is argued here that the decision making processes resemble evolutionary dynamics in terms of both selection and mutation, which can be expressed by the replicator-mutator equation. To support this argument, a mathematical model of decision making has been made from an analogy with evolutionary dynamics, in which there are three variables: initial support rate, business hierarchy, and power of persuasion. On the other hand, a survey of patterns in decision making in construction projects has also been performed through self-administered mail questionnaire to construction practitioners. Consequently, comparison between the numerical analysis of mathematical model and the statistical analysis of empirical data has shown a significant potential of the replicator-mutator equation as a tool to study dynamic properties of intentions in communication.

A revised Tobit procedure for mitigating bias in the presence of non-zero censoring with an application to milk-market participation in the Ethiopian highlands

Fixed transactions costs that prohibit exchange engender bias in supply analysis due to censoring of the sample observations. The associated bias in conventional regression procedures applied to censored data and the construction of robust methods for mitigating bias have been preoccupations of applied economists since Tobin [Econometrica 26 (1958) 24]. This literature assumes that the true point of censoring in the data is zero and, when this is not the case, imparts a bias to parameter estimates of the censored regression model. We conjecture that this bias can be significant; affirm this from experiments; and suggest techniques for mitigating this bias using Bayesian procedures. The bias-mitigating procedures are based on modifications of the key step that facilitates Bayesian estimation of the censored regression model; are easy to implement; work well in both small and large samples; and lead to significantly improved inference in the censored regression model. These findings are important in light of the widespread use of the zero-censored Tobit regression and we investigate their consequences using data on milk-market participation in the Ethiopian highlands. (C) 2004 Elsevier B.V. All rights reserved.

Time patterns in UK demand for alcohol and tobacco: an application of the EM algorithm

Capturing the pattern of structural change is a relevant task in applied demand analysis, as consumer preferences may vary significantly over time. Filtering and smoothing techniques have recently played an increasingly relevant role. A dynamic Almost Ideal Demand System with random walk parameters is estimated in order to detect modifications in consumer habits and preferences, as well as changes in the behavioural response to prices and income. Systemwise estimation, consistent with the underlying constraints from economic theory, is achieved through the EM algorithm. The proposed model is applied to UK aggregate consumption of alcohol and tobacco, using quarterly data from 1963 to 2003. Increased alcohol consumption is explained by a preference shift, addictive behaviour and a lower price elasticity. The dynamic and time-varying specification is consistent with the theoretical requirements imposed at each sample point. (c) 2005 Elsevier B.V. All rights reserved.

Application of real-time and multiplex PCR assays to study leaf blotch epidemics in barley

Leaf blotch, caused by Rhynchosporium secalis, was studied in a range of winter barley cultivars using a combination of traditional plant pathological techniques and newly developed multiplex and real-time polymerase chain reaction (PCR) assays. Using PCR, symptomless leaf blotch colonization was shown to occur throughout the growing season in the resistant winter barley cv. Leonie. The dynamics of colonization throughout the growing season were similar in both Leonie and Vertige, a susceptible cultivar. However, pathogen DNA levels were approximately 10-fold higher in the susceptible cultivar, which expressed symptoms throughout the growing season. Visual assessments and PCR also were used to determine levels of R. secalis colonization and infection in samples from a field experiment used to test a range of winter barley cultivars with different levels of leaf blotch resistance. The correlation between the PCR and visual assessment data was better at higher infection levels (R(2) = 0.81 for leaf samples with >0.3% disease). Although resistance ratings did not correlate well with levels of disease for all cultivars tested, low levels of infection were observed in the cultivar with the highest resistance rating and high levels of infection in the cultivar with the lowest resistance rating.

Modelling and validation of Lactobacillus plantarum fermentations in cereal-based media with different sugar concentrations and buffering capacities

An unstructured mathematical model is proposed to describe the fermentation kinetics of growth, lactic acid production, pH and sugar consumption by Lactobacillus plantarum as a function of the buffering capacity and initial glucose concentration of the culture media. Initially the experimental data of L plantarum fermentations in synthetic media with different buffering capacity and glucose were fitted to a set of primary models. Later the parameters obtained from these models were used to establish mathematical relationships with the independent variables tested. The models were validated with 6 fermentations of L. plantarum in different cereal-based media. In most cases the proposed models adequately describe the biochemical changes taking place during fermentation and are a promising approach for the formulation of cereal-based probiotic foods. (C) 2008 Elsevier B.V. All rights reserved.

Application of Gaussian multi-scale representation to feature tracking in meteorological satellite imagery

Gaussian multi-scale representation is a mathematical framework that allows to analyse images at different scales in a consistent manner, and to handle derivatives in a way deeply connected to scale. This paper uses Gaussian multi-scale representation to investigate several aspects of the derivation of atmospheric motion vectors (AMVs) from water vapour imagery. The contribution of different spatial frequencies to the tracking is studied, for a range of tracer sizes, and a number of tracer selection methods are presented and compared, using WV 6.2 images from the geostationary satellite MSG-2.

Algorithmic statistical analysis of electrophysiological data for the investigation of structure-activity relationship in single neurons

We are developing computational tools supporting the detailed analysis of the dependence of neural electrophysiological response on dendritic morphology. We approach this problem by combining simulations of faithful models of neurons (experimental real life morphological data with known models of channel kinetics) with algorithmic extraction of morphological and physiological parameters and statistical analysis. In this paper, we present the novel method for an automatic recognition of spike trains in voltage traces, which eliminates the need for human intervention. This enables classification of waveforms with consistent criteria across all the analyzed traces and so it amounts to reduction of the noise in the data. This method allows for an automatic extraction of relevant physiological parameters necessary for further statistical analysis. In order to illustrate the usefulness of this procedure to analyze voltage traces, we characterized the influence of the somatic current injection level on several electrophysiological parameters in a set of modeled neurons. This application suggests that such an algorithmic processing of physiological data extracts parameters in a suitable form for further investigation of structure-activity relationship in single neurons.

Impacts of pollution and climate change on ombrotrophic Sphagnum species in the UK: analysis of uncertainties in two empirical niche models

A significant challenge in the prediction of climate change impacts on ecosystems and biodiversity is quantifying the sources of uncertainty that emerge within and between different models. Statistical species niche models have grown in popularity, yet no single best technique has been identified reflecting differing performance in different situations. Our aim was to quantify uncertainties associated with the application of 2 complimentary modelling techniques. Generalised linear mixed models (GLMM) and generalised additive mixed models (GAMM) were used to model the realised niche of ombrotrophic Sphagnum species in British peatlands. These models were then used to predict changes in Sphagnum cover between 2020 and 2050 based on projections of climate change and atmospheric deposition of nitrogen and sulphur. Over 90% of the variation in the GLMM predictions was due to niche model parameter uncertainty, dropping to 14% for the GAMM. After having covaried out other factors, average variation in predicted values of Sphagnum cover across UK peatlands was the next largest source of variation (8% for the GLMM and 86% for the GAMM). The better performance of the GAMM needs to be weighed against its tendency to overfit the training data. While our niche models are only a first approximation, we used them to undertake a preliminary evaluation of the relative importance of climate change and nitrogen and sulphur deposition and the geographic locations of the largest expected changes in Sphagnum cover. Predicted changes in cover were all small (generally <1% in an average 4 m2 unit area) but also highly uncertain. Peatlands expected to be most affected by climate change in combination with atmospheric pollution were Dartmoor, Brecon Beacons and the western Lake District.

Application of interpolymer complexation to coat solid and soft surfaces

Water-soluble polymers are often capable of forming interpolymer complexes in solutions and at interfaces, which offers an excellent opportunity for surface modification. The complex formation may be driven by H-bonding between poly(carboxylic acids) and non-ionic polymers or by electrostatic attraction between oppositely-charged polyelectrolytes. In the present communication the following applications of interpolymer complexation in coating technologies will be considered: (1) Complexation between poly(acrylic acid) and non-ionic polymers via H-bonding was used to coat glass surfaces. It was realised using layer-by-layer deposition of IPC on glass surfaces with subsequent cross-linking of dry multilayers by thermal treatment. Depending on the glass surface functionality this complexation resulted in detachable and non-detachable hydrogel films; (2) Electrostatic layer-by-layer self-assembly between glycol chitosan and bovine serum albumin (BSA) was used to coat magnetic nanoparticles. It was demonstrated that the native structure of BSA remains unaffected by the self-assembling process.