any way to make correct prediction using SARIMA with multiple dependent variables in python?

Question

I have 14 years of monthly weather data , I intend to make a forecast on pollution_index for next 3 years. To do so, I used SARIMA model from