A novel method for the trace analysis of 17beta-estradiol (E2) in river water sample was developed, which involved stir bar sorptive extraction (SBSE) with in situ acylation (first derivatization) and thermal desorption (TD) with quartz wool assisted (QWA) in tube silylation (second derivatization), followed by gas chromatography-mass spectrometry (GC-MS), and is called the "dual derivatization method." The optimum conditions for SBSE with in situ acylation, such as the volume of acetic acid anhydride and the extraction time, were investigated. In addition, the optimum conditions for TD with QWA in tube silylation, such as the volume of N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA) and the TD temperature and hold time, were investigated as well. The detection limit (S/N = 3) and the quantitation limit (S/N>10) of E2 in the river water sample were 0.5 and 2 pg ml(-1) (ppt), respectively, by the dual derivatization method. In addition, the detection limit was 0.1 pg ml(-1) by using dual derivatization method with multi-shot mode. The calibration curve for E2 was linear in the range of 0.002-10 ng ml(-1) with correlation coefficients >0.999. The average recoveries of E2 (n = 6) at the concentrations of 0.05 and 1.0 ng ml(-1) from the river water sample were 93.1 (RSD: 1.4%) and 98.4% (RSD: 0.8%), respectively, with correction using the added surrogate standard, 17beta-estradiol-(13)C(4). This simple, accurate, sensitive and selective analytical method may be applicable to the determination of trace amounts of E2 in water samples.