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.

Mycolactone-dependent depletion of endothelial cell thrombomodulin is strongly associated with fibrin deposition in buruli ulcer lesions

A well-known histopathological feature of diseased skin in Buruli ulcer (BU) is coagulative necrosis caused by the Mycobacterium ulcerans macrolide exotoxin mycolactone. Since the underlying mechanism is not known, we have investigated the effect of mycolactone on endothelial cells, focussing on the expression of surface anticoagulant molecules involved in the protein C anticoagulant pathway. Congenital deficiencies in this natural anticoagulant pathway are known to induce thrombotic complications such as purpura fulimans and spontaneous necrosis. Mycolactone profoundly decreased thrombomodulin (TM) expression on the surface of human dermal microvascular endothelial cells (HDMVEC) at doses as low as 2ng/ml and as early as 8hrs after exposure. TM activates protein C by altering thrombin’s substrate specificity, and exposure of HDMVEC to mycolactone for 24 hours resulted in an almost complete loss of the cells’ ability to produce activated protein C. Loss of TM was shown to be due to a previously described mechanism involving mycolactone-dependent blockade of Sec61 translocation that results in proteasome-dependent degradation of newly synthesised ER-transiting proteins. Indeed, depletion from cells determined by live-cell imaging of cells stably expressing a recombinant TM-GFP fusion protein occurred at the known turnover rate. In order to determine the relevance of these findings to BU disease, immunohistochemistry of punch biopsies from 40 BU lesions (31 ulcers, nine plaques) was performed. TM abundance was profoundly reduced in the subcutis of 78% of biopsies. Furthermore, it was confirmed that fibrin deposition is a common feature of BU lesions, particularly in the necrotic areas. These findings indicate that there is decreased ability to control thrombin generation in BU skin. Mycolactone’s effects on normal endothelial cell function, including its ability to activate the protein C anticoagulant pathway are strongly associated with this. Fibrin-driven tissue ischemia could contribute to the development of the tissue necrosis seen in BU lesions.

Fibrin activates GPVI in human and mouse platelets

The glycoprotein VI (GPVI)-Fc receptor γ (FcRγ) chain is the major platelet signaling receptor for collagen. Paradoxically, in a FeCl3 injury model, occlusion, but not initiation of thrombus formation, is delayed in GPVI-deficient and GPVI-depleted mice. In this study, we demonstrate that GPVI is a receptor for fibrin and speculate that this contributes to development of an occlusive thrombus. We observed a marked increase in tyrosine phosphorylation, including the FcRγ chain and Syk, in human and mouse platelets induced by thrombin in the presence of fibrinogen and the αIIbβ3 blocker eptifibatide. This was not seen in platelets stimulated by a protease activated receptor (PAR)-4 peptide, which is unable to generate fibrin from fibrinogen. The pattern of tyrosine phosphorylation was similar to that induced by activation of GPVI. Consistent with this, thrombin did not induce tyrosine phosphorylation of Syk and the FcRγ chain in GPVI-deficient mouse platelets. Mouse platelets underwent full spreading on fibrin but not fibrinogen, which was blocked in the presence of a Src kinase inhibitor or in the absence of GPVI. Spreading on fibrin was associated with phosphatidylserine exposure (procoagulant activity), and this too was blocked in GPVI-deficient platelets. The ectodomain of GPVI was shown to bind to immobilized monomeric and polymerized fibrin. A marked increase in embolization was seen following FeCl3 injury in GPVI-deficient mice, likely contributing to the delay in occlusion in this model. These results demonstrate that GPVI is a receptor for fibrin and provide evidence that this interaction contributes to thrombus growth and stability.

Calibration of uncertain flood inundation models using remotely sensed water levels

A traditional method of validating the performance of a flood model when remotely sensed data of the flood extent are available is to compare the predicted flood extent to that observed. The performance measure employed often uses areal pattern-matching to assess the degree to which the two extents overlap. Recently, remote sensing of flood extents using synthetic aperture radar (SAR) and airborne scanning laser altimetry (LIDAR) has made more straightforward the synoptic measurement of water surface elevations along flood waterlines, and this has emphasised the possibility of using alternative performance measures based on height. This paper considers the advantages that can accrue from using a performance measure based on waterline elevations rather than one based on areal patterns of wet and dry pixels. The two measures were compared for their ability to estimate flood inundation uncertainty maps from a set of model runs carried out to span the acceptable model parameter range in a GLUE-based analysis. A 1 in 5-year flood on the Thames in 1992 was used as a test event. As is typical for UK floods, only a single SAR image of observed flood extent was available for model calibration and validation. A simple implementation of a two-dimensional flood model (LISFLOOD-FP) was used to generate model flood extents for comparison with that observed. The performance measure based on height differences of corresponding points along the observed and modelled waterlines was found to be significantly more sensitive to the channel friction parameter than the measure based on areal patterns of flood extent. The former was able to restrict the parameter range of acceptable model runs and hence reduce the number of runs necessary to generate an inundation uncertainty map. A result of this was that there was less uncertainty in the final flood risk map. The uncertainty analysis included the effects of uncertainties in the observed flood extent as well as in model parameters. The height-based measure was found to be more sensitive when increased heighting accuracy was achieved by requiring that observed waterline heights varied slowly along the reach. The technique allows for the decomposition of the reach into sections, with different effective channel friction parameters used in different sections, which in this case resulted in lower r.m.s. height differences between observed and modelled waterlines than those achieved by runs using a single friction parameter for the whole reach. However, a validation of the modelled inundation uncertainty using the calibration event showed a significant difference between the uncertainty map and the observed flood extent. While this was true for both measures, the difference was especially significant for the height-based one. This is likely to be due to the conceptually simple flood inundation model and the coarse application resolution employed in this case. The increased sensitivity of the height-based measure may lead to an increased onus being placed on the model developer in the production of a valid model

Uncertainty assessment of a process-based integrated catchment model of phosphorus

Despite the many models developed for phosphorus concentration prediction at differing spatial and temporal scales, there has been little effort to quantify uncertainty in their predictions. Model prediction uncertainty quantification is desirable, for informed decision-making in river-systems management. An uncertainty analysis of the process-based model, integrated catchment model of phosphorus (INCA-P), within the generalised likelihood uncertainty estimation (GLUE) framework is presented. The framework is applied to the Lugg catchment (1,077 km2), a River Wye tributary, on the England–Wales border. Daily discharge and monthly phosphorus (total reactive and total), for a limited number of reaches, are used to initially assess uncertainty and sensitivity of 44 model parameters, identified as being most important for discharge and phosphorus predictions. This study demonstrates that parameter homogeneity assumptions (spatial heterogeneity is treated as land use type fractional areas) can achieve higher model fits, than a previous expertly calibrated parameter set. The model is capable of reproducing the hydrology, but a threshold Nash-Sutcliffe co-efficient of determination (E or R 2) of 0.3 is not achieved when simulating observed total phosphorus (TP) data in the upland reaches or total reactive phosphorus (TRP) in any reach. Despite this, the model reproduces the general dynamics of TP and TRP, in point source dominated lower reaches. This paper discusses why this application of INCA-P fails to find any parameter sets, which simultaneously describe all observed data acceptably. The discussion focuses on uncertainty of readily available input data, and whether such process-based models should be used when there isn’t sufficient data to support the many parameters.

Driving saccade to pursuit using image motion

Within the context of active vision, scant attention has been paid to the execution of motion saccades—rapid re-adjustments of the direction of gaze to attend to moving objects. In this paper we first develop a methodology for, and give real-time demonstrations of, the use of motion detection and segmentation processes to initiate capture saccades towards a moving object. The saccade is driven by both position and velocity of the moving target under the assumption of constant target velocity, using prediction to overcome the delay introduced by visual processing. We next demonstrate the use of a first order approximation to the segmented motion field to compute bounds on the time-to-contact in the presence of looming motion. If the bound falls below a safe limit, a panic saccade is fired, moving the camera away from the approaching object. We then describe the use of image motion to realize smooth pursuit, tracking using velocity information alone, where the camera is moved so as to null a single constant image motion fitted within a central image region. Finally, we glue together capture saccades with smooth pursuit, thus effecting changes in both what is being attended to and how it is being attended to. To couple the different visual activities of waiting, saccading, pursuing and panicking, we use a finite state machine which provides inherent robustness outside of visual processing and provides a means of making repeated exploration. We demonstrate in repeated trials that the transition from saccadic motion to tracking is more likely to succeed using position and velocity control, than when using position alone.

Does increased hydrochemical model complexity decrease robustness?

The aim of this study was, within a sensitivity analysis framework, to determine if additional model complexity gives a better capability to model the hydrology and nitrogen dynamics of a small Mediterranean forested catchment or if the additional parameters cause over-fitting. Three nitrogen-models of varying hydrological complexity were considered. For each model, general sensitivity analysis (GSA) and Generalized Likelihood Uncertainty Estimation (GLUE) were applied, each based on 100,000 Monte Carlo simulations. The results highlighted the most complex structure as the most appropriate, providing the best representation of the non-linear patterns observed in the flow and streamwater nitrate concentrations between 1999 and 2002. Its 5% and 95% GLUE bounds, obtained considering a multi-objective approach, provide the narrowest band for streamwater nitrogen, which suggests increased model robustness, though all models exhibit periods of inconsistent good and poor fits between simulated outcomes and observed data. The results confirm the importance of the riparian zone in controlling the short-term (daily) streamwater nitrogen dynamics in this catchment but not the overall flux of nitrogen from the catchment. It was also shown that as the complexity of a hydrological model increases over-parameterisation occurs, but the converse is true for a water quality model where additional process representation leads to additional acceptable model simulations. Water quality data help constrain the hydrological representation in process-based models. Increased complexity was justifiable for modelling river-system hydrochemistry. Increased complexity was justifiable for modelling river-system hydrochemistry.

Genetic relatedness and phenotypic characteristics of Treponema associated with human periodontal tissues and ruminant foot disease

Treponema have been implicated recently in the pathogenesis of digital dermatitis (DID) and contagious ovine digital dermatitis (CODD) that are infectious diseases of bovine and ovine foot tissues, respectively. Previous analyses of treponemal 16S rDNA sequences, PCR-amplified directly from DID or CODD lesions, have suggested relatedness of animal Treponema to some human oral Treponema species isolated from periodontal tissues. In this study a range of adhesion and virulence-related properties of three animal Treponema isolates have been compared with representative human oral strains of Treponema denticola and Treponema vincentii. In adhesion assays using biotinylated treponemal cells, T denticola cells bound in consistently higher numbers to fibronectin, laminin, collagen type 1, gelatin, keratin and lactoferrin than did T. vincentii or animal Treponema isolates. However, animal DID strains adhered to fibrinogen at equivalent or greater levels than T denticola. All Treponema strains bound to the amino-terminal heparin l/fibrin I domain of fibronectin. 16S rDNA sequence analyses placed ovine strain UB1090 and bovine strain UB1467 within a cluster that was phylogenetically related to T vincentii, while ovine strain UB1466 appeared more closely related to T denticola. These observations correlated with phenotypic properties. Thus, T denticola ATCC 35405, GM-1, and Treponema UB1466 had similar outer-membrane protein profiles, produced chymotrypsin-like protease (CTLP), trypsin-like protease and high levels of proline iminopeptidase, and co-aggregated with human oral bacteria Porphyromonas gingivalis and Streptococcus crista. Conversely, T vincentii ATCC 35580, D2A-2, and animal strains UB1090 and UB1467 did not express CTLP or trypsin-like protease and did not co-aggregate with P. gingivalis or S. crista. Taken collectively, these results suggest that human oral-related Treponema have broad host specificity and that similar control or preventive strategies might be developed for human and animal Treponema-associated infections.

Binding properties and adhesion-mediating regions of the major sheath protein of Treponema denticola ATCC 35405

There is growing evidence that a number of oral Treponema species, in particular Treponema denticola, are associated with the progression of human periodontal disease. The major sheath (or surface) protein (Msp) of T. denticola is implicated in adhesion of bacteria to host cells and tissue proteins and is likely to be an important virulence factor. However, the binding regions of the Msp are not known. We have purified from Escherichia coli recombinant Msp (rMsp) polypeptides corresponding to the following: full-length Msp (rMsp) minus 13 N-terminal amino acid (aa) residues, an amino-terminal fragment (rN-Msp, 189 aa residues), a 57-aa residue segment from the central region (rV-Msp), and a C-terminal fragment (rC-Msp, 272 aa residues). rMsp (530 aa residues) bound to immobilized fibronectin, keratin, laminin, collagen type 1, fibrinogen, hyaluronic acid, and heparin. The N- and V-region polypeptides, but not rC-Msp, also bound to these substrates. Binding of rMsp to fibronectin was targeted to the N-terminal heparin I/fibrin I domain. Antibodies to the N-region or V-region polypeptides, but not antibodies to the rC-Msp fragment, blocked adhesion of T. denticola ATCC 35405 cells to a range of host protein molecules. These results suggest that the N-terminal half of Msp carries epitopes that are surface exposed and that are involved in mediating adhesion. Binding of rMsp onto the cell surface of low-level fibronectin-binding Treponema isolates conferred a 10-fold increase in fibronectin binding. This confirms that Msp functions autonomously as an adhesin and raises the possibility that phenotypic complementation of virulence functions might occur within mixed populations of Treponema species.

Platelet protein disulfide isomerase is required for thrombus formation but not essential for hemostasis in mice

Protein disulfide isomerase (PDI) derived from intravascular cells is required for thrombus formation. However, it remains unclear whether platelet PDI contributes to the process. Using platelet-specific PDI-deficient mice, we demonstrate that PDI-null platelets have defects in aggregation and ATP secretion induced by thrombin, collagen, and ADP. Such defects were rescued by exogenously-added wild-type but not mutant PDI, indicating that the isomerase activity of platelet surface PDI is critical for the regulatory effect. PDI-deficient platelets expressed increased levels of intracellular ERp57 and ERp72. Platelet PDI regulated αIIbβ3 integrin activation but not P-selectin exposure, Ca2+ mobilization, β3-talin interaction, and platelet spreading on immobilized fibrinogen. Inhibition of ERp57 further diminished αIIbβ3 integrin activation, aggregation and ATP secretion of activated PDI-deficient platelets, suggesting distinct roles of PDI and ERp57 in platelet functions. We found that platelet PDI is important for thrombus formation on collagen-coated surfaces under arteriolar shear. Intravital microscopy demonstrates that platelet PDI is important for platelet accumulation but not initial adhesion and fibrin generation following laser-induced arteriolar injury. Tail bleeding time and blood loss in platelet-specific PDI-deficient mice were not significantly increased. Our results provide important evidence that platelet PDI is essential for thrombus formation but not for hemostasis in mice.

Global runoff routing with the hydrological component of the ECMWF NWP system

A global river routing scheme coupled to the ECMWF land surface model is implemented and tested within the framework of the Global Soil Wetness Project II, to evaluate the feasibility of modelling global river runoff at a daily time scale. The exercise is designed to provide benchmark river runoff predictions needed to verify the land surface model. Ten years of daily runoff produced by the HTESSEL land surface scheme is input into the TRIP2 river routing scheme in order to generate daily river runoff. These are then compared to river runoff observations from the Global Runoff Data Centre (GRDC) in order to evaluate the potential and the limitations. A notable source of inaccuracy is bias between observed and modelled discharges which is not primarily due to the modelling system but instead of to the forcing and quality of observations and seems uncorrelated to the river catchment size. A global sensitivity analysis and Generalised Likelihood Uncertainty Estimation (GLUE) uncertainty analysis are applied to the global routing model. The ground water delay parameter is identified as being the most sensitive calibration parameter. Significant uncertainties are found in results, and those due to parameterisation of the routing model are quantified. The difficulty involved in parameterising global river discharge models is discussed. Detailed river runoff simulations are shown for the river Danube, which match well observed river runoff in upstream river transects. Results show that although there are errors in runoff predictions, model results are encouraging and certainly indicative of useful runoff predictions, particularly for the purpose of verifying the land surface scheme hydrologicly. Potential of this modelling system on future applications such as river runoff forecasting and climate impact studies is highlighted. Copyright © 2009 Royal Meteorological Society.

Calibration of structure in a distributed forecasting model for a semiarid flash flood: dynamic surface storage and channel roughness

Flash floods pose a significant danger for life and property. Unfortunately, in arid and semiarid environment the runoff generation shows a complex non-linear behavior with a strong spatial and temporal non-uniformity. As a result, the predictions made by physically-based simulations in semiarid areas are subject to great uncertainty, and a failure in the predictive behavior of existing models is common. Thus better descriptions of physical processes at the watershed scale need to be incorporated into the hydrological model structures. For example, terrain relief has been systematically considered static in flood modelling at the watershed scale. Here, we show that the integrated effect of small distributed relief variations originated through concurrent hydrological processes within a storm event was significant on the watershed scale hydrograph. We model these observations by introducing dynamic formulations of two relief-related parameters at diverse scales: maximum depression storage, and roughness coefficient in channels. In the final (a posteriori) model structure these parameters are allowed to be both time-constant or time-varying. The case under study is a convective storm in a semiarid Mediterranean watershed with ephemeral channels and high agricultural pressures (the Rambla del Albujón watershed; 556 km 2 ), which showed a complex multi-peak response. First, to obtain quasi-sensible simulations in the (a priori) model with time-constant relief-related parameters, a spatially distributed parameterization was strictly required. Second, a generalized likelihood uncertainty estimation (GLUE) inference applied to the improved model structure, and conditioned to observed nested hydrographs, showed that accounting for dynamic relief-related parameters led to improved simulations. The discussion is finally broadened by considering the use of the calibrated model both to analyze the sensitivity of the watershed to storm motion and to attempt the flood forecasting of a stratiform event with highly different behavior.

Tityus discrepans scorpion venom activates platelets through GPVI and a novel Src-dependent signaling pathway

In humans and other mammals, Tityus discrepans (Td) scorpion envenomation produces a variety of systemic effects including respiratory distress, a generalized inflammatory reaction, modulation of blood pressure, fibrin formation, and platelet activation. For many of these effects, the venom components and underlying mechanisms are not known. In the present study, we demonstrate that Td venom (TdV) stimulates integrin αIIbβ3-dependent aggregation of washed human and mouse platelets downstream of Src kinase activation. The pattern of increase in tyrosine phosphorylation induced by TdV in human platelets is similar to that induced by the collagen receptor GPVI, and includes FcR γ-chain, Syk, and PLC γ 2. Confirmation of GPVI activation by TdV was achieved by expression of human GPVI in chicken DT40 B cells and use of a reporter assay. To our surprise, TdV was able to activate mouse platelets deficient in the GPVI-FcR γ-chain complex through a pathway that was also dependent on Src kinases. TdV therefore activates platelets through GPVI and a second, as yet unidentified Src kinase-dependent pathway.