Analysis of low river flows provides important information for effective management of water resources in a region. Despite the critical importance of understanding low flow dynamics, there is a gap in the literature regarding the use of non-stationary models to analyze low flow data under climate change in Turkey. In this research, low flow series from 80 measuring stations in Turkey are investigated by employing both stationary and non-stationary models based on the Generalized Additive Models for Location, Scale and Shape (GAMLSS). For constructing non-stationary models, 31 explanatory variables consisting of time, precipitation, temperature and atmospheric oscillation indices were used to model the parameters of the chosen distributions. The results show that stationary models are more successful at 7 stations, while non-stationary models are more successful at 73 stations. Comparisons between non-stationary models showed that for most stations, the best performing models were non-stationary models with annual precipitation as covariates. In addition, successful results were obtained when Western Mediterranean Oscillation and North Atlantic Oscillation indices were used as explanatory variables. Additionally, this study investigated 20 and 50-year return levels by fitting the non-stationary frequency distribution models for low flows over historical and projection periods under SSP2-4.5 and SSP5-8.5 climate scenarios. GAMLSS incorporated annual total precipitation, which is the most effective explanatory variable for low flows, as a covariate, and thus changes in low flows were analyzed. The results show that decreases are expected in low flows, except for the stations in the upper Euphrates basin compared to the historical period.