Pith. sign in

REVIEW 2 major objections 4 minor 17 references

Invasion Fronts in Shifting Habitats and Competition Systems: A Hamilton-Jacobi Approach and Nonlocal Effects

T0 review · 2 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper surveys Hamilton–Jacobi methods for invasion spreads in shifting habitats and proves a new result: at intermediate habitat-shift speeds the invasion front is pinned to within bounded distance of the moving habitat edge, while a sl

desk verdict Useful survey with a genuinely new but modest front-location theorem; the proof has a real gap in Lemma 11.8 that should be fixed before publication. read the letter →

arxiv 2607.29001 v1 pith:WIWYV2DK submitted 2026-07-31 math.AP

King-Yeung Lam , Chang-Hong Wu , Xiao Yu This is my paper · ORCID
classification math.AP MSC 35K5792D2535F2135B4035Q92
keywords reaction-diffusionequationsinvasionfrontsshiftinghabitatsHamilton-JacobispreadingspeedslogarithmiccorrectionnonlocalpullingFisher-KPPequation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Drawing together recent work on spreading phenomena in reaction–diffusion models of ecological invasion, this paper develops and applies the Hamilton–Jacobi (WKB) approach to environments whose quality shifts in space and time. Its new contribution concerns a Fisher–KPP population in which the habitat ahead of the front is less favorable (growth rate 1−a) and the favorable region retreats along a curve X(t)=βt−η log(t+1). The authors prove that the δ-level set of the invasion front, ξδ(t), satisfies sharp two-sided bounds: for slow habitat shifts (β<2√(1−a)) the front lags by (3/(2λ_min)) log t, the classical logarithmic delay; for intermediate shifts (2√(1−a)<β<2) the front stays within O(1) of the habitat edge; and for any β<2 with η≥0 it never lags behind the habitat edge. This pins down the precise front location in a regime where only speeds were known.

What carries the argument

The key object is the invasion front's level set ξδ(t) together with the Hamilton–Jacobi rate function w(t,x)=lim_{ε→0} −ε log u(t/ε,x/ε), whose zero set gives the occupied region. The new proof works in the moving frame y=x−βt+η log(t+1), where the equation becomes v_t=v_yy+(β−η/(t+1))v_y+v(g(y)−v) with g(y)=1 for y≤0 and g(y)=1−a for y>0. The upper bound in (ii) is obtained by gluing a non-minimal forced wave for the shifted problem to an exponential tail, producing a generalized supersolution; the lower bound in (iii) uses a compactly supported sinusoidal subsolution in the region where the growth rate is the larger one. These comparisons convert the shifting-habitat problem into a sequen

What would settle it

Solve the moving-frame equation v_t=v_yy+(β−η/(t+1))v_y+v(g(y)−v) with g(y)=1 for y≤0 and g(y)=1−a for y>0, using compactly supported initial data satisfying (62), for a fixed β∈(0,2) and η≥0. Track ζ_{δ1}(t)−ζ_{δ2}(t) for two levels δ1,δ2∈(0,1−a) up to large t. If this difference is unbounded as t→∞, the level-set gap assertion fails and the lower bound in Theorem 11.3(iii) does not hold for all δ.

Watch

Extended reading notes

Core claim

The paper's central new result, Theorem 11.3, establishes that for the reaction–diffusion equation u_t = u_xx + u(r(t,x)−u) with r(t,x)=1 for x ≤ X(t) and r(t,x)=1−a for x > X(t), X(t)=βt−η log(t+1), the level set ξδ(t)=sup{x:u(t,x)≥δ} satisfies: (i) if β<2√(1−a), then ξδ(t)=2√(1−a)t − (3/(2λ_min)) log t + O(1) with λ_min=√(1−a); (ii) if 2√(1−a)<β<2, then ξδ(t)≤βt−η log(t+1)+O(1); and (iii) if 0<β<2 and η≥0, then ξδ(t)≥βt−η log(t+1)−O(1). In words, when the favorable region retreats slowly the classical logarithmic delay persists; when it retreats at an intermediate speed the invasion front is confined to a bounded neighbourhood of the moving habitat edge; and a logarithmic retarding of the

Load-bearing premise

The load-bearing premise is that, for the solution in the moving frame, the distance between any two δ-level sets ζ_{δ1}(t) and ζ_{δ2}(t) (with δ1,δ2∈(0,1−a)) remains bounded uniformly in time; this is stated without proof in Lemma 11.8 and is needed to pass from small δ to all δ in the lower-bound statement.

Editorial extensions

If this is right

  • If the habitat shift speed is below the minimal KPP speed of the less favorable region, the invasion front's position is 2√(1−a)t − (3/(2λ_min)) log t + O(1): the classical logarithmic delay survives the moving environment.
  • If the habitat shift speed lies strictly between 2√(1−a) and 2, the front cannot outrun the habitat edge: ξδ(t) is bounded above by X(t)+O(1).
  • If the habitat edge is retarding (η≥0) and shifts at any speed β<2, the front cannot lag behind it: ξδ(t)≥X(t)−O(1). Together with (ii) this pins the front to within O(1) of the edge when η=0 and 2√(1−a)≤β<2.
  • When the habitat shift is linear (η=0) and β∈[2√(1−a),2), the front position is βt+O(1): the invasion speed equals the environmental shift speed exactly.
  • These bounds refine Hamilton–Jacobi speed results into front-location asymptotics, connecting the spreading-speed picture to the logarithmic-correction literature.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural extension would be to prove profile convergence: in the intermediate regime the solution should approach a forced traveling wave with a bounded (possibly t-dependent) phase, analogous to the classical logarithmic convergence result.
  • The same comparison strategy may work for smooth monotone growth profiles r(x−c1t) instead of a step function, but the level-set gap assertion used in Lemma 11.8 would need to be re-established in that generality.
  • The paper's open case β=2, 1/2≤η<3/2 could be probed numerically: if the supercritical formula continues to hold, the front should exhibit the coefficient (1/λ∗)(3/2−√a η) log t; a direct simulation of the level set would test this.
  • If the unproved assertion that δ-level sets in the moving frame remain a bounded distance apart fails, Theorem 11.3(iii) would only be valid for small δ; this is checkable by computing ζ_{δ1}(t)−ζ_{δ2}(t) for two different levels.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. This manuscript is a survey of spreading phenomena in reaction-diffusion equations in shifting habitats, organized around the Hamilton-Jacobi approach. It reviews the Shigesada-Kawasaki conjecture, the reduction of competition systems to scalar shifting-habitat models, local versus nonlocal selection of spreading speeds, flux-limited viscosity solutions for non-monotone environments, spreading in competition and predator-prey systems, and entire solutions. The paper's original contribution is in Section 11.2: Theorem 11.3 gives Bramson-type logarithmic corrections for the level set xi_delta(t) when the region ahead of the shifting interface is less favorable, with Lemmas 11.6-11.8 providing the proof. The claimed result pins the level set to within O(1) of the interface in the intermediate regime 2*sqrt(1-a) < beta < 2 and gives the classical (3/(2*lambda_min)) log t delay when beta < 2*sqrt(1-a).

Significance. The survey is well-organized and serves a useful purpose in consolidating the recent literature, with clear statements of the Hamilton-Jacobi framework, nonlocal pulling, and open questions. If Theorem 11.3 is fully established, it is a significant new result: it extends Bramson's logarithmic correction to shifting habitats with a worse region ahead and shows that in the forced-wave regime the invasion front is pinned within bounded distance of X(t). The proof strategy is natural, combining comparison with homogeneous KPP dynamics, a forced-wave supersolution, and a stationary subsolution, and Lemmas 11.6 and 11.7 are convincing. However, the manuscript has a load-bearing gap in Lemma 11.8 that prevents Theorem 11.3(iii) and Corollary 11.4 from being fully supported as written.

major comments (2)
  1. [Section 11.2, Lemma 11.8 (final sentence)] The proof of Lemma 11.8 constructs a stationary subsolution v with amplitude kappa*phi(y*) and explicitly yields the lower bound zeta_delta(t) >= -C only for delta < kappa*phi(y*). The final sentence then asserts (84) for every delta in (0,1-a) 'due to the fact that zeta_{delta_1}(t)-zeta_{delta_2}(t)=O(1)', but no proof, citation, or argument is supplied. This fact is exactly what is needed to pass from the small-threshold construction to the full statement of Theorem 11.3(iii). It is not an immediate consequence of the comparison principle already used, because the solution has not been shown to converge to a single translate of a forced wave. Since (iii) is the core lower bound in the regime 2*sqrt(1-a)<beta<2, the theorem currently lacks support for all delta. Please provide a proof or a precise reference.
  2. [Section 11.2, Corollary 11.4] The proof of Corollary 11.4 says that 'the argument of (ii) in Theorem 11.3 can be suitably modified' to cover beta=2*sqrt(1-a), eta=0. This modification is not shown. Lemma 11.7 requires beta>2*sqrt(1-a) to obtain lambda_0>0 and a non-minimal forced wave U_beta with exponential decay; at beta=2*sqrt(1-a), lambda_0=0 and the construction as written breaks down. Since the corollary is the sharp 'bounded distance from X(t)' statement at the boundary of the forced-wave regime, the required modification needs to be spelled out in detail.
minor comments (4)
  1. [Title/header] The running title contains a typo: 'INV ASION' should be 'INVASION'.
  2. [Section 11.2, Lemma 11.8] The lower-case letter v is used both for the moving-frame solution v(t,y) and for the stationary subsolution v(y). This is confusing; consider renaming the subsolution (e.g., w or phi).
  3. [Section 11.2, Lemma 11.8] The support condition 'y in (-infinity, -L_0 + pi*ell) subset (-infinity, 0)' relies on L_0 > pi*ell; this is chosen, but the sentence 'By the choice of L_0' could be more explicit about the ordering of the constants.
  4. [Section 11.1, Theorem 11.2 and following display] The notation for the logarithmic correction in the transition case beta=2*(sqrt(a)+sqrt(1-a)) is not fully consistent: Theorem 11.2(iii) defines m_q(t) with q=-3/2+eta*sqrt(a), but the display after the theorem gives the special case eta=0 with coefficient -5/(4*lambda_min). Unifying these formulas would improve readability.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the new Theorem 11.3 derivation reduces to external classical results (Bramson [20,61], Berestycki–Fang [12]) and comparison arguments, not to the authors' own speed formulas; the only flagged weakness is an unproved δ-extension in Lemma 11.8, which is a rigor gap rather than circularity.

full rationale

The paper is a survey plus one new theorem (Theorem 11.3). The theorem's proof uses Lemmas 11.6–11.8. Lemma 11.6 compares with the homogeneous Fisher–KPP equation and uses classical Bramson asymptotics [20,61]. Lemma 11.7 uses the external forced-wave result [12, Theorem 1.3(iii)] to build a supersolution. Lemma 11.8 constructs an explicit stationary subsolution and uses the comparison principle directly; no fitted parameter is involved. The authors' own prior speed formulas, e.g., (65) and (70), are cited as background, not as inputs to the new proof. Self-citations in the survey, including [87,84,86,52,103], are attributions of published results and are not used to replace a proof in the new derivation. The one substantive weakness is in Lemma 11.8: after proving (84) only for δ∈(0,κφ(y*)), the text states without proof that the estimate 'still holds for any δ∈(0,1−a) ... due to the fact that ζ_{δ1}(t)−ζ_{δ2}(t)=O(1)'. That assertion is load-bearing for Theorem 11.3(iii) and Corollary 11.4, and it is not derived; however, it is an omitted argument, not a circular reduction. It should be assessed as a rigor/correctness risk. There is no construction by which an input parameter equals the output prediction, and no uniqueness claim from the authors' prior work is used to force the conclusion. Score 2 reflects minor self-citations in the survey and the flagged gap, not a circular derivation.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central new result introduces no fitted parameters and no new entities. It depends on standard comparison principles, Bramson's classical log-delay theorem, and the existence of non-minimal forced waves from [12]. One auxiliary fact about level-set differences is asserted without support.

assumptions (4)
  • standard math Generalized comparison principle for piecewise C^1 super/subsolutions of the shifted equation (74), as in [83, Remark 1.1.2].
    Used in Lemmas 11.7 and 11.8 to compare v(t,y) with the constructed gluing supersolution ⅕v and subsolution v. The principle is standard but is imported rather than proved here.
  • standard math Bramson's logarithmic correction for homogeneous Fisher–KPP level sets: ξδ(t)=c*t − (3/(2λ*)) log t + O(1) [20,61].
    Used in Lemma 11.6 for both lower and upper bounds when β<2√(1−a).
  • domain assumption Existence of a non-minimal forced wave Uβ for (78) with Uβ(−∞)=1, Uβ(+∞)=0, U′β<0, and Uβ(y)∼e^{−λ0 y} for β∈(2√(1−a),2), quoted from [12, Theorem 1.3(iii)].
    This profile and its exponential decay rate are load-bearing for the supersolution construction in Lemma 11.7; without it the upper bound in Theorem 11.3(ii) would fail.
  • domain assumption Level-set positions for nearby δ differ by O(1): ζδ1(t)−ζδ2(t)=O(1) for δ1,δ2∈(0,1−a).
    Used at the end of Lemma 11.8 to extend the lower bound from small δ to all δ∈(0,1−a). This is asserted without proof or citation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Invasion Fronts in Shifting Habitats and Competition Systems: A Hamilton-Jacobi Approach and Nonlocal Effects." pith.science (2026). https://pith.science/paper/WIWYV2DK

@misc{pith2026260729001,
  author       = {Pith},
  title        = {Pith review of: Invasion Fronts in Shifting Habitats and Competition Systems: A Hamilton-Jacobi Approach and Nonlocal Effects},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WIWYV2DK}},
  note         = {Machine review of arXiv:2607.29001}
}
read the original abstract

We review recent developments in the study of spreading phenomena in reaction--diffusion equations arising from ecological invasion models. Motivated by the conjecture of Shigesada and Kawasaki on staged invasions, we discuss how competition systems can lead to effective scalar models with shifting habitats. We present the Hamilton--Jacobi approach for determining spreading speeds and revisit the result of {Li--Bewick--Shang--Fagan} (2014) from this perspective. We then describe the emergence of nonlocally pulled fronts when the shifting habitat connects regions of distinct positive growth rates. Recent results including the works of Lam--Yu (2022) and Lam--Nadin--Yu (2025) are surveyed. We also discuss spreading phenomena in competition systems, predator-prey systems, and the existence of various classes of entire solutions. Finally, we discuss the logarithmic correction for invasion waves in moving environments and prove a new result.

Figures

Figures reproduced from arXiv: 2607.29001 by the authors.

Figure 1
Figure 1. A piecewise linear trajectory that first travels ahead of the front γ(s) = cnlps and then relaxes to the front point (t, cnlpt) In all cases, for compactly supported initial data, the spreading speed is always locally selected. This is the case for compactly supported initial data. A recent study on the nonlocal diffusion problem by Tao et al. [128] revealed a novel nonlocal speed enhancement mechanism. As we have s… view at source ↗
Figure 2
Figure 2. The hyperbolic limits of three families of entire solutions constructed in [85] in the form of diverging waves that connects E1 as x → −∞ and 0 as x → +∞. The second invasion front where E2 is destabilized is nonlinearly selected. Here, 0 = (0, 0), E1 = (1, 0), E2 = (0, 1) and E ∗ = ( 1−a 1−ab , 1−b 1−ab ). Open Question 4. Characterize the set of all entire solutions to (2). While different classes of entire soluti… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 4 linked inside Pith

  1. [5]

    in Math., 30 (1978), pp

    ,Multidimensional nonlinear diffusion arising in population genetics, Adv. in Math., 30 (1978), pp. 33–76. [6]M. A very, M. Holzer, and A. Scheel,Pulled fronts are not (just) pulled, 2025. Preprint, arXiv:2508.14864. [7]G. Barles and E. Chasseigne,On modern approaches of Hamilton-Jacobi equations and control problems with discontinuities—a guide to theory...

    arXiv 1978

  2. [19]

    ,The Bramson delay in the non-local Fisher-KPP equation, Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire, 37 (2020), pp. 51–77. [20]M. Bramson,Convergence of solutions of the Kolmogorov equation to travelling waves, Mem. Amer. Math. Soc., 44 (1983), pp. iv+190. [21]C. Br ¨andle and E. Chasseigne,Large deviations estimates for some non-local equations. Fas...

    2020

  3. [36]

    ,Spreading speeds for multidimensional reaction-diffusion systems of the prey-predator type, Calc. Var. Partial Differential Equations, 58 (2019). Paper No. 137, 34. [37]A. Ducrot and Z. Jin,Spreading speeds for time heterogeneous prey-predator systems with nonlocal diffusion on a lattice, J. Differential Equations, 396 (2024), pp. 257–313. [38]L. C. Evan...

    2019

  4. [45]

    ,Geometric optics approach to reaction-diffusion equations, SIAM J. Appl. Math., 46 (1986), pp. 222–232

    1986

  5. [46]

    Probab., 19 (1991), pp

    ,Coupled reaction-diffusion equations, Ann. Probab., 19 (1991), pp. 29–57. [47]J. Garnier, T. Giletti, F. Hamel, and L. Roques,Inside dynamics of pulled and pushed fronts, J. Math. Pures Appl. (9), 98 (2012), pp. 428–449. [48]J. Gertner and M. I. Fre ˘ıdlin,The propagation of concentration waves in periodic and random media, Dokl. Akad. Nauk SSSR, 249 (19...

    arXiv 1991

  6. [58]

    ,Reaction-diffusion problems in cylinders with no invariance by translation. II. Mono- tone perturbations, Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire, 14 (1997), pp. 555–596. [59]F. Hamel and N. Nadirashvili,Entire solutions of the KPP equation, Comm. Pure Appl. Math., 52 (1999), pp. 1255–1276. [60]F. Hamel and N. Nadirashvili,Travelling fronts and e...

    1997

  7. [67]

    Nonlinear Sci., 35 (2025)

    ,Global dynamics for a class of nonlocal evolution systems in a periodic shifting environment, J. Nonlinear Sci., 35 (2025). Paper No. 48, 35

    2025

  8. [68]

    Differential Equations, 432 (2025)

    ,Spatio-temporal dynamics for cooperative reaction-diffusion systems with asymptotic annihilation, J. Differential Equations, 432 (2025). Paper No. 113234, 25. [69]C. Hu, J. Shang, and B. Li,Spreading speeds for reaction-diffusion equations with a shifting habitat, J. Dynam. Differential Equations, 32 (2020), pp. 1941–1964. [70]M. Iida, R. Lui, and H. Nin...

    2025

Show all 17 references
  1. [72]

    ,Quasi-convex Hamilton-Jacobi equations posed on junctions: the multi-dimensional case, Discrete Contin. Dyn. Syst., 37 (2017), pp. 6405–6435. [73]H. Ishii,Hamilton-Jacobi equations with discontinuous Hamiltonians on arbitrary open sets, Bull. Fac. Sci. Engrg. Chuo Univ., 28 (...

    2017

  2. [74]

    2074 of Lecture Notes in Math., Springer, Heidelberg, 2013, pp

    ,A short introduction to viscosity solutions and the large time behavior of solutions of Hamilton-Jacobi equations, in Hamilton-Jacobi equations: approximations, numerical analysis and applications, vol. 2074 of Lecture Notes in Math., Springer, Heidelberg, 2013, pp. 111–249. ...

    arXiv

  3. [103]

    Dif- ferential Equations, 271 (2021), pp

    ,Stacked invasion waves in a competition-diffusion model with three species, J. Dif- ferential Equations, 271 (2021), pp. 665–718. [104]S. Liu, Q. Liu, and K.-Y. Lam,Asymptotic spreading of interacting species with multiple fronts II: Exponentially decaying initial data, J. Di...

    2021

  4. [112]

    ,Invasion speed of a predator-prey system, Appl. Math. Lett., 74 (2017), pp. 46–51. [113]R. Peng, C.-H. Wu, and M. Zhou,Sharp estimates for the spreading speeds of the Lotka- Volterra diffusion system with strong competition, Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire, 38...

    2017

  5. [129]

    Preprint, arXiv:2604.14573

    ,Propagation dynamics for nonlocal dispersal predator-prey systems in shifting habi- tats: A hamilton-jacobi approach, 2026. Preprint, arXiv:2604.14573. [130]W. Tao, W.-T. Li, and J.-W. Sun,A Hamilton-Jacobi approach to neural field equations, J. Differential Equations, 422 (2...

    2026 arXiv

  6. [149]

    ,Global dynamics of evolution systems with asymptotic annihilation, Discrete Contin. Dyn. Syst., 43 (2023), pp. 2693–2720

    2023

  7. [150]

    Differential Equations, 436 (2025)

    ,Spatio-temporal dynamics for non-monotone semiflows with limiting systems having spreading speeds, J. Differential Equations, 436 (2025). Paper No. 113267, 52

    2025

  8. [151]

    Differential Equations, 469 (2026)

    ,Spatio-temporal dynamics for evolution systems with asymptotic translation invari- ance, J. Differential Equations, 469 (2026). Paper No. 114349, 37. [152]Y. Yuan, Y. W ang, and X. Zou,Spatial dynamics of a lotka-volterra model with a shifting habitat, Discrete Contin. Dyn. S...

    2026

  9. [2026]

    Preprint, arXiv:2604.27556. [119]F. Rothe,Convergence to pushed fronts, Rocky Mountain J. Math., 11 (1981), pp. 617–633. 42 K.-Y. LAM, C.-H. WU AND X. YU [120]R. B. Salako,Invasion entire solutions for two-species diffusive monostable competitive systems, Nonlinear Anal. Real ...

    1981 arXiv

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.