The determination of trace rare-earth elements (REEs) can be used for the assessment of environmental pollution, and is of great significance to the study of toxicity and toxicology in animals and plants. N, N, N', N'-tetraoctyl diglycolamide (TODGA) is an environmental friendly extractant that is highly selective to REEs. In this study, an analytical method was developed for the simultaneous determination of 16 trace REEs in simulated water samples by inductively coupled plasma optical emission spectroscopy (ICP-OES). With this method, TODGA was used as the extractant to perform the liquid-liquid extraction (LLE) sample pretreatment procedure. All 16 REEs were extracted from a 3 M nitric acid medium into an organic phase by a 0.025 M TODGA petroleum ether solution. A 0.03 M ethylenediaminetetraacetic acid disodium salt (EDTA) solution was used for back-extraction to strip the REEs from the organic phase into the aqueous phase. The aqueous phase was concentrated using a vacuum rotary evaporator and the concentration of the 16 REEs was detected by ICP-OES. Under the optimum experimental conditions, the limits of detection (3σ, n = 7) for the REEs ranged from 0.0405 ng mL-1 (Nd) to 0.5038 ng mL-1 (Ho). The relative standard deviations (c = 100 ng mL-1, n = 7) were from 0.5% (Eu) to 4.0% (Tm) with a linear range of 4-1000 ng mL-1 (R2 > 0.999). The recoveries of 16 REEs ranged from 95% to 106%. The LLE-ICP-OES method established in this study has the advantages of simple operation, low detection limits, fast analysis speed and the ability to simultaneously determine 16 REEs, thereby acting as a viable alternative for the simultaneous detection of trace amounts of REEs in water samples.