As a semi-arid Mediterranean country, Tunisia is affected by the impacts of climate change, particularly the coastal regions like the Cap-Bon. Irregular rainfall, rising temperatures and the recurrence of extreme events are all indicators that affect ecosystems and populations and make them more vulnerable to the influence of climatic variables. Therefore, an analysis of the trends of climate variables can contribute to facilitating the development of effective adaptation strategies. In this matter, this study was conducted to assess the homogeneity and trends of minimum and maximum air temperature (Tmin and Tmax) and precipitation (P) in the Cap-Bon region. Daily data were collected from the meteorological station of Nabeul for the period of 1982–2020. Pettitt and SNHT tests for homogeneity were applied to identify the breakpoints in multi-time scales of Tmax, Tmin and P data series. The Mann–Kendall (MK) test was used to detect the change in the time-series trend. A modified Mann–Kendall (mMK) test was used to remove the autocorrelation effect from the data series. Both the MK and mMK tests were used at the 5% significant level. The magnitude of the climatic trend was estimated using the non-parametric Sen’s slope estimator. Contrary to Tmin and P, the results of the homogeneity tests revealed the existence of significant breakpoints in the annual, seasonal and monthly Tmax time series. For most cases, the breakpoint occurred around the year 2000. For Tmin, significant breakpoints were recorded in March and April, while a significant shift in the P time series was detected in December. The Mann–Kendall results show a significant warming trend in annual Tmax, with magnitudes equal to 0.065 and 0.045 °C/year before and after the breakpoint, respectively. Nevertheless, non-significant tendencies were observed in the annual Tmin and P time series. On the monthly time scale, Tmax exhibited a significant upward trend in June and August, before the observed breakpoints, with Sen’s slope values equal to 0.065 and 0.045 °C/year, respectively. Regarding the Tmin data, a significant positive trend was observed in July at a rate of 0.033 °C/year.

Systematic biases in general circulation models (GCM) and regional climate models (RCM) impede their direct use in climate change impact research. Hence, the bias correction of GCM-RCMs outputs is a primary step in such studies. This study compares the potential of two bias correction methods (the method from the third phase of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP3) and Detrended Quantile Matching (DQM)) applied to the raw outputs of daily data of minimum and maximum air temperatures and precipitation, in the Cap-Bon region, from eight GCM-RCM combinations. The outputs of GCM/RCM combinations were acquired from the European branch of the coordinated regional climate downscaling experiment (EURO-CORDEX) dataset for historical periods and under two representative concentration pathway (RCP4.5 and RCP8.5) scenarios. Furthermore, the best combination of bias correction/GCM-RCM was used to assess the impact of climate change on reference evapotranspiration (ET₀). Numerous statistical indicators were considered to evaluate the performance of the bias correction/historical GCM-RCMs compared to the observed data. Trends of the Hargreaves–Samani_ET₀ model during the historical and projected periods were determined using the TFPMK method. A comparison of the bias correction methods revealed that, for all the studied model combinations, ISIMIP3 performs better in reducing biases in monthly precipitation. However, for Tmax and Tmin, the biases are greatly removed when the DQM bias correction method is applied. In general, better results were obtained when the HadCCLM model was used. Before applying bias correction, the set of used GCM-RCMs projected reductions in precipitation for most of the months compared to the reference period (1982–2006). However, Tmin and Tmax are expected to increase in all months and for the three studied periods. Hargreaves–Samani ET₀ values obtained from the best combination (DQM/ HadCCLM) show that RCP8.5 (2075–2098) will exhibit the highest annual ET₀ increase compared to the RCP4.5 scenario and the other periods, with a change rate equal to 11.85% compared to the historical period. Regarding spring and summer seasons, the change rates of ET₀ are expected to reach 10.44 and 18.07%, respectively, under RCP8.5 (2075–2098). This study shows that the model can be used to determine long-term trends in ET₀ patterns for diverse purposes, such as water resources planning, agricultural crop management and irrigation scheduling in the Cap-Bon region.