Long-term aerosol climatology is derived using Terra MODIS (Collection 6) enhanced Deep Blue (DB) AOD retrieval algorithm to investigate decadal trend (2006–2015) in columnar aerosol loading, future scenarios and potential source fields over the Indo-Gangetic Plain (IGP), South Asia. Satellite based aerosol climatology was analyzed in two contexts: for the entire IGP considering area weighted mean AOD and for nine individual stations located at upper (Karachi, Multan, Lahore), central (Delhi, Kanpur, Varanasi, Patna) and lower IGP (Kolkata, Dhaka). A comparatively high aerosol loading (AOD: 0.50 ± 0.25) was evident over IGP with a statistically insignificant increasing trend of 0.002 year−1. Analysis highlights the existing spatial and temporal gradients in aerosol loading with stations over central IGP like Varanasi (decadal mean AOD±SD; 0.67 ± 0.28) and Patna (0.65 ± 0.30) exhibit the highest AOD, followed by stations over lower IGP (Kolkata: 0.58 ± 0.21; Dhaka: 0.60 ± 0.24), with a statistically significant increasing trend (0.0174–0.0206 year−1). In contrast, stations over upper IGP reveal a comparatively low aerosol loading, having an insignificant increasing trend. Variation in AOD across IGP is found to be mainly influenced by seasonality and topography. A distinct “aerosol pool” region over eastern part of Ganges plain is identified, where meteorology, topography, and aerosol sources favor the persistence of airborne particulates. A strong seasonality in aerosol loading and types is also witnessed, with high AOD and dominance of fine particulates over central to lower IGP, especially during post-monsoon and winter. The time series analyses by autoregressive integrated moving average (ARIMA) indicate contrasting patterns in randomness of AOD over individual stations with better performance especially over central IGP. Concentration weighted trajectory analyses identify the crucial contributions of western dry regions and partial contributions from central Highlands and north-eastern India, in regulating AOD over stations across IGP. Although our analyses provide some attributes to the observed changes in aerosol loading, we conclude that the spatial and temporal pattern of aerosol properties is highly complex and dynamic over IGP, and require further investigation in order to reduce uncertainty in aerosol-climate model.