The dispersion and mixing of passive scalars released from two concentrated sources into open-channel and homogeneous isotropic turbulent flows are studied using direct numerical simulation (DNS). The simulations are conducted using two fully-parallelized in-house codes developed using the FORTRAN 90/95 programming language. A comparative study has been conducted to investigate the effects of the source separation distance, Reynolds number, relative length scales of the plume and turbulent flow, and source elevation on the dispersion and mixing of two plumes. For both flow configurations, four distinct stages in the downwind development of the cross correlation between the fluctuating concentration fields have been identified which feature zero, destructive and constructive interferences and a complete mixing state. Differences between the exceedance probability of concentrations for the single and total plumes are highlighted and analyzed, and the effects of destructive and constructive interference on the exceedance probabilities for the total plume are used to explain these differences. It is found that the relationship between the third- and fourth-order concentration moments and the second-order concentration moment can be well predicted using a clipped-gamma model. This leads to an interesting conclusion that all the higher-order (third-order and above) moments of the total concentration can be inferred from a knowledge of only the first- and second-order concentration moments of each single plume and of the cross correlation coefficient. From a spectral analysis, it is observed that there exists a range of `leading scales' at which the rate of turbulent mixing of the two plumes becomes the most efficient and the coherency spectrum of the plumes approaches the asymptotic value of unity quicker than at any other scales. / October 2016