<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "https://jats.nlm.nih.gov/nlm-dtd/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" specific-use="SMUR" dtd-version="3.0" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="publisher">EGUsphere</journal-id>
<journal-title-group>
<journal-title>EGUsphere</journal-title>
<abbrev-journal-title abbrev-type="publisher">EGUsphere</abbrev-journal-title>
<abbrev-journal-title abbrev-type="nlm-ta">EGUsphere</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub"></issn>
<publisher><publisher-name>Copernicus Publications</publisher-name>
<publisher-loc>Göttingen, Germany</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.5194/egusphere-2026-2295</article-id>
<title-group>
<article-title>On the origin of the twilight color index maximum and its application to cloud-height retrieval</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Toledo</surname>
<given-names>Daniel</given-names>
<ext-link>https://orcid.org/0000-0002-0103-1891</ext-link>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
</contrib-group><aff id="aff1">
<label>1</label>
<addr-line>Instituto Nacional de Técnica Aeroespacial INTA, 28850 Torrejon de Ardoz, Spain</addr-line>
</aff>
<pub-date pub-type="epub">
<day>26</day>
<month>05</month>
<year>2026</year>
</pub-date>
<volume>2026</volume>
<fpage>1</fpage>
<lpage>28</lpage>
<permissions>
<copyright-statement>Copyright: &#x000a9; 2026 Daniel Toledo</copyright-statement>
<copyright-year>2026</copyright-year>
<license license-type="open-access">
<license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri"  xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p>
</license>
</permissions>
<self-uri xlink:href="https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2295/">This article is available from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2295/</self-uri>
<self-uri xlink:href="https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2295/egusphere-2026-2295.pdf">The full text article is available as a PDF file from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2295/egusphere-2026-2295.pdf</self-uri>
<abstract>
<p>A number of previous studies have demonstrated the capability of detecting high-altitude clouds during twilight using the color index (CI), defined as the ratio of zenith intensities at two different wavelengths, typically selected in the visible range or near infrared (NIR). When high clouds are present, a maximum or minimum (depending on the wavelengths selection) is observed in the CI signal (Sarkissian et al., 1991; Toledo et al., 2016). These studies also showed that the solar zenith angle (SZA) at which the CI maximum or minimum occurs (SZA&lt;sub&gt;max&lt;/sub&gt;) strongly depends on cloud altitude, enabling cloud-height retrieval through comparison with radiative transfer (RT) simulations. Twilight conditions require RT simulations in spherical geometry, which are computationally expensive. In this work, we introduce a single-scattering formulation of the CI that provides a physically transparent framework for identifying the mechanisms that determine the SZA of the CI maximum and, consequently, the inferred cloud altitude. The simplified formulation is explicitly compared with Monte Carlo RT simulations in spherical geometry and is shown to accurately reproduce the behavior of SZA&lt;sub&gt;max&lt;/sub&gt; over a wide range of conditions relevant for high-altitude clouds. In particular, the model provides reliable cloud-height estimates for cloud optical depths up to &amp;tau;&lt;sub&gt;&lt;em&gt;C&lt;/em&gt;&lt;/sub&gt; ≲ 0.3. Within this single-scattering formulation, we demonstrate that SZA&lt;sub&gt;max&lt;/sub&gt; occurs at the SZA for which the relative SZA-variations of the zenith intensity at the two selected wavelengths become equal, thereby explaining the emergence of the extremum in differential terms. However, achieving a well-defined CI extremum requires selecting two wavelengths with sufficient spectral separation, typically spanning distinct regions of the visible--NIR spectrum. To overcome this spectral dependence, we introduce a Rayleigh-referenced color index (CI&lt;sup&gt;&lt;em&gt;R&lt;/em&gt;&lt;/sup&gt;), defined as the ratio between the measured zenith intensity and the corresponding intensity expected for a purely Rayleigh-scattering atmosphere at the same wavelength. This index reproduces the characteristic extrema associated with high-altitude clouds while requiring simulations and observations at only a single wavelength. The proposed formulation facilitates extensive sensitivity studies and provides greater flexibility in spectral selection, particularly in regions affected by gas absorption.</p>
</abstract>
<counts><page-count count="28"/></counts>
<funding-group>
<award-group id="gs1">
<funding-source>Ministerio de Ciencia e Innovación</funding-source>
<award-id>PID2022-139386OA-I00</award-id>
</award-group>
</funding-group>
</article-meta>
</front>
<body/>
<back>
</back>
</article>