On inversion for mass distribution from global (time-variable) gravity field

Size: px
Start display at page:

Download "On inversion for mass distribution from global (time-variable) gravity field"

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. 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 information

Impact of short period, non-tidal, temporal mass variability on GRACE gravity estimates

Impact 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 information

is the coefficient of degree 2, order 0 of the non-dimensional spherical harmonic

is 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 information

Physics 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 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 information

Satellite 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 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 information

CHAPTER 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 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 (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 information

Non-isotropic filtering of GRACE temporal gravity for geophysical signal enhancement

Non-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 information

Scattering cross-section (µm 2 )

Scattering 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 information

The 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.

The 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 information

Co-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data

Co-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 information

Navier-Stokes Equation: Principle of Conservation of Momentum

Navier-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 information

The Rotational and Gravitational Signature of Recent Great Earthquakes

The 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 information

Geophysical Journal International

Geophysical 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 information

CHAPTER 6 SOLID EARTH TIDES

CHAPTER 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 information

Gravity 3. Gravity 3. Gravitational Potential and the Geoid. Chuck Connor, Laura Connor. Potential Fields Geophysics: Week 2.

Gravity 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 information

CHALLENGES 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 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 information

12 Gravity and Topography

12 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 information

Notes: 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.

Notes: 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 information

CHAPTER X. Second Half Review 2017

CHAPTER 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 information

d 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.

d 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 information

Time-variable gravity from SLR and DORIS tracking

Time-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 information

l=0 The expansion coefficients can be determined, for example, by finding the potential on the z-axis and expanding that result in z.

l=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 information

Physics 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 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 information

Practical in Numerical Astronomy, SS 2012 LECTURE 9

Practical 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 information

Classical Field Theory: Electrostatics-Magnetostatics

Classical 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 information

Vortex stretching in incompressible and compressible fluids

Vortex 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 information

Multipole moments. November 9, 2015

Multipole 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 information

Consistent Histories. Chapter Chain Operators and Weights

Consistent 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 information

Geophysics Departmental Exam: 2004 Part 1

Geophysics 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 information

Tidal Effects on Earth s Surface

Tidal 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 information

221B Lecture Notes Scattering Theory II

221B 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 information

2 Voltage Potential Due to an Arbitrary Charge Distribution

2 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 information

Supporting Information

Supporting 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 information

Short 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 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 information

Physics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.

Physics 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 information

4 Power Series Solutions: Frobenius Method

4 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 information

Incompatibility Paradoxes

Incompatibility 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 information

Jackson 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 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 information

Spherical Harmonics and Related Topics. David Randall 2 S = 0, r 2 r r S 2. S = r n Y n

Spherical 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 information

14. Rotational Kinematics and Moment of Inertia

14. 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 information

The Basic Properties of Surface Waves

The 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 information

3 Chapter. Gauss s Law

3 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 information

r( θ) = cos2 θ ω rotation rate θ g geographic latitude - - θ geocentric latitude - - Reference Earth Model - WGS84 (Copyright 2002, David T.

r( θ) = 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 information

Chapter 15 Appendix Moment of Inertia of a Spherical Shell

Chapter 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 information

Physics 221A Fall 1996 Notes 12 Orbital Angular Momentum and Spherical Harmonics

Physics 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 information

Introduction to Group Theory

Introduction 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 information

Numerical Investigation on Spherical Harmonic Synthesis and Analysis

Numerical 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 information

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature

Chapter 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 information

Solutions to Laplace s Equations- II

Solutions 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 information

Expansion of 1/r potential in Legendre polynomials

Expansion 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 information

A Mission to Planet Mars Gravity Field Determination

A 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 information

Information on internal structure from shape, gravity field and rotation

Information 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 information

Fundamentals 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 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 information

Checking Consistency. Chapter Introduction Support of a Consistent Family

Checking 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 information

Solutions: Homework 2

Solutions: 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 information

1 Potential due to a charged wire/sheet

1 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 information

Quasi-geostrophic ocean models

Quasi-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 information

Bohr & Wheeler Fission Theory Calculation 4 March 2009

Bohr & 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 information

Geomagnetic 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 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 information

Gravity Recovery and Climate Experiment (GRACE) alias error from ocean tides

Gravity 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 information

Lecture 6 - Introduction to Electricity

Lecture 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

### 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 information

Caltech Ph106 Fall 2001

Caltech 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 information

Mean Vertical Gradient of Gravity

Mean 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 information

10 Least-squares collocation

10 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 information

Physics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top

Physics 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 information

Physics 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 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 information

An Axiomatic Approach to Quantum Mechanics

An 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 information

Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations

Lecture 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 information

Solutions to PS 2 Physics 201

Solutions 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 information

Stochastic Processes

Stochastic 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 information

Electrodynamics Exam Solutions

Electrodynamics 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 information

1 Exercise: Linear, incompressible Stokes flow with FE

1 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 information

Gravity variations induced by core flows

Gravity 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 information

Principles of the Global Positioning System Lecture 18" Mathematical models in GPS" Mathematical models used in GPS"

Principles 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 information

Foundations of Geomagnetism

Foundations 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 information

Control Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017

Control 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 information

Lecture 1: Introduction and Review

Lecture 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 information

Stability of Thick Spherical Shells

Stability 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 information

The 3 dimensional Schrödinger Equation

The 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 information

Cosmic Variance of the Three-Point Correlation Function of the Cosmic Microwave Background

Cosmic 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 information

What 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

What 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 information

Consistency of Earth Rotation, Gravity, and Shape Measurements

Consistency 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 information

Shallow Water Gravity Waves: A Note on the Particle Orbits

Shallow 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 information

A wavenumber approach to characterizing the diffuse field conditions in reverberation rooms

A 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 information

Elastic 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. 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 information

Structural 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 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 information

Ocean Bottom Pressure Changes Lead to a Decreasing Length-of-Day in a Warming Climate

Ocean 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 information

HOMEWORK SET 1 SOLUTIONS MATH 456, Spring 2011 Bruce Turkington

HOMEWORK 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 information

1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018

1 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. 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 information

has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity.

has 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 information

Mathematics 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 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 information

Page 52. Lecture 3: Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 2008/10/03 Date Given: 2008/10/03

Page 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 information

arxiv: v1 [gr-qc] 26 Jun 2008

arxiv: 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 information

Chapter 4 Multipole model of the Earth's magnetic field

Chapter 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 information

Quantum Mechanics in Three Dimensions

Quantum 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 information

Chapter 1. Governing Equations of GFD. 1.1 Mass continuity

Chapter 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 information

Dielectrics. Lecture 20: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay

Dielectrics. 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