In the current study, we introduce a model for the temporal and spatial interactions between prey and predator. The model incorporates a nonlinear refuge mechanism for prey, along with linear harvesting of prey and nonlinear harvesting for predators. Initially, we examined the well-posed nature of the model by analyzing the presence of all feasible equilibria and investigating the corresponding dynamics. Following that, we delve into the dynamics of the temporal model, focusing specifically on aspects such as uniform boundedness, permanence, and stability of viable equilibria. We demonstrate analytically that the proposed model experiences transcritical, saddle–node, Hopf, and Bogdanov–Takens bifurcations. It shows a variety of intricate dynamics involving Generalized Hopf and double Hopf. Then, discrete-time delay effects arising from the gestation of predator species have been incorporated into the temporal system. Hopf bifurcation for the delay parameter was detected in this investigation. Subsequently, we established conditions for self-diffusion instability and Turing instability in the spatiotemporal model, both with and without delay, employing the homogeneous Neumann boundary condition. Moreover, the discussion of sensitivity analysis (PRCC) serves to illustrate how crucial parameters influence the dynamics of the system. In addition, we conduct numerical simulations aiming to corroborate and validate the analytical results obtained.