Trends in spectrally resolved outgoing longwave radiation from 10 years of satellite measurements

In recent years, the interest has grown in satellite-derived hyperspectral radiance measurements for assessing the individual impact of climate drivers and their cascade of feedbacks on the outgoing longwave radiation (OLR). In this paper, we use 10 years (2008–2017) of reprocessed radiances from the infrared atmospheric sounding interferometer (IASI) to evaluate the linear trends in clear-sky spectrally resolved OLR (SOLR) in the range [645–2300] cm−1. Spatial inhomogeneities are observed in most of the analyzed spectral regions. These mostly reflected the natural variability of the atmospheric temperature and composition but long-term changes in greenhouse gases concentrations are also highlighted. In particular, the increase of atmospheric CO2 and CH4 led to significant negative trends in the SOLR of −0.05 to −0.3% per year in the spectral region corresponding to the ν2 and the ν3 CO2 and in the ν4 CH4 band. Most of the trends associated with the natural variability of the OLR can be related to the El Niño/Southern Oscillation activity and its teleconnections in the studied period. This is the case for the channels most affected by the temperature variations of the surface and the first layers of the atmosphere but also for the channels corresponding to the ν2 H2O and the ν3 O3 bands.


INTRODUCTION
The Earth's climate is closely linked to the flow of energy in and out of the Earth-atmosphere system. At equilibrium, the outgoing longwave radiation (OLR) and the reflected shortwave radiation (RSW) (W ⋅ m −2 ) at the top of the atmosphere (TOA) compensate the Incoming Solar Radiation (ISR) on average. Any perturbation of this balance due to a variation of the climate drivers (e.g. greenhouse gases), known as a radiative forcing, leads to a climate feedback which attempts to bring the radiation budget back to equilibrium 1 . A good understanding of the Earth-atmosphere system and of its long-term changes therefore requires an accurate evaluation of the radiative effects of the climate drivers and feedbacks on the OLR. Considerable improvements have been achieved in the last four decades through the exploitation of spectrally integrated OLR from dedicated broadband instruments (e.g. refs. [1][2][3][4][5][6]. Different methods, such as the radiative kernel technique 7,8 and the partial radiative perturbation 9 , have been developed to evaluate the contribution of individual variables to the changes in the OLR from broadband measurements. While these techniques are very powerful, they rely on an accurate knowledge of the different parameters of interest. For instance, to study the impact of O 3 on the OLR, these techniques exploit infrared satellite retrievals or model outputs, which can both be subject to strong regional biases and inaccuracies. Independent constraints can, in principle, be obtained from spectrally resolved OLR (i.e. the integrand of broadband OLR, in units of W ⋅ m −2 ⋅ (cm −1 ) −1 ) derived from infrared hyperspectral sounders, which give access to the spectral signature of individual climate processes and feedbacks. While the potential for using spectrally resolved OLR (hereafter abbreviated as SOLR) for the study of climate change has been known for some time (e.g. refs. 10,11 ), the interest for these measurements has only grown relatively recently with the advent of stable multi-year observations 12 . Among the latest achievements, important results were obtained by refs. 13,14 and 15,16 from the measurements of the Atmospheric Infrared Sounder (AIRS, aboard NASA's Aqua satellite 17 ). In particular, they identified a series of biases in climate models that compensated each other partially, and did not show up in comparisons with broadband fluxes. More recently, based on an analysis of 10 years of AIRS global-mean radiances, Pan et al. 18 calculated significant cooling trends in the carbon dioxide (CO 2 ) ν 2 band for channels sensitive to the lower and the mid-stratosphere. Clear differences were highlighted in trends derived from synthetic AIRS radiances computed by considering the output of a free-running general circulation model (GCM) and the European Centre for Medium-Range Weather Forecasts (ECWMF) reanalysis data set (ERA-Interim), confirming the challenges in accurately modeling and assimilating the stratospheric climate.
The infrared atmospheric sounding interferometer (IASI 19,20 ) on board the Metop satellites is a hyperspectral sounder like AIRS, which has the advantage of covering the 645-2760 cm −1 thermal infrared region without any gaps, which makes it well suited for detecting changes in the Earth's SOLR related to trace gases variations. Because of their different overpass time (9.30 versus 13.30 a.m. and p.m. local time for IASI and AIRS, respectively), the joint exploitation of the measurements from the different infrared instruments (also including the Cross-track Infrared Sounder (CrIS) on Suomi NPP which was launched in 2011 21 but whose measurements were never used to date to derived SOLR) could provide meaningful information on, for example, the diurnal cycle of the SOLR or for the confirmation of the variations observed in the SOLR. Since its launch in 2006, IASI has shown a very good stability over time and its measurements provide a good fundamental climate data record 22,23 . Among others, IASI is used as a reference for the inter-calibration of infrared sensors by the Global Space-Based Inter-Calibration System 24,25 . Brindley et al. 26 were the first to exploit the IASI measurements to study interannual variability in the radiance based on 5 years of observations. Variations were found to be small enough to allow the identification of robust changes in regions affected by feedback processes. Right after, Bantges et al. 27 focused on the detection of changes in the Earth's OLR spectrum by comparing radiance spectra of the interferometric monitor for greenhouse gases (IMG), the interferometric infrared spectrometer (IRIS) and the IASI sounder. Among others, clear signatures were identified in the spectra related to increases in well-mixed greenhouse gases. Here, we analyze the trends in 10 years (2008-2017) of IASI clearsky SOLR retrieved from a dedicated algorithm recently developed by Whitburn et al. 28 . Considering the relatively short time period, this study can be seen as a proof of concept in which one of the primary goals is to evaluate the sensitivity of the IASI-derived SOLR product to detect and quantify small changes in the SOLR and to show how these can be linked to changes in surface and atmospheric conditions. In particular, we address the following questions: • Can statistically significant SOLR trends be detected from IASI channels that are sensitive at different altitudes in the troposphere and in the stratosphere? What is the magnitude of the changes?
• Can these statistically significant trends be related to known changes in climate processes?
• Can the effect of increasing greenhouse gases be detected on the SOLR in the 10 years of IASI measurements?
The data and the method for the trends estimation are detailed in section 4. Figure 1a shows the zonal SOLR trends at a resolution of 0.25 cm −1 . For the sake of clarity, these are expressed as a percentage of the average SOLR in the year 2008 (percentage per year). Indeed, as the radiance decreases with the wavenumber in the range [645 2300] cm −1 , trends presented in absolute values become lower for increasing wavenumber. Note that the choice of the reference year matters very little for the calculation of trends, as the relative changes in SOLR are very small. Stippling indicate trends not significantly different from zero. An upper bound was set at 75°North and South of the Equator since we observed that the cloud filter used to identify clear-sky scenes in the SOLR retrieval algorithm becomes less accurate at high latitudes and also because the surface occupied by seas is lower at the poles. The area-weighted mean of the zonal LT for the tropical region (24°S-24°N) and the Northern and Southern high latitudes (60°-75°N and S) are shown as well (Fig. 1b, d, respectively). These allow assessing the trends on a more global scale and to identify small changes in the SOLR that might not appear at the zonal scale. The associated confidence limit (shaded area) were calculated in the same way, i.e. averaged from the 95% lower and upper limit of the zonal LT. From a radiation budget perspective, To help the interpretation of the results, we also evaluated at what altitudes the largest changes in radiation occur (in absorption or emission) by calculating the relative change (in percentage) of the Earth's radiation after passage through 1-kmthick layers of the atmosphere: (L i+1 − L i )/(L i ) × 100, with L i the Earth's radiance (W ⋅ m −2 ⋅ (cm −1 ) −1 ⋅ sr −1 ) at the top of the layer i. The change in each layer (%) was determined for a tropical (Fig.  1c) and a subarctic (Fig. 1e) standard atmosphere 29 from a set of IASI synthetic spectra calculated with the Atmosphit line-by-line radiative transfer code 30 .

Figures
In addition to the zonal trends, we also calculated the global LT on a 2°× 2°grid for a set of selected channels sensitive to different altitudes, as well as global temperature LT derived from the ERA5 reanalysis dataset over the same time-period 31 at different levels of pressure. The calculation of the temperature LT was made from 24-hourly data that were averaged per day. Note that considering temperatures interpolated at 9.30 a.m. instead does not affect the results. These global LT distributions of SOLR and temperature are shown in Fig. 2 and are used to support our analysis of Fig. 1.

Surface and lower troposphere
In the atmospheric window regions (795-970 cm −1 , 1070-1230 cm −1 , and 2090-2170 cm −1 ), where the atmosphere is mostly transparent, the changes in the SOLR reflect principally the changes in the surface temperature (T s , compare Fig. 2a with for the ERA5 T s (Fig. 2g) where typical El Niño/Southern Oscillation (ENSO) and PDO patterns can be clearly identified over the Pacific Ocean (boxes 1 and 3) 35 . Other warming trends that could not be related to a particular climate phenomenon are observed in the tropics along the west coast of the Pacific Ocean and contribute also to the zonal mean positive LT. Note that it can be verified that a change in surface temperature results in a relative change in radiation flux that increases quasi-linearly with wavenumber as illustrated with the straight line drawn on Fig. 1b, d.
In the northern hemisphere between 45°and 60°, LT do not differ significantly from zero in the window channels (Fig. 1a). However, when looking at Fig. 2a, we find an important spatial heterogeneity with both cooling and warming regions. Two large areas that exhibit strong negative trends are identified. One is located in the North Atlantic (about −0.15% per year, extending up to~68°N) and can be related to a particular phenomenon called the North Atlantic Warming Hole (NAWH) (Fig. 2a, box 2), likely due to the melting of the Greenland ice sheet and affecting the Atlantic meridional overturning circulation (AMOC) in the same region 36 . The second one (left part of box 1) is found in the Western Pacific and is part of the PDO pattern 35 . These two cooling areas are counterbalanced by the strong PDO warming pattern in the Eastern Pacific. Above 60°N, the zonal LT (Fig. 1a) become negative (between −0.1 and −0.4% per year, about −5 × 10 −4 to −2.5 × 10 −4 W ⋅ m −2 per year). Up to 68°N, these are driven by the cooling in the North Atlantic because of the relative importance of the region in the total area occupied by seas at these latitudes (~65%). Oddly, above 68°N, the zonal trends remain negative while Fig. 2a shows a dominance of positive LT. However, the zonal LT at these latitudes are not trustworthy (uncertainty close to 100%) due to a low number of observations (partly because of a lower accuracy of the cloud filter) and a higher seasonality and local variations of the SOLR.
In the Southern Hemisphere (below 45°), different warming and cooling areas are also observed. Over the South Atlantic and Pacific (box 4, Fig. 2a), these likely reflect the activity of the  Antarctic Dipole (ADP) which corresponds to a high-latitude climate phenomenon in the air-sea-ice system that is strongly influenced by the ENSO variations 37 . During El Niño events, positive (negative) sea surface temperature anomalies generally occur in the South Pacific (Atlantic). The patterns are similar but opposite during La Niña events. Because of the two La Niña events in the beginning and the strong El Niño in the end of the studied period, the dipole exhibits here positive LT in the South Pacific and negative LT in the South Atlantic. Strong positive and negative LT are also observed close to the Antarctic shores in the Indian Ocean (between 0°and 90°E) and in the Western Pacific Ocean (between 100°and 150°E), respectively (box 5, Fig. 2a). The reason for these cooling and warming trends is unclear and reflects a variation in the air-sea-ice system, possibly linked to the activity of the ADP as well. As for the Northern Hemisphere, the positive and negative LT compensate each other between 50°and 60°S, resulting in zonal LT not significantly different from zero (Fig. 1a). Again, surprisingly, the zonal LT become negative below 60°S while positive LT dominate at these latitudes on the global distribution (Fig. 2a) and should therefore be considered with care. The strong regional disparities in the LT evidenced here for both the Northern and the Southern Hemisphere that do not show up in the zonal mean highlight well the limits of the approach if it is not completed by an analysis based on global distributions. With the follow-up of the IASI mission, longer time series will become available in the near future and allow to better evaluate and characterize the climate processes at high latitudes in the Northern and the Southern Hemisphere.
When analyzing the LT in the window regions on Fig. 1a for each channel individually, we find a fine spectral structure (LT slightly less positive or negative than on average) corresponding to low tropospheric CO 2 and H 2 O absorption lines which reflect the changes in their atmospheric concentrations. This becomes more evident on the area-weighted average trends (Fig. 1b, d). For CO 2 , the increasing concentrations in the atmosphere increase the amount of Earth's radiation absorbed and thus reduce the OLR 10 . This is confirmed when looking at the global LT distribution for the strong 2093.5 cm −1 CO 2 channel (Fig. 2b) mostly sensitive to the first layers of the atmosphere in the tropics and at mid-latitudes (Fig. 2b, left panel). In contrast to the temperature changes close to the surface (Fig. 2g), the LT are mainly negative, except over the Northern Pacific Ocean, where a positive PDO pattern remains (Fig. 2b, orange box). For the H 2 O window channels, the change in SOLR is closely linked with the changes of the T s and atmospheric temperature (T a ) as a consequence of the Clausius-Clapeyron process 38 . In the tropics (at high latitudes), the increase (decrease) of T s leads to an increase (decrease) of the specific humidity and, therefore, to an increase (decrease) of the absorbed radiations (Fig. 1b, d) 10,38 . In addition to CO 2

Mid-, upper troposphere and stratosphere
For the channels peaking in the mid-and upper troposphere, and in the stratosphere (Fig. 1c, e), corresponding mainly to the regions of absorption of H 2 O, CO 2 , CH 4 , and O 3 , the LT result from both the increase of greenhouse gases concentrations and the atmospheric temperature changes. The latter are determined by dynamical processes that originate mostly from changes in SST in the tropics. A rise of SST (mainly related to ENSO cycle here) is known to enhance the deep convection and to cause a strengthening of the Brewer-Dobson circulation (BDC) 40 . In the tropics, this translates to a warming of the troposphere but a cooling of the lower stratosphere [41][42][43][44] . This cooling is also reinforced by the impact of the increase in CO 2 concentrations 45 . Conversely, at high latitudes (mostly in the Southern Hemisphere), the strengthening of the BDC induces a stronger downward motion which leads to a warming of the stratosphere through adiabatic heating 40,43,44,46 . This warming is also reinforced by the effect of the recovery of the ozone hole 47 . These effects are observed in the sequence of Fig. 2h-l, where a progressive inversion of the sign of the LT can be seen on the atmospheric temperature distributions when moving from the troposphere to the stratosphere (starting at 17-18 km of altitude in the tropics and at 8-9 km in the Arctic regions).
We first focus on the changes in the tropical and mid-latitude regions. Immediately obvious on Fig. 1a, b are the systematic negative trends observed in the ν 2 (below 795 cm −1 ) and the ν 3 (after 2200 cm −1 ) CO 2 absorption bands and, to a lesser extent, in the ν 4 CH 4 band (around 1306 cm −1 ) which range from about −0.05 to −0.3% per year. Similar negative patterns were also reported by 48 in the LT derived from AIRS radiance measurements (expressed in brightness temperature) and by 12,38 based on model simulations. At the center of the ν 2 (below 700 cm −1 ) and ν 3 (after 2250 cm −1 ) CO 2 bands (and for the 791.75 cm −1 line), all the radiation emitted by the surface is absorbed by the atmospheric CO 2 (with the largest percentual changes occurring close to the tropopause around the 80-90 hPa levels (16-17 km, Fig. 1c)). The SOLR observed at the top of the atmosphere corresponds therefore only to the emission of radiations from above in the stratosphere. In these spectral ranges, the increase in CO 2 concentrations would thus result in increased SOLR emissions which translates into a positive LT. However, the cooling trend observed indicates that the effect of the stratospheric temperature decrease (Fig. 2l) dominates here. Conversely, in the wings of the CO 2 bands and in the ν 4 CH 4 band, characterized by lower maximum sensitivity levels (~350-500 hPa), the negative LT in the tropics and at Northern hemispheric mid-latitudes (between 45°N and 30°S) reflect the increase in the CO 2 and CH 4 IR absorption which shifts the effective emitting level of radiations to higher (and therefore substantially colder) altitudes of the troposphere 38,48,49 . The effect of the tropospheric warming is thus here negated by the CO 2 and CH 4 increase. This is evidenced on the LT distribution for the tropospheric CO 2 channel at 740.75 cm −1 (Fig.  2c) where homogeneous negative changes are observed between 45°N and S, while the temperature LT at the corresponding altitude of sensitivity (Fig. 2i, j) show a more important variability. This indicates that the effect of the CO 2 increase is the dominant factor driving the LT. For N 2 O, as it was also the case for the ν 1 band in the window region (around 2224 cm −1 ), the effect of increasing concentrations is again not visible in the LT in the ν 3 N 2 O band (around 1230-1320 cm −1 ) due to the overlapping atmospheric variations of CH 4 .
In the ν 2 H 2 O band (around 1595 cm −1 ), the LT also reflect the changes in the state of the atmosphere at different altitudes. For the mid-tropospheric channels (~500 hPa, Fig. 1c), a pronounced latitudinal structure is found on Fig. 1a between 40°N and S with trends being either positive or negative. The global distribution of the LT at 1441.25 cm −1 (Fig. 2d) reveals a strong regional pattern. The latter likely reflects the changes in the main convergence (upward motions) and subsidence (sinking motions) zones associated with large-scale atmospheric circulation which carries heat and moisture 50 . Convergence zones are typically associated with a strong local increase of water vapor which increases the amount of radiation absorbed. Its effect on the SOLR (known as super greenhouse effect (SGE) 14,50 ) exceeds that of the warming of the mid-troposphere and results in a decrease of the SOLR in total.
Conversely, subsidence zones are characterized by a drying of the air masses which translates into a decrease of the absorption (an increase of the SOLR) in the ν 2 H 2 O band. In particular, this effect is illustrated by the large crescent-shaped positive LT observed in the tropical Western Pacific (Fig. 2d, box 6) which corresponds to a well-known region affected by severe droughts during ENSO events 51,52 . Similar patterns were already observed for the longwave cloud feedbacks 53 , supporting our analysis on the convective processes.
In contrast to the mid-tropospheric channels, the very strong H 2 O absorption channels that have a maximum sensitivity higher in the troposphere (typically~150-250 hPa, Fig. 1c) show systematically negative LT in the tropics and at mid-latitudes (Fig. 1a). For those channels, as it was also the case for the stratospheric CO 2 channels, the atmosphere is almost opaque to the surface emissions and a significant fraction of the SOLR emerges from the lower stratosphere. Here, the SOLR decrease are likely the result of two combined effects 38 : (1) the cooling of the stratosphere in link with the strengthening of the BDC and the CO 2 increase, and (2) the increase of H 2 O concentrations in the upper troposphere that goes along with the rise in tropospheric temperature as a consequence of the Clausius-Clapeyron process. This is confirmed by the LT distribution of the strong H 2 O absorption line at 1507 cm −1 (Fig. 2e) which shows negative trends over most of the tropical and mid-latitude regions. Two regions (boxes in Fig. 2e), though, exhibits a significant positive trend. The first one is located in the Northern Pacific between 30 and 45°N and is likely influenced by the warming of the lower stratosphere in the area (Fig. 2k). The second one, found in the Southern Pacific, is not correlated with a temperature increase in the stratosphere and might instead reflect the drying of air masses in the upper troposphere, possibly related to a local subsidence zone.
The ν 3 O 3 band (around 1042 cm −1 ) is characterized by strong positive zonal LT comprised between +0.05 and +0.12% per year (+0.8 to +1.5 × 10 −4 W ⋅ m −2 per year), with the strongest changes occurring below 30°N and S (Fig. 1a, b). Between 30°a nd 45°N and S, close to the center of the band, LT generally do not differ significantly from zero. When looking at the global distribution of the LT for the strong O 3 channel at 1057.25 cm −1 (Fig. 2f), we find an important spatial heterogeneity. Between 30°N and S, most of the positive LT are observed over the Pacific Ocean. Over the Atlantic, LT are mainly negative while they are close to zero (or not significantly different from zero) above the Indian Ocean. Conversely, the [30°-45°N and S] regions show large negative patterns in the Pacific, which are compensated in the zonal mean by the positive LT elsewhere. In general, these LT in the tropical and mid-latitude regions do not correspond to the changes observed in the atmospheric temperature at~50 hPa (Fig. 2l), where the O 3 band is mostly sensitive, and this therefore suggests a variation in the O 3 concentrations. Over the Pacific Ocean (between 45°N and S), the observed LT can be explained by the ENSO cycles. O 3 is mainly produced in the tropics in the lower and mid stratosphere and its abundance is directly modulated by the strength of the upwelling motions in the tropical lower stratosphere. During an El Niño event (in particular the 2015-2016 event here), the strengthening of the upwelling of the BDC in the tropics brings more air generally poor in O 3 from the troposphere into the stratosphere, leading to a decrease in O 3 concentrations in the lower stratosphere (and thus a positive LT for the period considered here) 54,55 . Conversely, over the mid-latitude Pacific (between 30°and 45°N and S), the enhanced downward motions of the BDC lead to an increase of the O 3 concentrations which translates into negative LT 55,56 .
At high latitudes, the interpretation of the LT is easier because most of the changes can be explained by a change in temperature at different altitudes. In the center of the ν 2 and the ν 3 CO 2 band, and for the strong ν 4 CH 4 and ν 2 H 2 O absorption lines (Fig. 1a, d), almost opaque to the surface emissions, the LT are strongly positive (up to 0.2-0.4% per year) above 40°S and 55°N and reflect the warming of the stratosphere caused by a strengthening of the BDC. In the wings of the CO 2 bands, the LT become progressively negative as the atmosphere is more transparent to the surface at these wavenumbers and also because of the CO 2 increase. For H 2 O, apart from the strong lines in the center of the ν 2 band, the LT appear mostly not significantly different from zero. Finally, in the ν 3 O 3 band, the LT observed probably reflect the opposite effects of the stratospheric warming and the O 3 increase, the latter being likely due to both the strengthening of the BDC 55,56 and the stratospheric ozone recovery at 20-25 km of altitude 44,57,58 . This is especially evident in the center of the band, where both impacts almost cancel each other out (LT close to zero) (Fig. 1d). When moving away from the center, the LT become more negative because of the larger fraction of the signal coming from the surface.

DISCUSSION
Making use of 10 years (2008-2017) of IASI measurements, we have estimated the linear trends (LT) in zonally averaged spectrally resolved OLR (SOLR) at a spectral sampling of 0.25 cm −1 (6621 channels in total). The LT result from both the atmospheric composition and temperature changes. Despite the limited time period and the difficulty of disentangling the effect of the different parameters affecting the SOLR, clear spectral signatures of the long-term changes in greenhouse gases concentrations (in particular CO 2 , CH 4 but also CFC-11 and CFC-12) could be identified unambiguously, especially in channels sensitive to the mid-and upper troposphere. For CO 2 and CH 4 , the increase in the atmospheric concentrations resulted in an increase of the fraction of SOLR absorbed which translated into a negative LT of about −0.05 to −0.3% per year in the ν 2 and ν 3 CO 2 and in the ν 4 CH 4 absorption band. For the LT mainly reflecting the changes in atmospheric temperature, most of them may have been connected to changes in the atmospheric circulation originating from an increase of the SST driven by the ENSO and PDO activity. In particular, the LT for the ν 2 H 2 O mid-tropospheric sensitive channels revealed a strong regional pattern reflecting the changes in the main convergence and subsidence zones. In the three atmospheric window regions, the LT at high latitudes (above 60°N and S) were associated with high uncertainties because of the low number of observations and possible issues with the cloud filter. All these results show the potential of the IASI-derived spectrally resolved OLR in improving our understanding of the current Earth's climate and of its long-term changes. In particular, such high spectral samplingnever reached before with any satellitederived OLRis of great interest for the identification of errors and biases of opposite sign in the output of climate models that cancel each other in total OLR. With the launch of IASI on Metop-B in 2012 and on Metop-C in 2018, which L1C product is consistent with the reprocessed IASI-A radiance dataset 22 , and the future launch of IASI-NG (New Generation, with a spectral sampling of 0.125 cm −1 ) on the Metop-SG suite of satellites, the time series of spectrally resolved OLR will be extended by at least another 25 years of measurements.

METHOD Data
The clear-sky SOLR are derived from a complete reprocessed dataset of IASI/Metop-A radiance measurements (with the latest version of the L1C) 22 in the range 645-2300 cm −1 at the 0.25 cm −1 native spectral sampling of the L1C spectra. The algorithm for the conversion of the spectra to the fluxes is detailed in Whitburn et al. 28 . The SOLR are derived with a good precision generally within ±0.005 W m −2 (cm −1 ) −1 depending on the wavenumber. The 2300 cm −1 upper bound was chosen to avoid solar reflectance contamination. The cloud-free scenes were selected based on the cloud information from the Advanced Very High Resolution Radiometer (AVHRR) flying on board the MetOp satellites as well. Only IASI observations associated with a 0% cloud coverage were considered. The cloud flag is rather conservative, keeping only~14% of all observations. The distribution of the fraction of clear-sky scenes is shown in Supplementary Fig. 2. The raw dataset consists of 10 years of daily global SOLR (2008-2017) averaged on a 2°× 2°grid, separately for day and night observations. For our analysis, we focused on daytime measurements. We restricted the analysis to observations over sea to ensure capturing trends related to global climate phenomena. Indeed, measurements over land are subject to stronger local variations in space and time (especially for the channels that are significantly influenced by the lower layers of the atmosphere and the surface) because of the smaller heat capacity compared to seas and the larger heterogeneity of the surface parameters (e.g. orography and land cover). These variations make it more difficult to identify changes in the SOLR on larger scales, in particular when focusing on relatively short time periods, such as those considered here.

Trends estimation
For each IASI channel, linear trends (LT) were computed from the daily SOLR, averaged zonally by 2°latitude, using the approach described in Gardiner et al. 59 . The method relies on fitting a low order Fourier series (here n = 3) capturing the intra-annual variability in the SOLR plus a linear term (the annual trend) to the dataset. The confidence limits associated with the trend estimates are determined using a bootstrap resampling. A LT is considered significantly different from zero if the 95% confidence interval does not contain zero. Even though land measurements were filtered out, their influence on the calculated zonal trends remains, especially near the coasts, as changes that occur over land can propagate over oceans. In addition to this, another effect that likely influences the zonal trends is the heterogeneity in the distribution of the cloud-free observations (see Supplementary Fig. 2) which modulates the weight of each region when averaging the SOLR by band of latitude. The impact of the zonal variability in the SOLR and of the heterogeneity in the cloud distribution are reflected in the uncertainties in the zonal linear trend. Nevertheless, these limitations have to be kept in mind for the analysis of the results.

DATA AVAILABILITY
Monthly means of IASI-derived spectrally resolved OLR 60 are freely available for all users through the IASI-FT website: https://iasi-ft.eu/data-access/OLR/. Daily IASI-derived spectrally resolved OLR can be made available upon request to the corresponding author.