ABSTRACTPrecipitation changes dynamically in the Mediterranean region. Therefore, the projection of future precipitation and its historical distribution mechanism is essential for climate mitigation and adaptation. In this study, a stepwise clustered precipitation downscaling method (SCPD) was developed and adopted in the Mediterranean region to reveal the inherent variation rules and trends over the future 100 years under two SSP scenarios. A cutting and merging multivariate process is introduced to build a cluster tree for supporting further downscaling and projecting steps. The ensemble average from the global climate model (GCM) dataset is used for precipitation projections. The precipitation performance of SCPD, evaluated by R2, is fairly decent. The precipitation projections vary with the original rainfall patterns over the gauge stations. Dry places tend to become comparably drier in the future. Precipitation in the northern Mediterranean region shows a drier winter–spring and wetter summer–autumn. Opposite trends emerged in the southern part, with increasing winter precipitation and decreasing summer rainfall. The rising carbon dioxide concentration will further intensify the decrease in rainfall. However, the centres of these two EOFs are not identical. The contributions of NAO (positive) and Niño 3.4 (negative) to PC1 are relatively high. Accordingly, the strongest positive correlation with PC2 is SCAND, as well as negative correlations with AO, NAO and EAWR. Positive anomaly precipitation is attributed to PC1, whereas PC2 is responsible for most of the negative variance precipitation.