Abstract. A compilation of data from several cruises between 1998 and 2013 was used to derive polynomial fits that estimate total alkalinity (A_T) and total dissolved inorganic carbon (C_T) from measurements of salinity and temperature in the Mediterranean Sea surface waters. The optimal equations were chosen based on the 10-fold cross-validation results and revealed that second- and third-order polynomials fit the A_T and C_T data respectively. The A_T surface fit yielded a root mean square error (RMSE) of ± 10.6 μmol kg⁻¹, and salinity and temperature contribute to 96 % of the variability. Furthermore, we present the first annual mean C_T parameterization for the Mediterranean Sea surface waters with a RMSE of ± 14.3 μmol kg⁻¹. Excluding the marginal seas of the Adriatic and the Aegean, these equations can be used to estimate A_T and C_T in case of the lack of measurements. The identified empirical equations were applied on the 0.25° climatologies of temperature and salinity, available from the World Ocean Atlas 2013. The 7-year averages (2005–2012) showed that A_T and C_T have similar patterns with an increasing eastward gradient. The variability is influenced by the inflow of cold Atlantic waters through the Strait of Gibraltar and by the oligotrophic and thermohaline gradient that characterize the Mediterranean Sea. The summer–winter seasonality was also mapped and showed different patterns for A_T and C_T. During the winter, the A_T and C_T concentrations were higher in the western than in the eastern basin. The opposite was observed in the summer where the eastern basin was marked by higher A_T and C_T concentrations than in winter. The strong evaporation that takes place in this season along with the ultra-oligotrophy of the eastern basin determines the increase of both A_T and C_T concentrations.