The solubility of tebuconazole in nine organic solvents including methanol, ethanol, isopropanol, n-propanol, ethyl acetate, toluene, acetone, 2-butanone and acetonitrile was determined experimentally by using the isothermal saturation method over a temperature range from (278.15 to 313.15) K under 101.2 kPa. For the temperature range investigated, the mole fraction solubility of tebuconazole in the solvents increased with a rise of temperature. On the whole, they obeyed the following order from high to low in different solvents: 2-butanone > acetone > (ethyl acetate, toluene) > methanol > ethanol > isopropanol > n-propanol > acetonitrile. The acquired solubility data of tebuconazole in the studied solvents were correlated by using the modified Apelblat equation, lambda h equation, Wilson and NRTL models. The maximum values of root-mean-square deviation (RMSD) Was 12.35 x 10(-4), and the relative average deviation (RAD), 1.56%. Generally, the four thermodynamic models could all be employed to describe the solubility behaviour of tebuconazole in these solvents, and the modified Apelblat equation correlated the experimental data best according to the result of Akaike Information Criterion (AIC) analysis. Furthermore, the mixing Gibbs energy, mixing enthalpy, Mixing entropy, activity coefficient at infinitesimal concentration (gamma(infinity)(1)) and reduced excess enthalpy (H-1(E,infinity)) were computed. The solubility obtained and thermodynamic studies should be very helpful for optimizing the purification process of tebuconazole. (C) 2016 Elsevier Ltd.