Can the ensemble smoother with multiple data assimilation be used to predict discharge in an Alpine karst aquifer? The answer is yes, at least, for the Bossea aquifer studied. The ensemble smoother is used to fit a unit hydrograph simultaneously with other parameters in a hydrologic model, such as base flow, infiltration coefficient, or snow melting contribution. The fitting uses observed discharge flow rates, daily precipitations, and temperatures to define the model's parameters. The data assimilation approach gives excellent results for fitting individual events. After the analysis of 27 such events, two average models are defined to be used to predict flow discharge from precipitation and temperature, one model for prediction during spring (when snow melting has an impact) and another one during autumn, yielding acceptable results, particularly for the fall rainfall events. The lesser performance for the spring events may indicate that the snow melting approximation needs to be revised. The results also show that the parameterization of the infiltration 1Pansa, Butera, Gómez-Hernández, Vigna coefficient needs further exploration. Overall, the main conclusion is that the ensemble smoother could be used to define a characteristic "signature" of a karst aquifer to be used in forecast analyses.The reasons for using the ensemble smoother instead of other stochastic approaches are that it is easy to use and explain and provides an estimation of the uncertainty about the predictions.