Skip to content

Commit adea0a7

Browse files
L08 review fixes: Saturn layer scale + Garnet v_P
Adversarial physics review findings folded in: - fig:saturn-schematic: layer boundaries shifted from diffuse core 0.30, metallic-H 0.55, He-rain 0.62 (which placed the molecular/metallic H transition deeper than published models) to diffuse core 0.55, metallic-H 0.70, He-rain 0.78, matching Mankovich 2021 ring seismology + Guillot 1999 + Nettelmann 2013 + Helled 2020. Annotation anchors updated to point at the correct mid-radii. - birch_minerals.json: Garnet v_P 8.8 -> 9.1 km/s. Pyrope (the mantle-relevant species) measures at 9.0-9.3 km/s in Bass 1995 / Karki 2001; the previous 8.8 placed Garnet 0.5 km/s below the Birch trend, making it look like a near-outlier.
1 parent 89990b5 commit adea0a7

4 files changed

Lines changed: 9 additions & 6 deletions

File tree

-37 Bytes
Binary file not shown.
4.43 KB
Binary file not shown.

scripts/figures/L08_interiors/data/birch_minerals.json

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -8,7 +8,7 @@
88
{"name": "Forsterite", "density": 3.22, "v_P": 8.6, "outlier": false},
99
{"name": "Diopside", "density": 3.28, "v_P": 7.7, "outlier": false},
1010
{"name": "Periclase", "density": 3.58, "v_P": 9.7, "outlier": false},
11-
{"name": "Garnet", "density": 3.60, "v_P": 8.8, "outlier": false},
11+
{"name": "Garnet", "density": 3.60, "v_P": 9.1, "outlier": false},
1212
{"name": "Ringwoodite", "density": 3.72, "v_P": 9.8, "outlier": false},
1313
{"name": "Corundum", "density": 3.99, "v_P": 10.8, "outlier": false},
1414
{"name": "Bridgmanite", "density": 4.10, "v_P": 10.5, "outlier": false},

scripts/figures/L08_interiors/fig_saturn_schematic.py

Lines changed: 8 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -22,12 +22,15 @@
2222
REPO_ROOT = Path(__file__).resolve().parents[3]
2323
OUT_AVIF = REPO_ROOT / "book/08_interiors/figures/saturn_interior_schematic.avif"
2424

25-
# Concentric layers (outer to inner): label, outer fractional radius, color
25+
# Concentric layers (outer to inner): label, outer fractional radius, color.
26+
# Boundaries follow the Mankovich 2021 ring-seismology fuzzy core extending
27+
# to ~0.60 R_Sat, with the molecular/metallic-H transition at ~0.70 R_Sat
28+
# bracketing the He-rain layer (Guillot 1999; Nettelmann 2013; Helled 2020).
2629
LAYERS = [
2730
("Molecular H$_2$\ngas to liquid", 1.00, "#bca3d6"),
28-
("Helium-rain layer\nHe droplets settle", 0.62, "#f0e3a8"),
29-
("Metallic hydrogen\n(conducting; dynamo source)", 0.55, "#7a559e"),
30-
("Diffuse core (gradational)\n~17 M$_\\oplus$ heavy elements", 0.30, "#c08e57"),
31+
("Helium-rain layer\nHe droplets settle", 0.78, "#f0e3a8"),
32+
("Metallic hydrogen\n(conducting; dynamo source)", 0.70, "#7a559e"),
33+
("Diffuse core (gradational)\n~17 M$_\\oplus$ heavy elements", 0.55, "#c08e57"),
3134
]
3235

3336

@@ -42,7 +45,7 @@ def make_plot() -> Path:
4245
# Annotation arrows (right side)
4346
label_xs = [1.2, 1.2, 1.2, 1.2]
4447
label_ys = [0.78, 0.30, -0.20, -0.78]
45-
layer_anchor_rs = [0.85, 0.59, 0.45, 0.18]
48+
layer_anchor_rs = [0.89, 0.74, 0.625, 0.27]
4649
for (label, _, _), x_l, y_l, r_anchor in zip(LAYERS, label_xs, label_ys,
4750
layer_anchor_rs):
4851
ax.annotate(label, xy=(r_anchor * 0.5, r_anchor * 0.5),

0 commit comments

Comments
 (0)