# Birch's law

Birch's law, discovered by the geophysicist Francis Birch, establishes a linear relation between compressional wave velocity v$_p$ and density $\rho$ of rocks and minerals:
    $v_{\mathrm {p} }=a({\bar {M}})+b\rho$
where $\,{\bar {M}}\,$ is the mean atomic mass in formula units and $\,a(x)\,$ is an empirical function determined by experiment.

## Example

The mean atomic mass of forsterite (Mg$_2$SiO$_4$) is equal to the sum of the atomic masses divided by the number of atoms in the formula:
    ${\bar {M}}={\frac {2\times (24.3)\;+\;1\times (28.3)\;+\;4\times (16.0)}{2\qquad \;+\;\qquad 1\qquad \,+\;\qquad 4}}=20.13~.$
Typical oxides and silicates in the mantle have values close to 20, while in the Earth's core it is close to 50.

## Applications

Birch's law applies to rocks that are under pressures of a few tens of gigapascals, enough for most cracks to close. It can be used in the discussion of geophysical data. The law is used in forming compositional and mineralogical models of the mantle by using the change in the velocity of the seismic wave and its relationship with a change in density of the material the wave is moving in. Birch's law is used in determining chemical similarities in the mantle as well as the discontinuities of the transition zones. Birch's Law can also be employed in the calculation of an increase of velocity due to an increase in the density of material.

## Shortcomings

It had been previously assumed that the velocity-density relationship is constant. That is, that Birch's law will hold true in any case, but the relationship does not hold true deeper in the mantle for the increased pressures near the transition zone. In cases where Birch's law was applied beyond the transition zone, parts of the formula need to be revised. For higher pressure regimes, different laws may be needed to determine wave velocities.

## Determining Birch's law experimentally

The relationship between the density of a material and the velocity of a P wave moving through the material was noted when research was conducted on waves in different materials.
In the experiment, a pulse of voltage is applied to a circular plate of polarized barium titanate ceramic (the transducer) which is attached to the near end of the material sample. The added voltage creates vibrations in the sample. Those vibrations travel through the sample to a second transducer on the far end. The vibrations are then converted into an electrical wave which is viewed on an oscilloscope to determine the travel time. The velocity is the lender of the damper decided by the wave's travel time.
The resulting relationship between the density of the material and the discovered velocity is known as Birch's law.

## Velocity of compressional waves in rocks

The table shows the velocities for different rocks ranging in pressure from 10 bars to 10,000 bars. It represents how the change in density, as given in the second column, is related to the velocity of the P wave moving in the material. An increase in the density of the material leads to an increase in the velocity which can be determined using Birch's Law.
    Velocities of Compressional Waves in Rocks Rock Type Rock Location Rock Density Wave Velocity (in km/s) at Pressure: 10 bar500 bar2000 bar10,000 bar SerpentiniteThetford, Quebec2.6015.6—5.736.00 SerpentiniteLudlow, VT2.6144.76.336.596.82 Granite, “G.I.”Westerly, RI2.6194.15.635.976.23 GraniteQuincy, MA2.6215.16.046.206.45 GraniteRockport, MA2.6245.05.966.296.51 GraniteStone Mt., GA2.6253.75.426.166.40 GraniteChelmsford, MA2.6264.25.646.096.35 GneissPelham, MA2.6433.45.676.066.31 Quartz monzonitePorterville, CA2.6445.1—6.076.37 QuartziteMT2.6475.6—6.156.35 GraniteHyderabad, India2.6545.46.266.386.56 GraniteBarre, VT2.6555.15.866.156.39 SandstoneNY2.6593.95.05.445.85 Pyrophyllite graniteSacred Heart, MN2.6625.9—6.286.45 GraniteBarriefield, Ontario2.6725.76.216.356.51 GneissHell Gate, NY2.6755.16.066.236.50 GraniteHyderabad, India2.6765.7—6.466.61 “Granite”Englehart, Ontario2.6796.16.286.376.57 GreywackeNew Zealand2.6795.45.635.876.13 “Granite”Larchford, Ontario2.6835.76.136.256.41 AlbiteSylmar, PA2.6876.40—6.656.76 GranodioriteButte, MT2.7054.4—6.356.56 GraywackeQuebec2.7055.4—6.046.28 SerpentiniteCA2.7105.8—6.086.31 SlateMedford, MA2.7345.49—5.916.22 “Charnockite”Pallavaram, India2.7406.15—6.306.46 Granodiorite gneissNH2.7584.4—6.076.30 TonaliteVal Verde, CA2.7635.1—6.436.60 AnorthositeTahawus, NY2.7686.73—6.907.02 AnorthositeStillwater Complex, MT2.7706.5—7.017.10 Augite syeniteOntario2.7805.7—6.636.79 Mica schistWoodsville, VT2.7975.7—6.486.64 SerpentiniteLudlow, VT2.7986.4—6.576.84 Quartz dioriteSan Luis Rey quad., CA2.7985.1—6.526.71 AnorthositeBushveld Complex2.8075.76.927.057.21 Chlorite schistChester Quarry, VT2.8414.8—6.827.07 Quartz dioriteDedham, MA2.9065.5—6.536.71 Talc schistChester, VT2.9144.9—6.506.97 GabbroMellen, WI2.9316.87.047.097.21 DiabaseCenterville, VA2.9766.14—6.766.93 DiabaseHolyoke, MA2.9776.256.406.476.63 NoritePertoria, Transvaal2.9786.67.027.117.28 DuniteWebster, NC2.9806.0—6.466.79 DiabaseSudbury, Ontario3.0036.46.676.766.91 DiabaseFrederick, MD3.0126.76—6.806.92 GabbroFrench Creek, PA3.0545.86.747.027.23 AmphiboliteMadison Co., MT3.1206.89—7.127.35 JadeiteJapan3.1807.6—8.228.28 Actinoliter schistChester, VT3.1946.61—7.207.54 DuniteWebster, NC3.2447.0—7.597.78 PyroxeniteSonoma Co., CA3.2476.8—7.798.01 DuniteMt. Dun, New Zealand3.2587.57.697.808.00 DuniteBalsam Gap, NC3.2677.07.828.018.28 BronzititeStillwater Complex, MT3.2797.42—7.657.83 DuniteAddie, NC3.3047.70—8.058.28 DuniteTwin Sisters Peaks, WA3.3127.78.118.278.42 EclogiteTanganyika3.3286.647.307.467.71 JadeiteBurma3.3318.45—8.698.78 HarzburgiteBushveld Complex3.3696.97.747.817.95 EclogiteKimberley3.3767.177.657.737.87 EclogiteSunnmore, Norway3.3765.2—7.307.69 EclogiteHealdsburg, CA3.4417.31—7.818.01 GarnetCT3.5616.3—8.558.99 DuniteMoonihoek Mine, Transvaal3.7446.77.137.217.36