{"id":"7964adda-513a-46f4-b315-abc122ab209a","arxiv_id":"2506.19152","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A single-parameter sublimation model reproduces most large-scale observed patterns of Io's SO2 atmosphere and predicts a seasonal hemispheric asymmetry tied to Io's tilt.","lead":"This paper models Io's thin sulfur dioxide atmosphere as produced entirely by sunlight-driven sublimation of surface frost, using a simplified heat-diffusion model with one adjustable parameter. Its results support a largely sublimation-controlled atmosphere and predict a seasonal north-south asymmetry that future observations can test.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (8), the core T-to-column mapping, is inconsistent as printed: at the model's stated 110-116 K dayside temperatures it yields 1e19-1e20 cm^-2, not the reported 3.7e16-8.5e16 cm^-2; the prefactor appears ~10^3 too large.","rationale":"The reader's conditional verdict is appropriate and I do not move it, but the sharpest weakly supported link is more specific than 'vapor-pressure equilibrium.' The equilibrium caveat in Appendix A3 is a known and honestly stated limitation; winds and patchy frost could smooth the T-N relation, yet the observed eclipse collapse and the Walker et al. frost-covered thermal inertia estimates give that assumption some independent support. By contrast, Eq. (8) is directly checkable and appears inconsistent with the paper's own outputs: at the dayside temperatures the text reports (about 110-116 K), the printed formula gives columns around 1e19-1e20 cm^-2, two to three orders above the 1e16-1e17 cm^-2 values used in all comparison figures. This makes it impossible to audit the main claim from the manuscript alone. The proposed test is decisive but cheap, because the code is public. If the code constant differs from the printed one, the correction is a text-level fix and the seasonal out-of-sample prediction and the eclipse asymmetry argument remain valuable; if not, the central comparison fails. That uncertainty is exactly why 'conditional' is the right verdict. The reader already flagged Eq. (8) as unverifiable, but did not make it the weakest assumption, hence partial agreement.","tokens_in":29448,"tokens_out":16494,"duration_ms":161522,"concrete_test":"Run the public Zenodo/GitHub code for the default case at the sub-Jovian point around perihelion; print the surface temperature and the saturation-column value at that location. Verify whether the code constant is 1.1516e13 or 1.1516e16 (or unit-equivalent) and compare the code's dayside maxima with Eq. (8) evaluated at the printed T. If the code uses a corrected prefactor and reproduces Figure 5, the physics argument survives but Eq. (8) must be corrected; if the code uses the printed prefactor, the modeled columns are 2-3 orders of magnitude too high and the central quantitative comparison cannot be supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2's Eq. (8) is the only link between the simulated surface temperature and the sublimation column density that the central claim rests on, so its numerical content is load-bearing. Plugging in the paper's own values—m_SO2 = 1.064e-25 kg, g0 = 1.81 m s^-2, T = 110 K—gives N_sub = 9.4e18 cm^-2; at T = 116.5 K, the model's near-perihelion dayside maximum discussed in Section 3.5, it gives roughly 9e19 cm^-2. Figure 5 reports maximum dayside columns of 3.7e16 and 8.5e16 cm^-2. To reach 8.5e16 cm^-2 from Eq. (8) as printed, T would have to be about 99 K, not the 110-116 K stated for the dayside. Thus either the code uses a prefactor near 1e13 Pa rather than the printed 1.1516e16 Pa, or the reported column densities and temperatures cannot both follow from Eq. (8). All of the paper's quantitative support for a purely sublimation-driven atmosphere—the latitudinal falloff, the sub-anti-Jovian asymmetry, the seasonal cycle—flows through this equation, so this is not a cosmetic typo. The Appendix A3 caveat about failure of equilibrium at eclipse ingress/egress and at terminators does not address this; the discrepancy occurs at the dayside maximum.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a time-dependent 1-D heat-diffusion model of Io's surface and subsurface temperatures, driven by solar insolation, Jovian thermal and reflected radiation, and internal tidal heating, with SPICE-based exact celestial geometry including Io's inclination and eclipses. Surface temperatures are converted to SO2 column densities under the assumption of instantaneous local vapor-pressure equilibrium with a uniform frost layer (Eq. 8). A parameter study over thermal diffusivity selects a default value alpha = 2.41e-7 m^2/s (thermal inertia 298 MKS) by matching observed diurnal temperature/column variations and the latitudinal column falloff. The resulting model reproduces an equator-to-pole column decrease of more than an order of magnitude, a sub-Jovian versus anti-Jovian asymmetry with dayside maxima of 3.7e16 and 8.5e16 cm^-2, a diurnal collapse of nearly two orders of magnitude, seasonal hemispheric asymmetries tied to Io's inclination, and requires an enhanced polar conductive heat flux of about 1.2-2.1 W/m^2 to explain the warm poles.","tokens_in":29745,"tokens_out":15156,"duration_ms":133959,"significance":"If the quantitative mapping in Eq. (8) is corrected, the model provides a simple, transparent, and falsifiable baseline for how much of Io's large-scale atmospheric structure can be explained by a purely sublimation-driven atmosphere. Strengths include the use of exact celestial geometry, a publicly available code, and an unusually candid Appendix A3 listing model limitations. The main quantitative claims, however, are conditional on the correctness of Eq. (8), which as printed cannot produce the reported column densities, and the headline latitudinal comparison is partly in-sample because the free thermal diffusivity is tuned against that same latitudinal profile. The paper is within the scope of JGR: Planets and offers a useful reference framework for future observations, provided the central equation and the framing of the parameter fit are fixed.","major_comments":[{"comment":"Equation (8) as printed is numerically inconsistent with the reported results. For the stated values m_SO2 = 1.064e-25 kg and g0 = 1.81 m/s^2, the equation gives N_sub(T=110 K) ~ 9e18 cm^-2 and N_sub(T=116.5 K) ~ 9e19 cm^-2, whereas the manuscript reports maximum dayside column densities of 3.7e16 and 8.5e16 cm^-2. The prefactor appears to be about 10^3 too large; replacing 1.1516e16 by 1.1516e13 Pa would yield values consistent with the reported columns at the stated dayside temperatures. Because Eq. (8) is the sole link between the simulated surface temperatures and every atmospheric column density in Sections 3.3-3.6, this is a load-bearing issue that must be corrected and verified against the public code output before publication.","section":"Section 2, Eq. (8); Figures 5 and 9"},{"comment":"The default thermal diffusivity alpha = 0.1 alpha_0 is selected by 'quantitative comparison' with the observed latitudinal column-density variation (Strobel and Wolven 2001) and with the observed diurnal temperature/column variation (Section 3.1 and 3.3). Therefore the agreement shown in Figures 5-7 for the equator-to-pole falloff and the diurnal collapse is a consistency check, not an independent prediction, and the abstract's wording 'We find that a thermal diffusivity ... yields ... decreasing ... in accordance with the observed spatial variations' should be reframed to state explicitly that this parameter was tuned to those observations. An out-of-sample test, for example fixing alpha using only the diurnal temperature data and then comparing to the latitudinal column profile, would materially strengthen the central claim.","section":"Section 3.3 and Abstract"},{"comment":"The central T-to-column mapping assumes instantaneous vapor-pressure equilibrium with a uniform, frost-covered surface everywhere. The authors candidly acknowledge in Appendix A3 that this equilibrium fails near eclipse ingress/egress and at terminators, but the paper's headline quantitative comparisons are made at the dayside maximum, where the equilibrium assumption is most favorable. Still, the model also assumes no winds, no lateral transport, and no partial frost coverage; each of these could change columns by factors of several. I ask the authors to add a short quantitative sensitivity statement, or at least a discussion of how partial frost coverage or finite sublimation/condensation rates would alter the reported column densities, so that the reader can gauge how robust the 'purely sublimation-driven' conclusion is to violations of the local-equilibrium assumption.","section":"Appendix A3 and Eq. (8)"}],"minor_comments":[{"comment":"Please define all symbols in Eq. (8) explicitly, including the units of the prefactor (Pa), the mass of the SO2 molecule in kg, and the explicit statement that T is in Kelvin; the current notation 'm_SO2[kg]·g0[m s^-2]' is confusing.","section":"Eq. (8)"},{"comment":"The parameter-study range is printed as '1.61−96.4×10^-7' which is ambiguous; use consistent scientific notation such as 1.61e-7 to 9.64e-6 m^2/s.","section":"Table 1"},{"comment":"The text says the model 'underestimates these column density values by a facor of 2 - 6' and cites '0.75−1.1×10^17 cm^-2'; there are typographical errors ('facor') and inconsistent spacing that should be corrected.","section":"Section 3.5.1"},{"comment":"The caption does not explain all symbol types (stars, diamonds, triangles) or their color coding; please add a legend or a complete caption so the multi-epoch observational comparisons are self-contained.","section":"Figure 9"},{"comment":"The comparison in Figure 6 is made separately for northern and southern hemispheres because of Io's inclination, but the text does not quantify the uncertainty in the modeled latitudinal profiles; adding a spread of alpha values within the stated 272-333 MKS range would help the reader judge whether the equator-to-pole decrease is robust within the fitted parameter range.","section":"Section 3.3 and Figure 6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well suited to JGR: Planets and the authors have provided public code and a candid limitations appendix, which I view as strengths. The critical issue is Eq. (8): as printed, it cannot produce the reported column densities, and the discrepancy is large enough that the paper's quantitative results are not reproducible without correction. I believe the prefactor is a typo rather than a fundamental error, because the reported temperatures and columns are mutually consistent with a prefactor of ~1e13 Pa, so I recommend major revision rather than rejection. I would also encourage the editor to ask the authors to verify the corrected equation against the released code output."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a serious, transparent paper that has a plausible physics story, but the printed version has a load-bearing numerical inconsistency. Eq. (8), the only link between surface temperature and sublimation column density, produces about 1e19-1e20 cm^-2 at the stated dayside temperatures, not the reported 3.7e16-8.5e16 cm^-2. To get 8.5e16 from the printed equation you'd need ~99 K. Either the code uses a prefactor near 1e13 Pa instead of the printed 1.1516e16, or the reported columns and temperatures don't come from the same model. I checked the stress-test note; the arithmetic holds. That has to be fixed before any quantitative claim can be trusted.\n\nWhat is genuinely new: the seasonal hemispheric asymmetry tied to Io's inclination (northern summer near perihelion, ~4 K hemispheric surface-temperature shift), the night-side contribution from Jupiter's radiation raising column densities by a factor ~3.5 when solar insolation vanishes, the 0.6 m depth-of-solar-influence estimate, and the polar conductive heat-flux constraints. The heat-equation and geometry parts don't depend on the Eq. (8) prefactor. The paper also does real work: exact SPICE geometry, public Python code, a transparent parameter study, and a candid limitations appendix. The comparison against Strobel and Wolven's Lyman-alpha latitudinal profile targets the right observable, and the citation pattern is fair, including the prior DSMC work of Walker and the long TEXES/IRTF time series.\n\nWhere the paper is soft: the main free parameter alpha is fitted to the same latitudinal and diurnal observations that are then presented as reproduced. That is in-sample. The authors say this explicitly in Section 3.3, but the abstract and conclusions don't carry the caveat. The \"purely sublimation\" framing is also stronger than the model's own results: the anti-Jovian comparison is short by a factor of 2-6, and Giles et al.'s best fit needs a volcanic component of 0.74e17 cm^-2. The authors acknowledge both, but the title and abstract overstate. The physical soft spot is the equilibrium assumption: local column density equals vapor-pressure equilibrium over uniform frost. Winds, patchy frost, and finite sublimation/condensation times would decouple local columns from local surface temperature. For dayside, global-scale structure it is a reasonable starting assumption, not a proven closure; the authors list it in A3 but it's the load-bearing model choice.\n\nThis is a paper for Io observers and planetary atmospheric modelers. It deserves a serious referee, not a desk reject. I would ask for the Eq. (8) prefactor to be reconciled with the code, for residual plots, and for explicit in-sample versus out-of-sample labeling. If the code matches the reported numbers, this is a good revision; if not, the central claim collapses.","headline":"A transparent one-parameter model for Io's sublimation atmosphere with a real seasonal prediction, but Eq. (8) as printed is off by ~10^3 and the headline numbers can't be reproduced until it's fixed.","tokens_in":30348,"tokens_out":5418,"would_cite":true,"duration_ms":52740,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper shows that a purely sublimation-driven atmosphere, with a single tuned thermal diffusivity value, can reproduce the observed latitudinal, longitudinal, and temporal variability of Io's SO2 atmosphere.","keywords":["Io","sublimation atmosphere","SO2 column density","thermal inertia","thermal diffusivity","eclipse collapse","seasonal variability","tidal heating"],"falsifier":"A spatially resolved, minute-cadence observation of Io's SO2 column at the dawn terminator during an eclipse ingress would decide the issue: the model requires the column to drop by more than an order of magnitude within about ten minutes, tracking the computed surface temperature. If the observed column either lags the surface-temperature drop by substantially more than ten minutes or remains high over frost regions that have already cooled, the instantaneous-equilibrium assumption of Equation (8) is falsified.","tokens_in":29229,"feed_emoji":"🪐","tokens_out":8951,"duration_ms":79651,"temperature":0.7,"pith_summary":"This paper claims that the global structure of Io's SO2 atmosphere—its latitude bands, day–night asymmetry, sub-Jovian/anti-Jovian difference, eclipse collapse, and seasonal cycle—can all be produced by solar-driven sublimation of surface frost alone, without needing volcanoes to shape the large-scale pattern. The claim matters because it would settle a long-standing debate: volcanic outgassing then adds only local perturbations, while sublimation sets where and when the atmosphere is densest. The model's only free parameter is the subsurface thermal diffusivity, which it fixes at $\\alpha = 2.41\\times10^{-7}\\,\\mathrm{m^2\\,s^{-1}}$ (thermal inertia 298 MKS) by fitting the observed diurnal temperature range and latitudinal column-density decrease. On that single setting the model matches observed dayside column densities of $3.7\\times10^{16}\\,\\mathrm{cm^{-2}}$ (sub-Jovian) and $8.5\\times10^{16}\\,\\mathrm{cm^{-2}}$ (anti-Jovian), an equator-to-pole drop of over an order of magnitude, and a column-density variation of nearly an order of magnitude over a Jovian year.","feed_headline":"Io's atmosphere matches a purely sublimation-driven model","feed_subtitle":"A single thermal-inertia value reproduces its latitude, longitude, and seasonal structure.","key_machinery":"The central object is a one-dimensional, time-dependent heat-diffusion equation for Io's subsurface temperature, $\\partial T/\\partial t = \\alpha\\,\\partial^2 T/\\partial z^2 + (\\rho c_p)^{-1}(P - L)$, with solar insolation, Jupiter's thermal and reflected radiation, internal tidal heating at the lower boundary, and surface radiation loss as source and sink terms, all evaluated with SPICE-derived exact celestial geometry. The free parameter is the thermal diffusivity $\\alpha$, which controls how much of the diurnal heating is stored in the subsurface and re-emitted at night. The second step maps surface temperature to column density through the vapor-pressure equilibrium formula $N_{\\mathrm{sub}} = 1.1516\\times10^{16}\\,\\exp(-4510/T)/(m_{\\mathrm{SO_2}} g_0)$ (Wagman, 1979), which is what makes the model testable against observations. This machinery converts orbital geometry and one material property into a global, time-varying map of the atmosphere.","core_discovery":"The paper argues that, on global scales, Io's SO2 atmosphere is a purely sublimation-driven atmosphere in the quantitative sense that most of its observed spatial and temporal variability is reproduced by mapping the local surface temperature, computed from a 1-D heat-diffusion equation with exact orbital geometry, to the vapor-pressure equilibrium column density of SO2. With the best-fit thermal diffusivity $\\alpha = 2.41\\times10^{-7}\\,\\mathrm{m^2\\,s^{-1}}$ ($\\Gamma = 298$ MKS), the model produces dayside maxima of $3.7\\times10^{16}\\,\\mathrm{cm^{-2}}$ on the sub-Jovian and $8.5\\times10^{16}\\,\\mathrm{cm^{-2}}$ on the anti-Jovian hemisphere, a day–night contrast of almost two orders of magnitude at the equator, and a latitudinal decrease of more than one order of magnitude that matches Lyman-$\\alpha$ observations. The anti-Jovian excess by a factor of about four is a direct geometric consequence: because of Io's locked rotation, the anti-Jovian hemisphere is never in Jupiter's shadow and receives about two more hours of sunlight each day. Seasonal effects appear as hemispheric summers, with northern summer near perihelion and northern winter near aphelion, and the anomalous warm poles require an additional conductive heat flux of at least 1.2 W/$m^{2}$ at high latitudes.","pith_inferences":["If the equilibrium mapping is correct, eclipse collapse and regrowth become a remote probe of Io's surface thermal diffusivity at centimeter depths, giving a testable prediction for future spacecraft or high-cadence ground-based observations.","The same surface-temperature-to-column mapping could be applied to other sublimation-dominated tenuous atmospheres, such as Triton's nitrogen or Pluto's, where the equilibrium assumption is shakier; the model's success on Io would make it a useful null hypothesis there.","The polar heat-flux requirement (1.2–2.1 W/m^2) is essentially independent of insolation and albedo at high latitudes, so a future measurement of Io's polar heat flow would either confirm the deep-mantle heating scenario or falsify the model's polar extension."],"forward_implications":["If the claim holds, volcanoes control only local and transient patches of Io's atmosphere; the global pattern of SO2 is set by insolation and thermal inertia.","The anti-Jovian hemisphere should be denser than the sub-Jovian by about a factor of four at all seasons, as a direct consequence of Io's eclipse geometry.","Io's atmosphere should collapse by roughly an order of magnitude within about ten minutes of Jupiter's eclipse ingress, and recover similarly at egress.","Over a Jovian year the equatorial column density should change by about a factor of seven, with the northern and southern hemispheres alternately warmer by about 4 K.","Io's polar temperatures of 75–90 K require an internal conductive heat flux that increases toward the poles (roughly 1.2 W/m^2 in summer, 2.1 W/m^2 in winter), which favors deep-mantle tidal heating over a shallow magma ocean."],"supporting_citations":[{"why":"Supplies the vapor-pressure equilibrium relation (Equation 8) that converts surface temperature to SO2 column density.","marker":"Wagman (1979)"},{"why":"Provides the reference thermal diffusivity, mass density, and specific heat capacity that set the model's baseline.","marker":"Leone et al. (2011)"},{"why":"Provides the Lyman-alpha latitudinal column-density decrease that the model is designed to reproduce.","marker":"Strobel and Wolven (2001)"},{"why":"Supplies the observed surface temperature maps, including the anomalously warm poles that drive the polar heat-flux analysis.","marker":"Rathbun et al. (2004)"},{"why":"Provides seasonal column-density variations from TEXES that constrain the allowed thermal inertia range and the Jovian-year cycle.","marker":"Tsang et al. (2012)"},{"why":"Provides an independently determined best-fit thermal inertia (200±50 MKS) for a frost-covered surface used to validate the default value.","marker":"A. C. Walker et al. (2012)"},{"why":"Supplies observed sub- and anti-Jovian hemispheric column densities and the atmosphere's latitudinal extension.","marker":"Feaga et al. (2009)"},{"why":"Provides ALMA eclipse observations that constrain the collapse/recovery timescale and the depth-dependent thermal inertia.","marker":"de Pater et al. (2020b)"},{"why":"Supplies the SPICE kernel toolkit that computes the exact celestial geometry used in all simulations.","marker":"Acton (1996)"}],"fun_headline_variants":["Sublimation alone explains Io's atmospheric variability","Io's thin air: frost sublimation drives its entire structure","One thermal parameter recreates Io's observed atmosphere","Sublimation-only model fits Io's atmosphere data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the atmosphere is everywhere in instantaneous vapor-pressure equilibrium with a uniform, frost-covered surface, so that local column density is a pure function of local surface temperature; if the frost is patchy, if winds move the gas, or if sublimation and condensation lag behind temperature changes, the mapping in Equation (8) breaks, and the authors concede it fails at eclipse ingress/egress and at the terminators.","fun_headline_variants_meta":{"raw":{"variants":["Sublimation alone explains Io's atmospheric variability","Io's thin air: frost sublimation drives its entire structure","One thermal parameter recreates Io's observed atmosphere","Sublimation-only model fits Io's atmosphere data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000853,"raw_usage":{"total_tokens":3810,"prompt_tokens":1154,"completion_tokens":2656,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":770,"completion_tokens_details":{"reasoning_tokens":2593}},"tokens_in":770,"tokens_out":2656,"duration_ms":22128,"temperature":1.0,"reasoning_tokens":2593,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:36:35.097759+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A spatially resolved, minute-cadence observation of Io's SO2 column at the dawn terminator during an eclipse ingress would decide the issue: the model requires the column to drop by more than an order of magnitude within about ten minutes, tracking the computed surface temperature. If the observed column either lags the surface-temperature drop by substantially more than ten minutes or remains high over frost regions that have already cooled, the instantaneous-equilibrium assumption of Equation (8) is falsified.","supporting_citations":[{"cited_title":"APACrefauthors \\ 1979","cited_arxiv_id":null,"evidence_quote":"Supplies the vapor-pressure equilibrium relation (Equation 8) that converts surface temperature to SO2 column density."},{"cited_title":", Wilson, L","cited_arxiv_id":null,"evidence_quote":"Provides the reference thermal diffusivity, mass density, and specific heat capacity that set the model's baseline."},{"cited_title":", Spencer , J R","cited_arxiv_id":null,"evidence_quote":"Provides seasonal column-density variations from TEXES that constrain the allowed thermal inertia range and the Jovian-year cycle."},{"cited_title":", McGrath, M","cited_arxiv_id":null,"evidence_quote":"Supplies observed sub- and anti-Jovian hemispheric column densities and the atmosphere's latitudinal extension."}],"review_version":1}