Geophysics: Calculating Inner Core Radius

In summary: Geometrical optics approach is sufficient. So we don't have to talk about diffractions of the wave, just refractions at the interfaces.This is correct. The interfaces are simply points where the two media meet. The wave does not diffract at these points, it just refracts. 3) I'm going to use terms from geometrical optics, so I'll call the path a P-wave took to get from where it originated to where it was detected a "beam".This is correct. The terms you are using are related to geometrical optics.
  • #1
s7b
26
0
I have a question about some geophysics...


When we want to calculate the radius of an inner core how do you do it. If I were to be given the radius of the whole planet as well as density and the velocity of the Pwave, where do I begin to solve. Also let's say the angular distances where P-waves are observed are also given.
 
Physics news on Phys.org
  • #2
If you have the P-wave velocity profile, you will see the velocity gradually increases with a few little jumps in it, then it suddenly drops down dramatically and then steadily begins to increase, and then there is another giant leap up in the velocity. That last giant leap up is the velocity going from the outer core to the inner core -- the radius of the inner core is radius of the Earth minus the depth of this giant leap.

Incidentally, if you didnt have the velocity profile, but you had observations of P-waves from all over the world, you could invert the traveltime data to find the velocity profile which would then show you the inner core depth in the way I described above.

You might also want to look at the moment of inertia of Earth, as that constrains density distribution which might be useful in finding inner core radius.
 
  • #3
The actual wave solution is incredibly complex, but I think you can get a descent estimate based on the following assumptions.

1) There is only one sharp change in density, and that defines the core. (The crust is not going to make much of a difference for this approach, as it is sufficiently thin.)

2) Geometrical optics approach is sufficient. So we don't have to talk about diffractions of the wave, just refractions at the interfaces.

I'm going to use terms from geometrical optics, so I'll call the path a P-wave took to get from where it originated to where it was detected a "beam".

The problem is cylindrically symmetric, so we can consider a 2D problem only. So picture Earth as two concentric circles with radii R1<R2. R2 is Earth's radius, and R1 is what you are trying to find. Inner circle is the core with some density giving speed of P-wave propagation c1, and the outside of that circle the P-wave propagates at c2. I will denote the center of the two circles as point O.

Consider a beam that leaves point A on the outer circle and hits the core at point B. Call angle OAB θ2 and AOB θ1. Distance from A to B is d. Then:

[tex]R_1 sin(\theta_1) = d sin(\theta_2)[/tex]

[tex]R_1 cos(\theta_1) + d cos(\theta_2) = R1 + R2[/tex]

After entering the core, the beam is refracted and exits the core at point C. Refraction will depend on wave speed propagation. If we call angle of incidence α and angle of refraction β, we can define the following two relations.

[tex]\alpha = \theta_1 + \theta_2[/tex]

[tex]\frac{sin(\beta)}{c_1} = \frac{sin(\alpha)}{c_2}[/tex]

First is purely geometric, second is due to Snell's Law of refraction, which holds for P-waves in this approximation.

After exiting core at C, the beam continues to point D on outer circle, where it is detected. Note that angle OBC is the angle of refraction β, and it is equal to angle OCB. Angle BOC is then 180°-2β. That gives us angle AOD = AOB + BOC + COD = BOC + 2*AOD. And that's the angle from epicenter to the point where the wave was detected. Let me denote it as γ.

[tex]\gamma = 180° - 2 (\beta - \theta_1)[/tex]

So you can find point of exit as function of θ1. Though, not a trivial one. Now you need the time. You already know that the beam will travel distance d, defined by an earlier equation, before it hits the core. It will travel the same distance after leaving the core. What's left is to compute distance within the core. Call that distance b. Again, using triangle BOC for reference, we can see that the following holds.

[tex]b = 2 R_1 cos(\beta)[/tex]

And this finally gives us time.

[tex]t = 2 \frac{d}{c_2} + \frac{b}{c_1}[/tex]

Now you can take your data and fit it to the unknown R1 and θ1 on known γ and t. It's going to be an ugly fit, but presumably you'd have a computer program do it, so there shouldn't be any problems. If you have good data from opposite side of the planet with precise times, you should be able to get a very good estimate for the radius.
 
  • #4
True that it is a complex problem. In fact it is more complex than you make out.

K^2 said:
1) There is only one sharp change in density, and that defines the core. (The crust is not going to make much of a difference for this approach, as it is sufficiently thin.)

There are two profound interfaces. The first you correctly identified is the core/mantle boundary; the second is the outer/inner core boundary, which I believe is the boundary the questioner is interested in.

2) Geometrical optics approach is sufficient. So we don't have to talk about diffractions of the wave, just refractions at the interfaces.
Yes, geometrical ray theory gives a pretty good approximation to the true finite-frequency wavefield, which is useful in determining a useful background velocity model of the Earth.

Inner circle is the core with some density giving speed of P-wave propagation c1, and the outside of that circle the P-wave propagates at c2.
This is much too simple. In fact the velocity increases smoothly with depth in the Earth, this will have profound implications to the paths that rays take, and hence the observations you make at the surface.
 
  • #5
Well, generalizing to two interfaces is simple enough, but if you introduce variable density, you are pretty much stuck doing finite element analysis on the thing. You can still get away with a 2D case, and if we assume that density varies with radius only, this can be solved easily enough, but it's a purely numerical problem now. You make a guess for the density, construct the response, and then start adjusting the density function until you hit the correct response.
 

1. What is geophysics?

Geophysics is a branch of science that studies the physical properties and processes of the Earth, including its structure, composition, and dynamics.

2. How is the inner core radius calculated?

The inner core radius is calculated using a combination of seismic data and mathematical models. Seismic waves travel at different speeds through different materials, allowing scientists to determine the boundaries between the layers of the Earth. By studying these waves, scientists can estimate the radius and composition of the inner core.

3. Why is calculating the inner core radius important?

Calculating the inner core radius is important for understanding the structure and evolution of the Earth. It can also provide valuable information for predicting earthquakes and other geological events.

4. How accurate are the calculations of the inner core radius?

The calculations of the inner core radius are constantly being refined and improved, but they can still be affected by uncertainties and limitations in data and models. The current estimated radius of the inner core is around 1,220 kilometers with an uncertainty of about ±10 kilometers.

5. Are there any factors that can affect the accuracy of the inner core radius calculations?

Yes, there are several factors that can affect the accuracy of the inner core radius calculations, including the quality and quantity of seismic data, the complexity of the Earth's structure, and the assumptions and limitations of the mathematical models used. Ongoing research and advancements in technology and techniques are helping to improve the accuracy of these calculations.

Similar threads

  • Other Physics Topics
Replies
9
Views
2K
Replies
1
Views
828
  • Electromagnetism
2
Replies
42
Views
706
  • Introductory Physics Homework Help
Replies
2
Views
823
  • Introductory Physics Homework Help
Replies
6
Views
967
  • Astronomy and Astrophysics
Replies
3
Views
1K
  • Other Physics Topics
Replies
2
Views
1K
  • Special and General Relativity
Replies
21
Views
709
Replies
2
Views
1K
Back
Top