Lynch, J. T., Schneider, M. T. Y., Perriman, D. M., Scarvell, J. M., Pickering, M. R., Asikuzzaman, M., Galvin, C. R., Besier, T. F., & Smith, P. N. (2019). Statistical shape modelling reveals large and distinct subchondral bony differences in osteoarthritic knees. Journal of Biomechanics, 93, 177–184. https://doi.org/10.1016/j.jbiomech.2019.07.003
-or you can try to limit the area of alignment to some landmark more in the subchondral area, for instance you can do that manually putting a few landmarks on the subchondral area of all meshes and do a landmark alignment, or you can just try to make sure all the landmark correspondence that you're making in the ICP alignment phase are not somewhere on the length of the bone.
Best Wishes,
saeed