High-temperature anelasticity and elasticity of mantle peridotite

An apparatus has been devised which allows precise creep and relaxation measurements to be made on minerals and rocks at temperatures up to 1600°C and at very low deviatoric stresses (1 < σ < 300 bar). This paper is concerned with measurements on mantle peridotite (lherzolite) from Balmuccia (Zone of Ivrea, Italy).The reaction of the sample to a step-like increase in stress is called its “creep function”. It is shown that the creep function contains all the necessary information to derive the spectra of the quality factor Q(ω) and of Young's modulus E(ω), within the seismic range of frequencies, provided the material behaves as a linear system. This has been proven up to a strain of 5 × 10^?5.The Q^?1-spectra at 1200 and 1300°C, obtained by Fourier inversion from the creep function, show no pronounced peak in the frequency band 0.01 < tf < 1 Hz and exhibit a general tendency to decrease slightly with frequency. The creep function: ?(t) = ?_u · [1 + 3.7 · q · {(1 + 50t)^0.27 ? }], where q is related to Q, satisfactorily describes the data at high temperatures and leads to $\text{Q}{}^{\mathrm{?1}}\text{(ω, T) = 3 × 10}{}^3\text{· ω}{}^{\mathrm{<mtext>?0.27\; ·\; </mtext><mtext>exp</mtext>}}\text{(}\text{?30}\text{RT}\text{)}$E(ω) is related to Q(ω) by the material dispersion equation. Above 1100°C the unrelaxed Young's modulus decreases rapidly with temperature according to an activation energy of about 20 kcal/mole. A lowering of short period S-wave velocity by 40% and P-wave velocity by 10% occurs below the solidus. Therefore, no partial melting is required in the asthenosphere.Steady-state creep at low axial stresses (20 < σ < 100 bar), obtained from the same rock, follows the relation $\text{?}\text{?}\text{= 3 × 10}{}^7\text{· δ}{}^{1.4}\text{·}\text{exp}\text{(}\text{?125}\text{RT}\text{)}$ indicative of grain boundary diffusion or superplasticity. At higher stresses a power law $\text{?}\text{?}\text{= 45 · δ}{}^4\text{·}\text{exp}\text{(}\text{?125}\text{RT}\text{)}$ typical of dislocation creep, is found.The frequency dependence of Q and the ratio of the activation energies of Q and are indicative of so called “high-temperature background absorption”, as the dominant mechanism, and of a diffusion-controlled dislocation mobility common to both absorption and creep. From a, b, and c, relations between the effective viscosity η_f and Q of the form: $\text{log}\text{η}{}_{\mathrm{e??}}\text{=}\text{1}\text{α}\text{·}\text{log}\text{Q ? (n ? 1) ·}\text{log}\text{ω +}\text{log}\text{D}$ are derived, where α ～ 0.25, n is the power of σ, and D is a constant.     

Upper mantle plasticity from laboratory experiments

On the basis of two assumptions i.e. (1) plastic and anelastic behaviour of the upper mantle can be approximated by the behaviour of the dominant mineral olivine, and (2) the behaviour of natural olivine and synthetic forsterite are similar, we have investigated the flow laws and the flow microstructures of forsterite single crystals. The results obtained between 1400–1650°C and 10–100 MPa suggest a model of climb controlled creep in which the a edge dislocations are dominant. The activation energy measured in that regime is 4.7 eV, close to that of Si self-diffusion and the flow law is $\text{?}\text{?}\text{=10}{}^6\text{σ}{}^{2.6}\text{exp}\text{(?4.7}\text{eV/k}\text{T)}$, where σ is in MPa. Extrapolation of these results to the upper mantle would imply very low stresses (i.e. ?10 MPa) in the asthenosphere. However the effect of pressure and grain size are unknown and extrapolation to very low stresses is not straightforward.     

High-pressure recovery of olivine: implications for creep mechanisms and creep activation volume

High-temperature and high-pressure recovery experiments were made on experimentally deformed olivines at temperatures of 1613–1788 K and pressures of 0.1 MPa to 2.0 GPa. In the high-pressure experiments, a piston cylinder apparatus was used with BN and NaCl powder as the pressure medium, and the hydrostatic condition of the pressure was checked by test runs with low dislocation density samples. No dislocation multiplication was observed. The kinetics of the dislocation annihilation process were examined by different initial dislocation density runs and shown to be of second order, i.e.

$\text{d}\text{ρ}\text{d}\text{t}\text{= ?p}{}^2\text{K}{}_0\text{exp}\text{[?(}\text{E}{}^1\text{+PV}{}^1\text{RT}\text{]}$

where ρ is the dislocation density, k₀ is a constant, $\text{E}{}^1\text{and}\text{V}{}^1$ are the activation energy and volume respectively, and P, R and T are pressure, gas constant and temperature, respectively. Activation energy and volume were estimated from the temperature and pressure dependence of the dislocation annihilation rate as $\text{E}{}^1\text{=389±59}\text{kJ mol}{}^{\mathrm{?1}}$ and $\text{V}{}^1\text{=14±2}\text{cm}{}^3\text{mol}{}^{\mathrm{?1}}$, respectively.The diffusion constants relevant to the dislocation annihilation process were estimated from a theoretical relation k=αD where $\text{k=k}{}_0\text{exp}\text{[?(E}{}^1\text{+ PV}{}^1\text{)/RT], D}$ is the diffusion constant and α is a non-dimensional constant of ca. 300. The results agree well with the self-diffusion constant of oxygen in olivine. This suggests that the dislocation annihilation is rate-controlled by the (oxygen) diffusion-controlled dislocation climb.The mechanisms of creep in olivine and dry dunite are examined by using the experimental data of static recovery. It is suggested that the creep of dry dunite is rate-controlled by recovery at cell walls or at grain boundaries which is rate-controlled by oxygen diffusion. Creep activation volume is estimated to be 16±3 cm³ mol^?1.     

Geothermal effects in the formation of electrically conducting zones and temperature distribution in the earth

Pollack and Chapman have shown that the surface heat flow in continental regions is dependent not only on the earth's crust below the observation site, but also on the upper mantle there. Therefore heat flow can be used to investigate the role of the thermal conditions in the creation of the electrically conductive zones in both the crust and mantle.Empirical exponential formulas describe the depth to the conductivity increase in the crust corresponding to granitization, the depth to the conductive zone at the top of the asthenosphere (SLVZ), as a function of heat flow. Comparing the latter with temperature estimations in the asthenosphere it is concluded that partial melting of the upper mantle occurs only where $\text{q ? 42}\text{m W m}{}^{\mathrm{?2}}\text{? 1}\text{HFU}$.The depth to the conductivity increase corresponding to the mineralogic phase transition in the upper mantle is increased with high temperatures. Such a conductive zone shows that the maximum temperature difference between stable platform areas and active zones is about 1000°C.     

Constraints on the origins of lunar magnetism from electron reflection measurements of surface magnetic fields

Apollo 15 and 16 subsatellite measurements of lunar surface magnetic fields by the electron reflection method are summarized. Patches of strong surface fields ranging from less than $\text{1}\text{4}$° to tens of degrees in size are found distributed over the lunar surface, but in general no obvious correlation is observed between field anomalies and surface geology. In lunar mare regions a positive statistical correlation is found between the surface field strength and the geologic age of the surface as determined from crater erosion studies. However, there is a lack of correlation of surface field with impact craters in the mare, implying that mare do not have a strong large-scale uniform magnetization as might be expected from an ancient lunar dynamo. This lack of correlation also indicates that mare impact processes do not generate strong magnetization coherent over ～ 10 km scale size. In the lunar highlands fields of >100 nT are found in a region of order 10 km wide and >300 km long centered on and paralleling the long linear rille, Rima Sirsalis. These fields imply that the rille has a strong magnetization (>5 × 10^?6 gauss cm³ gm^?1 associated with it, either in the form of intrusive, magnetized rock or as a gap in a uniformly magnetic layer of rock. However, a survey of seven lunar farside magnetic anomalies observed by the Apollo 16 subsatellite suggests a correlation with inner ejecta material from large impact basins. The implications of these results for the origin of lunar magnetism are discussed.     

Global electromagnetic induction in the Moon and planets

A summary of experiments and analyses concerning electromagnetic induction in the Moon and other extraterrestrial bodies is presented. Magnetic step-transient measurements made on the lunar dark side show the eddy current response to be the dominant induction mode of the Moon. Analysis of the poloidal field decay of the eddy currents has yielded a range of monotonic conductivity profiles for the lunar interior: the conductivity rises from 3·10^?4 mho/m at a depth of 170 km to 10^?2 mho/m at 1000 km depth. The static magnetization field induction has been measured and the whole-Moon relative magnetic permeability has been calculated to be $\text{μ}\text{μ}{}_0\text{= 1.01 ± 0.06}$. The remanent magnetic fields, measured at Apollo landing sites, range from 3 to 327 γ. Simultaneous magnetometer and solar wind spectrometer measurements show that the 38-γ remanent field at the Apollo 12 site is compressed to 54 γ by a solar wind pressure increase of 7·10^?8 dyn/cm². The solar wind confines the induced lunar poloidal field; the field is compressed to the surface on the lunar subsolar side and extends out into a cylindrical cavity on the lunar antisolar side. This solar wind confinement is modeled in the laboratory by a magnetic dipole enclosed in a superconducting lead cylinder; results show that the induced poloidal field geometry is modified in a manner similar to that measured on the Moon. Induction concepts developed for the Moon are extended to estimate the electromagnetic response of other bodies in the solar system.     

Fracture toughness and subcritical crack growth during high-temperature tensile deformation of Westerly granite and Black gabbro

The double torsion testing method has been used to determine catastrophic and subcritical crack propagation parameters for pre-cracked specimens of Westerly granite and Black gabbro under a number of environmental conditions.The critical stress intensity factor for catastrophic crack propagation (fracture toughness) of granite and gabbro has been measured at temperatures from 20 to 400°C, in a vacuum. At 20°C, the fracture toughness of Westerly granite was $\text{1.79 ± 0.02}\text{MPa}\text{·}\text{m}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, and for two blocks of Black gabbro it was $\text{3.03 ± 0.08}\text{MPa}\text{·}\text{m}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ and $\text{2.71 ± 0.15}\text{MPa}\text{·}\text{m}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, respectively. These values are very close to those reported by other investigators for tests conducted in air of ambient humidity at room temperature. For both rocks, fracture toughness at first increased slightly, and then decreased steadily on raising the temperature above ambient conditions. This behaviour is explained in terms of the density and distribution of thermally induced microcracks, as determined by quantitative optical microscopy.Subcritical crack growth behaviour has been studied at temperatures up to 300°C, and under water vapour at pressures of 0.6 to 15 kPa. Both the load relaxation and incremental constant displacement rate forms of the double torsion testing method were utilised to generate stress intensity factor/crack velocity diagrams. Crack growth was measured over the velocity range 5 × 10^?3 to 10^?7 m · s^?1. Increasing both temperature and water vapour pressure resulted in substantially higher crack growth rates. The overall effect of raising the temperature over the range studied here (20–300°C) was to increase the crack growth rate in granite and gabbro by ～5 and 7 orders of magnitude, respectively, at constant stress intensity factor and vapour pressure of water. For both rocks, the slopes of stress intensity factor/crack velocity curves were sensitive to changes in both temperature and water vapour pressure at low values of the latter parameter. Slopes fell substantially on raising the water vapour pressure, but were relatively insensitive to changes in temperature at these higher pressures. No subcritical crack growth limit was encountered.Estimates of the uncertainty in our experimental data are given. From the results of multiple load relaxation experiments on Westerly granite specimens, we estimate the uncertainty in position of stress intensity factor/crack velocity curves along the stress intensity axis to be c. 10% of the fracture toughness, and the uncertainty in slope of such curves to be c. 12%.Problems associated with the extrapolation of our experimental data to regions of higher effective confining pressure in the Earth's crust are discussed.     

The cooling Earth: A reappraisal

By treating the lithosphere as a diffusive boundary layer to mantle convection, the convective speed or mantle creep rate, $\text{?}\text{?}$, can be related to the mantle-derived heat flux, $\text{Q}\text{?}$. If cell size is independent of $\text{Q}\text{?}{}^2$ then $\text{?}\text{?}\text{∝}\text{Q}\text{?}$. (If cell size varies with $\text{Q}\text{?}$, then a different power law prevails, but the essential conclusions are unaffected.) Then the factthat for constant thermodynamic efficiency of mantle convection, the mechanical power dissipation is proportionalto $\text{Q}\text{?}$, gives convective stress $\text{σ ∝}\text{Q}\text{?}{}^{\mathrm{?1}}$, i.e. the stress increases as the convection slows. This means an increasing viscosityor stiffness of the mantle which can be identified with a cooling rate in terms of a temperature-dependent creep law. If we suppose that the mantle was at or close to its melting point within 1 or 2 × 10⁸ years of accretionof the Earth, the whole scale of cooling is fixed. The present rate of cooling is estimated to be about 4.6 × 10^?8 deg y^?1 for the average mantle temperature, assumed to be 2500 K, but this very slow cooling rate represents a loss ofresidual mantle heat of 7 × 10¹² W, about 30% of the total mantle-derived heat flux. This conclusion requires theEarth to be deficient in radioactive heat, relative to carbonaceous chondrites. A consideration of mantle outgassing and atmospheric argon leads to the conclusion that the deficiency is due to depletion of potassium, and that the K/U ratio of the mantle is only about 2500, much less than either the crustal or carbonaceous chondritic values. Thetotal terrestrial potassium is estimated to be about 6 × 10²⁰ kg. Acceptance of the cooling of the Earth removes the necessity for potassium in the core.     

An experimental study of magnetite-titanomagnetite interdiffusion

Interdiffusion experiments were performed between Fe₃O₄ (single crystal) and Fe_2.8Ti_0.2O₄ (powder), under self-buffering conditions (temperature range 600–1034°C), and for various oxygen potentials at 1400°C. Profiles of Fe and Ti were obtained by electronprobe microanalysis, and the interdiffusion coefficient $\text{D}$ was calculated by the Boltzmann-Matano method. Low-temperature data at 3 mole% Ti could be described by $\text{D}\text{= (3.85}{}_{\mathrm{?1.11}}{}^{\mathrm{+1.68}}\text{) × 10}{}^{\mathrm{?3}}\text{exp}\text{(2.23 ± 0.04}\text{eV}\text{/kT)}\text{cm}{}^2\text{/}\text{s}$. An estimate is given for the time to interdiffuse 2μm at various temperatures, and the results compared with recent experiments.     

Rupture process of the Miyagi-Oki,Japan, earthquake of June 12, 1978

The faulting mechanism and multiple rupture process of the M = 7.4 Miyagi-Oki earthquake are studied using surface and body wave data from local and worldwide stations. The main results are as follows. (1) P-wave first motion data and radiation patterns of long-period surface waves indicate a predominantly thrust mechanism with strike N10° E, dip 20°W, and slip angle 76°. The seismic moment is 3.1 × 10²⁷ dyne-cm. (2) Farfield SH waveforms and local seismograms suggest that the rupture occurred in two stages, being concordant with the two zones of aftershock activity revealed by the microearthquake network of Tohoku University. The upper and lower zones, located along the westward-dipping plate interface, are separated by a gap at a depth of 35 km and have dimensions of 37 × 34 and 24 × 34 km², respectively. Rupture initiated at the southern end of the upper aftershock zone and propagated at N20°W subparallel to the trench axis. About 11 s later, the second shock, which was located 30 km landward (westward) of the first, initiated at the upper corner of the lower aftershock zone and propagated down-dip N80°W. Using Haskell modelling for this rupture process, synthetic seismograms were computed for teleseismic SH waves and nearfield body waves. Other parameters determined are: seismic moment $\text{M}{}_0\text{= 1.7 × 10}{}^{27}\text{dyne-cm, slip dislocation}\text{u}\text{= 1.9}\text{m}\text{, Δσ = 95}\text{bar, rupture velocity}\text{ν = 3.2}\text{km s}{}^{\mathrm{?1}}\text{,}\text{rise time}\text{τ = 2}\text{s}$, for the first event; $\text{M}{}_0\text{= 1.4 × 10}{}^{27}\text{dyne-cm}\text{,}\text{u}\text{= 2.4}\text{m}\text{, Δσ = 145}\text{bar}$, for the second event; and time separation between the two shocks ΔT = 11 s. The above two-segment model does not explain well the sharp onsets of the body waves at near-source stations. An initial break of a small subsegment on the upper zone, which propagated down-dip, was hypothesized to explain the observed near-source seismograms. (3) The multiple rupture of the event and the absence of aftershocks between the two fault zones suggests that the frictional and/or sliding characteristics along the plate interface are not uniform. The rupture of the first event was arrested, presumably by a region of high fracture strength between the two zones. The fracture energy of the barrier was estimated to be 10¹⁰ erg cm^?2. (4) The possible occurrence of a large earthquake has been noted for the region adjacent to and seaward of the area that ruptured during the 1978 event. The 1978 event does not appear to reduce the likelihood of occurrence of this expected earthquake.     

Physical constraints for the analysis of the geomagnetic secular variation

Slow changes in the magnetic field are believed to originate in the core of the Earth. Interpretation of these changes requires knowledge both of the vertical component of the field and of its rate of change at the core-mantle boundary (CMB). While various spherical harmonic models show some agreement for the field at the CMB, those for secular variation (SV) do not. SV models depend heavily on annual means at relatively few and poorly distributed magnetic observatories. In this paper, the SV at the CMB is modelled by fitting 15-year differences in the annual means of the X, Y and Z components (from 1959 to 1974). The model is made unique by imposing the constraint that $\text{?}{}_{\mathrm{CMB}}\text{B}\text{?}{}_{\mathrm{r}}{}^2\text{d}\text{S}$ be a minimum, using the method of Shure et al. (1982). If SV is attributed to motions of core fluid, then this model will yield, in some sense, the slowest core motions. The null space is determined by the distribution of observations, and therefore, to be consistent, only those observatories have been retained which recorded almost continuously throughout the interval 1959–1974.The method allows misfit between the model and the observations. The best value for the misfit can be derived from estimates of errors in the data, or alternatively, because larger misfit leads to smoother models (i.e., smaller $\text{?}\text{B}\text{?}{}_{\mathrm{r}}{}^2\text{d}\text{S}$), the best value can be estimated subjectively from the final appearance of the model. Both procedures have their counterparts in the conventional spherical harmonic expansion approach, when smoothing is achieved by lowering the truncation level. The new proposal made in this paper is to use objective criteria for determining the misfit, based on the assumption that diffusion is negligible, in which event all integrals $\text{∫}\text{B}\text{?}{}_{\mathrm{r}}{}^2\text{d}\text{S}$ will vanish when S_i is a region on the CMB bounded by a contour of zero vertical component of field. For the 1965 definitive model which is adopted here, and for most other contemporary models, there are six such areas, giving five independent integrals (the integrals over the six regions must sum to zero if ? · B = 0). Tabulating these integrals for various choices of the misfit gives minimum values near 2 nT y^?1. It is impossible to achieve this good a fit to the data using a reasonable model derived by truncating the spherical harmonic expansion. The value 2 nT y^?1 corresponds to errors of ～ 20 nT in individual annual means, which is rather larger than expected from the scatter in the data.     

Splitting of dislocations in olivine,cross-slip-controlled creep and mantle rheology

The splitting of the [100] dislocations of forsterite Mg₂SiO₄ is investigated in a hard-sphere model. Glissile splittings exist in (001) and the most energetically favorable is the one that does not involve cutting of SiO₄ tetrahedra: $\text{[100] → [}\text{1}\text{6}\text{,}\text{1}\text{9}\text{,}\text{1}\text{6}\text{+ [}\text{2}\text{3}\text{, 0, 0] + [}\text{1}\text{6}\text{,}\text{1}\text{9}\text{,}\text{1}\text{6}$ [100] screw dislocations are shown to be able to split simultaneously on (001) and (010) according to the reaction: $\text{[100] = [}\text{1}\text{6}\text{,}\text{1}\text{36}\text{,}\text{1}\text{4}\text{] + [}\text{1}\text{6}\text{,}\text{1}\text{9}\text{,}\text{1}\text{6}\text{] + [}\text{1}\text{3}\text{, 0, θ] + [}\text{1}\text{6}\text{,}\text{1}\text{36}\text{,}\text{1}\text{4}\text{] + [}\text{1}\text{6}\text{,}\text{1}\text{6}\text{,}\text{1}\text{6}\text{]}$ This sessile configuration is analogous to the one found for screw dislocations of body-centered cubic metals. It accounts forthe long rectilinear sessile screw segments commonly observed, in experimentally and naturally deformed olivine. A new creep model for the upper mantle is proposed, where recovery is controlled by cross-slip of screw dislocations instead of climb of edge dislocations. The creep law, fitted on the experimental results of the literature is: $\text{?}\text{?}\text{= 1.2 · 10}{}^4\text{σ}{}^2\text{exp}\text{? [(125000 ? 11.7σ + PV}{}^1\text{/RT]>}$ (σ in bars) the activation volume for cross-slip is estimated and viscosity-depth curves are plotted. The proposed creep mechanism is grain-size independent and less non-Newtonian than the climb-controlled one; it is found to be dominant over the latter at stresses smaller than 100 bar.     

Effect of oxygen partial pressure on the dislocation recovery in olivine: a new constraint on creep mechanisms

The dislocation annihilation rate in experimentally deformed olivine single crystals was measured as a function of oxygen partial pressure (P_O2). It was shown that the dislocation annihilation rate decreased with increasing P_O2. This result is inconsistent with the reported P_O2 dependence of creep rate () in single olivine crystals, thus indicating that the creep in single olivine crystals is not rate-controlled by recovery, under the experimentally investigated conditions.     

Variation of the magnetic properties in a basaltic dyke with concentric cooling zones

The effects of the variation of magnetic grain size on the magnetic properties of rocks have been studied throughout a reversely magnetized basaltic dyke with concentric cooling zones.Except in a few tachylites in which the magnetic mineral is a Ti-rich titanomagnetite, in the bulk of the dyke the magnetization is carried by almost pure magnetite grains. Although the percentage p of these magnetic oxides varies slightly, the large changes in the various magnetic parameters observed across the dyke are essentially attributable to large variations in the grain size of the magnetic particles.From the outer scoria region, where the magnetic grains are a mixture of single-domain (SD) and superparamagnetic (SP) grains, to the tachylite zone with finely crystallized basaltic glass containing interacting elongated SD particles, one observes an increase of both the ratio of the saturation remanent magnetization and the saturation induced magnetization J_rs/J_is, the bulk coercive force H_c, the median destructive field MDF, the intensity of the remanent magnetization J_r, and the Koenigsberger ratio Q. In the tachylites these parameters reach unusually high values, for subaerial basalts:

$\text{〈}\text{J}{}_{\mathrm{rs}}\text{J}{}_{\mathrm{is}}\text{〉 = 0.3, 〈H}{}_{\mathrm{c}}\text{〉 = 460 O}{}_{\mathrm{e}}\text{, 〈MDF〉 = 620 O}{}_{\mathrm{e}}\text{r.m.s., 〈J}{}_{\mathrm{r}}\text{〉 = 2.7 · 10}{}^{\mathrm{?2}}\text{e.m.u. cm}{}^{\mathrm{?3}}\text{〈Q〉 = 24}$

These parameters decrease in the basalt toward the centre of the dyke where pseudo-single-domain (pseudo-SD) particles coexist together with multidomain (MD) grains. The susceptibility remains approximately constant from the inner basalt to the tachylite, but increases in the scoria up to values 10 times higher owing to the presence of SP particles. The magnetic viscosity increases also drastically toward the margin of the dyke due to an increase of the fraction of the SD particles just above the superparamagnetic threshold.     

Earthquake similarity laws

Source parameters of 27 major shallow earthquakes in the magnitude range 7.0–8.6, which occurred during 1906–1969, are used to establish dimensionless invariants involving the fault dimensions, average slip, and the rise time. It is found that these entities are expressible as simple functions of the subsonic shear Mach number (M) and the cube root of the seismic potency, $\text{(}\text{U}\text{S)}{}^{\mathrm{<mtext>1</mtext><mtext>3</mtext>}}$. Moreover, a new principle is suggested according to which all dimensionless numbers which can be constructed from the basic fault elements are simple powers of the contraction factor $\text{(1-M}{}^2\text{)}{}^{\mathrm{<mtext>?</mtext><mtext>1</mtext><mtext>2</mtext>}}$ with coefficients of the order unity. The laws of dynamical similarity thus found are those appropriate for subsonic rupture in which the Mach number is very close to unity and the radiation efficiency is between $\text{1}\text{6}$ and $\text{1}\text{3}$. The empirical similarity laws are shown to be compatible with a source model in which the fault plane is simulated by a flexible membrane with additional restoring stiffness forces provided by an elastic medium attached to it on one side. Results suggest the possibility that earthquake rupture, together with the radiation of seismic waves, terminates at the moment that Mach 1 is reached.     

General relationships among sound speeds: II. Theory and discussion

The dependence of bulk sound speed V_φ upon mean atomic weight $\text{m}$ and density ρ can be expressed in a single equation:

$\text{V}{}_{\mathrm{\phi }}\text{=Bρ}{}^{\mathrm{\lambda }}\text{(}\text{m}{}_0\text{m}{}^{\mathrm{<mtext>[</mtext><mtext>1</mtext><mtext>2</mtext><mtext>+\lambda (1?c)]</mtext>}}\text{(}\text{km/sec}\text{)}$

Here B is an empirically determined “universal” parameter equal to 1.42, $\text{m}{}_0\text{= 20.2}$, a reference mean atomic weight for which well-determined elastic properties exist, and λ = 1.25 is a semi empirical parameter equal to $\text{γ ?}\text{1}\text{3}$ where γ is a Grüneisen parameter. The constant c = (? ln V_M/? ln $\text{m}\text{)}{}_{\mathrm{X}}$, where V_M is molar volume, is in general different for different crystal structure series and different cation substitutions. However, it is possible to use c_Fe = 0.14 for Fe²⁺Mg²⁺ and GeSi substitutions and c_Ca ? 1.3 for CaMg substitutional series. With these values it is pos to deduce from the above equation Birch's law, its modifications introduced by Simmons to account for Ca-bearing minerals, variations in the seismic equation of state observed by D.L. Anderson, and the apparent proportionality of bulk modulus K to V_M^?4.     

Spatial power spectra of the crustal geomagnetic field and core geomagnetic field

Lowes (1966, 1974) has introduced the function R_n defined by $\text{R}{}^{\mathrm{n}}\text{=(n + 1)}\text{∑}\text{m=0}\text{∞}\text{[(g}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{+ (h}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{]}\text{where}\text{g}{}_{\mathrm{n}}{}^{\mathrm{m}}\text{and}\text{h}{}_{\mathrm{n}}{}^{\mathrm{m}}$ are the coefficients of a spherical harmonic expansion of the scalar potential of the geomagnetic field at the Earth's surface. The mean squared value of the magnetic field B = ??V on a sphere of radius r > α is given by $\text{〈}\text{B}\text{·〉 =}\text{∑}\text{n=1}\text{∞}\text{R}{}^{\mathrm{n}}\text{(}\text{a/}\text{r}\text{)}{}^{\mathrm{2n=4}}\text{where}\text{a}$ is the Earth's radius. We refer to R_n as the spherical harmonic spatial power spectrum of the geomagnetic field.In this paper it is shown that Rⁿ = R^M_n = R^C_n where the components R_n^M due to the main (or core) field and R_n^C due to the crustal field are given approximately by $\text{R}{}^{\mathrm{M}}{}_{\mathrm{n}}\text{= [}\text{(n =1)/}\text{(n + 2)}\text{](1.142 × 10}{}^9\text{)(0.288}{}^{\mathrm{n}}\text{Λ}{}^2\text{R}{}^{\mathrm{C}}{}_{\mathrm{n}}\text{= [}\text{(n =1){[1 — exp(-n/290)]/}\text{(n/290)} 0.52 Λ}{}^2\text{where}\text{Iγ = 1 nT}$. The two components are approximately equal for n = 15.Lowes has given equations for the core and crustal field spectra. His equation for the crustal field spectrum is significantly different from the one given here. The equation given in this paper is in better agreement with data obtained on the POGO spacecraft and with data for the crustal field given by Alldredge et al. (1963).The equations for the main and crustal geomagnetic field spectra are consistent with data for the core field given by Peddie and Fabiano (1976) and data for the crustal field given by Alldredge et al. The equations are based on a statistical model that makes use of the principle of equipartition of energy and predicts the shape of both the crustal and core spectra. The model also predicts the core radius accurately. The numerical values given by the equations are not strongly dependent on the model.Equations relating average great circle power spectra of the geomagnetic field components to R_n are derived. The three field components are in the radial direction, along the great circle track, and perpendicular to the first two. These equations can, in principle, be inverted to compute the R_n for celestial bodies from average great circle power spectra of the magnetic field components.     

Flow at high Reynolds numbers through anisotropic porous media

An approximate expression is developed for the relationship between the hydraulic gradient (J), the specific discharge (q) and fluid and porous matrix properties in the case of saturated, steady and uniform (macroscopic) flow of a Newtonian liquid at high Reynolds numbers through a homogeneous anisotropic porous medium:

$\text{g}\text{J}\text{=(v}\text{w}{}^{\mathrm{(2)}}\text{+}\text{B}{}^{\mathrm{(4)}}\text{:}\text{qq}\text{/q+}\text{C}{}^{\mathrm{(3)}}\text{·}\text{q}\text{)·}\text{q}$

In this expression, the tensors w⁽²⁾, B⁽⁴⁾ and C⁽³⁾ denote properties of the solid matrix only. The tensors W⁽²⁾, and C⁽³⁾ are symmetrical; the tensor B⁽⁴⁾ is symmetrical only in the first and last pairs of indices. It seems that no mathematical expression with a finite number of parameters exists, which can serve as a universal exact expression for the sought relationship between J and q.     

Source mechanism of the 1906 San Francisco earthquake

Surface-wave amplitudes in the period range 50–100 s at eight European and North American stations, horizontal slip profiles along the rupture zone and the timing of certain events along the fault during rupture time are all engaged in unison to reconstruct the motion at the source. A modified source model is used to accommodate a moving rupture with variable dislocation in the direction of propagation.It is inferred that the rupture started at about 13 h 11 m 55 s GMT near San Juan Bautista and propagated unilaterally northwestward along N35°W over 400 km with an average rupture velocity of 3.5 km/s. At 13 h 12 m 12 s, the dynamic shear front, moving with the rupture speed, hit the Lick Observatory. Then, at 13 h 12 m 18 s, the rupture arrived to the vicinity of the epicenter in the Santa Cruz Mountains given by B. Bolt. There the slip changed sharply from an average of 0.5 m to a high value of 3 m causing extensive landslides and avalanches. At 13 h 12 m 32.5 s two railroad clocks at San Rafael were stopped. Finally, at 13 h 12 m 36 s the offset front hit the Naval Observatory at Mare Island and stopped the astronomical clocks there. Conspicuous surface waves, visible on Wiechert seismograms in Europe in the period range 55–65 s, reflect the true rupture time.The seismic data inversion yields an effective radiation source some 240 km long with an average vertical extent of some 34 km over a total fault length of 400 km ($\text{U}\text{d}\text{S ? 29,000}\text{m km}{}^2\text{or}\text{μ}\text{U}\text{d}\text{S ? 9 · 10}{}^{27}\text{dyn cm}$). It began at the Santa Cruz Mountains and ended some 20 km off coast Point Arena. Thus, due to the nonuniform slip profile, only $\text{3}\text{5}$ of the total fracture length contributed to the far radiation field.Although the product of the average source displacement (over the entire fault) and the vertical extent appears to be fairly well determined from the surface-wave spectrums, the separate values of these entities cannot be uniquely determined. If the average surface displacements (～ 3.2 m) are diagnostic of the entire fault, a vertical extent of H = 34 km is required.Finally, a new analysis of surface waves from the Alaska earthquake of July 10, 1958, the Queen Charlotte Islands earthquake of August 22, 1949 and the Kern County shock of July 21, 1952, enables us to draw parallels between the three biggest major events which occurred along the NE Pacific coast during 1906–1958. A common feature of all of these earthquakes is that vertical failure extents of 30–40 km are implied.     

A confirmation of the correlation between hydrocarbons and chlorophyll a in the upper euphotic zone

The average hydrocarbon content found in 14 water samples from the euphotic zone off Northwest Africa was 4.4 μg l.^?1 with no extreme values. The average chlorophyll $\text{a}$ content was 1.9 μg l.^?1. The data fit a model proposed in a previous paper (Zsolnay, 1977). The resulting equation was $\text{HC}\text{= ?2.7 + 4.82}\text{Ch}\text{?0.942}\text{Ch}{}^2$ and could explain 59% of the variance in the hydrocarbon distribution. This indicates a correlation that is significant at the 0.002 level.