Hydrological studies depend heavily on environmental variables, such as precipitation and resulting runoff, which exhibit highly seasonal and intermittent behaviour in semiarid basins. In these basins, the use of traditional methods to adjust biases in time series projections can lead to inaccurate results regarding the impacts of climate change. This study introduces a non-stationary bias adjustment methodology (NS) specifically designed for environmental variables characterized by sporadic events and substantial intensity variability, such as precipitation. The methodology involves establishing a probability threshold to adapt the occurrence of precipitation events and utilizes a quantile mapping method based on a non-stationary theoretical and parametric distribution to adjust biases associated with precipitation intensity. The NS method is applied to daily precipitation projections from seven regional climatic models under the RCP 8.5 scenario spanning 2006–2100, alongside historical concurrent data from 1970 to 2005. The present method is compared to the widely used quantile delta mapping approach (QDM), revealing significant differences in performance related to the distribution of precipitation events throughout the year and the behaviour of mean and extreme intensity values. Both approaches show reliable performance, with root mean square errors between the empirical distribution functions of the corrected hindcast time series and the observations being lower than or close to 1 mm for the percentiles smaller than 50th. The error increases with the percentiles, particularly at stations located at higher altitudes. In these locations, QDM doubles or triples the error obtained with the non-stationary approach for percentiles higher than 75th, reaching up to 34 mm. The proposed methodology demonstrates promising potential in reducing uncertainties associated with systematic errors in inter-annual precipitation variability. It is a first step to jointly apply a similar approach to all involved driving variables. Such methods are key to assessing hydrological responses and associated impacts in semi-arid mountainous basins all around the Globe.