The aim of this work is to provide the strong convergence results of numerical approximations of a general second order non-autonomous semilinear stochastic partial differential equation (SPDE) driven simultaneously by an additive fractional Brownian motion (fBm) with Hurst parameter H ∈ ( 1 2 , 1) and a Poisson random measure, more realistic in modelling real world phenomena. Approximations in space are performed by the standard finite element method and in time by the stochastic Magnus-type integrator or the linear semi-implicit Euler method. We investigate the mean-square errors estimates of our fully discrete schemes and the results show how the convergence orders depend on the regularity of the initial data and the driven processes. To the best of our knowledge, these two schemes are the first numerical methods to approximate the non-autonomous semilinear stochastic partial differential equation (SPDE) driven simultaneously by an additive fractional Brownian motion with Hurst parameter H and a Poisson random measure.