The complexity of boron species in aqueous solution becomes the bottleneck in development of boron resources from brine. The concentrations of different boron species and transformation mechanism among boron species are mainly affected with the total boron concentration, temperatures, pH of the solution, and cations in the solution. The change trends for concentration of total boron in the complicated aqueous solution can be obtained with the phase equilibrium results, but the accurate concentration of boron species in the solution cannot be obtained. The metastable zone width (MZW) of borates and H3BO3 changed with the concentrations of impurity ions, the stirring rate, cooling rate, and original temperature for the saturated solution. The Pitzer model is the effective method to calculate the concentration of boron species. More thermodynamic data at temperatures below and above 298.15 K are needed for model parameterization. The main problems existed in the systems containing boron species were summarized, and the future development direction is also discussed. The progresses on boron species in the aqueous solution can provide fundamental data in the comprehensive exploitation of the brine boron resource.