We use 3D cloud-resolving model (CRM) simulations of two mesoscale convective systems at midlatitudes and a simple statistical ensemble method to diagnose the scale dependency of convective momentum transport (CMT) and CMT-related properties and evaluate a parameterization scheme for the convection-induced pressure gradient (CIPG) developed by Gregory et al. Gregory et al. relate CIPG to a constant coefficient multiplied by mass flux and vertical mean wind shear. CRM results show that mass fluxes and CMT exhibit strong scale dependency in temporal evolution and vertical structure. The upgradient–downgradient CMT characteristics for updrafts are generally similar between small and large grid spacings, which is consistent with previous understanding, but they can be different for downdrafts across wide-ranging grid spacings. For the small to medium grid spacings (4–64 km), Gregory et al. reproduce some aspects of CIPG scale dependency except for underestimating the variations of CIPG as grid spacing decreases. However, for large grid spacings (128–512 km), Gregory et al. might even less adequately parameterize CIPG because it omits the contribution from either the nonlinear-shear or the buoyancy forcings. Further diagnosis of CRM results suggests that inclusion of nonlinear-shear forcing in Gregory et al. is needed for the large grid spacings. For the small to median grid spacings, a modified Gregory et al. with the three-updraft approach help better capture the variations of CIPG as grid spacing decreases compared to the single updraft approach. Further, the optimal coefficients used in Gregory et al. seem insensitive to grid spacings, but they might be different for updrafts and downdrafts, for different MCS types, and for zonal and meridional components.