Over the past few decades the horizontal resolution of regional climate models (RCMs) has steadily increased, leading to a better representation of small-scale topographic features and more details ...in simulating dynamical aspects, especially in coastal regions and over complex terrain. Due to its complex terrain, the broader Adriatic region represents a major challenge to state-of-the-art RCMs in simulating local wind systems realistically. The objective of this study is to identify the added value in near-surface wind due to the refined grid spacing of RCMs. For this purpose, we use a multi-model ensemble composed of CORDEX regional climate simulations at 0.11° and 0.44° grid spacing, forced by the ERA-Interim reanalysis, a COSMO convection-parameterizing simulation at 0.11° and a COSMO convection-resolving simulation at 0.02° grid spacing. Surface station observations from this region and satellite QuikSCAT data over the Adriatic Sea have been compared against daily output obtained from the available simulations. Both day-to-day wind and its frequency distribution are examined. The results indicate that the 0.44° RCMs rarely outperform ERA-Interim reanalysis, while the performance of the high-resolution simulations surpasses that of ERA-Interim. We also disclose that refining the grid spacing to a few km is needed to properly capture the small-scale wind systems. Finally, we show that the simulations frequently yield the accurate angle of local wind regimes, such as for the Bora flow, but overestimate the associated wind magnitude. Finally, spectral analysis shows good agreement between measurements and simulations, indicating the correct temporal variability of the wind speed.
In this study we investigate the scaling of precipitation extremes with temperature in the Mediterranean region by assessing against observations the present day and future regional climate ...simulations performed in the frame of the HyMeX and MED-CORDEX programs. Over the 1979–2008 period, despite differences in quantitative precipitation simulation across the various models, the change in precipitation extremes with respect to temperature is robust and consistent. The spatial variability of the temperature–precipitation extremes relationship displays a hook shape across the Mediterranean, with negative slope at high temperatures and a slope following Clausius–Clapeyron (CC)-scaling at low temperatures. The temperature at which the slope of the temperature–precipitation extreme relation sharply changes (or temperature break), ranges from about 20 °C in the western Mediterranean to <10 °C in Greece. In addition, this slope is always negative in the arid regions of the Mediterranean. The scaling of the simulated precipitation extremes is insensitive to ocean–atmosphere coupling, while it depends very weakly on the resolution at high temperatures for short precipitation accumulation times. In future climate scenario simulations covering the 2070–2100 period, the temperature break shifts to higher temperatures by a value which is on average the mean regional temperature change due to global warming. The slope of the simulated future temperature–precipitation extremes relationship is close to CC-scaling at temperatures below the temperature break, while at high temperatures, the negative slope is close, but somewhat flatter or steeper, than in the current climate depending on the model. Overall, models predict more intense precipitation extremes in the future. Adjusting the temperature–precipitation extremes relationship in the present climate using the CC law and the temperature shift in the future allows the recovery of the temperature–precipitation extremes relationship in the future climate. This implies negligible regional changes of relative humidity in the future despite the large warming and drying over the Mediterranean. This suggests that the Mediterranean Sea is the primary source of moisture which counteracts the drying and warming impacts on relative humidity in parts of the Mediterranean region.
Flooding of the Adriatic coastline is predominantly caused by storm surges induced by winds from the south-eastern sector. This phenomenon in Venice is known as
acqua alta
. We present a study of ...wind fields favouring storm-surge setups in the Adriatic, including their characteristics in the present climate and their expected characteristics in future scenarios. Analysis is based on (i) measured sea levels in Venice and Bakar (1984–2014), (ii) near-surface wind from ERA5 reanalysis, and (iii) simulations of wind fields with three regional climate models (ALADIN52, RCA4, and RegCM4) forced with several global models (CNRM-CM, MPI-ESM-MR/LR, HadGEM2-ES, EC-EARTH, and IPSL-CM5). For future climates, we considered two scenarios (RCP4.5 and RCP8.5) and two future periods (2041–2070 and 2071–2100) with respect to the historical 1971–2000 period. It was found that the probability that the frequency, intensity, annual cycle, and spatial structure of the wind inducing the Adriatic storm surges will change in future climates is small. The result is robust and consistent according to all considered criteria—it does not depend on the analysed regional climate models, boundary conditions, climate scenarios, or future time interval.
The main goal of this analysis is to provide an assessment of the fine grid regional climate model (RCM) ensemble over the broader Adriatic region, which is the domain with orographically complex ...terrain and a developed coastline. The results of this study are a beneficial complement to earlier work on using modified bivariate bias correction methods. This analysis was based on an ensemble of 12 combinations of three RCMs and four global climate models (GCMs) from 1971 to 2004, divided into calibration and validation parts. A comparison of one univariate and three versions of bivariate bias correction methods of summer and winter monthly temperature and precipitation was performed. To examine the effect of various methods, we used different marginal distributions and distributions of dependence (copula). Bivariate bias correction was conducted with parametric and empirical marginal distributions. In the quantile mapping (univariate) and the parametric bivariate method, gamma (precipitation) and normal (temperature) probability distributions were used. We documented the bias and the influence of bias correction methods on precipitation and temperature, including their climate change signals in historical data. The considered methods maintained the spatial distribution of trends from the native ensemble. If the relationship between variables was not distinct and the correlation coefficient was low, bivariate bias correction methods tended to weaken the correlation coefficient from the model. Two experiments on how the bias correction methods influence statistical measures of the considered variables and their relationships showed that the bivariate method with empirical distributions (eeG) performed slightly better than other methods. Although the eeG method performed better than the other methods, none of the considered methods was able to fully preserve the correlation coefficient between variables detected from observations. Considering two copulas (Gaussian and Student's t), we conclude that the effect of different copula types on method performance is small.
Bivariate bias correction of the regional climate model ensemble over the Adriatic region Renata Sokol Jurković*, Ivan Güttler, Zoran Pasarić Bias in RCM ensemble over the Adriatic region was examined, and several bias correction methods and their influence on statistical measures, especially on the correlation coefficient, were evaluated. A comparison of one univariate and three versions of bivariate bias correction methods of summer and winter monthly temperature and precipitation was made. The eeG method maintains a positive correlation coefficient along the coast and lowers negative correlations detected from the native ensemble, and performs slightly better than other bias correction methods.
The ability of a large ensemble of regional climate models to accurately simulate heat waves at the regional scale of Europe was evaluated. Within the EURO-CORDEX project, several state-of-the art ...models, including non-hydrostatic meso-scale models, were run for an extended time period (20 years) at high resolution (12 km), over a large domain allowing for the first time the simultaneous representation of atmospheric phenomena over a large range of spatial scales. Eight models were run in this configuration, and thirteen models were run at a classical resolution of 50 km. The models were driven with the same boundary conditions, the ERA-Interim re-analysis, and except for one simulation, no observations were assimilated in the inner domain. Results, which are compared with daily temperature and precipitation observations (ECA&D and E-OBS data sets) show that, even forced by the same re-analysis, the ensemble exhibits a large spread. A preliminary analysis of the sources of spread, using in particular simulations of the same model with different parameterizations, shows that the simulation of hot temperature is primarily sensitive to the convection and the microphysics schemes, which affect incoming energy and the Bowen ratio. Further, most models exhibit an overestimation of summertime temperature extremes in Mediterranean regions and an underestimation over Scandinavia. Even after bias removal, the simulated heat wave events were found to be too persistent, but a higher resolution reduced this deficiency. The amplitude of events as well as the variability beyond the 90th percentile threshold were found to be too strong in almost all simulations and increasing resolution did not generally improve this deficiency. Resolution increase was also shown to induce large-scale 90th percentile warming or cooling for some models, with beneficial or detrimental effects on the overall biases. Even though full causality cannot be established on the basis of this evaluation work, the drivers of such regional differences were shown to be linked to changes in precipitation due to resolution changes, affecting the energy partitioning. Finally, the inter-annual sequence of hot summers over central/southern Europe was found to be fairly well simulated in most experiments despite an overestimation of the number of hot days and of the variability. The accurate simulation of inter-annual variability for a few models is independent of the model bias. This indicates that internal variability of high summer temperatures should not play a major role in controlling inter-annual variability. Despite some improvements, especially along coastlines, the analyses conducted here did not allow us to generally conclude that a higher resolution is clearly beneficial for a correct representation of heat waves by regional climate models. Even though local-scale feedbacks should be better represented at high resolution, combinations of parameterizations have to be improved or adapted accordingly.
This study performs a systematic analysis of the recent and future changes of dry spells (DS) in Croatia. DS are defined as consecutive sequences of days with daily precipitation less than 5 mm of ...the precipitation-per-day threshold (DS5). Daily precipitation data come from a dense national rain gauge network (covering seven regions) and span the period 1961–2015. The spatial and temporal changes of the observed mean (MDS5) and maximum (MxDS5) seasonal and annual dry spells were analysed by means of the Kendall tau method and the partial trend method. Future changes of DS5 were assessed by employing the three regional climate models (RegCM4, CCLM4 and RCA4) covering the EURO-CORDEX domain with a 12.5-km horizontal resolution, resulting in a realistic orography and land–sea border over Croatia. The models were forced at their boundaries by the four CMIP5 global climate models. For the reference period 1971–2000, the observed, as well as modelled, DS5 were analysed, and the systematic model errors were assessed. Finally, the projections and future changes of the DS5 statistics based on simulations under the high and medium greenhouse gases concentration scenarios (i.e., RCP8.5 and RCP4.5) with a focus on the climate change signal between 1971–2000 and two future periods, 2011–2040 and 2041–2070, were examined. A prevailing increasing trend of MDS5 was found in the warm part of the year, being significant in the mountainous littoral and North Adriatic coastal region. An increasing trend of MxDS5 was also found in the warm part of the year (both the spring and summer), and it was particularly pronounced along the Adriatic coast, while a coherent negative trend pattern was found in the autumn. By applying the partial trend methodology, an increase was found in the very long DS5 (above the 90th percentile) in the recent half of the analysed 55-year period in all seasons, except in the autumn when shortening in the DS5 was detected. The climate change signal during the two analysed future periods was positive for the summer in all regions, weakly negative for the winter and not conclusive for the spring, autumn and year. It was found that no RCM-GCM combination is the best in all cases, since the most successful model combinations depend on the season and location.
The 2m temperature (T2m) and precipitation from five regional climate models (RCMs), which participated in the ENSEMBLES project and were integrated at a 25-km horizontal resolution, are compared ...with observed climatological data from 13 stations located in the Croatian coastal zone. The twentieth century climate was simulated by forcing RCMs with identical boundary conditions from the ERA-40 reanalysis and the ECHAM5/MPI-OM global climate model (GCM); climate change in the twenty-first century is based on the A1B scenario and assessed from the GCM-forced RCMs’ integrations. When forced by ERA-40, most RCMs exhibit cold bias in winter which contributes to an overestimation of the T2m annual cycle amplitude and the errors in interannual variability are in all RCMs smaller than those in the climatological mean. All models underestimate observed warming trends in the period 1951–2010. The largest precipitation biases coincide with locations/seasons with small observed amounts but large precipitation amounts near high orography are relatively well reproduced. When forced by the same GCM all RCMs exhibit a warming in the cold half-year and a cooling (or weak warming) in the warm period, implying a strong impact of GCM boundary forcing. The future eastern Adriatic climate is characterised by a warming, up to +5 °C towards the end of the twenty-first century; for precipitation, no clear signal is evident in the first half of the twenty-first century, but a reduction in precipitation during summer prevails in the second half. It is argued that land-sea contrast and complex coastal configuration of the Croatian coast, i.e. multitude of island and well indented coastline, have a major impact on small-scale variability. Orography plays important role only at small number of coastal locations. We hypothesise that the parameterisations related to land surface processes and soil hydrology have relatively stronger impact on variability than orography at those locations that include a relatively large fraction of land (most coastal stations), but affecting less strongly locations at the Adriatic islands.
Since changes in temperature and precipitation have different effects on (a) all developmental stages of grapevines in most of the wine regions worldwide (i.e., on their phenological characteristics) ...and (b) different varieties, a comprehensive database of bioclimatic indices has been calculated and analysed for Croatian wine producing regions. The database consists of the average growing season temperature, growing degree‐days, Huglin index, dryness index and cool night index that are based on all available meteorological measurements as well as the outputs of regional climate models (RCMs) from the EURO‐CORDEX database. The horizontal grid spacing of 0.11° from the RCM ensembles enabled a fine‐scale determination of bioclimatic indices for the present and future climate in Croatia. In addition, statistical analyses (standard statistical parameters and Bayesian method) were carried out to examine trends in sugar content, total acidity and date of harvest. Calculations were performed for the present and future climate on the basis of data from seven selected vineyards/wineries and four varieties (‘Graševina’, ‘Plavac mali’, ‘Chardonnay’ and ‘Merlot’). The results show whether the part of Croatia that is suitable for grape cultivation in the present climate will continue to be favourable in the future within the Mediterranean area. In general, projections suggest further warming and drying of the climate in Croatia and an earlier harvest, with some variations among varieties that show latitude dependence. Projections for the future climate also suggest that the existing viticultural zoning will be much less adequate for the Croatian territory because it reduces the economically sustainable production of wine in certain areas.
Spatial distributions of temperature‐based bioclimatic indices (e.g., growing degree‐days GDD) for high‐end climate change scenarios are pointing to further warming in the period 2041–2070. Differences between the two periods P2 (2041–2070) and P0 (1971–2000) clearly show that some indices in certain regions will likely cross into higher (warmer) classes within the index scale. All grape varieties examined indicated an increase in the number of earlier harvests and a reduction in the number of later harvests, regardless of location.
The main goal of this study is to present a recently developed classification method for weather types based on the vorticity and the location of the synoptic centers relative to the Adriatic region. ...The basis of the present objective classification, applied to the Adriatic region, is the subjective classification developed by Poje. Our algorithm considered daily mean sea-level pressure and 500 hPa geopotential height to define one out of 17 possible weather types. We applied the algorithm to identify which weather type was relevant in the generation of the two typical near-surface winds over the Adriatic region, namely Bora and Sirocco. Two high-resolution (0.11°) EURO-CORDEX regional climate models were used, SMHI-RCA4 and DHMZ-RegCM4, forced by several CMIP5 global climate models and analyzed for two 30-year periods: near-present day and mid-21st century climate conditions under the high-end Representative Concentration Pathway (RCP8.5) scenario. Bora and Sirocco days were extracted for each weather type and a distribution over the 30-year period was presented. Our results suggest that in the winter season, climate model projections indicate a reduction in the main cyclonic types relevant in the formation of Bora over the entire Adriatic region and an increase in the number of anticyclonic types relevant in Sirocco events. In contrast, for the summer season, an increase in the main anticyclonic Bora-related weather types is found in the ensemble over the northern Adriatic region.