The gross features of the lithosphere-asthenosphere system in Europe from seismic surface waves and body waves

Long-period recordings of dispersive Rayleigh waves along numerous station lines, or profiles, in Europe have for the first time permitted a uniform inversion of these observations based on a new method of phase velocity regionalization.Regional dispersion relations obtained by this method have then been subjected to a complete inversion procedure commonly known as the hedgehog method. The results are presented in a map outlining the thickness of the lower lithosphere (lid) and the shear (S) velocities in both the lid and the asthenosphere channel.A comparison of these results with the minimum compressional (P) wave velocities in the asthenosphere and their corresponding depths provides an estimate of theV _p /V _s ratio for the asthenosphere in the European area.Contribution No. 314, Institute of Geophysics, ETH-Zürich, Switzerland.     

Large-scale structure of magnetospheric plasma

Recent investigations of magnetospheric plasma structure are summarized under the broad categories of empirical models, transport across boundaries, formation, and dynamics of the plasma sheet. This report reviews work in these areas during the period 1991 to 1993. Fully three-dimensional empirical models and simulations have become important contributors to our understanding of the magnetospheric system. Some new structural concepts have appeared in the literature: the entry boundary and geopause, the plasma sheet region 1 vortices, the low-energy layer, the adiabaticity boundary or wall region, and a region in the tail to which we refer as the injection port. Traditional structural concepts have also been the subject of recent study, notably the plasmapause, the magnetopause, and the plasma sheet. Significant progress has been made in understanding the nature of plasma sheet formation and dynamics, but the acceleration of electrons to high energy remains somewhat mysterious.     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

Seismic signals in geothermal areas of active volcanism: a case study from ‘La Fossa’, Vulcano (Italy)

Since the last eruption (1888–1890) volcanism at Vulcano, Aeolian Archipelago, southern Tyrrhenian Sea, has taken the form of persistent fumarolic activity. The gas-vapour phases of the geothermal systems are mainly discharged within two restricted areas about 1 km apart from each other, in the northern part of the island. These areas are La Fossa crater, and the beach fumaroles of the Baia di Levante. Fluids released at the two main fumarolic fields display quite different chemical and temperature characteristics, implying different origins. The local seismicity essentially takes the form of discrete shocks of shallow origin (depth1 km) at La Fossa, usually with energy < 10¹³ ergs. They are thought to be related to the uprise of pressurized hot gases and vapours discharged at the crater fumaroles. The present investigation points to the existence of two principal categories of seismic events (called M-shocks and N-shocks). These are short events (normally < 10 s). M-type shocks are thought to be due to resonance vibrations within the interior of the volcano, probably driven by the excitation of shock-waves within cavities deeply affected by deposition and alteration of self-sealant hydrothermal minerals. N-type events display features that resemble those of volcano-tectonic earthquakes, but have no recognizable S-phases. Here they are tentatively attributed to microfracturing of rocks which have been extensively hydrothermally altered. Results of the present study permit a preliminary conceptual model of the local shallow seismic processes in the framework of geochemical modelling of fumarolic activity and geological inferences from geothermal drilling.     

Results of concentration and mobility measurements of positively and negatively charged particles taken by a rocket-borne parachuted aspiration (Gerdien) probe in the height region from 72 to 39 km

Summary With a self-neutralizing, rocket-borne parachute aspiration probe concentration and mobility measurements of positively and negatively charged carriers were made between 72 and 40 km over Sardinia. Between 72 km and 60 km a heavy and a light group of positively charged particles were found. The heavy group was more abundant than the hight group in the greater heights. For the negatively charged particles only one mobility group was found which roughly corresponded to the mobility of the light group of the positively charged carriers. The probe was designed so that no electrons could enter the analyzer section of the aspiration system. By this the electron concentration was given by the difference in concentration between the positively and negatively charged particles.This paper was read byH. Dolezalek, in an abbreviated form supplied by the authors.     

Intercomparison of New Zealand filter ozonometer and Dobson spectrophotometer total ozone measurements

A commercial version of the prototype New Zealand narrowband interference filter instrument has been run at Garmisch-Partenkirchen since autumn 1977. The vicinity of the Dobson stations Hohenpeissenberg and Arosa permitted a check of the filter instrument's utility, reliability and long-term stability in total ozone measurement by (i) station to station-or short-range intercomparisons of daily means and (ii) a several days' direct intercomparison at Arosa in spring 1978.The regression analysis with the Hohenpeissenberg data-covering the five month autumn-winter 1977/78 period with highly variable ozone-shows that the filter instrument's direct-sun modeX _AD total ozone values were systematically too low by a conversion factor of 0.93. The excellent proportionality between the instruments is indicated by the standard deviation of only 1,3%. The (–7.0±1.3%) bias perfectly agrees with the result reported byBasher [1] for a similar, yet direct, intercomparison.During the direct Arosa intercomparison the filter instrument also tracked down short-time variations in a remarkable, though biased, parallelism with the Dobson instruments, and theX _AD daily means fit well into the prognostic relationship derived from the preceding autumn-winter short-range intercomparison.A creeping deterioration of the characteristics of the 305.5 nm double filter then prevented any furtherX _AD total ozone intercomparison. No final decision is possible whether this behaviour after a half year's run is typical or just an unusual exception. The otherwise excellent long-term stability is demonstrated by the comparison of theX _CD daily means with the HohenpeissenbergX _AD data for a 1 1/2 year's period (autumn 1977–spring 1979), resulting in the same bias relationship as for theX _AD total ozone measurements.     

Numerical modeling of the high-latitude F-layer anomalies

The results of a detailed numerical study of the behaviour of the convective polar ionosphere are presented. The developed theoretical model produces three-dimensional distribution of electron density, electron and ion temperature. The effects of auroral particle precipitation on the density and temperature structures are studied for winter and low solar activity conditions. The high-latitude ionospheric features, such as the tongue of ionization, the main trough the polar ionization peak, the auroral ionization peak, the high-latitude ionization hole, the tongue of electron temperature, the high latitude minima of electron, temperature, and the ion temperature hot spot are obtained from calculations.Numerically obtained results are used for determining the HF propagation paths in the polar ionosphere. The effects of ionospheric irregularities on high frequency ionospheric radio waves are investigated by using a three-dimensional ray-tracing computer program. Ray-path trajectories are presented for different values of the elevation angle of transmission. From our study, it was found that large-scale irregularity structures of the high-latitude ionosphere, in the presence of the earth's magnetic field, significantly affect high frequency radio wave propagation.     

Far-field Gravity and Tilt Signals by Large Earthquakes: Real or Instrumental Effects?

A wide set of dynamics phenomena (i.e., Geodynamics, Post Glacial Rebound, seismicity and volcanic activity) can produce time gravity changes, which spectrum varies from short to long (more than 1 year) periods. The amplitude of the gravity variations is generally in the order of consequently their detection requires instruments with high sensitivity and stability: then, high quality experimental data. Spring and superconducting gravimeters are intensively used with this target and they are frequently jointed with tiltmeters recording stations in order to measure the elasto-gravitational perturbation of the Earth. The far-field effects produced by large earthquakes on records collected by spring gravimeters and tiltmeters are investigated here. Gravity and tilt records were analyzed on time windows spanning the occurrence of large worldwide earthquakes; the gravity records have been collected on two stations approximately 600 km distant. The background noise level at the stations was characterized, in each season, in order to detect a possible seasonal dependence and the presence of spectral components which could hide or mask other geophysical signals, such as, for instance, the highest mode of the Seismic Free Oscillation (SFO) of the Earth. Some spectral components (6.5; 8; 9; 14, 20, 51) have been detected in gravity and tilt records on the occasion of large earthquakes and the effect of the SFO has been hypothesized. A quite different spectral content of the EW and NS tiltmeter components has been detected and interpreted as a consequence of the radiation pattern of the disturbances due to the earthquakes. Through the analysis of the instrumental sensitivity, instrumental effects have been detected for gravity meters at very low frequency.F.S.E. (Fondo Sociale Europeo -European Community -)     

Plasma drift estimates from the Dynasonde: comparison with EISCAT measurements

Modern ionosondes make almost simultaneous measurements of the time rate of change of phase path in different directions and at different heights. By combining these Doppler measurements and angles of arrival of many such radar echoes it is possible to derive reliable estimates of plasma drift velocity for a defined scattering volume. Results from both multifrequency and kinesonde-mode soundings at 3-min resolution show that the Dynasonde-derived F-region drift velocity is in good agreement with EISCAT, despite data loss during intervals of blanketing by intense E-region ionisation. It is clear that the Tromsø Dynasonde, employing standard operating modes, gives a reliable indication of overall convection patterns during quiet to moderately active conditions.     

Variation of continental mantle heat flow with age. Possibility of discriminating between thermal models of the lithosphere

Continental mantle heat flow values are obtained by subtracting the radiogenic heat produced in the lower crust and lithosphere beneath the crust from reduced heat flow values reported for various heat flow provinces. The significance of continental mantle heat flow values thus obtained is that they can be considered essentially as representing the residual heat of cooling of the continental lithosphere. A plot of these mantle heat flow values against 1/t where t is the geologic age of the last thermal event suggests a linear trend. It is also found that the recently proposed relationQ=500 (1/t) for the variation of oceanic heat flowQ (in mW/M²) with aget (in million years) provides a reasonably good fit to the mantle heat flow data. The constant thickness plate model however, is found to be unsatisfactory in explaining the variation of continental mantle heat flow with age.     

Sealing law of far-field spectra based on observed parameters of the specific barrier model

A set of acceleration source spectra is constructed using the observed parameters of the specific barrier model of Papageorgiou and Aki. The spectra show a significant departure from the ²-model at the high frequency range. Specifically, the high frequency spectral amplitudes of seismic excitation are higher as compared to the level predicted by the ²-model. This is also supported by other observational evidence. The high frequency amplitudes of acceleration scale proportionally to the square root of the rupture areaS, to the rupture spreading velocityv, and to the local strain drop (/) (=strain drop in between barriers). The local strain drop in between barriers is not related in a simple way to the global strain drop, which is the strain drop estimated by assuming that it is uniform over the entire rupture area. Consequently, the similarity law does not apply. Using the source spectra which we constructed, we derive expressions for high frequency amplitudes of acceleration such asa _rms anda _max. Close to the fault both are independent of fault dimensions and scale as (/µ)(f)^1/2, while away from the fault plane they scale asW ^1/2(/µ)(f)^1/2, whereW is the width of the fault and f is the effective bandwidth of the spectra.     

Generators in atmospheric electricity

Summary The definition of a generator in atmospheric electricity is considered, and various phenomena are discussed as to whether they can be described as generators.     

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 role of wavegroup velocity in the blocking problem

Considering the blocking problem as a baroclinic instability problem in a dispersive wave system with diabatic heating effects, it is of great interest to investigate the role of wavegroup velocityv _gr in blocking processes, becausev _gr controls the energy transfer in the wave field. Using a Newtonian Cooling —type of forcing with a phase differencek to the main field and taking the linearized version of a two-level model, the phase speedc _r, the group velocityv _gr and the growth ratekc _i have been obtained as analytical functions of the mean zonal windU, the thermal windU _T, the coefficient of diabatic heating x, the phase differencek and the wavelengthL. Now the hypothesis is introduced, that a blocking should be expected, ifv _gr has a maximum value in the vicinity ofL _o, for whichc _r vanishes and thee-folding timet=1/kc _i (kc _i>0) is smaller than 6 days (see condition (20) in the text). One finds, that for special parameter combinations (U _T, U, ), where 15 m/secU _T25m/sec,U=10m/sec, 0.8·10^–51.5·10^–5 [sec^–1], certain valuesL _o with an appropriate phase differencek exist, which satisfy the conditions mentioned above (for values see Table 2). TherebyL _o varies within the range 8500 km <L _o<11000 km corresponding to the preferred planetary blocking wavenumber 2 in middle latitudes 50°<<70° N.     

Comparison of short- and long-periodP-wave travel times and arrival angles

Summary Using about 100 earthquakes recorded on both the short-period and long-period matched three-component seismograph systems at Umeå (UME) in northern Sweden, relativeP-wave travel times (LP-SP) are measured and arrival angles (incidence and azimuth) are determined from relative amplitudes ofP. Long-period onsets appear to register consistently later than short-period ones by a couple of tenths of a second. This figure is of the same order as the uncertainty in T _LS (difference between LP and SP onset times), but still is shown to be of some statistical significance. The effect is believed due, at least in part, to the relative inability of the long-period seismometer to visibly register an initial shorter-period lower-amplitude focal phase from some earthquakes.The dip on crustal interfaces implied by the short-period azimuth and incidence-angle anomalies is 7 to 15 times greater than that implied by the long-period anomalies. This effect is believed due to undulation or irregularity of, foremost, the Moho, as well as other crustal boundaries. The wavelength of, say, Moho undulations would seem to be limited above and below by considerations of the dominant long-period wavelengths (100 km) and of the diameter of Moho sampled by the short-period waves (50 km). The possible undesirable effect of the crustal layering on the components of surface motion, and thus on the arrival angles, is discussed and cautiously dismissed, due to the absence of sedimentary layering, in favour of the above interpretation.     

Review of Ray Theory Applications in Modelling and Imaging of Seismic Data

Throughout the last twenty years, 3D seismic ray modelling has developed from a research tool to a more operational tool that has gained growing interest in the petroleum industry. Various areas of application have been established and new ones are under development. Many of these applications require a modelling system with flexible, robust and efficient modelling algorithms in the core. The present paper reviews the basic elements of such a system, based on the open model concept and the wavefront construction technique. In the latter, ervený's dynamic ray tracing is an intrinsic part. The modelling system can be used for generating ray attributes and synthetic seismograms for realistic 3D surveys with tens of thousands of shots and receivers. Moreover, some other types of application areas are illustrated: Production of Green's functions for prestack depth migration and hybrid modelling (combined ray and finite-difference modelling), attribute mapping and illumination analysis, both for survey planning and interpretation. Finally, the concepts of isochron rays and velocity rays related to seismic isochrons have been introduced recently, with very interesting future applications.     

The Time-Varying Gravitational Potential Field of a Massive,Deformable Body

The internal as well as the external gravitational field of the Earth is computed under the assumption that (i) the mass distribution of the deformable body as well as (ii) the heights of the topographic surface and the depths of the interfaces (dynamic isostasy) vary in time. In order to represent those shape variations properly, the topographic masses as well as the interface masses are regarded as condensed following a proposal by F.R. Helmert (1884, p. 149-163). Accordingly time-varying simple layer mass densities are generated. Basic results are collected in 4 boxes.     

Low Free-Field Accelerations of the 1999 Kocaeli Earthquake?

The August 17, 1999 Kocaeli earthquake ruptured the earths surface along 145 km and produced a magnitude of M_w=7.4. As expected for such an event Modified Mercalli intensities of typically IX and X in the vicinity of the fault were determined. Yet the observed accelerations at the five near-fault sites remained amazingly small with horizontal PGA values of 0.14 g to 0.4 g. We attempt to resolve the enigma by modeling surface ground motion with a finite-difference algorithm, utilizing two different rupture and slip histories derived from the strong-motion observations and translate the computed horizontal motion in intensity values. We can show that (a) in a given simple crustal velocity model different slip distributions result in significantly different ground motion distributions in the vicinity of the fault even though both slip distributions fit the observed accelerometer data quite well. (b) Both slip distributions project high ground motion into areas adjacent to the fault where no accelerometer was located. (c) Both slip distributions are not fully compatible with observed intensity observations around the fault, although this could be partly attributed to the lack of knowledge regarding to the crustal structure. In the light of our results it would thus be misleading if the few strong-motion observations around the Kocaeli earthquake fault were taken as typical or representative for the entire area and for potential future events.     

Some exact solutions for the flow of fluid through tubes with uniformly porous walls

Summary This paper considers an incompressible fluid flowing through a straight, circular tube whose walls are uniformly porous. The flow is steady and one dimensional. The loss of fluid through the wall is proportional to the mean static pressure in the tube. Several formulations of the wall shear stress are considered; these formulations were motivated by the results from Hamel's radial flow problem, boundary layer flows/and boundary layer suction profiles. For each of these formulations exact solutions for the mean axial velocity and the mean static pressure of the fluid are obtained. Sample results are plotted on graphs. For the constant wall shear stress problem, the theoretical solutions compare favorably with some experimental results.Notations A, B, D, E constant parameters - a, b constant parameters - Ai(z), Bi(z) Airy functions - Ai, Bi derivatives of Airy functions - k constant of proportionality betweenV andp - L length of pores - p,p mean static pressure - p ₀ static pressure outside the tube - p ₀ value ofp atx=0 - Q constant exponent - R inside radius of the tube - T wall shear stress - T ₀ shear parameter - t wall thickness - U free stream velocity - ,u mean axial velocity - u ₀ value ofu atx=0 - V,V mean seepage velocity through the wall - v ₀ mean seepage velocity - x,x axial distance along the tube - z transformed axial distance - z ₀ value ofz atx=0 - mean outflow angle through the wall - cos - density of the fluid - wall shear stress - dynamic viscosity of the fluid - over-bar dimensional terms - no bar nondimensional terms The National Center for Atmospheric Research is sponsored by the National Science Foundation     

Comparison of electric and elastic parameters of indian precambrian rocks

Summary Examination of the relation between electric resistivity, compressional wave velocity and density of Indian Precambrian rocks is presented.

---

¶rt;m auumu ¶rt; muu nmuu, m n¶rt; uu u nmm ¶rt;uu n¶rt; u ¶rt;uu.