The ion-specific surface tension of aqueous electrolyte solutions is of fundamental importance in physical chemistry, solution chemistry, electrochemistry, and biochemistry, yet it remains a challenge to be qualitatively predicted. In this work, an analytical theory of ion-specific surface tension is developed. By modeling the solution to a restricted primitive model (RPM), the surface tension increment of the electrolyte solution reduces to the surface tension of the RPM, which is further determined analytically by using the integral equation theory and the morphological thermodynamics theory. According to our formula, the surface tension increment of the electrolyte solution consists of a dominant and positive hard sphere contribution, which depends on the salt concentration, and a secondary electrostatic contribution, which depends on the inverse Debye length and the salt concentration. Our theory is applied to 1:1, 2:2, 1:2, 2:1, 3:1, and 3:2 electrolyte solutions with typical salt concentrations up to several mol/L and compared with experimental data. Without introducing any adjustable parameters, our theory leads to a good prediction of the surface tension increment of more than 50 kinds of aqueous solutions. Such a good agreement demonstrates the great potential of our theory for a fundamental understanding of specific ion effects in a variety of electrolyte solutions.