An interpretation of the Nusselt-Rayleigh number relationship for convection in a spherical shell

Recent studies on the relationship between the Nusselt ( Nu ) and Rayleigh ( Ra ) numbers for base-heated convection in a spherical shell with a constant viscosity show that the power-law index is around 1/4, which is different from the value of 1/3 predicted by a simple boundary layer theory. We show that such a difference may be caused by the flow pattern due to the geometry. The flow pattern of the convection in a spherical shell at relatively low Ra , at least, less than 10⁶, is characterized by narrow upwelling and broad downwelling, which is similar to the opposite flow pattern of internally heated convection. Convection in the internally heated case predicts the power-law index of 1/4. We demonstrate this relationship based on the concept of 'local' Rayleigh ( Ra1 ) and Nusselt ( Nu1 ) numbers     

A boundary layer model for mantle convection with surface plates

Summary. A simple, analytical model for mantle convection with mobile surface plates is presented. Our aim is to determine under what conditions free convection can account for the observed plate motions, and to evaluate the thermal structure of the mantle existing under these conditions. Boundary layer methods are used to represent two-dimensional cellular convection at large Rayleigh and infinite Prandtl numbers. The steady-state structure consists of cells with isentropic interiors enclosed by thermal boundary layers. Lithospheric plates are represented as upper surfaces on each cell free to move at a uniform speed. Buoyancy forces are concentrated in narrow rising and decending thermal plumes; torques imparted by these plumes drive both the deformable mantle and overlying plate. Solutions are found for a comprehensive range of cell sizes. We derive an expression for the plate speed as a function of its length, the mantle viscosity and surface heat flux. Using mean values for these parameters, we find that thermal convection extending to 700 km depth can move plates at 1 cm yr^-1, while convection through the whole mantle can move plates at 4–5 cm yr^-1. Analysis of the steady-state temperature field, for the case of heating from below, shows that the upper thermal boundary layer develops a complex structure, including an 'asthenosphere' defined by a local maximum in the geotherm occurring at depths of 50–150 km.     

Three-dimensional Rayleigh hysteresis of oriented core samples from the German Continental Deep Drilling Program: susceptibility tensor, Rayleigh tensor, three-dimensional Rayleigh law

Rayleigh hysteresis, as defined by the well-known Rayleigh relations, has been observed not only when magnetization of pyrrhotite-bearing KTB-samples is measured in parallel to a weak dc magnetic field, but also in experiments where field and measuring directions have been adjusted strictly perpendicularly to each other. Nine-tupels of independent Rayleigh hysteresis loops could thus be compiled. Their characteristic coefficients X ^ijk of initial susceptibility together with the Rayleigh loss coefficients α^jk have been proved to determine completely the samples' weak-field magnetic anisotropy. Interpreting the coefficient matrices ( X ^ijk) and (α^jk) as the tensor of initial susceptibility and the Rayleigh tensor, respectively, generalization of the isotropic Rayleigh relations in terms of corresponding tensor relationships has been suggested for the anisotropic case. Application to the KTB samples showed 3-D Rayleigh hysteresis measurements to be an excellent tool for rock magnetic analysis in terms of ore content and crystalline texture. In particular, a magnetocrystalline double texture of the basal planes of pyrrhotite precipitates and their [1120] directions of easy magnetization have been clearly detected. Surprisingly, the welt-known theorem α= const. X ²_I, formulated by Néel (1942) for the isotropic case, has been found to hold true even in tensor generalization (α_jk) = const ( X ²_jk). To reach sufficient resolution for the measurements performed, a sensitive vibrating coil magnetometer (VCM) has been developed.     

Three-dimensional VP and VP /VS models of the upper crust in the Friuli area (northeastern Italy)

3-D images of P velocity and P - to S -velocity ratio have been produced for the upper crust of the Friuli area (northeastern Italy) using local earthquake tomography. The data consist of 2565 P and 930 S arrival times of high quality. The best-fitting V _P and V _P / V _S 1-D models were computed before the 3-D inversion. V _P was measured on two rock samples representative of the investigated upper layers of the Friuli crust. The tomographic V _P model was used for modelling the gravity anomalies, by converting the velocity values into densities along three vertical cross-sections. The computed gravity anomalies were optimized with respect to the observed gravity anomalies. The crust investigated is characterized by sharp lateral and deep V _P and V _P / V _S anomalies that are associated with the complex geological structure. High V _P / V _S values are associated with highly fractured zones related to the main faulting pattern. The relocated seismicity is generally associated with sharp variations in the V _P / V _S anomalies. The V _P images show a high-velocity body below 6 km depth in the central part of the Friuli area, marked also by strong V _P / V _S heterogeneities, and this is interpreted as a tectonic wedge. Comparison with the distribution of earthquakes supports the hypothesis that the tectonic wedge controls most of the seismicity and can be considered to be the main seismogenic zone in the Friuli area.     

Estimation of uncertainties in eigenspectral estimates from decaying geophysical time series

Summary. The response of many dynamical systems to an impulse is a linear combination of decaying cosines. The frequencies of the cosines have generally been estimated in geophysics by periodogram analysis and little formal indication of uncertainty has been provided. This work presents an estimation procedure by the methods of complex demodulation and nonlinear regression that specifically incorporates in the basic model the decaying aspect of the cosines (periodogram analysis does not). The use of plots of the instantaneous phase as a function of time is shown to greatly enhance resolution. Expressions for the variances of eigenfrequencies, amplitudes, phases and damping constants Q are derived by non-linear least-squares. The results are illustrated, for the problem of the free oscillations of the Earth, by computations with the record made at Trieste of the Chilean earthquake of 1960 May 22. Sample values are periods and standard errors of 737.79 ± 0.13 s, 506.25 ± 0.13 s and 429.60 ± 0.14 s for ₀ T ₈, ₀ T ₁₃ and ₀ T ₁₆ with Q values and standard errors of 200 ± 14, 230 ± 28 and 215 ± 30, respectively.     

The detection of single-domain greigite (Fe3S4) using rotational remanent magnetization (RRM) and the effective gyro field (Bg): mineral magnetic and palaeomagnetic applications

The intensity of rotational remanent magnetization (RRM) acquired by single-domain greigite at a rotation frequency of 5 rps was combined with measurements of anhysteretic remanent magnetization (ARM) to calculate the effective biasing field ( Bg ) that produced the RRM. Samples of single-domain greigite had Bg values between -137 and -84 μT, and a MDF_RRM of c. 80 mT. By contrast, a suite of natural and synthetic ferrimagnetic iron oxide samples, including single-domain magnetite and y Fe₃O₄ tape particles, acquired Bg values between -3 and -14 μT, and MDF_RRM ranged between 43 and 68 mT (when RRM was acquired). Multidomain magnetite did not acquire a RRM. Bg values at 5 rps were calculated from previously published data for magnetite particles of different grain sizes, which revealed a minimum Bg value of -24 μT and a MDF_RRM of 57 mT for the finest fraction (0.2-0.8 μm in diameter). In a geological example, measurements of Bg and MDF_RRM were used to detect the presence of greigite in a 4 m long Late Weichselian sediment core. Variations in inclination, declination and the intensity of the natural remanent magnetization (NRM) correlate with changes in magnetic mineralogy.     

Waveform inversion of mantle Love waves:the Born seismogram approach

Summary. Normal mode theory, extended to the slightly laterally heterogeneous earth by the first-order Born approximation, is applied to the waveform inversion of mantle Love wave (200–500 s) for the Earth's lateral heterogeneity at l = 2 and a spherically symmétric anelasticity ( Q _μ) structure. The data are from the Global Digital Seismograph Network (GDSN). The l =2 pattern is very similar to the results of other studies that used either different méthods, such as phase velocity measurements and multiplet location measurements, or a different data set, such as mantle Rayleigh waves from different instruments. The results are carefully analysed for variance reduction and are most naturally explained by heterogeneity in the upper 420 km. Because of the poor resolution of the data set for the deep interior, however, a fairly large heterogeneity in the transition zones, of the order of up to 3.5 per cent in shear wave velocity, is allowed. It is noteworthy that Love waves of this period range cannot constrain the structure below 420 km and thus any model presented by similar studies below this depth are likely to be constrained by Rayleigh waves (spheroidal modes) only.
The calculated modal Q values for the obtained Q _μ model fall within the error bars of the observations. The result demonstrates the discrepancy of Rayleigh wave Q and Love wave Q and indicates that care must be taken when both Rayleigh and Love wave data, including amplitude information, are inverted simultaneously.
Anomalous amplitude inversions of G2 and G3, for example, are observed for some source-receiver pairs. This is due to multipathing effects. One example near the epicentral region, which is modelled by the obtained l = 2 heterogeneity, is shown.     

Seismic amplitude tomography for crustal attenuation beneath China

Amplitude tomography reconstructs seismic attenuation directly from recorded wave amplitudes. We have applied the tomography to amplitude data reported in the 'Annual Bulletin of Chinese Earthquakes' and interpreted the regionally varying crustal attenuation in terms of tectonics. The seismic amplitudes were originally recorded for determining the M _L and M _S magnitudes. They generally correspond to the maximum amplitudes of the horizontal components of the short-period S waves and intermediate-period Rayleigh waves. Both sets of measurements are sensitive to crustal structure. The peak amplitudes from M _L amplitudes spread spherically with significant dispersion and scattering. M _S amplitudes show cylindrical spreading with little dispersion. Average crustal Q values for attenuation at 1 Hz are 737 and 505 for M _L and M _S, respectively, with substantial regional variations. Frequency dependence in the attenuation is also indicated. Regions with the lowest attenuation (high Q values) are beneath the south China Block, Sichuan Basin, Ordos Platform, the Daxinganling and the Korea Craton. These tend to be tectonically inactive regions, which are generally dominated by intrusive and cratonic rocks in the upper crust. Regions with the highest attenuation (low Q values) are beneath Bohai Basin, Yunnan, eastern Songpan-Ganzi Terrain, margins of the Ordos platform and the Qilian Shan. These are predominantly active basins, grabens and fold belts. The continental margin also highly attenuates both S and surface waves.     

Numerical experiments on the onset of convective instability in the Earth's mantle

Summary. A simplified model of convection in the mantle is used to investigate the transient effect of cooling a fluid layer from above, The model, representing the mantle overlain by the lithosphere, consists of a two-dimensional fluid layer overlain by a solid conducting lid. The initial temperature of both layers is the same, with the top surface of the lid kept at 0°C throughout. We observe the onset of small-scale flow in the model. In the absence of internal heating the behaviour of the system is controlled by the Rayleigh number, R , and the ratio of the thicknesses of the two layers, a . The onset time of convection as defined by reference to conduction temperature profiles is related simply to a boundary layer critical Rayleigh number. The mean temperature profiles for the convection model are also compared with the observed depth—age relation for oceanic lithosphere and the results are used to estimate the viscosity of the mantle.     

Lateral variation and frequency dependence of coda-Q in the southern part of Iberia

b
A method based on the coda attenuation law: Q = Q ₀( f/f ₀)^v leads to the determination of the lateral variation of coda- Q in the southern part of the Iberian Peninsula using seismograms belonging to the seismological network of the Cartuja Observatory, located in Granada. The lateral variation of Q ₀ ( Q value corresponding to a reference frequency f ₀ of 1 Hz) and its frequency dependence for the 1 to 5 Hz frequency range are, in general, in agreement with coda- Q values for frequencies less than about 1 Hz, previously determined in the region under study.
To determine the coda- Q values analytical functions have been used to fit the magnification curves of the vertical component short-period seismographs belonging to the Cartuja network. The problem is solved by using least-squares techniques and non-linear inversion. The determined coda- Q ₀ values and its frequency dependence correlate well with several known geophysical parameters in the southern part of the Iberian Peninsula.     

Multimode surface waveform tomography of the Pacific Ocean: a closer look at the lithospheric cooling signature

We present a regional surface waveform tomography of the Pacific upper mantle, obtained using an automated multimode surface waveform inversion technique on fundamental and higher mode Rayleigh waves, to constrain the   V_SV   structure down to ∼400 km depth. We have improved on previous implementations of this technique by robustly accounting for the effects of uncertainties in earthquake source parameters in the tomographic inversion. We have furthermore improved path coverage in the South Pacific region by including Rayleigh wave observations from the French Polynesian Pacific Lithosphere and Upper Mantle Experiment deployment. This improvement has led to imaging of vertical low-velocity structures associated with hotspots within the South Pacific Super-Swell region. We have produced an age-dependent average cross-section for the Pacific Ocean lithosphere and found that the increase in   V_SV   with age is broadly compatible with a half-space cooling model of oceanic lithosphere formation. We cannot confirm evidence for a Pacific-wide reheating event. Our synthetic tests show that detailed interpretation of average   V_SV   trends across the Pacific Ocean may be misleading unless lateral resolution and amplitude recovery are uniform across the region, a condition that is difficult to achieve in such a large oceanic basin with current seismic stations.     

Anelasticity of the Crust and Upper Mantle of South America From the Inversion of Observed Surface Wave Attenuation

Fundamental-mode Rayleigh and Love waves generated by several earthquakes situated along great-circle paths between pairs of seismograph stations have been analysed to obtain coefficients of attenuation, group velocities, phase velocities, and specific quality factors in the period range 18–80s in two regions of the South American continent. One set of paths crosses the shield region which lies on the eastern coast and another set traverses the mountainous region inland. the average attenuation coefficient values are clearly higher in the tectonically active western region throughout the entire period range than in the eastern or shield region.
Inversion of the attenuation data yielded shear wave internal friction ( Q ^-1_β) models as a function of depth in the crust and upper mantle in both regions. A low- Q zone below the lithosphere is prominent in both regions. the results show that substantial variations of Q _β occur in the two regions of South America. the Q_β values were found to be inversely related to the heat flow values or to the temperature.     

Mantle convection under simulated plates: effects of heating modes and ridge and trench migration, and implications for the core—mantle boundary, bathymetry, the geoid and Benioff zones

Summary. Numerical convection models are presented in which plates are simulated by imposing piecewise constant horizontal velocities on the upper boundary. A 4 × 1 box of constant viscosity fluid and two-dimensional (2-D) flow is assumed. Four heating modes are compared: the four combinations of internal or bottom heating and prescribed bottom temperature or heat flux. The case with internal heating and an isothermal base is relevant to lower mantle or whole mantle convection, and it yields a lower thermal boundary layer which is laterally variable and can be locally reversed, corresponding to heat flowing back into the core locally. When scaled to the whole mantle, the surface deflections and gravity and geoid perturbations calculated from the models are comparable to those observed at the Earth's surface. For models with migrating ridges and trenches, the flow structure lags well behind the changing surface 'plate'configurations. This may help to explain the poor correlation between the main geoid features and plate boundaries. Trench migration substantially affects the dip of the cool descending fluid because of induced horizontal shear in the vicinity of the trench. Such shear is small for whole mantle convection, but is large for upper mantle convection, and would probably result in the Tonga Benioff zone dipping to the SE, opposite to the observed dip, for the case of upper mantle convection.     

The traveltime perturbations for seismic body waves in factorized anisotropic inhomogeneous media

The traveltime perturbation equations for the quasi-compressional and the two quasi-shear waves propagating in a factorized anisotropic inhomogeneous (FAI) media are derived. The concept of FAI media simplifies considerably these equations. In the FAI medium, the density normalized elastic parameters a _ijkl ( X _i ) can be described by the relation a _ijkl ( X _i) = f ²( x _i ) A _ijkl, where A _ijkl are constants, independent of coordinates x _i and f ²( x _i) is a continuous smooth function of x _i . The types of anisotropy ( A _ijkl ) and inhomogeneity [ f ( x _i)] are not restricted. The traveltime perturbations of individual seismic body waves ( q ^P , qS 1 and qS 2) propagating in the FAI medium depend, of course, both on the structural pertubations [δ f ²( x _i)] and on the anisotropy perturbations (δ A _ijkl ), but both these effects are fully separated. The perturbation equations for the time delay between the two qS -waves propagating in the FAI medium are simplified even more. If the unperturbed (background) medium is isotropic, the perturbation of the time delay does not depend on the structural perturbations (δ f ²( x _i) at all. This striking result, valid of course only in the framework of first-order perturbation theory, will simplify considerably the interpretation of the time delay between the two split qS -waves in inhomogeneous anisotropic media. Numerical examples are presented.     

Upper crustal shear velocity models from higher mode Rayleigh wave dispersion in Scotland

Summary. Group velocities for first and second higher mode Rayleigh waves, in the frequency range 0.8–4.8 Hz, generated from a local earthquake of magnitude 3.7 M _L in western Scotland, are measured at stations along the 1974 LISPB line. These provide detailed information about the crustal structure west of the line. The data divide the region into seven apparently homogeneous provinces. Averaged higher mode velocity dispersion curves for each province are analysed simultaneously using a linearized inversion technique, yielding regionalized shear velocity profiles down to a depth of 17 km into the upper crust. Shear wave velocity is between 3.0 and 3.4 km s⁻¹ in the upper 2 km, with a slow increase to around 3.8 km s⁻¹. P -wave models computed using these results agree with profiles from the LISPB and LUST refraction experiments.     

Localization of gravity and topography: constraints on the tectonics and mantle dynamics of Venus

We develop a method for spatio-spectral localization of harmonic data on a sphere and use it to interpret recent high-resolution global estimates of the gravity and topography of Venus in the context of geodynamical models. Our approach applies equally to the simple spatial windowing of harmonic data and to variable-length-scale analyses, which are analogous to a wavelet transform in the Cartesian domain. Using the variable-length-scale approach, we calculate the localized RMS amplitudes of gravity and topography, as well as the spectral admittance between the two fields, as functions of position and wavelength. The observed admittances over 10 per cent of the surface of Venus (highland plateaus and tessera regions) are consistent with isostatic compensation of topography by variations in crustal thickness, while admittances over the remaining 90 per cent of the surface (rises, plains and lowlands) indicate that long-wavelength topography is dominantly the result of vertical convective tractions at the base of the lithosphere. The global average crustal thickness is less than 30 km, but can reach values as large as 40 km beneath tesserae and highland plateaus. We also note that an Earth-like radial viscosity structure cannot be rejected by the gravity and topography data and that, without a mechanical model of the lithosphere, admittance values cannot constrain the thickness of the thermal boundary layer of Venus. Modelling the lithosphere as a thin elastic plate indicates that at the time of formation of relief in highland plateaus and tesserae, the effective elastic plate thickness, T_e , was less than 20 km. Estimates of T_e at highland rises are consistently less than 30 km. Our inability to find regions with T_e > 30 km is inconsistent with predictions made by a class of catastrophic resurfacing models.     

Nature and origin of magnetic minerals within the Middle Jurassic shallow-water carbonate rocks of the Paris Basin, France: implications for magnetostratigraphic dating

The Middle and Upper Jurassic Bathonian-Oxfordian shallow-water carbonate rocks from the Paris Basin, France, consist mainly of oolitic and bioclastic limestones that are hydrocarbon reservoirs in the subsurface. Despite a preliminary positive study, these deposits have been considered to be largely remagnetized (Rochette, private communication), and hence not amenable to palaeomagnetic dating. To establish their magnetic mineralogy and test this remagnetization hypothesis, we have used an integrated investigation combining petrographic, geochemical, rock-magnetic and palaeomagnetic measurements on samples extracted from five cores from the Paris Basin and from outcrops in Burgundy. Magnetic minerals in the Bathonian-Oxfordian carbonates include: (1) primary biogenic single-domain magnetite and detrital multi-domain Ti-magnetite and their oxidized form, maghemite; (2) authigenic spheres of magnetite probably related to hydrocarbons; and (3) goethite, either restricted to ferruginous ooid layers or resulting from surficial alteration, notably replacement of pyrite framboids. Rock-magnetic experiments carried out on 68 samples reveal H _cr/ H _c and M _rs/ M _s ratios ranging from 1.88 to 5.58 and 0.017 to 0.314, respectively. These values are clearly distinct from diagnostic values for a chemical remagnetization. Pyrrhotite was not identified within these sediments. Moreover, the average H _cr/ H _c ratio of 3.14 is significantly different from the value of 1.333 for natural pyrrhotite (Dekkers 1988). These results have a direct implication for the preservation of the primary magnetization; consequently, these deposits are selectively amenable for magnetostratigraphic dating and possible regional correlations.     

Rayleigh-wave amplitudes from earthquakes in the range 0–150°

Summary . Vertical component Rayleigh-wave amplitudes from 1461 shallow earthquakes recorded in the distance range 0–150° are analysed to separate the effects of earthquake size, epicentral distance (Δ) and recording station.
The estimated decay of amplitude with distance has the form of a theoretical curve for the decay of Rayleigh waves with distance if the assumption is made that the decay due to dispersion for the data analysed is that of an Airy phase. Writing the decay due to anelastic attenuation as exp (- k Δ), k is estimated to be 0.676/rad over the whole range of distance. If the distance effects are represented by a straight line of the form h log Δ+ constant, h is estimated to be 1.15. The calibration function for computing M _s derived from the estimated distance effects is very similar to that of Marshall & Basham.
Station effects on Rayleigh-wave amplitudes though statistically significant are small, and can probably be ignored in the computation of M _s.
Comparing the estimated surface-wave magnitudes (earthquake size) obtained in this study with the long and short period body-wave magnitudes ( m ^LP_b and m ^SP_b respectively) obtained by Booth, Marshall & Young for the same earthquake shows that m ^LP_b is about equal to M _s over the magnitude range of interest (˜4.0–7.0). The m ^LP_b and M_s relationship shows that the greater the long-period energy radiated by an earthquake the smaller proportionately is the short-period energy.     

Strain accumulation along an oblique plate boundary: the Reykjanes Peninsula, southwest Iceland

We use annual GPS observations on the Reykjanes Peninsula (RP) from 2000 to 2006 to generate maps of surface velocities and strain rates across the active plate boundary. We find that the surface deformation on the RP is consistent with oblique plate boundary motion on a regional scale, although considerable temporal and spatial strain rate variations are observed within the plate boundary zone. A small, but consistent increase in eastward velocity is observed at several stations on the southern part of the peninsula, compared to the 1993–1998 time period. The 2000–2006 velocities can be modelled by approximating the plate boundary as a series of vertical dislocations with left-lateral motion and opening. For the RP plate boundary we estimate left-lateral motion  18⁺⁴₋₃ mm yr⁻¹  and opening of  7⁺³₋₂ mm yr⁻¹  below a locking depth of  7⁺¹₋₂ km  . The resulting deep motion of  20⁺⁴₋₃ mm yr⁻¹  in the direction of  N(100⁺⁸₋₆)°E  agrees well with the predicted relative North America–Eurasia rate. We calculate the areal and shear strain rates using velocities from two periods: 1993–1998 and 2000–2006. The deep motion along the plate boundary results in left-lateral shear strain rates, which are perturbed by shallow deformation due to the 1994–1998 inflation and elevated seismicity in the Hengill–Hrómundartindur volcanic system, geothermal fluid extraction at the Svartsengi power plant, and possibly earthquake activity on the central part of the peninsula.