Grasslands provide a wide range of ecosystem services (ESs). However, there is currently no method for easily diagnosing the level of ESs produced. Our aim was to develop ES indicators based on botanical surveys, which are readily available data and integrative of grassland spatiotemporal variability. Methods Based on academic knowledge and expertise, we identified several simple vegetation criteria that we aggregated using a multicriteria analysis tool to construct indicators of the level of ESs provided by grasslands. In this study, the indicators were calculated from over 2000 botanical surveys spread over a wide biogeographical gradient. Results Analyses of correlation between the various indicators show that “forage supply” and “diversity conservation” were not correlated. “Forage availability” and “nitrogen availability for the vegetation” were positively linked together and negatively linked to the robustness of the plant community to extreme events. A temporal approach highlights that the “biodiversity conservation” score decreased from 1970 to 2010 and that “nitrogen availability for the vegetation” was lower in 1970 and 1980 than in 2000 and 2010. Conclusions These results show that our aggregation method based on a large data set of botanical surveys could be appropriate for studying temporal dynamics of ESs.