Интервал в метке оси при использовании expression(paste(…))

Рассмотрим следующий пример:

plot(c(2,4,6)~c(1,2,3),xlab="x",
ylab=expression(paste('flux (g ',CO[2]~m^{-2}~h^{-1},')')))

Очевидно, я хочу полный пробел между "g" и "CO", но почему-то получаю меньший (с некоторыми метками даже нулевой ) пробел в метке графика.

Проблема еще более очевидна, если я сделаю это так:

plot(c(2,4,6)~c(1,2,3),xlab="x",
ylab=expression(paste('flux (g C',O[2]~m^{-2}~h^{-1},')')))

Я делаю что-то не так? Есть ли способ исправить интервал или даже лучший способ создать метки с большим количеством подстрочных/надстрочных индексов и греческих букв?