Improving approximate matrix factorizations for implicit time integration in air pollution modelling