The Saharan heat low (SHL) is a key component of the West African climate system and an important 10 driver of the West African Monsoon across a range of timescales of variability. The physical mechanisms driving the variability in the SHL remain uncertain, although water vapour has been implicated as of primary importance. Here, we quantify the independent effects of variability in dust and water vapour on the radiation budget and atmospheric heating of the region using a radiative transfer model configured with observational input data from the Fennec field campaign at the location of Bordj Badji Mokhtar (BBM) in southern Algeria (0.9E, 21.4N), close 15 to the SHL core, for June 2011. Overall, we find dust aerosol and water vapour to be of similar importance in driving variability in the top of atmosphere (TOA) radiation budget and therefore the column integrated heating over the SHL (~7 W m-2 per standard deviation of dust AOD). As such we infer that SHL intensity is likely to be similarly enhanced by the effects of dust and water vapour surge events. However, the details of the processes differ. Dust generates substantial radiative cooling at the surface (~11 W m-2 per standard deviation of dust AOD), 20 presumably leading to reduced sensible heat flux into the boundary layer, which is more than compensated by direct radiative heating from SW absorption by dust in the dusty boundary layer. In contrast water vapour invokes a radiative warming at the surface of ~6 W m-2 per standard deviation of column integrated water vapour in Kg m-2. Net effects involve a pronounced net atmospheric radiative convergence with heating rates on average of 0.5 K day-1 and up to 6 K day-1 during synoptic/meso-scale dust events from monsoon surges and convective cold 25 pool outflows (‘haboobs’). On this basis we make inferences on the processes driving variability in the SHL associated with radiative and advective heating/cooling. Depending on the synoptic context over the region processes driving variability involve both independent effects of water vapour and dust and compensating events in which dust and water vapour are co-varying. Forecast models typically have biases of up to 2 kg m-2 in column integrated water vapour (equivalent to a change in 2.6 W m-2 TOA net flux) and typically lack variability in dust, 30 and so are expected to poorly represent these couplings. An improved representation dust and water vapour and quantification of associated radiative impact is thus imperative in quest for the answer to what remains to be uncertain related with the climate system of the SHL region.