A new implicit-explicit local differential transform method (IELDTM) is derived here for time integration of the nonlinear (2 + 1)-dimensional advection-diffusion-reaction (ADR) equations. The IELDTM is adaptively constructed as a stability preserved and high order time integrator for spatially discretized ADR equations. For spatial discretization of the model equation, the Chebyshev spectral collocation method (ChCM) is utilized. A robust stability analysis and global error analysis of the IELDTM are presented with respect to the direction parameter theta. With the help of the global error analysis, adaptivity equations are derived to minimize the computational costs of the algorithms. The produced method is shown to eliminate the accuracy disadvantage of the classical theta-method and the stability disadvantages of the existing differential transform-based methods. Two examples of the Burgers equation in one and two space dimensions and the Chapman oxygen-ozone ADR model are solved via the ChCM-IELDTM hybridization. The present time integrator is proven to provide more efficient numerical characteristics than the various multi-step and multi-stage time integration methods. The IELDTM is extensively compared with the widely used MATLAB solvers, ode45 and ode15s. The adaptive IELDTM has been proven to integrate the stiff Chapman ADR equations with optimum costs over relatively long-time intervals.