{"id":"7cac3eef-6f81-4b32-a559-88013bf24632","arxiv_id":"2507.01441","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A pebble-based thermophysical model reproduces roughly the water outgassing of comet 67P through southern-hemisphere crust blow-off, but emits too much dust, CO2, and CO.","lead":"This paper tests a model in which gas pressure builds up under a comet's surface and blows off pebble-sized dust, using Rosetta data from comet 67P. It finds that the mechanism explains southern outgassing but overshoots observed dust and CO2 production, placing tight constraints on cometary material properties.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central mechanism requires b at or below the lower edge of the DSMC-derived range; preferred b=0.1Dp is below that range, so the activity claim hinges on an unvalidated microstructure.","rationale":"The reader's weakest assumption is indeed the low diffusivity. I agree this is the most load-bearing: unlike other limitations (overproduction, missing small dust, no northern activity), a shift in b within the quoted DSMC range changes the mechanism from active to inactive, not just from over- to under-producing. The paper's own numbers bracket the threshold tightly: activity at b = 0.3 Dp, cessation around 0.6 Dp. Since the literature range is 0.3–3 Dp, choosing 0.3 Dp selects the extreme lower boundary, and the preference for 0.1 Dp lies outside the range. The Fulle & Blum argument is suggestive but unquantified. I therefore keep the reader's CONDITIONAL verdict: the claim is plausible and internally consistent, but it rests on an unresolved microphysical parameter. A DSMC measurement of b for the actual proposed structure would settle the question. I do not escalate to REJECT because the mechanism is not contradicted by any internal inconsistency, and b = 0.3 Dp is within the cited range. My additional point about the preferred b = 0.1 Dp being below the range strengthens rather than changes the reader's concern.","tokens_in":21772,"tokens_out":8694,"duration_ms":105443,"concrete_test":"Run DSMC simulations of gas flow through the specific Fulle & Blum (2017) structure—random packings of 1 cm pebbles with fractal dust filling the interstices at the model's 0.76 total porosity—over several realizations, and measure b. Then rerun the fixed-heat-capacity model with the measured b and also with b = 0.4 Dp as a sensitivity case. If the measured b is ≥ 0.3 Dp, or if b = 0.4 Dp suppresses repeating southern-hemisphere ejections, the load-bearing assumption fails; if measured b ≈ 0.1 Dp and b = 0.4 Dp still activates the south, the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 defines permeability via the half-transmission thickness b and cites a theory/DSMC range b ≈ 0.3–3 Dp for effective porosities 0.4–0.8, with large scatter. The nominal model is run at b = 0.3 Dp, and Section 3.3 reports that for b ≳ 0.6 Dp only sporadic ejections occur and for b = 1 Dp there is no repeating activity. Section 5 further states a preference for b = 0.1 Dp, which is below the cited physical range. The only argument for allowing b < 0.3 Dp is the qualitative Fulle & Blum fractal-dust-filling narrative in Section 4.1, but no DSMC or laboratory measurement of b for that specific microstructure is provided. Thus the central mechanism operates only at the very edge of, or outside, the parameter range currently supported by gas-transport calculations. If the true b for the pebble-plus-fractal-dust structure is 0.4–0.5 Dp, within the quoted scatter, pressure build-up becomes insufficient and the entire ejecting-crust mechanism would not operate. The quantitative support is also strained: in the preferred b = 0.1 Dp, fixed-heat-capacity run (Table 1), CO2 is overproduced by a factor of 5.1 and dust by 7.7 relative to Rosetta estimates, so the 'roughly match' claim in the strongest claim is not a robust quantitative result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses a one-dimensional, pebble-based thermophysical model with gas-pressure build-up to simulate dust and volatile outgassing from comet 67P over three full orbits, with 19 latitude bins and separate H2O, CO2, and CO ice species. The authors show that for low gas diffusivity (half-transmission thickness b=0.3Dp or 0.1Dp) and favorable choices of dust-to-ice ratio, heat capacity, and tensile strength, southern latitudes undergo repeated ejection of pebble-to-decimetre-sized chunks, while northern latitudes drain volatiles deep and remain inactive. The resulting water production is in broad temporal agreement with Rosetta and produces a southern-dominated outgassing pattern qualitatively consistent with non-gravitational acceleration and torque modelling. However, the same runs produce total CO2 and dust masses 5–18 times the Rosetta estimates in the nominal cases, and the preferred low-diffusivity regime lies at or below the lower edge of the DSMC-derived range, so the quantitative claim of 'roughly matching' Rosetta is not yet established.","tokens_in":22103,"tokens_out":7139,"duration_ms":83628,"significance":"If the mechanism is ultimately validated, it would resolve the cohesion bottleneck and provide a natural explanation for 67P's southern 'blow-off' water emission pattern inferred from dynamics, making it an important contribution to cometary activity modelling. The paper has clear strengths: a full-orbit multi-latitude treatment, inclusion of three volatiles, systematic sensitivity tests over ice fractions, heat capacity, diffusivity, and strength, and comparison with multiple Rosetta datasets. It also reports numerical resolution checks and explicitly identifies the regimes where the model fails (northern activity, small dust, pre-perihelion CO/CO2). The central significance is currently conditional because the activity depends on a diffusivity parameter that is not independently validated and because the claimed match to observed dust and CO2 masses is not within a factor of ~2.","major_comments":[{"comment":"The conclusion that the model produces outgassing rates that 'roughly match those observed by Rosetta' is not supported by the totals in Table 1: the preferred fixed-c, b=0.1Dp run gives ΔM_CO2/observed=5.1 and ΔM_dust/observed=7.7, while the nominal b=0.3Dp run gives 10.7, 3.4, and 18.2 for CO2, CO, and dust, respectively. Since the abstract and conclusions present the mechanism as reproducing the global emission rates, this discrepancy is load-bearing for the paper's main claim; please either restrict the 'match' claim to water, or demonstrate a physical mechanism (e.g., partial fallback, reduced active area, or an ice-dependent strength law) that brings the other species within a factor of ~2 without destroying the water match.","section":"Section 5 and Table 1"},{"comment":"The mechanism operates only in a diffusivity regime that is at or below the low end of the quoted physical range. The paper states b≈0.3–3Dp for effective porosities 0.4–0.8, takes b=0.3Dp as the nominal value, reports no repeating activity for b≳0.6Dp, and prefers b=0.1Dp. The Fulle & Blum fractal-dust argument in Section 4.1 is qualitative and does not provide a DSMC or laboratory measurement of b for that specific microstructure. Given that a b value of 0.4–0.5Dp lies within the quoted scatter and would shut off the activity in this model, the central claim depends on an unvalidated parameter estimate; a direct DSMC or experimental determination of b for the pebble-plus-fractal-dust structure (or an observational constraint from Rosetta data) is needed.","section":"Sections 2 and 3.3 and 4.1"},{"comment":"The dust-mass results rest on the assumption that ejected particle size equals the depth of the ejection layer, but the paper offers no physical derivation for this relation. Since this assumption sets both the ejected mass and the size distribution, its sensitivity should be explored; otherwise the factor-of-seven dust overproduction in the preferred run cannot be reliably attributed to the physics of the ejecting crust.","section":"Section 3.4 and Figure 4"},{"comment":"The abstract claims reproduction of 'global emission rates of dust,' but the model only ejects particles at or above pebble size and cannot generate the ubiquitous sub-millimetre dust coma observed by OSIRIS. This is a model limitation acknowledged in Section 4.2, but it means the claim should be explicitly limited to the pebble/chunk component, and the implications for the 'dominant mechanism' claim should be discussed.","section":"Section 4.2 and Abstract"}],"minor_comments":[{"comment":"The supplied manuscript text has many missing spaces between words and some garbled inline notation (e.g., the first paragraph of the Introduction); please ensure the final submission is properly typeset.","section":"Throughout"},{"comment":"Table 1 does not include the water-content-dependent strength run discussed in Section 4.1; adding its ΔM_H2O and ΔM_CO values would make the comparison of all species possible.","section":"Table 1 and Section 4.1"},{"comment":"The phrase 'Strong constraints are placed' overstates the status of the parameters, which are varied to best match the Rosetta data rather than independently measured; consider rephrasing to 'best-fit parameters' or discussing the resulting degeneracies.","section":"Abstract and Section 5"},{"comment":"The data availability statement says data will be shared 'on reasonable request'; given the many free parameters and fitted curves, making the simulation code or tabulated output available would strengthen reproducibility.","section":"Data Availability"},{"comment":"Figure 1 would benefit from indicating the temperature range over which 67P's thermal inertia has actually been measured, so the reader can judge which pebble size and heat-capacity prescription is most relevant.","section":"Figure 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is squarely within scope for a cometary physics journal. The main risk is that the preferred b=0.1Dp lies below the DSMC-derived range quoted in the paper; if the authors can obtain or cite a direct calculation or measurement for the fractal-dust microstructure, the manuscript would be much stronger. The abstract and conclusions should also be made quantitatively consistent with Table 1. I do not see a novelty-disclosure issue; the work builds on prior papers in a transparent way."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it is a genuine step forward for the ejecting-crust framework: full orbit, 19 latitudes, CO added as a separate species, and a clear demonstration that the southern 'blow-off' can reproduce the spatial pattern of water outgassing inferred from non-gravitational accelerations and torques. Second, the fit to the data is not as good as the abstract implies. The nominal model exceeds observed CO2 by 10.7x, CO by 3.4x, and dust by 18.2x; even the best-fitting case in Table 1 leaves dust at 7.7x and CO2 at 5.1x. The authors are honest about this in the text, but the abstract's 'roughly match' is doing a lot of work.\n\nThe main load-bearing soft spot is the gas diffusivity. The model requires a half-transmission thickness b below about 0.3 Dp, with a stated preference for 0.1 Dp, while the cited DSMC and theoretical range is 0.3–3 Dp. The paper justifies the low end with the Fulle & Blum fractal-dust narrative, but gives no lab or DSMC measurement of b for that microstructure. If the true value is 0.4–0.5 Dp, within the quoted scatter, pressure build-up is insufficient and the ejecting-crust mechanism stops working entirely. That is a real fragility, not a nitpick.\n\nGive credit where it's due. The parameter exploration is systematic, the authors test heat capacity, diffusivity, ice fractions, strength laws, and numerical resolution. They explicitly flag the missing small-dust component and the failure to generate northern activity. They also discuss fallback and re-condensation as possible remedies. The paper is transparent about its own limitations to an unusual degree.\n\nThe 'constraints' on diffusivity, heat capacity, and tensile strength are really calibration results: parameters are adjusted until the output best matches the same Rosetta data used for comparison. That is not circular in a damning way, but it means the title's 'constraints' should be read as 'model requirements' rather than independent measurements. No code or data are released, which makes the calibration hard to audit.\n\nWho is this for? Anyone working on cometary thermophysical models or the Rosetta dataset. The southern blow-off mechanism is an interesting and plausible explanation for the NGA/torque results, and the paper lays out clear observational tests (e.g., OSIRIS images to distinguish homogeneous small-dust ejection from localized WEBs). It deserves a serious referee, but the referee should push hard on the diffusivity range and ask for a much more careful abstract.\n\nMy recommendation: send it to review, but expect major revisions. The honest text already acknowledges most of what a referee would say; the authors need to either bring the diffusivity into the validated range, find a way to measure b for fractal dust, or reframe the paper as a scenario study rather than a constrained model.","headline":"A serious, transparent attempt to make the ejecting-crust model work for all of 67P, but the mechanism only operates at the edge of the physically plausible diffusivity range and the abstract overstates the fit quality.","tokens_in":22714,"tokens_out":1609,"would_cite":true,"duration_ms":20799,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Subsurface gas pressure can eject centimetre-sized chunks and drive comet 67P's activity.","keywords":["comet 67P","thermophysical model","gas pressure build-up","dust ejection","tensile strength","gas diffusivity","outgassing","pebble structure"],"falsifier":"Measure the gas permeability of a realistic pebble packing with fractal dust filling the interstices—the structure assumed for 67P—using direct-simulation Monte Carlo or laboratory flow experiments. If the resulting half-transmission thickness comes out near one pebble diameter or larger, the pressure build-up that powers the model cannot occur; the paper itself reports that $b = 1\\,D_p$ yields only sporadic ejections and no repeating activity.","tokens_in":21526,"feed_emoji":"☄️","tokens_out":9262,"duration_ms":96448,"temperature":0.7,"pith_summary":"Cometary activity has long faced what is called the cohesion bottleneck: sublimating gas should be too weak to lift dust off a nucleus, yet comets vigorously shed dust. This paper argues that a pebble-structured surface with very low gas permeability can break the bottleneck, because sublimation pressure builds up under a thin crust and ejects pebbles and chunks a few millimetres to about a decimetre in size. Simulating 19 latitudes across comet 67P over three orbits, the authors obtain water, CO2, and CO emission rates that roughly match Rosetta's observations, with essentially all ejections coming from the southern hemisphere during perihelion. They argue this southern 'blow-off' naturally explains the strong southern water outgassing previously inferred from non-gravitational acceleration and torque modelling. The mechanism only works if the subsurface is far less gas-permeable than standard estimates suggest, and even then the model tends to overproduce dust, CO2, and CO.","feed_headline":"Gas pressure can blow off comet 67P's crust","feed_subtitle":"A pebble-scale model reproduces Rosetta's outgassing rates, but only if the subsurface is nearly gas-tight.","key_machinery":"The load-bearing object is the half-transmission thickness $b$ (the diffusion-scale parameter), defined as the number of particle layers that reduces the outgassing flux by half; it controls how easily sublimating gas escapes through the dust crust. For gas flow between pebbles, $b$ is proportional to pebble diameter $D_p$, with theory and direct-simulation Monte Carlo giving $b \\approx 0.3$–$3\\,D_p$. The paper shows that only the low end ($b \\le 0.3\\,D_p$, best fit $0.1\\,D_p$) lets pressure accumulate enough to overcome the layer strength, i.e., the sum of gravitational pressure and a depth- and size-dependent tensile strength. When pressure wins, one or more pebble layers are ejected, resetting a thin dust crust and keeping volatile ices within the top centimetre in the active south. Nearly all subsequent behaviour—the outgassing rates, the latitudinal asymmetry, and the parameter constraints—follows from the frequency and timing of these ejection events.","core_discovery":"The paper's central claim is that the ejecting-crust mechanism—gas pressure in the interstitial space between pebbles exceeding the low tensile strength of the pebble aggregate—can account for the global activity of comet 67P, provided the gas diffusivity is low. With a half-transmission thickness of $b = 0.3\\,D_p$ or less (best at $0.1\\,D_p$), repeated ejection cycles are established in the southern hemisphere near perihelion, keeping water ice within the top centimetre and letting CO2 actually reach the surface. The modelled total water emission is 1.6 times the observed value in the nominal run, while CO2, CO, and dust totals exceed observations by about 10.7, 3.4, and 18 times, respectively; the shape of the southern activation and the post-perihelion decline match well. The resulting active fraction is a few percent in the north and rises to roughly 20–30 percent in the south, matching the pattern inferred from the comet's non-gravitational acceleration and torque. The paper concludes that low gas diffusivity, large heat capacity, and a steeply depth- or ice-content-dependent tensile strength are required for the mechanism to work, and that the location and nature of erosion is the critical unknown for cometary activity.","pith_inferences":["One consequence the authors leave implicit: on any other comet with the same pebble fabric, the same blow-off mechanism should concentrate activity on whichever hemisphere is strongly illuminated near perihelion, so the latitudinal activity pattern of other Jupiter-family comets would be a direct test of the model.","The difficulty of suppressing CO2-driven ejection suggests an upper bound on the CO2 abundance of the non-water-enriched surface material: if the CO2-to-water ice fraction were much above roughly ten percent, the southern hemisphere would erode too rapidly and emit too much CO2, so the paper indirectly tightens constraints on 67P's primordial volatile inventory.","Because the model cannot produce the continuous small-dust coma seen at all latitudes, the real activity mechanism may be bimodal—large pebble and chunk ejection during southern summer plus a distributed small-dust removal process—and re-analysing coma images for whether small dust comes from localized patches or the whole disk would discriminate between these."],"forward_implications":["If the paper's central claim is right, the bulk of 67P's outgassing of all species, and essentially all ejection of pebbles and chunks, is confined to the southern hemisphere during the few months around perihelion.","The model reproduces the pattern inferred from non-gravitational acceleration and torque fits: water emission from the north is weak and roughly constant, while southern emission rises sharply at perihelion, with effective active fractions of a few percent in the north and about 20–30 percent in the south.","Ejected material in the model is limited to pebbles and chunks from a few millimetres up to about 15 centimetres, with a power-law size distribution that is shallower than small-particle observations but steeper than some large-chunk estimates.","Matching the Rosetta data requires a narrow parameter window—dust-to-ice mass ratio near 2, CO2-to-water fraction around 0.03–0.1, CO fraction near 0.01, low diffusivity, and either high heat capacity or a steep tensile-strength gradient—so the model constrains the subsurface structure as much as it explains the activity."],"supporting_citations":[{"why":"Supplies the pebble-based thermophysical model, the sublimation-rate and permeability equations, and the outgassing-area factor that the present simulations use.","marker":"Gundlach et al. 2020"},{"why":"Extended the model to a full cometary orbit on a single surface patch; its difficulty generating ejections motivates the multi-latitude, low-diffusivity setup here.","marker":"Bischoff et al. 2023"},{"why":"Previous development of the model showing that low gas diffusivity is essential and that internal-pebble pressure alone does not explain bulk activity; this paper continues that model for inter-pebble pressure.","marker":"Attree et al. 2024a"},{"why":"Provides the depth/size-dependent tensile-strength relation used to decide when gas pressure ejects a layer.","marker":"Skorov & Blum 2012"},{"why":"Comet-nucleus model of densely packed pebbles with fractal dust filling voids that supports the low gas-diffusivity assumption.","marker":"Fulle & Blum 2017"},{"why":"Rosetta water-production dataset the modelled water curve is compared against.","marker":"Läuter et al. 2020"},{"why":"Rosetta CO2 and CO production data used to normalize total volatile mass-loss and compare emission curves.","marker":"Läuter et al. 2019"},{"why":"Updated dust mass-loss estimate used to normalize and compare modelled dust ejection totals.","marker":"Laurent-Varin et al. 2024"},{"why":"Non-gravitational acceleration and torque modelling that inferred strong southern water emission, which the blow-off mechanism naturally explains.","marker":"Attree et al. 2024b"},{"why":"Theoretical and DSMC-based relation between half-transmission thickness b and particle size D_p that fixes the diffusivity range tested.","marker":"Macher et al. 2024"}],"fun_headline_variants":["Gas pressure can blow off comet 67P's crust in south","Comet 67P's southern blow-off explained by gas pressure","Pebble gas-tightness drives comet 67P's eruptions","Tight subsurface propels comet 67P's dust ejection","Comet 67P's activity pinned to low-diffusivity crust"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The subsurface must be much less gas-permeable than the microphysical simulations suggest: the half-transmission thickness has to be at or below about 0.3 pebble diameters (best around 0.1), whereas the literature-derived range spans 0.3 to 3 pebble diameters; if the real value is near 1, the model produces only sporadic ejections and no repeating activity.","fun_headline_variants_meta":{"raw":{"variants":["Gas pressure can blow off comet 67P's crust in south","Comet 67P's southern blow-off explained by gas pressure","Pebble gas-tightness drives comet 67P's eruptions","Tight subsurface propels comet 67P's dust ejection","Comet 67P's activity pinned to low-diffusivity crust"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001213,"raw_usage":{"total_tokens":5103,"prompt_tokens":1168,"completion_tokens":3935,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":784,"completion_tokens_details":{"reasoning_tokens":3842}},"tokens_in":784,"tokens_out":3935,"duration_ms":35651,"temperature":1.0,"reasoning_tokens":3842,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:51:00.845791+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the gas permeability of a realistic pebble packing with fractal dust filling the interstices—the structure assumed for 67P—using direct-simulation Monte Carlo or laboratory flow experiments. If the resulting half-transmission thickness comes out near one pebble diameter or larger, the pressure build-up that powers the model cannot occur; the paper itself reports that $b = 1\\,D_p$ yields only sporadic ejections and no repeating activity.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the pebble-based thermophysical model, the sublimation-rate and permeability equations, and the outgassing-area factor that the present simulations use."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extended the model to a full cometary orbit on a single surface patch; its difficulty generating ejections motivates the multi-latitude, low-diffusivity setup here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the depth/size-dependent tensile-strength relation used to decide when gas pressure ejects a layer."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Updated dust mass-loss estimate used to normalize and compare modelled dust ejection totals."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Theoretical and DSMC-based relation between half-transmission thickness b and particle size D_p that fixes the diffusivity range tested."}],"review_version":1}