# The two rings of (50000) Quaoar

C. L. Pereira, B. Sicardy, B. E. Morgado, F. Braga-Ribas, E.

Fernández-Valenzuela, D. Souami, B. J. Holler, R. C. Boufleur, G. Margoti,

M. Assafin, et al.

► To cite this version:

C. L. Pereira, B. Sicardy, B. E. Morgado, F. Braga-Ribas, E. Fernández-Valenzuela, et al. The two rings of (50000) Quaoar. Astronomy & Astrophysics, 2023, 673, 10.1051/0004-6361/202346365. insu-04479090

HAL Id: insu-04479090

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

Submitted on 5 Mar 2024

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
LETTER TO THE EDITOR

The two rings of (50000) Quaoar

C. L. Pereira¹,² ID, B. Sicardy³, B. E. Morgado⁴,¹,² ID, F. Braga-Ribas⁵,¹,² ID, E. Fernández-Valenzuela⁶,⁷
D. Souami³,⁸,⁹,¹⁰,\* ID, B. J. Holler¹¹ ID, R. C. Boufleur¹,², G. Margoti⁵ ID, M. Assafin⁴,² ID, J. L. Ortiz⁷ ID,
P. Santos-Sanz⁷ ID, B. Epinat¹²,¹³ ID, P. Kervella³ ID, J. Desmars¹⁴,¹⁵ ID, R. Vieira-Martins¹,² ID, Y. Kilic¹⁶,¹⁷ ID,
A. R. Gomes Júnior¹⁸,¹⁹,² ID, J. I. B. Camargo¹,² ID, M. Emilio²⁰,¹,⁵, M. Vara-Lubiano⁷, M. Kretlow⁷,²¹,²²,
L. Albert²³,²⁴ ID, C. Alcock²⁵ ID, J. G. Ball²⁶, K. Bender²⁷, M. W. Buie²⁸, K. Butterfield²⁹, M. Camarca³⁰,
J. H. Castro-Chacón³¹ ID, R. Dunford²⁹, R. S. Fisher³² ID, D. Gamble²⁹,³³, J. C. Geary²⁵, C. L. Gnilka³⁴,
K. D. Green³⁵ ID, Z. D. Hartman²⁶, C.-K. Huang³⁶ ID, H. Januszewski¹³, J. Johnston³⁷, M. Kagitani³⁸, R. Kamin²⁹,
J. J. Kavelaars³⁹, J. M. Keller³⁷, K. R. de Kleer³⁰ ID, M. J. Lehner³⁶ ID, A. Luken³², F. Marchis⁴⁰,⁴¹, T. Marlin³⁰,
K. McGregor⁴² ID, V. Nikitin²⁹,⁴³, R. Nolthenius²⁷, C. Patrick²⁹, S. Redfield⁴² ID, A. W. Rengstorf⁴⁴ ID,
M. Reyes-Ruiz³¹ ID, T. Seccull²⁶ ID, M. F. Skrutskie⁴⁵, A. B. Smith²⁶, M. Sproul⁴⁶ ID, A. W. Stephens²⁶ ID,
A. Szentgyorgyi²⁵ ID, S. Sánchez-Sanjuán³¹ ID, E. Tatsumi⁴⁷,⁴⁸ ID, A. Verbiscer⁴⁵ ID, S.-Y. Wang³⁶, F. Yoshida⁴⁹,⁵⁰ ID,
R. Young⁵¹, and Z.-W. Zhang³⁶

(Affiliations can be found after the references)

Received 9 March 2023 / Accepted 17 April 2023

ABSTRACT

**Context.** Quaoar is a classical trans-Neptunian object (TNO) with an area-equivalent diameter of 1100 km and an orbital semi-major axis of 43.3 astronomical units. Based on stellar occultations observed between 2018 and 2021, an inhomogeneous ring (Q1R, i.e., Quaoar's first ring) has been detected around this body.

**Aims.** A new stellar occultation by Quaoar was observed on August 9, 2022, with the aim of improving Quaoar's shape models and the physical parameters of Q1R, while searching for additional material around the body.

**Methods.** The occultation provided nine effective chords across Quaoar, pinning down its size, shape, and astrometric position. Large facilities, such as Gemini North and the Canada-France-Hawaii Telescope (CFHT), were used to obtain high acquisition rates and signal-to-noise ratios. The light curves were also used to characterize the Q1R ring (radial profiles and orbital elements).

**Results.** Quaoar's elliptical fit to the occultation chords yields the limb with an apparent semi-major axis of 579.5 ± 4.0 km, apparent oblateness of 0.12 ± 0.01, and area-equivalent radius of 543 ± 2 km. Quaoar's limb orientation is consistent with Q1R and Weywot orbiting in Quaoar's equatorial plane. The orbital radius of Q1R is refined to a value of 4057 ± 6 km. The radial opacity profile of the more opaque ring profile follows a Lorentzian shape that extends over 60 km, with a full width at half maximum (FWHM) of ~5 km and a peak normal optical depth of 0.4. Besides the secondary events related to the already reported rings, new secondary events detected during the August 2022 occultation in three different data sets are consistent with another ring around Quaoar with a radius of 2520 ± 20 km, assuming the ring is circular and co-planar with Q1R. This new ring has a typical width of 10 km and a normal optical depth of ~0.004. Just as Q1R, it also lies outside Quaoar's classical Roche limit.

**Key words.** methods: data analysis – methods: observational – techniques: photometric – Kuiper belt objects: individual: Quaoar – planets and satellites: rings

1. Introduction

In the last decade, three ring systems have been discovered
around minor bodies in the outer Solar System: the Centaur
Chariklo (Braga-Ribas et al. 2014), the dwarf planet Haumea
(Ortiz et al. 2017), and the trans-Neptunian object (TNO)
(50000) Quaoar (Morgado et al. 2023). Dense material has also
been detected around the Centaur Chiron (Ruprecht et al. 2015;
Ortiz et al. 2015; Sickafoose et al. 2020). However, the nature of
this material, namely, whether it is a permanent or transient ring
or a dust shell, is still a matter of debate.

Quaoar's ring, referred to as Q1R hereafter, was detected
during several stellar occultations observed between 2018 and

2021 (Morgado et al. 2023). Q1R has a radius of about 4100 km with significant azimuthal variations in the optical depth, ranging between 0.004 and 0.1–0.7, and in width, ranging from 5 km to 300 km. Like Chariklo’s and Haumea’s rings, Quaoar’s Q1R ring orbits close to the 1/3 spin-orbit resonance (SOR) with the central body, suggesting a link between this resonance and the ring (Salo et al. 2021; Sicardy et al. 2021; Morgado et al. 2023). Meanwhile, a unique property of Q1R is its location, which is far outside Quaoar’s classical Roche limit. This limit is estimated to be at 1780 km from the body center, assuming particles with a bulk density of ρ = 0.4 g cm⁻³. Outside the Roche limit, rings should accrete into satellites over timescales of less than 100 years (Kokubo et al. 2000; Takeda & Ida 2001). However, collisions more elastic than previously considered for Saturn’s ring may maintain a ring unaccreted at distances greater than the Roche limit (Morgado et al. 2023). The 6/1

<footnote>* Fulbright Visiting Scholar (2022–2023) at University of California, Berkeley.</footnote>
mean-motion resonance (MMR) with Quaoar's satellite, Weywot, may contribute to the confinement of the ring since the satellite's eccentricity creates an equilibrium region capable of concentrating ring material in an arc over a longitude interval (Morgado et al. 2023). The confinement of the rings may also be due to the presence of putative "shepherd" satellites.

This Letter presents observations of a stellar occultation by Quaoar that occurred on August 9, 2022. This event was observed in Hawaii, the continental US, and Mexico. The high image acquisition rate and high signal-to-noise-ratios (S/Ns) obtained at Mauna Kea with the Alopeke camera at the Gemini North 8.1-m telescope and Wide-field InfraRed Camera (WIRCam) at the Canada-France-Hawaii 3.6-m Telescope (CFHT) allowed the dense part of Q1R to be radially resolved and its optical depth to be probed in the *r', z' and Ks* bands. Furthermore, this event revealed a new ring around Quaoar, hereafter called the Q2R (Quaoar's second ring).

## 2. Prediction and observation

The occultation of August 9, 2022 was predicted within the framework of the European Research Council (ERC) Lucky Star project<fnref n="1" />. The campaign was managed by the Occultation Portal website<fnref n="2" /> described in Kilic et al. (2022). The occultation shadow crossed the continental United States, Mexico, and the Hawaiian archipelago with a geocentric shadow velocity of 17.6 km s⁻<fnref n="1" />. The shadow path and the sites involved in the observation campaign are shown in Appendix C.1. More details on the occulted star and Quaoar are provided in Table 1.

Table B.1 summarizes the observational circumstances. Among the stations, CFHT was equipped with the WIRCam (Puget et al. 2004) at 2.15 µm (Ks-band) with an 8.8 Hz acquisition rate. The nearby Gemini North telescope used the Alopeke instrument (Scott et al. 2021) to simultaneously record the event with the *z'(947 nm) and r'(620 nm) bands at a 10 Hz cadence. The images at the Tohoku University Haleakala Observatory (TUHO) in Hawaii and at the Transneptunian Automated Occultation Survey (TAOS II) in Baja California (Mexico) were taken with no filter to maximize the S/N. Q2R was not detected in these observations, likely due to insufficient S/N values. At the University of California, Santa Cruz (UCSC - California, US), the Dunrhomin Observatory (Colorado, US), Sommers-Bausch Observatory (Colorado, US), Nederland (Colorado, US), and a mobile station in Bonny Doon Eco Reserve (California, US), only the main body was recorded due to a low S/N.

## 3. Data analysis

The data sets were analyzed using the standard photometry procedures of the Platform for Reduction of Astronomical Images Automatically (PRAIA, Assafin et al. 2011). The flux of the occulted star was corrected by the flux of nearby reference stars to account for sky transparency variations and was normalized to the total star+Quaoar flux outside the occultation.

The Stellar Occultation Reduction and Analysis package (SORA, Gomes-Júnior et al. 2022) was used to derive the ingress and egress instants of the occultation by the main body, assuming an opaque body without atmosphere. More details are presented in Appendix A. The extremities of the resulting chords were finally used to fit an elliptical model to Quaoar's limb at the moment of the occultation.

<table><caption>**Table 1.** More details about the occulted star and Quaoar.</caption><thead><tr><th></th><th>Occulted star</th></tr></thead><tbody><tr><td>Epoch</td><td>2022-08-09 06:34:04 UTC</td></tr><tr><td>Source ID</td><td>Gaia DR3 4098214367441486592</td></tr><tr><td>Star position at epoch <sup>(1)</sup></td><td>α<sub>★</sub> = 18h21m42s86965 ± 0.2404 mas<br>δ<sub>★</sub> = -15°12′45′′9639 ± 0.2191 mas</td></tr><tr><td>Magnitudes <sup>(2)</sup></td><td>G = 15.3; RP = 14.1; BP = 16.8;<br>J = 12.0; H = 11.0; K = 10.7</td></tr><tr><td>Apparent diameter <sup>(3)</sup></td><td>0.44 mas/1.33 km</td></tr><tr><td>RUWE <sup>(4)</sup></td><td>0.959</td></tr><tr><td></td><td>(50000) Quaoar</td></tr><tr><td>Ephemeris version</td><td>NIMAv16 <sup>(5)</sup></td></tr><tr><td>Geocentric distance</td><td>41.983157 au</td></tr><tr><td>Apparent magnitude</td><td>V = 18.8</td></tr><tr><td>Mass <sup>(6)</sup></td><td>(1.2 ± 0.05) × 10<sup>21</sup> kg</td></tr><tr><td>Rotation period <sup>(7)</sup></td><td>17.6788 ± 0.0004 h</td></tr><tr><td></td><td>Weywot's orbital pole <sup>(6)</sup></td></tr><tr><td>RA</td><td>17h44m ± 40m</td></tr><tr><td>Dec</td><td>+50° ± 6°</td></tr></tbody></table>

**Notes.** (¹)The star position was taken from the Gaia Data Release 3 (GDR3) star catalog (Gaia Collaboration 2023) and is propagated to the event epoch using the formalism of Butkevich & Lindegren (2014) applied with SORA (Gomes-Júnior et al. 2022). (²)J, H, and K from the NOMAD catalog (Zacharias et al. 2004). (³)Limb-darkened angular diameter estimated using a fit of the spectral energy distribution of the star with Castelli & Kurucz (2003) atmosphere models, with the reddening definition of Fitzpatrick (1999). (⁴)Renormalized unit weight error (Lindegren et al. 2021). (⁵)Obtained from the Numerical Integration of the Motion of an Asteroid (NIMA, Desmars et al. 2015), based on astrometric positions derived from previous stellar occultations and the Minor Planet Center (MPC) database, available at https://lesia.obspm.fr/lucky-star/obj.php?p=958 (⁶)Obtained from the method of Vachier et al. (2012). (⁷)Assuming the double-peaked light curve from Ortiz et al. (2003).

In the case of rings, the events were fitted using the square box model described in Appendix E and pipelines based on the SORA package. This provides the radial width of the ring $W_r$ (assuming a circular ring), its normal opacity, $p_N$, and normal optical depth, $\tau_N$, from which the equivalent width, $E_p = W_r p_N$, and equivalent depth, $A_\tau = W_r \tau_N$, are derived.

In the case of the CFHT and Gemini data, the acquisition rate and S/N were high enough to provide resolved profiles of the dense part of Q1R. Then, $E_p$ and $A_\tau$ can be obtained directly from the numerical calculations of the integrals $A_\tau = \int_{W_r} (v_r \tau_N) dt$ and $E_p = \int_{W_r} (v_r p_N) dt$ (i.e., without using the square model described above), $v_r$ is the star velocity perpendicular to the ring, measured in the ring plane.

The dense part of Q1R has a radial profile reminiscent of the Saturnian F-ring core (Bosh et al. 2002). More precisely, it is consistent with a Lorentzian shape used here to provide the full width at half maximum (FWHM) of the profile and the event's timing, from which its location relative to Quaoar is obtained.

## 4. Results

### 4.1. Main body

The event on August 9, 2022 provided ten chords across Quaoar's main body. However, Gemini used a dual camera at two wavelengths, resulting in nine effective chords across the

<footnote n="1">¹ https://lesia.obspm.fr/lucky-star</footnote>

<footnote n="2">² https://occultation.tug.tubitak.gov.tr</footnote>
<table><caption>**Table 2.** Retrieved parameters of Quaoar's body and its rings.</caption><thead><tr><th colspan="2">Quaoar's limb fitting</th></tr></thead><tbody><tr><td>Apparent semi-major axis</td><td>a' = 579.5 ± 4.0 km</td></tr><tr><td>Apparent oblateness</td><td>ε' = 0.12 ± 0.01</td></tr><tr><td>Position angle</td><td>P = 345.2 ± 1.2 deg</td></tr><tr><td rowspan="2">Ephemeris offset</td><td>f<sub>c</sub> = -24.7 ± 2.0 km</td></tr><tr><td>g<sub>c</sub> = -13.9 ± 2.1 km</td></tr><tr><td>R<sub>equiv</sub></td><td>543 ± 2 km</td></tr><tr><td>χ<sup>2</sup> per degree of freedom</td><td>0.986</td></tr><tr><td>Quaoar's geocentric position (ICRS)</td><td>α = 18<sup>h</sup>21<sup>m</sup>42<sup>s</sup>8677703 ± 0.249 mas<br>δ = -15°12'45"829691 ± 0.230 mas</td></tr><tr><td>Q1R's preferred pole position (ICRS)<sup>(1)</sup></td><td>α<sub>P</sub> = 17<sup>h</sup>19<sup>m</sup>16<sup>s</sup>8 ± 55<sup>s</sup>2<br>δ<sub>P</sub> = +53°27' ± 18'</td></tr><tr><td colspan="2">Q1R's geometry</td></tr><tr><td>Radius</td><td>4,057 ± 6 km</td></tr><tr><td>Opening angle</td><td>B = -20.0 ± 0.3 deg</td></tr><tr><td>Position angle</td><td>P = 350.2 ± 0.2 deg</td></tr><tr><td>χ<sup>2</sup> per degree of freedom</td><td>1.633</td></tr><tr><td colspan="2">Q2R's geometry<sup>(2)</sup></td></tr><tr><td>Radius</td><td>2520 ± 20 km</td></tr></tbody></table>

**Notes.** Quaoar's limb, Quaoar's geocentric position, and Q1R's geometry are given at epoch August 9th, 2022, 06:34:03.560 UTC. The fits use the star position and the NIMAv16 ephemeris for Quaoar (Table 1). (1) The preferred solution is co-planar with Quaoar's equatorial plane and Weywot's orbit within 5 ± 7 deg. The complementary position α<sub>P</sub> = 10<sup>h</sup>00<sup>m</sup> ± 05<sup>m</sup> and δ<sub>P</sub> = +79°21' ± 08', is not considered here as it provides a less satisfactory fit to the Q1R detections and yields an inclination of 45 deg with respect to Weywot's orbit, see Sect. 4.2. (2) Assuming a circular ring co-planar with Q1R. The error bar for the Q2R ring radius may be underestimated since we only have two effective detections, and several parameters were assumed.

body, that is, N = 18 chord extremities with the timings listed in Table A.1.

We used the SORA package to fit an ellipse to those extremities projected in the sky plane. This ellipse is defined by M = 5 parameters and the fit has N - M = 13 degrees of freedom. The fitted parameters are the ephemeris offsets f<sub>c</sub> and g<sub>c</sub> to apply to Quaoar's center in right ascension and declination, respectively, the apparent semi-major axis a' of the limb, its apparent oblateness ε' = (a' - b')/a', where b' is the apparent semi-minor axis of the limb, and P is the position angle of b'. The elliptical fit also considers a σ<sub>model</sub>, which is the uncertainty of the ellipse associated with the existence of putative topographic features on the surface of Quaoar. From the methodology presented by (Johnson & McGetchin 1973), we estimated that Quaoar might support topographic features of about 5 km, given Quaoar's density ρ = 1.99 g cm<sup>-3</sup> and strength S = 0.0303 × 10<sup>9</sup> dynes cm<sup>-2</sup>, consistent with an ice-rich composition. The standard deviation of the observed radial residuals is ~8 km, providing an upper limit for topographic features consistent with this prediction. Considering the equivalent radius in area as 542 km, this corresponds to 1.5%. The fit results are listed in Table 2. The chords and the best-fitting elliptical limb are plotted in Fig. 1, with more details given in Fig. A.2.

## 4.2. Ring Q1R

The Q1R ring was detected on both sides of the main body in the Gemini and CFHT light curves, and their physical properties

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

**Fig. 1.** Representation of our results on Quaoar's shape (center) and the detection of the two rings Q1R (outer ring) and Q2R (inner ring). The red segments correspond to the 1-σ error bars on the particular events. The orbit of Q1R is determined from a simultaneous fit using the present work and previous detections of 2018, 2019, 2020, and 2021 reported by Morgado et al. (2023) (see Sect. 4.2). The solution for the orbit of the new Q2R ring assumes that this ring is co-planar and concentric with Q1R. The central part of the plot (occultation by the solid body) is enlarged in Fig. A.2. In yellow, we show the 1/3 SOR resonance with Quaoar, and in teal is the 5/7 SOR resonance (considering the double-peaked rotation period). The purple ellipse represents 6/1 MMR with Weywot, and the green ellipse presents the expected Roche limit, considering particles with a bulk density of ρ = 0.4 g cm<sup>-3</sup>. The arrow shows the star's motion relative to Quaoar. Note: the orbital radius of Weywot is about three times larger than that of Q1R, and thus it is not shown in this representation.

were determined as described in Sect. 3. Although the detections of the dense part of the Q1R in the Gemini and CFHT light curves indicate the presence of diffuse material around the ring, the quality of the TAOS II and TUHO light curves only allows for the narrow and dense part of the ring to be detected. Therefore, the parameters E<sub>p</sub> and A<sub>τ</sub> were determined by the square box model for these data.

Due to the low optical depth of the Q1R segment intercepted before the closest approach, only the light curves obtained at Gemini and CFHT have sufficient S/N for detection (see an example in Fig. 2 and more details in Figs. 3, D.1, E.1, and E.2). The Gemini detection in the r' bandpass seems more sharply defined than its counterpart in z', but this effect remains marginal considering the S/N. The CFHT light curve shows a drop simultaneously, but it may be affected by a seeing deterioration that compromises an accurate determination of the ring boundaries. Even so, the central times of the Gemini and CFHT detections are all consistent at the 1-σ level. The parameters of Q1R for all the detections are presented in Table E.1.

We have combined the data from this work and the previous observations to improve the Q1R orbital parameters, assuming a fixed ring pole orientation between 2018 and 2022, as per Morgado et al. (2023). We tested a range of pole orientations and ring radii using a χ<sup>2</sup> statistic, resulting in two complementary solutions presented in Table 2. The preferred solution has a χ<sup>2</sup> per degree of freedom χ<sup>2</sup><sub>pdf</sub> = 1.6 and shows a better agreement with observations than the mirror solution with χ<sup>2</sup><sub>pdf</sub> = 4.1. Moreover, the preferred solution is co-planar with Weywot's orbit to within 5 ± 7 deg, as expected from a primordial disk surrounding Quaoar that evolved into a ring and formed its satellite, while the mirror solution is inclined by 45 ± 7 deg with respect to Weywot's orbit.

The position angles of Quaoar's projected pole (345.2 ± 1.2 deg) and Q1R (350.2 ± 0.2 deg) are misaligned by ~5 deg (Table 2). Thus, our results are consistent with Q1R orbiting
<figure><img src="image_3.png" /></figure>

**Fig. 2.** Detection of Q1R, Q2R, and the main body in the Gemini (z') light curve. The normalized flux is plotted in black vs. time (UTC), while the best-fitting models are over-plotted in red. The shallow events caused by Q1R (ingress) and Q2R (ingress and egress) are fitted by a square model. A Lorentzian function fits the dense part of Q1R at egress; see Sect. 3. The spikes observed during the disappearance and the re-appearance of the star behind Quaoar stem from the diffraction by the sharp opaque limb of the body. Note: there is a variation in the flux close to -110 and 290 s, caused by sudden variations in the seeing, and they can be neglected. More detailed views of Q1R and Q2R obtained at other telescopes are displayed in Figs. 3, D.1, E.1, and E.2.

close to Quaoar's equatorial plane, assuming that the body is an oblate spheroid. Part of this misalignment could stem from the fact that Quaoar is a triaxial ellipsoid (or a body with a more complex shape). Hence, the position angle of Quaoar's limb does not necessarily coincide with the position angle of Quaoar's pole.

### 4.3. The discovery of a new ring around Quaoar

The unique photometric quality of the Gemini and CFHT data allowed for the detection of additional material around Quaoar (Fig. 2). These data sets reveal additional secondary events symmetrically located with respect to Quaoar (Fig. 3). Two such events are simultaneously detected in the Gemini z'-band and CFHT Ks-band light curves before the closest approach, with detections reaching around $5.5\sigma$ and $5.2\sigma$, respectively. Conversely, the light curves display simultaneous events after the closest approach with significant detections, standing at $5.7\sigma$, $3.7\sigma$, and $4.7\sigma$ for the Gemini z', Gemini r', and CFHT data sets, respectively. Assuming that the light curves have a normal distribution, the probability that individual points of equivalent width $E_p(3\sigma) > 12$ km occur randomly in each light curve is $p \approx 1.4 \times 10^{-3}$, with p approaching zero for values larger than $E_p$. While the two light curves of the Gemini instrument may be correlated as they were taken at the same telescope, the Gemini and CFHT data are independent in terms of fast-seeing fluctuations. Using Poisson statistics, we estimate that the probability that the simultaneous events in the Gemini and CFHT data occur randomly due to the seeing fluctuations is very low, with $p \approx 10^{-6}$.

Some dips in flux were observed in the vicinity of Q2R detection in the region after the closest approach (Fig. 3). These detections are at the detection limit of our data, so we cannot say whether these are due to seeing variations or additional structures. Furthermore, we do not observe counterparts of these dips in the region before the closest approach.

All these detections are consistent with a new circular ring (Q2R) co-planar and concentric with Q1R, orbiting at 2520 ±

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

**Fig. 3.** Detection of the Q2R ring in the Gemini and CFHT light curves. The data points are plotted in black, and the red lines are the square box model fits derived from the quantities listed in Table E.1. An arbitrary offset was applied in the vertical direction for clarity. The time axis is relative to August 9, 2022 at 06:34:49.26 UTC, the time of the closest approach of Mauna Kea to Quaoar's shadow center.

20 km from Quaoar (Table 2). The radial width, optical depth, equivalent width, and equivalent depth of Q2R are listed in Table E.1. They were obtained as described in Appendix E.

## 5. Discussion and conclusions

The stellar occultation of August 9th, 2022 provided nine effective occultation chords obtained in Hawaii, Mexico, and the continental United States. They constrained Quaoar's shape, providing an apparent semi-major axis of $a' = 579.5 \pm 4.0$ km, an apparent oblateness $\epsilon' = 0.12 \pm 0.01$, and an area-equivalent radius of $R_{\text{equiv}} = a'\sqrt{1-\epsilon'} = 543 \pm 2$ km. Using an absolute magnitude of $H_V = 2.73\pm0.06$ (Fornasier et al. 2013), this yields
a geometric albedo of $p_V = 0.124 \pm 0.006$. Our value of $R_{\text{equiv}}$ differs from that of Braga-Ribas et al. (2013), $555.0 \pm 2.5$ km, by about 12 km, that is, at the $4-\sigma$ level. This difference could be evidence that Quaoar is not an oblate spheroid, but rather a triaxial ellipsoid or a body with a more complex shape; alternatively, it may be caused by a change in the aspect-angle since 2011, as observed from Earth.

The continuous region and the dense part of Q1R were detected during this event. The dense part was radially resolved and showed a Lorentzian profile reminiscent of Saturn's F ring (Bosh et al. 2002) or Neptune's ring arcs (Nicholson et al. 1990; Sicardy et al. 1991). This contrasts with the sharp edges observed for Chariklo's main ring C1R (Bérard et al. 2017; Morgado et al. 2021). The dense part of Q1R is detected over a radial width of ~60 km with a peak optical depth of $\tau_N \approx 0.4$, an FWHM Lorentzian width of 5 km and an equivalent width $E_p$ of around 2 km (see Table E.1). The values presented here are more precise than (but consistent with) those published by Morgado et al. (2023).

The detections of the dense part of Q1R in 2021 imply a minimum arc length of 365 km, corresponding to an azimuthal extension of about 5.1 deg (Morgado et al. 2023). The detections from CFHT, Gemini, TUHO, and TAOS II in 2022 (Figs. E.1, E.2) suggest a minimum arc-length of 226 km or ~3.2 deg. Since 2011, we have obtained 19 cuts of the Q1R with a sufficient S/N to detect the densest region. The limited azimuthal extent of the two detections of the dense regions means that they could both be parts of the same arc-like structure. If this structure was detected two times among the 19 cuts, then its extent can be estimated using Poisson statistics. This analysis yields an arc length with a 70% chance of falling between 18 and 72 deg. In this case, all of the detections in either 2021 or 2022 would each be correlated samples of one part of this arc.

The more tenuous component of Q1R is radially resolved and shows no marked structures, being consistent with a square model within our S/N limits. The best light curves (Gemini and CFHT) provide a typical width of 80–100 km and a typical normal optical depth of 0.003 (and, thus, an equivalent width of ~0.3 km) for that component at the longitude where it was detected (Table E.1).

Our detections, combined with those reported by Morgado et al. (2023), improved the orbital elements of Q1R. They are consistent with a circular ring of radius $4057.2 \pm 5.8$ km (Table 2) corresponding to about $7.5 \times$ Quaoar's area-equivalent radius $R_{\text{equiv}}$. This value coincides with the 6/1 MMR with Weywot (4020 ± 60 km) and within $3-\sigma$ of the 1/3 SOR with Quaoar (4200 ± 60 km), considering the double-peaked rotation period of $17.6788 \pm 0.0004$ h (Ortiz et al. 2003).

Our preferred solution for Q1R's orbital pole is consistent with this ring being co-planar with Weywot's orbit. Moreover, the apparent semi-major axes of Quaoar's limb and Q1R's orbit are aligned at the $4-\sigma$ level, suggesting that Q1R lies in Quaoar's equatorial plane, as expected from a colliding ring system (Kokubo et al. 2000).

Our data reveal a new ring (Q2R) around Quaoar. The detections are consistent with a circular ring of radius $2520 \pm 20$ km co-planar with Q1R (Table 2, Figs. 3 and 1). It has a typical width of 10 km, an optical depth of about 0.004, and an equivalent depth of about 0.04 km (see accurate values in Table E.1). Although closer to Quaoar than Q1R, Q2R is at $4.6 \times$ Quaoar's radius, also outside Quaoar's Roche limit, which is estimated to be around 1780 km, assuming the ring particle density as $\rho = 0.4 \text{ g cm}^{-3}$ (Morgado et al. 2023).

Using previously determined values for Quaoar's mass and rotation period (Table 1), we derived a 5/7 SOR radius of

$2525 \pm 35$ km. This coincides with Q2R's radius, $2520 \pm 20$ km, to within the $1-\sigma$ error bars. Like the 1/3 SOR, the 5/7 SOR is a second-order resonance and, as such, may play an essential role in the confinement of Q2R. However, more solid determinations of Q2R's orbit and Quaoar's shape are required.

Table E.1 shows that differences in the equivalent widths of Q1R and Q2R are observed between the $z'$, $r'$, and $Ks$ filters. The significance and interpretation of these differences require more analysis and will be discussed in a forthcoming publication. Moreover, the high S/N obtained at Gemini and CFHT will be used to detect or put a stringent upper limit on a putative atmosphere. Finally, the comparison of past (and possibly future) results derived from multi-chord occultations by Quaoar's main body will be used to constrain its shape better. This will be important for better understanding the dynamics of Quaoar's ring system, particularly under the effect of spin-orbit resonances with the body.

*Acknowledgements.* C.L.P is thankful for the support of the CAPES and FAPERJ/DSC-10 (E26/204.141/2022). This work was carried out within the "Lucky Star" umbrella that agglomerates the efforts of the Paris, Granada, and Rio teams, funded by the European Research Council under the European Community's H2020 (ERC Grant Agreement No. 669416). This study was financed in part by the National Institute of Science and Technology of the e-Universe project (INCT do e-Universo, CNPq grant 465376/2014-2). This study was financed in part by CAPES - Finance Code 001. The following authors acknowledge the respective CNPq grants: B.E.M. 150612/2020-6; F.B.R. 314772/2020-0; R.V.M. 307368/2021-1; M.A. 427700/2018-3, 310683/2017-3, 473002/2013-2; J.I.B.C. 308150/2016-3, 305917/2019-6. R.C.B acknowledges the FAPERJ grant E26/202/125/2020. E.F.-V. acknowledges financial support from the Florida Space Institute and the Space Research Initiative. J.L.O., P.S.-S., M.V-L, and M.K. acknowledge financial support from the grant CEX2021-001131-S funded by MCIN/AEI/ 10.13039/501100011033, they also acknowledge the financial support by the Spanish grants PID2020-112789GB-100 from AEI and Proyecto de Excelencia de la Junta de Andalucía PY20-01309. Funding for RECON was provided by grants from USA: NSF AST-1413287, AST-1413072, AST-1848621, and AST-1212159. We thank RECON observers Doug Thompson, Ken Conway, Dorey Conway, Terry Miller, David Schulz, Michael von Schalscha, and Matt Christensen for their efforts in collecting data. Based on observations obtained with WIRCam, a joint project of CFHT, Taiwan, Korea, Canada, France, at the Canada-France-Hawaii Telescope (CFHT) which is operated from the summit of Maunakea by the National Research Council of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii. The observations at the Canada-France-Hawaii Telescope were performed with care and respect from the summit of Maunakea which is a significant cultural and historic site. We thank Marc Baril and Tom Vermeulen for their time dedicated to the observation performed at Canada-France-Hawaii Telescope (CFHT). Based on observations obtained at the international Gemini Observatory, a program of NSF's NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation on behalf of the Gemini Observatory partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Ministerio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovações e Comunicações (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). This work made use of data from GN-2022B-DD-101 observing program and were obtained with the High-Resolution Imaging instrument(s) 'Alopeke (and/or Zorro). 'Alopeke (and/or Zorro) was funded by the NASA Exoplanet Exploration Program and built at the NASA Ames Research Center by Steve B. Howell, Nic Scott, Elliott P. Horch, and Emmett Quigley. 'Alopeke (and/or Zorro) was mounted on the Gemini North (and/or South) telescope of the international Gemini Observatory, a program of NSF's NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. This work was enabled by observations made from the Gemini North telescope, located within the Maunakea Science Reserve and adjacent to the summit of Maunakea. We are grateful for the privilege of observing the Universe from a place that is unique in both its astronomical quality and its cultural significance. This research used SORA, a python package for stellar occultations reduction and analysis, developed with the support of ERC Lucky Star and LIneA/Brazil, within the collaboration of Rio-Paris-Granada teams. This work profited from unpublished occultations by Quaoar made at SOAR
(SO2019A-003) and the Pico dos Dias Observatory (OP2019A-004) to improve
the accuracy of the ephemeris NIMAv16.

References

Assafin, M., Vieira Martins, R., Camargo, J. I. B., et al. 2011, Gaia Follow-up Network for Solar System Objects, 85

Bérard, D., Sicardy, B., Camargo, J. I. B., et al. 2017, AJ, 154, 144

Bosh, A. S., Olkin, C. B., French, R. G., & Nicholson, P. D. 2002, Icarus, 157, 57

Braga-Ribas, F., Sicardy, B., Ortiz, J. L., et al. 2013, ApJ, 773, 26

Braga-Ribas, F., Sicardy, B., Ortiz, J. L., et al. 2014, Nature, 508, 72

Butkevich, A. G., & Lindegren, L. 2014, A&A, 570, A62

Castelli, F., & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, eds. N. Piskunov, W. W. Weiss, & D. F. Gray, 210, A20

Cuzzi, J. N. 1985, Icarus, 63, 312

Desmars, J., Camargo, J. I. B., Braga-Ribas, F., et al. 2015, A&A, 584, A96

Elliott, J. L., French, R. G., Meech, K. J., & Elias, J. H. 1984, AJ, 89, 1587

Fitzpatrick, E. L. 1999, PASP, 111, 63

Fornasier, S., Lellouch, E., Müller, T., et al. 2013, A&A, 555, A15

Gaia Collaboration (Vallenari, A., et al.) 2023, https://doi.org/10.1051/0004-6361/202243940

Gomes-Júnior, A. R., Morgado, B. E., Benedetti-Rossi, G., et al. 2022, MNRAS, 511, 1167

Johnson, T., & McGetchin, T. 1973, Icarus, 18, 612

Kilic, Y., Braga-Ribas, F., Kaplan, M., et al. 2022, MNRAS, 515, 1346

Kokubo, E., Ida, S., & Makino, J. 2000, Icarus, 148, 419

Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2

Morgado, B., Benedetti-Rossi, G., Gomes-Júnior, A. R., et al. 2019, A&A, 626, L4

Morgado, B. E., Sicardy, B., Braga-Ribas, F., et al. 2021, A&A, 652, A141

Morgado, B. E., Sicardy, B., Braga-Ribas, F., et al. 2023, Nature, 614, 239

Nicholson, P. D., Cooke, M. L., Matthews, K., Elias, J. H., & Gilmore, G. 1990, Icarus, 87, 1

Ortiz, J. L., Gutiérrez, P. J., Sota, A., Casanova, V., & Teixeira, V. R. 2003, A&A, 409, L13

Ortiz, J. L., Duffard, R., Pinilla-Alonso, N., et al. 2015, A&A, 576, A18

Ortiz, J. L., Santos-Sanz, P., Sicardy, B., et al. 2017, Nature, 550, 219

Puget, P., Stadler, E., Doyon, R., et al. 2004, in Ground-based Instrumentation for Astronomy, eds. A. F. M. Moorwood, & M. Iye, SPIE Conf. Ser., 5492, 978

Roques, F., Moncuquet, M., & Sicardy, B. 1987, AJ, 93, 1549

Ruprecht, J. D., Bosh, A. S., Person, M. J., et al. 2015, Icarus, 252, 271

Salo, H., Sicardy, B., Mondino-Llermanos, A., et al. 2021, European Planetary Science Congress, EPSC2021-338

Scott, N. J., Howell, S. B., Gnilka, C. L., et al. 2021, Front. Astron. Space Sci., 8, 138

Sicardy, B., Roques, F., & Brahic, A. 1991, Icarus, 89, 220

Sicardy, B., Salo, H., Souami, D., et al. 2021, European Planetary Science Congress, EPSC2021-91

Sickafoose, A. A., Bosh, A. S., Emery, J. P., et al. 2020, MNRAS, 491, 3643

Takeda, T., & Ida, S. 2001, ApJ, 560, 514

Vachier, F., Berthier, J., & Marchis, F. 2012, A&A, 543, A68

Zacharias, N., Monet, D. G., Levine, S. E., et al. 2004, Am. Astron. Soc. Meet. Abstr., 205, 48.15

7 Instituto de Astrofísica de Andalucía – Consejo Superior de Investigaciones Científicas, Glorieta de la Astronomía S/N, 18008 Granada, Spain

8 Departments of Astronomy, and of Earth and Planetary Science, 501, Campbell Hall, University of California, Berkeley, CA 94720, USA

9 naXys, Department of Mathematics, University of Namur, Rue de Bruxelles 61, 5000 Namur, Belgium

10 Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Laboratoire Lagrange, Bd de l'Observatoire, CS 34229, 06304 Nice Cedex 4, France

11 Space Telescope Science Institute, Baltimore, MD, USA

12 Aix Marseille Université, CNRS, CNES, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, 13388 Marseille, France

13 Canada-France-Hawaii Telescope, 65-1238 Mamalahoa Highway, Kamuela, HI 96743, USA

14 Institut Polytechnique des Sciences Avancées IPSA, 63 boulevard de Brandebourg, 94200 Ivry-sur-Seine, France

15 Institut de Mécanique Céleste et de Calcul des Éphémérides, IMCCE, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC Univ Paris 06, Univ. Lille, France

16 Akdeniz University, Faculty of Sciences, Department of Space Sciences and Technologies, 07058 Antalya, Turkey

17 TÜBİTAK National Observatory, Akdeniz University Campus, 07058 Antalya, Turkey

18 Institute of Physics, Federal University of Uberlândia, Uberlândia-MG, Brazil

19 UNESP – São Paulo State University, Grupo de Dinâmica Orbital e Planetologia, CEP 12516-410, Guaratinguetá, SP, Brazil

20 Universidade Estadual de Ponta Grossa, O.A. – DEGEO, Ponta Grossa, PR, Brazil

21 Internationale Amateursternwarte (IAS) e. V., Mittelstr. 6, 15749 Mittenwalde, Germany

22 International Occultation Timing Association – European Section (IOTA/ES), Am Brombeerhag 13, 30459 Hannover, Germany

23 Département de Physique et Observatoire du Mont-Mégantic, Université de Montréal, C.P. 6128, Succ. Centre-ville, Montréal, H3C 3J7 Québec, Canada

24 Institut Trottier de Recherche sur les exoplanètes, Université de Montréal, Montreal, Canada

25 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA

26 Gemini Observatory/NSF's NOIRLab, Hilo, HI, USA

27 Cabrillo College Astronomy, Aptos, CA, USA

28 Southwest Research Institute, 1050 Walnut Street, Suite 300, Boulder, CO 80302, USA

29 International Occultation Timing Association (IOTA), Topeka, KS, USA

30 California Institute of Technology, 1200 E California Blvd, Pasadena, CA 91125, USA

31 Universidad Nacional Autónoma de México, Instituto de Astronomía, AP 106, Ensenada 22800, BC, Mexico

32 University of Oregon, Eugene, OR, USA

33 CanCON - Canadian Collaborative Occultation Network, Canada

34 NASA Ames Research Center, Moffett Field, California, USA & NASA Exoplanet Science Institute, Caltech/IPAC, Mail Code 100-22, Pasadena, CA, USA

35 University of New Haven, Dept. of Mathematics and Physics, 300 Boston Post Road, West Haven, CT 06477, USA

36 Academia Sinica Institute of Astronomy and Astrophysics, 11F of AS/NTU Astronomy-Mathematics Building, No.1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan, ROC

37 University of Colorado, 2000 Colorado Avenue, Boulder, CO 80309, USA

38 Planetary Plasma and Atmospheric Research Center, Graduate School of Science, Tohoku University, 6-3 Aramaki -aza-aoba, Aoba-ku, Sendai 980-8578, Japan

39 Department of Physics and Astronomy, University of Victoria, Building, 3800 Finnerty Road, Victoria, BC V8P 5C2, Canada

<footnote n="1">¹ Observatório Nacional/MCTI, R. General José Cristino 77, CEP 20921-400 Rio de Janeiro, RJ, Brazil
e-mail: chrystianpereira@on.br</footnote>

<footnote n="2">² Laboratório Interinstitucional de e-Astronomia – LIneA, Rio de Janeiro, RJ, Brazil</footnote>

<footnote n="3">³ LESIA, Observatoire de Paris, Université PSL, Sorbonne Université, Université de Paris, CNRS, 92190 Meudon, France</footnote>

<footnote n="4">⁴ Universidade Federal do Rio de Janeiro – Observatório do Valongo, Ladeira Pedro Antônio 43, CEP 20.080-090 Rio de Janeiro, RJ, Brazil</footnote>

<footnote n="5">⁵ Federal University of Technology – Paraná (UTFPR-Curitiba), Rua Sete de Setembro, 3165, CEP 80230-901 Curitiba, PR, Brazil</footnote>

<footnote n="6">⁶ Florida Space Institute, UCF, 12354 Research Parkway, Partnership 1 Building, Room 211, Orlando, FL, USA</footnote>
40 Unistellar, 5 allée Marcel Leclerc, bâtiment B, 13008 Marseille, France

41 SETI Institute, Carl Sagan Center, Suite 200, 339 Bernardo Avenue, Mountain View, CA 94043, USA

42 Astronomy Department and Van Vleck Observatory, Wesleyan University, Middletown, CT 06459, USA

43 Research and Education Collaborative Occultation Network, USA

44 Purdue University Northwest, Department of Chemistry and Physics, Hammond, IN, USA

45 University of Virginia, Department of Astronomy, PO Box 400325, Charlottesville, VA 22904, USA

46 Private Observatory, PA, USA

47 Instituto de Astrofísica de Canarias, University of La Laguna, Tenerife, Spain

48 Dept. Earth and Planetary Science, The University of Tokyo, Tokyo, Japan

49 School of Medicine Department of Basic Sciences University of Occupational and Environmental Health, 1-1 Iseigaoka, Yahata, Kitakyusyu 807-8555, Japan

50 Planetary Exploration Research Center, Chiba Institute of Technology, 2-17-1 Tsudanuma, Narashino, Chiba 275-0016, Japan

51 Naylor Observatory, Lewisberry, PA, USA
# Appendix A: Main body occultation and elliptical fit

The light curves obtained with the Gemini and CFHT telescopes present prominent diffraction signatures in the main body occultation mainly due to the high acquisition rate, the large wavelength, and high S/N of the light curves. During the occultation, Quaoar was at a geocentric distance of 41.983157 au, implying a Fresnel scale, $F = \sqrt{\lambda D/2}$, from 1.5 to 2.6 km, for visible to near-infrared (CFHT Ks-band) wavelengths, respectively. Considering the apparent velocity of the event as 17.57 km s⁻¹ and the acquisition rate of the Gemini and CFHT instruments, we obtained a spatial resolution of these light curves of ~1.8 km. These high-S/N light curves have a Fresnel diffraction effect on the same order of magnitude as the instrumental response (e.g., the integration time).

<figure><img src="image_5.png" /><figcaption>**Fig. A.1.** Gemini z' light curve (black dots) and the modeled light curve (red) considering the local normal star velocity, the wavelength and apparent star size. The geometric model is in blue. The normalized flux is plotted as a function of the seconds before and after the local closest approach (2022-08-09 06:34:49.26 UT).</figcaption></figure>

With the limb-darkened angular diameter of the star estimated and considering the observation band pass and its FWHM, we also fit the local normal velocity ($v_⊥$). By minimizing the $\chi^2$ statistic, we found the velocity and the ingress (egress) time that best explained the observed light curve. This was applied to the Gemini z'-band data since this was our best light curve. The radial velocities obtained are $v_⊥ = 10.2$ km s⁻¹ and $v_⊥ = 16.6$ km s⁻¹ for ingress and egress instants, respectively. The separation between the chords obtained from the Gemini and CFHT when projected onto the sky plane was about 160 meters, a value smaller than the star's apparent diameter. Therefore, the velocities were considered the same at both sites. This procedure resulted in an optimal fit for both diffraction spikes and the baseline occultation by the main body. Figure A.1 presents an example of this fit with the Gemini z' light curve plotted with the modeled curve.

Ingress and egress times for the main body occultation are presented in Table A.1. The instant uncertainties of the Gemini

<table><caption>**Table A.1.** Ingress and egress times for Quaoar's main body.</caption><thead><tr><th>Site</th><th>Ingress<br>(hh:mm:ss.s)</th><th>Egress<br>(hh:mm:ss.s)</th></tr></thead><tbody><tr><td>CFHT</td><td>06:34:23.34 (0.02)</td><td>06:35:11.42 (0.01)</td></tr><tr><td>Gemini (z')</td><td>06:34:23.29 (0.01)</td><td>06:35:11.29 (0.01)</td></tr><tr><td>Gemini (r')</td><td>06:34:23.27 (0.03)</td><td>06:35:11.28 (0.02)</td></tr><tr><td>TUHO</td><td>06:34:23.5 (0.2)</td><td>06:35:18.1 (0.2)</td></tr><tr><td>TAOS II</td><td>06:30:38.3 (0.2)</td><td>06:31:35.1 (0.1)</td></tr><tr><td>Dunrhomin</td><td>06:30:03.3 (0.4)</td><td>06:31:02.8 (0.9)</td></tr><tr><td>Sommers-Bausch</td><td>06:30:02.6 (0.2)</td><td>06:31:03.8 (0.5)</td></tr><tr><td>UCSC</td><td>06:31:20.3 (0.1)</td><td>06:32:00.9 (0.1)</td></tr><tr><td>Nederland</td><td>06:30:03.4 (0.5)</td><td>06:31:04.4 (0.8)</td></tr><tr><td>Bonny Doon</td><td>06:31:20.7 (0.4)</td><td>06:32:01.3 (0.4)</td></tr></tbody></table>

**Notes.** Times are in UT and the error bars in parentheses are in seconds and given at the 1-$\sigma$ level. The significant differences in the CFHT and Gemini instants are due to the length of the chords, with the CFHT being ~1 km longer than the Gemini, as expected from the limb.

<figure><img src="image_6.png" /><figcaption>**Fig. A.2.** Close-up view of Fig. 1 shows the best elliptical fit (black ellipse) and ellipses within 1-$\sigma$ region (grey) to the ten chords derived from the timings in Table A.1. The grey 1-$\sigma$ fits differ only slightly from the black best-fit. The red segments represent the 1-$\sigma$ error bars in the ingress and egress times. At this scale, the Gemini and CFHT chords are superimposed, and their error bars would be too small to be seen (i.e., they are here represented by red squares, much bigger than their actual sizes).</figcaption></figure>

and CFHT light curves are limited by occultation modeling, not data quality. The elliptical fit parameters to the chord extremities are presented in Fig. A.2.
## Appendix B: Observational circumstances

The observational circumstances for this occultation campaign are presented in Table B.1. We present in this work the analysis of ten data sets with positive detections for the main body, of which four are detections of the Q1R ring and three are detections of both the Q1R and Q2R rings. In addition, we have more

eighteen sites that participated in the observation campaign, but reported cloudy skies. The target star was too faint for the setup in Woodside (CA, US) and Summerland BC, CA, even with 4 and 3 seconds of exposure, respectively; in addition, the Summerland data was acquired while clouds were drifting through the field.

<table><caption>Table B.1. Observational circumstances.</caption><thead><tr><th>Site</th><th>Longitude [° '"]<br>Latitude [° '"]<br>Altitude [m]</th><th>Observers</th><th>Aperture (mm)<br>Detector<br>Filter</th><th>Exp. Time<br>Cycle<br>(s)</th><th>Status<br>Light curve<br>Dispersion</th></tr></thead><tbody><tr><td>CFHT</td><td>-155 28 07.93</td><td>H. Januszewski</td><td>3,580</td><td>0.113</td><td>Positive</td></tr><tr><td>MaunaKea</td><td>19 49 31.01</td><td>M. Baril, B. Epinat</td><td>WIRCam</td><td>0.113</td><td>0.013</td></tr><tr><td>HI, US</td><td>4,206.1</td><td>T. Vermeulen</td><td>Ks</td><td></td><td></td></tr><tr><td>Gemini-North</td><td>-155 28 08.6</td><td>T. Seccull</td><td>8,100</td><td>0.100</td><td>Positive</td></tr><tr><td>MaunaKea</td><td>19 49 26</td><td>A. Stephens</td><td>'Alopeke<br>r' & z'</td><td>0.101</td><td>z': 0.012<br>r': 0.022</td></tr><tr><td>TUHO</td><td>-156 15 27</td><td>E. Tatsumi</td><td>600</td><td>1.500</td><td>Positive</td></tr><tr><td>Haleakala</td><td>20 42 25.14</td><td>F. Yoshida</td><td>ASI178MM</td><td>1.634</td><td>0.098</td></tr><tr><td>HI, US</td><td>3,051.8</td><td>M. Kagitani</td><td>None</td><td></td><td></td></tr><tr><td>TAOS II</td><td>-115 27 47.88</td><td>J. H. Castro-Chacón</td><td>1,300</td><td>0.200</td><td>Positive</td></tr><tr><td>San Pedro Martir</td><td>31 02 36.26</td><td>S. Sánchez-Sanjuán</td><td>Andor Ixon 888</td><td>0.205</td><td>0.274</td></tr><tr><td>BC, MX</td><td>2,824.8</td><td></td><td>None</td><td></td><td></td></tr><tr><td>Dunrhomin Obs.</td><td>-105 09 46.57</td><td></td><td>280</td><td>2</td><td>Positive</td></tr><tr><td>Longmont</td><td>40 15 08.46</td><td>M. Buie</td><td>QHY174M-GPS</td><td>2</td><td>0.229</td></tr><tr><td>CO, US</td><td>1,592.8</td><td></td><td>None</td><td></td><td></td></tr><tr><td>Sommers-Bausch Obs.</td><td>-105 15 47.09</td><td>J. Keller</td><td>508</td><td>0.9</td><td>Positive</td></tr><tr><td>Boulder</td><td>40 00 13.32</td><td>J. Johnston</td><td>QHY174M-GPS</td><td>0.9</td><td>0.259</td></tr><tr><td>CO, US</td><td>1,645</td><td></td><td>None</td><td></td><td></td></tr><tr><td>Nederland</td><td>-105 26 44.05</td><td>M. Skrutskie</td><td>200</td><td>2</td><td>Positive</td></tr><tr><td>CO, US</td><td>39 59 13.95</td><td>A. Verbicer</td><td>QHY174M-GPS<br>570 nm</td><td>2</td><td>0.313</td></tr><tr><td>Bonny Doon Eco Reserve</td><td>-122 08 18.62</td><td></td><td>203</td><td>0.260</td><td>Positive</td></tr><tr><td>Santa Cruz</td><td>37 03 03.12</td><td>K. Bender</td><td>Watec 910HX/RC</td><td>0.260</td><td>0.538</td></tr><tr><td>CA, US</td><td>488.4</td><td></td><td>None</td><td></td><td></td></tr><tr><td>UCSC</td><td>-122 04 46.20</td><td></td><td>203</td><td>0.500</td><td>Positive</td></tr><tr><td>Santa Cruz</td><td>37 01 03.47</td><td>R. Nolthenius</td><td>Watec 910HX/RC</td><td>0.534</td><td>0.556</td></tr><tr><td>CA, US</td><td>344.0</td><td></td><td>None</td><td></td><td></td></tr><tr><td>Summerland</td><td>-119 40 12</td><td></td><td>457</td><td>3</td><td>Negative</td></tr><tr><td>BC, CA</td><td>49 36 00</td><td>D. Gamble</td><td>QHY174M-GPS<br>None</td><td>3</td><td></td></tr><tr><td>Woodside</td><td>485</td><td></td><td>114</td><td>4</td><td>Inconclusive</td></tr><tr><td>CA, US</td><td>37 23 36.2</td><td>F. Marchis</td><td>eVscope/IMX347</td><td>4</td><td></td></tr><tr><td></td><td>367</td><td></td><td>None</td><td></td><td></td></tr><tr><td>Naylor Obs.</td><td>-76 53 43.48</td><td>R. Kamin</td><td>355.6</td><td>0.5</td><td>Overcast</td></tr><tr><td>Harrisburg</td><td>40 08 55.60</td><td>R. Young</td><td>QHY174M-GPS</td><td>0.5</td><td></td></tr><tr><td>PA, US</td><td>181.6</td><td></td><td></td><td></td><td></td></tr><tr><td>Nikitin Home</td><td>-105 11 00.96</td><td></td><td>279.4</td><td>6</td><td>Overcast</td></tr><tr><td>Gunbarrel</td><td>40 04 08.61</td><td>V. Nikitin</td><td>QHY174-GPS<br>None</td><td>6</td><td></td></tr><tr><td>CO, US</td><td>1,594</td><td></td><td></td><td></td><td></td></tr><tr><td>Van Vleck Obs.</td><td>-72 39 33.12</td><td>K. McGregor</td><td>610</td><td></td><td>Overcast</td></tr><tr><td>Middletown</td><td>41 33 18.62</td><td>S. Redfield</td><td>Apogee E2V<br>Luminance</td><td>No data</td><td></td></tr><tr><td>CT, US</td><td>68.9</td><td></td><td></td><td></td><td></td></tr><tr><td>Skychariot</td><td>-74 58 49.28</td><td></td><td>406.4</td><td></td><td>Overcast</td></tr><tr><td>Shohola</td><td>41 21 39.94</td><td>M. Sproul</td><td>ASI1600MC<br>None</td><td>No data</td><td></td></tr><tr><td>PA, US</td><td>414.6</td><td></td><td></td><td></td><td></td></tr><tr><td>NIRo - CAC</td><td>-87 22 31.08</td><td></td><td>508</td><td></td><td>Overcast</td></tr><tr><td>Lowell</td><td>41 16 15.03</td><td>A. W. Rengstorf</td><td>FLI PL09000<br>I</td><td>No data</td><td></td></tr><tr><td>IN, US</td><td>195.7</td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>**Table B.1.** [continued.] Observational circumstances.</caption><thead><tr><th>Site</th><th>Longitude [° '"]<br>Latitude [° '"]<br>Altitude [m]</th><th>Observers</th><th>Aperture (mm)<br>Detector<br>Filter</th><th>Exp. Time<br>Cycle<br>(s)</th><th>Status<br>Light curve<br>Dispersion</th></tr></thead><tbody><tr><td>Jimginny Obs.<br>Naperville<br>IL, US</td><td>-88 07 0.00<br>41 45 32.40<br>230</td><td>R. Dunford</td><td>356<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Pine Mountain Obs.<br>Bend<br>OR, US</td><td>-120 56 28.8<br>43 47 30.3<br>1,920</td><td>S. Fisher<br>A. Luken</td><td>610<br>Andor Xyla<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Palomar Obs.<br>San Diego<br>CA, US</td><td>-116 51 53.5<br>33 21 21.5<br>1,712</td><td>E. Fernandez-Valenzuela<br>K. de Kleer<br>T. Marlin, M. Camarca</td><td>5,100<br>CHIMERA</td><td>No data</td><td>Overcast</td></tr><tr><td>Westport Astronomical Obs.<br>Westport<br>CT, US</td><td>-73 19 39.0<br>41 10 16.0<br>87</td><td>K. D. Green</td><td>356<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Yerington<br>NV, US</td><td>-119 09 39<br>38 59 28<br>1,340</td><td>M. Christensen</td><td>280<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Bend<br>OR, US</td><td>-121 18 55<br>44 03 29<br>1,105</td><td>A-M. Eklund</td><td>280<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Yuma<br>CA, US</td><td>-114 37 40<br>32 41 34<br>62</td><td>D. Thompson<br>K. Conway<br>D. Conway</td><td>280<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Reno<br>NV, US</td><td>-119 49 06<br>39 32 42<br>1,387</td><td>T. Stoffel</td><td>280<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Cedarville<br>CA, US</td><td>-120 10 24<br>41 31 45<br>1,420</td><td>T. Miller<br>D. Schulz</td><td>280<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Kingman<br>AZ, US</td><td>-113 51 46<br>35 24 43<br>966</td><td>C. Patrick<br>K. Butterfield</td><td>305<br>QHY174M-GPS<br>0.5 focal reducer</td><td>No data</td><td>Overcast</td></tr><tr><td>Beatty<br>CA, US</td><td>-116 45 29<br>36 54 04<br>999</td><td>J. Heller</td><td>280<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Burney<br>CA, US</td><td>-121 23 56<br>41 02 45<br>1,012</td><td>M. von Schalscha</td><td>280<br>QHY174M-GPS<br>None</td><td>No data</td><td>Overcast</td></tr><tr><td>Gardnerville<br>NV, US</td><td>-119 44 59<br>38 56 29<br>1,449</td><td>J. Bardecker</td><td>304<br>Watec 910-HX<br>None</td><td>No data</td><td>Overcast</td></tr></tbody></table>

**Notes.** TAOS II: The Transneptunian Automated Occultation Survey. CFHT: Canada-France-Hawaii Telescope. TUHO: Tohoku University Haleakala Observatory. UCSC: University of California, Santa Cruz. NIRO: Northwest Indiana Robotic - Calumet Astronomy Center.
**Appendix C: Occultation map**

Figure C.1 presents the stellar occultation map with all the sites involved in this campaign. The green dots mark the positive detections and orange dots are stations that reported cloudy

weather. The white dot represents the Woodside observation with the eVscope where the event was not detected due to the low S/N of the data set. The continuous line represents the limit of the Quaoar's shadow projected over the Earth. The dashed-lines represent the Q1R and Q2R projections, as annotated in the image.

<figure><img src="image_7.png" /><figcaption>**Fig. C.1.** Occultation map with positive detections (black dots) and sites that reported cloudy weather (orange dots). The white dot represents no data due to the instrumental limits. In red we display a site with data but no occultation detected. The black solid lines mark Quaoar's shadow limits and the black dashed lines mark the Q1R and Q2R projections. The black arrow in the bottom right indicates the direction of motion of the shadow.</figcaption></figure>
## Appendix D: Light curves

Figure D.1 displays all the normalized light curves obtained during the stellar occultation by (50000) Quaoar of August 9, 2022. The stellar flux has been corrected for sky transparency fluctuations using reference stars. The time scale is relative to the closest approach time at each site. Besides the occultation by Quaoar and Q1R, the Gemini and CFHT light curves revealed the new Q2R ring. The Tohoku (TUHO) and TAOS II data allowed for the detection of the densest region of the Q1R ring. Due to the S/N limitations, the Nederland, Dunrhomin, Sommers-Bausch, UCSC, and Bonny Doon stations detected only the occultation

by the main body. Some data sets presented dropped frames during the photometric process, such as the Durhomin and Nederland. This occurred due to the degradation of the sky quality shortly after the occultation of the main body, in addition to clouds crossing the field of view of the telescope.

We noticed that the Gemini and CFHT light curves were misaligned in time, probably caused by an offset in the Gemini data. As we have a reliable time source for the CFHT data, we aligned the centers of these chords by applying an offset of +0.27 seconds on the Gemini chord, obtaining values of $\chi^2_{pdf}$ closer to 1. Although very close, this result is slightly better than applying an offset to the CFHT chord.

<figure><img src="image_8.png" /><figcaption>**Fig. D.1.** All the positive light curves obtained during the August 9, 2022 stellar occultation by Quaoar. The black dots represent the data points. These light curves are plotted as a function of the time in seconds relative to the closest approach for each site. The green and gray vertical lines stand for Q2R and Q1R, respectively. When detected, these lines stand for calculated detection times. For light curves where these secondary structures were not detected, the green and gray lines mark theoretical times expected from our best-fit solution for a circular ring. The horizontal gray dashed lines represent the baseline and minima of the stellar fluxes.</figcaption></figure>
## Appendix E: Modeled ring detections

The ring events were fit using the models described in Elliot et al. (1984) and Bérard et al. (2017), that is, square boxes with uniform opacity and sharp edges. As with the main body occultation modeling, the modeled curve was convolved with observation bandpass, apparent star diameter, and instrumental response, but now fitting an opacity for the square box. The outputs of a given fit are the width, W, of the event, its apparent opacity, p' (corresponding to the depth of the observed stellar flux drop), and its apparent optical depth, τ' = -ln(1 - p'). The geometry of the ring (assumed here to be circular) is defined by its opening angle, B, and the position angle, P, of its apparent semi-minor axis, counted positively from local celestial north to the celestial east direction. This provides the radial width, W<sub>r</sub>, the normal opacity, p<sub>N</sub> = |sin B| (1 - √1 - p'), and the normal optical depth, τ<sub>N</sub> = τ' |sin B|/2, of the ring from each event. The factor of 2 in the formula for τ<sub>N</sub> is due to the fact that Airy diffraction by individual ring particles results in a loss of light, resulting in an observed optical depth to be twice as large as would be measured at the ring level; see more details in Cuzzi (1985), Roques et al. (1987). We note that the formulae for p<sub>N</sub> and τ<sub>N</sub> are valid only for mono-layer and poly-layer rings, respectively. The equivalent width (respectively, depth)

E<sub>p</sub> = W<sub>r</sub>p<sub>N</sub> (respectively, A<sub>τ</sub> = W<sub>r</sub>τ<sub>N</sub>) measures the amount of ring material that blocked the stellar rays in the mono-layer (resp. poly-layer) case.

A Lorentzian function is fitted to the dense part of the Q1R in the Gemini and CFHT light curves, converted from flux to normal optical depth as a function of radial distance in the ring plane. The area under the curve equals the equivalent width, A<sub>τ</sub>, of the ring, and the full-width half minimum (FWHM) gives us the approximate width of the ring's core. The position of the function's valley gives us the average radial distance between the ring and Quaoar's center. The E<sub>p</sub> value was obtained from the integral of the ring profile in the curve of the normal opacity, p<sub>N</sub>, as a function of the radial distance in the sky plane. This modeling considers the stellar apparent diameter of 1.33 km, which has negligible influence on the ring profile.

Figs. E.1 and E.2 present all light curves in which the Q1R and Q2R rings were detected. The TAOS II light curve, used for detecting of the densest region of the Q1R ring, was obtained by applying aperture photometry to stacked images in order to increase the S/N. Stacking was performed using Python routines built from the `astropy` library as described in Morgado et al. (2019), where each new image is the average of six original images, resulting in a temporal resolution of 1.2 seconds for each stacked image.

<table><caption>**Table E.1.** Physical parameters of rings Q1R and Q2R.</caption><thead><tr><th>Ring</th><th>Detection</th><th>Mid-time<br>August 9<sup>th</sup>, 2022</th><th>r<br>(km)</th><th>W<sub>r</sub><br>(km)</th><th>p<sub>N</sub></th><th>τ<sub>N</sub></th><th>E<sub>p</sub><br>(km)</th><th>A<sub>τ</sub><br>(km)</th></tr></thead><tbody><tr><td rowspan="3">Q1R-ing</td><td>CFHT</td><td>06:31:15.8 (0.2)</td><td>4,009.9 (1.8)</td><td>105.5 (5.6)</td><td>0.0030 (0.0003)</td><td>0.0030 (0.0003)</td><td>0.31 (0.02)</td><td>0.31 (0.02)</td></tr><tr><td>Gemini (z')</td><td>06:31:15.72 (0.04)</td><td>3,995.4 (0.3)</td><td>76.4 (0.8)</td><td>0.0026 (0.0001)</td><td>0.0026 (0.0001)</td><td>0.20 (0.01)</td><td>0.20 (0.01)</td></tr><tr><td>Gemini (r')</td><td>06:31:15.7 (0.1)</td><td>3,995.9 (0.8)</td><td>85.8 (1.1)</td><td>0.0032 (0.003)</td><td>0.0032 (0.0003)</td><td>0.28 (0.02)</td><td>0.28 (0.02)</td></tr><tr><td rowspan="5">Q1R-egr</td><td>CFHT<sup>a</sup></td><td>06:38:32.581 (0.003)</td><td>4,123.69 (0.05)</td><td>6.1 (0.1)</td><td>0.330 (0.007)</td><td>0.395 (0.009)</td><td>2.01 (0.03)</td><td>2.41 (0.04)</td></tr><tr><td>Gemini (z')<sup>a</sup></td><td>06:38:32.444 (0.003)</td><td>4,123.11 (0.02)</td><td>5.3 (0.1)</td><td>0.353 (0.008)</td><td>0.417 (0.007)</td><td>1.87 (0.02)</td><td>2.21 (0.02)</td></tr><tr><td>Gemini (r')<sup>a</sup></td><td>06:38:32.429 (0.005)</td><td>4,122.85 (0.04)</td><td>5.1 (0.1)</td><td>0.341 (0.009)</td><td>0.398 (0.008)</td><td>1.74 (0.03)</td><td>2.03 (0.03)</td></tr><tr><td>TUHO</td><td>06:38:37.4 (0.7)</td><td>4,107.2 (7.6)</td><td>71 (17)</td><td>0.03 (0.01)</td><td>0.03 (0.01)</td><td>2.3 (0.9)</td><td>2.4 (1.1)</td></tr><tr><td>TAOS II</td><td>06:34:55.8 (0.5)</td><td>4,125.7 (4.4)</td><td>34 (16)</td><td>0.05 (0.02)</td><td>0.05 (0.03)</td><td>1.6 (1.1)</td><td>1.7 (1.3)</td></tr><tr><td rowspan="3">Q2R-ing</td><td>CFHT</td><td>06:32:43.9 (0.2)</td><td>2,490.4 (1.6)</td><td>16.5 (3.0)</td><td>0.0022 (0.0004)</td><td>0.0022 (0.0004)</td><td>0.036 (0.006)</td><td>0.036 (0.006)</td></tr><tr><td>Gemini (z')</td><td>06:32:43.61 (0.06)</td><td>2,493.8 (0.6)</td><td>11.3 (0.8)</td><td>0.0034 (0.0005)</td><td>0.0034 (0.0005)</td><td>0.038 (0.006)</td><td>0.039 (0.006)</td></tr><tr><td>Gemini (r')</td><td>n.d.<sup>b</sup></td><td>n.d.<sup>b</sup></td><td>n.d.<sup>b</sup></td><td>n.d.<sup>b</sup></td><td>n.d.<sup>b</sup></td><td>n.d.<sup>b</sup></td><td>n.d.<sup>b</sup></td></tr><tr><td rowspan="3">Q2R-egr</td><td>CFHT</td><td>06:37:00.20 (0.01)</td><td>2,540.7 (0.1)</td><td>10.8 (0.2)</td><td>0.0049 (0.0003)</td><td>0.0049 (0.0003)</td><td>0.052 (0.001)</td><td>0.053 (0.001)</td></tr><tr><td>Gemini (z')</td><td>06:37:00.09 (0.02)</td><td>2,540.45 (0.1)</td><td>6.3 (0.3)</td><td>0.0047 (0.0005)</td><td>0.0048 (0.0006)</td><td>0.030 (0.002)</td><td>0.030 (0.002)</td></tr><tr><td>Gemini (r')</td><td>06:36:59.93 (0.02)</td><td>2,528.3 (1.1)</td><td>9.1 (2.5)</td><td>0.003 (0.001)</td><td>0.003 (0.001)</td><td>0.030 (0.008)</td><td>0.030 (0.008)</td></tr></tbody></table>

**Notes.** The normal opacity, p<sub>N</sub>, and normal optical depth, τ<sub>N</sub>, were calculated from the ring opening angle, B, and position angle, P, on August 9, 2022, derived from the orientation of the body obtained from the 2018-2022 data. The other parameters are: Mid-time in UT and the error bars in parentheses given in seconds; r, the radial distance from Quaoar's center in kilometers; W<sub>r</sub>, the radial width in kilometers; E<sub>p</sub> and A<sub>τ</sub>, the equivalent width and equivalent depth in kilometers respectively. The terms "ing" and "egr" stand for ingress and egress, respectively, and refer to the fact that the detection occurred before and after the occultation by Quaoar's main body, respectively. The error bars in parentheses are at the 1-σ level. <sup>a</sup>From Lorentzian fit. The width W<sub>r</sub> is defined as the FWHM of the τ<sub>N</sub> profile. The E<sub>p</sub> and A<sub>τ</sub> values where obtained from respective integrals in the ring profile (see Sec. 3); <sup>b</sup>Not detected.
<figure><img src="image_9.png" /><figcaption>**Fig. E.1.** Fits to the Q1R and Q2R light curve data taken with Gemini (z'), Gemini (r'), and CFHT (Ks) (the corresponding filters are indicated in parentheses). These light curves are plotted as a function of the time in seconds relative to the local closest approach (C/A). The blue dots represent residuals with an arbitrary vertical offset for clarity. Note: the y-axis scale is different in the rightmost panel in each row (dense part of Q1R).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. E.2.** Fits to the Q1R dense ring in light curve data taken with TUHO and TAOS II telescopes, from top to bottom. TAOS II data are stacked every six images. These light curves are plotted as a function of the time in seconds relative to the local closest approach (C/A), with a different y-axis scale in the rightmost panel in each row (dense part of Q1R). The green vertical dashed lines represent the theoretical times for the Q1R and Q2R rings. Note: although marginal, the two data points that detect the Q1R ring in the TAOS II light curve are statistically significant at the 4.2-σ level.</figcaption></figure>