Abstract. This paper presents ensemble simulations with the global climate model developed at the A. M. Obukhov Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS CM). These simulations were forced by historical reconstruction of external forcings for 850–2005 AD and by the Representative Concentration Pathways (RCP) scenarios till year 2300. Different ensemble members were constructed by varying the governing parameters of the IAP RAS CM module to simulate natural fires. These members are constrained by the GFED–3.1 observational data set and further subjected to Bayesian averaging. This approach allows to select only changes in fire characteristics which are robust within the constrained ensemble. In our simulations, the present-day (1998–2011 AD) global area burnt due to natural fires is (2.1 ± 0.4) × 106 km2 yr−1 (ensemble means and intra-ensemble standard deviations are presented), and the respective CO2 emissions in the atmosphere are (1.4 ± 0.2) PgC yr−1. The latter value is in agreement with the corresponding observational estimates. Regionally, the model underestimates CO2 emissions in the tropics; in the extra-tropics, it underestimates these emissions in north-east Eurasia and overestimates them in Europe. In the 21st century, the ensemble mean global burnt area is increased by 13% (28%, 36%, 51%) under scenario RCP 2.6 (RCP 4.5, RCP 6.0, RCP 8.5). The corresponding global emissions increase is 14% (29%, 37%, 42%). In the 22nd–23rd centuries, under the mitigation scenario RCP 2.6 the ensemble mean global burnt area and respective CO2 emissions slightly decrease, both by 5% relative to their values in year 2100. Under other RCP scenarios, these variables continue to increase. Under scenario RCP 8.5 (RCP 6.0, RCP 4.5) the ensemble mean burnt area in year 2300 is higher by 83% (44%, 15%) than its value in year 2100, and the ensemble mean CO2 emissions are correspondingly higher by 31% (19%, 9%). All changes of natural fire characteristics in the 21st–23rd centuries are associated mostly with the corresponding changes in boreal regions of Eurasia and North America. However, under the RCP 8.5 scenario, increase of the burnt area and CO2 emissions in boreal regions during the 22nd–23rd centuries are accompanied by the respective decreases in the tropics and subtropics.