Aqueous speciation of rare earth elements (REE) controls their mobilization, fractionation, and enrichment in the natural waters. Geochemical modeling of their speciation is key to improve our understanding of the formation of economic mineral deposits, for developing mineral separation and mine tailing recovery technologies, and for characterizing the geochemistry of thermal water. However, our ability to predict the fate of REE in a wide pH and temperature range is limited by the scarcity of thermodynamic data for the REE hydroxyl complexes. In part I of this study (H. J. Han and A. P. Gysi,