{"id":"6979e43d-a575-4e3d-83fd-ee838f75235d","arxiv_id":"1908.05507","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":6,"one_line_summary":"In the Rosenzweig-MacArthur predator-prey model with a slowly declining carrying capacity, a rate-induced critical transition causes temporary collapse of both populations, with a canard trajectory as the tipping threshold.","lead":"Slow environmental decline can still cause a sudden collapse in a classical predator-prey model, provided the decline is faster than the ecosystem's balance point can move. The collapse occurs even though the system never loses stability, and the study maps exactly which starting conditions crash and which survive.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The canard-threshold analysis is only verified at kappa=0.01, the lower edge of the cited ecological range kappa in [0.01,0.08]; without a robustness check the quantitative collapse predictions may not generalize.","rationale":"The paper's novel contribution is not simply that ramping can tip the Rosenzweig-MacArthur model (Siteur et al. already demonstrated this for a time-varying growth rate), but that the tipping threshold is a folded-saddle canard and that this mechanism explains why collapse can occur for environmental rates slower than the slowest internal timescale. That mechanism is inherently a singular-limit (kappa -> 0) construction. The text invokes geometric singular perturbation theory to argue that for 0 < kappa << 1 the stable and unstable slow manifolds intersect along a maximal canard near the singular canard, but no theorem or numerical study quantifies the radius of validity in kappa. All simulations in the paper use kappa = 0.01, the smallest and most favorable value in the ecological range cited. At kappa = 0.08, the fast-slow separation is only about a factor of 12, and the O(kappa) splitting of the slow manifolds could shift the canard threshold by an amount comparable to the collapse-prone region; the critical-rate curves in Figs. D.1B and F.1 are singular-limit predictions. Without a check, the central quantitative conclusion is not established for the upper part of the claimed parameter range. The other two weaknesses are real but less central: the sign inconsistency in the desingularization rescaling is contradicted by the correct equations and does not affect the equilibrium locations, and initial-condition dependence is explicitly analyzed in Section 4. Thus the kappa-robustness is the load-bearing concern. This supports the reader's CONDITIONAL verdict; no verdict change is needed.","tokens_in":23986,"tokens_out":21011,"duration_ms":191984,"concrete_test":"Run the ramped system (9)-(11) with kappa = 0.04, 0.06, and 0.08 (keeping eta = 0.8, phi0 = 0.1, epsilon = 1e-6, and the same initial states as in Fig. 3B and Fig. 6), and compute the critical rate separating tracking from collapse. Compare this full-system critical rate with the singular-canard prediction obtained from Eqs. (E.17)-(E.18) and Fig. F.1. If the full-system threshold deviates from the singular-canard prediction by more than the expected O(kappa) amount (e.g., more than 10-20% at kappa = 0.08), or if the sharp threshold structure disappears, then the paper should either restrict its ecological claims to kappa < 0.01 or provide a finite-kappa correction to the threshold.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central construction is the folded-saddle canard threshold, derived in the singular limit kappa -> 0 (Appendix E) and then used to predict critical rates and collapse-prone initial conditions (Figs. 5-7, D.1, F.1). Every numerical verification in the paper fixes kappa = 0.01, which is the lower bound of the ecological range kappa in [0.01,0.08] cited in Section 3. In fast-slow systems with two slow variables, the maximal canard is indeed generic for sufficiently small kappa, but the distance between the singular canard and the full maximal canard is O(kappa) (with additional singular effects near the fold). At kappa = 0.08 the time-scale separation is only about a factor of 12, so the threshold can shift by an amount comparable to the width of the collapse-prone region, and the normal-hyperbolicity assumptions behind the singular perturbation picture may begin to fail. Because the headline conclusion ('very slow is too fast') is meant to apply to natural zooplankton-phytoplankton systems across this whole kappa range, the absence of any test at kappa > 0.01 leaves the quantitative predictions -- critical rates and the basin of collapse-prone states -- unvalidated in exactly the parameter range where ecological relevance is claimed. A secondary, less central issue is the sign inconsistency in the desingularization rescaling in Appendix E; however, the desingularized equations themselves appear to have the correct sign, so this is likely a typographical error rather than the cause of the numerical results.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the three-dimensional slow-fast Rosenzweig-MacArthur model (Eqs. 9–11) in which the carrying capacity declines linearly in time (φ increases at rate r). For an initial condition on the coexistence equilibrium e3(φ0) close to the fold, the authors show numerically that there is a critical ramp rate between r = 0.005 and r = 0.006 at which the system switches from tracking the moving equilibrium to a rate-induced critical transition in which the prey density temporarily collapses. They explain the mechanism using geometric singular perturbation theory: the reduced dynamics on the folded critical manifold has a folded saddle singularity, and the singular canard emanating from it separates tracking from tipping. They derive the folded saddle coordinates analytically (Appendix E), compute the singular canard as the tipping threshold (Figs. 5–7), characterize the set of collapse-prone initial conditions, and study the dependence of the critical rate on the predator's handling time η (Appendices D and F). Numerical simulations of the original three-dimensional system confirm the predicted threshold at κ = 0.01.","tokens_in":24373,"tokens_out":11893,"duration_ms":112783,"significance":"The result is significant because it transfers the concept of rate-induced critical transitions to a canonical ecological model with an ecologically motivated forcing (declining resource/carrying capacity), and it produces an analytic, parameter-free tipping threshold (the singular canard) rather than a fitted quantity. The claim that collapse can occur for ramp rates slower than the slowest intrinsic timescale is counterintuitive and is well explained through the fast-moving equilibrium argument. The paper also gives falsifiable predictions: the critical rate as a function of η and the location of collapse-prone initial states. Strengths include the explicit derivation of the desingularized system, the folded-saddle analysis, and the direct numerical verification in the original system. The main weakness is the lack of validation across the ecological timescale-separation range cited in the paper.","major_comments":[{"comment":"All numerical simulations fix κ = 0.01, which is the lower bound of the ecological range κ ∈ [0.01, 0.08] cited from Edwards and Brindley (1999). The singular canard is a κ → 0 construction, and the distance between the singular and maximal canards is O(κ), with additional effects near the fold; at κ = 0.08 the time-scale separation is only about a factor of 12 and the normal-hyperbolicity assumptions weaken. Since the paper's quantitative predictions (critical rates, collapse-prone initial-state region) are used to draw conclusions for natural plankton systems across this whole range, the authors should simulate the original system (9)–(11) for at least two larger κ values (e.g., κ = 0.04 and κ = 0.08) and report how the tipping threshold and critical rates deviate from the singular-canard prediction. This is a load-bearing missing check for the ecological claim.","section":"Section 3 and Figs. 3–7, D.1, F.1"},{"comment":"The paper identifies the maximal canard as the tipping threshold but computes only the singular canard (κ = 0), and the verification in Fig. 6 is visual rather than quantitative. Because the threshold underpins the collapse-prone sets in Figs. 5–7 and the critical rates in Fig. F.1, I ask for a quantitative check at κ = 0.01, for example by bisection in u0 for fixed v0 and φ0 to locate the actual threshold in the original system and comparing it with the singular canard curve. This would also provide a baseline for the κ-robustness test requested above.","section":"Section 4 and Fig. 6"}],"minor_comments":[{"comment":"The stated time rescaling t = −2φη(u_F − u)s has the wrong sign; with that choice one would obtain du/ds = −Λ, whereas Eqs. (E.17)–(E.18) follow from t = +2φη(u_F − u)s. The desingularized equations appear correct, but the text should be fixed.","section":"Appendix E, text after Eq. (E.16)"},{"comment":"The phrase 'collapse of a predator-prey system' is stronger than the demonstrated dynamics, in which the prey density temporarily collapses and then the system recovers to the coexistence equilibrium (Section 3, Fig. 1C). The authors correctly state the temporary nature in the body and discuss the role of noise in converting it to extinction, but the title and abstract should be tempered (e.g., 'temporary collapse') to avoid overstating the deterministic result.","section":"Title and Abstract"},{"comment":"The text near line 102 reads 'whitin a bounded φ-interval'; this should be 'within a bounded φ-interval'. Also, the citation '(S. Sakar and P.S. Dutta, unpublished manuscript)' is not listed in the reference list and should either be removed or given a full reference.","section":"Section 1"},{"comment":"There are two typos: 'the back dashed line' should be 'the black dashed line', and the text uses 'black solid line' for two different curves in the same paragraph; please re-read the caption and surrounding text of Fig. F.1 for consistency.","section":"Appendix F and Fig. F.1"},{"comment":"The caption states u0 = (1−η)^−1 and v0 = (1−φ0u0)(1 + ηu0); it may help the reader to note explicitly that this initial condition is exactly the coexistence equilibrium e3(φ0) on the critical manifold.","section":"Fig. 3B caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good fit for the journal and the core mechanism is sound. The main substantive issue is the lack of robustness tests over the cited κ range, which I believe the authors can address with additional numerical experiments. The sign typo in Appendix E and the overstatement in the title/abstract should be fixed. I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, well-explained application of the rate-induced tipping framework to the Rosenzweig-MacArthur model with a linearly declining carrying capacity. The central claim holds: the system can collapse even when the environmental ramp is slower than the slowest internal timescale, and the tipping threshold is the singular canard associated with a folded saddle. The reduced-system derivation in Appendix E checks out, and the numerical integrations of the full three-dimensional system back the singular-limit analysis.\n\nWhat is genuinely new: Siteur et al. (2016) ramped the prey growth rate and used a steady-lag approximation; this paper ramps the carrying capacity - a more ecologically accessible parameter - and gives the full geometric singular perturbation treatment: folded saddle, singular canard as threshold, and critical rates as a function of initial state and handling time. Appendix F also verifies the eigenvector approximation of the critical rate against direct simulation, so that part is solid.\n\nWhere the soft spots are. The quantitative predictions - critical rates and the collapse-prone basin - are verified only at kappa = 0.01, the lower edge of the ecological range kappa in [0.01, 0.08] quoted from Edwards and Brindley. The full canard lies O(kappa) from the singular canard, and at kappa = 0.08 the time-scale separation is only a factor of twelve. I don't think the mechanism disappears there, but the paper should demonstrate that the threshold and critical rates do not shift enough to alter the ecological conclusion. That is a real gap, not a nitpick. A smaller issue: Appendix E writes the desingularizing rescaling as t = -2 phi eta (u_F - u) s and then gives du/ds = Lambda; with that rescaling the equation should have -Lambda. The desingularized system itself appears to have the correct sign, so this is a typo, not a substantive error. Finally, the headline result applies to initial states near the fold; the paper is transparent about this, but the abstract could be read as claiming more.\n\nThis paper deserves a serious referee. It is careful, the math is honest, and the ecological message is worth taking seriously. With the kappa robustness check and the sign typo fixed, it is publishable in a good applied-dynamical-systems or ecology journal. I would cite it.","headline":"A careful, convincing application of canard theory to R-tipping in a classic predator-prey model; the main gap is that all numerics sit at the lower edge of the cited time-scale separation range.","tokens_in":24868,"tokens_out":4736,"would_cite":true,"duration_ms":39877,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34E15","37N25","92D25"],"pacs":[],"model":"deepseek-v4-flash","headline":"A smooth, slow decline in habitat quality can still tip a predator-prey system into sudden collapse, and the tipping threshold is a canard trajectory, not a bifurcation.","keywords":["rate-induced critical transition","R-tipping","population collapse","Rosenzweig-MacArthur model","canard trajectory","folded saddle singularity","carrying capacity decline","fast-slow dynamical systems"],"falsifier":"Numerically integrate the ramped system (9)-(11) at $\\kappa = 0.08$, the upper end of the ecological range the paper cites, with $\\eta=0.8$ and the initial state on the moving equilibrium $e_3(\\varphi_0)$, for ramping rates on both sides of the predicted canard threshold. If the predicted collapse does not occur for $r$ above the threshold, or occurs for $r$ below it, the canard prediction does not carry over to weaker timescale separation.","tokens_in":23825,"feed_emoji":"🌿","tokens_out":10245,"duration_ms":94466,"temperature":0.7,"pith_summary":"The paper studies a well-known predator-prey model in which the prey's habitat quality declines smoothly while predators reproduce much more slowly than prey. It claims that this system can undergo a rate-induced critical transition—a sudden, temporary collapse of the prey population—even when the environmental change is slower than the slowest intrinsic ecosystem timescale. The collapse is not preceded by a bifurcation or by loss of linear stability; a unique stable coexistence equilibrium exists for every fixed habitat quality, but the system can fail to track it as the equilibrium moves. The threshold separating tracking from collapse is a maximal canard trajectory associated with a folded saddle singularity. The practical stakes: ecosystems may be tipped by how fast conditions change, not only by how bad they become.","feed_headline":"Even very slow habitat decline can collapse a predator-prey system","feed_subtitle":"A canard trajectory, not a loss of stability, marks the tipping threshold in the classic model.","key_machinery":"The key object is the folded saddle singularity on the fold curve of the critical manifold, together with the singular canard trajectory through it. Desingularization removes the zero denominator of the reduced slow flow at the fold, turning the folded singularity into a regular equilibrium of a rescaled system; the stable manifold of that equilibrium is the tipping threshold. The speed of the moving equilibrium, $|\\dot e_3| = r/(1-\\eta)^2$, carries the explanation of why a slow ramp can be too fast: for handling time $\\eta$ close to one, even small $r$ makes the equilibrium sweep through phase space faster than the system can follow.","core_discovery":"The central discovery is that the ramped Rosenzweig-MacArthur system—fast prey, slower predator, and a third equation that linearly increases $\\phi$ (proportional to the inverse carrying capacity)—has a critical rate of environmental change below the slowest ecosystem timescale. When the ramping rate exceeds that critical rate, an initial state on the stable part of the folded critical manifold is drawn to the fold and jumps off in the fast direction, sending prey density to extremely low values while predator density remains high; the resulting overconsumption produces a temporary collapse of the prey population. This happens in a parameter range in which the coexistence equilibrium is linearly stable for every fixed $\\phi$, so standard stability analysis would see no danger. The paper computes the boundary between tracking and tipping as the singular/maximal canard trajectory through the folded saddle singularity; it is generic because two slow variables (predator density and $\\phi$) are present. The same threshold can be converted into a critical ramping rate and into a description of all collapse-prone initial states.","pith_inferences":["If the folded-saddle/canard mechanism is generic for systems with two slow variables, the same rate-induced collapse should appear in other population and ecosystem models whose critical manifold has a fold and whose environment is ramped; this could be tested in model families before field confirmation exists.","The mathematical critical rate (where the trajectory crosses the fold) and the ecological critical rate (where prey density drops below a conservation threshold) can differ substantially, so management targets should state which threshold they are using.","Adding noise to the ramped system should make the extinction probability depend on how long the prey spends at low density during the collapse; measuring that duration would quantify the extinction risk.","Because tipping here occurs without any loss of linear stability, variance-based early-warning signals would not precede it; monitoring the rate of environmental change itself may be a more practical warning indicator."],"forward_implications":["An environmental change that never leaves the range where a stable coexistence equilibrium exists can still collapse the prey population, provided the rate of change crosses the canard threshold.","The critical rate is not fixed by the ecosystem alone: it is smaller for longer predator handling times (inefficient predators) and larger for shorter handling times, following a highly nonlinear curve.","For a fixed ramping rate, the paper gives the full set of initial predator and prey densities that are collapse-prone; increasing the rate expands that set and can turn a tracking initial state into a tipping one.","The collapse is temporary in the deterministic model, but it leaves the prey at very low density, so with demographic or environmental noise the likely endpoint is extinction and breakdown of the system.","Observed regime shifts may be rate-induced rather than bifurcation-induced, so classifying a transition by stability analysis alone can miss the actual mechanism."],"supporting_citations":[{"why":"Develops the rate-induced tipping framework for slow-fast systems and shows that the singular canard from a folded saddle separates tracking from tipping; this paper applies that method.","marker":"Wieczorek et al. (2011)"},{"why":"Defines R-tipping and contrasts rate-induced tipping with bifurcation- and noise-induced tipping, giving the terminology the paper uses.","marker":"Ashwin et al. (2012)"},{"why":"The earlier demonstration of a rate-induced critical transition in a time-scaled Rosenzweig-MacArthur model with changing growth rate; this paper extends the analysis to a decline in carrying capacity.","marker":"Siteur et al. (2016)"},{"why":"Establishes the generic existence of canards and folded saddle singularities in systems with two slow variables, which underlies the threshold computation.","marker":"Szmolyan and Wechselberger (2001)"},{"why":"Provides the geometric singular perturbation theory used to define the critical manifold and the perturbed slow manifolds.","marker":"Fenichel (1979)"},{"why":"Source of the predator-prey model whose ramped version is studied in the paper.","marker":"Rosenzweig and MacArthur (1963)"},{"why":"Supplies the plankton parameter range $\\kappa \\in [0.01,0.08]$ used to argue that the computed critical rates are ecologically plausible.","marker":"Edwards and Brindley (1999)"},{"why":"Shows how folded singularities give rise to tipping thresholds in multi-scale systems and supports the claim that canard thresholds are not exceptional.","marker":"Perryman and Wieczorek (2014)"}],"fun_headline_variants":["Rate-induced collapse: slow environmental change can crash predator-prey systems","Even slow habitat decline can trigger a sudden predator-prey collapse","Rate, not level: slow habitat decline tips predator-prey to collapse","Slow habitat decline can still collapse predator-prey, no stability loss"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on a strong timescale separation, $\\kappa = 0.01$ in all simulations, because the canard threshold is an asymptotic object whose accuracy at finite $\\kappa$ is not tested across the cited ecological range $\\kappa \\in [0.01,0.08]$.","fun_headline_variants_meta":{"raw":{"variants":["Rate-induced collapse: slow environmental change can crash predator-prey systems","Even slow habitat decline can trigger a sudden predator-prey collapse","Rate, not level: slow habitat decline tips predator-prey to collapse","Slow habitat decline can still collapse predator-prey, no stability loss"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001135,"raw_usage":{"total_tokens":4758,"prompt_tokens":1029,"completion_tokens":3729,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":645,"completion_tokens_details":{"reasoning_tokens":3654}},"tokens_in":645,"tokens_out":3729,"duration_ms":26459,"temperature":1.0,"reasoning_tokens":3654,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:12:25.673889+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the ramped system (9)-(11) at $\\kappa = 0.08$, the upper end of the ecological range the paper cites, with $\\eta=0.8$ and the initial state on the moving equilibrium $e_3(\\varphi_0)$, for ramping rates on both sides of the predicted canard threshold. If the predicted collapse does not occur for $r$ above the threshold, or occurs for $r$ below it, the canard prediction does not carry over to weaker timescale separation.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Develops the rate-induced tipping framework for slow-fast systems and shows that the singular canard from a folded saddle separates tracking from tipping; this paper applies that method."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines R-tipping and contrasts rate-induced tipping with bifurcation- and noise-induced tipping, giving the terminology the paper uses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The earlier demonstration of a rate-induced critical transition in a time-scaled Rosenzweig-MacArthur model with changing growth rate; this paper extends the analysis to a decline in carrying capacity."},{"cited_title":"and Wechselberger, M","cited_arxiv_id":null,"evidence_quote":"Establishes the generic existence of canards and folded saddle singularities in systems with two slow variables, which underlies the threshold computation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the geometric singular perturbation theory used to define the critical manifold and the perturbed slow manifolds."},{"cited_title":"and MacArthur, R","cited_arxiv_id":null,"evidence_quote":"Source of the predator-prey model whose ramped version is studied in the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the plankton parameter range $\\kappa \\in [0.01,0.08]$ used to argue that the computed critical rates are ecologically plausible."},{"cited_title":"and Wieczorek, S","cited_arxiv_id":null,"evidence_quote":"Shows how folded singularities give rise to tipping thresholds in multi-scale systems and supports the claim that canard thresholds are not exceptional."}],"review_version":1}