Understanding and predicting forest fires have proved a highly difficult endeavour, which requires extending and adapting complex models used in different fields. Here, we apply a marked point process approach, commonly used in ecology, which uses multiple Gaussian random fields to represent dynamics of Mediterranean forest fires in a spatio-temporal distribution model. Inference is carried out using Integrated Nested Laplace Approximation (INLA) with inlabru, an accessible and computationally efficient approach for Bayesian hierarchical modelling, which is not yet widely used in species distribution models. Using the marked point process approach, intensity of forest fires and dispersion were predicted using socioeconomic factors and environmental and fire-related variables. This demonstrates the advantage of complex model components in accounting for spatio-temporal dynamics that are not explained by environmental variables. Introduction of spatio-temporal marked point process can provide a more realistic perspective of a system, which is of particular importance for a practical and impact-focused worldwide problem such as forest fires.