The isotope effects of gadolinium in Gd-EDTA ligand exchange system were studied by means of ion exchange chromatography. The separation coefficients of gadolinium isotopes, epsilon, and the local enrichment factors, beta, were calculated from the observed isotopic ratios at the front and rear boundaries of the gadolinium adsorption band. Clear mass independent anomalies were observed in the isotope effects of 155Gd and 157Gd. The relation between the isotope effects of gadolinium isotopes, studied by the three-isotope plot and the separation coefficient methods, and the mass of gadolinium isotopes was found to be related to the change in the mean square radius of the nuclear charge distribution parameter, (r2), of these isotopes, which suggests that the nucleus shape and size highly affect the gadolinium isotope effects in chemical exchange reactions.