Orbits, Gravity and Planetary Interiors
**Guiding question:** How do gravitational measurements reveal motions, shapes and hidden structures? Orbital geometry constrains dynamical relationships while gravity, shape, rotation and tidal response constrain internal mass distribution. An analytic two-body orbit is a benchmark, not a complete ephemeris. Interior
§2.1 Learning outcomes
By the end of this unit, you should be able to describe the dominant physical interactions without collapsing all mechanisms into a single cause; translate an observable into a measurable quantity with units; design a synthetic comparison that varies one quantity at a time; identify at least two independent failure modes of an inference; and communicate why a simple classroom model is not an operational planetary prediction.
§2.2 Five lessons
2.1 Keplerian orbits and the two-body benchmark
Central concept. For two isolated point masses the orbital period, semimajor axis and total gravitational parameter obey Kepler’s third law. Real solar-system ephemerides include perturbations, nonspherical gravity and changing reference frames. In this part of orbits, gravity and planetary interiors, the important skill is to connect a plausible physical mechanism to a particular, imperfect record. The phrase explains the evidence should never mean merely that a sketch resembles a photograph. A mechanism must be compatible with constraints on size, time, composition and energetics. It should predict something that another observation could disprove. The expected strength of the inference depends on the quality of the instrument and the uniqueness of the prediction, not on the visual appeal of a diagram.
What would actually be measured? Compute circular-orbit periods around a specified gravitating mass and contrast with an eccentric test orbit. This is a constructed teaching example; no row should be mistaken for a spacecraft measurement. Define the measurable quantities before choosing a fitting equation, and preserve calibration units, coordinate frame and acquisition conditions. For a remote world, resolution may change over the scene and the hidden part of a terrain or atmosphere may have fundamentally different characteristics. Report not only a central estimate but also which conditions remained unobserved and which parameters were held fixed for convenience.
Failure mode and competing account. An AU-year scaling assumes solar mass and a heliocentric, negligible-test-mass approximation. The strongest competing explanation is not an arbitrary alternative chosen to look weak. It should be physically plausible under the same observing circumstances. Ask which observable, spatial relationship or laboratory property differs between the explanations. If both models predict nearly the same signal at the available resolution, the honest result is non-identifiability. A mission proposal can still be useful in this situation by explicitly targeting that degeneracy. Do not increase certainty simply because repeated measurements use the same processing assumptions.
Worked quantitative case — Circular orbital period. Consider the fictional parameter value 0.6 AU. Under the stated classroom approximation, Keplerian two-body relation around a solar-mass central body. Substituting the input gives 0.464758 years. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Check your understanding. (1) Which assumption most strongly controls the numerical estimate for this lesson? (2) Describe how to distinguish the chosen physical process from its most plausible confounder using a new measurement. (3) Which claim would remain unjustified even if the calculation were arithmetically perfect? (4) Rewrite the result as a testable statement with units and explicit uncertainty. (5) State the advantage and the limitation of transferring an Earth-based analogy to this setting.
2.2 Surface gravity and escape speed
Central concept. Surface acceleration for a spherical mass distribution is GM/R² and ideal escape speed is the square root of 2GM/R. Escape speed characterizes energy, not the long-term persistence of an atmosphere under irradiation. In this part of orbits, gravity and planetary interiors, the important skill is to connect a plausible physical mechanism to a particular, imperfect record. The phrase explains the evidence should never mean merely that a sketch resembles a photograph. A mechanism must be compatible with constraints on size, time, composition and energetics. It should predict something that another observation could disprove. The expected strength of the inference depends on the quality of the instrument and the uniqueness of the prediction, not on the visual appeal of a diagram.
What would actually be measured? Vary a hypothetical rocky body’s mass while keeping radius fixed and compute both gravity and ideal escape speed. This is a constructed teaching example; no row should be mistaken for a spacecraft measurement. Define the measurable quantities before choosing a fitting equation, and preserve calibration units, coordinate frame and acquisition conditions. For a remote world, resolution may change over the scene and the hidden part of a terrain or atmosphere may have fundamentally different characteristics. Report not only a central estimate but also which conditions remained unobserved and which parameters were held fixed for convenience.
Failure mode and competing account. A rotating irregular body may have variable effective gravity, and high exobase temperatures permit loss processes beyond a single cutoff. The strongest competing explanation is not an arbitrary alternative chosen to look weak. It should be physically plausible under the same observing circumstances. Ask which observable, spatial relationship or laboratory property differs between the explanations. If both models predict nearly the same signal at the available resolution, the honest result is non-identifiability. A mission proposal can still be useful in this situation by explicitly targeting that degeneracy. Do not increase certainty simply because repeated measurements use the same processing assumptions.
Worked quantitative case — Ideal escape speed. Consider the fictional parameter value 0.8 Earth masses. Under the stated classroom approximation, Spherical nonrotating model at fixed Earth radius. Substituting the input gives 10.0051 km/s. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Check your understanding. (1) Which assumption most strongly controls the numerical estimate for this lesson? (2) Describe how to distinguish the chosen physical process from its most plausible confounder using a new measurement. (3) Which claim would remain unjustified even if the calculation were arithmetically perfect? (4) Rewrite the result as a testable statement with units and explicit uncertainty. (5) State the advantage and the limitation of transferring an Earth-based analogy to this setting.
2.3 Spin, obliquity and tidal torques
Central concept. Rotation produces equatorial acceleration and can redistribute oceans, atmospheres and interior mass. Obliquity modifies seasonal forcing; tidal torques exchange angular momentum between spin and orbit. In this part of orbits, gravity and planetary interiors, the important skill is to connect a plausible physical mechanism to a particular, imperfect record. The phrase explains the evidence should never mean merely that a sketch resembles a photograph. A mechanism must be compatible with constraints on size, time, composition and energetics. It should predict something that another observation could disprove. The expected strength of the inference depends on the quality of the instrument and the uniqueness of the prediction, not on the visual appeal of a diagram.
What would actually be measured? Construct illumination geometry for a tilted synthetic planet and estimate centrifugal acceleration at a chosen latitude. This is a constructed teaching example; no row should be mistaken for a spacecraft measurement. Define the measurable quantities before choosing a fitting equation, and preserve calibration units, coordinate frame and acquisition conditions. For a remote world, resolution may change over the scene and the hidden part of a terrain or atmosphere may have fundamentally different characteristics. Report not only a central estimate but also which conditions remained unobserved and which parameters were held fixed for convenience.
Failure mode and competing account. Spin evolution depends on internal dissipation and orbital companions; synchronous rotation is not guaranteed at all distances. The strongest competing explanation is not an arbitrary alternative chosen to look weak. It should be physically plausible under the same observing circumstances. Ask which observable, spatial relationship or laboratory property differs between the explanations. If both models predict nearly the same signal at the available resolution, the honest result is non-identifiability. A mission proposal can still be useful in this situation by explicitly targeting that degeneracy. Do not increase certainty simply because repeated measurements use the same processing assumptions.
Worked quantitative case — Equatorial centrifugal acceleration. Consider the fictional parameter value 12 hours per rotation. Under the stated classroom approximation, a = ω²R with R=6371 km; not latitude-dependent effective gravity. Substituting the input gives 0.134772 m/s². This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Check your understanding. (1) Which assumption most strongly controls the numerical estimate for this lesson? (2) Describe how to distinguish the chosen physical process from its most plausible confounder using a new measurement. (3) Which claim would remain unjustified even if the calculation were arithmetically perfect? (4) Rewrite the result as a testable statement with units and explicit uncertainty. (5) State the advantage and the limitation of transferring an Earth-based analogy to this setting.
2.4 Differentiation, density and moments of inertia
Central concept. Rock and metal segregation, melting and pressure-dependent phase changes make layered interiors possible. Bulk density constrains total mass per volume but cannot locate dense layers independently. In this part of orbits, gravity and planetary interiors, the important skill is to connect a plausible physical mechanism to a particular, imperfect record. The phrase explains the evidence should never mean merely that a sketch resembles a photograph. A mechanism must be compatible with constraints on size, time, composition and energetics. It should predict something that another observation could disprove. The expected strength of the inference depends on the quality of the instrument and the uniqueness of the prediction, not on the visual appeal of a diagram.
What would actually be measured? Create two mass-balanced two-layer spheres with the same mean density but different core radius and compare predicted moment-of-inertia factors. This is a constructed teaching example; no row should be mistaken for a spacecraft measurement. Define the measurable quantities before choosing a fitting equation, and preserve calibration units, coordinate frame and acquisition conditions. For a remote world, resolution may change over the scene and the hidden part of a terrain or atmosphere may have fundamentally different characteristics. Report not only a central estimate but also which conditions remained unobserved and which parameters were held fixed for convenience.
Failure mode and competing account. Treating density as constant within each layer masks compressibility, temperature and composition gradients. The strongest competing explanation is not an arbitrary alternative chosen to look weak. It should be physically plausible under the same observing circumstances. Ask which observable, spatial relationship or laboratory property differs between the explanations. If both models predict nearly the same signal at the available resolution, the honest result is non-identifiability. A mission proposal can still be useful in this situation by explicitly targeting that degeneracy. Do not increase certainty simply because repeated measurements use the same processing assumptions.
Worked quantitative case — Bulk density of a modeled sphere. Consider the fictional parameter value 0.66 10²⁴ kg. Under the stated classroom approximation, Uniform bulk density from a 3000 km radius, not an internal density profile. Substituting the input gives 5,835.68 kg/m³. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Check your understanding. (1) Which assumption most strongly controls the numerical estimate for this lesson? (2) Describe how to distinguish the chosen physical process from its most plausible confounder using a new measurement. (3) Which claim would remain unjustified even if the calculation were arithmetically perfect? (4) Rewrite the result as a testable statement with units and explicit uncertainty. (5) State the advantage and the limitation of transferring an Earth-based analogy to this setting.
2.5 Interior heat and long-term cooling
Central concept. Primordial accretional heat, radioactive decay, secular cooling and tides can all contribute to planetary energy budgets. Transport occurs by conduction, convection and advection where material properties allow it. In this part of orbits, gravity and planetary interiors, the important skill is to connect a plausible physical mechanism to a particular, imperfect record. The phrase explains the evidence should never mean merely that a sketch resembles a photograph. A mechanism must be compatible with constraints on size, time, composition and energetics. It should predict something that another observation could disprove. The expected strength of the inference depends on the quality of the instrument and the uniqueness of the prediction, not on the visual appeal of a diagram.
What would actually be measured? Evaluate a synthetic conductive heat flow through a crust for two thermal conductivity assumptions. This is a constructed teaching example; no row should be mistaken for a spacecraft measurement. Define the measurable quantities before choosing a fitting equation, and preserve calibration units, coordinate frame and acquisition conditions. For a remote world, resolution may change over the scene and the hidden part of a terrain or atmosphere may have fundamentally different characteristics. Report not only a central estimate but also which conditions remained unobserved and which parameters were held fixed for convenience.
Failure mode and competing account. A present surface heat flow does not uniquely reveal heat production because insulation and mantle transport change with time. The strongest competing explanation is not an arbitrary alternative chosen to look weak. It should be physically plausible under the same observing circumstances. Ask which observable, spatial relationship or laboratory property differs between the explanations. If both models predict nearly the same signal at the available resolution, the honest result is non-identifiability. A mission proposal can still be useful in this situation by explicitly targeting that degeneracy. Do not increase certainty simply because repeated measurements use the same processing assumptions.
Worked quantitative case — Crustal conductive flux. Consider the fictional parameter value 26 km conductive layer. Under the stated classroom approximation, Fourier q = k ΔT / L, k=2.5 W/(m K), ΔT=700 K. Substituting the input gives 0.0673077 W/m². This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Check your understanding. (1) Which assumption most strongly controls the numerical estimate for this lesson? (2) Describe how to distinguish the chosen physical process from its most plausible confounder using a new measurement. (3) Which claim would remain unjustified even if the calculation were arithmetically perfect? (4) Rewrite the result as a testable statement with units and explicit uncertainty. (5) State the advantage and the limitation of transferring an Earth-based analogy to this setting.
§2.3 Integrated investigation: Orbits, Gravity and Planetary Interiors
Research question. How do gravitational measurements reveal motions, shapes and hidden structures? Use the fabricated evidence packet generated for this unit, not a downloaded mission observation, to test the issue. The packet provides five numerical inputs corresponding to five distinct lesson calculations. Prepare an analysis log before computing results: specify which data are invented, which constants are adopted from reference physics, which quantities are exact in the model, and which are uncertain in an actual measurement. Do not call synthetic demonstration values planetary facts.
Workflow A — observations and provenance. Build a small metadata table with columns for source, timestamp or hypothetical epoch, spatial footprint, units, calibration state, quality flag and an explicit 'synthetic' indicator. Check that every row belongs to the same stated physical or geometrical reference when the comparison requires that. If different observations cover different regions or epochs, say what assumption permits their comparison. A useful practice is to write down the interpretation you expect before opening a calculated output; doing so can reveal a tendency to favor a result after seeing it.
Workflow B — model construction. Start with one conservation, geometry or rate relationship from the unit. Derive the form in symbols, then set parameter values. Write an extreme-case check: an input tending to zero, a very long orbital or cooling time, or a hypothetical measurement at the edge of the instrument's range. Not every extreme has physical meaning, but the exercise can reveal sign errors or a model's restricted domain. Compare at least one dimensionless ratio where possible. Record any calibration or geological assumptions that cannot be derived from the formula itself.
Workflow C — alternatives and sensitivity. Construct a second explanation able to reproduce at least one headline feature of the mock dataset. Change a single uncertain parameter by plus or minus ten percent and tabulate the output. Then alter the model structure itself: for example, replace uniform material properties with a layered assumption, permit a nonstationary flux, or introduce a time-varying forcing condition. When the inferred history changes more from model structure than from stated measurement error, the report must say that the quantity is weakly identified. A precise fitted value is not automatically a precise scientific claim.
Workflow D — cross-unit integration. Connect this unit to two other themes in the book. For example, relate a dynamical interpretation to geological chronology, a spectral inference to atmospheric chemistry, or an internal energy source to an icy-surface observation. Each connection should contain one shared variable and one new source of uncertainty. This prevents a detailed but isolated calculation from being treated as a full planetary system model. If a comparison relies on terrestrial experiments, identify the gravitational, thermal or compositional regime where the analogy ceases to be straightforward.
Workflow E — reporting and evaluation. Submit a one-page methods statement, a table of computed values, one plot with labeled axes, an uncertainty inventory, and a claim supported by the data packet. Include an explicit alternative interpretation and a feasible next observation. The teacher should assess whether units, underlying equations, provenance and falsifiability are correct; artistic polish alone is not evidence of scientific rigor. Use language such as consistent with, constrained by and not distinguished by these data when warranted. Do not advertise an invented training result as mission validation.
§2.4 Unit synthesis and research discussion
The scientific difficulty of orbits, gravity and planetary interiors is not simply locating one more fact. It is learning what the available observations exclude. A strong synthesis distinguishes a process being physically possible from that process being historically demonstrated. It also distinguishes a local interpretation from a global planetary pattern and a present state from the history that formed it. Return to the unit's guiding question and state three claims at different strengths: one direct observation, one explanation that can be tested, and one conjecture for which available information is insufficient. For each, name the likely effect of improved instrumental precision, better geographic coverage and stronger physical modeling. These improvements address different errors and should not be treated as interchangeable. Explore whether the most valuable next mission measurement is an imaging survey, an in-situ measurement, a laboratory calibration or repeated monitoring.
Capstone prompt. You have resources for one additional investigation only. Present the competing hypotheses in a table and estimate which candidate observation would most clearly separate their predictions per unit cost and risk. Your justification must acknowledge a scenario in which your chosen observation gives an ambiguous result. A responsible conclusion is an explicit account of what remains unknown, not an invented certainty.
Step-by-Step Statistical Simulation and Monte Carlo Solutions
Three original questions connect calculations, definitions and interpretation. Open each solution to follow the reasoning.
Consider the fictional parameter value 0.6 AU. Under the stated classroom approximation, Keplerian two-body relation around a solar-mass central body. Substituting the input gives 0.464758 years. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 0.6 AU. Under the stated classroom approximation, Keplerian two-body relation around a solar-mass central body. Substituting the input gives 0.464758 years. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 0.8 Earth masses. Under the stated classroom approximation, Spherical nonrotating model at fixed Earth radius. Substituting the input gives 10.0051 km/s. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 0.8 Earth masses. Under the stated classroom approximation, Spherical nonrotating model at fixed Earth radius. Substituting the input gives 10.0051 km/s. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 12 hours per rotation. Under the stated classroom approximation, a = ω²R with R=6371 km; not latitude-dependent effective gravity. Substituting the input gives 0.134772 m/s². This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 12 hours per rotation. Under the stated classroom approximation, a = ω²R with R=6371 km; not latitude-dependent effective gravity. Substituting the input gives 0.134772 m/s². This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 0.66 10²⁴ kg. Under the stated classroom approximation, Uniform bulk density from a 3000 km radius, not an internal density profile. Substituting the input gives 5,835.68 kg/m³. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 0.66 10²⁴ kg. Under the stated classroom approximation, Uniform bulk density from a 3000 km radius, not an internal density profile. Substituting the input gives 5,835.68 kg/m³. This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 26 km conductive layer. Under the stated classroom approximation, Fourier q = k ΔT / L, k=2.5 W/(m K), ΔT=700 K. Substituting the input gives 0.0673077 W/m². This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.
Consider the fictional parameter value 26 km conductive layer. Under the stated classroom approximation, Fourier q = k ΔT / L, k=2.5 W/(m K), ΔT=700 K. Substituting the input gives 0.0673077 W/m². This number is neither a retrieved physical parameter nor a validated forecast. Its role is to expose scale and sensitivity. First confirm that the units on both sides agree; next estimate how the result changes if the input moves by ten percent while all other conditions are fixed. Finally ask whether the model has an appropriate domain for extrapolation. A numerical answer without those checks can be exact arithmetic attached to an incorrect physical story.