<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "https://jats.nlm.nih.gov/nlm-dtd/publishing/3.0/journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" xml:lang="en" dtd-version="3.0" article-type="research-article">
  <front>
    <journal-meta><journal-id journal-id-type="publisher">SOIL</journal-id><journal-title-group>
    <journal-title>SOIL</journal-title>
    <abbrev-journal-title abbrev-type="publisher">SOIL</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">SOIL</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">2199-398X</issn><publisher>
    <publisher-name>Copernicus Publications</publisher-name>
    <publisher-loc>Göttingen, Germany</publisher-loc>
  </publisher></journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5194/soil-12-715-2026</article-id><title-group><article-title>Plastic film residues on cropland: monitoring soil contamination through optical remote sensing</article-title><alt-title>Plastic film residues on cropland: monitoring soil contamination</alt-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Fabrizi</surname><given-names>Alessandro</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Fiener</surname><given-names>Peter</given-names></name>
          
        <ext-link>https://orcid.org/0000-0001-6244-4705</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff2">
          <name><surname>Van Oost</surname><given-names>Kristof</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Wilken</surname><given-names>Florian</given-names></name>
          <email>florian.wilken@uni-a.de</email>
        </contrib>
        <aff id="aff1"><label>1</label><institution>Institute of Geography, University of Augsburg, Augsburg, Germany</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>Earth and Life Institute, Université Catholique de Louvain, Louvain-la-Neuve, Belgium</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Florian Wilken (florian.wilken@uni-a.de)</corresp></author-notes><pub-date><day>30</day><month>June</month><year>2026</year></pub-date>
      
      <volume>12</volume>
      <issue>1</issue>
      <fpage>715</fpage><lpage>731</lpage>
      <history>
        <date date-type="received"><day>5</day><month>August</month><year>2025</year></date>
           <date date-type="rev-request"><day>1</day><month>December</month><year>2025</year></date>
           <date date-type="rev-recd"><day>28</day><month>May</month><year>2026</year></date>
           <date date-type="accepted"><day>8</day><month>June</month><year>2026</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2026 Alessandro Fabrizi et al.</copyright-statement>
        <copyright-year>2026</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026.html">This article is available from https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026.html</self-uri><self-uri xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026.pdf">The full text article is available as a PDF file from https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026.pdf</self-uri>
      <abstract><title>Abstract</title>

      <p id="d2e113">Plastic films have been improving agricultural production and covering an increasing surface area of cropland in the last decades. Yet their use has been connected to the generation of plastic residues, potentially acting as a main secondary microplastic source in agricultural soils. Monitoring the generation of plastic film residues is crucial for identifying good management practices and assessing the risk of plastic use in agriculture. Remote sensing has been qualified as a valuable tool for monitoring macroplastic residues mainly on waters, while its use on agricultural soils is mostly unexplored. Our study combined proximal and remote sensing techniques to lay the foundations of UAV (Unmanned Aerial Vehicle) use for monitoring macroplastic film residues on cropland.</p>

      <p id="d2e116">Through proximal and UAV acquisitions of five-band multispectral data (i.e., blue, green, red, red edge, near infrared), we highlighted the potential of off-the-shelf miniaturised sensors and identified possible workflows for detecting macroplastic film residues. Our findings highlight a greater, or at least similar, efficacy of spatial resolution compared to spectral resolution, encouraging the use of high-resolution RGB cameras over multispectral cameras. Through proximal acquisitions of hyperspectral data on 8 different agricultural plastic films, we built spectral libraries and located absorption features at wavelengths that are not commonly available on off-the-shelf multispectral cameras. We highlighted that these absorption features unambiguously identify plastic films on cropland and offer the potential to distinguish plastic types, encouraging the development of sensors tailored for plastic detection.</p>
  </abstract>
    
<funding-group>
<award-group id="gs1">
<funding-source>H2020 Marie Skłodowska-Curie Actions</funding-source>
<award-id>955334</award-id>
</award-group>
</funding-group>
</article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d2e128">Soils are estimated to be a greater plastic sink than the ocean surface (Kedzierski et al., 2023; Lofty et al., 2022) and to potentially act as a plastic source for other environmental matrices through water and wind erosion (Rehm et al., 2021; Rezaei et al., 2022). In experimental setups, high to very high concentrations of plastic have been shown to degrade soil physicochemical properties and affect plant health (Zhang et al., 2022), but the threshold concentrations at which these effects occur are still poorly quantified. Most likely, the source and physicochemical properties of the plastic are important drivers (Wang et al., 2023a, b; Liu et al., 2023), highlighting the necessity of tracing the origin of plastic contamination. The knowledge about the effects of current plastic contamination is still limited (Landrigan et al., 2023), while the ubiquitous presence, persistence, bioaccumulation, and possible unknown effects of plastic threaten soil health and call for plastic contamination monitoring plans.</p>
      <p id="d2e131">Plastic in soil may span a wide range of polymers, shapes, sizes, and sources. One of the major sources and sinks of plastic are agricultural soils (Horton et al., 2017; Kawecki and Nowack, 2019), since plastic is extensively used in various ways to enhance productivity, such as mulching and greenhouse films, or seed and fertiliser coatings (FAO, 2021; EIP-AGRI Focus Group, 2021). Around 80 % of plastic use in European vegetable production is attributed to plastic films used as crop covers (Agriculture Plastic Environment Europe, 2019). While plastic films have been associated with the generation of macro- and microplastic residues (Steinmetz et al., 2022) (<inline-formula><mml:math id="M1" display="inline"><mml:mrow><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">5</mml:mn></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M2" display="inline"><mml:mrow><mml:mo>&lt;</mml:mo><mml:mn mathvariant="normal">5</mml:mn></mml:mrow></mml:math></inline-formula> mm, respectively; Thompson et al., 2004; Frias and Nash, 2019), their use can increase crop production and reduce the use of water and agrochemicals (Lamont, 2017; Espí et al., 2006). Mulching films can reduce evaporation from the soil surface and hence reduce the demand for water in dry regions or increase soil temperature in cold regions to improve the control of harvesting time. Additionally, mulching films can reduce weed growth and can be used for soil disinfestation to reduce the use of agrochemicals. Greenhouse films typically create a controlled environment for the growth of fruits and vegetables outside of typical growing seasons. Quantifying the residues generated after plastic film use and establishing monitoring tools is fundamental not only for assessing the environmental risk of plastic contamination, but also for identifying good management practices, aiming at a sustainable use of plastic films.</p>
      <p id="d2e154">Plastic covers used in agriculture are mainly LDPE (low-density polyethylene) films (EIP-AGRI Focus Group, 2021) with thickness ranging from 15 to 200 microns and lifetime ranging from a single cropping season to multiple years (e.g., 4 years) (Scarascia-Mugnozza et al., 2011). Thicker films are generally more resistant to mechanical stress and less likely to produce residues (Steinmetz et al., 2022). The usage beyond lifetime and other management practices, like application and removal modes, may also affect the amount of residues generated (FAO, 2021; EIP-AGRI Focus Group, 2021). Overall, earlier studies, such as the ones cited above, improved our understanding of plastic residue generation and concentration in soils, but quantitative analyses are still hampered by the methodologies used for mapping plastic residues. Plastic is typically quantified by combining soil sampling and ex situ analyses, making comprehensive data collection and analysis on multiple fields slow. However, large-scale data collection is needed to cover the variety of plastic films used and management practices adopted.</p>
      <p id="d2e157">The broad area coverage combined with very high spatial resolution offered by UAVs (Unmanned Aerial Vehicles) (Colomina and Molina, 2014) potentially enables data collection on multiple fields and detection of objects in the size range of macroplastic. UAV use has been explored to detect macroplastics on waters and coastal areas first and foremost (Veenstra and Churnside, 2012; Martin et al., 2018), while applications on land are still limited.</p>
      <p id="d2e161">Plastic materials are known to have unique absorption features in the SWIR (short-wave infrared) (Kühn et al., 2004), and hyperspectral imaging has already been used for detecting floating macroplastics (Garaba and Dierssen, 2018). On land, the use of SWIR absorption features has been investigated for non-destructive microplastic detection (Pahlawan et al., 2025). Hyperspectral imaging has been used for the detection of plastic-covered surfaces (Zhou et al., 2021; Garaba, 2025; Despini et al., 2025) and pushed up to the identification of their functionality (e.g., greenhouses or photovoltaic panels) (Zhou et al., 2023). Bonifazi et al. (2023) have first shown the potential of using absorption features and hyperspectral imaging for the ex-situ recognition of macroplastic residues found on agricultural soils. Still, reference spectral libraries are currently not available for agricultural plastics, and their establishment has been shown to push forward remote sensing algorithms for plastic detection (Garaba et al., 2018, 2021). Despite the high potential of SWIR-based plastic detection, it must be considered that the absorption bands are typically available on hyperspectral sensors only, and their limited availability may hamper large-scale monitoring. Multispectral cameras offer a good compromise between data volume and spectral resolution, and the availability of off-the-shelf miniaturised sensors has favoured their use on low-payload UAVs. To the best of our knowledge, multispectral cameras and the non-visible part of the spectrum in general have not been explored for remote sensing detection of plastic residues on agricultural soils. RGB cameras further decrease the complexity of data collection and analysis, as additional data processing, such as radiometric calibration, which can be more or less challenging depending on light conditions, is not required. Moreover, RGB cameras increase spatial resolution at the cost of spectral resolution. This allows for the acquisition of data with a higher ground sampling distance at the same flight height, or to fly higher and cover larger areas in shorter times to acquire data with the same ground sampling distance. The use of RGB cameras, coupled with deep learning, has been shown to result in reasonable plastic film residue mapping on cropland, but it has been limited to the detection of transparent film residues on cotton fields (Zhai et al., 2022; Qiu et al., 2022; Yang et al., 2024). Different film and soil types may require seeing beyond the visible.</p>
      <p id="d2e164">Our work aimed at a comprehensive assessment of UAV optical remote sensing use to detect different plastic film residues on cropland. We combined hyperspectral and multispectral proximal sensing with multispectral remote sensing to (i) build spectral libraries for the most common plastic films used in agriculture; (ii) define spectral changes occurring when moving from pristine plastic film to residue in soil; (iii) define possible workflows and sensors for detecting different plastic film residues using UAV technology.</p>

      <fig id="F1" specific-use="star"><label>Figure 1</label><caption><p id="d2e169">Experimental setup used for acquiring reflectance measurements of plastic films from a spectroradiometer (i.e., ASD) <bold>(a)</bold> and a multispectral camera <bold>(b)</bold>, and workflow used <bold>(c)</bold>. The same batch of samples is represented on the left <bold>(a)</bold> and on the right part of the image <bold>(b)</bold>. On the left side <bold>(a)</bold>, columns were labelled according to the treatment; different rows show different film types, from top to bottom: Black – biodegradable, Black – type 1, Black&amp;White – small (black side), Black&amp;White – small (white side), Transparent – large, Transparent – with holes. Film names refer to Table 1.</p></caption>
        <graphic xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026-f01.jpg"/>

      </fig>

</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Methods</title>
<sec id="Ch1.S2.SS1">
  <label>2.1</label><title>Proximal sensing</title>
<sec id="Ch1.S2.SS1.SSS1">
  <label>2.1.1</label><title>Experimental site and measuring devices</title>
      <p id="d2e219">Outdoor measurements of plastic film reflectance were carried out on a rooftop in Louvain-la-Neuve, Belgium, between 10:00 and 18:00 LT on 24 July 2022, with almost clear sky conditions. The measurements were performed through the parallel use of a spectroradiometer (ASD fieldspec 3; Fig. 1a; Malvern Panalytical, Ltd.), operated via an acquisition software in reflectance mode (RS<sup>3</sup>; Malvern Panalytical, Ltd.), and a multispectral camera (Micasense RedEdge-MX RX02; Fig. 1b; AgEagle Aerial Systems, Inc.), using a setup similar to Crucil et al. (2019).</p>
      <p id="d2e231">The spectroradiometer was used with a fibreoptic cable mounted on a pistol grip (Fig. 1a), which does not affect the 25° field of view of the bare fibreoptic cable. The device acquired hyperspectral data with 3 nm spectral resolution in the 350–1000 nm region, and 10 nm spectral resolution in the 1000–2500 nm region. The software automatically performed the conversion to reflectance measurements and output data at a 1 nm interval. A Spectralon panel with <inline-formula><mml:math id="M4" display="inline"><mml:mrow><mml:mn mathvariant="normal">12.7</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">12.7</mml:mn></mml:mrow></mml:math></inline-formula> cm<sup>2</sup> of reflective area (Fig. 1a), provided by the same company manufacturing the spectroradiometer, was used as a white reference to optimise and calibrate the instrument about every 4 min, and all the spectra were obtained from the average of 75 spectra, each with 34 ms of integration time. Thus, around 2.5 s were needed to perform each measurement. The fibreoptic tip was placed 11 cm above the targets at the nadir position, providing a shadow-free conical field of view with a diameter of 4.9 cm.</p>
      <p id="d2e255">The multispectral camera used five imaging sensors to acquire spectral data centred at 475, 560, 668, 717, 842 nm wavelengths – namely blue, green, red, red edge, and NIR (near infrared) bands – with bandwidths of 32, 27, 16, 12, 57 nm, respectively. Images of a reflectance panel were acquired before every measurement and used to calibrate and convert data to reflectance values in a later stage. The reflectance panel used had the following reflectance values for RedEdge-MX cameras with serial numbers RX02: blue <inline-formula><mml:math id="M6" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.49 %; green <inline-formula><mml:math id="M7" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.73 %; red <inline-formula><mml:math id="M8" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.81 %; red edge <inline-formula><mml:math id="M9" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.82 %; near infrared <inline-formula><mml:math id="M10" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.72 %. The camera was placed 70 cm above the targets at the nadir position, providing around 0.5 mm per pixel images.</p>
</sec>
<sec id="Ch1.S2.SS1.SSS2">
  <label>2.1.2</label><title>Plastic films, soil backgrounds, and experimental workflow</title>
      <p id="d2e301">The measurements were performed on eight different LDPE films (Table 1), where double-sided black and white films were used on both sides, and virtually increased the number of films to 10. Each plastic film was analysed in four different conditions: pristine, crumpled, dirty, and crumpled and dirty. Crumpling was obtained by manually rubbing the films to obtain a homogenous and randomised crease pattern (Fig. 1a), while dirtiness was obtained by rolling the films into a dry soil volume, thus creating a layer of soil particles attached to the film (Fig. 1a). A field soil (i.e., directly collected from a site), coming from arable land in the Belgian loam belt, where soils are characterised by an organic carbon content of <inline-formula><mml:math id="M11" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.26</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.37</mml:mn></mml:mrow></mml:math></inline-formula> % C and classified as silt loam according to the USDA classification (Zhang et al., 2021), and a standard soil (LUFA 2.2), characterised by organic carbon content of <inline-formula><mml:math id="M12" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.72</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.54</mml:mn></mml:mrow></mml:math></inline-formula> % C and classified as sandy loam according to USDA classification, were used for the experiment. The two soils were used to make the films dirty and as backgrounds for the measurements on plastic films. Both soils were sieved down to 250 <inline-formula><mml:math id="M13" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi></mml:mrow></mml:math></inline-formula>m, to remove coarser organic and inorganic matter, enabling the acquisition of stable soil spectra and ensuring that the pristine film's surface would not be influenced by soil roughness. Around 40 g of soil were placed into 87 mm petri dishes. The plastic films were cut to fit into the petri dishes, then the four treatments were applied, and finally, the films were placed on the soil surface in the petri dishes (Fig. 1). Additionally, some petri dishes were left with soil only and used as a reference spectrum of the soil backgrounds (Fig. 1). The use of black nitrile rubber O-rings (Fig. 1) with 85 mm of outer diameter and 75 mm of inner diameter was necessary to fix the films in the petri dishes and prevent wind from blowing them away.</p>

<table-wrap id="T1" specific-use="star"><label>Table 1</label><caption><p id="d2e339">Specifications of the LDPE plastic films used for the experiment. Thickness was obtained from the average of 10 measurements performed with a micrometre. The characteristics of typical use in plasticulture practices are obtained from the authors' knowledge and information provided by the manufacturing companies.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="6">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="center"/>
     <oasis:colspec colnum="3" colname="col3" align="left" colsep="1"/>
     <oasis:colspec colnum="4" colname="col4" align="left"/>
     <oasis:colspec colnum="5" colname="col5" align="left"/>
     <oasis:colspec colnum="6" colname="col6" align="left"/>
     <oasis:thead>
       <oasis:row>

         <oasis:entry colname="col1"/>

         <oasis:entry colname="col2"/>

         <oasis:entry colname="col3"/>

         <oasis:entry namest="col4" nameend="col6" align="center">Characteristics of typical use in </oasis:entry>

       </oasis:row>
       <oasis:row rowsep="1">

         <oasis:entry colname="col1"/>

         <oasis:entry namest="col2" nameend="col3" colsep="1">Physical characteristics </oasis:entry>

         <oasis:entry namest="col4" nameend="col6" align="center">plasticulture practices </oasis:entry>

       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>

         <oasis:entry colname="col1">Name</oasis:entry>

         <oasis:entry colname="col2">Thickness</oasis:entry>

         <oasis:entry colname="col3">Colour</oasis:entry>

         <oasis:entry colname="col4">Application</oasis:entry>

         <oasis:entry colname="col5">Crop</oasis:entry>

         <oasis:entry colname="col6">Duration</oasis:entry>

       </oasis:row>
       <oasis:row rowsep="1">

         <oasis:entry colname="col1"/>

         <oasis:entry colname="col2">(<inline-formula><mml:math id="M14" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi></mml:mrow></mml:math></inline-formula>m)</oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4">mode</oasis:entry>

         <oasis:entry colname="col5">type</oasis:entry>

         <oasis:entry colname="col6">use</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry rowsep="1" colname="col1">Black – type 1</oasis:entry>

         <oasis:entry rowsep="1" colname="col2"><inline-formula><mml:math id="M15" display="inline"><mml:mrow><mml:mn mathvariant="normal">18</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4"/>

         <oasis:entry colname="col5">vegetables</oasis:entry>

         <oasis:entry colname="col6">one</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry rowsep="1" colname="col1">Black – type 2</oasis:entry>

         <oasis:entry rowsep="1" colname="col2"><inline-formula><mml:math id="M16" display="inline"><mml:mrow><mml:mn mathvariant="normal">18</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>

         <oasis:entry colname="col3">black</oasis:entry>

         <oasis:entry colname="col4"/>

         <oasis:entry colname="col5">or fruits</oasis:entry>

         <oasis:entry colname="col6">growing</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry rowsep="1" colname="col1">Black – biodegradable</oasis:entry>

         <oasis:entry rowsep="1" colname="col2"><inline-formula><mml:math id="M17" display="inline"><mml:mrow><mml:mn mathvariant="normal">15</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>

         <oasis:entry rowsep="1" colname="col3"/>

         <oasis:entry colname="col4"/>

         <oasis:entry rowsep="1" colname="col5"/>

         <oasis:entry rowsep="1" colname="col6">season</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col1">Black&amp;White – small</oasis:entry>

         <oasis:entry rowsep="1" colname="col2" morerows="3"><inline-formula><mml:math id="M18" display="inline"><mml:mrow><mml:mn mathvariant="normal">99</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4"/>

         <oasis:entry colname="col5"/>

         <oasis:entry colname="col6"/>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col1">(black side)</oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4">On the soil</oasis:entry>

         <oasis:entry colname="col5"/>

         <oasis:entry colname="col6"/>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col1">Black&amp;White – small</oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4">surface</oasis:entry>

         <oasis:entry colname="col5"/>

         <oasis:entry colname="col6"/>

       </oasis:row>
       <oasis:row>

         <oasis:entry rowsep="1" colname="col1">(white side)</oasis:entry>

         <oasis:entry colname="col3">black/white</oasis:entry>

         <oasis:entry colname="col4"/>

         <oasis:entry colname="col5">asparagus</oasis:entry>

         <oasis:entry colname="col6"/>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col1">Black&amp;White – large</oasis:entry>

         <oasis:entry rowsep="1" colname="col2" morerows="3"><inline-formula><mml:math id="M19" display="inline"><mml:mrow><mml:mn mathvariant="normal">143</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">5</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4"/>

         <oasis:entry colname="col5"/>

         <oasis:entry colname="col6"/>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col1">(black side)</oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4"/>

         <oasis:entry colname="col5"/>

         <oasis:entry colname="col6">multiple</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col1">Black&amp;White – large</oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4"/>

         <oasis:entry colname="col5"/>

         <oasis:entry colname="col6"/>

       </oasis:row>
       <oasis:row>

         <oasis:entry rowsep="1" colname="col1">(white side)</oasis:entry>

         <oasis:entry rowsep="1" colname="col3"/>

         <oasis:entry colname="col4"/>

         <oasis:entry rowsep="1" colname="col5"/>

         <oasis:entry colname="col6">growing</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry rowsep="1" colname="col1">Transparent – with holes</oasis:entry>

         <oasis:entry rowsep="1" colname="col2"><inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:mn mathvariant="normal">44</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry rowsep="1" colname="col4">large tunnel</oasis:entry>

         <oasis:entry colname="col5">vegetables</oasis:entry>

         <oasis:entry colname="col6">seasons</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry rowsep="1" colname="col1">Transparent – small</oasis:entry>

         <oasis:entry rowsep="1" colname="col2"><inline-formula><mml:math id="M21" display="inline"><mml:mrow><mml:mn mathvariant="normal">157</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>

         <oasis:entry colname="col3">transparent</oasis:entry>

         <oasis:entry colname="col4">film for</oasis:entry>

         <oasis:entry colname="col5">or fruits</oasis:entry>

         <oasis:entry colname="col6"/>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col1">Transparent – large</oasis:entry>

         <oasis:entry colname="col2"><inline-formula><mml:math id="M22" display="inline"><mml:mrow><mml:mn mathvariant="normal">185</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4">greenhouses</oasis:entry>

         <oasis:entry colname="col5"/>

         <oasis:entry colname="col6"/>

       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap>

      <p id="d2e768">The measurements were performed on a table covered with a black plastic film, and the samples were placed on a black painted wood board (Fig. 1) to minimise background reflectance. However, background measurements were not performed. The petri dishes and the board were marked to allow for positioning the samples always at the same angle, minimising the influence of sample displacement on optical phenomena. The measurements were carried out in batches of a maximum of 25 samples (Fig. 1). First, hyperspectral measurements were performed on individual samples (Fig. 1a). At the end of the batch, one multispectral image was acquired over all the samples in the batch (Fig. 1b) Overall, 82 samples were analysed, i.e., (10 films <inline-formula><mml:math id="M23" display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> 4 treatments <inline-formula><mml:math id="M24" display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> 2 soil backgrounds) <inline-formula><mml:math id="M25" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula> 2 soils. For each sample, five replicate measurements were carried out throughout the day.</p>
</sec>
<sec id="Ch1.S2.SS1.SSS3">
  <label>2.1.3</label><title>Data processing and analysis</title>
      <p id="d2e801">Hyperspectral data were processed in R through the package prospectr (Stevens and Ramirez-Lopez, 2025). First, data at wavelengths below 400 nm and above 2400 nm were removed because of excessive noise in the outdoor measurements (Crucil et al., 2019; Crucil and Van Oost, 2021). Then, a Savitzky-Golay smoothing (Savitzky and Golay, 1964) was applied with a second order polynomial fit and a window size of 11 nm to remove random noise that can arise from outdoor measurements. Finally, water vapour absorption regions, between 1350 and 1450 nm and between 1800 and 1950 nm, were removed. Two spectra, corresponding to one replicate of pristine Black&amp;White – small (black side) on sandy loam and one replicate of pristine Black&amp;White – large (black side) on sandy loam (Table 1), were removed by visual assessment because of excessive noise (Fig. S1 in the Supplement), probably due to rapidly changing light conditions. An additional splice correction step could be applied to remove jumps in the spectra caused by sensor transitions in the spectroradiometer at 1000 and 1800 nm. However, as these spectral regions were not analysed for fine-scale spectral variations, splice correction was not applied.</p>
      <p id="d2e805">The multispectral images were pre-processed for radiometric referencing, vignetting correction, black level compensation and conversion to reflectance, following the procedure given by the original manufacturer (MicaSense, 2025) and using the script provided by Crucil (2022). Following, mean reflectance values were extracted over circular areas of around 5 cm in diameter, centred in each petri dish. This was done to relate the reflectance of an imaging sensor to a non-imaging sensor having a field of view with 4.9 cm of diameter.</p>
      <p id="d2e808">Mean values and standard deviation of the five replicates were calculated for each sample (<inline-formula><mml:math id="M26" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula>) and for each wavelength (<inline-formula><mml:math id="M27" display="inline"><mml:mi mathvariant="italic">λ</mml:mi></mml:math></inline-formula>) to account for variability in plastic film reflectance under changing light conditions. Mean spectra are presented for pristine films only. After calculating the average mean (<inline-formula><mml:math id="M28" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">μ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) and the average standard deviation (<inline-formula><mml:math id="M29" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">σ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) over the whole spectrum, the following coefficient of variation (CV) was used to describe the variation of measurements between replicates of the same sample:

              <disp-formula id="Ch1.E1" content-type="numbered"><label>1</label><mml:math id="M30" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="normal">CV</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">σ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">μ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>

            To compare the variability across specific groups of film colours and treatments, further aggregated means of the coefficient of variation were also calculated.</p>
      <p id="d2e888">To further compare the treatments, the difference in mean reflectance (<inline-formula><mml:math id="M31" display="inline"><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover></mml:math></inline-formula>) between the pristine and the treated film were calculated. For each treatment and for each wavelength, the difference (<inline-formula><mml:math id="M32" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi mathvariant="italic">ρ</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) is expressed as:

              <disp-formula id="Ch1.E2" content-type="numbered"><label>2</label><mml:math id="M33" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="normal">Δ</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mrow><mml:mi mathvariant="normal">treatment</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">λ</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mrow><mml:mi mathvariant="normal">pristine</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">λ</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mrow><mml:mi mathvariant="normal">treatment</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">λ</mml:mi></mml:mrow></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>

            As similarities between spectra of the same film colour groups were observed, <inline-formula><mml:math id="M34" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi mathvariant="italic">ρ</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> were aggregated for black, white, and transparent films through mean values.</p>
      <p id="d2e976">Lastly, the following three plastic indexes were calculated for hyperspectral data only, according to their original equations:

                  <disp-formula specific-use="align"><mml:math id="M35" display="block"><mml:mtable displaystyle="true"><mml:mtr><mml:mtd><mml:mstyle class="stylechange" displaystyle="true"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mi mathvariant="normal">HI</mml:mi><mml:mi mathvariant="normal">_</mml:mi><mml:mn mathvariant="normal">1215</mml:mn><mml:mo>=</mml:mo><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1216</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1197</mml:mn><mml:mo>)</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mn mathvariant="normal">1235</mml:mn></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mn mathvariant="normal">1197</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mn mathvariant="normal">1235</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1197</mml:mn></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mn mathvariant="normal">1197</mml:mn></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mn mathvariant="normal">1216</mml:mn></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mo>(</mml:mo><mml:mi mathvariant="normal">Garaba</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">and</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi mathvariant="normal">Dierssen</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mn mathvariant="normal">2018</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi mathvariant="normal">HI</mml:mi><mml:mi mathvariant="normal">_</mml:mi><mml:mn mathvariant="normal">1732</mml:mn><mml:mo>=</mml:mo><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1729</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1705</mml:mn><mml:mo>)</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mn mathvariant="normal">1741</mml:mn></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mn mathvariant="normal">1705</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mn mathvariant="normal">1741</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1705</mml:mn></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mn mathvariant="normal">1705</mml:mn></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mn mathvariant="normal">1729</mml:mn></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle class="stylechange" displaystyle="true"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mo>(</mml:mo><mml:mi mathvariant="normal">K</mml:mi><mml:mtext>ü</mml:mtext><mml:mi mathvariant="normal">hn</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi mathvariant="normal">et</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi mathvariant="normal">al</mml:mi><mml:mo>.</mml:mo><mml:mo>,</mml:mo><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mn mathvariant="normal">2004</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mi mathvariant="normal">ND</mml:mi><mml:mi mathvariant="normal">_</mml:mi><mml:mn mathvariant="normal">1715</mml:mn><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mrow><mml:mn mathvariant="normal">1590</mml:mn><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">to</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mn mathvariant="normal">1630</mml:mn></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mrow><mml:mn mathvariant="normal">1695</mml:mn><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">to</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mn mathvariant="normal">1735</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mrow><mml:mn mathvariant="normal">1590</mml:mn><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">to</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mn mathvariant="normal">1630</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mrow><mml:mn mathvariant="normal">1695</mml:mn><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">to</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mn mathvariant="normal">1735</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle class="stylechange" displaystyle="true"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mo>(</mml:mo><mml:mi mathvariant="normal">Castagna</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">et</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi mathvariant="normal">al</mml:mi><mml:mo>.</mml:mo><mml:mo>,</mml:mo><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mn mathvariant="normal">2023</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>

            where <inline-formula><mml:math id="M36" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the reflectance at the wavelength <inline-formula><mml:math id="M37" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula>, and <inline-formula><mml:math id="M38" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover><mml:mrow><mml:mi>x</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">to</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the average reflectance from wavelength <inline-formula><mml:math id="M39" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> to wavelength <inline-formula><mml:math id="M40" display="inline"><mml:mi>y</mml:mi></mml:math></inline-formula>. HI (hydrocarbon indexes) are narrow-band indexes involving three bands, where the middle band represents the wavelength at the absorption feature, and the other two are approximately located before and after the absorption region. ND_1715 captures the same absorption region represented by HI_1732, but it is a broadband index calculated through a normalised difference between the reflectance before the absorption region and the reflectance at the absorption region. ND_1715 was calculated to assess better the potential of a multispectral imaging sensor on a moving platform, which may require broader bandwidths to increase signal-to-noise ratio. To differentiate index values of plastic films from other possible targets (i.e. soil and crop residues), we have extracted index values of soils and crop residues from open spectral libraries (Hively et al., 2021; Kokaly et al., 2017).</p>
</sec>
</sec>
<sec id="Ch1.S2.SS2">
  <label>2.2</label><title>Remote sensing</title>
<sec id="Ch1.S2.SS2.SSS1">
  <label>2.2.1</label><title>Test site and equipment</title>
      <p id="d2e1400">An agricultural field located in the province of Córdoba, Spain (37°59<sup>′</sup>42.30<sup>′′</sup> N, 4°27<sup>′</sup>40.41<sup>′′</sup> W) was selected to test UAV capabilities in detecting plastic film residues (Fig. 2). The field is located next to the Guadalquivir River. The typical soil at this location is a calcaric Luvisol soil according to FAO classification (Junta de Andalucía, 2005), with soil organic carbon content of <inline-formula><mml:math id="M45" display="inline"><mml:mrow><mml:mn mathvariant="normal">0.77</mml:mn><mml:mo>±</mml:mo><mml:mn mathvariant="normal">0.64</mml:mn></mml:mrow></mml:math></inline-formula> % C (Muñoz-Rojas et al., 2012). Despite management practices not including the use of plastic directly on the field, a few plastic items were found on the field, probably resulting from littering, plastic contamination from neighboring fields, or plastic parts of machines used on the field.</p>

      <fig id="F2" specific-use="star"><label>Figure 2</label><caption><p id="d2e1459">Location of the study area, distribution of plastic films in the study area <bold>(a)</bold>, field picture <bold>(b)</bold>, and examples of black <bold>(c)</bold>, transparent <bold>(d)</bold>, and white <bold>(e)</bold> films placed on the field. On the bottom, workflow of the UAV experiment <bold>(f)</bold>. The geometry of the study area was chosen to maximise the number of plastic films, that are represented by coloured dots <bold>(a)</bold>. On the upper-left corner, the orthomosaic was obtained from the UAV images.</p></caption>
            <graphic xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026-f02.jpg"/>

          </fig>

      <p id="d2e1490">UAV images were acquired with a Phantom 4 pro (SZ DJI Technology Co., Ltd.) equipped with a downwelling light sensor and the same multispectral camera used in the proximal sensing experiment (Micasense RedEdge-MX RX02; Fig. 1b; AgEagle Aerial Systems, Inc.). The images were acquired on 28 September 2022, between 11:30 and 12:00 LT, with clear sky conditions. Images of a reflectance panel were acquired before the flight and used for the radiometric referencing of the images in a later stage. The reflectance panel used had the following reflectance values for RedEdge-MX cameras with serial numbers RX02: blue <inline-formula><mml:math id="M46" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.49 %; green <inline-formula><mml:math id="M47" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.73 %; red <inline-formula><mml:math id="M48" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.81 %; red edge <inline-formula><mml:math id="M49" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.82 %; near infrared <inline-formula><mml:math id="M50" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 47.72 %. The flight height was set to 7 m (Fig. 2f). The camera used had <inline-formula><mml:math id="M51" display="inline"><mml:mrow><mml:mn mathvariant="normal">1280</mml:mn><mml:mo>×</mml:mo><mml:mn mathvariant="normal">960</mml:mn></mml:mrow></mml:math></inline-formula> pixels of image resolution, 4.8 cm of sensor width, and 5.4 cm of focal length thus providing images with a spatial resolution of around 0.5 cm per pixel at 7 m of flight height. Additionally, two Reach RS2 RTK GNSS (global navigation satellite system) receivers (Emlid Tech Kft.) were used for placing the plastic films on the field and acquiring their coordinates.</p>
</sec>
<sec id="Ch1.S2.SS2.SSS2">
  <label>2.2.2</label><title>Plastic films and experimental workflow</title>
      <p id="d2e1550">Three of the films used in the proximal sensing experiment were used for the remote sensing experiment. With reference to Table 1, the plastic films used were: Black – type 1, Black&amp;White – small (white side), and Transparent – with holes. Films were selected to have one black, one white, and one transparent film. The films were cut into 5 cm <inline-formula><mml:math id="M52" display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> 10 cm pieces, and two small holes on the long edges of the films were made to facilitate the penetration of metal spikes that were used to fix the plastic films into the soil. The plastic films were placed on the field using the GNSS receivers, following a random distribution of points previously generated with QGIS (Free Software Foundation, Inc.) (Fig. 2a and f). After fixing the films into the soil, the metal spikes were covered with soil to avoid their influence on film reflectance. This left a film surface of around 5 cm <inline-formula><mml:math id="M53" display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> 5 cm uncovered and detectable (Fig. 2c–e). In total, 21 black films, 19 transparent films, and 13 white films were placed on the field (Fig. 2a and f). The experiment was originally designed to place a higher number of films, equally distributed among film colours, on a larger area. However, it was not possible to acquire images of the entire area, and the study area was restricted to that shown in Fig. 2. This resulted in an uneven distribution of films among the different colours.</p>
</sec>
<sec id="Ch1.S2.SS2.SSS3">
  <label>2.2.3</label><title>Data processing and analysis</title>
      <p id="d2e1575">Pix4D mapper (Pix4D SA) was used to create orthophoto mosaics through a structure-from-motion algorithm. The same software was used to perform the conversion to reflectance measurements, using data from the downwelling light sensor and the reflectance panel. Reflectance of saturated pixels, acquired as no data by the multispectral camera, was set to the maximum reflectance of 1. The outputs of the pre-processing – five orthophoto mosaics, each representing one of the five bands of the multispectral camera – were used for further analyses in ArcGIS Pro (Environmental Systems Research Institute, Inc.).</p>
      <p id="d2e1578">The centre of each plastic film was found on the images and labelled according to film colour (Fig. 2f). Then, 312 random points were generated and labelled by visual interpretation. All the randomly generated points were labelled as “soil” and, together with the plastic films, represented the ground observations. The ground observations were used for training and testing a pixel-based image classification algorithm with a 5-fold cross-validation (Fig. 2f). At each fold, 80 % of the observations is used for training and 20 % of the observations is used for validation. A buffer of 1 cm was applied on training points to increase training features (Fig. 2f). The buffer was applied after the division between training and validation data, to avoid the selection of validation data within the buffer. The subclasses “soil”, “shadow”, and “other” – mainly representing vegetation residues or unidentified classes – were created after the classification to evaluate spectral features and allow a better interpretation of the results.</p>
      <p id="d2e1581">The classification was performed through a random forest (Breiman, 2001), with the maximum number of trees set to 50 and the maximum depth of each tree set to 30. After the classification, the images were post-processed using a majority filter on 4-connected pixels (Fig. 2f), increasing the minimum mapping unit to pixel neighbourhoods but minimising isolated pixel artefacts. Four different datasets were created to compare the results of different spectral and spatial resolutions (Fig. 2f): <list list-type="order"><list-item>
      <p id="d2e1586">Multispectral 0.5 cm: it contains the original five multispectral bands, plus the following vegetation index:<disp-formula specific-use="align"><mml:math id="M54" display="block"><mml:mtable displaystyle="true"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mi mathvariant="normal">NDVI</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mo>(</mml:mo><mml:mi mathvariant="normal">Normalized</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">Difference</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi mathvariant="normal">Vegetation</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">Index</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle class="stylechange" displaystyle="true"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="normal">NIR</mml:mi><mml:mo>-</mml:mo><mml:mi mathvariant="normal">Red</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">NIR</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Red</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mo>(</mml:mo><mml:mi mathvariant="normal">Tucker</mml:mi><mml:mo>,</mml:mo><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mn mathvariant="normal">1979</mml:mn><mml:mo>)</mml:mo><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula></p></list-item><list-item>
      <p id="d2e1669">Multispectral 1.2 cm: it contains the same six bands as the multispectral 0.5 cm dataset. The spatial resolution was reduced to 1.2 cm by resampling the raster through nearest neighbour. The resolution was chosen to simulate the resolution obtained from the same multispectral camera at 20 m of flight altitude. This flight altitude is needed to acquire RGB data at 0.5 cm of spatial resolution, following the specifications of the RGB camera mounted on a Phantom 4 Pro.</p></list-item><list-item>
      <p id="d2e1673">RGB 0.5 cm: it contains the three original visible bands (i.e., red, green, blue).</p></list-item><list-item>
      <p id="d2e1677">RGB <inline-formula><mml:math id="M55" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula> indexes 0.5 cm: in addition to the three original bands, three band ratios were calculated as follows:<disp-formula id="Ch1.E3" specific-use="align" content-type="subnumberedsingle"><mml:math id="M56" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E3.4"><mml:mtd><mml:mtext>3a</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi mathvariant="normal">Red</mml:mi><mml:mrow><mml:mi mathvariant="normal">Green</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Blue</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E3.5"><mml:mtd><mml:mtext>3b</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi mathvariant="normal">Green</mml:mi><mml:mrow><mml:mi mathvariant="normal">Red</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Blue</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E3.6"><mml:mtd><mml:mtext>3c</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi mathvariant="normal">Blue</mml:mi><mml:mrow><mml:mi mathvariant="normal">Red</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Green</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>The addition of the band ratios in dataset 4 was done to ensure that at least one RGB dataset had the same number of bands as the multispectral datasets, hence providing the same number of input features into the random forest algorithm.</p></list-item></list> The detection of the plastic films placed on the field was evaluated against the ground observations by calculating the producer accuracy as the ratio between the number of plastic films correctly detected and the number of plastic films placed on the field (Fig. 2f). The presence of false positives instead (i.e., pixels classified as plastic but not covered by plastic) was checked by comparing the actual and the expected number of plastic-covered pixels. First, the expected number of plastic-covered pixels was calculated and weighted by the producer accuracy to account for the performance of the classifier in detecting the known plastic films. Then, the ratio between the actual number of pixels classified as plastic and the expected number of plastic-covered pixels was used to calculate the potential overestimation of plastic on the field (Fig. 2f). For each film, the expected number of plastic-covered pixels (<inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi mathvariant="normal">pixels</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) was calculated as follows:

              <disp-formula id="Ch1.E7" content-type="numbered"><label>4</label><mml:math id="M58" display="block"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi mathvariant="normal">pixels</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi mathvariant="normal">films</mml:mi></mml:msub><mml:mo>×</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mi mathvariant="normal">film</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mi mathvariant="normal">pixel</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>×</mml:mo><mml:mi mathvariant="normal">PA</mml:mi></mml:mrow></mml:math></disp-formula>

            where <inline-formula><mml:math id="M59" display="inline"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi mathvariant="normal">films</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the number of films placed on the field, <inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mi mathvariant="normal">film</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the area of a single film, <inline-formula><mml:math id="M61" display="inline"><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mi mathvariant="normal">pixel</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the area of a single pixel, and PA is the producer accuracy of the classification. This approach was used to account for the low rate of plastic detection compared to soil, avoiding the inflation of user accuracy due to bias in the ground observations. An example of the confusion matrix for the multispectral 0.5 cm dataset is provided in Table S1 in the Supplement, but is not included in the results.</p>
</sec>
</sec>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Results</title>
<sec id="Ch1.S3.SS1">
  <label>3.1</label><title>Spectra of pristine films and plastic indexes</title>
      <p id="d2e1854">On both soils, the spectra of the films are divided into three groups (Fig. 3): black (Black – type 1, Black – type 2, Black – biodegradable, Black&amp;White – small (black side), Black&amp;White – large (black side)), white (Black&amp;White – small (white side), Black&amp;White – large (white side)), and transparent (Transparent – with holes, Transparent – small, Transparent – large). White and black film reflectance is not particularly influenced by the background soil (Fig. 3). In contrast, transparent film reflectance changes according to the soil, with thinner films being closer to the soil spectra than thicker films (Fig. 3). However, the distance between the plastic film and the soil spectra changes across the spectral region. For instance, Transparent – with holes is closer to the spectra of silt loam in the visible, while Transparent – small gets closer in the NIR and SWIR (Fig. 3). On sandy loam, Transparent – with holes and Transparent – small spectra are very close in the visible, while Transparent – large gets closer to the soil spectra in the NIR and SWIR (Fig. 3). These relationships could be due to different additives in the plastic films, targeting different uses and influencing the transmittance of the films. Transmittance measurements could further elucidate how plastic film transmittance changes across the spectral region, and could be a relevant topic for crop production (Jones et al., 2021).</p>

      <fig id="F3" specific-use="star"><label>Figure 3</label><caption><p id="d2e1859">Spectra of pristine plastic films and of soils used as background. On the left side, spectra acquired on silt loam; on the right side, spectra acquired on sandy loam. Film colours are represented by different colours, and variations within film colours are represented by line shapes. Regions of absorption features are highlighted in light grey. Film names refer to Table 1.</p></caption>
          <graphic xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026-f03.png"/>

        </fig>

      <p id="d2e1868">Overall, the reflectance of pristine plastic films showed a high coefficient of variation of 40 % on average, while both soils had a similar coefficient of variation of around 4 %. The reflectance of black films is low across the entire spectrum, and the coefficient of variation is 71 % on average, ranging from 10 % for the Black – type 2 on silt loam to 135 % for the Black&amp;White – small (black side) on silt loam. White films have the highest mean reflectance and the lowest coefficient of variation, ranging from 6 % to 17 %. The two different white films present a similar spectral shape, being extremely reflective in the visible spectra, and decreasingly reflective through the SWIR, while the thicker film had higher reflectance across the whole spectra on both soils (Fig. 3). The transparent films generally follow the spectrum of the soil background, and transparent films on a less reflective soil result in lower reflectance of plastic films, and vice versa. The characteristic soil-shaped spectrum is interrupted by absorption features in the SWIR, which are shown by white films at the same wavelengths (Fig. 3). Specifically, the absorption features were found around the following wavelengths: 1215, 1730, 1765, 2312, and 2352 nm. As a result, HI_1215 and HI_1732 allow distinguishing all non-black films, in both pristine and treated conditions, from soils and mixtures, and from crop residues (Fig. 4a and b). While both indices guarantee a complete separability between plastic and non-plastic, HI_1732 values are one order of magnitude greater than HI_1215. ND_1715 has values similar to HI_1732, but the broader spectral interval results in outliers for soils and mixtures, represented by rare minerals (e.g., Erionite, Xenotime) overlaying with ND_1715 values of plastic (Fig. 4c). For both transparent and white films, the indexes generally have higher values for films with a higher thickness (Fig. 4).</p>

      <fig id="F4" specific-use="star"><label>Figure 4</label><caption><p id="d2e1874">Boxplots of plastic index values <bold>(a–c)</bold> for plastic films, soils, and crop residues. Index values include pristine, crumpled, dirty, and crumpled and dirty films. Values for all black films (i.e., Black – type 1, Black – type 2, Black&amp;White – small (black side), Black&amp;White – large (black side), Black – biodegradable) were aggregated, as they did not represent any relevant value. An axis break was used in <bold>(c)</bold> to enhance the visualisation of the plot. </p></caption>
          <graphic xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026-f04.png"/>

        </fig>

</sec>
<sec id="Ch1.S3.SS2">
  <label>3.2</label><title>Reflectance of plastic film as residues on soil</title>
      <p id="d2e1897">Compared to pristine films, crumpled and dirty films decreased the reflectance variability between replicates of hyperspectral measurements. The average coefficient of variation was reduced to 30 % for crumpled films, 9 % for crumpled and dirty films, and 9 % for dirty films. Overall, crumpling did not influence plastic film reflectance, as the changes in reflectance are within the range of the coefficient of variation for all the films (Figs. 5 and S2–S4). On the contrary, the presence of soil on the film surface influenced plastic film reflectance (Figs. 5 and S2–S4). This is particularly evident in the case of white films, where the reflectance of pristine films is substantially higher compared to dirty films in the visible spectra, while the differences between the spectra are reduced for higher wavelengths (Figs. 5 and S2–S4). Black film reflectance also changes when covered by soil, but the negative values of reflectance difference in Fig. 5 indicate that the presence of soil increased the reflectance of the plastic films. The differences between spectra increase for higher wavelengths, as soil reflectance increases (Fig. 5). Transparent films are the ones affected the least by treatments, with some changes induced by the presence of soil in the case of sandy loam only (Fig. 5). All these differences were consistently observed for the multispectral measurements in proximal sensing, despite generally higher coefficients of variation (e.g., 45 % for pristine films and 15 % for soils, on average).</p>

      <fig id="F5" specific-use="star"><label>Figure 5</label><caption><p id="d2e1902">Mean difference in reflectance between pristine films and treatments, presented for two different soils used as background in the proximal sensing experiment. Different data acquisition modes are represented by shapes, while different film colours are represented by the colour scale.</p></caption>
          <graphic xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026-f05.png"/>

        </fig>

</sec>
<sec id="Ch1.S3.SS3">
  <label>3.3</label><title>Plastic film detection from UAV</title>
      <p id="d2e1919">The spectral separability between films and soil observed in proximal sensing is confirmed by the multispectral acquisitions from the UAV (Fig. 6). White films have the highest spectral separability from non-plastic classes, especially in the bands of the visible spectra, while the reflectance slightly overlaps with other non-plastic elements when approaching the infrared (Fig. 6). Transparent film reflectance largely overlaps with the reflectance of soil and other non-plastic elements, while black films are quite distinct from soils, especially in the infrared bands, where soil reflectance increases (Fig. 6). Overall, all the 0.5 cm datasets had similar performances in detecting the plastic films placed on the field (Fig. 7a), while the multispectral 1.2 cm dataset had higher omission errors. The difference between the 1.2 cm dataset and the 0.5 cm datasets mainly consists in the missed detection of transparent films (Fig. 7a). White films placed on the field were detected with all four datasets in almost every fold, while black films had a generally higher and more fluctuating omission error (Fig. 7a).</p>

      <fig id="F6"><label>Figure 6</label><caption><p id="d2e1924">Reflectance of plastic films divided by colour (i.e., black, transparent, white), compared to soil reflectance, shadow, and other non-plastic objects found on the experimental site.</p></caption>
          <graphic xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026-f06.png"/>

        </fig>

      <fig id="F7" specific-use="star"><label>Figure 7</label><caption><p id="d2e1935">On the left side, accuracy (i.e., producer accuracy) of plastic detection with four different datasets <bold>(a)</bold>, separated by film and aggregated (i.e., “all”). Error bars represent the standard deviation of the producer accuracy within the 5-fold cross-validation; the numbers below the bars (<inline-formula><mml:math id="M62" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula>) represent the average number of points used for validation. On the right side, factor of plastic overestimation of four different datasets, separated by film and aggregated <bold>(b)</bold>. The factor of plastic overestimation was calculated as the ratio between the detected number and the expected number of plastic-covered pixels. Error bars represent the standard deviation of the overestimation within the 5-fold cross-validation.</p></caption>
          <graphic xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026-f07.png"/>

        </fig>

      <p id="d2e1958">Despite the good performances in detecting the plastic films placed on the field, all the datasets overestimated the presence of plastic on the field (Fig. 7b). Overall, plastic overestimation ranges from an average of 21 times for multispectral 1.2 cm dataset to 44 times for RGB 0.5 cm dataset. Black films encountered the highest overestimation compared to the other films, with RGB datasets having a higher overestimation compared to multispectral datasets (Fig. 7b). After black films, transparent films showed the second highest overestimation, finally followed by white films (Fig. 7b). The overestimation was generally related to the confusion of black films with shadows, white films with highly reflective pixels, and transparent films with soil (Fig. 8).</p>

      <fig id="F8" specific-use="star"><label>Figure 8</label><caption><p id="d2e1963">Detail of the true colour images <bold>(a–c)</bold> and associated classification results <bold>(d–f)</bold> obtained with the multispectral 0.5 cm dataset. The classified map is shown in panel <bold>(g)</bold>. In images <bold>(a–f)</bold>, the exact location of the plastic films is highlighted. In image <bold>(a)</bold>, the example of a highly reflective soil region is highlighted. The location of the plastic films placed on the field is highlighted in Fig. 2a.</p></caption>
          <graphic xlink:href="https://soil.copernicus.org/articles/12/715/2026/soil-12-715-2026-f08.png"/>

        </fig>

</sec>
</sec>
<sec id="Ch1.S4">
  <label>4</label><title>Discussion</title>
<sec id="Ch1.S4.SS1">
  <label>4.1</label><title>Spectral reflectance of plastic film residues</title>
      <p id="d2e2003">A spectral library of white, black, and transparent LDPE films was built, nearly covering all plastic film colours used in agriculture as crop covers. This work introduces an open-access spectral library for agricultural plastic films and explores non-visible regions of the electromagnetic spectrum for identifying plastic films as residues in agricultural soils. Within white and transparent films, differences related to thickness and transparency of the films mainly affected the mean reflectance, while the resulting spectral shapes were highly comparable (Fig. 3). Transparent film spectra are highly influenced by soil spectra (Fig. 3), and our experimental setup allowed us to acquire the resulting reflectance of soil and plastic films. An experimental setup such as the one used by Jones et al. (2021) would allow a complete characterisation of spectral radiative properties – reflectivity, transmissivity, and absorptivity – and provide insights into spectral changes of transparent films with different backgrounds. Differences between black films can hardly be observed, as spectra are flat and the reflectance is low, producing a low signal-to-noise ratio and high coefficient of variations (Fig. S5).</p>
      <p id="d2e2006">The mean reflectance of white and transparent films showed much higher variability across replicate measurements compared to soils, resulting in higher coefficients of variation. This is related to both the experimental design and the optical properties of the films. In fact, the measurements were acquired around 4 hours before and 4 h after solar noon, with solar zenith angle ranging from 32 to 59°, and solar azimuth angle ranging from 103 to 262° (SunCalc.Org., 2025). Moreover, the plastic films have a smooth surface, which leads away from the assumption of Lambertian reflector – typically adopted in remote sensing – and causes specular reflection to take over diffuse reflection phenomena (Goddijn-Murphy and Dufaur, 2018; Goddijn-Murphy et al., 2017). This is particularly evident in the case of pristine films, which had the highest coefficient of variation, while the presence of soil on the film surface increases the roughness of the films, favouring diffuse reflection and decreasing the variability of the reflectance between replicates. The implications for plastic residue detection are mostly positive, as particle deposition is expected on land, while specular reflection may be a limiting factor in other applications, such as the detection of clean plastic films with satellite data.</p>
      <p id="d2e2009">Despite a decrease in the coefficient of variation, we did not find any relevant influence of crumpling on film reflectance, while soil was the main driver of changes in reflectance between pristine and treated films (Fig. 5). Plastic films retained soil particles differently after the dirty treatment, and the comparability between films may be limited in our experimental setup. However, black films appeared to have higher retention of soil particles compared to white films (Fig. 1), while changes in plastic film reflectance were particularly evident for white films in the visible spectra, where the distance between the spectra of plastic films and soil is the highest (Figs. 3 and 5). The presence of soil did not have a particular influence on plastic absorption features, enabling plastic identification through spectral indexes, also on soil-covered plastic films (Fig. 4). However, as the amount of soil on the film increases, a progressive deterioration of plastic spectra is expected, until reaching a soil cover at which plastic films will not be detectable anymore. Identifying this threshold of soil cover will be helpful in delineating the boundaries of remote sensing use for detecting plastic film residues on cropland. Our experimental setup did not allow for controlling the amount of soil placed on the films, but aimed at representing field conditions, where soil particle retention of the film surface varies based on film properties. Moreover, indoor spectroscopy would allow exploiting a higher signal-to-noise ratio and better control over the amount of soil placed on the films compared to outdoor measurements, where wind can remove soil particles from the film surface. We also advise against using a contact probe in indoor spectroscopy, as we experienced plastic film surface melting in contact with the light source.</p>
</sec>
<sec id="Ch1.S4.SS2">
  <label>4.2</label><title>Use of multispectral and true colour UAV images for plastic film detection</title>
      <p id="d2e2021">Film spectra acquired from UAV presented similar spectral shapes to the proximal sensing acquisitions (Figs. 3 and 6). White films had very high reflectance for shorter wavelengths and decreasing reflectance for longer wavelengths, transparent films had reflectance similar to soil, and black films had a very low reflectance. However, it must be considered that the reflectance of white films in the UAV acquisitions was highly influenced by assigning the maximum reflectance of 1 to no data values, corresponding to saturated pixels. This particularly influenced the blue and green bands, where the number of saturated pixels was more than half, while the red and red edge bands had around 20 % and 10 % of saturated pixels, respectively, and the NIR had no saturated pixels.</p>
      <p id="d2e2024">Good performances were achieved in detecting true positives of plastic films (Fig. 7a). Despite the presence of other plastic residues on the field, plastic films were overestimated compared to the expected number of plastic-covered pixels (Fig. 7b). The results of our study might have been influenced by structural flaws of the survey, such as the absence of ground control points or the low overlap of images during flight, which can fall below the recommended threshold of 70 % when flying at 7 m above ground with the camera used. This could have resulted, for example, in band misalignment and halo artefacts that are visible in the RGB composite (Fig. 2a). However, general considerations can be made on how detection accuracy changes depending on the film colour and on the region of the spectrum used.</p>
      <p id="d2e2027">Specifically, black films had the highest overestimation, which was due to the presence of shadows on the field (Fig. 8a). In our study area, recent tillage practices induced soil clod formation and emphasised the impact of shadows (Fig. 2). However, surface microtopography is always expected on arable land, and we identify shadows as the main limiting factor in black film detection, similarly to previous works dealing with black plastic detection (Iordache et al., 2022; Shan et al., 2018). White films were slightly overestimated in correspondence with highly reflective soil regions (Fig. 8a). As with shadows, the topography of our study area might have emphasised the issue, creating soil surface angles that result in specular reflection and produce high reflectance values, comparable to those of white films. While techniques like brightness thresholding could improve the results of the classification (Iordache et al., 2022), their implementation and effectiveness may be site-specific, and a few measures could be adopted during data collection to address the challenges related to the illumination geometry. Possible solutions are multi-temporal flights, high-resolution DSM (Digital Surface Model), or flying with a homogenous cloud cover and the highest solar elevation angle possible.</p>
      <p id="d2e2030">While transparent films had a lower overestimation compared to black films (Fig. 7b), the source of error is related to a strong spectral overlay with soil in the visible and near infrared (Figs. 3 and 6), rather than a structural flaw of the survey. In this case, the only solution to substantially improve plastic detection may be limited to having sensors with bands located in the SWIR. The additional bands on the multispectral camera did not show a substantial improvement in transparent film detection, especially when compared with RGB<inline-formula><mml:math id="M63" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>indexes 0.5 cm dataset (Fig. 7). For black films, the datasets showed similar producer accuracies, but the multispectral datasets showed a lower overestimation compared to the RGB datasets (Fig. 7). However, the difference in spatial resolution must be considered, as coarser pixel size negatively affects both the areal estimates and the minimum mapping unit. At the same time, for equal pixel size, RGB cameras allow for increasing flight height, consequently reducing flight time and increasing area coverage. Overall, the results of our study suggest that the increased spatial resolution of an RGB camera should be favoured over the higher spectral resolution of multispectral cameras when aiming at the detection of macroplastic residues of around 5 cm.</p>
      <p id="d2e2041">Regarding the algorithm used and the experimental setup, it needs to be highlighted that using more sophisticated algorithms like deep learning and image segmentation techniques would likely increase detection accuracy. However, our experimental setup might have introduced bias when using textural features (e.g., fixed size and shape of plastic films, placement of plastic films under soil clods) and finding the best image classification algorithm goes beyond the scope of our study. Instead, a random forest allowed us to relate the results to spectral features, and to define which sensors and survey workflow are best suited for plastic films detection on soil. Lastly, a high surface occupation of plastic film residues in agricultural fields is rarely documented, and our experiment was designed to replicate a realistic scenario. Future studies could consider increasing the size of the observations related to plastic films to have more robust results.</p>
</sec>
<sec id="Ch1.S4.SS3">
  <label>4.3</label><title>A new generation of sensors for monitoring plastic on soil</title>
      <p id="d2e2052">We have shown potential workflows for detecting plastic film residues with RGB or multispectral cameras. However, the accuracy of these approaches will be site and film-specific. Machine learning approaches could be limited by the reduced size of training samples. At the same time, the absence of distinct spectral features on current sensor technology hampers the application of traditional methodologies such as maximum likelihood or threshold selection based on spectral indexes. A coarse-textured, clear soil may represent the optimal site for detecting black films, providing a good contrast between the film and the background, and reducing the complexity of the topography by limiting the formation of soil clods. Such soils are also more likely to induce sensor saturation and a lower contrast with clear films, making the detection of white films harder. On the contrary, a fine-textured dark soil may provide a good environment for white film detection, while black film detection would be limited by a decreased contrast and eventually hindered by the presence of soil clods and direct sunlight. Our study investigated a fraction of all possible soil–plastic film combinations, and further research is needed to comprehensively assess the capabilities of UAV-based detection with current technologies. However, plastic films do not share any absorption feature in the visible and near infrared (Fig. 3), and their identification with currently available multispectral broadband sensors is mostly colour-driven and influenced by spectral shape ambiguities (Garaba, 2025).</p>
      <p id="d2e2055">While the use of plastic absorption bands in the aquatic environment may be limited by water absorption features (Knaeps et al., 2021; Moshtaghi et al., 2021; Garaba and Dierssen, 2018), their use is extremely promising for soils. Plastic indexes enable unambiguous identification of plastic film residues on soil (Fig. 4), allowing their detection even when spectral features overlap with other regions of the spectrum, such as the visible range. We showcased the use of plastic indexes available in the literature, but we provided open access to the spectral library used in this study to promote the development of new spectral indexes and detection algorithms (Garaba, 2025). Traditional detection algorithms, such as threshold selection on spectral indexes, potentially enable overcoming the limitations of machine-learning-based methodologies related to the reduced size of training datasets. Within currently available indexes, even the broadband index ND_1715 resulted very effective in identifying plastic films from other elements that could possibly be found on cropland – except for some rare minerals that will not be found on arable land (Fig. 4c) – demonstrating the effectiveness of broad bands located at plastic absorption features and encouraging the development of multispectral cameras tailored to plastic detection.</p>
      <p id="d2e2058">Plastic residues rarely exceed 100 cm<sup>2</sup> on agricultural land (Piehl et al., 2018; Stefano and Pleissner, 2022; Wang et al., 2022), placing high requirements on spatial resolution, which is traded off with spectral resolution and areal coverage in sensor technology (Shaw and Burke, 2003). This strongly limits the use of technologies with large areal coverage, such as satellite data. In addition, the use of technologies with high spectral resolution on airborne platforms, such as hyperspectral sensors, would require low flight altitudes, resulting in increased flight time and associated costs, as well as the collection, storage, and processing of large volumes of data. The availability of broad bands related to SWIR absorption features on miniaturised sensors, together with high-resolution visible bands, would reduce costs, data volume, and complexity of current technologies, representing a milestone for monitoring any non-black plastic residue on land.</p>
      <p id="d2e2070">Additionally, according to Castagna et al. (2023), remote sensing offers the potential to distinguish at least three plastic types, exploiting the different locations of absorption features in the SWIR. PE, PP, and PVC – the most produced plastic types in Europe and the most used polymers in agriculture (Plastics Europe AISBL, 2024) – have similar absorption features, and their distinction might be challenging (Castagna et al., 2023). PET instead – the fourth most produced plastic type in Europe (Plastics Europe AISBL, 2024) – has different absorption features, and it can be distinguished from PE, PP, and PVC by using narrow-band indexes (Castagna et al., 2023). On agricultural land, this could mean distinguishing plastic residues generated by agricultural management from plastic residues generated by littering (e.g., PET bottles). Building spectral libraries and characterising the absorption features of the most common plastic residues is key to direct future efforts in developing new sensors, which we identify as a necessary technological development to monitor soil plastic contamination. The introduction of legal limits on plastic concentration in soil may not be far off, as initial regulatory examples begin to emerge (Meixner et al., 2020). This will likely drive demand for faster and standardised monitoring and increase the appeal of a sensor for on-soil plastic detection. Moreover, as macroplastic residues can fragment into microplastic (Yang et al., 2022; Song et al., 2017; Julienne et al., 2019), reliable macroplastic monitoring will also support better microplastic assessments.</p>
</sec>
</sec>
<sec id="Ch1.S5" sec-type="conclusions">
  <label>5</label><title>Conclusion</title>
      <p id="d2e2083">The introduction of miniaturised multispectral imaging sensors has boosted the use of UAV remote sensing. Compared to RGB cameras, multispectral sensors provide access to a few additional bands, typically at the cost of reduced spatial resolution and increased complexity of use. The additional bands available on commercial multispectral sensors are mostly designed for vegetation monitoring, while the detection of plastic residues still relies on the identification of their colour. Within current technologies, our results support using RGB cameras for monitoring plastic film residues on agricultural soils with UAVs, favouring spatial resolution over spectral resolution. However, the accuracy of these techniques will be highly dependent on the type of soil and on the colour of the residue.</p>
      <p id="d2e2086">The use of plastic absorption features provides an unambiguous identification of the plastic residue, especially on soils. Currently, these features are available only on expensive and complex sensors, hampering their use for monitoring large scales, at which remote sensing is most needed. Plastics are ubiquitous, persistent, and bioaccumulative contaminants that require monitoring plans, and the development of miniaturised sensors targeting plastic detection may represent a pivotal moment for assessing the risk of plastic contamination on land. We expect this work to further promote the exploration of remote sensing potential for monitoring soil plastic contamination and to encourage building and sharing spectral libraries of plastic residues commonly found on soils.</p>
</sec>

      
      </body>
    <back><notes notes-type="dataavailability"><title>Data availability</title>

      <p id="d2e2093">The spectral libraries used in this research are available at: <ext-link xlink:href="https://doi.org/10.5281/zenodo.14336252" ext-link-type="DOI">10.5281/zenodo.14336252</ext-link> (Fabrizi et al., 2024).</p>
  </notes><app-group>
        <supplementary-material position="anchor"><p id="d2e2099">The supplement related to this article is available online at <inline-supplementary-material xlink:href="https://doi.org/10.5194/soil-12-715-2026-supplement" xlink:title="pdf">https://doi.org/10.5194/soil-12-715-2026-supplement</inline-supplementary-material>.</p></supplementary-material>
        </app-group><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d2e2108">All authors: writing – review and editing, conceptualisation, methodology, resources. AF, FW, PF: visualisation. FW, PF: supervision, project administration, funding acquisition. AF: validation, investigation, software, formal analysis, data curation, writing – original draft.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d2e2114">At least one of the (co-)authors is a member of the editorial board of <italic>SOIL</italic>. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.</p>
  </notes><notes notes-type="disclaimer"><title>Disclaimer</title>

      <p id="d2e2123">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. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.</p>
  </notes><ack><title>Acknowledgements</title><p id="d2e2129">We thank the Institute of Sustainable Agriculture (IAS-CSIC, Córdoba) for the access provided to the research infrastructures. We thank Marco Bravin and Ana Carolina Cugler Moreira for their contribution to the data collection. We thank the members of the Water and Soil Resource Research group of the University of Augsburg for their constant and valuable feedback. Grammarly was used for automated grammar and language proofreading only.</p></ack><notes notes-type="financialsupport"><title>Financial support</title>

      <p id="d2e2134">This project has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement no. 955334.</p>
  </notes><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d2e2140">This paper was edited by Rafael Clemente and reviewed by two anonymous referees.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bib1"><label>1</label><mixed-citation>Agriculture Plastic Environment Europe: Statistics – plasticulture in Europe, <uri>https://apeeurope.eu/statistics/</uri> (last access: 22 December 2019), 2019.</mixed-citation></ref>
      <ref id="bib1.bib2"><label>2</label><mixed-citation>Bonifazi, G., Francesconi, E., Gasbarrone, R., Palmieri, R., and Serranti, S.: A preliminary study on the utilization of hyperspectral imaging for the on-soil recognition of plastic waste resulting from agricultural activities, Land, 12, <ext-link xlink:href="https://doi.org/10.3390/land12101934" ext-link-type="DOI">10.3390/land12101934</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib3"><label>3</label><mixed-citation>Breiman, L.: Random forests, Mach. Learn., 45, 5–32, <ext-link xlink:href="https://doi.org/10.1023/A:1010933404324" ext-link-type="DOI">10.1023/A:1010933404324</ext-link>, 2001.</mixed-citation></ref>
      <ref id="bib1.bib4"><label>4</label><mixed-citation>Castagna, A., Dierssen, H. M., Devriese, L. I., Everaert, G., Knaeps, E., and Sterckx, S.: Evaluation of historic and new detection algorithms for different types of plastics over land and water from hyperspectral data and imagery, Remote Sens. Environ., 298, <ext-link xlink:href="https://doi.org/10.1016/j.rse.2023.113834" ext-link-type="DOI">10.1016/j.rse.2023.113834</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib5"><label>5</label><mixed-citation>Colomina, I. and Molina, P.: Unmanned aerial systems for photogrammetry and remote sensing: a review, ISPRS J. Photogram. Remote Sens., 92, 79–97, <ext-link xlink:href="https://doi.org/10.1016/j.isprsjprs.2014.02.013" ext-link-type="DOI">10.1016/j.isprsjprs.2014.02.013</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bib6"><label>6</label><mixed-citation>Crucil, G.: Application of the radiometric calibration model and image correction: from raw Rededge images to calibrated reflectanc, GitHub [code], <uri>https://github.com/Debrisman/Micasense_Rededge_correction/blob/main/FUN_RededgeM_processing.R</uri> (last access: 15 February 2022), 2022.</mixed-citation></ref>
      <ref id="bib1.bib7"><label>7</label><mixed-citation>Crucil, G. and Van Oost, K.: Towards mapping of soil crust using multispectral imaging, Sensors, 21, <ext-link xlink:href="https://doi.org/10.3390/s21051850" ext-link-type="DOI">10.3390/s21051850</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib8"><label>8</label><mixed-citation>Crucil, G., Castaldi, F., Aldana-Jague, E., van Wesemael, B., Macdonald, A., and Van Oost, K.: Assessing the performance of UAS-compatible multispectral and hyperspectral sensors for soil organic carbon prediction, Sustainability, 11, <ext-link xlink:href="https://doi.org/10.3390/su11071889" ext-link-type="DOI">10.3390/su11071889</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib9"><label>9</label><mixed-citation>Despini, F., Costanzini, S., Parmeggiani, D., and la Cecilia, D.: Analysis of hyperspectral data to define multispectral sensor specifications for enhancing plastic greenhouse detection, Remote Sens. Appl. Soc. Environ., 101802, <ext-link xlink:href="https://doi.org/10.1016/j.rsase.2025.101802" ext-link-type="DOI">10.1016/j.rsase.2025.101802</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bib10"><label>10</label><mixed-citation>EIP-AGRI Focus Group: Reducing the plastic footprint of agriculture – final report, European Commission, <uri>https://eu-cap-network.ec.europa.eu/publications/eip-agri-focus-group-plastic-footprint-final-report_en</uri> (last access: 29 December 2021), 2021.</mixed-citation></ref>
      <ref id="bib1.bib11"><label>11</label><mixed-citation>Espí, E., Salmerón, A., Fontecha, A., García, Y., and Real, A. I.: Plastic films for agricultural applications, J. Plast. Film Sheet., 22, 85–102, <ext-link xlink:href="https://doi.org/10.1177/8756087906064220" ext-link-type="DOI">10.1177/8756087906064220</ext-link>, 2006.</mixed-citation></ref>
      <ref id="bib1.bib12"><label>12</label><mixed-citation>Fabrizi, A., Fiener, P., Van Oost, K., and Wilken, F.: Hyperspectral and multispectral reflectance of agricultural plastic films, Zenodo [data set], <ext-link xlink:href="https://doi.org/10.5281/zenodo.14336252" ext-link-type="DOI">10.5281/zenodo.14336252</ext-link>, 2024. </mixed-citation></ref>
      <ref id="bib1.bib13"><label>13</label><mixed-citation>FAO: Assessment of agricultural plastics and their sustainability. A call for action, <ext-link xlink:href="https://doi.org/10.4060/cb7856en" ext-link-type="DOI">10.4060/cb7856en</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib14"><label>14</label><mixed-citation>Frias, J. P. G. L. and Nash, R.: Microplastics: finding a consensus on the definition, Mar. Pollut. Bull., 138, 145–147, <ext-link xlink:href="https://doi.org/10.1016/j.marpolbul.2018.11.022" ext-link-type="DOI">10.1016/j.marpolbul.2018.11.022</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib15"><label>15</label><mixed-citation>Garaba, S. P.: Towards advanced mapping of plastic greenhouses from EMIT, EnMAP and PRISMA hyperspectral missions, Environ. Chall., 22, 101398, <ext-link xlink:href="https://doi.org/10.1016/j.envc.2025.101398" ext-link-type="DOI">10.1016/j.envc.2025.101398</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bib16"><label>16</label><mixed-citation>Garaba, S. P. and Dierssen, H. M.: An airborne remote sensing case study of synthetic hydrocarbon detection using short wave infrared absorption features identified from marine-harvested macro- and microplastics, Remote Sens. Environ., 205, 224–235, <ext-link xlink:href="https://doi.org/10.1016/j.rse.2017.11.023" ext-link-type="DOI">10.1016/j.rse.2017.11.023</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib17"><label>17</label><mixed-citation>Garaba, S. P., Aitken, J., Slat, B., Dierssen, H. M., Lebreton, L., Zielinski, O., and Reisser, J.: Sensing ocean plastics with an airborne hyperspectral shortwave infrared omager, Environ. Sci. Technol., 52, 11699–11707, <ext-link xlink:href="https://doi.org/10.1021/acs.est.8b02855" ext-link-type="DOI">10.1021/acs.est.8b02855</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib18"><label>18</label><mixed-citation>Garaba, S. P., Arias, M., Corradi, P., Harmel, T., de Vries, R., and Lebreton, L.: Concentration, anisotropic and apparent colour effects on optical reflectance properties of virgin and ocean-harvested plastics, J. Hazard. Mater., 406, 124290, <ext-link xlink:href="https://doi.org/10.1016/j.jhazmat.2020.124290" ext-link-type="DOI">10.1016/j.jhazmat.2020.124290</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib19"><label>19</label><mixed-citation>Goddijn-Murphy, L. and Dufaur, J.: Proof of concept for a model of light reflectance of plastics floating on natural waters, Mar. Pollut. Bull., 135, 1145–1157, <ext-link xlink:href="https://doi.org/10.1016/j.marpolbul.2018.08.044" ext-link-type="DOI">10.1016/j.marpolbul.2018.08.044</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib20"><label>20</label><mixed-citation>Goddijn-Murphy, L., Peters, S., van Sebille, E., James, N. A., and Gibb, S.: Concept for a hyperspectral remote sensing algorithm for floating marine macro plastics, Mar. Pollut. Bull., 126, 255–262, <ext-link xlink:href="https://doi.org/10.1016/j.marpolbul.2017.11.011" ext-link-type="DOI">10.1016/j.marpolbul.2017.11.011</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib21"><label>21</label><mixed-citation>Hively, W. D., Lamb, B. T., Daughtry, C. S. T., Quemada, M., Hunt Jr., R., Serbin, G., and Dennison, P.: Reflectance spectra of agricultural field conditions supporting remote-sensing evaluation of non-photosynthetic vegetation cover (ver. 1.1, November 2022), US Geological Survey data release [data set], <ext-link xlink:href="https://doi.org/10.5066/P9XK3867" ext-link-type="DOI">10.5066/P9XK3867</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib22"><label>22</label><mixed-citation>Horton, A. A., Walton, A., Spurgeon, D. J., Lahive, E., and Svendsen, C.: Microplastics in freshwater and terrestrial environments: Evaluating the current understanding to identify the knowledge gaps and future research priorities, Sci. Total Environ., 586, 127–141, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2017.01.190" ext-link-type="DOI">10.1016/j.scitotenv.2017.01.190</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib23"><label>23</label><mixed-citation>Iordache, M.-D., De Keukelaere, L., Moelans, R., Landuyt, L., Moshtaghi, M., Corradi, P., and Knaeps, E.: Targeting plastics: machine learning applied to litter detection in aerial multispectral images, Remote Sens., 14, <ext-link xlink:href="https://doi.org/10.3390/rs14225820" ext-link-type="DOI">10.3390/rs14225820</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib24"><label>24</label><mixed-citation>Jones, H., Black, T. A., Jassal, R. S., Nesic, Z., Johnson, M. S., and Smukler, S.: Characterization of shortwave and longwave properties of several plastic film mulches and their impact on the surface energy balance and soil temperature, Sol. Energy, 214, 457–470, <ext-link xlink:href="https://doi.org/10.1016/j.solener.2020.11.058" ext-link-type="DOI">10.1016/j.solener.2020.11.058</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib25"><label>25</label><mixed-citation>Julienne, F., Delorme, N., and Lagarde, F.: From macroplastics to microplastics: role of water in the fragmentation of polyethylene, Chemosphere, 236, 124409, <ext-link xlink:href="https://doi.org/10.1016/j.chemosphere.2019.124409" ext-link-type="DOI">10.1016/j.chemosphere.2019.124409</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib26"><label>26</label><mixed-citation>Junta de Andalucía: Mapa de Suelos de Andalucía, <uri>https://portalrediam.cica.es/VisorRediam/?lyr.add.wms=http://medioambiente/mapwms/REDIAM_Suelos_Andalucia</uri> (last access: 9 April 2025), 2005.</mixed-citation></ref>
      <ref id="bib1.bib27"><label>27</label><mixed-citation>Kawecki, D. and Nowack, B.: Polymer-specific modeling of the environmental emissions of seven commodity plastics as macro- and microplastics, Environ. Sci. Technol., 53, 9664–9676, <ext-link xlink:href="https://doi.org/10.1021/acs.est.9b02900" ext-link-type="DOI">10.1021/acs.est.9b02900</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib28"><label>28</label><mixed-citation>Kedzierski, M., Cirederf-Boulant, D., Palazot, M., Yvin, M., and Bruzaud, S.: Continents of plastics: an estimate of the stock of microplastics in agricultural soils, Sci. Total Environ., 880, 163294, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2023.163294" ext-link-type="DOI">10.1016/j.scitotenv.2023.163294</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib29"><label>29</label><mixed-citation>Knaeps, E., Sterckx, S., Strackx, G., Mijnendonckx, J., Moshtaghi, M., Garaba, S. P., and Meire, D.: Hyperspectral-reflectance dataset of dry, wet and submerged marine litter, Earth Syst. Sci. Data, 13, 713–730, <ext-link xlink:href="https://doi.org/10.5194/essd-13-713-2021" ext-link-type="DOI">10.5194/essd-13-713-2021</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib30"><label>30</label><mixed-citation>Kokaly, R. F., Clark, R. N., Swayze, G. A., Livo, K. E., Hoefen, T. M., Pearson, N. C., Wise, R. A., Benzel, W. M., Lowers, H. A., Driscoll, R. L., and Klein, A. J.: USGS Spectral Library Version 7 Data, US Geological Survey data release [data set], <ext-link xlink:href="https://doi.org/10.5066/F7RR1WDJ" ext-link-type="DOI">10.5066/F7RR1WDJ</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib31"><label>31</label><mixed-citation>Kühn, F., Oppermann, K., and Hörig, B.: Hydrocarbon Index – an algorithm for hyperspectral detection of hydrocarbons, Int. J. Remote Sens., 25, 2467–2473, <ext-link xlink:href="https://doi.org/10.1080/01431160310001642287" ext-link-type="DOI">10.1080/01431160310001642287</ext-link>, 2004.</mixed-citation></ref>
      <ref id="bib1.bib32"><label>32</label><mixed-citation>Lamont, W. J.: Plastic mulches for the production of vegetable crops, in: A guide to the manufacture, performance, and potential of plastics in agriculture, in: Plastic design library, edited by: Orzolek, M. D., Elsevier, 45–60, <ext-link xlink:href="https://doi.org/10.1016/b978-0-08-102170-5.00003-8" ext-link-type="DOI">10.1016/b978-0-08-102170-5.00003-8</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib33"><label>33</label><mixed-citation>Landrigan, P. J., Raps, H., Cropper, M., Bald, C., Brunner, M., Canonizado, E. M., Charles, D., Chiles, T. C., Donohue, M. J., Enck, J., Fenichel, P., Fleming, L. E., Ferrier-Pages, C., Fordham, R., Gozt, A., Griffin, C., Hahn, M. E., Haryanto, B., Hixson, R., Ianelli, H., James, B. D., Kumar, P., Laborde, A., Law, K. L., Martin, K., Mu, J., Mulders, Y., Mustapha, A., Niu, J., Pahl, S., Park, Y., Pedrotti, M.-L., Pitt, J. A., Ruchirawat, M., Seewoo, B. J., Spring, M., Stegeman, J. J., Suk, W., Symeonides, C., Takada, H., Thompson, R. C., Vicini, A., Wang, Z., Whitman, E., Wirth, D., Wolff, M., Yousuf, A. K., and Dunlop, S.: The Minderoo-Monaco commission on plastics and human health, Ann. Global Health, 89, 23, <ext-link xlink:href="https://doi.org/10.5334/aogh.4056" ext-link-type="DOI">10.5334/aogh.4056</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib34"><label>34</label><mixed-citation>Liu, M., Feng, J., Shen, Y., and Zhu, B.: Microplastics effects on soil biota are dependent on their properties: a meta-analysis, Soil Biol. Biochem., 178, 108940, <ext-link xlink:href="https://doi.org/10.1016/j.soilbio.2023.108940" ext-link-type="DOI">10.1016/j.soilbio.2023.108940</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib35"><label>35</label><mixed-citation>Lofty, J., Muhawenimana, V., Wilson, C. A. M. E., and Ouro, P.: Microplastics removal from a primary settler tank in a wastewater treatment plant and estimations of contamination onto European agricultural land via sewage sludge recycling, Environ. Pollut., 304, 119198, <ext-link xlink:href="https://doi.org/10.1016/j.envpol.2022.119198" ext-link-type="DOI">10.1016/j.envpol.2022.119198</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib36"><label>36</label><mixed-citation>Martin, C., Parkes, S., Zhang, Q., Zhang, X., McCabe, M. F., and Duarte, C. M.: Use of unmanned aerial vehicles for efficient beach litter monitoring, Mar. Pollut. Bull., 131, 662–673, <ext-link xlink:href="https://doi.org/10.1016/j.marpolbul.2018.04.045" ext-link-type="DOI">10.1016/j.marpolbul.2018.04.045</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib37"><label>37</label><mixed-citation>Meixner, K., Kubiczek, M., and Fritz, I.: Microplastic in soil–current status in Europe with special focus on method tests with Austrian samples, AIMS Environ. Sci., 7, 174–191, <ext-link xlink:href="https://doi.org/10.3934/environsci.2020011" ext-link-type="DOI">10.3934/environsci.2020011</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib38"><label>38</label><mixed-citation>MicaSense: Radiometric calibration model for MicaSense sensors, <ext-link xlink:href="https://support.micasense.com/hc/en-us/articles/115000351194-Radiometric-Calibration-Model-for-MicaSense-Sensors">https://support.micasense.com/hc/en-us/articles/115000351194-Radiometric-Calibration-Model-for-MicaSense- Sensors</ext-link> (last access: 11 April 2025), 2025.</mixed-citation></ref>
      <ref id="bib1.bib39"><label>39</label><mixed-citation>Moshtaghi, M., Knaeps, E., Sterckx, S., Garaba, S. P., and Meire, D.: Spectral reflectance of marine macroplastics in the VNIR and SWIR measured in a controlled environment, Sci. Rep., 11, 5436, <ext-link xlink:href="https://doi.org/10.1038/s41598-021-84867-6" ext-link-type="DOI">10.1038/s41598-021-84867-6</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib40"><label>40</label><mixed-citation>Muñoz-Rojas, M., Jordán, A., Zavala, L., De la Rosa, D., Abd-Elmabod, S., and Anaya-Romero, M.: Organic carbon stocks in Mediterranean soil types under different land uses (southern Spain), Solid Earth, 3, 375–386, <ext-link xlink:href="https://doi.org/10.5194/se-3-375-2012" ext-link-type="DOI">10.5194/se-3-375-2012</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bib41"><label>41</label><mixed-citation>Pahlawan, M. F. R., Kim, Y., Aline, U., Zahroh, A., Masithoh, R. E., Kim, M. S., Baek, I., and Cho, B.-K.: Non-destructive identification of microplastics in soil using spectroscopy and hyperspectral imaging, Trends Anal. Chem., 118216, <ext-link xlink:href="https://doi.org/10.1016/j.trac.2025.118216" ext-link-type="DOI">10.1016/j.trac.2025.118216</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bib42"><label>42</label><mixed-citation>Piehl, S., Leibner, A., Löder, M. G., Dris, R., Bogner, C., and Laforsch, C.: Identification and quantification of macro-and microplastics on an agricultural farmland, Sci. Rep., 8, 17950, <ext-link xlink:href="https://doi.org/10.1038/s41598-018-36172-y" ext-link-type="DOI">10.1038/s41598-018-36172-y</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib43"><label>43</label><mixed-citation>Plastics Europe AISBL: The circular economy for plastics – a European analysis, <uri>https://plasticseurope.org/wp-content/uploads/2024/03/CEreport_executivesummary_2024.pdf</uri> (last access: 11 April 2025), 2024.</mixed-citation></ref>
      <ref id="bib1.bib44"><label>44</label><mixed-citation>Qiu, F., Zhai, Z., Li, Y., Yang, J., Wang, H., and Zhang, R.: UAV imaging and deep learning based method for predicting residual film in cotton field plough layer, Front. Plant Sci., 13, 1010474, <ext-link xlink:href="https://doi.org/10.3389/fpls.2022.1010474" ext-link-type="DOI">10.3389/fpls.2022.1010474</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib45"><label>45</label><mixed-citation>Rehm, R., Zeyer, T., Schmidt, A., and Fiener, P.: Soil erosion as transport pathway of microplastic from agriculture soils to aquatic ecosystems, Sci. Total Environ., 795, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2021.148774" ext-link-type="DOI">10.1016/j.scitotenv.2021.148774</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib46"><label>46</label><mixed-citation>Rezaei, M., Abbasi, S., Pourmahmood, H., Oleszczuk, P., Ritsema, C., and Turner, A.: Microplastics in agricultural soils from a semi-arid region and their transport by wind erosion, Environ. Res., 212, <ext-link xlink:href="https://doi.org/10.1016/j.envres.2022.113213" ext-link-type="DOI">10.1016/j.envres.2022.113213</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib47"><label>47</label><mixed-citation>Savitzky, A. and Golay, M. J.: Smoothing and differentiation of data by simplified least squares procedures, Anal. Chem., 36, 1627–1639, <ext-link xlink:href="https://doi.org/10.1021/ac60214a047" ext-link-type="DOI">10.1021/ac60214a047</ext-link>, 1964.</mixed-citation></ref>
      <ref id="bib1.bib48"><label>48</label><mixed-citation>Scarascia-Mugnozza, G., Sica, C., and Russo, G.: Plastic materials in European agriculture: actual use and perspectives, J. Agric. Eng., 42, 15–28, <ext-link xlink:href="https://doi.org/10.4081/jae.2011.3.1" ext-link-type="DOI">10.4081/jae.2011.3.1</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bib49"><label>49</label><mixed-citation>Shan, J., Zhao, J., Liu, L., Zhang, Y., Wang, X., and Wu, F.: A novel way to rapidly monitor microplastics in soil by hyperspectral imaging technology and chemometrics, Environ. Pollut., 238, 121–129, <ext-link xlink:href="https://doi.org/10.1016/j.envpol.2018.03.026" ext-link-type="DOI">10.1016/j.envpol.2018.03.026</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib50"><label>50</label><mixed-citation> Shaw, G. A. and Burke, H. K.: Spectral imaging for remote sensing, Lincoln Lab. J., 14, 3–28, 2003.</mixed-citation></ref>
      <ref id="bib1.bib51"><label>51</label><mixed-citation>Song, Y. K., Hong, S. H., Jang, M., Han, G. M., Jung, S. W., and Shim, W. J.: Combined effects of UV exposure duration and mechanical abrasion on microplastic fragmentation by polymer type, Environ. Sci. Technol., 51, 4368–4376, <ext-link xlink:href="https://doi.org/10.1021/acs.est.6b06155" ext-link-type="DOI">10.1021/acs.est.6b06155</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib52"><label>52</label><mixed-citation>Stefano, N. and Pleissner, D.: Quantification and analysis of surface macroplastic contamination on arable areas, J. Soils Sediments, 22, 757–768, <ext-link xlink:href="https://doi.org/10.1007/s11368-021-03114-7" ext-link-type="DOI">10.1007/s11368-021-03114-7</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib53"><label>53</label><mixed-citation>Steinmetz, Z., Löffler, P., Eichhöfer, S., David, J., Muñoz, K., and Schaumann, G. E.: Are agricultural plastic covers a source of plastic debris in soil? A first screening study, Soil, 8, 31–47, <ext-link xlink:href="https://doi.org/10.5194/soil-8-31-2022" ext-link-type="DOI">10.5194/soil-8-31-2022</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib54"><label>54</label><mixed-citation>Stevens, A. and Ramirez-Lopez, L.: An introduction to the prospectr package, <uri>https://cran.r-project.org/web/packages/prospectr/vignettes/prospectr.html</uri> (last access: 10 April 2025), 2025.</mixed-citation></ref>
      <ref id="bib1.bib55"><label>55</label><mixed-citation>SunCalc.org.: <uri>https://www.suncalc.org/#/50.6693,4.62,17/2022.07.24/10:00/1/3</uri> (last access: 11 April 2025), 2025.</mixed-citation></ref>
      <ref id="bib1.bib56"><label>56</label><mixed-citation>Thompson, R. C., Olsen, Y., Mitchell, R. P., Davis, A., Rowland, S. J., John, A. W., McGonigle, D., and Russell, A. E.: Lost at sea: where is all the plastic?, Science, 304, 838, <ext-link xlink:href="https://doi.org/10.1126/science.1094559" ext-link-type="DOI">10.1126/science.1094559</ext-link>, 2004.</mixed-citation></ref>
      <ref id="bib1.bib57"><label>57</label><mixed-citation>Tucker, C. J.: Red and photographic infrared linear combinations for monitoring vegetation, Remote Sens. Environ., 8, 127–150, <ext-link xlink:href="https://doi.org/10.1016/0034-4257(79)90013-0" ext-link-type="DOI">10.1016/0034-4257(79)90013-0</ext-link>, 1979.</mixed-citation></ref>
      <ref id="bib1.bib58"><label>58</label><mixed-citation>Veenstra, T. S. and Churnside, J. H.: Airborne sensors for detecting large marine debris at sea, Mar. Pollut. Bull., 65, 63–68, <ext-link xlink:href="https://doi.org/10.1016/j.marpolbul.2010.11.018" ext-link-type="DOI">10.1016/j.marpolbul.2010.11.018</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bib59"><label>59</label><mixed-citation>Wang, P.-Y., Zhao, Z.-Y., Xiong, X.-B., Wang, N., Zhou, R., Zhang, Z.-M., Ding, F., Hao, M., Wang, S., Ma, Y., Uzamurera, A. G., Xiao, K.-W., Khan, A., Tao, X.-P., Wang, W.-Y., and Tao, H.-Y.: Microplastics affect soil bacterial community assembly more by their shapes rather than the concentrations, Water Res., 245, 120581, <ext-link xlink:href="https://doi.org/10.1016/j.watres.2023.120581" ext-link-type="DOI">10.1016/j.watres.2023.120581</ext-link>, 2023a.</mixed-citation></ref>
      <ref id="bib1.bib60"><label>60</label><mixed-citation>Wang, S., Fan, T., Cheng, W., Wang, L., Zhao, G., Li, S., Dang, Y., and Zhang, J.: Occurrence of macroplastic debris in the long-term plastic film-mulched agricultural soil: a case study of Northwest China, Sci. Total Environ., 831, 154881, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2022.154881" ext-link-type="DOI">10.1016/j.scitotenv.2022.154881</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib61"><label>61</label><mixed-citation>Wang, Y., Xiang, L., Amelung, W., Elsner, M., Gan, J., Kueppers, S., Christian, L., Jiang, X., Adu-Gyamfi, J., Heng, L., Sik Ok, Y., Ivleva, N. P., Luo, Y., Barceló, D., Schäffer, A., and Wang, F.: Micro-and nanoplastics in soil ecosystems: Analytical methods, fate, and effects, Trends Anal. Chem., 169, 117309, <ext-link xlink:href="https://doi.org/10.1016/j.trac.2023.117309" ext-link-type="DOI">10.1016/j.trac.2023.117309</ext-link>, 2023b. </mixed-citation></ref>
      <ref id="bib1.bib62"><label>62</label><mixed-citation>Yang, J., Zhai, Z., Li, Y., Duan, H., Cai, F., Lv, J., and Zhang, R.: Design and research of residual film pollution monitoring system based on UAV, Comput. Electron. Agr., 217, <ext-link xlink:href="https://doi.org/10.1016/j.compag.2023.108608" ext-link-type="DOI">10.1016/j.compag.2023.108608</ext-link>, 2024.</mixed-citation></ref>
      <ref id="bib1.bib63"><label>63</label><mixed-citation>Yang, Y., Li, Z., Yan, C., Chadwick, D., Jones, D. L., Liu, E., Liu, Q., Bai, R., and He, W.: Kinetics of microplastic generation from different types of mulch films in agricultural soil, Sci. Total Environ., 814, 152572, <ext-link xlink:href="https://doi.org/10.1016/j.scitotenv.2021.152572" ext-link-type="DOI">10.1016/j.scitotenv.2021.152572</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib64"><label>64</label><mixed-citation>Zhai, Z., Chen, X., Zhang, R., Qiu, F., Meng, Q., Yang, J., and Wang, H.: Evaluation of residual plastic film pollution in pre-sowing cotton field using UAV imaging and semantic segmentation, Front. Plant Sci., 13, 991191, <ext-link xlink:href="https://doi.org/10.3389/fpls.2022.991191" ext-link-type="DOI">10.3389/fpls.2022.991191</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib65"><label>65</label><mixed-citation>Zhang, H., Shi, P., Crucil, G., Wesemael, B., Limbourg, Q., and Van Oost, K.: Evaluating the capability of a UAV-borne spectrometer for soil organic carbon mapping in bare croplands, Land Degrad. Dev., 32, 4375–4389, <ext-link xlink:href="https://doi.org/10.1002/ldr.4043" ext-link-type="DOI">10.1002/ldr.4043</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib66"><label>66</label><mixed-citation>Zhang, J., Ren, S., Xu, W., Liang, C., Li, J., Zhang, H., Li, Y., Liu, X., Jones, D. L., Chadwick, D. R., Zhang, F., and Wang, K.: Effects of plastic residues and microplastics on soil ecosystems: a global meta-analysis, J. Hazard. Mater., 435, <ext-link xlink:href="https://doi.org/10.1016/j.jhazmat.2022.129065" ext-link-type="DOI">10.1016/j.jhazmat.2022.129065</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib67"><label>67</label><mixed-citation>Zhou, S., Kuester, T., Bochow, M., Bohn, N., Brell, M., and Kaufmann, H.: A knowledge-based, validated classifier for the identification of aliphatic and aromatic plastics by WorldView-3 satellite data, Remote Sens. Environ., 264, <ext-link xlink:href="https://doi.org/10.1016/j.rse.2021.112598" ext-link-type="DOI">10.1016/j.rse.2021.112598</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib68"><label>68</label><mixed-citation>Zhou, S., Mou, L., Hua, Y., Zhang, L., Kaufmann, H., and Zhu, X. X.: Can we use deep learning models to identify the functionality of plastics from space?, Int. J. Appl. Earth Obs. Geoinf., 123, <ext-link xlink:href="https://doi.org/10.1016/j.jag.2023.103491" ext-link-type="DOI">10.1016/j.jag.2023.103491</ext-link>, 2023.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>Plastic film residues on cropland: monitoring soil contamination through optical remote sensing</article-title-html>
<abstract-html/>
<ref-html id="bib1.bib1"><label>1</label><mixed-citation>
      
Agriculture Plastic Environment Europe: Statistics – plasticulture in Europe, <a href="https://apeeurope.eu/statistics/" target="_blank"/> (last access: 22 December 2019), 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>2</label><mixed-citation>
      
Bonifazi, G., Francesconi, E., Gasbarrone, R., Palmieri, R., and Serranti, S.: A preliminary study on the utilization of hyperspectral imaging for the
on-soil recognition of plastic waste resulting from agricultural activities,
Land, 12, <a href="https://doi.org/10.3390/land12101934" target="_blank">https://doi.org/10.3390/land12101934</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>3</label><mixed-citation>
      
Breiman, L.: Random forests, Mach. Learn., 45, 5–32, <a href="https://doi.org/10.1023/A:1010933404324" target="_blank">https://doi.org/10.1023/A:1010933404324</a>, 2001.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>4</label><mixed-citation>
      
Castagna, A., Dierssen, H. M., Devriese, L. I., Everaert, G., Knaeps, E.,
and Sterckx, S.: Evaluation of historic and new detection algorithms for
different types of plastics over land and water from hyperspectral data and
imagery, Remote Sens. Environ., 298, <a href="https://doi.org/10.1016/j.rse.2023.113834" target="_blank">https://doi.org/10.1016/j.rse.2023.113834</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>5</label><mixed-citation>
      
Colomina, I. and Molina, P.: Unmanned aerial systems for photogrammetry and
remote sensing: a review, ISPRS J. Photogram. Remote Sens., 92, 79–97,
<a href="https://doi.org/10.1016/j.isprsjprs.2014.02.013" target="_blank">https://doi.org/10.1016/j.isprsjprs.2014.02.013</a>, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>6</label><mixed-citation>
      
Crucil, G.: Application of the radiometric calibration model and image
correction: from raw Rededge images to calibrated reflectanc, GitHub [code],
<a href="https://github.com/Debrisman/Micasense_Rededge_correction/blob/main/FUN_RededgeM_processing.R" target="_blank"/> (last access: 15 February 2022), 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>7</label><mixed-citation>
      
Crucil, G. and Van Oost, K.: Towards mapping of soil crust using multispectral imaging, Sensors, 21, <a href="https://doi.org/10.3390/s21051850" target="_blank">https://doi.org/10.3390/s21051850</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>8</label><mixed-citation>
      
Crucil, G., Castaldi, F., Aldana-Jague, E., van Wesemael, B., Macdonald, A.,
and Van Oost, K.: Assessing the performance of UAS-compatible multispectral
and hyperspectral sensors for soil organic carbon prediction, Sustainability, 11, <a href="https://doi.org/10.3390/su11071889" target="_blank">https://doi.org/10.3390/su11071889</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>9</label><mixed-citation>
      
Despini, F., Costanzini, S., Parmeggiani, D., and la Cecilia, D.: Analysis
of hyperspectral data to define multispectral sensor specifications for
enhancing plastic greenhouse detection, Remote Sens. Appl. Soc. Environ.,
101802, <a href="https://doi.org/10.1016/j.rsase.2025.101802" target="_blank">https://doi.org/10.1016/j.rsase.2025.101802</a>, 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>10</label><mixed-citation>
      
EIP-AGRI Focus Group: Reducing the plastic footprint of agriculture – final
report, European Commission, <a href="https://eu-cap-network.ec.europa.eu/publications/eip-agri-focus-group-plastic-footprint-final-report_en" target="_blank"/>
(last access: 29 December 2021), 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>11</label><mixed-citation>
      
Espí, E., Salmerón, A., Fontecha, A., García, Y., and Real, A. I.: Plastic films for agricultural applications, J. Plast. Film Sheet., 22, 85–102, <a href="https://doi.org/10.1177/8756087906064220" target="_blank">https://doi.org/10.1177/8756087906064220</a>, 2006.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>12</label><mixed-citation>
      
Fabrizi, A., Fiener, P., Van Oost, K., and Wilken, F.: Hyperspectral and
multispectral reflectance of agricultural plastic films, Zenodo [data set],
<a href="https://doi.org/10.5281/zenodo.14336252" target="_blank">https://doi.org/10.5281/zenodo.14336252</a>, 2024.


    </mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>13</label><mixed-citation>
      
FAO: Assessment of agricultural plastics and their sustainability. A call
for action, <a href="https://doi.org/10.4060/cb7856en" target="_blank">https://doi.org/10.4060/cb7856en</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>14</label><mixed-citation>
      
Frias, J. P. G. L. and Nash, R.: Microplastics: finding a consensus on the
definition, Mar. Pollut. Bull., 138, 145–147, <a href="https://doi.org/10.1016/j.marpolbul.2018.11.022" target="_blank">https://doi.org/10.1016/j.marpolbul.2018.11.022</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>15</label><mixed-citation>
      
Garaba, S. P.: Towards advanced mapping of plastic greenhouses from EMIT,
EnMAP and PRISMA hyperspectral missions, Environ. Chall., 22, 101398,
<a href="https://doi.org/10.1016/j.envc.2025.101398" target="_blank">https://doi.org/10.1016/j.envc.2025.101398</a>, 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>16</label><mixed-citation>
      
Garaba, S. P. and Dierssen, H. M.: An airborne remote sensing case study of
synthetic hydrocarbon detection using short wave infrared absorption features identified from marine-harvested macro- and microplastics, Remote Sens. Environ., 205, 224–235, <a href="https://doi.org/10.1016/j.rse.2017.11.023" target="_blank">https://doi.org/10.1016/j.rse.2017.11.023</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>17</label><mixed-citation>
      
Garaba, S. P., Aitken, J., Slat, B., Dierssen, H. M., Lebreton, L., Zielinski, O., and Reisser, J.: Sensing ocean plastics with an airborne hyperspectral shortwave infrared omager, Environ. Sci. Technol., 52,
11699–11707, <a href="https://doi.org/10.1021/acs.est.8b02855" target="_blank">https://doi.org/10.1021/acs.est.8b02855</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>18</label><mixed-citation>
      
Garaba, S. P., Arias, M., Corradi, P., Harmel, T., de Vries, R., and Lebreton, L.: Concentration, anisotropic and apparent colour effects on
optical reflectance properties of virgin and ocean-harvested plastics, J.
Hazard. Mater., 406, 124290, <a href="https://doi.org/10.1016/j.jhazmat.2020.124290" target="_blank">https://doi.org/10.1016/j.jhazmat.2020.124290</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>19</label><mixed-citation>
      
Goddijn-Murphy, L. and Dufaur, J.: Proof of concept for a model of light
reflectance of plastics floating on natural waters, Mar. Pollut. Bull., 135,
1145–1157, <a href="https://doi.org/10.1016/j.marpolbul.2018.08.044" target="_blank">https://doi.org/10.1016/j.marpolbul.2018.08.044</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>20</label><mixed-citation>
      
Goddijn-Murphy, L., Peters, S., van Sebille, E., James, N. A., and Gibb, S.:
Concept for a hyperspectral remote sensing algorithm for floating marine macro plastics, Mar. Pollut. Bull., 126, 255–262,
<a href="https://doi.org/10.1016/j.marpolbul.2017.11.011" target="_blank">https://doi.org/10.1016/j.marpolbul.2017.11.011</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>21</label><mixed-citation>
      
Hively, W. D., Lamb, B. T., Daughtry, C. S. T., Quemada, M., Hunt Jr., R., Serbin, G., and Dennison, P.: Reflectance spectra of agricultural field
conditions supporting remote-sensing evaluation of non-photosynthetic
vegetation cover (ver. 1.1, November 2022), US Geological Survey data release [data set], <a href="https://doi.org/10.5066/P9XK3867" target="_blank">https://doi.org/10.5066/P9XK3867</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>22</label><mixed-citation>
      
Horton, A. A., Walton, A., Spurgeon, D. J., Lahive, E., and Svendsen, C.:
Microplastics in freshwater and terrestrial environments: Evaluating the current understanding to identify the knowledge gaps and future research
priorities, Sci. Total Environ., 586, 127–141, <a href="https://doi.org/10.1016/j.scitotenv.2017.01.190" target="_blank">https://doi.org/10.1016/j.scitotenv.2017.01.190</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>23</label><mixed-citation>
      
Iordache, M.-D., De Keukelaere, L., Moelans, R., Landuyt, L., Moshtaghi, M.,
Corradi, P., and Knaeps, E.: Targeting plastics: machine learning applied to
litter detection in aerial multispectral images, Remote Sens., 14,
<a href="https://doi.org/10.3390/rs14225820" target="_blank">https://doi.org/10.3390/rs14225820</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>24</label><mixed-citation>
      
Jones, H., Black, T. A., Jassal, R. S., Nesic, Z., Johnson, M. S., and Smukler, S.: Characterization of shortwave and longwave properties of several plastic film mulches and their impact on the surface energy balance and soil temperature, Sol. Energy, 214, 457–470, <a href="https://doi.org/10.1016/j.solener.2020.11.058" target="_blank">https://doi.org/10.1016/j.solener.2020.11.058</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>25</label><mixed-citation>
      
Julienne, F., Delorme, N., and Lagarde, F.: From macroplastics to microplastics: role of water in the fragmentation of polyethylene, Chemosphere, 236, 124409, <a href="https://doi.org/10.1016/j.chemosphere.2019.124409" target="_blank">https://doi.org/10.1016/j.chemosphere.2019.124409</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>26</label><mixed-citation>
      
Junta de Andalucía: Mapa de Suelos de Andalucía,
<a href="https://portalrediam.cica.es/VisorRediam/?lyr.add.wms=http://medioambiente/mapwms/REDIAM_Suelos_Andalucia" target="_blank"/> (last access: 9 April 2025), 2005.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>27</label><mixed-citation>
      
Kawecki, D. and Nowack, B.: Polymer-specific modeling of the environmental
emissions of seven commodity plastics as macro- and microplastics, Environ.
Sci. Technol., 53, 9664–9676, <a href="https://doi.org/10.1021/acs.est.9b02900" target="_blank">https://doi.org/10.1021/acs.est.9b02900</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib28"><label>28</label><mixed-citation>
      
Kedzierski, M., Cirederf-Boulant, D., Palazot, M., Yvin, M., and Bruzaud, S.: Continents of plastics: an estimate of the stock of microplastics in agricultural soils, Sci. Total Environ., 880, 163294, <a href="https://doi.org/10.1016/j.scitotenv.2023.163294" target="_blank">https://doi.org/10.1016/j.scitotenv.2023.163294</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib29"><label>29</label><mixed-citation>
      
Knaeps, E., Sterckx, S., Strackx, G., Mijnendonckx, J., Moshtaghi, M., Garaba, S. P., and Meire, D.: Hyperspectral-reflectance dataset of dry, wet and submerged marine litter, Earth Syst. Sci. Data, 13, 713–730,
<a href="https://doi.org/10.5194/essd-13-713-2021" target="_blank">https://doi.org/10.5194/essd-13-713-2021</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib30"><label>30</label><mixed-citation>
      
Kokaly, R. F., Clark, R. N., Swayze, G. A., Livo, K. E., Hoefen, T. M., Pearson, N. C., Wise, R. A., Benzel, W. M., Lowers, H. A., Driscoll, R. L., and Klein, A. J.: USGS Spectral Library Version 7 Data, US Geological Survey data release [data set], <a href="https://doi.org/10.5066/F7RR1WDJ" target="_blank">https://doi.org/10.5066/F7RR1WDJ</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib31"><label>31</label><mixed-citation>
      
Kühn, F., Oppermann, K., and Hörig, B.: Hydrocarbon Index – an
algorithm for hyperspectral detection of hydrocarbons, Int. J. Remote Sens.,
25, 2467–2473, <a href="https://doi.org/10.1080/01431160310001642287" target="_blank">https://doi.org/10.1080/01431160310001642287</a>, 2004.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib32"><label>32</label><mixed-citation>
      
Lamont, W. J.: Plastic mulches for the production of vegetable crops, in: A
guide to the manufacture, performance, and potential of plastics in
agriculture, in: Plastic design library, edited by: Orzolek, M. D., Elsevier,
45–60, <a href="https://doi.org/10.1016/b978-0-08-102170-5.00003-8" target="_blank">https://doi.org/10.1016/b978-0-08-102170-5.00003-8</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib33"><label>33</label><mixed-citation>
      
Landrigan, P. J., Raps, H., Cropper, M., Bald, C., Brunner, M., Canonizado,
E. M., Charles, D., Chiles, T. C., Donohue, M. J., Enck, J., Fenichel, P.,
Fleming, L. E., Ferrier-Pages, C., Fordham, R., Gozt, A., Griffin, C., Hahn, M. E., Haryanto, B., Hixson, R., Ianelli, H., James, B. D., Kumar, P.,
Laborde, A., Law, K. L., Martin, K., Mu, J., Mulders, Y., Mustapha, A., Niu,
J., Pahl, S., Park, Y., Pedrotti, M.-L., Pitt, J. A., Ruchirawat, M., Seewoo, B. J., Spring, M., Stegeman, J. J., Suk, W., Symeonides, C., Takada, H., Thompson, R. C., Vicini, A., Wang, Z., Whitman, E., Wirth, D., Wolff, M., Yousuf, A. K., and Dunlop, S.: The Minderoo-Monaco commission on plastics and human health, Ann. Global Health, 89, 23, <a href="https://doi.org/10.5334/aogh.4056" target="_blank">https://doi.org/10.5334/aogh.4056</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib34"><label>34</label><mixed-citation>
      
Liu, M., Feng, J., Shen, Y., and Zhu, B.: Microplastics effects on soil biota are dependent on their properties: a meta-analysis, Soil Biol. Biochem., 178, 108940, <a href="https://doi.org/10.1016/j.soilbio.2023.108940" target="_blank">https://doi.org/10.1016/j.soilbio.2023.108940</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib35"><label>35</label><mixed-citation>
      
Lofty, J., Muhawenimana, V., Wilson, C. A. M. E., and Ouro, P.: Microplastics removal from a primary settler tank in a wastewater treatment plant and estimations of contamination onto European agricultural land via sewage sludge recycling, Environ. Pollut., 304, 119198, <a href="https://doi.org/10.1016/j.envpol.2022.119198" target="_blank">https://doi.org/10.1016/j.envpol.2022.119198</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib36"><label>36</label><mixed-citation>
      
Martin, C., Parkes, S., Zhang, Q., Zhang, X., McCabe, M. F., and Duarte, C. M.: Use of unmanned aerial vehicles for efficient beach litter monitoring, Mar. Pollut. Bull., 131, 662–673, <a href="https://doi.org/10.1016/j.marpolbul.2018.04.045" target="_blank">https://doi.org/10.1016/j.marpolbul.2018.04.045</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib37"><label>37</label><mixed-citation>
      
Meixner, K., Kubiczek, M., and Fritz, I.: Microplastic in soil–current status in Europe with special focus on method tests with Austrian samples, AIMS Environ. Sci., 7, 174–191, <a href="https://doi.org/10.3934/environsci.2020011" target="_blank">https://doi.org/10.3934/environsci.2020011</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib38"><label>38</label><mixed-citation>
      
MicaSense: Radiometric calibration model for MicaSense sensors,
<a href="https://support.micasense.com/hc/en-us/articles/115000351194-Radiometric-Calibration-Model-for-MicaSense-Sensors" target="_blank">https://support.micasense.com/hc/en-us/articles/115000351194-Radiometric-Calibration-Model-for-MicaSense-
Sensors</a> (last access: 11 April 2025), 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib39"><label>39</label><mixed-citation>
      
Moshtaghi, M., Knaeps, E., Sterckx, S., Garaba, S. P., and Meire, D.: Spectral reflectance of marine macroplastics in the VNIR and SWIR measured
in a controlled environment, Sci. Rep., 11, 5436, <a href="https://doi.org/10.1038/s41598-021-84867-6" target="_blank">https://doi.org/10.1038/s41598-021-84867-6</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib40"><label>40</label><mixed-citation>
      
Muñoz-Rojas, M., Jordán, A., Zavala, L., De la Rosa, D., Abd-Elmabod, S., and Anaya-Romero, M.: Organic carbon stocks in Mediterranean soil types under different land uses (southern Spain), Solid Earth, 3, 375–386, <a href="https://doi.org/10.5194/se-3-375-2012" target="_blank">https://doi.org/10.5194/se-3-375-2012</a>, 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib41"><label>41</label><mixed-citation>
      
Pahlawan, M. F. R., Kim, Y., Aline, U., Zahroh, A., Masithoh, R. E., Kim, M.
S., Baek, I., and Cho, B.-K.: Non-destructive identification of microplastics in soil using spectroscopy and hyperspectral imaging, Trends Anal. Chem., 118216, <a href="https://doi.org/10.1016/j.trac.2025.118216" target="_blank">https://doi.org/10.1016/j.trac.2025.118216</a>, 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib42"><label>42</label><mixed-citation>
      
Piehl, S., Leibner, A., Löder, M. G., Dris, R., Bogner, C., and Laforsch, C.: Identification and quantification of macro-and microplastics on an agricultural farmland, Sci. Rep., 8, 17950, <a href="https://doi.org/10.1038/s41598-018-36172-y" target="_blank">https://doi.org/10.1038/s41598-018-36172-y</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib43"><label>43</label><mixed-citation>
      
Plastics Europe AISBL: The circular economy for plastics – a European
analysis, <a href="https://plasticseurope.org/wp-content/uploads/2024/03/CEreport_executivesummary_2024.pdf" target="_blank"/> (last access: 11 April 2025), 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib44"><label>44</label><mixed-citation>
      
Qiu, F., Zhai, Z., Li, Y., Yang, J., Wang, H., and Zhang, R.: UAV imaging and deep learning based method for predicting residual film in cotton field plough layer, Front. Plant Sci., 13, 1010474, <a href="https://doi.org/10.3389/fpls.2022.1010474" target="_blank">https://doi.org/10.3389/fpls.2022.1010474</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib45"><label>45</label><mixed-citation>
      
Rehm, R., Zeyer, T., Schmidt, A., and Fiener, P.: Soil erosion as transport
pathway of microplastic from agriculture soils to aquatic ecosystems, Sci.
Total Environ., 795, <a href="https://doi.org/10.1016/j.scitotenv.2021.148774" target="_blank">https://doi.org/10.1016/j.scitotenv.2021.148774</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib46"><label>46</label><mixed-citation>
      
Rezaei, M., Abbasi, S., Pourmahmood, H., Oleszczuk, P., Ritsema, C., and Turner, A.: Microplastics in agricultural soils from a semi-arid region and
their transport by wind erosion, Environ. Res., 212,
<a href="https://doi.org/10.1016/j.envres.2022.113213" target="_blank">https://doi.org/10.1016/j.envres.2022.113213</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib47"><label>47</label><mixed-citation>
      
Savitzky, A. and Golay, M. J.: Smoothing and differentiation of data by
simplified least squares procedures, Anal. Chem., 36, 1627–1639,
<a href="https://doi.org/10.1021/ac60214a047" target="_blank">https://doi.org/10.1021/ac60214a047</a>, 1964.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib48"><label>48</label><mixed-citation>
      
Scarascia-Mugnozza, G., Sica, C., and Russo, G.: Plastic materials in European agriculture: actual use and perspectives, J. Agric. Eng., 42, 15–28, <a href="https://doi.org/10.4081/jae.2011.3.1" target="_blank">https://doi.org/10.4081/jae.2011.3.1</a>, 2011.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib49"><label>49</label><mixed-citation>
      
Shan, J., Zhao, J., Liu, L., Zhang, Y., Wang, X., and Wu, F.: A novel way to
rapidly monitor microplastics in soil by hyperspectral imaging technology and chemometrics, Environ. Pollut., 238, 121–129, <a href="https://doi.org/10.1016/j.envpol.2018.03.026" target="_blank">https://doi.org/10.1016/j.envpol.2018.03.026</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib50"><label>50</label><mixed-citation>
      
Shaw, G. A. and Burke, H. K.: Spectral imaging for remote sensing, Lincoln Lab. J., 14, 3–28, 2003.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib51"><label>51</label><mixed-citation>
      
Song, Y. K., Hong, S. H., Jang, M., Han, G. M., Jung, S. W., and Shim, W. J.: Combined effects of UV exposure duration and mechanical abrasion on microplastic fragmentation by polymer type, Environ. Sci. Technol., 51,
4368–4376, <a href="https://doi.org/10.1021/acs.est.6b06155" target="_blank">https://doi.org/10.1021/acs.est.6b06155</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib52"><label>52</label><mixed-citation>
      
Stefano, N. and Pleissner, D.: Quantification and analysis of surface macroplastic contamination on arable areas, J. Soils Sediments, 22, 757–768,
<a href="https://doi.org/10.1007/s11368-021-03114-7" target="_blank">https://doi.org/10.1007/s11368-021-03114-7</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib53"><label>53</label><mixed-citation>
      
Steinmetz, Z., Löffler, P., Eichhöfer, S., David, J., Muñoz, K.,
and Schaumann, G. E.: Are agricultural plastic covers a source of plastic
debris in soil? A first screening study, Soil, 8, 31–47, <a href="https://doi.org/10.5194/soil-8-31-2022" target="_blank">https://doi.org/10.5194/soil-8-31-2022</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib54"><label>54</label><mixed-citation>
      
Stevens, A. and Ramirez-Lopez, L.: An introduction to the prospectr package,
<a href="https://cran.r-project.org/web/packages/prospectr/vignettes/prospectr.html" target="_blank"/>
(last access: 10 April 2025), 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib55"><label>55</label><mixed-citation>
      
SunCalc.org.:
<a href="https://www.suncalc.org/#/50.6693,4.62,17/2022.07.24/10:00/1/3" target="_blank"/> (last
access: 11 April 2025), 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib56"><label>56</label><mixed-citation>
      
Thompson, R. C., Olsen, Y., Mitchell, R. P., Davis, A., Rowland, S. J., John, A. W., McGonigle, D., and Russell, A. E.: Lost at sea: where is all the plastic?, Science, 304, 838, <a href="https://doi.org/10.1126/science.1094559" target="_blank">https://doi.org/10.1126/science.1094559</a>, 2004.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib57"><label>57</label><mixed-citation>
      
Tucker, C. J.: Red and photographic infrared linear combinations for monitoring vegetation, Remote Sens. Environ., 8, 127–150,
<a href="https://doi.org/10.1016/0034-4257(79)90013-0" target="_blank">https://doi.org/10.1016/0034-4257(79)90013-0</a>, 1979.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib58"><label>58</label><mixed-citation>
      
Veenstra, T. S. and Churnside, J. H.: Airborne sensors for detecting large
marine debris at sea, Mar. Pollut. Bull., 65, 63–68,
<a href="https://doi.org/10.1016/j.marpolbul.2010.11.018" target="_blank">https://doi.org/10.1016/j.marpolbul.2010.11.018</a>, 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib59"><label>59</label><mixed-citation>
      
Wang, P.-Y., Zhao, Z.-Y., Xiong, X.-B., Wang, N., Zhou, R., Zhang, Z.-M., Ding, F., Hao, M., Wang, S., Ma, Y., Uzamurera, A. G., Xiao, K.-W., Khan, A., Tao, X.-P., Wang, W.-Y., and Tao, H.-Y.: Microplastics affect soil bacterial community assembly more by their shapes rather than the concentrations, Water Res., 245, 120581, <a href="https://doi.org/10.1016/j.watres.2023.120581" target="_blank">https://doi.org/10.1016/j.watres.2023.120581</a>, 2023a.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib60"><label>60</label><mixed-citation>
      
Wang, S., Fan, T., Cheng, W., Wang, L., Zhao, G., Li, S., Dang, Y., and Zhang, J.: Occurrence of macroplastic debris in the long-term plastic film-mulched agricultural soil: a case study of Northwest China, Sci. Total
Environ., 831, 154881, <a href="https://doi.org/10.1016/j.scitotenv.2022.154881" target="_blank">https://doi.org/10.1016/j.scitotenv.2022.154881</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib61"><label>61</label><mixed-citation>
      
Wang, Y., Xiang, L., Amelung, W., Elsner, M., Gan, J., Kueppers, S., Christian, L., Jiang, X., Adu-Gyamfi, J., Heng, L., Sik Ok, Y., Ivleva, N.
P., Luo, Y., Barceló, D., Schäffer, A., and Wang, F.: Micro-and
nanoplastics in soil ecosystems: Analytical methods, fate, and effects,
Trends Anal. Chem., 169, 117309, <a href="https://doi.org/10.1016/j.trac.2023.117309" target="_blank">https://doi.org/10.1016/j.trac.2023.117309</a>, 2023b.


    </mixed-citation></ref-html>
<ref-html id="bib1.bib62"><label>62</label><mixed-citation>
      
Yang, J., Zhai, Z., Li, Y., Duan, H., Cai, F., Lv, J., and Zhang, R.: Design
and research of residual film pollution monitoring system based on UAV, Comput. Electron. Agr., 217, <a href="https://doi.org/10.1016/j.compag.2023.108608" target="_blank">https://doi.org/10.1016/j.compag.2023.108608</a>, 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib63"><label>63</label><mixed-citation>
      
Yang, Y., Li, Z., Yan, C., Chadwick, D., Jones, D. L., Liu, E., Liu, Q., Bai, R., and He, W.: Kinetics of microplastic generation from different types of mulch films in agricultural soil, Sci. Total Environ., 814, 152572,
<a href="https://doi.org/10.1016/j.scitotenv.2021.152572" target="_blank">https://doi.org/10.1016/j.scitotenv.2021.152572</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib64"><label>64</label><mixed-citation>
      
Zhai, Z., Chen, X., Zhang, R., Qiu, F., Meng, Q., Yang, J., and Wang, H.:
Evaluation of residual plastic film pollution in pre-sowing cotton field
using UAV imaging and semantic segmentation, Front. Plant Sci., 13, 991191,
<a href="https://doi.org/10.3389/fpls.2022.991191" target="_blank">https://doi.org/10.3389/fpls.2022.991191</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib65"><label>65</label><mixed-citation>
      
Zhang, H., Shi, P., Crucil, G., Wesemael, B., Limbourg, Q., and Van Oost, K.: Evaluating the capability of a UAV-borne spectrometer for soil organic carbon mapping in bare croplands, Land Degrad. Dev., 32, 4375–4389,
<a href="https://doi.org/10.1002/ldr.4043" target="_blank">https://doi.org/10.1002/ldr.4043</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib66"><label>66</label><mixed-citation>
      
Zhang, J., Ren, S., Xu, W., Liang, C., Li, J., Zhang, H., Li, Y., Liu, X.,
Jones, D. L., Chadwick, D. R., Zhang, F., and Wang, K.: Effects of plastic
residues and microplastics on soil ecosystems: a global meta-analysis, J.
Hazard. Mater., 435, <a href="https://doi.org/10.1016/j.jhazmat.2022.129065" target="_blank">https://doi.org/10.1016/j.jhazmat.2022.129065</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib67"><label>67</label><mixed-citation>
      
Zhou, S., Kuester, T., Bochow, M., Bohn, N., Brell, M., and Kaufmann, H.: A
knowledge-based, validated classifier for the identification of aliphatic
and aromatic plastics by WorldView-3 satellite data, Remote Sens. Environ.,
264, <a href="https://doi.org/10.1016/j.rse.2021.112598" target="_blank">https://doi.org/10.1016/j.rse.2021.112598</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib68"><label>68</label><mixed-citation>
      
Zhou, S., Mou, L., Hua, Y., Zhang, L., Kaufmann, H., and Zhu, X. X.: Can we
use deep learning models to identify the functionality of plastics from
space?, Int. J. Appl. Earth Obs. Geoinf., 123, <a href="https://doi.org/10.1016/j.jag.2023.103491" target="_blank">https://doi.org/10.1016/j.jag.2023.103491</a>, 2023.

    </mixed-citation></ref-html>--></article>
