{"id":"e0ce68c8-f0ad-44ac-a3aa-bbe4d19300f9","arxiv_id":"2508.12289","paper_version":4,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Horizontal convection with a water-like density maximum reorganizes into a single roll with full-depth plumes, and heat transport scales as Nu ~ Ra^{1/4} to Ra^{1/3}, faster than the classical Ra^{1/5} law.","lead":"This paper simulates a fluid heated and cooled along the same boundary and shows that water's density maximum near 4°C can change the whole flow pattern and let heat move much faster than classical theory predicts. The study matters for lakes, polar seas, and industrial models that assume water density changes smoothly with temperature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Nu-enhancement mechanism hinges on a Phi_i2 vs. z-hat scaling that is only 'suggested' and diagnosed from the same DNS runs it explains; until that relation is derived or tested out-of-sample, the claimed 1/4-to-1/3 exponent is not established.","rationale":"The reader's UNVERDICTED verdict was forced by the text mismatch; my concern is independent of that pipeline issue and targets the scientific argument as stated in the abstract. The weakest link is not the existence of the reorganization or the measured Nu slope—those are DNS findings—but the causal account: a single diagnostic quantity, plume height, is asserted to control the energy-budget term that produces the new scaling. Because the same DNS runs provide both z-hat and Nu, the agreement between the scaling model and the data is not an out-of-sample validation. I also note that 1/4-to-1/3 is a range of exponents, so the paper must show which exponent applies where and why; otherwise the 'enhanced scaling' claim is under-specified. The proposed test is feasible: Phi_i2 is an integral of known fields, and two additional Ra values or a shifted extremum temperature would break the calibration loop. I do not see grounds for rejection—the DNS evidence for flow reorganization and enhanced Nu could be correct—but acceptance should be conditional on the Phi_i2-z-hat relation being shown to be causal rather than diagnostic.","tokens_in":28055,"tokens_out":5203,"duration_ms":59264,"concrete_test":"Run two additional EXT-NELT DNS cases at Ra = 10^9 and 10^11 (or with the density-extremum temperature shifted by ±1 C), computing Phi_i2 directly from the stored fields as the nonlinear buoyancy-flux integral rather than as the residual of the OB dissipation balance. Then test whether the measured Phi_i2(z-hat) scaling predicts the DNS Nu exponent; if Phi_i2 deviates from the model's required Ra scaling even when z-hat ~ H, the mechanism is a restatement, not an explanation.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that in EXT-NELT the flow reorganizes to a single roll with full-depth plumes (z-hat ~ H), and that this changes the global kinetic-energy dissipation balance through the nonlinear-EOS term Phi_i2, raising Nu from the Rossby 1/5 scaling to 1/4-1/3. The load-bearing step is the assertion that Phi_i2 is controlled by z-hat. This step is not independently derived in the abstract: it rests on a 'scaling argument' whose key variable z-hat is diagnosed from the same DNS runs whose Nu(Ra) behavior the model then reproduces. If z-hat is extracted by post-processing the flow fields and Phi_i2 is inferred as the residual needed to close the OB energy budget, the agreement is partly built in. The quoted exponent range 1/4-to-1/3 also spans two different power laws, consistent with a non-asymptotic transition rather than a new dissipation closure. Since the full text with the Phi_i2 definition, equations, and DNS setup was not available, this mechanism could not be checked directly; as presented, the explanation is calibrated to the data rather than predictive.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript as identified by its title and abstract reports a two-dimensional direct numerical simulation study of horizontal convection with a nonlinear equation of state exhibiting a density extremum near 4°C. The abstract claims that in the EXT-NELT configuration the flow reorganizes from a bicellular structure into a single-roll circulation driven by full-depth mixing plumes, that heat transport is enhanced from the Rossby scaling Nu ~ Ra^{1/5} to Nu ~ Ra^{1/4}–Ra^{1/3}, and that this enhancement arises from an additional potential-energy transfer term Phi_{i2} introduced by the nonlinear equation of state, with the term's magnitude controlled by the characteristic plume height z-hat. The submitted full text, however, is not this paper: it is an unrelated manuscript on multi-modal quality analysis of weather radar forecasts (RadarQA, arXiv:2508.12291), containing no governing equations, simulation setup, energy-budget analysis, or numerical results for horizontal convection. The abstract-level claims are therefore the only reviewable content of the physics submission.","tokens_in":28192,"tokens_out":4421,"duration_ms":45053,"significance":"If the abstract's claims were fully substantiated, the result would be of genuine interest to the horizontal-convection and geophysical fluid dynamics communities: it would indicate that the standard Oberbeck-Boussinesq energy closure underpredicts heat transport in density-extremum configurations once plumes penetrate the full cavity depth, and it would identify a specific energy-budget term, Phi_{i2}, responsible for the change. The reported reorganization from a bicellular to a single-roll circulation is a plausible and potentially impactful observation for flows near the density maximum of water. However, the paper as submitted offers no machine-checked proofs, reproducible code, derivations, figures, or data tables for any of these claims; the actual manuscript body is a different work. The only quantitative statement available, that the model 'captures the main trends of the numerical data,' is too weak to constitute validation, and the claimed exponent range 1/4–1/3 is presented without fitting statistics or a stated Ra window.","major_comments":[{"comment":"The body of the submission is a different manuscript: the entire 'Full Text' section is the RadarQA paper on weather radar forecast quality analysis, with different authors, different subject matter, and no equations or results relating to horizontal convection, the density extremum, Phi_{i2}, z-hat, or Nu-Ra scaling. None of the central claims of the title and abstract can be checked because the supporting derivation, DNS setup, boundary conditions, equation-of-state form, and data-analysis procedure are absent. This is a load-bearing defect: the manuscript is internally inconsistent and cannot be reviewed as a physics paper.","section":"Full text"},{"comment":"The explanatory mechanism appears circular: the new potential-energy transfer term Phi_{i2} is asserted to be controlled by the characteristic plume height z-hat, and z-hat is diagnosed from the same EXT-NELT DNS runs whose heat-transport scaling the model then reproduces. The abstract states this only as a 'scaling argument' with 'suggests,' and no independent derivation of the Phi_{i2}(z-hat) relation is provided. As presented, the model is calibrated to the data it explains rather than offering an out-of-sample prediction; a test on a different Rayleigh-number range or an alternative equation-of-state parametrization would be needed to establish the claimed 1/4-to-1/3 exponent.","section":"Abstract, second paragraph"},{"comment":"The claim of 'enhanced heat transport scaling ranging from Nu ~ Ra^{1/4} to Nu ~ Ra^{1/3}' spans two distinct power-law exponents, and the abstract does not say whether this is a single fitted exponent that drifts with Rayleigh number, a crossover between two asymptotic regimes, or a fit with substantial uncertainty. Without the figures and fitting details that would appear in the full text, the reader cannot assess whether the data support a new dissipation closure or merely a transitional, non-asymptotic effect.","section":"Abstract, first paragraph"}],"minor_comments":[{"comment":"The term Phi_{i2} is introduced without a definition or an equation number, so its physical content is not accessible from the abstract alone.","section":"Abstract, second paragraph"},{"comment":"The abbreviations EXT, MON, LENT, and NELT are used without expansion, and the nonlinear equation-of-state form is never stated.","section":"Abstract, first paragraph"},{"comment":"The author list and reference list of the full text belong to the RadarQA manuscript, which is unrelated to the title and abstract of the submission; this mismatch should be resolved by the authors before any further review.","section":"Full text"}],"recommendation":"reject","confidential_remarks":"The submission contains a mismatched full text: the abstract describes a horizontal-convection physics study, while the body is the RadarQA weather-radar paper. This is not a fixable technical error within the paper's scope as submitted; it makes the manuscript non-reviewable. I recommend returning the submission to the authors to upload the correct full text. If the correct full text is provided, the circularity concern about z-hat raised in the stress-test note should be addressed directly by deriving Phi_{i2}(z-hat) or by performing an out-of-sample test, rather than only reporting that the model captures the main trends of the same data."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the practical problem: the \"full text\" attached to this record is not this paper. It's the RadarQA paper on weather-radar forecast quality analysis. So no governing equations, no grid resolution, no boundary conditions, no derivation of the Phi_i2 term, no actual Nu(Ra) data. Everything I say below is from the abstract alone.\n\nWhat the abstract describes is a plausible and genuinely new result for 2D horizontal convection: with a nonlinear equation of state and a 4°C density extremum (EXT-NELT), the large-scale flow reorganizes from two cells to a single roll with full-depth mixing plumes, and the heat-transport scaling moves from the classical Rossby Nu ~ Ra^{1/5} to somewhere between Ra^{1/4} and Ra^{1/3}. The proposed mechanism—an additional potential-energy transfer term Phi_i2 that modifies the OB kinetic-energy dissipation balance—is a real conceptual addition. The paper also does a clean-looking job of isolating the effects by crossing extremum/monotonic boundary conditions with linear/nonlinear EOS. The abstract is well written and honestly hedged: it says \"scaling argument suggests\" and \"possible energy budget interpretation,\" not \"we prove.\"\n\nThe soft spots are real. The load-bearing step is that Phi_i2 is controlled by the characteristic plume height z-hat, and z-hat is diagnosed from the same DNS runs whose Nu(Ra) the model reproduces. That makes the explanation partly calibrated to the data, and the abstract itself does not present an independent prediction or an out-of-sample test. The quoted exponent range is also uncomfortably broad: 1/4 to 1/3 spans two different power laws, which could just be a non-asymptotic transition rather than a new closure. And it's 2D, so the geophysical relevance is suggestive, not direct. These would be fixable concerns if the real manuscript provides a derivation or a prediction that goes beyond fitting.\n\nNet: this is exactly the kind of paper that deserves a serious referee—if the actual full text matches the abstract. The idea is subfield-scale but meaningful, and the mechanism is worth a careful look. But I cannot endorse the result on the current record, and I would not cite it until I've seen the real paper and checked whether the z-hat scaling is derived or just diagnosed.","headline":"The attached full text is a different paper (RadarQA), so this verdict rests entirely on the abstract: the density-extremum horizontal convection result is plausible and likely novel, but the central z-hat scaling argument is too under-supported to accept as established.","tokens_in":28891,"tokens_out":2588,"would_cite":false,"duration_ms":24704,"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":"In two-dimensional horizontal convection, combining a density extremum with a nonlinear equation of state reorganizes the flow into a single-roll circulation and, when mixing plumes reach full depth, enhances heat transport from the…","keywords":["horizontal convection","density extremum","nonlinear equation of state","mixing plumes","heat transport scaling","energy budget","direct numerical simulation","Rossby scaling"],"falsifier":"Run the EXT-NELT configuration in a taller cavity, or at a higher Rayleigh number, where $\\hat{z}$ no longer reaches $H$; if $Nu$ still follows a $Ra^{1/4}$–$Ra^{1/3}$ scaling, the plume-depth link fails. Conversely, a laboratory experiment with water near $4^\\circ\\mathrm{C}$ that directly measures plume penetration depth and Nusselt number would test whether $\\hat{z} \\sim H$ is necessary for the enhanced scaling.","tokens_in":27700,"feed_emoji":"🌊","tokens_out":5055,"duration_ms":50544,"temperature":0.7,"pith_summary":"This paper asks whether the standard Oberbeck-Boussinesq picture of horizontal convection survives when the fluid's density is not a linear function of temperature, as in water near $4^\\circ\\mathrm{C}$. Using two-dimensional direct numerical simulations from $Ra=10^6$ to $5\\times 10^{10}$, it compares linear and nonlinear equations of state with monotonic and density-extremum boundary conditions, and finds that only the extremum-plus-nonlinear case reorganizes the flow: a bicellular pattern gives way to a single roll driven by central mixing plumes. When those plumes penetrate the full cavity depth, heat transport follows a steeper scaling, $Nu \\sim Ra^{1/4}$ up to $Nu \\sim Ra^{1/3}$, instead of the usual Rossby $Nu \\sim Ra^{1/5}$. The paper traces this to an additional potential-energy transfer term in the total energy budget that the nonlinear equation of state introduces, whose size is set by the plume height. If correct, standard energy closures under-predict heat transport in density-extremum horizontal convection once plumes reach the full depth.","feed_headline":"Water's density extremum reshapes convection and speeds heat flow","feed_subtitle":"2D simulations show full-depth plumes push heat transport beyond Rossby's classic 1/5 power law.","key_machinery":"The load-bearing object is $\\Phi_{i2}$, an additional potential-energy transfer term appearing in the total energy budget when the equation of state is nonlinear, alongside the usual Oberbeck-Boussinesq terms. Its magnitude is argued to scale with the characteristic mixing-plume height $\\hat{z}$; full-depth plumes ($\\hat{z} \\sim H$) make the standard OB horizontal-convection kinetic-energy dissipation closure incomplete, which is how the paper moves the heat-transport exponent from $1/5$ to $1/4$–$1/3$.","core_discovery":"The central claim is that the nonlinearity of the equation of state, not the density-extremum boundary condition alone, is what changes horizontal convection. In the EXT-NELT configuration, the large-scale circulation shifts from two cells to a single roll, with central mixing plumes carrying the transport; when these plumes span the whole height $H$ ($\\hat{z} \\sim H$), the Nusselt number grows as $Ra^{1/4}$ to $Ra^{1/3}$, distinctly steeper than the Rossby $Ra^{1/5}$ seen in the other three configurations. The paper's mechanism is a new term $\\Phi_{i2}$ in the global kinetic-energy balance, a potential-energy transfer produced by the nonlinear equation of state, whose magnitude is controlled by $\\hat{z}$. With $\\hat{z} \\sim H$, this term changes the dissipation closure, and a scaling model built from it reproduces the main trends of the DNS data.","pith_inferences":["If the plume-height control holds, the transition from $Ra^{1/5}$ to $Ra^{1/3}$ should be continuous in $\\hat{z}/H$; a testable prediction is that the effective Nusselt exponent is a monotone function of plume penetration depth.","The single-roll reorganization suggests that three-dimensional simulations may show a qualitatively different large-scale flow, since two-dimensional central plumes are often sensitive to confinement; the scaling claim should be checked in three dimensions.","One can use the $\\Phi_{i2}$ budget term to design controlled experiments: changing the temperature of the cold boundary relative to $4^\\circ\\mathrm{C}$ should tune $\\hat{z}$ and produce a predictable shift in the Nusselt scaling.","Analogous potential-energy transfer terms may appear for other non-Oberbeck-Boussinesq nonlinearities, such as salinity or compressibility effects, so the same energy-budget analysis could identify enhanced transport in those settings."],"forward_implications":["For density-extremum horizontal convection with full-depth plumes, the standard Oberbeck-Boussinesq closure under-predicts heat transport; global energy budgets for such flows should include $\\Phi_{i2}$.","The flow structure itself changes: the bicellular circulation becomes a single-roll, central-plume regime, with transitional anomalies in the Reynolds-number scaling.","The heat-transport exponent is not fixed: it can lie between $1/4$ and $1/3$ depending on plume depth, so a single power law may not describe the full parameter range.","Models of lakes, ice-ocean settings, or industrial systems where water near $4^\\circ\\mathrm{C}$ is the working fluid should not assume the Rossby $1/5$ scaling once plumes become depth-penetrating.","The energy-budget model offers a diagnostic route: computing $\\Phi_{i2}$ and $\\hat{z}$ from DNS or experiments could predict when enhanced transport begins."],"supporting_citations":[],"fun_headline_variants":["Nonlinear water density spurs faster convection heat transport","Density extremum rewires convection, speeds heat flow","Beyond Rossby scaling: nonlinear density boosts convection","Mixing plumes drive single roll, boost heat transport in water","Nonlinear equation of state reshapes convection, speeds heat"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The explanation assumes that the extra energy term $\\Phi_{i2}$ is controlled by the plume height $\\hat{z}$ measured from the same simulations whose heat-transport scaling it is then used to reproduce, so the mechanism is calibrated to the data rather than independently predicting them.","fun_headline_variants_meta":{"raw":{"variants":["Nonlinear water density spurs faster convection heat transport","Density extremum rewires convection, speeds heat flow","Beyond Rossby scaling: nonlinear density boosts convection","Mixing plumes drive single roll, boost heat transport in water","Nonlinear equation of state reshapes convection, speeds heat"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000824,"raw_usage":{"total_tokens":3684,"prompt_tokens":1104,"completion_tokens":2580,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":720,"completion_tokens_details":{"reasoning_tokens":2500}},"tokens_in":720,"tokens_out":2580,"duration_ms":22188,"temperature":1.0,"reasoning_tokens":2500,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:25:12.704989+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the EXT-NELT configuration in a taller cavity, or at a higher Rayleigh number, where $\\hat{z}$ no longer reaches $H$; if $Nu$ still follows a $Ra^{1/4}$–$Ra^{1/3}$ scaling, the plume-depth link fails. Conversely, a laboratory experiment with water near $4^\\circ\\mathrm{C}$ that directly measures plume penetration depth and Nusselt number would test whether $\\hat{z} \\sim H$ is necessary for the enhanced scaling.","supporting_citations":[],"review_version":1}