The present study examines severe weather proxies including CAPE, convective inhibition (CIN), 0-1 km and 0-6 km wind shear in 12 GCMs from the Coupled Model Intercomparison Project phase 6 (CMIP6) and evaluate their skill relative to two independent reanalysis products – ERA5 and MERRA2. We also analyze composite parameters such as significant tornado parameter (STP) and supercell composite parameter (SCP) to evaluate the co-variability of component parameters. The analysis reveals that while most models show a reasonable skill in capturing the annual and seasonal distribution of observed parameters, significant bias exists within the ensemble. For instance, several models tend to overestimate CAPE values considerably, specifically over the oceanic regions. This characteristic is very likely associated with the errors arising from the boundary layer and convective parameterization schemes employed in the respective models. The mixed biases among the ensemble in component parameters (CAPE, CIN, and wind shear) consequently impact the skill of composite parameters. Additionally, we compare the distributions of vertical profiles of temperature, humidity, and wind in the models with rawinsonde profiles from Integrated Global Radiosonde Archive version 2 (IGRA2) during convective and non-convective periods to identify the biases in basic model variables. The findings from this study will offer valuable information to choose the right set of models for downscaling applications, and for investigating the future changes in severe storm environments over the US in general.

