Abstract. Temperature fields in evaporating liquids are measured by simultaneous use of Background Oriented Schlieren (BOS) technique for the side view and IR thermal imaging for the surface distribution. Good agreement between the two methods is obtained with typical measurement error less than 0.1 K. Two configurations of surface layer are observed: thermocapillary convection state with moving liquid surface and small thermal cells, associated with Marangoni convection, and "cool skin" with negligible velocity at the surface, larger cells and dramatic increase of velocity within 0.1 mm layer beneath the surface. These configurations are shown to be formed in various liquids (water with various degrees of purification, ethanol, butanol, decane, kerosene, glycerine) depending rather on initial conditions and ambient parameters than on the liquid. Water, which has been considered as the liquid without observable Marangoni convection, actually can exhibit both kinds of behavior during the same experimental run. Evaporation is also studied by means of numerical simulations. Separate problems in air and liquid are considered, with thermal imaging data of surface temperature making the separation possible. It is shown that evaporation rate can be predicted by numerical simulation of the air side with appropriate boundary conditions. Comparison is made with known empirical correlations for Sherwood-Rayleigh relationship. Numerical simulations of water-side problem reveal the issue of velocity boundary conditions at the free surface, determining the structure of surface layer. Flow field similar to observed in the experiments is obtained with special boundary conditions of third kind, presenting a combination of no-slip and surface tension boundary conditions.