Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

Geothermal palaeogradients and metamorphic zonation from the conodont colour alteration index (CAI)

The conodont colour alteration index (CAI) is potentially valuable in thermal history investigations, but there are problems in its use. An approach is proposed herein to overcome the problem of establishing reliable geothermal palaeogradients in areas with burial metamorphism, based on knowledge of the stratigraphic succession thickness, the boundaries of chronostratigraphic units, and the age and duration of deformation episodes. Although also problematical that a correlation between CAI and Kübler index (KI) values does not exist, there are areas where CAI and KI values are notably consistent for the anchizone boundaries; however, anomalous correlations between CAI and KI values have also been found in some areas. A CAI anchizone (or ancaizone) is defined in order to enable CAI values to be used independently to establish a metamorphic zonation.     

Shale basins,sulfur-deficient ore brines and the formation of exhalative base metal deposits

The massive sulfide deposits of the Iberian Pyrite Belt are interbedded with felsic volcanic rocks and shale, and underlain by several thousand meters of siliciclastic sedimentary rocks known as the PQ Group. Isotope geochemistry and regional geology are both consistent with equilibration of the ore-forming fluids with the PQ Group, prior to ore deposition near the former seafloor. The average Cu:Zn:Pb ratio of the PQ Group rocks (ca. 26:55:19) is similar to the weighted average of all the massive sulfide orebodies combined (ca. 25:52:23).The genetic relationship between massive sulfide deposits and a siliciclastic sedimentary metal source is explained here by a thermodynamic model, proposing that mildly reducing redox conditions imposed by equilibration with the sedimentary rocks are most critical for the formation of an effective ore-forming fluid. Relatively metal-rich but organic-poor pyrite-bearing shale undergoing dewatering of saline pore fluids is an effective source for the generation of sulfur-deficient but relatively iron and base metal-rich brines. Thus, we propose that the giant deposits of the Iberian Pyrite Belt owe their existence not to exceptionally metal-enriched (e.g., magmatic) fluids, but to the existence of a fairly ordinary but large metal source in reactive siliciclastic sediments, combined with an underlying igneous heat source and a particularly efficient mechanism of sulfide precipitation by mixing with H₂S-rich fluids at or near the seafloor.Essentially similar mineral equilibria are imposed when saline fluids are buffered by typical continental basement rocks. Leaching of retrograde minerals and possibly residual salts from their magmatic or metamorphic prehistory is expected to generate similar, variably metal-rich but relatively sulfide-deficient fluids. Thus, the existence of mildly reducing rocks can be the dominant chemical control in the source of fluids generating many volcanogenic, Irish-type or sedex deposits, many of which are known to precipitate their metal load in response to biogenic sulfide addition at the ore deposition site.     

Burial and scour around short cylinder under progressive shoaling waves

The results of a laboratory experimental program aimed at better understanding the scour around and burial of heavy cylindrical objects under oscillating flow on a sandy bed are described. This study was motivated by its application to the dynamics of isolated cobbles/mines on a sandy floor under nonlinear progressive waves, such as that occur in shallow coastal waters beyond the wave-breaking region. In the experiments, nonlinear progressive waves were generated in a long wave tank of rectangular cross-section with a bottom slope. Model mines (short cylinders) were placed on the sandy bottom and the temporal evolution of the bed profile and the velocity field in the near field of the object were observed. Experiments were conducted at relatively high Reynolds numbers for a range of flow conditions, which can be characterized by the Keulegan–Carpenter number and Shields parameter. Depending on the values of these parameters, four different scour regimes around the cylinder including periodical burial of cylinder under migrating sand ripples were observed; they were classified as: (i) no scour/burial, (ii) initial scour, (iii) expanded scour, and (iv) periodic burial cases. A scour regime diagram was developed and the demarcation criteria between different regimes were deduced. Semi-empirical formulae that permit estimation of the scour depth with time, the equilibrium maximum scour depth and length, and conditions necessary for the burial of the cylinder as a function of main external parameters are also proposed.     

A probability distribution for depth-limited extreme wave heights in a sea state

A model for the depth-limited distribution of the highest wave in a sea state is presented. The distribution for the extreme wave height is based on a probability density function (pdf) for depth-limited wave height distribution for individual waves [Méndez, F.J., Losada, I.J., Medina, R. 2004. Transformation model of wave height distribution. Coastal Eng, Vol. 50, 97:115.] and considers the correlation between consecutive waves. The model is validated using field data showing a good representation of the extreme wave heights in the surf zone. Some important statistical wave heights are parameterized obtaining useful expressions that can be used in further calculations.     

The center-to-limb variations of four Ca i lines in the photospheric spectrum at λ6500

We study the center-limb (CL) variation of the average profiles of four Ca i lines near 6500 and compare these observations with synthetic data obtained from several line formation models having different thermal structures, line parameters, LTE and non-LTE conditions, and micro and macroturbulence values, to assess the formation characteristics of our Ca i lines in the solar photosphere.Comparison of numerical results with observations indicates that non-LTE is indispensable to fit the CL variation of the central residual intensity for the line 6493, and anisotropic microturbulence is indispensable to improve the CL behavior of the equivalent widths for all lines. The Ca i line analysis favors a cool photospheric model, but this cannot be disentangled clearly from the effects of non-LTE and small-scale velocity fields on the grounds of the present line formation models.     

Accuracy assessment of lidar-derived digital elevation models

Despite the relatively high cost of airborne lidar-derived digital elevation models (DEMs), such products are usually presented without a satisfactory associated estimate of accuracy. For the most part, DEM accuracy estimates are typically provided by comparing lidar heights against a finite sample of check point coordinates from an independent source of higher accuracy, supposing a normal distribution of the derived height differences or errors. This paper proposes a new methodology to assess the vertical accuracy of lidar DEMs using confidence intervals constructed from a finite sample of errors computed at check points. A non-parametric approach has been tested where no particular error distribution is assumed, making the proposed methodology especially applicable to non-normal error distributions of the type usually found in DEMs derived from lidar. The performance of the proposed model was experimentally validated using Monte Carlo simulation on 18 vertical error data-sets. Fifteen of these data-sets were computed from original lidar data provided by the International Society for Photogrammetry and Remote Sensing Working Group III/3, using their respective filtered reference data as ground truth. The three remaining data-sets were provided by the Natural Environment Research Council's Airborne Research and Survey Facility lidar system, together with check points acquired using high precision kinematic GPS. The results proved promising, the proposed models reproducing the statistical behaviour of vertical errors of lidar using a favourable number of check points, even in the cases of data-sets with non-normally distributed residuals. This research can therefore be considered as a potentially important step towards improving the quality control of lidar-derived DEMs.     

The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem

This paper deals with the problem of detecting and correcting cycle-slips in Global Navigation Satellite System (GNSS) phase data by exploiting the Bayesian theory. The method is here applied to undifferenced observations, because repairing cycle-slips already at this stage could be a useful pre-processing tool, especially for a network of permanent GNSS stations. If a dual frequency receiver is available, the cycle-slips can be easily detected by combining two phase observations or phase and range observations from a single satellite to a single receiver. These combinations, expressed in a distance unit form, are completely free from the geometry and depend only on the ionospheric effect, on the electronic biases and on the initial integer ambiguities; since these terms are expected to be smooth in time, at least in a short period, a cycle-slip in one or both the two carriers can be modelled as a discontinuity in a polynomial regression. The proposed method consists in applying the Bayesian theory to compute the marginal posterior distribution of the discontinuity epoch and to detect it as a maximum a posteriori (MAP) in a very accurate way. Concerning the cycle-slip correction, a couple of simultaneous integer slips in the two carriers is chosen by maximazing the conditional posterior distribution of the discontinuity amplitude given the detected epoch. Numerical experiments on simulated and real data show that the discontinuities with an amplitude 2 or 3 times larger than the noise standard deviation are successfully identified. This means that the Bayesian approach is able to detect and correct cycle-slips using undifferenced GNSS observations even if the slip occurs by one cycle. A comparison with the scientific software BERNESE 5.0 confirms the good performance of the proposed method, especially when data sampled at high frequency (e.g. every 1 s or every 5 s) are available.