Modelling uncertainty of flood quantile estimation at ungauged sites by Bayesian networks

Prediction at ungauged sites is essential for water resources planning and management. Ungauged sites have no observations about the magnitude of floods, but some site and basin characteristics are known. Regression models relate physiographic and climatic basin characteristics to flood quantiles, which can be estimated from observed data at gauged sites. However, these models assume linear relationships between variables Prediction intervals are estimated by the variance of the residuals in the estimated model. Furthermore, the effect of the uncertainties in the explanatory variables on the dependent variable cannot be assessed. This paper presents a methodology to propagate the uncertainties that arise in the process of predicting flood quantiles at ungauged basins by a regression model. In addition, Bayesian networks were explored as a feasible tool for predicting flood quantiles at ungauged sites. Bayesian networks benefit from taking into account uncertainties thanks to their probabilistic nature. They are able to capture non-linear relationships between variables and they give a probability distribution of discharges as result. The methodology was applied to a case study in the Tagus basin in Spain.

Statistical distributions obtained from the compression of monodisperse, soft and frictionless particles

The cyclic compression of several granular systems has been simulated with a molecular dynamics code. All the samples consisted of bidimensional, soft, frictionless and equal-sized particles that were initially arranged according to a squared lattice and were compressed by randomly generated irregular walls. The compression protocols can be described by some control variables (volume or external force acting on the walls) and by some dimensionless factors, that relate stiffness, density, diameter, damping ratio and water surface tension to the external forces, displacements and periods. Each protocol, that is associated to a dynamic process, results in an arrangement with its own macroscopic features: volume (or packing ratio), coordination number, and stress; and the differences between packings can be highly significant. The statistical distribution of the force-moment state of the particles (i.e. the equivalent average stress multiplied by the volume) is analyzed. In spite of the lack of a theoretical framework based on statistical mechanics specific for these protocols, it is shown how the obtained distributions of mean and relative deviatoric force-moment are. Then it is discussed on the nature of these distributions and on their relation to specific protocols.

Univariate and bivariate truncated von Mises distributions

In this article we study the univariate and bivariate truncated von Mises distribution, as a generalization of the von Mises distribution (\cite{jupp1989}), (\cite{mardia2000directional}). This implies the addition of two or four new truncation parameters in the univariate and, bivariate cases, respectively. The results include the definition, properties of the distribution and maximum likelihood estimators for the univariate and bivariate cases. Additionally, the analysis of the bivariate case shows how the conditional distribution is a truncated von Mises distribution, whereas the marginal distribution that generalizes the distribution introduced in \cite{repe}. From the viewpoint of applications, we test the distribution with simulated data, as well as with data regarding leaf inclination angles (\cite{safari}) and dihedral angles in protein chains (\cite{prote}). This research aims to assert this probability distribution as a potential option for modelling or simulating any kind of phenomena where circular distributions are applicable.\par

Computer simulation of the interplay between fractal structures and surrounding heterogeneous multifractal distributions. Applications

In a large number of physical, biological and environmental processes interfaces with high irregular geometry appear separating media (phases) in which the heterogeneity of constituents is present. In this work the quantification of the interplay between irregular structures and surrounding heterogeneous distributions in the plane is made For a geometric set image and a mass distribution (measure) image supported in image, being image, the mass image gives account of the interplay between the geometric structure and the surrounding distribution. A computation method is developed for the estimation and corresponding scaling analysis of image, being image a fractal plane set of Minkowski dimension image and image a multifractal measure produced by random multiplicative cascades. The method is applied to natural and mathematical fractal structures in order to study the influence of both, the irregularity of the geometric structure and the heterogeneity of the distribution, in the scaling of image. Applications to the analysis and modeling of interplay of phases in environmental scenarios are given.

Computer simulation of random parckings for self-similar particle size distributions in soil and granular materials: porosity and pore size distribution

A 2D computer simulation method of random packings is applied to sets of particles generated by a self-similar uniparametric model for particle size distributions (PSDs) in granular media. The parameter p which controls the model is the proportion of mass of particles corresponding to the left half of the normalized size interval [0,1]. First the influence on the total porosity of the parameter p is analyzed and interpreted. It is shown that such parameter, and the fractal exponent of the associated power scaling, are efficient packing parameters, but this last one is not in the way predicted in a former published work addressing an analogous research in artificial granular materials. The total porosity reaches the minimum value for p = 0.6. Limited information on the pore size distribution is obtained from the packing simulations and by means of morphological analysis methods. Results show that the range of pore sizes increases for decreasing values of p showing also different shape in the volume pore size distribution. Further research including simulations with a greater number of particles and image resolution are required to obtain finer results on the hierarchical structure of pore space.

Computer Simulation of Packing of Particles with Size Distributions Produced by Fragmentation Processes

Fragmentation schemes inspired by theoretical results and conjectures of Kolmogorov are applied to produce particle size distributions of different natures, depending on fragmentation parameters. A two-dimensional computer simulation method of packing is applied to the resulting distributions and the void fraction is evaluated. The relationship between the void fraction and characteristic parameters of the fragmentation process is studied.

Computer simulation of packing of particles with size distributions produced by fragmentation processes

Fragmentation schemes inspired by theoretical results and conjectures of Kolmogorov are applied to produce particle size distributions of different natures, depending on fragmentation parameters. A two-dimensional computer simulation method of packing is applied to the resulting distributions and the void fraction is evaluated. The relationship between the void fraction and characteristic parameters of the fragmentation process is studied.

Intermittent pluri-sink model and the emergence of complex heterogeneity patterns: a simple paradigm for explaining complexity in soil chemical distributions

The spatial complexity of the distribution of organic matter, chemicals, nutrients, pollutants has been demonstrated to have multifractal nature (Kravchenco et al. [1]). This fact supports the possibility of existence of some emergent heterogeneity structure built under the evolution of the system. The aim of this note is providing a consistent explanation to the mentioned results via an extremely simple model.

Mapping species distributions : A comparison of skilled naturalist and lay citizen science recording

Acknowledgements We are grateful to Elaine O’Mahony, Imogen Pearce, Richard Comont, Anthony McCluskey and other BBCT staff for the many hours of BeeWatch species identification and for all people who submitted sightings to BeeWatch, OPAL, BWARS and the various local recording schemes and societies. We thank the NBN for allowing us to download the bumblebee records without strings attached, and the Essex, Greater London, Cumbria and Sussex based recording centres for providing records upon request. Finally, we are indebted to Tom August and two anonymous reviewers for their valuable critique on an earlier version of this work.

Equilibrium distributions of topological states in circular DNA: Interplay of supercoiling and knotting

Two variables define the topological state of closed double-stranded DNA: the knot type, K, and ΔLk, the linking number difference from relaxed DNA. The equilibrium distribution of probabilities of these states, P(ΔLk, K), is related to two conditional distributions: P(ΔLk|K), the distribution of ΔLk for a particular K, and P(K|ΔLk) and also to two simple distributions: P(ΔLk), the distribution of ΔLk irrespective of K, and P(K). We explored the relationships between these distributions. P(ΔLk, K), P(ΔLk), and P(K|ΔLk) were calculated from the simulated distributions of P(ΔLk|K) and of P(K). The calculated distributions agreed with previous experimental and theoretical results and greatly advanced on them. Our major focus was on P(K|ΔLk), the distribution of knot types for a particular value of ΔLk, which had not been evaluated previously. We found that unknotted circular DNA is not the most probable state beyond small values of ΔLk. Highly chiral knotted DNA has a lower free energy because it has less torsional deformation. Surprisingly, even at |ΔLk| > 12, only one or two knot types dominate the P(K|ΔLk) distribution despite the huge number of knots of comparable complexity. A large fraction of the knots found belong to the small family of torus knots. The relationship between supercoiling and knotting in vivo is discussed.

Invariant size–frequency distributions along a latitudinal gradient in marine bivalves

In the most extensive analysis of body size in marine invertebrates to date, we show that the size–frequency distributions of northeastern Pacific bivalves at the provincial level are surprisingly invariant in modal and median size as well as size range, despite a 4-fold change in species richness from the tropics to the Arctic. The modal sizes and shapes of these size–frequency distributions are consistent with the predictions of an energetic model previously applied to terrestrial mammals and birds. However, analyses of the Miocene–Recent history of body sizes within 82 molluscan genera show little support for the expectation that the modal size is an evolutionary attractor over geological time.

Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations

We describe and test a Markov chain model of microsatellite evolution that can explain the different distributions of microsatellite lengths across different organisms and repeat motifs. Two key features of this model are the dependence of mutation rates on microsatellite length and a mutation process that includes both strand slippage and point mutation events. We compute the stationary distribution of allele lengths under this model and use it to fit DNA data for di-, tri-, and tetranucleotide repeats in humans, mice, fruit flies, and yeast. The best fit results lead to slippage rate estimates that are highest in mice, followed by humans, then yeast, and then fruit flies. Within each organism, the estimates are highest in di-, then tri-, and then tetranucleotide repeats. Our estimates are consistent with experimentally determined mutation rates from other studies. The results suggest that the different length distributions among organisms and repeat motifs can be explained by a simple difference in slippage rates and that selective constraints on length need not be imposed.

The estimation of statistical parameters for local alignment score distributions

The distribution of optimal local alignment scores of random sequences plays a vital role in evaluating the statistical significance of sequence alignments. These scores can be well described by an extreme-value distribution. The distribution’s parameters depend upon the scoring system employed and the random letter frequencies; in general they cannot be derived analytically, but must be estimated by curve fitting. For obtaining accurate parameter estimates, a form of the recently described ‘island’ method has several advantages. We describe this method in detail, and use it to investigate the functional dependence of these parameters on finite-length edge effects.

Characteristics of the deep ocean carbon system during the past 150,000 years: ΣCO2 distributions, deep water flow patterns, and abrupt climate change

Studies of carbon isotopes and cadmium in bottom-dwelling foraminifera from ocean sediment cores have advanced our knowledge of ocean chemical distributions during the late Pleistocene. Last Glacial Maximum data are consistent with a persistent high-ΣCO2 state for eastern Pacific deep water. Both tracers indicate that the mid-depth North and tropical Atlantic Ocean almost always has lower ΣCO2 levels than those in the Pacific. Upper waters of the Last Glacial Maximum Atlantic are more ΣCO2-depleted and deep waters are ΣCO2-enriched compared with the waters of the present. In the northern Indian Ocean, δ13C and Cd data are consistent with upper water ΣCO2 depletion relative to the present. There is no evident proximate source of this ΣCO2-depleted water, so I suggest that ΣCO2-depleted North Atlantic intermediate/deep water turns northward around the southern tip of Africa and moves toward the equator as a western boundary current. At long periods (>15,000 years), Milankovitch cycle variability is evident in paleochemical time series. But rapid millennial-scale variability can be seen in cores from high accumulation rate series. Atlantic deep water chemical properties are seen to change in as little as a few hundred years or less. An extraordinary new 52.7-m-long core from the Bermuda Rise contains a faithful record of climate variability with century-scale resolution. Sediment composition can be linked in detail with the isotope stage 3 interstadials recorded in Greenland ice cores. This new record shows at least 12 major climate fluctuations within marine isotope stage 5 (about 70,000–130,000 years before the present).

Distinct desmocollin isoforms occur in the same desmosomes and show reciprocally graded distributions in bovine nasal epidermis.

The adhesive core of the desmosome is composed of cadherin-like glycoproteins of two families, desmocollins and desmogleins. Three isoforms of each are expressed in a tissue-specific and developmentally regulated pattern. In bovine nasal epidermis, the three desmocollin (Dsc) isoforms are expressed in overlapping domains; Dsc3 expression is strongest in the basal layer, while Dsc2 and Dsc1 are strongly expressed in the suprabasal layers. Herein we have investigated whether different isoforms are assembled into the same or distinct desmosomes by performing double immunogold labeling using isoform-specific antibodies directed against Dsc1 and Dsc3. The results show that individual desmosomes harbor both isoforms in regions where their expression territories overlap. Quantification showed that the ratio of the proteins in each desmosome altered gradually from basal to immediately suprabasal and upper suprabasal layers, labeling for Dsc1 increasing and Dsc3 decreasing. Thus desmosomes are constantly modified as cells move up the epidermis, with continuing turnover of the desmosomal glycoproteins. Statistical analysis of the quantitative data showed a possible relationship between the distributions of the two isoforms. This gradual change in desmosomal composition may constitute a vertical adhesive gradient within the epidermis, having important consequences for cell positioning and differentiation.