Analyzing and forecasting the Indian Summer Monsoon Rainfall (ISMR) is vital for South Asia’s socio-economic stability. Using 35 climate models from the latest generation of the Coupled Model Intercomparison Project (CMIP6) to simulate and project ISMR, we integrated statistical methods, such as Taylor diagrams, comprehensive rating indicators, and interannual variability scores, to compare performance differences between various models and analyze influencing mechanisms. The results show that the majority of models effectively simulate the climatology of the ISMR. However, they exhibit limitations in accurately capturing its interannual variability. Importantly, we observed no significant correlation between a model’s ability to simulate ISMR’s general climatology and its accuracy in representing annual variability. After a comprehensive assessment, models, like BCC-ESM1, EC-Earth3-Veg, GFDL-CM4, INM-CM5-0, and SAM0-UNICON were identified as part of the prime model mean ensemble (pMME), demonstrating superior performance in spatiotemporal simulations. The pMME can accurately simulate the sea surface temperature changes in the North Indian Ocean and the atmospheric circulation characteristics of South Asia. This accuracy is pivotal for CMIP6’s prime models to precisely simulate ISMR climatic variations. CMIP6 projections suggest that, by the end of the 21st century, ISMR will increase under low, medium, and high emission scenarios, with a significant rise in rainfall under the high emission scenario, especially in the western and northern parts of India. Among the pMME, the projected increase in rainfall across India is more moderate, with an estimated increase of 30%. The findings of this study suggest that selecting the best models for regional climate downscaling research will project regional climate changes more accurately. This provides valuable recommendations for model improvements in the Indian region.