*2.2. Specular Point Positions and Ray Points*

The specular point positions and the ray tracing of the direct, incident, and reflected signals are calculated based on the methodology presented in [11,20]. A geometrical model is employed to characterize specular reflections and determine the specular point position, considering the Earth's surface curvature. For this model, the transmitter (*Tx*) and receiver (*Rx*) positions are needed in an Earth-centered Earth-fixed (ECEF) frame. The *Rx* position is extracted from the Spire data files. To obtain this position, the Lemur-2 satellites are equipped with a zenith dual-frequency (L1 and L2) antenna, which facilitates precise orbit

determination (POD). The *Tx* position is derived from the broadcasted GPS ephemeris. The Earth's curvature is modeled with an osculating spherical surface with respect to the WGS-84 ellipsoid at a reference specular point. An iterative solution is employed to find the best-fitting sphere that satisfies the condition of equal incident and reflected angles (specular reflection) [20]. The specular point positions are calculated at 10 s intervals on the receiver trajectory. A ray-tracing module is set to compute ray points every 10 km along the three ray paths: *Tx* to *SP* (incident), *SP* to *Rx* (reflected), and *Tx* to *Rx* (direct). The positions of the ray points (latitude, longitude, and ellipsoidal height) are subsequently utilized to obtain the electron density from the ionospheric electron density models. Figure 3 illustrates an example of the electron density retrieval from the NEDM2020 model depicting the change along the specular point tracks every 10 s (blue stars), and the ray points change every 10 km (red dots) along the incident (in), reflected (re), and direct (dr) ray paths.

**Figure 3.** Representation of the specular point change along track (blue stars) for the SVN 79 and PRN 30 every 10 s (~45 km) and the ray points every 10 km along the direct, incident, and reflected ray paths (red dots).

Following the ray tracing, a total of 28,790 reflection events are obtained. The total number of reflection events by region is depicted in Figure 4a. Additionally, Figure 4b illustrates the distribution of reflection events concerning the elevation angle by region. Notably, the south pole region exhibits a higher number of events; however, all regions show similar behavior, with a higher concentration of events in the elevation range between 5◦ and 20◦.

**Figure 4.** (**a**) Total number of reflection events on 1 March 2021 by region. (**b**) The number of reflection events across grazing angles, ranging from 5◦ to 30◦.
