The mathematical model for the calculation of theoxidation-reduction species distribution in a flowing glass melt under non-isothermalconditions was formulated. The results of calculations show that the initialoxidation-reduction state of a glass melt influences both the rate of bubbleremoving from the melt and the bubble composition development. The very high orvery low values of the redox state accelerate the bubble removing process inglasses refined by sulphates while only high redox values are favourable inglasses refined by oxides. The concentration ratios CO2/N2as well as N2/Ar in analysed bubbles, used for the identification ofbubble sources in a glass melting space, generally grow when the bubble removingprocess is accelerated.