Interactive comment on “ Metrology of the Solar Spectral Irradiance at the Top Of Atmosphere in the Near Infrared Measured at Mauna Loa Observatory : The PYR-ILIOS campaign ” by Nuno

The manuscript describes measurements of the near infrared solar irradiance from the surface and its extrapolation to top-of-the-atmosphere (TOA) solar irradiance using the Langley-plot technique. The experiment was performed at the high-altitude Observatory in Mauna-Loa, Hawaii and is a follow-up experiment to a similar one performed at Izaña, Tenerife, Canary Islands. The resulting TOA irradiances are compared to space-based measurements of solar irradiance, and a rudimentary uncertainty analysis is provided. The paper provides an important dataset to try to resolve observed discrepancies in space-based solar irradiances in the near infrared range. The advan-


Introduction
An accurate knowledge of solar spectral irradiance (SSI) remains central to the study of the climate on Earth.The variability in the ultraviolet (UV) part of the spectrum and its influence on climate via the mechanisms of solar-terrestrial interactions, simulated by chemistry-climate models (Gray et al., 2010;Ermolli et al., 2013), constitutes most of the research in SSI measurements.Despite its extremely low variability, < 0.05 % over a solar cycle (Lean, 1991;Harder et al., 2009), the near-infrared (NIR) part of the spectrum plays a major role in the Earth's radiative budget due to its quasi-total absorption by water vapour (Collins et al., 2006).The determination of its absolute level remains challenging (Meftah et al., 2017): the measurement of the topof-atmosphere (TOA) SSI started nearly 50 years ago and evolved both with ground-based and space-borne instruments, and a consensus on the absolute level in the NIR part is still to be achieved (Bolsée et al., 2014;Hilbig et al., 2018).
Aircraft-borne instrumentation at an altitude of 12 km provided the first TOA SSI measurement dataset in 1969 (Arvesen et al., 1969) with an on-board standard of spectral irradiance.
-At Mauna Loa Observatory (MLO), Shaw (1982) conducted a campaign with a 10-channel (UV, visible and NIR) filter radiometer and Gröbner and Kerr (2001) Published by Copernicus Publications on behalf of the European Geosciences Union.
N. Pereira et al.: Metrology of solar spectral irradiance with a double Brewer spectrophotometer measuring in the range 300-355 nm.Kindel et al. (2001) provided TOA SSI in the range 350 to 2500 nm, measured with a spectroradiometer.All of these measurement campaigns used different types of 1000 W lamps, traceable to National Institute of Standards and Technology (NIST) standards as calibration sources.
Finally, the CAVIAR (Menang et al., 2013) and CAVIAR2 (Elsey et al., 2017) spectra were obtained with an infrared Fourier spectrometer (FTIR) calibrated against National Physical Laboratory (NPL) standards at the UK Met Office observation site in Camborne, in the range 1-2.5 µm.
TOA SSI values from all the above-mentioned groundbased campaigns were obtained using the Langley plot technique that permits extrapolation to the TOA irradiance in atmospheric windows chosen according to criteria detailed in Sect.2.2.The monitoring of the absolute spectral calibrations is secured through comparisons with relative stable secondary standards.The reliability of the traceability to primary irradiance standards is an advantage for groundbased measurement.Performing these measurements based on world reference sites for the determination of TOA physical quantities, such as IZO and MLO, on days with often pristine conditions, ensures a high accuracy of the TOA extrapolations (Shaw, 1975(Shaw, , 1976;;Kiedron and Michalsky, 2016;Toledano et al., 2018).
On the other hand, space-borne SSI measurements covering the NIR range started in the 1990s, though these were limited to wavelengths shorter than 2.4 µm.From the SOL-SPEC instrument family, the instrument SOSP (SOlar SPectrum) on board EURECA (Thuillier et al., 1981) that pioneered the space-borne NIR absolute solar spectroscopy released the ATLAS3 reference spectrum (Thuillier et al., 2003).An upgraded version of the SOLSPEC instrument, SOLAR/SOLSPEC, including a fully refurbished NIR channel, readout electronics and extended wavelength range up to 3 µm of SOLSPEC, flew from 2008 to 2017 on board the International Space Station (ISS) (Thuillier et al., 2009), releasing the SOLAR2 (Bolsée, 2012;Thuillier et al., 2014) and SOLAR-ISS(IR) (Meftah et al., 2017); SOLSPEC is the space-borne instrument that measured SSI farther in the NIR.The instrument providing the longest time series of SSI measurements in the NIR is the SIM (Spectral Irradiance Monitor) prism spectrometer on SORCE (Solar Radiation and Climate Experiment) launched in 2003 (Harder et al., 2000(Harder et al., , 2005) ) and still on orbit but with infrequent operational time, due to the end of battery life.Another instrument contributing to NIR SSI measurements is SCIAMACHY (Scanning Imaging Absorption Spectrometer for Atmospheric Chartography) (Noël et al., 1998;Burrows et al., 1995), a remote sensing spectrometer adapted to measure SSI.The latest data release is SCIAMACHY V9 (Hilbig et al., 2018).
All above-mentioned NIR datasets reasonably agree up to 1.3 µm.When comparing SORCE and ATLAS3, the differ-ence between both does not exceed 2 % in the NIR range, which is a consequence of SORCE being scaled up to AT-LAS3, due to incompatibilities of fractional TSI (total solar irradiance) between both datasets (Harder et al., 2010).
At 1.6 µm, corresponding to the minimum opacity value of the solar photosphere, differences up to 8 % (reaching 10 % at 2.2 µm) were observed between ATLAS3 and SO-LAR2 (Thuillier et al., 2014).This bias motivated the development of new ground-based instrumentation measuring the SSI NIR: CAVIAR and IRSPERAD (Bolsée et al., 2014;Menang et al., 2013).The data of both experiments confirmed this bias, both showing a level closer to that of SO-LAR2.Posteriorly, SOLSPEC and SCIAMACHY data reprocessing processes tend to intermediate values between ATLAS3 and SOLAR2 (Meftah et al., 2017;Hilbig et al., 2018).
In this paper we present a rerun of the IRSPERAD experiment, named PYR-ILIOS, carried out in July 2016.While still using the Langley plot technique and calibration against the PTB black body, this new experiment differs from IRSPERAD in three aspects: first, the observation site is MLO instead of IZO; second, possible sources of systematic uncertainties have been identified and fixed (see Sect. 2.1); third, the traceability of the calibration to the primary standard was improved (see Sect. 2.6).A detailed estimation of the uncertainty budget will be presented in Sect.3, followed by the presentation of the obtained spectrum and its comparison with space-borne and ground-based spectra described in this section, along with a discussion on the status of the NIR SSI measurement.

Instrumentation
The core of the direct Sun measurement instrumentation is a Bentham NIR spectrometer: it consists of a double monochromator placed inside a thermally stabilized container, with light detection by a PbS cell.An optical fiber guides the sunlight between the entrance slit of the spectrometer and the diffusor of a 7.2 • field-ofview (FOV) sunlight-collecting optics (telescope).The telescope is connected to an EKO Sun tracker that provides a tracking accuracy (https://eko-eu.com/products/solar-energy/sun-trackers/str-22g-sun-trackers, last access: 12 December 2018) of 0.01 • .The working wavelength range is from 0.6 to 2.3 µm, with a nominal 10 nm bandpass.The instrument characteristics are given in depth in Bolsée et al. (2014) and have remained unchanged since.No modifications have been made either to the telescope or to the spectrometer.Nevertheless, a factory defect in the assemblage of the components was detected and rectified: the lens focusing the light collected in the optic fiber into the spectrometer entrance slit was properly fixed into its barrel support for the PYR-ILIOS campaign, which was not previously the case for the IRSPERAD campaign at IZO.Another change relative to the IRSPERAD campaign was that the thermally stabilized spectrometer container was placed indoors in a thermally stabilized environment, which reduced thermal stress due to outdoor exposure and improved the stability of the spectrometer's response.

Langley plot method
The wavelength-dependent direct transmitted solar irradiance in the atmosphere is described by the Beer-Bouguer-Lambert (BBL) law.For spectral regions where molecular absorption is negligible and only Rayleigh and aerosol scattering are present, the BBL law is written in the following form: where E 0 is the irradiance at the top of the atmosphere (TOA), m is the air mass factor (AMF) as a function of the solar zenith angle (SZA) θ and τ is the optical depth that depends on λ.D is the ratio between the Earth-Sun distance at the moment of the measurement and the mean Earth-Sun distance; subscripts R and A stand for Rayleigh and aerosol, respectively.Because the aerosol vertical profile over the measurement site at the moment of the measurement is unknown, aerosol AMF is approximated to Rayleigh AMF (Schmid and Wehrli, 1995); considering m A ≈ m R ≈ m, defining τ = τ R +τ A and taking the logarithm of Eq. ( 1), it can be rewritten as Provided that τ (λ) remains constant for a series of measurements of E 0 (λ) taken over a given range of m (θ ) (spreading over a half day), the TOA value of E 0 (λ) is thus the intercept at the origin (m = 0) of the least-squares regression to the data series E (λ) as a function of m(θ ).

Atmospheric windows
The wavelength domains for which the Langley plot method described in Sect.2.2 is valid, i.e. atmospheric windows, were determined through the model using a procedure developed in Kindel et al. (2001) and also used in Bolsée et al. (2014): using a TOA reference spectrum as input, the MOD-TRAN (MODerate resolution atmospheric TRANsmission) (Berk et al., 2014) RTM (radiative transfer model) was used to simulate irradiances measured at the ground, as a function of the measurement site parameters, for a series of AMFs.The Langley plot method was applied to these simulated irradiances, and the wavelengths for which the synthetic E 0 recreated the input TOA within 0.5 % were kept as valuable wavelengths for the Langley plot; these set of wavelengths were grouped in contiguous windows called atmospheric windows.

Absolute calibration
The absolute calibration was performed against a primary standard of spectral irradiance, the BB3200pg black body of the PTB.It has been extensively described in Sapritsky et al. (1997) and Sperfeld et al. (1998Sperfeld et al. ( , 2000)).The spectral irradiance equation describing the black body emission is calculated using Planck's law: where ε BB and A BB stand, respectively, for the effective emissivity and the aperture of the BB3200pg, d BB for the distance between the black body aperture and the optic centre of the telescope and n for the refractive index of air; c 1 and c 2 are the first and second radiation constants.The fundamental parameter, the temperature of the cavity T BB , is known with a standard uncertainty of 0.5 K (∼ 0.02 % for a nominal temperature of 3000 K) with a drift lower than 0.5 Kh − 1 (Friedrich et al., 1995;Werner et al., 2000;Taubert et al., 2003).The uncertainties on ε BB and A BB are 1×10 −4 (0.01 %) and 0.04 mm (0.03 %), respectively (Woolliams et al., 2006).The distance between the black body aperture and the telescope optical active surface, the diffuser d BB , is the sum of two distances: d BB = d s + d T , where d S is the distance between the black body and the first optical surface of the telescope, the quartz plate, and d T is the distance between the quartz plate and the diffuser.The uncertainties on d s and d T are 0.05 mm (Woolliams et al., 2006) and 0.5 mm, respectively; the combined uncertainty on d BB is 0.5 mm, 0.04 % at the nominal distance of 1384.05 mm.
The absolute calibration coefficient R, that converts the spectrometer signal into irradiance, is given by Eq. ( 4): with S BB being the signal recorded by the spectrometer and E BB , the emission of the black body, given by Eq. (3).During the calibration campaign at PTB, two different temperature set points, 3016.5 and 2847.6 K, were used to build the response curve, R BB .The distance d BB was kept fixed at 1384.05 mm so that the black body aperture was seen by the entrance optics with an angular extension of 0.5 • .
N. Pereira et al.: Metrology of solar spectral irradiance

Radiometric characterization
The spectrometer was characterized at the laboratory of the Belgian Institute for Space Aeronomy (BIRA-IASB) for the uncertainty on the measured signal, the detector sensitivity to temperature and for the wavelength scale.The flat field of the detector was measured during the ground-based campaign at MLO and the linearity was verified during the calibrations at the PTB laboratory: -The flat field of the entrance optics was measured during the ground-based campaign.The telescope was angularly displaced from the normal Sun direction thanks to an angular fine-tuning mechanism, for a series of angular positions for two orthogonal directions.The agreement between both directions' data curves allows an insensitivity of the signal to solar depointing better than 0.05 • to be estimated, although a finer angular sampling would be necessary to accurately determine the angular limits of this insensitivity.Given the 0.01 • pointing accuracy of the Sun tracker, the response of the instrument is considered to be insensitive to pointing during the campaign.
-The temperature sensitivity of the spectrometer was determined in the laboratory (Bolsée et al., 2014).During the campaign, the spectrometer box was placed indoors with its temperature being constant within 0.1 • C, equivalent to the resolution of the temperature probe readout; no temperature correction on the signal was thus applied.
-For the verification of the linearity of the detector, the telescope was placed at several different distances from a stable 200 W lamp.The measured signal as a function of distance was successfully fitted to an inverse square law function, demonstrating the detector linearity within a 2-decade dynamic range.

Relative calibration
A set of six FEL lamps (F102, F104, F417, F418, F545, F546) were used as relative calibration standards, to monitor a possible change of response of the spectrometer during the measurement campaign.Taking as reference the lamps' signal measured at the PTB (27 April 2016), S PTB Fj (λ), four additional relative calibrations were performed: -Immediately before the start of the measurement campaign on 29 June (i = 1), the signal of the six lamps, S MLO1 Fj , was measured on site.This first MLO relative calibration was valuable to monitor the spectrometers' response change between the calibration at PTB and the beginning of the field measurements.During this 2-month period that included the transportation of the equipment, a decrease of response varying between 1 % and 3 % in the 1000 to 2200 nm range was detected.
-During the 20-day measurement campaign, three relative calibrations were performed: on 7 July (i = 2), 14 July (i = 3) and 19 July (i = 4).The cumulated loss of response between 29 June and July varied from 1.5 % to 0.5 % in the 800 nm to 1.8 µm domain.
The corresponding correction factor for each relative calibration is where N stands for the total number of lamps, j for the lamp number and i for the calibration day index.K(λ) was obtained by linear interpolation for all days of the campaign.

Ground-based campaign
The for the measurement of the concentration of atmospheric carbon dioxide.It is considered a world reference site to accurately determine extraterrestrial constants via the Langley plot method (Shaw, 1975(Shaw, , 1976;;Kiedron and Michalsky, 2016).

Data selection and analysis
From the 20-day campaign, 12 high-quality half-days, all during morning time, were kept for analysis.The selection criteria were verification of cloudless clear skies and a Langley plot correlation coefficient R 2 > 0.9.The morning data of the days 2, 3, 5, 7, 8, 9, 10, 11, 13, 14, 16 and 17 July 2016 were kept for analysis; a subset of these selected Langley plots is shown in Fig. 1, for four different wavelengths.
3 Uncertainty budget

Uncertainty on the spectrometer signal
The raw uncertainty of a spectrometer measured signal, S raw x , regardless of its source, either solar (S S ), black body (S BB ) or lamp S PTB Fj , S MLO Fj signal, is a function of the intrinsic noise of the measured physical signal convolved by the spectrometer's transmission and detector's response.The uncertainty on a measured signal u(S raw x ) was determined in the laboratory by calculating the standard deviation for a sample of measured signals at several intensities from a 1000 W stable lamp (Bolsée et al., 2014).This uncertainty is shown in Fig. A1.
Additionally, all measured signals, S raw x (λ), are affected by an uncertainty term due to the finite bandpass of the instrument, u(C λ ), and the uncertainty on the determination of the true wavelength scale, u(C λ ) (Obaton et al., 2007).
where δ(λ) stands for the maximum deviation in the determination of the real wavelength scale of the spectrometer.δ(λ) was determined in the laboratory by measuring the deviation between the measured and the corresponding nominal peak values of a series of well-known emission rays of Xe, Ar and Kr lamps as well as of lasers and Pen-Ray lamps, δ(λ) < 0.2 nm for the working wavelength range.BW stands for the spectrometer bandpass of 10.63 nm, measured in the 3.2 Langley plot sensitivity to aerosol daily variation

Uncertainty on a calibrated direct Sun measurement
The expression for a calibrated solar measurement, E(λ) is with S S (λ), R(λ) and K(λ) being expressed by Eqs. ( 6), ( 4) and ( 5), respectively.The uncertainties associated with the factors in Eq. ( 7) were calculated using the law of propagation of uncertainties (LPU) and are represented in Fig. 2. The similarity of shapes of the curves of the individual uncertainties reflects the convolution of the measured signals by the spectrometer's response.The largest contribution to the calibrated solar signal comes from the uncertainty on the absolute calibration which is dominated by the uncertainty u(S BB ) of the measured signal, S BB , of the black body, whereas the uncertainty on the emission of the black body, u(E BB ), is known within 0.2 % for the totality of the wavelength range.

Uncertainty on the determination of the TOA irradiance
The uncertainty in the determination of the TOA irradiance via the Langley plot method, u(E LP 0 ), corresponds to the unwww.atmos-meas-tech.net/11/6605/2018/Table 1.List of relative uncertainties terms expressed as percentages.The coverage factor is k = 1 for all terms.A and B stand respectively for type A and type B uncertainties, while C stands for combined uncertainty according to GUM (2008).u(C λ ) and u(C λ ) are calculated for a solar signal.The prefix u, for uncertainty, is omitted for each term of the first row, for the sake of clarity. AMF 0.04 870 0.11 1.29 0.18 0.62 0.52 0.43 1.29 1.46 0.66 0.06 0.67 4 0.19 1020 0.02 0.04 0.01 0.04 0.10 0.33 0.04 0.07 0.14 0.14 0.34 0.41 0.17 0.11 0.20 8 0.79 1640 0.09 0.17 0.06 0.06 0.13 0.08 0.19 0.26 0.12 0.06 0.14 certainty on the determination of the intercept at origin, P 0 , when applying a linear regression on Eq. ( 2).The uncertainty on the measured E in the Langley plot method logarithmic space, u(log(E)), and the uncertainty in the u(E LP 0 ) irradiance value are given by where E LP 0 = exp(P 0 ) gives the irradiance TOA value.The uncertainty in P 0 was estimated using two independent methods.
-A Monte Carlo method was employed.Given a measured Langley plot dataset consisting of (m i , log ) is created, where each log(E * i ) is affected by a random normal distributed quantity, with a standard uncertainty given by Eq. ( 8), and each m i is affected by an uncertainty defined in Sect.3.6.The standard deviation in the distribution of the N >>1 retrieved P 0 values corresponds to u(P 0 ), with u(E LP 0 ) given by Eq. ( 9).
-The weighted total least-squares algorithm developed by Krystek and Anton (2007) was used.It computes the uncertainty in the determination for both linear regression parameters using the uncertainties on the measured quantities as inputs, i.e. the uncertainties on E (Sect.3.3) and AMF (Sect.3.6).The uncertainty on the determination of the TOA irradiance, u(E 0 ), matches perfectly for both methods; it is below 1 % for the central wavelength range of 0.9 to 2.2 µm.Figure 2 shows the contribution of all the uncertainty terms detailed in Sect.3. In Table 1 a list of the uncertainty types and values at key wavelengths is presented.

Quantification of the circumsolar radiation
An ideal sunlight-collecting optic device should ideally have an acceptance angle equal to that of the solar disk seen on Earth, ∼ 0.5 • .In practice the FOV is much larger than 0.5 • such that Sun-and sky-scattered radiation enters the FOV of the sunlight-collecting optics, affecting the direct normal Sun measurement.Circumsolar radiation is strongly dependent on aerosols' size and their abundance, increasing with AMF and decreasing with wavelength due to Rayleigh scattering (Blanc et al., 2014).The estimation of circumsolar radiation was done with the aid of the LibRadtran (Mayer and Kylling, 2005) RTM.LibRadtran computes the radiance field of the Sun-and sky-scattered radiation.The integral of this radiance field over the solid angle of the acceptance cone of the entrance optics is the amount of circumsolar irradiance (CSI) measured by the spectrometer in excess of the normal direct Sun irradiance (DNI) (Gueymard, 2001).For standard clear-sky atmospheric conditions observed at MLO and for typical aerosol charges values measured during the mission, the quantification of CSI is shown in Fig. A2.Given the un-certainty budget, the impact of the circumsolar radiation can be considered negligible.

Estimation of air mass factors' uncertainty
As referred to in Sect.2.2, the absence of knowledge of the vertical profile of the relevant species, namely aerosols, is a limiting factor for accurately calculating the AMF.The uncertainty in the AMF calculation is based on the approach of Schmid and Wehrli (1995), who considered that m A , due to the presence of stratospheric aerosols, could take the form , to be used as input for the determination of the Langley plot parameters' uncertainty (Sect.3.4).
The possible bias introduced at the Langley plot's intercept at origin by a realistic non-constant aerosol concentration during the measurement was estimated considering a measured aerosol optical depth (AOD) profile.For a given measured Langley plot consisting of (m i , log(E i )) and regression parameters E LP 0 and τ , a synthetic Langley plot ) stands for the real diurnal aerosol optical depth profile measured with AERONET (available at λ = 870, 1020, 1640 nm) and τ R (λ, t i ) the Rayleigh optical depth calculated according N. Pereira et al.: Metrology of solar spectral irradiance to Bodhaine et al. (1999).This bias at the intercept at origin, expressed as a ratio, E LP 0 /E * 0 , averaged over the selected days is −0.2 %, +0.4 % and +0.1 % for 870, 1020 and 1640 nm, respectively.The signal of the bias replicates the signal of the AOD morning trend measured at MLO, and the larger negative bias at 1020 nm relative to 1640 nm is due to the more pronounced AOD negative trend at 1640 nm.Assuming that the interval | E LP 0 −E * 0 | comprises the true value of the intercept at origin, E 0 , within a rectangular distribution, the corresponding uncertainty u(E amounts to 0.06 % at 870 and 1640 nm and 0.1 % at 1020 nm, which is added quadratically to the uncertainty on E LP 0 (Sect.3.4) to determine the uncertainty on E 0 .u(E AOD 0 ) is interpolated linearly to the working wavelength range.

Results
The PYR-ILIOS TOA SSI results are obtained by averaging the E λ (0) obtained by the Langley plot method for the 12 half-days that satisfied the data selection criteria detailed in Sect.2.8.PYR-ILIOS and other space-borne and groundbased instruments' datasets described in the Introduction are compared to the SOLAR-ISS(IR) from Meftah et al. (2017) in Fig. 3.
The mismatch between the PYR-ILIOS and IRSPERAD dataset varies between 2 % and 4.5 % in the central wavelength range between 1.0 and 1.8 µm, attaining 5 % in the 2.1 µm window and peaking to a maximum of 6 % in the 1.5 and 2.2 µm windows.Except for the shorter wavelengths' (λ < 900 nm) region, uncertainties do not explain the observed mismatch between both.The higher disagreement is observed in the far end of the spectrum, with discrepancies of up to 13 % between CAVIAR2 and ATLAS3 and SORCE.Below 1.3 µm all the datasets are compatible within the uncertainties bars.

Discussion
The difference observed between IRSPERAD and PYR-ILIOS is not explained by the uncertainties of both datasets.An atmospheric bias is not considered because MLO and IZO are world reference sites for the determination of extraterrestrial constants (Shaw, 1976;Kiedron and Michalsky, 2016;Toledano et al., 2018) and the atmospheric perturbations in ground-based SSI measurements are negligible (Elsey et al., 2017;Bolsée et al., 2016;Weber, 2015).By carrying out the new PYR-ILIOS experiment, we unveiled a defect of fixation of the focusing lens.Due to the fact that the instrument was moved between the IRSPERAD precampaign relative calibration (31 May 2011) and the start of the Sun measurement campaign (1 June 2011 onwards), the effect of the lens' eventual movement was not considered and therefore not monitored; this defect likely biased the SSI obtained during the IRSPERAD campaign in a nonreproducible way.This defect was detected and corrected for the PYR-ILIOS campaign and the relative calibration strategy adapted to identify possible similar issues: the instrument was installed and powered on and the lamps were measured; the solar measurements began immediately afterwards, without displacing or powering off the instrument.The PYR-ILIOS relative calibration procedure highlights the importance of monitoring ground-based pre-campaign instruments' response with secondary standards.Additionally it justifies the choice of PYR-ILIOS as a more reliable measurement than IRSPERAD, due to the higher confidence in the traceability of the instrument's calibration to the black body primary standard.
In the higher disagreement region around 1.6 µm, the most recent data versions of SOLAR/SOLSPEC and SCIA-MACHY instruments, SOLSPEC-ISS and SCIAMACHY V9, respectively, as well as PYR-ILIOS converge to an intermediate level between SOLAR2 and ATLAS3.This convergence is also observed for longer wavelengths: in the 2 µm region PYR-ILIOS and Kindel et al. are in reasonable agreement, while the level of the two SCIAMACHY V9 adjacent bands (1.9-2.05 and 2.2-2.4 µm) suggests that it is also in agreement with the two ground-based datasets; on the other hand, in this region, both data versions of the SOL-SPEC/SOLAR still retain the 8 % difference to ATLAS3 and SORCE.
A rerun of the measurement campaign at IZO would be crucial to understand the observed discrepancy between PYR-ILIOS and IRSPERAD datasets.Data from the SORCE successor, TSIS, which has been on board ISS since December 2017, are expected to further increase the understanding of SSI in the NIR.

Figure 1 .
Figure 1.Measured irradiance and respective Langley plot fits for the four AERONET wavelengths, 870, 1020 and 1640 and 2065 nm, shown for the morning data of 2, 9, 10 and 13 July 2016.

Figure 2 .
Figure 2. Individual uncertainties contributing to the combined uncertainty in the TOA SSI, u(E 0 ).Black-body-associated quantities (u(E BB ), u(S BB ) and u(R)) and lamp-associated quantities (u(K)) are plotted for the full wavelength working range, while solarmeasurement-associated quantities (u(S s ), u(E) and u(E 0 )) are plotted in the atmospheric windows' wavelengths.

Figure 3 .
Figure 3. Ratio of ground-based and space-borne spectra relative to SOLSPEC-ISS(IR).Uncertainty at ±1σ is represented by the shaded areas.
1 and m O 3 standing for the ozone air mass.Assuming a rectangular distribution of m A delimited by k 1 = 1 and a k 1 = 0.2, the standard deviation of m A can be calculated as u PYR-ILIOS campaign took place during the first 20 days of July 2016 at the Mauna Loa Observatory (MLO) on the island of Hawaii.The MLO (19.53 • N, 155.58 • W) is situated at 3397 m a.s.l.; it is the leading long-term atmospheric monitoring facility on Earth, a primary calibration site for the AErosol Robotic NETwork (AERONET; https: //aeronet.gsfc.nasa.gov/,last access: 12 December 2018), a global station for the Global Atmosphere Watch (GAW) of the World Meteorological Organization (WMO) and the premier site (https://www.esrl.noaa.gov/gmd/obop/mlo/programs/esrl/co2/co2.html,last access: 12 December 2018)