{"id":"8350756c-1e45-4a43-9019-9fdd38717504","arxiv_id":"2606.10836","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"Introduces an OBBT framework that discretizes Darcy's law via finite volumes, applies McCormick relaxations with flow sign and irrotationality constraints, and yields guaranteed outer bounds on all variables for non-identifiable groundwater models.","lead":"The paper proposes using optimization-based bound tightening to compute guaranteed interval bounds on uncertain parameters and states in non-identifiable groundwater models, instead of relying on sampled ensembles. This deterministic approach could reduce the risk of underestimating uncertainty ranges in subsurface flow predictions.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Validity and tightness of bounds after adding flow sign prescription and irrotationality constraints to McCormick relaxations of the FV discretization","rationale":"The reader's weakest_assumption pinpoints exactly the same technical hinge. Because the review was performed on the abstract alone, the soundness of the remedies cannot be checked; the full manuscript would be required to perform the concrete test above. No other internal inconsistency (e.g., in the OBBT formulation itself) is visible from the given material.","tokens_in":1776,"tokens_out":339,"duration_ms":23152,"concrete_test":"From the full text, extract the precise algebraic form of the irrotationality constraint (discrete curl) and the flow sign prescription; confirm that every solution of the original nonlinear Darcy's law satisfies them. Then re-optimize the 2D steady-state example both with and without the remedies and measure the change in bound width on at least one state variable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (guaranteed outer bounds without sampling) requires that the McCormick-relaxed finite-volume system, once augmented by the two remedies, remains a valid outer approximation whose optima are both sound (original feasible set is contained) and not excessively loose. The abstract explicitly states that plain McCormick breaks sign coupling between fluxes and head gradients, allowing non-physical rotational flow; the remedies are therefore load-bearing. If either remedy is not a valid inequality for the original nonlinear problem, or if the residual relaxation gap stays large on the target models, the computed intervals cease to be reliable outer bounds on the physically admissible null space.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript proposes an interval-based uncertainty quantification method for non-identifiable groundwater models via Optimization-based Bound Tightening (OBBT). It discretizes Darcy's law with a finite-volume scheme, relaxes the resulting bilinear terms using McCormick envelopes, and augments the formulation with flow sign prescription and irrotationality constraints to restore physical consistency. The framework is demonstrated on a 1D steady-state model, a 2D steady-state model across four configurations, and a 2D transient model on a hexagonal grid, with the central claim being guaranteed outer bounds on all uncertain variables without sampling.","tokens_in":1904,"tokens_out":480,"duration_ms":23690,"significance":"If the added constraints produce a valid outer approximation whose optima remain both sound and sufficiently tight, the approach supplies a deterministic, conservative alternative to ensemble UQ that directly encodes the null space of admissible solutions. This would be a useful complement to sampling-based methods in data assimilation contexts, particularly where exhaustive exploration is impractical.","major_comments":[{"comment":"Abstract (paragraph on remedies): The assertion that flow sign prescription and irrotationality constraints are effective remedies is load-bearing for the guarantee of outer bounds. The manuscript must show that these inequalities are valid for the original nonlinear problem (i.e., every feasible point of the unrelaxed finite-volume system remains feasible after augmentation) and does not merely tighten the relaxation gap at the cost of excluding admissible solutions.","section":"Abstract (remedies paragraph)"},{"comment":"Numerical examples section: The demonstrations on the three models report no quantitative metrics of bound tightness (e.g., ratio of obtained interval widths to independently computed ground-truth ranges) or residual relaxation gap. Without such measures or comparison against a reference method, it is not possible to confirm that the remedies leave the bounds both valid and practically useful rather than excessively loose.","section":"Numerical examples"}],"minor_comments":[{"comment":"The abstract mentions characterizing the 'strengths and limitations' of the two remedies; the main text should supply explicit statements of these (e.g., which remedy is more effective on transient vs. steady problems) with supporting numerical evidence.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their thoughtful review and constructive feedback on our manuscript. We address each major comment below, indicating revisions where appropriate to strengthen the presentation of validity and tightness.","responses":[{"response":"We agree that an explicit demonstration of validity is required to support the outer-bound guarantee. The original manuscript characterizes the strengths and limitations of these constraints but does not contain a dedicated proof. In the revised version we will insert a new subsection (Section 3.4) that derives the flow-sign and irrotationality inequalities directly from the continuous Darcy and continuity equations, shows that every solution of the unrelaxed finite-volume system satisfies them, and therefore confirms that the augmented relaxation remains a valid outer approximation. This addition will be accompanied by a short remark on the conditions under which each constraint is active.","revision_made":"yes","referee_comment":"Abstract (paragraph on remedies): The assertion that flow sign prescription and irrotationality constraints are effective remedies is load-bearing for the guarantee of outer bounds. The manuscript must show that these inequalities are valid for the original nonlinear problem (i.e., every feasible point of the unrelaxed finite-volume system remains feasible after augmentation) and does not merely tighten the relaxation gap at the cost of excluding admissible solutions."},{"response":"We concur that quantitative assessment of tightness is necessary for readers to judge practical utility. For the 1-D and 2-D steady-state examples, where the admissible set can be enumerated or solved analytically, we will add tables reporting (i) the ratio of our interval widths to the true null-space ranges and (ii) the residual McCormick gap at the obtained optima. For the transient hexagonal-grid case we will report the same gap metric together with a comparison against a modest Monte-Carlo ensemble (10^4 realizations) to illustrate conservatism. These metrics and the associated discussion will be inserted into Section 4.","revision_made":"yes","referee_comment":"Numerical examples section: The demonstrations on the three models report no quantitative metrics of bound tightness (e.g., ratio of obtained interval widths to independently computed ground-truth ranges) or residual relaxation gap. Without such measures or comparison against a reference method, it is not possible to confirm that the remedies leave the bounds both valid and practically useful rather than excessively loose."}],"tokens_in":1476,"tokens_out":499,"duration_ms":17784,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"This paper gives a sampling-free OBBT interval method for bounding the null space in non-identifiable groundwater models, with concrete fixes for McCormick sign-coupling failures.\n\nIt discretizes Darcy's law on a finite-volume grid, relaxes the bilinear flux-gradient products with McCormick envelopes, then adds flow sign prescription and irrotationality constraints to restore physical coupling. The authors test the pipeline on a 1D steady case, a 2D steady case in four configurations, and a 2D transient hexagonal-grid case. The approach is new as a deterministic outer-bound alternative to ensemble UQ and the remedies are a practical, targeted contribution.\n\nThe central claim of guaranteed outer bounds holds on paper once the remedies are in place, and the connection to null-space theory is clean. What is missing is any reported measure of bound tightness, residual relaxation gap, or comparison against ground-truth ranges from exhaustive sampling. Without those numbers the practical value remains unclear.\n\nThe added constraints are load-bearing; if they are not shown to be valid inequalities for the original nonlinear system, the bounds lose their guarantee. The abstract does not supply that verification.\n\nThe work is for people who need conservative, deterministic UQ in subsurface flow rather than statistical coverage. It is coherent enough to deserve referee time, though the numerical validation section will need strengthening before acceptance.","headline":"The paper gives a sampling-free OBBT interval method for bounding the null space in non-identifiable groundwater models, with concrete fixes for McCormick sign-coupling failures.","tokens_in":2395,"tokens_out":356,"would_cite":false,"duration_ms":16357,"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":"Optimization-based bound tightening gives guaranteed outer bounds on uncertain variables in non-identifiable groundwater models without sampling.","keywords":["groundwater modeling","uncertainty quantification","optimization-based bound tightening","non-identifiable models","Darcy's law","finite volume discretization","McCormick relaxations","interval methods"],"falsifier":"An admissible solution found by exhaustive sampling or another method that lies outside the interval bounds computed by the OBBT procedure, or a case in which the bounds remain too loose to be useful after the constraints are applied.","tokens_in":2668,"feed_emoji":"💧","tokens_out":681,"duration_ms":25389,"temperature":0.7,"pith_summary":"The paper seeks to establish that uncertainty in groundwater models, where sparse data leave many parameter-state combinations equally plausible, can be quantified directly as tightened intervals rather than by exploring discrete realizations. It does so by casting the problem as an optimization task that extremizes each variable subject to a constraint system encoding Darcy's law, observations, and selected physical restrictions. Finite-volume discretization produces bilinear terms that are relaxed via McCormick envelopes; flow sign prescription and irrotationality constraints are added to restore physical fidelity and enable effective tightening. This approach is shown to work on 1D steady, 2D steady, and 2D transient test cases. A sympathetic reader would care because it replaces the risk of incomplete ensemble exploration with deterministic, conservative bounds that connect directly to the model's null space.","feed_headline":"Optimization tightens intervals to bound groundwater model uncertainty","feed_subtitle":"Guaranteed outer bounds replace sampling ensembles by extremizing variables over physical constraints and observations.","key_machinery":"Optimization-based Bound Tightening (OBBT) that represents uncertainty as intervals and tightens them by solving optimization problems over the discretized physical constraints and data.","core_discovery":"Optimization-based bound tightening applied to a finite-volume discretization of Darcy's law, with McCormick relaxations augmented by flow sign prescription and irrotationality constraints, produces guaranteed outer bounds on all uncertain parameters, states, and boundary conditions for non-identifiable groundwater models.","pith_inferences":["The same interval formulation could be tested on other under-determined inverse problems such as atmospheric transport or reservoir history matching.","The computed outer bounds could serve as a quantitative guide for choosing new observation locations that most reduce the null-space dimension.","Scalability on larger three-dimensional grids would depend on the tightness of the chosen relaxations and the cost of the successive optimization solves."],"forward_implications":["The method supplies conservative bounds on every uncertain variable without any sampling step.","It applies unchanged to steady-state and transient problems on Cartesian and hexagonal grids.","Flow sign prescription and irrotationality constraints each restore sign coupling and rotational consistency that plain McCormick relaxations lose.","The resulting bounds are naturally linked to the null-space structure of the under-determined inverse problem."],"fun_headline_variants":["OBBT tightens intervals to bound non-identifiable groundwater models","Guaranteed interval bounds replace sampling in groundwater models","OBBT bounds null space uncertainty for non-identifiable groundwater models","Optimization-based bounds for non-identifiable groundwater uncertainty"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The flow sign prescription and irrotationality constraints, when added to the McCormick relaxations, must keep the resulting bounds both valid and sufficiently tight for the target models.","fun_headline_variants_meta":{"raw":{"variants":["OBBT tightens intervals to bound non-identifiable groundwater models","Guaranteed interval bounds replace sampling in groundwater models","OBBT bounds null space uncertainty for non-identifiable groundwater models","Optimization-based bounds for non-identifiable groundwater uncertainty"]},"model":"grok-4.3","cost_usd":0.010301,"raw_usage":{"total_tokens":4575,"prompt_tokens":694,"num_sources_used":0,"completion_tokens":52,"cost_in_usd_ticks":103012000,"prompt_tokens_details":{"text_tokens":694,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3829,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":694,"tokens_out":52,"duration_ms":30352,"temperature":1.0,"reasoning_tokens":3829,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-27T10:57:21.178288+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"An admissible solution found by exhaustive sampling or another method that lies outside the interval bounds computed by the OBBT procedure, or a case in which the bounds remain too loose to be useful after the constraints are applied.","supporting_citations":[],"review_version":1}