There is a substantial lack of tacrolimus pharmacokinetic information in pediatric hematopoietic stem cell transplant (HSCT) patients. This study aimed to develop population pharmacokinetics (PopPK) of tacrolimus in pediatric HSCT patients and to devise model-guided dosage regimens. A retrospective analysis was performed on 86 pediatric HSCT patients who received tacrolimus intravenously or orally. A total of 578 tacrolimus trough concentrations (C) were available for pharmacokinetic analysis using a non-linear mixed-effects modeling method. Demographic and clinical data were included and assessed as covariates the stepwise method. Bayesian estimators were used to devise pediatric dosage regimens that targeted C of 5-15 ng mL. A one-compartment model with first-order absorption adequately described the tacrolimus pharmacokinetics. Clearance (CL), volume of distribution (V), and typical bioavailability (F) in this study were estimated to be 2.42 L h (10.84%), 79.6 L (16.51%), and 19% (13.01%), respectively. Body weight, hematocrit, post-transplantation days, and caspofungin and azoles concomitant therapy were considered significant covariates for tacrolimus CL. Hematocrit had a significant impact on the V of tacrolimus. In the subgroup cohort of children ( = 24) with genotype, the clearance was 1.38-fold higher in CYP3A5 expressers than in non-expressers. Simulation indicated that the initial dosage optimation of tacrolimus for intravenous and oral administration was recommended as 0.025 and 0.1 mg kg d (q12h), respectively. A PopPK model for tacrolimus in pediatric HSCT patients was developed, showing good predictive performance. Model-devised dosage regimens with trough tacrolimus concentrations provide a practical strategy for achieving the therapeutic range.