Seismic energy partitioning and scattering in laterally heterogeneous ocean crust

We present finite difference forward models of elastic wave propagation through laterally heterogeneous upper oceanic crust. The finite difference formulation is a 2-D solution to the elastic wave equation for heterogeneous media and implicitly calculatesP andSV propagation, compressional to shear conversion, interference effects and interface phenomena. Random velocity perturbations with Gaussian and self-similar autocorrelation functions and different correlation lengths (a) are presented which show different characteristics of secondary scattering. Heterogeneities scatter primary energy into secondary body waves and secondary Stoneley waves along the water-solid interface. The presence of a water-solid interface in the model allows for the existence of secondary Stoneley waves which account for much of the seafloor noise seen in the synthetic seismograms for the laterally heterogeneous models.Random incoherent secondary scattering generally increases aska (wavenumber,k, and correlation length,a) approaches one. Deterministic secondary scattering from larger heterogeneities is the dominant effect in the models aska increases above one. Secondary scattering also shows up as incoherence in the primary traces of the seisograms when compared to the laterally homogeneous case. Cross-correlation analysis of the initialP-diving wave arrival shows that, in general, the correlation between traces decreases aska approaches one. Also, because many different wave types exist for these marine models, the correlation between traces is range dependent, even for the laterally homogeneous case.     

Scattering ofP andS waves by a spherically symmetric inclusion

Scattering of an arbitrary elastic wave incident upon a spherically symmetric inclusion is considered and solutions are developed in terms of the spherical vector system of Petrashen, which produces results in terms of displacements rather than displacement potentials and in a form suitable for accurate numerical computations. Analytical expressions for canonical scattering coefficients are obtained for both the cases of incidentP waves and incidentS waves. Calculations of energy flux in the scattered waves lead to elastic optical theorems for bothP andS waves, which relate the scattering cross sections to the amplitude of the scattered fields in the forward direction. The properties of the solutions for a homogeneous elastic sphere, a sphere filled by fluid, and a spherical cavity are illustrated with scattering cross sections that demonstrate important differences between these types of obstacles. A general result is that the frequency dependence of the scattering is defined by the wavelength of the scattered wave rather than the wavelength of the incident wave. This is consistent with the finding that the intensity of thePS scattering is generally much stronger than theSP scattering. When averaged over all scattering angles, the mean intensity of thePS converted waves is2V _p ² /V _s ⁴ times the mean intensity of theSP converted waves, and this ratio is independent of frequency. The exact solutions reduce to simple and easily used expressions in the case of the low frequency (Rayleigh) approximation and the low contrast (Rayleigh-Born) approximation. The case of energy absorbing inclusions can also be obtained by assigning complex values to the elastic parameters, which leads to the result that an increase in attenuation within the inclusion causes an increased scattering cross section with a marked preference for scatteredS waves. The complete generality of the results is demonstrated by showing waves scattered by the earth's core in the time domain, an example of high-frequency scattering that reveals a very complex relationship between geometrical arrivals and diffracted waves.     

On the non-tidal secular acceleration of the Earth's rotation

Summary The secular positive acceleration of the Earth's rotation has been computed on the basis of the observed secular decrease of the second zonal harmonic[5]. It corresponds to the observed secular deceleration of the Earth's rotation which should be greater because of oceanic tides.

---

¶rt; u a¶rt;a u m aauu [5], u ¶rt;a num u au. m mmmu a¶rt;a u mu au u, m.. n¶rt; n m uu , ¶rt; m n auu nuua.

     

On Seasonal Variation of the 162 kHz Radio Wave Reflection Height: Model Calculations and Their Comparison with Experiments

The seasonal variation of the 162 kHz radio wave reflection height was measured indirectly at Panská Ves using the IPHA method described in Fier and Matys (1992). The fact that these waves are reflected at the height where the electron concentration is about 3.5 × 10 ⁸ m ^–3 was used to model this height using a 1-D model of the lower ionosphere electron concentration described by Ondráková (1993). The comparison of the measurements and model results indicates that: 1. The annual mean of the reflection height derived from measurements is several km higher than that derived from the model; 2. the character of the modelled seasonal variation is different from the real variation shown in Fier and Latovika (1992). The differences are discussed.     

Iterative geophysical tomography

Summary The algorithm of iterative geophysical tomography is presented. The medium is approximated smoothly by means of B-splines. The tww-point problem of ray computation is solved with the aid of paraxial approximation. The parameters of the medium are obtained from the iterative algorithm of minimizing the quadratic form. Two numerical 2-D examples are given.

---

u¶rt; au umamuuu mauu. ¶rt;a annuuaa n nu nu -na. ma na aa a nu nu naaua annuauu. aam ¶rt; n a umamu aua uuauauu a¶rt;amu . am nu¶rt; ¶rt;a 2-D u nua.

     

Secular changes of the geoid in the Fennoscandian land uplift region

Summary Two alternative methods are presented of determining the secular changes of the geoid in the Fennoscandian uplift region. Both are based on the Airy-Heiskanen hypothesis of isostatic compensation. The first alternative is the usual Stokes integration of changes of free-air anomalies in the given region. Another alternative was derived, which is based on calculating the gravitational potential of the surface uplift layer and of the same layer at the depth of compensation. The calculated changes of the shape of the geoid were then compared with the results of Bjerhammar's geophysical interpretation of the expansion of the Earth's gravity potential (regional depression of the geoid in the Fennoscandian uplift region). Regression analysis indicates a strong correlation between the two quantities (r=0.88). The time required to level this depression of the geoid is estimated. A relation was derived for improving the accuracy of the observed values of Fennoscandian uplifts, which takes into account the elastic reaction of the Earth's body to the change of gravitational potential and perturbations, caused by the secular decrease of the depth of the northern part of the Baltic Sea.

---

¶rt;ma ¶rt;a auama m¶rt;a n¶rt;u uuu¶rt;a amu n¶rt;a a¶rt;uu. a u¶rt;m uunm uamu nauu u-aa. auam — umuau ma uu aau ¶rt; ¶rt; ¶rt;a amu. ¶rt; m auam, a a uuuaumau nmuaa nm n¶rt;a, u a au nauu. am u uu u¶rt;a auam mamau aaa nuu umnmauu au nmuaa u mmu u (ua ¶rt;nuuu¶rt;a amu a¶rt;uauu). u aau aam a u au u uu (r=0,88). uam , ¶rt;u ¶rt; auau a ¶rt;nuuu¶rt;a. ¶rt; mu ¶rt; mu u au n¶rt; a¶rt;uauu, a n au u a uuaumau nmuaa u u, mau u uu amu amu .

     

A note on the general concept of wave breaking for Rossby and gravity waves

A recently proposed general definition of wave breaking is further discussed, in order to deal with some points on which misunderstanding appears to have arisen. As with surface and internal gravity waves, the classification of Rossby waves into breaking and not breaking is a generic classification based on dynamical considerations, and not a statement about any unique signature or automatically recognizable shape. Nor is it a statement about passive tracers uncorrelated with potential vorticity on isentropic surfaces. A strong motivation for the definition is that proofs of the nonacceleration theorem of wave, mean-flow interaction theory rely, explicitly or implicitly, on a hypothesis that the waves do not break in the sense envisaged.The general definition refers to the qualitative behaviour of a certain set of material contours, namely those, and only those, which would undulate reversibly, with small slopes, under the influence of the waves' restoring mechanism, in those circumstances for which linearized, nondissipative wave theory is a self-consistent approximation to nonlinear reality. The waves' restoring mechanism depends upon the basic-state vertical potential density gradient in the case of gravity waves, and upon the basic-state isentropic gradient of potential vorticity in the case of Rossby waves. In the usual linearized theory of planetary scale Rossby waves on a zonal shear flow, the relevant material contours lie along latitude circles when undisturbed.     

The gravitational effect of bodies with variable density

Summary Determination of the gravitational effect of some bodies, the density of which is supposed to be variable in the vertical direction, possibly in the horizontal direction, too.

---

n¶rt;uaumau ma m m, nmm m um mua uu muma anau.

     

Approximate solutions of the surface mass loading of a spherical Earth model

Summary Simple linear representation of the components of an approximate plane solution of point mass loading of the Earth's surface in a conveniently chosen coordinate system leads to selection of a 2nd-degree curve which is the best fit of the spherical solution for the given Earth model. The new approximate solution, which, analogously to the plane solution, can be called a parabolic solution, enables the simple input parameters of the plane solution to be used also for substantially larger angular distances. The comparison with the spherical solution is carried out by computing the effects of the M₂-wave of ocean tides. The results of the computations for the tidal station Brussels prove the two solutions to be in sufficient agreement for global problems as well.

---

m nu uau aa nuuum nm u m m au nmu u n¶rt;¶rt; um ¶rt;uam ¶rt;m ma u m mnu, ma auu a n¶rt;¶rt;um u u ¶rt; ¶rt;a ¶rt;u u. nuuum u, m n aauu nm u aam naauu u, nm unam nm ¶rt; ¶rt;a nm u ¶rt;a ¶rt; m u amu. au uu u m uu uu ₂ u nuu. mam uu ¶rt; nuu mauu ¶rt;aam ¶rt;mum au u u ¶rt;a ¶rt;a a¶rt;a.

     

Curvature of a group of difeomorphisms and the instability of Kármán's vortex street

a auu uu nmama uau u¶rt;a u¶rt;mu (nn ¶rt;uu au m a) ua na ááa m a muau¶rt;¶rt;uaua uma. mumu mmmm uu aa nu¶rt;umu n au uu ma a u ma mu.     

Time patterns of polarization currents generated in measuring the electrical conductivity of rocks

amamuu aau m¶rt;a uu mn¶rt;mu n¶rt; nu u mnama, m num , naa, m uau nuau nmu um a mam uu.     

The intensity of the main phase of geomagnetic storms in relation to the parameters of the solar wind

Summary On the basis of investigating 10 storms (1965–1967) good correlation was found between the density of the solar wind energy (₂=1/2mNv²) and the intensity of the main phase of the geomagnetic storms, expressed in terms of the maximum decrease of the horizontal intensity (B=H/cos). The relation between ₂, or Nv², and B could then be used to determine the quantities and ₀ ( is the factor expressing the increase in energy density in the magnetosphere, ₀ is the energy density of the particles in a quiet magnetosphere). A comparison with the directly observed distribution of the energy density of the particles in the magnetosphere indicates that the computed value of ₀ seems to be realistic. The magnitude of the factor will have to be checked again.     

Precession-nutation torque in terms of the stokes constants

um nu-mau m nuuau u auumu m nm ma u. m mu u¶rt;uu¶rt;a m ma nuua u uu.     

Friction constitutive law with rate and state dependences

A general formula for the Dieterich-Ruina friction constitutive law with rate and state (n-state variables,n=1, 2,...) dependences has been obtained and discussed under the assumption that the slip acceleration a varies ion a linearly with the slip displacement , namelya = a ₀ + (-₀). Wherea ₀, ₀ are initial constants, is the acceleration rate and constant.a ₀ and may be arbitrary constants (positive, negative or zero).The extreme value of frictional resistance and the existence condition of the extreme value, which are very important and govern to some degree the motion process of a frictionally slipping mechanical system, have been analyzed. A critical value _c which is the measure of the velocity weakening and velocity strengthening of the mechanical system, and its properties and the relationship to the extreme problem have been studied. Again, according to the critical value _c, the concepts of light or strong velocity weakening (or strengthening) are introduced.A possibly new phenomenon that frictional resistance may vary in some kind of decayed oscillation is found. Finally, the condition for the smallest frictional resistance for a slipping mechanical system with nonuniform acceleration has been obtained.     

Regularization of downward heat flow continuation

Summary Some regularization approaches to downward heat flow continuation from the Earth's surface to the depth of interest are presented, provided the heat transfer is described by the homogeneous steady-state equation. The Fourier transform of horizontal coordinates is used and the relation between the regularization of the heat flow spectrum and the regularization of the heat flow in the spatial domain are discussed. The utilization of an a priori information about the subsurface heat flow is mentioned.

---

¶rt; m zuau n¶rt;¶rt; n¶rt;u u nmz mnz nma n¶rt;nuu, m n mna nuam ¶rt;¶rt; maua au (au anaa). na naau zuma ¶rt;uam u ama ¶rt; zuau mnz nma u uau z nma. nma m nuu anu u n¶rt;nmz mnz nma.

     

A note on the stress-dilatancy relation for simulated fault gouge

Theoretical constraints on the stress-dilation relation for a deforming Coulomb material requirev ifC=0 andv sin^-1( _m / _m ) always, wherev is the dilation angle, is the friction angle,C is cohesion, _m is the maximum shear stress, and _m is the mean effective stress. Recent laboratory measurements of friction and dilatancy of simulated fault gouge show that small amplitude shear-load cycling causes compaction and consolidation. Comparison of the data with theory indicates that such load cycling produces: (1) increased coefficient of friction (or friction angle), (2) increased cohesion, and (3) increased dilatancy rate (or dilation angle). Under certain conditions of load cycling without significant plastic shear strain accumulation ( ^p <0.005) we find thatv exceeds both and, in contrast to theory, sin^-1( _m / _m ). This result is interpreted in terms of enhanced cohesion and overconsolidation, which lead to residual stresses within the gouge. An analogy is drawn between these special loading conditions and those extant on natural faults. In particular, our results imply that jostling and minor stress variations associated with microearthquakes may produce strengthening of fault gouge and changes in the fault zone's stress-dilatancy relation. Hence, compaction associated with microseismicity may lead to subsequent dilation of fault gouge, even for faults with large displacement rates and large net offsets (e.g., San Andreas). In regions where such dilation persists over sufficient displacements (on the order of the critical slip distance for seismic faulting) it may tend to inhibit unstable slip.     

Optimum determination of Euler Angles

a m m n¶rt; uua ¶rt; au anauu a nmu n¶rt;nuu, m u ¶rt;a a ¶rt;u -um, ma u amu um ¶rt;uam. a namu unm ¶rt;a anau, nu a¶rt;u m m aua n-a. nuam n¶rt; n¶rt;u au — D-nmuamu — naa uu anau ¶rt; n¶rt;u a.     

Explosive source and its influence on ground motion dynamics

Summary In adjusting measured values in sets A(r*), v(r*) and f(r*) by means of a power function in the form of P=Kr^* a region of discontinuity of the approximating curves was found at the distance r*11.5 m kg ^–1/3. It is assumed that this discontinuity was caused by the varying character of the source of seismic waves. For scaled distances r*>11.5 m kg ^–1/3 the explosion was considered to be a spherical source from the point of view of the charge geometry and of the distance of the pick-up from the centre of the charge, whereas if r*<11.5 m kg ^–1/3 the explosion in the borehole had the character of a cylindrical source. The difference of the two types of sources was reflected in the exponent with both the functions A(r*) and v(r*), so that for r*>11.5 m kg ^–1/3 –4.0 and–2.4, and for r*<11.5 m kg ^–1/3 –2.5 and–1.5. For the same intervals of scaled distance in the set f(r*)1.4 and1.2.     

On reflection and diffraction of pressure waves from floating ice

Summary The diffraction of a pressure wave arising from the edge of the floating ice due to an incident wave of the formH[t–(xcos/)+ (ysin/)], whereH(t) is the Heaviside step function, has been studied. The ice sheet has been taken to be moving with a velocity much less than the velocity of sound in the liquid. The problem has been solved using an approximate method of solving Wiener Hopf problems. Finally, the pressure field at a large distance from the edge has been derived in a closed form evaluating the integrals by the saddle point method.     

Constraints to the three-dimensional non-hydrostatic density distribution in the earth

Summary The paper is concerned with the properties of a density distribution within the Earth. A system of density parameter constraints involving Stokes' coefficients of the gravity field and the parameters describing the Earth's figure is derived. A density model, whose parameters fit these constraints, accounts for the fine structure of the gravity field and Earth's figure. Additional condition imposed on the average spherical density model are derived; they guarantee that the average spherical model is compatible with the 3-D density model.

---

aamuam ma an¶rt;u nmmu mu u. ¶rt;a uma u, auuau au naam nmm ¶rt;u; uma am m naamaumau n u naam nuau u u. mma ¶rt;, naam m ¶rt;mm mu auuau u, nm n¶rt;um uu umnmau m mmaumau n u u u. ¶rt; ma ¶rt;num u, aa¶rt;au a ¶rt; u ¶rt; nmmu, m nuam ¶rt;uu ¶rt; u ¶rt;u m nmm ¶rt;.