Cavity expansion in cross‐anisotropic rock

Schistous rock can be considered—in a first approximation—as cross‐anisotropic linear elastic material. The determination of the corresponding material constants on the basis of the laboratory investigation of rock samples often fails, as the extraction of appropriate cores proves to be unfeasible (the cores disintegrate if the schistosity is pronounced). In this paper a new method is presented to determine the material constants of a linear elastic cross‐anisotropic rock on the basis of cavity expansion field tests, e.g. with a radial jack. To this purpose, an analytic approximation for the deformation of a hydrostatically loaded cylindrical cavity in cross‐anisotropic rock is derived which serves to the inverse analysis of the material parameters. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison of ten notable meteorological drought indices on tracking the effect of drought on streamflow

ABSTRACT

Ten notable meteorological drought indices were compared on tracking the effect of drought on streamflow. A 730-month dataset of precipitation, temperature and evapotranspiration for 88 catchments in Oregon, USA, representing pristine conditions, was used to compute the drought indices. These indices were correlated with the monthly streamflow datasets of the minimum, maximum and mean discharge, and the discharge monthly fluctuation; it was revealed that the 3-month Z-score drought index (Z3) has the best association with the four streamflow variables. The Mann-Kendall trend detection test applied to the latter index time series mainly highlighted a downward trend in the autumn and winter drought magnitude (DM) and an upward trend in the spring and summer DM (p = 0.05). Finally, the Pettitt test indicated an abrupt decline in the annual and autumn DM, which began in 1984 and 1986, respectively.     

The temperature structure in a stably stratified internal boundary layer over a cold sea

Data from the Öresund experiment are used to investigate the structure of the stably stratified internal boundary layer (SIBL) which develops when warm air is advected from a heated land surface over a cooler sea. The present study is based on a theory developed by Stull (1983a, b, c). He proposed that the turbulence and the mean structure of the nocturnal boundary layer is controlled by the time-integrated value of surface heat flux and that the instantaneous heat flux is of less importance.Dimensional arguments are used to define simple, physically consistent, temperature, velocity and length scales. The dimensionless surface heat flux has a high value immediately downwind of the shoreline and it decreases rapidly in magnitude with increasing distance from the coast. Farther away, it is essentially constant. The dimensionless potential temperature change exhibits an exponential profile. It is estimated that turbulence accounts for 71% of boundary-layer cooling while clear-air radiational cooling is responsible for the remaining 29%.Finally it is found that theoretical predictions for the height of the SIBL are in a good agreement with observations.     

Determination of S-wave velocity structure using microtremors and spac method applied in Thessaloniki (Greece)

Array measurements of microtremors at 16 sites in the city of Thessaloniki were performed to estimate the V_s velocity of soil formations for site effect analysis. The spatial autocorrelation method was used to determine phase velocity dispersion curves in the frequency range from 0.8–1.5 to 6–7 Hz. A Rayleigh wave inversion technique (stochastic method) was subsequently applied to determine the Vs profiles at all the examined sites. The determination of V_s profiles reached a depth of 320 m. Comparisons with V_s values from cross-hole tests at the same sites proved the reliability of the SPAC method. The accuracy of the V_s profiles, the ability to reach large penetration depths in densely populated urban areas and its low cost compared to conventional geophysical prospecting, make Mictrotremor Exploration Method very attractive and useful for microzonation and site effects studies. An example of its application for the site characterization in Thessaloniki is presented herein.     

Analysis of fault rupture propagation through uniform soil cover

This paper presents results from numerical simulations of the propagation of an active dip–slip fault rupture through a uniform soil layer covering the rigid bedrock. Following verification of the numerical methodology against field evidence, a parametric study is performed for loose and dense sand, for normally consolidated and overconsolidated clay, as well as for different fault dip angles (normal and reverse faults) and for different thicknesses of the soil cover. The soil is modeled as an elasto-plastic, strain-softening material that obeys the Mohr–Coulomb failure criterion. The study aims at establishing criteria for the approximate depiction of the location and the width of the zone with significant ground surface distortion, where the differential ground displacements induced by the fault rupture may threaten the integrity of man-mad structures.     

Adaptive filtering of GOCE-derived gravity gradients of the disturbing potential in the context of the space-wise approach

Filtering and signal processing techniques have been widely used in the processing of satellite gravity observations to reduce measurement noise and correlation errors. The parameters and types of filters used depend on the statistical and spectral properties of the signal under investigation. Filtering is usually applied in a non-real-time environment. The present work focuses on the implementation of an adaptive filtering technique to process satellite gravity gradiometry data for gravity field modeling. Adaptive filtering algorithms are commonly used in communication systems, noise and echo cancellation, and biomedical applications. Two independent studies have been performed to introduce adaptive signal processing techniques and test the performance of the least mean-squared (LMS) adaptive algorithm for filtering satellite measurements obtained by the gravity field and steady-state ocean circulation explorer (GOCE) mission. In the first study, a Monte Carlo simulation is performed in order to gain insights about the implementation of the LMS algorithm on data with spectral behavior close to that of real GOCE data. In the second study, the LMS algorithm is implemented on real GOCE data. Experiments are also performed to determine suitable filtering parameters. Only the four accurate components of the full GOCE gravity gradient tensor of the disturbing potential are used. The characteristics of the filtered gravity gradients are examined in the time and spectral domain. The obtained filtered GOCE gravity gradients show an agreement of 63–84 mEötvös (depending on the gravity gradient component), in terms of RMS error, when compared to the gravity gradients derived from the EGM2008 geopotential model. Spectral-domain analysis of the filtered gradients shows that the adaptive filters slightly suppress frequencies in the bandwidth of approximately 10–30 mHz. The limitations of the adaptive LMS algorithm are also discussed. The tested filtering algorithm can be connected to and employed in the first computational steps of the space-wise approach, where a time-wise Wiener filter is applied at the first stage of GOCE gravity gradient filtering. The results of this work can be extended to using other adaptive filtering algorithms, such as the recursive least-squares and recursive least-squares lattice filters.     

Incompatible deformation in rock mechanics

Compatibility, i.e., continuous displacement field, is often taken for granted but experience teaches that rock deformation is in most cases discontinuous. The relation between incompatible deformation and continuous distribution of dislocations is established since long and is recalled in this paper. Contrary to several applications of Cosserat theory to already laminated (jointed) rock, the main goal of this paper is to give an explanation for the formation of discontinuities (joints) in rock as a result of a particular type of buckling which leads the body outside of the Euclidean space of its original configuration.     

Recursive algorithms for the computation of the potential harmonic coefficients of a constant density polyhedron

The gravitational potential of a constant density general polyhedron can be expressed both in terms of a closed analytical expression and as a series expansion involving the corresponding spherical harmonic coefficients. The latter can be obtained from two independent algorithms, which differ not only in their algorithmic architecture but in their efficiency and overall performance, especially when computing the coefficients of higher degree and order. In the present paper a comparative study of all these three approaches is carried out focusing on the numerical implementation of the recursive relations appearing in the two algorithms for the computation of the polyhedral potential harmonic coefficients. The performed numerical investigations show that the linear algorithm proposed by Jamet and Thomas (Proceedings of the second international GOCE user workshop, ‘GOCE, The Geoid and Oceanography’, ESA-ESRIN, Frascati, Italy, 8–10 March 2004, ESA SP-569, 2004), but so far not implemented, achieves a reasonable accuracy at a computational expense that opens to practical applications, for instance in the field of satellite gravimetry/gradiometry interpretation. The convergence behavior of the linear recursion algorithm is studied thoroughly and a computational procedure is proposed that enables the stable computation of potential harmonic coefficients up to degree 60 when referring to an arbitrarily shaped polyhedral body.