Investigation on potential and limitations of ERA5 Reanalysis downscaled on Italy by a convection-permitting model

This study presents a comprehensive assessment of a dynamical downscaling of ERA5 Reanalysis recently performed over Italy through the COSMO-CLM model at a convection-permitting scale (0.02°) over the period 1989–2020. Results are analysed against several independent observational datasets and reanalysis products. The capability of the downscaling to realistically represent the climatology for 2 m temperature and precipitation is analysed over the whole peninsula and subdomains. Hourly precipitation patterns, orography effects, and urban climate dynamics are also investigated, highlighting the weaknesses and strengths of the convection-permitting model. In particular, gains in performances are achieved in mountainous areas where the climate characteristics are correctly represented, as are the hourly precipitation characteristics. Losses in performances occur in coastal and flat areas of the Italian peninsula, where the convection-permitting model performance does not seem to be satisfactory, as opposed to complex orographic areas. The adopted urban parameterisation is demonstrated to simulate heat detection for two Italian cities: Rome and Milan. Finally, a subset of extreme climate indicators is evaluated, finding: (i) a region-dependent response, (ii) a notable performance of the convection-permitting model over mountainous areas and (iii) discrepancies in the South, Central and Insular subdomains. Climate indicators detect extreme events at a detailed scale, becoming an important tool for turning climate data into information.


Introduction
In recent years, Italy, like other countries in the Mediterranean area, has been observing an increase in extreme weather events (i.e., heavy rainfall and increase in temperature), causing huge impacts (i.e., floods, droughts, heat waves) with consequences to assets and people (Amponsah et al. 2018;Gaume et al. 2009;Llasat et al. 2010;2013).These occurrences are specifically related to peculiar characteristics of the region, such as demographic vulnerability (i.e., high-density population), orographic features, and the complex terrain close to the coast (Miglietta and Davolio 2022).Moreover, their occurrence in multiple climatic impact-drivers perspectives is expected to increase (with high confidence) since global warming is expected to increase faster over this area than the global mean temperature change (IPCC 2021).Therefore, there is an increasing need for more reliable and detailed climate information on the Italian peninsula to improve the assessment of climate hazards and risk management.
Convection-permitting regional climate models (CP-RCMs) are increasingly considered a relevant step forward in studies for understanding past climate and climate change at local scales and extreme weather events that most impact society (Kendon et al. 2021), holding promising potentials (Prein et al. 2013).CP-RCMs have indeed highlighted added-value concerning the representation of precipitation, whether it concerns their spatial or temporal characteristics, particularly on the sub-daily time scale, especially in summer (Keller et al. 2016;Lind et al. 2016;Fowler et al. 2021).Several authors (Adinolfi et al. 2021;Ban et al. 2021;Berthou et al. 2020;Reder et al. 2020) demonstrated these improvements in complex orographic contexts, such as the Alpine area, also in the frame of different European projects and initiative as the European Climate Prediction system-Horizon 2020 project (https:// www.eucp-proje ct.eu/) and the CORDEX Flagship Pilot Studies (Coppola et al. 2020) on Convective phenomena at high resolution over Europe and the Mediterranean (COR-DEX FPS-CONV).The performances of CP-RCMs with "event-based" approaches for short-term statistics were testes in other studies (Coppola et al. 2020;Raffa et al. 2021b).Coppola et al. (2020) focused on three case studies indicating that large-scale driven events are captured in climate-mode even if the inter-model spread increases.Moreover, they argued the relevance of an ensemblebased approach to investigating high-impact convective processes.Raffa et al. (2021b) provided a reliable eventbased analysis over two summer precipitation events in August 2007 and 2010 based on a direct downscaling of ERA5 (one-way nesting) at km-scale and a two-step nested one.The one-way nesting recognised the timing and intensity of the rainfall events.At the same time, the two-step nested one completely fails in their detection, possibly due to the influence of the lateral boundary condition by a freely evolving (i.e., not nudged) intermediate simulation allowing internal variability to develop.
Due to the previous promising outcomes of CP-RCMs, Raffa et al. (2021a) have released VHR-REA_IT, a gridded hourly dataset over Italy for recent climate (January 1989 -December 2020) integrated into the CMCC Data Delivery System (DDS-http:// dds.cmcc. it accessed on 4 August 2021).The VHR-REA_IT dataset derives from the homonymous hindcast simulation that is the dynamical downscaling of ERA5 Reanalysis (Hersbach et al. 2020) at a convectionpermitting scale (0.02°, i.e., ≃ 2.2 km).The dynamical downscaling has been performed with the COSMO-CLM (Consortium for Small-Scale Modeling in Climate Mode) Regional Climate Model (Rockel et al. 2008), switching on a specific module (i.e., TERRA-URB, Wouters et al. 2016) for modelling urban areas.VHR-REA_IT exploits the added value of CP-RCMs for short-duration precipitation extremes, avoiding error-prone deep convection parameterisation schemes, which are sometimes responsible for misinterpreting precipitation patterns and trends (Kendon et al. 2021).It also brings with it the advantages of using urban parameterisation.This latter acts in synergy with the enhanced horizontal resolution, allowing for a better description of land-use dynamics related to the presence of cities (Gutowski et al. 2020), as well as increasing the relevance of activating specific urban canopy parameterisations (Wouters et al. 2016).Preliminary experiments (Garbero et al. 2021) have demonstrated the capability of the COSMO-CLM with TERRA-URB to reasonably detect the Urban Heat Island (UHI) effect and improve air temperature forecasts in different urban areas with specific morphological characteristics.Such capability provides relevant perspectives for enhancing urban climate modelling activities.
This paper presents an in-depth evaluation of the VHR-REA_IT, highlighting the weaknesses and strengths of this dynamical downscaling of ERA5 at the convection-permitting scale.Furthermore, it aims to explore the capability of the hindcast simulation to represent the climatology of average and extreme patterns of temperature and precipitation in a manner consistent with different independent observation datasets.The overall goal of the paper is to highlight some steps forward in comparison with state-of-the-art: (i) providing climate information with increased space (≃ 2.2 km) and time (1 h) resolutions; (ii) exploiting the potential of convection-permitting models and urban parameterisation to represent the past and current climate synergically; (iii) characterising on a detailed scale extreme events having the potential for significant negative impacts on lives, infrastructure or the environment (Fowler and Ali 2022) Specific analyses focus on the hourly precipitation features, orography effects, and urban dynamics.Moreover, extreme events are investigated considering a subset of indices for precipitation and temperature of the World Meteorological Organization (WMO) Expert Team on Climate Change Detection and Indices (ETCCDI).
The paper is organized as follows.First, Sect. 2 introduces the regional climate model with initial and boundary conditions, VHR-REA_IT setup details, the observations and reanalysis products assumed as reference for evaluation, and the evaluation method considered for analysing extremes.Section 3 presents and discusses the evaluation results for temperature and precipitation, as regard mean characteristics (i.e., spatial patterns, orography effects) and extremes, giving evidence of specific peculiarities in urban dynamics and hourly precipitation patterns.Section 4 draws conclusions.

Regional Climate Model COSMO-CLM and set-up
VHR-REA_IT has been performed by dynamically downscaling ERA5 Reanalysis at 0.02° (≃ 2.2 km) with the regional climate model COSMO-CLM (Rockel et al. 2008).It is a non-hydrostatic and limited-area model designed for dynamically downscaling simulations at different horizontal resolutions varying from the meso-β (horizontal scales in-between 20-200 km) to the meso-γ (horizontal scales in-between 2-20 km).It exploits finite difference methods to solve fully compressible governing equations of fluid dynamics on a structured grid using finite difference methods.
Horizontal advection is calculated using a fifth-order upwind scheme, and vertical advection is computed with an implicit Crank-Nicholson scheme (Förstner and Doms 2004).Time integration is performed with a split-explicit third-order Runge-Kutta discretisation (Baldauf et al. 2011).Cloud microphysics is represented with a single-moment scheme (Doms and Baldauf 2011) using five hydrometeors (cloud water, rain, ice crystals, snow, and graupel).The radiation scheme is based on a δ two-stream approach described by Ritter and Geleyn (Ritter and Geleyn 1992).
Turbulent fluxes within the planetary boundary layer are parameterised with a 1.5-order turbulent kinetic energy (TKE)-based scheme (Mellor and Yamada, 1982;Raschendorfer 2001).The standard parameterisation of convection used in the model is based on the Tiedtke scheme (Tiedtke 1989), a mass-flux closure approach used to parameterise modifications to the vertical structure of the atmosphere due to deep, midlevel and shallow convection.In convectionpermitting mode, the deep convection is explicitly resolved; therefore, only the shallow convection part of the scheme is active.
The representation of soil moisture is performed using a 7-layer soil model (TERRA_ML) with a formulation for water runoff dependent on orography (Schrodin and Heise 2001).COSMO-CLM also allows turning on an intrinsic representation of urban physics by employing TERRA-URB parameterisation (Wouters et al. 2016) with a bulk scheme to discern each grid cell between urban canopy and natural land cover and computes adjusted soil and water fluxes considering the urban environment features.
The VHR-REA_IT dynamical downscaling relies on the COSMO-DE model configuration developed by the Deutscher Wetterdienst for numerical weather prediction application.Such a configuration has been modified adequately in Raffa et al. (2021a) to account for the urban effects, switching on the TERRA-URB parameterisation, and to solve convection explicitly as in the frame of the CORDEX FPS-CONV initiative.
The simulation was performed on the supercomputer cluster GALILEO available at the Consorzio Interuniversitario del Nord-Est per il Calcolo Automatico (CINECA) using the version COSMO-CLM 5.00_clm9 associated with the interpolator version int2lm_2.00_clm4.It has been initialised and driven by the ERA5 reanalysis (Hersbach et al. 2020) of the European Centre for Medium-Range Weather Forecast (ECMWF) using a direct downscaling strategy as described in Raffa et al. 2021b.ERA5 currently represents the most plausible climate descriptor suitable for several applications, such as monitoring climate change, research, education, policy-making, and business.It has a global coverage with a spatial resolution of 0.28° (≈ 31 km) and provides outputs at an hourly scale from 1950 to the present.Recently, several dynamical downscaling experiments (Reder et al. 2022;Wang et al. 2021;Pilguj et al. 2022) relied on initial and boundary conditions from ERA5, being assumed as the reference reanalysis for the experiment design of the dynamical downscaling of the Coupled Model Intercomparison Project Phase 6 (CMIP6) global models generation over the EURO-CORDEX domain.
The computational domain (Lon = 5°W-20°E; Lat = 36°N-48°N; see Fig. 1a) consists of 585 x 730 grid points, 50 vertical levels and 7 soil levels.The analysis domain is obtained by discarding the relaxation zone from the computational domain (25 points for each border).The covered period is 1989-2020, assuming the year 1988 as spin-up.Its time step for integration based on a third-order Runge-Kutta scheme is 20 s.
The external datasets adopted in VHR-REA_IT are GLC2000 (Bartholomé and Belward 2005) for land use, GLOBE for surface elevation, and FAO Digital Soil Map of the World for soil types.Moreover, additional external datasets are the impervious surface area obtained from the European Environment Agency (Maucha et al. 2010) and the anthropogenic heat flux for the year 2006 derived by Flanner (2009), as TERRA-URB is activated.

Observations and reanalysis products
As the resolution of climate models increases, the evaluation of CP-RCM simulations becomes challenging due to the lack of high-quality and high-resolution observational gridded datasets (Piazza et al. 2019).Moreover, despite being all based on meteorological stations, a large discrepancy between the different observed gridded datasets has been highlighted (Zittis et al. 2017).This is especially true in areas with low population density (i.e., islands) and mountainous regions where weather stations are sparsely distributed and precipitation undercatches snowfall in winter.
In the previous study (Raffa et al. 2021a), the VHR-REA_IT dataset was validated against the observational dataset E-OBS (Cornes et al. 2018).The main findings of this evaluation are also reported in the following figures and discussion of the results in order to depict a comprehensive evaluation of the hindcast simulation.In this work, VHR-REA_IT has been evaluated against different independent observational datasets available for Italy.Specifically: • E-OBS (Cornes et al. 2018;Haylock et al. 2008) In addition, in order to overcome the discrepancies among different observational gridded datasets, VHR-REA_IT is also compared with reanalysis products as ERA5 (the driven reanalysis of VHR-REA_IT as detailed in §2.1) and UERRA.Specifically:

Methods for the evaluation
This study investigates the performances of VHR-REA_IT for temperature and total precipitation.To this aim, independent observational datasets such as SCIA-ISPRA, E-OBS and GRIPHO (see §2.2) have been assumed as references.SCIA-ISPRA and E-OBS are used as a baseline for temperature and precipitation on the daily scale across the period 1989-2020.For SCIA-ISPRA, the daily mean temperature is the average value between daily minima and maxima values (Allen et al. 1998).The daily mean temperature thus obtained for SCIA-ISPRA only provides a rough estimate, as the variable of daily average temperature is used directly for other datasets.GRIPHO is adopted to evaluate precipitation at the hourly scale for the entire available period (2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016).
For both variables, the evaluation is first conducted by comparing observations and simulation in terms of spatial distributions through maps; the qualitative comparison is then turned into quantitative information using specific metrics (i.e., bias and variability) computed at annual and seasonal scales to corroborate the analyses.Moreover, a comparison with reanalysis products such as ERA5 (the driven reanalysis of VHR-REA_IT) and UERRA (see §2.2) have been performed thorough maps and metrics able to assess the added value.In particular, the DAV metric (Soares and Cardoso 2018) provides an objective and normalized measure of the added value in terms of potential gain in the performance of climate models due to the usage of a higher resolution, comparing higher-and coarser-resolution simulation probability density function (PDFs) to the observational PDF (Cardoso and Soares 2022;Careto et al. 2022;Tölle and Churiulin 2021).DAV accounts for the difference in skill scores between high resolution (hr) and low resolution (lr), assuming the observations (obs) as reference:(1) where S hr and S lr are the Perkins skill scores for high and low resolution respectively; n is the number of bins considered to obtain the PDF; Z hr , Z lr and Z obs are the frequencies of values in each bin for high resolution, low resolution, and observations respectively.Specifically, the Perkins skill score is a quantitative measure of how well each simulation resembles the observed probability density functions by measuring the common area between two probability density functions (Perkins et al. 2007).In general, DAV = 0 indicates that no gain is found; DAV <0 points out a loss associated with using a higher resolution relative to the lower resolution; DAV >0 expresses the beneficial impact of increasing the grid spacing.
In a second step, more detailed information is provided, assuming the Nomenclature of Territorial Units for Statistics (NUTS1) classification 2021 (Fig. 1b) as the spatial reference.NUTS1 for Italy splits the peninsula into five sub-areas (i.e., Northwest Italy, Northeast Italy, Central Italy, South Italy, and Insular Italy).For each sub-area, the evolution of the annual bias, the orography role, and the assessment of extremes through climate indicators are investigated.To address the orography role, the 30-year mean temperature and precipitation values are first grouped into 300-metre orographic bands and then averaged.Such an operation has been performed for observations and simulation independently.The evaluation of extremes relies on climate indicators for temperature and precipitation (see Table 1), selected among those provided by the experts of the CCl/CLIVAR/ JCOMM Team on Climate Change Detection and Indices (ETCCDI; Karl et al. 1999).These indicators generally turn climate data into information highlighting various characteristics of extremes, including frequency, amplitude, and persistence.In this work, some indicators are used following the precise definition provided by ETCCDI, although others are slightly modified from the originals.Operatively, some frequency indicators (i.e., FD, SU, TR, R20mm, CWD and CDD) are calculated yearly and then averaged over 30-year; others, such as the percentiles and the absolute indicators (i.e.TX90th, TN10th, RX1d and R99th), are instead calculated considering the entire distribution of values for the 30-year.Once calculated, each indicator is spatially aggregated at the NUTS1 level (see Fig. 1b), deriving spatial mean and spread between the 5th and 95th percentile of the spatial distribution.The outcomes are displayed using a "radar chart" representation.
Then, further separate analyses are performed for temperature and precipitation, focusing on the summer season (JJA).Maps of minimum temperature, normalised using the Min-Max procedure (Patro and Sahu 2015), are proposed for the two biggest Italian cities (i.e., Milan and Rome) to preliminary appreciate the model's ability to capture the UHI phenomena.In this application, measurements from LWS (on which SCIA-ISPRA relies) are processed with the Ordinary Kriging technique (Wackernagel 1995) and assumed as the baseline.On the other hand, for precipitation, hourly maps of intensity, frequency, and heavy precipitation are presented to verify the potential benefits of the CP-RCM.Specifically, intensity and frequency are assessed on wet hours (an hour with precipitation ≥ 0.1 mm/h).In contrast, heavy hourly precipitation is investigated through the amount of hourly precipitation above the 99.9th percentile computed from all data (wet and dry events), according to Schär et al. (2016).The amount of hourly precipitation is also investigated using whiskers and box plot representations to look for differences in the statistical distributions synthetically: the halfway mark indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively.Moreover, diurnal cycles averaged over NUTS1 are assessed during summer for wet-hour frequency and intensity.
It figures out a general good ability of the climate simulation in reproducing observational spatial patterns.In particular, the simulated 2m temperature agrees with observations (both SCIA-ISPRA and E-OBS) over the Alps and Apennine chains and the mountainous areas in southern Italy ( as Basilicata, Calabria, and the Etna volcano).On the other hand, the simulation slightly overestimates temperature over flat terrains such as the Po valley, northwestern Italy (as South Piemonte), and coastal areas, especially if compared with E-OBS.Several studies found a temperature overestimation in the Po valley due to a typical atmospheric dynamic acknowledged as surface-based temperature inversions (Cacciamani et al. 1995;Bucchignani et al. 2016) not reproduced by most climate models.
Yearly statistics (i.e.mean values and standard deviations) shown in Figure 2 and reported in Table 2, confirm the good agreement between observations and VHR-REA_ IT.Specifically, the yearly mean value for observations is 13.2 °C for SCIA-ISPRA and 13.1°C for E-OBS, with a spatial variability of ± 3.8 °C for both observational datasets.In comparison, it is 13.7°C for VHR-REA_IT with a spatial variability of ± 4.2 °C.UERRA and ERA5 show a yearly mean value of 12.5°C and 12.9°C, respectively, with a spatial variability of ± 3.7°C and ± 3.9°C, respectively.By looking at seasonal values (see Table 2), VHR-REA_IT against SCIA-ISPRA returns a cold bias (-1.2 °C) during winter and a warm tendency (+1.8 °C) during summer.In contrast, the biases are attenuated during the transition seasons (i.e., autumn and spring, with values = +0.1 °C and +0.5 °C, respectively).Such a tendency also emerged in the seasonal comparisons against E-OBS (also reported in Raffa et al. 2021a), highlighting cold winter bias and warm summer bias.Similar tendencies occur in the comparison of VHR-REA_IT against UERRA and ERA5.The overestimation of VHR-REA_IT in summer temperature may be related to the soil moisture-precipitation feedback modification (Pielke 2001;Koster et al. 2004) acting when an explicit treatment of deep-convection processes is turned on (Hohenegger et al. 2009;Taylor et al. 2013).Indeed, a different soil moisture-precipitation feedback influences the heat flux partitioning, affecting summer temperature.Such behaviour seems to be revealed by several CP-RCMs, as highlighted in Sangelantoni et al. (2023).
Moreover, the spatial patterns in Figure 2 show that UERRA reanalysis slightly underestimates 2m temperature over the whole domain if compared with SCIA-ISPRA observations fitting in a good way E-OBS, whereas ERA5 patterns are in better agreement with SCIA-ISPRA than E-OBS, especially over flat terrain.Such qualitative analyses about the reanalysis products performance are also confirmed by the DAV metric reported in Table 3, where SCIA-ISPRA is used as the reference.Generally, DAV assessed over Italy indicates an added value in VHR-REA_IT of 14% if compared with UERRA reanalysis and a slight loss of -4% if compared with ERA5.The DAV over NUTS1 indicates that an added value in VHR-REA_IT is found considering both UERRA and ERA5 over Insular, Central and Southern areas.In contrast, slight losses are quantified over Northern areas, mainly due to temperature overestimation of VHR-REA_IT in the Po Valley.
In addition, Figure 3 shows the yearly evolution of 2m temperature bias over NUTS1 subdomains.Such elaboration confirms a warm bias for each NUTS with values between 0.4 and 1.4 °C.The lower limit of this bias is in the Northwest area.It is also interesting to note the evolution of the bias in the Insular area, where the trend over time is in line with the other regions from 2003 onwards.In contrast, the evolution appears to be inconsistent in previous years.It may be due to the low potential reliability of LWS observations on which the SCIA-ISPRA grid relies.This point would need further investigation.

Mapping of urban heat during summer nights
The increase in built surfaces constitutes the main reason for the formation of urban heat, as urban canyons preclude the release of reflected radiation.The main contribution to the formation of urban heat is the missing night-cooling of horizontal surfaces, together with a cloudless sky and light  winds (Oke et al. 2017).Of course, there is also a contribution from indoor heating, vehicle presence, and waste heat from air conditioning and refrigeration systems.The COSMO-CLM model can cope with this effect by modelling urban environments with the bulk parameterisation scheme (i.e., TERRA-URB), exploiting prescribed anthropogenic heat flux, impervious surface area, and other external urban canopy parameters (Garbero et al. 2021).In general, cities store a high amount of heat during the daylight hours warming the town elements (i.e., roads, walls and roofs) that release this trapped heat during the night hours, warming the air within the town.Therefore, the released heat is maximised during the summer night hours (Reder et al. 2018).
So the analyses on the minimum temperatures are proposed in the following associated with tropical night occurrence.
Figure 5 shows the maps of minimum temperature, normalised using the Min-Max procedure, for the cities of Milan (Fig. 5a for observation SCIA_ISPRA LWS, Fig. 5b for reanalysis product ERA5 and Fig. 5c for the simulation) and Rome (Fig. 5d for observations SCIA_ISPRA LWS, Fig. 5e for reanalysis product ERA5, and Fig. 5f for the simulation) to track the urban heat during summer nights.The Min-Max procedure is a normalisation method by performing linear transformations of the original data (Kappal 2019).Specifically, the minimum value of data is transformed into 0 while the maximum to 1 so that other values gets transformed into a decimal between 0 and 1.
Observations capture the urban heat island effect, gradually disappearing as it moves away from the more urbanised area.ERA5, also due to its coarse resolution, is not able to reproduce the urban heat island effect while the results of the dynamical downscaling consistently highlight the ability of urban parameterisation to represent urban dynamics and thermal peculiarities compared to observations.This result is particularly evident for Milan.Finally, it is interesting to note that the warm blob detected by VHR-REA_IT (Fig. 5f) and identified with the triangular symbol, is due to the presence of a body of water (i.e. the Lake Bracciano) which reduces the diurnal air temperature variation (i.e., maximum and minimum temperature close to the mean value), returning higher minimum temperature than that of the areas surrounding the city.

Extreme temperature indicators
The last section for temperature evaluation focuses on analysing extreme temperature indicators introduced in §2 and listed in Table 1. Figure 6 shows the temperature-related ETCCDI climate indicators for 1989-2020 represented through the radar chart for each NUTS1.The plots provide mean (continuous line), 5th percentile (filled circle) and 95th percentile (filled square) values of indicators allowing the comparison between the VHR-REA_IT and SCIA-ISPRA datasets.The values are also provided in Table 4 as mean value differences between VHR-REA_IT and SCIA-ISPRA.
VHR-REA_IT generally overestimates the 90th percentile of daily maximum temperature (TX90th, see Fig. 6a) against SCIA-ISPRA observations with almost uniform To sum up, the indicators TX90th and TN10th are uniform and well reproduced by VHR-REA_IT.The same is for summer days and tropical nights.A response dependent on the region, especially in topographically complex areas, is performed by frost days magnifying differences between Northern and Southern areas, resulting in a non-uniform radar plot.
Although covering different periods, the patterns of precipitation from SCIA-ISPRA and GRIPHO observational datasets (Fig. 7a and c) depict a coherent and uniform spatial distribution.Specifically, mean precipitation is below 4 mm/ day, almost the entire peninsula except for the Alpine chain and regions, such as the northern Tyrrhenian coast, featuring more pronounced mean values.Moreover, GRIPHO highlights higher values than SCIA-ISPRA over the southern Mediterranean area (as the Neapolitan coast and coastal areas in the southern Campania region).E-OBS (Fig. 7b) reports lower values than the previous one, especially over the Alps chain and Insular areas.VHR-REA_IT (Fig. 7f) agrees with observations, SCIA-ISPRA and GRIPHO in particular, reproducing the spatial pattern of precipitation over the whole domain, except over the Alpine area, where the model tends to produce higher daily precipitation values.This behaviour may be partly due to the forcing data (i.e.ERA5 depicted in Fig. 7e) and partly to the snow effect (detailed in §4.2.2). Figure 7e proves that ERA5 tends to increase precipitation across Italy, especially if compared with E-OBS gridded dataset.The spatial refinement of ERA5 in VHR-REA_IT at 2.2 km reduces this wet bias.Moreover, the CP-RCMs multi-model study by Ban et al. (2021) also found these overestimations over the Alps, highlighting that daily precipitation has higher values than observations in large-scale precipitation from mid-latitude storms in winter and convective precipitation in summer.
Table 5 reports the yearly and seasonal mean precipitation values derived from SCIA-ISPRA, E-OBS, GRIPHO, UERRA, ERA5 and VHR-REA_IT.Yearly metrics (i.e.mean values and standard deviations) shown in Table 5 confirm the coherence between observations and the model.The mean value is 2.6 mm/d for SCIA-ISPRA and UERRA and 2.7 mm/d for GRIPHO and ERA5, close to VHR-REA_IT (2.5 mm/d); instead, it is 2.2 mm/d for E-OBS.In addition, spatial variability is around ± 0.8 mm/d for all observational datasets, while it is slightly higher for UERRA and ERA5 (± 1.0 mm/d) and even higher for VHR-REA_IT (± 1.2 mm/d).
By looking at the seasonal mean precipitation values, the observational references indicate high precipitation in autumn and winter, with values falling in spring and attaining a minimum in summer.In comparison with observations (mainly SCIA-ISPRA and GRIPHO), VHR-REA_IT reduces precipitation in autumn and winter (3.0 mm/d and 2.3 mm/d, respectively), increasing it in spring (2.7 mm/d) and summer (2.0 mm/d), when convection plays an increased role resulting in higher precipitation value (Ban et al. 2021).
Moreover, the comparison in Figure 7 with reanalysis products highlights that the spatial patterns depicted by UERRA reanalysis slightly overestimate total precipitation in some parts of central and southern Italy if compared with SCIA-ISPRA observations, whereas ERA5 patterns are in good agreement with it.Such qualitative analyses are also confirmed by the DAV metric reported in Table 6, where SCIA-ISPRA is used as the reference.Generally, DAV assessed over all Italy indicates an added value in VHR-REA_IT of 7% if compared with UERRA reanalysis and a loss of −5% if compared with ERA5.The DAV over NUTS1 indicates that an added value in VHR-REA_IT is found considering both UERRA and ERA5 over Northern and Southern areas, whereas losses are quantified over Central and Insular areas.These negative DAV values may be related to peculiar circulation patterns occurring, especially along coasts (i.e., where the sea mainly regulates the lateral boundaries, and CCLM is not able to reproduce these dynamics as it is uncoupled from an ocean model; Reder et al. 2022).
Figure 8 draws the yearly evolution of precipitation bias over NUTS1 subdomains.The time-series return a precise spatial repartition of the bias along the Italian peninsula.Specifically, northern Italy (i.e., Northwest NUTS1 and Northeast NUTS1) reveals a slight increase in the simulated total precipitation against observations; such bias is on average 0.5 mm/day with a temporal standard deviation inbetween 0.1-0.2mm/day.On the other hand, the rest of the peninsula points out a slight decrease in the simulated total precipitation; this decrease is limited in the South NUTS1 (≈ −0.3 mm/day) growing in the Central NUTS1 and Insular NUTS1 (≈ −0.6 mm/day).Finally, the temporal standard deviation for these areas is almost the same (≈ 0.2 mm/day).

Orography effects
Figure 9 displays the altitudinal profiles of mean total precipitation obtained for each NUTS1 (Fig. 9a-e) using SCIA-ISPRA and VHR-REA_IT.The altitudinal profiles obtained by interpreting the VHR-REA_IT total precipitation (red line) follow the observation patterns (black line) up to altitudes of around 1000 m and then diverge towards significantly higher values beyond this altitude.Between the bands 1000 m and 4000 m, snow accumulation is likely responsible for the wider spread (Monteiro et al. 2022).Indeed, COSMO-CLM total precipitation includes rainfall and snowfall (i.e., liquid and solid precipitation), while the SCIA-ISPRA observations predominantly incorporate rainfall based on rain gauge measurements.To account for such a difference in the quantification of the total precipitation between the simulation and observation, this work proposes a simplified method relying on temperature data to discern the snowfall component of the modelled total precipitation as the orography increases.Such a method is a modified version of the approach used at the monthly scale in the Thornthwaite water balance model (Thornthwaite and Mather 1955;McCabe and Markstrom 2007) to derive rainfall and snowfall from total precipitation according to mean temperature.It is assumed that the total precipitation is totally snowfall (hence rainfall contribution null) when  both maximum and minimum daily 2m-temperatures are lower than zero.At the same time, it is totally rainfall (thus snowfall contribution null) when both are greater than zero.
In the case of maximum daily 2m-temperature greater than zero and minimum one lower than zero, a weighted average law for total precipitation -based on mean, maximum and minimum temperatures -is defined to account for both rainfall and snowfall contributions.More details can be found in Appendix, where the procedure is outlined and validated using ERA5 data.
The method provides the snowfall profiles in Fig. 9a-e (grey line) and, consequently, the rainfall ones (blue line) as the difference between the total precipitation and snowfall profiles.The Figure shows snowfall has a negligible contribution for low elevations (less than 1000 m) and increases with altitude, exceeding rainfall contributions for altitude values greater than 2500 m.On the other hand, VHR-REA_IT rainfall returns altitudinal profiles that match the observations satisfactorily.Specifically, rainfall increases with altitudes up to 1000 m.It then reverses the altitudinal growth gradient beyond this altitude when part of the precipitation occurs as snowfall.This behaviour also clarifies the precipitation overestimates shown in Fig. 7f over high orography (i.e., the Alps) and is even more evident in Fig. 8 in terms of positive mean annual bias for Northwest Italy and Northeast Italy.

Analysis of the summer hourly precipitation
Several studies (Ban et al. 2021;Adinolfi et al. 2021;Kendon et al. 2021;Lucas-Picher et al. 2021) recognised the improved representation of hourly precipitation as an addedvalue of CP-RCMs.In this view, this Section analyses hourly precipitation features of VHR-REA_IT, comparing with the driven reanalysis ERA5 and assuming the available hourly dataset (i.e., GRIPHO) as a reference.
VHR-REA_IT overestimates intensity compared to observations (Fig. 10g against Fig. 10a) over Italy, except for the Alps, where it well-reproduces hourly precipitation intensity.The modelled frequency (Fig. 10h) agrees with the spatial patterns of observations (Fig. 10b), with higher values over the Alps than the rest of the peninsula.The 99.9th percentile of hourly precipitation (representing heavy precipitation) is fairly modelled by VHR-REA_IT (Fig. 10i) against observation (Fig. 10c) over the Alps.Still, it is overestimated in the Southern part and the insular areas.Ban et al. (2021) and Adinolfi et al. (2021) returned similar findings by comparing observations with CP-RCMs over a domain covering Northern and Central Italy.Specifically, Ban et al. (2021) argued that overestimations in intensity and 99.9th percentile of hourly precipitation are not necessarily signs of a bad model performance, as a part of these differences can be explained by observational uncertainties due to the sparse observational networks over higher altitudes that receive more precipitation and interpolation methods used to produce gridded datasets.Moreover, Figure 10 depicts that ERA5 (Fig. 10d-f) strongly overestimates hourly precipitation frequency, whereas it underestimates precipitation intensity and heavy precipitation, in comparison with GRIPHO (Fig. 10a-c) and VHR-REA_IT (Fig. 10g-i).Reder et al. (2022) returned similar findings analysing precipitation extremes for different cities at CP-RCMs.They stated that each grid point of ERA5 is informative for an area, at least equal to 31 × 31 km 2 , and then, for each grid cell, weather variables (i.e., for precipitation, Reder and Rianna 2021) can be affected in their magnitude by an averaging effect which produces a smoothing, especially for hourly extreme precipitation values.Such a smoothing effect decreases with spatial resolution enhancement.This aspect returns evidence of benefits induced by the convection-permitting model in representing hourly precipitation features.
In addition, heavy precipitation is statistically analysed as whiskers box plots in Figure 11 at the NUTS1 level.It is worth noting that VHR-REA_IT data is remapped over the GRIPHO grid, so the remapping procedure could induce artificial errors due to smoothing or representativity errors.The box plots from GRIPHO return Northwest and Northeast median values of around 11 mm/h.VHR-REA_IT represents with a good agreement the Northwest and overestimates Northeast NUTS1 with a shift of around 3 mm/h.The median values in the Insular area are close to 5 mm/h for both observation and model.The main differences between model and observation occur over the South area, in representing the median and 95th percentile, and in the 5th percentile for the Central NUTS1.To sum up, VHR-REA_IT increases the distribution of heavy hourly precipitation in median and percentiles.
The findings highlight that the spatial refinement of VHR-REA_IT is able to reproduce the hourly precipitation features more in the complex orographic areas (i.e., Northwest and Northeast NUTS1) than in other territories (i.e.South, Central and Insular NUTS1).They also agree with previous studies (Berthou et al. 2020;Adinolfi et al. 2021) and multi-model ensemble (Ban et al. 2021) at the convectionpermitting scale over the Alpine area.
Moreover, diurnal cycles of precipitation intensity and wet-hour frequency are averaged in summer (2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016) over NUTS1, as shown in Figures 12 and 13, respectively.The km-scale patterns improve the representation of diurnal cycles both for precipitation intensity and wet-hour frequency.For intensity (Fig. 12), the values and timing of GRIPHO over all NUTS1 are much better captured by VHR-REA_IT than ERA5.The latter systematically underestimates the hourly intensity.When looking into the diurnal cycle of wet-hour frequency (Fig. 13), VHR-REA_IT figures out a good match with GRIPHO, especially over Northeast and Northwest NUTS1.ERA5 systematically overestimates the wet-hour frequency over all NUTS1.The diurnal cycles of summer intensity and wet-hour frequency reveal that VHR-REA_IT shows superior performance to the coarse resolution ERA5 Reanalyses.It is a further confirmation that the longstanding problem of incorrect timing with too frequent and too weak precipitation is greatly improved by switching on the parameterisation of convection (Ban et al. 2021).However, some deficiencies persist in the CP-RCMs modelling system that need to be addressed in the future.

Extreme precipitation indices
Figure 14 displays the precipitation-based ETCCDI climate indicators, introduced in §2 and listed in Table 1.Data are represented for 1989-2020 using the radar chart for each NUTS1.The plots provide mean (continuous line), 5th percentile (filled circle) and 95th percentile (filled square) values of indicators allowing the comparison between the VHR-REA_IT and SCIA-ISPRA datasets.The values are also provided in Table 7 as mean value differences between VHR-REA_IT and SCIA-ISPRA.VHR-REA_IT fairly and uniformly detects the observed distribution of RX1d (Fig. 14a) and R99th (Fig. 14b) as for mean as for percentiles, although for R99th slight overestimations over South and Insular areas occur.Moreover, the R20mm (Fig. 14c) distribution is quite well represented by VHR-REA_IT compared with observations over the Insular, Southern and Central areas, and some overestimations are represented over the Northeast and Northwest areas.On the other hand, for CWD (Fig. 14d) and CDD (Fig. 14e), simulation differs from observations.Specifically, VHR-REA_IT returns a shifted distribution for CWD compared to observations over all NUTS1, but Northeast.It means that the simulation underestimates both mean and percentile values.Consequently, VHR-REA_IT generally overestimates CDD against observations.This is mainly evident in the Insular, South and Central areas.
To sum up, maximum daily precipitation and 99 th percentile of daily precipitation are uniform and well reproduced by VHR-REA_IT.The same is for the number of days with precipitation ≥ 20 mm.A response dependent on the region, especially in topographically complex areas, is performed by the maximum number of consecutive wet and dry days, highlighting a better performance of VHR-REA_IT in  mountainous Northern NUTS1 and differences in the South, Central and Insular NUTS1 with non-uniform radar plot.Zollo et al. (2016) performed a similar analysis for precipitation extremes over Italy (in particular in four regions, i.e., Veneto, Tuscany, Calabria, and Sardinia, which can be assumed as indicative of Northeast, Central, Southern and Insular NUTS1).They exploited a simulation performed with the RCM COSMO-CLM at a horizontal resolution of ~ 8 km driven by the ERA-Interim reanalysis (Dee et al. 2011).In general, the authors highlighted a general tendency of the model towards underestimation, except for the number of CDD.Compared with Zollo et al. (2016), this work improves the representation of RX1d, R99th, and R20mm by reducing the bias between the model and observations.However, some deficiencies still need improvement in representing the frequency of wet/dry days.As done in this work, adopting a threshold of 1 mm to discriminate between wet and dry days may partly mitigate such an issue.Although CP-RCMs largely improve the tendency (Figs. 12,13) to trigger excessively frequent and weak precipitation events (the so-called "drizzle effect" as in Gutowski et al., 2003) in the coarser resolution models (Pichelli et al. 2021), CP-RCMs present drier conditions than the corresponding RCMs (Sangelantoni et al. 2023), and it is not possible to fully disentangle the role of resolution increase against the role of improved physics, which are inextricably linked.Both the land-atmosphere interactions and microphysical parameterisations could be regarded as the key factors behind these modifications introduced by the km-scale modelling (Milovac et al. 2016).For example, large uncertainty affects the representation of the soil moisture-precipitation feedback, which remains a challenge (Yang et al. 2018;Kendon et al. 2021;Leutwyler et al. 2013).

Conclusions
The paper analysed VHR-REA_IT, a dynamical downscaling of ERA5 recently performed over Italy through the COSMO-CLM model at a convection-permitting scale, turning on the urban scheme TERRA-URB.The atmospheric variables coming from VHR-REA_IT (introduced in Raffa et al. 2021a) provide accurate information on the actual local synoptic weather situation.Although Raffa et al. (2021a) is a data descriptor paper, it presented some preliminary and seasonal analyses (i.e.spatial patterns and statistics) with respect to E-OBS and ERA5.The current study integrated Raffa et al. (2021a) analyses with further statistical evaluations, metrics and reference products, focusing on 2m temperature and precipitation, assumed as the two variables commonly adopted in climate assessments.The capability of VHR-REA_IT in realistically representing the climatology over the whole domain and subdomains (NUTS1) as well as capturing specific features, has been tested.The main results of the study are summarised in the following: 1.The first achievement is the assessment of climate data over Italy with increased space (2.2 km) and time (1 h) resolutions for a 30-year long period .Such continuous reanalysis-driven downscaling has the advantage of overcoming the common time-slice approach (i.e. based on 10-year period of simulation) limitations, usually used due to the computational burden, therefore it goes some steps forward compared with the current state-of-the-art.2), especially over the mountainous areas.On the other hand, it returns a warm temperature bias over flat terrain mainly due to the soil moisture-precipitation feedback modification acting when an explicit treatment of deep-convection processes is turned on (Hohenegger et al. 2009;Taylor et al. 2013).This issue is a well-known and an independent-model behaviour investigated in CP-RCMs multi-model studies (Sangelantoni et al. 2023;Soares et al. 2022), revealing the need to assess the applicability of parametrisations, designed for coarser resolutions to kilometre scales (as the land-atmosphere interactions and microphysical parameterisation).Furthermore, VHR-REA_IT consistently highlights the ability of urban parameterisation to represent urban heat as the local detection of minimum temperature during summer.This result is particularly evident in the case of Milan city (see Fig. 5).Such a capacity makes the urban parameterization used for the dynamical downscaling VHR-REA_IT a powerful tool to improve regional climate models, support local adap-tation actions and impact studies in the urban environment.
The spatial distribution of daily mean precipitation reveals that VHR-REA_IT accurately represents past and current precipitation patterns and agrees with observations over the whole domain, except the Alpine area, where the model returns higher values.It could be possibly due to the forcing data and partly to the snow effect.An interestingalthough simplified -approach is proposed to account for the snowfall component of the total precipitation at high altitudes: VHR-REA_IT rainfall returns vertical profiles that match the observations satisfactorily and increase with altitudes up to 1000 m.Beyond this altitude, part of the precipitation occurs as snowfall.This behaviour also clarifies the precipitation overestimates over the Northern part of the peninsula.Moreover, hourly analyses show that the modelled intensity of wet hours in summer and the heavy precipitation events, in terms of the 99.9th percentile, are overestimated if compared with reference products, especially over the South part of the peninsula.In contrast, over the Alps, the agreement is satisfactory.The modelled frequency of wet hours in summer agrees with observations.Both simulation and observations depict higher values over the Alps than in the remaining part of the peninsula.The ability of VHR-REA_IT is mainly in the representation of hourly intensity, frequency and heavy precipitation if compared with ERA5 patterns (see Figure 10), clearly showing the benefits of explicitly-solved convection.Such capability is also proved by the diurnal cycles (Figs. 12,13), where ERA5 systematically underestimates the hourly intensity and overestimates the wet-hour frequency over all NUTS1, but VHR-REA_IT figures out a good match with GRIPHO in the representation of values and timing.Therefore, this study corroborates the advantages of the use of CP-RCMs, as Prein et al. (2013), Kendon et al. (2021), Adinolfi et al. (2021) and multi-model studies (Ban et al. 2021;Müller et al. 2022).

The third achievement consists in assessing extreme
events on a detailed scale through climate indicators.It is worth noting that climate indicators are useful tools for the assessment of extreme events in terms of magnitude, occurrence and impacts (Kjellström et al. 2010;Devkota and Bhattarai 2018).For instance, the number of consecutive dry days constitutes essential information in drought studies.At the same time, the maximum daily precipitation may provide information supporting, as an example, flooding and landslides studies.Finally, the number of tropical nights is closely connected to the occurrence of heatwaves.In recent years, assessing the added value of CP-RCMs in the representation of precipitation patterns, compared to coarser resolution RCMs, has been one of the main topics of climate studies, as evidenced by recent initiatives such as the CORDEX FPS-CONV.This paper aimed to support the research activities with this goal, proving the limitations and abilities of VHR-REA_IT in representing past and current 2m temperature and precipitation patterns over complex contexts.This aim was addressed by analysing a 30-years long simulation against ERA5 and other reference products (i.e.observational gridded datasets and reanalyses).The main impact of this study is mainly intended for research purposes in the field of climate studies.Thus, the findings cannot properly support decision-making and adaptation assessments, as results are derived from one-single dynamical downscaling and affected by model uncertainties.To cope with this issue, multi-model ensembles are necessary to span the full range of uncertainties (Fosser et al. 2015;Ban et al. 2021).However, a multi-model long-time CP-RCMs ensemble derived from a coordinated initiative over Italy is currently unavailable.The key scientific results presented in this paper highlighted the benefits and limitations of CP-RCMs, contributing further to the current stateof-the-art on this topic.In particular, benefits are gained in mountainous areas where the CP-RCMs are proven to fairly represent the 2m temperature and hourly precipitation features.Over urban contexts, the here adopted urban parameterisation is proved to simulate city heat detection.Losses in performances occur over coastal and flat areas of the Italian peninsula, where the performance of the CP-RCM does not seem to be satisfactory.The presented potentialities and limitations are affected by the peculiarities of the Italian peninsula, characterised by orographic features, terrain close to the coast and complex urban development.
Moreover, the proposed findings would support the scientific activities in the climate community to find a proper trade-off between the improvements and obstacles of running CP-RCMs.On the one hand, the fine spatial and temporal scales and the demonstrated benefits make CP-RCMs attractive tools to represent climate features.They open up a wide field of research applications related to the impacts of climate change.On the other hand, the limitations of the CP-RCMs over some territories raise the question of whether the additional computational cost required by CP-RCMs is justified.Specifically, VHR_ REA-IT required 2160 cores and employed about 61 h to perform a 1-year simulation (Raffa et al. 2021a), producing a large amount of data (around 8 TB of output data and around 70 TB of forcing data in terms of 3-dimensional boundary data needed for the downscaling).
As future developments, the role of the data assimilation procedure on downscaling activities (not included in VHR-REA_IT) could be addressed by comparing it with regional reanalyses over Italy driven by ERA5 (Bonanno et al. 2019, Capecchi et al. 2022, Giordani et al. 2023and Cerenzia et al. 2022).
Given the coefficient k, it is possible to determine: The introduced approach is validated considering the snowfall data from the ERA5 reanalysis.For this purpose, 13 years of data (i.e., 2002-2014) have been processed.The validation consists of comparing the snowfall returned by the approach with the snowfall returned by ERA5.The analysis has been performed on the entire Italian domain and the results are shown both as monthly time series (Fig. 15 in Appendix) and scatter plots (Fig. 16 in Appendix).
The evolutions obtained with the simplified approach agree with those returned in output by ERA5.In particular, the rainfall variable shows a coefficient of determination R2 practically close to unity (0.98).At the same time, very satisfactory performances are returned by the snowfall variable, characterised by a coefficient of determination R2 of approximately 0.88.This simple comparison proves the validity of the simplified approach adopted, which makes it possible to easily discriminate the rainfall and snowfall components from the total precipitation on a daily scale.
Acknowledgements The authors sincerely thank the HIGHLANDER (HIGH performance computing to support smart LAND sERvices) project (https:// highl ander proje ct.eu/), which is funded by the Connecting European Facility (CEF) Telecommunications sector under agreement number INEA/CEF/ICT/A2018/1815462.The experiments have been performed using the COSMO model in Climate Mode (COSMO-CLM).COSMO-CLM is the community model of the German regional climate research jointly further developed by the CLM-Community.The authors acknowledge the members of the community for their common efforts to develop the model and to find the right setups.The authors would like to thank Emanuela Pichelli and Erika Coppola from UNESCO ICTP for allowing access to the Italian database of precipitation GRIPHO and the WCRP-CORDEX-FPS on Convective phenomena at high resolution over Europe and the Mediterranean.Acknowledgment is also made for the Italian Institute for Environmental Protection and Research (ISPRA) that made the (A.4) SCIA-ISPRA dataset freely available by through public web access (scia.isprambiente.it),exploiting the national system for climate data collection, processing and dissemination.UERRA and ERA5 have been obtained using Copernicus Climate Change Service Information and downloaded from the C3S CDS.Finally, we acknowledge the E-OBS dataset and the data providers in the ECA&D project ( https:// www.ecad.eu).The authors are also grateful to the anonymous reviews for their comments and suggestions, which enabled them to improve the manuscript.

Funding
The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.

Fig. 1 a
Fig. 1 a Surface height (m a.s.l.) of the computational domain; b area considered for validation defined according to the first level of Nomenclature of Territorial Units for Statistics (NUTS 1) for Italy

Fig. 3 Fig. 4
Fig. 3 Evolution of the annual bias of VHR-REA_IT against SCIA-ISPRA dataset of 2 m mean temperature over NUTS1

Fig. 5
Fig. 5 Spatial patterns of the minimum temperature during summer (representative of night-time conditions) over Milan urban context retrieved by SCIA-ISPRA local weather stations (a), ERA5 (b) and VHR-REA_IT (c).Respectively subplots (d), (e) and (d) refers to Rome urban context.Minimum temperature in JJA have been normal-

Fig. 6
Fig. 6 Radar chart for temperature-related ETCCDI climate indicators for 1989-2020 for each NUTS1.The plots provide mean (continuous line), 5th percentile (filled circle) and 95th percentile (filled square) values of indicators allowing the comparison between the

Fig. 8 Fig. 9
Fig. 8 Evolution of the annual bias of VHR-REA_IT against SCIA-ISPRA dataset of mean total precipitation over NUTS1

Fig. 11
Fig. 11 Box plots over NUTS1 for the period 2001-2016 in terms of summer heavy precipitation assessed as the 99.9 percentile of hourly precipitation from GRIPHO and VHR-REA_IT

Fig. 14
Fig. 14 Radar chart for precipitation-related ETCCDI climate indicators for 1989-2020 for each NUTS1.The plots provide mean (continuous line), 5th percentile (filled circle) and 95th percentile (filled square) values of indicators allowing the comparison between the VHR-REA_IT (red lines and fills) and SCIA-ISPRA (black lines and All authors contributed to the study conception and design.Conceptualization and methodology, AR, MA, MR and PM; software, MR; validation, MR, MA and AR; formal analysis, MR, AR; investigation, PM and AR; data curation, MR, MA and AR; writing-original draft preparation, MA; writing review and editing, MA, AR and PM; visualization, MR; supervision, PM All authors have read and agreed to the published version of the manuscript.

Table 1
List of indicators chosen for temperature and precipitation Number of days with precipitation greater or equal to 20 mmR20mm day/y Maximum number of consecutive wet days with daily precipitation greater than 1 mm CWD day/y Maximum number of consecutive dry days with daily precipitation lower than 1 mm CDD day/y

Table 3
DAV metric for 2 m temperature assessed over Italy (whole domain) and over NUTS1 from VHR-REA_IT (as high resolution product), SCIA (as observation), UERRA and ERA5 (as low resolu-

Table 5
Metrics (mean value and standard deviation) over the whole domain from SCIA-ISPRA, E-OBS, GRIPHO, UERRA, ERA5 and VHR-REA_IT, related to the spatial precipitation distribution for the period 1989-2020 (the period 2001-2016 is considered for GRIPHO) Metrics are yearly and seasonally (SON September-October-November, DJF December-January-February, MAM March-April-May, JJA June-

Table 7
VHR-REA_IT minus SCIA-ISPRA mean values of the selected climate indicators for precipitation over NUTS1