Rémy Ballot
email: remy.ballot@inrae.fr

Nicolas Guilpart

Marie-Hélène Jeuffroy

The first map of crop sequence types in Europe over

Crop diversification is considered a key element of agroecological transition, whereas current dominant cropping systems are known to rely on only a few crop species -like cereals in Europe. To assess the benefits of crop diversification at a large scale, an accurate description of current crop sequences is required as a baseline. However, such a description is lacking at the scale of Europe. Here, we developed the first map of dominant crop sequences in Europe for the period 2012-2018. We used the Land Use Cover Area frame statistical Survey (LUCAS) dataset that provides temporally incomplete (2012, 2015 and 2018) land cover information from a stable grid of points covering Europe. Eight crop sequence types were identified using hierarchical clustering implemented on LUCAS data and mapped over Europe. We show, in France, that the relative importance of these eight crop sequence types (as estimated from LUCAS data) was highly consistent with those derived from an almost spatially exhaustive temporally complete national dataset (the French Land Parcel Identification System) for the same period, thus validating the method and typology for this country. Land use (i.e. crop production area) derived from our map of dominant crop sequences was also highly consistent with land use reported by official statistics at both national and European levels, validating the approach at a European scale. This first map of dominant crop sequences in Europe should be useful for future studies dealing with agricultural issues that are sensitive to crop rotations. The map of dominant crop sequence types in Europe derived from

Introduction

Crop diversification -the increase in crop diversity from the field to national scales -is considered a key element of agroecological transition (Kremen et al., 2012;[START_REF] Lechenet | Profiling farming management strategies with contrasting pesticide use in France[END_REF][START_REF] Tilman | Cultivate biodiversity to harvest food security and sustainability[END_REF]Wezel et al., 2014), whereas current crop sequences are short and specialised in many parts of the world (Meynard et al., 2018;Salembier et al., 2016;[START_REF] Schott | Les oléoprotéagineux dans les systèmes de culture: évolution des assolements et des successions culturales depuis les années 1970 dans le bassin de la Seine[END_REF][START_REF] Stein | Identifying crop rotation practice by the typification of crop sequence patterns for arable farming systems -A case study from Central Europe[END_REF]. Indeed, crop diversification supports a number of ecosystem services, with positive effects on crop yield, soil fertility, nutrient cycling, carbon sequestration, climate and water regulation, pest control, biodiversity and pollination (Beillouin et al., 2019(Beillouin et al., , 2021a;;[START_REF] Tamburini | Agricultural diversification promotes multiple ecosystem services without compromising yield[END_REF]. Strategies of crop diversification practices include crop rotation, cultivar mixture, cover crops, intercropping and agroforestry. Among these strategies, diversifying the crop sequence is of partic-ular interest in order to (i) increase nutrient availability and limit fertiliser requirements (e.g. by the inclusion of legume crops), (ii) favour soil protection and conservation by enhancing soil cover, (iii) promote natural regulation of pests and diseases by avoiding the presence of successive host crops and (iv) reduce weed infestation (Wezel et al., 2014). The term "crop sequence" refers to the sequence of crops grown in succession in the same field over a given period of time (Dury et al., 2012). A crop sequence is then defined by the nature of its crops and their order of succession. Based on this definition, the temporal frequencies of crops are key features of a crop sequence. The term "crop rotation" is also commonly used to refer to the sequence of crops grown in succession in the same field [START_REF] Bullock | Crop rotation[END_REF] but includes a notion of cyclicality (e.g. rotation length) -at least to some degree [START_REF] Castellazzi | A systematic representation of crop rotations[END_REF]. Hereafter, we use the term Published by Copernicus Publications. crop sequence rather than crop rotation because we consider a fixed period of time (i.e. 2012-2018), focus on temporal frequencies of crops and do not analyse cyclicality in crop sequences.

Aiming at leveraging benefits from diversifying crop sequences, several foresight studies have assessed scenarios for the future of agriculture in which crop sequences are modified to increase the services they provide [START_REF] Billen | Reshaping the European agro-food system and closing its nitrogen cycle: The potential of combining dietary change, agroecology, and circularity[END_REF][START_REF] Poux | An agroecological Europe in 2050: multifunctional agriculture for healthy eating, Findings from the Ten Years For Agroecology (TYFA) modelling exercise[END_REF] or scenarios of organic farming expansion where significant differences in crop rotations between conventional and organic agriculture are acknowledged (Barbieri et al., 2017(Barbieri et al., , 2019(Barbieri et al., , 2021)). A robust assessment of benefits expected from crop diversification at a large scale requires an accurate description of current crop sequences as a baseline. Indeed, several ecosystem services provided by diversified farming systems depend on the precrop effect and thus need to know precisely the nature of the crops successively grown on the same field [START_REF] Bennett | Meeting the demand for crop production: the challenge of yield decline in crops grown in short rotations[END_REF]. However, such a description is still lacking for Europe for at least two reasons: (i) existing datasets at the European scale, e.g. the Land Use Cover Area frame statistical Survey (LUCAS) dataset, provide information about land use categories and crop species cultivated on agricultural land but no information about crop sequences because data are not available every year (e.g. the LUCAS dataset provides information every 3 years); (ii) even if national and sub-national datasets describing crop sequences are available in some European Union (EU) state members, e.g. the Land Parcel Identification System in France [START_REF] Levavasseur | RPG Explorer: A new tool to ease the analysis of agricultural landscape dynamics with the Land Parcel Identification System[END_REF], a lack of harmonisation (e.g. spatial and temporal resolution or nomenclature) between these datasets makes them difficult to use in a cross-analysis. This situation has resulted in the production of regional to national studies of crop sequences [START_REF] Levavasseur | RPG Explorer: A new tool to ease the analysis of agricultural landscape dynamics with the Land Parcel Identification System[END_REF]Peltonen-Sainio and Jauhiainen, 2019;[START_REF] Stein | Identifying crop rotation practice by the typification of crop sequence patterns for arable farming systems -A case study from Central Europe[END_REF][START_REF] Xiao | Modeling the spatial distribution of crop sequences at a large regional scale using land-cover survey data: A case from France[END_REF] but has hampered any analysis at the scale of the EU. Such an analysis would be useful, especially in the context of the EU's Farm to Fork strategy, in which the adoption of more diverse crop rotations is encouraged (European Commission, 2022). To overcome this problem, we developed an original method that combines European-level and national-level datasets to create a map of current dominant crop sequences at the European level. We show that the temporally incomplete (i.e. every 3 years) information on crop sequences provided by the LUCAS dataset can be used to derive robust estimates of crop frequencies in the sequence when compared to crop frequencies derived from a temporally complete national-level dataset. We acknowledge that this approach does not allow us to capture the exact order of crops in the succession. Nonetheless, it allows us to capture crop frequencies, which are a key feature of crop sequences [START_REF] Castellazzi | A systematic representation of crop rotations[END_REF]Peltonen-Sainio and Jauhiainen, 2019).

Materials and methods

Method overview

The purpose of this work is to map current dominant crop sequences from the European LUCAS dataset. For this study, the multi-year harmonised data from d 'Andrimont et al. (2020) were used. As this dataset is temporally incomplete (observations of land use on fixed points only every 3 years), we proceeded in three steps (Fig. 1) to assess how this incomplete information could be used to describe the diversity and localisation of major crop sequences across Europe. First, we filtered the LUCAS data by selecting points under non-perennial agricultural land cover. For each of these points, we calculated frequencies of eight crops or groups of crops across 2012, 2015 and 2018 (the last 3 years of observation). For example, a point identified as wheat in 2012 and as corn in 2015 and 2018 was converted into cereal frequency = 0.33 and corn frequency = 0.67. Second, based on this set of eight variables, LUCAS points were classified into eight groups (hereafter referred to as crop sequence types) by combining principal component analysis (PCA) with hierarchical clustering. Third, we assessed the consistency of this classification against crop sequences derived with the quasiexhaustive yearly French Land Parcel Identification System (LPIS) and with the Eurostat crop production data. All data sources used in this study are detailed in Table 1.

Data sources

The LUCAS dataset

The LUCAS dataset is an in situ land cover (i.e. physical cover observed at the earth's surface, e.g. cereals, root crops, fodder crops) and land use (i.e. socio-economic function of the observed earth's surface, e.g. agriculture, industry, residential) observation carried out on approximately 300 000 points sampled from a stable grid of around 1 100 000 points across Europe every 3 years since 2006. The last observation campaign was carried out in 2018. The work presented here is based on the harmonised LUCAS database developed by d 'Andrimont et al. (2020), which gathered, in a unique dataset, all the information collected from the beginning of this survey (Table 1). All 28 EU member states are represented in this dataset, except Croatia and Malta. Montenegro, which is not a member state, is also represented, as well as the United Kingdom, which is no longer a member state. Thus, "Europe" will be used hereafter to refer to the spatial scope of this study. 2009, 2012, 2015 and 2018. This choice would also have limited the analysis to 11 of the 27 EU countries (i.e. Belgium, Czechia, Germany, Spain, France, Hungary, Italy, Luxembourg, the Netherlands, Poland and Slovakia). We discarded points labelled as non-agricultural use or permanent agricultural use (e.g. orchards, vineyards) in at least 1 year among 2012, 2015 and 2018 (i.e. land cover (LC) not included in land cover classifications B11 to B55). As an exception, we also conserved 2609 points identified as permanent grasslands in 2018 (i.e. LC E20) but identified as non-permanent agricultural uses in 2012 or 2015, as they may be the result of a confusion between temporary and permanent grasslands during observation. This resulted in a dataset of 21 620 points with information about crops cultivated in 2012, 2015 and 2018. Thus, each point is associated with a single (temporally incomplete) crop sequence.

The French LPIS dataset

Depending on the purpose of an analysis, crop sequences can be described in many different ways, including length and flexibility of rotation [START_REF] Castellazzi | A systematic representation of crop rotations[END_REF], nature and function of crops or crop groups (Barbieri et al., 2017;[START_REF] Videla-Mensegue | Functional crop types are more important than diversity for the productivity, profit and risk of crop sequences in the inner Argentinean Pampas[END_REF], diversity of crops or crop groups (Beillouin et al., 2021b;[START_REF] Tamburini | Agricultural diversification promotes multiple ecosystem services without compromising yield[END_REF], diversity of sowing dates [START_REF] Weisberger | Does diversifying crop rotations suppress weeds? A meta-analysis[END_REF], order of succession of crops within the sequence (Peltonen-Sainio et al., 2019) and temporal frequency or return time of crops within the sequence [START_REF] Nowak | Assessment of the diversity of crop rotations based on network analysis indicators[END_REF]. Given the characteristics of the LUCAS dataset, the exact order of crop succession within the sequences could not be described, nor could the variety of sowing dates (e.g. the LUCAS dataset provides no information about whether the crop is sown in autumn or spring). Therefore, we chose to describe current dominant crop sequences by the temporal frequencies, over 3 years (i.e. 2012, 2015 and 2018), of eight crops or groups of crops described below. This choice is a compromise between addressing common issues related to crop rotations (e.g. nitrogen management and the role of legumes) and avoiding unnecessary complexity. These eight groups were defined based on two criteria: (i) their agronomic relevance and (ii) their importance in terms of cultivated area in Europe. The groups considered were (i) cereals (i.e. wheat, barley, oat, triticale and rye, corresponding to LC B11 to B19 except B16 in the LUCAS dataset nomenclature), (ii) corn (LC B16), (iii) rapeseed (LC B32), (iv) sunflower (LC B31), (v) pulses (i.e. dry pulses and soybean, corresponding to LC B33 and B41), (vi) root crops (i.e. beets and potatoes, corresponding to LC B21 to B23), (vii) forage legumes (i.e. alfalfa and clover, corresponding to LC B51 and B52) and (viii) temporary grassland (LC B53, B55 and E20, Table 2). As these eight groups did not encompass all the land cover categories, the sum of (groups of) crop frequencies may be lower than 1 for each point.

In order to serve the quality assessment step, the French LPIS dataset was preprocessed in the same way. First, fields under perennial use were discarded. Then, (groups of) crop frequencies were calculated for each field.

Identification of crop sequence types with hierarchical clustering

PCA was performed with the PCA() function of the R package FactoMineR v2.1 [START_REF] Lê | FactoMineR: an R package for multivariate analysis[END_REF] on all eight variables describing temporal frequencies of cereals, rapeseed, sunflower, pulses, corn, root crops, forage legumes and temporary grassland. Seven components displayed eigenvalues higher than or close to 1, totalling more than 99 % of variance (Fig. 2a). These seven components were used as inputs for a hierarchical clustering performed with the HCPC() function from the same R package. Inertia gain displayed a clear break from eight to nine (Fig. 2b), and thus eight groups of crop sequences (crop sequence types) were considered. A crop sequence type has thus been assigned to all points considered. To analyse the geographical distribution of dominant crop sequences, observed locations of these eight crop sequence types were mapped across Europe.

Quality assessment

The crop sequence types derived from incomplete temporal sequences (LUCAS dataset) were compared to crop sequence types derived from complete temporal sequences based on the French LPIS dataset. The French LPIS dataset represents the best available data for spatial distributions of crops in France in terms of spatial and temporal resolutions, spatial coverage and disaggregation by crop type, with a cov- erage higher than 98 % for all field crops or 92 % for temporary grasslands [START_REF] Cantelaube | Construction d'une Base de Données Géographiques à Échelle Fine Exhaustive Sur l'occupation Agricole Du Sol: Le " RPG Complété[END_REF][START_REF] Guilpart | How much agricultural land is there close to residential areas? An assessment at the national scale in France[END_REF]. With this aim, a three-step procedure was followed. First, a random forest (RF) model was trained on LUCAS data to predict crop sequence types from crop (or crop group) frequencies (i.e. a total of eight predictors). The RF model was fitted using the randomForest() function of the R package randomForest v4.6.14 [START_REF] Liaw | Classification and regression by random-Forest[END_REF], with default settings (i.e. 500 trees, two variables randomly sampled as candidates at each split). The RF model showed good performances, as indicated by an out-of-bag (OOB) error rate lower than 0.1 % (Table 4). Each tree of our RF model is constructed based on a random sample of the observations generated by bootstrap. The observations that are not part of the bootstrap sample are referred to as OOB observations, which are used to estimate the prediction error (the so-called OOB error). The OOB error is considered a good measure of the true prediction error [START_REF] Matthew | Bias of the random forest out-of-bag (OOB) error for certain input parameters[END_REF][START_REF] Janitza | On the overestimation of random forest's out-of-bag error[END_REF]. Second, the RF model was applied to French LPIS data to classify observed crop sequences into the eight crop sequence types. Third, the crop sequence type distributions derived from the LUCAS dataset and predicted from the French LPIS dataset were compared at the national and regional levels in France. This allowed us to check whether the lower temporal and spatial resolutions of the LUCAS dataset led to some bias in the identification of crop sequence types and their geographical distribution.

To ensure that the map of crop sequence types presented in this paper could be used as a baseline in studies assessing scenarios for the future of agriculture in Europe, it appears essential to assess its consistency with observed land use. With this aim, average crop or crop group harvested areas between 2012 and 2018 were calculated from (i) their temporal frequencies within the crop sequence types derived from the LUCAS dataset, (ii) the relative importance (i.e. spatial frequency) of crop sequence types and (iii) the total arable land area following Eq. ( 1):

harvested_area i = j F j,i , F j total_area.
(1)

F j is the spatial frequency of crop sequence type j , and F i,j is the temporal frequency of (group of) crop i within crop sequence type j . This calculation was applied for each country in Europe, at the national level, and compared to crop or crop group harvested areas reported by Eurostat (2023).

Relative importances of crop sequence types were calculated according to the number of points and without consideration of the field area. Indeed, the harmonised LUCAS dataset provides information about the sizes of the surveyed parcels in hectares, but this information is limited to four categories (i.e. < 0.5, 0.5-1, 1-10 and > 10 ha) that were not relevant to weight relative importance.

Results

Eight types of crop sequences identified based on the LUCAS dataset

Eight different types of crop sequences were identified based on the LUCAS dataset ( and cereals (17 %). The remaining three crop sequence types can be considered minor as they individually account for less than 10 % of LUCAS points. This includes a crop sequence type dominated by sunflower (42 %) and cereals (39 %), a second one dominated by pulses (37 %) and cereals (36 %) and a third one dominated by forage legumes (47 %) and cereals (23 %). Two (groups of) crops are present in all crop sequence types: cereals and corn (Table 2). This highlights the central role of these crops in crop rotations in Europe. On the other hand, pulses, and to a lesser extent sunflower, root crops and forage legumes, appear in only a few (i.e. one to three) crop sequence types. Some crop sequence types show a high degree of diversity, like the "pulses and cereals" in which all (groups of) crops are represented, while some others are much less diverse, like the "cereals" and "corn and cereals" types. These two crop sequence types account for almost half (45 %) of all LUCAS points and are composed of cereals and corn by more than 90 %.

If all the LUCAS campaigns available were considered (i.e. 2006-2018), the classification would have resulted in the same eight types of crop sequences. The eight types would not have ranked in the same order regarding their relative importance. Their relative importance would have changed slightly: 26 instead of 35 % for a cereal-dominated crop sequence type (same ranking), 17 instead of 8 % for sunflower and cereals (second instead of sixth place), 16 instead of 10 % for corn and cereals (third instead of fifth place), 13 instead of 11 % for a temporary grassland-dominated crop sequence type (fourth instead of third place), 8 instead of 13 % for rapeseed and cereals (fifth instead of second place), 8 instead of 10 % for root crops and cereals (sixth instead of fourth place), 7 % for pulses and cereals (relative importance and ranking unchanged) and 6 instead of 7 % for forage legumes and cereals (ranking unchanged). These differences were consistent as 16 of the 27 countries were discarded when we only considered points with observations for all the 2006-2018 LUCAS campaigns.

3.2 Quality assessment 3.2.1 Incomplete temporal crop sequences are good proxies of complete temporal crop sequences at the national and regional scales in France

On average at the national scale in France, the relative importance of crop sequence types estimated from the LUCAS dataset is in good agreement (R 2 = 0.75, root mean square error RMSE = 0.04) with estimates based on the French LPIS dataset, with no systematic bias (Figs. 3 and S1 in the Supplement). However, estimates based on the LUCAS dataset underestimate the proportion of the "grasslands" crop sequence type (-9 %) and slightly overestimate the proportion of the "corn and cereals" (+3 %), "cereals" (+2.5 %), "pulses and cereals" (+2 %) and "sunflower and cereals" (-2 %) ones. At the regional scale, the agreement between the relative importance of crop sequence types derived from the LUCAS or French LPIS datasets remains good but varies by crop group (Fig. 4). The crop sequence types "sunflower and cereals", "root crops and cereals", "rapeseed and cereals", "forage legumes and cereals" and "corn and cereals" show the best agreement, with high R 2 (higher than 0.8), low RMSE (lower than 0.1) and no systematic bias, whereas the crop sequence type "pulses and cereals" shows the lowest agreement (R 2 = 0.37, RMSE = 0.04). Underestimation of the crop sequence type "grasslands", especially for regions where relative importance is high, and overestimation of the crop sequence type "cereals" are confirmed but for regions where relative importance is low. Overall, regional differences in crop sequence types are well captured with the LUCAS dataset, as shown by two conclusions. First, any given crop sequence type may have a low relative importance in one region and a high one in another. This range of relative importance across regions is well captured by estimates derived from the LUCAS dataset for all crop sequence types (Figs. 4,S2). Second, regional specificities in terms of dominant crop sequence types are also well captured by LUCAS data (Fig. S2). For example, the Picardie region is dominated by the three crop sequence types "root crops and cereals", "cereals" and "rapeseed and cereals", while the Aquitaine region is dominated by the crop sequence types "corn and cereals" and "sunflower and cereals".

Land use derived from crop sequence types is in good agreement with land use reported by official statistics at national and EU levels

Comparison of crop harvested areas derived either from crop sequence types or from official statistics (EUROSTAT, 2022) shows good agreement at the European scale (Fig. 5), with R 2 higher than 0.99 and with no bias. Comparison for each country of the EU (Fig. 6) reveals good levels of correlation between areas derived from crop sequence types and areas reported by Eurostat, with R 2 ranging from 0.71 for temporary grasslands to 0.98 for cereals, rapeseed or pulses. Nevertheless, some crop-country combinations display bad agreement between estimated and observed harvested areas (Fig. S3).

Spatial distribution of crop sequence types in Europe

The maps of crop sequence types presented in Figs. 7 and8 show a strong spatial pattern of crop sequence type distribution in the EU. The crop sequence types "grasslands", "corn and cereals", "cereals" and, to a lesser extent, "pulses and cereals" and "forage legumes and cereals" are present in all EU countries, whereas the other crop sequence types are concentrated in specific parts of Europe. Indeed, "rapeseed and cereals" and "root crops and cereals" are mainly found in the north of France, Belgium, the Netherlands, Ger-Earth Syst. Sci. Data, 15, 5651-5666, 2023 https://doi.org/10.5194/essd-15-5651-2023 many, Czechia and Poland, whereas "sunflower and cereals" is mainly found in the south of France, Spain, Hungary, Bulgaria and Romania. These differences are best highlighted when looking at the latitudinal distribution of crop sequence types shown in Fig. 9. Indeed, "grasslands", "pulses and cereals" and "cereals" are found in a wide range of latitudes (from 40 to 60 • N), while other crop sequence types are limited to narrower ranges of latitudes, with some in the north, like "root crops and cereals", concentrated between 50 and 55 • N and some in the south, like "sunflower and cereals", concentrated between 40 and 50 • N. Most crop sequence types display regions in which they are more concentrated, like the north of France and Belgium for "root crops and cereals" or central Italy for "forage legumes and cereals". However, the "pulses and cereals" crop sequence type seems more regularly distributed from 40 to 60 • N. Several crop sequence types coexist in most regions. For example, all crop sequence types except "sunflower and cereals" and "forage legumes and cereals" are found in the very north of France (Fig. 8). It is well known that temporal and spatial crop diversities are not independent of each other (Aramburu Merlos and Hijmans, 2020). Therefore, it is expected that a crop sequence type characterised by a high temporal frequency of a given (group of) crop(s) will be frequent where this (group of) crop(s) is widely cultivated. For example, crop sequences including corn can only be observed where corn is grown. However, knowing where corn is grown does not say any-thing about the crop sequence in which corn is cultivated. Of course, knowing which other crops are grown in the same area as corn can inform us about possible crop sequences, but this is not sufficient. Our results provide a good example: many crop sequence types coexist in the very north of France (temporary grasslands, corn and cereals, root crops and cereals, rapeseed and cereals, pulses and cereals, and cereals). As a consequence, cereals (e.g. wheat) can be found in very specialised crop sequences (e.g. the "cereals" crop sequence type), moderately diversified crop sequences (e.g. the "root crops and cereals" crop sequence type) and diversified crop sequences (e.g. the "pulses and cereals" crop sequence type). This demonstrates that specialised crop sequences can still occur in areas where a substantial diversity of crops is cultivated, and this cannot be inferred from land use (e.g. individual crop maps of "rotation heads") only.

Data availability 6 Discussion

On the quality of the map of crop sequence types

The maps of crop sequence types derived from our study fill an important void in our knowledge of the spatial distribution of crop sequences in Europe. Despite them being based on incomplete data in both time and space (the LUCAS dataset), these maps have been shown to represent well the distribution of crop sequence types in France at both national and regional scales when compared to an almost spatially exhaustive dataset of temporally complete crop sequences. This consistency shows that the temporally incomplete information from LUCAS (i.e. only three crops known -2012, 2015 and 2018 -from a 7-year crop sequence) is a good proxy to temporally complete crop sequences at regional to national scales, although individual sequence classification at the field scale may be prone to some bias. At least three different types of bias at the individual field level can be identified. First, let us consider a 3-year cereal-beet-potato rotation, which is quite common in north-western Europe. This rotation may appear in the LUCAS dataset as cereal-cerealcereal, beet-beet-beet or potato-potato-potato, depending on the crop at the time of observation, and may therefore be classified as the "cereals" or "root crops and cereals" crop https://doi.org/10.5194/essd-15-5651-2023

Earth Syst. Sci. Data, 15, 5651-5666, 2023 sequence type. However, we believe that this individual field classification bias does not affect our results at larger scales, as all survey points are not likely to be synchronised in terms of crop sequence.

Second, the comparison with the French LPIS dataset revealed an underestimation of grassland crop sequences, which led to an overestimation of other crop sequence types (Fig. 3). A possible explanation for this underestimation could be the inability of the 3 sample years available in the Earth Syst. Sci. Data, 15, 5651-5666, 2023 https://doi.org/10.5194/essd-15-5651-2023 Note that crop sequence types have been given short names that reflect dominant crops or crop groups within the sequence. For example, the crop sequence type "grasslands" also includes cereals and corn, but these crops are less frequent in the sequence. The composition of each crop sequence type in terms of temporal frequencies of (groups of) crops is presented in Table 2.

LUCAS dataset (i.e. 2012, 2015 and 2018) to capture the full diversity of longer crop rotations, such as rotations including temporary grassland. For example, let us consider a cyclical crop rotation starting in 2012 with 3 years of consecutive grassland followed by wheat, wheat, maize and barley. Then, observations in 2012, 2015 and 2018 would be grassland, wheat and barley respectively, yielding a proportion of grassland of one-third instead of one-half. Another possible explanation could be that temporary grassland may sometimes be confused with permanent grassland during LUCAS observations as well as during farmer declarations in the French LPIS. In the LUCAS dataset, we decided to consider as temporary grasslands points observed as permanent grassland in 2018 but as non-perennial agricultural cover in 2012 and 2015. These points are more likely temporary grasslands confused with permanent ones. However, they could also be actual permanent grasslands after a perennial change in land use between 2015 and 2012. This may contribute to the underestimation of crop sequences with temporary grassland importance. As highlighted by previous research (e.g. Martin et al., 2020), crop sequences with temporary grasslands provide specific services such as carbon sequestration, bio-diversity conservation or weed control. Future users of our dataset should therefore be aware that this underestimation of the importance of crop sequences with temporary grasslands may result in underestimation of such services as well as overestimation of the production from other overestimated crops.

Third, the reliability of our validation dataset also needs to be discussed. The French LPIS dataset is based on farmers' declarations, which may not be 100 % correct. The RPG Explorer algorithm, which compiles LPIS annual raw data into multi-annual crop sequence data, has been validated. However, it cannot be 100 % correct when identifying the most likely crop sequence and when more than one crop is declared per block. A robust estimation of farmer declaration error and how it could propagate to the crop sequence is challenging and could not be done in this study. However, to date, the dataset we used has been the most complete one regarding crop sequences in France.

The comparison between land use derived from our map of crop sequence types and land use derived from official statistics revealed good consistency at the European level, except for some crop-country combinations. This could be partly short names that reflect dominant crops or groups of crops within the sequence. For example, the crop sequence type "grasslands" also includes cereals and corn, but these crops are less frequent in the sequence. The composition of each crop sequence type in terms of temporal frequencies of (groups of) crops is presented in Table 2.

explained by the spatial sampling effort of LUCAS data, which can be quantified as the number of LUCAS points per unit area of cropland by country. Analysis of this variable revealed important variation across countries, with the spatial sampling effort ranging from 0.10 (Hungary) to 1.42 (Montenegro) points per thousand hectares of cropland, with a median of 0.22 (Table 5). The accuracy of the map of crop sequence types should therefore be analysed in the light of this varying spatial sampling effort. For example, our results pointed to harvested areas overestimated or underestimated for most (groups of) crops for Romania, which displays a sampling effort lower than the median (0.12). On the other hand, Montenegro displays the highest sampling effort (1.42) but also an important discrepancy between estimated and observed harvested areas. Indeed, this high sampling rate hides a low number of observations (12) for small arable land. Thus, the map accuracy may be limited either by a low sampling effort or by a low number of observations in small-sized countries.

Lastly, we note that the quality assessment of our map of crop sequence types in the EU would benefit from comparisons with other datasets at national or sub-national levels in other countries. Nevertheless, we highlight that all crop sequence types are present in France, allowing for a complete quality assessment over the eight crop sequence types.

Perspectives

We assume our maps of crop sequences to be useful for future studies dealing with agricultural issues that are sensitive to crop rotations, including nature-based pest control [START_REF] Lechenet | Profiling farming management strategies with contrasting pesticide use in France[END_REF], pesticide use intensity [START_REF] Jacquet | An economic analysis of the possibility of reducing pesticides in French field crops[END_REF], nitrogen management and cycling [START_REF] De Notaris | Nitrogen leaching: A crop rotation perspective on the effect of N surplus, field management and use of catch crops[END_REF], biodiversity [START_REF] Li | Crop rotation history constrains soil biodiversity and multifunctionality relationships[END_REF], soil carbon sequestration (King and Blesh, 2018), water resource management [START_REF] Yang | Effect of diversified crop rotations on groundwater levels and crop water productivity in the North China Plain[END_REF] and crop yield [START_REF] Bennett | Meeting the demand for crop production: the challenge of yield decline in crops grown in short rotations[END_REF]. We also highlight that the methodology presented in this paper that leverages temporally incomplete information about crop sequences based on the LUCAS dataset can be adapted to spe-cific objectives of future studies. For example, looking at the specific place of soybean within crop sequences in the EU would probably require us to consider this crop specifically instead of as part of the group "pulses". Different indicators could also be calculated to describe different facets of crop sequences relevant to different ecosystem services. For example, the diversity of crop sowing dates is relevant to weed control [START_REF] Weisberger | Does diversifying crop rotations suppress weeds? A meta-analysis[END_REF], the crop functional diversity is relevant to nitrogen cycling when considering legume and non-legume crops [START_REF] Zhao | Nitrogen fixation and transfer between legumes and cereals under various cropping regimes[END_REF], the phylogenetic crop species diversity might be relevant or not to the control of pests and diseases depending on whether considered crop species are hosts of the same pests and diseases or not (e.g. wheat and barley vs. wheat and soybean) (Beillouin et al., 2021a), and the share of summer and winter crops is relevant to water management as water demand is usually higher during the summer. Here we chose the eight crops or groups of crops for their relevance to major agronomic issues (see Methods, Sect. 2.3.1). Cereals (wheat, barley, oat, triticale, rye) are all species of the Poaceae family, sharing similar features regarding their effects on crop sequences. For example, they all have straws with a relatively high C : N ra- short names that reflect dominant crops or groups of crops within the sequence. For example, the crop sequence type "grasslands" also includes cereals and corn, but these crops are less frequent in the sequence. The composition of each crop sequence type in terms of temporal frequencies of (groups of) crops is presented in Table 2. tio, which contributes to limiting nitrogen availability for the subsequent crop when not exported from the field (Justes et al., 2009). Pulses (e.g. pea, faba bean, soybean) are all from the Fabaceae family and have the ability to fix nitrogen from the atmosphere, so that they are usually not fertilised with nitrogen and contribute to increasing nitrogen availability to subsequent crops thanks to the relatively low C : N ratio of their crop residues (Guinet et al., 2020). Root crops (beets and potatoes) may have a specific effect on the topsoil and subsoil structure as belowground organs are harvested [START_REF] Gabarron-Galeote | BEETSOIL: a decision support tool for forecasting the impact of soil conditions on sugar beet harvest[END_REF]. Corn, sunflower and rapeseed are three widely grown crops and have specific characteristics, including low pesticide use for sunflower, high insecticide use for rapeseed [START_REF] Urruty | Impacts of agricultural land use changes on pesticide use in French agriculture[END_REF] and high irrigation requirements for maize, especially in the south of Europe [START_REF] Senthilkumar | Can farmers use maize earliness choice and sowing dates to cope with future water scarcity? A modelling approach applied to south-western France[END_REF]. Forage legumes and temporary grasslands are commonly cultivated for at least 2 or 3 years and present numerous benefits, including controlling weeds (Martin et al., 2020).

We hope the methodology developed here will stimulate a more detailed description of crop sequences, in both time and space, and their effect on agricultural production and sustainability in future studies. Review statement. This paper was edited by Chaoqun Lu and reviewed by three anonymous referees.

Figure 1 .

 1 Figure 1. Workflow diagram for developing spatially explicit maps of dominant crop sequences across Europe.

Figure 2 .

 2 Figure 2. Identification of crop sequence types based on the LU-CAS dataset. (a) Eigenvalues of the components from the principal component analysis. (b) Inertia gain according to the level of cutting. The hierarchical clustering is performed on 43 291 points from the harmonised LUCAS dataset for the years 2012, 2015 and 2018, based on the frequency of eight crop groups (see the text for details).

  rapeseed (43 %) and cereals (43 %); a second one dominated by grasslands (74 %), cereals (13 %) and corn (11 %); a third one dominated by root crops (43 %) and cereals (35 %); and a fourth one dominated by corn (82 %)

Figure 3 .

 3 Figure 3. Comparison of the relative importance of crop sequence types derived from the LUCAS and French LPIS datasets at the national scale in France. The dotted line denotes the 1 : 1 line. Also shown are the fitted linear regression model and the associated slope and intercept. RMSE: root mean square error.

Figure 4 .

 4 Figure 4. Comparison of the relative importance of crop sequence types derived from the LUCAS and French LPIS datasets for the 22 regions in France. The dotted line denotes the 1 : 1 line. Also shown are the fitted linear regression models and the associated slopes and intercepts. Labels identify the five regions with the highest difference between predicted and observed areas (the font size is proportional to the area difference relative to the observed area). AL: Alsace, AQ: Aquitaine, AU: Auvergne, BN: Basse-Normandie, BO: Bourgogne, BR: Bretagne, CE: Centre, CA: Champagne-Ardenne, CO: Corse, FC: Franche-Comté, HN: Haute-Normandie, IF: Île-de-France, LR: Languedoc-Roussillon, LI: Limousin, LO: Lorraine, MP: Midi-Pyrénées, NP: Nord-Pas-de-Calais, PL: Pays de la Loire, PI: Picardie, PC: Poitou-Charentes, PA: Provence-Alpes-Côte d'Azur, RA: Rhône-Alpes.

Figure 5 .

 5 Figure 5. Comparison of harvested areas reported by Eurostat and derived from crop sequence types based on the LUCAS dataset at the European scale. The dotted line denotes the 1 : 1 line. Also shown are the fitted linear regression model and the associated slope and intercept.

Figure 6 .

 6 Figure 6. Comparison of harvested areas reported by Eurostat and derived from the crop sequence types based on the LUCAS dataset for the various countries in Europe. The dotted line denotes the 1 : 1 line. Also shown are the fitted linear regression models and the associated slopes and intercepts. Labels identify the 10 countries with the highest difference between predicted and observed areas (the font size is proportional to the area difference relative to the observed area). AT: Austria, BE: Belgium, BG: Bulgaria, CY: Cyprus, CZ: Czechia, DE: Germany, DK: Denmark, EE: Estonia, EL: Greece, ES: Spain, FI: Finland, FR: France, HU: Hungary, IE: Ireland, IT: Italy, LT: Lithuania, LU: Luxembourg, LV: Latvia, MT: Montenegro, NL: Netherlands, PL: Poland, PT: Portugal, RO: Romania, SE: Sweden, SI: Slovenia, SK: Slovakia, UK: United Kingdom.

Figure 7 .

 7 Figure7. Map of crop sequence types in Europe based on the LUCAS dataset. Note that crop sequence types have been given short names that reflect dominant crops or crop groups within the sequence. For example, the crop sequence type "grasslands" also includes cereals and corn, but these crops are less frequent in the sequence. The composition of each crop sequence type in terms of temporal frequencies of (groups of) crops is presented in Table2.

Figure 8 .

 8 Figure 8. Individual maps of crop sequence types in Europe based on the LUCAS dataset. Note that crop sequence types have been given short names that reflect dominant crops or groups of crops within the sequence. For example, the crop sequence type "grasslands" also includes cereals and corn, but these crops are less frequent in the sequence. The composition of each crop sequence type in terms of temporal frequencies of (groups of) crops is presented in Table2.

Figure 9 .

 9 Figure 9. Latitudinal distribution of crop sequence types based on the LUCAS dataset. Note that crop sequence types have been given short names that reflect dominant crops or groups of crops within the sequence. For example, the crop sequence type "grasslands" also includes cereals and corn, but these crops are less frequent in the sequence. The composition of each crop sequence type in terms of temporal frequencies of (groups of) crops is presented in Table2.

Supplement.

  The supplement related to this article is available online at: https://doi.org/10.5194/essd-15-5651-2023-supplement.Author contributions. RB and MHJ initiated the research. RB, MHJ and NG designed the research. RB performed the analysis with the help of NG. RB and NG wrote the manuscript with contributions from MHJ.Competing interests. The contact author has declared that none of the authors has any competing interests. Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Acknowledgements. This work was supported by the LegValue project funded by the European Union's Horizon 2020 research and innovation programme under grant agreement no. 727672. Access to some confidential data, on which this work is based, has been made possible within a secure environment offered by CASD -Centre d'accès sécurisé aux données (ref. 10.34724/CASD).

Financial support .

 support This research has been supported by the European Union's Horizon 2020 framework programme (LegValue grant no. 727672).

  The French LPIS dataset is based on yearly declarations made by farmers in compliance with Common Agricultural Policy subsidies. Starting from 2008, it has provided annual almost spatially exhaustive information on land use for agri-

	Earth

Syst. Sci. Data, 15, 5651-5666, 2023 https://doi.org/10.5194/essd-15-5651-2023

Table 1 .

 1 Summary of the data sources used in this study.

	Name	Dataset description		References
	Harmonised LUCAS	Harmonised LUCAS in situ land cover	d'Andrimont et al. (2020)
	dataset	and use database for field surveys from	Available for download at https://data.jrc.ec.europa.eu/
		2006 to 2018 in the European Union		dataset/f85907ae-d123-471f-a44a-8cca993485a2
				(last access: 13 August 2020)
	French LPIS dataset	RPG Explorer Crop successions France	Martin et al. (2021)
		2007-2014, 2007-2019, 2015-2019		Available for download at https://data.inrae.fr/dataset.
				xhtml?persistentId=doi:10.15454/XH84QB
				(last access: 28 April 2021)
	Eurostat dataset	Annual crop acreage at the national	Available for download at
		scale for EU countries from 2000		https://ec.europa.eu/eurostat/
				estat-navtree-portlet-prod/BulkDownloadListing?
				file=data/apro_cpsh1.tsv.gz (last access: 27 April 2021)
	cultural land, detailing 28 crop categories until 2014 and ap-	al. (2021, Table 1). This dataset compiles all annual LPIS
	proximately 300 from 2015. Until 2014, information was col-	data for France into a unique dataset of crop sequences for
	lected at the block scale. Each block can enclose one or more	more than 18 million fields over the period 2012-2018 and
	agricultural parcels and thus one or more crops with a de-	associated areas. It also relies on an algorithm which identi-
	clared area for each one but no geolocalisation within the	fies the most likely crop sequence when more than one crop
	block. Due to parcel reconfiguration from one year to an-	is declared for one given block. This dataset represents the
	other, it is not straightforward to know pluriannual crop se-	most exhaustive dataset of crop sequences in France in terms
	quences from the French LPIS dataset for years older than	of both spatial and temporal resolution.
	2015. For our study, we used the RPG Explorer Crop suc-
	cessions France version 2.0 dataset developed by Martin et

Table 2 .

 2 Grouping of land cover categories considered for crop sequence classification from the harmonised LUCAS dataset and correspondence to the French LPIS and EUROSTAT datasets used for validation.

		Examples of crops	Harmonised LUCAS dataset	French LPIS dataset EUROSTAT dataset
	Cereals	Wheat, barley, oat, triticale, rye LC B11 to B15 and B17 to B19 ble, orge, cer_2nd	C1110 to C1420 and C1600 to C1900
	Corn		LC B16	mais	C1500, G3000
	Rapeseed		LC B32	colza	I1110
	Sunflower		LC B31	TRN	I1120
	Pulses	Dry pulses, soybean	LC B33 and B41	olea, prot	I1130, P110 to P9000
	Root crops	Beets, potatoes	LC B21 to B23	leg_fl, a_indus	R1000 to R3000
	Forage legumes	Alfalfa, clover	LC B51 and B52	fou	G2000
	Temporary grasslands		LC B53, B55 and E20	pt, pp	G1000

  Table 3). One crop sequence type largely dominates in terms of relative importance within Eu-

rope: it accounts for 35 % of LUCAS points and is dominated by cereals (79 %), followed by corn (8 %) and grasslands (6 %). Then comes a group of four crop sequence types of moderate relative importance, each one representing 10 % to 13 % of LUCAS points. This group includes a crop sequence https://doi.org/10.5194/essd-15-5651-2023

Earth Syst. Sci. Data, 15, 5651-5666, 2023

Table 3 .

 3 Average temporal frequency of crops or crop groups within each crop sequence type. Values in bold indicate the dominant (or two dominant) crop(s) or group(s) of crops in the crop sequence type. Values in brackets indicate standard deviations.

			Temporal frequency of crops or crop groups over the 2012-2018 period		Relative importance
	Crop sequence type (short name)	Cereals	Corn Rapeseed Sunflower Pulses	Root Grassland crops	Forage Total legumes	of crop sequence types across the EU
	1 -Temporary grassland	0.13	0.11	0	0	0	0	0.74	0	0.98 11 %
		(0.16) (0.22)	(0)	(0)	(0)	(0)	(0.20)	(0)
	2 -Forage legumes and cereals	0.23	0.11	0.02	0.02	0	0	0.13	0.47	0.98 7 %
		(0.25) (0.19)	(0.09)	(0.08)	(0) (0.02)	(0.22)	(0.20)
	3 -Corn and cereals	0.17	0.82	0	0	0	0	0	0	0.99 10 %
		(0.17) (0.17)	(0)	(0)	(0)	(0)	(0.10)	(0)
	4 -Root crops and cereals	0.35	0.09	0.05	0	0	0.43	0.04	0.01	0.97 10 %
		(0.26) (0.18)	(0.13)	(0)	(0) (0.19)	(0.13)	(0.05)
	5 -Rapeseed and cereals	0.43	0.09	0.43	0	0	0	0.03	0	0.98 13 %
		(0.25) (0.18)	(0.19)	(0)	(0)	(0)	(0.12)	(0)
	6 -Sunflower and cereals	0.39	0.11	0.03	0.42	0	0.01	0.02	0	0.98 8 %
		(0.25) (0.20)	(0.11)	(0.16)	(0) (0.05)	(0.10)	(0)
	7 -Pulses and cereals	0.36	0.09	0.04	0.03	0.37	0.03	0.04	0.01	0.97 7 %
		(0.26) (0.19)	(0.13)	(0.10) (0.12) (0.10)	(0.12)	(0.07)
	8 -Cereals	0.79	0.08	0	0	0	0	0.06	0	0.93 35 %
		(0.16) (0.13)	(0)	(0)	(0)	(0)	(0.11)	(0)

Table 4 .

 4 Confusion matrix of the random forest model.

	Observed class	Predicted class	Error rate (%)

Table 5 .

 5 Spatial sampling effort of LUCAS data by country. The sampling effort is quantified as the number of LUCAS points for which observations are available in 2012, 2015 and 2018 per unit area of arable land in 2018 by country as reported by EUROSTAT (2023).

	Country	Number of Arable land area	Number of LUCAS
		LUCAS points	(×1000 ha) points per agricultural area
	Montenegro	12	8	1.42
	Cyprus	84	94	0.89
	Slovenia	83	167	0.50
	Luxembourg	25	62	0.40
	Austria	544	1408	0.39
	Belgium	332	976	0.34
	Czechia	690	2069	0.33
	Italy	1833	6142	0.30
	Latvia	309	1075	0.29
	The Netherlands	332	1269	0.26
	Finland	565	2171	0.26
	Denmark	601	2347	0.26
	Germany	3007	11 772	0.26
	Poland	2448	10 908	0.22
	France	3988	19 092	0.21
	Sweden	507	2461	0.21
	Bulgaria	668	3289	0.20
	Estonia	118	661	0.18
	Spain	2140	12 396	0.17
	Ireland	79	476	0.17
	Lithuania	362	2209	0.16
	Greece	320	1972	0.16
	Portugal	169	1158	0.15
	Slovakia	178	1358	0.13
	Romania	1048	8946	0.12
	United Kingdom	749	7023	0.11
	Hungary	429	4228	0.10

https://doi.org/10.5194/essd-15-5651-2023Earth Syst. Sci. Data, 15, 5651-5666, 2023

Earth Syst. Sci.Data, 15, 5651-5666, 2023 https://doi.org/10.5194/essd-15-5651-2023

Data described in this paper can be accessed at https://doi.org/10.5281/zenodo.7016986 (Ballot et al., 2022).

Code availability

Code used to produce the data described in this paper as well as to create the figures and tables can be accessed at https://doi.org/10.5281/zenodo.7018283 (Ballot and Guilpart, 2022).