We present the results of downscaling CMIP6 global climate projections to local scales for the Mediterranean and Italian regions, aiming to produce high-resolution climate information for the assessment of climate change signals, with a focus on precipitation extreme events. We performed hindcast (i.e. ERA5-driven) and historical simulations (driven by the MPI-ESM1-2-HR model) to simulate the present (1980–2014) and future (2015–2100) climate under three different emission scenarios (SSP1-2.6, SSP2-4.5, SSP5-8.5).</p> <p>For each experiment, a double nesting approach is adopted to dynamically downscale global data to the regional domain of interest, firstly over the Europe (EURO) CORDEX domain, at a spatial resolution of 15 km, and then further refined (second nesting) over Italy and north-western Mediterranean, at a resolution of 5 km, i.e. in the so-called gray-zone (5–10 km), close to the convection permitting (CP) limit. Besides validating the experimental protocol, this work potentially questions the need for climate simulations to always resort to deep-convection parameterizations when spatial refinement is increased up to the limit of the CP scale, yet convective processes are still not explicitly resolved. Analyses of the most relevant Essential Climate Variables (ECVs) are presented, with a focus on the spatial distribution of precipitation, its probability density function, and the statistics of extreme events, for both current climate and far-end scenarios. By the end of the century for all the scenarios and seasons there is a projected general warming along with an intensification of the hydrological cycle over most of the continental EU and mean precipitation reduction over the Mediterranean region accompanied, over Italian Peninsula, by a strong increase in the intensity of extreme precipitation events, particularly relevant for the SSP5-8.5 scenario during autumn.