Three dimensional interstellar extinction maps provide a powerful tool for stellar population analysis. However, until now, these 3D maps were rather limited by sensitivity and spatial resolution. We use data from the VISTA Variables in the Via Lactea survey together with the Besancon stellar population synthesis model of the Galaxy to determine interstellar extinction as a function of distance in the Galactic bulge covering -10<l<10 and -10<b<5. We adopted a recently developed method to calculate the colour excess. First we constructed the H-Ks vs. Ks and J-Ks vs. Ks colour-magnitude diagrams based on the VVV catalogues that matched 2MASS. Then, based on the temperature-colour relation for M giants and the distance-colour relations, we derived the extinction as a function of distance. The observed colours were shifted to match the intrinsic colours in the Besancon model as a function of distance iteratively. This created an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. We present a 3D extinction map that covers the whole VVV area with a resolution of 6'x6', for J-Ks and H-Ks using distance bins of 0.5-1.0kpc. The high resolution and depth of the photometry allows us to derive extinction maps for a range of distances up to 10kpc and up to 30 magnitudes of extinction in AV (3.0mag in AKs). Integrated maps show the same dust features and consistent values as other 2D maps. We discuss the spatial distribution of dust features in the line of sight, which suggests that there is much material in front of the Galactic bar, specifically between 5-7kpc. We compare our dust extinction map with the high-resolution ^12^CO maps (NANTEN2) towards the Galactic bulge, where we find a good correlation between ^12^CO and A_V_. We determine the X factor by combining the CO map and our dust extinction map. Our derived average value X=2.5+/-0.47x10^20^/(cm^2^.K.km/s) is consistent with the canonical value of the Milky Way. The X-factor decreases with increasing extinction.

Cone search capability for table J/A+A/566/A120/table1jk (3D Interstellar extinction as a function of distance and E(J-K))

Cone search capability for table J/A+A/566/A120/table1hk (3D Interstellar extinction as a function of distance and E(H-K))