Explaining and predicting the distribution of ions at the air water interface has been a central challenge of physical chemistry for nearly a century. In essence, the problem amounts to calculating the change in the solvation energy of an ion as it approaches the interface. Here, we generalize our recently developed model of ionic solvation energies to calculate this interaction. The change in the Born energy as well as the static polarization response of the ion is included by using the conductor-like screening model (COSMO), which treats the ions quantum mechanically. Approximate expressions for the dispersion repulsion, cavity attraction, and surface potential contributions are also included. This model reproduces the surface tensions of electrolyte solutions and is consistent with ab initio molecular dynamics (MD) simulation. The model provides clear physical insight into iodide's adsorption. Unlike alternative models, no parameters are deliberately adjusted to reproduce surface tensions, and all of the important contributions to the interactions are included. Solving this problem has important direct implications for atmospheric chemistry and bubble properties. It also has important indirect implications for the more complex interactions of ions with protein and mineral surfaces. These play a fundamental role in a vast number of biological and industrial processes. The model is conceptually simple and has low computational demand, which facilitates its extension to these important applications.