The slow deformation of terrestrial orbits in the medium range, subject to lunisolar resonances, is well approximated by a family of Hamiltonian flow with 2.5 degree-of-freedom. The action variables of the system may experience chaotic variations and large drift that we may quantify. Using variational chaos indicators, we compute high-resolution portraits of the action space. Such refined meshes allow to reveal the existence of tori and structures filling chaotic regions. Our elaborate computations allow us to isolate precise initial conditions near specific zones of interest and study their asymptotic behaviour in time. Borrowing classical techniques of phase-space visualization, we highlight how the drift is mediated by the complement of the numerically detected KAM tori.