# Dataset about the adoption of winter cover crops at the municipality level for mainland France

Benjamin Nowak, Audrey A. Michaud, Gaëlle Marliac

<figure><img src="image_1.png" /></figure>

To cite this version:

Benjamin Nowak, Audrey A. Michaud, Gaëlle Marliac. Dataset about the adoption of winter cover crops at the municipality level for mainland France. Data in Brief, 2022, 44, pp.108544. 10.1016/j.dib.2022.108544. hal-03783670

HAL Id: hal-03783670

https://hal.inrae.fr/hal-03783670

Submitted on 22 Sep 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Data Article

# Dataset about the adoption of winter cover crops at the municipality level for mainland France

<figure><img src="image_4.png" /></figure>

Benjamin Nowak<sup>a,*</sup>, Audrey Michaud<sup>b</sup>, Gaëlle Marliac<sup>c</sup>

ARTICLE INFO

Article history:

Received 9 May 2022

Revised 10 August 2022

Accepted 11 August 2022

Available online 17 August 2022

Dataset link: Dataset about the adoption of winter cover crops at the municipality level for mainland France (Original data)

Keywords:

Cover crops

Nutrient leaching

Agriculture

Carbon sequestration

Soil erosion

Farming

ABSTRACT

Winter soil cover by vegetation is associated with multiple benefits, such as increasing soil carbon storage and reducing erosion and nutrient leaching. This dataset provides an estimate of winter soil cover before spring-sown crops at municipality level for mainland France for two years (2018 and 2019). These estimates were obtained through the monitoring of all plots with spring-sown crops, declared within the context of the European Common Agricultural Policy. Detection of plots with winter soil cover was achieved through the analysis of Normalized Difference Vegetation Index (NDVI) time series, computed from Sentinel-2 multispectral images. For this dataset, it was considered that soil cover had to exceed 50% for a plot to be considered as covered by vegetation. Based on the literature, this corresponds to a threshold NDVI value between 0.45 and 0.59. To allow for sensitivity and uncertainty analyses for future studies that may be conducted using these data, three estimates of winter soil cover are given: minimum (based on the cultivated area exceeding the upper NDVI threshold of 0.59), maximum (considering the lower NDVI threshold of 0.45) and best estimate (mean NDVI threshold of 0.52). This dataset may be useful primarily to researchers working on biogeochemical cycle modeling or to government agencies, as several public policies (such as the Nitrates Directive) aim at developing winter cover crops.

Dataset link: Dataset about the adoption of winter cover crops at the municipality level for mainland France (Original data)

<footnote><sup>a</sup> Université Clermont Auvergne, AgroParisTech, INRAE, VetAgro Sup, UMR Territoires, 63370, Lempdes, France</footnote>

<footnote><sup>b</sup> Université Clermont Auvergne, INRAE, VetAgro Sup, UMR Herbivores, 63122, Saint-Genès-Champanelle, France</footnote>

<footnote><sup>c</sup> Université Clermont Auvergne, INRAE, UMR GDEC, 63000, Clermont-Ferrand, France</footnote>

<footnote>* Corresponding author.</footnote>

<footnote>E-mail address: bjn.nowak@vetagro-sup.fr (B. Nowak).</footnote>
<table><caption>Specifications Table</caption><tr><td>Subject</td><td>Agronomy and Crop Science</td></tr><tr><td>Specific subject area</td><td>Evaluation of the adoption of winter cover crops before spring-sown for mainland France</td></tr><tr><td>Type of data</td><td>Table (csv format)</td></tr><tr><td>How the data were acquired</td><td>Large scale analysis of Normalized Difference Vegetation Index (NDVI) computed from Sentinel-2 multispectral images with the Google Earth Engine platform</td></tr><tr><td>Data format</td><td>Analyzed</td></tr><tr><td>Description of data collection</td><td>Selection of all plots with spring crops for 2018 and 2019 with the French Registre Parcellaire Graphique. Computation of Normalized Difference Vegetation Index (NDVI) time series for each plot. Detection of plots with winter soil cover using thresholds defined in the literature.</td></tr><tr><td>Data source location</td><td>Country scale: mainland France<br>Raw data sources used:<br>- French Registre Parcellaire Graphique to defined plot boundaries (available here: https://www.data.gouv.fr/fr/datasets/registre-parcellaire-graphique-rpg-contours-des-parcelles-et-ilots-culturaux-et-leur-groupe-de-cultures-majoritaire/)<br>- Sentinel-2 multispectral images to compute Normalized Difference Vegetation Index (NDVI) time series for each plot. (available here: https://scihub.copernicus.eu/dhus/#/home)</td></tr><tr><td>Data accessibility</td><td>Repository name: WinterSoilCover (https://github.com/BjnNowak/WinterSoilCover)<br>Data identification number (DOI): 10.5281/zenodo.6822857<br>Direct URL to data: https://raw.githubusercontent.com/BjnNowak/WinterSoilCover/main/dataset/dataset_winter_soilCover.csv<br>Direct URL to code: https://github.com/BjnNowak/WinterSoilCover/blob/main/code/gee_code_example<br>Instructions for accessing these data: None (open access related research article: B. Nowak, G. Marliac, A. Michaud, Estimation of winter soil cover by vegetation before spring-sown crops for mainland France using multispectral satellite imagery, Environ. Res. Lett. 16:064,024 (2021). https://doi.org/10.1088/1748-9326/ac007c)</td></tr><tr><td>Related research article</td><td>B. Nowak, G. Marliac, A. Michaud, Estimation of winter soil cover by vegetation before spring-sown crops for mainland France using multispectral satellite imagery, Environ. Res. Lett. 16:064,024 (2021). https://doi.org/10.1088/1748-9326/ac007c</td></tr></table>

Value of the Data

* Winter soil cover by vegetation is associated with multiple benefits, such as increasing soil carbon storage and reducing erosion and nutrient leaching.

* This dataset may be useful primarily to researchers working on biogeochemical cycle modeling or to government agencies.

* Estimating the adoption of winter cover crops can help assess the additional soil carbon storage provided by these crops.
* As the winter period is a critical period for leaching, the estimation of winter soil cover is essential for models aiming to evaluate the transfer of mineral elements to waterways.

* As several public policies (such as the Nitrates Directive) aim at developing winter cover crops, this dataset will allow to evaluate the efficiency of these policies.

# 1. Data Description

Dataset For each municipality and for each spring crop of mainland France, this dataset provides the total area occupied and an estimate of the area covered by vegetation during the winter prior to sowing the spring crop. The detailed description of the variables is listed below:

* YEAR: Spring crop harvest year (2018 or 2019).

* INSEE: INSEE code of the municipality.

* CODE_CULTU: Code identification of the following spring crop as given in the French Registre Parcellaire Graphique.

* FOLLOWING_CROP: Translation of the CODE_CULTU code in English.

* TOTAL_AREA: Total area (in hectares) cultivated with this spring crop.

* COVERED_AREA: "Best" estimation of the area with winter cover crops before sowing this spring crop (in hectares).

* COVERED_AREA_MIN: Minimum estimation of the area with winter cover crops before sowing this spring crop (in hectares).

* COVERED_AREA_MAX: Maximum estimation of the area with winter cover crops before sowing this spring crop (in hectares).

* DEP: Department number.

* REG: Region number.

# 2. Experimental Design, Materials and Methods

This dataset provides an estimate of winter soil cover before spring crops at municipality level for mainland France for two years (2018 and 2019). It has been used to evaluate the effect of crop rotation [1] or soil and climatic conditions [2] on the adoption of winter cover crops in France.

The procedure for assessing winter soil cover was as follows:

## 1. Definition of the borders of plots with spring crops

For France, the limits of the plots declared within the context of the European Common Agricultural Policy are stored in a vector file known as the "Registre Parcellaire Graphique" [3,4]. For each year, this file indicates the type of crop grown on each plot. It was thus possible to select only the plots with spring crops with the variable "CODE_CULTU" of the attributes table, which gives the identification of each crop with a three letter code. To facilitate future use of this dataset, a "FOLLOWING_CROP" column has been added here to give the detailed name of each crop in English. Compared to the original "Registre Parcellaire Graphique" files, a 20 m negative buffer was applied to the borders of each plot to avoid edge effects due to sensor resolution or ground geolocation uncertainty during the soil cover detection stage.

## 2. Estimation of winter soil cover using Sentinel-2 multi-spectral images

Soil cover rate during the winter prior to sowing spring crops was estimated from the computation of the Normalized Difference Vegetation Index (NDVI) from multi-spectral images. The NDVI is calculated as follows:

$$NDVI = \frac{\rho_{NIR} - \rho_{Red}}{\rho_{NIR} + \rho_{Red}} \qquad (1)$$

Where $\rho_{NIR}$ and $\rho_{Red}$ stand for near infrared and red reflectance, respectively.
<figure><img src="image_5.png" /><figcaption>**Fig. 1.** Estimation of soil cover by vegetation from the Normalized Difference Vegetation Index (NDVI). Based on the literature, a 50% soil cover corresponds to a threshold NDVI value between 0.45 and 0.59.</figcaption></figure>

Because of the way it is calculated (Eq. (1)), the NDVI values are between -1 and 1. A bare soil has an NDVI value of about 0.2, and this value will increase with the development of the vegetation cover (up to NDVI values that can reach 0.9 for crops with high biomass). The NDVI is especially appropriate for monitoring winter soil cover because soil cover has a linear relationship with NDVI (Fig. 1), whereas it is more difficult to distinguish between high and very high crop biomass because the NDVI tends to saturate once the soil is completely covered by vegetation.

Yet soil characteristics, such as color or moisture, and crop residues influence NDVI measurements, especially in the early stages of crop development, when the soil is still poorly covered. In order to limit these effects, it was therefore considered that soil cover had to exceed 50% for a plot to be considered as covered by vegetation.

The linear relationship between NDVI and soil cover is well established in the literature and the NDVI values for 50% soil cover are similar across species [5,6], such as wheat [7,8] or soybeans [9] for example. The values reported by the studies are also similar for any source of multispectral images, from field sensor [10] to satellite images [11,12]. Overall, the NDVI values corresponding to a soil cover of 50% are between 0.45 and 0.59 (Fig. 1). Thus, to allow for sensitivity and uncertainty analyses for future studies that may be conducted using these data, three estimates are given: best estimate (column "COVERED_AREA", corresponding to the NDVI threshold of 0.52), minimum ("COVERED_AREA_MIN", for a NDVI threshold of 0.59) and maximum ("COVERED_AREA_MAX", for a NDVI threshold of 0.45).

The NDVI computation was carried out through the Google Earth Engine platform [13], using Sentinel-2 multispectral images at 10 m spatial resolution [14], corrected to surface reflectance using Sen2Cor [15]. Two levels of filters were applied to remove invalid observations [16]. First, for the study period, only the least cloudy images (20% threshold) were selected. Then, a second filter at the pixel scale was applied to remove observations identified as clouds, shadows or snow (using the Scene Classification map provided with Sentinel 2 observations). An example summarizing the code used can be found in the GitHub repository associated with this dataset [16].
<figure><img src="image_6.png" /><figcaption>Fig. 2. Timeline for cover crop detection.</figcaption></figure>

For each year, winter soil cover monitoring was carried out for two months (December and January) during the winter before sowing the spring-sown crop (Fig. 2). December was chosen as the beginning of the study period to limit the risk of detecting unharvested spring-sown crops on the plots, such as grain maize or sugar beet that can be harvested late in the year. January was chosen as the end of the study period because some spring-sown crops, such as peas, can be sown as early as February. Furthermore, if a cover crop was present on the plot, it must have been already detected in December or January. For most plots, several NDVI values could be calculated for the study period. In this case, the maximum NDVI value was compared to the thresholds used to assess a soil cover rate above 50%.

Finally, results were aggregated at the municipality level. In the dataset, the “INSEE” column gives the INSEE code that identifies each municipality. For each crop and for each municipality, the total area occupied by the crop (in hectares) is specified, as well as three estimates (also in hectares) of the area covered, corresponding to the three NDVI thresholds defined above. For each crop and municipality combination, the given year corresponds to the year of harvest of the following spring crops (e.g. the year 2019 refers to the soil cover in December 2018 and January 2019).

## Ethics Statements

This work does not involve any type of human studies, animal studies, or data gathered using social media.

This manuscript adheres to ethics in publishing standards.

## Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

## Data Availability

Dataset about the adoption of winter cover crops at the municipality level for mainland France (Original data) (GitHub).

## CRediT Author Statement

**Benjamin Nowak:** Conceptualization, Methodology, Data curation, Writing – original draft, Writing – review & editing; **Audrey Michaud:** Conceptualization, Writing – review & editing; **Gaëlle Marliac:** Conceptualization, Writing – review & editing.
Acknowledgments

This research did not receive any specific grant from funding agencies in the public, commer-
cial, or not-for-profit sectors.

References

[1] B. Nowak, G. Marliac, A. Michaud, Estimation of winter soil cover by vegetation before spring-sown crops for mainland France using multispectral satellite imagery, Environ. Res. Lett. 16 (2021) 064024, doi:10.1088/1748-9326/ac007c.

[2] B. Nowak, A. Michaud, G. Marliac, Soil-climate factors have a greater influence on the presence of winter cover crops than regulatory constraints in France, Agron. Sustain. Dev. 42 (2022) 28, doi:10.1007/s13593-022-00770-y.

[3] data.gouv.fr, Registre parcellaire graphique (RPG) : contours des parcelles et îlots culturaux et leur groupe de cultures majoritaire, (n.d.). https://www.data.gouv.fr/fr/datasets/registre-parcellaire- graphique-rpg-contours-des- parcelles-et- îlots- culturaux- et-leur- groupe-de- cultures-majoritaire/. Accessed January, 2022.

[4] Géoservices I.G.N., Accéder au téléchargement des données libres IGN: https://geoservices.ign.fr/documentation/diffusion/telechargement-donnees-libres.html#rpg, (n.d.). https://geoservices.ign.fr/documentation/diffusion/telechargement-donnees-libres.html#rpg. Accessed March 24, 2021.

[5] J.C. Jiménez-Muñoz, J.A. Sobrino, A. Plaza, L. Guanter, J. Moreno, P. Martínez, Comparison between fractional vegetation cover retrievals from vegetation indices and spectral mixture analysis: case study of PROBA/CHRIS data over an agricultural area, Sensors 9 (2009) 768-793, doi:10.3390/s90200768.

[6] T.J. Trout, L.F. Johnson, J. Gartung, Remote sensing of canopy cover in horticultural crops, HortScience 43 (2008) 333-337, doi:10.21273/HORTSCI.43.2.333.

[7] S. Er-Raki, A. Chehbouni, N. Guemouria, B. Duchemin, J. Ezzahar, R. Hadria, Combining FAO-56 model and ground-
based remote sensing to estimate water consumptions of wheat crops in a semi-arid region, Agric. Water Manag.
87 (2007) 41–54, doi:10.1016/j.agwat.2006.02.004.

[8] R. López-Urrea, A. Montoro, J. González-Piqueras, P. López-Fuster, E. Fereres, Water use of spring wheat to raise water productivity, Agric. Water Manag. 96 (2009) 1305-1310, doi:10.1016/j.agwat.2009.04.015.

[9] A. de la Casa, G. Ovando, L. Bressanini, J. Martínez, G. Díaz, C. Miranda, Soybean crop coverage estimation from NDVI images with different spatial resolution to evaluate yield variability in a plot, ISPRS J. Photogramm. Remote Sens. 146 (2018) 531–547, doi:10.1016/j.isprsjprs.2018.10.018.

[10] K. Prabhakara, W.D. Hively, G.W. McCarty, Evaluating the relationship between biomass, percent groundcover and remote sensing indices across six winter cover crop fields in Maryland, United States, Int. J. Appl. Earth Obs. Geoinf. 39 (2015) 88–102, doi:10.1016/j.jag.2015.03.002.

[11] K. Imukova, J. Ingwersen, T. Streck, Determining the spatial and temporal dynamics of the green vegetation fraction of croplands using high-resolution RapidEye satellite images, Agric. For. Meteorol. 206 (2015) 113–123, doi:10.1016/j.agrformet.2015.03.003.

[12] L.F. Johnson, T.J. Trout, Satellite NDVI assisted monitoring of vegetable crop evapotranspiration in California's San Joaquin Valley, Remote Sens. 4 (2012) 439–455 (Basel), doi:10.3390/rs4020439.

[13] N. Gorelick, M. Hancher, M. Dixon, S. Ilyushchenko, D. Thau, R. Moore, Google earth engine: planetary-scale geospa-
tial analysis for everyone, Remote Sens. Environ. 202 (2017) 18–27, doi:10.1016/j.rse.2017.06.031.

[14] M. Drusch, U. Del Bello, S. Carlier, O. Colin, V. Fernandez, F. Gascon, B. Hoersch, C. Isola, P. Laberinti, P. Martimort, A. Meygret, F. Spoto, O. Sy, F. Marchese, P. Bargellini, Sentinel-2: ESA'S optical high-resolution mission for GMES operational services, Remote Sens. Environ. 120 (2012) 25–36, doi:10.1016/j.rse.2011.11.026.

[15] J. Louis, V. Debaecker, B. Pflug, M. Main-Korn, J. Bieniarz, U. Mueller-Wilm, E. Cadau, F. Gascon, L. Ouwehand, Sentinel-2 Sen2Cor: L2A processor for users, in: L. Ouwehand (Ed.), Proceedings of the Living Planet Sympo- sium 2016, Spacebooks Online, Prague, Czech Republic, 2016, pp. 1–8. http://esamultimedia.esa.int/multimedia/ publications/SP-740/SP-740_toc.pdf.. Accessed May 10, 2021.

[16] B. Nowak, A. Michaud, G. Marliac, Dataset about the adoption of winter cover crops at the municipality level for mainland France, (2022). doi:10.5281/ZENODO.6513633.