Abstract. The Bulgarian Geomagnetic Reference Field (Bul-GRF) for 2015.0 epoch and its secular variation model prediction up to 2020.0 is produced and presented in this paper. The main field model is based on the well-known polynomial approximation in latitude and longitude of the geomagnetic field elements. The challenge in our modelling strategy was to update the absolute field geomagnetic data from 1980.0 up to 2015.0 using secular measurements unevenly distributed in time and space. As a result, our model gives a set of six coefficients for the horizontal H , vertical Z, total field F , and declination D elements of the geomagnetic field. The extrapolation of BulGRF to 2020 is based on an autoregressive forecasting of the Panagyurishte observatory annual means. Comparison of the field values predicted by the model with Panagyurishte (PAG) observatory annual mean data and two vector field measurements performed in 2015 shows a close match with IGRF-12 values and some difference with the real (measured) values, which is probably due to the influence of crustal sources. BulGRF proves to be a reliable alternative to the global geomagnetic field models which together with its simplicity makes it a useful tool for reducing magnetic surveys to a common epoch carried out over the Bulgarian territory up to 2020.