Secular variation in numerical geodynamo models with lateral variations of boundary heat flow

We study magnetic field variations in numerical models of the geodynamo, with convection driven by nonuniform heat flow imposed at the outer boundary. We concentrate on cases with a boundary heat flow pattern derived from seismic anomalies in the lower mantle. At a Rayleigh number of about 100 times critical with respect to the onset of convection, the magnetic field is dominated by the axial dipole component and has a similar spectral distribution as Earth’s historical magnetic field on the core-mantle boundary (CMB). The time scales of variation of the low-order Gauss coefficients in the model agree within a factor of two with observed values. We have determined the averaging time interval needed to delineate deviations from the axial dipole field caused by the boundary heterogeneity. An average over 2000 years (the archeomagnetic time scale) is barely sufficient to reveal the long-term nondipole field. The model shows reduced scatter in virtual geomagnetic pole positions (VGPs) in the central Pacific, consistent with the weak secular variation observed in the historical field. Longitudinal drift of magnetic field structures is episodic and differs between regions. Westward magnetic drift is most pronounced beneath the Atlantic in our model. Although frozen flux advection by the large-scale flow is generally insufficient to explain the magnetic drift rates, there are some exceptions. In particular, equatorial flux spot pairs produced by expulsion of toroidal magnetic field are rapidly advected westward in localized equatorial jets which we interpret as thermal winds.     

Inner-core conductivity in numerical dynamo simulations

Speculation about its possible super-rotation has drawn the attention of many geophysical researchers to the Earth’s inner core. An issue of special interest for geodynamo modelling is the influence of the inner-core conductivity. It has been suggested that the finite magnetic diffusivity of the inner core prevents more frequent reversals of the Earth’s magnetic field. We explore the possible influence of the inner-core conductivity by comparing convection-driven 3D dynamo simulations with insulating or conducting inner cores (CIC) at various parameters. The influence on the field structure in the outer core is only marginal. The time behaviour of dipole-dominated non-reversing dynamos is also little affected. Concerning reversing dynamos, the inner-core conductivity reduces the number of short dipole-polarity intervals with a typical length of a few thousand years. Reversals are always correlated with low dipole strength and these short intervals are found in periods where the dipole moment stays low. Polarity intervals longer than about 10,000 years, where the dipole moment has time recover in strength, are equally likely in insulating and CIC models. Since these latter intervals are of more geophysical relevance, we conclude that the influence of the inner-core conductivity on Earth-like reversal sequences is insignificant for the dynamo model employed here.     

Detecting thermal boundary control in surface flows from numerical dynamos

The geomagnetic field and secular variation exhibit asymmetrical spatial features which are possibly originating from an heterogeneous thermal control of the Earth's lower mantle on the core. The identification of this control in magnetic data is subject to several difficulties, some of which can be alleviated by the use of core surface flow models. Using numerical dynamos driven by heterogeneous boundary heat flux, we confirm that within the parameter space accessible to simulations, time average surface flows obey a simple thermal wind equilibrium between the Coriolis and buoyancy forces, the Lorentz, inertial and viscous forces playing only a secondary role, even for Elsasser numbers significantly larger than 1. Furthermore, we average the models over the duration of three vortex turnovers, and correlate them with a longer time average which fully reveals the signature of boundary heterogeneity. This allows us to quantify the possibility of observing mantle control in core surface flows averaged over a short time period. A scaling analysis is performed in order to apply the results to the Earth's core. We find that three vortex turnovers could represent between 100 and 360 years of Earth time, and that the heat flux heterogeneity at the core-mantle boundary could be large enough to yield an observable signature of thermal mantle control in a time average core surface flow within reach of the available geomagnetic data.     

Numerical simulations of MHD dynamos

The generation of magnetic fields in space plasmas and in astrophysics is usually described within the framework of magnetohydrodynamics. Turbulent helical flows produce magnetic fields very efficiently, with correlation length scales larger than those characterizing the flow. Within the context of the solar magnetic cycle, a turbulent dynamo is responsible for the so-called alpha effect, while the Omega effect is associated to the differential rotation of the Sun.We present direct numerical simulations of turbulent magnetohydrodynamic dynamos including two-fluid effects such as the Hall current. More specifically, we study the evolution of an initially weak and small-scale magnetic field in a system maintained in a stationary regime of hydrodynamic turbulence, and explore the conditions for exponential growth of the magnetic energy. In all the cases considered, we find that the dynamo saturates at the equipartition level between kinetic and magnetic energy, and the total energy reaches a Kolmogorov power spectrum.     

A theoretical analysis of the observed variability of the geomagnetic dipole field

We present a detailed analysis of the Sint-800 virtual axial dipole moment (VADM) data in terms of an Ω mean field model of the geodynamo that features a non-steady generation of poloidal from toroidal magnetic field. The result is a variable excitation of the dipole mode and the overtones, and there are occasional dipole reversals. The model permits a theoretical evaluation of the statistical properties of the dipole mode. We show that the model correctly predicts the distribution of the VADM and the autocorrelation function inferred from the Sint-800 data. The autocorrelation technique allows us to determine the turbulent diffusion time τ_d=R²/β of the geodynamo. We find that τ_d is about 10–15 kyr. The model is able to reproduce the observed secular variation of the dipole mode, and the mean time between successive dipole reversals. On the other hand, the duration of a reversal is a factor 2 too long. This could be due to imperfections in the model or to unknown systematics in the Sint-800 data. The use of mean field theory is shown to be selfconsistent.     

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.     

The geodynamo as a bistable oscillator

Abstract

Our intent is to provide a simple and quantitative understanding of the variability of the axial dipole component of the geomagnetic field on both short and long time scales. To this end we study the statistical properties of a prototype nonlinear mean field model. An azimuthal average is employed, so that (1) we address only the axisymmetric component of the field, and (2) the dynamo parameters have a random component that fluctuates on the (fast) eddy turnover time scale. Numerical solutions with a rapidly fluctuating α reproduce several features of the geomagnetic field: (1) a variable, dominantly dipolar field with additional fine structure due to excited overtones, and sudden reversals during which the field becomes almost quadrupolar, (2) aborted reversals and excursions, (3) intervals between reversals having a Poisson distribution. These properties are robust, and appear regardless of the type of nonlinearity and the model parameters. A technique is presented for analysing the statistical properties of dynamo models of this type. The Fokker-Planck equation for the amplitude a of the fundamental dipole mode shows that a behaves as the position of a heavily damped particle in a bistable potential ∝(1 ? a ²)², subject to random forcing. The dipole amplitude oscillates near the bottom of one well and makes occasional jumps to the other. These reversals are induced solely by the overtones. Theoretical expressions are derived for the statistical distribution of the dipole amplitude, the variance of the dipole amplitude between reversals, and the mean reversal rate. The model explains why the reversal rate increases with increasing secular variation, as observed. Moreover, the present reversal rate of the geodynamo, once per (2?3) × 10⁵ year, is shown to imply a secular variation of the axial dipole moment of ～ 15% (about the current value). The theoretical dipole amplitude distribution agrees well with the Sint-800 data.     

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.     

Meteorite paleomagnetism - From magnetic domains to planetary fields and core dynamos

Meteorites represent the earliest records of the evolution of the solar system, providing information on the conditions, processes and chronology for formation of first solids, planetesimals and differentiated bodies. Evidence on the nature of magnetic fields in the early solar system has been derived from chondritic meteorites. Chondrules, which are millimeter sized silicate spherules formed by rapid melting and cooling, have been shown to retain remanent magnetization records dating from the time of chondrule formation and accretion of planetesimals. Studies on different meteorite classes, including ordinary and carbonaceous chondrites, have however provided contrasting results with wide ranges for protoplanetary disk magnetic fields. Developments on instrumentation and techniques for rock magnetic and paleointensity analyses are allowing increased precision. Micromagnetic and an array of geochemical, petrographic and electronic microscopy analyses provide unprecedented resolution, characterizing rock magnetic properties at magnetic domain scales. We review studies on chondrules from the Allende meteorite that reveal relationships among hysteresis parameters and physical properties. Coercivity, remanent and saturation remanence parameters correlate with chondrule size and density; in turn related to internal chondrule structure, mineralogy and morphology. Compound, fragmented and rimmed chondrules show distinct hysteresis properties, related to mineral composition and microstructures. The remanent magnetization record and paleointensity estimates derived from the Allende and other chondrites support remanent acquisition under influence of internal magnetic fields within parent planetesimals. Results support that rapid differentiation following formation of calcium-aluminum inclusions and chondrules gave rise to differentiated planetesimals with iron cores, capable of generating and sustaining dynamo action for million year periods. The Allende chondrite may have derived from a partly differentiated planetesimal which sustained an internal magnetic field.     

A numerical method to compute Euler potentials for non dipolar magnetic fields

The magnetospheric magnetic field may be conveniently described by two scalar functions (α, β), known as the Euler potentials. They are not uniquely defined, and they may be difficult to derive for configuration more complex than a simple dipole. We propose here a simple numerical method to compute one possible pair (α, β). In magnetospheric regions of closed field lines, α can be chosen as a function of the tube volume of unit magnetic flux. The method can be applied to a wide class of magnetic fields which describe the magnetospheric domain of closed field lines and the conjugated ionosphere. Here, it is used with the T87 Tsyganenko model. The results coincide with the dipolar potentials at close distances from the Earth. At larger distances, they display an increasing distortion with the radial distance (or the invariant latitude in the ionosphere) and the magnetic activity. In the magnetosphere, the contours of α and β are stretched towards the nightside. In the ionosphere, they also extend towards the nightside and present major distortions in a narrow ring at the polar cap boundary, which mapsdistant boundary layers in the magnetosphere.     

Control volume method for hydromagnetic dynamos in rotating spherical shells: Testing the code against the numerical dynamo benchmark

Hydromagnetic dynamos in rotating spherical shells are investigated using the control volume method. We present a validation of our code against the numerical dynamo benchmark. It is successfully benchmarked and we are able to conclude that the control volume method is another numerical method available for numerical modelling of self-consistent dynamos. In addition, the efficiency of our numerical code is tested. Computations provide conclusions that dynamo codes based on the spectral methods are much more efficient than our code based on the control volume method at the study of global fields on small and medium size parallel computers. However, our code could be much more efficient than codes based on the spectral methods on very large parallel computers, especially at the study of turbulence.     

A two-step unconditionally stable explicit method with controllable numerical dissipations

A family of unconditionally stable direct integration algorithm with controllable numerical dissipations is proposed. The numerical properties of the new algorithms are controlled by three parameters α, β and γ. By the consistent and stability analysis, the proposed algorithms achieve the second-order accuracy and are unconditionally stable under the condition that α≥-0.5, β≤ 0.5 and γ≥-(1+α)/2. Compared with other unconditionally stable algorithms, such as Chang's algorithms and CR algorithm, the proposed algorithms are found to be superior in terms of the controllable numerical damping ratios. The unconditional stability and numerical damping ratios of the proposed algorithms are examined by three numerical examples. The results demonstrate that the proposed algorithms have a superior performance and can be used expediently in solving linear elastic dynamics problems.     

Block-effective macrodispersion for numerical simulations of sorbing solute transport in heterogeneous porous formations

Heterogeneity is prevalent in aquifers and has an enormous impact on contaminant transport in groundwater. Numerical simulations are an effective way to deal with heterogeneity directly by assigning different hydraulic property values to each numerical grid block. Because hydraulic properties vary on different scales, but they cannot be sampled exhaustively and the number of numerical grid blocks is limited by computational considerations, the dispersive effects of unmodeled heterogeneity need to be accounted for. Dispersion tensors can be used to model the dispersion caused by unmodeled heterogeneity. The concept of block-effective macrodispersion tensors for modeling the effects of small-scale variability on solute transport introduced by Rubin et al. [Rubin Y, Sun A, Maxwell R, Bellin A. The concept of block-effective macrodispersivity and a unified approach for grid-scale- and plume-scale-dependent transport. J Fluid Mech 1999;395:161–80] is extended in this paper for use with reactive solutes. The tensors are derived for reactive solutes with spatially variable retardation factors and for solutes experiencing spatially uniform rate-limited sorption. The longitudinal block-effective macrodispersion coefficient is largest for perfect negative correlation between the log-hydraulic conductivity and the retardation factor. Because dispersion tensors, as they are usually implemented in numerical simulations, produce symmetric spreading, the applicability of the concept depends on the portion of the plume asymmetry caused by small-scale variability. The presented results show that the concept is applicable for rate-limited sorption for block sizes of one and two integral scales.     

Experimental and numerical modelling of sedimentation in a rectangular shallow basin

Numerical simulation of flows in shallow reservoirs has to be checked for its consistency in predicting real flow conditions and sedimentation patterns. Typical flow patterns may exhibit flow separation at the inlet, accompanied by several recirculation and stagnation areas all over the reservoir surface. The aim of the present research project is to study the influence of the geometry of a reservoir on sediment transport and deposition numerically and experimentally, focusing on a prototype reservoir depth between 5 and 15 m as well as suspended sediment transport.
A series of numerical simulations is presented and compared with scaled laboratory experiments, with the objective of testing the sensitivity to different flow and sediment parameters and different turbulence closure schemes. Different scenarios are analyzed and a detailed comparison of preliminary laboratory tests and some selected simulations are presented.
The laboratory experiments show that suspended sediment transport and deposition are determined by the initial flow pattern and by the upstream and downstream boundary conditions. In the experiments, deposition in the rectangular basin systematically developed along the left bank, although inflow and outflow were positioned symmetrically along the centre of the basin. Three major horizontal eddies developed influencing the sediment deposition pattern. Although asymmetric flow patterns are privileged, a symmetric pattern can appear from time to time. This particular behaviour could also be reproduced by a two-dimensional depth-averaged flow and sediment transport model （CCHE2D）. The paper presents numerical simulations using different turbulence closure schemes （k-ε and eddy viscosity models）. In spite of the symmetric setup, these generally produced an asymmetric flow pattern that can easily switch sides depending on the assumptions made for the initial and boundary conditions. When using the laboratory experiment as a reference, the most reliable numerical results have been obtai     

Sensitivity analysis of influencing parameters on slit-type barrier performance against debris flow using 3D-based numerical approach

Extreme rainfall-induced debris flow can be catastrophic to an urban area,and installation of slit-type barriers can prevent such damage while minimizing negative impact on environments.However,the performance of slit-type barriers against debris flows remains poorly identified partly due to the innate complexity in interactions between debris flow and solid structure.This paper investigated the flow behaviors of debris affected by slit-type barriers using the computational fluid dynamics(CFD)method,in which the numerical model based on the volume of fluid method was verified using the physical modeling results.The sensitivity analysis was performed by building metamodels to determine the primary parameters influencing the barrier performance against debris flows among various variables,in which the effect of input properties and design parameters,particularly the soil concentration in fluidized debris,initial velocity and volume of debris,the barrier height,and the opening ratio,was evaluated from the perspectives of the flow energy reduction and debris trapping.The initial velocity and volume of debris were found to play a significant role in determining the debris flow characteristics.A decrease in the opening ratio in the channel primarily facilitated the energy reduction and trapping due to the reduced opening size.However,the barrier height exhibited a limited effect when the height was sufficiently high to block the debris flow volume.In addition,it was observed that the double barrier system effectively increased the energy reduction while keeping the benefit of open-type barrier.The developed simulation method and obtained results provide an effective tool and an insight that can contribute to an optimum design of the debris-flow barrier.     

Two-dimensional numerical modeling of dam-break flows over natural terrain using a central explicit scheme

A two-dimensional (2D) numerical model has been developed to solve shallow water equations for simulation of dam-break flows. The spatial derivatives are discretized using a well-balanced explicit central upwind conservative scheme. The scheme is Riemann solver free and guarantees the positivity of the flow depth over complex topography if the Courant number is kept less than 0.25. The time integration is performed by Euler’s scheme. The model is verified against analytical results for water surface elevation and discharge for three benchmark test cases. A good agreement between analytical solutions and computed results is observed. The property of well-balancing in still water over an uneven bottom is also confirmed. The model is then validated by simulating a laboratory experiment in which a dam break flow propagates over a triangular obstacle. The model performance was found to be satisfactory. A dam break laboratory experimental test case on a frictionless horizontal bottom is also simulated for 2D validation of the model, and good agreement between simulation and the experimental data is observed. The suitability of the proposed model for real life applications is demonstrated by simulating the Malpasset dam-break event, which occurred in 1959 in France. The computed arrival time of the flood wave front and the maximum flow depths at various observation points matched well with the measurements on a 1/400 scale physical model. The overall performance indicates that this model can be applied for simulation of dam-break waves in real life cases.     

Experimental and numerical analysis of MASW tests for detection of buried timber trestles

This paper presents results from multi-channel analysis surface waves (MASW) tests conducted to locate buried timber trestles in two different sections (A and B) of an earth embankment. For section A, the location of the trestles is known as well as the soil properties; thus, it is used for calibration purposes. Conversely, the location of the timber trestles is unknown for section B. A seismic array of 24 geophones with a spacing of 0.5 m is used for surface-wave measurements. Different signal processing techniques are used for data analysis: dispersion curves, power spectral functions, frequency-spectra contour plots, two-dimensional-Fourier transform, and the wavelet transform. A new procedure is proposed for the surface location of buried trestles, which is based on the use of the normalized average power energy plot. Finite-element numerical simulations of MASW tests on layered and homogeneous media with and without buried trestles are performed to support experimental results and to explain the wave–trestle interaction. The numerical simulations show that a buried trestle induces vibration amplifications at the surface in front of its location but vibration attenuation immediately after. These effects, however, are observed only if the embedment depth of timber trestle is smaller than one third of the wavelength. Experimental and numerical results show that the results from MASW tests can be successfully used to define the surface location of decayed buried trestles.     

Development of a family of unconditionally stable explicit direct integration algorithms with controllable numerical energy dissipation

The implicit dissipative generalized‐ α method is analyzed using discrete control theory. Based on this analysis, a one‐parameter family of explicit direct integration algorithms with controllable numerical energy dissipation, referred to as the explicit KR‐α method, is developed for linear and nonlinear structural dynamic numerical analysis applications. Stability, numerical dispersion, and energy dissipation characteristics of the proposed algorithms are studied. It is shown that the algorithms are unconditionally stable for linear elastic and stiffness softening‐type nonlinear systems, where the latter indicates a reduction in post yield stiffness in the force–deformation response. The amount of numerical damping is controlled by a single parameter, which provides a measure of the numerical energy dissipation at higher frequencies. Thus, for a specific value of this parameter, the resulting algorithm is shown to produce no numerical energy dissipation. Furthermore, it is shown that the influence of the numerical damping on the lower mode response is negligible. It is further shown that the numerical dispersion and energy dissipation characteristics of the proposed explicit algorithms are the same as that of the implicit generalized‐ α method. A numerical example is presented to demonstrate the potential of the proposed algorithms in reducing participation of undesired higher modes by using numerical energy dissipation to damp out these modes. Copyright © 2014 John Wiley & Sons, Ltd.     

The effect of intensity of turbulence in umbrella cloud on tephra dispersion during explosive volcanic eruptions: Experimental and numerical approaches

During an explosive volcanic eruption, tephra fall out from the umbrella region of the eruption cloud to the ground surface. We investigated the effect of the intensity of turbulence in the umbrella cloud on dispersion and sedimentation of tephra by performing a series of laboratory experiments and three dimensional (3-D) numerical simulations. In the laboratory experiments, spherical glass-bead particles are mixed in stirred water with various intensities of turbulence, and the spatial distribution and the temporal evolution of the particle concentration are measured. The experimental results show that, when the root-mean-square of velocity fluctuation in the fluid (W_rms) is much greater than the particle terminal velocity (v_t), the particles are homogeneously distributed in the fluid, and settle at their terminal velocities at the base of the fluid where turbulence diminishes. On the other hand, when W_rms is as small as or smaller than v_t, the particle concentration increases toward the base of the fluid during settling, which substantially increases the rate of particle settling. The results of the 3-D simulations of eruption cloud indicate that W_rms is up to 40 m/s in most of the umbrella cloud even during a large scale plinian eruption with a magma discharge rate of 10⁹ kg/s. These results suggest that relatively coarse pyroclasts (more than a few mm in diameter) tend to concentrate around the base of the umbrella cloud, whereas fine pyroclasts (less than 1/8 mm in diameter) may be distributed homogeneously throughout the umbrella cloud during tephra dispersion. The effect of the gradient of particle concentration in the umbrella cloud explains the granulometric data of the Pinatubo 1991 plinian deposits.