The nonlinear coupled reaction-diffusion (NCRD) systems are important in the formation of spatiotemporal patterns in many scientific and engineering fields, including physical and chemical processes, biology, electrochemical processes, fractals, viscoelastic materials, porous media, and many others. In this study, a mixed-type modal discontinuous Galerkin approach is developed for one-and two-dimensional NCRD systems, including linear, Gray-Scott, Brusselator, isothermal chemical, and Schnakenberg models to yield the spatiotemporal patterns. These models essentially represent a variety of complicated natural spatiotemporal patterns such as spots, spot replication, stripes, hexagons, and so on. In this approach, a mixed-type formulation is presented to address the second-order derivatives emerging in the diffusion terms. For spatial discretization, hierarchical modal basis functions premised on the orthogonal scaled Legendre polynomials are used. Moreover, a novel reaction term treatment is proposed for the NCRD systems, demonstrating an intrinsic feature of the new DG scheme and preventing erroneous solutions due to extremely nonlinear reaction terms. The proposed approach reduces the NCRD systems into a framework of ordinary differential equations in time, which are addressed by an explicit third-order TVD Runge-Kutta algorithm. The spatiotemporal patterns generated with the present approach are very comparable to those found in the literature. This approach can readily be expanded to