Quasi-isotropic approximation of ray theory for anisotropic media

Wave propagation in weakly anisotropic inhomogeneous media is studied by the quasi-isotropic approximation of ray theory. The approach is based on the ray-tracing and dynamic ray-tracing differential equations for an isotropic background medium. In addition, it requires the integration of a system of two complex coupled differential equations along the isotropic ray.
The interference of the qS waves is described by traveltime and polarization corrections of interacting isotropic S waves. For qP waves the approach leads to a correction of the traveltime of the P wave in the isotropic background medium.
Seismograms and particle-motion diagrams obtained from numerical computations are presented for models with different strengths of anisotropy.
The equivalence of the quasi-isotropic approximation and the quasi-shear-wave coupling theory is demonstrated. The quasi-isotropic approximation allows for a consideration of the limit from weak anisotropy to isotropy, especially in the case of qS waves, where the usual ray theory for anisotropic media fails.     

Forecasting Aftershock Activity: 4. Estimating the Maximum Magnitude of Future Aftershocks

Izvestiya, Physics of the Solid Earth - In this paper, we consider the problem of forecasting the magnitude of the strongest aftershock starting from a certain instant of time in the future. This...     

Electrical structure of the upper mantle beneath Central Europe: Results of the CEMES project

In the years 2001–2003, we accomplished the experimental phase of the project CEMES by collecting long-period magnetotelluric data at positions of eleven permanent geomagnetic observatories situated within few hundreds kilometers along the south-west margin of the East European Craton. Five teams were engaged in estimating independently the magnetotelluric responses by using different data processing procedures. The conductance distributions at the depths of the upper mantle have been derived individually beneath each observatory. By averaging the individual cross-sections, we have designed the final model of the geoelectrical structure of the upper mantle beneath the CEMES region. The results indicate systematic trends in the deep electrical structure of the two European tectonic plates and give evidence that the electrical structure of the upper mantle differs between the East European Craton and the Phanerozoic plate of west Europe, with a separating transition zone that generally coincides with the Trans-European Suture Zone.     

The significance of isotope specific diffusion coefficients for reaction-transport models of sulfate reduction in marine sediments

Modeling isotopic signatures in systems affected by diffusion, advection, and a reaction which modifies the isotopic abundance of a given species, is a discipline in its infancy. Traditionally, much emphasis has been placed on kinetic isotope effects during biochemical reactions, while isotope effects caused by isotope specific diffusion coefficients have been neglected. A recent study by Donahue et al. (2008) suggested that transport related isotope effects may be of similar magnitude as microbially mediated isotope effects. Although it was later shown that the assumed differences in the isotope specific diffusion coefficients were probably overstated by one or two orders of magnitude (Bourg, 2008), this study raises several important issues: (1) Is it possible to directly calculate isotopic enrichment factors from measured concentration data without modeling the respective system? (2) Do changes in porosity and advection velocity modulate the influence of isotope specific diffusion coefficients on the fractionation factor α? (3) If one has no a priori knowledge whether diffusion coefficients are isotope specific or not, what is the nature and magnitude of the error introduced by either assumption? Here we argue (A) That the direct substitution of measured data into a differential equation is problematic and cannot be used as a replacement for a reaction-transport model; (B) That the transport related fractionation scales linearly with the difference between the respective diffusion coefficients of a given isotope system, but depends in a complex non-linear way on the interplay between advection velocity, and downcore changes of temperature and porosity. Last but not least, we argue that the influence of isotope specific diffusion coefficients on microbially mediated sulfate reduction in typical marine sediments is considerably smaller than the error associated with the determination of the fractionation factor.     

Modelling of global ocean tides with allowance for island effects

Summary Parameterization schemes for single islands, island chains and archipelagos are described. It is shown that their utilization within the framework of a global tidal model leads not only to local peculiarities in the spatial structure of tidal characteristics, but also to cardinal restructuring of the tidal pattern far beyond the islands and an appreciable increase (about 0.5·10¹²W) of the global tidal energy dissipation. The latter is commensurable by the order of magnitude with the imbalance between estimates of tidal energy dissipation in the World Ocean as a whole and on the continental shelf.

---

Modellierung der Gezeiten der weltmeere unter Berücksichtigung von Inseleffekten

Zusammenfassung Es werden parametrisierte Darstellungen von einzelnen Inseln, Inselketten und Inselgruppen beschrieben. Es wird gezeigt, daß ihre Anwendung bei einem globalen Gezeitenmodell nicht nur zu lokalen Besonderheiten in der geographischen Struktur der Gezeiten führt, sondern auch zur grundsätzlichen Restrukturierung des Gezeitenmusters weit über die Inseln hinaus und zu einer deutlichen Zunahme (etwa 0,5·10¹² W) des globalen Gezeiten-Energieverlusts für die M₂-Tide. Dieser hat etwa die Größenordnung des Unterschiedes zwischen dem Gezeiten-Energieverlust des gesamten Weltmeeres und dem auf dem Kontinentalschelf.

     

Application of the wavelet transform to automatic seismic signal detection

The paper is devoted to such a priority direction of digital seismology as the detection of seismic signals in a noisy medium. Various formulations of the problem of seismic event detection are analyzed. The reasons for the ineffectiveness of the classical formulation of the detection problem are revealed, and an adequate general formulation of the problem in the frequency-time space is proposed. A criterion for the automatic determination of the length of a seismic signal is proposed. The time-scale STA/LTA (Short Time Average to Long Time Average) detector, utilizing the expansion of a seismogram in a packet of wavelets, is proposed. The method eliminates the shortcomings of the standard scheme and enables the detection of seismic events with a low signal-to-noise ratio. The problems of the choice of the mother wavelet and the decomposition level of seismograms are discussed from the standpoint of physical properties of seismic wave fields.     

On the nature of short-period oscillations of the main Black Sea pycnocline, submesoscale eddies, and response of the marine environment to the catastrophic shower of 2012

Field studies performed at the Shirshov Institute of Oceanology, Russian Academy of Sciences (SIO RAS), Black Sea hydrophysical polygon in 2012 are illustrated. The variations in the vertical distribution of the hydrophysical characteristics (water temperature, salinity, and density, as well as current velocity) in the upper 200-m layer of the Black Sea above the continental slope in the cold season, obtained using an Aqualog autonomous profiler on a moored buoy station, have been analyzed. It has been established that the position of the permanent pycno-halocline and the hydrosulphuric zone upper boundary intensively oscillate with a characteristic period of 5–10 days. These oscillations cause short-period variations in the thickness of the oxigenated layer by 20–40 m, which reaches one-third of the total thickness of the layer. Measurements performed with autonomous stations (bottom ADCP, thermochain) at the experimental subsatellite polygon in the Gelendzhik coastal zone, as well as meteorological, ship, and satellite data obtained during the catastrophic rains and flooding on July 6–7, 2012, and afterward, have been simultaneously analyzed. It has been established that a catastrophic flow of turbid fresh water into the sea caused the formation of a belt of freshened (by 1.0–2.7 psu) less dense water with a high suspension concentration on the shelf and the upper continental slope. This water formed a quasi-geostrophic northwestward along-shore current, the velocity of which reached 40–50 cm/s. Therefore, the freshened and turbid water mostly escaped from the Gelendzhik region northwestward for two days after the flood, and the remaining water became free of suspension owing to its settlement during approximately the same period. The fields of the current velocity and suspension concentration in a submesoscale cyclonic eddy, identified on the satellite image, were measured at the hydrophysical polygon. It has been established that a high (when compared to the background values) suspension concentration in the surface-water layer in an eddy is related to intense upwelling at the eddy center and the rising of suspension (apparently phytoplankton) from the thermocline layer, where the suspension concentration is maximal.     

An iteratively reweighted least-squares approach to adaptive robust adjustment of parameters in linear regression models with autoregressive and <Emphasis Type="Italic">t</Emphasis>-distributed deviations

In this paper, we investigate a linear regression time series model of possibly outlier-afflicted observations and autocorrelated random deviations. This colored noise is represented by a covariance-stationary autoregressive (AR) process, in which the independent error components follow a scaled (Student’s) t-distribution. This error model allows for the stochastic modeling of multiple outliers and for an adaptive robust maximum likelihood (ML) estimation of the unknown regression and AR coefficients, the scale parameter, and the degree of freedom of the t-distribution. This approach is meant to be an extension of known estimators, which tend to focus only on the regression model, or on the AR error model, or on normally distributed errors. For the purpose of ML estimation, we derive an expectation conditional maximization either algorithm, which leads to an easy-to-implement version of iteratively reweighted least squares. The estimation performance of the algorithm is evaluated via Monte Carlo simulations for a Fourier as well as a spline model in connection with AR colored noise models of different orders and with three different sampling distributions generating the white noise components. We apply the algorithm to a vibration dataset recorded by a high-accuracy, single-axis accelerometer, focusing on the evaluation of the estimated AR colored noise model.