<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.3 20210610//EN" "JATS-journalpublishing1-3-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="1.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Remote Sens.</journal-id>
<journal-title-group>
<journal-title>Frontiers in Remote Sensing</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Remote Sens.</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">2673-6187</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1640320</article-id>
<article-id pub-id-type="doi">10.3389/frsen.2025.1640320</article-id>
<article-version article-version-type="Version of Record" vocab="NISO-RP-8-2008"/>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Original Research</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Gaia&#x2019;s Crown: a deep space mirage seen from DSCOVR/EPIC during lunar transit</article-title>
<alt-title alt-title-type="left-running-head">Blank et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/frsen.2025.1640320">10.3389/frsen.2025.1640320</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Blank</surname>
<given-names>Karin</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1183950"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Methodology" vocab-term-identifier="https://credit.niso.org/contributor-roles/methodology/">Methodology</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Software" vocab-term-identifier="https://credit.niso.org/contributor-roles/software/">Software</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Visualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/visualization/">Visualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing &#x2013; review and editing</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Herman</surname>
<given-names>Jay</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1367750"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Validation" vocab-term-identifier="https://credit.niso.org/contributor-roles/validation/">Validation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing &#x2013; review and editing</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Dangelo</surname>
<given-names>Sarah</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Software" vocab-term-identifier="https://credit.niso.org/contributor-roles/software/">Software</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Validation" vocab-term-identifier="https://credit.niso.org/contributor-roles/validation/">Validation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing &#x2013; review and editing</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Marshak</surname>
<given-names>Alexander</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1121340"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing &#x2013; review and editing</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tennenbaum</surname>
<given-names>Andrew</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x26; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/">Writing &#x2013; review and editing</role>
</contrib>
</contrib-group>
<aff id="aff1">
<label>1</label>
<institution>Goddard Space Flight Center, NASA</institution>, <city>Greenbelt</city>, <state>MD</state>, <country country="US">United States</country>
</aff>
<aff id="aff2">
<label>2</label>
<institution>University of Maryland Baltimore County</institution>, <city>Baltimore</city>, <state>MD</state>, <country country="US">United States</country>
</aff>
<author-notes>
<corresp id="c001">
<label>&#x2a;</label>Correspondence: Karin Blank, <email xlink:href="mailto:karin.b.blank@nasa.gov">karin.b.blank@nasa.gov</email>
</corresp>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="2026-01-12">
<day>12</day>
<month>01</month>
<year>2026</year>
</pub-date>
<pub-date publication-format="electronic" date-type="collection">
<year>2025</year>
</pub-date>
<volume>6</volume>
<elocation-id>1640320</elocation-id>
<history>
<date date-type="received">
<day>03</day>
<month>06</month>
<year>2025</year>
</date>
<date date-type="rev-recd">
<day>18</day>
<month>11</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>24</day>
<month>11</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2026 Blank, Herman, Dangelo, Marshak and Tennenbaum.</copyright-statement>
<copyright-year>2026</copyright-year>
<copyright-holder>Blank, Herman, Dangelo, Marshak and Tennenbaum</copyright-holder>
<license>
<ali:license_ref start_date="2026-01-12">https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License (CC BY)</ext-link>. The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</license-p>
</license>
</permissions>
<abstract>
<p>The Earth Polychromatic Imaging Camera (EPIC), onboard the Deep Space Climate Observatory (DSCOVR) spacecraft, located at the Earth-Sun Lagrange 1 point, has captured a unique optical effect during lunar occultation named &#x201c;Gaia&#x2019;s Crown.&#x201d; In EPIC images, the phenomenon appears as a small &#x201c;flange&#x201d; at the Earth&#x2013;Moon contact when the Moon is roughly half below Earth&#x2019;s limb; it is present in the visible and near-infrared channels but absent in the ultraviolet. Using atmospheric data and 3D, voxel-based ray tracing models, this effect was identified as a combination of atmospheric distortion and a complex mirage caused by variations in the Earth&#x2019;s atmosphere. Additionally, it is shown that while satellites closer to the Earth can see a similar phenomenon, Gaia&#x2019;s Crown presents unique distortion effects that demonstrate how EPIC&#x2019;s vantage point at 1.5 million kilometers from Earth provides a different perspective on atmospheric optics.</p>
</abstract>
<kwd-group>
<kwd>DSCOVR</kwd>
<kwd>EPIC</kwd>
<kwd>atmospheric optics</kwd>
<kwd>mirage</kwd>
<kwd>GOES (geostationary operational environmental satellite)</kwd>
<kwd>Himawari-8</kwd>
</kwd-group>
<funding-group>
<funding-statement>The authors declare that financial support was received for the research and/or publication of this article. The work utilized the DSCOVR/EPIC data, which utilizes computational resources and services provided by the NASA Center for Climate Simulation (NCCS) at the Goddard Space Flight Center.</funding-statement>
</funding-group>
<counts>
<fig-count count="15"/>
<table-count count="3"/>
<equation-count count="19"/>
<ref-count count="28"/>
<page-count count="16"/>
</counts>
<custom-meta-group>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Multi- and Hyper-Spectral Imaging</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<label>1</label>
<title>Introduction</title>
<p>The Deep Space Climate Observatory (DSCOVR) is a NOAA-operated satellite that orbits the Earth-Sun Lagrange-1 (L<sub>1</sub>) point. This location, at approximately 1.5 million kilometers away from the Earth, offers a unique view of our planet where the disk is almost fully illuminated. The DSCOVR spacecraft serves as a host to both Sun- and Earth-observing instruments.</p>
<p>The Earth Polychromatic Imaging Camera (EPIC), a NASA instrument onboard DSCOVR, utilizes a 30&#xa0;cm Cassegrain telescope with a 0.62&#xb0; field of view to photograph the Earth from this vantage point. Equipped with a 2048 &#xd7; 2048 charge-coupled device (CCD) sensor, it takes 13 or 22 image sets daily. Each set consists of 10 narrowband images covering wavelengths between 317 and 780&#xa0;nm, with the more frequent sets occurring during the Northern Hemisphere summer (<xref ref-type="bibr" rid="B14">Marshak et al., 2018</xref>). This wavelength range permits a unique view of the Earth in ultraviolet, visible, and near infrared, which enables detection of ozone, sulfur dioxide, aerosols, vegetation, ocean, cloud properties, and other science applications (<xref ref-type="bibr" rid="B14">Marshak et al., 2018</xref>).</p>
<p>This brief report documents an observational phenomenon in EPIC imagery and demonstrates a first-order, physically plausible mechanism. We do not attempt a full radiative-transfer or global occurrence analysis. Our aim is identification and plausibility, leaving comprehensive modeling for follow-up work.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Observations</title>
<p>The EPIC camera nominally observes just the sunlit Earth, but its 0.62&#xb0; field of view allows it to occasionally capture the Moon as well, as seen in the left portion of <xref ref-type="fig" rid="F1">Figure 1</xref>. The instrument observes these events roughly four times a year - two events with the Moon passing in front of the Earth and the remainder with it passing behind.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Left: Color image taken by EPIC of the Moon passing behind the Earth taken on 27 September 2015. Right: A close-up of the 680&#xa0;nm band that shows the lunar limb phenomenon where the Moon intersects with the Earth. Image has been gamma corrected to improve the visibiliy of the phenomenon.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g001.tif">
<alt-text content-type="machine-generated">A composite image shows the Earth with cloud patterns and land masses visible, and the Moon partially peeking from behind. On the right, a grayscale, pixelated close-up possibly depicts an element of Earth&#x27;s surface or atmosphere against a dark background.</alt-text>
</graphic>
</fig>
<p>In 2015, a special imaging was performed where the EPIC camera was trained upon the Moon, following it as it transited behind the Earth. Upon more detailed inspection of the images, an interesting effect was seen as the Moon passed behind the Earth. In this sequence, as seen in the right, zoomed-in portion of <xref ref-type="fig" rid="F1">Figure 1</xref>, a small &#x201c;flange&#x201d; is seen, appearing as almost a mount for the Moon against the planet&#x2019;s surface. It was not known what this phenomenon was or how it originated. This paper will discuss the origin of the phenomenon and the underlying optics.</p>
<p>A review of the original, unprocessed level 0 data, which contains the raw counts as received from the instrument, also revealed the feature, eliminating possible level 1a processing calibration artifacts as the cause. As part of this investigation, we determined that pixels in the feature are not saturated, indicating that the cause is not CCD bloom, nor does it mimic any known optical artifact of the telescope including stray light (<xref ref-type="bibr" rid="B4">Cede et al., 2021</xref>). This was true for all images investigated as part of this paper.</p>
<p>Because the sequence was part of a public outreach imaging to acquire color photos of the Moon passing behind the Earth&#x2019;s disk, the instrument recorded only the red (680&#xa0;nm), green (551&#xa0;nm), and blue (443&#xa0;nm) channels. The phenomenon was seen in all 3 wavelengths.</p>
<p>A review of the EPIC archive was conducted and all images where the Moon passed behind the Earth were obtained. Out of the 16 image sets obtained where the Moon intersected the Earth&#x2019;s horizon, only two additional sets contained the phenomenon. In the set of data collected over 10&#xa0;years, it occurred only 16.7% of the time. As seen in <xref ref-type="table" rid="T1">Table 1</xref>, the only common elements between times it is seen is that it occurred when the limb was dark, and the Moon was halfway below the horizon. The orientation of the Moon on the CCD, or the physical location the lunar limb intersects with the Earth appear to be independent. The phenomenon is not seen when the limb is illuminated where the Moon intersects.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Table contains observed metadata regarding the EPIC observations collected.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Date/Time UTC</th>
<th align="left">Phenomenon visible?</th>
<th align="left">Moon intersection: Earth&#x2019;s limb lighting conditions</th>
<th align="left">Imaging type</th>
<th align="left">Moon orientation on CCD</th>
<th align="left">Moon intersection: Location on earth</th>
<th align="left">Moon more than halfway above horizon</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">2015-09-27 10:35</td>
<td align="left">Yes</td>
<td align="left">Dark</td>
<td align="left">443, 551, 680</td>
<td align="left">Center</td>
<td align="left">Pacific</td>
<td align="left">No</td>
</tr>
<tr>
<td align="left">2015-11-25 20:29</td>
<td align="left">No</td>
<td align="left">Dark</td>
<td align="left">All 10 bands</td>
<td align="left">Top</td>
<td align="left">Pacific</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2016-03-22 19:55</td>
<td align="left">No</td>
<td align="left">Dark</td>
<td align="left">All 10 bands</td>
<td align="left">Top</td>
<td align="left">Pacific</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2016-09-16 04:38</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">All 10 bands</td>
<td align="left">Bottom</td>
<td align="left">South America</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2017-03-11 18:54</td>
<td align="left">No</td>
<td align="left">Dark</td>
<td align="left">All 10 bands</td>
<td align="left">Left</td>
<td align="left">Antarctica</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2017-10-05 00:45</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">All 10 bands</td>
<td align="left">Bottom</td>
<td align="left">Pacific</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2018-11-24 02:42</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">All 10 bands</td>
<td align="left">Bottom</td>
<td align="left">Pacific</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2019-06-17 17:16</td>
<td align="left">No</td>
<td align="left">Dark</td>
<td align="left">All 10 bands</td>
<td align="left">Top</td>
<td align="left">Africa</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2019-06-17 21:38</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">All 10 bands</td>
<td align="left">Bottom</td>
<td align="left">Pacific</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2020-10-02 06:31</td>
<td align="left">No</td>
<td align="left">Dark</td>
<td align="left">443, 551, 680</td>
<td align="left">Top</td>
<td align="left">Pacific</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2020-10-02 07:01</td>
<td align="left">Yes</td>
<td align="left">Dark</td>
<td align="left">443, 551, 680</td>
<td align="left">Top</td>
<td align="left">Pacific</td>
<td align="left">No</td>
</tr>
<tr>
<td align="left">2020-10-02 11:01</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">443, 551, 680</td>
<td align="left">Bottom left</td>
<td align="left">Atlantic</td>
<td align="left">No</td>
</tr>
<tr>
<td align="left">2020-10-02 11:31</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">443, 551, 680</td>
<td align="left">Bottom left</td>
<td align="left">North America</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2020-11-30 17:02</td>
<td align="left">Yes</td>
<td align="left">Dark</td>
<td align="left">All 10 bands</td>
<td align="left">Top right</td>
<td align="left">Atlantic</td>
<td align="left">No</td>
</tr>
<tr>
<td align="left">2021-12-18 09:41</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">All 10 bands</td>
<td align="left">Bottom left</td>
<td align="left">Pacific</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2021-12-18 15:05</td>
<td align="left">No</td>
<td align="left">Dark</td>
<td align="left">All 10 bands</td>
<td align="left">Top</td>
<td align="left">North America</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2022-02-16 12:49</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">All 10 bands</td>
<td align="left">Bottom left</td>
<td align="left">Atlantic</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2022-02-16 16:25</td>
<td align="left">No</td>
<td align="left">Illuminated</td>
<td align="left">All 10 bands</td>
<td align="left">Bottom left</td>
<td align="left">Atlantic</td>
<td align="left">Yes</td>
</tr>
<tr>
<td align="left">2023-02-05 20:48</td>
<td align="left">No</td>
<td align="left">Dark</td>
<td align="left">All 10 bands</td>
<td align="left">Top right</td>
<td align="left">Antarctica</td>
<td align="left">Yes</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Definitions of table items are as follows. &#x201c;Date/Time UTC&#x201d;: Date in YYYY-MM-DD HH: MM UTC, format. &#x201c;Phenomenon Visible&#x201d;: If the phenomenon is seen in the image. &#x201c;Moon Intersection: Earth&#x2019;s Limb Lighting Conditions&#x201d; if the Moon is on the illuminated or dark limb of the Earth. &#x201c;Imaging Type&#x201d;: Which bands were imaged. &#x201c;Moon Orientation CCD&#x201d;: Where the Moon is located on the CCD. &#x201c;Moon Intersection: Location on Earth&#x201d; Geophysical location where the Moon&#x2019;s intersection occurs.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>In an additional dataset, all 10 bands were obtained (<xref ref-type="fig" rid="F2">Figure 2</xref>), enabling a view of the limb phenomenon in different wavelengths. The phenomenon is mostly invisible in ultraviolet (UV) bands 317, 325, 340, and 388&#xa0;nm and becomes increasingly clear as the bands progress through 443, 551, 680, 688, 764, and 780&#xa0;nm. This indicates that the phenomenon is not likely due to atmospheric scattering, which is smaller in the visible and near-IR bands (<xref ref-type="fig" rid="F2">Figure 2</xref>).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Phenomenon in different wavelengths. From data taken on 30 November 2020.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g002.tif">
<alt-text content-type="machine-generated">Grayscale images showing a phenomenon at various wavelengths: 317, 325, 340, 388, 443, 551, 680, 688, 764, and 780 nanometers. Each image shows light and shadow patterns, with differences in contrast and texture between images.</alt-text>
</graphic>
</fig>
<p>Initial inspection of the Moon passing behind the Earth from other spacecraft, such as GOES and the ISS, did not yield the same phenomenon (as seen in <xref ref-type="fig" rid="F3">Figure 3</xref>). Furthermore, crew members on board the ISS who had viewed comparable Earth-Moon transits reported that no flange-like feature was visible to the naked eye, corroborating the absence of the phenomenon at low-Earth viewing geometry (<xref ref-type="bibr" rid="B20">Pettit and Hague, 2025</xref>). Both images from the International Space Station (ISS), at an altitude of &#x223c;400&#xa0;km, and Geostationary Operational Environmental Satellite (GOES) Advanced Baseline Imager (ABI), at &#x223c;35,800&#xa0;km, display expected optical compression from refraction through the atmosphere. While this image is not included in this paper, the optical compression can also be seen in images from Himawari-8 (<xref ref-type="bibr" rid="B25">Universe Space Tech, 2023</xref>). However, while this paper was in the review process, another, newly taken picture from Himawari-8 revealed a similar phenomenon to the flange. The origin of this and the situations in which it appears is discussed in <xref ref-type="sec" rid="s7">Section 7</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Views of the Moon passing behind the Earth. Left: Astronaut photography from International Space Station (ISS) (NASA via <xref ref-type="bibr" rid="B2">Astronomy.com, 2023</xref>). Right: From GOES Advanced Baseline Imager (ABI) (<xref ref-type="bibr" rid="B18">NOAA NESDIS, 2020</xref>). Both exhibit typical optical compression from atmospheric distortion.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g003.tif">
<alt-text content-type="machine-generated">Side-by-side images comparing the view of the Moon from two perspectives. The left shows the Moon above Earth&#x27;s atmosphere as seen from the International Space Station (ISS), with a blue gradient. The right displays the Moon partially above Earth&#x2019;s horizon with a clearer black sky, captured by the GOES satellite, featuring the NOAA logo.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s3">
<label>3</label>
<title>Mirage review</title>
<p>Based on the known mirage in <xref ref-type="fig" rid="F4">Figure 4</xref>, a further review of atmospheric phenomenon was conducted, and from a catalogue of possible features (<xref ref-type="bibr" rid="B16">Naylor, 2002</xref>), the following observation was made: that the appearance was very similar to a feature known as an &#x201c;omega mirage.&#x201d;</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Example of Sun &#x201c;omega mirage.&#x201d; This inferior mirage occurs when the ocean creates a surface of low-density hot air (<xref ref-type="bibr" rid="B27">Inaglory, 2007</xref>).</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g004.tif">
<alt-text content-type="machine-generated">A large, distorted image of the sun near the horizon, displaying vibrant red, orange, and yellow hues. The sun&#x27;s shape appears stretched and compressed, creating an unusual optical illusion against a dark background.</alt-text>
</graphic>
</fig>
<p>Mirages are an optical phenomenon where light is bent by different, contrasting combinations of air temperatures and humidity layers. This bending, produced by different refractive indices of the layers, causes displacement and distortion of distant objects in the sky.</p>
<p>Mirages are separated into roughly three classes: inferior, superior, and &#x201c;complex.&#x201d;</p>
<p>In inferior mirages, as seen in <xref ref-type="fig" rid="F5">Figure 5</xref>, an inversion of the light rays causes a mirror reflection of an object to appear below it. This is caused by scenarios where the air has high temperatures near the ground followed by layers of cooler temperatures as the height increases. A common example of this is the &#x201c;highway mirage&#x201d; where a hot road appears as if it has a reflective puddle, but instead it is caused by the temperature differences above the road (<xref ref-type="bibr" rid="B7">Hecht, 2002</xref>).</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Examples of superior and inferior mirages. Superior mirages, caused by layers of increasingly warmer air cause an object to appear to float. Inferior mirages, caused by layers of hot air that increasingly cools causes a mirror-like reflection of the object (<xref ref-type="bibr" rid="B11">Lorenzelli, 2014</xref>).</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g005.tif">
<alt-text content-type="machine-generated">Diagram comparing superior and inferior mirages. On the left, a superior mirage shows light bending downward as it passes through layers of warm and cold air, making the mirage appear above the real object. On the right, an inferior mirage shows light bending upward through layers of cool and hot air, making the mirage appear below the real object. Both diagrams depict observers perceiving light traveling in a straight line, illustrating the optical illusion caused by atmospheric refraction.</alt-text>
</graphic>
</fig>
<p>In superior mirages, as seen in <xref ref-type="fig" rid="F5">Figure 5</xref>, light bends in such a way to cause an object to appear to float above its original location. In this situation, layers of cool air followed by warm cause an expansion effect, which makes the objects appear higher. Superior mirages often cause objects that are visually below the line of sight to appear on the horizon when they would not in normal conditions (<xref ref-type="bibr" rid="B12">Lynch and Livingston, 1995</xref>).</p>
<p>The final mirage, a &#x201c;complex mirage,&#x201d; occurs when there are alternating layers between hot and cold air. In these mirages, known as fata morgana, the optical effects consist of inverted images, compressed and stretched layers, as well as objects &#x201c;floating&#x201d; on the horizon (<xref ref-type="bibr" rid="B26">van der Werf, 2022</xref>).</p>
</sec>
<sec id="s4">
<label>4</label>
<title>Theory development</title>
<p>Terrestrial mirages occur when different temperature layers cause changes in the air density, which leads to changes in the refractive index. Near the ground, these thermal layers form due to local conditions: for example, heat rising from the ocean or the sequence of cool air overlaid by warm air creates distinctive optical distortions.</p>
<p>It was considered that the phenomenon was perhaps based on temporary atmospheric conditions. To investigate, we used the Goddard Earth Observing System Forward Processing (GEOS-FP) data set, a product that consists of near-real-time, global weather maps. GEOS-FP merges millions of satellite and ground observations with a numerical weather model every 6&#xa0;h; this data-assimilation step produces hourly fields of temperature, humidity, winds, and other variables on a grid of roughly 25&#xa0;km resolution. (<xref ref-type="bibr" rid="B6">Global Modeling Assimilation Office, 2017</xref>). An example of GEOS data can be seen in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Vertical sample from a GEOS dataset (GEOS.fp.asm.inst3_3d_asm_Np.20230205_2100.V01.nc4) located at Goddard Space Flight Center (38.99N, 76.85W) for example purposes.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">GEOS layer &#x23;</th>
<th align="left">Altitude (km)</th>
<th align="left">Pressure (hPa)</th>
<th align="left">Temperature (K)</th>
<th align="left">Relative humidity (%)</th>
<th align="left">Air density (kg/m<sup>3</sup>)</th>
<th align="left">Refractive index <italic>n</italic> (unitless)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1</td>
<td align="left">0.17</td>
<td align="left">N/A</td>
<td align="left">N/A</td>
<td align="left">N/A</td>
<td align="left">N/A</td>
<td align="left">N/A</td>
</tr>
<tr>
<td align="left">8</td>
<td align="left">1.24</td>
<td align="left">825</td>
<td align="left">267.4004</td>
<td align="left">96.06</td>
<td align="left">1.07</td>
<td align="left">1.000241</td>
</tr>
<tr>
<td align="left">15</td>
<td align="left">3.67</td>
<td align="left">600</td>
<td align="left">253.4794</td>
<td align="left">80.99</td>
<td align="left">0.82</td>
<td align="left">1.000185</td>
</tr>
<tr>
<td align="left">22</td>
<td align="left">9.75</td>
<td align="left">250</td>
<td align="left">234.2164</td>
<td align="left">3.10</td>
<td align="left">0.37</td>
<td align="left">1.000084</td>
</tr>
<tr>
<td align="left">28</td>
<td align="left">21.97</td>
<td align="left">40</td>
<td align="left">227.2038</td>
<td align="left">0.29</td>
<td align="left">0.06</td>
<td align="left">1.000014</td>
</tr>
<tr>
<td align="left">35</td>
<td align="left">40.18</td>
<td align="left">3</td>
<td align="left">263.3156</td>
<td align="left">4.67E-04</td>
<td align="left">0.0039</td>
<td align="left">1.000001</td>
</tr>
<tr>
<td align="left">42</td>
<td align="left">65.61</td>
<td align="left">0.1</td>
<td align="left">224.6964</td>
<td align="left">1.11E-03</td>
<td align="left">0.00015</td>
<td align="left">1</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Cells outlined with a solid line were directly from the GEOS, dataset; cells outlined with a dash were derived from the formula in <xref ref-type="sec" rid="s5-2">Section 5.2</xref>. Some values are N/A because the location is lower than the mean altitude of the layer (i.e., it would be underground).</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Atmospheric parameters from the GEOS product, including temperature and humidity, were obtained and reprojected into the EPIC projection and field of view at the times of EPIC partial lunar occultation events. However, the results revealed no evidence of a relationship between temporal atmospheric conditions and whether the phenomenon appeared.</p>
<p>Further review of the data was conducted, and it was observed that all three times the phenomenon was observed when the Moon was halfway below the horizon. Although this is not statistically significant, but because no other explanation fits the observation as well, we adopt a working hypothesis: that this phenomenon was not fleeting but occurred every time the Moon passed behind the Earth and occurred while the Moon was at least halfway below the horizon. This is not as unlikely as it seems, as for the majority of the sets the imagery cadence was planned for normal, Earth-imaging operations and not intended for capturing the Moon or this event. Two out of the three datasets where the phenomenon was observed were fast-cadence imaging sequences specifically targeted for imaging the Moon and therefore more likely to catch it in the right physical position. Also, importantly, the imagery serves only as an initial clue; the hypothesis is tested chiefly by a physics-based ray tracing model, rather than imagery alone.</p>
<p>Under the assumption that this is a consistent event, we propose a new working hypothesis. The atmospheric layers, while generally characterized by a decrease in temperatures with an increase of altitude, do not follow a linear gradient. Instead, there is a switching back and forth of the temperature as altitude increases. The troposphere, stratosphere, mesosphere, and thermosphere are characterized by alternating cooling and rising temperatures. These different layers create an environment ripe for mirage-type optical distortion by providing contrasting refractive indices.</p>
</sec>
<sec sec-type="methods" id="s5">
<label>5</label>
<title>Methods</title>
<p>There are many ways to model the effect of atmospheric distortion on light. Common astronomical methods use heuristics as a function of viewing angle. Many, including the EPIC geolocation algorithm (<xref ref-type="bibr" rid="B3">Blank, 2019</xref>), use layer-based ray tracing models that treat the atmosphere as a linear, 2-dimensional pathway. However, these common methods had not revealed any type of optical phenomenon, other than compression, at distance. If the phenomenon is a mirage, resolving it requires a fully three-dimensional model rather than the simpler 2D approximations.</p>
<sec id="s5-1">
<label>5.1</label>
<title>Ray tracer model overview</title>
<p>To this end, a custom, voxel-based ray tracer utilizing vectors was developed. Ray tracers are software that estimate the path of viewing rays in a virtual scene. Voxels, a portmanteau of &#x201c;volumetric pixel,&#x201d; are essentially 3D pixels, where each voxel is populated with a property. In this case, the property is the refractive index of a physical location. Using vector math permits a complete 3D calculation of the ray paths and interactions with the Earth&#x2019;s atmosphere, determining how the rays are warped and where they intersect with the Moon. Furthermore, unlike most ray tracers that discard the intermediary products, a secondary voxel model can store all computed vectors and distortion levels, permitting inspection of the actual phenomenon in post-processing.</p>
<p>A basic renderer (<xref ref-type="fig" rid="F6">Figure 6</xref>) would consist of GEOS-derived refractive indices and an Earth and Moon model which would be used to populate a voxel model. The vector-based ray tracer would then compute on top of the model, producing vector voxel maps of the resulting view rays.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Voxel ray tracer pipeline. GEOS data is used to derive the refractive indices, which are then combined with geophysical models and used to populate the voxel model. The 3D vector-based ray tracer then calculates optical paths and stores the results the ray vector voxel map.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g006.tif">
<alt-text content-type="machine-generated">Diagram showing a workflow for generating a Ray Vector Voxel Map. Inputs include temperature, humidity, and refractive index through GEOS-derived products, and Earth/Moon models. The algorithm involves a Voxel Model and a 3D Vector Ray Tracer, leading to the final product, the Ray Vector Voxel Map.</alt-text>
</graphic>
</fig>
<p>However, this is not without challenge. A full model of the EPIC, Earth, and Moon space, consisting of these refractive voxels, has a computational and storage complexity of O(n<sup>3</sup>) where n is a function of resolution. EPIC has a resolution of &#x223c;8&#xa0;km, but obtaining clear ray paths for a potential mirage, which might contain ray inversions, requires a voxel space that is oversampled by a factor of 3&#x2013;4&#xd7; (<xref ref-type="bibr" rid="B22">Shannon, 1949</xref>). A basic calculation of a 2&#xa0;km voxel model of the Earth/Moon scene from EPIC would have a volume of &#x223c;10.8 trillion km<sup>3</sup> or &#x223c;5.4 trillion voxels. On a computer, using double precision floating point for all calculations and stored values, would require over 4.3 terabytes of memory, an impractical amount.</p>
<p>To improve resource utilization, a &#x201c;just-in-time&#x201d; renderer was developed. In this paradigm, only the portions of the scene that are immediately being computed upon are rendered. A further reduction in computation is obtained by clipping the rendered scene to a subset of the optical pathways needed for the optical effect. To ensure continuity of the calculation, the computation retains both the current and previous slice information. Every value that is needed to ensure a complete calculation is propagated.</p>
<p>The renderer is prepared by developing abstract models that contain the values needed, but not in the voxel format. Data from GEOS, including pressure, temperature, and humidity for each atmospheric layer, is ingested in equirectangular projection and the air density and subsequent refractive index are calculated. An abstract 3D model of the Earth and the Moon, in Earth-Centered, Earth-Fixed (ECEF) (<xref ref-type="bibr" rid="B17">Nishihama et al., 1997</xref>) XYZ coordinates is also separately computed.</p>
<p>During computation, instead of having the entire voxel model precomputed, one slice along the X-axis is rendered at a time. The abstract models are used to generate look-up indices into the equirectangular refraction datasets, and these values are pulled into a 2D, YZ-dimension slice. The vector renderer then calculates the resulting angular distortions from refraction and the new ray-vectors, along with distorted YZ index offsets, which are then stored in an HDF (hierarchical data format) file. If an opaque surface is encountered, or the calculation yields a reflection, the ray is &#x201c;snipped,&#x201d; and no more computations occur on it. The renderer marches through the model space, one slice at a time, until it reaches the end.</p>
<p>In sections where there is only the vacuum of space and no refractive voxels, such as the area between the Earth and the Moon, the renderer will skip drawing ray slices and compute new YZ index offsets based on the last computed ray angles and the distances covered.</p>
<p>The product of the ray tracer is a full 3D map of the viewing rays&#x2019; vectors and the amount of distortion per ray, stored in an HDF file. Making this data useful requires additional software that can pull the relative data and display in a viewable format.</p>
<p>A ray tracer, which only models the path the rays take through a scene, is not a full simulator as it does not implement a full lighting model and instrument model that would be required to duplicate the original image. A full lighting model would need to include Rayleigh scattering, which describes the wavelength-dependent way that air molecules scatter light; a radiative-transfer model, which tracks absorption, emission, scattering and the resulting gain or loss of radiance along the ray; and a physics-based shader, which determines surface brightness values based on specular and angular reflectance. In addition, to fully replicate the original image would require an optical and sensor model for the EPIC instrument. The working theory is that the distorted ray-paths are the primary contributor to the phenomenon, and that it is unnecessary to model the other physical aspects to understand the basic underlying behavior. The computational complexity of the ray tracer is very high, requiring millions of physics-based calculations, extremely reducing the likelihood of accidentally replicating the phenomenon. If the DSCOVR EPIC phenomenon and the GOES ABI compressive effect are both sourced from the refractive aspect of atmospheric optics, then a ray tracer will be able to produce both images, with the only change being the initial vectors set to imitate the instrument&#x2019;s field of view.</p>
<p>For the voxel model, while the minimum for imitating EPIC&#x2019;s resolution would require &#x223c;10&#xa0;km, this would not be adequate for interpreting the ray paths. According to Nyquist-Shannon sampling theorem (<xref ref-type="bibr" rid="B22">Shannon, 1949</xref>), the minimum sampling to resolve a signal is at least twice the resolution, which means that the voxel model needs to be at least 5&#xa0;km or better to interpret the results.</p>
</sec>
<sec id="s5-2">
<label>5.2</label>
<title>Refractive index computation</title>
<p>To calculate refraction, the GEOS-FP model, which contains 3D assimilated states at various pressure levels, is used. There are 42 layers in this set, each at a different pressure level, covering altitudes from sea level to &#x223c;65&#xa0;km. Each layer contains, in equirectangular projection, values for temperature, humidity, and pressure from which the air density and refractive index can be derived.</p>
<p>The process for calculating air density is essentially to derive the amount of dry air versus water vapor at pressure (<xref ref-type="bibr" rid="B21">Picard et al., 2008</xref>; <xref ref-type="bibr" rid="B19">OmniCalculator, 2024</xref>). The first step is to calculate saturation vapor pressure, <italic>svp</italic>, in kPa via Tetens equation (<xref ref-type="disp-formula" rid="e1">Equation 1</xref>) (<xref ref-type="bibr" rid="B15">Murray, 1967</xref>; <xref ref-type="bibr" rid="B28">Wikipedia, 2025</xref>). <italic>T</italic> is the temperature in Celsius, and <italic>RH</italic> is the relative humidity, both obtained from the GEOS model.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>v</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.61078</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>17.27</mml:mn>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>237.2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>Actual vapor pressure (<italic>avp), in kPa,</italic> is then:<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>v</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>v</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>R</mml:mi>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>Pressure (P), in units of kPa, is then used to estimate the dry air pressure, <italic>dap</italic>, (kPa) via:<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>P</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>v</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>Using 287.058&#xa0;J/(kg K) as the specific gas constant for dry air, and 461.495&#xa0;J/(kg K) as the specific gas constant for water vapor, the air density, in kg/m<sup>3</sup>, can be calculated via:<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>s</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1000</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>287.058</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1000</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>v</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>461.495</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>The refractive index (<italic>refindex</italic>) for visible light (<xref ref-type="bibr" rid="B24">Thormahlen, 1985</xref>) is then calculated (<xref ref-type="disp-formula" rid="e5">Equation 5</xref>) (<xref ref-type="bibr" rid="B13">MadSci, 1998</xref>) as follows, where 1.29&#xa0;kg/m<sup>3</sup> is the air density at room temperature and pressure and c is the speed of light in m/s. 0.00029 is based on the Edl&#xe9;n dispersion formula for standard air pressure (<xref ref-type="bibr" rid="B10">Jones, 1981</xref>) and the refractive index for air:<disp-formula id="e5">
<mml:math id="m5">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>&#x2a;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.00029</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mn>1.29</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>This calculation is done for each layer of the GEOS dataset.</p>
</sec>
<sec id="s5-3">
<label>5.3</label>
<title>Vector ray tracer</title>
<p>Refraction is modelled with Snell&#x2019;s law (<xref ref-type="disp-formula" rid="e6">Equation 6</xref>) (<xref ref-type="bibr" rid="B1">Angel, 2000</xref>), which describes that when a ray encounters a surface with a different index of refraction, the angle that the light gets transmitted through the surface is dependent on the ratio between the current and new refractive indices, as well as the angle of incidence. Essentially:<disp-formula id="e6">
<mml:math id="m6">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>Where <inline-formula id="inf1">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf2">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the incident and transmitted angle and <inline-formula id="inf3">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf4">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the refractive indices of the current and new surface.</p>
<p>To apply this property to vectors requires some adaptation (<xref ref-type="bibr" rid="B1">Angel, 2000</xref>). Two vectors are needed, one of which is the vector representing the current ray, defined as v<sub>r</sub>
<italic>.</italic> The second is for the voxel coordinate of the atmospheric surface with which the ray vector is colliding. This vector is from the Earth&#x2019;s center coordinate to the current voxel coordinate and is defined as v<sub>a</sub>.</p>
<p>The intersection angle, <inline-formula id="inf5">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, between the two vectors is computed:<disp-formula id="e7">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:msup>
<mml:mi>cos</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mo>&#xb7;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>The ratio between refractive indices is:<disp-formula id="e8">
<mml:math id="m13">
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>And the angle the light is refracted is:<disp-formula id="e9">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mi>cos</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2a;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>0.5</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>To calculate the refracted vector, v<sub>t</sub>, is then:<disp-formula id="e10">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>cos</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>If <inline-formula id="inf6">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x398;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is a complex number, the ray is encountering the critical angle and being reflected instead of transmitted. In this case, the program snips the ray as it is no longer propagated through the model. Rays are also snipped when they encounter a non-transparent surface, such as the surface of the Earth or the Moon.</p>
<p>The amount of distortion for each ray is tracked and updated through each propagation.</p>
<p>The magnitude factor <italic>f</italic> of each ray is calculated by normalizing the x component transmitted vector, since the algorithm can only advance the rays on discrete slices:<disp-formula id="e11">
<mml:math id="m17">
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>Where the new coordinates are:<disp-formula id="e12">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
<disp-formula id="equ1">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="equ2">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>z</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>z</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>z</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>Where <italic>s</italic> refers to the current slice, and <italic>s&#x2b;</italic>1 indicates results should be the new coordinates in the next slice to be computed on. Note that the x calculation can be skipped since it will always advance by 1 due to the discrete nature of the voxel model; it is shown here for completeness.</p>
<p>Before execution of the ray tracer, it is necessary to set the state of the rays based on the optical configuration of the instrument. The initial ray vector slice that starts the process is calculated based on the field of view (FOV) of the instrument. Where the FOV angle is <inline-formula id="inf7">
<mml:math id="m21">
<mml:mrow>
<mml:mi>&#x3c6;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, F is the maximum vector.<disp-formula id="e13">
<mml:math id="m22">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3c6;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>The rate of change across the slice is calculated, where vol is the voxel model dimensions in the y and z-axis:<disp-formula id="e14">
<mml:math id="m23">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
<disp-formula id="equ3">
<mml:math id="m24">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>z</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
<mml:mi>z</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>Then each ray vector for the initialization slice is calculated, where i is the column and j is the row, and <italic>init</italic> is the initial ray slice:<disp-formula id="e15">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>F</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
<disp-formula id="equ4">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>z</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>F</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:mi>z</mml:mi>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>The result is an array where the center vector has an offset angle close to zero, and from there the offset angles symmetrically spread as they approach the edge of the volumetric space.</p>
</sec>
<sec id="s5-4">
<label>5.4</label>
<title>Full ray tracing process</title>
<p>The ray tracer takes the above algorithms and pulls it all together in the following steps:<list list-type="order">
<list-item>
<p>The initial parameters for the voxel model are set up according to the scene, and dimensions are determined by the locations of the Earth and the Moon and the voxel resolution.</p>
</list-item>
<list-item>
<p>The GEOS-FP dataset is ingested and the refractive indices for the atmospheric layers are calculated (<xref ref-type="disp-formula" rid="e1">Equations 1</xref>&#x2013;<xref ref-type="disp-formula" rid="e5">5</xref>). If necessary, neighboring layers can be averaged together to meet the resolution of the voxel model.</p>
</list-item>
<list-item>
<p>The initial ray vector states are calculated (<xref ref-type="disp-formula" rid="e13">Equations 13</xref>&#x2013;<xref ref-type="disp-formula" rid="e15">15</xref>), based on the FOV angle. This becomes the current ray vector.</p>
</list-item>
<list-item>
<p>The ray tracer then runs iteratively through the voxel model, advancing forward slice by slice until it reaches the end. This consists of the following steps:</p>
<list list-type="alpha-lower">
<list-item>
<p>The current slice is rendered by pulling the relevant refractive indices from the layer model into their proper location in the slice.</p>
</list-item>
<list-item>
<p>The Earth&#x2019;s center coordinate to the current slice coordinate vector is built and the algorithm for calculating the transmitted angle is executed (<xref ref-type="disp-formula" rid="e6">Equations 6</xref>&#x2013;<xref ref-type="disp-formula" rid="e10">10</xref>). Any rays that are reflected or encounter a body are snipped. If a body is encountered, the type is stored (i.e. Earth vs. Moon).</p>
</list-item>
<list-item>
<p>The amount of distortion for each ray is calculated (<xref ref-type="disp-formula" rid="e11">Equations 11</xref> and <xref ref-type="disp-formula" rid="e12">12</xref>).</p>
</list-item>
<list-item>
<p>Data concerning the ray vectors, refraction, angle of incidence, and transmitted angles are stored in a file for future analysis.</p>
</list-item>
<list-item>
<p>Ray information is propagated to the next slice calculation to ensure continuity between slices.</p>
</list-item>
</list>
</list-item>
<list-item>
<p>The final simulated image is calculated based on the states of the ray vectors and the calculated level of distortion.</p>
</list-item>
</list>
</p>
</sec>
</sec>
<sec sec-type="results" id="s6">
<label>6</label>
<title>Results</title>
<p>The ray tracer was run at 2&#xa0;km resolution, resulting in <xref ref-type="fig" rid="F7">Figures 7</xref>, <xref ref-type="fig" rid="F8">8</xref>. As seen in the raytraced image compared to the actual image from EPIC, the same phenomenon is seen. A vertical sampling of the rays reveals a maximum 0.45&#xb0; distortion that decreases with distance away from the Earth. This is not unexpected, as the refractive indices decrease in relation to the air density, but it is interesting to note that the descent has several &#x201c;bumps&#x201d; where it increases instead of decreases.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Left: EPIC image of the Moon limb phenomenon. Right: Results from ray tracer. Because this model uses discrete rays, small black dots occur where there is no ray present due to diversion from refraction.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g007.tif">
<alt-text content-type="machine-generated">Two-panel image comparing features. Left: a grayscale image with a dark and light gradient transitioning to a textured surface. Right: a color-coded version with black, gray, and blue sections corresponding to the gradient and surface for clearer differentiation.</alt-text>
</graphic>
</fig>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Above: Diagrams that show where the data was obtained from the simulation. Left: Per ray magnitude of distortion along vertical axis. Right: Per ray magnitude of distortion along horizontal axis.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g008.tif">
<alt-text content-type="machine-generated">Two-section graphic showing ray distortion data. The top displays two rectangular plots with color gradients, possibly illustrating distortion patterns. Bottom left is a bar chart titled &#x22;DSCOVR Magnitude of Ray Distortion along Vertical (Y) Axis,&#x22; showing a decrease in distortion as distance increases from Earth&#x27;s Y-axis center. Bottom right features a bar chart titled &#x22;EPIC Magnitude of Ray Distortion along Horizontal (Z) Axis,&#x22; indicating variation in distortion magnitude with distance from Earth&#x27;s Z-axis center.</alt-text>
</graphic>
</fig>
<p>An inspection of the rays as they traverse through the atmosphere (<xref ref-type="fig" rid="F9">Figure 9</xref>) shows that the rays are primarily traveling through the stratosphere (&#x223c;20&#x2013;50&#xa0;km altitude) and lower mesosphere (&#x223c;50&#x2013;85&#xa0;km).</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Illustration of the ray paths. This shows how the rays travel through the atmosphere, and how they are affected by the atmosphere. The rays are bent very slightly in this image; the effect seen in phenomenon is due to a combination of this slight angular distortion and the distance between the Earth and the Moon. This view helps account for what happens to UV light, such as the EPIC 317, 325, and 340&#xa0;nm bands, as they pass through the stratosphere and would be filtered out by the long path alone the ozone layers. The rays that hit the troposphere approach the critical angle and are reflected.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g009.tif">
<alt-text content-type="machine-generated">Graph showing EPIC ray path along the Y-axis of the atmosphere, with atmospheric height on the Y-axis and distance from the spacecraft on the X-axis. Layers indicated are the troposphere, stratosphere, and mesosphere. Blue shaded areas represent different atmospheric layers.</alt-text>
</graphic>
</fig>
<p>This would account for the apparent lack of the phenomenon in UV light, as the light passes almost horizontally through the lower stratosphere (15&#x2013;30&#xa0;km) and it would be filtered out by the ozone layer. Although computing the atmospheric loss for scattering and absorption is outside of intent of this paper, for reference the loss at 41.8&#xb0; above the horizon is approximately 99.8% at 300&#xa0;nm and 47.8% at 350&#xa0;nm (<xref ref-type="bibr" rid="B23">Spectra-AM1.5, 2025</xref>). The rays that hit the troposphere approach the critical angle and are reflected.</p>
<p>To ensure that the model&#x2019;s result is not affected by the snipping of reflected rays, a review of a number of calculations that results in reflected rays was conducted. Of the 245,592,557 refraction calculations that were conducted, only 776 resulted in a reflected ray, or 0.0003%. This indicates that reflection does not play a statistically significant role in the construction of the final image.</p>
<p>The model was run on a very small segment at 1&#xa0;km resolution; the smaller voxel size did not alter the image morphology.</p>
<p>Viewing the full ray paths (<xref ref-type="fig" rid="F10">Figure 10</xref>) shows the degree of distortion and where the pixels of the phenomenon originate. The small slice of atmosphere and the distortion caused over the 384,400&#xa0;km distance between the Earth and the Moon generates over a 3,000&#xa0;km bending of the rays.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Top: Slice from ray tracer results that shows the full path of viewing rays and the level of divergence that is caused by atmospheric refraction and the distance between the Earth and the Moon. The location of the Moon is shown in green on the right. Colors of the rays have no meaning other than to distinguish them from each other. Bottom: Zoomed-in view of rays above. This magnified view reveals several ray inversions (circled in red), indicating a complex mirage. Color of the rays have no meaning other than to distinguish them from each other. Slice is taken from same location as seen in <xref ref-type="fig" rid="F8">Figure 8</xref>, top left.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g010.tif">
<alt-text content-type="machine-generated">Graphical representation of the ray paths from EPIC (Earth Polychromatic Imaging Camera) along the Y-Axis. The top graph shows the paths from Earth&#x27;s atmosphere to the Moon&#x27;s surface, while the bottom zooms into a segment closer to the atmosphere. Colorful lines represent different ray paths with Earth&#x2019;s outline depicted on the left. The distance axes are labeled in kilometers, and red circles emphasize specific areas in the lower zoomed graph.</alt-text>
</graphic>
</fig>
<p>Zooming in on the rays provides more detail in <xref ref-type="fig" rid="F11">Figure 11</xref>. Here, several inversions of the light rays can be seen, showing that this phenomenon is not just optical distortion, but a complex mirage. It is interesting to note that while the appearance of the phenomenon is similar to an &#x201c;omega mirage,&#x201d; the effect is obtained through different means.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Zoomed-in view of Z-axis slice. Many ray inversions (circled in red) are seen, indicating aspects of a complex mirage. Color of the rays have no meaning other than to distinguish them from each other. Slice is taken from same location as seen in <xref ref-type="fig" rid="F8">Figure 8</xref>, top right.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g011.tif">
<alt-text content-type="machine-generated">Diagram showing ray paths through Earth&#x27;s atmosphere with multiple colored lines crossing from the Earth on the left to the direction of the Moon on the right. The Earth is represented by a blue shaded area labeled &#x22;Not to scale.&#x22; Red circles highlight intersections of the lines, and annotations describe the &#x22;Length of atmosphere.&#x22; The axes measure distances from Earth&#x27;s center in kilometers.</alt-text>
</graphic>
</fig>
<p>Inspection of the rays along the Z-axis yields a different picture. Here, the magnitude of distortion is less, but with much greater variance along the slice.</p>
<p>Viewing the rays shows a much more complex distortion in the horizontal vs. vertical axis. A zoomed-in view shows much more mirage-like ray inversions.</p>
<p>While the analysis shows that the phenomenon has aspects of a complex mirage, it does not demonstrate the differences that occur between the views seen by EPIC versus geostationary imagery.</p>
</sec>
<sec id="s7">
<label>7</label>
<title>Comparative analysis</title>
<p>In order to determine how the GOES scene differs from EPIC, it is necessary to run the ray tracer from the GOES point of view. The GOES ABI has a 17.76&#xb0; field of view, and by initializing the first vectors with corresponding offsets, it is possible to create a GOES model of the scene.</p>
<p>Running the ray tracer for GOES produced the matching optical compressive effect, as seen in <xref ref-type="fig" rid="F12">Figure 12</xref>, <italic>top left</italic>. When viewing the magnitudes of distortion (<xref ref-type="fig" rid="F13">Figure 13</xref>), the degree of distortion is in correlation to the density of the atmosphere according to altitude. Unlike that seen in the EPIC data, there are no &#x201c;bumps&#x201d; or aberrations that cause the ray inversions. This is true of both the vertical and horizontal slices. Unlike the EPIC phenomenon, this is only a straightforward, atmospheric lensing effect, and has no evidence of being a complex mirage.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Top left: GOES ABI image showing atmospheric distortion. In this image, the optical properties of the atmosphere &#x201c;compress&#x201d; the Moon. Top right: Simulated scene from ray tracer using GOES ABI field of view, shows similar effect. Bottom left: Himawari-8 image showing the limb phenomenon. In both simulated images, because the voxel model uses discrete rays, small black dots occur where there is no ray present due to diversion from refraction.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g012.tif">
<alt-text content-type="machine-generated">Two images on the left show the moon rising over a blurry horizon with a gradient from gray to blue. On the right are corresponding silhouette representations against a black and blue background.</alt-text>
</graphic>
</fig>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>GOES ABI/Himawari-8 simulated magnitudes of refractive distortion. Left: Angular distortion of rays on horizontal axis. Right: Angular distortion of rays on vertical axis. Distances are negative since they are left of the center coordinate. In these charts, the rays bend in a more regular decline than EPIC, showing no evidence of ray inversions.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g013.tif">
<alt-text content-type="machine-generated">Two bar graphs illustrating the GOES magnitude of ray distortion in degrees. The left graph shows distortion along the vertical (Y) axis, peaking at 0.15 degrees at around 6255 km. The right graph shows distortion along the horizontal (Z) axis, gradually increasing to a peak of about 0.12 degrees near -3980 km. Both graphs use kilometers as the unit for distance from Earth&#x27;s axis center reference.</alt-text>
</graphic>
</fig>
<p>However, when looking at the magnitude of the rays for geostationary, as seen in <xref ref-type="fig" rid="F13">Figure 13</xref>, there is a substantial difference between what is seen for GOES, Himawari-8, and EPIC. The magnitude of distortion is far less, plus there is no evidence of the magnitude changes which would indicate the ray inversions as seen in <xref ref-type="fig" rid="F11">Figure 11</xref>. (This <xref ref-type="fig" rid="F11">Figure 11</xref> plot is not duplicated for GOES/Himawari-8 because the angular distortion due to these instruments&#x2019; wide field of view renders the relatively small effect of atmospheric disturbance invisible. The addition of this plot would exceed the number of figures permitted by the publication.)</p>
<p>It was initially believed that this limb effect was unique to EPIC. However, in an image taken by Himawari-8 on 11 August 2025, after an initial draft of this paper had been submitted for publication, showed otherwise, as seen in <xref ref-type="fig" rid="F12">Figure 12</xref>, bottom left. We were able to replicate this, without modification to the ray tracer, after setting the Earth and Moon configuration to represent the Himawari-8 scene. The reason it was originally missed was that the time interval at which the geostationary version of the phenomenon occurs is very short. EPIC&#x2019;s phenomenon occurs for over a course of &#x223c;30&#xa0;min. Empirical analysis of Himawari-8&#x2019;s phenomenon using the ray tracer indicates that it is likely visible for less than 2&#xa0;min. The shorter time span is partially due to the lesser amount of distortion seen at geostationary versus L<sub>1</sub>, approximately 2.5 times less. Additionally, the Moon&#x2019;s relative motion is faster from a geostationary platform. Locked to the Earth&#x2019;s rotation, the spacecraft is moving at a rate of 15.04&#xb0; per hour compared to L<sub>1</sub>&#x2019;s 0.041&#xb0; per hour.</p>
<p>What is happening in the EPIC image that is different from the geostationary view? How can two pictures of the same apparent scene cause a different optical phenomenon? The answer lies in the fact that although the images appear similar, that similarity is itself an illusion.</p>
<p>DSCOVR orbits the Earth-Sun Lagrange-1 point at over 1.5 million kilometers away from the Earth. Taking the pictures requires a telescope with a 0.62&#xb0; field of view. GOES, on the other hand, is in geostationary orbit, a relatively cozy 35,000&#xa0;km away from the Earth. Capturing the entire Earth in a single frame requires the instrument to image at a wide angle of 17.76&#xb0;. (Himawari-8 is very similar at 17.8&#xb0; field of view) (<xref ref-type="bibr" rid="B9">Japan Meteorological Agency, 2015</xref>). This influences the scene; when the Moon is visible to EPIC, it is obscured by the Earth for GOES.</p>
<p>When considering the angles and slice of the atmosphere through which the two instruments view the scene, EPIC&#x2019;s view is across a much wider swath of the atmosphere, while GOES&#x2019;s is a narrow section near the center (<xref ref-type="fig" rid="F14">Figure 14</xref> bottom left and right). This changes the entry angles of the rays; the more extreme angle for EPIC causes a corresponding increase in magnitude of refraction, which causes a more extreme reaction when the ray collides with the different refractive indices. This is why the ray inversions occur with EPIC, but not in GOES or other platforms closer to Earth, such as ISS or other lower orbit satellites.</p>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Top Image: Difference in perspective between DSCOVR and GOES. The EPIC instrument has a narrow field of view, so the rays are more parallel, and the GOES ABI has a wide field of view. Although the images appear similar, they are obtained through different means. Bottom left: Horizontal view angles without atmospheric distortion. EPIC has a 0.62&#xb0; field of view at 1.5 million km vs. GOES ABI&#x2019;s 17.76&#xb0; at 35,000&#xa0;km. Bottom right: Horizontal view angles with atmospheric refraction distortion. EPIC&#x2019;s rays enter closer to the edge of the atmospheric &#x201c;lens&#x201d;, which causes more bending of the rays than what is seen with GOES ABI.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g014.tif">
<alt-text content-type="machine-generated">Comparison of DSCOVR and GOES perspective angles. The upper diagram shows DSCOVR with a 0.62-degree field of view and GOES with a 17.76-degree field of view, both observing Earth and the Moon. The lower section illustrates view angles with and without distortion. On the left, naive angles for DSCOVR and GOES are shown. On the right, refracted view angles with intersection angles for both DSCOVR and GOES are depicted. Earth and the Moon are labeled in each diagram.</alt-text>
</graphic>
</fig>
<p>When looking at the GOES results in <xref ref-type="fig" rid="F13">Figure 13</xref> vs. EPIC&#x2019;s in <xref ref-type="fig" rid="F8">Figure 8</xref> it can be seen that the magnitude of distortion has much more consistent delineation. There is no evidence of ray inversions, only optical distortion from the atmosphere. Ray slice figures, such as seen for EPIC in <xref ref-type="fig" rid="F10">Figure 10</xref>, do not appear for GOES because the ray paths are so heavily dominated by the wide-angle view of the instrument that it hides the effect of refraction.</p>
<p>In the GOES and Himawari-8 versions of distortion, the vertical distortion is &#x223c;0.2&#xb0;, but the horizontal distortion is comparatively weak at &#x223c;0.02&#xb0;. The vertical distortion dominates the resulting images. In EPIC, the vertical distortion reaches a maximum of &#x223c;0.45&#xb0;, and the horizontal distortion, while at &#x223c;0.11&#xb0; still has a substantial impact to the image. Although the general shapes of the distortion look similar for both spacecraft, on close inspection of the Himawari-8 image, the craters on the Moon can be clearly made out, although warped. Unfortunately, the resolution for the EPIC data is far less, however at a resolution similar to Himawari-8, the images would likely look very different due to the additional horizontal distortion seen from EPIC. The mirage-like components caused by the ray inversions, seen only in the EPIC results, would create additional distortion in the instrument&#x2019;s images, possibly creating a more fata morgana-like appearance. The end result is two phenomena that are structurally very different from each other.</p>
</sec>
<sec id="s8">
<label>8</label>
<title>Model validation</title>
<p>The above results were presented to the DSCOVR/EPIC operations group. They agreed to conduct a &#x201c;special imaging&#x201d; activity with the instrument in order to capture the phenomenon. The operations group advised imaging a single wavelength band, as that would permit the shortest interval between exposures. For the activity, it was decided to image the Moon as it passed behind the Earth at 3-min intervals in the 780&#xa0;nm wavelength. The short interval between exposures would permit capturing the phenomenon in different stages; 780&#xa0;nm, as seen in <xref ref-type="fig" rid="F2">Figure 2</xref>, would provide the clearest view. In this configuration, along with the limits of onboard storage and downlink, would permit 40 images to be collected over a period of 2&#xa0;h.</p>
<p>Based on the orbital dynamics, a favorable configuration of the Earth and Moon, permitting a transit, was found for 11 May 2025, from 21:30 to 23:30 UTC.</p>
<p>The ray tracer, as seen in <xref ref-type="fig" rid="F15">Figure 15</xref>, was able to correctly predict the shape the flange-like shape as it begins broad and then shrinks as the Moon rises, until it eventually disappears once it is more than halfway above the limb. This demonstrates that the working hypothesis, that the phenomenon disappears when the Moon is halfway above the limb, was correct. The changes in shape are due to the movement of the Moon from the path of the highly bent rays that pass through the lower atmospheric layers, into lesser bent rays of the upper atmospheric layers.</p>
<fig id="F15" position="float">
<label>FIGURE 15</label>
<caption>
<p>Lunar transit from 11 May 2025. The ray tracer correctly predicts the behavior of the mirage, a wide swath that shrinks as the Moon rises, until it eventually disappears once the Moon is greater than halfway above the limb. The rays are &#x201c;longer&#x201d; in the simulation, as it does not include an atmospheric absorption or an instrument sensitivity model, which would attenuate dimmer rays.</p>
</caption>
<graphic xlink:href="frsen-06-1640320-g015.tif">
<alt-text content-type="machine-generated"> A four-panel graphic, each panel split into a top and bottom section. In the top sections, a gray circular shape partially appears from the right against a black background, bordered on the right by a bright blue curved strip. Across the panels from left to right, more of the gray circle becomes visible and the blue strip becomes narrower. In the bottom sections, grayscale images show the Moon rising from behind a darker curved foreground limb. The Moon's edge appears distorted and stretched in the first panel and becomes more normally rounded as it rises in the later panels</alt-text>
</graphic>
</fig>
<p>
<xref ref-type="table" rid="T3">Table 3</xref> contains the differences in length between the phenomenon as predicted in the simulation and the measured. While the level of difference may seem to be large, it is not unexpected as predicted rays become more spread out at the edges and the phenomenon becomes very dim at the edges in the original level 1a data, averaging about 400 counts/second. For contrast, the middle of the phenomenon is &#x223c;3,000 counts/second, the mean count value of the Earth is &#x223c;14,000 counts/second, and the mean count for the dark space around Earth is &#x223c;80 counts/second. The quantum efficiency of EPIC at 780&#xa0;nm is 58%, which would attenuate the signal more (<xref ref-type="bibr" rid="B5">Epic.gsfc.nasa.gov, 2024</xref>). Atmospheric loss from absorption and scattering would contribute as well, although not as significantly. Although computing the atmospheric loss from scattering and absorption is outside of the intent of this paper, for reference it is &#x223c;4.65% at 780&#xa0;nm at 41.8&#xb0; elevation (<xref ref-type="bibr" rid="B23">Spectra-AM1.5, 2025</xref>).</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Comparison of lunar limb phenomenon length between simulated and EPIC images.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Date/Time (UTC)</th>
<th align="left">Predicted limb length (pixels)</th>
<th align="left">Actual limb length (pixels)</th>
<th align="left">Limb length difference (pixels)</th>
<th align="left">Percent length match</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">2025-05-11 21:45</td>
<td align="left">243</td>
<td align="left">210</td>
<td align="left">33</td>
<td align="left">86%</td>
</tr>
<tr>
<td align="left">2025-05-11 21:57</td>
<td align="left">103</td>
<td align="left">50</td>
<td align="left">53</td>
<td align="left">48%</td>
</tr>
<tr>
<td align="left">2025-05-11 22:09</td>
<td align="left">52</td>
<td align="left">14</td>
<td align="left">38</td>
<td align="left">26%</td>
</tr>
<tr>
<td align="left">2025-05-11 22:21</td>
<td align="left">0</td>
<td align="left">0</td>
<td align="left">0</td>
<td align="left">100%</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Measurement is from the edge of the Moon to the visible end of the limb. This level of attenuation is not unexpected due to quantum efficiency of EPIC, at 780&#xa0;nm wavelength.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>The full set of images are available on the EPIC Gallery (<ext-link ext-link-type="uri" xlink:href="https://epic.gsfc.nasa.gov/galleries">https://epic.gsfc.nasa.gov/galleries</ext-link>). Many thanks to NOAA and NASA for conducting operations and permitting use of the instrument to investigate atmospheric optical properties.</p>
</sec>
<sec sec-type="conclusion" id="s9">
<label>9</label>
<title>Conclusion</title>
<p>The findings demonstrate that the phenomenon seen in the EPIC image when the Moon is at the horizon is a combination of unique atmospheric distortion effects including refraction and ray inversions, making it similar to a complex mirage caused by regular differences in temperature in the atmospheric layers. Although a visually similar effect is seen at geostationary, caused solely by atmospheric refraction, it lacks the ray inversions and distortions common with mirages. EPIC&#x2019;s phenomenon, &#x201c;Gaia&#x2019;s Crown&#x201d;, is an optical effect between the Earth and the Moon when the two visually intersect and is seen only in deep space due to the unique viewing geometry causing greater distortion from atmospheric refraction. The demonstrated sensitivity to refraction and high-altitude temperature inversions hints that L<sub>1</sub> observations could provide a coarse probe of atmospheric structure. The technique could be applied to other planets with significant atmosphere to help understand their composition.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s10">
<title>Data availability statement</title>
<p>Public input datasets used in this study are available at: DSCOVR/EPIC browse (color composites): <ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="https://epic.gsfc.nasa.gov">https://epic.gsfc.nasa.gov</ext-link>; DSCOVR/EPIC science data: <ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="https://search.earthdata.nasa.gov/">https://search.earthdata.nasa.gov/</ext-link>; DSCOVR/EPIC lunar-transit browse images: <ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="https://epic.gsfc.nasa.gov/galleries">https://epic.gsfc.nasa.gov/galleries</ext-link>; GEOS (e.g., GEOS-FP) meteorological fields: <ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="https://gmao.gsfc.nasa.gov">https://gmao.gsfc.nasa.gov</ext-link>; Derived materials from this study (e.g., the ray tracer, configuration files, and rendered outputs) are subject to U.S. Government release procedures and require prior approval by the sponsoring agency; release timelines are outside the authors&#x2019; control. The algorithmic approach and parameter settings necessary to understand and reproduce the identification are described in Methods. Requests to access controlled materials may be directed to the corresponding author (who will forward them for agency review).</p>
</sec>
<sec sec-type="author-contributions" id="s11">
<title>Author contributions</title>
<p>KB: Conceptualization, Investigation, Methodology, Software, Visualization, Writing &#x2013; original draft, Writing &#x2013; review and editing. JH: Conceptualization, Validation, Writing &#x2013; review and editing. SD: Software, Validation, Writing &#x2013; review and editing. AM: Writing &#x2013; review and editing. AT: Data curation, Writing &#x2013; review and editing.</p>
</sec>
<sec sec-type="COI-statement" id="s13">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="ai-statement" id="s14">
<title>Generative AI statement</title>
<p>The authors declare that no Generative AI was used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p>
</sec>
<sec sec-type="disclaimer" id="s15">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<fn-group>
<fn fn-type="custom" custom-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/101895/overview">Soo Chin Liew</ext-link>, National University of Singapore, Singapore</p>
</fn>
<fn fn-type="custom" custom-type="reviewed-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1066065/overview">Santo V. Salinas</ext-link>, National University of Singapore, Singapore</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1183221/overview">V&#xed;ctor Molina Garc&#xed;a</ext-link>, German Aerospace Center (DLR), Germany</p>
</fn>
</fn-group>
<ref-list>
<title>References</title>
<ref id="B1">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Angel</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2000</year>). <source>Interactive computer graphics: a top-down approach with OpenGL</source>. <edition>2nd ed</edition>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Addison Wesley Longman, Inc</publisher-name>.</mixed-citation>
</ref>
<ref id="B2">
<mixed-citation publication-type="book">
<collab>Astronomy.com</collab>. (<year>2023</year>). <source>A distorted view of a full moon intersecting earth&#x27;s horizon was photographed from the international space station</source>. <publisher-name>NASA</publisher-name>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.astronomy.com/observing/everything-you-need-to-know-about-the-Moon/">https://www.astronomy.com/observing/everything-you-need-to-know-about-the-Moon/</ext-link> (Accessed September 7, 2023)</comment>.</mixed-citation>
</ref>
<ref id="B3">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Blank</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2019</year>). <source>EPIC geolocation and color imagery algorithm, revision 6</source>. <publisher-name>NASA Langley Research Center, Atmospheric Science Data Center</publisher-name>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://asdc.larc.nasa.gov/documents/dscovr/DSCOVR_EPIC_Geolocation_V03.pdf">https://asdc.larc.nasa.gov/documents/dscovr/DSCOVR_EPIC_Geolocation_V03.pdf</ext-link> (Accessed August 1, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B4">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cede</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>McCauley</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Herman</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Blank</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Kowalewski</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Raw EPIC data calibration</article-title>. <source>Front. Remote Sens.</source> <volume>2</volume>, <fpage>702275</fpage>. <pub-id pub-id-type="doi">10.3389/frsen.2021.702275</pub-id>
</mixed-citation>
</ref>
<ref id="B5">
<mixed-citation publication-type="web">
<collab>Epic.gsfc.nasa.gov</collab> (<year>2024</year>). <article-title>What is EPIC</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://epic.gsfc.nasa.gov/about/epic">https://epic.gsfc.nasa.gov/about/epic</ext-link> (Accessed August 1, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B6">
<mixed-citation publication-type="web">
<collab>Global Modeling Assimilation Office</collab> (<year>2017</year>). <article-title>File specification for GEOS-5 FP (forward processing)</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://gmao.gsfc.nasa.gov/pubs/docs/Lucchesi1202.pdf">https://gmao.gsfc.nasa.gov/pubs/docs/Lucchesi1202.pdf</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B7">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Hecht</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2002</year>). <source>Optics</source>. <edition>Fourth Edition</edition>. <publisher-loc>San Francisco, CA</publisher-loc>: <publisher-name>Pearson Education Inc</publisher-name>.</mixed-citation>
</ref>
<ref id="B8">
<mixed-citation publication-type="web">
<person-group person-group-type="author">
<name>
<surname>Inaglory</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Schematic diagram explaining fata Morgana</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://commons.wikimedia.org/wiki/File:Fada_morgana_graphnn.JPG">https://commons.wikimedia.org/wiki/File:Fada_morgana_graphnn.JPG</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B27">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Inaglory</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2007</year>). <source>Sunset inferior mirage</source>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://commons.wikimedia.org/wiki/File:Sunset_inferio_mirage.jpg">https://commons.wikimedia.org/wiki/File:Sunset_inferio_mirage.jpg</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B9">
<mixed-citation publication-type="web">
<collab>Japan Meteorological Agency</collab> (<year>2015</year>). <article-title>Himawari-8/9 Himawari standard data user&#x2019;s guide, version 1.2</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.data.jma.go.jp/mscweb/en/himawari89/space_segment/hsd_sample/HS_D_users_guide_en_v12.pdf">https://www.data.jma.go.jp/mscweb/en/himawari89/space_segment/hsd_sample/HS_D_users_guide_en_v12.pdf</ext-link> (Accessed September 18, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B10">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jones</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>1981</year>). <article-title>The refractivity of air</article-title>. <source>J. Res. Natl. Bureau Stand.</source> <volume>86</volume> (<issue>1</issue>), <fpage>27</fpage>. <pub-id pub-id-type="doi">10.6028/jres.086.002</pub-id>
<pub-id pub-id-type="pmid">34566037</pub-id>
</mixed-citation>
</ref>
<ref id="B11">
<mixed-citation publication-type="web">
<person-group person-group-type="author">
<name>
<surname>Lorenzelli</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Basic diagram of the optical illusion</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://commons.wikimedia.org/w/index.php?curid=37081455">https://commons.wikimedia.org/w/index.php?curid&#x3d;37081455</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B12">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Lynch</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Livingston</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>1995</year>). <source>Color and light in nature</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>.</mixed-citation>
</ref>
<ref id="B13">
<mixed-citation publication-type="web">
<collab>MadSci</collab> (<year>1998</year>). <article-title>What is the speed of light as a function of air density?</article-title> <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.madsci.org/posts/archives/feb98/888690999.Ph.r.html">https://www.madsci.org/posts/archives/feb98/888690999.Ph.r.html</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B14">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Marshak</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Herman</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Szabo</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Blank</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Cede</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Carn</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Earth observations from DSCOVR/EPIC instrument</article-title>. <source>Bull. Am. Meteor. Soc.</source> <volume>99</volume>, <fpage>1829</fpage>&#x2013;<lpage>1850</lpage>. <pub-id pub-id-type="doi">10.1175/BAMS-D-17-0223.1</pub-id>
<pub-id pub-id-type="pmid">30393385</pub-id>
</mixed-citation>
</ref>
<ref id="B15">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Murray</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>1967</year>). <article-title>On the computation of saturation pressure</article-title>. <source>Am. Meteorological Soc</source>.</mixed-citation>
</ref>
<ref id="B16">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Naylor</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2002</year>). <source>Out of the blue: a 24-hour skywatcher&#x2019;s guide</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>.</mixed-citation>
</ref>
<ref id="B17">
<mixed-citation publication-type="web">
<person-group person-group-type="author">
<name>
<surname>Nishihama</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wolfe</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Solomon</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>MODIS level 1A Earth location: algorithm theoretical basis document version 3.0</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://modis.gsfc.nasa.gov/data/atbd/atbd_mod28_v3.pdf">https://modis.gsfc.nasa.gov/data/atbd/atbd_mod28_v3.pdf</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B18">
<mixed-citation publication-type="web">
<collab>NOAA NESDIS</collab> (<year>2020</year>). <article-title>An ode to the moon: how NOAA satellites capture earth&#x2019;s satellite</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.nesdis.noaa.gov/news/ode-the-Moon-how-noaa-satellites-capture-earths-satellite">https://www.nesdis.noaa.gov/news/ode-the-Moon-how-noaa-satellites-capture-earths-satellite</ext-link> (Accessed August 1, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B19">
<mixed-citation publication-type="web">
<collab>OmniCalculator</collab> (<year>2024</year>). <article-title>Air density calculator</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.omnicalculator.com/physics/air-density">https://www.omnicalculator.com/physics/air-density</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B20">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Pettit</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Hague</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2025</year>). <source>Expedition 72 postflight presentation</source>. <publisher-name>NASA Goddard Space Flight Center Video</publisher-name>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://video.ibm.com/recorded/134509377">https://video.ibm.com/recorded/134509377</ext-link> (Accessed September 18, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B21">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Picard</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Davis</surname>
<given-names>R. S.</given-names>
</name>
<name>
<surname>Gl&#xe4;ser</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Fujii</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Revised formula for the density of moist air (CIPM-2007)</article-title>. <source>Metrologia</source> <volume>45</volume>, <fpage>149</fpage>&#x2013;<lpage>155</lpage>. <pub-id pub-id-type="doi">10.1088/0026-1394/45/2/004</pub-id>
</mixed-citation>
</ref>
<ref id="B22">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Shannon</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>1949</year>). &#x201c;<article-title>Communication in the presence of noise</article-title>,&#x201d; in <source>Proceedings in the institute of radio engineers</source>. <pub-id pub-id-type="doi">10.1109/JRPROC.1949.232969</pub-id>
</mixed-citation>
</ref>
<ref id="B23">
<mixed-citation publication-type="book">
<collab>Spectra-am1.5</collab> (<year>2025</year>). <source>Reference air mass 1.5 spectra</source>. <publisher-name>National Renewable Energy Laboratory</publisher-name>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.nrel.gov/grid/solar-resource/spectra-am1.5">https://www.nrel.gov/grid/solar-resource/spectra-am1.5</ext-link>.</comment>
</mixed-citation>
</ref>
<ref id="B24">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thormahlen</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Straub</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Grigull</surname>
<given-names>U.</given-names>
</name>
</person-group> (<year>1985</year>). <article-title>Refractive index of water and its dependence on wavelength, temperature, and density</article-title>. <source>J. Phys. Chem. Reference Data</source> <volume>14</volume>, <fpage>933</fpage>&#x2013;<lpage>945</lpage>. <pub-id pub-id-type="doi">10.1063/1.555743</pub-id>
</mixed-citation>
</ref>
<ref id="B25">
<mixed-citation publication-type="web">
<collab>Universe Space Tech</collab> (<year>2023</year>). <article-title>The moon &#x201c;photobombs&#x201d; the Earth in the image of a Japanese satellite</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://universemagazine.com/en/the-Moon-photobombs-the-earth-in-the-image-of-a-japanese-satellite/">https://universemagazine.com/en/the-Moon-photobombs-the-earth-in-the-image-of-a-japanese-satellite/</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
<ref id="B26">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>van der Werf</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Novaya zemlya effect and fata morgana. Raytracing in a spherically non-symmetric atmosphere</article-title>. <source>Comptes Rendus Phys.</source> <volume>23</volume>, <fpage>365</fpage>&#x2013;<lpage>389</lpage>. <pub-id pub-id-type="doi">10.5802/crphys.102</pub-id>
</mixed-citation>
</ref>
<ref id="B28">
<mixed-citation publication-type="web">
<collab>Wikipedia</collab> (<year>2025</year>). <article-title>Tetens equation</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://en.wikipedia.org/wiki/Tetens_equation">https://en.wikipedia.org/wiki/Tetens_equation</ext-link> (Accessed April 8, 2025)</comment>.</mixed-citation>
</ref>
</ref-list>
</back>
</article>