On inversion for mass distribution from global (time-variable) gravity field
|
|
- Emmeline Wilkinson
- 5 years ago
- Views:
Transcription
1 Journal of Geodynamics 39 (2005) On inversion for mass distribution from global (time-variable) gravity field Benjamin F. Chao Space Geodesy Branch, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA Received 20 July 2004; received in revised form 9 November 2004; accepted 10 November 2004 Abstract The well-known non-uniqueness of the gravitational inverse problem states that the external gravity field, even if completely and exactly known, cannot uniquely determine the density distribution of the body that produces the gravity field. In this paper, we provide conceptual insight by examining the problem in terms of spherical harmonic expansion of the global gravity field. By comparing the multipoles and the moments of the density function, we show that in 3-D the degree of knowledge deficiency in trying to inversely recover the density distribution from an external gravity field solution is (n + 1)(n + 2)/2 (2n +1)=n(n 1)/2 for each harmonic degree n. On the other hand, on a 2-D spherical shell we show via a simple relationship that the inverse solution of the surface density distribution is unique. The latter applies quite readily in the inversion of time-variable gravity signals (such as those observed by the GRACE space mission) where the sources largely come from the Earth s surface over a wide range of timescales Elsevier Ltd. All rights reserved. Keywords: Inverse problem; Solution uniqueness; Gravity; Time-variable gravity; Spherical harmonics 1. Introduction The well-known non-uniqueness of the gravitational inverse problem states the following: The external gravity field, even if completely and exactly known, cannot uniquely determine the density distribution of the body that produces the gravity field. This is an intrinsic property of a field that obeys the Laplace equation, as already treated in mathematical as well as geophysical literature (e.g., Menke, 1989; Zidarov, 1990; Parker, 1994; Backus et al., 1996). In a simple example, imagine two uniform concentric spherical address: benjamin.f.chao@nasa.gov /$ see front matter 2005 Elsevier Ltd. All rights reserved. doi: /j.jog
2 224 B.F. Chao / Journal of Geodynamics 39 (2005) shells, one with mass m and the other m, superimposing onto an arbitrary body. It is obvious that the two shells collectively have zero effect on the external gravity field, and hence their existence cannot be inferred from outside gravitationally. The same assertion applies to the time-variable gravity (TVG) produced by a mass distribution varying with time. Global TVG signals detected from space contain information about mass redistributions that occur in the Earth system. Hence, it is becoming a new data type for monitoring climatic and geophysical changes. For example, the GRACE mission (Tapley et al., 2004) is to obtain a high-precision, global gravity map every month or so, at a spatial resolution of several hundred km. The month-to-month difference then signifies mass redistributions in or on the Earth during that given month. What exactly would such TVG data tell us about mass redistribution? For one thing, the TVG is the sum of the gravitational signals originating from all geophysical sources at work at any given time. Sorting out different geophysical signals in the data is a challenge (as well as an art ), but in principle can be facilitated by recognizing the different temporal and spatial characteristics of different geophysical phenomena (e.g., Chao, 1994). This is a problem depending on the real-world geophysics it is not the subject of the present paper. What we want to examine, and we do it analytically, is something more basic conceptually the uniqueness of the general gravitational inversion in terms of spherical harmonic functions for the global field. We shall assume that we have the complete and exact knowledge of the external gravity field (which in fact only constitutes such knowledge on a closed boundary thanks to the Laplace Equation), and we will not deal with practical issues such as solution stability, error propagation, or optimization schemes. The mathematical essence of various aspects of the relevant theory has appeared in variant forms in the literature; in that sense we here cast the theory in geophysical perspective, aiming towards a comprehensive appreciation and understanding of the problem. We shall also treat some situations from a practical viewpoint of TVG, particularly if the latter originates from the surface of the Earth. 2. Gravitational inverse problem in 3-D According to Newton s gravitational law, the external gravity potential field U produced by a (finite) body of internal density distribution ρ over volume is given by the volume integral: ρ(r 0 ) U(r) = G r r 0 d (1) where G is the universal gravitational constant. In the spherical coordinates, r = (radius r, co-latitude θ, longitude λ)=(r, solid angle Ω), if we expand into spherical harmonic functions by the addition theorem the function 1/ r r 0 = n n=0 m= n (2n + 1) 1 (r0 n/rn+1 )Ynm (Ω)Y nm(ω 0 ), we get the so-called multipole expansion of U (e.g., Jackson, 1975): U(r) = G n n=0 m= n 1 (2n + 1)r n+1 [ ρ(r 0 )r n 0 Y nm(ω 0 )d ] Ynm (Ω) (2) for a singly connected, where Y nm is the 4 -normalized surface spherical harmonic function of degree n and order m, and the volume integral in the bracket is the corresponding (complex-valued) multipole of the density distribution ρ(r) where dv =r 2 dω dr =r 2 sinθ dθ dλ dr.
3 B.F. Chao / Journal of Geodynamics 39 (2005) At the same time, because U satisfies the Laplace equation exterior of, it has a closed-form solution customarily expressed as (e.g., Kaula, 1966): U(r, θ, λ) = GM n ( a ) n+1pnm (cos θ)(c nm cos mλ + S nm sin mλ) (3) a r n=0 m=0 where M is the Earth s mass, P nm is the 4 -normalized Legendre function: Y nm (Ω)=P nm (cosθ) exp(imλ). The dimensionless coefficients C nm and S nm are known as the (normalized) Stokes coefficients of degree n (= 0,1,2,..., ) and order m (= 0,1,2,...,n); they constitute the values given to specify a gravity field model. The value of the Stokes coefficients actually refers to the value of some convenient length parameter a; in the case of the Earth a is often chosen to be the mean radius of the Earth along the Equator. It really does not matter what a is as long as r = a is conveniently exterior of the body, because once the Stokes coefficients are known with respect to the chosen a, the entire external gravity field is completely determined all one needs to do is to upward or downward continue, as the case may be, via the radial transfer function (a/r) n +1 for each degree n according to Eq. (3). Comparing Eqs. (2) and (3), one gets 1 C nm + is nm = ρ(r)r n Y (2n + 1)Ma n nm (Ω)dV (4) (where for brevity the subscript 0 in the integrand is omitted). That is, the Stokes coefficients are simply normalized multipoles of the density distribution of the body (e.g., Chao and Gross, 1987; Chao, 1994). The gravitational inverse problem in the present context, Eq. (4), is a canonical linear inversion: given a complete and exact knowledge of the gravity field in the form of the (infinite) set of Stokes coefficients, how much can we learn about ρ(r) in by way of Eq. (4)? All we know are the multipoles, i.e., the projections of ρ(r) onto the solid spherical harmonics functions r n Y nm (Ω) in a Hilbert space H that contains all possible ρ(r) functions. However, the (infinite) set of solid spherical harmonics does not form a complete set of basis function for H it only spans the harmonic sub-space of H, leaving an orthogonal or null sub-space wherefrom any element when added to any solution forms a new solution but having the same Stokes coefficients. To proceed we note that in the Cartesian coordinates one can invoke another expansion for 1/ r r 0, namely the Taylor expansion: 1 r r 0 = ( 1) n [ x 0 n! x + y 0 y + z 0 z n=0 ] n ( ) 1 r Now invoking the multinomial theorem on Eq. (5) and substituting it into Eq. (1), we get another expression for the solution of U (Morse and Feshbach, 1953, pp ): U(r) = G n=0 α+β+γ=n α,β,γ 0 ( 1) n [ α!β!γ! x α 0 yβ 0 zγ 0 ρ(r 0)d ] n ( ) 1 x α y β z γ r We shall call Eq. (6) the moment expansion for U, because the quantity in the bracket is recognized as the (3-D) moment of order α + β + γ = n for the density distribution ρ(r) in. The multipole expansion (2) and the moment expansion (6) are, of course, equivalent. Equating the degree of multipole with the order of moment (we have actually used the same symbol (5) (6)
4 226 B.F. Chao / Journal of Geodynamics 39 (2005) n), we note the following: For each n, there are 2n + 1 independent multipoles while there are (n + 1)(n + 2)/2 independent moments. [The latter is because the form x α y β z γ where α + β + γ = n has (n +1)=(n + 1)(n + 2)/2 terms.] This seeming paradox can be resolved by realizing that certain (in fact the majority of combination of) terms in (6) actually have zero contribution to U (see Morse and Feshbach, 1953, pp ); that is, they belong to the null sub-space in H with respect to U. On the other hand, it is of central importance in the present discussion that, as far as ρ(r) is concerned, the moments uniquely determine ρ(r). In other words, knowing the infinite set of moments of ρ(r) is equivalent to knowing ρ(r) itself. This is known as the Moment Theorem in statistics literature (e.g., Papoulis, 1965, p. 214), and is closely related to the Weierstrass Approximation Theorem (e.g., Backus et al., 1996). In the present context, the moments constitute of projections onto a complete basis function set that spans H. Thus, in trying to determine the moments but only knowing the multipoles, our degree of deficiency in knowledge is (n + 1)(n + 2)/2 (2n +1)=n(n 1)/2 for each n, where n = 0,1,2,...,. Conceptually, one can decompose the Hilbert space H into disjoint sub-spaces, H = H n, each one having a null sub-space of dimension n(n 1)/2. For n = 0, there is only one term in both multipole and moment expansions, i.e., the monopole = zeroth moment, corresponding to the total mass of the Earth. A complete determination is also assured for n = 1, where there are three terms for both expansions: the three degree-1, or dipole, Stokes coefficients correspond to the three first moments of ρ(r) which give the three components of the position vector of the body s center of mass, or the so-called geocenter in the case of the Earth. For n = 2, we have five degree-2, or quadrupole, Stokes coefficients but six second moments, i.e. the elements of the inertia tensor. The multipoles are no longer sufficient to determine the corresponding moments (e.g., Chao and Gross, 1987). The degree of deficiency grows rapidly as n increases to infinity; its numerability provides a quantitative assessment of how grossly non-unique the 3-D gravitational inversion is. The situation is summarized in Table 1. This degree of deficiency is associated with the freedom in determining the radial dependence of ρ(r), as will be shown in the next section. We should mention that, the multipoles (Eq. (2)) and the moments (Eq. (6)) being linear in ρ(r), conceptually we can replace ρ(r) by just the lateral density anomaly, i.e., the laterally heterogeneous part of ρ(r) after removing the dominant, spherically symmetric part (which only translates into the n =0 term). This scheme is useful in practical applications, because then under certain extra constraints the inversion of Eq. (4) can become unique so as to yield unique solutions for ρ(r). For example, applying Table 1 The degree of deficiency as a function of spherical harmonic degree n in the 3-D gravitational inversion Degree n No. of multipoles (2n + 1) No. of moments (n + 1)(n + 2)/2 Degree of deficiency n(n 1)/2 0 1 (monopole) 1 (total mass) (dipole) 3 (center of mass) (quadrupole) 6 (inertia tensor) (octupole) 10 (3rd moment)
5 B.F. Chao / Journal of Geodynamics 39 (2005) the physical condition of minimum shear energy in the mantle, Kaula (1963) determined a unique lateral density distribution for the Earth. Rubincam (1982) achieved the unique solution by employing the mathematical constraint of maximum entropy. Likewise, it can be shown that minimizing the norm-2 variance for the lateral distribution also leads to unique solutions. However, the variability of the lateral density obtained in the latter two cases is highly concentrated towards the surface as a consequence of the imposed mathematical conditions that do not necessarily reflect the reality. 3. Gravitational inverse problem on a 2-D spherical shell Anticipating the TVG case below, let us now consider a special case a gravitating body that is a spherical shell with infinitesimal thickness, called S 0. Let the (2-D) surface mass density distribution over S 0 be σ(ω). How much of the σ(ω) function can we recover given a complete knowledge of external Stokes coefficients as above? Intuitively, we no longer need to recover quite as much as in the 3-D case, because no radial dependence exists. We now show that such σ(ω) can be uniquely determined, and that the above non-uniqueness in 3-D is thus associated with the radial dependence. The proof of this uniqueness for the gravitational inversion over S 0 is rather simple, as follows. Eq. (4) reduces to C nm + is nm = a 2 (2n + 1)M S 0 σ(ω)y nm (Ω)dΩ (7) on S 0, where we have let a simply be the radius of the shell. Now note that the integral on the right side of Eq. (7), thanks to the orthogonality of spherical harmonics, is simply the coefficient of the (n,m)th harmonic component of σ(ω), that is, the σ nm (apart from a multiplicative constant) in the spherical harmonic expansion: σ(ω) = Σ n,m σ nm Ynm (Ω). Eq. (7) can then be written as σ nm = (2n + 1)M 4πa 2 (C nm + is nm ) (8) Eq. (8) states that each of the harmonic components of σ(ω) is determined by the corresponding harmonic component of the external gravity field (i.e. Stokes coefficient), by way of the transfer function 2n + 1 (and a multiplicative constant), which increases with increasing n. This seemingly trivial fact that Eq. (7) is invertible into (8) is actually quite profound: knowing the set of the gravitational Stokes coefficients is equivalent to knowing the spherical harmonic expansion of σ(ω), which in turn is equivalent to knowing σ(ω) itself (because the spherical harmonic functions form a complete set for S 0 ). An equivalent way of stating the above is that the surface spherical harmonics span the Hilbert space containing all 2-D σ(ω) functions defined on S 0. Alternatively, one can state that on a 2-D spherical surface the moment expansion has 2n + 1 terms, the same as the number of independent multipoles for any given n. Therefore, the 2-D spherical σ(ω) can be uniquely determined gravitationally. A corollary is that one can always construct a σ(ω)ons 0 to mimic the external field created by some arbitrary mass distribution in (e.g., Menke, 1989, p. 270). We shall elaborate on this point when discussing the TVG below.
6 228 B.F. Chao / Journal of Geodynamics 39 (2005) Time-variable gravity problem: 3-D body versus 2-D spherical shell When the mass distribution in the body changes with time, the external gravity field would change accordingly. This is the forward problem. The inverse problem is to recover to the extent possible the mass redistribution in from the observed TVG field based on Eqs. (4) and (7). There are two ways to describe a mass transport in a continuum: the Eulerian approach and the Lagrangian approach. They are conceptually equivalent to each other. Here for the sake of argument, we shall adopt the Eulerian description because of its simpler mathematics. In the Eulerian description for a 3-D body, the TVG is expressed in terms of time-variable Stokes coefficients (e.g., Chao, 1994): 1 C nm (t) + i S nm (t) = ρ(r; t)r n Y (2n + 1)Ma n nm (Ω)dV (9) The corresponding 2-D spherical shell case reduces to a 2 C nm (t) + i S nm (t) = σ(ω; t)y nm (Ω)dΩ (10) (2n + 1)M S 0 or equivalently for σ(ω; t) = Σ n,m σ nm (t)ynm (Ω), as Eq. (8), (2n + 1)M σ nm (t) = [ C 4πa 2 nm (t) + i S nm (t)] (11a) As before, Eq. (9) is not invertible and cannot uniquely determine the bodily density change ρ(r;t)in. In contrast, without the radial dependence Eq. (10) is completely invertible to yield a unique surface density change σ(ω;t)ons 0 as specified in Eqs. (11a) or (11b) below. A caveat should be injected at this point. In applying Eq. (11a), one should ascertain that σ truly represents the surface value. If, for example, σ is evaluated from an apparent surface water load in terms of water depth, then σ will be accompanied by an additional (and in general non-superficial) deformation in the solid Earth due to sold Earth s elastic yielding, which is proportional to the load itself. The factor of the proportionality, different for each n, is a negative fractional number known as the load Love number k n. The net effect is that σ is somewhat reduced by a factor of (1 + k n ), about 0.70 for n = 2 and monotonically approaching 1 as n increases. Thus, to undo this loading effect, the apparent σ becomes (2n + 1)M σ nm (t) = 4πa 2 (1 + k n )[ C nm(t) + i S nm (t)] (11b) This formula has been employed by, e.g., Chao et al. (1987), Chao (1994), and Wahr et al. (1998). Even though only applicable to a 2-D spherical shell, this invertibility is actually quite useful with respect to our real Earth. To be more specific, suppose we observed the TVG in terms of a set of changing Stokes coefficients. We can convert the latter readily into σ(ω) onr = a according to Eqs. (11a) or (11b) whichever is appropriate. In doing so, we construct a σ(ω) ons 0 that mimics the external field produced by the original mass redistributions, whether occurring on r = a or elsewhere. In principle that constitutes the entire information about the density redistribution that can be recovered from the Stokes coefficients alone, while nothing more can be inferred about the depth or altitude (i.e., radial distance) at which the mass redistributions actually take place. However, for the Earth, this is already a fairly accurate representation, because a majority of the important mass transports do take place near the Earth s surface.
7 B.F. Chao / Journal of Geodynamics 39 (2005) Notable departures from a spherical surface arise from the ellipticity of solid Earth s shape (less than half of 1/300a or 10 km), land topography (within a few km), atmospheric thickness (scale height about 10 km), depth of the oceans (within about 5 km), or crustal thickness (within say, 30 km). Mass transports do take place in the deeper interior of the Earth, but have typically much longer timescales which make the signal separation possible. A final point to make: suppose there is a near-surface 2-D mass redistributions, which occurs on a surface at some radial distance from r = a. How do we map the nominal density redistribution σ on r = a obtained above into its true value at its true location and depth? In the trivial case where the new surface is simply another concentric sphere, all one needs to do is to upward or downward continue σ(ω) by multiplying with the radial transfer function (a/r) n+1 just as one would for the gravity field as mentioned earlier (Eq. (2)). For an analytical departure such as the elliptical shape of the Earth, it is conceivable that some analytical expressions involving Wigner 3-j symbols can be obtained for the mapping, at least to first order of small amplitude of the radial departure. If so, then the mapping for an arbitrary bumpy surface slightly deviated from a sphere, whose amplitude can be expanded into spherical harmonics, can also be worked out. This is beyond our present interest, however. 5. Discussions and conclusions We examined the uniqueness problem of the gravitational inversion in terms of spherical harmonics. We showed that the reason for the non-uniqueness in the 3-D inversion is because the set of 2n + 1 multipoles are inadequate to recover the (n + 1)(n + 2)/2 moments for the (bodily) density function, and hence the degree of deficiency in knowledge is n(n 1)/2, which grows rapidly with the degree n. In complete contrast, in the case of a 2-D spherical shell without the radial dependence, the gravitational inversion for the (surface) density function proves to be unique. This has important implications for actual geophysical situations where the time-variable gravity signals often do originate from Earth s surface which can be approximated by a spherical shell. The simplicity in the 2-D inversion (8), (11a) and (11b) owes its existence to our use of the spherical harmonics as basis functions. The reason is traced to the fact that the solid spherical harmonics are themselves solutions to the Laplace equation which the external gravity field obeys. A consequence is the following. Since the spherical harmonics are global by definition and hence inefficient in expanding regional gravity or density variabilities, there have been efforts in adopting other types of localized basis function sets such as spatial wavelets. However, although potentially more efficient in numerical representations, these basis functions in general do not enjoy the simple relationship between the (timevariable) gravity and the density (re)distribution, as do spherical harmonics in Eqs. (8), (11a) and (11b). Hence their applicability to the gravitational inversion, particularly for the 2-D spherical shell, may not be straightforward. This of course awaits further investigations. Another advantage of treating the problem in spherical harmonics lies in the natural wavelength scale of the spherical harmonics depending on the degree n. In this paper, we have made our theoretical formulas as if we had complete knowledge of the gravity field with n going to infinity. In practice there is always a limit as to how high n can be, and one simply truncates the expansions at a certain maximum degree N that corresponds to some highest possible spatial resolution. For example, in space gravity missions, a rule of thumb is that the nominal spatial resolution 20,000/N km is comparable to the spacecraft altitude from which the spacecraft senses the gravity. Within N, the entire argument presented in the above holds
8 230 B.F. Chao / Journal of Geodynamics 39 (2005) valid perfectly. We can truncate at any N; the lack of information beyond N simply need not enter into the consideration. Acknowledgments This study is supported by NASA s Solid Earth and Natural Hazards Program. References Backus, G., Parker, R., Constable, C., Foundations of Geomagnetism. Cambridge University Press, New York. Chao, B.F., Gross, R.S., Changes in the Earth s rotation and low-degree gravitational field induced by earthquakes. Geophys. J. Roy. Astron. Soc. 91, Chao, B.F., O Connor, W.P., Chang, A.T.C., Hall, D.K., Foster, J.L., Snow-load effect on the Earth s rotation and gravitational field, J. Geophys. Res. 92, Chao, B.F., The geoid and Earth rotation. In: Vanicek, P., Christou, N. (Eds.), Geophysical Interpretations of Geoid. CRC Press, Boca Raton. Jackson, J.D., Classical Electrodynamics, 2nd ed. Wiley, New York. Kaula, W.M., Elastic models of the mantle corresponding to variations in the external gravity field. J. Geophys. Res. 68, Kaula, W.M., Theory of Satellite Geodesy. Blaisdell, Waltham, MA. Menke, W., Geophysical Data Analysis Discrete Inverse Theory. Academic Press, San Diego. Morse, P.M., Feshbach, H., Methods of Theoretical Physics. McGraw-Hill, New York. Papoulis, A., Probability, Random Variables, and Stochastic Processes. McGraw-Hill, New York. Parker, R.L., Geophysical Inverse Theory. Princeton University Press, Princeton, NJ. Rubincam, D.P., Information theory lateral density distribution for Earth inferred from global gravity field. J. Geophys. Res. 87, Tapley, B.D., Battadpur, S., Watkins, M., Reigber, C., The gravity and climate recovery experiment: mission overview and early results. Geophys. Res. Lett. 31, L09607, doi: /2004gl Wahr, J., Molenaar, M., Bryan, F., Time variability of the Earth s gravity field: Hydrological and oceanic effects and their possible detection using GRACE. J. Geophys. Res. 103, Zidarov, D., Inverse Gravimetric Problem in Geoprospecting and Geodesy. Elsevier, New York.
Copyright 2004 American Geophysical Union. Further reproduction or electronic distribution is not permitted.
Copyright 2004 American Geophysical Union. Further reproduction or electronic distribution is not permitted. Citation: Thompson, P. F., S. V. Bettadpur, and B. D. Tapley (2004), Impact of short period,
More informationImpact of short period, non-tidal, temporal mass variability on GRACE gravity estimates
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L06619, doi:10.1029/2003gl019285, 2004 Impact of short period, non-tidal, temporal mass variability on GRACE gravity estimates P. F. Thompson, S. V. Bettadpur, and
More informationis the coefficient of degree 2, order 0 of the non-dimensional spherical harmonic
Materials and Methods J is the coefficient of degree, order 0 of the non-dimensional spherical harmonic representation of the mass distribution of the Earth system. It is directly related to the diagonal
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationSatellite Gravimetry and its Application to Glaciology by Anthony Arendt for the UAF Summer School in Glaciology, June post-glacial rebound
Satellite Gravimetry and its Application to Glaciology by Anthony Arendt for the UAF Summer School in Glaciology, June 2010 1 Overview The Earth is a dynamic system in which components of the core, surface
More informationCHAPTER 3 POTENTIALS 10/13/2016. Outlines. 1. Laplace s equation. 2. The Method of Images. 3. Separation of Variables. 4. Multipole Expansion
CHAPTER 3 POTENTIALS Lee Chow Department of Physics University of Central Florida Orlando, FL 32816 Outlines 1. Laplace s equation 2. The Method of Images 3. Separation of Variables 4. Multipole Expansion
More information(ii) Observational Geomagnetism. Lecture 5: Spherical harmonic field models
(ii) Observational Geomagnetism Lecture 5: Spherical harmonic field models Lecture 5: Spherical harmonic field models 5.1 Introduction 5.2 How to represent functions on a spherical surface 5.3 Spherical
More informationNon-isotropic filtering of GRACE temporal gravity for geophysical signal enhancement
Geophys. J. Int. (2005) 163, 18 25 doi: 10.1111/j.1365-246X.2005.02756.x GJI Geodesy, potential field and applied geophysics Non-isotropic filtering of GRACE temporal gravity for geophysical signal enhancement
More informationScattering cross-section (µm 2 )
Supplementary Figures Scattering cross-section (µm 2 ).16.14.12.1.8.6.4.2 Total scattering Electric dipole, a E (1,1) Magnetic dipole, a M (1,1) Magnetic quardupole, a M (2,1). 44 48 52 56 Wavelength (nm)
More informationThe changes in normalized second degree geopotential coefficients. - f P 20 (sin4) J.), (la) /5 'GM*ft r] (lb) Ac 21 -ia* 21 = i n ^ S i f RI J^ GM.
CHAPTER 7 SOLID EARTH TIDES The solid Earth tide model is based on an abbreviated form of the Wahr model (Wahr, 98) using the Earth model 66A of Gilbert and Dziewonski (975). The Love numbers for the induced
More informationCo-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data
Chapter 2 Co-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data W.Sun,G.Fu,andSh.Okubo Abstract Dislocation theories were developed conventionally for a deformed earth
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationThe Rotational and Gravitational Signature of Recent Great Earthquakes
The Rotational and Gravitational Signature of Recent Great Earthquakes Richard S. Gross Jet Propulsion Laboratory California Institute of Technology Pasadena, CA 91109 8099, USA 7th IVS General Meeting
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2014) 199, 800 804 GJI Geodynamics and tectonics doi: 10.1093/gji/ggu301 On the gravitational energy associated with Earth s changing oblateness B. F.
More informationCHAPTER 6 SOLID EARTH TIDES
CHAPTER 6 SOLID EARTH TIDES The solid Earth tide model is based on an abbreviated form of the Wahr model (Wahr, 98) using the Earth model 66A of Gilbert and Dziewonski (975). The Love numbers for the induced
More informationGravity 3. Gravity 3. Gravitational Potential and the Geoid. Chuck Connor, Laura Connor. Potential Fields Geophysics: Week 2.
Gravitational Potential and the Geoid Chuck Connor, Laura Connor Potential Fields Geophysics: Week 2 Objectives for Week 1 Gravity as a vector Gravitational Potential The Geoid Gravity as a vector We can
More informationCHALLENGES TO THE SWARM MISSION: ON DIFFERENT INTERNAL SHA MAGNETIC FIELD MODELS OF THE EARTH IN DEPENDENCE ON SATELLITE ALTITUDES
CHALLENGES TO THE SWARM MISSION: ON DIFFERENT INTERNAL SHA MAGNETIC FIELD MODELS OF THE EARTH IN DEPENDENCE ON SATELLITE ALTITUDES Wigor A. Webers Helmholtz- Zentrum Potsdam, Deutsches GeoForschungsZentrum,
More information12 Gravity and Topography
12. Gravity and Topography 126 12 Gravity and Topography In the limit of perfect hydrostatic equilibrium, there is an exact relationship between gravity and topography... and no new information emerges:
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationCHAPTER X. Second Half Review 2017
CHAPTER X Second Half Review 217 Here is a quick overview of what we covered in the second half of the class. Remember that the final covers the whole course but there will naturally be a bias towards
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationTime-variable gravity from SLR and DORIS tracking
Time-variable gravity from SLR and DORIS tracking Frank G. Lemoine 1, Steven M. Klosko 2, Christopher M. Cox 3, Thomas J. Johnson 4 1. Planetary Geodynamics Laboratory, NASA Goddard Space Flight Center,
More informationl=0 The expansion coefficients can be determined, for example, by finding the potential on the z-axis and expanding that result in z.
Electrodynamics I Exam - Part A - Closed Book KSU 15/11/6 Name Electrodynamic Score = 14 / 14 points Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try
More informationPhysics 342 Lecture 4. Probability Density. Lecture 4. Physics 342 Quantum Mechanics I
Physics 342 Lecture 4 Probability Density Lecture 4 Physics 342 Quantum Mechanics I Monday, February 4th, 28 We review the basic notions of probability, in particular the role of probability density in
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationVortex stretching in incompressible and compressible fluids
Vortex stretching in incompressible and compressible fluids Esteban G. Tabak, Fluid Dynamics II, Spring 00 1 Introduction The primitive form of the incompressible Euler equations is given by ( ) du P =
More informationMultipole moments. November 9, 2015
Multipole moments November 9, 5 The far field expansion Suppose we have a localized charge distribution, confined to a region near the origin with r < R. Then for values of r > R, the electric field must
More informationConsistent Histories. Chapter Chain Operators and Weights
Chapter 10 Consistent Histories 10.1 Chain Operators and Weights The previous chapter showed how the Born rule can be used to assign probabilities to a sample space of histories based upon an initial state
More informationGeophysics Departmental Exam: 2004 Part 1
2004 Geophysics Departmental Exam: 2004 Part 1 This section is 90 minutes, closed book, and consists of questions designed to test your knowledge of facts and figures in the geosciences. The focus will
More informationTidal Effects on Earth s Surface
Tidal Effects on Earth s Surface Tom Murphy February, 1 This treatment follows the conventions of F. D. Stacey s Physics of the Earth, and is largely an elaboration on this work. 1 Tidal Potential The
More information221B Lecture Notes Scattering Theory II
22B Lecture Notes Scattering Theory II Born Approximation Lippmann Schwinger equation ψ = φ + V ψ, () E H 0 + iɛ is an exact equation for the scattering problem, but it still is an equation to be solved
More information2 Voltage Potential Due to an Arbitrary Charge Distribution
Solution to the Static Charge Distribution on a Thin Wire Using the Method of Moments James R Nagel Department of Electrical and Computer Engineering University of Utah, Salt Lake City, Utah April 2, 202
More informationSupporting Information
Supporting Information Harig and Simons 1.173/pnas.178519 SI Text Determination of Noise. Gravity Recovery and Climate Experiment (GRACE) data are released as spherical harmonic coefficients along with
More informationShort Wire Antennas: A Simplified Approach Part I: Scaling Arguments. Dan Dobkin version 1.0 July 8, 2005
Short Wire Antennas: A Simplified Approach Part I: Scaling Arguments Dan Dobkin version 1.0 July 8, 2005 0. Introduction: How does a wire dipole antenna work? How do we find the resistance and the reactance?
More informationPhysics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 374 (2010) 1063 1067 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Macroscopic far-field observation of the sub-wavelength near-field dipole
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationIncompatibility Paradoxes
Chapter 22 Incompatibility Paradoxes 22.1 Simultaneous Values There is never any difficulty in supposing that a classical mechanical system possesses, at a particular instant of time, precise values of
More informationJackson 6.4 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 6.4 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell PROBLEM: A uniformly magnetized and conducting sphere of radius R and total magnetic moment m = 4πMR 3
More informationSpherical Harmonics and Related Topics. David Randall 2 S = 0, r 2 r r S 2. S = r n Y n
! Revised April 10, 2017 1:32 PM! 1 Spherical Harmonics and Related Topics David Randall The spherical surface harmonics are convenient functions for representing the distribution of geophysical quantities
More information14. Rotational Kinematics and Moment of Inertia
14. Rotational Kinematics and Moment of nertia A) Overview n this unit we will introduce rotational motion. n particular, we will introduce the angular kinematic variables that are used to describe the
More informationThe Basic Properties of Surface Waves
The Basic Properties of Surface Waves Lapo Boschi lapo@erdw.ethz.ch April 24, 202 Love and Rayleigh Waves Whenever an elastic medium is bounded by a free surface, coherent waves arise that travel along
More information3 Chapter. Gauss s Law
3 Chapter Gauss s Law 3.1 Electric Flux... 3-2 3.2 Gauss s Law (see also Gauss s Law Simulation in Section 3.10)... 3-4 Example 3.1: Infinitely Long Rod of Uniform Charge Density... 3-9 Example 3.2: Infinite
More informationr( θ) = cos2 θ ω rotation rate θ g geographic latitude - - θ geocentric latitude - - Reference Earth Model - WGS84 (Copyright 2002, David T.
1 Reference Earth Model - WGS84 (Copyright 22, David T. Sandwell) ω spheroid c θ θ g a parameter description formula value/unit GM e (WGS84) 3.9864418 x 1 14 m 3 s 2 M e mass of earth - 5.98 x 1 24 kg
More informationChapter 15 Appendix Moment of Inertia of a Spherical Shell
Chapter 1 Appendix Moment of nertia of a Spherical Shell t is common to regard the rotation of a rigid object with spherical symmetry; we live on one such object (not entirely uniform and not entirely
More informationPhysics 221A Fall 1996 Notes 12 Orbital Angular Momentum and Spherical Harmonics
Physics 221A Fall 1996 Notes 12 Orbital Angular Momentum and Spherical Harmonics We now consider the spatial degrees of freedom of a particle moving in 3-dimensional space, which of course is an important
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
More informationNumerical Investigation on Spherical Harmonic Synthesis and Analysis
Numerical Investigation on Spherical Harmonic Synthesis and Analysis Johnny Bärlund Master of Science Thesis in Geodesy No. 3137 TRITA-GIT EX 15-006 School of Architecture and the Built Environment Royal
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationSolutions to Laplace s Equations- II
Solutions to Laplace s Equations- II Lecture 15: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Laplace s Equation in Spherical Coordinates : In spherical coordinates
More informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationA Mission to Planet Mars Gravity Field Determination
A Mission to Planet Mars Gravity Field Determination Department for Theoretical Geodesy Graz University of Technology and Space Research Institute Austrian Academy of Sciences Gravity field CHAMP GRACE
More informationInformation on internal structure from shape, gravity field and rotation
Information on internal structure from shape, gravity field and rotation Seismological information is available only for the Earth and in limited amounts for the Moon. Various geodetic data put constraints
More informationFundamentals on light scattering, absorption and thermal radiation, and its relation to the vector radiative transfer equation
Fundamentals on light scattering, absorption and thermal radiation, and its relation to the vector radiative transfer equation Klaus Jockers November 11, 2014 Max-Planck-Institut für Sonnensystemforschung
More informationChecking Consistency. Chapter Introduction Support of a Consistent Family
Chapter 11 Checking Consistency 11.1 Introduction The conditions which define a consistent family of histories were stated in Ch. 10. The sample space must consist of a collection of mutually orthogonal
More informationSolutions: Homework 2
Solutions: Homework 2 Ex 2.1: Particle moving in Magnetic Field (a) If the initial velocity is perpendicular to the magnetic field B, the magnitude of the Lorentz Force (Eq 1.52) is F = 1 c q qu B ub (1)
More information1 Potential due to a charged wire/sheet
Lecture XXX Renormalization, Regularization and Electrostatics Let us calculate the potential due to an infinitely large object, e.g. a uniformly charged wire or a uniformly charged sheet. Our main interest
More informationQuasi-geostrophic ocean models
Quasi-geostrophic ocean models March 19, 2002 1 Introduction The starting point for theoretical and numerical study of the three dimensional large-scale circulation of the atmosphere and ocean is a vorticity
More informationBohr & Wheeler Fission Theory Calculation 4 March 2009
Bohr & Wheeler Fission Theory Calculation 4 March 9 () Introduction The goal here is to reproduce the calculation of the limiting Z /A against spontaneous fission Z A lim a S. (.) a C as first done by
More informationGeomagnetic Field Modeling Lessons learned from Ørsted and CHAMP and prospects for Swarm
Geomagnetic Field Modeling Lessons learned from Ørsted and CHAMP and prospects for Swarm Nils Olsen RAS Discussion Meeting on Swarm October 9 th 2009 Nils Olsen (DTU Space) Ørsted, CHAMP, and Swarm 1 /
More informationGravity Recovery and Climate Experiment (GRACE) alias error from ocean tides
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2006jb004747, 2008 Gravity Recovery and Climate Experiment (GRACE) alias error from ocean tides K. W. Seo, 1,2 C. R. Wilson, 3 S. C. Han, 4,5 and
More informationLecture 6 - Introduction to Electricity
Lecture 6 - Introduction to Electricity A Puzzle... We are all familiar with visualizing an integral as the area under a curve. For example, a b f[x] dx equals the sum of the areas of the rectangles of
More information### dv where ρ is density, R is distance from rotation axis, and dv is
Comments This is one of my favorite problem assignments in our Solid Earth Geophysics class, typically taken by junior and senior concentrators and by first-year graduate students. I encourage students
More informationCaltech Ph106 Fall 2001
Caltech h106 Fall 2001 ath for physicists: differential forms Disclaimer: this is a first draft, so a few signs might be off. 1 Basic properties Differential forms come up in various parts of theoretical
More informationMean Vertical Gradient of Gravity
Mean Vertical Gradient of Gravity P. Vaníek, J. Janák Department of Geodesy and Geomatics Engineering, University of New Brunswick, P.O.Box 4400, Fredericton, New Brunswick, Canada, E3B 2M8 J. Huang Geodetic
More information10 Least-squares collocation
10 Least-squares collocation 10.1 Principles of least-squares collocation The principle of collocation is very simple. The anomalous potential T outside the earth is a harmonic function, that is, it satisfies
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationPhysics 221A Fall 1996 Notes 21 Hyperfine Structure in Hydrogen and Alkali Atoms
Physics 221A Fall 1996 Notes 21 Hyperfine Structure in Hydrogen and Alkali Atoms Hyperfine effects in atomic physics are due to the interaction of the atomic electrons with the electric and magnetic multipole
More informationAn Axiomatic Approach to Quantum Mechanics
Available online at www.worldscientificnews.com WSN 116 (019) 09-1 EISSN 39-19 An Axiomatic Approach to Quantum Mechanics M. A. I. S. Rajakaruna 1,a, K. A. I. L. Wijewardena Gamalath,b 1 178/A, Sri Wimalarama
More informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More informationSolutions to PS 2 Physics 201
Solutions to PS Physics 1 1. ke dq E = i (1) r = i = i k eλ = i k eλ = i k eλ k e λ xdx () (x x) (x x )dx (x x ) + x dx () (x x ) x ln + x x + x x (4) x + x ln + x (5) x + x To find the field for x, we
More informationStochastic Processes
qmc082.tex. Version of 30 September 2010. Lecture Notes on Quantum Mechanics No. 8 R. B. Griffiths References: Stochastic Processes CQT = R. B. Griffiths, Consistent Quantum Theory (Cambridge, 2002) DeGroot
More informationElectrodynamics Exam Solutions
Electrodynamics Exam Solutions Name: FS 215 Prof. C. Anastasiou Student number: Exercise 1 2 3 4 Total Max. points 15 15 15 15 6 Points Visum 1 Visum 2 The exam lasts 18 minutes. Start every new exercise
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationGravity variations induced by core flows
Geophys. J. Int. (2010) 180, 635 650 doi: 10.1111/j.1365-246X.2009.04437.x Gravity variations induced by core flows Mathieu Dumberry Department of Physics, University of Alberta, Edmonton, Canada T6G 2G7.
More informationPrinciples of the Global Positioning System Lecture 18" Mathematical models in GPS" Mathematical models used in GPS"
12.540 Principles of the Global Positioning System Lecture 18" Prof. Thomas Herring" Room 54-820A; 253-5941" tah@mit.edu" http://geoweb.mit.edu/~tah/12.540 " Mathematical models in GPS" Review assignment
More informationFoundations of Geomagnetism
Foundations of Geomagnetism GEORGE BACKUS University of California, San Diego ROBERT PARKER University of California, San Diego CATHERINE CONSTABLE University of California, San Diego m.m CAMBRIDGE UNIVERSITY
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationStability of Thick Spherical Shells
Continuum Mech. Thermodyn. (1995) 7: 249-258 Stability of Thick Spherical Shells I-Shih Liu 1 Instituto de Matemática, Universidade Federal do Rio de Janeiro Caixa Postal 68530, Rio de Janeiro 21945-970,
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationCosmic Variance of the Three-Point Correlation Function of the Cosmic Microwave Background
CfPA 93 th 18 astro-ph/9306012 June 1993 REVISED arxiv:astro-ph/9306012v2 14 Jul 1993 Cosmic Variance of the Three-Point Correlation Function of the Cosmic Microwave Background Mark Srednicki* Center for
More informationWhat is the force of gravity on (pointlike) M 2 caused by (pointlike) M 1? M 2 M 1 = +G M 1M 2. ˆ r. F grav. r 2 = G M 1 M 2. r 3
GRAVITY 2-1 GRAVITY - Newton s law of gravity - Gravitational force (direct integration) - gravitational field - Symmetry and Invariance arguments - gravitational potential, and PE 2-2 M 1 r A) B) C) D)
More informationConsistency of Earth Rotation, Gravity, and Shape Measurements
Consistency of Earth Rotation, Gravity, and Shape Measurements Richard S. Gross*, David A. Lavallée, Geoffrey Blewitt, and Peter J. Clarke *Jet Propulsion Laboratory, California Institute of Technology
More informationShallow Water Gravity Waves: A Note on the Particle Orbits
Journal of Oceanography Vol. 5, pp. 353 to 357. 1996 Shallow Water Gravity Waves: A Note on the Particle Orbits KERN E. KENYON 463 North Lane, Del Mar, CA 9014-4134, U.S.A. (Received 4 July 1995; in revised
More informationA wavenumber approach to characterizing the diffuse field conditions in reverberation rooms
PROCEEDINGS of the 22 nd International Congress on Acoustics Isotropy and Diffuseness in Room Acoustics: Paper ICA2016-578 A wavenumber approach to characterizing the diffuse field conditions in reverberation
More informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
More informationStructural Dynamics Prof. P. Banerji Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 1 Introduction
Structural Dynamics Prof. P. Banerji Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 1 Introduction Hello, I am Pradipta Banerji from the department of civil engineering,
More informationOcean Bottom Pressure Changes Lead to a Decreasing Length-of-Day in a Warming Climate
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:1.129/, Ocean Bottom Pressure Changes Lead to a Decreasing Length-of-Day in a Warming Climate Felix W. Landerer, Johann H. Jungclaus, Jochem Marotzke Max
More informationHOMEWORK SET 1 SOLUTIONS MATH 456, Spring 2011 Bruce Turkington
HOMEWORK SET SOLUTIONS MATH 456, Spring Bruce Turkington. Consider a roll of paper, like a toilet roll. Its cross section is very nearly an annulus of inner radius R and outer radius R. The thickness of
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
More information/. 3D. a very abstract sort of problem, this research is a step on the way to applicable knowledge.
/. 3D Corwin Atwood-Stone Physics Department, The College of Wooster, Wooster, Ohio 691, USA (Dated: May 7, 009) A simulation of the mutual orbital motion of a point and a stick of equal masses in three
More informationhas a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity.
http://preposterousuniverse.com/grnotes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been framed
More informationMathematics that Every Physicist should Know: Scalar, Vector, and Tensor Fields in the Space of Real n- Dimensional Independent Variable with Metric
Mathematics that Every Physicist should Know: Scalar, Vector, and Tensor Fields in the Space of Real n- Dimensional Independent Variable with Metric By Y. N. Keilman AltSci@basicisp.net Every physicist
More informationPage 52. Lecture 3: Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 2008/10/03 Date Given: 2008/10/03
Page 5 Lecture : Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 008/10/0 Date Given: 008/10/0 Inner Product Spaces: Definitions Section. Mathematical Preliminaries: Inner
More informationarxiv: v1 [gr-qc] 26 Jun 2008
Chaotic Motion Around Stellar Objects with Octupolar Deformation: Newtonian and Post Newtonian Approaches Javier Ramos-Caro, Framsol López-Suspes, and Guillermo A. González Grupo de Investigación en Relatividad
More informationChapter 4 Multipole model of the Earth's magnetic field
Chapter 4 Multipole model of the Earth's magnetic field 1 Previously A measurement of the geomagnetic field at any given point and time consists of a superposition of fields from different sources: Internal
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationDielectrics. Lecture 20: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay
What are dielectrics? Dielectrics Lecture 20: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay So far we have been discussing electrostatics in either vacuum or in a conductor.
More information