A system (100), method and computer program product for optimization of crop protection, for a particular crop in a particular field. The system receives a predetermined future time interval and real-time weather forecast data and can retrieve corresponding predicted target pressure data (122) representing expected biotic stress factor probabilities for said crop from a crop state prediction model. An optimizer (150) receives the predicted target pressure data (122) for one or more targets and determines an initial target pressure score for the particular field by aggregating all target pressure data in the predetermined future time interval. A protection product recommender (270) provides a plurality of product recommendations with each product recommendation including a mix of protection products suitable for treatment of the one or more targets, and each protection product having protection product characteristics (265) describing the treatment effect of the respective protection product by target. The optimizer determines further target pressure scores for a plurality of combinations of product recommendations and product application time slots by applying the respective protection product characteristics (265) to the target pressure data, and selects the combination with the target pressure score indicating the lowest expected biotic stress for the crop in the field as recommended combination.