pith. sign in

arxiv: 1907.01983 · v1 · pith:4PJRRVZHnew · submitted 2019-07-03 · 🌊 nlin.PS · q-bio.PE

Model of pattern formation in marsh ecosystems with nonlocal interactions

Sofya Zaytseva , Junping Shi , Leah B Shaw This is my paper

Pith reviewed 2026-05-25 09:31 UTC · model grok-4.3

classification 🌊 nlin.PS q-bio.PE
keywords pattern formationnonlocal interactionsMexican-hat kernelreaction-diffusion equationsmarsh ecosystemsscale-dependent feedbackbiharmonic approximationSpartina alterniflora
0
0 comments X

The pith

Spatial patterns emerge in marsh models when the Mexican-hat kernel width and amplitude exceed thresholds for scale-dependent feedback.

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

The paper builds a reaction-diffusion model of grass and sediment volume in tidal marshes that includes a nonlocal interaction term. A Mexican-hat kernel encodes short-range positive feedback, in which grass traps sediment and promotes growth, together with long-range negative feedback, in which diverted flow erodes troughs farther away. A steady-state biharmonic approximation of the system yields explicit conditions under which uniform states lose stability and spatial patterns appear. These patterns are proposed to account for the irregular shoreline shapes seen in real Spartina marshes. The central result is that pattern formation occurs only for restricted ranges of the kernel's width and amplitude.

Core claim

We propose a mathematical framework to model grass-sediment dynamics as a system of reaction-diffusion equations with an additional nonlocal term quantifying the short-range positive and long-range negative grass-sediment interactions. We use a Mexican-hat kernel function to model this scale-dependent feedback. We perform a steady state biharmonic approximation of our system and derive conditions for the emergence of spatial patterns, corresponding to a spatially varying marsh shoreline. We find that the emergence of such patterns depends on the spatial scale and strength of the scale-dependent feedback, specified by the width and amplitude of the Mexican-hat kernel function.

What carries the argument

Mexican-hat kernel function placed inside the nonlocal term of the reaction-diffusion equations to represent short-range positive and long-range negative grass-sediment interactions.

If this is right

  • Patterns appear only when the kernel width and amplitude satisfy derived inequalities from the biharmonic analysis.
  • The wavelength of emerging patterns is set by the kernel width parameter.
  • Uniform marsh states remain stable for kernels whose amplitude is too small or whose width is outside the critical range.
  • The model links a concrete mathematical form of scale-dependent feedback to the irregular shoreline geometry observed in the field.

Where Pith is reading between the lines

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

  • If the kernel form is realistic, then controlled removal or addition of grass patches at specific distances could be used to test whether erosion troughs form at the predicted length scale.
  • The same nonlocal structure might be applied to other sediment-vegetation systems whose positive and negative feedbacks operate at different spatial ranges.
  • Time-dependent numerical solutions of the original integro-differential system could reveal whether patterns coarsen or stabilize after the initial instability.

Load-bearing premise

The short-range positive and long-range negative grass-sediment interactions can be accurately represented by a single Mexican-hat kernel function inside the reaction-diffusion system.

What would settle it

Direct measurement of sediment accretion and erosion distances around Spartina patches that yields an interaction kernel whose shape deviates substantially from a Mexican-hat profile, or observation of uniform shorelines when the measured kernel width and amplitude lie inside the model's predicted patterning regime.

Figures

Figures reproduced from arXiv: 1907.01983 by Junping Shi, Leah B Shaw, Sofya Zaytseva.

Figure 1
Figure 1. Figure 1: a) Self-organization on the marsh edge in the York River, a tributary of the Chesapeake [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Diagram of grass-sediment interactions adapted from (Bertness 1984). [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Illustration of a) the cross-section of marsh edge used to model the marsh dynamics [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Schematics representations of parameter regimes for the positive coexistence steady state [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Bifurcation diagrams plotted using MatCont (Dhooge et al. 2008) for Case I with a saddle [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Parameter space of Turing-like instability satisfying conditions (3.34) for various values [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Spatial patterns produced through simulations of the biharmonic system (3.39) (panels [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Panel a) displays the final steady states of grass after 1800 time units in the biharmonic [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: We numerically integrate the fourth order biharmonic system (3.39) for different kernel [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Functions a) F(S) and b) L(G) from (2.5) . a) c) d) b) [PITH_FULL_IMAGE:figures/full_fig_p024_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Spatial patterns produced through simulations of the biharmonic system (3.39) (panels [PITH_FULL_IMAGE:figures/full_fig_p024_11.png] view at source ↗
read the original abstract

Smooth cordgrass Spartina alterniflora is a grass species commonly found in tidal marshes. It is an ecosystem engineer, capable of modifying the structure of its surrounding environment through various feedbacks. The scale-dependent feedback between marsh grass and sediment volume is particularly of interest. Locally, the marsh vegetation attenuates hydrodynamic energy, enhancing sediment accretion and promoting further vegetation growth. In turn, the diverted water flow promotes the formation of erosion troughs over longer distances. This scale-dependent feedback may explain the characteristic spatially varying marsh shoreline, commonly observed in nature. We propose a mathematical framework to model grass-sediment dynamics as a system of reaction-diffusion equations with an additional nonlocal term quantifying the short-range positive and long-range negative grass-sediment interactions. We use a Mexican-hat kernel function to model this scale-dependent feedback. We perform a steady state biharmonic approximation of our system and derive conditions for the emergence of spatial patterns, corresponding to a spatially varying marsh shoreline. We find that the emergence of such patterns depends on the spatial scale and strength of the scale-dependent feedback, specified by the width and amplitude of the Mexican-hat kernel function.

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

0 major / 2 minor

Summary. The manuscript proposes a system of reaction-diffusion equations with a nonlocal Mexican-hat kernel to capture short-range positive and long-range negative grass-sediment interactions in tidal marsh ecosystems. A steady-state biharmonic approximation is applied to derive conditions under which spatial patterns emerge, with the onset depending on the width and amplitude of the kernel.

Significance. If the derivations are correct, the work supplies an explicit mathematical link between kernel parameters and pattern formation, offering a mechanistic explanation for observed marsh shoreline heterogeneity. The framework is internally consistent as a modeling exercise and generates falsifiable predictions tied to measurable spatial scales of feedback.

minor comments (2)
  1. The abstract states the biharmonic approximation is performed but does not indicate the order of the expansion or the regime of validity; adding one sentence on the truncation error would improve clarity for readers.
  2. Notation for the kernel function (width and amplitude parameters) should be introduced once in the model section and used consistently thereafter to avoid any ambiguity when the conditions are stated.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary and significance assessment of our manuscript, as well as the recommendation for minor revision. No specific major comments were listed in the report, so we have no points to address point-by-point. We are pleased that the internal consistency and falsifiable predictions were recognized.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper constructs a reaction-diffusion model augmented by a chosen Mexican-hat kernel for scale-dependent feedback, then applies a steady-state biharmonic approximation to obtain explicit conditions on pattern onset in terms of kernel width and amplitude. This is a direct mathematical consequence of the stated model equations and approximation; the kernel shape is introduced as an explicit modeling assumption rather than derived or fitted from the target pattern result. No self-definitional steps, fitted inputs renamed as predictions, load-bearing self-citations, or imported uniqueness theorems appear in the provided abstract or description. The derivation remains self-contained against the model's own equations.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on the assumption that a Mexican-hat kernel captures the biological and physical interactions; no free parameters are fitted in the abstract, but the kernel width and amplitude function as tunable scales.

free parameters (2)
  • Mexican-hat kernel width
    Sets the spatial scale separating positive and negative feedback; chosen to match observed marsh pattern wavelengths.
  • Mexican-hat kernel amplitude
    Sets the relative strength of local facilitation versus long-range inhibition.
axioms (2)
  • domain assumption The scale-dependent feedback between vegetation and sediment can be represented by a single nonlocal integral term with Mexican-hat shape.
    Invoked when the authors replace local interactions with the nonlocal kernel in the reaction-diffusion system.
  • domain assumption The steady-state biharmonic approximation preserves the pattern-forming instability of the original system.
    Used to derive the explicit conditions on kernel parameters.

pith-pipeline@v0.9.0 · 5729 in / 1408 out tokens · 31011 ms · 2026-05-25T09:31:30.001285+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

13 extracted references · 13 canonical work pages · 1 internal anchor

  1. [1]

    & Ayres, D

    Adams, J., Grobler, A., Rowe, C., Riddin, T., Bornman, T. & Ayres, D. (2012), ‘Plant traits and spread of the invasive salt marsh grass, spartina alterniflora loisel., in the great brak estuary, 25 south africa’, African Journal of Marine Science 34(3), 313–322. Altieri, A. H., Silliman, B. R. & Bertness, M. D. (2007), ‘Hierarchical organization via a faci...

    work page 2012

  2. [2]

    Bates, P. W. & Ren, X. (1996), ‘Transition layer solutions of a higher order equation in an infinite tube’, Communications in Partial Differential Equations 21(1-2), 109–145. Bates, P. W. & Ren, X. (1997), ‘Heteroclinic orbits for a higher order phase transition problem’, European Journal of Applied Mathematics 8(02), 149–163. Bayliss, A. & Volpert, V. (201...

    work page 1996

  3. [3]

    & Ridolfi, L

    Borgogno, F., D’Odorico, P., Laio, F. & Ridolfi, L. (2009), ‘Mathematical models of vegetation pattern formation in ecohydrology’, Reviews of Geophysics 47(1). Bouma, T., Friedrichs, M., Van Wesenbeeck, B., Temmerman, S., Graf, G. & Herman, P. (2009), ‘Density-dependent linkage of scale-dependent feedbacks: A flume study on the intertidal macro- phyte spart...

    work page 2009

  4. [4]

    & Gai, C

    Chen, Y., Kolokolnikov, T., Tzou, J. & Gai, C. (2015), ‘Patterned vegetation, tipping points, and the rate of climate change’, European Journal of Applied Mathematics 26(6), 945–958. Couteron, P. & Lejeune, O. (2001), ‘Periodic spotted patterns in semi-arid vegetation explained by a propagation-inhibition model’, Journal of Ecology 89(4), 616–628. Dakos, ...

    work page 2015

  5. [5]

    Green, J. B. & Sharpe, J. (2015), ‘Positional information and reaction-diffusion: two big ideas in developmental biology combine’, Development 142(7), 1203–1211. Halpern, B. S., Silliman, B. R., Olden, J. D., Bruno, J. P. & Bertness, M. D. (2007), ‘Incorporat- ing positive interactions in aquatic restoration and conservation’, Frontiers in Ecology and the ...

    work page 2015

  6. [6]

    J., Herman, P

    Liu, Q.-X., Weerman, E. J., Herman, P. M., Olff, H. & van de Koppel, J. (2012), ‘Alternative mechanisms alter the emergent properties of self-organization in mussel beds’, Proceedings of the Royal Society of London B: Biological Sciences p. rspb20120157. Madzvamuse, A., Ndakwo, H. S. & Barreira, R. (2015), ‘Cross-diffusion-driven instability for reaction-di...

    work page 2012

  7. [7]

    From scale-dependent feedbacks to long-range competition alone: a short review on pattern-forming mechanisms in arid ecosystems

    Mariotti, G. & Fagherazzi, S. (2010), ‘A numerical model for the coupled long-term evolution of salt marshes and tidal flats’, Journal of Geophysical Research: Earth Surface 115(F1). Mart´ ınez-Garc´ ıa, R., Calabrese, J. M., Hern´ andez-Garc´ ıa, E. & L´ opez, C. (2013), ‘Vegetation pattern formation in semiarid systems without facilitative mechanisms’, G...

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  8. [8]

    Murray, J. D. (2001), Mathematical Biology. II Spatial Models and Biomedical Applications {Interdisciplinary Applied Mathematics V. 18}, Springer-Verlag New York Incorporated. Nakamasu, A., Takahashi, G., Kanbe, A. & Kondo, S. (2009), ‘Interactions between zebrafish pigment cells responsible for the generation of turing patterns’, Proceedings of the Nation...

    work page 2001

  9. [9]

    Schile, Lisa M., e. a. (2014), ‘”modeling tidal marsh distribution with sea-level rise: Evaluating the role of vegetation, sediment, and upland habitat in marsh resiliency.”’, PloS one 9.2 . Schwarz, C., Bouma, T., Zhang, L., Temmerman, S., Ysebaert, T. & Herman, P. (2015), ‘Inter- actions between plant traits and sediment characteristics influencing speci...

    work page 2014

  10. [10]

    Vandenbruwaene, W., Temmerman, S., Bouma, T., Klaassen, P., De Vries, M., Callaghan, D., Van Steeg, P., Dekker, F., Van Duren, L., Martini, E. et al. (2011), ‘Flow interaction with dynamic vegetation patches: Implications for biogeomorphic evolution of a tidal landscape’, Journal of Geophysical Research: Earth Surface 116(F1). Watt, C., Garbary, D. J. & L...

    work page 2011

  11. [11]

    (1998), ‘Spatial heterogeneity in three species, plant-parasite-hyperparasite, systems’, Philosophical Transactions of the Royal Society B: Biological Sciences 353(1368),

    White, K. (1998), ‘Spatial heterogeneity in three species, plant-parasite-hyperparasite, systems’, Philosophical Transactions of the Royal Society B: Biological Sciences 353(1368),

    work page 1998

  12. [12]

    31 Yang, W., Wang, Q., Pan, X., Li, B. et al. (2014), ‘Estimation of the probability of long-distance dispersal: Stratified diffusion of spartina alterniflora in the yangtze river estuary’, American Journal of Plant Sciences 5(24),

    work page 2014

  13. [13]

    Ysebaert, T., Yang, S.-L., Zhang, L., He, Q., Bouma, T. J. & Herman, P. M. (2011), ‘Wave attenuation by two contrasting ecosystem engineering salt marsh macrophytes in the intertidal pioneer zone’, Wetlands 31(6), 1043–1054. Zaytseva, S., Shaw, L., Lipcius, R., Shi, J. & Kirwan, M. (2018), ‘Pattern formation in marsh ecosys- tems modeled through the inter...

    work page 2011