A numerical study of three-dimensional solute transportat fracture intersections by using a particle tracking technique is presented.Two models of orthogonal fracture intersection are considered, namely, twoparallel-walled channels and two rough-walled Gaussian fractures. The fluidvelocity is calculated by solving the three-dimensional Stokes equation withno-slip boundary condition at the solid wall. Examples of individual trajectoriesof particles are first given in order to illustrate the main features of thephenomenon. Solute mass partitioning between outgoing fracture branches isconsidered for various transport regimes, characterized by the local P¿cletnumber, and for various ratios of the flow rates in the intersecting channels.Generally speaking, it can be said that at dominant diffusion the influenceof the flow rates ratio is weak, while it is important in the opposite situation.Validity of the classical models of solute mixing, stream tube routing, andperfect mixing is analyzed by comparing their predictions with the numericaldata. Preliminary recommendations are made for the use of these results inlarge-scale modeling.