Geoid Determination Through Ellipsoidal Stokes Boundary-Value Problem by Splitting Its Solution to the Low-Degree and the High-Degree Parts

The ellipsoidal Stokes boundary-value problem is used to compute the geoidal heights. The low degree part of the geoidal heights can be represented more accurately by Global Geopotential Models (GGM). So the disturbing potential is splitted into a low-degree reference potential and a higher-degree potential. To compute the low-degree part, the global geopotential model is used, and for the high-degree part, the solution of the ellipsoidal Stokes boundary-value problem in the form of the surface integral is used. We present an effective method to remove the singularity of the high-degree of the spherical and ellipsoidal Stokes functions around the computational point. Finally, the numerical results of solving the ellipsoidal Stokes boundary-value problem and the difference between the high-degree part of the solution of the ellipsoidal Stokes boundary-value problem and that of the spherical Stokes boundary-value problem is presented.     

Ellipsoidal Stokes Boundary-Value Problem with Ellipsoidal Corrections in the Boundary Condition

We would like to solve the Stokes boundary-value problem taking into consideration the ellipsoidal corrections in the boundary condition in ellipsoidal coordinates The original problem, i.e., the ellipsoidal Stokes boundary-value problem has been solved by Martinec and Grafarend (1997) We use the same philosophy expressed by Martinec (1998) to solve the spherical Stokes boundary-value problem with ellipsoidal corrections in the boundary condition We wish to show the magnitude of the integration kernel describing the effect of the ellipsoidal corrections in the boundary condition in a cap around the computational point.     

Solution to the Stokes Boundary-Value Problem on an Ellipsoid of Revolution

We have constructed Green's function to Stokes's boundary-value problem with the gravity data distributed over an ellipsoid of revolution. We show that the problem has a unique solution provided that the first eccentricity e₀ of the ellipsoid of revolution is less than 0·65041. The ellipsoidal Stokes function describing the effect of ellipticity of the boundary is expressed in the E-approximation as a finite sum of elementary functions which describe analytically the behaviour of the ellipsoidal Stokes function at the singular point = 0. We prove that the degree of singularity of the ellipsoidal Stokes function in the vicinity of its singular point is the same as that of the spherical Stokes function.     

Far-Zone Contribution in Ellipsoidal Stokes Boundary-Value Problem

In the first attempt to solve the Stokes boundary-value problem in ellipsoidal coordinates numerically (Ardestani and Martinec, 2000), we focused on the near-zone contribution since the effect of the ellipsoidal Stokes function in the far-zone contribution is not considered. We present a method for solving the ellipsoidal Stokes integral in far-zone contribution. The numerical results of computing the magnitude of this term for an area in north of Canada are presented.     

Maximum Entropy Density Determination for the Earth and Terrestrial Planets

The maximum entropy principle of the information theory gives rise to a general regularization strategy for ill-posed inverse problems. The methods based on this principle have become standard in various branches of engineering sciences. Of course, ill-posed problems frequently appear in Earth sciences, too. Nonetheless, the concept of maximum entropy is not very popular here. Therefore, we review the basic approaches employing the principle of maximum entropy in one way or another. We can distinguish at least three different approaches, partly yielding coincident results. One possible area of application is the determination of Earth and planetary models, although the paper cannot treat this in its practical complexity. Most of the discussion is restricted to the determination of the Earth's mass density function from various sources of data. Three sample problems are solved using the principle of maximum entropy: a spherical and an ellipsoidal problem related to the Earth and an ellipsoidal problem related to Mars. This illustrates the numerical procedure, which is non-trivial in many cases. It also shows some results, partly compared to standard solutions. The pros and cons of the approaches are discussed.     

Reproducing kernel and Galerkin��s matrix for the exterior of an ellipsoid: Application in gravity field studies

In the introductory part of the paper the importance of the topic for gravity field studies is outlined. Some concepts and tools often used for the representation of the solution of the respective boundary-value problems are mentioned. Subsequently a weak formulation of Neumann??s problem is considered with emphasis on a particular choice of function basis generated by the reproducing kernel of the respective Hilbert space of functions. The paper then focuses on the construction of the reproducing kernel for the solution domain given by the exterior of an oblate ellipsoid of revolution. First its exact structure is derived by means of the apparatus of ellipsoidal harmonics. In this case the structure of the kernel, similarly as of the entries of Galerkin??s matrix, becomes rather complex. Therefore, an approximation of ellipsoidal harmonics (limit layer approach), based on an approximation version of Legendre??s ordinary differential equation, resulting from the method of separation of variables in solving Laplace??s equation, is used. The kernel thus obtained shows some similar features, which the reproducing kernel has in the spherical case, i.e. for the solution domain represented by the exterior of a sphere. A numerical implementation of the exact structure of the reproducing kernel is mentioned as a driving impulse of running investigations.     

Atmospheric effects in the derivation of geoid-generated gravity anomalies

Parameters of the gravity field harmonics outside the geoid are sought in solving the Stokes boundary-value problem while harmonics outside the Earth in solving the Molodensky boundary-value problem. The gravitational field generated by the atmosphere is subtracted from the Earth’s gravity field in solving either the Stokes or Molodensky problem. The computation of the atmospheric effect on the ground gravity anomaly is of a particular interest in this study. In this paper in particular the effect of atmospheric masses is discussed for the Stokes problem. In this case the effect comprises two components, specifically the direct and secondary indirect atmospheric effects. The numerical investigation is conducted at the territory of Canada. Numerical results reveal that the complete effect of atmosphere on the ground gravity anomaly varies between 1.75 and 1.81 mGal. The error propagation indicates that precise determination of the atmospheric effect on the gravity anomaly depends mainly on the accuracy of the atmospheric mass density distribution model used for the computation.     

Magnetic instabilities in rapidly rotating spherical geometries I. from cylinders to spheres

Abstract

The solution of the full nonlinear hydromagnetic dynamo problem is a major numerical undertaking. While efforts continue, supplementary studies into various aspects of the dynamo process can greatly improve our understanding of the mechanisms involved. In the present study, the linear stability of an electrically conducting fluid in a rigid, electrically insulating spherical container in the presence of a toroidal magnetic field B_o(r,θ)lø and toroidal velocity field U_o(r,θ)lø, [where (r,θ,ø) are polar coordinates] is investigated. The system, a model for the Earth's fluid core, is rapidly rotating, the magnetostrophic approximation is used and thermal effects are excluded. Earlier studies have adopted a cylindrical geometry in order to simplify the numerical analysis. Although the cylindrical geometry retains the fundamental physics, a spherical geometry is a more appropriate model for the Earth. Here, we use the results which have been found for cylindrical systems as guidelines for the more realistic spherical case. This is achieved by restricting attention to basic states depending only on the distance from the rotation axis and by concentrating on the field gradient instability. We then find that our calculations for the sphere are in very good qualitative agreement both with a local analysis and with the predictions from the results of the cylindrical geometry. We have thus established the existence of field gradient modes in a realistic (spherical) model and found a sound basis for the study of various other, more complicated, classes of magnetically driven instabilities which will be comprehensively investigated in future work.     

A refined method of recovering potential coefficients from surface gravity data

Summary A new method for computing the potential coefficients of the Earth's external gravity field is presented. The gravimetric boundary-value problem with a free boundary is reduced to the problem with a fixed known telluroid. The main idea of the derivation consists in a continuation of the quantities from the physical surface to the telluroid by means of Taylor's series expansion in such a way that the terms whose magnitudes are comparable with the accuracy of today's gravity measurements are retained. Thus not only linear, but also non-linear terms are taken into account. Explicitly, the terms up to the order of the third power of the Earth's flattening are retained. The non-linear boundary-value problem on the telluroid is solved by an iteration procedure with successive approximations. In each iteration step the solution of the non-linear problem is estimated by the solutions of two linear problems utilizing the fact that the non-linear boundary condition may be split into two parts; the linear spherical approximation of the gravity anomaly whose magnitude is significantly greater than the others and the non-linear ellipsoidal corrections. Finally, in order to solve the problem in terms of spherical harmonics, the transform method composed of the fast Fourier transform and Gauss Legendre quadrature is theoretically outlined. Immediate data processing of gravity data measured on the physical Earth's surface without any continuation of gravity measurements to a reference level surface belongs to the main advantage of the presented method. This implies that no preliminary data handling is needed and that the error data propagation is, consequently, maximally suppressed.     

Use of Green’s function for determining the disturbing potential of an ellipsoidal Earth

A spherical approximation makes the basis for a majority of formulas in physical geodesy. However, the present-day accuracy in determining the disturbing potential requires an ellipsoidal approximation. The paper deals with constructing Green’s function for an ellipsoidal Earth by an ellipsoidal harmonic expansion and using it for determining the disturbing potential. From the result obtained the part that corresponds to the spherical approximation has been extracted. Green’s function is known to depend just on the geometry of the surface where boundary values are given. Thus, it can be calculated irrespective of the gravity data completeness. No changes of gravity data have an effect on Green’s function and they can be easily taken into account if the function has already been constructed. Such a method, therefore, can be useful in determining the disturbing potential of an ellipsoidal Earth.     

A note on three-dimensional scattering and diffraction by a hemispherical canyon–I: Vertically incident plane P-wave

The three-dimensional scattering by a hemi-spherical canyon in an elastic half-space subjected to seismic plane and spherical waves has long been a challenging boundary-value problem. It has been studied by earthquake engineers and strong-motion seismologists to understand the amplification effects caused by surface topography. The scattered and diffracted waves will, in all cases, consist of both longitudinal (P-) and shear (S-) shear waves. Together, at the half-space surface, these waves are not orthogonal over the infinite plane boundary of the half-space. Thus, to simultaneously satisfy both zero normal and shear stresses on the plane boundary numerical approximation of the geometry and/or wave functions were required, or in some cases, relaxed (disregarded). This paper re-examines this boundary-value problem from the applied mathematics point of view, and aims to redefine the proper form of the orthogonal spherical-wave functions for both the longitudinal and shear waves, so that they can together simultaneously satisfy the zero-stress boundary conditions at the half-space surface. With the zero-stress boundary conditions satisfied at the half-space surface, the most difficult part of the problem will be solved, and the remaining boundary conditions at the finite canyon surface will be easy to satisfy.     

Electron transport and energy degradation in the ionosphere: evaluation of the numerical solution,comparison with laboratory experiments and auroral observations

Auroral electron transport calculations are a critical part of auroral models. We evaluate a numerical solution to the transport and energy degradation problem. The numerical solution is verified by reproducing simplified problems to which analytic solutions exist, internal self-consistency tests, comparison with laboratory experiments of electron beams penetrating a collision chamber, and by comparison with auroral observations, particularly the emission ratio of the N₂ second positive to N⁺ ₂ first negative emissions. Our numerical solutions agree with range measurements in collision chambers. The calculated N₂2P to N⁺ ₂1N emission ratio is independent of the spectral characteristics of the incident electrons, and agrees with the value observed in aurora. Using different sets of energy loss cross sections and different functions to describe the energy distribution of secondary electrons that emerge from ionization collisions, we discuss the uncertainties of the solutions to the electron transport equation resulting from the uncertainties of these input parameters.     

Evolution of non-linear α 2-dynamos and taylor's constraint

A key non-linear mechanism in a strong-field geodynamo is that a finite amplitude magnetic field drives a flow through the Lorentz force in the momentum equation and this flow feeds back on the field-generation process in the magnetic induction equation, equilibrating the field. We make use of a simpler non-linear?α?²-dynamo to investigate this mechanism in a rapidly rotating fluid spherical shell. Neglecting inertia, we use a pseudo-spectral time-stepping procedure to solve the induction equation and the momentum equation with no-slip velocity boundary conditions for a finitely conducting inner core and an insulating mantle. We present calculations for Ekman numbers (E) in the range 2.5× 10^?3 to 5.0× 10^?5, for?α?=α ₀cos?θ?sin?π?(r?r_i ) (which vanishes on both inner and outer boundaries). Solutions are steady except at lower E and higher values of?α?₀. Then they are periodic with a reversing field and a characteristic rapid increase then equally rapid decrease in magnetic energy. We have investigated the mechanism for this and shown the influence of Taylor's constraint. We comment on the application of our findings to numerical hydrodynamic dynamos.     

Modeling of Physical Processes by Analysis of Hard X-Ray and Microwave Radiations in the Solar Flare of November 10, 2002

For electron acceleration during solar flares, it is very important to determine the pitch-angle and energy dependences of the electron distribution function. At present, this cannot be done directly from observations. Therefore, it is necessary to perform a numerical simulation of the propagation of accelerated electrons in the magnetic field of the flare loop (loops) and calculate the X-ray and radio emissions. For the solar flare of November 10, 2002, we have obtained qualitative and quantitative agreements of modeled X-ray and radio maps with the RHESSI satellite and Nobeyama Radioheliograph data. We have determined the flare model parameters that agree with observations. The pitch-angle anisotropy of electrons determined by highly directional functions of the S(α) = cos⁸(α) type, the energy spectrum consist of two electron populations, the low-energy part of the spectrum up to an energy of break of 350 keV is characterized by a power law with the exponent δ₁ = 2.7–2.9, and the energy spectrum is more rigid above 420 keV (δ₂ = 2–2.3).     

The boundary value problem of Poisson’s equation with Stokes’ boundary condition

The following Poisson’s equation with the Stokes’ boundary condition is dealt with $$\left\{ \begin{gathered} \nabla ^2 T = - 4\pi Gp outside S, \hfill \\ \left. {\frac{{\partial T}}{{\partial h}} = \frac{1}{\gamma }\frac{{\partial y}}{{\partial h}}T} \right|_s = - \Delta g, \hfill \\ T = O\left( {r^{ - 3} } \right) at infinity, \hfill \\ \end{gathered} \right.$$ whereS is reference ellipsord. Under spherical approximation transformation, the ellipsoidal correction terms about the boundary condition, the equation and the density in the above BVP are respectively given. Therefore, the disturbing potentialT can he obtained if the magnitudes aboveO(ε⁴) are neglected.     

The characteristics of magma reservoir failure beneath a volcanic edifice

Eruptions fed from subsurface reservoirs commonly construct volcanic edifices at the surface, and the growth of an edifice will in turn modify the subsurface stress state that dictates the conditions under which subsequent rupture of the inflating reservoir can occur. We re-examine this problem using axisymmetric finite element models of ellipsoidal reservoirs beneath conical edifices, explicitly incorporating factors (e.g., full gravitational loading conditions, an elastic edifice instead of a surface load, reservoir pressures sufficient to induce tensile rupture) that compromise previous solutions to illustrate why variations in rupture behavior can occur. Relative to half-space model results, the presence of an edifice generally rotates rupture toward the crest of a spherical reservoir, with increasing flank slope (for an edifice of constant volume) and larger edifices (or greater reservoir scaled depths) normally serving to enhance this trend. When non-spherical reservoirs are considered, the presence of an edifice amplifies previously identified half-space failure characteristics, shifting rupture to the crest more rapidly for prolate reservoirs while forcing rupture closer to the midpoint of oblate reservoirs. Rupture is always observed to occur in the σ_t orientation, and depending on where initial failure occurs rupture favors the initial emplacement of either lateral sills, circumferential intrusions or vertically ascending dikes. Ultimately, integration of our numerical model results with other information, for instance the sequence of intrusion/eruption events observed at a given volcano, can provide useful new insight into how a volcano's subsurface magma plumbing system evolved. We demonstrate this process through application of our model to Summer Coon, a well-studied stratocone on Earth, and Ilithyia Mons, a large conical shield volcano on Venus.     

Layer-Based Modelling of the Earth’s Gravitational Potential up to 10-km Scale in Spherical Harmonics in Spherical and Ellipsoidal Approximation

Global forward modelling of the Earth’s gravitational potential, a classical problem in geophysics and geodesy, is relevant for a range of applications such as gravity interpretation, isostatic hypothesis testing or combined gravity field modelling with high and ultra-high resolution. This study presents spectral forward modelling with volumetric mass layers to degree 2190 for the first time based on two different levels of approximation. In spherical approximation, the mass layers are referred to a sphere, yielding the spherical topographic potential. In ellipsoidal approximation where an ellipsoid of revolution provides the reference, the ellipsoidal topographic potential (ETP) is obtained. For both types of approximation, we derive a mass layer concept and study it with layered data from the Earth2014 topography model at 5-arc-min resolution. We show that the layer concept can be applied with either actual layer density or density contrasts w.r.t. a reference density, without discernible differences in the computed gravity functionals. To avoid aliasing and truncation errors, we carefully account for increased sampling requirements due to the exponentiation of the boundary functions and consider all numerically relevant terms of the involved binominal series expansions. The main outcome of our work is a set of new spectral models of the Earth’s topographic potential relying on mass layer modelling in spherical and in ellipsoidal approximation. We compare both levels of approximations geometrically, spectrally and numerically and quantify the benefits over the frequently used rock-equivalent topography (RET) method. We show that by using the ETP it is possible to avoid any displacement of masses and quantify also the benefit of mapping-free modelling. The layer-based forward modelling is corroborated by GOCE satellite gradiometry, by in-situ gravity observations from recently released Antarctic gravity anomaly grids and degree correlations with spectral models of the Earth’s observed geopotential. As the main conclusion of this work, the mass layer approach allows more accurate modelling of the topographic potential because it avoids 10–20-mGal approximation errors associated with RET techniques. The spherical approximation is suited for a range of geophysical applications, while the ellipsoidal approximation is preferable for applications requiring high accuracy or high resolution.     

Numerically quantifying energy loss caused by squirt flow

In interconnected microcracks, or in microcracks connected to spherical pores, the deformation associated with the passage of mechanical waves can induce fluid flow parallel to the crack walls, which is known as squirt flow. This phenomenon can also occur at larger scales in hydraulically interconnected mesoscopic cracks or fractures. The associated viscous friction causes the waves to experience attenuation and velocity dispersion. We present a simple hydromechanical numerical scheme, based on the interface-coupled Lamé–Navier and Navier–Stokes equations, to simulate squirt flow in the frequency domain. The linearized, quasi-static Navier–Stokes equations describe the laminar flow of a compressible viscous fluid in conduits embedded in a linear elastic solid background described by the quasi-static Lamé–Navier equations. Assuming that the heterogeneous model behaves effectively like a homogeneous viscoelastic medium at a larger spatial scale, the resulting attenuation and stiffness modulus dispersion are computed from spatial averages of the complex-valued, frequency-dependent stress and strain fields. An energy-based approach is implemented to calculate the local contributions to attenuation that, when integrated over the entire model, yield results that are identical to those based on the viscoelastic assumption. In addition to thus validating this assumption, the energy-based approach allows for analyses of the spatial dissipation patterns in squirt flow models. We perform simulations for a series of numerical models to illustrate the viability and versatility of the proposed method. For a 3D model consisting of a spherical crack embedded in a solid background, the characteristic frequency of the resulting P-wave attenuation agrees with that of a corresponding analytical solution, indicating that the dissipative viscous flow problem is appropriately handled in our numerical solution of the linearized, quasi-static Navier–Stokes equations. For 2D models containing either interconnected cracks or cracks connected to a circular pore, the results are compared with those based on Biot's poroelastic equations of consolidation, which are solved through an equivalent approach. Overall, our numerical simulations and the associated analyses demonstrate the suitability of the coupled Lamé–Navier and Navier–Stokes equations and of Biot's equations for quantifying attenuation and dispersion for a range of squirt flow scenarios. These analyses also allow for delineating numerical and physical limitations associated with each set of equations.     

A Brovar-type solution of the fixed geodetic boundary-value problem

For more than 150 years gravity anomalies have been used for the determination of geoidal heights, height anomalies and the external gravity field. Due to the fact that precise ellipsoidal heights could not be observed directly, traditionally a free geodetic boundary-value problem (GBVP) had to be formulated which after linearisation is related to gravity anomalies. Since nowadays the three-dimensional positions of gravity points can be determined by global navigation satellite systems very precisely, the modern formulation of the GBVP can be based on gravity disturbances which are related to a fixed GBVP using the known topographical surface of the Earth as boundary surface. The paper discusses various approaches into the solution of the fixed GBVP which after linearization corresponds to an oblique-derivative boundary-value problem for the Laplace equation. Among the analytical solution approaches a Brovar-type solution is worked out in detail, showing many similarities with respect to the classical solution of the scalar free GBVP.     

Estimation of the gravitational potential energy of the earth based on different density models

The estimation of the Earth’s gravitational potential energy E was obtained for different density distributions and rests on the expression E = − (W_min + ΔW) derived from the conventional relationship for E. The first component W_min expresses minimum amount of the work W and the second component ΔW represents a deviation from W_min interpreted in terms of Dirichlet’s integral applied on the internal potential. Relationships between the internal potential and E were developed for continuous and piecewise continuous density distributions. The global 3D density model inside an ellipsoid of revolution was chosen as a combined solution of the 3D continuous distribution and the reference PREM radial piecewise continuous profile. All the estimates of E were obtained for the spherical Earth since the estimated (from error propagation rule) accuracy σ_E of the energy E is at least two orders greater than the ellipsoidal reduction and the contribution of lateral density inhomogeneities of the 3D global density model. The energy E contained in the 2nd degree Stokes coefficients was determined. A good agreement between E = E_Gauss derived from Gaussian distribution and other E, in particular for E = E_PREM based on the PREM piecewise continuous density model and E-estimates derived from simplest Legendre-Laplace, Roche, Bullard and Gauss models separated into core and mantle only, suggests the Gaussian distribution as a basic radial model when information about density jumps is absent or incomplete.