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.     

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.     

Differentiation of the sensitivity to copper and cadmium in different life stages of a copepod

Copper is more toxic to all life stages of the copepod Tisbe than cadmium. The most sensitive life stage of Tisbe to both copper and cadmium is the one-day-old nauplius. The resistance of larval stages of Tisbe increases with age (one-day-old nauplii 48h $\text{lc}{}_{50}\text{=0.3142}$ mg Cu l^?1. and 0.5384 mg Cd l^?1, 0.3415 mg Cu l^?1. and 0.645 mg Cd l^?1. for five-days-old nauplii and 0.5289 mg Cu l^?1. and 0.9061 mg Cd l^?1. for ten-days-old nauplii. The two reproductive stages of Tisbe tested (females with ovigerous bands and females bearing the first ovigerous sac) demonstrated an increased sensitivity to metals and proved more sensitive than the ten-days-old copepodids (only females with ovigerous bands had a similar sensitivity to copper with the ten-days-old copepodids).     

The thermal boundary-layer interpretation of D″ and its role as a plume source

The “anomalous” layer in the lowermost mantle, identified as D″ in the notation of K.E. Bullen, appears in the PREM Earth model as a 150 km-thick zone in which the gradient of incompressibility with pressure, $\text{d}\text{K}\text{d}\text{P}$, is almost 1.6, instead of 3.2 as in the overlying mantle. Since PREM shows no accompanying change in density or density gradient, we identify D″ as a thermal boundary layer and not as a chemically distinct zone. The anomaly in $\text{d}\text{K}\text{d}\text{P}$ is related to the temperature gradient by the temperature dependence of K_s, for which we present a thermodynamic identity in terms of accessible quantities. This gives the numerical result (?K_s/?T)_P=?1.6×10⁷ Pa K^?1 for D″ material. The corresponding temperature increment over the D″ range is 840 K. Such a layer cannot be a static feature, but must be maintained by a downward motion of the lower mantle toward the core-mantle boundary with a strong horizontal flow near the base of D″. Assuming a core heat flux of 1.6 × 10¹² W, the downward speed is 0.07 mm y^?1 and the temperature profile in D″, scaled to match PREM data, is approximately exponential with a scale height of 73 km. The inferred thermal conductivity is 1.2 W m^?1 K^?1. Using these values we develop a new analytical model of D″ which is dynamically and thermally consistent. In this model, the lower-mantle material is heated and softened as it moves down into D″ where the strong temperature dependence of viscosity concentrates the horizontal flow in a layer ～ 12 km thick and similarly ensures its removal via narrow plumes.     

The Rn content of groundwater in strata of the Plio-Pleistocene Osaka Group in the Sennan area, Japan

Studies on the ²²²Rn content of groundwater were conducted to obtain knowledge of its distribution in the geosphere just below ground surface and below the depth of 100 m in the Osaka Group (Plio-Pleistocene). Samples analyzed from ten bore holes within 10 m depth below land surface showed that a sandy layer contained groundwater with consistent radon concentrations except for a few locations where the variations in the water level were large due to precipitation and/or the water level was within 2 m of ground surface. The average ²²²Rn concentration in groundwater within 10 m of land surface was found to be approximately equal to $\text{1}\text{3}$ to $\text{3}\text{5}$ of the highest concentration found in groundwater from a well penetrating 100–200 m (420 pCi/l^?1).     

How strong is the magnetic field in the Earth's liquid core?

A crucial step in the investigation of the energetics of motions in the Earth's core and the generation of the geomagnetic field by the hydromagnetic dynamo process is the estimation of the average strength $\text{B}$ of the magnetic field B = B_p + B_T in the core. Owing to the probability that the toroidal field B_T in the core, which has no radial component, is a good deal stronger than the poloidal field B_p, direct downward extrapolation of the surface field to the core-mantle interface gives no more than an extreme lower limit to $\text{B}$. This paper outlines the indirect methods by which $\text{B}$ can be estimated, arguing that $\text{B}$ is probably about 10^?2 T (100 Γ) but might be as low as 10^?3 T (10 Γ) or as high as 5 × 10^?2 T (500 Γ).     

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.     

Grain boundary migration in olivine at atmospheric pressure

Annealing experiments in order to study grain boundary migration (GBM) were carried out at temperatures of 1513–1773 K from 10 min to 100 hours at atmospheric pressure. Grain growth due to GBM is observed in the formation of margins of neoblastic grains which display very different structures of dislocations from that of consumed porphyroclastic and initial neoblastic grains.The velocity of GBM obtained here is approximated to be $\text{c=}\text{gh}$ in which k is 1.15×10^?9 cm³ s^?1, Q is the activation energy for GBM in olivine at 210±20 kJ mol^?1 and ρ is the dislocation density of consumed olivine (cm^?2.Inasmuch as GBM in static annealing reduces stored strain energy in olivines, it is one of the softening processes counteracting work-hardening by dislocation multiplication and tangling as well as dislocation annihilation. GBM-softening is dominant in low temperature annealing but in high temperatures dislocation recovery predominantly takes place.     

Wide-band feedback seismometers

The effects of feedback in seismometers are summarised and the responsivities and detectivities of various open-and closed-loop instruments are compared and discussed. It is concluded that the most satisfactory system is displacement force-feedback employing a capacitive displacement transducer.Details are given of a Willmore Mk III(a) short-period seismometer which has been converted to a wide-band force-feedback instrument by the addition of a capacitive transducer. It provides a wide-band response flat to acceleration from d.c. to 25 Hz, with a theoretical detection level of 10^?10 m s^?2$\text{Hz}{}^{\mathrm{<mtext>?</mtext><mtext>1</mtext><mtext>2</mtext>}}$ and a dynamic range of 140 dB. Wide-band velocity, long-period narrow-band or short-period responses are easily obtained by filters outside the loop. The instrument has been compared with standard seismometers over a period of several months, and found to compare satisfactorily when enclosed in a small pressure vessel.     

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.     

Acute toxicity of copper to sea catfish

Sea catfish (Arius felis) were exposed to aqueous solutions of reagent grade cupric chloride in artificial seawater ($\text{30.0±2.0‰}$, 21–23°C) in four static bioassays. The 24, 48, 72 and 96 h LC₅₀ were calculated and found to be 5.43, 4.17, 3.57 and 2.40 mg l.^?1 copper, respectively. Experimental concentrations of copper producing subtle behavioral changes in this species correspond to less than 0.3% of the 72 h LC₅₀. Based on this comparison with literature values, a new, maximum ‘safe’ concentration for copper in marine waters of 0.01 mg l.^?1 is proposed.     

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.     

Transmission electron microscope observation of α, β and γ (Mg,Fe)2SiO4 in shocked meteorites: planar defects and polymorphic transitions

High pressure polymorphs of olivine (Mg, Fe)₂SiO₄ have been observed by transmission electron microscopy in the shocked Tenham and Catherwood meteorites. Planar defects are characterized in β and γ polymorphs and their possible role in polymorphic transitions is assessed and discussed in relation with the published literature. The stacking fault $\text{1}\text{2}\text{[}\text{1}\text{?}\text{01](010)}$ observed in the β-phase can produce a shear transformation to the γ-phase if it occurs on every other (010) plane. Conversely, the stacking fault $\text{1}\text{2}\text{[1}\text{1}\text{?}\text{2](110)}$ observed in the γ-phase can be transformed to the β-phase if it is repeated on every other (110) plane and to an intermediate, yet unobserved, ε1-phase if it is repeated on every (110) plane. Shear transformation between olivine and γ-spinel could involve the ε1-phase as an intermediate stage.     

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.     

Some models of the geomagnetic field in western Europe from 1960 to 1980

Achache et al.'s averages over 10 west European observatories of the annual mean magnetic components X, Y, Z are studied from 1953 to 1979. A 1970 jerk is fitted to the data by least squares, and the fit is quantitatively compared with that of a polynomial having the same number of free parameters as the jerk (a quintic). A crude correction for the 11-y sunspot cycle is also attempted. The jerk and quintic are equally good fits to X, the jerk is clearly better for Y, and the quintic is modestly better for Z. The jerk amplitudes of X and Z but not Y depend heavily on whether a sunspot correction is made, and such a correction may be required to estimate the energy of the internal part of the jerk in the manner of Malin and Hodder. The sunspot correction to Y probably cannot be brought above the noise because it is small and highly correlated with both the jerk and the quintic. The sunspot corrections to X and Z are statistically significant but depend on whether the core signal is taken as a jerk or a quintic. Courtillot and Le Mouël's physical model for the jerk predicts that the jerk amplitudes of X, Y, Z will be proportional to the longitudinal derivatives of X, Y and Z. This prediction is roughly verified for the magnitudes but not the signs; it is just possible that the signs in X and Z are not statistically significant.     

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.     

Velocity-density systematics,mean atomic weight,and the interior of the moon

Mean atomic weight profiles for the lunar mantle have been calculated from velocity-density systematic relations using lunar density and seismic velocity models. Despite large variability among the models, the calculation including Poisson's ratio yields a range of mean atomic weight values between 22 and 23 g mol^?1 below 150 km. A similar calculation for the Earth's mantle produces a mean atomic weight of 21.1 ±0.4 g mol^?1. This suggests that the Moon cannot be derived directly from the Earth's mantle, or that it has had a differentiation history different from the Earth's. The lunar $\text{m}\text{'}\text{s}$ require an Fe mole fraction between 0.25 and 0.33 for a pure olivine mantle, or between 0.33 and 0.45 for pure pyroxene.The present profiles are 0.5–3.0 g mol^?1 higher than those calculated from lunar compositional models based on lunar rock compositions and petrology and assumed lunar histories, indicating inadequacies in either the seismic or compositional models, or in both. The mean atomic weight approach provides a method of comparing the consistency of seismic and compositional models of planetary interiors.     

Shear modulus and Q of forsterite and dunite near partial melting from forced-oscillation experiments

An apparatus designed to determine the complex shear modulus of rock samples by forced torsion oscillations at high temperature and in the seismic frequency band 0.003–30 Hz is briefly described. Measurements were performed on natural dunite from Åheim (Norway) up to 1400°C and on polycrystalline forsterite up to 1500°C at 1 atm pressure. The two materials were chosen to study, by comparison, the effect of melt on the elasticity and anelasticity of mantle rocks.Between 1000 and 1200°C the absolute values of the shear modulus G are almost equal for both materials. Above 1200°C G for natural dunite decreases progressively with temperature and at 1400°C and 1 Hz reaches $\text{～}\text{1}\text{3}$ of its value at 1100°C. In contrast, G of pure forsterite depends little on temperature. For petrological reasons, supported by simultaneous measurements of the electric resistivity, there is strong evidence that the decrease of G in dunite above 1200°C is due to melt from the lower melting components of the dunite. Based on different models estimates of the melt fraction are made.At high temperature, in both materials Q^?1 is characterized by a monotonic decrease with frequency according to ω^?α, with α ≈ 0.25. An apparent activation energy of 38±5 kcal mol^?1 for forsterite and 48±8 kcal mol^?1 for dunite was found with no significant change in the regime of partial melting. From this it is concluded that Q^?1, even at partial melting, is dominated by solid state high temperature background absorption. There is no indication from these experiments for a constant-Q-band at low seismic frequencies or an increase of Q proportional to frequency as suggested by some seismologists. The present results are in good qualitative agreement with those for Young's modulus obtained previously by strain retardation experiments.     

Creep response of the lunar crust in mare regions from an analysis of crater deformation

The settling trends of 318 lunar mare craters are compared with predictions of numerical finite-element models in order to determine the creep response of the upper lunar mare crust. No settling is evident in craters smaller than 5 km in diameter. Settling rates of larger craters increase as function of crater size in a manner suggesting a non-linear lunar creep response corresponding to the power law $\text{ε}\text{?}\text{= 8.3 · 10}{}^{\mathrm{?34}}\text{σ}{}^2$ where ε&#x0301; is the strain rate and σ is the differential stress. However, the observed nonlinearity is probably an apparent nonlinearity resulting from the temperature induced viscosity decrease with depth due to a lunar crustal temperature gradient of 3° C/km and a creep activation energy of 20 kcal/mole. It is concluded that creep in the lunar medium is essentially Newtonian, and that the effective viscosity of the upper lunar mare crust is (1.6 ± 0.3) · 10²⁵ poise.     

The effect of solar activity on the secular variation of the geomagnetic field in Romania

The 11-year solar cycle effect in the geomagnetic components H and Z is made clear for Surlari Observatory and 19 repeat stations for the interval 1952–1974. The correlation with Wolf number and its time derivative is discussed in terms of the effects of the external and induced current systems.The $\text{H}\text{?}$ data available for solar cycle 20 (1964–1976) were processed to give the geographical distribution of the secular variation impulse for epoch 1969.5 in Romania. It is suggested that this distribution might reflect the deep internal structure of the area considered.A qualitative correlation is noted between long-period solar activity and variation of the horizontal component of the geomagnetic field at some repeat stations.