Biomagnification of hydrophobic organic compounds (HOCs) increases the eco-environmental risks they pose. Here, we gained mechanistic insights into biomagnification of deuterated polycyclic aromatic hydrocarbons (PAHs- d10) in zebrafish with carefully controlled water (ng L-1) by a passive dosing method and dietary exposures using pre-exposed Daphnia magna and fish food. A new bioaccumulation kinetic model for fish was established to take into account discrete dietary uptake, while the frequently used model regards dietary uptake as a continuous process. We found that when freely dissolved concentrations of the PAHs- d10 were constant in water, the intake amount of the PAHs- d10 played an important role in affecting their steady-state concentrations in zebrafish, and there was a peak concentration in zebrafish after each dietary uptake. Moreover, considering the randomness of predation, the Monte Carlo simulation results showed that the probabilities of biomagnification of the PAHs- d10 in zebrafish increased with their dietary uptake amount and frequency. This study indicates that in addition to the well-known lipid-water partitioning, the bioaccumulation of HOCs in fish is also a discontinuous kinetic process caused by the fluctuation of HOC concentration in the gastrointestinal tract as a result of the discrete food ingestion. The discontinuity and randomness of dietary uptake can partly explain the differences among aquatic ecosystems with respect to biomagnification for species at similar trophic levels and provides new insight for future analysis of experimental and field bioaccumulation data for fish.