A toxicokinetic model is proposed to predict the time evolution of malathion and its metabolites, mono- and dicarboxylic acids (MCA, DCA) and phosphoric derivatives (dimethyl dithiophosphate [DMDTP], dimethyl thiophosphate [DMTP], and dimethyl phosphate [DMP]) in the human body and excreta, under a variety of exposure routes and scenarios. The biological determinants of the kinetics were established from published data on the in vivo time profiles of malathion and its metabolites in the blood and urine of human volunteers exposed by intravenous, oral, or dermal routes. In the model, body and excreta compartments were used to represent the time varying amounts of each of the following: malathion, MCA, DCA, DMDTP, DMTP, and DMP. The dynamic of intercompartment exchanges was described mathematically by a differential equation system that ensured conservation of mass at all times. The model parameters were determined by statistically adjusting the explicit solution of the differential equations to the experimental human data. Simulations provide a close approximation to kinetic data available in the published literature. When simulating a dermal exposure to malathion, the main route of entry for workers, the model predicts that it takes an average of 11.8 h to recover half of the absorbed dose of malathion eventually excreted in urine as metabolites, compared to 3.2 h following an intravenous injection and 4.0 h after oral administration. This shows that following a dermal exposure, the absorption rate governs the urinary excretion rate of malathion metabolites because the dermal absorption rate is much slower than biotransformation and renal clearance processes. The model served to establish biological reference values for malathion metabolites in urine since it allows links to be made between the absorbed dose of malathion and the time course of cumulative amounts of metabolites excreted in urine. From the no-observed-effect level (NOEL) of 0.61 micromol/kg/day derived from the data of Moeller and Rider (1962), the model predicts corresponding biological reference values for MCA, DCA, and phosphoric derivatives of 44, 13, and 62 nmol/kg, respectively, in 24-h urine samples. The latter were used to assess the health risk of workers exposed to malathion in botanical greenhouses, starting from urinary measurements of MCA and DCA metabolites.