This paper proposes an L1/LDG algorithm for a time–space fractional convection equation on a bounded rectangular domain
. Here, the temporal partial derivative is in the sense of Caputo with order
, and the spatial partial derivatives are of Riesz type with orders
in
‐ and
‐directions, respectively. L1 approximation is utilized to evaluate the temporal Caputo derivative, while the local discontinuous Galerkin (LDG) method is adopted to deal with the spatial Riesz derivative. The rigorous numerical stability analysis and error estimate are presented which are supported by illustrative numerical examples.