{"id":"80297fa0-e715-41a3-8f22-7ad1a8556f0c","arxiv_id":"2411.13238","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Stochastic noise blurs the deterministic Busse balloon boundary of the 1D Klausmeier model, with pattern stability depending strongly on noise intensity and position inside the balloon.","lead":"This paper studies how random noise changes the stability of periodic vegetation patterns in a mathematical model of drylands, showing that noise blurs the boundary of the region where such patterns are stable. It introduces two measures, average first exit time and local wave number, to describe when and how patterns break down or reorganize.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central exit-time reversal in Fig. 8c may depend on the heuristic pulse-counting parameters used to define Texit; the blurring claim needs a robustness check.","rationale":"Reader's verdict is already CONDITIONAL, and the weakest assumption identified there is the same pulse-counting heuristic. I agree that this is the most load-bearing point: the quantitative support for 'blurring the Busse balloon' is a single reversal in average exit times, and every Texit value depends on a classification pipeline whose parameters are chosen once and never varied. The paper has independent strengths: the deterministic continuation, the new combination of local Fourier diagnostics, and the reproducible repository listed in the data availability statement. I am not claiming the reversal is wrong; it may well be robust. But the authors' own admission that pulse counting is heuristic means a parameter-sweep robustness check is necessary before the central claim can be accepted. If the ordering survives the sweep, the concern is resolved and the conditional verdict can stand; if it does not, the strongest quantitative evidence for the paper's headline conclusion would disappear. Hence no change to the reader's CONDITIONAL recommendation is needed; the required check is the natural next step.","tokens_in":14685,"tokens_out":8567,"duration_ms":100376,"concrete_test":"Recompute the a=2, sigma=0.25 panel of Figure 8 with the same 50 seeds while varying the Appendix A.3 pulse-counting parameters: Gaussian smoothing widths 32 and 128 gridpoints crossed with MinProminence values 0.15 and 0.6 (a 2x2 sweep). If the ordering E[Texit](38) > E[Texit](23) and E[Texit](38) > E[Texit](24) does not persist under every one of the four settings, the blurring reversal depends on the heuristic detection rule and the central claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest evidence for the blurring claim is the reversal in Fig. 8c: at a=2, sigma=0.25, E[Texit(k=38)] exceeds E[Texit(k=23)] and E[Texit(k=24)], although 38 lies outside the deterministic Busse balloon and 23,24 lie inside. Texit is defined as the first time the number of pulses changes (Section 2.3), and pulse number is computed by smoothing the solution with a Gaussian over 64 gridpoints and counting extrema with MinProminence 0.3 (Appendix A.3). The authors explicitly call this procedure 'in some sense heuristic' (Section 2.3). The concern is that this detector has a wave-number-dependent sensitivity: a fixed smoothing width and a fixed prominence threshold do not respond in the same way to a decaying pulse in a short-wavelength, high-k pattern (k=38, wavelength about 13 spatial units) as in a longer-wavelength pattern near the lower edge (k=23,24). If the threshold keeps registering a dying pulse as present longer for high k, the measured ordering Texit(38)>Texit(23),Texit(24) would reflect the classification rule rather than the stochastic dynamics. No test in the paper shows that the reversal survives changes to the smoothing width or MinProminence, and the GitHub scripts alone do not establish invariance.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the one-dimensional stochastic Klausmeier model (1.1) with multiplicative noise in the vegetation mortality. It defines the first exit time Texit as the first time the number of pulses changes, and it introduces local Fourier transforms to assign local wavenumbers to transient and noisy patterns. Numerical experiments show that average exit times vary strongly with rainfall a, wavenumber k, and noise intensity sigma, and that for a=2 and sigma=0.25 the mean exit time at k=38 exceeds that at k=23 and 24 even though k=38 lies outside the deterministic Busse balloon. The authors also compute long-time local-wavenumber distributions, which concentrate around the middle of the balloon, and conclude that stochasticity blurs the boundary of the deterministic Busse balloon while the boundary becomes sharp again in the limit sigma -> 0 and Tmax -> infinity.","tokens_in":14900,"tokens_out":8237,"duration_ms":83813,"significance":"The manuscript addresses a timely question: whether the deterministic Busse balloon remains a useful stability concept under stochastic forcing. The proposed tools, average first exit times and local wavenumber distributions, are natural and well motivated, and the numerical experiments are direct simulations rather than fitted outputs, which avoids circularity. If the central reversal in Figure 8c is robust, the paper gives a concrete, falsifiable statement about stochastic destabilization in a canonical ecological model, and the open-source scripts on GitHub are a strength for reproducibility. The comparison with empirical Busse balloons from Bastiaansen et al. is a useful sanity check. However, the robustness of the pulse-counting definition and the unsupported stationarity label for sigma=0.2 are load-bearing gaps that need to be addressed before the claims can be accepted.","major_comments":[{"comment":"The central reversal in Figure 8c depends on the pulse-counting detector: Texit is the first time the pulse number changes, and in Appendix A.3 the pulse number is obtained by smoothing u with a Gaussian of width 64 gridpoints and counting extrema with MinProminence 0.3, a procedure the authors themselves call heuristic in Section 2.3. Because k=38 has a much shorter wavelength than k=23 or 24, a fixed smoothing width and prominence threshold can have wavenumber-dependent detection lag: a dying short-wavelength pulse at k=38 may be registered as present for a different length of time than a dying long-wavelength pulse, purely because of the detector. The paper provides no test that the ordering E[Texit(38)] > E[Texit(23)], E[Texit(24)] survives variation of the smoothing width or MinProminence, or comparison with a different destabilization criterion. Please add a robustness study of this specific reversal, for example by recomputing Texit with smoothing widths of 32 and 128 gridpoints and MinProminence values of 0.15 and 0.6, and by checking the same ordering with a criterion based on the local wavenumber distribution or on amplitude decay.","section":"Section 2.3, Appendix A.3, Figure 8c"},{"comment":"The claim that Figures 9b and 10c show stationary distributions is not supported by the paper's own convergence test. Section 3.2 says that Appendix A.4 validates stationarity, but Appendix A.4 explicitly states that for sigma=0.2 two different initial conditions lead to different distributions at T=Tmax=10^4 and that the simulations have not yet reached the stationary distribution. Nevertheless, Figure 9b is labeled as a stationary distribution for sigma=0.2, and Figure 10c uses the same object to compare stability measures. This internal inconsistency affects the observability conclusions. Please either restrict the stationarity claim to sigma=0.25, or present the sigma=0.2 results as finite-time long-time distributions with a quantitative convergence assessment, and state whether the comparison in Figure 10 changes when only verified stationary cases are used.","section":"Section 3.2, Appendix A.4, Figures 9b and 10c"},{"comment":"Figure 8c reports E[Texit] over 50 realizations and displays the standard deviation of the exit times, but the central claim that the mean exit time at k=38 exceeds those at k=23 and 24 is not accompanied by a standard error, a confidence interval, or a pairwise significance test. With heavy-tailed exit-time data and only 50 runs, the observed reversal may be within sampling error, and the caption's note that standard deviations are added to highlight the high variability makes this concern concrete. Additionally, Texit is right-censored at Tmax, so the fraction of runs that reached Tmax for each k must be reported, since censoring biases the estimated mean differently across wavenumbers. Please provide mean uncertainties and censoring counts for the bars in Figure 8c, and test the ordering of the means.","section":"Figure 8c and Section 3.1"}],"minor_comments":[{"comment":"The caption of Figure 9 states Tmax=2500, while Section 3.2 states Tmax=10^4 for the same computations; please reconcile the two values.","section":"Figure 9 caption"},{"comment":"Figure A.4 claims a power-law relation Texit ~ sigma^alpha with alpha approximately 10 based on a visual log-log fit, and Figure 7b claims an exponential relation between a and the maximum exit time from a straight-line fit; report the number of data points, regression standard errors, and residual diagnostics for both fits, and state how the a=0.40 outlier in Figure 7b is treated.","section":"Figure A.4 and Figure 7b"},{"comment":"Appendix A.3 mentions a moving median in time to filter short fluctuations but does not specify its window length; since the pulse number and thus Texit can depend on this parameter, please give the value used.","section":"Appendix A.3"},{"comment":"The local Fourier transform window width ell=50 is fixed a priori; a brief discussion of the sensitivity of the local wavenumber statistics to ell would help readers assess the observability claims.","section":"Section 2.2"}],"recommendation":"major_revision","confidential_remarks":"The paper's main quantitative claim is plausible but rests on a detector whose parameters are not varied, and the stationarity overclaim in Section 3.2 versus Appendix A.4 is the most serious internal inconsistency. I would like the revision to include the robustness tests described in Major Comment 1 and to correct or qualify the sigma=0.2 stationarity statements. The paper is otherwise within scope for a nonlinear-science journal and the numerical study is reproducible in principle."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper asks a genuinely good question: what happens to a deterministic Busse balloon when you add noise, and it is the first to answer it for the Klausmeier model with the tools it brings in. The local wave number machinery is well motivated and does work that the global Fourier mode cannot do during transients, and the comparison of three stability notions in Fig 10 is a useful conceptual contribution. The central qualitative claim, that noise blurs the deterministic boundary, is credible and supported by the exit-time data in Fig 8. The scaling laws in Fig 7b and A.4 are plausible but fitted by eye to numerical output; the paper would be stronger with error bars or at least an honest statement that the exponents are rough. The numerics are reproducible in principle: a GitHub repository is cited and the integration scheme is standard.\n\nThe weakest point is the definition of Texit via pulse counting. The smoothing width and MinProminence threshold are fixed by hand, and a high-k pattern near k=38 has a much shorter wavelength than the k=23-24 patterns it is compared against in Fig 8c. The reversal that drives the blurring claim could in principle be an artifact of the detector: a dying pulse in a short-wavelength pattern might be registered as present for longer by a fixed-width smoother. The paper itself calls the procedure 'in some sense heuristic', but it never tests whether the reversal survives changes to the smoothing or threshold parameters. That is a real gap, and it is the main reason my verdict is conditional rather than accept. The 'stationary' label in Figs 9-10 is also an overstatement: Appendix A.4 shows that for sigma=0.2 the distribution has not actually converged, so the term means 'visually stationary on this timescale'.\n\nThe authors are honest about the limits of their study, and the self-citations to their earlier continuation work are appropriate. The paper is a serious numerical study of a real question; it deserves a referee. I would send it out with a request to add a robustness check on the pulse-counting detector and to soften the 'stationary' language. If the robustness check confirms the Fig 8c reversal, the blurring result is solid and this will be a useful reference for anyone working on stochastic pattern-forming systems.","headline":"Noise blurs the Busse balloon in the Klausmeier model, but the pulse-counting definition of exit time needs a robustness check before I'd trust the quantitative claims.","tokens_in":15510,"tokens_out":1121,"would_cite":true,"duration_ms":11054,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B36","35R60","60H15"],"pacs":[],"model":"deepseek-v4-flash","headline":"Stochastic forcing blurs the boundary of the deterministic Busse balloon in a dryland vegetation model.","keywords":["stochastic reaction-diffusion equations","Klausmeier model","Busse balloon","first exit time","local wave number","periodic patterns","dryland vegetation","multiplicative noise"],"falsifier":"Run the $a=2.0$, $\\sigma=0.25$ stochastic simulations with the pulse-counting thresholds varied (smoothing width and prominence) and check whether wave number 38 still has a longer average first exit time than wave numbers 23 and 24; if the ordering reverses, the claimed blurring is an artifact of the pulse-counting definition.","tokens_in":14347,"feed_emoji":"🌾","tokens_out":6437,"duration_ms":61609,"temperature":0.7,"pith_summary":"This paper asks what happens to the Busse balloon—the set of wave numbers that are stable periodic patterns in a deterministic reaction-diffusion model—when the model is forced by noise. Using the one-dimensional Klausmeier model for dryland vegetation, it claims that stochastic forcing blurs the balloon's boundary: a wave number outside the deterministic balloon can have a longer average first exit time than wave numbers inside it. The authors introduce two numerical tools, the average first exit time and the local wave number, and show that stability under noise depends strongly on the rainfall parameter, the noise intensity, and the wave number's position inside the balloon. If true, this means deterministic Busse balloons are only a sharp predictor of observable patterns in the zero-noise limit; with realistic noise, pattern persistence is graded rather than binary.","feed_headline":"Noise blurs the Busse balloon's edge in a dryland vegetation model","feed_subtitle":"Deterministically unstable wave numbers can outlive stable ones when mortality noise is added to the Klausmeier model.","key_machinery":"Busse balloon: the region in parameter space of deterministically stable periodic wave numbers, traced by numerical continuation. First exit time: the first time the pulse number changes, estimated by averaging over many noise realizations; it mimics deterministic stability in the stochastic setting. Local wave number: the predominant mode of a Gaussian-windowed Fourier transform with width $\\ell=50$; it captures both the number of pulses and their spatial distribution. The pulse number itself is counted by smoothing the solution and locating extrema with a prominence threshold, a heuristic that defines when destabilization is said to occur.","core_discovery":"The paper studies the one-dimensional stochastic Klausmeier model (1.1), a reaction-diffusion system for dryland vegetation with multiplicative noise in the mortality rate. It claims that stochastic forcing does not simply destroy the deterministic Busse balloon but blurs its boundary: for fixed rainfall $a=2.0$ and noise $\\sigma=0.25$, the average first exit time of the deterministically unstable wave number 38 exceeds that of the deterministically stable wave numbers 23 and 24, so the observed boundary of the balloon is no longer sharp. The blurring is controlled by position within the balloon and noise level: patterns in the centre survive longest, exit times grow exponentially in the rainfall parameter $a$ and obey a power law in $\\sigma$, and the boundary sharpens again in the limit $\\sigma\\to 0$ and $T_{\\max}\\to\\infty$. The paper also introduces local wave numbers to describe what solutions look like after destabilization, showing that typical stochastic solutions are patchworks of different local wave numbers rather than a single global mode.","pith_inferences":["If the blurring mechanism is generic, the same pulse-deletion route should blur Busse balloons in other stochastic reaction-diffusion systems, such as Swift-Hohenberg or Gray-Scott; a direct test would be to repeat the exit-time scan for those models.","The measured exponential-in-$a$ and power-law-in-$\\sigma$ relations for exit times could be turned into a scaling law for pattern persistence in dryland vegetation, once dimensional rainfall and disturbance intensity are calibrated.","A Markov chain on pulse numbers with transition rates estimated from short stochastic runs could convert the blurred balloon into a quantitative prediction of how fast patterns degrade, although the paper notes the state definition is delicate.","In two-dimensional stripe patterns, the noise is correlated along the stripe direction in the present one-dimensional setup; in true two dimensions, local defects may nucleate differently, so the blurring could appear at lower noise levels."],"forward_implications":["Within the stochastic model, patterns near the centre of the deterministic Busse balloon are the most persistent, but the deterministic boundary loses its predictive power for observability once noise is present.","Reducing the noise intensity increases exit times most strongly in the balloon's interior, so the sharp deterministic boundary is recovered only in the double limit $\\sigma\\to 0$ and $T_{\\max}\\to\\infty$.","Typical stochastic solutions are not close to a single periodic pattern: the stationary distribution of local wave numbers is spread across the balloon and shifts with $a$ and $\\sigma$.","Because exit times are exponentially sensitive to rainfall and power-law sensitive to noise, the relevant pattern persistence timescale can change by orders of magnitude within the parameter range studied.","The three measures of stability peak at different wave numbers (31 for exit time, 30 for deterministic selection, 28 for stationary distribution), so 'most stable' depends on which observable is used."],"supporting_citations":[{"why":"Introduces the deterministic Klausmeier model whose stochastic version is studied here.","marker":"[20]"},{"why":"Supplies the continuation methods used to compute the stable and unstable periodic patterns used as initial conditions.","marker":"[33, 34]"},{"why":"Supplies the continuation software used to trace the deterministic Busse balloon boundary in Figure 1.","marker":"[10]"},{"why":"Provides the numerical integration scheme and the algorithm for generating spatially coloured noise used in the simulations.","marker":"[22]"},{"why":"Introduces the local Fourier transform and local wave number used to classify transient and noisy patterns.","marker":"[38]"},{"why":"Establishes the average first exit time as a measure of ecological resilience, motivating its use as a stochastic stability measure.","marker":"[1]"},{"why":"Provides empirical dryland pattern data whose wave-number distributions are compared with the model's stochastic Busse balloon.","marker":"[6]"}],"fun_headline_variants":["Noise blurs desert vegetation pattern edges","Stochastic forcing smears the Busse balloon's rim","Unstable waves outlive stable ones under noise","Rainfall noise erases sharp pattern stability zones","Klausmeier noise fuzzes the Busse balloon boundary"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The definition of first exit time rests on a heuristic pulse-counting procedure (Gaussian smoothing width 64 gridpoints, prominence threshold 0.3); if a different threshold is used, measured exit times and the reported blurring of the balloon boundary could change.","fun_headline_variants_meta":{"raw":{"variants":["Noise blurs desert vegetation pattern edges","Stochastic forcing smears the Busse balloon's rim","Unstable waves outlive stable ones under noise","Rainfall noise erases sharp pattern stability zones","Klausmeier noise fuzzes the Busse balloon boundary"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000624,"raw_usage":{"total_tokens":2842,"prompt_tokens":852,"completion_tokens":1990,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":468,"completion_tokens_details":{"reasoning_tokens":1925}},"tokens_in":468,"tokens_out":1990,"duration_ms":16690,"temperature":1.0,"reasoning_tokens":1925,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:39:39.717753+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the $a=2.0$, $\\sigma=0.25$ stochastic simulations with the pulse-counting thresholds varied (smoothing width and prominence) and check whether wave number 38 still has a longer average first exit time than wave numbers 23 and 24; if the ordering reverses, the claimed blurring is an artifact of the pulse-counting definition.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the deterministic Klausmeier model whose stochastic version is studied here."},{"cited_title":"Doedel, A","cited_arxiv_id":null,"evidence_quote":"Supplies the continuation software used to trace the deterministic Busse balloon boundary in Figure 1."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the numerical integration scheme and the algorithm for generating spatially coloured noise used in the simulations."},{"cited_title":"Vi˜ nals, E","cited_arxiv_id":null,"evidence_quote":"Introduces the local Fourier transform and local wave number used to classify transient and noisy patterns."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the average first exit time as a measure of ecological resilience, motivating its use as a stochastic stability measure."},{"cited_title":"Bastiaansen, O","cited_arxiv_id":null,"evidence_quote":"Provides empirical dryland pattern data whose wave-number distributions are compared with the model's stochastic Busse balloon."}],"review_version":1}