{"id":"1375a122-ba31-4300-add5-bc3135b762eb","arxiv_id":"2606.26042","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":2,"one_line_summary":"Top-heavy thermochemical geodynamo simulations show heterogeneous outer boundary heat flux reproduces observed longitudinal inclination anomaly variations, while elevated chemical driving reduces longitudinal structures and promotes polar radial field minima.","lead":"This paper runs top-heavy geodynamo simulations varying the balance of thermal and chemical buoyancy and compares the resulting time-averaged magnetic fields to paleomagnetic observations. The simulations indicate that heterogeneous mantle heat flux produces longitudinal inclination variations matching data, while chemical driving reduces those variations but creates polar field minima.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"Central claim hinges on unverified fidelity of top-heavy thermochemical setup to Earth's actual buoyancy sources and boundary conditions","rationale":"The reader's weakest_assumption is precisely the load-bearing point; the abstract alone supplies no additional internal checks (equations, tables, or convergence tests) that would allow a more technical objection. Therefore the UNVERDICTED verdict and low confidence remain appropriate.","tokens_in":1730,"tokens_out":335,"duration_ms":17582,"concrete_test":"Re-run the heterogeneous-boundary case at the reported chemical-driving strength using an independent geodynamo code (e.g., MagIC or Rayleigh) with identical non-dimensional parameters; if the time-averaged inclination anomaly loses its longitudinal variation or the amplitude differs by >30 % from the original suite, the claimed observational match is not robust to numerical implementation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline result—that heterogeneous outer-boundary heat flux produces longitudinal inclination-anomaly variations matching paleomagnetic constraints while homogeneous flux does not, and that elevated chemical driving only partially suppresses them—requires that the chosen ranges of chemical Rayleigh number, thermal driving strength, and imposed heat-flux heterogeneity correctly dominate the buoyancy budget and boundary conditions inside the outer core. The abstract provides no quantitative justification for these ranges, no resolution or benchmark tests, and no sensitivity analysis showing that the reported distinction survives changes in the relative strength of chemical versus thermal buoyancy or in the spatial pattern of the heterogeneous flux. If the model’s buoyancy partitioning or boundary implementation deviates from Earth, the contrast between heterogeneous and homogeneous cases could be an artifact rather than a robust signature.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper presents a suite of top-heavy thermochemical geodynamo simulations that vary the relative strength of chemical versus thermal buoyancy while imposing either homogeneous or heterogeneous outer-boundary heat flux. Spatially averaged field morphology metrics fail to distinguish the cases, but simulations with heterogeneous boundary forcing reproduce the longitudinal variations in time-averaged inclination anomaly required by paleomagnetic field models spanning centuries to tens of kyr; elevated chemical driving reduces but does not eliminate these structures and deepens polar radial-field minima. The authors conclude that both strong mantle heat-flux heterogeneity and chemical driving are likely present in Earth’s core.","tokens_in":1885,"tokens_out":498,"duration_ms":17690,"significance":"If the chosen buoyancy partitioning and boundary implementation are representative of the outer core, the work supplies a concrete mechanistic link between core-mantle thermal interactions and a specific, observationally testable paleomagnetic signature, moving beyond generic buoyancy assumptions common in earlier dynamo studies.","major_comments":[{"comment":"The central claim—that heterogeneous outer-boundary heat flux produces longitudinal inclination-anomaly variations matching observations while homogeneous flux does not—rests on the fidelity of the selected ranges of chemical Rayleigh number, thermal driving strength, and imposed heat-flux heterogeneity. No quantitative justification for these ranges, no resolution or benchmark tests, and no sensitivity sweeps demonstrating that the heterogeneous-versus-homogeneous distinction survives changes in buoyancy partitioning are provided.","section":"Methods / Simulation Setup"},{"comment":"The abstract states that observational constraints on longitudinal inclination-anomaly variations are 'readily matched' by the heterogeneous cases. Without explicit description of how the time-averaged inclination anomaly is extracted from the simulations, which paleomagnetic field models are used, data exclusion criteria, or error quantification in the comparison, it is impossible to judge whether the reported match is robust or an artifact of post-processing choices.","section":"Results / Comparison to Observations"}],"minor_comments":[{"comment":"Notation for the relative chemical driving strength should be defined once in the text and used consistently; the abstract introduces the concept without a symbol or equation reference.","section":null},{"comment":"Figure captions should state the exact time-averaging interval and number of realizations used for each reported field morphology metric.","section":null}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their constructive comments, which help clarify the presentation of our results. We respond to each major comment below.","responses":[{"response":"We agree that the manuscript would be strengthened by explicit justification and additional tests. The chosen ranges are motivated by estimates of chemical versus thermal buoyancy in Earth's core from prior literature, but we will add a new subsection in Methods providing quantitative rationale with references, along with resolution and benchmark tests. We will also include a sensitivity analysis (based on additional runs) showing the heterogeneous-versus-homogeneous distinction in longitudinal inclination anomalies persists across variations in buoyancy partitioning. These will be incorporated in the revised version.","revision_made":"yes","referee_comment":"[Methods / Simulation Setup] The central claim—that heterogeneous outer-boundary heat flux produces longitudinal inclination-anomaly variations matching observations while homogeneous flux does not—rests on the fidelity of the selected ranges of chemical Rayleigh number, thermal driving strength, and imposed heat-flux heterogeneity. No quantitative justification for these ranges, no resolution or benchmark tests, and no sensitivity sweeps demonstrating that the heterogeneous-versus-homogeneous distinction survives changes in buoyancy partitioning are provided."},{"response":"We acknowledge that the comparison procedure requires more detail for reproducibility. In the revised manuscript we will expand the relevant Results subsection to describe exactly how the time-averaged inclination anomaly is extracted (including averaging intervals and spatial binning), name the specific paleomagnetic field models used, state any exclusion criteria, and report quantitative metrics such as RMS misfit and correlation coefficients between simulated and observed longitudinal variations. This will permit direct evaluation of the robustness of the match.","revision_made":"yes","referee_comment":"[Results / Comparison to Observations] The abstract states that observational constraints on longitudinal inclination-anomaly variations are 'readily matched' by the heterogeneous cases. Without explicit description of how the time-averaged inclination anomaly is extracted from the simulations, which paleomagnetic field models are used, data exclusion criteria, or error quantification in the comparison, it is impossible to judge whether the reported match is robust or an artifact of post-processing choices."}],"tokens_in":1418,"tokens_out":453,"duration_ms":31910,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main point is that the new suite of top-heavy thermochemical geodynamo runs finds heterogeneous mantle heat flux reproduces the longitudinal variations in time-averaged inclination anomaly required by observational field models, whereas homogeneous flux does not. Elevated chemical driving reduces those structures without erasing them and also deepens polar minima in the radial field. Spatially averaged measures of morphology and variability do not separate the cases.\n\nWhat is new is the systematic variation of relative thermal versus chemical driving combined with direct tests against multi-timescale paleomagnetic constraints on that specific longitudinal feature. The paper does a clean job of isolating the boundary condition effect and showing it survives moderate chemical driving.\n\nThe soft spot is that the abstract gives no numbers on how close the matches are, no error quantification, and no sensitivity tests on the chosen ranges of chemical Rayleigh number or the spatial pattern of heat flux heterogeneity. If those ranges do not reflect the actual buoyancy budget in the core, the contrast between heterogeneous and homogeneous runs could shift. The stress-test concern about model fidelity is reasonable until the full methods and figures are checked.\n\nThis is for people working on core-mantle thermal coupling and paleomagnetic field morphology. A reader who wants concrete model predictions to compare with data will find it useful. It deserves a serious referee because the central comparison is grounded in independent observations and the setup is reproducible enough to test.\n\nRecommendation: send it to peer review.","headline":"These top-heavy simulations show heterogeneous outer-boundary heat flux produces longitudinal inclination anomalies matching paleomagnetic data while homogeneous cases do not, with chemical driving weakening but not removing the signal.","tokens_in":2389,"tokens_out":368,"would_cite":false,"duration_ms":17391,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Simulations with heterogeneous outer boundary heat flux match observed longitudinal variations in geomagnetic inclination anomaly, unlike homogeneous cases.","keywords":["geodynamo","paleomagnetism","core-mantle interactions","thermochemical convection","inclination anomaly","heterogeneous heat flux","geomagnetic field","chemical buoyancy"],"falsifier":"A paleomagnetic compilation covering multiple intervals that shows no longitudinal variations in time-averaged inclination anomaly, even though independent evidence indicates heterogeneous mantle heat flux, would falsify the reported match.","tokens_in":2639,"feed_emoji":"🧲","tokens_out":510,"duration_ms":30093,"temperature":0.7,"pith_summary":"This paper runs top-heavy thermochemical geodynamo simulations that vary the relative strength of thermal and chemical buoyancy to see what signatures appear in the time-averaged magnetic field. When these runs are compared against paleomagnetic field models spanning centuries to tens of thousands of years, spatially averaged measures of morphology and variability fail to separate different driving strengths or boundary conditions. Longitudinal variations in the time-averaged inclination anomaly, however, appear in the heterogeneous-heat-flux cases and are absent from the homogeneous cases. A sympathetic reader would care because the result ties surface magnetic records directly to the pattern of heat leaving the core and to the balance of chemical versus thermal buoyancy inside it.","feed_headline":"Heterogeneous heat flux matches geomagnetic inclination patterns","feed_subtitle":"Simulations show that only models with varied mantle heat flux reproduce observed longitudinal variations in time-averaged field anomaly.","key_machinery":"Top-heavy thermochemical geodynamo simulations with variable chemical driving strength and heterogeneous outer boundary heat flux, which separate equatorial thermal effects from polar chemical effects.","core_discovery":"Top-heavy geodynamo simulations that include heterogeneous outer boundary thermal forcing produce longitudinal variations in the time-averaged inclination anomaly that match observational constraints from paleomagnetic records, whereas simulations with homogeneous mantle heat flux do not. Increasing chemical driving reduces these longitudinal structures without erasing them and promotes deeper polar minima in the radial magnetic field. The results imply that Earth's outer core experiences both strong heat flux heterogeneity at the core-mantle boundary and significant chemical driving, resulting in small but persistent departures from the geocentric axial dipole approximation.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Mantle heat variations match inclination anomalies","Chemical buoyancy reduces field longitude variations","Thermal heterogeneity aligns with paleomagnetic records","Chemical driving deepens polar magnetic field minima"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The top-heavy thermochemical geodynamo model with the chosen ranges of chemical driving strength and boundary heat flux heterogeneity accurately represents the dominant buoyancy sources and boundary conditions in Earth's outer core.","fun_headline_variants_meta":{"raw":{"variants":["Mantle heat variations match inclination anomalies","Chemical buoyancy reduces field longitude variations","Thermal heterogeneity aligns with paleomagnetic records","Chemical driving deepens polar magnetic field minima"]},"model":"grok-4.3","cost_usd":0.00866,"raw_usage":{"total_tokens":3842,"prompt_tokens":702,"num_sources_used":0,"completion_tokens":50,"cost_in_usd_ticks":86603000,"prompt_tokens_details":{"text_tokens":702,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3090,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":702,"tokens_out":50,"duration_ms":30662,"temperature":1.0,"reasoning_tokens":3090,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-25T19:05:22.368687+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A paleomagnetic compilation covering multiple intervals that shows no longitudinal variations in time-averaged inclination anomaly, even though independent evidence indicates heterogeneous mantle heat flux, would falsify the reported match.","supporting_citations":[],"review_version":1}