Adams–Williamson equation
The Adams–Williamson equation, named after L. H. Adams and E. D. Williamson, is a relation between the velocities of seismic waves and the density of the Earth's interior. Given the average density of rocks at the Earth's surface and profiles of the Pwave and Swave speeds as function of depth, it can predict how density increases with depth. It assumes that the compression is adiabatic and that the Earth is spherically symmetric, homogeneous, and in hydrostatic equilibrium. It can also be applied to spherical shells with that property. It is an important part of models of the Earth's interior such as the Preliminary Reference Earth Model (PREM).^{[1]}^{[2]}
Contents
History
Williamson and Adams first developed the theory in 1923. They concluded that "It is therefore impossible to explain the high density of the Earth on the basis of compression alone. The dense interior cannot consist of ordinary rocks compressed to a small volume; we must therefore fall back on the only reasonable alternative, namely, the presence of a heavier material, presumably some metal, which, to judge from its abundance in the Earth's crust, in meteorites and in the Sun, is probably iron."^{[1]}
Theory
The two types of seismic body waves are compressional waves (Pwaves) and shear waves (Swaves). Both have speeds that are determined by the elastic properties of the medium they travel through, in particular the bulk modulus K, the shear modulus μ, and the density ρ. In terms of these parameters, the Pwave speed v_{p} and the Swave speed v_{s} are
\[ \begin{align} v_p &= \sqrt{\frac{K+(4/3)\mu}{\rho}} \\ v_s &= \sqrt{\frac{\mu}{\rho}}. \end{align}\]
These two speeds can be combined in a seismic parameter
{{{Script error}}}

</p> 
({{{3}}}) 
The definition of the bulk modulus,
\[K = V\frac{dP}{dV},\]
is equivalent to {{{Script error}}}

</p> 
({{{3}}}) 
Suppose a region at a distance r from the Earth's center can be considered a fluid in hydrostatic equilibrium, it is acted on by gravitational attraction from the part of the Earth that is below it and pressure from the part above it. Also suppose that the compression is adiabatic (so thermal expansion does not contribute to density variations). The pressure P(r) varies with r as
{{{Script error}}}

</p> 
({{{3}}}) 
where g(r) is the gravitational acceleration at radius r.^{[1]}
If Equations 1,2 and 3 are combined, we get the Adams–Williamson equation:
\[ \frac{d\rho}{dr} = \frac{\rho(r)g(r)}{\Phi(r)}.\]
This equation can be integrated to obtain
\[ \ln\left(\frac{\rho}{\rho_0}\right) = \int_{r_0}^r \frac{g(r)}{\Phi(r)}dr, \]
where r_{0} is the radius at the Earth's surface and ρ_{0} is the density at the surface. Given ρ_{0} and profiles of the P and Swave speeds, the radial dependence of the density can be determined by numerical integration.^{[1]}
References
 ↑ ^{1.0} ^{1.1} ^{1.2} ^{1.3} Poirier 2000
 ↑ Dziewonski & Anderson 1981
References
 Script error
 Script error