Pattern Dynamics Analysis of Host–Parasite Models with Aggregation Effect Based on Coupled Map Lattices

Host–parasitoid systems are an essential area of research in ecology and evolutionary biology due to their widespread occurrence in nature and significant impact on species evolution, population dynamics, and ecosystem stability. In such systems, the host is the organism being attacked by the parasitoid, while the parasitoid depends on the host to complete its life cycle. This paper investigates the effect of parasitoid aggregation attacks on a host in a host–parasitoid model with self-diffusion on two-dimensional coupled map lattices. We assume that in the simulation of biological populations on a plane, the interactions between individuals follow periodic boundary conditions. The primary objective is to analyze the complex dynamics of the host–parasitoid interaction model induced by an aggregation effect and diffusion in a discrete setting. Using the aggregation coefficient k as the bifurcating parameter and applying central manifold and normal form analysis, it has been shown that the system is capable of Neimark–Sacker and flip bifurcations even without diffusion. Furthermore, with the influence of diffusion, the system exhibits pure Turing instability, Neimark–Sacker–Turing instability, and Flip–Turing instability. The numerical simulation section explores the path from bifurcation to chaos through calculations of the maximum Lyapunov exponent and the construction of a bifurcation diagram. The interconversion between different Turing instabilities is simulated by adjusting the timestep and self-diffusion coefficient values, which is based on pattern dynamics in ecological modeling. This contributes to a deeper understanding of the dynamic behaviors driven by aggregation effects in the host–parasitoid model.