New induction sounding tested in Central Europe

The hourly data of nine geomagnetic observatories situated in Central Europe have been analyzed using the generalized magnetovariation (GMV) method designed recently for induction soundings of inhomogeneous media. In this method, impedance is one of transfer functions in the differential relation between spectra of the magnetic components and their derivatives. The peculiarity of this impedance is its correspondence to the magnetotelluric one estimated from the linear relations. Three transfer functions have been estimated simultaneously for data of geomagnetic observatories, using three different routines working in the period range from three hours up to two days. Noises in the source field components have been compared with noise in the estimated plane field divergence. The multivariate errors-in-variables method was used to extract spatially and temporally coherent geomagnetic field structure from the partially incoherent geomagnetic variations. This method allows estimating reliably impedances and gradient tippers for each observatory, taking into consideration the Earth’s sphericity. The obtained responses have been used for induction soundings and for detecting a deep inhomogeneity in the region.     

Application of horizontal-to-vertical (H/V) Fourier spectral ratio for analysis of site effect on rock (NEHRP-class B) sites in Taiwan

The frequency-dependent amplification for rock (NEHRP-class B) sites was studied using earthquake ground-motion database collected in Taiwan during implementation of the Taiwan Strong Motion Instrumentation Program. The database used includes several hundred records from earthquakes of M_L 4.0–7.3 occurred between 1993 and 2004. The characteristics of amplification were evaluated using the well-known technique of horizontal-to-vertical Fourier spectral ratio (H/V) of the S-wave phase [Lermo J, Chavez-Garcia FJ. Site effect evaluation using spectral ratios with only one station. Bull Seism Soc Am 1993;83:1574–94]. The study allows us to analyze peculiarities of rock sites amplification in Northern and Eastern Taiwan. It was suggested to divide the NEHRP-class B site amplification into four types based on frequency of maximum amplification and the shape of amplification function. The applicability of the technique was also checked for a few stiff and soft soil sites (NEHRP-classes D and E).     

A statistically significant signature of multi-decadal solar activity changes in atmospheric temperatures at three European stations

We study the effect of solar variability on temperatures recorded in three European stations with the longest gapless series available (Prague, Bologna and Uccle). Following a pattern recognition approach, we partition daily temperature “indices” (minimum, maximum and range) in two separate classes with respect to the level of solar activity (high H vs low L 11 year cycles). Using the two-sample Kolmogorov–Smirnov statistics, multiple shuffles of data, and other partitions, we demonstrate that the separation between the probability distribution functions of H and L temperatures is statistically significant and robust. We find that average annual variations for the H and L classes display common and site-specific patterns. For practically all series considered, differences between graphs of annual change for the two classes H and L are large (～1 °C). Solar activity accounts, at least in part, for the multi-decadal variations in temperature observed at these European sites in the past two centuries.     

QVOA techniques for fracture characterization

New computational techniques of QVOA analysis (Quality factor Versus Offset and Azimuth) for fracture characterization are developed. The techniques are applied to synthetic surface data of reflection with noise.     

Diamondiferous eclogites from Yakutia,Siberia: evidence for a diversity of protoliths

Major-element and REE compositions of 14 diamondiferous eclogites from the Udachnaya kimberlite in Yakutia, Siberia have been determined by electron microprobe and secondary ion mass spectrometer (SIMS). Based on previous clinopyroxene classification schemes (e.g., Taylor and Neal 1989), all of these eclogite xenoliths belong to Group B/C, although some of the garnet compositions and mineral REE abundances are inconsistent with the indicated groups. This demonstrates the inadequacy of the classification scheme based on African eclogites for application to Siberian samples. Because of the coarse grain size of the Udachnaya nodules, meaningful modal abundances could not be obtained. However, reconstructed REE compositions using various garnet: clinopyroxene ratios demonstrate relative insensitivity to changes in mode for common eclogitic assemblages. Many of these reconstructed REE compositions show LREE depletions. Some depletions are consistent with an origin (either directly or through partial melting) as normal or Type-I ocean floor basalt. Others, however, require material of eclogitic or pyroxenitic affinities to undergo partial melting; this facilitates the depletion of LREE while leaving the HREE at nearly original levels. Many of the eclogites of South Africa are consistent with a protolith of anomalous or Type II ocean floor basalt. This fundamental difference between the two regions is the likely cause of the inconsistencies with the chemicallybased classification.     

P‐wave stacking‐velocity tomography for VTI media

A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity V_P0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces.     

Bianchi-II Rotating World

Dynamical and kinematic properties of Bianchi-II cosmological models with rotation and expansion are investigated. Exact solutions of Einstein field equations are obtained which describe the evolution of a rotating Universe. Exact solutions of null, timelike and spacelike geodesics are constructed. Two new cosmological tests for rotating universes are discussed: cosmological lens effect and cosmological mirror effect.     

Impact of waves on the sea drag: measurements in the Baltic sea and a model interpretation

Measurements from the Baltic Sea and a wind-over-wave coupled model are used to study the wave impact on the sea drag. The study has been carried out for different wave conditions, namely a pure wind-sea, following-swell/ mixed sea and cross-swell/ mixed sea. Measurements reveal the fact that the sea drag is dependent on the sea-state. In stationary conditions and in the absence of severe cross-swell, swell reduces drag compared to wind-sea at the same wind speed. The cross-swell enhances the drag as compared to the following-swell case and the magnitude of the drag coefficient is increased with increasing the angle of swell propagation to the wind. It is shown that the agreement between the model results and measurements is good for pure wind-sea and stationary mixed-sea cases. Discrepancies occur at light winds, where most of the data represent pure swell conditions. During these pure swell conditions the data are characterized by a large variation of the drag coefficient. The variation is caused by mesoscale variability in the stress co-spectra, wind-cross-swell effects and nonstationarity in the wave and wind fields not represented in the model.     

Analysis of Thomsen parameters for finely layered VTI media

Since the work of Postma and Backus, much has been learned about elastic constants in vertical transversely isotropic (VTI) media when the anisotropy is due to fine layering of isotropic elastic materials. Nevertheless, there has continued to be some uncertainty about the possible range of Thomsen's anisotropy parameters ε and δ for such media. We use both Monte Carlo studies and detailed analysis of Backus' equations for both two- and three-component layered media to establish the results presented. We show that ε lies in the range ?3/8 ε ½[〈v²_p〉〈v^?2_p〉?1], for finely layered media; smaller positive and all negative values of ε occur for media with large fluctuations in the Lamé parameter λ in the component layers. We show that δ can also be either positive or negative, and that for constant density media, sign (δ) = sign (〈v^?2_p〉 ? 〈v^?2_s〉〈v²_s/v²_p〉). Monte Carlo simulations show that among all theoretically possible random media, positive and negative δ are equally likely in finely layered media. (Of course, the δs associated with real earth materials may span some smaller subset of those that are theoretically possible, but answering this important question is beyond our present scope.) Layered media having large fluctuations in λ are those most likely to have positive δ. This is somewhat surprising since ε is often negative or a small positive number for such media, and we have the general constraint that ε ? δ > 0 for layered VTI media. Since Gassmann's results for fluid-saturated porous media show that the mechanical effects of fluids influence only the Lamé parameter λ, not the shear modulus μ, these results suggest that small positive δ occurring together with small positive ε (but somewhat larger than δ) may be indicative of changing fluid content in a layered earth.     

Determination of Porosity and Permeability out of Dispersion Analyses of Borehole Wall Displacements

Waves carried in the borehole can be measured by the use of controlled sources inside the borehole and receivers located in the same borehole as well. This article gives the theoretical background and develops a new method to determine rock-physical parameters out of such controlled measurements of waves carried inside the borehole. Theoretical dispersion curves of the group velocities are matched to dispersion analyses of registrations obtained in the borehole. This inversion process intends to determine the in situ porosity and permeability of the rock.