Control interventions in sustainable pest management schemes are set according to the phenology and the population abundance of the pests. This information can be obtained using suitable mathematical models that describe the population dynamics based on individual life history responses to environmental conditions and resource availability. These responses are described by development, fecundity and survival rate functions, which can be estimated from laboratory experiments. If experimental data are not available, data on field population dynamics can be used for their estimation. This is the case of the extrinsic mortality term that appears in the mortality rate function due to biotic factors. We propose a Bayesian approach to estimate the probability density functions of the parameters in the extrinsic mortality rate function, starting from data on population abundance. The method investigates the time variability in the mortality parameters by comparing simulated and observed trajectories. The grape berry moth, a pest of great importance in European vineyards, has been considered as a case study. Simulated data have been considered to evaluate the convergence of the algorithm, while field data have been used to obtain estimates of the mortality for the grape berry moth.