A new airborne broadband radiometer system and an efﬁcient method to correct dynamic thermal offsets

. The instrumentation of the High Altitude and Long Range (HALO) research aircraft is extended by the new Broadband AirCrAft RaDiometer Instrumentation (BAC-ARDI) to quantify the radiative energy budget. Two sets of pyranometers and pyrgeometers are mounted to measure up-ward and downward solar (0.3–3 µm) and thermal–infrared (3–100 µm) irradiances. The radiometers are installed in a passively ventilated fairing to reduce the effects of the dynamic environment, e.g., fast changes in altitude and temperature. The remaining thermal effects range up to 20 W m − 2 for the pyranometers and 10 W m − 2 for the pyrgeometers. Using data collected by BACARDI during a night ﬂight, it is demonstrated that the dynamic components of the offsets can be parameterized by the rate of change of the radiometer sensor temperatures, providing a greatly simplifying correction of the dynamic thermal effects. The parameterization provides a linear correction function (200–500 W m − 2 K − 1 s) that depends on the radiometer type and the mounting position of the radiometer on HALO. Furthermore, BACARDI measurements from the EUREC 4 A (Elucidating the Role of Clouds—Circulation Coupling in Climate) ﬁeld campaign are analyzed to characterize the performance of the radiome-ters and to evaluate all corrections applied in the data processing. Vertical proﬁles of irradiance measurements up to 10 km altitude show that the thermal offset correction limits the bias due to temperature changes to values below 10 W m − 2 . Measurements with BACARDI during horizontal, circular ﬂight patterns in cloud-free conditions demonstrate that the common geometric attitude correction of the solar downward irradiance provides reliable measurements in this typical ﬂight section of EUREC 4 A, even without active stabilization of the radiometer.


Introduction
Measurements of solar and thermal-infrared irradiance are important to quantify the radiative impact of atmospheric components and surface properties on the Earth's radiative energy budget and to quantify their relevance for climate change. Ground-based observations of the broadband upward and downward irradiances are routinely performed within the Baseline Surface Radiation Network (BSRN) at locations distributed over the entire globe (Driemel et al., 2018). These observations were used in a variety of studies, e.g., characterizing the climatology of cloud radiative effects by Shupe and Intrieri (2004). However, BSRN observations are limited to fixed land locations representing a local environment, e.g., surface albedo or temperature regime. Observations over ocean are obtained only from a few research ships and buoys (Kalisch and Macke, 2012;Colbo and Weller, 2009). Instead, airborne or spaceborne observations resolve the spatial distribution of the radiative energy budget, which is strongly affected by the heterogeneity of the surface albedo, surface temperature, and clouds (Stapf et al., 2020(Stapf et al., , 2021. While satellite estimates of the irradiances at the top of the atmosphere require radiative transfer simulations, airborne observations provide direct measurements of the upward and downward as well as solar and thermal-infrared irradiance. Furthermore, radiative processes such as cloud-top cooling or aerosol layer warming need to be quantified to understand the influence of radiative processes on atmospheric dynamics (e.g., Wendisch et al., 2008;Simpfendoerfer et al., 2019). These quantities are derived from profiles of net (downward minus upward) irradiances, which can be measured directly only by airborne observations (Bucholtz et al., 2010) or from balloon and helicopter platforms (Egerer et al., 2019;Siebert et al., 2021).
Broadband irradiances F are measured by radiometers, in particular pyranometers (solar, 0.3-3 µm) and pyrgeometers (thermal-infrared, 3-100 µm). The measurement principle of most common radiometers, as discussed here, is based on thermopile sensors. Some radiometers use photo-diode sensors, which are sensitive only to a limited spectral range, while thermopile sensors in general detect the entire spectral range of electromagnetic radiation. To define the wavelength selectivity of a thermopile radiometer and to protect the sensor from environmental impacts, the sensor is capped by a dome. Special materials, e.g., quartz glass, or silicon, and filter coatings guarantee a relatively constant sensitivity of the instrument over the desired spectral range (Gröbner and Los, 2007). However, an unshaded pyrgeometer may suffer from leakage effects when solar radiation is transmitted above the cut-on wavelength of the pyrgeometer dome interference filter (Marty, 2000;Meloni et al., 2012). The overall performance of broadband radiometers is determined by the radiometric calibration accuracy, dome spectral transmissivity, angular response, direct solar heating, dome temperature effects, and long-term measurement stability (e.g., Philipona et al., 1995Philipona et al., , 2001Wendisch and Brenguier, 2013;Gröbner et al., 2014).
The combination of a thermopile sensor and an optical filter dome can affect the thermal equilibrium of the entire instrument and thus bias the measurements, especially when operating the radiometer on aircraft, where fast radiation and temperature changes may occur (e.g., Albrecht et al., 1974;Saunders et al., 1992). Thermopile radiometers typically have an inertia of a few seconds to changes in radiation. However, Curry and Herman (1985) and Foot (1986) showed that thermal equilibrium, especially during rapid ascents and descents, may be reached only after several minutes. This effect is caused by a differential change in the temperatures between the dome and the thermopile sensor, which was estimated by Philipona et al. (1995) at up to ±1.5 • C. To minimize these uncertainties, Philipona et al. (1995) suggested adding two additional temperature sensors in the dome and parameterizing the irradiance bias. However, commercially available broadband radiometers, which are built for groundbased operation, do not include these temperature sensors and require careful post-processing (Ehrlich and Wendisch, 2015).
Airborne measurements, especially of the downward solar irradiance, are also affected by the aircraft attitude when the radiometers are fixed to the aircraft fuselage (Wendisch et al., 2001). By definition, the atmospheric irradiance refers to a horizontally aligned surface, which is not maintained by the radiometer sensor during pitch and roll aircraft movements. Depending on the solar zenith angle, Wendisch et al. (2001) calculated that a misalignment of ±1 • already results in an offset of up to 3 % in the downward solar irradiance for a solar zenith angle of 60 • . Actively stabilized pyranometers, such as those proposed by Wendisch et al. (2001) and Bucholtz et al. (2008), can minimize such uncertainties, but these techniques are complex, expensive, and not applicable to all aircraft installations. A post-correction as suggested by, e.g., Bannehr and Schwiesow (1993) and Boers et al. (1998) is limited to the direct solar component of the incoming radiation (cloud-free conditions) and depends on the accuracy of the estimation of the fraction of direct solar radiation, the characterization of the pyranometer mounting, and the measurement of the aircraft attitude.
This attitude correction requires synchronized pyranometer and aircraft attitude measurements, which may not be given due to the slow response of the broadband radiometer (Freese and Kottmeier, 1998). As shown by Ehrlich and Wendisch (2015), characterizing the radiometer time response and reconstructing the measurement time series significantly improve the performance of airborne radiometers and help to analyze the radiation field in complex cloud and surface conditions (Egerer et al., 2019;Stapf et al., 2020).
Given these known issues, airborne measurements of broadband solar and thermal-infrared irradiance are delicate and require a proper setup of the radiometers on the aircraft as well as careful post-processing aimed at correcting positional and thermal biases. Here, a new radiometer package, the Broadband AirCrAft RaDiometer Instrumentation (BACARDI) installed on the High Altitude And Long Range (HALO) research aircraft operated by the German Aerospace Center (Deutsches Luft und Raumfahrtzentrum, DLR), is introduced. BACARDI is comprised of a set of two Kipp & Zonen pyranometers (CMP22) and pyrgeometers (CGR4) that are mounted in a fixed position to the aircraft fuselage. The housing and mounting is constructed to minimize thermal effects. However, thermal offsets remain, and therefore a novel approach to correct for them is developed. To illustrate the basis of the correction, Sect. 2 gives a review of the theory of the broadband radiometer radiative budget. The radiometer characteristics, data acquisition specification, and the instrument design, including the aircraft specific instrument mounting and shielding, are described in Sect. 3. All basic corrections including the radiometric calibration, the reconstruction of the time response, and the attitude correction of the solar downward irradiance are specified in Sect. 4. Based on a dedicated calibration Figure 1. Simplified radiative energy budget of a broadband radiometer composed of (1) the transmitted solar irradiance, (2) the transmitted thermal-infrared irradiance, (3) the irradiance emitted by the dome, and (4) the irradiance emitted by the sensor. Blue arrows indicate solar radiation, and red arrows are used for thermalinfrared radiation.
flight, a novel approach to correct the dynamic thermal offset of the radiometer in rapidly changing temperature conditions is developed. Section 5 outlines the correction approach and the application to measured vertical profiles of solar and thermal-infrared irradiance. The overall performance of BACARDI is tested using measurements during the EUREC 4 A (Elucidating the Role of Clouds-Circulation Coupling in Climate) field campaign (Bony et al., 2017;Stevens et al., 2021). In Sect. 6 measurements of heating rate profiles and the consistency of measurements during circular flight pattern are analyzed by comparison to radiative transfer simulations. The key benefits of the new system are summarized in Sect. 7.
2 The radiative energy budget of broadband radiometers

Basics
Broadband radiometers, which are based on thermopile sensors, use the temperature increase in the illuminated receiver area compared to a shaded reference area as the primary measure. A simplified drawing and the radiative budget of a broadband radiometer are illustrated in Fig. 1. Based on the Seebeck effect, which describes the thermopile sensitivity α (unit: V K −1 ), the resulting temperature difference T between the sensor surface temperature T s and the reference area temperature T ref generates a detectable voltage U th , which is used to compute the irradiance: In most cases, the temperature of the reference area is measured with a standard temperature sensor, e.g., Pt-100 or ther-mistor. The sensor itself absorbs the incoming irradiance F in and also emits radiation F out in the thermal-infrared wavelength range. In the case of pyrgeometers, the emission is a major issue. From a simplified energy budget of the thermopile sensor, the total effect of the net irradiance F net,dyn on the sensor can be described by with C the heat capacity of the sensor surface (unit: J m −2 K −1 ) and K the thermopile thermal conductance (unit: W m −2 K −1 ).

Thermal equilibrium
Assuming thermal equilibrium, ∂T s /∂t = 0, Eq. (2) reduces to the net irradiance F net,stat in static conditions: Following the simplified radiative budget of a broadband radiometer as shown in Fig. 1, the standard formulation of the calibration equation of broadband radiometers in radiative equilibrium is postulated (e.g., Fairall et al., 1998;Ji and Tsay, 2000). The incoming and outgoing irradiances F in and F out at the sensor surface can be expressed independently by with τ d,sol the solar transmissivity of the dome, T d the temperature of the dome, and τ d,ir and d,ir the thermal-infrared transmissivity and emissivity of the dome. The sensor is characterized by the sensor emissivity s,ir and temperature T s . σ is the Stefan-Boltzmann constant. The sensor reflects the incoming irradiance F in with the reflectivity ρ s , while the dome can reflect the outgoing irradiance F out with a reflectivity ρ d . For simplicity, all transmissivities, emissivities, and reflectivities are broadband quantities; spectral dependencies of the materials are not considered. Assuming that ρ s · ρ d 1 and using the following assumption for the polynomial of fourth order to replace T s in Eq. (5) by T ref , the equations finally can be resolved for the thermal-infrared and solar incident irradiance:

1566
A. Ehrlich et al.: A new airborne broadband radiometer system The last term in Eqs. (7) and (8), the so-called longwave and shortwave leakage, is only a function of the incident irradiance and the ratio between solar and thermal-infrared transmissivity of the dome, which are determined by material properties of the dome and/or filter coating. There is evidence for the existence of such errors due to spectral imperfections of the dome, and possible corrections were developed (e.g., Pascal and Josey, 2000). By careful selection of the dome material and coating (low τ d,sol and high τ d,ir for pyrgeometers and high τ d,sol and low τ d,ir for pyranometers), this error can be minimized or even neglected as confirmed by long-term comparison of different shaded and illuminated pyrgeometers (Meloni et al., 2012;Kipp & Zonen, 2014). For pyranometers, the longwave leakage is correlated with the net thermal-infrared irradiance measured by a pyrgeometer (third term in Eq. 8) and can thus hardly be distinguished from the thermal dome offset (second term in Eq. 8).
Neglecting this leakage effect, Eqs. (7) and (8) can be reduced to the commonly known formulas (Philipona et al., 1995). For the thermal-infrared irradiance measured by pyrgeometer it is with the parameters A 1 , A 2 , and A 3 summarizing the instrument characteristics. If the temperature dependence of the thermopile sensitivity, A 2 · σ · T 3 ref , is compensated for electronically within the radiometer, the first term of Eq. (9) can be further reduced, leading to the formulation by Albrecht et al. (1974): with a ir the adjusted pyrgeometer thermopile sensitivity (unit: V (W m −2 ) −1 ) and b ir the pyrgeometer dome factor A 3 . The last term, also known as the window heating offset, corrects for a thermal imbalance between the dome and sensor surface mainly caused by solar radiative heating of the dome in static conditions. The dome factor b ir in Eq. (10) theoretically defines the ratio of thermal-infrared emissivity d,ir to transmissivity of the dome τ d,ir . Ji and Tsay (2000) showed that the dome factor, experimentally determined from a black-body calibration of the instrument, yields significantly higher values than expected from theory. Only by using data obtained in thermal equilibrium is the theory fulfilled. This indicates that the commonly used higher dome factor implies non-equilibrium effects. Optimizing the thermal design of the radiometer can reduce the window heating offset such that a dome temperature measurement can be omitted and no dome factor is needed (Meloni et al., 2012;Gröbner et al., 2014;Kipp & Zonen, 2014). Applying a similar transformation, the solar irradiance measured by pyranometers in Eq. (8) reduces to The adjusted pyranometer thermopile sensitivity a sol (unit: V (W m −2 ) −1 ) includes the weak temperature dependence of the thermopile as defined in theory by Eq. (8), which can often be compensated for by the construction of the radiometer or determined in extended laboratory calibrations. The static pyranometer thermal dome effect is scaled with the dome factor b sol and the temperature difference between the dome and sensor. This effect is often called the zero or dark offset since it is mainly caused by radiative cooling of the dome and is best visualized as a negative offset during night measurements in the absence of solar irradiance. A second dome with high thermal conductivity, e.g., quartz, in good thermal contact with the instrument housing can reduce this error to a few watts per square meter (W m −2 ) (Philipona, 2002;Reda et al., 2005;Kipp & Zonen, 2016). Ventilation of the dome can further reduce the zero offset (Michalsky et al., 2017). If the thermal dome effect cannot be neglected, available correction methods are applied. These have been developed based on either an additional dome temperature measurement or the simultaneous measurement of the net thermal-infrared irradiance by a pyrgeometer (e.g., Bush et al., 2000;Haeffelin et al., 2001;Dutton et al., 2001).

Dynamic environment -no thermal equilibrium
The assumption of thermal equilibrium is valid for standard ground-based measurements with slowly varying environ-mental conditions. However, if the radiometers are subject to fast temperature changes like during airborne measurements, e.g., during ascents and descents, the slow adjustment of the sensor temperature (first term in Eq. 2) needs to be considered. An offset voltage U s and offset irradiance F s will be generated by the thermal lag between the reference and the sensor, which is initiated by the thermal conductance and capacity of the sensor.
Replacing the sensor temperature T s in Eq.
(2) by the reference temperature, T s = T ref + T , the thermal reaction of the sensor to an outside temperature change can be described by The assumption that T changes much less than T ref leads to the dynamic sensor thermal offset F s defined as the difference between the net irradiance in static F net,stat and dynamic conditions F net,dyn : This error correction term for dynamic temperature changes is proportional to the time derivative of the reference temperature. F s often is called "zero offset B" or "zero offset due to temperature change" and is mostly specified by the instrument manufacturer for a fixed temperature change of 5 K h −1 . However, during airborne observation, especially during ascents and descents, faster temperature changes of the order of Kelvin per minute (K min −1 ) occur. A similar behavior is expected for the dome, leading to a slow adjustment of the dome temperature T d . Due to the different thermal properties of the dome and sensor, the dynamic thermal offsets in both parts do not compensate. The dynamic dome effect F d can then be expressed as As indicated by Bush et al. (2000), the temperature difference of the dome T d depends linearly on the temporal change in T ref : with the coefficient γ (unit: s) characterizing the relationship. Assuming T d T ref and approximating the fourthorder polynomial similar to Eq. (6), the dome effect reduces to Adding both effects (Eqs. 13 and 16), the total dynamic thermal offset F is described by Based on the initial assumption that T d and T s can be related linearly to T ref , Eq. (17) indicates that the thermal offset can be linearly parameterized by the change rate of the sensor reference temperature, providing the dynamic thermal offset correction coefficient β (unit: W m −2 K −1 s). During data post-processing, such a parameterization can be applied to correct the irradiance measurements in high dynamic conditions.
3 Design of BACARDI for operation on HALO

Broadband radiometers
For the measurements of upward and downward broadband irradiance, F ↓ and F ↑ , separated into the solar and thermal-infrared spectral range, BACARDI combines two sets of Kipp & Zonen pyranometers (CMP22) and pyrgeometers (CGR4). The CMP22 pyranometers detect radiation in the wavelength range of 0.2-3.6 µm, which covers almost the entire solar spectral range (Kipp & Zonen, 2016). The CGR4 pyrgeometer is sensitive to wavelengths between 4.5 and 42 µm, covering a large fraction of thermalinfrared radiation (Kipp & Zonen, 2014). Both radiometers use thermopile sensors, providing a sensitivity in the range of 10 µV (W m −2 ) −1 . The radiometric calibration of the radiometers, which refers to the entire solar and thermalinfrared spectral range is repeated regularly by the manufacturer a few months in advance of a HALO measurement campaign. The radiometers are calibrated as a secondary standard (Class A) through comparison with a reference instrument traceable to the World Radiation Center. For the pyranometers, this comparison is done in the laboratory, and for the pyrgeometers, the comparison is performed outside under mainly cloud-free conditions during nighttime. Additionally, for both radiometers the temperature dependence of the thermopile sensitivity is determined within a climate chamber for the temperature range of −40 to 50 • C. Calibration uncertainties typically range below 1 % for the CMP22 pyranometers and 4 % for the CGR4 pyrgeometers (2 times the standard deviation confidence level). The temperature dependence of the thermopiles does not exceed 0.5 % for a wide temperature range (−30 to 50 • C). To track the sensor temperature, each radiometer is equipped by the manufacturer with a platinum (Pt-100) resistance thermometer. The respective quartz and silicon domes function as wavelength bandpass filters and are characterized by a cosine response, which is less than 1 % off from theory over the entire 180 • field of view (Kipp & Zonen, 2014. The optimized thermal design of both radiometers reduces the window heating offset to less than 4 W m −2 and makes them suited for aircraft operation. However, the time response of the radiometers needs to be considered for airborne measurements. As specified by the manufacturer, the CMP22 pyranometer typically reacts quicker with a 1/e (63 % adjust-ment) response time of about 2 s, while the CGR4 pyrgeometer is characterized by a response time on the order of 6 s.

Electronics and data acquisition
The CMP22 and CGR4 radiometers do not contain any internal signal conditioning and only provide a low-voltage (in the range of 10 µV (W m −2 ) −1 ) thermopile signal and a fourwire Pt-100 temperature signal.
To mitigate the effects of electromagnetic noise, the wiring of the low-voltage signal is as short as possible and a signal conditioning unit is used. This unit is placed on the radiometer mounting plate inside the fuselage, where it is protected by an electromagnetic compatibility shielding metal box. The signal conditioning is based on isolated Dataforth 8B modules that are plugged into a backplane. The 8B30-02 module with an input range of ±50 mV is used for amplification of the thermopile signals up to ±5 V, with an accuracy of 0.05 %. The four-wire Pt-100 resistance is translated into a voltage (0-5 V) by the 8B35-01 module, which covers the temperature range ±100 • C with an uncertainty of ±0.2 • C.
The output voltage signals have a bandwidth of 3 Hz and are recorded by the HALO basic data acquisition system BAHAMAS (BAsic HAlo Measurement And Sensor system; Krautstrunk and Giez, 2012) with a 10 Hz data rate and 18bit resolution. A signal path calibration is performed after aircraft installation of the radiometers, which includes all wiring, connectors, electronics, and the data acquisition. For the calibration, the radiometers are replaced by either a highprecision constant voltage source (Burster 4463) to simulate the thermopile output or a high-precision resistance decade (Burster 1427) to simulate the Pt-100. Both calibration references are set to values covering the operating range of the radiometer and Pt-100 thermometer by a computer-controlled calibration routine. The calibration factors are implemented in the first post-processing step of the BACARDI raw data.

Mounting on HALO and fairing
The integration of BACARDI on HALO uses the standard 10 in. × 7 in. fuselage apertures. Because HALO is equipped with four upper and six lower central apertures, some flexibility in installation depending on the layout of the actual scientific instrumentation is given. A drawing and visualization of one BACARDI sensor package is shown in Fig. 2. The mounting plates, to which the radiometers are attached, compensate for the mean pitch angle of the HALO aircraft, which amounts to about −3 • in normal flight conditions. To reduce the cable length between the radiometers and the electronics (amplifier and Pt-100 conditioner), the electronics housing is attached to the mounting plate on the opposite side of the radiometers inside the fuselage.
The radiometers of BACARDI are in an aerodynamic fairing to minimize the environmental influence on the radiation measurement by, e.g., ice aggregation or water droplet impact and heating by solar radiation. To minimize aerodynamically induced temperature gradients across the instrument, a passive ventilation of the fairing is implemented to keep the instruments close to thermal equilibrium with its surrounding environment. The ventilation is designed to divert the main airflow containing droplets or particles around the radiometer housings. The fairing exhaust acts as a water drain and avoids entrapment of water inside the fairing. Thus, the design of the fairing for the upward-looking radiometers slightly differs from that for the downward-looking radiometers. Figure 3 shows measurements of all four sensor temperatures compared to the ambient temperature measured on HALO. In general, the sensor temperatures are higher, especially in cold conditions, due to low heat transfer in the rather low-density air and the heat conduction from the cabin. Temperature adjustments to changes in ambient temperature (change in altitude) significantly lag in time and may lead to thermal offsets as discussed in Sect. 2.3. This is most prominently indicated by the hysteresis between ascent (upper branch) and descent (lower branch) in Fig. 3. However, comparing only the sensor temperatures, the differences are larger between the pyranometers and pyrgeometers than between the upper and lower setup. This indicates that temperature adjustments are rather a matter of the internal sensor housing of CGR4 and CMP22, their internal heat transfer, and the mounting order (CGR4 mounted in front of CMP22). The ventilation within the fairing is similar in the upper and lower sensor package.
To enable maintenance work, e.g., changing desiccant cartridges and signal calibration, easy access to the radiometers is considered necessary. Therefore, the upward-looking radiometers can be detached from inside the cabin without removing the fairing, whereas for the downward-looking radiometers, it is sufficient to dismount the fairing. The calibration of the pyranometers and pyrgeometers provided by the manufacturer includes tracking changes in the thermopile sensitivity with changing instrument temperatures. For the CMP22 pyranometers, the change in the sensitivity is estimated in the range of ±0.3 % for the temperature range between −20 and 50 • C. Significantly lower sensitivities of up to −2 % are registered when temperatures reach −40 • C. For the CGR4 pyrgeometers, lower differences in the sensitivity are reported. Here, deviations do not exceed ±0.5 %, with the largest positive biases observed for the lowest and highest temperatures and slight negative offsets in between. Figure 4a shows the sensor temperatures of all four radiometers and the ambient static air temperature for the  (1) with ventilation inlet (2) and exhaust (3), pyrgeometer (4), pyranometer (5), desiccant cartridge (6), and electronic box (7). The red arrow indicates the flight direction. EUREC 4 A flight on 22 January 2020 (flight ID HALO-0122). Except for takeoff and landing, the flight altitude and, thus, the ambient temperature changed only for one flight section when HALO climbed to 13.5 km altitude. At cruising altitude, minimum ambient temperatures down to −60 • C were observed. However, due to thermal conduction from the aircraft, the sensor temperatures remained significantly higher and did not reach −40 • C, which is the lower boundary of the calibration certificate. The same holds for all other flights during EUREC 4 A. In other environments, e.g., Arctic conditions, in which low temperatures are reached at lower altitudes with higher air densities, extending the calibration to lower temperatures needs to be considered as demonstrated by Su et al. (2008). The effect of the temperature dependence on the sensor sensitivities is shown in Fig. 4b and c. The changes in the sensor temperature are well documented, and the radiometric calibration is adjusted by up to 1 % (0.1 V (W m −2 ) −1 ). Converted into irradiance, this corresponds to a maximum correction of 5 W m −2 for the downward solar irradiance during local solar noon (16:00 UTC) when the Sun is high. Fluctuations in the time series are caused by the presence of clouds, which reduce or enhance the irradiance and thus the cor-rected bias. Due to the lower thermal-infrared irradiances, the differences here are 1 order of magnitude lower.

Correction of sensor response time
The response times τ r of the CMP22 and CGR4 radiometers provided by the manufacturer are evaluated by measurements during a test flight in cloud-free conditions. The response time of the upward-looking pyranometers is determined by the cross-correlation between the measured downward irradiance and the aircraft attitude angles, assuming that the aircraft attitude is recorded instantaneously by the GPS-aided inertial navigation system. A 1/e (63 % adjustment) response time of τ r = 1.2 s, which is slightly lower than reported by the manufacturer, is obtained. The same response time is assumed for the downward-looking pyranometer.
The response times of the CGR4 pyrgeometers are extensively characterized by Ehrlich and Wendisch (2015) in a laboratory study with a reported τ r of about 3 s. EUREC 4 A measurements of flight sections with sharp turns are used to validate the τ r of the BACARDI pyrgeometers. During the turns, the upward-looking radiometer partly observed the warmer lower hemisphere, which caused a sudden increase in the upward irradiance. Based on a detailed analysis of this systematic change, the response time is adjusted to 3.3 s.
The inertia of the measured irradiances caused by these response times are corrected following the deconvolution method proposed by Ehrlich and Wendisch (2015). To minimize numerical effects of the deconvolution at sharp gradients and the sensor noise, a cut-off frequency of 0.6 Hz and a moving average filter with 0.5 s window length are applied to the reconstruction of the pyranometer measurements. For the pyrgeometers, a slightly lower cut-off frequency of 0.5 Hz and a longer window length of 2 s are chosen.

Attitude correction of downward solar irradiance
BACARDI is fixed to the aircraft fuselage and does not actively align with the horizontal plane. Therefore, the measurements are affected by the aircraft attitude. Except for turns, changes in the roll and pitch angles of HALO typically do not exceed ±1 • , limiting the alignment error (Wendisch et al., 2001). Changes in the pitch angle are mostly related to flight altitude and true air speed and are up to ±3 • . For the downward solar irradiance, a post-correction following the approach by Bannehr and Schwiesow (1993) and Boers et al. (1998) is applied. This correction is valid only for the downward direct solar irradiance. Therefore, the relative fractions of direct and diffuse solar radiation in cloudfree conditions are estimated using radiative transfer simulations. The simulations are updated continuously based on available in-flight observations of temperature and humidity profiles. The one-dimensional (1D) plane-parallel radiative transfer solver DIScrete ORdinaTe DISORT 2.0 embedded in the library for radiative transfer is applied (libRadtran; Emde et al., 2016;Laszlo et al., 2016). For the conditions during ACLOUD, a 5 % uncertainty of the simulated fraction of direct radiation amounts to less than 1 % uncertainty of the corrected downward irradiance. In cloudy conditions with 100 % diffuse radiation no correction can be applied. Therefore, final BACARDI data include both uncorrected F ↓ sol to be used for cloudy conditions and a corrected product to be used in cloud-free conditions. A basic cloud mask that is based on a comparison with the expected cloud-free irradiance and the identification of enhanced variability of the downward solar irradiance within a 20 s running window is provided in the published data set.
For the correction, the offset angles of BACARDI, 0 for the roll and 0 for the pitch angle, characterizing the relative alignment of the radiometer with respect to the inertial navigation system of HALO are determined from measurements during test flights. In cloud-free conditions, flight sections in different flight directions are compared to simulations of the theoretical downward solar irradiances. By minimizing the differences between corrected and simulated F ↓ sol , the bestfitting pair of 0 and 0 is derived. For the installation of BACARDI during EUREC 4 A, two test flights are analyzed, one performed before the campaign in the vicinity of Oberpfaffenhofen, Germany, and one during the campaign based in Barbados. In both cases, 0 = +0.3 • and 0 = +2.5 • are obtained, so it can be assumed that the offset angles are stable once BACARDI is installed on HALO. To account for the limitations of the attitude correction, the downward solar irradiance is filtered before publishing the data set. Data are assumed to be valid when the attitude correction factors are less than 25 %. For larger correction factors, roll and pitch angles need to be smaller than 5 • . The excluded data correspond to turns with large roll angles or conditions with low Sun.

Dynamic thermal offset correction
As discussed in Sect. 2.3 and indicated by the sensor temperatures shown in Fig. 3, dynamic thermal offsets need to be considered if the radiometers are exposed to fast temperature changes. F is expected to be proportional to the time derivative of the sensor reference temperature (Eq. 17). To quantify and finally correct this effect for BACARDI as operated on HALO, an exemplary night flight was performed on 15 May 2019. The flight was about 1.5 h long and represents a typical HALO ascent and descent profile including a few level steps before reaching a maximum height of about 13 km. The static air temperature varied between −55 and +20 • C, and the takeoff time was more than 1.5 h after sunset, therefore ensuring that no solar radiation was present.
During the night flight, it is assumed that the solar irradiances measured by the pyranometer are zero: F sol = 0 W m −2 . Thus, deviations from zero are used to quantify the dynamic thermal offset F . According to the theory, F mainly depends on the derivative of the sensor reference temperature ∂T ref /∂t. For pyranometer measurements from the night flight, this relation is shown in Fig. 5. Data measured shortly after start and before landing (gray symbols) were not used to determine the thermal offsets. Since the calculation of ∂T ref /∂t amplifies the measurement noise, the signal is smoothed with a 10 s running mean filter before and after applying the derivative function. Through this treatment, no significant additional noise is added to the thermopile measurement when applying the thermal correction to the raw data. To remove long-term trends of the ambient temperature and instrument performance, the data are additionally detrended with a high-pass filter. This also removes potential static thermal offsets as described in Sect. 2.2. For the pyranometers, two averaging times are applied in the highpass filter: 100 s, displaying only very fast sensor responses, and 1000 s, which also includes slower adjustments of the thermal equilibrium. Both filters result in an almost identical trend that indicates that the pyranometers respond similarly to fast and slow temperature changes.
To quantify the dynamic thermal offset correction, two different fit approaches are selected. A simple linear fit (see Fig 5), which neglects the thermal dome effect, provides the correction coefficient β (unit: W m −2 K −1 s). A more complex multivariable fit (not shown here) including the absolute value of the sensor reference temperature T ref following Eq. (17) is applied but did not show significant improvement as a correction. This result shows that the dynamic dome effect can hardly be discriminated from the dynamic thermal offset of the thermopile, and the simple linear fit sufficiently corrects for F . All correction coefficients derived for the upper and lower pyranometer are listed in Table 1. For detrending the data with the 100 s high-pass filter, the coefficients of the upper pyranometer, β = 235 W m −2 K −1 s, and lower pyranometer, β = 439 W m −2 K −1 s, significantly differ by almost a factor of 2. This indicates that the lower radiometer dome is more strongly exposed to the airflow (slight negative pitch angle of HALO). Assuming similar changes in the internal sensor temperatures (see Sect. 3.3), the lower radiometer is affected by a stronger dynamic thermal offset. Therefore, the coefficients reported here for BACARDI operated on HALO cannot reliably be transferred to other broadband radiometers on other research aircraft. For the thermal-infrared irradiance measured by the pyrgeometers, the assumption of F ir = 0 W m −2 does not apply. On timescales of several minutes, F ir also varies with changing atmospheric conditions and altitude and cannot be assumed to be constant. Therefore, only the detrending with the 100 s high-pass filter is applied. Additionally, only selected flight segments are used to determine F . These sections are characterized by small variations in the thermal-infrared irradiance that match strong variations in temperature. The selected data are shown in Fig. 5b and d for both pyrgeometers. The remaining fluctuations of the pyrgeometers show an excellent correlation with ∂T ref /∂t, which is inverse to the correlation of the pyranometers. The dynamic thermal offset correction coefficient amounts to β = −491 W m −2 K −1 s for the upper pyrgeometer and β = −404 W m −2 K −1 s for the lower pyrgeometer. Compared to the CMP22 pyranometers, β values of both CGR4 pyrgeometers show only a small difference. This might be a consequence of the less exposed domes of the CGR4 compared to the CMP22 in combination with the more efficient ventilation of the CGR4 inside the BACARDI sensor mounting where the CGR4 is placed in front of CMP22 with respect to the flight direction.
Applying these parameterizations, all irradiances measured during the night flight were corrected. In Fig. 6, the upward and downward solar irradiances are compared to the uncorrected measurements for the entire flight. The data are presented as vertical profiles to compare the ascent and descent, which should agree after correction in the absence of any solar radiation. It is obvious that the uncorrected data show a similar pattern of fluctuations for upward and downward irradiance originating from the dynamic thermal effects. The dynamic thermal offset correction reduces this thermal error in both pyranometers from up to 20 W m −2 to values below 10 W m −2 . For the downward irradiance, the best agreement is found using the 100 s high-pass filter, while for the upward irradiance both filter options agree. The remaining bias to F sol = 0 W m −2 is caused by potential static thermal offsets as described in Sect. 2.2 and other uncertainties such as the radiometric calibration of the pyranometer.

Measurement examples
Measurements of BACARDI during the EUREC 4 A field campaign  are used to demonstrate how the applied corrections affect typical analysis of broadband radiation measurements.

Irradiance and heating rate profile
The dynamic thermal offset correction is most relevant when the temperature environment changes rapidly, such as during ascents and descents. Also, the aircraft flight velocity and the air density change the airflow around the sensors and control the adjustment of the thermal equilibrium. Figure 7 shows corrected and uncorrected profiles of all four irradiance components for an ascent up to 10 km altitude measured right after the start of the research flight on 7 February 2020 (flight ID HALO-0207). To interpret this profile, it needs to be considered that during such an ascent, HALO also covers a horizontal distance of about 200 km during which the atmospheric conditions may change. However, to estimate the effect of the dynamic thermal offset correction on the measurements, this case is well suited. Flight sections that do not comply at all with the required conditions, e.g., flight maneuvers of HALO, have been removed from the corrected data.
The ascent is characterized by an apparent cloud layer with cloud top at about 2 km as indicated by the increase in F ↓ sol and the decrease in F ↓ ir at this altitude. Above this cloud, cloud-free conditions above the aircraft prevail. The upward irradiances, solar and thermal-infrared, are both affected by the changing cloud situations below HALO. The general increase in reflected solar radiation above the low-level cloud layer is covered by F ↑ sol , while F ↑ ir drops only for a limited period. Afterwards, the low-level cloud layer likely became thinner along the flight track, resulting in a cloud-top temperature that is more similar to the surface temperature.
This general pattern is shown by both uncorrected and corrected data. Differences, as provided in Fig. 7c and f, increase with altitude and are related to the temperature profile. Between 2 and 3 km altitude, a temperature inversion was present and the thermal offsets became smaller. In general, the uncorrected data underestimate the solar irradiance and overestimate the thermal-infrared irradiance. This is due to the inverse correlation of temperature change and dynamic thermal offset for the CGR4 and CMP22 radiometers as discussed in Sect. 5. While both CGR4 pyrgeometers show an almost synchronized pattern (almost identical β), the dynamic thermal offset correction differs for both pyranometers (higher β for F ↑ sol ). Therefore, the upward solar irradiance is more affected than the downward irradiance. As the dynamic thermal offset is independent of the absolute magnitude of the irradiance, this behavior might also be valid for other conditions with higher surface reflectivity or the presence of more reflective clouds.
Profiles of broadband solar and thermal-infrared irradiance are often used to study and quantify the impact of water vapor, clouds, and aerosol particles on atmospheric heating rates. To calculate atmospheric heating rate profiles, the upward and downward irradiances are combined into the net irradiance F net , independently for both spectral ranges: For the measurement case shown above, the net irradiance profiles for corrected and uncorrected data and their differences are shown in Fig. 8. Differences between corrected and uncorrected data are below 8 W m −2 for the solar irradiance and below 5 W m −2 for the thermal-infrared irradiance. As the upward and downward radiometers are almost equally affected by the temperature change during the ascent, the dynamic thermal correction mostly cancels out for F net . Only the upper and lower pyranometer show slight differ- ences. This implies that the uncorrected irradiances can also be used to estimate F net . The dynamic thermal offset correction only becomes relevant for the profiles of F net,sol at higher altitudes, where temperature changes are quicker. The net thermal-infrared irradiance, F net,ir , significantly differs only for low altitudes below the cloud layer.
Consequently, the atmospheric heating rates, defined as the vertical change in net irradiance, also show only a minor impact of the radiometer dynamic thermal offsets. In Eq. (19) ρ represents the air density and c p the specific heat capacity of the air. From the example profile, heating rates are calculated for a 50 m layer thickness, showing the strongest heating rates of down to −4 K h −1 at the top of the low-level cloud layer. However, the differences between corrected and uncorrected data are less than ±0.2 K h −1 for the entire profile. These results demonstrate that for this specific application of BACARDI, which is based on differences of upper and lower broadband radiometer measurements, a dynamic thermal offset correction could be neglected. This might not be valid if the mounting position of BACARDI on HALO changes for other missions, which deviate from the instrument configurations presented by Stevens et al. (2021).

Impact of solar radiation on downward thermal-infrared irradiance
The HALO flights of EUREC 4 A have mostly been performed at flight altitudes above 10 km and often under cloudfree conditions above HALO, causing high values of downward solar irradiance. In such conditions the solar leakage of the pyrgeometer dome interference filter can produce an overestimation of the thermal-infrared irradiance (Philipona  , 1995;Marty, 2000;Meloni et al., 2012). For cloudfree ground-based measurements, Meloni et al. (2012) identified an overestimation of up to 10 W m −2 depending on the amount of downward solar irradiance. This bias was investigated for BACARDI using radiative transfer simulations of F ↓ ir , which are reliable and can serve as a benchmark for measurement above 10 km and under cloud-free conditions. The simulations have been performed along the HALO track, considering the time of day, the geographical position, and flight altitude of HALO with a temporal resolution of at least 30 s. The radiative transfer solver DISORT 2.0 and the lowtran parameterization of molecular absorption embedded in libRadtran are applied (Emde et al., 2016;Laszlo et al., 2016;Ricchiazzi and Gautier, 1998). In the simulations, the cloud-free atmosphere is defined by merged temperature and humidity profiles from the Barbados Cloud Observatory radiosondes (BCO; Stevens et al., 2016;Stephan et al., 2021) and the frequent dropsonde measurements from HALO (George, 2021).
Filtered for cloud-free conditions, Fig. 9 shows the difference between measured and simulated F ↓ ir as a function of the measured downward solar irradiance F ↓ sol . The data indicate a trend to overestimate F ↓ ir for increasing F ↓ sol . For values of F ↓ sol above 1000 W m −2 , typical for times around solar noon, the bias ranges up to 10 W m −2 , comparable to the findings of Meloni et al. (2012). A linear regression suggests an increase in the bias by 1 W m −2 for each 100 W m −2 increase in F ↓ sol . However, the data show a large variability and the regression suggests a negative bias for the absence of solar radiation. This may be attributed to remaining uncertainties of the radiative transfer simulations and the pyrgeometer sensitivity due to changes in water vapor concentrations above HALO (Nyeki et al., 2017), a permanent bias of the radiometer calibration, and a static thermal offset. Figure 9. Difference between measured and simulated downward thermal-infrared irradiance F ↓ ir for 12 EUREC 4 A flights filtered for cloud-free conditions and flight altitudes above 10 km. The differences are plotted and fitted as a function of the measured downward solar irradiance F ↓ sol .

Solar irradiance during horizontal, circular flight pattern
During EUREC 4 A, HALO frequently flew a circular flight pattern that aimed to quantify the large-scale vertical motion, an eminent parameter characterizing the dynamic state of the atmosphere (Bony et al., 2017). The typical circle had a diameter of roughly 220 km, which corresponds to a permanent roll angle of HALO between 2 • and 3 • . Therefore, the correction of the downward solar irradiance F ↓ sol for horizontal misalignment, as described in Sect. 4.3, becomes more important. At the same time, a circular flight pattern provides observations over the full range of relative solar azimuth angles and is thus an ideal test bed for evaluating the performance of the solar irradiance measurements.
The accuracy of the attitude correction is tested against radiative transfer simulations that are introduced in Sect. 6.2. For the simulations of F ↓ sol the absorption by ozone, which becomes relevant for the typical flight altitude of 10 km, is determined by satellite estimates of the atmospheric ozone column. The sea surface albedo is parameterized on the basis of Cox and Munk (1954) using the 10 m wind speed obtained from the lowermost wind speed value of the HALO dropsondes. For the high flight altitudes above 10 km and often under cloud-free conditions above HALO, the simulations of F ↓ sol are reliable and can serve as a benchmark due to the implementation of frequent radiosonde and dropsonde observations. Figure 10 compares downward and upward solar irradiance, F ↓ sol and F ↑ sol , measured by BACARDI during the entire flight of 7 February 2020 (flight ID HALO-0207) with alongtrack simulations for cloud-free conditions. To illustrate the effect of the attitude correction for the downward irradiance, data with (black line) and without attitude correction (blue line) are plotted. The uncorrected F ↓ sol shows oscillations of different frequency that are superposed to the diurnal cycle. The slow oscillations (between 0.5 and 1 h) are associated with the circular flight pattern and caused by a combination of a permanent roll angle of about 3 • and changes in latitude (solar zenith angle). Oscillations with higher frequencies between 0.5 and 1 Hz, e.g., most obvious between 17:00 and 19:30 UTC, result from variations of the roll and pitch angle due to turbulence and the aircraft autopilot. The postcorrection of F ↓ sol does remove most of these fast and slow oscillations. This confirms that the roll and pitch angle offsets are determined with sufficient accuracy and that the sensor response time of BACARDI is corrected so that the oscillations are synchronized in time with the aircraft attitude. Subsequently, the remaining slow oscillations are in phase with the simulations and are only caused by the changes in the local solar zenith angle (latitude).
A statistically more robust comparison of measured and simulated F ↓ sol is performed by merging 12 EUREC 4 A flights and filtering the data for cloud-free conditions and reliable attitude corrections. Data are assumed to be valid when the attitude correction factors are less than 25 %. For larger correction factors, roll and pitch angles need to be smaller than 5 • . A one-to-one comparison is shown in Fig. 11a. From the total number of almost 3 million individual measurement samples, 97.6 % agree with the simulations within an uncertainty range of less than 5 %. This indicates that the general performance of BACARDI including the thermal and attitude correction is stable over the entire campaign. A linear regression of all reliable filtered data shows only a slight deviation from the 1 : 1 slope with a correlation coefficient of 0.999. The absolute differences are limited by the applied filter but illustrate that for high solar irradiances, the outliers of the measurements tend to be lower than the simulations. This might be caused by a remaining contamination of the filtered data by clouds above the aircraft, which are not considered in the cloud-free simulations. For low values of F ↓ sol the measurements are slightly overestimated. These measurements correspond to conditions of high solar zenith angle, when the attitude correction becomes more critical. At the same time, the angular response of the CMP22 pyranometer is known to slightly deviate from an ideal cosine response at high solar zenith angles.
Making use of the circular flight pattern, a potential asymmetrical cosine response of the pyranometer inlet is investigated in Fig. 11b. The ratio of the corrected observations and simulations is analyzed as a function of solar zenith angle θ 0 and relative heading of HALO with respect to the solar azimuth. Only a subset of seven flights are used. Other flights are excluded because they either did not contain the circle flight pattern or show evidence of contamination by higher clouds like cirrus. Up to solar zenith angles of approximately 75 • , the observed F ↓ sol is within 5 % of the simulated values. The good agreement of the majority of the data points is regarded as an indicator that the attitude correction is independent of the flight direction over a wider range of illumination conditions (0-1200 W m −2 ) and solar zenith angles (0-75 • ).
For solar zenith angles larger than 75 • a slight directional dependence, relative to the position of the Sun with respect to the orientation of BACARDI (HALO), is obvious. F ↓ sol is overestimated by BACARDI between 30 and 210 • relative solar azimuth and underestimated if the Sun is in the opposite directions. These effects may result from different factors, which cannot be disentangled here. It might indicate slightly incorrect offset angles determined for the attitude correction, an azimuthal dependence of the cosine response of the pyranometer, or reflection by the aircraft fuselage and the tailplane fin at high zenith angles. Therefore, it is advisable to use the data at θ 0 > 75 • with some amount of caution.
The upward solar irradiance as well as the upward and downward terrestrial irradiance cannot be corrected for the aircraft attitude. However, these components are characterized by a nearly isotropic radiation field compared to the downward radiation, and the effects of a misalignment are minimal for a nearly level sensor (Bucholtz et al., 2008). To limit the remaining uncertainties due to the aircraft movement, measurements with roll and pitch angles exceeding ±4 • were removed from the data set. The time series of F ↑ sol shown in Fig. 10b indicates that the flight track covers an area with a generally low cloud cover and some patches of low-level stratiform clouds. The cloud-free areas correspond to the low values of F ↑ sol , which form a baseline at about 80 W m −2 . Only these cloud-free measurements can be compared to radiative transfer simulations (red line). The measurements match the simulated baseline and also follow their slight diurnal change with higher values observed at solar noon. The agreement of observed and simulated F ↑ sol indicates that the measurements in conditions like EUREC 4 A are reliable, even without any attitude correction. For observations over higher reflecting surfaces like sea ice, this needs to be confirmed.

Conclusions
A new radiometer package, the Broadband AirCrAft Ra-Diometer Instrumentation (BACARDI) for the HALO research aircraft, is introduced and characterized. BACARDI comprises two sets of upward-and downward-looking broadband radiometers covering the solar and thermal-infrared spectral ranges. The operation of broadband pyranometers and pyrgeometers as mounted on the HALO research aircraft is investigated in this paper. Especially for a fast and highflying aircraft such as HALO, for which the environmental conditions such as air temperature and density can change rapidly, a minimization of the related dynamic effects is required for the data to fulfill their scientific potential. Three basic corrections are applied to the measurements of BAC-ARDI.
-The post-processing of BACARDI measurements accounts for the temperature dependence of the sensor thermopile sensitivity. Due to the large range of environmental temperatures under which HALO operates (from the surface to the lower stratosphere), this correction amounts to about 5 W m −2 for the pyranometers (1 % change in sensitivity), while the pyrgeometer sensitivity is more stable with less than 1 W m −2 correction (0.5 % change in sensitivity).
-The corrections of the sensor response time make use of the 10 Hz sampling frequency and account for the fast change in irradiance, e.g., in the case of crossing cloud or sea ice edges. The deconvolution method by Ehrlich and Wendisch (2015) with a response time of 1.2 and 3.3 s for the pyranometers and pyrgeometers, respectively, is applied to reconstruct the high-frequency changes in irradiance.
-For the rather smooth changes in the HALO attitude (roll and pitch angle), the common correction method by Bannehr and Schwiesow (1993) is successfully applied to the downward solar irradiance as evaluated during circular flight patterns.
It is shown that known dynamic thermal effects occur for BACARDI when the sensor and dome temperatures do not change simultaneously, such as during ascents and descents into other temperature regimes. To correct for these dynamic thermal offsets, a new method is introduced. Historically, such effects were monitored and corrected with additional measurements of the dome temperature. The approach presented here is based on a simple parameterization that combines the dynamic dome effect and the dynamic thermal offset of the thermopile and, therefore, does not require measurements of the dome temperature. For the radiometers of BACARDI, the dynamic thermal offsets are found to correlate with the rate of change of the sensor temperature, which is expected from theory (see Eq. 17, Sect. 5). Using the sensor temperature as the proxy to determine the dynamic thermal offsets makes the post-processing straightforward as the sensor temperature is measured by the radiometers by default.
The parameterization of the dynamic thermal offset of BACARDI is derived from an exemplary calibration flight in nighttime conditions, in which the pyranometer measurement can be assumed to be zero. For the pyrgeometers, selected flight sections with strong temperature changes are analyzed. The magnitudes of the correction coefficients of the individual radiometers are in the range of 200-500 W m −2 K −1 s and depend on the radiometer type, the mounting position of the radiometer, and the aircraft angle of attack. As the radiometer position and environmental conditions might change between HALO missions, the coefficients should be determined regularly. It also has to be noted that the coefficients reported for BACARDI operated on HALO cannot be transferred to other broadband radiometers on other research aircraft.
The performance of BACARDI was evaluated by measurement examples from the EUREC 4 A field campaign . BACARDI was implemented on HALO for the first time during EUREC 4 A. The system extends the existing suite of active and passive remote sensing instruments on HALO, which lacked instrumentation to observe the solar and thermal-infrared radiative energy budget. BACARDI measurements during an ascent up to 10 km altitude demonstrate how strongly the new dynamic thermal offset affects the single irradiance components in fast-changing environmental conditions. In general, without thermal offset correction, the solar irradiance is underestimated, while the thermal-infrared irradiance is overestimated by up to 20 W m −2 . The exact offset correction depends on radiometer type, the mounting position of the radiometer, and the airflow around the aircraft but is independent of the magnitude of irradiance.
It is shown that net irradiances and atmospheric heating rates calculated from the upward and downward irradiances are less affected by the dynamic thermal effect. As upper and lower radiometers show a similar magnitude of the thermal offset, the thermal effects cancel out to a large extent. In contrast to ascents and descents, for straight flight legs maintaining constant flight levels, which are more typical for HALO observations, the temperature changes are small (below 5 K h −1 ), and potential dynamic thermal offsets range below 1 W m −2 for all broadband irradiances, which appears negligible compared to the uncertainties of the sensor sensitivities (1 % for the CMP22 pyranometer and 4 % for the CGR4 pyrgeometer). Nevertheless, temperature variations and sudden temperature gradients can appear along constant height levels, e.g., at upper-level frontal systems or tropopause disturbances. In conditions with high F ↓ sol , the pyrgeometer shows a slight bias due to leakage of solar radiation above the cut-on wavelength of the CGR4 interference filter. This bias correlates with F ↓ sol and is up to 10 W m −2 during solar noon.
The circular flight patterns frequently performed during EUREC 4 A are used to evaluate the attitude correction of F ↓ sol . Comparisons of the measurements to cloud-free radiative transfer simulations indicate the effectiveness of the corrections. The remaining biases after applying the attitude correction are significant only for solar zenith angles larger than 75 • , which were present during EUREC 4 A only briefly during early or late flights.
The processed broadband irradiances measured by BAC-ARDI during EUREC 4 A are published at the AERIS atmosphere Data and Services Centre . The data are used by Luebke et al. (2022) to assess the cloud radiative forcing with regard to the cloud life cycle and the cloud's temporal evolution, both of which are targets of EUREC 4 A. As shown by Su et al. (2008) the operation of radiometers on high-flying aircraft can exceed the range of environmental conditions for which the radiometer performance is typically certified by the manufacturer. Thus, an additional laboratory characterization of the radiometer might become relevant if operating in high-altitude or Arctic conditions.
Data availability. Processed data of BACARDI are published at the AERIS atmosphere Data and Services Centre (https://doi.org/10.25326/160, . Raw data can be obtained from the authors on request.
Author contributions. AE and MZ equally contributed to the paper. AE compiled the paper, performed the basic corrections of the BAC-ARDI measurements, and evaluated the thermal offset correction. MZ invented the concept of the thermal offset correction, which is presented in this paper, and is responsible for the design BACARDI project and the integration into the aircraft. AG, VN, and CM supported MZ and performed the thermal correction. RM and TR designed the instruments and provided the construction drawing. BS and MW were responsible for the scientific guidance of the project and designed and coordinated the flight strategy. AE, AEL, and KW provided the scientific background and the data analysis for the presented measurement examples.
Competing interests. At least one of the (co-)authors is a member of the editorial board of Atmospheric Measurement Techniques. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.