Aerosol Optical Thickness (AOT) is one of the critical factors for global atmospheric conditions, climate change, and air pollution. AOT has been exposed as a major component of air pollution in Bangladesh. This paper aims to map the seasonal distribution of AOT from 2002-2022 and to explore the internal relationship between AOT and ten air pollutants using remote sensing and machine learning tools. These ten air pollutants are Particulate matter (PM2.5), Methane (CH4), Carbon monoxide (CO), Nitrogen dioxide (NO2), Formaldehyde (HCHO), Ozone (O3), Sulfur dioxide (SO2), Aerosol Particulate Radius (APR), Nitrogen oxide (NOx) and Black carbon (BC). The results show that the concentrations of AOT were higher in December-January-February (mean value 0.50) and March-April-May (mean value 0.50) seasons, mostly in the central, western, and southern parts of Dhaka, Narayanganj, and Munshiganj districts. AOT was a bit less in June-July-August (mean value 0.33) and September-October-November (mean value 0.37). This paper also revealed that the AOT was correlated positively with PM2.5 (0.60), CH4 (0.80), NO (0.76), and BC (0.83) while correlated negatively with CO (-0.66), HCHO (-0.16), SO2 (-0.41), APR (-0.48), and NOx (-0.20). From the machine learning, the Rational quadratic GPR (RME-0.0024, MAE-0.0015, R2-0.96), Matern 5/2 GPR (RMSE-0.0023, MAE-0.0015, R2-0.96), and Squared Exponential GPR (RMSE-0.0015, MAE-0.0015, R2-0.96) were found good classifiers to predict AOT. UN agencies, government line departments, and local and regional development councils for air pollution mitigation and long-term protective measures may use the paper's key results.