Keiichi Umetsu, Michele Pizzardo, Antonaldo Diaferio, Margaret J. Geller

We present CLUMI+, a self-consistent, multi-probe methodology for reconstructing the mass distribution in and around galaxy clusters by combining gravitational lensing and dynamical observations. Building on the joint-likelihood framework of Umetsu (2013), CLUMI+ integrates weak-lensing shear and magnification data with projected escape velocity measurements in the cluster infall region, yielding tighter constraints on the gravitational potential without relying on equilibrium assumptions. The mass distribution is modeled using a flexible, piecewise-defined convergence profile that characterizes the azimuthally averaged surface mass density within the lensing field, transitioning to a projected power-law form at larger radii where phase-space constraints complement lensing. Additional strong-lensing constraints are incorporated via central aperture-mass measurements, enabling full-scale mass reconstruction from the cluster core to the outskirts. We validate CLUMI+ using synthetic weak-lensing and phase-space data for a massive cluster from the IllustrisTNG simulations, demonstrating unbiased recovery of projected and three-dimensional mass profiles and achieving 10%--30% improvement in precision at large radii. As a case study, we apply CLUMI+ to Abell 2261, combining Subaru and Hubble Space Telescope weak+strong lensing data with spectroscopic measurements from the Hectospec Cluster Survey. This analysis demonstrates the power of multi-probe, equilibrium-free modeling for robust cluster mass reconstruction.