We report magnetic birefringence measurements up to high fields (17.5 T) of dilute aqueous Suspensions of rod-like cellulose nanocrystals With well characterized distributions of lengths, widths and thicknesses. We compare these data with three Models, one with colinear (1), one with perpendicular cylindrically symmetric tensors for diamagnetic susceptibility and refractive index (2) and one with biaxial diamagnetic anisotropy (3). We find that taking into account polydispersities of length, width, and thickness is essential for accurate fitting and that model 1 is the most appropriate, presumably because of the twisting of the suspeiided nanocrystal along their long axis. The best-fitted susceptibility anisotropy was Delta chi(z(xy)) = chi(zz)-(chi(xx)+chi(yy)),)/2 = -2.44 x 10(-6) when considering only the crystalline core of nanocrystals and, more appropriately, Azz(xy) = -0.05 x 10 when including crystalline core and skin. The latter -value is slightly higher than Delta chi(z(xy)) = -0.68(5) x 10(-6) deduced from estimations using Pascal's additivitylaw. The specific birefringence of the nanocrystals in water was found to be delta n(0) = +0.120(2), which is well accounted for by the iutrinsic birefringence of crystalline cellulose (delta n(0)(intr) = n parallel to-n(1) = +0.0744) and the birefringence arising from the slender shape of nanocrystals.