“…However, it is not sufficient to quantify concurrent extremes by conducting pairs of univariate analyses (i.e., only analyzing each variable separately) as by doing so could lead to under-or overestimation of risk if the variables of interest are respectively positively or negatively related to each other. There are relatively few climate studies considering extremes in a multivariate setting despite a large body of work in the statistical community having been dedicated to modeling multivariate and spatial extremes (Tawn, 1988(Tawn, , 1990Smith, 1990;Coles and Tawn, 1991;Tawn, 1996, 1997;Coles et al, 1999;Heffernan and Tawn, 2004;Cooley et al, 2006;Naveau et al, 2009;Davison et al, 2012;Wadsworth and Tawn, 2012a;Huser and Davison, 2014;Wadsworth and Tawn, 2018;Huang et al, 2019a;Wadsworth and Tawn, 2019;Huser and Wadsworth, 2019;Cooley et al, 2019;Beranger et al, 2019;Bopp et al, 2020). The existing methods for modeling multivariate (including spatial) extremes mostly focus on "component-wise extremes", in which extreme values for each component (e.g., climate variable) are first extracted separately and then combined to create a new extremal data vector.…”