Total dissolved solids (TDS) concentration, as an essential variable in evaluating the quality of drinking and agricultural water, represents water body salinity. In a river system, high TDS concentration has negative impacts on human health and crops consuming water. Since anions and cations affect the TDS value, identifying a relationship between these variables and TDS can help predict and monitor river quality. This research investigates the Gaussian process regression (GPR) model capabilities as a data-driven model to capture the relationship and estimates the TDS value in the Tajan River watershed in Northern Iran. Monthly anions and cations measured over 16 years including bicarbonate (HCO3-), carbonate (CO32-), sulfate (SO42-), chloride (Cl-), calcium (Ca2+), magnesium (Mg2+), sodium (Na+), and potassium (K+) are considered as the predictor variables. Five GPR kernel functions are applied in the modeling, and their efficiency is evaluated using four statistics: the coefficient of determination (R-2), Mean absolute error (MAE), Mean squared error (MSE), and Nash-Sutcliffe efficiency (NSE). Also, the performance of the proposed method is assessed by comparing it to the Artificial neural network (ANN) model, as an efficient and popular prediction model. The results reveal that the GPR model with the rational quadratic kernel function performed better in terms of performance criteria (R-2 = 0.9836).