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.     

Gravity torques for surface waves

The effects of the gravity torques acting on the angular momentum of surface gravity waves are calculated theoretically. For short crested waves the gravity torque is caused by the force of gravity on the orbiting fluid particles acting down the slopes of the crests and troughs and in the direction parallel to the crests and troughs. The gravity torque tries to rotate the angular momentum vectors, and thus the waves themselves, counterclockwise in the horizontal plane, as viewed from above, in both hemispheres. The amount of rotation per unit time is computed to be significant assuming reasonable values for the along-crest and trough slopes for waves in a storm area. The gravity torque has a frequency which is double the frequency of the waves. For long crested waves the gravity torque acts in the vertical plane of the orbit and tries to decelerate the particles when they rise and accelerate them when they fall. By disrupting the horizontal cyclostrophic balance of forces on the fluid particles (centrifugal force versus pressure force) the gravity torque accounts qualitatively for the three characteristics of breaking waves: that they break at the surface, that they break at the crest, and that the crest breaks in the direction of wave propagation.     

SII型海洋重力仪稳定平台颤动故障与排除方法

S型海洋重力仪是目前海洋重力测量中较为常用的一种仪器，其产品性能也逐步得到提高。介绍了进行全自动化控制升级后，最新引进的SⅡ型海洋重力仪使用中出现的稳定平台颤动故障现象及排除办法。     

任意曲线曲面上重磁位场转换的B样条函数法

提出用B样条函数求解曲线、曲面上重磁位场的向上延拓,水平、垂向导数计算,磁异常分量互换的方法。该方法的特点是:原理简明,程序通用性强,计算精度高。     

Peak discharge of a Pleistocene lava-dam outburst flood in Grand Canyon, Arizona, USA

The failure of a lava dam 165,000 yr ago produced the largest known flood on the Colorado River in Grand Canyon. The Hyaloclastite Dam was up to 366 m high, and geochemical evidence linked this structure to outburst-flood deposits that occurred for 32 km downstream. Using the Hyaloclastite outburst-flood deposits as paleostage indicators, we used dam-failure and unsteady flow modeling to estimate a peak discharge and flow hydrograph. Failure of the Hyaloclastite Dam released a maximum 11 × 10⁹ m³ of water in 31 h. Peak discharges, estimated from uncertainty in channel geometry, dam height, and hydraulic characteristics, ranged from 2.3 to 5.3 × 10⁵ m³ s⁻¹ for the Hyaloclastite outburst flood. This discharge is an order of magnitude greater than the largest known discharge on the Colorado River (1.4 × 10⁴ m³ s⁻¹) and the largest peak discharge resulting from failure of a constructed dam in the USA (6.5 × 10⁴ m³ s⁻¹). Moreover, the Hyaloclastite outburst flood is the oldest documented Quaternary flood and one of the largest to have occurred in the continental USA. The peak discharge for this flood ranks in the top 30 floods (>10⁵ m³ s⁻¹) known worldwide and in the top ten largest floods in North America.     

Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

A parametric study on the impact of satellite attitude errors on GOCE gravity field recovery

In the recent design of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission, the gravity gradients are defined in the gradiometer reference frame (GRF), which deviates from the actual flight direction (local orbit reference frame, LORF) by up to 3–4°. The main objective of this paper is to investigate the effect of uncertainties in the knowledge of the gradiometer orientation due to attitude reconstitution errors on the gravity field solution. In the framework of several numerical simulations, which are based on a realistic mission configuration, different scenarios are investigated, to provide the accuracy requirements of the orientation information. It turns out that orientation errors have to be seriously considered, because they may represent a significant error component of the gravity field solution. While in a realistic mission scenario (colored gradiometer noise) the gravity field solutions are quite insensitive to small orientation biases, random noise applied to the attitude information can have a considerable impact on the accuracy of the resolved gravity field models.