Open Access
Issue
A&A
Volume 680, December 2023
Article Number A9
Number of page(s) 10
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202347477
Published online 05 December 2023

© The Authors 2023

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

Pluto’s atmosphere was discovered during the 1985 stellar occultation (Brosch 1995), and since then, stellar occultations have played a crucial role in studying its structure, composition, and evolution over time (Hubbard et al. 1988; Elliot et al. 1989, 2003; Yelle & Elliot 1997; Sicardy et al. 2003, 2011, 2016, 2021; Pasachoff et al. 2005, 2017; Young et al. 2008, 2021; Rannou & Durry 2009; Person et al. 2013, 2021; Olkin et al. 2015; Bosh et al. 2015; Gulbis et al. 2015; Dias-Oliveira et al. 2015; Meza et al. 2019; Arimatsu et al. 2020). A compilation of 12 occultations observed between 1988 and 2016 revealed a three-fold monotonic increase in the atmospheric pressure of Pluto during that period (Meza et al. 2019). This increase can be explained by the volatile transport model (VTM) of the Laboratoire de Météorologie Dynamique (LMD; Bertrand & Forget 2016; Forget et al. 2017; Bertrand et al. 2018, 2019), which was subsequently fine-tuned by Meza et al. (2019). This model provides a framework for simulating the volatile cycles on Pluto over both seasonal and astronomical timescales, allowing us to explore the long-term evolution of Pluto’s atmosphere and its response to seasonal variations over its 248-yr heliocentric orbital period (Meza et al. 2019). According to the LMD VTM in Meza et al. (2019, VTM19 hereafter), Pluto’s atmospheric pressure is expected to have reached its peak around the year 2020. The pressure increase is attributed to the progression of summer over the northern hemisphere of Pluto, exposing Sputnik Planitia (SP)1 to solar radiation. The surface of SP, which is composed of nitrogen (N2), methane (CH4), and carbon monoxide (CO) ices, is believed to sublimate and release volatile gases into the atmosphere during this period, leading to a pressure increase. After reaching its peak, the model predicts a gradual decline in pressure over the next two centuries under the combined effects of Pluto’s recession from the Sun and the prevalence of the winter season over SP.

On one hand, the VTM19 remains consistent with the analysis of Sicardy et al. (2021) of the 6 June 2020 occultation observed at Devasthal, where two colocated telescopes were used. This latter analysis suggests that Pluto’s atmosphere has been in a plateau phase since mid-2015, which aligns with the model predictions that the atmospheric pressure reached its peak around 2020.

On the other hand, the Arimatsu et al. (2020) analysis of the 17 July 2019 occultation observed by a single telescope (TUHO) suggests a rapid pressure decrease between 2016 and 2019. These authors detected a significant pressure drop at the 2.4σ level. However, it is worth noting that the geometry of this occultation is grazing. This may have introduced larger correlations between the pressure and the geocentric closest approach distance to Pluto’s shadow axis, leading to insufficient precision to confidently support the claim of a large pressure decrease followed by a return in 2020 to a pressure level close to that of 2015 (Sicardy et al. 2021).

These contrasting results highlight the need for occultation observations between 2019 and 2020 in order to better understand the behavior and evolution of Pluto’s atmosphere during this time period. Furthermore, while Young et al. (2021) support the presence of a pressure drop based on their analysis of the 15 August 2018 occultation, Sicardy et al. (2021) suggest that careful comparisons between measurements by independent teams should be made before drawing any conclusions on the pressure evolution.

Observations of the 5 September 2019 occultation, which have not been reported by other teams, are presented in Sect. 2, followed by a description of the light-curve fitting methods in Sect. 3. These unique observations allow us to track the changes in Pluto’s atmosphere during the time period between the events studied by Arimatsu et al. (2020) and Sicardy et al. (2021). Results are detailed in Sect. 4, and the pressure evolution is discussed in Sect. 5, including comparisons with the reanalyzed 15 August 2018 and 17 July 2019 events. Conclusions and recommendations are provided in Sect. 6.

2 Occultation observations

Two observation campaigns were organized in China for occultations in 2019 (see Appendix A). One occurred on 17 July 2019, which was studied by Arimatsu et al. (2020), and the other on 5 September 2019, which is reported in the present paper for the first time. Due to bad weather conditions in many areas, no effective light curves were observed by our stations for the first occultation, and only two light curves were obtained for the second.

Table 1 lists the circumstances of the 5 September 2019 event. Figure 1 presents all the observation stations and the reconstructed path of the shadow of Pluto2 during this event. Table A.2 lists the circumstances of stations with positive detections. Their station codes are DWM and HNU.

To ensure accurate and precise timing in stellar occultations, some stations (e.g., DWM as shown in Table A.2) were equipped with QHY174GPS cameras. These cameras, manufactured by QHYCCD3, offer precise recording of observation time and location for each frame using a GPS-based function, and have been used in many stellar occultation studies (e.g., Buie et al. 2020a,b; Morgado et al. 2021, 2022; Pereira et al. 2023). In the light-curve fitting procedures described in Sect. 3.2, the time-recording offsets of the QHY174GPS cameras are fixed to zero, considering their reliability and accuracy as time references.

All observational data were captured in the FITS format. These data were processed using the Tangra occultation photometric tool4 (Pavlov 2020) and our data-reduction code (see Appendix B). It was ensured that the targets and reference stars in all the images we used were not overexposed. The resulting light curves from the observations, after being normalized, are presented in Fig. 2. Each data point on the light curves is represented by fi(t) ± σi(t), where i indicates the quantities associated with a specific station, t represents the recorded timing of each frame, f the normalized total observed flux of the occulted star and the Pluto’s system, and σ the measurement error associated with each data point.

Table 1

Circumstances and light-curve fitting results of the 5 September 2019 event.

thumbnail Fig. 1

Reconstructed occultation map of the 5 September 2019 event.

3 Light-curve fitting methods

3.1 Light-curve model

In order to simulate observed light curves, we implemented a light-curve model, ϕ(t; A, s, Δt, Aτ, Δρ, p0), which is described in Appendix C and is consistent with DO15 (Dias-Oliveira et al. 2015; Sicardy et al. 2016, 2021; Meza et al. 2019). As a function of model parameters, its time-dependent Jacobian matrix was also implemented to represent the sensitivity of the model to the corresponding parameters to be estimated through fitting procedures.

The light-curve model of a given station can be formally written as ϕi(t)=Ai(siψi(t;Δti,Δτ,Δρ,p0)+(1si)),${\phi _i}(t) = {A_i} \cdot \left( {{s_i} \cdot {\psi _i}\left( {t;\Delta {t_i},\Delta \tau ,\Delta \rho ,{p_0}} \right) + \left( {1 - {s_i}} \right)} \right),$(1)

where i indicates the quantities associated with the station; for further details, the reader is referred to Appendix C. Here, the reference ephemerides we use are the NIMAv95 asteroidal ephemeris (Desmars et al. 2015, 2019) for the orbit of the Pluto system barycenter with respect to the Sun, the PLU0586 satellite ephemerides (Brozovic et al. 2015; Jacobson et al. 2019) for the orbit of Pluto with respect to the Pluto system barycenter, and the DE4408 planetary ephemerides (Park et al. 2021) for the orbits of the Earth and the Sun with respect to the Solar System barycenter. The reference star catalog where the data of the occulted star are obtained is Gaia DR3.

3.2 Fitting procedure

The light-curve model was fitted to the normalized observed light curves simultaneously by nonlinear least squares, returning a χ2-type value of goodness-of-fit. The goal is to minimize the objective function given by χobs2=i,j(ϕi(tij)fi(tij))2σi2(tij),$\chi _{{\rm{obs}}}^2 = \sum\limits_{i,j} {{{{{\left( {{\phi _i}\left( {{t_{ij}}} \right) - {f_i}\left( {{t_{ij}}} \right)} \right)}^2}} \over {\sigma _i^2\left( {{t_{ij}}} \right)}}}, $(2)

where tij represents the mid-exposure time of the jth observation of the station i.

In addition, with the used reference ephemerides and star catalog, some a priori information on Δρ can be obtained: Δρ=0 km±σρ,$\Delta \rho = 0{\rm{km}} \pm {\sigma _\rho },$(3)

where the uncertainty σρ is set to 72 km using the positional uncertainties listed in the “orbit quality” table of NIMAv9 and in Gaia DR3. This σρ value corresponds to about 3 mas on the sky at the geocentric distance of Pluto. The a priori information can be treated as independent observational data and used in the model fitting, with the objective function modified as: χapr2=χobs2+Δρ2σρ2.$\chi _{{\rm{apr}}}^2 = \chi _{{\rm{obs}}}^2 + {{\Delta {\rho ^2}} \over {\sigma _\rho ^2}}.$(4)

The fitting steps are as follows:

  • In order to find all local minima at which a nonlinear least-squares fitting could potentially get stuck, we explored the two-parameter space (Δρ, p0) by generating the variation of χ2 as a function of them. Figure 3 presents such two χ2 maps, labeled (a) and (b), which are analyzed in Sect. 4. The maps are generated by minimizing χobs2$\chi _{{\rm{obs}}}^2$ or χapr2$\chi _{{\rm{apr}}}^2$ at each fixed (Δρ, p0) point on a regularly spaced grid. The Levenberg-Marquardt (LM) method, which is implemented in the LMFIT package8, was used in each fitting procedure. The free parameters to be adjusted are Δτ9, Δti of any station with no reliable time reference system like QHY1 4GPS, and si and Ai of each station.

  • For a more accurate best-fitting solution for (Δρ, p0), the LM method is used again, with Δρ and p0 adjusted with initial guesses located at all known local minima of each χ2 map.

  • Each χ2 map, which provides information about the quality of the fit, is used to define confidence limits based on constant χ2 boundaries (Press et al. 2007).

4 Results

Figure 3a shows the χobs2$\chi _{{\rm{obs}}}^2$ map for the 5 September 2019 occultation. Two local minima are observed. However, considering the significant χ2 difference of 9 between the two local minima, the global minimum is more likely to be the correct solution. In addition, the Δρ value at the global minimum is more consistent with the NIMAv9 solution, Δρ = 0 km, at the 0.16 σρ level, compared with the other local one at the 2.44 σρ level.

In an effort to mitigate or at least further weaken the presence of multiple local minima, we calculated the χapr2$\chi _{{\rm{apr}}}^2$ map by adding the χ2-type value of the a priori information, (Δρ/σρ)2, into the χobs2$\chi _{{\rm{obs}}}^2$ map. Figure 3b presents the results, which show that two local minima are still present, but with a χ2 difference of about 14.5, which is larger than that of the χobs2$\chi _{{\rm{obs}}}^2$ map. Therefore, the global minimum is confidently accepted as the solution for (ρcag, psurf), as provided in Table 1.

Moreover, Fig. 3 presents the consistency of our derived psurf across the two different local minima. Our findings demonstrate that the specific choice of local minima does not significantly affect the value of psurf, further supporting the reliability of our solution for psurf.

thumbnail Fig. 2

Occultation observations and the best-fitting light-curve model of the 5 September 2019 event. Panel a: observed and simultaneously fitted light curves. Panels b and c: reconstructed stellar paths seen by DWM and HNU, respectively.

5 Pressure evolution

5.1 Comparisons and necessary reanalyses of historical events

In Fig. 4, the red plot represents our psurf measurement from the 5 September 2019 occultation. We also include other published measurements (Hinson et al. 2017; Meza et al. 2019; Arimatsu et al. 2020; Young et al. 2021; Sicardy et al. 2021) and the pressure evolution predicted by the VTM19 in order to provide a comprehensive view of the pressure variations on Pluto.

To avoid potential inconsistencies arising from different analysis methods, as discussed by Sicardy et al. (2021), we reanalyzed the 15 August 2018 event studied by Young et al. (2021) using the IXON observational data of Silva-Cabrera et al. (2022). The derived pressure measurement presented in Appendix D.1 is 12.0270.08+0.09$12.027_{ - 0.08}^{ + 0.09}$ µbar. In addition, we also reanalyze the 17 July 2019 event in Appendix D.2, deriving a pressure of psurf=9.4210.75+0.68${p_{{\rm{surf}}}} = 9.421_{ - 0.75}^{ + 0.68}$ µbar, which is similar to that of Arimatsu et al. (2020), of 9.560.34+0.52$9.56_{ - 0.34}^{ + 0.52}$ µbar. This similarity is expected because the same DO15 method is used. As this same method is used by Meza et al. (2019) and Sicardy et al. (2021), their pressure measurements, along with that of Arimatsu et al. (2020), can be fully compared with our new ones. Both the remeasurements are plotted in black in Fig. 4.

The pressure measurement from the 5 September 2019 event shows alignments with those from the 19 July 2016, 15 August 2018, and 6 June 2020 events within their combined 1σ levels. Our new measurement from the 15 August 2018 event does not show the significant pressure drop previously reported by Young et al. (2021). The previously reported pressure drop between the 19 July 2016 and 17 July 2019 events is still detected at the same level as in Arimatsu et al. (2020).

5.2 Discussion on pressure variations

While the VTM19 remains, on average, applicable and capable of predicting the main atmospheric behavior during the observed years, there are also two probable V-shaped pressure variations observed from 2010 to 2015 and from 2015 to 2020, especially when considering the measurements from the 4 June 2011 and 17 July 2019 events. These V-shaped variations suggest the presence of additional factors that have not been accounted for. Specifically, short-term changes in Pluto’s surface ices and their interaction with the atmosphere are likely contributing to the variation. Moreover, spectral monitoring of the surface composition has revealed some short-term changes in the ices over several Earth years (e.g., Grundy et al. 2014; Lellouch et al. 2022; Holler et al. 2022).

However, the validity of the V-shaped variations is debatable given the stellar faintness of the 4 June 2011 event and the grazing single-chord geometry of the 17 July 2019 event. If the debatable measurement from the 17 July 2019 event were discarded, no significant changes would be observed between 2016 and 2020. This more likely supports the plateau phase since 2015 predicted by the VTM19. In order to better understand the relationship between these factors, further observations using multiple observational techniques (occultation, spectroscopy, and photometry) are required, as well as simulations with a refined VTM.

thumbnail Fig. 3

The χ2 maps of the 5 September 2019 event. Panels a and b: The χobs2$\chi _{{\rm{obs}}}^2$ and χapr2$\chi _{{\rm{apr}}}^2$ maps, respectively. The χobs2$\chi _{{\rm{obs}}}^2$ denotes the goodness-of-fit only using observational data, while χapr2$\chi _{{\rm{apr}}}^2$ denotes the goodness-of-flt with the additional a priori information represented as the χ2 type value, (Δρρ)2. These maps are used to derive the best-fitting atmospheric pressure p0 at the reference radius r0 of 1215 km and the cross-track correction Δρ to the ephemerides, of which the a priori uncertainty σρ is 72 km. The surface pressure psurf and the geocentric closest approach distance ρcag to the shadow center are obtained by linear transformations of Δρ and p0, respectively. The best-fitting χobs2$\chi _{{\rm{obs}}}^2$ and χapr2$\chi _{{\rm{apr}}}^2$ values per degree of freedom are 1.160 and 1.155, respectively.

thumbnail Fig. 4

Pressure evolution of Pluto over a 248-yr heliocentric orbital period predicted by the VTM19, along with the measured pressures with 1σ error bars. The 3σ error bars of our new measurements and of the previous measurements of the 19 July 2016, 17 July 2019, and 6 June 2020 events are also presented using the published χ2 maps. AN2${A_{{{\rm{N}}_2}}}$ denotes the albedo of nitrogen ice.

6 Conclusions

The unique observations of the 5 September 2019 occultation provide a surface pressure of psurf = 11.478 ± 0.55 µbar. In order to avoid potential method inconsistencies in comparing with historical pressure measurements (Sicardy et al. 2021), we also reanalyzed the 15 August 2018 and 17 July 2019 events based on publicly available data (Silva-Cabrera et al. 2022; Arimatsu et al. 2020). All measurements are presented in Fig. 4.

The VTM 19 remains applicable on average. In addition, we also observed unaccounted-for V-shaped pressure variations with the previously reported pressure drop being a part of these variations; however, these variations are debatable. To better understand all significant pressure variations of Pluto, continuous pressure monitoring through occultation observations is essential where possible. Also, simultaneous and frequent spectroscopic and photometric monitoring of changes to its surface ice are important, as such comprehensive monitoring will provide more short-term and long-term evolution constraints of Pluto’s interacting atmosphere and surface.

Acknowledgements

We acknowledge Bruno Sicardy for his useful comments that helped improving this manuscript. This work has been supported by the National Natural Science Foundation of China (Grant Nos. 12203105 and 12103091). We acknowledge the science research grants from the China Manned Space Project with NO.CMS-CSST-2021-A12 and No. CMS-CSST-2021-B10. This research has made use of data from the 40 cm DOB telescope at the Dawei Mountain Observatory of Hunan Astronomical Association and from the 50 cm RC telescope at the Observatory of Hebei Normal University. We ackonwledge the support of Chinese amateur astronomers from Hunan Astronomical Association, Nanjing Amateur Astronomers Association, and Shenzhen Astronomical Observatory Team.

Appendix A Two occultation campaigns in 2019

Table A.1

Stations that encountered weather problems on 17 July 2019

Table A.2

Stations with positive detections on 5 September 2019

Table A.3

Stations that encountered weather problems on 5 September 2019

The 17 July 201910 and the 5 September 201911 stellar occultations by Pluto were originally predicted by the ERC Lucky Star project12.

A.1 The 17 July 2019 occultation campaign

Figure A.1a presents the reconstructed path of the shadow of Pluto13 during the 17 July 2019 event, which was studied by Arimatsu et al. (2020). We also organized an observation campaign involving stations detailed in Table A.1. All the above-mentioned stations are also presented in Figure A.1a. Unfortunately, all our stations encountered weather problems.

A.2 The 5 September 2019 occultation campaign

The occultation map of the 5 September 2019 event is given by Figure 1. Table A.2 lists the circumstances of stations with positive detections. Table A.3 lists the circumstances of stations that encountered weather problems.

thumbnail Fig. A.1

Reanalysis of the 17 July 2019 event based on the TUHO observational data from Figure 1 of Arimatsu et al. (2020) (credit: Arimatsu et al, A&A, 638, L5, 2020, reproduced with permission © ESO). Panel (a): Reconstructed occultation map. Panel (b): Observed and simultaneously fitted light curves. Panel (c): Reconstructed stellar paths seen by TUHO. Panel (d): The χ2 map, where χ2 denotes the goodness-of-fit value using these data. The best-fitting χ2 value per degree of freedom is 0.860.

Appendix B Data processing

The observational data were initially obtained in the FITS format, a commonly used format for astronomical data. To process and analyze these data, we used the Tangra occultation photometric tool developed by Pavlov (2020). This tool offers various functionalities and algorithms specifically designed for analyzing occultation events. When calibration images including bias, dark, and flat-field frames are available, they are used for image correction. The default measurement type and tracking method provided by Tangra were applied. The reduction method was aperture photometry with median background subtraction. The combined image of the occulted star and Pluto’s system is circled as the target.

For each station, a selection of suitable guiding and nearby reference stars (typically three) was made. All images of target and reference stars are confirmed as unaffected by overexposure. The signal flux counts, denoted as Isignal, and the background flux counts, denoted as Ibkg, were measured for the target and each reference star in all FITS files. The timing information and measurement results were exported to a CSV file specific to each station, from which the observed light curves would be generated.

Further data processing and analysis were performed using our data reduction code implemented in Python. Each flux measurement for the target or a reference star was calculated as I = IsignalIbkg. The measurement error, σI, was modeled using the classical form given by: σI=σbkg2+Igeff1,${\sigma _I} = \sqrt {\sigma _{{\rm{bkg}}}^2 + I \cdot g_{{\rm{eff}}}^{ - 1}} ,$(B.1)

where σbkg is the standard deviation of the background (Ibkg) and geff is the effective gain that represents the number of electrons per flux count.

The relative photometric result of the target is defined as F=ITIRS,$F = {{{I_T}} \over {{I_{RS}}}},$(B.2)

where IT is the I of target and IRS the sum of those of all reference stars. This step helps eliminate low-frequency variations caused by atmospheric and instrumental effects. The corresponding measurement error, σF, was derived using the usual formula for propagation of errors (Press et al. 2007).

If the effective gain geff was not precisely known, geff1$g_{eff}^{ - 1}$ was adjusted with the non-negative constraint to fit the derived σF model (σF(geff1)${\sigma _F}\left( {g_{{\rm{eff}}}^{ - 1}} \right)$ where geff10$g_{{\rm{eff}}}^{ - 1} \ge 0$) to the standard deviation of F outside the occultation part.

Finally, the normalized total observed flux, denoted as f, and its error, σ, were derived by dividing F and σF, respectively, by the median of F outside the occultation part. As a result, f outside the occultation part approximates unity.

Appendix C Synthetic light curve model for a stellar occultation by Pluto

For a given station, the light-curve model is represented as: ϕ(t)=A(sψ(t)+(1s)),$\phi (t) = A \cdot (s \cdot \psi (t) + (1 - s)),$(C.1)

where t is the recorded timing of the observation and ϕ is the received flux, as a function of t; A is the total flux of the star and Pluto’s system when the star is not occulted; s is the flux ratio of the occulted star to A; and ψ is the normalized (between zero and unity) flux of the occulted star, which represents the total flux of the primary (sometimes called the near-limb) and secondary (or far-limb, when available) images produced by Pluto’s spherical planetary atmosphere (Sicardy 2023) as shown in Figure 2.

To calculate ψ, the following equation is used: ψ(t)=| r+z || dr+dz |+| rz || drdz |,(z>0,r±>Rp)$\psi (t) = \left| {{{{r_ + }} \over z}} \right|\left| {{{d{r_ + }} \over {dz}}} \right| + \left| {{{{r_ - }} \over z}} \right|\left| {{{d{r_ - }} \over {dz}}} \right|,\quad \left( {z > 0,{r_ \pm } > {R_{\rm{p}}}} \right)$(C.2)

where z is the distance of the station to the shadow axis, that is, the line passing through the center of Pluto and parallel to the light ray from the star; r+ and r represent the closest approach distances of the primary and secondary images to the shadow axis, respectively; and Rp is the radius of Pluto’s body. The dividers, |r±/z|, and derivatives, |dr±/dz|, account for the compression and stretching effects produced by the spherical planetary atmosphere, respectively, caused by limb curvature and differential refraction.

Given reference ephemerides and a star catalog, the distance z is modeled in the International Celestial Reference Frame (ICRF) using the following equation: z(t)=| (ρp(t+Δti+Δτ)+Δρn^cagρi(t+Δti))×^s |,$z(t) = \left| {\left( {{\rho _{\rm{p}}}\left( {t + \Delta {t_i} + \Delta \tau } \right) + \Delta \rho \cdot {{\widehat {\bf{n}}}_{{\rm{cag}}}} - {{\bf{\rho }}_i}\left( {t + \Delta {t_i}} \right)} \right) \times {{\widehat \ell }_{\rm{s}}}} \right|$(C.3)

where t + Δti represents the real observation epoch corrected for the camera time recording offset Δti; ^S${\hat \ell _{\rm{s}}}$ is the unit vector in the direction of the geocentric astrometric position of the occulted star; ρi is the geocentric geometric position of i; ρp is the geocentric astrometric position of Pluto obtained at t + Δti + Δτ, with Δτ accounting for the ephemeris offset along the geocentric motion of Pluto; n^cag${\hat n_{{\rm{cag}}}}$ is the unit vector in the northern direction along the predicted geocentric closest approach to the shadow axis, multiplied by Δρ to account for the ephemeris offset across the geocentric motion of Pluto. The two ephemeris offset parameters, Δτ and Δρ, represent the corrections to the epoch tcag and distance ρcag of the geoncetric closest approach predicted with the given reference ephemerides and star catalog (see Section 3.1).

Then, the values of r± can be numerically solved using the relationship between z and r±: z=±(r±+Dω(r±)),$z = \pm \left( {{r_ \pm } + D \cdot \omega \left( {{r_ \pm }} \right)} \right),$(C.4)

where D represents the light travel distance, which can be approximated by the geocentric distance of Pluto, and ω(r) is the total deviation angle at r.

As detailed in Dias-Oliveira et al. (2015) and Sicardy (2023), ω(r) can be obtained by a ray-tracing code that follows these steps:

  • set the temperature profile T(r) of Pluto using Equation (4) from Dias-Oliveira et al. (2015), with parameters obtained from Sicardy et al. (2016);

  • set the gas molecular refractivity K corresponding to a stellar wavelength of λµm;

  • set the atmospheric pressure p0 at the reference radius r0 of 1215 km either directly, or determined from a surface pressure psurf (at Rp = 1187 km) with the ratio psurf/p0 = 1.837 used by previous studies (Meza et al. 2019; Sicardy et al. 2021);

  • set the needed boundary condition for the gas molecular density, n0 = p0/(kB · T(r0)), where kB is Boltzmann’s constant;

  • integrate the first-order differential equation (Equations (5) and (6) of Dias-Oliveira et al. (2015)) to obtain the gas molecular density profile n(r);

  • transform n(r) to the gas refractivity v(r) using v(r) = K · n(r) (Equation (7) from Dias-Oliveira et al. (2015));

  • calculate the total deviation based on the straight line approximation (Equation (16) of Sicardy (2023), with the upper limit of r set as 2600 km).

Finally, an interpolator of ω(r) for specific p0 and K values is built in advance. This allows efficient computation of ω(r) during the analysis of the occultation observations.

Appendix D Results of the reanalyzed events necessary for comparisons

D.1 The 15 August 2018 event

Based on the IXON observational data from the online supplemental data file14 of Silva-Cabrera et al. (2022), we reanalyzed the 15 August 2018 event15 using our fitting procedure. As the raw data were normalized, with the flux contribution of the Pluto system carefully removed by Silva-Cabrera et al. (2022), all Ai and si were fixed as unity in our fitting procedure. Moreover, as the uncertainties used for weighting are not given, we used the standard deviation of the data outside the occultation part for each station. Results are presented in Figure D.1. Our surface pressure remeasurement is 12.0270.08+0.09$12.027_{ - 0.08}^{ + 0.09}$ µbar.

D.2 The 17 July 2019 event

Based on the TUHO observational data and uncertainties from Figure 1 of Arimatsu et al. (2020) (credit: Arimatsu et al, A&A, 638, L5, 2020, reproduced with permission © ESO), we reanalyzed the 17 July 2019 event using our fitting procedure. Results are presented in Figure A.1. Our surface pressure remeasurement is 9.4210.75+0.68$9.421_{ - 0.75}^{ + 0.68}$ µbar, similar to the result obtained by Arimatsu et al. (2020), 9.560.34+0.52$9.56_{ - 0.34}^{ + 0.52}$ µbar.

thumbnail Fig. D.1

Reanalysis of the 15 August 2018 event based on the IXON observational data from the online supplemental data file of Silva-Cabrera et al. (2022). The raw data were normalized, with the flux contribution of the Pluto system carefully removed by Silva-Cabrera et al. (2022). The uncertainties used for weighting are the standard deviation of the data outside the occultation part for each station. Panel (a): Reconstructed occultation map. Panel (b): Observed and simultaneously fitted light curves. Panel (c): Reconstructed stellar paths seen by IXON. Panel (d): The χ2 map, where χ2 denotes the goodness-of-fit value using these data. The best-fitting χ2 value per degree of freedom is 1.046.

References

  1. Arimatsu, K., Hashimoto, G. L., Kagitani, M., et al. 2020, A&A, 638, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Bertrand, T., & Forget, F. 2016, Nature, 540, 86 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bertrand, T., Forget, F., Umurhan, O. M., et al. 2018, Icarus, 309, 277 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bertrand, T., Forget, F., Umurhan, O. M., et al. 2019, Icarus, 329, 148 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bosh, A. S., Person, M. J., Levine, S. E., et al. 2015, Icarus, 246, 237 [NASA ADS] [CrossRef] [Google Scholar]
  6. Brosch, N. 1995, MNRAS, 276, 571 [NASA ADS] [Google Scholar]
  7. Brozovic, M., Showalter, M. R., Jacobson, R. A., & Buie, M. W. 2015, Icarus, 246, 317 [CrossRef] [Google Scholar]
  8. Buie, M. W., Leiva, R., Keller, J. M., et al. 2020a, AJ, 159, 230 [CrossRef] [Google Scholar]
  9. Buie, M. W., Porter, S. B., Tamblyn, P., et al. 2020b, AJ, 159, 130 [NASA ADS] [CrossRef] [Google Scholar]
  10. Desmars, J., Camargo, J. I. B., Braga-Ribas, F., et al. 2015, A&A, 584, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Desmars, J., Meza, E., Sicardy, B., et al. 2019, A&A 625, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Dias-Oliveira, A., Sicardy, B., Lellouch, E., et al. 2015, ApJ, 811, 53 [NASA ADS] [CrossRef] [Google Scholar]
  13. Elliot, J. L., Dunham, E. W., Bosh, A. S., et al. 1989, Icarus, 77, 148 [NASA ADS] [CrossRef] [Google Scholar]
  14. Elliot, J. L., Ates, A., Babcock, B. A., et al. 2003, Nature, 424, 165 [NASA ADS] [CrossRef] [Google Scholar]
  15. Forget, F., Bertrand, T., Vangvichith, M., et al. 2017, Icarus, 287, 54 [NASA ADS] [CrossRef] [Google Scholar]
  16. Gaia Collaboration 2022, VizieR Online Data Catalog: I/355 [Google Scholar]
  17. Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Grundy, W. M., Olkin, C. B., Young, L. A., & Holler, B. J. 2014, Icarus, 235, 220 [NASA ADS] [CrossRef] [Google Scholar]
  19. Gulbis, A. A. S., Emery, J. P., Person, M. J., et al. 2015, Icarus, 246, 226 [CrossRef] [Google Scholar]
  20. Hinson, D. P., Linscott, I. R., Young, L. A., et al. 2017, Icarus, 290, 96 [CrossRef] [Google Scholar]
  21. Holler, B. J., Yanez, M. D., Protopapa, S., et al. 2022, Icarus, 373, 114729 [NASA ADS] [CrossRef] [Google Scholar]
  22. Hubbard, W. B., Hunten, D. M., Dieters, S. W., Hill, K. M., & Watson, R. D. 1988, Nature, 336, 452 [NASA ADS] [CrossRef] [Google Scholar]
  23. Jacobson, R. A., Brozovic, M., Showalter, M., et al. 2019, in Pluto System After New Horizons, 2133, 7031 [NASA ADS] [Google Scholar]
  24. Lellouch, E., Butler, B., Moreno, R., et al. 2022, Icarus, 372, 114722 [NASA ADS] [CrossRef] [Google Scholar]
  25. Meza, E., Sicardy, B., Assafin, M., et al. 2019, A&A 625, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Morgado, B. E., Sicardy, B., Braga-Ribas, F., et al. 2021, A&A, 652, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Morgado, B. E., Gomes-Júnior, A. R., Braga-Ribas, F., et al. 2022, AJ, 163, 240 [NASA ADS] [CrossRef] [Google Scholar]
  28. Newell, D., & Tiesinga, E. 2019, in The International System of Units (SI) Special Publication (NIST SP), National Institute of Standards and Technology, Gaithersburg, MD, https://doi.org/10.6028/NIST.SP.330-2019 [CrossRef] [Google Scholar]
  29. Olkin, C. B., Young, L. A., Borncamp, D., et al. 2015, Icarus, 246, 220 [NASA ADS] [CrossRef] [Google Scholar]
  30. Park, R. S., Folkner, W. M., Williams, J. G., & Boggs, D. H. 2021, AJ, 161, 105 [NASA ADS] [CrossRef] [Google Scholar]
  31. Pasachoff, J. M., Souza, S. P., Babcock, B. A., et al. 2005, AJ, 129, 1718 [CrossRef] [Google Scholar]
  32. Pasachoff, J. M., Babcock, B. A., Durst, R. F., et al. 2017, Icarus, 296, 305 [CrossRef] [Google Scholar]
  33. Pavlov, H. 2020, Astrophysics Source Code Library [record ascl:2004.002] [Google Scholar]
  34. Pereira, C. L., Sicardy, B., Morgado, B. E., et al. 2023, A&A, 673, A4 [Google Scholar]
  35. Person, M. J., Dunham, E. W., Bosh, A. S., et al. 2013, AJ, 146, 83 [NASA ADS] [CrossRef] [Google Scholar]
  36. Person, M. J., Bosh, A. S., Zuluaga, C. A., et al. 2021, Icarus, 356, 113572 [NASA ADS] [CrossRef] [Google Scholar]
  37. Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 2007, Numerical Recipes in C++: The Art of Scientific Computing, 3rd edn. (Cambridge University Press) [Google Scholar]
  38. Rannou, P., & Durry, G. 2009, J. Geophys. Res. Planets, 114, E11013 [NASA ADS] [CrossRef] [Google Scholar]
  39. Sicardy, B. 2023, Comptes Rendus Physique, 23, 213 [NASA ADS] [CrossRef] [Google Scholar]
  40. Sicardy, B., Widemann, T., Lellouch, E., et al. 2003, Nature, 424, 168 [NASA ADS] [CrossRef] [Google Scholar]
  41. Sicardy, B., Bolt, G., Broughton, J., et al. 2011, AJ, 141, 67 [CrossRef] [Google Scholar]
  42. Sicardy, B., Talbot, J., Meza, E., et al. 2016, ApJ, 819, L38 [NASA ADS] [CrossRef] [Google Scholar]
  43. Sicardy, B., Ashok, N. M., Tej, A., et al. 2021, ApJ, 923, L31 [NASA ADS] [CrossRef] [Google Scholar]
  44. Silva-Cabrera, J. S., Castro-Chacón, J. H., Reyes-Ruiz, M., et al. 2022, MNRAS, 511, 5550 [NASA ADS] [CrossRef] [Google Scholar]
  45. Stern, S. A., Bagenal, F., Ennico, K., et al. 2015, Science, 350, aad1815 [Google Scholar]
  46. Washburn, E. W., West, C. J., Dorsey, N. E., et al. 1930, International Critical Tables of Numerical Data, Physics, Chemistry and Technology (Washington, DC: The National Academies Press) [Google Scholar]
  47. Yelle, R. V., & Elliot, J. L. 1997, in Pluto and Charon, eds. S. A. Stern, & D. J. Tholen (Kluwer Academic Publishers), 347 [Google Scholar]
  48. Young, E. F., French, R. G., Young, L. A., et al. 2008, AJ, 136, 1757 [NASA ADS] [CrossRef] [Google Scholar]
  49. Young, E., Young, L. A., Johnson, P. E., & PHOT Team 2021, AAS Meet. Abstr., 53, 307.06 [NASA ADS] [Google Scholar]

1

Sputnik Planitia is a large plain or basin on Pluto discovered by the New Horizons spacecraft (Stern et al. 2015).

2

The occulted star is Gaia DR3 6771712487062767488, of which the astrometric and photometric parameters are obtained from VizieR (Gaia Collaboration 2022).

9

As a technical note, if no equipment for reliable time references is used, Δτ should be fixed as zero, avoiding fitting problems caused by its linear relationship with ti.

13

The occulted star is Gaia DR3 6772059623498733952, of which the astrometric and photometric parameters are obtained from VizieR.

14

“stac401_Supplement_File - zip file”, via the paper link https://doi.org/10.1093/mnras/stac401.

15

The occulted star is Gaia DR3 6772629170525258240, of which the astrometric and photometric parameters are obtained from VizieR.

All Tables

Table 1

Circumstances and light-curve fitting results of the 5 September 2019 event.

Table A.1

Stations that encountered weather problems on 17 July 2019

Table A.2

Stations with positive detections on 5 September 2019

Table A.3

Stations that encountered weather problems on 5 September 2019

All Figures

thumbnail Fig. 1

Reconstructed occultation map of the 5 September 2019 event.

In the text
thumbnail Fig. 2

Occultation observations and the best-fitting light-curve model of the 5 September 2019 event. Panel a: observed and simultaneously fitted light curves. Panels b and c: reconstructed stellar paths seen by DWM and HNU, respectively.

In the text
thumbnail Fig. 3

The χ2 maps of the 5 September 2019 event. Panels a and b: The χobs2$\chi _{{\rm{obs}}}^2$ and χapr2$\chi _{{\rm{apr}}}^2$ maps, respectively. The χobs2$\chi _{{\rm{obs}}}^2$ denotes the goodness-of-fit only using observational data, while χapr2$\chi _{{\rm{apr}}}^2$ denotes the goodness-of-flt with the additional a priori information represented as the χ2 type value, (Δρρ)2. These maps are used to derive the best-fitting atmospheric pressure p0 at the reference radius r0 of 1215 km and the cross-track correction Δρ to the ephemerides, of which the a priori uncertainty σρ is 72 km. The surface pressure psurf and the geocentric closest approach distance ρcag to the shadow center are obtained by linear transformations of Δρ and p0, respectively. The best-fitting χobs2$\chi _{{\rm{obs}}}^2$ and χapr2$\chi _{{\rm{apr}}}^2$ values per degree of freedom are 1.160 and 1.155, respectively.

In the text
thumbnail Fig. 4

Pressure evolution of Pluto over a 248-yr heliocentric orbital period predicted by the VTM19, along with the measured pressures with 1σ error bars. The 3σ error bars of our new measurements and of the previous measurements of the 19 July 2016, 17 July 2019, and 6 June 2020 events are also presented using the published χ2 maps. AN2${A_{{{\rm{N}}_2}}}$ denotes the albedo of nitrogen ice.

In the text
thumbnail Fig. A.1

Reanalysis of the 17 July 2019 event based on the TUHO observational data from Figure 1 of Arimatsu et al. (2020) (credit: Arimatsu et al, A&A, 638, L5, 2020, reproduced with permission © ESO). Panel (a): Reconstructed occultation map. Panel (b): Observed and simultaneously fitted light curves. Panel (c): Reconstructed stellar paths seen by TUHO. Panel (d): The χ2 map, where χ2 denotes the goodness-of-fit value using these data. The best-fitting χ2 value per degree of freedom is 0.860.

In the text
thumbnail Fig. D.1

Reanalysis of the 15 August 2018 event based on the IXON observational data from the online supplemental data file of Silva-Cabrera et al. (2022). The raw data were normalized, with the flux contribution of the Pluto system carefully removed by Silva-Cabrera et al. (2022). The uncertainties used for weighting are the standard deviation of the data outside the occultation part for each station. Panel (a): Reconstructed occultation map. Panel (b): Observed and simultaneously fitted light curves. Panel (c): Reconstructed stellar paths seen by IXON. Panel (d): The χ2 map, where χ2 denotes the goodness-of-fit value using these data. The best-fitting χ2 value per degree of freedom is 1.046.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.