In this study, an ensemble Bayesian inference-based copula approach (EBICA) has been proposed to evaluate the joint risk of extreme precipitation indices under climate change. EBICA incorporates GCMs, quantile delta mappingbased bias correction method, Monte Carlo Markov Chain (MCMC)-based Bayesian inference method, and copula functions into an integration, where 17 GCMs and 3 copulas are explored EBICA can reflect the inherent uncertainties in GCM and copula selections, as well as deal with the uncertainty caused by copula parameters when modelling the dependent structures.Ensemble results of GCMs and copulas are also analysed to promote the projection accuracy of dynamic features of extreme precipitation (characterized by precipitation intensity, very heavy precipitation amount and frequency, and maximum 1 or 5-day precipitation amount). EBICA is then applied to Fujian Province to investigate climate change impacts on extreme precipitation during 2021-2099. Major findings can be summarized as (a) GCM is the primary source of uncertainties that impacts extreme precipitation projection, and the ensemble of GCMs suggests that climate change would exacerbate future precipitation intensity, amount, and frequency over Fujian Province; (b) mountainous area in Fujian is more easily encountered with very heavy precipitation amount and coastal area is more easily faced with maximum 1-day or consecutive 5-day precipitation under climate change; (c) Fujian Province would averagely encounter with shorter co-occurrence period of [SDII, R95p] and [R20mm, R95p], corresponding to higher risk level in very heavy precipitation amount and frequency, especially in the north part. Generally, these findings could help lay the foundation for forecasting and warning of extreme precipitation events and reducing associated risk losses.