Large diffusion-weighted brain MRI (dMRI) studies in neonates are crucial for developmental neuroscience. Our aim was to investigate the utility of ComBat, and empirical Bayes tool for multisite harmonization, for removing site effects from white matter (WM) dMRI measures in healthy infants born 37-42+6 weeks from the Theirworld Edinburgh Birth Cohort (n=86) and Developing Human Connectome Project (n=287). Skeletonized fractional anisotropy (FA), mean, axial and radial diffusivity (MD, AD, RD) maps were harmonized. The differences between voxel-wise metrics, skeleton means and histogram widths (5th-95th percentile) were assessed before and after harmonization, as well as variance associated with gestational age at birth. Before harmonization, large cohort differences were observed in all measures. Harmonization removed all voxel-wise differences from MD maps and all metric means and histogram widths, however small voxel-wise differences (<1.5% of voxels) remained in FA, AD and RD. We detected significant relationships between GA at birth and all metrics. When comparing single site and multi-site harmonized datasets of equal sample sizes, harmonized data resulted in smaller standardized regression coefficients. ComBat will enable unprecedented sample sizes in developmental neuroscience, offering new horizons for biomarker discovery and validation, understanding typical and atypical brain development, and assessment of neuroprotective therapies.