The effect of anisotropic diffusivity on rotating magnetoconvection in the Earth's core

A thermal diffusive process in the Earth's core is principally enhanced by small-scale flows that are highly anisotropic because of the Earth's rapid rotation and a strong magnetic field. This means that a thermal eddy diffusivity should not be a scalar but a tensor. The effect of such anisotropic tensor diffusivity, which is to be prescribed, on dynamics in the Earth's core is investigated through numerical simulations of magnetoconvection in a rapidly rotating system. A certain degree of anisotropy has an insignificant effect on the character, like kinetic and magnetic energies, of magnetoconvection in a small region with periodic boundaries in the three directions. However, in a region with top and bottom rigid boundary surfaces, kinetic and magnetic energies of magnetoconvection can be altered by the same degree of anisotropy. This implies that anisotropic tensor diffusivity affects on dynamics in the core, in particular near the boundary surfaces.     

The effect of anisotropic heat transport in the earth's core on the geodynamo

Intermediate dynamos are axisymmetric, spherical models that evade Cowling's theorem by invoking an α-effect to create the meridional magnetic field from the zonal. Usually the energy source maintaining the motions is a specified thermal wind, but here the dynamo is driven by the buoyancy created by a uniform distribution of heat sources. It has been argued by Braginsky and Meytlis (this journal, vol. 55, 1990) that, in a rapidly rotating, strongly magnetic system such as the Earth's core, heat is transported principally by a small-scale turbulence that is highly anisotropic. They conclude that the diffusion of heat parallel to the rotation axis is then significantly greater than it is in directions away from that axis. A preliminary study of the consequences of this idea is reported here. Solutions are derived numerically using both isotropic and non-isotropic thermal diffusivity tensors, and the results are compared. It is shown that even a small degree of anisotropy can materially alter the character of the dynamo.     

气煤弹性各向异性系数实验测试

气煤裂隙丰富,具有弹性各向异性,确定气煤弹性各向异性类型与大小,对煤田各向异性研究具有重要意义. 本文在对气煤煤样孔隙率测试的基础上,测量了顺煤层面、垂直于煤层面、与煤层面成45°三个方向上纵波速度（V_P）和横波速度（V_SH、V_SV）,计算弹性常数. 对于具有垂向对称轴横向各向同性模型,各向异性系数平均值|ε|≤0.2、|δ|≤0.2和|γ|≤0.2,表明气煤是弱各向异性介质;对于横向对称轴横向各向异性模型,各向异性系数平均值ε(V)≈δ(V),按照Thomsen观点气煤各向异性是由椭圆裂隙引起弹性各向异性.本文还讨论了各向异性系数与孔隙率的关系,表明各向异性的大小与孔隙率之间有一定的关系. 因此当煤层厚度与其顶底板岩性不变或变化很小,可以认为气煤各向异性主要由裂隙引起的,气煤各向异性大小和方向能够反映煤的裂缝密度和方位.     

Anisotropic attenuation of stratified viscoelastic media

An important cause of seismic anisotropic attenuation is the interbedding of thin viscoelastic layers. However, much less attention has been devoted to layer‐induced anisotropic attenuation. Here, we derive a group of unified weighted average forms for effective attenuation from a binary isotropic, transversely isotropic‐ and orthorhombic‐layered medium in the zero‐frequency limit by using the Backus averaging/upscaling method and analyse the influence of interval parameters on effective attenuation. Besides the corresponding interval attenuation and the real part of stiffness, the contrast in the real part of the complex stiffness is also a key factor influencing effective attenuation. A simple linear approximation can be obtained to calculate effective attenuation if the contrast in the real part of stiffness is very small. In a viscoelastic medium, attenuation anisotropy and velocity anisotropy may have different orientations of symmetry planes, and the symmetry class of the former is not lower than that of the latter. We define a group of more general attenuation‐anisotropy parameters to characterize not only the anisotropic attenuation with different symmetry classes from the anisotropic velocity but also the elastic case. Numerical tests reveal the influence of interval attenuation anisotropy, interval velocity anisotropy and the contrast in the real part of stiffness on effective attenuation anisotropy. Types of effective attenuation anisotropy for interval orthorhombic attenuation and interval transversely isotropic attenuation with a vertical symmetry (vertical transversely isotropic attenuation) are controlled only by the interval attenuation anisotropy. A type of effective attenuation anisotropy for interval TI attenuation with a horizontal symmetry (horizontal transversely isotropic attenuation) is controlled by the interval attenuation anisotropy and the contrast in the real part of stiffness. The type of effective attenuation anisotropy for interval isotropic attenuation is controlled by all three factors. The magnitude of effective attenuation anisotropy is positively correlated with the contrast in the real part of the stiffness. Effective attenuation even in isotropic layers with identical isotropic attenuation is anisotropic if the contrast in the real part of stiffness is non‐zero. In addition, if the contrast in the real part of stiffness is very small, a simple linear approximation also can be performed to calculate effective attenuation‐anisotropy parameters for interval anisotropic attenuation.     

Reference ellipsoids for anisotropic media

Perturbation methods are common tools for describing wave propagation in weakly anisotropic media. The anisotropic medium is replaced by an average isotropic medium where wave propagation can be treated analytically and the correction for the effect of anisotropy is computed by perturbation techniques. This works well for anisotropies of up to 10%. Some materials (e.g. shales), however, can exhibit a much stronger anisotropy. In this case a background is required which still can be treated analytically but is applicable to stronger P-wave anisotropy. We present an averaging technique to compute a best-fitting ellipsoidal medium to an arbitrary anisotropic medium. Ellipsoidal media are sufficiently simple for analytical expressions to be available for many applications and allow consideration of strong P-wave anisotropy. The averaging of the arbitrary anisotropic medium can be carried out globally (i.e. for the whole sphere) or sectorially (e.g. for seismic waves propagating predominantly in the vertical direction). We derive linear relationships for the coefficients of the ellipsoid which depend on the elastic coefficients of the anisotropic medium. We also provide specifications for best-fitting elliptical and best-fitting isotropic media. Numerical examples for different rocks demonstrate the improved approximation of the anisotropic model obtained using the formulae derived, compared with the conventionally used average isotropic medium.     

Local turbulence in the Earth's core

Abstract

It is shown that in the Earth's core, where the geodynamo is at work (and is supplied with energy by the prevailing unstable density stratification), a buoyancy instability of a local character exists which is highly supercritical. This instability results in fully developed turbulence dominated by small scale vortices. The influence of the Earth's rotation and of the magnetic field produced by the geodynamo makes this small scale turbulence highly anisotropic. A qualitative picture of this local anisotropic turbulence is devised and the main parameters characterizing it are estimated. Expressions for the turbulent diffusivity are developed and discussed.     

A dynamo model with double diffusive convection for Mercury's core

A recent dynamo model for Mercury assumes that the upper part of the planet's fluid core is thermally stably stratified because the temperature gradient at the core–mantle boundary is subadiabatic. Vigorous convection driven by a superadiabatic temperature gradient at the boundary of a growing solid inner core and by the associated release of light constituents takes place in a deep sub-layer and powers a dynamo. These models have been successful at explaining the observed weak global magnetic field at Mercury's surface. They have been based on the concept of codensity, which combines thermal and compositional sources of buoyancy into a single variable by assuming the same diffusivity for both components. Actual diffusivities in planetary cores differ by a large factor. To overcome the limitation of the codensity model, we solve two separate transport equations with different diffusivities in a double diffusive dynamo model for Mercury. When temperature and composition contribute comparable amounts to the buoyancy force, we find significant differences to the codensity model. In the double diffusive case convection penetrates the upper layer with a net stable density stratification in the form of finger convection. Compared to the codensity model, this enhances the poloidal magnetic field in the nominally stable layer and outside the core, where it becomes too strong compared to observation. Intense azimuthal flow in the stable layer generates a strong axisymmetric toroidal field. We find in double diffusive models a surface magnetic field of the observed strength when compositional buoyancy plays an inferior role for driving the dynamo, which is the case when the sulphur concentration in Mercury's core is only a fraction of a percent.     

Anisotropic turbulent thermal diffusion and thermal convection in a rapidly rotating fluid sphere

Estimates of the molecular values of magnetic, viscous and thermal diffusion suggest that the state of the Earth’s core is turbulent and that complete numerical simulation of the geodynamo is not realizable at present. Large eddy simulation of the geodynamo with modelling of the sub-grid scale turbulence must be used. Current geodynamo models effectively model the sub-grid scale turbulence with isotropic diffusivities larger than the molecular values appropriate for the core. In the Braginsky and Meytlis (1990) picture of core turbulence the thermal and viscous diffusivities are enhanced up to the molecular magnetic diffusivity in the directions of the rotation axis and mean magnetic field. We neglect the mean magnetic field herein to isolate the effects of anisotropic thermal diffusion, enhanced or diminished along the rotation axis, and explore the instability of a steady conductive basic state with zero mean flow in the Boussinesq approximation. This state is found to be more stable (less stable) as the thermal diffusion parallel to the rotation axis is increased (decreased), if the transverse thermal diffusion is fixed. To examine the effect of simultaneously varying the diffusion along and transverse to the rotation axis, the Frobenius norm is used to control for the total thermal diffusion. When the Frobenius norm of the thermal diffusion tensor is fixed, it is found that increasing the thermal diffusion parallel to the rotation axis is destabilising. This result suggests that, for a fixed total thermal diffusion, geodynamo codes with anisotropic thermal diffusion may operate at lower modified Rayleigh numbers.     

A multi‐azimuth seismic refraction study for a horizontal transverse isotropic medium: physical modelling results

To investigate the characteristics of the anisotropic stratum, a multi‐azimuth seismic refraction technique is proposed in this study since the travel time anomaly of the refraction wave induced by this anisotropic stratum will be large for a far offset receiver. To simplify the problem, a two‐layer (isotropy–horizontal transverse isotropy) model is considered. A new travel time equation of the refracted P‐wave propagation in this two‐layer model is derived, which is the function of the phase and group velocities of the horizontal transverse isotropic stratum. In addition, the measured refraction wave velocity in the physical model experiment is the group velocity. The isotropic intercept time equation of a refraction wave can be directly used to estimate the thickness of the top (isotropic) layer of the two‐layer model because the contrast between the phase and group velocities of the horizontal transverse isotropic medium is seldom greater than 10% in the Earth. If the contrast between the phase and group velocities of an anisotropic medium is small, the approximated travel time equation of a refraction wave is obtained. This equation is only dependent on the group velocity of the horizontal transverse isotropic stratum. The elastic constants A₁₁, A₁₃, and A₃₃ and the Thomsen anisotropic parameter ε of the horizontal transverse isotropic stratum can be estimated using this multi‐azimuth seismic refraction technique. Furthermore, under a condition of weak anisotropy, the Thomsen anisotropic parameter δ of the horizontal transverse isotropic stratum can be estimated by this technique as well.     

Seismic anisotropy analysis at the Low-Noise Underground Laboratory (LSBB) of Rustrel (France)

Seismic anisotropy of a fractured karstic limestone massif in sub-parallel underground galleries is studied. As the fractures are mostly vertically oriented, the seismic properties of the massif are approximated by horizontal transverse isotropy (HTI). Several data inversion methods were applied to a seismic dataset of arrival-times of P and S-waves.The applied methods include: isotropic tomography, simple cosine function fit, homogeneous Monte-Carlo anisotropic inversion for the parameters of horizontal transverse isotropy and anisotropic tomography for tilted transversely isotropic bodies. All methods lead to the conclusion that there is indeed an anisotropy present in the rock massif and confirm the direction of maximum velocity parallel to the direction of fracturing. Strong anisotropy of about 15% is found in the studied area. Repeated measurements show variations of the P-wave parameters, but not of the S-wave parameters, which is reflecting a change in water saturation.     

Microseismic Full Waveform Modeling in Anisotropic Media with Moment Tensor Implementation

Seismic anisotropy which is common in shale and fractured rocks will cause travel-time and amplitude discrepancy in different propagation directions. For microseismic monitoring which is often implemented in shale or fractured rocks, seismic anisotropy needs to be carefully accounted for in source location and mechanism determination. We have developed an efficient finite-difference full waveform modeling tool with an arbitrary moment tensor source. The modeling tool is suitable for simulating wave propagation in anisotropic media for microseismic monitoring. As both dislocation and non-double-couple source are often observed in microseismic monitoring, an arbitrary moment tensor source is implemented in our forward modeling tool. The increments of shear stress are equally distributed on the staggered grid to implement an accurate and symmetric moment tensor source. Our modeling tool provides an efficient way to obtain the Green’s function in anisotropic media, which is the key of anisotropic moment tensor inversion and source mechanism characterization in microseismic monitoring. In our research, wavefields in anisotropic media have been carefully simulated and analyzed in both surface array and downhole array. The variation characteristics of travel-time and amplitude of direct P- and S-wave in vertical transverse isotropic media and horizontal transverse isotropic media are distinct, thus providing a feasible way to distinguish and identify the anisotropic type of the subsurface. Analyzing the travel-times and amplitudes of the microseismic data is a feasible way to estimate the orientation and density of the induced cracks in hydraulic fracturing. Our anisotropic modeling tool can be used to generate and analyze microseismic full wavefield with full moment tensor source in anisotropic media, which can help promote the anisotropic interpretation and inversion of field data.     

One-Dimensional qP-Wave Velocity Model of the Upper Crust for the West Bohemia/Vogtland Earthquake Swarm Region

The western part of the Bohemian Massif (West Bohemia/Vogtland region) is characteristic in the relatively frequent recurrence of intraplate earthquake swarms and in other manifestations of past-to-recent geodynamic activity. In this study we derived 1D anisotropic qP-wave model of the upper crust in the seismogenic West Bohemia/Vogtland region by means of joint inversion of two independent data sets - travel times from controlled shots and arrival times from local earthquakes extracted from the WEBNET seismograms. We derived also simple 1-D P-wave and S-wave isotropic models. Reasons for deriving these models were: (a) only simplified crustal velocity models, homogeneous half-space or 1D isotropic layered models of this region, have been derived up to now and (b) a significant effective anisotropy of the upper crust in the region which was indicated recently by S-wave splitting. Both our anisotropic qP-wave and isotropic P-and S-wave velocity models are constrained by four layers with the constant velocity gradient. Weak anisotropy for P-waves is assumed. The isotropic model is represented by 9 parameters and the anisotropic one is represented by 24 parameters. A new robust and effective optimization algorithm - isometric algorithm - was used for the joint inversion. A two-step inversion algorithm was used. During the first step the isotropic P- and S-wave velocity model was derived. In the second step, it was used as a background model and the parameters of anisotropy were sought. Our 1D models are adequate for the upper crust in the West Bohemia/Vogtland swarm region up to a depth of 15 km. The qP-wave velocity model shows 5% anisotropy, the minimum velocity in the horizontal direction corresponds to an azimuth of 170°. The isotropic model indicates the V_P/V_S ratio variation with depth. The difference between the hypocentre locations based on the derived isotropic and anisotropic models was found to be several hundreds of meters.     

Bayesian Markov chain Monte Carlo inversion for anisotropy of PP- and PS-wave in weakly anisotropic and heterogeneous media

A single set of vertically aligned cracks embedded in a purely isotropic background may be considered as a long-wavelength effective transversely isotropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly anisotropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Markov chain Monte Carlo inversion method of PP- and PS-wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inversion scheme taking only the PP-wave into account.     

Anisotropic velocity analysis1

Results from walkaway VSP and shale laboratory experiments show that shale anisotropy can be significantly anelliptic. Heterogeneity and anellipticity both lead to non-hyperbolic moveout curves and the resulting ambiguity in velocity analysis is investigated for the case of a factorizable anisotropic medium with a linear dependence of velocity on depth. More information can be obtained if there are several reflectors. The method of Dellinger et al. for anisotropic velocity analysis in layered transversely isotropic media is examined and is shown to be restricted to media having relatively small anellipticity. A new scheme, based on an expansion of the inverse-squared group velocity in spherical harmonics, is presented. This scheme can be used for larger anellipticity, and is applicable for horizontal layers having monoclinic symmetry with the symmetry plane parallel to the layers. The method is applied to invert the results of anisotropic ray tracing on a model Sand/shale sequence. For transversely isotropic media with small anisotropy, the scheme reduces to the method of Byun et al. and Byun and Corrigan. The expansion in spherical harmonics allows the P-phase slowness surface of each layer to be determined in analytic form from the layer parameters obtained by inversion without the need to assume that the anisotropy is weak.     

Comparison of DC sensitivity patterns for anisotropic and isotropic media

Many rocks possess electrical properties with a clearly expressed anisotropy. The anisotropic character of the rocks is often overlooked in forming the Fréchet derivatives or sensitivity functions for parameter updating during the inversion of DC resistivity data. In this study we have compared the sensitivity patterns for an isotropic, homogeneous model with that for a transversely isotropic (i.e. anisotropic) model having a tilted axis of symmetry using a pole–pole array. The sensitivity functions are expressed in terms of the derivatives of the electric potential U with respect to the average conductivity σ_m (geometric mean of the longitudinal and transverse conductivities) and the coefficient of anisotropy λ. Results are plotted in both cross-section form and plan view for various dip and strike angles of the axis of symmetry. The derivative dU/dλ decreases more rapidly than the isotropic value dU/dσ, and shows pronounced asymmetry and weakening of magnitude with increasing dip of the plane of symmetry. The derivative dU/dσ_m also exhibits the asymmetric pattern (except for vertical and horizontal dip cases). The positive region between the electrodes only extends to a small depth compared to the isotropic derivative, even in the case of a vertical axis of symmetry (VTI medium). The ratio of this anisotropic derivative to the isotropic derivative, when plotted as a function of position and depth shows prominent differences in both the sign and the magnitude of the sensitivities, especially for steep dips and for strongly anisotropic rocks. The plot highlights the dangers of an isotropic assumption. Even for mildly anisotropic rocks (λ < 1.2) the possibility for error in interpretation is considerable. Combined borehole and surface measurements are needed to diagnose anisotropy. Further work is needed to design optimal electrode configurations in anisotropic situations.     

Mean-field concept and direct numerical simulations of rotating magnetoconvection and the geodynamo

Mean-field theory describes magnetohydrodynamic processes leading to large-scale magnetic fields in various cosmic objects. In this study magnetoconvection and dynamo processes in a rotating spherical shell are considered. Mean fields are defined by azimuthal averaging. In the framework of mean-field theory, the coefficients which determine the traditional representation of the mean electromotive force, including derivatives of the mean magnetic field up to the first order, are crucial for analyzing and simulating dynamo action. Two methods are developed to extract mean-field coefficients from direct numerical simulations of the mentioned processes. While the first method does not use intrinsic approximations, the second one is based on the second-order correlation approximation. There is satisfying agreement of the results of both methods for sufficiently slow fluid motions. Both methods are applied to simulations of rotating magnetoconvection and a quasi-stationary geodynamo. The mean-field induction effects described by these coefficients, e.g., the α-effect, are highly anisotropic in both examples. An α²-mechanism is suggested along with a strong γ-effect operating outside the inner core tangent cylinder. The turbulent diffusivity exceeds the molecular one by at least one order of magnitude in the geodynamo example. With the aim to compare mean-field simulations with corresponding direct numerical simulations, a two-dimensional mean-field model involving all previously determined mean-field coefficients was constructed. Various tests with different sets of mean-field coefficients reveal their action and significance. In the magnetoconvection and geodynamo examples considered here, the match between direct numerical simulations and mean-field simulations is only satisfying if a large number of mean-field coefficients are involved. In the magnetoconvection example, the azimuthally averaged magnetic field resulting from the numerical simulation is in good agreement with its counterpart in the mean-field model. However, this match is not completely satisfactory in the geodynamo case anymore. Here the traditional representation of the mean electromotive force ignoring higher than first-order spatial derivatives of the mean magnetic field is no longer a good approximation.     

The effects of different parameter regimes in geodynamo simulations

Over the past 10 years, geodynamo simulations have grown rapidly in sophistication. However, it is still necessary to make certain approximations in order to maintain numerical stability. In addition, models are forced to make assumptions about poorly known parameters for the Earth's core. Different magnetic Prandtl numbers have been used and different assumptions about the presence of radiogenic heating have been made. This study examines some of the consequences of different approximations and assumptions using the Glatzmaier–Roberts geodynamo model. Here, we show that the choice of magnetic Prandtl number has a greater influence on the character of the magnetic field produced than the addition of a plausible amount of radiogenic heating. In particular, we find that prescribing a magnetic Prandtl number of unity with Ekman number limited by current computing resources, results in magnetic fields with significantly smaller intensities and variabilities compared with the much more Earth-like results obtained from simulations with large magnetic Prandtl numbers. A magnetic Prandtl number of unity, with both the viscous and magnetic diffusivities set to the Earth's magnetic diffusivity, requires a rotation rate much smaller than that of the Earth for currently reachable Ekman numbers. This results in a reduced dominance of the Coriolis forces relative to the buoyancy forces, and therefore, a reduction in the magnetic field intensity and the variability compared to the large Prandtl number cases.     

各向异性介质中的AVO

分析了横向各向同性和方位各向异性介质的本构关系,由此讨论弹性波在两种各向异性介质中的传播特点,提出可表征这两种介质各向异性程度的广义参数．以此为基础讨论了两种各向异性介质中存在水平界面时的反射系数近似式,将Dely等人推导的横向各向同性介质中的反射系数公式推广到方位各向异性介质的主轴方向上．根据算例讨论修正的Banik和Thomsen的近似式,着重分析两种各向异性介质中的AVO关系及其对实际勘探的影响和指导意义．