Focus in quality assessment of iron ores is the content of total iron (TFe). Laser-induced breakdown spectroscopy (LIBS) technology possesses the merits of rapid, in-situ, in real-time multi-element analysis for iron ores, but its application for quantitative TFe content is subject to interference of the iron matrix effect and the lack of suitable data mining tools. Herein, a new method of LIBS-based variable importance-back propagation artificial neural network (VI-BP-ANN) for quantitative TFe content in iron ores was first proposed. After the LIBS spectra of 80 representative iron samples were obtained, random forest (RF) was optimized by out-of-bag error and then used to measure and rank variable importance. The variable importance thresholds and the number of neurons were optimized by five-fold cross-validation (5-CV) with correlation coefficient (R2) and root mean square error (RMSE). With using only 1.40% of full spectral variables to construct BP-ANN model, the resulted R2, the root mean squared error of prediction (RMSEP) and the modeling time of the final VI-BP-ANN model was 0.9450, 0.3174 wt% and 24 s, respectively. Compared with full spectrum based-model (BP-ANN, RF, SVM, and PLS) and VI-RF model, the VI-BP-ANN model reduced over fitting and obtained the highest R2 and the lowest RMSE both for calibration and prediction. Meanwhile, characteristics of variables selected by variable importance were analyzed. In addition to the elemental emission lines of Ca, Al, Na, K, Mn, Si, Mg, Ti, Zr, and Li, partial spectral baselines of 540-610 nm and 820-970 nm are also selected as characteristic variables, which indicate that variable importance can take into full consideration the elemental interactions and the spectral baselines. Our approach shows that LIBS combined with VI-BP-ANN is able to rapidly and accurately quantify TFe content in iron ores.