Evidence on joint association of a phthalate mixture with thyroid function among children and its underlying mechanism is largely unknown. We aimed to explore the associations of 10 urinary phthalate metabolites (mPAEs), either as individuals or as a mixture, with thyroid function indicators [free thyroxine, free triiodothyronine (FT3), and thyroid-stimulating hormone (TSH)] in 144 children aged 4-12 years with up to 3 repeated visits across 3 seasons. Significant and positive associations were observed for mono-(2-ethylhexyl) phthalate (MEHP), mono-iso-butyl phthalate (MiBP), and mono-n-butyl phthalate (MnBP) with TSH, as well as monobenzyl phthalate (MBzP) with FT3 in dose-response manners. The relationship between MEHP and TSH remained robust in multiple-phthalate models. Bayesian kernel machine regression (BKMR) models revealed overall linear associations of the 10 mPAE mixture with higher TSH and FT3 levels, and MEHP and MBzP were major contributors. Meanwhile, MEHP, MiBP, and MnBP were linked to the elevation of multiple cytokines including CCL 27, CCL3, CXCL1, and IL-16. Among them, IL-16 mediated the relationships of MEHP and MiBP with TSH, and the mediated proportions were 24.16% and 24.27%, respectively. Our findings suggested that mPAEs dominated by MEHP were dose-responsively associated with elevated TSH among healthy children and mediated by IL-16.