In this paper, two time fractional coupled systems of differential equations have been investigated
by using the Lie symmetry analysis. Vector field generators of these variants of Boussinesq equations are
extracted and corresponding reductions to the systems of fractional ordinary differential equations (FODEs)
with Erd´elyi-Kober fractional derivatives are obtained. Conservation laws of these fractional systems are
extracted. Finally, the invariant subspace method is applied to find the corresponding exact solutions.