Numerical climate models are complex and combine a large number of physical processes. They are key tools in quantifying the relative contribution of potential anthropogenic causes (e.g., the current increase in greenhouse gases) on high-impact atmospheric variables like heavy rainfall or temperatures. These so- called climate extreme event attribution problems are particularly challenging in a multivariate context, that is, when the atmospheric variables are measured on a possibly high-dimensional grid. In addition, global climate models like any in sillico numerical experiments are affected by different types of bias. In this talk, I will discuss about how to combine to two statistical theories to assess causality in the context of extreme event attribution. In addition, the question of uncertainties quantification that remains a challenge in any climate attribution analysis will be explored from various directions. In particular, a simple model bias correction step for records will described in details. To illustrate our approach, we infer emergence times in precipitation from the CMIP5 and CMIP6 archives. Joint work with Anna Kiriliouk, Paula Gonzalez Soulivanh Thao and Julien Worms Biblio: - Naveau P. and S. Thao Multi-model errors and emergence times in climate attribution studies, journal of climate, (2022) - Worms J. and P. Naveau. Record events attribution in climate studies. Envirmetrics (2022, in press). ⟨hal-02938596⟩ -Kiriliouk, A., and P. Naveau, 2020: Climate extreme event attribution using multivariate peaks- over-thresholds modeling and counterfactual theory. Ann. Appl. Stat., 14 (3), 1342–1358