<?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">AMT</journal-id><journal-title-group>
    <journal-title>Atmospheric Measurement Techniques</journal-title>
    <abbrev-journal-title abbrev-type="publisher">AMT</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Atmos. Meas. Tech.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">1867-8548</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/amt-19-4141-2026</article-id><title-group><article-title>Assessment of the RFI environment in key passive microwave bands for Earth observation</article-title><alt-title>Systematic assessment of the RFI environment in key passive microwave bands</alt-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Onrubia</surname><given-names>Raul</given-names></name>
          <email>onrubia@zenithalblue.com</email>
        <ext-link>https://orcid.org/0000-0001-6636-1587</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Oliva</surname><given-names>Roger</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="yes" rid="aff2">
          <name><surname>Duncan</surname><given-names>David</given-names></name>
          <email>david.duncan@ecmwf.int</email>
        <ext-link>https://orcid.org/0000-0002-1955-4391</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff2">
          <name><surname>Bormann</surname><given-names>Niels</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff3">
          <name><surname>Barbosa</surname><given-names>Jose</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff3">
          <name><surname>Nestoras</surname><given-names>Ioannis</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff3">
          <name><surname>Jordão</surname><given-names>Adriano</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff4">
          <name><surname>Jorge</surname><given-names>Flavio</given-names></name>
          
        <ext-link>https://orcid.org/0000-0002-1656-7824</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff4">
          <name><surname>Challot</surname><given-names>Juliette</given-names></name>
          
        <ext-link>https://orcid.org/0009-0008-4350-4272</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff4">
          <name><surname>Soldo</surname><given-names>Yan</given-names></name>
          
        </contrib>
        <aff id="aff1"><label>1</label><institution>ZenithalBlue Technologies S.L.U., Barcelona, Spain</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>European Centre for Medium-Range Weather Forecasts (ECMWF), Reading, United Kingdom</institution>
        </aff>
        <aff id="aff3"><label>3</label><institution>Research and Development for Aerospace S.L.U., Zürich, Switzerland</institution>
        </aff>
        <aff id="aff4"><label>4</label><institution>European Space Agency, ESTEC, 2201 Noordwijk, the Netherlands</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Raul Onrubia (onrubia@zenithalblue.com) and David Duncan (david.duncan@ecmwf.int)</corresp></author-notes><pub-date><day>25</day><month>June</month><year>2026</year></pub-date>
      
      <volume>19</volume>
      <issue>12</issue>
      <fpage>4141</fpage><lpage>4163</lpage>
      <history>
        <date date-type="received"><day>30</day><month>September</month><year>2025</year></date>
           <date date-type="rev-request"><day>4</day><month>December</month><year>2025</year></date>
           <date date-type="rev-recd"><day>28</day><month>April</month><year>2026</year></date>
           <date date-type="accepted"><day>25</day><month>May</month><year>2026</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2026 Raul Onrubia 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://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026.html">This article is available from https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026.html</self-uri><self-uri xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026.pdf">The full text article is available as a PDF file from https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026.pdf</self-uri>
      <abstract><title>Abstract</title>

      <p id="d2e187">Radio Frequency Interference (RFI) is spreading worldwide, affecting numerous Earth Observation (EO) instruments. Among these, microwave radiometers play an essential role, providing critical measurements for climate monitoring, weather forecasting, and numerous other applications. In order to plan for future satellite missions and fully exploit currently available measurements, it is crucial to study the contamination levels at bands where radiometers operate. This work presents the Earth Observation RFI Scanner (EORFIScan), an RFI detection system for EO products that combines multiple RFI detection techniques in order to reduce missed detections. This software has been used to assess several passive microwave bands from 6 up to 200 GHz, including both exclusive and shared bands. Analysis and validation of this method is presented for the year 2022. The resulting RFI probability maps show significant contamination in the bands up to and including 18.7 GHz; thus for the assessment, we mainly focus on bands from 6 to 18.7 GHz, where the most widespread contamination is found. A few brightness temperatures in the range of 350–400 K have been observed at 23.8 GHz and one at 36.5 GHz, which suggest the presence of man-made emissions. At higher frequencies, RFI contamination is not clearly visible in the analysed data. Comparisons with simulated radiances from a numerical weather prediction model are presented as a way to evaluate the RFI detection, finding that flagged observations are typically warmer than model simulations, as would be expected for RFI. It is clear from the results presented that RFI is already a concern for users at lower frequency passive microwave bands, and it is recommended that real-time monitoring systems are developed to keep an eye on the evolving threat of RFI in EO bands.</p>
  </abstract>
    
<funding-group>
<award-group id="gs1">
<funding-source>European Space Agency</funding-source>
<award-id>4000141035/23/NL/SD</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="d2e199">Radio spectrum monitoring and enforcement actions are key to deliver radio protection and ensure that the different radio services operate in compliance with the in-force radio regulations and to secure an efficient and effective spectrum use, free from Radio Frequency Interference (RFI) <xref ref-type="bibr" rid="bib1.bibx37" id="paren.1"/>. Notwithstanding these essential actions, RFI is an inevitable consequence of the operation of radio systems, and it is a growing threat that has affected numerous Earth Observation (EO) satellite missions in recent years, especially passive microwave sensors <xref ref-type="bibr" rid="bib1.bibx10 bib1.bibx32 bib1.bibx38" id="paren.2"/>. These interferences vary widely in power level, radiation pattern, polarization, bandwidth, time duration, duty cycle, etc., and as a consequence there is no single algorithm that can optimally detect all levels and types of radiometric contamination <xref ref-type="bibr" rid="bib1.bibx40" id="paren.3"/>.</p>
      <p id="d2e211">The European Space Agency (ESA) and other space agencies rely on microwave radiometers for a wide variety of Earth Observation (EO) missions, operating in frequency bands from 1 to 700 GHz. RFI can negatively impact these radiometers in three different ways: High power RFIs can potentially damage the instrument, but have low occurrence probability; medium power RFIs cause data loss and gaps; and low power RFIs, mostly undetected, can distort scientific data, potentially compromising geophysical retrievals and weather forecasts.</p>
      <p id="d2e214">It is thus of utmost importance to detect such RFI occurrences, not only to address their impact on the users' side, but also to support respective radio protection actions. However, instances of interference can vary widely not just in magnitude as mentioned above, but also in whether it is possible or appropriate to assign responsibility. Albeit the role of the International Telecommunications Union (ITU) in facilitating international coordination of radio protection actions, national administrations, by means of their National Regulatory Authorities (NRA), are responsible for overseeing the spectrum use and enforce the law in the territory under their jurisdiction. The bands typically used in EO missions fall into EESS (Earth Exploration-Satellite Service), including passive allocations with either exclusive access or shared access with active services. Service disruptions should always be reported to NRA who are responsible for conducting the required investigations and actions on the field to secure integrity of the spectrum so it can be effectively and efficiently used. In the case of exclusive access bands, clear violations, from the law enforcement standpoint, are easier to act upon. In case of shared access bands, the regulatory framework is still to be enforced. If the regulatory framework is enforced and service disruptions still occur, those findings are relevant in closing the loop for future enhanced spectrum policy, planning and assignment conditions. Comprehensive details about radio protection policies and practices can be found in <xref ref-type="bibr" rid="bib1.bibx37" id="paren.4"/>. There are numerous instances of the process to locate and shut down illegal transmissions ongoing, but new RFI sources appear constantly and the procedure to detect, attribute, and shut down illegal operators takes time.</p>
      <p id="d2e220">In addition to clear violations of ITU regulations, many instances of apparent RFI come from legal operations in nearby bands or simply other users in a non-exclusive band. Furthermore, the spectral response of a radiometer can vary significantly from its specified bandwidth, and in this case apparent violations may be caused by radiometer sensitivity outside of the exclusive-use band due to being manufactured outside specified tolerances. Many space agencies now publicly provide measured spectral response functions (SRFs) for MW radiometers <xref ref-type="bibr" rid="bib1.bibx42" id="paren.5"/> and it is a recognized best practice <xref ref-type="bibr" rid="bib1.bibx14" id="paren.6"/>, but it is rare for heritage sensors and indeed most of the instruments studied here do not have measured SRFs publicly available. In other words, the cause of apparent RFI may be due to instrument design, illegal operators, or be a simple consequence of sharing a crowded part of the spectrum. No matter their source, satellite operators and downstream users need to address the inevitable presence of unwanted emissions.</p>
      <p id="d2e230">Regardless of the source of RFI, and although the RF environment is constantly changing, it remains critically important for the planning of future EO missions to know the current status of contamination at bands of interest. For example, the L-band has been widely studied <xref ref-type="bibr" rid="bib1.bibx32 bib1.bibx3 bib1.bibx36" id="paren.7"/>, and this has led to tangible changes in mission design and requirements for subsequent L-band missions. Nevertheless, the ESA RFI Monitoring &amp; Information Tool (ERMIT) showed that the decreasing trend in worldwide contamination suddenly changed in 2024 when the Russo-Ukrainian war and the Myanmar civil war led to massive contamination in all observations over Europe and most of the South and South East Asia <xref ref-type="bibr" rid="bib1.bibx15" id="paren.8"/>. C-band contamination was first discovered in AMSR-E measurements <xref ref-type="bibr" rid="bib1.bibx28" id="paren.9"/>. The presence of RFI in this band led to inclusion of an extra C-band channel on the follow-on AMSR2 mission for the purpose of RFI mitigation <xref ref-type="bibr" rid="bib1.bibx22" id="paren.10"/>. RFI surveys conducted with different methodologies have only been conducted up to and including the 23.8 GHz band <xref ref-type="bibr" rid="bib1.bibx11 bib1.bibx47 bib1.bibx6 bib1.bibx17" id="paren.11"/>, showing contamination in the 18.7 GHz band but not at higher frequencies. <xref ref-type="bibr" rid="bib1.bibx39" id="text.12"/> reports RFI up to the 36.5 GHz band, albeit not widespread. Despite the best efforts of space agencies to predict the potential contamination issues of future EO missions, the RFI environment could change rapidly as services require more and more bandwidth, technology becomes more accessible and thus widespread, active services are aiming at higher frequencies, and the spectrum becomes more crowded with different users. Spectrum use, in the past, present and future, by Earth observation missions has been reviewed in <xref ref-type="bibr" rid="bib1.bibx24" id="text.13"/>.</p>
      <p id="d2e255">The objective of this study is to survey RFI contamination in several key passive microwave bands from 6 up to 200 GHz using the same methodology, independent on the sensor characteristics (conical scan, across-track scan, …), in all bands commonly used by microwave radiometers using software developed under an ESA contract: the Earth Observation RFI Scanner (EORFIScan). This system is designed to be placed between the data generated from EO missions and the end users, such as the Numerical Weather Prediction (NWP) centres, with the goal of detecting and filtering RFI in order to avoid the assimilation of RFI-contaminated data. These new flags are appended to the original EO products. For this analysis, a comprehensive library of RFI detection algorithms is used, supported by external information and RFI databases. As the system analyses the incoming EO products, the results are stored in internal RFI databases that are used to improve future detection of RFI signals in the products. These databases show a picture of the RFI environment in each frequency band analysed. L-band is not analysed since EORFIScan is an upgrade of the Ground Detection RFI System (GRDS), which already surveyed this frequency band <xref ref-type="bibr" rid="bib1.bibx33" id="paren.14"/>, and frequency bands beyond 200 GHz the number of sensors working at these bands is reduced and RFIs are less likely to happen.</p>
      <p id="d2e261">It is worth emphasizing that proper validation of RFI detection methods is difficult because no source of truth is available. Traditional metrics like a false alarm ratio are then hard to define, in contrast to a parameter like precipitation for which rain gauges offer a local measure of truth. In this study, RFI detection is assessed using simulations from an NWP model, following the idea that NWP model simulations provide information on the expected natural variability, which can be used as reference for RFI detection. This allows some validation of the methods, both instantaneously and in aggregate, but does not lead to specific conclusions on aspects like detection limits or false alarms. Furthermore, even an excellent RFI detection method will miss contamination that is within instrument or geophysical noise levels. This type of low-level RFI would be effectively undetectable on an instantaneous basis. Therefore, the lack of detection does not preclude the existence of RFI in a given band. This is especially true at higher frequencies where natural variability is larger due to more significant atmospheric signals.</p>
      <p id="d2e264">This paper is structured as follows. Section <xref ref-type="sec" rid="Ch1.S2"/> describes the system architecture, the RFI detection techniques utilized, and the threshold computation procedure, Sect. <xref ref-type="sec" rid="Ch1.S3"/> shows the RFI maps of the passive microwave bands from 6 up to 200 GHz, which are discussed in Sect. <xref ref-type="sec" rid="Ch1.S4"/>. Section <xref ref-type="sec" rid="Ch1.S5"/> describes the results validation, and Sect. <xref ref-type="sec" rid="Ch1.S6"/> draws some conclusions.</p>
</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Materials and Methods</title>
      <p id="d2e285">This section presents a brief overview of the system architecture, describes which instruments were analysed during which time period, describes the RFI detection techniques applied in the analysis, and explains the procedure to compute the thresholds of each detection technique.</p>
<sec id="Ch1.S2.SS1">
  <label>2.1</label><title>Instruments Used and Temporal Scope of the Analysis</title>
      <p id="d2e295">In order to minimize the number of different instruments used to cover as many bands as possible, AMSR2 <xref ref-type="bibr" rid="bib1.bibx22" id="paren.15"/>, AMSU-A <xref ref-type="bibr" rid="bib1.bibx30" id="paren.16"/>, MWHS-2 <xref ref-type="bibr" rid="bib1.bibx49" id="paren.17"/>) and AMR-C <xref ref-type="bibr" rid="bib1.bibx29" id="paren.18"/> were chosen. Table <xref ref-type="table" rid="T1"/> shows the passive microwave bands under study and the bands and polarizations covered by each instrument. The numbers inside parentheses state which channels or bands each instrument has in each polarization in that band. AMSR2 provides observations at H (horizontal) and V (vertical) polarizations for all frequency bands, whereas the other instruments only measure one polarization. Cross-track scanners such as AMSU-A and MWHS-2 have polarization that shifts as a function of scan angle, with the polarization typically referred to as quasi-vertical or -horizontal (QV or QH) depending on the polarization at nadir. AMSR2, AMSU-A and MWHS-2 provide 2-D <inline-formula><mml:math id="M1" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> observations per frequency band and polarization. AMR-C is, instead, a single-beam radiometer; therefore, it provides 1-D observations. Using these 4 instruments, all bands above L-band up to 200 GHz were covered. Note that the 18.7 GHz band is observed by two sensors (AMR-C and AMSR2), the 23.8 GHz is observed by 3 (AMR-C, AMSR2 and AMSU-A), and 89 GHz is also observed by 3 (AMSR2, AMSU-A and MWHS-2). This allows cross-validation between instruments. AMSR2 is the sole payload on GCOM-W, and is a one-off instrument. In contrast, the other instruments studied have flown or will fly on multiple satellite platforms. The analysis here uses AMSU-A from Metop-C, MWHS-2 from FY-3E, and AMR-C from Sentinel-6A.</p>

<table-wrap id="T1" specific-use="star"><label>Table 1</label><caption><p id="d2e327">Passive microwave bands covered by this study, and instruments observing them. Parentheses state the channel number(s) that the instrument has in that band. Correspondence between channel numbers and band characteristics can be found in <xref ref-type="bibr" rid="bib1.bibx45" id="text.19"/>.</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="left"/>
     <oasis:colspec colnum="3" colname="col3" align="left"/>
     <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 rowsep="1">
         <oasis:entry colname="col1">Frequency Band</oasis:entry>
         <oasis:entry colname="col2">Band Status</oasis:entry>
         <oasis:entry colname="col3">AMSR2</oasis:entry>
         <oasis:entry colname="col4">AMSU-A</oasis:entry>
         <oasis:entry colname="col5">MWHS-2</oasis:entry>
         <oasis:entry colname="col6">AMR-C</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">6.4–7.7 GHz</oasis:entry>
         <oasis:entry colname="col2">Shared Services</oasis:entry>
         <oasis:entry colname="col3">H+V (6.9, 7.3)</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">10.6–10.7 GHz</oasis:entry>
         <oasis:entry colname="col2">Shared Services</oasis:entry>
         <oasis:entry colname="col3">H+V (10.65)</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">18.6–18.8 GHz</oasis:entry>
         <oasis:entry colname="col2">Shared Services</oasis:entry>
         <oasis:entry colname="col3">H+V (18.7)</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6">H<sup>∗</sup> (Ch. 1)</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">23.6–24 GHz</oasis:entry>
         <oasis:entry colname="col2">Exclusive</oasis:entry>
         <oasis:entry colname="col3">H+V (23.8)</oasis:entry>
         <oasis:entry colname="col4">QV (Ch.1)</oasis:entry>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6">H<sup>∗</sup> (Ch. 2)</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">31.3–31.8 GHz</oasis:entry>
         <oasis:entry colname="col2">Mostly Exclusive in Region 2</oasis:entry>
         <oasis:entry colname="col3"/>
         <oasis:entry colname="col4">QV (Ch. 2)</oasis:entry>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">33.5–34.5 GHz</oasis:entry>
         <oasis:entry colname="col2">Shared Services</oasis:entry>
         <oasis:entry colname="col3"/>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6">H<sup>∗</sup> (Ch. 3)</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">36–37 GHz</oasis:entry>
         <oasis:entry colname="col2">Shared Services</oasis:entry>
         <oasis:entry colname="col3">H+V (36.5)</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">50.2–50.4 GHz</oasis:entry>
         <oasis:entry colname="col2">Exclusive</oasis:entry>
         <oasis:entry colname="col3"/>
         <oasis:entry colname="col4">QV (Ch. 3)</oasis:entry>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">52.6–59.3 GHz</oasis:entry>
         <oasis:entry colname="col2">Exclusive up to 54.25 GHz</oasis:entry>
         <oasis:entry colname="col3"/>
         <oasis:entry colname="col4">QV (Ch. 4, 7)</oasis:entry>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1"/>
         <oasis:entry colname="col2"/>
         <oasis:entry colname="col3"/>
         <oasis:entry colname="col4">QH (Ch. 5, 6, 8–14)</oasis:entry>
         <oasis:entry colname="col5"/>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">86–92 GHz</oasis:entry>
         <oasis:entry colname="col2">Exclusive</oasis:entry>
         <oasis:entry colname="col3">H+V (89)</oasis:entry>
         <oasis:entry colname="col4">QV (Ch. 15)</oasis:entry>
         <oasis:entry colname="col5">QH (Ch. 1)</oasis:entry>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">114.25–122.25 GHz</oasis:entry>
         <oasis:entry colname="col2">Exclusive up to 116 GHz</oasis:entry>
         <oasis:entry colname="col3"/>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">QV (Ch. 2–9)</oasis:entry>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">164–167 GHz</oasis:entry>
         <oasis:entry colname="col2">Shared Services</oasis:entry>
         <oasis:entry colname="col3"/>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">QH (Ch. 10)</oasis:entry>
         <oasis:entry colname="col6"/>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">183 GHz</oasis:entry>
         <oasis:entry colname="col2">Exclusive</oasis:entry>
         <oasis:entry colname="col3"/>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">QV (Ch. 11–15)</oasis:entry>
         <oasis:entry colname="col6"/>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table><table-wrap-foot><p id="d2e333"><sup>∗</sup> AMR-C observes at nadir using H polarization. V polarization is present in the instrument as back-up, since there is no difference between H and V observations when observing at nadir.</p></table-wrap-foot></table-wrap>

      <p id="d2e688">If the product contains the water fraction of each observation, pixels are labeled as “land” (water fraction lower than 5 %), “sea” (water fraction higher than 95 %), and “coast” (water fraction between 5 % and 95 %, and pixels 50 km up to away from the coastline according to CIA World Data Bank II; <xref ref-type="bibr" rid="bib1.bibx8" id="altparen.20"/>). If the product only distinguishes pixels as land or sea, coast pixels are deduced as pixels up to 50 km away from the same coastlines database. Labels “sea-ice” and “storm/clouds” are added afterwards. In the case of AMSR2, sea-ice is detected using the ARTIST sea-ice algorithm <xref ref-type="bibr" rid="bib1.bibx43" id="paren.21"/>, which to the authors knowledge performs well in the Arctic marginal zone <xref ref-type="bibr" rid="bib1.bibx26" id="paren.22"/>, and storms are obtained using the procedure from <xref ref-type="bibr" rid="bib1.bibx9" id="text.23"/>. For AMSU-A, MWHS-2 and AMR-C, the <inline-formula><mml:math id="M6" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula>-means algorithm is used to separate “sea” into “sea-only”, “sea-ice”, and “storm/clouds”.</p>
      <p id="d2e711">Table <xref ref-type="table" rid="T2"/> shows the time span covered by the analysis, which tried to cover all 2022. However, MWHS-2 only provides observations for the 166 GHz band from October 2022 onward (in the case of satellite FY-3E, the first of the series FY that covers this band). In order to cover a similar period of time to be able to observe any seasonality, the time span analyzed from FY-3E covers 2023 instead of 2022.</p>

<table-wrap id="T2"><label>Table 2</label><caption><p id="d2e719">Time span for each satellite. The specified months are included in the analysis.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="3">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="left"/>
     <oasis:colspec colnum="3" colname="col3" align="left"/>
     <oasis:thead>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">Instrument</oasis:entry>
         <oasis:entry colname="col2">T. Start</oasis:entry>
         <oasis:entry colname="col3">T. Stop</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">AMSR2</oasis:entry>
         <oasis:entry colname="col2">January 2022</oasis:entry>
         <oasis:entry colname="col3">December 2022</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">AMSU-A</oasis:entry>
         <oasis:entry colname="col2">January 2022</oasis:entry>
         <oasis:entry colname="col3">December 2022</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">MWHS-2</oasis:entry>
         <oasis:entry colname="col2">October 2022</oasis:entry>
         <oasis:entry colname="col3">November 2023</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">AMR-C</oasis:entry>
         <oasis:entry colname="col2">January 2022</oasis:entry>
         <oasis:entry colname="col3">December 2022</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap>

</sec>
<sec id="Ch1.S2.SS2">
  <label>2.2</label><title>RFI Detection techniques</title>
      <p id="d2e807">Multiple RFI detection techniques are applied to the observations with three progressive threshold levels. Different thresholds will therefore correspond to how strict the detection is, and how confident we are that the detections picked up actual RFI (referred later as “confidence level”). The flags for each frequency band are combined with a logical OR, so they keep the highest confidence level. By doing so, observations have 4 possible flags that can be encoded in 2 bits: no RFI flag, low confidence flag, medium confidence flag, and high confidence flag. The RFI techniques used in the analysis are the following: <list list-type="bullet"><list-item>
      <p id="d2e812"><italic>Intensity algorithm</italic>: An observation with brightness temperature <inline-formula><mml:math id="M7" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is flagged if <inline-formula><mml:math id="M8" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub><mml:mo>&gt;</mml:mo><mml:mi mathvariant="italic">α</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">λ</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, where <inline-formula><mml:math id="M9" display="inline"><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">λ</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> represents the threshold as a function of the latitude <inline-formula><mml:math id="M10" display="inline"><mml:mi mathvariant="italic">λ</mml:mi></mml:math></inline-formula>. We apply this detection algorithm to all assessed instruments.</p></list-item><list-item>
      <p id="d2e871"><italic>Polarization ratio</italic>: An observation with brightness temperatures <inline-formula><mml:math id="M11" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">H</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M12" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">V</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> for H and V polarizations, respectively, is flagged if the polarization ratio <inline-formula><mml:math id="M13" display="inline"><mml:mi>P</mml:mi></mml:math></inline-formula> exceeds a threshold <inline-formula><mml:math id="M14" display="inline"><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">λ</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> that depends on latitude <inline-formula><mml:math id="M15" display="inline"><mml:mi mathvariant="italic">λ</mml:mi></mml:math></inline-formula>. We apply this detection algorithm only to AMSR2. The Polarization ratio is defined by <xref ref-type="bibr" rid="bib1.bibx7" id="text.24"/> as:<disp-formula id="Ch1.E1" content-type="numbered"><label>1</label><mml:math id="M16" display="block"><mml:mrow><mml:mi>P</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">V</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">H</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">V</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">H</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:math></disp-formula></p></list-item><list-item>
      <p id="d2e968"><italic>Spatial Variability</italic>: An observation is flagged if the modulus of the components of the spatial variability <xref ref-type="bibr" rid="bib1.bibx21" id="paren.25"/> is greater than a threshold <inline-formula><mml:math id="M17" display="inline"><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">λ</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> that depends on the latitude. For 2-D imaging, both components of the spatial variability are computed by doing the convolution of the data in the instrument reference frame (field of view number and scan number) with the matrices in Eqs. (<xref ref-type="disp-formula" rid="Ch1.E2"/>) and (<xref ref-type="disp-formula" rid="Ch1.E3"/>). We apply this detection algorithm to all assessed instruments.<disp-formula id="Ch1.E2" content-type="numbered"><label>2</label><mml:math id="M18" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="bold">H</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">SV</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfenced open="(" close=")"><mml:mtable class="array" columnalign="center center center"><mml:mtr><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">1</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula></p>
      <p id="d2e1045"><disp-formula id="Ch1.E3" content-type="numbered"><label>3</label><mml:math id="M19" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="bold">H</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">SV</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfenced close=")" open="("><mml:mtable class="array" columnalign="center center center"><mml:mtr><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">1</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>In the case of AMR-C, that is a 1-D instrument, the convolution is conducted with a vector composed of ten <inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>, a <inline-formula><mml:math id="M21" display="inline"><mml:mn mathvariant="normal">0</mml:mn></mml:math></inline-formula>, and ten <inline-formula><mml:math id="M22" display="inline"><mml:mn mathvariant="normal">1</mml:mn></mml:math></inline-formula>: <inline-formula><mml:math id="M23" display="inline"><mml:mrow><mml:mo>〈</mml:mo><mml:mtable class="array" columnalign="center center center"><mml:mtr><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi mathvariant="normal">…</mml:mi><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mn mathvariant="normal">0</mml:mn></mml:mtd><mml:mtd><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mi mathvariant="normal">…</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>〉</mml:mo></mml:mrow></mml:math></inline-formula></p></list-item><list-item>
      <p id="d2e1159"><italic>Image Enhancement</italic>: An observation is flagged if the modulus of the result of applying a high pass filter to the image <xref ref-type="bibr" rid="bib1.bibx21" id="paren.26"/> is greater than a threshold <inline-formula><mml:math id="M24" display="inline"><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">λ</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> that depends on the latitude <inline-formula><mml:math id="M25" display="inline"><mml:mi mathvariant="italic">λ</mml:mi></mml:math></inline-formula>. To apply the high pass filter in 2-D images, the data in the instrument reference frame (field of view number and scan number) is convoluted with the matrix:<disp-formula id="Ch1.E4" content-type="numbered"><label>4</label><mml:math id="M26" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="bold">H</mml:mi><mml:mi mathvariant="normal">HPF</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfenced close=")" open="("><mml:mtable class="array" columnalign="center center center"><mml:mtr><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1.5</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1.5</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mn mathvariant="normal">8</mml:mn></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1.5</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1.5</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>No useful equivalent has been found for 1-D sensors. We apply this detection algorithm to all assessed instruments except AMR-C.</p></list-item><list-item>
      <p id="d2e1258"><italic>Generalized RFI Index</italic>: An observation is flagged if the generalized RFI Index <inline-formula><mml:math id="M27" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:math></inline-formula> at a given frequency channel <inline-formula><mml:math id="M28" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula> is positive and larger than a given threshold <inline-formula><mml:math id="M29" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula>. The generalized RFI Index <xref ref-type="bibr" rid="bib1.bibx11" id="paren.27"/> is defined as the difference between observation <inline-formula><mml:math id="M30" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:math></inline-formula> and the expected value <inline-formula><mml:math id="M31" display="inline"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="true" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:math></inline-formula>):<disp-formula specific-use="gather" content-type="numbered"><mml:math id="M32" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E5"><mml:mtd><mml:mtext>5</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo><mml:mo>-</mml:mo><mml:mover accent="true"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow><mml:mo mathvariant="normal" stretchy="true">^</mml:mo></mml:mover><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E6"><mml:mtd><mml:mtext>6</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:mover accent="true"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow><mml:mo mathvariant="normal" stretchy="true">^</mml:mo></mml:mover><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>+</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>∈</mml:mo><mml:mfenced close="}" open="{"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>j</mml:mi><mml:mo>)</mml:mo><mml:mo>≠</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>i</mml:mi><mml:mo>)</mml:mo><mml:mo>)</mml:mo><mml:mo>)</mml:mo></mml:mrow></mml:mfenced></mml:mrow></mml:munder><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub><mml:mo>[</mml:mo><mml:mi>j</mml:mi><mml:mo>]</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:msubsup><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>[</mml:mo><mml:mi>j</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>The expected brightness temperature <inline-formula><mml:math id="M33" display="inline"><mml:mover accent="true"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow><mml:mo mathvariant="normal" stretchy="true">^</mml:mo></mml:mover></mml:math></inline-formula> is estimated using the brightness temperatures and their squared values of the other frequency channels <inline-formula><mml:math id="M34" display="inline"><mml:mrow><mml:mi>j</mml:mi><mml:mo>≠</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:math></inline-formula> weighted by the coefficients <inline-formula><mml:math id="M35" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> an <inline-formula><mml:math id="M36" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula>, which are determined statistically by minimizing <inline-formula><mml:math id="M37" display="inline"><mml:mrow><mml:mo>min⁡</mml:mo><mml:mfenced open="{" close="}"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo><mml:mo>-</mml:mo><mml:mover accent="true"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow><mml:mo mathvariant="normal" stretchy="true">^</mml:mo></mml:mover><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:mfenced></mml:mrow></mml:math></inline-formula> for a large number of observations. If RFI is present, the observation <inline-formula><mml:math id="M38" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:math></inline-formula> will be larger than the estimation <inline-formula><mml:math id="M39" display="inline"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="true" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>[</mml:mo><mml:mi>i</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:math></inline-formula>. Note that this technique assumes that the observations from other frequency channels are RFI-clean. Nonetheless, note that this is not always the case. This could lead to RFI contamination from other frequency channels leaking as false signatures in the desired channel. We apply this detection algorithm to all assessed instruments.</p></list-item></list></p>
</sec>
<sec id="Ch1.S2.SS3">
  <label>2.3</label><title>Threshold computation procedure</title>
      <p id="d2e1650">The EORFIScan has three threshold levels. Each is defined through a False Alarm Probability <inline-formula><mml:math id="M40" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. The higher <inline-formula><mml:math id="M41" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, the higher the number of detections, but their certainty decreases since the number of possible missed detections increases. In contrast, the lower <inline-formula><mml:math id="M42" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, the lower the number of detections, but the higher the certainty of detection. From now on, the first case will be referred to as less strict and the latter as more strict in terms on how certain we are that a detection is not a false alarm. Thresholds might depend on geometrical or geophysical parameters such as the latitude or the incidence angle, the polarization or the terrain type. For instance, the physical temperature tends to be higher at low latitudes, and colder at higher ones, so it is expected the brightness temperature to follow this tendency, or in a similar way, the natural range of brightness temperatures over sea is different than the range over land.</p>
      <p id="d2e1686">Since one of the applications is to use the cleaned data to feed geophysical models, thresholds must be determined statistically and be model independent. Otherwise, the data used to feed models would only contain the data that already fit the models, preventing thus any improvement in these same models.</p>
      <p id="d2e1689">The time period used for the threshold computation has been determined empirically. For AMSR2, one week of each consecutive month of 2021 was used, to take into account the variability between seasons. For AMSU-A and MWHS-2, which have less dense data sampling, two weeks of each month were used. AMR-C data is even more scarcely sampled, so two full years (2020 and 2021) were used. Each instrument's thresholds have been computed separately using uniquely its own data.</p>
      <p id="d2e1692">The ideal procedure for setting the RFI detection technique thresholds is to first apply the RFI detection techniques to a large number of clean observations. Then, to estimate the probability density function <inline-formula><mml:math id="M43" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> for each frequency band, polarization, and surface type (e.g. land, sea, coast). Then, a threshold <inline-formula><mml:math id="M44" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> defined by a given a priori probability of false alarm <inline-formula><mml:math id="M45" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is found using Eq. (<xref ref-type="disp-formula" rid="Ch1.E7"/>) <xref ref-type="bibr" rid="bib1.bibx40" id="paren.28"/>:

            <disp-formula id="Ch1.E7" content-type="numbered"><label>7</label><mml:math id="M46" display="block"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="normal">∞</mml:mi></mml:munderover><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></disp-formula></p>
      <p id="d2e1780">The a priori probability of false alarm is defined as erroneously choosing hypothesis <inline-formula><mml:math id="M47" display="inline"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> (RFI is present) when the true hypothesis is <inline-formula><mml:math id="M48" display="inline"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> (no RFI is present), which can be formulated as <inline-formula><mml:math id="M49" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>P</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>H</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>|</mml:mo><mml:msub><mml:mi>H</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>. The only exception to this procedure is the Generalized RFI Index, which first estimates the coefficients to compute the index, and then the thresholds are set to the PDF of the resulting error. In reality, it is impossible to exclude all samples contaminated with RFI, since this is actually the goal of the study. These samples modify the PDF used in Eq. (<xref ref-type="disp-formula" rid="Ch1.E7"/>). To minimize this effect, the strongest and most obvious contaminated samples and the nearby ones are excluded with a preliminary RFI screening method.</p>
      <p id="d2e1838">The function <inline-formula><mml:math id="M50" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> might not only depend on the random variable <inline-formula><mml:math id="M51" display="inline"><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, but also on the geometrical variables <inline-formula><mml:math id="M52" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>N</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> (latitude, longitude, incidence angle, azimuth angle, etc.). Ideally, the (<inline-formula><mml:math id="M53" display="inline"><mml:mrow><mml:mi>N</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>)-dimensional histogram of <inline-formula><mml:math id="M54" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>N</mml:mi></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> would be computed and the threshold <inline-formula><mml:math id="M55" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> would be computed for the desired <inline-formula><mml:math id="M56" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> for each <inline-formula><mml:math id="M57" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula>-dimensional bin <inline-formula><mml:math id="M58" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mi>N</mml:mi></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>:

            <disp-formula id="Ch1.E8" content-type="numbered"><label>8</label><mml:math id="M59" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub><mml:mo>=</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>Y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mi>N</mml:mi></mml:msub></mml:mrow></mml:mfenced></mml:mrow><mml:mi mathvariant="normal">∞</mml:mi></mml:munderover><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>N</mml:mi></mml:msub><mml:mo>|</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mi>N</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mi>N</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula></p>
      <p id="d2e2115">Estimating small <inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> values for all bins would require a huge amount of data, or else the threshold values for low <inline-formula><mml:math id="M61" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">fa</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> values would be quite noisy. Instead, to simplify, it will be assumed that <inline-formula><mml:math id="M62" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> depends at most on a single variable <inline-formula><mml:math id="M63" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> that stands for the most relevant geometrical variable. This variable is determined by inspecting the 2D histograms of <inline-formula><mml:math id="M64" display="inline"><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> against all possible <inline-formula><mml:math id="M65" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, taking into account that there might be no dependency. It also will be assumed that <inline-formula><mml:math id="M66" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> only changes the mean value of <inline-formula><mml:math id="M67" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> with a deterministic function <inline-formula><mml:math id="M68" display="inline"><mml:mrow><mml:mi>g</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>:

            <disp-formula id="Ch1.E9" content-type="numbered"><label>9</label><mml:math id="M69" display="block"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mi>g</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>

      <fig id="F1" specific-use="star"><label>Figure 1</label><caption><p id="d2e2282"><bold>(a)</bold> AMSR2 RFI pre-filtering map at 6.9 GHz and H polarization. Color axis shows the standard deviation of the high-pass filter value. <bold>(b)</bold> RFI Contamination maps observed by <xref ref-type="bibr" rid="bib1.bibx11" id="text.29"/>.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f01.jpg"/>

        </fig>

      <fig id="F2"><label>Figure 2</label><caption><p id="d2e2301">Biased 2-D histogram of the brightness temperature at 6.9 GHz V polarization in AMSR2 as function of the latitude.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f02.png"/>

        </fig>

      <fig id="F3" specific-use="star"><label>Figure 3</label><caption><p id="d2e2313"><bold>(a)</bold> Unbiased 2-D histogram of the brightness temperature at 6.9 GHz V polarization in AMSR2 as function of the latitude. The histogram is normalized in bins of 0.25° latitude. <bold>(b)</bold> Integrated 1-D histogram for the unbiased 2-D histogram in panel <bold>(a)</bold>. The figure has been rotated and vertically aligned to match the integration of panel <bold>(a)</bold> along the latitudinal axis aligning their peaks.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f03.png"/>

        </fig>

      <p id="d2e2333">The procedure for determining thresholds consists in first estimating <inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:mi>g</mml:mi><mml:mo>(</mml:mo><mml:mi>y</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> so that <inline-formula><mml:math id="M71" display="inline"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> becomes the probability density function of a unidimensional random variable with deterministic dependency on the function <inline-formula><mml:math id="M72" display="inline"><mml:mrow><mml:mi>g</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, then remove this dependency, and last compute the thresholds with respect this unidimensional random variable. The steps for this procedure are S.1 to S.9 described below.
          <def-list>
            <def-item><term>S.1</term><def>

      <p id="d2e2384">First, the data from AMSR2, AMSU-A and MWHS-2 was pre-filtered for strong RFI contamination using the high pass filter from the Image Enhancement algorithm. AMR-C data are used without RFI prefiltering. To do so, the filter is applied and the results are collocated in cells of 0.5° <inline-formula><mml:math id="M73" display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> 0.5° in latitude and longitude. The standard deviation of each cell is computed, and the cells are classified according to their surface type (“land”, “sea-only”, “coast”, “sea-ice”, “cloud/storm”). For a given terrain type “land” and “sea-only”, cells are discarded if they are considered an outlier following the criteria <inline-formula><mml:math id="M74" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">σ</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>≥</mml:mo><mml:mi mathvariant="normal">median</mml:mi><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:mn mathvariant="normal">4</mml:mn><mml:mo>⋅</mml:mo><mml:mi>c</mml:mi><mml:mo>⋅</mml:mo><mml:mi mathvariant="normal">MAD</mml:mi><mml:mo mathvariant="italic">{</mml:mo><mml:mi>X</mml:mi><mml:mo mathvariant="italic">}</mml:mo></mml:mrow></mml:math></inline-formula>, where <inline-formula><mml:math id="M75" display="inline"><mml:mrow><mml:mi mathvariant="normal">MAD</mml:mi><mml:mo mathvariant="italic">{</mml:mo><mml:mi>X</mml:mi><mml:mo mathvariant="italic">}</mml:mo><mml:mo>=</mml:mo><mml:mi mathvariant="normal">median</mml:mi><mml:mo>(</mml:mo><mml:mo>|</mml:mo><mml:mi>X</mml:mi><mml:mo>-</mml:mo><mml:mi mathvariant="normal">median</mml:mi><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mo>)</mml:mo><mml:mo>|</mml:mo><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is the Median Absolute Deviation, <inline-formula><mml:math id="M76" display="inline"><mml:mrow><mml:mi>c</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mrow><mml:msqrt><mml:mn mathvariant="normal">2</mml:mn></mml:msqrt><mml:mo>⋅</mml:mo><mml:msup><mml:mi mathvariant="normal">erfc</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup><mml:mo>(</mml:mo><mml:mn mathvariant="normal">3</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:math></inline-formula> <xref ref-type="bibr" rid="bib1.bibx27" id="paren.30"/>, <inline-formula><mml:math id="M77" display="inline"><mml:mrow><mml:mi mathvariant="normal">erfc</mml:mi><mml:mo>(</mml:mo><mml:mi>z</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mi mathvariant="normal">erf</mml:mi><mml:mo>(</mml:mo><mml:mi>z</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is the complimentary error function where <inline-formula><mml:math id="M78" display="inline"><mml:mrow><mml:mi mathvariant="normal">erf</mml:mi><mml:mo>(</mml:mo><mml:mi>z</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mn mathvariant="normal">2</mml:mn><mml:msqrt><mml:mi mathvariant="italic">π</mml:mi></mml:msqrt></mml:mfrac></mml:mstyle><mml:msubsup><mml:mo>∫</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mi>z</mml:mi></mml:msubsup><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:msup><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:math></inline-formula>. Figure <xref ref-type="fig" rid="F1"/>a shows the standard deviation of the high-pass filter in each cell for AMSR2 at 6.9 GHz H Polarization. Note how it matches the literature in Fig. <xref ref-type="fig" rid="F1"/>b <xref ref-type="bibr" rid="bib1.bibx11" id="paren.31"/>. Likely false positives due to sea-ice on the southern coast of Alaska, the Baltic Sea, and the Russian-Pacific coast are ignored, since they are “coast” pixels.</p>
            </def></def-item>
            <def-item><term>S.2</term><def>

      <p id="d2e2604">2-D histograms of each RFI detection technique (except the Generalized RFI Index), frequency band, and polarization are computed as a function of the available geometric parameters (latitude, incidence angle, etc.) to visualize which one has more influence. For instance, Fig. <xref ref-type="fig" rid="F2"/> shows the brightness temperature as function of the latitude. The dependencies found between thresholds and geometric parameters for each RFI detection techniques is mentioned in the list of RFI detection techniques in Sect. <xref ref-type="sec" rid="Ch1.S2.SS2"/>. The histograms are normalized in bins of the variable <inline-formula><mml:math id="M79" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> (the <inline-formula><mml:math id="M80" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula>-axis). In the case of the brightness temperature for ASMR2, the bins are 0.25° wide.</p>

      <p id="d2e2629">In the case of the Generalized RFI Index, once the coefficients <inline-formula><mml:math id="M81" display="inline"><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M82" display="inline"><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are estimated in Eq. (<xref ref-type="disp-formula" rid="Ch1.E6"/>), the 1-D histogram of the error in Eq. (<xref ref-type="disp-formula" rid="Ch1.E5"/>) is computed. In the case no geometrical dependencies are found, also the 1-D histogram is computed. Both cases jump directly to S.5.</p>
            </def></def-item>
            <def-item><term>S.3</term><def>

      <p id="d2e2664">The contribution <inline-formula><mml:math id="M83" display="inline"><mml:mrow><mml:mi>g</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is estimated as <inline-formula><mml:math id="M84" display="inline"><mml:mrow><mml:mover accent="true"><mml:mi>g</mml:mi><mml:mo mathvariant="normal" stretchy="true">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mo mathvariant="italic">{</mml:mo><mml:mi>x</mml:mi><mml:mo>:</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mtext> is maximum</mml:mtext><mml:mo mathvariant="italic">}</mml:mo></mml:mrow></mml:math></inline-formula>, and the unbiased 2-D histogram is computed as <inline-formula><mml:math id="M85" display="inline"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub></mml:mrow><mml:mo mathvariant="normal" stretchy="true">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>-</mml:mo><mml:mover accent="true"><mml:mi>g</mml:mi><mml:mo stretchy="true" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> as shown in Fig. <xref ref-type="fig" rid="F3"/>a.</p>
            </def></def-item>
            <def-item><term>S.4</term><def>

      <p id="d2e2807">The 2-D histogram in Fig. <xref ref-type="fig" rid="F3"/>a is integrated along the <inline-formula><mml:math id="M86" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> dimension (latitude) and divided by the total area to obtain a 1-D Probability Density Function as shown in Fig. <xref ref-type="fig" rid="F3"/>b.</p>
            </def></def-item>
            <def-item><term>S.5</term><def>

      <p id="d2e2827">The three desired thresholds <inline-formula><mml:math id="M87" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M88" display="inline"><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi mathvariant="normal">…</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula> (dashed lines in Fig. <xref ref-type="fig" rid="F3"/>b) and a reference threshold <inline-formula><mml:math id="M89" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> (continuous line in Fig. <xref ref-type="fig" rid="F3"/>b) are computed for a given set of <inline-formula><mml:math id="M90" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">fa</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M91" display="inline"><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi mathvariant="normal">…</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M92" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">fa</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> using Eq. (<xref ref-type="disp-formula" rid="Ch1.E7"/>). The reference false alarm probability is a high probability <inline-formula><mml:math id="M93" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">fa</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> used to smooth the threshold function as will be shown later. The other false alarm probabilities are set to <inline-formula><mml:math id="M94" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">fa</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">4</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M95" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">fa</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M96" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">fa</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">2.5</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, respectively. In the case of 1-D histograms only the three thresholds <inline-formula><mml:math id="M97" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M98" display="inline"><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mi mathvariant="normal">…</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula> are computed, and the procedure finishes here.</p>
            </def></def-item>
            <def-item><term>S.6</term><def>

      <p id="d2e3068">The distances between the thresholds to the reference threshold <inline-formula><mml:math id="M99" display="inline"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are calculated.</p>
            </def></def-item>
            <def-item><term>S.7</term><def>

      <p id="d2e3102">The reference threshold <inline-formula><mml:math id="M100" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> as function <inline-formula><mml:math id="M101" display="inline"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is computed in the 2-D biased histogram using Eq. (<xref ref-type="disp-formula" rid="Ch1.E10"/>), and is shown with a continuous red line in Fig. <xref ref-type="fig" rid="F4"/>.

                  <disp-formula id="Ch1.E10" content-type="numbered"><label>10</label><mml:math id="M102" display="block"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">fa</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub></mml:mrow></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mi mathvariant="normal">∞</mml:mi></mml:munderover><mml:msub><mml:mi>f</mml:mi><mml:mi>X</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:msub><mml:mi>x</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></disp-formula></p>
            </def></def-item>
            <def-item><term>S.8</term><def>

      <p id="d2e3220">The threshold approximation function for the reference false alarm probability <inline-formula><mml:math id="M103" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is approximated with a <inline-formula><mml:math id="M104" display="inline"><mml:mi>M</mml:mi></mml:math></inline-formula> order polynomial, shown in Fig. <xref ref-type="fig" rid="F4"/> as a continuous green line:

                  <disp-formula id="Ch1.E11" content-type="numbered"><label>11</label><mml:math id="M105" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>≈</mml:mo><mml:mi>p</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:msubsup><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn><mml:mi>M</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>⋅</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:math></disp-formula></p>
            </def></def-item>
            <def-item><term>S.9</term><def>

      <p id="d2e3332">The threshold approximation function for each threshold level is defined as <inline-formula><mml:math id="M106" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mi>p</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and is shown in Fig. <xref ref-type="fig" rid="F4"/> as a dashed green lines.</p>
            </def></def-item>
          </def-list></p>

      <fig id="F4"><label>Figure 4</label><caption><p id="d2e3381">Biased 2-D histogram of the brightness temperature at 6.9 GHz V polarization in AMSR2 as function of the latitude, with the reference threshold <inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="normal">ref</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> in red, the threshold approximation function for the reference false alarm probability in continuous green line, and the threshold approximation function for each threshold level in green dashed lines. The histogram is normalized in bins of 0.25° latitude.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f04.png"/>

        </fig>

      <p id="d2e3410">If Eqs. (<xref ref-type="disp-formula" rid="Ch1.E10"/>) and (<xref ref-type="disp-formula" rid="Ch1.E11"/>) were used to directly compute <inline-formula><mml:math id="M108" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, a large dataset would be required since <inline-formula><mml:math id="M109" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="normal">fa</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> are small, and small datasets will imply excessively noisy estimations. Still, large datasets will not ensure that two polynomial approximation functions <inline-formula><mml:math id="M110" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">α</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:msub><mml:mi>y</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> do not cross each other. The procedure described here reduces the amount of data needed, and ensures that threshold approximation functions do not cross each other, and therefore increasing false alarm probabilities always ensures a less strict threshold as defined at the beginning of this subsection.</p>
      <p id="d2e3472">Note that this subsection focused on the probability of false alarm and does not mention the a priori probability of missed detections <inline-formula><mml:math id="M111" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">md</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. Calculating <inline-formula><mml:math id="M112" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">md</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> requires knowing the statistics on the RFI sources. These statistics are scarce, vary over time, have been conducted in particular frequency bands (not necessarily generalized to other bands), and barely include low power RFIs, the ones more likely to be miss detected. As a consequence, this computation is not possible or may be inaccurate.</p>
      <p id="d2e3498">Note that since no ground truth is available, it is not possible to estimate the a posteriori probability of false alarms or the a posteriori probability of missed detections.</p>
</sec>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Results</title>
      <p id="d2e3511">In order to evaluate the status of the RFI environment, RFI probability maps are used. Neither AMSR2, AMSU-A or MWHS-2 data is gridded, so EASE2 grid with 25 km spacing has been used for the database. Observations are assigned to the closest grid point. Then, the probability of RFI is calculated at each grid point as <inline-formula><mml:math id="M113" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">RFI</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>D</mml:mi><mml:mo>/</mml:mo><mml:mi>O</mml:mi></mml:mrow></mml:math></inline-formula>, where <inline-formula><mml:math id="M114" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> is the number of RFI detections, and <inline-formula><mml:math id="M115" display="inline"><mml:mi>O</mml:mi></mml:math></inline-formula> is the total number of observations.</p>
      <p id="d2e3547">The maps are computed for each month along the temporal scope of each sensor analysis using both ascending and descending data, and both H and V polarization (if available). In the following, we present maps of mean RFI probabilities generated from these monthly maps, with the aim to highlight areas that are consistently flagged as RFI-affected. While these average maps may not capture more transient events that are only present during certain times, they provide a representative overview of the main affected areas. Average maps are a good trade-off between showing RFI contamination and mitigating false alarms.</p>
      <p id="d2e3550">Furthermore, after discarding measurements potentially affected by sun glint, observations over 350 K were stored along with their latitude, longitude and date per each frequency band. Then, they were grouped into clusters using <inline-formula><mml:math id="M116" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula>-means method.</p>
      <p id="d2e3560">Figures <xref ref-type="fig" rid="F5"/> to <xref ref-type="fig" rid="F9"/> show mean RFI probability maps for selected channels from the four instruments considered, covering the frequency range 6.9 to 183 GHz. Maps from AMSR2 have been prioritized over AMSU-A since two polarizations are available, and consequently the resulting maps are less noisy. Furthermore, the higher the frequency, the noisier the maps, since the spatial resolution is higher, the atmospheric effects on the <inline-formula><mml:math id="M117" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi mathvariant="normal">B</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> become more intense, and the atmospheric contributions are less smooth than the land or sea dynamics particularly because cloud and precipitation sensitivity generally increases with frequency.</p>

      <fig id="F5" specific-use="star"><label>Figure 5</label><caption><p id="d2e3581">RFI Probability maps for: <bold>(a)</bold> 6.9 GHz (AMSR2), <bold>(b)</bold> 7.3 GHz (AMSR2), and <bold>(c)</bold> 10.7 GHz (AMSR2).</p></caption>
        <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f05.jpg"/>

      </fig>

      <fig id="F6" specific-use="star"><label>Figure 6</label><caption><p id="d2e3601">RFI Probability maps for: <bold>(a)</bold> 18.7 GHz (AMSR2), <bold>(b)</bold> 23.8 GHz (AMSR2) and <bold>(c)</bold> 31.4 GHz (AMSU-A).</p></caption>
        <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f06.jpg"/>

      </fig>

      <fig id="F7" specific-use="star"><label>Figure 7</label><caption><p id="d2e3621">RFI Probability maps for: <bold>(a)</bold> 34 GHz (AMR-C), <bold>(b)</bold> 36.5 GHz (AMSR2), and <bold>(c)</bold> 50.3 GHz (AMSU-A).</p></caption>
        <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f07.jpg"/>

      </fig>

      <fig id="F8" specific-use="star"><label>Figure 8</label><caption><p id="d2e3641">RFI Probability maps for: <bold>(a)</bold> 52.6–59.3 GHz (AMSU-A), <bold>(b)</bold> 89 GHz (AMSR2), and <bold>(c)</bold> 118 GHz (MWHS-2).</p></caption>
        <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f08.jpg"/>

      </fig>

      <fig id="F9" specific-use="star"><label>Figure 9</label><caption><p id="d2e3662">RFI Probability maps for: <bold>(a)</bold> 166 GHz (MWHS-2), and <bold>(b)</bold> 183 GHz (MWHS-2).</p></caption>
        <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f09.png"/>

      </fig>

</sec>
<sec id="Ch1.S4">
  <label>4</label><title>Discussion of the Identified RFI Environment</title>
      <p id="d2e3685">The 6.4–7.7 GHz band is not exclusively allocated for passive applications and therefore passive sensors cannot claim protection. This regulatory framework can be linked to the significant contamination documented previously <xref ref-type="bibr" rid="bib1.bibx11 bib1.bibx46 bib1.bibx47" id="paren.32"/>. It is separated into the 6.9 and 7.3 GHz sub-bands. At the 6.9 GHz sub-band contamination is widespread in the US, and India. Some contamination is scattered in South America (with a denser detection area in Brazil), in Africa (with a denser detection area in South Africa), Europe (note the dense detection area over Ukraine and the detection from possible oil rigs in the North Sea), and the sea contamination in the Yellow Sea in front of China probably due to marine traffic. Point sources of 6.9 GHz contamination over sea had previously not been identified in the literature. In the 7.3 GHz sub-band contamination is present in all populated continents. Both North and South America show a few scattered detections, with strong detections in Central America. Europe also shows a few RFI sources, mostly in the North Sea due to oil rigs, and a large contaminated area over Ukraine. The Middle East is heavily contaminated in Turkey and Syria, with some contamination seen in Saudi Arabia. Indonesia and Thailand show strong signatures across their entire areas due to a telecommunication service using this band. Some isolated RFI sources can be seen in the rest of Oceania. Also note the reflection of geostationary satellite signals over the sea near Japan, in the Mediterranean Sea and the Atlantic Ocean near Europe, and in the Pacific Ocean between Hawaii and the US.</p>
      <p id="d2e3691">Only a part of the 10.7 GHz band is exclusive for passive services (10.68 to 10.70 GHz), whereas AMSR2 has a nominal bandwidth from 10.60 to 10.70 GHz. RFI contamination is mostly present in Europe due to the reflection of geostationary TV broadcast satellites over the sea <xref ref-type="bibr" rid="bib1.bibx41" id="paren.33"/>. The reflections are seen as stripes following the descending orbit direction. The same signatures can be seen as descending stripes in the Bay of Bengal near the coast, which were not observed in previous analyses. Some strong signatures can be seen in the UK, and some contamination is visible along the Nile River, in Nigeria, on the east coast of China, Mexico, Middle East, Central and Western Europe, Australia and widespread in Japan and Turkey.</p>
      <p id="d2e3697">The 18.7 GHz frequency band is not exclusive for Earth Observation, but despite this the RFI contamination maps are relatively clean. RFI contamination is found on the US coast, also including Hawaii. No other regions seem to be contaminated. However, the analysis focused on 2022; since then, strong contamination has been observed by the recent Compact Ocean Wind Vector Receiver onboard the International Space Station (COWVR-ISS) instrument <xref ref-type="bibr" rid="bib1.bibx4" id="paren.34"/>. It is true that COWVR-ISS has a larger bandwidth than AMSR2, and the detected RFIs could be out of the AMSR2 bandwidth. This fact plus a growing presence of contaminated areas makes an RFI detection and mitigation system increasingly attractive at this frequency band.</p>

      <fig id="F10" specific-use="star"><label>Figure 10</label><caption><p id="d2e3706">AMSR-2 RFI Probability maps centred at the Norwegian coastline for: <bold>(a)</bold> 6.9 GHz, <bold>(b)</bold> 7.3 GHz, <bold>(d)</bold> 18.7 GHz, and <bold>(e)</bold> 36.5 GHz; and <bold>(c)</bold> Oil Rig distribution near the Norwegian coast in 2022 <xref ref-type="bibr" rid="bib1.bibx31" id="paren.35"/>.</p></caption>
        <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f10.jpg"/>

      </fig>

      <p id="d2e3734">One area where we can zoom in is near the Norwegian coast in the North Sea and Arctic Ocean, where several spots of contamination at different bands can be observed. They closely match the location of the 2022 oil rig maps observed in Fig. <xref ref-type="fig" rid="F10"/> and it seems reasonable to view this as causal. These signatures are only observed by AMSR2. They are also visible at 6.9, 7.3 and 36.5 GHz, but not at 10.7 GHz, and 89 GHz (see Sect. <xref ref-type="sec" rid="Ch1.S3"/>). A possible explanation could be that the RFI Index technique is causing false positives due to the coefficients associated with certain bands. Bands where RFI is erroneously detected (tentatively 18.7 GHz) are explained by the coefficients in <inline-formula><mml:math id="M118" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M119" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula> corresponding to the <inline-formula><mml:math id="M120" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> of the channels where the RFI is present being negative, contributing positively to the RFI Index (see Eq. <xref ref-type="disp-formula" rid="Ch1.E5"/>). On the contrary, in bands that do not show this “cross-contamination” effect, the coefficients in <inline-formula><mml:math id="M121" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M122" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula> corresponding to the <inline-formula><mml:math id="M123" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> of the channels where the RFI is present is positive, contributing negatively to the RFI Index. Figure <xref ref-type="fig" rid="F10"/> shows that the upper RFI detected area at 7.3 GHz and the lower RFI detected area at 6.9 GHz seem to be “cross-contaminating” the 18.7 and 36.5 GHz maps, since patterns have similar shapes. It could be the opposite way, but it is less likely that an oil rig is transmitting at 36.5 GHz, since commercial radio links are common at 6.9 and 7.3 GHz bands.</p>

      <fig id="F11" specific-use="star"><label>Figure 11</label><caption><p id="d2e3798">Map showing the location of all RFI at 23.8 GHz (red), and at 36.5 GHz (blue).</p></caption>
        <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f11.png"/>

      </fig>

      <p id="d2e3807">The 23.8 GHz frequency band is exclusive for passive services, but the deployment of communication services for 5G services in the nearby band 24.25–27.5 GHz might cause RFI contamination in this band. This band was identified for use by 5G in 2019 <xref ref-type="bibr" rid="bib1.bibx35" id="paren.36"/>. Analysis of the AMSR2 measurements at 23.8 GHz highlighted instances in which the measured brightness temperature was between 350 and 400 K which indicates presence of RFI. Some of those observations are suspected to be caused by sun glint. These measurements were scattered along a constant latitude that changed from month to month, and matched the regions prone to sun glint shown in <xref ref-type="bibr" rid="bib1.bibx44" id="paren.37"/> and Sect. 3 of <xref ref-type="bibr" rid="bib1.bibx25" id="paren.38"/> whose latitude also change along the year. Some of these observations even showed abnormally high values compared to the surrounding area in the higher frequency bands. These observations were discarded from the analysis. The rest of those observations over 350 K were grouped in clusters and are listed in Table <xref ref-type="table" rid="T3"/>. The largest cluster of RFI signatures is located in Yuzhong County, Lanzhou, Gansu, China. These observations are, to our knowledge, the first confirmed detection of RFI in the 23.8 GHz band. The COWVR-ISS also reported RFI in this frequency band <xref ref-type="bibr" rid="bib1.bibx4" id="paren.39"/>, but since the instrument's bandwidth is 75 MHz larger than the allocated band for passive observations, it is not possible to ensured if it is an in-band or out-of-band emission. Adjacent channels include mobile and fixed communications, and amateur communications. All RFI observations have been located in a World map in Fig. <xref ref-type="fig" rid="F11"/>. These detections are not observed in the RFI Probability maps since the number of contaminated observations is proportionally low to the total number of observations collocated to a cell. Furthermore, the smaller antenna swath makes adjacent cells not to contain any contaminated observations, and therefore they do not show the typical “bump” shape. As consequence of both, the effect of these observations in the maps can not be distinguished from the detection noise.</p>

<table-wrap id="T3"><label>Table 3</label><caption><p id="d2e3830">Centroid coordinates of each cluster of RFI signatures at 23.8 GHz band, sorted by number of observations.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="2">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="justify" colwidth="4.7cm"/>
     <oasis:thead>
       <oasis:row>
         <oasis:entry colname="col1">Coordinates</oasis:entry>
         <oasis:entry colname="col2" align="left">Notes</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">(Lat, Lon)</oasis:entry>
         <oasis:entry colname="col2" align="left"/>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">35.9° N, 104.2° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Yuzhong County, Lanzhou, Gansu, China</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">32.8° N, 35.0° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Haifa, Israel</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">27.9° N, 113.0° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Yuetang District, Xiangtan, Hunan, China</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">41.0° N, 28.7° E</oasis:entry>
         <oasis:entry colname="col2" align="left">West Istanbul, Turkey</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">39.4° N, 75.7° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Shufu County, Prefectura de Kasgar, Sinkiang, China</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1"><inline-formula><mml:math id="M124" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>6.2° N, 106.8° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Yakarta, Indonesia</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">52.1° N, 21.1° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Warsav, Poland</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">41.0° N, 29.2° E</oasis:entry>
         <oasis:entry colname="col2" align="left">East side of Istanbul, Turkey</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">47.3° N, 38.9° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Outskirts of Taganrog, Russia</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1"><inline-formula><mml:math id="M125" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>4.0° N, <inline-formula><mml:math id="M126" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>79.2° E</oasis:entry>
         <oasis:entry colname="col2" align="left">San Cayetano, Ecuador</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">57.0° N, 24.0° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Riga, Latvia</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1"><inline-formula><mml:math id="M127" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>14.1° N, <inline-formula><mml:math id="M128" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>72.3° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Choqueca, Peru</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">49.4° N, 11.3° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Leinburg, Germany</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">55.7° N, 37.6° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Moscow, Russia</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1"><inline-formula><mml:math id="M129" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>23.5° N, <inline-formula><mml:math id="M130" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>46.5° E</oasis:entry>
         <oasis:entry colname="col2" align="left">São Paulo, Brasil</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">52.4° N, 19.8° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Nowe Wymyśle, Polonia</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">40.0° N, 116.3° E</oasis:entry>
         <oasis:entry colname="col2" align="left">Pekin, China</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap>

      <p id="d2e4066">The 31.3–31.8 GHz is exclusive for passive services in region 2, but no clear signs of RFI contamination can be seen except than false positives in already previously mentioned areas such as sea-ice missed classifications and coastal up-welling areas. The band from 33.5–34.5 GHz is shared but results do not suggest man-made RFI contamination. RFI has been found in this frequency band with COWVR <xref ref-type="bibr" rid="bib1.bibx4" id="paren.40"/>, but the larger bandwidth of COWVR (1975 MHz) also covers part of the adjacent bands. Therefore, these reported sources cannot be confirmed as in-band.</p>
      <p id="d2e4072">The 36–37 GHz band is shared with other services. RFIs have been reported in this frequency band in November 2018 caused by a radar facility in Roi Namur in the Kwajalein atoll, Marshall Islands. The RFI signal only affected the gain of the Sentinel-3 instrument, and the radiometric data quality was not affected <xref ref-type="bibr" rid="bib1.bibx39" id="paren.41"/>. No clear RFI signatures were detected in the considered dataset. AMSR2 data showed 1 single observations in 2022 over 350 K which could not be linked to sun glint effects. The only remaining observation is located at 37.41° N, 37.61° E in Araban, Gaziantep province, Turkey. This RFI observation has been located in a World map in Fig. <xref ref-type="fig" rid="F11"/>.</p>
      <p id="d2e4080">At higher frequencies, there were no clear detections of RFI signals in the considered dataset, regardless of whether the bands are exclusive or shared. The signatures shown in the maps seem to respond to geophysical signatures rather than RFI contamination due to their large extension, location, and variability throughout the year. Most bands show false detections at the sea-ice borders as a result of missed classifications in both hemispheres, even after using a sea ice concentration algorithm that performs well in the ice marginal zone <xref ref-type="bibr" rid="bib1.bibx26" id="paren.42"/>. Threshold could have been adjusted to reduce the false positives in the ice edges, but it would had come at the expense of misdetections in this area. Although not many RFI are expected here, some ships cross this area and their presence should not be neglected. False detections can also be seen at the borders of the Greenland ice sheet, or in the Himalayan mountains. Many coastlines, such as the Mauritanian coast, the Angola-Namibia coast, the Chilean-Peruvian coast, the north-west coast of Australia, the Yemen-Oman coast, the Somalia coast, the Baltic Sea coastline, the Yellow Sea coastline, and the California coast show a high RFI probability in many frequency bands (such as 18.7, 23.8, 36.5, or 89 GHz) and instruments. <xref ref-type="bibr" rid="bib1.bibx5" id="text.43"/> shows correspondence between these areas and the so-called upwelling regions of the ocean, of which the four main ones are the Chilean-Peruvian coast, the California coast, the Angola-Namibia coast, and the Mauritanian coast. The water anomaly observed by this effect in the Canary and California currents can reach 4 K after long time averaging and could explain the detection in these areas <xref ref-type="bibr" rid="bib1.bibx20 bib1.bibx2" id="paren.44"/>. For analogous reasons to the sea ice edge false alarms, softening the thresholds might result in misdetections in marien regions such as the Yellow sea of the seas around Japan. Other false positives are observed in South America and Africa can be seen in all frequency bands such as the Lakes Titicaca and Poopó. The RFI probability changes over the year with a stable decreasing trend until mid-year, when it starts to increase again. This could be caused by a reduction on the water level during summer until the rainy period in winter.</p>
</sec>
<sec id="Ch1.S5">
  <label>5</label><title>Results Validation</title>
      <p id="d2e4101">NWP models compare millions of observed radiances with model-simulated radiances every day as part of the data assimilation procedure. For microwave radiances at ECMWF, this process employs the RTTOV-SCATT radiative transfer model <xref ref-type="bibr" rid="bib1.bibx18" id="paren.45"/>, simulating radiative effects from the surface, atmospheric gases, clouds, and precipitation – a so-called all-sky simulation. The observed radiance (<inline-formula><mml:math id="M131" display="inline"><mml:mi>O</mml:mi></mml:math></inline-formula>) is compared to the model's short-range forecast simulation of the observation, known as the model background (<inline-formula><mml:math id="M132" display="inline"><mml:mi>B</mml:mi></mml:math></inline-formula>), interpolated in time and space to the observation location and with the same viewing geometry. This is the first step of the assimilation of the observations, which will subsequently lead to updated model fields. In addition to offering information to update the model state, analysis of departures between observations and model (<inline-formula><mml:math id="M133" display="inline"><mml:mrow><mml:mi>O</mml:mi><mml:mo>-</mml:mo><mml:mi>B</mml:mi></mml:mrow></mml:math></inline-formula>) can be instructive for monitoring the quality of observations themselves. It is worth noting that in this study the full spatial resolution of radiances are analyzed (i.e. with no averaging applied), in contrast to typical pre-processing of observations before assimilation at ECMWF that includes averaging (superobbing) for most microwave data.</p>
      <p id="d2e4133">Simulations of the NWP model background do not include artificial sources. The background thus represents a realization of a pristine RF environment, one with natural sources of emission only. Analysis of departures can therefore aid in the identification and quantification of RFI sources, as the natural range of <inline-formula><mml:math id="M134" display="inline"><mml:mrow><mml:mi>O</mml:mi><mml:mo>-</mml:mo><mml:mi>B</mml:mi></mml:mrow></mml:math></inline-formula> variability is known very well in most cases. For example, over the ocean at low frequencies that are primarily sensitive to SST, <inline-formula><mml:math id="M135" display="inline"><mml:mrow><mml:mi>O</mml:mi><mml:mo>-</mml:mo><mml:mi>B</mml:mi></mml:mrow></mml:math></inline-formula> variations are typically within about 1 K. This provides a tight constraint on expected behavior, allowing abnormal radiative signals to stand out as likely contamination. Instantaneous departures can show these anomalous scenes, but the real power of NWP-based monitoring is in the aggregate, with signals able to emerge even in noisy sampling over many days and weeks. Crucially, RFI should manifest as a positive <inline-formula><mml:math id="M136" display="inline"><mml:mrow><mml:mi>O</mml:mi><mml:mo>-</mml:mo><mml:mi>B</mml:mi></mml:mrow></mml:math></inline-formula> signal only, whereas other mismatches between the model and observations should be distributed between positive and negative departures. For example, a misplaced cloud in the model background will show up at a cloud-sensitive frequency band as a positive blob next to a negative one, whereas a region of RFI should be overwhelmingly positive. Even if RFI is coincident with cloud and precipitation, the large amount of sampling afforded by NWP-based analysis means that these signals should still be clear in time-averaged maps and histograms.</p>
      <p id="d2e4172">In this validation analysis, only results from AMSR2 low frequency channels are presented, covering March and October of 2022. This is to focus on frequency bands in which significant contamination was found. Fuller analysis can be found in the technical report from <xref ref-type="bibr" rid="bib1.bibx12" id="paren.46"/>, and comprehensive analysis of AMSR2 performance in the ECMWF system can be found in <xref ref-type="bibr" rid="bib1.bibx25" id="text.47"/>.</p>
<sec id="Ch1.S5.SS1">
  <label>5.1</label><title>Definitions, Screening, and Instantaneous Departures</title>
      <p id="d2e4188">Here we define the departure with inclusion of a bias correction term, as calculated by the variational bias correction scheme used within the data assimilation system <xref ref-type="bibr" rid="bib1.bibx1" id="paren.48"/>. This is a method of removing persistent sources of bias between observations and model, including overall instrument bias against the model plus terms for angular and atmospheric thickness components among others. Inclusion of the channel-dependent BC term here is important for analyzing AMSR2, which has roughly 2–5 K mean biases against the ECMWF model at several channels <xref ref-type="bibr" rid="bib1.bibx13" id="paren.49"/>. This should also aid in isolating signals of contamination rather than bias. The departure is defined:

            <disp-formula id="Ch1.E12" content-type="numbered"><label>12</label><mml:math id="M137" display="block"><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">BC</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo>(</mml:mo><mml:mi>O</mml:mi><mml:mo>-</mml:mo><mml:mi>B</mml:mi><mml:mo>)</mml:mo><mml:mo>-</mml:mo><mml:mi mathvariant="normal">BC</mml:mi></mml:mrow></mml:math></disp-formula></p>
      <p id="d2e4223">To illustrate what can be seen in the departures, Fig. <xref ref-type="fig" rid="F12"/> shows AMSR2 observations at 18.7V, the model background simulated values, and the difference. East of the United States' east coast is a long frontal region of cloud and precipitation, resulting in positive and negative departures side-by-side. In contrast, the area around Florida shows relatively weak but persistent positive departures of roughly 5–10 K, consistent with a reflected geostationary broadcast signal. In this way, RFI contamination can be identified by eye in <inline-formula><mml:math id="M138" display="inline"><mml:mrow><mml:mi>O</mml:mi><mml:mo>-</mml:mo><mml:mi>B</mml:mi></mml:mrow></mml:math></inline-formula> fields.</p>

      <fig id="F12"><label>Figure 12</label><caption><p id="d2e4242">Observed radiances (top), model background equivalents (middle), and background departures (<inline-formula><mml:math id="M139" display="inline"><mml:mrow><mml:mi>O</mml:mi><mml:mo>-</mml:mo><mml:mi>B</mml:mi></mml:mrow></mml:math></inline-formula>; bottom) from AMSR2 18.7V on 2 October 2022 over North America.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f12.jpg"/>

        </fig>

      <p id="d2e4264">However, the fidelity of model-simulated radiances is not the same everywhere and it is therefore necessary to introduce some screening of departures for statistical analysis. For example, microwave radiances over ice-free oceans have very well-modeled surface emissivities at the frequencies discussed, whereas more challenging surfaces such as sea-ice, high orography, and snow-covered ground have much poorer emissivity modeling. Certain regions are thus screened out of the validation analysis in order to exclude areas of known model bias or deficiencies in the forward modeling. This is why frozen surfaces of Canada, high orography of the Rocky Mountains, and some coastal scenes are not shown in Fig. <xref ref-type="fig" rid="F12"/>. This data selection essentially follows the screening rules for microwave imagers assimilated at ECMWF <xref ref-type="bibr" rid="bib1.bibx19" id="paren.50"/>. For more details on the data selection, see <xref ref-type="bibr" rid="bib1.bibx12" id="paren.51"/>.</p>
      <p id="d2e4275">RFI sources are spectrally narrow in most cases and should therefore be distinct between radiometry bands. This is important for departure-based RFI detection because model assumptions and radiative transfer modeling will be very similar at nearby frequencies. Consequently, a forward model error such as an erroneous surface emissivity will manifest similar at nearby frequencies whereas RFI will not. Figure <xref ref-type="fig" rid="F13"/> shows departures from the three lowest frequency bands on AMSR2 at V over Western Europe. The left panels show all departures, while the right panels show only points that EORFIScan flagged as likely RFI. As the figure presents descending nodes of the GCOM-W orbit and AMSR2 scans in the forward direction, it is possible to see reflected RFI from geostationary sources, particularly over sea at 10.65V in the Bay of Biscay, the North Sea, and in parts of the Mediterranean. It is also clear that some surfaces over land are simulated better than others, with larger model biases apparent over complex surfaces and some artefacts along coastlines. Regardless, several point sources of RFI are evident through localized positive departures over both land and sea, most evident at 6.9V in the North Sea, Luxembourg, and southern Spain, at 7.3V at different sites in the North Sea and Italy, and then very strong sources at 10.65V in the UK and Italy. Most of the clear examples of RFI in the left panels are identified by EORFIScan in the right panels. One possible exception is 7.3V contamination in Italy where despite both 6.9V and 10.65V RFI in the vicinity, the much stronger departures at 7.3 GHz suggest likely RFI here.</p>

      <fig id="F13"><label>Figure 13</label><caption><p id="d2e4282">Departures for AMSR2 channels 1, 3, and 5 (all V) on 2 March 2022 over western Europe. The right column only shows departures identified by EORFIScan as likely contamination.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f13.jpg"/>

        </fig>

</sec>
<sec id="Ch1.S5.SS2">
  <label>5.2</label><title>Statistical Evaluation</title>
      <p id="d2e4300">Moving from qualitative to quantitative evaluation of flagged departures, it is best to focus on areas in which there is known contamination to improve the signal to noise ratio of the analysis. For example, it is not useful to analyze the Southern Ocean for RFI flagging when it is clear from Figs. <xref ref-type="fig" rid="F5"/> and <xref ref-type="fig" rid="F6"/> that little to no real contamination is found there. The focus here is on a region encompassing North America, Europe, North Africa, and the Middle East (20–70° N, 140° W–50° E), as this region includes a variety of sources at different frequency bands. Figure <xref ref-type="fig" rid="F14"/> shows mean departure maps of this region for March 2022, gridded at a 1° resolution. The four lowest-frequency V channels from AMSR2 are shown, now with all screened departures in the top panels and flagged departures only in the bottom ones. In this analysis, the two higher-confidence flag levels (introduced in Sect. <xref ref-type="sec" rid="Ch1.S2.SS3"/>) are used together as simply “flagged data” for simplicity.</p>

      <fig id="F14" specific-use="star"><label>Figure 14</label><caption><p id="d2e4313">Mean departures from AMSR2 channels 6.9V, 7.3V, 10.65V, and 18.7V comprising data from March 2022. Plots alternate between all observations that passed the screening criteria (top) and RFI-flagged observations only (bottom). At least 1000 observations need to exist in a grid box for it to be plotted, and at least 10 flagged observations for the bottom panels.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f14.png"/>

        </fig>

      <p id="d2e4322">Due to atmospheric and surface variability, well-identified RFI will not exclusively manifest as positive departures, but will generally produce positive departures in the aggregate. Strong and/or persistent contamination should thus be red in the figure. The monthly mean departures for all data do indicate some regions for which frequent RFI can be picked out by eye without any filtering: e.g. Morocco and Ukraine at 7.3V or England and Italy at 10.65V. These are clearly regions with persistent contamination. For flagged points only, the expectation is that skillful detection would identify positive departures in the aggregate – lighter red for weak or infrequent true RFI, and dark red for consistently strong RFI. Most of the flagged points do match with positive mean departures, with larger magnitudes typical over land and weaker but more widespread RFI over sea. Some of the main features found by EORFIScan in Fig. <xref ref-type="fig" rid="F5"/>  –  e.g. reflected contamination at 10.65V around Europe and 18.7V around the USA  –  are verified by the departure analysis. Many smaller features that are either not persistent or not strong enough to be clear in Fig. <xref ref-type="fig" rid="F5"/> are visible in Fig. <xref ref-type="fig" rid="F14"/>, including the Nile Valley at 10.65V, Southern California at 6.9V, Spain at 7.3V, and so on.</p>

      <fig id="F15" specific-use="star"><label>Figure 15</label><caption><p id="d2e4334">Histograms of departures for screened observations over land (red; <inline-formula><mml:math id="M140" display="inline"><mml:mo>&gt;</mml:mo></mml:math></inline-formula> 95 % land fraction) and sea (blue; <inline-formula><mml:math id="M141" display="inline"><mml:mo>&lt;</mml:mo></mml:math></inline-formula> 1 % land fraction), comprising data from March 2022. The histograms of all points are given by solid lines, with the flagged histograms given by dotted lines. Region is the same as for the previous figure.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f15.png"/>

        </fig>

      <p id="d2e4357">Along the same lines, aggregate departures can be analysed across this region by separating points by land and sea. Figure <xref ref-type="fig" rid="F15"/> provides histograms of departures from this region, showing all data and flagged data from land and sea separately. This allows investigation of relative thresholds of detection and viewing the symmetry of the departures. It is clear that EORFIScan flags a large majority of very positive departures at 6.9 and 7.3 GHz over both land and sea, with nearly all departures larger than <inline-formula><mml:math id="M142" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>30 K flagged as likely contamination at 6.9V, for example. In these cases, it is quite likely that these are indeed RFI, as the departure PDFs have long positive tails. In general, EORFIScan flags very few negative departures, especially over sea, and this indicates a low rate of false positives. It is interesting to note that the histograms are significantly wider for H channels than for V at all frequencies, as might be expected from the larger dynamic range at H. Possibly as a result of this, EORFIScan is more cautious in flagging for H channels. This means that detection thresholds appear larger at H channels, as can be seen in the histogram peaks for flagged data over sea especially. One further feature to note is that EORFIScan thinks that most contamination over sea is fairly weak, with most cases being between about <inline-formula><mml:math id="M143" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>5 and <inline-formula><mml:math id="M144" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>20 K departures; this is in contrast to much of the contamination over land, where <inline-formula><mml:math id="M145" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>20 K signals are common and 7.3 GHz contamination is typically a <inline-formula><mml:math id="M146" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>25 K signal or more.</p>
</sec>
<sec id="Ch1.S5.SS3">
  <label>5.3</label><title>Validation Discussion and Critique</title>
      <p id="d2e4406">One aspect that is very important for data assimilation is that RFI screening does not remove real, natural signals of clouds and precipitation. This is more of a concern for higher frequency channels due to hydrometeor sensitivity generally increasing with frequency. The challenge for good RFI identification is to retain the natural signals while having a detection threshold that is low enough for the application of interest. Judging the performance of EORFIScan is not quite as easy as simply ensuring that the departure standard deviation decreases globally – if the flagging removes too many points with real geophysical signals, then it is undesirable.</p>

      <fig id="F16" specific-use="star"><label>Figure 16</label><caption><p id="d2e4411">Standard deviation of departures for all screened observations (top) and just observations with no flag (bottom) from AMSR2 channels 6.9V and 10.65V, comprising data from October 2022. At least 1000 (screened) observations need to exist in a grid box for it to be plotted.</p></caption>
          <graphic xlink:href="https://amt.copernicus.org/articles/19/4141/2026/amt-19-4141-2026-f16.png"/>

        </fig>

      <p id="d2e4420">With this in mind, Fig. <xref ref-type="fig" rid="F16"/> shows the standard deviation of departures before and after flagged data have been removed at 6.9V and 10.65V, now over October 2022 rather than March. Comparing the top and bottom panels, there are regions in which likely contamination appears almost entirely removed by EORFIScan flagging, for example over the UK at 10.65V or the northeast US at 6.9V. Some regions like southeast Ukraine do appear much improved after screening at 6.9V, but significant contamination remains. Outside of the known areas of contamination, the standard deviations of departures appear almost entirely unchanged; examples include the northwestern Atlantic in the Gulf Stream and the Sahara and Arabian deserts, in which large variability is expected and indeed retained by EORFIScan. This indicates that false positives are not significantly dampening the natural variability of departures.</p>
      <p id="d2e4426">Balancing detection against false positives is difficult to assess for RFI because there is not a truth dataset to compare against, and assessing the detection threshold is similarly challenging because low-level RFI is particularly hard to identify. In our analysis, negative departures represent likely (but not certain) false positives if they are flagged. With this definition of false positives, EORFIScan has a false positive rate of between 0.05 % to 0.2 %, depending on the frequency. Given this low rate of false positives, it could be suggested that the flagging is primarily identifying larger instances of RFI. Going back to the swath-level view of Fig. <xref ref-type="fig" rid="F13"/>, it does appear that flagging is rare near coastlines, and some points of likely RFI remain in the un-flagged maps of Fig. <xref ref-type="fig" rid="F16"/> along coastlines.</p>
      <p id="d2e4433">This balance is more difficult for H channels where the range of natural variability is larger because of more variable emissivity. The histograms of Fig. <xref ref-type="fig" rid="F15"/> suggest that detection limits of EORFIScan are perhaps twice or more what they are at V, though 10.65 GHz departures over sea are an interesting outlier with most flagged departures lying between about <inline-formula><mml:math id="M147" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>8 to 15 K at both polarizations. It is hard to say anything definitive about the strength of expected contamination because these could conceivably be quite different at V and H at the same frequency band. But if the assumption is that RFI magnitudes are similar across polarizations then the detection skill is stronger at V over both land and sea, especially at C-band.</p>
</sec>
</sec>
<sec id="Ch1.S6" sec-type="conclusions">
  <label>6</label><title>Conclusions</title>
      <p id="d2e4454">This paper described an assessment of key passive microwave bands exploited for Earth Observation from the 6.4–7.3 GHz band to 174–191 GHz. First the EORFIScan software has been described, which applies multiple RFI detection techniques to EO brightness temperature products. This part includes a description of the RFI detection techniques and the methodology to compute the thresholds. The software has then been used to analyze a year of AMSR2, AMSU-A, MWHS-2, and AMR-C data in order to cover all main EO passive microwave bands beyond the L-band. The methodology has been assessed for AMSR2 and AMSU-A via comparisons with NWP model equivalents, with results presented for the 6–19 GHz bands on AMSR2, where the large majority of corroborated RFI was found. RFI probability maps for 2022 have been computed as the number of detections over the number of observations accounting for all polarizations and orbital passes.</p>
      <p id="d2e4457">The results have shown heaviest contamination in the 6.4–7.3 GHz band, more moderate contamination at 10.7 and 18.7 GHz bands. The identified areas of systematic contamination have been assessed using all-sky NWP departures, and in most cases are found to correlate strongly with positive departures as would be expected for skilful flagging of RFI. Several new sources are observed in these bands with respect to the literature and corroborated by the NWP-based analysis. The emergence of 6.9 GHz contamination over sea is particularly troubling, as this is a key band for all-weather SST and sea-ice measurements <xref ref-type="bibr" rid="bib1.bibx41" id="paren.52"/>, and will be flown on the upcoming CIMR mission from ESA as well as on the successor instrument to AMSR2, AMSR3. First RFI signatures have been found at 23.8 and 36.5 GHz in AMSR2.</p>
      <p id="d2e4463">Consideration of RFI detection systems in these bands is strongly advised as part of future mission planning, as the contamination observed already has an impact on users and may be expected to worsen in the upcoming years. At bands from 50 GHz and up, no significant contamination has been detected in the dataset under consideration. However, the constantly evolving threat of RFI in EO bands would benefit greatly from monitoring in real time, especially for bands of consequence for NWP and other operational applications. It is clearly suboptimal to discover new sources of contamination a year or two after those radiances have been ingested into a variety of operational models and retrievals, and thus real-time monitoring would provide significant benefit to a wide community of users.</p>
      <p id="d2e4466">It is reasonable to expect that high-frequency technology will become more accessible, low-frequency bands will become more crowded, and new services will require greater bandwidths in the near future. As a result, higher frequency bands may start showing contamination and RFI at lower frequency bands could become more widespread than it already is, which will threaten their use for Earth observation. Given that new missions take years from design until the end of their operational lifetime, new EO missions should include RFI detection and mitigation systems to be prepared for the expected worsening of the RFI environment.</p>
</sec>

      
      </body>
    <back><notes notes-type="codedataavailability"><title>Code and data availability</title>

      <p id="d2e4473">The threshold coefficients, the RFI Index coefficients, the functions to compute the thresholds and the RFI Index, and an example script to generate them can be found in <ext-link xlink:href="https://doi.org/10.5281/zenodo.19727331" ext-link-type="DOI">10.5281/zenodo.19727331</ext-link> <xref ref-type="bibr" rid="bib1.bibx34" id="paren.53"/>. Data used for AMSR2, AMSU-A and MWHS-II are in BUFR format for convenience of ECMWF. The data in this format is not publicly available, however it is publicly available in other data formats such as HDF5, which containing additional auxiliary data not included in the BUFR ones. AMSR2 data can be found in <uri>https://doi.org/10.57746/EO.01GS73ANS548QGHAKNZDJYXD2H</uri> <xref ref-type="bibr" rid="bib1.bibx23" id="paren.54"/>, AMSU-A data can be found in <ext-link xlink:href="https://doi.org/10.7289/V5X63JT2" ext-link-type="DOI">10.7289/V5X63JT2</ext-link> <xref ref-type="bibr" rid="bib1.bibx50" id="paren.55"/>, MWHS-II can be found in <xref ref-type="bibr" rid="bib1.bibx48" id="text.56"/>, and AMR-C data can be downloaded from <ext-link xlink:href="https://doi.org/10.15770/EUM_SEC_CLM_0116" ext-link-type="DOI">10.15770/EUM_SEC_CLM_0116</ext-link> <xref ref-type="bibr" rid="bib1.bibx16" id="paren.57"/>.</p>
  </notes><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d2e4507">Ra.O., J.B., I.N., and A.J. wrote the software; Ra.O. processed the data; D.D. and N.B. validated the results; Ra.O., Ro.O., D.D., and N.B. interpreted the results; Ra.O. and D.D. graphic results generation; Ra.O., Ro.O., J.B., D.D., and N.B. wrote the resulting technical reports; Y.S., J.C., and F.J. reviewed the technical reports; Ra.O., Ro.O., D.D., N.B. and J.B. requested the funding; Ra.O., Ro.O., J.B. and N.B. managed the project; Y.S., J.C., and F.J. supervised the project; Ra.O. and D.D. wrote and edited the manuscript; all authors reviewed the manuscript.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d2e4513">The contact author has declared that none of the authors has any competing interests.</p>
  </notes><notes notes-type="disclaimer"><title>Disclaimer</title>

      <p id="d2e4519">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="d2e4525">The authors acknowledge the Japan Aerospace Exploration Agency (JAXA) for the GCOM-W BUFR data provided through ECMWF, to European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) for the AMSU-A and AMR-C data, and to the Chineese National Satellite Meteorological Center (NSMC) for the MWHS-II data. The authors also acknowledge Antonio Martelucci and Stephen English for his help and comments during the project.</p></ack><notes notes-type="financialsupport"><title>Financial support</title>

      <p id="d2e4530">This project was developed under ESA contract no. 4000141035/23/NL/SD.</p>
  </notes><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d2e4536">This paper was edited by Laura Bianco and reviewed by five anonymous referees.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bibx1"><label>Auligné et al.(2007)</label><mixed-citation>Auligné, T., McNally, A., and Dee, D.: Adaptive bias correction for satellite data in a numerical weather prediction system, Quart. J. Roy. Meteor. Soc., 133, 631–642, <ext-link xlink:href="https://doi.org/10.1002/qj.56" ext-link-type="DOI">10.1002/qj.56</ext-link>, 2007.</mixed-citation></ref>
      <ref id="bib1.bibx2"><label>Benazzouz et al.(2014)</label><mixed-citation>Benazzouz, A., Mordane, S., Orbi, A., Chagdali, M., Hilmi, K., Atillah, A., Lluís Pelegrí, J., and Hervé, D.: An improved coastal upwelling index from sea surface temperature using satellite-based approach – The case of the Canary Current upwelling system, Continental Shelf Research, 81, 38–54, <ext-link xlink:href="https://doi.org/10.1016/j.csr.2014.03.012" ext-link-type="DOI">10.1016/j.csr.2014.03.012</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bibx3"><label>Bringer et al.(2018)</label><mixed-citation>Bringer, A., Daehn, M., Johnson, J. T., Soldo, Y., Le Vine, D. M., de Matthaeis, P., Piepmeier, J. R., and Mohammed, P.: SMAP Mission: Changes in the RFI Environment, IEEE, 3754–3757, <ext-link xlink:href="https://doi.org/10.1109/igarss.2018.8519067" ext-link-type="DOI">10.1109/igarss.2018.8519067</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bibx4"><label>Brown and Morris(2024)</label><mixed-citation>Brown, S. and Morris, M.: Compact Ocean Wind Vector Radiometer (COWVR) Environmental Data Record (EDR) Quick-start User's Guide, Jet Propoulsion Laboratory (JPL) from California Institute of Technology (CALTECH), Pasadena, CA, version 1.02, <uri>https://archive.podaac.earthdata.nasa.gov/podaac-ops-cumulus-docs/cowvr-tempest/open/docs/STP_H8_COWVR_User-Guide_v6.pdf</uri> (last access: 22 June 2026), 2024.</mixed-citation></ref>
      <ref id="bib1.bibx5"><label>Capone and Hutchins(2013)</label><mixed-citation>Capone, D. G. and Hutchins, D. A.: Microbial biogeochemistry of coastal upwelling regimes in a changing ocean, Nat. Geosci., 6, 711–717, <ext-link xlink:href="https://doi.org/10.1038/ngeo1916" ext-link-type="DOI">10.1038/ngeo1916</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx6"><label>Casey et al.(2022)</label><mixed-citation>Casey, A., Gagnon, N., Bedard, J., Beaulne, A., Shahabadi, M. B., and Aparicio, J.: Potential impact of RFI on NWP: a Canadian perspective, in: RFI, ECMWF, <uri>https://events.ecmwf.int/event/258/contributions/2884/attachments/1564/2803/RFI2022_Casey_Beaulne.pdf</uri> (last access: 24 April 2024), 2022.</mixed-citation></ref>
      <ref id="bib1.bibx7"><label>Cavalieri et al.(1984)</label><mixed-citation>Cavalieri, D. J., Gloersen, P., and Campbell, W. J.: Determination of sea ice parameters with the NIMBUS 7 SMMR, J. Geophys. Res.-Atmos., 89, 5355–5369, <ext-link xlink:href="https://doi.org/10.1029/JD089iD04p05355" ext-link-type="DOI">10.1029/JD089iD04p05355</ext-link>, 1984.</mixed-citation></ref>
      <ref id="bib1.bibx8"><label>CIA(1985)</label><mixed-citation>CIA: World Data Bank II: North America, South America, Europe, Africa, Asia, Inter-university Consortium for Political and Social Research, <ext-link xlink:href="https://doi.org/10.3886/ICPSR08376.V1" ext-link-type="DOI">10.3886/ICPSR08376.V1</ext-link>, 1985.</mixed-citation></ref>
      <ref id="bib1.bibx9"><label>Draper(2016a)</label><mixed-citation>Draper, D. W.: Report on GMI Special Study #15: Radio Frequency Interference, Technical Report 20160003316, Goddard Space Flight Center, Greenbelt, USA, <uri>https://ntrs.nasa.gov/api/citations/20160003316/downloads/20160003316.pdf</uri> (last access: 22 June 2026), 2016a.</mixed-citation></ref>
      <ref id="bib1.bibx10"><label>Draper(2016b)</label><mixed-citation>Draper, D. W.: Terrestrial and space-based RFI observed by the GPM microwave imager (GMI) within NTIA semi-protected passive earth exploration bands at 10.65 and 18.7 GHz, in: 2016 Radio Frequency Interference (RFI), 26–30, <ext-link xlink:href="https://doi.org/10.1109/RFINT.2016.7833526" ext-link-type="DOI">10.1109/RFINT.2016.7833526</ext-link>, 2016b.</mixed-citation></ref>
      <ref id="bib1.bibx11"><label>Draper(2018)</label><mixed-citation>Draper, D. W.: Radio Frequency Environment for Earth-Observing Passive Microwave Imagers, IEEE J. Sel. Top. Appl., 11, 1913–1922, <ext-link xlink:href="https://doi.org/10.1109/JSTARS.2018.2801019" ext-link-type="DOI">10.1109/JSTARS.2018.2801019</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bibx12"><label>Duncan and Bormann(2024)</label><mixed-citation>Duncan, D. and Bormann, N.: Assessing RFI flags at passive microwave bands with an NWP model, Tech. rep., ESA Contract Report, Reading, <ext-link xlink:href="https://doi.org/10.21957/eefd6f0954" ext-link-type="DOI">10.21957/eefd6f0954</ext-link>, 2024.</mixed-citation></ref>
      <ref id="bib1.bibx13"><label>Duncan et al.(2024)</label><mixed-citation>Duncan, D. I., Geer, A., Bormann, N., and Dahoui, M.: Vicarious calibration monitoring of MWI and ICI using NWP fields, Tech. rep., EUMETSAT Contract Report, <ext-link xlink:href="https://doi.org/10.21957/7c2d18d2e1" ext-link-type="DOI">10.21957/7c2d18d2e1</ext-link>, 2024.</mixed-citation></ref>
      <ref id="bib1.bibx14"><label>English et al.(2020)</label><mixed-citation>English, S., Duncan, D., and Turner, E.: Assessment of Value of Provision of Spectral Response Functions from CGMS Agencies, in: CGMS-48 WG II, <uri>https://www.researchgate.net/publication/351358361_Assessment_of_Value_of_Provision_of_Spectral_Response_Functions_from_CGMS_Agencies</uri> (last access: 22 June 2026), 2020.</mixed-citation></ref>
      <ref id="bib1.bibx15"><label>ERMIT(2025)</label><mixed-citation>ERMIT: ESA RFI Monitoring &amp; Information Tool, <uri>https://rfi.smos.eo.esa.int/rfimanager/publicindex.jsp</uri>, (last access: 15 April 2025), 2025.</mixed-citation></ref>
      <ref id="bib1.bibx16"><label>EUMETSAT for Copernicus(2025)</label><mixed-citation>EUMETSAT for Copernicus: Advanced Microwave Radiometer (AMR-C) Level 2 Products (baseline version G01) – Sentinel-6 – Reprocessed, V. G01, EUMETSAT [data set], <ext-link xlink:href="https://doi.org/10.15770/EUM_SEC_CLM_0116" ext-link-type="DOI">10.15770/EUM_SEC_CLM_0116</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bibx17"><label>Farrar and Jones(2016)</label><mixed-citation>Farrar, S. and Jones, L.: Corruption of the TRMM microwave imager cold sky mirror due to RFI, in: 2016 14th Specialist Meeting on Microwave Radiometry and Remote Sensing of the Environment (MicroRad), 76–81, <ext-link xlink:href="https://doi.org/10.1109/MICRORAD.2016.7530508" ext-link-type="DOI">10.1109/MICRORAD.2016.7530508</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bibx18"><label>Geer et al.(2021)</label><mixed-citation>Geer, A. J., Bauer, P., Lonitz, K., Barlakas, V., Eriksson, P., Mendrok, J., Doherty, A., Hocking, J., and Chambon, P.: Bulk hydrometeor optical properties for microwave and sub-millimetre radiative transfer in RTTOV-SCATT v13.0, Geosci. Model Dev., 14, 7497–7526, <ext-link xlink:href="https://doi.org/10.5194/gmd-14-7497-2021" ext-link-type="DOI">10.5194/gmd-14-7497-2021</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bibx19"><label>Geer et al.(2022)</label><mixed-citation>Geer, A. J., Lonitz, K., Duncan, D. I., and Bormann, N.: Improved Surface Treatment for All-sky Microwave Observations, Tech. Rep. 894, ECMWF Tech. Memo., Shinfield Park, Reading, <ext-link xlink:href="https://doi.org/10.21957/zi7q6hau" ext-link-type="DOI">10.21957/zi7q6hau</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bibx20"><label>Gómez-Ocampo et al.(2018)</label><mixed-citation>Gómez-Ocampo, E., Gaxiola-Castro, G., Durazo, R., and Beier, E.: Effects of the 2013-2016 warm anomalies on the California Current phytoplankton, Deep-Sea Res. Pt. II, 151, 64–76, <ext-link xlink:href="https://doi.org/10.1016/j.dsr2.2017.01.005" ext-link-type="DOI">10.1016/j.dsr2.2017.01.005</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bibx21"><label>Gonzalez and Woods(2018)</label><mixed-citation> Gonzalez, R. C. and Woods, R. E.: Intensity Transformations and Spatial Filtering, Chap. 3, 4th global edn., Pearson, ISBN 978-1-292-22307-0, 2018.</mixed-citation></ref>
      <ref id="bib1.bibx22"><label>Imaoka et al.(2010)</label><mixed-citation> Imaoka, K., Kachi, M., Kasahara, M., Ito, N., Nakagawa, K., and Oki, T.: Instrument performance and calibration of AMSR-E and AMSR2, Int. Arch. Photogramm., 38, 13–18, 2010.</mixed-citation></ref>
      <ref id="bib1.bibx23"><label>JAXA(2012)</label><mixed-citation>Japan Aerospace Exploration Agency (JAXA): GCOM-W/AMSR2 L1B Brightness Temperature, Japan Aerospace Exploration Agency [data set], <uri>https://doi.org/10.57746/EO.01GS73ANS548QGHAKNZDJYXD2H</uri>, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx24"><label>Jorge et al.(2023)</label><mixed-citation> Jorge, F., Pedro, L., and Mendonça, S.: Emerging prospects for Earth observation in addressing sustainable development goals, ITU News Magazine, 5, 50–54, 2023.</mixed-citation></ref>
      <ref id="bib1.bibx25"><label>Kazumori et al.(2016)</label><mixed-citation>Kazumori, M., Geer, A. J., and English, S. J.: Effects of all-sky assimilation of GCOM-W/AMSR2 radiances in the ECMWF numerical weather prediction system, Quart. J. Roy. Meteor. Soc., <ext-link xlink:href="https://doi.org/10.1002/qj.2669" ext-link-type="DOI">10.1002/qj.2669</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bibx26"><label>Kern et al.(2019)</label><mixed-citation>Kern, S., Lavergne, T., Notz, D., Pedersen, L. T., Tonboe, R. T., Saldo, R., and Sørensen, A. M.: Satellite passive microwave sea-ice concentration data set intercomparison: closed ice and ship-based observations, The Cryosphere, 13, 3261–3307, <ext-link xlink:href="https://doi.org/10.5194/tc-13-3261-2019" ext-link-type="DOI">10.5194/tc-13-3261-2019</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bibx27"><label>Leys et al.(2013)</label><mixed-citation>Leys, C., Ley, C., Klein, O., Bernard, P., and Licata, L.: Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median, J. Exp. Soc. Psychol., 49, 764–766, <ext-link xlink:href="https://doi.org/10.1016/j.jesp.2013.03.013" ext-link-type="DOI">10.1016/j.jesp.2013.03.013</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx28"><label>Li et al.(2004)</label><mixed-citation>Li, L., Njoku, E., Im, E., Chang, P., and Germain, K.: A preliminary survey of radio-frequency interference over the U.S. in Aqua AMSR-E data, IEEE T. Geosci. Remote, 42, 380–390, <ext-link xlink:href="https://doi.org/10.1109/tgrs.2003.817195" ext-link-type="DOI">10.1109/tgrs.2003.817195</ext-link>, 2004.</mixed-citation></ref>
      <ref id="bib1.bibx29"><label>Maiwald et al.(2020)</label><mixed-citation>Maiwald, F., Brown, S. T., Koch, T., Milligan, L., Kangaslahti, P., Schlecht, E., Skalare, A., Bloom, M., Torossian, V., Kanis, J., Statham, S., Kang, S., and Vaze, P.: Completion of the AMR-C Instrument for Sentinel-6, IEEE J. Sel. Top. Appl., 13, 1811–1818, <ext-link xlink:href="https://doi.org/10.1109/JSTARS.2020.2991175" ext-link-type="DOI">10.1109/JSTARS.2020.2991175</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bibx30"><label>Mo(1996)</label><mixed-citation>Mo, T.: Prelaunch calibration of the advanced microwave sounding unit-A for NOAA-K, IEEE T. Microw. Theory, 44, 1460–1469, <ext-link xlink:href="https://doi.org/10.1109/22.536029" ext-link-type="DOI">10.1109/22.536029</ext-link>, 1996.</mixed-citation></ref>
      <ref id="bib1.bibx31"><label>NOD(2024)</label><mixed-citation>Norwegian Offshore Directorate (NOD): Maps of the Norweigean Continental Shelf, <uri>https://www.sodir.no/en/whats-new/publications/Maps-of-the-ncs</uri> (last access: 14 March 2025), 2024.</mixed-citation></ref>
      <ref id="bib1.bibx32"><label>Oliva et al.(2016)</label><mixed-citation>Oliva, R., Daganzo, E., Richaume, P., Kerr, Y., Cabot, F., Soldo, Y., Anterrieu, E., Reul, N., Gutierrez, A., Barbosa, J., and Lopes, G.: Status of Radio Frequency Interference (RFI) in the 1400–1427 MHz passive band based on six years of SMOS mission, Remote Sens. Environ., 180, 64–75, <ext-link xlink:href="https://doi.org/10.1016/j.rse.2016.01.013" ext-link-type="DOI">10.1016/j.rse.2016.01.013</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bibx33"><label>Oliva et al.(2021)</label><mixed-citation>Oliva, R., Onrubia, R., Martellucci, A., Daganzo-Eusebio, E., Jorge, F., Soldo, Y., English, S., De Rosnay, P., Weston, P., Barbosa, J., and Nestoras, I.: Results from the Ground RFI Detection System for Passive Microwave Earth Observation Data, in: 2021 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), IEEE, 1827–1830, <ext-link xlink:href="https://doi.org/10.1109/igarss47720.2021.9553636" ext-link-type="DOI">10.1109/igarss47720.2021.9553636</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bibx34"><label>Onrubia and Oliva(2026)</label><mixed-citation>Onrubia, R. and Oliva, R.: Threshold Generation Coefficients for the EORFIScan Software, Zenodo [data set], <ext-link xlink:href="https://doi.org/10.5281/zenodo.19727331" ext-link-type="DOI">10.5281/zenodo.19727331</ext-link>, 2026.</mixed-citation></ref>
      <ref id="bib1.bibx35"><label>Ookla(2023)</label><mixed-citation>Ookla: 5G in asia pacific: Deployment momentum continues, GSMA (GSM Association), <uri>https://www.gsma.com/get-involved/gsma-membership/gsma_resources/5g-in-asia-pacific-deployment-momentum-continues</uri> (last access: 22 June 2026), 2023.</mixed-citation></ref>
      <ref id="bib1.bibx36"><label>Parde et al.(2011)</label><mixed-citation>Parde, M., Zribi, M., Fanise, P., and Dechambre, M.: Analysis of RFI Issue Using the CAROLS L-Band Experiment, IEEE T. Geosci. Remote, 49, 1063–1070, <ext-link xlink:href="https://doi.org/10.1109/tgrs.2010.2069101" ext-link-type="DOI">10.1109/tgrs.2010.2069101</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bibx37"><label>Pedro et al.(2022)</label><mixed-citation>Pedro, L., Sá, M., Fernandes, R., Jorge, F., and Mendonça, S.: Protection of Earth Observation Satellites From Radio-Frequency Interference: Policies and Practices [Perspectives], IEEE Geoscience and Remote Sensing Magazine, 10, 278–288, <ext-link xlink:href="https://doi.org/10.1109/MGRS.2022.3221824" ext-link-type="DOI">10.1109/MGRS.2022.3221824</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bibx38"><label>Piepmeier et al.(2014)</label><mixed-citation>Piepmeier, J. R., Johnson, J. T., Mohammed, P. N., Bradley, D., Ruf, C., Aksoy, M., Garcia, R., Hudson, D., Miles, L., and Wong, M.: Radio-Frequency Interference Mitigation for the Soil Moisture Active Passive Microwave Radiometer, IEEE T. Geosci. Remote, 52, 761–775, <ext-link xlink:href="https://doi.org/10.1109/TGRS.2013.2281266" ext-link-type="DOI">10.1109/TGRS.2013.2281266</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bibx39"><label>Quartly et al.(2020)</label><mixed-citation>Quartly, G. D., Nencioli, F., Raynal, M., Bonnefond, P., Nilo Garcia, P., Garcia-Mondéjar, A., Flores de la Cruz, A., Crétaux, J.-F., Taburet, N., Frery, M.-L., Cancet, M., Muir, A., Brockley, D., McMillan, M., Abdalla, S., Fleury, S., Cadier, E., Gao, Q., Escorihuela, M. J., Roca, M., Bergé-Nguyen, M., Laurain, O., Bruniquel, J., Féménias, P., and Lucas, B.: The Roles of the S3MPC: Monitoring, Validation and Evolution of Sentinel-3 Altimetry Observations, Remote Sens., 12, <ext-link xlink:href="https://doi.org/10.3390/rs12111763" ext-link-type="DOI">10.3390/rs12111763</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bibx40"><label>Querol et al.(2017)</label><mixed-citation>Querol, J., Onrubia, R., Alonso-Arroyo, A., Pascual, D., Park, H., and Camps, A.: Performance Assessment of Time–Frequency RFI Mitigation Techniques in Microwave Radiometry, IEEE J. Sel. Top. Appl., 10, 3096–3106, <ext-link xlink:href="https://doi.org/10.1109/JSTARS.2017.2654541" ext-link-type="DOI">10.1109/JSTARS.2017.2654541</ext-link>, 2017. </mixed-citation></ref>
      <ref id="bib1.bibx41"><label>Scanlon et al.(2025)</label><mixed-citation>Scanlon, T., Geer, A. J., Bormann, N., and Duncan, D. I.: Using 6.9 and 10.65 GHz From the AMSR2 and GMI Microwave Imagers in the ECMWF NWP System (IFS), IEEE T. Geosci. Remote, 63, 1–19, <ext-link xlink:href="https://doi.org/10.1109/tgrs.2025.3627448" ext-link-type="DOI">10.1109/tgrs.2025.3627448</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bibx42"><label>SFCG(2014)</label><mixed-citation>SFCG: Report 32-1R2: Passive Sensor Filter Characteristics, Tech. rep., Space Frequency Coordination Group, <uri>https://sfcgonline.org/resources/reports/</uri> (last access: 22 June 2026), 2014.</mixed-citation></ref>
      <ref id="bib1.bibx43"><label>Spreen et al.(2008)</label><mixed-citation>Spreen, G., Kaleschke, L., and Heygster, G.: Sea ice remote sensing using AMSR-E 89-GHz channels, J. Geophys. Res.-Oceans, 113, <ext-link xlink:href="https://doi.org/10.1029/2005JC003384" ext-link-type="DOI">10.1029/2005JC003384</ext-link>, 2008.</mixed-citation></ref>
      <ref id="bib1.bibx44"><label>Wentz et al.(2021)</label><mixed-citation>Wentz, F., Meissner, T., Gentemann, C., Hilburn, K., and Scott, J.: RSS/AMSR-2 Data Images/Daily, Remote Sensing Systems, <uri>https://images.remss.com/amsr/amsr2_data_daily.html</uri> (last access: 23 April 2026), 2021.</mixed-citation></ref>
      <ref id="bib1.bibx45"><label>World Meteorological Organization (WMO)(2025)</label><mixed-citation>World Meteorological Organization (WMO): List of all Instruments, observing Systems Capability Analysis and Review Tool (OSCAR), <uri>https://space.oscar.wmo.int/instruments</uri> (last access: 22 June 2026), 2025.</mixed-citation></ref>
      <ref id="bib1.bibx46"><label>Wu et al.(2019)</label><mixed-citation>Wu, Y., Qian, B., Bao, Y., Li, M., Petropoulos, G. P., Liu, X., and Li, L.: Detection and Analysis of C-Band Radio Frequency Interference in AMSR2 Data over Land, Remote Sens., 11, <ext-link xlink:href="https://doi.org/10.3390/rs11101228" ext-link-type="DOI">10.3390/rs11101228</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bibx47"><label>Wu et al.(2020)</label><mixed-citation>Wu, Y., Li, M., Bao, Y., and Petropoulos, G. P.: Cross-Validation of Radio-Frequency-Interference Signature in Satellite Microwave Radiometer Observations over the Ocean, Remote Sensing, 12, <ext-link xlink:href="https://doi.org/10.3390/rs12203433" ext-link-type="DOI">10.3390/rs12203433</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bibx48"><label>Xian et al.(2021)</label><mixed-citation>Xian, D., Zhang, P., Gao, L., Sun, R., Zhang, H., and Jia, X.: Fengyun Meteorological Satellite Products for Earth System Science Applications, Adv. Atmos. Sci., 38, 1267–1284, <ext-link xlink:href="https://doi.org/10.1007/s00376-021-0425-3" ext-link-type="DOI">10.1007/s00376-021-0425-3</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bibx49"><label>Zhang et al.(2012)</label><mixed-citation>Zhang, S., Li, J., Wang, Z., Wang, H., Sun, M., Jiang, J., and He, J.: Design of the second generation microwave humidity sounder (MWHS-II) for Chinese meteorological satellite FY-3, in: 2012 IEEE International Geoscience and Remote Sensing Symposium, 4672–4675, <ext-link xlink:href="https://doi.org/10.1109/IGARSS.2012.6350423" ext-link-type="DOI">10.1109/IGARSS.2012.6350423</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx50"><label>Zou et al.(2013)</label><mixed-citation>Zou, C.-Z., Wang, W., and NOAA CDR Program: NOAA Fundamental Climate Data Record (FCDR) of AMSU-A Level 1c Brightness Temperature, Version 1.0, NOAA National Climatic Data Center [data set], <ext-link xlink:href="https://doi.org/10.7289/V5X63JT2" ext-link-type="DOI">10.7289/V5X63JT2</ext-link>, 2013.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>Assessment of the RFI environment in key passive microwave bands for Earth observation</article-title-html>
<abstract-html/>
<ref-html id="bib1.bib1"><label>Auligné et al.(2007)</label><mixed-citation>
      
Auligné, T., McNally, A., and Dee, D.: Adaptive bias correction for
satellite data in a numerical weather prediction system, Quart. J. Roy.
Meteor. Soc., 133, 631–642, <a href="https://doi.org/10.1002/qj.56" target="_blank">https://doi.org/10.1002/qj.56</a>, 2007.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>Benazzouz et al.(2014)</label><mixed-citation>
      
Benazzouz, A., Mordane, S., Orbi, A., Chagdali, M., Hilmi, K., Atillah, A.,
Lluís Pelegrí, J., and Hervé, D.: An improved coastal upwelling index from
sea surface temperature using satellite-based approach – The case of the
Canary Current upwelling system, Continental Shelf Research, 81, 38–54,
<a href="https://doi.org/10.1016/j.csr.2014.03.012" target="_blank">https://doi.org/10.1016/j.csr.2014.03.012</a>, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>Bringer et al.(2018)</label><mixed-citation>
      
Bringer, A., Daehn, M., Johnson, J. T., Soldo, Y., Le Vine, D. M., de Matthaeis, P., Piepmeier, J. R., and Mohammed, P.: SMAP Mission: Changes in the RFI Environment, IEEE, 3754–3757,
<a href="https://doi.org/10.1109/igarss.2018.8519067" target="_blank">https://doi.org/10.1109/igarss.2018.8519067</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>Brown and Morris(2024)</label><mixed-citation>
      
Brown, S. and Morris, M.: Compact Ocean Wind Vector Radiometer (COWVR) Environmental Data Record (EDR) Quick-start User's Guide, Jet Propoulsion Laboratory (JPL) from California Institute of Technology (CALTECH), Pasadena, CA, version 1.02, <a href="https://archive.podaac.earthdata.nasa.gov/podaac-ops-cumulus-docs/cowvr-tempest/open/docs/STP_H8_COWVR_User-Guide_v6.pdf" target="_blank"/> (last access: 22 June 2026), 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>Capone and Hutchins(2013)</label><mixed-citation>
      
Capone, D. G. and Hutchins, D. A.: Microbial biogeochemistry of coastal
upwelling regimes in a changing ocean, Nat. Geosci., 6, 711–717,
<a href="https://doi.org/10.1038/ngeo1916" target="_blank">https://doi.org/10.1038/ngeo1916</a>, 2013.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>Casey et al.(2022)</label><mixed-citation>
      
Casey, A., Gagnon, N., Bedard, J., Beaulne, A., Shahabadi, M. B., and Aparicio, J.: Potential impact of RFI on NWP: a Canadian perspective, in: RFI, ECMWF, <a href="https://events.ecmwf.int/event/258/contributions/2884/attachments/1564/2803/RFI2022_Casey_Beaulne.pdf" target="_blank"/> (last access: 24 April 2024), 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>Cavalieri et al.(1984)</label><mixed-citation>
      
Cavalieri, D. J., Gloersen, P., and Campbell, W. J.: Determination of sea ice parameters with the NIMBUS 7 SMMR, J. Geophys. Res.-Atmos., 89, 5355–5369, <a href="https://doi.org/10.1029/JD089iD04p05355" target="_blank">https://doi.org/10.1029/JD089iD04p05355</a>, 1984.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>CIA(1985)</label><mixed-citation>
      
CIA: World Data Bank II: North America, South America, Europe, Africa, Asia, Inter-university Consortium for Political and Social Research,
<a href="https://doi.org/10.3886/ICPSR08376.V1" target="_blank">https://doi.org/10.3886/ICPSR08376.V1</a>, 1985.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>Draper(2016a)</label><mixed-citation>
      
Draper, D. W.: Report on GMI Special Study #15: Radio Frequency Interference, Technical Report 20160003316, Goddard Space Flight Center, Greenbelt, USA, <a href="https://ntrs.nasa.gov/api/citations/20160003316/downloads/20160003316.pdf" target="_blank"/> (last access: 22 June 2026), 2016a.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>Draper(2016b)</label><mixed-citation>
      
Draper, D. W.: Terrestrial and space-based RFI observed by the GPM microwave imager (GMI) within NTIA semi-protected passive earth exploration bands at 10.65 and 18.7&thinsp;GHz, in: 2016 Radio Frequency Interference (RFI), 26–30,
<a href="https://doi.org/10.1109/RFINT.2016.7833526" target="_blank">https://doi.org/10.1109/RFINT.2016.7833526</a>, 2016b.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>Draper(2018)</label><mixed-citation>
      
Draper, D. W.: Radio Frequency Environment for Earth-Observing Passive Microwave Imagers, IEEE J. Sel. Top. Appl., 11, 1913–1922, <a href="https://doi.org/10.1109/JSTARS.2018.2801019" target="_blank">https://doi.org/10.1109/JSTARS.2018.2801019</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>Duncan and Bormann(2024)</label><mixed-citation>
      
Duncan, D. and Bormann, N.: Assessing RFI flags at passive microwave bands with an NWP model, Tech. rep., ESA Contract Report, Reading,
<a href="https://doi.org/10.21957/eefd6f0954" target="_blank">https://doi.org/10.21957/eefd6f0954</a>, 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>Duncan et al.(2024)</label><mixed-citation>
      
Duncan, D. I., Geer, A., Bormann, N., and Dahoui, M.: Vicarious calibration
monitoring of MWI and ICI using NWP fields, Tech. rep., EUMETSAT
Contract Report, <a href="https://doi.org/10.21957/7c2d18d2e1" target="_blank">https://doi.org/10.21957/7c2d18d2e1</a>, 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>English et al.(2020)</label><mixed-citation>
      
English, S., Duncan, D., and Turner, E.: Assessment of Value of Provision of Spectral Response Functions from CGMS Agencies, in: CGMS-48 WG II, <a href="https://www.researchgate.net/publication/351358361_Assessment_of_Value_of_Provision_of_Spectral_Response_Functions_from_CGMS_Agencies" target="_blank"/> (last access: 22 June 2026), 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>ERMIT(2025)</label><mixed-citation>
      
ERMIT: ESA RFI Monitoring &amp; Information Tool, <a href="https://rfi.smos.eo.esa.int/rfimanager/publicindex.jsp" target="_blank"/>, (last access: 15 April 2025), 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>EUMETSAT for Copernicus(2025)</label><mixed-citation>
      
EUMETSAT for Copernicus: Advanced Microwave Radiometer (AMR-C) Level 2 Products (baseline version G01) – Sentinel-6 – Reprocessed, V. G01, EUMETSAT [data set], <a href="https://doi.org/10.15770/EUM_SEC_CLM_0116" target="_blank">https://doi.org/10.15770/EUM_SEC_CLM_0116</a>, 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>Farrar and Jones(2016)</label><mixed-citation>
      
Farrar, S. and Jones, L.: Corruption of the TRMM microwave imager cold sky mirror due to RFI, in: 2016 14th Specialist Meeting on Microwave Radiometry and Remote Sensing of the Environment (MicroRad), 76–81,
<a href="https://doi.org/10.1109/MICRORAD.2016.7530508" target="_blank">https://doi.org/10.1109/MICRORAD.2016.7530508</a>, 2016.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>Geer et al.(2021)</label><mixed-citation>
      
Geer, A. J., Bauer, P., Lonitz, K., Barlakas, V., Eriksson, P., Mendrok, J., Doherty, A., Hocking, J., and Chambon, P.: Bulk hydrometeor optical properties for microwave and sub-millimetre radiative transfer in RTTOV-SCATT v13.0, Geosci. Model Dev., 14, 7497–7526, <a href="https://doi.org/10.5194/gmd-14-7497-2021" target="_blank">https://doi.org/10.5194/gmd-14-7497-2021</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>Geer et al.(2022)</label><mixed-citation>
      
Geer, A. J., Lonitz, K., Duncan, D. I., and Bormann, N.: Improved Surface Treatment for All-sky Microwave Observations, Tech. Rep. 894, ECMWF Tech. Memo., Shinfield Park, Reading, <a href="https://doi.org/10.21957/zi7q6hau" target="_blank">https://doi.org/10.21957/zi7q6hau</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>Gómez-Ocampo et al.(2018)</label><mixed-citation>
      
Gómez-Ocampo, E., Gaxiola-Castro, G., Durazo, R., and Beier, E.: Effects of the 2013-2016 warm anomalies on the California Current phytoplankton, Deep-Sea Res. Pt. II, 151, 64–76, <a href="https://doi.org/10.1016/j.dsr2.2017.01.005" target="_blank">https://doi.org/10.1016/j.dsr2.2017.01.005</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>Gonzalez and Woods(2018)</label><mixed-citation>
      
Gonzalez, R. C. and Woods, R. E.: Intensity Transformations and Spatial Filtering, Chap. 3, 4th global edn., Pearson, ISBN 978-1-292-22307-0, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>Imaoka et al.(2010)</label><mixed-citation>
      
Imaoka, K., Kachi, M., Kasahara, M., Ito, N., Nakagawa, K., and Oki, T.: Instrument performance and calibration of AMSR-E and AMSR2, Int. Arch. Photogramm., 38, 13–18, 2010.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>JAXA(2012)</label><mixed-citation>
      
Japan Aerospace Exploration Agency (JAXA): GCOM-W/AMSR2 L1B Brightness Temperature, Japan Aerospace Exploration Agency [data set], <a href="https://doi.org/10.57746/EO.01GS73ANS548QGHAKNZDJYXD2H" target="_blank"/>, 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>Jorge et al.(2023)</label><mixed-citation>
      
Jorge, F., Pedro, L., and Mendonça, S.: Emerging prospects for Earth observation in addressing sustainable development goals, ITU News Magazine, 5, 50–54, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>Kazumori et al.(2016)</label><mixed-citation>
      
Kazumori, M., Geer, A. J., and English, S. J.: Effects of all-sky assimilation of GCOM-W/AMSR2 radiances in the ECMWF numerical weather prediction system, Quart. J. Roy. Meteor. Soc., <a href="https://doi.org/10.1002/qj.2669" target="_blank">https://doi.org/10.1002/qj.2669</a>, 2016.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>Kern et al.(2019)</label><mixed-citation>
      
Kern, S., Lavergne, T., Notz, D., Pedersen, L. T., Tonboe, R. T., Saldo, R., and Sørensen, A. M.: Satellite passive microwave sea-ice concentration data set intercomparison: closed ice and ship-based observations, The Cryosphere, 13, 3261–3307, <a href="https://doi.org/10.5194/tc-13-3261-2019" target="_blank">https://doi.org/10.5194/tc-13-3261-2019</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>Leys et al.(2013)</label><mixed-citation>
      
Leys, C., Ley, C., Klein, O., Bernard, P., and Licata, L.: Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median, J. Exp. Soc. Psychol., 49, 764–766, <a href="https://doi.org/10.1016/j.jesp.2013.03.013" target="_blank">https://doi.org/10.1016/j.jesp.2013.03.013</a>, 2013.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib28"><label>Li et al.(2004)</label><mixed-citation>
      
Li, L., Njoku, E., Im, E., Chang, P., and Germain, K.: A preliminary survey of radio-frequency interference over the U.S. in Aqua AMSR-E data, IEEE T. Geosci. Remote, 42, 380–390, <a href="https://doi.org/10.1109/tgrs.2003.817195" target="_blank">https://doi.org/10.1109/tgrs.2003.817195</a>, 2004.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib29"><label>Maiwald et al.(2020)</label><mixed-citation>
      
Maiwald, F., Brown, S. T., Koch, T., Milligan, L., Kangaslahti, P., Schlecht, E., Skalare, A., Bloom, M., Torossian, V., Kanis, J., Statham, S., Kang, S., and Vaze, P.: Completion of the AMR-C Instrument for Sentinel-6, IEEE J. Sel. Top. Appl., 13, 1811–1818, <a href="https://doi.org/10.1109/JSTARS.2020.2991175" target="_blank">https://doi.org/10.1109/JSTARS.2020.2991175</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib30"><label>Mo(1996)</label><mixed-citation>
      
Mo, T.: Prelaunch calibration of the advanced microwave sounding unit-A for NOAA-K, IEEE T. Microw. Theory, 44, 1460–1469, <a href="https://doi.org/10.1109/22.536029" target="_blank">https://doi.org/10.1109/22.536029</a>, 1996.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib31"><label>NOD(2024)</label><mixed-citation>
      
Norwegian Offshore Directorate (NOD): Maps of the Norweigean Continental Shelf, <a href="https://www.sodir.no/en/whats-new/publications/Maps-of-the-ncs" target="_blank"/> (last access: 14 March 2025), 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib32"><label>Oliva et al.(2016)</label><mixed-citation>
      
Oliva, R., Daganzo, E., Richaume, P., Kerr, Y., Cabot, F., Soldo, Y., Anterrieu, E., Reul, N., Gutierrez, A., Barbosa, J., and Lopes, G.: Status of Radio Frequency Interference (RFI) in the 1400–1427&thinsp;MHz passive band based on six years of SMOS mission, Remote Sens. Environ., 180, 64–75, <a href="https://doi.org/10.1016/j.rse.2016.01.013" target="_blank">https://doi.org/10.1016/j.rse.2016.01.013</a>, 2016.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib33"><label>Oliva et al.(2021)</label><mixed-citation>
      
Oliva, R., Onrubia, R., Martellucci, A., Daganzo-Eusebio, E., Jorge, F., Soldo, Y., English, S., De Rosnay, P., Weston, P., Barbosa, J., and Nestoras, I.: Results from the Ground RFI Detection System for Passive Microwave Earth Observation Data, in: 2021 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), IEEE, 1827–1830, <a href="https://doi.org/10.1109/igarss47720.2021.9553636" target="_blank">https://doi.org/10.1109/igarss47720.2021.9553636</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib34"><label>Onrubia and Oliva(2026)</label><mixed-citation>
      
Onrubia, R. and Oliva, R.: Threshold Generation Coefficients for the EORFIScan Software, Zenodo [data set], <a href="https://doi.org/10.5281/zenodo.19727331" target="_blank">https://doi.org/10.5281/zenodo.19727331</a>, 2026.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib35"><label>Ookla(2023)</label><mixed-citation>
      
Ookla: 5G in asia pacific: Deployment momentum continues, GSMA (GSM Association), <a href="https://www.gsma.com/get-involved/gsma-membership/gsma_resources/5g-in-asia-pacific-deployment-momentum-continues" target="_blank"/> (last access: 22 June 2026), 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib36"><label>Parde et al.(2011)</label><mixed-citation>
      
Parde, M., Zribi, M., Fanise, P., and Dechambre, M.: Analysis of RFI Issue Using the CAROLS L-Band Experiment, IEEE T. Geosci. Remote, 49, 1063–1070, <a href="https://doi.org/10.1109/tgrs.2010.2069101" target="_blank">https://doi.org/10.1109/tgrs.2010.2069101</a>, 2011.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib37"><label>Pedro et al.(2022)</label><mixed-citation>
      
Pedro, L., Sá, M., Fernandes, R., Jorge, F., and Mendonça, S.: Protection of Earth Observation Satellites From Radio-Frequency Interference: Policies and Practices [Perspectives], IEEE Geoscience and Remote Sensing Magazine, 10,
278–288, <a href="https://doi.org/10.1109/MGRS.2022.3221824" target="_blank">https://doi.org/10.1109/MGRS.2022.3221824</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib38"><label>Piepmeier et al.(2014)</label><mixed-citation>
      
Piepmeier, J. R., Johnson, J. T., Mohammed, P. N., Bradley, D., Ruf, C., Aksoy, M., Garcia, R., Hudson, D., Miles, L., and Wong, M.: Radio-Frequency Interference Mitigation for the Soil Moisture Active Passive Microwave
Radiometer, IEEE T. Geosci. Remote, 52, 761–775,
<a href="https://doi.org/10.1109/TGRS.2013.2281266" target="_blank">https://doi.org/10.1109/TGRS.2013.2281266</a>, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib39"><label>Quartly et al.(2020)</label><mixed-citation>
      
Quartly, G. D., Nencioli, F., Raynal, M., Bonnefond, P., Nilo Garcia, P., Garcia-Mondéjar, A., Flores de la Cruz, A., Crétaux, J.-F., Taburet, N., Frery, M.-L., Cancet, M., Muir, A., Brockley, D., McMillan, M., Abdalla, S., Fleury, S., Cadier, E., Gao, Q., Escorihuela, M. J., Roca, M., Bergé-Nguyen, M., Laurain, O., Bruniquel, J., Féménias, P., and Lucas, B.: The Roles of the S3MPC: Monitoring, Validation and Evolution of Sentinel-3 Altimetry Observations, Remote Sens., 12, <a href="https://doi.org/10.3390/rs12111763" target="_blank">https://doi.org/10.3390/rs12111763</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib40"><label>Querol et al.(2017)</label><mixed-citation>
      
Querol, J., Onrubia, R., Alonso-Arroyo, A., Pascual, D., Park, H., and Camps, A.: Performance Assessment of Time–Frequency RFI Mitigation Techniques in Microwave Radiometry, IEEE J. Sel. Top. Appl., 10, 3096–3106,
<a href="https://doi.org/10.1109/JSTARS.2017.2654541" target="_blank">https://doi.org/10.1109/JSTARS.2017.2654541</a>, 2017.


    </mixed-citation></ref-html>
<ref-html id="bib1.bib41"><label>Scanlon et al.(2025)</label><mixed-citation>
      
Scanlon, T., Geer, A. J., Bormann, N., and Duncan, D. I.: Using 6.9 and 10.65&thinsp;GHz From the AMSR2 and GMI Microwave Imagers in the ECMWF NWP System (IFS), IEEE T. Geosci. Remote, 63, 1–19, <a href="https://doi.org/10.1109/tgrs.2025.3627448" target="_blank">https://doi.org/10.1109/tgrs.2025.3627448</a>, 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib42"><label>SFCG(2014)</label><mixed-citation>
      
SFCG: Report 32-1R2: Passive Sensor Filter Characteristics, Tech. rep., Space Frequency Coordination Group, <a href="https://sfcgonline.org/resources/reports/" target="_blank"/> (last access: 22 June 2026), 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib43"><label>Spreen et al.(2008)</label><mixed-citation>
      
Spreen, G., Kaleschke, L., and Heygster, G.: Sea ice remote sensing using
AMSR-E 89-GHz channels, J. Geophys. Res.-Oceans, 113,
<a href="https://doi.org/10.1029/2005JC003384" target="_blank">https://doi.org/10.1029/2005JC003384</a>, 2008.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib44"><label>Wentz et al.(2021)</label><mixed-citation>
      
Wentz, F., Meissner, T., Gentemann, C., Hilburn, K., and Scott, J.: RSS/AMSR-2 Data Images/Daily, Remote Sensing Systems, <a href="https://images.remss.com/amsr/amsr2_data_daily.html" target="_blank"/> (last access: 23 April 2026), 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib45"><label>World Meteorological Organization (WMO)(2025)</label><mixed-citation>
      
World Meteorological Organization (WMO): List of all Instruments, observing Systems Capability Analysis and Review Tool (OSCAR), <a href="https://space.oscar.wmo.int/instruments" target="_blank"/> (last access: 22 June 2026), 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib46"><label>Wu et al.(2019)</label><mixed-citation>
      
Wu, Y., Qian, B., Bao, Y., Li, M., Petropoulos, G. P., Liu, X., and Li, L.: Detection and Analysis of C-Band Radio Frequency Interference in AMSR2
Data over Land, Remote Sens., 11, <a href="https://doi.org/10.3390/rs11101228" target="_blank">https://doi.org/10.3390/rs11101228</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib47"><label>Wu et al.(2020)</label><mixed-citation>
      
Wu, Y., Li, M., Bao, Y., and Petropoulos, G. P.: Cross-Validation of Radio-Frequency-Interference Signature in Satellite Microwave Radiometer Observations over the Ocean, Remote Sensing, 12, <a href="https://doi.org/10.3390/rs12203433" target="_blank">https://doi.org/10.3390/rs12203433</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib48"><label>Xian et al.(2021)</label><mixed-citation>
      
Xian, D., Zhang, P., Gao, L., Sun, R., Zhang, H., and Jia, X.: Fengyun Meteorological Satellite Products for Earth System Science Applications, Adv. Atmos. Sci., 38, 1267–1284, <a href="https://doi.org/10.1007/s00376-021-0425-3" target="_blank">https://doi.org/10.1007/s00376-021-0425-3</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib49"><label>Zhang et al.(2012)</label><mixed-citation>
      
Zhang, S., Li, J., Wang, Z., Wang, H., Sun, M., Jiang, J., and He, J.: Design of the second generation microwave humidity sounder (MWHS-II) for Chinese meteorological satellite FY-3, in: 2012 IEEE International Geoscience and Remote Sensing Symposium, 4672–4675, <a href="https://doi.org/10.1109/IGARSS.2012.6350423" target="_blank">https://doi.org/10.1109/IGARSS.2012.6350423</a>, 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib50"><label>Zou et al.(2013)</label><mixed-citation>
      
Zou, C.-Z., Wang, W., and NOAA CDR Program: NOAA Fundamental Climate Data Record (FCDR) of AMSU-A Level 1c Brightness Temperature, Version 1.0, NOAA National Climatic Data Center [data set], <a href="https://doi.org/10.7289/V5X63JT2" target="_blank">https://doi.org/10.7289/V5X63JT2</a>, 2013.

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