We study the long time asymptotic behavior for the Cauchy problem of the modified Camassa-Holm (mCH) equation in the solitonic regions mt + (m (u2 - u2 )) x x + & UKappa;ux = 0, m = u - uxx, u(x, 0) = u0(x) & ISIN; H4,2(R), where & UKappa; is a positive constant characterizing the effect of the linear dispersion. Our main technical tool is the representation of the Cauchy problem with an associated matrix Riemann-Hilbert (RH) problem and the consequent asymptotic analysis of this RH problem. Based on the spectral analysis of the Lax pair associated with the mCH equation and scattering matrix, the solution of the Cauchy problem is characterized via the solution of a RH problem in the new scale (y, t). Further using the partial differential generalization of the Deift-Zhou steepest descent method, we derive different long time asymptotic expansions of the solution u(y, t) in different space-time solitonic regions of xi = y/t. We divide the half-plane {(y, t) : -& INFIN; < y < & INFIN;, t > 0} into four asymptotic regions: The phase function theta(z) has no stationary phase point on the jump contour in the space-time solitonic regions xi & ISIN; (-& INFIN;, -1/4) & OR; (2, +& INFIN;), corresponding asymptotic approximations can be characterized with an N(& ULambda;)-solitons with diverse residual error order O(t-1+2 rho); The phase function theta(z) has four phase points and eight phase points on the jump contour in the space-time solitonic regions xi & ISIN; (0, 2) and xi & ISIN; (-1/4, 0), respectively. The corresponding asymptotic approximations can be characterized with an N(& ULambda;)-soliton as well as an interaction term between soliton solutions and the dispersion term with diverse residual error order O(t-3/4). Our results also confirm the soliton resolution conjecture and asymptotically stability of the N-soliton solutions for the mCH equation.(c) 2022 Elsevier Inc. All rights reserved.