# Incorporating EarthCARE observations into a multi-lidar cloud climate record: the ATLID (Atmospheric Lidar) cloud climate product

Artem G. Feofilov, Hélène Chepfer, Vincent Noël, Frederic Szczap

► To cite this version:

Artem G. Feofilov, Hélène Chepfer, Vincent Noël, Frederic Szczap. Incorporating EarthCARE observations into a multi-lidar cloud climate record: the ATLID (Atmospheric Lidar) cloud climate product. Atmospheric Measurement Techniques, 2023, 16, pp.3363-3390. 10.5194/amt-16-3363-2023. insu-04195502

HAL Id: insu-04195502

https://insu.hal.science/insu-04195502

Submitted on 4 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Incorporating EarthCARE observations into a multi-lidar cloud climate record: the ATLID (Atmospheric Lidar) cloud climate product

Artem G. Feofilov<fnref n="1" />, Hélène Chepfer<fnref n="1" />, Vincent Noël<fnref n="2" />, and Frederic Szczap<fnref n="3" />

**Correspondence:** Artem G. Feofilov (artem.feofilov@lmd.polytechnique.fr)

Received: 31 October 2022 – Discussion started: 18 November 2022

Revised: 28 February 2023 – Accepted: 15 May 2023 – Published: 5 July 2023

**Abstract.** Despite significant advances in atmospheric measurements and modeling, clouds' response to human-induced climate warming remains the largest source of uncertainty in model predictions of climate. The launch of the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite in 2006 started the era of long-term spaceborne optical active sounding of Earth's atmosphere, which continued with the CATS (Cloud-Aerosol Transport System) lidar on board the International Space Station (ISS) in 2015 and the Atmospheric Laser Doppler Instrument (ALADIN) lidar on board Aeolus in 2018. The next important step is the Atmospheric Lidar (ATLID) instrument from the EarthCARE (Earth Clouds, Aerosols and Radiation Explorer) mission, expected to launch in 2024.

In this article, we define the ATLID Climate Product, Short-Term (CLIMP-ST) and ATLID Climate Product, Long-Term (CLIMP-LT). The purpose of CLIMP-ST is to help evaluate the description of cloud processes in climate models, beyond what is already done with existing space lidar observations, thanks to ATLID's new capabilities. The CLIMP-LT product will merge the ATLID cloud observations with previous space lidar observations to build a long-term cloud lidar record useful to evaluate the cloud climate variability predicted by climate models.

We start with comparing the cloud detection capabilities of ATLID and CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) in day- and nighttime, on a profile-to-profile basis in analyzing virtual ATLID (355 nm) and CALIOP

(532 nm) measurements over synthetic cirrus and stratocumulus cloud scenes. We show that solar background noise affects the cloud detectability in daytime conditions differently for ATLID and CALIPSO.

We found that the simulated daytime ATLID measurements have lower noise than the simulated daytime CALIOP measurements. This allows for lowering the cloud detection thresholds for ATLID compared to CALIOP and enables ATLID to better detect optically thinner clouds than CALIOP in daytime at high horizontal resolution without false cloud detection. These lower threshold values will be used to build the CLIMP-ST (Short-Term, related only to the ATLID observational period) product. This product should provide the ability to evaluate optically thin clouds like cirrus in climate models compared to the current existing capability.

We also found that ATLID and CALIPSO may detect similar clouds if we convert ATLID 355 nm profiles to 532 nm profiles and apply the same cloud detection thresholds as the ones used in GOCCP (GCM-Oriented CALIPSO Cloud Product; general circulation model). Therefore, this approach will be used to build the CLIMP-LT product. The CLIMP-LT data will be merged with the GOCCP data to get a long-term (2006–2030s) cloud climate record. Finally, we investigate the detectability of cloud changes induced by human-caused climate warming within a virtual long-term cloud monthly gridded lidar dataset over the 2008–2034 period that we obtained from two ocean–atmosphere coupled climate models coupled with a lidar simulator. We found that a long-

<footnote n="1">¹LMD/IPSL, Sorbonne Université, ENS, PSL Research University, École polytechnique,
Institut Polytechnique de Paris, CNRS, Paris, France</footnote>

<footnote n="2">²Laboratoire d'Aérologie (Laero), Observatoire Midi-Pyrénées, Université Toulouse 3, CNRS, IRD, Toulouse, France</footnote>

<footnote n="3">³Laboratoire de Météorologie Physique (LaMP), UMR 6016, CNRS, Aubière, France</footnote>
term trend of opaque cloud cover should emerge from short-
term natural climate variability after 4 years (possible life-
time) to 7 years (best-case scenario) for ATLID merged with
CALIPSO measurements according to predictions from the
considered climate models. We conclude that a long-term li-
dar cloud record built from the merging of the actual ATLID-
LT data with CALIPSO-GOCCP data will be a useful tool for
monitoring cloud changes and evaluating the realism of the
cloud changes predicted by climate models.

# 1 Introduction

Clouds play an important role in the radiative energy bud-
get of Earth. The radiative effect of clouds is twofold: on the
one hand, clouds reflect some of the Sun's radiance during
the day, thus preventing surface warming. On the other hand,
high thin clouds trap some of the outgoing infrared radiation
emitted by the surface and re-emit it back to the ground, thus
contributing to its heating. Overall, at a global scale, clouds
contribute to cooling Earth radiatively, but quantifying pre-
cisely this global effect as well as the influence of clouds
on Earth's radiative budget everywhere requires knowing the
coverage of clouds, as well as their geographical and verti-
cal distributions, temperature, and optical properties. Cloud
properties are expected to change under the influence of cli-
mate warming, leading to changes in the amplitude of the
overall cloud radiative cooling. But how cloud properties
change as climate warms is uncertain (e.g., Zelinka et al.,
2012, 2016; Chepfer et al., 2014; Vaillant de Guélis et al.,
2018; Perpina et al., 2021). Cloud feedback uncertainties are
an important contributor to climate sensitivity uncertainty
and therefore limit our ability to predict the future evolution
of climate for a given CO₂ emission scenario (e.g., Winker et
al., 2017; Zelinka et al., 2020).

Global-scale round-the-clock satellite observations of
Earth's atmosphere provide invaluable information that im-
proves our knowledge of current clouds' properties and helps
us to evaluate the cloud description in climate models in cur-
rent climate simulations. Among the remote sensing tech-
niques, active sounding plays a special role because of its
high vertical and horizontal resolution and high sensitiv-
ity. The launch in 2006 of the Cloud-Aerosol Lidar and In-
frared Pathfinder Satellite Observation (CALIPSO; Winker
et al., 2010) satellite started the era of operational spaceborne
optical active sounding of Earth's atmosphere for clouds
and aerosols. It was followed by the CATS (Cloud-Aerosol
Transport System) lidar on board the International Space Sta-
tion (ISS) in 2015 (McGill et al., 2015) and the Atmospheric
Laser Doppler Instrument (ALADIN) lidar on board Aeo-
lus in 2018 (Reitebuch et al., 2020; Straume et al., 2020).
The next important step is the Atmospheric Lidar (ATLID)
instrument (do Carmo et al., 2021), from the EarthCARE
(Earth Clouds, Aerosols and Radiation Explorer) mission

(e.g., Héliere et al., 2017; Illingworth et al., 2015), expected
to launch in 2024. With this lidar, the scientific community
will continue receiving invaluable vertically resolved infor-
mation of atmospheric optical properties needed for the esti-
mation of cloud occurrence frequency, thickness, and height.
Cloud profiles deduced from CALIOP (Cloud-Aerosol Lidar
with Orthogonal Polarization) observations have been widely
used to evaluate the cloud description in climate models (e.g.,
Nam et al., 2012; Cesana et al., 2019) and have provided
leads to improve this description (e.g., Konsta et al., 2012).
To avoid any discrepancy in cloud definition between models
and observations and to allow for consistent comparisons be-
tween clouds simulated by climate models and observed by
satellite, the Cloud Feedback Model Intercomparison Project
(CFMIP) has developed the CFMIP Observation Simulator
Package (COSP1; Bodas-Salcedo et al., 2011), which was
followed by COSP2 (Swales et al., 2018). These packages in-
clude a lidar simulator (Chepfer et al., 2008; Reverdy et al.,
2015; Guzman et al., 2017; Cesana et al., 2019) that mim-
ics the measurements that would be obtained by spaceborne
lidars if they were overflying the atmosphere simulated by
a climate model. In parallel to the COSP lidar simulator,
a Level 2 and 3 cloud product named CALIPSO-GOCCP
(GCM-Oriented CALIPSO Cloud Product; general circula-
tion model; Chepfer et al., 2008, 2010, 2013; Guzman et
al., 2017; Cesana et al., 2019) was designed to ensure scale-
aware and definition-aware comparison between simulated
and observed clouds.

Despite the similarity of the measuring principle of ATLID
and CALIOP lidars – the emitter sends a brief pulse of laser
radiation to the atmosphere, and the receiver registers a time-
resolved backscatter signal collected through its telescope
– the sensitivity of both lidars to the same clouds is differ-
ent. This is explained by differences in observational geome-
try; in wavelength, pulse energy, and repetition frequency; in
telescope diameter and detector type; in the capability of de-
tecting molecular backscatter separately from the particulate
one; in vertical and horizontal resolution and averaging; and
so on. Since the CALIPSO-GOCCP algorithm cannot be ap-
plied directly to ATLID data, a specific algorithm had to be
developed which generates the ATLID cloud product CLIMP
(Climate Product).

The present paper describes the design of CLIMP and its
associated algorithm, developed with the following two goals
in mind.

1. On short timescales, such as the period of the ATLID operation, CLIMP should help improve the current evaluation of cloud description in climate models beyond CALIPSO. From this perspective, CLIMP should take advantage of ATLID capabilities compared to CALIPSO from the point of view of evaluation of clouds in climate models while maintaining compliance with the COSP lidar framework.
2. On long timescales, CLIMP should enable building a merged CALIPSO-ATLID long-term lidar cloud product, in which the same clouds are detected despite the instrumental and orbital differences between ATLID and CALIOP. From this point of view, CLIMP should maximize consistency with GOCCP. The GOCCP-CLIMP long-term dataset should describe more than 20 years of cloud profiles at a global scale, which will enable the study and evaluation in climate models of inter-annual variability in cloud profiles due to multi-annual climate variations (e.g., El Niño, North Atlantic Oscillation, Madden-Julian Oscillation). Its analysis will moreover make possible the detection of cloud changes because of human-induced climate warming, as well as their evaluation in climate model simulations.

Therefore, CLIMP will be composed of two datasets named CLIMP-ST (Short-Term) and CLIMP-LT (Long-Term). They will mainly differ in their cloud detection threshold, as we will see later in the text. This threshold is parameterized in COSP lidar and can easily be changed when comparing simulated data to CLIMP-ST and CLIMP-LT.

The CLIMP product and algorithm inherit from the approach developed for CALIPSO-GOCCP. This algorithm processes Level 1 (L1) data in exactly the same way as the COSP lidar simulator does. GOCCP is part of the CFMIP-OBS (observations) database included in Obs4MIPS (Observations for Model Intercomparisons Project; Waliser et al., 2020) for model evaluation. Differences between GOCCP, NASA, and JAXA (Japan Aerospace Exploration Agency) CALIOP cloud products were documented in Chepfer et al. (2013) and Cesana et al. (2019).

The three key elements of the GOCCP algorithm, which need to be kept when developing CLIMP, are the following.

i. Lidar profiles are not averaged horizontally before cloud detection to (1) keep consistency with the subgrid module SCOPS (Subgrid Cloud Overlap Profile Sampler; Klein and Jacob, 1999) included in COSP that is required to respect the Eulerian framework of climate model simulations and (2) avoid overestimation of the cloud fraction in shallow clouds (e.g., Chepfer et al., 2008, 2013; Feofilov et al., 2022).

ii. Lidar measurements are averaged vertically every 480 m to improve the signal-to-noise ratio (SNR) while maintaining consistency with CloudSat data used for comparison with COSP radar outputs (Marchand et al., 2009; Haynes et al., 2007). This value of 480 m can be different in CLIMP as it can be changed in COSP lidar, but averaging the lidar signal vertically before cloud detection should remain the way to increase ATLID SNR when needed for climate mode evaluation.

iii. Cloud detection thresholds are chosen for consistency with COSP lidar and to prevent false cloud detections

in CALIOP L1 daytime data at full horizontal resolution and 480 m averaged vertical resolution. The cloud detection threshold can be modified in CLIMP but then should also be modified in COSP lidar. This threshold needs to be constant over a full dataset and cannot be scene-dependent.

We would like to stress that the main two purposes of this article are (a) to compare two spaceborne lidars in terms of cloud detection and the signal-to-noise ratio for given observational conditions and (b) to develop a method for merging the data from several spaceborne lidars into a continuous cloud record to detect long-term changes and get a seamless cloud climatology. We assume that the calibration of the instruments is performed dynamically on board the satellites and that the calibration coefficients and crosstalk parameters are known with high accuracy. In this case, we can study the theoretically achievable cloud detection for a given experimental setup, which is defined by a number of parameters like telescope diameter, transmission of the system, solar-noise filtering, detector type, and so on. For the sake of simplicity, we do not discuss the depolarized component of the radiation backscattered by particles, assuming that it is backscattered the same way at these wavelengths and that one can always consider a sum of parallel and perpendicular backscatter for cloud detection.

The structure of the article is as follows. In Sect. 2, we briefly describe the differences and similarities between ATLID and CALIOP, the formalism necessary to understand the analysis presented in the next sections, and the cloud variables used in this study. Section 3 describes the physical elements that matter for the development of CLIMP-ST, using synthetic cloud scenes (Sect. 3.1) and a numerical chain which simulates lidar profiles observed by CALIPSO and ATLID over the cloud scenes at full spatial resolution and instantaneous timescales (Sect. 3.2). In this section, we also pay specific attention to the estimates of lidar signal noise. Then we define the cloud detection scheme of CLIMP-ST (Sect. 3.3), and we try to answer whether ATLID might better observe optically thinner clouds in daytime than CALIOP at full horizontal resolution, a useful capability for evaluating the description of cirrus in climate models. Section 4 describes the physical elements that matter for the development of CLIMP-LT. Section 4.1 presents the cloud detection scheme used in CLIMP-LT to detect the same cloud as CALIPSO-GOCCP, despite the instrumental differences between ATLID and CALIOP. Then we analyze a long-term (multi-decadal, monthly averaged), global-scale space lidar virtual dataset built from climate models and COSP lidar simulation (Sect. 4.2) to illustrate how a merged CLIMP-LT-CALIPSO-GOCCP dataset could help evaluate climate models' predictions of multi-decadal cloud changes (Sect. 4.3). We conclude in Sect. 5.
## 2 Two spaceborne lidars, lidar equation, and cloud detection

### 2.1 Differences between CALIOP (CALIPSO) and ATLID (EarthCARE) spaceborne lidars

CALIOP, a two-wavelength polarization-sensitive near-nadir-viewing lidar, provides high-resolution vertical profiles of aerosols and clouds (Winker et al., 2010). Its initial orbital altitude was 705 km (now 688 km to match that of Cloud-Sat), and its orbit is inclined at 98.05°. The lidar overpasses the Equator at 01:30 and 13:30 LST (local solar time). It uses three receiver channels: one channel measuring the 1064 nm backscatter intensity and two channels measuring orthogonally polarized components of the 532 nm backscattered signal. Cloud and aerosol layers are detected by comparing the measured 532 nm signal return with the return expected from a molecular atmosphere (see the definitions later). The other instrumental parameters of this lidar are described in Table 1 (see also Fig. 1 of Hunt et al., 2009, for a block diagram of CALIOP).

The goals of the EarthCARE mission are “to retrieve vertical profiles of clouds and aerosols, and the characteristics of their radiative and microphysical properties to determine flux gradients within the atmosphere and fluxes at the Earth’s surface, as well as to measure directly the fluxes at the top of the atmosphere and also to clarify the processes involved in aerosol-cloud and cloud-precipitation-convection interactions” (Héliere et al., 2017; Illingworth et al., 2015). The ATLID instrument on board the EarthCARE satellite will measure the attenuated atmospheric backscatter with a vertical resolution of ~100 and ~500 m in the altitude ranges of 0–20 km and 20–40 km, respectively. ATLID is a polarization-sensitive, high-spectral-resolution lidar (HSRL), which can separate the thermally broadened molecular backscatter (Rayleigh) from the unbroadened backscatter from atmospheric particles (Mie) (Durand et al., 2007; see also Fig. 2 of do Carmo et al., 2021). This helps ATLID retrieve extinction and backscatter vertical profiles without assuming the extinction-to-backscatter ratio (as in CALIOP retrievals), which is poorly known, especially for aerosols (e.g., Rogers et al., 2014).

When considering signal quality and performance, some parameters are in favor of CALIOP (telescope diameter, energy per pulse, solar-filter bandwidth), whereas others favor ATLID (altitude, noise level). In the next section, we show how these differences affect the detectability of clouds. We excluded the multiple scattering coefficient from the table since it is an important and complex parameter of lidar instrument which depends on its several parameters. Instead, we discuss it in a dedicated paragraph below.

### 2.2 Lidar equation

The formalism used in this work was described in Feofilov et al. (2022). In this section, we repeat only the basic definitions needed for understanding the material presented below. Since we will discuss both conventional (non-HSRL) and HSRL lidars, we will introduce necessary quantities in parallel and label them correspondingly: the molecular, particulate, and total components will get the indices “mol”, “part”, and “tot”, respectively.

An atmospheric lidar sends a brief pulse of laser radiation towards the atmosphere. The lidar optics collect the backscattered photons and drive them to a detector. The detected signal is time-resolved: supposing each photon traveled straight forward and back, each time bin corresponds to a fixed distance from the lidar to the atmospheric layer where backscattering occurred. The propagation of laser light through the atmosphere and backwards to the detector is described by the following lidar equation:

$$ \text{ATB}(\lambda, z) = (\beta_{\text{mol}}(\lambda, z) + \beta_{\text{part}}(\lambda z)) \times e^{-2\int_{z_{\text{sat}}}^{z} (\alpha_{\text{mol}}(\lambda, z') + \eta\alpha_{\text{part}}(\lambda, z')) dz'} \quad (1) $$

where ATB stands for attenuated total backscatter (m⁻¹ sr⁻¹), βₘol (λ, z) and βₚart (λ, z) are the wavelength-dependent molecular and particulate backscatter coefficients (m⁻¹ sr⁻¹), αₘol (λ, z) and αₚart (λ, z) are the extinction coefficients (m⁻¹), Zₚsat is the altitude of the satellite, λ is the wavelength, and η is a multiple scattering coefficient (e.g., Platt, 1973; Garnier et al., 2015; Donovan, 2016).

For the HSRL lidar, one can write similar equations for the attenuated radiance backscattered from atmospheric particles and molecules (APBs and AMBs), respectively:

$$ \text{APB}(\lambda, z) = \beta_{\text{part}}(\lambda z) \times e^{-2\int_{z_{\text{sat}}}^{z} (\alpha_{\text{mol}}(\lambda, z') + \eta\alpha_{\text{part}}(\lambda, z')) dz'} \quad (2) $$

$$ \text{AMB}(\lambda, z) = \beta_{\text{mol}}(\lambda, z) \times e^{-2\int_{z_{\text{sat}}}^{z} (\alpha_{\text{mol}}(\lambda, z') + \eta\alpha_{\text{part}}(\lambda, z')) dz'} \quad (3) $$

For cloud definition, we will also need to define the attenuated molecular backscatter for clear-sky conditions:

$$ \text{ATB}_{\text{mol}}(\lambda, z) = \beta_{\text{mol}}(\lambda, z) \times e^{-2\int_{z_{\text{sat}}}^{z} \alpha_{\text{mol}}(\lambda, z') dz'} \quad (4) $$

The physical meaning of η in Eqs. (1)–(3) is an increase in the number of photons remaining in the lidar receiver field of view (FOV) besides the ones directly backscattered by the layer, and its value depends on the type of scattering media, FOV of the telescope, and laser beam divergence. The typical value of η varies between 0.5 and 0.8 for commonly used lidars (Chiriaco et al., 2006; Chepfer et al., 2008, 2013; Garnier et al., 2015; Donovan, 2016; see also Appendix B of Reverdy et al., 2015). Setting η to 1 means no multiple scattering and would correspond to an infinitely narrow FOV telescope combined with an infinitely small laser beam divergence. In CALIOP cloud products up to version 3, the η was set to 0.6 for ice clouds, whereas for version 4.10 a
<table><caption>**Table 1.** Specifications of the CALIOP and ATLID spaceborne lidars considered in this article. We gathered specifications from Hunt et al. (2009) for CALIOP and from do Carmo et al. (2021) for ATLID. PMT: photomultiplier tube, APD: avalanche photodiode.</caption><thead><tr><th>Parameter</th><th>Symbol</th><th>CALIOP</th><th>ATLID</th></tr></thead><tbody><tr><td>Altitude (km)</td><td>Z</td><td>688<sup>a</sup></td><td>393</td></tr><tr><td>Orbital inclination (°)</td><td>I</td><td>98.05</td><td>97.050</td></tr><tr><td>Wavelength (nm)</td><td>λ</td><td>532, 1064</td><td>355</td></tr><tr><td>Pulse repetition frequency (Hz)</td><td>PRF</td><td>20</td><td>51 (25.5)<sup>b</sup></td></tr><tr><td>Horizontal distance between profiles (m)</td><td>Δx</td><td>333</td><td>285</td></tr><tr><td>Finest vertical resolution (troposphere) (m)</td><td>Δz</td><td>30</td><td>100</td></tr><tr><td>Telescope diameter (m)</td><td>d<sub>tel</sub></td><td>1.0</td><td>0.6</td></tr><tr><td>Telescope field of view (μrad)</td><td>φ</td><td>130</td><td>64</td></tr><tr><td>Energy per pulse (mJ)</td><td>E<sub>pulse</sub></td><td>110</td><td>35 (70)<sup>b</sup></td></tr><tr><td>Footprint (m)</td><td>d<sub>fp</sub></td><td>90</td><td>29</td></tr><tr><td>Laser beam divergence (μrad)</td><td>θ</td><td>100</td><td>45</td></tr><tr><td>Solar-filter bandwidth (nm)</td><td>Δλ</td><td>0.04, 0.475</td><td>0.71 (0.35)<sup>c</sup></td></tr><tr><td>Solar-filter transmission</td><td>ξ<sub>filter</sub></td><td>0.85</td><td>0.87</td></tr><tr><td>Total optical-system transmission</td><td>ξ<sub>rec</sub></td><td>0.67, 0.68</td><td>0.62</td></tr><tr><td>Detector type</td><td></td><td>PMT, APD</td><td>CCD</td></tr><tr><td>Detector efficiency</td><td>γ</td><td>0.11, 0.4</td><td>0.79 mol, 0.75 part</td></tr><tr><td>Excess noise factor</td><td>ENF</td><td>1.46</td><td>1.44</td></tr><tr><td>Single-shot noise scale factor (square root of photoelectrons)</td><td>NSF</td><td>5.14</td><td>1.0</td></tr><tr><td>Dark current (photoelectrons per second)</td><td>N<sub>dark</sub></td><td>1331, 1.85 × 10<sup>7</sup></td><td>153</td></tr><tr><td>Readout noise (photoelectrons)</td><td>RON</td><td>3–5</td><td>&lt; 3</td></tr></tbody></table>

temperature-dependent coefficient was used, which varied in between 0.46 and 0.78 (Young et al., 2018). For water clouds, the η values are derived from the relationship developed in Hu et al. (2006) (also see Table 4 in Young et al., 2018). A detailed modeling of η for different cloud types observed by CALIOP and ATLID (Shcherbakov et al., 2022) shows that η depends on the cloud thickness and type and that the ATLID values are somewhat higher than those of CALIOP. Based on these works, we set a fixed value of η to 0.6 for CALIOP and to 0.75 for ATLID. This is an approximation, and a more complex approach might be required for processing real data, but our tests show that the conclusions of the present work do not change if we vary η within ±0.1 either for CALIOP or for ATLID (but not for both).

## 2.3 Cloud detection and cloud variables

To characterize the scattering properties of the atmosphere, it would be convenient to use some ratio of attenuated backscatter values (Eqs. 1–4), which would have a clear physical interpretation. Due to attenuation of AMBs below the clouds, using it in the denominator is counterproductive, so the ATB<sub>mol</sub>(λ, z) is used instead, and the scattering ratio (SR) is defined as

$$SR(532 \text{ nm}, z) = \frac{\text{ATB}(532 \text{ nm}, z)}{\text{ATB}_{\text{mol}}(532 \text{ nm}, z)}. \qquad (5)$$

Considering a single-pulse profile measurement, we define a layer as cloudy if the following two conditions are met:

$$SR(532 \text{ nm}, z) > 5$$

and

$$\text{ATB}(532 \text{ nm}, z) - \text{ATB}_{\text{mol}}(532 \text{ nm}, z) > 2.5 \times 10^{-6} \text{ m}^{-1} \text{ sr}^{-1}. \quad (6)$$

The second condition in Eq. (6) comes from the fact that the molecular backscatter in the upper troposphere is weak and the fluctuation in ATB might cause a false cloud detection if only SR is used. With the second condition, the cloud detection is more robust. This definition is used in CALIPSO-GOCCP (e.g., Chepfer et al., 2008, 2010, 2013), and we suggest keeping it for other lidars to ensure consistency between cloud products as discussed later.

In application to ATLID, this will mean using the values recalculated to 532 nm of ATB, which will be estimated from Eq. (1): β<sub>part</sub>(355 nm, z) and α<sub>part</sub>(355 nm, z) retrieved from the measurements (Eqs. 2 and 3) and β<sub>mol</sub>(532 nm, z) and α<sub>mol</sub>(532 nm, z) retrieved or estimated from pressure-temperature profiles from reanalysis. In the numerical experiment below, we calculated ATB<sub>mol</sub>(532 nm, z) = β<sub>mol</sub>(532 nm, z) × e<sup>-2∫z<sub>sat</sub>α<sub>mol</sub>(532, z')dz'</sup> using the available pressure-temperature profiles and the formalism provided in Feofilov et al. (2022). Here, we reproduce Eq. (8) of this paper.

<footnote><sup>a</sup> The nominal orbit altitude at launch was 705 km but was lowered to 688 km in September 2018 to maintain formation flying with CloudSat. <sup>b</sup> The original pulse repetition frequency of ATLID laser is 51 Hz at the energy of 35 mJ per pulse, but the measurements will be doubled on board the satellite (do Carmo et al., 2021), so one can consider the effective frequency and energy per pulse to be equal to 25.5 and 70 mJ, respectively. <sup>c</sup> The solar-filter bandwidth of ATLID is 0.71 nm, but the transmission function of the Mie channel is approximately half of that, so one should calculate the solar noise in this channel with a narrower effective filter width.</footnote>
per:

$$SR'(532 \text{ nm}, z) = \frac{e^{-2\int_{z_{\text{sat}}}^{z} (\alpha_{\text{mol}}(532, z') + \eta_{355} \alpha_{\text{part}}(535, z')) dz'}}{\beta_{\text{mol}}(532, z) \times e^{-2\int_{z_{\text{sat}}}^{z} \alpha_{\text{mol}}(532, z') dz'}} \quad (7)$$

In this conversion, we assume that the spectral dependence of particulate backscatter ($\beta_{part}(\lambda)$ and $\alpha_{part}(\lambda)$) is weak at the wavelengths used in this study. In Beyerle et al. (2001) it is stated that this is generally true for cirrus. In Voudouri et al. (2020), the values at two wavelengths agree within (relatively large) error bars. Therefore, we do not attempt to compensate for the spectral dependence. The only area where we noticed that this approach does not work for real data is the polar stratospheric region, where a direct application of Eq. (7) leads to an overestimation of polar stratospheric clouds (PSCs) (Fig. 8b of Feofilov et al., 2022).

For the cloud properties, we use the same variables as in CALIPSO-GOCCP (Chepfer et al., 2010): cloud fraction CF(z), opaque cloud cover $C_{opaque}$, and opaque cloud height $Z_{opaque}$. If a given atmospheric layer was observed multiple times or if it was sampled vertically at several points, we define the cloud fraction profile CF(z) in a usual way:

$$CF(z) = \frac{N_{\text{cld}}(z)}{N_{\text{tot}}(z)}, \quad (8)$$

where $N_{cld}(z)$ is the number of times the conditions of Eq. (6) is met and $N_{tot}(z)$ is the total number of measurements in this layer. The opaque cloud cover $C_{opaque}$ is used in long time series and is defined over the 2° × 2° latitude–longitude gridded data as follows:

$$SR(z) < 0.06; C_{\text{opaque}} = \frac{N_{\text{opaque\_prof}}}{N_{\text{total\_prof}}}, \quad (9)$$

where the first condition triggers the opaque cloud detection (Guzman et al., 2017), $N_{opaque\_prof}$ is the number of vertical profiles for which an attenuation corresponding to a presence of opaque cloud was found, and $N_{total\_prof}$ is the total number of measurements in a 2° × 2° grid box. For an individual lidar profile, $Z_{opaque}$ corresponds to an altitude of full attenuation of backscattered signal, whereas for gridded data, $Z_{opaque}$ is an opaque-cloud-cover-weighted sum (Guzman et al., 2017).

## 3 The CLIMP short-term dataset

In this section, we search for useful cloud information regarding model evaluation that can be retrieved from the ATLID but cannot be obtained from the CALIPSO data. For this purpose, we use high-resolution cloud scenes (Sect. 3.1), simulate how they would be observed by ATLID and CALIPSO (Sect. 3.2), and compare the SR(z) profiles seen by the two lidars (Sect. 3.3) and the clouds detected by the two instruments (Sect. 3.4). To address the comparability of

clouds observed by two spaceborne lidars, we used the existing methodology (Reverdy et al., 2015; Feofilov et al., 2022) but with a much finer-scaled cloud model, updated instrumental parameters of ATLID, and a new simulation chain which estimates noise at the detector level and propagates it to the cloud product level (the details are provided in Sect. 3.2.2 below). The main question we sought to answer in this section was whether ATLID can better observe optically thinner clouds than CALIPSO in daytime, a useful capability for evaluating thin cirrus clouds in climate models (e.g., Berry et al., 2019). At the same time, we checked whether the chosen cloud detection parameters and instrumental properties affect the detection of highly inhomogeneous low-level thick clouds.

### 3.1 Cloud-generating model

The 3DCLOUD model (Szczap et al., 2014) generates three-dimensional (3D) spatial structures of stratocumulus and fair-weather cumulus and cirrus that share some statistical properties observed in real clouds such as the inhomogeneity parameter $\rho$ (standard deviation normalized by the mean of the water content) and the Fourier spectral slope $\hat{\beta}$ close to -5/3 between the smallest scale of the simulation to the outer scale $L_{out}$ (where the spectrum becomes more flat). We assume that water content follows a gamma distribution. 3DCLOUD_V2 presented in Alkasem et al. (2017) is based on wavelet framework instead of the Fourier framework. First, 3DCLOUD assimilates meteorological profiles (humidity, pressure, temperature, and wind velocity) and solves drastically simplified basic atmospheric equations in order to simulate 3D water content. Second, the Fourier filtering method is used to constrain the intensity of mean water content, $\rho$, $\hat{\beta}$, and $L_{out}$, which are values provided by the user (Hogan and Illingworth, 2003; Kärcher et al., 2018).

Conditions of simulations to generate the stratocumulus in this study (see Fig. 1) are identical to those used in Szczap et al. (2014) for the DYCOMS2-RF01 case (the first research flight of the second Dynamics and Chemistry of Marine Stratocumulus) for the marine stratocumulus regime (Stevens et al., 2005). We have only changed the number of voxels in the x, y, and z directions to $N_x = N_y = 1000$ and $N_z = 50$, respectively. The corresponding spatial resolutions were set to $\Delta_x = \Delta_y = 100$ m and $\Delta_z = 24$ m, respectively. The vertical extension of the simulated area is still $L_z = 1200$ m, but the horizontal extensions for this study are $L_x = L_y = 100$ km.

If the number of voxels is large, the 3DCLOUD and 3DCLOUD_V2 are very time-consuming (see Table 1 in Szczap et al., 2014) and cannot assimilate the fractional coverage for cirrus clouds. Therefore, we developed 3DCLOUD_V3 to overcome these two drawbacks for a cirrus cloud. This model will be published elsewhere. Here, we present only an outline of the 3DCLOUD_V3 algorithm.

To increase the calculation speed in 3DCLOUD_V3, we generate clouds using modified statistical tools developed as
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Examples of the stratocumulus generated with 3DCLOUD: (a) 2D ice water path and (b) its volume rendering.</figcaption></figure>

part of stage 2 of 3DCLOUD. The first stage of 3DCLOUD (i.e., the step of solving simplified basic atmospheric equations, which is very time-consuming) is no longer carried out in 3DCLOUD_V3. Thereby, 3DCLOUD_V3 can be seen as a purely stochastic cirrus cloud generator. The user has to provide, in addition to $N_x, N_y, N_z, \Delta_x, \Delta_y$, and $\Delta_z$, the mean ice water path (IWP); $L_{\text{out}}$; the shape of the vertical profile of ice water content (IWC), $\rho(z)$, and $\hat{\beta}(z)$ and of horizontal wind velocity components $u(z)$ and $v(z)$; and finally the cloud fraction CF. The shape of the vertical profile of IWC can also be stipulated (rectangular, upper triangle, lower triangle, and isosceles trapezoid; Feofilov et al., 2015). The algorithm works as follows:

1. The 3D isotropic field with a Gaussian probability density function (PDF) is generated from a 3D inverse Fourier transform assuming a random phase for each Fourier amplitude and a 3D spectral energy density with 1D spectral slope $\hat{\beta}$ close to -5/3 between the smallest scale of $L_{\text{out}}$.

2. The 2D Gaussian PDF is transformed into a 2D Gamma PDF at each z level, satisfying the values of IWC(z), $\rho(z)$, and $\hat{\beta}(z)$.

3. There is horizontal displacement, at each z level, of 2D IWC (to simulate fall streaks) computed from $u(z)$ and $v(z)$, based on the model of sedimentation proposed by Hogan and Kew (2005). In 3DCLOUD_V3, the user can choose the value of the sedimentation velocity: either constant or as a function of IWC (see formula in Fig. 12 in Heymsfield et al., 2017). Alternatively, the wind velocity vertical profile can be computed from a constant value of the vertical wind shear prescribed by the user; in this case, the user also has to provide the “generated-level height” as explained in Hogan and Kew (2005).

4. The vertical profile of the cloud cover is iteratively modified in order to obtain the CF value prescribed by the user.

Figure 2 demonstrates the examples of 2D IWP and the 3D IWC volume rendering of the cirrus generated with 3DCLOUD_V3, where $N_x = N_y = 1000$, $N_z = 100$, $\Delta_x = \Delta_y = 100$ m, and $\Delta_z = 20$ m. The mean IWP is set to $1 \text{ g m}^{-2}$. The IWC vertical profile shape is “rectangular”. The geometric depth is 2 km. The outer scale is $L_{\text{out}} = 20$ km. We set the constant vertical wind shear to $5 \text{ m s}^{-1} \text{ km}^{-1}$ in the x and y directions, and the generated-level height is 400 m under the cloud top. The inhomogeneity parameter of IWC is $\rho = 0.4$. The spectral slope $\beta$ is equal to -5/3. Figure 2c shows the gamma-like PDF of the IWC (we ignored null values), and Fig. 2d shows the mean power spectra of IWP (and IWC) along x and y directions (and z direction), with the 1D spectral slope close to -2.0 (-1.3) between $L_{\text{out}} = 20$ km and the finest spatial resolution. As expected, values of the spectral slope of IWP are smaller than those of IWC (i.e., the IWP signal is “smoother” than the IWC signal) because IWP is the vertically integral quantity of IWC. One can note that the IWC spectral slope is slightly smaller than the prescribed theoretical value $\hat{\beta} = -5/3$ because of the many null values of the IWC; we plan to remove this bias in the final version of 3DCLOUD_V3.

## 3.2 Numerical chain to simulate cloud observations by CALIOP and ATLID at high resolution

### 3.2.1 Creating pseudo orbits

We performed the following numerical experiment, outlined in the flowchart in Fig. 3. First, we created a gridded global atmosphere from the output of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) atmosphere model version 1 (EAMv1; Rasch et al., 2019) for the conditions of autumn equinox in the Northern Hemisphere. Since we wanted to address both high- and low-level cloud detection, we picked up only the tropical part of the orbit between 5° S and 5° N and used this data as a set of smooth “background” profiles. Since this model does not provide the small-scale variability needed for our experiment, we used
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Examples of cirrus generated with 3DCLOUD_V3: (a) ice water path and (b) its volume rendering; (c) IWC PDF; (d) mean 1D power spectrum of IWP (red curves) and of IWC (blue curve) following the x, y and z direction (solid, dash-dot and dashed line, respectively). A theoretical power spectrum with spectral slope $\hat{\beta} = -5/3$ is added (dashed black line). A dotted vertical black line indicates the outer scale $L_{out} = 20$ km.</figcaption></figure>

the subgrid model described in Sect. 3.1, which generates realistic cloud profiles at a grid comparable or finer than the distance between two consecutive footprints of studied lidars. To address the most challenging observation conditions, we picked up two cloud types: (1) thin cirrus clouds with optical depths ($\tau$) of about 0.03–0.1 per layer (Sassen and Comstock, 2001) and (2) stratocumulus clouds with their high horizontal variability and large optical depths (up to 30 but with about one-third of semi-transparent clouds). These clouds were simulated using an updated 3DCLOUD_V3 model (see Sect. 3.3) and provided as gridded sets of ice water content (IWC) and liquid water content (LWC) values for cirrus and stratocumulus clouds, respectively. We do not consider another challenging case, a thin-cloud layer above a highly reflective cloud, but the daytime noise estimated for the stratocumulus scene will give an idea of what background noise will be interfering with the useful cloud signal in this case.

These gridded sets were converted to pseudo orbits by slicing them along the diagonal lines and arranging the slices into “lidar curtains”, each comprising 20 000 individual profiles and split into daytime and nighttime parts with 10 000 profiles each. This way we got almost seamless cloud distributions, which followed the variability prescribed by the 3DCLOUD_V3 model and at the same time resembled parts

of real lidar orbits. We show the most representative parts of these pseudo orbits in Figs. 4 and 5 for cirrus and stratocumulus clouds, respectively, and we discuss them below.

With these two datasets covering both the daytime and the nighttime scenes, we performed a full series of simulations, explained in Fig. 3. Namely, we fed the high-resolution atmospheric inputs described above to the CALIOP and ATLID simulators (Chepfer et al., 2008; Reverdy et al., 2015) included in the Cloud Feedback Model Intercomparison Project Observational Simulator Package version 2 (COSP2) simulator (Swales et al., 2018). These simulators do not account for instrumental noise effects, so their outputs were processed by a third part of the simulation chain (Fig. 3), which estimates noise and its propagation in the lidar system.

### 3.2.2 Estimating lidar signals and noise

As mentioned above, the outputs of COSP2 simulator are the noise-free APB (355 nm, z) and AMB (355 nm, z) profiles for ATLID and noise-free ATB (532 nm, z) profiles for CALIOP, both calculated at a horizontal resolution of 300 m and vertical resolution of 20 m. To estimate the noise for these profiles and to propagate it further to SR (532 nm, z) for CALIOP and to recalculated SR'(532 nm, z) for ATLID,
<figure><img src="image_5.png" /><figcaption>Figure 3. Flowchart explaining the numerical experiment on comparing clouds retrieved from CALIOP and ATLID observations. Green boxes list the input and output data. Black text between boxes describes actions performed on each dataset. Blue text in the boxes marks the datasets used in the estimation. White text in square brackets in the boxes indicates horizontal (H) and vertical (V) resolutions of the datasets. Note that the ATLID SR' values are estimated at 532 nm (see Sect. 2.3).</figcaption></figure>

we calculated the signals from scratch using the corresponding instrumental parameters introducing measurement noise as follows.

We start with a laser emission and estimate the number of emitted photons per sounding pulse, measured at the output of the sounding unit:

$$N_{\text{em}} = \frac{E_{\text{pulse}}}{hc/\lambda} \times \zeta, \qquad (10)$$

where $E_{\text{pulse}}$ is the energy per laser pulse, $h$ is Planck's constant, $c$ is the speed of light, $\lambda$ is the wavelength, and $\zeta$ is an effective coefficient of optical throughput of the emission path of the lidar. In the present work, it is assumed that there is no optical loss in the emission path and $\zeta$ is equal to 1. The numerical solution of the Eqs. (1)–(3) yields the number of photons per range gate ($t_i t_i + \Delta t_i$), which we will denote as $t_i$, coming through the CALIOP receiver or through the ATLID receiver before splitting into molecular and particulate components in the HSRL module (compare to Eqs. 2 and 3):

$$N_{\text{rec}}^{\text{tot}}(t_i) = N_{\text{rec}}^{\text{mol}}(t_i) + N_{\text{rec}}^{\text{par}}(t_i), \qquad (11)$$

$$N_{\text{rec}}^{\text{mol}}(t_i) = N_{\text{em}} \times \beta_{\text{mol}}(i) \times \Delta z_i \times \Omega(z_i) \times \xi_{\text{rec}} \times e^{-2\sum_{j=0}^{i}(\alpha_{\text{mol}}(j) + \eta\alpha_{\text{part}}(j))\Delta z_j}, \qquad (12)$$

$$N_{\text{rec}}^{\text{par}}(t_i) = N_{\text{em}} \times \beta_{\text{par}}(i) \times \Delta z_i \times \Omega(z_i) \times \xi_{\text{rec}} \times e^{-2\sum_{j=0}^{i}(\alpha_{\text{mol}}(j) + \eta\alpha_{\text{part}}(j))\Delta z_j}, \qquad (13)$$

$$\Omega(z_i) = \pi \times \frac{(d_{\text{tel}}/2)^2}{z_i^2}, \qquad (14)$$

where $N_{\text{rec}}^{\text{mol}}(t_i)$ and $N_{\text{rec}}^{\text{par}}(t_i)$ are the photons backscattered by molecules and particles, respectively; $\Omega(z_i)$ is an altitude-dependent solid angle with $z_i$ corresponding to the time of flight $t_i$ between the satellite and the measured layer $i$; and $\xi_{\text{rec}}$ is the receiver's transmission.

For an HSRL lidar, the molecular and particulate components are supposed to be registered individually, but this separation is not ideal because of the crosstalk between the channels: a part of molecular backscatter comes at the same wavelength as the original laser radiance and it “contaminates” the particulate channel, which is centered at this wavelength. Overall, the HSRL system is characterized by four crosstalk coefficients, $C_{\text{mm}}$, $C_{\text{pp}}$, $C_{\text{mp}}$, and $C_{\text{pm}}$. The first two show a contribution of the molecular and particulate channels to themselves, and in the ideal HSRL they should be equal to 1. The second pair shows how much energy “leaks” from a molecular channel to a particulate one and vice versa. In an ideal HSRL, these coefficients would be equal to 0. In the operational retrieval, these coefficients will be determined through a continuous calibration procedure performed on the orbit. For this exercise, we estimated these coefficients from the Fabry-Pérot interferometer spectral curves (Cheng
<figure><img src="image_6.png" /><figcaption>Figure 4. Example of cirrus cloud input data from the 3DCLOUD model used in the simulation: (a) ice water content (IWC), with night corresponding to one piece of the orbit; (b) IWC, with day corresponding to another piece of the orbit; (c) accumulated optical depth starting from the cloud top at night; (d) same as (c) during the day; (e) cloud mask at night; and (f) cloud mask during the day. We set the cloud mask to 1 whenever $IWC > 0$. The cloud masks presented here are called the “reference dataset” in the rest of the paper.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 5. Same as Fig. 4 but for stratocumulus cloud scenes. Note the color scale difference between Figs. 4 and 5.</figcaption></figure>

et al., 2013): $C_{mm} = 0.815$, $C_{pp} = 0.60$, $C_{mp} = 0.185$, and $C_{pm} = 0.40$. In the case of a non-ideal HSRL, the number of photoelectrons produced by each detector per range gate is as follows:

$$N_{\text{det}}^{\text{mol}}(t_i) = \gamma \times \xi_{\text{rec}} \times (N_{\text{rec}}^{\text{mol}}(t_i) \times c_{\text{mm}} + N_{\text{rec}}^{\text{par}}(t_i) \times c_{\text{pm}}), \quad (15)$$

$$N_{\text{det}}^{\text{par}}(t_i) = \gamma \times \xi_{\text{rec}} \times (N_{\text{rec}}^{\text{mol}}(t_i) \times c_{\text{mp}} + N_{\text{rec}}^{\text{par}}(t_i) \times c_{\text{pp}}), \quad (16)$$

where $\gamma$ and $\xi_{\text{rec}}$ are the detector’s quantum efficiency and transmittance of the optical path, respectively. To come back to “pure” $N_{\text{rec}}^{\text{mol}}(t_i)$ and $N_{\text{rec}}^{\text{par}}(t_i)$ used in the retrieval, one has

to solve this system:

$$N_{\text{rec}}^{\text{mol}}(t_i) = N_{\text{det}}^{\text{mol}}(t_i) \times k_a + N_{\text{det}}^{\text{par}}(t_i) \times k_b, \quad (17)$$

$$N_{\text{rec}}^{\text{par}}(t_i) = N_{\text{det}}^{\text{mol}}(t_i) \times k_c + N_{\text{det}}^{\text{par}}(t_i) \times k_d, \quad (18)$$

$$k_a = \frac{c_{\text{pp}}}{\kappa}; k_b = \frac{-c_{\text{pm}}}{\kappa}; k_c = \frac{-c_{\text{mp}}}{\kappa}; k_d = \frac{c_{\text{mm}}}{\kappa}; \\ = \gamma \times \xi_{\text{rec}} \times (c_{\text{mm}}c_{\text{pp}} - c_{\text{pm}}c_{\text{mp}}). \quad (19)$$

Besides the components related to atmospheric backscatter properties, the $N_{\text{det}}^{\text{mol}}(t_i)$ and $N_{\text{det}}^{\text{par}}(t_i)$ are affected by “parasite” solar backscattered photons during the daytime, which are not correlated with the laser shots. To estimate the solar background add-on to $N_{\text{det}}^{\text{mol}}(t_i)$ and $N_{\text{det}}^{\text{par}}(t_i)$, one has to solve the radiative transfer equation for the radiation emitted by the Sun, backscattered by air and particles in the atmosphere and by the surface in the direction of the spaceborne lidar, and attenuated by the atmospheric layers:

$$N_{\text{sol}}(t_i) = \Delta t_i \times N_{\text{sol}}^{\text{TOA}}(\lambda) \times [R_{\text{sol}}^{\text{atm}} + R_{\text{sol}}^{\text{surf}}], \quad (20)$$

$$R_{\text{sol}}^{\text{atm}} = \int_{z_{\text{sat}}}^{0} (\beta_{\text{m}}(z) \times \phi_{\text{m}}(\text{SZA}) + \beta_{\text{p}}(z) \times \phi_{\text{p}}(\text{SZA})) \times \cos(\text{SZA})^{-1} \\ \times \exp \left\{ -2 \int_{z_{\text{sat}}}^{z} (\alpha_{\text{m}}(z') + \alpha_{\text{p}}(z')) \cos(\text{SZA})^{-1} dz' \right\} dz, \quad (21)$$

$$R_{\text{sol}}^{\text{surf}} = A_{\text{surf}} \times \phi_{\text{surf}}(\text{SZA}) \times \exp \left\{ - \int_{z_{\text{sat}}}^{z} (\alpha_{\text{m}}(z) + \alpha_{\text{p}}(z)) \cos(\text{SZA})^{-1} dz \right\} \times \int_{z_{\text{sat}}}^{0} (\alpha_{\text{m}}(z) + \alpha_{\text{p}}(z)) dz, \quad (22)$$

where $N_{\text{sol}}^{\text{TOA}}(\lambda)$ is the top-of-atmosphere solar flux at wavelength $\lambda$ and for filter width $\Delta\lambda$; $R_{\text{sol}}^{\text{atm}}$ and $R_{\text{sol}}^{\text{surf}}$ represent the proportion of the incoming solar radiance reflected in the direction of lidar; $\phi_{\text{m}}(\text{SZA})$, $\phi_{\text{p}}(\text{SZA})$, and $\phi_{\text{surf}}(\text{SZA})$ are the scatter plots for the angle between the Sun and the nadir view of lidar for molecular scattering, scattering on particles, and scattering from the surface, respectively; $z_{\text{sat}}$ is the altitude of a satellite; and $A_{\text{surf}}$ is the surface albedo. We assumed Lambertian scattering from the surface with albedo equal to 0.08 for ocean and 0.15 for land (arbitrary values), we used Rayleigh scattering phase function for the molecular component, and we used the geometric-optics phase function approximation for particulate scattering.

The solar photons pass through the optical system and HSRL, hit the detectors, and produce the “solar-noise photoelectrons”:

$$N_{\text{det.sol.}}^{\text{mol}}(t_i) = \gamma \times \xi_{\text{rec}} \times \xi_{\text{mol.sol.}} \times N_{\text{sol}}(t_i), \quad (23)$$

$$N_{\text{det.sol.}}^{\text{par}}(t_i) = \gamma \times \xi_{\text{rec}} \times \xi_{\text{par.sol.}} \times N_{\text{sol}}(t_i), \quad (24)$$

where $\xi_{\text{mol.sol.}}$ and $\xi_{\text{par.sol.}}$ represent the convolution of the solar-filter spectral curve with the interferometric spectral
curve for a given channel (see the comment on solar-filter bandwidth in Table 1). In addition to solar noise, there is always a dark current of the detector ($N_{dark}$) and readout noise (RON) which are added to the signal. Since the $N_{det.sol.}^{mol}$ ($t_i$) and $N_{det.sol.}^{par}$ ($t_i$), $N_{dark}$, and RON are registered along with $N_{det}^{mol}$ ($t_i$) and $N_{det}^{par}$ ($t_i$) during daytime, they enter Eqs. (15) and (16) and affect the retrieval. For the non-HSRL lidar:

$$N_{det.sol.}(t_i) = \gamma \times \xi_{rec} \times \xi_{tot.sol.} \times N_{sol}(t_i), \quad (25)$$

where $\gamma$ and $\xi_{rec}$ stand for the corresponding parameters of non-HSRL lidar and $\xi_{tot.sol.}$ represents the transmission coefficient of a solar rejection filter, which is equal to a ratio of an integral of the spectral transmission curve of the filter to a full spectral width of the filter. A quick back-of-the-envelope estimate of the ratio of solar photons coming to the particulate detector of ATLID to the number of solar photons reaching the surface of CALIOP's detector per same sampling interval is

$$\frac{N_{sol}^{TOA} (355 \text{ nm}) \times d_{tel.ATLID}^2 \times d_{fp.ATLID}^2 \times Z_{CALIOP}^2 \times \xi_{rec.ATLID} \times \xi_{par.sol.}}{N_{sol}^{TOA} (532 \text{ nm}) \times d_{tel.CALIOP}^2 \times d_{fp.CALIOP}^2 \times Z_{ATLID}^2 \times \xi_{rec.CALIOP} \times \xi_{tot.sol.}} = \\ \frac{(1162.8 \times 355) \times 0.6^2 \times 690^2 \times 29^2 \times 0.62 \times (0.35 \times 0.87)}{(1900.0 \times 532) \times 1.0^2 \times 393^2 \times 90^2 \times 0.67 \times (0.04 \times 0.85)} = 0.38. \quad (26)$$

So, at first sight, the ATLID retrieval should be less solar-contaminated than CALIOP. But, this ratio alone is not enough for such a conclusion because the solar photons should be compared to the useful signal. Below, we show the results of simulations for two atmospheric scenarios which consider two-way radiative transfer both for solar radiance and for lidar sounding radiance and add the noise of the remaining detection path.

Now, when all the components of the signal are known, we can estimate the daytime and nighttime signal and noise and propagate them to the retrieved parameters. It is important to mention that the instruments compared in this work use the detectors of different types. Namely, the CALIOP lidar uses a photomultiplier tube (PMT), whereas the ATLID lidar detects the backscatter with the help of a charge-coupled device (CCD). Besides different characteristics like gain or dark current (see Table 1), these detectors are not the same in terms of applicable noise statistics (Liu et al., 2006). Even though the incoming photon flux distributions for both instruments are Poisson, the photoelectrons produced by the PMT do not follow a strict Poisson distribution. It is known that for Poisson-distributed signals, a one-to-one relationship exists between the mean and the variance of the photocurrent. As Liu et al. (2006) show, the mean and the variance of the PMT photocurrent are also proportional but not one to one, and the corresponding noise scale factor (NSF) has to be applied to estimate random errors for lidar systems using PMTs or avalanche detectors. The NSF is linked to an excess noise factor (ENF), but it is not equal to it. For the PMTs with identical gain factors $m$ for each dynode, the ENF is given by (Kingston, 1978; Liu and Sugimoto, 2002)

$$\text{ENF} = \frac{m}{m-1}. \quad (27)$$

For the analog detection, the NSF in the multiplied photoelectron domain can be either calculated from the detector's ENF and gain or estimated from the solar-noise-dominated signals (Liu et al., 2006):

$$\text{NSF} = \frac{\sigma(N_{\text{det}})}{\sqrt{(N_{\text{det}})}} = \frac{\sqrt{\text{var}(N_{\text{det}})}}{\sqrt{(N_{\text{det}})}}, \quad (28)$$

where $\sigma(N_{det})$, $\text{var}(N_{det})$, and $\langle N_{det} \rangle$ are the standard deviation, variance, and mean of the signal, respectively. The NSF value provided in the CALIPSO L1 version 4.10 files is equal to 5.14. However, using this value in synthetic-noise calculations leads to an overestimation of the daytime noise, so for the calculations below we took a more conservative value of NSF = 3.16. One can write the expressions for the variances of CALIOP and ATLID signals in the analog detection domain through the number of photoelectrons calculated for each channel:

$$\text{var}(N_{\text{det}}^{\text{mol}}) = \text{ENF} \times (N_{\text{det}}^{\text{mol}} + N_{\text{dark}} + N_{\text{det.sol.}}^{\text{mol}}) + \text{RON}, \quad (29)$$

$$\text{var}(N_{\text{det}}^{\text{par}}) = \text{ENF} \times (N_{\text{det}}^{\text{par}} + N_{\text{dark}} + N_{\text{det.sol.}}^{\text{par}}) + \text{RON}, \quad (30)$$

$$\text{var}(N_{\text{det}}^{\text{tot}}) = \text{NSF}^2 \times (N_{\text{det}}^{\text{tot}} + N_{\text{dark}} + N_{\text{det.sol.}}^{\text{tot}}) + \text{RON}. \quad (31)$$

We draw the reader's attention to the fact that the detector's parameters in Eq. (31) are not equal to those in Eqs. (29) and (30) and that for real calculations one has to use the values from Table 1 or a similar source. The variances of molecular, particulate, and total incoming photon fluxes, which are finally used in the optical property retrievals, are estimated in accordance with the standard error propagation formulae applied to the equations above:

$$\text{var}(\text{AMB}) = k_a^2 \times \text{var}(N_{\text{det}}^{\text{mol}}) + k_b^2 \times \text{var}(N_{\text{det}}^{\text{par}}) + 2 \times k_a \times k_b \times \text{cov}(N_{\text{det}}^{\text{mol}} N_{\text{det}}^{\text{par}}), \quad (32)$$

$$\text{var}(\text{APB}) = k_c^2 \times \text{var}(N_{\text{det}}^{\text{mol}}) + k_c^2 \times \text{var}(N_{\text{det}}^{\text{par}}) + 2 \times k_c \times k_d \times \text{cov}(N_{\text{det}}^{\text{mol}} N_{\text{det}}^{\text{par}}), \quad (33)$$

$$\text{var}(\text{ATB}) = \frac{1}{\kappa^2} \times \text{var}(N_{\text{det}}^{\text{tot}}), \quad (34)$$

where the $\text{cov}(N_{\text{det}}^{\text{mol}} N_{\text{det}}^{\text{par}})$ represents the covariance of molecular and particulate channels. This term is required because the signals in the channels are coupled through non-zero crosstalk.
When the variances are known, the original noise-free AMB, APB, and ATB profiles are modified by random noise, which is modulated by the standard deviation calculated from the variances, and the results are saved. If horizontal or vertical signal averaging is involved, the noise is scaled inversely proportional to a square root of the number of samples within the averaging interval.

### 3.2.3 Useful lidar signals and their SNRs

To address the information content of the backscattered radiance, it makes sense to define a useful signal and to estimate the SNR for this signal. For CALIOP the useful signal is represented by ATB (λ, z) (see Eq. 1), whereas ATLID can measure the molecular and particulate backscattered radiances separately, so it would be logical to call the APB (λ, z) (see Eq. 2) a signal which carries the information about the cloud and look at its SNR. For the sake of simplicity, we do not discuss here the perpendicular channels of these two space lidars assuming that the backscattered depolarized radiance is detected the same way, and adding the processing of this component to the formalism above would not change the conclusions of this work. Another aspect that we do not discuss here is the change in cloud microphysics, which can also affect the cloud detection and cloud radiative effects. We consider and model only the cloud occurrence, cloud cover, and cloud detectability.

For the simulated CALIOP signals, we estimate SR(z) at 532 nm and CF(z) according to Eqs. (6) and (7). The simulated ATLID signals are converted to equivalent 532 nm SR'(z) (see Sect. 2.2 and Feofilov et al., 2022). Then we calculate CF(z) for ATLID using Eqs. (6) and (7) with the same thresholds, and then we analyze the resulting cloud fraction.

To quantify the lidar cloud detection agreement and disagreement regarding the reference cloud dataset, we distinguish four cases: (1) when the lidar detects the actually cloudy layer as cloudy (YESSiYES case); (2) when there is no cloud and the lidar does not detect a cloud (NO_NO); (3) when the lidar does not detect an existing cloudy layer (YES_NO or false negative); and (4) when the lidar detects a cloud, whereas the layer does not contain a cloud (NOSiYES or false positive). We will define their occurrence ratios as

$$ 
\begin{aligned}
R_{YES\_YES}(z) &= \frac{N_{YES\_YES}(z)}{N_{tot}(z)}; \\
R_{NO\_NO}(z) &= \frac{N_{NO\_NO}(z)}{N_{tot}(z)}; \\
R_{YES\_NO}(z) &= \frac{N_{YES\_NO}(z)}{N_{tot}(z)}; \\
R_{NO\_YES}(z) &= \frac{N_{NO\_YES}(z)}{N_{tot}(z)}.
\end{aligned} 
\quad (35) $$

The sum of all four ratios in Eq. (35) yields unity. A perfect match between the cloud distribution in the atmosphere and the product retrieved from the measurement would be

when $R_{YES\_YES}(z) + R_{NO\_NO}(z) = 1$ and $R_{YES\_NO}(z) = R_{NO\_YES}(z) = 0$.

### 3.3 Simulated ATLID and CALIPSO lidar profiles over cirrus and stratocumulus scenes

The most representative parts of pseudo orbits generated with the help of the 3D_CLOUDV3 model (Sect. 3.3) are shown in Figs. 4 and 5 for cirrus and stratocumulus clouds, respectively. We arbitrarily split the “cloud curtain” generated from the output of this model (Sect. 3.2) to “daytime” and “nighttime” by setting the solar zenith angle (SZA) to 45 and 120°, respectively. These values are not linked with the cloud formation mechanisms in the 3D_CLOUDV3 model; they are just needed for a second half of the simulator chain (see noise-related boxes in Fig. 3). In Fig. 4a and b, one can see a fine structure of modeled cirrus clouds. Looking at Fig. 4c and d, one can say that the clouds are optically thin. This combination makes the detection of the clouds marked in Fig. 4e and f challenging.

The stratocumulus clouds shown in Fig. 5 belong to another category of challenging observations. The clouds are closely spaced along the horizontal axis, and at the same time they are optically thick: about two-thirds of the clouds have optical thickness larger than 3 (Fig. 5c, d), but the scene also contains about one-third of semi-transparent clouds like the ones that were reported in Leahy et al. (2022). From Fig. 5c and d, one can conclude that at present there is no space-based measurement that can retrieve all the optical properties of cloud layers shown in Fig. 5e and f. Another problem of these clouds is that their horizontal averaging might bias the estimated cloud fraction (see, e.g., Fig. 4 and discussion of Feofilov et al., 2022).

In Figs. 6 and 7, we demonstrate the differences between two lidars for four scenes (cirrus/stratocumulus clouds, day/night) using the simulated backscatter signal. For the cirrus cloud scene (Fig. 6), both the ATB (532 nm, z) of CALIOP and the APB (355 nm, z) of ATLID show a detectable signal in the areas marked by a cloud mask in Fig. 6e and f. But if one defines the signal detection level as $3\sigma$, one will see that a part of thin clouds will be missing. This is not surprising since we compare a “pure” modeled cloud with its noisy representation in the measuring system. What can be estimated from the image is the potential reliability of cloud detection from ATLID and CALIOP: according to the SNR values (Fig. 6g, h vs. c, d), the APB (355 nm, z) signal from ATLID (Fig. 6e, f) reaches higher SNR values than the ATB (532 nm, z) signal from CALIOP (Fig. 6a, b). This hints at the fact that the cloud detection from this instrument might be somewhat better than that from CALIOP and that one can lower the detection threshold and still get the cloud instead of noise. This is a subject of one of the experiments described below. As for the daytime vs. nighttime difference, we do not see a big change between the left-hand side and right-hand side panels for ATLID (Fig. 6e–h), whereas the
<figure><img src="image_8.png" /><figcaption>Figure 6. Signals and signal-to-noise ratio for the cirrus cloud scene. CALIOP: (a) ATB (532 nm, z) at night for one piece of the orbit, (b) ATB (532 nm, z) during the day for another piece of the pseudo orbit, (c) SNR at night, and (d) SNR during the day. ATLID: (e) APBs (355 nm, z) at night, (f) APBs (355 nm, z) during the day, (g) SNR at night, and (h) SNR during the day. Note that the scene contains only these clouds and a clear sky below. For the reflective clouds beneath the cirrus layer, the daytime noise will be higher (see the right-hand side panels of Fig. 7).</figcaption></figure>

CALIOP shows higher noise in Fig. 6b and d. We note here that the calculations were performed for the cases when only a thin cirrus cloud was present in the atmospheric column, whereas the rest of it corresponded to clear-sky conditions. In the real world, though, the second cloud layer beneath cirrus might increase the solar noise (see the right-hand side panels of Fig. 7), and this will adversely affect the thin-cloud detection, especially from the CALIOP measurements. This is explained by a larger field of view of the CALIOP lidar (see Table 1). In our exercise, we wanted to estimate the best achievable results for a given cloud scene for each instrument and to compare the lidar performances. This way, the conclusions made below for the daytime scenes refer to the minimal differences between the two instruments. As for the stratocumulus clouds (Fig. 7), both the signals and SNRs are strong for both lidars, day and night. The altitudes beneath these clouds correspond to areas without a useful signal: at these heights, the signal is already attenuated by a cloud above, and the attenuation is so strong that even the cloud base is not visible at optical wavelengths (e.g., Guzman et al., 2017). Another remarkable feature shown in this plot is higher daytime noise for CALIOP (Fig. 7b, d). Even though this high noise level does not affect the stratocumulus cloud detection itself, it might affect the aforementioned higher-level cloud detection, and, from this point of view, ATLID has an advantage over CALIOP.

Summarizing, one can say that the ATB (532 nm, z) signals of CALIOP and the APB (355 nm, z) signals of ATLID carry similar type information for the same cloud scenes, but their SNRs suggest (a) that the daytime cloud detection from ATLID should be more reliable and (b) that one can lower the

<figure><img src="image_9.png" /><figcaption>Figure 7. Same as Fig. 6 but for stratocumulus cloud scenes.</figcaption></figure>

detection threshold for this instrument without admixing numerous noise-triggered clouds. Let us now see how the signal quality transforms into the product quality and, in particular, to cloud detection quality.

## 3.4 Capability of ATLID to better detect optically thinner clouds than CALIPSO

Here, we describe the test we performed seeking to answer whether the cloud detection limits (Eq. 6) defined in Chepfer et al. (2010) could be lowered to detect thinner clouds. For this test, we followed the second half of the flowchart (Fig. 3) and calculated the SR (532 nm, z) for CALIOP and the CALIOP-like SR (532 nm, z) for ATLID (Eqs. 2–3), but we changed the cloud detection thresholds of Eq. (6) to the following ones:

$$ \begin{array}{l} \text{SR}(532\,\text{nm}, z) > 3 \\ \text{and} \\ \text{ATB}(532\,\text{nm}, z) - \text{ATB}_{\text{mol}}(532\,\text{nm}, z) > 1.5 \times 10^{-6} \, \text{m}^{-1} \, \text{sr}^{-1}. \quad (36) \end{array} $$

Then we estimated the cloud fractions and statistical agreement with the source cloud data (Eqs. 6, 7). The threshold in the left-hand side of Eq. (36) implies that the particulate backscatter in a layer, which we call a cloudy one, is twice the molecular one.

The threshold in the right-hand side of Eq. (36) corresponds to the absolute values of ATB (532 nm, z) recalculated for SR (532 nm, z) = 3 at the height of 8 km (Chepfer et al., 2010), but overall the rationale for selecting these very values is based on the SNR values levels we observed in the test simulations. Further lowering the threshold will lead to an increased number of false-positive cloud detections in ATLID.

Since the “native” CALIOP profiles are averaged over three points above 8 km, we applied an averaging procedure over ~1 km distance to all simulated signals and repeated the analysis. To compare apples to apples in terms of signal statistics, we averaged over four CALIPSO shots and over two effective ATLID shots, yielding the actual average over 1330 and 1140 m, respectively. To reduce the number of
<table><caption>**Table 2.** Cloud detection statistics for CALIOP and ATLID when the cloud definition corresponds to SR (532 nm, z) > 3 and ATB (λ, z) – ATB<sub>mol</sub> (λ, z) > 1.5 × 10<sup>-6</sup> m<sup>-1</sup> sr<sup>-1</sup> (Eq. 36). The bias and RMS values are defined for the clouds detected in the columns (see text), and we define the total score in percent as 100% × (1 – (YES_NO + NO<sub>YES</sub>)/(YESSiYES + NO_NO)). Ci: cirrus clouds, Sc: stratocumulus clouds.</caption><thead><tr><th colspan="2">Lidar</th><th colspan="4">CALIOP</th><th colspan="4">ATLID</th></tr><tr><th colspan="2">Day/night</th><th colspan="2">Night</th><th colspan="2">Day</th><th colspan="2">Night</th><th colspan="2">Day</th></tr><tr><th>Cloud</th><th>Ci</th><th>Sc</th><th>Ci</th><th>Sc</th><th>Ci</th><th>Sc</th><th>Ci</th><th>Sc</th><th>Ci</th><th>Sc</th></tr><tr><th>Averaged</th><th>N</th><th>Y</th><th>N</th><th>Y</th><th>N</th><th>Y</th><th>N</th><th>Y</th><th>N</th><th>Y</th></tr></thead><tbody><tr><td>YESSiYES</td><td>7</td><td>7</td><td>8</td><td>8</td><td>8</td><td>7</td><td>9</td><td>8</td><td>8</td><td>8</td></tr><tr><td>NO_NO</td><td>80</td><td>78</td><td>82</td><td>86</td><td>69</td><td>72</td><td>68</td><td>80</td><td>83</td><td>81</td></tr><tr><td>YES_NO</td><td>10</td><td>14</td><td>6</td><td>6</td><td>13</td><td>18</td><td>4</td><td>5</td><td>9</td><td>11</td></tr><tr><td>NOSiYES</td><td>2</td><td>0</td><td>4</td><td>0</td><td>10</td><td>3</td><td>18</td><td>6</td><td>0</td><td>0</td></tr><tr><td>Total score</td><td>85</td><td>85</td><td>88</td><td>93</td><td>56</td><td>71</td><td>69</td><td>87</td><td>91</td><td>89</td></tr><tr><td>Bias</td><td>11</td><td>-4</td><td>-1</td><td>-6</td><td>21</td><td>11</td><td>16</td><td>2</td><td>-5</td><td>-8</td></tr><tr><td>RMS</td><td>18</td><td>13</td><td>8</td><td>6</td><td>23</td><td>21</td><td>9</td><td>6</td><td>10</td><td>9</td></tr></tbody></table>

Of course, the examples considered in this section do not cover the whole range of high-, middle-, and low-level clouds, but they draw a line between the threshold values that can be used for cloud definition for CALIOP and ATLID and show that the difference is linked to noise characteristics of the instruments. This result suggests that ATLID should be able to better observe optically thinner clouds than CALIOP in daytime at full horizontal resolution.

To illustrate this point, we used the available dataset for cirrus and estimated the minimal detectable backscatter (MDB) for ATLID in terms of equivalent ATB (532 nm, z) for comparison with CALIOP values obtained for 5 km horizontal averaging of cirrus measured at 15 km height (McGill et al., 2007). For this numerical experiment, we used noisy APBs and noise-free AMBs to keep the consistency with our approach of cloud detection using only one noisy component, the particulate one. For this horizontal averaging, we obtained MDB = 3.0 ± 1.0 × 10<sup>-7</sup> m<sup>-1</sup> sr<sup>-1</sup> for the nighttime and MDB = 4.0 ± 1.0 × 10<sup>-7</sup> m<sup>-1</sup> sr<sup>-1</sup> for the daytime in equivalent ATB (532 nm) values, whereas for CALIOP we obtained MDB = 4.0 ± 2.0 × 10<sup>-7</sup> m<sup>-1</sup> sr<sup>-1</sup> for the nighttime and MDB = 1.3 ± 0.2 × 10<sup>-6</sup> m<sup>-1</sup> sr<sup>-1</sup> for the daytime in its native ATB (532 nm). The daytime value estimated for CALIOP is in good agreement with the measured one (McGill et al., 2007), whereas the estimated nighttime value is somewhat lower than the measured MDB = 8.0 ± 1.0 × 10<sup>-7</sup> m<sup>-1</sup> sr<sup>-1</sup>. From this comparison, we cannot conclude that the ATLID will provide better sensitivity to thin clouds during nighttime, but we can conclude that its daytime thin-cloud detection at 5 km averaging capacity should be comparable to that of CALIOP for the nighttime, and this will be an important achievement for daytime vs. nighttime cloud comparison. Using the cloud detection thresholds defined by Eq. (36) and refined for the real data flow using the methodology outlined above, the CLIMP short-term product will be produced.

# 4 The CLIMP long-term dataset

## 4.1 Capability of CLIMP and CALIPSO-GOCCP to detect the same clouds

One of the overarching goals of our study is to develop a method for merging the data from several spaceborne lidars into a continuous cloud record to detect long-term changes and get a seamless cloud climatology. Since the low threshold tested in the previous section revealed the sensitivity mismatch between the two instruments, we had to test whether the cloud detection thresholds developed for CALIOP (Chepfer et al., 2010) are applicable to ATLID and whether the clouds retrieved using these thresholds are consistent between the two lidars. For this exercise, we followed the same scheme as in the previous section, but this time the clouds were defined in Eq. (6) as in Chepfer et al. (2010) and the follow-up works (e.g., Cesana et al., 2019; Guzman et al., 2017).

Figure 10 demonstrates the daytime and nighttime scattering ratios above the detection thresholds (Eq. 6) and the corresponding cirrus cloud detection statistics for CALIOP. The SR (532 nm, z) in Fig. 10a and b demonstrates the same patterns as the ATB (532 nm, z) signals in Fig. 6a and b. As expected, this time the daytime noise is less pronounced (compare Fig. 10b to Fig. 8b). Still, the daytime scene contains a certain number of false detections marked in red in Fig. 10d. The same analysis performed for ATLID (Fig. 11) also shows somewhat less noise in daytime (compare Fig. 11b to Fig. 9b). The cloud detection quality of ATLID is like that of CALIOP (see Table 3). In this setup, the ATLID is just slightly better than CALIOP with its somewhat higher rate of false detections during the day (compare panels c and d of Figs. 8, 9, 10, and 11 and the corresponding columns in Table 3). For stratocumulus clouds (Figs. 1 and 5), with their strong signals, the agreement between CALIOP and ATLID
<figure><img src="image_10.png" /><figcaption>**Figure 10.** Scattering ratio and cloud detection statistics estimated for cirrus clouds observed by CALIOP using Eq. (6): (a) scattering ratio at night, (b) cloud detection at night, (c) scattering ratio during the day, and (d) cloud detection during the day. Note the color scale difference for (a) and (c) compared to (b) and (d).</figcaption></figure>

is also better than for the clouds defined by Eq. (36) (compare panels c and d of Figs. A1, A2, A3, and A4). The 1 km averaging further improves the agreement between the datasets (Table 3).

Summarizing, using the thresholds (Eq. 6) to define the clouds makes the cloud datasets from CALIOP and ATLID comparable. Further adjustment will be needed for real ATLID data to compensate the effects of the diurnal cycle (Noel et al., 2018; Chepfer et al., 2019; Feofilov and Stubenrauch, 2019). Other compensations might be required when the real ATLID data become available. Since there is a high chance that there will be no overlapping period for these two satellite instruments, an intercalibration procedure will be required. For this, one can use the average cloud amount for low, middle, and high clouds in different zones (tropics, middle latitudes, and polar) to track the changes and to introduce feedback to the cloud detection algorithm. This way, the number of cases measured for each zone will be high, and the uncertainty will be low. The daytime and nighttime observations should be considered separately to address the diurnal cycle and daytime noise issues. In the sections below, we assume that the intercalibration has been performed and that the cloud datasets agree.

## 4.2 Numerical chain to simulate long-term lidar record and method to quantify time of emergence

The previous section shows that ATLID and CALIOP data may be merged to build a long-term dataset, even though their instrumental or orbital differences might necessitate further reconciliation. Here we suppose perfect reconciliation will eventually be reached, and we build a long-term space lidar synthetic dataset spanning more than 30 years to examine when a change in cloud properties attributable to human-induced warming would be detectable in the lidar

cloud record according to climate model simulations. This approach is directly inspired by the one pioneered in Chepfer et al. (2018) and later expanded in Perpina et al. (2021).

We use climate predictions from IPSL-CM6 (Institut Pierre Simon Laplace climate model; Boucher et al., 2020) and CESM2 (Community Earth System Model; Hurrell et al., 2013), two ocean-atmosphere coupled GCMs which took part in Climate Model Intercomparison Project (CMIP) Phase 6 (Eyring et al., 2016). We use predictions that start in 2008 and end in 2034 and follow the RCP8.5 (Representative Concentration Pathway) scenario, which tracks the observed CO₂ emissions closely (Schwalm et al., 2020). Predictions are provided as monthly grids with spatial resolutions of 1.27° × 2.5° on 79 vertical levels (IPSL-CM6) and 1.25° × 0.94° on 40 vertical levels (CESM). On these predictions of atmospheric conditions, we apply the COSP1.4 lidar simulator (Sect. 3.2), which generates on similar spatial grids the monthly averaged cloud properties that would be observed by a spaceborne lidar flying over the simulated atmosphere. In addition to the simulation steps described in Sect. 3.2, here, as a first step of the simulation, for each grid box of the GCM-created atmosphere an ensemble of subgrid-scale profiles are stochastically generated by the Subgrid Cloud Overlap Profile Sampler (Klein and Jakob, 1999). Each of these profiles is fed to the COSP simulator, which generates a synthetic lidar profile, on which cloud detection is performed. All subgrid-scale cloud detection profiles are eventually averaged to generate a single vertical profile for each grid box (see Chepfer et al., 2008, for details).

From the synthetic cloud properties, we considered two climate diagnostics whose trend should be related to climate change: first the fraction of opaque clouds $C_{opaque}$, defined as the number of lidar profiles in which an opaque cloud is detected in a given latitude-longitude grid box divided by the total number of profiles sampled in the same grid box. Opaque clouds are responsible for the majority of cloud radiative effects in the tropics (Vaillant de Guélis et al., 2017), and the cloud amount has been identified as one of the main drivers of cloud feedbacks on climate (Zelinka et al., 2016); thus the fraction of opaque clouds should be closely tied to climate change. Second, we considered the altitude of full attenuation $Z_{opaque$ (Guzman et al., 2017) averaged over all opaque profiles in every grid box. The vertical distribution of clouds is closely linked to their longwave radiative impact and to climate change (Vaillant de Guélis et al., 2018), and their altitude is expected to increase by several hundred meters per century (Richardson et al., 2022). Altitude is among the cloud properties whose change is expected to be detectable the earliest using active remote sensing (Chepfer et al., 2014; Takahashi et al., 2019; Aerenson et al., 2022).

From the GCM predictions, the COSP lidar simulator generates monthly grids of $C_{opaque}$ and $Z_{opaque$ that we spatially average over the tropics (30° S–30° N) to get monthly time series. We des seasonalize those time series to get their monthly anomalies over the 2008–2034 period. For any time
<table><caption>**Table 3.** Cloud detection statistics for CALIOP and ATLID in the case when the cloud is defined as SR (532 nm, z) > 5 and ATB (λ, z) – ATB<sub>mol</sub> (λ, z) > 2.5 × 10<sup>-6</sup> m<sup>-1</sup> sr<sup>-1</sup> (Eq. 6).</caption><thead><tr><th rowspan="3">Lidar</th><th colspan="8">CALIOP</th><th colspan="8">ATLID</th></tr><tr><th colspan="2">Day/night</th><th colspan="2">Night</th><th colspan="2">Day</th><th colspan="2">Night</th><th colspan="2">Day</th><th colspan="2">Night</th><th colspan="2">Day</th></tr><tr><th>Cloud</th><th>Ci</th><th>Y</th><th>Ci</th><th>Y</th><th>Ci</th><th>Y</th><th>Ci</th><th>Y</th><th>Ci</th><th>Y</th><th>Ci</th><th>Y</th><th>Ci</th><th>Y</th></tr></thead><tbody><tr><td>Averaged</td><td>N</td><td>Y</td><td>N</td><td>Y</td><td>N</td><td>Y</td><td>N</td><td>Y</td><td>N</td><td>Y</td><td>N</td><td>Y</td><td>N</td><td>Y</td></tr><tr><td>YESSi</td><td>7</td><td>6</td><td>7</td><td>8</td><td>7</td><td>6</td><td>9</td><td>8</td><td>7</td><td>6</td><td>7</td><td>7</td><td>6</td><td>6</td></tr><tr><td>NO_NO</td><td>82</td><td>77</td><td>86</td><td>86</td><td>72</td><td>73</td><td>79</td><td>85</td><td>82</td><td>80</td><td>86</td><td>86</td><td>79</td><td>77</td></tr><tr><td>YES_NO</td><td>11</td><td>17</td><td>6</td><td>7</td><td>15</td><td>20</td><td>5</td><td>6</td><td>10</td><td>14</td><td>7</td><td>7</td><td>15</td><td>18</td></tr><tr><td>NOSi</td><td>0</td><td>0</td><td>0</td><td>0</td><td>6</td><td>1</td><td>7</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Total score</td><td>88</td><td>81</td><td>93</td><td>93</td><td>67</td><td>75</td><td>85</td><td>93</td><td>89</td><td>85</td><td>93</td><td>93</td><td>83</td><td>80</td></tr><tr><td>Bias</td><td>-2</td><td>-12</td><td>-6</td><td>-7</td><td>13</td><td>-2</td><td>3</td><td>-5</td><td>-8</td><td>-9</td><td>-7</td><td>-7</td><td>-4</td><td>-10</td></tr><tr><td>RMS</td><td>14</td><td>11</td><td>6</td><td>6</td><td>22</td><td>20</td><td>7</td><td>5</td><td>9</td><td>9</td><td>6</td><td>6</td><td>18</td><td>17</td></tr></tbody></table>

t along these time series, the record length is equivalent to the period between 1 January 2008 and t, and we computed the trend $w(t)$ as the linear regression of the time series of anomalies over that period. The uncertainty $\sigma_w(t)$ in the trend $w(t)$ at time t was computed, as in Chepfer et al. (2018), as $\sigma_w(t) = \sigma_N \sqrt{\frac{1+\phi}{1-\phi}} n^{-\frac{3}{2}}$, with n being the number of years in the record at time t, $\phi$ being the lag-1 autocorrelation coefficient of the series between 0 and t, and $\sigma_N$ being the standard deviation of the noise remaining in the series between 0 and t once it has been deseasonalized and the autocorrelated part has been removed.

The following analysis focuses on the tropical regions (30° S–30° N), where the atmospheric circulation will be impacted by the weakening of the Hadley and Walker circulations expected in the upcoming century by most climate predictions (Davis and Rosenlof, 2012; Su et al., 2014; Kjellsson, 2015; Chemke, 2021). These changes will have important effects on the spatial distribution of tropical clouds (Su et al., 2014), which provide the basis for our climate diagnostics. Cloud opacity is one of the cloud properties most closely linked to their radiative impact (Zelinka et al., 2012), which explains why our diagnostics are based on the properties of opaque clouds (as in Perpina et al., 2021). The results below assume it will be possible to process ATLID measurements in such a way that CLIMP and GOCCP cloud properties are consistent.

## 4.3 How many years of ATLID observation are required in addition to CALIPSO to evaluate the climate model prediction of cloud changes?

Figure 12 shows how the uncertainty in the retrieved trend for $C_{opaque}$ changes with the length of the record of lidar-based cloud properties, starting in 2008, according to predictions from IPSL-CM6 (blue) and CESM1 (orange). The uncertainty is generally the largest and fluctuates most when the record is short and decreases and stabilizes as the record gets longer. At any time t if we require a 95% confidence

<figure><img src="image_11.png" /><figcaption>**Figure 11.** Same as Fig. 10 but for ATLID.</figcaption></figure>

level in the prediction and assume trends are normally distributed, the real trend will lie in the $w(t) \pm 2\sigma_w(t)$ interval. The sign of the trend will be robust once $|\frac{w(t)}{\sigma_w(t)}| > 2$. This is when the uncertainty in the trend becomes small compared to the trend itself and marks the time of emergence of cloud change induced by anthropogenic warming. This occurs earlier for strong, stable trends and might never occur for very small trends or trends whose sign changes over time. Times of emergence in the $C_{opaque}$ time series are indicated in Fig. 12 with triangles for three confidence levels (50%, 70%, and 95%). Reaching a reliable sign requires a longer record if the required confidence level is strong.

According to predictions from IPSL-CM6 (blue), a reliable trend should emerge from the natural variability at a 50% to 70% confidence level between 2030 and 2032. In other words, IPSL-CM6 predicts that revealing a reliable long-term trend in the fraction of opaque clouds would require an uninterrupted spaceborne lidar record of 22 years, which would be achievable if EarthCARE operates for at least 7 years. Reaching 95% confidence levels on the retrieved trend would require extending the record beyond
<figure><img src="image_12.png" /></figure>

**Figure 12.** Evolution of the uncertainty in the $C_{opaque}$ trend as a function of the length of the spaceborne lidar record, according to atmospheric conditions predicted by IPSL-CM6 (blue; Laboratoire de Météorologie Dynamique, LMD) and CESM (orange) in the period between 2008 and 2034 following the RCP8.5 scenario. The first 2 years of the record (2008–2010) are considered in the analysis, but trend uncertainties during that period are very large and are masked in the figure to improve the legibility of later years. CALIPSO’s (CAL) planned end of operation (2023) is marked by a vertical blue line. Supposing EarthCARE (EC) begins operation right afterward, its nominal 3-year operation point is marked by a vertical purple line, and an optimistic 10-year operation point is marked by a vertical black line.

25 years, most probably through another spaceborne lidar mission further in time. CESM1, meanwhile, predicts that a reliable long-term trend in the fraction of opaque clouds (at similar confidence levels between 50% and 70%) would be reached between 2025 and 2027, requiring 2 to 4 years of EarthCARE operation. A highly reliable trend (95% confidence levels) would be detectable in 2029, after 6 years of EarthCARE operation. In summary, if a 50% confidence level is acceptable, detecting a reliable trend would either be possible within the EarthCARE nominal operation time frame (2 years after launch), according to CESM1, or require EarthCARE to operate 4 years beyond its planned lifetime, according to IPSL-CM6.

If we consider the $Z_{opaque}$ diagnostic (Fig. 13), the IPSL-CM6 model now predicts a trend will be detectable at high 95% confidence levels in 2024, i.e., 1 year into EarthCARE’s nominal operation period. Meanwhile, according to CESM1 predictions, detecting a reliable trend (even at a modest 50% confidence level) would require EarthCARE operating for 8 years, 5 years beyond its nominal operation time frame. This very fast detection of a reliable $Z_{opaque}$ trend predicted by IPSL-CM6 is consistent with how this model expects important and fast changes in the vertical distribution of opaque clouds in the tropics (Perpina et al., 2021).

We sum these results up in Table 4, which in addition provides similar record lengths to detect reliable trends when considering grid boxes dominated by either low or high clouds. Tropical low opaque clouds include sparse shallow cumulus (Konsta et al., 2012) and optically thicker stratocumulus along the western coasts of continents (Guzman et al., 2017), both confined to the boundary layer and most frequent in subsidence regions. By contrast, tropical high opaque clouds are more localized and strongly correlated

<figure><img src="image_13.png" /></figure>

**Figure 13.** Same as Fig. 12 but for the altitude of opacity $Z_{opaque}$ instead of $C_{opaque}$.

with deep convection. Since both kinds of clouds are driven by very different processes, it is not unreasonable to assume they will probably evolve differently in the upcoming century, which justifies their separate studies. In practice, we identified grid boxes dominated by low clouds as those where $Z_{opaque$ was below 3 km and high-cloud grid boxes as those where $Z_{opaque$ was above 3 km. The results, shown in Table 4, suggest that the nominal ATLID (EarthCARE) operation will be enough to validate or invalidate the trends in opaque tropical low clouds predicted by CESM. It will be possible to validate or invalidate other model-based cloud predictions only if EarthCARE performs beyond its nominal lifetime (which is not impossible, as CALIPSO demonstrated) or if measurements from a follow-up spaceborne lidar mission after ATLID are included in the cloud profile record. These results are consistent with the trends, uncertainties, and times of emergence found when conducting a relatively simpler comparison of HadGEM2-A (Hadley Centre Global Environment Model) predictions in current vs. +4 K conditions (Chepfer et al., 2018). Needless to say that the treatment of any follow-up mission (e.g., Atmosphere Observing System, AOS; Aeolus-2) will require the compensation for all the differences between the lidars like it is done in this work.

As stated up front, these results depend on rather strong hypotheses of perfect continuity and perfect intercalibration between the consecutive spaceborne lidars that provide the measurements from which the cloud properties are derived. Imperfect continuity would occur if, for instance, EarthCARE starts operation later than CALIOP stops. The missing years in the record would delay the detection of a reliable trend by at least the same time period (Chepfer et al., 2018). Perfect intercalibration supposes the effects of instrumental differences in technical specifications (wavelengths, pulse energy, field of view, etc.) and orbital characteristics (local time of overpass, altitude) on lidar measurements are reconciled somehow. For instance, ATLID operates at 355 nm and CALIOP operates at 532 nm, and the impact this has on measurements can be reconciled by converting the ATLID signal at 532 nm as done in the current study, but the costs of this conversion are not completely understood and will require re-examination when actual ATLID data are available. Imperfect intercalibration could lead to offsets in one spaceborne lidar’s record compared to the other and would increase the uncertainties in the retrieved trends. Increased delays be-
<table><caption>**Table 4.** When will a spaceborne lidar record starting in 2008 be long enough to enable a reliable detection (at 70% confidence level) of $C_{opaque}$ or $Z_{opaque}$ trends according to predictions from IPSL-CM6 or CESM? The required years of EarthCARE operation are shown in parentheses, supposing they begin in 2023. The monthly evolution of the trend uncertainties for low and high clouds are provided in Figs. B1 and B2 of Appendix B.</caption><thead><tr><th rowspan="2"></th><th colspan="2">C<sub>opaque</sub></th><th colspan="2">Z<sub>opaque</sub></th></tr><tr><th>IPSL-CM6</th><th>CESM</th><th>IPSL-CM6</th><th>CESM</th></tr></thead><tbody><tr><td>All clouds</td><td>2030 (7 years)</td><td>2027 (4 years)</td><td>2021</td><td>2034 (11 years)</td></tr><tr><td>Low clouds only (&lt; 3 km)</td><td>No trend</td><td>2024 (1 year)</td><td>No trend</td><td>2025 (2 years)</td></tr><tr><td>High clouds only (&gt; 3 km)</td><td>2027 (4 years)</td><td>2031 (8 years)</td><td>2018</td><td>No trend</td></tr></tbody></table>

tween the operation of both instruments would complicate their intercalibration. The different local times of overpass (01:30 and 13:30 local solar time, LST, for CALIPSO; 06:00 and 18:00 LST for EarthCARE) are also quite problematic, since each instrument will sample clouds at a different phase of their diurnal cycle (Noel et al., 2018; Chepfer et al., 2019; Feofilov and Stubenrauch, 2019). In particular, this will impact high clouds related to deep convection that exhibit a marked diurnal cycle. It is out of the scope of the present work to evaluate how this change could bias the retrieved long-term trends. The same applies to a follow-up lidar mission, which may or may not operate at the same orbit with the same local solar time of overpass and may or may not measure the depolarized backscatter.

Finally, the times of emergence presented here must not be understood as definite but as predictions by climate models. It is worth noting, for instance, that, according to predictions from IPSL-CM6, a reliable trend should already be readily detectable in the existing record of $Z_{opaque}$ that is today only built on CALIOP (CALIPSO) (Table 4). Such a trend has not been identified yet. This is consistent with the fact that in current climate conditions IPSL-CM6 overestimates the altitude of opaque clouds in tropical convective regions and brings them significantly higher (+2 km) near the end of 21st century (Perpina et al., 2021). Such rapid changes are not present in CESM predictions. These important model differences highlight the crucial need for continued long-term cloud lidar observations able to monitor the actual cloud changes and disambiguate model predictions.

## 5 Conclusions

This study presents the physical basis for the ATLID cloud climate product named CLIMP. This product builds on previous work on CALIPSO, a space lidar dedicated to cloud and aerosol observations like ATLID. CALIPSO data have been used for 16 years to evaluate the description of clouds in climate models using a dedicated product named GOCCP and a dedicated lidar simulator named COSP lidar. The present work also builds on recent work on AEOLUS, a space lidar with HSRL capability operating in the UV, like ATLID.

Based on this legacy, we have defined the CLIMP-ST (Short-Term) and CLIMP-LT (Long-Term) products, both dedicated to cloud climate studies. Both contain the same variables as GOCCP (see Table D1 in Appendix D) on the same horizontal and vertical resolutions, but CLIMP-ST and CLIMP-LT have different cloud detection thresholds because they aim to tackle slightly different science objectives.

The CLIMP-ST product is designed to make full use of the ATLID capability for evaluating cloud description in climate models. CLIMP-ST is expected to contain optically thin cloud detected in daytime conditions at full resolution that were not observed by former space lidars at such high spatial resolutions during daytime. This new information, if confirmed in actual data, will help make progress regarding our current understanding of processes tied to thin ice clouds in the climate system. It will help evaluate the description in climate models of optically thin clouds in regions where they are frequent and important for climate, for example in the tropics and polar regions.

The CLIMP-LT product is designed to detect the same clouds as CALIPSO-GOCCP. Merging CLIMP-LT with GOCCP will allow for building a multi-decadal cloud profile record, useful to monitor the cloud inter-annual natural variability and cloud changes induced by human-caused climate warming. This record, if quality is sufficient, will be useful for evaluating climate prediction of cloud changes and for helping reduce uncertainties in model-based climate feedbacks and climate sensitivity.

To design CLIMP-ST and CLIMP-LT, we examined the differences between CALIOP and ATLID, space lidars that operate at different wavelengths and use different observation techniques and detectors. We sought to answer two questions. (1) Can the HSRL capability of ATLID help reconcile its cloud retrievals with the CALIOP record? (2) Does the cloud product retrieved from ATLID observations compare well with the one retrieved from CALIOP observations, and if so, how many years of ATLID observations are needed to detect trends in opaque cloud cover or altitude of opaque clouds, assuming ATLID operation will follow CALIOP without a gap?

To answer these questions, we coupled the outputs of the 3DCLOUD model with the COSP2 simulator and added in-
strumental noise for two cloud scenes, thin cirrus clouds at ~15 km in the tropics and stratocumulus clouds at ~1 km height. CALIOP and ATLID orbits over these cloud scenes were simulated both for nighttime and daytime conditions, at full vertical and horizontal (1/3 km) resolution and at 1 km horizontal resolution. Then, we applied a wavelength conversion algorithm to ATLID observations to convert UV lidar profiles into 532 nm lidar profiles and added synthetic noise generated for each instrument in accordance with its characteristics.

We addressed the first question for CLIMP-ST. We showed that the lower daytime noise of ATLID allows for applying more sensitive thresholds for cloud detection (SR (532 nm, z) > 3; ATB (532 nm, z) -ATB<sub>mol</sub> (532 nm, z) > 1.5 × 10<sup>-6</sup> m<sup>-1</sup> sr<sup>-1</sup>) than for CALIPSO at full spatial resolution in daytime without introducing a bias. This suggests that ATLID may provide new information on optically thin clouds at daytime conditions at full spatial resolution.

We addressed the second question for CLIMP-LT. We search for consistency between ATLID and CALIPSO-GOCCP in cloud detection; therefore we applied the same cloud detection threshold (SR (532 nm, z) > 5; ATB (532 nm, z) -ATB<sub>mol</sub> (532 nm, z) > 2.5 × 10<sup>-6</sup> m<sup>-1</sup> sr<sup>-1</sup>) to both instruments, then their nighttime cloud products are comparable, whereas the daytime CALIOP clouds are characterized by a somewhat higher false-detection rate. This suggests ATLID and CALIPSO might observe the same clouds, with some adjustment in the cloud detection scheme. Then we analyzed 24 years of predictions from two general circulation models (IPSL-CM6 and CESM2) in the RCP8.5 scenario, coupled with the COSP lidar simulator. We show that IPSL-CM6 predicts the opaque-cloud-cover trend detection will require 7 years of ATLID operation besides the existing CALIOP cloud dataset, whereas CESM2 predicts the opaque-cloud-cover trend can be detected in 4 years. For the clouds above 3 km altitude, these numbers change to 4 and 8 years, respectively, and for the altitudes below 3 km the IPSL-CM6 clouds indicate no trend and CESM cloud trend detection will require 1 year of ATLID operation. These differences in climate predictions highlight the need for a multi-decadal cloud lidar record.

The current results rely on a comparison of exactly the same atmospheric scenes “virtually observed” by two space lidars, and they were obtained in the framework of comparing the cloud detection capabilities of these two instruments. However, the comparison of the actual ATLID measurements with actual CALIOP ones will face with an uncompensated difference linked to the sampling of local solar time by CALIOP and ATLID. The difference in the diurnal cycle will bias the detected cloud amount and height. This is a separate issue that should be compensated for, and this should be a subject of a separate work. Moreover, the comparison of actual ATLID measurements with CALIOP ones will probably face unexpected differences other than the ones foreseen in this paper. Therefore, the CLIMP algorithm will require an adjustment after ATLID launch to take those into account.

That being said, this study suggests that it is likely that ATLID will provide new information useful to help evaluate cloud description in climate models beyond the existing space lidar observations. Moreover, merging the ATLID data with the CALIOP data will probably provide important information on cloud response to climate warming.
Appendix A

<figure><img src="image_14.png" /></figure>

**Figure A1.** Scattering ratio and cloud detection statistics estimated for stratocumulus clouds observed by CALIOP using Eq. (36): (a) scattering ratio at night, (b) cloud detection at night, (c) scattering ratio during the day, and (d) cloud detection during the day.

<figure><img src="image_15.png" /></figure>

**Figure A2.** Same as Fig. A1 but for ATLID.

<figure><img src="image_16.png" /></figure>

**Figure A3.** Scattering ratio and cloud detection statistics estimated for stratocumulus clouds observed by CALIOP using Eq. (6): (a) scattering ratio at night, (b) cloud detection at night, (c) scattering ratio during the day, and (d) cloud detection during the day.

<figure><img src="image_17.png" /></figure>

**Figure A4.** Same as Fig. A3 but for ATLID.
Appendix B

<figure><img src="image_18.png" /><figcaption>**Figure B1.** Same as Fig. 12 but with a separate analysis of the altitude of attenuation (ZFA) of (a) high-level clouds and (b) low-level clouds.</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>**Figure B2.** Same as Fig. 13 but with a separate analysis of (a) high-level clouds and (b) low-level clouds.</figcaption></figure>
Appendix C

<table><caption>**Table C1.** Variable definitions for the ATLID cloud product. Variables are of the type real (float64) unless specified otherwise. Shaded variables are used as dimensions. L1b: Level 1b.</caption><thead><tr><th>Variable name</th><th>Unit</th><th>Dimension</th><th>Remarks</th></tr></thead><tbody><tr><td>Time_UTC</td><td>s</td><td>Ntime</td><td>Same unit as in the ATLID L1b file</td></tr><tr><td>Altitude</td><td>m</td><td>NZ</td><td></td></tr><tr><td>Levels</td><td>m</td><td>Nlev (4)</td><td></td></tr><tr><td>Flags</td><td>Unitless</td><td>Nflags (6)</td><td></td></tr><tr><td>Lon</td><td>°</td><td>Ntime</td><td></td></tr><tr><td>Lat</td><td>°</td><td>Ntime</td><td></td></tr><tr><td>Surface_elevation</td><td>m</td><td>Ntime</td><td>From DEM and/or lidar ground return</td></tr><tr><td>Temperature</td><td>°C</td><td>Ntime × NZ</td><td>From ECMWF in ATLID L1b</td></tr><tr><td>Pressure</td><td>hPa</td><td>Ntime × NZ</td><td>From ECMWF in ATLID L1b</td></tr><tr><td>Scattering_ratio</td><td>Unitless</td><td>Ntime × NZ</td><td></td></tr><tr><td>Layer_identification_mask</td><td>Unitless (int8)</td><td>Ntime × NZ</td><td>See Table D2</td></tr><tr><td>Quality_flags</td><td>0/1 (int8)</td><td>Ntime × NZ × Nflags</td><td>See Table D3</td></tr><tr><td>Cloud_presence</td><td>0/1 (int8)</td><td>Ntime × Nlev</td><td>Nlev cloud flag at specific vertical levels<br>Nlev = 0 – anywhere in the profile<br>Nlev = 1 – at low levels<br>Nlev = 2 – at mid levels<br>Nlev = 3 – at high levels</td></tr></tbody></table>

<table><caption>**Table C2.** Layer identification mask description.</caption><thead><tr><th>Bin</th><th>Corresponding SR values</th></tr></thead><tbody><tr><td>0</td><td>Fully attenuated region:<br>SR &lt; SR_bins [0] (default 0.01)</td></tr><tr><td>1</td><td>Clear-sky region: SR_bins [0] &lt;<br>SR &lt; SR_bins [1] (default 1.2)</td></tr><tr><td>2</td><td>Unclassified region: SR_bins [1] &lt;<br>SR &lt; SR_bins [2] (default 3.0)</td></tr><tr><td>3 to 11</td><td>Cloud region: SR &gt; SR_bins [2]<br>The actual bin number provides information on<br>SR intensity within the cloud, with 3 being the weakest<br>signal and 11 being the strongest signal. Defaults are as follows.</td></tr><tr><td>3</td><td>4</td><td>5</td><td>6</td><td>7</td><td>8</td><td>9</td><td>10</td><td>11</td></tr><tr><td>5</td><td>7</td><td>10</td><td>15</td><td>20</td><td>25</td><td>30</td><td>40</td><td>50</td></tr></tbody></table>
<table><caption>**Table C3.** Quality flag indicator.</caption><thead><tr><td>Flag value</td><td>Explanation</td></tr></thead><tbody><tr><td>0</td><td>Missing or unreliable data, according to crosstalk information from ATLID Level 1b; if Mie, Rayleigh, geolocalization, or atmospheric quality are not good enough, the profile will be rejected and be considered missing or unreliable</td></tr><tr><td>1</td><td>Data located below the surface elevation</td></tr><tr><td>2</td><td>Noisy data, according to molecular calibration; if the calibration <i>R</i> is not within range, the entire profile is flagged as noisy</td></tr><tr><td>3</td><td>Conflicting cloud detection indicators in the upper troposphere: SR &lt; 3 and ΔATB &gt; 1.5 × 10<sup>-6</sup> m<sup>-1</sup> sr<sup>-1</sup></td></tr><tr><td>4</td><td>Presence of very bright clouds (SR &gt; 50) anywhere in the profile</td></tr><tr><td>5</td><td>Negative SR (SR &lt; 0); can appear in fully attenuated cloud mask (SR &lt; 0.01)</td></tr></tbody></table>
*Code and data availability.* We do not present the code in this work. The "lidar-consistent opaque cloud fractions and altitude of attenuation" dataset used in this work is available for download at https://doi.org/10.5281/zenodo.8059374 (Noel et al., 2023).

*Author contributions.* HC, VN, and AF: conceptualization, investigation, methodology, and validation. FS: calculation of cloud distributions. AF: data curation and formal analysis. AF: writing of original draft. AF, HC, and VN: review and editing.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

*Special issue statement.* This article is part of the special issue "EarthCARE Level 2 algorithms and data products". It is not associated with a conference.

*Acknowledgements.* This work was supported by the European Space Agency (ESA) (contract no. 4000xxxxxx/221/I-DT-lr). This work was also supported by the Centre National de la Recherche Scientifique (CNRS) and by the Centre National d'Études Spatiales (CNES) through the Expecting EarthCARE, Learning from A-Train (EECLAT) project as well as by the Office National d'Études et de Recherches Aérospatiales (ONERA, the French Aerospace Lab). The authors want to thank Institut Pascal at Université Paris-Saclay and the organizers of the 4-week Atmospheric Cycle Workshop for inspiring this study (French investissements d'avenir grant no. ANR-11-IDEX-0003-01). Special thanks go to Valery Shcherbakov (Laboratoire de Météorologie Physique, Université Clermont Auvergne, Aubière, France) for discussing the multiple scattering coefficient in relation to CALIOP and ATLID. Lastly, the authors express their gratitude to the reviewers, Sebastian Bley (ESA ESRIN, LITR) and Mark Vaughan (NASA Langley), for their thorough reviews and valuable comments. The authors would like to extend their sincere appreciation to the late Tobias Wehr, who served as the EarthCARE mission scientist at ESA, for his unwavering support of their work in the course of last 10 years.

*Financial support.* This work was supported by the European Space Agency (ESA) (contract no. 4000xxxxxx/221/I-DT-lr). This work was also supported by the Centre National de la Recherche Scientifique (CNRS) and by the Centre National d'Études Spatiales (CNES) through the Expecting EarthCARE, Learning from A-Train (EECLAT) project as well as by the Office National d'Études et de Recherches Aérospatiales (ONERA, the French Aerospace Lab).

*Review statement.* This paper was edited by Hajime Okamoto and reviewed by Mark Vaughan and Sebastian Bley.

## References

Aerenson, T., Marchand, R., Chepfer, H., and Medeiros, B.: When Will MISR Detect Rising High Clouds? *J. Geophys. Res.-Atmos.*, 127, e2021JD035865, https://doi.org/10.1029/2021JD035865, 2022.

Alkasem A., Szczap, F., Cornet, C., Shcherbakov, V., Gour, Y., Jourdan, O., Labonnnote, L. C., and Mioche, G.: Effects of cirrus heterogeneity on lidar CALIOP/CALIPSO data, JQSRT, 202, 38–49, https://doi.org/10.1016/j.jqsrt.2017.07.005, 2017.

Berry, E., Mace, G. G., and Gettelman, A.: Using A-Train Observations to Evaluate Cloud Occurrence and Radiative Effects in the Community Atmosphere Model during the Southeast Asia Summer Monsoon, J. Climate, 32, 4145–4165, https://doi.org/10.1175/JCLI-D-18-0693.1, 2019.

Beyerle, G., Gross, M. R., Haner, D. A., Kjome, N. T., McDermid, I. S., McGee, T. J., Rosen, J. M., Schäfer, H.-J., and Schrems, O.: A Lidar and Backscatter Sonde Measurement Campaign at Table Mountain during February–March 1997: Observations of Cirrus Clouds. J. Appl. Meteor., 40, 1275–1287, https://doi.org/10.1175/1520-0469(2001)058<1275:ALABSM>2.0.CO;2, 2001.

Bodas-Salcedo, A., Webb, M. J., Bony, S., Chepfer, H., J.-l. Dufresne, J.-L., Klein, S. A., Zhang, Y., Marchand, R., Haynes, J. M., Pincus, R., and John, V. O.: COSP, Satellite simulation software for model assessment, B. Am. Meteorol. Soc., 92, 1023–1043, https://doi.org/10.1175/2011BAMS2856.1, 2011.

Boucher, O., Servonnat, J., Albright, A. L., Aumont, O., Balkanski, Y., Bastrikov, V., Bekki, S., Bonnet, R., Bony, S., Bopp, L., Braconnot, P., Brockmann, P., Cadule, P., Caubel, A., Cheruy, F., Codron, F., Cozic, A., Cugnet, D., D'Andrea, F., Davini, P., de Lavergne, C., Denvil, S., Deshayes, J., Devilliers, M., Ducharne, A., Dufresne, J. L., Dupont, E., Éthé, C., Fairhead, L., Falletti, L., Flavoni, S., Fojols, M. A., Gardoll, S., Gastineau, G., Ghattas, J., Grandpeix, J. Y., Guenet, B., Guez, L. E., Guilyardi, E., Guimberteau, M., Hauglustaine, D., Hourdin, F., Idelkadi, A., Joussaume, S., Kageyama, M., Khodri, M., Krinner, G., Lebas, N., Levavasseur, G., Lévy, C., Li, L., Lott, F., Lurton, T., Luyssaert, S., Madec, G., Madeleine, J. B., Maignan, F., Marchand, M., Marti, O., Mellul, L., Meurdesoif, Y., Mignot, J., Musat, I., Ottlé, C., Peylin, P., Planton, Y., Polcher, J., Rio, C., Rochetin, N., Rousset, C., Sepulchre, P., Sima, A., Swingedouw, D., Thiéblemont, R., Traore, A. K., Vancoppenolle, M., Vial, J., Vialard, J., Viovy, N., and Vuichard, N.: Presentation and Evaluation of the IPSL-CM6A-LR Climate Model, J. Adv. Model. Earth Syst., 12, e2019MS002010, https://doi.org/10.1029/2019MS002010, 2020.

Cesana, G., Del Genio, A. D., and Chepfer, H.: The Cumulus And Stratocumulus CloudSat-CALIPSO Dataset (CASCAD), Earth Syst. Sci. Data, 11, 1745–1764, https://doi.org/10.5194/essd-11-1745-2019, 2019.

Chemke, R.: Future changes in the Hadley circulation: The role of ocean heat transport, Geophys. Res. Lett., 48, e2020GL091372. https://doi.org/10.1029/2020GL091372, 2021.
Cheng, Z., Liu, D., Yang, Y., Yang, L., and Huang, H.: Interferometric filters for spectral discrimination in high-spectral-resolution lidar: performance comparisons between Fabry-Perot interferometer and field-widened Michelson interferometer, Appl. Opt., 52, 7838–7850, https://doi.org/10.1364/AO.52.007838, 2013.

Chepfer, H., Bony, S., Winker, D., Chiriaco, M., Dufresne, J.-L., and Sèze, G.: Use of CALIPSO lidar observations to evaluate the cloudiness simulated by a climate model, Geophys. Res. Lett., 35, L15704, https://doi.org/10.1029/2008GL034207, 2008.

Chepfer, H., Bony, S., Winker, D., Cesana, G., Dufresne, J.-L., Minnis, P., Stubenrauch, C. J., and Zeng, S.: The GCM Oriented Calipso Cloud Product (CALIPSO-GOCCP), J. Geophys. Res., 115, D00H16, https://doi.org/10.1029/2009JD012251, 2010.

Chepfer, H., Cesana, G., Winker, D., Getzewich, B., and Vaughan, M.: Comparison of two different cloud climatologies derived from CALIOP Level 1 observations: the CALIPSO-ST and the CALIPSO-GOCCP, J. Atmos. Ocean. Technol., 30, 725–744 https://doi.org/10.1175/JTECH-D-12-00057.1, 2013.

Chepfer, H., Noël, V., Winker, D., and Chiriaco, M.: Where and when will we observe cloud changes due to climate warming?, Geophys. Res. Lett., 41, 8387–8395, https://doi.org/10.1002/2014GL061792, 2014.

Chepfer, H., Noël, V., Chiriaco, M., Wielicki, B., Winker, D., Loeb, N., and Wood, R.: The potential of multi-decades space-born lidar to constrain cloud feedbacks, J. Geophys. Res.-Atmos., 123, 5433–5454, https://doi.org/10.1002/2017JD027742, 2018.

Chepfer, H., Brogniez, H., and Noël, V.: Diurnal variations of cloud and relative humidity profiles across the tropics, Sci. Rep., 9, 16045, https://doi.org/10.1038/s41598-019-52437-6, 2019.

Chiriaco, M., Vautard, R., Chepfer, H., Haeffelin, M., Dudhia, J., Wanherdrick, Y., Morille, Y., and Protat, A.: The Ability of MM5 to Simulate Ice Clouds: Systematic Comparison between Simulated and Measured Fluxes and Lidar/Radar Profiles at the SIRTA Atmospheric Observatory, Mon. Weather Rev., 134, 897–918, https://doi.org/10.1175/MWR3102.1, 2006.

Davis, S. M. and Rosenlof, K. H.: A multidiagnostic intercomparison of tropical-width time series using reanalyses and satellite observations, J. Climate, 25, 1061–1078, 2012.

do Carmo, J. P., Villele, de V., Wallace, K., Lefebvre, A., Ghose, K., Kanitz, T., Chassat, F., Corselle, B., Belhadj, T., and Bravetti, P.: ATmospheric LIDar(ATLID): Pre-launch testing and calibration of the European Space Agency instrument that will measure aerosols and thin clouds in the atmosphere, Atmosphere, 12, 76, https://doi.org/10.3390/atmos12010076, 2021.

Donovan, D. P.: The Expected Impact of Multiple Scattering on ATLID Signals, EPJ Conf., 119, 01006, https://doi.org/10.1051/epjconf/201611901006, 2016.

Durand, Y., Hélière, A., Bézy, J.-L., and Meynart, R.: The ESA EarthCARE mission: Results of the ATLID instrument pre-developments, Proc. SPIE, 6750, 675015, https://doi.org/10.1117/12.737932, 2007.

Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958, https://doi.org/10.5194/gmd-9-1937-2016, 2016.

Feofilov, A. G. and Stubenrauch, C. J.: Diurnal variation of high-level clouds from the synergy of AIRS and IASI space-

borne infrared sounders, Atmos. Chem. Phys., 19, 13957–13972, https://doi.org/10.5194/acp-19-13957-2019, 2019.

Feofilov, A. G., Stubenrauch, C. J., and Delanoë, J.: Ice water content vertical profiles of high-level clouds: classification and impact on radiative fluxes, Atmos. Chem. Phys., 15, 12327–12344, https://doi.org/10.5194/acp-15-12327-2015, 2015.

Feofilov, A. G., Chepfer, H., Noël, V., Guzman, R., Gindre, C., Ma, P.-L., and Chiriaco, M.: Comparison of scattering ratio profiles retrieved from ALADIN/Aeolus and CALIOP/CALIPSO observations and preliminary estimates of cloud fraction profiles, Atmos. Meas. Tech., 15, 1055–1074, https://doi.org/10.5194/amt-15-1055-2022, 2022.

Garnier, A., Pelon, J., Vaughan, M. A., Winker, D. M., Trepte, C. R., and Dubuisson, P.: Lidar multiple scattering factors inferred from CALIPSO lidar and IIR retrievals of semi-transparent cirrus cloud optical depths over oceans, Atmos. Meas. Tech., 8, 2759–2774, https://doi.org/10.5194/amt-8-2759-2015, 2015.

Guzman, R., Chepfer, H., Noël, V., Vaillant de Guelis, T., Kay, J. E., Raberanto, P., Cesana, G., Vaughan, M. A., and Winker, D. M.: Direct atmosphere opacity observations from CALIPSO provide new constraints on cloud-radiation interactions, J. Geophys. Res.-Atmos., 122, 1066–1085, https://doi.org/10.1002/2016JD025946, 2017.

Haynes, J. M., Marchand, R. T., Luo, Z., Bodas-Salcedo, A., and Stephens, G. L.: A multipurpose radar simulation package: QuickBeam, B. Am. Meteorol. Soc., 88, 1723–1727, https://doi.org/10.1175/BAMS-88-11-1723, 2007.

Hélière, A., Gelsthorpe, R., Le Hors, L., and Toulemont, Y.: ATLID, the Atmospheric Lidar on board the EarthCARE Satellite, Proc. SPIE 10564, International Conference on Space Optics – ICSO 2012, 105642D (20 November 2017), https://doi.org/10.1117/12.2309095, 2017.

Heymsfield, A., Krämer, M., Wood, N. B., Gettelman, A., Field, P. R., and Liu, G.: Dependence of the Ice Water Content and Snowfall Rate on Temperature, Globally: Comparison of in Situ Observations, Satellite Active Remote Sensing Retrievals, and Global Climate Model Simulations, J. App. Meteorol. Climatol., 56, 189–215, https://doi.org/10.1175/JAMC-D-16-0230.1, 2017.

Heymsfield, A., Bansemer, A., Wood, N. B., Liu, G., Tanelli, S., Sy, O. O., Poellot, M., and Liu, C.: Toward Improving Ice Water Content and Snow-Rate Retrievals from Radars. Part II: Results from Three Wavelength Radar–Collocated In Situ Measurements and CloudSat–GPM–TRMM Radar Data, J. Appl. Met. Clim., 57, 365–389, https://doi.org/10.1175/JAMC-D-17-0164.1, 2018.

Hogan, R. J. and Illingworth, A. J.: Parameterizing Ice Cloud Inhomogeneity and the Overlap of Inhomogeneities Using Cloud Radar Data, J. Atmos. Sci., 60, 756–767, https://doi.org/10.1175/1520-0469(2003)060<0756:PICIAT>2.0.CO;2, 2003.

Hogan, R. J. and Kew, S.F.: A 3D stochastic cloud model for investigating the radiative properties of inhomogeneous cirrus clouds, Q. J. Roy. Meteorol. Soc., 131, 2585–2608, https://doi.org/10.1256/qj.04.144, 2005.

Hu, Y., Liu, Z., Winker, D., Vaughan, M., Noel, V., Bissonnette, L., Roy, G., and McGill, M.: Simple relation between lidar multiple scattering and depolarization for water clouds, Opt. Lett., 31, 1809–1811, 2006.

Hunt, W. H., Winker, D. M., Vaughan, M. A., Powell, K. A., Lucker, P. L., and Weimer, C.: CALIPSO Lidar Description and Perfor-
mance Assessment, J. Atmos. Ocean. Technol., 26, 1214–1228,
https://doi.org/10.1175/2009JTECHA1223.1, 2009.

Hurrell, J. W., Holland, M. M., Gent, P. R., Ghan, S., Kay, J.E., Kushner, P. J., Lamarque, J.-F., Large, W. G., Lawrence, D., Lindsay, K., Lipscomb, W. H., Long, M. C., Mahowald, N., Marsh, D. R., Neale, R. B., Rasch, P., Vavrus, S., Vertenstein, M., Bader, D., Collins, W. D., Hack, J. J., Kiehl, J., and Marshall, S.: The Community Earth System Model: A Framework for Collaborative Research, B. Am. Meteorol. Soc., 94, 1339–1360, 2013.

Illingworth, A. J., Barker, H. W., Beljaars, A. et al.: The Earth-
CARE satellite: the next step forward in global measurements
of clouds, aerosols, precipitation, and radiation, B. Am. Meteo-
rol. Soc., 96, 1311–1332, https://doi.org/10.1175/BAMS-D-12-
00227.1, 2015.

Kärcher, B., Thornberry, T. D., Krämer, M., and Jensen,
E. J.: On the statistical distribution of total water in
cirrus clouds, Geophys. Res. Lett., 45, 9963–9971,
https://doi.org/10.1029/2018GL079780, 2018.

Kingston, R. H.: Detection of Optical and Infrared Radiation, Springer Series in Optical Sciences, v10, Springer, https://doi.org/10.1007/978-3-540-35948-7, 142 pp., 1978.

Kjellsson, J.: Weakening of the global atmospheric circula-
tion with global warming, Clim. Dynam., 45, 975–988,
https://doi.org/10.1007/s00382-014-2337-8, 2015.

Klein, S. A. and Jakob, C.: Validation and Sensitivities of Frontal Clouds Simulated by the ECMWF Model, Mon. Weather Rev., 127, 2514–2531, https://doi.org/10.1175/1520-0493(1999)127<2514:VASOFC>2.0.CO;2, 1999.

Konsta, D., Chepfer, H., and Dufresne, J.-L.: A process ori-
ented characterization of tropical oceanic clouds for climate
model evaluation, based on a statistical analysis of day-
time A-Train observations, Clim. Dynam., 39, 2091–2108,
https://doi.org/10.1007/s00382012-1533-7, 2012.

Leahy, L. V., Wood, R., Charlson, R. J., Hostetler, C. A., Rogers, R. R., Vaughan, M. A., and Winker, D. M.: On the nature and extent of optically thin marine low clouds, J. Geophys. Res.-Atmos., 117, D22201, https://doi.org/10.1029/2012JD017929, 2022.

Liu, Z. and Sugimoto, N.: Simulation study for cloud detection with space lidars using analog detection photomultiplier tubes, Appl. Opt., 41, 1750–1759, 2002.

Liu, Z., Hunt, W. H., Vaughan, M., Hostetler, C., McGill, M., Pow-
ell, K., Winker, D., and Hu, Y.: Estimating random errors due
to shot noise in backscatter lidar observations, Appl. Opt., 45,
4437–4447, https://doi.org/10.1364/AO.45.004437, 2006.

Marchand, R., Haynes, J., Mace, G. G., Ackerman, T., and Stephens, G.: A comparison of simulated cloud radar output from the multiscale modeling framework global climate model with CloudSat cloud radar observations, J. Geophys. Res., 114, D00A20, https://doi.org/10.1029/2008JD009790, 2009.

McGill, M. J., Vaughan, M. A., Trepte, C. R., Hart, W. D., Hlavka, D. L., Winker, D. M., and Kuehn, R.: Air-borne validation of spatial properties measured by the CALIPSO lidar, J. Geophys. Res.-Atmos., 112, D20201, https://doi.org/10.1029/2007JD008768, 2007.

McGill, M. J., Yorks, J. E., Scott, V. S., Kupchock, A. W., and Selmer, P. A.: The Cloud-Aerosol Transport System (CATS): A technology demonstration on the International Space Station, Proc. Spie., 9612, https://doi.org/10.1117/12.2190841, 2015.

Nam C., Bony, S., Dufresne, J. L., and Chepfer, H.: The ‘too few, too bright’ tropical low-cloud problem in CMIP5 models, Geophys. Res. Lett., 39, 21, https://doi.org/10.1029/2012GL053421, 2012.

Noel, V.: Monthly global grids of lidar-consistent opaque cloud fractions and altitude of attenuation 2008–2034 according to predictions from CESM2 and IPSL-CM6 GCMs, Zenodo [data set], https://doi.org/10.5281/zenodo.8059374, 2023.

Noel, V., Chepfer, H., Chiriaco, M., and Yorks, J.: The diur-
nal cycle of cloud profiles over land and ocean between 51°
S and 51°N, seen by the CATS spaceborne lidar from the In-
ternational Space Station, Atmos. Chem. Phys., 18, 9457–9473,
https://doi.org/10.5194/acp-18-9457-2018, 2018.

Perpina, M., Noel, V., Chepfer, H., Guzman, R., and Feofilov, A. G.: Link between opaque cloud properties and atmospheric dynamics in observations and simulations of current climate in the Tropics, and impact on future predictions. J. Geophys. Res.-Atmos., 126, e2020JD033899. https://doi.org/10.1029/2020JD033899, 2021.

Platt, C. M. R.: Lidar and radiometric observations of cirrus clouds, J. Atmos. Sci., 30, 1191–1204, 1973.

Rasch, P. J., Xie, S., Ma, P.-L., Lin, W., Wang, H., Tang, Q., Burrows, S. M., Caldwell, P., Zhang, K., Easter, R. C., Cameron-Smith, P., Singh, B., Wan, H., Golaz, J.-C., Harrop, B. E., Roesler, E., Bacmeister, J., Larson, V. E., Evans, K. J., Qian, Y., Taylor, M., Leung, L. R., Zhang, Y., Brent, L., Branstetter, M., Hannay, C., Mahajan, S., Mametjanov, A., Neale, R., Richter, J. H., Yoon, J.-H., Zender, C. S., Bader, D., Flanner, M., Foucar, J. G., Jacob, R., Keen, N., Klein, S. A., Liu, X., Salinger, A. G., Shrivastava, M., and Yang, Y.: An Overview of the Atmospheric Component of the Energy Exascale Earth System Model, J. Adv. Model. Earth Syst., 11, 2377–2411, https://doi.org/10.1029/2019MS001629, 2019.

Reitebuch, O., Lemmerz, C., Lux, O., Marksteiner, U., Rahm, S., Weiler, F., Witschas, B., Meringer, M., Schmidt, K., Huber, D., Nikolaus, I., Geiss, A., Vaughan, M., Dabas, A., Flament, T., Stieglitz, H., Isaksen, L., Rennie, M., de Kloe, J., Marseille, G.-J., Stoffelen, A., Wernham, D., Kanitz, T., Straume, A.-G., Fehr, T., von Bismarck, J., Floberghagen, R., and Parrinello, T.: Initial Assessment of the Performance of the First Wind Lidar in Space on Aeolus, EPJ Web Conf., 237, 01010, https://doi.org/10.1051/epjconf/202023701010, 2020.

Reverdy M., Chepfer, H., Donovan, D., Noël, V., Cesana,
G., Hoareau, C., Chiriaco, M., and Bastin, S.: An Earth-
CARE/ATLID simulator to evaluate cloud description in cli-
mate models, J. Geophys. Res.-Atmos., 120, 11090–11113,
https://doi.org/10.1002/2015JD023919, 2015.

Richardson, M. T., Roy, R. J., and Lebsock, M. D.: Satel-
lites Suggest Rising Tropical High Cloud Altitude:
2002–2021, Geophys. Res. Lett., 49, e2022GL098160,
https://doi.org/10.1029/2022GL098160, 2022.

Rogers, R. R., Vaughan, M. A., Hostetler, C. A., Burton, S. P., Fer-
rare, R. A., Young, S. A., Hair, J. W., Obland, M. D., Harper, D.
B., Cook, A. L., and Winker, D. M.: Looking through the haze:
evaluating the CALIPSO level 2 aerosol optical depth using air-
borne high spectral resolution lidar data, Atmos. Meas. Tech.,
7, 4317–4340, https://doi.org/10.5194/amt-7-4317-2014, 2014.

Sassen, K. and Comstock, J. M.: A midlatitude cirrus
cloud climatology from the facility for atmospheric
remote sensing. Part III: Radiative properties, J. At-
mos. Sci., 58, 2113-2127, https://doi.org/10.1175/1520-0469(2001)058<2113:Amcccf>2.0.Co;2, 2001.

Schwalm, C. R., Glendon, S., and Duffy, P. B.: RCP85 tracks cumulative CO₂ emissions, P. Natl. Acad. Sci. USA, 117, 19656-19657, https://doi.org/10.1073/pnas.2007117117, 2020.

Shcherbakov, V., Szczap, F., Alkasem, A., Mioche, G., and Cornet, C.: Empirical model of multiple-scattering effect on single-wavelength lidar data of aerosols and clouds, Atmos. Meas. Tech., 15, 1729-1754, https://doi.org/10.5194/amt-15-1729-2022, 2022.

Stevens, B., H-Moeng, C., Ackerman A. S., Bretherton, C. S., Chlond, A., de Roode, S., Edwards, J., Golaz, J. C., Jiang, H., Khairoutdinov, M., Kirkpatrick, M. P., Lewellen, D. C., Lock, A. L., Müller, F., Stevens, D. E., Whelan, E., and Zhu, P.: Evaluation of large-eddy simulations via observations of nocturnal marine stratocumulus, Mon. Weather Rev., 133, 1443-1462, 2005.

Straume, A. G., Rennie, M., Isaksen, L., de Kloe, J., Marseille, G.-J., Stoffelen, A., Flament, T., Stieglitz, H., Dabas, A., Huber, D., Reitebuch, O., Lemmerz, C., Lux, O., Marksteiner, U., Weiler, F., Witschas, B., Meringer, M., Schmidt, K., Nikolaus, I., Geiss, A., Flamant, P., Kanitz, T., Wernham, D., von Bismarck, J., Bley, S., Fehr, T., Floberghagen, R., and Parrinello, T.: ESA's Space-Based Doppler Wind Lidar Mission Aeolus - First Wind and Aerosol Product Assessment Results, EPJ Web Conf. 237, 01007, https://doi.org/10.1051/epjconf/202023701007, 2020.

Su, H., Jiang, J. H., Zhai, C., Shen, T. J., Neelin, J. D., Stephens, G. L., and Yung, Y. L.: Weakening and strengthening structures in the Hadley Circulation change under global warming and implications for cloud response and climate sensitivity: Circulation, Clouds, Climate Sensitivity, J. Geophys. Res.-Atmos., 119, 5787-5805. https://doi.org/10.1002/2014JD021642, 2014.

Swales, D. J., Pincus, R., and Bodas-Salcedo, A.: The Cloud Feedback Model Intercomparison Project Observational Simulator Package: Version 2, Geosci. Model Dev., 11, 77–81, https://doi.org/10.5194/gmd-11-77-2018, 2018.

Szczap, F., Gour, Y., Fauchez, T., Cornet, C., Faure, T., Jourdan, O., Penide, G., and Dubuisson, P.: A flexible three-dimensional stratocumulus, cumulus and cirrus cloud generator (3DCLOUD) based on drastically simplified atmospheric equations and the Fourier transform framework, Geosci. Model Dev., 7, 1779–1801, https://doi.org/10.5194/gmd-7-1779-2014, 2014.

Takahashi, H., Lebsock, M. D., Richardson, M., Marchand, R., and Kay, J. E.: When Will Spaceborne Cloud Radar Detect Upward Shifts in Cloud Heights?, J. Geophys. Res.-Atmos., 124, 2018JD030242, https://doi.org/10.1029/2018JD030242, 2019.

Vaillant de Guélis, T., Chepfer, H., Noel, V., Guzman, R., Dubuisson, P., Winker, D. M., and Kato, S.: The link between outgoing longwave radiation and the altitude at which a spaceborne lidar beam is fully attenuated, Atmos. Meas. Tech., 10, 4659–4685, https://doi.org/10.5194/amt-10-4659-2017, 2017.

Vaillant de Guélis T., H. Chepfer, Noël V., Guzman R., Bonazola M., and Winker D. M.: Space lidar observations constrain longwave cloud feedback, Nature Sci. Rep., 8, 16570, https://doi.org/10.1038/s41598-018-34943-1, 2018.

Voudouri, K. A., Giannakaki, E., Komppula, M., and Balis, D.: Variability in cirrus cloud properties using a PollyXT Raman lidar over high and tropical latitudes, Atmos. Chem. Phys., 20, 4427–4444, https://doi.org/10.5194/acp-20-4427-2020, 2020.

Waliser, D., Gleckler, P. J., Ferraro, R., Taylor, K. E., Ames, S., Biard, J., Bosilovich, M. G., Brown, O., Chepfer, H., Cinquini, L., Durack, P. J., Eyring, V., Mathieu, P.-P., Lee, T., Pinnock, S., Potter, G. L., Rixen, M., Saunders, R., Schulz, J., Thépaut, J.-N., and Tuma, M.: Observations for Model Intercomparison Project (Obs4MIPs): status for CMIP6, Geosci. Model Dev., 13, 2945–2958, https://doi.org/10.5194/gmd-13-2945-2020, 2020.

Winker, D. M., Pelon, J., Coakley Jr., J. A., Ackerman, S. A., Charlson, R. J., Colarco, P. R., Flamant, P., Fu, Q., Hoff, R. M., Kittaka, C., Kubar, T. L., Le Treut, H., Mccormick, M. P., Mégie, G., Poole, L., Powell, K., Trepte, C., Vaughan, M. A., and Wielicki, B. A.: The CALIPSO Mission, A Global 3D View of Aerosols and Clouds, B. Am. Meteor. Soc., 91, 1211–1230, https://doi.org/10.1175/2010BAMS3009.1, 2010.

Winker, D., Chepfer, H., Noël V., and Cai, X.: Observational constraints on cloud feedbacks: The role of active satellite sensors, Surv. Geophys., 38, 1483–1508, https://doi.org/10.1007/s10712-017-9452-0, 2017.

Young, S. A., Vaughan, M. A., Garnier, A., Tackett, J. L., Lambeth, J. D., and Powell, K. A.: Extinction and optical depth retrievals for CALIPSO’s Version 4 data release, Atmos. Meas. Tech., 11, 5701–5727, https://doi.org/10.5194/amt-11-5701-2018, 2018.

Zelinka, M. D., Klein, S. A., and Hartmann, D. L.: Computing and Partitioning Cloud Feedbacks Using Cloud Property Histograms. Part II: Attribution to Changes in Cloud Amount, Altitude, and Optical Depth, J. Climate, 25, 3736–3754, https://doi.org/10.1175/JCLI-D-11-00249.1, 2012.

Zelinka, M. D., Zhou, C., and Klein, S. A.: Insights from a refined decomposition of cloud feedbacks, Geophys. Res. Lett., 43, 9259–9269, https://doi.org/10.1002/2016GL069917, 2016.

Zelinka, M. D., Myers, T. A., McCoy, D. T., Po-Chedley, S., Caldwell, P. M., Ceppi, P., Klein, S. A., and Taylor, K. E.: Causes of higher climate sensitivity in CMIP6 models, Geophys. Res. Lett., 47, e2019GL085782, https://doi.org/10.1029/2019GL085782, 2020.