A free boundary problem for spreading under shifting climate

Xianfa Song; Xinan Hao; Yihong Du; Yuanyang Hu

A free boundary problem for spreading under shifting climate

T0 review · 0 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read A single critical speed $c_0$ completely determines whether a species spreading into a shifting, improving climate vanishes or spreads, and what profile it takes.

desk verdict Solid extension of the Du-Wei-Zhou shifting-climate free boundary model to the favourable-shift case; the forced-speed semi-wave selection rule is new and the proofs essentially check out, modulo a repairable sign slip in Lemma 3.8. read the letter →

arxiv 1908.04041 v1 pith:M4DL7CPJ submitted 2019-08-12 math.AP

Yuanyang Hu , Xinan Hao , Xianfa Song , Yihong Du This is my paper · ORCID

classification math.AP MSC 35K5735R3592D25

keywords diffusivelogisticequationfreeboundaryproblemspreadingandvanishingdichotomyshiftingclimatesemi-wavecriticalspeedasymptoticprofileinvasivespecies

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

The paper studies a species whose habitat interval expands at a moving front while the climate shifts from lethal to livable at speed $c$. It claims that the eventual outcome is always one of two possibilities—extinction or unbounded spread—and that the boundary between them is a single critical speed $c_0$ computed from the homogeneous favourable habitat. When spread occurs below $c_0$, the front settles into a travelling profile that lags the climate edge by a fixed distance $L_0$; at or above $c_0$, the front moves at the classical speed $c_0$ with the classical semi-wave profile. It also gives a sharp initial-size rule: populations starting in an interval at least $\frac{\pi}{2}\sqrt{d/a}$ long always spread.

What carries the argument

The auxiliary semi-wave problem (1.5), $$-dv''-cv'=A(x)v-$bv^{2}$,\quad -\infty<x<L,\ v(L)=0,$$ has a unique positive solution $v_L$ for each $L\ge 0$, with $v_L(-\infty)=a/b$ and $v_L'<0$. The map $L\mapsto -\mu(A(L))v_L'(L)$ is strictly increasing, so there is a unique $L_0$ with $-\mu(A(L_0))v_{L_0}'(L_0)=c$; $L_0=0$ exactly when $c=c_0$, where $c_0$ is the classical spreading speed and $v_0=q_{c_0}$. This monotone family supplies the forced profiles and the comparison barriers used to control $h(t)-ct$ and to extract convergent limits along time shifts.

What would settle it

Solve (1.5) numerically for a function $A$ satisfying (1.3) but with a non-monotone transition on $[0,l_0]$, and plot $-\mu(A(L))v_L'(L)$; if this curve crosses the level $c$ more than once for some $0<c<c_0$, the uniqueness of $L_0$ fails. Alternatively, simulate (1.4) in that case and check whether two different initial data produce different limits of $h(t)-ct$.

Watch

Extended reading notes

Core claim

The central discovery is a complete classification of the long-time behaviour of the free-boundary logistic equation (1.4), where the environment changes from unfavourable ($A=a<0$) to favourable ($A=a_0>0$) and the front moves by the Stefan-type condition $h'(t)=-\mu(A(h(t)-ct))u_x(h(t),t)$. The paper proves a spreading–vanishing dichotomy: every solution either has $h(t)\to h_\infty<\infty$ with $u\to 0$, or $h(t)\to\infty$. In the spreading case the asymptotic profile depends on a critical speed $c_0$ from the homogeneous problem. If $0<c<c_0$, then $h(t)-ct\to L_0$ and $u(\cdot,t)-v_{L_0}(\cdot+L_0-h(t))\to 0$ in $L^\infty$, with $v_{L_0}$ the unique semi-wave of (1.5); if $c=c_0$ or $c>c_0$, then $h(t)-c_0t$ converges to a constant and $u(\cdot,t)-q_{c_0}(\cdot-h(t))\to 0$, where $q_{c_0}$ is the classical semi-wave. The paper also establishes a sharp initial-range threshold for the spread-vanish alternative.

Load-bearing premise

The whole classification rests on the auxiliary semi-wave problem having a unique solution $v_L$ and a unique $L_0$ where $-\mu(A(L))v_L'(L)=c$; this uniqueness relies on $A$ and $\mu$ being monotone in the assumed directions.

Editorial extensions

If this is right

For $0<c<c_0$, a spreading population's front converges to $h(t)-ct\to L_0$, so the invasion lags the shifting habitat edge by a fixed distance and the density approaches the forced semi-wave $v_{L_0}$.
For $c\ge c_0$, the front asymptotically moves at the homogeneous-environment speed $c_0$ rather than the climate speed, with profile $q_{c_0}$; in particular a faster climate shift does not speed up the invasion.
Every solution either vanishes or spreads; there is no intermediate state, and vanishing implies the population density decays to zero uniformly over the shrinking range.
If the initial range is at least $\frac{\pi}{2}\sqrt{d/a}$, vanishing is impossible, so a sufficiently large starting habitat guarantees spread regardless of the climate speed.
For smaller initial ranges with initial density $\sigma\varphi$, there is a threshold $\sigma_0$: densities below it vanish, above it spread, with $\sigma_0=+\infty$ left as an open possibility.

Reading between the lines

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

Because $L_0$ is characterized by a single monotone equation, the model predicts a directly measurable lag: a spreading front below the critical speed should trail the moving habitat edge by roughly $L_0$ units; time series of range edges could test this without any parameter fitting beyond $c$ and $c_0$.
The monotone dependence in Proposition 1.1 suggests that as $c$ rises to $c_0$, $L_0$ shrinks to zero and the forced semi-wave $v_{L_0}$ degenerates into $q_{c_0}$, so the two regimes connect continuously; the paper does not prove this convergence explicitly.
If the same classification were attempted with a non-monotone $A$ or $\mu$, multiple $L_0$ values could appear and the front could in principle select different lags from different initial data, so the dichotomy may be a special property of monotone environmental gradients.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Worth knowing: this is a competent, careful follow-up to Du-Wei-Zhou [9] that covers the opposite climate-shift sign (a>0>a0) and adds a climate-dependent Stefan coefficient mu(A(h-ct)). The main novelty is the asymptotic selection rule: for 0L0; for c>=c0 the classical c0 semi-wave wins. That classification is genuinely new and is not present in [9] or the other cited works. The paper also delivers the auxiliary semi-wave theory in Proposition 1.1, which is the load-bearing piece, and it is proved, not imported: existence, uniqueness, monotonicity of v'_L(L), and existence/uniqueness of L0 are all established. The comparison arguments in The proofs of Theorems 3.7 and 3.9 are structurally sound, and the paper is honest about the open threshold question in Remark 1.6. Delegating several standard estimates to earlier published papers is acceptable practice here. The citation pattern is not a problem; the paper builds on the authors' own prior work, but the cited items are published tools, not unpublished claims. Soft spots: the monotonicity assumptions on A and mu are explicit and are doing real work in Proposition 1.1, but the paper acknowledges this in the model statement. Lemma 3.8 contains a sign slip: after uniform convergence of u(0,t) to a/b, the text asserts u(0,t)>a/b+epsilon for large t, which cannot be right. The intended lower-comparison goes through with a/b-epsilon (or a/(b+epsilon)), so the lemma and its use in Theorem 3.9 are not endangered. Proposition 1.1 also asserts continuous dependence of v'_L(L) on L without proof; standard elliptic regularity covers it, but the paper should say so. These are minor blemishes in a paper whose central argument holds up. This is for mathematical biology and PDE readers working on Stefan-type free boundaries and shifting environments. It deserves a serious referee; the classification is a solid advance and the errors are repairable. I would send it out.

Referee Report

0 major / 4 minor

Summary. This paper analyzes a free boundary problem for a diffusive logistic equation in a habitat that shifts to a favorable environment at speed c>0. The model (1.4) is the Du-Wei-Zhou model but with the free boundary condition h'(t) = -mu(A(h(t)-ct)) u_x(h(t),t), where mu depends monotonically on the resource function A. The main theorems (Theorems 1.2-1.4) establish a spreading-vanishing dichotomy and, in the spreading case, sharp asymptotics: for 0<c<c0 the front position satisfies h(t)-ct -> L0 and the solution converges to the forced semi-wave v_{L0}(.+L0-h(t)); for c>=c0 the profile is the classical semi-wave q_{c0} with speed c0. Theorem 1.5 gives initial-data threshold criteria for vanishing and spreading. The proofs rely on an auxiliary elliptic problem (Proposition 1.1) that constructs a unique forced semi-wave and a unique L0 satisfying -mu(A(L0))v'_L0(L0)=c, plus comparison arguments and compactness/limiting arguments.

Significance. The paper gives a complete asymptotic classification for a natural variant of the shifting-climate free boundary model, complementing earlier works that treated the unfavorable-shift case. The identification of the forced semi-wave as the spreading profile when c<c0 is a new structural result, and the sharp convergence statements are strong. The proofs are detailed and largely self-contained, using standard tools such as comparison principles, upper/lower solutions, and parabolic compactness. The explicit monotonicity assumptions on A and mu are exactly what make the auxiliary profile unique, so the classification is coherent and non-circular. The main results make precise, falsifiable predictions about the asymptotic position and shape of the front, and the manuscript ships with no ad hoc parameters fitted to data.

minor comments (4)

[Section 3.2, Lemma 3.8] In the proof of Lemma 3.8, the inequality 'u(0,t) > a/b + epsilon for t >= T' is incompatible with the uniform convergence (3.5) to a/b; the correct lower bound is, for example, u(0,t) > a/(b+epsilon) for large t (or u(0,t) > a/b - epsilon), which suffices for the subsequent comparison with u_epsilon. With this correction, the lower-solution argument and the conclusion of the lemma go through unchanged.
[Section 3.1, Lemma 3.6] In the proof of Lemma 3.6, the displayed chain '-mu(A(H))V'_H(H) = -mu(A(H*))V'_H*(H*) > mu(A(L0))V'_L0(L0) = c' is missing a minus sign on the right-hand side; it should read '-mu(A(H*))V'_H*(H*) > -mu(A(L0))V'_L0(L0) = c'. The intended comparison is clear and the argument is otherwise correct.
[Section 2.4, Proposition 1.1(iii)] In Proposition 1.1(iii), the existence of L0 invokes 'continuous dependence of mu(A(L))v'_L(L) on L' without proof; since this is a standard elliptic regularity consequence, please add a brief justification (as is done for the convergence in Lemma 3.1) so the proof is fully rigorous.
[Throughout] There are several typographical errors: 'Liptschitz' should be 'Lipschitz' (Introduction); in Lemma 3.3 the definition 'B(t) := min_{t in [0,+infty)} {h(t), ct}' should use a different variable inside the min (e.g., B(t) := min{h(t), ct}); and in several chain inequalities the quantity mu(A(L))v'_L(L) appears without the leading minus sign, which can confuse the reader even though the context makes the intended sign clear.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the classification is proved from Proposition 1.1 and published comparison/spreading results; no prediction is fitted to data or defined by the target limit.

full rationale

This paper does not exhibit circular reasoning. The critical speed c0 and semi-wave q_c0 are defined from the homogeneous Stefan-type problem via published results [3,4,7], and Proposition 1.1 independently constructs the forced semi-wave family v_L, proves its uniqueness and monotonicity, and defines L0 as the unique point where -mu(A(L))v'_L(L)=c. The asymptotic limits h(t)-ct -> L0 and u(.,t) -> v_L0(.+L0-h(t)) are then proved in Theorems 3.7 and 3.9 by compactness, comparison, and limiting-profile arguments; L0 is not fitted from the solution, and the limit is shown to satisfy the defining equation only after the proof, not by construction. The paper relies heavily on earlier works by the same group (notably [4], [6], [7], [8]) for standard existence, comparison, and homogeneous-spreading results, but these are published mathematical tools external to the present classification claim and are not used to smuggle in the target conclusion. The one concrete blemish, the inequality u(0,t)>a/b+epsilon in Lemma 3.8, is a repairable sign error and does not turn any prediction into an input by construction. Therefore no circular step is identified.

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

No free parameters are fitted: c, d, a, b, a0, l0, h0, mu, and the initial data are specified inputs, and c0 is derived from the homogeneous semi-wave problem. The paper introduces no new physical entities; the forced semi-wave v_L is a mathematical construction used to describe the asymptotic profile. The main proofs rely on established PDE comparison and compactness arguments and on prior free-boundary results by the authors and others.

assumptions (5)

standard math Local existence, uniqueness, and uniform estimates for free-boundary problems of type (1.4), stated as Theorem 2.1, Lemma 2.2, and Theorem 2.3.
The paper proves these by arguments 'similar to' [4] with corrections from [20], and does not reproduce the full proofs.
standard math Comparison principles for the free boundary problem and for logistic equations on unbounded domains, including Lemma 2.1 of [6].
These are used throughout Sections 2-6 for all supersolution and subsolution arguments.
standard math Homogeneous free-boundary spreading speed and profile results from [4] and [7], including the existence of c0 and q_c0 and the limit h(t)-c0t tending to a constant.
Theorem 1.2 of [7] supplies the homogeneous asymptotic used in Sections 4-5.
standard math Classification of entire solutions with moving front for the homogeneous equation, taken from Section 4.2 of [8].
Used in Theorem 1.3, case (III), to identify the limiting entire solution as q_c0.
domain assumption The structural assumptions on A and mu: A is Lipschitz, equals a for nonpositive argument and a0 for argument at least l0, and is strictly monotone on [0,l0] with a>0>a0; mu is continuous and increasing on [a0,a].
This monotonicity ensures the unique forced semi-wave L0 and the ordering arguments in Lemmas 3.1, 3.3, 3.5, and 3.6.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A free boundary problem for spreading under shifting climate." pith.science (2026). https://pith.science/paper/M4DL7CPJ

@misc{pith2026190804041,
  author       = {Pith},
  title        = {Pith review of: A free boundary problem for spreading under shifting climate},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/M4DL7CPJ}},
  note         = {Machine review of arXiv:1908.04041}
}

abstract

In this paper we consider a free boundary problem which models the spreading of an invasive species whose spreading is enhanced by the changing climate. We assume that the climate is shifting with speed c and obtain a complete classification of the long-time dynamical behaviour of the species. The model is similar to that in [9] with a slight refinement in the free boundary condition. While [9], like many works in the literature, investigates the case that unfavourable environment is shifting into the favourable habitat of the concerned species, here we examine the situation that the unfavourable habitat of an invasive species is replaced by a favourable environment with a shifting speed c. We show that a spreading-vanishing dichotomy holds, and there exists a critical speed$c_0$ such that when spreading happens in the case $c < c_0$, the spreading profile is determined by a semi-wave with forced speed c, but when $c \geq c_0$, the spreading profile is determined by the usual semi-wave with speed $c_0$.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 20 canonical work pages

[9]

Y. Du, L. Wei, and L. Zhou, Spreading in a shifting environ ment modeled by the diﬀusive logistic equation with a free boundary, J. Dyn. Diﬀ. Equat., 30(2018), 1389-14 26

work page 2018
[1]

Berestycki, O

H. Berestycki, O. Diekmann, C. J. Nagelkerke, and P. A. Ze geling, Can a species keep pace with a shifting climate?, Bull. Math. Biol. 71 (2009), no. 2, 399-429

work page 2009
[2]

Berestycki and Jian Fang, Forced waves of the Fisher-K PP equation in a shifting environment, J

H. Berestycki and Jian Fang, Forced waves of the Fisher-K PP equation in a shifting environment, J. Diﬀerential Equations 264 (2018), no. 3, 2157-2183

work page 2018
[3]

Bunting, Y

G. Bunting, Y. Du, and K. Krakowski, Spreading speed revi sited: Analysis of a free boundary model, Networks and Heterogeneous Media, 7(2012), 583-603

work page 2012
[4]

Du and Z

Y. Du and Z. Lin, Spreading-vanishing dichotomy in the di ﬀusive logistic model with a free boundary, SIAM J. Math. Anal., 42(2010) 377-405. Erratum: SIAM J. Math . Anal., 45(2013), 1995-1996

work page 2010
[5]

Du and B

Y. Du and B. Lou, Spreading and vanishing in nonlinear diﬀ usion problems with free boundaries, J. Eur. Math. Soc., 17(2015) 2673-2724

work page 2015
[6]

Du and L

Y. Du and L. Ma, Logistic type equations on RN by a squeezing method involving boundary blow-up solutions, J. London Math. Soc., 64(2001), 107-124

work page 2001
[7]

Y. Du, H. Matsuzawa, and M. Zhou, Sharp estimate of the spr eading speed determined by nonlinear free boundary problems, SIAM J. Math. Anal., 46(2014), 375-396

work page 2014

Show all 20 references

[8]

Y. Du, H. Matsuzawa, and M. Zhou, Spreading speed and proﬁ le for nonlinear Stefan problems in high space dimensions, J. Math. Pures Appl., 103(2015), 741-787

2015
[10]

J. Fang, Y. Lou, and J. Wu, Can pathogen spread keep pace w ith its host invasion? SIAM J. Appl. Math. 76 (2016), no. 4, 1633-1657

2016
[11]

Fey and C

S.B. Fey and C. M. Herren. Temperature-mediated biotic interactions inﬂuence enemy release of nonnative species in warming environments, Ecology, 95(2014), no. 8, 2246-2256

2014
[12]

Gan and P

W. Gan and P. Zhou, A revisit to the diﬀusive logistic mod el with free boundary condition, Discrete Contin. Dyn. Syst. B, 21 (2016), 837-847

2016
[13]

Hu and X

H. Hu and X. Zou, Existence of an extinction wave in the Fi sher equation with a shifting habitat, Proc. Amer. Math. Soc. 145 (2017), no. 11, 4763-4771

2017
[14]

Lei and Y

C. Lei and Y. Du, Asymptotic proﬁle of the solution to a fr ee boundary problem arising in a shifting climate model, Discrete Contin. Dyn. Syst. B, 22 (2017), 895 -911

2017
[15]

M. A. Lewis, N. G. Marculis, and Z. Shen, Integrodiﬀeren ce equations in the presence of climate change: persistence criterion, travelling waves and inside dynami cs, J. Math. Biol. 77 (2018), no. 6-7, 1649-1687

2018
[16]

B. Li, S. Bewick, M. R. Barnard, and W. F. Fagan, Persiste nce and spreading speeds of integro-diﬀerence equations with an expanding or contracting habitat, Bull. M ath. Biol. 78 (2016), no. 7, 1337-1379

2016
[17]

B. Li, S. Bewick, J. Shang, and W. F. Fagan, Persistence a nd spread of species with a shifting habitat edge, SIAM J. Appl. Math. 5 (2014), 1397-1417; Erratum: 75 (2 015), 2379-2380

2014
[18]

W.T. Li, J.B. Wang, and X.Q. Zhao, Spatial dynamics of a n onlocal dispersal population model in a shifting environment, J. Nonlinear Sci. 28 (2018), no. 4, 11 89-1219

2018
[19]

Wang and X.Q

J.B. Wang and X.Q. Zhao, Uniqueness and global stabilit y of forced waves in a shifting environment, Proc. Amer. Math. Soc. 147 (2019), no. 4, 1467-1481

2019
[20]

Wang, Existence and uniqueness of solutions of free b oundary problems in heterogeneous environ- ments, Discrete Contin

M. Wang, Existence and uniqueness of solutions of free b oundary problems in heterogeneous environ- ments, Discrete Contin. Dyn. Syst. B, 24(2019), 415-421

2019

Pith tools

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

[1] [9]

Y. Du, L. Wei, and L. Zhou, Spreading in a shifting environ ment modeled by the diﬀusive logistic equation with a free boundary, J. Dyn. Diﬀ. Equat., 30(2018), 1389-14 26

work page 2018

[2] [1]

Berestycki, O

H. Berestycki, O. Diekmann, C. J. Nagelkerke, and P. A. Ze geling, Can a species keep pace with a shifting climate?, Bull. Math. Biol. 71 (2009), no. 2, 399-429

work page 2009

[3] [2]

Berestycki and Jian Fang, Forced waves of the Fisher-K PP equation in a shifting environment, J

H. Berestycki and Jian Fang, Forced waves of the Fisher-K PP equation in a shifting environment, J. Diﬀerential Equations 264 (2018), no. 3, 2157-2183

work page 2018

[4] [3]

Bunting, Y

G. Bunting, Y. Du, and K. Krakowski, Spreading speed revi sited: Analysis of a free boundary model, Networks and Heterogeneous Media, 7(2012), 583-603

work page 2012

[5] [4]

Du and Z

Y. Du and Z. Lin, Spreading-vanishing dichotomy in the di ﬀusive logistic model with a free boundary, SIAM J. Math. Anal., 42(2010) 377-405. Erratum: SIAM J. Math . Anal., 45(2013), 1995-1996

work page 2010

[6] [5]

Du and B

Y. Du and B. Lou, Spreading and vanishing in nonlinear diﬀ usion problems with free boundaries, J. Eur. Math. Soc., 17(2015) 2673-2724

work page 2015

[7] [6]

Du and L

Y. Du and L. Ma, Logistic type equations on RN by a squeezing method involving boundary blow-up solutions, J. London Math. Soc., 64(2001), 107-124

work page 2001

[8] [7]

Y. Du, H. Matsuzawa, and M. Zhou, Sharp estimate of the spr eading speed determined by nonlinear free boundary problems, SIAM J. Math. Anal., 46(2014), 375-396

work page 2014