Ambiguous Moment Tensors and Radiation Patterns in Anisotropic Media with Applications to the Modeling of Earthquake Mechanisms in W-Bohemia

Anisotropic material properties are usually neglected during inversions for source parameters of earthquakes. In general anisotropic media, however, moment tensors for pure-shear sources can exhibit significant non-double-couple components. Such effects may be erroneously interpreted as an indication for volumetric changes at the source. Here we investigate effects of anisotropy on seismic moment tensors and radiation patterns for pure-shear and tensile-type sources. Anisotropy can significantly influence the interpretation of the source mechanisms. For example, the orientation of the slip within the fault plane may affect the total seismic moment. Also, moment tensors due to pure-shear and tensile faulting can have similar characteristics depending on the orientation of the elastic tensor. Furthermore, the tensile nature of an earthquake can be obscured by near-source anisotropic properties. As an application, we consider effects of inhomogeneous anisotropic properties on the seismic moment tensor and the radiation patterns of a selected type of micro-earthquakes observed in W-Bohemia. The combined effects of near-source and along-path anisotropy cause characteristic amplitude distortions of the P, S1 and S2 waves. However, the modeling suggests that neither homogeneous nor inhomogeneous anisotropic properties alone can explain the observed large non-double-couple components.The results also indicate that a correct analysis of the source mechanism, in principle, is achievable by application of anisotropic moment tensor inversion.     

Are open fractures necessarily aligned with maximum horizontal stress?

Fluid flow in fractured rock is an increasingly central issue in recovering water and hydrocarbon supplies and geothermal energy, in predicting flow of pollutants underground, in engineering structures, and in understanding large-scale crustal behaviour. Conventional wisdom assumes that fluids prefer to flow along fractures oriented parallel or nearly parallel to modern-day maximum horizontal compressive stress, or S_Hmax. The reasoning is that these fractures have the lowest normal stresses across them and therefore provide the least resistance to flow. For example, this view governs how geophysicists design and interpret seismic experiments to probe fracture fluid pathways in the deep subsurface. Contrary to these widely held views, here we use core, stress measurement, and fluid flow data to show that S_Hmax does not necessarily coincide with the direction of open natural fractures in the subsurface (>3 km depth). Consequently, in situ stress direction cannot be considered to predict or control the direction of maximum permeability in rock. Where effective stress is compressive and fractures are expected to be closed, chemical alteration dictates location of open conduits, either preserving or destroying fracture flow pathways no matter their orientation.     

Ar/Ar ages of tephras intercalated in astronomically tuned Neogene sedimentary sequences in the eastern Mediterranean

We present new ⁴⁰Ar/³⁹Ar data for sanidine and biotite derived from volcanic ash layers that are intercalated in Pliocene and late Miocene astronomically dated sequences in the Mediterranean with the aim to solve existing inconsistencies in the intercalibration between the two independent absolute dating methods. ⁴⁰Ar/³⁹Ar sanidine ages are systematically younger by 0.7-2.3% than the astronomical ages for the same ash layers. The significance of the discrepancy disappears except for the upper Ptolemais ashes, which reveal the largest difference, if an improved full error propagation method is applied to calculate the absolute error in the ⁴⁰Ar/³⁹Ar ages. The total variance is dominated by that of the activity of the decay of ⁴⁰K to ⁴⁰Ar (∼70%) and that the amount of radiogenic ⁴⁰Ar_p in the primary standard GA1550 biotite (∼15%). If the ⁴⁰Ar/³⁹Ar ages are calculated relative to an astronomically dated standard, the influence of these parameters is greatly reduced, resulting in a more reliable age and in a significant reduction of the error in ⁴⁰Ar/³⁹Ar dating.Astronomically calibrated ages for Taylor Creek Rhyolite (TCR) and Fish Canyon Tuff (FCT) sanidine are 28.53±0.02 and 28.21±0.04 Ma (±1 S.E.), respectively, if we start from the more reliable results of the Cretan A1 ash layer. The most likely explanation for the large discrepancy found for the younger Ptolemais ash layers (equivalent to FCT of 28.61 Ma) is an error in the tuning of this part of the sequence.     

Using Eigenpattern Analysis to Constrain Seasonal Signals in Southern California

— Earthquake fault systems are now thought to be an example of a complex nonlinear system (Bak, et al., 1987; Rundle and Klein, 1995). The spatial and temporal complexity of this system translates into a similar complexity in the surface expression of the underlying physics, including deformation and seismicity. Here we show that a new pattern dynamic methodology can be used to define a unique, finite set of deformation patterns for the Southern California Integrated GPS Network (SCIGN). Similar in nature to the empirical orthogonal functions historically employed in the analysis of atmospheric and oceanographic phenomena (Preisendorfer, 1988), the method derives the eigenvalues and eigenstates from the diagonalization of the correlation matrix using a Karhunen-Loeve expansion (KLE) (Fukunaga, 1970; Rundle et al., 2000; Tiampo et al., 2002). This KLE technique may be used to determine the important modes in both time and space for the southern California GPS data, modes that potentially include such time-dependent signals as plate velocities, viscoelasticity, and seasonal effects. Here we attempt to characterize several of the seasonal vertical signals on various spatial scales. These, in turn, can be used to better model geophysical signals of interest such as coseismic deformation, viscoelastic effects, and creep, as well as provide data assimilation and model verification for large-scale numerical simulations of southern California.     

Dynamic performance of twin curved cable‐stayed bridges

The dynamic behaviour of two curved cable‐stayed bridges, recently constructed in northern Italy, has been investigated by full‐scale testing and theoretical models. Two different excitation techniques were employed in the dynamic tests: traffic‐induced ambient vibrations and free vibrations. Since the modal behaviour identified from the two types of test are very well correlated and a greater number of normal modes was detected during ambient vibration tests, the validity of the ambient vibration survey is assessed in view of future monitoring. For both bridges, 11 vibration modes were identified in the frequency range of 0ndash;10Hz, being a one‐to‐one correspondence between the observed modes of the two bridges. Successively, the information obtained from the field tests was used to validate and improve 3D finite elements so that the dynamic performance of the two systems were assessed and compared based on both the experimental results and the updated theoretical models. Copyright © 2003 John Wiley & Sons, Ltd.     

Seismic tests on reinforced concrete and steel frames retrofitted with dissipative braces

Seismic tests have been conducted on two 3‐storey structures protected with pressurized fluid‐viscous spring damper devices. One of the structures was a reinforced concrete frame with clay elements in the slabs, while the other one was a steel frame with steel/concrete composite slabs. The spring dampers were installed through K bracing in between the floors. The tests were performed by means of the pseudodynamic method, which allowed the use of large and full‐size specimens, and by implementing a specific compensation strategy for the strain‐rate effect at the devices. The test results allowed the verification of the adequacy of the attachment system as well as the comparison of the behaviour of the unprotected buildings with several protected configurations, showing the benefits of the application of the devices and the characteristics of their performance. The response of the protected structures was always safer than that of the unprotected ones mainly due to a significant increase of equivalent damping. The increase in the damping ratio depends on the level of deformation. Copyright © 2004 John Wiley & Sons, Ltd.     

Bayesian algorithm for the estimation of the dielectric constant from active and passive remotely sensed data

An inversion technique based on the merging of microwave remotely sensed data is applied to ground-based radiometer and scatterometer data acquired for the same area. The purpose of this technique is to retrieve the dielectric constant of bare soils. The algorithm is based on a Bayesian approach and combines prior information on the dielectric constant and surface roughness with observed data, in order to obtain a marginal posterior probability density function. The function describes how the probability is distributed within the range of the dielectric constant values, given the measured values of emissivity and backscattering coefficient. The algorithm allows for the incorporation of all the available sources of information, such as multipolarization and multifrequency data. Several criteria, which have been used to compare the predicted and the observed values, show that for dielectric constant values higher than 10 the best performance is achieved when data with one polarization and one or two frequencies are exploited. For dielectric constant values of less than 10, the configuration with two polarizations produces the best estimates.