Abstract:We study steady states in a reaction-diffusion system for a benthic bacteria-nutrient model in a marine sediment over 1D and 2D domains by using Landau reductions and numerical path following methods. We point out how the system reacts to changes of the strength of food supply and ingestion. We find that the system has a stable homogeneous steady state for relatively large rates of food supply and ingestion, while this state becomes unstable if one of these rates decreases and Turing patterns such as hexagons and stripes start to exist. One of the main results of the present work is a global bifurcation diagram for solutions over a bounded 2D domain. This bifurcation diagram includes branches of stripes, hexagons, and mixed modes. Furthermore, we find a number of snaking branches of stationary states, which are spatial connections between homogeneous states and hexagons, homogeneous states and stripes as well as stripes and hexagons in parameter ranges, where both corresponding states are stable. The system under consideration originally contains some spatially varying coefficients and with these exhibits layerings of patterns. The existence of spatial connections between different steady states in bistable ranges shows that spatially varying patterns are not necessarily due to spatially varying coefficients. The present work gives another example, where these effects arise and shows how the analytical and numerical observations can be used to detect signs that a marine bacteria population is in danger to die out or on its way to recovery, respectively. We find a type of hexagon patches on a homogeneous background, which seems to be new discovery. We show the first numerically calculated solution-branch, which connects two different types of hexagons in parameter space. We check numerically for bounded domains whether the stability changes for hexagons and stripes, which are extended homogeneously into the third dimension. We find that stripes and one type of hexagons have the same stable range over bounded 2D and 3D domains. This does not hold for the other type of hexagons. Their stable range is shorter for the bounded 3D domain, which we used here. We find a snaking branch, which bifurcates when the hexagonal prisms loose their stability. Solutions on this branch connects spatially between hexagonal prisms and a genuine 3D pattern (balls).

Numerical results for snaking of patterns over patterns in some 2D Selkov-Schnakenberg Reaction-Diffusion systems

Pattern formation with pde2path -- a tutorial

Homoclinic snaking of contact defects in reaction-diffusion equations

Unified framework for localized patterns in reaction–diffusion systems; the Gray–Scott and Gierer–Meinhardt cases

Pattern analysis in a benthic bacteria-nutrient system

Dissecting the snake: the transition from localized patterns to isolated spikes in pattern formation systems

Effects of cross-diffusion on Turing patterns in a reaction-diffusion Schnakenberg model

Spatial pattern formation in reaction–diffusion models: a computational approach

Time-dependent localized patterns in a predator-prey model

Stability of patterns on thick curved surfaces

Pattern formation of bulk-surface reaction-diffusion systems in a ball

Bifurcation Analysis of Reaction Diffusion Systems on Arbitrary Surfaces

Spot Patterns in the 2-D Schnakenberg Model with Localized Heterogeneities

Localized and extended patterns in the cubic-quintic Swift-Hohenberg equation on a disk

Localized patterns and semi-strong interaction, a unifying framework for reaction–diffusion systems

Variational approximations to homoclinic snaking

Differential geometric bifurcation problems in pde2path -- algorithms and tutorial examples

Oscillatory Instabilities of a One-Spot Pattern in the Schnakenberg Reaction-Diffusion System in $3$-D Domains

Pattern formation in nonlocal Kondo model

Pattern formations and optimal packing