<p>I dug around quite a bit and was unable to find anyone making measurements of the type that would lead one to determine the moment of inertia of Uranus or Neptune without first making assumptions about their composition.  The way to determine the moment of inertia without such assumptions is to use the gravitational quadrupole moment, <span class="math-container">$J_2$</span>, and oblateness (see below).</p>
<p><a href="https://ui.adsabs.harvard.edu/abs/1928PASP...40..234M/abstract" rel="noreferrer">Moore &amp; Menzel (1928)</a> report on the rotation of Neptune.  They state that &quot;it is possible to calculate the oblateness and period of rotation, first making certain assumptions as to the planet's internal constitution. Thus Jackson,
assuming a constitution for Neptune similar to that of Jupiter, derives an [oblateness] of <span class="math-container">$1/65.7$</span>...&quot;</p>
<p>The same duo also reported on Uranus's rotation <a href="https://ui.adsabs.harvard.edu/abs/1930PASP...42..330M/abstract" rel="noreferrer">(Moore &amp; Menzel 1930)</a>, but they make no mention of the oblateness.</p>
<p><a href="https://ntrs.nasa.gov/citations/19810011459" rel="noreferrer">Goody (1981)</a> gives the most useful historical account of Uranian observations I could find.  They state:</p>
<blockquote>
<p>The first report of a rotation rate close to <span class="math-container">$11$</span> h, based on a theoretical analysis of the planet's figure, was by Berstrand (1909). <strong>At that time, there were no measurements of <span class="math-container">$J_2$</span> and no reliable data on the oblateness</strong>. The crucial early work was the spectrographic determination of <span class="math-container">$10.8 \pm 0.3$</span> h by Lowell and Slipher (1912).</p>
</blockquote>
<p>Based on this I would infer that any knowledge of the moments of inertia of Uranus and Neptune relied on hydrodynamical models of the planets' interiors, particularly the <a href="https://en.wikipedia.org/wiki/Darwin%E2%80%93Radau_equation" rel="noreferrer">Darwin-Radau Approximation</a>, which dates back to the 19th century.</p>
<p>Goody references a text, Alexander (1965), as having even more historical details.  If I can track that down, I'll update this.</p>
<hr />
<p>If you do want to determine a planet's moment of inertia remotely without assumptions, you'll need to precisely measure its gravitational field.</p>
<p>Any gravitational field can be expressed as a sum of <a href="https://en.wikipedia.org/wiki/Spherical_harmonics" rel="noreferrer">spherical harmonics</a>.  This was first worked out by <a href="https://en.wikipedia.org/wiki/Pierre-Simon_Laplace#On_the_figure_of_the_Earth" rel="noreferrer">Laplace in the late 18th century</a>.  The terms in this sum are often referred to as the <a href="https://en.wikipedia.org/wiki/Multipole_expansion" rel="noreferrer">multipole moments</a> of mass.  The first <span class="math-container">$\ell=0$</span> term is the mass monopole, which is proportional to the object's total mass.  The <span class="math-container">$\ell=1$</span> dipole terms vanish for the right choice of coordinates.  The <span class="math-container">$\ell=2$</span> quadrupole terms are related to the moments of inertia of the object.  The expansion continues to higher and higher terms.</p>
<p>As you suggest, to map out the multipoles of an object you could use geodesy.  Just watch something orbit the object in question and carefully map out the motion of the orbiter.  If the orbiter follows a perfect elliptical orbit, then it is moving in a purely monopolar field of spherical object.  If the orbiter deviates from a perfect ellipse, then there must be higher order corrections to describe the field.  The first non-zero correction to a monopole is the quadrupole, which encodes the moment of inertia.</p>
<p>For an axially symmetric object of mass <span class="math-container">$M$</span> and equatorial radius <span class="math-container">$R$</span> the quadrupole term is:
<span class="math-container">$$J_2 = \frac{I_z - I_x}{M\,R^2},$$</span>
where <span class="math-container">$I_z$</span> is the moment of inertia about the symmetric axis (i.e. rotational), and <span class="math-container">$I_x=I_y$</span> is the moment of inertia about an axis perpendicular to the symmetric axis (i.e. polar).  This vanishes for a spherically symmetric object, because <span class="math-container">$I_z=I_x$</span>.  To fully disentangle <span class="math-container">$I_z$</span> from <span class="math-container">$J_2$</span> you'd also need the oblateness.</p>
<p>Luckily Uranus and Neptune have moons, so you can do geodesy without sending a satellite flyby.  The earliest reports of <span class="math-container">$J_2$</span> and oblateness I could find are from the 1970s (e.g., <a href="https://ui.adsabs.harvard.edu/abs/1972Obs....92...84C/abstract" rel="noreferrer">Cook 1972</a>, <a href="https://ui.adsabs.harvard.edu/abs/1978ApJ...219.1079T/abstract" rel="noreferrer">Trauger et al 1978</a>).  This predates the <a href="https://en.wikipedia.org/wiki/Voyager_2#Encounter_with_Uranus" rel="noreferrer">Voyager 2 flyby in 1986</a>.</p>
<p>These days we can do better than <span class="math-container">$\ell=2$</span>.  <a href="https://ui.adsabs.harvard.edu/abs/2014AJ....148...76J/abstract" rel="noreferrer">Jacobson (2014)</a> used a wide range of historical and more recent data, including observations as far back as 1912 and the Voyager 2 flyby data, to determine many parameters of the Uranian system, including multipole moments up to <span class="math-container">$\ell=6$</span>.</p>
<p>With dedicated orbiters we can do even better.  For example the <a href="https://en.wikipedia.org/wiki/Earth_Gravitational_Model#EGM2008" rel="noreferrer">EGM2008 model</a> for the Earth's gravitational field (based on data from the <a href="https://en.wikipedia.org/wiki/GRACE_and_GRACE-FO" rel="noreferrer">GRACE</a> mission) goes up to <span class="math-container">$\ell=2159$</span>.</p>

CC BY-SA 4.0
2023-10-31T08:03:12.787Z
Doesn't that fail for a spherical object (which also have moment of inertia). My guess is that you must measure changes od rotation.
Leos Ondra
CC BY-SA 4.0
2023-11-02T19:11:09.663Z
I've edited the answer to incorporate updates based on the comments.
Paul T.
CC BY-SA 4.0
2023-11-04T10:44:41.607Z
I've upvoted you answer as I provides a lot of useful information and corresponding context.  Unfortunately I'm still in search of what could have been accomplished in the 1930s.  I may reach out to S&T and ask if they have a reference as in the article the context is to use the MOI to estimate the composition.
GrapefruitIsAwesome
CC BY-SA 4.0
2023-10-30T18:47:32.640Z
2023-11-02T19:12:33.423Z
2023-11-02T19:12:33.423Z
Paul T.
<p>In the <a href="https://skyandtelescope.org/sky-and-telescope-magazine/inside-the-july-2023-issue/" rel="nofollow noreferrer">July 2023 issue</a> of <a href="https://skyandtelescope.org" rel="nofollow noreferrer">Sky &amp; Telescope</a> there is an article titled <a href="https://www.nxtbook.com/nxtbooks/aas/st_202307/index.php#/p/14" rel="nofollow noreferrer">Sights Set on Uranus</a>.  In that article the following statement is made:</p>
<blockquote>
<p>Since the 1930's we've suspected that Uranus and Neptune are made
mostly of ice. (Ice refers to materials that are typically liquids or
gases on Earth but are frozen in the outer solar system, including
water and other lightweight molecular compounds like methane and
ammonia.) But that time, astronomers had measured each planet's mass,
volume, and <em>moment of inertia</em>, a measure of how concentrated the
mass is towards the planet's center.</p>
</blockquote>
<p>I can easily imagine how one would estimate the mass and volume of a planet remotely in the 1930s, but how was the moment of inertia estimated?  I would have thought this would have required an orbiter or a fly-by to measure the shape of the gravitational potential.</p>
<hr>
<p>Information below is more background information that I was able to dig up:</p>
<p>Wildt, R., 1939. <a href="https://ui.adsabs.harvard.edu/abs/1947MNRAS.107...84W/abstract" rel="nofollow noreferrer">The Constitution of the Planets</a>. <em>Proceedings of the American Philosophical Society, 81</em>(2), pp.135-152.</p>
<blockquote>
<p>The greatest moment of inertia, <span class="math-container">$C$</span>, of a planet can be computed
from its oblateness on the Radau-Darwin theory of the
ellipticity of the&gt; Earth.<sup>*</sup>  For a homogeneous body
the dimensionless quantity, <span class="math-container">$C/Ma^2$</span>, has the maximum value,
2/5. The values of <span class="math-container">$C/Ma^2$</span> given in Table I for the
terrestrial planets have been taken from H. Jeffreys <span class="math-container">$\dagger$</span>;
they correspond to his models of these bodies which will be
discussed in Section II. The other values have been calculated
from the dynamical ellipticities of the giant planets
determined by N. Lvoff.<span class="math-container">$\ddagger$</span> Since the deviations of the
individual values from Lvoff’s mean, 0.24, are not well
established at present, only this mean will be used in the
further discussions; but there is evidence that the quantity
<span class="math-container">$C/Ma^2$</span> is decidedly less for Saturn than for Jupiter.</p>
</blockquote>
<p>Note: <span class="math-container">$C$</span> - greatest moment of inertia; <span class="math-container">$M$</span> - mass, unit 5.966x10<sup>27</sup> g; <span class="math-container">$a$</span> - equatorial radius, unit 6278.388 km</p>
<p>* G. H. Darwin, <em>M.N.</em>, <strong>60</strong>, 82, 1900; H. Jeffreys, <em>M.N.</em>, <strong>84</strong>, 534, 1924.<br>
<span class="math-container">$\dagger$</span> H. Jeffreys, <em>M.N., Geophys. Suppl.</em>, <strong>4</strong>, 62, 1940.<br>
<span class="math-container">$\ddagger$</span> N. Lvoff, <em>Russ. Astr. J.</em>, <strong>9</strong>, 68, 1932.<br></p>
<p>There three references I am currently struggling to track down.</p>

CC BY-SA 4.0
2023-11-11T13:03:03.253Z
Inspired by [Paul T.](https://astronomy.stackexchange.com/users/43713/paul-t)'s answer I started down the rabbit hole of what was done in in and prior to the 1930s.  So far I've arrived at [The Constitution of the planets](https://ui.adsabs.harvard.edu/abs/1947MNRAS.107...84W/abstract) from 1947, but its hard to chase down the older references from that paper.  Still working...
GrapefruitIsAwesome
<planet><solar-system><planetary-science>
How is a Planet's Moment of Inertia Measured Remotely?
GrapefruitIsAwesome
# How is a Planet's Moment of Inertia Measured Remotely?

(Asked by: username_0 on 2023-10-28)

In the [July 2023 issue](https://skyandtelescope.org/sky-and-telescope-magazine/inside-the-july-2023-issue/) of [Sky & Telescope](https://skyandtelescope.org) there is an article titled [Sights Set on Uranus](https://www.nxtbook.com/nxtbooks/aas/st_202307/index.php#/p/14). In that article the following statement is made:

> Since the 1930's we've suspected that Uranus and Neptune are made
> mostly of ice. (Ice refers to materials that are typically liquids or
> gases on Earth but are frozen in the outer solar system, including
> water and other lightweight molecular compounds like methane and
> ammonia.) But that time, astronomers had measured each planet's mass,
> volume, and *moment of inertia*, a measure of how concentrated the
> mass is towards the planet's center.

I can easily imagine how one would estimate the mass and volume of a planet remotely in the 1930s, but how was the moment of inertia estimated? I would have thought this would have required an orbiter or a fly-by to measure the shape of the gravitational potential.

---

Information below is more background information that I was able to dig up:

Wildt, R., 1939. [The Constitution of the Planets](https://ui.adsabs.harvard.edu/abs/1947MNRAS.107...84W/abstract). *Proceedings of the American Philosophical Society, 81*(2), pp.135-152.

> The greatest moment of inertia, $C$, of a planet can be computed
> from its oblateness on the Radau-Darwin theory of the
> ellipticity of the> Earth.\* For a homogeneous body
> the dimensionless quantity, $C/Ma^2$, has the maximum value,
> 2/5. The values of $C/Ma^2$ given in Table I for the
> terrestrial planets have been taken from H. Jeffreys $\dagger$;
> they correspond to his models of these bodies which will be
> discussed in Section II. The other values have been calculated
> from the dynamical ellipticities of the giant planets
> determined by N. Lvoff.$\ddagger$ Since the deviations of the
> individual values from Lvoff’s mean, 0.24, are not well
> established at present, only this mean will be used in the
> further discussions; but there is evidence that the quantity
> $C/Ma^2$ is decidedly less for Saturn than for Jupiter.

Note: $C$ - greatest moment of inertia; $M$ - mass, unit 5.966x1027 g; $a$ - equatorial radius, unit 6278.388 km

\* G. H. Darwin, *M.N.*, **60**, 82, 1900; H. Jeffreys, *M.N.*, **84**, 534, 1924.  
$\dagger$ H. Jeffreys, *M.N., Geophys. Suppl.*, **4**, 62, 1940.  
$\ddagger$ N. Lvoff, *Russ. Astr. J.*, **9**, 68, 1932.

There three references I am currently struggling to track down.

**Tags:** planet, solar-system, planetary-science
> **Comments:**
> **username_0**: Inspired by [Paul T.](https://astronomy.stackexchange.com/users/43713/paul-t)'s answer I started down the rabbit hole of what was done in in and prior to the 1930s.  So far I've arrived at [The Constitution of the planets](https://ui.adsabs.harvard.edu/abs/1947MNRAS.107...84W/abstract) from 1947, but its hard to chase down the older references from that paper.  Still working...  

---

## 1 answers

### Answer by username_1 on 2023-10-30. Score: 9

I dug around quite a bit and was unable to find anyone making measurements of the type that would lead one to determine the moment of inertia of Uranus or Neptune without first making assumptions about their composition. The way to determine the moment of inertia without such assumptions is to use the gravitational quadrupole moment, $J\_2$, and oblateness (see below).

[Moore & Menzel (1928)](https://ui.adsabs.harvard.edu/abs/1928PASP...40..234M/abstract) report on the rotation of Neptune. They state that "it is possible to calculate the oblateness and period of rotation, first making certain assumptions as to the planet's internal constitution. Thus Jackson,
assuming a constitution for Neptune similar to that of Jupiter, derives an [oblateness] of $1/65.7$..."

The same duo also reported on Uranus's rotation [(Moore & Menzel 1930)](https://ui.adsabs.harvard.edu/abs/1930PASP...42..330M/abstract), but they make no mention of the oblateness.

[Goody (1981)](https://ntrs.nasa.gov/citations/19810011459) gives the most useful historical account of Uranian observations I could find. They state:

> The first report of a rotation rate close to $11$ h, based on a theoretical analysis of the planet's figure, was by Berstrand (1909). **At that time, there were no measurements of $J\_2$ and no reliable data on the oblateness**. The crucial early work was the spectrographic determination of $10.8 \pm 0.3$ h by Lowell and Slipher (1912).

Based on this I would infer that any knowledge of the moments of inertia of Uranus and Neptune relied on hydrodynamical models of the planets' interiors, particularly the [Darwin-Radau Approximation](https://en.wikipedia.org/wiki/Darwin%E2%80%93Radau_equation), which dates back to the 19th century.

Goody references a text, Alexander (1965), as having even more historical details. If I can track that down, I'll update this.

---

If you do want to determine a planet's moment of inertia remotely without assumptions, you'll need to precisely measure its gravitational field.

Any gravitational field can be expressed as a sum of [spherical harmonics](https://en.wikipedia.org/wiki/Spherical_harmonics). This was first worked out by [Laplace in the late 18th century](https://en.wikipedia.org/wiki/Pierre-Simon_Laplace#On_the_figure_of_the_Earth). The terms in this sum are often referred to as the [multipole moments](https://en.wikipedia.org/wiki/Multipole_expansion) of mass. The first $\ell=0$ term is the mass monopole, which is proportional to the object's total mass. The $\ell=1$ dipole terms vanish for the right choice of coordinates. The $\ell=2$ quadrupole terms are related to the moments of inertia of the object. The expansion continues to higher and higher terms.

As you suggest, to map out the multipoles of an object you could use geodesy. Just watch something orbit the object in question and carefully map out the motion of the orbiter. If the orbiter follows a perfect elliptical orbit, then it is moving in a purely monopolar field of spherical object. If the orbiter deviates from a perfect ellipse, then there must be higher order corrections to describe the field. The first non-zero correction to a monopole is the quadrupole, which encodes the moment of inertia.

For an axially symmetric object of mass $M$ and equatorial radius $R$ the quadrupole term is:
$$J\_2 = \frac{I\_z - I\_x}{M\,R^2},$$
where $I\_z$ is the moment of inertia about the symmetric axis (i.e. rotational), and $I\_x=I\_y$ is the moment of inertia about an axis perpendicular to the symmetric axis (i.e. polar). This vanishes for a spherically symmetric object, because $I\_z=I\_x$. To fully disentangle $I\_z$ from $J\_2$ you'd also need the oblateness.

Luckily Uranus and Neptune have moons, so you can do geodesy without sending a satellite flyby. The earliest reports of $J\_2$ and oblateness I could find are from the 1970s (e.g., [Cook 1972](https://ui.adsabs.harvard.edu/abs/1972Obs....92...84C/abstract), [Trauger et al 1978](https://ui.adsabs.harvard.edu/abs/1978ApJ...219.1079T/abstract)). This predates the [Voyager 2 flyby in 1986](https://en.wikipedia.org/wiki/Voyager_2#Encounter_with_Uranus).

These days we can do better than $\ell=2$. [Jacobson (2014)](https://ui.adsabs.harvard.edu/abs/2014AJ....148...76J/abstract) used a wide range of historical and more recent data, including observations as far back as 1912 and the Voyager 2 flyby data, to determine many parameters of the Uranian system, including multipole moments up to $\ell=6$.

With dedicated orbiters we can do even better. For example the [EGM2008 model](https://en.wikipedia.org/wiki/Earth_Gravitational_Model#EGM2008) for the Earth's gravitational field (based on data from the [GRACE](https://en.wikipedia.org/wiki/GRACE_and_GRACE-FO) mission) goes up to $\ell=2159$.

> **Comments:**
> **username_2**: Doesn't that fail for a spherical object (which also have moment of inertia). My guess is that you must measure changes od rotation.  
> **username_1**: I've edited the answer to incorporate updates based on the comments.  
> **username_0**: I've upvoted you answer as I provides a lot of useful information and corresponding context.  Unfortunately I'm still in search of what could have been accomplished in the 1930s.  I may reach out to S&T and ask if they have a reference as in the article the context is to use the MOI to estimate the composition.  

---

