Model averaging provides the kinetic model taking into account model uncertainty.
ETBE is formed faster than TAEE in the simultaneous etherification.
Reactions mechanism is deduced from the proposed kinetic model.
Active sites are mainly occupied by adsorbed ethanol, ETBE and TAEE.
Two active sites participate in etherification reactions, and one in isomerization.