(reported in December, 2001)
We have tested the impact of lateral viscosity variations in the top 200 km of the mantle on the interpretation of Fennoscandian postglacial uplift data (Martinec, Cadek and Fleitout, 2001) and presented an analytical form of the layer propagator matrix for the response of a locally incompressible, layered, linear-viscoelastic sphere assuming that the initial density stratification is parametrized by Darwin's law (Martinec, Thoma and Wolf, 2001). We have also studied the numerical problems associated with the viscoelastic response of the Earth to external loads in the time domain by means of the method of lines, (Hanyk et al., in press).
We have continued in dealing with the geodetic boundary-value problems. We particularly discussed the effect of the planar and spherical terrain models (Vanicek et al, 2001., Novak et al. 2001) and the ellipsoidal model of the geoid (Ardestani and Martinec, 2001) on precise geoid determination. The satellite gradiometric mission GOCE (to be launched in 2004) motivated us to solve the spherical gradiometric boundary-value problems to find the external Earth's gravitational potential (Martinec, 2001).
We have simulated subduction of slabs and concentrated on the dynamic effect of slab weakening due to grain size reduction at the phase boundary under the presence of trench migration in Cizkova et al. (submitted).
We have developed novel approaches to electromagnetic induction modeling for axisymmetric 2-D models of electrical conductivity (Martinec, Everett and Velimsky, submitted), and fully 3-D models (Velimsky and Martinec, to be submitted till the end of 2001), based on direct integration in the time domain.
Classical Boussinesq approximation of thermal convection was used by Vecsey and Matyska (2001) to demonstrate the ability of wavelet spectra to reveal the multiscale temporal dynamics of convecting systems. We have also employed the extended-Boussinesq convection models to extract the Bullen parameter (Matyska and Yuen, 2001; Matyska and Yuen, in press), which revealed non-adiabaticity of plumes and similarities between the Bullen parameter magnitudes obtained from tomographic models and the values yielded by thermal convection. The effects of adiabatic heating/cooling and viscous dissipation on 3-D rapidly rotating thermal convection was studied by Mistr et al. (in press). The applicability of new multidimensional wavelet analysis to various geophysical fields was studied by Yuen et al. (2001, and the paper in press) and by Vecsey et al.(2001).
References:Development of the algorithms for two-point ray tracing and travel-time interpolation in 3-D continued. The equations derived for the estimation of the average Lyapunov exponents, describing the ray chaos due to heterogeneities in the velocity model, have been applied to the construction of velocity models suitable for ray tracing and other high-frequency asymptotic methods. The designed algorithm of constructing velocity models have been tested on various 2-D and 3-D synthetic structures. Particular attention have been devoted to the resolution of seismic inversion techniques, model fitting, conversion and smoothing with the application of medium correlation functions and Sobolev scalar products.
The equations for calculation of the third and higher partial derivatives of travel time in both isotropic and anisotropic media have been derived. The equations for the linear paraxial approximation of the polarization vectors and for the variation of the polarization vectors with a velocity perturbation have been derived. Fermat's variational principle for anisotropic inhomogeneous media and its relation to the Finslerian metric and Hamiltonian have been studied.
Explicit equations for approximate linearized reflection/transmission coefficients at a generally oriented weak-contrast interface separating generally anisotropic media have been derived. A comparison of the reflection/transmission coefficients for dissipative and nondissipative media has been also performed. Various kinds of the coupling ray theory for weakly anisotropic models have been studied and compared with the exact solution derived for the "twisted crystal" model.
A new algorithm of the optimization of the shape of Gaussian beams and packets for prestack migrations in complex models have been designed, including smoothing of the beam parameters over the Hamiltonian hypersurface in the phase space.
The contemporary state of the seismic ray theory has been described by V. Cerveny in his monograph. Professor V. Cerveny received the Ernst Mach Award from the Academy of Sciences of the Czech Republic. K. Zacek received the Award of the Minister of Education, Youth and Sports, awarded annually to the five best university students in the Czech Republic.
References:Records of several Turkey earthquakes of 1999, obtained at broad-band seismic stations jointly operated by the Charles University and Patras University since 1997 in the Corinth Gulf, Greece, have been used to study dispersion of Love and Rayleigh waves. As compared with previous structural models for the territory of Greece, the dispersion data require significantly lower velocities in the uppermost crust, and smaller crustal thickness (Novotny et al., 2001).
In conjunction with an active litospheric structural experiment in Europe, CELEBRATION, four additional short profiles (of about 50 km long) were measured and interpreted, with a special emphasis on the upper-crustal structure in the earthquake swarm region of Western Bohemia (Malek et al., in press-a). The previously gathered controlled-source data (1989-91), and microearthquake data (1996-99) from Western Bohemia were re-interpreted with the aim to reveal anisotropy of the upper crust (Malek et al., in press-b). Methodical and practical aspects of the hypocentral location by deterministic grid search were solved (Jansky, 2000).
Focal mechanism retrieval based on amplitude spectra of complete waveforms and first-motion polarities, ASPO (Zahradnik et al., 2001), has been successfully validated on the 1999 Athens and the 2001 Skyros, Greece earthquakes (Zahradnik, submitted-a,b). Classical method of empirical Green's functions has been generalized for the practically important case of unequal focal mechanisms of the mainshock and aftershock, and the innovated method has been applied to inversion of the rupture nucleation point of the Athens 1999 earthquake (Plicka and Zahradnik, submitted). The finite-extent model of the Athens earthquake was then used to synthesize strong-ground motions by the newly developed combined deterministic-stochastic method, PEXT; the modeling results have been found in agreement with the observed peak acceleration values and empirical attenuation relations (Zahradnik and Tselentis, submitted).
A new 3D finite-difference code for the elastic wavefield modeling on spatially irregular (rectangular) grids was written and successfully tested. Attention was focused on efficiency; to that aim, the source- and path-effects are treated separately of the site effects by means of a hybrid fully 3D approach (Oprsal and Zahradnik, in press).
References: