ABSTRACT The assessment of extreme hydrological events relies on the assumption that time series are independent and identically distributed. However, empirical evidence contradicts this assumption, indicating the presence of nonstationarity, including trends, step changes, or both, in the river records of the Argentinean Central Andes. This study proposes the use of generalized additive models in location, scale and shape (GAMLSS) to model the annual maximum and minimum flow of seven rivers. Firstly, the presence of trends and step changes was assessed. Subsequently, the models were formulated considering time and various climatic indices as covariates. The findings reveal declining trends and negative jumps over the past 50 to 70 years. The location parameters of the models exhibit linear or smoothing dependence with time or climate indices covariates. Different modes of climate variability are associated with hydrological extremes. The nonstationary models provide a novel complementary framework for water resource planning in the region.