A hybrid method for the simulation of radionuclide contaminant plumes in heterogeneous unsaturated formations