A dynamic model has been developed to describe the performance of an anoxic biotrickling filter for biogas desulfurization. The model considers the most relevant phenomena involved in the biotrickling filter operation: convection, absorption, diffusion and biodegradation. The model also considers that a fraction of the liquid phase is stagnant -an assumption that increases the importance of diffusion phenomena for low liquid flow rates. The model was calibrated and validated using experimental data from a pilot-scale plant installed in a WWTP. In the calibration stage a set of periods with a wide range of operating conditions was used; i.e., biogas flow rate in the range 1-5 m 3 h −1 , recirculation flow rate in the range 1-3 m 3 h −1 , and nitrate concentration in the range 1-423 gN−NO3 -m −3 . The predictions obtained on using the model were consistent with experimental data and the divergence was less than 2%.The model was subsequently validated using two faultless periods (recirculation flow rates of 1.5 and 3 m 3 h −1, biogas flow rate in the range 1-5.2 m 3 h −1 , and inlet H2Sconcentration steps in the range 3600-5500 ppmV). An ANOVA study was carried out in order to quantify the suitability of the predictions. The results indicated that the differences between experimental and simulated outlet H2S concentrations were not statistically significant. The model was also able to predict simultaneously the dynamic concentrations of sulfide and nitrate in the liquid phase. Once the model had been validated, six control strategies were analyzed for different scenarios and purposes:i.e., to minimize the nitrate consumption and/or to maximize the H2S removal efficiency. The developed model is a potential tool to enhance and optimize the performance of biotrickling filters for the anoxic removal of H2S.3