We investigate the impact of injection mode (flux and resident injection) and heterogeneity in hydraulic properties on dispersion of advecting particles in two-dimensional discrete fracture network models, using a Monte Carlo method. We find that the injection mode has a significant effect on dispersion: The resident injection mode exhibits anomalous features of transport whereas the flux injection mode tends to Gaussian transport; this observation is easily understood by considering phase diagrams where a limited number of particles entering low velocity fractures greatly increase macrodispersion. In spite of a sizeable portion of negative longitudinal velocities, it is shown that multiple crossings are negligible when quantifying longitudinal macrodispersion. A simple probabilistic expression of particle mass balance is shown to predict well the spatial distribution of advecting particles.