pith. sign in

arxiv: 2604.16291 · v1 · submitted 2026-04-17 · 🧮 math.DS

Global dynamics and regime shifts in a resource-consumer model with facilitation and habitat loss

Teodoro Mayayo , Josep Sardany\'es , Joan Torregrosa This is my paper

Pith reviewed 2026-05-10 07:06 UTC · model grok-4.3

classification 🧮 math.DS
keywords resource-consumer modelfacilitationhabitat lossheteroclinic bifurcationlimit cyclePoincaré diskpiecewise-linear approximationstochastic dynamics
0
0 comments X

The pith

A resource-consumer model with facilitation has a unique stable limit cycle whose coexistence region shrinks via an analytically characterized heteroclinic bifurcation as habitat loss increases.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies a cubic planar dynamical system for resource and consumer populations in which positive plant-plant interactions are progressively weakened by habitat destruction. All phase portraits are enumerated in the Poincaré disk under biologically plausible parameter ranges, revealing a single stable limit cycle in the positive quadrant. The authors derive the exact curve of the heteroclinic bifurcation that links saddle points at infinity and drives simultaneous extinction of both populations, thereby showing how the oscillatory coexistence domain contracts with rising habitat loss. A piecewise-linear approximation supplies closed-form expressions for the same bifurcation threshold, and stochastic perturbations to the resource are shown to move the system across this threshold at lower habitat-loss values than in the deterministic case.

Core claim

The cubic system possesses a unique stable limit cycle. The heteroclinic bifurcation curve, expressed as a function of the habitat-loss parameter, is obtained analytically and bounds the region of parameter space in which interior equilibria and cycles exist; beyond this curve both populations collapse to zero. The piecewise-linear reduction preserves the qualitative dynamics and yields explicit formulas confirming the location of the bifurcation. Stochastic forcing applied to the resource accelerates crossing into the co-extinction regime.

What carries the argument

The analytically derived heteroclinic bifurcation curve in the Poincaré compactification of the cubic system, which separates the basin of the stable limit cycle from the extinction states under increasing habitat loss.

If this is right

  • The interval of parameters supporting stable oscillations narrows monotonically as the habitat-loss intensity rises.
  • The piecewise-linear model supplies explicit algebraic expressions for the critical habitat-loss value at which the limit cycle disappears.
  • Extrinsic noise on the resource population shifts the effective bifurcation point to smaller habitat-loss values, causing earlier co-extinctions.
  • Global phase-portrait enumeration shows that the only attractors in the ecological domain are the interior limit cycle and the extinction equilibria.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Preserving facilitation networks in real ecosystems could widen the buffer of oscillatory persistence against moderate habitat reduction.
  • Time-series observations of declining cycle amplitude followed by abrupt joint extinction would provide a testable signature of the predicted regime shift.
  • The same cubic structure and bifurcation-tracking approach may apply directly to other facilitative or mutualistic consumer-resource systems.
  • Incorporating spatial heterogeneity would test whether the heteroclinic mechanism survives in spatially extended versions of the model.

Load-bearing premise

The chosen cubic polynomial terms for facilitation and the linear scaling for habitat loss accurately capture the dominant ecological interactions without introducing spurious dynamical behaviors.

What would settle it

Numerical continuation or direct integration that reveals either multiple coexisting limit cycles or the lack of a heteroclinic orbit at the parameter values predicted by the analytic curve would falsify the uniqueness and bifurcation results.

Figures

Figures reproduced from arXiv: 2604.16291 by Joan Torregrosa, Josep Sardany\'es, Teodoro Mayayo.

Figure 1
Figure 1. Figure 1: x0 x0+x1 2 x1 x0+x1 (b+1)x0x1 [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Transcritical bifurcation diagrams around x0 = xe (left) and xe = x1 (right). In black we represent the equilibria of saddle type, while in gray we represent the node type. 3.7. Slow-Fast phenomena. For 0 < ε ≪ 1, systems of the form x˙ = f(x, y, ε), y˙ = ε g(x, y, ε), (8) are said to exhibit a slow-fast structure. The variable x evolves on a fast time scale, with dynamics of order O(1), whereas y evolves … view at source ↗
Figure 3
Figure 3. Figure 3: Parameter space in (xe, F) showing the different ecological regimes. The grey region ΩNLC corresponds to the no-limit-cycle regime (see Corollary 3.26); it includes both the collapse region, where the resource and the consumer unavoidably go extinct, and the static coexistence regime, where (xe, ye) is a stable equilibrium. The central white region represents the potential persistent oscillation regime. Po… view at source ↗
Figure 4
Figure 4. Figure 4: Parameter space (xe, F). The orange curve of the form xe(F) corresponds to the heteroclinic bifurcation curve. The green curve corresponds to the Hopf bifurcation. (a) Ω1. (b) Ω2. (c) Ω3. (d) Ω4 [PITH_FULL_IMAGE:figures/full_fig_p023_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Phase portraits of all the regions in the parameter space (xe, F), where F > 0 and x0 < xe < x1. Stable and unstable manifolds are colored in blue and red, respectively. 3.10. Piecewise linear model. In this section we introduce a piecewise linear (PWL) approximation of the original model. The motivation for this approach is twofold. First, the PWL system allows us to write the solutions of both vector fie… view at source ↗
Figure 6
Figure 6. Figure 6: Phase portraits of all the regions in the parameter space (xe, F), where F > 0 and x0 < xe < x1. Stable manifolds are coloured blue and unstable manifolds coloured in red. x0 √ x0x1 x0+x1 2 x1 Ω1 Ω4 Ω2 Ω6 Ω7 [PITH_FULL_IMAGE:figures/full_fig_p027_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Parameter space (xe, F). The orange curve of the form xe(F) corresponds to the heteroclinic bifurcation curve (region Ω5). The green curve corresponds to the pseudo-Hopf bifurcation (region Ω3). The grey curves separate the no return regions in the parameter space from the ones that have return. They are analogous to the focus-node transition curves separating the node and focus. 3.10.3. Effect of habitat … view at source ↗
Figure 8
Figure 8. Figure 8: Survival probability fixing a certain parameter value F over xe when increas￾ing the noise intensity σ. The transition of survival probability from one to zero for σ = 0 corresponds to the heteroclinic bifurcation. So for σ > 0, there is that the bifurcation anticipates. The stability region between the Hopf and the heteroclinic bifurcation corre￾sponds to the vertical length of xe, given an F and a certai… view at source ↗
Figure 10
Figure 10. Figure 10: Extinction times for fixed parameters x0 = 1, x1 = 3, and F = 1, with xe ∈ [x1,(x0+x1)/2] shown in linear-log scale. Each curve corresponds to a different noise intensity: green for σ = 5, yellow for σ = 2, blue for σ = 1, and black for σ = 0. Solid lines represent the mean extinction time, and the shaded areas represent the standard deviation. Each time series ends at the heteroclinic bifurcation point, … view at source ↗
Figure 9
Figure 9. Figure 9: Time series showing the effect of noise on extinction and persistence. Fixed F = 1; we vary both xe and σ. The choice of xe distinguishes between values near the Hopf point and values near the heteroclinic bifurcation. Black denotes x(t), and grey denotes y(t). For values of xe near the Hopf bifurcation and under high noise intensity, the limit cycle is both small and surrounded by a narrow basin of attrac… view at source ↗
Figure 11
Figure 11. Figure 11: Local charts (Uk, ϕk), k = 1, 2, 3, of la es￾fera de Poincar´e. Here, d = max{deg(P), deg(Q)}. Each chart corresponds to a local parametrization of the vector field in the sphere, allowing us to describe the dynamics near the equator (i.e., at infinity in the plane). Local charts (U1, ϕ1) and (U2, ϕ2) study the points at infinity. Definition A.1 (Finite singular points). Singular points lying in S 2 \ S 1… view at source ↗
Figure 12
Figure 12. Figure 12: Return map near the heteroclinic connection. The limit cycle occurs when r1(ρ) = r2(ρ). D.1.1. Map S1. From equation (14), we obtain (A1,1ρ − A1,0)B1 C1 = r1D1(E0,1 − E1,1r1) υ F1 , or, equivalently, (A1,1ρ − A1,0)z1 = − r1D1(E0,1 − E1,1r1) υ , (16) where the constants are chosen to be positive. The exponent υ is given by υ := − (x1 − x0)(x0x1 − x 2 e ) (x0 + x1 − 2xe)(x1 − xe)x0 < 0, and we introduce the… view at source ↗
read the original abstract

Modelling how populations respond to habitat loss is crucial for understanding ecosystem stability, especially when positive interactions among resource species, such as plant-plant facilitation, play a key role. Habitat loss not only reduces available organic nutrients and space for primary producers but also disrupts the positive feedbacks that sustain resource populations, thereby affecting consumer persistence and the overall system's stability. We analyse a cubic planar model describing resource-consumer dynamics with facilitation under progressive habitat loss. Our study characterizes the parameter space and enumerates all the phase portraits within the Poincar\'e disk under ecologically relevant conditions. We show that the system has a unique stable limit cycle and characterize analytically the heteroclinic bifurcation curve involving the collapse of the resource and the consumer, enabling us to determine how the parameter region sustaining coexistence oscillations narrows under habitat destruction. To further explore these dynamics, we construct a piecewise-linear (PWL) approximation that preserves the system's qualitative behaviour, allowing us to obtain an explicit expression for the heteroclinic bifurcation. Finally, we investigate how extrinsic noise affecting the resource species impacts the overall dynamics, showing that stochasticity can anticipate the onset of the heteroclinic bifurcation causing earlier co-extinctions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript studies the global dynamics of a cubic planar resource-consumer model incorporating plant facilitation and habitat loss. It enumerates all possible phase portraits in the Poincaré disk, proves the existence of a unique stable limit cycle, and provides an analytical description of the heteroclinic bifurcation curve that governs the collapse of both species. The analysis shows how habitat destruction reduces the parameter region for oscillatory coexistence. A piecewise linear approximation is developed to obtain an explicit expression for the bifurcation, and the influence of extrinsic noise on the resource is explored, indicating that noise can cause the system to reach the extinction bifurcation earlier.

Significance. Should the claims be verified, this paper makes a valuable contribution to mathematical ecology by delivering a complete qualitative analysis of an important model class. The enumeration of phase portraits and the analytical bifurcation curve are notable strengths, providing a parameter-free understanding of regime shifts. The PWL approximation and stochastic extension add practical value by enabling explicit calculations and highlighting the role of noise in accelerating extinctions. These results could inform models of ecosystem resilience under environmental change.

major comments (2)
  1. [§3] §3: The analytical characterization of the heteroclinic bifurcation is a key result; however, the transition to the PWL approximation in deriving the explicit curve should include a rigorous justification that the qualitative dynamics, including the uniqueness of the limit cycle, are preserved without introducing new bifurcations.
  2. [§5] §5: In the noise analysis, the claim that stochasticity anticipates the onset of the heteroclinic bifurcation leading to earlier co-extinctions requires specification of the noise type and intensity range to ensure the result is not sensitive to particular choices of parameters.
minor comments (2)
  1. Figure captions should explicitly state the parameter values used to generate the phase portraits for reproducibility.
  2. The abstract mentions 'all the phase portraits' but the main text could benefit from a table summarizing the conditions for each portrait.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading, positive assessment of the contribution, and recommendation for minor revision. We address the two major comments point by point below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [§3] §3: The analytical characterization of the heteroclinic bifurcation is a key result; however, the transition to the PWL approximation in deriving the explicit curve should include a rigorous justification that the qualitative dynamics, including the uniqueness of the limit cycle, are preserved without introducing new bifurcations.

    Authors: We agree that a more explicit justification is desirable. The manuscript states that the PWL approximation preserves the qualitative behaviour of the original cubic system, but we acknowledge that the transition step would benefit from additional rigor. In the revised version we will add a dedicated paragraph (or short appendix) that (i) compares the nullclines and vector fields of the two systems, (ii) verifies that the PWL model retains the same number and stability of equilibria, and (iii) confirms, via a combination of index theory and numerical continuation, that no new bifurcations appear and that the unique stable limit cycle persists for the same parameter region. This will directly address the concern while keeping the explicit bifurcation curve as the main practical output. revision: yes

  2. Referee: [§5] §5: In the noise analysis, the claim that stochasticity anticipates the onset of the heteroclinic bifurcation leading to earlier co-extinctions requires specification of the noise type and intensity range to ensure the result is not sensitive to particular choices of parameters.

    Authors: We accept the need for greater precision. The stochastic extension employs additive white noise acting on the resource equation. In the revision we will explicitly state the noise type, the range of intensities examined (small to moderate values for which the anticipation effect is observed), and include a brief sensitivity check showing that the earlier crossing of the heteroclinic threshold remains consistent across this interval. These clarifications will be inserted in §5 together with the existing simulation results. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper conducts a direct mathematical analysis of a cubic planar ODE system, enumerating all phase portraits in the Poincaré disk and deriving the existence of a unique stable limit cycle plus the heteroclinic bifurcation curve from the model equations and their equilibria under stated ecological constraints. The subsequent PWL approximation is applied only after the qualitative results to obtain an explicit formula while preserving the already-established dynamics; it does not redefine or substitute for the original-system claims. No data-fitting occurs, no predictions reduce to fitted inputs by construction, and no load-bearing steps rely on self-citations or imported uniqueness theorems. The chain from vector field to phase-portrait classification to bifurcation loci is independent of the target results.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available; no explicit free parameters, axioms, or invented entities are listed in the provided text.

pith-pipeline@v0.9.0 · 5507 in / 1125 out tokens · 54836 ms · 2026-05-10T07:06:19.030139+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

57 extracted references · 57 canonical work pages

  1. [1]

    A. A. Andronov, A. A. Vitt, and S. E. Khaikin.Theory of oscillators. Pergamon Press, Oxford- New York-Toronto, Ont., 1966. Translated from the Russian by F. Immirzi; translation edited and abridged by W. Fishwick

    work page 1966

  2. [2]

    A. D. Bazykin.Nonlinear dynamics of interacting populations, volume 11 ofWorld Scientific Series on Nonlinear Science. Series A: Monographs and Treatises. World Scientific Publishing Co., Inc., River Edge, NJ, 1998

    work page 1998

  3. [3]

    A. D. Bazykin and F. S. Berezovskaya. Allee effect, lower critical population density and dynamics of predator-prey system.Problemy ekologicheskogo monitoringa i modelirovaniya ekosistem, 2:161–175, 1979

    work page 1979

  4. [4]

    Benoˆ ıt

    E. Benoˆ ıt. Chasse au canard. II. Tunnels—entonnoirs—peignes.Collect. Math., 32(2):77–97, 1981

    work page 1981

  5. [5]

    Benoˆ ıt and J.-L

    E. Benoˆ ıt and J.-L. Callot. Chasse au canard. IV. Annexe num´ erique.Collect. Math., 32(2):115–119, 1981

    work page 1981

  6. [6]

    Berdugo, S

    M. Berdugo, S. K´ efi, S. Soliveres, and F. Maestre. Plant spatial patterns identify alternative ecosys- tem multifunctionality states in global drylands.Nature Ecol. Evol., 1:003, 2017

    work page 2017

  7. [7]

    M. D. Bertness and R. Callaway. Positive interactions in communities.Trends in Ecology & Evolution, 9(5):191–193, 1994

    work page 1994

  8. [8]

    Borsotti and M

    J. Borsotti and M. Sensi. A geometric analysis of the Bazykin–Berezovskaya predator-prey model with Allee effect in an economic framework.Nonlinear Analysis: Real World Applications, 89:104534, June 2026

    work page 2026

  9. [9]

    R. W. Brooker, F. T. Maestre, R. M. Callaway, C. J. Lortie, L. A. Cavieres, G. Kunstler, P. Liancourt, K. Tielb¨ orger, J. M. Travis, F. Anthelme, et al. Facilitation in plant communities: The past, the present, and the future.Journal of Ecology, 96(1):18–34, 2008

    work page 2008

  10. [10]

    T. M. Brooks, R. A. Mittermeier, C. G. Mittermeier, G. A. B. Da Fonseca, A. B. Rylands, W. R. Konstant, P. Flick, J. Pilgrim, S. Oldfield, G. Magin, and C. Hilton-Taylor. Habitat loss and extinc- tion in the hotspots of biodiversity.Conservation Biology, 16(4):909–923, 2002

    work page 2002

  11. [11]

    J.-L. Callot. Chasse au canard. III. Les canards ont la vie br` eve.Collect. Math., 32(2):99–114, 1981

    work page 1981

  12. [12]

    W. Coppel. A survey of quadratic systems.J. Differ. Equations, 2(3):293–304, 1966

    work page 1966

  13. [13]

    R. F. Costantino, R. A. Desharnais, J. M. Cushing, and B. Dennis. Chaotic dynamics in an insect population.Science, 275:389–439, 1997

    work page 1997

  14. [14]

    Dakos, B

    V. Dakos, B. Matthews, A. Hendry, J. Levine, N. Loeuille, J. Norberg, P. Nosil, M. Scheffer, and L. D. Meester. Ecosystem tipping points in an evolving world.Nature Ecology Evolution, 3:355–362, 2019

    work page 2019

  15. [15]

    Dennis, R

    B. Dennis, R. A. Desharnais, J. M. Cushing, and R. F. Costantino. Estimating chaos and complex dynamics in an insect population.Journal of Animal Ecology, 66:704–729, 1997

    work page 1997

  16. [16]

    Diener and M

    F. Diener and M. Diener. Chasse au canard. I. Les canards.Collect. Math., 32(1):37–74, 1981

    work page 1981

  17. [17]

    Diz-Pita, J

    E. Diz-Pita, J. Llibre, and M. V. Otero-Espinar.Study of the nonelementary singular points and the dynamics near the infinity in predator-prey systems, volume 38 ofNonlinear Syst. Complex., pages 103–128. Springer, Cham, 2023

    work page 2023

  18. [18]

    G. Duff. Limit-cycles and rotated vector fields.Annals of Mathematics, 57:15, 1953

    work page 1953

  19. [19]

    Dumortier, J

    F. Dumortier, J. Llibre, and J. C. Art´ es.Qualitative theory of planar differential systems. Universi- text. Springer-Verlag, Berlin, 2006

    work page 2006

  20. [20]

    Dumortier and R

    F. Dumortier and R. Roussarie.Canard Cycles and Center Manifolds. American Mathematical So- ciety: Memoirs of the American Mathematical Society. American Mathematical Society, 1996

    work page 1996

  21. [21]

    C. S. Elton. Fluctuations in the numbers of animals: Their causes and effects.British Journal of Experimental Biology, 2:119–163, 1924

    work page 1924

  22. [22]

    C. S. Elton and M. Nicholson. The 10-year cycle in numbers of the lynx in Canada.Journal of Animal Ecology, 11:215–244, 1924

    work page 1924

  23. [23]

    L. Fahrig. Effects of habitat fragmentation on biodiversity.Annu. Rev. Ecol. Evol. Syst., 34:487–515, 2003

    work page 2003

  24. [24]

    Ferragut and A

    A. Ferragut and A. Gasull. Non-algebraic oscillations for predator-prey models.Publ. Mat, 195:207, 2014

    work page 2014

  25. [25]

    A. F. Filippov.Differential equations with discontinuous righthand sides, volume 18 ofMathematics and its Applications (Soviet Series). Kluwer Academic Publishers Group, Dordrecht, 1988. Trans- lated from the Russian

    work page 1988

  26. [26]

    A. B. Franklin et al. Range-wide declines of northern spotted owl populations in the pacific northwest: A meta-analysis.Biological Conservation, 259:109168, 2021. GLOBAL ANALYSIS OF A RESOURCE-CONSUMER MODEL 35

    work page 2021

  27. [27]

    Gonzalez, J

    A. Gonzalez, J. Lawton, F. Gilbert, T. Blackburn, and I. Evans-Freke. Metapopulation dynamics, abundance, and distribution in a microecosystem.Science, 281:2045–2047, 1998

    work page 2045

  28. [28]

    N. M. Haddad, L. A. Brudvig, J. Clobert, K. F. Davies, A. Gonzalez, et al. Habitat fragmentation and its lasting impact on earth’s ecosystems.Sci. Adv., 1(2):e1500052, 2015

    work page 2015

  29. [29]

    Hanski.The Shrinking World: Ecological Consequences of Habitat Loss

    I. Hanski.The Shrinking World: Ecological Consequences of Habitat Loss. 14 (Excellence in Ecology). International Ecology Institute, 2005

    work page 2005

  30. [30]

    I. Hanski. Habitat loss, the dynamics of biodiversity, and a perspective on conservation.Ambio, 40:248–255, 2011

    work page 2011

  31. [31]

    Hanski and M

    I. Hanski and M. Gilpin. Metapopulation dynamics: brief history and conceptual domain.Biological Journal of the Linnean Society, 42:3–16, 1991

    work page 1991

  32. [32]

    Keitt, M

    T. Keitt, M. Lewis, and R. Holt. Allee effects, invasion pinning, and species’ borders.Am. Nat., 157:203–216, 1997

    work page 1997

  33. [33]

    Krupa and P

    M. Krupa and P. Szmolyan. Extending geometric singular perturbation theory to nonhyper- bolic points-fold and canard points in two dimensions.SIAM Journal on Mathematical Analysis, 33(2):286–314, 2001

    work page 2001

  34. [34]

    Kuznetsov.Elements of Applied Bifurcation Theory

    Y. Kuznetsov.Elements of Applied Bifurcation Theory. Applied Mathematical Sciences. Springer New York, 2008

    work page 2008

  35. [35]

    Y. A. Kuznetsov, S. Rinaldi, and A. Gragnani. One-parameter bifurcations in planar Filippov sys- tems.Int. J. Bifurcation Chaos Appl. Sci. Eng., 13(8):2157–2188, 2003

    work page 2003

  36. [36]

    W. F. Laurance. Habitat destruction: death by a thousand cuts. InConservation Biology for All. Oxford University Press, 2010

    work page 2010

  37. [37]

    R. Levins. Some demographic and genetic consequences of environmental heterogeneity for biological control.Bulletin of the Entomological Society of America, 15:237–240, 1969

    work page 1969

  38. [38]

    Li and Y

    B. Li and Y. Kuang. Heteroclinic bifurcation in the Michaelis–Menten-type ratio-dependent predator- prey system.SIAM J. Appl. Math., 67(5):1453–1464, 2007

    work page 2007

  39. [39]

    Ovaskainen and I

    O. Ovaskainen and I. Hanski. From individual behavior to metapopulation dynamics: unifying the patchy population and classic metapopulation models.Am. Nat., 164:346–377, 2004

    work page 2004

  40. [40]

    L. Perko. Rotated vector fields and the global behavior of limit cycles for a class of quadratic systems in the plane.J. Differ. Equations, 18(1):63–86, 1975

    work page 1975

  41. [41]

    L. M. Perko. A global analysis of the Bogdanov–Takens system.SIAM Journal on Applied Mathe- matics, 52(4):1172–1192, 1992

    work page 1992

  42. [42]

    Sardany´ es and E

    J. Sardany´ es and E. Fontich. On the metapopulation dynamics of autocatalysis: extinction transients related to ghosts.Int. J. Bifurc. Chaos, 120:1–8, 2010

    work page 2010

  43. [43]

    Sardany´ es, F

    J. Sardany´ es, F. Iv´ an, and B. Vidiella. Identifying regime shifts, transients and late warning signals for proactive ecosystem management.Biological Conservation, 290:110433, 2024

    work page 2024

  44. [44]

    Sardany´ es, J

    J. Sardany´ es, J. Piero, and R. Sol´ e. Habitat loss-induced tipping points in metapopulations with facilitation.Pop. Ecol., 61:436–449, 2019

    work page 2019

  45. [45]

    W. M. Schaffer. Stretching and folding in lynx fur returns: Evidence for a strange attractor in nature?American Naturalist, 124:798–820, 1984

    work page 1984

  46. [46]

    Scheffer.Critical Transitions in Nature and Society

    M. Scheffer.Critical Transitions in Nature and Society. Princeton University Press, Princeton, NJ, USA, 2009

    work page 2009

  47. [47]

    Scheffer, S

    M. Scheffer, S. Carpenter, J. A. Foley, C. Folke, and B. Walker. Catastrophic shifts in ecosystems. Nature, 413:591–596, 2001

    work page 2001

  48. [48]

    Seenivasana and P

    R. Seenivasana and P. Paul. Slow–fast dynamics in small trophic chains with habitat loss.Eur. Phys. J. B., 98:105, 2025

    work page 2025

  49. [49]

    Sigmund.Kolmogorov and population dynamics, pages 177–186

    K. Sigmund.Kolmogorov and population dynamics, pages 177–186. Springer Berlin Heidelberg, Berlin, Heidelberg, 2007

    work page 2007

  50. [50]

    D. L. Silvano and M. V. Segalla. Conservation of brazilian amphibians.Conservation Biology, 19(3):653–658, 2005

    work page 2005

  51. [51]

    Soliveres, C

    S. Soliveres, C. Smit, and F. T. Maestre. Moving forward on facilitation research: response to changing environments and effects on the diversity, functioning and evolution of plant communities. Biological Reviews, 90(1):297–313, 2015

    work page 2015

  52. [52]

    R. Sol´ e. Nonequilibrium dynamics in conservation biology: Scales, attractors and critical points. Biological Conservation, 294:110601, 2024

    work page 2024

  53. [53]

    Tilman, R

    D. Tilman, R. May, C. Lehman, and M. Nowak. Habitat destruction and the extinction debt.Nature, 371:65–66, 1994

    work page 1994

  54. [54]

    Vidiella, E

    B. Vidiella, E. Fontich, S. Valverde, and J. Sardany´ es. Habitat loss causes long extinction transients in small trophic chains.Theoretical Ecology, 14(4):Paper No. 641, 661, 2021. 36 T. MAYAYO, J. SARDANY ´ES, AND J. TORREGROSA

    work page 2021

  55. [55]

    J. Wang, J. Shi, and J. Wei. Predator-prey system with strong Allee effect in prey.J. Math. Biol., 62(3):291–331, 2011

    work page 2011

  56. [56]

    With and A

    K. With and A. King. Dispersal success on fractal landscapes: a consequence of lacunarity thresholds. Landscape Ecology, 14:73–82, 1999

    work page 1999

  57. [57]

    Zegeling, , H

    A. Zegeling, , H. Wang, and G. Zhu. Uniqueness of limit cycles in a predator-prey model with sigmoid functional response.Journal of Nonlinear Modeling and Analysis, 5(4):790–802, 2023. AppendixA.Poincar ´e compactification The Poincar´ e compactification is a common tool for studying the limit sets of planar polynomial differential systems. To carry out t...

    work page 2023