The main objective of this study was to develop simple models for the prediction of bromate formation in ozonated bottled waters, using rapidly and practically measurable raw water quality and/or operational parameters. A total of 6 multi-linear regression ( MLR) with or without principal component analysis (PCA) and 2 artificial neural networks ( ANN) models with multilayer perceptron architecture were developed for the prediction of bromate formation. PCA was employed to better identify relations between variables and reduce the number of variables. Experimental data used in modeling was provided from the ozonation of samples from 5 groundwater sources at various applied ozone dose and contact time. MLR models#1 and #2 well-predicted bromate formation although correlations ( i.e., the signs of regression constants) among pH ( as input variable) and bromate concentrations did not agree with the chemistry. MLR model#6, containing practical input parameters that are measured on-line in full-scale treatment plants, adequately predicted bromate formation and agreed with the chemistry, although fewer input parameters were used compared to MLR#1 and #2. Although both of the ANN models exhibited high regression coefficients ( R 2) ( 0.97 for both) ANN#1 was found to provide better prediction of bromate formation based on mean square error (MSE) values. However, since ANN#2 included easily measurable input parameters it may be practically used by water companies employing ozonation. Results overall indicated that ANN models have stronger prediction capabilities of bromate formation than MLR models. ANN modeling appears to be a strong tool in situations where the relations between variables are nonlinear, interactive and complex, as in the bromate formation by ozonation.